Blame view

README.md 3.57 KB
d08e41f0   Jalabert   feat files
1
2
## Requirement

9b18e070   Jalabert   feat add and corr...
3
4
5
6
7
```
pip install numpy  
pip install astropy  
pip install matplotlib  
```
d08e41f0   Jalabert   feat files
8

9b18e070   Jalabert   feat add and corr...
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
## Programs

**radiation_fields.py** :  
*Description* :  
This program calculates the scaling factor of the radiation field G0 relative 
to the data 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

**three_levels_model.py** :  
*Description* :  
This program seeks to calculate how a neutral gas subjected to a radiation field 
is heated, involving PAHs in the ionization of the gas. The expressions used are based on the model with 3 energy levels of the molecule, going from the neutral molecule (charge state Z=0) 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,   
512b446c   Jalabert   addition of an au...
49
            electron density in cm-3 (n_e = n_h * 1.6e-4, n_h : hydrogen density)  
9b18e070   Jalabert   feat add and corr...
50
51
52
53
54
55
56
57
58
59
        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)  
            
17fb87a4   Jalabert   feat correction file
60
**four_levels_model.py** :  
9b18e070   Jalabert   feat add and corr...
61
*Description* :  
512b446c   Jalabert   addition of an au...
62
This program does the same thing and has the same parameters as the program three_levels_model.py, the difference being that we add expressions related to the 4 energy levels model. This model consists in taking into account the anions, i.e. the PAHs in a charge state Z = -1.
9b18e070   Jalabert   feat add and corr...
63
64
65

These last two programs use the parameters returned by radiations_fields.py in its calculations

384ac2cd   Jalabert   feat correction
66
## Useful files
9b18e070   Jalabert   feat add and corr...
67
68
69

2 files categories are usefull to run the programs :

a08dc7ee   Jalabert   add commentaries
70
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.  
9b18e070   Jalabert   feat add and corr...
71
72
73
74
75
76
77
78
79
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
0917f1d1   Jalabert   feat add commenta...
80
81
82
83
    2. stars file

## Authors

512b446c   Jalabert   addition of an au...
84
The authors of the programs are O. Berné, C. Joblin, S. Foschino and F. Jalabert