We address the expressions for the rates of the Keplerian orbital elements within a two-body problem perturbed by the tides in both partners. The formulae for these rates have appeared in the literature in various forms, at times with errors. We reconsider, from scratch, the derivation of these rates and arrive at the Lagrange-type equations which, in some details, differ from the corresponding equations obtained previously by Kaula (1964). We also write down detailed expressions for $da/dt$, $de/dt$ and $di/dt$, to order $e^4$. They differ from Kaula's expressions which contain a redundant factor of $M/(M+M^{\prime}),$ with $M$ and $M^{\prime}$ being the masses of the primary and the secondary. As Kaula was interested in the Earth-Moon system, this redundant factor was close to unity and was unimportant in his developments. This factor, however, must be reinstated when Kaula's theory is applied to a binary composed of partners of comparable masses. We have found that, while it is legitimate to simply sum the primary's and secondary's inputs in $da/dt$ or $de/dt$, this is not the case for $di/dt$. So our expression for $di/dt$ differs from that of Kaula in two regards. First, the contribution due to the dissipation in the secondary averages out when the apsidal precession is uniform. Second, we have obtained an additional term which emerges owing to the conservation of the angular momentum: a change in the inclination of the orbit causes a change of the primary's plane of equator.

Observational evidence of white dwarf planetary systems is dominated by the remains of exo-asteroids through accreted metals, debris discs, and orbiting planetesimals. However, exo-planets in these systems play crucial roles as perturbing agents, and can themselves be perturbed close to the white dwarf Roche radius. Here, we illustrate a procedure for computing the tidal interaction between a white dwarf and a near-spherical solid planet. This method determines the planet's inward and/or outward drift, and whether the planet will reach the Roche radius and be destroyed. We avoid constant tidal lag formulations and instead employ the self-consistent secular Darwin-Kaula expansions from Bou\'{e} & Efroimsky (2019), which feature an arbitrary frequency dependence on the quality functions. We adopt wide ranges of dynamic viscosities and spin rates for the planet in order to straddle many possible outcomes, and provide a foundation for the future study of individual systems with known or assumed rheologies. We find that: (i) massive Super-Earths are destroyed more readily than minor planets (such as the ones orbiting WD 1145+017 and SDSS J1228+1040), (ii) low-viscosity planets are destroyed more easily than high-viscosity planets, and (iii) the boundary between survival and destruction is likely to be fractal and chaotic.

We perform numerical simulations of the TRAPPIST-1 system of seven exoplanets orbiting a nearby M dwarf, starting with a previously suggested stable configuration. The long-term stability of this configuration is confirmed, but the motion of planets is found to be chaotic. The eccentricity values are found to vary within finite ranges. The rates of tidal dissipation and tidal evolution of orbits are estimated, assuming an Earth-like rheology for the planets. We find that under this assumption the planets b, d, e were captured in the 3:2 or higher spin-orbit resonances during the initial spin-down but slipped further down into the 1:1 resonance. Dependent on its rheology, the innermost planet b may be captured in a stable pseudosynchronous rotation. Non-synchronous rotation ensures higher levels of tidal dissipation and internal heating. The positive feedback between the viscosity and the dissipation rate -- and the ensuing runaway heating -- are terminated by a few self-regulation processes. When the temperature is high and the viscosity is low enough, the planet spontaneously leaves the 3:2 resonance. Further heating is stopped either by passing the peak dissipation or by the emergence of partial melt in the mantle. In the post-solidus state, the tidal dissipation is limited to the levels supported by the heat transfer efficiency. The tides on the host star are unlikely to have had a significant dynamical impact. The tides on the synchronized inner planets tend to reduce these planets' orbital eccentricity, possibly contributing thereby to the system's stability.

Highly volcanic exoplanets, which can be variously characterized as 'lava worlds', 'magma ocean worlds', or 'super-Ios' are high priority targets for investigation. The term 'lava world' may refer to any planet with extensive surface lava lakes, while the term 'magma ocean world' refers to planets with global or hemispherical magma oceans at their surface. 'Highly volcanic planets', including super-Ios, may simply have large, or large numbers of, active explosive or extrusive volcanoes of any form. They are plausibly highly diverse, with magmatic processes across a wide range of compositions, temperatures, activity rates, volcanic eruption styles, and background gravitational force magnitudes. Worlds in all these classes are likely to be the most characterizable rocky exoplanets in the near future due to observational advantages that stem from their preferential occurrence in short orbital periods and their bright day-side flux in the infrared. Transit techniques should enable a level of characterization of these worlds analogous to hot Jupiters. Understanding processes on highly volcanic worlds is critical to interpret imminent observations. The physical states of these worlds are likely to inform not just geodynamic processes, but also planet formation, and phenomena crucial to habitability. Volcanic and magmatic activity uniquely allows chemical investigation of otherwise spectroscopically inaccessible interior compositions. These worlds will be vital to assess the degree to which planetary interior element abundances compare to their stellar hosts, and may also offer pathways to study both the very young Earth, and the very early form of many silicate planets where magma oceans and surface lava lakes are expected to be more prevalent. We suggest that highly volcanic worlds may become second only to habitable worlds in terms of both scientific and public long-term interest.

The geophysics of extrasolar planets is a scientific topic often regarded as standing largely beyond the reach of near-term observations. This reality in no way diminishes the central role of geophysical phenomena in shaping planetary outcomes, from formation, to thermal and chemical evolution, to numerous issues of surface and near-surface habitability. We emphasize that for a balanced understanding of extrasolar planets, it is important to look beyond the natural biases of current observing tools, and actively seek unique pathways to understand exoplanet interiors as best as possible during the long interim prior to a time when internal components are more directly accessible. Such pathways include but are not limited to: (a) enhanced theoretical and numerical modeling, (b) laboratory research on critical material properties, (c) measurement of geophysical properties by indirect inference from imprints left on atmospheric and orbital properties, and (d) the purpose-driven use of Solar System object exploration expressly for its value in comparative planetology toward exoplanet-analogs. Breaking down barriers that envision local Solar System exploration, including the study of Earth's own deep interior, as separate from and in financial competition with extrasolar planet research, may greatly improve the rate of needed scientific progress for exoplanet geophysics. As the number of known rocky and icy exoplanets grows in the years ahead, we expect demand for expertise in 'exogeoscience' will expand at a commensurately intense pace. We highlight key topics, including: how water oceans below ice shells may dominate the total habitability of our galaxy by volume, how free-floating nomad planets may often attain habitable subsurface oceans supported by radionuclide decay, and how deep interiors may critically interact with atmospheric mass loss via dynamo-driven magnetic fields.

Internal dissipation in a tidally perturbed librating body differs from the tidal dissipation in a steadily spinning rotator. First, libration changes the spectral distribution of tidal damping across the tidal modes, as compared to the case of steady spin. This changes both the tidal heating rate and the tidal torque. Second, while a non-librating rotator experiences alternating deformation only due to the potential force exerted on it by the perturber, a librating body is also subject to a toroidal force proportional to the angular acceleration. Third, while the centrifugal force in a steadily spinning body renders only a permanent deformation, in a librating body this force contains two alternating components $-$ one radial, another a degree-2 potential force. Both contribute to heating, as well as to the tidal torque and potential. We build a formalism to describe dissipation in a homogeneous terrestrial body performing small-amplitude libration in longitude. This formalism incorporates a linear rheological law defining the response of the material to forcing. While the formalism can work with an arbitrary linear rheology, we consider a simple example of a Maxwell material. We show that, independent of rheology, the forced libration in longitude can provide a considerable and even leading input in the tidal heating. Based on the observed parameters, this input amounts to 52% in Phobos, 33% in Mimas, 23% in Enceladus, and 96% in Epimetheus. This supports the hypothesis by Makarov & Efroimsky (2014) that the additional damping due to forced libration may have participated in the early heating up of some moons. As one possibility, a moon could have been chipped by collisions $-$ whereby it acquired a higher triaxiality and a higher forced-libration magnitude and, consequently, a higher heating rate. After the moon warms up, its triaxiality reduces, and so does the tidal heating.

In the preceding paper (Efroimsky 2017), we derived an expression for the tidal dissipation rate in a homogeneous near-spherical Maxwell body librating in longitude. Now, by equating this expression to the outgoing energy flux due to the vapour plumes, we estimate the mean tidal viscosity of Enceladus, under the assumption that the Enceladean mantle behaviour is Maxwell. This method yields a value of $\,0.24\times 10^{14}\;\mbox{Pa~s}\,$ for the mean tidal viscosity, which is very close to the viscosity of ice near the melting point.

Perturbations of all sorts destabilise the rotation of a small body and leave it in a non-principal spin state. In such a state, the body experiences alternating stresses generated by the inertial forces. This yields nutation relaxation, i.e., evolution of the spin towards the principal rotation about the maximal-inertia axis. Knowledge of the timescales needed to damp the nutation is crucial in studies of small bodies' dynamics. In the literature hitherto, nutation relaxation has always been described with aid of an empirical quality factor $\,Q\,$ introduced to parameterise the energy dissipation rate. Among the drawbacks of this approach was its inability to describe the dependence of the relaxation rate upon the current nutation angle. This inability stemmed from our lack of knowledge of the quality factor's dependence on the forcing frequency. In this article, we derive our description of nutation damping directly from the rheological law obeyed by the material. This renders us the nutation damping rate as a function of the current nutation angle, as well as of the shape and the rheological parameters of the body. In contradistinction from the approach based on an empirical $\,Q\,$-factor, our development gives a zero damping rate in the spherical-shape limit. Our method is generic and applicable to any shape and to any linear rheological law. However, to simplify the developments, here we consider a dynamically oblate rotator with a Maxwell rheology.

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 non-librating 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.

We use a damped mass-spring model within an N-body code to simulate the tidal evolution of the spin and orbit of a self-gravitating viscoelastic spherical body moving around a point-mass perturber. The damped mass-spring model represents a Kelvin-Voigt viscoelastic solid. We measure the tidal quality function (the dynamical Love number $\,k_2\,$ divided by the tidal quality factor $\,Q\,$) from the numerically computed tidal drift of the semimajor axis of the binary. The shape of $\,k_2/Q\,$, as a function of the principal tidal frequency, reproduces the kink shape predicted by Efroimsky (2012a; CeMDA 112$\,:\,$283) for the tidal response of near-spherical homogeneous viscoelastic rotators. We demonstrate that we can directly simulate the tidal evolution of spinning viscoelastic objects. In future, the mass-spring N-body model can be generalised to inhomogeneous and/or non-spherical bodies.

A formula for the tidal dissipation rate in a spherical body is derived from first principles, to correct some mathematical inaccuracies found in the literature. The development is combined with the Darwin-Kaula formalism for tides. Our intermediate results are compared with those by Zschau (1978) and Platzman (1984). When restricted to the special case of an incompressible spherical planet spinning synchronously without libration, our final formula can be compared with the commonly used expression from Peale & Cassen (1978, Eqn. 31). The two turn out to differ. In our expression, the contributions from all Fourier modes are positive-definite, this not being the case of the formula from Ibid. (The presence of negative terms in their formula was noticed by Makarov 2013.) Examples of application of our expression for the tidal damping rate are provided in the work by Makarov and Efroimsky (2014).

The rate of tidal evolution of asteroidal binaries is defined by the dynamical Love numbers divided by quality factors. Common is the (often illegitimate) approximation of the dynamical Love numbers with their static counterparts. As the static Love numbers are, approximately, proportional to the inverse rigidity, this renders a popular fallacy that the tidal evolution rate is determined by the product of the rigidity by the quality factor: $\,k_l/Q\propto 1/(\mu Q)\,$. In reality, the dynamical Love numbers depend on the tidal frequency and all rheological parameters of the tidally perturbed body (not just rigidity). We demonstrate that in asteroidal binaries the rigidity of their components plays virtually no role in tidal friction and tidal lagging, and thereby has almost no influence on the intensity of tidal interactions (tidal torques, tidal dissipation, tidally induced changes of the orbit). A key quantity that determines the tidal evolution is a product of the effective viscosity $\,\eta\,$ by the tidal frequency $\,\chi\,$. The functional form of the torque's dependence on this product depends on who wins in the competition between viscosity and self-gravitation. Hence a quantitative criterion, to distinguish between two regimes. For higher values of $\,\eta\chi\,$ we get $\,k_l/Q\propto 1/(\eta\chi)\;$; $\,$while for lower values we obtain $\,k_l/Q\propto \eta\chi\,$. Our study rests on an assumption that asteroids can be treated as Maxwell bodies. Applicable to rigid rocks at low frequencies, this approximation is used here also for rubble piles, due to the lack of a better model. In the future, as we learn more about mechanics of granular mixtures in a weak gravity field, we may have to amend the tidal theory with other rheological parameters, ones that do not show up in the description of viscoelastic bodies.

In Efroimsky & Makarov (2014), we derived from the first principles a formula for the tidal heating rate in a tidally perturbed homogeneous sphere. We compared it with the formulae used in the literature, and pointed out the differences. Using this result, we now present three case studies - Mercury, Kepler-10b, and a triaxial Io. A very sharp frequency-dependence of k2/Q near spin-orbit resonances yields a similarly sharp dependence of k2/Q on the spin rate. This indicates that physical libration may play a major role in tidal heating of synchronously rotating bodies. The magnitude of libration in the spin rate being defined by the planet's triaxiality, the latter should be a factor determining the dissipation rate. Other parameters equal, a synchronously rotating body with a stronger triaxiality should generate more heat than a similar body of a more symmetrical shape. Further in the paper, we discuss scenarios where initially triaxial objects melt and lose their triaxiality. Thereafter, dissipation in them becomes less intensive; so the bodies freeze. The tidal bulge becomes a new permanent figure, with a new triaxiality lower than the original. In the paper, we also derive simplified, approximate expressions for dissipation rate in a rocky planet of the Maxwell rheology, with a not too small Maxwell time. The three expressions derived pertain to the cases of a synchronous spin, a 3:2 resonance, and a nonresonant rotation; so they can be applied to most close-in super-Earth exoplanets detected thus far. In such bodies, the rate of tidal heating outside of synchronous rotation is weakly dependent on the eccentricity and obliquity, provided both these parameters are small or moderate. According to our calculation, Kepler-10b could hardly survive the great amount of tidal heating without being synchronised, circularised and also reshaped through a complete or partial melt-down.

While it is accepted that the eccentricity of Mercury (0.206) favours entrapment into the 3:2 spin-orbit resonance, open is the question how and when the capture took place. A recent work by Makarov (2012) has demonstrated that trapping into this resonance is certain if the eccentricity is larger than 0.2, provided that we use a realistic tidal model, the one which is based on the Darwin-Kaula expansion of the tidal torque. The physics-based tidal model changes dramatically the statistics of the possible final spin states. First, we discover that after only one encounter with the spin-orbit 3:2 resonance this resonance becomes the most probable end-state. Second, if a capture into this (or any other) resonance takes place, the capture becomes final, several crossings of the same state being forbidden by our model. Third, within our model the trapping of Mercury happens much faster than previously believed: for most histories, 10 - 20 Myr are sufficient. Fourth, even a weak laminar friction between the solid mantle and a molten core would most likely result in a capture in the 2:1 or even higher resonance. So the principal novelty of our paper is that the 3:2 end-state is more ancient than the same end-state obtained when the constant time lag model is employed. The swift capture justifies our treatment of Mercury as a homogeneous, unstratified body whose liquid core had not yet formed by the time of trapping. We also provide a critical analysis of the hypothesis by Wieczorek et al. (2012) that the early Mercury might had been retrograde, whereafter it synchronised its spin and then accelerated it to the 3:2 resonance. Accurate processing of the available data on cratering does not support that hypothesis, while the employment of a realistic rheology invalidates a key element of the hypothesis, an intermediate pseudosynchronous state needed to spin-up to the 3:2 resonance.

GJ 581d is a potentially habitable super-Earth in the multiple system of exoplanets orbiting a nearby M dwarf. We investigate this planet's long-term dynamics, with an emphasis on its probable final rotation states acquired via tidal interaction with the host. The published radial velocities for the star are re-analysed with a benchmark planet detection algorithm, to confirm that there is no evidence for the recently proposed two additional planets (f and g). Limiting the scope to the four originally detected planets, we assess the dynamical stability of the system and find bounded chaos in the orbital motion. For the planet d, the characteristic Lyapunov time is 38 yr. Long-term numerical integration reveals that the system of four planets is stable, with the eccentricity of the planet d changing quasi-periodically in a tight range around 0.27, and with its semimajor axis varying only a little. The spin-orbit interaction of GJ 581d with its host star is dominated by the tides exerted by the star on the planet. We model this interaction, assuming a terrestrial composition of the mantle. Besides the customarily included secular parts of the triaxiality-caused and tidal torques, we also include these torques' oscillating components. It turns out that, dependent on the mantle temperature, the planet gets trapped into the 2:1 or an even higher spin-orbit resonance. It is very improbable that the planet could have reached the 1:1 resonance. This enhances the possibility of the planet being suitable for sustained life.

Tidal torques play a key role in rotational dynamics of celestial bodies. They govern these bodies' tidal despinning, and also participate in the subtle process of entrapment of these bodies into spin-orbit resonances. This makes tidal torques directly relevant to the studies of habitability of planets and their moons. Our work begins with an explanation of how friction and lagging should be built into the theory of bodily tides. Although much of this material can be found in various publications, a short but self-consistent summary on the topic has been lacking in the hitherto literature, and we are filling the gap. After these preparations, we address a popular concise formula for the tidal torque, which is often used in the literature, for planets or stars.We explain why the derivation of this expression, offered in the paper by Goldreich (1966; AJ 71, 1 - 7) and in the books by Kaula (1968, eqn. 4.5.29), and Murray & Dermott (1999, eqn. 4.159), implicitly sets the time lag to be frequency independent. Accordingly, the ensuing expression for the torque can be applied only to bodies having a very special (and very hypothetical) rheology which makes the time lag frequency independent, i.e, the same for all Fourier modes in the spectrum of tide. This expression for the torque should not be used for bodies of other rheologies. Specifically, the expression cannot be combined with an extra assertion of the geometric lag (or the phase lag) being constant, because at finite eccentricities the said assumption is incompatible with the constant-time-lag condition.

We reexamine the popular belief that a telluric planet or satellite on an eccentric orbit can, outside a spin-orbit resonance, be captured in a quasi-static tidal equilibrium called pseudosynchronous rotation. The existence of such configurations was deduced from oversimplified tidal models assuming either a constant tidal torque or a torque linear in the tidal frequency. A more accurate treatment requires that the torque be decomposed into the Darwin-Kaula series over the tidal modes, and that this decomposition be combined with a realistic choice of rheological properties of the mantle. This development demonstrates that there exist no stable equilibrium states for solid planets and moons, other than spin-orbit resonances.

Spin-orbit coupling is often described in the "MacDonald torque" approach which has become the textbook standard. Within this method, a concise expression for the additional tidal potential, derived by MacDonald (1964; Rev. Geophys. 2, 467), is combined with an assumption that the Q factor is frequency-independent (i.e., that the geometric lag angle is constant in time). This makes the approach unphysical because MacDonald's derivation of the said formula was implicitly based on keeping the time lag frequency-independent, which is equivalent to setting Q to scale as the inverse tidal frequency. The contradiction requires the MacDonald treatment of both non-resonant and resonant rotation to be rewritten. The non-resonant case was reconsidered by Efroimsky & Williams (2009; CMDA 104, 257), in application to spin modes distant from the major commensurabilities. We continue this work by introducing the necessary alterations into the MacDonald-torque-based model of falling into a 1:1 resonance. (For the original version of the model, see Goldreich 1966; AJ 71, 1.) We also study the effect of the triaxiality on both circulating and librating rotation near the synchronous state. Circulating rotation may evolve toward the libration region or toward a spin rate larger than synchronous (pseudosynchronous spin). Which behaviour depends on the eccentricity, the triaxiality of the primary, and the mass ratio of the secondary and primary bodies. The spin evolution will always stall for the oblate case. For small-amplitude librations, expressions are derived for the libration frequency, damping rate, and average orientation. However, the stability of pseudosynchronous spin hinges upon the dissipation model. Makarov and Efroimsky (2012; arXiv:1209.1616) have found that a more realistic dissipation model than the corrected MacDonald torque makes pseudosynchronous spin unstable.

While the seismic quality factor and phase lag are defined solely by the bulk properties of the mantle, their tidal counterparts are determined both by the bulk properties and self-gravitation of a body as a whole. For a qualitative estimate, we model the body with a homogeneous sphere and express the tidal phase lag through the lag in a sample of material. Although simplistic, our model is sufficient to understand that the lags are not identical. The difference emerges because self-gravitation pulls the tidal bulge down. At low frequencies, this reduces strain and makes tidal damping less efficient in larger bodies. At high frequencies, competition between self-gravitation and rheology becomes more complex, though for sufficiently large superearths the same rule works: the larger the body, the weaker tidal damping in it. Being negligible for small terrestrial planets and moons, the difference between the seismic and tidal lagging (and likewise between the seismic and tidal damping) becomes very considerable for superearths. In those, it is much lower than what one might expect from using a seismic quality factor. The tidal damping rate deviates from the seismic damping rate especially in the zero-frequency limit, and this difference takes place for bodies of any size. So the equal in magnitude but opposite in sign tidal torques, exerted on one another by the primary and the secondary, go smoothly through zero as the secondary crosses the synchronous orbit. We describe the mantle rheology with the Andrade model, allowing it to lean towards the Maxwell model at the lowest frequencies. To implement this additional flexibility, we reformulate the Andrade model by endowing it with a free parameter which is the ratio of the anelastic timescale to the viscoelastic Maxwell time of the mantle. Some uncertainty in this parameter's frequency-dependence does not influence our principal conclusions.

Description of tides is based on the form of dependence of the geometric lag on the tidal frequency. Some authors assume the lag angle to be constant, others set it to be linear in the frequency. The actual dependence of the lag on the frequency is complicated and is determined by the planet's rheology. A particular form of this dependence will fix the form of the frequency dependence of the tidal quality factor Q. Since at present we know the frequency-dependence of the quality factor, we can reverse our line of reasoning and obtain the appropriate frequency-dependence of the lag. Employment of a realistic frequency-dependence for Q renders considerable changes in timescales defined by tidal dynamics.

We point out that the MacDonald formula for body-tide torques is valid only in the zeroth order of e/Q, while its time-average is valid in the first order. So the formula cannot be used for analysis in higher orders of e/Q. This necessitates corrections in the theory of tidal despinning and libration damping. We prove that when the inclination is low and phase lags are linear in frequency, the Kaula series is equivalent to a corrected version of the MacDonald method. The correction to MacDonald's approach would be to set the phase lag of the integral bulge proportional to the instantaneous frequency. The equivalence of descriptions gets violated by a nonlinear frequency-dependence of the lag. We explain that both the MacDonald- and Darwin-torque-based derivations of the popular formula for the tidal despinning rate are limited to low inclinations and to the phase lags being linear in frequency. The Darwin-torque-based derivation, though, is general enough to accommodate both a finite inclination and the actual rheology. Although rheologies with Q scaling as the frequency to a positive power make the torque diverge at a zero frequency, this reveals not the impossible nature of the rheology, but a flaw in mathematics, i.e., a common misassumption that damping merely provides lags to the terms of the Fourier series for the tidal potential. A hydrodynamical treatment (Darwin 1879) had demonstrated that the magnitudes of the terms, too, get changed. Reinstating of this detail tames the infinities and rehabilitates the "impossible" scaling law (which happens to be the actual law the terrestrial planets obey at low frequencies).

Spin-orbit coupling can be described in two approaches. The method known as "the MacDonald torque" is often combined with an assumption that the quality factor Q is frequency-independent. This makes the method inconsistent, because the MacDonald theory tacitly fixes the rheology by making Q scale as the inverse tidal frequency. Spin-orbit coupling can be treated also in an approach called "the Darwin torque". While this theory is general enough to accommodate an arbitrary frequency-dependence of Q, this advantage has not yet been exploited in the literature, where Q is assumed constant or is set to scale as inverse tidal frequency, the latter assertion making the Darwin torque equivalent to a corrected version of the MacDonald torque. However neither a constant nor an inverse-frequency Q reflect the properties of realistic mantles and crusts, because the actual frequency-dependence is more complex. Hence the necessity to enrich the theory of spin-orbit interaction with the right frequency-dependence. We accomplish this programme for the Darwin-torque-based model near resonances. We derive the frequency-dependence of the tidal torque from the first principles, i.e., from the expression for the mantle's compliance in the time domain. We also explain that the tidal torque includes not only the secular part, but also an oscillating part. We demonstrate that the lmpq term of the Darwin-Kaula expansion for the tidal torque smoothly goes through zero, when the secondary traverses the lmpq resonance (e.g., the principal tidal torque smoothly goes through nil as the secondary crosses the synchronous orbit). We also offer a possible explanation for the unexpected frequency-dependence of the tidal dissipation rate in the Moon, discovered by LLR.

Mathematical modeling of bodily tides can be carried out in various ways. Most straightforward is the method of complex amplitudes, which is often used in the planetary science. Another method, employed both in planetary science and astrophysics, is based on decomposition of each harmonic of the tide into two bulges oriented orthogonally to one another. We prove that the two methods are equivalent. Specifically, we demonstrate that the two-bulge method is not a separate approximation, but ensues directly from the Fourier expansion of a linear tidal theory equipped with an arbitrary rheological model involving a departure from elasticity. To this end, we use the most general mathematical formalism applicable to linear bodily tides. To express the tidal amendment to the potential of the perturbed primary, we act on the tide-raising potential of the perturbing secondary with a convolution operator. This enables us to interconnect a complex Fourier component of the tidally generated potential of the perturbed primary with the appropriate complex Fourier component of the tide-raising potential of the secondary. Then we demonstrate how this interrelation entails the two-bulge description. While less economical mathematically, the two-bulge approach has a good illustrative power, and may be employed on a par with a more concise method of complex amplitudes. At the same time, there exist situations where the two-bulge method becomes more practical for technical calculations.

In geophysics and seismology, it is a common knowledge that the quality factors Q of the mantle and crust materials scale as the tidal frequency to a positive fractional power (Karato 2007, Efroimsky and Lainey 2007). In astronomy, there exists an equally common belief that such rheological models introduce discontinuities into the equations and thus are unrealistic at low frequencies. We demonstrate that, while such models indeed make the conventional expressions for the tidal torque diverge for vanishing frequencies, the emerging infinities reveal not the impossible nature of one or another rheology, but a subtle flaw in the underlying mathematical model of friction. Flawed is the common misassumption that the tidal force and torque are inversely proportional to the quality factor. In reality, they are proportional to the sine of the tidal phase lag, while the inverse quality factor is commonly identified with the tangent of the lag. The sine and tangent of the lag are close everywhere {\it{except in the vicinity of the zero frequency}}. Reinstating of this detail tames the fake infinities and rehabilitates the "impossible" scaling law (which happens to be the actual law the mantles obey). This preprint is a pilot paper. A more comprehensive treatise on tidal torques is to be published (Efroimsky and Williams 2009).

It was believed until very recently that a near-equatorial satellite would always keep up with the planet's equator (with oscillations in inclination, but without a secular drift). As explained in Efroimsky and Goldreich (2004), this opinion originated from a wrong interpretation of a (mathematically correct) result obtained in terms of non-osculating orbital elements. A similar analysis carried out in the language of osculating elements will endow the planetary equations with some extra terms caused by the planet's obliquity change. Some of these terms will be nontrivial, in that they will not be amendments to the disturbing function. Due to the extra terms, the variations of a planet's obliquity may cause a secular drift of its satellite orbit inclination. In this article we set out the analytical formalism for our study of this drift. We demonstrate that, in the case of uniform precession, the drift will be extremely slow, because the first-order terms responsible for the drift will be short-period and, thus, will have vanishing orbital averages (as anticipated 40 years ago by Peter Goldreich), while the secular terms will be of the second order only. However, it turns out that variations of the planetary precession make the first-order terms secular. For example, the planetary nutations will resonate with the satellite's orbital frequency and, thereby, may instigate a secular drift. A detailed study of this process will be offered in the subsequent publication, while here we work out the required mathematical formalism and point out the key aspects of the dynamics.

Given the impressive investment by the nation in observational Astronomy and Astrophysics facilities coming on line now and in the near future, we advocate for an increased investment in applied and fundamental research on Astrophysical and Celestial Dynamics (ACD). Specifically we call for a) continued and expanded support for applied research in ACD, b) creation of support for fundamental research in ACD and its subfields, and c) the creation of a unified program to help scientists coordinate and collaborate in their research in these fields. The benefits of this proposal are threefold. First, it will enable researchers to interpret and understand the implications of newly observed phenomena that will invariably arise from new facilities and surveys. Second, research on fundamentals will foster connections between specialists, leveraging advances found in one sub-field and making them available to others. Third, a coordinated approach for applied and fundamental research in ACD will help academic institutions in the United States to produce future researchers trained and knowledgeable in essential subfields such as Mathematical Celestial Mechanics and able to continue its advancement in conjunction with the increase in observations.

We revisit the Lagrange and Delaunay systems of equations for the orbital elements, and point out a previously neglected aspect of these equations: in both cases the orbit resides on a certain 9-dimensional submanifold of the 12-dimensional space spanned by the orbital elements and their time derivatives. We demonstrate that there exists a vast freedom in choosing this submanifold. This freedom of choice (=freedom of gauge fixing) reveals a symmetry hiding behind Lagrange's and Delaunay's systems, which is, mathematically, analogous to the gauge invariance in electrodynamics. Just like a convenient choice of gauge simplifies calculations in electrodynamics, so the freedom of choice of the submanifold may, potentially, be used to create simpler schemes of orbit integration. On the other hand, the presence of this feature may be a previously unrecognised source of numerical instability.

Any model of tides is based on a specific hypothesis of how lagging depends on the tidal-flexure frequency. For example, Gerstenkorn (1955), MacDonald (1964), and Kaula (1964) assumed constancy of the geometric lag angle, while Singer (1968) and Mignard (1979, 1980) asserted constancy of the time lag. Thus, each of these two models was based on a certain law of scaling of the geometric lag. The actual dependence of the geometric lag on the frequency is more complicated and is determined by the rheology of the planet. Besides, each particular functional form of this dependence will unambiguously fix the appropriate form of the frequency dependence of the tidal quality factor, Q. Since at present we know the shape of the dependence of Q upon the frequency, we can reverse our line of reasoning and single out the appropriate actual frequency-dependence of the angular lag. This dependence turns out to be different from those employed hitherto, and it entails considerable alterations in the time scales of the tide-generated dynamical evolution. Phobos' fall on Mars is an example we consider.

Construction of a theory of orbits about a precessing oblate planet, in terms of osculating elements defined in a frame of the equator of date, was started in Efroimsky and Goldreich (2004) and Efroimsky (2005, 2006). We now combine that analytical machinery with numerics. The resulting semianalytical theory is then applied to Deimos over long time scales. In parallel, we carry out a purely numerical integration in an inertial Cartesian frame. The results agree to within a small margin, for over 10 Myr, demonstrating the applicability of our semianalytical model over long timescales. This will enable us to employ it at the further steps of the project, enriching the model with the tides, the pull of the Sun, and the planet's triaxiality. Another goal of our work was to check if the equinoctial precession predicted for a rigid Mars could have been sufficient to repel the orbits away from the equator. We show that for low initial inclinations, the orbit inclination reckoned from the precessing equator of date is subject only to small variations. This is an extension, to non-uniform precession given by the Colombo model, of an old result obtained by Goldreich (1965) for the case of uniform precession and a low initial inclination. However, near-polar initial inclinations may exhibit considerable variations for up to +/- 10 deg in magnitude. Nevertheless, the analysis confirms that an oblate planet can, indeed, afford large variations of the equinoctial precession over hundreds of millions of years, without repelling its near-equatorial satellites away from the equator of date: the satellite inclination oscillates but does not show a secular increase. Nor does it show secular decrease, a fact that is relevant to the discussion of the possibility of high-inclination capture of Phobos and Deimos.

The Hamiltonian theory of Earth rotation, known as the Kinoshita-Souchay theory, operates with nonosculating Andoyer elements. This situation parallels a similar phenomenon that often happens (but seldom gets noticed) in orbital dynamics, when the standard Lagrange-type or Delaunay-type planetary equations unexpectedly render nonosculating orbital elements. In orbital mechanics, osculation loss happens when a velocity-dependent perturbation is plugged into the standard planetary equations. In attitude mechanics, osculation is lost when an angular-velocity-dependent disturbance is plugged in the standard dynamical equations for the Andoyer elements. We encounter exactly this situation in the theory of Earth rotation, because this theory contains an angular-velocity-dependent perturbation (the switch from an inertial frame to that associated with the precessing ecliptic of date). While the osculation loss does not influence the predictions for the figure axis of the planet, it considerably alters the predictions for the instantaneous spin-axis' orientation. We explore this issue in great detail.

This paper reviews the Serret-Andoyer (SA) canonical formalism in rigid-body dynamics and presents some new results. As is well known, the problem of unsupported and unperturbed rigid rotator can be reduced. The availability of this reduction is offered by the underlying symmetry, which stems from conservation of the angular momentum and rotational kinetic energy. When a perturbation is turned on, these quantities are no longer preserved. Nonetheless, the language of reduced description remains extremely instrumental even in the perturbed case. We describe the canonical reduction performed by the Serret-Andoyer (SA) method, and discuss its applications to attitude dynamics and to the theory of planetary rotation. Specifically, we consider the case of angular-velocity-dependent torques, and discuss the variation-of-parameters-inherent antinomy between canonicity and osculation. Finally, we address the transformation of the Andoyer variables into the action-angle ones, using the method of Sadov.

We have designed and constructed a ``dispersed Fourier Transform Spectrometer'' (dFTS), consisting of a conventional FTS followed by a grating spectrometer. By combining these two devices, we negate a substantial fraction of the sensitivity disadvantage of a conventional FTS for high resolution, broadband, optical spectroscopy, while preserving many of the advantages inherent to interferometric spectrometers. In addition, we have implemented a simple and inexpensive laser metrology system, which enables very precise calibration of the interferometer wavelength scale. The fusion of interferometric and dispersive technologies with a laser metrology system yields an instrument well-suited to stellar spectroscopy, velocimetry, and extrasolar planet detection, which is competitive with existing high-resolution, high accuracy stellar spectrometers. In this paper, we describe the design of our prototype dFTS, explain the algorithm we use to efficiently reconstruct a broadband spectrum from a sequence of narrowband interferograms, and present initial observations and resulting velocimetry of stellar targets.

We continue the study undertaken in Efroimsky (2005a) where we explored the influence of spin-axis variations of an oblate planet on satellite orbits. Near-equatorial satellites had long been believed to keep up with the oblate primary's equator in the cause of its spin-axis variations. As demonstrated by Efroimsky and Goldreich (2004), this opinion had stemmed from an inexact interpretation of a correct result by Goldreich (1965). Though Goldreich (1965) mentioned that his result (preservation of the initial inclination, up to small oscillations about the moving equatorial plane) was obtained for non-osculating inclination, his admonition has been persistently ignored for forty years. It was explained in Efroimsky and Goldreich (2004) that the equator precession influences the osculating inclination of a satellite orbit already in the first order over the perturbation caused by a transition from an inertial to an equatorial coordinate system. It was later shown in Efroimsky (2005a) that the secular part of the inclination is affected only in the second order. This fact, anticipated by Goldreich (1965), remains valid for a constant rate of the precession. It turns out that non-uniform variations of the planetary spin state generate changes in the osculating elements, that are linear in the planetary equator's total precession rate, rate that includes the equinoctial precession, nutation, the Chandler wobble, and the polar wander. We work out a formalism which will help us to determine if these factors cause a drift of a satellite orbit away from the evolving planetary equator.

Both orbital and rotational dynamics employ the method of variation of parameters. We express, in a non-perturbed setting, the coordinates (Cartesian, in the orbital case, or Eulerian in the rotation case) via the time and six adjustable constants called elements (orbital elements or rotational elements). If, under disturbance, we use this expression as ansatz and endow the "constants" with time dependence, then the perturbed velocity (Cartesian or angular) will consist of a partial derivative with respect to time and a so-called convective term, one that includes the time derivatives of the variable "constants." Out of sheer convenience, the so-called Lagrange constraint is often imposed. It nullifies the convective term and, thereby, guarantees that the functional dependence of the velocity upon the time and "constants" stays, under perturbation, the same as it used to be in the undisturbed setting. When the dynamical equations, written in terms of the "constants," are demanded to be symplectic (and the "constants" make conjugated pairs $ Q, P$), these "constants" are called Delaunay elements, in the orbital case, or Andoyer elements, in the rotational case. The Andoyer and Delaunay sets of elements share a feature not visible with a naked eye: in certain cases, the standard equations render these elements non-osculating. Hence, even though the Andoyer variables in the Kinoshita-Souchay theory are introduced in a precessing frame of the Earth orbit, they nevertheless return the angular velocity relative to an inertial frame.

In orbital and attitude dynamics the coordinates and the Euler angles are expressed as functions of the time and six constants called elements. Under disturbance, the constants are endowed with time dependence. The Lagrange constraint is then imposed to guarantee that the functional dependence of the perturbed velocity on the time and constants stays the same as in the undisturbed case. Constants obeying this condition are called osculating elements. The constants chosen to be canonical are called Delaunay elements, in the orbital case, or Andoyer elements, in the spin case. (As some Andoyer elements are time dependent even in the free-spin case, the role of constants is played by their initial values.) The Andoyer and Delaunay sets of elements share a feature not readily apparent: in certain cases the standard equations render them non-osculating. In orbital mechanics, elements furnished by the standard planetary equations are non-osculating when perturbations depend on velocities. To preserve osculation, the equations must be amended with extra terms that are not parts of the disturbing function. In the case of Delaunay parameterisation, these terms destroy canonicity. So under velocity-dependent disturbances, osculation and canonicity are incompatible. (Efroimsky and Goldreich 2003, 2004) Similarly, the Andoyer elements turn out to be non-osculating under angular-velocity-dependent perturbation. Amendment of only the Hamiltonian makes the equations render nonosculating elements. To make them osculating, more terms must enter the equations (and the equations will no longer be canonical). In practical calculations, is often convenient to deliberately deviate from osculation by substituting the Lagrange constraint with a condition that gives birth to a family of nonosculating elements.

In most books the Delaunay and Lagrange equations for the orbital elements are derived by the Hamilton-Jacobi method: one begins with the 2-body Hamilton equations, performs a canonical transformation to the orbital elements, and obtains the Delaunay system. A standard trick is then used to generalise the approach to the N-body case. We re-examine this step and demonstrate that it contains an implicit condition which restricts the dynamics to a 9(N-1)-dimensional submanifold of the 12(N-1)-dimensional space spanned by the elements and their time derivatives. The tacit condition is equivalent to the so-called Lagrange constraint. It is the condition of the orbital elements being osculating, i.e., of the instantaneous ellipse (or hyperbola) being always tangential to the physical velocity. Imposure of any condition different from the Lagrange constraint (but compatible with the equations of motion) is legitimate and will not alter the physical trajectory or velocity (though will alter the mathematical form of the planetary equations). This freedom of nomination of the supplementary constraint reveals a gauge-type internal symmetry of the celestial-mechanics equations and has consequences for the stability of numerical integrators. Another important aspect of this freedom is that any gauge different from that of Lagrange makes the Delaunay system non-canonical. In a more general setting, when the disturbance depends not only upon positions but also upon velocities, there is a "generalised Lagrange gauge" wherein the Delaunay system is symplectic. This gauge renders orbital elements that are osculating in the phase space. It coincides with the regular Lagrange gauge when the perturbation is velocity-independent.

We summarise research reported in (Efroimsky 2002, 2003; Efroimsky and Goldreich 2003a,b) and develop its application to planetary equations in non-inertial frames. We provide a practical example illustrating how the gauge formalism considerably simplifies the calculation of satellite motion about an oblate precessing planet.

Neutron stars, asteroids, comets, cosmic-dust granules, spacecraft, as well as whatever other freely spinning body dissipate energy when they rotate about any axis different from principal. We discuss the internal-dissipation-caused relaxation of a freely precessing rotator towards its minimal-energy mode (mode that corresponds to the spin about the maximal-inertia axis). While the body nutates at some rate, the internal stresses and strains within the body oscillate at frequencies both higher and lower than this rate. The internal dissipation takes place mostly the second and higher harmonics. We discuss the application of our findings to asteroids. Regarding the comets, estimates show that the currently available angular resolution of spacecraft-based instruments makes it possible to observe wobble damping within year- or maybe even month-long spans of time. We also discuss cosmic-dust astrophysics; in particular, the role played by precession damping in the dust alignment. We show that this damping provides coupling of the grain's rotational and vibrational degrees of freedom; this entails occasional flipping of dust grains due to thermal fluctuations. During such a flip, grain preserves its angular momentum, but the direction of torques arising from H2 formation reverses. As a result, flipping grain will not rotate fast in spite of the action of uncompensated H2 formation torques. The grains get ``thermally trapped,'' and their alignment is marginal.

We develop a comprehensive quantitative description of the cross-section mechanism discovered several years ago by Lazarian. This is one of the processes that determine grain orientation in clouds of suprathermal cosmic dust. The cross-section mechanism manifests itself when an ensemble of suprathermal paramagnetic granules is placed in a magnetic field and is subject to ultrasonic gas bombardment. The mechanism yields dust alignment whose efficiency depends upon two factors: the geometric shape of the granules, and the angle Phi between the magnetic line and the gas flow. We calculate the quantitative measure of this alignment, and study its dependence upon the said factors. It turns out that, irrelevant of the grain shape, the action of a flux does not lead to alignment if Phi = arccos(1/sqrt{3}).

Whenever a freely spinning body is found in a complex rotational state, this means that either the body is a recent victim of an impact or a tidal interaction, or is a fragment of a recently disrupted progenitor. Another factor (relevant for comets) is outgassing. Due to impacts, tidal forces and outgassing, the asteroidal and cometary precession must be a generic phenomenon: while some rotators are in the state of visible tumbling, a much larger amount of objects must be performing narrow-cone precession not so easily observable from the Earth. The internal dissipation in a freely precessing top leads to relaxation (gradual damping of the precession) and sometimes to spontaneous changes in the rotation axis. Recently developed theory of dissipative precession of a rigid body reveals that this is a highly nonlinear process: while the body is precessing at an angular rate $ \omega$, the precession-caused stresses and strains in the body contain components oscillating at other frequencies. Dependent upon the spin state, those frequencies may be higher or, most remarkably, lower than the precession rate. In many states dissipation at the harmonics is comparable to or even exceeds that at the principal frequency. For this and other reasons, in many spin states the damping of asteroidal and cometary wobble happens faster, by several orders, than believed previously. This makes it possible to measure the precession-damping rate. The narrowing of the precession cone through the period of about a year can be registered by the currently available spacecraft-based observational means. However, in the near-separatrix spin states a precessing rotator can considerably slow down its relaxation.

When a solid body is freely rotating at an angular velocity ${\bf \Omega}$, the ellipsoid of constant angular momentum, in the space $\Omega_1, \Omega_2, \Omega_3$, has poles corresponding to spinning about the minimal-inertia and maximal-inertia axes. The first pole may be considered stable if we neglect the inner dissipation, but becomes unstable if the dissipation is taken into account. This happens because the bodies dissipate energy when they rotate about any axis different from principal. In the case of an oblate symmetrical body, the angular velocity describes a circular cone about the vector of (conserved) angular momentum. In the course of relaxation, the angle of this cone decreases, so that both the angular velocity and the maximal-inertia axis of the body align along the angular momentum. The generic case of an asymmetric body is far more involved. Even the symmetrical prolate body exhibits a sophisticated behaviour, because an infinitesimally small deviation of the body's shape from a rotational symmetry (i.e., a small difference between the largest and second largest moments of inertia) yields libration: the precession trajectory is not a circle but an ellipse. In this article we show that often the most effective internal dissipation takes place at twice the frequency of the body's precession. Applications to precessing asteroids, cosmic-dust alignment, and rotating satellites are discussed.

Asteroids and comets dissipate energy when they rotate about the axis different from the axis of the maximal moment of inertia. We show that the most efficient internal relaxation happens at the double frequency of body's tumbling. Therefore the earlier estimates that ignore double frequency input underestimate the internal relaxation in asteroids and comets. We show that the Earth seismological data may poorly represent acoustic properties of asteroids and comet as internal relaxation increases in the presence of moisture. At the same time owing to non-linearlity of inelastic relaxation small angle nutations can persist for very long time spans, but our ability to detect such precessions is limited by the resolution of the radar-generated images. Wobbling may provide valuable information on the composition and structure of asteroids and on their recent history of external impacts.

A body dissipates energy when it freely rotates about any axis different from principal. This entails relaxation, i.e., decrease of the rotational energy, with the angular momentum preserved. The spin about the major-inertia axis corresponds to the minimal kinetic energy, for a fixed angular momentum. Thence one may expect comets and asteroids (as well as spacecraft or cosmic-dust granules) stay in this, so-called principal, state of rotation, unless they are forced out of this state by a collision, or a tidal interaction, or cometary jetting, or by whatever other reason. As is well known, comet P/Halley, asteroid 4179 Toutatis, and some other small bodies exhibit very complex rotational motions attributed to these objects being in non-principal states of spin. Most probably, the asteroid and cometary wobble is quite a generic phenomenon. The theory of wobble with internal dissipation has not been fully developed as yet. In this article we demonstrate that in some spin states the effectiveness of the inelastic-dissipation process is several orders of magnitude higher than believed previously, and can be measured, by the presently available observational instruments, within approximately a year span. We also show that in some other spin states both the precession and precession-relaxation processes slow down considerably. (We call it near-separatrix lingering effect.) Such spin states may evolve so slowly that they can mimic the principal-rotation state.

Motivated by a recent study by Lazarian and Draine, which showed that a high degree of grain alignment of the paramagnetic dust is achievable if the rates of internal relaxation are controlled by the Barnett relaxation process, we undertake a study of an alternative mechanism of internal dissipation, namely, the inelastic dissipation of energy in oblate dust grains. We find that deformations at double frequency that were disregarded in earlier studies dominate the inelastic relaxation. Our results indicate that for grains with 4:1 axis ratio, or/and grains formed via agglomeration inelastic relaxation dominates the Barnett relaxation within large (a>0.1 micron) grains even when they rotate at thermal rate. For grains with axis ratio less than 1:2 the inelastic relaxation is dominant when grains rotate suprathermally.

This paper provides a quantitative account of a recently introduced mechanism of mechanical alignment of suprathermally rotating grains. These rapidly rotating grains are essentially not susceptible to random torques arising from gas-grain collisions, as the timescales for such torques to have significant effect are orders of magnitude greater than the mean time between crossovers. Such grains can be aligned by gaseous torques during the short periods of crossovers and/or due to the difference in the rate at which atoms arrive at grain surface. The latter is a result of the difference in orientation of a grain in respect to the supersonic flow. This process, which we call cross-section alignment, is the subject of our present paper. We derive expressions for the measure of cross-section alignment for oblate grains and study how this measure depends upon the angle between the interstellar magnetic field and the gaseous flow and upon the grain shape.