thermodynamics.h

Go to the documentation of this file.

 
#ifndef __THERMODYNAMICS__
#define __THERMODYNAMICS__
 
#include <vector>
 
#include "background.h"
//#include "arrays.h"
//#include "helium.h"
//#include "hydrogen.h"
 
enum recombination_algorithm { recfast, hyrec };
 
enum reionization_parametrization {
  reio_none,      
  reio_camb,      
  reio_bins_tanh, 
  reio_half_tanh, 
  reio_many_tanh, 
  reio_inter      
};
 
enum reionization_z_or_tau {
  reio_z,  
  reio_tau 
};
 
#define f1(x) \
  (-0.75 * x * (x * x / 3. - 1.) + 0.5)     
#define f2(x) (x * x * (0.5 - x / 3.) * 6.) 
#define _BBN_ -1
 
struct thermo {
 
  double YHe = _BBN_; 
  enum recombination_algorithm recombination = recfast; 
  enum reionization_parametrization reio_parametrization = reio_camb; 
  enum reionization_z_or_tau reio_z_or_tau =
      reio_z; 
  double tau_reio = 0.0925; 
  double z_reio = 11.357; 
  bool compute_cb2_derivatives =
      false; 
  bool compute_damping_scale =
      false; 
  double reionization_width = 0.5; 
  double reionization_exponent = 1.5; 
  double helium_fullreio_redshift = 3.5; 
  double helium_fullreio_width = 0.5; 
  int binned_reio_num = 0; 
  std::vector<double> binned_reio_z; 
  std::vector<double> binned_reio_xe; 
  double binned_reio_step_sharpness =
      0.3; 
  int many_tanh_num = 0; 
  std::vector<double> many_tanh_z; 
  std::vector<double>
      many_tanh_xe; 
  double many_tanh_width = 0.5; 
  int reio_inter_num = 0; 
  std::vector<double> reio_inter_z; 
  std::vector<double> reio_inter_xe; 
  double annihilation =
      0.; 
  bool has_on_the_spot =
      true; 
  double decay = 0.; 
  double annihilation_variation = 0.; 
  double annihilation_z = 1000.; 
  double annihilation_zmax = 2500.; 
  double annihilation_zmin = 30.; 
  double annihilation_f_halo =
      0.; 
  double annihilation_z_halo = 30.; 
 
 
  short thermodynamics_verbose =
      0; 
};
 
struct recombination {
 
  int index_re_z;          
  int index_re_xe;         
  int index_re_Tb;         
  int index_re_wb;         
  int index_re_cb2;        
  int index_re_dkappadtau; 
  int re_size;             
 
 
  int rt_size; 
  std::vector<double>
      recombination_table; 
 
 
  double CDB;     
  double CR;      
  double CK;      
  double CL;      
  double CT;      
  double fHe;     
  double CDB_He;  
  double CK_He;   
  double CL_He;   
  double fu;      
  double H_frac;  
  double Tnow;    
  double Nnow;    
  double Bfact;   
  double CB1;     
  double CB1_He1; 
  double CB1_He2; 
  double H0;      
  double YHe;     
  /* parameters for energy injection */
 
  double
      annihilation; 
  bool has_on_the_spot; 
  double decay; 
  double annihilation_variation; 
  double annihilation_z; 
  double annihilation_zmax; 
  double annihilation_zmin; 
  double annihilation_f_halo; 
  double annihilation_z_halo; 
};
 
struct reionization {
 
  int index_re_z;          
  int index_re_xe;         
  int index_re_Tb;         
  int index_re_wb;         
  int index_re_cb2;        
  int index_re_dkappadtau; 
  int index_re_dkappadz; 
  int index_re_d3kappadz3; 
  int re_size;             
 
 
  int rt_size; 
  std::vector<double>
      reionization_table; 
 
 
  double
      reionization_optical_depth; 
 
 
  /* parameters used by reio_camb */
 
  int index_reio_redshift; 
  int index_reio_exponent; 
  int index_reio_width; 
  int index_reio_xe_before; 
  int index_reio_xe_after;  
  int index_helium_fullreio_fraction; 
  int index_helium_fullreio_redshift; 
  int index_helium_fullreio_width; 
  /* parameters used by reio_bins_tanh, reio_many_tanh, reio_inter */
 
  int reio_num_z;         
  int index_reio_first_z; 
  int index_reio_first_xe; 
  int index_reio_step_sharpness; 
  /* parameters used by all schemes */
 
  int index_reio_start; 
 
  std::vector<double>
      reionization_parameters; 
  int reio_num_params; 
 
 
  int index_reco_when_reio_start; 
};
 
 
 
 
#define _m_e_ 9.10938215e-31  
#define _m_p_ 1.672621637e-27 
#define _m_H_ 1.673575e-27    
#define _not4_ 3.9715         
#define _sigma_ 6.6524616e-29 
 
 
#define _RECFAST_INTEG_SIZE_ 3
 
#define _Lambda_ 8.2245809
#define _Lambda_He_ 51.3
#define _L_H_ion_ 1.096787737e7
#define _L_H_alpha_ 8.225916453e6
#define _L_He1_ion_ 1.98310772e7
#define _L_He2_ion_ 4.389088863e7
#define _L_He_2s_ 1.66277434e7
#define _L_He_2p_ 1.71134891e7
#define _A2P_s_ 1.798287e9           /*updated like in recfast 1.4*/
#define _A2P_t_ 177.58e0             /*updated like in recfast 1.4*/
#define _L_He_2Pt_ 1.690871466e7     /*updated like in recfast 1.4*/
#define _L_He_2St_ 1.5985597526e7    /*updated like in recfast 1.4*/
#define _L_He2St_ion_ 3.8454693845e6 /*updated like in recfast 1.4*/
#define _sigma_He_2Ps_ 1.436289e-22  /*updated like in recfast 1.4*/
#define _sigma_He_2Pt_ 1.484872e-22  /*updated like in recfast 1.4*/
 
 
 
#define _a_PPB_ 4.309
#define _b_PPB_ -0.6166
#define _c_PPB_ 0.6703
#define _d_PPB_ 0.5300
#define _T_0_ pow(10., 0.477121) /* from recfast 1.4 */
#define _a_VF_ pow(10., -16.744)
#define _b_VF_ 0.711
#define _T_1_ pow(10., 5.114)
#define _a_trip_ pow(10., -16.306) /* from recfast 1.4 */
#define _b_trip_ 0.761             /* from recfast 1.4 */
 
 
 
#define _YHE_BIG_ 0.5    
#define _YHE_SMALL_ 0.01 
#define _Z_REC_MAX_ 2000.
#define _Z_REC_MIN_ 500.
 
 
#endif