Usage

Emulchem is a python wrapper for astrochemical emulators. The module contains emulators for astrochemcial models of molecules. These can be accessed by instantiating a python object associated to the molecule as follows:

import emulchem
CS = emulchem.ChemistryEmulator(specie="CS")

To get a list of all of the molecules included in the package, the following command can be run:

emulchem.molecule_list()

Predictions can be obtained through the instantiated object as follows:

radfield = 10
zeta = 10
density = 10**5
temperature =  200
metallicity = 1
av = 10
CS.get_prediction(radfield,zeta,density,av,temperature,metallicity)

Notice that it is possible to see the units expected by the emulator (the units are consistent between molecules) by calling:

help(CS.get_prediction)

The module also contains emulators for the radiative transfer module RADEX. These are useful for prototyping but are not perfectly accurate and are run over a limited grid. As such they should not be run in final analysis!!!

Here is a more complicated example of using the module that involves running an MCMC grid using the emcee python package:

import emcee

def something():
    pass

The associated posteriors are then as follows