Search

Supernova explosion and the associated neutron star natal kicks are important events on a pathway of a binary to become a gravitational wave source, an X-ray binary or a millisecond radio pulsar. Weak natal kicks often lead to binary survival, while strong kicks frequently disrupt the binary. In this article, we aim to further constrain neutron star natal kicks in binaries. We explore binary population synthesis models by varying prescription for natal kick, remnant mass and mass accretion efficiency. We introduce a robust statistical technique to analyse combined observations of different nature. Using this technique, we further test different models using parallax and proper motion measurements for young isolated radio pulsars and similar measurements for Galactic Be X-ray binaries. Our best model for natal kicks is consistent with both measurements and contains a fraction of $w=0.2\pm 0.1$ weak natal kicks with $\sigma_1 = 45^{+25}_{-15}$ km/s, the remaining natal kicks are drawn from the high-velocity component, same as in previous works: $\sigma_2 = 336$ km/s. We found that currently used models for natal kicks of neutron stars produced by electron capture supernova (combination of maxwellian $\sigma=265$ km/s and $\sigma = 30$ km/s for electron capture) are inconsistent or marginally consistent with parallaxes and proper motions measured for isolated radio pulsars. We suggest a new model for natal kicks of ecSN, which satisfy both observations of isolated radio pulsars and Be X-ray binaries., Comment: 21 page, 36 figures, accepted to MNRAS on September 21

Observations found that star clusters contain a large fraction of binaries. Tight binaries are an important heating source that influences the long-term dynamical evolution of star clusters. However, due to the limitation of $N$-body tool, previous theoretical modelling for globular clusters (GCs) by using direct $N$-body simulations have not investigated how a large fraction of primordial binaries affect their long-term evolution. In this work, by using the high-performance $N$-body code, PeTar, we carry out star-by-star models for intermediate massive GCs ($N=100000$) with the primordial binary fraction varying from 0 to 1. We find that when a stellar-mass black hole (BH) subsystem exists, the structural evolution of GCs (core and half-mass radii) only depends on the properties of massive primordial binaries, because they affect the number of BH binaries (BBHs), which dominate the binary heating process. Low-mass binaries including double white dwarf binaries (BWDs) have almost no influence on the dynamics. Meanwhile, only gravitational wave (GW) mergers from BBHs are strongly affected by dynamical interactions, while low-mass mergers from BWDs show no difference in the isolated environment (field) and in GCs. Low-mass binaries become important only after most BHs escape and the core collapse of light stars occurs. Our result suggests that for $N$-body modelling of GCs with a black hole subsystem dominating binary heating, it is not necessary to include low-mass binaries. These binaries can be studied separately by using standalone binary stellar evolution codes. This way can significantly reduce the computing cost., Comment: 10 pages, 6 figures, submitted to MNRAS

We present results from the B-type Binaries Characterisation (BBC) programme, a multi-epoch spectroscopic study of 88 early B-type binary candidates in the 30 Doradus region of the Large Magellanic Cloud (LMC). From radial-velocity analysis of 29 observational epochs we confirm the binary status of 64 of our targets, comprising 50 SB1 and 14 SB2 B-type binaries. A further 20 systems (classified as SB1*) show clear signs of periodicity but with more tentative periods. Orbital solutions are presented for these 84 systems, providing the largest homogeneous sample to date of the binary properties of early B-type stars. Our derived orbital-period distribution is generally similar to those for samples of more massive (O-type) binaries in both the LMC and the Galaxy. This similarity with the properties of the more massive O-type binaries is important as early B-type stars are expected to account for the majority of core-collapse supernovae. Differences in the period distributions of the different samples start to increase above 4 d, and are also present between the earliest (B0-0.7) and later-type (B1-2.5) systems within the BBC sample, although further study is required to understand if this is an observational bias or a real physical effect. We have examined the semi-amplitude velocities and orbital periods of our sample to identify potential candidates that could hide compact companions. Comparing with probability distributions of finding black hole companions to OB-type stars from a recent theoretical study, we have found 16 binaries in the higher probability region that warrant further study., 42 pages (main body: 28 pages, supplementary material: 14 pages), 20 figures. Accepted for publication in MNRAS

We study gravitational lensing of gravitational waves from compact object binaries as a probe of compact dark matter (DM) objects such as primordial black holes. Assuming a point mass lens, we perform parameter estimation of lensed gravitational wave signals from compact object binaries to determine the detectability of the lens with ground based laser interferometers. Then, considering binary populations that LIGO-Virgo has been probing, we derive a constraint on the abundance of compact DM from non-observation of lensed events. We find that the LIGO-Virgo observations imply that compact objects heavier than $M_l = 200M_\odot$ can not constitute all DM and less than $40\%$ of DM can be in compact objects heavier than $M_l = 400M_\odot$. We also show that the DM fraction in compact objects can be probed by LIGO in its final sensitivity for $M_l > 40M_\odot$ reaching $2\%$ of the DM abundance at $M_l > 200M_\odot$, and by ET for $M_l > 1M_\odot$ reaching DM fraction as low as $7\times 10^{-5}$ at $M_l > 40M_\odot$., 8 pages, 4 figures. published version

As gravitational-wave (GW) interferometers become more sensitive and probe ever more distant reaches, the number of detected binary neutron star mergers will increase. However, detecting more events farther away with GWs does not guarantee corresponding increase in the number of electromagnetic counterparts of these events. Current and upcoming wide-field surveys that participate in GW follow-up operations will have to contend with distinguishing the kilonova from the ever increasing number of transients they detect, many of which will be consistent with the GW sky-localization. We have developed a novel tool based on a temporal convolutional neural network architecture, trained on sparse early-time photometry and contextual information for Electromagnetic Counterpart Identification (El-CID). The overarching goal for El-CID is to slice through list of new transient candidates that are consistent with the GW sky localization, and determine which sources are consistent with kilonovae, allowing limited target-of-opportunity resources to be used judiciously. In addition to verifying the performance of our algorithm on an extensive testing sample, we validate it on AT2017gfo - the only EM counterpart of a binary neutron star merger discovered to date - and AT2019npv - a supernova that was initially suspected as a counterpart of the gravitational-wave event, GW190814, but was later ruled out after further analysis., Comment: 17 pages, 10 figures, accepted in MNRAS

Many young and intermediate age massive stellar clusters host bimodal distributions in the rotation rates of their stellar populations, with a dominant peak of rapidly rotating stars and a secondary peak of slow rotators. The origin of this bimodal rotational distribution is currently debated and two main theories have been put forward in the literature. The first is that all/most stars are born as rapid rotators and that interacting binaries brake a fraction of the stars, resulting in two populations. The second is that the rotational distribution is a reflection of the early evolution of pre-main sequence stars, in particular, whether they are able to retain or lose their protoplanetary discs during the first few Myr. Here, we test the binary channel by exploiting multi-epoch VLT/MUSE observations of NGC 1850, a 100Myr massive cluster in the LMC, to search for differences in the binary fractions of the slow and fast rotating populations. If binarity is the cause of the rotational bimodality, we would expect that the slowly rotating population should have a much larger binary fraction than the rapid rotators. However, in our data we detect similar fractions of binary stars in the slow and rapidly rotating populations (5.9+/-1.1% and 4.5+/-0.6%, respectively).Hence, we conclude that binarity is not a dominant mechanism in the formation of the observed bimodal rotational distributions., Comment: 5 pages, 2 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

Some surviving companions of type Ia supernovae (SNe Ia) from the white-dwarf + main-sequence (WD + MS) channel may evolve to hot subdwarfs. In this paper, we preformed stellar evolution calculations for surviving companions of close WD + MS systems in the spin-up/spin-down model and the canonical non-rotating model to map out the initial parameter spaces in the orbital period - secondary mass plane in which the surviving companions can evolve to hot subdwarfs. Based on these results, we carried out a series of binary population synthesis calculation to obtain the Galactic birth rate of the hot subdwarfs from the WD + MS channel, which is $2.3-6\times10^{\rm -4}\,{\rm yr}^{\rm -1}$ for the spin-up/spin-down model and $0.7-3\times10^{\rm -4}\,{\rm yr}^{\rm -1}$ for the canonical non-rotating model. We also show the distributions of some integral properties of the hot subdwarfs, e.g. the mass and the space velocity, for different models. In addition, comparing our results with the observations of the intermediate helium-rich (iHe-rich) hot subdwarfs, the hot subdwarfs from the WD + MS channel may explain some observational features of the iHe-rich hot subdwarfs, especially for those from the spin-up/spin-down model. Although we expect that the SN Ia channel can only contribute a small fraction of the iHe-rich hot subdwarf population, some of these may help to explain cases with unusual kinematics., Comment: 20 figures. Published by MNRAS

We present the discovery and characterisation of an eclipsing binary identified by the Next Generation Transit Survey in the $\sim$115 Myr old Blanco 1 open cluster. NGTS J0002-29 comprises three M dwarfs: a short-period binary and a companion in a wider orbit. This system is the first well-characterised, low-mass eclipsing binary in Blanco 1. With a low mass ratio, a tertiary companion and binary components that straddle the fully convective boundary, it is an important benchmark system, and one of only two well-characterised, low-mass eclipsing binaries at this age. We simultaneously model light curves from NGTS, TESS, SPECULOOS and SAAO, radial velocities from VLT/UVES and Keck/HIRES, and the system's spectral energy distribution. We find that the binary components travel on circular orbits around their common centre of mass in $P_{\rm orb} = 1.09800524 \pm 0.00000038$ days, and have masses $M_{\rm pri}=0.3978\pm 0.0033$ M$_{\odot}$ and $M_{\rm sec}=0.2245\pm 0.0018$ M$_{\odot}$, radii $R_{\rm pri}=0.4037\pm 0.0048$ R$_{\odot}$ and $R_{\rm sec}=0.2759\pm 0.0055$ R$_{\odot}$, and effective temperatures $T_{\rm pri}=3372\,^{+44}_{-37}$ K and $T_{\rm sec}=3231\,^{+38}_{-31}$ K. We compare these properties to the predictions of seven stellar evolution models, which typically imply an inflated primary. The system joins a list of 19 well-characterised, low-mass, sub-Gyr, stellar-mass eclipsing binaries, which constitute some of the strongest observational tests of stellar evolution theory at low masses and young ages., 21 pages, 9 figures. Accepted for publication in MNRAS

Accreting magnetic white dwarfs offer an opportunity to understand the interplay between spin-up and spin-down torques in binary systems. Monitoring of the white dwarf spin may reveal whether the white dwarf spin is currently in a state of near-equilibrium, or of uni-directional evolution towards longer or shorter periods, reflecting the recent history of the system and providing constraints for evolutionary models. This makes the monitoring of the spin history of magnetic white dwarfs of high interest. In this paper we report the results of a campaign of follow-up optical photometry to detect and track the 39 sec white dwarf spin pulses recently discovered in Hubble Space Telescope data of the cataclysmic variable V1460 Her. We find the spin pulsations to be present in g-band photometry at a typical amplitude of 0.4%. Under favourable observing conditions, the spin signal is detectable using 2-meter class telescopes. We measured pulse-arrival times for all our observations, which allowed us to derive a precise ephemeris for the white dwarf spin. We have also derived an orbital modulation correction that can be applied to the measurements. With our limited baseline of just over four years, we detect no evidence yet for spin-up or spin-down of the white dwarf, obtaining a lower limit of |P/Pdot|> 4e7 years, which is already 4 to 8 times longer than the timescales measured in two other cataclysmic variable systems containing rapidly rotating white dwarfs, AE Aqr and AR Sco., 8 pages, 6 figures. Accepted for publication in MNRAS

We investigate the low-luminosity supernova SN 2016bkv and its peculiar early-time interaction. For that, we compute radiative transfer models using the CMFGEN code. Because SN 2016bkv shows signs of interaction with material expelled by its progenitor, it offers a great opportunity to constrain the uncertain evolutionary channels leading to low-luminosity supernovae. Our models indicate that the progenitor had a mass-loss rate of (6.0 +- 2.0) x 1e-4 Msun/yr (assuming a velocity of 150 km/s). The surface abundances of the progenitor are consistent with solar contents of He and CNO. If SN 2016bkv's progenitor evolved as a single star, it was an odd red supergiant that did not undergo the expected dredge up for some reason. We propose that the progenitor more likely evolved through binary interaction. One possibility is that the primary star accreted unprocessed material from a companion and avoided further rotational and convective mixing until the SN explosion. Another possibility is a merger with a lower mass star, with the primary remaining with low N abundance until core collapse. Given the available merger models, we can only put a loose constraint on the pre-explosion mass around 10-20 Msun, with lower values being favored based on previous observational constraints from the nebular phase., 10 pages, 8 figures, MNRAS accepted

Photometric observations of the nearly equal-mass binary near-Earth asteroid 2017 YE5 were carried out at the Observatório Astronômico do Sertão de Itaparica (OASI, Brazil) and at the Blue Mountains Observatory (BMO, Australia) between 2018 July and August, shortly after it made a close approach to Earth in 2018 June. These observations allowed to determine an orbital period of the system in good agreement with the result of the radar observations. Our results also indicate that the 2017 YE5 system is not fully locked in a synchronous spin–orbit resonance, as there is a possible asynchronous component or tumbling rotator in the system. Additional data obtained at the Observatorio Astronómico Nacional de San Pedro Mártir (OAN-SPM, Mexico) in 2018 August allowed to derive the colour indices and the low-resolution spectrum, which indicate that the object has a very reddish surface similar to the outer Solar system objetcs. This is consistent with a dark albedo of about 3 per cent, derived from an effective diameter of the combined components and an absolute magnitude of 19.2. Infrared data obtained at the NASA’s Infrared Telescope Facility exhibit a thermal emission in the 2.5 $\mu$m range for which a low albedo of 2–4 per cent was adjusted by applying a thermal model. Regarding to the taxonomic classification, we found that the 2017 YE5 is a D-type in the Bus-DeMeo taxonomy. Therefore, according to our results and considering that the 2017 YE5 system has a typical comet orbit (TJ = 2.87), we suggest that it is a dormant Jupiter-family binary comet.

One possibility for detecting low-amplitude pulsational variations is through gravitational microlensing. During a microlensing event, the temporary brightness increase leads to improvement in the signal-to-noise ratio, and thereby better detectability of pulsational signatures in light curves. We explore this possibility under two primary considerations. The first is when the standard point-source and point-lens approximation applies. In this scenario, dividing the observed light curve by the best-fitted microlensing model leads to residuals that result in pulsational features with improved uncertainties. The second is for transit events (single lens) or caustic crossing (binary lens). The point-source approximation breaks down, and residuals relative to a simple best-fitted microlensing model display more complex behavior. We employ a Monte-Carlo simulation of microlensing of pulsating variables toward the Galactic bulge for the surveys of OGLE and of KMTNet. We demonstrate that the efficiency for detecting pulsational signatures with intrinsic amplitudes of $350$, is $\sim 50-60\%$. The maximum efficiency occurs for pulsational periods $P \simeq 0.1-0.3$ days. We also study the possibility that high-magnification microlensing events of non-radially pulsating stars (NRPs) could be misinterpreted as planetary or binary microlensing events. We conclude that small asymmetric features around lightcurve peaks due to stellar pulsations could be misdiagnosed with crossing (or passing close to) small caustic curves., Comment: 11 pages, 4 Figures, 4 Tables

Analysis of timing variation of the eclipsing binary RR Dra implies the existence of an additional object with mass of no less than 3.5(±0.2) M⊙. This third object is orbiting around the central eclipsing pair once every 98(±1) yr, making the whole a hierarchical triple system. However, reliable photometric solutions based on light curves from Transiting Exoplanet Survey Satellite and WASP show that contribution of the third-light takes only about 2 per cent of the total luminosity. It could photometrically be attributed to other unrelated star located within the aperture. The tertiary component is thus a massive object but invisible in optical wavelengths. Besides, evidence of misalignment between the barycenter and the centre of light of the system is also found. This strengthen the existence of a hidden black hole candidate in the form of the third body. The potential black hole may play an essential role in extracting angular momentum from the central binary pair, forming the current state. As a good laboratory, RR Dra is an interesting system that can help to study stellar-mass black hole under the circumbinary case.

Young star clusters are dynamically active stellar systems and are a common birthplace for massive stars. Low-mass star clusters ($\sim{}300-10^3$ M$_\odot$) are more numerous than massive systems and are characterized by a two-body relaxation time scale of a few Myr: the most massive stars sink to the cluster core and dynamically interact with each other even before they give birth to compact objects. Here, we explore the properties of black holes (BHs) and binary black holes (BBHs) formed in low-mass young star clusters, by means of a suite of $10^5$ direct $N$-body simulations with a high original binary fraction (100 % for stars with mass $>5$ M$_\odot$). Most BHs are ejected in the first $\sim{}20$ Myr by dynamical interactions. Dynamical exchanges are the main formation channel of BBHs, accounting for $\sim{}40-80$ % of all the systems. Most BBH mergers in low-mass young star clusters involve primary BHs with mass $, 14 pages, 13 figures, comments welcome

Common envelope evolution (CEE) is believed to be an important stage in the evolution of binary/multiple stellar systems. Following this stage, the CE is thought to be ejected, leaving behind a compact binary (or a merger product). Although extensively studied, the CEE process is still little understood, and although most binaries have non-negligible eccentricity, the effect of initial eccentricity on the CEE has been little explored. Moreover, most studies assume a complete circularization of the orbit by the CE onset, while observationally such eccentricities are detected in many post-CE binaries. Here we use smoothed particle hydro-dynamical simulations (SPH) to study the evolution of initially eccentric ($0\le e\le0.95$) CE-systems. We find that initially eccentric binaries only partially circularize. In addition, higher initial eccentricity leads to a higher eccentricity following the end of the inspiral phase, with eccentricities as high as 0.18 in the most eccentric cases, and even higher if the initial peri-center of the orbit is located inside the star (e.g. following a kick into an eccentric orbit, rather than a smooth transition). CEE of more eccentric binaries leads to enhanced dynamical mass-loss of the CE compared with more circular binaries, and depends on the initial closest approach of the binary. We show that our results and the observed eccentricities of post-CE binaries suggest that the typical assumptions of circular orbits following CEE might potentially be revised. We expect post-CE eccentricities to affect the delay time distributions of various transients such as supernovae, gamma-ray bursts and gravitational-wave sources by up to tens of percents., Accepted to MNRAS

The vast majority of binaries containing a compact object and a regular star spend most of their time in a quiescent state where no strong interactions occur between components. Detection of these binaries is extremely challenging and only few candidates have been detected through optical spectroscopy. Self-lensing represents a new means of detecting compact objects in binaries, where gravitational lensing of the light from the visible component by the compact object produces periodic optical flares. Here we show that current and planned large-area optical surveys can detect a significant number ($\sim 100$-$10,000$s) of these self-lensing binaries and provide insights into the properties of the compact lenses. We show that many of the predicted population of observable self-lensing binaries will be observed with multiple self-lensing flares; this both improves the chances of detection and also immediately distinguishes them from chance-alignment micro-lensing events. Through self-lensing we can investigate long - but previously hidden - stages of binary evolution and consequently provide new constraints on evolutionary models which impact on the number and nature of double compact object mergers., 11 pages, 6 figures, 3 tables

Over the past 3 yr, the fading non-thermal emission from the GW170817 remained generally consistent with the afterglow powered by synchrotron radiation produced by the interaction of the structured jet with the ambient medium. Recent observations by Hajela et al. indicate the change in temporal and spectral behaviour in the X-ray band. We show that the new observations are compatible with the emergence of a new component due to non-thermal emission from the fast tail of the dynamical ejecta of ab-initio binary neutron star merger simulations. This provides a new avenue to constrain binary parameters. Specifically, we find that equal mass models with soft equations of state (EOSs) and high-mass ratio models with stiff EOSs are disfavoured as they typically predict afterglows that peak too early to explain the recent observations. Moderate stiffness and mass ratio models, instead, tend to be in good overall agreement with the data.

New multicolour photometric observations were carried out on 22 nights in three observation missions between 2015 October and 2016 February for eclipsing binary HIP 7666. High- and low-resolution spectroscopic observations were also carried out in the winters of 2015, 2016, and 2019, respectively. The fully phase-covered light curves and radial velocity curves are presented. All times of light minima are used to calculate the orbital period [2.372 2200(4) d] and new ephemerides. The photometric solution and stellar physical parameters are derived, showing that HIP 7666 is a detached binary with the absolute parameters M1 = 1.53(3) M⊙, R1 = 2.08(2) R⊙, log L1/L⊙ = 0.99(3), log g1 = 3.98(2), Mbol1 = 2.26(8), Teff,1 = 7100(100) K for the primary, and M2 = 1.23(3) M⊙, R2 = 1.06(2) R⊙, log L2/L⊙ = 0.13(4), log g2 = 4.47(2), Mbol2 = 4.4(1), Teff,2 = 6029(67) K for the secondary. The age of 1.69 Gyr and 1.755 Gyr are estimated from parsec isochrones and mesa evolutionary tracks, respectively. In Herzsprung–Russell (H–R) diagram, the primary component has evolved to the late stages, while the secondary still locates at the early stages of the main sequence. Frequency analysis yields three frequencies of f1 = 24.631(4) cd−1, f2 = 21.193(1) cd−1, and f3 = 28.07(7) cd−1. Comparisons between models and observations suggest that the primary component is most likely a p-mode non-radial δ Scuti oscillator.

We revisit the well-studied outburst of the low mass X-ray binary (LMXB) system V404 Cygni, and claims of gamma-ray excesses observed with the Fermi-LAT instrument. Upon analysing an 11.5 year dataset with the 8-year LAT point source catalogue and 8-year background models, we find no evidence to suggest that there is high energy gamma-ray emission during the outburst period (or at any other time) from V404 Cygni. This is due to the proximity of V404 Cygni to the gamma-ray emitting blazar B2023+336, a luminous source approximately 0.3 degrees away, which causes source confusion at the position of V404 Cygni, the luminous gamma-ray background, and the use of older background models and catalogues in previous studies., Comment: 11 pages, 6 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

We study the enrichment of the interstellar medium with rapid neutron capture (r-process) elements produced in binary neutron star (BNS) mergers. We use a semi-analytic model to describe galactic evolution, with merger rates and time delay distributions of BNS mergers consistent with the latest population synthesis models. In order to study the dispersion of the relative abundances of r-process elements and iron, we applied a turbulent mixing scheme, where the freshly synthesized elements are gradually dispersed in the interstellar medium. We show that within our model the abundances observed in Milky-Way stars, in particular the scatter at low metallicities, can be entirely explained by BNS mergers. Our results suggest that binary neutron star mergers could be the dominant source of r-process elements in the Galaxy., Comment: 16 pages, 19 figures. Modified to match the version accepted by MNRAS (some changes to main results and conclusions, new section added)

The parameter space for modelling stellar systems is vast and complicated. To find best-fitting models for a star one needs a statistically robust way of exploring this space. We present a new machine-learning approach to predict the modelling parameters for detached double-lined eclipsing binary systems, including the system age, based on observable quantities. Our method allows for the estimation of the importance of several physical effects which are included in a parametrized form in stellar models, such as convective core overshoot or stellar spot coverage. The method yields probability distribution functions for the predicted parameters which take into account the statistical and, to a certain extent, the systematic errors which is very difficult to do using other methods. We employ two different approaches to investigate the two components of the system either independently or in a combined manner. Furthermore, two different grids are used as training data. We apply the method to 26 selected objects and test the predicted best solutions with an on-the-fly optimization routine which generates full hydrostatic models. While we do encounter failures of the predictions, our method can serve as a rapid estimate for stellar ages of detached eclipsing binaries taking full account of the uncertainties in the observables.

Θ Mus = HD 113904 is a massive multiple system containing the WC5/6 + O6/7V binary WR48 in a $19.1\, \mathrm{ d}$ circular orbit. Previous attempts to constrain the variable photometric properties of this binary subsystem have been thwarted by the dominating stochastically variable light from a 10-times brighter blue supergiant (BSG), located only 46 mas away. Even now, with extensive optical space-based photometry from one of the BRITE-Constellation satellites, we were unable to beat down the intrinsic stochastic variability from the BSG enough to provide a convincing detection of a low-level atmospheric eclipse of the WC + O system, as often seen in other short-period WR + O systems. We explore the variability of the dominating BSG and find that its behaviour is similar to that of other BSGs, with a forest of low-frequency Fourier peaks likely from stochastic gravity waves reaching the stellar surface. Then, by adopting an orbital inclination from another more reliable source, we obtain a clumping-independent, linear-density-dependent upper limit of the mass-loss rate for the WR component of $(6.5 \pm 0.3) \times 10^{-5}\, \mathrm{ M}_{\odot }\,\mathrm{ yr}^{ -1}$, which is consistent with values of other WC5/6 subtypes. This corresponds to an upper limit of 5.0 ± 0.2 mmag for the depth of the atmospheric eclipse in the WR48 subsystem when observed together with the BSG.

Extremely low-mass white dwarfs (ELM WDs) are the result of binary evolution in which a low-mass donor star is stripped by its companion leaving behind a helium-core white dwarf. We explore the formation of ELM WDs in binary systems considering the Convection And Rotation Boosted magnetic braking treatment. Our evolutionary sequences were calculated using the MESA code, with initial masses of 1.0 and 1.2 Msun (donor), and 1.4 (accretor), compatible with low mass X-ray binaries (LMXB) systems. We obtain ELM models in the range 0.15 to 0.27 Msun from a broad range of initial orbital periods, 1 to 25 d. The bifurcation period, where the initial period is equal to the final period, ranges from 20 to 25 days. In addition to LMXBs, we show that ultra-compact X-ray binaries (UCXB) and wide-orbit binary millisecond pulsars can also be formed. The relation between mass and orbital period obtained is compatible with the observational data from He white dwarf companions to pulsars., Comment: 16 pages, 6 figures. Accepted for publication in MNRAS

The Legacy Survey of Space and Time (LSST) by the Vera C. Rubin Observatory is expected to discover tens of millions of quasars. A significant fraction of these could be powered by coalescing massive black hole (MBH) binaries, since many quasars are believed to be triggered by mergers. We show that under plausible assumptions about the luminosity functions, lifetimes, and binary fractions of quasars, we expect the full LSST quasar catalogue to contain between 20 and 100 million compact MBH binaries with masses M = 105–9M⊙, redshifts z = 0–6, and orbital periods P = 1–70 d. Their light-curves are expected to be distinctly periodic, which can be confidently distinguished from stochastic red-noise variability, because LSST will cover dozens, or even hundreds of cycles. A very small subset of 10–150 ultracompact (P ≲ 1 d) binary quasars among these will, over ∼5–15 yr, evolve into the mHz gravitational-wave frequency band and can be detected by LISA. They can therefore be regarded as ‘LISA verification binaries’, analogous to short-period Galactic compact-object binaries. The practical question is how to find these handful of ‘needles in the haystack’ among the large number of quasars: this will likely require a tailored co-adding analysis optimized for this purpose.

Peters' formula is an analytical estimate of the time-scale of gravitational wave (GW)-induced coalescence of binary systems. It is used in countless applications, where the convenience of a simple formula outweighs the need for precision. However, many promising sources of the Laser Interferometer Space Antenna (LISA), such as supermassive black hole binaries and extreme mass-ratio inspirals (EMRIs), are expected to enter the LISA band with highly eccentric ($e \gtrsim$ 0.9) and highly relativistic orbits. These are exactly the two limits in which Peters' estimate performs the worst. In this work, we expand upon previous results and give simple analytical fits to quantify how the inspiral time-scale is affected by the relative 1.5 post-Newtonian (PN) hereditary fluxes and spin-orbit couplings. We discuss several cases that demand a more accurate GW time-scale. We show how this can have a major influence on quantities that are relevant for LISA event-rate estimates, such as the EMRI critical semi-major axis. We further discuss two types of environmental perturbations that can play a role in the inspiral phase: the gravitational interaction with a third massive body and the energy loss due to dynamical friction and torques from a surrounding gas medium ubiquitous in galactic nuclei. With the aid of PN corrections to the time-scale in vacuum, we find simple analytical expressions for the regions of phase space in which environmental perturbations are of comparable strength to the effects of any particular PN order, being able to qualitatively reproduce the results of much more sophisticated analyses., Accepted for publication in MNRAS. Comments welcome!

The detection of gravitational waves from a neutron star merger, GW170817, marked the dawn of a new era in time-domain astronomy. Monitoring of the radio emission produced by the merger, including high-resolution radio imaging, enabled measurements of merger properties including the energetics and inclination angle. In this work, we compare the capabilities of current and future gravitational wave facilities to the sensitivity of radio facilities to quantify the prospects for detecting the radio afterglows of gravitational wave events. We consider three observing strategies to identify future mergers – wide field follow-up, targeting galaxies within the merger localization and deep monitoring of known counterparts. We find that while planned radio facilities like the Square Kilometre Array will be capable of detecting mergers at gigaparsec distances, no facilities are sufficiently sensitive to detect mergers at the range of proposed third-generation gravitational wave detectors that would operate starting in the 2030s.

We carry out a systematic study of the response of companion stars in massive binaries after being impacted by supernova ejecta. A total of 720 1D stellar evolution calculations are performed to follow the inflation and contraction of the star in response to the energy injection and how it depends on various parameters. We find that the maximum luminosity achieved during the inflated phase is only dependent on the stellar mass and we derive an analytic formula to describe the relation. There is also a tight correlation between the duration of expansion and the intersected energy. These correlations will be useful to constrain pre-supernova binary parameters from future detections of inflated companions. We also discuss the possible outcomes of the binary system when the companion inflation is taken into account. Based on simple binary population synthesis, we estimate that $\sim$1-3% of stripped-envelope supernovae may have observable inflated companions. Finally, we apply our models to the observed companion of SN2006jc and place strong constraints on the possible pre-supernova binary parameters., 15 pages, 15 figures, Accepted for publication in MNRAS

We show that the odds of the mass-gap (secondary) object in GW190814 being a neutron star (NS) improve if one allows for a stiff high-density equation of state (EoS) or a large spin. Since its mass is $\in (2.50,2.67) M_{\odot}$, establishing its true nature will make it either the heaviest NS or the lightest black hole (BH), and can have far-reaching implications on NS EoS and compact object formation channels. When limiting oneself to the NS hypothesis, we deduce the secondary's properties by using a Bayesian framework with a hybrid EoS formulation that employs a parabolic expansion-based nuclear empirical parameterization around the nuclear saturation density augmented by a generic 3-segment piecewise polytrope (PP) model at higher densities and combining a variety of astrophysical observations. For the slow-rotation scenario, GW190814 implies a very stiff EoS and a stringent constraint on the EoS specially in the high-density region. On the other hand assuming the secondary object is a rapidly rotating NS, we constrain its rotational frequency to be $f=1170^{+389}_{-495}$ Hz, within a $90\%$ confidence interval. In this scenario, the secondary object in GW190814 would qualify as the fastest rotating NS ever observed. However, for this scenario to be viable, rotational instabilities would have to be suppressed both during formation and the subsequent evolution until merger, otherwise the secondary of GW190814 is more likely to be a BH., Comment: Matches with the MNRAS published version

We explore chemically homogeneous evolution (CHE) as a formation channel for massive merging binary black holes (BBHs). We develop methods to include CHE in a rapid binary population synthesis code, Compact Object Mergers: Population Astrophysics and Statistics (COMPAS), which combines realistic models of binary evolution with cosmological models of the star-formation history of the Universe. For the first time, we simultaneously explore conventional isolated binary star evolution under the same set of assumptions. This approach allows us to constrain population properties and make simultaneous predictions about the gravitational-wave detection rates of BBH mergers for the CHE and conventional formation channels. The overall mass distribution of detectable BBHs is consistent with existing gravitational-wave observations. We find that the CHE channel may yield up to ~70% of all gravitational-wave detections of BBH mergers coming from isolated binary evolution., 15 pages, 16 figures

Trans-Neptunian binaries are among the objects carrying the most relevant information on the conditions that prevailed in the primitive outer Solar system. However, their present orbital architecture not necessarily reflects those conditions, because it is the result of complex evolutionary processes. There is evidence suggesting that binary systems dominated the primitive trans-Neptunian populations. Therefore, binary–binary interactions could be an important evolutionary path. In this paper, a series of numerical experiments of binary–binary close encounters were carried out, varying specific parameters related to their physical and orbital characteristics. The relative speeds with which the encounters occur play an important role in the relative variation of orbital elements, showing a good agreement with impulsive approximation models involving a binary and a single intruder. However, simulations varying the mass, semimajor axis, and closest approach distance between pairs of binaries, showed marked deviations from this model.

We study the detached eclipsing binary, KIC 5359678, with starspot modulation using the high-quality Kepler photometry and LAMOST spectroscopy. The phoebe model, optimal for this binary, reveals that this system is a circular detached binary, composed of two F-type main-sequence stars. The masses and radii of the primary and the secondary are M1 = 1.31 ± 0.05 M⊙, R1 = 1.52 ± 0.04 R⊙, M2 = 1.12 ± 0.04 M⊙, and R2 = 1.05 ± 0.06 R⊙, respectively. The age of this binary is estimated to be about 2 Gyr, a value much longer than the synchronization time-scale of 17.8 Myr. The residuals of light curves show quasi-sinusoidal signals, which could be induced by starspots. We apply autocorrelation function analysis on the out-of-eclipse residuals and find that the spot with rotational period close to the orbital period, while, the decay time-scale of starspots is longer than that on the single stars with the same temperature, period range, and rms scatter. A two-starspot model is adopted to fit the signals with two-dip pattern, whose result shows that the longitude decreases with time.

We present the results of a medium resolution optical spectroscopic survey of 92 cool ($3,000 \lesssim T_{\rm eff} \lesssim 4,500\,$K) southern TESS candidate planet hosts, and describe our spectral fitting methodology used to recover stellar parameters. We quantify model deficiencies at predicting optical fluxes, and while our technique works well for $T_{\rm eff}$, further improvements are needed for [Fe/H]. To this end, we developed an updated photometric [Fe/H] calibration for isolated main sequence stars built upon a calibration sample of 69 cool dwarfs in binary systems, precise to $\pm0.19\,$dex, from super-solar to metal poor, over $1.51 < {\rm Gaia}~(B_P-R_P) < 3.3$. Our fitted $T_{\rm eff}$ and $R_\star$ have median precisions of 0.8% and 1.7%, respectively and are consistent with our sample of standard stars. We use these to model the transit light curves and determine exoplanet radii for 100 candidate planets to 3.5% precision and see evidence that the planet-radius gap is also present for cool dwarfs. Our results are consistent with the sample of confirmed TESS planets, with this survey representing one of the largest uniform analyses of cool TESS candidate planet hosts to date., Accepted version. 33 Pages, 12 Figures, 9 Tables

Novae are some of the most commonly detected optical transients and have the potential to provide valuable information about binary evolution. Binary population synthesis codes have emerged as the most effective tool for modelling populations of binary systems, but such codes have traditionally employed greatly simplified nova physics, precluding detailed study. In this work, we implement a model treating H and He novae as individual events into the binary population synthesis code \binaryc. This treatment of novae represents a significant improvement on the `averaging' treatment currently employed in modern population synthesis codes. We discuss the evolutionary pathways leading to these phenomena and present nova event rates and distributions of several important physical parameters. Most novae are produced on massive white dwarfs, with approximately 70 and 55 per cent of nova events occurring on O/Ne white dwarfs for H and He novae respectively. Only 15 per cent of H-nova systems undergo a common-envelope phase, but these systems are responsible for the majority of H nova events. All He-accreting He-nova systems are considered post-common-envelope systems, and almost all will merge with their donor star in a gravitational-wave driven inspiral. We estimate the current annual rate of novae in M31 (Andromeda) to be approximately $41 \pm 4$ for H novae, underpredicting the current observational estimate of $65^{+15}_{-16}$, and $0.14\pm0.015$ for He novae. When varying common-envelope parameters, the H nova rate varies between 20 and 80 events per year., Accepted, MNRAS. 7 Jun 2020: Minor correction regarding AM CVn masses at period bounce, courtesy of P. Neuteufel

Many stars are in binaries or higher-order multiple stellar systems. Although in recent years a large number of binaries have been proven to host exoplanets, how planet formation proceeds in multiple stellar systems has not been studied much yet from the theoretical standpoint. In this paper we focus on the evolution of the dust grains in planet-forming discs in binaries. We take into account the dynamics of gas and dust in discs around each component of a binary system under the hypothesis that the evolution of the circumprimary and the circumsecondary discs is independent. It is known from previous studies that the secular evolution of the gas in binary discs is hastened due to the tidal interactions with their hosting stars. Here we prove that binarity affects dust dynamics too, possibly in a more dramatic way than the gas. In particular, the presence of a stellar companion significantly reduces the amount of solids retained in binary discs because of a faster, more efficient radial drift, ultimately shortening their lifetime. We prove that how rapidly discs disperse depends both on the binary separation, with discs in wider binaries living longer, and on the disc viscosity. Although the less-viscous discs lose high amounts of solids in the earliest stages of their evolution, they are dissipated slowly, while those with higher viscosities show an opposite behaviour. The faster radial migration of dust in binary discs has a striking impact on planet formation, which seems to be inhibited in this hostile environment, unless other disc substructures halt radial drift further in. We conclude that if planetesimal formation were viable in binary discs, this process would take place on very short time scales., 18 pages, 11 figures; accepted for publication in MNRAS

Swift J004929.5-733107 is an X-ray source in the Small Magellanic Cloud (SMC) that has been reported several times, but the optical counterpart has been unclear due to source confusion in a crowded region of the SMC. Previous works proposed [MA93] 302 as the counterpart, however we show here, using data obtained from the S-CUBED project, that the X-ray positio is inconsistent with that object. Instead we propose a previously unclassified object which has all the indications of being a newly identified Be star exhibiting strong HU emission. Evidence for the presence of significant I-band variability strongly suggest that this is, in fact, a Be type star with a large circumstellar disk. Over 18 years worth of optical monitoring by the OGLE project reveal a periodic modulation at a period of 413d, probably the binary period of the system. A SALT optical spectrum shows strong Balmer emission and supports a proposed spectral classification of B1-3 III-IVe. The X-ray data obtained from the S-CUBED project reveal a time-averaged spectrum well fitted by a photon index = 0.93 pm 0.16. Assuming the known distance to the SMC the flux corresponds to a luminosity 10E35 erg/s. All of these observational facts suggest that this is confirmed as a Be star-neutron star X-ray binary (BeXRB) in the SMC, albeit one with an unusually long binary period at the limits of the Corbet Diagram., 9 pages 12 figures

We report on the Insight-HXMT observations of the new black hole X-ray binary MAXI J1820+070 during its 2018 outburst. Detailed spectral analysis via the continuum fitting method shows an evolution of the inferred spin during its high soft sate. Moreover, the hardness ratio, the non-thermal luminosity and the reflection fraction also undergo an evolution, exactly coincident to the period when the inferred spin transition takes place. The unphysical evolution of the spin is attributed to the evolution of the inner disc, which is caused by the collapse of a hot corona due to condensation mechanism or may be related to the deceleration of a jet-like corona. The studies of the inner disc radius and the relation between the disc luminosity and the inner disc radius suggest that, only at a particular epoch, did the inner edge of the disc reach the innermost stable circular orbit and the spin measurement is reliable. We then constrain the spin of MAXI J1820+070 to be a*=0.2^{+0.2}_{-0.3}. Such a slowly spinning black hole possessing a strong jet suggests that its jet activity is driven mainly by the accretion disc rather than by the black hole spin., Comment: 14 pages, 13 figures, 5 tables, accepted for publication in MNRAS

High-resolution spectra ($R\sim 48\, 000$) of the massive binary system HH Carinae have been analysed. Precise absolute parameters were derived from a simultaneous solution of the radial velocities and the light curves. The primary component is found to be an O9-type main-sequence star with a temperature of 33 500 K, while the secondary component is a B0-type giant/subgiant star with a temperature of 27 500 K. An analysis of the spectroscopic and photometric data has shown that the primary component rotates at a speed of vrot1 = 220 km s−1, which is three times faster than the synchronous rotation, while the secondary component synchronously rotates with the orbit at a speed of vrot2 = 150 km s−1. The distance to the system and the velocity of the centre of mass are determined as d = 4.6 ± 0.8 kpc and Vγ = –16 km s−1, respectively. The distance of the system is in agreement with the most probable Gaia distance of $4.9^{+0.9}_{-0.7}$ kpc and the distance in the latest data release (DR3) of 4.4 ± 0.3. Emission structures at the Hα region were modelled using the code shellspec where the derived absolute parameters of the components have been considered. Because the components are massive stars, mass loss as a result of stellar winds is expected. Produced models confirm that the components do indeed have strong stellar winds and there is mass transfer from the secondary to the primary. Stellar winds and the gas stream between the components have been modelled as a hot shell around the system, with a temperature of ∼22 000 K. Models also indicate that the interaction between the wind and the gas stream causes the formation of a high-temperature (100 000 K) impact region.

We present an analysis of the most massive white dwarf candidates in the Montreal White Dwarf Database 100 pc sample. We identify 25 objects that would be more massive than $1.3~M_{\odot}$ if they had pure H atmospheres and CO cores, including two outliers with unusually high photometric mass estimates near the Chandrasekhar limit. We provide follow-up spectroscopy of these two white dwarfs and show that they are indeed significantly below this limit. We expand our model calculations for CO core white dwarfs up to $M=1.334\ M_\odot$, which corresponds to the high-density limit of our equation-of-state tables, $\rho = 10^9$ g cm$^{-3}$. We find many objects close to this maximum mass of our CO core models. A significant fraction of ultramassive white dwarfs are predicted to form through binary mergers. Merger populations can reveal themselves through their kinematics, magnetism, or rapid rotation rates. We identify four outliers in transverse velocity, four likely magnetic white dwarfs (one of which is also an outlier in transverse velocity), and one with rapid rotation, indicating that at least 8 of the 25 ultramassive white dwarfs in our sample are likely merger products., Comment: MNRAS, in press

BG Ind is a well studied, bright, nearby binary consisting of a pair of F stars in a 1.46-day orbit. We have discovered in the TESS lightcurve for TIC 229804573 (aka BG Ind) a second eclipsing binary in the system with a 0.53-day. Our subsequent analyses of the recent TESS and archival ground-based photometric and radial velocity data, reveal that the two binaries are gravitationally bound in a 721-day period, moderately eccentric orbit. We present the results of a joint spectro-photodynamical analysis of the eclipse timing variation curves of both binaries based on TESS and ground-based archival data, the TESS lightcurve, archival radial velocity data and the spectral energy distribution, coupled with the use of PARSEC stellar isochrones. We confirm prior studies of BG Ind which found that the brighter binary A consists of slightly evolved F-type stars with refined masses of 1.32 and 1.43 $M_\odot$, and radii of 1.59 and 2.34 $R_\odot$. The previously unknown binary B has two less massive stars of 0.69 and 0.64 $M_\odot$ and radii of 0.64 and 0.61 $R_\odot$. Based on a number of different arguments which we discuss, we conclude that the three orbital planes are likely aligned to within 17$^\circ$., 17 pages. Accepted for publication in MNRAS

New spectrometric data on V Pup are combined with satellite photometry (HIPPARCOS and recent TESS) to allow a revision of the absolute parameters with increased precision. We find: $M_1$ = 14.0$\pm$0.5, $M_2$ = 7.3$\pm$0.3 (M$_\odot$); $R_{1}$ = 5.48$\pm$0.18, $R_2$ = 4.59$\pm$0.15 (R$_\odot$); $T_{1}$ 26000$\pm 1000$, $T_2$ 24000 $\pm$1000 (K), age 5 $\pm$1 (Myr), photometric distance 320 $\pm$10 (pc). The TESS photometry reveals low-amplitude ($\sim$0.002 mag) variations of the $\beta$ Cep kind, consistent with the deduced evolutionary condition and age of the optical primary. This fact provides independent support to our understanding of the system as in a process of Case A type interactive evolution that can be compared with $\mu^1$ Sco. The $\sim$10 M$_{\odot}$ amount of matter shed by the over-luminous present secondary must have been mostly ejected from the system rather than transferred, thus taking angular momentum out of the orbit and keeping the pair in relative close proximity. New times of minima for V Pup have been studied and the results compared with previous analyses. The implied variation of period is consistent with the Case A evolutionary model, though we offer only a tentative sketch of the original arrangement of this massive system. We are not able to confirm the previously reported cyclical variations having a 5.47 yr period with the new data, though a direct comparison between the HIPPARCOS and TESS photometry points to the presence of third light from a star that is cooler than those of the close binary, as mentioned in previous literature., Comment: 13 pages, 12 figures, 8 tables. Accepted by MNRAS

We introduce here our new approach to modeling particle cloud evolution off surface of small bodies (asteroids and comets), following the evolution of ejected particles requires dealing with various time and spatial scales, in an efficient, accurate and modular way. In order to improve computational efficiency and accuracy of such calculations, we created an N-body modeling package as an extension to the increasingly popular orbital dynamics N-body integrator Rebound. Our code is currently a stand-alone variant of the Rebound code and is aimed at advancing a comprehensive understanding of individual particle trajectories, external forcing, and interactions, at the scale which is otherwise overlooked by other modeling approaches. The package we developed -- Rebound Ejecta Dynamics (RED) -- is a Python-based implementation with no additional dependencies. It incorporates several major mechanisms that affect the evolution of particles in low-gravity environments and enables a more straightforward simulation of combined effects. We include variable size and velocity distributions, solar radiation pressure, ellipsoidal gravitational potential, binary or triple asteroid systems, and particle-particle interactions. In this paper, we present a sample of the RED package capabilities. These are applied to a small asteroid binary system (characterized following the Didymos/Dimorphos system, which is the target for NASA's Double Asteroid Redirection Test mission), Comment: 13 pages, 11 figures, 4 tables, 2 appendices. Published (advance access) on 2/15/2021, see https://academic.oup.com/mnras/advance-article-abstract/doi/10.1093/mnras/stab408/6136266?redirectedFrom=fulltext

We introduce a phenomenological, physically motivated, model for the effective tidal deformability of a neutron star, adding the frequency dependence (associated with the star's fundamental mode of oscillation) that comes into play during the late stages of the binary inspiral. Testing the model against alternative descriptions, we demonstrate that it provides an accurate representation of the dynamical tide up to close to merger. The simplicity of the prescription makes it an attractive alternative for a gravitational-wave data analysis implementation, facilitating an inexpensive construction of a large number of templates covering the relevant parameter space., Revised version with added details. To appear in MNRAS

The metallicity dependence of the wide-binary fraction in stellar populations plays a critical role in resolving the open question of wide binary formation. In this paper, we investigate the metallicity ([Fe/H]) and age dependence of the wide-binary fraction (binary separations between $10^3$ and $10^4$ AU) for field F and G dwarfs within 500 pc by combining their metallicity and radial velocity measurements from LAMOST DR5 with the astrometric information from Gaia DR2. We show that the wide-binary fraction strongly depends on the metallicity: as metallicity increases, the wide-binary fraction first increases, peaks at [Fe/H]$\simeq 0$, and then decreases at the high metallicity end. The wide-binary fraction at [Fe/H]$=0$ is about two times larger than that at [Fe/H]$=-1$ and [Fe/H]$=+0.5$. This metallicity dependence is dominated by the thin-disk stars. Using stellar kinematics as a proxy of stellar age, we show that younger stars have a higher wide-binary fraction at fixed metallicity close to solar. We propose that multiple formation channels are responsible for the metallicity and age dependence. In particular, the positive metallicity correlation at [Fe/H]$0$ can be inherited from the similar metallicity dependence of close binaries, and radial migration may play a role in enhancing the wide-binary fraction around the solar metallicity., Accepted by MNRAS. Fig 3 and 7 are the key results

The recently discovered colliding-wind binary (CWB) Apep has been shown to emit luminously from radio to X-rays, with the emission driven by a binary composed of two Wolf-Rayet (WR) stars of one carbon-sequence (WC8) and one nitrogen-sequence (WN4-6b). Mid-infrared imaging revealed a giant spiral dust plume that is reminiscent of a pinwheel nebula but with additional features that suggest Apep is a unique system. We have conducted observations with the Australian Long Baseline Array to resolve Apep's radio emission on milliarcsecond scales, allowing us to relate the geometry of the wind-collision region to that of the spiral plume. The observed radio emission shows a bow-shaped structure, confirming its origin as a wind-collision region. The shape and orientation of this region is consistent with being originated by the two stars and with being likely dominated by the stronger wind of the WN4-6b star. This shape allowed us to provide a rough estimation of the opening angle of $\sim 150^\circ$ assuming ideal conditions. The orientation and opening angle of the emission also confirms it as the basis for the spiral dust plume. We also provide estimations for the two stars in the system to milliarcsecond precision. The observed radio emission, one order of magnitude brighter and more luminous than any other known non-thermal radio-emitting CWB, confirms it is produced by an extremely powerful wind collision. Such a powerful wind-collision region is consistent with Apep being a binary composed of two WR stars, so far the first unambiguously confirmed system of its kind., Comment: 9 pages, 3 figures, accepted for publication in MNRAS

We present in this paper a first step toward a semi-analytical treatment of the accretion process in wind-fed neutron star supergiant X-ray binaries with eccentric orbits. We consider the case of a spherically symmetric wind for the supergiant star and a simplified model for the accretion onto the compact object. A self-consistent calculation of the photoionization of the stellar wind by the X-rays from the accreting neutron star is included. This effect is convolved with the modulation of the mass accretion rate induced by the eccentric orbit to obtain the expected X-ray luminosity of a system along the orbit. As part of our results, we first show that the bi-modality of low and high X-ray luminosity solutions for supergiant X-ray binaries reported in previous papers is likely to result from the effect of the neutron star approaching first and then moving away from the companion (without coexisting simultaneously). We propose that episodes of strong wind photoionization can give rise to off-states of the sources. Our calculations are applied to the case of a few classical supergiant X-ray binary systems with known eccentricities (Vela X-1, 4U 1907+09, GX 301-2) and to the case of the only supergiant fast X-ray transient with a confirmed eccentric orbit, IGR J08408-4503. The results are compared with observational findings on these sources. We also discuss the next steps needed to expand the calculations toward a more comprehensive treatment in future publications., Comment: Accepted for publication in MNRAS

The intermediate mass Herbig Ae star V1787 Ori is a member of the L1641 star-forming region in the Orion A molecular cloud. We report the detection of an M-type companion to V1787 Ori at a projected separation of 6.66" (corresponding to 2577 au), from the analysis of VLT/NACO adaptive optics $K_s$-band image. Using astrometric data from Gaia DR2, we show that V1787 Ori A and B share similar distance ($d$ $\sim$ 387 pc) and proper motion, indicating that they are physically associated. We estimate the spectral type of V1787 Ori B to be M5 $\pm$ 2 from color--spectral type calibration tables and template matching using SpeX spectral library. By fitting PARSEC models in the Pan-STARRS color-magnitude diagram, we find that V1787 Ori B has an age of 8.1$^{+1.7}_{-1.5}$ Myr and a mass of 0.39$^{+0.02}_{-0.05}$ $M_\odot$. We show that V1787 Ori is a pre-main sequence wide binary system with a mass ratio of 0.23. Such a low mass ratio system is rarely identified in Herbig Ae/Be binary systems. We conclude this work with a discussion on possible mechanisms for the formation of V1787 Ori wide binary system., 10 pages,3 tables, 5 figures, accepted for publication in MNRAS

While the majority of massive stars have a stellar companion, most pulsars appear to be isolated. Taken at face value, this suggests that most massive binaries break apart due to strong natal kicks received in supernova explosions. However, the observed binary fraction can still be subject to strong selection effects, as monitoring of newly discovered pulsars is rarely carried out for long enough to conclusively rule out multiplicity. Here, we use the second Gaia Data Release (DR2) to search for companions to 1534 rotation-powered pulsars with positions known to better than 0.5 arcseconds. We find 22 matches to known pulsars, including one not reported elsewhere, and 8 new possible companions to young pulsars. We examine the photometric and kinematic properties of these systems and provide empirical relations for identifying Gaia sources with potential millisecond pulsar companions. Our results confirm that the observed multiplicity fraction is small. However, we show that the number of binaries below the sensitivity of Gaia and radio timing in our sample could still be significantly higher. We constrain the binary fraction of young pulsars to be $f_{\rm young}^{\rm true}\leq 5.3(8.3)\%$ under realistic(conservative) assumptions for the binary properties and current sensitivity thresholds. For massive stars ($\geq 10$ M$_{\odot}$) in particular, we find $f_{\rm OB}^{\rm true}\leq 3.7\%$ which sets a firm independent upper limit on the galactic neutron-star merger rate, $\leq 7.2\times 10^{-4}$ yr$^{-1}$. Ongoing and future projects such as the CHIME/pulsar program, MeerTime, HIRAX and ultimately the SKA, will significantly improve these constraints in the future., Comment: Accepted for publication in MNRAS; comments are welcome

We modelled the Chandra and RXTE X-ray spectra of the massive binary WR 140 in the framework of the standard colliding stellar wind (CSW) picture. Models with partial electron heating at the shock fronts are a better representation of the X-ray data than those with complete temperature equalization. Emission measure of the X-ray plasma in the CSW region exhibits a considerable decrease at orbital phases near periastron. This is equivalent to variable effective mass-loss rates over the binary orbit. At orbital phases near periastron, a considerable X-ray absorption in excess to that from the stellar winds in WR 140 is present. The standard CSW model provides line profiles that in general do not match well the observed line profiles of the strong line features in the X-ray spectrum of WR 140. The variable effective mass-loss rate could be understood qualitatively in CSW picture of clumpy stellar winds where clumps are efficiently dissolved in the CSW region near apastron but not at periastron. However, future development of CSW models with non-spherically-symmetric stellar winds might be needed to get a better correspondence between theory and observations., Accepted for publication in MNRAS; 13 pages, 10 figures, 4 tables

We use a combined binary evolution code including dynamical effects to study nova eruptions in a symbiotic system. Following the evolution, over $\sim10^5$ years, of multiple consecutive nova eruptions on the surface of a $1.25M_\odot$ white dwarf (WD) accretor, we present a comparison between simulations of two types of systems. The first is the common, well known, cataclysmic variable (CV) system in which a main sequence donor star transfers mass to its WD companion via Roche-lobe overflow. The second is a detached, widely separated, symbiotic system in which an asymptotic giant branch donor star transfers mass to its WD companion via strong winds. For the latter we use the Bondi-Hoyle-Lyttleton prescription along with orbital dynamics to calculate the accretion rate. We use the combined stellar evolution code to follow the nova eruptions of both simulations including changes in mass, accretion rate and orbital features. We find that while the average accretion rate for the CV remains fairly constant, the symbiotic system experiences distinct epochs of high and low accretion rates. The examination of epochs for which the accretion rates of both simulations are similar, shows that the evolutionary behaviors are identical. We obtain that for a given WD mass, the rate that mass is accreted ultimately determines the development, and that the stellar class of the donor is of no significance to the development of novae. We discuss several observed systems and find that our results are consistent with estimated parameters of novae in widely separated symbiotic systems., Comment: Accepted for publication in MNRAS

Fast radio bursts (FRBs) are millisecond-duration radio transients from cosmological distances. Their isotropic energies follow a power-law distribution with a possible exponential cutoff, but their intrinsic redshift distribution, which contains information about the FRB sources, is not well understood. We attempt to constrain both distributions by means of Monte Carlo simulations and comparing the simulations results with the available FRB specific fluence distribution, dispersion measure (DM) distribution, and the estimated energy distribution data. Two redshift distribution models, one tracking the star formation history of the universe and another tracking compact binary mergers, are tested. For the latter model, we consider three merger delay timescale distribution (Gaussian, log-normal, and power law) models. Two FRB samples detected by Parkes and the Australian Square Kilometre Array Pathfinder (ASKAP), respectively, are used to confront the simulation results. We confirm the $\sim -1.8$ power law index for the energy distribution but the exponential cutoff energy of the distribution, if any, is unconstrained. For the best energy distribution model, none of the redshift distributions we considered are rejected by the data. A future, larger, uniform FRB sample (such as the one collected by the Canadian Hydrogen Intensity Mapping Experiment, CHIME) can provide better constraints on the intrinsic FRB redshift distribution using the methodology presented in this paper., MNRAS, in press