Source author record

Cornelius Rampf

Cornelius Rampf appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2026arXiv

Complex-time singular structure of the 1D Hou-Luo model

Starting from smooth initial data, we investigate the complex-time analytic structure of the one-dimensional Hou--Luo (HL) model, a wall approximation of the three-dimensional axisymmetric Euler equations. While the finite-time blow-up in this setting has been already established, here we chart the entire singular landscape. This analysis is enabled by a novel formulation of the HL model in Lagrangian coordinates, in which the time-Taylor coefficients of the flow fields are evaluated symbolically to high truncation order. Our results are threefold. First, we show that the Lagrangian series for the vorticity converges within the complex-time disc of radius~$t_\star >0$ and is free from (early-time) resonances that impede the Eulerian formulation. Second, applying asymptotic analysis on the series, we recover both the blow-up time and the singularity exponent with high accuracy. This also enables a quantitative assessment of the Beale--Kato--Majda criterion, which we find correctly identifies the blow-up time, but washes out the local singularity exponent, as it relies on a spatial supremum. Third, and most importantly, we develop a Lagrangian singularity theory that predicts the eye-shaped singularity profile observed in Eulerian coordinates by exploiting the driving mechanism of the blow-up: The accumulation of multiple fluid particles at the same Eulerian position. The employed techniques extend recently introduced methods for the inviscid Burgers equation [C. Rampf et al., Phys. Rev. Fluids 7 (2022) 10, 104610], and can be further adapted to higher spatial dimensions or other hydrodynamical equations.

preprint2023arXiv

Renormalization group and UV completion of cosmological perturbations: Gravitational collapse as a critical phenomenon

Cosmological perturbation theory is known to converge poorly for predicting the spherical collapse and void evolution of collisionless matter. Using the exact parametric solution as a testing ground, we develop two asymptotic methods in spherical symmetry that resolve the gravitational evolution to much higher accuracy than Lagrangian perturbation theory (LPT), which is the current gold standard in the literature. One of the methods selects a stable fixed-point solution of the renormalization-group flow equation, thereby predicting already at the leading order the critical exponent of the phase transition to collapsed structures. The other method completes the truncated LPT series far into the UV regime, by adding a non-analytic term that captures the critical nature of the gravitational collapse. We find that the UV method most accurately resolves the evolution of the nonlinear density as well as its one-point probability distribution function. Similarly accurate predictions are achieved with the renormalization-group method, especially when paired with Padé approximants. Further, our results yield new, very accurate, formulas to relate linear and nonlinear density contrasts. Finally, we chart possible ways on how to adapt our methods to the case of cosmological random field initial conditions.

preprint2022arXiv

Analytical growth functions for cosmic structures in a $Λ$CDM Universe

The cosmological fluid equations describe the early gravitational dynamics of cold dark matter (CDM), exposed to a uniform component of dark energy, the cosmological constant $Λ$. Perturbative predictions for the fluid equations typically assume that the impact of $Λ$ on CDM can be encapsulated by a refined growing mode $D$ of linear density fluctuations. Here we solve, to arbitrary high perturbative orders, the nonlinear fluid equations with an {\it Ansatz} for the fluid variables in increasing powers of $D$. We show that $Λ$ begins to populate the solutions starting at the fifth order in this strict $D$-expansion. By applying suitable resummation techniques, we recast these solutions to a standard perturbative series where not $D$, but essentially the initial gravitational potential serves as the bookkeeping parameter within the expansion. Then, by using the refined growth functions at second and third order in standard perturbation theory, we determine the matter power spectrum to one-loop accuracy as well as the leading-order contribution to the matter bispectrum. We find that employing our refined growth functions impacts the total power- and bispectra at a precision that is below one percent at late times. However, for the power spectrum, we find a characteristic scale-dependent suppression that is fairly similar to what is observed in massive neutrino cosmologies. Therefore, we recommend employing our refined growth functions in order to reduce theoretical uncertainties for analysing data in related pipelines.

preprint2016arXiv

Cosmological $N$-body simulations including radiation perturbations

Cosmological $N$-body simulations are the standard tool to study the emergence of the observed large-scale structure of the Universe. Such simulations usually solve for the gravitational dynamics of matter within the Newtonian approximation, thus discarding general relativistic effects such as the coupling between matter and radiation ($\equiv$ photons and neutrinos). In this paper we investigate novel hybrid simulations which incorporate interactions between radiation and matter to the leading order in General Relativity, whilst evolving the matter dynamics in full non-linearity according to Newtonian theory. Our hybrid simulations come with a relativistic space-time and make it possible to investigate structure formation in a unified framework. In the present work we focus on simulations initialized at $z=99$, and show that the extracted matter power spectrum receives up to $3\%$ corrections on very large scales through radiation. Our numerical findings compare favourably with linear analytical results from Fidler et al. (2016), from which we deduce that there cannot be any significant non-linear mode-coupling induced through linear radiation corrections.

preprint2016arXiv

Effects of Critical Collapse on Primordial Black-Hole Mass Spectra

Certain inflationary models as well as realisations of phase transitions in the early Universe predict the formation of primordial black holes. For most mass ranges, the fraction of matter in the form of primordial black holes is limited by many different observations on various scales. Primordial black holes are assumed to be formed when overdensities that cross the horizon have Schwarzschild radii larger than the horizon. Traditionally it was therefore assumed that primordial black-hole masses were equal to the horizon mass at their time of formation. However, detailed calculations of their collapse show that primordial black holes formed at each point in time should rather form a spectrum of different masses, obeying critical scaling. Though this has been known for more than fifteen years, the effect of this scaling behaviour is largely ignored when considering predictions for primordial black hole mass spectra. In this paper we consider the critical collapse scaling for a variety of models which produce primordial black holes, and find that it generally leads to a shift, broadening and an overall decrease of the mass contained in primordial black holes. This effect is model and parameter dependent and cannot be contained by a constant rescaling of the spectrum; it can become important and should be taken into account when comparing to observational constraints.

preprint2016arXiv

Lagrangian theory for cosmic structure formation with vorticity: Newtonian and post-Friedmann approximations

We study the nonlinear gravitational dynamics of a universe filled with a pressureless fluid and a cosmological constant $Λ$ in the context of Newtonian gravity, and in the relativistic post-Friedmann approach proposed in paper I [I. Milillo et al., Phys. Rev. D 92, 023519 (2015).]. The post-Friedmann approximation scheme is based on the $1/c$ expansion of the space-time metric and the energy-momentum tensor, and includes nonlinear Newtonian cosmology. Here we establish the nonlinear post-Friedmann framework in the Lagrangian-coordinates approach for structure formation. For this we first identify a Lagrangian gauge which is suitable for incorporating nonzero vorticity. We analyze our results in two limits: at the leading order we recover the fully nonlinear Newtonian cosmological equations in the Lagrangian formulation, and we provide a space-time metric consistent from the perspective of general relativity. We then linearize our expressions and recover the relativistic results at first order in cosmological perturbation theory. Therefore, the introduced approximation scheme provides a unified treatment for the two leading-order regimes, from the small scales described by Newtonian gravity to the large linear scale, where first-order relativistic cosmological perturbation theory gives a very good description of structure formation.

preprint2016arXiv

Relativistic Interpretation of Newtonian Simulations for Cosmic Structure Formation

The standard numerical tools for studying non-linear collapse of matter are Newtonian $N$-body simulations. Previous work has shown that these simulations are in accordance with General Relativity (GR) up to first order in perturbation theory, provided that the effects from radiation can be neglected. In this paper we show that the present day matter density receives more than 1$\%$ corrections from radiation on large scales if Newtonian simulations are initialised before $z=50$. We provide a relativistic framework in which \emph{unmodified} Newtonian simulations are compatible with linear GR even in the presence of radiation. Our idea is to use GR perturbation theory to keep track of the evolution of relativistic species and the relativistic space-time consistent with the Newtonian trajectories computed in $N$-body simulations. If metric potentials are sufficiently small, they can be computed using a first-order Einstein--Boltzmann code such as CLASS. We make this idea rigorous by defining a class of GR gauges, the \emph{Newtonian motion} gauges, which are defined such that matter particles follow Newtonian trajectories. We construct a simple example of a relativistic space-time within which unmodified Newtonian simulations can be interpreted.

preprint2015arXiv

Boltzmann hierarchy for interacting neutrinos I: formalism

Starting from the collisional Boltzmann equation, we derive for the first time and from first principles the Boltzmann hierarchy for neutrinos including interactions with a scalar particle. Such interactions appear, for example, in majoron-like models of neutrino mass generation. We study two limits of the scalar mass: (i) An extremely massive scalar whose only role is to mediate an effective 4-fermion neutrino-neutrino interaction, and (ii) a massless scalar that can be produced in abundance and thus demands its own Boltzmann hierarchy. In contrast to, e.g., the first-order Boltzmann hierarchy for Thomson-scattering photons, our interacting neutrino/scalar Boltzmann hierarchies contain additional momentum-dependent collision terms arising from a non-negligible energy transfer in the neutrino-neutrino and neutrino-scalar interactions. This necessitates that we track each momentum mode of the phase space distributions individually, even if the particles were massless. Comparing our hierarchy with the commonly used $(c_{\rm eff}^2,c_{\rm vis}^2)$-parameterisation, we find no formal correspondence between the two approaches, which raises the question of whether the latter parameterisation even has an interpretation in terms of particle scattering. Lastly, although we have invoked majoron-like models as a motivation for our study, our treatment is in fact generally applicable to all scenarios in which the neutrino and/or other ultrarelativistic fermions interact with scalar particles.

preprint2015arXiv

General relativistic corrections to $N$-body simulations and the Zel'dovich approximation

The initial conditions for Newtonian $N$-body simulations are usually generated by applying the Zel'dovich approximation to the initial displacements of the particles using an initial power spectrum of density fluctuations generated by an Einstein-Boltzmann solver. We show that in most gauges the initial displacements generated in this way receive a first-order relativistic correction. We define a new gauge, the $N$-body gauge, in which this relativistic correction vanishes and show that a conventional Newtonian $N$-body simulation includes all first-order relativistic contributions (in the absence of radiation) if we identify the coordinates in Newtonian simulations with those in the relativistic $N$-body gauge.

preprint2015arXiv

How smooth are particle trajectories in a $Λ$CDM Universe?

It is shown here that in a flat, cold dark matter (CDM) dominated Universe with positive cosmological constant ($Λ$), modelled in terms of a Newtonian and collisionless fluid, particle trajectories are analytical in time (representable by a convergent Taylor series) until at least a finite time after decoupling. The time variable used for this statement is the cosmic scale factor, i.e., the "$a$-time", and not the cosmic time. For this, a Lagrangian-coordinates formulation of the Euler-Poisson equations is employed, originally used by Cauchy for 3-D incompressible flow. Temporal analyticity for $Λ$CDM is found to be a consequence of novel explicit all-order recursion relations for the $a$-time Taylor coefficients of the Lagrangian displacement field, from which we derive the convergence of the $a$-time Taylor series. A lower bound for the $a$-time where analyticity is guaranteed and shell-crossing is ruled out is obtained, whose value depends only on $Λ$ and on the initial spatial smoothness of the density field. The largest time interval is achieved when $Λ$ vanishes, i.e., for an Einstein-de Sitter universe. Analyticity holds also if, instead of the $a$-time, one uses the linear structure growth $D$-time, but no simple recursion relations are then obtained. The analyticity result also holds when a curvature term is included in the Friedmann equation for the background, but inclusion of a radiation term arising from the primordial era spoils analyticity.

preprint2015arXiv

Second-order cosmological perturbation theory and initial conditions for $N$-body simulations

We use gauge-invariant cosmological perturbation theory to calculate the displacement field that sets the initial conditions for $N$-body simulations. Using first and second-order fully relativistic perturbation theory in the synchronous-comoving gauge, allows us to go beyond the Newtonian predictions and to calculate relativistic corrections to it. We use an Einstein--de Sitter model, including both growing and decaying modes in our solutions. The impact of our results should be assessed through the implementation of the featured displacement in cosmological $N$-body simulations.

preprint2014arXiv

A Relativistic view on large scale N-body simulations

We discuss the relation between the output of Newtonian N-body simulations on scales that approach or exceed the particle horizon to the description of General Relativity. At leading order, the Zeldovich approximation is correct on large scales, coinciding with the General Relativistic result. At second order in the initial metric potential, the trajectories of particles deviate from the second order Newtonian result and hence the validity of 2LPT initial conditions should be reassessed when used in very large simulations. We also advocate using the expression for the synchronous gauge density as a well behaved measure of density fluctuations on such scales.

preprint2014arXiv

Astrophysical Bose-Einstein Condensates and Superradiance

We investigate gravitational analogue models to describe slowly rotating objects (e.g., dark-matter halos, or boson stars) in terms of Bose-Einstein condensates, trapped in their own gravitational potentials. We begin with a modified Gross-Pitaevskii equation, and show that the resulting background equations of motion are stable, as long as the rotational component is treated as a small perturbation. The dynamics of the fluctuations of the velocity potential are effectively governed by the Klein-Gordon equation of a "Eulerian metric," where we derive the latter by the use of a relativistic Lagrangian extrapolation. Superradiant scattering on such objects is studied. We derive conditions for its occurence and estimate its strength. Our investigations might give an observational handle to phenomenologically constrain Bose-Einstein condensates.

preprint2014arXiv

Frame dragging and Eulerian frames in General Relativity

The physical interpretation of cold dark matter perturbations is clarified by associating Bertschinger's Poisson gauge with a Eulerian/observer's frame of reference. We obtain such an association by using a Lagrangian approach to relativistic cosmological structure formation. Explicitly, we begin with the second-order solution of the Einstein equations in a synchronous/comoving coordinate system---which defines the Lagrangian frame, and transform it to a Poissonian coordinate system. The generating vector of this coordinate/gauge transformation is found to be the relativistic displacement field. The metric perturbations in the Poissonian coordinate system contain known results from standard/Eulerian Newtonian perturbation theory, but contain also purely relativistic corrections. On sub-horizon scales these relativistic corrections are dominated by the Newtonian bulk part. These corrections however set up non-linear constraints for the density and for the velocity which become important on scales close to the horizon. Furthermore, we report the occurence of a transverse component in the displacement field, and find that it induces a non-linear frame dragging as seen in the observer's frame, which is sub-dominant at late-times and sub-horizon scales. Finally, we find two other gauges which can be associated with a Eulerian frame. We argue that the Poisson gauge is to be preferred because it comes with the simplest physical interpretation.

preprint2014arXiv

Relativistic Lagrangian displacement field and tensor perturbations

We investigate the purely spatial Lagrangian coordinate transformation from the Lagrangian to the basic Eulerian frame. We demonstrate three techniques for extracting the relativistic displacement field from a given solution in the Lagrangian frame. These techniques are (a) from defining a local set of Eulerian coordinates embedded into the Lagrangian frame; (b) from performing a specific gauge transformation; and (c) from a fully non-perturbative approach based on the ADM split. The latter approach shows that this decomposition is not tied to a specific perturbative formulation for the solution of the Einstein equations. Rather, it can be defined at the level of the non-perturbative coordinate change from the Lagrangian to the Eulerian description. Studying such different techniques is useful because it allows us to compare and develop further the various approximation techniques available in the Lagrangian formulation. We find that one has to solve the gravitational wave equation in the relativistic analysis, otherwise the corresponding Newtonian limit will necessarily contain spurious non-propagating tensor artefacts at second order in the Eulerian frame. We also derive the magnetic part of the Weyl tensor in the Lagrangian frame, and find that it is not only excited by gravitational waves but also by tensor perturbations which are induced through the non-linear frame-dragging. We apply our findings to calculate for the first time the relativistic displacement field, up to second order, for a $Λ$CDM Universe in the presence of a local primordial non-Gaussian component. Finally, we also comment on recent claims about whether mass conservation in the Lagrangian frame is violated.

preprint2014arXiv

The recursion relation in Lagrangian perturbation theory

We derive a recursion relation in the framework of Lagrangian perturbation theory, appropriate for studying the inhomogeneities of the large scale structure of the universe. We use the fact that the perturbative expansion of the matter density contrast is in one-to-one correspondence with standard perturbation theory (SPT) at any order. This correspondence has been recently shown to be valid up to fourth order for a non-relativistic, irrotational and dust-like component. Assuming it to be valid at arbitrary (higher) order, we express the Lagrangian displacement field in terms of the perturbative kernels of SPT, which are itself given by their own and well-known recursion relation. We argue that the Lagrangian solution always contains more non-linear information in comparison with the SPT solution, (mainly) if the non-perturbative density contrast is restored after the displacement field is obtained.

preprint2013arXiv

Initial conditions for cold dark matter particles and General Relativity

We describe the irrotational dust component of the universe in terms of a relativistic gradient expansion and transform the resulting synchronous metric to a Newtonian coordinate system. The two metrics are connected via a space-like displacement field and a time-like perturbation, providing a relativistic generalization of the transformation from Lagrangian to Eulerian coordinates. The relativistic part of the displacement field generates already at initial time a non-local density perturbation at second order. This is a purely relativistic effect since it originates from space-time mixing. We give two options, the passive and the active approach, on how to include the relativistic corrections for example in N-body simulations. In the passive approach we treat the corrections as a non-Gaussian modification of the initial Gaussian field (primordial non-Gaussianity could be incorporated as well). The induced non-Gaussianity depends on scale and the redshift at which initial conditions are set, with f_NL ~ few for small enough scales and redshifts. In the active approach we show how to use the relativistic trajectory to obtain the initial displacement and velocity of particles for N-body simulations without modifying the initial Gaussian field.

preprint2012arXiv

Lagrangian perturbations and the matter bispectrum I: fourth-order model for non-linear clustering

We investigate the Lagrangian perturbation theory of a homogeneous and isotropic universe in the non-relativistic limit, and derive the solutions up to the fourth order. These solutions are needed for example for the next-to-leading order correction of the (resummed) Lagrangian matter bispectrum, which we study in an accompanying paper. We focus on flat cosmologies with a vanishing cosmological constant, and provide an in-depth description of two complementary approaches used in the current literature. Both approaches are solved with two different sets of initial conditions---both appropriate for modelling the large-scale structure. Afterwards we consider only the fastest growing mode solution, which is not affected by either of these choices of initial conditions. Under the reasonable approximation that the linear density contrast is evaluated at the initial Lagrangian position of the fluid particle, we obtain the nth-order displacement field in the so-called initial position limit: the nth order displacement field consists of 3(n-1) integrals over n linear density contrasts, and obeys self-similarity. Then, we find exact relations between the series in Lagrangian and Eulerian perturbation theory, leading to identical predictions for the density contrast and the peculiar-velocity divergence up to the fourth order.

preprint2012arXiv

Lagrangian perturbations and the matter bispectrum II: the resummed one-loop correction to the matter bispectrum

This is part two in a series of papers in which we investigate an approach based on Lagrangian perturbation theory (LPT) to study the non-linear evolution of the large-scale structure distribution in the universe. Firstly, we compute the matter bispectrum in real space using LPT up one-loop order, for both Gaussian and non-Gaussian initial conditions. In the initial position limit, we find that the one-loop bispectrum computed in this manner is identical to its counterpart obtained from standard Eulerian perturbation theory (SPT). Furthermore, the LPT formalism allows for a simple reorganisation of the perturbative series corresponding to the resummation of an infinite series of perturbations in SPT. Applying this method, we find a resummed one-loop bispectrum that compares favourably with results from N-body simulations. We generalise the resummation method also to the computation of the redshift-space bispectrum up to one loop.

preprint2012arXiv

Zel'dovich approximation and General Relativity

We show how the Zel'dovich approximation and the second order displacement field of Lagrangian perturbation theory can be obtained from a general relativistic gradient expansion in ΛCDM cosmology. The displacement field arises as a result of a second order non-local coordinate transformation which brings the synchronous/comoving metric into a Newtonian form. We find that, with a small modification, the Zel'dovich approximation holds even on scales comparable to the horizon. The corresponding density perturbation is not related to the Newtonian potential via the usual Poisson equation but via a modified Helmholtz equation. This is a consequence of causality not present in the Newtonian theory. The second order displacement field receives relativistic corrections that are subdominant on short scales but are comparable to the second order Newtonian result on scales approaching the horizon. The corrections are easy to include when setting up initial conditions in large N-body simulations.

preprint2010arXiv

Cosmological parameters from large scale structure - geometric versus shape information

The matter power spectrum as derived from large scale structure (LSS) surveys contains two important and distinct pieces of information: an overall smooth shape and the imprint of baryon acoustic oscillations (BAO). We investigate the separate impact of these two types of information on cosmological parameter estimation, and show that for the simplest cosmological models, the broad-band shape information currently contained in the SDSS DR7 halo power spectrum (HPS) is by far superseded by geometric information derived from the baryonic features. An immediate corollary is that contrary to popular beliefs, the upper limit on the neutrino mass m_νpresently derived from LSS combined with cosmic microwave background (CMB) data does not in fact arise from the possible small-scale power suppression due to neutrino free-streaming, if we limit the model framework to minimal LambdaCDM+m_ν. However, in more complicated models, such as those extended with extra light degrees of freedom and a dark energy equation of state parameter w differing from -1, shape information becomes crucial for the resolution of parameter degeneracies. This conclusion will remain true even when data from the Planck surveyor become available. In the course of our analysis, we introduce a new dewiggling procedure that allows us to extend consistently the use of the SDSS HPS to models with an arbitrary sound horizon at decoupling. All the cases considered here are compatible with the conservative 95%-bounds \sum m_ν< 1.16 eV, N_eff = 4.8 \pm 2.0.