Evolution Of Close Binary Systems Research Articles

Orbital evolution of close binary systems: comparing viscous and wind-driven circumbinary disc models

ABSTRACT Previous work has shown that interactions between a central binary system and a circumbinary disc (CBD) can lead to the binary orbit either shrinking or expanding, depending on the properties of the disc. In this work, we perform two-dimensional hydrodynamical simulations of CBDs surrounding equal mass binary systems that are on fixed circular orbits, using the athena++ code in Cartesian coordinates. Previous studies have focused on discs where viscosity drives angular momentum transport. The aim of this work is to examine how the evolution of a binary system changes when angular momentum is extracted from the disc by a magnetized wind. In this proof-of-concept study, we mimic the effects of a magnetic field by applying an external torque that results in a prescribed radial mass flux through the disc. For three different values of the radial mass flux, we compare how the binary system evolves when the disc is either viscous or wind driven. In all cases considered, our simulations predict that the binary orbit should shrink faster by a factor of a few when surrounded by a wind-driven CBD compared to a corresponding viscous CBD. In-spiral time-scales of ∼106–107 yr are obtained for circular binaries surrounded by CBDs with masses typical of protoplanetary discs, indicating that significant orbital shrinkage can occur through binary–disc interactions during Class I/II pre-main-sequence phases.

On the non-dissipative tidal evolution of the misalignment between spin and orbital angular momenta

ABSTRACT We extend our previous work on the evolution of close binary systems with misaligned orbital and spin angular momenta resulting from non-dissipative tidal interaction to include all physical effects contributing to apsidal motion. In addition to tidal distortion of the primary by the compact secondary, these include relativistic Einstein precession and the rotational distortion of the primary. The influence of the precession of the line of nodes is included. The dependence of the tidal torque on the apsidal angle $\hat{\varpi }$ couples the apsidal motion to the rate of evolution of the misalignment angle β which is found to oscillate. We provide analytical estimates for the oscillation amplitude Δβ over a wide range of parameter space confirmed by numerical integrations. This is found to be more significant near critical curves on which ${\mathrm{d}}{\hat{\varpi }} /{\mathrm{d}}t=0$ for a specified β. We find that to obtain 0.1 < Δβ < ∼1, the mass ratio q > ∼1, the initial eccentricity should be modest, $\cos \beta \lt 1/\sqrt{5},$ with cos β < 0 corresponding to retrograde rotation, initially, and the primary rotation rate should be sufficiently large. The extended discussion of apsidal motion and its coupled evolution to the misalignment angle given here has potential applications to close binaries with anomalous apsidal motion as well as transiting exoplanets such as warm Jupiters.

X-RAY ASTRONOMY AND CLOSE BINARY SYSTEMS

The discovery in 1962 of the compact source Sco X-1, the first X-ray source located outside the Solar System, marked a new (“golden”) era in the study of close binary systems (CBS). Accreting neutron stars and first black hole candidates in close binary systems were discovered. The ability to “weigh” neutron stars and black holes in close binary systems made it possible to distinguish accreting black holes from neutron stars. The theory of accretion onto relativistic objects in close binary systems has been developed, as well as the theory of the evolution of close binary systems with mass exchange up to the latest stages, including binary black holes and neutron stars. Close binary systems have become the cutting edge of astrophysics. Subsequent observations of gravitation waves and observations on the EHT intercontinental radio interferometer with an angular resolution of ~10–5 arcsec made it possible to finally prove the existence of black holes in the Universe. Thus, the modern triumph of black holes has largely been due to the development of the science of close binary systems. The paper is based on a talk presented at the astrophysical memorial seminar “Novelties in Understanding the Evolution of Binary Stars”, dedicated to the 90th anniversary of Professor M.A. Svechnikov.

The role of stellar expansion on the formation of gravitational wave sources

ABSTRACT Massive stars are the progenitors of black holes and neutron stars, the mergers of which can be detected with gravitational waves (GW). The expansion of massive stars is one of the key factors affecting their evolution in close binary systems, but it remains subject to large uncertainties in stellar astrophysics. For population studies and predictions of GW sources, the stellar expansion is often simulated with the analytic formulae from Hurley et al. (2000). These formulae need to be extrapolated and are often considered outdated. In this work, we present five different prescriptions developed from 1D stellar models to constrain the maximum expansion of massive stars. We adopt these prescriptions to investigate how stellar expansion affects mass transfer interactions and in turn the formation of GW sources. We show that limiting radial expansion with updated 1D stellar models, when compared to the use of Hurley et al. (2000) radial expansion formulae, does not significantly affect GW source properties (rates and masses). This is because most mass transfer events leading to GW sources are initialized in our models before the donor star reaches its maximum expansion. The only significant difference was found for the mass distribution of massive binary black hole mergers (Mtot > 50 M⊙) formed from stars that may evolve beyond the Humphreys–Davidson limit, whose radial expansion is the most uncertain. We conclude that understanding the expansion of massive stars and the origin of the Humphrey–Davidson limit is a key factor for the study of GW sources.

Close Binary Stars. VII: On the Evolution of Close Binary Systems

Close Binary Stars. VII: On the Evolution of Close Binary Systems

Influence of tidal dissipation on outcomes of binary–single encounters between stars and black holes in stellar clusters

ABSTRACT In the cores of dense stellar clusters, close gravitational encounters between binary and single stars can frequently occur. Using the tsunami code, we computed the outcome of a large number of binary–single interactions involving two black holes (BHs) and a star to check how the inclusion of orbital energy losses due to tidal dissipation can change the outcome of these chaotic interactions. Each interaction was first simulated without any dissipative processes and then we systematically added orbital energy losses due to gravitational wave emission [using post-Newtonian (PN) corrections] and dynamical tides and recomputed the interactions. We find that the inclusion of tides increases the number of BH–star mergers by up to 75 per cent; however, it does not affect the number of BH–BH mergers. These results highlight the importance of including orbital energy dissipation due to dynamical tides during few-body encounters and evolution of close binary systems within stellar cluster simulations. Consistent with previous studies, we find that the inclusion of PN terms increases the number of BH–BH mergers during binary–single encounters. However, BH–star mergers are largely unaffected by the inclusion of these terms.

Effect of tides on the orbital evolution of the redback system PSR 1723–2837

ABSTRACT The standard model of stellar evolution in close binary systems assumes that during mass transfer episodes, the system is in a synchronized and circularized state. Remarkably, the redback system PSR J1723–2837 has an orbital period derivative $\dot{P}_{\rm orb}$ too large to be explained by this model. Motivated by this fact, we investigate the action of tidal forces in between two consecutive mass transfer episodes for a system under irradiation feedback, which is a plausible progenitor for PSR J1723–2837. We base our analysis on Hut’s treatment of equilibrium tidal evolution, generalized by considering the donor as a two layers object that may not rotate as a rigid body. We also analyse three different relations for the viscosity with the tidal forcing frequency. We found that the large value measured for $\dot{P}_{\rm orb}$ can be reached by systems where the donor star rotates slower (by few per cent) than the orbit just after mass transfer episodes. Van Staden & Antoniadis have observed this object and reported a lack of synchronism, opposite to that required by the Hut’s theory to account for the observed $\dot{P}_{\rm orb}$. Motivated by this discrepancy, we analyse photometric data obtained by the spacecraft Kepler second mission K2, with the purpose of identifying the periods present in PSR J1723–2837. We notice several periods close to those of the orbit and the rotation. The obtained periods pattern reveals the presence of a more complex phenomenology, which would not be well described in the frame of the weak friction model of equilibrium tides.

The Birth Function for Black Holes and Neutron Stars in Close Binaries

Abstract The mass function for black holes and neutron stars at birth is explored for mass-losing helium stars. These should resemble, more closely than similar studies of single hydrogen-rich stars, the results of evolution in close binary systems. The effects of varying the mass-loss rate and metallicity are calculated using a simple semi-analytic approach to stellar evolution that is tuned to reproduce detailed numerical calculations. Though the total fraction of black holes made in stellar collapse events varies considerably with metallicity, mass-loss rate, and mass cutoff, from 5% to 30%, the shapes of their birth functions are very similar for all reasonable variations in these quantities. Median neutron star masses are in the range 1.32–1.37 regardless of metallicity. The median black hole mass for solar metallicity is typically 8–9 if only initial helium cores below 40 (ZAMS mass less than 80 ) are counted, and 9–13 , in most cases, if helium cores with initial masses up to 150 (ZAMS mass less than 300 ) contribute. As long as the mass-loss rate as a function of mass exhibits no strong nonlinearities, the black hole birth function from 15 to 35 has a slope that depends mostly on the initial mass function for main-sequence stars. These findings imply the possibility of constraining the initial mass function and the properties of mass loss in close binaries using ongoing measurements of gravitational-wave radiation. The expected rotation rates of the black holes are briefly discussed.

Binary evolution leading to the formation of the very massive neutron star in the J0740+6620 binary system

ABSTRACT We study the evolution of close binary systems in order to account for the existence of the recently observed binary system containing the most massive millisecond pulsar ever detected, PSR J0740+6620, and its ultra-cool helium white dwarf companion. In order to find a progenitor for this object we compute the evolution of several binary systems composed by a neutron star and a normal donor star employing our stellar code. We assume conservative mass transfer. We also explore the effects of irradiation feedback on the system. We find that irradiated models also provide adequate models for the millisecond pulsar and its companion, so both irradiated and non irradiated systems are good progenitors for PSR J0740+6620. Finally, we obtain a binary system that evolves and accounts for the observational data of the system composed by PSR J0740+6620 (i.e. orbital period, mass, effective temperature and inferred metallicity of the companion, and mass of the neutron star) in a time scale smaller than the age of the Universe. In order to reach an effective temperature as low as observed, the donor star should have an helium envelope as demanded by observations.

Identifying the formation mechanism of redback pulsars

ABSTRACT We analyse the evolution of close binary systems containing a neutron star that lead to the formation of redback pulsars. Recently, there has been some debate on the origin of such systems and the formation mechanism of redbacks may still be considered as an open problem. We show that the operation of a strong evaporation mechanism, starting from the moment when the donor star becomes fully convective (or alternatively since the formation of the neutron star by accretion-induced collapse), produces systems with donor masses and orbital periods in the range corresponding to redbacks with donors appreciably smaller than their Roche lobes, i.e. they have low filling factors (lower than 0.75). Models of redback pulsars can be constructed assuming the occurrence of irradiation feedback. They have been shown to undergo cyclic mass transfer during the epoch at which they attain donor masses and orbital periods corresponding to redbacks, and stay in quasi-Roche lobe overflow conditions with high filling factors. We show that, if irradiation feedback occurs and radio ejection inhibits further accretion on to the neutron star after the first mass transfer cycle, the redback systems feature high filling factors. We suggest that the filling factor should be considered as a useful tool for discriminating among those redback formation mechanisms. We compare theoretical results with available observations and conclude that observations tend to favour models with high filling factors.

What information can we derive from historical Far Eastern guest stars for modern research on novae and cataclysmic variables?

ABSTRACTRecently, there have been several studies on the evolution of binary systems using historical data that are treated as facts in the chain of argument. This paper discusses six case studies of modern dwarf novae with suggested historical counterparts from the historical point of view, as well as the derived consequences for the physics of close binary systems (the dwarf novae Z Cam and AT Cnc, the nebula in M22, and the possible Nova 101, Nova 483 and Nova 1437). I consider the historical Far Eastern reports and, after a careful re-reading of the text, map the given information on to the sky. In some cases, the positions given in modern lists of classical nova–guest star pairs turn out to be wrong, or they have to be considered highly approximate: the historical position should, in most cases, be transformed into areas of the celestial sphere and not into point coordinates. Based on the correct information, I consider the consequences concerning the evolution of close binary systems. The result is that none of the cases of cataclysmic variables suggested to have a historical counterpart can be (fully) supported. Because the identification of the historical record of observation with the cataclysmic variables known today turns out to be always uncertain, a potential historical observation alone cannot be relied on to draw conclusions on the evolution of binaries. Evolution scenarios should be derived from astrophysical observations and modelling only.

Correlation of the rate of Type Ia supernovae with the parent galaxy properties: Light and shadows

Context. The identification of the progenitors of Type Ia supernovae (SNIa) is extremely important in several astrophysical contexts, ranging from stellar evolution in close binary systems to evaluating cosmological parameters. Determining the distribution of the delay times (DTD) of SNIa progenitors can shed light on their nature. The DTD can be constrained by analysing the correlation between the SNIa rate and those properties of the parent galaxy which trace the average age of their stellar populations. Aims. We investigate the diagnostic capabilities of this correlation by examining its systematics with the various parameters at play: simple stellar population models, the adopted description for the star formation history (SFH) in galaxies, and the way in which the masses of the galaxies are evaluated. Methods. We computed models for the diagnostic correlations for a variety of input ingredients and for a few astrophysically motivated DTD laws appropriate for a wide range of possibilities for the SNIa progenitors. The models are compared to the results of three independent observational surveys. Results. The scaling of the SNIa rate with the properties of the parent galaxy is sensitive to all input ingredients mentioned above. This is a severe limitation on the possibility to discriminate alternative DTDs. In addition, current surveys show some discrepancies for the reddest and bluest galaxies, likely because of limited statistics and the inhomogeneity of the observations. For galaxies with intermediate colours the rates are in agreement, leading to a robust determination of the productivity of SNIa from stellar populations of ≃0.8 events per 1000 M⊙. Conclusions. Large stastistics of SNIa events along with accurate measurements of the SFH in the galaxies are required to derive firm constraints on the DTD. The LSST will achieve these results by providing a homogeneous, unbiased, and vast database on both SNIa and galaxies.

Progenitors of gravitational wave mergers: binary evolution with the stellar grid-based code ComBinE

The first gravitational wave detections of mergers between black holes and neutron stars represent a remarkable new regime of high-energy transient astrophysics. The signals observed with LIGO-Virgo detectors come from mergers of extreme physical objects which are the end products of stellar evolution in close binary systems. To better understand their origin and merger rates, we have performed binary population syntheses at different metallicities using the new grid-based binary population synthesis code ComBinE. Starting from newborn pairs of stars, we follow their evolution including mass loss, mass transfer and accretion, common envelopes and supernova explosions. We apply the binding energies of common envelopes based on dense grids of detailed stellar structure models, make use of improved investigations of the subsequent Case BB Roche-lobe overflow and scale supernova kicks according to the stripping of the exploding stars. We demonstrate that all the double black hole mergers, GW150914, LVT151012, GW151226, GW170104, GW170608 and GW170814, as well as the double neutron star merger GW170817, are accounted for in our models in the appropriate metallicity regime. Our binary interaction parameters are calibrated to match the accurately determined properties of Galactic double neutron star systems, and we discuss their masses and types of supernova origin. Using our default values for the input physics parameters, we find a double neutron star merger rate of about 3.0 Myr^-1 for Milky-Way equivalent galaxies. Our upper limit to the merger-rate density of double neutron stars is R=400 yr^-1 Gpc^-3 in the local Universe (z=0).

Evolution of redback radio pulsars in globular clusters

We study the evolution of close binary systems composed of a normal, intermediate mass star and a neutron star considering a chemical composition typical of that present in globular clusters (Z = 0.001). We look for similarities and differences with respect to solar composition donor stars, which we have extensively studied in the past. As a definite example, we perform an application on one of the redbacks located in a globular cluster. We performed a detailed grid of models in order to find systems that represent the so-called redback binary radio pulsar systems with donor star masses between 0.6 and 2.0 solar masses and orbital periods in the range 0.2 - 0.9 days. We find that the evolution of these binary systems is rather similar to those corresponding to solar composition objects, allowing us to account for the occurrence of redbacks in globular clusters, as the main physical ingredient is the irradiation feedback. Redback systems are in the quasi-RLOF state, that is, almost filling their corresponding Roche lobe. During the irradiation cycle the system alternates between semi-detached and detached states. While detached the system appears as a binary millisecond pulsar, called a redback. Circumstellar material, as seen in redbacks, is left behind after the previous semi-detached phase. The evolution of binary radio pulsar systems considering irradiation successfully accounts for, and provides a way for, the occurrence of redback pulsars in low-metallicity environments such as globular clusters. This is the case despite possible effects of the low metal content of the donor star that could drive systems away from redback configuration.

The evolution of close binary stars

A review of our current understanding of the physics and evolution of close binary stars with various masses under the influence of the nuclear evolution of their components and their magnetic stellar winds is presented. The role of gravitational-wave radiation by close binaries on their evolution and the loss of their orbital angular momentum is also considered. The final stages in the evolution of close binary systems are described. The review also notes the main remaining tasks related to studies of the physics and evolution of various classes of close binaries, including analyses of collisions of close binaries and supermassive black holes in galactic nuclei. Such a collision could lead to the capture of one of the components by the black hole and the acceleration of the remaining component to relativistic speeds.

BINARY CENTRAL STARS OF PLANETARY NEBULAE DISCOVERED THROUGH PHOTOMETRIC VARIABILITY. III. THE CENTRAL STAR OF ABELL 65

A growing number of close binary stars are being discovered among central stars of planetary nebulae. Recent and ongoing surveys are finding new systems and contributing to our knowledge of the evolution of close binary systems. The push to find more systems was largely based on early discoveries which suggested that 10 to 15% of all central stars are close binaries. One goal of this series of papers is confirmation and classification of these systems as close binaries and determination of binary system parameters. Here we provide time-resolved multi-wavelength photometry of the central star of Abell 65 as well as further analysis of the nebula and discussion of possible binary--nebula connections. Our results for Abell 65 confirm recent work showing that it has a close, cool binary companion, though several of our model parameters disagree with the recently published values. With our longer time baseline of photometric observations from 1989--2009 we also provide a more precise orbital period of 1.0037577 days.

THE QUASI-ROCHE LOBE OVERFLOW STATE IN THE EVOLUTION OF CLOSE BINARY SYSTEMS CONTAINING A RADIO PULSAR

We study the evolution of close binary systems formed by a normal (solar composition), intermediate mass donor star together with a neutron star. We consider models including irradiation feedback and evaporation. These non-standard ingredients deeply modify the mass transfer stages of these binaries. While models that neglect irradiation feedback undergo continuous, long standing mass transfer episodes, models including these effect suffer a number cycles of mass transfer and detachment. During mass transfer the systems should reveal themselves as low-mass X-ray binaries (LMXBs), whereas when detached they behave as a binary radio pulsars. We show that at these stages irradiated models are in a Roche lobe overflow (RLOF) state or in a quasi-RLOF state. Quasi-RLOF stars have a radius slightly smaller than its Roche lobe. Remarkably, these conditions are attained for orbital period and donor mass values in the range corresponding to a family of binary radio pulsars known as "redbacks". Thus, redback companions should be quasi-RLOF stars. We show that the characteristics of the redback system PSR J1723-2837 are accounted for by these models. In each mass transfer cycle these systems should switch from LMXB to binary radio pulsar states with a timescale of \sim million years. However, there is recent and fast growing evidence of systems switching on far shorter, human timescales. This should be related to instabilities in the accretion disc surrounding the neutron star and/or radio ejection, still to be included in the model having the quasi-RLOF state as a general condition.

UNDERSTANDING THE EVOLUTION OF CLOSE BINARY SYSTEMS WITH RADIO PULSARS

We calculate the evolution of close binary systems (CBSs) formed by a neutron star (behaving as a radio pulsar) and a normal donor star, evolving either to helium white dwarf (HeWD) or ultra short orbital period systems. We consider X-ray irradiation feedback and evaporation due to radio pulsar irradiation. We show that irradiation feedback leads to cyclic mass transfer episodes, allowing CBSs to be observed in-between as binary radio pulsars under conditions in which standard, non-irradiated models predict the occurrence of a low mass X-ray binary. This behavior accounts for the existence of a family of eclipsing binary systems known as redbacks. We predict that redback companions should almost fill their Roche lobe, as observed in PSR J1723-2837. This state is also possible for systems evolving with larger orbital periods. Therefore, binary radio pulsars with companion star masses usually interpreted as larger than expected to produce HeWDs may also result in such {\it quasi - Roche Lobe Overflow} states, rather than hosting a carbon-oxygen WD. We found that CBSs with initial orbital periods $\mathrm{P_{i}<1}$ day evolve into redbacks. Some of them produce low mass HeWDs, and a subgroup with shorter $\mathrm{P_{i}}$ become black widows (BWs). Thus, BWs descent from redbacks, although not all redbacks evolve into BWs. There is mounting observational evidence favoring that BW pulsars are very massive ($\mathrm{\gtrsim 2\; M_{\odot}}$). As they should be redback descendants, redback pulsars should also be very massive, since most of the mass is transferred before this stage.

Wolf-Rayet stars with relativistic companions

The evolution of close binary systems containing Wolf-Rayet (WR) stars and black holes (BHs) is analyzed numerically. Both the stellar wind from the donor star itself and the induced stellar wind due to irradiation of the donor with hard radiation arising during accretion onto the relativistic component are considered. The mass and angular momentum losses due to the stellar wind are also taken into account at phases when the WR star fills its Roche lobe. It is shown that, if a WR star with a mass higher than ∼10M ⊙ fills its Roche lobe in an initial evolutionary phase, the donor star will eventually lose contact with the Roche lobe as the binary loses mass and angular momentum via the stellar wind, suggesting that the semi-detached binary will become detached. The star will remain a bright X-ray source, since the stellar wind that is captured by the black hole ensures a near-Eddington accretion rate. If the initial mass of the helium donor is below ∼5M ⊙, the donor may only temporarily detach from its Roche lobe. Induced stellar wind plays a significant role in the evolution of binaries containing helium donors with initial masses of ∼2M ⊙. We compute the evolution of three observed WR-BH binaries: Cyg X-3, IC 10 X-1, and NGC 300 X-1, as well as the evolution of the SS 433 binary system, which is a progenitor of such systems, under the assumption that this binary will avoid a common-envelope stage in its further evolution, as it does in its current evolutionary phase.

EFFECTS OF ROTATIONALLY INDUCED MIXING IN COMPACT BINARY SYSTEMS WITH LOW-MASS SECONDARIES AND IN SINGLE SOLAR-TYPE STARS

Many population synthesis and stellar evolution studies have addressed the evolution of close binary systems in which the primary is a compact remnant and the secondary is filling its Roche lobe, thus triggering mass transfer. Although tidal locking is expected in such systems, most studies have neglected the rotationally-induced mixing that may occur. Here we study the possible effects of mixing in the mass-losing stars for a range in secondary star masses and metallicities. We find that tidal locking can induce rotational mixing prior to contact and thus affect the evolution of the secondary star if the effects of the Spruit-Tayler dynamo are included both for angular momentum and chemical transport. Once contact is made, the effect of mass transfer tends to be more rapid than the evolutionary time scale, so the effects of mixing are no longer directly important, but the mass transfer strips matter to inner layers that may have been affected by the mixing. These effects are enhanced for secondaries of 1-1.2 Msun and for lower metallicities. We discuss the possible implications for the paucity of carbon in the secondaries of the cataclysmic variable SS Cyg and the black hole candidate XTE J1118+480 and for the progenitor evolution of Type Ia supernovae. We also address the issue of the origin of blue straggler stars in globular and open clusters. We find that for models that include rotation consistent with that observed for some blue straggler stars, evolution is chemically homogeneous. This leads to tracks in the HR diagram that are brighter and bluer than the non-rotating main-sequence turn-off point. Rotational mixing could thus be one of the factors that contribute to the formation of blue stragglers.