Characterizing Variable Stars with APOGEE + TESS

Characterizing Variable Stars with APOGEE + TESS#

Learning Goals#

By the end of this tutorial, you will:

  • Understand how to search the MAST Archive and download SDSS APOGEE data using astroquery.mast

  • Plot a stellar spectrum and an HR diagram and using APOGEE data

  • Understand how to search for TESS light curves complementing the APOGEE observations

  • Explore different types of variable stars using their APOGEE parameters and TESS light curves

Table of Contents#

  • Introduction

  • Imports

  • Accessing APOGEE data at MAST

    • Querying all APOGEE data

    • Searching for a specific star

    • Downloading APOGEE data products

    • Plotting an APOGEE spectrum

    • Downloading the APOGEE allStar catalog

    • Plotting an HR diagram from APOGEE

  • Searching for TESS data of this star

    • Coordinate search using astroquery.mast

    • Plotting a light curve from TESS

  • Combining APOGEE and TESS data

    • Plotting APOGEE and TESS data together

    • Exploring Different types of Variable Stars

      • RR-Lyrae Variable: 2M11361176+8117369

      • Rotationally Variable - 2M19181706+5141323

      • Eclipsing Binary - 2M19203184+3830492

  • Additional Resources

    • How to Cite

    • About This Notebook

Introduction#

The Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey provides infrared-wavelength high-resolution spectroscopy for over 650,000 unique stars from the Milky Way and nearby dwarf galaxies. APOGEE collected data between 2011 - 2020 as part of the Sloan Digital Sky Survey (SDSS-IV) project. APOGEE data is now available at the Mikulski Archive for Space Telescopes (MAST) through the SDSS Legacy Archive at MAST.

In this notebook tutorial, we will demonstrate how to access APOGEE data at MAST using Python. One APOGEE star, a Cepheid Variable named V1154-Cyg will be used to demonstrate the basics of how to download and plot APOGEE data. We will then combine this APOGEE spectrum with light curves from the TESS mission, also accessible from MAST, to study this variable star.

Imports#

The main packages we’re using for this notebook and their use-cases are:

  • astroquery.mast Observations for searching the MAST archive

  • astropy.io fits for accessing FITS files

  • astropy.table Table for accessing FITS tables

  • astropy.units for working with astronomical units

  • astropy.coordinates for handling astronomical coordinates

  • matplotlib.pyplot for plotting data

  • numpy to handle array functions

%matplotlib inline

from astroquery.mast import Observations
import astropy.io.fits as fits
from astropy.table import Table
from astropy.coordinates import SkyCoord

import astropy.units as u
import matplotlib.pyplot as plt
import numpy as np

This cell updates some of the settings in matplotlib to use larger font sizes in the figures:

#Update Plotting Parameters
params = {'axes.labelsize': 12, 'xtick.labelsize': 12, 'ytick.labelsize': 12, 
          'text.usetex': False, 'lines.linewidth': 1,
          'axes.titlesize': 18, 'font.family': 'serif', 'font.size': 12}
plt.rcParams.update(params)

Accessing APOGEE data at MAST#

The SDSS Legacy Archive at MAST hosts all of the science-ready data products from the SDSS-IV APOGEE Survey, which includes spectra for over 650,000 unique stars in the Milky Way. APOGEE acquired observations in both the northern and southern hemispheres, targeting the disk, bulge, bar, and halo components of the Milky Way, as well as several hundred stars in nearby satellite galaxies (including the LMC and SMC). This notebook will demonstrate how to search and download APOGEE data using MAST!

Querying all APOGEE data#

Searching for APOGEE data is straightforward with astroquery.mast. In this example, we use Observations.query_criteria and search for provenance_name = 'APOGEE'. This will return a table describing all of the APOGEE data hosted by the MAST archive.

Other useful search parameters for APOGEE data might include:

  • obs_collection = 'SDSS': searches for all SDSS data

  • To specify telescope, use the instrument_name field. For example, instrument_name = 'apogee-lco25m' limits the search to the stars that were observed with the Las Campanas Observatory 2.5-m telescope in the southern hemisphere.

  • Use target_name to search for stars using their APOGEE ID (usually the 2MASS ID - for example '2M05320041-0017041')

  • obs_id can help search for specific targets or fields. Note that wild cards (*) are allowed in the search fields, for example:

    • obs_id='*n7789*' for anything in the “N7789” field

    • obs_id='*2M05320041-0017041*' for a specific star name (2MASS ID)

# Search for APOGEE data
# Using the pagesize and page parameters to only return first 10 results
apogee_obs_list = Observations.query_criteria(provenance_name='APOGEE', pagesize=10, page=1)

# Display First Ten Entries in Table
apogee_obs_list
Table masked=True length=10
intentTypeobs_collectionprovenance_nameinstrument_nameprojectfilterswavelength_regiontarget_nametarget_classificationobs_ids_ras_decdataproduct_typeproposal_picalib_levelt_mint_maxt_exptimeem_minem_maxobs_titlet_obs_releaseproposal_idproposal_typesequence_numbers_regionjpegURLdataURLdataRightsmtFlagsrcDenobsidobjID
str7str4str6str13str4str4str8str18str4str51float64float64str8str18int64float64float64float64float64float64str63float64str3str1int64str51str1str92str6boolfloat64str9str9
scienceSDSSAPOGEEapogee-lco25mSDSSNoneINFRARED2M06015966-6511312STARsdss_apogee_lco25m_275-29-o_tess_2m06015966-651131290.498613-65.192024spectrumSDSS Collaboration358803.1698020254658803.169802025464000.01514.01696.0Apache Point Observatory Galactic Evolution Experiment (APOGEE)59554.0N/A----CIRCLE 90.498613 -65.192024 0.00018194444444444445--mast:SDSS/apogee/lco25m/275-29-O_TESS/2M06015966-6511312/asStar-dr17-2M06015966-6511312.fitsPUBLICFalsenan234526617675892862
scienceSDSSAPOGEEapogee-lco25mSDSSNoneINFRARED2M19040099-3105425STARsdss_apogee_lco25m_sgrc-3_2m19040099-3105425286.004128-31.095142spectrumSDSS Collaboration358272.4046338773158386.0883347569520000.01514.01696.0Apache Point Observatory Galactic Evolution Experiment (APOGEE)59554.0N/A----CIRCLE 286.004128 -31.095142 0.00018194444444444445--mast:SDSS/apogee/lco25m/SGRC-3/2M19040099-3105425/asStar-dr17-2M19040099-3105425.fitsPUBLICFalsenan234602474676189155
scienceSDSSAPOGEEapogee-apo25mSDSSNoneINFRARED2M04164154+6021333STARsdss_apogee_apo25m_146+08_2m04164154+602133364.17308660.359257spectrumSDSS Collaboration357375.19290077546657405.10632836805612000.01514.01696.0Apache Point Observatory Galactic Evolution Experiment (APOGEE)59554.0N/A----CIRCLE 64.173086 60.359257 0.0002777777777777778--mast:SDSS/apogee/apo25m/146+08/2M04164154+6021333/apStar-dr17-2M04164154+6021333.fitsPUBLICFalsenan234102062674246652
scienceSDSSAPOGEEapogee-apo25mSDSSNoneINFRARED2M01302583+6121267STARsdss_apogee_apo25m_128+00_2m01302583+612126722.60765161.35743spectrumSDSS Collaboration355927.0503363657456234.1844604513912000.01514.01696.0Apache Point Observatory Galactic Evolution Experiment (APOGEE)59554.0N/A----CIRCLE 22.607651 61.35743 0.0002777777777777778--mast:SDSS/apogee/apo25m/128+00/2M01302583+6121267/apStar-dr17-2M01302583+6121267.fitsPUBLICFalsenan234163418674578448
scienceSDSSAPOGEEapogee-apo25mSDSSNoneINFRARED2M00230856+5909407STARsdss_apogee_apo25m_120-04_2m00230856+59094075.78567459.161308spectrumSDSS Collaboration355850.19970814814655876.1034505208312000.01514.01696.0Apache Point Observatory Galactic Evolution Experiment (APOGEE)59554.0N/A----CIRCLE 5.785674 59.161308 0.0002777777777777778--mast:SDSS/apogee/apo25m/120-04/2M00230856+5909407/apStar-dr17-2M00230856+5909407.fitsPUBLICFalsenan233850987673171837
scienceSDSSAPOGEEapogee-lco25mSDSSNoneINFRARED2M06532565-5934359STARsdss_apogee_lco25m_270-23-o_2m06532565-5934359103.356886-59.576645spectrumSDSS Collaboration358831.288909687558858.2394033101916000.01514.01696.0Apache Point Observatory Galactic Evolution Experiment (APOGEE)59554.0N/A----CIRCLE 103.356886 -59.576645 0.00018194444444444445--mast:SDSS/apogee/lco25m/270-23-O/2M06532565-5934359/asStar-dr17-2M06532565-5934359.fitsPUBLICFalsenan234508082675834494
scienceSDSSAPOGEEapogee-lco25mSDSSNoneINFRARED2M10374795-6005195STARsdss_apogee_lco25m_287-01-c_2m10374795-6005195159.449817-60.08876spectrumSDSS Collaboration358913.99908027777659187.2700110300948000.01514.01696.0Apache Point Observatory Galactic Evolution Experiment (APOGEE)59554.0N/A----CIRCLE 159.449817 -60.08876 0.00018194444444444445--mast:SDSS/apogee/lco25m/287-01-C/2M10374795-6005195/asStar-dr17-2M10374795-6005195.fitsPUBLICFalsenan234568700676053921
scienceSDSSAPOGEEapogee-apo25mSDSSNoneINFRARED2M03032813+0107147STARsdss_apogee_apo25m_178-48_mga_2m03032813+010714745.867221.120754spectrumSDSS Collaboration357779.1181316782457783.09716141203516000.01514.01696.0Apache Point Observatory Galactic Evolution Experiment (APOGEE)59554.0N/A----CIRCLE 45.86722 1.120754 0.0002777777777777778--mast:SDSS/apogee/apo25m/178-48_MGA/2M03032813+0107147/apStar-dr17-2M03032813+0107147.fitsPUBLICFalsenan234193034674700500
scienceSDSSAPOGEEapogee-lco25mSDSSNoneINFRARED2M06224728-7539285STARsdss_apogee_lco25m_lmc16_2m06224728-753928595.697025-75.657944spectrumSDSS Collaboration358096.2277165509358473.1985585879736000.01514.01696.0Apache Point Observatory Galactic Evolution Experiment (APOGEE)59554.0N/A----CIRCLE 95.697025 -75.657944 0.00018194444444444445--mast:SDSS/apogee/lco25m/LMC16/2M06224728-7539285/asStar-dr17-2M06224728-7539285.fitsPUBLICFalsenan234493306675790568
scienceSDSSAPOGEEapogee-apo25mSDSSNoneINFRARED2M00381127+5351499STARsdss_apogee_apo25m_120-08_2m00381127+53514999.54698753.863873spectrumSDSS Collaboration356226.21209560185556254.1438172222212000.01514.01696.0Apache Point Observatory Galactic Evolution Experiment (APOGEE)59554.0N/A----CIRCLE 9.546987 53.863873 0.0002777777777777778--mast:SDSS/apogee/apo25m/120-08/2M00381127+5351499/apStar-dr17-2M00381127+5351499.fitsPUBLICFalsenan234132807674410076

Searching for a specific star#

Let’s narrow down the search to one particular star: cepheid variable “V1154-Cyg”. Cepheid Variables are pulsating stars, which physically grow larger and smaller, and as a result brighter and dimmer, on a regular rhythm of few days. Cepheid variables in particular are useful to study in astronomy because their pulsation period correlates with their intrinsic brightness, making it easy to determine how far away the star is! Later in this notebook, we will look at some other kinds of variable stars.

We can search for this star in particular using the target_name keyword:

# Search for APOGEE star V1154_Cyg
apogee_obs_list = Observations.query_criteria(provenance_name='APOGEE',
                                              target_name='V1154_Cyg')

# Display results
apogee_obs_list
Table masked=True length=1
intentTypeobs_collectionprovenance_nameinstrument_nameprojectfilterswavelength_regiontarget_nametarget_classificationobs_ids_ras_decdataproduct_typeproposal_picalib_levelt_mint_maxt_exptimeem_minem_maxobs_titlet_obs_releaseproposal_idproposal_typesequence_numbers_regionjpegURLdataURLdataRightsmtFlagsrcDenobsidobjID
str7str4str6str12str4str4str8str9str4str35float64float64str8str18int64float64float64float64float64float64str63float64str3str1int64str55str1str67str6boolfloat64str9str9
scienceSDSSAPOGEEapogee-apo1mSDSSNoneINFRAREDV1154_CygSTARsdss_apogee_apo1m_cepheid_v1154_cyg297.0643920943.1268806458spectrumSDSS Collaboration358002.32424549768658069.12407793981424000.01514.01696.0Apache Point Observatory Galactic Evolution Experiment (APOGEE)59554.0N/A----CIRCLE 297.06439209 43.1268806458 0.0002777777777777778--mast:SDSS/apogee/apo1m/cepheid/V1154_Cyg/apStar-dr17-V1154_Cyg.fitsPUBLICFalsenan233753893675403444

From this results table, we can see some basic metadata related to this observation:

  • This star was observed with the APO 1-M telescope (instrument_name)

  • It’s coordinates are in the s_ra and s_dec columns

  • From the t_min column, we can see that this star was first observed on the date of MJD 58002 (Correpsonding to 2017-09-06) and last observed on MJD 58069 (2017-11-12 )

  • APOGEE provides infrared-wavelength (wavelength_region) spectra (dataproduct_type) with wavelength range of 1514 - 1696 nanometers (em_min, em_max)

Downloading APOGEE data products#

List all of the data products available for this star using Observations.get_product_list().

There are 8 total files available for this star, which include the individual spectra from each visit (APVISIT; apVisit-dr17-*-V1154_Cyg.fits), the combined spectrum (APSTAR; apStar-dr17-V1154_Cyg.fits), the processed and anaylzed spectrum (ASPCAPSTAR; aspcapStar-dr17-V1154_Cyg.fits) which contains the best-fit model and stellar parameters from the APOGEE Stellar Parameters and Chemical Abundance Pipeline (ASPCAP). Only the apStar and aspcapStar files are tagged as “Minimum Recommended Products’

More information on the APOGEE data products available at MAST can be found on the APOGEE Data Products in the Archive Manual, and more information on all of these products can be seen in the search results table:

# List all products available for this observation
products = Observations.get_product_list(apogee_obs_list)

# Show table
products
Table masked=True length=10
obsIDobs_collectiondataproduct_typeobs_iddescriptiontypedataURIproductTypeproductGroupDescriptionproductSubGroupDescriptionproductDocumentationURLprojectprvversionproposal_idproductFilenamesizeparent_obsiddataRightscalib_levelfilters
str9str4str12str35str209str1str74str7str28str15str1str6str4str3str33int64str9str6int64str4
233753893SDSSspectrumsdss_apogee_apo1m_cepheid_v1154_cygSDSS APOGEE single-visit spectrum from one night of observation.Smast:SDSS/apogee/apo1m/cepheid/V1154_Cyg/apVisit-dr17-58002-V1154_Cyg.fitsSCIENCE--APVISIT--APOGEEdr17N/AapVisit-dr17-58002-V1154_Cyg.fits498240233753893PUBLIC3None
233753893SDSSspectrumsdss_apogee_apo1m_cepheid_v1154_cygSDSS APOGEE single-visit spectrum from one night of observation.Smast:SDSS/apogee/apo1m/cepheid/V1154_Cyg/apVisit-dr17-58004-V1154_Cyg.fitsSCIENCE--APVISIT--APOGEEdr17N/AapVisit-dr17-58004-V1154_Cyg.fits498240233753893PUBLIC3None
233753893SDSSspectrumsdss_apogee_apo1m_cepheid_v1154_cygSDSS APOGEE single-visit spectrum from one night of observation.Smast:SDSS/apogee/apo1m/cepheid/V1154_Cyg/apVisit-dr17-58013-V1154_Cyg.fitsSCIENCE--APVISIT--APOGEEdr17N/AapVisit-dr17-58013-V1154_Cyg.fits498240233753893PUBLIC3None
233753893SDSSspectrumsdss_apogee_apo1m_cepheid_v1154_cygSDSS APOGEE single-visit spectrum from one night of observation.Smast:SDSS/apogee/apo1m/cepheid/V1154_Cyg/apVisit-dr17-58041-V1154_Cyg.fitsSCIENCE--APVISIT--APOGEEdr17N/AapVisit-dr17-58041-V1154_Cyg.fits498240233753893PUBLIC3None
233753893SDSSspectrumsdss_apogee_apo1m_cepheid_v1154_cygSDSS APOGEE single-visit spectrum from one night of observation.Smast:SDSS/apogee/apo1m/cepheid/V1154_Cyg/apVisit-dr17-58050-V1154_Cyg.fitsSCIENCE--APVISIT--APOGEEdr17N/AapVisit-dr17-58050-V1154_Cyg.fits498240233753893PUBLIC3None
233753893SDSSspectrumsdss_apogee_apo1m_cepheid_v1154_cygSDSS APOGEE single-visit spectrum from one night of observation.Smast:SDSS/apogee/apo1m/cepheid/V1154_Cyg/apVisit-dr17-58069-V1154_Cyg.fitsSCIENCE--APVISIT--APOGEEdr17N/AapVisit-dr17-58069-V1154_Cyg.fits498240233753893PUBLIC3None
233753893SDSSspectrumsdss_apogee_apo1m_cepheid_v1154_cygSDSS APOGEE combined, RV-corrected spectrum from all nights of observation for a single star.Smast:SDSS/apogee/apo1m/cepheid/V1154_Cyg/apStar-dr17-V1154_Cyg.fitsSCIENCEMinimum Recommended ProductsAPSTAR--APOGEEdr17N/AapStar-dr17-V1154_Cyg.fits4057920233753893PUBLIC3None
233753893SDSSspectrumsdss_apogee_apo1m_cepheid_v1154_cygSDSS APOGEE post-analysis pipeline (ASPCAP) spectrum. The combined spectrum is pseudo-continuum normalized, and includes the measured parameter estimates and best matching synthetic ASPCAP spectrum.Smast:SDSS/apogee/apo1m/cepheid/V1154_Cyg/aspcapStar-dr17-V1154_Cyg.fitsSCIENCEMinimum Recommended ProductsASPCAPSTAR--APOGEEdr17N/AaspcapStar-dr17-V1154_Cyg.fits247680233753893PUBLIC3None
234628801SDSSmeasurementssdss_apogee_allstarSDSS APOGEE allStar catalog. This summary table contains stellar parameters, radial velocities, chemical compositions and observation metadata for each star observed in APOGEE, combined across multiple visits.Dmast:SDSS/apogee/allStar-dr17-synspec_rev1.fitsSCIENCE--APOGEE Catalogs--APOGEEdr17N/AallStar-dr17-synspec_rev1.fits3970621440233753893PUBLIC3--
234630899SDSSmeasurementssdss_apogee_allvisitSDSS APOGEE allVisit catalog. This summary table contains target information, photometry, and observation summaries for each individual visit observed in APOGEE.Dmast:SDSS/apogee/allVisit-dr17-synspec_rev1.fitsSCIENCE--APOGEE Catalogs--APOGEEdr17N/AallVisit-dr17-synspec_rev1.fits2800627200233753893PUBLIC3--

Now we will download the aspcapStar spectrum for this star using Observations.download_products(). The download will print a status message when completed.

manifest = Observations.download_products(products, productSubGroupDescription='ASPCAPSTAR')

Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:SDSS/apogee/apo1m/cepheid/V1154_Cyg/aspcapStar-dr17-V1154_Cyg.fits to ./mastDownload/SDSS/sdss_apogee_apo1m_cepheid_v1154_cyg/aspcapStar-dr17-V1154_Cyg.fits ...

 [Done]

Plotting an APOGEE Spectrum#

Now let’s take a look at the file and plot the spectrum.

Based on the descriptions in the aspcapStar Data Model, this file has four extensions:

  • HDU0: The Primary Header and file metadata

  • HDU1: The array containing the observed combined spectrum

  • HDU2: The error array for the spectrum

  • HDU3: The best fit model spectrum from ASPCAP

  • HDU4: The ASPCAP Data Table, containing the best fit stellar parameter values and other information for this star

# Open file
aspcap_star = fits.open(manifest['Local Path'][0])

# Display file information
aspcap_star.info()

Filename: ./mastDownload/SDSS/sdss_apogee_apo1m_cepheid_v1154_cyg/aspcapStar-dr17-V1154_Cyg.fits
No.    Name      Ver    Type      Cards   Dimensions   Format
  0  PRIMARY       1 PrimaryHDU       5   ()      
  1                1 ImageHDU        11   (8575,)   float64   
  2                1 ImageHDU        11   (8575,)   float64   
  3                1 ImageHDU        11   (8575,)   float64   
  4                1 BinTableHDU    270   1R x 123C   [64A, 30A, 6A, J, 20A, 30A, D, D, D, D, E, E, E, E, E, E, 16A, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, J, J, E, E, 16A, E, 32A, E, E, J, J, J, J, J, J, 132A, 32A, 32A, J, J, J, J, E, E, K, 132A, K, 132A, E, E, E, E, E, E, E, E, E, E, E, J, J, E, E, E, E, E, E, K, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, E, 8A, 189E, 21E, 9E, 81E, E, 9E, 81E, 9K, K, 256A, E, E, E, 27E, 27D, 27E, 27E, 27E, 27E, 27E, 27E, 27K]

We can plot the spectrum using this information! The wavelength information can be found in the file hearer (HUD0), and observed flux is contained in the first extension (HDU1).

plt.figure(figsize=(15, 5))

# Get wavelength information from the header
# CRVAL1 is the minimum wavlength, CDELT1 is the step size (in log space) and NAXIS1 is the number of pixels
apogee_wls = np.logspace(aspcap_star[1].header['CRVAL1'],
                         aspcap_star[1].header['CRVAL1']+aspcap_star[1].header['CDELT1']*aspcap_star[3].header['NAXIS1'],
                         aspcap_star[1].header['NAXIS1'])


# Get the observed and model flux from extensions 1 and 3
observed_flux = aspcap_star[1].data
model_flux = aspcap_star[3].data

# Mask the data using the error array in extenion 2
pixel_mask = (aspcap_star[2].data < 0.1) # use only pixels with small errors (10% or less)

# Plot the observed spectrum
plt.plot(apogee_wls[pixel_mask], observed_flux[pixel_mask], c='k', label='Observed Spectrum')

# Plot the model spectrum
plt.plot(apogee_wls[pixel_mask], model_flux[pixel_mask], c='darkred', label='Best Fit Model Spectrum')

# Set axes limits
plt.ylim(0.6, 1.1)
plt.xlim(np.min(apogee_wls), np.max(apogee_wls))

# Set labels
apogee_id = aspcap_star[4].data['APOGEE_ID'][0] # String containing star's name
plt.title(f"APOGEE Spectrum - {apogee_id}")
plt.xlabel(r'Wavelength [$\AA$]')
plt.ylabel('Flux')
plt.grid()
plt.legend()

plt.show()
../../../_images/e3804fdec567d6284654b451c3d7178432f5f6ad4ada73e1ec05bccfc69e9361.png

This is the APOGEE spectrum for this star! It looks great. We can tell it must be a hot star because of the large Hydrogen Absorption lines (for example around 15346, 15443, 15560 and 16113 angstroms - this is the Hydrogen Brackett Series!) The gaps in the spectrum around 15800 and 16400 Angstroms are due to the APOGEE instrumentation: APOGEE observes across 3 ccds, which have a small gap in wavelength coverage between them.

Downloading the APOGEE allStar Catalog#

We can also download the APOGEE allStar catalog, which contains a lot of useful information for the full APOGEE sample, including stellar parameters, chemical abundances, positions, and observation data for all of the stars in APOGEE.

We can do this with the download_file function in astroquery, knowing the MAST data URL for the allStar catalog is mast:SDSS/apogee/allStar-dr17-synspec_rev1.fits.

This is a fairly large file (3.9 GB), so this download may take a few minutes!

# Download allStar file - this may take a few minutes!
Observations.download_file('mast:SDSS/apogee/allStar-dr17-synspec_rev1.fits')

Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:SDSS/apogee/allStar-dr17-synspec_rev1.fits to allStar-dr17-synspec_rev1.fits ...

 [Done]

('COMPLETE', None, None)

# Open file 
allstar = Table.read('allStar-dr17-synspec_rev1.fits', hdu=1)

# Display table
allstar[:10]
Table length=10
FILEAPOGEE_IDTARGET_IDAPSTAR_IDASPCAP_IDTELESCOPELOCATION_IDFIELDALT_IDRADECGLONGLATJJ_ERRHH_ERRKK_ERRSRC_HWASH_MWASH_M_ERRWASH_T2WASH_T2_ERRDDO51DDO51_ERRIRAC_3_6IRAC_3_6_ERRIRAC_4_5IRAC_4_5_ERRIRAC_5_8IRAC_5_8_ERRIRAC_8_0IRAC_8_0_ERRWISE_4_5WISE_4_5_ERRTARG_4_5TARG_4_5_ERRWASH_DDO51_GIANT_FLAGWASH_DDO51_STAR_FLAGTARG_PMRATARG_PMDECTARG_PM_SRCAK_TARGAK_TARG_METHODAK_WISESFD_EBVAPOGEE_TARGET1APOGEE_TARGET2APOGEE2_TARGET1APOGEE2_TARGET2APOGEE2_TARGET3APOGEE2_TARGET4TARGFLAGSSURVEYPROGRAMNAMENVISITSSNRSNREVSTARFLAGSTARFLAGSANDFLAGANDFLAGSVHELIO_AVGVSCATTERVERRRV_TEFFRV_LOGGRV_FEHRV_ALPHARV_CARBRV_CHI2RV_CCFWHMRV_AUTOFWHMRV_FLAGN_COMPONENTSMEANFIBSIGFIBMIN_HMAX_HMIN_JKMAX_JKGAIAEDR3_SOURCE_IDGAIAEDR3_PARALLAXGAIAEDR3_PARALLAX_ERRORGAIAEDR3_PMRAGAIAEDR3_PMRA_ERRORGAIAEDR3_PMDECGAIAEDR3_PMDEC_ERRORGAIAEDR3_PHOT_G_MEAN_MAGGAIAEDR3_PHOT_BP_MEAN_MAGGAIAEDR3_PHOT_RP_MEAN_MAGGAIAEDR3_DR2_RADIAL_VELOCITYGAIAEDR3_DR2_RADIAL_VELOCITY_ERRORGAIAEDR3_R_MED_GEOGAIAEDR3_R_LO_GEOGAIAEDR3_R_HI_GEOGAIAEDR3_R_MED_PHOTOGEOGAIAEDR3_R_LO_PHOTOGEOGAIAEDR3_R_HI_PHOTOGEOASPCAP_GRIDFPARAM_GRIDCHI2_GRIDFPARAMFPARAM_COVASPCAP_CHI2PARAMPARAM_COVPARAMFLAGASPCAPFLAGASPCAPFLAGSFRAC_BADPIXFRAC_LOWSNRFRAC_SIGSKYFELEMFELEM_ERRX_HX_H_ERRX_MX_M_ERRELEM_CHI2ELEMFRACELEMFLAGEXTRATARGMEMBERFLAGMEMBERX_H_SPECX_M_SPECTEFFTEFF_ERRLOGGLOGG_ERRM_HM_H_ERRALPHA_MALPHA_M_ERRVMICROVMACROVSINITEFF_SPECLOGG_SPECC_FEC_FE_SPECC_FE_ERRC_FE_FLAGCI_FECI_FE_SPECCI_FE_ERRCI_FE_FLAGN_FEN_FE_SPECN_FE_ERRN_FE_FLAGO_FEO_FE_SPECO_FE_ERRO_FE_FLAGNA_FENA_FE_SPECNA_FE_ERRNA_FE_FLAGMG_FEMG_FE_SPECMG_FE_ERRMG_FE_FLAGAL_FEAL_FE_SPECAL_FE_ERRAL_FE_FLAGSI_FESI_FE_SPECSI_FE_ERRSI_FE_FLAGP_FEP_FE_SPECP_FE_ERRP_FE_FLAGS_FES_FE_SPECS_FE_ERRS_FE_FLAGK_FEK_FE_SPECK_FE_ERRK_FE_FLAGCA_FECA_FE_SPECCA_FE_ERRCA_FE_FLAGTI_FETI_FE_SPECTI_FE_ERRTI_FE_FLAGTIII_FETIII_FE_SPECTIII_FE_ERRTIII_FE_FLAGV_FEV_FE_SPECV_FE_ERRV_FE_FLAGCR_FECR_FE_SPECCR_FE_ERRCR_FE_FLAGMN_FEMN_FE_SPECMN_FE_ERRMN_FE_FLAGFE_HFE_H_SPECFE_H_ERRFE_H_FLAGCO_FECO_FE_SPECCO_FE_ERRCO_FE_FLAGNI_FENI_FE_SPECNI_FE_ERRNI_FE_FLAGCU_FECU_FE_SPECCU_FE_ERRCU_FE_FLAGCE_FECE_FE_SPECCE_FE_ERRCE_FE_FLAGYB_FEYB_FE_SPECYB_FE_ERRYB_FE_FLAGVISIT_PK
bytes64bytes30bytes58bytes71bytes77bytes6int32bytes20bytes30float64float64float64float64float32float32float32float32float32float32bytes16float32float32float32float32float32float32float32float32float32float32float32float32float32float32float32float32float32float32int32int32float32float32bytes16float32bytes32float32float32int32int32int32int32int32int32bytes132bytes32bytes32int32float32float32int64bytes132int64bytes132float32float32float32float32float32float32float32float32float32float32float32int32int32float32float32float32float32float32float32int64float32float32float32float32float32float32float32float32float32float32float32float32float32float32float32float32float32bytes8float32[21,9]float32[21]float32[9]float32[9,9]float32float32[9]float32[9,9]int64[9]int64bytes256float32float32float32float32[27]float64[27]float32[27]float32[27]float32[27]float32[27]float32[27]float32[27]int64[27]int32int64bytes10float32[27]float32[27]float32float32float32float32float32float32float32float32float32float32float32float32float32float32float32float32int32float32float32float32int32float32float32float32int32float32float32float32int32float32float32float32int32float32float32float32int32float32float32float32int32float32float32float32int32float32float32float32int32float32float32float32int32float32float32float32int32float32float32float32int32float32float32float32int32float32float32float32int32float32float32float32int32float32float32float32int32float32float32float32int32float32float32float32int32float32float32float32int32float32float32float32int32float32float32float32int32float32float32float32int32float32float32float32int32int32[100]
apStar-dr17-VESTA.fitsVESTAapo1m.calibration.VESTAapogee.apo1m.stars.calibration.VESTAapogee.apo1m.synspec_fix.calibration.VESTAapo1m1calibration----292.21913091040517-30.6029193636988199.9990.099.9990.099.9990.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.0000.00.0------041943040419430400APOGEE_1MTARGETAPOGEE2_1M_TARGETapo1m2418.9276441.157440--0--9.25731210.1855150.0240308985980.21344.643887-0.203651370.00.016.5416925.0923920.430801223.000020.0-9999.999999.99-9999.999999.990----------------------------------GKd_b-- .. ---- .. --5766.8 .. 0.0222.64 .. 0.02.5654295694.8984 .. 0.0644.19434 .. 1.0256 .. 0133TEFF_WARN,VMICRO_WARN,STAR_WARN0.00226244350.0259515570.0228906050.002652 .. -0.444710.007184497080743313 .. 0.0117932185530662540.0083094 .. --0.029537989 .. 0.156177790.002652 .. --0.02910603 .. 0.156177791.7949402 .. 0.97651520.9752284 .. 1.00 .. 25890--0.0083094 .. --0.002652 .. --5694.898425.3809834.42851970.0199736320.00565740.00503307250.00667910.00651770040.412648230.04.77892265766.84.38320.00484679970.00484679970.02953798900.01127990.01127990.02859264600.124264810.124264810.0230221500.11493770.11024480.0493361200.1466680.06820440.1113361900.0351473-0.03214120.01149605400.0752214040.06279540.01362693100.0671891-0.00595390.0118011360----0.085440852-0.0894567-0.0623892020.041454110-0.0944697-0.18929260.0326718870-0.01682920.05245480.0108100160-0.23147619-0.184545190.0561924840------640.094247950.12911740.0960142540-0.0257743020.13676740.04869532600.04086991-0.11181260.01146729800.00346260.00346260.00540092170------640.0512783970.04892840.0102655040----0.0594029722------64----1.020 .. -2147483648
apStar-dr17-2M00000002+7417074.fits2M00000002+7417074apo25m.120+12.2M00000002+7417074apogee.apo25m.stars.120+12.2M00000002+7417074apogee.apo25m.synspec_fix.120+12.2M00000002+7417074apo25m5046120+12none0.00010374.285408119.4018072243182211.7674135414862428.5970.0397.6670.0297.3140.0182MASS----------------------------7.3530.027.3530.02-11-2.82.4UCAC40.24235167RJCE_WISE0.242351670.43400200-2147465196000APOGEE2_TWOBIN_GT_0_8,APOGEE2_WISE_DERED,APOGEE2_SHORT,APOGEE2_NORMAL_SAMPLEapogee2disk3827.1565669.41653072PERSIST_MED,PERSIST_LOW0---51.7314530.204533160.0205957353671.80151.0904876-0.40140440.00.082.418517.40865912.7905980123.95515330.171737.012.20.89999.995380282167077157120.285441670.0194750.048581990.02403150.52224850.022098111.698813.272210.4956-51.92470.3656463234.57473027.02953464.55643268.6013115.6933507.7324Mg_d-- .. ---- .. --3621.7 .. 0.00.64348 .. 0.015.2932583723.9111 .. 0.016.48606 .. 1.00 .. 00--0.00.00.0-0.0037053 .. 0.224820.00035680527798831463 .. 0.0001349740632576868-0.15138529 .. --0.008765366 .. 0.023297967-0.0037053 .. --0.005022976 .. 0.02329796710.580349 .. 11.9453811.0 .. 1.00 .. 200---0.15138529 .. ---0.0037053 .. --3723.91114.06030270.90459830.021523084-0.147680.00718340840.0369220.00473469261.71218233.223484--3621.71.02450.0092947040.0092947040.0087653660-0.0068259984-0.00682599840.0176710600.151220.151220.01019773700.083401990.0267390010.00820233500.0501119940.0307099970.03325760700.030429006-0.0330940040.010857150------64-0.028965399-0.0235660080.0121813950----0.086512072------640.0862190.0164600010.0330202360-0.098511994-0.080507010.0122871830------64------64-0.10815802-0.153120010.031854160------64------64-0.16068-0.160680.00925129700.075089990.0270199920.02017369700.0076825023-0.0357100070.0106090360----0.00122116722------64----1.022 .. -2147483648
apStar-dr17-2M00000019-1924498.fits2M00000019-1924498apo25m.060-75.2M00000019-1924498apogee.apo25m.stars.060-75.2M00000019-1924498apogee.apo25m.synspec_fix.060-75.2M00000019-1924498apo25m5071060-75none0.000832-19.41385163.39412165971994-75.9063967700245311.0740.02210.740.02610.670.0232MASS----------------------------10.6660.02-----1119.6-6.5UCAC40.006158138SFD0.0220314530.02039118700-2147399648000APOGEE2_SFD_DERED,APOGEE2_SHORT,APOGEE2_NORMAL_SAMPLE,APOGEE2_ONEBIN_GT_0_3apogee2halo3229.87506233.314590--0--19.0738620.087965140.0385987465854.6874.653328-0.436138960.00.05.79853724.57305119.39943701166.00.07.012.20.39999.9924139298125874590720.082092860.31508120.2900220.302588-10.341260.27387112.288212.514911.659717.88042.120882294.90011648.59553922.614888.34845797.92316970.81085GKd_b-- .. ---- .. --5522.4 .. 0.0236.7 .. 0.02.80995925501.773 .. 0.0652.9846 .. 1.00 .. 068VMICRO_WARN,N_M_WARN0.00013308490.0248868780.0247537920.054958 .. 0.424430.01083282008767128 .. 0.005686562974005938-0.21379201 .. --0.033176932 .. 0.147667810.054958 .. --0.032482434 .. 0.147667812.0648096 .. 1.10576030.97642267 .. 0.909090940 .. 25800---0.21379201 .. --0.054958 .. --5501.77325.5535644.3041150.024554765-0.268750.006752810.09097830.0073728020.409694370.06.27364545522.44.1950.0617380.0617380.03317693200.070630.070630.03355040-0.1119-0.11190.0326301052560.235342890.230650.0462971140-0.19082639-0.269290.1457230900.16523850.097950.013435200.296303720.28387770.01554261300.096592990.0234499870.0134200520----0.0927759820.09333250.12040.0511001700.1155729140.0207500160.0355027600.097265990.166550.01375008200.112659010.159590.0587675460------640.0505005570.0853700040.101880510-0.130981710.0315600040.0604799840-0.07802749-0.230710.0171355290-0.27553-0.275530.00725272670------640.0139300230.011580020.01383546650----0.1025938242------64----1.025 .. -2147483648
apStar-dr17-2M00000032+5737103.fits2M00000032+5737103apo25m.116-04.2M00000032+5737103apogee.apo25m.stars.116-04.2M00000032+5737103apogee.apo25m.synspec_fix.116-04.2M00000032+5737103apo25m4424116-04none0.00133557.61953116.06537094825998-4.56476751407882110.9050.02310.6350.0310.4830.0222MASS----------------------------10.4470.0210.430.0-110.00.0NOMAD0.14228994RJCE_WISE_PARTSKY0.12668462---2147483584-21474831360000APOGEE_NO_DERED,APOGEE_TELLURICapogeeapogee3121.3259152.67575512PERSIST_HIGH512PERSIST_HIGH-20.5451640.130395950.0653062766399.0024.2702193-0.70242210.00.03.577794628.49852616.8539640124.6530926.664583-9999.999999.99-9999.999999.994225966799645137921.29853120.00952245-0.410274060.00872931-3.22660140.0089724312.220912.6511.6189-19.1961.35642762.9458758.2913767.3334761.33756.7403766.5322Fd_d-- .. 0.0-- .. 1.92779126210.7 .. 0.0568.35 .. 0.01.92779126099.781 .. 0.02101.9775 .. 1.00 .. 00--0.00.0260846430.0258184730.10486 .. 0.499420.011564601212739944 .. 0.009124253876507282-0.14024001 .. --0.06806532 .. 0.26746480.10486 .. --0.06770368 .. 0.26746481.7376808 .. 1.73600110.97735685 .. 1.00 .. 259210---0.14024001 .. --0.10486 .. --6099.78145.847333.67397050.030661006-0.24510.0070071280.02915930.01175542552.14343330.014.9658016210.73.52590.112730.112730.0680653200.1244400140.1244400140.0466395600.114560010.114560010.3859263320.0450689050.0403760080.145615770-0.17934635-0.257809970.253778900.038494498-0.028793990.01549950500.0595960170.0471700130.01845200700.1382180.065075010.0170270520----0.176573142-0.09011951-0.0630520.0596198400.14228290.0474600050.0550385940-0.12913498-0.059850990.01790468600.170519010.217450010.116481510------64-0.00546944140.0294000060.222812950-1.1124718-0.94993010.1023763900.06363252-0.0890499950.0190900150-0.25297-0.252970.0076459660------64-0.108749986-0.111099990.0165547320----0.123838582------64----1.028 .. -2147483648
apStar-dr17-2M00000032+5737103.fits2M00000032+5737103apo25m.N7789.2M00000032+5737103apogee.apo25m.stars.N7789.2M00000032+5737103apogee.apo25m.synspec_fix.N7789.2M00000032+5737103apo25m4264N7789none0.00133557.61953116.06537094825998-4.56476751407882110.9050.02310.6350.0310.4830.0222MASS----------------------------10.4470.0210.430.0-110.00.0NOMAD0.10556992RJCE_WISE_PARTSKY0.12668462---214748160000000APOGEE_SHORTapogeeapogee3221.82442218.282410--0---20.434650.163896650.062999316401.96344.253023-0.72442070.00.04.1942330.61907817.23286401263.637337.50555375.012.50.09999.994225966799645137921.29853120.00952245-0.410274060.00872931-3.22660140.0089724312.220912.6511.6189-19.1961.35642762.9458758.2913767.3334761.33756.7403766.5322Fd_a-- .. ---- .. --6284.1 .. 0.0409.34 .. 0.02.56974486162.0303 .. 0.01743.381 .. 1.00 .. 00--0.00.0239552840.0234229430.039991 .. 0.136470.010498571209609509 .. 0.009732984937727451-0.181519 .. --0.058728103 .. 0.251368940.039991 .. --0.058453497 .. 0.251368942.2316186 .. 2.14501330.976819 .. 1.00 .. 200---0.181519 .. --0.039991 .. --6162.030341.7538153.71556070.025935415-0.221510.0056726380.04806630.00985560751.95789780.014.9933956284.13.55850.0326509920.0326509920.05872810300.070725010.070725010.03793010-0.48906004-0.489060040.67714272880.14057290.135880010.12401360-1.1662664-1.244730.214010540-0.0010954887-0.068384010.01289359800.0130859910.000659987330.0139255800.1300760.0569330.0135523750----0.1431851220.12813250.15520.0475980160-0.00018709898-0.095010.046923150-0.170857-0.1015730050.0139221160------1------64-0.059629455-0.0247600080.17689870-0.7412817-0.578740.082906820-0.15302749-0.305710020.0150238220-0.21417-0.214170.0061514160------64-0.06943999-0.0717899950.0126378440----0.107103192------64----1.028 .. -2147483648
asStar-dr17-2M00000035-7323394.fits2M00000035-7323394lco25m.SMC12.2M00000035-7323394apogee.lco25m.stars.SMC12.2M00000035-7323394apogee.lco25m.synspec_fix.SMC12.2M00000035-7323394lco25m7218SMC12none0.001467-73.394287307.93944096852624-43.2303046496328215.0080.04514.3520.05814.2770.0682MASS----------------------------14.1010.035-----110.2-1.1GAIA-DR20.0109211225SFD0.184518370.0361626600-2139095040-214741811200APOGEE2_MAGCLOUD_CANDIDATE,APOGEE2_GAIA_OVERLAPapogee2smagclouds335.32554636.7778470--0--165.671940.111918340.188125154667.1592.1917586-1.35303750.00.01.918990119.57644714.38868101186.957373.4641016-9999.999999.99-9999.999999.994689447878791422208-0.023941280.04500720.255408050.0535888-1.26917550.0557216.747817.386115.9896----12516.22958956.69917395.0377310.75155323.14211760.366GKg_a-- .. ---- .. --4378.4 .. 0.0451.09 .. 0.02.41318154555.4043 .. 0.0451.09 .. 1.032 .. 00--0.0058557360.105003990.099813685-0.27891 .. -1.49730.01720406860113144 .. 0.21653637290000916-1.46211 .. --0.05615539 .. 0.21653637-0.27891 .. --0.052546624 .. 0.216536372.9201825 .. 0.594524860.89734745 .. 1.00 .. 3510---1.46211 .. ---0.27891 .. --4555.404321.238881.49885110.08263534-1.18320.019806060.0084420010.0259100371.38140185.9081445--4378.41.4375-0.2907101-0.29071010.056155390------10.47552990.47552990.0537535300.152532940.09586990.051480500.548311950.52890990.3106469800.0068899393-0.0566329960.0376912240-0.4556191-0.314300060.06000968400.037651540.0430510040.043763850----0.334372072-0.38900805-0.412790060.209710170-0.92324126-0.993000150.110225220-0.33387506-0.315870050.061497510-0.31694686-0.31652010.065656270------10.186761920.141799930.2647038700.64053390.638119940.19102240-0.15649807-0.273400070.0604968440-1.1714-1.17140.02140348400.020969868-0.0271000860.257513170-0.03820753-0.081600070.053996870----0.27679532-0.15516007-0.099500060.146836860----1.0214 .. -2147483648
apStar-dr17-2M00000068+5710233.fits2M00000068+5710233apo25m.N7789.2M00000068+5710233apogee.apo25m.stars.N7789.2M00000068+5710233apogee.apo25m.synspec_fix.N7789.2M00000068+5710233apo25m4264N7789none0.0028557.173164115.97715352490295-5.00239236972543710.6640.02310.1320.0310.0180.022MASS----------------------------10.0320.0210.0360.0-118.0-18.0NOMAD0.0055077374RJCE_WISE_PARTSKY0.045900350.52296805-214748160000000APOGEE_SHORTapogeeapogee3294.87552282.745940--0---12.6737870.12024520.029277045315.28273.8085988-0.287583740.00.09.71490321.64597517.97247101263.851874.5825765.012.50.09999.994210775972675511041.38635560.01756725.79172470.0161734-12.7355890.015647312.275912.831111.5693-12.73590.843714718.71344709.5037727.03937715.0294706.5003723.68005GKg_a-- .. ---- .. --4991.4 .. 0.050.409 .. 0.04.35080865031.2637 .. 0.0102.92217 .. 1.00 .. 00--0.00.0167686980.016768698-0.025931 .. 0.143310.002885030349716544 .. 0.0019492049468681216-0.186011 .. --0.015217459 .. 0.084045514-0.025931 .. --0.013979203 .. 0.0840455142.8813765 .. 1.5712670.9859233 .. 1.00 .. 200---0.186011 .. ---0.025931 .. --5031.263710.1450573.4561320.023554068-0.160080.00601273030.0557040.0060627560.67197023.2469554--4991.43.4848-0.025730997-0.0257309970.0152174590-0.048318997-0.0483189970.02504017200.0841760040.0841760040.01859093600.1136980060.0570350.02343307200.1164910050.097089010.07013439400.092060010.0285370050.01130474600.0141167940.15543580.0183022200.0148876010.0202870070.0136431160----0.1443554920.02548790.00170589980.034081970-0.027150989-0.0969099850.03666611400.0382410060.0562460050.013792179500.0212841780.0217110070.0223428230-0.18249099-0.140739990.089618100.20249780.15753580.0946040600.0289040060.0264900030.042318540-0.029748008-0.146650.0145151680-0.16028-0.160280.00688997840-0.26304-0.311110020.0660599800.0590725020.015680.0116009210----0.0044014412-0.0555800058.0004334e-050.086572440----1.0226 .. -2147483648
apStar-dr17-2M00000103+1525513.fits2M00000103+1525513apo25m.107-46_MGA.2M00000103+1525513apogee.apo25m.stars.107-46_MGA.2M00000103+1525513apogee.apo25m.synspec_fix.107-46_MGA.2M00000103+1525513apo25m2897107-46_MGAnone0.00432215.430942105.06543952239785-45.6493480677228811.2780.02311.0060.02410.9470.0192MASS12.5353490.00847058711.6990410.01381888612.512720.0060115126----------------10.9370.0210.9370.0201-2.2-6.3URAT10.017441347RJCE_WISE0.0174413470.0508918900-2147383280000APOGEE2_WISE_DERED,APOGEE2_SHORT,APOGEE2_MANGA_LED,APOGEE2_ONEBIN_GT_0_3manga-apogee2manga2126.0876125.954730--0---12.0879240.073868450.0598753436144.2124.5084305-0.62963920.00.02.965346822.9536518.26930401222.999980.07.011.50.39999.9927720976194176087042.0844050.0180373-3.2428450.0171125-8.7514420.010061612.282112.581211.8179-13.931.53426473.27496469.28427477.08853473.27005468.36456477.71826Fd_b-- .. ---- .. --6032.1 .. 0.01158.2 .. 0.01.79110075945.751 .. 0.01957.7793 .. 1.00 .. 068VMICRO_WARN,N_M_WARN0.00359329260.0488421620.0440511030.055067 .. 0.497740.017340991646051407 .. 0.012145369313657284-0.21197301 .. --0.06668289 .. 0.251714650.055067 .. --0.06618328 .. 0.251714651.6491114 .. 1.53567740.9537347 .. 1.00 .. 25900---0.21197301 .. --0.055067 .. --5945.75144.24684.12319230.033373106-0.267040.0081474990.01240879950.0123255970.45948360.06.75227366032.14.13320.0499169830.0499169830.0666828900.089093980.089093980.0512538660-0.48443002-0.484430022.7495032880.0164178760.0117249790.13338370------10.046251476-0.0210370120.01732564200.1528260.140399990.02206321600.06266858-0.0104744140.0193041130----0.183835122-0.06725854-0.0401910250.0685782200.24792290.153099980.05728067500.051335990.120619980.02055560800.079758990.126689990.114116470------65-0.2820894-0.247219980.22983770-0.0741516950.088389990.1089736400.05873251-0.093950.0223818850-0.26189-0.261890.0089061940------64-0.037420005-0.0397700070.0196000470----0.133414032------64----1.0229 .. -2147483648
apStar-dr17-2M00000133+5721163.fits2M00000133+5721163apo25m.NGC7789_btx.2M00000133+5721163apogee.apo25m.stars.NGC7789_btx.2M00000133+5721163apogee.apo25m.synspec_fix.NGC7789_btx.2M00000133+5721163apo25m5922NGC7789_btxnone0.00555857.354549116.0147761731873-4.82491659463106113.2610.02912.6050.03612.3830.0252MASS----------------------------12.2850.02312.2850.023-11-2.3-1.7URAT10.2478597RJCE_WISE0.2478597--00-2147463150000APOGEE2_TWOBIN_0_5_TO_0_8,APOGEE2_WISE_DERED,APOGEE2_MEDIUM,APOGEE2_NORMAL_SAMPLEapogee2-mangacal_btx7116.9411120.2954850331652BRIGHT_NEIGHBOR,MTPFLUX_LT_75,MTPFLUX_LT_500---101.7134860.1242127640.066519525121.1662.971028-0.693380360.00.03.004735519.88148915.87781301209.9607446.42505612.212.80.50.84210863633054368000.127609940.0295789-1.1467650.02521950.364161520.02601515.55116.426714.6271----5883.3785068.65627047.50345923.9555205.1376869.6465GKg_b-- .. ---- .. --4897.1 .. 0.0129.53 .. 0.02.17355184974.601 .. 0.0234.4223 .. 1.00 .. 00--0.00.0266169810.025419217-0.065567 .. -0.36550.005662243347615004 .. 0.010711676441133022-0.57870704 .. --0.028721372 .. 0.13720466-0.065567 .. --0.02698478 .. 0.137204661.7831587 .. 3.15456870.97238356 .. 1.00 .. 200---0.57870704 .. ---0.065567 .. --4974.60115.3108562.38825060.035517264-0.513140.0098355930.107150.0106554521.6445993.991985--4897.12.4374-0.06917703-0.069177030.0287213720-0.12475002-0.124750020.0454574800.243079990.243079990.03148038300.185942980.129279970.0393129900.015072018-0.00432997940.139961600.133488980.069965990.01890180300.158510980.299830020.0312654700.092230590.0976299940.0230380330----0.2163851920.211851980.188070.0725127100.154339020.0845800040.0592240700.0919950.1099999850.02630023100.0795971450.0800239740.0382046660-0.088947-0.047196030.138037960-0.09750801-0.142470.166298910-0.20686603-0.209280010.0847979860-0.107198-0.22410.0267691260-0.50953-0.509530.01101191800.084710.036640020.1377790900.022392541-0.0209999680.0226617730----0.0251306762-0.25649-0.200829980.13238880----1.0231 .. -2147483648
apStar-dr17-2M00000211+6327470.fits2M00000211+6327470apo25m.117+01.2M00000211+6327470apogee.apo25m.stars.117+01.2M00000211+6327470apogee.apo25m.synspec_fix.117+01.2M00000211+6327470apo25m4591117+01none0.00880263.463078117.227985339705781.162166785239817611.880.02411.1070.0310.8610.0222MASS----------------------------10.7960.0210.7790.02-110.012.0NOMAD0.25520405RJCE_WISE_ALLSKY0.23959874---214748158400000APOGEE_WISE_DERED,APOGEE_SHORTapogeeapogee4111.05336169.020421536PERSIST_HIGH,PERSIST_MED0---106.631860.053484740.0290859794880.472.722735-0.16598170.00.06.985388819.06305715.88991450171.306621.266177.012.20.59999.994315949800534227200.069583740.0354776-0.923787950.0366333-0.90800610.035702214.38315.350413.3986----8216.6386704.69149909.0396924.2865814.23148930.6GKg_c-- .. ---- .. --4597.5 .. 0.022.344 .. 0.03.46282084681.4634 .. 0.076.84869 .. 1.00 .. 00--0.00013308490.0291455950.025286132-0.18688 .. -0.145790.0018862661672756076 .. 0.0017332916613668203-0.255647 .. --0.013141024 .. 0.059426643-0.18688 .. --0.011139332 .. 0.0594266432.789973 .. 1.407570.970987 .. 1.00 .. 200---0.255647 .. ---0.18688 .. --4681.46348.7663382.2455270.024334472-0.0687670.00697149760.00769400040.0060165881.48251813.078036--4597.52.3951-0.180636-0.1806360.0131410240-0.27499598-0.274995980.02485163700.3455940.3455940.0160572400.027782995-0.028880.0165598600.1255750.1061730.0603687100.012056999-0.0514660030.0120063140-0.024401010.1169180050.0221654330-0.014249399-0.0088500010.0146067750----0.1305191520.1055330.0817510.0384314320-0.07762-0.1473790.038437560-0.00176250190.0162424970.0146273820-0.065623835-0.0651970060.0211846510-0.15478702-0.113036010.07854650-0.40178698-0.446748970.0760388150-0.046804994-0.0492189970.03753106300.010702997-0.1061990.0153161230-0.075011-0.0750110.0083958290-0.164809-0.212878990.0548548330-0.025256492-0.0686489940.0127329750------2590.194040980.249700980.070006490----1.0238 .. -2147483648

Plotting an HR Diagram (or Kiel Diagram) with APOGEE#

When characterizing stars, a useful plot is the HR diagram, which helps classify stars by their different types (for example, a star can be a Sun-like main sequence star or a massive Red Giant star). Let’s plot an HR Diagram using APOGEE data.

Traditional HR diagrams plot a star’s color on the x-axis and its brightness (or magnitude) on the y-axis. For this excerise, we’ll be using the star’s temperature (\(T_{eff}\)) as a substitute for color and the surface gravity (\(\log g\)) as a proxy for brightness. This variation on the HR diagram is technically called a Kiel Diagram, but they are very similar in appearance and interpretation.

# Plot HR Diagram form APOGEE catalog
plt.figure(figsize=(10, 10))

# Plot Temperature (TEFF) against surface gravity (LOGG)
# And color points by their metallicity (FE_H)
im = plt.scatter(allstar['TEFF'], allstar['LOGG'], c=allstar['FE_H'],
                 marker='.', s=1, zorder=1,
                 cmap='jet', vmin=-1, vmax=0.5)
        
# Add colorbar legend to plot
plt.colorbar(im, location='bottom', label='Metallicity [Fe/H]')

        
# Labels and titles
plt.xlim(7000, 3000)
plt.ylim(5.5, -0.5)
plt.xlabel(r'Effective Temperature [$T_{\rm{eff}}$, K]')
plt.ylabel(r'Surface Gravity [$\log(g)$]')
plt.grid()
plt.title('APOGEE Stellar Parameters - allStar')

plt.show()
../../../_images/d80dd6a45a3e611ca25319d59219c4dc42261495d80779d5227a600236f80a02.png

This is an HR Diagram made with APOGEE data! The Main Sequence stars are along the bottom, and the Red Giant Branch is in the upper-right portion. The Horizontal Branch stars are the upper-left group of dark blue (metal poor) points.

Searching for TESS Data of this Star#

Now let’s search MAST for a light curve of our Cepheid Variable star.

Coordinate search using astroquery.mast#

# Retrieve RA and DEC of APOGEE star
ra = apogee_obs_list['s_ra'][0]
dec = apogee_obs_list['s_dec'][0]
# make a SkyCoord object from these coordinates
coord = SkyCoord(ra=ra*u.deg, dec=dec*u.deg)
print(coord)

# Search for TESS observations based on coordinates
tess_obs = Observations.query_criteria(coordinates=coord, # Search by coordinates
                                       radius=(1/60/60)*u.deg, # search within 1 arcsecond
                                       # Search for TESS observations
                                       obs_collection='TESS',
                                       # Select only Science observations (not calibration files)
                                       intentType='science',
                                       # Search for time series data (light curves)
                                       dataproduct_type='timeseries',
                                       dataURL='*_lc.fits'
                                       )

# Display Results
tess_obs

<SkyCoord (ICRS): (ra, dec) in deg
    (297.06439209, 43.12688065)>
Table masked=True length=4
intentTypeobs_collectionprovenance_nameinstrument_nameprojectfilterswavelength_regiontarget_nametarget_classificationobs_ids_ras_decdataproduct_typeproposal_picalib_levelt_mint_maxt_exptimeem_minem_maxobs_titlet_obs_releaseproposal_idproposal_typesequence_numbers_regionjpegURLdataURLdataRightsmtFlagsrcDenobsidobjIDobjID1distance
str7str4str4str10str4str4str7str9str1str47float64float64str10str14int64float64float64float64float64float64str1float64str7str1int64str47str1str73str6boolfloat64str8str9str9float64
scienceTESSSPOCPhotometerTESSTESSOptical272943843--tess2019226182529-s0015-0000000272943843-0151-s297.06444299685143.1269331469927timeseriesRicker, George358710.8662031458736.91161649120.0600.01000.0--58756.3333334G022062--15CIRCLE ICRS 297.06444300 43.12693315 0.00138889--mast:TESS/product/tess2019226182529-s0015-0000000272943843-0151-s_lc.fitsPUBLICFalsenan2748376470551762705517620.0
scienceTESSSPOCPhotometerTESSTESSOptical272943843--tess2021204101404-s0041-0000000272943843-0212-s297.06444299685143.1269331469927timeseriesRicker, George359419.490829687559446.08111446759120.0600.01000.0--59477.0G04028--41CIRCLE 297.064443 43.12693315 0.00138889--mast:TESS/product/tess2021204101404-s0041-0000000272943843-0212-s_lc.fitsPUBLICFalsenan628034641202463081202463080.0
scienceTESSSPOCPhotometerTESSTESSOptical272943843--tess2022190063128-s0054-0000000272943843-0227-s297.06444299685143.1269331469927timeseriesRicker, George359769.4001347569459795.63100395833120.0600.01000.0--59828.0G04028--54CIRCLE 297.064443 43.12693315 0.00138889--mast:TESS/product/tess2022190063128-s0054-0000000272943843-0227-s_lc.fitsPUBLICFalsenan915573751752899171752899170.0
scienceTESSSPOCPhotometerTESSTESSOptical272943843--tess2022217014003-s0055-0000000272943843-0242-s297.06444299685143.1269331469927timeseriesRicker, George359796.6018387268559823.76548444445120.0600.01000.0--59841.0G04028--55CIRCLE 297.064443 43.12693315 0.00138889--mast:TESS/product/tess2022217014003-s0055-0000000272943843-0242-s_lc.fitsPUBLICFalsenan937845111822515391822515390.0

There are a handful of results, but we can find the closest match to this star using the “distance” column which returns the distance (in arcsecs) of each match to our cone search. There are several results with distance of 0, (this star was observed in multiple TESS sectors), so we will just use the first result.

tess_obs[np.argmin(tess_obs['distance'])]
Row index=0 masked=True
intentTypeobs_collectionprovenance_nameinstrument_nameprojectfilterswavelength_regiontarget_nametarget_classificationobs_ids_ras_decdataproduct_typeproposal_picalib_levelt_mint_maxt_exptimeem_minem_maxobs_titlet_obs_releaseproposal_idproposal_typesequence_numbers_regionjpegURLdataURLdataRightsmtFlagsrcDenobsidobjIDobjID1distance
str7str4str4str10str4str4str7str9str1str47float64float64str10str14int64float64float64float64float64float64str1float64str7str1int64str47str1str73str6boolfloat64str8str9str9float64
scienceTESSSPOCPhotometerTESSTESSOptical272943843--tess2019226182529-s0015-0000000272943843-0151-s297.06444299685143.1269331469927timeseriesRicker, George358710.8662031458736.91161649120.0600.01000.0--58756.3333334G022062--15CIRCLE ICRS 297.06444300 43.12693315 0.00138889--mast:TESS/product/tess2019226182529-s0015-0000000272943843-0151-s_lc.fitsPUBLICFalsenan2748376470551762705517620.0

And download this file similar to what we did with APOGEE before:

# Get product list
tess_products = Observations.get_product_list(tess_obs[np.argmin(tess_obs['distance'])])

# Show products
tess_products
Table masked=True length=2
obsIDobs_collectiondataproduct_typeobs_iddescriptiontypedataURIproductTypeproductGroupDescriptionproductSubGroupDescriptionproductDocumentationURLprojectprvversionproposal_idproductFilenamesizeparent_obsiddataRightscalib_levelfilters
str8str4str10str47str18str1str73str7str28str2str1str4str20str7str55int64str8str6int64str4
27483764TESStimeseriestess2019226182529-s0015-0000000272943843-0151-sLight curvesSmast:TESS/product/tess2019226182529-s0015-0000000272943843-0151-s_lc.fitsSCIENCEMinimum Recommended ProductsLC--SPOCspoc-4.0.28-20200407G022062tess2019226182529-s0015-0000000272943843-0151-s_lc.fits190656027483764PUBLIC3TESS
27483764TESStimeseriestess2019226182529-s0015-0000000272943843-0151-sTarget pixel filesSmast:TESS/product/tess2019226182529-s0015-0000000272943843-0151-s_tp.fitsSCIENCEMinimum Recommended ProductsTP--SPOCspoc-4.0.28-20200407G022062tess2019226182529-s0015-0000000272943843-0151-s_tp.fits4595616027483764PUBLIC2TESS
# Download Products
manifest = Observations.download_products(tess_products, description='Light curves')

Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:TESS/product/tess2019226182529-s0015-0000000272943843-0151-s_lc.fits to ./mastDownload/TESS/tess2019226182529-s0015-0000000272943843-0151-s/tess2019226182529-s0015-0000000272943843-0151-s_lc.fits ...

 [Done]

Plotting a light curve from TESS#

The TESS light curve file has three extensions:

  • The primary header (HDU0)

  • The light curve data (HDU1)

  • aperture Information (HDU2)

# Open file
tess_lc = fits.open(manifest['Local Path'][0])

# Print file information
tess_lc.info()

Filename: ./mastDownload/TESS/tess2019226182529-s0015-0000000272943843-0151-s/tess2019226182529-s0015-0000000272943843-0151-s_lc.fits
No.    Name      Ver    Type      Cards   Dimensions   Format
  0  PRIMARY       1 PrimaryHDU      44   ()      
  1  LIGHTCURVE    1 BinTableHDU    161   18757R x 20C   [D, E, J, E, E, E, E, E, E, J, D, E, D, E, D, E, D, E, E, E]   
  2  APERTURE      1 ImageHDU        49   (11, 11)   int32

We can plot the light curve by using the TIME and SAP_FLUX data from the first extension:

plt.figure(figsize=(15, 5))

# Plot the observed spectrum
plt.plot(tess_lc['LIGHTCURVE'].data['TIME'], tess_lc['LIGHTCURVE'].data['SAP_FLUX'], c='k', label='TESS Light Curve')

# Plot a vertical line for every day
for i in range(int(np.nanmin(tess_lc['LIGHTCURVE'].data['TIME'])), int(np.nanmax(tess_lc['LIGHTCURVE'].data['TIME']))):
    plt.axvline(i, c='lightgrey')
    
# Set labels
plt.title(f"TESS-SPOC Light Curve - {tess_lc[0].header['OBJECT']}")
plt.xlabel('Time [BTJD]')
plt.ylabel('Flux')
plt.grid()

plt.show()
../../../_images/c793696a513df6a892cc2012ecaace5e037593d79a4238295b7e24bed2890064.png

This is a very typical light curve for a Cepheid Variable star, and TESS observed the variations in flux really well! This star cycles over a period of about 6 days.

Combining APOGEE and TESS data#

Now we are ready to combine everything: the HR Diagram, the spectrum, and the TESS light curve into a single plot!

Plotting APOGEE and TESS data together#

In the below cell, we define a helper function to plot all of the data together for our Cepheid Variable star:

def plot_spectrum_and_lc(apogee_spec, tess_lightcurve):
    """Helper function for plotting the HR Diagram, APOGEE Spectrum, and TESS Light Curve for a star"""
    # Set up axes for subplots
    plt.figure(figsize=(20, 10))
    ax1 = plt.subplot2grid((3, 5), (0, 0), colspan=2, rowspan=3)
    ax2 = plt.subplot2grid((3, 5), (0, 2), colspan=3, rowspan=1)
    ax3 = plt.subplot2grid((3, 5), (1, 2), colspan=3, rowspan=1)
    
    # ===================
    # AX1 - HR Diagram
    # ===================
    # Plot HR Diagram form APOGEE catalog
    # Only plot 1 out of every 10 points with [::10] to reduce plotting time
    im = ax1.scatter(allstar['TEFF'][::10], allstar['LOGG'][::10], c=allstar['FE_H'][::10],
                     marker='.', s=1, cmap='jet', vmin=-1, vmax=0.5, zorder=1)
            
    plt.colorbar(im, ax=ax1, location='bottom', label='Metallicity [Fe/H]')
            
    # Plot the specific star
    indx = np.where(allstar['APOGEE_ID'] == apogee_spec[4].data['APOGEE_ID'][0])[0][0]
    if allstar['TEFF'][indx]: # if the star has valid apogee parameters
        ax1.scatter(allstar['TEFF'][indx], allstar['LOGG'][indx], c=allstar['FE_H'][indx],
                    marker='*', s=1000, cmap='jet', edgecolor='k', lw=2,
                    vmin=-1, vmax=0.5, zorder=5, label='APOGEE Parameters')
    else: # params from aspcapstar (uncalibrated)
        print('WARNING: Quality Warning Flags on ASPCAP star parameters. ASPCAP parameters may not be reliable for this star.')
        ax1.scatter(apogee_spec[4].data['FPARAM'][0][0], apogee_spec[4].data['FPARAM'][0][1], c=apogee_spec[4].data['FPARAM'][0][3],
                    marker='*', s=1000, cmap='jet', edgecolor='k', lw=2, vmin=-1, vmax=0.5, zorder=5, label='APOGEE Parameters')
                
    # Labels and titles
    ax1.set_xlim(7000, 3000)
    ax1.set_ylim(5.5, -0.5)
    ax1.set_xlabel(r'Effective Temperature [$T_{\rm{eff}}$, K]')
    ax1.set_ylabel(r'Surface Gravity [$\log(g)$]')
    ax1.grid()
    ax1.set_title('APOGEE Stellar Parameters - allStar')
            
    # ===================
    # AX2 - APOGEE Spectrum
    # ===================        
    # Get the observed and model flux from extensions 1 and 3
    observed_flux = apogee_spec[1].data
    model_flux = apogee_spec[3].data
    
    # Mask the data using the error array in extenion 2
    pixel_mask = (apogee_spec[2].data < 0.1) # use only pixels with small errors (10% or less)
    
    # Plot the observed spectrum
    ax2.plot(apogee_wls[pixel_mask], observed_flux[pixel_mask], c='k', label='Observed Spectrum')
    
    # Plot the model spectrum
    ax2.plot(apogee_wls[pixel_mask], model_flux[pixel_mask], c='darkred', label='Best Fit Model Spectrum')
    
    # Set axes limits
    ax2.set_ylim(0.6, 1.1)
    ax2.set_xlim(np.min(apogee_wls), np.max(apogee_wls))
    
    apogee_id = apogee_spec[4].data['APOGEE_ID'][0] # String containing star's name
    ax2.set_title(f"APOGEE Spectrum - {apogee_id}")
    ax2.set_xlabel(r'Wavelength [$\AA$]')
    ax2.set_ylabel('Flux')
    ax2.grid()
    ax2.legend()
    
    # ===================
    # AX3 - TESS Light Curve
    # ===================
    # Plot Light Curve
    ax3.plot(tess_lightcurve['LIGHTCURVE'].data['TIME'], tess_lightcurve['LIGHTCURVE'].data['SAP_FLUX'], c='k', label='TESS Light Curve')
    
    # Plot a vertical line for every day
    for i in range(int(np.nanmin(tess_lightcurve['LIGHTCURVE'].data['TIME'])), int(np.nanmax(tess_lightcurve['LIGHTCURVE'].data['TIME']))):
        ax3.axvline(i, c='lightgrey')
        
    # Set labels
    ax3.set_title(f"TESS-SPOC Light Curve - {tess_lightcurve[0].header['OBJECT']}")
    ax3.set_xlabel('Time [BTJD]')
    ax3.set_ylabel('Flux')
    ax3.grid()
            
    plt.tight_layout()
    save_name = f"{apogee_id}_APOGEE_spec_TESS_lightcurve.png"
    plt.savefig(save_name, bbox_inches='tight')
    print(f"Saved to {save_name}")
    plt.show()

aspcap_star = fits.open('mastDownload/SDSS/sdss_apogee_apo1m_cepheid_v1154_cyg/aspcapStar-dr17-V1154_Cyg.fits')
tess_lc = fits.open(manifest['Local Path'][0])

plot_spectrum_and_lc(aspcap_star, tess_lc)

Saved to V1154_Cyg_APOGEE_spec_TESS_lightcurve.png
../../../_images/24b903482b2ab0097a7a2b425ecadb07a5f1f7cb3583315dc49c40988b30dbd5.png

Like we saw before, this star is a metal-poor horiztonal branch star, and has a period of about 6 days. A very cool example of a Cepheid Variable!

Exploring different types of Variable Stars#

So far in this notebook, we’ve plotted the stellar parameters, spectrum, and light curve of a Cepheid variable star. But what about other types of variables? Here are a few more stars in APOGEE that we will explore and plot in this section!

  • 2M11361176+8117369: a RR Lyrae Variable Star, a type of quickly-pulsating old, low-mass star

  • 2M19181706+5141323: a “Rotationally Variable” Star, a main sequence dwarf star with large starspots blocking out some of the flux on one side of the star as it rotates

  • 2M19223768+5044541: an Eclipsing Binary Star, a type of close binary system with two stars rotating around each other

Since we’re going to be going through the same process for each star (search for APOGEE data -> search for TESS data -> plot results), Let’s write a helper function:

def plot_variable_star(apogee_id: str) -> None:
    """
    This function takes a star name from APOGEE and downloads the APOGEE and TESS
    data for this star. It then plots the APOGEE spectrum and TESS Lightcurve together.

    Parameters:
    ============
    apogee_id: str
        name of star we want to plot
    """
    print("="*50)
    print(f"Plotting Spectrum and Light Curve for {apogee_id}")
    print("="*50)

    # Search for APOGEE Data
    print("Searching For APOGEE Data...")
    apogee_obs_list = Observations.query_criteria(provenance_name='APOGEE',
                                                  target_name=apogee_id)
    print("Downloading APOGEE Data...")
    # Download APOGEE Spectrum
    product_list = Observations.get_product_list(apogee_obs_list)
    manifest = Observations.download_products(product_list, productSubGroupDescription='ASPCAPSTAR')
    # Open APOGEE Spectrum
    aspcapstar = fits.open(manifest['Local Path'][0])

    # Coordinates of Star
    coord = SkyCoord(ra=apogee_obs_list['s_ra'][0]*u.deg, dec=apogee_obs_list['s_dec'][0]*u.deg)

    # Search for TESS observations based on coordinates
    print("Searching For TESS Data...")
    tess_obs = Observations.query_criteria(coordinates=coord, # Search by coordinates
                                           radius=(1/60/60)*u.deg, # search within 1 arcsecond
                                           obs_collection='TESS', # Search for TESS observations
                                           intentType='science', # Select only Science observations (not calibration files)
                                           dataproduct_type='timeseries', dataURL="*lc.fits") # Search for time series data (light curves)

    # Pick closest match by distance
    tess_obs = tess_obs[np.argmin(tess_obs['distance'])]

    # Download TESS Light Curve
    print("Downloading TESS Data...")
    tess_products = Observations.get_product_list(tess_obs)
    manifest = Observations.download_products(tess_products, description='Light curves')

    # Open light curve
    tess_lc = fits.open(manifest['Local Path'][0])

    # Plot star
    print("Making Plot...")
    plot_spectrum_and_lc(aspcapstar, tess_lc)

RR Lyrae Variable - 2M11361176+8117369#

# Define the star name we want to search for
star_name = '2M11361176+8117369'
# Call our helper function to download data and plot the results
plot_variable_star(star_name)

==================================================
Plotting Spectrum and Light Curve for 2M11361176+8117369
==================================================
Searching For APOGEE Data...
Downloading APOGEE Data...

Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:SDSS/apogee/apo1m/rrlyr/2M11361176%2B8117369/aspcapStar-dr17-2M11361176%2B8117369.fits to ./mastDownload/SDSS/sdss_apogee_apo1m_rrlyr_2m11361176+8117369/aspcapStar-dr17-2M11361176+8117369.fits ...

 [Done]

Searching For TESS Data...

Downloading TESS Data...

Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:TESS/product/tess2021175071901-s0040-0000000372758482-0211-s_lc.fits to ./mastDownload/TESS/tess2021175071901-s0040-0000000372758482-0211-s/tess2021175071901-s0040-0000000372758482-0211-s_lc.fits ...

 [Done]

Making Plot...

Saved to 2M11361176+8117369_APOGEE_spec_TESS_lightcurve.png
../../../_images/36e56ab90c936d73ff58027ca10fc7ab4e7120c8e8daac84575141c5af05beb1.png

This star is also a horizontal branch star, but it is hotter and more metal-poor than the Cepheid Variable was. It pulsates a lot more quickly too - from the light curve, the pulsation period is about 12 hours!

Rotationally Variable - 2M19181706+5141323#

star_name = '2M19181706+5141323'
plot_variable_star(star_name)

==================================================
Plotting Spectrum and Light Curve for 2M19181706+5141323
==================================================
Searching For APOGEE Data...
Downloading APOGEE Data...

Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:SDSS/apogee/apo25m/K01_082%2B17/2M19181706%2B5141323/aspcapStar-dr17-2M19181706%2B5141323.fits to ./mastDownload/SDSS/sdss_apogee_apo25m_k01_082+17_2m19181706+5141323/aspcapStar-dr17-2M19181706+5141323.fits ...

 [Done]

Searching For TESS Data...

Downloading TESS Data...

Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:TESS/product/tess2019226182529-s0015-0000000416267333-0151-s_lc.fits to ./mastDownload/TESS/tess2019226182529-s0015-0000000416267333-0151-s/tess2019226182529-s0015-0000000416267333-0151-s_lc.fits ...

 [Done]

Making Plot...

Saved to 2M19181706+5141323_APOGEE_spec_TESS_lightcurve.png
../../../_images/3cd4c8e95ee4226130bdc0bc056e28fe8b5a321dd4110ed6f970f81a4f2ccf49.png

This star is a main-sequence dwarf star, a lot like the Sun but a little bit cooler in temperature. Its light curve is not as periodic as the other two, because the variation in flux is caused by the stars rotation, not pulsation! Large star spots likely block out the flux on part of this star, causing it to dip in brightness roughly every 14 days.

Eclipsing Binary Star - 2M19203184+3830492#

star_name = '2M19203184+3830492'
plot_variable_star(star_name)

==================================================
Plotting Spectrum and Light Curve for 2M19203184+3830492
==================================================
Searching For APOGEE Data...
Downloading APOGEE Data...

Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:SDSS/apogee/apo25m/K21_071%2B10/2M19203184%2B3830492/aspcapStar-dr17-2M19203184%2B3830492.fits to ./mastDownload/SDSS/sdss_apogee_apo25m_k21_071+10_2m19203184+3830492/aspcapStar-dr17-2M19203184+3830492.fits ...

 [Done]

Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:SDSS/apogee/apo25m/N6791/2M19203184%2B3830492/aspcapStar-dr17-2M19203184%2B3830492.fits to ./mastDownload/SDSS/sdss_apogee_apo25m_n6791_2m19203184+3830492/aspcapStar-dr17-2M19203184+3830492.fits ...

 [Done]

Searching For TESS Data...

Downloading TESS Data...

Downloading URL https://mast.stsci.edu/api/v0.1/Download/file?uri=mast:TESS/product/tess2021175071901-s0040-0000000122304494-0211-s_lc.fits to ./mastDownload/TESS/tess2021175071901-s0040-0000000122304494-0211-s/tess2021175071901-s0040-0000000122304494-0211-s_lc.fits ...

 [Done]

Making Plot...
WARNING: Quality Warning Flags on ASPCAP star parameters. ASPCAP parameters may not be reliable for this star.

Saved to 2M19203184+3830492_APOGEE_spec_TESS_lightcurve.png
../../../_images/641be8f3eb42a8e90f098c7060b9b731ff3404ca0a95e8cd0f95197e775bd722.png

Every few days, this star’s flux dips dramatically in brightness by over 25%! This system is an eclipsing binary pair - a pair of stars orbiting around each other, and eclipsing one another along our line-of-sight, causing the fluctuations in brightness seen in the light curve.

Because this is an eclipsing binary pair, the ASPCAP parameters have a quality warning and may not be as reliable as the parameter estimates for other stars (which explains why the star is slightly off of the main sequence in the HR Diagram!). The single visit spectra (apVisit files) are probably more useful for most science cases involving binary stars than the combined spectrum (aspcapStar file) we plotted here. The apVisit files can be dowloaded using the productSubGroupDescription='APVISIT' keyword when downloading the products.

Congratulations! You have reached the end of this tutorial notebook. You have learned how to access and download APOGEE data from MAST, and combine it with TESS light curves to characterize different types of variable stars.

Additional Resources#

Additional resources are linked below:

Citations#

If you use MaNGA data for published research, see the following links for information on which citations to include in your paper:

About this Notebook#

Author(s): Julie Imig (jimig@stsci.edu)
Keyword(s): Tutorial, SDSS, APOGEE, TESS, stars
First published: January 2025
Last updated: January 2025

Top of Page Space Telescope Logo

Contents