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

Robert Feldmann

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

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

Catalog footprint

What is connected

32works
6topics
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

32 published item(s)

preprint2022arXiv

Black hole -- galaxy scaling relations in FIRE: the importance of black hole location and mergers

The concurrent growth of supermassive black holes (SMBHs) and their host galaxies remains to be fully explored, especially at high redshift. While often understood as a consequence of self-regulation via AGN feedback, it can also be explained by alternative SMBH accretion models. Here, we expand on previous work by studying the growth of SMBHs with the help of a large suite of cosmological zoom-in simulations (MassiveFIRE) that are part of the Feedback in Realistic Environments (FIRE) project. The growth of SMBHs is modelled in post-processing with different black hole accretion models, placements, and merger treatments, and validated by comparing to on-the-fly calculations. Scaling relations predicted by the gravitational torque driven accretion (GTDA) model agree with observations at low redshift without the need for AGN feedback, in contrast to models in which the accretion rate depends strongly on SMBH mass. At high redshift, we find deviations from the local scaling relations in line with previous theoretical results. In particular, SMBHs are under-massive, presumably due to stellar feedback, but start to grow efficiently once their host galaxies reach $M_* \sim 10^{10} M_{\odot}$. We analyse and explain these findings in the context of a simple analytic model. Finally, we show that the predicted scaling relations depend sensitively on the SMBH location and the efficiency of SMBH merging, particularly in low-mass systems. These findings highlight the relevance of understanding the evolution of SMBH-galaxy scaling relations to predict the rate of gravitational wave signals from SMBH mergers across cosmic history.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

From EMBER to FIRE: predicting high resolution baryon fields from dark matter simulations with Deep Learning

Hydrodynamic simulations provide a powerful, but computationally expensive, approach to study the interplay of dark matter and baryons in cosmological structure formation. Here we introduce the EMulating Baryonic EnRichment (EMBER) Deep Learning framework to predict baryon fields based on dark-matter-only simulations thereby reducing computational cost. EMBER comprises two network architectures, U-Net and Wasserstein Generative Adversarial Networks (WGANs), to predict two-dimensional gas and HI densities from dark matter fields. We design the conditional WGANs as stochastic emulators, such that multiple target fields can be sampled from the same dark matter input. For training we combine cosmological volume and zoom-in hydrodynamical simulations from the Feedback in Realistic Environments (FIRE) project to represent a large range of scales. Our fiducial WGAN model reproduces the gas and HI power spectra within 10% accuracy down to ~10 kpc scales. Furthermore, we investigate the capability of EMBER to predict high resolution baryon fields from low resolution dark matter inputs through upsampling techniques. As a practical application, we use this methodology to emulate high-resolution HI maps for a dark matter simulation of a L=100 Mpc/h comoving cosmological box. The gas content of dark matter haloes and the HI column density distributions predicted by EMBER agree well with results of large volume cosmological simulations and abundance matching models. Our method provides a computationally efficient, stochastic emulator for augmenting dark matter only simulations with physically consistent maps of baryon fields.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Galaxies lacking dark matter produced by close encounters in a cosmological simulation

The standard cold dark matter plus cosmological constant model predicts that galaxies form within dark-matter haloes, and that low-mass galaxies are more dark-matter dominated than massive ones. The unexpected discovery of two low-mass galaxies lacking dark matter immediately provoked concerns about the standard cosmology and ignited explorations of alternatives, including self-interacting dark matter and modified gravity. Apprehension grew after several cosmological simulations using the conventional model failed to form adequate numerical analogues with comparable internal characteristics (stellar masses, sizes, velocity dispersions and morphologies). Here we show that the standard paradigm naturally produces galaxies lacking dark matter with internal characteristics in agreement with observations. Using a state-of-the-art cosmological simulation and a meticulous galaxy-identification technique, we find that extreme close encounters with massive neighbours can be responsible for this. We predict that approximately 30 percent of massive central galaxies (with at least 1e11 solar masses in stars) harbour at least one dark-matter-deficient satellite (with 1e8 - 1e9 solar masses in stars). This distinctive class of galaxies provides an additional layer in our understanding of the role of interactions in shaping galactic properties. Future observations surveying galaxies in the aforementioned regime will provide a crucial test of this scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Hot-mode accretion and the physics of thin-disk galaxy formation

We use FIRE simulations to study disk formation in z~0, Milky Way-mass galaxies, and conclude that a key ingredient for the formation of thin stellar disks is the ability for accreting gas to develop an aligned angular momentum distribution via internal cancellation *prior* to joining the galaxy. Among galaxies with a high fraction (>70%) of their young stars in a thin disk (h/R~0.1) we find that: (i) hot, virial-temperature gas dominates the inflowing gas mass on halo scales (>~20 kpc), with radiative losses offset by compression heating; (ii) this hot accretion proceeds until angular momentum support slows inward motion, at which point the gas cools to T~10^4 K or less; (iii) prior to cooling, the accreting gas develops an angular momentum distribution that is aligned with the galaxy disk, and while cooling transitions from a quasi-spherical spatial configuration to a more flattened, disk-like configuration. We show that the existence of this "rotating cooling flow" accretion mode is strongly correlated with the fraction of stars forming in a thin disk among a sample of 17 z~0 galaxies spanning a halo mass range of 10^10.5 solar masses to 10^12 solar masses, or a stellar mass range 10^8 solar masses to 10^11 solar masses. Notably, galaxies with a thick disk or irregular morphology do not undergo significant angular momentum alignment of gas prior to accretion and show no correspondence between halo gas cooling and flattening. Our results suggest that rotating cooling flows (or, more generally, rotating subsonic flows) that become coherent and angular momentum-supported prior to accretion onto the galaxy are likely a necessary condition for the formation of thin, star-forming disk galaxies in a LambdaCDM universe.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quenching and the UVJ diagram in the SIMBA cosmological simulation

Over the past decade, rest-frame color-color diagrams have become popular tools for selecting quiescent galaxies at high redshift, breaking the color degeneracy between quiescent and dust-reddened star-forming galaxies. In this work, we study one such color-color selection tool -- the rest-frame $U-V$ vs. $V-J$ diagram -- by employing mock observations of cosmological galaxy formation simulations. In particular, we conduct numerical experiments assessing both trends in galaxy properties in UVJ space and the color-color evolution of massive galaxies as they quench at redshifts $z\sim 1$--$2$. We find that our models broadly reproduce the observed UVJ diagram at $z=1$--$2$, including (for the first time in a cosmological simulation) reproducing the population of extremely dust-reddened galaxies in the top right of the UVJ diagram. However, our models primarily populate this region with low-mass galaxies and do not produce as clear a bimodality between star-forming and quiescent galaxies as is seen in observations. The former issue is due to an excess of dust in low-mass galaxies and relatively gray attenuation curves in high-mass galaxies, while the latter is due to the overpopulation of the green valley in SIMBA. When investigating the time evolution of galaxies on the UVJ diagram, we find that the quenching pathway on the UVJ diagram is independent of the quenching timescale, and instead dependent primarily on the average specific star formation rate in the 1 Gyr prior to the onset of quenching. Our results support the interpretation of different quenching pathways as corresponding to the divergent evolution of post-starburst and green valley galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Star Formation Suppresion by Tidal Removal of Cold Molecular Gas from an Intermediate-Redshift Massive Post-Starburst Galaxy

Observations and simulations have demonstrated that star formation in galaxies must be actively suppressed to prevent the formation of over-massive galaxies. Galactic outflows driven by stellar feedback or supermassive black hole accretion are often invoked to regulate the amount of cold molecular gas available for future star formation, but may not be the only relevant quenching processes in all galaxies. We present the discovery of vast molecular tidal features extending up to 64 kpc outside of a massive z=0.646 post-starburst galaxy that recently concluded its primary star-forming episode. The tidal tails contain (1.2 +/- 0.1)x10^10 Msun of molecular gas, 47 +/- 5 % of the total cold gas reservoir of the system. Both the scale and magnitude of the molecular tidal features are unprecedented compared to all known nearby or high-redshift merging systems. We infer that the cold gas was stripped from the host galaxies during the merger, which is most likely responsible for triggering the initial burst phase and the subsequent suppression of star formation. While only a single example, this result shows that galaxy mergers can regulate the cold gas contents in distant galaxies by directly removing a large fraction of the molecular gas fuel, and plausibly suppress star formation directly, a qualitatively different physical mechanism than feedback-driven outflows.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Compact Structures of Massive $z\sim0.7$ Post-Starburst Galaxies in the SQuIGG$\vec{L}$E Sample

We present structural measurements of 145 spectroscopically selected intermediate-redshift (z$\sim$0.7), massive ($M_\star \sim 10^{11} \ M_\odot$) post-starburst galaxies from the SQuIGG$\vec{L}$E Sample measured using wide-depth Hyper Suprime-Cam i-band imaging. This deep imaging allows us to probe the sizes and structures of these galaxies, which we compare to a control sample of star forming and quiescent galaxies drawn from the LEGA-C Survey. We find that post-starburst galaxies systematically lie $\sim0.1$ dex below the quiescent mass-size (half-light radius) relation, with a scatter of $\sim0.2$ dex. This finding is bolstered by non-parametric measures, such as the Gini coefficient and the concentration, which also reveal these galaxies to have more compact light profiles than both quiescent and star-forming populations at similar mass and redshift. The sizes of post-starburst galaxies show either negative or no correlation with the time since quenching, such that more recently quenched galaxies are larger or similarly sized. This empirical finding disfavors the formation of post-starburst galaxies via a purely central burst of star formation that simultaneously shrinks the galaxy and shuts off star formation. We show that the central densities of post-starburst and quiescent galaxies at this epoch are very similar, in contrast with their effective radii. The structural properties of z$\sim$0.7 post-starburst galaxies match those of quiescent galaxies that formed in the early universe, suggesting that rapid quenching in the present epoch is driven by a similar mechanism to the one at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The observability of galaxy merger signatures in nearby gas-rich spirals

Galaxy mergers are crucial to understanding galaxy evolution, therefore we must determine their observational signatures to select them from large IFU galaxy samples such as MUSE and SAMI. We employ 24 high-resolution idealised hydrodynamical galaxy merger simulations based on the "Feedback In Realistic Environment" (FIRE-2) model to determine the observability of mergers to various configurations and stages using synthetic images and velocity maps. Our mergers cover a range of orbital configurations at fixed 1:2.5 stellar mass ratio for two gas rich spirals at low redshift. Morphological and kinematic asymmetries are computed for synthetic images and velocity maps spanning each interaction. We divide the interaction sequence into three: (1) the pair phase; (2) the merging phase; and (3) the post-coalescence phase. We correctly identify mergers between first pericentre passage and 500 Myr after coalescence using kinematic asymmetry with 66% completeness, depending upon merger phase and the field-of-view of the observation. We detect fewer mergers in the pair phase (40%) and many more in the merging and post-coalescence phases (97%). We find that merger detectability decreases with field-of-view, except in retrograde mergers, where centrally concentrated asymmetric kinematic features enhances their detectability. Using a cut-off derived from a combination of photometric and kinematic asymmetry, we increase these detections to 89% overall, 79% in pairs, and close to 100% in the merging and post-coalescent phases. By using this combined asymmetry cut-off we mitigate some of the effects caused by smaller fields-of-view subtended by massively multiplexed integral field spectroscopy programmes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The rotation of planet-hosting stars

Understanding the distribution of angular momentum during the formation of planetary systems is a key topic in astrophysics. Data from the $\textit{Kepler}$ and $\textit{Gaia}$ missions allow to investigate whether stellar rotation is correlated with the presence of planets around Sun-like stars. Here, we perform a statistical analysis of the rotation period of 493 planet-hosting stars. These are matched to a control sample, without detected planets, with similar effective temperatures, masses, radii, metallicities, and ages. We find that planet-hosting stars rotate on average $1.63 \pm 0.40$ days slower. The difference in rotation is statistically significant both in samples including and not including planets confirmed by radial velocity follow-up observations. We also analyse the dependence of rotation distribution on various stellar and planetary properties. Our results could potentially be explained by planet detection biases depending on the rotation period of their host stars in both RV and transit methods. Alternatively, they could point to a physical link between the existence of planets and stellar rotation, emphasising the need to understand the role of angular momentum in the formation and evolution planetary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Now you see it, now you don't: Star formation truncation precedes the loss of molecular gas by ~100 Myr in massive post-starburst galaxies at z~0.6

We use ALMA observations of CO(2-1) in 13 massive ($M_{\star}\gtrsim 10^{11} M_{\odot}$) post-starburst galaxies at $z\sim0.6$ to constrain the molecular gas content in galaxies shortly after they quench their major star-forming episode. The post-starburst galaxies in this study are selected from the Sloan Digital Sky Survey spectroscopic samples (DR14) based on their spectral shapes, as part of the SQuIGGLE program. Early results showed that two post-starburst galaxies host large H$_2$ reservoirs despite their low inferred star formation rates. Here we expand this analysis to a larger statistical sample of 13 galaxies. Six of the primary targets (45%) are detected, with $M_{H_2}\gtrsim10^9 M_{\odot}$. Given their high stellar masses, this mass limit corresponds to an average gas fraction of $\langle f_{H_2} \equiv M_{H_2}/M_{\star} \rangle \sim7\%$, or ${\sim}14\%$ using lower stellar masses estimates derived from analytic, exponentially declining star formation histories. The gas fraction correlates with the $D_n4000$ spectral index, suggesting that the cold gas reservoirs decrease with time since burst, as found in local K+A galaxies. Star formation histories derived from flexible stellar population synthesis modeling support this empirical finding: galaxies that quenched $\lesssim 150$ Myr prior to observation host detectable CO(2-1) emission, while older post-starburst galaxies are undetected. The large $\mathrm{H_2}$ reservoirs and low star formation rates in the sample imply that the quenching of star formation precedes the disappearance of the cold gas reservoirs. However, within the following 100-200 Myrs, the SQuIGGLE galaxies require the additional and efficient heating or removal of cold gas to bring their low star formation rates in line with standard $\mathrm{H_2}$ scaling relations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Out of sight, out of mind? The impact of correlated clustering in substructure lensing

A promising route for revealing the existence of dark matter structures on mass scales smaller than the faintest galaxies is through their effect on strong gravitational lenses. We examine the role of local, lens-proximate clustering in boosting the lensing probability relative to contributions from substructure and unclustered line-of-sight (LOS) halos. Using two cosmological simulations that can resolve halo masses of $M_{\rm halo} \simeq 10^{9}\ M_{\odot}$ (in a simulation box of length $L_{\rm box}{\sim}100\,{\rm Mpc}$) and $10^{7}\ M_{\odot}$ ($L_{\rm box}\sim20\,{\rm Mpc}$), we demonstrate that clustering in the vicinity of the lens host produces a clear enhancement relative to an assumption of unclustered halos that persists to $> 20\,R_{\rm vir}$. This enhancement exceeds estimates that use a two-halo term to account for clustering, particularly within $2-5\,R_{\rm vir}$. We provide an analytic expression for this excess, clustered contribution. We find that local clustering boosts the expected count of $10^9 \ M_\odot$ perturbing halos by ${\sim}35\%$ compared to substructure alone, a result that will significantly enhance expected signals for low-redshift ($z_l \simeq 0.2$) lenses, where substructure contributes substantially compared to LOS halos. We also find that the orientation of the lens with respect to the line of sight (e.g., whether the line of sight passes through the major axis of the lens) can also have a significant effect on the lensing signal, boosting counts by an additional $\sim 50\%$ compared to a random orientations. This could be important if discovered lenses are biased to be oriented along their principal axis.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

SQuIGGLE: Studying Quenching in Intermediate-z Galaxies -- Gas, AnguLar Momentum, and Evolution

We describe the SQuIGGLE survey of intermediate-redshift post-starburst galaxies. We leverage the large sky coverage of the SDSS to select ~1300 recently-quenched galaxies at 0.5<z<~0.9 based on their unique spectral shapes. These bright, intermediate-redshift galaxies are ideal laboratories to study the physics responsible for the rapid quenching of star formation: they are distant enough to be useful analogs for high-redshift quenching galaxies, but low enough redshift that multi-wavelength follow-up observations are feasible with modest telescope investments. We use the Prospector code to infer the stellar population properties and non-parametric star formation histories of all galaxies in the sample. We find that SQuIGGLE galaxies are both very massive (M* ~ 10^11.25 Msun) and quenched, with inferred star formation rates <~1Msun/yr, more than an order of magnitude below the star-forming main sequence. The best-fit star formation histories confirm that these galaxies recently quenched a major burst of star formation: >75% of SQuIGGLE galaxies formed at least a quarter of their total stellar mass in the recent burst, which ended just ~200Myr before observation. We find that SQuIGGLE galaxies are on average younger and more burst-dominated than most other z<~1 post-starburst samples. This large sample of bright post-starburst galaxies at intermediate redshift opens a wide range of studies into the quenching process. In particular, the full SQuIGGLE survey will investigate the molecular gas reservoirs, morphologies, kinematics, resolved stellar populations, AGN incidence, and infrared properties of this unique sample of galaxies in order to place definitive constraints on the quenching process.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Measuring dynamical masses from gas kinematics in simulated high-redshift galaxies

Advances in instrumentation have recently extended detailed measurements of gas kinematics to large samples of high-redshift galaxies. Relative to most nearby, thin disk galaxies, in which gas rotation accurately traces the gravitational potential, the interstellar medium (ISM) of z>1 galaxies is typically more dynamic and exhibits elevated turbulence. If not properly modeled, these effects can strongly bias dynamical mass measurements. We use high-resolution FIRE-2 cosmological zoom-in simulations to analyze the physical effects that must be considered to correctly infer dynamical masses from gas kinematics. Our analysis covers a range of galaxy properties from low-redshift Milky-Way-mass galaxies to massive high-redshift galaxies (M_* > 10^11 M_sun at z=1). Selecting only snapshots where a disk is present, we calculate the rotational profile v_phi(r) of the cool (10^3.5 K < T < 10^4.5 K) gas and compare it to the circular velocity v_c=sqrt(GM/r). In the simulated galaxies, the gas rotation traces the circular velocity at intermediate radii, but the two quantities diverge significantly in the center and in the outer disk. Our simulations appear to over-predict observed rotational velocities in the centers of massive galaxies (likely from a lack of black hole feedback), so we focus on larger radii. Gradients in the turbulent pressure at these radii can provide additional radial support and bias dynamical mass measurements low by up to 40%. In both the interior and exterior, the gas' motion can be significantly non-circular due to e.g. bars, satellites, and inflows/outflows. We discuss the accuracy of commonly-used analytic models for pressure gradients (or "asymmetric drift") in the ISM of high-redshift galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Positive feedback at the disc-halo interface

The flat star formation (SF) history of the Milky Way (MW) requires gas in the Galactic disc to be replenished, most likely from a reservoir outside the Galaxy. Such a replenishment may be achieved by a form of `positive' feedback, whereby SF feedback creates a Galactic fountain cycle that collects and cools additional gas from the hot halo surrounding the Galaxy. In this paper we present a model of this process for the MW. A section of the Galactic disc is allowed to form stars which subsequently explode as supernovae (SNe) and send gas out into the hot halo. The gas that is sent out is colder than the hot halo gas and, as it mixes, the halo gas is cooled, providing fuel for further SF as the mixture falls back onto the Galactic disc. We find that this process can be sufficient to maintain a roughly-constant cold gas mass in the MW over at least 3 Gyr. Our results further suggest that there is a positive feedback trend whereby increasing SF leads to an increase in the cold gas budget at average SF rates below 0.5 solar masses per year but above which becomes negative, where further increasing the SFR causes the cold gas budget to decrease. We have constructed an analytical model for this that reproduces the data well and could have profound implications for galaxy evolution in feedback-dominated regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The AGORA high-resolution galaxy simulations comparison project: Public data release

As part of the AGORA High-resolution Galaxy Simulations Comparison Project (Kim et al. 2014, 2016) we have generated a suite of isolated Milky Way-mass galaxy simulations using 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical galaxy formation community. In these simulations we adopted identical galactic disk initial conditions, and common physics models (e.g., radiative cooling and ultraviolet background by a standardized package). Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production were carefully constrained. Here we release the simulation data to be freely used by the community. In this release we include the disk snapshots at 0 and 500Myr of evolution per each code as used in Kim et al. (2016), from simulations with and without star formation and feedback. We encourage any member of the numerical galaxy formation community to make use of these resources for their research - for example, compare their own simulations with the AGORA galaxies, with the common analysis yt scripts used to obtain the plots shown in our papers, also available in this release.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The IRX-$β$ relation of high-redshift galaxies

The relation between infrared excess (IRX) and UV spectral slope ($β_{\rm UV}$) is an empirical probe of dust properties of galaxies. The shape, scatter, and redshift evolution of this relation are not well understood, however, leading to uncertainties in estimating the dust content and star formation rates (SFRs) of galaxies at high redshift. In this study, we explore the nature and properties of the IRX-$β_{\rm UV}$ relation with a sample of $z=2-6$ galaxies ($M_*\approx 10^9-10^{12}\,M_\odot$) extracted from high-resolution cosmological simulations (MassiveFIRE) of the Feedback in Realistic Environments (FIRE) project. The galaxies in our sample show an IRX-$β_{\rm UV}$ relation that is in good agreement with the observed relation in nearby galaxies. IRX is tightly coupled to the UV optical depth, and is mainly determined by the dust-to-star geometry instead of total dust mass, while $β_{\rm UV}$ is set both by stellar properties, UV optical depth, and the dust extinction law. Overall, much of the scatter in the IRX-$β_{\rm UV}$ relation of our sample is found to be driven by variations of the intrinsic UV spectral slope. We further assess how the IRX-$β_{\rm UV}$ relation depends on viewing direction, dust-to-metal ratio, birth-cloud structures, and the dust extinction law and we present a simple model that encapsulates most of the found dependencies. Consequently, we argue that the reported `deficit' of the infrared/sub-millimetre bright objects at $z>5$ does not necessarily imply a non-standard dust extinction law at those epochs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Role of Active Galactic Nuclei in the Quenching of Massive Galaxies in the SQuiGGLE Survey

We study the incidence of nuclear activity in a large sample of massive post-starburst galaxies at z~0.7 selected from the Sloan Digital Sky Survey, and identify active galactic nuclei based on radio continuum and optical emission lines. Over our mass range of 10^10.6-10^11.5 Msun, the incidence of radio activity is weakly dependent on stellar mass and independent of stellar age, while radio luminosity depends strongly on stellar mass. Optical nuclear activity incidence depends most strongly on the Dn4000 line index, a proxy for stellar age, with an active fraction that is ~ten times higher in the youngest versus oldest post-starburst galaxies. Since a similar trend is seen between age and molecular gas fractions, we argue that, like in local galaxies, the age trend reflects a peak in available fueling rather than feedback from the central black hole on the surrounding galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Predicting dark matter halo formation in N-body simulations with deep regression networks

Dark matter haloes play a fundamental role in cosmological structure formation. The most common approach to model their assembly mechanisms is through N-body simulations. In this work we present an innovative pathway to predict dark matter halo formation from the initial density field using a Deep Learning algorithm. We implement and train a Deep Convolutional Neural Network (DCNN) to solve the task of retrieving Lagrangian patches from which dark matter halos will condense. The volumetric multi-label classification task is turned into a regression problem by means of the euclidean distance transformation. The network is complemented by an adaptive version of the watershed algorithm to form the entire protohalo identification pipeline. We show that splitting the segmentation problem into two distinct sub-tasks allows for training smaller and faster networks, while the predictive power of the pipeline remains the same. The model is trained on synthetic data derived from a single full N-body simulation and achieves deviations of ~10% when reconstructing the dark matter halo mass function at z=0. This approach represents a promising framework for learning highly non-linear relations in the primordial density field. As a practical application, our method can be used to produce mock dark matter halo catalogues directly from the initial conditions of N-body simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Giant clumps in the FIRE simulations: a case study of a massive high-redshift galaxy

The morphology of massive star-forming galaxies at high redshift is often dominated by giant clumps of mass ~10^8-10^9 Msun and size ~100-1000 pc. Previous studies have proposed that giant clumps might have an important role in the evolution of their host galaxy, particularly in building the central bulge. However, this depends on whether clumps live long enough to migrate from their original location in the disc or whether they get disrupted by their own stellar feedback before reaching the centre of the galaxy. We use cosmological hydrodynamical simulations from the FIRE (Feedback in Realistic Environments) project that implement explicit treatments of stellar feedback and ISM physics to study the properties of these clumps. We follow the evolution of giant clumps in a massive (stellar mass ~10^10.8 Msun at z=1), discy, gas-rich galaxy from redshift z>2 to z=1. Even though the clumpy phase of this galaxy lasts over a gigayear, individual gas clumps are short-lived, with mean lifetime of massive clumps of ~20 Myr. During that time, they turn between 0.1% and 20% of their gas into stars before being disrupted, similar to local GMCs. Clumps with M>10^7 Msun account for ~20% of the total star formation in the galaxy during the clumpy phase, producing ~10^10 Msun of stars. We do not find evidence for net inward migration of clumps within the galaxy. The number of giant clumps and their mass decrease at lower redshifts, following the decrease in the overall gas fraction and star-formation rate.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The formation of massive, quiescent galaxies at cosmic noon

The cosmic noon (z~1.5-3) marked a period of vigorous star formation for most galaxies. However, about a third of the more massive galaxies at those times were quiescent in the sense that their observed stellar populations are inconsistent with rapid star formation. The reduced star formation activity is often attributed to gaseous outflows driven by feedback from supermassive black holes, but the impact of black hole feedback on galaxies in the young Universe is not yet definitively established. We analyze the origin of quiescent galaxies with the help of ultra-high resolution, cosmological simulations that include feedback from stars but do not model the uncertain consequences of black hole feedback. We show that dark matter halos with specific accretion rates below ~0.25-0.4 per Gyr preferentially host galaxies with reduced star formation rates and red broad-band colors. The fraction of such halos in large dark matter only simulations matches the observed fraction of massive quiescent galaxies (~10^10-10^11 Msun). This strongly suggests that halo accretion rate is the key parameter determining which massive galaxies at z~1.5-3 become quiescent. Empirical models that connect galaxy and halo evolution, such as halo occupation distribution or abundance matching models, assume a tight link between galaxy properties and the masses of their parent halos. These models will benefit from adding the specific accretion rate of halos as a second model parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The impact of stellar feedback on hot gas in galaxy haloes: the Sunyaev-Zel'dovich effect and soft X-ray emission

The thermal Sunyaev-Zel'dovich (SZ) effect and soft X-ray emission are routinely observed around massive galaxies and in galaxy groups and clusters. We study these observational diagnostics of galaxy haloes for a suite of cosmological `zoom-in' simulations from the `Feedback In Realistic Environments' project, which spans a large range in halo mass 10^10-10^13 Msun). We explore the effect of stellar feedback on the hot gas observables. The properties of our simulated groups, such as baryon fractions, SZ flux, and X-ray luminosities (L_X), are broadly consistent with existing observations, even though feedback from active galactic nuclei is not included. We make predictions for future observations of lower-mass objects for both SZ and diffuse X-ray measurements, finding that they are not just scaled-down versions of massive galaxies, but more strongly affected by galactic winds driven by star formation. Low-mass haloes (<~10^11 Msun) retain a low fraction of their baryons, which results in a strong suppression of the SZ signal. Our simulations therefore predict a scaling with halo mass that is steeper than self-similar for haloes less massive than 10^13 Msun. For halo masses <~10^12 Msun, L_X is time-variable and correlated primarily with the star formation rate (SFR). For these objects, the diffuse X-ray emission is powered mostly by galactic winds and the gas dominating the X-ray emission is flowing out with radial velocities close to the halo's circular velocity. For halo masses >~10^13 Msun, on the other hand, L_X is much less variable and not correlated with the SFR, because the emission originates from the quasi-hydrostatic, virialized halo gas.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The equilibrium view on dust and metals in galaxies: Galactic outflows drive low dust-to-metal ratios in dwarf galaxies

Most galaxy evolution simulations as well as a variety of observational methods assume a linear scaling between the (galaxy-averaged) dust-to-gas ratio D and metallicity Z of the interstellar medium (ISM). Indeed, nearby galaxies with solar or moderately sub-solar metallicities clearly follow this trend albeit with significant scatter. However, a growing number of observations show that the linear scaling breaks down for metal-poor galaxies (Z<0.2 Z_sun), highlighting the need for a more sophisticated modeling of the dust-to-metal ratio of galaxies. Here we study the co-evolution of dust and metal abundances in galaxies with the help of a dynamical, one-zone model that incorporates dust formation and destruction processes in addition to gas inflows, outflows, and metal enrichment. The dynamical model is consistent with various observational constraints, including the stellar mass -- metallicity relation, the stellar mass -- halo mass relation, and the observed Z -- D relation for both metal-poor and metal-rich galaxies. The functional form of the Z -- D relation follows from a basic equilibrium ansatz, similar to the ideas used previously to model the stellar mass -- metallicity relation. Galactic outflows regulate the inflow rate of gas from the cosmic web for galaxies of a given star formation rate. The mass loading factor of outflows thus dictates the rate at which the dust and metal content of the ISM is diluted. The stellar mass dependence of the mass loading factor drives the evolution of metallicities, dust-to-gas ratios, and dust-to-metal ratios in galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Formation of Submillimetre-Bright Galaxies from Gas Infall over a Billion Years

Submillimetre-luminous galaxies at high-redshift are the most luminous, heavily star-forming galaxies in the Universe, and are characterised by prodigious emission in the far-infrared at 850 microns (S850 > 5 mJy). They reside in halos ~ 10^13Msun, have low gas fractions compared to main sequence disks at a comparable redshift, trace complex environments, and are not easily observable at optical wavelengths. Their physical origin remains unclear. Simulations have been able to form galaxies with the requisite luminosities, but have otherwise been unable to simultaneously match the stellar masses, star formation rates, gas fractions and environments. Here we report a cosmological hydrodynamic galaxy formation simulation that is able to form a submillimetre galaxy which simultaneously satisfies the broad range of observed physical constraints. We find that groups of galaxies residing in massive dark matter halos have rising star formation histories that peak at collective rates ~ 500-1000 Msun/yr at z=2-3, by which time the interstellar medium is sufficiently enriched with metals that the region may be observed as a submillimetre-selected system. The intense star formation rates are fueled in part by a reservoir gas supply enabled by stellar feedback at earlier times, not through major mergers. With a duty cycle of nearly a gigayear, our simulations show that the submillimetre-luminous phase of high-z galaxies is a drawn out one that is associated with significant mass buildup in early Universe proto-clusters, and that many submillimetre-luminous galaxies are actually composed of numerous unresolved components (for which there is some observational evidence).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Detecting Dark Matter Substructures around the Milky Way with Gaia

Cold Dark Matter (CDM) theory, a pillar of modern cosmology and astrophysics, predicts the existence of a large number of starless dark matter halos surrounding the Milky Way (MW). However, clear observational evidence of these "dark" substructures remains elusive. Here, we present a detection method based on the small, but detectable, velocity changes that an orbiting substructure imposes on the stars in the MW disk. Using high-resolution numerical simulations we estimate that the new space telescope Gaia should detect the kinematic signatures of a few starless substructures provided the CDM paradigm holds. Such a measurement will provide unprecedented constraints on the primordial matter power spectrum at low-mass scales and offer a new handle onto the particle physics properties of dark matter.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Argo Simulation: I. Quenching of Massive Galaxies at High Redshift as a Result of Cosmological Starvation

Observations show a prevalence of high redshift galaxies with large stellar masses and predominantly passive stellar populations. A variety of processes have been suggested that could reduce the star formation in such galaxies to observed levels, including quasar mode feedback, virial shock heating, or galactic winds driven by stellar feedback. However, the main quenching mechanisms have yet to be identified. Here we study the origin of star formation quenching using Argo, a cosmological, hydrodynamical zoom-in simulation that follows the evolution of a massive galaxy at $z\geq{}2$. This simulation adopts the same sub-grid recipes of the Eris simulations, which have been shown to form realistic disk galaxies, and, in one version, adopts also a mass and spatial resolution identical to Eris. The resulting galaxy has properties consistent with those of observed, massive (M_* ~ 1e11 M_sun) galaxies at z~2 and with abundance matching predictions. Our models do not include AGN feedback indicating that supermassive black holes likely play a subordinate role in determining masses and sizes of massive galaxies at high z. The specific star formation rate (sSFR) of the simulated galaxy matches the observed M_* - sSFR relation at early times. This period of smooth stellar mass growth comes to a sudden halt at z=3.5 when the sSFR drops by almost an order of magnitude within a few hundred Myr. The suppression is initiated by a leveling off and a subsequent reduction of the cool gas accretion rate onto the galaxy, and not by feedback processes. This "cosmological starvation" occurs as the parent dark matter halo switches from a fast collapsing mode to a slow accretion mode. Additional mechanisms, such as perhaps radio mode feedback from an AGN, are needed to quench any residual star formation of the galaxy and to maintain a low sSFR until the present time.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Argo Simulation: II. The Early Build-up of the Hubble Sequence

The Hubble sequence is a common classification scheme for the structure of galaxies. Despite the tremendous usefulness of this diagnostic, we still do not fully understand when, where, and how this morphological ordering was put in place. Here, we investigate the morphological evolution of a sample of 22 high redshift ($z\geq3$) galaxies extracted from the Argo simulation. Argo is a cosmological zoom-in simulation of a group-sized halo and its environment. It adopts the same high resolution ($\sim10^4$ M$_\odot$, $\sim100$ pc) and sub-grid physical model that was used in the Eris simulation but probes a sub-volume almost ten times bigger with as many as 45 million gas and star particles in the zoom-in region. Argo follows the early assembly of galaxies with a broad range of stellar masses ($\log M_{\star}/{\rm M}_{\odot}\sim8-11$ at $z\simeq3$), while resolving properly their structural properties. We recover a diversity of morphologies, including late-type/irregular disc galaxies with flat rotation curves, spheroid dominated early-type discs, and a massive elliptical galaxy, already established at $z\sim3$. We identify major mergers as the main trigger for the formation of bulges and the steepening of the circular velocity curves. Minor mergers and non-axisymmetric perturbations (stellar bars) drive the bulge growth in some cases. The specific angular momenta of the simulated disc components fairly match the values inferred from nearby galaxies of similar $M_{\star}$ once the expected redshift evolution of disc sizes is accounted for. We conclude that morphological transformations of high redshift galaxies of intermediate mass are likely triggered by processes similar to those at low redshift and result in an early build-up of the Hubble sequence.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Circum-Galactic Gas and the Isotropic Gamma Ray Background

Interactions of cosmic rays with the interstellar gas and radiation fields of the Milky Way provide the majority of the gamma rays observed by the Fermi Gamma Ray Space Telescope. In addition to the gas which is densely concentrated along the Galactic Disk, hydrodynamical simulations and observational evidence favor the presence of a halo of hot (T~10^6 K) ionized hydrogen (H_II), extending with non-negligible densities out to the virial radius of the Milky Way. We show that cosmic ray collisions with this circum-galactic gas should be expected to provide a significant flux of gamma rays, on the order of 10% of the observed isotopic gamma ray background at energies above 1 GeV. In addition, gamma rays originating from the extended H_II halos of other galaxies along a given line-of-sight should contribute to this background at a similar level.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The fate of high redshift massive compact galaxies in dense environments

Massive compact galaxies seem to be more common at high redshift than in the local universe, especially in denser environments. To investigate the fate of such massive galaxies identified at z~2 we analyse the evolution of their properties in three cosmological hydrodynamical simulations that form virialised galaxy groups of mass ~10^13 Msun hosting a central massive elliptical/S0 galaxy by redshift zero. We find that at redshift ~2 the population of galaxies with M_*> 2 10^10 Msun is diverse in terms of mass, velocity dispersion, star formation and effective radius, containing both very compact and relatively extended objects. In each simulation all the compact satellite galaxies have merged into the central galaxy by redshift 0 (with the exception of one simulation where one of such satellite galaxy survives). Satellites of similar mass at z = 0 are all less compact than their high redshift counterparts. They form later than the galaxies in the z = 2 sample and enter the group potential at z < 1, when dynamical friction times are longer than the Hubble time. Also, by z = 0 the central galaxies have increased substantially their characteristic radius via a combination of in situ star formation and mergers. Hence in a group environment descendants of compact galaxies either evolve towards larger sizes or they disappear before the present time as a result of the environment in which they evolve. Since the group-sized halos that we consider are representative of dense environments in the LambdaCDM cosmology, we conclude that the majority of high redshift compact massive galaxies do not survive until today as a result of the environment.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The relation between mid-plane pressure and molecular hydrogen in galaxies: Environmental dependence

Molecular hydrogen (H2) is the primary component of the reservoirs of cold, dense gas that fuel star formation in our galaxy. While the H2 abundance is ultimately regulated by physical processes operating on small scales in the interstellar medium (ISM), observations have revealed a tight correlation between the ratio of molecular to atomic hydrogen in nearby spiral galaxies and the pressure in the mid-plane of their disks. This empirical relation has been used to predict H2 abundances in galaxies with potentially very different ISM conditions, such as metal-deficient galaxies at high redshifts. Here, we test the validity of this approach by studying the dependence of the pressure -- H2 relation on environmental parameters of the ISM. To this end, we follow the formation and destruction of H2 explicitly in a suite of hydrodynamical simulations of galaxies with different ISM parameters. We find that a pressure -- H2 relation arises naturally in our simulations for a variety of dust-to-gas ratios or strengths of the interstellar radiation field in the ISM. Fixing the dust-to-gas ratio and the UV radiation field to values measured in the solar neighborhood results in fair agreement with the relation observed in nearby galaxies with roughly solar metallicity. However, the parameters (slope and normalization) of the pressure -- H2 relation vary in a systematical way with ISM properties. A particularly strong trend is the decrease of the normalization of the relation with a lowering of the dust-to-gas ratio of the ISM. We show that this trend and other properties of the pressure -- H2 relation are natural consequences of the transition from atomic to molecular hydrogen with gas surface density.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The X-factor in Galaxies: II. The molecular hydrogen -- star formation relation

There is ample observational evidence that the star formation rate (SFR) surface density, Sigma_SFR, is closely correlated with the surface density of molecular hydrogen, Sigma_H2. This empirical relation holds both for galaxy-wide averages and for individual >=kpc sized patches of the interstellar medium (ISM), but appears to degrade substantially at a sub-kpc scale. Identifying the physical mechanisms that determine the scale-dependent properties of the observed Sigma_H2-Sigma_SFR relation remains a challenge from a theoretical perspective. To address this question, we analyze the slope and scatter of the Sigma_H2-Sigma_SFR relation using a set of cosmological, galaxy formation simulations with a peak resolution of ~100 pc. These simulations include a chemical network for molecular hydrogen, a model for the CO emission, and a simple, stochastic prescription for star formation that operates on ~100 pc scales. Specifically, star formation is modeled as a Poisson process in which the average SFR is directly proportional to the present mass of H2. The predictions of our numerical model are in good agreement with the observed Kennicutt-Schmidt and Sigma_H2-Sigma_SFR relations. We show that observations based on CO emission are ill suited to reliably measure the slope of the latter relation at low (<20 M_sun pc^-2) H2 surface densities on sub-kpc scales. Our models also predict that the inferred Sigma_H2-Sigma_SFR relation steepens at high H2 surface densities as a result of the surface density dependence of the CO/H2 conversion factor. Finally, we show that on sub-kpc scales most of the scatter in the relation is a consequence of discreteness effects in the star formation process. In contrast, variations of the CO/H2 conversion factor are responsible for most of the scatter measured on super-kpc scales.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Intrinsic Alignment of Cluster Galaxies: the Redshift Evolution

We present measurements of two types of cluster galaxy alignments based on a volume limited and highly pure ($\ge$ 90%) sample of clusters from the GMBCG catalog derived from SDSS DR7. We detect a clear BCG alignment (the alignment of major axis of the BCG toward the distribution of cluster satellite galaxies). We find that the BCG alignment signal becomes stronger as the redshift and BCG absolute magnitude decrease, and becomes weaker as BCG stellar mass decreases. No dependence of the BCG alignment on cluster richness is found. We can detect a statistically significant ($\ge$ 3 sigma) satellite alignment (the alignment of the major axes of the cluster satellite galaxies toward the BCG) only when we use the isophotal fit position angles (PAs, hereafter), and the satellite alignment depends on the apparent magnitudes rather than the absolute magnitudes of the BCGs. This suggests the detected satellite alignment based on isophotoal PAs from the SDSS pipeline is possibly due to the contamination from the diffuse light of nearby BCGs. We caution that this should not be simply interpreted as non-existence of the satellite alignment, but rather that we cannot detect them with our current photometric SDSS data. We perform our measurements on both SDSS $r$ band and $i$ band data, but did not observe a passband dependence of the alignments.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The X-factor in Galaxies: I. Dependence on Environment and Scale

Characterizing the conversion factor between CO emission and column density of molecular hydrogen, X_CO, is crucial in studying the gaseous content of galaxies, its evolution, and relation to star formation. In most cases the conversion factor is assumed to be close to that of giant molecular clouds (GMCs) in the Milky Way, except possibly for mergers and star-bursting galaxies. However, there are physical grounds to expect that it should also depend on the gas metallicity, surface density, and strength of the interstellar radiation field. The XCO factor may also depend on the scale on which CO emission is averaged due to effects of limited resolution. We study the dependence of X_CO on gas properties and averaging scale using a model that is based on a combination of results of sub-pc scale magneto-hydrodynamic simulations and on the gas distribution from self-consistent cosmological simulations of galaxy formation. Our model predicts a value of X_CO that is consistent with the Galactic value for interstellar medium conditions typical for the Milky Way. For such conditions the predicted X_CO varies by only a factor of two for gas surfaced densities in the range \sim 50 - 500 M_sun / pc^2. However, the model also predicts that more generally on the scale of GMCs, X_CO is a strong function of metallicity, and depends on the column density and the interstellar UV flux. We show explicitly that neglecting these dependencies in observational estimates can strongly bias the inferred distribution of H2 column densities of molecular clouds to have a narrower and offset range compared to the true distribution. We find that when averaged on \sim kpc scales the X-factor depends only weakly on radiation field and column density, but is still a strong function of metallicity. The predicted metallicity dependence can be approximated as X_CO \sim Z^{-γ} with γ ~ 0.5 - 0.8.

0 citations0 reviewsNo code linkedOpen access