Modeling the luminosity function of galactic low-mass X-ray binaries

The evolution of the family of binaries with a low-mass star and a compact neutron star companion (low-mass X-ray binaries (LMXBs) with neutron stars) ismodeled by the method of population synthesis. Continuous Roche-lobe filling by the optical star in LMXBs is assumed to be maintained by the removal of orbital angular momentum from the binary by a magnetic stellar wind from the optical star and the radiation of gravitational waves by the binary. The developed model of LMXB evolution has the following significant distinctions: (1) allowance for the effect of the rotational evolution of a magnetized compact remnant on themass transfer scenario in the binary, (2) amore accurate allowance for the response of the donor star to mass loss at the Roche-lobe filling stage. The results of theoretical calculations are shown to be in good agreement with the observed orbital period-X-ray luminosity diagrams for persistent Galactic LMXBs and their X-ray luminosity function. This suggests that the main elements of binary evolution, on the whole, are correctly reflected in the developed code. It is shown that most of the Galactic bulge LMXBs at luminosities L _x > 10³⁷ erg s^?1 should have a post-main-sequence Roche-lobe-filling secondary component (low-mass giants). Almost all of the models considered predict a deficit of LMXBs at X-ray luminosities near ～10^36.5 erg s^?1 due to the transition of the binary from the regime of angular momentum removal by a magnetic stellar wind to the regime of gravitational waves (analogous to the widely known period gap in cataclysmic variables, accreting white dwarfs). At low luminosities, the shape of the model luminosity function for LMXBs is affected significantly by their transient behavior-the accretion rate onto the compact companion is not always equal to the mass transfer rate due to instabilities in the accretion disk around the compact object. The best agreement with observed binaries is achieved in the models suggesting that heavy neutron stars with masses 1.4–1.9M _⊙ can be born.     

A model of low-mass neutron stars with a quark core

We consider an equation of state that leads to a first-order phase transition from the nucleon state to the quark state with a transition parameter λ>3/2 (λ=ρ_Q/(ρ_N+P₀/c²)) in superdense nuclear matter. Our calculations of integrated parameters for superdense stars using this equation of state show that on the stable branch of the dependence of stellar mass on central pressure dM/dP_c>0) in the range of low masses, a new local maximum with M_max=0.082_⊙ and R=1251 km appears after the formation of a toothlike kink (M=0.08M_⊙, R=205 km) attributable to quark production. For such a star, the mass and radius of the quark core are M_core=0.005M_⊙ and R_core=1.73 km, respectively. In the model under consideration, mass accretion can result in two successive transitions to a quark-core neutron star with energy release similar to a supernova explosion: initially, a low-mass star with a quark core is formed; the subsequent accretion leads to configurations with a radius of ～1000 km; and, finally, the second catastrophic restructuring gives rise to a star with a radius of ～100 km.     

The gas content in galactic disks: Correlation with kinematics

We consider the relationship between the total HI mass in late-type galaxies and the kinematic properties of their disks. The mass M_HI for galaxies with a wide variety of properties, from dwarf dIrr galaxies with active star formation to giant low-brightness galaxies, is shown to correlate with the product V_cR₀ (V_c is the rotational velocity, and R₀ is the radial photometric disks cale length), which characterizes the specific angular momentum of the disk. This correlation, along with the decrease in the relative mass of the gas in a galaxy with increasing V_c, can be explained in terms of the previous assumption that the gas density in the disks of most galaxies is maintained at a level close to the threshold (marginal) stability of a gaseous layer to local gravitational perturbations. In this case, the regulation mechanism of the star formation rate associated with the growth of local gravitational instability in the gaseous layer must play a crucial role in the evolution of the gas content in the galactic disk.     

Radius of the neutron star magnetosphere during disk accretion

The dependence of the spin frequency derivative \(\dot \nu \) of accreting neutron stars with a strongmagnetic field (X-ray pulsars) on the mass accretion rate (bolometric luminosity, L_bol) has been investigated for eight transient pulsars in binary systems with Be stars. Using data from the Fermi/GBM and Swift/BAT telescopes, we have shown that for seven of the eight systems the dependence \(\dot \nu \) (L_bol) can be fitted by the model of angular momentum transfer through an accretion disk, which predicts the relation \(\dot \nu \) ～ L^6/7_bol. Hysteresis in the dependence \(\dot \nu \) (L_bol) has been confirmed in the system V 0332+53 and has been detected for the first time in the systems KS 1947+300, GRO J1008-57, and 1A 0535+26. Estimates for the radius of the neutron star magnetosphere in all of the investigated systems have been obtained. We show that this quantity varies from pulsar to pulsar and depends strongly on the analytical model and the estimates for the neutron star and binary system parameters.     

Evolution of low-mass helium stars in semidetached binaries

We systematically investigate the evolution of low-mass (0.35, 0.40, and 0.65M _⊙) helium donors in semidetached binaries with white-dwarf accretors. The initial periods of the binaries are chosen in such a way that the helium abundance in the center of the models at the time of Roche lobe overflow varies between Y _c = 0.98 and Y _c ? 0.1. The results of our calculations can be used to analyze the formation scenarios and evolutionary status of AM CVn stars. We show that the minimum orbital periods of the semidetached binaries depend weakly on the total mass of the components and the evolutionary phase of the donor at the time of Roche lobe overflow and are 9–10 min. The differences in the mass transfer rates after P _orb reaches its minimum in the range P _orb ≈ 10–40 min do not exceed a factor of ～2.5. For P _orb ? 20 min, the mass-losing stars are weakly degenerate homogeneous cooling objects; the He, C, N, O, and Ne abundances depend on the evolutionary phase at which Roche lobe overflow occurred. For the binaries that are currently believed to be the most probable candidates for AM CVn stars with helium donors, Y ? 0.4, X _C ? 0.3, X _O ? 0.25, and X _N ? 0.5 × 10^?2. In the binaries under consideration, once P _orb ≈ 40 min has been reached, the mass loss time scale begins to exceed the thermal time scale of the donors, the latter begin to contract, their matter becomes degenerate, and the populations of AMCVn stars with white-dwarf and helium-star progenitors of their donors probably merge together.     

A Model of the Mira–Type Star T UMi

Stellar evolution calculations were carried out from the main sequence to the final stage of the asymptotic giant branch for stars with initial masses 1 M_⊙ ≤ M_ZAMS ≤ 2 M_⊙ and metallicity Z = 0.01. Selected models of evolutionary sequences were used as initial conditions for solution of the equations of radiation hydrodynamics and time–dependent convection describing radial stellar pulsations. The study was aimed to construct the hydrodynamic models of Mira–type stars that show the secular decrease in the pulsation period Π commenced in 1970th at Π = 315 day. We show that such a condition for the period change is satisfied with evolutionary sequences 1 M_⊙ ≤ M_ZAMS ≤ 1.2 M_⊙ and the best agreement with observations is obtained for M_ZAMS = 1.2 M_⊙. The pulsation period reduction is due to both the stellar radius decrease during the thermal pulse of the helium burning shell and mode switch from the fundamental mode to the first overtone. Theoretical estimates of the fundament parameters of the star at the onset of pulsation period reduction are as follows: the mass is M = 0.93 M_⊙, the luminosity is L = 4080 L_⊙, and the radius is R = 220 R_⊙. The mode switch occurs 35 years after the onset of period reduction.     

Einstein–Maxwell Field Equation in Isotropic Coordinates: An Application to Modeling Superdense Star

We present a charged analogue of Pant et al. (2010, Astrophys. Space Sci., 330, 353) solution of the general relativistic field equations in isotropic coordinates by using simple form of electric intensity E that involve charge parameter K. Our solution is well behaved in all respects for all values of X lying in the range 0 <X≤ 0.11, K lying in the range 4 <K≤ 6.2 and Schwarzschild compactness parameter u lying in the range 0 <u≤ 0.247. Since our solution is well behaved for wide ranges of the parameters, we can model many different types of ultra-cold compact stars like quark stars and neutron stars. We have shown that corresponding to X = 0.077 and K = 6.13 for which u = 0.2051 and by assuming surface density ρ _b =4.6888×10¹⁴ g cm ^?3 the mass and radius are found to be 1.509M _⊙, 10.906 km respectively which match with the observed values of mass 1.51M _⊙ and radius 10.90 km of the quark star XTE J1739-217. The well behaved class of relativistic stellar models obtained in this work might have astrophysical significance in the study of more realistic internal structures of compact stars.     

Flux-Vortex Pinning and Neutron Star Evolution

G. Srinivasan et al. (1990) proposed a simple and elegant explanation for the reduction of the neutron star magnetic dipole moment during binary evolution leading to low mass X-ray binaries and eventually to millisecond pulsars: Quantized vortex lines in the neutron star core superfluid will pin against the quantized flux lines of the proton superconductor. As the neutron star spins down in the wind accretion phase of binary evolution, outward motion of vortex lines will reduce the dipole magnetic moment in proportion to the rotation rate. The presence of a toroidal array of flux lines makes this mechanism inevitable and independent of the angle between the rotation and magnetic axes. The incompressibility of the flux-line array (Abrikosov lattice) determines the epoch when the mechanism will be effective throughout the neutron star. Flux vortex pinning will not be effective during the initial young radio pulsar phase. It will, however, be effective and reduce the dipole moment in proportion with the rotation rate during the epoch of spindown by wind accretion as proposed by Srinivasan et al. The mechanism operates also in the presence of vortex creep.     

Magnetic and spin evolution of neutron stars in close binaries

The evolution of neutron stars in close binary systems with a low-mass companion is considered, assuming the magnetic field to be confined within the solid crust. We adopt the standard scenario for the evolution in a close binary system, in which the neutron star passes through four evolutionary phases ('isolated pulsar'–'propeller'– accretion from the wind of a companion – accretion resulting from Roche-lobe overflow). Calculations have been performed for a great variety of parameters characterizing the properties of both the neutron star and the low-mass companion. We find that neutron stars with more or less standard magnetic field and spin period that are processed in low-mass binaries can evolve to low-field rapidly rotating pulsars. Even if the main-sequence life of a companion is as long as 10¹⁰ yr, the neutron star can maintain a relatively strong magnetic field to the end of the accretion phase. The model that is considered can account well for the origin of millisecond pulsars.     

On a family of well behaved perfect fluid balls as astrophysical objects in general relativity

A family of well behaved perfect fluid balls has been derived starting with the metric potential g ₄₄=B(1+Cr ²) ⁿ for all positive integral values of n. For n≥4, the members of this family are seen to satisfy the various physical conditions e.g. c ² ρ≥p≥0,dp/dr<0,dρ/dr<0, along with the velocity of sound \((\sqrt{dp/c^{2}d\rho} )< 1\) and the adiabatic index ((p+c ² ρ)/p)(dp/(c ² dρ))>1. Also the pressure, energy density, velocity of sound and ratio of pressure and energy density are of monotonically decreasing towards the pressure free interface (r=a). The fluid balls join smoothly with the Schwarzschild exterior model at r=a. The well behaved perfect fluid balls so obtained are utilised to construct the superdense star models with their surface density 2×10¹⁴  gm/cm³. We have found that the maximum mass of the fluid balls corresponding to various values of n are decreasing with the increasing values of n. Over all maximum mass for the whole family turns out to be 4.1848M _Θ and the corresponding radius as 19.4144 km while the red shift at the centre and red shift at surface as Z ₀=1.6459 and Z _a =0.6538 respectively this all happens for n=4. It is interesting to note that for higher values of n viz n≥170, the physical data start merging with that of Kuchowicz superdense star models and hence the family of fluid models tends to the Kuchowicz fluid models as n→∞. Consequently the maximum mass of the family of solution can not be less than 1.6096 M _Θ which is the maximum mass occupied by the Kuchowicz superdense ball. Hence each member of the family for n≥4 provides the astrophysical objects like White dwarfs, Quark star, typical neutron star.     

Instability of radial modes at the Wolf-Rayet evolutionary stage

The evolution of a Population-I star with an initial mass M _ZAMS = 60 M _⊙ has been calculated. At the stage when a red giant turns into an early-type helium star, the vast bulk of the stellar mass is concentrated in a compact core surrounded by an extended envelope that is unstable with respect to radial oscillations. The range of effective temperatures within which the instability arises extends to T _eff ? 10⁵ K. For the models corresponding to the Wolf-Rayet evolutionary stage (5 × 10⁴ K ≤ T _eff ≤ 1.05 × 10⁵ K), hydrodynamic calculations of self-exciting radial stellar pulsations have been performed. The pulsational instability develops in a time interval comparable to the dynamic timescale. Once the amplitude has ceased to grow, the pulsational motions are nonlinear traveling waves propagating from the core boundary to the stellar surface. The velocity amplitude of the outer layers is 500 km s^?1 < ΔU < 10³ km s^?1, depending on the effective temperature. During the evolution of a helium star, the mean ratio of the maximum expansion velocity of the outer layers to the local escape velocity decreases and lies within the range 0.25 < U _max/v _esc < 0.6 for the models considered. The nonlinearity of the stellar pulsations is responsible for the increase in the mean radius \(\bar r\) of the Lagrangian layers compared to the equilibrium radius r _eq. The effect of the increase in mean radius decreases with rising effective temperature from\(\bar r\)/r ～ 10 at T _eff = 7 × 10⁴ K to \(\bar r\)/r ≈ 2 at T _eff = 10⁵ K. The radial pulsation periods for the models considered lie within the range 0.1 day ≤ Π ≤ 1.6 day and the amplitude of the bolometric magnitude variations does not exceed 0 _. ^m 2.     

Bimodality of anomalous pulsars?

At present, it is widely believed that anomalous X-ray pulsars (AXPs), soft gamma-ray repeaters (SGRs), rotational radio transients (RRATs), compact central objects (CCOs) in supernova remnants, and X-ray dim isolated neutron stars (XDINSs) belong to different classes of anomalous objects in which the central bodies are isolated neutron stars. Previously, we have shown that AXPs and SGRs can be described in terms of the drift model for parameters of the central neutron star typical of radio pulsars (rotation periods P ～ 0.1–1 s and surface magnetic fields B ～ 10¹¹–10¹³ G). Here, we show that some of the peculiarities of the sources under consideration can be explained by their geometry (in particular, by the angle β between the rotation axis and the magnetic moment). If β ? 10° (an aligned rotator), the drift waves in the outer layers of the neutron star magnetosphere can account for the observed periodicity in the radiation. For large β (a nearly orthogonal rotator), the observed modulation of the radiation and its short bursts can be explained by mass accretion from the ambient medium (e.g., a relic disk).     

Outskirts of Galaxy Clusters A1139, A1314, A1656, A2040, A 2052, A2107: Star-Formation Rate

We investigate the variation of the fraction of galaxies with suppressed star formation (M_K < ?21 _. ^m 5) and early-type galaxies (frac_E) of the “red sequence” along the projected radius in six galaxy clusters:Coma (A1656), A1139, and A1314 in the Leo supercluster region (z ≈ 0.037) and A2040, A2052, A2107 in the Hercules supercluster region (z ≈ 0.036). According to SDSS (DR10) data, frac_E is the highest in the central regions of galaxy clusters and it is, on the average, equal to 0.62 ± 0.03, whereas in the 2–3R/R_200c interval and beyond the R_sp ≈ 0.95 ± 0.04 R_200m radius that we inferred from the observed profile frac_E is minimal and equal to 0.25 ± 0.02. This value coincides with the estimate frac_E = 0.24 ± 0.01 that we inferred for field galaxies located between the Hercules and Leo superclusters at the same redshifts. We show that the fraction of galaxies with suppressed star formation decreases continuously with cluster radius from 0.87 ± 0.02 in central regions down to 0.43 ± 0.03 in the 2–3 R/R_200c interval and beyond R_sp, but remains, on the average, higher than 26% than the corresponding fraction for field objects. This decrease is especially conspicuous in the galaxy mass interval log M_* [M_⊙] = 9.5–10. We found that galaxies with ongoing star formation have average clustercentric distances 1.5–2.5 R/R_200c and that their radial-velocity dispersions are higher than those of galaxies with suppressed star formation.     

Change in the orbital period of a binary system due to dynamical tides for main-sequence stars

We investigate the change in the orbital period of a binary system due to dynamical tides by taking into account the evolution of a main-sequence star. Three stars with masses of one, one and a half, and two solar masses are considered. A star of one solar mass at lifetimes t = 4.57 × 10⁹ yr closely corresponds to our Sun. We show that a planet of one Jupiter mass revolving around a star of one solar mass will fall onto the star in the main-sequence lifetime of the star due to dynamical tides if the initial orbital period of the planet is less than P _orb ≈ 2.8 days. Planets of one Jupiter mass with an orbital period P _orb ≈ 2 days or shorter will fall onto a star of one and a half and two solar masses in the mainsequence lifetime of the star.     

Type-Ia supernovae in semidetached binaries

We consider the evolutionary scenarios for close binaries that lead to the formation of semidetached systems in which a white dwarf can accumulate the Chandrasekhar mass through mass accretion from its companion, a main sequence star or a subgiant of mass M ～ 2M_⊙. Such dwarfs probably explode as type-Ia supernovae or collapse to form a neutron star. The population synthesis method is used to analyze the dependence of the model rate of these events in the Galaxy on the common envelope parameter, the mass transfer rate, and the response of a main-sequence star to helium accretion at an intermediate evolutionary stage. The rate of explosions in semidetached systems of this type in the Galaxy was found to be no higher than ?0.2×10^?3 yr^?1, which is less than 10% of the lower level for the empirically estimated SNe Ia rate.     

Symbiotic halo star LT Del: Phase variations of the emission spectrum and parameters of the cool component

Long-term photometric and spectroscopic observations of the yellow symbiotic star LT Del are analyzed. UBV light curves are presented. Based on the observations of 20 cycles, we have refined the orbital period of the star, P = 476 _· ^d 0 ± 1 _· ^d 0. The brightness has been found to be unstable at some orbital phases with an amplitude up to 0 _· ^m 3. We have measured the fluxes in hydrogen and helium emission lines and in continuum and investigated their relationship to the orbital period. The fluxes in hydrogen and HeI lines follow the UBV light curves in phase; the He II 4686 Å flux does not depend on the phase and is constant within the accuracy of our measurements. The intensity ratio of the 4686 Å andHβ lines changes from 0.2 to 0.9 over the period. We interpret the spectroscopic observations based on the hypothesis of heating and ionization of the stellar wind from a cool component by high-frequency radiation from a hot star with a temperature of 105 K. We have estimated the spectral type of the cool star from our photometry and its continuum energy distribution as a bright K2–4 red giant branch halo star. The bolometric luminosity and mass loss rate have been estimated for the K component to be L _bol ～ 700L _⊙ and \(\dot{M}\) ～ 10^?8 M _⊙ yr^?1, respectively.     

A hydrodynamic model for asymmetric explosions of rapidly rotating collapsing supernovae with a toroidal atmosphere

We numerically solved the two-dimensional axisymmetric hydrodynamic problem of the explosion of a low-mass neutron star in a circular orbit. In the initial conditions, we assumed a nonuniform density distribution in the space surrounding the collapsed iron core in the form of a stationary toroidal atmosphere that was previously predicted analytically and computed numerically. The configuration of the exploded neutron star itself was modeled by a torus with a circular cross section whose central line almost coincided with its circular orbit. Using an equation of state for the stellar matter and the toroidal atmosphere in which the nuclear statistical equilibrium conditions were satisfied, we performed a series of numerical calculations that showed the propagation of a strong divergent shock wave with a total energy of ～0.2×10⁵¹ erg at initial explosion energy release of ～1.0×10⁵¹ erg. In our calculations, we rigorously took into account the gravitational interaction, including the attraction from a higher-mass (1.9M_⊙) neutron star located at the coordinate origin, in accordance with the rotational explosion mechanism for collapsing supernovae. We compared in detail our results with previous similar results of asymmetric supernova explosion simulations and concluded that we found a lower limit for the total explosion energy.     

Effect of partial mixing of matter on the hydrodynamic angular momentum transport processes in massive main-sequence stars

We consider the evolution of a rotating star with a mass of 16M _⊙ and an angular momentum of 3.25 × 10⁵² g cm² s^?1, along with the hydrodynamic transport of angular momentum and chemical elements in its interiors. When the partial mixing of matter of the turbulent radiative envelope and the convective core is taken into account, the efficiency of the angular momentum transport by meridional circulation in the stellar interiors and the duration of the hydrogen burning phase increase. Depending on the Schmidt number in the turbulent radiative stellar envelope, the ratio of the equatorial rotational velocity to the circular one increases with time in the process of stellar evolution and can become typical of early-type Be stars during an additional evolution time of the star on the main sequence. Partial mixing of matter is a necessary condition under which the hydrodynamic transport processes can increase the angular momentum of the outer stellar layer to an extent that the equatorial rotational velocity begins to increase during the second half of the evolutionary phase of the star on the main sequence, as shown by observations of the brightest stars in open star clusters with ages of 10–25 Myr. When the turbulent Schmidt number is 0.4, the equatorial rotational velocity of the star increases during the second half of the hydrogen burning phase in the convective core from 330 to 450 km s^?1.     

The formation and evolution of low-mass close binaries with compact components

We discuss the formation and evolution of interacting low-mass close binaries with a He-1CO- or ONe-dwarf neutron star or a black hole as a compact component. Mass exchange leads to cataclysmic events in such systems. The rate of semidetached low-mass close binary formation is 5×10^–3 yr^–1 if the accreting component is a He degenerate dwarf, 5×10^–3 yr^–1 if it is a CO-dwarf and 3×10^–8 yr^–1 if it is a neutron star. Systems with compact accretors arise as the result of the common envelope phase of close binary evolution or due to collisions of single neutron stars or dwarfs with low-mass single stars in dense stellar clusters. Evolution of LMCB to the contact phase in semi-detached stages is determined mainly by the angular momentum losses by a magnetic stellar wind and radiation of gravitational waves. Numerical computations of evolution with momentum loss explain observed mass exchange rates in such systems, the absence of cataclysmic variables with orbital periods 2^h–3^h, the low number and the evolutionary status of systems with orbital periods shorter than 80^m. In conclusion we list unsolved problems related to magnetic stellar wind, the distribution of young close binaries over main initial parameters, stability of mass exchange.Paper presented at the IAU Colloquium No. 93 on Cataclysmic Variables. Recent Multi-Frequency Observations and Theoretical Developments, held at Dr. Remeis-Sternwarte Bamberg, F.R.G., 16–19 June, 1986.