Dynamical tides in close binary systems

The aim of the present paper will be to investigate the circumstances under which an irreversible dissipation of the kinetic energy into heat is generated by the dynamical tides in close binary systems if (a) their orbit is eccentric; (b) the axial rotation of the components is not synchronized with the revolution; or (c) the equatorial planes are inclined to that of the orbit.In Section 2 the explicit form of the viscous dissipation function will be set up in terms of the velocity-components of spheroidal deformation arising from the tides; in Section 3, the principal partial tides contributing to the dissipation will be detailed; Section 4 will be devoted to a determination of the extent of stellar viscosity — both gas and radiative; while in the concluding Section 5 quantitative estimates will be given of the actual rate at which the kinetic energy of dynamical tides gets dissipated into heat by viscous friction in stellar plasma.The results disclose that the amount of heat produced per unit time by tidal interaction between components of actual close binaries equals only about 10^–10th part of their nuclear energy production; and cannot, therefore, affect the internal structure of evolution of the constituent stars to any appreciable extent. Moreover, it is shown that the kinetic energy of their axial rotation can be influenced by tidal friction only on a nuclear, rather than gravitational (Kelvin) time-scale — as long as plasma or radiative viscosity constitute the sole sources of dissipation. However, the emergence of turbulent viscosity in secondary components of late spectral types, which have evolved away from the Main Sequence, can accelerate the dissipation 10⁵–10⁶ times, and thus give rise to appreciable changes in the elements of the system (particularly, in the orbital periods) over time intervals of the order of 10⁵–10⁶ years. Lastly, it is pointed out that, in close binary systems consisting of a pair of white dwarfs, a dissipation of the kinetic energy through viscous tides in degenerate fermion-gas could produce enough heat to account, by itself, for the observed luminosity of such objects.     

Angular momentum and kinetic energy of rotating deformable bodies

The general equations of angular momentum and kinetic energy of a rotating deformable (or not rigid) body are discussed for a fixed and a rotating coordinate system. A new system of equations is developed for a deformable body of arbitrary form using the Lagrangian (vector) cisplacement up to the first order terms. The equations are, then, illustrated for a self-gravitating ceformable body perturbed by tides.     

The precession and nutation of deformable bodies,III

In preceding papers of this series (Kopal, 1968; 1969) the Eulerian equations have been set up which govern the precession and nutation of self-gravitating fluid globes of arbitrary structures in inertial coordinates (space-axes) as well as with respect to the rotating body axes; with due account being taken of the effects arising from equilibrium as well as dynamical tides.In Section 1 of the present paper, the explicit form of these equations is recapitulated for subsequent solations. Section 2 contains then a detailed discussion of the coplanar case (in which the equation of the rotating configuration and the plane of its orbit coincide with the invariable plane of the system); and small fluctuations in the angular velocity of axial rotation arising from the tidal breathing in eccentric binary systems are investigated.In Section 3, we consider the angular velocity of rotation about theZ-axis to be constant, but allow for finite inclination of the equator to the orbital plane. The differential equations governing such a problem are set up exactly in terms of the time-dependent Eulerian angles and , and their coefficients averaged over a cycle. In Section 4, these equations are linearized by the assumption that the inclinations of the equator and the orbit to the invariable plane of the system are small enough for their squares to be negligible; and the equations of motion reduced to their canonical form.The solution of these equations — giving the periods of precession and nutation of rotating components of close binary systems, as well as the rate of nodal regression which is synchronised with precession — are expressed in terms of the physical properties of the respective system and of its constituent components; while the concluding Section 6 contains a discussion of the results, in which the differences between the precession and nutation of rigid and fluid bodies are pointed out.     

Tidal dissipation at arbitrary eccentricity and obliquity

Expressions for tidal dissipation in a body in synchronous rotation at arbitrary orbital eccentricity and obliquity are derived. The rate of tidal dissipation for a synchronously rotating body is compared to that in a body in asymptotic nonsynchronous rotation.     

Tides in a body librating about a spin–orbit resonance: generalisation of the Darwin–Kaula theory

The Darwin-Kaula theory of bodily tides is intended for celestial bodies rotating without libration. We demonstrate that this theory, in its customary form, is inapplicable to a librating body. Specifically, in the presence of libration in longitude, the actual spectrum of Fourier tidal modes differs from the conventional spectrum rendered by the Darwin–Kaula theory for a nonlibrating celestial object. This necessitates derivation of formulae for the tidal torque and the tidal heating rate, that are applicable under libration. We derive the tidal spectrum for longitudinal forced libration with one and two main frequencies, generalisation to more main frequencies being straightforward. (By main frequencies we understand those emerging due to the triaxiality of the librating body.) Separately, we consider a case of free libration at one frequency (once again, generalisation to more frequencies being straightforward). We also calculate the tidal torque. This torque provides correction to the triaxiality-caused physical libration. Our theory is not self-consistent: we assume that the tidal torque is much smaller than the permanent-triaxiality-caused torque, so the additional libration due to tides is much weaker than the main libration due to the permanent triaxiality. Finally, we calculate the tidal dissipation rate in a body experiencing forced libration at the main mode, or free libration at one frequency, or superimposed forced and free librations.     

Dynamical tides in close binary systems

Differential equations governing the dynamical tides in close binary systems consisting of centrally condensed components of viscous gas are split up (Section 2) in their real and imaginary parts, the ratio of which defines the tidal lag. In Sections 3 and 4 these equations will be particularized to a case in which the central mass-point of each star is surrounded by an evanescent envelope the density of which decreases as the inverse square of the central distance. It is shown that self-gravitating configurations built up in accordance with this model are incapable of performing free nonradial oscillations with a frequency comprised between 0 ² ; but explicit expressions for forced oscillations representing dynamical tides are given for an arbitrary form of the external field of force. Equations for the imaginary components of the displacement, constructed for the same model in Section 4, disclose that if the viscosity of stellar material is identified with that of hydrogen plasma, the tidal lag due to a viscous dissipation of kinetic energy may produce dynamical effects, the cumulative outcome of which becomes appreciable on the Kelvin time-scale, but over short intervals of time their stationary photometric effects should be negligible. The latter can become observationally significant only for stars in which turbulent viscosity under near-adiabatic conditions becomes and important factor.     

On secular stability of tidally distorted stars of arbitrary structure

The aim of the present paper will be to develop a theory which should make it possible to investigate secular stability of close binary systems, consisting of tidally-distorted components of arbitrary internal structure, by a minimization of the potential energy of the system as a whole. In the second section which follows brief introductory remarks, appropriate expressions for the total potential energy of a close binary will be formulated. Section 3 will be concerned mainly with the nature of the tide-generating potential, and its effects on the shape of each star. In Section 4, the amplitudes of partial tides raised by this potential will be specified, for stars of arbitrary structure, correctly to terms of second order in superficial distortion; and in Section 5 we shall investigate the effects of interaction between rotation and tides to the same degree of approximation. The concluding Section 6 will then contain an explicit formulation of different constituents adding up to the total potential energy of the system, which can be used as a basis for its secular stability by the methods outlined already in our previous investigation (Kopal, 1973).     

The effects of viscous friction on the precession and nutation of celestial bodies

The aim of the present study has been to set the system of differential equations which govern the precession and nutation of self-gravitating globes of compressible viscous fluid, due to the attraction exerted on the rotating configuration by its companion; and to construct their approximate solution which are correct to terms of the second order in small dependent variables of the problem. Section 2 contains an explicit formulation of the effects of viscosity arising in this connection, given exactly as far as the viscosity remains a function of radial distancer only; but irrespective of its magnitude. In Section 3 the equations of motion will be linearized for the case of near-circular orbits and small inclinations andi of the equator of the rotating configuration, and of its orbital plane, to the invariable plane of the system; while in Section 4 further simplifications will be introduced which are legitimate for studies of secular (or long-periodic) motions of the nodes and inclinations. The actual solutions of so simplified a system of equations are constructed in Section 5; and these represent a generalization of the results obtained in our previous investigation (Kopal, 1969) of the inviscid case.The physical significance of the new results will be discussed in the concluding Section 6. It is demonstrated that the axes of rotation of deformable components in close binary systems are initially inclined to the orbital plane, viscous dissipation produced by dynamical tides will tend secularly to rectify their positions until perpendicularity to the orbital plane has been established, and the equators as well as orbit made to coincide with the invariable plane of the system-in a similar manner as other effects of tidal friction are bound eventually to synchronize the velocity of axial rotation with that of orbital revolution in the course of time.An application of the results of the present study to the dynamics of the Earth-Moon system discloses that the observed inclination of 1°.5 of the lunar equator to the ecliptic cannot be regarded as being secularly constant, but representing the present deviations from perpendicularity of oscillatory motion of very long period.The Lunar Science Institute is operated by the Universities Space Research Association under Contract No. NSR-09-051-001 with the National Aeronautics and Space Administration. This paper constitutes the Lunar Science Institute Contribution No. 85.     

The precession and nutation of deformable bodies

The aim of the present paper will be to derive from the fundamental equations of hydrodynamics the explicit form of the Eulerian equations which govern the motion about the centre of gravity of self-gravitating bodies, consisting of compressible fluid of arbitrary viscosity, in an arbitrary external field of force. If the problem is particularized so that the external field of force represents the attaction of the sun and the moon, this motion would represent the luni-solar precession and nutation of a fluid viscous earth; if, on the other hand, the external field of force were governed by the earth (and the sun), the motion would define the physical librations of the moon regarded as a deformable body. The same equations are, moreover, equally applicable to the phenomena of precession and nutation of rotating fluid components in close binary systems, distorted by mutual tidal action; and the present paper contains the first formulation of the effects of viscosity on such phenomena.Investigation supported in part by the U.S. National Aeronautics and Space Administration under Contract No. NASW-1470.     

An observational test for the origin of the Titan-Hyperion orbital resonance

If Hyperion's radius is near the upper limit of recent estimates, and tidal dissipation in Hyperion is reasonably well represented by a frequency-independent Q ? 2–300, finding Hyperion rotating in the 3:2 spin-orbit resonance like Mercury would imply a primordial origin for the Titan-Hyperion 4:3 orbital resonance. Independent of this test, observation of Hyperion's spin rate will place an upper bound on the average tidal effective Q for the satellite as a function of its assumed initial angular velocity.     

The generalized non-conservative model of a 1-planet system revisited

We study the long-term dynamics of a planetary system composed of a star and a planet. Both bodies are considered as extended, non-spherical, rotating objects. There are no assumptions made on the relative angles between the orbital angular momentum and the spin vectors of the bodies. Thus, we analyze full, spatial model of the planetary system. Both objects are assumed to be deformed due to their own rotations, as well as due to the mutual tidal interactions. The general relativity corrections are considered in terms of the post-Newtonian approximation. Besides the conservative contributions to the perturbing forces, there are also taken into account non-conservative effects, i.e., the dissipation of the mechanical energy. This dissipation is a result of the tidal perturbation on the velocity field in the internal zones with non-zero turbulent viscosity (convective zones). Our main goal is to derive the equations of the orbital motion as well as the equations governing time-evolution of the spin vectors (angular velocities). We derive the Lagrangian equations of the second kind for systems which do not conserve the mechanical energy. Next, the equations of motion are averaged out over all fast angles with respect to time-scales characteristic for conservative perturbations. The final equations of motion are then used to study the dynamics of the non-conservative model over time scales of the order of the age of the star. We analyze the final state of the system as a function of the initial conditions. Equilibria states of the averaged system are finally discussed.     

Slowly rotating relativistic stars

Equations are given which determine the moment of inertia of a rotating relativistic fluid star to second order in the angular velocity with no other approximation being made. The equations also determine the moment of inertia of matter located between surfaces of constant density in a rotationally distorted star; for example, the moments of inertia of the crust and core of a rotationally distorted neutron star can be calculated in this way. The method is applied ton=3/2 relativistic polytropes and to neutron star models constructed from the Baym-Bethe-Pethick-Sutherland-Pandharipande equation of state. Supported in part by the National Science Foundation. Alfred P. Sloan Research Fellow.     

Tidal evolution in close binary systems

The aim of the present paper will be to give a mathematical outline of the theory of tidal evolution in close binary systems of secularly constant total momentum — an evolution activated by viscous friction of dynamical tides raised by the two components on each other. The first section contains a general outline of the problem; and in Section 2 we shall establish the basic expressions for the energy and momenta of close binaries consisting of components of arbitrary internal structure. In Section 3 we shall investigate the maximum and minimum values of the energy (kinetic and potential) which such systems can attain for given amount of total momentum; while in Section 4 we shall compare these results with the actual facts encountered in binaries with components whose internal structure (and, therefore, rotational momenta) are known to us from evidence furnished by the observed rates of apsidal advance.The results show that all such systems — be these of detached or semi-detached type — disclose that more than 99% of their total momenta are stored in the orbital momentum. The sum of the rotational momenta of the constituent components amounts to less than a percent of the total — a situation characteristic of a state close to the minimum energy for given total momentum. This appears, moreover, to be true not only of the systems with both components on the Main Sequence, but also of those possessing evolved components in contact with their Roche limits.Under such conditions, a synchronism between rotation and revolution (characteristic of both extreme states of maximum and minimum energy) is not only possible, but appears to have been actually approached — if not attained — in the majority of cases. In other words, it would appear that — in at least a large majority of known cases — the existing close binaries have already attained orbits of maximum distension consistent with their momenta; and tidal evolution alone can no longer increase the present separations of the components to any appreciable extent.The virtual absence, in the sky, of binary systems intermediate between the stages of maximum and minimum energy for given momentum leads us to conjecture that the process of dynamical evolution activated by viscous tides may enroll on a time-scale which is relatively short in comparison with their total age — even for systems like Y Cygni or AG Persei, whose total age can scarcely exceed 10⁷ yr. A secular increase of the semi-major axes of relative orbits is dynamically coupled with a corresponding variation in the velocity of axial rotation of both components through the tidal lag arising from the viscosity of stellar material. The differential equations of so coupled a system are given in Section 5; but their solution still constitutes a task for the future.The Lunar Science Institute Contribution No. 90. The Lunar Science Institute is operated by the Universities Space Research Association under Contract No. NSR 09-051-001 with the National Aeronautics and Space Administration.     

Tidal dissipation in rotating fluid bodies: a simplified model

We study the tidal forcing, propagation and dissipation of linear inertial waves in a rotating fluid body. The intentionally simplified model involves a perfectly rigid core surrounded by a deep ocean consisting of a homogeneous incompressible fluid. Centrifugal effects are neglected, but the Coriolis force is considered in full, and dissipation occurs through viscous or frictional forces. The dissipation rate exhibits a complicated dependence on the tidal frequency and generally increases with the size of the core. In certain intervals of frequency, efficient dissipation is found to occur even for very small values of the coefficient of viscosity or friction. We discuss the results with reference to wave attractors, critical latitudes and other features of the propagation of inertial waves within the fluid, and comment on their relevance for tidal dissipation in planets and stars.     

Effects of rotation on internal structure of the stars

The aim of the present paper will be to establish the explicit form of the equations which govern the internal structure of stars rotating with constant angular velocity formulated in terms of Clairaut coordinates (cf. Kopal, 1980) in which the radial coordinate is replaced by the total potential, which for equilibrium configurations remains constant over distorted level surfaces. The introductory Section 1 contains an account of previous work on rotating stars, commencing with Milne (1923), von Zeipel (1924) and Chandrasekhar (1933), who all employed orthogonal coordinates for their analysis. In Section 2 we shall apply to this end the curvilinear Clairaut coordinates introduced already in our previous work (cf. Kopal, 1980, 1981); and although these are not orthogonal, this disadvantage is more than offset by the fact that, in their terms, the fundamental equation of our problem will assume the form of ordinary differential equations, subject to very simple boundary conditions. The explicit form of these equations — exact to terms of fourth order in surficial distortion caused by centrifugal force—will be obtained in Section 3; while in the concluding Section 4 these will be particularized (for the sake of comparison with work of previous investigators) to stars of initially polytropic structure. These will prove to be much simpler in Clairaut coordinates than they were in any previously used frame of reference. Lastly, in Appendix A we shall present the explicit forms, in Clairaut coordinates, of the differential operators which were needed to establish the results given in Sections 3–4; while Appendix B will summarize other auxiliary algebraic relations of which use was made to formulate our fourth-order theory developed in Section 3.     

Vibrational stability of the components of close binary systems

The aim of the present paper will be to extend our previous investigation of the vibrational stability of rotating configurations (Kopal, 1981) to a similar investigation of the stability of the components of close binary systems which not only rotate, but also distort each other by tidal action. To this end, differential equations which govern first-order oscillations of arbitrary spherical-harmonic symmetry will be set up in Clairaut coordinates in which the radial coordinate is replaced by the potential which remains constant over level surfaces of equilibrium configurations; introduced by us in an earlier paper (Kopal, 1980), and their form detailed for surface distorted by second-, third-, and fourth-harmonic tides raised by the external mass; and their boundary conditions established. A solution of such differential boundary-value problems arising in connection with the stars of arbitrary structure remains, of course, a task for automatic computers. It may only be added that the tide-generating potential Ψ _T established in this paper should enable us to study, by the same method, not only free, but also forced oscillations of the components of close binary systems, arising from orbital eccentricity of the respective couples, dynamical tides, or other causes likely to be operative in such systems.     

潮汐能量耗散与地月系统演化

对地月系统而言, 在很大程度上角动量守恒是正确的. 地月距离的变化主要是受到月球引起的潮汐能量耗散的影响. 根据月球的平均运动和它的长期加速度, 就可以计算出月潮能量耗散的数值. 海洋是潮汐能量耗散的主要区域. 由于潮汐的高度正比于月球对潮汐隆起的万有引力, 由此可导出总的月球潮汐摩擦力正比于月球平均运动的平方. 如果采用月球平均加速度数值-20.72$''\cdot$cy^-2, 就可以推算出35亿年来地月之间的距离以及回归年日数和朔望月日数的演化. 此理论结果与古生物钟的数据进行比对, 两者符合较好.     

On the dissipation rates for Alfvén waves in the solar transition region

The present paper concerns Alfvén waves, in a resistive and viscous atmosphere, under a steep temperature gradient (Section 1). The dissipative Alfvén wave equation is deduced assuming uniform vertical background magnetic field, and allowing for arbitrary profiles of Alfvén speed, and viscous and resistive diffusivities as functions of altitude (Section 2). A three-parameter family of temperature profiles, allowing for independent choice of initial and asymptotic temperature, and of initial temperature gradient, is used to re-write the wave equation, with the temperature as the independent variable, instead of altitude (Section 3). It is shown that, for the conditions prevailing in the solar transition region between the chromosphere and corona, two approximations of the dissipative wave equations may be considered, the simplest leading to solution in terms of Gaussian hypergeometric functions (Section 4). The exact analytical solution allows calculation of the (i) velocity and (ii) magnetic field perturbations, (iii) kinetic, (iv) magnetic and (v) total energy density, (vi) energy flux, (vii) rate-of-strain and (viii) electric current, and (ix) viscous, (x) resistive and (xi) total rate of dissipation (Section 5). These are plotted versus temperature, across the transition region from the chromosphere to the corona, for the quiet and active Sun (Section 6). The feasibility of heating of the transition region by dissipation of Alfvén waves is discussed (Section 7), by comparing empirical heating rates, with theoretical values for a range of physical conditions, including initial velocity perturbations 5 to 15 km s ^–1, background magnetic field 12 to 120 G, wave periods 60 to 300 s, thickness of the transition region 100 to 300 km, resistive and anomalous diffusivities to 100 and viscous and turbulent diffusivities to 100 . The conclusion is that dissipation of Alfvén waves is not an effective heating mechanism for the transition region and corona, although it may be for the chromosphere (see Campos and Mendes, 1995, and references therein).     

Effects of rotation and tidal distortions on the shapes of radial velocity curves of polytropic models of pulsating variable stars

Anharmonic oscillations of rotating stars have been studied by various authors in literature to explain the observed features of certain variable stars. However, there is no study available in literature that has discussed the combined effect of rotation and tidal distortions on the anharmonic oscillations of stars. In this paper, we have created a model to determine the effect of rotation and tidal distortions on the anharmonic radial oscillations associated with various polytropic models of pulsating variable stars. For this study we have used the theory of Rosseland to obtain the anharmonic pulsation equation for rotationally and tidally distorted polytropic models of pulsating variable stars. The main objective of this study is to investigate the effect of rotation and tidal distortions on the shapes of the radial velocity curves for rotationally and tidally distorted polytropic models of pulsating variable stars. The results of the present study show that the rotational effects cause more deviations in the shapes of radial velocity curves of pulsating variable stars as compared to tidal effects.