Markarian galaxies in UV and multiwavelength studies

The UV properties of 1152 Markarian galaxies have been investigated based on GALEX data. These objects have been investigated also in other available wavelengths using multi-wavelength data from X-ray to radio. Using our classification for activity types for 779 Markarian galaxies based on SDSS spectroscopy, we have investigated these objects on the GALEX, 2MASS and WISE color-magnitude and color-color diagrams by the location of objects of different activity types and have revealed a number of loci. UV contours overplotted on the optical images revealed additional structures, particularly spiral arms of a number of Markarian galaxies. UV (FUV and NUV) and optical absolute magnitudes and luminosities have been calculated showing graduate transition from AGN to Composites, HIIs and Absorption line galaxies from (average \(M\)) \(-17.56^{m}\) to \(-15.20^{m}\) in FUV, from \(-18.07^{m}\) to \(-15.71^{m}\) in NUV and from AGN to Composites, Absorption line galaxies and HII from \(-21.14^{m}\) to \(-19.42^{m}\) in optical wavelengths and from (average \(L\)) \(7\times10^{9}\) to \(4 \times 10^{8}\) in FUV, from \(1\times 10^{10}\) to \(5\times10^{8}\) in NUV and from AGN to Composites, Absorption line galaxies and HII from \(7\times10^{10}\) to \(1\times10^{10}\) in optical wavelengths.     

Averaged model to study long-term dynamics of a probe about Mercury

This paper provides a method for finding initial conditions of frozen orbits for a probe around Mercury. Frozen orbits are those whose orbital elements remain constant on average. Thus, at the same point in each orbit, the satellite always passes at the same altitude. This is very interesting for scientific missions that require close inspection of any celestial body. The orbital dynamics of an artificial satellite about Mercury is governed by the potential attraction of the main body. Besides the Keplerian attraction, we consider the inhomogeneities of the potential of the central body. We include secondary terms of Mercury gravity field from \(J_2\) up to \(J_6\), and the tesseral harmonics \(\overline{C}_{22}\) that is of the same magnitude than zonal \(J_2\). In the case of science missions about Mercury, it is also important to consider third-body perturbation (Sun). Circular restricted three body problem can not be applied to Mercury–Sun system due to its non-negligible orbital eccentricity. Besides the harmonics coefficients of Mercury’s gravitational potential, and the Sun gravitational perturbation, our average model also includes Solar acceleration pressure. This simplified model captures the majority of the dynamics of low and high orbits about Mercury. In order to capture the dominant characteristics of the dynamics, short-period terms of the system are removed applying a double-averaging technique. This algorithm is a two-fold process which firstly averages over the period of the satellite, and secondly averages with respect to the period of the third body. This simplified Hamiltonian model is introduced in the Lagrange Planetary equations. Thus, frozen orbits are characterized by a surface depending on three variables: the orbital semimajor axis, eccentricity and inclination. We find frozen orbits for an average altitude of 400 and 1000 km, which are the predicted values for the BepiColombo mission. Finally, the paper delves into the orbital stability of frozen orbits and the temporal evolution of the eccentricity of these orbits.     

The Poynting-Robertson effect in the Newtonian potential with a Yukawa correction

We consider a Yukawa-type gravitational potential combined with the Poynting-Robertson effect. Dust particles originating within the asteroid belt and moving on circular and elliptic trajectories are studied and expressions for the time rate of change of their orbital radii and semimajor axes, respectively, are obtained. These expressions are written in terms of basic particle parameters, namely their density and diameter. Then, they are applied to produce expressions for the time required by the dust particles to reach the orbit of Earth. For the Yukawa gravitational potential, dust particles of diameter \(10^{ - 3}\) m in circular orbits require times of the order of \(8.557 \times 10^{6}\) yr and for elliptic orbits of eccentricities \(e =0.1, 0.5\) require times of \(9.396 \times 10^{6}\) and \(2.129 \times 10^{6}\) yr respectively to reach Earth’s orbit. Finally, various cases of the Yukawa potential are studied and the corresponding particle times to reach Earth’s are derived per case along with numerical results for circular and various elliptical orbits.     

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.     

Dust color temperature distribution of two FIR cavities at IRIS and AKARI maps

By systematically searching the region of far infrared loops, we found a number of huge cavity-like dust structures at \(60\,\mu \hbox {m}\) and \(100\,\mu \hbox {m}\) IRIS maps. By checking these with AKARI maps (\(90\,\mu \hbox {m}\) and \(140\,\mu \hbox {m}\)), two new cavity-like structures (sizes \(\sim \) \( 2.7\,\hbox {pc} \times 0.8\,\hbox {pc}\) and \(\sim \) \( 1.8\,\hbox {pc} \times 1\,\hbox {pc}\)) located at R.A. (\(\hbox {J}2000)=14^{h}41^{m}23^{s}\) and Dec. \((\hbox {J}2000)=-64^{\circ }04^{\prime }17^{{\prime }{\prime }}\) and R.A. \((\hbox {J}2000)=05^{h}05^{m}35^{s}\) and Dec. \((\hbox {J}2000)=-\,69^{\circ }35^{\prime } 25^{{\prime }{\prime }}\) were selected for the study. The difference in the average dust color temperatures calculated using IRIS and AKARI maps of the cavity candidates were found to be \(3.2\pm 0.9\,\hbox {K}\) and \(4.1\pm 1.2\,\hbox {K}\), respectively. Interestingly, the longer wavelength AKARI map gives larger values of dust color temperature than that of the shorter wavelength IRIS maps. Possible explanation of the results will be presented.     

Motion of the fast exoplanets

It is shown that a number of superfast, with periods \(< 2\) d, exoplanets revolve around parent stars with periods, near-commensurate with \(P_{E}\) and/or \(2 P_{E} / \pi\), where the exoplanet resonance timescale \(P_{E}=9603(85)\) s agrees fairly well with the period \(P_{0}= 9600.606(12)\) s of the so-called “cosmic oscillation” (the probability that the two timescales would coincide by chance is near \(3 \times10^{-4}\); the \(P_{0}\) period was discovered first in the Sun, and later on—in other objects of Cosmos). True nature of the exoplanet \(P_{0}\) resonance is unknown.     

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}\).     

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.     

Modeling the Global Coronal Field with Simulated Synoptic Magnetograms from Earth and the Lagrange Points $L_{3}$, L4$L_{4}$, and L5$L_{5}$

The solar photospheric magnetic flux distribution is key to structuring the global solar corona and heliosphere. Regular full-disk photospheric magnetogram data are therefore essential to our ability to model and forecast heliospheric phenomena such as space weather. However, our spatio-temporal coverage of the photospheric field is currently limited by our single vantage point at/near Earth. In particular, the polar fields play a leading role in structuring the large-scale corona and heliosphere, but each pole is unobservable for \({>}\,6\) months per year. Here we model the possible effect of full-disk magnetogram data from the Lagrange points \(L_{4}\) and \(L_{5}\), each extending longitude coverage by \(60^{\circ}\). Adding data also from the more distant point \(L_{3}\) extends the longitudinal coverage much further. The additional vantage points also improve the visibility of the globally influential polar fields. Using a flux-transport model for the solar photospheric field, we model full-disk observations from Earth/\(L_{1}\), \(L_{3}\), \(L_{4}\), and \(L_{5}\) over a solar cycle, construct synoptic maps using a novel weighting scheme adapted for merging magnetogram data from multiple viewpoints, and compute potential-field models for the global coronal field. Each additional viewpoint brings the maps and models into closer agreement with the reference field from the flux-transport simulation, with particular improvement at polar latitudes, the main source of the fast solar wind.     

On the photo-gravitational restricted four-body problem with variable mass

This paper deals with the photo-gravitational restricted four-body problem (PR4BP) with variable mass. Following the procedure given by Gascheau (C. R. 16:393–394, 1843) and Routh (Proc. Lond. Math. Soc. 6:86–97, 1875), the conditions of linear stability of Lagrange triangle solution in the PR4BP are determined. The three radiating primaries having masses \(m_{1}\), \(m_{2}\) and \(m_{3}\) in an equilateral triangle with \(m_{2}=m_{3}\) will be stable as long as they satisfy the linear stability condition of the Lagrangian triangle solution. We have derived the equations of motion of the mentioned problem and observed that there exist eight libration points for a fixed value of parameters \(\gamma (\frac{m \ \text{at time} \ t}{m \ \text{at initial time}}, 0<\gamma\leq1 )\), \(\alpha\) (the proportionality constant in Jeans’ law (Astronomy and Cosmogony, Cambridge University Press, Cambridge, 1928), \(0\leq\alpha\leq2.2\)), the mass parameter \(\mu=0.005\) and radiation parameters \(q_{i}, (0< q_{i}\leq1, i=1, 2, 3)\). All the libration points are non-collinear if \(q_{2}\neq q_{3}\). It has been observed that the collinear and out-of-plane libration points also exist for \(q_{2}=q_{3}\). In all the cases, each libration point is found to be unstable. Further, zero velocity curves (ZVCs) and Newton–Raphson basins of attraction are also discussed.     

Higher harmonics electrostatic ion cyclotron parallel flow velocity shear instability with inhomogeneous DC electric field in the magnetosphere of Saturn

In present paper higher harmonic electrostatic ion-cyclotron (EIC) parallel flow velocity shear instability in presence of perpendicular inhomogeneous DC electric field with the ambient magnetic field has been studied, in different regions of the magnetosphere of Saturn. Dimensionless growth rate variation of EIC waves has been observed with respect to \(k_{ \bot } \rho _{i}\) for various plasma parameters. Effect of velocity shear scale length (\(A_{i}\)), temperature anisotropy (\(T_{ \bot } /T_{\|}\)), magnetic field (\(B\)), electric field (\(E\)), inhomogeneity (\(P/a\)), angle of propagation (\(\theta \)), ratio of electron to ion temperature (\(T_{e}/T_{i}\)) and density gradient (\(\varepsilon _{n}\rho _{i}\)) on the growth of EIC waves in the inner magnetosphere of Saturn has been studied and analyzed. The mathematical formulation for dispersion relation and growth rate has been done by using the method of characteristic solution and kinetic approach. This theoretical analysis has been done taking the data from the Cassini in the inner magnetosphere of Saturn in the extended region where ion cyclotron waves have been observed. The change in the growth of these waves due to the presence of Enceladus has been analyzed.     

Investigation of the Geoeffectiveness of Disk-Centre Full-Halo Coronal Mass Ejections

We studied the occurrence and characteristics of geomagnetic storms associated with disk-centre full-halo coronal mass ejections (DC-FH-CMEs). Such coronal mass ejections (CMEs) can be considered as the most plausible cause of geomagnetic storms. We selected front-side full-halo coronal mass ejections detected by the Large Angle and Spectrometric Coronagraph onboard the Solar and Heliospheric Observatory (SOHO/LASCO) from the beginning of 1996 till the end of 2015 with source locations between solar longitudes E10 and W10 and latitudes N20 and S20. The number of selected CMEs was 66 of which 33 (50%) were deduced to be the cause of 30 geomagnetic storms with \(\mathrm{Dst} \leq- 50~\mbox{nT}\). Of the 30 geomagnetic storms, 26 were associated with single disk-centre full-halo CMEs, while four storms were associated, in addition to at least one disk-centre full-halo CME, also with other halo or wide CMEs from the same active region. Thirteen of the 66 CMEs (20%) were associated with 13 storms with \(-100~\mbox{nT} < \mbox{Dst} \leq- 50~\mbox{nT}\), and 20 (30%) were associated with 17 storms with \(\mbox{Dst}\leq- 100~\mbox{nT}\). We investigated the distributions and average values of parameters describing the DC-FH-CMEs and their interplanetary counterparts encountering Earth. These parameters included the CME sky-plane speed and direction parameter, associated solar soft X-ray flux, interplanetary magnetic field strength, \(B_{t}\), southward component of the interplanetary magnetic field, \(B_{s}\), solar wind speed, \(V_{sw}\), and the \(y\)-component of the solar wind electric field, \(E_{y}\). We found only a weak correlation between the Dst of the geomagnetic storms associated with DC-FH-CMEs and the CME sky-plane speed and the CME direction parameter, while the correlation was strong between the Dst and all the solar wind parameters (\(B_{t}\), \(B_{s}\), \(V_{sw}\), \(E_{y}\)) measured at 1 AU. We investigated the dependences of the properties of DC-FH-CMEs and the associated geomagnetic storms on different phases of solar cycles and the differences between Solar Cycles 23 and 24. In the rise phase of Solar Cycle 23 (SC23), five out of eight DC-FH-CMEs were geoeffective (\(\mbox{Dst} \leq- 50~\mbox{nT}\)). In the corresponding phase of SC24, only four DC-FH-CMEs were observed, three of which were nongeoeffective (\(\mbox{Dst} > - 50~\mbox{nT}\)). The largest number of DC-FH-CMEs occurred at the maximum phases of the cycles (21 and 17, respectively). Most of the storms with \(\mbox{Dst}\leq- 100~\mbox{nT}\) occurred at or close to the maximum phases of the cycles. When comparing the storms during epochs of corresponding lengths in Solar Cycles 23 and 24, we found that during the first 85 months of Cycle 23 the geoeffectiveness rate of the disk-centre full-halo CMEs was 58% with an average minimum value of the Dst index of \(- 146~\mbox{nT}\). During the corresponding epoch of Cycle 24, only 35% of the disk-centre full-halo CMEs were geoeffective with an average value of Dst of \(- 97~\mbox{nT}\).     

Fast Spinning of Planets

Spin periods of Jupiter, Saturn, Uranus and Neptune are specified by the analysis of the resonant motion of large satellites: \(P = 0.445(2)\,\hbox {d}\), 0.448(1) d, 0.673(9) d and 0.561(7) d, respectively. They occur to be near-commensurate with \(P_0=9600.606(12)\,\hbox {s}\), the period of the “cosmic” oscillation, discovered first in the Sun, then in other variable objects of the Universe. The like analysis of spin rates of the total set of the largest and fastest rotators of the Solar system (with mean diameters \(\ge 500\,\hbox {km}\) and \(P < 2\,\hbox {d}\),—of planets, asteroids and satellites) resulted in the best commensurable, or “synchronizing”, timescale 9594(65) s, coinciding fairly well with \(P_0\) too (the probability that the two timescales could agree by chance, is less than \(10^{-5}\)). True origin of this odd common resonance of our planetary system is unknown.     

Analysis of the exoplanet containing system Kepler-13

We have applied the close binary system analysis program WinFitter, with its physically detailed fitting function, to an intensive study of the complex multiple system Kepler-13 using photometry data from all 13 short cadence quarters downloaded from the NASA Exoplanet Archive (NEA) (http://exoplanetarchive.ipac.caltech.edu). The data-point error of our normalized, phase-sequenced and binned (380 points per bin: 0.00025 phase interval) flux values, at 14 ppm, allows the model’s specification for the mean reference flux level of the system to a precision better than 1 ppm. Our photometrically derived values for the mass and radius of KOI13.01 are \(6.8\pm0.6~\mbox{M}_{\mathrm{J}}\) and \(1.44\pm0.04~\mbox{R}_{\mathrm{J}}\). The star has a radius of \(1.67\pm0.05~\mbox{R}_{\odot}\). Our modelling sets the mean of the orbital inclination \(i\) at \(94.35\pm0.14^{\circ}\), with the star’s mean precession angle \(\phi_{p}\)—\(49.1\pm5.0^{\circ}\) and obliquity \(\theta_{o}\)\(67.9 \pm 3.0^{\circ}\), though there are known ambiguities about the sense in which such angles are measured.Our findings did not confirm secular variation in the transit modelling parameters greater than their full correlated errors, as argued by previous authors, when each quarter’s data was best-fitted with a determinable parameter set without prejudice. However, if we accept that most of the parameters remain the same for each transit, then we could confirm a small but steady diminution in the cosine of the orbital inclination over the 17 quarter timespan. This is accompanied by a slight increase of the star’s precession angle (less negative), but with no significant change in the obliquity of its spin axis. There are suggestions of a history of strong dynamical interaction with a highly distorted planet rotating in a 3:2 resonance with its revolution, together with a tidal lag of \(\sim30~\mbox{deg}\). The mean precessional period is derived to be about 1000 y, but at the present time the motion of the star’s rotation axis appears to be supporting the gravitational torque, rather than providing the balance against it that would be expected over long periods of time.The planet has a small but detectable backwarming effect on the star, which helps to explain the difference in brightness just after transit and just before occultation eclipses. In assessing these findings it is recognized that sources of uncertainty remain, notably with possible inherent micropulsational effects, variations from other components of the multiple star, stellar activity, differential rotation and the neglect of higher order terms (than \(r_{1}^{5}\)) in the fitting function, where \(r_{1}\) is the ratio of the radius of the star to the mean orbital separation of planet and host star.     

The Na <Emphasis Type="SmallCaps">i</Emphasis> and Sr <Emphasis Type="SmallCaps">ii</Emphasis> Resonance Lines in Solar Prominences

We estimate the electron density, \(n_{\mathrm{e}}\), and its spatial variation in quiescent prominences from the observed emission ratio of the resonance lines Na?i?5890 Å (D₂) and Sr?ii?4078 Å. For a bright prominence (\(\tau_{\alpha}\approx25\)) we obtain a mean \(n_{\mathrm{e}}\approx2\times10^{10}~\mbox{cm}^{-3}\); for a faint one (\(\tau _{\alpha }\approx4\)) \(n_{\mathrm{e}}\approx4\times10^{10}~\mbox{cm}^{-3}\) on two consecutive days with moderate internal fluctuation and no systematic variation with height above the solar limb. The thermal and non-thermal contributions to the line broadening, \(T_{\mathrm{kin}}\) and \(V_{\mathrm{nth}}\), required to deduce \(n_{\mathrm{e}}\) from the emission ratio Na?i/Sr?ii cannot be unambiguously determined from observed widths of lines from atoms of different mass. The reduced widths, \(\Delta\lambda_{\mathrm{D}}/\lambda_{0}\), of Sr?ii?4078 Å show an excess over those from Na?D₂ and \(\mbox{H}\delta\,4101\) Å, assuming the same \(T_{\mathrm{kin}}\) and \(V_{\mathrm{nth}}\). We attribute this excess broadening to higher non-thermal broadening induced by interaction of ions with the prominence magnetic field. This is suggested by the finding of higher macro-shifts of Sr?ii?4078 Å as compared to those from Na?D₂.     

Local stellar kinematics from RAVE data—VIII. Effects of the Galactic disc perturbations on stellar orbits of red clump stars

We aim to probe the dynamic structure of the extended Solar neighborhood by calculating the radial metallicity gradients from orbit properties, which are obtained for axisymmetric and non-axisymmetric potential models, of red clump (RC) stars selected from the RAdial Velocity Experiment’s Fourth Data Release. Distances are obtained by assuming a single absolute magnitude value in near-infrared, i.e. \(M_{Ks}=-1.54\pm0.04\) mag, for each RC star. Stellar orbit parameters are calculated by using the potential functions: (i) for the MWPotential2014 potential, (ii) for the same potential with perturbation functions of the Galactic bar and transient spiral arms. The stellar age is calculated with a method based on Bayesian statistics. The radial metallicity gradients are evaluated based on the maximum vertical distance (\(z_{max}\)) from the Galactic plane and the planar eccentricity (\(e_{p}\)) of RC stars for both of the potential models. The largest radial metallicity gradient in the \(0< z_{max} \leq0.5\) kpc distance interval is \(-0.065\pm0.005~\mbox{dex}\,\mbox{kpc}^{-1}\) for a subsample with \(e_{p}\leq0.1\), while the lowest value is \(-0.014\pm0.006~\mbox{dex}\,\mbox{kpc}^{-1}\) for the subsample with \(e_{p}\leq0.5\). We find that at \(z_{max}>1\) kpc, the radial metallicity gradients have zero or positive values and they do not depend on \(e_{p}\) subsamples. There is a large radial metallicity gradient for thin disc, but no radial gradient found for thick disc. Moreover, the largest radial metallicity gradients are obtained where the outer Lindblad resonance region is effective. We claim that this apparent change in radial metallicity gradients in the thin disc is a result of orbital perturbation originating from the existing resonance regions.     

Addendum to: Strength of the Solar Coronal Magnetic Field – A Comparison of Independent Estimates Using Contemporaneous Radio and White-Light Observations

This addendum uses an alternate fit for the electron density distribution \(N(r)\) (see Figure 1) and estimates the coronal magnetic field using the new model. We find that the estimates of the magnetic field are in close agreement using both the models.

We have fit the \(N(r)\) distribution obtained from STEREO-A/COR1 and SOHO/LASCO-C2 using a fifth-order polynomial (see Figure 1). The expression can be written as

$$\begin{aligned} N_{\text{cor}}(r) &= 1.43 \times 10^{9} r^{-5} - 1.91 \times 10^{9} r^{-4} + 1.07 \times 10^{9} r^{-3} - 2.87 \times 10^{8} r^{-2} \\ &\quad {} + 3.76 \times 10^{7} r^{-1} - 1.91 \times 10^{6} , \end{aligned}$$

(1)

where \(N_{\text{cor}}(r)\) is in units of cm^?3 and \(r\) is in units of \(\mathrm{R}_{\odot}\). The background coronal electron density is enhanced by a factor of 5.5 at 2.63 \(\mathrm{R}_{\odot}\) during the coronal mass ejection (CME). The estimated coronal magnetic field strength (\(B\)) using radio data indicates that \(B(r) \approx(0.51\text{\,--\,}0.48) \pm 0.02\ \mathrm{G}\) in the range \(r \approx2.65\text{\, --\,}2.82\ \mathrm{R}_{\odot}\). The field strengths for STEREO-A/COR1 and SOHO/LASCO-C2 are ≈?0.32 G at \(r \approx 3.11\ \mathrm{R}_{\odot}\) and ≈?0.12 G at \(r \approx 4.40\ \mathrm{R}_{\odot}\), respectively.

     

Hill stability of the satellite in the elliptic restricted four-body problem

We consider an elliptic restricted four-body system including three primaries and a massless particle. The orbits of the primaries are elliptic, and the massless particle moves under the mutual gravitational attraction. From the dynamic equations, a quasi-integral is obtained, which is similar to the Jacobi integral in the circular restricted three-body problem (CRTBP). The energy constant \(C\) determines the topology of zero velocity surfaces, which bifurcate at the equilibrium point. We define the concept of Hill stability in this problem, and a criterion for stability is deduced. If the actual energy constant \(C_{\mathrm{ac}}\ ( {>} 0 ) \) is bigger than or equal to the critical energy constant \(C_{\mathrm{cr}}\), the particle will be Hill stable. The critical energy constant is determined by the mass and orbits of the primaries. The criterion provides a way to capture an asteroid into the Earth–Moon system.     

Prediction of the Length of Upcoming Solar Cycles

The forecast of solar cycle (SC) characteristics is crucial particularly for several space-based missions. In the present study, we propose a new model for predicting the length of the SC. The model uses the information of the width of an autocorrelation function that is derived from the daily sunspot data for each SC. We tested the model on Versions 1 and 2 of the daily international sunspot number data for SCs 10?–?24. We found that the autocorrelation width \(A_{\mathrm{w}} ^{n}\) of SC \(n\) during the second half of its ascending phase correlates well with the modified length that is defined as \(T_{\mathrm{cy}}^{n+2} - T_{\mathrm{a}}^{n}\). Here \(T_{\mathrm{cy}}^{n+2}\) and \(T_{ \mathrm{a}}^{n}\) are the length and ascent time of SCs \(n+2\) and \(n\), respectively. The estimated correlation coefficient between the model parameters is 0.93 (0.91) for Version 1 (Version 2) sunspot series. The standard errors in the observed and predicted lengths of the SCs for Version 1 and Version 2 data are 0.38 and 0.44 years, respectively. The advantage of the proposed model is that the predictions of the length of the upcoming two SCs (i.e., \(n+1\), \(n+2\)) are readily available at the time of the peak of SC \(n\). The present model gives a forecast of 11.01, 10.52, and 11.91 years (11.01, 12.20, and 11.68 years) for the length of SCs 24, 25, and 26, respectively, for Version 1 (Version 2).