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.     

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.     

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.     

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.     

<Emphasis Type="Italic">GALEX</Emphasis> far-ultraviolet observations of metal-poor subdwarfs

Far-ultraviolet photometry derived from the GALEX satellite observatory has been compiled for a sample of metal-poor subdwarfs with \(\mathrm{[Fe/H]} < -1.0\). The FUV properties of these subdwarfs are compared with those of a set of Population I dwarfs that are known to have low levels of chromospheric activity. Comparisons are made via a number of photometric plots, including an absolute FUV magnitude versus \((V-K_{s})\) diagram, two-colour diagrams involving both \((m_{ \mathrm{FUV}}-B)\) and \((m_{\mathrm{FUV}}-V)\) versus \(B-V\), and a two-colour diagram composed of \((m_{\mathrm{FUV}}-V)\) versus \((V-K_{s})\). The warmest subdwarfs with \((V-K_{s}) \sim1.2\mbox{--}1.4\) show FUV excesses ranging from \(\sim2\mbox{--}3~\mbox{mag}\) relative to the Population I dwarfs, with the amount of FUV enhancement decreasing among subdwarfs of decreasing effective temperature. The coolest dwarfs that are compared have \((V-K_{s}) \sim1.8\), and among these stars the subdwarfs with \(-2.0 \leq{\mathrm{[Fe/H]}} \leq-1.0\) approach the locus of low activity Population I dwarfs in the \((m_{\mathrm{FUV}}-V, V-K_{s})\) diagram. In the \((m_{\mathrm{FUV}}-B, B-V)\) diagram the subdwarfs in this metallicity range overlap the Population I dwarf sequence for \((B-V) > 0.6\). The behaviour of the subdwarfs is consistent with their FUV fluxes being determined by a combination of a photospheric FUV spectrum, the strength of which diminishes towards cooler effective temperatures, and a spectrum of emission lines arising from a chromosphere and/or transition region which are of comparable strength between the coolest dwarfs and subdwarfs.     

CCD <Emphasis Type="Italic">UBV</Emphasis> photometric study of five open clusters—Dolidze 36, NGC 6728, NGC 6800, NGC 7209, and Platais 1

In this study, we present CCD UBV photometry of poorly studied open star clusters, Dolidze 36, NGC 6728, NGC 6800, NGC 7209, and Platais 1, located in the first and second Galactic quadrants. Observations were obtained with T100, the 1-m telescope of the TÜB?TAK National Observatory. Using photometric data, we determined several astrophysical parameters such as reddening, distance, metallicity and ages and from them, initial mass functions, integrated magnitudes and colours. We took into account the proper motions of the observed stars to calculate the membership probabilities. The colour excesses and metallicities were determined independently using two-colour diagrams. After obtaining the colour excesses of the clusters Dolidze 36, NGC 6728, NGC 6800, NGC 7209, and Platais 1 as \(0.19\pm0.06\), \(0.15\pm0.05\), \(0.32\pm0.05\), \(0.12\pm 0.04\), and \(0.43\pm0.06\) mag, respectively, the metallicities are found to be \(0.00\pm0.09\), \(0.02\pm0.11\), \(0.03\pm0.07\), \(0.01\pm0.08\), and \(0.01\pm0.08\) dex, respectively. Furthermore, using these parameters, distance moduli and age of the clusters were also calculated from colour-magnitude diagrams simultaneously using PARSEC theoretical models. The distances to the clusters Dolidze 36, NGC 6728, NGC 6800, NGC 7209, and Platais 1 are \(1050\pm90\), \(1610\pm190\), \(1210\pm150\), \(1060\pm90\), and \(1710\pm250\) pc, respectively, while corresponding ages are \(400\pm100\), \(750\pm150\), \(400\pm100\), \(600\pm100\), and \(175\pm50\) Myr, respectively. Our results are compatible with those found in previous studies. The mass function of each cluster is derived. The slopes of the mass functions of the open clusters range from 1.31 to 1.58, which are in agreement with Salpeter’s initial mass function. We also found integrated absolute magnitudes varying from ?4.08 to ?3.40 for the clusters.     

Corotating Shock Waves and the Solar-wind Source of Energetic Ion bundances: Power Laws in $$/Q$Q$

We find that element abundances in energetic ions accelerated by shock waves formed at corotating interaction regions (CIRs) mirror the abundances of the solar wind modified by a decreasing power-law dependence on the mass-to-charge ratio \(A\)/\(Q\) of the ions. This behavior is similar in character to the well-known power-law dependence on \(A\)/\(Q\) of abundances in large gradual solar energetic particles (SEP). The CIR ions reflect the pattern of \(A\)/\(Q\), with \(Q\) values of the source plasma temperature or freezing-in temperature of 1.0?–?1.2 MK typical of the fast solar wind in this case. Thus the relative ion abundances in CIRs are of the form \((A\mbox{/}Q)^{a}\), where \(a\) is nearly always negative and evidently decreases with distance from the shocks, which usually begin beyond 1 AU. For one unusual historic CIR event where \(a \approx 0\), the reverse shock wave of the CIR seems to occur at 1 AU, and these abundances of the energetic ions become a direct proxy for the abundances of the fast solar wind.     

Correlation Between the Magnetic Field and Plasma Parameters at 1 AU

The physical parameters of the solar wind observed in-situ near 1 AU have been studied for several decades, and relationships between them, such as the positive correlation between the solar wind plasma temperature, \(T\), and velocity, \(V\), and the negative correlation between density, \(N\), and velocity, \(V\), are well known. However, the magnetic field intensity, \(B\), does not appear to be well correlated with any individual plasma parameter. In this article, we discuss previously under-reported correlations between \(B\) and the combined plasma parameters \(\sqrt{N V^{2}} \) as well as between \(B\) and \(\sqrt{NT}\). These two correlations are strong during periods of corotating interaction regions and high-speed streams, and moderate during intervals of slow solar wind. The results indicate that the magnetic pressure in the solar wind is well correlated both with the plasma dynamic pressure and the thermal pressure.     

Role of primordial black holes in the direct collapse scenario of supermassive black hole formation at high redshifts

In this paper, we explore the possibility of accreting primordial black holes as the source of heating for the collapsing gas in the context of the direct collapse black hole scenario for the formation of super-massive black holes (SMBHs) at high redshifts, \(z\sim \) 6–7. One of the essential requirements for the direct collapse model to work is to maintain the temperature of the in-falling gas at \(\approx \)10\(^4\) K. We show that even under the existing abundance limits, the primordial black holes of masses \(\gtrsim \)10\(^{-2}M_\odot \), can heat the collapsing gas to an extent that the \(\mathrm{H}_2\) formation is inhibited. The collapsing gas can maintain its temperature at \(10^4\) K till the gas reaches a critical density \(n_{{c}} \,{\approx }\, 10^3~\hbox {cm}^{-3}\), at which the roto-vibrational states of \(\mathrm{H}_2\) approaches local thermodynamic equilibrium and \(\mathrm{H}_2\) cooling becomes inefficient. In the absence of \(\mathrm{H}_2\) cooling, the temperature of the collapsing gas stays at \(\approx \)10\(^4\) K even as it collapses further. We discuss scenarios of subsequent angular momentum removal and the route to find collapse through either a supermassive star or a supermassive disk.     

Relation Between Coronal Hole Areas and Solar Wind Speeds Derived from Interplanetary Scintillation Measurements

We investigate the relation between coronal hole (CH) areas and solar wind speeds during 1995?–?2011 using the potential field (PF) model analysis of magnetograph observations and interplanetary scintillation (IPS) observations by the Institute for Space-Earth Environmental Research (formerly Solar-Terrestrial Environment Laboratory) of Nagoya University. We obtained a significant positive correlation between the CH areas (\(A\)) derived from the PF model calculations and solar wind speeds (\(V\)) derived from the IPS observations. The correlation coefficients between them are usually high, but they drop significantly in solar maxima. The slopes of the \(A\)?–?\(V\) relation are roughly constant except for the period around solar maximum, when flatter or steeper slopes are observed. The excursion of the correlation coefficients and slopes at solar maxima is ascribed partly to the effect of rapid structural changes in the coronal magnetic field and solar wind, and partly to the predominance of small CHs. It is also demonstrated that \(V\) is inversely related to the flux expansion factor (\(f\)) and that \(f\) is closely related to \(A^{-1/2}\); hence, \(V \propto A^{1/2}\). A better correlation coefficient is obtained from the \(A^{1/2}\)?–?\(V\) relation, and this fact is useful for improving space weather predictions. We compare the CH areas derived from the PF model calculations with He i 1083 nm observations and show that the PF model calculations provide reliable estimates of the CH area, particularly for large \(A\).     

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.     

The Cadmium Zinc Telluride Imager on AstroSat

The Cadmium Zinc Telluride Imager (CZTI) is a high energy, wide-field imaging instrument on AstroSat. CZTI’s namesake Cadmium Zinc Telluride detectors cover an energy range from 20 keV to \(>200\) keV, with 11% energy resolution at 60 keV. The coded aperture mask attains an angular resolution of 17\(^\prime \) over a 4.6\(^\circ \) \(\times \) 4.6\(^\circ \)  (FWHM) field-of-view. CZTI functions as an open detector above 100 keV, continuously sensitive to GRBs and other transients in about 30% of the sky. The pixellated detectors are sensitive to polarization above \(\sim \)100 keV, with exciting possibilities for polarization studies of transients and bright persistent sources. In this paper, we provide details of the complete CZTI instrument, detectors, coded aperture mask, mechanical and electronic configuration, as well as data and products.     

The Pleiades apex and its kinematical structure

A study of cluster characteristics and internal kinematical structure of the middle-aged Pleiades open star cluster is presented. The individual star apexes and various cluster kinematical parameters including the velocity ellipsoid parameters are determined using both Hipparcos and Gaia data. Modern astrometric parameters were taken from the Gaia Data Release 1 (DR1) in combination with the Radial Velocity Experiment Fifth Data Release (DR5). The necessary set of parameters including parallaxes, proper motions and radial velocities are used for \(n=17\) stars from Gaia DR1+RAVE DR5 and for \(n=19\) stars from the Hipparcos catalog using SIMBAD data base. Single stars are used to improve accuracy by eliminating orbital movements. RAVE DR5 measurements were taken only for the stars with the radial velocity errors not exceeding \(2~\mbox{km}/\mbox{s}\). For the Pleiades stars taken from Gaia, we found mean heliocentric distance as \(136.8 \pm 6.4\) pc, and the apex position is calculated as: \(A_{CP}=92^{\circ }.52\pm 1^{\circ }.72\), \(D_{CP}=-42^{\circ }.28\pm 2^{\circ }.56\) by the convergent point method and \(A_{0}=95^{\circ }.59\pm 2^{\circ }.30\) and \(D_{0}=-50^{\circ }.90\pm 2^{\circ }.04\) using AD-diagram method (\(n=17\) in both cases). The results are compared with those obtained historically before the Gaia mission era.     

Fourier Power Spectra of Solar Noise Storms

We analyzed three noise storms recorded on 200?–?400 MHz Trieste Callisto radio spectra on 2 July 2012, 8 July 2012, and 16 July 2012 by the Fourier method. We divided intervals of the noise storms into five-minute intervals, and in these intervals we computed the mean Fourier spectra as a function of the wave numbers in the frequency and height-scale spaces. We found that these Fourier spectra, where the spectrum from the quiet-activity interval was subtracted, are power-law spectra. The mean power-law index of these spectra in the range \(\ln(k_{z}) = [1.8, 2.9]\) (where \(k_{z}\) is the wave number in the height-scale space) is \(-1.7\pm0.14\), \(-1.6\pm0.14\), and \(-1.5 \pm0.12\) for the 2 July 2012, the 8 July 2012, and the 16 July 2012 noise storms, respectively. It appears that as the number of Type-I bursts in the studied interval increases, the power-law index becomes closer to \(-5/3\); this is known as the Kolmogorov spectral index. The power-law index of the noise storms is very similar to that of the narrowband dm-spikes found in our previous studies. Furthermore, we found a break in the power spectra at \(\ln(k_{z}) \approx2.9\), and the mean power-law index values above this break are \(-2.9\pm0.46\), \(-3.1\pm0.65\), and \(-3.4\pm0.98\), respectively.     

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.     

Relating Solar Energetic Particle Event Fluences to Peak Intensities

Recently we (Kahler and Ling, Solar Phys.292, 59, 2017: KL) have shown that time–intensity profiles [\(I(t)\)] of 14 large solar energetic particle (SEP) events can be fitted with a simple two-parameter fit, the modified Weibull function, which is characterized by shape and scaling parameters [\(\alpha\) and \(\beta\)]. We now look for a simple correlation between an event peak energy intensity [\(I_{\mathrm{p}}\)] and the time integral of \(I(t)\) over the event duration: the fluence [\(F\)]. We first ask how the ratio of \(F/I_{\mathrm{p}}\) varies for the fits of the 14 KL events and then examine that ratio for three separate published statistical studies of SEP events in which both \(F\) and \(I_{\mathrm{p}}\) were measured for comparisons of those parameters with various solar-flare and coronal mass ejection (CME) parameters. The three studies included SEP energies from a 4?–?13 MeV band to \(E > 100~\mbox{MeV}\). Within each group of SEP events, we find a very robust correlation (\(\mathrm{CC} > 0.90\)) in log–log plots of \(F\)versus\(I_{\mathrm{p}}\) over four decades of \(I_{\mathrm{p}}\). The ratio increases from western to eastern longitudes. From the value of \(I_{\mathrm{p}}\) for a given event, \(F\) can be estimated to within a standard deviation of a factor of \({\leq}\,2\). Log–log plots of two studies are consistent with slopes of unity, but the third study shows plot slopes of \({<}\,1\) and decreasing with increasing energy for their four energy ranges from \(E > 10~\mbox{MeV}\) to \({>}\,100~\mbox{MeV}\). This difference is not explained.     

Intermittency of the Solar Magnetic Field and Solar Magnetic Activity Cycle

Small-scale solar magnetic fields demonstrate features of fractal intermittent behavior, which requires quantification. For this purpose we investigate how the observational estimate of the solar magnetic flux density \(B\) depends on resolution \(D\) in order to obtain the scaling \(\ln B_{D} = - k \ln D +a\) in a reasonably wide range. The quantity \(k\) demonstrates cyclic variations typical of a solar activity cycle. In addition, \(k\) depends on the magnetic flux density, i.e. the ratio of the magnetic flux to the area over which the flux is calculated, at a given instant. The quantity \(a\) demonstrates some cyclic variation, but it is much weaker than in the case of \(k\). The scaling obtained generalizes previous scalings found for the particular cycle phases. The scaling is typical of fractal structures. In our opinion, the results obtained trace small-scale action in the solar convective zone and its coexistence with the conventional large-scale solar dynamo based on differential rotation and mirror-asymmetric convection.     

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.     

Theoretical scaling law of coronal magnetic field and electron power-law index in solar microwave burst sources

It is first proposed a theoretical scaling law respectively for the coronal magnetic field strength \(B\) and electron power-law index \(\delta\) versus frequency and coronal height in solar microwave burst sources. Based on the non-thermal gyro-synchrotron radiation model (Ramaty in Astrophys. J. 158:753, 1969), \(B\) and \(\delta\) are uniquely solved by the observable optically-thin spectral index and turnover (peak) frequency, the other parameters (plasma density, temperature, view angle, low and high energy cutoffs, etc.) are relatively insensitive to the calculations, thus taken as some typical values. Both of \(B\) and \(\delta\) increase with increasing of radio frequency but with decreasing of coronal height above photosphere, and well satisfy a square or cubic logarithmic fitting.