Showing 1–50 of 214 results for author: Gnedin, N

On the minimum number of radiation field parameters to specify gas cooling and heating functions

Authors: David Robinson, Camille Avestruz, Nickolay Y. Gnedin

Abstract: Fast and accurate approximations of gas cooling and heating functions are needed for hydrodynamic galaxy simulations. We use machine learning to analyze atomic gas cooling and heating functions in the presence of a generalized incident local radiation field computed by Cloudy. We characterize the radiation field through binned radiation field intensities instead of the photoionization rates used i… ▽ More Fast and accurate approximations of gas cooling and heating functions are needed for hydrodynamic galaxy simulations. We use machine learning to analyze atomic gas cooling and heating functions in the presence of a generalized incident local radiation field computed by Cloudy. We characterize the radiation field through binned radiation field intensities instead of the photoionization rates used in our previous work. We find a set of 6 energy bins whose intensities exhibit relatively low correlation. We use these bins as features to train machine learning models to predict Cloudy cooling and heating functions at fixed metallicity. We compare the relative SHAP importance of the features. From the SHAP analysis, we identify a feature subset of 3 energy bins ($0.5-1, 1-4$, and $13-16 \, \mathrm{Ry}$) with the largest importance and train additional models on this subset. We compare the mean squared errors and distribution of errors on both the entire training data table and a randomly selected 20% test set withheld from model training. The machine learning models trained with 3 and 6 bins, as well as 3 and 4 photoionization rates, have comparable accuracy everywhere. We conclude that 3 energy bins (or 3 analogous photoionization rates: molecular hydrogen photodissociation, neutral hydrogen HI, and fully ionized carbon CVI) are sufficient to characterize the dependence of the gas cooling and heating functions on our assumed incident radiation field model. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: 12 + 5 pages, 7 + 4 figures. To be submitted to OJA. Comments welcome!

On the Physical Nature of Ly$α$ Transmission Spikes in High Redshift Quasar Spectra

Authors: Hanjue Zhu, Nickolay Gnedin, Camille Avestruz

Abstract: We investigate Lyman-alpha (Ly$α$) transmission spikes at $5.2 < z < 6.8$ using synthetic quasar spectra from the ``Cosmic Reionization On Computers" simulations. We focus on understanding the relationship between these spikes and the properties of the intergalactic medium (IGM). Disentangling the complex interplay between IGM physics and the influence of galaxies on the generation of these spikes… ▽ More We investigate Lyman-alpha (Ly$α$) transmission spikes at $5.2 < z < 6.8$ using synthetic quasar spectra from the ``Cosmic Reionization On Computers" simulations. We focus on understanding the relationship between these spikes and the properties of the intergalactic medium (IGM). Disentangling the complex interplay between IGM physics and the influence of galaxies on the generation of these spikes presents a significant challenge. To address this, we employ Explainable Boosting Machines, an interpretable machine learning algorithm, to quantify the relative impact of various IGM properties on the Ly$α$ flux. Our findings reveal that gas density is the primary factor influencing absorption strength, followed by the intensity of background radiation and the temperature of the IGM. Ionizing radiation from local sources (i.e. galaxies) appears to have a minimal effect on Ly$α$ flux. The simulations show that transmission spikes predominantly occur in regions of low gas density. Our results challenge recent observational studies suggesting the origin of these spikes in regions with enhanced radiation. We demonstrate that Ly$α$ transmission spikes are largely a product of the large-scale structure, of which galaxies are biased tracers. △ Less

Submitted 9 January, 2024; originally announced January 2024.

Comments: 10 pages, 8 figures, comments are welcome

Do Mini-halos Affect Cosmic Reionization?

Authors: Nickolay Y. Gnedin

Abstract: The role of unresolved structures ("mini-halos") in determining the consumption of ionizing photons during cosmic reionization remains an unsolved problem in modeling cosmic reionization, despite recent extensive studies with small-box high-resolution simulations by Park et al. and Chan et al., because the small-box studies are not able to fully sample all environments. In this paper these simulat… ▽ More The role of unresolved structures ("mini-halos") in determining the consumption of ionizing photons during cosmic reionization remains an unsolved problem in modeling cosmic reionization, despite recent extensive studies with small-box high-resolution simulations by Park et al. and Chan et al., because the small-box studies are not able to fully sample all environments. In this paper these simulations are combined with large-box simulations from the "Cosmic Reionization On Computers" (CROC) project, allowing one to account for the full range of environments and to produce an estimate for the number of recombinations per hydrogen atom that are missed in large-scale simulations like CROC or Thesan. I find that recombinations in unresolved mini-halos are completely negligible compared to recombinations produced in large-scale cosmic structures and inside more massive, fully resolved halos. Since both Park et al. and Chan et al. studies have severe limitations, the conclusions of this paper may need to be verified with more representative sets of small-box high-resolution simulations. △ Less

Submitted 1 December, 2023; originally announced December 2023.

Comments: Submitted to ApJ

Emergence of the Temperature-Density Relation in the Low Density Intergalactic Medium

Authors: Alexandra Wells, David Robinson, Camille Avestruz, Nickolay Y. Gnedin

Abstract: We examine the evolution of the phase diagram of the low-density intergalactic medium (IGM) during the Epoch of Reionization in simulation boxes with varying reionization histories from the Cosmic Reionization on Computers project. The PDF of gas temperature at fixed density exhibits two clear modes: a warm and cold temperature mode, corresponding to the gas inside and outside of ionized bubbles.… ▽ More We examine the evolution of the phase diagram of the low-density intergalactic medium (IGM) during the Epoch of Reionization in simulation boxes with varying reionization histories from the Cosmic Reionization on Computers project. The PDF of gas temperature at fixed density exhibits two clear modes: a warm and cold temperature mode, corresponding to the gas inside and outside of ionized bubbles. We find that the transition between the two modes is "universal" in the sense that its timing is accurately parameterized by the value of the volume-weighted neutral fraction for any reionization history. This "universality" is more complex than just a reflection of the fact that ionized gas is warm and neutral gas is cold: it holds for the transition at a fixed value of gas density, and gas at different densities transitions from the cold to the warm mode at different values of the neutral fraction, reflecting a non-trivial relationship between the ionization history, the evolving gas density PDF, and the spectrum of ionizing radiation. Furthermore, the "emergence" of the tight temperature-density relation in the warm mode is also approximately "universally" controlled by the volume-weighted neutral fraction for any reionization history. In particular, the "emergence" of the temperature-density relation (as quantified by the rapid decrease in its width) occurs when the neutral fraction is $10^{-4}\lesssim X_\mathrm{HI} \lesssim10^{-3}$ for any reionization history. Our results indicate that the neutral fraction is a primary quantity controlling the various properties of the temperature-density relation, regardless of reionization history. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 7 pages, 10 figures. Submitted to MNRAS

Modeling molecular hydrogen in low metallicity galaxies

Authors: Ava Polzin, Andrey V. Kravtsov, Vadim A. Semenov, Nickolay Y. Gnedin

Abstract: We use a suite of hydrodynamics simulations of the interstellar medium (ISM) within a galactic disk, which include radiative transfer, a non-equilibrium model of molecular hydrogen, and a realistic model for star formation and feedback, to study the structure of the ISM and H$_2$ abundance as a function of local ISM properties. We show that the star formation rate and structure of the ISM are sens… ▽ More We use a suite of hydrodynamics simulations of the interstellar medium (ISM) within a galactic disk, which include radiative transfer, a non-equilibrium model of molecular hydrogen, and a realistic model for star formation and feedback, to study the structure of the ISM and H$_2$ abundance as a function of local ISM properties. We show that the star formation rate and structure of the ISM are sensitive to the metallicity of the gas with a progressively smoother density distribution with decreasing metallicity. In addition to the well-known trend of the HI-H$_2$ transition shifting to higher densities with decreasing metallicity, the maximum achieved molecular fraction in the interstellar medium drops drastically at $Z \lesssim 0.2 \, Z_\odot$ as the formation time of H$_2$ becomes much longer than a typical lifetime of dense regions of the ISM. We present accurate fitting formulae for both volumetric and projected $f_\mathrm{H_2}$ measured on different scales as a function of gas metallicity, UV radiation field, and gas density. We show that when the formulae are applied to the patches in the simulated galaxy the overall molecular gas mass is reproduced to better than a factor of $\lesssim 1.5$ across the entire range of metallicities and scales. We also show that the presented fit is considerably more accurate than any of the previous $f_\mathrm{H_2}$ models and fitting formulae in the low-metallicity regime. The fit can thus be used for modeling molecular gas in low-resolution simulations and semi-analytic models of galaxy formation in the dwarf and high-redshift regimes. △ Less

Submitted 16 October, 2023; originally announced October 2023.

Comments: 17 pages, 12 figures; submitted to ApJ

Exploring the Dependence of Gas Cooling and Heating Functions on the Incident Radiation Field with Machine Learning

Authors: David Robinson, Camille Avestruz, Nickolay Y. Gnedin

Abstract: Gas cooling and heating functions play a crucial role in galaxy formation. But, it is computationally expensive to exactly compute these functions in the presence of an incident radiation field. These computations can be greatly sped up by using interpolation tables of pre-computed values, at the expense of making significant and sometimes even unjustified approximations. Here, we explore the capa… ▽ More Gas cooling and heating functions play a crucial role in galaxy formation. But, it is computationally expensive to exactly compute these functions in the presence of an incident radiation field. These computations can be greatly sped up by using interpolation tables of pre-computed values, at the expense of making significant and sometimes even unjustified approximations. Here, we explore the capacity of machine learning to approximate cooling and heating functions with a generalized radiation field. Specifically, we use the machine learning algorithm XGBoost to predict cooling and heating functions calculated with the photoionization code Cloudy at fixed metallicity, using different combinations of photoionization rates as features. We perform a constrained quadratic fit in metallicity to enable a fair comparison with traditional interpolation methods at arbitrary metallicity. We consider the relative importance of various photoionization rates through both a principal component analysis (PCA) and calculation of SHapley Additive exPlanation (SHAP) values for our XGBoost models. We use feature importance information to select different subsets of rates to use in model training. Our XGBoost models outperform a traditional interpolation approach at each fixed metallicity, regardless of feature selection. At arbitrary metallicity, we are able to reduce the frequency of the largest cooling and heating function errors compared to an interpolation table. We find that the primary bottleneck to increasing accuracy lies in accurately capturing the metallicity dependence. This study demonstrates the potential of machine learning methods such as XGBoost to capture the non-linear behavior of cooling and heating functions. △ Less

Submitted 18 January, 2024; v1 submitted 13 October, 2023; originally announced October 2023.

Comments: 12+4 pages, 6+2 figures. Published in MNRAS. Updated to match MNRAS accepted version

Journal ref: MNRAS, Volume 528, Issue 1, February 2024

Modeling Dust Production, Growth, and Destruction in Reionization-Era Galaxies with the CROC Simulations II: Predicting the Dust Content of High-Redshift Galaxies

Authors: Clarke J. Esmerian, Nickolay Y. Gnedin

Abstract: We model the interstellar dust content of the reionization era with a suite of cosmological, fluid-dynamical simulations of galaxies with stellar masses ranging from $\sim 10^5 - 10^9 M_{\odot}$ in the first $1.2$ billion years of the universe. We use a post-processing method that accounts for dust creation and destruction processes, allowing us to systematically vary the parameters of these proce… ▽ More We model the interstellar dust content of the reionization era with a suite of cosmological, fluid-dynamical simulations of galaxies with stellar masses ranging from $\sim 10^5 - 10^9 M_{\odot}$ in the first $1.2$ billion years of the universe. We use a post-processing method that accounts for dust creation and destruction processes, allowing us to systematically vary the parameters of these processes to test whether dust-dependent observable quantities of galaxies at these epochs could be useful for placing constraints on dust physics. We then forward model observable properties of these galaxies to compare to existing data. We find that we are unable to simultaneously match existing observational constraints with any one set of model parameters. Specifically, the models which predict the largest dust masses $D/Z \gtrsim 0.1$ at $z = 5$ -- because of high assumed production yields and/or efficient growth via accretion in the interstellar medium -- are preferred by constraints on total dust mass and infrared luminosities, but these models produce far too much extinction in the ultraviolet, preventing them from matching observations of $β_{\rm UV}$. To investigate this discrepancy, we analyze the relative spatial distribution of stars and dust as probed by infrared (IR) and ultraviolet (UV) emission, which appear to exhibit overly symmetric morphologies compared to existing data, likely due to the limitations of the stellar feedback model used in the simulations. Our results indicate that the observable properties of the dust distribution in high redshift galaxies are a particularly strong test of stellar feedback. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: Submitted to ApJ, comments welcome

Cosmic Reionization On Computers: Statistics, Physical Properties and Environment of Lyman Limit Systems at $z\sim6$

Authors: Jiawen Fan, Hanjue Zhu, Camille Avestruz, Nickolay Y. Gnedin

Abstract: Lyman limit systems (LLSs) are dense hydrogen clouds with high enough HI column densities to absorb Lyman continuum photons emitted from distant quasars. Their high column densities imply an origin in dense environments; however, the statistics and distribution of LLSs at high redshifts still remain uncertain. In this paper, we use self-consistent radiative transfer cosmological simulations from t… ▽ More Lyman limit systems (LLSs) are dense hydrogen clouds with high enough HI column densities to absorb Lyman continuum photons emitted from distant quasars. Their high column densities imply an origin in dense environments; however, the statistics and distribution of LLSs at high redshifts still remain uncertain. In this paper, we use self-consistent radiative transfer cosmological simulations from the "Cosmic Reionization On Computers" (CROC) project to study the physical properties of LLSs at the tail end of cosmic reionization at $z\sim6$. We generate 3000 synthetic quasar sightlines to obtain a large number of LLS samples in the simulations. In addition, with the high physical fidelity and resolution of CROC, we are able to quantify the association between these LLS samples and nearby galaxies. Our results show that the fraction LLSs spatially associated with nearby galaxies is increasing with the HI column density. Moreover, we find that LLSs that are not near any galaxy typically reside in filamentary structures connecting neighboring galaxies in the intergalactic medium (IGM). This quantification of the distribution and associations of LLSs to large scale structures informs our understanding of the IGM-galaxy connection during the Epoch of Reionization, and provides a theoretical basis for interpreting future observations. △ Less

Submitted 1 March, 2024; v1 submitted 14 December, 2022; originally announced December 2022.

Comments: 11 pages, 10 figures, submitted to ApJ

Journal ref: 2024 ApJ 963 45

Mass-Metallicity Relation during the Epoch of Reionization in the CROC Simulations

Authors: Isaac Noel, Hanjue Zhu, Nickolay Gnedin

Abstract: The low-redshift mass-metallicity relation (MZR) is well studied, but the high-redshift MZR remains difficult to observe. To study the early MZR further, we analyze the Cosmic Reionization on Computers (CROC) simulations with a focus on the MZR from redshifts 5 to 10. We find that, across all redshifts, CROC galaxies exhibit similar stellar-phase and gas-phase MZRs that flatten with higher stellar… ▽ More The low-redshift mass-metallicity relation (MZR) is well studied, but the high-redshift MZR remains difficult to observe. To study the early MZR further, we analyze the Cosmic Reionization on Computers (CROC) simulations with a focus on the MZR from redshifts 5 to 10. We find that, across all redshifts, CROC galaxies exhibit similar stellar-phase and gas-phase MZRs that flatten with higher stellar mass. We attribute this flattening to the inaccurate star formation and feedback modeling in CROC (star formation is overly suppressed in massive CROC galaxies). In addition, we show that the ratio between stellar metallicity and gas metallicity ($Z_*/Z_{gas}$) decreases as stellar age increases, meaning that in CROC galaxies, gas accretion rate is lower than metal production rate. With JWST we will be able to compare our predictions to observations of the Epoch of Reionization and understand better early galaxy formation. △ Less

Submitted 30 October, 2022; originally announced October 2022.

Comments: 5 pages, 5 figures, submitted to ApJ

Cosmic Reionization On Computers: Baryonic Effects on Halo Concentrations During the Epoch of Reionization

Authors: Hanjue Zhu, Nickolay Y. Gnedin

Abstract: Baryons both increase halo concentration through adiabatic contraction and expel mass through feedback processes. However, it is not well understood how the radiation fields prevalent during the epoch of reionization affect the evolution of concentration in dark matter halos. We investigate how baryonic physics during the epoch of reionization modify the structure of dark matter halos in the Cosmi… ▽ More Baryons both increase halo concentration through adiabatic contraction and expel mass through feedback processes. However, it is not well understood how the radiation fields prevalent during the epoch of reionization affect the evolution of concentration in dark matter halos. We investigate how baryonic physics during the epoch of reionization modify the structure of dark matter halos in the Cosmic Reionization On Computers (CROC) simulations. We use two different measures of halo concentration to quantify the effects. We compare concentrations of halos matched between full physics simulations and dark-matter-only simulations with identical initial conditions between $5 \leq z \leq 9$. Baryons in full physics simulations do pull matter towards the center, increasing the maximum circular velocity compared to dark-matter-only simulations. However, their overall effects are much less than if all the baryons were simply centrally concentrated indicating that heating processes efficiently counteract cooling effects. Finally, we show that the baryonic effects on halo concentrations at $z\approx5$ are relatively insensitive to environmental variations of reionization history. These results are pertinent to models of galaxy-halo connection during the epoch of reionization. △ Less

Submitted 29 September, 2022; originally announced September 2022.

Comments: 9 pages, 7 figures, comments are welcome

Compressing the cosmological information in one-dimensional correlations of the Lyman-$α$ forest

Authors: Christian Pedersen, Andreu Font-Ribera, Nickolay Y. Gnedin

Abstract: Observations of the Lyman-$α$ (Ly$α$) forest from spectroscopic surveys such as BOSS/eBOSS, or the ongoing DESI, offer a unique window to study the growth of structure on megaparsec scales. Interpretation of these measurements is a complicated task, requiring hydrodynamical simulations to model and marginalise over the thermal and ionisation state of the intergalactic medium. This complexity has l… ▽ More Observations of the Lyman-$α$ (Ly$α$) forest from spectroscopic surveys such as BOSS/eBOSS, or the ongoing DESI, offer a unique window to study the growth of structure on megaparsec scales. Interpretation of these measurements is a complicated task, requiring hydrodynamical simulations to model and marginalise over the thermal and ionisation state of the intergalactic medium. This complexity has limited the use of Ly$α$ clustering measurements in joint cosmological analyses. In this work we show that the cosmological information content of the 1D power spectrum ($P_\mathrm{1D}$) of the Ly$α$ forest can be compressed into a simple two-parameter likelihood without any significant loss of constraining power. We simulate $P_\mathrm{1D}$ measurements from DESI using hydrodynamical simulations and show that the compressed likelihood is model independent and lossless, recovering unbiased results even in the presence of massive neutrinos or running of the primordial power spectrum. △ Less

Submitted 7 May, 2023; v1 submitted 20 September, 2022; originally announced September 2022.

Comments: Published in ApJ. 17 pages, 7 figures

Journal ref: The Astrophysical Journal, Volume 944, Number 2, 2023

Reconstructing Large-scale Temperature Profiles around $z\sim 6$ Quasars

Authors: Huanqing Chen, Rupert Croft, Nickolay Y. Gnedin

Abstract: High-redshift quasars ionize HeII into HeIII around them, heating the IGM in the process and creating large regions with elevated temperature. In this work, we demonstrate a method based on a convolutional neural network (CNN) to recover the spatial profile for $T_0$, the temperature at the mean cosmic density, in quasar proximity zones. We train the neural network with synthetic spectra drawn fro… ▽ More High-redshift quasars ionize HeII into HeIII around them, heating the IGM in the process and creating large regions with elevated temperature. In this work, we demonstrate a method based on a convolutional neural network (CNN) to recover the spatial profile for $T_0$, the temperature at the mean cosmic density, in quasar proximity zones. We train the neural network with synthetic spectra drawn from a Cosmic Reionization on Computers simulation. We discover that the simple CNN is able to recover the temperature profile with an accuracy of $\approx 1400$ K in an idealized case of negligible observational uncertainties. We test the robustness of the CNN and discover that it is robust against the uncertainties in quasar host halo mass, quasar continuum and ionizing flux. We also find that the CNN has good generality with regard to the hardness of quasar spectra. Saturated pixels pose a bigger problem for accuracy and may downgrade the accuracy to $1700$ K in the outer parts of the proximity zones. Using our method, one could distinguish whether gas is inside or outside the HeIII region created by the quasar. Because the size of the HeIII region is closely related to the total quasar lifetime, this method has great potential in constraining the quasar lifetime on $\sim $Myr timescales. △ Less

Submitted 29 August, 2022; originally announced August 2022.

Comments: 10 pages, 11 figures, comments welcome

Modeling Dust Production, Growth, and Destruction in Reionization-Era Galaxies with the CROC Simulations: Methods and Parameter Exploration

Authors: Clarke J. Esmerian, Nickolay Y. Gnedin

Abstract: We introduce a model for the explicit evolution of interstellar dust in a cosmological galaxy formation simulation. We post-process a simulation from the Cosmic Reionization on Computers project (CROC, Gnedin 2014), integrating an ordinary differential equation for the evolution of the dust-to-gas ratio along pathlines in the simulation sampled with a tracer particle technique. This model incorpor… ▽ More We introduce a model for the explicit evolution of interstellar dust in a cosmological galaxy formation simulation. We post-process a simulation from the Cosmic Reionization on Computers project (CROC, Gnedin 2014), integrating an ordinary differential equation for the evolution of the dust-to-gas ratio along pathlines in the simulation sampled with a tracer particle technique. This model incorporates the effects of dust grain production in asymptotic giant branch (AGB) star winds and supernovae (SN), grain growth due to the accretion of heavy elements from the gas phase of the interstellar medium (ISM), and grain destruction due to thermal sputtering in the high temperature gas of supernova remnants (SNRs). A main conclusion of our analysis is the importance of a carefully chosen dust destruction model, for which different reasonable parameterizations can predict very different values at the $\sim 100$ pc resolution of the ISM in our simulations. We run this dust model on the single most massive galaxy in a 10$h^{-1}$ co-moving Mpc box, which attains a stellar mass of $\sim 2\times10^9 M_{\odot}$ by $z=5$. We find that the model is capable of reproducing dust masses and dust-sensitive observable quantities broadly consistent with existing data from high-redshift galaxies. The total dust mass in the simulated galaxy is somewhat sensitive to parameter choices for the dust model, especially the timescale for grain growth due to accretion in the ISM. Consequently, observable quantities that can constrain galaxy dust masses at these epochs are potentially useful for placing constraints on dust physics. △ Less

Submitted 3 August, 2022; originally announced August 2022.

Comments: Submitted to ApJ, comments welcome!

Modeling Cosmic Reionization

Authors: Nickolay Y. Gnedin, Piero Madau

Abstract: The transformation of cold neutral intergalactic hydrogen into a highly ionized warm plasma marks the end of the cosmic dark ages and the beginning of the age of galaxies. The details of this process reflect the nature of the early sources of radiation and heat, the statistical characteristics of the large-scale structure of the Universe, the thermodynamics and chemistry of cosmic baryons, and the… ▽ More The transformation of cold neutral intergalactic hydrogen into a highly ionized warm plasma marks the end of the cosmic dark ages and the beginning of the age of galaxies. The details of this process reflect the nature of the early sources of radiation and heat, the statistical characteristics of the large-scale structure of the Universe, the thermodynamics and chemistry of cosmic baryons, and the histories of star formation and black hole accretion. A number of massive data sets from new ground- and space-based instruments and facilities over the next decade are poised to revolutionize our understanding of primeval galaxies, the reionization photon budget, the physics of the intergalactic medium (IGM), and the fine-grained properties of hydrogen gas in the "cosmic web". In this review we survey the physics and key aspects of reionization-era modeling and describe the diverse range of computational techniques and tools currently available in this field. △ Less

Submitted 5 October, 2022; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: Accepted to Living Reviews in Computational Astrophysics. This is a "living" review and will be updated regularly. If we missed or misrepresented your recent work, we offer sincere apology and welcome a friendly comment - we will fix our error in the next revision

Revealing the Galaxy-Halo Connection Through Machine Learning

Authors: Ryan Hausen, Brant E. Robertson, Hanjue Zhu, Nickolay Y. Gnedin, Piero Madau, Evan E. Schneider, Bruno Villasenor, Nicole E. Drakos

Abstract: Understanding the connections between galaxy stellar mass, star formation rate, and dark matter halo mass represents a key goal of the theory of galaxy formation. Cosmological simulations that include hydrodynamics, physical treatments of star formation, feedback from supernovae, and the radiative transfer of ionizing photons can capture the processes relevant for establishing these connections. T… ▽ More Understanding the connections between galaxy stellar mass, star formation rate, and dark matter halo mass represents a key goal of the theory of galaxy formation. Cosmological simulations that include hydrodynamics, physical treatments of star formation, feedback from supernovae, and the radiative transfer of ionizing photons can capture the processes relevant for establishing these connections. The complexity of these physics can prove difficult to disentangle and obfuscate how mass-dependent trends in the galaxy population originate. Here, we train a machine learning method called Explainable Boosting Machines (EBMs) to infer how the stellar mass and star formation rate of nearly 6 million galaxies simulated by the Cosmic Reionization on Computers (CROC) project depend on the physical properties of halo mass, the peak circular velocity of the galaxy during its formation history $v_\mathrm{peak}$, cosmic environment, and redshift. The resulting EBM models reveal the relative importance of these properties in setting galaxy stellar mass and star formation rate, with $v_\mathrm{peak}$ providing the most dominant contribution. Environmental properties provide substantial improvements for modeling the stellar mass and star formation rate in only $\lesssim10\%$ of the simulated galaxies. We also provide alternative formulations of EBM models that enable low-resolution simulations, which cannot track the interior structure of dark matter halos, to predict the stellar mass and star formation rate of galaxies computed by high-resolution simulations with detailed baryonic physics. △ Less

Submitted 21 April, 2022; originally announced April 2022.

Comments: 27 pages, 18 figures, Submitted to AAS Journals

Cosmic Reionization on Computers: Physical Origin of Long Dark Gaps in Quasar Absorption Spectra

Authors: Nickolay Y. Gnedin

Abstract: I explore the properties of "dark gaps" - regions in quasar absorption spectra without significant transmission - with several simulations from the Cosmic Reionization On Computers (CROC) project. CROC simulations in largest available boxes (120 cMpc) come close to matching both the distribution of mean opacities and the frequency of dark gaps, but alas not in the same model: the run that matches… ▽ More I explore the properties of "dark gaps" - regions in quasar absorption spectra without significant transmission - with several simulations from the Cosmic Reionization On Computers (CROC) project. CROC simulations in largest available boxes (120 cMpc) come close to matching both the distribution of mean opacities and the frequency of dark gaps, but alas not in the same model: the run that matches the mean opacities fails to contain enough dark gaps and vice versa.:( Never-the-less, the run that matches the dark gap distributions serves as a counter-example to claims in the literature that the dark gap statistics requires a late end to reionization - in that run reionization ends at z=6.7 (likely too early). While multiple factors contribute to the frequency of large dark gaps in the simulations, the primary factor that controls the overall shape of the dark gap distribution is the ionization level in voids - the lowest density regions produce the highest transmission spikes that terminate long gaps. As the result, the dark gap distribution correlates strongly with the fraction of the spectrum above the gap detection threshold, the observed distribution is matched by the simulation in which this fraction is 2%. Hence, the gap distribution by itself does not constrain the timing of reionization. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: Submitted to ApJ

Testing multiflavored ULDM models with SPARC

Authors: Lauren Street, Nickolay Y. Gnedin, L. C. R. Wijewardhana

Abstract: We perform maximum likelihood estimates (MLEs) for single and double flavor ultralight dark matter (ULDM) models using the Spitzer Photometry and Accurate Rotation Curves (SPARC) database. These estimates are compared to MLEs for several commonly used cold dark matter (CDM) models. By comparing various CDM models we find, in agreement with previous studies, that the Burkert and Einasto models tend… ▽ More We perform maximum likelihood estimates (MLEs) for single and double flavor ultralight dark matter (ULDM) models using the Spitzer Photometry and Accurate Rotation Curves (SPARC) database. These estimates are compared to MLEs for several commonly used cold dark matter (CDM) models. By comparing various CDM models we find, in agreement with previous studies, that the Burkert and Einasto models tend to perform better than other commonly used CDM models. We focus on comparisons between the Einasto and ULDM models and analyze cases for which the ULDM particle masses are: free to vary; and fixed. For each of these analyses, we perform fits assuming the soliton and halo profiles are: summed together; and matched at a given radius. When we let the particle masses vary, we find a negligible preference for any particular range of particle masses, within $10^{-25}\,\text{eV}\leq m\leq10^{-19}\,\text{eV}$, when assuming the summed models. For the matched models, however, we find that almost all galaxies prefer particles masses in the range $10^{-23}\,\text{eV}\lesssim m\lesssim10^{-20}\,\text{eV}$. For both double flavor models we find that most galaxies prefer approximately equal particle masses. We find that the summed models give much larger variances with respect to the soliton-halo (SH) relation than the matched models. When the particle masses are fixed, the matched models give median and mean soliton and halo values that fall within the SH relation bounds, for most masses scanned. When the particle masses are fixed in the fitting procedure, we find the best fit results for the particle mass $m=10^{-20.5}\,\text{eV}$ (for the single flavor models) and $m_1=10^{-20.5}\,\text{eV}$, $m_2=10^{-20.2}\,\text{eV}$ for the double flavor, matched model. We discuss how our study will be furthered using a reinforcement learning algorithm. △ Less

Submitted 4 April, 2022; originally announced April 2022.

Comments: 36 pages, 25 figures, 2 appendices

Journal ref: Phys. Rev. D 106, 043007 (2022)

Measuring the Density Fields around Bright Quasars at $z\sim 6$ with XQR-30 Spectra

Authors: Huanqing Chen, Anna-Christina Eilers, Sarah E. I. Bosman, Nickolay Y. Gnedin, Xiaohui Fan, Feige Wang, Jinyi Yang, Valentina D'Odorico, George D. Becker, Manuela Bischetti, Chiara Mazzucchelli, Andrei Mesinger, Andrea Pallottini

Abstract: Measuring the density of the intergalactic medium using quasar sightlines in the epoch of reionization is challenging due to the saturation of Lyman-$α$ absorption. Near a luminous quasar, however, the enhanced radiation creates a proximity zone observable in the quasar spectra where the Lyman-$α$ absorption is not saturated. In this study, we use $10$ high-resolution ($R\gtrsim 10,000$)… ▽ More Measuring the density of the intergalactic medium using quasar sightlines in the epoch of reionization is challenging due to the saturation of Lyman-$α$ absorption. Near a luminous quasar, however, the enhanced radiation creates a proximity zone observable in the quasar spectra where the Lyman-$α$ absorption is not saturated. In this study, we use $10$ high-resolution ($R\gtrsim 10,000$) $z\sim 6$ quasar spectra from the extended XQR-30 sample to measure the density field in the quasar proximity zones. We find a variety of environments within $3$ pMpc distance from the quasars. We compare the observed density cumulative distribution function (CDF) with models from the $\textit{Cosmic Reionization on Computers}$ simulation, and find a good agreement between $1.5$ to $3$pMpc from the quasar. This region is far away from the quasar hosts and hence approaching the mean density of the universe, which allows us to use the CDF to set constraints on the cosmological parameter $σ_8=0.6\pm0.3$. The uncertainty is mainly due to the limited number of high-quality quasar sightlines currently available. Utilizing the more than $>200$ known quasars at $z\gtrsim 6$, this method will allow us in the future to tighten the constraint on $σ_8$ to the percent level. In the region closer to the quasar within $1.5$ pMpc, we find the density is higher than predicted in the simulation by $1.23 \pm 0.17$, suggesting the typical host dark matter halo mass of a bright quasar ($M_{\rm 1450}<-26.5$) at $z\sim 6$ is $\log_{\rm 10} (M/M_\odot)=12.5^{+0.4}_{-0.7}$. △ Less

Submitted 26 October, 2021; originally announced October 2021.

Comments: 16 pages, 11 figures, submitted to ApJ, comments are welcome

Cosmic Reionization on Computers: Evolution of the Flux Power Spectrum

Authors: Nishant Mishra, Nickolay Gnedin

Abstract: We explore the evolution of the flux power spectrum in the Cosmic Reionization On Computers (CROC) simulations. We find that, contrary to some previous studies, the shape of the flux power spectrum is rather insensitive to the timing of reionization. However, the amplitude of the flux power spectrum does strongly evolve with time, and that evolution is almost perfectly correlated with the timing o… ▽ More We explore the evolution of the flux power spectrum in the Cosmic Reionization On Computers (CROC) simulations. We find that, contrary to some previous studies, the shape of the flux power spectrum is rather insensitive to the timing of reionization. However, the amplitude of the flux power spectrum does strongly evolve with time, and that evolution is almost perfectly correlated with the timing of reionization. We show how such correlation can be used in a (futuristic) measurement to determine the redshift of overlap of ionized bubbles. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: Submitted to ApJ

Approximating Density Probability Distribution Functions Across Cosmologies

Authors: Huanqing Chen, Nickolay Y. Gnedin, Philip Mansfield

Abstract: Using a suite of self-similar cosmological simulations, we measure the probability distribution functions (PDFs) of real-space density, redshift-space density, and their geometric mean. We find that the real-space density PDF is well-described by a function of two parameters: $n_s$, the spectral slope, and $σ_L$, the linear rms density fluctuation. For redshift-space density and the geometric mean… ▽ More Using a suite of self-similar cosmological simulations, we measure the probability distribution functions (PDFs) of real-space density, redshift-space density, and their geometric mean. We find that the real-space density PDF is well-described by a function of two parameters: $n_s$, the spectral slope, and $σ_L$, the linear rms density fluctuation. For redshift-space density and the geometric mean of real- and redshift-space densities, we introduce a third parameter, $s_L={\sqrt{\langle(dv^L_{\rm pec}/dr)^2\rangle}}/{H}$. We find that density PDFs for the LCDM cosmology is also well-parameterized by these three parameters. As a result, we are able to use a suite of self-similar cosmological simulations to approximate density PDFs for a range of cosmologies. We make the density PDFs publicly available and provide an analytical fitting formula for them. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: 9 pages, 10 figures, submitted to ApJ, comments welcome

Can Cooling and Heating Functions be Modeled with Homogeneous Radiation Fields?

Authors: David Robinson, Camille Avestruz, Nickolay Y. Gnedin

Abstract: Cooling and heating functions describe how radiative processes impact the thermal state of a gas as a function of its temperature and other physical properties. In a most general case the functions depend on the detailed distributions of ionic species and on the radiation spectrum. Hence, these functions may vary on a very wide range of spatial and temporal scales. In this paper, we explore coolin… ▽ More Cooling and heating functions describe how radiative processes impact the thermal state of a gas as a function of its temperature and other physical properties. In a most general case the functions depend on the detailed distributions of ionic species and on the radiation spectrum. Hence, these functions may vary on a very wide range of spatial and temporal scales. In this paper, we explore cooling and heating functions between $5\leq z \leq10$ in simulated galaxies from the Cosmic Reionization On Computers (CROC) project. We compare three functions. First, the actual cooling and heating rates of hydrodynamic cells as a function of cell temperature. Second, the median cooling and heating functions computed using median interstellar medium (ISM) properties (median ISM). Last, the median of the cooling and heating functions of all gas cells (instantaneous). We find that the median ISM and instantaneous approaches to finding a median cooling and heating function give identical results within the spread due to cell-to-cell variation. However, the actual cooling (heating) rates experienced by the gas at different temperatures in the simulations do not correspond to either summarized cooling (heating) functions. In other words, the thermodynamics of the gas in the simulations cannot be described by a single set of a cooling plus a heating function with a spatially constant radiation field that could be computed with common tools, such as Cloudy. △ Less

Submitted 1 August, 2022; v1 submitted 3 September, 2021; originally announced September 2021.

Comments: 12 pages, 7 figures. Key figure: figure 2. Accepted to ApJ

Journal ref: 2022, ApJ, 936, 50

Reconstruction of the Density Power Spectrum from Quasar Spectra using Machine Learning

Authors: Maria Han Veiga, Xi Meng, Oleg Y. Gnedin, Nickolay Y. Gnedin, Xun Huan

Abstract: We describe a novel end-to-end approach using Machine Learning to reconstruct the power spectrum of cosmological density perturbations at high redshift from observed quasar spectra. State-of-the-art cosmological simulations of structure formation are used to generate a large synthetic dataset of line-of-sight absorption spectra paired with 1-dimensional fluid quantities along the same line-of-sigh… ▽ More We describe a novel end-to-end approach using Machine Learning to reconstruct the power spectrum of cosmological density perturbations at high redshift from observed quasar spectra. State-of-the-art cosmological simulations of structure formation are used to generate a large synthetic dataset of line-of-sight absorption spectra paired with 1-dimensional fluid quantities along the same line-of-sight, such as the total density of matter and the density of neutral atomic hydrogen. With this dataset, we build a series of data-driven models to predict the power spectrum of total matter density. We are able to produce models which yield reconstruction to accuracy of about 1% for wavelengths $k \leq 2 h Mpc^{-1}$, while the error increases at larger $k$. We show the size of data sample required to reach a particular error rate, giving a sense of how much data is necessary to reach a desired accuracy. This work provides a foundation for developing methods to analyse very large upcoming datasets with the next-generation observational facilities. △ Less

Submitted 19 July, 2021; originally announced July 2021.

Comments: 10 pages, 9 figures

A Hydro-Particle-Mesh Code for Efficient and Rapid Simulations of the Intracluster Medium

Authors: Yizhou He, Hy Trac, Nickolay Y. Gnedin

Abstract: We introduce the cosmological HYPER code based on an innovative hydro-particle-mesh (HPM) algorithm for efficient and rapid simulations of gas and dark matter. For the HPM algorithm, we update the approach of Gnedin & Hui (1998) to expand the scope of its application from the lower-density intergalactic medium (IGM) to the higher-density intracluster medium (ICM). While the original algorithm trac… ▽ More We introduce the cosmological HYPER code based on an innovative hydro-particle-mesh (HPM) algorithm for efficient and rapid simulations of gas and dark matter. For the HPM algorithm, we update the approach of Gnedin & Hui (1998) to expand the scope of its application from the lower-density intergalactic medium (IGM) to the higher-density intracluster medium (ICM). While the original algorithm tracks only one effective particle species, the updated version separately tracks the gas and dark matter particles as they do not exactly trace each other on small scales. For the approximate hydrodynamics solver, the pressure term in the gas equations of motion is calculated using robust physical models. In particular, we use a dark matter halo model, ICM pressure profile, and IGM temperature-density relation, all of which can be systematically varied for parameter-space studies. We show that the HYPER simulation results are in good agreement with the halo model expectations for the density, temperature, and pressure radial profiles. Simulated galaxy cluster scaling relations for Sunyaev-Zel'dovich (SZ) and X-ray observables are also in good agreement with mean predictions, with scatter comparable to that found in hydrodynamic simulations. HYPER also produces lightcone catalogs of dark matter halos and full-sky tomographic maps of the lensing convergence, SZ effect, and X-ray emission. These simulation products are useful for testing data analysis pipelines, generating training data for machine learning, understanding selection and systematic effects, and for interpreting astrophysical and cosmological constraints. △ Less

Submitted 9 July, 2021; originally announced July 2021.

Comments: 24 pages, 14 figures

Gravitational Self-force Errors of Poisson Solvers on Adaptively Refined Meshes

Authors: Hanjue Zhu, Nickolay Y. Gnedin

Abstract: An error in the gravitational force that the source of gravity induces on itself (a self-force error) violates both the conservation of linear momentum and the conservation of energy. If such errors are present in a self-gravitating system and are not sufficiently random to average out, the obtained numerical solution will become progressively more unphysical with time: the system will acquire or… ▽ More An error in the gravitational force that the source of gravity induces on itself (a self-force error) violates both the conservation of linear momentum and the conservation of energy. If such errors are present in a self-gravitating system and are not sufficiently random to average out, the obtained numerical solution will become progressively more unphysical with time: the system will acquire or lose momentum and energy due to numerical effects. In this paper, we demonstrate how self-force errors can arise in the case where self-gravity is solved on an adaptively refined mesh when the refinement is nonuniform. We provide the analytical expression for the self-force error and numerical examples that demonstrate such self-force errors in idealized settings. We also show how these errors can be corrected to an arbitrary order by straightforward addition of correction terms at the refinement boundaries. △ Less

Submitted 20 April, 2021; originally announced April 2021.

Comments: 8 pages, 5 figures, ApJS accepted

Spatial Decorrelation of Young Stars and Dense Gas as a Probe of the Star Formation-Feedback Cycle in Galaxies

Authors: Vadim A. Semenov, Andrey V. Kravtsov, Nickolay Y. Gnedin

Abstract: The spatial decorrelation of dense molecular gas and young stars observed on $\lesssim 1$ kiloparsec scales in nearby galaxies indicates rapid dispersal of star-forming regions by stellar feedback. We explore the sensitivity of this decorrelation to different processes controlling the structure of the interstellar medium, the abundance of molecular gas, star formation, and feedback in a suite of s… ▽ More The spatial decorrelation of dense molecular gas and young stars observed on $\lesssim 1$ kiloparsec scales in nearby galaxies indicates rapid dispersal of star-forming regions by stellar feedback. We explore the sensitivity of this decorrelation to different processes controlling the structure of the interstellar medium, the abundance of molecular gas, star formation, and feedback in a suite of simulations of an isolated dwarf galaxy with structural properties similar to NGC300 that self-consistently model radiative transfer and molecular chemistry. Our fiducial simulation reproduces the magnitude of decorrelation and its scale dependence measured in NGC300, and we show that this agreement is due to different aspects of feedback, including H$_2$ dissociation, gas heating by the locally variable UV field, early mechanical feedback, and supernovae. In particular, early radiative and mechanical feedback affects the correlation on $\lesssim 100$ pc scales, while supernovae play a significant role on $\gtrsim 100$ pc scales. The correlation is also sensitive to the choice of the local star formation efficiency per freefall time, $ε_{\rm ff}$, which provides a strong observational constraint on $ε_{\rm ff}$ when the global star formation rate is independent of its value. Finally, we explicitly show that the degree of correlation between the peaks of molecular gas and star formation density is directly related to the distribution of the lifetimes of star-forming regions. △ Less

Submitted 13 September, 2021; v1 submitted 24 March, 2021; originally announced March 2021.

Comments: 21 pages + appendix, 17 figures; published in ApJ; the animated version of Figure 1 is available in the online version of the article and at https://vadimsemenov.com/research/ngc300

Journal ref: ApJ 918, 13 (2021)

Recovering Density Fields inside Quasar Proximity Zones at $z\sim 6$

Authors: Huanqing Chen, Nickolay Y. Gnedin

Abstract: The matter density field at $z\sim 6$ is very challenging to probe. One of the traditional probes of the low density IGM that works successfully at lower redshifts is the Lyman-alpha forest in quasar spectra. However, at the end of reionization, the residual neutral hydrogen usually creates saturated absorption, thus much of the information about the gas density is lost. Luckily, in a quasar proxi… ▽ More The matter density field at $z\sim 6$ is very challenging to probe. One of the traditional probes of the low density IGM that works successfully at lower redshifts is the Lyman-alpha forest in quasar spectra. However, at the end of reionization, the residual neutral hydrogen usually creates saturated absorption, thus much of the information about the gas density is lost. Luckily, in a quasar proximity zone, the ionizing radiation is exceptionally intense, thus creating a large region with non-zero transmitted flux. In this study we use the synthetic spectra from simulations to investigate how to recover the density fluctuations inside the quasar proximity zones. We show that, under ideal conditions, the density can be recovered accurately with a small scatter. We also discuss how systematics such as the quasar continuum fitting and reionization models affect the results. This study shows that by analyzing the absorption features inside quasar proximity zones we can potentially constrain quasar properties and the environments they reside in. △ Less

Submitted 26 September, 2021; v1 submitted 27 January, 2021; originally announced January 2021.

Comments: 12 pages, 12 figures, comments welcome

Natively Periodic Fast Multipole Method: Approximating the Optimal Green Function

Authors: Nickolay Y. Gnedin

Abstract: The Fast Multipole Method (FMM) obeys periodic boundary conditions "natively" if it uses a periodic Green function for computing the multipole expansion in the interaction zone of each FMM oct-tree node. One can define the "optimal" Green function for such a method that results in the numerical solution that converges to the equivalent Particle-Mesh solution in the limit of sufficiently high order… ▽ More The Fast Multipole Method (FMM) obeys periodic boundary conditions "natively" if it uses a periodic Green function for computing the multipole expansion in the interaction zone of each FMM oct-tree node. One can define the "optimal" Green function for such a method that results in the numerical solution that converges to the equivalent Particle-Mesh solution in the limit of sufficiently high order of multipoles. A discrete functional equation for the optimal Green function can be derived, but is not practically useful as methods for its solution are not known. Instead, this paper presents an approximation for the optimal Green function that is accurate to better than 1e-3 in LMAX norm and 1e-4 in L2 norm for practically useful multipole counts. Such an approximately optimal Green function offers a practical way for implementing FMM with periodic boundary conditions "natively", without the need to compute lattice sums or to rely on hybrid FMM-PM approaches. △ Less

Submitted 11 January, 2021; v1 submitted 13 November, 2020; originally announced November 2020.

Comments: Accepted for publication by ApJ

The Distribution and Evolution of Quasar Proximity Zone Sizes

Authors: Huanqing Chen, Nickolay Y. Gnedin

Abstract: In this paper, we study the sizes of quasar proximity zones with synthetic quasar absorption spectra obtained by post-processing a Cosmic Reionization On Computers (CROC) simulation. CROC simulations have both relatively large box sizes and high spacial resolution, allowing us to resolve Lyman limit systems, which are crucial for modeling the quasar absorption spectra. We find that before reioniza… ▽ More In this paper, we study the sizes of quasar proximity zones with synthetic quasar absorption spectra obtained by post-processing a Cosmic Reionization On Computers (CROC) simulation. CROC simulations have both relatively large box sizes and high spacial resolution, allowing us to resolve Lyman limit systems, which are crucial for modeling the quasar absorption spectra. We find that before reionization most quasar proximity zone sizes grow steadily for $\sim 10$ Myr, while after reionization they grow rapidly but only for $\sim 0.1$ Myr. We also find a slow growth of $R_{\rm obs}$ with decreasing turn-on redshift. In addition, we find that $\sim 1-2\%$ of old quasars ($30$ Myr old) display extremely small proximity zone sizes ($<1$ proper Mpc), of which the vast majority are due to the occurrence of a damped Ly$α$ absorber (DLA) or a Lyman limit system (LLS) along the line of sight. These DLAs and LLSs are contaminated with metal, which offers a way to distinguish them from the normal proximity zones of young quasars. △ Less

Submitted 25 March, 2021; v1 submitted 11 August, 2020; originally announced August 2020.

Comments: 18 pages, 15 figures, comments welcome!

Reionization with galaxies and active galactic nuclei

Authors: Pratika Dayal, Marta Volonteri, Tirthankar Roy Choudhury, Raffaella Schneider, Maxime Trebitsch, Nickolay Y. Gnedin, Hakim Atek, Michaela Hirschmann, Amy Reines

Abstract: In this work we investigate the properties of the sources that reionized the intergalactic medium (IGM) in the high-redshift Universe. Using a semi-analytical model aimed at reproducing galaxies and black holes in the first 1.5 Gyr of the Universe, we revisit the relative role of star formation and black hole accretion in producing ionizing photons that can escape into the IGM. Both star formation… ▽ More In this work we investigate the properties of the sources that reionized the intergalactic medium (IGM) in the high-redshift Universe. Using a semi-analytical model aimed at reproducing galaxies and black holes in the first 1.5 Gyr of the Universe, we revisit the relative role of star formation and black hole accretion in producing ionizing photons that can escape into the IGM. Both star formation and black hole accretion are regulated by supernova feedback, resulting in black hole accretion being stunted in low-mass halos. We explore a wide range of combinations for the escape fraction of ionizing photons (redshift-dependent, constant and scaling with stellar mass) from both star formation ($\langle f_{\rm esc}^{\rm sf} \rangle$) and AGN ($f_{\rm esc}^{\rm bh}$) to find: (i) the ionizing budget is dominated by stellar radiation from low stellar mass ($M_*<10^9 {\rm M_\odot}$ ) galaxies at $z>6$ with the AGN contribution (driven by $M_{bh}>10^6 {\rm M_\odot}$ black holes in $M_* > 10^9 {\rm M_\odot}$ galaxies) dominating at lower redshifts; (ii) AGN only contribute $10-25\%$ to the cumulative ionizing emissivity by $z=4$ for the models that match the observed reionization constraints; (iii) if the stellar mass dependence of $\langle f_{\rm esc}^{\rm sf} \rangle$ is shallower than $f_{\rm esc}^{\rm bh}$, at $z<7$ a transition stellar mass exists above which AGN dominate the escaping ionizing photon production rate; (iv) the transition stellar mass decreases with decreasing redshift. While AGN dominate the escaping emissivity above the knee of the stellar mass function at $z \sim 6.8$, they take-over at stellar masses that are a tenth of the knee mass by $z=4$. △ Less

Submitted 21 April, 2020; v1 submitted 16 January, 2020; originally announced January 2020.

Comments: Accepted to MNRAS

Cosmic Reionization On Computers: The Galaxy-Halo Connection between $5 \leq z \leq10$

Authors: Hanjue Zhu, Camille Avestruz, Nickolay Y. Gnedin

Abstract: We explore the connection between the stellar component of galaxies and their host halos during the epoch of reionization ($5 \leq z\leq10$) using the CROC (Cosmic Reionization on Computers) simulations. We compare simulated galaxies with observations and find that CROC underpredicts the abundance of luminous galaxies when compared to observed UV luminosity functions, and analogously the most mass… ▽ More We explore the connection between the stellar component of galaxies and their host halos during the epoch of reionization ($5 \leq z\leq10$) using the CROC (Cosmic Reionization on Computers) simulations. We compare simulated galaxies with observations and find that CROC underpredicts the abundance of luminous galaxies when compared to observed UV luminosity functions, and analogously the most massive galaxies when compared to observed stellar mass functions. We can trace the deficit of star formation to high redshifts, where the slope of the star formation rate to stellar mass relation is consistent with observations, but the normalization is systematically low. This results in a star formation rate density and stellar mass density that is systematically offset from observations. However, the less luminous or lower stellar mass objects have luminosities and stellar masses that agree fairly well with observational data. We explore the stellar-to-halo mass ratio, a key quantity that is difficult to measure at high redshifts and that models do not consistently predict. In CROC, the stellar-to-halo mass ratio {\it decreases} with redshift, a trend opposite to some abundance matching studies. These discrepancies uncover where future effort should be focused in order to improve the fidelity of modeling cosmic reionization. We also compare the CROC galaxy bias with observational measurements using Lyman-Break Galaxy (LBG) samples. The good agreement of simulation and data shows that the clustering of dark matter halos is properly captured in CROC. △ Less

Submitted 16 April, 2021; v1 submitted 7 January, 2020; originally announced January 2020.

Comments: 11 pages, 9 figures, comments are welcome

The Effect of Warm Dark Matter on Early Star Formation Histories of Massive Galaxies: Predictions from the CROC Simulations

Authors: Clarke J. Esmerian, Nickolay Y. Gnedin

Abstract: Several massive ($M_{*} > 10^8 M_{\odot}$), high-redshift ($z = 8-10$) galaxies have recently been discovered to contain stars with ages of several hundred million years, pushing the onset of star formation in these galaxies back to $z\sim15$. The very existence of stars formed so early may serve as a test for cosmological models with suppressed small-scale power (and, hence, late formation of cos… ▽ More Several massive ($M_{*} > 10^8 M_{\odot}$), high-redshift ($z = 8-10$) galaxies have recently been discovered to contain stars with ages of several hundred million years, pushing the onset of star formation in these galaxies back to $z\sim15$. The very existence of stars formed so early may serve as a test for cosmological models with suppressed small-scale power (and, hence, late formation of cosmic structure). We explore the ages of the oldest stars in numerical simulations from the Cosmic Reionization On Computers (CROC) project with cold dark matter (CDM) and two warm dark matter (WDM) cosmologies with 3 and 6 keV particles. There are statistically significant differences of $\sim 5\;{\rm Myr}$ between average stellar ages of massive galaxies in CDM and 3 keV WDM, while CDM and 6 keV WDM are statistically indistinguishable. Even this 5 Myr difference, however, is much less than current observational uncertainties on the stellar population properties of high-redshift galaxies. The age distributions of all galaxies in all cosmologies fail to produce a substantial Balmer break, although uncertainties in dust attenuation are a potentially significant factor. Finally, we assess the convergence of our simulation predictions and find that the systematic uncertainties on individual galaxy properties are comparable to the differences between cosmologies, suggesting these differences may not be numerically robust. △ Less

Submitted 31 May, 2021; v1 submitted 22 September, 2019; originally announced September 2019.

Comments: Accepted in ApJ, update to match published version

Journal ref: ApJ 910, 117 (2021)

Hierarchical Particle-Mesh: an FFT-accelerated Fast Multipole Method

Authors: Nickolay Y. Gnedin

Abstract: I describe a modification to the original Fast Multipole Method (FMM) of Greengard & Rokhlin that approximates the gravitation field of an FMM cell as a small uniform grid (a "gridlet") of effective masses. The effective masses on a gridlet are set from the requirement that the multipole moments of the FMM cells are reproduced exactly, hence preserving the accuracy of the gravitational field repre… ▽ More I describe a modification to the original Fast Multipole Method (FMM) of Greengard & Rokhlin that approximates the gravitation field of an FMM cell as a small uniform grid (a "gridlet") of effective masses. The effective masses on a gridlet are set from the requirement that the multipole moments of the FMM cells are reproduced exactly, hence preserving the accuracy of the gravitational field representation. The calculation of the gravitational field from a multipole expansion can then be computed for all multipole orders simultaneously, with a single Fast Fourier Transform, significantly reducing the computational cost at a given value of the required accuracy. The described approach belongs to the class of "kernel independent" variants of the FMM method and works with any Green function. △ Less

Submitted 25 June, 2019; originally announced June 2019.

Comments: Accepted for publication in the ApJ

Cosmology with the Highly Redshifted 21cm Line

Authors: Adrian Liu, James Aguirre, Yacine Ali-Haimoud, Marcelo Alvarez, Adam Beardsley, George Becker, Judd Bowman, Patrick Breysse, Volker Bromm, Philip Bull, Jack Burns, Isabella P. Carucci, Tzu-Ching Chang, Xuelei Chen, Hsin Chiang, Joanne Cohn, David DeBoer, Joshua Dillon, Olivier Doré, Cora Dvorkin, Anastasia Fialkov, Steven Furlanetto, Nick Gnedin, Bryna Hazelton, Jacqueline Hewitt , et al. (25 additional authors not shown)

Abstract: In addition to being a probe of Cosmic Dawn and Epoch of Reionization astrophysics, the 21cm line at $z>6$ is also a powerful way to constrain cosmology. Its power derives from several unique capabilities. First, the 21cm line is sensitive to energy injections into the intergalactic medium at high redshifts. It also increases the number of measurable modes compared to existing cosmological probes… ▽ More In addition to being a probe of Cosmic Dawn and Epoch of Reionization astrophysics, the 21cm line at $z>6$ is also a powerful way to constrain cosmology. Its power derives from several unique capabilities. First, the 21cm line is sensitive to energy injections into the intergalactic medium at high redshifts. It also increases the number of measurable modes compared to existing cosmological probes by orders of magnitude. Many of these modes are on smaller scales than are accessible via the CMB, and moreover have the advantage of being firmly in the linear regime (making them easy to model theoretically). Finally, the 21cm line provides access to redshifts prior to the formation of luminous objects. Together, these features of 21cm cosmology at $z>6$ provide multiple pathways toward precise cosmological constraints. These include the "marginalizing out" of astrophysical effects, the utilization of redshift space distortions, the breaking of CMB degeneracies, the identification of signatures of relative velocities between baryons and dark matter, and the discovery of unexpected signs of physics beyond the $Λ$CDM paradigm at high redshifts. △ Less

Submitted 14 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to Decadal 2020 survey

Astro2020 Science White Paper: First Stars and Black Holes at Cosmic Dawn with Redshifted 21-cm Observations

Authors: Jordan Mirocha, Daniel Jacobs, Josh Dillon, Steve Furlanetto, Jonathan Pober, Adrian Liu, James Aguirre, Yacine Ali-Haïmoud, Marcelo Alvarez, Adam Beardsley, George Becker, Judd Bowman, Patrick Breysse, Volker Bromm, Jack Burns, Xuelei Chen, Tzu-Ching Chang, Hsin Chiang, Joanne Cohn, David DeBoer, Cora Dvorkin, Anastasia Fialkov, Nick Gnedin, Bryna Hazelton, Masui Kiyoshi , et al. (17 additional authors not shown)

Abstract: The "cosmic dawn" refers to the period of the Universe's history when stars and black holes first formed and began heating and ionizing hydrogen in the intergalactic medium (IGM). Though exceedingly difficult to detect directly, the first stars and black holes can be constrained indirectly through measurements of the cosmic 21-cm background, which traces the ionization state and temperature of int… ▽ More The "cosmic dawn" refers to the period of the Universe's history when stars and black holes first formed and began heating and ionizing hydrogen in the intergalactic medium (IGM). Though exceedingly difficult to detect directly, the first stars and black holes can be constrained indirectly through measurements of the cosmic 21-cm background, which traces the ionization state and temperature of intergalactic hydrogen gas. In this white paper, we focus on the science case for such observations, in particular those targeting redshifts z $\gtrsim$ 10 when the IGM is expected to be mostly neutral. 21-cm observations provide a unique window into this epoch and are thus critical to advancing first star and black hole science in the next decade. △ Less

Submitted 14 March, 2019; originally announced March 2019.

Comments: White paper submitted to the Astro2020 Decadal Survey

Astro2020 Science White Paper: Insights Into the Epoch of Reionization with the Highly-Redshifted 21-cm Line

Authors: Steven Furlanetto, Chris L. Carilli, Jordan Mirocha, James Aguirre, Yacine Ali-Haimoud, Marcelo Alvarez, Adam Beardsley, George Becker, Judd D. Bowman, Patrick Breysse, Volker Bromm, Philip Bull, Jack Burns, Isabella P. Carucci, Tzu-Ching Chang, Xuelei Chen, Hsin Chiang, Joanne Cohn, Frederick Davies, David DeBoer, Joshua Dillon, Olivier Doré, Cora Dvorkin, Anastasia Fialkov, Nick Gnedin , et al. (25 additional authors not shown)

Abstract: The epoch of reionization, when photons from early galaxies ionized the intergalactic medium about a billion years after the Big Bang, is the last major phase transition in the Universe's history. Measuring the characteristics of the transition is important for understanding early galaxies and the cosmic web and for modeling dwarf galaxies in the later Universe. But such measurements require probe… ▽ More The epoch of reionization, when photons from early galaxies ionized the intergalactic medium about a billion years after the Big Bang, is the last major phase transition in the Universe's history. Measuring the characteristics of the transition is important for understanding early galaxies and the cosmic web and for modeling dwarf galaxies in the later Universe. But such measurements require probes of the intergalactic medium itself. Here we describe how the 21-cm line of neutral hydrogen provides a powerful probe of the reionization process and therefore important constraints on both the galaxies and intergalactic absorbers at that time. While existing experiments will make precise statistical measurements over the next decade, we argue that improved 21-cm analysis techniques - allowing imaging of the neutral gas itself - as well as improved theoretical models, are crucial for testing our understanding of this important era. △ Less

Submitted 18 March, 2019; v1 submitted 14 March, 2019; originally announced March 2019.

Comments: 8 pages, 4 figures

Astro2020 Science White Paper: Synergies Between Galaxy Surveys and Reionization Measurements

Authors: Steven Furlanetto, Adam Beardsley, Chris L. Carilli, Jordan Mirocha, James Aguirre, Yacine Ali-Haimoud, Marcelo Alvarez, George Becker, Judd D. Bowman, Patrick Breysse, Volker Bromm, Philip Bull, Jack Burns, Isabella P. Carucci, Tzu-Ching Chang, Hsin Chiang, Joanne Cohn, Frederick Davies, David DeBoer, Mark Dickinson, Joshua Dillon, Olivier Doré, Cora Dvorkin, Anastasia Fialkov, Steven Finkelstein , et al. (25 additional authors not shown)

Abstract: The early phases of galaxy formation constitute one of the most exciting frontiers in astrophysics. It is during this era that the first luminous sources reionize the intergalactic medium - the moment when structure formation affects every baryon in the Universe. Here we argue that we will obtain a complete picture of this era by combining observations of galaxies with direct measurements of the r… ▽ More The early phases of galaxy formation constitute one of the most exciting frontiers in astrophysics. It is during this era that the first luminous sources reionize the intergalactic medium - the moment when structure formation affects every baryon in the Universe. Here we argue that we will obtain a complete picture of this era by combining observations of galaxies with direct measurements of the reionization process: the former will provide a detailed understanding of bright sources, while the latter will constrain the (substantial) faint source population. We further describe how optimizing the comparison of these two measurements requires near-infrared galaxy surveys covering large volumes and retaining redshift information and also improvements in 21-cm analysis, moving those experiments into the imaging regime. △ Less

Submitted 18 March, 2019; v1 submitted 14 March, 2019; originally announced March 2019.

Comments: 9 pages, 3 figures

Cosmic Reionization On Computers: Reionization Histories of Present-day Galaxies

Authors: Hanjue Zhu, Camille Avestruz, Nickolay Y. Gnedin

Abstract: We examine the reionization history of present-day galaxies by explicitly tracing the building blocks of halos from the Cosmic Reionization On Computers project. We track dark matter particles that belong to $z=0$ halos to trace the neutral fractions at corresponding positions during rapid global reionization. The resulting particle reionization histories allow us to explore different definitions… ▽ More We examine the reionization history of present-day galaxies by explicitly tracing the building blocks of halos from the Cosmic Reionization On Computers project. We track dark matter particles that belong to $z=0$ halos to trace the neutral fractions at corresponding positions during rapid global reionization. The resulting particle reionization histories allow us to explore different definitions of a halo's reionization redshift and to account for the neutral content of the interstellar medium. Consistent with previous work, we find a systematic trend of reionization redshift with mass - present day halos with higher masses have earlier reionization times. Finally, we quantify the spread of reionization times within each halo, which also has a mass dependence. △ Less

Submitted 28 November, 2019; v1 submitted 14 March, 2019; originally announced March 2019.

Comments: 7 pages, 5 figures, Accepted to ApJ

Constraining the Tail-End of Reionization Using Lyman-$α$ Transmission Spikes

Authors: Enrico Garaldi, Nickolay Gnedin, Piero Madau

Abstract: We investigate Lyman-$α$ transmission spikes at $z > 5$ in synthetic quasar spectra, and discuss their connection to the properties of the intergalactic medium and their ability to constrain reionization models. We use state-of-the-art radiation-hydrodynamic simulations from the Cosmic Reionization On Computers series to predict the number of transmission spikes as a function of redshift, both in… ▽ More We investigate Lyman-$α$ transmission spikes at $z > 5$ in synthetic quasar spectra, and discuss their connection to the properties of the intergalactic medium and their ability to constrain reionization models. We use state-of-the-art radiation-hydrodynamic simulations from the Cosmic Reionization On Computers series to predict the number of transmission spikes as a function of redshift, both in the ideal case of infinite spectral resolution and in a realistic observational setting. Transmission spikes are produced in highly-ionized underdense regions located in the vicinity of UV sources. We find that most of the predicted spikes are unresolved by current observations, and show that our mock spectra are consistent with observations of the quasar ULAS J1120+0641 in about 15% of the realizations. The spike height correlates with both the gas density and the ionized fraction, but the former link is erased when synthetic spectra are smoothed to realistically achievable spectral resolutions. There exists a linear relationship between spike width and the extent of the associated underdense region, with a slope that is redshift-dependent. In agreement with observations, the spike transmitted flux is suppressed at small distance from bright galaxies as these reside in overdense regions. We argue that this anti-correlation can be used to constrain large-scale density modes. △ Less

Submitted 22 March, 2019; v1 submitted 20 February, 2019; originally announced February 2019.

Comments: 16 pages, 12 figures, accepted by ApJ

What Sets the Slope of the Molecular Kennicutt-Schmidt Relation?

Authors: Vadim A. Semenov, Andrey V. Kravtsov, Nickolay Y. Gnedin

Abstract: The surface densities of molecular gas, $Σ_{\rm H_2}$, and the star formation rate (SFR), $\dotΣ_\star$, correlate almost linearly on kiloparsec scales in observed star-forming (non-starburst) galaxies. We explore the origin of the linear slope of this correlation using a suite of isolated $L_\star$ galaxy simulations. We show that in simulations with efficient feedback, the slope of the… ▽ More The surface densities of molecular gas, $Σ_{\rm H_2}$, and the star formation rate (SFR), $\dotΣ_\star$, correlate almost linearly on kiloparsec scales in observed star-forming (non-starburst) galaxies. We explore the origin of the linear slope of this correlation using a suite of isolated $L_\star$ galaxy simulations. We show that in simulations with efficient feedback, the slope of the $\dotΣ_\star$-$Σ_{\rm H_2}$ relation on kiloparsec scales is insensitive to the slope of the $\dotρ_\star$-$ρ$ relation assumed at the resolution scale. We also find that the slope on kiloparsec scales depends on the criteria used to identify star-forming gas, with a linear slope arising in simulations that identify star-forming gas using a virial parameter threshold. This behavior can be understood using a simple theoretical model based on conservation of interstellar gas mass as the gas cycles between atomic, molecular, and star-forming states under the influence of feedback and dynamical processes. In particular, we show that the linear slope emerges when feedback efficiently regulates and stirs the evolution of dense, molecular gas. We show that the model also provides insights into the likely origin of the relation between the SFR and molecular gas in real galaxies on different scales. △ Less

Submitted 18 December, 2018; v1 submitted 19 September, 2018; originally announced September 2018.

Comments: 14 pages + appendix, 15 figures; accepted for publication in the Astrophysical Journal

Journal ref: ApJ 870, 79 (2019)

Constraints On Duty Cycles Of Quasars At Z $\sim$ 6

Authors: Huanqing Chen, Nickolay Y. Gnedin

Abstract: We study the mass of quasar-hosting dark matter halos at z $\sim$ 6 and further constrain the fraction of dark matter halos hosting an active quasar $f_{on}$ and the quasar beaming angle $i_{\rm max}$ using observations of CII lines in the literature. We make assumptions that (1) more massive halos host brighter quasars, (2) a fraction of the halos host active quasars with a certain beaming angle,… ▽ More We study the mass of quasar-hosting dark matter halos at z $\sim$ 6 and further constrain the fraction of dark matter halos hosting an active quasar $f_{on}$ and the quasar beaming angle $i_{\rm max}$ using observations of CII lines in the literature. We make assumptions that (1) more massive halos host brighter quasars, (2) a fraction of the halos host active quasars with a certain beaming angle, (3) cold gas in galaxies has rotational velocity $V_{\rm circ}=αV_{\rm max}$, and that (4) quasars point randomly on the sky. We find that for a choice of specific $α\gtrsim 1$, the most likely solution has $f_{\rm on} < 0.01$, corresponding to a small duty cycle of quasar activity. However, if we marginalize over $α$, for some choices of a prior a second solution with $f_{\rm on}=1$ appears. Overall, our the constraints are highly sensitive to $α$ and hence inconclusive. Stronger constraints on $f_{\rm on}$ can be made if we better understand the dynamics of cold gas in these galaxies. △ Less

Submitted 5 September, 2018; originally announced September 2018.

Criterions for retrograde rotation of accreting black holes

Authors: A. G. Mikhailov, M. Yu. Piotrovich, Yu. N. Gnedin, T. M. Natsvlishvili, S. D. Buliga

Abstract: Rotating supermassive black holes produce jets and their origin is connected to magnetic field that is generated by accreting matter flow. There is a point of view that electromagnetic fields around rotating black holes are brought to the hole by accretion. In this situation the prograde accreting disks produce weaker large-scale black hole threading magnetic fields, implying weaker jets that in r… ▽ More Rotating supermassive black holes produce jets and their origin is connected to magnetic field that is generated by accreting matter flow. There is a point of view that electromagnetic fields around rotating black holes are brought to the hole by accretion. In this situation the prograde accreting disks produce weaker large-scale black hole threading magnetic fields, implying weaker jets that in retrograde regimes. The basic goal of this paper is to find the best candidates for retrograde accreting systems in observed active galactic nuclei. We show that active galactic nuclei with low Eddington ratio are really the best candidates for retrograde systems. This conclusion is obtained for kinetically dominated FRII radio galaxies, flat spectrum radio loud narrow line Seyfert I galaxies and a number of nearby galaxies. Our conclusion is that the best candidates for retrograde systems are the noticeable population of active galactic nuclei in the Universe. This result corresponds to the conclusion that in the merging process the interaction of merging black holes with a retrograde circumbinary disk is considerably more effective for shrinking the binary system. △ Less

Submitted 9 March, 2018; originally announced March 2018.

Comments: 5 pages

How galaxies form stars: the connection between local and global star formation in galaxy simulations

Authors: Vadim A. Semenov, Andrey V. Kravtsov, Nickolay Y. Gnedin

Abstract: Using a suite of isolated $L_\star$ galaxy simulations, we show that global depletion times and star-forming gas mass fractions in simulated galaxies exhibit systematic and well-defined trends as a function of the local star formation efficiency per freefall time, $ε_{\rm ff}$, strength of stellar feedback, and star formation threshold. We demonstrate that these trends can be reproduced and explai… ▽ More Using a suite of isolated $L_\star$ galaxy simulations, we show that global depletion times and star-forming gas mass fractions in simulated galaxies exhibit systematic and well-defined trends as a function of the local star formation efficiency per freefall time, $ε_{\rm ff}$, strength of stellar feedback, and star formation threshold. We demonstrate that these trends can be reproduced and explained by a simple physical model of global star formation in galaxies. Our model is based on mass conservation and the idea of gas cycling between star-forming and non-star-forming states on certain characteristic time scales under the influence of dynamical and feedback processes. Both the simulation results and our model predictions exhibit two limiting regimes with rather different dependencies of global galactic properties on the local parameters. When $ε_{\rm ff}$ is small and feedback is inefficient, the total star-forming mass fraction, $f_{\rm sf}$, is independent of $ε_{\rm ff}$ and the global depletion time, $τ_{\rm dep}$, scales inversely with $ε_{\rm ff}$. When $ε_{\rm ff}$ is large or feedback is very efficient, these trends are reversed: $f_{\rm sf} \propto ε_{\rm ff}^{-1}$ and $τ_{\rm dep}$ is independent of $ε_{\rm ff}$ but scales linearly with the feedback strength. We also compare our results with the observed depletion times and mass fractions of star-forming and molecular gas and show that they provide complementary constraints on $ε_{\rm ff}$ and the feedback strength. We show that useful constraints on $ε_{\rm ff}$ can also be obtained using measurements of the depletion time and its scatter on different spatial scales. △ Less

Submitted 13 June, 2018; v1 submitted 28 February, 2018; originally announced March 2018.

Comments: 17 pages + appendix, 13 figures, 1 table; accepted for publication in ApJ

Journal ref: ApJ 861, 4 (2018)

Enforcing the Courant-Friedrichs-Lewy Condition in Explicitly Conservative Local Time Stepping Schemes

Authors: Nickolay Y. Gnedin, Vadim A. Semenov, Andrey V. Kravtsov

Abstract: An optimally efficient explicit numerical scheme for solving fluid dynamics equations, or any other parabolic or hyperbolic system of partial differential equations, should allow local regions to advance in time with their own, locally constrained time steps. However, such a scheme can result in violation of the Courant-Friedrichs-Lewy (CFL) condition, which is manifestly non-local. Although the v… ▽ More An optimally efficient explicit numerical scheme for solving fluid dynamics equations, or any other parabolic or hyperbolic system of partial differential equations, should allow local regions to advance in time with their own, locally constrained time steps. However, such a scheme can result in violation of the Courant-Friedrichs-Lewy (CFL) condition, which is manifestly non-local. Although the violations can be considered to be "weak" in a certain sense and the corresponding numerical solution may be stable, such calculation does not guarantee the correct propagation speed for arbitrary waves. We use an experimental fluid dynamics code that allows cubic "patches" of grid cells to step with independent, locally constrained time steps to demonstrate how the CFL condition can be enforced by imposing a condition on the time steps of neighboring patches. We perform several numerical tests that illustrate errors introduced in the numerical solutions by weak CFL condition violations and show how strict enforcement of the CFL condition eliminates these errors. In all our tests the strict enforcement of the CFL condition does not impose a significant performance penalty. △ Less

Submitted 9 January, 2018; originally announced January 2018.

Comments: Accepted for publication in the Journal of Computational Physics

Was there an early reionization component in our universe?

Authors: Pablo Villanueva-Domingo, Stefano Gariazzo, Nickolay Y. Gnedin, Olga Mena

Abstract: A deep understanding of the Epoch of Reionization is still missing in our knowledge of the universe. While future probes will allow us to test the precise evolution of the free electron fraction from redshifts between $z\simeq 6$ and $z\simeq 20$, at present one could ask what kind of reionization processes are allowed by present Cosmic Microwave Background temperature and polarization measurement… ▽ More A deep understanding of the Epoch of Reionization is still missing in our knowledge of the universe. While future probes will allow us to test the precise evolution of the free electron fraction from redshifts between $z\simeq 6$ and $z\simeq 20$, at present one could ask what kind of reionization processes are allowed by present Cosmic Microwave Background temperature and polarization measurements. An early contribution to reionization could imply a departure from the standard picture where star formation determines the reionization onset. BBy considering a broad class of possible reionization parameterizations, we find that current data do not require an early reionization component in our universe and that only one marginal class of models, based on a particular realization of reionization, may point to that. In addition, the frequentist Akaike Information Criterion (AIC) provides strong evidence against alternative reionization histories, favoring the most simple reionization scenario, which describes reionization by means of only one (constant) reionization optical depth $τ$. △ Less

Submitted 27 March, 2018; v1 submitted 7 December, 2017; originally announced December 2017.

Comments: 16 pages, 3 figures. Version accepted for publication in JCAP

Star Cluster Formation in Cosmological Simulations. II. Effects of Star Formation Efficiency and Stellar Feedback

Authors: Hui Li, Oleg Y. Gnedin, Nickolay Y. Gnedin

Abstract: The implementation of star formation and stellar feedback in cosmological simulations plays a critical role in shaping galaxy properties. In the first paper of the series, we presented a new method to model star formation as a collection of star clusters. In this paper, we improve the algorithm by eliminating accretion gaps, boosting momentum feedback, and introducing a subgrid initial bound fract… ▽ More The implementation of star formation and stellar feedback in cosmological simulations plays a critical role in shaping galaxy properties. In the first paper of the series, we presented a new method to model star formation as a collection of star clusters. In this paper, we improve the algorithm by eliminating accretion gaps, boosting momentum feedback, and introducing a subgrid initial bound fraction, $f_i$, that distinguishes cluster mass from stellar particle mass. We perform a suite of simulations with different star formation efficiency per freefall time $ε_{\rm ff}$ and supernova momentum feedback intensity $f_{\rm boost}$. We find that the star formation history of a Milky Way-sized galaxy is sensitive to $f_{\rm boost}$, which allows us to constrain its value, $f_{\rm boost}\approx5$, in the current simulation setup. Changing $ε_{\rm ff}$ from a few percent to 200\% has little effect on global galaxy properties. However, on smaller scales, the properties of star clusters are very sensitive to $ε_{\rm ff}$. We find that $f_i$ increases with $ε_{\rm ff}$ and cluster mass. Through the dependence on $f_i$, the shape of the cluster initial mass function varies strongly with $ε_{\rm ff}$. The fraction of clustered star formation and maximum cluster mass increase with the star formation rate surface density, with the normalization of both relations dependent on $ε_{\rm ff}$. The cluster formation timescale systematically decreases with increasing $ε_{\rm ff}$. Local variations in the gas accretion history lead to a 0.25~dex scatter for the integral cluster formation efficiency. Joint constraints from all the observables prefer the runs that produce a median integral efficiency of 16\%. △ Less

Submitted 25 June, 2018; v1 submitted 4 December, 2017; originally announced December 2017.

Comments: 18 pages, 14 figure., ApJ in press

Constraints on Spin of a Supermassive Black Hole in Quasars with Big Blue Bump

Authors: M. Yu. Piotrovich, Yu. N. Gnedin, T. M. Natsvlishvili, S. D. Buliga

Abstract: We determined the spin value of supermassive black hole (SMBH) in active galactic nuclei (AGN) with investigated ultraviolet-to-optical spectral energy distribution, presented in the sample of Shang et al. (2005). The estimates of the spin values have been produced at the base of the standard geometrically thin accretion disk model and with using the results of the polarimetric observations. The p… ▽ More We determined the spin value of supermassive black hole (SMBH) in active galactic nuclei (AGN) with investigated ultraviolet-to-optical spectral energy distribution, presented in the sample of Shang et al. (2005). The estimates of the spin values have been produced at the base of the standard geometrically thin accretion disk model and with using the results of the polarimetric observations. The polarimetric observations are very important for determining the inclination angle of AGN disk. We presented the results of our determinations of the radiation efficiency of the accretion flow and values of the spins of SMBHs, that derives the coefficient of radiation efficiency. The majority of SMBHs of AGNs from Shang et al. (2005) sample are to be the Kerr black holes with the high spin value. △ Less

Submitted 20 November, 2017; originally announced November 2017.

Comments: 5 pages, 1 figure

Warm Dark Matter and Cosmic Reionization

Authors: Pablo Villanueva-Domingo, Nickolay Y. Gnedin, Olga Mena

Abstract: In models with dark matter made of particles with keV masses, such as a sterile neutrino, small-scale density perturbations are suppressed, delaying the period at which the lowest mass galaxies are formed and therefore shifting the reionization processes to later epochs. In this study, focusing on Warm Dark Matter (WDM) with masses close to its present lower bound, i.e. around the $3$ keV region,… ▽ More In models with dark matter made of particles with keV masses, such as a sterile neutrino, small-scale density perturbations are suppressed, delaying the period at which the lowest mass galaxies are formed and therefore shifting the reionization processes to later epochs. In this study, focusing on Warm Dark Matter (WDM) with masses close to its present lower bound, i.e. around the $3$ keV region, we derive constraints from galaxy luminosy functions, the ionization history and the Gunn-Peterson effect. We show that even if star formation efficiency in the simulations is adjusted to match the observed UV galaxy luminosity functions in both CDM and WDM models, the full distribution of Gunn-Peterson optical depth retains the strong signature of delayed reionization in the WDM model. However, until the star formation and stellar feedback model used in modern galaxy formation simulations is constrained better, any conclusions on the nature of dark matter derived from reionization observables remain model-dependent. △ Less

Submitted 28 August, 2017; originally announced August 2017.

Comments: 7 pages, 6 figures

Baryonic effects in cosmic shear tomography: PCA parametrization and importance of extreme baryonic models

Authors: Irshad Mohammed, Nickolay Y. Gnedin

Abstract: Baryonic effects are amongst the most severe systematics to the tomographic analysis of weak lensing data which is the principal probe in many future generations of cosmological surveys like LSST, Euclid etc.. Modeling or parameterizing these effects is essential in order to extract valuable constraints on cosmological parameters. In a recent paper, Eifler et al. (2015) suggested a reduction techn… ▽ More Baryonic effects are amongst the most severe systematics to the tomographic analysis of weak lensing data which is the principal probe in many future generations of cosmological surveys like LSST, Euclid etc.. Modeling or parameterizing these effects is essential in order to extract valuable constraints on cosmological parameters. In a recent paper, Eifler et al. (2015) suggested a reduction technique for baryonic effects by conducting a principal component analysis (PCA) and removing the largest baryonic eigenmodes from the data. In this article, we conducted the investigation further and addressed two critical aspects. Firstly, we performed the analysis by separating the simulations into training and test sets, computing a minimal set of principle components from the training set and examining the fits on the test set. We found that using only four parameters, corresponding to the four largest eigenmodes of the training set, the test sets can be fitted thoroughly with an RMS $\sim 0.0011$. Secondly, we explored the significance of outliers, the most exotic/extreme baryonic scenarios, in this method. We found that excluding the outliers from the training set results in a relatively bad fit and degraded the RMS by nearly a factor of 3. Therefore, for a direct employment of this method to the tomographic analysis of the weak lensing data, the principle components should be derived from a training set that comprises adequately exotic but reasonable models such that the reality is included inside the parameter domain sampled by the training set. The baryonic effects can be parameterized as the coefficients of these principle components and should be marginalized over the cosmological parameter space. △ Less

Submitted 7 July, 2017; originally announced July 2017.

Comments: 7 pages, 8 figures, 1 table. Submitted to ApJ

Report number: FERMILAB-PUB-17-223-A

Constraining Recoiling Velocities of Black Holes Ejected by Gravitational Radiation in Galaxy Mergers

Authors: M. Yu. Piotrovich, Yu. N. Gnedin, T. M. Natsvlishvili, S. D. Buliga

Abstract: Recent general relativistic simulations have shown that the coalescence of two spinning black holes (BH) can lead to recoiling speeds of the BH remnant of up to thousands of km/s as a result of the gravitational radiation emission. It is important that the accretion disc remains bound to ejected BH within the region where the gas orbital velocity is larger than the ejection speed. We considered th… ▽ More Recent general relativistic simulations have shown that the coalescence of two spinning black holes (BH) can lead to recoiling speeds of the BH remnant of up to thousands of km/s as a result of the gravitational radiation emission. It is important that the accretion disc remains bound to ejected BH within the region where the gas orbital velocity is larger than the ejection speed. We considered the situation when the recoiling kick radius coincides with the radius of the broad line region (BLR). We show that in this situation the observed polarization data of accretion disk emission allow to determine the value of the recoil velocity. We present the estimates of the kick velocity for AGN with determined polarization data. △ Less

Submitted 5 May, 2017; originally announced May 2017.

Comments: 7 pages, 1 figures, 2 tables. Accepted for publication in Astrophysics and Space Sciences

The physical origin of long gas depletion times in galaxies

Authors: Vadim A. Semenov, Andrey V. Kravtsov, Nickolay Y. Gnedin

Abstract: We present a model that explains why galaxies form stars on a time scale significantly longer than the time scales of processes governing the evolution of interstellar gas. We show that gas evolves from a non-star-forming to a star-forming state on a relatively short time scale and thus the rate of this evolution does not limit the star formation rate. Instead, the star formation rate is limited b… ▽ More We present a model that explains why galaxies form stars on a time scale significantly longer than the time scales of processes governing the evolution of interstellar gas. We show that gas evolves from a non-star-forming to a star-forming state on a relatively short time scale and thus the rate of this evolution does not limit the star formation rate. Instead, the star formation rate is limited because only a small fraction of star-forming gas is converted into stars before star-forming regions are dispersed by feedback and dynamical processes. Thus, gas cycles into and out of star-forming state multiple times, which results in a long time scale on which galaxies convert gas into stars. Our model does not rely on the assumption of equilibrium and can be used to interpret trends of depletion times with the properties of observed galaxies and the parameters of star formation and feedback recipes in simulations. In particular, the model explains how feedback self-regulates the star formation rate in simulations and makes it insensitive to the local star formation efficiency. We illustrate our model using the results of an isolated $L_*$-sized galaxy simulation that reproduces the observed Kennicutt-Schmidt relation for both molecular and atomic gas. Interestingly, the relation for molecular gas is almost linear on kiloparsec scales, although a nonlinear relation is adopted in simulation cells. We discuss how a linear relation emerges from non-self-similar scaling of the gas density PDF with the average gas surface density. △ Less

Submitted 8 August, 2017; v1 submitted 13 April, 2017; originally announced April 2017.

Comments: 16 pages, 9 figures; accepted for publication in ApJ

Journal ref: ApJ 845, 133 (2017)