The merging rate of cosmic structures is computed, relying on the Ansatz that they can be predicted in the initial linear density field from the coalescence of critical points with increasing smoothing scale, used here as a proxy for cosmic time. Beyond the mergers of peaks with saddle points (a proxy for halo mergers), we consider the coalescence and nucleation of all sets of critical points, including wall-saddle to filament-saddle and wall-saddle to minima (a proxy for filament and void mergers respectively), as they impact the geometry of galactic infall, and in particular filament disconnection. Analytical predictions of the one-point statistics are validated against multiscale measurements in 2D and 3D realisations of Gaussian random fields (the corresponding code being available upon request) and compared qualitatively to cosmological $N$-body simulations at early times ($z\geq 10$) and large scales ($\geq 5\, \mathrm{Mpc}/h$). The rate of filament coalescence is compared to the merger rate of haloes and the two-point clustering of these events is computed, along with their cross-correlations with critical points. These correlations are qualitatively consistent with the preservation of the connectivity of dark matter haloes, and the impact of the large scale structures on assembly bias. The destruction rate of haloes and voids as a function of mass and redshift is quantified down to $z=0$ for a $\Lambda$CDM cosmology. The one-point statistics in higher dimensions are also presented, together with consistency relations between critical point and critical event counts.

We consider the application of peaks theory to the calculation of the number density of peaks relevant for primordial black hole (PBH) formation. For PBHs, the final mass is related to the amplitude and scale of the perturbation from which it forms, where the scale is defined as the scale at which the compaction function peaks. We therefore extend peaks theory to calculate not only the abundance of peaks of a given amplitude, but peaks of a given amplitude and scale. A simple fitting formula is given in the high-peak limit relevant for PBH formation. We also adapt the calculation to use a Gaussian smoothing function, ensuring convergence regardless of the choice of power spectrum.

We investigate the impact of the number of filaments connected to the nodes of the cosmic web on the physical properties of their galaxies using the Sloan Digital Sky Survey. We compare these measurements to the cosmological hydrodynamical simulations Horizon-(no)AGN and Simba. We find that more massive galaxies are more connected, in qualitative agreement with theoretical predictions and measurements in dark matter only simulation. The star formation activity and morphology of observed galaxies both display some dependence on the connectivity of the cosmic web at fixed stellar mass: less star forming and less rotation supported galaxies also tend to have higher connectivity. These results qualitatively hold both for observed and virtual galaxies, and can be understood given that the cosmic web is the main source of fuel for galaxy growth. The simulations show the same trends at fixed halo mass, suggesting that the geometry of filamentary infall impacts galaxy properties beyond the depth of the local potential well. Based on simulations, it is also found that AGN feedback is key in reversing the relationship between stellar mass and connectivity at fixed halo mass. Technically, connectivity is a practical observational proxy for past and present accretion (minor mergers or diffuse infall).

The simplest models of dark matter halo formation are based on the heuristic assumption, motivated by spherical collapse, that virialized haloes originate from initial regions that are maxima of the smoothed density field. Here, we replace this notion with the dynamical requirement that protohalo patches be regions where the local gravitational flow converges to a point. For this purpose, we look for spheres whose gravitational acceleration at the boundary -- relative to their center of mass -- points towards their geometric center: that is, spheres with null dipole moment. We show that these configurations are minima of the total energy, i.e. the most energetically bound spheres. For this reason, we study peaks of the energy overdensity field, and argue that the approach shows considerable promise. This change simply requires that one modify the standard top-hat filter, with the added important benefit that, for power spectra of cosmological interest, the resulting model is no longer plagued by divergences. Although the formalism is no more complicated than the overdensity based approach, the model is richer in the sense that it naturally predicts scatter in the overdensities of protohalo patches that are destined to form haloes of the same mass, in qualitative agreement with simulations of halo formation.

We report on the confirmation and mass determination of a transiting planet orbiting the old and inactive G7 dwarf star HD219666 (Mstar = 0.92 +/- 0.03 MSun, Rstar = 1.03 +/- 0.03 RSun, tau_star = 10 +/- 2 Gyr). With a mass of Mb = 16.6 +/- 1.3 MEarth, a radius of Rb = 4.71 +/- 0.17 REarth, and an orbital period of P ~ 6 days, HD219666b is a new member of a rare class of exoplanets: the hot-Neptunes. The Transiting Exoplanet Survey Satellite (TESS) observed HD219666 (also known as TOI-118) in its Sector 1 and the light curve shows four transit-like events, equally spaced in time. We confirmed the planetary nature of the candidate by gathering precise radial-velocity measurements with HARPS@ESO3.6m. We used the co-added HARPS spectrum to derive the host star fundamental parameters (Teff = 5527 +/- 65 K, log g = 4.40 +/- 0.11 (cgs), [Fe/H]= 0.04 +/- 0.04 dex, log R'HK = -5.07 +/- 0.03), as well as the abundances of many volatile and refractory elements. The host star brightness (V = 9.9) makes it suitable for further characterisation by means of in-transit spectroscopy. The determination of the planet orbital obliquity, along with the atmospheric metal-to-hydrogen content and thermal structure could provide us with important clues on the formation mechanisms of this class of objects.

The strikingly anisotropic large-scale distribution of matter made of an extended network of voids delimited by sheets, themselves segmented by filaments, within which matter flows towards compact nodes where they intersect, imprints its geometry on the dynamics of cosmic flows, ultimately shaping the distribution of galaxies and the redshift evolution of their properties. The (filament-type) saddle points of this cosmic web provide a local frame in which to quantify the induced physical and morphological evolution of galaxies on large scales. The properties of virtual galaxies within the Horizon-AGN simulation are stacked in such a frame. The iso-contours of the galactic number density, mass, specific star formation rate (sSFR), kinematics and age are clearly aligned with the filament axis with steep gradients perpendicular to the filaments. A comparison to a simulation without feedback from active galactic nuclei (AGN) illustrates its impact on quenching star formation of centrals away from the saddles. The redshift evolution of the properties of galaxies and their age distribution are consistent with the geometry of the bulk flow within that frame. They compare well with expectations from constrained Gaussian random fields and the scaling with the mass of non-linearity, modulo the redshift dependent impact of feedback processes. Physical properties such as sSFR and kinematics seem not to depend only on mean halo mass and density: the residuals trace the geometry of the saddle, which could point to other environment-sensitive physical processes, such as spin advection, and AGN feedback at high mass.

The role of the cosmic web in shaping galaxy properties is investigated in the GAMA spectroscopic survey in the redshift range $0.03 \leq z \leq 0.25$. The stellar mass, $u - r$ dust corrected colour and specific star formation rate (sSFR) of galaxies are analysed as a function of their distances to the 3D cosmic web features, such as nodes, filaments and walls, as reconstructed by DisPerSE. Significant mass and type/colour gradients are found for the whole population, with more massive and/or passive galaxies being located closer to the filament and wall than their less massive and/or star-forming counterparts. Mass segregation persists among the star-forming population alone. The red fraction of galaxies increases when closing in on nodes, and on filaments regardless of the distance to nodes. Similarly, the star-forming population reddens (or lowers its sSFR) at fixed mass when closing in on filament, implying that some quenching takes place. Comparable trends are also found in the state-of-the-art hydrodynamical simulation Horizon-AGN. These results suggest that on top of stellar mass and large-scale density, the traceless component of the tides from the anisotropic large-scale environment also shapes galactic properties. An extension of excursion theory accounting for filamentary tides provides a qualitative explanation in terms of anisotropic assembly bias: at a given mass, the accretion rate varies with the orientation and distance to filaments. It also explains the absence of type/colour gradients in the data on smaller, non-linear scales.

The mass, accretion rate and formation time of dark matter haloes near proto-filaments (identified as saddle points of the potential) are analytically predicted using a conditional version of the excursion set approach in its so-called "upcrossing" approximation. The model predicts that at fixed mass, mass accretion rate and formation time vary with orientation and distance from the saddle, demonstrating that assembly bias is indeed influenced by the tides imposed by the cosmic web. Starved, early forming haloes of smaller mass lie preferentially along the main axis of filaments, while more massive and younger haloes are found closer to the nodes. Distinct gradients for distinct tracers such as typical mass and accretion rate occur because the saddle condition is anisotropic, and because the statistics of these observables depend on both the conditional means and their covariances. The theory is extended to other critical points of the potential field. The response of the mass function to variations of the matter density field (the so-called large scale bias) is computed, and its trend with accretion rate is shown to invert along the filament. The signature of this model should correspond at low redshift to an excess of reddened galactic hosts at fixed mass along preferred directions, as recently reported in spectroscopic and photometric surveys and in hydrodynamical simulations. The anisotropy of the cosmic web emerges therefore as a significant ingredient to describe jointly the dynamics and physics of galaxies, e.g. in the context of intrinsic alignments or morphological diversity.

We present precise measurements of the assembly bias of dark matter halos, i.e. the dependence of halo bias on other properties than the mass, using curved "separate universe" N-body simulations which effectively incorporate an infinite-wavelength matter overdensity into the background density. This method measures the LIMD (local-in-matter-density) bias parameters $b_n$ in the large-scale limit. We focus on the dependence of the first two Eulerian biases $b^E_1$ and $b^E_2$ on four halo properties: the concentration, spin, mass accretion rate, and ellipticity. We quantitatively compare our results with previous works in which assembly bias was measured on fairly small scales. Despite this difference, our findings are in good agreement with previous results. We also look at the joint dependence of bias on two halo properties in addition to the mass. Finally, using the excursion set peaks model, we attempt to shed new insights on how assembly bias arises in this analytical model.

The dark matter halos that host galaxies and clusters form out of initial high-density patches, providing a biased tracer of the linear matter density field. In the simplest local bias approximation, the halo field is treated as a perturbative series in the average overdensity of the Lagrangian patch. In more realistic models, however, additional quantities will affect the clustering of halo-patches, and this expansion becomes a function of several stochastic variables. In this paper, we present a general multivariate expansion scheme that can parametrize the clustering of any biased Lagrangian tracer, given only the variables involved and their symmetry (in our case rotational invariance). This approach is based on an expansion in the orthonormal polynomials associated with the relevant variables, so that no renormalization of the coefficients ever occurs. We provide explicit expression for the series coefficients, or Lagrangian bias parameters, in the case of peaks of the linear density field. As an application of our formalism, we present a simple derivation of the original BBKS formula, and compute the non-Gaussian bias in the presence of a primordial trispectrum of the local shape.

The excursion set approach uses the statistics of the density field smoothed on a wide range of scales, to gain insight into a number of interesting processes in nonlinear structure formation, such as cluster assembly, merging and clustering. The approach treats the curve defined by the height of the overdensity fluctuation field when changing the smoothing scale as a random walk. The steps of the walks are often assumed to be uncorrelated, so that the walk heights are a Markov process, even though this assumption is known to be inaccurate for physically relevant filters. We develop a model in which the walk steps, rather than heights, are a Markov process, and correlations between steps arise because of nearest neighbour interactions. This model is a particular case of a general class, which we call Markov Velocity models. We show how these can approximate the walks generated by arbitrary power spectra and filters, and, unlike walks with Markov heights, provide a very good approximation to physically relevant models. We define a Markov Velocity Monte Carlo algorithm to generate walks whose first crossing distribution is very similar to that of TopHat-smoothed LCDM walks. Finally, we demonstrate that Markov Velocity walks generically exhibit a simple but realistic form of assembly bias, so we expect them to be useful in the construction of more realistic merger history trees.

The simplest stochastic halo formation models assume that the traceless part of the shear field acts to increase the initial overdensity (or decrease the underdensity) that a protohalo (or protovoid) must have if it is to form by the present time. Equivalently, it is the difference between the overdensity and (the square root of the) shear that must be larger than a threshold value. To estimate the effect this has on halo abundances using the excursion set approach, we must solve for the first crossing distribution of a barrier of constant height by the random walks associated with the difference, which is now (even for Gaussian initial conditions) a non-Gaussian variate. The correlation properties of such non-Gaussian walks are inherited from those of the density and the shear, and, since they are independent processes, the solution is in fact remarkably simple. We show that this provides an easy way to understand why earlier heuristic arguments about the nature of the solution worked so well. In addition to modelling halos and voids, this potentially simplifies models of the abundance and spatial distribution of filaments and sheets in the cosmic web.

Insight into a number of interesting questions in cosmology can be obtained from the first crossing distributions of physically motivated barriers by random walks with correlated steps. We write the first crossing distribution as a formal series, ordered by the number of times a walk upcrosses the barrier. Since the fraction of walks with many upcrossings is negligible if the walk has not taken many steps, the leading order term in this series is the most relevant for understanding the massive objects of most interest in cosmology. This first term only requires knowledge of the bivariate distribution of the walk height and slope, and provides an excellent approximation to the first crossing distribution for all barriers and smoothing filters of current interest. We show that this simplicity survives when extending the approach to the case of non-Gaussian random fields. For non-Gaussian fields which are obtained by deterministic transformations of a Gaussian, the first crossing distribution is simply related to that for Gaussian walks crossing a suitably rescaled barrier. Our analysis shows that this is a useful way to think of the generic case as well. Although our study is motivated by the possibility that the primordial fluctuation field was non-Gaussian, our results are general. In particular, they do not assume the non-Gaussianity is small, so they may be viewed as the solution to an excursion set analysis of the late-time, nonlinear fluctuation field rather than the initial one. They are also useful for models in which the barrier height is determined by quantities other than the initial density, since most other physically motivated variables (such as the shear) are usually stochastic and non-Gaussian. We use the Lognormal transformation to illustrate some of our arguments.

In this paper we provide an alternative method to compute correlation functions in the in-in formalism, with a modified set of Feynman rules to compute loop corrections. The diagrammatic expansion is based on an iterative solution of the equation of motion for the quantum operators with only retarded propagators, which makes each diagram intrinsically local (whereas in the standard case locality is the result of several cancellations) and endowed with a straightforward physical interpretation. While the final result is strictly equivalent, as a bonus the formulation presented here also contains less graphs than other diagrammatic approaches to in-in correlation functions. Our method is particularly suitable for applications to cosmology.

We prove an exact relation between the tensor and the scalar primordial power spectra generated during inflation. Such a mapping considerably simplifies the derivation of any power spectra as they can be obtained from the study of the tensor modes only, which are much easier to solve. As an illustration, starting from the second order slow-roll tensor power spectrum, we derive in a few lines the next-to-next-to-leading order power spectrum of the comoving curvature perturbation in generalized single field inflation with a varying speed of sound.

Recently, we provided a simple but accurate formula which closely approximates the first crossing distribution associated with random walks having correlated steps. The approximation is accurate for the wide range of barrier shapes of current interest and is based on the requirement that, in addition to having the right height, the walk must cross the barrier going upwards. Therefore, it only requires knowledge of the bivariate distribution of the walk height and slope, and is particularly useful for excursion set models of the massive end of the halo mass function. However, it diverges at lower masses. We show how to cure this divergence by using a formulation which requires knowledge of just one other variable. While our analysis is general, we use examples based on Gaussian initial conditions to illustrate our results. Our formulation, which is simple and fast, yields excellent agreement with the considerably more computationally expensive Monte-Carlo solution of the first crossing distribution, for a wide variety of moving barriers, even at very low masses.

If one accounts for correlations between scales, then nonlocal, k-dependent halo bias is part and parcel of the excursion set approach, and hence of halo model predictions for galaxy bias. We present an analysis that distinguishes between a number of different effects, each one of which contributes to scale-dependent bias in real space. We show how to isolate these effects and remove the scale dependence, order by order, by cross-correlating the halo field with suitably transformed versions of the mass field. These transformations may be thought of as simple one-point, two-scale measurements that allow one to estimate quantities which are usually constrained using n-point statistics. As part of our analysis, we present a simple analytic approximation for the first crossing distribution of walks with correlated steps which are constrained to pass through a specified point, and demonstrate its accuracy. Although we concentrate on nonlinear, nonlocal bias with respect to a Gaussian random field, we show how to generalize our analysis to more general fields.

We provide a simple formula that accurately approximates the first crossing distribution of barriers having a wide variety of shapes, by random walks with a wide range of correlations between steps. Special cases of it are useful for estimating halo abundances, evolution, and bias, as well as the nonlinear counts in cells distribution. We discuss how it can be extended to allow for the dependence of the barrier on quantities other than overdensity, to construct an excursion set model for peaks, and to show why assembly and scale dependent bias are generic even at the linear level.

We prove that, in a generic single-field model, the consistency relation for the 3-point function in the squeezed limit receives corrections that vanish quadratically in the ratio of the momenta, i.e. as (k_L/k_S)^2. This implies that a detection of a bispectrum signal going as 1/k_L^2 in the squeezed limit, that is suppressed only by one power of k_L compared with the local shape, would rule out all single-field models. The absence of this kind of terms in the bispectrum holds also for multifield models, but only if all the fields have a mass much smaller than H. The detection of any scale dependence of the bias, for scales much larger than the size of the haloes, would disprove all single-field models. We comment on the regime of squeezing that can be probed by realistic surveys.

We study the effects of primordial non-Gaussianity on the large scale structure in the excursion set approach, accounting for correlations between steps of the random walks in the smoothed initial density field. These correlations are induced by realistic smoothing filters (as opposed to a filter that is sharp in k-space), but have been ignored by many analyses to date. We present analytical arguments -- building on existing results for Gaussian initial conditions -- which suggest that the effect of the filter at large smoothing scales is remarkably simple, and is in fact identical to what happens in the Gaussian case: the non-Gaussian walks behave as if they were smooth and deterministic, or "completely correlated". As a result, the first crossing distribution (which determines, e.g., halo abundances) follows from the single-scale statistics of the non-Gaussian density field -- the so-called "cloud-in-cloud" problem does not exist for completely correlated walks. Also, the answer from single-scale statistics is simply one half that for sharp-k walks. We explicitly test these arguments using Monte Carlo simulations of non-Gaussian walks, showing that the resulting first crossing distributions, and in particular the factor 1/2 argument, are remarkably insensitive to variations in the power spectrum and the defining non-Gaussian process. We also use our Monte Carlo walks to test some of the existing prescriptions for the non-Gaussian first crossing distribution. Since the factor 1/2 holds for both Gaussian and non-Gaussian initial conditions, it provides a theoretical motivation (the first, to our knowledge) for the common practice of analytically prescribing a ratio of non-Gaussian to Gaussian halo abundances.

The abundance of collapsed objects in the universe, or halo mass function, is an important theoretical tool in studying the effects of primordially generated non-Gaussianities on the large scale structure. The non-Gaussian mass function has been calculated by several authors in different ways, typically by exploiting the smallness of certain parameters which naturally appear in the calculation, to set up a perturbative expansion. We improve upon the existing results for the mass function by combining path integral methods and saddle point techniques (which have been separately applied in previous approaches). Additionally, we carefully account for the various scale dependent combinations of small parameters which appear. Some of these combinations in fact become of order unity for large mass scales and at high redshifts, and must therefore be treated non-perturbatively. Our approach allows us to do this, and to also account for multi-scale density correlations which appear in the calculation. We thus derive an accurate expression for the mass function which is based on approximations that are valid over a larger range of mass scales and redshifts than those of other authors. By tracking the terms ignored in the analysis, we estimate theoretical errors for our result and also for the results of others. We also discuss the complications introduced by the choice of smoothing filter function, which we take to be a top-hat in real space, and which leads to the dominant errors in our expression. Finally, we present a detailed comparison between the various expressions for the mass functions, exploring the accuracy and range of validity of each.

Primordial non-Gaussianity (NG) affects the large scale structure (LSS) of the universe by leaving an imprint on the distribution of matter at late times. Much attention has been focused on using the distribution of collapsed objects (i.e. dark matter halos and the galaxies and galaxy clusters that reside in them) to probe primordial NG. An equally interesting and complementary probe however is the abundance of extended underdense regions or voids in the LSS. The calculation of the abundance of voids using the excursion set formalism in the presence of primordial NG is subject to the same technical issues as the one for halos, which were discussed e.g. in arXiv:1005.1203. However, unlike the excursion set problem for halos which involved random walks in the presence of one barrier $\delta_c$, the void excursion set problem involves two barriers $\delta_v$ and $\delta_c$. This leads to a new complication introduced by what is called the "void-in-cloud" effect discussed in the literature, which is unique to the case of voids. We explore a path integral approach which allows us to carefully account for all these issues, leading to a rigorous derivation of the effects of primordial NG on void abundances. The void-in-cloud issue in particular makes the calculation conceptually rather different from the one for halos. However, we show that its final effect can be described by a simple yet accurate approximation. Our final void abundance function is valid on larger scales than the expressions of other authors, while being broadly in agreement with those expressions on smaller scales.

We study the consequences of imposing an approximate Galilean symmetry on the Effective Theory of Inflation, the theory of small perturbations around the inflationary background. This approach allows us to study the effect of operators with two derivatives on each field, which can be the leading interactions due to non-renormalization properties of the Galilean Lagrangian. In this case cubic non-Gaussianities are given by three independent operators, containing up to six derivatives, two with a shape close to equilateral and one peaking on flattened isosceles triangles. The four-point function is larger than in models with small speed of sound and potentially observable with the Planck satellite.

A method to estimate the reliability of a perturbative expansion of the stochastic inflationary Langevin equation is presented and discussed. The method is applied to various inflationary scenarios, as large field, small field and running mass models. It is demonstrated that the perturbative approach is more reliable than could be naively suspected and, in general, only breaks down at the very end of inflation.

A perturbative method for solving the Langevin equation of inflationary cosmology in presence of backreaction is presented. In the Gaussian approximation, the method permits an explicit calculation of the probability distribution of the inflaton field for an arbitrary potential, with or without the volume effects taken into account. The perturbative method is then applied to various concrete models namely large field, small field, hybrid and running mass inflation. New results on the stochastic behavior of the inflaton field in those models are obtained. In particular, it is confirmed that the stochastic effects can be important in new inflation while it is demonstrated they are negligible in (vacuum dominated) hybrid inflation. The case of stochastic running mass inflation is discussed in some details and it is argued that quantum effects blur the distinction between the four classical versions of this model. It is also shown that the self-reproducing regime is likely to be important in this case.

The behavior of the quintessence field is studied during inflation. In order to have a satisfactory model of dark energy, the quintessence field value today should be as insensible to the initial conditions as possible. Usually, only the dependence on the initial conditions specified at the end of inflation or, equivalently, at the beginning of the radiation dominated era, is considered. Provided the quintessence field is initially within a large but, crucially, finite interval, its present value becomes independent of the initial value it started from. The question as to whether inflation naturally drives the quintessence field to the above-mentioned interval is addressed. Since the quantum effects turn out to be important, the formalism of stochastic inflation is used in order to calculate the evolution of the quintessence field. Moreover, the quantum effects originating from the inflaton field are also taken into account and are proved to be sub-dominant in most cases. Finally, the requirement that the quintessence field is on tracks today is shown to imply quite tight constraints on the initial values of the quintessence and inflaton fields at the beginning of inflation. In particular, the initial value of the inflaton field cannot be too large which indicates that the quintessential scenario seems to be compatible with inflation only if the total number of e-folds is quite small.

We generalize the treatment of inflationary perturbations to deal with the non-Markovian colored noise emerging from any realistic approach to stochastic inflation. We provide a calculation of the power-spectrum of the gauge-invariant comoving curvature perturbation to first order in the slow-roll parameters. Properly accounting for the constraint that our local patch of the Universe is homogeneous on scales just above the present Hubble radius, we find a blue tilt of the power-spectrum on the largest observable scales, in agreement with the WMAP data which show an unexpected suppression of the low multipoles of the CMB anisotropy. Our explanation of the anomalous behaviour of the lower multipoles of the CMB anisotropies does not invoke any ad-hoc introduction of new physical ingredients in the theory.

Using the techniques of out-of-equilibrium field theory, we study the influence on the properties of cosmological perturbations generated during inflation on observable scales coming from fluctuations corresponding today to scales much bigger than the present Hubble radius. We write the effective action for the coarse-grained inflaton perturbations integrating out the sub-horizon modes, which manifest themselves as a colored noise and lead to memory effects. Using the simple model of a scalar field with cubic self-interactions evolving in a fixed de Sitter background, we evaluate the two- and three-point correlation function on observable scales. Our basic procedure shows that perturbations do preserve some memory of the super-horizon-scale dynamics, in the form of scale-dependent imprints in the statistical moments. In particular, we find a blue tilt of the power-spectrum on large scales, in agreement with the recent results of the WMAP collaboration which show a suppression of the lower multipoles in the Cosmic Microwave Background anisotropies, and a substantial enhancement of the intrinsic non-Gaussianity on large scales