Cthulhu: An Open Source Molecular and Atomic Cross Section Computation Code for Substellar Atmospheres

Arnav Agrawal,Ryan J. MacDonald

DOI: https://doi.org/10.21105/joss.06894

2024-10-18

Abstract:Atmospheric studies of exoplanets and brown dwarfs are a cutting-edge and rapidly evolving area of astrophysics research. Calculating models of exoplanet or brown dwarf spectra requires knowledge of the wavelength-dependent absorption of light (cross sections) by the molecules and atoms in the atmosphere. Here we introduce Cthulhu, a pure Python package that rapidly calculates cross sections from atomic and molecular line lists. Cthulhu includes modules to automatically download molecular line lists from online databases (e.g. ExoMol and HITRAN) and compute cross sections on a user-specified temperature, pressure, and wavenumber grid. Cthulhu requires only CPUs and can run on a user's laptop (for smaller line lists with < 100 million lines) or on a large cluster in parallel (for many billion lines). Cthulhu includes in-depth Jupyter tutorials in the online documentation. Finally, Cthulhu can be used as an educational tool to demystify the process of making cross sections for atmospheric models.

Instrumentation and Methods for Astrophysics,Earth and Planetary Astrophysics,Solar and Stellar Astrophysics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to simplify and improve the accessibility and efficiency of calculating molecular and atomic absorption cross - sections in order to support the research of exoplanet and brown dwarf atmospheric models. Specifically:

1. **Background problems**:
- With the emergence of powerful new telescopes (such as the James Webb Space Telescope (JWST) and the upcoming Extremely Large Telescopes (ELTs)), scientists can obtain spectral data of exoplanets and brown dwarfs in more detail, thereby studying their chemical compositions and physical properties.
- Calculating the atmospheric spectral models of these celestial bodies requires accurate wavelength - dependent absorption cross - section data. If these data are inaccurate, the measured chemical abundances and other atmospheric properties may be biased.

2. **Existing challenges**:
- Absorption cross - sections are usually pre - calculated from line lists in large quantum - mechanical transition databases (such as ExoMol, HITRAN, HITEMP, and VALD). However, this process is computationally intensive and complex and requires specialized tools.
- For many exoplanet and brown dwarf researchers, calculating molecular and atomic absorption cross - sections is a non - specialized area with a relatively high technical threshold.

3. **Paper solutions**:
- **Cthulhu** is a pure - Python package aimed at rapidly calculating absorption cross - sections from atomic and molecular line lists. It has the following features:
- Automatically download molecular line lists from online databases and re - format them into efficiently stored HDF5 files.
- Support users to specify temperature, pressure, and wavenumber grids to calculate absorption cross - sections.
- It can run only with a CPU, can handle smaller line lists (< 100 million) on an ordinary laptop, and can also process billions of line lists in parallel on a large cluster.
- Provide detailed Jupyter tutorials to help users learn how to use this tool.

4. **Target users**:
- Cthulhu is not only for scientific researchers but also serves as an educational tool, aiming to lower the threshold for beginners to enter the field of absorption cross - section calculation, making it easier for them to understand and operate relevant databases and tools.

By providing such a user - friendly tool, Cthulhu hopes to promote more researchers to accurately and efficiently conduct research on exoplanet and brown dwarf atmospheric models.