Cosmological Model

EMRI event-rate cosmological model and H₀ evaluation orchestration.

This module defines the cosmological event-rate model used for simulation: Model1CrossCheck samples EMRI events from a cosmological rate model, while LamCDMScenario and DarkEnergyScenario define the parameter spaces.

The Hubble-constant posterior evaluation has been extracted to master_thesis_code.bayesian_inference.bayesian_statistics (BayesianStatistics) and is invoked via main.py:evaluate() / --evaluate.

class master_thesis_code.cosmological_model.CosmologicalParameter(symbol: str, unit: str, lower_limit: float, upper_limit: float, value: float = 0.0, derivative_epsilon: float = 1e-06, is_fixed: bool = False, randomize_by_distribution: collections.abc.Callable[[float, float, numpy.random._generator.Generator], float] = <function uniform at 0x7f90ca6916c0>, fiducial_value: float = 1.0)[source]

Bases: Parameter

Parameters:

  • symbol (str)

  • unit (str)

  • lower_limit (float)

  • upper_limit (float)

  • value (float)

  • derivative_epsilon (float)

  • is_fixed (bool)

  • randomize_by_distribution (Callable[[float, float, Generator], float])

  • fiducial_value (float)

fiducial_value: float = 1.0
class master_thesis_code.cosmological_model.DarkEnergyScenario[source]

Bases: object

Omega_DE: float = 0.7274
Omega_m: float = 0.2726
h: float = 0.73
w_0: CosmologicalParameter
w_a: CosmologicalParameter
class master_thesis_code.cosmological_model.LamCDMScenario[source]

Bases: object

https://arxiv.org/pdf/2102.01708.pdf

Omega_m: CosmologicalParameter
h: CosmologicalParameter
w_0: float = -1.0
w_a: float = 0.0
master_thesis_code.cosmological_model.MBH_spin_distribution(lower_limit, upper_limit)[source]

https://iopscience.iop.org/article/10.1088/0004-637X/762/2/68/pdf

Parameters:

  • lower_limit (float)

  • upper_limit (float)

Return type:

float

class master_thesis_code.cosmological_model.Model1CrossCheck(rng=None, max_redshift_override=None)[source]

Bases: object

cross check of Model M1 in PHYSICAL REVIEW D 95, 103012 (2017)

Parameters:

  • rng (Generator | None)

  • max_redshift_override (float | None)

static R_emri(M)[source]
Parameters:

M (float)

Return type:

float

static dN_dz_of_mass(mass, redshift)[source]
Parameters:

  • mass (float)

  • redshift (float)

Return type:

float

emri_distribution(M, redshift)[source]
Parameters:

  • M (float)

  • redshift (float)

Return type:

float

emri_rate: int = 294
parameter_space: ParameterSpace
sample_emri_events(number_of_samples)[source]
Parameters:

number_of_samples (int)

Return type:

list[ParameterSample]

setup_emri_events_sampler()[source]
Return type:

None

setup_simplified_event_sampler()[source]
Return type:

None

simplified_event_mass_dependency(mass)[source]
Parameters:

mass (float)

Return type:

float

snr_threshold: float = 20
master_thesis_code.cosmological_model.gaussian(x, mu, sigma, a)[source]
Parameters:

  • x (float | ndarray[tuple[Any, ...], dtype[float64]])

  • mu (float)

  • sigma (float)

  • a (float)

Return type:

float | ndarray[tuple[Any, …], dtype[float64]]

master_thesis_code.cosmological_model.polynomial(x, a, b, c, d, e, f, g, h, i)[source]
Parameters:

  • x (float | ndarray[tuple[Any, ...], dtype[float64]])

  • a (float)

  • b (float)

  • c (float)

  • d (float)

  • e (float)

  • f (float)

  • g (float)

  • h (float)

  • i (float)

Return type:

float