Requirement

pip install numpy  
pip install astropy  
pip install matplotlib  

Programs

radiation_fields.py :
Description :
This program calculates the scaling factor of the radiation field G0 relative to the spectrum of an object or a set of objects radiating for given distances A file containing my G0 values with the associated distance will then be created

Input parameters :
filename : str,
name of the .txt containing intensity (in erg cm-2 s-1 nm-1 sr-1) and wavelength (in nm) data
parsec : float,
distance in parsec from the star
star_radius : float,
star radius, in unit of solar radius
ISRF : bool,
Interstellar Radiation Field; if true, then the filename is a file for ISRF
if false, the radiation field of a star is studied
RF_list: bool,
If false, the calculation of G0 is done for a given distance, if true, it is done for a hundred distances whose input distance will correspond to a G0 of the order of 1e6

For the interstellar medium, the parameters star_radius, parsec and RF_list are no longer important; the part of the program executed will be independent of them

four_levels_model.py :
Description :
This main program computes the photoelectric heating rate and efficiency by PAHs. The expressions used are based on the model with 4 charge states of the molecule, going from the anion molecule (charge state Z=-1) to the dication (charge state Z=2) by ionization.

Input parameters :
filename : str,
name of the file containing wavelength in nm and intensity of a star or of the interstellar medium
star_radius : float,
star radius, in unit of solar radius
t_gas : float,
gas temperature
n_e : float,
electron density in cm-3 (n_e = n_h * 1.6e-4, n_h : hydrogen density)
n_c : float,
number of carbon atoms in pah molecules (size of the PAH)
parsec : float,
distance in parsec from the star
fc_pah : float,
fraction of cosmic carbon locked in PAHs (default: 0.1)
ISRF : bool,
Interstellar Radiation Field; if true, then the filename is a file for ISRF if false, the radiation field of a star is studied (default: False)

This program uses radiations_fields.py.

Useful files

2 files categories are usefull to run the programs :

1. First, cross-section data (in Mb) of each PAH for sizes ranging from minimum 24 C atoms to 66 C atoms with associated energy (in eV) for each charge state of the molecule.
   Here you can find :

   1. the anions file
   2. the neutrals file
   3. the cations file
   4. the dications file

2. Then you have the files containing the data of the stars or the interstellar medium for different approximations. Each file has a column with wavelengths (in nm) associated with an intensity (in erg cm-2 s-1 nm-1 sr-1).
   Here you can find :

   1. ISRF file
   2. stars file

Authors

The authors of the programs are O. Berné, C. Joblin, S. Foschino and F. Jalabert