Sledgehamr (ScaLar fiEld Dynamics Getting solvEd witH Adaptive Mesh Refinement) simulates the dynamics of coupled scalar fields on a 3-dimensional mesh. Adaptive mesh refinement (AMR) can boost performance if spatially localized regions of the scalar field require high resolution. sledgehamr is compatible with both GPU and CPU clusters, and, because it is AMReX-based (ascl:2409.012), offers a flexible and customizable framework. This framework enables various applications, such as the generation of gravitational wave spectra.

cosmocnc evaluates the number count likelihood of galaxy cluster catalogs. Fast Fourier Transform (FFT) convolutions are used to evaluate some of the likelihood integrals. The code supports three types of likelihoods (unbinned, binned, and an extreme value likelihood); it also supports the addition of stacked cluster data (e.g., stacked lensing profiles), which is modeled in a consistent way with the cluster catalog. The package produce mass estimates for each cluster in the sample, which are derived assuming the hierarchical model that is used to model the mass observables, and generates synthetic cluster catalogs for a given observational set-up. cosmocnc interfaces with the Markov chain Monte Carlo (MCMC) code Cobaya (ascl:1910.019), allowing for easy-to-run MCMC parameter estimation.

SZiFi (pronounced "sci-fi") implements the iterative multi-frequency matched filter (iMMF) galaxy cluster finding method. It can be used to detect galaxy clusters with mm intensity maps through their thermal Sunyaev-Zeldovich (tSZ) signal. As a novel feature, SZiFi can perform foreground deprojection via a spectrally constrained MMF or sciMMF, and can also be used for point source detection.

hmvec is a pure Python/numpy vectorized general halo model and HOD code. It includes support for 3d power spectra involving NFW, Battaglia electron density profiles and galaxy HODs. It also supports 2d power spectra including tSZ, cosmic shear, galaxy-galaxy lensing and CMB lensing. hmvec calculates a vectorized FFT for a given profile over all points in mass and redshift, using one double loop over mass and redshift to interpolate the profile Fourier transforms to the target wavenumbers; every other part of the code is vectorized.

Spectool is a toolkit designed for processing astronomical spectral data, offering a collection of common spectral analysis algorithms. The package includes functions for spectral resampling, spectral flattening, radial velocity measurements, spectral convolution broadening, and more. Each function in the package is implemented independently, allowing users to select and utilize the desired features as needed. The functions are designed with simple and intuitive interfaces, ensuring ease of use for various data sets and analysis tasks.

Hierarchical Semi-Sparse Cube (HiSS-Cube) framework provides highly parallel processing of combined multi-modal multi-dimensional big data. The package builds a database on top of the HDF5 framework which supports parallel queries. A database index on top of HDF5 can be easily constructed in parallel, and the code supports efficient multi-modal big data combinations. The performance of HiSS-Cube is bounded by the I/O bandwidth and I/O operations per second of the underlying parallel file system; it scales linearly with the number of I/O nodes and can be extended to any kind of multidimensional data combination and information retrieval.

NbodyGradient computes gradients of N-body integrations for Newtonian gravity and arbitrary N-body hierarchies. Developed for transit-timing analyses and written in Julia, NbodyGradient gives derivatives of the transit times with respect to the initial conditions, either masses and Cartesian coordinates/velocities or orbital elements.

The spectools_ir suite analyzes medium/high-resolution IR molecular astronomical spectra. It has three main sub-modules (flux_calculator, slabspec, and slab_fitter) and also offers a sub-module (utils) with a few additional functions. Written with infrared medium/high-resolution molecular spectroscopy in mind, spectools_ir generally assumes spectra are in units of Jy and microns and uses information from the HITRAN molecular database. Some routines are more general, but users interested in other applications should proceed with caution.

ROCKE-3D (Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics) models the atmospheres and oceans of solar system and exoplanetary terrestrial planets. Written in Fortran, it is a three-dimensional General Circulation Model (GCM). ROCKE-3D requires Panoply, the SOCRATES radiation code and spectral files, and has several additional dependencies.

Deep-Transit detects transits using a deep learning based 2D object detection algorithm. The code determines the light curve and outputs the transiting candidates' bounding boxes and confidence scores. It has been trained for Kepler and TESS data, and can be extended to other photometric surveys and even ground-based observations. Deep-Transit also provides an interface for training new datasets.

Spinifex is a pure Python tooling for ionospheric corrections in radio astronomy, e.g. getting total electron content and rotation measures.

This repository implements an optimized XGBoost-based framework for photometric classification of Type Ia supernovae, addressing class imbalance through PR-AUC and F1-score prioritization. The approach is designed for scalability in large-scale astronomical surveys such as LSST and ensures improved classification robustness compared to traditional metrics like ROC-AUC.

exoscene simulates direct images of exoplanetary systems. Written in Python, the package has three modules. These modules can determine a planet's relative astrometry ephemeris, its phase function, and flux ratio, compute the band-integrated irradiance of a star, and accurately resample an image model array to a detector array. exoscene also offers modeling and mapping functions and has additional capabilities.

AESTRA (Auto-Encoding STellar Radial-velocity and Activity) uses deep learning for precise radial velocity measurements in the presence of stellar activity noise. The architecture combines a convolutional radial-velocity estimator and a spectrum auto-encoder. The input consists of a collection of hundreds or more of spectra of a single star, which span a variety of activity states and orbital motion phases of any potential planets. AESTRA does not require any prior knowledge about the star.

APPLESOSS (A Producer of ProfiLEs for SOSS) builds 2D spatial profiles for the first, second, and third diffraction orders for a NIRISS/SOSS GR700XD/CLEAR observation. The profiles are entirely data driven, retain a high level of fidelity to the original observations, and can be used as the specprofile reference file for ATOCA (ascl:2502.016). They can also be used as a PSF weighting for optimal extractions.

ATMOSPHERIX reads t.fits files from the Canada-France-Hawaii Telescope's near-infrared spectropolarimeter SPIRou, processes the data to remove telluric/stella contributions, and performs the correlation analysis for a given planet atmosphere template. The correlation function computes the correlation between the data and model for a grid of planet velocimetric semi-amplitude and systemic velocity. ATMOSPHERIX takes transmission spectroscopy into account and allows the user to inject a synthetic planet if desired.

Bioverse assesses the diagnostic power of a statistical exoplanet survey of the properties of nearby terrestrial exoplanets via direct imaging or transit spectroscopy. 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. The code contains a versatile module for population-level hypothesis testing supporting trade studies and survey optimization. Bioverse supports direct imaging or transit missions, and its modularity makes it adaptable to any mission concept that makes measurements on a sample of exoplanets.

CROCODILE (CROss-COrrelation retrievals of Directly-Imaged self-Luminous Exoplanets) runs atmospheric retrievals of directly observed gas giant exoplanets by adopting adequate likelihood functions. The code makes use of petitRADTRANS (ascl:2207.014) and PyMultiNest (ascl:1606.005) and provides a statistical framework to interpret the photometry, low-resolution spectroscopy, and medium (and higher) resolution cross-correlation spectroscopy.

DIA (Delta function Difference Imaging Code) provides a difference image analysis pipeline that employs a delta-function kernel; this is useful for reducing TESS Full Frame Images. DIA's scripts are available in both Python and IDL and are nearly identical in their outputs. Together, the scripts make a pipeline that cleans and aligns images, generates a master frame by combining all available images, performs image subtraction, generates light curves, and does a basic detrending to the light curves based on magnitude. DIA can also apply bias subtraction, flat fielding, background subtraction and align images to the first image in the list.

Given mass, radius, and equilibrium temperature, ExoMDN can deliver a full posterior distribution of mass fractions and thicknesses of each planetary layer. A machine-learning model for the interior characterization of exoplanets based on Mixture Density Networks (MDN), ExoMDN is trained on a large dataset of more than 5.6 million synthetic planets below 25 Earth masses. These synthetic planets consist of an iron core, a silicate mantle, a water and high-pressure ice layer, and a H/He atmosphere. ExoMDN uses log-ratio transformations to convert the interior structure data into a form that the MDN can easily handle.

GEOCLIM.jl, written in Julia, replicates some features of the original GEOCLIM model written in Fortran. It also extends the original weathering equations WHAK and MAC, which ignore direct dependence on pCO2 and include direct pCO2 dependence respectively. The code estimates global silicate weathering rates from gridded climatology. GEOCLIM.jl estimates weathering during periods of Earth history when the continental configuration was radically different, typically more than 100 million years ago, and includes functions to compute, for example, land/ocean fraction, area-weighted average, area-weighted sum, and land mass perimeter, among other values.

Gollum performs spectral visualization and analysis. It offers both a programmatic interface and a visual interface that help users analyze stellar and substellar spectra, with support included for a set of precomputed synthetic spectral model grids.

GPS (Genesis Population Synthesis) develops population synthesis models. The code suite uses the Genesis database of planet formation models for small exoplanets (super-Earths and Mini-Neptunes). Although the codebase focuses on the Genesis models, aother models can easily be integrated with GPS. It computes the bulk compositions of the planets and simulates atmospheric loss and evolution to find the final states of the planets that can be observationally verified. GPS also offers tools to process and analyze the data from recent observations of small exoplanets in order to compare them with the models.

The IGRINS_transit data reduction pipeline takes high-resolution observations of transiting exoplanets with Gemini-S/IGRINS and produces cross-correlation detections of molecules in the exoplanet's atmosphere. IGRINS_transit removes low signal-to-noise orders, performs a secondary wavelength calibration, and uses a singular value decomposition (SVD) to separate out the signature of the transiting planet from the host star and telluric contamination.

IsoFATE (Isotopic Fractionation via ATmospheric Escape) models mass fractionation resulting from diffusive separation in escaping planetary atmospheres and numerically computes atmospheric species abundance over time. The model is tuned to sub-Neptune sized planets with rocky cores of Earth-like bulk composition and primordial H/He atmospheres. F, G, K, and M type stellar fluxes are readily implemented. IsoFATE has two versions, the first of which simulates a ternary mixture of H, He, and D (deuterium); the second version is coupled to the magma ocean-atmosphere equilibrium chemistry model Atmodeller.

kpic_pipeline reduces data taken with the Keck Planet Imager and Characterizer (KPIC). Written in Python, the code processes high resolution spectroscopy data taken with KPIC to study exoplanet atmospheres; it processes and calibrate the data to enable spectroscopic model fitting. kpic_pipeline can reduce the observed data into 1D spectra for one given science target or can be used to reduce the full nightly data.

luas builds Gaussian processes (GPs) primarily for two-dimensional data sets. It uses different optimizations to make the application of GPs to 2D data sets possible within a reasonable timeframe. The code is implemented using Jax (ascl:2111.002), which helps calculate derivatives of the log-likelihood as well as permitting the code to be easily run on either CPU or GPU. luas can be used with popular inference frameworks such as NumPyro and PyMC. The package makes it easier to account for systematics correlated across two dimensions in data sets, in addition to being helpful for any other applications (e.g., interpolation).

Mini-chem solves chemical kinetics for gas giant atmospheric modeling. It is pared down from large chemical networks to make use of "net forward reaction tables"; this significantly reduces the number of reactions and species required to be evolved in the ODE solvers. The code's NCHO network consists of only 12 species with 10 reactions, making it a lightweight and easy to couple network to large scale 3D GCM models, or other models of interest (such as 1D or 2D kinetic modelling efforts). Mini-chem is written in Fortran and has three main parts: the input routine, the chemistry routines, and the still ODE solver.

NcorpiON integrates collisional and fragmenting systems of planetesimals or moonlets orbiting a central mass. It features a fragmentation model, based on crater scaling and ejecta models, that realistically simulates a violent impact. Written in C, the code detects collisions, computes mutual gravity, and can resolve a collision by fragmentation. The fast multipole expansions are implemented up to order six to allow for a high precision in mutual gravity computation.

This library of synthetic X-ray spectra provides a tabulated version of the slim disk model for fitting tidal disruption events (TDEs). The library is created by ray-tracing stationary, general relativistic slim disks and consistently incorporating gravitational redshift, Doppler, and lensing effects.

S³Fit fits spectrum and multi-band photometric Spectral Energy Distribution (SED) simultaneously for analyzing observational data of galaxies. It improves the moderate constraints on properties of continuum models in a pure spectral fitting due to the limited wavelength coverage. The code supports multiple models with multiple components, and can handle complex systems with a mixed contribution of Active Galactic Nucleus (AGN) and its host galaxy in both of continua and emission lines (e.g., narrow lines and broad outflow lines). The fitting strategy is optimized to enable an efficient solution of the best-fit results for several tens of parameters and model components. S³Fit is also extensible for adding functions and components by users such as new band filters, star formation history functions, emission lines, and also types of models.

LESSPayne performs semi-automatic analysis for echelle spectra of stars. It uses a neural network emulator to do a full spectrum fit to estimate stellar parameters and performs automatic continuum and equivalent width fits normalization with theoretical masks. The code uses MOOG (ascl:1202.009) for spectrum synthesis fitting, ATLAS model atmosphere interpolation, and equivalent width abundance determination. LESSPayne can also perform automatic abundance uncertainty analysis with error propagation and summary tables, and should be viewed as providing a high-quality initialization for an smhr file that reduces the time for a standard analysis.

StellarSpecModel interpolates the stellar spectral grid; provided with stellar parameters (Teff, FeH, logg), the package will return the corresponding stellar spectrum. It also generates and analyzes theoretical stellar spectral energy distributions (SEDs). StellarSpecModel includes functionality for both single and binary star systems, incorporating extinction models and the ability to handle photometric data in various filter bands.

GaMPEN (Galaxy Morphology Posterior Estimation Network) estimates robust posteriors (i.e., values + uncertainties) for structural parameters of galaxies using a Bayesian machine learning framework. The code also automatically crops input images to an optimal size before structural parameter estimation. The package produces extremely well-calibrated (less than 5% deviation) predicted posteriors; these have been shown to be up to 60% more accurate compared to the uncertainties predicted by many light-profile fitting algorithms. Once trained, it takes GaMPEN less than a millisecond to perform a single model evaluation on a CPU. Thus, GaMPEN’s posterior prediction capabilities are ready for large galaxy samples expected from upcoming large imaging surveys, such as Rubin-LSST, Euclid, and NGRST.

GaMorNet classifies galaxies morphologically using a Convolutional Neural Network. The code does not need a large amount of training data, as it is trained on simulations and then transfer-learned on a small portion of real data, and can be applied on multiple datasets. The software has a misclassification rate of less than 5%. GaMorNet is written in Python and uses the Keras and TFLearn deep learning libraries to perform all of the machine learning operations. Both these aforementioned libraries in turn use TensorFlow for their underlying tensor operations.

Superbubble Finding Algorithm identifies superbubbles in HI column density maps of both observed and simulated galaxies that has only two adjustable parameters. The algorithm takes an input column density galaxy image and returns the labels and basic measurements of the detected bubbles so plotting and further analysis can be done. The algorithm includes an automated galaxy-background separation step to focus the analysis on the galactic disk. Functions to solve for the superbubble radii, superbubble galactic radii location, and an external function to plot the detected bubbles are also packaged in Superbubble Finding Algorithm.

FELINE (Find Emission LINEs) combines a fully parallelized galaxy line template matching with the matched filter approach for individual emission features of LSDcat (ascl:1612.002). For the 3D matched filtering, the complete data cube is first median filtered to remove all continuum sources, and then cross-correlated with a template of an isolated emission feature in two spatial and one spectral dimension. FELINE then evaluates the likelihood in each spectrum of the cube for emission lines at the positions provided by a given redshift and a certain combination of typical emission features.

ExoSim 2 (Exoplanet Observation Simulator 2) makes spectro-photometric observations of transiting exoplanets from space, ground, and sub-orbital platforms. It is a complete rewrite of ExoSim (ascl:2002.008); it is implemented in Python 3 and uses object-oriented design principles. The package follows a three-step workflow: the creation of focal planes, the production of Sub-Exposure blocks, and the generation of non-destructive reads (NDRs). ExoSim 2 has demonstrated consistency in estimating photon conversion efficiency, saturation time, and signal generation. The simulator has also been validated independently for instantaneous read-out and jitter simulation, and for astronomical signal representation.

The Blooming Tree (BT) algorithm identifies clusters, groups, and substructures from galaxy redshift surveys. Based on the hierarchical clustering method, it takes the projected binding energy as the linking length and provides three main analysis approaches to trim the hierarchical tree; 1.) the direct trimming (binding energy, velocity disperion, or eta); 2.) the sigma plateau, when no trimming threshold is appointed; and 3.) the blooming tree. The Blooming Tree algorithm tool works only in the terminal.

XGPaint, written in Julia, generates maps of extragalactic foregrounds, using astrophysical models designed to replicate the statistics of the millimeter sky. The code computes simulated galaxies from the Cosmic Infrared Background (CIB), radio galaxies, and contributions and distortions from the Sunyaev-Zeldovich (SZ) effect. XGPaint is multithreaded, and supports both HEALPix and Plate Carrée pixelizations.

Colume (COLUMn to vOLUME) uses the statistical and spatial distribution of a column density map to infer a likely volume density distribution along each line of sight. The Python package incorporates all pre-processing (in particular re-sampling) functions needed to efficiently work on the column density maps. Colume's outputs are saved in Numpy format.

Gradus.jl traces geodesics and calculates observational signatures of accreting compact objects. The code requires only a specification of the non-zero metric components of a chosen spacetime in order to solve the geodesic equation and compute a wide variety of trajectories and orbits. Gradus includes algorithms for calculating physical quantities are implemented generically, so they may be used with different classes of spacetime with minimal effort.

NRPyElliptic sets up initial data (ID) for numerical relativity (NR) using the same numerical methods employed for solving hyperbolic evolution equations. The code implements a hyperbolic relaxation method to solve complex nonlinear elliptic PDEs for NR ID. The hyperbolic PDEs are evolved forward in (pseudo)time, resulting in an exponential relaxation of the arbitrary initial guess to a steady state that coincides with the solution of the elliptic system. The package solves these equations on highly efficient numerical grids exploiting underlying symmetries in the physical scenario. NRPyElliptic is built in the NRPy+ (ascl:1807.025) framework, which facilitates the solution of hyperbolic PDEs on Cartesian-like, spherical-like, cylindrical-like, or bispherical-like numerical grids.

SCONE (Supernova Classification with a Convolutional Neural Network) classifies supernovae (SNe) by type using multi-band photometry data (lightcurves) using a convolutional neural networks. SCONE takes in supernova (SN) photometry data in the format output by SNANA simulations, separated into two types of files: metadata and observation data. Photometric data is pre-processed via 2D Gaussian process regression, which smooths over irregular sampling rates between filters and also allows SCONE to be independent of the filter set on which it was trained.

R2D2 (Residual-to-Residual DNN series for high-Dynamic range imaging) performs synthesis imaging for radio interferometry. The R2D2 algorithm takes a hybrid structure between a Plug-and-Play (PnP) algorithm and a learned version of the well-known Matching Pursuit algorithm. Its reconstruction is formed as a series of residual images, iteratively estimated as outputs of iteration-specific Deep Neural Networks (DNNs), each taking the previous iteration’s image estimate and associated back-projected data residual as inputs. The primary application of the R2D2 algorithm is to solve large-scale high-resolution high-dynamic range inverse problems in radio astronomy, more specifically 2D planar monochromatic intensity imaging.

LeR calculates detectable rates of gravitational waves events (both lensed and un-lensed events). Written in Python, it performs statistical simulation and forecasting of gravitational wave (GW) events and their rates. The code samples gravitational wave source properties and lens galaxies attributes and source redshifts, and can generate image properties such as source position, magnification, and time delay. The package also calculates detectable merger rates per year. Key features of LeR include efficient sampling, optimized SNR calculations, and systematic archiving of results. LeR is tailored to support both GW population study groups and GW lensing research groups by providing a comprehensive suite of tools for GW event analysis.

A Python package for helping with detecting exoplanet transit dips in TESS light curves.

Features:
Load and clean .fits light curve or target pixel files
Visualize light curves and transit candidates
Designed for students, researchers, and citizen scientists
Easily download and process raw data in bulk
Easily shortlist potential exoplanet candidates

nuance uses linear models and Gaussian processes to simultaneously search for planetary transits while modeling correlated noises (e.g., stellar variability). The code computes the likelihood of a transit being present in some correlated noise without disentangling the two; it searches the transit signal while simultaneously modeling correlated noise, assuming that the light curve can be modeled as a Gaussian process. nuance detects single or periodic transits and can find transits in light curves from multiple instruments, whether space-based or ground-based; it can also run in parallel on CPUs or GPUs.

PBjam analyzes the oscillation spectra of solar-like oscillators. The code performs two main tasks: identifying a set of modes of interest in a spectrum of oscillations, and accurately modeling those modes to measure their frequencies. Mode identification relies on a large set of previous observations of the model parameters, which are then used to construct a prior distribution to inform the sampling. PJjam models the modes using a nested sampling or MCMC algorithm, where Lorentzian profiles are fit to each of the identified modes.

The 1D radiative-convective code PCM_LBL simulates the climates of diverse planetary atmospheres. The code is written in modular modern Fortran and uses a 'brute-force' spectral approach where absorption coefficients are computed on a fixed spectral grid directly from line data. This allows climate calculations to be performed more simply and at higher accuracy than in a correlated-k approach. PCM_LBL allows the user to iterate rapidly between fast, lower accuracy calculations and slow high accuracy calculations. By default, the model is set up to run fairly fast at moderate resolution; the accuracy of the code can be adjusted with a few (documented) changes.

photoevolver simulates the atmospheric escape of extrasolar planets and their evolution. The code evolves the gaseous atmosphere of a planet backwards and forwards in time, taking into account its internal structure and cooling rate, atmospheric mass loss processes, and the stellar emission history. photoevolver determines whether a palent's atmosphere survives or ise completely stripped by radiation from its host star.

plaNETic uses a Bayesian neural network-based to model small (masses between 0.5 and 15 M_earth) exoplanets. The code efficiently computes posteriors of a planet's internal structure based on its observed planetary and stellar parameters. It uses a full grid accept-reject sampling algorithm. plaNETic also allows for different choices in priors concerning the expected abundance of water (formation inside vs. outside of iceline) and the planetary Si/Mg/Fe ratios (stellar vs. iron-enriched vs. free).

pycdata imports datasets from various telescopes/instruments in pycheops (ascl:2312.034), thus providing the facility pycheops lacks to model transits, eclipses, phase curves from other telescopes/instruments and the PSF photometry produced by PIPE (ascl:2404.002). pycdata automatically puts resultant data products into the pycheops cache directory so that it can be directly readable from the pycheops command line.

RTModel models and interprets microlensing events. It uses photometric time series collected from ground and/or space telescopes to propose one or more of the following possible models:

- single-lens-single-source microlensing;
- single-lens-binary-source microlensing, with or without xallarap; and/or
- binary-lens-single-source microlensing, including planetary microlensing, parallax and orbital motion.

All models include the finite-size of the source(s). The modeling strategy is based on a grid search in the parameter space for single-lens models, whereas a template library for binary-lens models is used including all possible geometries of the source trajectory with respect to the caustics. In addition to this global search, planets are searched where maximal deviations from a Paczynski model occurs. The RTModel package also including subpackages for creating an immediate visualization of models and the possibility to review each individual fitting process as an animated GIF.

Spright predicts planetary masses, densities, and radial velocity semi-amplitudes given a small planet's radius or planetary radii given the small planet's mass. The package contains two relations: one for small planets orbiting M dwarfs and another for planets orbiting FGK stars. The radial velocity semi-amplitude can be predicted given the planet's radius, orbital period, orbital eccentricity (optional), and the host star mass. Spright offers both a command line script and a set of Python classes. The command line script can directly create publication-quality plots, and the classes offer a full access to the predicted numerical distributions.

SWAMPE models the dynamics of exoplanetary atmospheres; it 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 can be modified to model dissimilar space objects, from Brown Dwarfs to terrestrial, potentially habitable exoplanets. SWAMPE can be easily run on a personal laptop.

tglc (TESS-Gaia Light Curve) produces PSF-based TESS full-frame image (FFI) light curves for any sector and any star with custom options. Using Gaia DR3 as priors, the code has forward modeled FFIs with the effective point spread function to remove contamination from nearby stars. The resulting light curves show a photometric precision closely tracking the pre-launch prediction of the noise level: TGLC's photometric precision consistently reaches <2% at 16th TESS magnitude even in crowded fields, demonstrating excellent decontamination and deblending power.

TriArc (reTRIeval ARCturus) uses Bayesian statistics to determine the minimum abundance of an atmospheric species in a given model atmosphere (excluding the species of interest) and spectral noise profile. The code is configured for transmission spectroscopy and is built using the forward modeling capabilities of petitRADTRANS (ascl:2207.014); it also uses PandExo (ascl:1906.016). TriArc has been used to calculate prebiosignature detection thresholds for various potential JWST targets.

Turbospectrum_NLTE updates the spectral synthesis code Turbospectrum (ascl:1205.004) with NLTE capabilities. The code takes a 1D model atmosphere, one or several line lists, and computes the emergent spectrum (flux and/or intensities at various angles), with a prescribed chemical composition. Various parameters can be adjusted, such as microturbulence (vmicro), individual abundances, and isotopic ratios. Turbospectrum_NLTE can also handle the computation for a single chunk of spectrum with a constant wavelength step, or for a number of smaller windows, e.g., around lines of interest. Calculations can also be done for plane-parallel or spherically symmetric models.

VSPEC (Variable Star PhasE Curve) simulates exoplanet observations. It combines NASA’s Planetary Spectrum Generator (PSG) with a custom variable star model. Originally built to simulate the infrared excess of non-transiting planets, VSPEC supports transit, eclipse, phase curve geometries as well as spots, faculae, flares, granulation, and the transit light source effect.

CHIMERA (CaltecH Inverse ModEling and Retrieval Algorithms) retrieves exoplanet atmospheres, and can be used for both transmission and emission geometries with options for both the "free" and "chemically consistent" abundance retrievals. The code uses correlated-K opacities (R=100) with the random-overlap resort-rebin procedure and includes full multiple scattering in emission (both planetary and stellar reflected light) using a two stream approximation variant. CHIMERA includes multiple Bayesian samplers, including PyMultiNest (ascl:1606.005) and dynesty (ascl:1809.013).

DYNAMITE (DYNAmical Multi-planet Injection TEster) predicts the presence of unseen planets in multi-planet systems via population statistics. The code uses the specific (yet often incomplete) data on the orbital and physical parameters for the planets in any given system's architecture and combines it with detailed statistical population models and a dynamical stability criterion to predict the likelihood for the parameters of one additional planet in the system. DYNAMITE's predictions are given in the form of observable values (transit depth measurements, RV semi-amplitudes, or direct imaging separation and contrast), which can be tested by follow-up observations.

TROPF (Tidal Response Of Planetary Fluids) enables efficient terrestrial fluid tidal studies across a wide range of parameter space. The software includes several different solutions to the governing equations in classical tidal theory, and can calculate millions of such solutions on several-minute-long timescales. Written in MATLAB/Octave, TROPF can be ported to Python and other languages, as the instructions for building the operator matrices are described in detail and the coding of core TROPF routines adheres to generic sparse matrix operations and avoids functions specific to MATLAB.

Pytmosph3R computes transmission and emission spectra based on 3D atmospheric simulations, for example, performed with the LMDZ generic global climate model. It produces transmittance maps of the atmospheric limb at all wavelengths that can then be spatially integrated to yield the transmission spectrum. Pytmosph3R can use 3D time-varying atmospheric structures from a GCM as well as simpler, parameterized 1D or 2D structures, and can be used in notebooks or on the command line.

tpfplotter creates a TESS Target Pixel File of a source, overplotting the aperture mask used by the SPOC pipeline and the Gaia catalog to check for possible contaminations within the aperture. The software can create 1-column paper-ready figures, overplotting the Gaia DR2 catalog to the TESS Target Pixel Files, and can create plots for any target observed by TESS. tpfplotter can search by coordinates if the TIC number of the source is not known.

ExoInt devolatilizes stellar abundances to produce rocky exoplanetary bulk composition to constrain the modeling of the exoplanet interiors; the code uses Monte Carlo simulations that assume that each element’s abundance (within its uncertainty) follows a Gaussian distribution. ExoInt also contains a module to provide the mineralogy based on the stoichemitric output of mantle and core compositions, core mass fraction, along with the given mass and radius information.

ChromaStarPy computes the vertical structure of a static, plane-parallel, one-dimensional stellar atmosphere in local thermodynamic equilibrium (LTE); it also computes the emergent spectrum incorporating opacity computed with a comprehensive atomic line list from the NIST Atomic Spectra Database. The code provides post-processed data products that are ready to visualize in a Python IDE such as spyder. ChromaStarPy is a port of ChromaStarServer (ascl:1701.009); the code enables users to experiment with and develop a stellar astrophysical modeling code in a graphical IDE, and to compare observational data to ad hoc model output.

GalClass facilitates visual morphological classifications of large samples of galaxies taking advantage of multi-wavelength imaging and ancillary information. It offers a versatile Graphic User Interface (GUI), which adapts to the provided classification scheme. It displays a series of pre-prepared PDF files for classification, grouping by galaxy and filter, while also listing relevant metadata and displaying a color image of each source. GalClass enables easy navigation through the sample and continuously outputs classification results in a JSON file. Finally, it offers an analysis submodule which combines and processes output files of multiple classifications.

MultiREx generates synthetic transmission spectra of exoplanets. This tool extends the functionalities of the TauREx (ascl:2209.015) framework, enabling the mass production of spectra and observations with added noise. Though the package was originally conceived to train machine learning models in the identification of biosignatures in noisy spectra, it can also be used for other purposes.

AstroPT trains astronomical large observation models using imagery data. The code follows a similar saturating log-log scaling law to textual models and the models' performances on downstream tasks as measured by linear probing improves with model size up to the model parameter saturation point. Other modalities can be folded into the AstroPT model, and use of a causally trained autoregressive transformer model enables integration with the wider deep learning FOSS community.

PDQ predicts the positions on the sky of high-redshift quasars that should provide photons that are both acausal and uncorrelated. The predicted signal-to-noise ratios are calculated at framerate sufficient for random-number generation input to a loophole-free Bell test, and are calibrated against a public archival dataset of four pairs of highly-separated bright stars observed simultaneously (and serendipitously) at 17 Hz with that same instrumentation in 2019 to 2021.

RFIClean excises periodic RFI (broadband as well as narrow-band) in the Fourier domain, and then mitigates narrow-band spectral line RFI as well as broadband bursty time-domain RFI using robust statistics. Primarily designed to efficiently search and mitigate periodic RFI from GMRT time-domain data, RFIClean has evolved to mitigate any spiky (in time or frequency) RFI as well, and from any SIGPROC filterbank format data file. RFIClean uses several modules from SIGPROC (ascl:1107.016) to handle the filterbank format I/O.