Source author record

Piero Madau

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

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA astro-ph astro-ph.HE hep-ph astro-ph.IM astro-ph.SR gr-qc

Catalog footprint

What is connected

73works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2026arXiv

Dynamical Dark Energy Imprints in the Lyman-Alpha Forest

The nature of dark energy (DE) remains elusive, even though it constitutes the dominant energy-density component of the Universe and drives the late-time acceleration of cosmic expansion. By combining measurements of the expansion history from baryon acoustic oscillations, supernova surveys, and cosmic microwave background data, the Dark Energy Spectroscopic Instrument (DESI) Collaboration has inferred that the DE equation of state may evolve over time. The profound implications of a time-variable, ``dynamical" DE (DDE) that departs from a cosmological constant motivate the need for independent observational tests. In this work, we use cosmological hydrodynamical simulations of structure formation to investigate how DDE affects the properties of the Lyman-Alpha ``forest'' of absorption features produced by neutral hydrogen in the cosmic web. We find that DDE models consistent with the DESI constraints induce a spectral tilt in the forest transmitted flux power spectrum, imprinting a scale- and redshift-dependent signature relative to standard Lambda-CDM cosmologies. These models also yield higher intergalactic medium temperatures and reduced Lyman-Alpha opacity compared to Lambda-CDM. We discuss the observational implications of these trends as potential avenues for independent confirmation of DDE.

preprint2022arXiv

Inferring the Thermal History of the Intergalactic Medium from the Properties of the Hydrogen and Helium Lyman-alpha Forest

The filamentary network of intergalactic medium (IGM) gas that gives origin to the Lyman-alpha forest in the spectra of distant quasars encodes information on the physics of structure formation and the early thermodynamics of diffuse baryonic material. Here, we use a massive suite of more than 400 high-resolution cosmological hydrodynamical simulations run with the Graphics Processing Unit-accelerated code Cholla to study the IGM at high spatial resolution maintained over the entire computational volume. The simulations capture a wide range of possible thermal histories of intergalactic gas by varying the amplitude and timing of the photoheating and photoionizing background produced by star-forming galaxies and active galactic nuclei. A statistical comparison of synthetic spectra with the observed 1D flux power spectra of hydrogen in 14 redshift bins over the full range 2.2 <= z <= 5.0 and with the Lyman-alpha opacity of helium in 5 redshift bins over the range 2.4 < z < 2.9 tightly constrains the photoionization and photoheating history of the IGM. By leveraging the constraining power of the available Lyman-alpha forest data to break model degeneracies, we find that the IGM experienced two main reheating events from the non-equilibrium ionization of hydrogen and helium over 1.2 Gyr of cosmic time. For our best-fit model, hydrogen reionization completes by z~6.0 with a first IGM temperature peak T_0 ~ 1.3x10^4 K, and is followed by the reionization of HeII that completes by z~3.0 and yields a second temperature peak of T_0 ~ 1.4x10^4 K. We discuss how our results can be used to obtain information on the timing and the sources of hydrogen and helium reionization.

preprint2022arXiv

The Dawn of Disk Formation in a Milky Way-sized Galaxy Halo: Thin Stellar Disks at $z > 4$

We present results from \textsc{GigaEris}, a cosmological, $N$-body hydrodynamical "zoom-in" simulation of the formation of a Milky Way-sized galaxy halo with unprecedented resolution, encompassing of order a billion particles within the refined region. The simulation employs a modern implementation of smoothed-particle hydrodynamics, including metal-line cooling and metal and thermal diffusion. We focus on the early assembly of the galaxy, down to redshift $z=4.4$. The simulated galaxy has properties consistent with extrapolations of the main sequence of star-forming galaxies to higher redshifts and levels off to a star formation rate of $\sim$60$\, M_{\odot}$~yr$^{-1}$ at $z=4.4$. A compact, thin rotating stellar disk with properties analogous to those of low-redshift systems arises already at $z \sim 8$. The galaxy rapidly develops a multi-component structure, and the disk, at least at these early stages, does not grow "upside-down" as often reported in the literature. Rather, at any given time, newly born stars contribute to sustain a thin disk. The kinematics reflect the early, ubiquitous presence of a thin disk, as a stellar disk component with $v_ϕ/σ_R$ larger than unity is already present at $z \sim 9$--10. Our results suggest that high-resolution spectro-photometric observations of very high-redshift galaxies should find thin rotating disks, consistent with the recent discovery of cold rotating gas disks by ALMA. Finally, we present synthetic images for the JWST NIRCam camera, showing how the early disk would be easily detectable already at those early times.

preprint2021arXiv

Deep Realistic Extragalactic Model (DREaM) Galaxy Catalogs: Predictions for a Roman Ultra-Deep Field

In the next decade, deep galaxy surveys from telescopes such as the James Webb Space Telescope and Roman Space Telescope will provide transformational data sets that will greatly enhance the understanding of galaxy formation during the epoch of reionization (EoR). In this work, we present the Deep Realistic Extragalactic Model (DREaM) for creating synthetic galaxy catalogs. Our model combines dark matter simulations, subhalo abundance matching and empirical models, and includes galaxy positions, morphologies, and spectral energy distributions (SEDs). The resulting synthetic catalog extends to redshifts $z \sim 12 $, and galaxy masses $\log_{10}(M/M_{\odot}) = 5 $ covering an area of $1\, {\rm deg}^2$ on the sky. We use DREaM to explore the science returns of a $1\, {\rm deg}^2$ \emph{Roman} UDF, and to provide a resource for optimizing ultra-deep survey designs. We find that a \emph{Roman} UDF to $\sim 30\,m_{\rm AB}$ will potentially detect more than $10^6$ $M_{\rm UV}<-17$ galaxies, with more than $10^4$ at redshifts $z>7$, offering an unparalleled dataset for constraining galaxy properties during the EoR. Our synthetic catalogs and simulated images are made publicly available to provide the community with a tool to prepare for upcoming data.

preprint2020arXiv

Formation of the first stars and black holes

We review the current status of knowledge concerning the early phases of star formation during cosmic dawn. This includes the first generations of stars forming in the lowest mass dark matter halos in which cooling and condensation of gas with primordial composition is possible at very high redshift ($z > 20$), namely metal-free Population III stars, and the first generation of massive black holes forming at such early epochs, the so-called black hole seeds. The formation of black hole seeds as end states of the collapse of Population III stars, or via direct collapse scenarios, is discussed. In particular, special emphasis is given to the physics of supermassive stars as potential precursors of direct collapse black holes, in light of recent results of stellar evolution models, and of numerical simulations of the early stages of galaxy formation. Furthermore, we discuss the role of the cosmic radiation produced by the early generation of stars and black holes at high redshift in the process of reionization.

preprint2020arXiv

Global torques and stochasticity as the drivers of massive black hole pairing in the young Universe

The forthcoming Laser Interferometer Space Antenna (LISA) will probe the population of coalescing massive black hole (MBH) binaries up to the onset of structure formation. Here we simulate the galactic-scale pairing of $\sim10^6 M_\odot$ MBHs in a typical, non-clumpy main-sequence galaxy embedded in a cosmological environment at $z = 7-6$. In order to increase our statistical sample, we adopt a strategy that allows us to follow the evolution of six secondary MBHs concomitantly. We find that the magnitude of the dynamical-friction induced torques is significantly smaller than that of the large-scale, stochastic gravitational torques arising from the perturbed and morphologically evolving galactic disc, suggesting that the standard dynamical friction treatment is inadequate for realistic galaxies at high redshift. The dynamical evolution of MBHs is very stochastic, and a variation in the initial orbital phase can lead to a drastically different time-scale for the inspiral. Most remarkably, the development of a galactic bar in the host system either significantly accelerates the inspiral by dragging a secondary MBH into the centre, or ultimately hinders the orbital decay by scattering the MBH in the galaxy outskirts. The latter occurs more rarely, suggesting that galactic bars overall promote MBH inspiral and binary coalescence. The orbital decay time can be an order of magnitude shorter than what would be predicted relying on dynamical friction alone. The stochasticity, and the important role of global torques, have crucial implications for the rates of MBH coalescences in the early Universe: both have to be accounted for when making predictions for the upcoming LISA observatory.

preprint2020arXiv

Globular Cluster Formation from Colliding Substructure

We investigate a scenario where the formation of Globular Clusters (GCs) is triggered by high-speed collisions between infalling atomic-cooling subhalos during the assembly of the main galaxy host, a special dynamical mode of star formation that operates at high gas pressures and is intimately tied to LCDM hierarchical galaxy assembly. The proposed mechanism would give origin to "naked" globulars, as colliding dark matter subhalos and their stars will simply pass through one another while the warm gas within them clashes at highly supersonic speed and decouples from the collisionless component, in a process reminiscent of the Bullet galaxy cluster. We find that the resulting shock-compressed layer cools on a timescale that is typically shorter than the crossing time, first by atomic line emission and then via fine-structure metal-line emission, and is subject to gravitational instability and fragmentation. Through a combination of kinetic theory approximation and high-resolution $N$-body simulations, we show that this model may produce: (a) a GC number-halo mass relation that is linear down to dwarf galaxy scales and agrees with the trend observed over five orders of magnitude in galaxy mass; (b) a population of old globulars with a median age of 12 Gyr and an age spread similar to that observed; (c) a spatial distribution that is biased relative to the overall mass profile of the host; and (d) a bimodal metallicity distribution with a spread similar to that observed in massive galaxies.

preprint2020arXiv

The impact of Lyman-$α$ emission line heating and cooling on the cosmic dawn 21-cm signal

Allowing for enhanced Ly$α$ photon line emission from Population III dominated stellar systems in the first forming galaxies, we show the 21-cm cosmic dawn signal at $10<z<30$ may substantially differ from standard scenarios. Energy transfer by Ly$α$ photons emerging from galaxies may heat intergalactic gas if HII regions within galaxies are recombination bound, or cool the gas faster than by adiabatic expansion if reddened by winds internal to the haloes. In some cases, differential 21-cm antenna temperatures near $-500$ mK may be achieved at $15<z<25$, similar to the signature detected by the EDGES 21-cm cosmic dawn experiment.

preprint2019arXiv

Momentum Injection by Clustered Supernovae: Testing Subgrid Feedback Prescriptions

Using a 1D Lagrangian code specifically designed to assess the impact of multiple, time-resolved supernovae (SNe) from a single star cluster on the surrounding medium, we test three commonly used feedback recipes: delayed cooling (e.g., used in the GASOLINE-2 code), momentum-energy injection (a resolution-dependent transition between momentum-dominated feedback and energy-dominated feedback used, e.g., in the FIRE-2 code), and simultaneous energy injection (e.g., used in the EAGLE simulations). Our work provides an intermediary test for these recipes: we analyse a setting that is more complex than the simplified scenarios for which many were designed, but one more controlled than a full galactic simulation. In particular, we test how well these models reproduce the enhanced momentum efficiency seen for an 11 SN cluster simulated at high resolution (0.6 pc; a factor of 12 enhancement relative to the isolated SN case) when these subgrid recipes are implemented in low resolution (20 pc) runs. We find that: 1) the delayed cooling model performs well -- resulting in 9 times the momentum efficiency of the fiducial isolated SN value -- when SNe are clustered and $10^{51}$ erg are injected per SN, while clearly over-predicting the momentum efficiency in the single SN test case; 2) the momentum-energy model always achieves good results, with a factor of 5 boost in momentum efficiency; and 3) injecting the energy from all SNe simultaneously does little to prevent over-cooling and greatly under-produces the momentum deposited by clustered SNe, resulting in a factor of 3 decrease in momentum efficiency on the average.

preprint2016arXiv

Bar-driven evolution and quenching of spiral galaxies in cosmological simulations

We analyse the output of the hi-res cosmological zoom-in simulation ErisBH to study self-consistently the formation of a strong stellar bar in a Milky Way-type galaxy and its effect on the galactic structure, on the central gas distribution and on star formation. The simulation includes radiative cooling, star formation, SN feedback and a central massive black hole which is undergoing gas accretion and is heating the surroundings via thermal AGN feedback. A large central region in the ErisBH disk becomes bar-unstable after z~1.4, but a clear bar-like structure starts to grow significantly only after z~0.4, possibly triggered by the interaction with a massive satellite. At z~0.1 the bar reaches its maximum radial extent of l~2.2 kpc. As the bar grows, it becomes prone to buckling instability, which we quantify based on the anisotropy of the stellar velocity dispersion. The actual buckling event is observable at z~0.1, resulting in the formation of a boxy-peanut bulge clearly discernible in the edge-on view of the galaxy at z=0. The bar in ErisBH does not dissolve during the formation of the bulge but remains strongly non-axisymmetric down to the resolution limit of ~100 pc at z=0. During its early growth, the bar exerts a strong torque on the gas within its extent and drives gas inflows that enhance the nuclear star formation on sub-kpc scales. Later on the infalling gas is nearly all consumed into stars and, to a lesser extent, accreted onto the central black hole, leaving behind a gas-depleted region within the central ~2 kpc. Observations would more likely identify a prominent, large-scale bar at the stage when the galactic central region has already been quenched. Bar-driven quenching may play an important role in disk-dominated galaxies at all redshift. [Abridged]

preprint2016arXiv

Clumpy high-z galaxies as a testbed for feedback-regulated galaxy formation

We study the dependence of fragmentation in massive gas-rich galaxy disks at $z > 1$ on feedback model and hydrodynamical method, employing the GASOLINE2 SPH code and the lagrangian mesh-less code GIZMO in finite mass mode. We compare non-cosmological galaxy disk runs with standard blastwave supernovae (SN)feedback, which introduces delayed cooling in order to drive winds, and runs with the new superbubble SN feedback, which produces winds naturally by modelling the detailed physics of SN-driven bubbles and leads to efficient self-regulation of star formation. We find that, with blastwave feedback, massive star forming clumps form in comparable number and with very similar masses in GASOLINE2 and GIZMO. The typical masses are in the range $10^7-10^8 M_{\odot}$, lower than in most previous works, while giant clumps with masses above $10^9 M_{\odot}$ are exceedingly rare. With superbubble feedback, instead, massive bound star forming clumps do not form because galaxies never undergo a phase of violent disk instability. Only sporadic, unbound star forming overdensities lasting only a few tens of Myr can arise that are triggered by perturbations of massive satellite companions. We conclude that there is a severe tension between explaining massive star forming clumps observed at $z > 1$ primarily as the result of disk fragmentation driven by gravitational instability and the prevailing view of feedback-regulated galaxy formation. The link between disk stability and star formation efficiency should thus be regarded as a key testing ground for galaxy formation theory.

preprint2016arXiv

Cold Dark Matter Substructures in Early-Type Galaxy Halos

We present initial results from the "Ponos" zoom-in numerical simulations of dark matter substructures in massive ellipticals. Two very highly resolved dark matter halos with $M_{\rm vir}=1.2\times 10^{13}$ $M_{\odot}$ and $M_{\rm vir}=6.5\times 10^{12}$ $M_{\odot}$ and different ("violent" vs. "quiescent") assembly histories have been simulated down to $z=0$ in a $Λ$CDM cosmology with a total of 921,651,914 and 408,377,544 particles, respectively. Within the virial radius, the total mass fraction in self-bound $M_{\rm sub}>10^6$ $M_{\odot}$ subhalos at the present epoch is 15% for the violent host and 16.5% for the quiescent one. At $z=0.7$, these fractions increase to 19 and 33%, respectively, as more recently accreted satellites are less prone to tidal destruction. In projection, the average fraction of surface mass density in substructure at a distance of $R/R_{\rm vir}=0.02$ ($\sim 5-10$ kpc) from the two halo centers ranges from 0.6% to $\gtrsim 2$%, significantly higher than measured in simulations of Milky Way-sized halos. The contribution of subhalos with $M_{\rm sub} < 10^9$ $M_{\odot}$ to the projected mass fraction is between one fifth and one third of the total, with the smallest share found in the quiescent host. We assess the impact of baryonic effects via twin, lower-resolution hydrodynamical simulations that include metallicity-dependent gas cooling, star formation, and a delayed-radiative-cooling scheme for supernova feedback. Baryonic contraction produces a super-isothermal total density profile and increases the number of massive subhalos in the inner regions of the main host. The host density profiles and projected subhalo mass fractions appear to be broadly consistent with observations of gravitational lenses.

preprint2016arXiv

Diffuse coronae in cosmological simulations of Milky Way-sized galaxies

We investigate the properties of halo gas using three cosmological `zoom-in' simulations of realistic Milky Way-galaxy analogs with varying sub-grid physics. In all three cases, the mass of hot ($T > 10^6$ K) halo gas is $\sim 1\%$ of the host's virial mass. The X-ray luminosity of two of the runs is consistent with observations of the Milky Way, while the third simulation is X-ray bright and resembles more closely a very massive, star-forming spiral. Hot halos extend to 140 kpc from the galactic center and are surrounded by a bubble of warm-hot ($T = 10^5 - 10^6$K) gas that extends to the virial radius. Simulated halos agree well outside 20-30 kpc with the $β$-model of Miller & Bregman (2013) based on OVII absorption and OVIII emission measurements. Warm--hot and hot gas contribute up to $80\%$ of the total gas reservoir, and contain nearly the same amount of baryons as the stellar component. The mass of warm-hot and hot components falls into the range estimated for $L^*$ galaxies. With key observational constraints on the density of the Milky Way corona being satisfied, we show that concealing of the ubiquitous warm-hot baryons, along with the ejection of just $20-30 \%$ of the diffuse gas out of the potential wells by supernova-driven outflows, can solve the "missing baryon problem". The recovered baryon fraction within 3 virial radii is $90\%$ of the universal value. With a characteristic density of $\sim 10^{-4}$ cm$^{-3}$ at $50-80$ kpc, diffuse coronae meet the requirement for fast and complete ram--pressure stripping of the gas reservoirs in dwarf galaxy satellites.

preprint2016arXiv

Dispersal of tidal debris in a Milky-Way-sized dark matter halo

We simulate the tidal disruption of a collisionless N-body globular star cluster in a total of 300 different orbits selected to have galactocentric radii between 10 and 30 kpc in four dark matter halos: (a) a spherical halo with no subhalos, (b) a spherical halo with subhalos, (c) a realistic halo with no subhalos, and (d) a realistic halo with subhalos. This allows us to isolate and study how the halo's (lack of) dynamical symmetry and substructures affect the dispersal of tidal debris. The realistic halos are constructed from the snapshot of the Via Lactea II simulation at redshift zero. We find that the overall halo's lack of dynamical symmetry disperses tidal debris to make the streams fluffier, consistent with previous studies of tidal debris of dwarf galaxies in larger orbits than ours in this study. On the other hand, subhalos in realistic potentials can locally enhance the densities along streams, making streams denser than their counterparts in smooth potentials. We show that many long and thin streams can survive in a realistic and lumpy halo for a Hubble time. This suggests that upcoming stellar surveys will likely uncover more thin streams which may contain density gaps that have been shown to be promising probes for dark matter substructures.

preprint2016arXiv

Enhanced Momentum Feedback from Clustered Supernovae

Young stars typically form in star clusters, so the supernovae (SNe) they produce are clustered in space and time. This clustering of SNe may alter the momentum per SN deposited in the interstellar medium (ISM) by affecting the local ISM density, which in turn affects the cooling rate. We study the effect of multiple SNe using idealized 1D hydrodynamic simulations which explore a large parameter space of the number of SNe, and the background gas density and metallicity. The results are provided as a table and an analytic fitting formula. We find that for clusters with up to ~100 SNe the asymptotic momentum scales super-linearly with the number of SNe, resulting in a momentum per SN that can be an order of magnitude larger than for a single SN, with a maximum efficiency for clusters with 10-100 SNe. We argue that additional physical processes not included in our simulations -- self-gravity, breakout from a galactic disk, and galactic shear -- can slightly reduce the momentum enhancement from clustering, but the average momentum per SN still remains a factor of 4 larger than the isolated SN value when averaged over a realistic cluster mass function for a star-forming galaxy. We conclude with a discussion of the possible role of mixing between hot and cold gas, induced by multi-dimensional instabilities or preexisting density variations, as a limiting factor in the buildup of momentum by clustered SNe, and suggest future numerical experiments to explore these effects.

preprint2015arXiv

Black Hole Starvation and Bulge Evolution in a Milky Way-like Galaxy

We present a new zoom-in hydrodynamical simulation, "Erisbh", which follows the cosmological evolution and feedback effects of a supermassive black hole at the center of a Milky Way-type galaxy. ErisBH shares the same initial conditions, resolution, recipes of gas cooling, star formation and feedback, as the close Milky Way-analog "Eris", but it also includes prescriptions for the formation, growth and feedback of supermassive black holes. We find that the galaxy's central black hole grows mainly through mergers with other black holes coming from infalling satellite galaxies. The growth by gas accretion is minimal because very little gas reaches the sub-kiloparsec scales. The final black hole is, at z=0, about 2.6 million solar masses and it sits closely to the position of SgrA* on the MBH-MBulge and MBH-sigma planes, in a location consistent with what observed for pseudobulges. Given the limited growth due to gas accretion, we argue that the mass of the central black hole should be above 10^5 solar masses already at z~8. The effect of AGN feedback on the host galaxy is limited to the very central few hundreds of parsecs. Despite being weak, AGN feedback seems to be responsible for the limited growth of the central bulge with respect to the original Eris, which results in a significantly flatter rotation curve in the inner few kiloparsecs. Moreover, the disk of ErisBH is more prone to instabilities, as its bulge is smaller and its disk larger then Eris. As a result, the disk of ErisBH undergoes a stronger dynamical evolution relative to Eris and around z=0.3 a weak bar grows into a strong bar of a few disk scale lengths in size. The bar triggers a burst of star formation in the inner few hundred parsecs, provides a modest amount of new fuel to the central black hole, and causes the bulge of ErisBH to have, by z=0, a box/peanut morphology.(Abridged)

preprint2015arXiv

Cosmic Reionization after Planck: Could Quasars Do It All?

We assess a model of late cosmic reionization in which the ionizing background radiation arises entirely from high redshift quasars and other active galactic nuclei (AGNs). The low optical depth to Thomson scattering reported by the Planck Collaboration pushes the redshift of instantaneous reionization down to z=8.8^{+1.7}_{-1.4} and greatly reduces the need for significant Lyman-continuum emission at very early times. We show that, if recent claims of a numerous population of faint AGNs at z=4-6 are upheld, and the high inferred AGN comoving emissivity at these epochs persists to higher, z~10, redshifts, then active galaxies may drive the reionization of hydrogen and helium with little contribution from normal star-forming galaxies. We discuss an AGN-dominated scenario that satisfies a number of observational constraints: the HI photoionization rate is relatively flat over the range 2<z<5, hydrogen gets fully reionized by z=5.7, and the integrated Thomson scattering optical depth is tau=0.056, in agreement with measurements based on the Lya opacity of the intergalactic medium (IGM) and cosmic microwave background (CMB) polarization. It is a prediction of the model that helium gets doubly reionized before redshift 4, the heat input from helium reionization dominates the thermal balance of the IGM after hydrogen reionization, and z>5 AGNs provide a significant fraction of the unresolved X-ray background at 2 keV. Singly- and doubly-ionized helium contribute about 13% to tau, and the HeIII volume fraction is already 50% when hydrogen becomes fully reionized.

preprint2015arXiv

Evidence of patchy hydrogen reionization from an extreme Ly$α$ trough below redshift six

We report the discovery of an extremely long ($\sim$110 Mpc/$h$) and dark ($τ_{\rm eff} \gtrsim 7$) Ly$α$ trough extending down to $z \simeq 5.5$ towards the $z_{\rm em} \simeq 6.0$ quasar ULAS J0148+0600. We use these new data in combination with Ly$α$ forest measurements from 42 quasars at $4.5 \le z_{\rm em} \le 6.4$ to conduct an updated analysis of the line-of-sight variance in the intergalactic Ly$α$ opacity over $4 \le z \le 6$. We find that the scatter in transmission among lines of sight near $z \sim 6$ significantly exceeds theoretical expectations for either a uniform ultraviolet background (UVB) or simple fluctuating UVB models in which the mean free path to ionizing photons is spatially invariant. The data, particularly near $z \simeq 5.6$-5.8, instead require fluctuations in the volume-weighted hydrogen neutral fraction that are a factor of 3 or more beyond those expected from density variations alone. We argue that these fluctuations are most likely driven by large-scale variations in the mean free path, consistent with expectations for the final stages of inhomogeneous hydrogen reionization. Even by $z \simeq 5.6$, however, a large fraction of the data are consistent with a uniform UVB, and by $z \sim 5$ the data are fully consistent with opacity fluctuations arising solely from the density field. This suggests that while reionization may be ongoing at $z \sim 6$, it has fully completed by $z \sim 5$.

preprint2015arXiv

Signatures of Kinematic Substructure in the Galactic Stellar Halo

Tidal debris from infalling satellites can leave observable structure in the phase-space distribution of the Galactic halo. Such substructure can be manifest in the spatial and/or velocity distributions of the stars in the halo. This paper focuses on a class of substructure that is purely kinematic in nature, with no accompanying spatial features. To study its properties, we use a simulated stellar halo created by dynamically populating the Via Lactea II high-resolution N-body simulation with stars. A significant fraction of the stars in the inner halo of Via Lactea share a common speed and metallicity, despite the fact that they are spatially diffuse. We argue that this kinematic substructure is a generic feature of tidal debris from older mergers and may explain the detection of radial-velocity substructure in the inner halo made by the Sloan Extension for Galactic Understanding and Exploration. The GAIA satellite, which will provide the proper motions of an unprecedented number of stars, should further characterize the kinematic substructure in the inner halo. Our study of the Via Lactea simulation suggests that the stellar halo can be used to map the speed distribution of the local dark-matter halo, which has important consequences for dark-matter direct-detection experiments.

preprint2015arXiv

Simulating tidal streams in a high resolution dark matter halo

We simulate tidal streams in the presence and absence of substructures inside the zero redshift snapshot of the Via Lactea II (VL-2) simulation. A halo finder is used to remove and isolate the subhalos found inside the high resolution dark matter halo of VL-2, and the potentials for both the main halo and all the subhalos are constructed individually using the self-consistent field (SCF) method. This allows us to make direct comparison of tidal streams between a smooth halo and a lumpy halo without assuming idealized profiles or triaxial fits. We simulate the kinematics of a star cluster starting with the same orbital position but two different velocities. Although these two orbits are only moderately eccentric and have similar apo- and pericentric distances, we find that the two streams have very different morphologies. We conclude that our model of the potential of VL-2 can provide insights about tidal streams that have not been explored by previous studies using idealized or axisymmetric models.

preprint2015arXiv

The History of R-Process Enrichment in the Milky Way

We investigate the production sites and the enrichment history of $r$-process elements in the Galaxy, as traced by the [Eu/Fe] ratio, using the high resolution, cosmological zoom-in simulation `Eris'. At $z=0$, Eris represents a close analog to the Milky Way, making it the ideal laboratory to understand the chemical evolution of our Galaxy. Eris formally traces the production of oxygen and iron due to Type-Ia and Type-II supernovae. We include in post-processing the production of $r$-process elements from compact binary mergers. Unlike previous studies, we find that the nucleosynthetic products from compact binary mergers can be incorporated into stars of very low metallicity and at early times, even with a minimum delay time of 100 Myr. This conclusion is relatively insensitive to modest variations in the merger rate, minimum delay time, and the delay time distribution. By implementing a first-order prescription for metal-mixing, we can further improve the agreement between our model and the data for the chemical evolution of both [$α$/Fe] and [Eu/Fe].We argue that compact binary mergers could be the dominant source of $r$-process nucleosynthesis in the Galaxy.

preprint2015arXiv

The photoheating of the intergalactic medium in synthesis models of the UV background

We compare cosmological hydrodynamical simulations combined with the homogeneous metagalactic UV background (UVB) of Haardt & Madau (2012) (HM2012) to observations of the Lyman-alpha forest that are sensitive to the thermal and ionization state of the intergalactic medium (IGM). The transition from optically thick to thin photoheating predicted by the simple one-zone, radiative transfer model implemented by HM2012 predicts a thermal history that is in remarkably good agreement with the observed rise of the IGM temperature at z~3 if we account for the expected evolution of the volume filling factor of HeIII. Our simulations indicate that there may be, however, some tension between the observed peak in the temperature evolution and the rather slow evolution of the HeII opacities suggested by recent Hubble Space Telescope/COS measurements. The HM2012 UVB also underpredicts the metagalactic hydrogen photoionization rate required by our simulations to match the observed opacity of the forest at z>4 and z<2.

preprint2014arXiv

A cosmic web filament revealed in Lyman-alpha emission around a luminous high-redshift quasar

Simulations of structure formation in the Universe predict that galaxies are embedded in a "cosmic web", where the majority of baryons reside as rarefied and highly ionized gas. This material has been studied for decades in absorption against background sources, but the sparseness of these inherently one-dimensional probes preclude direct constraints on the three-dimensional morphology of the underlying web. Here we report observations of a cosmic web filament in Lyman-alpha emission, discovered during a survey for cosmic gas fluorescently "illuminated" by bright quasars at z=2.3. With a projected size of approximately 460 physical kpc, the Lyman-alpha emission surrounding the radio-quiet quasar UM287 extends well beyond the virial radius of any plausible associated dark matter halo. The estimated cold gas mass of the nebula from the observed emission is at least ten times larger than what is typically found by cosmological simulations, suggesting that a population of intergalactic gas clumps with sub-kpc sizes may be missing within current numerical models.

preprint2014arXiv

Building Late-Type Spiral Galaxies by In-Situ and Ex-Situ Star Formation

We analyze the formation and evolution of the stellar components in "Eris", a 120 pc-resolution cosmological hydrodynamic simulation of a late-type spiral galaxy. The simulation includes the effects of a uniform UV background, a delayed-radiative-cooling scheme for supernova feedback, and a star formation recipe based on a high gas density threshold. It allows a detailed study of the relative contributions of "in-situ" (within the main host) and "ex-situ" (within satellite galaxies) star formation to each major Galactic component in a close Milky Way analog. We investigate these two star-formation channels as a function of galactocentric distance, along different lines of sight above and along the disk plane, and as a function of cosmic time. We find that: 1) approximately 70 percent of today's stars formed in-situ; 2) more than two thirds of the ex-situ stars formed within satellites after infall; 3) the majority of ex-situ stars are found today in the disk and in the bulge; 4) the stellar halo is dominated by ex-situ stars, whereas in-situ stars dominate the mass profile at distances < 5 kpc from the center at high latitudes; and 5) approximately 25% of the inner, r < 20 kpc, halo is composed of in-situ stars that have been displaced from their original birth sites during Eris' early assembly history.

preprint2014arXiv

Carbon-Enhanced Metal-Poor Stars: Relics from the Dark Ages

We use detailed nucleosynthesis calculations and a realistic prescription for the environment of the first stars to explore the first episodes of chemical enrichment that occurred during the dark ages. Based on these calculations, we propose a novel explanation for the increased prevalence of carbon-enhanced metal-poor (CEMP) stars with decreasing Fe abundance: The observed chemistry for the most metal-poor Galactic halo stars is the result of an intimate link between the explosions of the first stars and their host minihalo's ability to retain its gas. Specifically, high-energy supernovae produce a near solar ratio of C/Fe, but are effective in evacuating the gas from their host minihalo, thereby suppressing the formation of a second generation of stars. On the other hand, minihalos that host low-energy supernovae are able to retain their gas and form a second stellar generation but, as a result, the second stars are born with a supersolar ratio of C/Fe. Our models are able to accurately reproduce the observed distributions of [C/Fe] and [Fe/H], as well as the fraction of CEMP stars relative to non-CEMP stars as a function of [Fe/H] without any free parameters. We propose that the present lack of chemical evidence for very massive stars (>140 Msun), that ended their lives as a highly energetic pair-instability supernova, does not imply that such stars were rare or did not exist; the chemical products of these very massive first stars may have escaped from their host minihalo, and were never incorporated into subsequent generations of stars. Finally, our models suggest that the most Fe-poor stars currently known may have seen the enrichment from a small multiple of metal-free stars, and need not have been exclusively enriched by a solitary first star. These calculations support the idea that some of the surviving dwarf satellite galaxies of the Milky Way are relics of the first galaxies.

preprint2014arXiv

Cosmic Star Formation History

Over the past two decades, an avalanche of data from multiwavelength imaging and spectroscopic surveys has revolutionized our view of galaxy formation and evolution. Here we review the range of complementary techniques and theoretical tools that allow astronomers to map the cosmic history of star formation, heavy element production, and reionization of the Universe from the cosmic "dark ages" to the present epoch. A consistent picture is emerging, whereby the star-formation rate density peaked approximately 3.5 Gyr after the Big Bang, at z~1.9, and declined exponentially at later times, with an e-folding timescale of 3.9 Gyr. Half of the stellar mass observed today was formed before a redshift z = 1.3. About 25% formed before the peak of the cosmic star-formation rate density, and another 25% formed after z = 0.7. Less than ~1% of today's stars formed during the epoch of reionization. Under the assumption of a universal initial mass function, the global stellar mass density inferred at any epoch matches reasonably well the time integral of all the preceding star-formation activity. The comoving rates of star formation and central black hole accretion follow a similar rise and fall, offering evidence for co-evolution of black holes and their host galaxies. The rise of the mean metallicity of the Universe to about 0.001 solar by z = 6, one Gyr after the Big Bang, appears to have been accompanied by the production of fewer than ten hydrogen Lyman-continuum photons per baryon, a rather tight budget for cosmological reionization.

preprint2014arXiv

Dark Matter Heating and Early Core Formation in Dwarf Galaxies

We present more results from a fully cosmological LCDM simulation of a group of isolated dwarf galaxies that has been shown to reproduce the observed stellar mass and cold gas content, resolved star formation histories, and metallicities of dwarfs in the Local Volume. Here we investigate the energetics and timetable of the cusp-core transformation. As suggested by previous work, supernova-driven gas outflows remove dark matter (DM) cusps and create kpc-size cores in all systems having a stellar mass above 1e6 Msun. The "DM core mass removal efficiency" -- dark mass ejected per unit stellar mass -- ranges today from a few to a dozen, and increases with decreasing host mass. Because dwarfs form the bulk of their stars prior to redshift 1 and the amount of work required for DM heating and core formation scales approximately as Mvir^{5/3}, the unbinding of the DM cusp starts early and the formation of cored profiles is not as energetically onerous as previously claimed. DM particles in the cusp typically migrate to 2-3 core radii after absorbing a few percent of the energy released by supernovae. The present-day slopes of the inner dark matter mass profiles, Gamma=d log M/dlog R=2.5-3, of the simulated "Bashful" and "Doc" dwarfs are similar to those measured in the luminous Fornax and Sculptor dwarf spheroidals. None of the simulated galaxies has a circular velocity profile exceeding 20 km/s in the inner 1 kpc, implying that supernova feedback is key to solve the "too-big-to-fail" problem for Milky Way subhalos.

preprint2014arXiv

Faint dwarfs as a test of DM models: WDM vs. CDM

We use high resolution Hydro$+$N-Body cosmological simulations to compare the assembly and evolution of a small field dwarf (stellar mass ~ 10$^{6-7}$ M$\odot$, total mass 10$^{10}$ M$\odot$ in $Λ$ dominated CDM and 2keV WDM cosmologies. We find that star formation (SF) in the WDM model is reduced and delayed by 1-2 Gyr relative to the CDM model, independently of the details of SF and feedback. Independent of the DM model, but proportionally to the SF efficiency, gas outflows lower the central mass density through `dynamical heating', such that all realizations have circular velocities $<$ 20kms at 500$~$pc, in agreement with local kinematic constraints. As a result of dynamical heating, older stars are less centrally concentrated than younger stars, similar to stellar population gradients observed in nearby dwarf galaxies. Introducing an important diagnostic of SF and feedback models, we translate our simulations into artificial color-magnitude diagrams and star formation histories in order to directly compare to available observations. The simulated galaxies formed most of their stars in many $\sim$10 Myr long bursts. The CDM galaxy has a global SFH, HI abundance and Fe/H and alpha-elements distribution well matched to current observations of dwarf galaxies. These results highlight the importance of directly including `baryon physics' in simulations when 1) comparing predictions of galaxy formation models with the kinematics and number density of local dwarf galaxies and 2) differentiating between CDM and non-standard models with different DM or power spectra.

preprint2014arXiv

Reversal of Fortune: Increased Star Formation Efficiencies in the Early Histories of Dwarf Galaxies?

On dwarf galaxy scales, the different shapes of the galaxy stellar mass function and the dark halo mass function require a star-formation efficiency (SFE) in these systems that is currently more than 1 dex lower than that of Milky Way-size halos. Here, we argue that this trend may actually be reversed at high redshift. Specifically, by combining the resolved star-formation histories of nearby isolated dwarfs with the simulated mass-growth rates of dark matter halos, we show that the assembly of these systems occurs in two phases: (1) an early, fast halo accretion phase with a rapidly deepening potential well, characterized by a high SFE; and (2) a late slow halo accretion phase where, perhaps as a consequence of reionization, the SFE is low. Nearby dwarfs have more old stars than predicted by assuming a constant or decreasing SFE with redshift, a behavior that appears to deviate qualitatively from the trends seen amongst more massive systems. Taken at face value, the data suggest that, at sufficiently early epochs, dwarf galaxy halos above the atomic cooling mass limit can be among the most efficient sites of star formation in the universe.

preprint2014arXiv

Super-Critical Growth of Massive Black Holes from Stellar-Mass Seeds

We consider super-critical accretion with angular momentum onto stellar-mass black holes as a possible mechanism for growing billion-solar-mass holes from light seeds at early times. We use the radiatively-inefficient "slim disk" solution -- advective, optically thick flows that generalize the standard geometrically thin disk model -- to show how mildly super-Eddington intermittent accretion may significantly ease the problem of assembling the first massive black holes when the Universe was less than 0.8 Gyr old. Because of the low radiative efficiencies of slim disks around non-rotating as well as rapidly rotating holes, the mass e-folding timescale in this regime is nearly independent of the spin parameter. The conditions that may lead to super-critical growth in the early Universe are briefly discussed.

preprint2014arXiv

The Baryon Cycle of Dwarf Galaxies: Dark, Bursty, Gas-Rich Polluters

We present results from a fully cosmological, very high-resolution, LCDM "zoom-in" simulation of a group of seven field dwarf galaxies with present-day virial masses in the range M_vir=4.4e8-3.6e10 Msun. The simulation includes a blastwave scheme for supernova feedback, a star formation recipe based on a high gas density threshold, metal-dependent radiative cooling, a scheme for the turbulent diffusion of metals and thermal energy, and a uniform UV background that modifies the ionization and excitation state of the gas. The properties of the simulated dwarfs are strongly modulated by the depth of the gravitational potential well. All three halos with M_vir < 1e9 Msun are devoid of stars, as they never reach the density threshold for star formation of 100 atoms/cc. The other four, M_vir > 1e9 Msun dwarfs have blue colors, low star formation efficiencies, high cold gas to stellar mass ratios, and low stellar metallicities. Their bursty star formation histories are characterized by peak specific star formation rates in excess of 50-100 1/Gyr, far outside the realm of normal, more massive galaxies, and in agreement with observations of extreme emission-line starbursting dwarfs by the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey. Metal-enriched galactic outflows produce sub-solar effective yields and pollute with heavy elements a Mpc-size region of the intergalactic medium, but are not sufficient to completely quench star formation activity and are not ubiquitous in our dwarfs. Within the limited size of the sample, our simulations appear to simultaneously reproduce the observed stellar mass and cold gas content, resolved star formation histories, stellar kinematics, and metallicities of field dwarfs in the Local Volume.

preprint2014arXiv

The Dark Matter Contribution to Galactic Diffuse Gamma Ray Emission

Observations of diffuse Galactic gamma ray emission (DGE) by the Fermi Large Area Telescope (LAT) allow a detailed study of cosmic rays and the interstellar medium. However, diffuse emission models of the inner Galaxy underpredict the Fermi-LAT data at energies above a few GeV and hint at possible non-astrophysical sources including dark matter (DM) annihilations or decays. We present a study of the possible emission components from DM using the high-resolution Via Lactea II N-body simulation of a Milky Way-sized DM halo. We generate full-sky maps of DM annihilation and decay signals that include modeling of the adiabatic contraction of the host density profile, Sommerfeld enhanced DM annihilations, $p$-wave annihilations, and decaying DM. We compare our results with the DGE models produced by the Fermi-LAT team over different sky regions, including the Galactic center, high Galactic latitudes, and the Galactic anti-center. This work provides possible templates to fit the observational data that includes the contribution of the subhalo population to DM gamma-ray emission, with the significance depending on the annihilation/decay channels and the Galactic regions being considered.

preprint2014arXiv

The Distribution of Dark Matter in the Milky Way's Disk

We present an analysis of the effects of dissipational baryonic physics on the local dark matter (DM) distribution at the location of the Sun, with an emphasis on the consequences for direct detection experiments. Our work is based on a comparative analysis of two cosmological simulations of a Milky Way halo: Eris, a full hydrodynamic simulation, and ErisDark, its DM-only realization. We find that two distinct processes lead in Eris to a 30% enhancement of DM in the disk plane at the location of the Sun: the accretion and disruption of satellites resulting in a DM component with net angular momentum and the contraction of baryons pulling DM into the disk plane without forcing it to co-rotate. Owing to its particularly quiescent merger history for dark halos of Milky Way mass, the co-rotating dark disk in Eris is less massive than what has been suggested by previous work, contributing only 9% of the local DM density. The speed distribution in Eris is broadened and shifted to higher speeds compared to its DM-only twin simulation ErisDark. At high speeds $f(v)$ falls more steeply in Eris than in ErisDark or the SHM, easing the tension between recent results from the CDMS-II and XENON100 experiments. The non-Maxwellian aspects of $f(v)$ are still present, but much less pronounced in Eris than in DM-only runs. The weak dark disk increases the time-averaged scattering rate by only a few percent at low recoil energies. On the high velocity tail, however, the increase in typical speeds due to baryonic contraction results in strongly enhanced mean scattering rates compared to ErisDark, although they are still suppressed compared to the SHM. Similar trends are seen regarding the amplitude of the annual modulation, while the modulated fraction is increased compared to the SHM and decreased compared to ErisDark.

preprint2014arXiv

The Photon Underproduction Crisis

We examine the statistics of the low-redshift Lyman-alpha forest from smoothed particle hydrodynamic simulations in light of recent improvements in the estimated evolution of the cosmic ultraviolet background (UVB) and recent observations from the Cosmic Origins Spectrograph (COS). We find that the value of the metagalactic photoionization rate required by our simulations to match the observed properties of the low-redshift Lyman-alpha forest is a factor of 5 larger than the value predicted by state-of-the art models for the evolution of this quantity. This mismatch results in the mean flux decrement of the Lyman-alpha forest being underpredicted by at least a factor of 2 (a 10-sigma discrepancy with observations) and a column density distribution of Lyman-alpha forest absorbers systematically and significantly elevated compared to observations over nearly two decades in column density. We examine potential resolutions to this mismatch and find that either conventional sources of ionizing photons (galaxies and quasars) must be significantly elevated relative to current observational estimates or our theoretical understanding of the low-redshift universe is in need of substantial revision.

preprint2014arXiv

Towards a Unified Description of the Intergalactic Medium at Redshift z~2.5

We examine recent measurements of the z~2.5 intergalactic medium (IGM) which constrain the HI frequency distribution f(N) and the mean free path MFP to ionizing radiation. We argue that line-blending and the clustering of strong absorption-line systems have led previous authors to systematically overestimate the effective Lyman limit opacity, yielding too small of a MFP for the IGM. We further show that recently published measurements of f(N) at log NHI ~ 16 lie in strong disagreement, implying underestimated uncertainty from sample variance and/or systematics like line-saturation. Allowing for a larger uncertainty in the f(N) measurements, we provide a new cubic Hermite spline model for f(N) which reasonably satisfies all of the observational constraints under the assumption of randomly distributed absorption systems. We caution, however, that this formalism is invalid in light of absorber clustering and use a toy model to estimate the effects. Future work must properly account for the non-Poissonian nature of the IGM.

preprint2013arXiv

A "Light," Centrally-Concentrated Milky Way Halo?

We discuss a novel approach to "weighing" the Milky Way dark matter halo, one that combines the latest samples of halo stars selected from the Sloan Digital Sky Survey (SDSS) with state-of-the-art numerical simulations of Milky Way analogs. The fully cosmological runs employed in the present study include "Eris", one of the highest-resolution hydrodynamical simulations of the formation of a M_vir=8e11 M_sun late-type spiral, and the dark-matter only M_vir=1.7e12 M_sun "Via Lactea II" simulation. Eris provides an excellent laboratory for creating mock SDSS samples of tracer halo stars, and we successfully compare their density, velocity anisotropy, and radial velocity dispersion profiles with the observational data. Most mock SDSS realizations show the same "cold veil" recently observed in the distant stellar halo of the Milky Way, with tracers as cold as sigma_los ~ 50 km/s between 100-150 kpc. Controlled experiments based on the integration of the spherical Jeans equation as well as a particle tagging technique applied to Via Lactea II show that a "heavy" M_vir 2e12 M_sun realistic host produces a poor fit to the kinematic SDSS data. We argue that these results offer added evidence for a "light," centrally-concentrated Milky Way halo.

preprint2013arXiv

A Population of Relic Intermediate-Mass Black Holes in the Halo of the Milky Way

If "seed" central black holes were common in the subgalactic building blocks that merged to form present-day massive galaxies, then relic intermediate-mass black holes (IMBHs) should be present in the Galactic bulge and halo. We use a particle tagging technique to dynamically populate the N-body Via Lactea II high-resolution simulation with black holes, and assess the size, properties, and detectability of the leftover population. The method assigns a black hole to the most tightly bound central particle of each subhalo at infall according to an extrapolation of the M_BH-sigma_* relation, and self-consistently follows the accretion and disruption of Milky Way progenitor dwarfs and their holes in a cosmological "live" host from high redshift to today. We show that, depending on the minimum stellar velocity dispersion, sigma_m, below which central black holes are assumed to be increasingly rare, as many as ~2000 (sigma_m=3 km/s) or as few as ~70 (sigma_m=12 km/s) IMBHs may be left wandering in the halo of the Milky Way today. The fraction of IMBHs kicked out of their host by gravitational recoil is < 20%. We identify two main Galactic subpopulations, "naked" IMBHs, whose host subhalos were totally destroyed after infall, and "clothed" IMBHs residing in dark matter satellites that survived tidal stripping. Naked IMBHs typically constitute 40-50% of the total and are more centrally concentrated. We show that, in the sigma_m=12 km/s scenario, the clusters of tightly bound stars that should accompany naked IMBHs would be fainter than m_V=16 mag, spatially resolvable, and have proper motions of 0.1-10 milliarcsec per year. Their detection may provide an observational tool to constrain the formation history of massive black holes in the early Universe.

preprint2013arXiv

An Off-center Density Peak in the Milky Way's Dark Matter Halo?

We show that the position of the central dark matter density peak may be expected to differ from the dynamical center of the Galaxy by several hundred parsec. In Eris, a high resolution cosmological hydrodynamics simulation of a realistic Milky-Way-analog disk galaxy, this offset is 300 - 400 pc (~3 gravitational softening lengths) after z=1. In its dissipationless dark-matter-only twin simulation ErisDark, as well as in the Via Lactea II and GHalo simulations, the offset remains below one softening length for most of its evolution. The growth of the DM offset coincides with a flattening of the central DM density profile in Eris inwards of ~1 kpc, and the direction from the dynamical center to the point of maximum DM density is correlated with the orientation of the stellar bar, suggesting a bar-halo interaction as a possible explanation. A dark matter density offset of several hundred parsec greatly affects expectations of the dark matter annihilation signals from the Galactic Center. It may also support a dark matter annihilation interpretation of recent reports by Weniger (2012) and Su & Finkbeiner (2012) of highly significant 130 GeV gamma-ray line emission from a region 1.5 degrees (~200 parsec projected) away from Sgr A* in the Galactic plane.

preprint2013arXiv

Direct Gravitational Imaging of Intermediate Mass Black Holes in Extragalactic Halos

A galaxy halo may contain a large number of intermediate mass black holes (IMBHs) with masses in the range of 10^{2-6} solar mass. We propose to directly detect these IMBHs by observing multiply imaged QSO-galaxy or galaxy-galaxy strong lens systems in the submillimeter bands with high angular resolution. The silhouette of an IMBH in the lensing galaxy halo would appear as either a monopole-like or a dipole-like variation at the scale of the Einstein radius against the Einstein ring of the dust-emitting region surrounding the QSO. We use a particle tagging technique to dynamically populate a Milky Way-sized dark matter halo with black holes, and show that the surface mass density and number density of IMBHs have power-law dependences on the distance from the center of the host halo if smoothed on a scale of ~ 1 kpc. Most of the black holes orbiting close to the center are freely roaming as they have lost their dark matter hosts during infall due to tidal stripping. Next generation submillimeter telescopes with high angular resolution (< 0.3 mas) will be capable of directly mapping such off-nuclear freely roaming IMBHs with a mass of ~ 10^6 solar mass in a lensing galaxy that harbours a O(10^9) solar mass supermassive black hole in its nucleus.

preprint2013arXiv

Dwarf Galaxy Formation with H2-Regulated Star Formation: II. Gas-Rich Dark Galaxies at Redshift 2.5

We present a cosmological hydrodynamic simulation of the formation of dwarf galaxies at redshifts z>~2.5 using a physically-motivated model for H2-regulated star formation. Our simulation, performed using the Enzo code and reaching a peak resolution of 109 proper parsecs at z=2.5, extends the results of Kuhlen et al. (2012) to significantly lower redshifts. We show that a star formation prescription regulated by the local H2 abundance leads to the suppression of star formation in dwarf galaxy halos with M_h <~ 10^10 Msun and to a large population of gas-rich "dark galaxies" at z=2.5 with low star formation efficiencies and gas depletion timescales >20 Gyr. The fraction of dark galaxies is 60% at M_h ~ 10^10 Msun and increases rapidly with decreasing halo mass. Dark galaxies form late and their gaseous disks never reach the surface densities, > ~5700 Msun / pc^2 (Z/10^-3 Zsun)^(-0.88), that are required to build a substantial molecular fraction. Despite this large population of dark galaxies, we show that our H2-regulated simulation is consistent with both the observed luminosity function of galaxies and the cosmological mass density of neutral gas at z>~2.5. Moreover, our results provide a theoretical explanation for the recent detection in fluorescent Ly-alpha emission of gaseous systems at high redshift with little or no associated star formation. We further propose that H2-regulation may offer a fresh solution to a number of outstanding "dwarf galaxy problems" in LambdaCDM. In particular, H2-regulation leads galaxy formation to become effectively stochastic on mass scales of M_h ~ 10^10 Msun, and thus these massive dwarfs are not "too big to fail".

preprint2013arXiv

Inside Out and Upside Down: Tracing the Assembly of a Simulated Disk Galaxy Using Mono-Age Stellar Populations

We analyze the present-day structure and assembly history of a high resolution hydrodynamic simulation of the formation of a Milky Way (MW)-like disk galaxy, from the "Eris" simulation suite, dissecting it into cohorts of stars formed at different epochs of cosmic history. At z=0, stars with t_form < 2 Gyr mainly occupy the stellar spheroid, with the oldest (earliest forming) stars having more centrally concentrated profiles. The younger age cohorts populate disks of progressively longer radial scale length and shorter vertical scale height. At a given radius, the vertical density profiles and velocity dispersions of stars vary smoothly as a function of age, and the superposition of old, vertically-extended and young, vertically-compact cohorts gives rise to a double-exponential profile like that observed in the MW. Turning to formation history, we find that the trends of spatial structure and kinematics with stellar age are largely imprinted at birth, or immediately thereafter. Stars that form during the active merger phase at z>3 are quickly scattered into rounded, kinematically hot configurations. The oldest disk cohorts form in structures that are radially compact and relatively thick, while subsequent cohorts form in progressively larger, thinner, colder configurations from gas with increasing levels of rotational support. The disk thus forms "inside-out" in a radial sense and "upside-down" in a vertical sense. Secular heating and radial migration influence the final state of each age cohort, but the changes they produce are small compared to the trends established at formation. The predicted correlations of stellar age with spatial and kinematic structure are in good qualitative agreement with the correlations observed for mono-abundance stellar populations in the MW.

preprint2013arXiv

The circumgalactic medium of massive galaxies at z~3: a test for stellar feedback, galactic outflows, and cold streams

We present new results on the kinematics, thermal and ionization state, and spatial distribution of metal-enriched gas in the circumgalactic medium (CGM) of massive galaxies at redshift 3, using the "Eris" suite of cosmological "zoom-in" simulations. The reference run adopts a blastwave scheme for supernova feedback that produces galactic outflows, a star formation recipe based on a high gas density threshold, metal-dependent radiative cooling, and a model for the diffusion of metals and thermal energy. Synthetic spectra through the multiphase CGM produce interstellar absorption line strengths of Lya, CII, CIV, SiII, and SiIV as a function of galactocentric impact parameter (scaled to the virial radius) that are in broad agreement with those observed at high-redshift by Steidel et al. (2010). Only about one third of all the gas within R_vir is outflowing. The fraction of sightlines within one virial radius that intercept optically thick material is 27%, in agreement with recent observations by Rudie et al. (2012). Such optically thick absorption is shown to trace inflowing "cold" streams that penetrate deep inside the virial radius. The streams, enriched to metallicities above 0.01 solar, give origin to strong (log N > 13) CII absorption with a covering factor of 22% within R_vir and 10% within 2 R_vir. Galactic outflows do not cause any substantial suppression of the cold accretion mode. The central galaxy is surrounded by a large OVI halo, with a typical column density log N>14 and a near unity covering factor maintained all the way out to 150 kpc. This matches the trends recently observed in star-forming galaxies at low redshift by Tumlinson et al. (2011). Our zoom-in simulations of this single system appear to reproduce quantitatively the complex baryonic processes that determine the exchange of matter, energy, and metals between galaxies and their surroundings. (Abridged)

preprint2013arXiv

The HST/ACS+WFC3 Survey for Lyman Limit Systems II: Science

We present the first science results from our Hubble Space Telescope Survey for Lyman limit absorption systems (LLS) using the low dispersion spectroscopic modes of the Advanced Camera for Surveys and the Wide Field Camera 3. Through an analysis of 71 quasars, we determine the incidence frequency of LLS per unit redshift and per unit path length, l(z) and l(x) respectively, over the redshift range 1 < z< 2.6, and find a weighted mean of l(x)=0.29 +/-0.05 for 2.0 < z < 2.5 through a joint analysis of our sample and that of Ribaudo et al. (2011). Through stacked spectrum analysis, we determine a median (mean) value of the mean free path to ionizing radiation at z=2.4 of lambda_mfp = 243(252)h^(-1) Mpc, with an error on the mean value of +/- 43h^(-1) Mpc. We also re-evaluate the estimates of lambda_mfp from Prochaska et al. (2009) and place constraints on the evolution of lambda_mfp with redshift, including an estimate of the "breakthrough" redshift of z = 1.6. Consistent with results at higher z, we find that a significant fraction of the opacity for absorption of ionizing photons comes from systems with N_HI <= 10^{17.5} cm^(-2) with a value for the total Lyman opacity of tau_lyman = 0.40 +/- 0.15. Finally, we determine that at minimum, a 5-parameter (4 power-law) model is needed to describe the column density distribution function f(N_HI, X) at z \sim 2.4, find that f(N_HI,X) undergoes no significant change in shape between z \sim 2.4 and z \sim 3.7, and provide our best fit model for f(N_HI,X).

preprint2012arXiv

Insight Into the Formation of the Milky Way Through Cold Halo Substructure. III. Statistical Chemical Tagging in the Smooth Halo

We find that the relative contribution of satellite galaxies accreted at high redshift to the stellar population of the Milky Way's smooth halo increases with distance, becoming observable relative to the classical smooth halo about 15 kpc from the Galactic center. In particular, we determine line-of-sight-averaged [Fe/H] and [alpha/Fe] in the metal-poor main-sequence turnoff (MPMSTO) population along every Sloan Extension for Galactic Understanding and Exploration (SEGUE) spectroscopic line of sight. Restricting our sample to those lines of sight along which we do not detect elements of cold halo substructure (ECHOS), we compile the largest spectroscopic sample of stars in the smooth component of the halo ever observed in situ beyond 10 kpc. We find significant spatial autocorrelation in [Fe/H] in the MPMSTO population in the distant half of our sample beyond about 15 kpc from the Galactic center. Inside of 15 kpc however, we find no significant spatial autocorrelation in [Fe/H]. At the same time, we perform SEGUE-like observations of N-body simulations of Milky Way analog formation. While we find that halos formed entirely by accreted satellite galaxies provide a poor match to our observations of the halo within 15 kpc of the Galactic center, we do observe spatial autocorrelation in [Fe/H] in the simulations at larger distances. This observation is an example of statistical chemical tagging and indicates that spatial autocorrelation in metallicity is a generic feature of stellar halos formed from accreted satellite galaxies.

preprint2012arXiv

Pseudobulge Formation as a Dynamical Rather than a Secular Process

We investigate the formation and evolution of the pseudobulge in "Eris", a high-resolution N-body + smoothed particle hydrodynamics (SPH) cosmological simulation that successfully reproduces a Milky Way-like massive late-type spiral in a cold dark matter (LCDM) universe. At the present epoch, Eris has a virial mass Mvir=8x10^11 Msun, a photometric stellar mass M*=3.2x10^10 Msun, a bulge-to-total ratio B/T = 0.26, and a weak nuclear bar. We find that the bulk of the pseudobulge forms quickly at high redshift via a combination of non-axisymmetric disk instabilities and tidal interactions or mergers both occurring on dynamical timescales, not through slow secular processes at lower redshift. Its subsequent evolution is not strictly secular either, and is closely intertwined with the evolution of the stellar bar. In fact, the structure that we recognize as a pseudobulge today evolves from a stellar bar that formed at high redshift, was destroyed by minor mergers at z~3, reformed shortly after, and weakened again following a steady gas inflow at z<1. The gradual dissolution of the bar ensues at z=1 and continues until the present without increasing the stellar velocity dispersion in the inner regions. In this scenario the pseudobulge is not a separate component from the inner disk in terms of formation path, rather it is the first step in the inside-out formation of the baryonic disk, in agreement with the fact that pseudobulges of massive spiral galaxies have typically a dominant old stellar population. If our simulations do indeed reproduce the formation mechanisms of massive spirals, then the progenitors of late-type galaxies should have strong bars and small photometric pseudobulges at high redshift.

preprint2012arXiv

The origin of metals in the circum-galactic medium of massive galaxies at z=3

We present a detailed study of the metal-enriched circum-galactic medium of a massive galaxy at z=3 using "ErisMC", a new cosmological hydrodynamic "zoom-in" simulations of a disk galaxy with mass comparable to the Milky Way. The run adopts a blastwave scheme for supernova feedback that generates galactic outflows without explicit wind particles, a star formation recipe based on a high gas density threshold, and high temperature metal cooling. ErisMC's main progenitor at z=3 resembles a "Lyman break" galaxy of mass M_vir=2.4e11 M_sun, virial radius R_vir=48 kpc, and star formation rate 18 M_sun/yr, and its metal-enriched CGM extends to 200 (physical) kpc from its center. Approximately 41, 9, and 50 percent of all gas-phase metals at z=3 are locked in a hot (T> 3e5 K), warm (3e5 >T> 3e4 K), and cold (T< 3e4 K) medium, respectively. We identify three sources of heavy elements: 1) the main host, responsible for 60% of all the metals found within 3R_vir; 2) its satellite progenitors, responsible for 28% of all the metals within 3R_vir, and for only 5% of those beyond 3R_vir; and nearby dwarfs, which give origin to 12% of all the metals within 3R_vir and 95% of those beyond 3R_vir. Late (z<5) galactic "superwinds" account for only 9% of all the metals observed beyond 2R_vir, the bulk having been released at redshifts 5< z < 8 by early star formation and outflows. In the CGM, lower overdensities are typically enriched by `older', colder metals. Heavy elements are accreted onto Eris along filaments via low-metallicity cold inflows, and are ejected hot via galactic outflows at a few hundred km/s. The outflow mass-loading factor is of order unity for the main halo, but can exceed 10 for nearby dwarfs. We stress that our "zoom-in" simulation focuses on the CGM of a single massive system and cannot describe the enrichment history of the IGM as a whole. (abridged)

preprint2011arXiv

A galaxy as the source of a Civ absorption system close to the epoch of reionization

We find a bright (L_{UV}=2.5 L*_{z=6}) Lyman alpha emitter at redshift z=5.719 at a projected distance of 79 physical kpc from a strong triply ionized carbon (Civ) absorption system at redshift z=5.7238 previously reported in the spectrum of the z_{em} = 6.309 QSO SDSS J1030+0524. This is the highest redshift galaxy-absorber pair detected to-date, supporting the idea that galaxy-wide outflows were already in place at the end of the epoch of reionization. The proximity of this object makes it the most likely source of metals, consistent with models of outflows at lower redshift where significant observational evidence relates metal absorption systems with galaxies hosting outflows. In a typical outflow scenario, a wind of 200 km/s, active since the universe was only 0.6 Gyr old (z ~8.4), could eject metals out to 79 kpc at z=5.719. Although the origin of metals in the intergalactic medium (IGM) is still under debate, our results are consistent with predictions from cosmological simulations which reproduce the evolution of the cosmic density of Civ, from z ~ 6 to the present day based on outflow-driven enrichment of the IGM. We also report two more Lyman alpha emitters in this field, at z=5.973\pm 0.002 and z=5.676\pm 0.002 respectively, the former confirming the original identification by Stiavelli et al. Our results suggest that the colour cut typically used to identify i-dropouts (i_{775}-z_{850}>1.3) misses a non-negligible fraction of blue galaxies with faint UV continuum at z \geq 5.7.

preprint2011arXiv

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey

The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) is designed to document the first third of galactic evolution, over the approximate redshift (z) range 8--1.5. It will image >250,000 distant galaxies using three separate cameras on the Hubble Space Telescope, from the mid-ultraviolet to the near-infrared, and will find and measure Type Ia supernovae at z>1.5 to test their accuracy as standardizable candles for cosmology. Five premier multi-wavelength sky regions are selected, each with extensive ancillary data. The use of five widely separated fields mitigates cosmic variance and yields statistically robust and complete samples of galaxies down to a stellar mass of 10^9 M_\odot to z \approx 2, reaching the knee of the ultraviolet luminosity function (UVLF) of galaxies to z \approx 8. The survey covers approximately 800 arcmin^2 and is divided into two parts. The CANDELS/Deep survey (5σ point-source limit H=27.7 mag) covers \sim 125 arcmin^2 within GOODS-N and GOODS-S. The CANDELS/Wide survey includes GOODS and three additional fields (EGS, COSMOS, and UDS) and covers the full area to a 5σ point-source limit of H \gtrsim 27.0 mag. Together with the Hubble Ultra Deep Fields, the strategy creates a three-tiered "wedding cake" approach that has proven efficient for extragalactic surveys. Data from the survey are nonproprietary and are useful for a wide variety of science investigations. In this paper, we describe the basic motivations for the survey, the CANDELS team science goals and the resulting observational requirements, the field selection and geometry, and the observing design. The Hubble data processing and products are described in a companion paper.

preprint2011arXiv

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey - The Hubble Space Telescope Observations, Imaging Data Products and Mosaics

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

preprint2011arXiv

Delayed Enrichment by Unseen Galaxies: Explaining the Rapid Rise in IGM CIV Absorption from z = 6-5

In the near future, measurements of metal absorption features in the intergalactic medium (IGM) will become an important constraint on models of the formation and evolution of the earliest galaxies, the properties of the first stars, and the reionization and enrichment of the IGM. The first measurement of a metal abundance in the IGM at a redshift approaching the epoch of reionization already offers intriguing hints. Between z=5.8 and 4.7 (a 0.3 Gyr interval only 1 Gyr after the big bang), the measured density of CIV absorbers in the IGM increased by a factor of ~ 3.5 (Ryan-Weber et al. 2009; Becker, Rauch & Sargent 2009). If these values prove to be accurate, they pose two puzzles: (1) The total amount of CIV at z=5.8 implies too little star formation to reionize the IGM by z=6 or to match the WMAP electron scattering optical depth (tau). (2) The rapid growth from z = 6-5 is faster than the buildup of stellar mass or the increase in the star formation rate density over the same interval. We show that a delay of ~ 0.4-0.7 Gyr between the instantaneous production of ionizing photons and the later production of metal absorption features (added to the delay due to stellar lifetimes) can provide the full explanation for both puzzles. We calculate the delay in metal production due to finite stellar lifetimes alone and find that it is too short to explain the rapid CIV density increase. The additional delay could naturally be explained as the result of ~ 200 km/s outflows carrying carbon to distances of ~ 100 kpc, the typical distance between galaxies and CIV absorbers in enrichment simulations, and the typical outflow or absorption region scale observed at z ~ 2-3.

preprint2011arXiv

Forming Realistic Late-Type Spirals in a LCDM Universe: The Eris Simulation

Simulations of the formation of late-type spiral galaxies in a cold dark matter LCDM universe have traditionally failed to yield realistic candidates. Here we report a new cosmological N-body/SPH simulation of extreme dynamic range in which a close analog of a Milky Way disk galaxy arises naturally. Termed Eris, the simulation follows the assembly of a galaxy halo of mass Mvir=7.9x10^11 Msun with a total of N=18.6 million particles (gas + dark matter + stars) within the final virial radius, and a force resolution of 120 pc. It includes radiative cooling, heating from a cosmic UV field and supernova explosions, a star formation recipe based on a high gas density threshold (nSF=5 atoms cm^-3 rather than the canonical nSF=0.1 atoms cm^-3), and neglects AGN feedback. At the present epoch, the simulated galaxy has an extended rotationally-supported disk with a radial scale length Rd=2.5 kpc, a gently falling rotation curve with circular velocity at 2.2 disk scale lenghts of V2.2=214 km/s, a bulge-to-disk ratio B/D=0.35, and a baryonic mass fraction that is 30% below the cosmic value. The disk is thin, is forming stars in the region of the Sigma_SFR - Sigma_HI plane occupied by spiral galaxies, and falls on the photometric Tully-Fisher and the stellar mass-halo virial mass relations. Hot (T>3x10^5 K), X-ray luminous halo gas makes only 26% of the universal baryon fraction and follows a flattened density profile proportional to r^-1.13 out to r=100 kpc. Eris appears then to be the first cosmological hydrodynamic simulation in which the galaxy structural properties, the mass budget in the various components, and the scaling relations between mass and luminosity are all consistent with a host of observational constraints. (Abridged)

preprint2011arXiv

New Population Synthesis Models of the Cosmic UV Background

We present improved synthesis models of the evolving spectrum of the UV/X-ray diffuse background, updating and extending our previous results. A "minimal cosmic reionization model" is presented in which the galaxy UV emissivity traces recent determinations of the cosmic history of star formation, the luminosity-weighted escape fraction of hydrogen-ionizing radiation increases rapidly with redshift, from 1% at z=3 to 50% at z=8.5, the clumping factor of the intergalactic medium is C=3, and Population III stars make a negligible contribution to the metagalactic flux. The model provides a good fit to the hydrogen-ionization rates inferred from flux decrement measurements, and yields an optical depth to Thomson scattering, tau_e=0.085, that is in agreement with WMAP results.

preprint2011arXiv

On the assembly of the Milky Way dwarf satellites and their common mass scale

We use a particle tagging technique to dynamically populate the N-body Via Lactea II high-resolution simulation with stars. The method is calibrated using the observed luminosity function of Milky Way satellites and the concentration of their stellar populations, and self-consistently follows the accretion and disruption of progenitor dwarfs and the build-up of the stellar halo in a cosmological "live host". Simple prescriptions for assigning stellar populations to collisionless particles are able to reproduce many properties of the observed Milky Way halo and its surviving dwarf satellites, like velocity dispersions, sizes, brightness profiles, metallicities, and spatial distribution. Our model predicts the existence of approximately 1,850 subhalos harboring "extremely faint" satellites (with mass-to-light ratios >5,000) lying beyond the Sloan Digital Sky Survey detection threshold. Of these, about 20 are "first galaxies", i.e. satellites that formed a stellar mass above 10 Msun before redshift 9. The ten most luminous satellites (L> 1e6 Lsun) in the simulation are hosted by subhalos with peak circular velocities today in the range V_max=10-40 km/s that have shed between 80% and 99% of their dark mass after being accreted at redshifts 1.7< z <4.6. The satellite maximum circular velocity and stellar line-of-sight velocity dispersion today follow the relation V_max=2.2 sigma_los. We apply a standard mass estimation algorithm based on Jeans modelling of the line-of-sight velocity dispersion profiles to the simulated dwarf spheroidals, and test the accuracy of this technique. The inner (within 300 pc) mass-luminosity relation for currently detectable satellites is nearly flat in our model, in qualitative agreement with the "common mass scale" found in Milky Way dwarfs. We do, however, predict a weak, but significant positive correlation for these objects: M_300 ~L^{0.088 \pm 0.024}.

preprint2011arXiv

Radiative transfer in a clumpy universe: IV. New synthesis models of the cosmic UV/X-ray background

We present improved synthesis models of the evolving spectrum of the UV/X-ray diffuse background, updating and extending our previous results. Five new main components are added to our radiative transfer code CUBA: (1) the sawtooth modulation of the background intensity from resonant line absorption in the Lyman series of cosmic hydrogen and helium; (2) the X-ray emission from obscured and unobscured quasars; (3) a piecewise parameterization of the distribution in redshift and column density of intergalactic absorbers that fits recent measurements of the mean free path of 1 ryd photons; (4) an accurate treatment of the photoionization structure of absorbers; and (5) the UV emission from star-forming galaxies at all redshifts. We provide tables of the predicted HI and HeII photoionization and photoheating rates for use, e.g., in cosmological hydrodynamics simulations of the Lya forest, and a new metallicity-dependent calibration to the UV luminosity density-star formation rate density relation. A "minimal cosmic reionization model" is also presented in which the galaxy UV emissivity traces recent determinations of the cosmic history of star formation, the luminosity-weighted escape fraction of hydrogen-ionizing radiation increases rapidly with lookback time, the clumping factor of the high-redshift intergalactic medium evolves following the results of recent hydrodynamic simulations, and Population III stars and miniquasars make a negligible contribution to the metagalactic flux. The model provides a good fit to the hydrogen-ionization rates inferred from flux decrement and proximity effect measurements, predicts that cosmological HII (HeIII) regions overlap at redshift 6.7 (2.8), and yields an optical depth to Thomson scattering that is in agreement with WMAP results. (Abridged)

preprint2010arXiv

Fermi-LAT Sensitivity to Dark Matter Annihilation in Via Lactea II Substructure

We present a study of the ability of the Fermi Gamma-ray Space Telescope to detect dark-matter annihilation signals from the Galactic subhalos predicted by the Via Lactea II N-body simulation. We implement an improved formalism for estimating the boost factor needed to account for the effect of dark-matter clumping on scales below the resolution of the simulation, and we incorporate a detailed Monte Carlo simulation of the response of the Fermi-LAT telescope, including a simulation of its all-sky observing mode integrated over a ten year mission. We find that for WIMP masses up to about 150 GeV in standard supersymmetric models with velocity-averaged cross section 3*10^-26 cm^3 s^-1, a few subhalos could be detectable with >5 standard deviations significance and would likely deviate significantly from the appearance of a point source.

preprint2010arXiv

Recoiling Massive Black Holes in Gas-Rich Galaxy Mergers

The asymmetric emission of gravitational waves produced during the coalescence of a massive black hole (MBH) binary imparts a velocity "kick" to the system that can displace the hole from the center of its host. Here we study the trajectories and observability of MBHs recoiling in three (one major, two minor) gas-rich galaxy merger remnants that were previously simulated at high resolution, and in which the pairing of the MBHs had been shown to be successful. We run new simulations of MBHs recoiling in the major merger remnant with Mach numbers in the range 1<M<6, and use simulation data to construct a semi-analytical model for the orbital evolution of MBHs in gas-rich systems. We show that: 1) in major merger remnants the energy deposited by the moving hole into the rotationally supported, turbulent medium makes a negligible contribution to the thermodynamics of the gas. This contribution becomes significant in minor merger remnants, potentially allowing for an electromagnetic signature of MBH recoil; 2) in major merger remnants, the combination of both deeper central potential well and drag from high-density gas confines even MBHs with kick velocities as high as 1200 km/s within 1 kpc from the host's center; 3) kinematically offset nuclei may be observable for timescales of a few Myr in major merger remnants in the case of recoil velocities in the range 700-1,000 km/s; 4) in minor mergers remnants the effect of gas drag is weaker, and MBHs with recoil speeds in the range 300-600 km/s will wander through the host halo for longer timescales. When accounting for the probability distribution of kick velocities, however, we find that the likelihood of observing recoiling MBHs in gas-rich galaxy mergers is very low, typically below 10^-5 - 10^-6.

preprint2010arXiv

Tidal stellar disruptions by massive black hole pairs: II. Decaying binaries

Tidal stellar disruptions have traditionally been discussed as a probe of the single, massive black holes (MBHs) that are dormant in the nuclei of galaxies. In Chen et al. (2009), we used numerical scattering experiments to show that three-body interactions between bound stars in a stellar cusp and a non-evolving "hard" MBH binary will also produce a burst of tidal disruptions, caused by a combination of the secular "Kozai effect" and by close resonant encounters with the secondary hole. Here we derive basic analytical scalings of the stellar disruption rates with the system parameters, assess the relative importance of the Kozai and resonant encounter mechanisms as a function of time, discuss the impact of general relativistic (GR) and extended stellar cusp effects, and develop a hybrid model to self-consistently follow the shrinking of an MBH binary in a stellar background, including slingshot ejections and tidal disruptions. In the case of a fiducial binary with primary hole mass M_1=10^7\msun and mass ratio q=M_2/M_1=1/81, embedded in an isothermal cusp, we derive a stellar disruption rate \dot{N_*} ~ 0.2/yr lasting ~ 3X10^5 yr. This rate is 3 orders of magnitude larger than the corresponding value for a single MBH fed by two-body relaxation, confirming our previous findings. For q<<0.01, the Kozai/chaotic effect could be quenched due to GR/cusp effects by an order of magnitude, but even in this case the stellar-disruption rate is still two orders of magnitude larger than that given by standard relaxation processes around a single MBH. Our results suggest that >~10% of the tidal-disruption events may originate in MBH binaries.

preprint2009arXiv

At the heart of the matter: the origin of bulgeless dwarf galaxies and Dark Matter cores

For almost two decades the properties of "dwarf" galaxies have challenged the Cold Dark Matter (CDM) paradigm of galaxy formation. Most observed dwarf galaxies consists of a rotating stellar disc embedded in a massive DM halo with a near constant-density core. Yet, models based on the CDM scenario invariably form galaxies with dense spheroidal stellar "bulges" and steep central DM profiles, as low angular momentum baryons and DM sink to the center of galaxies through accretion and repeated mergers. Processes that decrease the central density of CDM halos have been identified, but have not yet reconciled theory with observations of present day dwarfs. This failure is potentially catastrophic for the CDM model, possibly requiring a different DM particle candidate. This Letter presents new hydrodynamical simulations in a Lambda$CDM framework where analogues of dwarf galaxies, bulgeless and with a shallow central DM profile, are formed. This is achieved by resolving the inhomogeneous interstellar medium, resulting in strong outflows from supernovae explosions which remove low angular momentum gas. This inhibits the formation of bulges and decreases the dark-matter density to less than half within the central kiloparsec. Realistic dwarf galaxies are thus shown to be a natural outcome of galaxy formation in the CDM scenario.

preprint2009arXiv

GeV Gamma-Ray Attenuation and the High-Redshift UV Background

We present new calculations of the evolving UV background out to the epoch of cosmological reionization and make predictions for the amount of GeV gamma-ray attenuation by electron-positron pair production. Our results are based on recent semi-analytic models of galaxy formation, which provide predictions of the dust-extinguished UV radiation field due to starlight, and empirical estimates of the contribution due to quasars. We account for the reprocessing of ionizing photons by the intergalactic medium. We test whether our models can reproduce estimates of the ionizing background at high redshift from flux decrement analysis and proximity effect measurements from quasar spectra, and identify a range of models that can satisfy these constraints. Pair-production against soft diffuse photons leads to a spectral cutoff feature for gamma rays observed between 10 and 100 GeV. This cutoff varies with redshift and the assumed star formation and quasar evolution models. We find only negligible amounts of absorption for gamma rays observed below 10 GeV for any emission redshift. With observations of high-redshift sources in sufficient numbers by the Fermi Gamma-ray Space Telescope and new ground-based instruments it should be possible to constrain the extragalactic background light in the UV and optical portion of the spectrum.

preprint2009arXiv

He II absorption and the sawtooth spectrum of the cosmic far-UV background

Cosmic ultraviolet background radiation between 3 and 4 Ryd is reprocessed by resonant line absorption in the Lyman series of intergalactic He II. This process results in a sawtooth modulation of the radiation spectrum from the He II Lya frequency to the Lyman limit. The size of this modulation is a sensitive probe of the epoch of helium reionization and of the sources that keep the intergalactic medium (IGM) highly ionized. For large absorption opacities, the background intensity will peak at frequencies just above each resonance, go to zero at resonance, and fluctuate greatly just below resonance. The He II sawtooth modulation may be one of the missing ingredients needed in the modelling of the abundances of metal ions like C III and Si IV observed in the IGM at redshift 3.

preprint2009arXiv

Probing the Epoch of Reionization with Milky Way Satellites

While the connection between high-redshift star formation and the local universe has recently been used to understand the observed population of faint dwarf galaxies in the Milky Way (MW) halo, we explore how well these nearby objects can probe the epoch of first light. We construct a detailed, physically motivated model for the MW satellites based on the state-of-the-art Via Lactea II dark-matter simulations. Our model incorporates molecular hydrogen (H_2) cooling in low-mass systems and inhomogeneous photo-heating feedback during the internal reionization of our own galaxy. We find that the existence of MW satellites fainter than M_V ~ -5 is strong evidence for H_2 cooling in low-mass halos, while satellites with -5 > M_V > -9 were affected by hydrogen cooling and photoheating feedback. The age of stars in very low-luminosity systems and the minimum luminosity of these satellites are key predictions of our model. Most of the stars populating the brightest MW satellites could have formed after the epoch of reionization. Our models also predict a significantly larger dispersion in M_300 values than observed and a number of luminous satellites with M_300 as low as 10^6 solar masses.

preprint2008arXiv

Formation and evolution of galaxy dark matter halos and their substructure

We use the ``Via Lactea'' simulation to study the co-evolution of a Milky Way-size LambdaCDM halo and its subhalo population. While most of the host halo mass is accreted over the first 6 Gyr in a series of major mergers, the physical mass distribution [not M_vir(z)] remains practically constant since z=1. The same is true in a large sample of LambdaCDM galaxy halos. Subhalo mass loss peaks between the turnaround and virialization epochs of a given mass shell, and declines afterwards. 97% of the z=1 subhalos have a surviving bound remnant at the present epoch. The retained mass fraction is larger for initially lighter subhalos: satellites with maximum circular velocities Vmax=10 km/s at z=1 have today about 40% of their mass back then. At the first pericenter passage a larger average mass fraction is lost than during each following orbit. Tides remove mass in substructure from the outside in, leading to higher concentrations compared to field halos of the same mass. This effect, combined with the earlier formation epoch of the inner satellites, results in strongly increasing subhalo concentrations towards the Galactic center. We present individual evolutionary tracks and present-day properties of the likely hosts of the dwarf satellites around the Milky Way. The formation histories of ``field halos'' that lie today beyond the Via Lactea host are found to strongly depend on the density of their environment. This is caused by tidal mass loss that affects many field halos on eccentric orbits.

preprint2008arXiv

Quantifying the heart of darkness with GHALO - a multi-billion particle simulation of our galactic halo

We perform a series of simulations of a Galactic mass dark matter halo at different resolutions, our largest uses over three billion particles and has a mass resolution of 1000 M_sun. We quantify the structural properties of the inner dark matter distribution and study how they depend on numerical resolution. We can measure the density profile to a distance of 120 pc (0.05% of R_vir) where the logarithmic slope is -0.8 and -1.4 at (0.5% of R_vir). We propose a new two parameter fitting function that has a linearly varying logarithmic density gradient which fits the GHALO and VL2 density profiles extremely well. Convergence in the density profile and the halo shape scales as N^(-1/3), but the shape converges at a radius three times larger at which point the halo becomes more spherical due to numerical resolution. The six dimensional phase-space profile is dominated by the presence of the substructures and does not follow a power law, except in the smooth under-resolved inner few kpc.

preprint2007arXiv

Redefining the Missing Satellites Problem

Numerical simulations of Milky-Way size Cold Dark Matter (CDM) halos predict a steeply rising mass function of small dark matter subhalos and a substructure count that greatly outnumbers the observed satellites of the Milky Way. Several proposed explanations exist, but detailed comparison between theory and observation in terms of the maximum circular velocity (Vmax) of the subhalos is hampered by the fact that Vmax for satellite halos is poorly constrained. We present comprehensive mass models for the well-known Milky Way dwarf satellites, and derive likelihood functions to show that their masses within 0.6 kpc (M_0.6) are strongly constrained by the present data. We show that the M_0.6 mass function of luminous satellite halos is flat between ~ 10^7 and 10^8 M_\odot. We use the ``Via Lactea'' N-body simulation to show that the M_0.6 mass function of CDM subhalos is steeply rising over this range. We rule out the hypothesis that the 11 well-known satellites of the Milky Way are hosted by the 11 most massive subhalos. We show that models where the brightest satellites correspond to the earliest forming subhalos or the most massive accreted objects both reproduce the observed mass function. A similar analysis with the newly-discovered dwarf satellites will further test these scenarios and provide powerful constraints on the CDM small-scale power spectrum and warm dark matter models.

preprint2006arXiv

Dark matter substructure and gamma-ray annihilation in the Milky Way halo

We present initial results from ``Via Lactea'', the highest resolution simulation to date of Galactic CDM substructure. It follows the formation of a Milky Way-size halo with Mvir=1.8x10^12 Msun in a WMAP 3-year cosmology, using 234 million particles. Over 10,000 subhalos can be identified at z=0: Their cumulative mass function is well-fit by N(>Msub)= 0.0064 (Msub/Mvir)^(-1) down to Msun=4x10^6 Msun. The total mass fraction in subhalos is 5.3%, while the fraction of surface mass density in substructure within a projected distance of 10 kpc from the halo center is 0.3%. Because of the significant contribution from the smallest resolved subhalos, these fractions have not converged yet. Sub-substructure is apparent in all the larger satellites, and a few dark matter lumps are resolved even in the solar vicinity. The number of dark satellites with peak circular velocities above 10 km/s (5 km/s) is 124 (812): of these, 5 (26) are found within 0.1 Rvir, a region that appeared practically smooth in previous simulations. The neutralino self-annihilation gamma-ray emission from dark matter clumps is approximately constant per subhalo mass decade. Therefore, while in our run the contribution of substructure to the gamma-ray luminosity of the Galactic halo amounts to only 40% of the total spherically-averaged smooth signal, we expect this fraction to grow significantly as resolution is increased further. An all-sky map of the expected annihilation gamma-ray flux reaching a fiducial observer at 8 kpc from the Galactic center shows that at the current resolution a small number of subhalos start to be bright enough to be visible against the background from the smooth density field surrounding the observer.

preprint2006arXiv

Multi-scale simulations of merging galaxies with supermassive black holes

We present the results of the first multi-scale N-Body+SPH simulations of merging galaxies containing central supermassive black holes (SMBHs) and having a spatial resolution of only a few parsecs. Strong gas inflows associated with equal-mass mergers produce non-axisymmetric nuclear disks with masses of order $10^9 M_{\odot}$, resolved by about $10^6$ SPH particles. Such disks have sizes of several hundred parsecs but most of their mass is concentrated within less than $50$ pc. We find that a close SMBH pair forms after the merger. The separation of the two SMBHs then shrinks further owing to dynamical friction against the predominantly gaseous background. The orbits of the SMBHs decay down to the minimum resolvable scale in a few million years for an ambient gas temperature and density typical of a region undergoing a starburst. These results suggest the initial conditions necessary for the eventual coalescence of the two holes arise naturally from the merging of two equal-mass galaxies whose structure and orbits are consistent with the predictions of the $Λ$CDM model. Our findings have important implications for planned gravitational wave detection experiments such as {\it LISA}.

preprint2005arXiv

Fate of Supermassive Black Holes and Evolution of the M(BH)-sigma Relation in Merging Galaxies: The Effect of Gaseous Dissipation

We analyze the effect of dissipation on the orbital evolution of supermassive black holes (SMBHs) using high-resolution self-consistent gasdynamical simulations of binary equal- and unequal-mass mergers of disk galaxies. The galaxy models are consistent with the LCDM paradigm of structure formation and the simulations include the effects of radiative cooling and star formation. We find that equal-mass mergers always lead to the formation of a close SMBH pair at the center of the remnant with separations limited solely by the adopted force resolution of ~ 100 pc. Instead, the final SMBH separation in unequal-mass mergers depends sensitively on how the central structure of the merging galaxies is modified by dissipation. In the absence of dissipation, the satellite galaxy can be entirely disrupted before the merger is completed leaving its SMBH wandering at a distance too far from the center of the remnant for the formation of a close pair. In contrast, we show that gas cooling facilitates the pairing process by increasing the resilience of the companion galaxy to tidal disruption. Moreover, we demonstrate that merging disk galaxies constructed to obey the M(BH)-sigma relation, move relative to it depending on whether they undergo a dissipational or collisionless merger, regardless of the mass ratio of the merging systems. Collisionless simulations reveal that remnants tend to move away from the mean relation highlighting the role of gas-poor mergers as a possible source of scatter. In dissipational mergers, the interplay between strong gas inflows associated with the formation of massive nuclear disks and the consumption of gas by star formation provides the necessary fuel to the SMBHs and allows the merger remnants to satisfy the relation.

preprint2004arXiv

A New Non-Parametric Approach to Galaxy Morphological Classification

We present two new non-parametric methods for quantifying galaxy morphology: the relative distribution of the galaxy pixel flux values (the Gini coefficient or G) and the second-order moment of the brightest 20% of the galaxy's flux (M20). We test the robustness of G and M20 to decreasing signal-to-noise and spatial resolution, and find that both measures are reliable to within 10% at average signal-to-noise per pixel greater than 3 and resolutions better than 1000 pc and 500 pc, respectively. We have measured G and M20, as well as concentration (C), asymmetry (A), and clumpiness (S) in the rest-frame near-ultraviolet/optical wavelengths for 150 bright local "normal" Hubble type galaxies (E-Sd) galaxies and 104 0.05 < z < 0.25 ultra-luminous infrared galaxies (ULIRGs).We find that most local galaxies follow a tight sequence in G-M20-C, where early-types have high G and C and low M20 and late-type spirals have lower G and C and higher M20. The majority of ULIRGs lie above the normal galaxy G-M20 sequence, due to their high G and M20 values. Their high Gini coefficients arise from very bright nuclei, while the high second-order moments are produced by multiple nuclei and bright tidal tails. All of these features are signatures of recent and on-going mergers and interactions. We also find that in combination with A and S, G is more effective than C at distinguishing ULIRGs from the "normal" Hubble-types. Finally, we measure the morphologies of 45 1.7 < z < 3.8 galaxies from HST NICMOS observations of the Hubble Deep Field North. We find that many of the z $\sim$ 2 galaxies possess G and A higher than expected from degraded images of local elliptical and spiral galaxies, and have morphologies more like low-redshift single nucleus ULIRGs.

preprint2004arXiv

The distribution and cosmic evolution of massive black hole spins

We study the expected distribution of massive black hole (MBH) spins and its evolution with cosmic time in the context of hierarchical galaxy formation theories. Our model uses Monte Carlo realizations of the merger hierarchy in a LCDM cosmology, coupled to semi-analytical recipes, to follow the merger history of dark matter halos, the dynamics of the MBHs they host, and their growth via gas accretion and binary coalescences. The coalescence of comparable mass holes increases the spin of MBHs, while the capture of smaller companions in randomly-oriented orbits acts to spin holes down. We find that, given the distribution of MBH binary mass ratios in hierarchical models, binary coalescences alone do not lead to a systematic spin-up or spin-down of MBHs with time: the spin distribution retains memory of its initial conditions. By contrast, because of the Bardeen-Petterson effect, gas accretion via a thin disk tends to spin holes up even if the direction of the spin axis changes randomly in time. In our models, accretion dominates over black hole captures and efficiently spins holes up. The spin distribution is heavily skewed towards fast-rotating Kerr holes, is already in place at early epochs, and does not change much below redshift 5. If accretion is via a thin disk, about 70% of all MBHs are maximally rotating and have radiative efficiencies approaching 30% (assuming a "standard'' spin-efficiency conversion). Even in the conservative case where accretion is via a geometrically thick disk, about 80% of all MBHs have spin parameters a/m > 0.8 and accretion efficiencies > 12%. Rapidly spinning holes with high radiative efficiencies may satisfy constraints based on comparing the local MBH mass density with the mass density inferred from luminous quasars (Soltan's argument).

preprint2003arXiv

The Size Evolution of High Redshift Galaxies

Hubble Space Telescope images of high-redshift galaxies selected via color and photometric redshifts are used to examine the size and axial-ratio distribution of galaxies as a function of redshift at lookback times t > 8 Gyr. These parameters are measured at rest-frame UV wavelengths (1200-2000 Angstroms) on images with a rest-frame resolution of less than 0.8 kpc. Galaxy radii are found to scale with redshift approximately as the inverse of the Hubble parameter H(z). This is in accord with the theoretical expectation that the typical sizes of the luminous parts of galaxies should track the expected evolution in the virial radius of dark-matter halos. The mean ratio of semi-major to semi-minor axis for a bright well-resolved sample of galaxies at z ~ 4 is b/a = 0.65, suggesting that these Lyman break galaxies are not drawn from a spheroidal population. However the median concentration index of this sample is C = 3.5, which is closer to the typical concentration indices, C ~ 4, of nearby elliptical galaxies than to the values, C<2 for local disk galaxies of type Sb and later.

preprint2001arXiv

Compound gravitational lensing as a probe of dark matter substructure within galaxy halos

We show how observations of multiply-imaged quasars at high redshift can be used as a probe of dark matter clumps (subhalos with masses ~ 10^9 solar masses) within the virialized extent of more massive lensing halos. A large abundance of such satellites is predicted by numerical simulations of galaxy formation in cold dark matter (CDM) cosmogonies. Small-scale structure within galaxy halos affects the flux ratios of the images without appreciably changing their positions. We use numerical simulations to quantify the effect of dark matter substructure on the distribution of magnifications, and find that the magnification ratio of a typical image pair will deviate significantly from the value predicted by a smooth lensing potential if, near the Einstein radius, only a few percent of the lens surface density is contained in subhalos. The angular size of the continuum source dictates the range of subclump masses that can have a detectable effect: to avoid confusion with gravitational microlensing caused by stars in the lens galaxy, the background source must be larger than the optical continuum-emitting region of a QSO. We also find that substructure will cause distortions to images on milli-arcsecond scales and bias the distribution of QSO magnification ratios -- two other possible methods of detection.

preprint2001arXiv

Early Metal Enrichment of the Intergalactic Medium by Pregalactic Outflows

We assess supernova (SN)-driven pregalactic outflows as a mechanism for distributing the product of stellar nucleosynthesis over large cosmological volumes prior to the reionization epoch. SN ejecta will escape the grasp of halos with virial temperatures T_vir>10^{4.3} K (corresponding to masses M>10^8 h^{-1} M_sun at redshift z=9 when they collapse from 2-sigma fluctuations) if rapid cooling can take place, and a significant fraction of their baryonic mass is converted into stars over a dynamical timescale. We study the evolution of SN-driven bubbles as they blow out from subgalactic halos and propagate into the intergalactic medium (IGM), and show that to lift the halo gas out of the potential well the energy injection must continue at least until blow-away occurs. The collective explosive output of about ten thousands SNe per halo at these early epochs could pollute vast regions of intergalactic space to a mean metallicity <Z>=Omega_Z/Omega_b>0.003 (comparable to the levels observed in the Lyman-alpha forest at z=3) without much perturbing the IGM hydrodynamically, i.e. producing large variations of the baryons relative to the dark matter. The volume filling factor of the ejecta is higher than 20% if the star formation efficiency is of order 10%. When the filling factor of the ejecta becomes significant, enriched material will typically be at a higher adiabat than expected from photoionization.

preprint1996arXiv

High Redshift Galaxies in the Hubble Deep Field. Color Selection and Star Formation History to z=4

The Lyman decrement associated with the cumulative effect of HI in QSO absorption systems along the line of sight provides a distinctive feature for identifying galaxies at z>2.5. The Hubble Deep Field (HDF) observations offer the opportunity to exploit the ubiquitous effect of intergalactic absorption and obtain useful statistical constraints on the redshift distribution of galaxies considerably fainter than current spectroscopic limits. We model the HI cosmic opacity as a function of redshift, including scattering in resonant lines of the Lyman series and Lyman-continuum absorption, and use stellar population synthesis models with a wide variety of ages, metallicities, dust contents, and redshifts, to derive color selection criteria that provide a robust separation between high redshift and low redshift galaxies. From the HDF images we construct a sample of star-forming galaxies at 2<z<4.5. While none of the more than 60 objects in the HDF having known Keck/LRIS spectroscopic redshifts in the range 0<z<1.4 is found to contaminate our high-redshift sample, our color criteria are able to efficiently select the 2.6<z<3.2 galaxies identified by Steidel et al. (1996b). The ultraviolet (and blue) dropout technique opens up the possibility of investigating cosmic star and element formation in the early universe.

Piero Madau

What is connected

Connect this record

See the researcher in context

Building this map preview

73 published item(s)

Dynamical Dark Energy Imprints in the Lyman-Alpha Forest

Inferring the Thermal History of the Intergalactic Medium from the Properties of the Hydrogen and Helium Lyman-alpha Forest

The Dawn of Disk Formation in a Milky Way-sized Galaxy Halo: Thin Stellar Disks at $z > 4$

Deep Realistic Extragalactic Model (DREaM) Galaxy Catalogs: Predictions for a Roman Ultra-Deep Field

Formation of the first stars and black holes

Global torques and stochasticity as the drivers of massive black hole pairing in the young Universe

Globular Cluster Formation from Colliding Substructure

The impact of Lyman-$α$ emission line heating and cooling on the cosmic dawn 21-cm signal

Momentum Injection by Clustered Supernovae: Testing Subgrid Feedback Prescriptions

Bar-driven evolution and quenching of spiral galaxies in cosmological simulations

Clumpy high-z galaxies as a testbed for feedback-regulated galaxy formation

Cold Dark Matter Substructures in Early-Type Galaxy Halos

Diffuse coronae in cosmological simulations of Milky Way-sized galaxies

Dispersal of tidal debris in a Milky-Way-sized dark matter halo

Enhanced Momentum Feedback from Clustered Supernovae

Black Hole Starvation and Bulge Evolution in a Milky Way-like Galaxy

Cosmic Reionization after Planck: Could Quasars Do It All?

Evidence of patchy hydrogen reionization from an extreme Ly$α$ trough below redshift six

Signatures of Kinematic Substructure in the Galactic Stellar Halo

Simulating tidal streams in a high resolution dark matter halo

The History of R-Process Enrichment in the Milky Way

The photoheating of the intergalactic medium in synthesis models of the UV background

A cosmic web filament revealed in Lyman-alpha emission around a luminous high-redshift quasar

Building Late-Type Spiral Galaxies by In-Situ and Ex-Situ Star Formation

Carbon-Enhanced Metal-Poor Stars: Relics from the Dark Ages

Cosmic Star Formation History

Dark Matter Heating and Early Core Formation in Dwarf Galaxies

Faint dwarfs as a test of DM models: WDM vs. CDM

Reversal of Fortune: Increased Star Formation Efficiencies in the Early Histories of Dwarf Galaxies?

Super-Critical Growth of Massive Black Holes from Stellar-Mass Seeds

The Baryon Cycle of Dwarf Galaxies: Dark, Bursty, Gas-Rich Polluters

The Dark Matter Contribution to Galactic Diffuse Gamma Ray Emission

The Distribution of Dark Matter in the Milky Way's Disk

The Photon Underproduction Crisis

Towards a Unified Description of the Intergalactic Medium at Redshift z~2.5

A "Light," Centrally-Concentrated Milky Way Halo?

A Population of Relic Intermediate-Mass Black Holes in the Halo of the Milky Way

An Off-center Density Peak in the Milky Way's Dark Matter Halo?

Direct Gravitational Imaging of Intermediate Mass Black Holes in Extragalactic Halos

Dwarf Galaxy Formation with H2-Regulated Star Formation: II. Gas-Rich Dark Galaxies at Redshift 2.5

Inside Out and Upside Down: Tracing the Assembly of a Simulated Disk Galaxy Using Mono-Age Stellar Populations

The circumgalactic medium of massive galaxies at z~3: a test for stellar feedback, galactic outflows, and cold streams

The HST/ACS+WFC3 Survey for Lyman Limit Systems II: Science

Insight Into the Formation of the Milky Way Through Cold Halo Substructure. III. Statistical Chemical Tagging in the Smooth Halo

Pseudobulge Formation as a Dynamical Rather than a Secular Process

The origin of metals in the circum-galactic medium of massive galaxies at z=3

A galaxy as the source of a Civ absorption system close to the epoch of reionization

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey

CANDELS: The Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey - The Hubble Space Telescope Observations, Imaging Data Products and Mosaics

Delayed Enrichment by Unseen Galaxies: Explaining the Rapid Rise in IGM CIV Absorption from z = 6-5

Forming Realistic Late-Type Spirals in a LCDM Universe: The Eris Simulation

New Population Synthesis Models of the Cosmic UV Background

On the assembly of the Milky Way dwarf satellites and their common mass scale

Radiative transfer in a clumpy universe: IV. New synthesis models of the cosmic UV/X-ray background

Fermi-LAT Sensitivity to Dark Matter Annihilation in Via Lactea II Substructure

Recoiling Massive Black Holes in Gas-Rich Galaxy Mergers

Tidal stellar disruptions by massive black hole pairs: II. Decaying binaries

At the heart of the matter: the origin of bulgeless dwarf galaxies and Dark Matter cores

GeV Gamma-Ray Attenuation and the High-Redshift UV Background

He II absorption and the sawtooth spectrum of the cosmic far-UV background

Probing the Epoch of Reionization with Milky Way Satellites

Formation and evolution of galaxy dark matter halos and their substructure

Quantifying the heart of darkness with GHALO - a multi-billion particle simulation of our galactic halo

Redefining the Missing Satellites Problem

Dark matter substructure and gamma-ray annihilation in the Milky Way halo

Multi-scale simulations of merging galaxies with supermassive black holes

Fate of Supermassive Black Holes and Evolution of the M(BH)-sigma Relation in Merging Galaxies: The Effect of Gaseous Dissipation

A New Non-Parametric Approach to Galaxy Morphological Classification

The distribution and cosmic evolution of massive black hole spins

The Size Evolution of High Redshift Galaxies

Compound gravitational lensing as a probe of dark matter substructure within galaxy halos

Early Metal Enrichment of the Intergalactic Medium by Pregalactic Outflows

High Redshift Galaxies in the Hubble Deep Field. Color Selection and Star Formation History to z=4