Binary Evolution Research Articles

Triple trouble with PSR J1618−3921: Mass measurements and orbital dynamics of an eccentric millisecond pulsar

Context. PSR J1618−3921 is one of five known millisecond pulsars (MSPs) in eccentric orbits (eMPSs) located in the Galactic plane, whose formation is poorly understood. Earlier studies of these objects revealed significant discrepancies between observations and predictions from standard binary evolution scenarios of pulsar-helium white dwarf (HeWD) binaries, especially in the case of PSR J0955−6150, for which mass measurements ruled out most eMSP formation models. Aims. We aim to measure the masses of the pulsar and its companion, and constrain the orbital configuration of PSR J1618−3921. This facilitates understanding similarities among eMSPs and could offer hints as to their formation mechanism. Methods. We conducted observations with the L-band receiver of the MeerKAT radio telescope and the UWL receiver of the Parkes Murriyang radio telescope between 2019 and 2021. These data were added to archival Parkes and Nançay observations. We performed a full analysis on this joint data set with a timing baseline of 23 years. We also used the data from recent observations to give a brief account of the emission properties of J1618−3921, including a rotating vector model (RVM) fit of the linear polarisation position angle of the pulsar. Results. From the timing analysis, we measure a small but significant proper motion of the pulsar. The long timing baseline allowed for a highly significant measurement of the rate of advance of periastron of ω̇ = (0.00145±0.00010)°yr−1. Despite the tenfold improvement in timing precision from MeerKAT observations, we can only report a low-significance detection of the orthometric Shapiro delay parameters, h3 = 2.70−1.47+2.07 μs and ς = 0.68−0.09+0.13. Under the assumption of the validity of general relativity (GR), the self-consistent combination of these three parameters leads to mass estimates of the total and individual masses in the binary of Mtot = 1.42−0.19+0.20 M⊙, Mc = 0.20−0.03+0.11 M⊙, and Mp = 1.20−0.20+0.19 M⊙. We detect an unexpected change in the orbital period of Ṗb = −2.26−0.33+0.35 × 10−12, that is an order of magnitude larger and carries an opposite sign to what is expected from the Galactic acceleration and the Shklovskii effect, which are a priori the only non-negligible contributions expected for Ṗb. We also detect a significant second derivative of the spin frequency, f̈. The RVM fit reveals a viewing angle of ζ = (111 ± 1)°. Furthermore, we report an unexpected, abrupt change in the mean pulse profile in June 2021 of unknown origin. Conclusions. We propose that the anomalous Ṗb and f̈ we measure for J1618−3921 indicate an additional varying acceleration due to a nearby mass. The J1618−3921 binary system is likely part of a hierarchical triple, but with the third component much farther away than the outer component of the MSP in a triple star system, PSR J0337+1715. This finding suggests that at least some eMSPs might have formed in triple star systems. Although the uncertainties are large, the binary companion mass is consistent with the Pb − MWD relation, which has been verified for circular HeWD binaries and also for the two HeWDs in the PSR J0337+1715 system. Future regular observations with the MeerKAT telescope will, due to the further extension of the timing baseline, improve the measurement of Ṗb and f̈. This will help us further understand the nature of this system, and perhaps improve our understanding of eMSPs in general.

Binary Yellow Supergiants in the Magellanic Clouds. I. Photometric Candidate Identification

Recent works have constrained the binary fraction of evolved populations of massive stars in local galaxies such as red supergiants and Wolf–Rayet stars, but the binary fraction of yellow supergiants (YSGs) in the Hertzsprung gap remains unconstrained. Binary evolution theory predicts that the Hertzsprung gap is home to multiple populations of binary systems with varied evolutionary histories. In this paper, we develop a method to distinguish single YSGs from YSG plus O- or B-type main-sequence binaries using optical and ultraviolet photometry, and then apply this method to identify candidate YSG binaries in the Magellanic Clouds. After constructing a set of combined stellar atmosphere models, we find that optical photometry is, given typical measurement and reddening uncertainties, sufficient to discern single YSGs from YSG+OB binaries if the OB-star is at least ∼5M ⊙ for T eff,YSG ∼ 4000 K, but requires a ∼20M ⊙ OB star for YSGs up to T eff,YSG ∼ 9000 K. For these hotter YSG temperatures, ultraviolet photometry allows binaries with OB companions as small as ∼7M ⊙ to be identified. We use color–color spaces developed from these models to search for evidence of excess blue or ultraviolet light in a set of ∼1000 YSG candidates in the Magellanic Clouds. We identify hundreds of candidate YSG binary systems and report a preliminary fraction of YSGs that show a blue/UV color excess of 20%–60%. Spectroscopic follow-up is now required to confirm the true nature of this population.

An upper limit on the spins of merging binary black holes formed through isolated binary evolution

As the sensitivity of ground-based gravitational wave detectors progressively increases, observations of black hole mergers will provide us with the joint distribution of their masses and spins. This will be a critical benchmark to validate different formation scenarios. Merging binary black holes formed through the evolution of isolated binary systems require both components to be stripped of their hydrogen envelopes before core-collapse. The rotation rates of such stripped stars are constrained by the critical rotation limit at their surface, including its deviation from the Keplerian value owing to the outward force provided by radiation. This sets a restriction on their angular momentum content at core-collapse. We aim to determine if this restriction plays a role in the spins of binary black hole mergers. We used detailed calculations of stripped stars with the MESA code at low metallicities ($Z=Z_ $Z_ and $Z_ to determine the dimensionless spins of black holes produced by critically rotating stellar progenitors. To study how such progenitors can arise, we considered their formation through chemically homogeneous evolution (CHE) in binary stars. We used a semi-analytical model to study the physical processes that determine the final angular momentum of CHE binaries, and compared our results against available population synthesis models that rely on detailed binary evolution calculations. We find that above black hole masses of $ 25M_ the dimensionless spin parameter of critically rotating stripped stars ($a=Jc/(GM^2)$) is below unity. This results in an exclusion region at high chirp masses and effective spins that cannot be populated by isolated binary evolution. CHE can produce binaries where both black holes hit this limit, producing a pileup at the boundary of the excluded region. High-spin black holes arise from very low-metallicity CHE systems with short delay times, which merge at higher redshifts. On the other hand, the contribution of CHE to merging binary black holes detected in the third observing run of the LVK collaboration is expected to be dominated by systems with low spins ($ eff <0.5$) that merge near redshift zero. Owing to its higher projected sensitivity and runtime, the fourth observing run of the LVK collaboration can potentially place constraints on the high-spin population and the existence of a limit set by critical rotation.

Tidal Spin-up of Subdwarf B Stars

Hot subdwarf B (sdB) stars are stripped helium-burning stars that are often found in close binaries, where they experience strong tidal interactions. The dissipation of tidally excited gravity waves alters their rotational evolution throughout the sdB lifetime. While many sdB binaries have well-measured rotational and orbital frequencies, there have been few theoretical efforts to accurately calculate the tidal torque produced by gravity waves. In this work, we directly calculate the tidal excitation of internal gravity waves in realistic sdB stellar models and integrate the coupled spin–orbit evolution of sdB binaries. We find that for canonical sdB (M sdB = 0.47 M ⊙) binaries, the transitional orbital period below which they could reach tidal synchronization in the sdB lifetime is ∼0.2 day, with weak dependence on the companion masses. For low-mass sdBs (M sdB = 0.37 M ⊙) formed from more massive progenitor stars, the transitional orbital period becomes ∼0.15 day. These values are very similar to the tidal synchronization boundary (∼0.2 day) evident from observations. We discuss the dependence of tidal torques on stellar radii, and we make predictions for the rapidly rotating white dwarfs formed from synchronized sdB binaries.

Uncertainty of the white dwarf astrophysical gravitational wave background

The astrophysical gravitational wave background (AGWB) is a stochastic gravitational wave (GW) signal emitted by different populations of in-spiralling binary systems containing compact objects throughout the Universe . In the frequency range between 10$^ $ and 10$^ $ hertz (Hz), it will be detected by future space-based gravitational wave detectors, such as Laser Interferometer Space Antenna (LISA). In a recent work, we concluded that the white dwarf (WD) contribution to the AGWB dominates that of black holes (BHs) and neutron stars (NSs) We aim to investigate the uncertainties of the WD AGWB that arise from the use of different stellar metallicities, star formation rate density (SFRD) models, and binary evolution models. We used the code we previously developed to determine the WD component of the AGWB. We used a metallicity-dependent SFRD based on an earlier work to construct five different SFRD models. We used four different population models based on a range of common-envelope treatments and six different metallicities for each model. For all possible combinations, the WD component of the AGWB is dominant over other populations of compact objects. The effects of metallicity and population model are less significant than the effect of a (metallicity dependent) SFRD model. We find a range of about a factor of 5 in the level of the WD AGWB around a value of $ WD $ at 1 mHz and a shape that is weakly dependent on the model. We find the uncertainty for the WD component of the AGWB to be about a factor of 5. We note that there are other uncertainties that have an effect on this signal as well. We discuss whether the turnover of the WD AGWB at 10 mHz will be detectable by LISA and find it to be likely. We confirm our previous findings asserting that the WD component of the AGWB dominates over other populations, in particular, BHs.

Stable case BB/BC mass transfer to form GW190425-like massive binary neutron star mergers

On April 25, 2019, the LIGO-Virgo Collaboration discovered a gravitational-wave (GW) signal from a binary neutron star (BNS) merger, that is, GW190425. Due to the inferred large total mass, the origin of GW190425 remains unclear. Assuming GW190425 originated from the standard isolated binary evolution channel, its immediate progenitor is considered to be a close binary system, consisting of a He-rich star and a NS just after the common envelope phase. We aim to study the formation of GW190425 in a solar-like environment by using the detailed binary evolution code MESA We perform detailed stellar structure and binary evolution calculations that take into account mass loss, internal differential rotation, and tidal interactions between a He-rich star and a NS companion. We explore the parameter space of the initial binary properties, including initial NS and He-rich masses and initial orbital period. We find that the immediate post-common-envelope progenitor system, consisting of a primary $ NS and a secondary He-rich star with an initial mass of $ in a close binary with an initial period of $ days days $), that experiences stable Case BB/BC mass transfer (MT) during binary evolution, can reproduce the formation of GW190425-like BNS events. Our studies reveal that the secondary He-rich star of the GW190425's progenitor before its core collapse can be efficiently spun up through tidal interaction, finally remaining as a NS with rotational energy even reaching $ erg $, which is always much higher than the neutrino-driven energy of the supernova (SN) explosion. If the newborn secondary NS is a magnetar, we expect that GW190425 can be the remnant of a magnetar-driven SN, namely a magnetar-driven ultra-stripped SN, a superluminous SN, or a broad-line Type Ic SN. Our results show that GW190425 could be formed through the isolated binary evolution, which involves a stable Case BB/BC MT just after the common envelope phase. On top of that, we show the He-rich star can be tidally spun up, potentially forming a spinning magnetized NS (magnetar) during the second SN explosion.

Accretion mediated spin-eccentricity correlations in LISA massive black hole binaries

Abstract We examine expected effective spin (χeff, 1yr) and orbital eccentricity ($e_{1\rm yr}$) correlations for a population of observable equal-mass massive black hole binaries (MBHBs) with total redshifted mass Mz ∼ [104.5, 107.5] M⊙ embedded in a circumbinary disc (CBD) at redshifts z = 1 and z = 2, one-year before merging in the LISA band. We find a strong correlation between measurable eccentricity and negative effective spin for MBHBs that are carried to merger by retrograde accretion. This is due to the well-established eccentricity pumping of retrograde accretion and less-well-established formation of retrograde minidiscs coupled with a stable retrograde CBD throughout the binary evolution from the self-gravitating radius. Conversely, prograde accretion channels result in positive $\chi _{{\rm eff},1\rm yr}$ and non-measurable $e_{1\rm yr}$ except for almost unity Eddington ratio and Mz ≲ 105 M⊙ MBHBs at z = 1. This clear contrast between the two CBD orientations – and particularly the unique signature of retrograde configurations – provides a promising way to unlock the mysteries of MBHB formation channels in the LISA era.

From spherical stars to disk-like structures: 3D common-envelope evolution of massive binaries beyond inspiral

From spherical stars to disk-like structures: 3D common-envelope evolution of massive binaries beyond inspiral

Advanced classification of hot subdwarf binaries using artificial intelligence techniques and Gaia DR3 data

Hot subdwarf stars are compact blue evolved objects, burning helium in their cores surrounded by a tiny hydrogen envelope. In the Hertzsprung-Russell Diagram they are located by the blue end of the Horizontal Branch. Most models agree on a quite probable common envelope binary evolution scenario in the Red Giant phase. However, the current binarity rate for these objects is yet unsolved, but key, question in this field. This study aims to develop a novel classification method for identifying hot subdwarf binaries within large datasets using Artificial Intelligence techniques and data from the third Gaia data release (GDR3). The results will be compared with those obtained previously using Virtual Observatory techniques on coincident samples. The methods used for hot subdwarf binary classification include supervised and unsupervised machine learning techniques. Specifically, we have used Support Vector Machines (SVM) to classify 3084 hot subdwarf stars based on their colour-magnitude properties. Among these, 2815 objects have Gaia DR3 BP/RP spectra, which were classified using Self-Organizing Maps (SOM) and Convolutional Neural Networks (CNN). In order to ensure spectral quality, previously to SOM and CNN classification, our 2815 BP/RP set were pre-analysed with two different approaches: the cosine similarity technique and the Uniform Manifold Approximation and Projection (UMAP) technique. Additional analysis onto a golden sample of 88 well-defined objects, is also presented. The findings demonstrate a high agreement level (sim 70-90<!PCT!>) with the classifications from the Virtual Observatory Sed Analyzer (VOSA) tool. This shows that the SVM, SOM, and CNN methods effectively classify sources with an accuracy comparable to human inspection or non-AI techniques. Notably, SVM in a radial basis function achieves 70.97<!PCT!> reproducibility for binary targets using photometry, and CNN reaches 84.94<!PCT!> for binary detection using spectroscopy. We also found that the single--binary differences are especially observable on the infrared flux in our Gaia DR3 BP/BR spectra, at wavelengths larger than $ sim $700 nm. We find that all the methods used are in fairly good agreement and are particularly effective to discern between single and binary systems. The agreement is also consistent with the results previously obtained with VOSA. In global terms, considering all quality metrics, CNN is the method that provides the best accuracy. The methods also appear effective for detecting peculiarities in the spectra. While promising, challenges in dealing with uncertain compositions highlight the need for caution, suggesting further research is needed to refine techniques and enhance automated classification reliability, particularly for large-scale surveys.

Constraining the overcontact phase in massive binary evolution. III. Period stability of known B+B and O+B overcontact systems

Binary systems play a crucial role in massive star evolution. Systems composed of B-type and O-type stars are of particular interest due to their potential to lead to very energetic phenomena or the merging of exotic compact objects. We aim to determine the orbital period variations of a sample of B+B and O+B massive overcontact binaries, with the primary objectives of characterizing the evolutionary timescales of these systems and addressing the existing discrepancy between observational data and theoretical predictions derived from population synthesis models. We used Period04 to analyze archival photometric data going back a century for a sample of seven binary systems to measure their orbital periods. We then determine the period variations using a linear fit. We find that the period variation timescales of five truly overcontact binary systems align with the nuclear timescale, in agreement with previous findings for more massive overcontact binaries. Additionally, we noticed a clear distinction between the five systems that had been unambiguously classified as overcontact systems and both SV Cen and VFTS 066, which seem to be evolving on thermal timescales and might be misclassified as overcontact systems. In the case of the five overcontact binaries, our results indicate a noticeable mismatch between the observational data and the theoretical predictions derived from population synthesis models. Furthermore, our results suggest that additional physical mechanisms must be investigated to compare the observed variations more thoroughly with theoretical predictions.

The Sunburst Arc with JWST

We report the detection of a population of Wolf-Rayet (WR) stars in the Sunburst Arc, a strongly gravitationally lensed galaxy at redshift z = 2.37. As the brightest known lensed galaxy, the Sunburst Arc has become an important cosmic laboratory for studying star and cluster formation, Lyman α (Lyα) radiative transfer, and Lyman continuum (LyC) escape. Here, we present the first results of JWST/NIRSpec IFU observations of the Sunburst Arc, focusing on a stacked spectrum of the 12-fold imaged Sunburst LyC-emitting (LCE) cluster. In agreement with previous studies, we find that the Sunburst LCE cluster is a very massive, compact star cluster with Mdyn = (9 ± 1)×106 M⊙. Our age estimate of 4.2–4.5 Myr is much larger than the crossing time of tcross = 183 ± 9 kyr, indicating that the cluster is dynamically evolved and consistent with it being gravitationally bound. We find a significant nitrogen enhancement of the low ionization state interstellar medium (ISM), with log(N/O) = − 0.74 ± 0.09, which is ≈0.8 dex above typical values for H II regions of a similar metallicity in the local Universe. We find broad stellar emission complexes around He IIλ4686 and C IVλ5808 with associated nitrogen emission; this is the first time WR signatures have been directly observed at redshifts above ∼0.5. The strength of the WR signatures cannot be reproduced by stellar population models that only include single-star evolution. While models with binary evolution better match the WR features, they still struggle to reproduce the nitrogen-enhanced WR features. JWST reveals the Sunburst LCE cluster to be a highly ionized proto-globular cluster with low oxygen abundance and extreme nitrogen enhancement that hosts a population of WR stars, likely including a previously suggested population of very massive stars (VMSs), which together are rapidly enriching the surrounding medium.

The spectroscopic binary fraction of the young stellar cluster M17

Context. Significant progress has been made toward understanding the formation of massive (M > 8 M⊙) binaries in close orbits (with periods of less than a month). Some of the observational studies leading to this progress are the detection of a very low velocity dispersion among the massive stars in the young region M17 and the measurement of a positive trend of velocity dispersion with age in Galactic clusters. The velocity dispersion observed in M17 could be explained either by the lack of binaries among the stars in this region, which implies the highly unlikely scenario of a different formation mechanism for M17 than for other Galactic regions, or by larger binary separations than typically observed, but with a binary fraction similar to other young Galactic clusters. The latter implies that, over time, the binary components migrate toward each other. This is in agreement with the finding that the radial velocity dispersion of young Galactic clusters correlates positively with their age. Aims. We aim to determine the origin of the strikingly low velocity dispersion by determining the observed and intrinsic binary fraction of massive stars in M17 through multi-epoch spectroscopy. Methods. We performed a multi-epoch spectroscopic survey consisting of three epochs separated by days and months, respectively. We complemented this survey with existing data covering timescales of years. We determined the radial velocity of each star at each epoch by fitting the stellar absorption profiles. The velocity shifts between epochs were used to determine whether a close companion is present. Results. We determined an observed binary fraction of 27% and an intrinsic binary fraction of 87%, consistent with that of other Galactic clusters. We conclude that the low velocity dispersion is due to a large separation among the young massive binaries in M17. Our result is in agreement with a migration scenario in which massive stars are born in binaries or higher-order systems at large separation and harden within the first million years of evolution. Such an inward migration may either be driven by interaction with a remnant accretion disk or with other young stellar objects present in the system, or by dynamical interactions within the cluster. Our results imply that possibly both dynamical interactions and binary evolution are key processes in the formation of gravitational wave sources.

Irradiation effect and binary radio pulsars with giant companions

Context. Pulsars are neutron stars that rotate rapidly. Most of the pulsars in binary systems tend to spin faster than those in isolation. According to binary evolution theory, radio pulsars in binary systems can have various types of companion stars. However, binary radio pulsars with giant stars, helium stars, and black holes as companions are lacking so far. Aims. We aim to investigate the possible parameter space of binary radio pulsars with giant companions. Methods. We used the MESA stellar evolution code to consider the effect of irradiation in binary evolution and evolved a series of binary models. Results. We present the potential physical properties of binary radio pulsars with giant star companions in the framework of the classical recycling scenario and the irradiation model. We found that the parameter space of binary radio pulsars with giant companions can be greatly expanded by the effect of irradiation. Moreover, when irradiation is strong enough, submillisecond pulsars might be accompanied by giant companion stars. We also demonstrate a correlation between the timescale of the binary system as a low-mass X-ray binary and the irradiation efficiency. Conclusions. The birthrate problem between millisecond pulsars and low-mass X-ray binaries might not be resolved by the irradiation effect. A more detailed binary population synthesis is necessary. Our findings may offer guidance for observers that wish to locate binary radio pulsars with giant companions or submillisecond pulsars.

Five New Heartbeat Star Systems with Tidally Excited Oscillations Discovered Based on TESS Data

Heartbeat stars (HBSs) with tidally excited oscillations (TEOs) are ideal astrophysical laboratories for studying the internal properties of the systems. In this paper, five new HBSs exhibiting TEOs are discovered using TESS photometric data. The orbital parameters are derived using a corrected version of Kumar et al.'s model based on the Markov Chain Monte Carlo method. The TEOs in these objects are examined, and their pulsation phases and modes are identified. The pulsation phases of the TEOs in TIC 266809405, TIC 266894805, and TIC 412881444 are consistent with the dominant l = 2, m = 0, or ±2 spherical harmonic. For TIC 11619404, although the TEO phase is close to the m = +2 mode, the m = 0 mode cannot be excluded because of the low inclination in this system. The TEO phase in TIC 447927324 shows a large deviation (>2σ) from the adiabatic expectations, suggesting that it is expected to be a traveling wave or that the pulsations are nonadiabatic. In addition, these TEOs occur at relatively low orbital harmonics, and we cautiously suggest that this may be an observational bias. These objects are valuable sources for studying the structure and evolution of eccentricity orbit binaries and extending the TESS HBS catalog with TEOs.

The Influence of the Magnetic Braking Laws on the Evolution of Persistent and Transient Low-mass X-Ray Binaries

Swift J1858.6−0814 (hereafter J1858) is a transient neutron star (NS) low-mass X-ray binary (LMXB). There is controversy regarding its donor mass derived from observations and theoretical calculations. In this paper, we adopt seven magnetic braking (MB) prescriptions suggested in the literature and different metallicity Z to simulate the evolution of the LMXB. Our results show that, employing the MB model proposed by A. Reiners & S. Mohanty (“rm12”), the convection- and rotation-boosted (“carb”) model, and the intermediate (“inter”) and convection-boosted (“cboost”) models in K. X. Van et al. can match (part of) the observational parameters of J1858 well. We then apply our method to other observed LMXBs and find that the “rm12” and “inter” MB laws are most promising in explaining transient LMXBs. In comparison, the simulations with the “cboost” and “carb” MB laws are more inclined to reproduce persistent LMXBs and ultracompact X-ray binaries, respectively. Our results, though subject to computational and/or observational bias, show that it is challenging to find a unified MB law that applies to the NS LMXB subpopulations simultaneously, indicating our lack of understanding of the true MB law. In addition, we explore the influence of various MB laws on the magnitude of the bifurcation periods in LMXBs.

Formation of GW230529 from Isolated Binary Evolution

In this paper, we explore the formation of the mass-gap black hole-neutron star (mgBHNS) merger detected in gravitational wave (GW) event, i.e., GW230529, from the isolated binary evolution channel, and study potential signatures of its electromagnetic counterparts. By adopting the “delayed” supernova prescription and reasonable model realizations, our population synthesis simulation results can simultaneously match the rate densities of mgBHNS and total BHNS mergers inferred from the population analyses, along with the population distribution of the BH mass in BHNS mergers reported by the LIGO–Virgo–KAGRA Collaboration. Because GW230529 contributes significantly to the inferred mgBHNS rate densities, we suggest that GW230529 can be explained through the isolated binary evolution channel. Considering the AP4 (DD2) equation of state, the probability that GW230529 can make tidal disruption is 12.8% (63.2%). If GW230529 is a disrupted event, its kilonova peak apparent magnitude is predicted ∼23–24 mag, and hence, can be detected by the present survey projects and Large Synoptic Survey Telescope. Since GW230529 could be an off-axis event inferred from the GW observation, its associated gamma-ray burst (GRB) might be too dim to be observed by γ-ray detectors, interpreting the lack of GRB observations. Our study suggests the existence of mgBHNS mergers formed through the isolated binary evolution channel due to the discovery of GW230529, indicating that BHNS mergers are still likely to be multimessenger sources that emit GWs, GRBs, and kilonovae. Although mgBHNS mergers account for ∼50% of the cosmological BHNS population, we find that ≳90% of disrupted BHNS mergers are expected to originate from mgBHNS mergers.

On the Cosmic Variance of the Merger Rate Density of Binary Neutron Stars

The cosmic variance on the star formation history may lead to bias in the merger rate density estimation of binary neutron star (BNS) mergers by compact binary population synthesis. In this paper, we take advantage of the large box size of the Millennium Simulation combined with the semianalytic galaxy formation model GABE and the parameterized population binary star evolution model to examine how much effect the cosmic variance will introduce on the estimation of the merger rate density of BNS mergers. We find that for subbox sizes of 100 and 200 Mpc, the variance of merger rate density σ R /R at different redshifts is about 23%–35% and 13%–20%, respectively. On the one hand, as for the variance of the detection rate on BNS mergers with the current LIGO–Virgo–KAGRA (LVK) detector network, this value is very small at ≲10%, which indicates ignoring the cosmic variance is reasonable for estimating the merger rate density from current LVK observations. On the other hand, with next-generation gravitational wave detectors, it is possible to localize BNS mergers within subboxes possessing a length of 40 Mpc for a source redshift z s < 0.2. In such a small box, the cosmic variance of the merger rate density is significant, i.e., the value of σ R /R is about ∼55%. This hints that estimating the merger rate density of BNS in different sky areas may provide useful information on the cosmic variance.

Constraining the Initial Mass Function via Stellar Transients

The stellar initial mass function (IMF) represents a fundamental quantity in astrophysics and cosmology describing the mass distribution of stars from low mass all the way up to massive and very massive stars. It is intimately linked to a wide variety of topics, including stellar and binary evolution, galaxy evolution, chemical enrichment, and cosmological reionization. Nonetheless, the IMF still remains highly uncertain. In this work, we aim to determine the IMF with a novel approach based on the observed rates of transients of stellar origin. We parametrize the IMF with a simple but flexible Larson shape, and insert it into a parametric model for the cosmic UV luminosity density, local stellar mass density, type Ia supernova (SN Ia), core-collapse supernova (CCSN), and long gamma-ray burst (LGRB) rates as a function of redshift. We constrain our free parameters by matching the model predictions to a set of empirical determinations for the corresponding quantities via a Bayesian Markov Chain Monte Carlo method. Remarkably, we are able to provide an independent IMF determination with a characteristic mass mc=0.10−0.08+0.24M⊙ and high-mass slope ξ=−2.53−0.27+0.24 that are in accordance with the widely used IMF parameterizations (e.g., Salpeter, Kroupa, Chabrier). Moreover, the adoption of an up-to-date recipe for the cosmic metallicity evolution allows us to constrain the maximum metallicity of LGRB progenitors to Zmax=0.12−0.05+0.29Z⊙. We also find which progenitor fraction actually leads to SN Ia or LGRB emission (e.g., due to binary interaction or jet-launching conditions), put constraints on the CCSN and LGRB progenitor mass ranges, and test the IMF universality. These results show the potential of this kind of approach for studying the IMF, its putative evolution with the galactic environment and cosmic history, and the properties of SN Ia, CCSN, and LGRB progenitors, especially considering the wealth of data incoming in the future.

Contribution of bismuth melts to gold endowment in the Baolun gold deposit, Hainan Island, South China

Bismuth (Bi) melts and related polymetallic alloys can efficiently scavenge gold (Au) from Au-unsaturated aqueous solutions, contributing to the Au endowment in many hydrothermal deposits. However, original evidence for Au–Bi melts is often poorly preserved due to post-precipitation alteration. This complicates investigation of the liquid bismuth collector model in hydrothermal Au deposits. Here we present primary evidence for the occurrence of Au–Bi melt in hydrothermal systems through an examination of the Au–Bi phases and textures in the Baolun Au deposit (Hainan Island). This study found that maldonite, native bismuth, Au–Bi symplectite and Au–Bi melt blebs appeared successively following the precipitation of native gold. Notably, large Au–Bi melt blebs were well preserved in natural systems due to the rapid cooling of fluids. The mineral assemblages and their corresponding fluid physiochemical conditions suggest that the evolution of the ore fluids at Baolun was characterized by a continuous reduction in Au and Bi concentrations and Au/Bi ratios alongside decreasing fluid temperatures and sulfur fugacity (fS2). These observations offer direct evidence for Au scavenging by Bi melts in a mesothermal (275–450 °C), low-fS2 and reduced hydrothermal system, aligning well with the Au–Bi binary phase evolution established in metallurgy. As the fluid cooled further, the Au–Bi phases were subsequently overprinted by later Te–S-rich fluids, as evidenced by the formation of bismuthinite, jonassonite, and joséite-A around the Au–Bi phases. Importantly, our study reveals that the Baolun Au deposit ischaracterizedby a Au–Bi–Te hydrothermal system and that the metamorphic country rocks at Baolun are the major source of Au, Bi, and Te.

Exploring the boundary between stable mass transfer and L2 overflow in close binary evolution

The majority of massive stars resides in binary systems, which are expected to experience mass transfer during their evolution. However, the conditions under which mass transfer leads to a common envelope, and thus possibly to a merging of both stars, are currently only poorly understood. The main uncertainties arise from the possible swelling of the mass gainer and from angular momentum loss from the binary system during non-conservative mass transfer. We have computed a dense grid of detailed models of stars that accrete mass at constant rates to determine the radius increase that is due to their thermal disequilibrium. While we find that models with an accretion that is faster than the thermal timescale expand in general, this expansion remains quite limited in the intermediate-mass regime even for accretion rates that exceed the thermal timescale accretion rate by a factor of 100. Our models of massive stars expand to extreme radii under these conditions. When the accretion rate exceed the Eddington accretion rate, our models expand rapidly. We derived analytical fits to the radius evolution of our models and a prescription for the boundary between stable mass transfer and overflow for arbitrary accretion efficiencies. We then applied our results to grids of binary models adopting various constant mass-transfer efficiencies and angular momentum budgets. We find that the first parameter affects the outcome of the Roche-lobe overflow more strongly. Our results are consistent with detailed binary evolution models and often lead to a smaller initial parameter space for stable mass transfer than do other recipes in the literature. We used this method to investigate the origin of Wolf-Rayet stars with O star companions in the Small Magellanic Cloud, and we found that the efficiency of the mass transfer process that led to the formation of the Wolf-Rayet star was likely lower than 50&lt;!PCT!&gt;.