primpy.efolds package

Calculations for the primordial Universe w.r.t. e-folds N.

primpy.efolds.inflation module

Differential equations for inflation w.r.t. e-folds N.

class primpy.efolds.inflation.InflationEquationsN(K, potential, track_time=False, track_eta=False, verbose=False)[source]

Bases: InflationEquations

Background equations during inflation w.r.t. e-folds N.

Solves background variables with e-folds N of the scale factor as independent variable for curved and flat universes using the Klein-Gordon and Friedmann equations.

Independent variable:

_N: e-folds of the scale-factor (the underscore here means that this is the as of yet uncalibrated scale factor)

Dependent variables:

phi: inflaton field
dphidN: d(phi)/dN
t: time (optional)
eta: conformal time (optional)

Methods

`H2`(x, y)	Compute the square of the Hubble parameter using the Friedmann equation.
`__call__`(x, y)	System of coupled ODEs for underlying variables.
`inflating`(x, y)	Inflation diagnostic for event tracking.
`sol`(sol, **kwargs)	Post-processing of `scipy.integrate.solve_ivp()` solution.
`w`(x, y)	Compute the equation of state parameter.

__call__(x, y)[source]: System of coupled ODEs for underlying variables.

H2(x, y)[source]: Compute the square of the Hubble parameter using the Friedmann equation.

w(x, y)[source]: Compute the equation of state parameter.

inflating(x, y)[source]: Inflation diagnostic for event tracking.

sol(sol, **kwargs)[source]: Post-processing of scipy.integrate.solve_ivp() solution.

primpy.efolds.perturbations module

Curvature perturbations with respect to e-folds _N.

class primpy.efolds.perturbations.PerturbationN(background, k, **kwargs)[source]

Bases: Perturbation

Curvature perturbation for wavenumber k with respect to e-folds _N=ln(a).

Solves the Mukhanov–Sasaki equations w.r.t. number of e-folds _N of the scale factor a for curved universes.

Parameters:

backgroundBunch object same as returned by scipy.integrate.solve_ivp(): Background solution as returned by primpy.efolds.inflation.InflationEquationsN.sol().
kfloat: wavenumber

class primpy.efolds.perturbations.ScalarModeN(background, k, **kwargs)[source]

Bases: ScalarMode

Template for scalar modes.

Methods

`__call__`(x, y)	Vector of derivatives.
`get_vacuum_ic_RST`()	Get initial conditions for scalar modes for RST vacuum w.r.t.
`mukhanov_sasaki_frequency_damping`()	Frequency and damping term of the Mukhanov-Sasaki equations for scalar modes.

__call__(x, y)[source]: Vector of derivatives.

mukhanov_sasaki_frequency_damping()[source]

Frequency and damping term of the Mukhanov-Sasaki equations for scalar modes.

Frequency and damping term of the Mukhanov-Sasaki equations for the comoving curvature perturbations R w.r.t. e-folds _N, where the e.o.m. is written as ddR + 2 * damping * dR + frequency**2 * R = 0.

get_vacuum_ic_RST()[source]: Get initial conditions for scalar modes for RST vacuum w.r.t. e-folds _N.

class primpy.efolds.perturbations.TensorModeN(background, k, **kwargs)[source]

Bases: TensorMode

Template for tensor modes.

Methods

`__call__`(x, y)	Vector of derivatives.
`get_vacuum_ic_RST`()	Get initial conditions for tensor modes for RST vacuum w.r.t.
`mukhanov_sasaki_frequency_damping`()	Frequency and damping term of the Mukhanov-Sasaki equations for tensor modes.

__call__(x, y)[source]: Vector of derivatives.

mukhanov_sasaki_frequency_damping()[source]

Frequency and damping term of the Mukhanov-Sasaki equations for tensor modes.

Frequency and damping term of the Mukhanov-Sasaki equations for the tensor perturbations h w.r.t. e-folds _N, where the e.o.m. is written as ddh + 2 * damping * dh + frequency**2 * h = 0.

get_vacuum_ic_RST()[source]: Get initial conditions for tensor modes for RST vacuum w.r.t. e-folds _N.