Spherical Collapse Solvers

Class providing numerical solvers for the collapse of spherical density perturbations in a given cosmology—computing the critical linear overdensity \(\delta_\mathrm{c}\), the virial density contrast \(\Delta_\mathrm{vir}\), and the turnaround-to-virial radius ratio as functions of cosmic time by integrating the spherical collapse equations. These tabulated quantities are consumed by criticalOverdensityClass and virialDensityContrastClass implementations, allowing the halo collapse threshold and virial radius definition to vary with redshift and cosmological model.

Default implementation: sphericalCollapseSolverCllsnlssMttrCsmlgclCnstnt

Methods

criticalOverdensity

Returns a tabulation of the critical overdensity for collapse.

  • double precision time [in]

  • logical tableStore [in]

  • class(table1D) criticalOverdensity_ [inout]

virialDensityContrast

Returns a tabulation of the virial density contrast.

  • double precision time [in]

  • logical tableStore [in]

  • class(table1D) virialDensityContrast_ [inout]

radiusTurnaround

Returns a tabulation of the ratio of turnaround to virial radii.

  • double precision time [in]

  • logical tableStore [in]

  • class(table1D) radiusTurnaround_ [inout]

linearNonlinearMap

Returns a mapping of linear to nonlinear density overdensity.

  • double precision time [in]

  • class(table2DLinLinLin) linearNonlinearMap_ [inout]

sphericalCollapseSolverBaryonsDarkMatterDarkEnergy

A spherical collapse solver for universes consisting of baryons, collisionless dark matter, and dark energy, computing the evolution of a spherical top-hat overdensity through turnaround and virialization. Baryons can be treated as either clustering or non-clustering via [baryonsClustering], and the energy-fixing radius for dark energy is set by [radiusFixed].

Parameters

  • [baryonsCluster] (boolean) — If true baryons are assumed to cluster in the same way as collisionless matter. If false, baryons are assumed not to cluster at all.

  • [energyFixedAt] (string; one of undefined, turnaround, virialization; default turnaround) — The radius at which the energy of a spherical top-hat perturbation in a dark energy cosmology can be considered to be fixed.

  • [tablePointsPerOctave] (integer; default 300) — The number of points per octave of time at which to tabulate solutions.

  • [perturbationSmall] (real; default 1.0d-3) — The largest initial perturbation considered to be sufficiently small. Larger initial perturbations will trigger an error.

sphericalCollapseSolverCllsnlssMttrCsmlgclCnstnt

A spherical collapse solver for universes consisting of collisionless matter and a cosmological constant.

(Default implementation)

Methods

getTable

Get the requested table.

  • double precision time [in]

  • type(varying_string) fileName [in]

  • logical tableStore [in]

  • type(enumerationCllsnlssMttCsmlgclCnstntClcltnType) calculationType [in]

  • class(table1D) table_ [inout]

tabulate

Construct a tabulated solution.

  • double precision time [in]

  • class(table1D) sphericalCollapse_ [inout]

  • type(enumerationCllsnlssMttCsmlgclCnstntClcltnType) calculationType [in]

sphericalCollapseSolverCllsnlssMttrDarkEnergy

A spherical collapse solver for universes consisting of collisionless matter and dark energy.

Parameters

  • [energyFixedAt] (string; one of undefined, turnaround, virialization; default turnaround) — The radius at which the energy of a spherical top-hat perturbation in a dark energy cosmology can be considered to be fixed.