• We model dust formation in the core collapse supernova explosion SN 1987A by treating the gas-phase formation of dust grain nuclei as a chemical process. To compute the synthesis of fourteen species of grains we integrate a non-equilibrium network of nucleating and related chemical reactions and follow the growth of the nuclei into grains via accretion and coagulation. The effects of the radioactive cobalt, titanium, and sodium on the thermodynamics and chemistry of the ejecta are taken into account. The grain temperature, which we allow to differ from the gas temperature, affects the surface-tension-corrected evaporation rate. We also account for He$^+$, Ne$^+$, Ar$^+$, and O weathering. We combine our dust synthesis model with a crude prescription for anisotropic radioactive nickel dredge-up into the core ejecta, the so-called `nickel bubbles', to compute the total dust mass and molecular-species-specific grain size distribution. The total mass varies between $0.41\,M_\odot$ and $0.73\,M_\odot$, depending on the bubble shell density contrast. In the decreasing order of abundance, the grain species produced are: magnesia, silicon, forsterite, iron sulfide, carbon, silicon dioxide, alumina, and iron. The combined grain size distribution is a power law $dN/da\propto a^{-4.39}$. Early ejecta compaction by expanding radioactive nickel bubbles strongly enhances dust synthesis. This underscores the need for improved understanding of hydrodynamic transport and mixing over the entire pre-homologous expansion.
  • Due to its accuracy and generality, Monte Carlo radiative transfer (MCRT) has emerged as the prevalent method for Ly$\alpha$ radiative transfer in arbitrary geometries. The standard MCRT encounters a significant efficiency barrier in the high optical depth, diffusion regime. Multiple acceleration schemes have been developed to improve the efficiency of MCRT but the noise from photon packet discretization remains a challenge. The discrete diffusion Monte Carlo (DDMC) scheme has been successfully applied in state-of-the-art radiation hydrodynamics (RHD) simulations. Still, the established framework is not optimal for resonant line transfer. Inspired by the DDMC paradigm, we present a novel extension to resonant DDMC (rDDMC) in which diffusion in space and frequency are treated on equal footing. We explore the robustness of our new method and demonstrate a level of performance that justifies incorporating the method into existing Ly$\alpha$ codes. We present computational speedups of $\sim 10^2$-$10^6$ relative to contemporary MCRT implementations with schemes that skip scattering in the core of the line profile. This is because the rDDMC runtime scales with the spatial and frequency resolution rather than the number of scatterings - the latter is typically $\propto \tau_0$ for static media, or $\propto (a \tau_0)^{2/3}$ with core-skipping. We anticipate new frontiers in which on-the-fly Ly$\alpha$ radiative transfer calculations are feasible in 3D RHD. More generally, rDDMC is transferable to any computationally demanding problem amenable to a Fokker-Planck approximation of frequency redistribution.
  • The physics of star formation at its extreme, in the nuclei of the densest and the most massive star clusters in the universe - potential massive black hole nurseries - has for decades eluded scrutiny. Spectroscopy of these systems has been scarce, whereas theoretical arguments suggest that radiation pressure on dust grains somehow inhibits star formation. Here, we harness an accelerated Monte Carlo radiation transport scheme to report a radiation hydrodynamical simulation of super star cluster formation in turbulent clouds. We find that radiation pressure reduces the global star formation efficiency by 30-35%, and the star formation rate by 15-50%, both relative to a radiation-free control run. Overall, radiation pressure does not terminate the gas supply for star formation and the final stellar mass of the most massive cluster is $\sim 1.3\times10^6\,M_\odot$. The limited impact as compared to in idealized theoretical models is attributed to a radiation-matter anti-correlation in the supersonically turbulent, gravitationally collapsing medium. In isolated regions outside massive clusters, where the gas distribution is less disturbed, radiation pressure is more effective in limiting star formation. The resulting stellar density at the cluster core is $\ge 10^{8}\,M_{\odot}\,\textrm{pc}^{-3}$, with stellar velocity dispersion $\gtrsim 70\,\text{km}\,\text{s}^{-1}$. We conclude that the super star cluster nucleus is propitious to the formation of very massive stars via dynamical core collapse and stellar merging. We speculate that the very massive star may avoid the claimed catastrophic mass loss by continuing to accrete dense gas condensing from a gravitationally-confined ionized phase.
  • Extremely metal poor stars have been the focus of much recent attention owing to the expectation that their chemical abundances can shed light on the metal and dust yields of the earliest supernovae. We present our most realistic simulation to date of the astrophysical pathway to the first metal enriched stars. We simulate the radiative and supernova hydrodynamic feedback of a $60\,M_\odot$ Population III star starting from cosmological initial conditions realizing Gaussian density fluctuations. We follow the gravitational hydrodynamics of the supernova remnant at high spatial resolution through its freely-expanding, adiabatic, and radiative phases, until gas, now metal-enriched, has resumed runaway gravitational collapse. Our findings are surprising: while the Population III progenitor exploded with a low energy of $10^{51}\,\text{erg}$ and injected an ample metal mass of $6\,M_\odot$, the first cloud to collapse after the supernova explosion is a dense surviving primordial cloud on which the supernova blastwave deposited metals only superficially, in a thin, unresolved layer. The first metal-enriched stars can form at a very low metallicity, of only $2-5\times10^{-4}\,Z_\odot$, and can inherit the parent cloud's highly elliptical, radially extended orbit in the dark matter gravitational potential.
  • Recent observations have shown that the characteristic luminosity of the rest-frame ultraviolet (UV) luminosity function does not significantly evolve at 4 < z < 7 and is approximately M*_UV ~ -21. We investigate this apparent non-evolution by examining a sample of 178 bright, M_UV < -21 galaxies at z=4 to 7, analyzing their stellar populations and host halo masses. Including deep Spitzer/IRAC imaging to constrain the rest-frame optical light, we find that M*_UV galaxies at z=4-7 have similar stellar masses of log(M/Msol)=9.6-9.9 and are thus relatively massive for these high redshifts. However, bright galaxies at z=4-7 are less massive and have younger inferred ages than similarly bright galaxies at z=2-3, even though the two populations have similar star formation rates and levels of dust attenuation. We match the abundances of these bright z=4-7 galaxies to halo mass functions from the Bolshoi Lambda-CDM simulation to estimate the halo masses. We find that the typical halo masses in ~M*_UV galaxies decrease from log(M_h/Msol)=11.9 at z=4 to log(M_h/Msol)=11.4 at z=7. Thus, although we are studying galaxies at a similar mass across multiple redshifts, these galaxies live in lower mass halos at higher redshift. The stellar baryon fraction in units of the cosmic mean Omega_b/Omega_m rises from 5.1% at z=4 to 11.7% at z=7; this evolution is significant at the ~3-sigma level. This rise does not agree with simple expectations of how galaxies grow, and implies that some effect, perhaps a diminishing efficiency of feedback, is allowing a higher fraction of available baryons to be converted into stars at high redshifts.
  • The elegance of inflationary cosmology and cosmological perturbation theory ends with the formation of the first stars and galaxies, the initial sources of light that launched the phenomenologically rich process of cosmic reionization. Here we review the current understanding of early star formation, emphasizing unsolved problems and technical challenges. We begin with the first generation of stars to form after the Big Bang and trace how they influenced subsequent star formation. The onset of chemical enrichment coincided with a sharp increase in the overall physical complexity of star forming systems. Ab-initio computational treatments are just now entering the domain of the predictive and are establishing contact with local observations of the relics of this ancient epoch.
  • We investigate the formation of a galaxy reaching a virial mass of $~ 10^8$ solar mass at $z=10$ by carrying out a zoomed radiation-hydrodynamical cosmological simulation. This simulation traces Population~III (Pop~III) star formation, characterized by a modestly top-heavy initial mass function (IMF), and considers stellar feedback such as photoionization heating from Pop III and Population~II (Pop~II) stars, mechanical and chemical feedback from supernovae (SNe), and X-ray feedback from accreting black holes (BHs) and high-mass X-ray binaries (HMXBs). We self-consistently impose a transition in star formation mode from top-heavy Pop III to low-mass Pop~II, and find that the star formation rate in the computational box is dominated by Pop~III until $z=13$, and by Pop~II thereafter. The simulated galaxy experiences bursty star formation, with a substantially reduced gas content due to photoionization heating from Pop~III and Pop~II stars, together with SN feedback. All the gas within the simulated galaxy is metal-enriched above $10^{-5}$ solar, such that there are no remaining pockets of primordial gas. The simulated galaxy has an estimated observed flux of $~10^{-3} nJy$, which is too low to be detected by the James Webb Space Telescope (JWST) without strong lensing amplification. We also show that our simulated galaxy is similar in terms of stellar mass to Segue 2, the least luminous dwarf known in the Local Group.
  • Radiation pressure can be dynamically important in star-forming environments such as ultra-luminous infrared and submillimeter galaxies. Whether and how radiation drives turbulence and bulk outflows in star formation sites is still unclear. The uncertainty in part reflects the limitations of direct numerical schemes that are currently used to simulate radiation transfer and radiation-gas coupling. An idealized setup in which radiation is introduced at the base of a dusty atmosphere in a gravitational field has recently become the standard test for radiation-hydrodynamics methods in the context of star formation. To a series of treatments featuring the flux-limited-diffusion approximation as well as a short-characteristics tracing and M1 closure for the variable Eddington tensor approximation, we here add another treatment that is based on the Implicit Monte Carlo radiation transfer scheme. Consistent with all previous treatments, the atmosphere undergoes Rayleigh-Taylor instability and readjusts to a near-Eddington-limited state. We detect late-time net acceleration in which the turbulent velocity dispersion matches that reported previously with the short-characteristics-based radiation transport closure, the most accurate of the three preceding treatments. Our technical result demonstrates the importance of accurate radiation transfer in simulations of radiative feedback.
  • We present a simulation of the long-term evolution of a Population III supernova remnant in a cosmological minihalo. Employing passive Lagrangian tracer particles, we investigate how chemical stratification and anisotropy in the explosion can affect the abundances of the first low-mass, metal-enriched stars. We find that reverse shock heating can leave the inner mass shells at entropies too high to cool, leading to carbon-enhancement in the re-collapsing gas. This hydrodynamic selection effect could explain the observed incidence of carbon-enhanced metal-poor (CEMP) stars at low metallicity. We further explore how anisotropic ejecta distributions, recently seen in direct numerical simulations of core-collapse explosions, may translate to abundances in metal-poor stars. We find that some of the observed scatter in the Population II abundance ratios can be explained by an incomplete mixing of supernova ejecta, even in the case of only one contributing enrichment event. We demonstrate that the customary hypothesis of fully-mixed ejecta clearly fails if post-explosion hydrodynamics prefers the recycling of some nucleosynthetic products over others. Furthermore, to fully exploit the stellar-archaeological program of constraining the Pop III initial mass function from the observed Pop II abundances, considering these hydrodynamical transport effects is crucial. We discuss applications to the rich chemical structure of ultra-faint dwarf satellite galaxies, to be probed in unprecedented detail with upcoming spectroscopic surveys.
  • To constrain the properties of the first stars with the chemical abundance patterns observed in metal-poor stars, one must identify any non-trivial effects that the hydrodynamics of metal dispersal can imprint on the abundances. We use realistic cosmological hydrodynamic simulations to quantify the distribution of metals resulting from one Population III supernova and from a small number of such supernovae exploding in close succession. Overall, supernova ejecta are highly inhomogeneously dispersed throughout the simulations. When the supernova bubbles collapse, quasi-virialized metal-enriched clouds, fed by fallback from the bubbles and by streaming of metal-free gas from the cosmic web, grow in the centers of the dark matter halos. Partial turbulent homogenization on scales resolved in the simulation is observed only in the densest clouds where the vortical time scales are short enough to ensure true homogenization on subgrid scales. However, the abundances in the clouds differ from the gross yields of the supernovae. Continuing the simulations until the cloud have gone into gravitational collapse, we predict that the abundances in second-generation stars will be deficient in the innermost mass shells of the supernova (if only one has exploded) or in the ejecta of the latest supernovae (when multiple have exploded). This indicates that hydrodynamics gives rise to biases complicating the identification of nucleosynthetic sources in the chemical abundance spaces of the surviving stars.
  • It is widely recognized that nucleosynthetic output of the first, Population III supernovae was a catalyst defining the character of subsequent stellar generations. Most of the work on the earliest enrichment was carried out assuming that the first stars were extremely massive and that the associated supernovae were unusually energetic, enough to completely unbind the baryons in the host cosmic minihalo and disperse the synthesized metals into the intergalactic medium. Recent work, however, suggests that the first stars may in fact have been somewhat less massive, with a characteristic mass scale of a few tens of solar masses. We present a cosmological simulation following the transport of the metals synthesized in a Population III supernova assuming that it had an energy of 1e51 ergs, compatible with standard Type II supernovae. A young supernova remnant is inserted in the first star's relic HII region in the free expansion phase and is followed for 40 Myr employing adaptive mesh refinement and Lagrangian tracer particle techniques. The supernova remnant remains partially trapped within the minihalo and the thin snowplow shell develops pronounced instability and fingering. Roughly half of the ejecta turn around and fall back toward the center of the halo, with 1% of the ejecta reaching the center in 30 kyr and 10% in 10 Myr. The average metallicity of the combined returning ejecta and the pristine filaments feeding into the halo center from the cosmic web is 0.001 - 0.01 Z_sun, but the two remain unmixed until accreting onto the central hydrostatic core that is unresolved at the end of the simulation. We conclude that if Population III stars had less extreme masses, they promptly enriched the host minihalos with metals and triggered Population II star formation.
  • We present the Cosmic Lyman-$\alpha$ Transfer code (COLT), a massively parallel Monte-Carlo radiative transfer code, to simulate Lyman-$\alpha$ (Ly$\alpha$) resonant scattering through neutral hydrogen as a probe of the first galaxies. We explore the interaction of centrally produced Ly$\alpha$ radiation with the host galactic environment. Ly$\alpha$ photons emitted from the luminous starburst region escape with characteristic features in the line profile depending on the density distribution, ionization structure, and bulk velocity fields. For example, anisotropic ionization exhibits a tall peak close to line centre with a skewed tail that drops off gradually. Idealized models of first galaxies explore the effect of mass, anisotropic H II regions, and radiation pressure driven winds on Ly$\alpha$ observables. We employ mesh refinement to resolve critical structures. We also post-process an ab initio cosmological simulation and examine images captured at various escape distances within the 1 Mpc$^3$ comoving volume. Finally, we discuss the emergent spectra and surface brightness profiles of these objects in the context of high-$z$ observations. The first galaxies will likely be observed through the red damping wing of the Ly$\alpha$ line. Observations will be biased toward galaxies with an intrinsic red peak located far from line centre that reside in extensive H II super bubbles, which allows Hubble flow to sufficiently redshift photons away from line centre and facilitate transmission through the intergalactic medium (IGM). Even with gravitational lensing to boost the luminosity this preliminary work indicates that Ly$\alpha$ emission from stellar clusters within haloes of $M_{\rm vir}<10^9~{\rm M}_\odot$ is generally too faint to be detected by the James Webb Space Telescope (JWST).
  • We simulate the formation of a low metallicity (0.01 Zsun) stellar cluster in a dwarf galaxy at redshift z~14. Beginning with cosmological initial conditions, the simulation utilizes adaptive mesh refinement and sink particles to follow the collapse and evolution of gas past the opacity limit for fragmentation, thus resolving the formation of individual protostellar cores. A time- and location-dependent protostellar radiation field, which heats the gas by absorption on dust, is computed by integration of protostellar evolutionary tracks with the MESA code. The simulation also includes a robust non-equilibrium chemical network that self-consistently treats gas thermodynamics and dust-gas coupling. The system is evolved for 18 kyr after the first protostellar source has formed. In this time span, 30 sink particles representing protostellar cores form with a total mass of 81 Msun. Their masses range from ~0.1 Msun to 14.4 Msun with a median mass ~0.5-1 Msun. Massive protostars grow by competitive accretion while lower-mass protostars are stunted in growth by close encounters and many-body ejections. In the regime explored here, the characteristic mass scale is determined by the temperature floor set by the cosmic microwave background and by the onset of efficient dust-gas coupling. It seems unlikely that host galaxies of the first bursts of metal-enriched star formation will be detectable with the James Webb Space Telescope or other next-generation infrared observatories. Instead, the most promising access route to the dawn of cosmic star formation may lie in the scrutiny of metal-poor, ancient stellar populations in the Galactic neighborhood. The observable targets that correspond to the system simulated here are ultra-faint dwarf satellite galaxies such as Bootes II, Segue I and II, and Willman I.
  • We use cosmological simulations to assess how the explosion of the first stars in supernovae (SNe) influences early cosmic history. Specifically, we investigate the impact by SNe on the host systems for Population~III (Pop~III) star formation and explore its dependence on halo environment and Pop~III progenitor mass. We then trace the evolution of the enriched gas until conditions are met to trigger second-generation star formation. To this extent, we quantify the recovery timescale, which measures the time delay between a Pop~III SN explosion and the appearance of cold, dense gas, out of which second-generation stars can form. We find that this timescale is highly sensitive to the Pop~III progenitor mass, and less so to the halo environment. For more massive progenitors, including those exploding in pair instability SNe, second-generation star formation is delayed significantly, for up to a Hubble time. The dependence of the recovery time on the mass of the SN progenitor is mainly due to the ionizing impact of the progenitor star. Photoionization heating increases the gas pressure and initiates a hydrodynamical response that reduces the central gas density, an effect that is stronger in more massive. The gas around lower mass Pop~III stars remains denser and hence the SN remnants cool more rapidly, facilitating the subsequent re-condensation of the gas and formation of a second generation of stars. In most cases, the second-generation stars are already metal-enriched to ~2-5 X 10^{-4}\zsun, thus belonging to Population~II. The recovery timescale is a key quantity governing the nature of the first galaxies, able to host low-mass, long-lived stellar systems. These in turn are the target of future deep-field campaigns with the James Webb Space Telescope.
  • Recent work suggests that the first generation of stars, the so-called Population III (Pop III), could have formed primarily in binaries or as members of small multiple systems. Here we investigate the impact of X-ray feedback from High-Mass X-ray Binaries (HMXBs) left behind in stellar binary systems after the primary forms a black hole (BH), accreting gas at a high rate from the companion, a process that is thought to be favored at the low metallicities characteristic of high-redshift gas. Thanks to their large mean free path, X-rays are capable of preionizing and preheating the gas in the intergalactic medium (IGM) and in haloes long before the reionization of the Universe is complete, and thus could have strongly affected the formation of subsequent generations of stars as well as reionization. We have carried out zoomed hydrodynamical cosmological simulations of minihaloes, accounting for the formation of Pop III stars and their collapse into BHs and HMXBs, and the associated radiation-hydrodynamic feedback from UV and X-ray photons. We find no strong net feedback from HMXBs on the simulated star formation history. On the other hand, the preheating of the IGM by HMXBs leads to a strong suppression of small-scale structures and significantly lowers the recombination rate in the IGM, thus yielding a net positive feedback on reionization. We further show that X-ray feedback from HMXBs can augment the ionizing feedback from the Pop III progenitor stars to suppress gas accretion onto the first BHs, limiting their growth into supermassive BHs. Finally, we show that X-ray ionization by HMXBs leaves distinct signatures in the properties of the high-redshift hydrogen that may be probed in upcoming observations of the redshifted 21cm spin-flip line.
  • Dwarf spheroidal satellite galaxies have emerged a powerful probe of small-scale dark matter clustering and of cosmic reionization. They exhibit structural and chemical continuity with dwarf irregular galaxies in the field and with spheroidal galaxies in high-density environments. By combining empirical constraints derived for star formation at low gas column densities and metallicities in the local universe with a model for dark matter and baryonic mass assembly, we provide an analytical description of how the dwarf spheroidals acquired their stellar content. Their progenitors formed stars until the gas content, initially reduced from the cosmic average by the thermal pressure of the reionized intergalactic medium, was finally ram pressure stripped during the progenitors' accretion on to the host galaxy. Dwarf spheroidal satellites of differing luminosities seem to share very similar most massive progenitor histories that reach thresholds for gas cooling by atomic line emission at epochs at which the Lagrangian volume of the Local Group should have been reionized. We hypothesize that dwarf spheroidals formed the bulk of their stars in partially rotationally supported HI disks in a reionized universe. This model provides an explanation for the "common mass scale" relation and reproduces the empirical luminosity-size and luminosity-metallicity relations. Explosive feedback phenomena, such as outflows driven by the concerted action of supernovae, need not have been significant in the dwarf spheroidals' formation. We further speculate that the true pre-reionization fossils should exhibit a structure distinct from that of the dwarf spheroidals, e.g., in the form of dense isolated or nuclear star clusters.
  • We simulate the formation of a metal-poor (10^-2 Zsun) stellar cluster in one of the first galaxies to form in the early Universe, specifically a high-redshift atomic cooling halo (z~14). This is the first calculation that resolves the formation of individual metal-enriched stars in simulations starting from realistic cosmological initial conditions. We follow the evolution of a single dense clump among several in the parent halo. The clump forms a cluster of ~40 stars and sub-stellar objects within 7000 years and could continue forming stars ~5 times longer. Protostellar dust heating has a negligible effect on the star formation efficiency, at least during the early evolutionary stages, but it moderately suppresses gaseous fragmentation and brown dwarf formation. We observe fragmentation in thin gaseous filaments and sustained accretion in larger, rotating structures as well as ejections by binary interactions. The stellar initial mass function above 0.1 Msun, evaluated after ~10^4 years of fragmentation and accretion, seems in agreement with the recent measurement in ultra-faint dwarf spheroidal Galactic satellites of Geha et al. (2013).
  • Population III stars are believed to have been more massive than typical stars today and to have formed in relative isolation. The thermodynamic impact of metals is expected to induce a transition leading to clustered, low-mass Population II star formation. In this work, we present results from three cosmological simulations, only differing in gas metallicity, that focus on the impact of metal fine-structure line cooling on the formation of stellar clusters in a high-redshift atomic cooling halo. Introduction of sink particles allows us to follow the process of gas hydrodynamics and accretion onto cluster stars for 4 Myr corresponding to multiple local free-fall times. At metallicities at least $10^{-3}\, Z_{\odot}$, gas is able to reach the CMB temperature floor and fragment pervasively resulting in a stellar cluster of size $\sim1$ pc and total mass $\sim1000\, M_{\odot}$. The masses of individual sink particles vary, but are typically $\sim100\, M_{\odot}$, consistent with the Jeans mass when gas cools to the CMB temperature, though some solar mass fragments are also produced. At the low metallicity of $10^{-4}\, Z_{\odot}$, fragmentation is completely suppressed on scales greater than 0.01 pc and total stellar mass is lower by a factor of 3 than in the higher metallicity simulations. The sink particle accretion rates, and thus their masses, are determined by the mass of the gravitationally unstable gas cloud and the prolonged gas accretion over many Myr. The simulations thus exhibit features of both monolithic collapse and competitive accretion. Even considering possible dust induced fragmentation that would occur at higher densities, the formation of a bona fide stellar cluster seems to require metal line cooling and metallicities of at least $10^{-3}\, Z_{\odot}$.
  • We investigate how radiative feedback from the first stars affects the assembly of the first dwarf galaxies. We perform cosmological zoomed SPH simulations of a dwarf galaxy assembling inside a halo of virial mass 10^9 solar at z = 10. The simulations follow the non-equilibrium chemistry/cooling of primordial gas and the conversion of the gas into metal-free stars. To quantify the radiative feedback, we compare a simulation in which stars emit both molecular hydrogen dissociating and hydrogen/helium ionizing radiation with a simulation in which stars emit only dissociating radiation, and with a simulation in which stars remain dark. Photodissociation and -ionization exert a strong negative feedback on the assembly of the simulated galaxy. Gas condensation is strongly impeded, and star formation is strongly suppressed in comparison with the simulation in which stars remain dark. The feedback on the gas implies a suppression of the central dark matter densities in the minihalo progenitor by factors of up to a few, which is a significant deviation from the singular isothermal density profile characterizing the dark matter distribution in the absence of radiative feedback. The evolution of gas densities, star formation rates, and the distribution of dark matter becomes insensitive to the inclusion of dissociating radiation in the late stages of the minihalo assembly, and it becomes insensitive to the inclusion of ionizing radiation once the minihalo turns into an atomically cooling galaxy. The formation of an extended disk inside the dwarf galaxy is a robust outcome not affected by the inclusion of radiation. We estimate that dwarf galaxies such as simulated here will be among the faintest galaxies the upcoming James Webb Space Telescope will detect. Our conclusions are subject to our neglect of feedback from supernovae and chemical enrichment as well as to cosmic variance. [abridged]
  • We investigate guide-field magnetic reconnection and particle acceleration in relativistic pair plasmas with three-dimensional particle-in-cell (PIC) simulations of a kinetic-scale current sheet in a periodic geometry at low magnetizations. The tearing instability is the dominant mode in the current sheet for all guide field strengths, while the linear kink mode is less important even without guide field. Oblique modes seem to be suppressed entirely. In its nonlinear evolution, the reconnection layer develops a network of interconnected and interacting magnetic flux ropes. As smaller flux ropes merge into larger ones, the reconnection layer evolves toward a three-dimensional, disordered state in which the resulting flux rope segments contain magnetic substructure on plasma skin depth scales. Embedded in the flux ropes, we detect spatially and temporally intermittent sites of dissipation reflected in peaks in the parallel electric field. Magnetic dissipation and particle acceleration persist until the end of the simulations, with simulations with higher magnetization and lower guide field strength exhibiting greater and faster energy conversion and particle energization. At the end of our largest simulation, the particle energy spectrum attains a tail extending to high Lorentz factors that is best modeled with a combination of two additional thermal components. We confirm that the primary energization mechanism is acceleration by the electric field in the X-line region. We discuss the implications of our results for macroscopic reconnection sites, and which of our results may be expected to hold in systems with higher magnetizations.
  • We investigate the process of metal-free star formation in the first galaxies with a high-resolution cosmological simulation. We consider the cosmologically motivated scenario in which a strong molecule-destroying Lyman-Werner (LW) background inhibits effective cooling in low-mass haloes, delaying star formation until the collapse or more massive haloes. Only when molecular hydrogen (H2) can self-shield from LW radiation, which requires a halo capable of cooling by atomic line emission, will star formation be possible. To follow the formation of multiple gravitationally bound objects, at high gas densities we introduce sink particles which accrete gas directly from the computational grid. We find that in a 1 Mpc^3 (comoving) box, runaway collapse first occurs in a 3x10^7 M_sun dark matter halo at z~12 assuming a background intensity of J21=100. Due to a runaway increase in the H2 abundance and cooling rate, a self-shielding, supersonically turbulent core develops abruptly with ~10^4 M_sun in cold gas available for star formation. We analyze the formation of this self-shielding core, the character of turbulence, and the prospects for star formation. Due to a lack of fragmentation on scales we resolve, we argue that LW-delayed metal-free star formation in atomic cooling haloes is very similar to star formation in primordial minihaloes, although in making this conclusion we ignore internal stellar feedback. Finally, we briefly discuss the detectability of metal-free stellar clusters with the James Webb Space Telescope.
  • Theoretical models predict that some of the first stars ended their lives as extremely energetic pair-instability supernovae (PISNe). With energies approaching 10^53 ergs, these supernovae are expected to be within the detection limits of the upcoming James Webb Space telescope (JWST), allowing observational constraints to be placed on the properties of the first stars. We estimate the source density of PISNe using a semi-analytic halo mass function based approach, accounting for the effects of feedback from star formation on the PISN rates using cosmological simulations. We estimate an upper limit of ~0.2 PISNe per JWST field of view at any given time. Feedback can reduce this rate significantly, e.g., lowering it to as little as one PISN per 4000 JWST fields of view for the most pessimistic explosion models. We also find that the main obstacle to observing PISNe from the first stars is their scarcity, not their faintness; exposures longer than a few times 10^4s will do little to increase the number of PISNe found. Given this we suggest a mosaic style search strategy for detecting PISNe from the first stars. Even rather high redshift PISNe are unlikely to be missed by moderate exposures, and a large number of pointings will be required to ensure a detection.
  • We study how the first galaxies were assembled under feedback from the accretion onto a central black hole (BH) that is left behind by the first generation of metal-free stars through self-consistent, cosmological simulations. X-ray radiation from the accretion of gas onto BH remnants of Population III (Pop III) stars, or from high-mass X-ray binaries (HMXBs), again involving Pop III stars, influences the mode of second generation star formation. We track the evolution of the black hole accretion rate and the associated X-ray feedback starting with the death of the Pop III progenitor star inside a minihalo and following the subsequent evolution of the black hole as the minihalo grows to become an atomically cooling galaxy. We find that X-ray photoionization heating from a stellar-mass BH is able to quench further star formation in the host halo at all times before the halo enters the atomic cooling phase. X-ray radiation from a HMXB, assuming a luminosity close to the Eddington value, exerts an even stronger, and more diverse, feedback on star formation. It photoheats the gas inside the host halo, but also promotes the formation of molecular hydrogen and cooling of gas in the intergalactic medium and in nearby minihalos, leading to a net increase in the number of stars formed at early times. Our simulations further show that the radiative feedback from the first BHs may strongly suppress early BH growth, thus constraining models for the formation of supermassive BHs.
  • Observational evidence suggests a link between long duration gamma ray bursts (LGRBs) and Type Ic supernovae. Here, we propose a potential mechanism for Type Ic supernovae in LGRB progenitors powered solely by accretion energy. We present spherically-symmetric hydrodynamic simulations of the long-term accretion of a rotating gamma-ray burst progenitor star, a "collapsar," onto the central compact object, which we take to be a black hole. The simulations were carried out with the adaptive mesh refinement code FLASH in one spatial dimension and with rotation, an explicit shear viscosity, and convection in the mixing length theory approximation. Once the accretion flow becomes rotationally supported outside of the black hole, an accretion shock forms and traverses the stellar envelope. Energy is carried from the central geometrically thick accretion disk to the stellar envelope by convection. Energy losses through neutrino emission and nuclear photodisintegration are calculated but do not seem important following the rapid early drop of the accretion rate following circularization. We find that the shock velocity, energy, and unbound mass are sensitive to convective efficiency, effective viscosity, and initial stellar angular momentum. Our simulations show that given the appropriate combinations of stellar and physical parameters, explosions with energies ~5x10^50 ergs, velocities 3000 km/s, and unbound material masses >6 solar masses are possible in a rapidly rotating 16 solar mass main sequence progenitor star. Further work is needed to constrain the values of these parameters, to identify the likely outcomes in more plausible and massive LRGB progenitors, and to explore nucleosynthetic implications.
  • The association of long-duration gamma-ray bursts (LGRBs) with Type Ic supernovae presents a challenge to supernova explosion models. In the collapsar model for LGRBs, gamma rays are produced in an ultrarelativistic jet launching from the magnetosphere of the black hole that forms in the aftermath of the collapse of a rotating progenitor star. The jet is collimated along the star's rotation axis, but the concomitant luminous supernova should be relatively--though certainly not entirely--spherical, and should synthesize a substantial mass of 56Ni. Our goal is to provide a qualitative assessment of the possibility that accretion of the progenitor envelope onto the black hole, which powers the LGRB, could also deposit sufficient energy and nickel mass in the envelope to produce a luminous supernova. For this, the energy dissipated near the black hole during accretion must be transported outward, where it can drive a supernova-like shockwave. Here we suggest that the energy is transported by convection and develop an analytical toy model, relying on global mass and energy conservation, for the dynamics of stellar collapse. The model suggests that a ~10,000 km/s shock can be driven into the envelope and that ~10^51 erg explosions are possible. The efficiency with which the accretion energy is being transferred to the envelope is governed by the competition of advection and convection at distances ~100-1,000 km from the black hole and is sensitive to the values of the convective mixing length, the magnitude of the effective viscous stress, and the specific angular momentum of the infalling envelope. Substantial masses of 56Ni may be synthesized in the convective accretion flow over the course of tens of seconds from the initial circularization of the infalling envelope around the black hole. The synthesized nickel is convectively mixed with a much larger mass of unburned ejecta.