The era of large transient surveys, gravitational-wave observatories, and multi-messenger astronomy has opened up new possibilities for our understanding of the evolution and final fate of massive ...stars. Most massive stars are born in binary or higher-order multiple systems and exchange mass with a companion star during their lives. In particular, the progenitors of a large fraction of compact-object mergers, and Galactic neutron stars (NSs) and black holes (BHs) have been stripped of their envelopes by a binary companion. Here, we study the evolution of single and stripped binary stars up to core collapse with the stellar evolution code M
ESA
and their final fates with a parametric supernova (SN) model. We find that stripped binary stars can have systematically different pre-SN structures compared to genuine single stars and thus also different SN outcomes. These differences are already established by the end of core helium burning and are preserved up to core collapse. Consequently, we find that Case A and B stripped stars and single and Case C stripped stars develop qualitatively similar pre-SN core structures. We find a non-monotonic pattern of NS and BH formation as a function of CO core mass that is different in single and stripped binary stars. In terms of initial mass, single stars of ≳35
M
⊙
all form BHs, while this transition is only at about 70
M
⊙
in stripped stars. On average, stripped stars give rise to lower NS and BH masses, higher explosion energies, higher kick velocities, and higher nickel yields. Within a simplified population-synthesis model, we show that our results lead to a significant reduction in the rates of BH–NS and BH–BH mergers with respect to typical assumptions made on NS and BH formation. Therefore, our models predict lower detection rates of such merger events with for example the advanced Laser Interferometer Gravitational-Wave Observatory (LIGO) than is often considered. Further, we show how certain features in the NS–BH mass distribution of single and stripped stars relate to the chirp-mass distribution of compact object mergers. Further implications of our findings are discussed with respect to the missing red-supergiant problem, a possible mass gap between NSs and BHs, X-ray binaries, and observationally inferred nickel masses from Type Ib/c and IIP SNe.
The presence of a nearby companion alters the evolution of massive stars in binary systems, leading to phenomena such as stellar mergers, x-ray binaries, and gamma-ray bursts. Unambiguous constraints ...on the fraction of massive stars affected by binary interaction were lacking. We simultaneously measured all relevant binary characteristics in a sample of Galactic massive O stars and quantified the frequency and nature of binary interactions. More than 70% of all massive stars will exchange mass with a companion, leading to a binary merger in one-third of the cases. These numbers greatly exceed previous estimates and imply that binary interaction dominates the evolution of massive stars, with implications for populations of massive stars and their supernovae.
Because the majority of massive stars are born as members of close binary systems, populations of massive main-sequence stars contain stellar mergers and products of binary mass transfer. We simulate ...populations of massive stars accounting for all major binary evolution effects based on the most recent binary parameter statistics and extensively evaluate the effect of model uncertainties. Assuming constant star formation, we find that 8 super(+9) sub(-4)% of a sample of early-type stars are the products of a merger resulting from a close binary system. In total we find that 30 super(+10) sub(-15)% of massive main-sequence stars are the products of binary interaction. We show that the commonly adopted approach to minimize the effects of binaries on an observed sample by excluding systems detected as binaries through radial velocity campaigns can be counterproductive. Systems with significant radial velocity variations are mostly pre-interaction systems. Excluding them substantially enhances the relative incidence of mergers and binary products in the non-radial velocity variable sample. This poses a challenge for testing single stellar evolutionary models. It also raises the question of whether certain peculiar classes of stars, such as magnetic O stars, are the result of binary interaction and it emphasizes the need to further study the effect of binarity on the diagnostics that are used to derive the fundamental properties (star-formation history, initial mass function, mass-to-light ratio) of stellar populations nearby and at high redshift.
The relative importance of primordial molecular cloud fragmentation versus large-scale accretion still remains to be assessed in the context of massive core/star formation. Studying the kinematics of ...the dense gas surrounding massive-star progenitors can tell us the extent to which large-scale flow of material impacts the growth in mass of star-forming cores. Here we present a comprehensive dataset of the 5500(±800) M⊙ infrared dark cloud SDC335.579-0.272 (hereafter SDC335), which exhibits a network of cold, dense, parsec-long filaments. Atacama Large Millimeter Array (ALMA) Cycle 0 observations reveal two massive star-forming cores, MM1 and MM2, sitting at the centre of SDC335 where the filaments intersect. With a gas mass of 545(-385+770) M⊙ contained within a source diameter of 0.05 pc, MM1 is one of the most massive, compact protostellar cores ever observed in the Galaxy. As a whole, SDC335 could potentially form an OB cluster similar to the Trapezium cluster in Orion. ALMA and Mopra single-dish observations of the SDC335 dense gas furthermore reveal that the kinematics of this hub-filament system are consistent with a global collapse of the cloud. These molecular-line data point towards an infall velocity Vinf = 0.7( ± 0.2) km s-1, and a total mass infall rate Ṁinf ≃ 2.5(±1.0) × 10-3 M⊙ yr-1 towards the central pc-size region of SDC335. This infall rate brings 750(±300) M⊙ of gas to the centre of the cloud per free-fall time (tff = 3 × 105 yr). This is enough to double the mass already present in the central pc-size region in 3.5-1.0+2.2 × tff. These values suggest that the global collapse of SDC335 over the past million year resulted in the formation of an early O-type star progenitor at the centre of the cloud’s gravitational potential well.
Context. Recent studies with, for example, Spitzer and Herschel have suggested that star formation in dense molecular gas may be governed by essentially the same “law” in Galactic clouds and external ...galaxies. This conclusion remains controversial, however, in large part because different tracers have been used to probe the mass of dense molecular gas in Galactic and extragalactic studies. Aims. We aimed to calibrate the HCN and HCO+ lines commonly used as dense gas tracers in extragalactic studies and to test the possible universality of the star-formation efficiency in dense gas (≳104 cm-3), SFEdense. Methods. We conducted wide-field mapping of the Aquila, Ophiuchus, and Orion B clouds at ~0.04 pc resolution in the J = 1 − 0 transition of HCN, HCO+, and their isotopomers. For each cloud, we derived a reference estimate of the dense gas mass MHerschelAV > 8, as well as the strength of the local far-ultraviolet (FUV) radiation field, using Herschel Gould Belt survey data products, and estimated the star-formation rate from direct counting of the number of Spitzer young stellar objects. Results. The H13CO+(1–0) and H13CN(1–0) lines were observed to be good tracers of the dense star-forming filaments detected with Herschel. Comparing the luminosities LHCN and LHCO+ measured in the HCN and HCO+ lines with the reference masses MHerschelAV > 8, the empirical conversion factors αHerschel − HCN (=MHerschelAV > 8/LHCN) and αHerschel − HCO+ (=MHerschelAV > 8/LHCO+) were found to be significantly anti-correlated with the local FUV strength. In agreement with a recent independent study of Orion B by Pety et al., the HCN and HCO+ lines were found to trace gas down to AV ≳ 2. As a result, published extragalactic HCN studies must be tracing all of the moderate density gas down to nH2 ≲ 103 cm-3. Estimating the contribution of this moderate density gas from the typical column density probability distribution functions in nearby clouds, we obtained the following G0-dependent HCN conversion factor for external galaxies: αHerschel − HCNfit′ = 64 × G0-0.34. Re-estimating the dense gas masses in external galaxies with αHerschel − HCNfit′(G0), we found that SFEdense is remarkably constant, with a scatter of less than 1.5 orders of magnitude around 4.5 × 10-8 yr-1, over eight orders of magnitude in dense gas mass. Conclusions. Our results confirm that SFEdense of galaxies is quasi-universal on a wide range of scales from ~ 1–10 pc to > 10 kpc. Based on the tight link between star formation and filamentary structure found in Herschel studies of nearby clouds, we argue that SFEdense is primarily set by the “microphysics” of core and star formation along filaments.
ABSTRACT
We investigate the impact of strong initial magnetic fields in core-collapse supernovae of non-rotating progenitors by simulating the collapse and explosion of a $16.9\, \mathrm{M}_\odot$ ...star for a strong- and weak-field case assuming a twisted-torus field with initial central field strengths of ${\approx }10^{12}$ and ${\approx }10^{6}\, \mathrm{G}$. The strong-field model has been set up with a view to the fossil-field scenario for magnetar formation and emulates a pre-collapse field configuration that may occur in massive stars formed by a merger. This model undergoes shock revival already $100\, \mathrm{ms}$ after bounce and reaches an explosion energy of $9.3\times 10^{50}\, \mathrm{erg}$ at $310\, \mathrm{ms}$, in contrast to a more delayed and less energetic explosion in the weak-field model. The strong magnetic fields help trigger a neutrino-driven explosion early on, which results in a rapid rise and saturation of the explosion energy. Dynamically, the strong initial field leads to a fast build-up of magnetic fields in the gain region to 40 per cent of kinetic equipartition and also creates sizable pre-shock ram pressure perturbations that are known to be conducive to asymmetric shock expansion. For the strong-field model, we find an extrapolated neutron star kick of ${\approx }350\, \mathrm{km}\, \mathrm{s}^{-1}$, a spin period of ${\approx }70\, \mathrm{ms}$, and no spin-kick alignment. The dipole field strength of the proto-neutron star is $2\times 10^{14}\, \mathrm{G}$ by the end of the simulation with a declining trend. Surprisingly, the surface dipole field in the weak-field model is stronger, which argues against a straightforward connection between pre-collapse fields and the birth magnetic fields of neutron stars.
The 30 Doradus star-forming region in the Large Magellanic Cloud is a nearby analog of large star-formation events in the distant universe. We determined the recent formation history and the initial ...mass function (IMF) of massive stars in 30 Doradus on the basis of spectroscopic observations of 247 stars more massive than 15 solar masses (Formula: see text). The main episode of massive star formation began about 8 million years (My) ago, and the star-formation rate seems to have declined in the last 1 My. The IMF is densely sampled up to 200 Formula: see text and contains 32 ± 12% more stars above 30 Formula: see text than predicted by a standard Salpeter IMF. In the mass range of 15 to 200 Formula: see text, the IMF power-law exponent is Formula: see text, shallower than the Salpeter value of 2.35.
Pacific low‐frequency variability (timescale > 8 years) exhibits a well‐known El Niño‐like pattern of basin‐scale sea surface temperature, which is found in all the major modes of Pacific decadal ...climate. Using a set of climate model experiments and observations, we decompose the mechanisms contributing to the growth, peak, and decay of the Pacific low‐frequency spatial variance. We find that the El Niño‐like interdecadal pattern is established through the combined actions of Pacific meridional modes (MM) and the El Niño–Southern Oscillation (ENSO). Specifically, in the growth phase of the pattern, subtropical stochastic excitation of the MM energizes the tropical low‐frequency variance acting as a red noise process. Once in the tropics, this low‐frequency variance is amplified by ocean‐atmospheric feedbacks as the pattern reaches its peak phase. At the same time, atmospheric teleconnections distribute the variance from the tropics to the extratropics, where the pattern ultimately decays. In this stochastic red noise model of Pacific climate, the timescale of the extra‐tropical/tropical interactions (1–2 years) permits the stochastic excitation of the ENSO‐like pattern of decadal and interdecadal variance.
Key Points
Pacific climate variability arises from the interaction of ENSO and Pacific Meridional Modes
Context. The earliest phases of the formation of high-mass stars are not well known. It is unclear whether high-mass cores in monolithic collapse exist or not, and what the accretion process and ...origin of the material feeding the precursors of high-mass stars are. As outflows are natural consequences of the accretion process, they represent one of the few (indirect) tracers of accretion. Aims. We aim to search for individual outflows from high-mass cores in Cygnus X and to study the characteristics of the detected ejections. We compare these to what has been found for the low-mass protostars, to understand how ejection and accretion change and behave with final stellar mass. Methods. We used CO (2–1) PdBI observations towards six massive dense clumps, containing a total of 9 high-mass cores. We estimated the bolometric luminosities and masses of the 9 high-mass cores and measured the energetics of outflows. We compared our sample to low-mass objects studied in the literature and developed simple evolutionary models to reproduce the observables. Results. We find that 8 out of 9 high-mass cores are driving clear individual outflows. They are therefore true equivalents of Class 0 protostars in the high-mass regime. The remaining core, CygX-N53 MM2, has only a tentative outflow detection. It could be one of the first examples of a true individual high-mass prestellar core. We also find that the momentum flux of high-mass objects has a linear relation to the reservoir of mass in the envelope, as a scale up of the relations previously found for low-mass protostars. This suggests a fundamental proportionality between accretion rates and envelope masses. The linear dependency implies that the timescale for accretion is similar for high- and low-mass stars. Conclusions. The existence of strong outflows driven by high-mass cores in Cygnus X clearly indicates that high-mass Class 0 protostars exist. The collapsing envelopes of these Class 0 objects have similar sizes and a similar fragmentation scale to the low-mass equivalents, and have enough mass to directly form high-mass stars from a monolithic collapse. If the pre-collapse evolution is quasi-static, the fragmentation scale is expected to limit the size of the initial mass reservoirs for all masses leading to higher densities at birth and therefore shorter free-fall times for higher mass stars. However, we find the collapse timescales to be similar for both low- and high-mass objects. This implies that in a quasi-static view, we would require significant turbulent/magnetic support to slow down the collapse of the more massive envelopes. But with this support still to be discovered, and based on independent indications of large dynamics in pre-collapse gas for high-mass star formation, we propose that such an identical collapse timescale implies that the initial densities, which should set the duration of the collapse, should be similar for all masses. Since the fragmentation scale is identical for all masses, a lower initial density requires that the mass that incorporates massive stars has to have been accreted from larger scales than those of low-mass stars and in a dynamical way.
Context.
Dense molecular filaments are ubiquituous in the interstellar medium, yet their internal physical conditions and the role of gravity, turbulence, the magnetic field, radiation, and the ...ambient cloud during their evolution remain debated.
Aims.
We study the kinematics and physical conditions in the Musca filament, the ambient cloud, and the Chamaeleon-Musca complex to constrain the physics of filament formation.
Methods.
We produced CO(2–1) isotopologue maps with the APEX telescope that cut through the Musca filament. We further study a NANTEN2
12
CO(1–0) map of the full Musca cloud, H
I
emission of the Chamaeleon-Musca complex, a
Planck
polarisation map, line radiative transfer models,
Gaia
data, and synthetic observations from filament formation simulations.
Results.
The Musca cloud, with a size of ~3–6 pc, contains multiple velocity components. Radiative transfer modelling of the CO emission indicates that the Musca filament consists of a cold (~10 K), dense (
n
H
2
∼ 10
4
cm
−3
) crest, which is best described with a cylindrical geometry. Connected to the crest, a separate gas component at
T
~ 15 K and
n
H
2
∼ 10
3
cm
−3
is found, the so-called strands. The velocity-coherent filament crest has an organised transverse velocity gradient that is linked to the kinematics of the nearby ambient cloud. This velocity gradient has an angle ≥30° with respect to the local magnetic field orientation derived from
Planck
, and the magnitude of the velocity gradient is similar to the transonic linewidth of the filament crest. Studying the large scale kinematics, we find coherence of the asymmetric kinematics from the 50 pc H
I
cloud down to the Musca filament. We also report a strong C
18
O/
13
CO abundance drop by an order of magnitude from the filament crest to the strands over a distance <0.2 pc in a weak ambient far-ultraviolet (FUV) field.
Conclusions.
The dense Musca filament crest is a long-lived (several crossing times), dynamic structure that can form stars in the near future because of continuous mass accretion replenishing the filament. This mass accretion on the filament appears to be triggered by a H
I
cloud–cloud collision, which bends the magnetic field around dense filaments. This bending of the magnetic field is then responsible for the observed asymmetric accretion scenario of the Musca filament, which is, for instance, seen as a V-shape in the position–velocity (PV) diagram.