spino.pipeline.phase_tsm_calculator module

class spino.pipeline.phase_tsm_calculator.TSM_ESM_Calculator[source]

Bases: object

Transmission / Emission Spectroscopy Metrics (Kempton et al. 2018) PASP, 130, 114401 - doi:10.1088/1538-3873/aadf6f

TSM = ScaleFactor × (Rp³ Teq) / (Mp Rs²) × 10^(−mJ/5) [Eq. 1] ESM = 4.29×10⁶ × B₇.₅(Tday)/B₇.₅(T★) × (Rp/R★)² × 10^(−mK/5) [Eq. 4] Tday = 1.10 × Teq

SCALE = {'large_sub_neptune': 1.28, 'small_sub_neptune': 1.26, 'sub_jovian': 1.15, 'terrestrial': 0.19}
TSM_THR = {'large_sub_neptune': 90, 'small_sub_neptune': 90, 'sub_jovian': 90, 'terrestrial': 10}
ESM_THR = 7.5
R_EARTH_CM = 637100000.0
R_SUN_CM = 69570000000.0
RHO_EARTH = 5.515
RJUP2REARTH = 11.2
MJUP2MEARTH = 317.8
compute(p)[source]
Parameters:

  • p (dict, single self-contained planet dictionary.)

  • keys (Required)

  • -------------

  • name (str (planet identifier))

  • radius) (Rp_Rearth OR Rp_Rjup (planet)

  • Rs_Rsun (float (stellar radius, [R☉]))

  • Teff_star (float (stellar Teff, [K]))

  • mag_J (float (host star J mag))

  • mag_K (float (host star K mag))

  • of (Plus ONE) – Teq : float (equilibrium temperature [K]) a_AU : float (semi-major axis [AU]; Teq computed)

  • Optional

  • --------

  • mass) (Mp_Mearth OR Mp_Mjup (planet) – (if absent → Chen & Kipping 2017 estimate)

  • P_days (float (orbital period [d]))

  • a_AU (float (semi-major axis [AU]))

Return type:

dict with TSM, ESM and all intermediates.

spino.pipeline.phase_tsm_calculator.analyze(planets, output_csv='tsm_esm_results.csv')[source]