03. Explore Trans-Neptunian Objects (TNOs) in DP0.3 (Intermediate)#

Contact authors: Greg Madejski and Melissa Graham

Last verified to run: September 15, 2023

Targeted learning level: Intermediate

Introduction#

This tutorial demonstrates how to identify and explore a population of Trans-Neptunian Objects (TNOs) in the simulated DP0.3 catalogs. TNOs are defined by having orbits with semi-major axes beyond the orbit of Neputne (> 30.1 AU). As the semi-major axis (a) can be derived from the orbit’s ellipticiy (e) and perihelion distance (q) via a = q/(1. - e), and as both ellipticity and perihelion are available in the MPCORB table, a sample of TNOs can be identified in the DP0.3 data set (see Step 1). Their properties (specifically, the relationship between their semi-major axis and eccentricity, as well as the distribution of their derived diameters) will be explored in Step 2. Note that some of the objects might not be moving in elliptical orbits (e > 1 - meaning they are not bound to the Solar System, but moving on parabolic or hyperbolic orbits). Such objects will be excluded from this tutorial, as an application of the formula above would result in a negative value of a.

Compared to the Solar System objects closer to the Earth, such as Main Belt Asteroids or Near-Earth Objects (NEOs), TNOs move relatively slowly across the sky. This relatively slow movement means that TNOs that fall within an LSST Deep Drilling Field (DDF) can stay within that field, and LSST can accumulate thousands of observations of them. This tutorial explores the position on the sky of one such TNO (Step 3) and plots time-domain quantities such as magnitude and phase angle (Step 4). Finally, it provides a visualization of its trajectory projected into 2D (see Step 5).

More information about the LSST DDFs can be found on the LSST DDF webpage and in Section 2.6 of the Survey Cadence Optimization Committee’s Phase 2 Recommendations report (PSTN-055). Note that DP0.2 did not include DDF observations, so the ability to explore science with a DDF-like cadence is unique to the DP0.3 simulation.

This tutorial assumes the successful completion of the beginner-level DP0.3 Portal tutorials, and uses the Astronomy Data Query Language (ADQL), which is similar to SQL (Structured Query Language). For more information about the DP0.3 catalogs, tables, and columns, see the DP0.3 Data Products Definition Document (DPDD).

Step 1. Identify a population of TNOs#

1.1. Log into the Rubin Science Platform at data.lsst.cloud and select the Portal Aspect. At upper right, next to “TAP Services” choose to “Show”, and then select “LSST DP0.3 SSO” from the drop-down menu that appears at the top.

1.2. At upper right, next to “View” choose “Edit ADQL”. Enter the following ADQL statement into the ADQL Query box. It will return the ellipticity (e), perihelion distance (q), and inclination (incl) for a random subset of objects in the MPCORB table. For an explanation of why this constraint on ssObjectId returns a random sample, see Step 2 of DP0.3 Portal tutorial 01, “Introduction to DP0.3: the MPCORB and SSObject tables”.

SELECT e, q, incl
FROM dp03_catalogs_10yr.MPCORB
WHERE ssObjectId > 9000000000000000000

1.3. Set the “Row Limit” to be 200000 and click “Search”.

1.4. The default results view will show a heatmap plot of q vs. e at left, and the table view at right.

A screenshot of the default results view for the query.

The default results view for the query, with the table at right and the heatmap at left.#

1.5. Exclude the objects moving on unbound orbits. Note that a small fraction of the objects - roughly one in a thousand - have derived eccentricities > 1 meaning those are not bound to the Solar System. Those objects can be excluded from further analysis by entering < 1 in the box underneath the table heading e, and hitting “enter.” This will result in a slightly modified display as below.

A screenshot of the default results view for the modified query.

The view for the query with e < 1.#

1.6. Create a new column in the table, containing semi-major axis, a. In the upper right column of the table panel, click on the icon to add a column (a tall narrow rectangle to the left of a + sign). In the pop-up window to “Add a column”, set the “Name” to “a”, the “Expression” to “q/(1-e)”, the “Units” to “au”, and the “Description” to “semi-major axis”. Click “Add Column”, and see the new column appear in the table.

A screenshot of the pop-up window to add a column.

The “Add a column” pop-up window.#

1.7. Create a scatter plot of inclination vs. semi-major axis. In the plot panel, click the “Settings” icon (double gears), and select “Add New Chart”. Set the “Plot Type” to “Scatter”, the “X” to “a”, “Y” to “incl”. In the “Chart Options” dropdown menu, set the “X Min” to “0”, the “X Max” to 60, the “Y Min” to 0, and the “Y Max” to 80. Click “OK”.

A screenshot of the plot parameters pop-up window.

Create a new plot with these parameters.#

1.8. Delete the default plot by clicking on the blue cross in the upper right corner, so that only the newly-created plot appears (it should look like the plot below). TNOs appear as a distinct population with a > 30.1 au in this parameter space.

A screenshot of the inclination versus semi-major axis plot, showing a clear population of TNOs.

The population of TNOs has x-values greater than 30 au.#

1.9. Notice how in the plot above, the majority of objects returned by the query were closer to the Sun than 30.1 au. In fact, only about 600 of the moving objects from the query were TNOs. The total number of TNOs is estimated to be far less than Main Belt Asteroids, and due to their distance they move more slowly and are fainter, and so are harder to detect and characterize. In the next step, a revised query will be used to only retrieve objects with semi-major axis greater than 30.1 au.

Step 2. Explore the properties of a population of TNOs#

2.1. Now that the population of the Trans-Neptunian Objects has been identified, it is possible to further explore their properties. To study the properties of a larger sample of TNOs, return to the ADQL query interface by clicking on “RSP TAP Search” tab, and clicking on “Edit ADQL” button.

2.2. Clear the ADQL query, and execute a query below, simiar to the one in Step 1.2, but which includes only objects at a > 30.1 au. Also include the absolute H magnitude mpcH which will be used in the derivation of diameters of TNOs in the subsequent step (2.6) below. As TNOs aren’t the only solar system objects beyond Neptune, reject objects with mpcDesignation as Long Period Comets (LPC).

SELECT e, incl, q, mpcH, mpcDesignation
FROM dp03_catalogs_10yr.MPCORB
WHERE q / (1 - e) > 30.1
AND SUBSTRING(mpcDesignation, 1, 3) != 'LPC'

2.3. Keep the “Row limit” to 200000, and click “Search”. This query will return 62,961 objects. The default plot in the results view will be a heatmap of inclination vs. eccentricity.

2.4. Plot the eccentricity of the orbit e as a function of the semi-major axis a. This time (in contrast to Step 1.6 but accomplishing the same goal), calculate a from e and q via setting appropriate plot parameters rather than creating another column in the right-hand table. Start by clicking on the “plot settings” (two gears) on the left-hand panel, and click on “Add New Chart.”

2.5. In the “Plot Parameters” pop-up window, select “Heatmap” for the plot type. Enter “q/(1-e)” for the X-axis, and “e” for the y-axis. Choose the “YlOrRd” (yellow-orange-red) color scale. Increase the number of bins to 200 for both x and y to improve the resolution of the heatmap. Expand the “Chart Options” and set the titles and labels as below. Restrict the x-axis to 30 < a < 100 au to illustrate in more detail the region of 1 to three times Neptune’s semi-major axis.

A screenshot of the plot parameters for the eccentricity vs. semi-major axis plot

The plot parameters for the eccentricity vs. semi-major axis plot.#

2.6. Click on “OK” or “Apply” in the “Plot Parameters” window, and view the plot (see below). Delete the default plot of inclination vs. eccentricity as it is not needed.

A screenshot of the plot of the eccentricity vs. semi-major axis

The plot of eccentricity vs. semi-major axis of solar system objects with a semi-major axis > 30 au.#

2.7. Multiple populations are apparent in the above plot. The majority of the objects have low eccentricity and a semi-major axis of about 30 to about 50 au. There are several sub-populations of Trans-Neptunian Objects (TNOs), such as the classical, resonant, scattered-disk, and detached. A full review of all TNO populations is beyond the scope of this tutorial.

2.8. Estimate the diameters of the objects using their absolute H magnitudes. Where H is the absolute H magnitude (column mpcH), and A is the albedo, the diameter \(d\) in kilometers is \(d = 10^{(3.1236 - 0.5 \times log(A) - 0.2 \times H)}\). This tutorial adopts an albedo value of 0.15 (as is commonly adopted, e.g., Vilenius et al. 2012; the diameter depends only weakly on the albedo), with which the expression reduces to \(d = 10^{(3.536 - (0.2 \times H))}\) km. Click on the “add column” icon. Enter D in the “name” field, and power(10,(3.536 - 0.2 * mpcH)) in the expression field, as below. Click the “Add Column” button.

screenshot illustrating the expression needed to make the new column containing the diameter

How to add a new column containing the estimated diameter.#

2.9. Plot the distribution of estimated diameters in log-space. Open the “Plot Parameters” pop-up window, click on “Add New Chart,” select “Histogram”, and enter the parameters as below.

screenshot illustrating the plot parameters for displaying the distribution of estimated diameters

The plot parameters for the distribution of estimated diameters.#

2.10. Click on the “Apply” or “OK” button. This will result in the plot showing the distribution of estimated diameters.

screenshot illustrating the distribution of estimated diameters

The distribution of estimated diameters.#

2.11. Notice that the tail of the distribution extends to very small diameters. This is suprising, as detecting kilometer-sized objects at the distance of Neptune should be quite challenging. This illustrates, in part, some of the shortcomings of the assumptions (such as albedo) which were used to estimate the diameters. See also the first exercise for the learner in Step 6.

Step 3. Find and explore a well-observed TNO#

3.1. Return to the RSP TAP Search form by clicking on the RSP TAP Search button at upper left. Navigate to the ADQL interface by clicking on the “Edit ADQL” button.

3.2. To identify a distant Solar System object with a large number of observations, enter the query below. This query joins the MPCORB table with the DiaSource table in order to retrive the number of detections: the count of the number of DiaSource table rows for a given solar system object, each of which has a unique ssObjectId. This query also applies a constraint that the semi-major axis be between 30 and 100 AU, and a constraint on the ssObjectId to return a random subset (similar to Step 1.2).

SELECT mpc.ssObjectId, COUNT(ds.ssObjectId), mpc.e, mpc.q
FROM dp03_catalogs_10yr.MPCORB AS mpc
JOIN dp03_catalogs_10yr.DiaSource AS ds ON mpc.ssObjectId = ds.ssObjectId
WHERE mpc.ssObjectId < -700000000000000000
AND mpc.q > 30 * (1 - mpc.e)
AND mpc.q < 100 * (1 - mpc.e)
GROUP BY mpc.ssObjectId, mpc.e, mpc.q

3.3. Click on “Search”. This search might take up to a minute. The query returns 12,589 objects.

3.4. The default view of the search above is a plot the first two columns against each other, ssObjectId and COUNT, which is not a particularly useful plot aside from showing that the number of detections for the most oft-detected objects in the outer Solar System is in the thousands. Click twice on the COUNT column header to order the entries by descending count and identify the most oft-detected outer Solar System object.

A screenshot of the default results view with the table sorted by count.

The default results view from the ADQL query above.#

3.5. Continue with the object with the largest number of observations: ssObjectId = -735085100561880491, which was detected 12,103 times. Its modest eccentricity of 0.1512 implies that this is a TNO (unlikely to be a comet).

3.6. Return to the ADQL query interface and use the following statement to retrieve the sky coordinates, magnitudes, filter (band), and time of observations (midPointMjdTai) for the oft-observed TNO with ssObjectId as above.

SELECT ra, dec, mag, band, midPointMjdTai
FROM dp03_catalogs_10yr.DiaSource
WHERE ssObjectId = -735085100561880491

3.7. The default results view will show the “Coverage” map at upper left. In the future, with real LSST data, this map would have an underlay of the LSST deeply stacked image. Since DP0.3 has no images, the “Coverage” map only shows the overlay of RA vs. Dec, which is redundant with the default plot. At upper right, click on “Bi-view Tables” to hide the “Coverage” map and show only the table and plot.

The default results view after clicking on bi-view tables.

The “Bi-view Tables” results view for the query of ssObjectId = -735085100561880491.#

3.8. Set the color of individual points to represent the time of the observation to better illustrate how the object moves across the sky. In the plot panel, click on the “Settings” icon (double gears) to open the “Plot Parameters” pop-up window. Under “Trace Options”, for “Color Map” enter “midPointMjdTai” and for “Color Scale” enter “Rainbow”. Then click “Apply”.

A screenshot of the plot of sky coordinates colored as a function of time.

Purple color corresponds to earlier observtations, and the red color corresponds to later observations.#

3.9. In the plot above, the 10 loops in the object’s path on the sky is a result of Earth’s orbital period and the 10-year LSST duration. As described in the introduction, this particular TNO was detected by LSST over ten thousand times because it happened to be in a deep drilling field. This will not be the case for the majority of solar system objects.

Step 4. Plot the time-domain quantities for the TNO#

Note that no time domain evolution in object brightness was included in the DP0.3 simulation (e.g., rotation curves for non-spherical objects, outgassing events). All changes in the brightness of DP0.3 objects with time are due to changes in the distance and phase angle from Earth.

4.1. Return to the search form and execute the following ADQL query to retrieve the r-band magnitudes, phase angles, heliocentric and topocentric distances, and time of the observations for the TNO explored in Step 3.

SELECT ds.midPointMjdTai, ds.mag, ds.band,
ss.phaseAngle, ss.topocentricDist, ss.heliocentricDist
FROM dp03_catalogs_10yr.DiaSource AS ds
JOIN dp03_catalogs_10yr.SSSource AS ss ON ds.diaSourceId = ss.diaSourceId
WHERE ss.ssObjectId = -735085100561880491
AND ds.band = 'r'

4.2. The default plot will have the r-band magnitude as a function of time. Use the plot “Settings” function to add a scatter plot showing the phase angle as a function of time. For the x-axis, use midPointMjdTai - 60000 to show more clearly the timescales between observations.

4.3. As mentioned above, the simulated Solar System data does not include any time-varying features. The changes in apparent magnitude are due to the object changing in phase angle and distance from Earth as a function of time. Add two new scatter plots showing the r-band magnitude as a function of phase angle and as a function of topocentric (Earth-centered) distance, as is shown below. The results view for four plots automatically reconfigures to a two-by-two grid. Notice how the magnitude is a monotonic function of phase angle and distance, but not time.

A screenshot of four plots showing magnitude and phase angle are not correlated with time, and that magnitude is correlated with phase angle and distance from Earth.

Four plots demonstrating that the apparent magnitude depends on phase angle and distance from Earth.#

4.4. Plot the topocentric and heliocentric distances of the object as a function of time already retrieved in Step 4.1. First, delete all but one of the plots prepared in Step 4.3 by clicking on the blue X in the upper right-hand part of the plot panels to make space for new plots. Then add a pair of new scatter plots that show topocentricDist and heliocentricDist as a function of midPointMjdTai - 60000. Then delete the remaining old plot so that only the two new plots are displayed.

A screenshot of two plots showing the heliocentric and topocentric distance of the trans-Neptunian object as a function of time.

Heliocentric and topocentric distance of the TNO as a function of time.#

4.5. The left plot shows the periodic change of the topocentric distance with time resulting from the Earth’s motion around the Sun - a different view of the same effect seen in Step 3. The right plot shows that this object is on a slightly inbound trajectory with respect to the Sun.

Step 5. View the 2-D projection of the TNO’s orbit to visualize its 3-D trajectory#

5.1. The goal of Step 5 is to visualize the 3-D trajectory of the well-observed trans-Neptunian object, via viewing the projections of its 3-D helio- and topocentric distances as a function of time into 2-D. Navigate to the ADQL query interface. Execute the query below to extract the helio- and topocentric X, Y, and Z distances of the TNO - so you can visualize its trajectory.

SELECT heliocentricX, heliocentricY, heliocentricZ,
topocentricX, topocentricY, topocentricZ, ssObjectId
FROM dp03_catalogs_10yr.SSSource
WHERE ssObjectId = -735085100561880491

5.2. The default plot will be the heliocentric Y distance as a function of heliocentic X distance as in the screenshot below. Note that the object moves slowly in heliocentric coordinate X as well as in Y (by a comparison to, e.g., Earth’s motion), covering only a few au in 10 years. This is expected given its multi-au distance from the Sun.

A screenshot of a plot showing the heliocentric Y vs. heliocentric X distance of the trans-Neptunian object.

Heliocentric Y vs. X distance of the trans-Neptunian object.#

5.3. Now plot the heliocentric Z distance as a function of heliocentric X distance. Click on “Plot Settings” and click on “Add New Chart.” Select heliocentricZ for y and heliocentricX for x. Click on “Apply” or “OK.”

5.4. Observe that the object’s trajectory is not constant in Z - and that means that its orbit is not in the plane of the Ecliptic during the simulated Rubin observation, but the object does pass through the ecliptic plane when Z = 0.

A screenshot of plots showing the heliocentric Y and heliocentric Z vs. heliocentric X distance of the trans-Neptunian object as a function of time.

Heliocentric Y vs. X distance as well as helliocentric Z vs. X distance of the trans-Neptunian object as a function of time.#

5.5. Next, plot the topocentricY vs. topocentricX and topocentricZ vs. topocentricX distances. On the same screen where you generated the plots in previous two steps, click on “Plot Settings” and click on “Add New Chart.” First, select topocentricY for y and topocentricX for x. and click “Apply” or “OK.” Next, click on “Plot Settings” and click on “Add New Chart.” Select topocentricZ for y and topocentricX for x, and click “Apply” or “OK.” There, the effect of position of the TNO on the sky as a result of Earth’s orbital motion is clearly apparent.

A screenshot of four plots showing the heliocentric and topocentric distances of the trans-Neptunian object as a function of time.

Visualization of the 3-D TNO’s trajectory by viewing the 2-D projections of its trajectory as measured from the Sun (top two plots) and the Earth (bottom two plots).#

Step 6. Exercises for the learner#

6.1. In Step 2, some of the sizes of the TNOs were on order one kilometer, quite small for objects at the distance of Neptune. However, objects with high eccentricities could come closer to Earth, and be detected despite their small size. For the objects returned by the query in Step 2, plot the eccentricity vs. estimated diameter. Explore whether some of the smallest objects have large eccentricities.

6.2. Plot the histogram of the number of visits to the Solar System objects in the dp03_catalogs.SSObject for objects observed more than 1000 times.

6.3. Repeat the steps 4 and 5 for another object with a large number of observations (say another one with numObs > 2,000). Note that you already identified objects with large number of observations in Steps 3.1, 3.2, and 3.3.