Researcher profile

Cristiano Porciani

Cristiano Porciani contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 published item(s)

preprint2022arXiv

Mechanical feedback from stellar winds with an application to galaxy formation at high redshift

We compute different sets of stellar evolutionary tracks in order to quantify the energy, mass, and metals yielded by massive main-sequence and post-main-sequence winds. Our aim is to investigate the impact of binary systems and of a metallicity-dependent distribution of initial rotational velocities on the feedback by stellar winds. We find significant changes compared to the commonly used non-rotating, single-star scenario. The largest differences are noticeable at low metallicity, where the mechanical-energy budget is substantially increased. So as to establish the maximal (i.e. obtained by neglecting dissipation in the near circumstellar environment) influence of winds on the early stages of galaxy formation, we use our new feedback estimates to simulate the formation and evolution of a sub-$L_*$ galaxy at redshift 3 (hosted by a dark-matter halo with a mass of $1.8\times 10^{11}$ M$_\odot$) and compare the outcome with simulations in which only supernova (SN) feedback is considered. Accounting for the continuous energy injection by winds reduces the total stellar mass, the metal content, and the burstiness of the star-formation rate as well as of the outflowing gas mass. However, our numerical experiment suggests that the enhanced mechanical feedback from the winds of rotating and binary stars has a limited impact on the most relevant galactic properties compared to the non-rotating single-star scenario. Eventually, we look at the relative abundance between the metals entrained in winds and those ejected by SNe and find that it stays nearly constant within the simulated galaxy and its surrounding halo.

preprint2021arXiv

The reach of next-to-leading-order perturbation theory for the matter bispectrum

We provide a comparison between the matter bispectrum derived with different flavours of perturbation theory at next-to-leading order and measurements from an unprecedentedly large suite of $N$-body simulations. We use the $χ^2$ goodness-of-fit test to determine the range of accuracy of the models as a function of the volume covered by subsets of the simulations. We find that models based on the effective-field-theory (EFT) approach have the largest reach, standard perturbation theory has the shortest, and `classical' resummed schemes lie in between. The gain from EFT, however, is less than in previous studies. We show that the estimated range of accuracy of the EFT predictions is heavily influenced by the procedure adopted to fit the amplitude of the counterterms. For the volumes probed by galaxy redshift surveys, our results indicate that it is advantageous to set three counterterms of the EFT bispectrum to zero and measure the fourth from the power spectrum. We also find that large fluctuations in the estimated reach occur between different realisations. We conclude that it is difficult to unequivocally define a range of accuracy for the models containing free parameters. Finally, we approximately account for systematic effects introduced by the $N$-body technique either in terms of a scale- and shape-dependent bias or by boosting the statistical error bars of the measurements (as routinely done in the literature). We find that the latter approach artificially inflates the reach of EFT models due to the presence of tunable parameters.

preprint2020arXiv

A comparison of $\text{H}_2$ formation models at high redshift

Modelling the molecular gas that is routinely detected through CO observations of high-redshift galaxies constitutes a major challenge for ab initio simulations of galaxy formation. We carry out a suite of cosmological hydrodynamic simulations to compare three approximate methods that have been used in the literature to track the formation and evolution of the simplest and most abundant molecule, H$_2$. Namely, we consider: i) a semi-empirical procedure that associates H$_2$ to dark-matter haloes based on a series of scaling relations inferred from observations, ii) a model that assumes chemical equilibrium between the H$_2$ formation and destruction rates, and iii) a model that fully solves the out-of-equilibrium rate equations and accounts for the unresolved structure of molecular clouds. We study the impact of finite spatial resolution and show that robust H$_2$ masses at redshift $z\approx 4$ can only be obtained for galaxies that are sufficiently metal enriched in which H$_2$ formation is fast. This corresponds to H$_2$ reservoirs with masses $M_{\mathrm{H_2}}\gtrsim 6\times 10^9 \mathrm{M}_\odot$. In this range, equilibrium and non-equilibrium models predict similar molecular masses (but different galaxy morphologies) while the semi-empirical method produces less H$_2$. The star formation rates as well as the stellar and H$_2$ masses of the simulated galaxies are in line with those observed in actual galaxies at similar redshifts that are not massive starbursts. The H$_2$ mass functions extracted from the simulations at $z\approx 4$ agree well with recent observations that only sample the high-mass end. However, our results indicate that most molecular material at high $z$ lies yet undetected in reservoirs with $10^9<M_{\mathrm H_2}<10^{10} \mathrm{M}_\odot$.

preprint2020arXiv

Forecasts for Next Generation tSZ Surveys: the Impact of a Cosmology-Dependent Selection Function

The thermal Sunyaev-Zel&#39;dovich (tSZ) effect is one of the primary tools for finding and characterizing galaxy clusters. Several ground-based experiments are either underway or are being planned for mapping wide areas of the sky at $\sim 150$ GHz with large-aperture telescopes. We present cosmological forecasts for a &#39;straw man&#39; tSZ survey that will observe a sky area between $200$ and $10^4$ deg$^2$ to an rms noise level between 2.8 and 20.2 $μ$K-arcmin. The probes we consider are the cluster number counts (as a function of the integrated Compton-$Y$ parameter and redshift) and their angular clustering (as a function of redshift). At fixed observing time, we find that wider surveys constrain cosmology slightly better than deeper ones due to their increased ability to detect rare high-mass clusters. In all cases, we notice that adding the clustering information does not practically improve the constraints derived from the number counts. We compare forecasts obtained by sampling the posterior distribution with the Markov-chain-Monte-Carlo method against those derived using the Fisher-matrix formalism. We find that the latter produces slightly optimistic constraints where errors are underestimated at the 10 per cent level. Most importantly, we use an analytic method to estimate the selection function of the survey and account for its response to variations of the cosmological parameters in the likelihood function. Our analysis demonstrates that neglecting this effect (as routinely done in the literature) yields artificially tighter constraints by a factor of 2.2 and 1.7 for $σ_8$ and $Ω_\mathrm{M}$, respectively.

preprint2020arXiv

The $n$-point streaming model: how velocities shape correlation functions in redshift space

Starting from first principles, we derive the fundamental equations that relate the $n$-point correlation functions in real and redshift space. Our result generalises the so-called `streaming model&#39; to higher-order statistics: the full $n$-point correlation in redshift-space is obtained as an integral of its real-space counterpart times the joint probability density of $n-1$ relative line-of-sight peculiar velocities. Equations for the connected $n$-point correlation functions are obtained by recursively applying the generalised streaming model for decreasing $n$. Our results are exact within the distant-observer approximation and completely independent of the nature of the tracers for which the correlations are evaluated. Focusing on 3-point statistics, we use an $N$-body simulation to study the joint probability density function of the relative line-of-sight velocities of pairs of particles in a triplet. On large scales, we find that this distribution is approximately Gaussian and that its moments can be accurately computed with standard perturbation theory. We use this information to formulate a phenomenological 3-point Gaussian streaming model. A practical implementation is obtained by using perturbation theory at leading order to approximate several statistics in real space. In spite of this simplification, the resulting predictions for the matter 3-point correlation function in redshift space are in rather good agreement with measurements performed in the simulation. We discuss the limitations of the simplified model and suggest a number of possible improvements. Our results find direct applications in the analysis of galaxy clustering but also set the basis for studying 3-point statistics with future peculiar-velocity surveys and experiments based on the kinetic Sunyaev-Zel&#39;dovich effect.

preprint2020arXiv

Toward a robust inference method for the galaxy bispectrum: likelihood function and model selection

The forthcoming generation of galaxy redshift surveys will sample the large-scale structure of the Universe over unprecedented volumes with high-density tracers. This advancement will make robust measurements of three-point clustering statistics possible. In preparation for this improvement, we investigate how several methodological choices can influence inferences based on the bispectrum about galaxy bias and shot noise. We first measure the real-space bispectrum of dark-matter haloes extracted from 298 N-body simulations covering a volume of approximately $1000 h^{-3} \mathrm{Gpc}^3$. We then fit a series of theoretical models based on tree-level perturbation theory to the numerical data. To achieve this, we estimate the covariance matrix of the measurement errors by using 10,000 mock catalogues generated with the Pinocchio code. We study how the model constraints are influenced by the binning strategy for the bispectrum configurations and by the form of the likelihood function. We also use Bayesian model-selection techniques to single out the optimal theoretical description of our data. We find that a three-parameter bias model combined with Poissonian shot noise is necessary to model the halo bispectrum up to scales of $k_\mathrm{max}\lesssim 0.08 h \mathrm{Mpc}^{-1}$, although fitting formulae that relate the bias parameters can be helpful to reduce the freedom of the model without compromising accuracy. Our data clearly disfavour local Eulerian and local Lagrangian bias models and do not require corrections to Poissonian shot noise. We anticipate that model-selection diagnostics will be particularly useful to extend the analysis to smaller scales as, in this case, the number of model parameters will grow significantly large.

preprint2019arXiv

Renormalisation of linear halo bias in N-body simulations

The interpretation of redshift surveys requires modeling the relationship between large-scale fluctuations in the observed number density of tracers, $δ_\mathrm{h}$, and the underlying matter density, $δ$. Bias models often express $δ_\mathrm{h}$ as a truncated series of integro-differential operators acting on $δ$, each weighted by a bias parameter. Due to the presence of `composite operators&#39; (obtained by multiplying fields evaluated at the same spatial location), the linear bias parameter measured from clustering statistics does not coincide with that appearing in the bias expansion. This issue can be cured by re-writing the expansion in terms of `renormalised&#39; operators. After providing a pedagogical and comprehensive review of bias renormalisation in perturbation theory, we generalize the concept to non-perturbative dynamics and successfully apply it to dark-matter haloes extracted from a large suite of N-body simulations. When comparing numerical and perturbative results, we highlight the effect of the window function employed to smooth the random fields. We then measure the bias parameters as a function of halo mass by fitting a non-perturbative bias model (both before and after applying renormalisation) to the cross spectrum $P_{δ_\mathrm{h}δ}(k)$. Finally, we employ Bayesian model selection to determine the optimal operator set to describe $P_{δ_\mathrm{h}δ}(k)$ for $k<0.2\,h$ Mpc$^{-1}$ at redshift $z=0$. We find that it includes $δ, \nabla^2δ, δ^2$ and the square of the traceless tidal tensor, $s^2$. Considering higher-order terms (in $δ$) leads to overfitting as they cannot be precisely constrained by our data. We also notice that next-to-leading-order perturbative solutions are inaccurate for $k\gtrsim 0.1\,h$ Mpc$^{-1}$.

preprint2010arXiv

Structure formation from non-Gaussian initial conditions: multivariate biasing, statistics, and comparison with N-body simulations

We study structure formation in the presence of primordial non-Gaussianity of the local type with parameters f_NL and g_NL. We show that the distribution of dark-matter halos is naturally described by a multivariate bias scheme where the halo overdensity depends not only on the underlying matter density fluctuation delta, but also on the Gaussian part of the primordial gravitational potential phi. This corresponds to a non-local bias scheme in terms of delta only. We derive the coefficients of the bias expansion as a function of the halo mass by applying the peak-background split to common parametrizations for the halo mass function in the non-Gaussian scenario. We then compute the halo power spectrum and halo-matter cross spectrum in the framework of Eulerian perturbation theory up to third order. Comparing our results against N-body simulations, we find that our model accurately describes the numerical data for wavenumbers k < 0.1-0.3 h/Mpc depending on redshift and halo mass. In our multivariate approach, perturbations in the halo counts trace phi on large scales and this explains why the halo and matter power spectra show different asymptotic trends for k -> 0. This strongly scale-dependent bias originates from terms at leading order in our expansion. This is different from what happens using the standard univariate local bias where the scale-dependent terms come from badly behaved higher-order corrections. On the other hand, our biasing scheme reduces to the usual local bias on smaller scales where |phi| is typically much smaller than the density perturbations. We finally discuss the halo bispectrum in the context of multivariate biasing and show that, due to its strong scale and shape dependence, it is a powerful tool for the detection of primordial non-Gaussianity from future galaxy surveys.

preprint2009arXiv

Lagrangian Statistics of Dark Halos in a LCDM Cosmology

New statistical properties of dark matter halos in Lagrangian space are presented. Tracing back the dark matter particles constituting bound halos resolved in a series of N-body simulations, we measure quantitatively the correlations of the proto-halo&#39;s inertia tensors with the local tidal tensors and investigate how the correlation strength depends on the proto-halo&#39;s sphericity, local density and filtering scale. It is shown that the majority of the proto-halos exhibit strong correlations between the two tensors provided that the tidal field is smoothed on the proto-halo&#39;s mass scale. The correlation strength is found to increase as the proto-halo&#39;s sphericity increases, as the proto-halo&#39;s mass increases, and as the local density becomes close to the critical value, delta_{ec}. It is also found that those peculiar proto-halos which exhibit exceptionally weak correlations between the two tensors tend to acquire higher specific angular momentum in Eulerian space, which is consistent with the linear tidal torque theory. In the light of our results, it is intriguing to speculate a hypothesis that the low surface brightness galaxies observed at present epoch correspond to the peculiar proto-halos with extreme low-sphericity whose inertia tensors are weakly correlated with the local tidal tensors.

preprint2007arXiv

Strong MgII systems in quasar and gamma-ray burst spectra

The incidence of strong MgII systems in gamma-ray burst (GRB) spectra is a few times higher than in quasar (QSO) spectra. We investigate several possible explanations for this effect, including: dust obscuration bias, clustering of the absorbers, different beam sizes of the sources, multiband magnification bias of GRBs, association of the absorbers with the GRB event or the circumburst environment. We find that: i) the incidence rate of MgII systems in QSO spectra could be underestimated by a factor 1.3-2 due to dust obscuration; ii) the equivalent-width distribution of the MgII absorbers along GRBs is consistent with that observed along QSOs thus suggesting that the absorbers are more extended than the beam sizes of the sources; iii) on average, GRB afterglows showing more than one MgII system are a factor of 1.7 brighter than the others, suggesting a lensing origin of the observed discrepancy; iv) gravitational lensing (in different forms, from galaxy lensing to microlensing) can bias high the counts of MgII systems along GRBs if the luminosity functions of the prompt gamma-ray emission and of the optical afterglows have a mean faint-end slope approaching -5/3 -- -2; v) some of the absorbers can be associated with the circumburst environment or produced by supernova remnants unrelated to the GRB event itself but lying in the same star-forming region. With the possible exception of magnification bias, it is unlikely that one of these effects on its own can fully account for the observed counts. However, the combined action of some of them can substantially reduce the statistical significance of the discrepancy.

preprint2000arXiv

A universal angular momentum profile for galactic halos

[Abridged] We study the angular-momentum profiles of a statistical sample of halos drawn from a high-resolution N-body simulation of the LCDM cosmology. We find that the cumulative mass distribution of specific angular momentum, j, in a halo of mass Mv is well fit by a universal function, M(<j) = Mv μj/(j_0+j). This profile is defined by one shape parameter (μor j_0) in addition to the global spin parameter λ. It follows a power-law over most of the mass, and flattens at large j, with the flattening more pronounced for small values of μ. Compared to a uniform sphere in solid-body rotation, most halos have a higher fraction of their mass in the low- and high-j tails of the distribution. The spatial distribution of angular momentum in halos tends to be cylindrical and is well-aligned within each halo for ~80% of the halos. We investigate two ideas for the origin of this profile. The first is based on a revised version of linear tidal-torque theory combined with extended Press-Schechter mass accretion, and the second focuses on j transport in minor mergers. Finally, we briefly explore implications of the M(<j) profile on the formation of galactic disks assuming that j is conserved during an adiabatic baryonic infall. The implied gas density profile deviates from an exponential disk, with a higher density at small radii and a tail extending to large radii. The steep central density profiles may imply disk scale lengths that are smaller than observed. This is reminiscent of the &#34;angular-momentum problem&#34; seen in hydrodynamic simulations, even though we have assumed perfect j conservation. A possible solution is to associate the central excesses with bulge components and the outer regions with extended gaseous disks.