Search

The variations in metallicity and spatial patterns within star-forming regions of galaxies result from diverse physical processes unfolding throughout their evolutionary history, with a particular emphasis in recent events. Analysing MaNGA and \textsc{eagle} galaxies, we discovered an additional dependence of the mass-metallicity relation (MZR) on metallicity gradients ($\nabla_{{\rm (O/H)}}$). Two regimes emerged for low and high stellar mass galaxies, distinctly separated at approximately ${\rm M_{\star}} >10^{9.75}$. Low-mass galaxies with strong positive $\nabla_{{\rm (O/H)}}$ appear less enriched than the MZR median, while those with strong negative gradients are consistently more enriched in both simulated and observed samples. Interestingly, low-mass galaxies with strong negative $\nabla_{{\rm (O/H)}}$ exhibit high star-forming activity, regardless of stellar surface density or $\nabla_{{\rm (O/H)}}$. In contrast, a discrepancy arises for massive galaxies between MaNGA and \textsc{eagle} datasets. The latter exhibit a notable anticorrelation between specific star formation rate and stellar surface density, independent of $\nabla_{{\rm (O/H)}}$, while MaNGA galaxies show this trend mainly for strong positive $\nabla_{{\rm (O/H)}}$. Further investigation indicates that galaxies with strong negative gradients tend to host smaller central black holes in observed datasets, a trend not replicated in simulations. These findings suggest disparities in metallicity recycling and mixing history between observations and simulations, particularly in massive galaxies with varying metallicity gradients. These distinctions could contribute to a more comprehensive understanding of the underlying physics.

Observed and simulated galaxies exhibit correlations between stellar mass, metallicity and morphology. We use the EAGLE cosmological simulation to examine the origin of these correlations for galaxies in the stellar mass range $10^9~\rm{M_\odot} \leqslant\ M_\star \leqslant 10^{10}~\rm{M_\odot}$, and the extent to which they contribute to the scatter in the mass-metallicity relation. We find that rotationally supported disc galaxies have lower metallicity than dispersion supported spheroidal galaxies at a given mass, in agreement with previous findings. In EAGLE, this correlation arises because discs form stars at later times, redshift $z\leqslant 1$, from the accretion of low-metallicity gas, whereas spheroidal galaxies galaxies typically form stars earlier, mainly by consumption of their gas reservoir. The different behaviour reflects the growth of their host dark matter halo: at a given stellar mass, disc galaxies inhabit dark matter haloes with lower mass that formed later compared to the haloes of spheroidal galaxies. Halo concentration plays a secondary role., Comment: 18 pages, 12 figures. Accepted for publication in MNRAS

There is a severe tension between the observed star formation rate (SFR) - stellar mass (${\rm M}_{\star}$) relations reported by different authors at $z = 1-4$. In addition, the observations have not been successfully reproduced by state-of-the-art cosmological simulations which tend to predict a factor of 2-4 smaller SFRs at a fixed ${\rm M}_{\star}$. We examine the evolution of the SFR$-{\rm M}_{\star}$ relation of $z = 1-4 $ galaxies using the SKIRT simulated spectral energy distributions of galaxies sampled from the EAGLE simulations. We derive SFRs and stellar masses by mimicking different observational techniques. We find that the tension between observed and simulated SFR$-{\rm M}_{\star}$ relations is largely alleviated if similar methods are used to infer the galaxy properties. We find that relations relying on infrared wavelengths (e.g. 24 ${\rm \mu m}$, MIPS - 24, 70 and 160 ${\rm \mu m}$ or SPIRE - 250, 350, 500 ${\rm \mu m}$) have SFRs that exceed the intrinsic relation by 0.5 dex. Relations that rely on the spectral energy distribution fitting technique underpredict the SFRs at a fixed stellar mass by -0.5 dex at $z \sim 4$ but overpredict the measurements by 0.3 dex at $z \sim 1$. Relations relying on dust-corrected rest-frame UV luminosities, are flatter since they overpredict/underpredict SFRs for low/high star forming objects and yield deviations from the intrinsic relation from 0.10 dex to -0.13 dex at $z \sim 4$. We suggest that the severe tension between different observational studies can be broadly explained by the fact that different groups employ different techniques to infer their SFRs., Comment: 15 pages, 6 Figures, 4 tables, accepted for publication in MNRAS, Figures 1, 2, 3, 4 summarize results

We present an analysis of the gas dynamics of star-forming galaxies at z~1.5 using data from the KMOS Galaxy Evolution Survey (KGES). We quantify the morphology of the galaxies using $HST$ CANDELS imaging parametrically and non-parametrically. We combine the H$\alpha$ dynamics from KMOS with the high-resolution imaging to derive the relation between stellar mass (M$_{*}$) and stellar specific angular momentum (j$_{*}$). We show that high-redshift star-forming galaxies at z~1.5 follow a power-law trend in specific stellar angular momentum with stellar mass similar to that of local late-type galaxies of the form j$_*$$\propto$M$_*^{0.53 \pm 0.10}$. The highest specific angular momentum galaxies are mostly disc-like, although generally, both peculiar morphologies and disc-like systems are found across the sequence of specific angular momentum at a fixed stellar mass. We explore the scatter within the j$_{*}$-M$_{*}$ plane and its correlation with both the integrated dynamical properties of a galaxy (e.g. velocity dispersion, Toomre Q$_{\rm g}$, H$\alpha$ star formation rate surface density $\Sigma_{\rm SFR}$) and its parameterised rest-frame UV/optical morphology (e.g. S\'ersic index, bulge to total ratio, Clumpiness, Asymmetry and Concentration). We establish that the position in the j$_{*}$-M$_{*}$ plane is strongly correlated with the star-formation surface density and the Clumpiness of the stellar light distribution. Galaxies with peculiar rest-frame UV/optical morphologies have comparable specific angular momentum to disc-dominated galaxies of the same stellar mass, but are clumpier and have higher star-formation rate surface densities. We propose that the peculiar morphologies in high--redshift systems are driven by higher star formation rate surface densities and higher gas fractions leading to a more clumpy inter-stellar medium., Comment: Accepted for publication in MNRAS. Comments welcome. 25 pages, 14 figures. The full pdf, which includes an additional 17 pages of supplementary figures is available on this http://astro.dur.ac.uk/KROSS/publications.html

The metallicity of star-forming gas in galaxies from the EAGLE simulations increases with stellar mass. Here we investigate whether the scatter around this relation correlates with morphology and/or stellar kinematics. At redshift $z=0$, galaxies with more rotational support have lower metallicities on average when the stellar mass is below $M_\star\approx 10^{10}~{\rm M}_\odot$. This trend inverts at higher values of $M_\star$, when prolate galaxies show typically lower metallicity. At increasing redshifts, the trend between rotational support and metallicity becomes weaker at low stellar mass but more pronounced at high stellar mass. We argue that the secondary dependence of metallicity on stellar kinematics is another manifestation of the observed anti-correlation between metallicity and star formation rate at a given stellar mass. At low masses, such trends seem to be driven by the different star-formation histories of galaxies and stellar feedback. At high masses, feedback from active galactic nuclei and galaxy mergers play a dominant role., Comment: 5 pages, 4 figures. Accepted for publication in MNRAS

We present the goals, design, and first results of the MUSE Ultra Deep Field (MUDF) survey, a large programme using the Multi Unit Spectroscopic Explorer (MUSE) instrument at the ESO Very Large Telescope. The MUDF survey is collecting ~ 150 hours on-source of integral field optical spectroscopy in a 1.5 x 1.2 square arcmin region which hosts several astrophysical structures along the line of sight, including two bright z ~ 3.2 quasars with close separation (~ 500 kpc). Following the description of the data reduction procedures, we present the analysis of the galaxy environment and gaseous properties of seven groups detected at redshifts 0.5 < z < 1.5, spanning a large dynamic range in halo mass, log(Mh/Msun) ~ 11 - 13.5. For four of the groups, we find associated MgII absorbers tracing cool gas in high-resolution spectroscopy of the two quasars, including one case of correlated absorption in both sightlines at distance ~ 480 kpc. The absorption strength associated with the groups is higher than what has been reported for more isolated galaxies of comparable mass and impact parameters. We do not find evidence for widespread cool gas giving rise to strong absorption within these groups. Combining these results with the distribution of neutral and ionised gas seen in emission in lower-redshift groups, we conclude that gravitational interactions in the group environment strip gas from the galaxy haloes into the intragroup medium, boosting the cross section of cool gas and leading to the high fraction of strong MgII absorbers that we detect., Comment: 20 pages, 12 figures, 4 tables, Accepted for publication in MNRAS

We present adaptive optics assisted integral field spectroscopy of 34 star-forming galaxies at $z$ = 0.8-3.3 selected from the HiZELS narrow-band survey. We measure the kinematics of the ionised interstellar medium on $\sim$1 kpc scales, and show that the galaxies are turbulent, with a median ratio of rotational to dispersion support of $v$/$\sigma$=0.82$\pm$0.13. We combine the dynamics with high-resolution rest-frame optical imaging and extract emission line rotation curves. We show that high-redshift star-forming galaxies follow a similar power-law trend in specific angular momentum with stellar mass as that of local late type galaxies. We exploit the high resolution of our data and examine the radial distribution of angular momentum within each galaxy by constructing total angular momentum profiles. Although the stellar mass of a typical star-forming galaxy is expected to grow by a factor $\sim$8 in the $\sim$5 Gyrs between $z$$\sim$3.3 and $z$$\sim$0.8, we show that the internal distribution of angular momentum becomes less centrally concentrated in this period i.e the angular momentum grows outwards. To interpret our observations, we exploit the EAGLE simulation and trace the angular momentum evolution of star forming galaxies from $z$$\sim$3 to $z$$\sim$0, identifying a similar trend of decreasing angular momentum concentration. This change is attributed to a combination of gas accretion in the outer disk, and feedback that preferentially arises from the central regions of the galaxy. We discuss how the combination of the growing bulge and angular momentum stabilises the disk and gives rise to the Hubble sequence., Comment: 31 pages, 18 figures including appendices. Accepted for publication in MNRAS 2019 March 12. Updated text and figures

We present first results from Multi Unit Spectroscopic Explorer (MUSE) observations at the Very Large Telescope in the MUSE Ultra Deep Field (MUDF), a $\approx 1.2\times 1.4$ arcmin$^2$ region for which we are collecting $\approx$200 hours of integral field spectroscopy. The $\approx 40$-hour observation completed to date reveals the presence of a group of three Ly$\alpha$ nebulae associated with a bright quasar pair at $z\simeq3.23$ with projected separation of $\approx 500\rm~kpc$. Two of the nebulae are physically associated with the quasars which are likely powering the Ly$\alpha$ emission, and extend for $\gtrsim 100~\rm kpc$ at a surface brightness level of $\approx 6\times 10^{-19}~\rm erg~s^{-1}~cm^{-2}~arcsec^{-2}$. A third smaller ($\approx$35 kpc) nebula lies at a velocity offset of $\approx 1550$ km s$^{-1}$. Despite their clustered nature, the two large nebulae have properties similar to those observed in isolated quasars and exhibit no sharp decline in flux at the current depth, suggesting an even more extended distribution of gas around the quasars. We interpret the shape and the alignment of the two brighter nebulae as suggestive of the presence of an extended structure connecting the two quasar host galaxies, as seen for massive galaxies forming within gas-rich filaments in cosmological simulations., Comment: MNRAS Letters in press, 3 figures, revised version: minor changes following proofs

We analyse maps of the spatially-resolved nebular emission of $\approx$1500 star-forming galaxies at $z\approx0.6$-$2.2$ from deep KMOS and MUSE observations to measure the average shape of their rotation curves. We use these to test claims for declining rotation curves at large radii in galaxies at $z\approx1$-$2$ that have been interpreted as evidence for an absence of dark matter. We show that the shape of the average rotation curves, and the extent to which they decline beyond their peak velocities, depends upon the normalisation prescription used to construct the average curve. Normalising in size by the galaxy stellar disk-scale length after accounting for seeing effects ($R_{\rm{d}}^{\prime}$), we construct stacked position-velocity diagrams that trace the average galaxy rotation curve out to $6R_{\rm{d}}^{\prime}$ ($\approx$13 kpc, on average). Combining these curves with average HI rotation curves for local systems, we investigate how the shapes of galaxy rotation curves evolve over $\approx$10 Gyr. The average rotation curve for galaxies binned in stellar mass, stellar surface mass density and/or redshift is approximately flat, or continues to rise, out to at least $6R_{\rm{d}}^{\prime}$. We find a trend between the outer slopes of galaxies' rotation curves and their stellar mass surface densities, with the higher surface density systems exhibiting flatter rotation curves. Drawing comparisons with hydrodynamical simulations, we show that the average shapes of the rotation curves for our sample of massive, star-forming galaxies at $z\approx0$-$2.2$ are consistent with those expected from $\Lambda$CDM theory and imply dark matter fractions within $6R_{\rm{d}}$ of at least $\approx60$ percent., Comment: 32 pages, 21 figures including appendices. Accepted for publication in MNRAS after addressing final comments from referee. Additional appendices added, plus some updates and additions to plots and text

We present sensitive 870$\mu$m continuum measurements from our ALMA programmes of 114 X-ray selected AGN in the CDF-S and COSMOS fields. We use these observations in combination with data from Spitzer and Herschel to construct a sample of 86 X-ray selected AGN, 63 with ALMA constraints at $z=1.5-3.2$ with stellar mass $>2\times10^{10}M_{\odot}$. We constructed broad-band spectral energy distributions in the infrared band (8 - 1000$\mu$m) and constrain star-formation rates (SFRs) uncontaminated by the AGN. Using a hierarchical Bayesian method that takes into account the information from upper limits, we fit SFR and specific SFR (sSFR) distributions. We explore these distributions as a function of both X-ray luminosity and stellar mass. We compare our measurements to two versions of the EAGLE hydrodynamical simulations: the reference model with AGN feedback and the model without AGN. We find good agreement between the observations and that predicted by the EAGLE reference model for the modes and widths of the sSFR distributions as a function of both X-ray luminosity and stellar mass; however, we found that the EAGLE model without AGN feedback predicts a significantly narrower width when compared to the data. Overall, from the combination of the observations with the model predictions, we conclude that (1) even with AGN feedback, we expect no strong relationship between the sSFR distribution parameters and instantaneous AGN luminosity and (2) a signature of AGN feedback is a broad distribution of sSFRs for all galaxies (not just those hosting an AGN) with stellar masses above $\approx 10^{10}$M$_{\odot}$., Comment: Accepted for publication in MNRAS

This work presents a study of galactic outflows driven by stellar feedback. We extract main sequence disc galaxies with stellar mass $10^9\le$ M$_{\star}/$M$_{\odot} \le 5.7\times10^{10}$ at redshift $z=0$ from the highest resolution cosmological simulation of the Evolution and Assembly of GaLaxies and their Environments (EAGLE) set. Synthetic gas rotation velocity and velocity dispersion ($\sigma$) maps are created and compared to observations of disc galaxies obtained with the Sydney-AAO Multi-object Integral field spectrograph (SAMI), where $\sigma$-values greater than $150$ km s$^{-1}$ are most naturally explained by bipolar outflows powered by starburst activity. We find that the extension of the simulated edge-on (pixelated) velocity dispersion probability distribution depends on stellar mass and star formation rate surface density ($\Sigma_{\rm SFR}$), with low-M$_{\star}/$low-$\Sigma_{\rm SFR}$ galaxies showing a narrow peak at low $\sigma$ ($\sim30$ km s$^{-1}$) and more active, high-M$_{\star}/$high-$\Sigma_{\rm SFR}$ galaxies reaching $\sigma>150$ km s$^{-1}$. Although supernova-driven galactic winds in the EAGLE simulations may not entrain enough gas with T $<10^5$ K compared to observed galaxies, we find that gas temperature is a good proxy for the presence of outflows. There is a direct correlation between the thermal state of the gas and its state of motion as described by the $\sigma$-distribution. The following equivalence relations hold in EAGLE: $i)$ low-$\sigma$ peak $\,\Leftrightarrow\,$ disc of the galaxy $\,\Leftrightarrow\,$ gas with T $<10^5$ K; $ii)$ high-$\sigma$ tail $\,\Leftrightarrow\,$ galactic winds $\,\Leftrightarrow\,$ gas with T $\ge 10^5$ K., Comment: 18 pages, 12 figures, accepted for publication in MNRAS

We investigate the evolution of the galaxy Star Formation Rate Function (SFRF) and Cosmic Star Formation Rate Density (CSFRD) of $z\sim 0-8 $ galaxies in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulations. In addition, we present a compilation of UV, IR and H$\alpha$ SFRFs and compare these with the predictions from the EAGLE suite of cosmological hydrodynamic simulations. We find that the constraints implied by different indicators are inconsistent with each other for the highest star-forming objects at z < 2, a problem that is possibly related to selection biases and the uncertainties of dust attenuation effects. EAGLE's feedback parameters were calibrated to reproduce realistic galaxy sizes and stellar masses at z = 0.1. In this work we test if and why those choices yield realistic Star Formation Rates (SFRs) for $z \sim 0-8$ as well. We demonstrate that SNe feedback plays a major role at setting the abundance of galaxies at all star-forming regimes, especially at high redshifts. On the contrary, Active Galactic Nuclei (AGN) feedback becomes more prominent at lower redshifts and is a major mechanism that affects only the highest star-forming systems. Furthermore, we find that galaxies with SFR $\sim 1-10 \, {\rm M_{\odot} \, yr^{-1}}$ dominate the CSFRD at redshifts z < 5, while rare high star-forming galaxies (SFR $\sim 10-100 \,{\rm M_{\odot} \, yr^{-1}}$) contribute significantly only briefly around the peak era ($z \sim 2$) and then are quenched by AGN feedback. In the absence of this prescription objects with SFR $\sim 10-100 \,{\rm M_{\odot} \, yr^{-1}}$ would dominate the CSFRD, while the cosmic budget of star formation would be extremely high. Finally, we demonstrate that the majority of the cosmic star formation occurs in relatively rare high mass halos ($ {\rm M_{Halo}} \sim 10^{11-13} \, {\rm M_{\odot}}$) even at the earliest epochs., Comment: 21 pages, 12 figures, accepted for publication in MNRAS

We present the evolution of galaxy sizes, from redshift 2 to 0, for actively star forming and passive galaxies in the cosmological hydrodynamical 1003 cMpc3 simulation of the EAGLE project. We find that the sizes increase with stellar mass , but that the relation weakens with increasing redshift. Separating galaxies by their star formation activity, we find that passive galaxies are typically smaller than active galaxies at fixed stellar mass. These trends are consistent with those found in observations and the level of agreement between the predicted and observed size - mass relation is of order 0.1 dex for z < 1 and 0.2-0.3 dex from redshift 1 to 2. We use the simulation to compare the evolution of individual galaxies to that of the population as a whole. While the evolution of the size-stellar mass relation for active galaxies provides a good proxy for the evolution of individual galaxies, the evolution of individual passive galaxies is not well represented by the observed size - mass relation due to the evolving number density of passive galaxies. Observations of z \approx 2 galaxies have revealed an abundance of massive red compact galaxies, that depletes below z \approx 1. We find that a similar population forms naturally in the simulation. Comparing these galaxies to their z = 0 descendants, we find that all compact galaxies grow in size due to the high-redshift stars migrating outwards. Approximately 60% of the compact galaxies increase in size further due to renewed star formation and/or mergers., Comment: 14 pages, 10 figures

We study high-resolution hydrodynamic simulations of Milky Way type galaxies obtained within the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) project, and identify the those that best satisfy observational constraints on the Milky Way total stellar mass, rotation curve, and galaxy shape. Contrary to mock galaxies selected on the basis of their total virial mass, the Milky Way analogues so identified consistently exhibit very similar dark matter profiles inside the solar circle, therefore enabling more accurate predictions for indirect dark matter searches. We find in particular that high resolution simulated haloes satisfying observational constraints exhibit, within the inner few kiloparsecs, dark matter profiles shallower than those required to explain the so-called Fermi GeV excess via dark matter annihilation., Comment: 33 pages, 15 figures, 4 tables. Matches published version

We investigate the evolution of galaxy masses and star formation rates in the Evolution and Assembly of Galaxies and their Environment (EAGLE) simulations. These comprise a suite of hydrodynamical simulations in a $\Lambda$CDM cosmogony with subgrid models for radiative cooling, star formation, stellar mass loss, and feedback from stars and accreting black holes. The subgrid feedback was calibrated to reproduce the observed present-day galaxy stellar mass function and galaxy sizes. Here we demonstrate that the simulations reproduce the observed growth of the stellar mass density to within 20 per cent. The simulation also tracks the observed evolution of the galaxy stellar mass function out to redshift z = 7, with differences comparable to the plausible uncertainties in the interpretation of the data. Just as with observed galaxies, the specific star formation rates of simulated galaxies are bimodal, with distinct star forming and passive sequences. The specific star formation rates of star forming galaxies are typically 0.2 to 0.4 dex lower than observed, but the evolution of the rates track the observations closely. The unprecedented level of agreement between simulation and data makes EAGLE a powerful resource to understand the physical processes that govern galaxy formation., Comment: Accepted for publication by MNRAS; 18 pages, 7 figures; comments welcome

We use the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) cosmological hydrodynamic simulation at z=0 to study the distribution and environmental dependence of neutral hydrogen (HI) gas in the outskirts of simulated galaxies. This gas can currently be probed directly in, for example, Ly$\alpha$ absorption via the observation of background quasars. Radio facilities, such as the Square Kilometre Array, will provide a complementary probe of the diffuse HI in emission and will constrain the physics underpinning the complex interplay between accretion and feedback mechanisms which affect the intergalactic medium. We extract a sample of 488 galaxies from a re-simulation of the average cosmic density GIMIC region. We estimate the neutral hydrogen content of these galaxies and the surrounding intergalactic medium within which they reside. We investigate the average HI radial profiles by stacking the individual profiles according to both mass and environment. We find high HI column densities at large impact parameters in group environments and markedly lower HI densities for non-group galaxies. We suggest that these results likely arise from the combined effects of ram pressure stripping and tidal interactions present in group environments., Comment: 9 pages, 3 figures, MNRAS accepted

Feedback from energy liberated by gas accretion onto black holes (BHs) is an attractive mechanism to explain the exponential cut-off at the massive end of the galaxy stellar mass function (SMF). Semi-analytic models of galaxy formation in which this form of feedback is assumed to suppress cooling in haloes where the gas cooling time is large compared to the dynamical time do indeed achieve a good match to the observed SMF. Furthermore, hydrodynamic simulations of individual halos in which gas is assumed to accrete onto the central BH at the Bondi rate have shown that a self-regulating regime is established in which the BH grows just enough to liberate an amount of energy comparable to the thermal energy of the halo. However, this process is efficient at suppressing the growth not only of massive galaxies but also of galaxies like the Milky Way, leading to disagreement with the observed SMF. The Bondi accretion rate, however, is inappropriate when the accreting material has angular momentum. We present an improved accretion model that takes into account the circularisation and subsequent viscous transport of infalling material and include it as a "subgrid" model in hydrodynamic simulations of the evolution of halos with a wide range of masses. The resulting accretion rates are generally low in low mass ($\lsim 10^{11.5} \msun$) halos, but show outbursts of Eddington-limited accretion during galaxy mergers. During outbursts these objects strongly resemble quasars. In higher mass haloes, gas accretion occurs continuously, typically at $~10$ % of the Eddington rate, which is conducive to the formation of radio jets. The resulting dependence of the accretion behaviour on halo mass induces a break in the relation between galaxy stellar mass and halo mass in these simulations that matches observations., Comment: 24 pages, 14 figures. Accepted to MNRAS for publication. Updated version with typos corrected

We present a new parallel implementation of the PINpointing Orbit Crossing-Collapsed HIerarchical Objects (PINOCCHIO) algorithm, a quick tool, based on Lagrangian Perturbation Theory, for the hierarchical build-up of Dark Matter halos in cosmological volumes. To assess its ability to predict halo correlations on large scales, we compare its results with those of an N-body simulation of a 3 Gpc/h box sampled with 2048^3 particles taken from the MICE suite, matching the same seeds for the initial conditions. Thanks to the FFTW libraries and to the relatively simple design, the code shows very good scaling properties. The CPU time required by PINOCCHIO is a tiny fraction (~1/2000) of that required by the MICE simulation. Varying some of PINOCCHIO numerical parameters allows one to produce a universal mass function that lies in the range allowed by published fits, although it underestimates the MICE mass function of FoF halos in the high mass tail. We compare the matter-halo and the halo-halo power spectra with those of the MICE simulation and find that these 2-point statistics are well recovered on large scales. In particular, when catalogs are matched in number density, agreement within ten per cent is achieved for the halo power spectrum. At scales k>0.1 h/Mpc, the inaccuracy of the Zel'dovich approximation in locating halo positions causes an underestimate of the power spectrum that can be modeled as a Gaussian factor with a damping scale of d=3 Mpc/h at z=0, decreasing at higher redshift. Finally, a remarkable match is obtained for the reduced halo bispectrum, showing a good description of nonlinear halo bias. Our results demonstrate the potential of PINOCCHIO as an accurate and flexible tool for generating large ensembles of mock galaxy surveys, with interesting applications for the analysis of large galaxy redshift surveys., Comment: 14 pages, 9 PDF figures, accepted by Monthly Notices of the Royal Astronomical Society. The PINOCCHIO code can be downloaded here: http://adlibitum.oats.inaf.it/monaco/Homepage/Pinocchio/index.html

We use observations and simulation to study the relationship between star-forming galaxies and the intergalactic medium at z~3. The observed galaxy sample is based on spectroscopic redshift data from a combination of the VLT LBG Redshift Survey and Keck observations in fields centred on bright z>3 QSOs, whilst the simulation data is taken from GIMIC. In the simulation, we find that the dominant peculiar velocities are in the form of large-scale coherent motions of gas and galaxies. Gravitational infall of galaxies towards one another is also seen. At smaller scales, the peculiar velocities in the simulation over-predict the difference between the simulated real- and z-space galaxy correlation functions. Peculiar velocity pairs separated by <1Mpc/h have a smaller dispersion and explain the z-space correlation function better. The Ly{\alpha} auto- and cross-correlation functions in the GIMIC simulation show infall smaller than implied by previous work. This reduced infall may be due to the galaxy wide outflows implemented in the simulation. The main challenge in comparing these simulated results with the observed correlation functions comes from the presence of velocity errors for the observed LBGs which dominate at ~1Mpc/h scales. When these are taken into account, the observed LBG correlation function is well matched by a simulated $M_*>10^9M_\odot$ galaxy sample. The simulated cross-correlation shows similar neutral gas densities around galaxies as are seen in the observations. The simulated and observed Ly{\alpha} z-space autocorrelation functions agree well with each other. Our overall conclusion is that gas and galaxy peculiar velocities are towards the low end of expectation. Finally, little direct evidence is seen in either simulation or observations for high transmission near galaxies due to feedback. (Abridged), Comment: 24 pages, 22 figures, Published in MNRAS

We present a catalogue of 2135 galaxy redshifts from the VLT LBG Redshift Survey (VLRS), a spectroscopic survey of z ~ 3 galaxies in wide fields centred on background quasi-stellar objects. We have used deep optical imaging to select galaxies via the Lyman-break technique. Spectroscopy of the Lyman-break galaxies (LBGs) was then made using the Visible Multi-Object Spectrograph (VIMOS), giving a mean redshift of z=2.79. We analyse the clustering properties of the VLRS sample and also of the VLRS sample combined with the smaller area Keck-based survey of Steidel et al. From the semiprojected correlation function, wp({\sigma}) we find that the results are well fit with a single power-law model, with clustering scale lengths of r0=3.46+-0.41 and 3.83+-0.24 Mpc/h, respectively. We note that the corresponding combined {\xi}(r) slope is flatter than for local galaxies at {\gamma} = 1.5-1.6 rather than {\gamma}=1.8. This flat slope is confirmed by the z-space correlation function, {\xi}(s), and in the range 10

We compare the results of thirteen cosmological gasdynamical codes used to simulate the formation of a galaxy in the LCDM structure formation paradigm. The various runs differ in their hydrodynamical treatment (SPH, moving-mesh and AMR) but share the same initial conditions and adopt their latest published model of cooling, star formation and feedback. Despite the common halo assembly history, we find large code-to-code variations in the stellar mass, size, morphology and gas content of the galaxy at z=0, due mainly to the different implementations of feedback. Compared with observation, most codes tend to produce an overly massive galaxy, smaller and less gas-rich than typical spirals, with a massive bulge and a declining rotation curve. A stellar disk is discernible in most simulations, though its prominence varies widely from code to code. There is a well-defined trend between the effects of feedback and the severity of the disagreement with observation. Models that are more effective at limiting the baryonic mass of the galaxy come closer to matching observed galaxy scaling laws, but often to the detriment of the disk component. Our conclusions hold at two different numerical resolutions. Some differences can also be traced to the numerical techniques: more gas seems able to cool and become available for star formation in grid-based codes than in SPH. However, this effect is small compared to the variations induced by different feedback prescriptions. We conclude that state-of-the-art simulations cannot yet uniquely predict the properties of the baryonic component of a galaxy, even when the assembly history of its host halo is fully specified. Developing feedback algorithms that can effectively regulate the mass of a galaxy without hindering the formation of high-angular momentum stellar disks remains a challenge., Comment: 28 pages, 16 figures, MNRAS accepted. High resolution figures and additional images can be found at http://www.aip.de/People/cscannapieco/aquila/

The formation of galaxies is regulated by a balance between the supply of gas and the rate at which it is ejected. Traditional explanations of gas ejection equate the energy required to escape the galaxy or host halo to an estimate for the energy yield from supernovae. This yield is usually assumed to be a constant fraction of the total available from the supernova, or is derived from the assumption of a consistent momentum yield. By applying these ideas in the context of a cold dark matter cosmogony, we derive a 1st-order analytic connection between these working assumptions and the expected relationship between baryon content and galaxy circular velocity, and find that these quick predictions straddle recent observational estimates. To examine the premises behind these theories in more detail, we then explore their applicability to a set of gasdynamical simulations of idealised galaxies. We show that different premises dominate to differing degrees in the simulated outflow, depending on the mass of the system and the resolution with which it is simulated. Using this study to anticipate the emergent behaviour at arbitrarily high resolution, we motivate more comprehensive analytic model which allows for the range of velocities with which the gas may exit the system, and incorporates both momentum and energy-based constraints on the outflow. Using a trial exit velocity distribution, this is shown to be compatible with the observed baryon fractions in intermediate-mass systems, but implies that current estimates for low-mass systems can not be solely accounted for by supernova winds under commonly-held assumptions., Comment: 16 pages, 13 figures. Accepted to MNRAS

(Abridged) We study the phase-space distribution of satellite galaxies associated with late-type galaxies in the GIMIC suite of simulations. GIMIC consists of re-simulations of 5 cosmologically representative regions from the Millennium simulation, which have higher resolution and incorporate baryonic physics. Whilst the disc of the galaxy is well aligned with the inner regions (r ~ 0.1r_200) of the dark matter halo, both in shape and angular momentum, there can be substantial misalignments at larger radii (r ~ r_200). Misalignments of > 45 deg are seen in ~ 30% of our sample. We find that the satellite population aligns with the shape (and angular momentum) of the outer dark matter halo. However, the alignment with the galaxy is weak owing to the mismatch between the disc and dark matter halo. Roughly 20% of the satellite systems with ten bright galaxies within r_200 exhibit a polar spatial alignment with respect to the galaxy --- an orientation reminiscent of the classical satellites of the Milky Way. We find that a small fraction (~10%) of satellite systems show evidence for rotational support which we attribute to group infall. There is a bias towards satellites on prograde orbits relative to the spin of the dark matter halo (and to a lesser extent with the angular momentum of the disc). This preference towards co-rotation is stronger in the inner regions of the halo where the most massive satellites accreted at relatively early times are located. We attribute the anisotropic spatial distribution and angular momentum bias of the satellites at z=0 to their directional accretion along the major axes of the dark matter halo. The satellite galaxies have been accreted relatively recently compared to the dark matter mass and have experienced less phase-mixing and relaxation --- the memory of their accretion history can remain intact to z=0., Comment: 19 pages, 17 Figures. Accepted by MNRAS

We have measured redshifts for 243 z ~3 quasars in nine VLT VIMOS LBG redshift survey areas, each of which is centred on a known bright quasar. Using spectra of these quasars, we measure the cross-correlation between neutral hydrogen gas causing the Lya forest and 1020 Lyman-break galaxies at z ~3. We find an increase in neutral hydrogen absorption within 5 h^-1 Mpc of a galaxy in agreement with the results of Adelberger et al. (2003, 2005). The Lya-LBG cross-correlation can be described by a power-law on scales larger than 3 h^-1 Mpc. When galaxy velocity dispersions are taken into account our results at smaller scales (<2 h^-1 Mpc) are also in good agreement with the results of Adelberger et al. (2005). There is little immediate indication of a region with a transmission spike above the mean IGM value which might indicate the presence of star-formation feedback. To measure the galaxy velocity dispersions, which include both intrinsic LBG velocity dispersion and redshift errors, we have used the LBG-LBG redshift space distortion measurements of Bielby et al. (2010). We find that the redshift-space transmission spike implied in the results of Adelberger et al. (2003) is too narrow to be physical in the presence of the likely LBG velocity dispersion and is likely to be a statistical fluke. Nevertheless, neither our nor previous data can rule out the presence of a narrow, real-space transmission spike, given the evidence of the increased Lya absorption surrounding LBGs which can mask the spike's presence when convolved with a realistic LBG velocity dispersion. Finally, we identify 176 CIV systems in the quasar spectra and find an LBG-CIV correlation strength on scales of 10 h^-1 Mpc consistent with the relation measured at ~Mpc scales., Comment: 33 pages, 16 figures, accepted by MNRAS

We present the initial imaging and spectroscopic data acquired as part of the VLT VIMOS Lyman-break galaxy Survey. UBR (or UBVI) imaging covers five 36'x36' fields centred on bright z>3 QSOs, allowing ~21,000 22 LBGs and 10 z>2 QSOs from a total of 19 VIMOS pointings. From the galaxy spectra, we observe a 625+-510 km/s velocity offset between the ISM absorption and Ly-alpha emission line redshifts. Using the photometric and spectroscopic catalogues, we have analysed the galaxy clustering at z~3. In the photometric case, the angular correlation function, w(theta), is well fit by a double power-law with clustering scale-length, r_0 = 3.19+0.32-0.54 Mpc/h for r < 1 Mpc/h and r_0 = 4.59+0.31-0.33 Mpc/h at larger scales. Using the redshift sample we estimate the semi-projected correlation function, w_p(sigma) and find r_0 = 3.67+0.23-0.24 Mpc/h for the VLT sample and r_0 = 3.98+0.14-0.15 Mpc/h for a combined VLT+Keck sample. From the z-space correlation functions and assuming the above xi(r) models, we find that the combined VLT and Keck surveys require a galaxy velocity dispersion, ^1/2 ~ 700 km/s, higher than the ~400 km/s found by previous authors. We also measure a value for the gravitational growth rate parameter of beta(z=3) = 0.48+-0.17, implying a low value for the bias of b = 2.06+1.1-0.5. This value is consistent with the galaxy clustering amplitude which gives b = 2.22+-0.16, assuming the standard cosmology, implying that the evolution of the gravitational growth rate is also consistent with Einstein gravity. We have compared our LBG clustering amplitudes with lower redshift measurements and find that the clustering strength is not inconsistent with that of low-redshift L* spirals for simple 'long-lived' galaxy models., Comment: 27 pages, 25 figures, accepted in MNRAS

(Abridged) The thermal state of the intracluster medium results from a competition between gas cooling and heating. The heating comes from two distinct sources: gravitational heating from the collapse of the dark matter halo and thermal input from galaxy/black hole formation. However, a long standing problem has been that cosmological simulations based on smoothed particle hydrodynamics (SPH) and Eulerian mesh codes predict different results even when cooling and galaxy/black hole heating are switched off. Clusters formed in SPH simulations show near powerlaw entropy profiles, while those formed in mesh simulations develop a core and do not allow gas to reach such low entropies. Since the cooling rate is closely connected to the minimum entropy of the gas, the differences are of potentially key importance. In this paper, we investigate the origin of this discrepancy. By comparing simulations run using the GADGET-2 SPH code and the FLASH adaptive Eulerian mesh code, we show that the discrepancy arises during the idealised merger of two clusters. The difference is not sensitive to the resolution of our simulations, nor is it is due differences in the gravity solvers, Galilean non-invariance of the mesh code, or an effect of unsuitable artificial viscosity in the SPH code. Instead, we find that the difference is inherent to the treatment of eddies and fluid instabilities. These are suppressed in the SPH simulations, while the cluster mergers generate strong vortices in the mesh simulations that efficiently mix the fluid and erase the low entropy gas. Consequently, particles in the SPH simulations retain a close connection to their initial entropy, while this connection is much weaker in the mesh simulations. We discuss the potentially profound implications of these results., Comment: 19 pages, 10 figures, Accepted for publication in MNRAS

We describe a numerical implementation of star formation in disk galaxies, in which the conversion of cooling gas to stars in the multiphase interstellar medium is governed by the rate at which molecular clouds are formed and destroyed. In the model, clouds form from thermally unstable ambient gas and get destroyed by feedback from massive stars and thermal conduction. Feedback in the ambient phase cycles gas into a hot galactic fountain or wind. We model the ambient gas hydrodynamically using smoothed particle hydrodynamics (SPH). However, we cannot resolve the Jeans mass in the cold and dense molecular gas and, therefore, represent the cloud phase with ballistic particles that coagulate when colliding. We show that this naturally produces a multiphase medium with cold clouds, a warm disk, hot supernova bubbles and a hot, tenuous halo. Our implementation of this model is based on the Gadget N-Body code. We illustrate the model by evolving an isolated Milky Way-like galaxy and study the properties of a disk formed in a rotating spherical collapse. Many observed properties of disk galaxies are reproduced well, including the molecular cloud mass spectrum, the molecular fraction as a function of radius, the Schmidt law, the stellar density profile and the appearance of a galactic fountain., Comment: 26 pages including 22 figures. Accepted for publication in MNRAS. Version including full resolution figures available from http://star-www.dur.ac.uk/~dph3cmb/sp.ps

We investigate the effect of peculiar velocities on the redshift space distribution of z>~2 galaxies, and we focus in particular on Ly-alpha emitters. We generate catalogues of dark matter (DM) halos and identify emitters with halos of the same co-moving space density (M(Ly-alpha emitters) ~ 3x10^11 M_sun). We decompose the peculiar velocity field of halos into streaming, gradient and random components, and compute and analyse these as a function of scale. Streaming velocities are determined by fluctuations on very large scales, strongly affected by sample variance, but have a modest impact on the interpretation of observations. Gradient velocities are the most important as they distort structures in redshift space, changing the thickness and orientation of sheets and filaments. Random velocities are typically below or of the same order as the typical observational uncertainty on the redshift. We discuss the importance of these effects for the interpretation of data on the large-scale structure as traced by Ly-alpha emitters (or similar kinds of astrophysical high-redshift objects), focusing on the induced errors in the viewing angles of filaments. We compare our predictions of velocity patterns for Ly-alpha emitters to observations and find that redshift clumping of Ly-alpha emitters, as reported for instance in the fields of high-redshift radio galaxies, does not allow to infer whether an observed field is sampling an early galaxy overdensity., Comment: 10 pages, postscript figures included, uses aa.cls. In press on Astronomy & Astrophysics

We present a method for studying the proximity effect and the density structure around redshift z=2-3 quasars. It is based on the probability distribution of Lyman-alpha pixel optical depths and its evolution with redshift. We validate the method using mock spectra obtained from hydrodynamical simulations, and then apply it to a sample of 12 bright quasars at redshifts 2-3 observed with UVES at the VLT-UT2 Kueyen ESO telescope. These quasars do not show signatures of associated absorption and have a mean monochromatic luminosity of 5.4 10^{31} h^{-2} erg/s/Hz at the Lyman limit. The observed distribution of optical depth within 10 Mpc/h from the QSO is statistically different from that measured in the general intergalactic medium at the same redshift. Such a change will result from the combined effects of the increase in photoionisation rate above the mean UV-background due to the extra ionizing photons from the quasar radiation (proximity effect), and the higher density of the IGM if the quasars reside in overdense regions (as expected from biased galaxy formation). The first factor decreases the optical depth whereas the second one increases the optical depth, but our measurement cannot distinguish a high background from a low overdensity. An overdensity of the order of a few is required if we use the amplitude of the UV-background inferred from the mean Lyman-$\alpha$ opacity. If no overdensity is present, then we require the UV-background to be higher, and consistent with the existing measurements based on standard analysis of the proximity effect., Comment: 15 pages. 13 figures. revised version

The expansion of the Universe leads to a rapid drop in the mean hydrogen Ly-alpha opacity tau of the intergalactic medium (IGM), tau (1+z)^{3.8}, between redshifts 4 and 3. Measurements of the temperature evolution of the IGM and of the He II opacity both suggest that He II reionizes in this redshift range. We use hydrodynamical simulations to show that the temperature increase associated with He II reionization leads to a relatively sudden decrease in tau around the reionization epoch of ~ 10 per cent. We find clear evidence for such a feature in the evolution of tau determined from a sample of ~ 1100 quasars obtained from the SDSS. He II reionization starts at redshift ~ 3.4, and lasts for (Delta z) ~ 0.4. The increase in the IGM temperature also explains the widths of hydrogen absorption lines as measured in high-resolution spectra., Comment: 4 pages, 3 figures. Minor changes, to be published in ApJL

We study the ability of PINOCCHIO (PINpointing Orbit-Crossing Collapsed HIerarchical Objects) to predict the merging histories of dark matter (DM) haloes, comparing the PINOCCHIO predictions with the results of two large N-body simulations run from the same set of initial conditions. We focus our attention on quantities most relevant to galaxy formation and large-scale structure studies. PINOCCHIO is able to predict the statistics of merger trees with a typical accuracy of 20 per cent. Its validity extends to higher-order moments of the distribution of progenitors. The agreement is valid also at the object-by-object level, with 70-90 per cent of the progenitors cleanly recognised when the parent halo is cleanly recognised itself. Predictions are presented also for quantities that are usually not reproduced by semi-analytic codes, such as the two-point correlation function of the progenitors of massive haloes and the distribution of initial orbital parameters of merging haloes. For the accuracy of the prediction and for the facility with which merger histories are produced, PINOCCHIO provides a means to generate catalogues of DM haloes which is extremely competitive to large-scale N-body simulations, making it a suitable tool for galaxy formation and large-scale structure studies., Comment: 12 pages, 10 eps figures included, submitted to MNRAS. PINOCCHIO is available at http://www.daut.univ.trieste.it/pinocchio

PINOCCHIO (PINpointing Orbit-Crossing Collapsed Hierarchical Objects) is a new algorithm for identifying dark matter halos in a given numerical realisation of the linear density field in a hierarchical universe (Monaco et al. 2001). It is shown that Lagrangian perturbation theory, and in particular its ellipsoidal truncation, is able to predict accurately the collapse, in the orbit-crossing sense, of generic mass elements. Some points that have undergone orbit crossing are assigned to the network of filaments and sheets that connects the halos; it is demonstrated that this network resembles closely that found in N-body simulations. The code generates a catalogue of dark matter halos with known mass, position, velocity, merging history and angular momentum. It is shown that the predictions of the code are very accurate when compared with the results of large N-body simulations that cover a range of cosmological models, box sizes and numerical resolutions. The mass function is recovered with an accuracy of better than 10 per cent in number density for halos with at least 30-50 particles. A similar accuracy is reached in the estimate of the correlation length r_0. The good agreement is still valid on the object-by-object level, with 70-100 per cent of the objects with more than 50 particles in the simulations also identified by our algorithm. For these objects the masses are recovered with an error of 20-40 per cent, and positions and velocities with a root mean square error of ~1-2 Mpc (0.5-2 grid lengths) and ~100 km/s, respectively. The recovery of the angular momentum of halos is considerably noisier and accuracy at the statistical level is achieved only by introducing free parameters. The algorithm requires negligible computer time as compared with performing a numerical $N$-body simulation., Comment: 24 pages, 15 eps figures included. Revised to match the version accepted by MNRAS. New appendix added on N-body estimate of halo spins. PINOCCHIO is available at http://www.daut.univ.trieste.it/pinocchio

We present a new algorithm (PINOCCHIO, PINpointing Orbit-Crossing Collapsed HIerarchical objects) to predict accurately the formation and evolution of individual dark matter haloes in a given realization of an initial linear density field. Compared with the halo population formed in a large (360^3 particles) collisionless simulation of a CDM universe, our method is able to predict to better than 10 per cent statistical quantities such as the mass function, two-point correlation function and progenitor mass function of the haloes. Masses of individual haloes are estimated accurately as well, with errors typically of order 30 per cent in the mass range well resolved by the numerical simulation. These results show that the hierarchical formation of dark matter haloes can be accurately predicted using local approximations to the dynamics when the correlations in the initial density field are properly taken into account. The approach allows one to automatically generate a large ensemble of accurate merging histories of haloes with complete knowledge of their spatial distribution. The construction of the full merger tree for a 256^3 realisation requires a few hours of CPU-time on a personal computer, orders of magnitude faster than the corresponding N-body simulation would take, and without needing any extensive post-processing. The technique can be efficiently used, for instance, for generating the input for galaxy formation modeling., Comment: 9 pages, 6 eps figures included. ApJ, in press (1 January 2002) PINOCCHIO is available at http://www.daut.univ.trieste.it/pinocchio

Recent numerical simulations have lead to a paradigm shift in our understanding of the intergalactic medium, and the loss of a physical justification for Voigt profile fitting of the Lyman-alpha forest. Many individual lines seen in simulated spectra have significant departures from the Voigt profile, yet could be well fitted by a blend of two or more such lines. We discuss the expected effect on the line profiles due to ongoing gravitational structure formation and Hubble expansion. We develop a method to detect departures from Voigt profiles of the absorption lines in a statistical way and apply this method to simulated Lyman-alpha forest spectra, confirming that the profiles seen do statistically differ from Voigt profiles., Comment: Accepted for publication in ApJL. 10 pages, 3 figures

The low-redshift evolution of the intergalactic medium is investigated using hydrodynamic cosmological simulations. The assumed cosmological model is a critical density cold dark matter universe. The imposed uniform background of ionizing radiation has the amplitude, shape and redshift evolution as computed from the observed quasar luminosity function by Haardt & Madau. We have analysed simulated Lyman-alpha spectra using Voigt-profile fitting, mimicking the procedure with which quasar spectra are analysed. Our simulations reproduce the observed evolution of the number of Lyman-alpha absorption lines over the whole observed interval of z=0.5 to z=4. In particular, our simulations show that the decrease in the rate of evolution of Lyman-alpha absorption lines at z< 2, as observed by Hubble Space Telescope, can be explained by the steep decline in the photo-ionizing background resulting from the rapid decline in quasar numbers at low redshift., Comment: 4 pages, latex (mn.sty), 3 figures, submitted to MNRAS

A P3M (Particle-Particle, Particle-Mesh) algorithm to compute the gravitational force on a set of particles is described. The gravitational force is computed using Fast Fourier Transforms. This leads to an incorrect force when the distance between two particles is of the order of a grid cell. This incorrect force is subtracted exactly from all particles in parallel using convolution with the appropriate Green's function in real space in a time of order $N_T$, irrespective of the degree of clustering of particles. Next, the correct $1/r^2$ force is added for all neighbouring particles in parallel, leading to an accurate algorithm which runs efficiently on a highly parallel computer. A full force calculation for $128k$ particles on a $128^3$ grid in a mildly clustered situation requires approximately 196 seconds on a $8k$ Connection Machine 2 with 8MHz clock. This decreases to an estimated 9.8 seconds on a full-sized $64k$ CM200., Comment: plain TeX, 2 .ps Figs. available on request; SNS-ASTRO-93-0901