BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Adrianne Slyz

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

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-ph.IMastro-ph.HEastro-ph.SR

Catalog footprint

What is connected

47works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

47 published item(s)

preprint2022arXiv

Catalogues of voids as antihalos in the local Universe

A recently-proposed algorithm identifies voids in simulations as the regions associated with halos when the initial overdensity field is negated. We apply this method to the real Universe by running a suite of constrained simulations of the 2M++ volume with initial conditions inferred by the BORG algorithm, along with the corresponding inverted set. Our 101 inverted and uninverted simulations, spanning the BORG posterior, each identify ~150,000 "voids as antihalos" with mass exceeding $4.38\times10^{11} \: \mathrm{M_\odot}$ (100 particles) at $z=0$ in a full-sky sphere of radius 155 Mpc/h around the Milky Way. We calculate the size function, volume filling fraction, ellipticity, central density, specific angular momentum, clustering and stacked density profile of the voids, and cross-correlate them with those produced by VIDE on the same simulations. We make our antihalo and VIDE catalogues publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Impact of Radiation Feedback on the Formation of Globular Cluster Candidates during Cloud-Cloud Collisions

To understand the impact of radiation feedback during the formation of a globular cluster (GC), we simulate a head-on collision of two turbulent giant molecular clouds (GMCs). A series of idealized radiation-hydrodynamic simulations is performed, with and without stellar radiation or Type II supernovae. We find that a gravitationally bound, compact star cluster of mass $M_{\rm GC} \sim 10^5\,M_\odot$ forms within $\approx 3\,{\rm Myr}$ when two GMCs with mass $M_{\rm GMC}=3.6\times10^5\,M_\odot$ collide. The GC candidate does not form during a single collapsing event but emerges due to the mergers of local dense gas clumps and gas accretion. The momentum transfer due to the absorption of the ionizing radiation is the dominant feedback process that suppresses the gas collapse and photoionization becomes efficient once a sufficient number of stars form. The cluster mass is larger by a factor of $\sim 2$ when the radiation feedback is neglected, and the difference is slightly more pronounced (16\%) when extreme $\rm Lyα$ feedback is considered in the fiducial run. In the simulations with radiation feedback, supernovae explode after the star-forming clouds are dispersed, and their metal ejecta are not instantaneously recycled to form stars.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

MgII in the JWST Era: a Probe of Lyman Continuum Escape?

Limited constraints on the evolution of the Lyman Continuum (LyC) escape fraction represent one of the primary uncertainties in the theoretical determination of the reionization history. Due to the intervening intergalactic medium (IGM), the possibility of observing LyC photons directly in the epoch of reionization is highly unlikely. For this reason, multiple indirect probes of LyC escape have been identified, some of which are used to identify low-redshift LyC leakers (e.g. O32), while others are primarily useful at $z>6$ (e.g. [OIII]/[CII] far infrared emission). The flux ratio of the resonant MgII doublet emission at 2796$\dot{\rm A}$ and 2803$\dot{\rm A}$ as well as the MgII optical depth have recently been proposed as ideal diagnostics of LyC leakage that can be employed at $z>6$ with JWST. Using state-of-the-art cosmological radiation hydrodynamics simulations post-processed with CLOUDY and resonant-line radiative transfer, we test whether MgII is indeed a useful probe of LyC leakage. Our simulations indicate that the majority of bright, star-forming galaxies with high LyC escape fractions are expected to be MgII emitters rather than absorbers at $z=6$. However, we find that the MgII doublet flux ratio is a more sensitive indicator of dust rather than neutral hydrogen, limiting its use as a LyC leakage indicator to only galaxies in the optically thin regime. Given its resonant nature, we show that MgII will be an exciting probe of the complex kinematics in high-redshift galaxies in upcoming JWST observations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the Viability of Determining Galaxy Properties from Observations I: Star Formation Rates and Kinematics

We explore how observations relate to the physical properties of the emitting galaxies by post-processing a pair of merging $z\sim2$ galaxies from the cosmological, hydrodynamical simulation NewHorizon using LCARS (Light from Cloudy Added to RAMSES) to encode the physical properties of the simulated galaxy into H$α$ emission line. By carrying out mock observations and analysis on these data cubes we ascertain which physical properties of the galaxy will be recoverable with the HARMONI spectrograph on the European Extremely Large Telescope (ELT). We are able to estimate the galaxy's star formation rate and dynamical mass to a reasonable degree of accuracy, with values within a factor of $1.81$ and $1.38$ of the true value. The kinematic structure of the galaxy is also recovered in mock observations. Furthermore, we are able to recover radial profiles of the velocity dispersion and are therefore able to calculate how the dynamical ratio varies as a function of distance from the galaxy centre. Finally, we show that when calculated on galaxy scales the dynamical ratio does not always provide a reliable measure of a galaxy's stability against gravity or act as an indicator of a minor merger.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Simulating Jellyfish Galaxies: A Case Study for a Gas-Rich Dwarf Galaxy

We investigate the formation of jellyfish galaxies using radiation-hydrodynamic simulations of gas-rich dwarf galaxies with a multi-phase interstellar medium (ISM). We find that the ram-pressure-stripped (RPS) ISM is the dominant source of molecular clumps in the near wake within 10 kpc from the galactic plane, while in-situ formation is the major channel for dense gas in the distant tail of the gas-rich galaxy. Only 20% of the molecular clumps in the near wake originate from the intracluster medium (ICM); however, the fraction reaches 50% in the clumps located at $80\,{\rm kpc}$ from the galactic center since the cooling time of the RPS gas tends to be short due to the ISM--ICM mixing ($\lesssim$ 10 Myr). The tail region exhibits a star formation rate of $0.001-0.01\,{\rm M_{\odot}\,yr^{-1}}$, and most of the tail stars are born in the stripped wake within 10 kpc from the galactic plane. These stars induce bright H$α$ blobs in the tail, while H$α$ tails fainter than $6\times10^{38}\,{\rm erg\,s^{-1}\,kpc^{-2}}$ are mostly formed via collisional radiation and heating due to mixing. We also find that the stripped tails have intermediate X-ray to H$α$ surface brightness ratios (1.5$\lesssim F_{\rm X}/F_{\rm Hα}\lesssim$20), compared to the ISM ($\lesssim$1.5) or pure ICM ($\gg$20). Our results suggest that jellyfish features emerge when the ISM from gas-rich galaxies is stripped by strong ram pressure, mixes with the ICM, and enhances the cooling in the tail.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Nature of High ${\rm [OIII]}_{\rm 88μm}$/${\rm [CII]}_{\rm 158μm}$ Galaxies in the Epoch of Reionization: Low Carbon Abundance and a Top-Heavy IMF?

ALMA observations of $z>6$ galaxies have revealed abnormally high [OIII]$_{\rm 88μm}$/[CII]$_{\rm 158μm}$ ratios and [CII]$_{\rm 158μm}$ deficits compared to local galaxies. Numerous solutions have been proposed including differences in C and O abundance ratios, observational bias, and differences in ISM properties. In order to elucidate the underlying physics that drives this high-redshift phenomenon, we employ SPHINX$^{20}$, a state-of-the-art, cosmological radiation-hydrodynamics simulation, that resolves detailed ISM properties of thousands of galaxies in the epoch of reionization. We find that the observed $z>6$ [OIII]-SFR and [CII]-SFR relations can only be reproduced when the C/O abundance ratio is $\sim8\times$ lower than Solar and the total metal production is $\sim4\times$ higher than that of a Kroupa IMF. This implies that high-redshift galaxies are potentially primarily enriched by low-metallicity core-collapse supernovae with a more top-heavy IMF. As AGB stars and type-Ia supernova begin to contribute to the galaxy metallicity, both the [CII]-SFR and [CII] luminosity functions are predicted to converge to observed values at $z\sim4.5$. While we demonstrate that ionisation parameter, LyC escape fraction, ISM gas density, and CMB attenuation all drive galaxies towards higher [OIII]/[CII], observed values at $z>6$ can only be reproduced with substantially lower C/O abundances compared to Solar. The combination of [CII] and [OIII] can be used to predict the values of ionisation parameter, ISM gas density, and LyC escape fraction and we provide estimates of these quantities for nine observed $z>6$ galaxies. Finally, we demonstrate that [OI]$_{\rm 63μm}$ can be used as a replacement for [CII] in high-redshift galaxies where [CII] is unobserved and argue that more observation time should be used to target [OI] at $z>6$. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Obelisk simulation: galaxies contribute more than AGN to HI reionization of protoclusters

We present the Obelisk project, a cosmological radiation-hydrodynamics simulation following the assembly and reionization of a protocluster progenitor during the first two billions of years from the big bang, down to $z = 3.5$. The simulation resolves haloes down to the atomic cooling limit, and tracks the contribution of different sources of ionization: stars, active galactic nuclei, and collisions. The Obelisk project is designed specifically to study the coevolution of high redshift galaxies and quasars in an environment favouring black hole growth. In this paper, we establish the relative contribution of these two sources of radiation to reionization and their respective role in establishing and maintaining the high redshift ionizing background. Our volume is typical of an overdense region of the Universe and displays star formation rate and black hole accretion rate densities similar to high redshift protoclusters. We find that hydrogen reionization happens inside-out and is completed by $z \sim 6$ in our overdensity, and is predominantly driven by galaxies, while accreting black holes only play a role at $z \sim 4$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Towards convergence of turbulent dynamo amplification in cosmological simulations of galaxies

Our understanding of the process through which magnetic fields reached their observed strengths in present-day galaxies remains incomplete. One of the advocated solutions is a turbulent dynamo mechanism that rapidly amplifies weak magnetic field seeds to the order of ${\sim}μ$G. However, simulating the turbulent dynamo is a very challenging computational task due to the demanding span of spatial scales and the complexity of the required numerical methods. In particular, turbulent velocity and magnetic fields are extremely sensitive to the spatial discretisation of simulated domains. To explore how refinement schemes affect galactic turbulence and amplification of magnetic fields in cosmological simulations, we compare two refinement strategies. A traditional quasi-Lagrangian adaptive mesh refinement approach focusing spatial resolution on dense regions, and a new refinement method that resolves the entire galaxy with a high resolution quasi-uniform grid. Our new refinement strategy yields much faster magnetic energy amplification than the quasi-Lagrangian method, which is also significantly greater than the adiabatic compressional estimate indicating that the extra amplification is produced through stretching of magnetic field lines. Furthermore, with our new refinement the magnetic energy growth factor scales with resolution following $\propto \Dres^{-1/2}$, in much better agreement with small-scale turbulent box simulations. Finally, we find evidence suggesting most magnetic amplification in our simulated galaxies occurs in the warm phase of their interstellar medium, which has a better developed turbulent field with our new refinement strategy.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Momentum deposition of supernovae with cosmic rays

The cataclysmic explosions of massive stars as supernovae are one of the key ingredients of galaxy formation. However, their evolution is not well understood in the presence of magnetic fields or cosmic rays (CRs). We study the expansion of individual supernova remnants (SNRs) using our suite of 3D hydrodynamical (HD), magnetohydrodynamical (MHD) and CRMHD simulations generated using RAMSES. We explore multiple ambient densities, magnetic fields and fractions of supernova energy deposited as CRs ($χ_{\rm CR}$), accounting for cosmic ray anisotropic diffusion and streaming. All our runs have comparable evolutions until the end of the Sedov-Taylor phase. However, our CRMHD simulations experience an additional CR pressure-driven snowplough phase once the CR energy dominates inside the SNR. We present a model for the final momentum deposited by supernovae that captures this new phase: $p_{\rm SNR} = 2.87\times 10^{5} (χ_{\text{CR}} + 1)^{4.82}\left(\frac{n}{\text{cm}^{-3}}\right)^{-0.196} M_{\odot}$ km s$^{-1}$. Assuming a 10% fraction of SN energy in CRs leads to a 50% boost of the final momentum, with our model predicting even higher impacts at lower ambient densities. The anisotropic diffusion of CRs assuming an initially uniform magnetic field leads to extended gas and cosmic ray outflows escaping from the supernova poles. We also study a tangled initial configuration of the magnetic field, resulting instead in a quasi-isotropic diffusion of CRs and earlier momentum deposition. Finally, synthetic synchrotron observations of our simulations using the POLARIS code show that the local magnetic field configuration in the interstellar medium modifies the overall radio emission morphology and polarisation.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Predicting the Observability of Population III Stars with ELT-HARMONI via the Helium $1640{\rmÅ}$ emission line

Population III (Pop. III) stars, as of yet, have not been detected, however as we move into the era of extremely large telescopes this is likely to change. One likely tracer for Pop. III stars is the HeII$\lambda1640$ emission line, which will be detectable by the HARMONI spectrograph on the European Extremely Large Telescope (ELT) over a broad range of redshifts ($2\leq z\leq14$). By post-processing galaxies from the cosmological, AMR-hydrodynamical simulation NewHorizon with theoretical spectral energy distributions (SED) for Pop. III stars and radiative transfer (i.e. the Yggdrasil Models and CLOUDY look-up tables respectively) we are able to compute the flux of HeII$\lambda1640$ for individual galaxies. From mock 10 hour observations of these galaxies we show that HARMONI will be able to detect Pop. III stars in galaxies up to $z\sim10$ provided Pop. III stars have a top heavy Initial Mass Function (IMF). Furthermore, we find that should Pop. III stars instead have an IMF similar to those of the Pop. I stars, the HeII$\lambda1640$ line would only be observable for galaxies with Pop. III stellar masses in excess of $10^{7}{\,M}_\odot\,{\rm yr}^{-1}$, average stellar age $<1{\rm Myr}$ at $z=4$. Finally, we are able to determine the minimal intrinsic flux required for HARMONI to detect Pop. III stars in a galaxy up to $z=10$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Rivers of Gas I.: Unveiling The Properties of High Redshift Filaments

At high redshift, the cosmic web is widely expected to have a significant impact on the morphologies, dynamics and star formation rates of the galaxies embedded within it, underscoring the need for a comprehensive study of the properties of such a filamentary network. With this goal in mind, we perform an analysis of high-$z$ gas and dark matter (DM) filaments around a Milky Way-like progenitor simulated with the {\sc ramses} adaptive mesh refinement (AMR) code from cosmic scales ($\sim$1 Mpc) down to the virial radius of its DM halo host ($\sim$20 kpc at $z=4$). Radial density profiles of both gas and DM filaments are found to have the same functional form, namely a plummer-like profile modified to take into account the wall within which these filaments are embedded. Measurements of the typical filament core radius $r_0$ from the simulation are consistent with that of isothermal cylinders in hydrostatic equilibrium. Such an analytic model also predicts a redshift evolution for the core radius of filaments in fair agreement with the measured value for DM $(r_0 \propto (1+z)^{-3.18\pm 0.28})$. Gas filament cores grow as $(r_0 \propto (1+z)^{-2.72\pm 0.26})$. In both gas and DM, temperature and vorticity sharply drop at the edge of filaments, providing an excellent way to constrain the outer filament radius. When feedback is included the gas temperature and vorticity fields are strongly perturbed, hindering such a measurement in the vicinity of the galaxy. However, the core radius of the filaments as measured from the gas density field is largely unaffected by feedback, and the median central density is only reduced by about 20%.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Detecting the Cosmic Web: Lyα Emission from Simulated Filaments at z=3

The standard cosmological model ($Λ$CDM) predicts the existence of the cosmic web: a distribution of matter into sheets and filaments connecting massive halos. However, observational evidence has been elusive due to the low surface brightness of the filaments. Recent deep MUSE/VLT data and upcoming observations offer a promising avenue for Ly$α$ detection, motivating the development of modern theoretical predictions. We use hydrodynamical cosmological simulations run with the AREPO code to investigate the potential detectability of large-scale filaments, excluding contributions from the halos embedded in them. We focus on filaments connecting massive ($M_{200c}\sim(1-3)\times10^{12} M_\odot$) halos at z=3, and compare different simulation resolutions, feedback levels, and mock-image pixel sizes. We find increasing simulation resolution does not substantially improve detectability notwithstanding the intrinsic enhancement of internal filament structure. By contrast, for a MUSE integration of 31 hours, including feedback increases the detectable area by a factor of $\simeq$5.5 on average compared with simulations without feedback, implying that even the non-bound components of the filaments have substantial sensitivity to feedback. Degrading the image resolution from the native MUSE scale of (0.2")$^2$ per pixel to (5.3")$^2$ apertures has the strongest effect, increasing the detectable area by a median factor of $\simeq$200 and is most effective when the size of the pixel roughly matches the width of the filament. Finally, we find the majority of Ly$α$ emission is due to electron impact collisional excitations, as opposed to radiative recombination.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

How primordial magnetic fields shrink galaxies

As one of the prime contributors to the interstellar medium energy budget, magnetic fields naturally play a part in shaping the evolution of galaxies. Galactic magnetic fields can originate from strong primordial magnetic fields provided these latter remain below current observational upper limits. To understand how such magnetic fields would affect the global morphological and dynamical properties of galaxies, we use a suite of high-resolution constrained transport magneto-hydrodynamic cosmological zoom simulations where we vary the initial magnetic field strength and configuration along with the prescription for stellar feedback. We find that strong primordial magnetic fields delay the onset of star formation and drain the rotational support of the galaxy, diminishing the radial size of the galactic disk and driving a higher amount of gas towards the centre. This is also reflected in mock UVJ observations by an increase in the light profile concentration of the galaxy. We explore the possible mechanisms behind such a reduction in angular momentum, focusing on magnetic braking. Finally, noticing that the effects of primordial magnetic fields are amplified in the presence of stellar feedback, we briefly discuss whether the changes we measure would also be expected for galactic magnetic fields of non-primordial origin.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

How to Quench a Dwarf Galaxy: The Impact of Inhomogeneous Reionization on Dwarf Galaxies and Cosmic Filaments

We use the SPHINX suite of high-resolution cosmological radiation hydrodynamics simulations to study how spatially and temporally inhomogeneous reionization impacts the baryonic content of dwarf galaxies and cosmic filaments. The SPHINX simulations simultaneously model an inhomogeneous reionization, follow the escape of ionising radiation from thousands of galaxies, and resolve haloes well below the atomic cooling threshold. This makes them an ideal tool for examining how reionization impacts star formation and the gas content of dwarf galaxies. We compare simulations with and without stellar radiation to isolate the effects of radiation feedback from that of supernova, cosmic expansion, and numerical resolution. We find that the gas content of cosmic filaments can be reduced by more than 80% following reionization. The gas inflow rates into haloes with $M_{vir}<10^8M_{\odot}$ are strongly affected and are reduced by more than an order of magnitude compared to the simulation without reionization. A significant increase in gas outflow rates is found for halo masses $M_{vir}<7\times10^7M_{\odot}$. Our simulations show that inflow suppression, rather than photoevaporation, is the dominant mechanism by which the baryonic content of high-redshift dwarf galaxies is regulated. At fixed redshift and halo mass, there is a large scatter in the halo baryon fractions that is entirely dictated by the timing of reionization in the local region surrounding a halo which can change by $Δz>3$ at fixed mass. Finally, although the gas content of high-redshift dwarf galaxies is significantly impacted by reionization, we find that most haloes with $M_{vir}<10^8M_{\odot}$ can remain self-shielded and form stars long after reionization, until their local gas reservoir is depleted, suggesting that local group dwarf galaxies do not necessarily exhibit star formation histories that peak prior to $z=6$...

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

New Methods for Identifying Lyman Continuum Leakers and Reionization-Epoch Analogues

Identifying low-redshift galaxies that emit Lyman Continuum radiation (LyC leakers) is one of the primary, indirect methods of studying galaxy formation in the epoch of reionization. However, not only has it proved challenging to identify such systems, it also remains uncertain whether the low-redshift LyC leakers are truly "analogues" of the sources that reionized the Universe. Here, we use high-resolution cosmological radiation hydrodynamics simulations to examine whether simulated galaxies in the epoch of reionization share similar emission line properties to observed LyC leakers at $z\sim3$ and $z\sim0$. We find that the simulated galaxies with high LyC escape fractions ($f_{\rm esc}$) often exhibit high O32 and populate the same regions of the R23-O32 plane as $z\sim3$ LyC leakers. However, we show that viewing angle, metallicity, and ionisation parameter can all impact where a galaxy resides on the O32-$f_{\rm esc}$ plane. Based on emission line diagnostics and how they correlate with $f_{\rm esc}$, lower-metallicity LyC leakers at $z\sim3$ appear to be good analogues of reionization-era galaxies. In contrast, local [SII]-deficient galaxies do not overlap with the simulated high-redshift LyC leakers on the SII-BPT diagram; however, this diagnostic may still be useful for identifying leakers. We use our simulated galaxies to develop multiple new diagnostics to identify LyC leakers using IR and nebular emission lines. We show that our model using only [CII]$_{\rm 158μm}$ and [OIII]$_{\rm 88μm}$ can identify potential leakers from non-leakers from the local Dwarf Galaxy Survey. Finally, we apply this diagnostic to known high-redshift galaxies and find that MACS1149_JD1 at $z=9.1$ is the most likely galaxy to be actively contributing to the reionization of the Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The impact of AGN feedback on galaxy intrinsic alignments in the Horizon simulations

The intrinsic correlations of galaxy shapes and orientations across the large-scale structure of the Universe are a known contaminant to weak gravitational lensing. They are known to be dependent on galaxy properties, such as their mass and morphologies. The complex interplay between alignments and the physical processes that drive galaxy evolution remains vastly unexplored. We assess the sensitivity of intrinsic alignments (shapes and angular momenta) to Active Galactic Nuclei -AGN- feedback by comparing galaxy alignment in twin runs of the cosmological hydrodynamical Horizon simulation, which do and do not include AGN feedback respectively. We measure intrinsic alignments in three dimensions and in projection at z=0 and z=1. We find that the projected alignment signal of all galaxies with resolved shapes with respect to the density field in the simulation is robust to AGN feedback, thus giving similar predictions for contamination to weak lensing. The relative alignment of galaxy shapes around galaxy positions is however significantly impacted, especially when considering high-mass ellipsoids. Using a sample of galaxy "twins" across simulations, we determine that AGN changes both the galaxy selection and their actual alignments. Finally, we measure the alignments of angular momenta of galaxies with their nearest filament. Overall, these are more significant in the presence of AGN as a result of the higher abundance of massive pressure-supported galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The impact of the connectivity of the cosmic web on the physical properties of galaxies at its nodes

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).

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

When galaxies align: intrinsic alignments of the progenitors of elliptical galaxies in the Horizon-AGN simulation

Elliptical galaxies today appear aligned with the large-scale structure of the Universe, but it is still an open question when they acquire this alignment. Observational data is currently insufficient to provide constraints on the time evolution of intrinsic alignments, and hence existing models range from assuming that galaxies gain some primordial alignment at formation, to suggesting that they react instantaneously to tidal interactions with the large-scale structure. Using the cosmological hydrodynamical simulation Horizon-AGN, we measure the relative alignments between the major axes of galaxies and eigenvectors of the tidal field as a function of redshift. We focus on constraining the time evolution of the alignment of the main progenitors of massive $z=0$ elliptical galaxies, the main weak lensing contaminant at low redshift. We show that this population, which at $z=0$ has a stellar mass above $10^{10.4}$ M$_\odot$, transitions from having no alignment with the tidal field at $z=3$, to a significant alignment by $z=1$. From $z=0.5$ they preserve their alignment at an approximately constant level until $z=0$. We find a mass-dependence of the alignment signal of elliptical progenitors, whereby ellipticals that are less massive today ($10^{10.4}<M/{\rm M}_\odot<10^{10.7}$) do not become aligned till later redshifts ($z<2$), compared to more massive counterparts. We also present an extended study of progenitor alignments in the parameter space of stellar mass and galaxy dynamics, the impact of shape definition and tidal field smoothing.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Bursty star formation feedback and cooling outflows

We study how outflows of gas launched from a central galaxy undergoing repeated starbursts propagate through the circumgalactic medium (CGM), using the simulation code RAMSES. We assume that the outflow from the disk can be modelled as a rapidly moving bubble of hot gas at $\mathrm{\sim1\;kpc}$ above disk, then ask what happens as it moves out further into the halo around the galaxy on $\mathrm{\sim 100\;kpc}$ scales. To do this we run 60 two-dimensional simulations scanning over parameters of the outflow. Each of these is repeated with and without radiative cooling, assuming a primordial gas composition to give a lower bound on the importance of cooling. In a large fraction of radiative-cooling cases we are able to form rapidly outflowing cool gas from in situ cooling of the flow. We show that the amount of cool gas formed depends strongly on the 'burstiness' of energy injection; sharper, stronger bursts typically lead to a larger fraction of cool gas forming in the outflow. The abundance ratio of ions in the CGM may therefore change in response to the detailed historical pattern of star formation. For instance, outflows generated by star formation with short, intense bursts contain up to 60 per cent of their gas mass at temperatures $<5 \times 10^4\,\mathrm{K}$; for near-continuous star formation the figure is $\lesssim$ 5 per cent. Further study of cosmological simulations, and of idealised simulations with e.g., metal-cooling, magnetic fields and/or thermal conduction, will help to understand the precise signature of bursty outflows on observed ion abundances.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Comparing Simulations of AGN Feedback

We perform adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) cosmological zoom simulations of a region around a forming galaxy cluster, comparing the ability of the methods to handle successively more complex baryonic physics. In the simplest, non-radiative case, the two methods are in good agreement with each other, but the SPH simulations generate central cores with slightly lower entropies and virial shocks at slightly larger radii, consistent with what has been seen in previous studies. The inclusion of radiative cooling, star formation, and stellar feedback leads to much larger differences between the two methods. Most dramatically, at z=5, rapid cooling in the AMR case moves the accretion shock well within the virial radius, while this shock remains near the virial radius in the SPH case, due to excess heating, coupled with poorer capturing of the shock width. On the other hand, the addition of feedback from active galactic nuclei (AGN) to the simulations results in much better agreement between the methods. In this case both simulations display halo gas entropies of 100 keV cm^2, similar decrements in the star-formation rate, and a drop in the halo baryon content of roughly 30%. This is consistent with AGN growth being self-regulated, regardless of the numerical method. However, the simulations with AGN feedback continue to differ in aspects that are not self-regulated, such that in SPH a larger volume of gas is impacted by feedback, and the cluster still has a lower entropy central core.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Redshift and luminosity evolution of the intrinsic alignments of galaxies in Horizon-AGN

Intrinsic galaxy shape and angular momentum alignments can arise in cosmological large-scale structure due to tidal interactions or galaxy formation processes. Cosmological hydrodynamical simulations have recently come of age as a tool to study these alignments and their contamination to weak gravitational lensing. We probe the redshift and luminosity evolution of intrinsic alignments in Horizon-AGN between z=0 and z=3 for galaxies with an r-band absolute magnitude of <-20. Alignments transition from being radial at low redshifts and high luminosities, dominated by the contribution of ellipticals, to being tangential at high redshift and low luminosities, where discs dominate the signal. This cannot be explained by the evolution of the fraction of ellipticals and discs alone: intrinsic evolution in the amplitude of alignments is necessary. The alignment amplitude of elliptical galaxies alone is smaller in amplitude by a factor of ~2, but has similar luminosity and redshift evolution as in current observations and in the nonlinear tidal alignment model at projected separations of > 1 Mpc. Alignments of discs are null in projection and consistent with current low redshift observations. The combination of the two populations yields an overall amplitude a factor of ~4 lower than observed alignments of luminous red galaxies with a steeper luminosity dependence. The restriction on accurate galaxy shapes implies that the galaxy population in the simulation is complete only to an r-band absolute magnitude of <-20. Higher resolution simulations will be necessary to avoid extrapolation of the intrinsic alignment predictions to the range of luminosities probed by future surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Detailed Study of Feedback from a Massive Star

We present numerical simulations of a 15 solar mass star in a suite of idealised environments in order to quantify the amount of energy transmitted to the interstellar medium (ISM). We include models of stellar winds, UV photoionisation and the subsequent supernova based on theoretical models and observations of stellar evolution. The system is simulated in 3D using RAMSES-RT, an Adaptive Mesh Refinement Radiation Hydrodynamics code. We find that stellar winds have a negligible impact on the system owing to their relatively low luminosity compared to the other processes. The main impact of photoionisation is to reduce the density of the medium into which the supernova explodes, reducing the rate of radiative cooling of the subsequent supernova. Finally, we present a grid of models quantifying the energy and momentum of the system that can be used to motivate simulations of feedback in the ISM unable to fully resolve the processes discussed in this work.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Black hole evolution: I. Supernova-regulated black hole growth

The growth of a supermassive black hole (BH) is determined by how much gas the host galaxy is able to feed it, which in turn is controlled by the cosmic environment, through galaxy mergers and accretion of cosmic flows that time how galaxies obtain their gas, but also by internal processes in the galaxy, such as star formation and feedback from stars and the BH itself. In this paper, we study the growth of a 10^12 Msun halo at z=2, which is the progenitor of al group of galaxies at z=0, and of its central BH by means of a high-resolution zoomed cosmological simulation, the Seth simulation. We study the evolution of the BH driven by the accretion of cold gas in the galaxy, and explore the efficiency of the feedback from supernovae (SNe). For a relatively inefficient energy input from SNe, the BH grows at the Eddington rate from early times, and reaches self-regulation once it is massive enough. We find that at early cosmic times z>3.5, efficient feedback from SNe forbids the formation of a settled disc as well as the accumulation of dense cold gas in the vicinity of the BH and starves the central compact object. As the galaxy and its halo accumulate mass, they become able to confine the nuclear inflows provided by major mergers and the BH grows at a sustained near-to-Eddington accretion rate. We argue that this mechanism should be ubiquitous amongst low-mass galaxies, corresponding to galaxies with a stellar mass below <10^9 Msun in our simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Intrinsic alignment of simulated galaxies in the cosmic web: implications for weak lensing surveys

The intrinsic alignment of galaxy shapes (by means of their angular momentum) and their cross-correlation with the surrounding dark matter tidal field are investigated using the 160 000, z=1.2 synthetic galaxies extracted from the high-resolution cosmological hydrodynamical simulation Horizon-AGN. One- and two-point statistics of the spin of the stellar component are measured as a function of mass and colour. For the low-mass galaxies, this spin is locally aligned with the tidal field `filamentary' direction while, for the high-mass galaxies, it is perpendicular to both filaments and walls. The bluest galaxies of our synthetic catalog are more strongly correlated with the surrounding tidal field than the reddest galaxies, and this correlation extends up to 10 Mpc/h comoving distance. We also report a correlation of the projected ellipticities of blue, intermediate mass galaxies on a similar scale at a level of 10^(-4) which could be a concern for cosmic shear measurements. We do not report any measurable intrinsic alignments of the reddest galaxies of our sample. This work is a first step toward the use of very realistic catalog of synthetic galaxies to evaluate the contamination of weak lensing measurement by the intrinsic galactic alignments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Intrinsic alignments of galaxies in the Horizon-AGN cosmological hydrodynamical simulation

The intrinsic alignments of galaxies are recognised as a contaminant to weak gravitational lensing measurements. In this work, we study the alignment of galaxy shapes and spins at low redshift ($z\sim 0.5$) in Horizon-AGN, an adaptive-mesh-refinement hydrodynamical cosmological simulation box of 100 Mpc/h a side with AGN feedback implementation. We find that spheroidal galaxies in the simulation show a tendency to be aligned radially towards over-densities in the dark matter density field and other spheroidals. This trend is in agreement with observations, but the amplitude of the signal depends strongly on how shapes are measured and how galaxies are selected in the simulation. Disc galaxies show a tendency to be oriented tangentially around spheroidals in three-dimensions. While this signal seems suppressed in projection, this does not guarantee that disc alignments can be safely ignored in future weak lensing surveys. The shape alignments of luminous galaxies in Horizon-AGN are in agreement with observations and other simulation works, but we find less alignment for lower luminosity populations. We also characterize the systematics of galaxy shapes in the simulation and show that they can be safely neglected when measuring the correlation of the density field and galaxy ellipticities.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Towards simulating star formation in turbulent high-z galaxies with mechanical supernova feedback

Feedback from supernovae is essential to understanding the self-regulation of star formation in galaxies. However, the efficacy of the process in a cosmological context remains unclear due to excessive radiative losses during the shock propagation. To better understand the impact of SN explosions on the evolution of galaxies, we perform a suite of high-resolution (12 pc), zoom-in cosmological simulations of a Milky Way-like galaxy at z=3 with adaptive mesh refinement. We find that SN explosions can efficiently regulate star formation, leading to the stellar mass and metallicity consistent with the observed mass-metallicity relation and stellar mass-halo mass relation at z~3. This is achieved by making three important changes to the classical feedback scheme: i) the different phases of SN blast waves are modelled directly by injecting radial momentum expected at each stage, ii) the realistic time delay of SNe, commencing at as early as 3 Myr, is required to disperse very dense gas before a runaway collapse sets in at the galaxy centre via mergers of gas clumps, and iii) a non-uniform density distribution of the ISM is taken into account below the computational grid scale for the cell in which SN explodes. The last condition is motivated by the fact that our simulations still do not resolve the detailed structure of a turbulent ISM in which the fast outflows can propagate along low-density channels. The simulated galaxy with the SN feedback model shows strong outflows, which carry approximately ten times larger mass than star formation rate, as well as smoothly rising circular velocity. Other feedback models that do not meet the three conditions form too many stars, producing a peaked rotation curve. Our results suggest that understanding the structure of the turbulent ISM may be crucial to assess the role of SN and other feedback processes in galaxy formation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Black hole evolution: II. Spinning black holes in a supernova-driven turbulent interstellar medium

Supermassive black holes (BH) accrete gas from their surroundings and coalesce with companions during galaxy mergers, and both processes change the BH mass and spin. By means of high-resolution hydrodynamical simulations of galaxies, either idealised or embedded within the cosmic web, we explore the effects of interstellar gas dynamics and external perturbations on BH spin evolution. All these physical quantities were evolved on-the-fly in a self-consistent manner. We use a `maximal' model to describe the turbulence induced by stellar feedback to highlight its impact on the angular momentum of the gas accreted by the BH. Periods of intense star formation are followed by phases where stellar feedback drives large-scale outflows and hot bubbles. We find that BH accretion is synchronised with star formation, as only when gas is cold and dense do both processes take place. During such periods, gas motion is dominated by consistent rotation. On the other hand, when stellar feedback becomes substantial, turbulent motion randomises gas angular momentum. However BH accretion is strongly suppressed in that case, as cold and dense gas is lacking. In our cosmological simulation, at very early times (z>6), the galactic disc has not yet settled and no preferred direction exists for the angular momentum of the accreted gas, so the BH spin remains low. As the gas settles into a disc (6>z>3), the BH spin then rapidly reaches its maximal value. At lower redshifts (z<3), even when galaxy mergers flip the direction of the angular momentum of the accreted gas, causing it to counter-rotate, the BH spin magnitude only decreases modestly and temporarily. Should this be a typical evolution scenario for BH, it potentially has dramatic consequences regarding their origin and assembly, as accretion on maximally spinning BH embedded in thin Shakura-Sunyaev disc is significantly reduced.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dancing in the dark: galactic properties trace spin swings along the cosmic web

A large-scale hydrodynamical cosmological simulation, Horizon-AGN, is used to investigate the alignment between the spin of galaxies and the cosmic filaments above redshift 1.2. The analysis of more than 150 000 galaxies per time step in the redshift range 1.2<z<1.8 with morphological diversity shows that the spin of low-mass blue galaxies is preferentially aligned with their neighbouring filaments, while high-mass red galaxies tend to have a perpendicular spin. The reorientation of the spin of massive galaxies is provided by galaxy mergers, which are significant in their mass build-up. We find that the stellar mass transition from alignment to misalignment happens around 3.10^10 M_sun. Galaxies form in the vorticity-rich neighbourhood of filaments, and migrate towards the nodes of the cosmic web as they convert their orbital angular momentum into spin. The signature of this process can be traced to the properties of galaxies, as measured relative to the cosmic web. We argue that a strong source of feedback such as active galactic nuclei is mandatory to quench in situ star formation in massive galaxies and promote various morphologies. It allows mergers to play their key role by reducing post-merger gas inflows and, therefore, keeping spins misaligned with cosmic filaments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Swirling around filaments: are large-scale structure vortices spinning up dark halos?

The kinematic analysis of dark matter and hydrodynamical simulations suggests that the vorticity in large-scale structure is mostly confined to, and predominantly aligned with their filaments, with an excess of probability of 20 per cent to have the angle between vorticity and filaments direction lower than 60 degrees relative to random orientations. The cross sections of these filaments are typically partitioned into four quadrants with opposite vorticity sign, arising from multiple flows, originating from neighbouring walls. The spins of halos embedded within these filaments are consistently aligned with this vorticity for any halo mass, with a stronger alignment for the most massive structures up to an excess of probability of 165 per cent. On large scales, adiabatic/cooling hydrodynamical simulations display the same vorticity in the gas as in the dark matter. The global geometry of the flow within the cosmic web is therefore qualitatively consistent with a spin acquisition for smaller halos induced by this large-scale coherence, as argued in Codis et al. (2012). In effect, secondary anisotropic infall (originating from the vortex-rich filament within which these lower-mass halos form) dominates the angular momentum budget of these halos. The transition mass from alignment to orthogonality is related to the size of a given multi-flow region with a given polarity. This transition may be reconciled with the standard tidal torque theory if the latter is augmented so as to account for the larger scale anisotropic environment of walls and filaments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Angular momentum transfer to a Milky Way disk at high redshift

An Adaptive Mesh Refinement cosmological resimulation is analyzed in order to test whether filamentary flows of cold gas are responsible for the build-up of angular momentum within a Milky Way like disk at z>=3. A set of algorithms is presented that takes advantage of the high spatial resolution of the simulation (12 pc) to identify: (i) the central gas disk and its plane of orientation; (ii) the complex individual filament trajectories that connect to the disk, and; (iii) the infalling satellites. The results show that two filaments at z>5.5, which later merge to form a single filament at z<4, drive the angular momentum and mass budget of the disk throughout its evolution, whereas luminous satellite mergers make negligible fractional contributions. Combined with the ubiquitous presence of such filaments in all large-scale cosmological simulations that include hydrodynamics, these findings provide strong quantitative evidence that the growth of thin disks in haloes with masses below 10^{12} M_{sun}, which host the vast majority of galaxies, is supported via inflowing streams of cold gas at intermediate and high redshifts.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Blowing cold flows away: the impact of early AGN activity on the formation of a brightest cluster galaxy progenitor

Supermassive black holes (BH) are powerful sources of energy that are already in place at very early epochs of the Universe (by z=6). Using hydrodynamical simulations of the formation of a massive M_vir=5 10^11 M_sun halo by z=6 (the most massive progenitor of a cluster of M_vir=2 10^15 M_sun at z=0), we evaluate the impact of Active Galactic Nuclei (AGN) on galaxy mass content, BH self-regulation, and gas distribution inside this massive halo. We find that SN feedback has a marginal influence on the stellar structure, and no influence on the mass distribution on large scales. In contrast, AGN feedback alone is able to significantly alter the stellar-bulge mass content by quenching star formation when the BH is self-regulating, and by depleting the cold gas reservoir in the centre of the galaxy. The growth of the BH proceeds first by a rapid Eddington-limited period fed by direct cold filamentary infall. When the energy delivered by the AGN is sufficiently large to unbind the cold gas of the bulge, the accretion of gas onto the BH is maintained both by smooth gas inflow and clump migration through the galactic disc triggered by merger-induced torques. The feedback from the AGN has also a severe consequence on the baryon mass content within the halo, producing large-scale hot superwinds, able to blow away some of the cold filamentary material from the centre and reduce the baryon fraction by more than 30 per cent within the halo's virial radius. Thus in the very young universe, AGN feedback is likely to be a key process, shaping the properties of the most massive galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Connecting the cosmic web to the spin of dark halos: implications for galaxy formation

We investigate the alignment of the spin of dark matter halos relative (i) to the surrounding large-scale filamentary structure, and (ii) to the tidal tensor eigenvectors using the Horizon 4pi dark matter simulation which resolves over 43 million dark matter halos at redshift zero. We detect a clear mass transition: the spin of dark matter halos above a critical mass tends to be perpendicular to the closest filament, and aligned with the intermediate axis of the tidal tensor, whereas the spin of low-mass halos is more likely to be aligned with the closest filament. Furthermore, this critical mass of 5 10^12 is redshift-dependent and scales as (1+z)^-2.5. We propose an interpretation of this signal in terms of large-scale cosmic flows. In this picture, most low-mass halos are formed through the winding of flows embedded in misaligned walls; hence they acquire a spin parallel to the axis of the resulting filaments forming at the intersection of these walls. On the other hand, more massive halos are typically the products of later mergers along such filaments, and thus they acquire a spin perpendicular to this direction when their orbital angular momentum is converted into spin. We show that this scenario is consistent with both the measured excess probabilities of alignment w.r.t. the eigen-directions of the tidal tensor, and halo merger histories. On a more qualitative level, it also seems compatible with 3D visualization of the structure of the cosmic web as traced by "smoothed" dark matter simulations or gas tracer particles. Finally, it provides extra support to the disc forming paradigm presented by Pichon et al (2011) as it extends it by characterizing the geometry of secondary infall at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constraining stellar assembly and AGN feedback at the peak epoch of star formation

We study stellar assembly and feedback from active galactic nuclei (AGN) around the epoch of peak star formation (1<z<2), by comparing hydrodynamic simulations to rest-frame UV-optical galaxy colours from the Wide Field Camera 3 (WFC3) Early-Release Science (ERS) Programme. Our Adaptive Mesh Refinement simulations include metal-dependent radiative cooling, star formation, kinetic outflows due to supernova explosions, and feedback from supermassive black holes. Our model assumes that when gas accretes onto black holes, a fraction of the energy is used to form either thermal winds or sub-relativistic momentum-imparting collimated jets, depending on the accretion rate. We find that the predicted rest-frame UV-optical colours of galaxies in the model that includes AGN feedback is in broad agreement with the observed colours of the WFC3 ERS sample at 1<z<2. The predicted number of massive galaxies also matches well with observations in this redshift range. However, the massive galaxies are predicted to show higher levels of residual star formation activity than the observational estimates, suggesting the need for further suppression of star formation without significantly altering the stellar mass function. We discuss possible improvements, involving faster stellar assembly through enhanced star formation during galaxy mergers while star formation at the peak epoch is still modulated by the AGN feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Feeding compact bulges and supermassive black holes with low angular-momentum cosmic gas at high redshift

We use cosmological hydrodynamical simulations to show that a significant fraction of the gas in high redshift rare massive halos falls nearly radially to their very centre on extremely short timescales. This process results in the formation of very compact bulges with specific angular momentum a factor 5-30$smaller than the average angular momentum of the baryons in the whole halo. Such low angular momentum originates both from segregation and effective cancellation when the gas flows to the centre of the halo along well defined cold filamentary streams. These filaments penetrate deep inside the halo and connect to the bulge from multiple rapidly changing directions. Structures falling in along the filaments (satellite galaxies) or formed by gravitational instabilities triggered by the inflow (star clusters) further reduce the angular momentum of the gas in the bulge. Finally, the fraction of gas radially falling to the centre appears to increase with the mass of the halo; we argue that this is most likely due to an enhanced cancellation of angular momentum in rarer halos which are fed by more isotropically distributed cold streams. Such an increasingly efficient funnelling of low-angular momentum gas to the centre of very massive halos at high redshift may account for the rapid pace at which the most massive supermassive black holes grow to reach observed masses around $10^9$M$_\odot$ at an epoch when the Universe is barely 1 Gyr old.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Lyman-alpha emission properties of simulated galaxies: interstellar medium structure and inclination effects

[abridged] Aims. The aim of this paper is to assess the impact of the interstellar medium (ISM) physics on Lyman-alpha (Lya) radiation transfer and to quantify how galaxy orientation with respect to the line of sight alters observational signatures. Methods. We compare the results of Lya radiation transfer calculations through the ISM of a couple of idealized galaxy simulations with different ISM models. Results. First, the small-scale structuration of the ISM plays a determinant role in shaping a galaxys Lya properties.The artificially warm, and hence smooth, ISM of G1 yields an escape fraction of 50 percent at the Lya line center, and produces symmetrical double-peak profiles. On the contrary, in G2, most young stars are embedded in thick star-forming clouds, and the result is a 10 times lower escape fraction. G2 also displays a stronger outflowing velocity field, which favors the escape of red-shifted photons, resulting in an asymmetric Lya line. Second, the Lya properties of G2 strongly depend on the inclination at which it is observed: From edge-on to face-on, the line goes from a double-peak profile with an equivalent width of -5 Angstrom to a 15 times more luminous red-shifted asymmetric line with EW 90 Angstrom. Conclusions. Lya radiation transfer calculations can only lead to realistic properties in simulations where galaxies are resolved into giant molecular clouds, putting these calculations out of reach of current large scale cosmological simulations. Finally, we find inclination effects to be much stronger for Lya photons than for continuum radiation. This could potentially introduce severe biases in the selection function of narrow-band Lya emitter surveys, which could indeed miss a significant fraction of the high-z galaxy population.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Satellite Survival in Highly Resolved Milky Way Class Halos

Surprisingly little is known about the origin and evolution of the Milky Way's satellite galaxy companions. UV photoionisation, supernova feedback and interactions with the larger host halo are all thought to play a role in shaping the population of satellites that we observe today, but there is still no consensus as to which of these effects, if any, dominates. In this paper, we revisit the issue by re-simulating a Milky Way class dark matter (DM) halo with unprecedented resolution. Our set of cosmological hydrodynamic Adaptive Mesh Refinement (AMR) simulations, called the Nut suite, allows us to investigate the effect of supernova feedback and UV photoionisation at high redshift with sub-parsec resolution. We subsequently follow the effect of interactions with the Milky Way-like halo using a lower spatial resolution (50pc) version of the simulation down to z=0. This latter produces a population of simulated satellites that we compare to the observed satellites of the Milky Way and M31. We find that supernova feedback reduces star formation in the least massive satellites but enhances it in the more massive ones. Photoionisation appears to play a very minor role in suppressing star and galaxy formation in all progenitors of satellite halos. By far the largest effect on the satellite population is found to be the mass of the host and whether gas cooling is included in the simulation or not. Indeed, inclusion of gas cooling dramatically reduces the number of satellites captured at high redshift which survive down to z=0.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

AGN feedback using AMR cosmological simulations

Feedback processes are thought to solve some of the long-standing issues of the numerical modelling of galaxy formation: over-cooling, low angular momentum, massive blue galaxies, extra-galactic enrichment, etc. The accretion of gas onto super-massive black holes in the centre of massive galaxies can release tremendous amounts of energy to the surrounding medium. We show, with cosmological Adaptive Mesh Refinement simulations, how the growth of black holes is regulated by the feedback from Active Galactic Nuclei using a new dual jet/heating mechanism. We discuss how this large amount of feedback is able to modify the cold baryon content of galaxies, and perturb the properties of the hot plasma in their vicinity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Are cold flows detectable with metal absorption lines?

[Abridged] Cold gas flowing within the "cosmic web" is believed to be an important source of fuel for star formation at high redshift. However, the presence of such filamentary gas has never been observationally confirmed. In this work, we investigate in detail whether such cold gas is detectable using low-ionisation metal absorption lines, such as CII \lambda1334 as this technique has a proven observational record for detecting gaseous structures. Using a large statistical sample of galaxies from the Mare Nostrum N-body+AMR cosmological simulation, we find that the typical covering fraction of the dense, cold gas in 10^12 Msun haloes at z~2.5 is lower than expected (~5%). In addition, the absorption signal by the interstellar medium of the galaxy itself turns out to be so deep and so broad in velocity space that it completely drowns that of the filamentary gas. A detectable signal might be obtained from a cold filament exactly aligned with the line of sight, but this configuration is so unlikely that it would require surveying an overwhelmingly large number of candidate galaxies to tease it out. Finally, the predicted metallicity of the cold gas in filaments is extremely low (\leq 0.001 Zsun). Should this result persist when higher resolution runs are performed, it would significantly increase the difficulty of detecting filamentary gas inflows using metal lines. However, even if we assume that filaments are enriched to Zsun, the absorption signal that we compute is still weak. We are therefore led to conclude that it is extremely difficult to observationally prove or disprove the presence of cold filaments as the favorite accretion mode of galaxies using low-ionisation metal absorption lines. The Ly-alpha emission route looks more promising but due to the resonant nature of the line, radiative transfer simulations are required to fully characterize the observed signal.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

How AGN feedback and metal cooling shape cluster entropy profiles

Observed clusters of galaxies essentially come in two flavors: non cool core clusters characterized by an isothermal temperature profile and a central entropy floor, and cool-core clusters where temperature and entropy in the central region are increasing with radius. Using cosmological resimulations of a galaxy cluster, we study the evolution of its intracluster medium (ICM) gas properties, and through them we assess the effect of different (sub-grid) modelling of the physical processes at play, namely gas cooling, star formation, feedback from supernovae and active galactic nuclei (AGN). More specifically we show that AGN feedback plays a major role in the pre-heating of the proto-cluster as it prevents a high concentration of mass from collecting in the center of the future galaxy cluster at early times. However, AGN activity during the cluster's later evolution is also required to regulate the mass flow into its core and prevent runaway star formation in the central galaxy. Whereas the energy deposited by supernovae alone is insufficient to prevent an overcooling catastrophe, supernovae are responsible for spreading a large amount of metals at high redshift, enhancing the cooling efficiency of the ICM gas. As the AGN energy release depends on the accretion rate of gas onto its central black hole engine, the AGN responds to this supernova enhanced gas accretion by injecting more energy into the surrounding gas, and as a result increases the amount of early pre-heating. We demonstrate that the interaction between an AGN jet and the ICM gas that regulates the growth of the AGN's BH, can naturally produce cool core clusters if we neglect metals. However, as soon as metals are allowed to contribute to the radiative cooling, only the non cool core solution is produced.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

How Does Feedback Affect Milky Way Satellite Formation?

We use sub-parsec resolution hydrodynamic resimulations of a Milky Way (MW) like galaxy at high redshift to investigate the formation of the MW satellite galaxies. More specifically, we assess the impact of supernova feedback on the dwarf progenitors of these satellite, and the efficiency of a simple instantaneous reionisation scenario in suppressing star formation at the low-mass end of this dwarf distribution. Identifying galaxies in our high redshift simulation and tracking them to z=0 using a dark matter halo merger tree, we compare our results to present-day observations and determine the epoch at which we deem satellite galaxy formation must be completed. We find that only the low-mass end of the population of luminous subhalos of the Milky-Way like galaxy is not complete before redshift 8, and that although supernovae feedback reduces the stellar mass of the low-mass subhalos (log(M/Msolar) < 9), the number of surviving satellites around the Milky-Way like galaxy at z = 0 is the same in the run with or without supernova feedback. If a luminous halo is able to avoid accretion by the Milky-Way progenitor before redshift 3, then it is likely to survive as a MW satellite to redshift 0.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Self-regulated growth of supermassive black holes by a dual jet/heating AGN feedback mechanism: methods, tests and implications for cosmological simulations

We develop a new sub-grid model for the growth of supermassive Black Holes (BHs) and their associated Active Galactic Nuclei (AGN) feedback in hydrodynamical cosmological simulations. Assuming that BHs are created in the early stages of galaxy formation, they grow by mergers and accretion of gas at a Eddington-limited Bondi accretion rate. However this growth is regulated by AGN feedback which we model using two different modes: a quasar-heating mode when accretion rates onto the BHs are comparable to the Eddington rate, and a radio-jet mode at lower accretion rates. In other words, our feedback model deposits energy as a succession of thermal bursts and jet outflows depending on the properties of the gas surrounding the BHs. We assess the plausibility of such a model by comparing our results to observational measurements of the coevolution of BHs and their host galaxy properties, and check their robustness with respect to numerical resolution. We show that AGN feedback must be a crucial physical ingredient for the formation of massive galaxies as it appears to be the only physical mechanism able to efficiently prevent the accumulation of and/or expel cold gas out of halos/galaxies and significantly suppress star formation. Our model predicts that the relationship between BHs and their host galaxy mass evolves as a function of redshift, because of the vigorous accretion of cold material in the early Universe that drives Eddington-limited accretion onto BHs. Quasar activity is also enhanced at high redshift. However, as structures grow in mass and lose their cold material through star formation and efficient BH feedback ejection, the AGN activity in the low-redshift Universe becomes more and more dominated by the radio mode, which powers jets through the hot circum-galactic medium.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The angular momentum of baryons and dark matter halos revisited

Recent theoretical studies have shown that galaxies at high redshift are fed by cold, dense gas filaments, suggesting angular momentum transport by gas differs from that by dark matter. Revisiting this issue using high-resolution cosmological hydrodynamics simulations with adaptive mesh refinement, we find that at the time of accretion, gas and dark matter do carry a similar amount of specific angular momentum, but that it is systematically higher than that of the dark matter halo as a whole. At high redshift, freshly accreted gas rapidly streams into the central region of the halo, directly depositing this large amount of angular momentum within a sphere of radius r=0.1rvir. In contrast, dark matter particles pass through the central region unscathed, and a fraction of them ends up populating the outer regions of the halo (r/rvir>0.1), redistributing angular momentum in the process. As a result, large-scale motions of the cosmic web have to be considered as the origin of gas angular momentum rather than its virialised dark matter halo host. This generic result holds for halos of all masses at all redshifts, as radiative cooling ensures that a significant fraction of baryons remain trapped at the centre of the halos. Despite this injection of angular momentum enriched gas, we predict an amount for stellar discs which is in fair agreement with observations at z=0. This arises because the total specific angular momentum of the baryons remains close to that of dark matter halos. We propose a new scenario where gas efficiently carries the angular momentum generated by large-scale structure motions deep inside dark matter halos, redistributing it only in the vicinity of the disc.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The impact of ISM turbulence, clustered star formation and feedback on galaxy mass assembly through cold flows and mergers

Two of the dominant channels for galaxy mass assembly are cold flows (cold gas supplied via the filaments of the cosmic web) and mergers. How these processes combine in a cosmological setting, at both low and high redshift, to produce the whole zoo of galaxies we observe is largely unknown. Indeed there is still much to understand about the detailed physics of each process in isolation. While these formation channels have been studied using hydrodynamical simulations, here we study their impact on gas properties and star formation (SF) with some of the first simulations that capture the multiphase, cloudy nature of the interstellar medium (ISM), by virtue of their high spatial resolution (and corresponding low temperature threshold). In this regime, we examine the competition between cold flows and a supernovae(SNe)-driven outflow in a very high-redshift galaxy (z {\approx} 9) and study the evolution of equal-mass galaxy mergers at low and high redshift, focusing on the induced SF. We find that SNe-driven outflows cannot reduce the cold accretion at z {\approx} 9 and that SF is actually enhanced due to the ensuing metal enrichment. We demonstrate how several recent observational results on galaxy populations (e.g. enhanced HCN/CO ratios in ULIRGs, a separate Kennicutt Schmidt (KS) sequence for starbursts and the population of compact early type galaxies (ETGs) at high redshift) can be explained with mechanisms captured in galaxy merger simulations, provided that the multiphase nature of the ISM is resolved.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The impact of supernovae driven winds on stream-fed protogalaxies

SNe driven winds are widely thought to be very influential in the high-redshift Universe, shaping the properties of the circum-galactic medium, enriching the IGM with metals and driving the evolution of low-mass galaxies. However, it is not yet fully understood how SNe driven winds interact with their surroundings in a cosmological context, nor is it clear whether they are able to significantly impact the evolution of low-mass galaxies from which they originate by altering the amount of cold material these accrete from the cosmic web. We implement a standard Taylor-Sedov type solution, widely used in the community to depict the combined action of many SN explosions, in a cosmological resimulation of a low mass galaxy at z =9 from the 'Nut' suite. However, in contrast with previous work, we achieve a resolution high enough to capture individual SN remnants in the Taylor-Sedov phase, for which the solution provides an accurate description of the expansion. We report the development of a high-velocity, far-reaching galactic wind produced by the combined action of SNe in the main galaxy and its satellites, which are located in the same or a neighbouring dark matter halo. Despite this, we find that (i) this wind carries out very little mass (the measured outflow is of the order of a tenth of the inflow/star formation rate) and (ii) the cold gas inflow rate remains essentially unchanged from the run without SNe feedback. Moreover, there are epochs during which star formation is enhanced in the feedback run relative to its radiative cooling only counterpart. We attribute this 'positive' feedback to the metal enrichment that is present only in the former. We conclude that at very high redshift, efficient SNe feedback can drive large-scale galactic winds but does not prevent massive cold gas inflow from fuelling galaxies, resulting in long-lived episodes of intense star formation.(abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Galactic star formation in parsec-scale resolution simulations

The interstellar medium (ISM) in galaxies is multiphase and cloudy, with stars forming in the very dense, cold gas found in Giant Molecular Clouds (GMCs). Simulating the evolution of an entire galaxy, however, is a computational problem which covers many orders of magnitude, so many simulations cannot reach densities high enough or temperatures low enough to resolve this multiphase nature. Therefore, the formation of GMCs is not captured and the resulting gas distribution is smooth, contrary to observations. We investigate how star formation (SF) proceeds in simulated galaxies when we obtain parsec-scale resolution and more successfully capture the multiphase ISM. Both major mergers and the accretion of cold gas via filaments are dominant contributors to a galaxy's total stellar budget and we examine SF at high resolution in both of these contexts.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Jet-regulated cooling catastrophe

We present the first implementation of Active Galactic Nuclei (AGN) feedback in the form of momentum driven jets in an Adaptive Mesh Refinement (AMR) cosmological resimulation of a galaxy cluster. The jets are powered by gas accretion onto Super Massive Black Holes (SMBHs) which also grow by mergers. Throughout its formation, the cluster experiences different dynamical states: both a morphologically perturbed epoch at early times and a relaxed state at late times allowing us to study the different modes of BH growth and associated AGN jet feedback. BHs accrete gas efficiently at high redshift (z>2), significantly pre-heating proto-cluster halos. Gas-rich mergers at high redshift also fuel strong, episodic jet activity, which transports gas from the proto-cluster core to its outer regions. At later times, while the cluster relaxes, the supply of cold gas onto the BHs is reduced leading to lower jet activity. Although the cluster is still heated by this activity as sound waves propagate from the core to the virial radius, the jets inefficiently redistribute gas outwards and a small cooling flow develops, along with low-pressure cavities similar to those detected in X-ray observations. Overall, our jet implementation of AGN feedback quenches star formation quite efficiently, reducing the stellar content of the central cluster galaxy by a factor 3 compared to the no AGN case. It also dramatically alters the shape of the gas density profile, bringing it in close agreement with the beta model favoured by observations, producing quite an isothermal galaxy cluster for gigayears in the process. However, it still falls short in matching the lower than Universal baryon fractions which seem to be commonplace in observed galaxy clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Skeleton: Connecting Large Scale Structures to Galaxy Formation

We report on two quantitative, morphological estimators of the filamentary structure of the Cosmic Web, the so-called global and local skeletons. The first, based on a global study of the matter density gradient flow, allows us to study the connectivity between a density peak and its surroundings, with direct relevance to the anisotropic accretion via cold flows on galactic halos. From the second, based on a local constraint equation involving the derivatives of the field, we can derive predictions for powerful statistics, such as the differential length and the relative saddle to extrema counts of the Cosmic web as a function of density threshold (with application to percolation of structures and connectivity), as well as a theoretical framework to study their cosmic evolution through the onset of gravity-induced non-linearities.

0 citations0 reviewsNo code linkedOpen access