Source author record

Massimo Ricotti

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

32works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

32 published item(s)

preprint2026arXiv

PEARLS: 21 Transients Found in the Three-Epoch NIRCam Observations in the Continuous Viewing Zone of the James Webb Space Telescope

We present 21 transients from our three-epoch, four-band NIRCam observations covering 14.16 arcmin^2 in the Spitzer IRAC Dark Field (IDF), taken by the JWST Prime Extragalactic Areas for Reionization and Lensing Science program with a time cadence of ~6 months. A separate Hubble Space Telescope program provided Advanced Camera for Surveys optical imaging contemporaneous with the second and third epochs of the NIRCam observations. The NIRSpec spectroscopy on three transients confirmed a Type Ia supernova at z=1.63 and the host galaxies of the other two at z=2.64 and 1.90, respectively. Combining these with the photometric redshifts (z_ph) of the host galaxies in the rest of the sample, we find that the transients are in either a "mid-z" group at z>1.6 with M_V < -16.0 mag or a "low-z" group at z < 0.4 with M_H > -14.0 mag. The mid-z transients are consistent with supernovae. In contrast, the low-z transients' luminosities fall in the range of the so-called "gap transients" between supernovae and novae. However, this latter conclusion is only tentative due to possible catastrophic failures in z_ph that could bias them to low-z. Conversely, if they are indeed at z < 0.4, it would be worth studying similar transients in the future. Our work further demonstrates the power of NIRCam in transient science and also shows that it would be more fruitful to carry out a long-term monitoring program with more passbands, a higher cadence and prompt follw-up spectroscopy. Being in the continuous viewing zone of the JWST, the IDF is an ideal field for this purpose.

preprint2026arXiv

VENUS: Two Faint Little Red Dots Separated by $\sim70\,\mathrm{pc}$ Hidden in a Single Lensed Galaxy at $z\sim7$

We report the identification of a pair of faint little red dots (LRDs), dubbed Red Eyes, in a strongly-lensed galaxy at $z\sim7$ behind the PLCKG004.5-10.5 cluster, identified from the JWST Treasury program VENUS. Red Eyes are spatially resolved on the image plane with distinct colors, while the critical curve lies far north of Red Eyes, clearly requiring two different LRDs rather than a single LRD. Red Eyes is an extremely close pair of LRDs separated by $\sim70\,\mathrm{pc}$ in the source plane with a magnification of $μ\sim20$, which consistently explains another counter-image detected to the north-west. Red Eyes is hosted in a typical star-forming galaxy with $M_{\mathrm{UV,int}}\sim -19$, but its own UV emission is very faint ($M_{\mathrm{UV,int}} \gtrsim -16$). Moreover, Red Eyes does not reside at the galaxy center but lies at an offset position of approximately one effective radius $R_{\mathrm{e}}$ away from the galaxy center. If observed without lensing, Red Eyes would appear as a typical star-forming galaxy at $z\sim 7$ with $M_{\mathrm{UV}}\sim -19$, showing no apparent LRD signatures in either morphology or SED. These results suggest that multiple off-center LRDs, similar to Red Eyes, may be commonly hidden in a typical high-$z$ star-forming galaxy. In this case, various plausible scenarios may emerge, one of which is that intermediate-mass black holes (IMBHs) with $M_\mathrm{BH}\sim10^{4\text{--}6}\,M_\odot$ may form in star clusters on a stellar disk and contribute to the growth of the central supermassive black hole via mergers, with some IMBHs detectable as luminous LRDs in a sufficiently active and massive phase.

preprint2022arXiv

Ghostly Stellar Haloes and their Relationship to Ultra-faint Dwarfs

Ghostly stellar haloes are extended haloes of stars composed solely of debris of pre-reionization fossil galaxies and should exist in dwarf galaxies with total masses $<10^{10}$ M$_\odot$. Fossil galaxies are even smaller mass dwarf galaxies that stopped forming stars after the epoch of reionization and have been identified in the Local Group as the ultra-faint dwarf satellites. Using cosmological N-body simulations we present an empirical model for the shapes and masses of ghostly stellar haloes. We compare the model to available observations of stellar haloes in six isolated dwarf galaxies in the Local Group (Leo T, Leo A, IC 10, WLM, IC 1613, NGC 6822) to infer the star formation efficiency in dwarf galaxies at the epoch of reionization. We find an efficiency of star formation in dark matter haloes with masses $10^6 - 10^8$ M$_\odot$ at $z\sim7$ in rough agreement with independent methods using data on the luminosity function of ultra-faint dwarf galaxies but systematically higher by a factor of 3-5. The systematic uncertainty of our results is still large, mainly because available observations of stellar halo profiles do not extend over a sufficiently large distance from the center of the host dwarf galaxy. Additional observations, easily within reach of current telescopes, can significantly improve the accuracy of this method and can also be used to constrain the present day dark matter masses of dwarf galaxies in the Local Group. Our method is based on a set of observations never used before, hence it is a new independent test of models of hierarchical galaxy formation.

preprint2020arXiv

Simulating Star Clusters Across Cosmic Time: II. Escape Fraction of Ionizing Photons from Molecular Clouds

We calculate the hydrogen and helium-ionizing radiation escaping star-forming molecular clouds, as a function of the star cluster mass and compactness, using a set of high-resolution radiation-magneto-hydrodynamic simulations of star formation in self-gravitating, turbulent molecular clouds. In these simulations, presented in He, Ricotti and Geen (2019), the formation of individual massive stars are well resolved, and their UV radiation feedback and lifetime on the main sequence are modeled self-consistently. We find that the escape fraction of ionizing radiation from molecular clouds, $\langle f_{\rm esc}^{\scriptscriptstyle \rm MC}\rangle$, decreases with increasing mass of the star cluster and with decreasing compactness. Molecular clouds with densities typically found in the local Universe have negligible $\langle f_{\rm esc}^{\scriptscriptstyle \rm MC}\rangle$, ranging between $0.5\%$ to $5\%$. Ten times denser molecular clouds have $\langle f_{\rm esc}^{\scriptscriptstyle \rm MC}\rangle \approx 10\%-20\%$, while $100\times$ denser clouds, which produce globular cluster progenitors, have $\langle f_{\rm esc}^{\scriptscriptstyle \rm MC}\rangle \approx 20\%-60\%$. We find that $\langle f_{\rm esc}^{\scriptscriptstyle \rm MC}\rangle$ increases with decreasing gas metallicity, even when ignoring dust extinction, due to stronger radiation feedback. However, the total number of escaping ionizing photons decreases with decreasing metallicity because the star formation efficiency is reduced. We conclude that the sources of reionization at $z>6$ must have been very compact star clusters forming in molecular clouds about $100\times$ denser than in today's Universe, which leads to a significant production of old globular clusters progenitors.

preprint2020arXiv

Structure and Instability of the Ionization Fronts around Moving Black Holes

In this paper we focus on understanding the physical processes that lead to stable or unstable ionization fronts (I-fronts) observed in simulations of moving black holes (BHs). The front instability may trigger bursts of gas accretion, rendering the BH significantly more luminous than at steady-state. We perform a series of idealized three dimensional radiation hydrodynamics simulations resolving the I-fronts around BHs of mass $M_\mathrm{BH}$ and velocity $v_\infty$ accreting from a medium of density $n_\mathrm{H}$. The I-front, with radius $R_\mathrm{I}$, transitions from D-type to R-type as the BH velocity becomes larger than a critical value $v_\mathrm{R}\sim 40\,\mathrm{km/s}$. The D-type front is preceded by a bow-shock of thickness $ΔR_\mathrm{I}$ that decreases as $v_\infty$ approaches $v_\mathrm{R}$. We find that both D-type and R-type fronts can be unstable given the following two conditions: i) for D-type fronts the shell thickness must be $ΔR_\mathrm{I}/R_\mathrm{I}<0.05$ (i.e., $v_\infty \gtrsim 20\,\mathrm{km/s}$.), while no similar restriction holds for R-type fronts; ii) the temperature jump across the I-front must be $T_\mathrm{II}/T_\mathrm{I}>3$. This second condition is satisfied if $T_\mathrm{I}<5000\,\mathrm{K}$ or if $n_\mathrm{H}\,M_\mathrm{BH} \gtrsim 10^6\,M_\odot\,\mathrm{cm^{-3}}$. Due to X-ray pre-heating typically $T_\mathrm{I} \sim 10^4\,\mathrm{K}$, unless the D-type shell is optically thick to X-rays, which also happens when $n_\mathrm{H}\,M_\mathrm{BH}$ is greater than a metallicity-dependent critical value. We thus conclude that I-fronts around BHs are unstable only for relatively massive BHs moving trough very dense molecular clouds. We briefly discuss the observational consequences of the X-ray luminosity bursts likely associated with this instability.

preprint2016arXiv

A Common Origin for Globular Clusters and Ultra-faint Dwarfs in Simulations of the First Galaxies

In this paper, the first in a series on galaxy formation before reionization, we focus on understanding what determines the size and morphology of stellar objects in the first low mass galaxies, using parsec- scale cosmological simulations performed with an adaptive mesh hydrodynamics code. Although the dense gas in which stars are formed tends to have a disk structure, stars are found in spheroids with little rotation. Halos with masses between 10^6 M_sun and 5x10^8 M_sun form stars stochastically, with stellar masses in the range 10^4 M_sun to 2x10^6 M_sun. Nearly independent of stellar mass, we observe a large range of half-light radii for the stars, from a few parsecs to a few hundred parsecs and surface brightnesses and mass-to-light ratios ranging from those typical of globular clusters to ultra-faint dwarfs. In our simulations, stars form in dense stellar clusters with high gas-to-star conversion efficiencies and rather uniform metallicities. A fraction of these clusters remain bound after the gas is removed by feedback, but others are destroyed, and their stars, which typically have velocity dispersions of 20 to 40 km/s, expand until they become bound by the dark matter halo. We thus speculate that the stars in ultra-faint dwarf galaxies may show kinematic and chemical signatures consistent with their origin in a few distinct stellar clusters. On the other hand, some globular clusters may form at the center of primordial dwarf galaxies and may contain dark matter, perhaps detectable in the outer parts.

preprint2016arXiv

Modeling Reionization in a Bursty Universe

We present semi-analytic models of the epoch of reionization focusing on the differences between continuous and bursty star formation (SF). Our model utilizes physically motivated analytic fits to 1D radiative transfer simulations of HII regions around dark matter halos in a representative cosmic volume. Constraining our simulations with observed and extrapolated UV luminosity functions of high redshift galaxies, we find that for a fixed halo mass, stellar populations forming in bursty models produce larger HII regions which leave behind long-lived relic HII regions which are able to maintain partial ionization in the intergalactic medium (IGM) in a manner similar to an early X-ray background. The overall effect is a significant increase in the optical depth of the IGM, $τ_e$, and a milder increase of the redshift of reionization. To produce $τ_e =0.066$ observed by Planck and complete reionization by redshift $z_{\rm re}\sim 6$, models with bursty SF require an escape fraction $f_{\rm esc}\sim 2\%-10\%$ that is 2-10 times lower than $f_{\rm esc}\sim 17\%$ found assuming continuous SF and is consistent with upper limits on $f_{\rm esc}$ from observations at $z=0$ and $z\sim 1.3-6$. The ionizing photon budget needed to reproduce the observed $τ_e$ and $z_{\rm re}$ depends on the period and duty cycle of the bursts of SF and the temperature of the neutral IGM. These results suggest that the tension between observed and predicted ionizing photon budget for reionization can be alleviated if reionization is driven by short bursts of SF, perhaps relating to the formation of Population~III stars and compact star clusters such as proto-globular clusters.

preprint2016arXiv

The Lost Dwarfs of Centaurus A and the Formation of its Dark Globular Clusters

We present theoretical constraints for the formation of the newly discovered dark star clusters (DSCs) with high mass-to-light (M/L) ratios, from Taylor et al (2015). These compact stellar systems photometrically resemble globular clusters (GCs) but have dynamical M/L ratios of ~ 10 - 100, closer to the expectations for dwarf galaxies. The baryonic properties of the dark star clusters (DSCs) suggest their host dark matter halos likely virialized at high redshift with M > 10^8 M_sun. We use a new set of high-resolution N-body simulations of Centaurus A to determine if there is a set of z=0 subhalos whose properties are in line with these observations. While we find such a set of subhalos, when we extrapolate the dark matter density profiles into the inner 20 pc, no dark matter halo associated with Centaurus A in our simulations, at any redshift, can replicate the extremely high central mass densities of the DSCs. Among the most likely options for explaining 10^5 - 10^7 M_sun within 10 pc diameter subhalos is the presence of a central massive black hole. We, therefore, propose that the DSCs are remnant cusps of stellar systems surrounding the central black holes of dwarf galaxies which have been almost completely destroyed by interactions with Centaurus A.

preprint2016arXiv

X-ray Twinkles and Pop III Stars

Pop III stars are typically massive stars of primordial composition forming at the centers of the first collapsed dark matter structures. Here we estimate the optimal X-ray emission in the early universe for promoting the formation of Pop III stars. This is important in determining the number of dwarf galaxies formed before reionization and their fossils in the local universe, as well as the number of intermediate-mass seed black holes. A mean X-ray emission per source above the optimal level reduces the number of Pop III stars because of the increased Jeans mass of the intergalactic medium (IGM), while a lower emission suppresses the formation rate of H2 preventing or delaying star formation in dark matter minihalos above the Jeans mass. The build up of the H2 dissociating background is slower than the X-ray background due to the shielding effect of resonant hydrogen Lyman lines. Hence, the nearly unavoidable X-ray emission from supernova remnants of Pop III stars is sufficient to boost their number to few tens per comoving Mpc^3 by redshift z~15. We find that there is a critical X-ray to UV energy ratio emitted per source that produces a universe where the number of Pop III stars is largest: 400 per comoving- Mpc^3. This critical ratio is very close to the one provided by 20-40 M_sun Pop III stars exploding as hypernovae. High mass X-ray binaries in dwarf galaxies are far less effective at increasing the number of Pop III stars than normal supernova remnants, we thus conclude that supernovae drove the formation of Pop III stars.

preprint2015arXiv

Bulge-driven Fueling of Seed Black Holes

We examine radiation-regulated accretion onto intermediate-mass and massive black holes (BHs) embedded in a bulge component. Using spherically symmetric one-dimensional radiation-hydrodynamics simulations, we track the growth of BHs accreting from a cold, neutral gas reservoir with temperature T=10^4 K. We find that the accretion rate of BHs embedded in bulges is proportional to r_{B,eff}/r_B, where r_{B,eff} is the increased effective Bondi radius that includes the gravitational potential of the bulge, and r_B is the Bondi radius of the BH. The radiative feedback from the BH suppresses the cold accretion rate to ~1 percent of the Bondi rate when a bulge is not considered. However, we find that the BH fueling rate increases rapidly when the bulge mass M_bulge is greater than the critical value of 10^6 M_sun and is proportional to M_bulge. Since the critical bulge mass is independent of the central BH mass M_{BH}, the growth rate of BHs with masses of 10^2, 10^4, and 10^6 M_sun exhibits distinct dependencies on the bulge-to-BH mass ratio. Our results imply that light seed BHs (<= 10^2 M_sun) which might be the remnants of the Pop III stars, cannot grow through accretion coevally with the early assembly of the bulge of the host galaxies until the bulge reaches the critical mass. However, massive BH seeds (>= 10^5 M_sun) that may form via direct collapse, are more likely to be embedded in a supercritical bulge and thus can grow efficiently coupling to the host galaxies and driving the early evolution of the M_{BH}-$σ$ relationship.

preprint2015arXiv

On the Physical Requirements for a Pre-Reionization Origin of the Unresolved Near-Infrared Background

The study of the Cosmic Near-Infrared Background (CIB) light after subtraction of resolved sources can push the limits of current observations and infer the level of galaxy and black hole activity in the early universe. However, disentangling the relative contribution from low- and high-redshift sources is not trivial. Spatial fluctuations of the CIB exhibit a clustering excess at angular scales $\sim 1^\circ$ whose origin has not been conclusively identified. We explore the likelihood that this signal is dominated by emission from galaxies and accreting black holes in the early Universe. We find that, if the first small mass galaxies have a normal IMF, the light of their ageing stars (fossils) integrated over cosmic time contributes a comparable amount to the CIB as their pre-reionization progenitors. However, the measured fluctuation signal is too large to be produced by galaxies at redshifts $z>8$ unless their star formation efficiencies are much larger than those inferred from the observed Lyman-dropout population. In order to produce the observed level of CIB fluctuation without violating constraints from galaxy counts and the electron optical depth of the IGM, minihalos at $z>12$ must form stars with efficiency $f_\star \gtrsim 0.1$ and, although a top-heavy IMF is preferred, have a very low escape fraction of ionizing radiation, $f_{\rm esc}<0.01$. If instead the CIB fluctuations are produced by high-$z$ black holes, one requires vigorous accretion in the early universe reaching $ρ_{\rm acc} \gtrsim 10^5M_\odot{\rm Mpc^{-3}}$ by $z\simeq 10$. This growth must stop by $z \sim 6$ and be significantly obscured not to overproduce the soft cosmic X-ray background (CXB) and its observed coherence with the CIB. We therefore find the range of suitable possibilities at high-$z$ to be narrow, but could possibly be widened by including additional physics and evolution at those epochs.

preprint2014arXiv

Clues on the Missing Sources of Reionization from Self-consistent Modeling of Milky Way and Dwarf Galaxy Globular Clusters

Globular clusters are unique tracers of ancient star formation. We determine the formation efficiencies of globular clusters across cosmic time by modeling the formation and dynamical evolution of the globular cluster population of a Milky Way type galaxy in hierarchical cosmology, using the merger tree from the Via Lactea II simulation. All of the models are constrained to reproduce the observed specific frequency and initial mass function of globular clusters in isolated dwarfs. Globular cluster orbits are then computed in a time varying gravitational potential after they are either accreted from a satellite halo or formed in situ, within the Milky Way halo. We find that the Galactocentric distances and metallicity distribution of globular clusters are very sensitive to the formation efficiencies of globular clusters as a function of redshift and halo mass. Our most accurate models reveal two distinct peaks in the globular cluster formation efficiency at z~2 and z~7-12 and prefer a formation efficiency that is mildly increasing with decreasing halo mass, the opposite of what expected for feedback-regulated star formation. This model accurately reproduces the positions, velocities, mass function, metallicity distribution, and age distribution of globular clusters in the Milky Way and predicts that ~ 40% formed in situ, within the Milky Way halo, while the other ~ 60% were accreted from about 20 satellite dwarf galaxies with Vc > 30 km/s, and about 29% or all globular clusters formed at redshifts z > 7. These results further strengthen the notion that globular cluster formation was an important mode of star formation in high-redshift galaxies and likely played a significant role in the reionization of the intergalactic medium

preprint2014arXiv

The Contribution of z<6 Sources to the Spatial Coherence in the Unresolved Cosmic near-Infrared and X-ray Backgrounds

A spatial clustering signal has been established in Spitzer/IRAC measurements of the unresolved Cosmic near-Infrared Background (CIB) out to large angular scales, ~1 deg. This CIB signal, while significantly exceeding the contribution from the remaining known galaxies, was further found to be coherent at a highly statistically significant level with the unresolved soft Cosmic X-ray Background (CXB). This measurement probes the unresolved CXB to very faint source levels using deep near-IR source subtraction. We study contributions from extragalactic populations at low to intermediate redshifts to the measured positive cross-power signal of the CIB fluctuations with the CXB. We model the X-ray emission from AGN, normal galaxies and hot gas residing in virialized structures, calculating their CXB contribution including their spatial coherence with all infrared emitting counterparts. We use a halo model framework to calculate the auto and cross-power spectra of the unresolved fluctuations based on the latest constraints of the halo occupation distribution and the biasing of AGN, galaxies and diffuse emission. At small angular scales (<1'), the 4.5mic vs 0.5-2 keV coherence can be explained by shot noise from galaxies and AGN. However, at large angular scales (~10 arcmin) we find that the net contribution from the modeled populations is only able to account for ~3% of the measured CIBxCXB cross-power. The discrepancy suggests that the CIBxCXB signal originates from the same unknown source population producing the CIB clustering signal out to ~1 deg.

preprint2014arXiv

The role of Compton heating in radiation-regulated accretion on to black holes

We investigate the role of Compton heating in radiation-regulated accretion on to black holes from a neutral dense medium using 1D radiation-hydrodynamic simulations. We focus on the relative effects of Compton-heating and photo-heating as a function of the spectral slope α, assuming a power-law spectrum in the energy range of 13.6 eV--100 keV. While Compton heating is dominant only close to the black hole, it can reduce the accretion rate to 0.1 % ($l \propto \dot{m}^2$ model)--0.01 % ($l \propto \dot{m}$ model) of the Bondi accretion rate when the BH radiation is hard (α ~ 1), where $l$ and $\dot{m}$ are the luminosity and accretion rate normalised by Eddington rates, respectively. The oscillatory behaviour otherwise typically seen in simulations with α > 1, become suppressed when α ~ 1 only for the $l \propto \dot{m}$ model. The relative importance of the Compton heating over photo-heating decreases and the oscillatory behaviour becomes stronger as the spectrum softens. When the spectrum is soft (α > 1.5), photo-heating prevails regardless of models making the effect of Compton heating negligible. On the scale of the ionization front, where the gas supply into the Strömgren sphere from large scale is regulated, photo-heating dominates. Our simulations show consistent results with the advection-dominated accretion flow ($l \propto \dot{m}^2$) where the accretion is inefficient and the spectrum is hard (α ~ 1).

preprint2013arXiv

Accretion onto Black Holes from Large Scales Regulated by Radiative Feedback. III. Enhanced Luminosity of Intermediate Mass Black Holes Moving at Supersonic Speeds

In this third paper of a series, we study the growth and luminosity of black holes (BHs) in motion with respect to their surrounding medium. We run a large set of two-dimensional (2D) axis-symmetric simulations to explore a large parameter space of initial conditions and formulate an analytical model for the accretion. Contrary to the case without radiation feedback, the accretion rate increases with increasing BH velocity $\vbh$ reaching a maximum value at $\vbh=2c_{\rm s,in}\sim 50$ km/s, where $c_{\rm s,in}$ is the sound speed inside the "cometary-shaped" \hii region around the BH, before decreasing as $\vbh^{-3}$ when the I-front becomes {\it R}-type (rarefied) and the accretion rate approaches the classical Bondi-Hoyle-Lyttleton solution. The increase of the accretion rate with $\vbh$ is produced by the formation of a {\it D}-type (dense) ionization front (I-front) preceded by a standing bow-shock that reduces the downstream gas velocity to transonic values. There is a range of densities and velocities where the dense shell is unstable producing periodic accretion rate peaks which can significantly increase the detectability of intermediate mass BHs. We find that the mean accretion rate for a moving BH is larger than that of a stationary BH of the same mass if the medium temperature is $T_\infty<10^4$ K. This result could be important for the growth of seed BHs in the multi-phase medium of the first galaxies and for building and early X-ray background that may affect the formation of the first galaxies and the reionization process.

preprint2013arXiv

Rayleigh-Taylor Instability at Ionization Fronts: Perturbation Analysis

The linear growth rate of the Rayleigh-Taylor instability (RTI) at ionization fronts is investigated via perturbation analysis in the limit of incompressible fluids. In agreement with previous numerical studies is found that absorption of ionizing radiation inside the HII region due to hydrogen recombinations suppresses the growth of instabilities. In the limit of a large density contrast at the ionization front the RTI growth rate has the simple analytical solution n=-nur+(nur^2+gk)^(1/2), where nur is the hydrogen recombination rate inside the HII region, k is the perturbation's wavenumber and g is the effective acceleration in the frame of reference of the front. Therefore, the growth of surface perturbations with wavelengths lambda >> lambda_{cr} = 2πg/nur^2 is suppressed by a factor (lambda_{cr}/4lambda)^(1/2) with respect to the non-radiative incompressible RTI. Implications on stellar and black hole feedback are briefly discussed.

preprint2013arXiv

Rayleigh-Taylor instability of ionization front around black holes

We examine the role of ionizing radiation emitted from black holes (BHs) in suppressing the growth of the Rayleigh-Taylor instability (RTI) across the ionization front (I-front) that forms when the gas fueling the BH is neutral. We use radiation-hydrodynamic simulations to show that the RTI is suppressed for non-accelerating fronts on all scales resolved in our simulations. A necessary condition for the stability of the I-front is that the radius of the Strömgren sphere is larger than the Bondi radius. When this condition is violated the I-front collapses producing an accretion luminosity burst. Transient growth of the RTI occurs only during the accretion burst when the effective acceleration in the frame of reference of the I-front increases significantly due to the rapid expansion of the Strömgren sphere.

preprint2012arXiv

Accretion onto Black Holes from Large Scales Regulated by Radiative Feedback. II. Growth Rate and Duty Cycle

This paper, the second of a series on radiation-regulated accretion onto black holes (BHs) from galactic scales, focuses on the effects of radiation pressure and angular momentum of the accreting gas. We simulate accretion onto intermediate-mass black holes, but we derive general scaling relationships that are solutions of the Bondi problem with radiation feedback valid for any mass of the BH $M_{\rm bh}$. Thermal pressure of the ionized sphere around the BH regulates the accretion rate producing periodic and short-lived luminosity bursts. We find that for ambient gas densities exceeding $n^{\rm cr}_{\rm H,\infty} \propto M_{\rm bh}^{-1}$, the period of the oscillations decreases rapidly and the duty cycle increases from 6%, in agreement with observations of the fraction of active galactic nuclei at $z\sim 3$, to 50%. The mean accretion rate becomes Eddington limited for $n_{\rm H,\infty}>n_{\rm H,\infty}^{\rm Edd} \simeq n_{\rm H,\infty}^{\rm cr} T_{\infty,4}^{-1}$ where $T_{\infty,4}$ is the gas temperature in units of $10^4$ K. In the sub-Eddington regime, the mean accretion rate onto BHs is about $1% T_{\infty,4}^{2.5}$ of the Bondi rate, thus is proportional to the thermal pressure of the ambient medium. The period of the oscillations coincides with depletion timescale of the gas inside the ionized bubble surrounding the BH. Gas depletion is dominated by a pressure gradient pushing the gas outward if $n_{\rm H,\infty}<n^{\rm cr}_{\rm H,\infty}$ and by accretion onto the BH otherwise. Generally, for $n_{\rm H,\infty}<n^{\rm cr}_{\rm H,\infty}$ angular momentum does not affect significantly the accretion rate and period of the oscillations.

preprint2012arXiv

Can the Near-IR Fluctuations Arise from Known Galaxy Populations?

Spatial Fluctuations in the Cosmic Infrared Background have now been measured out to sub-degree scales showing a strong clustering signal from unresolved sources. We attempt to explain these measurement by considering faint galaxy populations at z<6 as the underlying sources for this signal using 233 measured UV, optical and NIR luminosity functions (LF) from a variety of surveys covering a wide range of redshifts. We populate the lightcone and calculate the total emission redshifted into the near-IR bands in the observer frame and recover the observed optical and near-IR galaxy counts to a good accuracy. Using a halo model for the clustering of galaxies with an underlying LCDM density field, we find that fluctuations from known galaxy populations are unable to account for the large scale CIB clustering signal seen by HST/NICMOS, Spitzer/IRAC and AKARI/IRC and continue to diverge out to larger angular scales. Our purely empirical reconstruction shows that known galaxy populations are not responsible for the bulk of the fluctuation signal seen in the measurements and suggests an unknown population of very faint and highly clustered sources dominating the signal.

preprint2012arXiv

Reconstructing the Near-IR Background Fluctuations from known Galaxy Populations using Multiband Measurements of Luminosity Functions

We model fluctuations in the Cosmic Infrared Background (CIB) arising from known galaxy populations using 233 measured UV, optical and NIR luminosity functions (LF) from a variety of surveys spanning a wide range of redshifts. We compare best-fit Schechter parameters across the literature and find clear indication of evolution with redshift. Providing fitting formulae for the multi-band evolution of the LFs out to z~5, we calculate the total emission redshifted into the near-IR bands in the observer frame and recover the observed optical and near-IR galaxy counts to a good accuracy. Our empirical approach, in conjunction with a halo model describing the clustering of galaxies, allows us to compute the fluctuations of the unresolved CIB and compare the models to current measurements. We find that fluctuations from known galaxy populations are unable to account for the large scale CIB clustering signal seen by Spitzer/IRAC and AKARI/IRC and continue to diverge out to larger angular scales. This holds true even if the LFs are extrapolated out to faint magnitudes with a steep faint-end slope all the way to z=8. We also show that removing resolved sources to progressively fainter magnitude limits, isolates CIB fluctuations to increasingly higher redshifts. Our empirical approach suggests that known galaxy populations are not responsible for the bulk of the fluctuation signal seen in the measurements and favors a very faint population of highly clustered sources.

preprint2012arXiv

Two Epochs of Globular Cluster Formation from Deep Fields Luminosity Functions: Implications for Reionization and the Milky Way Satellites

The ages of globular clusters in our own Milky Way are known with precision of about \pm 1 Gyr, hence their formation at redshifts z>~3 and their role in hierarchical cosmology and the reionization of the intergalactic medium remain relatively undetermined. Here we analyze the effect of globular cluster formation on the observed rest-frame UV luminosity functions (LFs) and UV continuum slopes of high redshift galaxies in the Hubble Ultra Deep Fields. We find that the majority of present day globular clusters have formed during two distinct epochs: at redshifts z ~ 2-3 and at redshifts z>~6. The birth of proto-GC systems produce the steep, faint-end slopes of the galaxy LFs, and because the brightness of proto-GCs fades 5 Myrs after their formation, their blue colors are in excellent agreement with observations. Our results suggest that: i) the bulk of the old globular cluster population with estimated ages >~12 Gyr (about 50% of the total population), formed in the relatively massive dwarf galaxies at redshifts z>~6; ii) proto-GC formation was an important mode of star formation in those dwarf galaxies, and likely dominated the reionization process. Another consequence of this scenario is that some of the most massive Milky Way satellites may be faint and yet undiscovered because tidal stripping of a dominant GC population precedes significant stripping of the dark matter halos of these satellites. This scenario may alleviate some remaining tensions between CDM simulations and observations.

preprint2011arXiv

Accretion onto Intermediate Mass Black Holes Regulated by Radiative Feedback I. Parametric Study for Spherically Symmetric Accretion

We study the effect of radiative feedback on accretion onto intermediate mass black holes (IMBHs) using the hydrodynamical code ZEUS-MP with a radiative transfer algorithm. In this paper, the first of a series, we assume accretion from a uniformly dense gas with zero angular momentum and extremely low metallicity. Our 1D and 2D simulations explore how X-ray and UV radiation emitted near the black hole regulates the gas supply from large scales. Both 1D and 2D simulations show similar accretion rate and period between peaks in accretion, meaning that the hydro-instabilities that develop in 2D simulations do not affect the mean flow properties. We present a suite of simulations exploring accretion across a large parameter space, including different radiative efficiencies and radiation spectra, black hole masses, density and temperature, $T_\infty$, of the neighboring gas. In agreement with previous studies we find regular oscillatory behavior of the accretion rate, with duty cycle $\sim 6%$, mean accretion rate 3% $(T_{\infty}/10^4 {\rm K})^{2.5}$ of the Bondi rate and peak accretion $\sim 10$ times the mean for $T_{\infty}$ ranging between 3000K and 15000K. We derive parametric formulas for the period between bursts, the mean accretion rate and the peak luminosity of the bursts and thus provide a formulation of how feedback regulated accretion operates. The temperature profile of the hot ionized gas is crucial in determining the accretion rate, while the period of the bursts is proportional to the mean size of the Strömgren sphere and we find qualitatively different modes of accretion in the high vs. low density regimes. We also find that softer spectrum of radiation produces higher mean accretion rate.

preprint2011arXiv

Constraints on the Dark Matter Particle Mass from the Number of Milky Way Satellites

We have conducted N-body simulations of the growth of Milky Way-sized halos in cold and warm dark matter cosmologies. The number of dark matter satellites in our simulated Milky Ways decreases with decreasing mass of the dark matter particle. Assuming that the number of dark matter satellites exceeds or equals the number of observed satellites of the Milky Way we derive lower limits on the dark matter particle mass. We find with 95% confidence m_s > 13.3 keV for a sterile neutrino produced by the Dodelson and Widrow mechanism, m_s > 8.9 keV for the Shi and Fuller mechanism, m_s > 3.0 keV for the Higgs decay mechanism, and m_{WDM} > 2.3 keV for a thermal dark matter particle. The recent discovery of many new dark matter dominated satellites of the Milky Way in the Sloan Digital Sky Survey allows us to set lower limits comparable to constraints from the complementary methods of Lyman-alpha forest modeling and X-ray observations of the unresolved cosmic X-ray background and of dark matter halos from dwarf galaxy to cluster scales. Future surveys like LSST, DES, PanSTARRS, and SkyMapper have the potential to discover many more satellites and further improve constraints on the dark matter particle mass.

preprint2011arXiv

The Halo Model of Large Scale Structure for Warm Dark Matter

We present a comprehensive analysis of the halo model of cosmological large to small-scale structure statistics in the case of warm dark matter (WDM) structure formation scenarios. We include the effects of WDM on the linear matter power spectrum, halo density profile, halo concentration relation, halo mass function, subhalo density profile, subhalo mass function and biasing of the smooth dark matter component. As expected, we find large differences at the smallest physical scales in the nonlinear matter power spectrum predicted in the halo model between WDM and cold dark matter even for reasonably high-scale WDM particle masses. We find that significant effects are contributed from the alteration of the halo density profile and concentration, as well as the halo mass function. We further find that the effects of WDM on the subhalo population are important but sub-dominant. Clustering effects of the biasing of the smooth component in WDM is not largely significant.

preprint2011arXiv

Where are the Fossils of the First Galaxies? I. Local Volume Maps and Properties of the Undetected Dwarfs

We present a new method for generating initial conditions for LCDM N-body simulations which provides the dynamical range necessary to follow the evolution and distribution of the fossils of the first galaxies on Local Volume, 5-10 Mpc, scales. The initial distribution of particles represents the position, velocity and mass distribution of the dark and luminous halos extracted from pre-reionization simulations. We confirm previous results that ultra-faint dwarfs have properties compatible with being well preserved fossils of the first galaxies. However, because the brightest pre-reionization dwarfs form preferentially in biased regions, they most likely merge into non-fossil halos with circular velocities >20-30 km/s. Hence, we find that the maximum luminosity of true-fossils in the Milky Way is L_V<10^5 L_solar, casting doubts on the interpretation that some classical dSphs are true-fossils. In addition, we argue that most ultra-faints at small galactocentric distance, R<50 kpc, had their stellar properties modified by tides, while a large population of fossils is still undetected due to their extremely low surface brightness log(Sigma_V) < -1.4. We estimate that the region outside R_50 (~ 400 kpc) up to 1 Mpc from the Milky Way contains about a hundred true fossils of the first galaxies with V-band luminosities 10^3 - 10^5 L_solar and half-light radii, r_hl ~ 100-1000 pc.

preprint2011arXiv

Where are the Fossils of the First Galaxies? II. True Fossils, Ghost Halos, and the Missing Bright Satellites

We use a new set of cold dark matter simulations of the local universe to investigate the distribution of fossils of primordial dwarf galaxies within, and around the Milky Way. Throughout, we build upon previous results showing agreement between the observed stellar properties of a subset of the ultra-faint dwarfs and our simulated fossils. Here, we show that fossils of the first galaxies have galactocentric distributions and cumulative luminosity functions consistent with observations. In our model there are ~ 300 luminous satellites orbiting the Milky Way, ~50-70% of which are well preserved fossils, with this fraction decreasing with galactocentric distance. Within the Milky Way virial radius, the majority of these fossils have luminosities L_V<10^5 L_solar. This work produces an overabundance of bright dwarf satellites (L_V > 10^4 L_solar) with respect to observations where observations are nearly complete. The "bright satellite problem" is most evident in the outer parts of the Milky Way. We estimate that, although relatively bright, the primordial stellar populations are very diffuse, producing a population with surface brightnesses below surveys` detection limits and are easily stripped by tidal forces. Although we cannot yet present unmistakable evidence for the existence of the fossils of first galaxies in the Local Group, the results of our studies suggest observational strategies that may demonstrate their existence. Primarily, the detection of "ghost halos" of primordial stars around isolated dwarfs would prove that stars formed in minihalos (M<10^8 M_solar) before reionization, and strongly suggest that at least a fraction of the ultra-faint dwarfs are fossils of the first galaxies.

preprint2009arXiv

Pre-reionization Fossils, Ultra-faint Dwarfs and the Missing Galactic Satellite Problem

We argue that, at least a fraction of the newly discovered population of ultra-faint dwarf spheroidal galaxies in the Local Group constitute the fossil relic of a once ubiquitous population of dwarf galaxies formed before reionization with circular velocities smaller than $v_{c}^{cr} \sim 20$ km/s. We present several arguments in support of this model. The number of luminous Milky Way satellites inferred from observations is larger than the estimated number of dark halos in the Galaxy that have, or had in the past, circular velocity $>v_{c}^{cr}$, as predicted by the "Via Lactea" simulation. This implies that some ultra-faint dwarfs are fossils. However, this argument is weakened by recent results from the "Aquarius" simulations showing that the number of Galactic dark matter satellites is 2.5 larger than previously believed. Secondly, the existence of a population of ultra-faint dwarfs was predicted by cosmological simulations in which star formation in the first minihalos is reduced -- but not suppressed -- by radiative feedback. Here, we show the statistical properties of the fossil galaxies in those simulations are consistent with observations of the new dwarf population and with the number and radial distribution of Milky Way satellites as a function of their luminosity. Finally, the observed Galactocentric distribution of dwarfs is consistent with a fraction of dSphs being fossils. To make our case more compelling, future work should determine whether stellar chemical abundances of simulated "fossils" can reproduce observations and whether the tidal scenarios for the formation of Local Group dwarf spheroidals are equally consistent with all available observations.

preprint2009arXiv

The First Galaxies and the Likely Discovery of their Fossils in the Local Group

In cold dark matter cosmologies, small mass halos outnumber larger mass halos at any redshift. However, the lower bound for the mass of a galaxy is unknown, as are the typical luminosity of the smallest galaxies and their numbers in the universe. The answers depend on the extent to which star formation in the first population of small mass halos may be suppressed by radiative feedback loops operating over cosmological distance scales. If early populations of dwarf galaxies did form in significant number, their relics should be found today in the Local Group. These relics have been termed "fossils of the first galaxies". This paper is a review that summarizes our ongoing efforts to simulate and identify these fossils around the Milky Way and Andromeda. It is widely believed that reionization of the intergalactic medium would have stopped star formation in the fossils of the first galaxies. Thus, they should be among the oldest objects in the Universe. However, here we dispute this idea and discuss a physical mechanism whereby relatively recent episodes of gas accretion and star formation would be produced in some fossils of the first galaxies. We argue that fossils may be characterized either by a single old population of stars or by a bimodal star formation history. We also propose that the same mechanism could turn small mass dark halos formed before reionization into gas-rich but starless "dark galaxies". We believe that current observational data support the thesis that a fraction of the new ultra-faint dwarfs recently discovered in the Local Group are in fact fossils of the first galaxies.

preprint2007arXiv

Bondi accretion in the early universe

This paper presents a study of quasi-steady spherical accretion in the early Universe, before the formation of the first stars and galaxies. The main motivation is to derive the basic formulas that will be used in a companion paper to calculate the accretion luminosity of primordial black holes and their effect on the cosmic ionization history. The following cosmological effects are investigated: the coupling of the gas to the CMB photon fluid (i.e., Compton drag), Hubble expansion, and the growth of the dark matter halo seeded by the gravitational potential of the central point mass. The gas equations of motion are solved assuming either a polytropic or an isothermal equation of state. We consider the cases in which the accreting object is a point mass or a spherical dark matter halo with power-law density profile, as predicted by the theory of "secondary infall''. Analytical solutions for the sonic radius and fitting formulas for the accretion rate are provided. Different accretion regimes exist depending on the mass of the accreting object. If the black hole mass is smaller than 50-100 Msun, gas accretion is unaffected by Compton drag. A point mass and an extended dark halo of equal mass accrete at the same rate if M>5000 Msun, while smaller mass dark halos accrete less efficiently than the equivalent point mass. For masses M>3 x 10^4 Msun, the viscous term due to the Hubble expansion becomes important and the assumption of quasi-steady flow fails. Hence, the steady Bondi solutions transition to the time-dependent self-similar solutions for "cold cosmological infall".

preprint2006arXiv

Growth of structure seeded by primordial black holes

We discuss the possibilities for primordial black holes (PBHs) to grow via the accretion of dark matter. In agreement with previous works, we find that accretion during the radiation-dominated era does not lead to a significant mass increase. However, during matter-domination, PBHs may grow by up to two orders of magnitude in mass through the acquisition of large dark matter halos. We discuss the possibility of PBHs being an important component in dark matter halos of galaxies as well as their potential to explain the ultra-luminous x-ray sources (ULXs) observed in nearby galactic disks. We point out that although PBHs are ruled out as the dominant component of dark matter, there is still a great deal of parameter space open to them playing a role in the modern-day universe. For example, a primordial halo population of PBHs each at $10^{2.5} M_\odot$ making up 0.1% of the dark matter grow to $10^{4.5} M_\odot$ via the accumulation of dark matter halos to account for $\sim 10%$ of the dark matter mass by a redshift of $z \approx 30$. These intermediate mass black holes may then ``light up'' when passing through molecular clouds, becoming visible as ULXs at the present day, or they may form the seeds for supermassive black holes at the centers of galaxies.

preprint2002arXiv

Dwarf Sphs/First-galaxies connection

I analyze the properties of the first galaxies in cosmological simulations with radiative feedback. Preliminary results indicate similarities with the observed properties of the bulk of dwarf spheroidal galaxies (dSphs) in the Local Group and Andromeda. I briefly discuss observational tests that could help in understanding the impact of a population of small primordial objects on the cosmic evolution.

preprint1997arXiv

Energy Dissipation in Interstellar Cloud Collisions

We present a study of the kinetic energy dissipation in interstellar cloud collisions. The main aim is to understand the dependence of the elasticity (defined as the ratio of the final to the initial kinetic energy of the clouds) on the velocity and mass ratio of the colliding clouds, magnetic field strength, and gas metallicity for head-on collisions. The problem has been studied both analytically and via numerical simulations. We have derived handy analytical relationships that well approximate the analogous numerical results. The main findings of this work are: (i) the kinetic energy dissipation in cloud collisions is minimum (i.e. the collision elasticity is maximum) for a cloud relative velocity $v_r \simeq 30 km s^{-1}$; (ii) the above minimum value is proportional $Z L_c^2$, where $Z$ is the metallicity and $L_c$ is the cloud size: the larger is $Z L_c^2$ the more dissipative (i.e. inelastic) the collision will be; (iii) in general, we find that the energy dissipation decreases when the magnetic field strength, and mass ratio of the clouds are increased and the metallicity is decreased, respectively. We briefly discuss the relevance of this study to the global structure of the interstellar medium and to galaxy formation and evolution.