Showing 1–50 of 87 results for author: Ricotti, M

Explanations for the two-component spectral energy distributions of gravitationally lensed stars at high redshifts

Authors: Armin Nabizadeh, Erik Zackrisson, Emma Lundqvist, Massimo Ricotti, Seyong Park, Brian Welch, Jose M. Diego

Abstract: Observations of gravitationally lensed, high-mass stars at redshifts $\gtrsim1$ occasionally reveal spectral energy distributions that contain two components with different effective temperatures. Given that two separate stars are involved, it suggests that both stars have simultaneously reached very high magnification, as expected for two stars in a binary system close to the caustic curve of the… ▽ More Observations of gravitationally lensed, high-mass stars at redshifts $\gtrsim1$ occasionally reveal spectral energy distributions that contain two components with different effective temperatures. Given that two separate stars are involved, it suggests that both stars have simultaneously reached very high magnification, as expected for two stars in a binary system close to the caustic curve of the foreground galaxy-cluster lens. The inferred effective temperatures and luminosities of these stars are, however, difficult to reconcile with known binaries, or even with isolated stars of the same age. Here, we explore three alternative explanations for these cases: circumstellar dust around the cooler of the two stars; age differences of a few Myr among stars in the same star cluster, and a scenario in which the stars originate in two separate star clusters of different age along the lensing caustic. While all of these scenarios are deemed plausible in principle, dust solutions would require more circumstellar extinction than seen in local observations of the relevant super/hypergiant stars. Hence, we argue that age differences between the two stars are the most likely scenario, given the current data. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: 7 pages, 7 figures, 1 table, submitted to A&A

JWST NIRSpec High-resolution Spectroscopy of MACS0647-JD at z=10.167: Resolved [OII] Doublet and Electron Density in an Early Galaxy

Authors: Abdurro'uf, Rebecca L. Larson, Dan Coe, Tiger Yu-Yang Hsiao, Javier Álvarez-Márquez, Alejandro Crespo Gómez, Angela Adamo, Rachana Bhatawdekar, Arjan Bik, Larry D. Bradley, Christopher J. Conselice, Pratika Dayal, Jose M. Diego, Seiji Fujimoto, Lukas J. Furtak, Taylor A. Hutchison, Intae Jung, Meghana Killi, Vasily Kokorev, Matilde Mingozzi, Colin Norman, Tom Resseguier, Massimo Ricotti, Jane R. Rigby, Eros Vanzella , et al. (4 additional authors not shown)

Abstract: We present JWST/NIRSpec high-resolution spectroscopy G395H/F290LP of MACS0647-JD, a gravitationally lensed galaxy merger at $z=10.167$. The new spectroscopy, which is acquired for the two lensed images (JD1 and JD2), detects and resolves emission lines in the rest-frame ultraviolet (UV) and blue optical, including the resolved [OII]3726,3729 doublet, [NeIII]3870, [HeI]3890, H$δでるた$, H$γがんま$, and [OIII]4… ▽ More We present JWST/NIRSpec high-resolution spectroscopy G395H/F290LP of MACS0647-JD, a gravitationally lensed galaxy merger at $z=10.167$. The new spectroscopy, which is acquired for the two lensed images (JD1 and JD2), detects and resolves emission lines in the rest-frame ultraviolet (UV) and blue optical, including the resolved [OII]3726,3729 doublet, [NeIII]3870, [HeI]3890, H$δでるた$, H$γがんま$, and [OIII]4363. This is the first observation of the resolved [OII]3726,3729 doublet for a galaxy at $z>8$. We measure a line flux ratio [OII]3729/3726 $= 0.9 \pm 0.3$, which corresponds to an estimated electron density of $\log(n_{e} / \rm{cm}^{-3}) = 2.9 \pm 0.5$. This is significantly higher than the electron densities of local galaxies reported in the literature. We compile the measurements from the literature and further analyze the redshift evolution of $n_{e}$. We find that the redshift evolution follows the power-law form of $n_{e} = A\times (1+z)^{p}$ with $A=54^{+31}_{-23}$ cm$^{-3}$ and $p=1.2^{+0.4}_{-0.4}$. This power-law form may be explained by a combination of metallicity and morphological evolution of galaxies, which become, on average, more metal-poor and more compact with increasing redshift. △ Less

Submitted 4 July, 2024; v1 submitted 24 April, 2024; originally announced April 2024.

Comments: 20 pages, 11 figures. Accepted for publication in ApJ

JWST MIRI detections of H$αあるふぁ$ and [O III] and direct metallicity measurement of the $z=10.17$ lensed galaxy MACS0647$-$JD

Authors: Tiger Yu-Yang Hsiao, Javier Álvarez-Márquez, Dan Coe, Alejandro Crespo Gómez, Abdurro'uf, Pratika Dayal, Rebecca L. Larson, Arjan Bik, Carmen Blanco-Prieto, Luis Colina, Pablo Guillermo Pérez-González, Luca Costantin, Carlota Prieto-Jiménez, Angela Adamo, Larry D. Bradley, Christopher J. Conselice, Seiji Fujimoto, Lukas J. Furtak, Taylor A. Hutchison, Bethan L. James, Yolanda Jiménez-Teja, Intae Jung, Vasily Kokorev, Matilde Mingozzi, Colin Norman , et al. (8 additional authors not shown)

Abstract: JWST spectroscopy has revolutionized our understanding of galaxies in the early universe. Covering wavelengths up to $5.3\,{\rm μみゅーm}$, NIRSpec can detect rest-frame optical emission lines H$αあるふぁ$ out to $z = 7$ and [O III] to $z = 9.5$. Observing these lines in more distant galaxies requires longer wavelength spectroscopy with MIRI. Here we present MIRI MRS IFU observations of the lensed galaxy merger… ▽ More JWST spectroscopy has revolutionized our understanding of galaxies in the early universe. Covering wavelengths up to $5.3\,{\rm μみゅーm}$, NIRSpec can detect rest-frame optical emission lines H$αあるふぁ$ out to $z = 7$ and [O III] to $z = 9.5$. Observing these lines in more distant galaxies requires longer wavelength spectroscopy with MIRI. Here we present MIRI MRS IFU observations of the lensed galaxy merger MACS0647$-$JD at $z = 10.165$. With exposure times of 4.2 hours in each of two bands, we detect H$αあるふぁ$ at $9σしぐま$, [O III]$\,\lambda5008$ at $11σしぐま$, and [O III]$\,\lambda4960$ at $3σしぐま$. Combined with previously reported NIRSpec spectroscopy that yields seven emission lines including the auroral line [O III]$\,\lambda4363$, we present the first direct metallicity measurement of a $z > 10$ galaxy: $12+{\rm log(O/H)}= 7.79\pm0.09$, or $0.13^{+0.02}_{-0.03}\,Z_{\odot}$. This is similar to galaxies at $z \sim 4 - 9$ with direct metallicity measurements, though higher than expected given the high specific star formation rate ${\rm log(sSFR / yr^{-1})} = -7.4 \pm 0.3$. We further constrain the ionization parameter ${\rm log}(U)$ = $-1.9 \pm 0.1$, ionizing photon production efficiency ${\rm log}(ξくしー_{\rm ion})$ = $25.3\pm0.1$, and star formation rate $5.0\pm0.6\,M_{\odot}/{\rm yr}$ within the past $10\,{\rm Myr}$. These observations demonstrate the combined power of JWST NIRSpec and MIRI for studying galaxies in the first $500$ million years. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: 15 pages, 8 figures, submitted to ApJ

Unveiling the Cosmic Gems Arc at $z\sim10.2$ with JWST

Authors: Larry D. Bradley, Angela Adamo, Eros Vanzella, Keren Sharon, Gabriel Brammer, Dan Coe, Jose M. Diego, Vasily Kokorev, Guillaume Mahler, Masamune Oguri, Abdurro'uf, Rachana Bhatawdekar, Lise Christensen, Seiji Fujimoto, Takuya Hashimoto, Tiger Y. -Y Hsiao, Akio K. Inoue, Yolanda Jiménez-Teja, Matteo Messa, Colin Norman, Massimo Ricotti, Yoichi Tamura, Rogier A. Windhorst, Xinfeng Xu, Adi Zitrin

Abstract: We present recent JWST NIRCam imaging observations of SPT0615-JD (also known as the Cosmic Gems Arc), lensed by the galaxy cluster SPT-CL J0615-5746. The 5-arcsec-long arc is the most highly magnified $z>10$ galaxy known, straddling the lensing critical curve and revealing five star clusters with radii $\sim 1$ pc or less. We measure the full arc to have F200W 24.5 AB mag, consisting of two mirror… ▽ More We present recent JWST NIRCam imaging observations of SPT0615-JD (also known as the Cosmic Gems Arc), lensed by the galaxy cluster SPT-CL J0615-5746. The 5-arcsec-long arc is the most highly magnified $z>10$ galaxy known, straddling the lensing critical curve and revealing five star clusters with radii $\sim 1$ pc or less. We measure the full arc to have F200W 24.5 AB mag, consisting of two mirror images, each 25.3 AB mag with a magnification $μみゅー\sim 60$ (delensed 29.7 AB mag, $M_{UV} = -17.8$). The galaxy has an extremely strong Lyman break F115W$-$F200W $>3.2$ mag ($2σしぐま$ lower limit), is undetected in all bluer filters ($< 2σしぐま$), and has a very blue continuum slope redward of the break ($βべーた= -2.7 \pm 0.1$), resulting in a photometric redshift $z_{phot} = 10.2 \pm 0.2$ (95% confidence) with no significant likelihood below $z < 9.8$. Based on SED fitting to the total photometry, we estimate an intrinsic stellar mass of $M_{*} \sim 2.4 - 5.6 \times 10^{7} M_{\odot}$, young mass-weighted age of $\sim 21 - 79$ Myr, low dust content ($A_V < 0.15$), and a low metallicity of $\lesssim 1\%~Z_{\odot}$. We identify a fainter third counterimage candidate within 2.2 arcsec of the predicted position, lensed to AB mag 28.4 and magnified by $μみゅー\sim 2$, suggesting the fold arc may only show $\sim60$% of the galaxy. SPT0615-JD is a unique laboratory to study star clusters observed within a galaxy just 460 Myr after the Big Bang. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: 22 pages, 8 figures, 4 tables, submitted to ApJ

Formation of Large Circumstellar Discs in Multi-scale, ideal-MHD Simulations of Magnetically Critical Pre-stellar Cores

Authors: Chong-Chong He, Massimo Ricotti

Abstract: The formation of circumstellar discs is a critical step in the formation of stars and planets. Magnetic fields can strongly affect the evolution of angular momentum during prestellar core collapse, potentially leading to the failure of protostellar disc formation. This phenomenon, known as the magnetic braking catastrophe, has been observed in ideal-MHD simulations. In this work, we present result… ▽ More The formation of circumstellar discs is a critical step in the formation of stars and planets. Magnetic fields can strongly affect the evolution of angular momentum during prestellar core collapse, potentially leading to the failure of protostellar disc formation. This phenomenon, known as the magnetic braking catastrophe, has been observed in ideal-MHD simulations. In this work, we present results from ideal-MHD simulations of circumstellar disc formation from realistic initial conditions of strongly magnetised, massive cores with masses between $30 ~{\rm M}_\odot$ and $300 ~{\rm M}_\odot$ resolved by zooming into Giant Molecular Clouds with masses $\sim 10^4 \ {\rm M}_\odot$ and initial mass-to-flux ratios $0.6 \le μみゅー_0 \le 3$. Due to the large turbulence caused by the non-axisymmetric gravitational collapse of the gas, the dominant vertical support of discs is turbulent motion, while magnetic and turbulent pressures contribute equally in the outer toroid. The magnetic field topology is extremely turbulent and incoherent, reducing the effect of magnetic braking by roughly one order of magnitude and leading to the formation of large Keplerian discs even in magnetically critical or near-critical cores. Only cores in GMCs with $μみゅー_0 < 1$ fail to form discs. Instead, they collapse into a sheet-like structure and produce numerous low-mass stars. We also discuss a universal $B-ρろー$ relation valid over a large range of scales from the GMC to massive cores, irrespective of the GMC magnetisation. This study differs from the vast literature on this topic which typically focus on smaller mass discs with idealised initial and boundary conditions, therefore providing insights into the initial conditions of massive prestellar core collapse and disc formation. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: 15 pages, 11 figures, submitted to MNRAS

Violent starbursts and quiescence induced by FUV radiation feedback in metal-poor galaxies at high-redshift

Authors: Kazuyuki Sugimura, Massimo Ricotti, Jongwon Park, Fred Angelo Batan Garcia, Hidenobu Yajima

Abstract: JWST observations of galaxies at $z\gtrsim 8$ suggest that they are more luminous and clumpier than predicted by most models, prompting several proposals on the physics of star formation and feedback in the first galaxies. In this paper, we focus on the role of ultraviolet (UV) radiation in regulating star formation by performing a set of cosmological radiation hydrodynamics simulations of one gal… ▽ More JWST observations of galaxies at $z\gtrsim 8$ suggest that they are more luminous and clumpier than predicted by most models, prompting several proposals on the physics of star formation and feedback in the first galaxies. In this paper, we focus on the role of ultraviolet (UV) radiation in regulating star formation by performing a set of cosmological radiation hydrodynamics simulations of one galaxy at sub-pc resolution with different radiative feedback models. We find that the suppression of cooling by far UV (FUV) radiation (i.e., $\mathrm{H_2}$ dissociating radiation) from Pop II stars is the main physical process triggering the formation of compact and massive star clusters and is responsible for the bursty star formation observed in metal-poor galaxies at $z\gtrsim 10$. Indeed, artificially suppressing FUV radiation leads to a less intense continuous mode of star formation distributed into numerous, but low-mass open star clusters. Due to the intense FUV field, low-metallicity clouds remain warm ($\sim 10^4\,\mathrm{K}$) until they reach a relatively high density ($\gtrsim 10^3\,\mathrm{cm^{-3}}$), before becoming self-shielded and transitioning to a colder ($\sim 100\,\mathrm{K}$), partially molecular phase. As a result, star formation is delayed until the clouds accumulate enough mass to become gravitationally unstable. At this point, the clouds undergo rapid star formation converting gas into stars with high efficiency. We, therefore, observe exceptionally bright galaxies (ten times brighter than for continuous star formation) and subsequent quenched "dead" galaxies that did not form stars for tens of Myrs. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: 23 pages, 18 figures, submitted to AAS journal

Bound star clusters observed in a lensed galaxy 460 Myr after the Big Bang

Authors: Angela Adamo, Larry D. Bradley, Eros Vanzella, Adélaïde Claeyssens, Brian Welch, Jose M Diego, Guillaume Mahler, Masamune Oguri, Keren Sharon, Abdurro'uf, Tiger Yu-Yang Hsiao, Xinfeng Xu, Matteo Messa, Augusto E. Lassen, Erik Zackrisson, Gabriel Brammer, Dan Coe, Vasily Kokorev, Massimo Ricotti, Adi Zitrin, Seiji Fujimoto, Akio K. Inoue, Tom Resseguier, Jane R. Rigby, Yolanda Jiménez-Teja , et al. (3 additional authors not shown)

Abstract: The Cosmic Gems arc is among the brightest and highly magnified galaxies observed at redshift $z\sim10.2$. However, it is an intrinsically UV faint galaxy, in the range of those now thought to drive the reionization of the Universe. Hitherto the smallest features resolved in a galaxy at a comparable redshift are between a few hundreds and a few tens of parsecs. Here we report JWST observations of… ▽ More The Cosmic Gems arc is among the brightest and highly magnified galaxies observed at redshift $z\sim10.2$. However, it is an intrinsically UV faint galaxy, in the range of those now thought to drive the reionization of the Universe. Hitherto the smallest features resolved in a galaxy at a comparable redshift are between a few hundreds and a few tens of parsecs. Here we report JWST observations of the Cosmic Gems. The light of the galaxy is resolved into five star clusters located in a region smaller than 70 parsec. They exhibit minimal dust attenuation and low metallicity, ages younger than 50 Myr and intrinsic masses of $\sim10^6$ M$_{\odot}$. Their lensing-corrected sizes are approximately 1 pc, resulting in stellar surface densities near $10^5$~M$_{\odot}$/pc$^2$, three orders of magnitude higher than typical young star clusters in the local universe. Despite the uncertainties inherent to the lensing model, they are consistent with being gravitationally bound stellar systems, i.e., proto-globular clusters. We conclude that star cluster formation and feedback likely contributed to shape the properties of galaxies during the epoch of reionization. [Abridged] △ Less

Submitted 12 June, 2024; v1 submitted 6 January, 2024; originally announced January 2024.

Comments: Accepted for publication

To be, or not to be: Balmer breaks in high-z galaxies with JWST

Authors: Anton Vikaeus, Erik Zackrisson, Stephen Wilkins, Armin Nabizadeh, Vasily Kokorev, Abdurrouf, Larry D. Bradley, Dan Coe, Pratika Dayal, Massimo Ricotti

Abstract: Standard models of structure formation allow us to predict the cosmic timescales relevant for the onset of star formation and the assembly history of galaxies at high redshifts ($z > 10$). The strength of the Balmer break represents a well-known diagnostic of the age and star formation history of galaxies, which enables us to compare observations with contemporary simulations - thus shedding light… ▽ More Standard models of structure formation allow us to predict the cosmic timescales relevant for the onset of star formation and the assembly history of galaxies at high redshifts ($z > 10$). The strength of the Balmer break represents a well-known diagnostic of the age and star formation history of galaxies, which enables us to compare observations with contemporary simulations - thus shedding light on the predictive power of our current models of star formation in the early universe. Here, we measure the Balmer break strength for 23 spectroscopically confirmed galaxies at redshifts $6 \lesssim z \lesssim 12$ using public JWST NIRSpec data from the cycle 1 GO 1433 and GO 2282 programs (PI Coe), as well as public spectroscopic data from the JWST Deep Extragalactic Survey (JADES). We find that the range of observed Balmer break strengths agree well with that of current simulations given our measurement uncertainties. No cases of anomalously strong Balmer breaks are detected, and therefore no severe departures from the predictions of contemporary models of star formation. However, there are indications that the number of outliers in the observed distribution, both in direction of strong and weak Balmer breaks, is higher than that predicted by simulations. △ Less

Submitted 7 September, 2023; v1 submitted 5 September, 2023; originally announced September 2023.

Comments: 9 pages, 9 figures, 1 table. v2: Edited metadata and format errors

Reaching for the stars -- JWST/NIRSpec spectroscopy of a lensed star candidate at $z=4.76$

Authors: Lukas J. Furtak, Ashish K. Meena, Erik Zackrisson, Adi Zitrin, Gabriel B. Brammer, Dan Coe, José M. Diego, Jan J. Eldridge, Yolanda Jiménez-Teja, Vasily Kokorev, Massimo Ricotti, Brian Welch, Rogier A. Windhorst, Abdurro'uf, Felipe Andrade-Santos, Rachana Bhatawdekar, Larry D. Bradley, Tom Broadhurst, Wenlei Chen, Christopher J. Conselice, Pratika Dayal, Brenda L. Frye, Seiji Fujimoto, Tiger Y. -Y. Hsiao, Patrick L. Kelly , et al. (9 additional authors not shown)

Abstract: We present JWST/NIRSpec observations of a highly magnified star candidate at a photometric redshift of $z_{\mathrm{phot}}\simeq4.8$, previously detected in JWST/NIRCam imaging of the strong lensing (SL) cluster MACS J0647+7015 ($z=0.591$). The spectroscopic observation allows us to precisely measure the redshift of the host arc at $z_{\mathrm{spec}}=4.758\pm0.004$, and the star's spectrum displays… ▽ More We present JWST/NIRSpec observations of a highly magnified star candidate at a photometric redshift of $z_{\mathrm{phot}}\simeq4.8$, previously detected in JWST/NIRCam imaging of the strong lensing (SL) cluster MACS J0647+7015 ($z=0.591$). The spectroscopic observation allows us to precisely measure the redshift of the host arc at $z_{\mathrm{spec}}=4.758\pm0.004$, and the star's spectrum displays clear Lyman- and Balmer-breaks commensurate with this redshift. A fit to the spectrum suggests a B-type super-giant star of surface temperature $T_{\mathrm{eff,B}}\simeq15000$ K with either a redder F-type companion ($T_{\mathrm{eff,F}}\simeq6250$K) or significant dust attenuation ($A_V\simeq0.82$) along the line of sight. We also investigate the possibility that this object is a magnified young globular cluster rather than a single star. We show that the spectrum is in principle consistent with a star cluster, which could also accommodate the lack of flux variability between the two epochs. However, the lack of a counter image and the strong upper limit on the size of the object from lensing symmetry, $r\lesssim0.5$ pc, could indicate that this scenario is somewhat less likely -- albeit not completely ruled out by the current data. The presented spectrum seen at a time when the Universe was only $\sim1.2$ Gyr old showcases the ability of JWST to study early stars through extreme lensing. △ Less

Submitted 25 September, 2023; v1 submitted 31 July, 2023; originally announced August 2023.

Comments: Accepted for publication in MNRAS letters. v2 updated to match the published version

On the origin of outward migration of Population III stars

Authors: Jongwon Park, Massimo Ricotti, Kazuyuki Sugimura

Abstract: Outward migration of massive binary stars or black holes in their circumbinary disc is often observed in simulations and it is key to the formation of wide black hole binaries. Using numerical simulations of Population III (Pop III) star formation, we study the angular momentum of Pop III binaries and the torques between stars and gas discs to understand the origin of outward migration and high el… ▽ More Outward migration of massive binary stars or black holes in their circumbinary disc is often observed in simulations and it is key to the formation of wide black hole binaries. Using numerical simulations of Population III (Pop III) star formation, we study the angular momentum of Pop III binaries and the torques between stars and gas discs to understand the origin of outward migration and high ellipticity. The outward migration of protostars is produced by gravitational torques exerted on them by their circumstellar minidiscs. The minidiscs, on the other hand, migrate outward mainly by gaining angular momentum by accreting gas from the circumbinary disc. The angular momentum transfer is most efficient for rapidly accreting equal-mass binaries, and weaker when the secondary mass is small or the massive companion evaporates the gas disc via radiative feedback. We conclude that outward migration and the formation of wide equal-mass massive binaries is common in metal-free/metal-poor star formation, mainly driven by their large accretion rates. We expect that the lower gas temperature and accretion rates in metal-enriched circumstellar discs would lead more often to inward migration and closer binary separations. We also observe inward migration for smaller mass Pop III protostars/fragments, leading to the rapid merging of sink particles and likely the formation of close binary black holes that, however, reach separations below the resolution of our simulations. We discuss the implications that Pop III separations and ellipticity may have on the interpretation that gravitational wave signals from merging intermediate-mass black holes come from Pop III remnants. △ Less

Submitted 15 February, 2024; v1 submitted 26 July, 2023; originally announced July 2023.

Comments: 20 pages, 17+4 figures. Accepted for publication in MNRAS

Optimizing Roman's High Latitude Wide Area Survey for Low Surface Brightness Astronomy

Authors: Mireia Montes, Francesca Annibali, Michele Bellazzini, Alejandro S. Borlaff, Sarah Brough, Fernando Buitrago, Nushkia Chamba, Chris Collins, Ian Dell'Antonio, Ivanna Escala, Anthony H. Gonzalez, Benne Holwerda, Sugata Kaviraj, Johan Knapen, Anton Koekemoer, Seppo Laine, Pamela Marcum, Garreth Martin, David Martinez-Delgado, Chris Mihos, Massimo Ricotti, Ignacio Trujillo, Aaron E. Watkins

Abstract: One of the last remaining frontiers in optical/near-infrared observational astronomy is the low surface brightness regime (LSB, V-band surface brightness, $μみゅー_V>$ 27 AB mag/arcsec$^2$). These are the structures at very low stellar surface densities, largely unseen by even current wide-field surveys such as the Legacy Survey. Studying this domain promises to be transformative for our understanding o… ▽ More One of the last remaining frontiers in optical/near-infrared observational astronomy is the low surface brightness regime (LSB, V-band surface brightness, $μみゅー_V>$ 27 AB mag/arcsec$^2$). These are the structures at very low stellar surface densities, largely unseen by even current wide-field surveys such as the Legacy Survey. Studying this domain promises to be transformative for our understanding of star formation in low-mass galaxies, the hierarchical assembly of galaxies and galaxy clusters, and the nature of dark matter. It is thus essential to reach depths beyond $μみゅー_V$ = 30 AB mag/arcsec$^2$ to detect the faintest extragalactic sources, such as dwarf galaxies and the stellar halos around galaxies and within galaxy clusters. The High Latitude Wide Area Survey offers a unique opportunity to statistically study the LSB universe at unprecedented depths in the IR over an area of $\sim$2000 square degrees. The high spatial resolution will minimize source confusion, allowing an unbiased characterization of LSB structures, including the identification of stars in nearby LSB galaxies and globular clusters. In addition, the combination of Roman with other upcoming deep imaging observatories (such as Rubin) will provide multi-wavelength coverage to derive photometric redshifts and infer the stellar populations of LSB objects. △ Less

Submitted 15 June, 2023; originally announced June 2023.

Comments: White paper submitted to the call for input for the Roman Space Telescope's Core Community Surveys

An extremely metal poor star complex in the reionization era: Approaching Population III stars with JWST

Authors: E. Vanzella, F. Loiacono, P. Bergamini, U. Mestric, M. Castellano, P. Rosati, M. Meneghetti, C. Grillo, F. Calura, M. Mignoli, M. Bradac, A. Adamo, G. Rihtarsic, M. Dickinson, M. Gronke, A. Zanella, F. Annibali, C. Willott, M. Messa, E. Sani, A. Acebron, A. Bolamperti, A. Comastri, R. Gilli, K. I. Caputi , et al. (9 additional authors not shown)

Abstract: We present JWST/NIRSpec integral field spectroscopy (IFS) of a lensed Population III candidate stellar complex (dubbed Lensed And Pristine 1, LAP1), with a lensing-corrected stellar mass ~<10^4 Msun, absolute luminosity M_UV > -11.2 (m_UV > 35.6), confirmed at redshift 6.639 +/- 0.004. The system is strongly amplified (μみゅー>~ 100) by straddling a critical line of the Hubble Frontier Field galaxy clus… ▽ More We present JWST/NIRSpec integral field spectroscopy (IFS) of a lensed Population III candidate stellar complex (dubbed Lensed And Pristine 1, LAP1), with a lensing-corrected stellar mass ~<10^4 Msun, absolute luminosity M_UV > -11.2 (m_UV > 35.6), confirmed at redshift 6.639 +/- 0.004. The system is strongly amplified (μみゅー>~ 100) by straddling a critical line of the Hubble Frontier Field galaxy cluster MACS J0416. Despite the stellar continuum is currently not detected in the Hubble and JWST/NIRCam and NIRISS imaging, arclet-like shapes of Lyman and Balmer lines, Lya, Hg, Hb and Ha are detected with NIRSpec IFS with signal-to-noise ratios SNR=5-13 and large equivalent widths (>300-2000A), along with a remarkably weak [OIII]4959-5007 at SNR ~ 4. LAP1 shows a large ionizing photon production efficiency, log(ξくしー_{ion}[erg~Hz^{-1}])>26. From the metallicity indexes R23 = ([OIII]4959-5007 + [OII]3727) / Hb ~< 0.74 and R3 = ([OIII]5007 / Hb) = 0.55 +/- 0.14, we derive an oxygen abundance 12+log(O/H) ~< 6.3. Intriguingly, the Ha emission is also measured in mirrored sub-components where no [OIII] is detected, providing even more stringent upper limits on the metallicity if in-situ star formation is ongoing in this region (12+log(O/H) < 6, or Z < 0.002 Zsun). The formal stellar mass limit of the sub-components would correspond to ~10^{3} Msun or M_UV fainter than -10. Alternatively, such a metal-free pure line emitting region could be the first case of a fluorescing HI gas region, induced by transverse escaping ionizing radiation from a nearby star-complex. The presence of large equivalent-width hydrogen lines and the deficiency of metal lines in such a small region, make LAP1 the most metal poor star-forming region currently known in the reionization era and a promising site that may host isolated, pristine stars. △ Less

Submitted 23 May, 2023; originally announced May 2023.

Comments: 15 pages, 10 figures, 1 table. Submitted to A&A. Comments welcome

Journal ref: A&A 678, A173 (2023)

JWST NIRSpec spectroscopy of the triply-lensed $z = 10.17$ galaxy MACS0647$-$JD

Authors: Tiger Yu-Yang Hsiao, Abdurro'uf, Dan Coe, Rebecca L. Larson, Intae Jung, Matilde Mingozzi, Pratika Dayal, Nimisha Kumari, Vasily Kokorev, Anton Vikaeus, Gabriel Brammer, Lukas J. Furtak, Angela Adamo, Felipe Andrade-Santos, Jacqueline Antwi-Danso, Marusa Bradac, Larry D. Bradley, Tom Broadhurst, Adam C. Carnall, Christopher J. Conselice, Jose M. Diego, Megan Donahue, Jan J. Eldridge, Seiji Fujimoto, Alaina Henry , et al. (16 additional authors not shown)

Abstract: We present JWST/NIRSpec prism spectroscopy of MACS0647-JD, the triply-lensed $z \sim 11$ candidate discovered in HST imaging and spatially resolved by JWST imaging into two components A and B. Spectroscopy of component A yields a spectroscopic redshift $z=10.17$ based on 7 detected emission lines: CIII] $λらむだλらむだ$1907,1909, [OII] $λらむだ$3727, [NeIII] $λらむだ$3869, [NeIII] $λらむだ$3968, H$δでるた$ $λらむだ$4101, H$γがんま$ $λらむだ$4340, and… ▽ More We present JWST/NIRSpec prism spectroscopy of MACS0647-JD, the triply-lensed $z \sim 11$ candidate discovered in HST imaging and spatially resolved by JWST imaging into two components A and B. Spectroscopy of component A yields a spectroscopic redshift $z=10.17$ based on 7 detected emission lines: CIII] $λらむだλらむだ$1907,1909, [OII] $λらむだ$3727, [NeIII] $λらむだ$3869, [NeIII] $λらむだ$3968, H$δでるた$ $λらむだ$4101, H$γがんま$ $λらむだ$4340, and [OIII] $λらむだ$4363. These are the second-most distant detections of these emission lines to date, in a galaxy observed just 460 million years after the Big Bang. Based on observed and extrapolated line flux ratios we derive a gas-phase metallicity $Z =$ log(O/H) = $7.5 - 8.0$, or $(0.06 - 0.2)$ $Z_\odot$, ionization parameter log($U$) $\sim -1.9\pm0.2$, and an ionizing photon production efficiency ${\rm log}(ξくしー_{\rm ion})=25.2\pm0.2\,$erg$^{-1}$ Hz. The spectrum has a softened Lyman-$αあるふぁ$ break, evidence for a strong Ly$αあるふぁ$ damping wing, suggesting that MACS0647-JD was unable to ionize its surroundings beyond its immediate vicinity ($R_{\text{HII}} \ll 1$ pMpc). The Ly$αあるふぁ$ damping wing also suppresses the F150W photometry, explaining the slightly overestimated photometric redshift $z = 10.6 \pm 0.3$. MACS0647-JD has a stellar mass log($M/M_\odot$) = $8.1 \pm 0.3$, including $\sim$ 6$\times 10^7 M_\odot$ in component A, most of which formed recently (within $\sim$ 20 Myr) with a star formation rate $2\pm1 M_\odot$ / yr, all within an effective radius $70\pm24\,$pc. The smaller component B ($r \sim 20$) pc is likely older ($\sim$100 Myr) with more dust ($A_V \sim 0.1$ mag), as found previously. Spectroscopy of a fainter companion galaxy C separated by a distance of \about\ 3$\,$kpc reveals a Lyman break consistent with $z = 10.17$. MACS0647-JD is likely the most distant galaxy merger known. △ Less

Submitted 2 July, 2024; v1 submitted 4 May, 2023; originally announced May 2023.

Comments: 22 pages, 6 figures, accepted for publication in ApJ

The Extended [CII] under Construction? Observation of the brightest high-z lensed star-forming galaxy at z = 6.2

Authors: Yoshinobu Fudamoto, Akio K. Inoue, Dan Coe, Brian Welch, Ana Acebron, Massimo Ricotti, Nir Mandelker, Rogier A. Windhorst, Xinfeng Xu, Yuma Sugahara, Franz E. Bauer, Maruša Bradač, Larry D. Bradley, Jose M. Diego, Michael Florian, Brenda Frye, Seiji Fujimoto, Takuya Hashimoto, Alaina Henry, Guillaume Mahler, Pascal A. Oesch, Swara Ravindranath, Jane Rigby, Victoria Strait, Yoichi Tamura , et al. (4 additional authors not shown)

Abstract: We present results of [CII]$\,158\,\rm{μみゅーm}$ emission line observations, and report the spectroscopic redshift confirmation of a strongly lensed ($μみゅー\sim20$) star-forming galaxy, MACS0308-zD1 at $z=6.2078\pm0.0002$. The [CII] emission line is detected with a signal-to-noise ratio $>6$ within the rest-frame UV bright clump of the lensed galaxy (zD1.1) and exhibits multiple velocity components; the na… ▽ More We present results of [CII]$\,158\,\rm{μみゅーm}$ emission line observations, and report the spectroscopic redshift confirmation of a strongly lensed ($μみゅー\sim20$) star-forming galaxy, MACS0308-zD1 at $z=6.2078\pm0.0002$. The [CII] emission line is detected with a signal-to-noise ratio $>6$ within the rest-frame UV bright clump of the lensed galaxy (zD1.1) and exhibits multiple velocity components; the narrow [CII] has a velocity full-width-half-maximum (FWHM) of $110\pm20\,\rm{km/s}$, while broader [CII] is seen with an FWHM of $230\pm20\,\rm{km/s}$. The broader [CII] component is blueshifted ($-80\pm20\,\rm{km/s}$) with respect to the narrow [CII] component, and has a morphology which extends beyond the UV-bright clump. We find that while the narrow [CII] emission is most likely associated with zD1.1, the broader component is possibly associated with outflowing gas. Based on the non-detection of $λらむだ_{\rm 158\,μみゅーm}$ dust continuum, we find that MACS0308-zD1's star-formation activity occurs in a dust-free environment with the stringent upper limit of infrared luminosity $\lesssim9\times10^{8}\,{\rm L_{\odot}}$. Targeting this strongly lensed faint galaxy for follow-up ALMA and JWST observations will be crucial to characterize the details of typical galaxy growth in the early Universe. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: 10 pages, 6 figures, 1 table; submitted to ApJ

Star Cluster Formation and Survival in the First Galaxies

Authors: Fred Angelo Batan Garcia, Massimo Ricotti, Kazuyuki Sugimura, Jongwon Park

Abstract: Using radiation-hydrodynamic cosmological simulations, we present a detailed ($0.1$ pc resolution), physically motivated portrait of a typical-mass dwarf galaxy before the epoch of reionization, resolving the formation and evolution of star clusters into individual $10\:\mathrm{M_{\odot}}$ star particles. In the rest-frame UV, the galaxy has an irregular morphology with no bulge or galactic disk,… ▽ More Using radiation-hydrodynamic cosmological simulations, we present a detailed ($0.1$ pc resolution), physically motivated portrait of a typical-mass dwarf galaxy before the epoch of reionization, resolving the formation and evolution of star clusters into individual $10\:\mathrm{M_{\odot}}$ star particles. In the rest-frame UV, the galaxy has an irregular morphology with no bulge or galactic disk, dominated by light emitted from numerous, compact, and gravitationally-bound star clusters. This is especially interesting in light of recent HST and JWST observations that -- aided by the magnifying power of gravitational lenses -- have imaged, at parsec-scale resolution, individual young star clusters in the process of forming in similar galaxies at $z>6$. Because of their low metallicities and high temperatures, star-forming gas clouds in this galaxy have densities $\sim 100$ times higher than typical giant molecular clouds; hence, their star formation efficiencies are high enough ($f_*\sim10-70$ per cent) to produce a sizeable population of potential globular cluster progenitors but typically smaller (between a few $100\:-\: 2\times10^4\:\mathrm{M_{\odot}}$, sizes of $0.1-3$ pc) and of lower metallicities ($10^{-3.5}-10^{-2.5}\:\mathrm{Z_{\odot}}$). The initial mass function of the star-forming clouds is log-normal while the bound star cluster mass function is a power-law with a slope that depends mainly on $f_*$ but also on the temporal proximity to a major starburst. We find slopes between $-0.5$ and $-2.5$ depending on the assumed sub-grid $f_*$. Star formation is self-regulated on galactic scales; however, the multi-modal metallicity distribution of the star clusters and the fraction of stars locked into surviving bound star clusters depends on $f_*$. △ Less

Submitted 28 December, 2022; originally announced December 2022.

Comments: 20 pages, 15 figures, submitted to MNRAS. See https://fred144.github.io/research.html for movie renders

Population III star formation in an X-ray background: III. Periodic radiative feedback and luminosity induced by elliptical orbits

Authors: Jongwon Park, Massimo Ricotti, Kazuyuki Sugimura

Abstract: We model Pop III star formation in different FUV and X-ray backgrounds, including radiation feedback from protostars. We confirm previous results that a moderate X-ray background increases the number of Pop III systems per unit cosmological volume, but masses and multiplicities of the system are reduced. The stellar mass function also agrees with previous results, and we confirm the outward migrat… ▽ More We model Pop III star formation in different FUV and X-ray backgrounds, including radiation feedback from protostars. We confirm previous results that a moderate X-ray background increases the number of Pop III systems per unit cosmological volume, but masses and multiplicities of the system are reduced. The stellar mass function also agrees with previous results, and we confirm the outward migration of the stars within the protostellar discs. We find that nearly all Pop III star systems are hierarchical, i.e., binaries of binaries. Typically, two equal-mass stars form near the centre of the protostellar disc and migrate outward. Around these stars, mini-discs fragment forming binaries that also migrate outward. Stars may also form at Lagrange points L4/L5 of the system. Afterward, star formation becomes more stochastic due to the large multiplicity, and zero-metallicity low-mass stars can form when rapidly ejected from the disc. Stars in the disc often have eccentric orbits, leading to a periodic modulation of their accretion rates and luminosities. At the pericenter, due to strong accretion, the star can enter a red-supergiant phase reaching nearly Eddington luminosity in the optical bands ($m_{\rm AB} \sim 34$ for a $100~M_{odot}$ star at $z=6$). During this phase, the star, rather than its nebular lines, can be observed directly by JWST, if sufficiently magnified by a gravitational lens. The $\sim 10,000$ AUえーゆー separations and high eccentricities of many Pop III star binaries in our simulations are favorable parameters for IMBH mergers - and gravitational waves emission - through orbital excitation by field stars. △ Less

Submitted 8 December, 2022; originally announced December 2022.

Comments: Submitted to MNRAS. Movies are available at https://jwpark5064.github.io/ or https://www.youtube.com/@jp_astro_simulation

Bimodal Star Formation in Simulations of Strongly Magnetized Giant Molecular Clouds

Authors: Ronan Hix, Chong-Chong He, Massimo Ricotti

Abstract: We present the results of a set of radiation magnetohydrodynamic simulations of turbulent molecular clouds in which we vary the initial strength of the magnetic field within a range ($1 \lesssim μみゅー\lesssim 5$) consistent with observations of local giant molecular clouds (GMCs). We find that as we increase the strength of the magnetic field, star formation transitions from unimodal (the baseline cas… ▽ More We present the results of a set of radiation magnetohydrodynamic simulations of turbulent molecular clouds in which we vary the initial strength of the magnetic field within a range ($1 \lesssim μみゅー\lesssim 5$) consistent with observations of local giant molecular clouds (GMCs). We find that as we increase the strength of the magnetic field, star formation transitions from unimodal (the baseline case, $μみゅー=5$, with a single burst of star formation and Salpeter IMF) to bimodal. This effect is clearest in the most strongly magnetized GMCs ($μみゅー=1$): a first burst of star formation with duration, intensity and IMF comparable to the baseline case, is followed by a second star formation episode in which only low-mass stars are formed. Overall, due to the second burst of star formation, the strongly magnetized case results in a longer star formation period and a higher efficiency of star formation. The second burst is produced by gas that is not expelled by radiative feedback, instead remaining trapped in the GMC by the large-scale B-field, producing a nearly one-dimensional flow of gas along the field lines. The trapped gas has a turbulent and magnetic topology that differs from that of the first phase and strongly suppresses gas accretion onto protostellar cores, reducing their masses. We speculate that this star formation bimodality may be an important ingredient to understand the origin of multiple stellar populations observed in massive globular clusters. △ Less

Submitted 1 September, 2023; v1 submitted 8 December, 2022; originally announced December 2022.

Comments: 16 Pages, 10 Figures, Published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 522, Issue 4, July 2023, Page 6203

Massive Prestellar Cores in Radiation-magneto-turbulent Simulations of Molecular Clouds

Authors: Chong-Chong He, Massimo Ricotti

Abstract: We simulate the formation and collapse of prestellar cores at few-AUえーゆー resolution in a set of radiation-magneto-hydrodynamic simulations of giant molecular clouds (GMCs) using the grid-based code RAMSES-RT. We adopt, for the first time to our best knowledge, realistic initial/boundary conditions by zooming-in onto individual massive prestellar cores within the GMC. We identify two distinct modes of… ▽ More We simulate the formation and collapse of prestellar cores at few-AUえーゆー resolution in a set of radiation-magneto-hydrodynamic simulations of giant molecular clouds (GMCs) using the grid-based code RAMSES-RT. We adopt, for the first time to our best knowledge, realistic initial/boundary conditions by zooming-in onto individual massive prestellar cores within the GMC. We identify two distinct modes of fragmentation: "quasi-spherical" and "filamentary". In both modes the fragments eventually become embedded in a quasi-steady accretion disk or toroid with radii ~ 500-5000 AUえーゆー and opening angles $H/R \sim 0.5-1$. The disks/toroids are Toomre stable but the accreted pre-existing fragments are found orbiting the outer disk, appearing as disk fragmentation. Each core converts nearly 100 percent of the gas mass into a few massive stars forming near the disk center. Large and massive disks around high-mass stars are supported by magnetic pressure in the outer disk, at radii >200-1000 AUえーゆー, and turbulent pressure in the inner disk. The most massive core accretes several times more mass than its initial mass, forming a (proto)star cluster of 8 massive stars enshrouded by a toroid, suggesting a competitive accretion scenario for ultra-high-mass star formation. We also find that the HII regions produced by a single massive star remain trapped in the dense circumstellar disks for a few hundred kiloyears, while the dynamic motions of massive stars in wide binaries or multiple systems displace the stars from the densest parts of the disk, allowing UV radiation to escape producing steady or pulsating bipolar HII regions. △ Less

Submitted 20 October, 2022; originally announced October 2022.

Comments: 20 pages, 13 figures, submitted to MNRAS. Comments welcome. Simulation movie available at https://chongchonghe.github.io/data/#cores

Kindling the First Stars: I. Dependence of Detectability of the First Stars with JWST on the Pop III Stellar Masses

Authors: Mia Sauda Bovill, Massimo Stiavelli, Alessa Ibrahim Wiggins, Massimo Ricotti, Michele Trenti

Abstract: The first Pop III stars formed out of primordial, metal free gas, in minihalos at z>20, and kickstarted the cosmic processes of reionizaton and enrichment. While these stars are likely more massive than their enriched counterparts, the current unknowns of their astrophysics include; when the first Pop III stars ignited, how massive they were, and when and how the era of the first stars ended. Inve… ▽ More The first Pop III stars formed out of primordial, metal free gas, in minihalos at z>20, and kickstarted the cosmic processes of reionizaton and enrichment. While these stars are likely more massive than their enriched counterparts, the current unknowns of their astrophysics include; when the first Pop III stars ignited, how massive they were, and when and how the era of the first stars ended. Investigating these questions requires an exploration of a multi-dimensional parameter space, including the slope of the Pop III stellar initial mass function (IMF) and the strength of the non-ionizing UV background. In this work, we present a novel model which treats both the slope and maximum mass of Pop III stars as truly free parameters while including the physics of the fragmentation of primordial gas. Our results also hint at a non-universal Pop III IMF which is dependent on the efficiency of primordial gas fragmentation. Our relatively simple model reproduces the results from hydrodynamic simulations, but with a computational efficiency which allows us to investigate the observable differences between a wide range of potential Pop III IMFs. In addition, the evolution of the number density of Pop III stars may provide insight into the evolution of the H2 dissociating background. While the slope of the Pop III IMF does not significantly affect the predicted number density of the first stars, more top heavy IMFs produce Pop III star clusters which are 2-3 magnitudes brighter than their more bottom heavy counterparts. While the Pop III star clusters are too dim for direct detection by JWST, we find they are within the reach of gravitational lensing. △ Less

Submitted 20 October, 2022; v1 submitted 18 October, 2022; originally announced October 2022.

Comments: 13 pages, 10 pages, submitted to ApJ

Devouring the Milky Way Satellites: Modeling Dwarf Galaxies with Galacticus

Authors: Sachi Weerasooriya, Mia Sauda Bovill, Andrew Benson, Alexi M. Musick, Massimo Ricotti

Abstract: Dwarf galaxies are ubiquitous throughout the universe and are extremely sensitive to various forms of internal and external feedback. Over the last two decades, the census of dwarf galaxies in the Local Group and beyond has increased markedly. While hydrodynamic simulations (e.g. FIRE II, MINT Justice League) have reproduced the observed dwarf properties down to the ultra-faints, such simulations… ▽ More Dwarf galaxies are ubiquitous throughout the universe and are extremely sensitive to various forms of internal and external feedback. Over the last two decades, the census of dwarf galaxies in the Local Group and beyond has increased markedly. While hydrodynamic simulations (e.g. FIRE II, MINT Justice League) have reproduced the observed dwarf properties down to the ultra-faints, such simulations require extensive computational resources to run. In this work, we constrain the standard physical implementations in the semi-analytic model Galacticus to reproduce the observed properties of the Milky Way satellites down to the ultra-faint dwarfs found in SDSS. We run Galacticus on merger trees from our high-resolution N-body simulation of a Milky Way analog. We determine the best-fit parameters by matching the cumulative luminosity function and luminosity-metallicity relation from both observations and hydrodynamic simulations. With the correct parameters, the standard physics in Galacticus can reproduce the observed luminosity function and luminosity-metallicity relation of the Milky Way dwarfs. In addition, we find a multi-dimensional match with half-light radii, velocity dispersions and mass-to-light ratios at z = 0 down to M_V <= -6 (L >= 10^4 L_solar). In addition to successfully reproducing the properties of the z = 0 Milky Way satellite population, our modeled dwarfs have star formation histories which are consistent with those of the Local Group dwarfs. △ Less

Submitted 27 September, 2022; originally announced September 2022.

Comments: 18 pages, 15 figures, submitted to ApJ

Ghostly Stellar Haloes and their Relationship to Ultra-faint Dwarfs

Authors: Massimo Ricotti, Emil Polisensky, Emily Cleland

Abstract: 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… ▽ More 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. △ Less

Submitted 6 September, 2022; v1 submitted 16 May, 2022; originally announced May 2022.

Comments: 18 pages, 17 figures, published on MNRAS. Movie of Fig.7 available on the online published version and at https://www.astro.umd.edu/~ricotti/NEWWEB/group.html

Journal ref: 2022MNRAS.515..302R

Population III Star Formation in an X-ray background: II. Protostellar Discs, Multiplicity and Mass Function of the Stars

Authors: Jongwon Park, Massimo Ricotti, Kazuyuki Sugimura

Abstract: Disc fragmentation plays an important role in determining the number of primordial stars (Pop III stars), their masses, and hence the initial mass function. In this second paper of a series, we explore the effect of uniform FUV H$_2$-photodissociating and X-ray radiation backgrounds on the formation of Pop~III stars using a grid of high-resolution zoom-in simulations. We find that, in an X-ray bac… ▽ More Disc fragmentation plays an important role in determining the number of primordial stars (Pop III stars), their masses, and hence the initial mass function. In this second paper of a series, we explore the effect of uniform FUV H$_2$-photodissociating and X-ray radiation backgrounds on the formation of Pop~III stars using a grid of high-resolution zoom-in simulations. We find that, in an X-ray background, protostellar discs have lower surface density and higher Toomre $Q$ parameter, so they are more stable. For this reason, X-ray irradiated discs undergo fewer fragmentations and typically produce either binary systems or low-multiplicity systems. In contrast, the cases with weak or no X-ray irradiation produce systems with a typical multiplicity of $6 \pm 3$. In addition, the most massive protostar in each system is smaller by roughly a factor of two when the disc is irradiated by X-rays, due to lower accretion rate. With these two effects combined, the initial mass function of fragments becomes more top-heavy in a strong X-ray background and is well described by a power-law with slope $1.53$ and high-mass cutoff of $61$ M$_\odot$. Without X-rays, we find a slope $0.49$ and cutoff mass of $229$ M$_\odot$. Finally, protostars migrate outward after their formation due to the accretion of high-angular momentum gas from outside and the migration is more frequent and significant in absence of X-ray irradiation. △ Less

Submitted 16 July, 2021; originally announced July 2021.

Population III Star Formation in an X-ray background: I. Critical Halo Mass of Formation and Total Mass in Stars

Authors: Jongwon Park, Massimo Ricotti, Kazuyuki Sugimura

Abstract: The first luminous objects forming in the universe produce radiation backgrounds in the FUV and X-ray bands that affect the formation of Population III stars. Using a grid of cosmological hydrodynamics zoom-in simulations, we explore the impact of the Lyman-Warner (LW) and X-ray radiation backgrounds on the critical dark matter halo mass for Population III star formation and the total mass in star… ▽ More The first luminous objects forming in the universe produce radiation backgrounds in the FUV and X-ray bands that affect the formation of Population III stars. Using a grid of cosmological hydrodynamics zoom-in simulations, we explore the impact of the Lyman-Warner (LW) and X-ray radiation backgrounds on the critical dark matter halo mass for Population III star formation and the total mass in stars per halo. We find that the LW radiation background lowers the H$_2$ fraction and delays the formation of the Population III stars. On the other hand, X-ray irradiation anticipates the redshift of collapse and reduces the critical halo mass, unless the X-ray background is too strong and gas heating shuts down gas collapse into the halos and prevents star formation. Therefore, an X-ray background can increase the number of dark matter halos forming Population III stars by about a factor of ten, but the total mass in stars forming in each halo is reduced. This is because X-ray radiation increases the molecular fraction and lowers the minimum temperature of the collapsing gas (or equivalently the mass of the quasi-hydrostatic core) and therefore slows down the accretion of the gas onto the central protostar. △ Less

Submitted 16 July, 2021; originally announced July 2021.

Multi-wavelength detectability of isolated black holes in the Milky Way

Authors: Francesca Scarcella, Daniele Gaggero, Riley Connors, Julien Manshanden, Massimo Ricotti, Gianfranco Bertone

Abstract: Isolated black holes in our Galaxy have eluded detection so far. We present here a comprehensive study on the detectability of isolated stellar-mass astrophysical black holes that accrete interstellar gas from molecular clouds in both the local region and the Central Molecular Zone. We adopt a state-of-the-art model for the accretion physics backed up by numerical simulations, and study the number… ▽ More Isolated black holes in our Galaxy have eluded detection so far. We present here a comprehensive study on the detectability of isolated stellar-mass astrophysical black holes that accrete interstellar gas from molecular clouds in both the local region and the Central Molecular Zone. We adopt a state-of-the-art model for the accretion physics backed up by numerical simulations, and study the number of observable sources in both the radio and X-ray band, as a function of a variety of parameters. We discuss in particular the impact of the astrophysical uncertainties on our prediction for the number of bright X-ray sources in the central region of the Galaxy. We finally consider future developments in the radio domain, and assess the potential of SKA to detect a population of astrophysical black holes accreting gas in our Galaxy. △ Less

Submitted 18 December, 2020; originally announced December 2020.

Biconical-dominated Accretion Flow onto Seed Black Holes in a Hyper-accretion Regime

Authors: KwangHo Park, John H. Wise, Tamara Bogdanović, Massimo Ricotti

Abstract: Hyperaccretion occurs when the gas inflow rate onto a black hole (BH) is so high that the radiative feedback cannot reverse the accretion flow. This extreme process is a promising mechanism for the rapid growth of seed BHs in the early universe, which can explain high-redshift quasars powered by billion solar mass BHs. In theoretical models, spherical symmetry is commonly adopted for hyperaccretio… ▽ More Hyperaccretion occurs when the gas inflow rate onto a black hole (BH) is so high that the radiative feedback cannot reverse the accretion flow. This extreme process is a promising mechanism for the rapid growth of seed BHs in the early universe, which can explain high-redshift quasars powered by billion solar mass BHs. In theoretical models, spherical symmetry is commonly adopted for hyperaccretion flows; however, the sustainability of such structures on timescales corresponding to the BH growth has not been addressed yet. Here we show that stochastic interactions between the ionizing radiation from the BH and nonuniform accretion flow can lead to the formation of a rotating gas disk around the BH. Once the disk forms, the supply of gas to the BH preferentially occurs via biconical-dominated accretion flow perpendicular to the disk, avoiding the centrifugal barrier of the disk. Biconical-dominated accretion flows from opposite directions collide in the vicinity of the BH supplying high-density, low angular momentum gas to the BH, whereas most of the gas with nonnegligible angular momentum is deflected to the rotationally supported outflowing decretion disk. The disk becomes reinforced progressively as more mass from the biconical flow transfers to the disk and some of the outflowing gas from the disk is redirected to the biconical accretion funnels through a meridional structure. This axisymmetric hydrodynamic structure of a biconical-dominated accretion flow and decretion disk continues to provide uninterrupted flow of high-density gas to the BH. △ Less

Submitted 25 November, 2020; v1 submitted 11 June, 2020; originally announced June 2020.

Comments: 14 pages, 10 figures, accepted for publication in ApJ

Structure and Instability of the Ionization Fronts around Moving Black Holes

Authors: Kazuyuki Sugimura, Massimo Ricotti

Abstract: 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… ▽ More 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. △ Less

Submitted 13 June, 2020; v1 submitted 12 March, 2020; originally announced March 2020.

Comments: 14 pages, 9 figures, published in MNRAS, movie available at https://www.youtube.com/watch?v=pj1dy-0AM0o

Journal ref: MNRAS 495, 2966-2978 (2020)

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

Authors: Chong-Chong He, Massimo Ricotti, Sam Geen

Abstract: 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 a… ▽ More 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. △ Less

Submitted 16 January, 2020; originally announced January 2020.

Comments: 17 pages, 13 figures, accepted for publication in MNRAS

Simulating Star Clusters Across Cosmic Time: I. Initial Mass Function, Star Formation Rates and Efficiencies

Authors: Chong-Chong He, Massimo Ricotti, Sam Geen

Abstract: We present radiation-magneto-hydrodynamic simulations of star formation in self-gravitating, turbulent molecular clouds, modeling the formation of individual massive stars, including their UV radiation feedback. The set of simulations have cloud masses between $m_{\rm gas}=10^3$~M$_\odot$ to $3 \times 10^5$~M$_\odot$ and gas densities typical of clouds in the local universe (… ▽ More We present radiation-magneto-hydrodynamic simulations of star formation in self-gravitating, turbulent molecular clouds, modeling the formation of individual massive stars, including their UV radiation feedback. The set of simulations have cloud masses between $m_{\rm gas}=10^3$~M$_\odot$ to $3 \times 10^5$~M$_\odot$ and gas densities typical of clouds in the local universe ($\overline n_{\rm gas} \sim 1.8\times 10^2$~cm$^{-3}$) and 10$\times$ and 100$\times$ denser, expected to exist in high-redshift galaxies. The main results are: {\it i}) The observed Salpeter power-law slope and normalisation of the stellar initial mass function at the high-mass end can be reproduced if we assume that each star-forming gas clump (sink particle) fragments into stars producing on average a maximum stellar mass about $40\%$ of the mass of the sink particle, while the remaining $60\%$ is distributed into smaller mass stars. Assuming that the sinks fragment according to a power-law mass function flatter than Salpeter, with log-slope $0.8$, satisfy this empirical prescription. {\it ii}) The star formation law that best describes our set of simulation is $dρろー_*/dt \propto ρろー_{gas}^{1.5}$ if $\overline n_{gas}<n_{cri}\approx 10^3$~cm$^{-3}$, and $dρろー_*/dt \propto ρろー_{\rm gas}^{2.5}$ otherwise. The duration of the star formation episode is roughly $6$ cloud's sound crossing times (with $c_s=10$~km/s). {\it iii}) The total star formation efficiency in the cloud is $f_*=2\% (m_{\rm gas}/10^4~M_\odot)^{0.4}(1+\overline n_{\rm gas}/n_{\rm cri})^{0.91}$, for gas at solar metallicity, while for metallicity $Z<0.1$~Z$_\odot$, based on our limited sample, $f_*$ is reduced by a factor of $\sim 5$. {\it iv)} The most compact and massive clouds appear to form globular cluster progenitors, in the sense that star clusters remain gravitationally bound after the gas has been expelled. △ Less

Submitted 26 September, 2019; v1 submitted 16 April, 2019; originally announced April 2019.

Comments: 21 pages, 18 figures, Published in MNRAS. References added; Fig. 5 and Fig. 15 added

Electromagnetic probes of primordial black holes as dark matter

Authors: Y. Ali-Haimoud, S. Clesse, J. Garcia-Bellido, A. Kashlinsky, L. Wyrzykowski, A. Achucarro, L. Amendola, J. Annis, A. Arbey, R. G. Arendt, F. Atrio-Barandela, N. Bellomo, K. Belotsky, J-L. Bernal, S. Bird, V. Bozza, C. Byrnes, S. Calchi Novati, F. Calore, B. J. Carr, J. Chluba, I. Cholis, A. Cieplak, P. Cole, I. Dalianis , et al. (69 additional authors not shown)

Abstract: The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence wit… ▽ More The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing. △ Less

Submitted 12 March, 2019; v1 submitted 11 March, 2019; originally announced March 2019.

Comments: Science whitepaper submitted to the Astro2020 Decadal Survey

Populations behind the source-subtracted cosmic infrared background anisotropies

Authors: A. Kashlinsky, R. G. Arendt, M. Ashby, F. Atrio-Barandela, V. Bromm, N. Cappelluti, S. Clesse, A. Comastri, J-G. Cuby, S. Driver, G. Fazio, A. Ferrara, A. Finoguenov, D. Fixsen, J. Garcia-Bellido, G. Hasinger, K. Helgason, R. J. Hill, R. Jansen, J. Kruk, J. Mather, P. Natarajan, N. Odegard, T. Reiprich, M. Ricotti , et al. (5 additional authors not shown)

Abstract: While the upcoming telescopes will reveal correspondingly fainter, more distant galaxies, a question will persist: what more is there that these telescopes cannot see? One answer is the source-subtracted Cosmic Infrared Background (CIB). The CIB is comprised of the collective light from all sources remaining after known, resolved sources are accounted for. Ever-more-sensitive surveys will identify… ▽ More While the upcoming telescopes will reveal correspondingly fainter, more distant galaxies, a question will persist: what more is there that these telescopes cannot see? One answer is the source-subtracted Cosmic Infrared Background (CIB). The CIB is comprised of the collective light from all sources remaining after known, resolved sources are accounted for. Ever-more-sensitive surveys will identify the brightest of these, allowing them to be removed, and - like peeling layers off an onion - reveal deeper layers of the CIB. In this way it is possible to measure the contributions from populations not accessible to direct telescopic observation. Measurement of fluctuations in the source-subtracted CIB, i.e., the spatial power spectrum of the CIB after subtracting resolved sources, provides a robust means of characterizing its faint, and potentially new, populations. Studies over the past 15 years have revealed source-subtracted CIB fluctuations on scales out to ~100' which cannot be explained by extrapolating from known galaxy populations. Moreover, they appear highly coherent with the unresolved Cosmic X-ray Background, hinting at a significant population of accreting black holes among the CIB sources. Characterizing the source-subtracted CIB with high accuracy, and thereby constraining the nature of the new populations, is feasible with upcoming instruments and would produce critically important cosmological information in the next decade. New coextensive deep and wide-area near-infrared, X-ray, and microwave surveys will bring decisive opportunities to examine, with high fidelity, the spatial spectrum and origin of the CIB fluctuations and their cross-correlations with cosmic microwave and X-ray backgrounds, and determine the formation epochs and the nature of the new sources (stellar nucleosynthetic or accreting black holes). △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: Science whitepaper submitted to the Astro2020 Decadal Survey

Multi-wavelength astronomical searches for primordial black holes

Authors: Julien Manshanden, Daniele Gaggero, Gianfranco Bertone, Riley M. T. Connors, Massimo Ricotti

Abstract: If primordial black holes of $\mathcal{O}(1-100) M_{\odot}$ constitute a significant portion of the dark matter in the Universe, they should be very abundant in our Galaxy. We present here a detailed analysis of the radio and X-ray emission that these objects are expected to produce due to the accretion of gas from the interstellar medium. With respect to previous studies, we relax the assumption… ▽ More If primordial black holes of $\mathcal{O}(1-100) M_{\odot}$ constitute a significant portion of the dark matter in the Universe, they should be very abundant in our Galaxy. We present here a detailed analysis of the radio and X-ray emission that these objects are expected to produce due to the accretion of gas from the interstellar medium. With respect to previous studies, we relax the assumption of a monochromatic mass function, and introduce an improved treatment of the physics of gas accretion onto isolated, moving compact objects, based on a set of state-of-the-art numerical simulations. By comparing our predictions with known radio and X-ray sources in the Galactic center region, we show that the maximum relic density of primordial black holes in the mass range of interest is $\sim 10^{-3}$ smaller than that of dark matter. The new upper bound is two orders of magnitude stronger with respect to previous results, based on a conservative phenomenological treatment of the accretion physics. We also provide a comprehensive critical discussion on the reliability of this bound, and on possible future developments in the field. We argue in particular that future multi-wavelength searches will soon start to probe the galactic population of astrophysical black holes. △ Less

Submitted 2 July, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: 21 pages, 5 figures, final version accepted for publication in JCAP

Ghostly Halos in Dwarf Galaxies: Constraints on the Star Formation Efficiency before Reionization

Authors: Hoyoung Kang, Massimo Ricotti

Abstract: Stellar halos observed around normal galaxies are extended and faint stellar structures formed by debris of tidally disrupted dwarf galaxies accreted over time by the host galaxy. Around dwarf galaxies, these stellar halos may not exist if all the accreted satellites are dark halos without stars. However, if a stellar halo is found in sufficiently small mass dwarfs, the whole stellar halo is compo… ▽ More Stellar halos observed around normal galaxies are extended and faint stellar structures formed by debris of tidally disrupted dwarf galaxies accreted over time by the host galaxy. Around dwarf galaxies, these stellar halos may not exist if all the accreted satellites are dark halos without stars. However, if a stellar halo is found in sufficiently small mass dwarfs, the whole stellar halo is composed of tidal debris of fossil galaxies, and we refer to it as ghostly halo. Fossil galaxies are called so because they formed most of their stars before the epoch of reionization, and have been identified as the ultra-faint dwarf galaxies found around the Milky Way and M31. In this paper we carry out semi-analytical simulations to characterize the sizes and stellar masses of ghostly stellar halos in dwarf galaxies as a function of their dark matter halo mass. By comparing the models to observations of six isolated dwarf galaxies in the Local Group showing evidence of extended stellar halos, we are able to constrain the star formation efficiency in fossil galaxies. We find that at redshift z~6, dark matter halos in the mass range 10^7-10^9 M_sun have a mean star formation efficiency f_* = M_*/M_{dm} ~ 0.05%-0.1% nearly constant as a function of the dark matter halo mass. There is a tentative indication of a sharp increase of f_* at M_{dm} ~ 10^6-10^7 M_sun, but it is very uncertain and based on only one dwarf galaxy (Leo~T). △ Less

Submitted 2 October, 2018; originally announced October 2018.

Comments: 15 pages, 7 figures. Submitted to MNRAS

ARC: Adaptive Ray-tracing with CUDA, a New Ray Tracing Code for Parallel GPUs

Authors: Blake Hartley, Massimo Ricotti

Abstract: We present the methodology of a photon-conserving, spatially-adaptive, ray-tracing radiative transfer algorithm, designed to run on multiple parallel Graphic Processing Units (GPUs). Each GPU has thousands computing cores, making them ideally suited to the task of tracing independent rays. This ray-tracing implementation has speed competitive with approximate momentum methods, even with thousands… ▽ More We present the methodology of a photon-conserving, spatially-adaptive, ray-tracing radiative transfer algorithm, designed to run on multiple parallel Graphic Processing Units (GPUs). Each GPU has thousands computing cores, making them ideally suited to the task of tracing independent rays. This ray-tracing implementation has speed competitive with approximate momentum methods, even with thousands of ionization sources, without sacrificing accuracy and resolution. Here, we validate our implementation with the selection of tests presented in the "cosmological radiative transfer codes comparison project," to demonstrate the correct behavior of the code. We also present a selection of benchmarks to demonstrate the performance and computational scaling of the code. As expected, our method scales linearly with the number of sources and with the square of the dimension of the 3D computational grid. Our current implementation is scalable to an arbitrary number of nodes possessing GPUs, but is limited to a uniform resolution 3D grid. Cosmological simulations of reionization with tens of thousands of radiation sources and intergalactic volumes sampled with 1024$^3$ grid points take about 30 days on 64 GPUs to reach complete reionization. △ Less

Submitted 18 July, 2018; originally announced July 2018.

Comments: 18 pages, 13 figures

Dusty Gas Accretion onto Massive Black Holes and Infrared Diagnosis of the Eddington Ratio

Authors: Hidenobu Yajima, Massimo Ricotti, KwangHo Park, Kazuyuki Sugimura

Abstract: Evidence for dust around supermassive black holes (SMBHs) in the early Universe is strongly suggested by recent observations. However, the accretion mechanism of SMBHs in dusty gas is not well understood yet. We investigate the growth of intermediate-mass black-holes (IMBHs) of $\sim 10^{5}~M_{\odot}$ in dusty clouds by using one-dimensional radiative-hydrodynamics simulations. We find that the ac… ▽ More Evidence for dust around supermassive black holes (SMBHs) in the early Universe is strongly suggested by recent observations. However, the accretion mechanism of SMBHs in dusty gas is not well understood yet. We investigate the growth of intermediate-mass black-holes (IMBHs) of $\sim 10^{5}~M_{\odot}$ in dusty clouds by using one-dimensional radiative-hydrodynamics simulations. We find that the accretion of dusty gas onto IMBHs proceeds gently with small fluctuations of the accretion rate, whereas that of pristine gas causes more violent periodic bursts. At dust-to-gas mass ratios similar to the solar neighborhood, the time averaged luminosity becomes smaller than that for primordial gas by one order of magnitude and the time-averaged Eddington ratio ranges from $\sim 10^{-4}$ to $\sim 10^{-2}$ in clouds with initial gas densities of $n_{\rm H} = 10 - 1000~\rm cm^{-3}$. Our calculations show that the effect of dust opacity alone is secondary compared to the radiation pressure on dust in regulating the BH growth. We also derive spectral energy distributions at IR bands by calculating dust thermal emission and show that the flux ratio between $λらむだ\lesssim 20~\rm μみゅーm$ and $\gtrsim 100~\rm μみゅーm$ is closely related to the Eddington ratio. Thermal emission from hot dust near the BH dominates only during the high accretion phase, producing higher flux density at $\lesssim 20~\rm μみゅーm$. Therefore, we suggest that the combinations of MIR observations by JWST and FIR observation by ALMA or Spitzer can be used to estimate the Eddington ratio of massive BHs. △ Less

Submitted 18 April, 2017; originally announced April 2017.

Comments: 13 pages, 12 figures, submitted to ApJ

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

Authors: Mia Sauda Bovill, Thomas H. Puzia, Massimo Ricotti, Matthew A. Taylor

Abstract: 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) sugge… ▽ More 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. △ Less

Submitted 24 August, 2016; originally announced August 2016.

Comments: 11 pages, 10 figure, 1 table. Accepted to ApJ

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

Authors: Massimo Ricotti, Owen H. Parry, Nickolay Y. Gnedin

Abstract: 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 spheroid… ▽ More 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. △ Less

Submitted 14 July, 2016; originally announced July 2016.

Comments: 15 pages, 12 figure, submitted to ApJ

X-ray Twinkles and Pop III Stars

Authors: Massimo Ricotti

Abstract: 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 numbe… ▽ More 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. △ Less

Submitted 14 July, 2016; originally announced July 2016.

Comments: 10 pages, 4 figures. Accepted for publication on MNRAS

Modeling Reionization in a Bursty Universe

Authors: Blake Hartley, Massimo Ricotti

Abstract: 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 re… ▽ More 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. △ Less

Submitted 19 February, 2016; originally announced February 2016.

Comments: 17 pages, 11 figures, 1 table. Submitted to MNRAS

Bulge-driven Fueling of Seed Black Holes

Authors: KwangHo Park, Massimo Ricotti, Priyamvada Natarajan, Tamara Bogdanović, John H. Wise

Abstract: 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, whe… ▽ More 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. △ Less

Submitted 10 December, 2015; originally announced December 2015.

Comments: 10 pages, 10 Figures, 1 Table, submitted to ApJ

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

Authors: Kári Helgason, Massimo Ricotti, Alexander Kashlinsky, Volker Bromm

Abstract: 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… ▽ More 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. △ Less

Submitted 23 September, 2015; v1 submitted 27 May, 2015; originally announced May 2015.

Comments: 13 pages, 11 figures. Accepted for publication in MNRAS

Fingerprints of the initial conditions on the density profiles of cold and warm dark matter haloes

Authors: E. Polisensky, M. Ricotti

Abstract: We use N-body simulations of dark matter haloes in cold dark matter (CDM) and a large set of different warm dark matter (WDM) cosmologies to demonstrate that the spherically averaged density profile of dark matter haloes has a shape that depends on the power spectrum of matter perturbations. Density profiles are steeper in WDM but become shallower at scales less than one percent of the virial radi… ▽ More We use N-body simulations of dark matter haloes in cold dark matter (CDM) and a large set of different warm dark matter (WDM) cosmologies to demonstrate that the spherically averaged density profile of dark matter haloes has a shape that depends on the power spectrum of matter perturbations. Density profiles are steeper in WDM but become shallower at scales less than one percent of the virial radius. Virialization isotropizes the velocity dispersion in the inner regions of the halo but does not erase the memory of the initial conditions in phase space. The location of the observed deviations from CDM in the density profile and in phase space can be directly related to the ratio between the halo mass and the filtering mass and are most evident in small mass haloes, even for a 34 keV thermal relic WDM. The rearrangement of mass within the haloes supports analytic models of halo structure that include angular momentum. We also find evidence of a dependence of the slope of the inner density profile in CDM cosmologies on the halo mass with more massive haloes exhibiting steeper profiles, in agreement with the model predictions and with previous simulation results. Our work complements recent studies of microhaloes near the filtering scale in CDM and strongly argue against a universal shape for the density profile. △ Less

Submitted 8 April, 2015; originally announced April 2015.

Comments: 14 pages, 11 figures. Submitted to MNRAS 05/2014, accepted for publication 04/2015

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

Authors: KwangHo Park, Massimo Ricotti, Tiziana Di Matteo, Christopher S. Reynolds

Abstract: 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 b… ▽ More 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). △ Less

Submitted 21 October, 2014; v1 submitted 7 August, 2014; originally announced August 2014.

Comments: 6 pages, 5 figures, accepted for publication in MNRAS

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

Authors: Harley Katz, Massimo Ricotti

Abstract: 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… ▽ More 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 △ Less

Submitted 23 June, 2014; originally announced June 2014.

Comments: 22 pages, 17 figures, submitted to MNRAS, comments welcome

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

Authors: Kari Helgason, Nico Cappelluti, Guenther Hasinger, Alexander Kashlinsky, Massimo Ricotti

Abstract: 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… ▽ More 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. △ Less

Submitted 26 February, 2014; v1 submitted 5 November, 2013; originally announced November 2013.

Comments: 18 pages, 12 figures. Accepted for publication in ApJ

Rayleigh-Taylor Instability at Ionization Fronts: Perturbation Analysis

Authors: Massimo Ricotti

Abstract: 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 i… ▽ More 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. △ Less

Submitted 5 October, 2013; originally announced October 2013.

Comments: Accepted for publication in MNRAS. 9 pages

Massive Milky Way Satellites in Cold and Warm Dark Matter: Dependence on Cosmology

Authors: E. Polisensky, M. Ricotti

Abstract: We investigate the claim that the largest subhaloes in high resolution dissipationless cold dark matter (CDM) simulations of the Milky Way are dynamically inconsistent with observations of its most luminous satellites. We find that the inconsistency is largely attributable to the large values of σしぐま_8 and n_s adopted in the discrepant simulations producing satellites that form too early and therefor… ▽ More We investigate the claim that the largest subhaloes in high resolution dissipationless cold dark matter (CDM) simulations of the Milky Way are dynamically inconsistent with observations of its most luminous satellites. We find that the inconsistency is largely attributable to the large values of σしぐま_8 and n_s adopted in the discrepant simulations producing satellites that form too early and therefore are too dense. We find the tension between observations and simulations adopting parameters consistent with WMAP9 is greatly diminished making the satellites a sensitive test of CDM. We find the Via Lactea II halo to be atypical for haloes in a WMAP3 cosmology, a discrepancy that we attribute to its earlier formation epoch than the mean for its mass. We also explore warm dark matter (WDM) cosmologies for 1--4 keV thermal relics. In 1 keV cosmologies subhaloes have circular velocities at kpc scales ~ 60% lower than their CDM counterparts, but are reduced by only 10% in 4 keV cosmologies. Since relic masses < 2-3 keV are ruled out by constraints from the number of Milky Way satellites and Lyman-αあるふぁ forest, WDM has a minor effect in reducing the densities of massive satellites. Given the uncertainties on the mass and formation epoch of the Milky Way, the need for reducing the satellite densities with baryonic effects or WDM is alleviated. △ Less

Submitted 1 October, 2013; originally announced October 2013.

Comments: 11 pages, 7 figures, submitted to MNRAS

Journal ref: 2014, MNRAS, 437, 2922

Rayleigh-Taylor instability of ionization front around black holes

Authors: KwangHo Park, Massimo Ricotti, Tiziana Di Matteo, Christopher S. Reynolds

Abstract: 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 f… ▽ More 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. △ Less

Submitted 13 November, 2013; v1 submitted 23 August, 2013; originally announced August 2013.

Comments: 9 pages, 8 figures, accepted for publication in MNRAS

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

Authors: Harley Katz, Massimo Ricotti

Abstract: 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 h… ▽ More 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. △ Less

Submitted 26 November, 2012; originally announced November 2012.

Comments: 12 pages, 6 figures, Submitted to MNRAS

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

Authors: KwangHo Park, Massimo Ricotti

Abstract: 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 w… ▽ More 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. △ Less

Submitted 27 March, 2013; v1 submitted 2 November, 2012; originally announced November 2012.

Comments: Accepted for publication in ApJ. 12 pages in emulated ApJ format. 11 figures. For associated MPEG files see http://www.andrew.cmu.edu/user/kwanghop/research.html

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

Authors: Kari Helgason, Massimo Ricotti, Alexander Kashlinsky

Abstract: 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 r… ▽ More 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. △ Less

Submitted 21 August, 2012; originally announced August 2012.

Comments: 3 pages, 2 figures. To appear in proceedings of First Stars IV meeting (Kyoto, Japan; 2012)