Search

Binary stars play a crucial role in our understanding of the formation and evolution of star clusters and their stellar populations. We use Gaia Data Release 3 to homogeneously analyze 78 Galactic open clusters and the unresolved binary systems they host, each composed of two main sequence (MS) stars. We first investigated the structural parameters of these clusters, such as the core radius and the central density, and determined the cluster mass function (MF) and total mass by interpolating the density profile of each cluster. We measured the fraction of binaries with a large mass ratio and the fraction of blue straggler stars (BSSs), and finally investigated possible connections between the populations of binary stars and BSSs with the main parameters of the host cluster. {Remarkably, we find that the MFs of 78 analyzed open clusters follow a similar trend and are well reproduced by two single power-law functions, with a change in slope around masses of 1$M_{\odot}$. The fraction of binary stars ranges from $\sim$15\% to more than $\sim$60\% without significant correlation with the mass and the age of the host cluster. Moreover, we detect hints of a correlation between the total fraction of binary stars and the central density of the host cluster. We compared the fraction of binary stars with that of BSSs, finding that clusters with high and low central density exhibit different trends. The fraction of binaries does not significantly change with the mass of the primary star and the mass ratio. The radial distribution of binary stars depends on cluster age. The binaries of clusters younger than $\sim$800\,Myr typically show a flat radial distribution, with some hints of a double peak. In contrast, the binaries of the remaining clusters are more centrally concentrated than the single stars, which is similar to what is observed in globular clusters., Comment: 27 pages, 27 figures, accepted for publication in Astronomy & Astrophysics

More than 100 millisecond pulsars (MSPs) have been discovered in radio observations of gamma-ray sources detected by the Fermi Large Area Telescope (LAT), but hundreds of pulsar-like sources remain unidentified. Here we present the first results from the targeted survey of Fermi-LAT sources being performed by the Transients and Pulsars with MeerKAT (TRAPUM) Large Survey Project. We observed 79 sources identified as possible gamma-ray pulsar candidates by a Random Forest classification of unassociated sources from the 4FGL catalogue. Each source was observed for 10 minutes on two separate epochs using MeerKAT's L-band receiver (856-1712 MHz), with typical pulsed flux density sensitivities of $\sim$100$\,\mu$Jy. Nine new MSPs were discovered, eight of which are in binary systems, including two eclipsing redbacks and one system, PSR J1526$-$2744, that appears to have a white dwarf companion in an unusually compact 5 hr orbit. We obtained phase-connected timing solutions for two of these MSPs, enabling the detection of gamma-ray pulsations in the Fermi-LAT data. A follow-up search for continuous gravitational waves from PSR J1526$-$2744 in Advanced LIGO data using the resulting Fermi-LAT timing ephemeris yielded no detection, but sets an upper limit on the neutron star ellipticity of $2.45\times10^{-8}$. We also detected X-ray emission from the redback PSR J1803$-$6707 in data from the first eROSITA all-sky survey, likely due to emission from an intra-binary shock., Comment: 17 pages, 8 figures, accepted for publication in MNRAS

M33 X-7 is the only known eclipsing black hole high mass X-ray binary. The system is reported to contain a very massive O supergiant donor and a massive black hole in a short orbit. The high X-ray luminosity and its location in the metal-poor galaxy M33 make it a unique laboratory for studying the winds of metal-poor donor stars with black hole companions and it helps us to understand the potential progenitors of black hole mergers. Using phase-resolved simultaneous HST- and XMM-Newton-observations, we traced the interaction of the stellar wind with the black hole. Our comprehensive spectroscopic investigation of the donor star (X-ray+UV+optical) yields new stellar and wind parameters for the system that differs significantly from previous estimates. In particular, the masses of the components are considerably reduced to 38 for the O-star donor and 11.4 for the black hole. The O giant is overfilling its Roche lobe and shows surface He enrichment. The donor shows a densely clumped wind with a mass-loss rate that matches theoretical predictions. We investigated the wind-driving contributions from different ions and the changes in the ionization structure due to X-ray illumination. Toward the black hole, the wind is strongly quenched due to strong X-ray illumination. For this system, the standard wind-fed accretion scenario alone cannot explain the observed X-ray luminosity, pointing toward an additional mass overflow, which is in line with our acceleration calculations. The X-ray photoionization creates an He II emission region emitting $10^{47}$ ph/s. We computed binary evolutionary tracks for the system using MESA. Currently, the system is transitioning toward an unstable mass transfer phase, resulting in a common envelope of the black hole and donor. Since the mass ratio is q~3.3 and the period is short, the system is unlikely to survive the common envelope, but will rather merge., Accepted for publication in Astronomy & Astrophysics

Without doubt, mass transfer in close binary systems contributes to the populations of Wolf-Rayet (WR) stars in the Milky Way and the Magellanic Clouds. However, the binary formation channel is so far not well explored. We want to remedy this by exploring large grids of detailed binary and single star evolution models computed with the publicly available MESA code, for a metallicity appropriate for the Large Magellanic Cloud (LMC). The binary models are calculated through Roche-lobe overflow and mass transfer, until the initially more massive star exhausts helium in its core. We distinguish models of WR and helium stars based on the estimated stellar wind optical depth. We use these models to build a synthetic WR population, assuming constant star formation. Our models can reproduce the WR population of the LMC to significant detail, including the number and luminosity functions of the main WR subtypes. We find that for binary fractions of 100% (50%), all LMC WR stars below $10^6\,L_{\odot}$ ($10^{5.7}\,L_{\odot}$) are stripped binary mass donors. We also identify several insightful mismatches. With a single star fraction of 50\%, our models produce too many yellow supergiants, calling either for a larger initial binary fraction, or for enhanced mass-loss near the Humphreys-Davidson limit. Our models predict more long-period WR binaries than observed, arguably due to an observational bias towards short periods. Our models also underpredict the shortest-period WR binaries, which may have implications for understanding the progenitors of double black hole mergers. The fraction of binary produced WR stars may be larger than often assumed, and outline the risk to mis-calibrate stellar physics when only single star models are used to reproduce the observed WR stars., 15 pages + 13 pages appendix, 14 figures, 2 tables

We present the analysis of the optical variability of the early, nitrogen-rich Wolf-Rayet (WR) star WR7. The analysis of multi-sector Transiting Exoplanet Survey Satellite (TESS) light curves and high-resolution spectroscopic observations confirm multi-periodic variability that is modulated on time-scales of years. We detect a dominant period of $2.6433 \pm 0.0005$ d in the TESS sectors 33 and 34 light curves in addition to the previously reported high-frequency features from sector 7. We discuss the plausible mechanisms that may be responsible for such variability in WR7, including pulsations, binarity, co-rotating interacting regions (CIRs) and clumpy winds. Given the lack of strong evidence for the presence of a stellar or compact companion, we suggest that WR7 may pulsate in quasi-coherent modes in addition to wind variability likely caused by CIRs on top of stochastic low-frequency variability. WR7 is certainly a worthy target for future monitoring in both spectroscopy and photometry to sample both the short ($\lesssim 1$ d) and long ($\gtrsim 1000$ d) variability time scales., 9 pages, 5 Figures, 1 Table; Accepted to MNRAS

New CCD photometric observations and their investigation of the W UMa-type binary, V870 Ara, are presented. Light curves of the system were taken through BVI filters from the Congarinni Observatory in Australia. The new ephemeris is calculated based on seven new determined minimum times, together with the TESS data and others compiled from the literature. Photometric solutions determined by the Wilson-Devinney (W-D) code are combined with the Monte Carlo simulation to determine the adjustable parameters' uncertainties. These solutions suggest that V870 Ara is a contact binary system with a mass ratio of 0.082, a fillout factor of 96+-4 percent, and an inclination of 73.60+-0.64 degrees. The absolute parameters of V870 Ara were determined by combining the Gaia EDR3 parallax and photometric elements., Comment: 5 tables, 6 figs, accepted in the Open Astronomy Journal

About a quarter of massive binary stars undergo mass transfer while both stars burn hydrogen at their cores, first on the thermal and then on the nuclear timescale. The nuclear timescale mass transfer leads to observable counterparts: the semi-detached so-called massive Algol binaries. However, comprehensive model predictions for these systems are sparse. We study them using a large grid of ~10,000 detailed binary evolution models calculated with the stellar evolution code MESA, covering initial donor masses between 10-40 M$_{\odot}$ and initial orbital periods above 1.4 d, at a metallicity suitable for the Large Magellanic Cloud (LMC). Our models imply ~30, or ~3% of the ~1,000 core hydrogen burning O-star binaries in the LMC to be currently in the semi-detached phase. Our donor models are up to 25-times more luminous than single stars of identical mass and effective temperature, which agrees with the observed Algols. A comparison of our models with the observed orbital periods and mass ratios implies rather conservative mass transfer in some systems, while very inefficient one in others. This is generally well reproduced by our spin-dependent mass transfer algorithm, except for the lowest considered masses. The observations reflect the slow increase of the surface nitrogen enrichment of the donors during the semi-detached phase all the way to CNO equilibrium. We also investigate the properties of our models after core hydrogen depletion of the donor star, when these models correspond to Wolf-Rayet/helium+OB star binaries. A dedicated spectroscopic survey of massive Algol systems may allow to derive the dependence of the efficiency of thermal timescale mass transfer on the binary parameters, as well as the efficiency of semiconvective mixing in the stellar interior. This would be a crucial step towards reliable binary models up to the formation of supernovae and compact objects., Accepted for publication in A&A (December 6, 2021)

Context. Many well-known bright stars have been observed by the ongoing Transiting Exoplanet Survey Satellite (TESS) space mission. For several of them, these new data reveal previously unobserved variability, such as tidally perturbed pulsations in close binary stars. Aims. Using newly detected gravity-mode (g-mode) pulsations in V456 Cyg, we aim to determine the global stellar properties of this short-period eclipsing binary and evaluate the interaction between these pulsations and the tides. Methods. We model the binary orbit and determine the physical properties of the component stars using the TESS photometry and published spectroscopy. We then measure the pulsation frequencies from the residuals of the light curve fit using iterative prewhitening, and analyse them to determine the global asteroseismic stellar parameters. We evaluate the pulsation parameters as a function of the orbital phase. Results. We find that the pulsations belong to the secondary component of V456 Cyg and that this star likely has a uniform radial rotation profile, synchronous ($\nu_{\rm rot} = 1.113(14) \rm d^{-1}$) with the binary orbit ($\nu_{\rm orb} = 1.122091(8) \rm d^{-1}$). The observed g modes are amplified by almost a factor three in the stellar hemisphere facing the primary. We present evidence that this is caused by tidal perturbation of the pulsations, with the mode coupling being strongly affected. Conclusions. V456 Cyg is only the second object for which tidally perturbed high-order g-mode pulsations are identified, after $\pi^5$ Ori. This opens up new opportunities for tidal g-mode asteroseismology, as it demonstrates another avenue in which g modes and tides can interact with each other., Comment: 8 pages, 6 figures, 3 tables. Accepted for publication in Astronomy & Astrophysics on 12 January 2022

Context. Young multiple stellar systems can host both circumstellar and circumbinary discs composed of gas and dust, and the orientations of circumbinary discs can be sculpted by the orientation and eccentricity of the central binaries. Studying multiple binary systems and their associated discs enables our understanding of the size and distribution of the planetary systems that subsequently form around them. Aims. A deep (~70%) and extended (~150 days) eclipse was seen towards the young multiple stellar system V773 Tau in 2010. We interpret it as being due to the passage of a circumbinary disc around the B components moving in front of the A components. Our aim is to characterise the orientation and structure of the disc, to refine the orbits of the sub-components, and to predict when the next eclipse will occur. Methods. We combined the photometry from several ground-based surveys, constructed a model for the light curve of the eclipse, and used high angular resolution imaging to refine the orbits of the three spatially resolved components of the system: A, B, and C. A frequency analysis of the light curves, including from the TESS satellite, enabled the characterisation of the rotational periods of the Aa and Ab stars. Results. A toy model of the circumbinary disc shows that it extends out to approximately 5 au around the B binary and has an inclination of 73° with respect to the orbital plane of AB, where the lower bound of the radius of the disc is constrained by the geometry of the AB orbit and the upper bound is set by the stability of the disc. We identify several frequencies in the photometric data that we attribute to rotational modulation of the Aa and Ab stellar companions. We produced the first determination of the orbit of the more distant C component around the AB system and limited its inclination to 93°. Conclusions. The high inclination and large diameter of the disc, together with the expected inclination of the disc from theory, suggest that B is an almost equal-mass, moderately eccentric binary. We identify the rotational periods of the Aa and Ab stars and a third frequency in the light curve that we attribute to the orbital period of the stars in the B binary. We predict that the next eclipse will occur around 2037, during which both detailed photometric and spectroscopic monitoring will characterise the disc in greater detail.

We investigate the aperiodic variability for a relatively large sample of accreting neutron stars and intermediate polars, focusing on the properties of the characteristic break commonly observed in power spectra of accreting objects. In particular, we investigate the relation of the break frequency and the magnetic field strength, both of which are connected to the size of the magnetosphere. We find that for the majority of objects in our sample the measured break frequency values indeed agree with estimated inner radii of the accretion disc, which allows to use observed break frequencies to independently assess the magnetic field strength and structure in accreting compact objects. As a special case, we focus on Hercules X-1 which is a persistent, medium-luminosity X-ray pulsar accreting from its low-mass companion. In the literature, it has been suggested that the complex pulse profiles, the spin-up behaviour and the luminosity-correlation of the cyclotron energy seen in Her X-1 can be explained with a complex magnetic field structure of the neutron star. Here, we connect the measured break frequency to the magnetospheric radius and show that the magnetic field strength derived assuming a dipole configuration is nearly an order of magnitude smaller than the magnetic field strength corresponding to the cyclotron energy. Accordingly, this discrepancy can be explained with the magnetic field having strong multipole components. The multipolar structure would also increase the accreting area on the neutron star surface, explaining why the critical luminosity for accretion column formation is puzzlingly high in this source., Comment: 10 pages, 3 figures

Accurate measurements of the masses of neutron stars are necessary to test binary evolution models, and to constrain the neutron star equation of state. In pulsar binaries with no measurable post-Keplerian parameters, this requires an accurate estimate of the binary system's inclination and the radial velocity of the companion star by other means than pulsar timing. In this paper, we present the results of a new method for measuring this radial velocity using the binary synthesis code Icarus. This method relies on constructing a model spectrum of a tidally distorted, irradiated star as viewed for a given binary configuration. This method is applied to optical spectra of the newly discovered black widow PSR J1555-2908. By modelling the optical spectroscopy alongside optical photometry, we find that the radial velocity of the companion star is $397\pm4$ km s$^{-1}$ (errors quoted at 95\% confidence interval), as well as a binary inclination of $>75^{\rm o}$. Combined with $\gamma$-ray pulsation timing information, this gives a neutron star mass of 1.67$^{+0.15}_{-0.09}$ M$_\odot$ and a companion mass of 0.060$^{+0.005}_{-0.003}$ M$_\odot$, placing PSR J1555-2908 at the observed upper limit of what is considered a black widow system., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 15 pages, 7 Figures. Underlying data available at https://zenodo.org/record/5653061

In the realm of massive stars, strong binary interaction is commonplace. One extreme case are overcontact systems, which is expected to be part of the evolution of all stars evolving towards a merger, and is hypothesized to play a role in the formation of binary black holes. However, important simplifications are made to model the evolution of overcontact binaries, namely that tidal or rotational deformation is frequently ignored. Yet, both observation and theory show that overcontact stars are heavily tidally deformed, leaving a potentially important effect on the outer layers unaccounted for in models. In this work we develop the methodology to represented tidally deformed stars. Using numerical methods, we compute the structure correction factors to the 1D spherical stellar structure equations due to the binary Roche potential, and compare them to existing results and the structure corrections of single rotating stars. We implement the new structure correction factors into the stellar evolution code MESA and explore several case studies. We compare the differences between our simulations when no rotation is included, when we treat rotation using single star corrections (i.e. only accounting for centrifugal deformation) or when we use tidal deformation. We find that ignoring rotation in deformed detached eclipsing binaries can produce a radius discrepancy of up to 5%. The difference between tidal and single star centrifugal distortion models is more benign at 1%, showing that single rotating star models are a suitable approximation of tidally deformed stars in a binary system. In overcontact configurations, we find a similar 5% variation in surface properties as a result of tidal distortion with respect to non-rotating models, showing that it is inappropriate to model binary stars that fill their Roche lobe significantly, as non-rotating., Comment: 12 pages, 9 figures, accepted to appear in Astronomy & Astrophysics

Given that mergers are often invoked to explain many exotic phenomena in massive star evolution, understanding the evolutionary phase directly preceding a merger, the overcontact phase, is of crucial importance. Despite its importance, large uncertainties exist in our understanding of the evolution of massive overcontact binaries. We aim to provide robust observational constraints on the future dynamical evolution of massive overcontact systems by measuring the rate at which the periods change for a sample of six such objects. Furthermore, we aim to investigate whether the periods of unequal mass systems show higher rates of change than their equal mass counterparts as theoretical models predict. Using archival photometric data from various ground- and space-based missions covering up to ~40 years, we measure the periods of each system over several smaller time spans. We then fit a linear regression through the measured periods to determine the rate at which the period is changing over the entire data set. We find that all of the stars in our sample have very small period changes and that there does not seem to be a correlation with the mass ratio. This implies that the orbital periods for these systems are stable on the nuclear timescale, and that the unequal mass systems may not equalize as expected. When comparing our results with population synthesis distributions, we find large discrepancies between the expected mass ratios and period stabilities. We find that these discrepancies can be mitigated to a degree by removing systems with shorter initial periods, suggesting that the observed sample of overcontact systems may originate from binary systems with longer initial orbital periods., Comment: 11 pages, 6 (+5) Figures, accepted for publication in A&A

Abbreviated. We aim to identify new low-mass companions to young stars using the astrometric measurements provided by the Gaia space mission and complementary VLT/SPHERE data. We identify companion candidates from a sample of K-type, pre-main sequence stars in the Scorpius Centaurus association using the early version of the third data release of the Gaia space mission. Based on the provided positions, proper motions, and magnitudes, we identify all objects within a predefined radius whose differential proper motions are consistent with a gravitationally bound system. We derive companion masses through comparison with evolutionary tracks. For seven identified companion candidates we use additional data collected with VLT/SPHERE and VLT/NACO to assess the accuracy of the properties of the companions based on Gaia photometry alone. We identify 110 comoving companions that have a companionship likelihood of more than $95\,\%$. We identify ten especially intriguing companions that have masses in the brown dwarf regime down to $20\,M_\mathrm{Jup}$. Our high-contrast imaging data confirm both astrometry and photometric masses derived from Gaia alone. We discover a new brown dwarf companion, TYC 8252-533-1 B, with a projected separation of approximately $570\,\mathrm{au}$ from its Sun-like primary. SED modeling provides a companion mass of $52^{+17}_{-11}\,M_\mathrm{Jup}$. We show that the Gaia database can identify low-mass companions at wide separations from their host stars. For K-type Sco-Cen members Gaia can detect sub-stellar objects at projected separations larger than $300\,\mathrm{au}$ and is sensitivity limited beyond $1,000\,\mathrm{au}$ with a lower mass limit down to $20\,M_\mathrm{Jup}$. A similar analysis of other star-forming regions could significantly enlarge the sample size of such objects and test formation and evolution theories of planetary systems., Accepted for publication in A&A (31 pages, 20 figures, 12 tables)

We present for the neutron-star low-mass X-ray binary 4U 1636–53, and for the first time for any source of kilohertz quasi-periodic oscillations (kHz QPOs), the two-dimensional behaviour of the fractional rms amplitude of the kHz QPOs in the parameter space defined by QPO frequency and photon energy. We find that the rms amplitude of the lower kHz QPO increases with energy up to $\sim \!\!12$ keV and then decreases at higher energies, while the rms amplitude of the upper kHz QPO either continues increasing or levels off at high energies. The rms amplitude of the lower kHz QPO increases and then decreases with frequency, peaking at $\sim \!\!760$ Hz, while the amplitude of the upper kHz QPO decreases with frequency, with a local maximum at around $\sim \!\!770$ Hz, and is consistent with becoming zero at the same QPO frequency, $\sim \!\!1400$ Hz, in all energy bands, thus constraining the neutron-star mass at $M_{NS} \le 1.6 \, \mathrm{M}_{\odot }$, under the assumption that this QPO reflects the Keplerian frequency at the inner edge of the accretion disc. We show that the slope of the rms energy spectrum is connected to the changing properties of the kHz QPOs in different energy bands as its frequencies change. Finally, we discuss a possible mechanism responsible for the radiative properties of the kHz QPOs and, based on a model in which the QPO arises from oscillations in a Comptonizing cloud of hot electrons, we show that the properties of the kHz QPOs can constrain the thermodynamic properties of the inner accretion flow.

We performed a detailed spectroscopic study of the SB2 system HD 60803 based on high-resolution spectra obtained with the different spectrographs. The analysis was done with two independent methods: a) the direct modelling of the observed binary spectrum by a sum of synthetic spectra varying a set of free parameters and minimizing a difference between the observed and theoretical spectra; b) spectrum disentangling and an independent modelling of the individual components. Being applied to binary spectra from different spectrographs both methods converge to a consistent solution for the fundamental parameters of the HD 60803 components: $T_{\rm eff}$=6055$\pm$70 K, $\log{g}$=4.08$\pm$0.12, $\zeta_{\rm RT}$=1.45$\pm$0.18 km s$^{-1}$, [M/H]=0.03$\pm$0.06 (primary), and $T_{\rm eff}$=6069$\pm$70 K, $\log{g}$=4.14$\pm$0.09, $\zeta_{\rm RT}$=1.48$\pm$0.18 km s$^{-1}$, [M/H]=0.03$\pm$0.06 (secondary). Differential abundance analysis of the components did not reveal any significant difference in their chemical composition. Besides Li both components have solar atmospheric abundances. Li abundance exceeds the solar one by $\sim$2 dex, but it agrees with Li abundance in main-sequence late F-stars. Relative-to-solar abundances in both components slightly correlate with the condensation temperature the same way as was found in the solar analogs with/without detected giant planets. The estimated age of the system is 5.5$\pm$0.5 Gyr., Comment: 10 pages, 8 figures, 4 tables, accepted in MNRAS

The accuracy of theoretical mass, radius and effective temperature values for M-dwarf stars is an active topic of debate. Differences between observed and theoretical values have raised the possibility that current theoretical stellar structure and evolution models are inaccurate towards the low-mass end of the main sequence. To explore this issue we use the CHEOPS satellite to obtain high-precision light curves of eclipsing binaries with low mass stellar companions. We use these light curves combined with the spectroscopic orbit for the solar-type companion to measure the mass, radius and effective temperature of the M-dwarf star. Here we present the analysis of three eclipsing binaries. We use the pycheops data analysis software to fit the observed transit and eclipse events of each system. Two of our systems were also observed by the TESS satellite -- we similarly analyse these light curves for comparison. We find consistent results between CHEOPS and TESS, presenting three stellar radii and two stellar effective temperature values of low-mass stellar objects. These initial results from our on-going observing programme with CHEOPS show that we can expect to have ~24 new mass, radius and effective temperature measurements for very low mass stars within the next few years., Comment: 12 pages, 8 figures, accepted for publication in MNRAS

Context. Eta Carinae is a highly eccentric, massive binary system (semimajor axis ~15.5 au) with powerful stellar winds and a phase-dependent wind-wind collision (WWC) zone. The primary star, η Car A, is a luminous blue variable (LBV); the secondary, η Car B, is a Wolf-Rayet or O star with a faster but less dense wind. Aperture-synthesis imaging allows us to study the mass loss from the enigmatic LBV η Car. Understanding LBVs is a crucial step toward improving our knowledge about massive stars and their evolution. Aims. Our aim is to study the intensity distribution and kinematics of η Car’s WWC zone. Methods. Using the VLTI-MATISSE mid-infrared interferometry instrument, we perform Brα imaging of η Car’s distorted wind. Results. We present the first VLTI-MATISSE aperture-synthesis images of η Car A’s stellar windin several spectral channels distributed across the Brα 4.052 μm line (spectral resolving power R ~ 960). Our observations were performed close to periastron passage in February 2020 (orbital phase ~ 14.0022). The reconstructed iso-velocity images show the dependence of the primary stellar wind on wavelength or line-of-sight (LOS) velocity with a spatial resolution of 6 mas (~14 au). The radius of the faintest outer wind regions is ~26 mas (~60 au). At several negative LOS velocities, the primary stellar wind is less extended to the northwest than in other directions. This asymmetry is most likely caused by the WWC. Therefore, we see both the velocity field of the undisturbed primary wind and the WWC cavity. In continuum spectral channels, the primary star wind is more compact than in line channels. A fit of the observed continuum visibilities with the visibilities of a stellar wind CMFGEN model (CMFGEN is an atmosphere code developed to model the spectra of a variety of objects) provides a full width at half maximum fit diameter of the primary stellar wind of 2.84 ± 0.06 mas (6.54 ± 0.14 au). We comparethe derived intensity distributions with the CMFGEN stellar wind model and hydrodynamic WWC models.

Coherent radio emission from stars can be used to constrain fundamental coronal plasma parameters, such as plasma density and magnetic field strength. It is also a probe of chromospheric and magnetospheric acceleration mechanisms. Close stellar binaries, such as RS Canum Venaticorum (RS CVn) systems, are particularly interesting as their heightened level of chromospheric activity and possible direct magnetic interaction make them a unique laboratory to study coronal and magnetospheric acceleration mechanisms. RS CVn binaries are known to be radio-bright but coherent radio emission has only conclusively been detected previously in one system. Here, we present a population of 14 coherent radio emitting RS CVn systems. We identified the population in the ongoing LOFAR Two Metre Sky Survey as circularly polarised sources at 144MHz that are astrometrically associated with known RS CVn binaries. We show that the observed emission is powered by the electron cyclotron maser instability. We use numerical calculations of the maser's beaming geometry to argue that the commonly invoked 'loss-cone' maser cannot generate the necessary brightness temperature in some of our detections and that a more efficient instability, such as the shell or horseshoe maser, must be invoked. Such maser configurations are known to operate in planetary magnetospheres. We also outline two acceleration mechanisms that could produce coherent radio emission, one where the acceleration occurs in the chromosphere and one where the acceleration is due to an electrodynamic interaction between the stars. We propose radio and optical monitoring observations that can differentiate between these two mechanisms., This article consists of 14 pages, 10 normal figures and 4 gifs. The gifs themselves each consist of 32 figures. This article is to be published in Astronomy and Astrophysics, section 8. 'Stellar Atmospheres'

We report new radio observations of a sample of thirty-six neutron star (NS) X-ray binaries, more than doubling the sample in the literature observed at current-day sensitivities. These sources include thirteen weakly-magnetised ($B, Accepted for publication in MNRAS; 31 pages, 8 figures, 8 tables (including online supplementary materials)

We perform a search for rapidly oscillating Ap stars in the Kepler long-cadence data, where true oscillations above the Nyquist limit of 283.21 {\mu}Hz can be reliably distinguished from aliases as a consequence of the barycentric time corrections applied to the Kepler data. We find evidence for rapid oscillations in six stars: KIC 6631188, KIC 7018170, KIC 10685175, KIC 11031749, KIC 11296437 and KIC 11409673, and identify each star as chemically peculiar through either pre-existing classifications or spectroscopic measurements. For each star, we identify the principal pulsation mode, and are able to observe several additional pulsation modes in KIC 7018170. We find that KIC 7018170 and KIC 11409673 both oscillate above their theoretical acoustic cutoff frequency, whilst KIC 11031749 oscillates at the cutoff frequency within uncertainty. All but KIC 11031749 exhibit strong amplitude modulation consistent with the oblique pulsator model, confirming their mode geometry and periods of rotation., Comment: 20 pages, 14 figures, accepted for publication in MNRAS

We present the broadband spectral analysis of all the six hard, intermediate and soft state NuSTAR observations of the recently discovered transient black hole X-ray binary MAXI J1348-630 during its first outburst in 2019. We first model the data with a combination of a multi-colour disc and a relativistic blurred reflection, and, whenever needed, a distant reflection. We find that this simple model scheme is inadequate in explaining the spectra, resulting in a very high iron abundance. We, therefore, explore the possibility of reflection from a high-density disc. We use two different sets of models to describe the high-density disc reflection: relxill-based reflection models, and reflionx-based ones. The reflionx-based high-density disc reflection models bring down the iron abundance to around the solar value, while the density is found to be $10^{20.3-21.4} \rm cm^{-3}$. We also find evidence of a high-velocity outflow in the form of $\sim$7.3 keV absorption lines. The consistency between the best-fit parameters for different epochs and the statistical significance of the corresponding model indicates the existence of high-density disc reflection in MAXI J1348-630., 15 pages, 12 figures, 5 tables; accepted for publication in MNRAS

New observational data obtained during the 2018, 2019, and 2020 speckle runs with the 4.1-m Southern Astrophysical Research (SOAR) telescope located at Cerro Pachón (Chile) allowed us to recalculate the orbits of the visual binaries: WDS 00277-1625 (YR 1 Aa,Ab), WDS 00462-2214 (RST4155), WDS 03124-4425 (JC 8), WDS 07427-3510 (HDS 1091), WDS 10093+2020 (A2145), WDS 10116+1321 (HU 874 AB), WDS 10217-0946 (BU 25), WDS 11585-2350 (RST 3767 AB), WDS 13117-2633 (FIN 305), WDS 13305+0729 (A 1789), and WDS 16458-0046 (A1141). In addition, we present the first orbits calculated for the following binaries: WDS 07303-5657 (FIN 105), WDS 09110-1929 (I 824), WDS 12111-5302 (HU 1604), WDS 14592-4206 (HDS 2116 Aa,Ab), WDS 15157-2736 (BU 350), WDS 15493+0503 (A 1126), WDS 16402-2800 (VOU 44 AB), and WDS 18126+1224 (HDS 2570). All of them are placed below Dec. + 21° and the majority are main-sequence stars. Except in cases with giant components or close triple systems, the ANAPAR method was used in order to obtain precise dynamical parallaxes and individual masses. These parallaxes were compared with those obtained by Gaia and/or Hipparcos satellites. In the case of FIN 305, we present two different orbital solutions. Also, using the dynamical parallaxes given by these orbits, we have been able to calculate the luminosity of these systems. Said luminosities allow us to indicate an approximate age for the components of these systems, situating them within the Hertzsprung-Russell (HR) diagram. In addition, a commentary for each binary about the physical and dynamical properties of the studied binaries has been included We are very grateful for the continuous support of the Chilean Telescope Allocation Committee: http://www.das.uchile.cl/das_cntac.html under programs CN2018A-1, CN2019A-2, CN2019B-13, and CN2020A-19 and 2020B-10. This paper was supported by the ‘Xunta de Galicia’ (Spain), under the ED431B 2020/38 grant. RAM and EC acknowledge support from ANID/FONDECYT grant 1190038 SI

Abbott, R., et al. (LIGO and Virgo Collaboration), We search for signatures of gravitational lensing in the gravitational-wave signals from compact binary coalescences detected by Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) and Advanced Virgo during O3a, the first half of their third observing run. We study: (1) the expected rate of lensing at current detector sensitivity and the implications of a non-observation of strong lensing or a stochastic gravitational-wave background on the merger-rate density at high redshift; (2) how the interpretation of individual high-mass events would change if they were found to be lensed; (3) the possibility of multiple images due to strong lensing by galaxies or galaxy clusters; and (4) possible wave-optics effects due to point-mass microlenses. Several pairs of signals in the multiple-image analysis show similar parameters and, in this sense, are nominally consistent with the strong lensing hypothesis. However, taking into account population priors, selection effects, and the prior odds against lensing, these events do not provide sufficient evidence for lensing. Overall, we find no compelling evidence for lensing in the observed gravitational-wave signals from any of these analyses., This material is based upon work supported by National Science Foundationʼs LIGO Laboratory, which is a major facility fully funded by the National Science Foundation. The authors also gratefully acknowledge the support of the Science and Technology Facilities Council (STFC) of the United Kingdom, the Max-Planck-Society (MPS), and the State of Niedersachsen/Germany for support of the construction of Advanced LIGO and construction and operation of the GEO600 detector. Additional support for Advanced LIGO was provided by the Australian Research Council. The authors gratefully acknowledge the Italian Istituto Nazionale di Fisica Nucleare (INFN), the French Centre National de la Recherche Scientifique (CNRS), and the Netherlands Organization for Scientific Research for the construction and operation of the Virgo detector and the creation and support of the European Gravitational Observatory (EGO) consortium. The authors also gratefully acknowledge research support from these agencies as well as by the Council of Scientific and Industrial Research of India; the Department of Science and Technology, India; the Science & Engineering Research Board (SERB), India; the Ministry of Human Resource Development, India; the Spanish Agencia Estatal de Investigación; the Vicepresidència i Conselleria d’Innovació Recerca i Turisme and the Conselleria d’Educació i Universitat del Govern de les Illes Balears; the Conselleria d’Innovació Universitats; Ciència i Societat Digital de la Generalitat Valenciana and the CERCA Programme Generalitat de Catalunya, Spain; the National Science Centre of Poland and the Foundation for Polish Science (FNP); the Swiss National Science Foundation (SNSF); the Russian Foundation for Basic Research; the Russian Science Foundation; the European Commission; the European Regional Development Funds (ERDF); the Royal Society; the Scottish Funding Council; the Scottish Universities Physics Alliance; the Hungarian Scientific Research Fund (OTKA); the French Lyon Institute of Origins (LIO); the Belgian Fonds de la Recherche Scientifique (FRS-FNRS); Actions de Recherche Concertées (ARC) and Fonds Wetenschappelijk Onderzoek—Vlaanderen (FWO), Belgium; the Paris Î le-de-France Region; the National Research, Development and Innovation Office Hungary (NKFIH); the National Research Foundation of Korea; the Natural Science and Engineering Research Council Canada; Canadian Foundation for Innovation (CFI); the Brazilian Ministry of Science, Technology, and Innovations; the International Center for Theoretical Physics South American Institute for Fundamental Research (ICTP-SAIFR); the Research Grants Council of Hong Kong; the National Natural Science Foundation of China (NSFC); the Leverhulme Trust; the Research Corporation; the Ministry of Science and Technology (MOST), Taiwan; the United States Department of Energy; the Kavli Foundation; and the Gordon and Betty Moore Foundation. The authors gratefully acknowledge the support of the NSF, STFC, INFN, and CNRS for provisions of the computational resources.

When a core-collapse supernova explodes in a binary star system, the ejecta might encounter an overdense shell, where the stellar winds of the two stars previously collided. In this work, we investigate effects of such interactions on supernova light curves on time-scales from the early flash ionization signatures to approximately one year after the explosion. We construct a model of the colliding-wind shell in an orbiting binary star system and we provide an analytical expression for the shell thickness and density, which we calibrate with three-dimensional adaptive mesh refinement hydrodynamical simulations probing different ratios of wind momenta and different regimes of radiative cooling efficiency. We model the angle-dependent interaction of supernova ejecta with the circumstellar medium and estimate the shock radiative efficiency with a realistic cooling function. We find that the radiated shock power exceeds typical Type IIP supernova luminosity only for double red supergiant binaries with mass ratios $q \gtrsim 0.9$, wind mass-loss rates $\dot{M} \gtrsim 10^{-4} M_\odot\,\text{yr}^{-1}$, and separations between about 50 and 1500 AU. The required $\dot{M}$ increases for binaries with smaller $q$ or primaries with faster wind. We estimate that $\ll 1\%$ of all collapsing massive stars satisfy the conditions on binary mass ratio and separation. Recombination luminosities due to colliding wind shells are at most a factor of 10 higher than for an otherwise unperturbed constant-velocity wind, but higher densities associated with wind acceleration close to the star provide much stronger signal., Comment: 15 pages, 13 figures, 1 table. Accepted to MNRAS

Sirius-like systems are wide binaries composed of a white dwarf (WD) and a companion of a spectral type earlier than M0. The WD progenitor evolves in isolation, but its wind during the AGB phase pollutes the companion surface and transfers some angular momentum. Within SHINE survey that uses SPHERE at the VLT, we acquired images of HD2133, HD114174, and CD-567708 and combined this data with high resolution spectra of the primaries, TESS, and literature data. We performed accurate abundance analyses for the MS. We found brighter J and K magnitudes for HD114174B than obtained previously and extended the photometry down to 0.95 micron. Our new data indicate a higher temperature and then shorter cooling age (5.57+/-0.02 Gyr) and larger mass (0.75+/-0.03 Mo) for this WD than previously assumed. This solved the discrepancy previously found with the age of the MS star. The two other WDs are less massive, indicating progenitors of ~1.3 Mo and 1.5-1.8 Mo for HD2133B and CD-56 7708B, respectively. We were able to derive constraints on the orbit for HD114174 and CD-56 7708. The composition of the MS stars agrees fairly well with expectations from pollution by the AGB progenitors of the WDs: HD2133A has a small enrichment of n-capture elements, which is as expected for pollution by an AGB star with a mass 3.0 Mo. On the other hand, none of these stars show the excesses of C that are expected to go along with those of n-capture elements. This might be related to the fact that these stars are at the edges of the mass range where we expect nucleosynthesis related to thermal pulses., Accepted on Astronomy and Astrophysics. 19 pages, 15 figures

Context. Relativistic jets are ubiquitous in astrophysics. High Mass Microquasars (HMMQs) are useful labs to study such jets, as they are relatively close and evolve over observable time scales. The ambient medium into which the jet propagates is, however, far from homogeneous. Corresponding simulation studies to date consider various forms of a wind-shaped ambient medium but typically neglect radiative cooling and relativistic effects. Aims. We investigate the dynamical and structural effects of radiative losses and system parameters on relativistic jets in HMMQs, from the jet launch to its propagation over several tens of orbital separations. Methods. We use 3D relativistic hydrodynamical simulations including parameterized radiative cooling derived from relativistic thermal plasma distribution to carry out parameter studies around two fiducial cases inspired by Cygnus X-1 and Cygnus X-3. Results. Radiative losses are found to be more relevant in Cygnus X-3 than Cygnus X-1. Varying jet power, jet temperature, or the wind of the donor star tends to have a larger impact at early times, when the jet forms and instabilities initially develop, than at later times when the jet has reached a turbulent state. Conclusions. Radiative losses may be dynamically and structurally relevant at least in Cygnus X-3 case, and thus should be examined in more detail., Comment: 28 pages, 25 figures, 8 tables. Accepted for publication in A&A

Even though tens of directly imaged companions have been discovered in the past decades, the number of directly confirmed multiplanet systems is still small. Dynamical analysis of these systems imposes important constraints on formation mechanisms of these wide-orbit companions. As part of the Young Suns Exoplanet Survey (YSES) we report the detection of a second planetary-mass companion around the 17 Myr-old, solar-type star TYC 8998-760-1 that is located in the Lower Centaurus Crux subgroup of the Scorpius-Centaurus association. The companion has a projected physical separation of 320 au and several individual photometric measurements from 1.1 to 3.8 microns constrain a companion mass of $6\pm1\,M_\mathrm{Jup}$, which is equivalent to a mass ratio of $q=0.57\pm0.10\%$ with respect to the primary. With the previously detected $14\pm3\,M_\mathrm{Jup}$ companion that is orbiting the primary at 160 au, TYC 8998-760-1 is the first directly imaged multiplanet system that is detected around a young, solar analog. We show that circular orbits are stable, but that mildly eccentric orbits for either/both components ($e > 0.1$) are chaotic on Gyr timescales, implying in-situ formation or a very specific ejection by an unseen third companion. Due to the wide separations of the companions TYC 8998-760-1 is an excellent system for spectroscopic and photometric follow-up with space-based observatories such as the James Webb Space Telescope., Accepted for publication in ApJL (12 pages, 8 figures, 2 tables)

Massive stars play an important role in the mechanical and chemical evolution of galaxies. Understanding the internal processes of these stars is vital to our understanding of their evolution and eventual end products. Deformations from spherical geometry are common for massive stars; however, the tools that are currently available for the study of these systems are almost exclusively one-dimensional. We present a new spectroscopic analysis tool tailored for massive stars that deviate from spherical symmetry. This code (entitled SPAMMS) is a spectroscopic patch model that takes the three-dimensional surface geometry of the system into account to produce spectral profiles at given phases and orientations. In using the Wilson-Devinney-like code PHOEBE in combination with the non local thermodynamic equilibrium (NLTE) radiative transfer code FASTWIND, we created a three-dimensional mesh that represents the surface geometry of our system and we assigned FASTWIND emergent intensity line profiles to each mesh triangle, which take the local parameters such as temperature, surface gravity, and radius into account. These line profiles were then integrated across the visible surface, where their flux contribution and radial velocity are taken into account, thus returning a final line profile for the visible surface of the system at a given phase. We demonstrate that SPAMMS can accurately reproduce the morphology of observed spectral line profiles for overcontact systems. Additionally, we show how line profiles of rapidly-rotating single stars differ when taking rotational distortion into account, and the effects that these can have on the determined parameters. Finally, we demonstrate the code's ability to reproduce the Rossiter-Mclaughlin and Struve-Sahade effects., 10 pages, 10 figures, accepted for publication in A&A

One of the proposed channels of binary black hole mergers involves dynamical interactions of three black holes. In such scenarios, it is possible that all three black holes merge in a so-called hierarchical merger chain, where two of the black holes merge first and then their remnant subsequently merges with the remaining single black hole. Depending on the dynamical environment, it is possible that both mergers will appear within the observable time window. Here we perform a search for such merger pairs in the public available LIGO and Virgo data from the O1/O2 runs. Using a frequentist p-value assignment statistics we do not find any significant merger pair candidates, the most significant being GW170809-GW151012 pair. Assuming no observed candidates in O3/O4, we derive upper limits on merger pairs to be $\sim11-110\ {\rm year^{-1}Gpc^{-3}}$, corresponding to a rate that relative to the total merger rate is $\sim 0.1-1.0$. From this we argue that both a detection and a non-detection within the next few years can be used to put useful constraints on some dynamical progenitor models., 7 pages, 4 figures

We investigate the role of binaries and chromospheric emission on HST photometry of globular clusters' stars. We quantify their respective effects on the position of stars in the chromosome map, especially among the first population. We computed atmosphere models and synthetic spectra for stars of different chemical compositions, based on isochrones produced by stellar evolution calculations with abundance variations representative of first and second populations in GCs. From this we built synthetic chromosome maps for a mixture of stars of different chemical compositions. We subsequently replaced a fraction of stars with binaries, or stars with chromospheric emission, using synthetic spectroscopy. We studied how the position of stars is affected in the chromosome map. Binaries can, in principle, explain the extension of the first population in the chromosome map. However, we find that given the binary fraction reported for GCs, the density of stars in the extended part is too small. Another difficulty of the binary explanation is that the shape of the distribution of the first population in the chromosome map is different in clusters with similar binary fractions. Additionally, we find that the contribution of chromospheric emission lines to the HST photometry is too small to have an observable impact on the shape of the chromosome map. Continuum chromospheric emission has an effect qualitatively similar to binaries. We conclude that binaries do have an impact on the morphology of the chromosome map of GCs, but they are unlikely to explain entirely the shape of the extended distribution of the first population stars. Uncertainties in the properties of continuum chromospheric emission of stars in GCs prevent any quantitative conclusion. Therefore, the origin of the extended first population remains unexplained., 11 pages, 10 figures. Accepted in Astronomy & Astrophysics

We present the novel algorithmically regularised integration method MSTAR for high accuracy ($|\Delta E/E| \gtrsim 10^{-14}$) integrations of N-body systems using minimum spanning tree coordinates. The two-fold parallelisation of the $\mathcal{O}(N_\mathrm{part}^2)$ force loops and the substep divisions of the extrapolation method allows for a parallel scaling up to $N_\mathrm{CPU} = 0.2 \times N_\mathrm{part}$. The efficient parallel scaling of MSTAR makes the accurate integration of much larger particle numbers possible compared to the traditional algorithmic regularisation chain (AR-CHAIN) methods, e.g. $N_\mathrm{part} = 5000$ particles on $400$ CPUs for $1$ Gyr in a few weeks of wall-clock time. We present applications of MSTAR on few particle systems, studying the Kozai mechanism and N-body systems like star clusters with up to $N_\mathrm{part} =10^4$ particles. Combined with a tree or a fast multipole based integrator the high performance of MSTAR removes a major computational bottleneck in simulations with regularised subsystems. It will enable the next generation galactic-scale simulations with up to $10^9$ stellar particles (e.g. $m_\star = 100 M_\odot$ for a $M_\star = 10^{11} M_\odot$ galaxy) including accurate collisional dynamics in the vicinity of nuclear supermassive black holes., Comment: 20 pages, 16 figures, accepted for publication in MNRAS

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I = 12 mag, and it was covered in great detail with almost 25 000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57 +/- 0.05 M-circle dot and 0.36 +/- 0.03 M-circle dot at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes., Polish National Sci-ence Centre (NCN): HARMONIA NCN grant [2015/18/M/ST9/00544]; Polish National Sci-ence Centre (NCN): OPUS NCN grant [2015/17/B/ST9/03167]; Polish National Sci-ence Centre (NCN): DAINA NCN grant [2017/27/L/ST9/03221]; European CommissionEuropean CommissionEuropean Commission Joint Research Centre [312430, 730890]; Polish Ministry of Higher EducationMinistry of Science and Higher Education, Poland [3040/7.PR/2014/2, DIR/WK/2018/12]; MAESTRO NCN grant [2014/14/A/ST9/00121]; DDT programmes [SW2016b12, A34DDT3]; UK Science and Technology Facilities CouncilUK Research & Innovation (UKRI)Science & Technology Facilities Council (STFC); Foundation for Polish Science (Program START); National Science Center, PolandNational Science Centre, Poland [ETIUDA 2018/28/T/ST9/00096]; Slovenian Research AgencySlovenian Research Agency - Slovenia [P1-0031, J1-8136]; COST Action GWverseEuropean Cooperation in Science and Technology (COST) [CA16104]; National Research Foundation of KoreaNational Research Foundation of Korea [2017R1A4A1015178]; STFCUK Research & Innovation (UKRI)Science & Technology Facilities Council (STFC); European Space Agency under the NELIOTA program [4000112943]; GROWTH project - National Science Foundation [1545949]; MINECO (Spanish Ministry of Economy) [ESP2016-80079-C2-1-R, ESP2014-55996-C2-1-R, RTI2018-095076B-C21]; ICCUB (Unidad de Excelencia a Maria de Maeztu) [MDM-2014-0369]; Agenzia Spaziale Italiana (ASI)Agenzia Spaziale Italiana (ASI) [ASI I/058/10/0, ASI 2014-025-R.1.2015]; Royal Society Research GrantRoyal Society of London [RG170230]; Ministry of Science and Technology (Taiwan)Ministry of Science and Technology, Taiwan [104-2112-M-008-012MY3, 104-2923-M-008-004-MY5]; European Research Council under the European Union's Seventh Framework Programme (FP/2007-2013)/ERC Grant [320964]; Gordon and Betty Moore FoundationGordon and Betty Moore Foundation [GBMF5490]; NSFNational Science Foundation (NSF) [AST-1515927, AST-0908816, AST-1412587]; Mt. Cuba Astronomical Foundation; Center for Cosmology and AstroParticle Physics at the Ohio State UniversityOhio State University; Chinese Academy of Sciences South America Center for Astronomy (CAS-SACA); Villum Foundation; MINECO under the Ramon y Cajal programme [RYC-2016-20254]; MINECO [AYA2017-86274-P]; AGAUR grantAgencia de Gestio D'Ajuts Universitaris de Recerca Agaur (AGAUR) [SGR-661/2017]; Science and Technology Facilities CouncilUK Research & Innovation (UKRI)Science & Technology Facilities Council (STFC) [ST/P000541/1]; STFC consolidated grant [ST/M001296/1]; Research Council of LithuaniaResearch Council of Lithuania (LMTLT) [S-LL-19-2]; TUB.ITAKTurkiye Bilimsel ve Teknolojik Arastirma Kurumu (TUBITAK); IKI; KFU; AST; National Science Fund (Bulgaria)National Science Fund of Bulgaria [DN 18/13-12.12.2017]; Institute of Astronomy via research project "Study of ICRF radio-sources and fast variable astronomical objects"; NAO Rozhen, BAS, via research project "Study of ICRF radio-sources and fast variable astronomical objects"; Ministry of Education, Science and Technological Development of the Republic of Serbia [176011, 176004, 176021]; DFGGerman Research Foundation (DFG)European Commission [SPP 1992, WA 1074/11-1]; NCN [2016/21/B/ST9/01126]; MINECO Ramon y Cajal programme [AYA2017-86274-P, RYJ-2016-20254]; Robert Martin Ayers Sciences Fund; National Aeronautics and Space AdministrationNational Aeronautics & Space Administration (NASA) [NNX08AR22G]; National Science FoundationNational Science Foundation (NSF) [AST-1238877]; W. M. Keck FoundationW.M. Keck Foundation; ICCUB (Unidad de Excelencia 'Maria de Maeztu') [MDM-2014-0369]; [3.6714.2017/8.9], This work relies on the results from the European Space Agency (ESA) space mission Gaia. Gaia data are being processed by the Gaia Data Processing and Analysis Consortium (DPAC). Funding for the DPAC is provided by national institutions, in particular the institutions participating in the Gaia Multi-Lateral Agreement (MLA). The Gaia mission website is https://www.cosmos.esa.int/gaia.In particular we acknowledge Gaia Photomet-ric Science Alerts Team, website http://gsaweb.ast.cam.ac.uk/alerts.We thank the members of the OGLE team for discussions and support. We also would like to thank the Polish Children Fund (KFnRD) for support of an internship of their pupils in Ostrowik Observatory of the Warsaw Univ-ersity, during which some of the data were collected, in particular we thank: Robert Nowicki, Michal Por.ebski and Karol Niczyj. The work presented here has been supported by the following grants from the Polish National Sci-ence Centre (NCN): HARMONIA NCN grant 2015/18/M/ST9/00544, OPUS NCN grant 2015/17/B/ST9/03167, DAINA NCN grant 2017/27/L/ST9/03221, as well as European Commission's FP7 and H2020 OPTICON grants (312430 and 730890), Polish Ministry of Higher Education support for OPTICON FP7, 3040/7.PR/2014/2, MNiSW grant DIR/WK/2018/12. PMr and JS acknowledge support from MAESTRO NCN grant 2014/14/A/ST9/00121 to Andrzej Udal-ski. We would like to thank the following members of the AAVSO for their amazing work with collecting vast amounts of data: Teofilo Arranz, James Boardman, Stephen Brincat, Geoff Chaplin, Emery Erdelyi, Rafael Farfan, William Goff, Franklin Guenther, Kevin Hills, Jens Jacobsen, Raymond Kneip, David Lane, Fernando Limon Martinez, Gianpiero Locatelli, Andrea Mantero, Attila Madai, Peter Meadows, Otmar Nickel, Arto Oksanen, Luis Perez, Roger Pieri, Ulisse Quadri, Diego Rodriguez Perez, Frank Schorr, George Sjoberg, Andras Timar, Ray Tomlin, Tonny Vanmunster, Klaus Wenzel, Thomas Wikander. We also thank the amateur observers from around the world, in particular, Pietro Capuozzo, Leone Trascianelli, Igor Zharkov from Ardingly College and Angelo Tomassini, Karl-Ludwig Bath. We also thank Roger Pickard from the British Astronomical Association and Matthias Penselin from the German Haus der Astronomie association for their contributions. KS thanks Dr. Dmitry Chulkov and Dr. Panagiotis Gavras for the interesting discussion of stellar multiplicity. We acknowledge support of DDT programmes SW2016b12 (WHT) and A34DDT3 (TNG). The INT, TNG and WHT are operated on the island of La Palma by the Isaac Newton Group of Telescopes in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias. The Liverpool Telescope is operated on the island of La Palma by Liverpool John Moores University in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias with financial support from the UK Science and Technology Facilities Council. SJF would like to thank the UCL students who assisted with the collection and checking of UCLO data for the observing campaign: Martina Aghopian, Ashleigh Arendt, Artem Barinov, Luke Barrett, Jasper Berry-Gair, Arjun Bhogal, Charles Bowesman, William Boyd, Andrei Cuceu, Michael Davies, Max Freedman, Gabriel Fu, Abirami Govindaraju, Iandeep Hothi, Clara Matthews Torres, Darius Modirrousta-Galian, Petru Neague, George Pattinson, Xiaoxi Song, and Brian Yu. P.Mr. acknowledges support from the Foundation for Polish Science (Program START) and the National Science Center, Poland (grant ETIUDA 2018/28/T/ST9/00096).; r r AC, AG and NI acknowledge the financial support from the Slovenian Research Agency (research core funding No. P1-0031 and project grant No. J1-8136)and networking support by the COST Action GWverse CA16104. Skinakas Observatory is a collaborative project of the University of Crete and the Foundation for Research and Technology-Hellas. Work by C.H. was supported by the grant (2017R1A4A1015178) of National Research Foundation of Korea. KW acknowledges funding from STFC, and thanks the University of Lei-cester for the investment in instrumentation. We gratefully acknowledge financial support by the European Space Agency under the NELIOTA program, contract No. 4000112943. This work has made use of data obtained with the Kryoneri Prime Focus Instrument, developed by the European Space Agency NELIOTA project on the 1.2 m Kryoneri telescope, which is operated by IAASARS, National Observatory of Athens, Greece. The Aristarchos tele-scope is operated on Helmos Observatory by the IAASARS of the National Observatory of Athens. This work was supported by the GROWTH project funded by the National Science Foundation under Grant No 1545949. This work was supported by the MINECO (Spanish Ministry of Economy) through grant ESP2016-80079-C2-1-R (MINECO/FEDER, UE) and ESP2014-55996-C2-1-R (MINECO/FEDER, UE) and MDM-2014-0369 of ICCUB (Unidad de Excelencia a Maria de Maeztu). This work was supported by the MINECO (Spanish Ministry of Economy) through grant ESP2016-80079-C2-1-R and RTI2018-095076B-C21 (MINECO/FEDER, UE), and MDM-2014-0369 of ICCUB (Unidad de Excelencia 'Maria de Maeztu'). The Joan Oro Telescope (TJO) of the Montsec Astronomical Observatory (OAdM) is owned by the Catalan Government and is operated by the Institute for Space Studies of Catalonia (IEEC). Support to this study has been provided by Agenzia Spaziale Italiana (ASI) through grants ASI I/058/10/0 and ASI 2014-025-R.1.2015. KW thanks Dipali Thanki and Ray McErlean for their technical support of the scientific programme of the University of Leicester observatory. This work was supported by Royal Society Research Grant RG170230. CCN thanks the funding from Ministry of Science and Technology (Taiwan) under the contracts 104-2112-M-008-012MY3 and 104-2923-M-008-004-MY5. The research leading to these results has received funding from the European Research Council under the European Union's Seventh Framework Programme (FP/2007-2013)/ERC Grant Agreement n. 320964 (WDTracer). We thank the Las Cumbres Observatory and its sta ff for its continuing support of the project. ASAS-SN is supported by the Gordon and Betty Moore Foundation through grant GBMF5490 to the Ohio State University and NSF grant AST-1515927. Development of ASAS-SN has been supported by NSF grant AST-0908816, the Mt. Cuba Astronomical Foundation, the Center for Cosmology and AstroParticle Physics at the Ohio State University, the Chinese Academy of Sciences South America Center for Astronomy (CAS-SACA), the Villum Foundation, and George Skestos. ARM acknowledges support from the MINECO under the Ramon y Cajal programme (RYC-2016-20254) and the AYA2017-86274-P grant, and the AGAUR grant SGR-661/2017. We acknowledge support from the Science and Technology Facilities Council (TB and RWW; ST/P000541/1). K.Horne acknowledges support from STFC consolidated grant ST/M001296/1. This work was partly supported by the Research Council of Lithuania, grant No. S-LL-19-2 Authors thank to TUB.ITAK, IKI, KFU, and AST for partial supports in using RTT150 (Russian-Turkish 1.5-m telescope in Antalya).; r r This work was partially funded by the subsidy 3.6714.2017/8.9 allocated to Kazan Federal University for the state assignment in the sphere of scientific activities. This research was partially supported by contract DN 18/13-12.12.2017 with the National Science Fund (Bulgaria). Work by YS was supported by an appointment to the NASA Postdoctoral Program at the Jet Propulsion Laboratory, California Institute of Technology, administered by Uni-versities Space Research Associ-ation through a contract with NASA. GD gratefully acknowledges the obser-ving grant support from the Institute of Astronomy and NAO Rozhen, BAS, via bilateral joint research project "Study of ICRF radio-sources and fast variable astronomical objects" (PI:G.Damljanovic). This work is a part of the Projects no. 176011 "Dynamics and kinematics of celestial bodies and systems", no. 176004 "Stellar physics" and no. 176021 "Visible and invisible matter in nearby galaxies: theory and observations" supported by the Ministry of Education, Science and Technological Development of the Republic of Serbia. YT acknowledges the support of DFG priority program SPP 1992 "Exploring the diversity of Extrasolar Planets" (WA 1074/11-1). This work of PMi, DM and ZK was supported by the NCN grant no. 2016/21/B/ST9/01126. ARM acknowledges support from the MINECO Ramon y Cajal programme RYJ-2016-20254 and grant AYA2017-86274-P and from the AGAUR grant SGR-661/2017. The work by C. R. was supported by an appointment to the NASA Postdoctoral Program at the Goddard Space Flight Center, administered by USRA through a contract with NASA. The Faulkes Telescope Project is an education partner of Las Cumbres Observatory (LCO). The Faulkes Telescopes are maintained and operated by LCO. This research was made possible through the use of the AAVSO Photomet-ric All-Sky Survey (APASS), funded by the Robert Martin Ayers Sciences Fund and NSF AST-1412587. The Pan-STARRS1 Surveys (PS1) and the PS1 public science archive have been made possible through contributions by the Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project O ffice, the Max-Planck Society and its participating institutes, the Max Planck Institute for Astronomy, Heidelberg and the Max Planck Institute for Extraterrestrial Physics, Garching, The Johns Hopkins University, Durham University, the University of Edinburgh, the Queen's University Belfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Observatory Global Telescope Network Incor-porated, the National Central University of Taiwan, the Space Telescope Sci-ence Institute, the National Aeronautics and Space Administration under Grant No. NNX08AR22G issued through the Planetary Science Division of the NASA Science Mission Directorate, the National Science Foundation Grant No. AST-1238877, the University of Maryland, Eotvos Lorand University (ELTE), the Los Alamos National Laboratory, and the Gordon and Betty Moore Foundation. Some of the data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Ins-titute of Technology, the University of California and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of the W. M. Keck Foundation.

The Laser Interferometer Space Antenna (LISA) will enable Galactic gravitational wave (GW) astronomy by individually resolving $ > 10^4$ signals from double white dwarf (DWD) binaries throughout the Milky Way. In this work we assess for the first time the potential of LISA data to map the Galactic stellar bar and spiral arms, since GWs are unaffected by stellar crowding and dust extinction unlike optical observations of the bulge region. To achieve this goal we combine a realistic population of Galactic DWDs with a high-resolution N-Body simulation a galaxy in good agreement with the Milky Way. We then model GW signals from our synthetic DWD population and reconstruct the structure of the simulated Galaxy from mock LISA observations. Our results show that while the low signal contrast between the background disc and the spiral arms hampers our ability to characterise the spiral structure, the stellar bar will instead clearly appear in the GW map of the bulge. The bar length and bar width derived from these synthetic observations are underestimated, respectively within $1\sigma$ and at a level greater than $2\sigma$, but the resulting axis ratio agrees to well within $1\sigma$, while the viewing angle is recovered to within one degree. These are competitive constraints compared to those from electromagnetic tracers, and they are obtained with a completely independent method. We therefore foresee that the synergistic use of GWs and electromagnetic tracers will be a powerful strategy to map the bar and the bulge of the Milky Way., Comment: Accepted by MNRAS, 14 pages, 11 figures, 1 table

We study the mutual evolution of the orbital properties of high mass ratio, circular, co-planar binaries and their surrounding discs, using 3D Smoothed Particle Hydrodynamics simulations. We investigate the evolution of binary and disc eccentricity, cavity structure and the formation of orbiting azimuthal over-dense features in the disc. Even with circular initial conditions, all discs with mass ratios $q>0.05$ develop eccentricity. We find that disc eccentricity grows abruptly after a relatively long time-scale ($\sim 400\textrm{--}700$ binary orbits), and is associated with a very small increase in the binary eccentricity. When disc eccentricity grows, the cavity semi-major axis reaches values $a_{\rm cav}\approx 3.5\, a_{\rm bin}$. We also find that the disc eccentricity correlates linearly with the cavity size. Viscosity and orbit crossing, appear to be responsible for halting the disc eccentricity growth -- eccentricity at the cavity edge in the range $e_{\rm cav}\sim 0.05\textrm{--} 0.35$. Our analysis shows that the current theoretical framework cannot fully explain the origin of these evolutionary features when the binary is almost circular ($e_{\rm bin}\lesssim 0.01$); we speculate about alternative explanations. As previously observed, we find that the disc develops an azimuthal over-dense feature in Keplerian motion at the edge of the cavity. A low contrast over-density still co-moves with the flow after 2000 binary orbits; such an over-density can in principle cause significant dust trapping, with important consequences for protoplanetary disc observations., 20 pages, 17 figures, 1 table, accepted for publication in MNRAS

Eclipsing binary stars allow derivation of accurate and precise masses and radii. When they reside in star clusters, properties of even higher precision, along with additional information, can be extracted. Asteroseismology of solar-like oscillations offers similar possibilities for single stars. We improve the previously established properties of the Hyades eclipsing binary HD27130 and re-assess the asteroseismic properties of the giant star $\epsilon$ Tau. The physical properties of these members of the Hyades are then used to constrain the helium content and age of the cluster. New multi-colour light curves were combined with multi-epoch radial velocities to yield masses and radii of HD27130. $T_{\rm eff}$ was derived from spectroscopy and photometry, and verified using the Gaia parallax. We estimate the cluster age from re-evaluated asteroseismic properties of $\epsilon$ Tau while using HD27130 to constrain the helium content. The masses and radii, and $T_{\rm eff}$ of HD 27130 were found to be $M=1.0245\pm0.0024 M_{\odot}$, $R=0.9226\pm0.015 R_{\odot}$, $T_{\rm eff}=5650\pm50$ K for the primary, and $M=0.7426\pm0.0016 M_{\odot}$, $R=0.7388\pm0.026 R_{\odot}$, $T_{\rm eff}=4300\pm100$ K for the secondary component. Our re-evaluation of $\epsilon$ Tau suggests that the previous literature estimates are trustworthy, and that the Hipparcos parallax is more reliable than the Gaia DR2 parallax. The helium content of HD27130 and thus of the Hyades is found to be $Y=0.27$ but with significant model dependence. Correlations with the adopted metallicity results in a robust helium enrichment law with $\frac{\Delta Y}{\Delta Z}$ close to 1.2. We estimate the age of the Hyades to be 0.9 $\pm$ 0.1 (stat) $\pm$ 0.1 (sys) Gyr in slight tension with recent age estimates based on the cluster white dwarfs. (abridged), Comment: 14 pages, 8 figures, accepted for publication in Astronomy & Astrophysics

The study of radiation emitted from black hole accretion discs represents a crucial way to understand the main physical properties of these sources, and in particular the black hole spin. Beside spectral analysis, polarimetry is becoming more and more important, motivated by the development of new techniques which will soon allow to perform measurements also in the X- and {\gamma}-rays. Photons emitted from black hole accretion discs in the soft state are indeed expected to be polarized, with an energy dependence which can provide an estimate of the black hole spin. Calculations performed so far, however, considered scattering as the only process to determine the polarization state of the emitted radiation, implicitly assuming that the temperatures involved are such that material in the disc is entirely ionized. In this work we generalize the problem by calculating the ionization structure of a surface layer of the disc with the public code CLOUDY , and then by determining the polarization properties of the emerging radiation using the Monte Carlo code STOKES . This allows us to account for absorption effects alongside scattering ones. We show that including absorption can deeply modify the polarization properties of the emerging radiation with respect to what is obtained in the pure-scattering limit. As a general rule, we find that the polarization degree is larger when absorption is more important, which occurs e.g. for low accretion rates and/or spins when the ionization of the matter in the innermost accretion disc regions is far from complete., Comment: 13 pages, 9 figures, accepted for publication in MNRAS

(Abridged) Eclipsing, spectroscopic double-lined binary star systems (SB2) are excellent laboratories for calibrating theories of stellar interior structure and evolution. We aim to investigate the mass discrepancy in binary stars. We study the effect of near-core mixing on the mass of the convective core of the stars and interpret the results in the context of the mass discrepancy. Two scenarios are considered, where individual stellar components of a binary system are treated independent of each other and where they are forced to have the same age and initial chemical composition. We find that the mass discrepancy is present in our sample and that it is anti-correlated with the surface gravity of the star. No correlations are found with other fundamental and atmospheric parameters, including the stellar mass. The mass discrepancy can be partially accounted for by increasing the amount of near-core mixing in stellar evolution models. We also find that ignoring the microturbulent velocity and turbulent pressure in stellar atmosphere models of hot evolved stars results in overestimation of their effective temperature by up to 8%. Together with enhanced near-core mixing, this can almost entirely account for the 30% mass discrepancy found for the evolved primary component of V380 Cyg. We find a strong link between the mass discrepancy and the convective core mass. The mass discrepancy can be solved by considering the combined effect of extra near-core boundary mixing and consistent treatment in the spectrum analysis of hot evolved stars. Our binary modelling results in convective core masses between 17 and 35% of the stellar mass, in excellent agreement with results from gravity-mode asteroseismology of single stars. This implies larger helium core masses near the end of the main sequence than anticipated so far., 22 pages (17 pages of main text + 5 pages of appendix), 14 figures. Accepted for publication in A&A

Observations of the star EPIC 220208795 (2MASS J01105556+0018507) reveal a single, deep and asymmetric eclipse, which we hypothesize is due to an eclipsing companion surrounded by a tilted and inclined opaque disk, similar to those seen around V928 Tau and EPIC 204376071. We aim to derive physical parameters of the disk and orbital parameters for the companion around the primary star. The modeling is carried out using a modified version of the python package pyPplusS, and optimization is done using emcee. The period analysis makes use of photometry from ground-based surveys, where we perform a period folding search for other possible eclipses by the disk. Parameters obtained by the best model fits are used to obtain the parameter space of the orbital parameters, while the most likely period obtained is used to constrain these parameters. The best model has an opaque disk with a radius of $1.14\pm0.03$ $R_{\odot}$, an impact parameter of $0.61\pm0.02$ $R_{\odot}$, an inclination of $77.01^{\circ}\pm0.03^{\circ}$, a tilt of $36.81^{\circ}\pm0.05^{\circ}$ and a transverse velocity of $77.45\pm0.05$ km s$^{-1}$. The two most likely periods are $\sim 290$ days and $\sim 236$ days, corresponding to an eccentricity of $\sim 0.7$, allowing us to make predictions for the epochs of the next eclipses. All models with tilted and inclined disks result in a minimum derived eccentricity of 0.3, which in combination with the two other known small transiting disk candidates V928 Tau and EPIC 204376071, suggest that there may be a common origin for their eccentric orbits., 13 pages, 10 figures, 10 tables, accepted for publication in Astronomy & Astrophysics. Reduced data and reduction script on GitHub at https://github.com/lizvdkamp/EPIC2202_disk

Gravitational-wave detections are now starting to probe the mass distribution of stellar-mass black holes (BHs). Robust predictions from stellar models are needed to interpret these. Theory predicts the existence of a gap in the BH mass distribution because of pair-instability supernova. The maximum BH mass below the gap is the result of pulsational mass loss. We evolve massive helium stars through their late hydrodynamical phases of evolution using the open-source MESA stellar evolution code. We find that the location of the lower edge of the mass gap at 45$M_\odot$ is remarkably robust against variations in the metallicity ($\approx 3M_\odot$), the treatment of internal mixing ($\approx 1M_\odot$), stellar wind mass loss ($\approx 4M_\odot$), making it the most robust predictions for the final stages of massive star evolution. The reason is that the onset of the instability is dictated by the near-final core mass, which in turn sets the resulting BH mass. However, varying $^{12}C\left(\alpha,\gamma\right)^{16}O$ reaction rate within its $1\sigma$ uncertainties shifts the location of the gap between $40M_\odot$ and $56M_\odot$. We provide updated analytic fits for population synthesis simulations. Our results imply that the detection of merging BHs can provide constraints on nuclear astrophysics. Furthermore, the robustness against metallicity suggests that there is a universal maximum for the location of the lower edge of the gap, which is insensitive to the formation environment and redshift for first-generation BHs. This is promising for the possibility to use the location of the gap as a "standard siren" across the Universe., Comment: 17 pages, 5 figures, 1 Table, Accepted ApJ

Millisecond pulsars are old, fast spinning neutron stars thought to have evolved from classical pulsars in binary systems, where the rapid rotation is caused by the accretion of matter and angular momentum from their companion. During this transition between classical and millisecond pulsars, there is a magnetic field reduction of around 4 orders of magnitude, which is not well understood. According to the standard scenario, the magnetic field is reduced as a consequence of accretion, either through ohmic dissipation or through screening by the accreted matter. We explored an alternative hypothesis in which the magnetic field is reduced through ambipolar diffusion before the accretion. This is particularly effective during the long epoch in which the pulsar has cooled, but has not yet started accreting. This makes the final magnetic field dependent on the evolution time of the companion star and thus its initial mass. We use observed binary systems to constrain the time available for the magnetic field decay based on the current pulsar companion: a helium white dwarf, a carbon-oxygen white dwarf, or another neutron star. Based on a simplified model without baryon pairing, we show that the proposed process agrees with the general distribution of observed magnetic field strengths in binaries, but is not able to explain some mildly recycled pulsars where no significant decay appears to have occurred. We discuss the possibility of other formation channels for these systems and the conditions under which the magnetic field evolution would be set by the neutron star crust rather than the core., Comment: 10 pages, 4 figures, 2 tables

In this paper, we address the question of how a central intermediate-mass black hole (IMBH) in a globular cluster (GC) affects dynamics, core collapse, and formation of the binary population. It is shown that the central IMBH forms a binary system that affects dynamics of stars in the cluster significantly. The presence of an intermediate-mass black hole with mass ≥ 1.0-1.7%of the total stellar mass in the cluster inhibits the formation of binary stars population.

In latest years, the high sensitivity of the instruments on-board the XMM-Newton and Chandra satellites allowed us to explore the properties of the X-ray emission from hot subdwarf stars. The small but growing sample of X-ray detected hot subdwarfs includes binary systems, in which the X-ray emission is due to wind accretion onto a compact companion (white dwarf or neutron star), as well as isolated sdO stars, in which X-rays are probably due to shock instabilities in the wind. X-ray observations of these low-mass stars provide information which can be useful for our understanding of the weak winds of this type of stars and can lead to the discovery of particularly interesting binary systems. Here we report the most recent results we have recently obtained in this research area., Comment: 8 pages, 3 figures. To appear in the Proceedings of the 8th Meeting on Hot Subdwarf Stars and Related Objects, 9-15 July 2017, Cracow, Poland. Eds. A. Baran, A. E. Lynas-Gray, Open Astronomy, in press

We present the detection of an additional pulsation mode (7.35 mHz) of a subdwarf B star, PG 1613+426, and periodic Observed minus Calculated (O-C) variations for two existing pulsations. PG 1613+426 is near the hot end of the sdB instability strip. One pulsation mode (6.94 mHz) was detected so far by Bonanno et al. (2002) and another pulsation mode candidate (7.05 mHz) was proposed with a confidence level above 90% by Kuassivi and Ferlet (2005). To constrain sdB star evolutional scenarios, this star was monitored in 2010, 2011, 2015, and 2017 as a part of a project for finding companions to sdB stars using the pulsation timing method. The photometric analysis of those data shows an additional 7.35 mHz pulsation mode as well as the previously detected 6.93 mHz mode. However the 7.05 mHz mode was not detected. Nightly amplitude changes of 7.35 mHz mode were observed in the 2011 data, however the 2017 data did not show nightly amplitude shifts. O-C variations were detected in both 6.93 mHz and 7.35 mHz pulsations, indicating that PG 1613+426 may have a low mass companion star. However, more observations are needed to confirm it.

One of the important issues regarding the final evolution of stars is the impact of binarity. A rich zoo of peculiar, evolved objects are born from the interaction between the loosely bound envelope of a giant, and the gravitational pull of a companion. However, binary interactions are not understood from first principles, and the theoretical models are subject to many assumptions. It is currently agreed upon that hot subdwarf stars can only be formed through binary interaction, either through common envelope ejection or stable Roche-lobe overflow (RLOF) near the tip of the red giant branch (RGB). These systems are therefore an ideal testing ground for binary interaction models. With our long term study of wide hot subdwarf (sdB) binaries we aim to improve our current understanding of stable RLOF on the RGB by comparing the results of binary population synthesis studies with the observed population. In this article we describe the current model and possible improvements, and which observables can be used to test different parts of the interaction model., 5 pages, 2 figures, accepted for publication in open astronomy as part of the conference proceedings for the 8th meeting on hot subdwarfs and related objects in Krakow

A radial velocity follow-up of the long-period sdOB+G1V type spectroscopic binary SB 744 revealed strong lines of fluorine and lead in the optical spectrum of the sdOB star and subsolar metallicity in the G1V companion. With high-quality observations and Gaia astrometric data, we aim at measuring the chemical composition and Galactic kinematics of the system to put it in context with known populations of hot subdwarfs. Such binary systems have high potential, as they give insights into the late stages of binary evolution as well as into the mysterious formation of stripped core helium-burning stars. We have analyzed the optical spectra with homogeneous atmospheric models to derive surface parameters of the binary members from a direct wavelength space decomposition and independently measured the atmospheric properties of the cool companion. The two independent methods reached consistent results, which, amended with constraints from spectral energy distributions provided a subdwarf mass. The Gaia astrometry allowed us to derive the Galactic kinematics of the system. SB 744 turned out to be an old, Population II system, that has gone through dramatic events. The hot subdwarf star belongs to the heavy-metal subclass of sdOB stars and we report super-solar abundances of lead, based on Pb III/IV lines. The He abundance of the hot subdwarf is the lowest among the known heavy-metal sdOB stars. The presence of fluorine implies that SB 744 was once a hierarchical triple system and the inner binary has merged in the near past. As an alternative scenario, single-star evolution through late core helium flash and atmospheric mixing can also produce the observed fluorine abundances. The atmospheric metal over-abundances currently observed are perhaps the results of a combination of mixing processes during formation and radiative support., Comment: Accepted for publication in Astronomy and Astrophysics

Our aim is to derive a self-consistent age, distance and composition for the globular cluster $47\,$Tucanae ($47\,$Tuc; NGC104). First, we reevaluate the reddening towards the cluster resulting in a nominal $E(B-V)=0.03\pm0.01$ as the best estimate. The $T_{\rm eff}$ of the components of the eclipsing binary member V69 is found to be $5900\pm72$ K from both photometric and spectroscopic evidence. This yields a true distance modulus $(m-M)_0=13.21\pm0.06$(random)$ \pm0.03 $(systematic) to $47\,$Tuc when combined with existing measurements of V69 radii and luminosity ratio. We then present a new completely self-consistent isochrone fitting method to ground based and $\textit{HST}$ cluster colour-magnitude diagrams and the eclipsing binary member V69. The analysis suggests that the composition of V69, and by extension one of the populations of $47\,$Tuc, is given by [Fe/H]$\sim-0.70$, [O/Fe]$\sim+0.60$, and $Y\sim0.250$ on the solar abundance scale of Asplund, Grevesse & Sauval. However, this depends on the accuracy of the model $T_{\rm eff}$ scale which is 50-75 K cooler than our best estimate but within measurement uncertainties. Our best estimate of the age of $47\,$Tuc is 11.8 Gyr, with firm ($3 ��$) lower and upper limits of 10.4 and 13.4 Gyr, respectively, in satisfactory agreement with the age derived from the white dwarf cooling sequence if our determination of the distance modulus is adopted., 19 pages, 8 figures, accepted for publication in MNRAS

The presence of an unseen ‘Planet 9’ on the outskirts of the Solar system has been invoked to explain the unexpected clustering of the orbits of several Edgeworth–Kuiper Belt Objects. We use N-body simulations to investigate the probability that Planet 9 was a free-floating planet (FFLOP) that was captured by the Sun in its birth star formation environment. We find that only 1–6 per cent of FFLOPs are ensnared by stars, even with the most optimal initial conditions for capture in star-forming regions (one FFLOP per star, and highly correlated stellar velocities to facilitate capture). Depending on the initial conditions of the star-forming regions, only 5–10 of 10 000 planets are captured on to orbits that lie within the constraints for Planet 9. When we apply an additional environmental constraint for Solar system formation – namely the injection of short-lived radioisotopes into the Sun’s protoplanetary disc from supernovae – we find the probability for the capture of Planet 9 to be almost zero., Monthly Notices of the Royal Astronomical Society. Letters, 472 (1), ISSN:1745-3933, ISSN:1745-3925

The current $super-active$ state of the recurrent nova T CrB has been observed with unprecedented detail. Previously published observations provide strong evidence that this state is due to an enhancement of the flow of material through the accretion disk, which increased the optical depth of its most internal region, the boundary layer. $NuSTAR$ and $Swift$ observed T CrB in 2015 September, roughly halfway through the rise to optical maximum. In our analysis of these data, we have found that: $i$) the UV emission, as observed with $Swift$/UVOT in 2015, was already as bright as it became in 2017, after the optical peak; $ii$) the soft X-ray emission (E $\lesssim$ 0.6 keV) observed in 2017 after the optical peak, on the other hand, had not yet developed during the rising phase in 2015; $iii$) the hard X-ray emitting plasma (E $\gtrsim$ 2 keV) had the same temperature and about half the flux of that observed during quiescence in 2006. This phenomenology is akin to that observed during dwarf novae in outburst, but with the changes in the spectral energy distribution happening on a far longer time scale., Accepted in ApJ

International audience; The stochastic gravitational wave background produced by supernovas, magnetars, and merger of binaries constituted of a pair of compact objects is reviewed and updated. The merger of systems composed of two black holes dominates by far the background signal, whose amplitude in the range 10–100 Hz is above the sensitivity of the planned Einstein laser interferometer (ET). The background signal at 25 Hz estimated by the laser interferometric gravitational observatory (LIGO)‐VIRGO collaboration, based on the available merger detection data, is in good agreement with the present theoretical predictions.