Top |
GEnum ╰── NcHIPertGWVars GInterface ╰── NcHIPertIGW GObject ╰── NcmHOAA ╰── NcHIPertGW
This object provides the computation of the gwatic mode for the cosmological perturbations. It solves the equation of motion for the gauge invariant variable (see Vitenti (2013) for notation and details) $$\zeta \equiv \Psi - \frac{2\bar{K}}{\kappa(\bar{\rho} + \bar{p})} + H\mathcal{V}.$$ Its conjugated momentum is give by \begin{split} P_\zeta &= \frac{2\bar{D}^2_\bar{K}\Psi}{x^3H}, \end{split}
The equations of motion in their first order form are
\begin{align}
\zeta^\prime &= \frac{P_\zeta}{m_\zeta}, \\
P_\zeta^\prime &= -m_\zeta\mu_\zeta^2\zeta.
\end{align}
The mass $m_\zeta$ and the frequency $\mu_\zeta$ are defined by
\begin{align}
m_\zeta &= \frac{3\Delta_\bar{K}(\bar{\rho} + \bar{p})}{\rho_\text{crit0} N x^3 c_s^2 E^2}, \\
\mu_\zeta^2 &= x^2N^2c_s^2k^2,
\end{align}
where $\bar{\rho} + \bar{p}$ is the background total energy density plus pressure,
$E^2 = H^2/H_0^2$ is the dimensionless Hubble function squared (nc_hicosmo_E2()
), $c_s^2$ the speed of sound,
$N$ is the lapse function that in this case (using $\alpha$ as time variable) is $N \equiv \vert{}E\vert^{-1}$, $\rho_\text{crit0}$
is the critical density today defined by $\rho_\text{crit0} \equiv 3H_0^2/\kappa$ and $$\Delta_\bar{K} \equiv \frac{k^2}{k^2 + \Omega_{k0}}.$$
gdouble nc_hipert_igw_eval_mnu (NcHIPertIGW *igw
,const gdouble tau
,const gdouble k
);
FIXME
gdouble nc_hipert_igw_eval_nu (NcHIPertIGW *igw
,const gdouble tau
,const gdouble k
);
FIXME
gdouble nc_hipert_igw_eval_dlnmnu (NcHIPertIGW *igw
,const gdouble tau
,const gdouble k
);
FIXME
void nc_hipert_igw_eval_system (NcHIPertIGW *igw
,const gdouble tau
,const gdouble k
,gdouble *nu
,gdouble *dlnmnu
);
FIXME
void nc_hipert_igw_get_sing_info (NcHIPertIGW *igw
,const gdouble k
,const guint sing
,gdouble *ts
,gdouble *dts_i
,gdouble *dts_f
,NcmHOAASingType *st
);
gdouble nc_hipert_igw_eval_sing_mnu (NcHIPertIGW *igw
,const gdouble tau_m_taus
,const gdouble k
,const guint sing
);
FIXME
gdouble nc_hipert_igw_eval_sing_dlnmnu (NcHIPertIGW *igw
,const gdouble tau_m_taus
,const gdouble k
,const guint sing
);
FIXME
void nc_hipert_igw_eval_sing_system (NcHIPertIGW *igw
,const gdouble tau_m_taus
,const gdouble k
,const guint sing
,gdouble *nu
,gdouble *dlnmnu
);
FIXME
NcHIPertGW *
nc_hipert_gw_ref (NcHIPertGW *pa
);
Increases the reference count of pa
.
void
nc_hipert_gw_clear (NcHIPertGW **pa
);
Decreases the reference count of *pa
and sets *pa
to NULL.