The stars in the Main Sequence and their dynamical parameters

The stars in the Main Sequence are seen as a hierarchy of objects with different massesM and effective dynamical radiiR _eff=R/α given by the stellar radii and the coefficients for the inner structure of the stars. As seen in a previous work (Paper I), during the lifetime in the Main SequenceR _eff(t) remains a near invariant when compared to the variation in the time ofR(t) and α(t). With such an effectiveR _eff one obtains the amounts of actionA _c(M), the effective densities ρ_eff(M)=ρ(M)α³(M), the densities of action and of energy (or mean presures in the stellar interior)a _c(M),e _c(M), and the potential energiesE _p(M). The amounts of action areA _c∝M ^k withk≈1.87 for the M stars,k≈5/3 for the KGF stars, andk≈1.83 for the A and earlier stars, representing very simples conditions for the other dynamical parameters. For instancek≈5/3 means a near invariant effective density α_eff for the KGF stars, while for such stars the mean densities and coefficients α present the strongest variations with masses ρ(M)∝M ^?1.81, α(M)∝M^0.6. The cases for the M stars (e _c(M)∝M ^?1) and for the A and earlier stars (betweena _c(M)=constant and α_eff(M)∝M ^?1) and also discussed. These conditions for the earlier stars also represent reasonable mean values for the whole stellar hierarchy in the range of masses 0.2M _⊙≤M≤25M _⊙. With all this, one can build ‘dynamical’ HR diagrams withA _c(M), E_p(M), α_eff M ^?p , etc., whose characteristics are analogous to these in the photometrical HR diagram. A comparison is made betweenA _c(M) from the models here and the HR diagram with the best known stars of luminosity classes IV, V, and white dwarfs. The comparison of the potential energiesE _p(M)∝M ^?p according to the stellar models used here and the observed frequency function ψ(M∝M _?q (number of stars in a given interval of masses) from different authors suggests the possibility that the productE _p(M)ψ(M) is a constant, but this must be confirmed with further studies of the function ψ(M) and its fine structure. There are analogies between the formulation used here for the stellar hierarchy and other physical processes, for instance, in modified forms of the Kolmogorov law of turbulence and in the formulation used for the hierarchy of molecular clouds in gravitational equilibrium. Besides, the function of actionA _c(M) for the stars has analogous properties to the relations of angular momenta and massesJ(M) for different types of objects. The cosmological implications of all this are discussed.     

Neutron stars in globular clusters: Formation and observational manifestations

Population synthesis is used to model the number of neutron stars in globular clusters that are observed as low-mass X-ray sources and millisecond radio pulsars. The dynamical interactions between binary and single stars in a cluster are assumed to take place only with a continuously replenished “background” of single stars whose properties keep track of the variations in parameters of the cluster as a whole and the evolution of single stars. We use the hypothesis that the neutron stars forming in binary systems from components with initial masses of ～8–12 M _⊙ during the collapse of degenerate O-Ne-Mg cores through electron captures do not acquire a high space velocity. The remaining neutron stars (from single stars with masses >8 M _⊙ or from binary components with masses >12 M _⊙) are assumed to be born with high space velocities. According to this hypothesis, a sizeable fraction of the forming neutron stars remain in globular clusters (about 1000 stars in a cluster with a mass of 5 × 10⁵ M _⊙). The number of millisecond radio pulsars forming in such a cluster in the case of accretion-driven spinup in binary systems is found to be ～10, in agreement with observations. Our modeling also reproduces the observed shape of the X-ray luminosity function for accreting neutron stars in binary systems with normal and degenerate components and the distribution of spin periods for millisecond pulsars.     

Stellar dynamics in young clusters: the formation of massive runaways and very massive runaway mergers

In the present paper we combine an N-body code that simulates the dynamics of young dense stellar systems with a massive star evolution handler that accounts in a realistic way for the effects of stellar wind mass loss. We discuss two topics.

1. The formation and the evolution of very massive stars (with masses >120 M_⊙) is followed in detail. These very massive stars are formed in the cluster core as a consequence of the successive (physical) collisions of the 10–20 most massive stars in the cluster (this process is known as ‘runaway merging’). The further evolution is governed by stellar wind mass loss during core hydrogen and core helium burning (the WR phase of very massive stars). Our simulations reveal that, as a consequence of runaway merging in clusters with solar and supersolar values, massive black holes can be formed, but with a maximum mass ≈70 M_⊙. In low-metallicity clusters, however, it cannot be excluded that the runaway-merging process is responsible for pair-instability supernovae or for the formation of intermediate-mass black holes with a mass of several 100 M_⊙.
2. Massive runaways can be formed via the supernova explosion of one of the components in a binary system (the Blaauw scenario), or via dynamical interaction of a single star and a binary or between two binaries in a star cluster. We explore the possibility that the most massive runaways (e.g. ζ Pup, λ Cep, BD+43°3654) are the product of the collision and merger of two or three massive stars.

     

Accurate masses and radii of normal stars

Summary Binary stars are the main source of fundamental data on stellar masses and radii (M, R). Considerable progress has been made in recent years in the quality and quantity of such data, and stellar masses and radii of high accuracy have led to a number of qualitatively new and interesting results on the properties and evolution of normal stars. This paper reviews the current status of fundamentalM andR determinations which (i) have errors 2%, the limit for non-trivial results in many applications, and (ii) can be presumed valid for single stars. These two conditions limit the discussion to data fromdetached, doublelined eclipsing binary systems.After a brief discussion (Sect. 2) of the main tests for accuracy and consistency which must be met for observational data to be included in the sample, data for 45 binary systems (90 single stars) are presented in Sect. 3 (Table 1 and Figs. 2–5). Spectral types are O8-M1 on the main sequence, with only two stars clearly in the red-giant region. From the review by Popper (1980), data for only 6 systems survive unchanged in the present list, while improved data are given for 18 systems; 21 systems are new additions. Broadband colours, effective temperatures, and luminosities are also given, but are scale-dependent and considerably less reliably determined thanM andR.The observed ranges inM andR for a given colour far exceed the observational errors, primarily due to evolutionary effects within the main sequence. For this reason, single-parameter relations used to predictM andR for single stars are limited to an accuracy of some ±15% inM and ±50% inR, basically independent of the number and accuracy of the data used to establish the relations. Two-parameter calibrations are discussed (Sect. 4) which can eventually reduce these errors to & 5% in bothM andR. At this level, abundance effects become significant and presumably account for the residual scatter.Comparison of the data with stellar evolution models is the topic of Sect. 5. Characteristic features of the data which are crucial in such work are emphasized, rather than attempts to prove the validity of any particular set of models. Already fromM andR alone, some significant constraints can be derived (Fig. 4). When bothM, R, andT _e are known, the initial helium abundanceY can be estimated if the metal-abundance parameter Z is assumed or determined. Studies in which binaries with accurate values ofM, R, and Z are fit by models calculated for the precise observed masses, and withY and mixing length constrained to solar values, provide the most stringent tests of the models. Probing further model refinements such as convective overshooting requires full use of the potential of the data. For example, models may yield general main-sequence limits which are consistent with the observations, but still be unable to fit any single system to the precision of the data. Conditions for critical, informative tests are discussed. Tidal effects in binaries are briefly discussed in Sect. 6. As tidal forces are extremely sensitive to the dimensions and internal structure of the stars, the present sample is well suited for such studies. Recent success in matching computed and observed apsidal-motion parameters for early-type binaries is mentioned. Finally, main priorities for future work are outlined.     

Role of the Rossiter-McLaughlin effect in the study of close binaries

This paper is devoted to binary stars belonging to the class of eclipsing-variable systems.Photometric and spectroscopic analysis of eclipses allows us to determine geometric parameters of the orbit and physical characteristics of stellar components as well as inclinations of stellar equators to the orbital plane. Estimations of inclinations can be obtained from measurement of the Rossiter-McLaughlin effect, which is discussed using examples of some eccentric binaries with an anomalous apsidal effect. Our task is to find the complete spectrum of solutions of the equation of apsidal motion, depending on the inclinations of the polar axes of the components to the orbital one for these systems, based on their individual spectroscopic and photometric observational data. The matrix of solutions allows us to select those pairs of polar inclinations that provide agreement with the observational apsidal period.     

Spectroscopic study of the optical counterpart to the fast X-ray transient IGR J17544-2619 based on observations at the 1.5-m RTT-150 telescope

We present the results of our long-term photometric and spectroscopic observations at the Russian–Turkish RTT-150 telescope for the optical counterpart to one of the best-known sources, representatives of the class of fast X-ray transients, IGR J17544-2619. Based on our optical data, we have determined for the first time the orbital and physical parameters of the binary system by the methods of Doppler spectroscopy.We have calculated theoretical spectra of the optical counterpart by applying non- LTE corrections for selected lines and obtained the parameters of the stellar atmosphere (T _eff = 33 000 K, log g = 3.85, R = 9.5 R _⊙, and M = 23 M _⊙). The latter suggest that the optical star is not a supergiant as has been thought previously.     

Properties of Wolf-Rayet stars in eclipsing binary systems

The results of investigations of a number of eclipsing Wolf-Rayet binaries are presented. The ‘core’ radiuses, the ‘core’ temperatures and masses of WR stars in the eclipsing WR+OB binary systems V 444 Cyg, CX Cep, CQ Cep, and CV Ser are obtained (see Table I). The results obtained from the light curves analysis of the V 444 Cyg in the range λλ2460 Å-3.5μ give strong evidence for the Beals (1944) model of WR phenomenon. The chromospheric-coronal effects in the WN5 extended atmosphere are not observed up to a distance ofr?20R _⊙. In the Hertzsprung—Russell diagram all the WR stars lie on the left side from the main sequence between the main sequence and the sequence of uniform helium stars (see Figure 9). Their locations are close to those of the helium remnants formed as a result of mass exchange in massive close binary systems. The period variations in the systems V 444 Cyg and CQ Cep have been discovered and a reliable value of the mass loss rateM=10^?5 M _⊙yr^?1 is obtained, for the two WR stars. The results of the photometric and spectroscopic investigations of the WR stars with low mass companions (post X-ray binary stage?) are presented too (see Table II). The masses of the companions are (1–2)M _⊙, their optical luminosity is ～10³⁶, erg s^?1 which implies that these companions cannot be the normal stars. It is possible that these companions are neutron stars accreting from the stellar wind of the WR stars. Low values of the X-ray luminosities of such WR stars with low mass companions imply that the accretion of matter in such systems is distinct from the accretion process in classical X-ray binary systems. It is noted also that the parameters of low massive companions coupled with WR stars are close to those of helium stars.     

Photometric elements and apsidal motion of the eclipsing binary NSV 18773

The photometric elements of the eclipsing binary NSV 18773 (HD 99898) have been determined for the first time by analyzing its V-and I-band light curves from the ASAS-2 and ASAS-3 catalogs. Based on these elements and using other published spectroscopic and photometric data, we constructed a consistent system of geometrical and physical parameters for the system that consists of two stars (M ₁ = 20M _⊙, Sp₁=B0V, R ₁ = 5.0R _⊙ and M ₂ = 14M _⊙, Sp₂ = B1V, R ₂ = 6.5R _⊙) in elliptical orbits (P = 5 _. ^d 049, e = 0.365, a = 40.1R _⊙). The distance to the system is d = 3.3 kpc, the interstellar extinction is A _V = 2 _. ^m 0, and the age is t = 2.8 × 10⁶ yr. NSV 18773 is a visual binary with components V _A = 9 _. ^m 9 and V _B = 10 _. ^m 3 separated by 0 _. ^" 8. The third light (L ₃ = 0.61) that we found by analyzing the light curves shows that the eclipsing binary is the system’s fainter component B. We confirmed the rapid apsidal motion of the star detected by Otero and Wils (2006) and refined its observed period: U _obs = 150 ± 6 yr. Our photometric elements and physical parameters allowed the apsidal parameter $\bar k_2^{obs} = 0.0135(14)$ , which reflects the density distribution along the radii of the component stars, to be determined. Within the error limits, the derived parameter agrees with its theoretically expected value, $\bar k_2^{th} = 0.0119(8)$ , from current evolutionary models of stars of the corresponding masses and ages.     

Graphite grain formation in the atmospheres of R coronae borealis stars

The results of calculations of graphite grain formation in the atmospheres of R CrB stars are given. The parameters for the models wereM=1M _⊙,M _bol=?6 mag. The effective temperature ranged from 5300K to 8300K. The chemical composition corresponded to the hydrogen-deficient carbon rich mixture:X=0,Y=0.9,Z _c=0.1. The results obtained show the existence of a critical mass loss rate which is ranged fromM _*≈10^?6 M _⊙yr^?1 forT _eff=5300 K toM _*≈10^?5 M _⊙ yr^?1 forT _eff=8300 K. As soon as the rate of mass loss exceedsM _* by 3–5 times the degree of condensation of carbon changes from 0 to 0.7. The finite radii of grains are about from 0.01 μm to 0.6 μm depending on the density near the condensation point, the velocity of matter outflow, and the stellar effective temperature. The duration of grain growth should amount to some dozens of days. It is supposed that the most probable explanation of dust-shell formation around R CrB stars is graphite condensation behind a shock wave arising from nonlinear stellar pulsation.     

Effects of stellar wind,dynamical friction,and star mergers on the dynamical evolution of multiple stars

The dynamical evolution of small stellar groups composed of N=6 components was numerically simulated within the framework of a gravitational N-body problem. The effects of stellar mass loss in the form of stellar wind, dynamical friction against the interstellar medium, and star mergers on the dynamical evolution of the groups were investigated. A comparison with a purely gravitational N-body problem was made. The state distributions at the time of 300 initial system crossing times were analyzed. The parameters of the forming binary and stable triple systems as well as the escaping single and binary stars were studied. The star-merger and dynamical-friction effects are more pronounced in close systems, while the stellar wind effects are more pronounced in wide systems. Star-mergers and stellar wind slow down the dynamical evolution. These factors cause the mean and median semimajor axes of the final binaries as well as the semimajor axes of the internal and external binaries in stable triple systems to increase. Star mergers and dynamical friction in close systems decrease the fraction of binary systems with highly eccentric orbits and the mean component mass ratios for the final binaries and the internal and external binaries in stable triple systems. Star mergers and dynamical friction in close systems increase the fraction of stable triple systems with prograde motions. Dynamical friction in close systems can both increase and decrease the mean velocities of the escaping single stars, depending on the density of the interstellar medium and the mean velocity of the stars in the system.     

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.     

Accuracy of abundance determinations by weighting functions method of curves of growth

Study of curve of growth abundance determinations is performed for the method of weighting functions. We evaluated abundances of iron for 15 stars changing their effective temperatures and gravities within the mean error of determination of these parameters, i.e.θ _eff ± 0.05θ _eff and logg ± 0.5. Wide range of spectral types of stars taken into account, permits analysis of accuracy of abundance determinations for large part of H-R diagram. The results prove, that in spite of using the most accurate method of curve of growth analysis, the accuracy of abundance determination (due to errors inT _eff and log g) is sometimes very small (in extreme cases not better than 0.4 dex) and the value of accuracy is connected with of spectral types of stars.     

Atmospheric Parameters and Abundances of sdB Stars

We summarize recent results of quantitative spectral analyses using NLTE and metal line-blanketed LTE model atmospheres. Temperatures and gravities derived for hundreds of sdB stars are now available and allow us to investigate systematic uncertainties of T _e ff, log g scales and to test the theory of stellar evolution and pulsations. Surface abundance patterns of about two dozen sdB stars are surprisingly homogenous. In particular the iron abundance is almost solar for most sdBs. We highlight one iron-deficient and three super metal-rich sdBs, a challenge to diffusion theory. sdB stars are slowly rotating stars unless they are in close binary systems, which is hard to understand if the sdB stars were formed in merger events. The only exception is the pulsator PG 1605+072 rotating at vsin i= 39 km/s. Signatures of stellar winds from sdB stars may have been found.     

FUSE Observations of EC14026 Stars

We present a FUSE abundance analysis of EC14026 stars. We compare the abundances of heavy elements in the atmospheres of EC14026 stars to non-pulsating stars with similar atmospheric parameters, and investigate whether weak stellar winds could explain the coexistence of variable and non-variable sdB stars in the log g – T_eff diagram. We also present preliminary results on time-dependent diffusion calculations of iron in presence of radiative levitation and mass loss, and show how weak stellar winds can affect the diffusive equilibrium between gravitational settling and radiative support.     

Determining the character of motion in quiet and active galaxies with a satellite companion

A galaxy model with a satellite companion is used to study the character of motion for stars moving in the x ‐y plane. It is observed that a large part of the phase plane is covered by chaotic orbits. The percentage of chaotic orbits increases when the galaxy has a dense nucleus of massM_n. The presence of the dense nucleus also increases the stellar velocities near the center of the galaxy. For small values of the distance R between the two bodies, low energy stars display a chaotic region near the centre of the galaxy, when the dense nucleus is present, while for larger values of R the motion in active galaxies is regular for low energy stars. Our results suggest that in galaxies with a satellite companion, the chaotic character of motion is not only a result of galactic interaction but also a result caused by the dense nucleus. Theoretical arguments are used to support the numerical outcomes. We follow the evolution of the galaxy, as mass is transported adiabatically from the disk to the nucleus. Our numerical results are in satisfactory agreement with observational data from M51‐type binary galaxies (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stellar adiabatic mass loss model and applications

Roche-lobe overflow and common envelope evolution are very important in binary evolution, which is believed to be the main evolutionary channel to hot subdwarf stars. The details of these processes are difficult to model, but adiabatic expansion provides an excellent approximation to the structure of a donor star undergoing dynamical time scale mass transfer. We can use this model to study the responses of stars of various masses and evolutionary stages as potential donor stars, with the urgent goal of obtaining more accurate stability criteria for dynamical mass transfer in binary population synthesis studies. As examples, we describe here several models with the initial masses equal to 1 M _⊙ and 10 M _⊙, and identify potential limitations to the use of our results for giant-branch stars.     

Study of the apsidal motion in the eclipsing binary MZ Lac

A series of highly accurate photoelectric observations of the eclipsing binary MZ Lac was obtained with a 48-cm AZT-14 reflector at the Tien-Shan High-Altitude Station of the Sternberg Astronomical Institute from 1985 to 2004 to study its apsidal motion. We constructed a consistent system of physical and geometrical parameters of the components and the binary’s orbit: we determined their masses (M₁ = 1.50M_⊙, M₂ = 1.29M_⊙), radii (R₁ = 1.86R_⊙, R₂ = 1.35R_⊙), luminosities (L₁ = 0.79L_⊙, L₂ = 0.45L_⊙), surface gravities (log_g1 = 4.06, log_g2 = 4.27), age (t = 1.9 × 10⁹ yr), and the distance to the binary (d = 510 pc). The binary exhibits apsidal motion with the period U_obs = 480 ± 40 yr, while its theoretically expected value is U_th = 450 ± 40 yr. Spectroscopic studies of MZ Lac and calculations of the absolute parameters of the components are required to test our conclusions.     

Period distributions in pulsating and binary stars

We analyze the hypothesis of quantization in bands for the angular momenta of binary systems and for the maount of actionA _c in stable and pulsating stars. This parameter isA _c=Mv _eff R _eff, where the effective velocity corresponds to the kinetic energy in the stellar interior and the effective radius corresponds to the potential energyGM ²/R _eff. Analogous parameters can be defined for a pulsating star withm=M where is the rate of the massm participating in the oscillation to the total massM andv _osc,R _osc the effective velocity and oscillation radius.From an elementary dimensional analysis one has thetA _c (energy x time) (period)^1/3 independently ifA _c corresponds to the angular momentum in a binary system, or to the oscillation in a pulsating star or the inner energy and its time-scaleP _eff in a stable star.From evolving stellar models one has that P _effP _eff(solar)1.22 hr a near-invariant for the Main Sequence and for the range of masses 0.6M <M<1.6M .With this one can give scalesn _k=kn ₁ withk integers andn ₁=(P/P ₁)^1/3 withP ₁=P _eff1.22 hr. In these scales proportional toA _c, one sees that the periods in binary and pulsating stars are clustered in discrete unitsn ₁,n ₂,n ₃, etc.This can be seen in pulsating Scuti, Cephei, RR Lyrae, W Virginis, Cephei, semi-regular variables, and Miras and in binary stars as cataclysmic binaries, W Ursa Majoris, Algols, and Lyrae with the corresponding subgroups in all these materials. Phase functions (n _k) in RR Lyrae and Cephei are also associated with discrete levelsn _k.the suggested scenario is that the potential energies and the amounts of actionE _p(t), A_c(t) are indeed time-dependent, but the stars remain more time in determinated most proble states. The Main Sequence itself is an example of this. These most probable states in binary systems, or pulsating or stable stars, must be associated with velocities sub-multiplesc/ _F , given by the velocity of light and the fine structure constant.Additional tests for such a hypothesis are suggested when the sufficient amount of observational data are available. They can made with oscillation velocities in pulsating stars and velocity differences of pairs of galaxies.     

Photometric elements,apsidal motion,and a third body in the eclipsing binary HP Aur

Highly accurate W BV R photometric measurements of the eclipsing binary HP Aur were performed in 2002–2003 with the 48-cm AZT-14 reflector at the Tien-Shan High-Altitude Observatory to determine the rate of apsidal motion. A consistent system of physical and geometrical parameters of the components and the binary as a whole has been constructed for the first time by analyzing these new measurements together with other published data: we determined their radii (R₁ = 1.05R_⊙, R₂ = 0.82R_⊙) and luminosities (L₁ = 1.10L_⊙, L₂ = 0.46L_⊙), spectral types (G2V + G8V) and surface gravities (log g₁ = 4.38, log g₂ = 4.51), age (t = 9.5 × 10⁹ yr), and the distance to the binary (d = 197 pc). We detected an ultraviolet excess in the spectra of both components, \(\Delta (W - B) \simeq - 0\mathop .\limits^m 25\), that is probably attributable to a metal deficiency in the atmospheres of these stars. In this system of two solar-type stars, we found a third body with the mass M₃ sin i ₃ ³ = 0.17_M⊙ that revolved with the period P₃ = 13.7 yr around the eclipsing binary in a highly eccentric elliptical orbit: e₃ = 0.70 and A₃ sin i₃ ? 7 AU. The orbit of the eclipsing binary itself was shown to be also elliptical, but with a low eccentricity (e = 0.0025(5)), while apsidal motion with a period U_obs > 80 yr was observed at a theoretically expected period U_th ≈ 92 yr. At least 20 to 30 more years of photoelectric measurements of this star will be required to reliably determine U_obs.     

Probability distribution of terrestrial planets in habitable zones around host stars

With more and more exoplanets being detected, it is paid closer attention to whether there are lives outside solar system. We try to obtain habitable zones and the probability distribution of terrestrial planets in habitable zones around host stars. Using Eggleton’s code, we calculate the evolution of stars with masses less than 4.00 M _⊙. We also use the fitting formulae of stellar luminosity and radius, the boundary flux of habitable zones, the distribution of semimajor axis and mass of planets and the initial mass function of stars. We obtain the luminosity and radius of stars with masses from 0.08 to 4.00 M _⊙, and calculate the habitable zones of host stars, affected by stellar effective temperature. We achieve the probability distribution of terrestrial planets in habitable zones around host stars. We also calculate that the number of terrestrial planets in habitable zones of host stars is 45.5 billion, and the number of terrestrial planets in habitable zones around K type stars is the most, in the Milky Way.