pahprops.c File Reference

#include "pahprops.h"
#include <gsl/gsl_math.h>

Include dependency graph for pahprops.c:

Functions
double	bndust_PAHNctoHc (double Nc)
double	bndust_PAHatoNC (double a)
double	bndust_PAHatoHc (double a)
void	dust_CPAHDrudes (double a, dust_ionisation state, dust_pahenhancements *e, dust_pahprofiles res)
double	bndust_PAHCutoff (double lambda, double a, dust_ionisation state)
	The cutoff of PAH absorption. See Eq A2 and A3 of the Li & Draine paper.
Variables
dust_drudeprofs	dust_DL_pahs_neutral [14]
dust_drudeprofs	dust_DL_pahs_ionised [14]

---

Detailed Description

Bojan Nikolic <bojan@bnikolic.co.uk>, <b.nikolic@mrao.cam.ac.uk>

---

Function Documentation

double bndust_PAHatoHc

(

double

a

)

Returns the H/C ratio given a grain radius (calls bndust_PAHNctoHc).

References bndust_PAHatoNC(), and bndust_PAHNctoHc().

Referenced by dust_CPAHDrudes(), and dust_nCHmodes().

Here is the call graph for this function:

double bndust_PAHatoNC

(

double

a

)

Returns the number of carbon atoms in a grain of radius a (a in in metres).

note a must be metres

Referenced by bndust_PAHatoHc(), bndust_PAHCutoff(), dust_nCCIPmodes(), dust_nCCOPmodes(), and Dust::GenPAHCC::GenPAHCC().

double bndust_PAHNctoHc

(

double

Nc

)

Returns an approximate ratio of number of hydrogren atoms to carbon atoms given the number of carbon atoms. See eq 4 of 2001ApJ...554..778L

Referenced by bndust_PAHatoHc().

void dust_CPAHDrudes	(	double	a,
		dust_ionisation	state,
		dust_pahenhancements *	e,
		dust_pahprofiles	res
	)

Calculate the appropriate drude profiles for a particular PAH. a is the grain radius

References bndust_PAHatoHc(), dust_pahenhancements::E62, dust_pahenhancements::E77, dust_pahenhancements::E86, dust_drudeprofs::gamma, dust_drudeprofs::sigma, and dust_drudeprofs::wavel.

Here is the call graph for this function:

---

Variable Documentation

dust_drudeprofs dust_DL_pahs_ionised[14]

Initial value:

{
  { .0722 ,     0.195,  7.97e7 },
  {.2175  ,     0.217 , 1.23e7 },
  { 3.3   ,     0.012,  44.7    },
  {6.2    ,     0.032,  157    },
  {7.7    ,     0.091,  548    },
  {8.6    ,     0.047,  242     },
  {11.3   ,     0.018,  400     },
  {11.9   ,     0.025,  61.4     },
  {12.7   ,     0.024,  149     },
  {16.4   ,     0.010,  5.52     },
  {18.3   ,     0.036,  6.04     },
  {21.2   ,     0.038,  10.8     },
  {23.1   ,     0.046,  2.78    },  
  {26.0   ,     0.69 ,  15.2    } 
}

dust_drudeprofs dust_DL_pahs_neutral[14]

Initial value:

 {
  { .0722 ,     0.195,  7.97e7},
  {.2175  ,     0.217 , 1.23e7},
  { 3.3   ,     0.012,  197 },
  {6.2    ,     0.032,  19.6 },
  {7.7    ,     0.091,  60.9 },
  {8.6    ,     0.047,  34.7 },
  {11.3   ,     0.018,  427 },
  {11.9   ,     0.025,  72.7 },
  {12.7   ,     0.024,  167 },
  {16.4   ,     0.010,  5.52 },
  {18.3   ,     0.036,  6.04 },
  {21.2   ,     0.038,  10.8 },
  {23.1   ,     0.046,  2.78},  
  {26.0   ,     0.69 ,  15.2 } 
}

These are the profiles as given in Table 1 of 2001ApJ...554..778L. Wavelengths in micron. Gamma is dimensionless. Cross section in 10^-20 cm per Carbon atom.