WR 140 (=V1687 Cyg): Infrared photometry, 2001–2010

We present the results of our infrared observations of WR 140 (=V1687 Cyg) in 2001–2010. Analysis of the observations has shown that the J brightness at maximum increased near the periastron by about 0 ^m .3; the M brightness increased by ∼2 ^m in less than 50 days. The minimum J brightness and the minimum L and M brightnesses were observed 550–600 and 1300–1400 days after the maximum, respectively. The JHKLM brightness minimum was observed in the range of orbital phases 0.7–0.9. The parameters of the primary O5 component of the binary have been estimated to be the following: R(O5) ≈ 24.7R _⊙, L(O5) ≈ 8 × 10⁵ L _⊙, and M _bol(O5) ≈ −10 ^m . At the infrared brightness minimum, T _g ∼ 820–880 K, R _g ≈ 2.6 × 10⁵ R _⊙, the optical depth of the shell at 3.5 μm is ∼5.3 × 10⁻⁶, and its mass is ≈1.4 × 10⁻⁸ M _⊙. At the maximum, the corresponding parameters are ∼1300 K, 8.6 × 10⁴ R _⊙, ∼2 × 10⁻⁴, and ∼6 × 10⁻⁸ M _⊙; the mean rate of dust inflow (condensation) into the dust structure is ∼3.3 × 10⁻⁸ M _⊙ yr⁻¹. The mean escape velocity of the shell from the heating source is ∼10³ km s⁻¹ and the mean dispersal rate of the shell is ∼1.1 × 10⁻⁸ M _⊙ yr⁻¹.     

A Model Atmosphere Analysis of the F6V Star π3 Ori

Model atmosphere analysis, based on Kurucz models has been applied to study the F6V star π³ Ori (=BS1543=HD30652). The following values of the effective temperature, surface gravity and microturbulence velocity were obtained: = 6270±200 K, log g = 3.80.2, ξ_t =3.5±0.5 km/s. The abundances of 10 elements were determined. The resulting element abundances for the π³ Ori were found to be about three times lower with respect to the Sun. From evolutionary calculations we derived a mass, radius and luminosity for π³ Ori of M =1.3 M_⊙, R =2.38 R_⊙, L =7.9 L_⊙. Hence this star should be classified F6IV instead of F6 V. This revised version was published online in July 2006 with corrections to the Cover Date.     

Pulsational instability of yellow hypergiants

Instability of population I (X = 0.7, Z = 0.02) massive stars against radial oscillations during the post-main-sequence gravitational contraction of the helium core is investigated. Initial stellar masses are in the range 65M _⊙ ≤ M _ZAMS ≤ 90M _⊙. In hydrodynamic computations of self-exciting stellar oscillations we assumed that energy transfer in the envelope of the pulsating star is due to radiative heat conduction and convection. The convective heat transfer was treated in the framework of the theory of time-dependent turbulent convection. During evolutionary expansion of outer layers after hydrogen exhaustion in the stellar core the star is shown to be unstable against radial oscillations while its effective temperature is T _eff > 6700 K for M _ZAMS = 65M _⊙ and T _eff > 7200 K for M _ZAMS = 90M _⊙. Pulsational instability is due to the κ-mechanism in helium ionization zones and at lower effective temperature oscillations decay because of significantly increasing convection. The upper limit of the period of radial pulsations on this stage of evolution does not exceed ≈200 day. Radial oscillations of the hypergiant resume during evolutionary contraction of outer layers when the effective temperature is T _eff > 7300 K for M _ZAMS = 65M _⊙ and T _eff > 7600 K for M _ZAMS = 90M _⊙. Initially radial oscillations are due to instability of the first overtone and transition to fundamental mode pulsations takes place at higher effective temperatures (T _eff > 7700 K for M _ZAMS = 65M _⊙ and T _eff > 8200 K for M _ZAMS = 90M _⊙). The upper limit of the period of radial oscillations of evolving blueward yellow hypergiants does not exceed ≈130 day. Thus, yellow hypergiants are stable against radial stellar pulsations during the major part of their evolutionary stage.     

Instability of LBV stars against radial oscillations

Hydrodynamic calculations of nonlinear radial oscillations of LBV stars with effective temperatures 1.5 × 10⁴ K ⩽ T _eff ⩽ 3 × 10⁴ K and luminosities 1.2 × 10⁶ L _⊙ ⩽ L ⩽ 1.9 × 10⁶ L _⊙ have been performed. Models for the evolutionary sequences of Population I stars (X = 0.7, Z = 0.02) with initial masses 70M _⊙ ⩽ M _ZAMS ⩽ 90M _⊙ at the initial helium burning stage have been used as the initial conditions. The radial oscillations develop on a dynamical time scale and are nonlinear traveling waves propagating from the core boundary to the stellar surface. The amplitude of the velocity variations for the outer layers is several hundred km s⁻¹, while the bolometric magnitude variations are within ΔM _bol ⩽ 0_· ^m 2. The onset of oscillations is not related to the κ-mechanism and is attributable to the instability of a self-gravitating envelope gas whose adiabatic index is close to its critical value of Γ₁ = 4/3 due to the dominant contribution of radiation in the internal energy and pressure. The interval of magnitude variation periods (6 days ≤ II ≤ 31 days) encompasses all currently available estimates of the microvariability periods for LBV stars, suggesting that this type of nonstationarity is pulsational in origin.     

Absolute dimensions of the close binary systems V453 Monocerotis and V523 Cassiopeiae

We present multi-colour CCD observations of the low-temperature contact binaries, V453 Mon and V523 Cas. Their light curves are modelled to determine a new set of stellar and orbital parameters. Analysis of mid-eclipse times yields a new linear ephemeris for both systems. A period decrease (dP/dt=2.3×10⁻⁷ days/yr) in V453 Mon is discovered. V523 Cas, however, is detected to show a period increase (dP/dt=9.8×10⁻⁸ days/yr) because of the mass transfer of a rate of 1.1×10⁻⁷ M_⊙ yr⁻¹, from a less massive donor. Using these findings we can determine the physical parameters of the components of V523 Cas to be M ₁=0.76 (3)M _⊙, M ₂=0.39 (2)M _⊙, R ₁=0.74 (2)R _⊙, R ₂=0.55 (2)R _⊙, L ₁=0.19 (3)L _⊙, L ₂=0.14 (3)L _⊙, and the distance of system as 46(9) pc.     

Aq Psc – Analysis of New Light Curves

New BV light curves of the A-type W UMa star AQ Psc (P = 0.476d) have been observed and are described. A few times of minimum light are obtained and the ephemeris is improved. The light curves are analyzed for the binary parameters with a light curve synthesis method. Combining the results with Lu and Rucinski’s spectroscopic mass ratio we determined the masses and radii of the components: M ₁ = 1.69M _⊙, M ₂ = 0.38M _⊙, R ₁ = 1.77R _⊙, and R ₂ = 0.89R _⊙.     

The evolution of BD+60°2522: with mass loss and overshooting

The ionizing star BD+60°2522 is known as the central star of Bubble Nebulae NGC 7635—wind-blown bubble created by the interaction of the stellar wind of BD+60°2522 (O6.5 IIIef, V=8.7 mag, mass loss rate 10^−5.76 M _⊙/year) with the ambient interstellar medium. From the evolutionary calculations for the star with mass loss and overshooting, we find that the initial mass of the star is 60M _⊙, its present age is 2.5×10⁶ years, and the present mass is 45M _⊙.     

Astrometric study of the relative motion of three stars with possible invisible companions based on homogeneous series obtained at Pulkovo with a 26-inch refractor

We investigate the relative motion of three stars, ADS 7446, 9346, and 9701, based on long-term observations with the Pulkovo 26-inch refractor. The relative motion of all three stars shows a perturbation that could be produced by the gravitational influence of an invisible companion. For ADS 7446, we have determined the orbit of the photocenter with a period of 7.9 yr; the mass of the companion is more than 0.4M _⊙. For ADS 9346, we have determined the radial velocities of the components: −14.60 km s⁻¹ for A and −13.94 km s⁻¹ for B. For ADS 9346 and 9701, we have determined the dynamical parallaxes, 24 and 20 mas, respectively, which are larger than those in the Hipparcos catalog by 5 mas, and calculated the orbits by the apparent motion parameter (AMP) method. The new orbit of ADS 9346 is: a = 5″.2, P = 2035 yr, and e = 0.46 at the system’s mass M = 2.5M _⊙. The new orbits of ADS 9701 are: (a = 2″.9, P = 829 yr, e = 0.54, M = 4.3M _⊙) and (a = 3″.8, P = 1157 yr, e = 0.53, M = 5.0M _⊙).     

Radial pulsations of massive main-sequence stars

The evolution of Population I stars (X = 0.7, Z = 0.02) with initial masses 40M _⊙ ≤ M _ZAMS ≤ 120M _⊙ until core hydrogen exhaustion has been computed. Models of evolutionary sequences have been used as the initial conditions in solving the equations of radiation hydrodynamics that describe the spherically symmetric motion of a self-gravitating gas. Stars with initial masses M _ZAMS ≥ 50M _⊙ are shown to become unstable against radial oscillations during the main-sequence evolution. The instability growth rate and the limit-cycle oscillation amplitude increase as the star evolves and as its initial mass increases. The pulsational instability is attributable to the iron Z-bump κ mechanism (T ∼ 2 × 10⁵ K). Convection that transfers from 20 to 50% of the total energy flux and, thus, reduces the efficiency of the κ mechanism emerges in the same layers. The periods of the radial oscillations of main-sequence stars lie within the range from 0.09 to 8 days. The boundaries of the instability region of radial pulsations in the Hertzsprung-Russell diagram have been determined and observational criteria for revealing pulsating variable main-sequence stars have been proposed.     

Two new LBV candidates in the M 33 galaxy

We present two new luminous blue variable (LBV) candidate stars discovered in the M33 galaxy. We identified these stars as massive star candidates at the final stages of evolution, presumably with a notable interstellar extinction. The candidates were selected from the Massey et al. catalog based on the following criteria: emission in H _α , V<18^./m 5 and 0_.^m 35 < (B - V) < 1_.^m 2. The spectra of both stars reveal a broad and strong H _α emission with extended wings (770 and 1000 kms⁻¹). Based on the spectra we estimated the main parameters of the stars. Object N45901 has a bolometric luminosity log(L/L_⊙) = 6.0–6.2 with the value of interstellar extinction A _V = 2.3 ± 0.1. The temperature of the star’s photosphere is estimated as T⋆ ∼ 13000–15000 K, its probable mass on the Zero Age Main Sequence is M∼ 60–80 M_⊙. The infrared excess in N 45901 corresponds to the emission of warm dust with the temperature Twarm ∼ 1000 K, and amounts to 0.1%of the bolometric luminosity. A comparison of stellar magnitude estimates from different catalogs points to the probable variability of the object N45901. Bolometric luminosity of the second object, N125093, is log(L/L_⊙) = 6.3 − 6.6, the value of interstellar extinction is A _V = 2.75 ± 0.15. We estimate its photosphere’s temperature as T⋆∼ 13000–16000K, the initial mass as M ∼ 90–120 M_⊙. The infrared excess in N125093 amounts to 5–6% of the bolometric luminosity. Its spectral energy distribution reveals two thermal components with the temperatures T_warm ∼ 1000K and T_cold ∼ 480 K. The [Ca II] λλ7291, 7323 lines, observed in LBV-like stars Var A and N93351 in M33 are also present in the spectrum of N 125093. These lines indicate relatively recent gas eruptions and dust activity linked with them. High bolometric luminosity of these stars and broad H α emissions allow classifying the studied objects as LBV candidates.     

On the Formation of Massive Primordial Stars

We investigate the formation by accretion of massive primordial protostars in the range 10 to 300 M _⊙. The high accretion rate used in the models (M _⊙ = 4.4 x 10^-3 M_⊙ yr^-1) causes the structure and evolution to differ significantly from those of both present-day protostars and primordial zero-age main sequence stars. The stellar surface is not visible throughout most of the main accretion phase, since a photosphere is formed in the in falling envelope. Significant nuclear burning does not take place until a protostellar mass of about 80 M _⊙. As the interior luminosity approaches the Eddington luminosity, the protostellar radius rapidly expands owing to the radiation pressure. Eventually, a final swelling occurs when the stellar mass reaches about 300 M _⊙. This expansion is likely to signal the end of the main accretion phase, thus setting an upper limit to the protostellar mass formed in these conditions. This revised version was published online in September 2006 with corrections to the Cover Date.     

Absolute parameters of the close binary BW Boo with a superasynchronous A2m primary component

Based on two high-dispersion spectra of the close binary BW Boo, we have detected lines of the secondary component whose contribution to the combined spectrum does not exceed 2%. We have determined the rotation velocities of the components and spectroscopic orbital elements. Numerous lines of neutral and ionized iron have been used to determine the effective temperature and surface gravity for the primary component. The photometric light curves for this binary have been solved for the first time. Its primary component is an A2Vm star with a mass of 2 ± 0.1M _⊙ and a radius of 1.9 ± 0.4R _⊙. Its rotation velocity is 2 km s⁻¹, which is a factor of 18 lower than the pseudo-synchronous velocity for this component. The G6 secondary component, a T Tau star, has a rotation velocity of 17 km s⁻¹, amass of 1.1M _⊙, and a radius of 1 R _⊙. The age of the binary has been estimated to be 10⁷ yr.     

Photometric study of the W UMa eclipsing binaries VW LMi and BX Dra

New ephemeris and the absolute parameters—masses, radii and luminosities—of the contact systems VW LMi and BX Dra have been obtained, by means of the analysis of the minima data available in the literature (for the determination of the ephemeris) and combining the previously published spectroscopic information and the results of the Wilson-Devinney method using photometric data (for the determination of the absolute parameters). The VW LMi O−C analysis confirms the multiplicity of the system detected previously from the spectroscopic data. Masses of the VW LMi contact system primary and secondary components are 1.67 ± 0.02M _⊙ and 0.70 ± 0.02M _⊙, respectively. The corresponding radii are 1.709 ± 0.007R _⊙ and 1.208 ± 0.006R _⊙, respectively. For the BX Dra contact system the masses are 2.19 ± 0.13M _⊙ and 0.63 ± 0.06M _⊙, and the radii, 2.13 ± 0.04R _⊙ and 1.26 ± 0.03R _⊙, for the primary and secondary, respectively. In both cases, the estimated luminosities seem to be slightly greater that the values derived from the Hipparcos distances.     

Spectroscopic Orbit of the Eclipsing Binary AI Draconis

New radial velocity measurements of the Algol-type eclipsing binary AI Dra, based on Reticon observations, are presented. The velocity measures themselves are based on fitting theoretical profiles, generated by a physical model of the binary, to the observed cross-correlation function (ccf). Such profiles match this function very well, much better in fact that Gaussian profiles which are generally used. Measuring the ccf's with Gaussian profiles yields following results: m_p sin³ i=2.55± 0.05m_⊙, m_s sin³ i = 1.14 ± 0.03m_⊙, (a_p + a_s) sin i=7.34 ±0.05R_⊙, and m_p/m_s =2.23± 0.05. Where as measuring the ccf's with theoretical profiles yields a mass ratio of 2.33 and following results: m_p sin³ i=2.84± 0.05m_⊙, m_s sin³ i=1.22 ± 0.03m_⊙, (a_p +a_s) sin i=7.56± 0.05R_⊙. The system comprising a semi-detached configuration. From the solution of a previously published light curved and combining it's results with the spectroscopic orbit, one can lead to the following physical parameters: m_p =2.99m_⊙, m_s =1.28m_⊙, > T_p < =9600 K, > T_s < =5400 K, > R_p < =2.35R_⊙, > R_s < =2.12R_⊙. The system comprising an AO primary and a secondary of G2 spectral type. This revised version was published online in July 2006 with corrections to the Cover Date.     

Mass-luminosity relation for massive stars

A catalog of massive (⩾10 M _⊙) stars in binary and multiple systems with well-known masses and luminosities has been compiled. The catalog is analyzed using a theoretical mass-luminosity relation. This relation allows both normal main-sequence stars and stars with peculiarities: with clear manifestations of mass transfer, mass accretion, and axial rotation, to be identified. Least-squares fitting of the observational data in the range of stellar masses 10M _⊙ ⩽ M ≲ 50 M _⊙ yields the relation L ∼ M ^2.76. An erratum to this article is available at .     

The bolometric light curve of SN 1993J and the nature of its progenitor

We have constructed the bolometric light curve of SN 1993J based on UBVRI(JHK) photometric data obtained from various sources and assumingA _V = 0 and a distance modulus of 27.6. Effective temperatures and photosphere radius at various times have been obtained from detailed blackbody fits. The bolometric light curve shows two maxima. The short rise time to the second maximum, and the luminosities at the minimum and the second maximum are used to constrain the properties of the progenitor star. The total mass of the hydrogen envelope M_H, in the star is found to be ≲ 0.2 M_⊙ at the time of explosion, and the explosion ejected about 0.05 M_⊙ of Ni⁵⁶. Thin hydrogen envelope combined with a sufficient presupernova luminosity suggest that the exploding star was in a binary with a probable period range of 5yr ≤P _orb ≤ 11yr.     

Changes in the Sun’s mass and gravitational constant estimated using modern observations of planets and spacecraft

More than 635 thousand positional observations of planets and spacecraft of various types (mostly radiotechnical ones, 1961–2010) were used to estimate possible changes in the gravitational constant, Sun’s mass, and semi-major axes of planetary orbits, as well as the associated value of the astronomical unit. The observations were analyzed based on the EPM2010 ephemerides constructed at the Institute of Applied Astronomy (Russian Academy of Sciences) in a post-Newtonian approximation as a result of simultanious numerical integration of the equations of motion of nine major planets, the Sun, the Moon, asteroids, and trans-Neptunian objects. The heliocentric gravitational constant GM _⊙ was found to vary with a rate of (GṀ _⊙/GM _⊙ = (−5.0 ± 4.1)) × 10⁻¹⁴ per year (at the 3σ level). The positive secular changes in the semimajor axes ȧ _i /a _i were found for Mercury, Venus, Mars, Jupiter, and Saturn provided by high-precision observations. These changes also correspond to the decrease in the heliocentric gravitational constant. The changing of GM _⊙, itself is probably caused by the loss of the mass M _⊙ of the Sun due to its radiation and solar wind; these effects are partly compensated by the material falling onto the Sun. Allowing for the maximum bounds on the possible change in the Sun’s mass M _⊙, it has been found from the change obtained in GM _⊙ that the annual change Ġ/G of the gravitational constant G falls within the interval −4.2 × 10⁻¹⁴ < ȧ/G < +7.5 × 10⁻¹⁴ with a 95% probability. The astronomical unit (AU) is connected by its definition only with the heliocentric gravitational constant. The decrease of GM _⊙ obtained in this paper should correspond to a secular decrease in the AU. It is shown, however, that the modern level of accuracy does not allow us to determine a change in the AU. The attained posibility of determining changes in GM _⊙ using high-accuracy observations encourages us to have a relation between GM _⊙ and the AU fixed for a certain moment in time, since it is inconvenient to have a time-dependent length for the AU.     

Extrapolating SMBH correlations down the mass scale: the case for IMBHs in globular clusters

Empirical evidence for both stellar mass black holes (M _•<10²  M _⊙) and supermassive black holes (SMBHs, M _•>10⁵  M _⊙) is well established. Moreover, every galaxy with a bulge appears to host a SMBH, whose mass is correlated with the bulge mass, and even more strongly with the central stellar velocity dispersion σ _c , the M _•–σ relation. On the other hand, evidence for “intermediate-mass” black holes (IMBHs, with masses in the range 100–10⁵ M _⊙) is relatively sparse, with only a few mass measurements reported in globular clusters (GCs), dwarf galaxies and low-mass AGNs. We explore the question of whether globular clusters extend the M _•–σ relationship for galaxies to lower black hole masses and find that available data for globular clusters are consistent with the extrapolation of this relationship. We use this extrapolated M _•–σ relationship to predict the putative black hole masses of those globular clusters where existence of central IMBH was proposed. We discuss how globular clusters can be used as a constraint on theories making specific predictions for the low-mass end of the M _•–σ relation.     

Orbital elements of the RS CVn eclipsing binary,SV Camelopardalis

The extensiveUBV observations of SV Camelopardalis by Patkos (1982) have been analysed to derive the orbital elements of the system. The data were corrected for the effect of third body (Sarma, Sarma & Abhyankar 1985) and for the ‘RS CVn’ distortion wave (Sarma, Vivekanandarao & Sarma 1988). The cleaned data were used to obtain a preliminary solution by a modified version of Wellmann method (Sarma & Abhyankar 1979) from which we concluded that the primary eclipse is a transit. The final orbital elements of SV Cam were obtained by the modified version (Sarma 1988; Sarmaet al. 1987) of WINK program by Wood (1972). The colour and median brightness variation are discussed. From the spectroscopic mass functionf(m) = 0.118 M_⊙ (Hiltner 1953), the absolute dimensions of the components are found to be 0.826 M_bd & 0.592 M_⊙ and 1.236 R_⊙ & 0.778 R_⊙ for the primary and secondary components, respectively. The age of the binary system is estimated to be 6.0 ± 1.0 × 10⁸ years     

Radial Velocity and Light Curve Analysis of the Eclipsing Binary Nn Vir

The eclipsing binary NN Vir is a short period system showing an EW-type light curve. Photometric observations of NN Vir were done by Gomez–Ferrellad and Garcia–Melendo (1997) at Esteve Duran Observatory. The first spectroscopic observations of this system were obtained by Rucinski and Lu (1999). The radial velocity and light curves analysis was made with the latest version of the Wilson program (1998), and the geometric and physical elements of the system are derived. From the simultaneous solutions of the system, we determined the masses and radii of the components: 1.89 M _⊙ and 1.65 R _⊙ for the primary component; 0.93 M _⊙ and 1.23 R _⊙ for the secondary component. We estimated effective temperatures of 7030 K for the primary and 6977 K for the secondary component.