BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Nickolay Y. Gnedin

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

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-phastro-ph.IMhep-phastro-ph.HEphysics.comp-ph

Catalog footprint

What is connected

65works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

65 published item(s)

preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Testing multiflavored ULDM models with SPARC

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Approximating Density Probability Distribution Functions Across Cosmologies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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}$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Natively Periodic Fast Multipole Method: Approximating the Optimal Green Function

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reionization with galaxies and active galactic nuclei

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$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cosmic Reionization On Computers: Numerical and Physical Convergence

In this paper I show that simulations of reionization performed under the Cosmic Reionization On Computers (CROC) project do converge in space and mass, albeit rather slowly. A fully converged solution (for a given star formation and feedback model) can be determined at a level of precision of about 20%, but such a solution is useless in practice, since achieving it in production-grade simulations would require a large set of runs at various mass and spatial resolutions, and computational resources for such an undertaking are not yet readily available. In order to make progress in the interim, I introduce a weak convergence correction factor in the star formation recipe, which allows one to approximate the fully converged solution with finite resolution simulations. The accuracy of weakly converged simulations approaches a comparable, ~20% level of precision for star formation histories of individual galactic halos and other galactic properties that are directly related to star formation rates, like stellar masses and metallicities. Yet other properties of model galaxies, for example, their HI masses, are recovered in the weakly converged runs only within a factor of two.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cosmic Reionization On Computers: The Faint End of the Galaxy Luminosity Function

Using numerical cosmological simulations completed under the "Cosmic Reionization On Computers" (CROC) project, I explore theoretical predictions for the faint end of the galaxy UV luminosity functions at $z\geq 6$. A commonly used Schechter function approximation with the magnitude cut at $M_{\rm CUT}\sim-13$ provides a reasonable fit to the actual luminosity function of simulated galaxies. When the Schechter functional form is forced on the luminosity functions from the simulations, the magnitude cut $M_{\rm CUT}$ is found to vary between $-12$ and $-14$ with a mild redshift dependence. An analytical model of reionization from Madau, Haardt & Rees (1997), as used by Robertson et al. (2015), provides a good description of the simulated results, but only if the redshift dependence of the effective escape fraction (induced by physical processes not captured by the Madau, Haardt & Rees model) is accounted for.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cosmic Reionization On Computers. Mean and Fluctuating Redshifted 21 cm Signal

We explore the mean and fluctuating redshifted 21 cm signal in numerical simulations of cosmic reionization from the Cosmic Reionization On Computers (CROC) project. We find that the mean signal varies between about $\pm20\rm{mK}$. Most significantly, we find that the negative pre-reionization dip at $z\sim10-15$ only extends to $\langleΔT_B\rangle\sim-20\rm{mK}$, in agreement with prior simulation results and in significant contrast to Pritchard & Loeb analytical model, requiring substantially higher sensitivity from global signal experiments that operate in this redshift range (EDGES-II, LEDA, SCI-HI, and DARE). We also explore the role of dense substructure (filaments and embedded galaxies) in the formation of 21 cm power spectrum. We find that by neglecting the semi-neutral substructure inside ionized bubbles, the power spectrum can be mis-estimated by 25-50\% at scales $k\sim 0.1-1h\rm{Mpc}^{-1}$. This scale range is of a particular interest, because the upcoming 21 cm experiments (MWA, PAPER, HERA) are expected to be most sensitive within it.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Cosmic Reionization On Computers. Ultraviolet Continuum Slopes and Dust Opacities in High Redshift Galaxies

We compare the properties of stellar populations of model galaxies from the Cosmic Reionization On Computers (CROC) project with the exiting UV and IR data. Since CROC simulations do not follow cosmic dust directly, we adopt two variants of the dust-follows-metals ansatz to populate model galaxies with dust. Using the dust radiative transfer code Hyperion, we compute synthetic stellar spectra, UV continuum slopes, and IR fluxes for simulated galaxies. We find that the simulation results generally match observational measurements, but, perhaps, not in full detail. The differences seem to indicate that our adopted dust-follows-metals ansatzes are not fully sufficient. While the discrepancies with the exiting data are marginal, the future JWST data will be of much higher precision, rendering highly significant any tentative difference between theory and observations. It is, therefore, likely, that in order to fully utilize the precision of JWST observations, fully dynamical modeling of dust formation, evolution, and destruction may be required.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Line Overlap and Self-Shielding of Molecular Hydrogen in Galaxies

The effect of line overlap in the Lyman and Werner bands, often ignored in galactic studies of the atomic-to-molecular transition, greatly enhances molecular hydrogen self-shielding in low metallicity environments, and dominates over dust shielding for metallicities below about 10% solar. We implement that effect in cosmological hydrodynamics simulations with an empirical model, calibrated against the observational data, and provide fitting formulae for the molecular hydrogen fraction as a function of gas density on various spatial scales and in environments with varied dust abundance and interstellar radiation field. We find that line overlap, while important for detailed radiative transfer in the Lyman and Werner bands, has only a minor effect on star formation on galactic scales, which, to a much larger degree, is regulated by stellar feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the Proper Use of the Reduced Speed of Light Approximation

I show that the Reduced Speed of Light (RSL) approximation, when used properly (i.e. as originally designed - only for the local sources but not for the cosmic background), remains a highly accurate numerical method for modeling cosmic reionization. Simulated ionization and star formation histories from the "Cosmic Reionization On Computers" (CROC) project are insensitive to the adopted value of the reduced speed of light for as long as that value does not fall below about 10% of the true speed of light. A recent claim of the failure of the RSL approximation in the Illustris reionization model appears to be due to the effective speed of light being reduced in the equation for the cosmic background too, and, hence, illustrates the importance of maintaining the correct speed of light in modeling the cosmic background.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cosmic Reionization On Computers III. The Clumping Factor

We use fully self-consistent numerical simulations of cosmic reionization, completed under the Cosmic Reionization On Computers (CROC) project, to explore how well the recombinations in the ionized IGM can be quantified by the effective "clumping factor." The density distribution in the simulations (and, presumably, in a real universe) is highly inhomogeneous and more-or-less smoothly varying in space. However, even in highly complex and dynamic environments the concept of the IGM remains reasonably well-defined; the largest ambiguity comes from the unvirialized regions around galaxies that are over-ionized by the local enhancement in the radiation field ("proximity zones"). That ambiguity precludes computing the IGM clumping factor to better than about 20\%. We also discuss a "local clumping factor," defined over a particular spatial scale, and quantify its scatter on a given scale and its variation as a function of scale.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Effects of Dark Matter Annihilation on Cosmic Reionization

We revisit the possibility of constraining the properties of dark matter (DM) by studying the epoch of cosmic reionization. Previous studies have shown that DM annihilation was unlikely to have provided a large fraction of the photons that ionized the universe, but instead played a subdominant role relative to stars and quasars. The DM, however, begins to efficiently annihilate with the formation of primordial microhalos at $z\sim100-200$, much earlier than the formation of the first stars. Therefore, if DM annihilation ionized the universe at even the percent level over the interval $z \sim 20-100$, it can leave a significant imprint on the global optical depth, $τ$. Moreover, we show that cosmic microwave background (CMB) polarization data and future 21 cm measurements will enable us to more directly probe the DM contribution to the optical depth. In order to compute the annihilation rate throughout the epoch of reionization, we adopt the latest results from structure formation studies and explore the impact of various free parameters on our results. We show that future measurements could make it possible to place constraints on the dark matter's annihilation cross section that are at a level comparable to those obtained from the observations of dwarf galaxies, cosmic ray measurements, and studies of recombination.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cosmic Reionization On Computers I. Design and Calibration of Simulations

Cosmic Reionization On Computers (CROC) is a long-term program of numerical simulations of cosmic reionization. Its goal is to model fully self-consistently (albeit not necessarily from the first principles) all relevant physics, from radiative transfer to gas dynamics and star formation, in simulation volumes of up to 100 comoving Mpc, and with spatial resolution approaching 100 pc in physical units. In this method paper we describe our numerical method, the design of simulations, and the calibration of numerical parameters. Using several sets (ensembles) of simulations in 20 Mpc/h and 40 Mpc/h boxes with spatial resolution reaching 125 pc at z=6, we are able to match the observed galaxy UV luminosity functions at all redshifts between 6 and 10, as well as obtain reasonable agreement with the observational measurements of the Gunn-Peterson optical depth at z<6.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cosmic Reionization On Computers II. Reionization History and Its Back Reaction on Early Galaxies

We compare the results from several sets of cosmological simulations of cosmic reionization, produced under Cosmic Reionization On Computers (CROC) project, with existing observational data on the high-redshift Ly-alpha forest and the abundance of Ly-alpha emitters. We find good consistency with the observational measurements and the previous simulation work. By virtue of having several independent realizations for each set of numerical parameters, we are able to explore the effect of cosmic variance on observable quantities. One unexpected conclusion we are forced into is that cosmic variance is unusually large at z>6, with both our simulations and, most likely, observational measurements are still not fully converged for even such basic quantities as the average Gunn-Peterson optical depth or the volume-weighted neutral fraction. We also find that reionization has little effect on the early galaxies or on global cosmic star formation history, because galaxies whose gas content is affected by photoionization contain no molecular (i.e. star-forming) gas in the first place. In particular, measurements of the faint end of the galaxy luminosity function by JWST are unlikely to provide a useful constraint on reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Emergence of the Kennicutt-Schmidt Relation from the Small-Scale SFR-Density Relation

We use simulations of isolated galaxies with a few parsec resolution to explore the connection between the small-scale star formation rate - gas density relation and the induced large-scale correlation between the star formation rate surface density and the surface density of the molecular gas (the Kennicutt-Schmidt relation). We find that, in the simulations, a power-law small-scale "star formation law" directly translates into an identical power-law Kennicutt-Schmidt relation. If this conclusion holds in the reality as well, it implies that the observed approximately linear Kennicutt-Schmidt relation must reflect the approximately linear small-scale "star formation law".

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Modeling Physical Processes at Galactic Scales and Above

What should these lectures be? The subject assigned to us is so broad that many books can be written about it. So, in planning these lectures I had several options. One would be to focus on a narrow subset of topics and to cover them in great detail. Such a subset necessarily would be highly personal and useful to a few read- ers at best. Another option would be to give a very shallow overview of the whole field, but then it won't be very much different from a highly compressed version of a university course (which anyone can take if they wish so). So, I decided to be selfish and to prepare these lectures as if I was teaching my own graduate student. Given my research interests, I selected what the student would need to know to be able to discuss science with me and to work on joint research projects. So, the story presented below is both personal and incomplete, but it does cover several subjects that are poorly represented in the existing textbooks (if at all). Some of topics I focus on below are closely connected, others are disjoint, some are just side detours on specific technical questions. There is an overlapping theme, however. Our goal is to follow the cosmic gas from large scales, low densities, (rel- atively) simple physics to progressively smaller scales, higher densities, closer rela- tion to galaxies, and more complex and uncertain physics. So, we (you - the reader, and me - the author) are going to follow a "yellow brick road" from the gas well be- yond any galaxy confines to the actual sites of star formation and stellar feedback. On the way we will stop at some places for a tour and run without looking back through some others. So, the road will be uneven, but I hope that some readers find it useful.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the mass of the Local Group

We use recent proper motion measurements of the tangential velocity of M31, along with its radial velocity and distance, to derive the likelihood of the sum of halo masses of the Milky Way and M31. This is done using a sample halo pairs in the Bolshoi cosmological simulation of $Λ$CDM cosmology selected to match properties and environment of the Local Group. The resulting likelihood gives estimate of the sum of masses of $M_{\rm MW,200}+M_{\rm M31,200}=$ $2.40_{-1.05}^{+1.95}\times10^{12}\,M_{\odot}$ ($90\%$ confidence interval). This estimate is consistent with individual mass estimates for the Milky Way and M31 and is consistent, albeit somewhat on the low side, with the mass estimated using the timing argument. We show that although the timing argument is unbiased on average for all pairs, for pairs constrained to have radial and tangential velocities similar to that of the Local Group the argument overestimates the sum of masses by a factor of $1.6$. Using similar technique we estimate the total dark matter mass enclosed within $1$ Mpc from the Local Group barycenter to be $M_{\rm LG}(r<1\, {\rm Mpc})=4.2_{-2.0}^{+3.4}\times10^{12}\,M_{\odot}$ ($90\%$ confidence interval).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Recombination clumping factor during cosmic reionization

We discuss the role of recombinations in the IGM, and the related concept of the clumping factor, during cosmic reionization. The clumping factor is, in general, a local quantity that depends on the local over-density and the scale below which the baryon density field can be assumed smooth. That scale, called the filtering scale, is itself depended on over-density and local thermal history. We present a method for building a self-consistent analytical model of inhomogeneous reionization assuming the linear growth rate of the density fluctuation, which accounts for these effects simultaneously. We show that taking into account the local clumping factor introduces significant corrections to the total recombination rate, comparing to the model with a globally uniform clumping factor.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Satellites in MW-like hosts: Environment dependence and close pairs

Previous studies showed that an estimate of the likelihood distribution of the Milky Way halo mass can be derived using the properties of the satellites similar to the Large and Small Magellanic Clouds (LMC and SMC). However, it would be straightforward to interpret such an estimate only if the properties of the Magellanic Clouds (MCs) are fairly typical and are not biased by the environment. In this study we explore whether the environment of the Milky Way affects the properties of the SMC and LMC such as their velocities. To test for the effect of the environment, we compare velocity distributions for MC-sized subhalos around Milky Way hosts in a sample selected simply by mass and in the second sample of such halos selected with additional restrictions on the distance to the nearest cluster and the local galaxy density, designed to mimic the environment of the Local Group (LG). We find that satellites in halos in the LG-like environments do have somewhat larger velocities, as compared to the halos of similar mass in the sample without environmental constraints. We derive the host halo likelihood distribution for the samples in the LG-like envirionment and in the control sample and find that the environment does not significantly affect the derived likelihood. We use the updated properties of the SMC and LMC to derive the constraint on the MW halo mass $\log{({\rm M}_{200} /\msol)}=12.06^{+0.31}_{-0.19}$ (90% confidence interval). We also explore the incidence of close pairs with relative velocities and separations similar to those of the LMC and SMC and find that such pairs are quite rare among $Λ$CDM halos. Taking into account the close separation of the MCs in the Busha et al.\ 2011 method results in the shift of the MW halo mass estimate to smaller masses, with the peak shifting approximately by a factor of two.[Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Simulations of disk galaxies with cosmic ray driven galactic winds

We present results from high-resolution hydrodynamic simulations of isolated SMC- and Milky Way-sized galaxies that include a model for feedback from galactic cosmic rays (CRs). We find that CRs are naturally able to drive winds with mass loading factors of up to ~10 in dwarf systems. The scaling of the mass loading factor with circular velocity between the two simulated systems is consistent with \propto v_c^{1-2} required to reproduce the faint end of the galaxy luminosity function. In addition, simulations with CR feedback reproduce both the normalization and the slope of the observed trend of wind velocity with galaxy circular velocity. We find that winds in simulations with CR feedback exhibit qualitatively different properties compared to SN driven winds, where most of the acceleration happens violently in situ near star forming sites. In contrast, the CR-driven winds are accelerated gently by the large-scale pressure gradient established by CRs diffusing from the star-forming galaxy disk out into the halo. The CR-driven winds also exhibit much cooler temperatures and, in the SMC-sized system, warm (T~10^4 K) gas dominates the outflow. The prevalence of warm gas in such outflows may provide a clue as to the origin of ubiquitous warm gas in the gaseous halos of galaxies detected via absorption lines in quasar spectra.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Circum-Galactic Gas and the Isotropic Gamma Ray Background

Interactions of cosmic rays with the interstellar gas and radiation fields of the Milky Way provide the majority of the gamma rays observed by the Fermi Gamma Ray Space Telescope. In addition to the gas which is densely concentrated along the Galactic Disk, hydrodynamical simulations and observational evidence favor the presence of a halo of hot (T~10^6 K) ionized hydrogen (H_II), extending with non-negligible densities out to the virial radius of the Milky Way. We show that cosmic ray collisions with this circum-galactic gas should be expected to provide a significant flux of gamma rays, on the order of 10% of the observed isotopic gamma ray background at energies above 1 GeV. In addition, gamma rays originating from the extended H_II halos of other galaxies along a given line-of-sight should contribute to this background at a similar level.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cooling and Heating Functions of Photoionized Gas

Cooling and heating functions of cosmic gas are a crucial ingredient for any study of gas dynamics and thermodynamics in the interstellar and intergalactic medium. As such, they have been studied extensively in the past under the assumption of collisional ionization equilibrium. However, for a wide range of applications, the local radiation field introduces a non-negligible, often dominant, modification to the cooling and heating functions. In the most general case, these modifications cannot be described in simple terms, and would require a detailed calculation with a large set of chemical species using a radiative transfer code (the well-known code Cloudy, for example). We show, however, that for a sufficiently general variation in the spectral shape and intensity of the incident radiation field, the cooling and heating functions can be approximated as depending only on several photoionization rates, which can be thought of as representative samples of the overall radiation field. This dependence is easy to tabulate and implement in cosmological or galactic-scale simulations, thus economically accounting for an important but rarely-included factor in the evolution of cosmic gas. We also show a few examples where the radiation environment has a large effect, the most spectacular of which is a quasar that suppresses gas cooling in its host halo without any mechanical or non-radiative thermal feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effect of halo bias and Lyman Limit Systems on the history of cosmic reionization

We extend the existing analytical model of reionization by Furlanetto et al. (2004) to include the biasing of reionization sources and additional absorption by Lyman Limit systems. Our model is, by construction, consistent with the observed evolution of the galaxy luminosity function at z<8 and with the observed evolution of Ly-α forest at z<6. We also find that, for a wide range of values for the relative escape fraction that we consider reasonable, and which are consistent with the observational constraints on the relative escape fraction from lower redshifts, our reionization model is consistent with the WMAP constraint on the Thompson optical depth and with the SPT and EDGES constraints on the duration of reionization. We, therefore, conclude that it is possible to develop physically realistic models of reionization that are consistent with all existing observational constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the Baryonic Contents of Low Mass Galaxies

The baryonic Tully-Fisher relation is an important observational constraint on cosmological and galactic models. However, it is critical to keep in mind that in observations only stars, molecular, and atomic gas are counted, while the contribution of the ionized gas is almost universally missed. The ionized gas is, however, expected to be present in the gaseous disks of dwarf galaxies simply because they are exposed to the cosmic ionizing background and to the stellar radiation that manages to escape from the central regions of the galactic disks into their outer layers. Such an expectation is, indeed, born out both by cosmological numerical simulations and by simple analytical models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the Origin of the High Column Density Turnover in the HI Column Density Distribution

We study the high column density regime of the HI column density distribution function and argue that there are two distinct features: a turnover at NHI ~ 10^21 cm^-2 which is present at both z=0 and z ~ 3, and a lack of systems above NHI ~ 10^22 cm^-2 at z=0. Using observations of the column density distribution, we argue that the HI-H2 transition does not cause the turnover at NHI ~ 10^21 cm^-2, but can plausibly explain the turnover at NHI > 10^22 cm^-2. We compute the HI column density distribution of individual galaxies in the THINGS sample and show that the turnover column density depends only weakly on metallicity. Furthermore, we show that the column density distribution of galaxies, corrected for inclination, is insensitive to the resolution of the HI map or to averaging in radial shells. Our results indicate that the similarity of HI column density distributions at z=3 and z=0 is due to the similarity of the maximum HI surface densities of high-z and low-z disks, set presumably by universal processes that shape properties of the gaseous disks of galaxies. Using fully cosmological simulations, we explore other candidate physical mechanisms that could produce a turnover in the column density distribution. We show that while turbulence within GMCs cannot affect the DLA column density distribution, stellar feedback can affect it significantly if the feedback is sufficiently effective in removing gas from the central 2-3 kpc of high-redshift galaxies. Finally, we argue that it is meaningful to compare column densities averaged over ~ kpc scales with those estimated from quasar spectra which probe sub-pc scales due to the steep power spectrum of HI column density fluctuations observed in nearby galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Reionization history and CMB parameter estimation

We study how uncertainty in the reionization history of the universe affects estimates of other cosmological parameters from the Cosmic Microwave Background. We analyze WMAP7 data and synthetic Planck-quality data generated using a realistic scenario for the reionization history of the universe obtained from high-resolution numerical simulation. We perform parameter estimation using a simple sudden reionization approximation, and using the Principal Component Analysis (PCA) technique proposed by Mortonson and Hu. We reach two main conclusions: (1) Adopting a simple sudden reionization model does not introduce measurable bias into values for other parameters, indicating that detailed modeling of reionization is not necessary for the purpose of parameter estimation from future CMB data sets such as Planck. (2) PCA analysis does not allow accurate reconstruction of the actual reionization history of the universe in a realistic case.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simulations of Early Baryonic Structure Formation with Stream Velocity: I. Halo Abundance

It has been recently shown that the relative velocity between the dark matter and the baryons (vbc) at the time of recombination can affect the structure formation in the early universe (Tseliakhovich & Hirata 2010). We statistically quantify this effect using large cosmological simulations. We use three different high resolution sets of simulations (with separate transfer functions for baryons and dark matter) that vary in box size, particle number, and the value of the relative velocity between dark matter and baryons. We show that the total number density of halos is suppressed by ~ 20% at z = 25 for vbc = 1σ(vbc), where σ(vbc) is the variance of the relative velocity, while for vbc = 3.4σ(vbc) the relative suppression at the same redshift reaches 50%, remaining at or above the 30% level all the way to z = 11. We also find high abundance of "empty halos", i.e., halos that have gas fraction below half of the cosmic mean baryonic fraction fb. Specifically we find that for vbc = 1σ(vbc) all halos below 10^5M\odot are empty at z \geq 19. The high abundance of empty halos results in significant delay in the formation of gas rich mini-halos and the first galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simulations of Early Baryonic Structure Formation with Stream Velocity: II. The Gas Fraction

Understanding the gas content of high redshift halos is crucial for studying the formation of the first generation of galaxies and reionization. Recently, Tseliakhovich & Hirata showed that the relative "stream" velocity between the dark matter and baryons at the time of recombination - formally a second order effect, but an unusually large one - can influence the later structure formation history of the Universe. We quantify the effect of the stream velocity on the so-called "characteristic mass" - the minimum mass of a dark matter halo capable of retaining most of its baryons throughout its formation epoch - using three different high-resolution sets of cosmological simulations (with separate transfer functions for baryons and dark matter) that vary in box size, particle number, and the value of the relative velocity between the dark matter and baryons. In order to understand this effect theoretically, we generalize the linear theory filtering mass to properly account for the difference between the dark matter and baryonic density fluctuation evolution induced by the stream velocity. We show that the new filtering mass provides an accurate estimate for the characteristic mass, while other theoretical ansatzes for the characteristic mass are substantially less precise.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The relation between mid-plane pressure and molecular hydrogen in galaxies: Environmental dependence

Molecular hydrogen (H2) is the primary component of the reservoirs of cold, dense gas that fuel star formation in our galaxy. While the H2 abundance is ultimately regulated by physical processes operating on small scales in the interstellar medium (ISM), observations have revealed a tight correlation between the ratio of molecular to atomic hydrogen in nearby spiral galaxies and the pressure in the mid-plane of their disks. This empirical relation has been used to predict H2 abundances in galaxies with potentially very different ISM conditions, such as metal-deficient galaxies at high redshifts. Here, we test the validity of this approach by studying the dependence of the pressure -- H2 relation on environmental parameters of the ISM. To this end, we follow the formation and destruction of H2 explicitly in a suite of hydrodynamical simulations of galaxies with different ISM parameters. We find that a pressure -- H2 relation arises naturally in our simulations for a variety of dust-to-gas ratios or strengths of the interstellar radiation field in the ISM. Fixing the dust-to-gas ratio and the UV radiation field to values measured in the solar neighborhood results in fair agreement with the relation observed in nearby galaxies with roughly solar metallicity. However, the parameters (slope and normalization) of the pressure -- H2 relation vary in a systematical way with ISM properties. A particularly strong trend is the decrease of the normalization of the relation with a lowering of the dust-to-gas ratio of the ISM. We show that this trend and other properties of the pressure -- H2 relation are natural consequences of the transition from atomic to molecular hydrogen with gas surface density.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The X-factor in Galaxies: II. The molecular hydrogen -- star formation relation

There is ample observational evidence that the star formation rate (SFR) surface density, Sigma_SFR, is closely correlated with the surface density of molecular hydrogen, Sigma_H2. This empirical relation holds both for galaxy-wide averages and for individual >=kpc sized patches of the interstellar medium (ISM), but appears to degrade substantially at a sub-kpc scale. Identifying the physical mechanisms that determine the scale-dependent properties of the observed Sigma_H2-Sigma_SFR relation remains a challenge from a theoretical perspective. To address this question, we analyze the slope and scatter of the Sigma_H2-Sigma_SFR relation using a set of cosmological, galaxy formation simulations with a peak resolution of ~100 pc. These simulations include a chemical network for molecular hydrogen, a model for the CO emission, and a simple, stochastic prescription for star formation that operates on ~100 pc scales. Specifically, star formation is modeled as a Poisson process in which the average SFR is directly proportional to the present mass of H2. The predictions of our numerical model are in good agreement with the observed Kennicutt-Schmidt and Sigma_H2-Sigma_SFR relations. We show that observations based on CO emission are ill suited to reliably measure the slope of the latter relation at low (<20 M_sun pc^-2) H2 surface densities on sub-kpc scales. Our models also predict that the inferred Sigma_H2-Sigma_SFR relation steepens at high H2 surface densities as a result of the surface density dependence of the CO/H2 conversion factor. Finally, we show that on sub-kpc scales most of the scatter in the relation is a consequence of discreteness effects in the star formation process. In contrast, variations of the CO/H2 conversion factor are responsible for most of the scatter measured on super-kpc scales.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Towards a complete accounting of energy and momentum from stellar feedback in galaxy formation simulations

Stellar feedback plays a key role in galaxy formation by regulating star formation, driving interstellar turbulence and generating galactic scale outflows. Although modern simulations of galaxy formation can resolve scales of 10-100 pc, star formation and feedback operate on smaller, "subgrid" scales. Great care should therefore be taken in order to properly account for the effect of feedback on global galaxy evolution. We investigate the momentum and energy budget of feedback during different stages of stellar evolution, and study its impact on the interstellar medium using simulations of local star forming regions and galactic disks at the resolution affordable in modern cosmological zoom-in simulations. In particular, we present a novel subgrid model for the momentum injection due to radiation pressure and stellar winds from massive stars during early, pre-supernova evolutionary stages of young star clusters. Early injection of momentum acts to clear out dense gas in star forming regions, hence limiting star formation. The reduced gas density mitigates radiative losses of thermal feedback energy from subsequent supernova explosions, leading to an increased overall efficiency of stellar feedback. The detailed impact of stellar feedback depends sensitively on the implementation and choice of parameters. Somewhat encouragingly, we find that implementations in which feedback is efficient lead to approximate self-regulation of global star formation efficiency. We compare simulation results using our feedback implementation to other phenomenological feedback methods, where thermal feedback energy is allowed to dissipate over time scales longer than the formal gas cooling time. We find that simulations with maximal momentum injection suppress star formation to a similar degree as is found in simulations adopting adiabatic thermal feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Comparison of Methods for Determining the Molecular Content of Model Galaxies

Recent observations indicate that star formation occurs only in the molecular phase of a galaxy's interstellar medium. A realistic treatment of star formation in simulations and analytic models of galaxies therefore requires that one determine where the transition from the atomic to molecular gas occurs. In this paper we compare two methods for making this determination in cosmological simulations where the internal structures of molecular clouds are unresolved: a complex time-dependent chemistry network coupled to a radiative transfer calculation of the dissociating ultraviolet (UV) radiation field, and a simple time-independent analytic approximation. We show that these two methods produce excellent agreement at all metallicities >~10^-2 of the Milky Way value across a very wide range of UV fields. At lower metallicities the agreement is worse, likely because time-dependent effects become important; however, there are no observational calibrations of molecular gas content at such low metallicities, so it is unclear if either method is accurate. The comparison suggests that, in many but not all applications, the analytic approximation provides a viable and nearly cost-free alternative to full time-dependent chemistry and radiative transfer.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Halo Contraction Effect in Hydrodynamic Simulations of Galaxy Formation

The condensation of gas and stars in the inner regions of dark matter halos leads to a more concentrated dark matter distribution. While this effect is based on simple gravitational physics, the question of its validity in hierarchical galaxy formation has led to an active debate in the literature. We use a collection of several state-of-the-art cosmological hydrodynamic simulations to study the halo contraction effect in systems ranging from dwarf galaxies to clusters of galaxies, at high and low redshift. The simulations are run by different groups with different codes and include hierarchical merging, gas cooling, star formation, and stellar feedback. We show that in all our cases the inner dark matter density increases relative to the matching simulation without baryon dissipation, at least by a factor of several. The strength of the contraction effect varies from system to system and cannot be reduced to a simple prescription. We present a revised analytical model that describes the contracted mass profile to an rms accuracy of about 10%. The model can be used to effectively bracket the response of the dark matter halo to baryon dissipation. The halo contraction effect is real and must be included in modeling of the mass distribution of galaxies and galaxy clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Implementing the DC Mode in Cosmological Simulations with Supercomoving Variables

As emphasized by previous studies, proper treatment of the density fluctuation on the fundamental scale of a cosmological simulation volume - the "DC mode" - is critical for accurate modeling of spatial correlations on scales > 10% of simulation box size. We provide further illustration of the effects of the DC mode on the abundance of halos in small boxes and show that it is straightforward to incorporate this mode in cosmological codes that use the "supercomoving" variables. The equations governing evolution of dark matter and baryons recast with these variables are particularly simple and include the expansion factor, and hence the effect of the DC mode, explicitly only in the Poisson equation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Non-Thermal Dark Matter Mimicking An Additional Neutrino Species In The Early Universe

The South Pole Telescope (SPT), Atacama Cosmology Telescope (ACT), and Wilkinson Microwave Anisotropy Probe (WMAP) have each reported measurements of the cosmic microwave background's (CMB) angular power spectrum which favor the existence of roughly one additional neutrino species, in addition to the three contained in the standard model of particle physics. Neutrinos influence the CMB by contributing to the radiation density, which alters the expansion rate of the universe during the epoch leading up to recombination. In this paper, we consider an alternative possibility that the excess kinetic energy implied by these measurements was possessed by dark matter particles that were produced through a non-thermal mechanism, such as late-time decays. In particular, we find that if a small fraction (<1%) of the dark matter in the universe today were produced through the decays of a heavy and relatively long-lived state, the expansion history of the universe can be indistinguishable from that predicted in the standard cosmological model with an additional neutrino. Furthermore, if these decays take place after the completion of big bang nucleosynthesis, this scenario can avoid tension with the value of three neutrino species preferred by measurements of the light element abundances.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On determining the shape of matter distributions

A basic property of objects, like galaxies and halos that form in cosmological structure formation simulations, is their shape. Here, we critically investigate shape determination methods that are commonly used in the literature. It is found that using an enclosed integration volume and weight factors r^{-2} and r_{ell}^{-2} (elliptical radius) for the contribution of each particle or volume element in the shape tensor leads to biased axis ratios and smoothing of details when calculating the local shape as a function of distance from the center. To determine the local shape of matter distributions as a function of distance for well resolved objects (typically more than O(10^4) particles), we advocate a method that (1) uses an ellipsoidal shell (homoeoid) as an integration volume without any weight factors in the shape tensor and (2) removes subhalos.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The impact of baryon physics on the structure of high-redshift galaxies

We study the detailed structure of galaxies at redshifts z > 2 using cosmological simulations with improved modeling of the interstellar medium and star formation. The simulations follow the formation and dissociation of molecular hydrogen, and include star formation only in cold molecular gas. The molecular gas is more concentrated towards the center of galaxies than the atomic gas, and as a consequence, the resulting stellar distribution is very compact. For halos with total mass above 10^{11} Mo, the median half-mass radius of the stellar disks is 0.8 kpc at z = 3. The vertical structure of the molecular disk is much thinner than that of the atomic neutral gas. Relative to the non-radiative run, the inner regions of the dark matter halo change shape from prolate to mildly oblate and align with the stellar disk. However, we do not find evidence for a significant dark disk of dark matter around the stellar disk. The outer halo regions retain the orientation acquired during accretion and mergers, and are significantly misaligned with the inner regions. The radial profile of the dark matter halo contracts in response to baryon dissipation, establishing an approximately isothermal profile throughout most of the halo. This effect can be accurately described by a modified model of halo contraction. The angular momentum of a fixed amount of inner dark matter is approximately conserved over time, while in the dissipationless case most of it is transferred outward during mergers. The conservation of the dark matter angular momentum provides supporting evidence for the validity of the halo contraction model in a hierarchical galaxy formation process.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The X-factor in Galaxies: I. Dependence on Environment and Scale

Characterizing the conversion factor between CO emission and column density of molecular hydrogen, X_CO, is crucial in studying the gaseous content of galaxies, its evolution, and relation to star formation. In most cases the conversion factor is assumed to be close to that of giant molecular clouds (GMCs) in the Milky Way, except possibly for mergers and star-bursting galaxies. However, there are physical grounds to expect that it should also depend on the gas metallicity, surface density, and strength of the interstellar radiation field. The XCO factor may also depend on the scale on which CO emission is averaged due to effects of limited resolution. We study the dependence of X_CO on gas properties and averaging scale using a model that is based on a combination of results of sub-pc scale magneto-hydrodynamic simulations and on the gas distribution from self-consistent cosmological simulations of galaxy formation. Our model predicts a value of X_CO that is consistent with the Galactic value for interstellar medium conditions typical for the Milky Way. For such conditions the predicted X_CO varies by only a factor of two for gas surfaced densities in the range \sim 50 - 500 M_sun / pc^2. However, the model also predicts that more generally on the scale of GMCs, X_CO is a strong function of metallicity, and depends on the column density and the interstellar UV flux. We show explicitly that neglecting these dependencies in observational estimates can strongly bias the inferred distribution of H2 column densities of molecular clouds to have a narrower and offset range compared to the true distribution. We find that when averaged on \sim kpc scales the X-factor depends only weakly on radiation field and column density, but is still a strong function of metallicity. The predicted metallicity dependence can be approximated as X_CO \sim Z^{-γ} with γ ~ 0.5 - 0.8.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Ultra-faint dwarf galaxies as a test of early enrichment and metallicity-dependent star formation

The tight relation of star formation with molecular gas indicated by observations and assumed in recent models implies that the efficiency with which galaxies convert their gas into stars depends on gas metallicity. This is because the abundance of molecular hydrogen is sensitive to the abundance of dust, which catalyzes the formation of H_2 and helps to shield it from dissociating radiation. In this study we point out that in the absence of significant pre-enrichment by Population III stars forming out of zero metallicity gas, such H_2-based star formation is expected to leave an imprint in the form of bi-modality in the metallicity distribution among dwarf galaxies and in the metallicity distribution of stars within individual galaxies. The bi-modality arises because when gas metallicity (and dust abundance) is low, formation of molecular gas is inefficient, the gas consumption time scale is long, and star formation and metal enrichment proceed slowly. When metallicity reaches a critical threshold value star formation and enrichment accelerate, which leads to rapid increase in both stellar mass and metallicity of galaxies. We demonstrate this process both using a simple analytical model and full cosmological simulations. In contrast, observed metallicity distributions of dwarf galaxies or stars within them are not bi-modal. We argue that this discrepancy points to substantial early stochastic pre-enrichment by population III stars to levels Z ~ 0.01 Z_sun in dense, star forming regions of early galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Analysis of methods for detecting the proximity effect in quasar spectra

Using numerical simulations of structure formation, we investigate multiple methods of determining the strength of the proximity effect in the HI Lyalpha forest. We analyze three high resolution (~10kpc) redshift snapshots (z=4,3,2.25) of a Hydro-Particle-Mesh simulation to obtain realistic absorption spectra of the HI Lyalpha forest. We begin our analysis investigating the intrinsic biases thought to arise in the widely adopted standard technique of combining multiple lines of sight when searching for the proximity effect. We confirm the existence of this biases. We then concentrate on the analysis of the proximity effect along individual lines of sight. We construct the proximity effect strength distribution (PESD) and confirm that the PESD inferred from a simple averaging technique accurately recovers the input strength of the proximity effect at all redshifts. Moreover, the PESD closely follows the behaviors found in observed samples of quasar spectra. However, the PESD obtained from our new simulated sight lines presents some differences to that of simple Monte Carlo simulations. After developing three new theoretical methods of recovering the strength of the proximity effect on individual lines of sight, we compare their accuracy to the PESD from the simple averaging technique. All our new approaches are based on the maximization of the likelihood function, albeit invoking some modifications. The new techniques presented here fail to recover the input proximity effect in an un-biased way. Thus, employing complex 3D simulations, we provide strong evidence in favor of the proximity effect strength distribution obtained from the simple averaging technique, as method of estimating the UV background intensity, free of any biases.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Effect of Cosmic UV Background on Star Formation in High Redshift Galaxies

The effect of the cosmic UV background on star formation in high redshift galaxies is explored by means of high resolutions cosmological simulations. The simulations include star formation, 3D radiative transfer, and a highly detailed ISM model, and reach spatial resolution sufficient to resolve formation sites for molecular clouds. In the simulations the local radiation field in the Lyman-Werner band around star-forming molecular clouds dominates over the cosmic UV background by a factor of 100, similarly to the interstellar radiation field in the Milky Way and in a few high redshift galaxies for which measurements exist. The cosmic UV background, therefore, is essentially irrelevant for star formation in normal galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Environmental Dependence of the Kennicutt-Schmidt Relation in Galaxies

We present a detailed description of a phenomenological H2 formation model and local star formation prescription based on the density of molecular (rather than total) gas. Such approach allows us to avoid the arbitrary density and temperature thresholds typically used in star formation recipes. We present results of the model based on realistic cosmological simulations of high-z galaxy formation for a grid of numerical models with varied dust-to-gas ratios and interstellar far UV (FUV) fluxes. Our results show that both the atomic-to-molecular transition on small, ~10 pc scales and the Kennicutt-Schmidt (KS) relation on ~kpc scales are sensititive to the dust-to-gas ratio and the FUV flux. The atomic-to-molecular transition as a function of gas density or column density has a large scatter but is rather sharp and shifts to higher densities with decreasing dust-to-gas ratio and/or increasing FUV flux. Consequently, star formation is concentrated to higher gas surface density regions, resulting in steeper slope and lower amplitude of the KS relation at a given gas surface density, in less dusty and/or higher FUV flux environments. These trends should have a particularly strong effect on the evolution of low-mass, low surface brightness galaxies which typically have low dust content and anemic star formation, but are also likely to be important for evolution of the Milky Way-sized systems. We parameterize the dependencies observed in our simulations in convenient fitting formulae, which can be used to model the dependence of the KS relation on the dust-to-gas ratio and FUV flux in semi-analytic models and in cosmological simulations that do not include radiative transfer and H2 formation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Inelastic Dark Matter As An Efficient Fuel For Compact Stars

Dark matter in the form of weakly interacting massive particles is predicted to become gravitationally captured and accumulate in stars. While the subsequent annihilations of such particles lead to the injection of energy into stellar cores, elastically scattering dark matter particles do not generally yield enough energy to observably impact stellar phenomenology. Dark matter particles which scatter inelastically with nuclei (such that they reconcile the annual modulation reported by DAMA with the null results of CDMS and other experiments), however, can be captured by and annihilate in compact stars at a much higher rate. As a result, old white dwarf stars residing in high dark matter density environments can be prevented from cooling below several thousand degrees Kelvin. Observations of old, cool white dwarfs in dwarf spheroidal galaxies, or in the inner kiloparsec of the Milky Way, can thus potentially provide a valuable test of the inelastic dark matter hypothesis.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Measuring Gas Accretion and Angular Momentum near Simulated Supermassive Black Holes

Using cosmological simulations with a dynamic range in excess of 10 million, we study the transport of gas mass and angular momentum through the circumnuclear region of a disk galaxy containing a supermassive black hole (SMBH). The simulations follow fueling over relatively quiescent phases of the galaxy's evolution (no mergers) and without feedback from active galactic nuclei (AGNs), as part of the first stage of using state-of-the-art, high-resolution cosmological simulations to model galaxy and black hole co-evolution. We present results from simulations at different redshifts (z=6, 4, and 3) and three different black hole masses (30 million, 90 million, and 300 million solar masses; at z=4), as well as a simulation including a prescription that approximates optically thick cooling in the densest regions. The interior gas mass throughout the circumnuclear disk shows transient and chaotic behavior as a function of time. The Fourier transform of the interior gas mass follows a power law with slope -1 throughout the region, indicating that, in the absence of the effects of galaxy mergers and AGN feedback, mass fluctuations are stochastic with no preferred timescale for accretion over the duration of each simulation (~ 1-2 Myr). The angular momentum of the gas disk changes direction relative to the disk on kiloparsec scales over timescales less than 1 Myr, reflecting the chaotic and transient gas dynamics of the circumnuclear region. Infalling clumps of gas, which are driven inward as a result of the dynamical state of the circumnuclear disk, may play an important role in determining the spin evolution of an SMBH, as has been suggested in stochastic accretion scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modeling The Large Scale Bias of Neutral Hydrogen

We present new analytical estimates of the large scale bias of neutral hydrogen (HI). We use a simple, non-parametric model which monotonically relates the total mass of a halo M_tot with its HI mass M_HI at zero redshift; for earlier times we assume limiting models for the HI density evolution consistent with the data presently available, as well as two main scenarios for the evolution of our M_HI - M_tot relation. We find that both the linear and the first nonlinear bias terms exhibit a strong evolution with redshift, regardless of the specific limiting model assumed for the H I density over time. These analytical predictions are then shown to be consistent with measurements performed on the Millennium Simulation. Additionally, we show that this strong bias evolution does not sensibly affect the measurement of the HI power spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Computational Eulerian Hydrodynamics and Galilean Invariance

Eulerian hydrodynamical simulations are a powerful and popular tool for modeling fluids in astrophysical systems. In this work, we critically examine recent claims that these methods violate Galilean invariance of the Euler equations. We demonstrate that Eulerian hydrodynamics methods do converge to a Galilean-invariant solution, provided a well-defined convergent solution exists. Specifically, we show that numerical diffusion, resulting from diffusion-like terms in the discretized hydrodynamical equations solved by Eulerian methods, accounts for the effects previously identified as evidence for the Galilean non-invariance of these methods. These velocity-dependent diffusive terms lead to different results for different bulk velocities when the spatial resolution of the simulation is kept fixed, but their effect becomes negligible as the resolution of the simulation is increased to obtain a converged solution. In particular, we find that Kelvin-Helmholtz instabilities develop properly in realistic Eulerian calculations regardless of the bulk velocity provided the problem is simulated with sufficient resolution (a factor of 2-4 increase compared to the case without bulk flows for realistic velocities). Our results reiterate that high-resolution Eulerian methods can perform well and obtain a convergent solution, even in the presence of highly supersonic bulk flows.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Computer Simulations of Cosmic Reionization

The cosmic reionization of hydrogen was the last major phase transition in the evolution of the universe, which drastically changed the ionization and thermal conditions in the cosmic gas. To the best of our knowledge today, this process was driven by the ultra-violet radiation from young, star-forming galaxies and from first quasars. We review the current observational constraints on cosmic reionization, as well as the dominant physical effects that control the ionization of intergalactic gas. We then focus on numerical modeling of this process with computer simulations. Over the past decade, significant progress has been made in solving the radiative transfer of ionizing photons from many sources through the highly inhomogeneous distribution of cosmic gas in the expanding universe. With modern simulations, we have finally converged on a general picture for the reionization process, but many unsolved problems still remain in this young and exciting field of numerical cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Cosmological Radiative Transfer Comparison Project II: The Radiation-Hydrodynamic Tests

The development of radiation hydrodynamical methods that are able to follow gas dynamics and radiative transfer self-consistently is key to the solution of many problems in numerical astrophysics. Such fluid flows are highly complex, rarely allowing even for approximate analytical solutions against which numerical codes can be tested. An alternative validation procedure is to compare different methods against each other on common problems, in order to assess the robustness of the results and establish a range of validity for the methods. Previously, we presented such a comparison for a set of pure radiative transfer tests (i.e. for fixed, non-evolving density fields). This is the second paper of the Cosmological Radiative Transfer (RT) Comparison Project, in which we compare 9 independent RT codes directly coupled to gasdynamics on 3 relatively simple astrophysical hydrodynamics problems: (5) the expansion of an H II region in a uniform medium; (6) an ionization front (I-front) in a 1/r^2 density profile with a flat core, and (7), the photoevaporation of a uniform dense clump. Results show a broad agreement between the different methods and no big failures, indicating that the participating codes have reached a certain level of maturity and reliability. However, many details still do differ, and virtually every code has showed some shortcomings and has disagreed, in one respect or another, with the majority of the results. This underscores the fact that no method is universal and all require careful testing of the particular features which are most relevant to the specific problem at hand.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Modeling Molecular Hydrogen and Star Formation in Cosmological Simulations

We describe a phenomenological model for molecular hydrogen formation suited for applications in galaxy formation simulations, which includes on-equilibrium formation of molecular hydrogen on dust and approximate treatment of both its self-shielding and shielding by dust from the dissociating UV radiation. The model is applicable in simulations in which individual star forming regions - the giant molecular complexes - can be identified (resolution of tens of pc) and their mean internal density estimated reliably, even if internal structure is not resolved. In agreement with previous studies, calculations based on our model show that the transition from atomic to fully molecular phase depends primarily on the metallicity, which we assume is directly related to the dust abundance, and clumpiness of the interstellar medium. The clumpiness simply boosts the formation rate of molecular hydrogen, while dust serves both as a catalyst of molecular hydrogen formation and as an additional shielding from dissociating UV radiation. The upshot is that it is difficult to form fully-shielded giant molecular clouds while gas metallicity is low. However, once the gas is enriched to Z ~ 0.01-0.1 solar, the subsequent star formation and enrichment can proceed at a much faster rate. This may keep star formation efficiency in the low-mass, low-metallicity progenitors of galaxies very low for a certain period of time with the effect similar to a strong "feedback" mechanism. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On the Kennicutt-Schmidt relation of low-metallicity high-redshift galaxies

We present results of self-consistent, high-resolution cosmological simulations of galaxy formation at z~3. The simulations employ recently developed recipe for star formation based on the local abundance of molecular hydrogen, which is tracked self-consistently during the course of simulation. The phenomenological H2 formation model accounts for the effects of dissociating UV radiation of stars in each galaxy, as well as self-shielding and shielding of H2 by dust, and therefore allows us to explore effects of lower metallicities and higher UV fluxes prevalent in high redshift galaxies on their star formation. We compare stellar masses, metallicities, and star formation rates of the simulated galaxies to available observations of the Lyman Break Galaxies (LBGs) and find a reasonable agreement. We find that the Kennicutt-Schmidt (KS) relation exhibited by our simulated galaxies at z~3 is substantially steeper and has a lower amplitude than the z=0 relation at Sigma_gas < 100 Msun/pc^2. The predicted relation, however, is consistent with existing observational constraints for the z~3 Damped Lyman $α$ (DLA) and LBGs. Our tests show that the main reason for the difference from the local KS relation is lower metallicity of the ISM in high redshift galaxies. We discuss several implications of the metallicity-dependence of the KS relation for galaxy evolution and interpretation of observations.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

A New Constraint on the Escape Fraction in Distant Galaxies Using Gamma-ray Burst Afterglow Spectroscopy

We describe a new method to measure the escape fraction fesc of ionizing radiation from distant star-forming galaxies using the afterglow spectra of long-duration gamma-ray bursts (GRBs). Optical spectra of GRB afterglows allow us to evaluate the optical depth of the host ISM, according to the neutral hydrogen column density N(HI) observed along the sightlines toward the star-forming regions where the GRBs are found. Different from previous effort in searching for faint, transmitted Lyman continuum photons, our method is not subject to background subtraction uncertainties and does not require prior knowledge of either the spectral shape of the host galaxy population or the IGM Lya forest absorption along these GRB sightlines. Because most GRBs occur in sub-L_* galaxies, our study also offers the first constraint on fesc for distant low-mass galaxies that dominate the cosmic luminosity density. We have compiled a sample of 27 GRBs at redshift z>2 for which the underlying N(HI) in the host ISM are known. These GRBs together offer a statistical sampling of the integrated optical depth to ionizing photons along random sightlines from star-forming regions in the host galaxies, and allow us to estimate the mean escape fraction <fesc> averaged over different viewing angles. We find <fesc>=0.02\pm 0.02 and place a 95% c.l. upper limit <fesc> <= 0.075 for these hosts. We discuss possible biases of our approach and implications of the result. Finally, we propose to extend this technique for measuring <fesc> at z~0.2 using spectra of core-collapse supernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Escape of Ionizing Radiation from High Redshift Galaxies

We model the escape of ionizing radiation from high-redshift galaxies using high-resolution Adaptive Mesh Refinement N-body + hydrodynamics simulations. Our simulations include time-dependent and spatially-resolved transfer of ionizing radiation in three dimensions, including effects of dust absorption. For galaxies of total mass M > 10^11 Msun and star formation rates SFR ~ 1-5 Msun/yr, we find angular averaged escape fractions of 0.01-0.03 over the entire redshift interval studied (3<z<9). In addition, we find that the escape fraction varies by more than an order of magnitude along different lines-of-sight within individual galaxies, from the largest values near galactic poles to the smallest along the galactic disk. The escape fraction declines steeply at lower masses and SFR. We show that the low values of escape fractions are due to a small fraction of young stars located just outside the edge of HI disk. We compare our predicted escape fraction of ionizing photons with previous results, and find a general agreement with both other simulation results and available direct detection measurements at z ~ 3. We also compare our simulations with a novel method to estimate the escape fraction in galaxies from the observed distribution of neutral hydrogen column densities along the lines of sights to long duration gamma-ray bursts. Using this method we find escape fractions of the GRB host galaxies of 2-3%, consistent with our theoretical predictions. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Scaling Relations of Dwarf Galaxies without Supernova-Driven Winds

Nearby dwarf galaxies exhibit tight correlations between their global stellar and dynamical properties, such as circular velocity, mass-to-light ratio, stellar mass, surface brightness, and metallicity. Such correlations have often been attributed to gas or metal-rich outflows driven by supernova energy feedback to the interstellar medium. We use high-resolution cosmological simulations of high-redshift galaxies with and without energy feedback, as well as analytic modeling, to investigate whether the observed correlations can arise without supernova-driven outflows. We find that the simulated dwarf galaxies exhibit correlations similar to those observed as early as z~10, regardless of whether supernova feedback is included. We also show that the correlations can be well reproduced by our analytic model that accounts for realistic gas inflow but assumes no outflows, and star formation rate obeying the Kennicutt-Schmidt law with a critical density threshold. We argue that correlations in simulated galaxies arise due to the increasingly inefficient conversion of gas into stars in low-mass dwarf galaxies rather than supernova-driven outflows. We also show that the decrease of the observed effective yield in low-mass objects, often used as an indicator of gas and metal outflows, can be reasonably reproduced in our simulations without outflows. We show that this trend can arise if a significant fraction of metals in small galaxies is spread to the outer regions of the halo outside the stellar extent via mixing. In this case the effective yield can be significantly underestimated if only metals within the stellar radius are taken into account. Measurements of gas metallicity in the outskirts of gaseous disks of dwarfs would thus provide a key test of such explanation.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Star Formation in a Cosmological Simulation of Reionization

We study the luminosity functions of high-redshift galaxies in detailed hydrodynamic simulations of cosmic reionization, which are designed to reproduce the evolution of the Lyman-alpha forest between z=5 and z=6. We find that the luminosity functions and total stellar mass densities are in agreement with observations when plausible assumptions about reddenning at z=6 are made. Our simulations support the conclusion that stars alone reionized the universe.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

AGN Outflows and the Matter Power Spectrum

We have investigated the effects of AGN outflows on the amplitude of the matter power spectrum in a simple model of spherically symmetric outflows around realistically clustered AGN population. We find that two competing effects influence the matter power spectrum in two opposite directions. First, AGN outflows move baryons from high to low density regions, decreasing the amplitude of the matter power spectrum by up to 20%. Second, high clustering of the AGN transfers the power from small to larger scales. The exact balance between these two effects depends on the details of outflows on small scales, and quantitative estimates will require much more sophisticated modeling than presented here.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Cosmic Reionization Redux

We show that numerical simulations of reionization that resolve the Lyman Limit systems (and, thus, correctly count absorptions of ionizing photons) have converged to about 10% level for 5<z<6.2 and are in reasonable agreement (within 10%) with the SDSS data in this redshift interval. The SDSS data thus constraint the redshift of overlap of cosmic HII regions to z_{OVL} = 6.1+-0.15. At higher redshifts, the simulations are far from convergence on the mean Gunn-Peterson optical depth, but achieve good convergence for the mean neutral hydrogen fraction. The simulations that fit the SDSS data, however, do not have nearly enough resolution to resolve the earliest episodes of star formation, and are very far from converging on the precise value of the optical depth to Thompson scattering - any value between 6 and 10% is possible, depending on the convergence rate of the simulations and the fractional contribution of PopIII stars. This is generally consistent with the third-year WMAP results, but much higher resolution simulation are required to come up with the sufficiently precise value for the Thompson optical depth that can be statistically compared with the WMAP data.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Incorporating Radiative Cooling into a Cosmological Hydrodynamic Code

A possible inconsistency arising when a radiative cooling term is incorporated in a finite resolution self-gravitating hydrodynamic code is discussed. The inconsistency appears when the heating-cooling balance within the cooling and collapsing gas cloud is broken near the resolution limit of a numerical code. As the result, the cooling time of a fluid element increases enormously leading to the unphysical conclusion that the fluid element does not cool and is therefore stable against the collapse. A special cooling consistency condition is introduced which approximately restores the heating-cooling balance and leads to a numerical solution that closely mimics the exact (infinite resolution) solution.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Equation of State of the Photoionized Intergalactic Medium

We develop an efficient method to study the effects of reionization history on the temperature-density relation of the intergalactic medium in the low density limit (overdensity less than 5). It is applied to the study of photo-reionization models in which the amplitude, spectrum and onset epoch of the ionizing flux, as well as the cosmology, are systematically varied. We find that the mean temperature-density relation at z=2-4 is well approximated by a power-law equation of state for uniform reionization models. We derive analytical expressions for its evolution and exhibit its asymptotic behavior: it is found that for sufficiently early reionization, imprints of reionization history prior to z=10 on the temperature-density relation are washed out. In this limit the temperature at cosmic mean density is proportional to (Ω_b h/\sqrtΩ_0)^{1/1.7}. While the amplitude of the radiation flux at the ionizing frequency of HI is found to have a negligible effect on the temperature-density relation as long as the universe reionizes before z=5, the spectrum can change the overall temperature by about 20%, through variations in the abundances of helium species. However the slope of the mean equation of state is found to lie within a narrow range for all reionization models we study, where reionization takes place before z=5. We discuss the implications of these findings for the observational properties of the Lyman-alpha forest. In particular, uncertainties in the temperature of the intergalactic medium, due to the uncertain reionization history of our universe, introduces a 30% scaling in the amplitude of the column density distribution while the the slope of the distribution is only affected by about 5%. Finally, we discuss how a fluctuating ionizing field affects the above results. We argue that under

0 citations0 reviewsNo code linkedOpen access