Advective accretion flow properties around rotating black holes – application to GRO J1655-40

We examine the properties of the viscous dissipative accretion flow around rotating black holes in the presence of mass loss. Considering the thin disc approximation, we self-consistently calculate the inflow-outflow solutions and observe that the mass outflow rates decrease with the increase in viscosity parameter (\(\alpha \)). Further, we carry out the model calculation of quasi-periodic oscillation frequency (\(\nu _{\mathrm{QPO}}\)) that is frequently observed in black hole sources and observe that \(\nu ^\mathrm{max}_{\mathrm{QPO}}\) increases with the increase of black hole spin (\(a_k\)). Then, we employ our model in order to explain the High Frequency Quasi-Periodic Oscillations (HFQPOs) observed in black hole source GRO J1655-40. While doing this, we attempt to constrain the range of \(a_k\) based on observed HFQPOs (\(\sim \)300 Hz and \(\sim \)450 Hz) for the black hole source GRO J1655-40.     

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.     

Existence of static wormhole solutions in $$ gravity

This work investigates some feasible regions for the existence of traversable wormhole geometries in \(f(R,G)\) gravity, where \(R\) and \(G\) represent the Ricci scalar and the Gauss-Bonnet invariant respectively. Three different matter contents anisotropic fluid, isotropic fluid and barotropic fluid have been considered for the analysis. Moreover, we split \(f(R,G)\) gravity model into Strobinsky like \(f(R)\) model and a power law \(f(G)\) model to explore wormhole geometries. We select red-shift and shape functions which are suitable for the existence of wormhole solutions for the chosen \(f(R,G)\) gravity model. It has been analyzed with the graphical evolution that the null energy and weak energy conditions for the effective energy-momentum tensor are usually violated for the ordinary matter content. However, some small feasible regions for the existence of wormhole solutions have been found where the energy conditions are not violated. The overall analysis confirms the existence of the wormhole geometries in \(f(R,G)\) gravity under some reasonable circumstances.     

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.     

An interacting new holographic dark energy in the framework of fractal cosmology

In this paper, we study an interacting holographic dark energy model in the framework of fractal cosmology. The features of fractal cosmology could pass ultraviolet divergencies and also make a better understanding of the universe in different dimensions. We discuss a fractal FRW universe filled with the dark energy and cold dark matter interacting with each other. It is observed that the Hubble parameter embraces the recent observational range while the deceleration parameter demonstrates an accelerating universe and a behavior similar to \(\Lambda \mbox{CDM}\). Plotting the equation of state shows that it lies in phantom region for interaction mode. We use \(\mathit{Om}\)-diagnostic tool and it shows a phantom behavior of dark energy which is a condition of avoiding the formation of black holes. Finally we execute the StateFinder diagnostic pair and all the trajectories for interacting and non-interacting state of the model meet the fixed point \(\Lambda \mbox{CDM}\) at the start of the evolution. A behavior similar to Chaplygin gas also can be observed in statefinder plane. We find that new holographic dark energy model (NHDE) in fractal cosmology expressed the consistent behavior with recent observational data and can be considered as a model to avoid the formation of black holes in comparison with the main model of NHDE in the simple FRW universe. It has also been observed that for the interaction term varying with matter density, the model generates asymptotic de-Sitter solution. However, if the interaction term varies with energy density, then the model shows Big-Rip singularity. Using our modified CAMB code, we observed that the interacting model suppresses the CMB spectrum at low multipoles \(l<50\) and enhances the acoustic peaks. Based on the observational data sets used in this paper and using Metropolis-Hastings method of MCMC numerical calculation, it seems that the best value with \(1\sigma \) and \(2\sigma \) confidence interval are \(\Omega _{m0}=0.278^{+0.008~+0.010} _{-0.007~-0.009}\), \(H_{0}=69.9^{+0.95~+1.57}_{-0.95~-1.57}\), \(r_{c}=0.08^{+0.02~+0.027}_{-0.002~-0.0027}\), \(\beta =0.496^{+0.005~+0.009} _{-0.005~-0.009}\), \(c= 0.691^{+0.024~+0.039}_{-0.025~-0.037}\) and \(b^{2}=0.035\) according to which we find that the proposed model in the presence of interaction is compatible with the recent observational data.     

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.     

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

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.     

$$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.     

An analysis of the convergence of Newton iterations for solving elliptic Kepler’s equation

In this note a study of the convergence properties of some starters \( E_0 = E_0(e,M)\) in the eccentricity–mean anomaly variables for solving the elliptic Kepler’s equation (KE) by Newton’s method is presented. By using a Wang Xinghua’s theorem (Xinghua in Math Comput 68(225):169–186, 1999) on best possible error bounds in the solution of nonlinear equations by Newton’s method, we obtain for each starter \( E_0(e,M)\) a set of values \( (e,M) \in [0, 1) \times [0, \pi ]\) that lead to the q-convergence in the sense that Newton’s sequence \( (E_n)_{n \ge 0}\) generated from \( E_0 = E_0(e,M)\) is well defined, converges to the exact solution \(E^* = E^*(e,M)\) of KE and further \( \vert E_n - E^* \vert \le q^{2^n -1}\; \vert E_0 - E^* \vert \) holds for all \( n \ge 0\). This study completes in some sense the results derived by Avendaño et al. (Celest Mech Dyn Astron 119:27–44, 2014) by using Smale’s \(\alpha \)-test with \(q=1/2\). Also since in KE the convergence rate of Newton’s method tends to zero as \( e \rightarrow 0\), we show that the error estimates given in the Wang Xinghua’s theorem for KE can also be used to determine sets of q-convergence with \( q = e^k \; \widetilde{q} \) for all \( e \in [0,1)\) and a fixed \( \widetilde{q} \le 1\). Some remarks on the use of this theorem to derive a priori estimates of the error \( \vert E_n - E^* \vert \) after n Kepler’s iterations are given. Finally, a posteriori bounds of this error that can be used to a dynamical estimation of the error are also obtained.     

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.     

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.     

A family of stacked central configurations in the planar five-body problem

We study planar central configurations of the five-body problem where three bodies, \(m_1, m_2\) and \(m_3\), are collinear and ordered from left to right, while the other two, \(m_4\) and \(m_5\), are placed symmetrically with respect to the line containing the three collinear bodies. We prove that when the collinear bodies form an Euler central configuration of the three-body problem with \(m_1=m_3\), there exists a new family, missed by Gidea and Llibre (Celest Mech Dyn Astron 106:89–107, 2010), of stacked five-body central configuration where the segments \(m_4m_5\) and \(m_1m_3\) do not intersect.     

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

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.     

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.     

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.     

Compressive and rarefactive double layers in non-uniform plasma with $$-nonextensive distributed electrons

Electrostatic solitary waves and double layers (DLs) formed by the coupled ion acoustic (IA) and drift waves have been investigated in non-uniform plasma using \(q\)-nonextensive distribution function for the electrons and assuming ions to be cold \(T_{i}< T_{e}\). It is found that both compressive and rarefactive nonlinear structures (solitary waves and DLs) are possible in such a system. The steeper gradients are supportive for compressive solitary (and double layers) and destructive for rarefactive ones. The \(q\)-nonextensivity parameter \(q\) and the magnitudes of gradient scale lengths of density and temperature have significant effects on the amplitude of the double layers (and double layers) as well as on the speed of these structures. This theoretical model is general which has been applied here to the \(F\)-region ionosphere for illustration.     

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.