Public repository of the Cosmic Linear Anisotropy Solving System (master for the most recent version of the standard code; GW_CLASS to include Cosmic Gravitational Wave Background anisotropies; classnet branch for acceleration with neutral networks; ExoCLASS branch for exotic energy injection; class_matter branch for FFTlog)
223
stars
292
forks
source link
Errors when varying HyRec's meR and fsR parameters #440
My collaborators and I have been trying to vary the parameters meR and fsR, which can be input directly to HyRec (the ratios of the electron mass and the fine structure constant at recombination compared to today).
For certain values of these parameters, in particular when they are less than 1, running CLASS produces an error and the recombination history is not successfully computed. There are two different kinds of errors that we have encountered.
Some examples:
with meR = 0.9, fsR = 0.98
**** ERROR HAS OCCURRED in HYREC-2 ****
in interpolate_rates: TR = 0.435027 is out of range.
called from rec_swift_hyrec_dxHIIdlna
called from rec_dxHIIdlna
with meR = 0.85, fsR = 0.98
**** ERROR HAS OCCURRED in HYREC-2 ****
xHeI = -5.523619E-05 < 0 at z = 2869.987666 in rec_helium_dxHeIIdlna.
You should try and extend the hydrogen post-saha phase by increasing DXHII_MAX in history.h
My collaborators and I have been trying to vary the parameters meR and fsR, which can be input directly to HyRec (the ratios of the electron mass and the fine structure constant at recombination compared to today).
For certain values of these parameters, in particular when they are less than 1, running CLASS produces an error and the recombination history is not successfully computed. There are two different kinds of errors that we have encountered.
Some examples:
with meR = 0.9, fsR = 0.98
with meR = 0.85, fsR = 0.98
Thank you all for your time!