Creating a HR Diagram with calculations

Let's create and compare a few color magnitude diagrams.

Note

This guide is slightly outdated, but everything still applies. The plots just look different, really. For example, all references to best here should be substituted with mwmlite.

To start, we'll load the mwmlite reductions, which contain SDSS-V's best parameter estimates for all stars. This is automatically done when you load the app.

image_20241211095950

We can create a new heatmap or scatter plot with the Add View button.

image_20241211100053

Let's add a heatmap. This will automatically create a basic HR (teff vs. logg; without flipped axes), but we can flip them very easily with the plot settings menu.

image_20241211100458

Heatmap with expanded settings menu. Both axes are set to be flipped.

Let's say we want to plot Gaia magnitudes, say BP - RP colour against absolute Gaia G magnitude. There isn't a bp_rp nor a G absolute magnitude column in the SDSS-V dataset by default, but there is the relevant information to make it.

To solve this, we can use the Virtual Calculations Menu in the side bar to create both a bp_rp and g_abs_mag column.

image_20241211101144

Note the virtual calculation menu on the bottom left in the sidebar

We'll say that bp_rp is simply:

\[ BP - RP \]

image_20241211100833

Adding a virtual column under the title bp_rp, which is just the BP - RP color index.

Likewise, the absolute G band magnitude is just:

\[ m_{G} - 5\log_{10}(\bar{\omega} / 1000) \]

since the parallax column plx is given in arcsec.

image_20241211101237

Both virtual columns now appear in the menu, and can be deleted as needed

When we go back to our plot settings menu, we can now actually see the new columns we added!

image_20241211101326

Look at them!

Setting the plot to both of these gives us our CMD.

image_20241211101944

Since the diagram was pretty uninformative, we'll change it to plot the logarithmic count of sources. We'll also filter down to just Milky Way Mapper (mwm) targets with the Targeting Filters.

image_20241211102143

If we restrict it to high S/N targets, and change the colorscale, we can clearly see the main tracks and the red clump. Note the effect of reddening on the color index.

image_20241211102258

We can additionally add reddening corrections, based on any of the \(E(B-V)\) columns, derived from a variety of dustmaps.

image_20241211102938

Using the Column Glossary, we can look up all the relevant reddening columns in our dataset.

Let's use the ebv_edenhofer_2023 column. Note that this limits us to bright sources within 1.25 kpc. Nonetheless, we get a much cleaner set of sources with the reddening corrected color. image_20241211102838

Pretty cool CMD, using the reddening-corrected color.


Written by Riley Thai, Dec 11 2024