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 0x7ff7e2ec7600>, fiducial_value: float = 1.0)[source]

Bases: Parameter

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

Bases: object

Omega_DE: float = 0.75
Omega_m: float = 0.25
de_equation(z)[source]
Parameters:

z (float)

Return type:

float

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:
Return type:

float

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

Bases: object

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

Parameters:

rng (Generator | None)

static R_emri(M)[source]
Parameters:

M (float)

Return type:

float

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

float

emri_distribution(M, redshift)[source]
Parameters:
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:
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:
Return type:

float