Recursive computation of mutual potential between two polyhedra

Recursive computation of mutual potential, force, and torque between two polyhedra is studied. Based on formulations by Werner and Scheeres (Celest Mech Dyn Astron 91:337–349, 2005) and Fahnestock and Scheeres (Celest Mech Dyn Astron 96:317–339, 2006) who applied the Legendre polynomial expansion to gravity interactions and expressed each order term by a shape-dependent part and a shape-independent part, this paper generalizes the computation of each order term, giving recursive relations of the shape-dependent part. To consider the potential, force, and torque, we introduce three tensors. This method is applicable to any multi-body systems. Finally, we implement this recursive computation to simulate the dynamics of a two rigid-body system that consists of two equal-sized parallelepipeds.     

The two-body interaction potential in the STF tensor formalism: an application to binary asteroids

The symmetric trace free (STF) tensor formalism, developed by Hartmann et al. (Celest Mech Dyn Astron 60:139–159. doi:10.1007/BF00693097, 1994), is a nice tool, not much used in Celestial Mechanics. It is fully equivalent to the usual spherical harmonics but permits more elegant and compact formulations. The coupling between the gravitational fields of extended bodies with this formalism has been used in Mathis and Le Poncin-Lafitte (Astron Astrophys 497:889–910. doi:10.1051/0004-6361/20079054, 2009) for binary stars or planetary systems, but not yet applied to binary asteroids. However, binary asteroids are common in the Solar System and usually their study requires a full two rigid body approach. The formulation of the two-body interaction potential in the STF formalism in the full two rigid body problem is detailed and completed in this article. An application to the binary asteroid (66391) 1999 KW4 is presented with a comparison of our results with other results of the literature for validation.     

Erratum to: Reports of the IAU Working Group on Cartographic Coordinates and Rotational Elements: 2006 & 2009

The primary poles for (243) Ida and (134340) Pluto and its satellite (134340) Pluto : I Charon were redefined in the IAU Working Group on Cartographic Coordinates and Rotational Elements (WGCCRE) 2006 report (Seidelmann et al. in Celest Mech Dyn Astr 98:155, 2007), and 2009 report (Archinal et al. in Celest Mech Dyn Astr 109:101, 2011), respectively, to be consistent with the primary poles of similar Solar System bodies. However, the WGCCRE failed to take into account the effect of the redefinition of the poles on the values of the rotation angle W at J2000.0. The revised relationships in Table 3 of Archinal et al. 2011) are $$\begin{array}{llll} W & = & 274^{\circ}.05 +1864^{\circ}.6280070\, d\;{\rm for\; (243)\,Ida} \\ W & = & 302^{\circ} .695 + 56^{\circ} .3625225\, d\;{\rm for\; (134340)\,Pluto,\; and}\\ W & = & 122^{\circ} .695 + 56^{\circ} .3625225\, d\;{\rm for\; (134340)\,Pluto : I \,Charon}\end{array}$$ where d is the time in TDB days from J2000.0 (JD2451545.0).     

Increase of the mean Sun–Earth distance caused by a secular mass accumulation

Based on many planetary observations between the years 1971 and 2003, Krasinsky and Brumberg (Celest. Mech. Dyn. Astron. 90:267–288, 2004) have estimated a rate of increase in the mean Sun-Earth distance of (15±4) m per century. Together with other anomalous observations in the solar system, this increase appears to be unexplained (Lämmerzahl et al. in Astrophys. Space Sci. Lib., vol. 349, pp. 75–101, 2008). We explain these findings by invoking a recently proposed gravitational impact model (Wilhelm et al. in Astrophys. Space Sci. 343:135–144, 2013) that implies a secular mass increase of all massive bodies. This allows us to formulate a quantitative understanding of the effect within the parameter range of the model with a mass accumulation rate of the Sun of (6.4±1.7)×10¹⁰ kg?s^?1.     

Position and velocity sensitivities at the triangular libration points in the restricted problem of three bodies when the bigger primary is an oblate body

In this paper we have examined the stability of triangular libration points in the restricted problem of three bodies when the bigger primary is an oblate spheroid. Here we followed the time limit and computational process of Tuckness (Celest. Mech. Dyn. Mech. 61, 1–19, 1995) on the stability criteria given by McKenzie and Szebehely (Celest. Mech. 23, 223–229, 1981). In this study it was found that in comparison to other studies the value of the critical mass μ _c has been reduced due to oblateness of the bigger primary, i.e. the range of stability of the equilateral triangular libration points reduced with the increase of the oblateness parameter I and hence the order of commensurability was increased.     

A new set of integrals of motion to propagate the perturbed two-body problem

A formulation of the perturbed two-body problem that relies on a new set of orbital elements is presented. The proposed method represents a generalization of the special perturbation method published by Peláez et al. (Celest Mech Dyn Astron 97(2):131–150, 2007) for the case of a perturbing force that is partially or totally derivable from a potential. We accomplish this result by employing a generalized Sundman time transformation in the framework of the projective decomposition, which is a known approach for transforming the two-body problem into a set of linear and regular differential equations of motion. Numerical tests, carried out with examples extensively used in the literature, show the remarkable improvement of the performance of the new method for different kinds of perturbations and eccentricities. In particular, one notable result is that the quadratic dependence of the position error on the time-like argument exhibited by Peláez’s method for near-circular motion under the $J_{2}$ perturbation is transformed into linear. Moreover, the method reveals to be competitive with two very popular element methods derived from the Kustaanheimo-Stiefel and Sperling-Burdet regularizations.     

On the convergence of an algorithm constructing the normal form for elliptic lower dimensional tori in planetary systems

We give a constructive proof of the existence of elliptic lower dimensional tori in nearly integrable Hamiltonian systems. In particular we adapt the classical Kolmogorov normalization algorithm to the case of planetary systems, for which elliptic tori may be used as replacements of elliptic Keplerian orbits in Lagrange-Laplace theory. With this paper we support with rigorous convergence estimates the semi-analytic work in our previous article (Sansottera et al., Celest Mech Dyn Astron 111:337–361, 2011), where an explicit calculation of an invariant torus for a planar model of the Sun-Jupiter-Saturn-Uranus system has been made. With respect to previous works on the same subject we exploit the characteristic of Lie series giving a precise control of all terms generated by our algorithm. This allows us to slightly relax the non-resonance conditions on the frequencies.     

Very Fast and Accurate Azimuth Disambiguation of Vector Magnetograms

We present a method for fast and accurate azimuth disambiguation of vector magnetogram data regardless of the location of the analyzed region on the solar disk. The direction of the transverse field is determined with the principle of minimum deviation of the field from the reference (potential) field. The new disambiguation (NDA) code is examined on the well-known models of Metcalf et al. (Solar Phys. 237, 267, 2006) and Leka et al. (Solar Phys. 260, 83, 2009), and on an artificial model based on the observed magnetic field of AR 10930 (Rudenko, Myshyakov, and Anfinogentov, Astron. Rep. 57, 622, 2013). We compare Hinode/SOT-SP vector magnetograms of AR 10930 disambiguated with three codes: the NDA code, the nonpotential magnetic-field calculation (NPFC: Georgoulis, Astrophys. J. Lett. 629, L69, 2005), and the spherical minimum-energy method (Rudenko, Myshyakov, and Anfinogentov, Astron. Rep. 57, 622, 2013). We then illustrate the performance of NDA on SDO/HMI full-disk magnetic-field observations. We show that our new algorithm is more than four times faster than the fastest algorithm that provides the disambiguation with a satisfactory accuracy (NPFC). At the same time, its accuracy is similar to that of the minimum-energy method (a very slow algorithm). In contrast to other codes, the NDA code maintains high accuracy when the region to be analyzed is very close to the limb.     

Dynamics of Kepler problem with linear drag

We study the dynamics of Kepler problem with linear drag. We prove that motions with nonzero angular momentum have no collisions and travel from infinity to the singularity. In the process, the energy takes all real values and the angular velocity becomes unbounded. We also prove that there are two types of linear motions: capture–collision and ejection–collision. The behaviour of solutions at collisions is the same as in the conservative case. Proofs are obtained using the geometric theory of ordinary differential equations and two regularizations for the singularity of Kepler problem equation. The first, already considered in Diacu (Celest Mech Dyn Astron 75:1–15, 1999), is mainly used for the study of the linear motions. The second, the well known Levi-Civita transformation, allows to complete the study of the asymptotic values of the energy and to prove the existence of collision solutions with arbitrary energy.     

Inclusion of sdBs in evolutionary population synthesis for binary stellar populations and the application: the determination of photo-z and galaxy morphology

Subdwarf B stars (sdBs) can significantly change the ultraviolet spectra of populations at age t～1 Gyr, and have been even included in the evolutionary population synthesis (EPS) models by Han et al. (Mon. Not. R. Astron. Soc. 380:1098, 2007). In this study we present the spectral energy distributions (SEDs) of binary stellar populations (BSPs) by combining the EPS models of Han et al. (Mon. Not. R. Astron. Soc. 380:1098, 2007) and those of the Yunnan group (Zhang et al. in Astron. Astrophys. 415:117, 2004; Mon. Not. R. Astron. Soc. 357:1088, 2005), which have included various binary interactions (except sdBs) in EPS models. This set of SEDs is available upon request from the authors. Using this set of SEDs of BSPs we build the spectra of Burst, E, S0–Sd and Irr types of galaxies by using the package of Bruzual and Charlot (Mon. Not. R. Astron. Soc. 344:1000, 2003). Combined with the photometric data (filters and magnitudes), we obtain the photometric redshifts and morphologies of 1502 galaxies by using the Hyperz code of Bolzonella et al. (Astron. Astrophys. 363:476, 2000). This sample of galaxies is obtained by removing those objects, mismatched with the SDSS/DR7 and GALEX/DR4, from the catalogue of Fukugita et al. (Astron. J. 134:579, 2007). By comparison the results with the SDSS spectroscopic redshifts and the morphological index of Fukugita et al. (Astron. J. 134:579, 2007), we find that the photo-z fluctuate with the SDSS spectroscopic redshifts, while the Sa–Sc galaxies in the catalogue of Fukugita et al. (Astron. J. 134:579, 2007) are classified earlier as Burst-E galaxies.     

Secondary resonances and the boundary of effective stability of Trojan motions

One of the most interesting features in the libration domain of co-orbital motions is the existence of secondary resonances. For some combinations of physical parameters, these resonances occupy a large fraction of the domain of stability and rule the dynamics within the stable tadpole region. In this work, we present an application of a recently introduced ‘basic Hamiltonian model’ \(H_\mathrm{b}\) for Trojan dynamics (Páez and Efthymiopoulos in Celest Mech Dyn Astron 121(2):139, 2015; Páez et al. in Celest Mech Dyn Astron 126:519, 2016): we show that the inner border of the secondary resonance of lowermost order, as defined by \(H_\mathrm{b}\), provides a good estimation of the region in phase space for which the orbits remain regular regardless of the orbital parameters of the system. The computation of this boundary is straightforward by combining a resonant normal form calculation in conjunction with an ‘asymmetric expansion’ of the Hamiltonian around the libration points, which speeds up convergence. Applications to the determination of the effective stability domain for exoplanetary Trojans (planet-sized objects or asteroids) which may accompany giant exoplanets are discussed.     

Spectral behavior of AM Her and QQ Vul in high and intermediate states in the UV

We present low resolution UV spectra of two polar systems, AM Her and QQ Vul from the observations taken by the IUE (International Ultraviolet Explorer) of the period between 1978–1996 and 1983–1996 for both systems respectively, to accomplish a large scale study of what happens to the ultraviolet flux of C IV 1550 Å spectral line during different orbital phases. Two spectra for both systems showing the variations in line fluxes and line widths at different orbital phases in high and intermediate states are presented. We concentrated on calculating the line fluxes and line widths of C IV 1550 Å emission line originating in the accretion stream. Our results show that there is spectral variability for the aformentioned physical parameters at different times, similar to that known for the light curve (Heise and Verbunt, Astron. Astrophys. 189:112, 1988; Gansicke et al., Astron. Astrophys. 303:127, 1995; Kafka and Honeycutt, Astron. J. 125:2188K, 2003). We attribute it to the variations of both density and temperature as a result of changing the mass transfer rate (Hutchings et al., Astron. J. 123:2841H, 2002; King and Lasota, Astron. Astrophys. 140L:16K, 1984) which is responsible for this spectral variability. Also we found that the line fluxes of AM Her are greater than the line fluxes of QQ Vul, while the line widths of both systems are approximately the same.     

Viscoelastic tides: models for use in Celestial Mechanics

In this note, by using Smale’s \(\alpha \)-theorem on the convergence of Newton’s method, the \(\alpha \)-sets of convergence of some starters of solving the elliptic Kepler’s equation are derived. For each starter we compute the exact \(\alpha \)-set in the eccentricity-main anomaly \((e,M)\in [0,1)\times [0,\pi ]\), showing that these sets are larger than those derived by Avendaño et al. (Celest Mech Dyn Astron 119:27–44, 2014). Further, new convergence tests based on the Newton–Kantorowitch theorem are given comparing with the derived from Smale’s \(\alpha \)-test.     

Long-Term Variations of Solar Differential Rotation and Sunspot Activity: Revisited

Long-term variations of solar differential rotation and sunspot activity are investigated through re-analyzing the data on parameters of the differential-rotation law obtained by Makarov, Tlatov, and Callebaut (Solar Phys. 170, 373, 1997), Javaraiah, Bertello, and Ulrich (Astrophys. J. 626, 579, 2005a; Solar Phys. 232, 25, 2005b), and Javaraiah et al. (Solar Phys. 257, 61, 2009). Our results indicate that the solar-surface-rotation rate at the Equator (indicated by the A-parameter of the standard solar-rotation law) shows a secular decrease since Cycle 12 onwards, given by about 1?–?1.5×10^?3 (deg?day^?1?year^?1). The B-parameter of the standard differential-rotation law seems to also show a secular decrease since Cycle 12 onwards, but of weak statistical significance. The rotation rate averaged over latitudes 0^°?–?40^° does not show a secular trend of statistical significance. Moreover, the average sunspot area shows a secular increase of statistical significance since Cycle 12 onwards, while a negative correlation is found between the level of sunspot activity (indicated by the average sunspot area) and the solar equatorial rotation on long-term scales.     

The effects of deformation inertia (kinetic energy) in the orbital and spin evolution of close-in bodies

The purpose of this work is to evaluate the effect of deformation inertia on tide dynamics, particularly within the context of the tide response equations proposed independently by Boué et al. (Celest Mech Dyn Astron 126:31–60, 2016) and Ragazzo and Ruiz (Celest Mech Dyn Astron 128(1):19–59, 2017). The singular limit as the inertia tends to zero is analyzed, and equations for the small inertia regime are proposed. The analysis of Love numbers shows that, independently of the rheology, deformation inertia can be neglected if the tide-forcing frequency is much smaller than the frequency of small oscillations of an ideal body made of a perfect (inviscid) fluid with the same inertial and gravitational properties of the original body. Finally, numerical integration of the full set of equations, which couples tide, spin and orbit, is used to evaluate the effect of inertia on the overall motion. The results are consistent with those obtained from the Love number analysis. The conclusion is that, from the point of view of orbital evolution of celestial bodies, deformation inertia can be safely neglected. (Exceptions may occur when a higher-order harmonic of the tide forcing has a high amplitude.)     

On the Origin of Pulsations of Sub-THz Emission from Solar Flares

We propose a model to explain fast pulsations in sub-THz emission from solar flares. The model is based on the approach of a flaring loop as an equivalent electric circuit and explains the pulse-repetition rate, the high-quality factor, Q≥10³, low modulation depth, pulse synchronism at different frequencies, and the dependence of the pulse-repetition rate on the emission flux, observed by Kaufmann et al. (Astrophys. J. 697, 420, 2009). We solved the nonlinear equation for electric current oscillations using a Van der Pol method and found the steady-state value for the amplitude of the current oscillations. Using the pulse rate variation during the flare on 4 November 2003, we found a decrease of the electric current from 1.7×10¹² A in the flare maximum to 4×10¹⁰ A just after the burst. Our model is consistent with the plasma mechanism of sub-THz emission suggested recently by Zaitsev, Stepanov, and Melnikov (Astron. Lett. 39, 650, 2013).     

Modified Coronal Index of the Solar Activity

The original coronal index of the solar activity (CI) has been constructed on the basis of ground-based measurements of the intensities of the coronal line of 530.3 nm (Rybanský in Bull. Astron. Inst. Czechoslov., 28, 367, 1975; Rybanský et al. in J. Geophys. Res., 110, A08106, 2005). In this paper, CI is compared with the EUV measurements on the CELIAS/SEM equipment based on the same idea as the original idea of the coronal index. The correlation is very good for the period 1996?–?2005 (r=0.94 for daily values). The principal result of this paper is the introduction of the modified coronal index (MCI) which in all uses and contexts can replace the existing CI index. Daily MCI values extend over a time period of six solar activity cycles. Future MCI measurements will be derived from more reliable measurements made by space-based observatories that are not influenced by the weather. MCI measurements are and will continue to be archived at the web site of the Slovak Central Observatory in Hurbanovo ( http://www.suh.sk/obs/vysl/MCI.htm ).     

Effects of Random Flows on the Solar f Mode: II. Horizontal and Vertical Flow

We study the influence of horizontal and vertical random flows on the solar f mode in a plane-parallel, incompressible model that includes a static atmosphere. The incompressible limit is an adequate approximation for f-mode type of surface waves that are highly incompressible. The paper revisits and extends the problem investigated earlier by Murawski and Roberts (Astron. Astrophys. 272, 601, 1993). We show that the consideration of the proposed velocity profile requires several restrictive assumptions to be made. These constraints were not recognised in previous studies. The impact of the inconsistencies in earlier modelling is analysed in detail. Corrections to the dispersion relation are derived and the relevance of these corrections is analysed. Finally, the importance of the obtained results is investigated in the context of recent helioseismological data. Detailed comparison with our complementary studies on random horizontal flows (Mole, Kerekes, and Erdélyi, Solar Phys., accepted, 2008) and the random magnetic model of Erdélyi, Kerekes, and Mole (Astron. Astrophys. 431, 1083, 2005) is also given. In particular, for realistic solar parameters we find significant frequency reduction and wave damping, both of which increase with the characteristic thickness of the random layer.     

Implementation and Comparison of Acoustic Travel-Time Measurement Procedures for the Solar Dynamics Observatory/Helioseismic and Magnetic Imager Time?–?Distance Helioseismology Pipeline

The Helioseismic and Magnetic Imager (HMI) instrument onboard the Solar Dynamics Observatory (SDO) satellite is designed to produce high-resolution Doppler-velocity maps of oscillations at the solar surface with high temporal cadence. To take advantage of these high-quality oscillation data, a?time?–?distance helioseismology pipeline (Zhao et al., Solar Phys. submitted, 2010) has been implemented at the Joint Science Operations Center (JSOC) at Stanford University. The aim of this pipeline is to generate maps of acoustic travel times from oscillations on the solar surface, and to infer subsurface 3D flow velocities and sound-speed perturbations. The wave travel times are measured from cross-covariances of the observed solar oscillation signals. For implementation into the pipeline we have investigated three different travel-time definitions developed in time?–?distance helioseismology: a Gabor-wavelet fitting (Kosovichev and Duvall, SCORE’96: Solar Convection and Oscillations and Their Relationship, ASSL, Dordrecht, 241, 1997), a?minimization relative to a reference cross-covariance function (Gizon and Birch, Astrophys. J. 571, 966, 2002), and a linearized version of the minimization method (Gizon and Birch, Astrophys. J. 614, 472, 2004). Using Doppler-velocity data from the Michelson Doppler Imager (MDI) instrument onboard SOHO, we tested and compared these definitions for the mean and difference travel-time perturbations measured from reciprocal signals. Although all three procedures return similar travel times in a quiet-Sun region, the method of Gizon and Birch (Astrophys. J. 614, 472, 2004) gives travel times that are significantly different from the others in a magnetic (active) region. Thus, for the pipeline implementation we chose the procedures of Kosovichev and Duvall (SCORE’96: Solar Convection and Oscillations and Their Relationship, ASSL, Dordrecht, 241, 1997) and Gizon and Birch (Astrophys. J. 571, 966, 2002). We investigated the relationships among these three travel-time definitions, their sensitivities to fitting parameters, and estimated the random errors that they produce.     

Fast numerics for the spin orbit equation with realistic tidal dissipation and constant eccentricity

We present an algorithm for the rapid numerical integration of a time-periodic ODE with a small dissipation term that is \(C^1\) in the velocity. Such an ODE arises as a model of spin–orbit coupling in a star/planet system, and the motivation for devising a fast algorithm for its solution comes from the desire to estimate probability of capture in various solutions, via Monte Carlo simulation: the integration times are very long, since we are interested in phenomena occurring on timescales of the order of \(10^6\)–\(10^7\) years. The proposed algorithm is based on the high-order Euler method which was described in Bartuccelli et al. (Celest Mech Dyn Astron 121(3):233–260, 2015), and it requires computer algebra to set up the code for its implementation. The payoff is an overall increase in speed by a factor of about 7.5 compared to standard numerical methods. Means for accelerating the purely numerical computation are also discussed.