Search

We present the continued analysis of polarization and Faraday rotation for the supernova remnants (SNRs) G46.8−0.3 and G39.2−0.3 in L -band (1–2 GHz) radio continuum in The H i /OH/Recombination line (THOR) survey. In this work, we present our investigation of Faraday depth fluctuations from angular scales comparable to the size of the SNRs down to scales less than our 16″ beam (≲0.7 pc) from Faraday dispersion (σ ϕ ). From THOR, we find median σ ϕ of 15.9 ± 3.2 rad m−2 for G46.8−0.3 and 17.6 ± 1.6 rad m−2 for G39.2−0.3. When comparing to polarization at 6 cm, we find evidence for σ ϕ ≳ 30 rad m−2 in localized regions where we detect no polarization in THOR. We combine Faraday depth dispersion with the rotation measure (RM) structure function (SF) and find evidence for a break in the SF on scales less than the THOR beam. We estimate the RM SF of the foreground interstellar medium using the SF of extragalactic radio sources and pulsars to find that the RM fluctuations we measure originate within the SNRs for all but the largest angular scales. [ABSTRACT FROM AUTHOR]

The semianalytical model a-sloth (Ancient Stars and Local Observables by Tracing Halos) is the first public code that connects the formation of the first stars and galaxies to observables. After several successful projects with this model, we publish the source code (https://gitlab.com/thartwig/asloth) and describe the public version in this paper. The model is based on dark matter merger trees that can either be generated based on Extended Pressâ€"Schechter theory or be imported from dark matter simulations. On top of these merger trees, a-sloth applies analytical recipes for baryonic physics to model the formation of both metal-free and metal-poor stars and the transition between them with unprecedented precision and fidelity. a-sloth samples individual stars and includes radiative, chemical, and mechanical feedback. It is calibrated based on six observables, such as the optical depth to Thomson scattering, the stellar mass of the Milky Way and its satellite galaxies, the number of extremely metal-poor stars, and the cosmic star formation rate density at high redshift. a-sloth has versatile applications with moderate computational requirements. It can be used to constrain the properties of the first stars and high- z galaxies based on local observables, predicts properties of the oldest and most metal-poor stars in the Milky Way, can serve as a subgrid model for larger cosmological simulations, and predicts next-generation observables of the early universe, such as supernova rates or gravitational wave events. [ABSTRACT FROM AUTHOR]

The lack of observations of abundance patterns originating in pair-instability supernovae has been a long-standing problem in relation to the first stars. This class of supernovae is expected to have an abundance pattern with a strong oddâ€"even effect, making it substantially different from present-day supernovae. In this study, we use a cosmological radiation hydrodynamics simulation to model such supernovae and the subsequent formation of the second generation of stars. We incorporate streaming velocities for the first time. There are 14 star-forming minihalos in our 1 cMpc h âˆ'1 box, leading to 14 supernovae occurring before redshift z = 19.5, where we start reducing the complexity of the simulation. Following the explosions, extremely metal-poor stars form in 10 halos via internal and external enrichment, which makes it the most common outcome. Only one halo does not recollapse during the simulations. This result is at variance with the current (lack of) observations of metal-poor stars with pair-instability supernova abundance patterns, suggesting that these very massive stars might be rare even in the early universe. The results from this simulation also give us insights into what drives different modes of recollapse and what determines the mixing behavior of metals after very energetic supernovae. [ABSTRACT FROM AUTHOR]

A unique filament is identified in the Herschel maps of the Orion A giant molecular cloud. The filament, which we name the Stick, is ruler-straight and at an early evolutionary stage. Transverse position–velocity diagrams show two velocity components closing in on the Stick. The filament shows consecutive rings/forks in C18O (1−0) channel maps, which is reminiscent of structures generated by magnetic reconnection. We propose that the Stick formed via collision-induced magnetic reconnection (CMR). We use the magnetohydrodynamics code Athena++ to simulate the collision between two diffuse molecular clumps, each carrying an antiparallel magnetic field. The clump collision produces a narrow, straight, dense filament with a factor of >200 increase in density. The production of the dense gas is seven times faster than freefall collapse. The dense filament shows ring/fork-like structures in radiative transfer maps. Cores in the filament are confined by surface magnetic pressure. CMR can be an important dense-gas-producing mechanism in the Galaxy and beyond. [ABSTRACT FROM AUTHOR]

Understanding the formation of stellar clusters requires following the interplay between gas and newly formed stars accurately. We therefore couple the magnetohydrodynamics code FLASH to the N-body code ph4 and the stellar evolution code SeBa using the Astrophysical Multipurpose Software Environment (AMUSE) to model stellar dynamics, evolution, and collisional N-body dynamics and the formation of binary and higher-order multiple systems, while implementing stellar feedback in the form of radiation, stellar winds, and supernovae in FLASH. We here describe the algorithms used for each of these processes. We denote this integrated package Torch. We then use this novel numerical method to simulate the formation and early evolution of several examples of open clusters of ∼1000 stars formed from clouds with a mass range of 103M⊙ to 105M⊙. Analyzing the effects of stellar feedback on the gas and stars of the natal clusters, we find that in these examples, the stellar clusters are resilient to disruption, even in the presence of intense feedback. This can even slightly increase the amount of dense, Jeans unstable gas by sweeping up shells; thus, a stellar wind strong enough to trap its own H ii region shows modest triggering of star formation. Our clusters are born moderately mass segregated, an effect enhanced by feedback, and retained after the ejection of their natal gas, in agreement with observations. [ABSTRACT FROM AUTHOR]

We present the first complete, velocity-resolved [C ii ] 158 μm image of the M51 grand-design spiral galaxy, observed with the upgraded German Receiver for Astronomy at Terahertz frequencies instrument on the Stratospheric Observatory for Infrared Astronomy. [C ii ] is an important tracer of various phases of the interstellar medium (ISM), including ionized gas, neutral atomic, and diffuse molecular regions. We combine the [C ii ] data with H i , CO, 24 μm dust continuum, FUV, and NIR K-band observations to study the evolution of the ISM across M51's spiral arms in both position–position and position–velocity space. Our data show strong velocity gradients in H i , 12CO, and [C ii ] at the locations of the stellar arms (traced by K-band data) with a clear offset in position–velocity space between upstream molecular gas (traced by 12CO) and downstream star formation (traced by [C ii ]). We compare the observed position–velocity maps across the spiral arms with synthetic observations from numerical simulations of galaxies with both dynamical and quasi-stationary steady spiral arms that predict both tangential and radial velocities at the location of the spiral arms. We find that our observations, based on the observed velocity gradients and associated offset between CO and [C ii ], are consistent with the presence of shocks in the spiral arms in the inner parts of M51 and in the arm connecting the companion galaxy, M51b, in the outer parts of M51. [ABSTRACT FROM AUTHOR]

The time evolution of the probability density function (PDF) of the mass density is formulated and solved for systems in free-fall using a simple approximate function for the collapse of a sphere. We demonstrate that a pressure-free collapse results in a power-law tail on the high-density side of the PDF. The slope quickly asymptotes to the functional form PV(ρ)∝ρ–1.54 for the (volume-weighted) PDF and PM(ρ)∝ρ–0.54 for the corresponding mass-weighted distribution. From the simple approximation of the PDF we derive analytic descriptions for mass accretion, finding that dynamically quiet systems with narrow density PDFs lead to retarded star formation and low star formation rates (SFRs). Conversely, strong turbulent motions that broaden the PDF accelerate the collapse causing a bursting mode of star formation. Finally, we compare our theoretical work with observations. The measured SFRs are consistent with our model during the early phases of the collapse. Comparison of observed column density PDFs with those derived from our model suggests that observed star-forming cores are roughly in free-fall. [ABSTRACT FROM AUTHOR]

Observations of atomic or molecular lines can provide important information about the physical state of star-forming regions. In order to investigate the line profiles from dynamical collapsing massive star-forming regions (MSFRs), we model the emission from hydrodynamic simulations of a collapsing cloud in the absence of outflows. By performing radiative transfer calculations, we compute the optically thick HCO+ and optically thin N2H+ line profiles from two collapsing regions at different epochs. Due to large-scale collapse, theMSFRs have large velocity gradients, reaching up to 20 km s-1 pc-1 across the central core. The optically thin lines typically contain multiple velocity components resulting from the superposition of numerous density peaks along the line of sight. The optically thick lines are only marginally shifted to the blue side of the optically thin line profiles, and frequently do not have a central depression in their profiles due to self-absorption. As the regions evolve, the lines become brighter and the optically thick lines become broader. The lower-order HCO+ (1-0) transitions are better indicators of collapse than the higher-order (4-3) transitions. We also investigate how the beam sizes affect profile shapes. Smaller beams lead to brighter and narrower lines that are more skewed to the blue in HCO+ relative to the true core velocity, but show multiple components in N2H+. High-resolution observations (e.g., with Atacama Large Millimeter Array) can test these predictions and provide insights into the nature of MSFRs. [ABSTRACT FROM AUTHOR]

The first stars to form in the universe are believed to have distribution of masses biased toward massive stars. This contrasts with the present-day initial mass function, which has a predominance of stars with masses lower than 1M☉. Therefore, the mode of star formation must have changed as the universe evolved. Such a transition is attributed to a more efficient cooling provided by increasing metallicity. Especially dust cooling can overcome the compressional heating, which lowers the gas temperature thus increasing its instability to fragmentation. The purpose of this paper is to verify if dust cooling can efficiently cool the gas, and enhance the fragmentation of gas clouds at the early stages of the universe. To confirm that, we calculate a set of hydrodynamic simulations that include sink particles, which represent contracting protostars. The thermal evolution of the gas during the collapse is followed by making use of a primordial chemical network and also a recipe for dust cooling. We model four clouds with different amounts of metals (10-4, 10-5, 10-6 Z☉, and 0), and analyze how this property affect the fragmentation of star-forming clouds. We find evidence for fragmentation in all four cases, and hence conclude that there is no critical metallicity below which fragmentation is impossible. Nevertheless, there is a clear change in the behavior of the clouds at Z ≲ 10-5 Z☉, caused by the fact that at this metallicity, fragmentation takes longer to occur than accretion, leading to a flat mass function at lower metallicities. [ABSTRACT FROM AUTHOR]

Observations are revealing the ubiquity of filamentary structures in molecular clouds. As cores are often embedded in filaments, it is important to understand how line profiles from such systems differ from those of isolated cores. We perform radiative transfer calculations on a hydrodynamic simulation of a molecular cloud in order to model line emission from collapsing cores embedded in filaments. We model two optically thick lines, CS(2-1) and HCN(1-0), and one optically thin line, N2H+ (1-0), from three embedded cores. In the hydrodynamic simulation, gas self-gravity, turbulence, and bulk flows create filamentary regions within which cores form. Though the filaments have large dispersions, the N2H+ (1-0) lines indicate subsonic velocities within the cores. We find that the observed optically thick line profiles of CS(2-1) and HCN(1-0) vary drastically with viewing angle. In over 50% of viewing angles, there is no sign of a blue asymmetry, an idealized signature of infall motions in an isolated spherical collapsing core. Profiles that primarily trace the cores, with little contribution from the surrounding filament, are characterized by a systematically higher HCN(1-0) peak intensity. The N2H+ (1-0) lines do not follow this trend. We demonstrate that red asymmetric profiles are also feasible in the optically thick lines, due to emission from the filament or one-sided accretion flows onto the core. We conclude that embedded cores may frequently undergo collapse without showing a blue asymmetric profile, and that observational surveys including filamentary regions may underestimate the number of collapsing cores if based solely on profile shapes of optically thick lines. [ABSTRACT FROM AUTHOR]

Turbulence plays a major role in the formation and evolution of molecular clouds. Observationally, turbulent velocities are convolved with the density of an observed region. To correct for this convolution, we investigate the relation between the turbulence spectrum of model clouds, and the statistics of their synthetic observations obtained from principal component analysis (PCA). We apply PCA to spectral maps generated from simulated density and velocity fields, obtained from hydrodynamic simulations of supersonic turbulence, and from fractional Brownian motion (fBm) fields with varying velocity, density spectra, and density dispersion. We examine the dependence of the slope of the PCA pseudo-structure function, αPCA on intermittency, on the turbulence velocity (βv) and density (βn) spectral indexes, and on density dispersion. We find that PCA is insensitive to β n and to the log-density dispersion σ s, provided σs ≤ 2. For σs > 2, αPCA increases with σs due to the intermittent sampling of the velocity field by the density field. The PCA calibration also depends on intermittency. We derive a PCA calibration based on fBm structures with σs ≤ 2 and apply it to 367 13CO spectral maps of molecular clouds in the Galactic Ring Survey. The average slope of the PCA structure function, <αPCA> = 0.62 ± 0.2, is consistent with the hydrodynamic simulations and leads to a turbulence velocity exponent of (βv = 2.06 ± 0.6 for a non-intermittent, low density dispersion flow. Accounting for intermittency and density dispersion, the coincidence between the PCA slope of the GRS clouds and the hydrodynamic simulations suggests βv ≃ 1.9, consistent with both Burgers and compressible intermittent turbulence. [ABSTRACT FROM AUTHOR]

The cosmic dark ages ended a few hundred million years after the big bang, when the first stars began to fill the universe with new light. It has generally been argued that these stars formed in isolation and were extremely massive--perhaps 100 times as massive as the Sun. In a recent study, Clark and collaborators showed that this picture requires revision. They demonstrated that the accretion disks that build up around Population III stars are strongly susceptible to fragmentation and that the first stars should therefore form in clusters rather than in isolation. We here use a series of high-resolution hydrodynamical simulations performed with the moving mesh code AREPO to follow up on this proposal and to study the influence of environmental parameters on the level of fragmentation. We model the collapse of five independent minihalos from cosmological initial conditions, through the runaway condensation of their central gas clouds, to the formation of the first protostar, and beyond for a further 1000 years. During this latter accretion phase, we represent the optically thick regions of protostars by sink particles. Gas accumulates rapidly in the circumstellar disk around the first protostar, fragmenting vigorously to produce a small group of protostars. After an initial burst, gravitational instability recurs periodically, forming additional protostars with masses ranging from ~0.1 to 10 M⊙. Although the shape, multiplicity, and normalization of the protostellar mass function depend on the details of the sink-particle algorithm, fragmentation into protostars with diverse masses occurs in all cases, confirming earlier reports of Population III stars forming in clusters. Depending on the efficiency of later accretion and merging, Population III stars may enter the main sequence in clusters and with much more diverse masses than are commonly assumed. [ABSTRACT FROM AUTHOR]

It has recently been demonstrated that coherent relative streaming velocities of order 30km s-1 between dark matter and gas permeated the universe on scales below a few Mpc directly after recombination. Here we use a series of high-resolution moving-mesh calculations to show that these supersonic motions significantly influence the virialization of the gas in minihalos and delay the formation of the first stars. As the gas streams into minihalos with bulk velocities around 1 km s-1 at z ~ 20, the additional momentum and energy input reduces the gas fractions and central densities of the halos, increasing the typical virial mass required for efficient cooling by a factor of three and delaying Population III star formation by Δz ≃ 4. Since the distribution of the magnitude of the streaming velocities is narrowly peaked around a non-negligible value, this effect is important in most regions of the universe. As a consequence, the increased minimum halo mass implies a reduction of the absolute number of minihalos that can be expected to cool and form Population III stars by up to an order of magnitude. We further find that the streaming velocities increase the turbulent velocity dispersion of the minihalo gas, which could affect its ability to fragment and hence alter the mass function of the first stars. [ABSTRACT FROM AUTHOR]

We present a new method for estimating the H2 cooling rate in the optically thick regime in simulations of primordial star formation. Our new approach is based on the TreeCol algorithm, which projects matter distributions onto a spherical grid to create maps of column densities for each fluid element in the computational domain. We have improved this algorithm by using the relative gas velocities to weight the individual matter contributions with the relative spectral line overlaps, in order to properly account for the Doppler effect. We compare our new method to the widely used Sobolev approximation, which yields an estimate for the column density based on the local velocity gradient and the thermal velocity. This approach generally underestimates the photon escape probability because it neglects the density gradient and the actual shape of the cloud. We present a correction factor for the true line overlap in the Sobolev approximation and a new method based on local quantities, which fits the exact results reasonably well during the collapse of the cloud, with the error in the cooling rates always being less than 10%. Analytical fitting formulae fail at determining the photon escape probability after formation of the first protostar (error of ∼40%) because they are based on the assumption of spherical symmetry and therefore break down once a protostellar accretion disk has formed. Our method yields lower temperatures and hence promotes fragmentation for densities above ∼1010 cm–3 at a distance of ∼200 AU from the first protostar. Since the overall accretion rates are hardly affected by the cooling implementation, we expect Pop III stars to have lower masses in our simulations, compared to the results of previous simulations that used the Sobolev approximation. [ABSTRACT FROM AUTHOR]

Powerful jets and outflows are launched from the protostellar disks around newborn stars. These outflows carry enough mass and momentum to transform the structure of their parent molecular cloud and to potentially control star formation itself. Despite their importance, we have not been able to fully quantify the impact of jets and outflows during the formation of a star cluster. The main problem lies in limited computing power. We would have to resolve the magnetic jet-launching mechanism close to the protostar and at the same time follow the evolution of a parsec-size cloud for a million years. Current computer power and codes fall orders of magnitude short of achieving this. In order to overcome this problem, we implement a subgrid-scale (SGS) model for launching jets and outflows, which demonstrably converges and reproduces the mass, linear and angular momentum transfer, and the speed of real jets, with ∼1000 times lower resolution than would be required without the SGS model. We apply the new SGS model to turbulent, magnetized star cluster formation and show that jets and outflows (1) eject about one-fourth of their parent molecular clump in high-speed jets, quickly reaching distances of more than a parsec, (2) reduce the star formation rate by about a factor of two, and (3) lead to the formation of ∼1.5 times as many stars compared to the no-outflow case. Most importantly, we find that jets and outflows reduce the average star mass by a factor of ∼ three and may thus be essential for understanding the characteristic mass of the stellar initial mass function. [ABSTRACT FROM AUTHOR]

The formation of high-mass stars is usually accompanied by powerful protostellar outflows. Such high-mass outflows are not simply scaled-up versions of their lower-mass counterparts, since observations suggest that the collimation degree degrades with stellar mass. Theoretically, the origins of massive outflows remain open to question because radiative feedback and fragmentation of the accretion flow around the most massive stars, with M > 15 M☼, may impede the driving of magnetic disk winds. We here present a three-dimensional simulation of the early stages of core fragmentation and massive star formation that includes a subgrid-scale model for protostellar outflows. We find that stars that form in a common accretion flow tend to have aligned outflow axes, so that the individual jets of multiple stars can combine to form a collective outflow. We compare our simulation to observations with synthetic H2 and CO observations and find that the morphology and kinematics of such a collective outflow resembles some observed massive outflows, such as Cepheus A and DR 21. We finally compare physical quantities derived from simulated observations of our models to the actual values in the models to examine the reliability of standard methods for deriving physical quantities, demonstrating that those methods indeed recover the actual values to within a factor of two to three. [ABSTRACT FROM AUTHOR]

We derive the CO-to-H2 conversion factor, XCO = N(H2)/ICO, across the Perseus molecular cloud on sub-parsec scales by combining the dust-based N(H2) data with the ICO data from the COMPLETE Survey. We estimate an average XCO ∼ 3 × 1019 cm–2 K–1 km–1 s and find a factor of ∼3 variations in XCO between the five sub-regions in Perseus. Within the individual regions, XCO varies by a factor of ∼100, suggesting that XCO strongly depends on local conditions in the interstellar medium. We find that XCO sharply decreases at AV ≲ 3 mag but gradually increases at AV ≳ 3 mag, with the transition occurring at AV where ICO becomes optically thick. We compare the N(H I), N(H2), ICO, and XCO distributions with two models of the formation of molecular gas, a one-dimensional photodissociation region (PDR) model and a three-dimensional magnetohydrodynamic (MHD) model, tracking both the dynamical and chemical evolution of gas. The PDR model based on the steady state and equilibrium chemistry reproduces our data very well but requires a diffuse halo to match the observed N(H I) and ICO distributions. The MHD model matches our data reasonably well, suggesting that time-dependent effects on H2 and CO formation are insignificant for an evolved molecular cloud like Perseus. However, we find interesting discrepancies, including a broader range of N(H I), likely underestimated ICO, and a large scatter of ICO at small AV. These discrepancies most likely result from strong compressions and rarefactions and density fluctuations in the MHD model. [ABSTRACT FROM AUTHOR]

The role of turbulence and magnetic fields is studied for star formation in molecular clouds. We derive and compare six theoretical models for the star formation rate (SFR)—the Krumholz & McKee (KM), Padoan & Nordlund (PN), and Hennebelle & Chabrier (HC) models, and three multi-freefall versions of these, suggested by HC—all based on integrals over the log-normal distribution of turbulent gas. We extend all theories to include magnetic fields and show that the SFR depends on four basic parameters: (1) virial parameter αvir; (2) sonic Mach number ; (3) turbulent forcing parameter b, which is a measure for the fraction of energy driven in compressive modes; and (4) plasma with the Alfvén Mach number . We compare all six theories with MHD simulations, covering cloud masses of 300 to 4 × 106M☼ and Mach numbers -50 and -∞, with solenoidal (b = 1/3), mixed (b = 0.4), and compressive turbulent (b = 1) forcings. We find that the SFR increases by a factor of four between and 50 for compressive turbulent forcing and αvir ∼ 1. Comparing forcing parameters, we see that the SFR is more than 10 times higher with compressive than solenoidal forcing for simulations. The SFR and fragmentation are both reduced by a factor of two in strongly magnetized, trans-Alfvénic turbulence compared to hydrodynamic turbulence. All simulations are fit simultaneously by the multi-freefall KM and multi-freefall PN theories within a factor of two over two orders of magnitude in SFR. The simulated SFRs cover the range and correlation of SFR column density with gas column density observed in Galactic clouds, and agree well for star formation efficiencies SFE = 1%-10% and local efficiencies ε = 0.3-0.7 due to feedback. We conclude that the SFR is primarily controlled by interstellar turbulence, with a secondary effect coming from magnetic fields. [ABSTRACT FROM AUTHOR]

We present the first three-dimensional simulations to include the effects of dark matter annihilation feedback during the collapse of primordial minihalos. We begin our simulations from cosmological initial conditions and account for dark matter annihilation in our treatment of the chemical and thermal evolution of the gas. The dark matter is modeled using an analytical density profile that responds to changes in the peak gas density. We find that the gas can collapse to high densities despite the additional energy input from the dark matter. No objects supported purely by dark matter annihilation heating are formed in our simulations. However, we find that dark matter annihilation heating has a large effect on the evolution of the gas following the formation of the first protostar. Previous simulations without dark matter annihilation found that protostellar disks around Population III stars rapidly fragmented, forming multiple protostars that underwent mergers or ejections. When dark matter annihilation is included, however, these disks become stable to radii of 1000 AU or more. In the cases where fragmentation does occur, it is a wide binary that is formed. [ABSTRACT FROM AUTHOR]

The formation of massive stars exceeding 10 M☼ usually results in large-scale molecular outflows. Numerical simulations, including ionization, of the formation of such stars show evidence for ionization-driven molecular outflows. Here we examine whether the outflows seen in these models reproduce the observations. We compute synthetic ALMA and CARMA maps of CO emission lines of the outflows, and compare their signatures to existing single-dish and interferometric data. We find that the ionization-driven models can only reproduce weak outflows around high-mass star-forming regions. We argue that expanding H II regions probably do not represent the dominant mechanism for driving observed outflows. We suggest instead that observed outflows are driven by the collective action of the outflows from the many lower-mass stars that inevitably form around young massive stars in a cluster. [ABSTRACT FROM AUTHOR]

We study the amplification of magnetic fields during the formation of primordial halos. The turbulence generated by gravitational infall motions during the formation of the first stars and galaxies can amplify magnetic fields very efficiently and on short timescales up to dynamically significant values. Using the Kazantsev theory, which describes the so-called small-scale dynamo—a magnetohydrodynamical process converting kinetic energy from turbulence into magnetic energy—we can then calculate the growth rate of the small-scale magnetic field. Our calculations are based on a detailed chemical network and we include non-ideal magnetohydrodynamical effects such as ambipolar diffusion and Ohmic dissipation. We follow the evolution of the magnetic field up to larger scales until saturation occurs on the Jeans scale. Assuming a weak magnetic seed field generated by the Biermann battery process, both Burgers and Kolmogorov turbulence lead to saturation within a rather small density range. Such fields are likely to become relevant after the formation of a protostellar disk and, thus, could influence the formation of the first stars and galaxies in the universe. [ABSTRACT FROM AUTHOR]

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 BH 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 BH 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 an 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. [ABSTRACT FROM AUTHOR]