The equilibrium and stability of magnetopolytropes

The theory of the oscillations of axisymmetric gaseous configurations with a prevalent magnetic field is presented. The virial tensor method is used to obtain the nine second harmonic modes of oscillations of the system. It is found that out of the nine modes, three are neutral, four are non-radial, and two are coupled. For the Prendergast spherical model it is found that one of the coupled modes is radial and the other non-radial. Both the radial and the non-radial modes obtained in this case agree with the corresponding formulae obtained byChandrasekhar andLimber (1954) andWoltjer (1962).The equilibrium structure of gaseous polytropes with toroidal magnetic fields is also investigated in detail for values of the polytropic indexn=1, 1.5, 2, 3 and 3.5. For this model the components of the moment of intertia and potential energy tensors together with the non-zero components of the supermatrix potential are obtained. The final results in terms of the effect of weak toroidal magnetic fields on the characteristic frequencies of distorted polytropes are presented in the form of tables.     

Dynamical tides in close binary systems,I

The aim of the present paper will be to develop from the fundamental equations of hydrodynamics a theory of dynamical tides in close binary systems, the components of which are regarded to consist of heterogeneous viscous fluid, and to revolve around their common centre of gravity in eccentric orbits; moreover, the equatorial planes of their axial rotation and the orbital plane need not be co-planar, but all may be inclined to the invariable plane of the system of arbitrary amounts. The changes in the pressure or density invoked by time-dependent deformation will be regarded as adiabatic; but, in the equilibrium state, both the density and viscosity of the material of our components may be arbitrary functions of the radial distance.Following a brief exposition in Section 2 of the fundamental equations linearized to small oscillations — be these free or forced — in Section 3 we shall particularize them to describe spheroidal deformations; with due regard to all terms arising from viscosity. Section 4 will contain a specification of the boundary conditions to be imposed upon such oscillations; and in Section 5 we shall solve the problem of non-radial oscillations of self-gravitating inviscid configurations in terms of hypergeometric series. The remaining Sections 6–8 will be devoted to a discussion of the phenomena arising from viscosity: in particular, we shall solve in a closed form the problem of non-radial oscillations of incompressible viscous globes in the terms of Bessel functions. It will be shown that the effect of viscosity — like those of compressibility — tend to de-stabilize all non-radial oscillations of homogeneous configurations.At the other extreme, a similar treatment of a mass-point model — as well as of one exhibiting high but finite degree of central condensation — is being postponed for a subsequent communication.     

Vibrational stability of rotating stars

The aim of the present paper will be to detail the explicit form of the equations which govern first-order oscillations of fast-rotating globes of self-gravitating fluids; with due account taken of the effects arising from the centrifugal as well as Coriolis force. As such configurations oscillate in general about distorted figures of equilibrium, the equations governing them can be conveniently expressed in terms of the Clairaut coordinates, associated with distorted spheroidal figures, and introduced in our previous paper (Kopal, 1980) for this purpose.In Section 2 which follows a brief outline of our problem, the equilibrium properties of fast-rotating configurations or arbitrary structure will be formulated. In Section 3 we shall carry out a separation of the variables in the equations of motion, and reduce the partial differential equations of the problem to an equivalent system of ordinary differential equations, by an expansion of expressions for the velocity componentsU, V, W in terms of tesseral harmonicsY _n ^m (, ). The explicit form of such a system, including the effects of all tesseral harmonics of orders up tom=n=4, will be specified in Section 3 for configurations whose equilibrium form is a sphere; while in Section 4 this latter condition will be relaxed to allow for the equilibrium configuration to become a rotational spheroid.In the concluding Section 5 we shall convert the complex form of our equations of motion into real terms, amenable to a solution-analytical or numerical-in terms of real variables; and shall establish the boundary conditions necessary for a specification of the characteristic frequencies of oscillation.     

A dynamical equivalent to the equilateral libration points of the earth-moon system

Consider the Earth-Moon-particle system as a Restricted Three Body Problem. There are two equilateral libration points. In the actual world system, those points are no longer relative equilibrium points mainly due to the effect of the Sun and to the noncircular motion of the Moon around the Earth. In this paper we present the problem as a perturbation of the RTBP and we look for the dynamical equivalent of L _4,5. It turns out to be a quasiperiodic orbit. It is obtained for a simplified model but the procedure to obtain it is general and can be carried out with an additional computational effort.     

On the oscillations in Mercury's obliquity

Mercury's capture into the 3:2 spin–orbit resonance can be explained as a result of its chaotic orbital dynamics. One major objective of MESSENGER and BepiColombo spatial missions is to accurately measure Mercury's rotation and its obliquity in order to obtain constraints on internal structure of the planet. Analytical approaches at the first-order level using the Cassini state assumptions give the obliquity constant or quasi-constant. Which is the obliquity's dynamical behavior deriving from a complete spin–orbit motion of Mercury simultaneously integrated with planetary interactions? We have used our SONYR model (acronym of Spin–Orbit N-bodY Relativistic model) integrating the spin–orbit N-body problem applied to the Solar System (Sun and planets). For lack of current accurate observations or ephemerides of Mercury's rotation, and therefore for lack of valid initial conditions for a numerical integration, we have built an original method for finding the libration center of the spin–orbit system and, as a consequence, for avoiding arbitrary amplitudes in librations of the spin–orbit motion as well as in Mercury's obliquity. The method has been carried out in two cases: (1) the spin–orbit motion of Mercury in the 2-body problem case (Sun–Mercury) where an uniform precession of the Keplerian orbital plane is kinematically added at a fixed inclination (S_2K case), (2) the spin–orbit motion of Mercury in the N-body problem case (Sun and planets) (S_n case). We find that the remaining amplitude of the oscillations in the S_n case is one order of magnitude larger than in the S_2K case, namely 4 versus 0.4 arcseconds (peak-to-peak). The mean obliquity is also larger, namely 1.98 versus 1.80 arcminutes, for a difference of 10.8 arcseconds. These theoretical results are in a good agreement with recent radar observations but it is not excluded that it should be possible to push farther the convergence process by drawing nearer still more precisely to the libration center. We note that the dynamically driven spin precession, which occurs when the planetary interactions are included, is more complex than the purely kinematic case. Nevertheless, in such a N-body problem, we find that the 3:2 spin–orbit resonance is really combined to a synchronism where the spin and orbit poles on average precess at the same rate while the orbit inclination and the spin axis orientation on average decrease at the same rate. As a consequence and whether it would turn out that there exists an irreducible minimum of the oscillation amplitude, quasi-periodic oscillations found in Mercury's obliquity should be to geometrically understood as librations related to these synchronisms that both follow a Cassini state. Whatever the open question on the minimal amplitude in the obliquity's oscillations and in spite of the planetary interactions indirectly acting by the solar torque on Mercury's rotation, Mercury remains therefore in a stable equilibrium state that proceeds from a 2-body Cassini state.     

Periodic and quasi-periodic motions of a solar sail close to SL 1 in the Earth–Sun system

Solar sails are a proposed form of spacecraft propulsion using large membrane mirrors to propel a satellite taking advantage of the solar radiation pressure. To model the dynamics of a solar sail we have considered the Earth–Sun Restricted Three Body Problem including the Solar radiation pressure (RTBPS). This model has a 2D surface of equilibrium points parametrised by the two angles that define the sail orientation. In this paper we study the non-linear dynamics close to an equilibrium point, with special interest in the bounded motion. We focus on the region of equilibria close to SL ₁, a collinear equilibrium point that lies between the Earth and the Sun when the sail is perpendicular to the Sun–sail direction. For different fixed sail orientations we find families of planar, vertical and Halo-type orbits. We have also computed the centre manifold around different equilibria and used it to describe the quasi-periodic motion around them. We also show how the geometry of the phase space varies with the sail orientation. These kind of studies can be very useful for future mission applications.     

Radial velocities of pulsating subdwarf B stars: KPD 2109+4401 and PB 8783

High-speed spectroscopy of two pulsating subdwarf B stars, KPD 2109+4401 and PB 8783, is presented. Radial motions are detected with the same frequencies as reported from photometric observations and with amplitudes of ∼2 km s⁻¹ in two or more independent modes. These represent the first direct observations of surface motion arising from multimode non-radial oscillations in subdwarf B stars. In the case of the sdB+F binary PB 8783, the velocities of both components are resolved; high-frequency oscillations are found only in the sdB star and not the F star. There also appears to be evidence for mutual motion of the binary components. If confirmed, it implies that the F-type companion is ≳1.2 times more massive than the sdB star, while the amplitude of the F-star acceleration over 4 h would constrain the orbital period to lie between 0.5 and 3.2 d.     

Perturbation of the radial and non-radial oscillations of a star by a magnetic field

A generalization of the perturbation method is applied to the problem of the radial and the non-radial oscillations of a gaseous star which is distorted by a magnetic field. An expression is derived for the perturbation of the oscillation frequencies due to the presence of a weak magnetic field when the equilibrium configuration is a spheroid. The particular application to the homogeneous model with a purely poloidal field inside, due to a current distribution proportional to the distance from the axis of symmetry, and a dipole type field outside is considered in detail. The main result is that the magnetic field has a large and almost stabilizing effect on unstableg-modes, particularly on higher order modes. With the considered magnetic field the surface layers appear to have a large weight.     

On the vertical families of two-dimensional tori near the triangular points of the Bicircular problem

This paper focuses on some aspects of the motion of a small particle moving near the Lagrangian points of the Earth–Moon system. The model for the motion of the particle is the so-called bicircular problem (BCP), that includes the effect of Earth and Moon as in the spatial restricted three body problem (RTBP), plus the effect of the Sun as a periodic time-dependent perturbation of the RTBP. Due to this periodic forcing coming from the Sun, the Lagrangian points are no longer equilibrium solutions for the BCP. On the other hand, the BCP has three periodic orbits (with the same period as the forcing) that can be seen as the dynamical equivalent of the Lagrangian points. In this work, we first discuss some numerical methods for the accurate computation of quasi-periodic solutions, and then we apply them to the BCP to obtain families of 2-D tori in an extended neighbourhood of the Lagrangian points. These families start on the three periodic orbits mentioned above and they are continued in the vertical (z and ż) direction up to a high distance. These (Cantor) families can be seen as the continuation, into the BCP, of the Lyapunov family of periodic orbits of the Lagrangian points that goes in the (z, ż) direction. These results are used in a forthcoming work [9] to find regions where trajectories remain confined for a very long time. It is remarkable that these regions seem to persist in the real system. This revised version was published online in July 2006 with corrections to the Cover Date.     

Hydromagnetic oscillations of a self-gravitating fluid spheroid

The oscillations of a homogeneous, compressible, self gravitating fluid spheroid in static equilibrium with a poloidal magnetic field inside and a dipole type field outside are studied using the second order tensor virial equations. It is found that for small values of the eccentricity, the equilibrium model is dynamically stable provided the usual criterion of pulsative stability in the absence of a magnetic field (>4/3) is satisfied. The magnetic field removes the accidental degeneracy of the radial and the non-radial modes of oscillation which exists for =1.6 in the absence of a magnetic field.     

A generalized time-dependent local convection theory

The time-dependent convection theory given in [1] is generalized to include non-radial oscillations. Under some reasonable assumptions, dynamic equations are derived of the auto-correlation and cross-correlation functions of the velocity and temperature fluctuations. Steady and pulsating components are separated and solved for. Our steady component agrees with Vitense's formulae. The pulsating component contains a phase-lag, due to the inertia of the convective motion.     

Effect of Moon perturbation on the energy curves and equilibrium points in the Sun–Earth–Moon system

In this paper, we have considered that the Moon motion around the Earth is a source of a perturbation for the infinitesimal body motion in the Sun–Earth system. The perturbation effect is analyzed by using the Sun–Earth–Moon bi–circular model (BCM). We have determined the effect of this perturbation on the Lagrangian points and zero velocity curves. We have obtained the motion of infinitesimal body in the neighborhood of the equivalent equilibria of the triangular equilibrium points. Moreover, to know the nature of the trajectory, we have estimated the first order Lyapunov characteristic exponents of the trajectory emanating from the vicinity of the triangular equilibrium point in the proposed system. It is noticed that due to the generated perturbation by the Moon motion, the results are affected significantly, and the Jacobian constant is fluctuated periodically as the Moon is moving around the Earth. Finally, we emphasize that this model could be applicable to send either satellite or telescope for deep space exploration.     

Non-Radial Pulsations in Components of Symbiotic Stars

Photometric observations of symbiotic stars in the blue and in the red spectral regions make it possible to reveal non-radial oscillations both of the cool and of the hot components. Light variations of red giants in the symbiotic systems CI Cyg and AG Peg show several periods in the 10–80^d range, interpreted as p-mode pulsations. These modes are excited by a bright spot produced by radiation flux from the hot component. The spot moves on the red giant’s photosphere at a velocity close to the sound speed. During the active phase of the symbiotic star CH Cyg, at least 25 frequencies of oscillations in the 150–6000 s range of periods were found in the light of the white dwarf. Their features correspond to non-radial g-modes. In the frame of 2D gas dynamical non-adiabatic models, the interaction between gas flows and the accretion disk leads to formation of a system of shock waves propagating towards the compact object, which is one of possible mechanisms to excite non-radial pulsations of white dwarfs in symbiotic systems.     

The vector approach to the problem of physical libration of the Moon: the linearized problem

A flexible and informative vector approach to the problem of physical libration of the rigid Moon has been developed in which three Euler differential equations are supplemented by 12 kinematic ones. A linearized system of equations can be split into an even and odd systems with respect to the reflection in the plane of the lunar equator, and rotational oscillations of the Moon are presented by superposition of librations in longitude and latitude. The former is described by three equations and consists of unrestricted oscillations with a period of T ₁ = 2.878 Julian years (amplitude of 1.855″) and forced oscillations with periods of T ₂ = 27.201 days (15.304″), one stellar year (0.008″), half a year (0.115″), and the third of a year (0.0003″) (five harmonics altogether). A zero frequency solution has also been obtained. The effect of the Sun on these oscillations is two orders of magnitude less than that of the Earth. The libration in latitude is presented by five equations and, at pertrubations from the Earth, is described by two harmonics of unrestricted oscillations (T ₅ ≈ 74.180 Julian years, T ₆ ≈ 27.347 days) and one harmonic of forced oscillations (T ₃ = 27.212 days). The motion of the true pole is presented by the same harmonics, with the maximum deviation from the Cassini pole being 45.3″. The fifth (zero) frequency yields a stationary solution with a conic precession of the rotation axis (previously unknown). The third Cassini law has been proved. The amplitudes of unrestricted oscillations have been determined from comparison with observations. For the ratio $ \frac{{\sin I}} {{\sin \left( {I + i} \right)}} \approx 0.2311 $ \frac{{\sin I}} {{\sin \left( {I + i} \right)}} \approx 0.2311 , the theory gives 0.2319, which confirms the adequacy of the approach. Some statements of the previous theory are revised. Poinsot’s method is shown to be irrelevant in describing librations of the Moon. The Moon does not have free (Euler) oscillations; it has oscillations with a period of T ₅ ≈ 74.180 Julian years rather than T ≈ 148.167 Julian years.     

Geodesic study of regular Hayward black hole

This paper is devoted to study the geodesic structure of regular Hayward black hole. The timelike and null geodesics have been studied explicitly for radial and non-radial motion. For timelike and null geodesic in radial motion there exists analytical solution, while for non-radial motion the effective potential has been plotted, which investigates the position and turning points of the particle. It has been found that massive particle moving along timelike geodesics path are dragged towards the BH and continues to move around black hole in particular orbits.     

Global Coronal Seismology

Following the observation and analysis of large-scale coronal-wave-like disturbances, we discuss the theoretical progress made in the field of global coronal seismology. Using simple mathematical techniques we determine average values for the magnetic field together with a magnetic map of the quiet Sun. The interaction between global coronal waves and coronal loops allows us to study loop oscillations in a much wider context, i.e. we connect global and local coronal oscillations.     

Radiation pressure effects in the oscillations of compressible rotating homogeneous spheroids

Earlier models of compressible, rotating, and homogeneous ellipsoids with gas pressure are generalized to include the presence of radiation pressure. Under the assumptions of a linear velocity field of the fluid and a bounded ellipsoidal surface, the dynamical behaviour of these models can be described by ordinary differential equations. These equations are used to study the finite oscillations of massive radiative models with masses 10M and 30M in which the effects of radiation pressure are expected to be important.Models with two different degrees of equilibrium are chosen: an equilibrium (i.e., dynamically stable) model with an initial asymmetric inward velocity, and a nonequilibrium model with a nonequilibrium central temperature and which falls inwards from rest. For each of these two degrees of equilibrium, two initial configurations are considered: rotating spheroidal and nonrotating spherical models.From the numerical integration of the differential equations for these models, we obtain the time evolution of their principal semi-diametersa ₁ anda ₃, and of their central temperatures, which are graphically displayed by making plots of the trajectories in the (a ₁,a ₃) phase space, and of botha ₁ and the total central pressureP _c against time.It is found that in all the equilibrium radiative models (in which radiation pressure is taken into account), the periods of the oscillations of botha ₁ andP _c are longer than those of the corresponding nonradiative models, while the reverse is true for the nonequilibrium radiative models. The envelopes of thea ₁ oscillations of the equilibrium radiative models also have much longer periods; this result also holds for the nonequilibrium models whenever the envelope is well defined. Further, as compared to the nonradiative models, almost all the radiative models collapse to smaller volumes before rebouncing, with the more massive model undergoing a larger collapse and attaining a correspondingly larger peakP _c.When the mass is increased, the dynamical behavior of the radiative model generally becomes more nonperiodic. The ratio of the central radiation pressure to the central gas pressure, which is small for low mass models, increases with mass, and at the center of the more massive model, the radiation pressure can be comparable in magnitude to the gas pressure. In all the radiative models, the average periods as well as the average amplitudes of both thea ₁ andP _c oscillations also increase with mass.When either rotation or radiation pressure effects or both are included in the equilibrium nonradiative model, the period of the envelope of thea ₁ oscillations is increased. The presence of rotation in the equilibrium radiative model, however, decreases this period.Some astrophysical implications of this work are briefly discussed.     

Motion of a Rigid Body in a Tidal Field

We investigate the motion, near the equilibrium configurations, of an initially spinless rigid body subject to an external tidal field. Two cases are considered: when the center of mass of the body is at rest at the equilibrium point of the field generated by a generic mass distribution, and when it is placed on a circular orbit subject to a spherically symmetric potential. A complete analysis of the equilibrium configurations is carried out for both cases. First, we derive the conditions for the stable equilibria, and then we analyze the frequencies of oscillations around the equilibrium positions. In view of these results, we consider the problem of alignment of galaxies in clusters. After estimating the period of the oscillations induced on the galaxies by the tidal field of the cluster, we discuss the possible effect of resonances between stellar orbits inside the galaxy and the oscillations of the galaxy as a whole; this may be a mechanism responsible for producing an intracluster stellar population. This revised version was published online in July 2006 with corrections to the Cover Date.     

Local Fractional Frequency Shifts Used as Tracers of Magnetic Activity

Using data from SOI-MDI (Haber et al., 2000), we compute the local frequencies of high-degree p modes and f modes. The frequencies are obtained through ring-diagram mode fitting. The Dense-Pack data set consists of a mosaic of 189 overlapping tiles, each tracked separately at the surface rotation rate over 1664-min time intervals during the Dynamics Programs. Each tile is 16° square and the tile centers are separated by 7.5° in latitude and longitude. For each sampling day and for each tile, we have computed the frequency shift measured relative to the temporal and spatial average of the entire set of frequencies. The motion of active regions as they rotate across the solar disk is vividly traced by these measurements. Active regions appear as locations of large positive frequency shifts. If the shifts are averaged over the solar disk and are scaled down to the appropriate wave number regime, the magnitude and frequency dependence of the shifts are consistent with the measured changes in global oscillation frequencies that occur over the solar cycle. As with the frequency shifts of low-degree global oscillations, the frequency dependence of the shifts indicates that the physical phenomena inducing the shifts is confined to the surface layers of the Sun.