Data behind Figures 1-3 in "The doubly librating Plutinos"

Data behind the figures for:

Article title: The doubly librating PlutinosAuthors: Renu Malhotra ([e-mail]) and Takashi Ito

Accepted for publication in Astrophysical Journal (Jan 2025)

This depository consists of three data files for data behind Figures 1--3.Each of these three files has two header rows, followed by data rows.

dbf1.txt: plain text file with six columns of data for Figure 1      column 1 ... designation of each Plutino      column 2 ... osculating semimajor axis (au) at initial time      column 3 ... osculating eccentricity at initial time      column 4 ... osculating inclination (degrees) at initial time      column 5 ... dsigma = sigma_max - sigma_min  (degrees)      column 6 ... dg = g_max - g_min  (degrees)

      Sources:       column 1 ... Volk K., Van Laerhoven C., 2024, Res. Notes AAS 8 36, data file rnaasad22d4t1_mrt.txt https://doi.org/10.3847/2515-5172/ad22d4      columns 2,3,4 ... osculating orbital elements for epoch A.D. 2024-Oct-17 00:00:00.0000 (JDTDB 2460600.5); based on cartesian heliocentric elements obtained from JPL Horizons System data service (retrieved on Nov 21, 2024), converted to barycentric elements, with solar system invariable plane as the reference plane      column 5,6 ... from output of 100 myr numerical integrations with SWIFT_RMVS; Plutinos that did not survive to the full 100 myr have dsigma and dg given the false value of 1.0e99

dbf2.txt: plain text file with seven columns of data for Figure 2      column 1 ... designation of each Plutino      column 2 ... osculating eccentricity at initial time      column 3 ... osculating inclination (degrees) at initial time      column 4 ... time averaged eccentricity (over 100 myr)      column 5 ... time averaged inclination (degrees) (over 100 myr)      column 6 ... dsigma = sigma_max - sigma_min  (degrees)      column 7 ... dg = g_max - g_min  (degrees)

      Sources:      column 1 ... Volk K., Van Laerhoven C., 2024, Res. Notes AAS 8 36, data file rnaasad22d4t1_mrt.txt https://doi.org/10.3847/2515-5172/ad22d4      columns 2,3 ... osculating orbital elements for epoch A.D. 2024-Oct-17 00:00:00.0000 (JDTDB 2460600.5); based on cartesian heliocentric elements obtained from JPL Horizons System data service (retrieved on Nov 21, 2024), converted to barycentric elements, with solar system invariable plane as the reference plane      column 4,5 ... time averaged values from output of 100 myr numerical integrations with SWIFT_RMVS      column 6,7 ... from output of 100 myr numerical integrations with SWIFT_RMVS; Plutinos that did not survive to the full 100 myr have dsigma and dg given the false value of 1.0e99

dbf3.txt: plain text file with seven columns of data for Figure 3      column 1 ... designation of each Plutino (only those identified as g librators)      column 2 ... time averaged eccentricity (over 100 myr)      column 3 ... time averaged inclination (degrees) (over 100 myr)

      column 4,5 ... eccentricity and inclination (degrees), respectively for the black dot-dashed curve      column 6,7 ... eccentricity and inclination (degrees), respectively for the red curve

      Sources:      column 1 ... Volk K., Van Laerhoven C., 2024, Res. Notes AAS 8 36, data file rnaasad22d4t1_mrt.txt https://doi.org/10.3847/2515-5172/ad22d4      column 2,3 ... time averaged values from output of 100 myr numerical integrations with SWIFT_RMVS      column 4,5 ... from numerical quadrature of the disturbing function of the CR3BP (Sun and Neptune in circular orbits, and a massless Plutino) at Neptune's exterior 3:2 mean motion resonance      column 6,7 ... from numerical quadrature of the disturbing function of the CR6BP (Sun, Jupiter, Saturn, Uranus and Neptune fixed in circular and coplanar orbits, and a massless Plutino-like object) at Neptune's exterior 3:2 mean motion resonance