Class ColumnDensity

Initializing an instance of the ColumnDensity class.

To be calculated:

• r_des: The minimum radius, ie where the dust species starts to exist. This is usually where the density of the species rises to 1% of the maximum density of the shell (in cm).
• t_des: The temperature at the destruction radius r_des of the species
• r_max: The maximum radius at which species exists, taken to be outer radius, or where the density of species drops below 10**(-10) times the maximum density in the shell (in cm).
• t_min: The temperature at the maximum radius r_max of the species
• coldens: The column density in g/cm2 of the dust species between r_des and r_max

Parameters:

• star (Star()) - The model parameter set

Overrides: object.__init__

Read all column densities, min/max temperatures and min/max radii for the species involved in the MCMax model.

Note that the self.coldens dictionary does not give real column densities! This dict merely gives column densities in a prescribed shell with given min and max radius, in order to compare with the H2 col density.

Calculate the full column density of a dust species in the shell.

This is NOT the value used for determining the abundance of the species

Parameters:

• species (string) - The dust species

Returns: float

The column density of the dust species in the full envelope in g/cm2

Calculate the full NUMBER column density of a dust species in the shell.

This is NOT the number used in determining the equivalent molecular abundance!

Parameters:

• species (string) - The dust species

Returns: float

The number column density of the dust species in the full envelope in cm-2

Calculate the molecular abundance of a dust species with respect to H2.

Note that the self.coldens dictionary does not give real column densities! This dict merely gives column densities in a prescribed shell with given min and max radius, in order to compare with the H2 col density.

Parameters:

• species (string) - the dust species (from Dust.dat)

Returns: float

The molecular abundance with respect to H2 of the dust species

Calculate the column number density of molecular hydrogen between two radial distances.

If a GASTRoNOoM model is calculated, the h2 density profile is taken from there. Otherwise, the value is derived from the total mass-loss rate, assuming everything is H2.

The use of a proper velocity profile can make a big difference for low mass-loss rate sources!

Parameters:

• species (string) - the dust species for which the H2 col dens is calculated. Required for the radial information.

Returns: float

The molecular hydrogen column number density between the given radii (cm-2)