Equilibrium Structures of Differentially Rotating Primary Components of Binary Stars

In this paper a method is proposed for computing the equilibrium structures and various other observable physical parameters of the primary components of stars in binary systems assuming that the primary is more massive than the secondary and is rotating differentially about its axis. Kippenhahn and Thomas averaging approach (1970) is used in a manner earlier used by Mohan, Saxena and Agarwal (1990) to incorporate the rotational and tidal effects in the equations of stellar structure. Explicit expressions for the distortional terms appearing in the stellar structure equations have been obtained by assuming a general law of differential rotation of the typeω² = b ₀+b ₁ s ²+b ₂ s ⁴, where ω is the angular velocity of rotation of a fluid element in the star at a distance s from the axis of rotation, and b ₀, b ₁, b ₂ are suitably chosen numerical constants. The expressions incorporate the effects of differential rotation and tidal distortions up to second order terms. The use of the proposed method has been illustrated by applying it to obtain the structures and observable parameters of certain differentially rotating primary components of the binary stars assuming the primary components to have polytropic structures. This revised version was published online in July 2006 with corrections to the Cover Date.     

Equilibrium structure of stars obeying a generalized differential rotation law

In the present paper we have considered the problem of determining the equilibrium structure of differentially rotating stars in which the angular velocity of rotation varies both along the axis of rotation and in directions perpendicular to it. For this purpose, a generalized law of differential rotation of the type ² =b ₀+b ₁ s ²+b ₂ s ⁴+b ₃ z ²+b ₄ z ⁴+b ₅ z ² s ² (here is a nondimensional measure of the angular velocity of a fluid element distants from the axis of rotation andz from the plane through the centre of the star perpendicular to the axis of rotation, andb's are suitably chosen parameters) has been used. Whereas Kippenhahn and Thomas averaging approach has been used to incorporate the rotational effects in the stellar structure equations, Kopal's results on Roche equipotentials have been used to obtain the explicit form of the stellar structure equations, which incorporate the rotational effects up to second order of smallness in the distortion parameters. The method has been used to compute the equilibrium structure of certain differentially rotating polytropes. Certain differentially rotating polytropes. Certain differentially rotating models of the Sun have also been computed by using this approach.     

Equilibrium Structures of Differentially Rotating and Tidally Distorted White Dwarf Models of Stars

In this paper we present a method for computing the equilibrium structures and various physical parameters of a primary component of the binary system assuming that the primary is more massive than the secondary and is rotating differentially according to the law of the w² = b₀ + b₁ × s² + b₂ × s⁴, w being the angular velocity of rotation of a fluid element distant s from the axis of rotation and b₀, b₁, b₂ suitably chosen numerical constants. This method utilizes the averaging approach of Kippenhahn and Thomas (1997) and the concept of Roche equipotentials in a manner earlier used by Mohan et al. (1997) to incorporate the effects of rotation and tidal distortions on the equilibrium structures of certain rotationally and tidally distorted stellar models. The use of the method has been illustrated by applying it to obtain the structures and some observable parameters of certain differentially rotating and tidally distorted binary systems whose primary component is assumed to be a white dwarf star.     

Kinematics of the galaxy from OB stars with proper motions from the Gaia DR1 catalogue

We consider two samples of OB stars with different distance scales that we have studied previously. The first and second samples consist of massive spectroscopic binaries with photometric distances and distances determined from interstellar calcium lines, respectively. The OB stars are located at heliocentric distances up to 7 kpc. We have identified them with the Gaia DR1 catalogue. Using the proper motions taken from the Gaia DR1 catalogue is shown to reduce the random errors in the Galactic rotation parameters compared to the previously known results. By analyzing the proper motions and parallaxes of 208 OB stars from the Gaia DR1 catalogue with a relative parallax error of less than 200%, we have found the following kinematic parameters: (U, V)_⊙ = (8.67, 6.63)± (0.88, 0.98) km s^?1, Ω₀ = 27.35 ± 0.77 km s^?1 kpc^?1, Ω′₀ = ?4.13 ± 0.13 km s^?1 kpc^?2, and Ω″₀ = 0.672 ± 0.070 km s^?1 kpc^?3, the Oort constants are A = ?16.53 ± 0.52 km s^?1 kpc^?1 and B = 10.82 ± 0.93 km s^?1 kpc^?1, and the linear circular rotation velocity of the local standard of rest around the Galactic rotation axis is V ₀ = 219 ± 8 km s^?1 for the adopted R ₀ = 8.0 ± 0.2 kpc. Based on the same stars, we have derived the rotation parameters only from their line-of-sight velocities. By comparing the estimated values of Ω′₀, we have found the distance scale factor for the Gaia DR1 catalogue to be close to unity: 0.96. Based on 238 OB stars of the combined sample with photometric distances for the stars of the first sample and distances in the calcium distance scale for the stars of the second sample, line-of-sight velocities, and proper motions from the Gaia DR1 catalogue, we have found the following kinematic parameters: (U, V, W)_⊙ = (8.19, 9.28, 8.79)± (0.74, 0.92, 0.74) km s^?1, Ω₀ = 31.53 ± 0.54 km s^?1 kpc^?1, Ω′₀ = ?4.44 ± 0.12 km s^?1 kpc^?2, and Ω″₀ = 0.706 ± 0.100 km s^?1 kpc^?3; here, A = ?17.77 ± 0.46 km s^?1 kpc^?1, B = 13.76 ± 0.71 km s^?1 kpc^?1, and V ₀ = 252 ± 8 km s^?1.     

Equilibrium models of differentially-rotating polytropic cylinders

The equilibrium structure of differentially rotating polytropic cylinders is determined numerically. We setn=3 and use a quadratic function for the law of differential rotation. We construct different models by varying the angular velocity at the axis and the ratio of the angular velocity at the surface to the angular velocity at the axis. By taking a decreasing function for the rotation law we are able to treat models with an angular velocity at the axis greater than the break-up velocity of uniformly rotating cylinders. We also determine whether a Richardson-like criterion for stability is violated in the models.     

Eigenfrequencies of small adiabatic barotropic modes of oscillations of rotationally and tidally-distorted stars

In this paper a method is proposed for computing the eigenfrequencies of small adiabatic barotropic modes of oscillations of rotationally and tidally-distorted stars. The method utilizes Kippenhahn and Thomas approach and concepts of Roche equipotentials to incorporate up to second-order the effects of rotation and tidal distortion terms on the eigenfrequencies. The proposed method has also been used to compute the eigenfrequencies of certain barotropic modes of oscillation of some rotationally and tidally distorted models of 10M , and 2.5M Main-Sequence stars.     

Relativistic polytropic spheres in general relativity

Spherical-symmetric massive structures with a relativistic polytropic equation of state have been studied in detail. Thevarious static parameters have been calculated under all possible variations of ₀ (=P ₀/E ₀0 pertaining to the central value) and the polytropic indexn. In particular, the values of the static parameters have been identified under the extreme relativistic conditions and for the most bound structures. The structures have been used to model neutron stars and the various neutron parameters have been calculated. Radial pulsation and rotation of these configurations have also been discussed.     

Hydrodynamic equations of differentially rotating stellar systems and the velocity ellipsoid

Anisotropic hydrodynamic equations for differentially rotating collisionless stellar systems are derived. These equations can describe the evolution of the systems in a time span longer than their rotation periods.As a by-product of derivation of hydrodynamic equations, the well-known relation that the ratio of the principal axes of the velocity ellipse in a differentially rotating stellar disk is [B/(B-A)]^1/2 is re-found if the system is in a purely circular rotation, whereA andB are the Oort's constants. In addition, we find a systematic mean motion superposed on a purely circular differential rotation makes the directions of axes of the velocity ellipse deviate from the radial and the transverse direction. The observed deviation of directions of axes in our neighbourhood in the Galaxy can be explained if in the mean motion superposed on a purely circular differential rotatin the gas of stars near us is compressed in the radial direction or rarefied in the transverse directions, with irregularities of the order of 5 km/sec in amplitude of velocity and 1 kpc in size. These magnitudes of irregularities agree with those actually observed or with those anticipated from other theoretical considerations.     

Testing the Distance Scale of the Gaia TGAS Catalogue by the Kinematic Method

We have studied the simultaneous and separate solutions of the basic kinematic equations obtained using the stellar velocities calculated on the basis of data from the Gaia TGAS and RAVE5 catalogues. By comparing the values of Ω'₀ found by separately analyzing only the line-of-sight velocities of stars and only their proper motions, we have determined the distance scale correction factor p to be close to unity, 0.97 ± 0.04. Based on the proper motions of stars from the Gaia TGAS catalogue with relative trigonometric parallax errors less than 10% (they are at a mean distance of 226 pc), we have found the components of the group velocity vector for the sample stars relative to the Sun (U, V,W)_⊙ = (9.28, 20.35, 7.36) ± (0.05, 0.07, 0.05) km s^?1, the angular velocity of Galactic rotation Ω⁰ = 27.24 ± 0.30 km s?1 kpc?1, and its first derivative Ω'₀ = ?3.77 ± 0.06 km s^?1 kpc^?2; here, the circular rotation velocity of the Sun around the Galactic center is V₀ = 218 ± 6 km s^?1 kpc (for the adopted distance R₀ = 8.0 ± 0.2 kpc), while the Oort constants are A = 15.07 ± 0.25 km s^?1 kpc^?1 and B = ?12.17 ± 0.39 km s^?1 kpc^?1, p = 0.98 ± 0.08. The kinematics of Gaia TGAS stars with parallax errors more than 10% has been studied by invoking the distances from a paper by Astraatmadja and Bailer-Jones that were corrected for the Lutz–Kelker bias. We show that the second derivative of the angular velocity of Galactic rotation Ω'₀ = 0.864 ± 0.021 km s^?1 kpc^?3 is well determined from stars at a mean distance of 537 pc. On the whole, we have found that the distances of stars from the Gaia TGAS catalogue calculated using their trigonometric parallaxes do not require any additional correction factor.     

Kinematics of the galaxy from Cepheids with proper motions from the Gaia DR1 catalogue

A sample of classical Cepheids with known distances and line-of-sight velocities has been supplemented with proper motions from the Gaia DR1 catalogue. Based on the velocities of 260 stars, we have found the components of the peculiar solar velocity vector (U, V, W)_⊙ = (7.90, 11.73, 7.39) ± (0.65, 0.77, 0.62) km s^?1 and the following parameters of the Galactic rotation curve: Ω₀ = 28.84 ± 0.33 km s^?1 kpc^?1, Ω′₀ = ?4.05 ± 0.10 km s^?1 kpc^?2, and Ω″₀ = 0.805 ± 0.067 km s^?1 kpc^?3 for the adopted solar Galactocentric distance R ₀ = 8 kpc; the linear rotation velocity of the local standard of rest is V ₀ = 231 ± 6 km s^?1.     

Kinematic analysis of solar-neighborhood stars based on RAVE4 data

We consider stars with radial velocities, proper motions, and distance estimates from the RAVE4 catalogue. Based on a sample of more than 145 000 stars at distances r < 0.5 kpc, we have found the following kinematic parameters: \({\left( {U,{\kern 1pt} V,{\kern 1pt} W} \right)_ \odot }\) = (9.12, 20.80, 7.66) ± (0.10, 0.10, 0.08) km s^?1, Ω₀ = 28.71 ± 0.63 km s^?1 kpc^?1, and Ω₀ ^’ = ?4.28 ± 0.11 km s^?1 kpc^?2. This gives the linear rotation velocity V ₀ = 230 ± 12 km s^?1 (for the adopted R ₀ = 8.0 ± 0.4 kpc) and the Oort constants A = 17.12 ± 0.45 km s^?1 kpc^?1 and B = ?11.60 ± 0.77 km s^?1 kpc^?1. The 2D velocity distributions in the UV, UW, and VW planes have been constructed using a local sample, r < 0.25 kpc, consisting of ~47 000 stars. A difference of the UV velocity distribution from the previously known ones constructed from a smaller amount of data has been revealed. It lies in the fact that our distribution has an extremely enhanced branch near the Wolf 630 peak. A previously unknown peak at (U, V) = (?96, ?10) km s^?1 and a separate new feature in the Wolf 630 stream, with the coordinates of its center being (U, V) = (30, ?40) km s^?1, have been detected.     

Correlation and variability of the HeI 4471 and MgII 4481 absorption lines in Be stars: A possible diagnostic tool for nonradial pulsations

In a total of 181 photographic spectrograms, obtained at the Manuel Foster Observatory of the Catholic University in Santiago (Chile), the equivalent widths of the absorption lines HeI 4471 and MgII 4481 have been determined for 10 selected southern Be stars and 4 normal B-type stars. The mean equivalent widths of each star are compared to published mean values from largers samples of B stars. On the average, Be stars coincide with normal B-type stars in their equivalent width. However, a striking difference between rapidly and slowly rotating stars in their line ratios (HeI/MgII) was found, which is probably an artifact due to a systematic misclassification of spectral types of rapidly rotating B stars (v sini>200 km s^–1). Six of the ten Be stars show evidence for variability in HeI and/or MgII with time-scales shorter than a few days and amplitudes up to a factor 2 in equivalent width. Line variability occurs in the entire range of projected rotation velocity (70 km s^–1<v sini<350 km s^–1), but seems to be restricted to spectral types B2-B4. In addition, variations at longer time-scales (years) were observed in two cases. Three of the variable stars (HR 4074, HR 4537, and Ara) present correlated variations in both lines, one ( Cen) anti-correlated variability. We interpret our results in terms of nonradial pulsations and stress the importance to study the HeI/MgII lines and their variability in a larger sample of Be, Bn, and normal B stars, including archival material in order to follow-up line variability for several decades in the past.     

Differential rotation in F stars

Differential rotation can be detected in single line profiles of stars rotating more rapidly than about v sin i = 10km s-1 with the Fourier transform technique. This allows to search for differential rotation in large samples to look for correlations between differential rotation and other stellar parameters. I analyze the fraction of differentially rotating stars as a function of color, rotation, and activity in a large sample of F-type stars. Color and rotation exhibit a correlation with differential rotation in the sense that more stars are rotating differentially in the cooler, less rapidly rotating stars. Effects of rotation and color, however, cannot be disentangled in the underlying sample. No trend with activity is found. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Pre-Main-Sequence evolution of rotating low-mass stars

The evolutionary behaviour of rotating low-mass stars in the mass range 0.2 and 0.9M has been investigated during the pre-Main-Sequence phase. The angular momentum is conserved locally in radiative regions and totally in convective regions, according to a predetermined angular velocity distribution depending on the structure of the star. As the stars contract toward the zero-age Main Sequence, they spin up under the assumption that the angular momentum is conserved during the evolution of the stars. When the stars have differential rotations, their inner regions rotate faster than the outer regions. The effective temperatures and luminosities of rotating low-mass stars are obtained lower than those of non-rotating stars. They have lower central temperature and density values compared to those of non-rotating stars.     

Rapidly rotating polytropes in the post-Newtonian approximation to general relativity

The study of uniformly polytropes with axial symmetry is extended to include all rotational terms of order ⁴, where is the angular velocity, consistently within the first post-Newtonian approximation to general relativity. The equilibrium structure is determined by treating the effects of rotation and post-Newtonian gravitation as independent perturbations on the classical polytropic structure. The perturbation effects are characterized by a rotation parameter = ²/2G _c and a relativity parameter, =p _c / _c C ² , wherep _c and _c are the central pressure and density respectively. The solution to the structural problem is obtained by following Chandrasekhar's series expansion technique and is complete to the post-Newtonian rotation terms of order ². The critical rotation parameterv _c , which characterizes the configuration with maximum uniform rotation, is accurately evaluated as a function of . Numerical values for all the structural parameters needed to determine the equilibrium configurations are presented for polytropes with indicesn=1, 1.5, 2, 2 5, 3, and 3.5.     

Two-and three-layer models of Uranus

We present simple two-layer models of Uranus with rocky core and polytropic envelope satisfying exactly the observed mass, radius and the gravitational moments. The models show that the value of the fourth order zonal harmonic isJ ₄ –38×10^–6, whileJ ₆ 10^–6. More elaborate threelayer models fail to satisfy the observational constraints of the ice/rock ratio and/or of the rotation period. We conclude that three-layer models with uniform chemical composition in each layer may be too restrictive. More realistic models should account for variable chemical composition within each layer.     

H12CO+ outflow and a rotating H13CO+ disk around L1551 IRS-5

High-velocity resolution (V=0.07 km s^–1) H¹²CO⁺ (J=1–0) and H¹³CO⁺ (J=1–0) observations have been carried out towards L1551 by use of the Metsähovi 14-m radio telescope. The observations reveal a bipolar H¹²CO⁺ outflow from the pre-Main-Sequence star IRS-5 which is centred on a flattened, 3 long H¹³CO⁺ cloud clump. This disk-like cloud structure has a velocity gradient ranging from 6.56 km s^–1 in the SE to 7.06 km s^–1 in the NW. It is noteworthy that the direction of the H¹²CO⁺ ion outflow is oriented E-W, and not along the NE-SW axis of the more extended CO outflow. In the disk area the H¹²CO⁺ spectra show to distinct velocity components. The right-hand H¹²CO⁺ velocity component agrees with the velocity of the H¹³CO⁺ disk. The left-hand H¹²CO⁺ component seems to belong to the outflow and the dense lobe material. The H¹²CO⁺ isovelocity contour map indicates that the dense lobe material is rotating (V _rot 0.6 km s^–1) in the same sense as the H¹³CO⁺ disk. This supports hydromagnetic outflow models.Paper presented at the 11th European Regional Astronomical Meetings of the IAU on New Windows to the Universe, held 3–8 July, 1989, Tenerife, Canary Islands, Spain     

Rapidly rotating supermassive stars in the first post-Newtonian approximation to general relativity

The structure and stability of rapidly uniformly rotating supermassive stars is investigated using the full post-Newtonian equations of hydrodynamics. The standard model of a supermassive star, a polytrope of index three, is adopted. All rotation terms up to and including those of order ⁴, where is the angular velocity, are retained. The effects of rotation and post-Newtonian gravitation on the classical configuration are explicitly evaluated and shown to be very small. The dynamical stability of the model is treated by using the binding energy approach. The most massive objects are found to be dynamically unstable when =1/c ².p _c / _c 2.2 × 10^–3, wherep _c and _c are the central pressure and density, respectively. Hence, the higher-order terms considered in this analysis do not appreciably alter the previously known stability limits.The maximum mass that can be stabilized by uniform rotation in the hydrogen-burning phase is found to be 2.9×10⁶ M , whereM is the solar mass. The corresponding nuclear-generated luminosity of 6×10⁴⁴ erg/sec^–1 is too small for the model to be applicable to the quasi-stellar objects. The maximum kinetic energy of a uniformly rotating supermassive star is found to be 3×10^–5 Mc ², whereM is the mass of the star. Masses in excess of 10¹⁰ M are required if an adequate store of kinetic energy is to be made available to a pulsar like QSO. However such large masses have rotation periods in excess of 100 yr and thus could not account for any short term periodic variability. It is concluded then that the uniformly rotating supermassive star does not provide a suitable base for a model of a QSO.     

Mass loss from coronae and its effect upon stellar rotation

The acoustic energy-generation rate from the convective zone was calculated for various models. Results show that chromosphere and corona can be expected around stars with temperature lower than 8000K at the main sequence, and lower than 6500K at logg=2.When a star is rotating rapidly, mass loss from its corona is large, and can be an effective mechanism of braking the stellar rotation. If this mechanism is effective, we can explain the slow rotation of stars later than F2 to be the result of the loss of the angular momentum through a stellar wind that is effective in their main sequence phase. Stars with massM>1.5M lose mass through a stellar wind during their contraction phase. The mass-loss rate is larger than the solar value because of the larger energy input into the chromosphere-corona system and because of the smaller gravitational potential at the surface. T Tauri stars may be the observational counterparts for such stars. As the duration of contraction phase is very short (less than 10⁷ years), the braking mechanism works only in the presence of a strong magnetic field (Ap) or in the presence of a companion (Am).Presented at the Trieste Colloquium on Mass Loss from Stars, September 12–16, 1968.     

Kinematics of the Galaxy from OB Stars with Data from the Gaia DR2 Catalogue

We have selected and analyzed a sample of OB stars with known line-of-sight velocities determined through ground-based observations and with trigonometric parallaxes and propermotions from the Gaia DR2 catalogue. Some of the stars in our sample have distance estimates made from calcium lines. A direct comparison with the trigonometric distance scale has shown that the calcium distance scale should be reduced by 13%. The following parameters of the Galactic rotation curve have been determined from 495 OB stars with relative parallax errors less than 30%: (U, V,W)_⊙ = (8.16, 11.19, 8.55)± (0.48, 0.56, 0.48) km s^?1, Ω₀ = 28.92 ± 0.39 km s^?1 kpc^?1, Ω'₀ = ?4.087 ± 0.083 km s^?1 kpc^?2, and Ω″ ₀ = 0.703 ± 0.067 km s?1 kpc^?3, where the circular velocity of the local standard of rest is V0 = 231 ± 5 km s^?1 (for the adopted R₀ = 8.0 ± 0.15 kpc). The parameters of the Galactic spiral density wave have been found from the series of radial, V_R, residual tangential, ΔV_circ, and vertical, W, velocities of OB stars by applying a periodogram analysis. The amplitudes of the radial, tangential, and vertical velocity perturbations are f_R = 7.1± 0.3 km s^?1, f_θ = 6.5 ± 0.4 km s^?1, and f_W = 4.8± 0.8 km s^?1, respectively; the perturbation wavelengths are λ_R = 3.3 ± 0.1 kpc, λ_θ = 2.3 ± 0.2 kpc, and λ_W = 2.6 ± 0.5 kpc; and the Sun’s radial phase in the spiral density wave is (χ_⊙)_R = ?135? ± 5?, (χ_⊙)θ = ?123? ± 8?, and (χ_⊙)_W = ?132? ± 21? for the adopted four-armed spiral pattern.