HyRex.hyrex module

class HyRex.hyrex.recomb_model(integration_spacing=0.0005, z0=8000.0, z1=0.0)[source]

Bases: Module

Complete recombination model implementation.

Combines helium and hydrogen recombination calculations with reionization modeling to compute full ionization history.

Attributes:

  • integration_spacing float

    Step size in log scale factor for integration

  • lna_axis_full array

    Full log scale factor grid for recombination evolution

  • z1 float

    Final redshift for evolution

  • twog_redshift float

    Redshift at which two-photon processes become negligible

  • He4equil_redshift float

    Redshift threshold for HeII+III equilibrium phase

  • idx_4He_equil array

    Indices for HeII+III equilibrium phase on full grid

Methods:

  • get_history :

    Compute complete recombination and reionization history (units: dimensionless)

get_history(h, omega_b, omega_cdm, Neff, YHe, z_reion=11, Delta_z_reion=0.5, z_reion_He=3.5, Delta_z_reion_He=0.5, exp_reion=1.5, rtol=1e-06, atol=1e-09, solver=Kvaerno3(), max_steps=1024)[source]

Compute complete recombination and reionization history.

Parameters:

  • h float

    Hubble parameter

  • omega_b float

    The baryon density Omega_b h^2

  • omega_cdm float

    The density of Cold Dark Matter Omega_cdm h^2

  • Neff float

    Effective number of neutrinos

  • YHe float

    Helium fraction

  • z_reion float, optional

    Reionization redshift of hydorgen and neutral helium (default: 11)

  • Delta_z_reion float, optional

    Reionization transition width for hydorgen and neutral helium (default: 0.5)

  • z_reion_He float, optional

    Reionization redshift of singly-ionized helium (default: 3.5)

  • Delta_z_reion_He float, optional

    Reionization transition width for singly-ionized helium (default: 0.5)

  • exp_reion float, optional

    Power of 1+z appearing in tanh argument during reionization (default: 3/2)

  • rtol float, optional

    Relative tolerance for ODE solver (default: 1e-6)

  • atol float, optional

    Absolute tolerance for ODE solver (default: 1e-9)

  • solver diffrax.Solver, optional

    ODE solver instance (default: Kvaerno3())

  • max_steps int, optional

    Maximum solver steps (default: 1024)

Returns:

tuple

(xe_full_reion, lna_full, Tm, lna_Tm) - complete ionization history with reionization, log scale factor, matter temperature, and temperature grid