petitRADTRANS.retrieval.chemistry

Module Contents

Functions

get_abundances(pressures, temperatures, line_species, ...)

This function takes in the C/O ratio, metallicity, and quench pressures and uses them

get_easychem_abundances(pressures, temperatures, ...)

update_atom_abundances(atom_abundances_list, ...)

get_exoatmos_abundances(pressures, temperatures, ...)

setup_exoatmos()

petitRADTRANS.retrieval.chemistry.get_abundances(pressures, temperatures, line_species, cloud_species, parameters, AMR=False)

This function takes in the C/O ratio, metallicity, and quench pressures and uses them to compute the gas phase and equilibrium condensate abundances from an interpolated table. Alternatively, a “free chemistry” approach can be used to set abundances. This function assumes a hydrogen-helium dominated atmosphere, and enforces <100% trace gas abundance by mass.

Args:
pressuresnumpy.ndarray

A log spaced pressure array. If AMR is on it should be the full high resolution grid.

temperaturesnumpy.ndarray

A temperature array with the same shape as pressures

line_speciesList(str)

A list of gas species that will contribute to the line-by-line opacity of the pRT atmosphere.

cloud_speciesList(str)

A list of condensate species that will contribute to the cloud opacity of the pRT atmosphere.

parametersdict

A dictionary of model parameters, in particular it must contain the names C/O, Fe/H and log_pquench. Additionally the cloud parameters log_X_cb_Fe(c) and MgSiO3(c) must be present.

AMRbool

Turn the adaptive mesh grid on or off. See fixed_length_amr for implementation.

Returns:
abundancesdict

Mass fraction abundances of all atmospheric species

MMWnumpy.ndarray

Array of the mean molecular weights in each pressure bin

small_indexnumpy.ndarray

The indices of the high resolution grid to use to define the adaptive grid.

PBasesdict

A dictionary of the cloud base pressures, either computed from equilibrium condensation or set by the user.

petitRADTRANS.retrieval.chemistry.get_easychem_abundances(pressures, temperatures, line_species, cloud_species, parameters, AMR=False)
petitRADTRANS.retrieval.chemistry.update_atom_abundances(atom_abundances_list, atom_name_list, parameters)
petitRADTRANS.retrieval.chemistry.get_exoatmos_abundances(pressures, temperatures, parameters, AMR=False)
petitRADTRANS.retrieval.chemistry.setup_exoatmos()