HITRAN/HITEMP

May (2021) Hajime Kawahara

See “Computing CO cross section using HITRAN (opacity calculator = LPF)” for tutorial to compute a cross section profile using HITRAN/HITEMP.

Molecular Database

An example to use HITRAN / HITEMP from exojax is like that.

>>> from exojax.utils.grids import wavenumber_grid
>>> from exojax.spec import moldb
>>> nus,wav,res=nugrid(22880.,23000.,1000,unit="AA")
>>> mdbCO=moldb.MdbHit('.database/05_HITEMP2019.par.bz2',nus)
bunziping

HITEMP H2O and CO2

For H2O and CO2, HITEMP provides multiple par files. To use those files, provide the directory path for moldb.MdbHit as follows.

>>> mdbH2O=moldb.MdbHit('.database/01_HITEMP2010',nus)
>>> mdbCO2=moldb.MdbHit('.database/02_HITEMP2010',nus)

extract option

extract=True in moldb.MdbHit extracts the opacity data in the wavenumber range of nus with margin. Theforefore it can reduce the use of DRAM. It may be useful for large databases such as CH4.

>>> mdbCH4=moldb.MdbHit('.database/06_HITEMP2020.par.bz2',nus,extract=True)

This creates a new directory (such as 6101.281269066504_6108.7354917532075_1.0 ) and the extracted data and header files (06_HITEMP2020.header and 06_HITEMP2020.par ) in it:

ls .database

06_HITEMP2020.par
6101.281269066504_6108.7354917532075_1.0:
06_HITEMP2020.header  06_HITEMP2020.par

Basic Quantities

These are the basic quantities of MdbHit.

quantity

instance

unit

np/jnp

line center

nu_lines

cm-1

np

line center

dev_nu_lines

cm-1

jnp

lower state energy

elower

cm-1

jnp

natural broadening

gamma_natural

cm-1

jnp

air pressure broadening

gamma_air

cm-1

jnp

self broadning

gamma_self

cm-1

jnp

Einstein coefficient

A

s-1

jnp

reference line strength

Sij0

cm

np

log_e Sij0

logsij0

jnp

statistical weight

gpp

jnp

temperature exponent

n_air

jnp

Collision Induced Absorption (CIA)

Collision Induced Absorption (CIA) is one of the important continuum in exoplanet/brown dwarfs atmosphere. HITRAN provides the CIA files . To load the CIA file, use contdb.

>>> from exojax.spec import contdb
>>> cdbH2H2=contdb.CdbCIA('.database/H2-H2_2011.cia',nus)

See “CIA coefficient” for tutorial.