The study of the multiplicity of massive stars gives hints of their formation processes and their evolution path. Optical interferometry is mandatory to fulfill our knowledge of their multiplicity by probing the separation gap between 1 and 50 mas. We demonstrated the capability of the new interferometric instrument MIRC-X, located at the CHARA array, to study a large sample of more than 120 (H < 7.5) O-type stars. We observed 29 O-type star systems, including a couple of systems in average atmospheric conditions around a magnitude of H = 7.5. Out of these 29 systems, we detected 18 companions in 16 different systems, resulting in a multiplicity fraction fm = 16/29 = 0.55, and a companion fraction of fc = 18/29 = 0.62. We observed for the first time 11 of these detected companions. This study concludes that a large survey on more than 120 Northern O-type stars is possible with MIRC-X.

Recent abrupt changes of CW Leonis may indicate that we are witnessing the moment that the central carbon star is evolving off the Asymptotic Giant Branch (AGB) and entering into the pre-planetary nebula (PPN) phase. The recent appearance of a red compact peak at the predicted stellar position is possibly an unveiling event of the star, and the radial beams emerging from the stellar position resemble the feature of the PPN Egg Nebula. The increase of light curve over two decades is also extraordinary, and it is possibly related to the phase transition. Decadal-period variations are further found in the residuals of light curves, in the relative brightness of radial beams, and in the extended halo brightness distribution. Further monitoring of the recent dramatic and decadal-scale changes of this most well-known carbon star CW Leonis at the tip of AGB is still highly essential, and will help us gain a more concrete understanding on the conditions for transition between the late stellar evolutionary phases.

We run numerical simulations of massive colliding wind binaries, and quantify the accretion onto the secondary under different conditions. We set 3D simulation of a LBV–WR system and vary the LBV mass loss rate to obtain different values of wind momentum ratio η. We show that the mean accretion rate for stationary systems fits a power law Macc∝ η–1.6 for a wide range of η, until for extremely small η saturation in the accretion is reached. We find that the stronger the primary wind, the smaller the opening angle of the colliding wind structure (CWS), and compare it with previous analytical estimates. We demonstrate the efficiency of clumpy wind in penetrating the CWS and inducing smaller scale clumps that can be accreted. We propose that simulations of colliding winds can reveal more relations as the ones we found, and can be used to constrain stellar parameters.

The population of black widows, binary systems containing a millisecond pulsar and a very low-mass companion star exposed to the high-energy pulsar wind, has grown exponentially in the past few years. The number of black widow candidates is now over 30 systems, but only 14 have been confirmed so far. Their relevance in analysing the extremes of the neutron stars properties led to multiwavelength dedicated studies that revealed a rich phenomenology. In this work, we provide a glimpse into the black widow class through modelling of high-cadence multi-band light curves of 6 systems, accounting for almost half of the confirmed population. A better understanding of the black widow population, which hosts some of the most massive and fastest spinning neutron stars, will ultimately benefit future modelling of compact object mergers.

During the late stages of a neutron star binary inspiral finite-size effects come into play, with the tidal deformability of the supranuclear density matter leaving an imprint on the gravitational-wave signal. As demonstrated in the case of GW170817—the first direct detection of gravitational waves from a neutron star binary—this can lead to strong constraints on the neutron star equation of state. As detectors become more sensitive, effects which may have a smaller influence on the neutron star tidal deformability need to be taken into consideration. Dynamical effects, such as oscillation mode resonances triggered by the orbital motion, have been shown to contribute to the tidal deformability, especially close to the neutron star coalesence. We calculate the contribution of the various stellar oscillation modes to the tidal deformability and demonstrate the (anticipated) dominance of the fundamental mode, showing what the impact of the matter composition is on the tidal deformability.

Knowing the late stages of the stellar evolution is crucial for understanding the fate of planets around subdwarfs and white dwarfs. Simulations by (Staff et al.2016) show, that exoplanets engulfed in the extending stellar envelope will quickly spiral down onto the parent star. Therefore, we do not expect to find planets on close by orbits to the subdwarfs (Blokesz et al.2019) or white dwarfs. However, the recent observation of planetary debris around WD 1145+017 white dwarf suggests, there might exists planets farther away from these stars. Using binarograms, O-C diagrams and Fourier transform for the Kepler space telescope data, we investigate a problem of missing planets around white dwarfs in binary systems, single white dwarfs and subdwarfs type B. The last ones, being the only stars which (due to the lack of hydrogen) go directly to the white dwarf cooling track after their red giant phase.

Formation of close double white dwarfs likely requires the initial binary system to evolve through two successive common envelope (CE) phases. A prominent method for describing CE outcomes involves defining an ejection efficiency, αeff, which quantifies the fraction of orbital energy available to unbind the envelope. Reproducing observed post-CE orbital parameters has proven difficult for numerical simulations, as the companion’s decaying orbit fails to eject the envelope. The ejection failure seen in numerical simulations may be resolved with a proper treatment of convection, whereby the binary orbit shrinks before energy can drive ejection. Where the orbital decay timescale exceeds the convective transport timescale, the energy released during inspiral is carried to the stellar surface and radiated away. By including convection, we produce sub-day post-CE orbital periods, a result consistent with observations. We comment on the effects of convection for the population of double white dwarfs that evolve through two CEs.

Stars are embedded in different environments of Giant Molecular Clouds during their formation phase. Despite this fact, it is common practice to assume an isolated spherical core as the initial condition for models of individual star formation. To avoid the uncertainties of initial and boundary conditions, we use an alternative approach of zoom-in simulations to account for the environment in which protostars form. Our models show that injections of 26Al from a close-by supernova into the young solar system were highly unlikely. Moreover, we find that the accretion process of protostars is heterogeneous and environment-dependent.

What are the formation channels of merging black holes and neutron stars? The first two observing runs of Advanced LIGO and Virgo give us invaluable insights to address this question, but a new approach to theoretical models is required, in order to match the challenges posed by the new data. In this review, I discuss the impact of stellar winds, core-collapse and pair instability supernovae on the formation of compact remnants in both isolated and dynamically formed binaries. Finally, I show that dynamical processes, such as the runaway collision scenario and the Kozai-Lidov mechanism, leave a clear imprint on the demography of merging systems.

. A fraction of high-mass X-ray binaries are supergiant fast X-ray transients. These systems have on average low X-ray luminosities, but display short flares during which their X-ray luminosity rises by a few orders of magnitude. The leading model for the physics governing this X-ray behaviour suggests that the winds of the donor OB supergiants are magnetized. In agreement with this model, the first spectropolarimetric observations of the SFXT IGR J11215-5952 using the FORS 2 instrument at the Very Large Telescope indicate the presence of a kG longitudinal magnetic field. Based on these results, it seems possible that the key difference between supergiant fast X-ray transients and other high-mass X-ray binaries are the properties of the supergiant’s stellar wind and the physics of the wind’s interaction with the neutron star magnetosphere.

Previous circular polarization observations obtained with the ESO FOcal Reducer low dispersion spectrograpgh at the VLT in 2007–2008 revealed the presence of a weak longitudinal magnetic field on the surface of the optical component of the X-ray binary Cyg X-1, which contains a black hole and an O9.7Iab supergiant on a 5.6 d orbit. In this contribution we report on recently acquired FORS 2 spectropolarimetric observations of Cyg X-1 along with measurements of a few additional high-mass X-ray binaries.

The SALT transient follow-up programme began in 2016 and will continue for 5 semesters (until 31 Oct 2018), with an expectation of renewal thereafter. It is currently the only SALT Large Science Programme, and was awarded ~250 ksec. per semester, with a significant fraction (60%) given for the highest priority target-of-opportunity time. The aim is to characterise and study transients across a wide range of classes, currently including (from closest to most distant) cataclysmic variables, novæ and other associated eruptive variables, low- and high-mass X-ray binaries, OGLE and Gaia transients (including tidal disruption events), super-luminous and unusual core-collapse supernovæ, kilonovæ and other candidate optical counterparts to gravitational-wave events, flaring blazars and AGN, and gamma-ray bursts. This programme currently involves four SALT partners, of which South Africa is the major contributor of time (74%) and resources and includes five institutions with over 30 co-investigators. This talk reviewed the nature of the programme and highlighted some of the results to date.

It is widely accepted that binary interactions are responsible for the shaping of planetary nebula. However, these binary interactions and evolutionary channels are poorly understood. Our recent study revealed a newly discovered population of low-luminosity, low-metallicity, likely binaries in the Magellanic Clouds: dusty post-RGB stars. They are likely to have evolved off the RGB via binary interaction. In this paper we present preliminary results of the first radial velocity monitoring of the post-RGB stars in the Small Magellanic Cloud (SMC) and the implications on stellar (binary) evolution. We also investigate their link, if any, to the planetary nebulae systems.

As we enter the era of gravitational wave astronomy, we are beginning to collect observations which will enable us to explore aspects of astrophysics of massive stellar binaries which were previously beyond reach. In this paper we describe COMPAS (Compact Object Mergers: Population Astrophysics and Statistics), a new platform to allow us to deepen our understanding of isolated binary evolution and the formation of gravitational-wave sources. We describe the computational challenges associated with their exploration, and present preliminary results on overcoming them using Gaussian process regression as a simulation emulation technique.

Using the Yunnan-II evolutionary population synthesis models comprising binary stars, we find that the inclusion of binary stars can raise the derived stellar metallicity Z* and/or age t (degeneracy problem), raise the stellar mass M*, lower the gaseous metallicity Zgas and star formation rate (SFR) of galaxies. This means that a few stars form recently in galaxies, while more stars form during the entire evolution process when considering binary stars. If the degeneracy between t and Z* can be broken, its effect on the feedback process and star formation history can be determined.

Perturbation theory in the three body problem has greatly advanced our ability to understand and model a variety of systems ranging from artificial satellites to stars and from extrasolar planets to asteroid-Jupiter interactions. In a series of papers, we developed an analytical technique for estimating the orbital eccentricity of the inner binary in hierarchical triple systems. The method combined the secular theory with calculations of short period terms. The derivation of the short term component was based on an expansion of the rate of change of the Runge-Lenz vector by using first order perturbation theory, while canonical perturbation theory was used to investigate the secular evolution of the system. In the present work we extend the calculation to the orbit of the outer binary. At the same time, we provide an improved version for some previous results. A post-Newtonian correction is included in our model. Our analytical estimates are compared with numerical and analytical results on the subject and applications to stellar triples and extrasolar planets are discussed.

Barium stars are extrinsic Asymptotic Giant Branch (AGB) stars. They present the s-enhancement characteristic for AGB and post-AGB stars, but are in an earlier evolutionary stage (main sequence dwarfs, subgiants, red giants). They are believed to form in binary systems, where a more massive companion evolved faster, produced the s-elements during its AGB phase, polluted the present barium star through stellar winds and became a white dwarf. The samples of barium stars of Allen & Barbuy (2006) and of Smiljanic et al. (2007) are analysed here. Spectra of both samples were obtained at high-resolution and high S/N. We compare these observations with AGB nucleosynthesis models using different initial masses and a spread of 13C-pocket efficiencies. Once a consistent solution is found for the whole elemental distribution of abundances, a proper dilution factor is applied. This dilution is explained by the fact that the s-rich material transferred from the AGB to the nowadays observed stars is mixed with the envelope of the accretor. We also analyse the mass transfer process, and obtain the wind velocity for giants and subgiants with known orbital period. We find evidence that thermohaline mixing is acting inside main sequence dwarfs and we present a method for estimating its depth.

Bipolar molecular outflows have been observed and studied extensively in the past, but some recent observations of periodic variations in maser intensity pose new challenges. Even quasi-periodic maser flares have been observed and reported in the literature. Motivated by these data, we have tried to study situations in binary systems with specific attention to the two observed features, i.e., the bipolar flows and the variabilities in the maser intensity. We have studied the evolution of spherically symmetric wind from one of the bodies in the binary system, in the plane of the binary. Our approach includes the analytical study of rotating flows with numerical computation of streamlines of fluid particles using PLUTO code. We present the results of our findings assuming simple configurations, and discuss the implications.

To date, there are several reported exoplanet detections within binary star systems. These findings are based on radial velocity data for the target star. However, the companion star could in turn have a companion of which we are not aware. We describe how this hidden binary system affects the radial velocity of the target star, mimicking a planet in some circumstances We also explain what can be done in practice in order to distinguish between these two effects.

To inspire and provoke lively discussions, I argue that the accuracy of the basic physical properties of stars, based on analyses of well-observed detached binaries, might be worse than usually believed. I offer some ways how to deal with the situation. I end with a few comments on the studies of extra-solar planets.