Motion Around The Triangular Equilibrium Points Of The Restricted Three-Body Problem Under Angular Velocity Variation

We study numerically the asymmetric periodic orbits which emanate from the triangular equilibrium points of the restricted three-body problem under the assumption that the angular velocity ω varies and for the Sun–Jupiter mass distribution. The symmetric periodic orbits emanating from the collinear Lagrangian point L ₃, which are related to them, are also examined. The analytic determination of the initial conditions of the long- and short-period Trojan families around the equilibrium points, is given. The corresponding families were examined, for a combination of the mass ratio and the angular velocity (case of equal eigenfrequencies), and also for the critical value ω = 2

The distribution law of stellarp-modes

We show that the overall densityg() of asymptotic acoustic frequencies of a star obeys a Weyl lawg() ^D–1, whereD is the dimensionality of the oscillating stellar configuration. For realistic stars with a finite non-zero surface sound speed,D is equal to the actual dimensionality of the star,D=3. For formal models with a vanishing sound velocity at the surface, heuristic arguments lead to a dimensionality parameterD=4.5. The empirical frequencies of Eddington's standard model are found to be consistent with the latter distribution, with reasonable agreement already occurring in the low-frequency range > _i 2× fundamental radial mode. We argue that real stars obey this 3.5-power law in some finite frequency interval _i << _f , _f being a very high frequency critically depending on the surface sound velocity, while the full asymptotic law, withD=3, holds for > _f .     

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.     

Macro-and micro-turbulent filter functions for weak lines in stellar atmospheres

Following a similar discussion given earlier for the solar case (De Jager, 1972) we compute in this paper spectral line profiles for the spatial wavelengths in which a stellar motion field can be decomposed, and thereafter the macro-and micro-turbulent filter functionsf _M(k) andf (k), where is the optical scale height andf ²(k) dk the fraction of the energy of the turbulent motions between wavenumbersk andk+dk of the spectrum of turbulence that contributes to either kind of turbulence. If micro-and macro-turbulent velocity components are known for a certain star, and if the spectrum of turbulence is sharp enough, the ratiof _M/f would enable one to derive the average size of the turbulent elements in the star's atmosphere. The computations apply to weak lines in idealized stellar atmospheres, and refer to two cases: isotropic turbulence, and radial pulsations. These filters can be suitably used in a diagnostic method for the analysis of the motion field in the solar and stellar atmospheres. Some examples of applications to stars of very different kinds are given.     

Apsidal motion in the eclipsing binary system of BW AQR

The first photoelectric light curve of the eclipsing binary system BW Aqr (F71V+F81V;P=6^d.7;V=10 ^m .31), discovered by Miss Leavitt at the beginning of the century, was obtained. The photometric elements were detemined. The components of this system are considerably evolved stars: the age of the system is about 2×10⁹ yr. It follows from the photometric data that the secondary component should have greater mass than the primary one The zero-age spectral classes of components were F2V and F1V. The system has an elliptical orbit with the eccentricitye=0.18. The angular rate of the apsidal motion (_obs = 0.070 deg yr^–1) and the corresponding value of the apsidal parameterk ₂=0.0090 (the relativistic term included) were found. The derived valuek ₂ exceeds by more than a factor of 2 the theoretical coefficient obtained from the modern theory of internal structure of stars with moderate masses .     

Invariant properties of families of moon-to-earth trajectories

Analytical techniques are employed to demonstrate certain invariant properties of families of moon-to-earth trajectories. The analytical expressions which demonstrate these properties have been derived from an earlier analytical solution of the restricted three-body problem which was developed by the method of matched asymptotic expansions. These expressions are given explicitly to orderµ ^1/2 where is the dimensionless mass of the moon. It is also shown that the inclusion of higher order corrections does not affect the nature of the invariant properties but only increases the accuracy of the analytic expressions.The results are compared with the work of Hoelker, Braud, and Herring who first discovered invariant properties of earth-to-moon trajectories by exact numerical integration of the equations of motion. (Similar properties for moon-to-earth trajectories follow from the principle of reflection). In each instance the analytical expressions result in properties which are equivalent, to orderµ ^1/2, with those found by numerical integration. Some quantitative comparisons are presented which show the analytical expressions to be quite accurate for calculating particular geometrical characteristics.

Nomenclature

Orbital Elements near the Moon energy - angular momentum - semi-major axis - eccentricity - inclination - argument of node - argument of pericynthion Orbital Elements near the Earth h _e energy - l _e angular momentum - i inclination - argument of node - argument of perigee - t _f time of flight Other symbols parameters used in matehing - U a function of the energy near the earth - a function of the angular momentum near the earth - r _p perigee radius - perincynthion radius - radius at node near moon - true anomaly of node near moon - initial angle between node near moon and earth-moon line - a function ofU, , andi - earth phase angle - dimensionless mass of the moon - U ₀, U₁ U=U ₀+U ₁ - i ₀, i_1/2, i₁ i=i ₀+µ ^1/2 i _1/2+µ i ₁ - ₀, _1/2, ₁ = ₀+µ ^1/2 i _1/2+µ i ₁ - _p longitude of vertex line - _n latitude of vertex line - R _o ,S _o ,N _o functions ofU ₀ and - a function ofU ₀, and      

Tidal torque constants as a critical test for an evolution of close eclipsing binaries

Individual tidal torque λ ₂,E ₂ and apsidal-motion k ₂ constants were calculated for 112 close eclipsing binaries (CEBs) with Detached components belonging to the Main Sequence (DMS-type) from the catalogue by Svechnikov and Perevozkina (Catalogue of orbital elements, masses and luminosities of variable stars of DMS-type and some results of its statistical treatment, Ural State University Press, Yekaterinburg, pp. 1–5, 1999) and for 95 detached binaries taken from the catalogue by Torres et al. (Astron. Astrophys. Rev. 18:67, 2010) on the base of theoretical evolutionary stellar models including tidal torque constants by Claret (Astron. Astrophys. 424:919, 2004). A method of the inversion of model track grid into isochrones was formulated as a complex interpolation procedure for DMS-binaries data. Sets of isochrones were computed in k ₂–M, k ₂–R, λ ₂–M, λ ₂–R, E ₂–M, and E ₂–R planes. Calculated tidal torque constants allow to test stellar structure theory by comparing observed and estimated values of apsidal motion period and analyzing the correlation between timescales of synchronization, circularization, magnetic braking, as well as nuclear burning of DMS-components.     

Third and fourth order resonances in Hamiltonian systems

The stability of the origin of an autonomous Hamiltonian system is investigated when the system possesses a third or fourth-order resonance.H ₂, the quadratic part ofH isH ₂=ⁿ _{i=1 i} J _i and the resonance condition is ⁿ _i=1 k _{i i} where thek 0,i = 1, 2, ...,n are the natural or fundamental frequencies. It is shown that the only case in which the origin can be unstable is ifk _i0,i=1,2,...,n. The condition for instability is then given in terms of the coefficients of the higher order terms in the Hamiltonian. The transfer of energy between modes is also investigated when a near-resonant condition exists.     

Angular momentum transport in pre-main-sequence stars of intermediate mass

Pre-Main-Sequence stars with masses between 2 and 5 M (Herbig Ae/Be stars) have radiative subphotospheric envelopes. However, they possess strong stellar winds and show definite signs of activity which could be linked to surface magnetic field. Therefore, they must lose angular momentum at a significant rate.We investigate the effect of such angular momentum losses on the internal structure of these stars, and on the distribution of angular velocity inside them. This paper presents a preliminary analysis guided by an analogy with laboratory and geophysical fluids. We propose that the friction exerted at the stellar surface by the angular momentum losses produces a mixed layer below the surface, separated from the unperturbed interior by an interface. Using scaling laws established by experimental studies of sheared stratified fluids, we discuss a simplified model for the evolution of the mixed layer.Although this model is still too preliminary to allow quantitative predictions, we show that for a reasonable choice of parameters, the mixed layer penetrates into the stellar interior on a time-scale of 10⁶ years, comparable to the Kelvin time-scale for the Herbig Ae/Be stars.     

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.     

Stellar analogs of solar magnetic activity

The techniques and principal results of observational studies of stellar activity are summarized. Both chromospheric and coronal emission clearly track surface magnetic field properties, but it is not well known how the detailed relation between the emission and surface magnetic fields varies with spectral type. For lower Main-Sequence stars of the same spectral type, there is clear evidence of a close relationship between mean activity level and rotation period P _rot. There is also less definitive evidence for a similar dependence on convective overturn time _c , such that activity depends on the single parameter Ro = P _rot/ _c . For single stars, stellar rotation, and magnetic activity both decline smoothly with age. This implies a feedback between angular momentum loss rate and activity level. Temporal variations in mean stellar activity level mimic the solar cycle only for old stars like the Sun, being much more irregular for younger stars. The characteristic timescale of the variations (the cycle period) appears to depend on Ro for old stars, but shows no clear dependence on either rotation rate or spectral type for younger stars. Further data on mean activity and its variation for a large number of lower Main-Sequence stars should contribute significantly to our understanding of the causes of stellar magnetic activity.     

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.     

The angular momentum loss of low-mass close binary systems

The detailed evolution of low-mass main-sequence stars (M < 1M ) with a compact companion is studied. For angular momentum loss associated with magnetic braking it is found that about 10^–11–10^–12 M yr^–1 in stellar wind loss would be required. This wind is 10²–10³ times stronger than the solar wind, so we believe here magnetic stellar wind is insufficient. It is well known that there is mass outflow in low-mass close binary systems. We believe here that these outflows are centrifugal driven winds from the outer parts of the accretion disks. The winds extract angular momentum from these systems and therefore drive secular evolution. Disk winds are preferred to winds from the secondary, because of the lower disk surface gravity.     

Existence of the continuations in theN-body problem

The method of obtaining the estimates of the maximalt-interval ( _–, ₊) on which the solution of theN-body problem exists and which is such that some fixed mutual distance (e. g. ₁₂) exceeds some fixed non-negative lower bound, for allt contained in ( _–, ₊), is considered. For given masses and initial data, the increasing sequences of the numbers _k , each of which provides the estimate ₊ > _k , are constructed. It appears that if ₊ = +, then .     

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.     

Models of pulsating low-massive yellow supergiants

The nonlinear self-excited oscillations of the envelopes of low-massive highly luminous stars are described. The parameters for these models wereM=0.8M ,M _bol=–5.5, –5.84 mag,T _eff=4500, 5000, 5500 K. The oscillations have been found to consist of the standing wave pulsation near the envelope bottom and running waves in outer layers. The ratio of the standing wave frequency _s to the average frequency of the running waves _r increases with the stellar luminosity: _s / _r =1.7 whenM _bol=–5.5 mag and _s / _r =2.4 whenM _bol=–5.84 mag. The frequency of oscillations near the photosphere is found to be in close agreement with the critical frequency for running waves. Mass loss from these stars is caused by shocks. It has been shown that agreement between FG Sge's period change observed during the last decade and the period-luminosity relation for double shell stars takes place when FG Sge's luminosity isM _bol=–5.96 mag.     

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 .     

Some applications of Jacobi dynamics to the stellar evolution

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 _cM ^k withk1.87 for the M stars,k5/3 for the KGF stars, andk1.83 for the A and earlier stars, representing very simples conditions for the other dynamical parameters. For instancek5/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 M25M .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 (MM _–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.     

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.