Showing 1–50 of 147 results for author: Shapiro, P

Multiple Beads-on-a-string: Dark Matter-Deficient Galaxy Formation in a Mini-bullet Satellite-satellite Galaxy Collision

Authors: Joohyun Lee, Eun-jin Shin, Ji-hoon Kim, Paul R. Shapiro, Eunwoo Chung

Abstract: Dark matter-deficient galaxies (DMDGs) discovered in the survey of ultra-diffuse galaxies (UDGs), in apparent conflict with standard CDM, may be produced by high-velocity galaxy-galaxy collisions, the $\textit{Mini-bullet}$ scenario. Recent observations of an aligned trail of $7-11$ UDGs near NGC1052, including DMDGs DF2 and DF4, suggesting a common formation event, $\sim8.9\pm1.5$ Gyr ago, provid… ▽ More Dark matter-deficient galaxies (DMDGs) discovered in the survey of ultra-diffuse galaxies (UDGs), in apparent conflict with standard CDM, may be produced by high-velocity galaxy-galaxy collisions, the $\textit{Mini-bullet}$ scenario. Recent observations of an aligned trail of $7-11$ UDGs near NGC1052, including DMDGs DF2 and DF4, suggesting a common formation event, $\sim8.9\pm1.5$ Gyr ago, provide a test. Hydro/N-body simulations, supplemented by galaxy orbit integrations, demonstrate that satellite-satellite collisions outside the host-galaxy virial radius can reproduce the observed UDGs in the NGC1052 group. A trail of $\sim10$ DMDGs is shown to form, including two massive ones that replicate the observed motions of DF2 and DF4. The linear relation, $v=Ax+v_{0}$, conjectured previously to relate positions ($x$) and velocities ($v$) of the aligned DMDGs as a signature of the collision event, is approximately obeyed, but individual DMDGs can deviate significantly from it. The progenitors whose collision spawned the trail of DMDGs survive the collision without, themselves, becoming DMDGs. We predict one progenitor is located at the end of the trail, testable by observing the difference between its stars, formed pre-collision, from those of the DMDGs, formed post-collision. By contrast, stellar ages and metallicities of the DMDGs are nearly identical. We further offer a hint that the tidal field of host NGC1052 may contribute to making DMDGs diffuse. $Λ$CDM simulation in a 100 cMpc box finds our required initial conditions $\sim10$ times at $z<3$. These results indicate current observations are consistent with the $\textit{Mini-bullet}$ scenario. △ Less

Submitted 9 February, 2024; v1 submitted 18 December, 2023; originally announced December 2023.

Comments: 27 pages, 14 figures, minor revisions for clarification only in response to referee report, + 2 new figures. For associated movies and information, see https://joohyun-lee.github.io/research/DMDG3/

Report number: UT-WI-11-2024

Journal ref: ApJ, 966, 72 (2024)

Impact of Self-shielding Minihalos on the Ly$α$ Forest at High Redshift

Authors: Hyunbae Park, Zarija Lukić, Jean Sexton, Marcelo Alvarez, Paul R. Shapiro

Abstract: Dense gas in minihalos with masses of $10^6-10^8~M_\odot$ can shield themselves from reionization for $\sim100$ Myr after being exposed to the UV background. These self-shielded systems, often unresolved in cosmological simulations, can introduce strong absorption in quasar spectra. This paper is the first systematic study on the impact of these systems on the Ly$α$ forest. We first derive the HI… ▽ More Dense gas in minihalos with masses of $10^6-10^8~M_\odot$ can shield themselves from reionization for $\sim100$ Myr after being exposed to the UV background. These self-shielded systems, often unresolved in cosmological simulations, can introduce strong absorption in quasar spectra. This paper is the first systematic study on the impact of these systems on the Ly$α$ forest. We first derive the HI column density profile of photoevaporating minihalos by conducting 1D radiation-hydrodynamics simulations. We utilize these results to estimate the Ly$α$ opacity from minihalos in a large-scale simulation that cannot resolve self-shielding. When the ionization rate of the background radiation is $0.03\times10^{-12}~{\rm s}^{-1}$, as expected near the end of reionization at $z\sim5.5$, we find that the incidence rate of damped Ly$α$ absorbers increases by a factor of $\sim2-4$ compared to at $z=4.5$. The Ly$α$ flux is, on average, suppressed by $\sim 3\%$ of its mean due to minihalos. The absorption features enhance the 1D power spectrum up to $\sim5\%$ at $k\sim0.1~h~{\rm Mpc}^{-1}~({\rm or}~10^{-3}~{\rm km}^{-1}~{\rm s})$, which is comparable to the enhancement caused by inhomogeneous reionization. The flux is particularly suppressed in the vicinity of large halos along the line-of-sight direction at separations of up to $10~h^{-1}~{\rm Mpc}$ at $r_\perp\lesssim2~h^{-1}~{\rm Mpc}$. However, these effects become much smaller for higher ionizing rates ($\gtrsim0.3\times10^{-12}~{\rm s}^{-1}$) expected in the post-reionization Universe. Our findings highlight the need to consider minihalo absorption when interpreting the Ly$α$ forest at $z\gtrsim5.5$. Moreover, the sensitivity of these quantities to the ionizing background intensity can be exploited to constrain the intensity itself. △ Less

Submitted 15 June, 2024; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 16 pages; 11 figures; Accepted for publication in ApJ

The Effect of Boundary Conditions on Structure Formation in Fuzzy Dark Matter

Authors: Iván Alvarez-Ríos, Francisco S. Guzmán, Paul R. Shapiro

Abstract: We illustrate the effect of boundary conditions on the evolution of structure in Fuzzy Dark Matter. Scenarios explored include the evolution of single, ground-state equilibrium solutions of the Schrödinger-Poisson system, the relaxation of a Gaussian density fluctuation, mergers of two equilibrium configurations, and the random merger of many solitons. For comparison, each scenario is evolved twic… ▽ More We illustrate the effect of boundary conditions on the evolution of structure in Fuzzy Dark Matter. Scenarios explored include the evolution of single, ground-state equilibrium solutions of the Schrödinger-Poisson system, the relaxation of a Gaussian density fluctuation, mergers of two equilibrium configurations, and the random merger of many solitons. For comparison, each scenario is evolved twice, with isolation boundary conditions and periodic boundary conditions, the two commonly used to simulate isolated systems and structure formation, respectively. Replacing isolation boundary conditions by periodic boundary conditions changes the domain topology and dynamics of each scenario, by affecting the outcome of gravitational cooling. With periodic boundary conditions, the ground-state equilibrium solution and Gaussian fluctuation each evolve toward the single equilibrium solitonic core of the isolated case, but surrounded by a tail, unlike the isolated versions. The case of head-on, binary mergers illustrates additional effects, caused by the pull suffered by the system due to the infinite network of periodic images along each direction of the domain. Binary merger with angular momentum is the first scenario we found in which the tail has a polynomial profile when using a periodic domain. Finally, the 3D merger of many, randomly-placed solitonic cores of different mass makes a solitonic core surrounded by a tail with power-law-like density profile, for periodic boundary conditions, while producing a core with a much sharper fall-off in the isolated case. This suggests that the conclusion of earlier work that the ground-state equilibrium solution is an attractor for the asymptotic state is true even in 3D and for general circumstances, but only if gravitational cooling is able to carry mass and energy off to infinity, which isolation boundary conditions allow, but periodic ones do not. △ Less

Submitted 1 June, 2023; v1 submitted 6 April, 2023; originally announced April 2023.

Comments: 17 pages, 18 figures, modifications applied in order to match the accepted version in Phys. Rev. D

Journal ref: Phys. Rev. D 107, 123524 (2023)

Cosmic Variance and the Inhomogeneous UV Luminosity Function of Galaxies During Reionization

Authors: Taha Dawoodbhoy, Paul R. Shapiro, Pierre Ocvirk, Joseph S. W. Lewis, Dominique Aubert, Jenny G. Sorce, Kyungjin Ahn, Ilian T. Iliev, Hyunbae Park, Romain Teyssier, Gustavo Yepes

Abstract: When the first galaxies formed and starlight escaped into the intergalactic medium to reionize it, galaxy formation and reionization were both highly inhomogeneous in time and space, and fully-coupled by mutual feedback. To show how this imprinted the UV luminosity function (UVLF) of reionization-era galaxies, we use our large-scale, radiation-hydrodynamics simulation CoDa II to derive the time- a… ▽ More When the first galaxies formed and starlight escaped into the intergalactic medium to reionize it, galaxy formation and reionization were both highly inhomogeneous in time and space, and fully-coupled by mutual feedback. To show how this imprinted the UV luminosity function (UVLF) of reionization-era galaxies, we use our large-scale, radiation-hydrodynamics simulation CoDa II to derive the time- and space-varying halo mass function and UVLF, from $z\simeq6$-15. That UVLF correlates strongly with local reionization redshift: earlier-reionizing regions have UVLFs that are higher, more extended to brighter magnitudes, and flatter at the faint end than later-reionizing regions observed at the same $z$. In general, as a region reionizes, the faint-end slope of its local UVLF flattens, and, by $z=6$ (when reionization ended), the global UVLF, too, exhibits a flattened faint-end slope, `rolling-over' at $M_\text{UV}\gtrsim-17$. CoDa II's UVLF is broadly consistent with cluster-lensed galaxy observations of the Hubble Frontier Fields at $z=6$-8, including the faint end, except for the faintest data point at $z=6$, based on one galaxy at $M_\text{UV}=-12.5$. According to CoDa II, the probability of observing the latter is $\sim5\%$. However, the effective volume searched at this magnitude is very small, and is thus subject to significant cosmic variance. We find that previous methods adopted to calculate the uncertainty due to cosmic variance underestimated it on such small scales by a factor of 2-4, primarily by underestimating the variance in halo abundance when the sample volume is small. △ Less

Submitted 8 August, 2023; v1 submitted 16 February, 2023; originally announced February 2023.

Comments: 16 pages, 16 figures, accepted by MNRAS 07/20/23, comments welcome

Reionization time of the Local Group and Local-Group-like halo pairs

Authors: Jenny G. Sorce, Pierre Ocvirk, Dominique Aubert, Stefan Gottloeber, Paul R. Shapiro, Taha Dawoodbhoy, Gustavo Yepes, Kyungjin Ahn, Ilian T. Iliev, Joseph S. W. Lewis

Abstract: Patchy cosmic reionization resulted in the ionizing UV background asynchronous rise across the Universe. The latter might have left imprints visible in present day observations. Several numerical simulation-based studies show correlations between reionization time and overdensities and object masses today. To remove the mass from the study, as it may not be the sole important parameter, this paper… ▽ More Patchy cosmic reionization resulted in the ionizing UV background asynchronous rise across the Universe. The latter might have left imprints visible in present day observations. Several numerical simulation-based studies show correlations between reionization time and overdensities and object masses today. To remove the mass from the study, as it may not be the sole important parameter, this paper focuses solely on the properties of paired halos within the same mass range as the Milky Way. For this purpose, it uses CoDaII, a fully-coupled radiation hydrodynamics reionization simulation of the local Universe. This simulation holds a halo pair representing the Local Group, in addition to other pairs, sharing similar mass, mass ratio, distance separation and isolation criteria but in other environments, alongside isolated halos within the same mass range. Investigations of the paired halo reionization histories reveal a wide diversity although always inside-out given our reionization model. Within this model, halos in a close pair tend to be reionized at the same time but being in a pair does not bring to an earlier time their mean reionization. The only significant trend is found between the total energy at z = 0 of the pairs and their mean reionization time: pairs with the smallest total energy (bound) are reionized up to 50 Myr earlier than others (unbound). Above all, this study reveals the variety of reionization histories undergone by halo pairs similar to the Local Group, that of the Local Group being far from an average one. In our model, its reionization time is ~625 Myr against 660+/-4 Myr (z~8.25 against 7.87+/-0.02) on average. △ Less

Submitted 26 July, 2022; originally announced July 2022.

Comments: Accepted for publication in MNRAS, 11 pages, 5 figures

The Physics of Light Relics

Authors: Cora Dvorkin, Joel Meyers, Peter Adshead, Mustafa Amin, Carlos A. Argüelles, Thejs Brinckmann, Emanuele Castorina, Timothy Cohen, Nathaniel Craig, David Curtin, Francis-Yan Cyr-Racine, Peizhi Du, Lloyd Knox, Bohua Li, Marilena Loverde, Kaloian Lozanov, Julian B. Muñoz, Katelin Schutz, Paul Shapiro, Benjamin Wallisch, Zachary J. Weiner, Weishuang Linda Xu

Abstract: Many well-motivated extensions of the Standard Model predict the existence of new light species that may have been produced in the early universe. Prominent examples include axions, sterile neutrinos, gravitinos, dark photons, and more. The gravitational influence of light relics leaves imprints in the cosmic microwave background fluctuations, the large-scale structure of the universe and the prim… ▽ More Many well-motivated extensions of the Standard Model predict the existence of new light species that may have been produced in the early universe. Prominent examples include axions, sterile neutrinos, gravitinos, dark photons, and more. The gravitational influence of light relics leaves imprints in the cosmic microwave background fluctuations, the large-scale structure of the universe and the primordial element abundances. In this paper, we detail the physics of cosmological light relics, and describe how measurements of their relic density and mass serve as probes of physics beyond the Standard Model. A measurement of the light relic density at the precision of upcoming cosmological surveys will point the way toward new physics or severely constrain the range of viable extensions to the Standard Model. △ Less

Submitted 15 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021. 59 pages, 12 figures

Scattering of Lyα Photons through the Reionizing Intergalactic Medium: I. Spectral Energy Distribution

Authors: Hyunbae Park, Hyo Jeong Kim, Kyungjin Ahn, Hyunmi Song, Intae Jung, Pierre Ocvirk, Paul R. Shapiro, Taha Dawoodbhoy, Jenny G. Sorce, Ilian T. Iliev

Abstract: During reionization, a fraction of galactic Ly$α$ emission is scattered in the intergalactic medium (IGM) and appears as a diffuse light extending megaparsecs from the source. We investigate how to probe the properties of the early galaxies and their surrounding IGM using this scattered light. We create a Monte Carlo algorithm to track individual photons and reproduce several test cases from previ… ▽ More During reionization, a fraction of galactic Ly$α$ emission is scattered in the intergalactic medium (IGM) and appears as a diffuse light extending megaparsecs from the source. We investigate how to probe the properties of the early galaxies and their surrounding IGM using this scattered light. We create a Monte Carlo algorithm to track individual photons and reproduce several test cases from previous literature. Then, we run our code on the simulated IGM of the CoDaII simulation. We find that the scattered light can leave an observable imprint on the emergent spectrum if collected over several square arcminutes. Scattering can redden the emission by increasing the path lengths of photons, but it can also make the photons bluer by upscattering them according to the peculiar motion of the scatterer. The photons emitted on the far blue side of the resonance appear more extended in both frequency and space compared to those emitted near the resonance. This provides a discriminating feature for the blueward emission, which cannot be constrained from the unscattered light coming directly from the source. The ionization state of the IGM also affects the scattered light spectrum. When the source is in a small HII region, the emission goes through more scatterings in the surrounding HI region regardless of the initial frequency and ends up more redshifted and spatially extended. This can result in a weakening of the scattered light toward high $z$ during reionization. Our results provide a framework for interpreting the scattered light to be measured by high-$z$ integral-field-unit surveys. △ Less

Submitted 16 May, 2022; v1 submitted 13 February, 2022; originally announced February 2022.

Comments: Accepted for publication in the ApJ

The short ionizing photon mean free path at z=6 in Cosmic Dawn III, a new fully-coupled radiation-hydrodynamical simulation of the Epoch of Reionization

Authors: Joseph S. W. Lewis, Pierre Ocvirk, Jenny G. Sorce, Yohan Dubois, Dominique Aubert, Luke Conaboy, Paul R. Shapiro, Taha Dawoodbhoy, Romain Teyssier, Gustavo Yepes, Stefan Gottlöber, Yann Rasera, Kyungjin Ahn, Ilian T. Iliev, Hyunbae Park, Émilie Thélie

Abstract: Recent determinations of the mean free path of ionising photons (mfp) in the intergalactic medium (IGM) at $\rm z=6$ are lower than many theoretical predictions. To gain insight into this issue, we investigate the evolution of the mfp in our new massive fully coupled radiation hydrodynamics cosmological simulation of reionization: Cosmic Dawn III (CoDaIII). CoDaIII's scale ($\rm 94^3 \, cMpc^3$) a… ▽ More Recent determinations of the mean free path of ionising photons (mfp) in the intergalactic medium (IGM) at $\rm z=6$ are lower than many theoretical predictions. To gain insight into this issue, we investigate the evolution of the mfp in our new massive fully coupled radiation hydrodynamics cosmological simulation of reionization: Cosmic Dawn III (CoDaIII). CoDaIII's scale ($\rm 94^3 \, cMpc^3$) and resolution ($\rm 8192^3$ grid) make it particularly suitable to study the evolution of the IGM during Reionization. The simulation was performed with RAMSES-CUDATON on Summit, and used 131072 processors coupled to 24576 GPUs, making it the largest Reionization simulation, and largest RAMSES simulation ever performed. A superior agreement with global constraints on Reionization is obtained in CoDaIII over CoDaII especially for the evolution of the neutral hydrogen fraction and the cosmic photo-ionization rate, thanks to an improved calibration, later end of reionization ($\rm z=5.6$), and higher spatial resolution. Analyzing the mfp, we find that CoDaIII reproduces the most recent observations very well, from $\rm z=6$ to $\rm z=4.6$. We show that the distribution of the mfp in CoDaIII is bimodal, with short (neutral) and long (ionized) mfp modes, respectively, due to the patchiness of reionization and the co-existence of neutral versus ionized regions during Reionization. The neutral mode peaks at sub-kpc to kpc scales of mfp, while the ionized mode peak evolves from $\rm 0.1 Mpc/h$ at $\rm z=7$ to $\sim 10$ Mpc/h at $\rm z=5.2$. Computing the mfp as the average of the ionized mode provides the best match to the recent observational determinations. The distribution reduces to a single neutral (ionized) mode at $\rm z>13$ ($\rm z<5$). △ Less

Submitted 22 August, 2022; v1 submitted 11 February, 2022; originally announced February 2022.

Comments: Submitted to MNRAS Accepted 2022 August 08. Received 2022 August 05; in original form 2022 February 11

Precision cosmology and the stiff-amplified gravitational-wave background from inflation: NANOGrav, Advanced LIGO-Virgo and the Hubble tension

Authors: Bohua Li, Paul R. Shapiro

Abstract: The recent NANOGrav finding of a common-spectrum process has invited interpretations as possible evidence of a primordial stochastic gravitational-wave background (SGWB) stronger than predicted by standard inflation+LCDM. Such an SGWB would contribute an extra radiation component to the background Universe which may affect its expansion history. As such, it may help alleviate the current Hubble te… ▽ More The recent NANOGrav finding of a common-spectrum process has invited interpretations as possible evidence of a primordial stochastic gravitational-wave background (SGWB) stronger than predicted by standard inflation+LCDM. Such an SGWB would contribute an extra radiation component to the background Universe which may affect its expansion history. As such, it may help alleviate the current Hubble tension, a novel connection between gravitational waves and cosmology. We demonstrate this by considering a cosmological model, the "standard inflation + stiff amplification" scenario, with two components added to the LCDM model: a stiff component (w=1) and the primordial SGWB. Previously, we showed that even for standard inflation, the SGWB may be detectable at the high frequencies probed by laser interferometers, if it is amplified by a possible early stiff era after reheating. Models that boost the SGWB enough to cause significant backreaction, however, must still preserve the well-measured radiation-matter equality, as precision cosmology demands. For that, we calculate the fully-coupled evolution of the SGWB and expansion history, sampling parameter space (tensor-to-scalar ratio, reheating temperature and temperature at stiff-to-radiation equality). We then perform a joint analysis of the NANOGrav results and latest upper bounds from Planck, big bang nucleosynthesis and Advanced LIGO-Virgo, to constrain the model. The resulting blue-tilted, stiff-amplified SGWB is still too small to explain the NANOGrav results. However, if someday, Advanced LIGO-Virgo detects the SGWB, our model can explain it within standard inflation (without requiring an initial blue tilt). Meanwhile, this model may bring current high-z measurements of the Hubble constant within 3.4 sigma of the low-z measurements by SH0ES (from 4.4 sigma) and within 2.6 sigma of those by H0LiCOW (from 3.1 sigma), reducing the tension. △ Less

Submitted 13 October, 2021; v1 submitted 26 July, 2021; originally announced July 2021.

Comments: 31 pages, 11 figures. Published in JCAP. (This version) Updated to match the published version

Journal ref: JCAP10 (2021) 024

Cosmological structure formation in scalar field dark matter with repulsive self-interaction: The Incredible Shrinking Jeans Mass

Authors: Paul R. Shapiro, Taha Dawoodbhoy, Tanja Rindler-Daller

Abstract: Scalar Field Dark Matter (SFDM) comprised of ultralight ($\gtrsim 10^{-22}$ eV) bosons is an alternative to standard, collisionless Cold Dark Matter (CDM) that is CDM-like on large scales but inhibits small-scale structure formation. As a Bose-Einstein condensate, its free-field ("fuzzy") limit (FDM) suppresses structure below the de Broglie wavelength, $λ_\text{deB}$, creating virialized haloes w… ▽ More Scalar Field Dark Matter (SFDM) comprised of ultralight ($\gtrsim 10^{-22}$ eV) bosons is an alternative to standard, collisionless Cold Dark Matter (CDM) that is CDM-like on large scales but inhibits small-scale structure formation. As a Bose-Einstein condensate, its free-field ("fuzzy") limit (FDM) suppresses structure below the de Broglie wavelength, $λ_\text{deB}$, creating virialized haloes with central cores of radius $\simλ_\text{deB}$, surrounded by CDM-like envelopes, and a halo mass function (HMF) with a sharp cut-off on small scales. With a strong enough repulsive self-interaction (SI), structure is inhibited, instead, below the Thomas-Fermi (TF) radius, $R_\text{TF}$ (the size of an SI-pressure-supported ($n=1$)-polytrope), when $R_\text{TF} > λ_\text{deB}$. Previously, we developed tools to describe SFDM dynamics on scales above $λ_\text{deB}$ and showed that SFDM-TF haloes formed by Jeans-unstable collapse from non-cosmological initial conditions have $R_\text{TF}$-sized cores, surrounded by CDM-like envelopes. Revisiting SFDM-TF in the cosmological context, we simulate halo formation by cosmological infall and collapse, and derive its transfer function from linear perturbation theory to produce cosmological initial conditions and predict statistical measures of structure formation, such as the HMF. Since FDM and SFDM-TF transfer functions both have small-scale cut-offs, we can align them to let observational constraints on FDM proxy for SFDM-TF, finding FDM with particle masses $1 \lesssim m/(10^{-22} \text{ eV}/c^2) \lesssim 30$ corresponds to SFDM-TF with $10 \gtrsim R_\text{TF}/(1 \text{ pc}) \gtrsim 1$, favoring sub-galactic (sub-kpc) core-size. The SFDM-TF HMF cuts off gradually, however, leaving more small-mass haloes: its Jeans mass shrinks so fast, scales filtered early can still recover and grow! △ Less

Submitted 3 November, 2021; v1 submitted 24 June, 2021; originally announced June 2021.

Comments: 30 pages, 10 figures, accepted by MNRAS 09/22/21, comments welcome

Crucial Factors for Lyα Transmission in the Reionizing Intergalactic Medium: Infall Motion, HII Bubble Size, and Self-shielded Systems

Authors: Hyunbae Park, Intae Jung, Hyunmi Song, Pierre Ocvirk, Paul R. Shapiro, Taha Dawoodbhoy, Ilian T. Iliev, Kyungjin Ahn, Michele Bianco, Hyo Jeong Kim

Abstract: Using the CoDa II simulation, we study the Ly$α$ transmissivity of the intergalactic medium (IGM) during reionization. At $z>6$, a typical galaxy without an active galactic nucleus fails to form a proximity zone around itself due to the overdensity of the surrounding IGM. The gravitational infall motion in the IGM makes the resonance absorption extend to the red side of Ly$α$, suppressing the tran… ▽ More Using the CoDa II simulation, we study the Ly$α$ transmissivity of the intergalactic medium (IGM) during reionization. At $z>6$, a typical galaxy without an active galactic nucleus fails to form a proximity zone around itself due to the overdensity of the surrounding IGM. The gravitational infall motion in the IGM makes the resonance absorption extend to the red side of Ly$α$, suppressing the transmission up to roughly the circular velocity of the galaxy. In some sight lines, an optically thin blob generated by a supernova in a neighboring galaxy results in a peak feature, which can be mistaken for a blue peak. Redward of the resonance absorption, the damping-wing opacity correlates with the global IGM neutral fraction and the UV magnitude of the source galaxy. Brighter galaxies tend to suffer lower opacity because they tend to reside in larger HII regions, and the surrounding IGM transmits redder photons, which are less susceptible to attenuation, owing to stronger infall velocity. The HII regions are highly nonspherical, causing both sight-line-to-sight-line and galaxy-to-galaxy variation in opacity. Also, self-shielded systems within HII regions strongly attenuate the emission for certain sight lines. All these factors add to the transmissivity variation, requiring a large sample size to constrain the average transmission. The variation is largest for fainter galaxies at higher redshift. The 68\% range of the transmissivity is similar to or greater than the median for galaxies with $M_{\rm UV}\ge-21$ at $z\ge7$, implying that more than a hundred galaxies would be needed to measure the transmission to 10\% accuracy. △ Less

Submitted 6 December, 2021; v1 submitted 22 May, 2021; originally announced May 2021.

Comments: 18 pages, 16 figures, 4 tables, Published in ApJ

Core-Envelope Haloes in Scalar Field Dark Matter with Repulsive Self-Interaction: Fluid Dynamics Beyond the de Broglie Wavelength

Authors: Taha Dawoodbhoy, Paul R. Shapiro, Tanja Rindler-Daller

Abstract: Scalar Field Dark Matter (SFDM) comprised of ultralight bosons has attracted great interest as an alternative to standard, collisionless Cold Dark Matter (CDM) because of its novel structure-formation dynamics, described by the coupled Schrödinger-Poisson equations. In the free-field ("fuzzy") limit of SFDM (FDM), structure is inhibited below the de Broglie wavelength, but resembles CDM on larger… ▽ More Scalar Field Dark Matter (SFDM) comprised of ultralight bosons has attracted great interest as an alternative to standard, collisionless Cold Dark Matter (CDM) because of its novel structure-formation dynamics, described by the coupled Schrödinger-Poisson equations. In the free-field ("fuzzy") limit of SFDM (FDM), structure is inhibited below the de Broglie wavelength, but resembles CDM on larger scales. Virialized haloes have "solitonic" cores of radius $\simλ_\text{deB}$, surrounded by CDM-like envelopes. When a strong enough repulsive self-interaction (SI) is also present, structure can be inhibited below a second length scale, $λ_\text{SI}$, with $λ_\text{SI}> λ_\text{deB}$ -- called the Thomas-Fermi (TF) regime. FDM dynamics differs from CDM because of quantum pressure, and SFDM-TF differs further by adding SI pressure. In the small-$λ_\text{deB}$ limit, however, we can model all three by fluid conservation equations for a compressible, $γ=5/3$ ideal gas, with ideal gas pressure sourced by internal velocity dispersion and, for the TF regime, an added SI pressure, $P_\text{SI}\propto ρ^2$. We use these fluid equations to simulate halo formation from gravitational collapse in 1D, spherical symmetry, demonstrating for the first time that SFDM-TF haloes form with cores the size of $R_\text{TF}$, the radius of an SI-pressure-supported $(n=1)$-polytrope, surrounded by CDM-like envelopes. In comparison with rotation curves of dwarf galaxies in the local Universe, SFDM-TF haloes pass the ["too-big-to-fail" + "cusp-core"]-test if $R_\text{TF}\gtrsim 1$ kpc. △ Less

Submitted 31 July, 2021; v1 submitted 14 April, 2021; originally announced April 2021.

Comments: 27 pages, 11 figures, accepted by MNRAS 06/28/21, comments welcome

Journal ref: MNRAS, 506, 2418-2444 (2021)

Angular Momentum and the Absence of Vortices in the Cores of Fuzzy Dark Matter Haloes

Authors: Sonja O. Schobesberger, Tanja Rindler-Daller, Paul R. Shapiro

Abstract: Scalar Field Dark Matter (SFDM), comprised of ultralight ($\gtrsim 10^{-22}$ eV) bosons, is distinguished from massive ($\gtrsim$ GeV), collisionless Cold Dark Matter (CDM) by its novel structure-formation dynamics as Bose-Einstein condensate (BEC) and quantum superfluid with wave-like properties, described by the Gross-Pitaevski and Poisson (GPP) equations. In the free-field (fuzzy) limit of SFDM… ▽ More Scalar Field Dark Matter (SFDM), comprised of ultralight ($\gtrsim 10^{-22}$ eV) bosons, is distinguished from massive ($\gtrsim$ GeV), collisionless Cold Dark Matter (CDM) by its novel structure-formation dynamics as Bose-Einstein condensate (BEC) and quantum superfluid with wave-like properties, described by the Gross-Pitaevski and Poisson (GPP) equations. In the free-field (fuzzy) limit of SFDM (FDM), structure is inhibited below the de Broglie wavelength $λ_{\text{deB}}$, but resembles CDM on larger scales. Virialized haloes have solitonic cores of radius $\sim λ_{\text{deB}}$ that follow the ground-state attractor solution of GPP, surrounded by CDM-like envelopes. As superfluid, SFDM is irrotational (vorticity-free) but can be unstable to vortex formation. We previously showed this can happen in halo cores, from angular momentum arising during structure formation, when repulsive self-interaction (SI) is present to support them out to a second length scale $λ_{\text{SI}}$ with $λ_{\text{SI}} > λ_{\text{deB}}$ (the Thomas-Fermi regime), but only if SI is strong enough. This suggested FDM cores (without SI) would not form vortices. FDM simulations later found vortices, but only outside halo cores, consistent with our previous suggestion based upon TF-regime analysis. We extend that analysis now to FDM, to show explicitly that vortices should not arise in solitonic cores from angular momentum, modelling them as either Gaussian spheres or compressible, ($n = 2$)-polytropic, irrotational Riemann-S ellipsoids. We find that, for typical halo spin parameters, angular momentum per particle is below $\hbar$, the minimum required even for one singly-quantized vortex in the centre. Even for larger angular momentum, however, vortex formation is not energetically favoured. △ Less

Submitted 22 April, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Comments: accepted by MNRAS, 29 pages, 9 figures; minor revisions

The impact of inhomogeneous subgrid clumping on cosmic reionization II: modelling stochasticity

Authors: Michele Bianco, Ilian T. Iliev, Kyungjin Ahn, Sambit K. Giri, Yi Mao, Hyunbae Park, Paul R. Shapiro

Abstract: Small-scale density fluctuations can significantly affect reionization but are typically modelled quite crudely. Unresolved fluctuations in numerical simulations and analytical calculations are included using a gas clumping factor, typically assumed to be independent of the local environment. In Paper I, we presented an improved, local density-dependent model for the sub-grid gas clumping. Here we… ▽ More Small-scale density fluctuations can significantly affect reionization but are typically modelled quite crudely. Unresolved fluctuations in numerical simulations and analytical calculations are included using a gas clumping factor, typically assumed to be independent of the local environment. In Paper I, we presented an improved, local density-dependent model for the sub-grid gas clumping. Here we extend this using an empirical stochastic model based on the results from high-resolution numerical simulations which fully resolve all relevant fluctuations. Our model reproduces well both the mean density-clumping relation and its scatter. We applied our stochastic model, along with the mean clumping one and the Paper I deterministic model, to create a large-volume realisation of the clumping field, and used these in radiative transfer simulations of cosmic reionization. Our results show that the simplistic mean clumping model delays reionization compared to local density-dependent models, despite producing fewer recombinations overall. This is due to the very different spatial distribution of clumping, resulting in much higher photoionization rates in the latter cases. The mean clumping model produces smaller HII regions throughout most of the reionization, but those percolate faster at late times. It also causes a significant delay in the 21-cm fluctuations peak and yields lower non-Gaussianity and many fewer bright pixels in the PDF distribution. The stochastic density-dependent model shows relatively minor differences from the deterministic one, mostly concentrated around overlap, where it significantly suppresses the 21-cm fluctuations, and at the bright tail of the 21-cm PDFs, where it produces noticeably more bright pixels. △ Less

Submitted 5 May, 2021; v1 submitted 5 January, 2021; originally announced January 2021.

Comments: 18 pages, 17 figures Minor change: enlarged the inset table in Fig. 1

Journal ref: Monthly Notices of the Royal Astronomical Society: 24 March 2021

Cosmic Reionization May Still Have Started Early and Ended Late: Confronting Early Onset with CMB Anisotropy and 21 cm Global Signals

Authors: Kyungjin Ahn, Paul R. Shapiro

Abstract: The global history of reionization was shaped by the relative amounts of starlight released by three halo mass groups: atomic-cooling halos (ACHs) with virial temperatures Tvir > 10^4 K, either (1) massive enough to form stars even after reionization (HMACHs, >~ 10^9 Msun) or (2) less-massive (LMACHs), subject to star formation suppression when overtaken by reionization, and (3) H2-cooling minihal… ▽ More The global history of reionization was shaped by the relative amounts of starlight released by three halo mass groups: atomic-cooling halos (ACHs) with virial temperatures Tvir > 10^4 K, either (1) massive enough to form stars even after reionization (HMACHs, >~ 10^9 Msun) or (2) less-massive (LMACHs), subject to star formation suppression when overtaken by reionization, and (3) H2-cooling minihalos (MHs) with Tvir < 10^4 K, whose star formation is predominantly suppressed by the H2-dissociating Lyman-Werner (LW) background. Our previous work showed that including MHs caused two-stage reionization - early rise to x ~ 0.1, driven by MHs, followed by a rapid rise, late, to x ~ 1, driven by ACHs - with a signature in CMB polarization anisotropy predicted to be detectable by the Planck satellite. Motivated by this prediction, we model global reionization semi-analytically for comparison with Planck CMB data and the EDGES global 21cm absorption feature, for models with: (1) ACHs, no feedback; (2) ACHs, self-regulated; and (3) ACHs and MHs, self-regulated. Model (3) agrees well with Planck E-mode polarization data, even with a substantial tail of high-redshift ionization, beyond the limit proposed by the Planck Collaboration (2018). No model reproduces the EDGES feature. For model (3), |dTb| <~ 60 mK across the EDGES trough, an order of magnitude too shallow, and absorption starts at higher z but is spectrally featureless. Early onset reionization by Population III stars in MHs is compatible with current constraints, but only if the EDGES interpretation is discounted or else other processes we did not include account for it. △ Less

Submitted 9 April, 2021; v1 submitted 6 November, 2020; originally announced November 2020.

Comments: Version accepted to ApJ with minor mod; added references; new figures for H-ionizing UV background prediction

Core-Halo Mass Relation in Scalar Field Dark Matter Models and its Consequences for the Formation of Supermassive Black Holes

Authors: Luis E. Padilla, Tanja Rindler-Daller, Paul R. Shapiro, Tonatiuh Matos, J. Alberto Vázquez

Abstract: Scalar-field dark matter (SFDM) halos exhibit a core-envelope structure with soliton-like cores and CDM-like envelopes. Simulations without self-interaction (free-field case) report a core-halo mass relation $M_c\propto M_{h}^β$, with either $β=1/3$ or $β=5/9$, which can be understood if core and halo obey certain energy or velocity scalings. We extend the core-halo mass relations to include SFDM… ▽ More Scalar-field dark matter (SFDM) halos exhibit a core-envelope structure with soliton-like cores and CDM-like envelopes. Simulations without self-interaction (free-field case) report a core-halo mass relation $M_c\propto M_{h}^β$, with either $β=1/3$ or $β=5/9$, which can be understood if core and halo obey certain energy or velocity scalings. We extend the core-halo mass relations to include SFDM with self-interaction (SI), either repulsive or attractive, and investigate its implications for the gravitational instability and collapse of solitonic cores, leading to supermassive black hole (SMBH) formation. For SFDM parameters that make $\sim$ Kpc-sized cores and CDM-like structure formation on large scales but suppressed on small scales, cores are stable for all galactic halos of interest, from the free-field to the repulsive SI limit. For attractive SI, however, halos masses $M_h\sim (10^{10}-10^{12}) M_\odot$ have cores that collapse to SMBHs with $M_{SMBH}\sim 10^{6}-10^8 M_\odot$, as observations seem to require, while smaller-mass halos have stable cores, for particle masses $m\simeq (2.14\times 10^{-22}-9.9\times 10^{-20})\rm{eV}/c^2$, if the free-field has $β=1/3$, or $m = 2.23\times 10^{-21}-1.7\times 10^{-18}\rm{eV}/c^2$, if $β=5/9$. For free-field and repulsive cases, however, if previous constraints on particle parameters are relaxed to allow much smaller (sub-galactic scale) cores, then halos can also form SMBHs, for the same range of halo and BH masses, as long as $β= 5/9$ is correct for the free-field. In that case, structure formation in SFDM would be largely indistinguishable from Cold Dark Matter (CDM). Such SFDM models might not resolve the small-scale structure problems of CDM, but they would explain the formation of SMBHs quite naturally. Since CDM, itself, has not yet been ruled out, such SFDM models must also be viable (Abbreviated). △ Less

Submitted 19 April, 2021; v1 submitted 23 October, 2020; originally announced October 2020.

Comments: 27 pages, 7 figures, minor updates to accord with the published version in Phys. Rev. D

Journal ref: Phys. Rev. D 103, 063012 (2021)

Large-scale variation in reionization history caused by Baryon-dark matter streaming velocity

Authors: Hyunbae Park, Paul R. Shapiro, Kyungjin Ahn, Naoki Yoshida, Shingo Hirano

Abstract: At cosmic recombination, there was supersonic relative motion between baryons and dark matter, which originated from the baryonic acoustic oscillations in the early universe. This motion has been considered to have a negligible impact on the late stage of cosmic reionization because the relative velocity quickly decreases. However, recent studies have suggested that the recombination in gas clouds… ▽ More At cosmic recombination, there was supersonic relative motion between baryons and dark matter, which originated from the baryonic acoustic oscillations in the early universe. This motion has been considered to have a negligible impact on the late stage of cosmic reionization because the relative velocity quickly decreases. However, recent studies have suggested that the recombination in gas clouds smaller than the local Jeans mass ($\lesssim$ $10^8~M_\odot$) can affect the reionization history by boosting the number of ultraviolet photons required for ionizing the intergalactic medium. Motivated by this, we performed a series of radiation-hydrodynamic simulations to investigate whether the streaming motion can generate variation in the local reionization history by smoothing out clumpy small-scale structures and lowering the ionizing photon budget. We found that the streaming velocity can add a variation of $Δz_e$ $\sim$ $0.05$ $-$ $0.5$ in the end-of-reionization redshift, depending on the level of X-ray preheating and the time evolution of ionizing sources. The variation tends to be larger when the ionizing efficiency of galaxies decreases toward later times. Given the long spatial fluctuation scales of the streaming motion ($\gtrsim 100$ Mpc), it can help to explain the Ly$α$ opacity variation observed from quasars and leave large-scale imprints on the ionization field of the intergalactic medium during the reionization. The pre-reionization heating by X-ray sources is another critical factor that can suppress small-scale gas clumping and can diminish the variation in $z_e$ introduced by the streaming motion. △ Less

Submitted 22 January, 2021; v1 submitted 23 October, 2020; originally announced October 2020.

Comments: 12 pages, 8 figures

Lyman-alpha transmission properties of the intergalactic medium in the CoDaII simulation

Authors: Max Gronke, Pierre Ocvirk, Charlotte Mason, Jorryt Matthee, Sarah E. I. Bosman, Jenny G. Sorce, Joseph Lewis, Kyungjin Ahn, Dominique Aubert, Taha Dawoodbhoy, Ilian T. Iliev, Paul R. Shapiro, Gustavo Yepes

Abstract: The decline in abundance of Lyman-$α$ (Ly$α$) emitting galaxies at $z \gtrsim 6$ is a powerful and commonly used probe to constrain the progress of cosmic reionization. We use the CoDaII simulation, which is a radiation hydrodynamic simulation featuring a box of $\sim 94$ comoving Mpc side length, to compute the Ly$α$ transmission properties of the intergalactic medium (IGM) at $z\sim 5.8$ to $7$.… ▽ More The decline in abundance of Lyman-$α$ (Ly$α$) emitting galaxies at $z \gtrsim 6$ is a powerful and commonly used probe to constrain the progress of cosmic reionization. We use the CoDaII simulation, which is a radiation hydrodynamic simulation featuring a box of $\sim 94$ comoving Mpc side length, to compute the Ly$α$ transmission properties of the intergalactic medium (IGM) at $z\sim 5.8$ to $7$. Our results mainly confirm previous studies, i.e., we find a declining Ly$α$ transmission with redshift and a large sightline-to-sightline variation. However, motivated by the recent discovery of blue Ly$α$ peaks at high redshift, we also analyze the IGM transmission on the blue side, which shows a rapid decline at $z\gtrsim 6$ of the blue transmission. This low transmission can be attributed not only to the presence of neutral regions but also to the residual neutral hydrogen within ionized regions, for which a density even as low as $n_{\rm HI}\sim 10^{-9}\,\mathrm{cm}^{-3}$ (sometimes combined with kinematic effects) leads to a significantly reduced visibility. Still, we find that $\sim 1\%$ of sightlines towards $M_{\mathrm{1600AB}}\sim -21$ galaxies at $z\sim 7$ are transparent enough to allow a transmission of a blue Ly$α$ peak. We discuss our results in the context of the interpretation of observations. △ Less

Submitted 21 October, 2021; v1 submitted 29 April, 2020; originally announced April 2020.

Comments: 13 pages, 12 figures. Published version

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 508, Issue 3, December 2021, Pages 3697-3709

Galactic ionising photon budget during the Epoch of Reionisation in the Cosmic Dawn II simulation

Authors: Joseph S. W. Lewis, Pierre Ocvirk, Dominique Aubert, Jenny G. Sorce, Paul R. Shapiro, Nicolas Deparis, Taha Dawoodbhoy, Romain Teyssier, Gustavo Yepes, Stefan Gottlöber, Kyungjin Ahn, Ilian T. Iliev, Jonathan Chardin

Abstract: Cosmic Dawn ("CoDa") II yields the first statistically-meaningful determination of the relative contribution to reionization by galaxies of different halo mass, from a fully-coupled radiation-hydrodynamics simulation of the epoch of reionization large enough ($\sim$ 100 Mpc) to model global reionization while resolving the formation of all galactic halos above $\sim 10^8 M_\odot$. Cell transmissio… ▽ More Cosmic Dawn ("CoDa") II yields the first statistically-meaningful determination of the relative contribution to reionization by galaxies of different halo mass, from a fully-coupled radiation-hydrodynamics simulation of the epoch of reionization large enough ($\sim$ 100 Mpc) to model global reionization while resolving the formation of all galactic halos above $\sim 10^8 M_\odot$. Cell transmission inside high-mass haloes is bi-modal -- ionized cells are transparent, while neutral cells absorb the photons their stars produce - and the halo escape fraction $f_{esc}$ reflects the balance of star formation rate ("SFR") between these modes. The latter is increasingly prevalent at higher halo mass, driving down $f_{esc}$ (we provide analytical fits to our results), whereas halo escape luminosity, proportional to $f_{esc} \times$SFR, increases with mass. Haloes with dark matter masses within $6.10^{8} M_\odot < M_h < 3.10^{10} M_\odot$ produce $\sim 80$% of the escaping photons at z=7, when the Universe is 50% ionized, making them the main drivers of cosmic reionization. Less massive haloes, though more numerous, have low SFRs and contribute less than 10% of the photon budget then, despite their high $f_{esc}$. High mass haloes are too few and too opaque, contributing $<10$% despite their high SFRs. The dominant mass range is lower (higher) at higher (lower) redshift, as mass function and reionization advance together (e.g. at z$=8.5$, x$_{\rm HI}=0.9$, $M_h < 5.10^9 M_\odot$ haloes contributed $\sim$80%). Galaxies with UV magnitudes $M_{AB1600}$ between $-12$ and $-19$ dominated reionization between z$=6$ and 8. △ Less

Submitted 18 June, 2020; v1 submitted 21 January, 2020; originally announced January 2020.

Comments: Submitted to MNRAS on 09/01/2020. Accepted on 12/06/2020

The impact of inhomogeneous subgrid clumping on cosmic reionization

Authors: Yi Mao, Jun Koda, Paul R. Shapiro, Ilian T. Iliev, Garrelt Mellema, Hyunbae Park, Kyungjin Ahn, Michele Bianco

Abstract: Cosmic reionization was driven by the imbalance between early sources and sinks of ionizing radiation, both of which were dominated by small-scale structure and are thus usually treated in cosmological reionization simulations by subgrid modelling. The recombination rate of intergalactic hydrogen is customarily boosted by a subgrid clumping factor, ${\left<n^2\right>/\left<n\right>^2}$, which corr… ▽ More Cosmic reionization was driven by the imbalance between early sources and sinks of ionizing radiation, both of which were dominated by small-scale structure and are thus usually treated in cosmological reionization simulations by subgrid modelling. The recombination rate of intergalactic hydrogen is customarily boosted by a subgrid clumping factor, ${\left<n^2\right>/\left<n\right>^2}$, which corrects for unresolved fluctuations in gas density ${n}$ on scales below the grid-spacing of coarse-grained simulations. We investigate in detail the impact of this inhomogeneous subgrid clumping on reionization and its observables, as follows: (1) Previous attempts generally underestimated the clumping factor because of insufficient mass resolution. We perform a high-resolution $N$-body simulation that resolves haloes down to the pre-reionization Jeans mass to derive the time-dependent, spatially-varying local clumping factor and a fitting formula for its correlation with local overdensity. (2) We then perform a large-scale $N$-body and radiative transfer simulation that accounts for this inhomogeneous subgrid clumping by applying this clumping factor-overdensity correlation. Boosting recombination significantly slows the expansion of ionized regions, which delays completion of reionization and suppresses 21 cm power spectra on large scales in the later stages of reionization. (3) We also consider a simplified prescription in which the globally-averaged, time-evolving clumping factor from the same high-resolution $N$-body simulation is applied uniformly to all cells in the reionization simulation, instead. Observables computed with this model agree fairly well with those from the inhomogeneous clumping model, e.g. predicting 21 cm power spectra to within 20% error, suggesting it may be a useful approximation. △ Less

Submitted 12 November, 2019; v1 submitted 6 June, 2019; originally announced June 2019.

Comments: 23 pages, 16 figures, 4 tables. MNRAS accepted

Journal ref: MNRAS, Volume 491, Issue 2, p.1600-1621 (2020)

Cosmic Dawn II (CoDa II): a new radiation-hydrodynamics simulation of the self-consistent coupling of galaxy formation and reionization

Authors: Pierre Ocvirk, Dominique Aubert, Jenny G. Sorce, Paul R. Shapiro, Nicolas Deparis, Taha Dawoodbhoy, Joseph Lewis, Romain Teyssier, Gustavo Yepes, Stefan Gottlöber, Kyungjin Ahn, Ilian T. Iliev, Yehuda Hoffman

Abstract: Cosmic Dawn II (CoDa II) is a new, fully-coupled radiation-hydrodynamics simulation of cosmic reionization and galaxy formation and their mutual impact, to redshift $z < 6$. With $4096^3$ particles and cells in a 94 Mpc box, it is large enough to model global reionization and its feedback on galaxy formation while resolving all haloes above $10^8$ M$_{\odot}$. Using the same hybrid CPU-GPU code RA… ▽ More Cosmic Dawn II (CoDa II) is a new, fully-coupled radiation-hydrodynamics simulation of cosmic reionization and galaxy formation and their mutual impact, to redshift $z < 6$. With $4096^3$ particles and cells in a 94 Mpc box, it is large enough to model global reionization and its feedback on galaxy formation while resolving all haloes above $10^8$ M$_{\odot}$. Using the same hybrid CPU-GPU code RAMSES-CUDATON as CoDa I in Ocvirk et al. (2016), CoDa II modified and re-calibrated the subgrid star-formation algorithm, making reionization end earlier, at $z \gtrsim 6$, thereby better matching the observations of intergalactic Lyman-alpha opacity from quasar spectra and electron-scattering optical depth from cosmic microwave background fluctuations. CoDa II predicts a UV continuum luminosity function in good agreement with observations of high-z galaxies, especially at $z = 6$. As in CoDa I, reionization feedback suppresses star formation in haloes below $\sim 2 \times 10^9$ M$_{\odot}$, though suppression here is less severe, a possible consequence of modifying the star-formation algorithm. Suppression is environment-dependent, occurring earlier (later) in overdense (underdense) regions, in response to their local reionization times. Using a constrained realization of $Λ$CDM constructed from galaxy survey data to reproduce the large-scale structure and major objects of the present-day Local Universe, CoDa II serves to model both global and local reionization. In CoDa II, the Milky Way and M31 appear as individual islands of reionization, i.e. they were not reionized by the progenitor of the Virgo cluster, nor by nearby groups, nor by each other. △ Less

Submitted 14 October, 2020; v1 submitted 27 November, 2018; originally announced November 2018.

Comments: 20 pages, 14 figures, accepted for publication in MNRAS

Suppression of Star Formation in Low-Mass Galaxies Caused by the Reionization of their Local Neighborhood

Authors: Taha Dawoodbhoy, Paul R. Shapiro, Pierre Ocvirk, Dominique Aubert, Nicolas Gillet, Jun-Hwan Choi, Ilian T. Iliev, Romain Teyssier, Gustavo Yepes, Stefan Gottlöber, Anson D'Aloisio, Hyunbae Park, Yehuda Hoffman

Abstract: Photoheating associated with reionization suppressed star formation in low-mass galaxies. Reionization was inhomogeneous, however, affecting different regions at different times. To establish the causal connection between reionization and suppression, we must take this local variation into account. We analyze the results of CoDa (`Cosmic Dawn') I, the first fully-coupled radiation-hydrodynamical s… ▽ More Photoheating associated with reionization suppressed star formation in low-mass galaxies. Reionization was inhomogeneous, however, affecting different regions at different times. To establish the causal connection between reionization and suppression, we must take this local variation into account. We analyze the results of CoDa (`Cosmic Dawn') I, the first fully-coupled radiation-hydrodynamical simulation of reionization and galaxy formation in the Local Universe, in a volume large enough to model reionization globally but with enough resolving power to follow all atomic-cooling galactic halos in that volume. For every halo identified at a given time, we find the redshift at which the surrounding IGM reionized, along with its instantaneous star formation rate (`SFR') and baryonic gas-to-dark matter ratio ($M_\text{gas}/M_\text{DM}$). The average SFR per halo with $M < 10^9 \text{ M}_\odot$ was steady in regions not yet reionized, but declined sharply following local reionization. For $M > 10^{10} \text{ M}_\odot$, this SFR continued through local reionization, increasing with time, instead. For $10^9 < M < 10^{10} \text{ M}_\odot$, the SFR generally increased modestly through reionization, followed by a modest decline. In general, halo SFRs were higher for regions that reionized earlier. A similar pattern was found for $M_\text{gas}/M_\text{DM}$, which declined sharply following local reionization for $M < 10^9 \text{ M}_\odot$. Local reionization time correlates with local matter overdensity, which determines the local rates of structure formation and ionizing photon consumption. The earliest patches to develop structure and reionize ultimately produced more stars than they needed to finish and maintain their own reionization, exporting their `surplus' starlight to help reionize regions that developed structure later. △ Less

Submitted 25 July, 2018; v1 submitted 14 May, 2018; originally announced May 2018.

Comments: 15 pages, 15 figures, accepted by MNRAS 07/17/18, comments welcome

Journal ref: MNRAS, 480, 1740-1753 (2018)

The Inhomogeneous Reionization Times of Present-day Galaxies

Authors: Dominique Aubert, Nicolas Deparis, Pierre Ocvirk, Paul R. Shapiro, Ilian T. Iliev, Gustavo Yepes, Stefan Gottloeber, Yehuda Hoffman, Romain Teyssier

Abstract: Today's galaxies experienced cosmic reionization at different times in different locations. For the first time, reionization ($50\%$ ionized) redshifts, $z_R$, at the location of their progenitors are derived from new, fully-coupled radiation-hydrodynamics simulation of galaxy formation and reionization at $z > 6$, matched to N-body simulation to z = 0. Constrained initial conditions were chosen t… ▽ More Today's galaxies experienced cosmic reionization at different times in different locations. For the first time, reionization ($50\%$ ionized) redshifts, $z_R$, at the location of their progenitors are derived from new, fully-coupled radiation-hydrodynamics simulation of galaxy formation and reionization at $z > 6$, matched to N-body simulation to z = 0. Constrained initial conditions were chosen to form the well-known structures of the local universe, including the Local Group and Virgo, in a (91 Mpc)$^3$ volume large enough to model both global and local reionization. Reionization simulation CoDa I-AMR, by CPU-GPU code EMMA, used (2048)$^3$ particles and (2048)$^3$ initial cells, adaptively-refined, while N-body simulation CoDa I-DM2048, by Gadget2, used (2048)$^3$ particles, to find reionization times for all galaxies at z = 0 with masses $M(z=0)\ge 10^8 M_\odot$. Galaxies with $M(z=0) \gtrsim 10^{11} M_\odot$ reionized earlier than the universe as a whole, by up to $\sim$ 500 Myrs, with significant scatter. For Milky-Way-like galaxies, $z_R$ ranged from 8 to 15. Galaxies with $M(z=0) \lesssim 10^{11} M_\odot$ typically reionized as late or later than globally-averaged $50\%$ reionization at $\langle z_R\rangle =7.8$, in neighborhoods where reionization was completed by external radiation. The spread of reionization times within galaxies was sometimes as large as the galaxy-to-galaxy scatter. The Milky Way and M31 reionized earlier than global reionization but later than typical for their mass, neither dominated by external radiation. Their most massive progenitors at $z>6$ had $z_R$ = 9.8 (MW) and 11 (M31), while their total masses had $z_R$ = 8.2 (both). △ Less

Submitted 20 February, 2018; v1 submitted 5 February, 2018; originally announced February 2018.

Comments: 8 pages, 5 figures, updated version accepted for publication in ApJL

Bose-Einstein-condensed scalar field dark matter and the gravitational wave background from inflation: new cosmological constraints and its detectability by LIGO

Authors: Bohua Li, Paul R. Shapiro, Tanja Rindler-Daller

Abstract: We consider an alternative cold dark matter candidate, ultralight bosons ($m>10^{-22}$eV) described by a complex scalar field (SFDM) with global U(1) symmetry, with comoving particle number density conserved after particle production during standard reheating. We allow for repulsive self-interaction. In a Lambda-SFDM universe, SFDM starts relativistic, evolving from stiff (w=1) to radiation-like (… ▽ More We consider an alternative cold dark matter candidate, ultralight bosons ($m>10^{-22}$eV) described by a complex scalar field (SFDM) with global U(1) symmetry, with comoving particle number density conserved after particle production during standard reheating. We allow for repulsive self-interaction. In a Lambda-SFDM universe, SFDM starts relativistic, evolving from stiff (w=1) to radiation-like (w=1/3), becoming nonrelativistic (w=0) at late times. Thus, a stiff-SFDM-dominated era precedes the familiar radiation-dominated era. SFDM particle mass $m$ and quartic self-interaction strength λ, are therefore constrained by cosmological observables, N_{eff}, the effective number of neutrino species during BBN, and z_{eq}, the matter-radiation equality redshift. Since the stochastic gravitational wave background (SGWB) from inflation is amplified during the stiff-SFDM-dominated era, it can also contribute a radiationlike component large enough to affect these observables. Remarkably, this amplification makes this SGWB detectable by current GW experiments, e.g., aLIGO/Virgo and LISA, for Lambda-SFDM models satisfying cosmological constraints, for a range of reheat temperatures T_{re} and currently allowed values of tensor-to-scalar ratio $r$. For given r and $λ/(mc^2)^2$, the marginally-allowed Lambda-SFDM model for each T_{re} has the smallest m that satisfies cosmological constraints. For example, for marginally-allowed models with r=0.01 and $λ/(mc^2)^2=10^{-18}$eV$^{-1}$cm$^3$, null detection by the aLIGO O1 run excludes 8.75*10^3<T_{re} (GeV)<1.7*10^5 at 95% confidence, demonstrating that GW experiments already place a new kind of cosmological constraint on SFDM. A wider parameter range should be accessible to aLIGO/Virgo O5, with potential to detect this signature of Lambda-SFDM. For this same illustrative family, 3-sigma detection is predicted for 600<T_{re} (GeV)<10^7. △ Less

Submitted 25 April, 2017; v1 submitted 23 November, 2016; originally announced November 2016.

Comments: 54 pages, 13 figures, 2 tables, revised and resubmitted to Physical Review D. This version expanded the manuscript to take into account the implications on the SGWB from the recent release of the Advanced LIGO O1 data run. Typos fixed. References added

Journal ref: Phys. Rev. D 96, 063505 (2017)

On the Contribution of Active Galactic Nuclei to the High-Redshift Metagalactic Ionizing Background

Authors: Anson D'Aloisio, Phoebe R. Upton Sanderbeck, Matthew McQuinn, Hy Trac, Paul R. Shapiro

Abstract: Motivated by the claimed detection of a large population of faint active galactic nuclei (AGN) at high redshift, recent studies have proposed models in which AGN contribute significantly to the z > 4 H I ionizing background. In some models, AGN are even the chief sources of reionization. If correct, these models would make necessary a complete revision to the standard view that galaxies dominated… ▽ More Motivated by the claimed detection of a large population of faint active galactic nuclei (AGN) at high redshift, recent studies have proposed models in which AGN contribute significantly to the z > 4 H I ionizing background. In some models, AGN are even the chief sources of reionization. If correct, these models would make necessary a complete revision to the standard view that galaxies dominated the high-redshift ionizing background. It has been suggested that AGN-dominated models can better account for two recent observations that appear to be in conflict with the standard view: (1) large opacity variations in the z ~ 5.5 H I Lyman-alpha forest, and (2) slow evolution in the mean opacity of the He II Lyman-alpha forest. Large spatial fluctuations in the ionizing background from the brightness and rarity of AGN may account for the former, while the earlier onset of He II reionization in these models may account for the latter. Here we show that models in which AGN emissions source >~ 50 % of the ionizing background generally provide a better fit to the observed H I Lyman-alpha forest opacity variations compared to standard galaxy-dominated models. However, we argue that these AGN-dominated models are in tension with constraints on the thermal history of the intergalactic medium (IGM). Under standard assumptions about the spectra of AGN, we show that the earlier onset of He II reionization heats up the IGM well above recent temperature measurements. We further argue that the slower evolution of the mean opacity of the He II Lyman-alpha forest relative to simulations may reflect deficiencies in current simulations rather than favor AGN-dominated models as has been suggested. △ Less

Submitted 26 February, 2018; v1 submitted 21 July, 2016; originally announced July 2016.

Comments: Matches version published by MNRAS with clarifications and expanded discussion

Journal ref: Monthly Notices of the Royal Astronomical Society, 2017, Volume 468, p.4691-4701

The Hydrodynamic Feedback of Cosmic Reionization on Small-Scale Structures and Its Impact on Photon Consumption during the Epoch of Reionization

Authors: Hyunbae Park, Paul R. Shapiro, Jun-hwan Choi, Naoki Yoshida, Shingo Hirano, Kyungjin Ahn

Abstract: Density inhomogeneity in the intergalactic medium (IGM) can boost the recombination rate of ionized gas substantially, affecting the growth of HII regions during reionization. Previous attempts to quantify this effect typically failed to resolve down to the Jeans scale in the pre-ionization IGM, which is important in establishing this effect, along with the hydrodynamical back-reaction of reioniza… ▽ More Density inhomogeneity in the intergalactic medium (IGM) can boost the recombination rate of ionized gas substantially, affecting the growth of HII regions during reionization. Previous attempts to quantify this effect typically failed to resolve down to the Jeans scale in the pre-ionization IGM, which is important in establishing this effect, along with the hydrodynamical back-reaction of reionization on it. Towards that end, we perform a set of fully-coupled, radiation-hydrodynamics simulations from cosmological initial conditions, extending the mass resolution of previous work to the scale of minihalos. Pre-reionization structure is evolved until a redshift $z_i$ at which the ionizing radiation from external sources arrives to sweep an R-type ionization front supersonically across the volume in a few Myr, until it is trapped on the surfaces of minihalos and converted to D-type, after which the minihalo gas is removed by photoevaporative winds. Small-scale density structures during this time lead to a high ($>$10) clumping factor for ionized gas, which hugely boosts the recombination rate until the structures are disrupted by the hydrodynamic feedback after $\sim 10-100~\rm{Myr}$. For incoming stellar radiation with intensity $J_{21}$ in a $200~h^{-1}~\rm{kpc}$ box with the mean density contrast $\barδ$, the number of extra recombinations per H atom, on top of what is expected from homogeneously distributed gas, is given by $0.32[J_{21}]^{0.12}[(1+z_i)/11]^{-1.7}[1+\barδ]^{2.5}$. In models in which most of the volume is ionized toward the end of reionization, this can add more than one recombination per H atom to the ionizing photon budget to achieve reionization. △ Less

Submitted 17 October, 2016; v1 submitted 20 February, 2016; originally announced February 2016.

Comments: 21 pages, 19 figure, Accepted to ApJ

The Effects of Lyman-Limit Systems on the Evolution and Observability of the Epoch of Reionization

Authors: Hemant Shukla, Garrelt Mellema, Ilian T. Iliev, Paul R. Shapiro

Abstract: We present the first large-scale, full radiative transfer simulations of the reionization of the intergalactic medium in the presence of Lyman-limit systems (LLSs). To illustrate the impact of LLS opacity, possibly missed by previous simulations, we add either a uniform or spatially-varying hydrogen bound-free opacity. This opacity, implemented as the mean free path (mfp) of the ionizing photons,… ▽ More We present the first large-scale, full radiative transfer simulations of the reionization of the intergalactic medium in the presence of Lyman-limit systems (LLSs). To illustrate the impact of LLS opacity, possibly missed by previous simulations, we add either a uniform or spatially-varying hydrogen bound-free opacity. This opacity, implemented as the mean free path (mfp) of the ionizing photons, extrapolates the observed, post-reionization redshift dependence into the epoch of reionization. In qualitative agreement with previous studies, we find that at late times the presence of LLSs slows down the ionization fronts, and alters the size distribution of H II regions. We quantitatively characterize the size distribution and morphological evolution of H II regions and examine the effects of the LLSs on the redshifted 21-cm signal from the patchy reionization. The presence of LLSs extends the ionization history by $Δz \sim 0.8$. The LLS absorbers significantly impede the late-time growth of the H II regions. The position dependent LLS distribution slows reionization further and additionally limits the late growth of the ionized regions. However, there is no "freeze out" of the H II regions and the largest regions grow to the size of the simulation volume. The 21-cm power spectra show that at large scales the power drops by a factor of 2 for 50% and 75% ionization stages (at $k = 0.1$ $\text{h} \, \text{Mpc}^{-1} $) reflecting the limiting effect of the LLSs on the growth of ionized patches. The statistical observables such as the RMS of the brightness temperature fluctuations and the peak amplitudes of the 21-cm power spectra at large-scales ($k = 0.05 - 0.1$ $\text{h} \, \text{Mpc}^{-1} $) are diminished by the presence of LLS. △ Less

Submitted 2 February, 2016; originally announced February 2016.

The Large-Scale Observational Signatures of Low-Mass Galaxies During Reionization

Authors: Keri L. Dixon, Ilian T. Iliev, Garrelt Mellema, Kyungjin Ahn, Paul R. Shapiro

Abstract: Observations of the epoch of reionization give us clues about the nature and evolution of the sources of ionizing photons, or early stars and galaxies. We present a new suite of structure formation and radiative transfer simulations from the PRACE4LOFAR project designed to investigate whether the mechanism of radiative feedback, or the suppression of star formation in ionized regions from UV radia… ▽ More Observations of the epoch of reionization give us clues about the nature and evolution of the sources of ionizing photons, or early stars and galaxies. We present a new suite of structure formation and radiative transfer simulations from the PRACE4LOFAR project designed to investigate whether the mechanism of radiative feedback, or the suppression of star formation in ionized regions from UV radiation, can be inferred from these observations. Our source halo mass extends down to $10^8 M_\odot$, with sources in the mass range $10^8$ to $10^9 M_\odot$ expected to be particularly susceptible to feedback from ionizing radiation, and we vary the aggressiveness and nature of this suppression. Not only do we have four distinct source models, we also include two box sizes (67 Mpc and 349 Mpc), each with two grid resolutions. This suite of simulations allows us to investigate the robustness of our results. All of our simulations are broadly consistent with the observed electron-scattering optical depth of the cosmic microwave background and the neutral fraction and photoionization rate of hydrogen at $z\sim6$. In particular, we investigate the redshifted 21-cm emission in anticipation of upcoming radio interferometer observations. We find that the overall shape of the 21-cm signal and various statistics are robust to the exact nature of source suppression, the box size, and the resolution. There are some promising model discriminators in the non-Gaussianity and small-scale power spectrum of the 21-cm signal. △ Less

Submitted 11 December, 2015; originally announced December 2015.

Comments: 21 pages, 15 figures, MNRAS accepted

Cosmic Dawn (CoDa): the First Radiation-Hydrodynamics Simulation of Reionization and Galaxy Formation in the Local Universe

Authors: Pierre Ocvirk, Nicolas Gillet, Paul R. Shapiro, Dominique Aubert, Ilian T. Iliev, Romain Teyssier, Gustavo Yepes, Jun-Hwan Choi, David Sullivan, Alexander Knebe, Stefan Gottloeber, Anson D'Aloisio, Hyunbae Park, Yehuda Hoffman, Timothy Stranex

Abstract: Cosmic reionization by starlight from early galaxies affected their evolution, thereby impacting reionization, itself. Star formation suppression, for example, may explain the observed underabundance of Local Group dwarfs relative to N-body predictions for Cold Dark Matter. Reionization modelling requires simulating volumes large enough [~(100Mpc)^3] to sample reionization "patchiness", while reso… ▽ More Cosmic reionization by starlight from early galaxies affected their evolution, thereby impacting reionization, itself. Star formation suppression, for example, may explain the observed underabundance of Local Group dwarfs relative to N-body predictions for Cold Dark Matter. Reionization modelling requires simulating volumes large enough [~(100Mpc)^3] to sample reionization "patchiness", while resolving millions of galaxy sources above ~10^8 Msun , combining gravitational and gas dynamics with radiative transfer. Modelling the Local Group requires initial cosmological density fluctuations pre-selected to form the well-known structures of the local universe today. Cosmic Dawn ("CoDa") is the first such fully-coupled, radiation-hydrodynamics simulation of reionization of the local universe. Our new hybrid CPU-GPU code, RAMSES-CUDATON, performs hundreds of radiative transfer and ionization rate-solver timesteps on the GPUs for each hydro-gravity timestep on the CPUs. CoDa simulated (91Mpc)^3 with 4096^3 particles and cells, to redshift 4.23, on ORNL supercomputer Titan, utilizing 8192 cores and 8192 GPUs. Global reionization ended slightly later than observed. However, a simple temporal rescaling which brings the evolution of ionized fraction into agreement with observations also reconciles ionizing flux density, cosmic star formation history, CMB electron scattering optical depth and galaxy UV luminosity function with their observed values. Photoionization heating suppressed the star formation of haloes below ~2 x 10^9 Msun , decreasing the abun- dance of faint galaxies around MAB_1600 = [-10,-12]. For most of reionization, star formation was dominated by haloes between 10^10 - 10^11 Msun , so low-mass halo suppression was not reflected by a distinct feature in the global star formation history. (Abridged) △ Less

Submitted 6 September, 2016; v1 submitted 30 October, 2015; originally announced November 2015.

Comments: 26 pages, 16 figures, accepted in MNRAS

How scalar-field dark matter may conspire to facilitate baryogenesis at the electroweak scale

Authors: Tanja Rindler-Daller, Bohua Li, Paul R. Shapiro, Marek Lewicki, James D. Wells

Abstract: The cosmic evolution of a dark matter model which behaves relativistically in the early Universe is explored. Dark matter is described as a complex scalar field, whose earliest evolution is characterized by a stiff equation of state ($p \simeq ρ$). In this phase, it is the dominant component in the Universe. We present constraints from Big Bang nucleosynthesis and primordial gravity waves from inf… ▽ More The cosmic evolution of a dark matter model which behaves relativistically in the early Universe is explored. Dark matter is described as a complex scalar field, whose earliest evolution is characterized by a stiff equation of state ($p \simeq ρ$). In this phase, it is the dominant component in the Universe. We present constraints from Big Bang nucleosynthesis and primordial gravity waves from inflation. Also, we study how the associated enhanced expansion rate due to the stiff phase might facilitate a first-order electroweak symmetry breaking phase transition, in light of the recently measured value of the Higgs boson mass. △ Less

Submitted 30 October, 2015; v1 submitted 28 October, 2015; originally announced October 2015.

Comments: Presentation at the DPF2015 Meeting of the American Physical Society Division of Particles and Fields, Ann Arbor, Michigan, August 4-8, 2015; 8 pages; 3 figures

Report number: MCTP-15-25; TCC-012-15

The Impact of Nonlinear Structure Formation on the Power Spectrum of Transverse Momentum Fluctuations and the Kinetic Sunyaev-Zel'dovich Effect

Authors: Hyunbae Park, Eiichiro Komatsu, Paul R. Shapiro, Jun Koda, Yi Mao

Abstract: Cosmological transverse momentum fields, whose directions are perpendicular to Fourier wave vectors, induce temperature anisotropies in the cosmic microwave background via the kinetic Sunyaev-Zeldovich (kSZ) effect. The transverse momentum power spectrum contains the four-point function of density and velocity fields, $\langleδδv v\rangle$. In the post-reionization epoch, nonlinear effects dominat… ▽ More Cosmological transverse momentum fields, whose directions are perpendicular to Fourier wave vectors, induce temperature anisotropies in the cosmic microwave background via the kinetic Sunyaev-Zeldovich (kSZ) effect. The transverse momentum power spectrum contains the four-point function of density and velocity fields, $\langleδδv v\rangle$. In the post-reionization epoch, nonlinear effects dominate in the power spectrum. We use perturbation theory and cosmological $N$-body simulations to calculate this nonlinearity. We derive the next-to-leading order expression for the power spectrum with a particular emphasis on the connected term that has been ignored in the literature. While the contribution from the connected term on small scales ($k>0.1\,h\,\rm{Mpc}^{-1}$) is subdominant relative to the unconnected term, we find that its contribution to the kSZ power spectrum at $\ell = 3000$ at $z<6$ can be as large as ten percent of the unconnected term, which would reduce the allowed contribution from the reionization epoch ($z>6$) by twenty percent. The power spectrum of transverse momentum on large scales is expected to scale as $k^2$ as a consequence of momentum conservation. We show that both the leading and the next-to-leading order terms satisfy this scaling. In particular, we find that both of the unconnected and connected terms are necessary to reproduce $k^2$. △ Less

Submitted 8 November, 2015; v1 submitted 16 June, 2015; originally announced June 2015.

Comments: 13 pages, 5 figures, Accepted to ApJ

Cosmology on the Largest Scales with the SKA

Authors: S. Camera, A. Raccanelli, P. Bull, D. Bertacca, X. Chen, P. G. Ferreira, M. Kunz, R. Maartens, Y. Mao, M. G. Santos, P. R. Shapiro, M. Viel, Y. Xu

Abstract: The study of the Universe on ultra-large scales is one of the major science cases for the Square Kilometre Array (SKA). The SKA will be able to probe a vast volume of the cosmos, thus representing a unique instrument, amongst next-generation cosmological experiments, for scrutinising the Universe's properties on the largest cosmic scales. Probing cosmic structures on extremely large scales will ha… ▽ More The study of the Universe on ultra-large scales is one of the major science cases for the Square Kilometre Array (SKA). The SKA will be able to probe a vast volume of the cosmos, thus representing a unique instrument, amongst next-generation cosmological experiments, for scrutinising the Universe's properties on the largest cosmic scales. Probing cosmic structures on extremely large scales will have many advantages. For instance, the growth of perturbations is well understood for those modes, since it falls fully within the linear regime. Also, such scales are unaffected by the poorly understood feedback of baryonic physics. On ultra-large cosmic scales, two key effects become significant: primordial non-Gaussianity and relativistic corrections to cosmological observables. Moreover, if late-time acceleration is driven not by dark energy but by modifications to general relativity, then such modifications should become apparent near and above the horizon scale. As a result, the SKA is forecast to deliver transformational constraints on non-Gaussianity and to probe gravity on super-horizon scales for the first time. △ Less

Submitted 19 January, 2015; v1 submitted 15 January, 2015; originally announced January 2015.

Comments: 18 pages, 3 figures; updated acknowledgments and references. This article is part of the 'SKA Cosmology Chapter, Advancing Astrophysics with the SKA (AASKA14) Conference, Giardini Naxos (Italy), June 9th-13th 2014'

The Linear Perturbation Theory of Reionization in Position-Space: Cosmological Radiative Transfer Along the Light-Cone

Authors: Yi Mao, Anson D'Aloisio, Benjamin D. Wandelt, Jun Zhang, Paul R. Shapiro

Abstract: The linear perturbation theory of inhomogeneous reionization (LPTR) has been developed as an analytical tool for predicting the global ionized fraction and large-scale power spectrum of ionized density fluctuations during reionization. In the original formulation of the LPTR, the ionization balance and radiative transfer equations are linearized and solved in Fourier space. However, the LPTR's app… ▽ More The linear perturbation theory of inhomogeneous reionization (LPTR) has been developed as an analytical tool for predicting the global ionized fraction and large-scale power spectrum of ionized density fluctuations during reionization. In the original formulation of the LPTR, the ionization balance and radiative transfer equations are linearized and solved in Fourier space. However, the LPTR's approximation to the full solution of the radiative transfer equation is not straightforward to interpret, since the latter is most intuitively conceptualized in position space. To bridge the gap between the LPTR and the language of numerical radiative transfer, we present a new, equivalent, position-space formulation of the LPTR that clarifies the approximations it makes and facilitates its interpretation. We offer a comparison between the LPTR and the excursion-set model of reionization (ESMR), and demonstrate the built-in capability of the LPTR to explore a wide range of reionization scenarios, and to go beyond the ESMR in exploring scenarios involving X-rays. △ Less

Submitted 14 April, 2015; v1 submitted 25 November, 2014; originally announced November 2014.

Comments: 10 pages, 1 figure, 1 table. Accepted by PRD with minor changes

Journal ref: Phys. Rev. D 91, 083015 (2015)

Nonlinear Bias of Cosmological Halo Formation in the Early Universe

Authors: Kyungjin Ahn, Ilian T. Iliev, Paul R. Shapiro, ChaiChalit Srisawat

Abstract: We present estimates of the nonlinear bias of cosmological halo formation, spanning a wide range in the halo mass from $\sim 10^{5} M_\odot$ to $\sim 10^{12} M_\odot$, based upon both a suite of high-resolution cosmological N-body simulations and theoretical predictions. The halo bias is expressed in terms of the mean bias and stochasticity as a function of local overdensity ($δ$), under different… ▽ More We present estimates of the nonlinear bias of cosmological halo formation, spanning a wide range in the halo mass from $\sim 10^{5} M_\odot$ to $\sim 10^{12} M_\odot$, based upon both a suite of high-resolution cosmological N-body simulations and theoretical predictions. The halo bias is expressed in terms of the mean bias and stochasticity as a function of local overdensity ($δ$), under different filtering scales, which is realized as the density of individual cells in uniform grids. The sampled overdensities span a range wide enough to provide the fully nonlinear bias effect on the formation of haloes. A strong correlation between $δ$ and halo population overdensity $δ_h$ is found, along with sizable stochasticity. We find that the empirical mean halo bias matches, with good accuracy, the prediction by the peak-background split method based on the excursion set formalism, as long as the empirical, globally-averaged halo mass function is used. Consequently, this bias formalism is insensitive to uncertainties caused by varying halo identification schemes, and can be applied generically. We also find that the probability distribution function of biased halo numbers has wider distribution than the pure Poisson shot noise, which is attributed to the sub-cell scale halo correlation. We explicitly calculate this correlation function and show that both overdense and underdense regions have positive correlation, leading to stochasticity larger than the Poisson shot noise in the range of haloes and halo-collapse epochs we study. △ Less

Submitted 1 April, 2015; v1 submitted 9 July, 2014; originally announced July 2014.

Comments: 18 pages, 8 figures, in press for publication in MNRAS; supplementary material (additional 16 figures) separately supplied (supplement.pdf) as a part of source files

Light cone effect on the reionization 21-cm signal II: Evolution, anisotropies and observational implications

Authors: Kanan K. Datta, Hannes Jensen, Suman Majumdar, Garrelt Mellema, Ilian T. Iliev, Yi Mao, Paul R. Shapiro, Kyungjin Ahn

Abstract: Measurements of the HI 21-cm power spectra from the reionization epoch will be influenced by the evolution of the signal along the line-of-sight direction of any observed volume. We use numerical as well as semi-numerical simulations of reionization in a cubic volume of 607 Mpc across to study this so-called light cone effect on the HI 21-cm power spectrum. We find that the light cone effect has t… ▽ More Measurements of the HI 21-cm power spectra from the reionization epoch will be influenced by the evolution of the signal along the line-of-sight direction of any observed volume. We use numerical as well as semi-numerical simulations of reionization in a cubic volume of 607 Mpc across to study this so-called light cone effect on the HI 21-cm power spectrum. We find that the light cone effect has the largest impact at two different stages of reionization: one when reionization is $\sim 20\%$ and other when it is $\sim 80\%$ completed. We find a factor of $\sim 4$ amplification of the power spectrum at the largest scale available in our simulations. We do not find any significant anisotropy in the 21-cm power spectrum due to the light cone effect. We argue that for the power spectrum to become anisotropic, the light cone effect would have to make the ionized bubbles significantly elongated or compressed along the line-of-sight, which would require extreme reionization scenarios. We also calculate the two-point correlation functions parallel and perpendicular to the line-of-sight and find them to differ. Finally, we calculate an optimum frequency bandwidth below which the light cone effect can be neglected when extracting power spectra from observations. We find that if one is willing to accept a $10 \%$ error due to the light cone effect, the optimum frequency bandwidth for $k= 0.056 \, \rm{Mpc}^{-1}$ is $\sim 7.5$ MHz. For $k = 0.15$ and $0.41 \, \rm{Mpc}^{-1}$ the optimum bandwidth is $\sim 11$ and $\sim 16$ MHz respectively. △ Less

Submitted 7 May, 2014; v1 submitted 3 February, 2014; originally announced February 2014.

Comments: 17 pages, accepted for publication in MNRAS, minor changes

Journal ref: MNRAS, 2014, Volume 442, Issue 2, p.1491-1506

Complex scalar field dark matter on galactic scales

Authors: Tanja Rindler-Daller, Paul R. Shapiro

Abstract: The nature of the cosmological dark matter remains elusive. Recent studies have advocated the possibility that dark matter could be composed of ultra-light, self-interacting bosons, forming a Bose-Einstein condensate in the very early Universe. We consider models which are charged under a global U(1)-symmetry such that the dark matter number is conserved. It can then be described as a classical co… ▽ More The nature of the cosmological dark matter remains elusive. Recent studies have advocated the possibility that dark matter could be composed of ultra-light, self-interacting bosons, forming a Bose-Einstein condensate in the very early Universe. We consider models which are charged under a global U(1)-symmetry such that the dark matter number is conserved. It can then be described as a classical complex scalar field which evolves in an expanding Universe. We present a brief review on the bounds on the model parameters from cosmological and galactic observations, along with the properties of galactic halos which result from such a dark matter candidate. △ Less

Submitted 17 April, 2014; v1 submitted 5 December, 2013; originally announced December 2013.

Comments: 17 pages; Minor corrections to reflect the published version

Journal ref: Mod. Phys. Lett. A, Vol. 29, No. 2 (2014) 1430002

The Effect of Large-Scale Structure on the Magnification of High-Redshift Sources by Cluster-Lenses

Authors: Anson D'Aloisio, Priyamvada Natarajan, Paul R. Shapiro

Abstract: Cluster gravitational lensing surveys like the Hubble Space Telescope Frontier Fields survey will detect distant galaxies 10-50 times fainter than any yet discovered. Using these surveys to measure the luminosity function of such faint, distant galaxies, however, requires that magnification maps built from the constraints of strongly-lensed images be accurate. For models that assume the cluster an… ▽ More Cluster gravitational lensing surveys like the Hubble Space Telescope Frontier Fields survey will detect distant galaxies 10-50 times fainter than any yet discovered. Using these surveys to measure the luminosity function of such faint, distant galaxies, however, requires that magnification maps built from the constraints of strongly-lensed images be accurate. For models that assume the cluster and nearby (correlated) structures are the only significant sources of lensing, a potential source of error in these maps comes from the fact that light rays also suffer weak deflections by uncorrelated large-scale structure along the line-of-sight, i.e. cosmic weak lensing (CWL). To demonstrate the magnitude of this effect, we calculate the magnification change which results when the same cluster-lens is placed along different lines of sight. Using a simple density profile for a cluster-lens at z~0.3-0.5 and the power spectrum of the matter density fluctuations responsible for CWL, we show that the typical magnifications of ~5(10) of sources at z=6-10 can differ by ~10-20(20-30)% from one line-of-sight to another. However, these fluctuations rise to greater than order unity near critical curves, indicating that CWL tends to make its greatest contribution to the most magnified images. We conclude that the neglect of CWL in determining the intrinsic luminosities of highly-magnified galaxies may introduce errors significant enough to warrant further effort to include this contribution in cluster-lens modeling. We suggest that methods of modeling CWL in galaxy-strong-lensing systems should be generalized to cluster-lensing systems. △ Less

Submitted 13 October, 2014; v1 submitted 7 November, 2013; originally announced November 2013.

Comments: Accepted by MNRAS. Changes to text for improved clarity, minor typos corrected, and references updated. All results remain the same

Simulating cosmic reionization: How large a volume is large enough?

Authors: Ilian T. Iliev, Garrelt Mellema, Kyungjin Ahn, Paul R. Shapiro, Yi Mao, Ue-Li Pen

Abstract: We present the largest-volume (425 Mpc/h=607 Mpc on a side) full radiative transfer simulation of cosmic reionization to date. We show that there is significant additional power in density fluctuations at very large scales. We systematically investigate the effects this additional power has on the progress, duration and features of reionization, as well as on selected reionization observables. We… ▽ More We present the largest-volume (425 Mpc/h=607 Mpc on a side) full radiative transfer simulation of cosmic reionization to date. We show that there is significant additional power in density fluctuations at very large scales. We systematically investigate the effects this additional power has on the progress, duration and features of reionization, as well as on selected reionization observables. We find that comoving simulation volume of ~100 Mpc/h per side is sufficient for deriving a convergent mean reionization history, but that the reionization patchiness is significantly underestimated. We use jackknife splitting to quantify the convergence of reionization properties with simulation volume for both mean-density and variable-density sub-regions. We find that sub-volumes of ~100 Mpc/h per side or larger yield convergent reionization histories, except for the earliest times, but smaller volumes of ~50 Mpc/h or less are not well converged at any redshift. Reionization history milestones show significant scatter between the sub-volumes, of Delta z=0.6-1 for ~50 Mpc/h volumes, decreasing to Delta z=0.3-0.5 for ~100 Mpc/h volumes, and $Δz$~0.1 for ~200 Mpc/h volumes. If we only consider mean-density sub-regions the scatter decreases, but remains at Delta z~0.1-0.2 for the different size sub-volumes. Consequently, many potential reionization observables like 21-cm rms, 21-cm PDF skewness and kurtosis all show good convergence for volumes of ~200 Mpc/h, but retain considerable scatter for smaller volumes. In contrast, the three-dimensional 21-cm power spectra at large scales (k<0.25 h/Mpc) do not fully converge for any sub-volume size. These additional large-scale fluctuations significantly enhance the 21-cm fluctuations, which should improve the prospects of detection considerably, given the lower foregrounds and greater interferometer sensitivity at higher frequencies. (abridged) △ Less

Submitted 28 October, 2013; originally announced October 2013.

Comments: 20 pages, 19 figures (most in color)

Cosmological Constraints on Bose-Einstein-Condensed Scalar Field Dark Matter

Authors: Bohua Li, Tanja Rindler-Daller, Paul R. Shapiro

Abstract: (Abridged) Despite the success of Cold Dark Matter (CDM) in explaining a wide range of observations, the microscopic nature of dark matter is still unknown. Attempts to detect WIMPs, the most commonly studied form, have not yet succeeded. Some predictions are even in apparent conflict with observations (e.g. cuspy-cored halos and `missing satellite problem'). Therefore, it is important to consider… ▽ More (Abridged) Despite the success of Cold Dark Matter (CDM) in explaining a wide range of observations, the microscopic nature of dark matter is still unknown. Attempts to detect WIMPs, the most commonly studied form, have not yet succeeded. Some predictions are even in apparent conflict with observations (e.g. cuspy-cored halos and `missing satellite problem'). Therefore, it is important to consider different forms of CDM. We focus on dark matter comprised of ultralight bosons that form a Bose-Einstein Condensate (BEC), described by a complex scalar field, and solve the Klein-Gordon and Einstein field equations to evolve the Friedmann-Robertson-Walker universe. We find, in addition to phases of radiation-, matter- and Λ-domination, an earlier phase of scalar-field-domination. Moreover, while WIMP CDM is non-relativistic at all times after it decouples, BEC scalar field dark matter (SFDM) is relativistic at early times, evolving from stiff to radiation-like, before it becomes non-relativistic (CDM-like) at late times. The timing of transitions between these phases yields fundamental constraints on SFDM model parameters, particle mass m and self-interaction coupling strength λ. We derive the range of particle parameters required to match observations of the evolving background universe, including the CMB and abundances of light elements produced by BBN, characterized by N_eff, the effective number of neutrino species, and the epoch of matter-radiation equality z_eq. This yields m >= 2.4*10^{-21} eV/c^2 and 9.5*10^{-19} eV^{-1}cm^3 <= λ/(mc^2)^2 <= 4*10^{-17} eV^{-1}cm^3. Our model accommodates current observations in which N_eff at BBN is higher than at z_eq, as probed by the CMB, otherwise unexplained by WIMP CDM. SFDM without self-interaction (`Fuzzy Dark Matter') cannot satisfy current BBN constraints within 68% confidence and is therefore disfavored. △ Less

Submitted 29 October, 2014; v1 submitted 22 October, 2013; originally announced October 2013.

Comments: 21 pages, 6 figures, minor updates to accord with the published version in Phys. Rev. D

Journal ref: Phys. Rev. D 89, 083536 (2014)

Primordial Non-Gaussianity Estimation using 21 cm Tomography from the Epoch of Reionization

Authors: Yi Mao, Anson D'Aloisio, Jun Zhang, Paul R. Shapiro

Abstract: Measuring the small primordial nonGaussianity (PNG) predicted by cosmic inflation theories may help diagnose them. The detectability of PNG by its imprint on the 21cm power spectrum from the epoch of reionization is reassessed here in terms of $f_{NL}$, the local nonlinearity parameter. We find that an optimum, multi-frequency observation by SKA can achieve $Δf_{NL} \sim 3$ (comparable to recent P… ▽ More Measuring the small primordial nonGaussianity (PNG) predicted by cosmic inflation theories may help diagnose them. The detectability of PNG by its imprint on the 21cm power spectrum from the epoch of reionization is reassessed here in terms of $f_{NL}$, the local nonlinearity parameter. We find that an optimum, multi-frequency observation by SKA can achieve $Δf_{NL} \sim 3$ (comparable to recent Planck CMB limits), while a cosmic-variance-limited array of this size like Omniscope can even detect $Δf_{NL} \sim 0.2$. This substantially revises the methods and results of previous work. △ Less

Submitted 24 October, 2013; v1 submitted 1 May, 2013; originally announced May 2013.

Comments: Accepted by PRD with minor changes. References added and updated

Journal ref: Phys. Rev. D 88, 081303 (R) (2013)

The scale-dependent signature of primordial non-Gaussianity in the large-scale structure of cosmic reionization

Authors: Anson D'Aloisio, Jun Zhang, Paul R. Shapiro, Yi Mao

Abstract: (ABRIDGED)The rise of cosmic structure depends upon the statistical distribution of initial density fluctuations generated by inflation. While the simplest models predict an almost perfectly Gaussian distribution, more-general models predict a level of primordial non-Gaussianity (PNG) that observations might yet be sensitive enough to detect. Recent Planck Collaboration measurements of the CMB tem… ▽ More (ABRIDGED)The rise of cosmic structure depends upon the statistical distribution of initial density fluctuations generated by inflation. While the simplest models predict an almost perfectly Gaussian distribution, more-general models predict a level of primordial non-Gaussianity (PNG) that observations might yet be sensitive enough to detect. Recent Planck Collaboration measurements of the CMB temperature anisotropy bispectrum significantly tighten the observational limits, but they are still far from the PNG level predicted by the simplest models of inflation. Probing levels below CMB sensitivities will require other methods, such as searching for the statistical imprint of PNG on galactic halo clustering. During the epoch of reionization (EoR), the first stars and galaxies released radiation into the intergalactic medium (IGM) that created ionized patches whose large-scale geometry and evolution reflected the underlying abundance and large-scale clustering of the star-forming galaxies. This statistical connection between ionized patches in the IGM and galactic halos suggests that observing reionization may be another way to constrain PNG. We employ the linear perturbation theory of reionization and semi-analytic models based on the excursion-set formalism to model the effects of PNG on the EoR. We quantify the effects of PNG on the large-scale structure of reionization by deriving the ionized density bias, i.e. ratio of ionized atomic to total matter overdensities in Fourier space, at small wavenumber. Just as previous studies found that PNG creates a scale-dependent signature in the halo bias, so, too, we find a scale-dependent signature in the ionized density bias. Our results, which differ significantly from previous attempts in the literature to characterize this PNG signature, will be applied elsewhere to predict its observable consequences, e.g. in the cosmic 21cm background. △ Less

Submitted 30 May, 2013; v1 submitted 23 April, 2013; originally announced April 2013.

Comments: Accepted by MNRAS with minor changes. Minor typos corrected and references added

Probing reionization with LOFAR using 21-cm redshift space distortions

Authors: Hannes Jensen, Kanan K. Datta, Garrelt Mellema, Emma Chapman, Filipe B. Abdalla, Ilian T. Iliev, Yi Mao, Mario G. Santos, Paul R. Shapiro, Saleem Zaroubi, G. Bernardi, M. A. Brentjens, A. G. de Bruyn, B. Ciardi, G. J. A. Harker, V. Jelić, S. Kazemi, L. V. E. Koopmans, P. Labropoulos, O. Martinez, A. R. Offringa, V. N. Pandey, J. Schaye, R. M. Thomas, V. Veligatla , et al. (2 additional authors not shown)

Abstract: One of the most promising ways to study the epoch of reionization (EoR) is through radio observations of the redshifted 21-cm line emission from neutral hydrogen. These observations are complicated by the fact that the mapping of redshifts to line-of-sight positions is distorted by the peculiar velocities of the gas. Such distortions can be a source of error if they are not properly understood, bu… ▽ More One of the most promising ways to study the epoch of reionization (EoR) is through radio observations of the redshifted 21-cm line emission from neutral hydrogen. These observations are complicated by the fact that the mapping of redshifts to line-of-sight positions is distorted by the peculiar velocities of the gas. Such distortions can be a source of error if they are not properly understood, but they also encode information about cosmology and astrophysics. We study the effects of redshift space distortions on the power spectrum of 21-cm radiation from the EoR using large scale $N$-body and radiative transfer simulations. We quantify the anisotropy introduced in the 21-cm power spectrum by redshift space distortions and show how it evolves as reionization progresses and how it relates to the underlying physics. We go on to study the effects of redshift space distortions on LOFAR observations, taking instrument noise and foreground subtraction into account. We find that LOFAR should be able to directly observe the power spectrum anisotropy due to redshift space distortions at spatial scales around $k \sim 0.1$ Mpc$^{-1}$ after $\gtrsim$ 1000 hours of integration time. At larger scales, sample errors become a limiting factor, while at smaller scales detector noise and foregrounds make the extraction of the signal problematic. Finally, we show how the astrophysical information contained in the evolution of the anisotropy of the 21-cm power spectrum can be extracted from LOFAR observations, and how it can be used to distinguish between different reionization scenarios. △ Less

Submitted 19 July, 2013; v1 submitted 22 March, 2013; originally announced March 2013.

Comments: Accepted for publication in MNRAS

The Kinetic Sunyaev-Zel'dovich effect as a probe of the physics of cosmic reionization: the effect of self-regulated reionization

Authors: Hyunbae Park, Paul R. Shapiro, Eiichiro Komatsu, Ilian T. Iliev, Kyungjin Ahn, Garrelt Mellema

Abstract: We calculate the angular power spectrum of the Cosmic Microwave Background (CMB) temperature fluctuations induced by the kinetic Sunyaev-Zel'dovich (kSZ) effect from the epoch of reionization (EOR). We use detailed N-body+radiative transfer simulations to follow inhomogeneous reionization of the intergalactic medium (IGM). For the first time we take into account the "self-regulation" of reionizati… ▽ More We calculate the angular power spectrum of the Cosmic Microwave Background (CMB) temperature fluctuations induced by the kinetic Sunyaev-Zel'dovich (kSZ) effect from the epoch of reionization (EOR). We use detailed N-body+radiative transfer simulations to follow inhomogeneous reionization of the intergalactic medium (IGM). For the first time we take into account the "self-regulation" of reionization: star formation in low-mass dwarf galaxies (10^8 M_\sun \lesssim M \lesssim 10^9 M_\sun) or minihalos (10^5 M_\sun \lesssim M \lesssim 10^8 M_\sun) is suppressed if these halos form in the regions that were already ionized or Lyman-Werner dissociated. Some previous work suggested that the amplitude of the kSZ power spectrum from the EOR can be described by a two-parameter family: the epoch of half ionization and the duration of reionization. However, we argue that this picture applies only to simple forms of the reionization history which are roughly symmetric about the half-ionization epoch. In self-regulated reionization, the universe begins to be ionized early, maintains a low level of ionization for an extended period, and then finishes reionization as soon as high-mass atomically-cooling halos dominate. While inclusion of self-regulation affects the amplitude of the kSZ power spectrum only modestly (\sim 10 %), it can change the duration of reionization by a factor of more than two. We conclude that the simple two-parameter family does not capture the effect of a physical, yet complex, reionization history caused by self-regulation. When added to the post-reionization kSZ contribution, our prediction for the total kSZ power spectrum is below the current upper bound from the South Pole Telescope. Therefore, the current upper bound on the kSZ effect from the EOR is consistent with our understanding of the physics of reionization. △ Less

Submitted 16 January, 2013; originally announced January 2013.

Comments: Submitted to ApJ; 16 pages, 8 figures

The halo mass function through the cosmic ages

Authors: William A. Watson, Ilian T. Iliev, Anson D'Aloisio, Alexander Knebe, Paul R. Shapiro, Gustavo Yepes

Abstract: In this paper we investigate how the halo mass function evolves with redshift, based on a suite of very large (with N_p = 3072^3 - 6000^3 particles) cosmological N-body simulations. Our halo catalogue data spans a redshift range of z = 0-30, allowing us to probe the mass function from the dark ages to the present. We utilise both the Friends-of-Friends (FOF) and Spherical Overdensity (SO) halofind… ▽ More In this paper we investigate how the halo mass function evolves with redshift, based on a suite of very large (with N_p = 3072^3 - 6000^3 particles) cosmological N-body simulations. Our halo catalogue data spans a redshift range of z = 0-30, allowing us to probe the mass function from the dark ages to the present. We utilise both the Friends-of-Friends (FOF) and Spherical Overdensity (SO) halofinding methods to directly compare the mass function derived using these commonly used halo definitions. The mass function from SO haloes exhibits a clear evolution with redshift, especially during the recent era of dark energy dominance (z < 1). We provide a redshift-parameterised fit for the SO mass function valid for the entire redshift range to within ~20% as well as a scheme to calculate the mass function for haloes with arbitrary overdensities. The FOF mass function displays a weaker evolution with redshift. We provide a `universal' fit for the FOF mass function, fitted to data across the entire redshift range simultaneously, and observe redshift evolution in our data versus this fit. The relative evolution of the mass functions derived via the two methods is compared and we find that the mass functions most closely match at z=0. The disparity at z=0 between the FOF and SO mass functions resides in their high mass tails where the collapsed fraction of mass in SO haloes is ~80% of that in FOF haloes. This difference grows with redshift so that, by z>20, the SO algorithm finds a ~50-80% lower collapsed fraction in high mass haloes than does the FOF algorithm, due in part to the significant over-linking effects known to affect the FOF method. △ Less

Submitted 9 November, 2013; v1 submitted 1 December, 2012; originally announced December 2012.

Comments: v4, 16 pages, 16 colour figures. Changed to match MNRAS print version. NOTE: v1 of this paper has a typo in the fitting function. Please ensure you use the latest version

Journal ref: 2013MNRAS.433.1230W

Will Nonlinear Peculiar Velocity and Inhomogeneous Reionization Spoil 21cm Cosmology from the Epoch of Reionization?

Authors: Paul R. Shapiro, Yi Mao, Ilian T. Iliev, Garrelt Mellema, Kanan K. Datta, Kyungjin Ahn, Jun Koda

Abstract: The 21cm background from the epoch of reionization is a promising cosmological probe: line-of-sight velocity fluctuations distort redshift, so brightness fluctuations in Fourier space depend upon angle, which linear theory shows can separate cosmological from astrophysical information. Nonlinear fluctuations in ionization, density and velocity change this, however. The validity and accuracy of the… ▽ More The 21cm background from the epoch of reionization is a promising cosmological probe: line-of-sight velocity fluctuations distort redshift, so brightness fluctuations in Fourier space depend upon angle, which linear theory shows can separate cosmological from astrophysical information. Nonlinear fluctuations in ionization, density and velocity change this, however. The validity and accuracy of the separation scheme are tested here for the first time, by detailed reionization simulations. The scheme works reasonably well early in reionization (< 40% ionized), but not late (> 80% ionized). △ Less

Submitted 28 February, 2013; v1 submitted 8 November, 2012; originally announced November 2012.

Comments: 2 figures, matches published PRL version

Simulating Cosmic Reionization and the Radiation Backgrounds from the Epoch of Reionization

Authors: Paul R. Shapiro, Ilian T. Iliev, Garrelt Mellema, Kyungjin Ahn, Yi Mao, Martina Friedrich, Kanan Datta, Hyunbae Park, Eiichiro Komatsu, Elizabeth Fernandez, Jun Koda, Mia Bovill, Ue-Li Pen

Abstract: Large-scale reionization simulations are described which combine the results of cosmological N-body simulations that model the evolving density and velocity fields and identify the galactic halo sources, with ray-tracing radiative transfer calculations which model the nonequilibrium ionization of the intergalactic medium. These simulations have been used to predict some of the signature effects of… ▽ More Large-scale reionization simulations are described which combine the results of cosmological N-body simulations that model the evolving density and velocity fields and identify the galactic halo sources, with ray-tracing radiative transfer calculations which model the nonequilibrium ionization of the intergalactic medium. These simulations have been used to predict some of the signature effects of reionization on cosmic radiation backgrounds, including the CMB, near-IR, and redshifted 21cm backgrounds. We summarize some of our recent progress in this work, and address the question of whether observations of such signature effects can be used to distinguish the relative contributions of galaxies of different masses to reionization. △ Less

Submitted 2 November, 2012; originally announced November 2012.

Comments: Based on invited review talk by Paul Shapiro at First Stars IV - From Hayashi to the Future, Kyoto, Japan, May 25, 2012; talk slides available at http://tpweb2.phys.konan-u.ac.jp/~FirstStar4/presentation_files/PShapiro.pdf

Journal ref: Shapiro, P.R. et al., First Stars IV - From Hayashi to the Future -, AIP Conf. Proc. 1480, 248-260 (2012); View online: http://dx.doi.org/10.1063/1.4754363

Finding new signature effects on galactic dynamics to constrain Bose-Einstein-condensed cold dark matter

Authors: Tanja Rindler-Daller, Paul R. Shapiro

Abstract: If cosmological cold dark matter (CDM) consists of light enough bosonic particles that their phase-space density exceeds unity, they will comprise a Bose-Einstein condensate (BEC). The nature of this BEC-CDM as a quantum fluid may then distinguish it dynamically from the standard form of CDM involving a collisionless gas of non-relativistic particles that interact purely gravitationally. We summar… ▽ More If cosmological cold dark matter (CDM) consists of light enough bosonic particles that their phase-space density exceeds unity, they will comprise a Bose-Einstein condensate (BEC). The nature of this BEC-CDM as a quantum fluid may then distinguish it dynamically from the standard form of CDM involving a collisionless gas of non-relativistic particles that interact purely gravitationally. We summarize some of the dynamical properties of BEC-CDM that may lead to observable signatures in galactic halos and present some of the bounds on particle mass and self-interaction coupling strength that result from a comparison with observed galaxies. △ Less

Submitted 17 April, 2014; v1 submitted 9 September, 2012; originally announced September 2012.

Comments: 10 pages, 5 figures; Minor corrections to reflect the published version

Journal ref: Astrophys. Space Sci. Proc. Vol.38 (2014) p.163-182

Detecting the Rise and Fall of the First Stars by Their Impact on Cosmic Reionization

Authors: Kyungjin Ahn, Ilian T. Iliev, Paul R. Shapiro, Garrelt Mellema, Jun Koda, Yi Mao

Abstract: The intergalactic medium was reionized before redshift z~6, most likely by starlight which escaped from early galaxies. The very first stars formed when hydrogen molecules (H2) cooled gas inside the smallest galaxies, minihalos of mass between 10^5 and 10^8 solar masses. Although the very first stars began forming inside these minihalos before redshift z~40, their contribution has, to date, been i… ▽ More The intergalactic medium was reionized before redshift z~6, most likely by starlight which escaped from early galaxies. The very first stars formed when hydrogen molecules (H2) cooled gas inside the smallest galaxies, minihalos of mass between 10^5 and 10^8 solar masses. Although the very first stars began forming inside these minihalos before redshift z~40, their contribution has, to date, been ignored in large-scale simulations of this cosmic reionization. Here we report results from the first reionization simulations to include these first stars and the radiative feedback that limited their formation, in a volume large enough to follow the crucial spatial variations that influenced the process and its observability. We show that, while minihalo stars stopped far short of fully ionizing the universe, reionization began much earlier with minihalo sources than without, and was greatly extended, which boosts the intergalactic electron-scattering optical depth and the large-angle polarization fluctuations of the cosmic microwave background significantly. Although within current WMAP uncertainties, this boost should be readily detectable by Planck. If reionization ended as late as z_ov<~7, as suggested by other observations, Planck will thereby see the signature of the first stars at high redshift, currently undetectable by other probes. △ Less

Submitted 2 August, 2012; v1 submitted 21 June, 2012; originally announced June 2012.

Comments: Accepted to ApJL; Some minor changes made and several new references added

On the Use of Ly-alpha Emitters as Probes of Reionization

Authors: Hannes Jensen, Peter Laursen, Garrelt Mellema, Ilian T. Iliev, Jesper Sommer-Larsen, Paul R. Shapiro

Abstract: We use numerical simulations to study the effects of the patchiness of a partly reionized intergalactic medium (IGM) on the observability of Ly-alpha emitters (LAEs) at high redshifts (z ~ 6). We present a new model that divides the Ly-alpha radiative transfer into a (circum-)galactic and an extragalactic (IGM) part, and investigate how the choice of intrinsic line model affects the IGM transmissi… ▽ More We use numerical simulations to study the effects of the patchiness of a partly reionized intergalactic medium (IGM) on the observability of Ly-alpha emitters (LAEs) at high redshifts (z ~ 6). We present a new model that divides the Ly-alpha radiative transfer into a (circum-)galactic and an extragalactic (IGM) part, and investigate how the choice of intrinsic line model affects the IGM transmission results. We use our model to study the impact of neutral hydrogen on statistical observables such as the Ly-alpha restframe equivalent width (REW) distribution, the LAE luminosity function and the two-point correlation function. We find that if the observed changes in LAE luminosity functions and equivalent width distributions between z ~ 6 and z ~ 7 are to be explained by an increased IGM neutral fraction alone, we require an extremely late and rapid reionization scenario, where the Universe was ~ 40 % ionized at z = 7, ~ 50 % ionized at z = 6.5 and ~ 100 % ionized at z = 6. This is in conflict with other observations, suggesting that intrinsic LAE evolution at z > 6 cannot be completely neglected. We show how the two-point correlation function can provide more robust constraints once future observations obtain larger LAE samples, and provide predictions for the sample sizes needed to tell different reionization scenarios apart. △ Less

Submitted 18 October, 2012; v1 submitted 18 June, 2012; originally announced June 2012.

Comments: Accepted for publication in MNRAS

Halo statistics in non-Gaussian cosmologies: the collapsed fraction, conditional mass function, and halo bias from the path-integral excursion set method

Authors: Anson D'Aloisio, Jun Zhang, Donghui Jeong, Paul R. Shapiro

Abstract: Characterizing the level of primordial non-Gaussianity (PNG) in the initial conditions for structure formation is one of the most promising ways to test inflation and differentiate among different scenarios. The scale-dependent imprint of PNG on the large-scale clustering of galaxies and quasars has already been used to place significant constraints on the level of PNG in our observed Universe. Su… ▽ More Characterizing the level of primordial non-Gaussianity (PNG) in the initial conditions for structure formation is one of the most promising ways to test inflation and differentiate among different scenarios. The scale-dependent imprint of PNG on the large-scale clustering of galaxies and quasars has already been used to place significant constraints on the level of PNG in our observed Universe. Such measurements depend upon an accurate and robust theory of how PNG affects the bias of galactic halos relative to the underlying matter density field. We improve upon previous work by employing a more general analytical method - the path-integral extension of the excursion set formalism - which is able to account for the non-Markovianity caused by PNG in the random-walk model used to identify halos in the initial density field. This non-Markovianity encodes information about environmental effects on halo formation which have so far not been taken into account in analytical bias calculations. We compute both scale-dependent and -independent corrections to the halo bias, along the way presenting an expression for the conditional collapsed fraction for the first time, and a new expression for the conditional halo mass function. To leading order in our perturbative calculation, we recover the halo bias results of Desjacques et. al. (2011), including the new scale-dependent correction reported there. However, we show that the non-Markovian dynamics from PNG can lead to marked differences in halo bias when next-to-leading order terms are included. We quantify these differences here. [abridged] △ Less

Submitted 22 October, 2012; v1 submitted 14 June, 2012; originally announced June 2012.

Comments: Accepted for publication in MNRAS. Includes minor revisions recommended by referee, slightly revised notation for clarity, and corrected typos