Skip to content

Commit

Permalink
Update portal-dp0-3-3.rst
Browse files Browse the repository at this point in the history
  • Loading branch information
sgreenstreet authored Feb 27, 2025
1 parent 4f5b063 commit 278fe1c
Showing 1 changed file with 18 additions and 19 deletions.
37 changes: 18 additions & 19 deletions tutorials-dp0-3/portal-dp0-3-3.rst
Original file line number Diff line number Diff line change
Expand Up @@ -111,7 +111,7 @@ This will result in a slightly modified display as below.

**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",
In the pop-up window to "Add a column", set the "Name" to "a", the "Expression" to "q/(1.0-e)", the "Units" to "au",
and the "Description" to "semimajor axis".
Click "Add Column", and see the new column appear in the table.

Expand Down Expand Up @@ -160,11 +160,11 @@ In the next step, a revised query will be used to only retrieve objects with sem
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.
**2.1.** Now isolate the population of Transneptunian Objects and further explore their properties.
To study the properties of a larger sample of TNOs, return to the ADQL query interface by clicking on "DP0.3 Catalogs" 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.
**2.2.** Clear the ADQL query, and execute the query below, which is simiar to the one in Step 1.2 but includes only objects at ``a`` > 30.1 au.
Also include the absolute H magnitude ``mpcH`` which will be used in the derivation of TNO diameters 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).

Expand All @@ -180,44 +180,43 @@ Long Period Comets (LPC).
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``.
**2.4.** Plot the eccentricity of the orbit ``e`` as a function of the semimajor 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 "+" sign on the left-hand panel to add a 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.
**2.5.** In the "Add New Chart" pop-up window, select "Scatter" for the plot type.
Enter "q/(1.0-e)" for the X-axis, and "e" for the y-axis.
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.
Restrict the x-axis to ``30 < a < 100`` au.

.. figure:: /_static/portal_tut03_step02a.png
:width: 400
:name: portal_tut03_step02a
:alt: A screenshot of the plot parameters for the eccentricity vs. semi-major axis plot
:alt: A screenshot of the plot parameters for the eccentricity vs. semimajor axis plot

Figure 6: The plot parameters for the eccentricity vs. semi-major axis plot.
Figure 6: The plot parameters for the eccentricity vs. semimajor axis plot.

**2.6.** Click on the "OK" button in the "Plot Parameters" window, and view the plot (see below).
**2.6.** Click on the "OK" button in the "Add New Chart" window, and view the plot (see below).
Delete the default plot of inclination vs. eccentricity as it is not needed.

.. figure:: /_static/portal_tut03_step02b.png
:width: 600
:name: portal_tut03_step02b
:alt: A screenshot of the plot of the eccentricity vs. semi-major axis

Figure 7: The plot of eccentricity vs. semi-major axis of solar system objects with a semi-major axis > 30 au.
Figure 7: The plot of eccentricity vs. semimajor axis of Transneptunian Objects (TNOs).

**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.7.** Multiple sub-populations are apparent in the above plot.
The majority of the objects have low eccentricity (``e``<0.3) and a semimajor axis of about 30 to about 50 au.
There are several sub-populations of Transneptunian Objects, such as the classical, resonant, scattered-disk, and detached sub-populations.
A full review of all TNO sub-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 :math:`d`
in kilometers is :math:`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 <https://arxiv.org/pdf/1204.0697.pdf>`_;
the diameter depends only weakly on the albedo),
This tutorial adopts an albedo value of 0.15 (as is commonly adopted, e.g., `Vilenius et al. 2012 <https://arxiv.org/pdf/1204.0697.pdf>`_),
with which the expression reduces to :math:`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.
Expand All @@ -240,7 +239,7 @@ Click on the "+" sign in the pop-up window, click on "Add New Chart," select "Hi

Figure 9: The plot parameters for the distribution of estimated diameters.

**2.10.** Click on the "Apply" button.
**2.10.** Click on the "OK" button.
This will result in the plot showing the distribution of estimated diameters.

.. figure:: /_static/portal_tut03_step02e.png
Expand Down

0 comments on commit 278fe1c

Please sign in to comment.