The Sitnikov problem for several primary bodies configurations

In this paper we address an \(n+1\)-body gravitational problem governed by the Newton’s laws, where n primary bodies orbit on a plane \(\varPi \) and an additional massless particle moves on the perpendicular line to \(\varPi \) passing through the center of mass of the primary bodies. We find a condition for the described configuration to be possible. In the case when the primaries are in a rigid motion, we classify all the motions of the massless particle. We study the situation when the massless particle has a periodic motion with the same minimal period as the primary bodies. We show that this fact is related to the existence of a certain pyramidal central configuration.     

Secular tidal changes in lunar orbit and Earth rotation

Small tidal forces in the Earth–Moon system cause detectable changes in the orbit. Tidal energy dissipation causes secular rates in the lunar mean motion n, semimajor axis a, and eccentricity e. Terrestrial dissipation causes most of the tidal change in n and a, but lunar dissipation decreases eccentricity rate. Terrestrial tidal dissipation also slows the rotation of the Earth and increases obliquity. A tidal acceleration model is used for integration of the lunar orbit. Analysis of lunar laser ranging (LLR) data provides two or three terrestrial and two lunar dissipation parameters. Additional parameters come from geophysical knowledge of terrestrial tides. When those parameters are converted to secular rates for orbit elements, one obtains dn/dt = \(-25.97\pm 0.05 ''/\)cent\(^{2}\), da/dt = 38.30 ± 0.08 mm/year, and di/dt = ?0.5 ± 0.1 \(\upmu \)as/year. Solving for two terrestrial time delays and an extra de/dt from unspecified causes gives \(\sim \) \(3\times 10^{-12}\)/year for the latter; solving for three LLR tidal time delays without the extra de/dt gives a larger phase lag of the N2 tide so that total de/dt = \((1.50 \pm 0.10)\times 10^{-11}\)/year. For total dn/dt, there is \(\le \)1 % difference between geophysical models of average tidal dissipation in oceans and solid Earth and LLR results, and most of that difference comes from diurnal tides. The geophysical model predicts that tidal deceleration of Earth rotation is \(-1316 ''\)/cent\(^{2}\) or 87.5 s/cent\(^{2}\) for UT1-AT, a 2.395 ms/cent increase in the length of day, and an obliquity rate of 9 \(\upmu \)as/year. For evolution during past times of slow recession, the eccentricity rate can be negative.     

Laplace series for the level ellipsoid of revolution

The outer gravitational potential V of the level ellipsoid of revolution T is uniquely determined by two quantities: the eccentricity \(\varepsilon \) of the ellipsoid and Clairaut parameter q, proportional to the angular velocity of rotation squared and inversely proportional to the mean density of the ellipsoid. Quantities \(\varepsilon \) and q are independent, though they lie in a rather strict two-dimensional domain. It follows that Stokes coefficients \(I_n\) of Laplace series representing the outer potential of T are uniquely determined by \(\varepsilon \) and q. In this paper, we have found explicit expressions for Stokes coefficients via \(\varepsilon \) and q, as well as their asymptotics when \(n\rightarrow \infty \). If T does not coincide with a Maclaurin ellipsoid, then \(|I_n|\sim B\varepsilon ^n/n\) with a certain constant B. Let us compare this asymptotics with one of \(I_n\) for ellipsoids constrained by the only condition of increasing (even nonstrict) of oblateness from the centre to the periphery: \(|I_n|\sim \bar{B}\varepsilon ^n/(n^2)\). Hence, level ellipsoids with ellipsoidal equidensites do not exist. The only exception represents Maclaurin ellipsoids. It should be recalled that we confine ourselves by ellipsoids of revolution.     

Counting relative equilibrium configurations of the full two-body problem

Consider a system of two rigid, massive bodies interacting according to their mutual gravitational attraction. In a relative equilibrium motion, the bodies rotate rigidly and uniformly about a fixed axis in \({\mathbb {R}}^3\). This is possible only for special positions and orientations of the bodies. After fixing the angular momentum, these relative equilibrium configurations can be characterized as critical points of a smooth function on configuration space. The goal of this paper is to use Morse theory and Lusternik–Schnirelmann category theory to give lower bounds for the number of critical points when the angular momentum is sufficiently large. In addition, the exact number of critical points and their Morse indices are found in the limit as the angular momentum tends to infinity.     

Contractivity-preserving explicit multistep Hermite–Obrechkoff series differential equation solvers

New optimal, contractivity-preserving (CP), explicit, d-derivative, k-step Hermite–Obrechkoff series methods of order p up to \(p=20\), denoted by CP HO(d, k, p), with nonnegative coefficients are constructed. These methods are used to solve nonstiff first-order initial value problems \(y'=f(t,y)\), \(y(t_0)=y_0\). The upper bound \(p_u\) of order p of HO(d, k, p) can reach, approximately, as high as 2.4 times the number of derivatives d. The stability regions of HO(d, k, p) have generally a good shape and grow with decreasing \(p-d\). We, first, note that three selected CP HO methods: 4-derivative 7-step HO of order 13, denoted by HO(4, 7, 13), 5-derivative 6-step HO of order 13, denoted by HO(5, 6, 13), and 9-derivative 2-step HO of order 13, denoted by CMDAHO(13) compare favorably with Adams–Cowell of order 13, denoted by AC(13), in solving standard N-body problems over an interval of 1000 periods on the basis of the relative error of energy as a function of the CPU time. Next, the three HO methods compare positively with AC(13) in solving standard N-body problems on the basis of the growth of relative positional error and relative energy error over 10, 000 periods of integration. Finally, these three methods compare also well with P-stable methods of Cash and Franco et al. on some quasi periodic, second-order linear and nonlinear problems. The coefficients of selected HO methods are listed in the appendix.     

Behaviour of physical parameters in extended gravity with hyperbolic function

In this paper, we have constructed the cosmological model of the universe in f(R, T) theory of gravity in a Bianchi type \(\mathrm{VI}_h\) universe for the functional f(R, T) in the form \(f(R,T)=\mu R+\mu T\), where R and T are respectively Ricci scalar and trace of energy momentum tensor and \(\mu \) is a constant. We have made use of the hyperbolic scale factor to find the physical parameters and metric potentials defined in the space-time. The physical parameters are constrained from different representative values to build up a realistic cosmological model aligned with the observational behaviour. The state finder diagnostic pair is found to be in the acceptable range. The energy conditions of the model are also studied.     

Average Magnetic Field Magnitude Profiles of <Emphasis Type="Italic">Wind</Emphasis> Magnetic Clouds as a Function of Closest Approach to the Clouds’ Axes and Comparison to Model

We examine the average magnetic field magnitude (\(| \boldsymbol{B} | \equiv B\)) within magnetic clouds (MCs) observed by the Wind spacecraft from 1995 to July 2015 to understand the difference between this \(B\) and the ideal \(B\)-profiles expected from using the static, constant-\(\alpha\), force-free, cylindrically symmetric model for MCs of Lepping, Jones, and Burlaga (J. Geophys. Res. 95, 11957, 1990, denoted here as the LJB model). We classify all MCs according to an assigned quality, \(Q_{0}\) (\(= 1, 2, 3\), for excellent, good, and poor). There are a total of 209 MCs and 124 when only \(Q_{0} = 1\), 2 cases are considered. The average normalized field with respect to the closest approach (\(\mathit{CA}\)) is stressed, where we separate cases into four \(\mathit{CA}\) sets centered at 12.5 %, 37.5 %, 62.5 %, and 87.5 % of the average radius; the averaging is done on a percentage-duration basis to treat all cases the same. Normalized \(B\) means that before averaging, the \(B\) for each MC at each point is divided by the LJB model-estimated \(B\) for the MC axis, \(B_{0}\). The actual averages for the 209 and 124 MC sets are compared to the LJB model, after an adjustment for MC expansion (e.g. Lepping et al. in Ann. Geophys. 26, 1919, 2008). This provides four separate difference-relationships, each fitted with a quadratic (Quad) curve of very small \(\sigma\). Interpreting these Quad formulae should provide a comprehensive view of the variation in normalized \(B\) throughout the average MC, where we expect external front and rear compression to be part of its explanation. These formulae are also being considered for modifying the LJB model. This modification will be used in a scheme for forecasting the timing and magnitude of magnetic storms caused by MCs. Extensive testing of the Quad formulae shows that the formulae are quite useful in correcting individual MC \(B\)-profiles, especially for the first \({\approx\,}1/3\) of these MCs. However, the use of this type of \(B\) correction constitutes a (slight) violation of the force-free assumption used in the original LJB MC model.     

Satellite de-orbiting via controlled solar radiation pressure

We analyze the families of central configurations of the spatial 5-body problem with four masses equal to 1 when the fifth mass m varies from 0 to \(+\infty \). In particular we continue numerically, taking m as a parameter, the central configurations (which all are symmetric) of the restricted spatial (\(4+1\))-body problem with four equal masses and \(m=0\) to the spatial 5-body problem with equal masses (i.e. \(m=1\)), and viceversa we continue the symmetric central configurations of the spatial 5-body problem with five equal masses to the restricted (\(4+1\))-body problem with four equal masses. Additionally we continue numerically the symmetric central configurations of the spatial 5-body problem with four equal masses starting with \(m=1\) and ending in \(m=+\infty \), improving the results of Alvarez-Ramírez et al. (Discrete Contin Dyn Syst Ser S 1: 505–518, 2008). We find four bifurcation values of m where the number of central configuration changes. We note that the central configurations of all continued families varying m from 0 to \(+\infty \) are symmetric.     

Ratio of the Core to the Extended Emissions in the Comoving Frame for Blazars

In a two-component jet model, the emissions are the sum of the core and extended emissions: \(S^{\mathrm{ob}}=S_{\mathrm{core}}^{\mathrm{ob}}+S_{\mathrm{ext}}^{\mathrm{ob}}\), with the core emissions, \(S_{\mathrm{core}}^{\mathrm{ob}}= f S_{\mathrm{ext}}^{\mathrm{ob}}\delta ^{q}\) being a function of the Doppler factor \(\delta \), the extended emission \(S_{\mathrm{ext}}^{\mathrm{ob}}\), the jet type dependent factor q, and the ratio of the core to the extended emissions in the comoving frame, f. The f is an unobservable but important parameter. Following our previous work, we collect 65 blazars with available Doppler factor \(\delta \), superluminal velocity \(\beta _{\mathrm{app}}\), and core-dominance parameter, R, and calculated the ratio, f, and performed statistical analyses. We found that the ratio, f, in BL Lacs is on average larger than that in FSRQs. We suggest that the difference of the ratio f between FSRQs and BL Lacs is one of the possible reasons that cause the difference of other observed properties between them. We also find some significant correlations between \(\log f\) and other parameters, including intrinsic (de-beamed) peak frequency, \(\log \nu _{\mathrm{p}}^{\mathrm{in}}\), intrinsic polarization, \(\log P^{\mathrm{in}}\), and core-dominance parameter, \(\log R\), for the whole sample. In addition, we show that the ratio, f, can be estimated by R.     

Dynamics of Large-Scale Solar-Wind Streams Obtained by the Double Superposed Epoch Analysis: 3. Deflection of the Velocity Vector

This work is a continuation of our previous articles (Yermolaev et al. in J. Geophys. Res.120, 7094, 2015 and Yermolaev et al. in Solar Phys.292, 193, 2017), which describe the average temporal profiles of interplanetary plasma and field parameters in large-scale solar-wind (SW) streams: corotating interaction regions (CIRs), interplanetary coronal mass ejections (ICMEs, including both magnetic clouds (MCs) and ejecta), and sheaths as well as interplanetary shocks (ISs). Changes in the longitude angle, \(\varphi\), in CIRs from ?2 to \(2^{\circ}\) agree with earlier results (e.g. Gosling and Pizzo, 1999). We have also analyzed the average temporal profiles of the bulk velocity angles in sheaths and ICMEs. We have found that the angle \(\varphi\) in ICMEs changes from 2 to \(-2^{\circ}\), while in sheaths it changes from ?2 to \(2^{\circ}\) (similar to the change in CIRs), i.e. the angle in CIRs and sheaths deflects in the opposite sense to ICMEs. When averaging the latitude angle \(\vartheta\) on all the intervals of the chosen SW types, the angle \(\vartheta\) is almost constant at \({\sim}\,1^{\circ}\). We made for the first time a selection of SW events with increasing and decreasing \(\vartheta\) and found that the average \(\vartheta\) temporal profiles in the selected events have the same “integral-like” shape as for \(\varphi\). The difference in \(\varphi\) and \(\vartheta\) average profiles is explained by the fact that most events have increasing profiles for the angle in the ecliptic plane as a result of solar rotation, while for the angle in the meridional plane, the numbers of events with increasing and decreasing profiles are equal.     

On the period of the periodic orbits of the restricted three body problem

We will show that the period T of a closed orbit of the planar circular restricted three body problem (viewed on rotating coordinates) depends on the region it encloses. Roughly speaking, we show that, \(2 T=k\pi +\int _\Omega g\) where k is an integer, \(\Omega \) is the region enclosed by the periodic orbit and \(g:{\mathbb {R}}^2\rightarrow {\mathbb {R}}\) is a function that only depends on the constant C known as the Jacobian constant; it does not depend on \(\Omega \). This theorem has a Keplerian flavor in the sense that it relates the period with the space “swept” by the orbit. As an application we prove that there is a neighborhood around \(L_4\) such that every periodic solution contained in this neighborhood must move clockwise. The same result holds true for \(L_5\).     

Relative equilibria of the restricted three-body problem in curved spaces

We use a formulation of the N-body problem in spaces of constant Gaussian curvature, \({\kappa }\in \mathbb {R}\), as widely used by A. Borisov, F. Diacu and their coworkers. We consider the restricted three-body problem in \(\mathbb {S}^2\) with arbitrary \({\kappa }>0\) (resp. \(\mathbb {H}^2\) with arbitrary \({\kappa }<0\)) in a formulation also valid for the case \({\kappa }=0\). For concreteness when \({\kappa }>0\) we restrict the study to the case of the three bodies at the upper hemisphere, to be denoted as \(\mathbb {S}^2_+\). The main goal is to obtain the totality of relative equilibria as depending on the parameters \({\kappa }\) and the mass ratio \(\mu \). Several general results concerning relative equilibria and its stability properties are proved analytically. The study is completed numerically using continuation from the \({\kappa }=0\) case and from other limit cases. In particular both bifurcations and spectral stability are also studied. The \(\mathbb {H}^2\) case is similar, in some sense, to the planar one, but in the \(\mathbb {S}^2_+\) case many differences have been found. Some surprising phenomena, like the coexistence of many triangular-like solutions for some values \(({\kappa },\mu )\) and many stability changes will be discussed.     

Spatial collinear restricted four-body problem with repulsive Manev potential

We outline some aspects of the dynamics of an infinitesimal mass under the Newtonian attraction of three point masses in a symmetric collinear relative equilibria configuration when a repulsive Manev potential (\(-1/r +e/r^{2}\)), \(e>0\), is applied to the central mass. We investigate the relative equilibria of the infinitesimal mass and their linear stability as a function of the mass parameter \(\beta \), the ratio of mass of the central body to the mass of one of two remaining bodies, and e. We also prove the nonexistence of binary collisions between the central body and the infinitesimal mass.     

A study of the geomagnetic indices asymmetry based on the interplanetary magnetic field polarities

Data of geomagnetic indices (aa, Kp, Ap, and Dst) recorded near 1 AU over the period 1967–2016, have been studied based on the asymmetry between the interplanetary magnetic field (IMF) directions above and below of the heliospheric current sheet (HCS). Our results led to the following conclusions: (i) Throughout the considered period, 31 random years (62%) showed apparent asymmetries between Toward (\(\mathbf{T}\)) and Away (\(\mathbf{A}\)) polarity days and 19 years (38%) exhibited nearly a symmetrical behavior. The days of \(\mathbf{A}\) polarity predominated over the \(\mathbf{T}\) polarity days by 4.3% during the positive magnetic polarity epoch (1991–1999). While the days of \(\mathbf{T}\) polarity exceeded the days of \(\mathbf{A}\) polarity by 5.8% during the negative magnetic polarity epoch (2001–2012). (ii) Considerable yearly North–South (N–S) asymmetries of geomagnetic indices observed throughout the considered period. (iii) The largest toward dominant peaks for \(aa\) and \(Ap\) indices occurred in 1995 near to minimum of solar activity. Moreover, the most substantial away dominant peaks for \(aa\) and \(Ap\) indices occurred in 2003 (during the descending phase of the solar cycle 23) and in 1991 (near the maximum of solar activity cycle) respectively. (iv) The N–S asymmetry of \(Kp\) index indicated a most significant away dominant peak occurred in 2003. (v) Four of the away dominant peaks of Dst index occurred at the maxima of solar activity in the years 1980, 1990, 2000, and 2013. The largest toward dominant peak occurred in 1991 (at the reversal of IMF polarity). (vi) The geomagnetic indices (aa, Ap, and \(Kp\)) all have northern dominance during positive magnetic polarity epoch (1971–1979), while the asymmetries shifts to the southern solar hemisphere during negative magnetic polarity epoch (2001–2012).     

$$M_{\bullet } - \sigma $$ relation in spherical systems

To investigate the \(M_\bullet -\sigma \) relation, we consider realistic elliptical galaxy profiles that are taken to follow a single power-law density profile given by \(\rho (r) = \rho _{0}(r/ r_{0})^{-\gamma }\) or the Nuker intensity profile. We calculate the density using Abel’s formula in the latter case by employing the derived stellar potential; in both cases. We derive the distribution function f(E) of the stars in the presence of the supermassive black hole (SMBH) at the center and hence compute the line-of-sight (LoS) velocity dispersion as a function of radius. For the typical range of values for masses of SMBH, we obtain \(M_{\bullet } \propto \sigma ^{p}\) for different profiles. An analytical relation \(p = (2\gamma + 6)/(2 + \gamma )\) is found which is in reasonable agreement with observations (for \(\gamma = 0.75{-}1.4\), \(p = 3.6{-}5.3\)). Assuming that a proportionality relation holds between the black hole mass and bulge mass, \(M_{\bullet } =f M_\mathrm{b}\), and applying this to several galaxies, we find the individual best fit values of p as a function of f; also by minimizing \(\chi ^{2}\), we find the best fit global p and f. For Nuker profiles, we find that \(p = 3.81 \pm 0.004\) and \(f = (1.23 \pm 0.09)\times 10^{-3}\) which are consistent with the observed ranges.     

No hanging out in neighborhoods of infinity in the three-body problem

Consider the spatial Newtonian three-body problem at fixed negative energy and fixed angular momentum. The moment of inertia I provides a measure of the overall size of a three-body system. We will prove that there is a positive number \(I_0\) depending on the energy and angular momentum levels as well as the masses such that every solution at these levels passes through \(I\le I_0\) at some instant of time. Motivation for this result comes from trying to prove the impossibility of realizing a certain syzygy sequence in the zero angular momentum problem.     

First integrals for the Kepler problem with linear drag

In this work we consider the Kepler problem with linear drag, and prove the existence of a continuous vector-valued first integral, obtained taking the limit as \(t\rightarrow +\infty \) of the Runge–Lenz vector. The norm of this first integral can be interpreted as an asymptotic eccentricity \(e_{\infty }\) with \(0\le e_{\infty } \le 1\). The orbits satisfying \(e_{\infty } <1\) approach the singularity by an elliptic spiral and the corresponding solutions \(x(t)=r(t)e^{i\theta (t)}\) have a norm r(t) that goes to zero like a negative exponential and an argument \(\theta (t)\) that goes to infinity like a positive exponential. In particular, the difference between consecutive times of passage through the pericenter, say \(T_{n+1} -T_n\), goes to zero as \(\frac{1}{n}\).     

Measurement of coronal properties of Seyfert galaxies from NuSTAR’s hard X-ray spectrum

Precise measurement of the coronal properties of Active Galactic Nuclei (AGN) requires the availability of high signal-to-noise ratio data covering a wide range of X-ray energies. The Nuclear Spectroscopic Telescope Array (NuSTAR) which is highly sensitive to earlier missions in its operational energy range of 3–79 keV, allows us to arrive at precise estimates of the coronal parameters such as cut-off energy (\(E_\mathrm{cut}\)), coronal temperature (\(\textit{kT}_e\)) and geometry of the corona at least for sources that have \(E_\mathrm{cut}\) within the energy range of NuSTAR. In this paper, we present our preliminary results on the spectral analysis of two Seyfert galaxies namely 3C 120 and NGC 4151 using NuSTAR observations in the 3–79 keV band. We investigated the continuum and coronal parameters, the photon index \(\Gamma \), \(E_\mathrm{cut}\) and \(\textit{kT}_{e}\). By fitting the X-ray spectrum of 3C 120 and NGC 4151 with a simple phenomenological model, we found that both the sources showed a clear cut-off in their spectrum.     

GONG p-Mode Parameters Through Two Solar Cycles

We investigate the parameters of global solar p-mode oscillations, namely damping width \(\Gamma\), amplitude \(A\), mean squared velocity \(\langle v^{2}\rangle\), energy \(E\), and energy supply rate \(\mathrm{d}E/\mathrm{d}t\), derived from two solar cycles’ worth (1996?–?2018) of Global Oscillation Network Group (GONG) time series for harmonic degrees \(l=0\,\mbox{--}\,150\). We correct for the effect of fill factor, apparent solar radius, and spurious jumps in the mode amplitudes. We find that the amplitude of the activity-related changes of \(\Gamma\) and \(A\) depends on both frequency and harmonic degree of the modes, with the largest variations of \(\Gamma\) for modes with \(2400~\upmu\mbox{Hz}\le\nu\le3300~\upmu\mbox{Hz}\) and \(31\le l \le60\) with a minimum-to-maximum variation of \(26.6\pm0.3\%\) and of \(A\) for modes with \(2400~\upmu\mbox{Hz}\le\nu\le 3300~\upmu\mbox{Hz}\) and \(61\le l \le100\) with a minimum-to-maximum variation of \(27.4\pm0.4\%\). The level of correlation between the solar radio flux \(F_{10.7}\) and mode parameters also depends on mode frequency and harmonic degree. As a function of mode frequency, the mode amplitudes are found to follow an asymmetric Voigt profile with \(\nu_{\text{max}}=3073.59\pm0.18~\upmu\mbox{Hz}\). From the mode parameters, we calculate physical mode quantities and average them over specific mode frequency ranges. In this way, we find that the mean squared velocities \(\langle v^{2}\rangle\) and energies \(E\) of p modes are anticorrelated with the level of activity, varying by \(14.7\pm0.3\%\) and \(18.4\pm0.3\%\), respectively, and that the mode energy supply rates show no significant correlation with activity. With this study we expand previously published results on the temporal variation of solar p-mode parameters. Our results will be helpful to future studies of the excitation and damping of p modes, i.e., the interplay between convection, magnetic field, and resonant acoustic oscillations.     

Secular resonances between bodies on close orbits: a case study of the Himalia prograde group of jovian irregular satellites

The gravitational interaction between two objects on similar orbits can effect noticeable changes in the orbital evolution even if the ratio of their masses to that of the central body is vanishingly small. Christou (Icarus 174:215–229, 2005) observed an occasional resonant lock in the differential node \(\varDelta \varOmega \) between two members in the Himalia irregular satellite group of Jupiter in the N-body simulations (corresponding mass ratio \(\sim 10^{-9}\)). Using a semianalytical approach, we have reproduced this phenomenon. We also demonstrate the existence of two additional types of resonance, involving angle differences \(\varDelta \omega \) and \(\varDelta (\varOmega +\varpi )\) between two group members. These resonances cause secular oscillations in eccentricity and/or inclination on timescales \(\sim \)1 Myr. We locate these resonances in (a, e, i) space and analyse their topological structure. In subsequent N-body simulations, we confirm these three resonances and find a fourth one involving \(\varDelta \varpi \). In addition, we study the occurrence rates and the stability of the four resonances from a statistical perspective by integrating 1000 test particles for 100 Myr. We find \(\sim \)10 to 30 librators for each of the resonances. Particularly, the nodal resonance found by Christou is the most stable: 2 particles are observed to stay in libration for the entire integration.