Search

A neutron star (NS) accreting matter from a companion star in a low-mass X-ray binary (LMXB) system can spin up to become a millisecond pulsar (MSP). Properties of many such MSP systems are known, which is excellent for probing fundamental aspects of NS physics when modelled using the theoretical computation of NS LMXB evolution. Here, we systematically compute the long-term evolution of NS, binary, and companion parameters for NS LMXBs using the stellar evolution code mesa. We consider the baryonic to gravitational mass conversion to calculate the NS mass evolution and show its cruciality for the realistic computation of some parameters. With computations using many combinations of parameter values, we find the general nature of the complex NS spin frequency (⁠|$\nu$|⁠) evolution, which depends on various parameters, including accretion rate, fractional mass-loss from the system, and companion star magnetic braking. Further, we utilize our results to precisely match some main observed parameters, such as |$\nu$|⁠ , orbital period (⁠|$P_{\rm orb}$|⁠), etc. of four accreting millisecond X-ray pulsars (AMXPs). By providing the |$\nu$|⁠ , |$P_{\rm orb}$|⁠ , and the companion mass spaces for NS LMXB evolution, we indicate the distribution and plausible evolution of a few other AMXPs. We also discuss the current challenges in explaining the parameters of AMXP sources with brown dwarf companions and indicate the importance of modelling the transient accretion in LMXBs as a possible solution. [ABSTRACT FROM AUTHOR]

Recently inaugurated telescopes, such as the MeerKAT radio telescope and the upcoming Rubin Observatory Legacy Survey of Space and Time, will be able to detect millions of transient events in the night sky. Stellar collisions between white dwarfs (WDs) and main-sequence (MS) stars may be detectable among such transients. Simulations will play a key role in characterizing these events and selecting targets for follow-up. We present 3D smoothed particle hydrodynamics models of dynamical interactions between a |$0.6\ {\rm M}_{\odot }$|  WD and 0.3, 0.6, and |$1.2\ {\rm M}_{\odot }$|  MS stars within globular cluster environments. Utilizing a 34-isotope nuclear network, we investigate the energetics, gas morphologies, and mass-loss properties of these collisions for different stellar mass ratios. Our models predict an overabundance of |$^{13}$| C, |$^{15}$| N, and |$^{17}$| O isotopes relative to solar abundances. Moreover, we find that the time-scale of the collisions is too short and maximum temperatures too low for any significant hydrogen burning or triple-alpha reactions to occur. This combined with a negligible production of elements heavier than neon may be key signatures in distinguishing these events from other transient events with similar peak bolometric luminosities (⁠|$\sim\!\! 10^{38}\!\!-\!\!10^{41}$|  erg s |$^{-1}$|⁠). [ABSTRACT FROM AUTHOR]

Star formation is a multiscale problem, and only global simulations that account for the connection from the molecular cloud-scale gas flow to the accreting protostar can reflect the observed complexity of protostellar systems. Star-forming regions are characterized by supersonic turbulence, and as a result, it is not possible to simultaneously design models that account for the larger environment and in detail reproduce observed stellar systems. Instead, the stellar inventories can be matched statistically, and the best matches found that approximate specific observations. Observationally, a combination of single-dish telescopes and interferometers are now able to resolve the nearest protostellar objects on all scales from the protostellar core to the inner |$10\, \mathrm{au}$|⁠. We present a new non-parametric methodology which uses high-resolution simulations and post-processing methods to match simulations and observations using deep learning. Our goal is to perform a down-selection from large data sets of synthetic images to a ranked list of best-matching candidates with respect to the observation. This is particularly useful for binary and multiple stellar systems that form in turbulent environments. The objective is to accelerate the rate at which we can do such comparisons, remove biases from hand-picking matches, and contribute to identifying the underlying physical processes that drive the creation and evolution of observed protostellar systems. [ABSTRACT FROM AUTHOR]

Protostellar binaries harbour complex environment morphologies. Observations represent a snapshot in time, and projection and optical depth effects impair our ability to interpret them. Careful comparison with high-resolution models that include the larger star-forming region can help isolate the driving physical processes and give context in the time domain to the observations. We carry out four zoom-in simulations with au scale resolution that result in three binaries and a single star. For the first time ever, we follow the detailed evolution of a protobinary in a full molecular cloud context until a circumbinary disc forms. We investigate the gas dynamics around the young stars and extract disc sizes. Using radiative transfer, we obtain the evolutionary tracer T bol of the binary systems. We find that the centrifugal radius in prestellar cores is a poor estimator of the resulting disc size due to angular momentum transport at all scales. For binaries, the disc sizes are regulated periodically by the binary orbit, having larger radii close to the apastron. The bolometric temperature differs systematically between edge-on and face-on views and shows a high-frequency time dependence correlated with the binary orbit and a low-frequency time dependence with larger episodic accretion events. These oscillations can cause the appearance of the system to change rapidly from class 0 to class I and, for short periods, even bring it to class II. The highly complex structure in early stages, as well as the binary orbit itself, affects the classical interpretation of protostellar classes, and the direct translation to evolutionary stages has to be done with caution and include other evolutionary indicators such as the extent of envelope material. [ABSTRACT FROM AUTHOR]

Context. 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. Aims. 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. Methods. 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. Results. The findings demonstrate a high agreement level (∼70–90%) 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% reproducibility for binary targets using photometry, and CNN reaches 84.94% 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 ∼700 nm. Conclusions. 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. [ABSTRACT FROM AUTHOR]

Context. Stars and planets form in regions of enhanced stellar density, subjecting protoplanetary discs to gravitational perturbations from neighbouring stars. Observations in the Taurus star-forming region have uncovered evidence of at least three recent, star-disc encounters that have truncated discs (HV/DO Tau, RW Aurigae, and UX Tau), raising questions about the frequency of such events. Aims. We aim to assess the probability of observing truncating star-disc encounters in Taurus. Methods. We generated a physically motivated dynamical model including binaries and a spatial-kinematic substructure to follow the historical dynamical evolution of the Taurus star-forming region. We used this model to track star-disc encounters and the resulting outer disc truncation over the lifetime of Taurus. Results. A quarter of discs are truncated below 30 au by dynamical encounters, but this truncation mostly occurs in binaries over the course of a few orbital periods, on a timescale ≲0.1 Myr. Nonetheless, some truncating encounters still occur up to the present age of Taurus. Strongly truncating encounters (ejecting ≳10 percent of the disc mass) occur at a rate ∼10 Myr−1, sufficient to explain the encounter between HV and DO Tau ∼0.1 Myr ago. If encounters that eject only ∼1 percent of the disc mass are responsible for RW Aurigae and UX Tau, then they are also expected with encounter rate Γenc ∼ 100–200 Myr−1. However, the observed sample of recent encounters is probably incomplete, since these examples occurred in systems that are not consistent with a random drawing from the mass function. One more observed example would statistically imply additional physics, such as replenishment of the outer disc material. Conclusions. The marginal consistency of the frequency of observed recent star-disc encounters with theoretical expectations underlines the value of future large surveys searching for external structures associated with recent encounters. The outcome of such a survey offers a highly constraining, novel probe of protoplanetary disc physics. [ABSTRACT FROM AUTHOR]

We present the first X-ray observation of the energetic millisecond pulsar binary PSR J1431−4715, performed with XMM-Newton and complemented with fast optical multi-band photometry acquired with the ULTRACAM instrument at ESO-NTT. It is found as a faint X-ray source without a significant orbital modulation. This contrasts with the majority of systems that instead display substantial X-ray orbital variability. The X-ray spectrum is dominated by non-thermal emission and, due to the lack of orbital modulation, does not favour an origin in an intrabinary shock between the pulsar and companion star wind. While thermal emission from the neutron star polar cap cannot be excluded in the soft X-rays, the dominance of synchrotron emission favours an origin in the pulsar magnetosphere that we describe at both X-ray and gamma-ray energies with a synchro-curvature model. The optical multi-colour light curve folded at the 10.8 h orbital period is double-humped and dominated by ellipsoidal effects, but also affected by irradiation. The ULTRACAM light curves are fit with several models encompassing direct heating and a cold spot, or heat redistribution after irradiation either through convection or convection plus diffusion. Despite the inability to constrain the best irradiation models, the fits provide consistent system parameters, giving an orbital inclination of 59 ± 6° and a distance of 3.1 ± 0.3 kpc. The companion is found to be an F-type star, underfilling its Roche lobe (fRL = 73 ± 4%) with a mass of 0.20 ± 0.04 M⊙, confirming the redback status, but hotter than the majority of redbacks. The stellar dayside and nightside temperatures of 7500 K and 7400 K, respectively, indicate a weak irradiation effect on the companion, likely due to its high intrinsic luminosity. Although the pulsar mass cannot be precisely derived, a heavy (1.8−2.2 M⊙) neutron star is favoured. [ABSTRACT FROM AUTHOR]

Context. To date, puzzlingly few bona fide RR Lyrae stars have been identified in binaries. Binarity in pulsating stars can be revealed through well-timed photometric data over sufficiently long time bases because of the light travel time effect (LTTE) on the pulsations, which manifests as variation in the timings of maximum light in the O−C (observed minus calculated) diagram. However, O−C variations can also have other causes, such as the Blazhko effect or a sudden or gradual change in the main pulsation period. Aims. We approach this challenge by disentangling the Blazhko effect and period changes from the potential LTTE on V1109 Cas, an RR Lyrae star suspected to be part of a binary system based on its O−C data in the GEOS database. In doing so, we aim to uncover the subtler signals indicative of a companion (LTTE). Methods. We analysed nine years of ground-based photometric data, using Fourier analysis to model the pulsation modulated by the Blazhko effect. From the fit to the observed light curves, we constructed a refined O−C diagram without the scatter caused by Blazhko modulation. Subsequently, we considered different possible scenarios, because distinguishing intrinsic period changes or breaks from the LTTE is challenging. Results. If the remaining O−C variation is due to a period break, refining the O−C diagram can almost entirely remove the trends. If we interpret the variation as an LTTE, we can consider possible configurations. While we currently favour a period break as the explanation for the residual O−C variations, more data on the star in the forthcoming years will help reveal whether the O−C variations are due to the LTTE or intrinsic period changes. Conclusions. Our study, based on a detailed light-curve fitting of V1109 Cas, offers insights into the discernment between Blazhko modulation, intrinsic period breaks, and the LTTE in binary systems. Our results highlight the complexity of determining binarity from O−C data traditionally used for detecting binary motion, calling for caution in this process. [ABSTRACT FROM AUTHOR]

Context. High-resolution observations reveal that the outflows of evolved low- and intermediate-mass stars harbour complex morphological structures that are linked to the presence of one or multiple companions. Hydrodynamical simulations provide a way to study the impact of a companion on the shaping of the asymptotic giant branch (AGB) star out-flow. Aims. Using smoothed particle hydrodynamic (SPH) simulations of an AGB star undergoing mass loss, which also has a binary companion, we study the impact of including H I atomic line cooling on the flow morphology. We also study how this affects the properties of the accretion disks that form around the companion. Methods. We used the PHANTOM code to perform high-resolution 3D SPH simulations of the interaction of a solar-mass companion with the outflow of an AGB star, using different wind velocities and eccentricities. We compared the model properties, computed with and without the inclusion of H I cooling. Results. The inclusion of H I cooling produces a sizeable decrease in the temperature, up to one order of magnitude, in the region closely surrounding the companion star. As a consequence, the morphological irregularities and relatively energetic (bipolar) outflows that were obtained without H I cooling no longer appear. In the case of an eccentric orbit and a low wind velocity, these morphologies are still highly asymmetric, but the same structures recur at every orbital period, making the morphology more regular. Flared accretion disks, with a (sub-)Keplerian velocity profile, are found to form around the companion in all our models with H I cooling, provided the accretion radius is small enough. The disks have radial sizes ranging from about 0.4 to 0.9 au and masses around 10−7−10−8 M⊙. For the considered wind velocities, mass accretion onto the companion is up to a factor of 2 higher than predicted by the standard Bondi Hoyle Littleton rate, ranging between ~4 to 21% of the AGB wind mass loss rate. The lower the wind velocity at the location of the companion, the larger and the more massive the disk and the higher the mass accretion efficiency. In eccentric systems, the disk size, disk mass, and mass accretion efficiency vary, depending on the orbital phase. Conclusions. H I cooling is an essential ingredient to properly model the medium around the companion where density-enhanced wind structures form and it favours the formation of an accretion disk. [ABSTRACT FROM AUTHOR]

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 L2 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%. [ABSTRACT FROM AUTHOR]

Context. Unresolved binaries have a strong influence on the observed parameters of stellar clusters (SCs). Aims. We quantify this influence and compute the resulting mass underestimates and stellar mass function (MF). Methods. N-body simulations of realistic SCs were used to investigate the evolution of the binary population in a SC and its tidal tails. Together with an empirically gauged stellar mass-luminosity relation, the results were then used to determine how the presence of binaries changes the photometric mass and MF of the SC and its tails as deduced from observations. Results. Tail 1 (T1), which is the tidal tail caused by gas expulsion, contains a larger fraction of binaries than both the SC and Tail 2 (T2), which forms after gas expulsion. Additionally, T1 has a larger velocity dispersion. Using the luminosity of an unresolved binary, an observer would underestimate its mass. This bias sensitively depends on the companion masses due to the structure of the stellar mass-luminosity relation. Combining the effect of all binaries in the simulation, the total photometric mass of the SC is underestimated by 15%. Dark objects (black holes and neutron stars) increase the difference between the real and observed mass of the SC further. For both the SC and the tails, the observed power-law index of the MF between a stellar mass of 0.3 and 0.7 M⊙ is smaller by up to 0.2 than the real one, the real initial mass function (IMF) being steeper by this amount. This difference is larger for stars with a larger velocity dispersion or binary fraction. Conclusions. Since the stars formed in SCs are the progenitors of the Galactic field stars, this work suggests that the binary fractions of different populations of stars in the Galactic disc will differ as a function of the velocity dispersion. However, the direction of this correlation is currently unclear, and a complete population synthesis will be needed to investigate this effect. Variations in the binary fractions of different clusters can lead to perceived variations of the deduced stellar MFs. [ABSTRACT FROM AUTHOR]

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

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]