Color-magnitude diagram: narrowband filters

In this tutorial, we will compute synthetic photometry for the narrowband H2 and H3 filters of VLT/SPHERE. We will use empirical data from the SpeX Prism Spectral Library and theoretical isochrones and spectra from the AMES-Cond and AMES-Dusty models. The synthetic colors and magnitudes are compared with the photometric data of directly imaged companions in a color-magnitude diagram.

Getting started

We start by importing the required Python modules.

[1]:
import urllib.request
import numpy as np
import species

Next, we initiate the species workflow with SpeciesInit. This will create the configuration file with default values and the HDF5 database.

[2]:
species.SpeciesInit()
Initiating species v0.4.0... [DONE]
Creating species_config.ini... [DONE]
Database: /Users/tomasstolker/applications/species/docs/tutorials/species_database.hdf5
Data folder: /Users/tomasstolker/applications/species/docs/tutorials/data
Working folder: /Users/tomasstolker/applications/species/docs/tutorials
Creating species_database.hdf5... [DONE]
Creating data folder... [DONE]
[2]:
<species.core.init.SpeciesInit at 0x14aa55a90>

Later on, we will interpolate the isochrone data at an age of 20 Myr for 100 logarithmically-spaced masses between 1 MJup and 1000 MJup.

[3]:
age = 20.  # (Myr)
masses = np.logspace(0., 3., 100)  # (Mjup)

Adding data and models

We will now add the required data to the database by first creating an instance of Database.

[4]:
database = species.Database()

The photometric data of directly imaged companions that are available in species is added by using the add_companion method with name=None. Alternative, add_object can be used to manually add photometric and spectroscopic data of an individual object.

[5]:
database.add_companion(name=None)
Getting GPI_YJHK spectrum of beta Pic b... [DONE]
IMPORTANT: Please cite Chilcote et al. 2017, AJ, 153, 182
           when making use of this spectrum in a publication
Getting GRAVITY spectrum of beta Pic b... [DONE]
IMPORTANT: Please cite Gravity Collaboration et al. 2020, A&A, 633, 110
           when making use of this spectrum in a publication
Downloading Vega spectrum (270 kB)... [DONE]
Adding Vega spectrum... [DONE]
Adding filter: Magellan/VisAO.Ys... [DONE]
Adding filter: Paranal/NACO.J... [DONE]
Adding filter: Gemini/NICI.ED286... [DONE]
Adding filter: Paranal/NACO.H... [DONE]
Adding filter: Paranal/NACO.Ks... [DONE]
Adding filter: Paranal/NACO.NB374... [DONE]
Adding filter: Paranal/NACO.Lp... [DONE]
Adding filter: Paranal/NACO.NB405... [DONE]
Adding filter: Paranal/NACO.Mp... [DONE]
Adding object: beta Pic b
   - Distance (pc) = 19.75 +/- 0.13
   - Magellan/VisAO.Ys:
      - Apparent magnitude = 15.53 +/- 0.34
      - Flux (W m-2 um-1) = 4.25e-15 +/- 1.35e-15
   - Paranal/NACO.J:
      - Apparent magnitude = 14.11 +/- 0.21
      - Flux (W m-2 um-1) = 6.87e-15 +/- 1.34e-15
   - Gemini/NICI.ED286:
      - Apparent magnitude = 13.18 +/- 0.15
      - Flux (W m-2 um-1) = 6.99e-15 +/- 9.69e-16
   - Paranal/NACO.H:
      - Apparent magnitude = 13.32 +/- 0.14
      - Flux (W m-2 um-1) = 5.47e-15 +/- 7.08e-16
   - Paranal/NACO.Ks:
      - Apparent magnitude = 12.64 +/- 0.11
      - Flux (W m-2 um-1) = 4.04e-15 +/- 4.10e-16
   - Paranal/NACO.NB374:
      - Apparent magnitude = 11.25 +/- 0.23
      - Flux (W m-2 um-1) = 1.69e-15 +/- 3.61e-16
   - Paranal/NACO.Lp:
      - Apparent magnitude = 11.30 +/- 0.06
      - Flux (W m-2 um-1) = 1.59e-15 +/- 8.79e-17
   - Paranal/NACO.NB405:
      - Apparent magnitude = 10.98 +/- 0.05
      - Flux (W m-2 um-1) = 1.61e-15 +/- 7.42e-17
   - Paranal/NACO.Mp:
      - Apparent magnitude = 11.10 +/- 0.12
      - Flux (W m-2 um-1) = 7.86e-16 +/- 8.70e-17
   - Spectrum:
      - Database tag: GPI_YJHK
      - Filename: ./data/companion_data/betapicb_gpi_yjhk.dat
      - Data shape: (185, 3)
      - Wavelength range (um): 0.95 - 2.40
      - Mean flux (W m-2 um-1): 5.40e-15
      - Mean error (W m-2 um-1): 6.44e-16
   - GRAVITY spectrum:
      - Object: Unknown
      - Database tag: GRAVITY
      - Filename: ./data/companion_data/BetaPictorisb_2018-09-22.fits
      - Data shape: (237, 3)
      - Wavelength range (um): 1.97 - 2.49
      - Mean flux (W m-2 um-1): 4.65e-15
      - Mean error (W m-2 um-1): 1.00e-16
   - GRAVITY covariance matrix:
      - Object: Unknown
      - Database tag: GRAVITY
      - Filename: ./data/companion_data/BetaPictorisb_2018-09-22.fits
      - Data shape: (237, 237)
   - Spectral resolution:
      - GPI_YJHK: 40.0
      - GRAVITY: 500.0
Adding filter: MKO/NSFCam.K... [DONE]
Adding object: beta Pic c
   - Distance (pc) = 19.75 +/- 0.13
   - MKO/NSFCam.K:
      - Apparent magnitude = 14.30 +/- 0.10
      - Flux (W m-2 um-1) = 7.93e-16 +/- 7.31e-17
Adding filter: Paranal/SPHERE.IRDIS_D_H23_2... [DONE]
Adding filter: Paranal/SPHERE.IRDIS_D_H23_3... [DONE]
Adding filter: Paranal/SPHERE.IRDIS_D_K12_1... [DONE]
Adding filter: Paranal/SPHERE.IRDIS_D_K12_2... [DONE]
Adding object: HIP 65426 b
   - Distance (pc) = 109.21 +/- 0.75
   - Paranal/SPHERE.IRDIS_D_H23_2:
      - Apparent magnitude = 17.94 +/- 0.05
      - Flux (W m-2 um-1) = 8.73e-17 +/- 4.02e-18
   - Paranal/SPHERE.IRDIS_D_H23_3:
      - Apparent magnitude = 17.58 +/- 0.06
      - Flux (W m-2 um-1) = 1.03e-16 +/- 5.70e-18
   - Paranal/SPHERE.IRDIS_D_K12_1:
      - Apparent magnitude = 17.01 +/- 0.09
      - Flux (W m-2 um-1) = 7.63e-17 +/- 6.33e-18
   - Paranal/SPHERE.IRDIS_D_K12_2:
      - Apparent magnitude = 16.79 +/- 0.09
      - Flux (W m-2 um-1) = 7.21e-17 +/- 5.98e-18
   - Paranal/NACO.Lp:
      - Apparent magnitude = 15.33 +/- 0.12
      - Flux (W m-2 um-1) = 3.88e-17 +/- 4.30e-18
   - Paranal/NACO.NB405:
      - Apparent magnitude = 15.23 +/- 0.22
      - Flux (W m-2 um-1) = 3.21e-17 +/- 6.55e-18
   - Paranal/NACO.Mp:
      - Apparent magnitude = 14.65 +/- 0.29
      - Flux (W m-2 um-1) = 2.99e-17 +/- 8.07e-18
Getting SPHERE_YJH spectrum of 51 Eri b... [DONE]
IMPORTANT: Please cite Samland et al. 2017, A&A, 603, 57
           when making use of this spectrum in a publication
Adding filter: MKO/NSFCam.J... [DONE]
Adding filter: MKO/NSFCam.H... [DONE]
Adding filter: Paranal/SPHERE.IRDIS_B_H... [DONE]
Adding filter: Keck/NIRC2.Lp... [DONE]
Adding filter: Keck/NIRC2.Ms... [DONE]
Adding object: 51 Eri b
   - Distance (pc) = 29.78 +/- 0.12
   - MKO/NSFCam.J:
      - Apparent magnitude = 19.04 +/- 0.40
      - Flux (W m-2 um-1) = 7.52e-17 +/- 2.83e-17
   - MKO/NSFCam.H:
      - Apparent magnitude = 18.99 +/- 0.21
      - Flux (W m-2 um-1) = 3.12e-17 +/- 6.07e-18
   - MKO/NSFCam.K:
      - Apparent magnitude = 18.67 +/- 0.19
      - Flux (W m-2 um-1) = 1.42e-17 +/- 2.49e-18
   - Paranal/SPHERE.IRDIS_B_H:
      - Apparent magnitude = 19.45 +/- 0.29
      - Flux (W m-2 um-1) = 2.06e-17 +/- 5.57e-18
   - Paranal/SPHERE.IRDIS_D_H23_2:
      - Apparent magnitude = 18.41 +/- 0.26
      - Flux (W m-2 um-1) = 5.66e-17 +/- 1.37e-17
   - Paranal/SPHERE.IRDIS_D_K12_1:
      - Apparent magnitude = 17.55 +/- 0.14
      - Flux (W m-2 um-1) = 4.64e-17 +/- 6.00e-18
   - Keck/NIRC2.Lp:
      - Apparent magnitude = 16.20 +/- 0.11
      - Flux (W m-2 um-1) = 1.79e-17 +/- 1.81e-18
   - Keck/NIRC2.Ms:
      - Apparent magnitude = 16.10 +/- 0.50
      - Flux (W m-2 um-1) = 8.47e-18 +/- 4.04e-18
   - Spectrum:
      - Database tag: SPHERE_YJH
      - Filename: ./data/companion_data/51erib_sphere_yjh.dat
      - Data shape: (33, 3)
      - Wavelength range (um): 1.01 - 1.64
      - Mean flux (W m-2 um-1): 3.50e-17
      - Mean error (W m-2 um-1): 1.06e-17
   - Spectral resolution:
      - SPHERE_YJH: 25.0
Adding filter: Subaru/CIAO.z...
/Users/tomasstolker/applications/species/species/data/filters.py:188: UserWarning: The minimum transmission value of Subaru/CIAO.z is smaller than zero (-1.80e-03). Wavelengths with negative transmission values will be removed.
  warnings.warn(f'The minimum transmission value of {filter_id} is smaller than zero '
 [DONE]
Adding filter: Paranal/SPHERE.IRDIS_B_J... [DONE]
Adding filter: Keck/NIRC2.H... [DONE]
Adding filter: Keck/NIRC2.Ks... [DONE]
Adding object: HR 8799 b
   - Distance (pc) = 41.29 +/- 0.15
   - Subaru/CIAO.z:
      - Apparent magnitude = 21.22 +/- 0.29
      - Flux (W m-2 um-1) = 1.93e-17 +/- 5.22e-18
   - Paranal/SPHERE.IRDIS_B_J:
      - Apparent magnitude = 19.78 +/- 0.09
      - Flux (W m-2 um-1) = 3.90e-17 +/- 3.24e-18
   - Keck/NIRC2.H:
      - Apparent magnitude = 18.05 +/- 0.09
      - Flux (W m-2 um-1) = 7.32e-17 +/- 6.08e-18
   - Paranal/SPHERE.IRDIS_D_H23_2:
      - Apparent magnitude = 18.08 +/- 0.14
      - Flux (W m-2 um-1) = 7.67e-17 +/- 9.92e-18
   - Paranal/SPHERE.IRDIS_D_H23_3:
      - Apparent magnitude = 17.78 +/- 0.10
      - Flux (W m-2 um-1) = 8.57e-17 +/- 7.91e-18
   - Keck/NIRC2.Ks:
      - Apparent magnitude = 17.03 +/- 0.08
      - Flux (W m-2 um-1) = 6.98e-17 +/- 5.15e-18
   - Paranal/SPHERE.IRDIS_D_K12_1:
      - Apparent magnitude = 17.15 +/- 0.06
      - Flux (W m-2 um-1) = 6.71e-17 +/- 3.71e-18
   - Paranal/SPHERE.IRDIS_D_K12_2:
      - Apparent magnitude = 16.97 +/- 0.09
      - Flux (W m-2 um-1) = 6.11e-17 +/- 5.07e-18
   - Paranal/NACO.Lp:
      - Apparent magnitude = 15.52 +/- 0.10
      - Flux (W m-2 um-1) = 3.26e-17 +/- 3.01e-18
   - Paranal/NACO.NB405:
      - Apparent magnitude = 14.82 +/- 0.18
      - Flux (W m-2 um-1) = 4.69e-17 +/- 7.80e-18
   - Keck/NIRC2.Ms:
      - Apparent magnitude = 16.05 +/- 0.30
      - Flux (W m-2 um-1) = 8.87e-18 +/- 2.48e-18
Adding object: HR 8799 c
   - Distance (pc) = 41.29 +/- 0.15
   - Paranal/SPHERE.IRDIS_B_J:
      - Apparent magnitude = 18.60 +/- 0.13
      - Flux (W m-2 um-1) = 1.16e-16 +/- 1.39e-17
   - Keck/NIRC2.H:
      - Apparent magnitude = 17.06 +/- 0.13
      - Flux (W m-2 um-1) = 1.82e-16 +/- 2.19e-17
   - Paranal/SPHERE.IRDIS_D_H23_2:
      - Apparent magnitude = 17.09 +/- 0.12
      - Flux (W m-2 um-1) = 1.91e-16 +/- 2.12e-17
   - Paranal/SPHERE.IRDIS_D_H23_3:
      - Apparent magnitude = 16.78 +/- 0.10
      - Flux (W m-2 um-1) = 2.15e-16 +/- 1.99e-17
   - Keck/NIRC2.Ks:
      - Apparent magnitude = 16.11 +/- 0.08
      - Flux (W m-2 um-1) = 1.63e-16 +/- 1.20e-17
   - Paranal/SPHERE.IRDIS_D_K12_1:
      - Apparent magnitude = 16.19 +/- 0.05
      - Flux (W m-2 um-1) = 1.62e-16 +/- 7.48e-18
   - Paranal/SPHERE.IRDIS_D_K12_2:
      - Apparent magnitude = 15.86 +/- 0.07
      - Flux (W m-2 um-1) = 1.70e-16 +/- 1.10e-17
   - Paranal/NACO.Lp:
      - Apparent magnitude = 14.65 +/- 0.11
      - Flux (W m-2 um-1) = 7.27e-17 +/- 7.38e-18
   - Paranal/NACO.NB405:
      - Apparent magnitude = 13.97 +/- 0.11
      - Flux (W m-2 um-1) = 1.03e-16 +/- 1.04e-17
   - Keck/NIRC2.Ms:
      - Apparent magnitude = 15.03 +/- 0.14
      - Flux (W m-2 um-1) = 2.27e-17 +/- 2.93e-18
Adding object: HR 8799 d
   - Distance (pc) = 41.29 +/- 0.15
   - Paranal/SPHERE.IRDIS_B_J:
      - Apparent magnitude = 18.59 +/- 0.37
      - Flux (W m-2 um-1) = 1.17e-16 +/- 4.06e-17
   - Keck/NIRC2.H:
      - Apparent magnitude = 16.71 +/- 0.24
      - Flux (W m-2 um-1) = 2.52e-16 +/- 5.61e-17
   - Paranal/SPHERE.IRDIS_D_H23_2:
      - Apparent magnitude = 17.02 +/- 0.17
      - Flux (W m-2 um-1) = 2.04e-16 +/- 3.20e-17
   - Paranal/SPHERE.IRDIS_D_H23_3:
      - Apparent magnitude = 16.85 +/- 0.16
      - Flux (W m-2 um-1) = 2.02e-16 +/- 2.99e-17
   - Keck/NIRC2.Ks:
      - Apparent magnitude = 16.09 +/- 0.12
      - Flux (W m-2 um-1) = 1.66e-16 +/- 1.84e-17
   - Paranal/SPHERE.IRDIS_D_K12_1:
      - Apparent magnitude = 16.20 +/- 0.07
      - Flux (W m-2 um-1) = 1.61e-16 +/- 1.04e-17
   - Paranal/SPHERE.IRDIS_D_K12_2:
      - Apparent magnitude = 15.84 +/- 0.10
      - Flux (W m-2 um-1) = 1.73e-16 +/- 1.60e-17
   - Paranal/NACO.Lp:
      - Apparent magnitude = 14.55 +/- 0.14
      - Flux (W m-2 um-1) = 7.97e-17 +/- 1.03e-17
   - Paranal/NACO.NB405:
      - Apparent magnitude = 13.87 +/- 0.15
      - Flux (W m-2 um-1) = 1.12e-16 +/- 1.56e-17
   - Keck/NIRC2.Ms:
      - Apparent magnitude = 14.65 +/- 0.35
      - Flux (W m-2 um-1) = 3.22e-17 +/- 1.06e-17
Getting SPHERE_YJH spectrum of HR 8799 e... [DONE]
IMPORTANT: Please cite Zurlo et al. 2016, A&A, 587, 57
           when making use of this spectrum in a publication
Adding object: HR 8799 e
   - Distance (pc) = 41.29 +/- 0.15
   - Paranal/SPHERE.IRDIS_B_J:
      - Apparent magnitude = 18.40 +/- 0.21
      - Flux (W m-2 um-1) = 1.39e-16 +/- 2.71e-17
   - Paranal/SPHERE.IRDIS_D_H23_2:
      - Apparent magnitude = 16.91 +/- 0.20
      - Flux (W m-2 um-1) = 2.25e-16 +/- 4.18e-17
   - Paranal/SPHERE.IRDIS_D_H23_3:
      - Apparent magnitude = 16.68 +/- 0.21
      - Flux (W m-2 um-1) = 2.36e-16 +/- 4.59e-17
   - Keck/NIRC2.Ks:
      - Apparent magnitude = 15.91 +/- 0.22
      - Flux (W m-2 um-1) = 1.96e-16 +/- 4.00e-17
   - Paranal/SPHERE.IRDIS_D_K12_1:
      - Apparent magnitude = 16.12 +/- 0.10
      - Flux (W m-2 um-1) = 1.73e-16 +/- 1.60e-17
   - Paranal/SPHERE.IRDIS_D_K12_2:
      - Apparent magnitude = 15.82 +/- 0.11
      - Flux (W m-2 um-1) = 1.76e-16 +/- 1.79e-17
   - Paranal/NACO.Lp:
      - Apparent magnitude = 14.49 +/- 0.21
      - Flux (W m-2 um-1) = 8.42e-17 +/- 1.64e-17
   - Paranal/NACO.NB405:
      - Apparent magnitude = 13.72 +/- 0.20
      - Flux (W m-2 um-1) = 1.29e-16 +/- 2.39e-17
   - Spectrum:
      - Database tag: SPHERE_YJH
      - Filename: ./data/companion_data/hr8799e_sphere_yjh.dat
      - Data shape: (39, 3)
      - Wavelength range (um): 0.94 - 1.65
      - Mean flux (W m-2 um-1): 9.61e-17
      - Mean error (W m-2 um-1): 3.20e-17
   - Spectral resolution:
      - SPHERE_YJH: 25.0
Adding filter: Gemini/GPI.H... [DONE]
Adding filter: Gemini/GPI.K1... [DONE]
Adding object: HD 95086 b
   - Distance (pc) = 86.44 +/- 0.24
   - Gemini/GPI.H:
      - Apparent magnitude = 20.51 +/- 0.25
      - Flux (W m-2 um-1) = 7.41e-18 +/- 1.72e-18
   - Gemini/GPI.K1:
      - Apparent magnitude = 18.99 +/- 0.20
      - Flux (W m-2 um-1) = 1.38e-17 +/- 2.56e-18
   - Paranal/NACO.Lp:
      - Apparent magnitude = 16.27 +/- 0.19
      - Flux (W m-2 um-1) = 1.63e-17 +/- 2.87e-18
Getting SPHERE_YJH spectrum of PDS 70 b... [DONE]
IMPORTANT: Please cite Müller et al. 2018, A&A, 617, 2
           when making use of this spectrum in a publication
Adding object: PDS 70 b
   - Distance (pc) = 113.43 +/- 0.52
   - Paranal/SPHERE.IRDIS_D_H23_2:
      - Apparent magnitude = 18.12 +/- 0.21
      - Flux (W m-2 um-1) = 7.40e-17 +/- 1.44e-17
   - Paranal/SPHERE.IRDIS_D_H23_3:
      - Apparent magnitude = 17.97 +/- 0.18
      - Flux (W m-2 um-1) = 7.20e-17 +/- 1.20e-17
   - Paranal/SPHERE.IRDIS_D_K12_1:
      - Apparent magnitude = 16.66 +/- 0.04
      - Flux (W m-2 um-1) = 1.05e-16 +/- 3.88e-18
   - Paranal/SPHERE.IRDIS_D_K12_2:
      - Apparent magnitude = 16.37 +/- 0.06
      - Flux (W m-2 um-1) = 1.06e-16 +/- 5.87e-18
   - MKO/NSFCam.J:
      - Apparent magnitude = 20.04 +/- 0.09
      - Flux (W m-2 um-1) = 2.99e-17 +/- 2.48e-18
   - MKO/NSFCam.H:
      - Apparent magnitude = 18.24 +/- 0.04
      - Flux (W m-2 um-1) = 6.22e-17 +/- 2.29e-18
   - Paranal/NACO.Lp:
      - Apparent magnitude = 14.68 +/- 0.22
      - Flux (W m-2 um-1) = 7.07e-17 +/- 1.44e-17
   - Paranal/NACO.NB405:
      - Apparent magnitude = 14.68 +/- 0.27
      - Flux (W m-2 um-1) = 5.33e-17 +/- 1.34e-17
   - Paranal/NACO.Mp:
      - Apparent magnitude = 13.80 +/- 0.27
      - Flux (W m-2 um-1) = 6.53e-17 +/- 1.64e-17
   - Keck/NIRC2.Lp:
      - Apparent magnitude = 14.64 +/- 0.18
      - Flux (W m-2 um-1) = 7.52e-17 +/- 1.25e-17
   - Spectrum:
      - Database tag: SPHERE_YJH
      - Filename: ./data/companion_data/pds70b_sphere_yjh.dat
      - Data shape: (39, 3)
      - Wavelength range (um): 0.96 - 1.64
      - Mean flux (W m-2 um-1): 3.37e-17
      - Mean error (W m-2 um-1): 1.09e-17
   - Spectral resolution:
      - SPHERE_YJH: 25.0
Adding object: PDS 70 c
   - Distance (pc) = 113.43 +/- 0.52
   - Paranal/NACO.NB405:
      - Apparent magnitude = 14.91 +/- 0.35
      - Flux (W m-2 um-1) = 4.31e-17 +/- 1.41e-17
   - Keck/NIRC2.Lp:
      - Apparent magnitude = 15.50 +/- 0.46
      - Flux (W m-2 um-1) = 3.40e-17 +/- 1.49e-17
Adding filter: HST/NICMOS1.F090M... [DONE]
Adding filter: HST/NICMOS1.F110M... [DONE]
Adding filter: HST/NICMOS1.F145M... [DONE]
Adding filter: HST/NICMOS1.F160W... [DONE]
Adding object: 2M1207 B
   - Distance (pc) = 64.42 +/- 0.65
   - HST/NICMOS1.F090M:
      - Apparent magnitude = 22.58 +/- 0.35
      - Flux (W m-2 um-1) = 7.90e-18 +/- 2.59e-18
   - HST/NICMOS1.F110M:
      - Apparent magnitude = 20.61 +/- 0.15
      - Flux (W m-2 um-1) = 2.71e-17 +/- 3.76e-18
   - HST/NICMOS1.F145M:
      - Apparent magnitude = 19.05 +/- 0.03
      - Flux (W m-2 um-1) = 4.33e-17 +/- 1.20e-18
   - HST/NICMOS1.F160W:
      - Apparent magnitude = 18.27 +/- 0.02
      - Flux (W m-2 um-1) = 6.37e-17 +/- 1.17e-18
   - Paranal/NACO.J:
      - Apparent magnitude = 20.00 +/- 0.20
      - Flux (W m-2 um-1) = 3.03e-17 +/- 5.61e-18
   - Paranal/NACO.H:
      - Apparent magnitude = 18.09 +/- 0.21
      - Flux (W m-2 um-1) = 6.76e-17 +/- 1.32e-17
   - Paranal/NACO.Ks:
      - Apparent magnitude = 16.93 +/- 0.11
      - Flux (W m-2 um-1) = 7.78e-17 +/- 7.89e-18
   - Paranal/NACO.Lp:
      - Apparent magnitude = 15.28 +/- 0.14
      - Flux (W m-2 um-1) = 4.07e-17 +/- 5.26e-18
Adding object: AB Pic B
   - Distance (pc) = 50.12 +/- 0.07
   - Paranal/NACO.J:
      - Apparent magnitude = 16.18 +/- 0.10
      - Flux (W m-2 um-1) = 1.02e-15 +/- 9.41e-17
   - Paranal/NACO.H:
      - Apparent magnitude = 14.69 +/- 0.10
      - Flux (W m-2 um-1) = 1.55e-15 +/- 1.43e-16
   - Paranal/NACO.Ks:
      - Apparent magnitude = 14.14 +/- 0.08
      - Flux (W m-2 um-1) = 1.02e-15 +/- 7.49e-17
Getting SPHERE_YJH spectrum of HD 206893 B... [DONE]
IMPORTANT: Please cite Delorme et al. 2017, A&A, 608, 79
           when making use of this spectrum in a publication
Adding object: HD 206893 B
   - Distance (pc) = 40.81 +/- 0.11
   - Paranal/SPHERE.IRDIS_B_H:
      - Apparent magnitude = 16.79 +/- 0.06
      - Flux (W m-2 um-1) = 2.39e-16 +/- 1.32e-17
   - Paranal/SPHERE.IRDIS_D_K12_1:
      - Apparent magnitude = 15.20 +/- 0.10
      - Flux (W m-2 um-1) = 4.04e-16 +/- 3.73e-17
   - Paranal/SPHERE.IRDIS_D_K12_2:
      - Apparent magnitude = 14.88 +/- 0.09
      - Flux (W m-2 um-1) = 4.19e-16 +/- 3.48e-17
   - Paranal/NACO.Lp:
      - Apparent magnitude = 13.79 +/- 0.31
      - Flux (W m-2 um-1) = 1.60e-16 +/- 4.64e-17
   - Paranal/NACO.NB405:
      - Apparent magnitude = 13.16 +/- 0.34
      - Flux (W m-2 um-1) = 2.16e-16 +/- 6.88e-17
   - Paranal/NACO.Mp:
      - Apparent magnitude = 12.77 +/- 0.27
      - Flux (W m-2 um-1) = 1.69e-16 +/- 4.24e-17
   - Spectrum:
      - Database tag: SPHERE_YJH
      - Filename: ./data/companion_data/hd206893b_sphere_yjh.dat
      - Data shape: (39, 3)
      - Wavelength range (um): 0.96 - 1.64
      - Mean flux (W m-2 um-1): 1.31e-16
      - Mean error (W m-2 um-1): 2.81e-17
   - Spectral resolution:
      - SPHERE_YJH: 25.0
Adding filter: Paranal/SPHERE.IRDIS_B_Ks... [DONE]
Adding object: RZ Psc B
   - Distance (pc) = 195.86 +/- 4.03
   - Paranal/SPHERE.IRDIS_B_H (2 values):
      - Apparent magnitude = 13.71 +/- 0.14
      - Flux (W m-2 um-1) = 4.07e-15 +/- 5.26e-16
      - Apparent magnitude = 13.85 +/- 0.26
      - Flux (W m-2 um-1) = 3.58e-15 +/- 8.65e-16
   - Paranal/SPHERE.IRDIS_B_Ks:
      - Apparent magnitude = 13.51 +/- 0.20
      - Flux (W m-2 um-1) = 1.69e-15 +/- 3.12e-16
Adding filter: HST/WFPC2-PC.F606W... [DONE]
Adding filter: HST/WFPC2-PC.F814W... [DONE]
Adding filter: HST/NICMOS2.F171M... [DONE]
Adding filter: HST/NICMOS2.F190N... [DONE]
Adding filter: HST/NICMOS2.F215N... [DONE]
Adding filter: Magellan/VisAO.ip... [DONE]
Adding filter: Magellan/VisAO.zp... [DONE]
Adding filter: Subaru/CIAO.CH4s... [DONE]
Adding filter: Subaru/CIAO.K... [DONE]
Adding filter: Subaru/CIAO.Lp... [DONE]
Adding object: GQ Lup B
   - Distance (pc) = 154.10 +/- 0.69
   - HST/WFPC2-PC.F606W:
      - Apparent magnitude = 19.19 +/- 0.07
      - Flux (W m-2 um-1) = 5.94e-16 +/- 3.83e-17
   - HST/WFPC2-PC.F814W:
      - Apparent magnitude = 17.67 +/- 0.05
      - Flux (W m-2 um-1) = 1.02e-15 +/- 4.68e-17
   - HST/NICMOS2.F171M:
      - Apparent magnitude = 13.84 +/- 0.13
      - Flux (W m-2 um-1) = 2.91e-15 +/- 3.49e-16
   - HST/NICMOS2.F190N:
      - Apparent magnitude = 14.08 +/- 0.20
      - Flux (W m-2 um-1) = 1.65e-15 +/- 3.06e-16
   - HST/NICMOS2.F215N:
      - Apparent magnitude = 13.40 +/- 0.15
      - Flux (W m-2 um-1) = 1.94e-15 +/- 2.69e-16
   - Magellan/VisAO.ip:
      - Apparent magnitude = 18.89 +/- 0.24
      - Flux (W m-2 um-1) = 3.74e-16 +/- 8.33e-17
   - Magellan/VisAO.zp:
      - Apparent magnitude = 16.40 +/- 0.10
      - Flux (W m-2 um-1) = 2.31e-15 +/- 2.13e-16
   - Magellan/VisAO.Ys:
      - Apparent magnitude = 15.88 +/- 0.10
      - Flux (W m-2 um-1) = 3.08e-15 +/- 2.84e-16
   - MKO/NSFCam.H:
      - Apparent magnitude = 14.01 +/- 0.13
      - Flux (W m-2 um-1) = 3.06e-15 +/- 3.67e-16
   - Paranal/NACO.NB405:
      - Apparent magnitude = 12.29 +/- 0.06
      - Flux (W m-2 um-1) = 4.82e-16 +/- 2.66e-17
   - Paranal/NACO.Mp:
      - Apparent magnitude = 11.97 +/- 0.08
      - Flux (W m-2 um-1) = 3.53e-16 +/- 2.60e-17
   - Paranal/NACO.Ks (4 values):
      - Apparent magnitude = 13.47 +/- 0.03
      - Flux (W m-2 um-1) = 1.88e-15 +/- 5.36e-17
      - Apparent magnitude = 13.39 +/- 0.03
      - Flux (W m-2 um-1) = 2.03e-15 +/- 6.00e-17
      - Apparent magnitude = 13.50 +/- 0.05
      - Flux (W m-2 um-1) = 1.84e-15 +/- 8.47e-17
      - Apparent magnitude = 13.50 +/- 0.03
      - Flux (W m-2 um-1) = 1.83e-15 +/- 4.72e-17
   - Subaru/CIAO.CH4s:
      - Apparent magnitude = 13.76 +/- 0.26
      - Flux (W m-2 um-1) = 3.96e-15 +/- 9.58e-16
   - Subaru/CIAO.K:
      - Apparent magnitude = 13.37 +/- 0.12
      - Flux (W m-2 um-1) = 1.87e-15 +/- 2.07e-16
   - Subaru/CIAO.Lp:
      - Apparent magnitude = 12.44 +/- 0.22
      - Flux (W m-2 um-1) = 5.78e-16 +/- 1.18e-16
Adding filter: Paranal/SPHERE.ZIMPOL_R_PRIM... [DONE]
Adding filter: Paranal/SPHERE.ZIMPOL_I_PRIM... [DONE]
Adding filter: Gemini/NIRI.H2S1v2-1-G0220... [DONE]
Adding object: PZ Tel B
   - Distance (pc) = 47.13 +/- 0.13
   - Paranal/SPHERE.ZIMPOL_R_PRIM:
      - Apparent magnitude = 17.84 +/- 0.31
      - Flux (W m-2 um-1) = 1.83e-15 +/- 5.29e-16
   - Paranal/SPHERE.ZIMPOL_I_PRIM:
      - Apparent magnitude = 15.16 +/- 0.12
      - Flux (W m-2 um-1) = 1.09e-14 +/- 1.20e-15
   - Paranal/SPHERE.IRDIS_D_H23_2:
      - Apparent magnitude = 11.78 +/- 0.19
      - Flux (W m-2 um-1) = 2.54e-14 +/- 4.47e-15
   - Paranal/SPHERE.IRDIS_D_H23_3:
      - Apparent magnitude = 11.65 +/- 0.19
      - Flux (W m-2 um-1) = 2.43e-14 +/- 4.27e-15
   - Paranal/SPHERE.IRDIS_D_K12_1:
      - Apparent magnitude = 11.56 +/- 0.09
      - Flux (W m-2 um-1) = 1.15e-14 +/- 9.58e-16
   - Paranal/SPHERE.IRDIS_D_K12_2:
      - Apparent magnitude = 11.29 +/- 0.10
      - Flux (W m-2 um-1) = 1.14e-14 +/- 1.05e-15
   - Paranal/NACO.J:
      - Apparent magnitude = 12.47 +/- 0.20
      - Flux (W m-2 um-1) = 3.11e-14 +/- 5.76e-15
   - Paranal/NACO.H:
      - Apparent magnitude = 11.93 +/- 0.14
      - Flux (W m-2 um-1) = 1.97e-14 +/- 2.55e-15
   - Paranal/NACO.Ks:
      - Apparent magnitude = 11.53 +/- 0.07
      - Flux (W m-2 um-1) = 1.12e-14 +/- 7.25e-16
   - Paranal/NACO.Lp:
      - Apparent magnitude = 11.04 +/- 0.22
      - Flux (W m-2 um-1) = 2.02e-15 +/- 4.12e-16
   - Paranal/NACO.NB405:
      - Apparent magnitude = 10.94 +/- 0.07
      - Flux (W m-2 um-1) = 1.67e-15 +/- 1.08e-16
   - Paranal/NACO.Mp:
      - Apparent magnitude = 10.93 +/- 0.03
      - Flux (W m-2 um-1) = 9.19e-16 +/- 2.54e-17
   - Gemini/NICI.ED286:
      - Apparent magnitude = 11.68 +/- 0.14
      - Flux (W m-2 um-1) = 2.78e-14 +/- 3.60e-15
   - Gemini/NIRI.H2S1v2-1-G0220:
      - Apparent magnitude = 11.39 +/- 0.14
      - Flux (W m-2 um-1) = 1.06e-14 +/- 1.37e-15
Adding filter: Subaru/CIAO.J... [DONE]
Adding filter: Subaru/CIAO.H... [DONE]
Adding filter: Subaru/CIAO.Ks... [DONE]
Adding filter: Keck/NIRC2.NB_4.05... [DONE]
Adding filter: LBT/LMIRCam.M_77K... [DONE]
Adding object: kappa And b
   - Distance (pc) = 50.06 +/- 0.87
   - Subaru/CIAO.J:
      - Apparent magnitude = 15.86 +/- 0.21
      - Flux (W m-2 um-1) = 1.41e-15 +/- 2.74e-16
   - Subaru/CIAO.H:
      - Apparent magnitude = 14.95 +/- 0.13
      - Flux (W m-2 um-1) = 1.29e-15 +/- 1.55e-16
   - Subaru/CIAO.Ks:
      - Apparent magnitude = 14.32 +/- 0.09
      - Flux (W m-2 um-1) = 8.53e-16 +/- 7.08e-17
   - Keck/NIRC2.Lp:
      - Apparent magnitude = 13.12 +/- 0.10
      - Flux (W m-2 um-1) = 3.05e-16 +/- 2.81e-17
   - Keck/NIRC2.NB_4.05:
      - Apparent magnitude = 13.00 +/- 0.20
      - Flux (W m-2 um-1) = 2.48e-16 +/- 4.60e-17
   - LBT/LMIRCam.M_77K:
      - Apparent magnitude = 13.30 +/- 0.30
      - Flux (W m-2 um-1) = 1.02e-16 +/- 2.86e-17
Adding filter: MKO/NSFCam.Ks... [DONE]
Adding object: HD 1160 B
   - Distance (pc) = 125.90 +/- 1.20
   - MKO/NSFCam.J:
      - Apparent magnitude = 14.69 +/- 0.05
      - Flux (W m-2 um-1) = 4.13e-15 +/- 1.90e-16
   - MKO/NSFCam.H:
      - Apparent magnitude = 14.21 +/- 0.02
      - Flux (W m-2 um-1) = 2.55e-15 +/- 4.69e-17
   - MKO/NSFCam.Ks:
      - Apparent magnitude = 14.12 +/- 0.05
      - Flux (W m-2 um-1) = 1.03e-15 +/- 4.73e-17
   - Paranal/NACO.Lp:
      - Apparent magnitude = 13.60 +/- 0.10
      - Flux (W m-2 um-1) = 1.91e-16 +/- 1.76e-17
   - Keck/NIRC2.Ms:
      - Apparent magnitude = 13.81 +/- 0.24
      - Flux (W m-2 um-1) = 6.98e-17 +/- 1.56e-17
Adding filter: MKO/NSFCam.Kp... [DONE]
Adding object: ROXs 12 B
   - Distance (pc) = 144.16 +/- 1.53
   - MKO/NSFCam.J:
      - Apparent magnitude = 15.82 +/- 0.03
      - Flux (W m-2 um-1) = 1.46e-15 +/- 4.03e-17
   - MKO/NSFCam.H:
      - Apparent magnitude = 14.83 +/- 0.03
      - Flux (W m-2 um-1) = 1.44e-15 +/- 3.98e-17
   - MKO/NSFCam.Kp:
      - Apparent magnitude = 14.14 +/- 0.03
      - Flux (W m-2 um-1) = 1.07e-15 +/- 2.96e-17
Adding filter: Keck/NIRC2.J... [DONE]
Adding object: ROXs 42 Bb
   - Distance (pc) = 144.16 +/- 1.53
   - Keck/NIRC2.J:
      - Apparent magnitude = 16.91 +/- 0.11
      - Flux (W m-2 um-1) = 5.32e-16 +/- 5.40e-17
   - Keck/NIRC2.H:
      - Apparent magnitude = 15.88 +/- 0.05
      - Flux (W m-2 um-1) = 5.40e-16 +/- 2.49e-17
   - Keck/NIRC2.Ks:
      - Apparent magnitude = 15.01 +/- 0.06
      - Flux (W m-2 um-1) = 4.49e-16 +/- 2.48e-17
   - Keck/NIRC2.Lp:
      - Apparent magnitude = 13.97 +/- 0.06
      - Flux (W m-2 um-1) = 1.39e-16 +/- 7.70e-18
   - Keck/NIRC2.NB_4.05:
      - Apparent magnitude = 13.90 +/- 0.08
      - Flux (W m-2 um-1) = 1.08e-16 +/- 7.99e-18
   - Keck/NIRC2.Ms:
      - Apparent magnitude = 14.01 +/- 0.23
      - Flux (W m-2 um-1) = 5.80e-17 +/- 1.24e-17
Adding filter: Paranal/SPHERE.IRDIS_D_Y23_2... [DONE]
Adding filter: Paranal/SPHERE.IRDIS_D_Y23_3... [DONE]
Adding filter: Paranal/SPHERE.IRDIS_D_J23_3... [DONE]
Adding filter: Subaru/IRCS.Lp... [DONE]
Adding object: GJ 504 b
   - Distance (pc) = 17.54 +/- 0.08
   - Paranal/SPHERE.IRDIS_D_Y23_2:
      - Apparent magnitude = 20.98 +/- 0.20
      - Flux (W m-2 um-1) = 2.43e-17 +/- 4.51e-18
   - Paranal/SPHERE.IRDIS_D_Y23_3:
      - Apparent magnitude = 20.14 +/- 0.09
      - Flux (W m-2 um-1) = 4.45e-17 +/- 3.69e-18
   - Paranal/SPHERE.IRDIS_D_J23_3:
      - Apparent magnitude = 19.01 +/- 0.17
      - Flux (W m-2 um-1) = 7.08e-17 +/- 1.11e-17
   - Paranal/SPHERE.IRDIS_D_H23_2:
      - Apparent magnitude = 18.95 +/- 0.30
      - Flux (W m-2 um-1) = 3.44e-17 +/- 9.64e-18
   - Paranal/SPHERE.IRDIS_D_H23_3:
      - Apparent magnitude = 21.81 +/- 0.35
      - Flux (W m-2 um-1) = 2.09e-18 +/- 6.87e-19
   - Paranal/SPHERE.IRDIS_D_K12_1:
      - Apparent magnitude = 18.77 +/- 0.20
      - Flux (W m-2 um-1) = 1.51e-17 +/- 2.79e-18
   - Subaru/CIAO.J:
      - Apparent magnitude = 19.78 +/- 0.10
      - Flux (W m-2 um-1) = 3.80e-17 +/- 3.51e-18
   - Subaru/CIAO.H:
      - Apparent magnitude = 20.01 +/- 0.14
      - Flux (W m-2 um-1) = 1.22e-17 +/- 1.58e-18
   - Subaru/CIAO.Ks:
      - Apparent magnitude = 19.38 +/- 0.11
      - Flux (W m-2 um-1) = 8.07e-18 +/- 8.19e-19
   - Subaru/CIAO.CH4s:
      - Apparent magnitude = 19.58 +/- 0.13
      - Flux (W m-2 um-1) = 1.86e-17 +/- 2.23e-18
   - Subaru/IRCS.Lp:
      - Apparent magnitude = 16.70 +/- 0.17
      - Flux (W m-2 um-1) = 1.14e-17 +/- 1.80e-18
Adding filter: Gemini/GMOS-S.z... [DONE]
Adding filter: CFHT/Wircam.Y... [DONE]
Adding filter: CFHT/Wircam.J... [DONE]
Adding filter: CFHT/Wircam.H... [DONE]
Adding filter: CFHT/Wircam.Ks... [DONE]
Adding filter: WISE/WISE.W1... [DONE]
Adding filter: WISE/WISE.W2... [DONE]
Adding object: GU Psc b
   - Distance (pc) = 47.61 +/- 0.16
   - Gemini/GMOS-S.z:
      - Apparent magnitude = 21.75 +/- 0.07
      - Flux (W m-2 um-1) = 1.63e-17 +/- 1.05e-18
   - CFHT/Wircam.Y:
      - Apparent magnitude = 19.40 +/- 0.05
      - Flux (W m-2 um-1) = 1.06e-16 +/- 4.87e-18
   - CFHT/Wircam.J:
      - Apparent magnitude = 18.12 +/- 0.03
      - Flux (W m-2 um-1) = 1.72e-16 +/- 4.76e-18
   - CFHT/Wircam.H:
      - Apparent magnitude = 17.70 +/- 0.03
      - Flux (W m-2 um-1) = 1.02e-16 +/- 2.82e-18
   - CFHT/Wircam.Ks:
      - Apparent magnitude = 17.40 +/- 0.03
      - Flux (W m-2 um-1) = 4.98e-17 +/- 1.38e-18
   - WISE/WISE.W1:
      - Apparent magnitude = 17.17 +/- 0.33
      - Flux (W m-2 um-1) = 1.14e-17 +/- 3.52e-18
   - WISE/WISE.W2:
      - Apparent magnitude = 15.41 +/- 0.22
      - Flux (W m-2 um-1) = 1.70e-17 +/- 3.48e-18
Adding object: 2M0103 ABb
   - Distance (pc) = 47.20 +/- 3.10
   - Paranal/NACO.J:
      - Apparent magnitude = 15.47 +/- 0.30
      - Flux (W m-2 um-1) = 1.96e-15 +/- 5.49e-16
   - Paranal/NACO.H:
      - Apparent magnitude = 14.27 +/- 0.20
      - Flux (W m-2 um-1) = 2.28e-15 +/- 4.23e-16
   - Paranal/NACO.Ks:
      - Apparent magnitude = 13.67 +/- 0.20
      - Flux (W m-2 um-1) = 1.57e-15 +/- 2.90e-16
   - Paranal/NACO.Lp:
      - Apparent magnitude = 12.67 +/- 0.10
      - Flux (W m-2 um-1) = 4.50e-16 +/- 4.15e-17
Adding filter: Gemini/NIRI.J-G0202w... [DONE]
Adding filter: Gemini/NIRI.H-G0203w... [DONE]
Adding filter: Gemini/NIRI.K-G0204w... [DONE]
Adding filter: Gemini/NIRI.Lprime-G0207w... [DONE]
Adding object: 1RXS 1609 B
   - Distance (pc) = 139.67 +/- 1.33
   - Gemini/NIRI.J-G0202w:
      - Apparent magnitude = 17.90 +/- 0.12
      - Flux (W m-2 um-1) = 2.12e-16 +/- 2.35e-17
   - Gemini/NIRI.H-G0203w:
      - Apparent magnitude = 16.87 +/- 0.07
      - Flux (W m-2 um-1) = 2.11e-16 +/- 1.36e-17
   - Gemini/NIRI.K-G0204w:
      - Apparent magnitude = 16.17 +/- 0.18
      - Flux (W m-2 um-1) = 1.38e-16 +/- 2.30e-17
   - Gemini/NIRI.Lprime-G0207w:
      - Apparent magnitude = 14.80 +/- 0.30
      - Flux (W m-2 um-1) = 6.56e-17 +/- 1.84e-17
Adding filter: MKO/NSFCam.Lp... [DONE]
Adding object: GSC 06214 B
   - Distance (pc) = 108.84 +/- 0.51
   - MKO/NSFCam.J:
      - Apparent magnitude = 16.24 +/- 0.04
      - Flux (W m-2 um-1) = 9.91e-16 +/- 3.65e-17
   - MKO/NSFCam.H:
      - Apparent magnitude = 15.55 +/- 0.04
      - Flux (W m-2 um-1) = 7.41e-16 +/- 2.73e-17
   - MKO/NSFCam.Kp:
      - Apparent magnitude = 14.95 +/- 0.05
      - Flux (W m-2 um-1) = 5.08e-16 +/- 2.34e-17
   - MKO/NSFCam.Lp:
      - Apparent magnitude = 13.75 +/- 0.07
      - Flux (W m-2 um-1) = 1.73e-16 +/- 1.12e-17
   - LBT/LMIRCam.M_77K:
      - Apparent magnitude = 13.75 +/- 0.30
      - Flux (W m-2 um-1) = 6.75e-17 +/- 1.89e-17
Adding object: HD 72946 B
   - Distance (pc) = 25.87 +/- 0.03
   - Paranal/SPHERE.IRDIS_D_H23_2:
      - Apparent magnitude = 14.56 +/- 0.07
      - Flux (W m-2 um-1) = 1.96e-15 +/- 1.27e-16
   - Paranal/SPHERE.IRDIS_D_H23_3:
      - Apparent magnitude = 14.40 +/- 0.07
      - Flux (W m-2 um-1) = 1.93e-15 +/- 1.24e-16
Adding object: HIP 64892 B
   - Distance (pc) = 125.20 +/- 1.42
   - Paranal/SPHERE.IRDIS_D_H23_2:
      - Apparent magnitude = 14.21 +/- 0.17
      - Flux (W m-2 um-1) = 2.71e-15 +/- 4.26e-16
   - Paranal/SPHERE.IRDIS_D_H23_3:
      - Apparent magnitude = 13.94 +/- 0.17
      - Flux (W m-2 um-1) = 2.94e-15 +/- 4.63e-16
   - Paranal/SPHERE.IRDIS_D_K12_1:
      - Apparent magnitude = 13.77 +/- 0.17
      - Flux (W m-2 um-1) = 1.51e-15 +/- 2.37e-16
   - Paranal/SPHERE.IRDIS_D_K12_2:
      - Apparent magnitude = 13.45 +/- 0.19
      - Flux (W m-2 um-1) = 1.56e-15 +/- 2.75e-16
   - Paranal/NACO.Lp:
      - Apparent magnitude = 13.09 +/- 0.17
      - Flux (W m-2 um-1) = 3.06e-16 +/- 4.81e-17
Adding filter: Paranal/SPHERE.IRDIS_D_J23_2... [DONE]
Adding object: HD 13724 B
   - Distance (pc) = 43.45 +/- 0.03
   - Paranal/SPHERE.IRDIS_D_J23_2:
      - Apparent magnitude = 17.09 +/- 0.16
      - Flux (W m-2 um-1) = 5.39e-16 +/- 7.97e-17
   - Paranal/SPHERE.IRDIS_D_J23_3:
      - Apparent magnitude = 17.82 +/- 0.32
      - Flux (W m-2 um-1) = 2.12e-16 +/- 6.33e-17
   - Paranal/SPHERE.IRDIS_D_H23_2:
      - Apparent magnitude = 18.23 +/- 0.40
      - Flux (W m-2 um-1) = 6.68e-17 +/- 2.52e-17
   - Paranal/SPHERE.IRDIS_D_H23_3:
      - Apparent magnitude = 17.10 +/- 0.05
      - Flux (W m-2 um-1) = 1.60e-16 +/- 7.39e-18
   - Paranal/SPHERE.IRDIS_D_K12_1:
      - Apparent magnitude = 16.67 +/- 0.18
      - Flux (W m-2 um-1) = 1.04e-16 +/- 1.74e-17
   - Paranal/SPHERE.IRDIS_D_K12_2:
      - Apparent magnitude = 17.48 +/- 0.46
      - Flux (W m-2 um-1) = 3.82e-17 +/- 1.67e-17
Adding object: TYC 8988 B
   - Distance (pc) = 94.60 +/- 0.30
   - Paranal/SPHERE.IRDIS_D_Y23_2:
      - Apparent magnitude = 17.03 +/- 0.21
      - Flux (W m-2 um-1) = 9.25e-16 +/- 1.80e-16
   - Paranal/SPHERE.IRDIS_D_Y23_3:
      - Apparent magnitude = 16.67 +/- 0.16
      - Flux (W m-2 um-1) = 1.09e-15 +/- 1.61e-16
   - Paranal/SPHERE.IRDIS_D_J23_2:
      - Apparent magnitude = 16.27 +/- 0.08
      - Flux (W m-2 um-1) = 1.15e-15 +/- 8.46e-17
   - Paranal/SPHERE.IRDIS_D_J23_3:
      - Apparent magnitude = 15.73 +/- 0.07
      - Flux (W m-2 um-1) = 1.45e-15 +/- 9.36e-17
   - Paranal/SPHERE.IRDIS_D_H23_2:
      - Apparent magnitude = 15.11 +/- 0.08
      - Flux (W m-2 um-1) = 1.18e-15 +/- 8.73e-17
   - Paranal/SPHERE.IRDIS_D_H23_3:
      - Apparent magnitude = 14.78 +/- 0.07
      - Flux (W m-2 um-1) = 1.36e-15 +/- 8.76e-17
   - Paranal/SPHERE.IRDIS_D_K12_1:
      - Apparent magnitude = 14.44 +/- 0.04
      - Flux (W m-2 um-1) = 8.14e-16 +/- 3.00e-17
   - Paranal/SPHERE.IRDIS_D_K12_2:
      - Apparent magnitude = 14.07 +/- 0.04
      - Flux (W m-2 um-1) = 8.83e-16 +/- 3.25e-17
   - Paranal/SPHERE.IRDIS_B_J:
      - Apparent magnitude = 15.73 +/- 0.38
      - Flux (W m-2 um-1) = 1.63e-15 +/- 5.81e-16
   - Paranal/SPHERE.IRDIS_B_H:
      - Apparent magnitude = 15.87 +/- 0.38
      - Flux (W m-2 um-1) = 5.57e-16 +/- 1.99e-16
   - Paranal/SPHERE.IRDIS_B_Ks:
      - Apparent magnitude = 14.70 +/- 0.14
      - Flux (W m-2 um-1) = 5.63e-16 +/- 7.29e-17
   - Paranal/NACO.Lp:
      - Apparent magnitude = 13.30 +/- 0.08
      - Flux (W m-2 um-1) = 2.52e-16 +/- 1.86e-17
   - Paranal/NACO.Mp:
      - Apparent magnitude = 13.08 +/- 0.20
      - Flux (W m-2 um-1) = 1.27e-16 +/- 2.35e-17
Adding object: TYC 8988 C
   - Distance (pc) = 94.60 +/- 0.30
   - Paranal/SPHERE.IRDIS_D_Y23_3:
      - Apparent magnitude = 22.37 +/- 0.31
      - Flux (W m-2 um-1) = 5.70e-18 +/- 1.65e-18
   - Paranal/SPHERE.IRDIS_D_J23_2:
      - Apparent magnitude = 21.81 +/- 0.22
      - Flux (W m-2 um-1) = 6.98e-18 +/- 1.42e-18
   - Paranal/SPHERE.IRDIS_D_J23_3:
      - Apparent magnitude = 21.17 +/- 0.15
      - Flux (W m-2 um-1) = 9.68e-18 +/- 1.34e-18
   - Paranal/SPHERE.IRDIS_D_H23_2:
      - Apparent magnitude = 19.78 +/- 0.08
      - Flux (W m-2 um-1) = 1.60e-17 +/- 1.18e-18
   - Paranal/SPHERE.IRDIS_D_H23_3:
      - Apparent magnitude = 19.32 +/- 0.06
      - Flux (W m-2 um-1) = 2.08e-17 +/- 1.15e-18
   - Paranal/SPHERE.IRDIS_D_K12_1:
      - Apparent magnitude = 18.34 +/- 0.04
      - Flux (W m-2 um-1) = 2.24e-17 +/- 8.26e-19
   - Paranal/SPHERE.IRDIS_D_K12_2:
      - Apparent magnitude = 17.85 +/- 0.09
      - Flux (W m-2 um-1) = 2.72e-17 +/- 2.25e-18
   - Paranal/SPHERE.IRDIS_B_H:
      - Apparent magnitude = 19.69 +/- 0.23
      - Flux (W m-2 um-1) = 1.65e-17 +/- 3.52e-18
   - Paranal/NACO.Lp:
      - Apparent magnitude = 16.29 +/- 0.21
      - Flux (W m-2 um-1) = 1.60e-17 +/- 3.12e-18
Adding object: HD 142527 B
   - Distance (pc) = 159.26 +/- 0.72
   - Paranal/NACO.J:
      - Apparent magnitude = 10.86 +/- 0.05
      - Flux (W m-2 um-1) = 1.37e-13 +/- 6.31e-15
   - Paranal/NACO.H (2 values):
      - Apparent magnitude = 10.50 +/- 0.20
      - Flux (W m-2 um-1) = 7.35e-14 +/- 1.36e-14
      - Apparent magnitude = 10.30 +/- 0.50
      - Flux (W m-2 um-1) = 8.83e-14 +/- 4.21e-14
   - Paranal/NACO.Ks (2 values):
      - Apparent magnitude = 10.00 +/- 0.30
      - Flux (W m-2 um-1) = 4.60e-14 +/- 1.29e-14
      - Apparent magnitude = 9.80 +/- 0.10
      - Flux (W m-2 um-1) = 5.53e-14 +/- 5.10e-15
   - Paranal/NACO.Lp (2 values):
      - Apparent magnitude = 9.10 +/- 0.10
      - Flux (W m-2 um-1) = 1.21e-14 +/- 1.11e-15
      - Apparent magnitude = 9.10 +/- 0.10
      - Flux (W m-2 um-1) = 1.21e-14 +/- 1.11e-15
   - Paranal/NACO.Mp:
      - Apparent magnitude = 9.20 +/- 0.20
      - Flux (W m-2 um-1) = 4.52e-15 +/- 8.37e-16
Adding object: CS Cha B
   - Distance (pc) = 168.77 +/- 1.92
   - Paranal/SPHERE.IRDIS_B_J:
      - Apparent magnitude = 19.16 +/- 0.21
      - Flux (W m-2 um-1) = 6.91e-17 +/- 1.34e-17
   - Paranal/SPHERE.IRDIS_B_H:
      - Apparent magnitude = 17.65 +/- 0.62
      - Flux (W m-2 um-1) = 1.08e-16 +/- 6.51e-17
   - Paranal/NACO.Ks:
      - Apparent magnitude = 17.40 +/- 0.16
      - Flux (W m-2 um-1) = 5.04e-17 +/- 7.46e-18
Adding object: CT Cha B
   - Distance (pc) = 189.95 +/- 0.42
   - Paranal/NACO.J:
      - Apparent magnitude = 16.61 +/- 0.30
      - Flux (W m-2 um-1) = 6.87e-16 +/- 1.92e-16
   - Paranal/NACO.Ks (2 values):
      - Apparent magnitude = 14.95 +/- 0.30
      - Flux (W m-2 um-1) = 4.82e-16 +/- 1.35e-16
      - Apparent magnitude = 14.89 +/- 0.30
      - Flux (W m-2 um-1) = 5.09e-16 +/- 1.42e-16
Adding object: SR 12 C
   - Distance (pc) = 125.00 +/- 25.00
   - MKO/NSFCam.J:
      - Apparent magnitude = 15.93 +/- 0.03
      - Flux (W m-2 um-1) = 1.32e-15 +/- 3.64e-17
   - MKO/NSFCam.H:
      - Apparent magnitude = 15.18 +/- 0.03
      - Flux (W m-2 um-1) = 1.04e-15 +/- 2.88e-17
   - MKO/NSFCam.Ks:
      - Apparent magnitude = 14.57 +/- 0.03
      - Flux (W m-2 um-1) = 6.79e-16 +/- 1.88e-17
   - MKO/NSFCam.Lp:
      - Apparent magnitude = 13.10 +/- 0.08
      - Flux (W m-2 um-1) = 3.15e-16 +/- 2.32e-17
Adding object: DH Tau B
   - Distance (pc) = 133.45 +/- 0.45
   - Subaru/CIAO.J:
      - Apparent magnitude = 15.71 +/- 0.05
      - Flux (W m-2 um-1) = 1.61e-15 +/- 7.44e-17
   - Subaru/CIAO.H:
      - Apparent magnitude = 14.96 +/- 0.04
      - Flux (W m-2 um-1) = 1.28e-15 +/- 4.70e-17
   - Subaru/CIAO.Ks:
      - Apparent magnitude = 14.19 +/- 0.02
      - Flux (W m-2 um-1) = 9.62e-16 +/- 1.77e-17

The spectra from the SpeX Prism Spectral Library are also downloaded and added to the database. For each spectrum, the SIMBAD Astronomical Database is queried for the SIMBAD identifier. The identifier is then used to extract the distance of the object (calculated from the parallax). A NaN value is stored for the distance if the object could not be identified in the SIMBAD database so these objects are not used in the color-magnitude diagram.

[6]:
database.add_spectra('spex')
Downloading SpeX Prism Spectral Library... [DONE]
Adding filter: 2MASS/2MASS.H... [DONE]
Adding spectra... 2MASS J04414825+2534304
/Users/tomasstolker/.pyenv/versions/3.9.6/envs/species3.9/lib/python3.9/site-packages/astroquery/simbad/core.py:135: UserWarning: Warning: The script line number 3 raised an error (recorded in the `errors` attribute of the result table): Identifier not found in the database : 2MASS J21225635+3656002
  warnings.warn("Warning: The script line number %i raised "
Adding spectra... 2MASSW J1239272+551537
/Users/tomasstolker/.pyenv/versions/3.9.6/envs/species3.9/lib/python3.9/site-packages/astroquery/simbad/core.py:135: UserWarning: Warning: The script line number 3 raised an error (recorded in the `errors` attribute of the result table): Identifier not found in the database : 2MASS J14044941-3159329
  warnings.warn("Warning: The script line number %i raised "
Adding spectra... SDSS J213240.36+102949.4
/Users/tomasstolker/.pyenv/versions/3.9.6/envs/species3.9/lib/python3.9/site-packages/astroquery/simbad/core.py:135: UserWarning: Warning: The script line number 3 raised an error (recorded in the `errors` attribute of the result table): Identifier not found in the database : 2MASS J05103520-4208140
  warnings.warn("Warning: The script line number %i raised "
Adding spectra... 2MASS J03382862+0001296
/Users/tomasstolker/.pyenv/versions/3.9.6/envs/species3.9/lib/python3.9/site-packages/astroquery/simbad/core.py:135: UserWarning: Warning: The script line number 3 raised an error (recorded in the `errors` attribute of the result table): Identifier not found in the database : 2MASS J14162394+1348363
  warnings.warn("Warning: The script line number %i raised "
Adding spectra... [DONE]

The AMES-Cond model spectra (i.e. for cloudless atmospheres) are also downloaded and imported into the database.

[7]:
database.add_model('ames-cond')
Downloading AMES-Cond model spectra (150 MB)... [DONE]
Unpacking AMES-Cond model spectra (150 MB)... [DONE]
Adding AMES-Cond model spectra... [DONE]
Grid points stored in the database:
   - Teff = [ 100.  200.  300.  400.  500.  600.  700.  800.  900. 1000. 1100. 1200.
 1300. 1400. 1500. 1600. 1700. 1800. 1900. 2000. 2100. 2200. 2300. 2400.
 2500. 2600. 2700. 2800. 2900. 3000. 3100. 3200. 3300. 3400. 3500. 3600.
 3700. 3800. 3900. 4000. 4100. 4200. 4300. 4400. 4500. 4600. 4700. 4800.
 5000. 5100. 5200. 5300. 5400. 5500. 5600. 5700. 5800. 5900. 6000. 6100.
 6200. 6300. 6400. 6500. 6600.]
   - log(g) = [2.5 3.  3.5 4.  4.5 5.  5.5]
Number of grid points per parameter:
   - teff: 65
   - logg: 7
/Users/tomasstolker/applications/species/species/util/data_util.py:282: RuntimeWarning: divide by zero encountered in log10
  flux = np.log10(flux)
Fix missing grid points with a linear interpolation:
   - teff = 200.0, logg = 5.5
   - teff = 900.0, logg = 2.5
   - teff = 4700.0, logg = 4.5
   - teff = 4700.0, logg = 5.0
   - teff = 4700.0, logg = 5.5
   - teff = 4800.0, logg = 2.5
   - teff = 4800.0, logg = 3.0
   - teff = 4800.0, logg = 4.5
   - teff = 6100.0, logg = 2.5
   - teff = 6100.0, logg = 3.0
   - teff = 6100.0, logg = 3.5
   - teff = 6100.0, logg = 4.0
   - teff = 6100.0, logg = 4.5
   - teff = 6100.0, logg = 5.5
   - teff = 6200.0, logg = 3.5
   - teff = 6200.0, logg = 4.0
Number of stored grid points: 455
Number of interpolated grid points: 16
Number of missing grid points: 0

And also the AMES-Dusty model spectra (i.e. for dusty atmospheres) are downloaded and imported into the database.

[8]:
database.add_model('ames-dusty')
Downloading AMES-Dusty model spectra (59 MB)... [DONE]
Unpacking AMES-Dusty model spectra (59 MB)... [DONE]
Adding AMES-Dusty model spectra... [DONE]
Grid points stored in the database:
   - Teff = [ 500.  600.  700.  800.  900. 1000. 1100. 1200. 1300. 1400. 1500. 1600.
 1700. 1800. 1900. 2000. 2100. 2200. 2300. 2400. 2500. 2600. 2700. 2800.
 2900. 3000. 3100. 3200. 3300. 3400. 3500. 3600. 3700. 3800. 3900. 4000.]
   - log(g) = [3.5 4.  4.5 5.  5.5 6. ]
Number of grid points per parameter:
   - teff: 36
   - logg: 6
Fix missing grid points with a linear interpolation:
   - teff = 900.0, logg = 6.0
   - teff = 1200.0, logg = 5.5
   - teff = 2100.0, logg = 3.5
   - teff = 2100.0, logg = 4.5
   - teff = 2200.0, logg = 3.5
   - teff = 2400.0, logg = 5.0
   - teff = 3100.0, logg = 3.5
   - teff = 3200.0, logg = 3.5
   - teff = 3300.0, logg = 3.5
   - teff = 3400.0, logg = 3.5
   - teff = 3500.0, logg = 3.5
   - teff = 3600.0, logg = 3.5
   - teff = 3700.0, logg = 3.5
   - teff = 3800.0, logg = 3.5
   - teff = 3900.0, logg = 3.5
   - teff = 3900.0, logg = 6.0
   - teff = 4000.0, logg = 3.5
   - teff = 4000.0, logg = 4.0
   - teff = 4000.0, logg = 4.5
   - teff = 4000.0, logg = 5.5
   - teff = 4000.0, logg = 6.0
Could not interpolate 15 grid points so storing zeros instead. [WARNING]
The grid points that are missing:
   - teff = 3100.0, logg = 3.5
   - teff = 3200.0, logg = 3.5
   - teff = 3300.0, logg = 3.5
   - teff = 3400.0, logg = 3.5
   - teff = 3500.0, logg = 3.5
   - teff = 3600.0, logg = 3.5
   - teff = 3700.0, logg = 3.5
   - teff = 3800.0, logg = 3.5
   - teff = 3900.0, logg = 3.5
   - teff = 3900.0, logg = 6.0
   - teff = 4000.0, logg = 3.5
   - teff = 4000.0, logg = 4.0
   - teff = 4000.0, logg = 4.5
   - teff = 4000.0, logg = 5.5
   - teff = 4000.0, logg = 6.0
Number of stored grid points: 216
Number of interpolated grid points: 6
Number of missing grid points: 15

The AMES-Cond and AMES-Dusty isochones are now downloaded. We note that the evolutionary data of these isochrones are identical and only the magnitudes are different, but for clarity we consider the isochrones from both AMES-Cond and AMES-Dusty. Since the SSL certificate of the website has expired, we need a workaround to use the urlretrieve function without causing an error.

[9]:
import ssl
ssl._create_default_https_context = ssl._create_unverified_context
[10]:
urllib.request.urlretrieve('https://phoenix.ens-lyon.fr/Grids/AMES-Cond/ISOCHRONES/model.AMES-Cond-2000.M-0.0.NaCo.Vega',
                           'data/model.AMES-Cond-2000.M-0.0.NaCo.Vega')
[10]:
('data/model.AMES-Cond-2000.M-0.0.NaCo.Vega',
 <http.client.HTTPMessage at 0x14ae390a0>)
[11]:
urllib.request.urlretrieve('https://phoenix.ens-lyon.fr/Grids/AMES-Dusty/ISOCHRONES/model.AMES-dusty.M-0.0.NaCo.Vega',
                           'data/model.AMES-dusty.M-0.0.NaCo.Vega')
[11]:
('data/model.AMES-dusty.M-0.0.NaCo.Vega',
 <http.client.HTTPMessage at 0x14ae39460>)

The isochrone data are added to the database with the add_isochrones method.

[12]:
database.add_isochrones(filename='data/model.AMES-Cond-2000.M-0.0.NaCo.Vega',
                        tag='iso_ames-cond')
Adding isochrones: iso_ames-cond... [DONE]
[13]:
database.add_isochrones(filename='data/model.AMES-dusty.M-0.0.NaCo.Vega',
                        tag='iso_ames-dusty')
Adding isochrones: iso_ames-dusty... [DONE]

Synthetic colors and magnitudes

For the narrowband filters of SPHERE, the synthetic colors and magnitudes are computed from the SpeX spectra. We start by creating an instance of ReadColorMagnitude.

[14]:
read_color = species.ReadColorMagnitude(library='spex',
                                        filters_color=('Paranal/SPHERE.IRDIS_D_H23_2', 'Paranal/SPHERE.IRDIS_D_H23_3'),
                                        filter_mag='Paranal/SPHERE.IRDIS_D_H23_2')

All spectra are used and the colors and magnitudes are stored in a ColorMagObject.

[15]:
colorbox = read_color.get_color_magnitude(object_type=None)
/Users/tomasstolker/applications/species/species/analysis/photometry.py:230: UserWarning: Spectum contains NaN so can not calculate the error.
  warnings.warn('Spectum contains NaN so can not calculate the error.')
/Users/tomasstolker/applications/species/species/analysis/photometry.py:304: UserWarning: Spectum contains NaN so can not calculate the error.
  warnings.warn('Spectum contains NaN so can not calculate the error.')

Next, the isochrone data is extracted by first creating and instance of ReadIsochrone.

[16]:
read_iso_cond = species.ReadIsochrone(tag='iso_ames-cond')
read_iso_dusty = species.ReadIsochrone(tag='iso_ames-dusty')

To interpolate the isochrones and compute synthetic photometry, we use the get_color_magnitude method and give the grid of model spectra (i.e. that was stored with add_model) as argument of model. The method will return a ColorMagBox object.

[17]:
modelcolor1 = read_iso_cond.get_color_magnitude(age=age,
                                                masses=masses,
                                                model='ames-cond',
                                                filters_color=('Paranal/SPHERE.IRDIS_D_H23_2', 'Paranal/SPHERE.IRDIS_D_H23_3'),
                                                filter_mag='Paranal/SPHERE.IRDIS_D_H23_2')
[18]:
modelcolor2 = read_iso_dusty.get_color_magnitude(age=age,
                                                 masses=masses,
                                                 model='ames-dusty',
                                                 filters_color=('Paranal/SPHERE.IRDIS_D_H23_2', 'Paranal/SPHERE.IRDIS_D_H23_3'),
                                                 filter_mag='Paranal/SPHERE.IRDIS_D_H23_2')
/Users/tomasstolker/applications/species/species/read/read_isochrone.py:247: UserWarning: The value of logg is 2.917784581524866, which is below the lower bound of the model grid (3.5). Setting the magnitudes to NaN for the following isochrone sample: {'teff': 501.3728693602355, 'logg': 2.917784581524866, 'mass': 1.0, 'distance': 10.0}.
  warnings.warn(f'The value of {item_bounds} is '
/Users/tomasstolker/applications/species/species/read/read_isochrone.py:247: UserWarning: The value of logg is 2.972779343287681, which is below the lower bound of the model grid (3.5). Setting the magnitudes to NaN for the following isochrone sample: {'teff': 502.51561905255875, 'logg': 2.972779343287681, 'mass': 1.0722672220103233, 'distance': 10.0}.
  warnings.warn(f'The value of {item_bounds} is '
/Users/tomasstolker/applications/species/species/read/read_isochrone.py:247: UserWarning: The value of logg is 3.035244115774989, which is below the lower bound of the model grid (3.5). Setting the magnitudes to NaN for the following isochrone sample: {'teff': 504.9110749694396, 'logg': 3.035244115774989, 'mass': 1.1497569953977358, 'distance': 10.0}.
  warnings.warn(f'The value of {item_bounds} is '
/Users/tomasstolker/applications/species/species/read/read_isochrone.py:247: UserWarning: The value of logg is 3.1022230438434617, which is below the lower bound of the model grid (3.5). Setting the magnitudes to NaN for the following isochrone sample: {'teff': 507.4796438308816, 'logg': 3.1022230438434617, 'mass': 1.2328467394420661, 'distance': 10.0}.
  warnings.warn(f'The value of {item_bounds} is '
/Users/tomasstolker/applications/species/species/read/read_isochrone.py:247: UserWarning: The value of logg is 3.2067743991749165, which is below the lower bound of the model grid (3.5). Setting the magnitudes to NaN for the following isochrone sample: {'teff': 519.9830583907551, 'logg': 3.2067743991749165, 'mass': 1.321941148466029, 'distance': 10.0}.
  warnings.warn(f'The value of {item_bounds} is '
/Users/tomasstolker/applications/species/species/read/read_isochrone.py:247: UserWarning: The value of logg is 3.3063735241310974, which is below the lower bound of the model grid (3.5). Setting the magnitudes to NaN for the following isochrone sample: {'teff': 533.9699983677936, 'logg': 3.3063735241310974, 'mass': 1.4174741629268053, 'distance': 10.0}.
  warnings.warn(f'The value of {item_bounds} is '
/Users/tomasstolker/applications/species/species/read/read_isochrone.py:247: UserWarning: The value of logg is 3.338631698699009, which is below the lower bound of the model grid (3.5). Setting the magnitudes to NaN for the following isochrone sample: {'teff': 543.9407068706029, 'logg': 3.338631698699009, 'mass': 1.5199110829529336, 'distance': 10.0}.
  warnings.warn(f'The value of {item_bounds} is '
/Users/tomasstolker/applications/species/species/read/read_isochrone.py:247: UserWarning: The value of logg is 3.3721167739343065, which is below the lower bound of the model grid (3.5). Setting the magnitudes to NaN for the following isochrone sample: {'teff': 556.1780019724513, 'logg': 3.3721167739343065, 'mass': 1.6297508346206442, 'distance': 10.0}.
  warnings.warn(f'The value of {item_bounds} is '
/Users/tomasstolker/applications/species/species/read/read_isochrone.py:247: UserWarning: The value of logg is 3.406958013871118, which is below the lower bound of the model grid (3.5). Setting the magnitudes to NaN for the following isochrone sample: {'teff': 570.788844526598, 'logg': 3.406958013871118, 'mass': 1.7475284000076838, 'distance': 10.0}.
  warnings.warn(f'The value of {item_bounds} is '
/Users/tomasstolker/applications/species/species/read/read_isochrone.py:247: UserWarning: The value of logg is 3.444317133429558, which is below the lower bound of the model grid (3.5). Setting the magnitudes to NaN for the following isochrone sample: {'teff': 586.455572083363, 'logg': 3.444317133429558, 'mass': 1.8738174228603839, 'distance': 10.0}.
  warnings.warn(f'The value of {item_bounds} is '
/Users/tomasstolker/applications/species/species/read/read_isochrone.py:247: UserWarning: The value of logg is 3.484376092775238, which is below the lower bound of the model grid (3.5). Setting the magnitudes to NaN for the following isochrone sample: {'teff': 603.2544905186481, 'logg': 3.484376092775238, 'mass': 2.0092330025650473, 'distance': 10.0}.
  warnings.warn(f'The value of {item_bounds} is '

Some warnings are printed for temperatures (i.e. masses) that are outside the grid of isochrone data and also when the log(g) from the isochrone data is outside the log(g) range of the grid with model spectra. In those case, the magnitudes and colors are set to NaN and will be ignored in the color-magnitude diagram.

Selecting directly imaged planets

Let’s have a look at the available companion data that were added to the database with add_companion.

[19]:
database.list_companions()
Object name = beta Pic b
Distance (pc) = 19.75 +/- 0.13
Magellan/VisAO.Ys (mag) = 15.53 +/- 0.34
Paranal/NACO.J (mag) = 14.11 +/- 0.21
Gemini/NICI.ED286 (mag) = 13.18 +/- 0.15
Paranal/NACO.H (mag) = 13.32 +/- 0.14
Paranal/NACO.Ks (mag) = 12.64 +/- 0.11
Paranal/NACO.NB374 (mag) = 11.25 +/- 0.23
Paranal/NACO.Lp (mag) = 11.3 +/- 0.06
Paranal/NACO.NB405 (mag) = 10.98 +/- 0.05
Paranal/NACO.Mp (mag) = 11.1 +/- 0.12
GPI_YJHK spectrum from Chilcote et al. 2017, AJ, 153, 182
GRAVITY spectrum from Gravity Collaboration et al. 2020, A&A, 633, 110

Object name = beta Pic c
Distance (pc) = 19.75 +/- 0.13
MKO/NSFCam.K (mag) = 14.3 +/- 0.1

Object name = HIP 65426 b
Distance (pc) = 109.21 +/- 0.75
Paranal/SPHERE.IRDIS_D_H23_2 (mag) = 17.94 +/- 0.05
Paranal/SPHERE.IRDIS_D_H23_3 (mag) = 17.58 +/- 0.06
Paranal/SPHERE.IRDIS_D_K12_1 (mag) = 17.01 +/- 0.09
Paranal/SPHERE.IRDIS_D_K12_2 (mag) = 16.79 +/- 0.09
Paranal/NACO.Lp (mag) = 15.33 +/- 0.12
Paranal/NACO.NB405 (mag) = 15.23 +/- 0.22
Paranal/NACO.Mp (mag) = 14.65 +/- 0.29

Object name = 51 Eri b
Distance (pc) = 29.78 +/- 0.12
MKO/NSFCam.J (mag) = 19.04 +/- 0.4
MKO/NSFCam.H (mag) = 18.99 +/- 0.21
MKO/NSFCam.K (mag) = 18.67 +/- 0.19
Paranal/SPHERE.IRDIS_B_H (mag) = 19.45 +/- 0.29
Paranal/SPHERE.IRDIS_D_H23_2 (mag) = 18.41 +/- 0.26
Paranal/SPHERE.IRDIS_D_K12_1 (mag) = 17.55 +/- 0.14
Keck/NIRC2.Lp (mag) = 16.2 +/- 0.11
Keck/NIRC2.Ms (mag) = 16.1 +/- 0.5
SPHERE_YJH spectrum from Samland et al. 2017, A&A, 603, 57

Object name = HR 8799 b
Distance (pc) = 41.29 +/- 0.15
Subaru/CIAO.z (mag) = 21.22 +/- 0.29
Paranal/SPHERE.IRDIS_B_J (mag) = 19.78 +/- 0.09
Keck/NIRC2.H (mag) = 18.05 +/- 0.09
Paranal/SPHERE.IRDIS_D_H23_2 (mag) = 18.08 +/- 0.14
Paranal/SPHERE.IRDIS_D_H23_3 (mag) = 17.78 +/- 0.1
Keck/NIRC2.Ks (mag) = 17.03 +/- 0.08
Paranal/SPHERE.IRDIS_D_K12_1 (mag) = 17.15 +/- 0.06
Paranal/SPHERE.IRDIS_D_K12_2 (mag) = 16.97 +/- 0.09
Paranal/NACO.Lp (mag) = 15.52 +/- 0.1
Paranal/NACO.NB405 (mag) = 14.82 +/- 0.18
Keck/NIRC2.Ms (mag) = 16.05 +/- 0.3

Object name = HR 8799 c
Distance (pc) = 41.29 +/- 0.15
Paranal/SPHERE.IRDIS_B_J (mag) = 18.6 +/- 0.13
Keck/NIRC2.H (mag) = 17.06 +/- 0.13
Paranal/SPHERE.IRDIS_D_H23_2 (mag) = 17.09 +/- 0.12
Paranal/SPHERE.IRDIS_D_H23_3 (mag) = 16.78 +/- 0.1
Keck/NIRC2.Ks (mag) = 16.11 +/- 0.08
Paranal/SPHERE.IRDIS_D_K12_1 (mag) = 16.19 +/- 0.05
Paranal/SPHERE.IRDIS_D_K12_2 (mag) = 15.86 +/- 0.07
Paranal/NACO.Lp (mag) = 14.65 +/- 0.11
Paranal/NACO.NB405 (mag) = 13.97 +/- 0.11
Keck/NIRC2.Ms (mag) = 15.03 +/- 0.14

Object name = HR 8799 d
Distance (pc) = 41.29 +/- 0.15
Paranal/SPHERE.IRDIS_B_J (mag) = 18.59 +/- 0.37
Keck/NIRC2.H (mag) = 16.71 +/- 0.24
Paranal/SPHERE.IRDIS_D_H23_2 (mag) = 17.02 +/- 0.17
Paranal/SPHERE.IRDIS_D_H23_3 (mag) = 16.85 +/- 0.16
Keck/NIRC2.Ks (mag) = 16.09 +/- 0.12
Paranal/SPHERE.IRDIS_D_K12_1 (mag) = 16.2 +/- 0.07
Paranal/SPHERE.IRDIS_D_K12_2 (mag) = 15.84 +/- 0.1
Paranal/NACO.Lp (mag) = 14.55 +/- 0.14
Paranal/NACO.NB405 (mag) = 13.87 +/- 0.15
Keck/NIRC2.Ms (mag) = 14.65 +/- 0.35

Object name = HR 8799 e
Distance (pc) = 41.29 +/- 0.15
Paranal/SPHERE.IRDIS_B_J (mag) = 18.4 +/- 0.21
Paranal/SPHERE.IRDIS_D_H23_2 (mag) = 16.91 +/- 0.2
Paranal/SPHERE.IRDIS_D_H23_3 (mag) = 16.68 +/- 0.21
Keck/NIRC2.Ks (mag) = 15.91 +/- 0.22
Paranal/SPHERE.IRDIS_D_K12_1 (mag) = 16.12 +/- 0.1
Paranal/SPHERE.IRDIS_D_K12_2 (mag) = 15.82 +/- 0.11
Paranal/NACO.Lp (mag) = 14.49 +/- 0.21
Paranal/NACO.NB405 (mag) = 13.72 +/- 0.2
SPHERE_YJH spectrum from Zurlo et al. 2016, A&A, 587, 57

Object name = HD 95086 b
Distance (pc) = 86.44 +/- 0.24
Gemini/GPI.H (mag) = 20.51 +/- 0.25
Gemini/GPI.K1 (mag) = 18.99 +/- 0.2
Paranal/NACO.Lp (mag) = 16.27 +/- 0.19

Object name = PDS 70 b
Distance (pc) = 113.43 +/- 0.52
Paranal/SPHERE.IRDIS_D_H23_2 (mag) = 18.12 +/- 0.21
Paranal/SPHERE.IRDIS_D_H23_3 (mag) = 17.97 +/- 0.18
Paranal/SPHERE.IRDIS_D_K12_1 (mag) = 16.66 +/- 0.04
Paranal/SPHERE.IRDIS_D_K12_2 (mag) = 16.37 +/- 0.06
MKO/NSFCam.J (mag) = 20.04 +/- 0.09
MKO/NSFCam.H (mag) = 18.24 +/- 0.04
Paranal/NACO.Lp (mag) = 14.68 +/- 0.22
Paranal/NACO.NB405 (mag) = 14.68 +/- 0.27
Paranal/NACO.Mp (mag) = 13.8 +/- 0.27
Keck/NIRC2.Lp (mag) = 14.64 +/- 0.18
SPHERE_YJH spectrum from Müller et al. 2018, A&A, 617, 2

Object name = PDS 70 c
Distance (pc) = 113.43 +/- 0.52
Paranal/NACO.NB405 (mag) = 14.91 +/- 0.35
Keck/NIRC2.Lp (mag) = 15.5 +/- 0.46

Object name = 2M1207 B
Distance (pc) = 64.42 +/- 0.65
HST/NICMOS1.F090M (mag) = 22.58 +/- 0.35
HST/NICMOS1.F110M (mag) = 20.61 +/- 0.15
HST/NICMOS1.F145M (mag) = 19.05 +/- 0.03
HST/NICMOS1.F160W (mag) = 18.27 +/- 0.02
Paranal/NACO.J (mag) = 20.0 +/- 0.2
Paranal/NACO.H (mag) = 18.09 +/- 0.21
Paranal/NACO.Ks (mag) = 16.93 +/- 0.11
Paranal/NACO.Lp (mag) = 15.28 +/- 0.14

Object name = AB Pic B
Distance (pc) = 50.12 +/- 0.07
Paranal/NACO.J (mag) = 16.18 +/- 0.1
Paranal/NACO.H (mag) = 14.69 +/- 0.1
Paranal/NACO.Ks (mag) = 14.14 +/- 0.08

Object name = HD 206893 B
Distance (pc) = 40.81 +/- 0.11
Paranal/SPHERE.IRDIS_B_H (mag) = 16.79 +/- 0.06
Paranal/SPHERE.IRDIS_D_K12_1 (mag) = 15.2 +/- 0.1
Paranal/SPHERE.IRDIS_D_K12_2 (mag) = 14.88 +/- 0.09
Paranal/NACO.Lp (mag) = 13.79 +/- 0.31
Paranal/NACO.NB405 (mag) = 13.16 +/- 0.34
Paranal/NACO.Mp (mag) = 12.77 +/- 0.27
SPHERE_YJH spectrum from Delorme et al. 2017, A&A, 608, 79

Object name = RZ Psc B
Distance (pc) = 195.86 +/- 4.03
Paranal/SPHERE.IRDIS_B_H (mag) = (13.71, 0.14) +/- (13.85, 0.26)
Paranal/SPHERE.IRDIS_B_Ks (mag) = 13.51 +/- 0.2

Object name = GQ Lup B
Distance (pc) = 154.1 +/- 0.69
HST/WFPC2-PC.F606W (mag) = 19.19 +/- 0.07
HST/WFPC2-PC.F814W (mag) = 17.67 +/- 0.05
HST/NICMOS2.F171M (mag) = 13.84 +/- 0.13
HST/NICMOS2.F190N (mag) = 14.08 +/- 0.2
HST/NICMOS2.F215N (mag) = 13.4 +/- 0.15
Magellan/VisAO.ip (mag) = 18.89 +/- 0.24
Magellan/VisAO.zp (mag) = 16.4 +/- 0.1
Magellan/VisAO.Ys (mag) = 15.88 +/- 0.1
MKO/NSFCam.H (mag) = 14.01 +/- 0.13
Paranal/NACO.NB405 (mag) = 12.29 +/- 0.06
Paranal/NACO.Mp (mag) = 11.97 +/- 0.08
Paranal/NACO.Ks (mag) = (13.474, 0.031) +/- (13.386, 0.032)
Subaru/CIAO.CH4s (mag) = 13.76 +/- 0.26
Subaru/CIAO.K (mag) = 13.37 +/- 0.12
Subaru/CIAO.Lp (mag) = 12.44 +/- 0.22

Object name = PZ Tel B
Distance (pc) = 47.13 +/- 0.13
Paranal/SPHERE.ZIMPOL_R_PRIM (mag) = 17.84 +/- 0.31
Paranal/SPHERE.ZIMPOL_I_PRIM (mag) = 15.16 +/- 0.12
Paranal/SPHERE.IRDIS_D_H23_2 (mag) = 11.78 +/- 0.19
Paranal/SPHERE.IRDIS_D_H23_3 (mag) = 11.65 +/- 0.19
Paranal/SPHERE.IRDIS_D_K12_1 (mag) = 11.56 +/- 0.09
Paranal/SPHERE.IRDIS_D_K12_2 (mag) = 11.29 +/- 0.1
Paranal/NACO.J (mag) = 12.47 +/- 0.2
Paranal/NACO.H (mag) = 11.93 +/- 0.14
Paranal/NACO.Ks (mag) = 11.53 +/- 0.07
Paranal/NACO.Lp (mag) = 11.04 +/- 0.22
Paranal/NACO.NB405 (mag) = 10.94 +/- 0.07
Paranal/NACO.Mp (mag) = 10.93 +/- 0.03
Gemini/NICI.ED286 (mag) = 11.68 +/- 0.14
Gemini/NIRI.H2S1v2-1-G0220 (mag) = 11.39 +/- 0.14

Object name = kappa And b
Distance (pc) = 50.06 +/- 0.87
Subaru/CIAO.J (mag) = 15.86 +/- 0.21
Subaru/CIAO.H (mag) = 14.95 +/- 0.13
Subaru/CIAO.Ks (mag) = 14.32 +/- 0.09
Keck/NIRC2.Lp (mag) = 13.12 +/- 0.1
Keck/NIRC2.NB_4.05 (mag) = 13.0 +/- 0.2
LBT/LMIRCam.M_77K (mag) = 13.3 +/- 0.3

Object name = HD 1160 B
Distance (pc) = 125.9 +/- 1.2
MKO/NSFCam.J (mag) = 14.69 +/- 0.05
MKO/NSFCam.H (mag) = 14.21 +/- 0.02
MKO/NSFCam.Ks (mag) = 14.12 +/- 0.05
Paranal/NACO.Lp (mag) = 13.6 +/- 0.1
Keck/NIRC2.Ms (mag) = 13.81 +/- 0.24

Object name = ROXs 12 B
Distance (pc) = 144.16 +/- 1.53
MKO/NSFCam.J (mag) = 15.82 +/- 0.03
MKO/NSFCam.H (mag) = 14.83 +/- 0.03
MKO/NSFCam.Kp (mag) = 14.14 +/- 0.03

Object name = ROXs 42 Bb
Distance (pc) = 144.16 +/- 1.53
Keck/NIRC2.J (mag) = 16.91 +/- 0.11
Keck/NIRC2.H (mag) = 15.88 +/- 0.05
Keck/NIRC2.Ks (mag) = 15.01 +/- 0.06
Keck/NIRC2.Lp (mag) = 13.97 +/- 0.06
Keck/NIRC2.NB_4.05 (mag) = 13.9 +/- 0.08
Keck/NIRC2.Ms (mag) = 14.01 +/- 0.23

Object name = GJ 504 b
Distance (pc) = 17.54 +/- 0.08
Paranal/SPHERE.IRDIS_D_Y23_2 (mag) = 20.98 +/- 0.2
Paranal/SPHERE.IRDIS_D_Y23_3 (mag) = 20.14 +/- 0.09
Paranal/SPHERE.IRDIS_D_J23_3 (mag) = 19.01 +/- 0.17
Paranal/SPHERE.IRDIS_D_H23_2 (mag) = 18.95 +/- 0.3
Paranal/SPHERE.IRDIS_D_H23_3 (mag) = 21.81 +/- 0.35
Paranal/SPHERE.IRDIS_D_K12_1 (mag) = 18.77 +/- 0.2
Subaru/CIAO.J (mag) = 19.78 +/- 0.1
Subaru/CIAO.H (mag) = 20.01 +/- 0.14
Subaru/CIAO.Ks (mag) = 19.38 +/- 0.11
Subaru/CIAO.CH4s (mag) = 19.58 +/- 0.13
Subaru/IRCS.Lp (mag) = 16.7 +/- 0.17

Object name = GU Psc b
Distance (pc) = 47.61 +/- 0.16
Gemini/GMOS-S.z (mag) = 21.75 +/- 0.07
CFHT/Wircam.Y (mag) = 19.4 +/- 0.05
CFHT/Wircam.J (mag) = 18.12 +/- 0.03
CFHT/Wircam.H (mag) = 17.7 +/- 0.03
CFHT/Wircam.Ks (mag) = 17.4 +/- 0.03
WISE/WISE.W1 (mag) = 17.17 +/- 0.33
WISE/WISE.W2 (mag) = 15.41 +/- 0.22

Object name = 2M0103 ABb
Distance (pc) = 47.2 +/- 3.1
Paranal/NACO.J (mag) = 15.47 +/- 0.3
Paranal/NACO.H (mag) = 14.27 +/- 0.2
Paranal/NACO.Ks (mag) = 13.67 +/- 0.2
Paranal/NACO.Lp (mag) = 12.67 +/- 0.1

Object name = 1RXS 1609 B
Distance (pc) = 139.67 +/- 1.33
Gemini/NIRI.J-G0202w (mag) = 17.9 +/- 0.12
Gemini/NIRI.H-G0203w (mag) = 16.87 +/- 0.07
Gemini/NIRI.K-G0204w (mag) = 16.17 +/- 0.18
Gemini/NIRI.Lprime-G0207w (mag) = 14.8 +/- 0.3

Object name = GSC 06214 B
Distance (pc) = 108.84 +/- 0.51
MKO/NSFCam.J (mag) = 16.24 +/- 0.04
MKO/NSFCam.H (mag) = 15.55 +/- 0.04
MKO/NSFCam.Kp (mag) = 14.95 +/- 0.05
MKO/NSFCam.Lp (mag) = 13.75 +/- 0.07
LBT/LMIRCam.M_77K (mag) = 13.75 +/- 0.3

Object name = HD 72946 B
Distance (pc) = 25.87 +/- 0.03
Paranal/SPHERE.IRDIS_D_H23_2 (mag) = 14.56 +/- 0.07
Paranal/SPHERE.IRDIS_D_H23_3 (mag) = 14.4 +/- 0.07

Object name = HIP 64892 B
Distance (pc) = 125.2 +/- 1.42
Paranal/SPHERE.IRDIS_D_H23_2 (mag) = 14.21 +/- 0.17
Paranal/SPHERE.IRDIS_D_H23_3 (mag) = 13.94 +/- 0.17
Paranal/SPHERE.IRDIS_D_K12_1 (mag) = 13.77 +/- 0.17
Paranal/SPHERE.IRDIS_D_K12_2 (mag) = 13.45 +/- 0.19
Paranal/NACO.Lp (mag) = 13.09 +/- 0.17

Object name = HD 13724 B
Distance (pc) = 43.45 +/- 0.03
Paranal/SPHERE.IRDIS_D_J23_2 (mag) = 17.09 +/- 0.16
Paranal/SPHERE.IRDIS_D_J23_3 (mag) = 17.82 +/- 0.32
Paranal/SPHERE.IRDIS_D_H23_2 (mag) = 18.23 +/- 0.4
Paranal/SPHERE.IRDIS_D_H23_3 (mag) = 17.1 +/- 0.05
Paranal/SPHERE.IRDIS_D_K12_1 (mag) = 16.67 +/- 0.18
Paranal/SPHERE.IRDIS_D_K12_2 (mag) = 17.48 +/- 0.46

Object name = TYC 8988 B
Distance (pc) = 94.6 +/- 0.3
Paranal/SPHERE.IRDIS_D_Y23_2 (mag) = 17.03 +/- 0.21
Paranal/SPHERE.IRDIS_D_Y23_3 (mag) = 16.67 +/- 0.16
Paranal/SPHERE.IRDIS_D_J23_2 (mag) = 16.27 +/- 0.08
Paranal/SPHERE.IRDIS_D_J23_3 (mag) = 15.73 +/- 0.07
Paranal/SPHERE.IRDIS_D_H23_2 (mag) = 15.11 +/- 0.08
Paranal/SPHERE.IRDIS_D_H23_3 (mag) = 14.78 +/- 0.07
Paranal/SPHERE.IRDIS_D_K12_1 (mag) = 14.44 +/- 0.04
Paranal/SPHERE.IRDIS_D_K12_2 (mag) = 14.07 +/- 0.04
Paranal/SPHERE.IRDIS_B_J (mag) = 15.73 +/- 0.38
Paranal/SPHERE.IRDIS_B_H (mag) = 15.87 +/- 0.38
Paranal/SPHERE.IRDIS_B_Ks (mag) = 14.7 +/- 0.14
Paranal/NACO.Lp (mag) = 13.3 +/- 0.08
Paranal/NACO.Mp (mag) = 13.08 +/- 0.2

Object name = TYC 8988 C
Distance (pc) = 94.6 +/- 0.3
Paranal/SPHERE.IRDIS_D_Y23_3 (mag) = 22.37 +/- 0.31
Paranal/SPHERE.IRDIS_D_J23_2 (mag) = 21.81 +/- 0.22
Paranal/SPHERE.IRDIS_D_J23_3 (mag) = 21.17 +/- 0.15
Paranal/SPHERE.IRDIS_D_H23_2 (mag) = 19.78 +/- 0.08
Paranal/SPHERE.IRDIS_D_H23_3 (mag) = 19.32 +/- 0.06
Paranal/SPHERE.IRDIS_D_K12_1 (mag) = 18.34 +/- 0.04
Paranal/SPHERE.IRDIS_D_K12_2 (mag) = 17.85 +/- 0.09
Paranal/SPHERE.IRDIS_B_H (mag) = 19.69 +/- 0.23
Paranal/NACO.Lp (mag) = 16.29 +/- 0.21

Object name = HD 142527 B
Distance (pc) = 159.26 +/- 0.72
Paranal/NACO.J (mag) = 10.86 +/- 0.05
Paranal/NACO.H (mag) = (10.5, 0.2) +/- (10.3, 0.5)
Paranal/NACO.Ks (mag) = (10.0, 0.3) +/- (9.8, 0.1)
Paranal/NACO.Lp (mag) = (9.1, 0.1) +/- (9.1, 0.1)
Paranal/NACO.Mp (mag) = 9.2 +/- 0.2

Object name = CS Cha B
Distance (pc) = 168.77 +/- 1.92
Paranal/SPHERE.IRDIS_B_J (mag) = 19.16 +/- 0.21
Paranal/SPHERE.IRDIS_B_H (mag) = 17.65 +/- 0.62
Paranal/NACO.Ks (mag) = 17.4 +/- 0.16

Object name = CT Cha B
Distance (pc) = 189.95 +/- 0.42
Paranal/NACO.J (mag) = 16.61 +/- 0.3
Paranal/NACO.Ks (mag) = (14.95, 0.3) +/- (14.89, 0.3)

Object name = SR 12 C
Distance (pc) = 125.0 +/- 25.0
MKO/NSFCam.J (mag) = 15.93 +/- 0.03
MKO/NSFCam.H (mag) = 15.18 +/- 0.03
MKO/NSFCam.Ks (mag) = 14.57 +/- 0.03
MKO/NSFCam.Lp (mag) = 13.1 +/- 0.08

Object name = DH Tau B
Distance (pc) = 133.45 +/- 0.45
Subaru/CIAO.J (mag) = 15.71 +/- 0.05
Subaru/CIAO.H (mag) = 14.96 +/- 0.04
Subaru/CIAO.Ks (mag) = 14.19 +/- 0.02

We select a few of these directly imaged planets and brown dwarfs for which SPHERE H2 and H3 photometry are available in the database. For each companion, a tuple needs to be provided with the name (as used in the database), the filter names for the color, and the filter name for the magnitude.

[20]:
objects = [('HIP 65426 b', 'Paranal/SPHERE.IRDIS_D_H23_2', 'Paranal/SPHERE.IRDIS_D_H23_3', 'Paranal/SPHERE.IRDIS_D_H23_2'),
           ('HR 8799 b', 'Paranal/SPHERE.IRDIS_D_H23_2', 'Paranal/SPHERE.IRDIS_D_H23_3', 'Paranal/SPHERE.IRDIS_D_H23_2'),
           ('HR 8799 d', 'Paranal/SPHERE.IRDIS_D_H23_2', 'Paranal/SPHERE.IRDIS_D_H23_3', 'Paranal/SPHERE.IRDIS_D_H23_2'),
           ('PZ Tel B', 'Paranal/SPHERE.IRDIS_D_H23_2', 'Paranal/SPHERE.IRDIS_D_H23_3', 'Paranal/SPHERE.IRDIS_D_H23_2')]

Plotting a color-magnitude diagram

Finally, we will create the color-magnitude diagram by providing the ColorMagBox objects as list to the boxes parameter of plot_color_magnitude. Further details on the various parameters of this function can be found in the API documentation.

[21]:
species.plot_color_magnitude(boxes=[colorbox, modelcolor1, modelcolor2],
                             objects=objects,
                             mass_labels=[(3., 'left'), (5., 'left'), (10., 'left'), (20., 'right'), (50., 'right'), (100., 'right')],
                             companion_labels=True,
                             field_range=('late M', 'late T'),
                             label_x=r'H2 $-$ H3',
                             label_y=r'M$_\mathregular{H2}$',
                             xlim=(-2.7, 1.8),
                             ylim=(17.3, 7),
                             offset=(-0.07, -0.1),
                             legend=(0.08, 0.8),
                             output='color_mag.png')
Plotting color-magnitude diagram: color_mag.png... [DONE]
[21]:
from IPython.display import Image
Image('color_mag.png')
[21]:
../_images/tutorials_color_magnitude_narrowband_46_0.png