The evolution of stars paired with supermassive black holes

A star located in the close vicinity of a supermassive black hole (SMBH) in a galactic nucleus or a globular-cluster core could form a close binary with the SMBH, with the star possibly filling its Roche lobe. The evolution of such binary systems is studied assuming that the SMBH mainly accretes matter from the companion star and that the presence of gas in the vicinity of the SMBH does not appreciably influence variations in the star’s orbit. The evolution of the star–SMBH system is mainly determined by the same processes as those determining the evolution of ordinary binaries. The main differences are that the star is subject to an incident flux of hard radiation arising during the accretion of matter by the SMBH, and, in detached systems, the SMBH captures virtually all the wind emitted by its stellar companion, which appreciably influences the evolution of the major axis of the orbit. Moreover, the exchange between the orbital angular momentum and the angular momentum of the overflowing matter may not be entirely standard in such systems. The computations assume that there will be no such exchange of angular momentum if the characteristic timescale for mass transfer is shorter than the thermal time scale of the star. The absorption of external radiation in the stellar envelope was computed using the same formalism applied when computing the opacity of the stellar matter. The numerical simulations show that, with the adopted assumptions, three types of evolution are possible for such a binary system, depending on the masses and the initial separation of the SMBH and star. Type I evolution leads to the complete destruction of the star. Only this type of evolution is realized for low-mass main-sequence (MS) stars, even those with large initial separations from their SMBHs. Massive MS stars will also be destroyed if the initial separation is sufficiently small. However, two other types of evolution are possible for massive stars, with a determining role in the time variations of the parameters of the star–SMBH system being played by the possible growth of the massive star into a red giant during the time it is located in the close vicinity of the SMBH. Type II evolution can be realized for massive MS stars that are initially farther from the SMBH than in the case of disruption. In this case, the massive star fills its Roche lobe during its expansion, but is not fully destroyed; the star retreats inside its Roche lobe after a period of intense mass loss. This type of evolution is characterized by an increase in the orbital period of the system with time. As a result, the remnant of the star (its former core) is preserved as a white dwarf, and can end up at a fairly large distance from the SMBH. Type III evolution can be realized formassiveMSstars that are initially located still farther from their SMBHs, and also for massive stars that are already evolved at the initial time. In these cases, the star moves away from the SMBH without filling its Roche lobe, due to its intense stellar wind. The remnants of such stars can also end up at a fairly large distances from their SMBHs.     

The spectrum of EG UMa

Spectra and radial-velocity curves for the precataclysmic variable EG UMa are analyzed. The system ephemeris has been improved and all system parameters determined. The parameters of the secondary are consistent with current evolutionary models for single main-sequence M stars with atmospheric metallicities exceeding the solar value by 0.5 dex. We have verified that the rotational velocity of the red dwarf exceeds the velocity corresponding to synchronous orbital motion by a factor of two to three. This suggests that the efficiency of tidal interactions between the components in the synchronization of their motion is low. The observed Ca II emission lines display reflection effects in a number of uniform spectra of EG UMa obtained during the quiescent state of the secondary.     

Changes in the parameters of the accretion disk around the dwarf nova SDSS J090350.73+330036.1 during an outburst

R-band photometric light curves of the eruptive eclipsing binary SDSS J090350.73+330036.1 obtained during a superoutburst in May 2010 (JD 2455341-2455347) are analyzed. Observations covering an interval near the outburst maximum and the post-maximum decrease by 0.7 ^m are presented. Oscillations (superhumps) whose period differs from the orbital period by several percent are observed in the light curve together with eclipses, suggesting that the studied system is a SU UMa dwarf nova. A ??spiral arm?? model is used to fit the light curves and determine the parameters of the accretion disk and other components of the binary system. Together with a hot line, this model takes into account, geometrical inhomogeneities on the surface of the accretion disk, namely, two thickenings at its outer edge that decrease exponentially in the vertical direction with approach toward the white dwarf. The increase in the R-band flux from the system during the superoutburst mainly results from the enhanced luminosity of the accretion disk due to the increase in its radius by up to ??0.44a ₀ at the outburst maximum (a ₀ is the component separation), as well as a shallower radial temperature decrease law than in the canonical case. As the superoutburst faded, the disk radius decreased smoothly at the end of our observation (to ??0.33a ₀), the thickness of its outer edge and temperature of its boundary layer decreased, and the parameter ?? _g approached its canonical value. Deviations from the mean brightness of the system as a function of the superhump period P _sh are analyzed for each out-of-eclipse set of observations. Various factors affecting the appearance and amplitudes of superhumps in the orbital light curves are considered.     

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.     

Light-curve synthesis for close binaries: An extended shock as an analog of a hot spot in cataclysmic variables

An algorithm is presented for the synthesis of the light curve of a close binary system consisting of a red dwarf that fills its Roche lobe and a spherical white dwarf. The spherical component is surrounded by an elliptical, geometrically thick accretion disk. The code models an extended shock located along the edge of the stream near the outer boundary of the disk. The observational manifestations of the shock show that it can be considered as an analog of a hot spot at the edge of the disk. Synthetic light curves for the SU UMa system OY Car at various phases of its activity indicate that the model can describe both typical and peculiar light curves for this cataclysmic variable reasonably well.     

Close binaries containing Supermassive Black Holes

We consider the evolution of binary systems formed by a Supermassive Black Hole (SMBH) residing in the center of a galaxy or a globular cluster and a star in its immediate vicinity. The star is assumed to fill its Roche lobe, and the SMBH accretes primarily the matter of this star. The evolution of such a system is mainly determined by the same processes as for an ordinary binary. The main differences are that the donor star is irradiated by hard radiation emitted during accretion onto the SMBH; in a detached system, nearly all the donor wind is captured by the black hole, which strongly affects the evolution of the semi-major axis; it is not possible for companions of the most massive SMBHs to fill their Roche lobes, since the corresponding orbital separations are smaller than the radius of the last stable orbit in the gravitational field of the SMBH. Moreover, there may not be efficient exchange between the orbital angular momentum and the angular momentum of the overflowing matter in such systems. Our computations assumed that, if the characteristic timescale for mass transfer is smaller than the thermal timescale of the star, no momentum exchange occurs. Absorption of incident external radiation in the stellar envelope was treated using the same formalism that was used when computing the radiative transfer in the stellar atmosphere. Numerical simulations show that Roche-lobe overflow is possible for a broad range of initial system parameters. The evolution of semi-detached systems containing a star and a SMBH nearly always ends with the dynamical disruption of the star. Stars with masses close to the solar mass are destroyed immediately after they fill their Roche lobes. During the accretion of matter of disrupted stars, the SMBH can achieve quasar luminosities. If the SMBH accretes ambient gas as well as gas stripped from stars, the star is subject to additional radiation in the detached phase of its evolution, strengthening its stellar wind. This leads to an increase of the semi-major axis and subsequent decrease of the probability of Roche-lobe overflow during the subsequent evolution of the system.     

Correlations between the development of a flare in the Blazar 3C 454.3 in the radio and optical

Radio and optical data are used to analyze the development of the flare in the blazar 3C 454.3 observed in 2004–2007. A detailed correspondance between the optical and radio flares is established, with a time delay that depends on the observing frequency. The variation of the delay of the radio flare relative to the optical flare is opposite to the dispersion delay expected for the propagation of radiation in the interstellar medium, testifying to an intrinsic origin for the observed outburst. Small-scale flux variations on time intervals of 5–10 days in the millimeter and optical are also correlated, with a time delay of about ten months. This may provide evidence for a single source generating the radiation at all wavelengths. Rapid flux fluctuations in the radio and optical that are correlated with the indicated time delays could be associated with inhomogeneities in the accretion disk. Detailed studies of the flux variations of Active Galactic Nuclei (AGN) can be used to analyze the structure of the accretion disk. A model for the energy release in AGN that is not associated purely with accretion onto supermassive black holes is proposed. As is the case for other active members of the AGN family, estimates of the lifetime of the binary black-hole system in 3C 454.3 suggest that this object is in a stage of its evolution that is fairly close to the coalescence of its black holes. The energy that is released as the companion of the central black hole loses orbital angular momentum is sufficient to explain the observed AGN phenomena. The source of primary energy release could be heating of the gas behind shock fronts that arise due to the friction between the companion black hole and the ambient gaseous medium. The orbit of the companion could be located at the periphery of the accretion disk of the central body at its apocenter and plunge more deeply into the accretion disk at its pericenter, inducing flares at all wavelengths. Energy-release parameters such as the temperature and density of the heated gas are estimated for 3C 454.3. The model considered assumes omnidirectional radiation of the medium in the presence of a magnetic field. The radiation corresponding to the minimum flux level (base level) could represent omnidirectional radiation due to the orbit of the moving companion. The fraction of the energy that is transferred to directed jets is small, comprising 1–2% of the total energy released due to the loss of orbital angular momentum by the companion.     

AE Aquarii represents a new subclass of Cataclysmic Variables

We analyze properties of the unique nova-like star AE Aquarii identified with a close binary system containing a red dwarf and a very fast rotating magnetized white dwarf. It cannot be assigned to any of the three commonly adopted sub-classes of Cataclysmic Variables: Polars, Intermediate Polars, and Accreting non-magnetized White Dwarfs. Our study has shown that the white dwarf in AE Aqr is in the ejector state and its dipole magnetic moment is ???1.5 × 10³⁴ G cm³. It switched into this state due to intensive mass exchange between the system components during a previous epoch. A high rate of disk accretion onto the white dwarf surface resulted in temporary screening of its magnetic field and spin-up of the white dwarf to its present spin period. Transition of the white dwarf to the ejector state had occurred at a final stage of the spin-up epoch as its magnetic field emerged from the accreted plasma due to diffusion. In the frame of this scenario AE Aqr represents a missing link in the chain of Polars evolution and the white dwarf resembles a recycled pulsar.     

Nonconservative evolution of cataclysmic binaries

We present a mechanism to take into account angular-momentum loss in binary systems with non-conservative mass transfer. In a number of cases, mass loss in the system can increase the orbital angular momentum of the stars. Including this mechanism in evolutionary models substantially expands the domain of stable mass transfer in binary systems. All observed cataclysmic binaries with known component masses fall within the calculated area for stable mass transfer.     

Formation of planetary systems and brown dwarfs around single stars

The conditions for the formation of planets and brown dwarfs around single main-sequence stars are considered in two scenarios. The formation of planets and brown dwarfs requires that the initial specific angular momentum of a solar-mass protostar be (0.32)×10¹⁸ cm²/s. The accreted matter of the protostar envelope forms a compact gas ring (disk) around the young star. If the viscosity of the matter in this ring (disk) is small, increasing its mass above a certain limit results in gravitational instability and the formation of a brown dwarf. If the viscosity of the gas is sufficiently large, the bulk of the protostar envelope material will be accreted by the young star, and the gas disk will grow considerably to the size of a protoplanetary dust disk due to the conservation of angular momentum. The formation of dust in the cool part of the extended disk and its subsequent collisional coalescence ultimately results in the formation of solar-type planetary systems.     

Formation of planets during the evolution of single and binary stars

Current views of the origin and evolution of single and binary stars suggest that the planets can form aroundmain-sequence single and binary stars, degenerate dwarfs, neutron stars, and stellarmass black holes according to several scenarios. Planets can arise during the formation of a star mainly due to excess angular momentum leading to the formation of an accretion-decretion disk of gas and dust around a single star or the components of a binary. It is the evolution of such disks that gives rise to planetary systems. A disk can arise around a star during its evolution due to the accretion of matter from dense interstellar clouds of gas and dust onto the star, the accretion of mass froma companion in a binary system, and the loss of matter during the contraction of a rapidly rotating star, in particular, if the star rotates as a rigid body and the rotation accelerates with its evolution along the main sequence. The fraction of stars with planetary systems is theoretically estimated as 30–40%, which is close to the current observational estimate of ∼34%.     

Evolution of Stars Paired with Intermediate-Mass Black Holes

Under certain conditions, stars close to intermediate-mass black holes (IMBHs) can form close binary systems with these objects, in which the Roche lobe can be filled by the star and intense accretion of the star’s matter onto the IMBH is possible. Recently, accreting IMBHs have been associated with hyperluminous X-ray sources (HLXs), whose X-ray luminosities can exceed 10⁴¹ erg/s. In this paper, the evolution of star—IMBH binary systems is investigated assuming that the IMBH mainly accretes the matter of its companion star, and that the presence of gas in the vicinity of the IMBH does not appreciably affect changes in the orbit of the star. The computations take into account all processes determining the evolution of ordinary binary systems, as well as the irradiation of a star by hard radiation during the accretion of its matter onto the IMBH. The absorption of external radiation in the stellar envelope was calculated applying the same formalism that is used to calculate the opacity of the stellar matter. The computations also assumed that, if the characteristic time for the mass transfer is less than the thermal time scale of the star, there is no exchange betwween the orbital angular momentum of the system and the angular momentum of the matter flowing onto the IMBH.

Numerical simulations have shown that, under these assumptions, three types of evolution are possible for such a binary system, depending on the mass of the IMBH and the star, as well as on the star’s initial distance from the IMBH. The first type ends with the destruction of the star. For low-mass main sequence (MS) stars, only this option is realized, even in the case of large initial distances from IMBH. For massive MS stars, the star is also destroyed if the mass of the IMBH is high and the initial distance of the star from the IMBH is sufficiently small.

The second type of evolution can occur for massive MS stars, which are initially located farther from the IMBH than in the first type of evolution. In this case, the massive star fills its Roche lobe during its evolutionary expansion, after which a stage of intense mass transfer begins. It is in this phase of the evolution that the star- IMBH system can manifest itself as a HLX, when its X-ray luminosity L_X exceeds 10⁴¹ erg/s for a fairly long time. Numerical simulations show that the initial mass of the donor star in systems with M_BH = (10³?10⁵)M_⊙ must be close to ～10 M_⊙ in this case. The characteristic duration of the HLX stage is 30 000–70 000 years. For smaller initial donor masses close to ～5M_⊙, L_X does not reach 10⁴¹ erg/s in the stage of intense mass transfer, but can exceed 10⁴⁰ erg/s. The duration of this stage of evolution is 300 000–800 000 years. A characteristic feature of this second type of evolution is an increase in the orbital period of the system over time. As a result, after a period of intense mass loss, the star “retreats” inside the Roche lobe. A remnant of the star in the form of a white dwarf is left behind, and can end up fairly far from the IMBH.

The third type of evolution can occur for massive MS stars that are initially even farther from the IMBH, as well as for massive stars that are already evolved at the initial time. In this case, conservative mass exchange in the presence of intense stellar wind leads to the star moving away from the IMBH, without filling its Roche lobe at all. For massive stars with sufficiently strong stellar winds (for example, stars with masses ～50M_⊙), the accretion rate of matter onto the IMBH in this case can reach values that are characteristic of HLXs. As in the case of the second type of evolution, the stellar remnant can remain at a fairly large distance from the IMBH.

     

Structure of the accretion disk in the cataclysmic variable UX UMa in the transition to its quiescent state

Our analysis of BV RI light curves for the cataclysmic variable UX UMa obtained in intermediate activity states, in the transition between the active and quiescent states of the system on March 12, 1997 and May 3, 2000, shows that the shapes of these light curves cannot be adequately described using the standard hot-spot model. A model with a “hot line” near the edge of the disk and a two-armed spiral structure on the disk surface reproduces much better out-of-eclipse variations in the light curves. The presence of an extended hot line can explain the anomalous shape of the I light curve on March 12, 1997. The decrease in the observed luminosity of the system between March 12, 1997 and May 3, 2000 could be due to a decrease in the disk luminosity by a factor of 2–2.5; the higher disk luminosity on the earlier date is associated with appreciable deviations of the radial temperature distribution of the disk material from that for the standard model. The phases and depths of dips in the out-of-eclipse sections of the UX UMa light curves are due primarily to the parameters of the complex shape of the accretion disk, which has a spiral structure located mainly near its outer edge. The contribution of the spiral arms in the V filter reaches 20–50% of the total disk radiation. The crest of the first spiral wave in our model maintains its approximate position in azimuth; this structure could represent a bulge in a halo at the outer edge of the disk near orbital phases φ ～ 0.7, in the direction of the continuation of the extended shock in the disk itself. The position of the crest of the second spiral arm changes with time. This structure may represent a one-armed spiral wave near the apastron of the weakly elliptical disk. Finally, the observations testify to the presence of another spiral arm that is les clearly manifest in terms of both its luminosity and its height above the unperturbed disk surface. Thus, in an intermediate activity state of UX UMa, the surface of the accretion disk is distorted by the action of a two-armed spiral structure in the outer regions of the disk, which is asymmetric in both its luminosity and dimensions, and a bulge at the disk edge in the region of its interaction with the inflow to the disk.     

Meridional circulation in young massive stars with laminar rotation

We study the rotation of a chemically homogeneous star with a mass of 16M_⊙, assuming that the angular-momentum distribution in its radiative envelope is determined by hydrodynamical processes—flows and turbulent diffusion. Meridional circulation and horizontal shear turbulence are the main hydrodynamical processes forming the radial distribution of the angular momentum in young massive stars in the absence of magnetic fields. The rotation of such stars is close to steady-state. The angular velocity of rotation of the convective core can be ～5–20% higher than the surface value. Under these conditions, the characteristic time for the radial transport of angular momentum by meridional flows and shear turbulence is comparable to the nuclear time scale.     

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.     

A photometric study of the eclipsing dwarf nova GY Cnc in quiescence and during an outburst

The results of photometric observations of the dwarf nova GY Cnc in the Rc filter acquired in 2013–2015 (~3900 orbital cycles, 19 nights in total) are presented, including observations during its outburst in April 2014. The binary’s orbital elements have been refined. The orbital period has changed only insignificantly during the ~30 000P_orb since the earlier observations; no systematic O–C variations were detected, only fluctuations within 0.004^d on time scales of 1500–2000P_orb. A “combined” model is used to solve for the parameters of GY Cnc during two states of the system. The flux from the white dwarf is negligible due to the star’s small size. The temperature of the donor star, T₂ ~ 3667 K (Sp M0.2 V), varies between 3440 and 3900 K (Sp K8.8–M1.7 V). The semi-major axis of the disk is a ~ 0.22a₀, on average. In quiescence, a varies within ~40%. The disk has a considerable eccentricity (e ~ 0.2?0.3) for a < 0.2a₀. The disk shape becomes more circular (e < 0.1) with increasing a. The outburst of GY Cnc was associated with increased luminosity of the disk due to the parameter α_g (related to the viscosity of the disk material) decreasing to 0.1–0.2 and the temperature in the inner parts of the disk increasing twofold, to T_in ~ 95 000 K. These changes were apparently due to the infall of matter onto the surface of the white dwarf as the outburst developed. All parameters of the accretion disk in quiescence display considerable variations about their mean values.     

The vorticity and angular momentum budgets of Asian summer monsoon

The study delineates the vorticity and angular momentum balances of Asian summer monsoon during the evolution and established phases. It also elucidates the differences between these balances in the National Centre for Environmental Prediction/National Centre for Atmospheric Research (NCEP/NCAR) reanalysis and the National Centre for Medium Range Weather Forecasts (NCMRWF) analysis fields. The NCEP/NCAR reanalysis for a 40 year period (1958-97) and the NCMRWF analysis for a three year (1994-96) period are made use of for the purpose. The time mean summer monsoon circulation is bifurcated into stable mean and transient eddy components and the mean component is elucidated. The generation of vorticity due to stretching of isobars balances most of the vorticity transported out of the monsoon domain during the evolution period. However, during the established period, the transportation by the relative and planetary vorticity components exceeds the generation due to stretching. The effective balancing mechanism is provided by vorticity generation due to sub-grid scale processes. The flux convergence of omega and relative momenta over the monsoon domain is effectively balanced by pressure torque during the evolution and established phases. Nevertheless, the balance is stronger during the established period due to the increase in the strength of circulation. Both the NCMRWF and NCEP fields indicate the mean features related to vorticity and angular momentum budgets realistically. Apart from the oceanic bias (strong circulation over oceans rather than continents), the summer monsoon circulation indicated by the NCEP is feeble compared to NCMRWF. The significant terms in the large-scale budgets of vorticity and angular momentum enunciate this aspect     

Photometry of the intermediate polar 1RXS J062518.2+733433: Characteristics of its periodic brightness variations

We present the results of our photometry of the recently discovered intermediate polar 1RXS J062518.2+733433. The observations were made using the 70 cm telescope of the Astronomical Observatory of the Ural State University with a multichannel photometer, and were carried out during 13 nights in March–May 2004. Our analysis reveals brightness variations with periods of 19.788 ± 0.003 min, corresponding to the rotational period of the white dwarf, 21.273 ± 0.003 min, corresponding to the orbital sideband, and 283.3 ± 0.5 min, which is the orbital period. The variations with the white dwarf’s rotational period show a stable amplitude and a quasi-sinusoidal pulse shape persistent over a long time. In contrast, the orbital-sideband variations have an unstable amplitude and a significantly nonsinusoidal pulse shape that varies with time. Variations of the amplitude and pulse shape of the orbital-sideband variations can be related to structural changes of the accretion disk surrounding the white dwarf.     

Accretion disk of the cataclysmic variable IY UMa in quiescence and its active state

We present the results of a study of the eclipsing binary IY UMa (type SU UMa) during its superoutburst of 2004 and in quiescence.We have refined the orbital period of the system. Light curves are presented for various states of activity.We estimate the parameters of IYUMa during the superoutburst and in quiescence for hot-line and spiral-wave models. The spiral-wave model, which takes into account the presence of vertical perturbations in outer parts of the disk, is able to reproduce the light-curve shapes and phases of dips in the out-of-eclipse parts of the binary’s light curves during the outburst both qualitatively and quantitatively.     

Stability of the multiple star system ι UMa (ADS 7114)

The physical and orbital parameters of the quadruple star system ι UMa (HD 76644 = ADS 7114) were determined earlier, when it was concluded based on modeling the system’s dynamics and applying theoretical stability criteria that the system was probably unstable. Here the stability of the ι UMa system is studied by calculating the Lyapunov characteristic exponents for representative sets of parameters and initial conditions. The conclusions on the system’s stability (or instability) based on various stability criteria and the calculated Lyapunov exponents are compared. The instability of the system as a whole is confirmed rigorously based on massive computations of the Lyapunov exponents. This system appears to be the only known multiple system whose instability has been rigorously established. The Lyapunov time-disruption time statistical relations are constructed, which show that the Hamiltonian intermittency of the second kind dominates. Typical disruption times are shorter than 1000 years, and the Lyapunov times are shorter than 100 years.