Search

We investigate the effect of the cutting-edge circumbinary disc (CBD) evolution models on massive black hole binary (MBHB) populations and the gravitational wave background (GWB). We show that CBD-driven evolution leaves a tell-tale signature in MBHB populations, by driving binaries towards an equilibrium eccentricity that depends on the binary mass ratio. We find high orbital eccentricities (⁠|$e_{\rm b} \sim 0.5$|⁠) as MBHBs enter multimessenger observable frequency bands. The CBD-induced eccentricity distribution of MBHB populations in observable bands is independent of the initial eccentricity distribution at binary formation, erasing any memory of eccentricities induced in the large-scale dynamics of merging galaxies. Our results suggest that eccentric MBHBs are the rule rather than the exception in upcoming transient surveys, provided that CBDs regularly form in MBHB systems. We show that the GWB amplitude is sensitive to CBD-driven preferential accretion onto the secondary, resulting in an increase in GWB amplitude |$A_{\rm yr^{-1}}$| by over 100 per cent with just 10 per cent Eddington accretion. As we self-consistently allow for binary hardening and softening, we show that CBD-driven orbital expansion does not diminish the GWB amplitude, and instead increases the amplitude by a small amount. We further present detection rates and population statistics of MBHBs with |$M_{\rm b} \gtrsim 10^6 \, {\rm M}_{\odot }$| in Laser Interferometer Space Antenna , showing that most binaries have equal mass ratios and can retain residual eccentricities up to |$e_{\rm b} \sim 10^{-3}$| due to CBD-driven evolution. [ABSTRACT FROM AUTHOR]

Galactic white dwarf binaries (WDBs) and black hole binaries (BHBs) will be gravitational wave (GW) sources for LISA. Their detection will provide insights into binary evolution and the evolution of our Galaxy through cosmic history. Here, we make predictions of the expected WDB and BHB population within our Galaxy. We combine predictions of the compact remnant binary populations expected by stellar evolution from the detailed Binary Population and Spectral Synthesis (BPASS) code, with a Milky Way analogue galaxy model from the Feedback in Realistic Environment (FIRE) simulations. We use PhenomA and legwork to simulate LISA observations. Both packages make similar predictions that on average four Galactic BHBs and 673 Galactic WDBs are above the signal-to-noise ratio (SNR) threshold of 7 after a four-year mission. We compare these predictions to earlier results using the binary star evolution (BSE) code with the same FIRE model galaxy. We find that BPASS predicts a few more LISA observable Galactic BHBs and a twentieth of the Galactic WDBs. The differences are due to the different physical assumptions that have gone into the binary evolution calculations. These results indicate that the expected population of compact binaries that LISA will detect depends very sensitively on the binary population synthesis models used and thus observations of the LISA population will provide tight constraints on our modelling of binary stars. Finally, from our synthetic populations, we have created mock LISA signals that can be used to test and refine data processing methods of the eventual LISA observations. [ABSTRACT FROM AUTHOR]

Millisecond pulsars with white dwarf companions have typical eccentricities |$e\sim 10^{-6}{\!-\!}10^{-3}$|⁠. The eccentricities of helium white dwarfs are explained well by applying the fluctuation–dissipation theorem to convective eddies in their red giant progenitors. We extend this theory to more massive carbon–oxygen (CO) white dwarfs with asymptotic giant branch (AGB) progenitors. Due to the radiation pressure in AGB stars, the dominant factor in determining the remnant white dwarf's eccentricity is the critical residual hydrogen envelope mass |$m_{\rm env}$| required to inflate the star to giant proportions. Using a suite of mesa stellar evolution simulations with |$\Delta m_{\rm c}=10^{-3}\, {\rm M}_{\odot }$| core-mass intervals, we resolved the AGB thermal pulses and found that the critical |$m_{\rm env}\propto m_{\rm c}^{-6}$|⁠. The resulting eccentricity |$e\sim 3\times 10^{-3}$| is almost independent of the remnant CO white dwarf's mass |$m_{\rm c}$|⁠. Nearly all of the measured eccentricities lie below this robust theoretical limit, indicating that the eccentricity is damped during the common-envelope inspiral that follows the unstable Roche lobe overflow of the AGB star. Specifically, we focused on white dwarfs with median masses |$m_{\rm c}\gt 0.6\, {\rm M}_{\odot }$|⁠. These massive white dwarfs begin their inspiral with practically identical orbital periods and eccentricities, eliminating any dependence on the initial conditions. For this sub-sample, we find an empirical relation |$e\propto P^{3/2}$| between the final period and eccentricity that is much tighter than previous studies – motivating theoretical work on the eccentricity evolution during the common envelope phase. The eccentricities of lower mass CO white dwarfs may be explained by alternative formation channels. [ABSTRACT FROM AUTHOR]

The Wide Angle Search for Planets (WASP) survey used transit photometry to discover nearly 200 gas-giant exoplanets and derive their planetary and stellar parameters. Reliable determination of the planetary density depends on accurate measurement of the planet's radius, obtained from the transit depth and photodynamical determination of the stellar radius. The stellar density and hence the stellar radius are typically determined in a model-independent way from the star's reflex orbital acceleration and the transit profile. Additional flux coming from the system due to a bright, undetected stellar binary companion can, however, potentially dilute the transit curve and radial velocity signal, leading to underestimation of the planet's mass and radius, and to overestimation of the planet's density. In this study, we cross-check the published radii of all the WASP planet-host stars, determined from their transit profiles and radial velocity curves, against radiometric measurements of stellar radii derived from their angular diameters (via the infrared flux method) and trigonometric parallaxes. We identify eight systems showing radiometric stellar radii significantly greater than their published photodynamical values: WASPs 20, 85, 86, 103, 105, 129, 144, and 171. We investigate these systems in more detail to establish plausible ranges of angular and radial velocity separations within which such 'stealth binaries' could evade detection, and deduce their likely orbital periods, mass ratios, and flux ratios. After accounting for the dilution of transit depth and radial velocity amplitude, we find that, on average, the planetary densities for the identified stealth binary systems should be reduced by a factor of 1.3. [ABSTRACT FROM AUTHOR]

Relativistic discrete ejecta launched by black hole X-ray binaries (BH XRBs) can be observed to propagate up to parsec-scales from the central object. Observing the final deceleration phase of these jets is crucial to estimate their physical parameters and to reconstruct their full trajectory, with implications for the jet powering mechanism, composition, and formation. In this paper, we present the results of the modelling of the motion of the ejecta from three BH XRBs: MAXI J1820 |$+$| 070, MAXI J1535–571, and XTE J1752–223, for which high-resolution radio and X-ray observations of jets propagating up to |$\sim$| 15 arcsec (⁠|$\sim$| 0.6 pc at 3 kpc) from the core have been published in the recent years. For each jet, we modelled its entire motion with a dynamical blast-wave model, inferring robust values for the jet Lorentz factor, inclination angle and ejection time. Under several assumptions associated to the ejection duration, the jet opening angle and the available accretion power, we are able to derive stringent constraints on the maximum jet kinetic energy for each source (between |$10^{43}$| and |$10^{44}$| erg, including also H1743–322), as well as placing interesting upper limits on the density of the ISM through which the jets are propagating (from |$n_{\rm ISM} \lesssim 0.4$| cm |$^{-3}$| down to |$n_{\rm ISM} \lesssim 10^{-4}$| cm |$^{-3}$|⁠). Overall, our results highlight the potential of applying models derived from gamma-ray bursts to the physics of jets from BH XRBs and support the emerging picture of these sources as preferentially embedded in low-density environments. [ABSTRACT FROM AUTHOR]

Pair-instability supernovae (PISNe) have long been predicted to be the final fates of near-zero-metallicity very massive stars (⁠|$Z \lt Z_\odot /3$|⁠ , M |$_\mathrm{ZAMS} \gtrsim 140\, \text{M}_\odot$|⁠). However, no definite PISN has been observed to date, leaving theoretical modelling validation open. To investigate the observability of these explosive transients, we combine detailed stellar evolution models for PISNe formation, computed from the binary population and spectral synthesis code suite, bpass , with the star formation history of all individual computational elements in the Illustris-TNG simulation. This allows us to compute comic PISN rates and predict their host galaxy properties. Of particular importance is that IllustrisTNG galaxies do not have uniform metallicities throughout, with metal-enriched galaxies often harbouring metal-poor pockets of gas where PISN progenitors may form. Accounting for the chemical inhomogeneities within these galaxies, we find that the peak redshift of PISNe formation is |$z=3.5$| instead of the value of |$z=6$| when ignoring chemical inhomogeneities within galaxies. Furthermore, the rate increases by an order of magnitude from 1.9 to 29 PISN Gpc |$^{-3}$| yr |$^{-1}$| at |$z=0$|⁠ , if the chemical inhomogeneities are considered. Using state-of-the-art theoretical PISN light curves, we find an observed rate of 13.8 (1.2) visible PISNe per year for the Euclid-Deep survey, or 83 (7.3) over the 6-yr lifetime of the mission when considering chemically inhomogeneous (homogenous) systems. Interestingly, only 12 per cent of helium PISN progenitors are sufficiently massive to power a superluminous supernova event, which can potentially explain why PISN identification in time-domain surveys remains elusive and progress requires dedicated strategies. [ABSTRACT FROM AUTHOR]

We present the Ultracool Dwarf Companion Catalogue of 278 multiple systems, 32 of which are newly discovered, each with at least one spectroscopically confirmed Ultracool Dwarf, within a 100 pc volume-limited sample. This catalogue is compiled using the Gaia Catalogue of Nearby Stars for stellar primaries and the Gaia Ultracool Dwarf Sample for low-mass companions and includes 241 doubles, 33 triples, and 4 higher-order systems established from positional, proper motion, and parallax constraints. The catalogue seeks to identify probable benchmark systems within 100 pc to obtain model-independent astrophysical parameters of Ultracool Dwarfs. Chance alignment probabilities are calculated to evaluate the physical nature of each system. Astrometric and photometric data from Gaia Data Release 3 and the Two Micron All Sky Survey are included for all objects. We identify potential unseen companions using a combination of the renormalized unit weight error, image parameter determination statistics, non-single star solutions, and photometric blending as provided by Gaia , identifying hierarchical Ultracool triple systems. Our catalogue includes 17 White Dwarf – Ultracool Dwarf systems, whose ages are determined using cooling models. We also use the Gaia FLAME results and the BANYAN Σ procedures to age 40 and 34 systems, respectively, and derive mass estimates from evolutionary models. [ABSTRACT FROM AUTHOR]

The correction map method means extended phase-space algorithm with correction map. In our research, we have developed a correction map method, specifically the dissipated correction map method with trapezoidal rule, for numerical simulations of gravitational waves from spinning compact binary systems. This new correction map method, denoted as |$CM3$|⁠ , has shown remarkable performance in various simulation results, such as phase-space distance, dissipated energy error, and gravitational waveform, closely resembling the high-order precision implicit Gaussian algorithm. When compared with the previously used mid-point map which is denoted as |$C_2$|⁠ , the |$CM3$| consistently exhibits a closer alignment with the highly accurate Gaussian algorithm in waveform evolution and orbital trajectory analysis. Through detailed comparisons and analyses, it is evident that |$CM3$| outperforms other algorithms, including |$CM2$| and |$C_2$| mentioned in this paper, in terms of accuracy and precision in simulating spinning compact binary systems. The incorporation of the trapezoidal rule and the optimization with a scale factor |$\gamma$| have significantly enhanced the performance of |$CM3$|⁠ , making it a promising method for future numerical simulations in astrophysics. With the groundbreaking detection of gravitational waves by the LIGO/VIRGO collaboration, interest in this research domain has soared. Our work contributes valuable insights for the application of matched filtering techniques in the analysis of gravitational wave signals, enhancing the precision and reliability of these detection. [ABSTRACT FROM AUTHOR]

It has been generally accepted that the originators of the double star astronomy were Christian Mayer and William Herschel. We recovered the memory of the poorly known Italian astronomer Giovanni Battista Hodierna, who published the first catalogue of stellar binaries over a century before Mayer and Herschel. We analysed the fourth section of 1654 G. B. Hodierna's book ' De systemate orbis cometici deque admirandis coeli characteribus '. There, Hodierna listed a dozen pairs of stars whose identification with modern star names had been lost for centuries. To identify the pairs, we used Hodierna's Latin descriptions of location in constellations for all primary stars, ecliptic coordinates and angular separations to companions for some, and the Washington Double Star, Hipparcos , and Gaia catalogues. We were able to identify the twelve primaries and eleven multiple systems with companions, of which nine were double and two were triple. Besides, with up-to-date data, we confirmed that four systems are physically bound: Atlas and Pleione, |$\alpha ^{1,2}$| Lib, |$\nu ^{1,2}$| Dra, and |$\theta ^{1}$| Ori A, C, and D. The other seven pairs are alignments of very bright stars at different distances. [ABSTRACT FROM AUTHOR]

Hypervelocity stars (HVSs) are stars which have been ejected from the Galactic Centre (GC) at velocities of up to a few thousand |$\text{km}\, \text{s}^{-1}$|⁠. They are tracers of the Galactic potential and can be used to infer properties of the GC, such as the initial mass function and assembly history. HVSs are rare, however, with only about a dozen promising candidates discovered so far. In this work, we make use of a novel, highly efficient method to identify new HVS candidates in Gaia. This method uses the nearly radial trajectories of HVSs to infer their distances and velocities based on their position and Gaia proper motion alone. Through comparison of inferred distances with Gaia parallaxes and photometry, we identified 600 HVS candidates with G < 20 including the previously discovered S5-HVS1, out of which we obtained ground-based follow-up observations for 196 stars. As we found no new HVSs based on their radial velocity, we used detailed HVS ejection simulations to significantly improve previous HVS ejection rate constraints. In particular, the ejection rate of HVSs more massive than 1 |$\mathrm{M_\odot }$| cannot be higher than |$10^{-5}$| yr |$^{-1}$| at |$2\sigma$| significance. Additionally, we predict that there are 5–45 unbound HVSs in the complete Gaia catalogue (⁠|$1\sigma$| interval), most of which will be main-sequence stars of a few M |$_\odot$| at heliocentric distances of tens to hundreds of kpc. By comparing our results to literature HVS candidates, we find an indication of either a time-dependent ejection rate of HVSs or a non-GC origin of previously identified HVS candidates. [ABSTRACT FROM AUTHOR]

Surveys in the Milky Way and Large Magellanic Cloud have revealed that the majority of massive stars will interact with companions during their lives. However, knowledge of the binary properties of massive stars at low metallicity, and therefore in conditions approaching those of the Early Universe, remain sparse. We present the Binarity at LOw Metallicity (BLOeM) campaign, an ESO large programme designed to obtain 25 epochs of spectroscopy for 929 massive stars in the Small Magellanic Cloud, allowing us to probe multiplicity in the lowest-metallicity conditions to date (Z = 0.2 Z⊙). BLOeM will provide (i) the binary fraction, (ii) the orbital configurations of systems with periods of P ≲ 3 yr, (iii) dormant black-hole binary candidates (OB+BH), and (iv) a legacy database of physical parameters of massive stars at low metallicity. Main sequence (OB-type) and evolved (OBAF-type) massive stars are observed with the LR02 setup of the GIRAFFE instrument of the Very Large Telescope (3960–4570 Å resolving power R = 6200; typical signal-to-noise ratio(S/N) ≈70–100). This paper utilises the first nine epochs obtained over a three-month time period. We describe the survey and data reduction, perform a spectral classification of the stacked spectra, and construct a Hertzsprung-Russell diagram of the sample via spectral-type and photometric calibrations. Our detailed classification reveals that the sample covers spectral types from O4 to F5, spanning the effective temperature and luminosity ranges 6.5 ≲ Teff/kK ≲ 45 and 3.7 < log L/L⊙ < 6.1 and initial masses of 8 ≲ Mini ≲ 80 M⊙. The sample comprises 159 O-type stars, 331 early B-type (B0–3) dwarfs and giants (luminosity classes V–III), 303 early B-type supergiants (II–I), and 136 late-type BAF supergiants. At least 82 stars are OBe stars: 20 O-type and 62 B-type (13% and 11% of the respective samples). In addition, the sample includes 4 high-mass X-ray binaries, 3 stars resembling luminous blue variables, 2 bloated stripped-star candidates, 2 candidate magnetic stars, and 74 eclipsing binaries. [ABSTRACT FROM AUTHOR]

This work aims to analyze some of the polluters proposed in the self-enrichment scenarios put forward to explain the multiple populations in globular clusters (GCs), extending previous studies. Three scenarios with different polluter stars were tested: asymptotic giant branch stars (AGBs), high-mass interacting binaries (IBs), and fast rotating massive stars (FRMSs). With abundance data available from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey and ΔY estimates from precise Hubble Space Telescope (HST) photometry, twenty-six clusters were studied, increasing the number of clusters in previous studies by more than a factor of three. We also included the study of the abundances of N, C, Mg, and Al, extending previous studies that mainly focused on the abundances of He, O, and Na. In addition, we constructed an empirical model to test whether one could explain the chemical signatures of the "enriched" population of GC stars with a fixed source and dilution process based on empirical data. In agreement with work by other authors, we found that the proposed polluters can generally predict the qualitative abundance patterns in GC stars and in some cases quantitatively predict some elements, but in most cases when we compare the model yields with the observations, we find that they cannot explain the entire set of observed abundance patterns. The empirical model succeeds in reproducing the abundances of Al for a given ΔY (and vice versa), showing that there is a direct relationship between Al and He, with one increasing proportionally to the other. However, the empirical model fails to reproduce the observed abundances of Na and N, in agreement with the results of previous works. The observed decoupling between the maximum abundances of CNO-cycle elements such as N and Na with those of Al and He provides new information and constraints for future models and could take us a step closer to understanding the origin of the peculiar abundance variations of GC stars. [ABSTRACT FROM AUTHOR]

Accreting neutron stars (NSs) are expected to emit a redshifted 2.2 MeV line due to the capture of neutrons produced through the spallation processes of 4He and heavier ions in their atmospheres. Detecting this emission would offer an independent method for constraining the equation of state of NSs and provide valuable insights into nuclear reactions occurring in extreme gravitational and magnetic environments. Typically, a higher mass accretion rate is expected to result in a higher 2.2 MeV line intensity. However, when the mass accretion rate approaches the critical threshold, the accretion flow is decelerated by the radiative force, leading to a less efficient production of free neutrons and a corresponding drop in the flux of the spectral line. This makes the brightest X-ray pulsars unsuitable candidates for γ-ray line detection. In this work, we present a theoretical framework for predicting the optimal X-ray luminosity required to detect a redshifted 2.2 MeV line in a strongly magnetized NS. As the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL) mission nears its conclusion, we have undertaken a thorough investigation of the SPectrometer on board INTEGRAL (SPI) data of this line in a representative sample of accreting NSs. No redshifted 2.2 MeV line was detected. For each spectrum, we have determined the 3σ upper limits of the line intensity, assuming different values of the line width. Although the current upper limits are still significantly above the expected line intensity, they offer valuable information for designing future gamma-ray telescopes aimed at bridging the observational MeV gap. Our findings suggest that advancing our understanding of the emission mechanism of the 2.2 MeV line, as well as the accretion flow responsible for it, will require a substantial increase in sensitivity from future MeV missions. For example, for a bright X-ray binary such as Sco X−1, we would need at least a 3σ line point source sensitivity of ≈10−6 ph cm−2 s−1, that is, about two orders of magnitude better than that currently achieved. [ABSTRACT FROM AUTHOR]

Context. High-precision and high-cadence photometric surveys such as Kepler or TESS are making huge progress not only in the detection of new extrasolar planets but also in the study of a great number of variable stars. This is the case for central stars of planetary nebulae (PNe), which have similarly benefited from the capabilities of these missions, increasing the number of known binary central stars and helping us to constrain the relationship between binarity and the complex morphologies of their host PNe. Aims. In this paper, we analyse the TESS light curves of a large sample of central stars of PNe with the aim of detecting signs of variability that may hint at the presence of short-period binary nuclei. This will have important implications in understanding PN formation evolution as well as the common envelope phase. Methods. We analysed 62 central stars of true, likely, or possible PNe and modelled the detected variability through an MCMC approach accounting for three effects: reflection, ellipsoidal modulations due to tidal forces, and the so-called Doppler beaming. Among the 62 central stars, only 38 are amenable for this study. The remaining 24 show large contamination from nearby sources preventing an optimal analysis. Also, eight targets are already known binary central stars, which we revisit here with the new high precision of the TESS data. Results. In addition to recovering the eight already known binaries in our sample, we find that 18 further central stars show clear signs of periodic variability in the TESS data, probably resulting from different physical effects compatible with the binary scenario. We propose them as new candidate binary central stars. We also discuss the origin of the detected variability in each particular case by using the TESS_localize algorithm. Finally, 12 targets show no or only weak evidence of variability at the sensitivity of TESS. Conclusions. Our study demonstrates the power of space-based photometric surveys in searching for close binary companions of central stars of PNe. Although our detections can only be catalogued as candidate binaries, we find a large percentage of possible stellar pairs associated with PNe, supporting the hypothesis that binarity plays a key role in shaping these celestial structures. [ABSTRACT FROM AUTHOR]

PSR J1227−6208 is a 34.53-ms recycled pulsar with a massive companion. This system has long been suspected to belong to the emerging class of massive recycled pulsar−ONeMg white dwarf systems such as PSR J2222−0137, PSR J1528−3146, and J1439−5501. Here, we present an updated emission and timing analysis with more than 11 years of combined Parkes and MeerKAT data, including 19 hours of high-frequency data from the newly installed MeerKAT S-band receivers. We measure a scattering timescale of 1.22 ms at 1 GHz with a flat scattering index of 3.33 < β < 3.62, and a mean flux density of 0.53 − 0.62 mJy at 1 GHz with a steep spectral index of 2.06 < α < 2.35. Around 15% of the emission is linearly and circularly polarised, but the polarisation angle does not follow the rotating vector model. Thanks to the sensitivity of MeerKAT, we successfully measure a rate of periastron advance of ω7 = 0.0171(11) deg yr−1, and a Shapiro delay with an orthometric amplitude of h3 = 3.6 ± 0.5 μs and an orthometric ratio of ς = 0.85 ± 0.05. The main source of uncertainty in our timing analysis is chromatic correlated dispersion measure noise, which we model as a power law in the Fourier space thanks to the large frequency coverage provided by the Parkes UWL receiver. Assuming general relativity and accounting for the measurements across all the implemented timing noise models, the total mass, companion mass, pulsar mass, and inclination angle are constrained at 2.3 < Mt/M⊙ < 3.2, 1.21 < Mc/M⊙ < 1.47, 1.16 < Mp/M⊙ < 1.69, and 77.5 < i/deg < 80.3. We also constrain the longitude of ascending node to either Ωa = 266 ± 78 deg or Ωa = 86 ± 78 deg. We argue against a neutron star nature of the companion based on the very low orbital eccentric of the system (e = 1.15 × 10−3), and instead classify the companion of PSR J1227−6208 as a rare, massive ONeMg white dwarf close to the Chandrasekhar limit. [ABSTRACT FROM AUTHOR]

There could be a significant population of double white dwarf binaries (DWDs) inside globular clusters (GCs); however, these binaries are often too faint to be individually observed. We have utilized a large number GC models evolved with the Monte Carlo Cluster Simulator (MOCCA) code to create a large statistical dataset of DWDs. These models include multiple-stellar populations, resulting in two distinct initial populations: one dense and the other less dense. Due to the lower density of one population, a large number of objects escape during the early GC evolution, leading to a high mass-loss rate. In this dataset we have analyzed three main groups of DWDs, namely in-cluster binaries, escaped binaries, and binaries formed from the isolated evolution of primordial binaries. We compared the properties of these groups to observations of close and wide binaries. We find that the number of escaping DWDs is significantly larger than the number of in-cluster binaries and those that form via the isolated evolution of all primordial binaries in our GC models. This suggests that dynamics play an important role in the formation of DWDs. For close binaries, we found a good agreement in the separations of escaped binaries and isolated binaries, but in-cluster binaries showed slight differences. We could not reproduce the observed extremely low mass WDs due to the limitations of our stellar and binary evolution prescriptions. For wide binaries, we also found a good agreement in the separations and masses, after accounting for observational selection effects. Even though the current observational samples of DWDs are extremely biased and incomplete, we conclude that our results compare reasonably well with observations. [ABSTRACT FROM AUTHOR]

Context. Many fast-rotating stars (rotation periods of < 2 days) are found to be unresolved binaries with separations of tens of AU. This correlation between fast rotators and binarity leads to the question of whether the formation of binary stars inherently produces fast rotators. Aims. Our goal is to understand the spin evolution of protostars and whether the formation of companions plays a role in spinning up stars. Methods. We used magneto-hydrodynamic simulations to study the formation of multiple star systems from turbulent and non-turbulent protostellar cores. We tracked the angular momentum accreted by individual star and inner disc systems by using a sink (star) particle technique. We ran a resolution study to extrapolate protostellar properties. Results. We find in all simulations that the primary star can experience a spin-up event correlated with the formation of companions, namely fragmentation into binaries or higher-order systems. The primary star can spin up by up to 84% of its pre-fragmentation angular momentum and by up to 18% of its pre-fragmentation mass-specific angular momentum. The mechanism for the spin-up is gravitational disc instabilities in the circumstellar disc around the primary star, which leads to the accretion of material with high specific angular momentum. The simulations that experience the strongest disc instabilities fragment to form companions. Simulations with weaker spin-up events experience disc instabilities triggered by a companion flyby, and the disc instability in these cases typically does not produce further fragments (i.e. they remain binary systems). Conclusions. The primary star in multiple star systems can end up with a higher spin than single stars. This is because gravitational instabilities in the circumstellar disc around the primary star can trigger a spin-up event. In the strongest spin-up events, the instability is likely to cause disc fragmentation and the formation of companions. This spin-up mechanism, coupled with shorter disc lifetimes due to truncated circumstellar discs (and thus short spin-down times), may help produce fast rotators. [ABSTRACT FROM AUTHOR]

Context. The detected population of "spiders", referring to millisecond pulsar binaries, has significantly grown in the past decade thanks to multiwavelength follow-up investigations of unidentified Fermi sources. These systems consist of low-mass stellar companions orbiting rotation-powered millisecond pulsars in short periods of a few hours up to day. Among them, a subset of intriguing objects called transitional millisecond pulsars (tMSPs) has been shown to exhibit a remarkable behavior, transitioning between pulsar-binary and faint low-mass X-ray binary states over a span of a few years. Aims. Our objective is to study the interaction of stellar winds in tMSPs in order to understand their observational properties. To this end we focus on the parameter range that places the system near Roche-lobe overflow. Methods. Employing the adaptative mesh refinement (AMR) AMRVAC 2.0 code, we performed 2D hydrodynamical (HD) simulations of the interaction between the flows from both stars, accounting for the effects of gravity and orbital motion. Results. By studying the mass loss and launch speed of the winds, we successfully recreated two phenomenologically distinct regimes: the accretion stream and the pulsar radio state. We also identified the tipping point that marks the sharp transition between these two states. In the accretion stream state, we discover a very strong variability induced by the pulsar wind. In the pulsar state, we reconstructed the corresponding X-ray light curves of the system that produces the characteristic double-peak pattern of these systems. The position of the peaks is shifted due to orbital motion and the leading peak is weaker due to eclipsing by the companion. Conclusions. This study highlights the importance of gravity and orbital motion in the interaction between the companion and pulsar winds. Our setup allows the study of the complex interaction between the pulsar wind and an accretion stream during mass transfer. We suggest that a smaller leading peak in X-rays is indicative of a nearly edge-on system. [ABSTRACT FROM AUTHOR]

We report the discovery of the new accreting millisecond X-ray pulsar SRGA J144459.2−604207 using data of the SRG/ART-XC. The source was observed twice in February 2024 during the declining phase of the outburst. The timing analysis revealed a coherent signal near 447.9 Hz modulated by the Doppler effect due to the orbital motion. The derived parameters for the binary system are consistent with a circular orbit with a period of ∼5.2 h. The pulse profiles of the persistent emission, showing a sine-like part during half a period with a plateau in between, can be well modeled by emission from two circular spots that are partially eclipsed by the accretion disk. Additionally, during our observations with an exposure of 133 ks, we detected 19 thermonuclear X-ray bursts. All bursts have similar shapes and energetics, and none show any signs of an expanding photospheric radius. The burst recurrence times decreases linearly from ∼1.6 h at the beginning of observations to ∼2.2 h at the end and anticorrelate with the persistent flux. The spectral evolution during the bursts is consistent with the models of the neutron star atmospheres that are heated by accretion and implies a neutron star radius of 11–12 km and a distance to the source of 8–9 kpc. We also detected coherent pulsations during the bursts and showed that the pulse profiles differ substantially from those observed in the persistent emission. However, we could not find a simple physical model explaining the pulse profiles detected during the bursts. [ABSTRACT FROM AUTHOR]

Gaia BH3 is the first observed dormant black hole (BH) with a mass of ≈30 M⊙, and it represents the first confirmation that such massive BHs are associated with metal-poor stars. Here, we explore the isolated binary formation channel for Gaia BH3, focusing on the old and metal-poor stellar population of the Milky Way halo. We used the MIST stellar models and our open-source population synthesis code SEVN to evolve 5.6 × 108 binaries, exploring 20 sets of parameters that encompass different natal kicks, metallicities, common envelope efficiencies and binding energies, and models for the Roche-lobe overflow. We find that systems such as Gaia BH3 form preferentially from binaries initially composed of a massive star (40–60 M⊙) and a low-mass companion (<1 M⊙) in a wide (P > 103 days) and eccentric orbit (e > 0.6). Such progenitor binary stars do not undergo any Roche-lobe overflow episode during their entire evolution, so the final orbital properties of the BH-star system are essentially determined at the core collapse of the primary star. Low natal kicks (≲ 10 km/s) significantly favour the formation of Gaia BH3-like systems, but high velocity kicks up to ≈220 km/s are also allowed. We estimated the formation efficiency for Gaia BH3-like systems in old (t >10 Gyr) and metal-poor (Z < 0.01) populations to be ∼4 × 10−8 M⊙−1 (for our fiducial model), representing ~3% of the whole simulated BH-star population. We expect up to ≈4000 BH-star systems in the Galactic halo formed through isolated evolution, of which ≈100 are compatible with Gaia BH3. Gaia BH3-like systems represent a common product of isolated binary evolution at low metallicity (Z < 0.01), but given the steep density profile of the Galactic halo, we do not expect more than one at the observed distance of Gaia BH3. Our models show that even if it was born inside a stellar cluster, Gaia BH3 is compatible with a primordial binary star that escaped from its parent cluster without experiencing significant dynamical interactions. [ABSTRACT FROM AUTHOR]

Products of stellar mergers are predicted to be common in stellar populations and can potentially explain stars with peculiar properties. When the merger occurs after the initially more massive star has evolved into the Hertzsprung gap, the merger product may remain in the blue part of the Hertzsprung–Russell diagram for millions of years. Such objects could, therefore, explain the overabundance of observed blue stars, such as blue supergiants. However, it is currently not straightforward to distinguish merger products from genuine single stars or other stars with similar surface diagnostics. In this work, we made detailed asteroseismic comparisons between models of massive post-main-sequence merger products and genuine single stars to identify which asteroseismic diagnostics can be used to distinguish them. In doing so, we developed tools for the relatively young field of merger seismology. Genuine single stars in the Hertzsprung gap are fully radiative, while merger products have a convective He-burning core and convective H-burning shell while occupying similar locations in the Hertzsprung–Russell diagram. These major structural differences are reflected in lower asymptotic period spacing values for merger products and the appearance of deep dips in their period spacing patterns. Our genuine single-star models with masses above roughly 11.4 solar masses develop short-lived intermediate convective zones during their Hertzsprung gap evolution. This also leads to deep dips in their period spacing patterns. Because of the lack of a convective core, merger products and genuine single stars can be distinguished based on their asymptotic period spacing value in this mass range. We performed the comparisons with and without the effects of slow rotation included in the pulsation equations and conclude that the two types of stars are seismically distinguishable in both cases. The observability of the distinguishing asteroseismic features of merger products can now be assessed and exploited in practice. [ABSTRACT FROM AUTHOR]

A consensus prevails with regard to star-disk systems accreting most of their mass and angular momentum during the collapse of a prestellar core. However, recent results have indicated that stars experience post-collapse or late infall, during which the star and its disk are refreshed with material from the protostellar environment through accretion streamers. Apart from adding mass to the star-disk system, infall potentially supplies a substantial amount of angular momentum, as the infalling material is initially not bound to the collapsing prestellar core. We investigate the orientation of infall on star-disk systems by analyzing the properties of accreting tracer particles in three-dimensional magnetohydrodynamical (3D MHD) simulations of a molecular cloud that is (4 pc)3 in volume. In contrast to the traditional picture, where the rotational axis is inherited from the collapse of a coherent pre-stellar core, the orientation of star-disk systems changes substantially throughout the accretion process, thereby extending the possibility of primordial misalignment as the source of large obliquities. In agreement with previous results that show larger contributions of late infall for increasing stellar masses, a misaligned infall is more likely to lead to a prolonged change in orientation for stars of higher final mass. On average, brown dwarfs and very low mass stars are more likely to form and accrete all of their mass as part of a multiple system, while stars with final masses above a few 0.1 M⊙ are more likely to accrete part of their mass as single stars. Finally, we find an overall trend among our sample: the post-collapse accretion phase is more anisotropic than the early collapse phase. This result is consistent with a scenario of Bondi-Hoyle-Littletlon accretion during the post-collapse phase, while the initial collapse is less anisotropic – despite the fact that material is funneled through accretion channels. [ABSTRACT FROM AUTHOR]

Massive binary evolution models are needed to predict massive star populations in star-forming galaxies, the supernova diversity, and the number and properties of gravitational wave sources. Such models are often computed using so-called rapid binary evolution codes, which approximate the evolution of the binary components based on detailed single star models. However, about one-third of the interacting massive binary stars undergo mass transfer during core hydrogen-burning (Case A mass transfer), whose outcome is difficult to derive from single star models. For this work, we used a large grid of detailed binary evolution models for primaries in the initial mass range 10–40 M⊙ with a Large and Small Magellanic Cloud composition, to derive analytic fits for the key quantities needed in rapid binary evolution codes, that is, the duration of core hydrogen-burning, and the resulting donor star mass. We find that systems with shorter orbital periods produce up to 50% lighter stripped donors and have a lifetime up to 30% larger than wider systems. Both quantities strongly depend on the initial binary orbital period, but the initial mass ratio and the mass-transfer efficiency of the binary have little impact on the outcome. Our results are easily parameterisable and can be used to capture the effects of Case A mass transfer more accurately in rapid binary evolution codes. [ABSTRACT FROM AUTHOR]

Context. The hot inner flow in black-hole X-ray binaries is not just a static corona rotating around the black hole: it must be partially outflowing. It is therefore a mildly relativistic "outflowing corona". We have developed a model in which Comptonization takes place in this outflowing corona. In all of our previous work, we assumed a rather high outflow speed of 0.8c. Aims. Here, we investigate whether an outflow with a significantly lower speed can also reproduce the observations. Thus, in this work we consider an outflow speed of 0.1c or less. Methods. As in all of our previous work, we used a Monte Carlo code to compute not only the emergent X-ray spectra, but also the time lags that are introduced to the higher-energy photons with respect to the lower-energy ones via multiple scatterings. We also record the angle (with respect to the symmetry axis of the outflow) and the height at which photons escape. Results. Our results are very similar to those of our previous work, with some small quantitative differences that can be easily explained. We are again able to quantitatively reproduce five observed correlations: (a) the time lag as a function of Fourier frequency, (b) the time lag as a function of photon energy, (c) the time lag as a function of Γ, (d) the time lag as a function of the cutoff energy in the spectrum, and (e) the long-standing radio–X-ray correlation – and all of them with only two parameters, which vary in the same ranges for all the correlations. Conclusions. Our model does not require a compact, narrow relativistic jet, although its presence does not affect the results. The essential ingredient of our model is the parabolic shape of the Comptonizing corona. The outflow speed plays a minor role. Furthermore, the bottom of the outflow, in the hard state, looks like a "slab" to the incoming soft photons from the disk, and this can explain the observed X-ray polarization, which is along the outflow. In the hard-intermediate state, we predict that the polarization of GX 339−4 will be perpendicular to the outflow. [ABSTRACT FROM AUTHOR]

Aims. In this study, we demonstrate the efficacy of the Ultraviolet Imaging Telescope (UVIT) in identifying and characterizing white dwarfs (WDs) within the Milky Way Galaxy. Methods. Leveraging the UVIT point-source catalogue towards the Small Magellanic Cloud and cross-matching it with Gaia DR3 data, we identified 43 single WDs (37 new detections), 13 new WD+main-sequence candidates, and 161 UV bright main-sequence stars by analysing their spectral energy distributions. Using the WD evolutionary models, we determined the masses, effective temperatures, and cooling ages of these identified WDs. Results. The masses of these WDs range from 0.2 to 1.3 M⊙ and the effective temperatures (Teff) lie between 10 000 K to 15 000 K, with cooling ages spanning 0.1–2 Gyr. Notably, we detect WDs that are hotter than reported in the literature, which we attribute to the sensitivity of UVIT. Furthermore, we report the detection of 20 new extremely low-mass candidates from our analysis. Future spectroscopic studies of the extremely low-mass candidates will help us understand the formation scenarios of these exotic objects. Despite limitations in Gaia DR3 distance measurements for optically faint WDs, we provide a crude estimate of the WD space density within 1kpc of 1.3 × 10−3 pc−3, which is higher than previous estimates in the literature. Conclusions. Our results underscore the instrumental capabilities of UVIT and anticipate forthcoming UV missions such as INSIST for systematic WD discovery. Our method sets a precedent for future analyses in other UVIT fields to find more WDs and perform spectroscopic studies to verify their candidacy. [ABSTRACT FROM AUTHOR]

Classical Cepheids provide valuable insights into the evolution of stellar multiplicity among intermediate-mass stars. Here, we present a systematic investigation of single-lined spectroscopic binaries (SB1s) based on high-precision velocities measured by the VELOcities of CEpheids (VELOCE) project. We detected 76 (29%) SB1 systems among the 258 Milky Way Cepheids in the first VELOCE data release, 32 (43%) of which were not previously known to be SB1 systems. We determined 30 precise and three tentative orbital solutions, 18 (53%) of which are reported for the first time. This large set of Cepheid orbits provides a detailed view of the eccentricity e and orbital period Porb distribution among evolved intermediate-mass stars, ranging from e ∈ [0.0, 0.8] and Porb ∈ [240, 9000] d. The orbital motion on timescales exceeding the 11 yr VELOCE baseline was investigated using a template-fitting technique applied to literature data. Particularly interesting objects include (a) R Cru, the Cepheid with the shortest orbital period in the Milky Way (∼238 d); (b) ASAS J103158−5814.7, a short-period overtone Cepheid exhibiting time-dependent pulsation amplitudes as well as orbital motion; and (c) 17 triple systems with outer visual companions, among other interesting objects. Most VELOCE Cepheids (21/23) that exhibit evidence of a companion based on a Gaia proper motion anomaly are also spectroscopic binaries, whereas the remaining do not exhibit significant (> 3σ) orbital radial velocity variations. Gaia quality flags, notably the renormalized unit weight error (RUWE), do not allow Cepheid binaries to be identified reliably although statistically the average RUWE of SB1 Cepheids is slightly higher than that of non-SB1 Cepheids. A comparison with Gaia photometric amplitudes in G-, Bp, and Rp also does not allow one to identify spectroscopic binaries among the full VELOCE sample, indicating that the photometric amplitudes in this wavelength range are not sufficiently informative of companion stars. [ABSTRACT FROM AUTHOR]

In this work we analyzed the long-term X-ray variation of Cir X-1 and its NICER spectra. The dip width and the peak flux after periastron are obtained. Our findings show a negative correlation between the dip width and the peak flux. A new scenario is introduced to explain the correlation. The disk height is higher at a larger radius in the standard model. When the neutron star (NS) moves away from periastron, the contraction of the equipotential surface causes the height of the outer disk to increase significantly. It blocks the X-ray emission near the NS. Subsequently, the height of the outer disk gradually decreases due to the viscosity. The viscous timescale is inversely proportional to the height of the disk, and the height of the outer disk is proportional to the mass transferred from the companion near periastron, thus the peak flux. Consequently, the duration of the obscuring, corresponding to the dip width, is inversely proportional to the peak flux. Moreover, we introduce a new ephemeris MJD(N) = 43075.0 + 16.5843N − 4.778 × 10−5N2 based on the dip-in phases (ingress phase of the dip). [ABSTRACT FROM AUTHOR]

Ultra-light dark matter perturbs the orbital motion of binary pulsars, in particular by causing peculiar time variations of a binary's orbital parameters, which then induce variations in the pulses' times of arrival. Binary pulsars have therefore been shown to be promising detectors of ultra-light dark matter. To date, the sensitivity of binary pulsars to ultra-light dark matter has only been studied for dark matter masses in a narrow resonance band around a multiple of the binary pulsar orbital frequency. In this study we devise a two-step, bayesian method that enables us to compute semi-analytically the sensitivity for all masses, also away from the resonance, and to combine several observed binaries into one global sensitivity curve. We then apply our method to the case of a universal, linearly-coupled, scalar ultra-light dark matter. We find that with next-generation radio observatories the sensitivity to the ultra-light dark matter coupling will surpass that of Solar-System constraints for a decade in mass around m ∼ 10−21 eV, even beyond resonance [ABSTRACT FROM AUTHOR]

Context. Energy and momentum feedback from stars is a key element in models of galaxy formation and interstellar medium (ISM) dynamics, but resolving the relevant length scales in order to directly include this feedback remains beyond the reach of current-generation simulations. Aims. We aim to constrain the energy feedback of winds, photoionisation, and supernovae (SNe) from massive stars. Methods. We measure the thermal and kinetic energy imparted to the ISM on various length scales, which we calculate from high-resolution 1D radiation-hydrodynamics simulations. Our grid of simulations covers a broad range of densities, metallicities, and state-of-the-art evolutionary models of single and binary stars. Results. A single star or binary system can carve a cavity of tens of parsecs (pc) in size into the surrounding medium. During the pre-SN phase, post-main sequence stellar winds and photoionisation dominate. While SN explosions dominate the total energy budget, the pre-SN feedback is of great importance by reducing the circumstellar gas density and delaying the onset of radiative losses in the SN remnant. Contrary to expectations, the metallicity dependence of the stellar wind has little effect on the cumulative energy imparted by feedback to the ISM; the only requirement is the existence of a sufficient level of pre-SN radiative and mechanical feedback. The ambient medium density determines how much and when feedback energy reaches distances of ≳10–20 pc and affects the division between kinetic and thermal feedback. Conclusions. Our results can be used as a subgrid model for feedback in large-scale simulations of galaxies. The results reinforce that the uncertain mapping of stellar evolution sequences to SN explosion energy is very important for determining the overall feedback energy from a stellar population. [ABSTRACT FROM AUTHOR]

Context. With the advent of large spectroscopic surveys, automated stellar parameter determination has become commonplace. Nevertheless, spectral classification still offers a quick and useful alternative for obtaining parameter estimates for large samples of spectra of varying quality. Aims. We present a new atlas of stellar spectra covering the B-type range, with the intention of providing detailed classification criteria valid for modern spectra and improving the grid of reliable standards. This new grid will be used in future works to provide classification criteria beyond the classical classification range and addressing, in particular, the use of Gaia/RVS spectra. Methods. We analysed historical standards by means of multiple high-resolution spectra, marking out problematic cases and complementing the grid with a new set of reliable comparators. We then elaborated on a new set of classification criteria based on high-quality R = 4000 spectra. Results. Our new classification grid is much thicker than any previous set of standards, presenting a high degree of self-consistency. Although it is based entirely on morphological criteria, the grid demonstrates a much better correlation with different physical parameters. Conclusions. The new grid can be used to study classification criteria in other spectral ranges, providing a valuable tool for the study of B-type stars that covers a very wide range of temperatures, luminosities, and stellar masses. The very process of classification also offers valuable insights into stellar evolution. [ABSTRACT FROM AUTHOR]

Context. Wolf–Rayet (WR) stars are massive evolved stars that exhibit particularly fast and dense stellar winds. Although they constitute a very short phase near the end of a massive star's life, they play a crucial role in the evolution of massive stars and have a substantial impact on their surrounding environment. Aims. We present the most comprehensive and deepest X-ray study to date of the properties of the richest Wolf–Rayet population observed in a single stellar cluster, Westerlund 1 (Wd1). By examining the X-ray signatures of WR stars, we aim to shed light on the hottest plasma in their stellar winds and gain insights into whether they exist as single stars or within binary systems. Methods. This work is based on 36 Chandra observations obtained from the "Extended Westerlund 1 and 2 Open Clusters Survey" (EWOCS) project, plus 8 archival Chandra observations. The overall exposure depth Ms) and baseline of the EWOCS observations extending over more than one year enable us to perform a detailed photometric, colour, and spectral analysis, as well as to search for short- and long-term periodicity. Results. In X-rays, we detect 20 out of the 24 known Wolf–Rayet stars in Wd1 down to an observed luminosity of ~7 × 1029 erg s−1 (assuming a distance of 4.23 kpc to Wd1), with 8 WR stars being detected in X-rays for the first time. Nine stars show clear evidence of variability over the year-long baseline, with clear signs of periodicity. The X-ray colours and spectral analysis reveal that the vast majority of the WR stars are hard X-ray sources (kT≥2.0 keV). The Fe XXV emission line at ~6.7 keV, which commonly originates from the wind–wind collision zone in binary systems, is detected for the first time in the spectra of 17 WR stars in Wd1. In addition the ~6.4 keV fluorescent line is observed in the spectra of three stars, which are among the very few massive stars exhibiting this line, indicating that dense cold material coexists with the hot gas in these systems. Overall, our X-ray results alone suggest a very high binary fraction (≥80%) for the WR star population in Wd1. When combining our results with properties of the WR population from other wavelengths, we estimate a binary fraction of ≥92%, which could even reach unity. This suggests that either all the most massive stars are found in binary systems within Wd1, or that binarity is essential for the formation of such a rich population of WR stars. [ABSTRACT FROM AUTHOR]

Context. The existence of low-mass giants with large amounts of lithium (Li) in their surfaces has challenged stellar evolution for decades. One of the possibilities usually discussed in the literature to explain these Li-rich giants involves the interaction with a close binary companion, a scenario that predicts that, when compared against their non-enriched counterparts, Li-rich giants should preferentially be found as part of binary systems. Aims. We aim to assemble the largest possible sample of low-mass giants with well-measured Li abundances, to determine with high statistical significance the close binary fractions of Li-rich and Li-normal giants, and thus test the binary interaction scenario for the emergence of Li-rich giants. Methods. We developed a method that uses radial velocities (RVs) at three different epochs to quantify the degree of RV variability, which we used as a proxy for the presence of a close binary companion. The method was tested and calibrated against samples of known RV standard stars and known spectroscopic binaries. We then assembled a sample of 1418 giants with available RVs from RAVE, GALAH, and Gaia, as well as stellar parameters and Li abundances from GALAH, to which we applied our variability classification. We could determine an evolutionary state for 1030 of these giants. We also compared the results of our RV variability analysis with binarity indicators from the Gaia mission. Results. When applying our methodology to the control samples, we found that the accuracy of the classification is controlled by the precision of the RVs used in the analysis. For the set of RVs available for the giants, this accuracy is 80–85%. Consistent with seismic studies, the resulting sample of giants contains a fraction of Li-rich objects in the red clump (RC) that is twice as large as that in the first ascent red giant branch (RGB). Among RC giants, the fractions of Li-rich objects with a high RV variability and with no RV variability are the same as those for Li-normal objects, but we find some evidence that these fractions may be different for giants in the first-ascent RGB. Analysis of binary indicators in Gaia DR3 shows a smaller fraction of binary giants than our criteria, but no relation can be seen between Li enrichment and binarity either. Conclusions. Our RV variability analysis indicates that there is no preference for Li-rich giants in the RC to be part of binary systems, thus arguing against a binary interaction scenario for the genesis of the bulk of Li-rich giants at that evolutionary stage. On the other hand, Li-rich giants in the RGB appear to have a small but measurable preference for having close companions, something that deserves further scrutiny with more and better data. Additional measurements of the RVs of these giants at a higher RV precision would greatly help in confirming and more robustly quantifying these results. [ABSTRACT FROM AUTHOR]

Context. The transition zone between the white dwarf (WD) envelope and a circumstellar accretion disk in classical novae, the boundary layer, is a region of strong dissipation and intense vorticity. In this strongly sheared layer, the hydrogen-rich accreted gas is expected to mix with the underlying WD outermost layers so the conditions for the onset of the thermonuclear runaway (TNR) in classical nova will be different from the standard treatment of the onset and subsequent mixing. Aims. We applied the critical layer instability (CLI) to the boundary between a disk-accreted H/He zone and the C/O- or O/Ne – rich outer layers of a mass-accreting WD in a cataclysmic binary and then used the resulting structure as input to one-dimensional nuclear-hydrodynamic simulations of the nova outburst. Methods. We simulated the subsonic mixing process in two dimensions for conditions appropriate for the inner disk and a CO 0.8 M⊙ and CO and ONe 1.25 M⊙ WDs using the compressible hydrodynamics code PLUTO. The resulting compositional profile was then imported into the one-dimensional nuclear-hydrodynamics code SHIVA to simulate the triggering and growth rate for the TNR and subsequent envelope ejection. Results. We find that the deep shear driven mixing changes the triggering and development of the TNR. In particular, the time to reach peak temperature is significantly shorter, and the ejected mass and maximum velocity of the ejecta substantially greater, than the current treatment. The 7Li yield is reduced by about an order of magnitude relative to the current treatments. [ABSTRACT FROM AUTHOR]

We present a complete numerical model of the afterglow of a laterally structured relativistic ejecta from the radio to very high energies (VHE). This includes a self-consistent calculation of the synchrotron radiation, with its maximum frequency, and of synchrotron self-Compton (SSC) scattering that takes the Klein-Nishina regime into account. Attenuation due to pair production is also included. This model is computationally efficient and allows multi-wavelength data fitting. As a validation test, the radiative model was used to fit the broad-band spectrum of GRB 190114C at 90 s up to the TeV range. The full model was then used to fit the afterglow of GW 170817 and predict its VHE emission. We find that the SSC flux at the peak was much dimmer than the upper limit from H.E.S.S. observations. However, we show that either a smaller viewing angle or a higher external density would make similar off-axis events detectable in the future at VHE, even above 100 Mpc with the sensitivity of the Cherenkov telescope array. High external densities are expected in the case of fast mergers, but the existence of a formation channel for these binary neutron stars is still uncertain. We highlight that VHE afterglow detections would help to efficiently probe systems like this. [ABSTRACT FROM AUTHOR]

Many circumbinary gas giant planets have been recently discovered. The formation mechanism of circumbinary planets on wide orbits is unclear. We investigate the formation of Delorme 1 (AB)b, a |$13 \pm 5 \ \mathrm{ M_J}$| planet, orbiting its host binary at 84 au. The planet is accreting while having an estimated age of 40 Myr, which is unexpected, as this process should have ceased due to the dissipation of the protoplanetary disc. Using the smoothed particle hydrodynamics code seren , we model three formation scenarios for this planet. In Scenario I, the planet forms in situ on a wide orbit in a massive disc (by gravitational instability), in Scenario II closer to the binary in a massive disc (by gravitational instability), and in Scenario III much closer to the binary in a less massive disc (by core accretion). Planets in Scenario I stay at the observed separation and have mass accretion rates consistent with observed value, but their final mass is too high. In Scenario II, the planet reaches the observed separation through outward migration or scattering by the binary, and has mass accretion rate comparable to the observed; however, the planet mass is above the observed value. In Scenario III, the planet's final mass and mass accretion rate are comparable to the observed ones, but the planet's separation is smaller. We conclude that all models may explain some features of the observations but not all of them, raising questions about how gas is accreted on to the planet from its circumplanetary disc, and the presumed age of the system. [ABSTRACT FROM AUTHOR]

We present the analysis of the microlensing event OGLE-2015-BLG-0845, which was affected by both the microlensing parallax and xallarap effects. The former was detected via the simultaneous observations from the ground and Spitzer , and the latter was caused by the orbital motion of the source star in a relatively close binary. The combination of these two effects led to a mass measurement of the lens object, revealing a low-mass (⁠|$0.14 \pm 0.05 \, \mathrm{ M}_{\odot }$|⁠) M dwarf at the bulge distance (⁠|$7.6 \pm 1.0$| kpc). The source binary consists of a late F-type subgiant and a K-type dwarf of |$\sim 1.2$| and |$\sim 0.9 \mathrm{ M}_{\odot }$|⁠ , respectively, and the orbital period is |$70 \pm 10$| d. OGLE-2015-BLG-0845 is the first single-lens event in which the lens mass is measured via the binarity of the source. Given the abundance of binary systems as potential microlensing sources, the xallarap effect may not be a rare phenomenon. Our work thus highlights the application of the xallarap effect in the mass determination of microlenses, and the same method can be used to identify isolated dark lenses. [ABSTRACT FROM AUTHOR]

We investigate the thermodynamics of close encounters between stellar mass black holes (BHs) in the gaseous discs of active galactic nuclei (AGNs), during which binary black holes (BBHs) may form. We consider a suite of 2D viscous hydrodynamical simulations within a shearing box prescription using the Eulerian grid code athena ++. We study formation scenarios where the fluid is either an isothermal gas or an adiabatic mixture of gas and radiation in local thermal equilibrium. We include the effects of viscous and shock heating, as well as optically thick cooling. We co-evolve the embedded BHs with the gas, keeping track of the energetic dissipation and torquing of the BBH by gas and inertial forces. We find that compared to the isothermal case, the minidiscs formed around each BH are significantly hotter and more diffuse, though BBH formation is still efficient. We observe massive blast waves arising from collisions between the radiative minidiscs during both the initial close encounter and subsequent periapsis periods for successfully bound BBHs. These 'disc novae' have a profound effect, depleting the BBH Hill sphere of gas and injecting energy into the surrounding medium. In analysing the thermal emission from these events, we observe periodic peaks in local luminosity associated with close encounters/periapses, with emission peaking in the optical/near-infrared (IR). In the AGN outskirts, these outbursts can reach 4 per cent of the AGN luminosity in the IR band, with flares rising over 0.5–1 yr. Collisions in different disc regions, or when treated in 3D with magnetism, may produce more prominent flares. [ABSTRACT FROM AUTHOR]

We report the discovery of five white dwarf + ultracool dwarf systems identified as common proper motion wide binaries in the Gaia Catalogue of Nearby Stars. The discoveries include a white dwarf + L subdwarf binary, VVV 1256−62AB, a gravitationally bound system located 75.6 |$^{+1.9}_{-1.8}$|  pc away with a projected separation of 1375 |$^{+35}_{-33}$|  au. The primary is a cool DC white dwarf with a hydrogen dominated atmosphere, and has a total age of |$10.5^{+3.3}_{-2.1}$|  Gyr, based on white dwarf model fitting. The secondary is an L subdwarf with a metallicity of [M/H] = |$-0.72^{+0.08}_{-0.10}$| (i.e. [Fe/H] = |$-0.81\pm 0.10$|⁠) and |$T_{\rm eff}$| = 2298 |$^{+45}_{-43}$|  K based on atmospheric model fitting of its optical to near infrared spectrum, and likely has a mass just above the stellar/substellar boundary. The subsolar metallicity of the L subdwarf and the system's total space velocity of 406 km s−1 indicates membership in the Galactic halo, and it has a flat eccentric Galactic orbit passing within 1 kpc of the centre of the Milky Way every |$\sim$| 0.4 Gyr and extending to 15–31 kpc at apogal. VVV 1256−62B is the first L subdwarf to have a well-constrained age, making it an ideal benchmark of metal-poor ultracool dwarf atmospheres and evolution. [ABSTRACT FROM AUTHOR]

We present a statistical analysis of the ages and metallicities of triple stellar systems that are known to host exoplanets. With controversial cases disregarded, so far 27 of those systems have been identified. Our analysis, based on an exploratory approach, shows that those systems are on average notably younger than stars situated in the solar neighborhood. Though the statistical significance of this result is not fully established, the most plausible explanation is a possible double selection effect due to the relatively high mass of planet-hosting stars of those systems (which spend less time on the main-sequence than low-mass stars) and that planets in triple stellar systems may be long-term orbitally unstable. The stellar metallicities are on average solar-like; however, owing to the limited number of data, this result is not inconsistent with the previous finding that stars with planets tend to be metal-rich as the deduced metallicity distribution is relatively broad. [ABSTRACT FROM AUTHOR]

The Rossiter–McLaughlin effect measures the misalignment between a planet's orbital plane and its host star's rotation plane. Around 10 per cent of planets exhibit misalignments in the approximate range 80°–125°, with their origin remaining a mystery. On the other hand, large misalignments may be common in eccentric circumbinary systems due to misaligned discs undergoing polar alignment. If the binary subsequently merges, a polar circumbinary disc – along with any planets that form within it – may remain inclined near 90 |$^{\circ }$| to the merged star's rotation. To test this hypothesis, we present N -body simulations of the evolution of a polar circumbinary debris disc comprised of test particles around an eccentric binary during a binary merger that is induced by tidal dissipation. After the merger, the disc particles remain on near-polar orbits. Interaction of the binary with the polar-aligned gas disc may be required to bring the binary to the small separations that trigger the merger by tides. Our findings imply that planets forming in discs that are polar-aligned to the orbit of a high-eccentricity binary may, following the merger of the binary, provide a possible origin for the population of near-polar planets around single stars. [ABSTRACT FROM AUTHOR]

Binary star systems are expected to follow Newtonian dynamics similarly to planetary systems. However, reports have been made of wide binary systems with separations around 0.01 pc and larger, showing potential deviations from standard Newtonian motion. This phenomenon, suggestive of the flattening of galactic rotation curves, calls for closer inspection. This study presents an analysis of wide binary stars using data from Gaia Data Release 3 (DR3), a space-based astrometry mission funded by the European Space Agency. The study compares different choices of selection criteria to examine the nature of these apparent anomalous kinematics within the solar neighbourhood. The Gaia data set furnishes detailed parameters, including radial velocity, mass, age, and binary probability in addition to standard parameters. A custom Python tool named bynary facilitated both data processing and analysis. This report reveals that the signs of any anomalous signals systematically diminish as the initial selection criteria are relaxed for degrouping while subsequent filtering remains stringent, leading to the complete elimination of any apparent non-Newtonian motion for binary separations within 0.5 pc. The study shows that any observed anomalous behaviour in solar neighbourhood wide binaries within 130 pc must be produced either by faint companion stars orbiting primary or secondary stars, or by flyby stars. The findings emphasize the importance of the choice of selection criteria in disentangling genuine binary dynamics from external influences. The conclusions align with the predictions of Newtonian mechanics and general relativity, though they do not exclude other phenomena at larger scales. [ABSTRACT FROM AUTHOR]

Over the last couple of years, the appearance of the third data release from the Gaia satellite has triggered various wide binary low acceleration gravity tests. Wide binaries with typical total masses |$\approx 1.0 - 1.6\,\mathrm{ M}_{\odot }$| and separations above a few thousand au probe the low acceleration |$a \lesssim a_{0}$| regime, where at galactic and larger scales gravitational anomalies typically attributed to the presence of an as yet undetected dark matter component appear, where |$a_{0} \approx 1.2\times 10^{-10}$| m s |$^{-2}$| is the acceleration scale of Modified Newtonian Dynamics (MOND). Thus, studies of the relative velocities and separations on the plane of the sky, |$v_{2\mathrm{ D}}$| and |$s_{2\mathrm{ D}}$|⁠ , respectively, of wide binary stars extending to separations above a few kau, provide an independent approach on the empirical study of gravity in the interesting |$a \lesssim a_{0}$| acceleration range. Two independent groups, through complementary approaches, have obtained evidence for a departure from Newtonian predictions in the low acceleration regime, in consistency with MOND expectations for wide binary orbits in the Solar Neighbourhood. Two other groups however, have instead reported results showing a clear preference for Newtonian gravity over various MOND alternatives tested, over the same low acceleration regime. We here take a critical look at the various studies in question, from sample selection to statistical treatment of the wide binary relative velocities obtained. We discover a couple of critical problems in the formal design and statistical implementation shared by the two latter groups, and show explicitly how these yield biased conclusions. [ABSTRACT FROM AUTHOR]

Disc warping, and possibly disc breaking, has been observed in protoplanetary discs around both single and multiple stars. Large warps can break the disc, producing multiple observational signatures. In this work, we use comparisons of disc time-scales to derive updated formulae for disc breaking, with better predictions as to when and where a disc is expected to break and how many breaks could occur. Disc breaking is more likely for discs with small inner cavities, cooler temperatures, and steeper power-law profiles, such that thin, polar-aligning discs are more likely to break. We test our analytical formulae using three-dimensional grid-based simulations of protoplanetary discs warped by the gravitational torque of an inner binary. We reproduce the expected warp behaviours in different viscosity regimes and observe disc breaking at locations in agreement with our derived equations. As our simulations only show disc breaking when disc viscosity is low, we also consider a viscous criterion for disc breaking, where rapid alignment to the precession vector can prevent a break by reducing the maximum misalignment between neighbouring rings. We apply these results to the GW Orionis circumtriple disc and find that the precession induced from the central stars can break the disc if it is relatively thin. We expect repeated or multiple disc breaking to occur for discs with sufficiently steep power-law profiles. We simulate a polar-aligning disc around an eccentric binary with steep power-law profiles and observe two separate breaking events at locations in rough agreement with our analytical predictions. [ABSTRACT FROM AUTHOR]

Radiation-driven outflows play a crucial role in extracting mass and angular momentum from binary systems undergoing rapid mass transfer at super-Eddington rates. To study the mass transfer process from a massive donor star to a stellar-mass black hole (BH), we perform multidimensional radiation-hydrodynamical simulations that follow accretion flows from the first Lagrange point down to about a hundred times the Schwarzschild radius of the accreting BH. Our simulations reveal that rapid mass transfer occurring at over a thousand times the Eddington rate leads to significant mass-loss from the accretion disc via radiation-driven outflows. Consequently, the inflow rates at the innermost radius are regulated by two orders of magnitude smaller than the transfer rates. We find that convective motions within the accretion disc drive outward energy and momentum transport, enhancing the radiation pressure in the outskirts of the disc and ultimately generating large-scale outflows with sufficient energy to leave the binary. Furthermore, we observe strong anisotropy in the outflows, which occur preferentially toward both the closest and furthest points from the donor star. However, when averaged over all directions, the specific angular momentum of the outflows is nearly comparable to the value predicted in the isotropic emission case. Based on our simulation results, we propose a formula that quantifies the mass growth rates on BHs and the mass-loss rates from binaries due to radiation-driven outflows. This formula provides important implications for the binary evolution and the formation of merging binary BHs. [ABSTRACT FROM AUTHOR]

How massive stars end their lives depends on the core mass, core angular momentum, and hydrogen envelopes at death. However, these key physical facets of stellar evolution can be severely affected by binary interactions. In turn, the effectiveness of binary interactions itself varies greatly depending on the initial conditions of the binaries, making the situation much more complex. We investigate systematically how binary interactions influence core–collapse progenitors and their fates. Binary evolution simulations are performed to survey the parameter space of supernova progenitors in solar metallicity binary systems and to delineate major evolutionary paths. We first study fixed binary mass ratios (⁠|$q=M_2/M_1$| = 0.5, 0.7, and 0.9) to elucidate the impacts of initial mass and initial separation on the outcomes, treating separately Type Ibc supernova, Type II supernova, accretion-induced collapse (AIC), rapidly rotating supernova (Ibc-R), black hole formation, and long gamma ray burst (long GRB). We then conduct 12 binary population synthesis model calculations, varying the initial condition distributions and binary evolution parameters, to estimate various supernova fractions. We obtain a Milky Way supernova rate |$R_{\rm SN} = (1.78$| – |$2.47) \times 10^{-2} \, {\rm yr}^{-1}$| which is consistent with observations. We find the rates of AIC, Ibc-R, and long GRB to be |$\sim 1/100$| the rate of regular supernovae. Our estimated long GRB rates are higher than the observed long GRB rate and close to the low luminosity GRB rate, although care must be taken considering our models are computed with solar metallicity. Furthering binary modelling and improving the inputs one by one will enable more detailed studies of these and other transients associated with massive stars. [ABSTRACT FROM AUTHOR]

Many of the most intriguing features, including spirals and cavities, in the current disc observations are found in binary systems like GG Tau, HD 142527, or HD 100453. Such features are evidence of the dynamic interactions between binary stars and the viscous disc. Understanding these dynamic interactions and how they result in the structure and growth of asymmetric circumbinary discs is a difficult problem, for which there is no complete analytical solution, that predicts the shape of the observed disc accurately. We use numerical simulations to evolve circumbinary discs with varying disc viscosities and investigate the size and shape of the inner cavities in such protoplanetary discs. We have simulated over 140 locally isothermal 2D grid-based disc models for |$\ge$| 3 × 104 binary orbits each and mapped out the parameter space relevant for protoplanetary discs. With this, it becomes possible to create parametrized profiles for individual discs to compare to observations and find limits to their binary eccentricity or internal viscosity from the simulation data. In the long-term simulations, larger cavity sizes than previously considered are possible within the parameter space (⁠|$\le$| 6 binary separations). As an example, we find that the eccentricity of the disc around HD 142527 suggests the impact of the binary dynamics on the disc. However, even considering the larger cavity sizes, the large size of the cavity in HD 142527 remains unexplained by the simulations considering the most recent orbital constraints. [ABSTRACT FROM AUTHOR]

We study the wide-binary eccentricity (e) distribution in young star clusters and the role of turbulence in setting the form of the e distribution using magnetohydrodynamical simulations of star cluster formation. The simulations incorporate gravity, turbulence, magnetic fields, protostellar heating, and jets/outflows. We find that (1) simulations that employ purely compressive turbulence driving produce binaries with a superthermal e distribution [ |$\alpha \gt 1$| in |$p(e) \propto e^\alpha$| ], while simulations with purely solenoidal driving or natural mixture of driving modes produce subthermal/thermal distributions (⁠|$\alpha \le$| 1), (2) the e distribution over the full range of binary separations in our simulations is set at the early stages of the star cluster formation process, (3) while binaries (separation of |$r_{\mathrm{pair}} \le 1000\, \mathrm{AU}$|⁠) have subthermal to thermal e distributions (⁠|$\alpha \sim 0.8$|⁠), wide binaries (⁠|$r_{\mathrm{pair}} \gt 1000\, \mathrm{AU}$|⁠) have a superthermal distribution (⁠|$\alpha \sim 1.8$|⁠), and (4) low-mass binary systems (system masses of |$M_{\mathrm{sys}} \le 0.8\, \mathrm{M_\odot }$|⁠) have a highly superthermal distribution (⁠|$\alpha \sim 2.4$|⁠), whereas high-mass systems (⁠|$M_{\mathrm{sys}} \gt 0.8\, \mathrm{M_\odot }$|⁠) exhibit a subthermal/thermal distribution (⁠|$\alpha \sim 0.8$|⁠). The binary eccentricity distribution is often modelled as a thermal distribution. However, our results suggest that the e distribution depends on the range of separation of the sampled binaries, which agrees with the findings from recent Gaia observations. We conclude that the dependence of the e distribution on the binary separation and mass is linked to the binary formation mechanism governed by the turbulent properties of the parent cloud. [ABSTRACT FROM AUTHOR]

We report results from 8 h of JWST /MIRI low resolution spectroscopic (LRS) monitoring directly followed by 7 h of JWST /NIRSpec prism spectroscopic monitoring of the benchmark binary brown dwarf WISE 1049AB, the closest, brightest brown dwarfs known. We find water, methane, and CO absorption features in both components, including the 3.3  μ m methane absorption feature and a tentative detection of small grain (⁠|$\lt $| 1 μ m) silicate absorption at |$\gt $| 8.5 μ m in WISE 1049A. Both components vary significantly (⁠|$\gt 1~{{\rm per\ cent}}$|⁠), with WISE 1049B displaying larger variations than WISE 1049A. Using K -means clustering, we find three main transition points in wavelength for both components of the binary: (1) change in behaviour at |$\sim$| 2.3 μ m coincident with a CO absorption bandhead, (2) change in behaviour at 4.2 μ m, close to the CO fundamental band at |$\lambda \gt $|  4.4 µm, and (3) change in behaviour at 8.3–8.5 µm, potentially corresponding to silicate absorption. We interpret the light curves observed with both NIRSpec and MIRI as likely stemming from (1) a deep pressure level driving the double-peaked variability seen in WISE 1049B at wavelengths |$\lt $| 2.3 and |$\gt $| 8.5 µm, (2) an intermediate pressure level shaping the light-curve morphology between 2.3 and 4.2 µm, and (3) a higher altitude pressure level producing single-peaked and plateaued light-curve behaviour between 4.2 and 8.5 µm. [ABSTRACT FROM AUTHOR]

General-relativistic equilibria of differentially rotating stars are expected in a number of astrophysical scenarios, from core-collapse supernovae to the remnant of binary neutron-star mergers. The latter, in particular, have been the subject of extensive studies where they were modelled with a variety of laws of differential rotation with varying degree of realism. Starting from accurate and fully general-relativistic simulations of binary neutron-star mergers with various equations of state and mass ratios, we establish the time when the merger remnant has reached a quasi-stationary equilibrium and extract in this way realistic profiles of differential rotation. This allows us to explore how well traditional laws reproduce such differential-rotation properties and to derive new laws of differential rotation that better match the numerical data in the low-density Keplerian regions of the remnant. In this way, we have obtained a novel and somewhat surprising result: the dynamical stability line to quasi-radial oscillations computed from the turning-point criterion can have a slope that is not necessarily negative with respect to the central rest-mass density, as previously found with traditional differential-rotation laws. Indeed, for stellar models reproducing well the properties of the merger remnants, the slope is actually positive, thus reflecting remnants with angular momentum at large distances from the rotation axis, and hence with cores having higher central rest-mass densities and slower rotation rates. [ABSTRACT FROM AUTHOR]