Engineering topics

Gnedin, Nickolay Y.

Publications and source records attributed to Gnedin, Nickolay Y..

8 records

DOE OSTI2024

On the Physical Nature of Lyα Transmission Spikes in High-redshift Quasar Spectra

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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Modeling Molecular Hydrogen in Low-metallicity Galaxies

We use a suite of hydrodynamics simulations of the interstellar medium (ISM) within a galactic disk, which includes radiative transfer, a nonequilibrium 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 H I–H 2 transition shifting to higher densities with decreasing metallicity, the maximum achieved molecular fraction in the ISM drops drastically at Z ≲ 0.2 Z ⊙ 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 ƒ 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 ≲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 ƒ 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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Do Minihalos Affect Cosmic Reionization?

The role of unresolved structures (“minihalos”) 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 minihalos are unimportant compared to recombinations produced in large-scale cosmic structures and inside more massive, fully resolved halos. Since both the 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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Cosmic Reionization on Computers: Statistics, Physical Properties, and Environments of Lyman Limit Systems at z ~ 6

Lyman limit systems (LLSs) are dense hydrogen clouds with high enough H i 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 ~ 6. We generate 3000 synthetic quasar sight lines 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 of LLSs spatially associated with nearby galaxies increases with H i 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 association of LLSs to large-scale structure informs our understanding of the IGM–galaxy connection during the "Epoch of Reionization," and provides a theoretical basis for interpreting future observations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Emergence of the temperature–density relation in the low-density intergalactic medium

We examine the evolution of the phase diagram of the low-density intergalactic medium during the Epoch of Reionization in simulation boxes with varying reionization histories from the Cosmic Reionization on Computers project. The probability density function (PDF) of gas temperature at fixed density exhibits two clear modes: a warm and a 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 parametrized 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 and the evolving gas density PDF. 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 ≲ X HI ≲ 10 –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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Exploring the dependence of gas cooling and heating functions on the incident radiation field with machine learning

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 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 behaviour of cooling and heating functions.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1995

Massive black holes and light-element nucleosynthesis in a baryonic universe

We reexamine the model proposed by Gnedin & Ostriker (1992) in which Jeans mass black holes (M(sub BH) approximately = 10(exp 6) solar mass) form shortly after decoupling. There is no nonbaryonic dark matter in this model, but we examine the possibility that Omega(sub b) is considerably larger than given by normal nucleosynthesis. Here we allow for the fact that much of the high baryon-to-photon ratio material will collapse leaving the universe of remaining material with light-element abundances more in accord with the residual baryonic density (approximately = 10(exp -2)) than with Omega(sub 0) and the initial baryonic density (approximately = 10(exp -1)). We find that no reasonable model can be made with random-phase density fluctuations, if the power on scales smaller than 10(exp 6) solar mass is as large as expected. However, phase-correlated models of the type that might occur in connection with topological singularities can be made with Omega(sub b) h(exp 2) = 0.013 +/- 0.001, 0.15 approximately less than Omega(sub 0) approximately less than 0.4, which are either flat (Omega(sub lambda) = 1 - Omega(sub 0)) or open (Omega(sub lambda) = 0) and which satisfy all the observational constraints which we apply, including the large baryon-to-total mass ratio found in the X-ray clusters. The remnant baryon density is thus close to that obtained in the standard picture (Omega(sub b) h(exp 2) = 0.0125 +/- 0.0025; Walker et al. 1991). The spectral index implied for fluctuations in the baryonic isocurvature scenario, -1 less than m less than 0, is in the range expected by other arguments based on large-scale structure and microwave fluctuation constraints. The dark matter in this picture is in the form of massive black holes. Accretion onto them at early epochs releases high-energy photons which significantly heat and reionize the universe. But photodissociation does not materially change light-element abundances. A typical model gives bar-y approximately = 1 x 10(exp -5), n(sub e)/n(sub H)(z = 30) approximately = 0.1, and a diffuse gamma-ray background at 100 keV near the Cosmic Background Explorer Satellite (COBE) limit of the order of 10% of that observed which originates from high-redshift quasars. Reionization in this model occurs at redshift 600 and reaches (H II/H(sub tot) approximately = 0.1-0.2.

Gnedin, Nickolay Y.↗

NASA NTRS1992

A tilted cold dark matter cosmological scenario

A new cosmological scenario based on CDM but with a power spectrum index of about 0.7-0.8 is suggested. This model is predicted by various inflationary models with no fine tuning. This tilted CDM model, if normalized to COBE, alleviates many problems of the standard CDM model related to both small-scale and large-scale power. A physical bias of galaxies over dark matter of about two is required to fit spatial observations.

Cen, Renyue↗