Source author record

Tom Abel

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

48works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

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

48 published item(s)

preprint2022arXiv

Detection of spatial clustering in the 1000 richest SDSS DR8 redMaPPer clusters with Nearest Neighbor distributions

Distances to the $k$-nearest-neighbor ($k$NN) data points from volume-filling query points are a sensitive probe of spatial clustering. Here we present the first application of $k$NN summary statistics to observational clustering measurement, using the 1000 richest redMaPPer clusters ($0.1\leqslant z\leqslant 0.3$) from the SDSS DR8 catalog. A clustering signal is defined as a difference in the cumulative distribution functions (CDFs) of $k$NN distances from fixed query points to the observed clusters versus a set of unclustered random points. We find that the $k=1,2$-NN CDFs of redMaPPer deviate significantly from the randoms' across scales of 35 to 155 Mpc, which is a robust signature of clustering. In addition to $k$NN, we also measure the two-point correlation function for the same set of redMaPPer clusters versus random points, which shows a noisier and less significant clustering signal within the same radial scales. Quantitatively, the $χ^2$ distribution for both the $k$NN-CDFs and the two-point correlation function measured on the randoms peak at $χ^2\sim 50$ (null hypothesis), whereas the $k$NN-CDFs ($χ^2\sim 300$, $p = 1.54\times 10^{-36}$) pick up a much more significant clustering signal than the two-point function ($χ^2\sim 100$, $p = 1.16\times 10^{-6}$) when measured on redMaPPer. Finally, the measured 3NN and 4NN CDFs deviate from the predicted $k=3, 4$-NN CDFs assuming an ideal Gaussian field, indicating a non-Gaussian clustering signal for redMaPPer clusters, although its origin might not be cosmological due to observational systematics. Therefore, $k$NN serves as a more sensitive probe of clustering complementary to the two point correlation function, providing a novel approach for constraining cosmology and galaxy-halo connection.

preprint2022arXiv

Modeling Nearest Neighbor distributions of biased tracers using Hybrid Effective Field Theory

We investigate the application of Hybrid Effective Field Theory (HEFT) -- which combines a Lagrangian bias expansion with subsequent particle dynamics from $N$-body simulations -- to the modeling of $k$-Nearest Neighbor Cumulative Distribution Functions ($k{\rm NN}$-${\rm CDF}$s) of biased tracers of the cosmological matter field. The $k{\rm NN}$-${\rm CDF}$s are sensitive to all higher order connected $N$-point functions in the data, but are computationally cheap to compute. We develop the formalism to predict the $k{\rm NN}$-${\rm CDF}$s of discrete tracers of a continuous field from the statistics of the continuous field itself. Using this formalism, we demonstrate how $k{\rm NN}$-${\rm CDF}$ statistics of a set of biased tracers, such as halos or galaxies, of the cosmological matter field can be modeled given a set of low-redshift HEFT component fields and bias parameter values. These are the same ingredients needed to predict the two-point clustering. For a specific sample of halos, we show that both the two-point clustering \textit{and} the $k{\rm NN}$-${\rm CDF}$s can be well-fit on quasi-linear scales ($\gtrsim 20 h^{-1}{\rm Mpc}$) by the second-order HEFT formalism with the \textit{same values} of the bias parameters, implying that joint modeling of the two is possible. Finally, using a Fisher matrix analysis, we show that including $k{\rm NN}$-${\rm CDF}$ measurements over the range of allowed scales in the HEFT framework can improve the constraints on $σ_8$ by roughly a factor of $3$, compared to the case where only two-point measurements are considered. Combining the statistical power of $k{\rm NN}$ measurements with the modeling power of HEFT, therefore, represents an exciting prospect for extracting greater information from small-scale cosmological clustering.

preprint2022arXiv

Tracer-Field Cross-Correlations with $k$-Nearest Neighbor Distributions

In astronomy and cosmology, significant effort is devoted to characterizing and understanding spatial cross-correlations between points - e.g. galaxy positions, high energy neutrino arrival directions, X-ray and AGN sources, and continuous field - e.g. weak lensing and Cosmic Microwave Background (CMB) maps. Recently, we introduced the $k$-nearest neighbor formalism to better characterize the clustering of discrete (point) datasets. Here we extend it to the point-field cross-correlation analysis. It combines $k$NN measurements of the point data set with measurements of the field smoothed on many scales. The resulting statistics are sensitive to all orders in the joint clustering of the points and the field. We demonstrate that this approach, unlike the 2-pt cross-correlation, can measure the statistical dependence of two datasets even when there are no linear (Gaussian) correlations. We further demonstrate that this framework is far more effective than the two-point function in detecting cross-correlations when the continuous field is contaminated by high levels of noise. For a particularly high level of noise, the cross-correlations between halos and the underlying matter field in a cosmological simulation, between $10h^{-1}{\rm Mpc}$ and $30h^{-1}{\rm Mpc}$, is detected at $>5σ$ significance using the technique presented here, when the two-point cross-correlation significance is $\sim 1σ$. Finally, we show that the $k$NN cross-correlations of halos and the matter field can be well-modeled on quasilinear scales by the Hybrid Effective Field Theory (HEFT) framework, with the same set of bias parameters as are used for the two-point cross-correlations. The substantial improvement in the statistical power of detecting cross-correlations with this method makes it a promising tool for various cosmological applications.

preprint2021arXiv

Nearest Neighbor distributions: new statistical measures for cosmological clustering

The use of summary statistics beyond the two-point correlation function to analyze the non-Gaussian clustering on small scales is an active field of research in cosmology. In this paper, we explore a set of new summary statistics -- the $k$-Nearest Neighbor Cumulative Distribution Functions ($k{\rm NN}$-${\rm CDF}$). This is the empirical cumulative distribution function of distances from a set of volume-filling, Poisson distributed random points to the $k$-nearest data points, and is sensitive to all connected $N$-point correlations in the data. The $k{\rm NN}$-${\rm CDF}$ can be used to measure counts in cell, void probability distributions and higher $N$-point correlation functions, all using the same formalism exploiting fast searches with spatial tree data structures. We demonstrate how it can be computed efficiently from various data sets - both discrete points, and the generalization for continuous fields. We use data from a large suite of $N$-body simulations to explore the sensitivity of this new statistic to various cosmological parameters, compared to the two-point correlation function, while using the same range of scales. We demonstrate that the use of $k{\rm NN}$-${\rm CDF}$ improves the constraints on the cosmological parameters by more than a factor of $2$ when applied to the clustering of dark matter in the range of scales between $10h^{-1}{\rm Mpc}$ and $40h^{-1}{\rm Mpc}$. We also show that relative improvement is even greater when applied on the same scales to the clustering of halos in the simulations at a fixed number density, both in real space, as well as in redshift space. Since the $k{\rm NN}$-${\rm CDF}$ are sensitive to all higher order connected correlation functions in the data, the gains over traditional two-point analyses are expected to grow as progressively smaller scales are included in the analysis of cosmological data.

preprint2021arXiv

On the single classical field description of interacting scalar fields

We test the degree to which interacting Bosonic systems can be approximated by a classical field as total occupation number is increased. This is done with our publicly available code repository, \href{https://github.com/andillio/QIBS}{QIBS}, a massively parallel solver for these systems. We use a number of toy models well studied in the literature and track when the classical field description admits quantum corrections, called the quantum breaktime. This allows us to test claims in the literature regarding the rate of convergence of these systems to the classical evolution. We test a number of initial conditions, including coherent states, number eigenstates, and field number states. We find that of these initial conditions, only number eigenstates do not converge to the classical evolution as occupation number is increased. We find that systems most similar to scalar field dark matter exhibit a logarithmic enhancement in the quantum breaktime with total occupation number. Systems with contact interactions or with field number state initial conditions, and linear dispersions, exhibit a power law enhancement. Finally, we find that the breaktime scaling depends on both model interactions and initial conditions.

preprint2020arXiv

Investigating the use of field solvers for simulating classical systems

We explore the use of field solvers as approximations of classical Vlasov-Poisson systems. This correspondence is investigated in both electrostatic and gravitational contexts. We demonstrate the ability of field solvers to be excellent approximations of problems with cold initial condition into the non linear regime. We also investigate extensions of the Schrödinger-Poisson system that employ multiple stacked cold streams, and the von Neumann-Poisson equation as methods that can successfully reproduce the classical evolution of warm initial conditions. We then discuss how appropriate simulation parameters need to be chosen to avoid interference terms, aliasing, and wave behavior in the field solver solutions. We present a series of criteria clarifying how parameters need to be chosen in order to effectively approximate classical solutions.

preprint2019arXiv

High-redshift Galaxy Formation with Self-consistently Modeled Stars and Massive Black Holes: Stellar Feedback and Quasar Growth

As computational resolution of modern cosmological simulations reach ever so close to resolving individual star-forming clumps in a galaxy, a need for "resolution-appropriate" physics for a galaxy-scale simulation has never been greater. To this end, we introduce a self-consistent numerical framework that includes explicit treatments of feedback from star-forming molecular clouds (SFMCs) and massive black holes (MBHs). In addition to the thermal supernovae feedback from SFMC particles, photoionizing radiation from both SFMCs and MBHs is tracked through full 3-dimensional ray tracing. A mechanical feedback channel from MBHs is also considered. Using our framework, we perform a state-of-the-art cosmological simulation of a quasar-host galaxy at z~7.5 for ~25 Myrs with all relevant galactic components such as dark matter, gas, SFMCs, and an embedded MBH seed of ~> 1e6 Ms. We find that feedback from SFMCs and an accreting MBH suppresses runaway star formation locally in the galactic core region. Newly included radiation feedback from SFMCs, combined with feedback from the MBH, helps the MBH grow faster by retaining gas that eventually accretes on to the MBH. Our experiment demonstrates that previously undiscussed types of interplay between gas, SFMCs, and a MBH may hold important clues about the growth and feedback of quasars and their host galaxies in the high-redshift Universe.

preprint2016arXiv

Non-Thermal Electron Energization from Magnetic Reconnection in Laser-Driven Plasmas

The possibility of studying non-thermal electron energization in laser-driven plasma experiments of magnetic reconnection is studied using two- and three-dimensional particle-in-cell simulations. It is demonstrated that non-thermal electrons with energies more than an order of magnitude larger than the initial thermal energy can be produced in plasma conditions currently accessible in the laboratory. Electrons are accelerated by the reconnection electric field, being injected at varied distances from the X-points, and in some cases trapped in plasmoids, before escaping the finite-sized system. Trapped electrons can be further energized by the electric field arising from the motion of the plasmoid. This acceleration gives rise to a non-thermal electron component that resembles a power-law spectrum, containing up to ~ 8% of the initial energy of the interacting electrons and ~ 24 % of the initial magnetic energy. Estimates of the maximum electron energy and of the plasma conditions required to observe suprathermal electron acceleration are provided, paving the way for a new platform for the experimental study of particle acceleration induced by reconnection.

preprint2015arXiv

An exact general remeshing scheme applied to physically conservative voxelization

We present an exact general remeshing scheme to compute analytic integrals of polynomial functions over the intersections between convex polyhedral cells of old and new meshes. In physics applications this allows one to ensure global mass, momentum, and energy conservation while applying higher-order polynomial interpolation. We elaborate on applications of our algorithm arising in the analysis of cosmological N-body data, computer graphics, and continuum mechanics problems. We focus on the particular case of remeshing tetrahedral cells onto a Cartesian grid such that the volume integral of the polynomial density function given on the input mesh is guaranteed to equal the corresponding integral over the output mesh. We refer to this as "physically conservative voxelization". At the core of our method is an algorithm for intersecting two convex polyhedra by successively clipping one against the faces of the other. This algorithm is an implementation of the ideas presented abstractly by Sugihara (1994), who suggests using the planar graph representations of convex polyhedra to ensure topological consistency of the output. This makes our implementation robust to geometric degeneracy in the input. We employ a simplicial decomposition to calculate moment integrals up to quadratic order over the resulting intersection domain. We also address practical issues arising in a software implementation, including numerical stability in geometric calculations, management of cancellation errors, and extension to two dimensions. In a comparison to recent work, we show substantial performance gains. We provide a C implementation intended to be a fast, accurate, and robust tool for geometric calculations on polyhedral mesh elements.

preprint2015arXiv

Cold accretion in early galaxy formation and its Lyman-alpha signatures

The Lyman-alpha (Lya) emission has played an important role in detecting high-redshift galaxies, including recently distant ones at redshift z > 7. It may also contain important information on the origin of these galaxies. Here, we investigate the formation of a typical L* galaxy and its observational signatures at the earliest stage, by combining a cosmological hydrodynamic simulation with three-dimensional radiative transfer calculations using the newly improved ART^2 code. Our cosmological simulation uses the Aquila initial condition which zooms in onto a Milky Way-like halo with high resolutions, and our radiative transfer couples multi-wavelength continuum, Lya line, and ionization of hydrogen. We find that the modeled galaxy starts to form at redshift z ~ 24 through efficient accretion of cold gas, which produces a strong Lya line with a luminosity of L(Lya) ~ 10^42 erg/s as early as z ~ 14. The Lya emission appears to trace the cold, dense gas. The lines exhibit asymmetric, single-peak profiles, and are shifted to the blue wing, a characteristic feature of gas inflow. Moreover, the contribution to the total Lya luminosity by excitation cooling increases with redshift, and it becomes dominant at z >~ 6. We predict that L* galaxies such as the modeled one may be detected at z <~ 8 by JWST and ALMA with a reasonable integration time. Beyond redshift 12, however, only Lya line may be observable by spectroscopic surveys. Our results suggest that Lya line is one of the most powerful tools to detect the first generation of galaxies, and to decipher their formation mechanism.

preprint2015arXiv

Simplex-in-Cell Technique for Collisionless Plasma Simulations

We extend the simplex-in-cell (SIC) technique recently introduced in the context of collisionless dark matter fluids (Abel et al. 2012; Hahn et al. 2012) to the case of collisionless plasmas. The six-dimensional phase space distribution function $f(\mathbf x,\mathbf v)$ is represented by an ensemble of three-dimensional manifolds, which we refer to as sheets. The electric potential field is obtained by solving the Poisson equation on a uniform mesh, where the charge density is evaluated by a spatial projection of the phase space sheets. The SIC representation of phase space density facilitates robust, high accuracy numerical evolution of the Vlasov-Poisson system using significantly fewer tracer particles than comparable particle-in-cell (PIC) approaches by reducing the numerical shot-noise associated with the latter. We introduce the SIC formulation and describe its implementation in a new code, which we validate using standard test problems including plasma oscillations, Landau damping, and two stream instabilities in one dimension. Merits of the new scheme are shown to include higher accuracy and faster convergence rates in the number of particles. We finally motivate and outline the efficient application of SIC to higher dimensional problems.

preprint2015arXiv

The Properties of Cosmic Velocity Fields

Understanding the velocity field is very important for modern cosmology: it gives insights to structure formation in general, and also its properties are crucial ingredients in modelling redshift-space distortions and in interpreting measurements of the kinetic Sunyaev-Zeldovich effect. Unfortunately, characterising the velocity field in cosmological N-body simulations is inherently complicated by two facts: i) The velocity field becomes manifestly multi-valued after shell-crossing and has discontinuities at caustics. This is due to the collisionless nature of dark matter. ii) N-body simulations sample the velocity field only at a set of discrete locations, with poor resolution in low-density regions. In this paper, we discuss how the associated problems can be circumvented by using a phase-space interpolation technique. This method provides extremely accurate estimates of the cosmic velocity fields and its derivatives, which can be properly defined without the need of the arbitrary "coarse-graining" procedure commonly used. We explore in detail the configuration-space properties of the cosmic velocity field on very large scales and in the highly nonlinear regime. In particular, we characterise the divergence and curl of the velocity field, present their one-point statistics, analyse the Fourier-space properties and provide fitting formulae for the velocity divergence bias relative to the non-linear matter power spectrum. We furthermore contrast some of the interesting differences in the velocity fields of warm and cold dark matter models. We anticipate that the high-precision measurements carried out here will help to understand in detail the dynamics of dark matter and the structures it forms.

preprint2014arXiv

Escape of Lyman-alpha and continuum photons from star-forming galaxies

A large number of high-redshift galaxies have been discovered via their narrow-band Lya line or broad-band continuum colors in recent years. The nature of the escaping process of photons from these early galaxies is crucial to understanding galaxy evolution and the cosmic reionization. Here, we investigate the escape of Lya, non-ionizing UV-continuum (l = 1300 - 1600 angstrom in rest frame), and ionizing photons (l < 912 angstrom) from galaxies by combining a cosmological hydrodynamic simulation with three-dimensional multi-wavelength radiative transfer calculations. The galaxies are simulated in a box of 5^3 h^-3 Mpc^3 with high resolutions using the Aquila initial condition which reproduces a Milky Way-like galaxy at redshift z=0. We find that the escape fraction (fesc) of these different photons shows a complex dependence on redshift and galaxy properties: fesc(Lya) and fesc(UV) appear to evolve with redshift, and they show similar, weak correlations with galaxy properties such as mass, star formation, metallicity, and dust content, while fesc(Ion) remains roughly constant at ~ 0.2 from z ~ 0 - 10, and it does not show clear dependence on galaxy properties. fesc(Lya) correlates more strongly with fesc(UV) than with fesc(Ion). In addition, we find a relation between the emergent Lya luminosity and the ionizing photon emissivity of Lyman Alpha Emitters (LAEs). By combining this relation with the observed luminosity functions of LAEs at different redshift, we estimate the contribution from LAEs to the reionization of intergalactic medium (IGM). Our result suggests that ionizing photons from LAEs alone are not sufficient to ionize IGM at z > 6, but they can maintain the ionization of IGM at z ~ 0 - 5.

preprint2014arXiv

The Birth of a Galaxy - III. Propelling reionisation with the faintest galaxies

Starlight from galaxies plays a pivotal role throughout the process of cosmic reionisation. We present the statistics of dwarf galaxy properties at z > 7 in haloes with masses up to 10^9 solar masses, using a cosmological radiation hydrodynamics simulation that follows their buildup starting with their Population III progenitors. We find that metal-enriched star formation is not restricted to atomic cooling ($T_{\rm vir} \ge 10^4$ K) haloes, but can occur in haloes down to masses ~10^6 solar masses, especially in neutral regions. Even though these smallest galaxies only host up to 10^4 solar masses of stars, they provide nearly 30 per cent of the ionising photon budget. We find that the galaxy luminosity function flattens above M_UV ~ -12 with a number density that is unchanged at z < 10. The fraction of ionising radiation escaping into the intergalactic medium is inversely dependent on halo mass, decreasing from 50 to 5 per cent in the mass range $\log M/M_\odot = 7.0-8.5$. Using our galaxy statistics in a semi-analytic reionisation model, we find a Thomson scattering optical depth consistent with the latest Planck results, while still being consistent with the UV emissivity constraints provided by Ly$α$ forest observations at z = 4-6.

preprint2013arXiv

A new approach to simulating collisionless dark matter fluids

Recently, we have shown how current cosmological N-body codes already follow the fine grained phase-space information of the dark matter fluid. Using a tetrahedral tesselation of the three-dimensional manifold that describes perfectly cold fluids in six-dimensional phase space, the phase-space distribution function can be followed throughout the simulation. This allows one to project the distribution function into configuration space to obtain highly accurate densities, velocities, and velocity dispersions. Here, we exploit this technique to show first steps on how to devise an improved particle-mesh technique. At its heart, the new method thus relies on a piecewise linear approximation of the phase space distribution function rather than the usual particle discretisation. We use pseudo-particles that approximate the masses of the tetrahedral cells up to quadrupolar order as the locations for cloud-in-cell (CIC) deposit instead of the particle locations themselves as in standard CIC deposit. We demonstrate that this modification already gives much improved stability and more accurate dynamics of the collisionless dark matter fluid at high force and low mass resolution. We demonstrate the validity and advantages of this method with various test problems as well as hot/warm-dark matter simulations which have been known to exhibit artificial fragmentation. This completely unphysical behaviour is much reduced in the new approach. The current limitations of our approach are discussed in detail and future improvements are outlined.

preprint2013arXiv

Enzo: An Adaptive Mesh Refinement Code for Astrophysics

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.

preprint2013arXiv

How closely do baryons follow dark matter on large scales?

We investigate the large-scale clustering and gravitational interaction of baryons and dark matter (DM) over cosmic time using a set of collisionless N-body simulations. Both components, baryons and DM, are evolved from distinct primordial density and velocity power spectra as predicted by early-universe physics. We first demonstrate that such two-component simulations require an unconventional match between force and mass resolution (i.e. force softening on at least the mean particle separation scale). Otherwise, the growth on any scale is not correctly recovered because of a spurious coupling between the two species at the smallest scales. With these simulations, we then demonstrate how the primordial differences in the clustering of baryons and DM are progressively diminished over time. In particular, we explicitly show how the BAO signature is damped in the spatial distribution of baryons and imprinted in that of DM. This is a rapid process, yet it is still not fully completed at low redshifts. On large scales, the overall shape of the correlation function of baryons and DM differs by 2% at z = 9 and by 0.2% at z = 0. The differences in the amplitude of the BAO peak are approximately a factor of 5 larger: 10% at z = 9 and 1% at z = 0. These discrepancies are, however, smaller than effects expected to be introduced by galaxy formation physics in both the shape of the power spectrum and in the BAO peak, and are thus unlikely to be detected given the precision of the next generation of galaxy surveys. Hence, our results validate the standard practice of modelling the observed galaxy distribution using predictions for the total mass clustering in the Universe.

preprint2013arXiv

Noiseless Gravitational Lensing Simulations

The microphysical properties of the DM particle can, in principle, be constrained by the properties and abundance of substructures in DM halos, as measured through strong gravitational lensing. Unfortunately, there is a lack of accurate theoretical predictions for the lensing signal of substructures, mainly because of the discreteness noise inherent to N-body simulations. Here we present Recursive-TCM, a method that is able to provide lensing predictions with an arbitrarily low discreteness noise, without any free parameters or smoothing scale. This solution is based on a novel way of interpreting the results of N-body simulations, where particles simply trace the evolution and distortion of Lagrangian phase-space volume elements. We discuss the advantages of this method over the widely used cloud-in-cells and adaptive-kernel smoothing density estimators. Applying the new method to a cluster-sized DM halo simulated in warm and cold DM scenarios, we show how the expected differences in their substructure population translate into differences in the convergence and magnification maps. We anticipate that our method will provide the high-precision theoretical predictions required to interpret and fully exploit strong gravitational lensing observations.

preprint2013arXiv

Reionization Histories of Milky Way Mass Halos

We investigate the connection between the epoch of reionization and the present day universe, by examining the extended mass reionization histories of dark matter halos identified at z=0. We combine an N-body dark matter simulation of a 600 Mpc volume with a three-dimensional, seminumerical reionization model. This provides reionization redshifts for each particle, which can then be connected with the properties of their halos at the present time. We find that the vast majority of present-day halos with masses larger than ~ few x 10^11 Msun reionize earlier than the rest of the universe. We also find significant halo-to-halo diversity in mass reionization histories, and find that in realistic inhomogenous models, the material within a given halo is not expected to reionize at the same time. In particular, the scatter in reionization times within individual halos is typically larger than the scatter among halos. From our fiducial reionization model, we find that the typical 68% scatter in reionization times within halos is ~ 115 Myr for 10^(12 \pm 0.25) Msun halos, decreasing slightly to ~ 95 Myr for 10^(15 \pm 0.25) Msun halos. We find a mild correlation between reionization history and environment: halos with shorter reionization histories are typically in more clustered environments, with the strongest trend on a scale of ~ 20 Mpc. Material in Milky Way mass halos with short reionization histories is preferentially reionized in relatively large HII regions, implying reionization mostly by sources external to the progenitors of the present-day halo. We investigate the impact on our results of varying the reionization model parameters, which span a range of reionization scenarios with varying timing and morphology.

preprint2013arXiv

The Warm DM halo mass function below the cut-off scale

Warm Dark Matter (WDM) cosmologies are a viable alternative to the Cold Dark Matter (CDM) scenario. Unfortunately, an accurate scrutiny of the WDM predictions with N-body simulations has proven difficult due to numerical artefacts. Here, we report on cosmological simulations that, for the first time, are devoid of those problems, and thus, are able to accurately resolve the WDM halo mass function well below the cut-off. We discover a complex picture, with perturbations at different evolutionary stages populating different ranges in the halo mass function. On the smallest mass scales we can resolve, identified objects are typically centres of filaments that are starting to collapse. On intermediate mass scales, objects typically correspond to fluctuations that have collapsed and are in the process of relaxation, whereas the high mass end is dominated by objects similar to haloes identified in CDM simulations. We then explicitly show how the formation of low-mass haloes is suppressed, which translates into a strong cut-off in the halo mass function. This disfavours some analytic formulations that predict a halo mass function that would extend well below the free streaming mass. We argue for a more detailed exploration of the formation of the smallest structures expected to form in a given cosmology, which, we foresee, will advance our overall understanding of structure formation.

preprint2012arXiv

A Novel Approach to Visualizing Dark Matter Simulations

In the last decades cosmological N-body dark matter simulations have enabled ab initio studies of the formation of structure in the Universe. Gravity amplified small density fluctuations generated shortly after the Big Bang, leading to the formation of galaxies in the cosmic web. These calculations have led to a growing demand for methods to analyze time-dependent particle based simulations. Rendering methods for such N-body simulation data usually employ some kind of splatting approach via point based rendering primitives and approximate the spatial distributions of physical quantities using kernel interpolation techniques, common in SPH (Smoothed Particle Hydrodynamics)-codes. This paper proposes three GPU-assisted rendering approaches, based on a new, more accurate method to compute the physical densities of dark matter simulation data. It uses full phase-space information to generate a tetrahedral tessellation of the computational domain, with mesh vertices defined by the simulation's dark matter particle positions. Over time the mesh is deformed by gravitational forces, causing the tetrahedral cells to warp and overlap. The new methods are well suited to visualize the cosmic web. In particular they preserve caustics, regions of high density that emerge, when several streams of dark matter particles share the same location in space, indicating the formation of structures like sheets, filaments and halos. We demonstrate the superior image quality of the new approaches in a comparison with three standard rendering techniques for N-body simulation data.

preprint2012arXiv

ART^2 : Coupling Lyman-alpha Line and Multi-wavelength Continuum Radiative Transfer

Narrow-band Lya line and broad-band continuum have played important roles in the discovery of high-redshift galaxies in recent years. Hence, it is crucial to study the radiative transfer of both Lya and continuum photons in the context of galaxy formation and evolution in order to understand the nature of distant galaxies. Here, we present a three-dimensional Monte Carlo radiative transfer code, All-wavelength Radiative Transfer with Adaptive Refinement Tree (ART^2), which couples Lya line and multi-wavelength continuum, for the study of panchromatic properties of galaxies and interstellar medium. This code is based on the original version of Li et al., and features three essential modules: continuum emission from X-ray to radio, Lya emission from both recombination and collisional excitation, and ionization of neutral hydrogen. The coupling of these three modules, together with an adaptive refinement grid, enables a self-consistent and accurate calculation of the Lya properties. As an example, we apply ART^2 to a cosmological simulation that includes both star formation and black hole growth, and study in detail a sample of massive galaxies at redshifts z=3.1 - 10.2. We find that these galaxies are Lya emitters (LAEs), whose Lya emission traces the dense gas region, and that their Lya lines show a shape characteristic of gas inflow. Furthermore, the Lya properties, including photon escape fraction, emergent luminosity, and equivalent width, change with time and environment. Our results suggest that LAEs evolve with redshift, and that early LAEs such as the most distant one detected at z ~ 8.6 may be dwarf galaxies with a high star formation rate fueled by infall of cold gas, and a low Lya escape fraction.

preprint2012arXiv

Single-Pass GPU-Raycasting for Structured Adaptive Mesh Refinement Data

Structured Adaptive Mesh Refinement (SAMR) is a popular numerical technique to study processes with high spatial and temporal dynamic range. It reduces computational requirements by adapting the lattice on which the underlying differential equations are solved to most efficiently represent the solution. Particularly in astrophysics and cosmology such simulations now can capture spatial scales ten orders of magnitude apart and more. The irregular locations and extensions of the refined regions in the SAMR scheme and the fact that different resolution levels partially overlap, poses a challenge for GPU-based direct volume rendering methods. kD-trees have proven to be advantageous to subdivide the data domain into non-overlapping blocks of equally sized cells, optimal for the texture units of current graphics hardware, but previous GPU-supported raycasting approaches for SAMR data using this data structure required a separate rendering pass for each node, preventing the application of many advanced lighting schemes that require simultaneous access to more than one block of cells. In this paper we present a single-pass GPU-raycasting algorithm for SAMR data that is based on a kD-tree. The tree is efficiently encoded by a set of 3D-textures, which allows to adaptively sample complete rays entirely on the GPU without any CPU interaction. We discuss two different data storage strategies to access the grid data on the GPU and apply them to several datasets to prove the benefits of the proposed method.

preprint2012arXiv

The Birth of a Galaxy. II. The Role of Radiation Pressure

Massive stars provide feedback that shapes the interstellar medium of galaxies at all redshifts and their resulting stellar populations. Here we present three adaptive mesh refinement radiation hydrodynamics simulations that illustrate the impact of momentum transfer from ionising radiation to the absorbing gas on star formation in high-redshift dwarf galaxies. Momentum transfer is calculated by solving the radiative transfer equation with a ray tracing algorithm that is adaptive in spatial and angular coordinates. We find that momentum input partially affects star formation by increasing the turbulent support to a three-dimensional rms velocity equal to the circular velocity of early haloes. Compared to a calculation that neglects radiation pressure, the star formation rate is decreased by a factor of five to 1.8 x 10^{-2} Msun/yr in a dwarf galaxy with a dark matter and stellar mass of 2.0 x 10^8 and 4.5 x 10^5 solar masses, respectively, when radiation pressure is included. Its mean metallicity of 10^{-2.1} Z_sun is consistent with the observed dwarf galaxy luminosity-metallicity relation. However, what one may naively expect from the calculation without radiation pressure, the central region of the galaxy overcools and produces a compact, metal-rich stellar population with an average metallicity of 0.3 Z_sun, indicative of an incorrect physical recipe. In addition to photo-heating in HII regions, radiation pressure further drives dense gas from star forming regions, so supernovae feedback occurs in a warmer and more diffuse medium, launching metal-rich outflows. Capturing this aspect and a temporal separation between the start of radiative and supernova feedback are numerically important in the modeling of galaxies to avoid the "overcooling problem". We estimate that dust in early low-mass galaxies is unlikely to aid in momentum transfer from radiation to the gas.

preprint2012arXiv

The Effect of Absorption Systems on Cosmic Reionization

We use large-scale simulations to investigate the morphology of reionization during the final, overlap phase. Our method uses an efficient three-dimensional smoothing technique which takes into account the finite mean free path due to absorption systems, lambda, by only smoothing over scales R_s<lambda. The large dynamic range of our calculations is necessary to resolve the neutral patches left at the end of reionization within a representative volume; we find that simulation volumes exceeding several hundred Mpc on a side are necessary in order to properly model reionization when the neutral fraction is ~0.01-0.3. Our results indicate a strong dependence of percolation morphology on a large and uncertain region of model parameter space. The single most important parameter is the mean free path to absorption systems, which serve as opaque barriers to ionizing radiation. If these absorption systems were as abundant as some realistic estimates indicate, the spatial structure of the overlap phase is considerably more complex than previously predicted. In view of the lack of constraints on the mean free path at the highest redshifts, current theories that do not include absorption by Lyman-limit systems, and in particular three-dimensional simulations, may underestimate the abundance of neutral clouds at the end of reionization. This affects predictions for the 21 cm signal associated with reionization, interpretation of absorption features in quasar spectra at z ~5-6, the connection between reionization and the local universe, and constraints on the patchiness and duration of reionization from temperature fluctuations measured in the cosmic microwave background arising from the kinetic Sunyaev-Zel'dovich effect.

preprint2012arXiv

Tracing the Dark Matter Sheet in Phase Space

The primordial velocity dispersion of dark matter is small compared to the velocities attained during structure formation. The initial density distribution is close to uniform and it occupies an initial sheet in phase space that is single valued in velocity space. Because of gravitational forces this three dimensional manifold evolves in phase space without ever tearing, conserving phase-space volume and preserving the connectivity of nearby points. N-body simulations already follow the motion of this sheet in phase space. This fact can be used to extract full fine-grained phase-space-structure information from existing cosmological N-body simulations. Particles are considered as the vertices of an unstructured three dimensional mesh, moving in six dimensional phase-space. On this mesh, mass density and momentum are uniquely defined. We show how to obtain the space density of the fluid, detect caustics, and count the number of streams as well as their individual contributions to any point in configuration-space. We calculate the bulk velocity, local velocity dispersions, and densities from the sheet - all without averaging over control volumes. This gives a wealth of new information about dark matter fluid flow which had previously been thought of as inaccessible to N-body simulations. We outline how this mapping may be used to create new accurate collisionless fluid simulation codes that may be able to overcome the sparse sampling and unphysical two-body effects that plague current N-body techniques.

preprint2012arXiv

Were progenitors of local L* galaxies Lyman-alpha emitters at high redshift?

The Lya emission has been observed from galaxies over a redshift span z ~ 0 - 8.6. However, the evolution of high-redshift Lya emitters (LAEs), and the link between these populations and local galaxies, remain poorly understood. Here, we investigate the Lya properties of progenitors of a local L* galaxy by combining cosmological hydrodynamic simulations with three-dimensional radiative transfer calculations using the new ART^2 code. We find that the main progenitor (the most massive one) of a Milky Way-like galaxy has a number of Lya properties close to those of observed LAEs at z ~ 2 - 6, but most of the fainter ones appear to fall below the detection limits of current surveys. The Lya photon escape fraction depends sensitively on a number of physical properties of the galaxy, such as mass, star formation rate, and metallicity, as well as galaxy morphology and orientation. Moreover, we find that high-redshift LAEs show blue-shifted Lya line profiles characteristic of gas inflow, and that the Lya emission by excitation cooling increases with redshift, and becomes dominant at z > 6. Our results suggest that some observed LAEs at z ~ 2-6 with luminosity of L_Lya ~ 10^{42-43} ergs/s may be similar to the main progenitor of the Milky Way at high redshift, and that they may evolve into present-day L* galaxies.

preprint2011arXiv

Comparing Numerical Methods for Isothermal Magnetized Supersonic Turbulence

We employ simulations of supersonic super-Alfvenic turbulence decay as a benchmark test problem to assess and compare the performance of nine astrophysical MHD methods actively used to model star formation. The set of nine codes includes: ENZO, FLASH, KT-MHD, LL-MHD, PLUTO, PPML, RAMSES, STAGGER, and ZEUS. We present a comprehensive set of statistical measures designed to quantify the effects of numerical dissipation in these MHD solvers. We compare power spectra for basic fields to determine the effective spectral bandwidth of the methods and rank them based on their relative effective Reynolds numbers. We also compare numerical dissipation for solenoidal and dilatational velocity components to check for possible impacts of the numerics on small-scale density statistics. Finally, we discuss convergence of various characteristics for the turbulence decay test and impacts of various components of numerical schemes on the accuracy of solutions. We show that the best performing codes employ a consistently high order of accuracy for spatial reconstruction of the evolved fields, transverse gradient interpolation, conservation law update step, and Lorentz force computation. The best results are achieved with divergence-free evolution of the magnetic field using the constrained transport method, and using little to no explicit artificial viscosity. Codes which fall short in one or more of these areas are still useful, but they must compensate higher numerical dissipation with higher numerical resolution. This paper is the largest, most comprehensive MHD code comparison on an application-like test problem to date. We hope this work will help developers improve their numerical algorithms while helping users to make informed choices in picking optimal applications for their specific astrophysical problems.

preprint2011arXiv

Enzo+Moray: Radiation Hydrodynamics Adaptive Mesh Refinement Simulations with Adaptive Ray Tracing

We describe a photon-conserving radiative transfer algorithm, using a spatially-adaptive ray tracing scheme, and its parallel implementation into the adaptive mesh refinement (AMR) cosmological hydrodynamics code, Enzo. By coupling the solver with the energy equation and non-equilibrium chemistry network, our radiation hydrodynamics framework can be utilised to study a broad range of astrophysical problems, such as stellar and black hole (BH) feedback. Inaccuracies can arise from large timesteps and poor sampling, therefore we devised an adaptive time-stepping scheme and a fast approximation of the optically-thin radiation field with multiple sources. We test the method with several radiative transfer and radiation hydrodynamics tests that are given in Iliev et al. (2006, 2009). We further test our method with more dynamical situations, for example, the propagation of an ionisation front through a Rayleigh-Taylor instability, time-varying luminosities, and collimated radiation. The test suite also includes an expanding H II region in a magnetised medium, utilising the newly implemented magnetohydrodynamics module in Enzo. This method linearly scales with the number of point sources and number of grid cells. Our implementation is scalable to 512 processors on distributed memory machines and can include radiation pressure and secondary ionisations from X-ray radiation. It is included in the newest public release of Enzo.

preprint2011arXiv

Magnetic Fields in Population III Star Formation

We study the buildup of magnetic fields during the formation of Population III star-forming regions, by conducting cosmological simulations from realistic initial conditions and varying the Jeans resolution. To investigate this in detail, we start simulations from identical initial conditions, mandating 16, 32 and 64 zones per Jeans length, and studied the variation in their magnetic field amplification. We find that, while compression results in some amplification, turbulent velocity fluctuations driven by the collapse can further amplify an initially weak seed field via dynamo action, provided there is sufficient numerical resolution to capture vortical motions (we find this requirement to be 64 zones per Jeans length, slightly larger than, but consistent with previous work run with more idealized collapse scenarios). We explore saturation of amplification of the magnetic field, which could potentially become dynamically important in subsequent, fully-resolved calculations. We have also identified a relatively surprising phenomena that is purely hydrodynamic: the higher-resolved simulations possess substantially different characteristics, including higher infall-velocity, increased temperatures inside 1000 AU, and decreased molecular hydrogen content in the innermost region. Furthermore, we find that disk formation is suppressed in higher-resolution calculations, at least at the times that we can follow the calculation. We discuss the effect this may have on the buildup of disks over the accretion history of the first clump to form as well as the potential for gravitational instabilities to develop and induce fragmentation.

preprint2011arXiv

Multi-scale initial conditions for cosmological simulations

We discuss a new algorithm to generate multi-scale initial conditions with multiple levels of refinements for cosmological "zoom-in" simulations. The method uses an adaptive convolution of Gaussian white noise with a real space transfer function kernel together with an adaptive multi-grid Poisson solver to generate displacements and velocities following first (1LPT) or second order Lagrangian perturbation theory (2LPT). The new algorithm achieves RMS relative errors of order 10^(-4) for displacements and velocities in the refinement region and thus improves in terms of errors by about two orders of magnitude over previous approaches. In addition, errors are localized at coarse-fine boundaries and do not suffer from Fourier-space induced interference ringing. An optional hybrid multi-grid and Fast Fourier Transform (FFT) based scheme is introduced which has identical Fourier space behaviour as traditional approaches. Using a suite of re-simulations of a galaxy cluster halo our real space based approach is found to reproduce correlation functions, density profiles, key halo properties and subhalo abundances with per cent level accuracy. Finally, we generalize our approach for two-component baryon and dark-matter simulations and demonstrate that the power spectrum evolution is in excellent agreement with linear perturbation theory. For initial baryon density fields, it is suggested to use the local Lagrangian approximation in order to generate a density field for mesh based codes that is consistent with Lagrangian perturbation theory instead of the current practice of using the Eulerian linearly scaled densities.

preprint2011arXiv

The Birth of a Galaxy: Primordial Metal Enrichment and Stellar Populations

By definition, Population III stars are metal-free, and their protostellar collapse is driven by molecular hydrogen cooling in the gas-phase, leading to large characteristic masses. Population II stars with lower characteristic masses form when the star-forming gas reaches a critical metallicity of 10^{-6} - 10^{-3.5} Z_\odot. We present an adaptive mesh refinement radiation hydrodynamics simulation that follows the transition from Population III to II star formation. The maximum spatial resolution of 1 comoving parsec allows for individual molecular clouds to be well-resolved and their stellar associations to be studied in detail. We model stellar radiative feedback with adaptive ray tracing. A top-heavy initial mass function for the Population III stars is considered, resulting in a plausible distribution of pair-instability supernovae and associated metal enrichment. We find that the gas fraction recovers from 5 percent to nearly the cosmic fraction in halos with merger histories rich in halos above 10^7 solar masses. A single pair-instability supernova is sufficient to enrich the host halo to a metallicity floor of 10^{-3} Z_\odot and to transition to Population II star formation. This provides a natural explanation for the observed floor on damped Lyman alpha (DLA) systems metallicities reported in the literature, which is of this order. We find that stellar metallicities do not necessarily trace stellar ages, as mergers of halos with established stellar populations can create superpositions of t-Z evolutionary tracks. A bimodal metallicity distribution is created after a starburst occurs when the halo can cool efficiently through atomic line cooling.

preprint2011arXiv

Weak Primordial Magnetic Fields and Anisotropies in the Cosmic Microwave Background Radiation

It is shown that small-scale magnetic fields present before recombination induce baryonic density inhomogeneities of appreciable magnitude. The presence of such inhomogeneities changes the ionization history of the Universe, which in turn decreases the angular scale of the Doppler peaks and increases Silk damping by photon diffusion. This unique signature could be used to (dis)prove the existence of primordial magnetic fields of strength as small as B~10^(-11) Gauss by upcoming cosmic microwave background observations.

preprint2010arXiv

A Multi-Code Analysis Toolkit for Astrophysical Simulation Data

The analysis of complex multiphysics astrophysical simulations presents a unique and rapidly growing set of challenges: reproducibility, parallelization, and vast increases in data size and complexity chief among them. In order to meet these challenges, and in order to open up new avenues for collaboration between users of multiple simulation platforms, we present yt (available at http://yt.enzotools.org/), an open source, community-developed astrophysical analysis and visualization toolkit. Analysis and visualization with yt are oriented around physically relevant quantities rather than quantities native to astrophysical simulation codes. While originally designed for handling Enzo's structure adaptive mesh refinement (AMR) data, yt has been extended to work with several different simulation methods and simulation codes including Orion, RAMSES, and FLASH. We report on its methods for reading, handling, and visualizing data, including projections, multivariate volume rendering, multi-dimensional histograms, halo finding, light cone generation and topologically-connected isocontour identification. Furthermore, we discuss the underlying algorithms yt uses for processing and visualizing data, and its mechanisms for parallelization of analysis tasks.

preprint2010arXiv

Effects of Varying the Three-Body Molecular Hydrogen Formation Rate in Primordial Star Formation

The transformation of atomic hydrogen to molecular hydrogen through three-body reactions is a crucial stage in the collapse of primordial, metal-free halos, where the first generation of stars (Population III stars) in the Universe are formed. However, in the published literature, the rate coefficient for this reaction is uncertain by nearly an order of magnitude. We report on the results of both adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) simulations of the collapse of metal-free halos as a function of the value of this rate coefficient. For each simulation method, we have simulated a single halo three times, using three different values of the rate coefficient. We find that while variation between halo realizations may be greater than that caused by the three-body rate coefficient being used, both the accretion physics onto Population III protostars as well as the long-term stability of the disk and any potential fragmentation may depend strongly on this rate coefficient.

preprint2010arXiv

High-Entropy Polar Regions Around the First Protostars

We report on simulations of the formation of the first stars in the Universe, where we identify regions of hot atomic gas (fH2 < 10-6) at densities above 10-14 g/cc, heated to temperatures ranging between 3000 and 8000 K. Within this temperature range atomic hydrogen is unable to cool effectively. We describe the kinetic and thermal characteristics of these regions and investigate their origin. We find that these regions, while small in total mass fraction of the cloud, may be dynamically important over the accretion timescale for the central clump in the cloud, particularly as a chemical, rather than radiative, mechanism for clearing the polar regions of the accretion disk of material and terminating accretion along these directions. These inherently three-dimensional effects stress the need for multi-dimensional calculations of protostellar accretion for reliable predictions of the masses of the very first stars.

preprint2010arXiv

Lowering the Characteristic Mass of Cluster Stars by Magnetic Fields and Outflow Feedback

Magnetic fields are generally expected to increase the characteristic mass of stars formed in stellar clusters, because they tend to increase the effective Jeans mass. We test this expectation using adaptive mesh refinement (AMR) magnetohydrodynamic simulations of cluster formation in turbulent magnetized clumps of molecular clouds, treating stars as accreting sink particles. We find that, contrary to the common expectation, a magnetic field of strength in the observed range decreases, rather than increases, the characteristic stellar mass. It (1) reduces the number of intermediate-mass stars that are formed through direct turbulent compression, because sub-regions of the clump with masses comparable to those of stars are typically magnetically subcritical and cannot be compressed directly into collapse, and (2) increases the number of low-mass stars that are produced from the fragmentation of dense filaments. The filaments result from mass accumulation along the field lines. In order to become magnetically supercritical and fragment, the filament must accumulate a large enough column density (proportional to the field strength), which yields a high volume density (and thus a small thermal Jeans mass) that is conducive to forming low-mass stars. We find, in addition, that the characteristic stellar mass is reduced further by outflow feedback. The conclusion is that both magnetic fields and outflow feedback are important in shaping the stellar initial mass function (IMF).

preprint2010arXiv

rpSPH: a novel Smoothed Particle Hydrodynamics Algorithm

We suggest a novel discretisation of the momentum equation for Smoothed Particle Hydrodynamics (SPH) and show that it significantly improves the accuracy of the obtained solutions. Our new formulation which we refer to as relative pressure SPH, rpSPH, evaluates the pressure force in respect to the local pressure. It respects Newtons first law of motion and applies forces to particles only when there is a net force acting upon them. This is in contrast to standard SPH which explicitly uses Newtons third law of motion continuously applying equal but opposite forces between particles. rpSPH does not show the unphysical particle noise, the clumping or banding instability, unphysical surface tension, and unphysical scattering of different mass particles found for standard SPH. At the same time it uses fewer computational operations. and only changes a single line in existing SPH codes. We demonstrate its performance on isobaric uniform density distributions, uniform density shearing flows, the Kelvin-Helmholtz and Rayleigh-Taylor instabilities, the Sod shock tube, the Sedov-Taylor blast wave and a cosmological integration of the Santa Barbara galaxy cluster formation test. rpSPH is an improvement these cases. The improvements come at the cost of giving up exact momentum conservation of the scheme. Consequently one can also obtain unphysical solutions particularly at low resolutions.

preprint2009arXiv

Adaptive Mesh Fluid Simulations on GPU

We describe an implementation of compressible inviscid fluid solvers with block-structured adaptive mesh refinement on Graphics Processing Units using NVIDIA's CUDA. We show that a class of high resolution shock capturing schemes can be mapped naturally on this architecture. Using the method of lines approach with the second order total variation diminishing Runge-Kutta time integration scheme, piecewise linear reconstruction, and a Harten-Lax-van Leer Riemann solver, we achieve an overall speedup of approximately 10 times faster execution on one graphics card as compared to a single core on the host computer. We attain this speedup in uniform grid runs as well as in problems with deep AMR hierarchies. Our framework can readily be applied to more general systems of conservation laws and extended to higher order shock capturing schemes. This is shown directly by an implementation of a magneto-hydrodynamic solver and comparing its performance to the pure hydrodynamic case. Finally, we also combined our CUDA parallel scheme with MPI to make the code run on GPU clusters. Close to ideal speedup is observed on up to four GPUs.

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.

preprint2009arXiv

Outflow Feedback Regulated Massive Star Formation in Parsec-Scale Cluster Forming Clumps

(Abridged) We investigate massive star formation in turbulent, magnetized, parsec-scale clumps of molecular clouds including protostellar outflow feedback using Enzo-based MHD simulations with accreting sink particles and effective resolution $2048^3$. We find that, in the absence of regulation by magnetic fields and outflow feedback, massive stars form readily in a turbulent, moderately condensed clump of $\sim 1,600$ solar masses, along with a cluster of hundreds of lower mass stars. The massive stars are fed at high rates by (1) transient dense filaments produced by large-scale turbulent compression at early times, and (2) by the clump-wide global collapse resulting from turbulence decay at late times. In both cases, the bulk of the massive star's mass is supplied from outside a 0.1 pc-sized "core" that surrounds the star. In our simulation, the massive star is clump-fed rather than core-fed. The need for large-scale feeding makes the massive star formation prone to regulation by outflow feedback, which directly opposes the feeding processes. The outflows reduce the mass accretion rates onto the massive stars by breaking up the dense filaments that feed the massive star formation at early times, and by collectively slowing down the global collapse that fuel the massive star formation at late times. The latter is aided by a moderate magnetic field of strength in the observed range. We conclude that the massive star formation in our simulated turbulent, magnetized, parsec-scale clump is outflow-regulated and clump-fed (ORCF for short). An important implication is that the formation of low-mass stars in a dense clump can affect the formation of massive stars in the same clump, through their outflow feedback on the clump dynamics.

preprint2009arXiv

The Impact of Inhomogeneous Reionization on the Satellite Galaxy Population of the Milky Way

We use the publicly available subhalo catalogs from the Via Lactea simulation along with a Gpc-scale N-body simulation to understand the impact of inhomogeneous reionization on the satellite galaxy population of the Milky Way. The large-volume simulation is combined with a model for reionization that allows us to predict the distribution of reionization times for Milky Way mass halos. Motivated by this distribution, we identify candidate satellite galaxies in the simulation by requiring that any subhalo must grow above a specified mass threshold before it is reionized; after this time the photoionizing background will suppress both the formation of stars and the accretion of gas. We show that varying the reionization time over the range expected for Milky Way mass halos can change the number of satellite galaxies by roughly two orders of magnitude. This conclusion is in contradiction with a number of studies in the literature, and we conclude that this is a result of inconsistent application of the results of Gnedin (2000). We compare our satellite galaxies to observations using both abundance matching and stellar population synthesis methods to assign luminosities to our subhalos and account for observational completeness effects. Additionally, if we assume that the mass threshold is set by the virial temperature Tvir = 8e3K we find that our model accurately matches the vmax distribution, radial distribution, and luminosity function of observed Milky Way satellites for a reionization time zreion = 9.6^{1.0}_{-2.1}, assuming that the Via Lacteasubhalo distribution is representative of the Milky Way. This results in the presence of 119^{+202}_{-50} satellite galaxies.

preprint2008arXiv

First Stars III Conference Summary

The understanding of the formation, life, and death of Population III stars, as well as the impact that these objects had on later generations of structure formation, is one of the foremost issues in modern cosmological research and has been an active area of research during the past several years. We summarize the results presented at "First Stars III," a conference sponsored by Los Alamos National Laboratory, the Kavli Institute for Particle Astrophysics and Cosmology, and the Joint Institute for Nuclear Astrophysics. This conference, the third in a series, took place in July 2007 at the La Fonda Hotel in Santa Fe, New Mexico, U.S.A.

preprint2006arXiv

Formation of Primordial Stars in a LCDM Universe

We study the formation of the first generation of stars in the standard cold dark matter model, using a very high-resolution hydordynamic simulations. Our simulation achieves a dynamic range of 10^{10} in length scale. With accurate treatment of atomic and molecular physics, it allows us to study the chemo-thermal evolution of primordial gas clouds to densities up to n = 10^{16}/cc without assuming any a priori equation of state; a six orders of magnitudes improvement over previous three-dimensional calculations. All the relevant atomic and molecular cooling and heating processes, including cooling by collision-induced continuum emission, are implemented. For calculating optically thick H2 cooling at high densities, we use the Sobolev method. To examine possible gas fragmentation owing to thermal instability, we compute explicitly the growth rate of isobaric perturbations. We show that the cloud core does not fragment in either the low-density or high-density regimes. We also show that the core remains stable against gravitational deformation and fragmentation. We obtain an accurate gas mass accretion rate within a 10 Msun innermost region around the protostar. The protostar is accreting the surrounding hot gas at a rate of 0.001-0.01 Msun/yr. From these findings we conclude that primordial stars formed in early minihalos are massive. We carry out proto-stellar evolution calculations using the obtained accretion rate. The resulting mass of the first star is M_ZAMS = 60-100 Msun, with the exact mass dependent on the actual accretion rate.

preprint2006arXiv

The Role of Primordial Kicks on Black Hole Merger Rates

Primordial stars are likely to be very massive $\geq30\Msun$, form in isolation, and will likely leave black holes as remnants in the centers of their host dark matter halos in the mass range $10^{6}-10^{10}\Ms$. Such early black holes, at redshifts z$\gtsim10$, could be the seed black holes for the many supermassive black holes found in galaxies in the local universe. If they exist, their mergers with nearby supermassive black holes may be a prime signal for long wavelength gravitational wave detectors. We simulate formation of black holes in the center of high redshift dark matter halos and explore implications of initial natal kick velocities conjectured by some formation models. The central concentration of early black holes in present day galaxies is reduced if they are born even with moderate kicks of tens of km/s. The modest kicks allow the black holes to leave their parent halo, which consequently leads to dynamical friction being less effective on the lower mass black holes as compared to those still embedded in their parent halos. Therefore, merger rates may be reduced by more than an order of magnitude. Using analytical and illustrative cosmological N--body simulations we quantify the role of natal kicks of black holes formed from massive metal free stars on their merger rates with supermassive black holes in present day galaxies. Our results also apply to black holes ejected by the gravitational slingshot mechanism.

preprint2005arXiv

The Formation of Primordial Luminous Objects

In these lecture notes we review the current knowledge about the formation of the first luminous objects. We start from the cosmological context of hierarchical models of structure formation, and discuss the main physical processes which are believed to lead to primordial star formation, i.e. the cooling processes and the chemistry of molecules (especially H2) in a metal-free gas. We then describe the techniques and results of numerical simulations, which indicate that the masses of the first luminous objects are likely to be much larger than that of present-day stars. Finally, we discuss the scenario presented above, exposing some of the most interesting problems which are currently being investigated, such as that of the feedback effects of these objects.

preprint2001arXiv

The Formation of the First Star in the Universe

We describe results from a fully self-consistent three dimensional hydrodynamical simulation of the formation of one of the first stars in the Universe. Dark matter dominated pre-galactic objects form because of gravitational instability from small initidal density perturbations. As they assemble via hierarchical merging, primordial gas cools through ro-vibrational lines of hydrogen molecules and sinks to the center of the dark matter potential well. The high redshift analog of a molecular cloud is formed. When the dense, central parts of the cold gas cloud become self-gravitating, a dense core of approximately 100 solar mass undergoes rapid contraction. At densities n>10^9 cm^-3 a one solar mass proto-stellar core becomes fully molecular due to three-body H_2 formation. Contrary to analytical expectations this process does not lead to renewed fragmentation and only one star is formed. The calculation is stopped when optical depth effects become important, leaving the final mass of the fully formed star somewhat uncertain. At this stage the protostar is acreting material very rapidly (~0.01 solar masses per year). Radiative feedback from the star will not only halt its growth but also inhibit the formation of other stars in the same pre-galactic object (at least until the first star ends its life, presumably as a supernova). We conclude that at most one massive (M >> 1 solar mass) metal free star forms per pre-galactic halo, consistent with recent abundance measurements of metal poor galactic halo stars.

preprint1999arXiv

The Radiative Feedback of the First Cosmological Objects

In hierarchical models of structure formation, an early cosmic UV background (UVB) is produced by the small (T_vir < 10^4 K) halos that collapse before reionization. The UVB at energies below 13.6eV suppresses the formation of stars or black holes inside small halos, by photo-dissociating their only cooling agent, molecular H2. We self-consistently compute the buildup of the early UVB in Press-Schechter models, coupled with H2 photo-dissociation both in the intergalactic medium (IGM), and inside virialized halos. We find that the intergalactic H2 has a negligible effect on the UVB, both because its initial optical depth is small (tau<0.1), and because it is photo-dissociated at an early stage. If the UV sources in the first collapsed halos are stars, then their UV flux suppresses further star-formation inside small halos. This results in a pause in the buildup of the UVB, and reionization is delayed until larger halos (T_vir> 10^4 K) collapse. If the small halos host mini-quasars with hard spectra extending to approximately 1 keV, then their X-rays balance the effects of the UVB, the negative feedback does not occur, and reionization can be caused by the small halos.