Welcome to the GSFC Exoplanet Modeling and Analysis Center (EMAC)
EMAC serves as a catalog, repository and integration platform for modeling and analysis resources focused on the study of exoplanet characteristics and environments. EMAC is a key project of the GSFC Sellers Exoplanet Environments Collaboration (SEEC).
If you've used EMAC in any part of your research, please cite our RNAAS paper either in your methods section or in the "Software used" portion of any manuscripts; see the FAQ for more information.
More Information on EMAC for first-time visitors...If you make use of tools linked or hosted on EMAC: please use the following statement in your publication acknowledgements: “This research made use of the NASA Exoplanet Modeling and Analysis Center (EMAC), which is funded by the NASA Planetary Science Division’s Internal Scientist Funding Model.”
Stay up to date with EMAC!- Subscribe to our monthly RSS messages on new updates and tools
- Check out the Twitter account @ExoplanetModels (not an official NASA account), where new tools and features are highlighted
- Email us with general feedback at and tell us what you’d change or improve.
- Click the
icon in a resource box to provide suggestions for an individual tool or tools.
- EMAC is intended as a clearinghouse for the whole research community interested in exoplanets, where any software or model developer can submit their tool/model or their model output as a contribution for others to use.
- EMAC provides a searchable and sortable database for available source code and data output files - both resources hosted locally by EMAC as well as existing external tools and repositories hosted elsewhere.
- The EMAC team also helps develop new web interfaces for tools that can be run “on-demand” or model grids that can be interpolated for more individualized results.
- If you would like to submit a new tool/model to EMAC, please visit our Submit a Resource page.
- For help with tutorials for select resources/tools use the “Demo” buttons below and subscribe to our YouTube channel.
- Watch this video for a walk-through of the whole EMAC site, including how to submit a new tool and how to access information for each resource.
The P.I. is Avi Mandell, and the Deputy P.I. is Eric Lopez; more information on EMAC staffing and organization can be found on Our Team page.

SWAMPE: A Shallow-Water Atmospheric Model in Python for Exoplanets
Landgren, E. and Nadeau, A.
SWAMPE is a Python package for modeling the dynamics of exoplanetary atmospheres. SWAMPE is an intermediate-complexity, two-dimensional shallow-water general circulation model. Benchmarked for synchronously rotating hot Jupiters and sub-Neptunes, the code is modular and could be easily modified to model dissimilar space objects, from Brown Dwarfs to terrestrial, potentially habitable exoplanets. SWAMPE can be easily run on a personal laptop.

Bioverse: Simulation framework for Bayesian hypothesis testing of statistical exoplanet missions
Alex Bixel et al.
Bioverse is a quantitative framework for assessing the diagnostic power of a statistical exoplanet survey. It combines Gaia-based stellar samples with Kepler-derived exoplanet demographics and a mission simulator that enables exploration of a variety of observing, follow-up, and characterization strategies. Uniquely, Bioverse contains a versatile module for population-level hypothesis testing supporting trade studies and survey optimization. It currently supports direct imaging or transit missions, but its modularity makes it adaptable to any mission concept that makes measurements on a sample of exoplanets.

TurbospectrumNLTE: Synthetic stellar spectra calculator LTE / NLTE
Plez, Bertrand; Gerber, Jeff; Magg, Ekaterina; Bergemann, Maria
Latest version of TS (Turbospectrum), with NLTE capabilities. Computation of stellar spectra (flux and intensities) in 1D or average <3D> stellar atmosphere models. In order to compute NLTE stellar spectra, additional data is needed, downloadable outside GitHub. See documentation in DOC folder Python wrappers are available at https://github.com/EkaterinaSe/TurboSpectrum-Wrapper/ and https://github.com/JGerbs13/TSFitPy They allow interpolation between models and fitting of spectra to derive stellar parameters.

Applefy: Robust detection limits for high-contrast imaging
Bonse, Markus J.; Gebhard, Timothy D.
Applefy calculates detection limits for exoplanet high contrast imaging (HCI) datasets. The package provides a number of features and functionalities to improve the accuracy and robustness of contrast curve calculations. Applefy implements the classical approach based on the t-test as well as the parametric boostrap test for non-Gaussian residual noise. Written in Python, it computes contrast curves and contrast grids.

Tiberius: Time series spectral extraction and transit light curve fitting
Kirk, J.
Tiberius is a Python library for reducing time series spectra and fitting exoplanet transit light curves. This can be used to extract spectra from JWST (all 4 instruments), along with ground-based long-slit spectrographs and Keck/NIRSPEC echelle spectra (beta). The light curve fitting routines can be used as as standalone to fit, for example, HST light curves extracted with other methods.

light-curve: Irregular time series analysis toolbox for Rust and Python
Malanchev, K.; Lavrukhina, A.
The "light-curve" project is a toolbox for analyzing irregular time-series data, consisting of two components: a feature extractor and a dm-dt mapper. Feature extraction is available as Rust and Python libraries, including various feature extractors like magnitude statistics, shape-based features, Lomb-Scargle periodogram peaks, and parametric fits. The dm-dt mapper represents observation pairs as 2-D points based on magnitude and time differences, available in the same Python library, the Rust library, and a binary executable for generating PNG images.

ExoCcycleGeo: Geophysical and geochemical controls on abiotic carbon cycling on Earth-like planets
Neveu, M. et al.
This model of geological carbon fluxes leverages widely used geochemical codes of solid-melt equilibria for silicate rocks (MELTS) and of equilibria and kinetics of water-rock interactions (PHREEQC). Coupled with a simple numerical computation of global thermal evolution, this model enables investigation of the effects of planet size (mass) and bulk, surface, and upper mantle composition on carbon cycling through geologic time. Its applicable size range (0.5 to 2 Earth masses) is limited by the fidelity of the geodynamic model. The applicable range of compositions is limited by those that can be handled by MELTS and PHREEQC.

MATRIX ToolKit: ToolKit for Multi-phAse Transits Recovery from Injected eXoplanets
Dévora-Pajares, Martín & Pozuelos, Francisco J.
The MATRIX ToolKit has been specially designed to establish detection limits of photometry data sets by performing robust injection-and-recovery analyses on a three dimensional grid of scenarios (orbital period vs planetary radius vs transit epoch). This kind of scientific detection threshold determination can now be done with a simple python command with the significant addition of taking into account different transit epochs, which helps to establish a more reliable detection rate for a given period and radius.

WATSON: Visual Vetting and Analysis of Transits from Space ObservatioNs
Dévora-Pajares, M.
WATSON (Visual Vetting and Analysis of Transits from Space ObservatioNs is a lightweight software package that enables a comfortable visual vetting of a transiting signal candidate from Kepler, K2 and TESS missions. WATSON looks for transit-like signals that could be generated by other sources or instrument artifacts. The code runs simplified tests on scenarios including:
- Transit shape model fit
- Odd-even transits checks
- Centroids shifts
- Optical ghost effects
- Transit source offsets
- and more...

PEPITA: Prediction of Exoplanet Transit Parameters Precisions using Information Analysis Techniques
Julio Hernandez Camero
PEPITA is a Python package that allows making predictions for the precision of exoplanet parameters using transit light-curves, without the need of performing a fit to the data. Behind scenes, it makes use of the Information Analysis techniques to predict the best precision that can be obtained by fitting a light-curve without actually needing to perform the fit.

PyMieScatt: The Python Mie Scattering package
Sumlin, B., Heinson, W., and Chakrabarty, R.
PyMieScatt is a comprehensive forward and inverse Mie theory solver for Python 3. This package calculates relevant parameters such as absorption, scattering, extinction, asymmetry, backscatter, and more. It also contains single-line functions to calculate optical coefficients (in Mm-1) of ensembles of particles in lognormal (with single or multiple modes) or custom size distributions. The inverse calculations retrieve the complex refractive index from laboratory measurements of scattering and absorption (or backscatter), useful for studying atmospheric organic aerosol of unknown composition. Read more in our JQSRT paper!

SHERLOCK: Searching for Hints of Exoplanets fRom Lightcurves Of spaCe-based seeKers
Dévora-Pajares, M. & Pozuelos, F. J. et al.
SHERLOCK is an end-to-end pipeline that allows the users to explore the data from space-based missions to search for planetary candidates. It can be used to recover alerted candidates by the automatic pipelines such as SPOC and the QLP, the so-called Kepler objects of interest (KOIs) and TESS objects of interest (TOIs), and to search for candidates that remain unnoticed due to detection thresholds, lack of data exploration or poor photometric quality.

PACMAN: A pipeline to reduce and analyze Hubble Wide Field Camera 3 IR Grism data
Zieba, Sebastian and Kreidberg, Laura
Here we present PACMAN, an end-to-end pipeline developed to reduce and analyze HST/WFC3 data. The pipeline includes both spectral extraction and light curve fitting. The foundation of PACMAN has been already used in numerous publications (e.g., Kreidberg et al., 2014; Kreidberg et al., 2018) and these papers have already accumulated hundreds of citations.

PHOTOe: An efficient Monte Carlo model for the slowing down of photoelectrons
Monte Carlo model for simulating the slowing down of photoelectrons in gases in the local deposition approximation.
Version v1 works with H and He atoms plus thermal electrons. The model is described in: García Muñoz, Icarus, Volume 392, 1 March 2023, 115373
The model is available on https://antoniogarciamunoz.wordpress.com/ and upon email request from the author.

SpecMatch-Empirical: Spectroscopic characterization of stars with an empirical spectral library
Yee, Samuel; Petigura, Erik; von Braun, Kaspar
We present SpecMatch-Empirical, a tool for measuring the fundamental properties of stars from their spectra by comparing them against an empirical spectral library of FGKM stars. The spectral library comprises high-resolution, high signal-to-noise observed spectra from Keck/HIRES for 404 touchstone stars with well-determined stellar parameters derived from interferometry, asteroseismology, and spectrophotometry. The code achieves accuracies of 100K, 15%, and 0.09 dex in Teff, Rstar, and [Fe/H] respectively for FGKM dwarfs.

Astroquery: Library with tools to query astronomical databases
Astroquery is a collection of tools for requesting data from databases hosted on remote servers with interfaces exposed on the internet, including those with web pages but without formal application program interfaces. These tools are built on the Python requests package, which is used to make HTTP requests, and astropy, which provides most of the data parsing functionality. Astroquery modules generally attempt to replicate the web page interface provided by a given service as closely as possible, making the transition from browser-based to command-line interaction easy. Astroquery enables the creation of fully reproducible workflows from data acquisition through publication.

spaceKLIP: Pipeline for reducing & analyzing JWST coronagraphy data with the help of pyKLIP
JWST ERS 1386 Team
The spaceKLIP pipeline enables to reduce & analyze JWST NIRCam & MIRI coronagraphy data. It provides functions to run the official jwst stage 1 and 2 data reduction pipelines with several modifications that were made to improve the quality of high-contrast imaging reductions. It then performs PSF subtraction based on the KLIP algorithm as implemented in the widely used pyKLIP package, outputs contrast curves, and enables forward model PSF fitting for any detected companions in order to extract their properties (offset and flux). The pipeline is still under heavy development.

Butterpy: realistic star spot evolution and light curves in Python
Claytor, Zachary R. et al.
Butterpy is a Python-based tool for simulating star spot emergence, evolution, and decay as well as stellar rotational light curves. It is adapted from the physically motivated model used by Aigrain et al. (2015, MNRAS, 450, 3211) to test the recovery of stellar rotation periods using different frequency analysis techniques. Butterpy allows the user to simulate light curves of stars with variable activity level, rotation period, spot lifetime, magnetic cycle duration and overlap, spot emergence latitudes, and latitudinal differential rotation shear. The name Butterpy is a portmanteau of "butterfly" (like the solar butterfly diagram) and "Python."

Kiauhoku: Python utilities for stellar model grid interpolation
Claytor, Zachary R. et al.
Kiauhoku is a Python package for interacting with, interpolating, and fitting stellar evolutionary tracks to observational data. It includes popular stellar model grids like MIST, Dartmouth, and GARSTEC, as well as a few custom YREC grids, with more being added over time.
From Hawaiian:
1. vt. To sense the span of a star's existence (i.e., its age).
2. n. The speed of a star (in this case, its rotational speed).
This name was created in partnership with Dr. Larry Kimura and Bruce Torres Fischer, a student participant in A Hua He Inoa, a program to bring Hawaiian naming practices to new astronomical discoveries. We are grateful for their collaboration.

rfast: A fast tool for planetary spectral forward and inverse modeling.
T. Robinson
The rfast tool is an ultra-quick planetary spectrum simulator and remote sensing tool, originally designed for rapid retrieval explorations for mission concept studies. Through a convenient runscript, users can generate a noise-free spectrum of a planetary environment, add instrumental noise, and perform inverse modeling. The rfast tool is capable of applications to simulated and real observations spanning reflected-light, thermal emission, and transit transmission.

Molecfit: A general tool for telluric absorption correction
Molecfit is a tool to correct for telluric absorption lines based on synthetic modelling of the Earth’s atmospheric transmission. It can be used with data obtained with various ground-based telescopes and instruments. It combines a publicly available radiative transfer code, a molecular line database, atmospheric profiles, and various kernels to model the instrument LSF. The atmospheric profiles are created by merging a standard atmospheric profile representative of a given observatory’s climate, of local meteorological data, and of dynamically retrieved altitude profiles for temperature, pressure, and humidity.

Pyshellspec: Binary systems with circumstellar matter (β Lyrae)
Brož, M., Nemravová, J.
PYSHELLSPEC is an astrophysical tool for modeling of binary systems with circumstellar matter (e.g. accretion disk, jet, shell), computation of interferometric observables |V2|, arg T3, |T3|, |dV|, arg dV, comparison of light curves, spectro-interferometry, spectra, and SED with observations, and both global and local optimisation of system parameters. It is based on Shellspec, a long-characteristic LTE radiation transfer code by Budaj & Richards (2004).

AnalyticLC: An Accurate 3D Analytic Modeling Tool for Exoplanetary Photometry, Radial Velocity and Astrometry
Yair Judkovsky, Aviv Ofir and Oded Aharonson
AnalyticLC is an accurate photometry, radial velocity, and astrometry modeling tool. It is based on a fourth-order expansion of the Disturbing Function, incorporating 3D orbital dynamics. The analytic approach of AnalyticLC elucidates the relation between orbital dynamics and observable quantities. In addition, it offers advantages for analyzing observations with a long time span, a scenario becoming increasingly common in this era of multiple space missions. AnalyticLC has been used to interpret Kepler data and obtain estimates of more than a hundred exoplanets' physical and orbital properties.

Optimal BLS: Optimally-efficient code for transit searches in long time series
Aviv Ofir
Explicitly including Keplerian dynamics in the transit search allows Optimal BLS to enhance transit detectability while allowing such searches to be done with much-reduced resources and time. By using the (standard) BLS, one is either fairly insensitive to long-period planets or less sensitive to short-period planets and computationally slower by a significant factor of ~330 (for a 3 yr long dataset). Physical system parameters, such as the host star's size and mass, directly affect transit search. This understanding can then be used to optimize the search for every star individually. The code is well-used by the community.

pyPplusS: Fast and precise light-curve model for transiting exoplanets with rings
Eden Rein and Aviv Ofir
The Polygon + Segments model allows modeling the light curve of an exoplanet with rings. This high-precision model includes full ring geometry as well as possible ring transparency and the host star’s limb darkening. Additionally, it can model oblate ringless planets as an opaque “ring” (same shape as a planet). pyPplusS is also computationally efficient, requiring just a 1D integration over a small range, making it faster than existing techniques. The algorithm at its core is further generalized to compute the light curve of any set of convex primitive shapes in transit (e.g. multiple planets, oblate planets, moons, rings, combination thereof, etc.) while accounting for their overlaps.

PySME: Stellar Spectral Synthesis and Parameter Fitting
Wehrhahn, A., Piskunov, N., Valenti, J.
Spectroscopy Made Easy (SME) is a software tool that fits an observed spectrum of a star with a model spectrum. Since its initial release in 1996, SME has been a suite of IDL routines that call a dynamically linked library, which is compiled from C++ and Fortran. This classic IDL version of SME is available for download. In 2018, we began began reimplementing the IDL part of SME in python 3, adopting an object oriented paradigm and continuous integration practices (code repository, build automation, self-testing, frequent builds).

ATOCA: Algorithm to Treat Order Contamination. Used to decontaminate and extract spectra from image.
ATOCA is used to extract and decontaminate spectroscopic images with multiple sources or diffraction orders. The inputs are, for all orders and sources: the wavelength solutions, the trace profiles, the throughputs and the spectral resolution kernels. From this, ATOCA can model simultaneously the detector and extract the spectra. See Darveau-Bernier et al. (2022) for more details.

FORECAST: Finely Optimised REtrieval of Companions of Accelerating STars
M. Bonavita, R. Gratton, C. Fontanive, A. Sozzetti
A significant proper motion difference between two catalogues for a given star is a good indication of the presence of a perturbing body. FORECAST allows you to identify the region where a companion compatible with the measured Δμ should appear if the star is directly imaged. It also provides an estimate of the mass of the companion compatible with the astrometric signal at each position in the allowed region. FORECAST maps can be used both to identify and confirm potential direct imaged sub-stellar candidates.

Exo-DMC: Exoplanet Detection Map Calculator
M. Bonavita, Silvano Desidera, Ernst de Moij, Arthur Vigan, Justine Lannier
Exo-DMC (Exoplanet Detection Map Calculator) is a Monte Carlo tool for the statistical analysis of exoplanet surveys results. It combines the information on the target stars with the instrument detection limits to estimate the probability of detection of a given synthetic planet population, ultimately generating detection probability maps. The Exo-DMC is the latest (although the first one in Python) rendition of the MESS (Multi-purpose Exoplanet Simulation System). Like MESS, the DMC allows for a high level of flexibility in terms of possible assumptions on the synthetic planet population to be used for the determination of the detection probability.

blase: An Interpretable Machine Learning Framework for Modeling High-Resolution Spectroscopic Data
Michael Gully-Santiago & Caroline V. Morley
Blasé introduces a powerful new approach to whole-spectrum fitting: clone 10,000+ spectral lines from a precomputed synthetic spectral model template, and then learn the perturbations to those lines through comparison to real data. Each spectral line has 4 parameters, yielding possibly 40,000+ parameters. The technique hinges on the magic of autodiff, the enabling technology behind Machine Learning, to tune all of those parameters precisely and quickly. The tool has conceivable extensions to Doppler imaging, Precision RV's, abundances, and more. It is built in PyTorch, with native GPU support.

SysSimPyPlots: Loading, analyzing, and plotting catalogs generated from the SysSim models
SysSimPyPlots is a Python package for loading, plotting, and otherwise visualizing the simulated catalogs generated from ExoplanetsSysSim, a comprehensive forward modeling framework for studying planetary systems based on the Kepler mission. In particular, it is designed to work with the SysSim clustered planetary system models (https://github.com/ExoJulia/SysSimExClusters) that characterize the underlying occurrence and intra-system correlations of multi-planet systems.

SysSimPyMMEN: Inferring Minimum-Mass Extrasolar Nebulae from the SysSim models
SysSimPyMMEN is a Python package for inferring the minimum-mass extrasolar nebula (MMEN), a power-law profile for the minimum mass in disk solids required to form the existing exoplanets if they formed in their present locations. It is designed to work with the SysSim clustered planetary system models (https://github.com/ExoJulia/SysSimExClusters) that characterize the underlying occurrence and intra-system correlations of multi-planet systems, but can be easily applied to any other planetary system by the user.

RAPOC: Rosseland And Planck Opacity Converter
Lorenzo V. Mugnai and Darius Modirrousta-Galian
RAPOC (Rosseland and Planck Opacity Converter) is a Python 3 code that calculates Rosseland and Planck mean opacities from wavelength-dependent opacities for a given temperature, pressure, and wavelength range. In addition to being user-friendly and rapid, RAPOC can interpolate between discrete data points, making it flexible and widely applicable to the astrophysical and Earth-sciences fields, as well as in engineering. For the input data, RAPOC can use ExoMol and DACE data, or any user-defined data, provided that it is in a readable format.

Kamodo: a CCMC tool for access, interpolation, and visualization of data in python.
The Community Coordinated Modeling Center at NASA GSFC
Kamodo allows model developers to represent simulation results as mathematical functions which may be manipulated directly by end users. Kamodo handles unit conversion transparently and supports interactive science discovery through jupyter notebooks with minimal coding and is accessible through python.

ECLIPS3D: Public code for linear wave and circulation calculations
F. Debras et al.
Public Fortran 90 code for linear wave and circulation calculations, developed originally for planetary atmospheres, with python scripts provided for data analysis.4 setups are provided: 2D_axi: eigenvector setup in spherical coordinates assuming axisymmetry around the axis of rotation. A longitudinal wavenumber, m, must therefore be provided. 2D_shallow: eigenvector setup for shallow water beta-plane. The latitude of the beta plane and characteristic height can be changed. 3D: eigenvector setup in full 3D, spherical coordinates. 3D_steady: linear circulation setup, hence matrix inversion. A forcing and a dissipation have to be implemented for a linear steady state to exist.

TESS-SIP: TESS Systematics Insensitive Periodogram
Hedges et al.
Tool for creating a Systematics-insensitive Periodogram (SIP) to detect long period rotation in NASA's TESS mission data. Read more in our published Research Note of the American Astronomical Society. SIP is a method of detrending telescope systematics (the TESS scattered light) simultaneously with calculating a Lomb-Scargle periodogram. This allows us to estimate of the rotation rate of variables with a period of >30days when there are multiple sectors. You can read a more in-depth work of how SIP is used in NASA's Kepler/K2 data here.

IGRINS RV: A Radial Velocity Pipeline for IGRINS
Stahl, A. et al. 2021; Tang, S-Y. et al. 2021
IGRINS RV is a python open source pipeline for extracting radial velocities (RVs) from spectra taken with the Immersion GRating INfrared Spectrometer (IGRINS). It uses a modified forward modeling technique that leverages telluric absorption lines as a common-path wavelength calibrator. IGRINS RV achieves an RV precision in the H and K bands of around 25-30 m/s for narrow-line stars, and it has successfully recovered the planet-induced RV signals of both HD 189733 and τ Boo A. Visit Stahl et al. 2021 to see the published paper.

AccretR: A planetary accretion and composition code in R
Mohit Melwani Daswani
AccretR is a planetary body accretion Monte Carlo code that calculates mass, radius and bulk composition along a specified growth track, for orderly/hierarchical, runaway, and random particle accretion models, optimized for icy ocean worlds in our Solar System (priors can be modified for other systems). Elements in the model include concentrations of: H, C, N, O, Na, Mg, Al, Si, S, Cl, K, Ca, and Fe. Maximal water is also computed, assuming all H goes into forming water. Accretional heat is also calculated.

FastChem: Ultra-fast Equilibrium Chemistry
Daniel Kitzmann, Joachim Stock
FastChem 2 is a new version of the established semi-analytical thermochemical equilibrium code FastChem. Whereas the original version of FastChem is limited to atmospheres containing a significant amount of the element hydrogen, FastChem 2 is now also applicable to chemical mixtures dominated by any other species such as CO2, N2, or Si for example. The code is written in object-oriented C++ and also offers an optional Python module.

Staralt: Calculating Object Visibility for Ground-based Telescopes
Isaac Newton Group of Telescopes
Staralt is a program that shows the observability of objects in various ways: either you can plot altitude against time for a particular night (Staralt), or plot the path of your objects across the sky for a particular night (Startrack), or plot how altitude changes over a year (Starobs), or get a table with the best observing date for each object (Starmult).

ExoAtmospheres: IAC community database for exoplanet atmospheric observations
The Exoplanets and Astrobiology group at IAC
The Exoplanet Atmospheres Database is built for the community and maintained by the community. Exoplanet atmospheres is an exciting and vibrant field of research, where new discoveries and publications occur at a very fast pace, and it is easy to miss many interesting results. The main purpose of this database is to become a quick and useful repository of all available exoplanet atmospheres observations, and also to help in the gathering of useful references for a given planet or planet types.

THAI: TRAPPIST Habitable Atmosphere Intercomparison GCM Data Repository
THAI Team (T. Fauchez et al.)
The TRAPPIST Habitable Atmosphere Intercomparison (THAI) project is a model inter-comparison effort between four GCMs: ExoCAM, LMD-G, ROCKE3D and the UM – examining a single interesting test case (TRAPPIST-1e) under several different atmosphere scenarios. The CKAN data repository provides NetCDF files for each case, allowing for examination and intercomparison of results from the different models. Scripts to process the data and plot them are available on our Github repository.

The TESS Triple-9 (TT9) Catalog: 999 uniformly vetted TESS candidate exoplanets
Luca Cacciapuoti et al.
The TESS Triple 9 Catalog (TT9, in short) contains dispositions for 999 candidate exoplanetary signals detected by the Transiting Exoplanet Survey Satellite (TESS). These signals were classified as Planet Candidate (PC), False Positive (FP) or Potential False Positive (PFP) based on TESS data (both light curves and images) as well as ancillary information such as stellar catalog. The data has been analyzed using the DAVE pipeline (https://github.com/exoplanetvetting/DAVE) and with the aid of Citizen Science, through the NASA funded Planet Patrol project (https://www.zooniverse.org/projects/marckuchner/planet-patrol). DAVE pdf products are available on https://exofop.ipac.caltech.edu/tess

Prose: A Python framework for modular astronomical images processing
Lionel Garcia
Built for astronomy, Prose is instrument-agnostic and allows the construction of data reduction pipelines using a wide range of building blocks, pre-implemented or user-defined. Using its modular architecture, it features basic reduction pipelines to deal with common tasks such as automatic reduction and photometric extraction.

exoVista: Planetary System Models for Survey Analyses
Christopher Stark
ExoVista quickly produces a large number of planetary system models that can serve as a standard set for simulated study with the direct imaging, transit, astrometric, interferometric, and astrometric methods at scattered light wavelengths. ExoVista distributes planets consistent with the Kepler occurrence rates around Hipparcos stars within 50 pc, assigns a mass, radius, and albedo to each planet, checks for stability of the orbits, evolves all objects with a gravitational n-body integrator, and generates a quasi-self-consistent debris disk for each system consistent with the LBTI HOSTS exozodi survey.

M_-M_K-: Estimating realistic stellar masses from magnitudes
Andrew Mann et al.
Code converts absolute 2MASS Ks-band magnitude (or a distance and a Ks-band magnitude) into an estimate of the stellar mass using the empirical relation derived from the resolved photometry and orbits of astrometric binaries. The code outputs errors based on the relationship's scatter and errors in the provided distance and Ks magnitude.

DustPy: A Python Package for Dust Evolution in Protoplanetary Disks
Sebastian Markus Stammler ; Tilman Birnstiel

ALMA: A Fortran program for computing the viscoelastic Love numbers of a spherically symmetric planet
Melini, D., Saliby, C., Spada, G.
ALMA is a Fortran code that computes loading and tidal Love numbers for a spherically symmetric, radially stratified, incompressible planet. ALMA can evaluate i) real (time-domain) Love numbers and their time derivatives for a Heaviside or ramp-shaped forcing time history, or ii) complex (frequency-domain) Love numbers for a periodic forcing. The planetary structure can include an arbitrary number of homogeneous layers, and each layer can have a different rheological law. ALMA can model the following linear rheologies: Elastic, Maxwell visco-elastic, Newtonian viscous fluid, Kelvin-Voigt solid, Burgers and Andrade transient rheologies. Additional rheological laws can be easily implemented.

VPLanet: The Virtual Planet Simulator
Rory Barnes et al.
VPLanet simulates planetary system evolution with a single executable: 1) thermal and magnetic evolution of terrestrial planets, 2) magma oceans, 3) radiogenic heating of interiors, 4) tidal effects, 5) rotational axis evolution, 6) stellar evolution, including pre-MS, XUV, and spin-down, 7) stellar flares, 8) climate via a 1-D EBM, 9) atmospheric escape, including water photolysis and H escape, 10) approximate orbital evolution, 11) exact orbital evolution, 12) circumbinary planet orbits, and 13) galactic perturbations on planetary systems. The code is validated by reproducing selected Solar System, exoplanet, and binary star properties. Documentation and numerous examples are provided.

Eureka!: An End-to-End Pipeline for JWST Time-Series Observations
Bell, T. J. et al.
Eureka! is a data reduction and analysis pipeline for exoplanet time-series observations, with a particular focus on James Webb Space Telescope (JWST) data. The goal of Eureka! is to provide an end-to-end pipeline that starts with raw, uncalibrated FITS files and ultimately yields precise exoplanet transmission and/or emission spectra. The pipeline has a modular structure with six stages, each with intermediate figures and outputs that allow users to compare Eureka!’s performance using different parameter settings or to compare Eureka! with an independent pipeline.

Pytmosph3R: Transmission/emission spectra from 3D atmospheric simulations (GCMs, ...)
Falco, Aurélien ; Pluriel, William ; Leconte, Jérémy ; Caldas, Anthony
Pytmosph3R is a Python-3 library that computes transmission and emission spectra based on 3D atmospheric simulations, for example performed with the LMDZ generic global climate model. Pytmosph3R can be used in notebooks or on the command line, using a configuration similar to that of TauREx. The library should include a feature to generate phase/light-curves in the next release.

celmech: A Python package for celestial mechanics
Hadden, Sam; Tamayo, Daniel
An open-source Python package designed to facilitate a wide variety of celestial mechanics calculations. The package allows users to formulate and integrate equations of motion by incorporating user-specified terms from the classical disturbing function expansion. The package is designed to interface seamlessly with the popular REBOUND N-body code.

pycheops: Python package for the analysis of light curves from the ESA CHEOPS mission
Maxted et al.
The pycheops python module is an open-source software package that has been developed to easily and efficiently analyse CHEOPS light curve data using state-of-the-art techniques. The models in the package can also be applied to other types of data. The package included a "cook book" and examples, plus a command-line tool that aids in the preparation of observing requests for CHEOPS observers (make_xml_files). For discussion and announcements, please join the pycheops google group.

HELIOS-K: A GPU opacity calculator for exoplanetary atmospheres
Simon Grimm, Kevin Heng
HELIOS-K calculates opacity functions for planetary atmopheres by using opacity line lists from different databases. Before the opacity functions can be calculated, the line lists need to be downloaded and preprocessed into binary files that can be read from HELIOS-K. HELIOS-K provides tools to automatically download and preprocess files from the ExoMol, HITRAN, HITEMP, NIST, Kurucz and VALD3 databases. HELIOS-K is running on GPUs and require a Nvidia GPU with compute capability of 3.0 or higher.

l1periodogram: Searching for planetary signatures in radial velocity time-series with a sparse recovery technique
Nathan C. Hara
The l1 periodogram is designed to search for periodicities in unevenly sampled time series. It can be used similarly as a Lomb-Scargle periodogram, and retrieves a figure which has a similar aspect but has fewer peaks due to aliasing. It is primarily designed for the search of exoplanets in radial velocity data, but can be also used for other purposes. It is based the sparse recovery technique called "Basis Pursuit" (Chen & Donoho 1998).

Magnitude-squared coherence: Frequency-domain view of the cross-correlation between RV and activity indicator time series
Dodson-Robinson, S. E., et al.
NWelch is a python3 implementation of Welch's method for estimating the magnitude-squared coherence (MSC) between contemporaneous RV and activity-indicator time series. While it's impossible to directly calculate the cross-correlation between two unevenly sampled time series, we can use a nonuniform fast Fourier transform to estimate the frequency-domain version of their cross-correlation - the MSC. Observers should be suspicious of planet candidates at frequencies with high MSC. For univariate time series (for example, RV only), Welch's method can deliver power spectrum estimates with lower variance than the Lomb-Scargle periodogram.

MARGE: A Python package to train and evaluate neural networks
Himes et al.
MARGE (Machine learning Algorithm for Radiative transfer of Generated Exoplanets) is an all-in-one package to generate exoplanet spectra across a defined parameter space, process the output, and train machine learning (ML) models as a fast approximation to radiative transfer (RT). Despite its backronym name, MARGE is a general package that can train neural networks on a provided data set of inputs and outputs. MARGE is an open-source project under the Reproducible Research Software License and welcomes improvements from the community to be submitted via pull requests on Github.

HOMER: A Bayesian inverse modeling code
Himes et al.
HOMER (the Helper Of My Eternal Retrievals) is a Bayesian inverse modeling code. Given some data and uncertainties, the posterior distribution is determined for some model. HOMER allows for both nested sampling and Markov chain Monte Carlo (MCMC) frameworks. HOMER's forward model is a neural network (NN) surrogate model trained by MARGE. For details on MARGE, see the MARGE User Manual at https://exosports.github.io/MARGE/doc/MARGE_User_Manual.html. HOMER is released under the Reproducible Research Software License and welcomes community contributions via pull requests on Github.

Aladin: Aladin Sky Atlas
CDS - Centre de Données astronomiques de Strasourg, Université de Strasbourg/CNRS
Aladin is an interactive sky atlas allowing the user to visualize digitized astronomical images or full surveys, superimpose entries from astronomical catalogues or databases, and interactively access related data and information from the Simbad database, the VizieR service and other archives for all known astronomical objects in the field. The Aladin sky atlas is available in two modes: Aladin Desktop, a regular application and Aladin Lite an HTML5 javascript web widget.

FitsMap: A Simple, Lightweight Tool For Displaying Interactive Astronomical Image and Catalog Data
Ryan Hausen and Brant E. Robertson
The visual inspection of image and catalog data continues to be a valuable aspect of astronomical data analysis. As the scale of astronomical image and catalog data continues to grow, visualizing the data becomes increasingly difficult. In this work, we introduce FitsMap, a simple, lightweight tool for visualizing astronomical image and catalog data. FitsMap only requires a simple web server and can scale to over gigapixel images with tens of millions of sources. Further, the web-based visualizations can be viewed performantly on mobile devices. FitsMap is implemented in Python and is open source (https://github.com/ryanhausen/fitsmap).

exocartographer: Forward-model for constraining exoplanet maps and orbital parameters from reflected lightcurves
Farr, B., Farr, W., Cowan, N., Haggard, H., and Robinson, T.

SteParSyn: A Bayesian code to infer stellar atmospheric parameters using spectral synthesis
Tabernero, H. M. et al.
SteParSyn is a Python code designed to infer the stellar atmospheric parameters Teff, log(g), and [Fe/H] of late-type stars (FGKM) using the spectral synthesis method. It uses a Markov chain Monte Carlo (MCMC) sampler to explore the parameter space by comparing synthetic spectra to the observations. The code has been employed to study stars in open clusters, cepheids, stars in the Magellanic clouds, exoplanet hosts observed with ESPRESSO and CARMENES, and to characterize the first super-AGB candidate in our Galaxy. The code is available to the community in a GitHub repository.

PyExoRaMa: An Interactive Tool in Python to Investigate the Radius-Mass Diagram for Exoplanets
Amadori, F.; Damasso, M.; Zeng, L.; Sozzetti, A.
This is the python version of the software originally developed with Mathematica by Zeng et al (https://doi.org/10.5281/zenodo.5899463). The code represents a useful tool for visualizing and manipulating data related to extrasolar planets and their host stars in a multi-dimensional parameter space. It can be used to identify possible interdependence among several physical parameters and to compare observables with theoretical models describing the exoplanet composition and structure. Our transposition to Python presents some new features with respect to the original version

special: SPEctral Characterization of ImAged Low-mass companions
Valentin Christiaens
special is a package developed for the spectral characterization of directly imaged low-mass companions (MLT dwarfs). Nonetheless, this toolkit can also be used in a more general way for the characterisation of any object with a measured spectrum, provided an input model or template grid. The available tools range from the calculation of spectral covariance matrices (e.g. for IFS datacubes) and empirical spectral indices to the Bayesian inference of atmospheric parameters provided an input grid of models. In the latter case, an MCMC or nested sampler can be used, and additional parameters such as (extra) black body component(s) and extinction can be considered.

GENGA: A GPU N-body integrator for planet formation and planetary system evolution
Simon Grimm, Joachim Stadel
GENGA is a hybrid symplectic N-body integrator, designed to integrate planet and planetesimal dynamics in the late stage of planet formation and stability analysis of planetary systems. GENGA is based on the integration scheme of the Mercury code (Chambers 1999), which handles close encounters with very good energy conservation. The GENGA code supports three simulation modes: Integration of up to 60000 - 100000 massive bodies, integration with up to a million test particles, or parallel integration of a large number of individual planetary systems. GENGA is written in CUDA C and runs on all NVidia GPUs with compute capability of at least 2.0.

SOAP 2.0: RV stellar activity simulation including spot and faculae
X. Dumusque, I. Boisse, N.Santos
SOAP 2.0 allows to model the RV effect (and the effect on the FWHM, CONTRAST and BIS SPAN of the CCF) induced by magnetic regions (spot and faculae) on the surface of a star. The temperature of the magnetic features, their size, their location on the stellar surface, the resolution of the instrument used, the stellar properties (radius, effective temperature, inclination), are all parameters that can be adjusted. The code is efficient, as all the backend computation are performed in C, and is friendly to use as the interface is in python.

NbodyGradient: Differentiable symplectic N-body code for arbitrary orbital architectures
Agol, E., Hernandez, D. & Langford, Z.
NbodyGradient is a symplectic integrator for Newtonian gravity and arbitrary N-body hierarchies. It computes the derivatives of the output with respect to the input coordinates with high numerical precision. It was developed for transit-timing analyses, and is the first code to give derivatives of the transit times with respect to the initial conditions (either masses & cartesian coordinates/velocities or orbital elements). It is written in the Julia language. It is being extended for modeling RV, astrometry, and photodynamics.

Roman Coronagraph Exposure Time Calculator: Estimates integration times for the Roman Coronagraph instrument
© 2021. Government sponsorship acknowledged. The research was carried out in part at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration.
The Roman Coronagraph Exposure Time Calculator (Roman_Coronagraph_ETC for short) is the public version of the exposure time calculator of the Coronagraph Instrument aboard the Nancy Grace Roman Space Telescope funded by NASA. Roman_Coronagraph_ETC methods are based upon peer reviewed research articles and a collection of instrumental and modeling parameters of both the Coronagraph Instrument and the Nancy Grace Roman Space Telescope. The values in these files do not contain any ITAR or export control information. Roman_Coronagraph_ETC is licensed under Apache v2.

RASSINE: A tool for stellar spectrum continuum fitting
Michael Cretignier
RASSINE is a Python3 code dedicated to fitting stellar continuum using an alpha-shape algorithm. The code is divided in a sequence of five consecutive steps (SNAKE) : 1) Smoothing 2) Neighbourhood local maxima 3) Alpha-shape 4) Killing outliers 5) Envelop interpolation The code contains a fully automatic mode in case of 1d spectrum is given as input but the stellar spectrum can be also fitted by hand using the GUI interface.

SISTER: Starshade Imaging Simulation Toolkit for Exoplanet Reconnaissance
S.R. Hildebrandt and S.B. Shaklan are the principal authors of SISTER
The Starshade Imaging Simulations tool is a versatile tool designed to provide enough accuracy and variety when predicting how an exoplanet system would look like in an instrument that utilizes an Starshade to block the light from the host star. The tool allows for controlling a set of parameters of the whole instrument that have to do with: (1) the Starshade design, (2) the exoplanetary system, (3) the optical system (telescope) and (4) the detector (camera). There is a built-in plotting software added, but the simulations may be stored on disk and be plotted with any other software. SISTER has an online tutorial and has been published in JATIS. Visit sister.caltech.edu for details.

spectools_ir: A suite of tools designed for analysis of medium/high-resolution IR molecular spectra
Colette Salyk
Spectools_ir is a small suite of tools designed for analysis of medium/high-resolution IR molecular astronomical spectra. It consists of three main sub-modules (flux_calculator, slabspec, and slab_fitter) as well as a 'utils' sub-module, with a few additional functions. Spectools_ir was written with infrared medium/high-resolution molecular spectroscopy in mind. It often assumes spectra are in units of Jy and microns, and it uses information from the HITRAN molecular database. Some routines are more general, but users interested in other applications should proceed with caution.

ROCKE-3D: A fully coupled ocean atmosphere 3-D General Circulation Model
Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics (ROCKE-3D) is a three-dimensional General Circulation Model (GCM) developed at the NASA Goddard Institute for Space Studies for the modeling of atmospheres and oceans of solar system and exoplanetary terrestrial planets.

BASTA: The BAyesian STellar Algorithm
The BASTA team

dips: Detrending strictly periodic signals
A. Prsa
dips is an algorithm for detrending timeseries of strictly periodic signals. It does not assume any functional form for the signal or the background or the noise; it disentangles the strictly periodic component from everything else. We use it in astronomy for detrending Kepler, K2 and TESS timeseries of periodic variable stars, eclipsing binary stars, exoplanets etc. The algorithm is described in detail in Prša, Zhang and Wells (2019), PASP 131, 8001. A new, generalized version of dips will be explained in Horvat and Prša (2022), currently in preparation.

RADIS: Fast Line-by-line code for infrared emission & absorption spectra at equilibrium & non-LTE
Pannier, E. ; Bekerom D. v. d ; Minesi N. and the RADIS contributors
RADIS uses a new algorithm that can resolve spectra with millions of lines within seconds on a single-CPU, and can be GPU-accelerated for almost-instant-computation (up to 5e14 lines*spectral points/s). It supports HITRAN, HITEMP and ExoMol out-of-the-box (auto-download), and therefore is particularly suitable to compute cross-sections or transmission spectra at high-temperature. It includes equilibrium calculations for all species, and non-LTE for CO2 and CO. The code is an open-source Python library (https://github.com/radis/radis) and can also be executed in an online environment with pre-configured HITEMP databases (https://radis.github.io/radis-lab/).