Search NASA⌕ Search

Engineering topics

Stacy, Athena

Publications and source records attributed to Stacy, Athena.

Magnetic fields in the formation of the first stars – II. Results

ABSTRACT Beginning with cosmological initial conditions at z = 100, we simulate the effects of magnetic fields on the formation of Population III stars and compare our results with the predictions of Paper I. We use gadget-2 to follow the evolution of the system while the field is weak. We introduce a new method for treating kinematic fields by tracking the evolution of the deformation tensor. The growth rate in this stage of the simulation is lower than expected for diffuse astrophysical plasmas, which have a very low resistivity (high magnetic Prandtl number); we attribute this to the large numerical resistivity in simulations, corresponding to a magnetic Prandtl number of order unity. When the magnetic field begins to be dynamically significant in the core of the minihalo at z = 27, we map it on to a uniform grid and follow the evolution in an adaptive mesh refinement, MHD simulation in orion2. The non-linear evolution of the field in the orion2 simulation violates flux-freezing and is consistent with the theory proposed by Xu & Lazarian. The fields approach equipartition with kinetic energy at densities ∼1010–1012 cm−3. When the same calculation is carried out in orion2 with no magnetic fields, several protostars form, ranging in mass from ∼1 to 30 M⊙; with magnetic fields, only a single ∼30 M⊙ protostar forms by the end of the simulation. Magnetic fields thus suppress the formation of low-mass Pop III stars, yielding a top-heavy Pop III initial mass function and contributing to the absence of observed Pop III stars.

79 ASTRONOMY AND ASTROPHYSICS↗

The First Stars: A Low-Mass Formation Mode

We perform numerical simulations of the growth of a Population III stellar system under photodissociating feedback. We start from cosmological initial conditions at z = 100, self-consistently following the formation of a minihalo at z = 15 and the subsequent collapse of its central gas to high densities. The simulations resolve scales as small as approx. 1 AU, corresponding to gas densities of 10(exp 16)/cu cm. Using sink particles to represent the growing protostars, we evolve the stellar system for the next 5000 yr. We find that this emerging stellar group accretes at an unusually low rate compared with minihalos which form at earlier times (z = 20-30), or with lower baryonic angular momentum. The stars in this unusual system will likely reach masses ranging from <1Stellar Mass to approx. 5 Stellar Mass by the end of their main-sequence lifetimes, placing them in the mass range for which stars will undergo an asymptotic giant branch (AGB) phase. Based upon the simulation, we predict the rare existence of Population III stars that have survived to the present day and have been enriched by mass overflow from a previous AGB companion.

Low-Mass↗

Constraining the Statistics of Population III Binaries

We perform a cosmological simulation in order to model the growth and evolution of Population III (Pop III) stellar systems in a range of host minihalo environments. A Pop III multiple system forms in each of the ten minihaloes, and the overall mass function is top-heavy compared to the currently observed initial mass function in the Milky Way. Using a sink particle to represent each growing protostar, we examine the binary characteristics of the multiple systems, resolving orbits on scales as small as 20 AU. We find a binary fraction of approx. 36, with semi-major axes as large as 3000 AU. The distribution of orbital periods is slightly peaked at approx. < 900 yr, while the distribution of mass ratios is relatively flat. Of all sink particles formed within the ten minihaloes, approx. 50 are lost to mergers with larger sinks, and 50 of the remaining sinks are ejected from their star-forming disks. The large binary fraction may have important implications for Pop III evolution and nucleosynthesis, as well as the final fate of the first stars.

Statistics↗

Probing Pre-Galactic Metal Enrichment with High-Redshift Gamma-Ray Bursts

We explore high-redshift gamma-ray bursts (GRBs) as promising tools to probe pre-galactic metal enrichment. We utilize the bright afterglow of a Population III (Pop III) GRB exploding in a primordial dwarf galaxy as a luminous background source, and calculate the strength of metal absorption lines that are imprinted by the first heavy elements in the intergalactic medium (IGM). To derive the GRB absorption line diagnostics, we use an existing highly resolved simulation of the formation of a first galaxy which is characterized by the onset of atomic hydrogen cooling in a halo with virial temperature approximately greater than10(exp 4) K.We explore the unusual circumburst environment inside the systems that hosted Pop III stars, modeling the density evolution with the self-similar solution for a champagne flow. For minihalos close to the cooling threshold, the circumburst density is roughly proportional to (1 + z) with values of about a few cm(exp −3). In more massive halos, corresponding to the first galaxies, the density may be larger, n approximately greater than100 cm(exp −3). The resulting afterglow fluxes are weakly dependent on redshift at a fixed observed time, and may be detectable with the James Webb Space Telescope and Very Large Array in the near-IR and radio wavebands, respectively, out to redshift z approximately greater than 20. We predict that the maximum of the afterglow emission shifts from near-IR to millimeter bands with peak fluxes from mJy to Jy at different observed times. The metal absorption line signature is expected to be detectable in the near future. GRBs are ideal tools for probing the metal enrichment in the early IGM, due to their high luminosities and featureless power-law spectra. The metals in the first galaxies produced by the first supernova (SN) explosions are likely to reside in low-ionization stages (C II, O I, Si II and Fe II). We show that, if the afterglow can be observed sufficiently early, analysis of the metal lines may distinguish whether the first heavy elements were produced in a pair-instability supernova or a core-collapse (Type II) SN, thus constraining the initial mass function of the first stars.

Metal↗

Effect of Population III Multiplicity on Dark Star Formation

We numerically study the mutual interaction between dark matter (DM) and Population III (Pop III) stellar systems in order to explore the possibility of Pop III dark stars within this physical scenario. We perform a cosmological simulation, initialized at z approx. 100, which follows the evolution of gas and DM. We analyze the formation of the first mini halo at z approx. 20 and the subsequent collapse of the gas to densities of 10(exp 12)/cu cm. We then use this simulation to initialize a set of smaller-scale 'cut-out' simulations in which we further refine the DM to have spatial resolution similar to that of the gas. We test multiple DM density profiles, and we employ the sink particle method to represent the accreting star-forming region. We find that, for a range of DM configurations, the motion of the Pop III star-disk system serves to separate the positions of the protostars with respect to the DM density peak, such that there is insufficient DM to influence the formation and evolution of the protostars for more than approx. 5000 years. In addition, the star-disk system causes gravitational scattering of the central DM to lower densities, further decreasing the influence of DM over time. Any DM-powered phase of Pop III stars will thus be very short-lived for the typical multiple system, and DM will not serve to significantly prolong the life of Pop III stars.

Stacy, Athena↗