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

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

Establishing the Galactic Centre distance using VVV Bulge RR Lyrae variables

This study’s objective was to exploit infrared VVV (VISTA Variables in the Via Lactea) photometry for high latitude RRab stars to establish an accurate Galactic Centre distance. RRab candidates were discovered and reaffirmed (\(n=4194\)) by matching \(K_{s}\) photometry with templates via \(\chi ^{2}\) minimization, and contaminants were reduced by ensuring targets adhered to a strict period-amplitude (\(\Delta K_{s}\)) trend and passed the Elorietta et al. classifier. The distance to the Galactic Centre was determined from a high latitude Bulge subsample (\(|b|>4^{\circ}\), \(R_{\mathit{GC}}=8.30 \pm 0.36\) kpc, random uncertainty is relatively negligible), and importantly, the comparatively low color-excess and uncrowded location mitigated uncertainties tied to the extinction law, the magnitude-limited nature of the analysis, and photometric contamination. Circumventing those problems resulted in a key uncertainty being the \(M_{K_{s}}\) relation, which was derived using LMC RRab stars (\(M_{K_{s}}=-(2.66\pm 0.06) \log {P}-(1.03\pm 0.06)\), \((J-K_{s})_{0}=(0.31\pm 0.04) \log {P} + (0.35\pm 0.02)\), assuming \(\mu _{0,\mathit{LMC}}=18.43\)). The Galactic Centre distance was not corrected for the cone-effect. Lastly, a new distance indicator emerged as brighter overdensities in the period-magnitude-amplitude diagrams analyzed, which arise from blended RRab and red clump stars. Blending may thrust faint extragalactic variables into the range of detectability.     

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.     

Stellar parameters of the two binary systems: HIP 14075 and HIP 14230

We present the stellar parameters of the individual components of the two old close binary systems HIP 14075 and HIP 14230 using synthetic photometric analysis. These parameters are accurately calculated based on the best match between the synthetic photometric results within three different photometric systems with the observed photometry of the entire system. From the synthetic photometry, we derive the masses and radii of HIP 14075 as \({\mathcal {M}}^A=0.99\pm 0.19 \mathcal {M_\odot }\), \(R_{A}=0.877\pm 0.08 R_\odot \) for the primary and \({\mathcal {M}}^B=0.96\pm 0.15 \mathcal {M_\odot }\), \(R_{B}=0.821\pm 0.07 R_\odot \) for the secondary, and of HIP 14230 as \({\mathcal {M}}^A=1.18\pm 0.22 \mathcal {M_\odot }\), \(R_{A}=1.234\pm 0.05 R_\odot \) for the primary and \({\mathcal {M}}^B=0.84\pm 0.12 \mathcal {M_\odot }\) , \(R_{B}=0.820\pm 0.05 R_\odot \) for the secondary. Both systems depend on Gaia parallaxes. Based on the positions of the components of the two systems on a theoretical Hertzsprung–Russell diagram, we find that the age of HIP 14075 is \(11.5\pm 2.0\) Gyr and of HIP 14230 is \(3.5\pm 1.5\) Gyr. Our analysis reveals that both systems are old close binary systems (\(\approx > 4\) Gyr). Finally, the positions of the components of both the systems on the stellar evolutionary tracks and isochrones are discussed.     

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.

     

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

A new holographic dark energy model in Brans-Dicke theory with logarithmic scalar field

We study a holographic dark energy model in the framework of Brans-Dicke (BD) theory with taking into account the interaction between dark matter and holographic dark energy. We use the recent observational data sets, namely SN Ia compressed Joint Light-Analysis (cJLA) compilation, Baryon Acoustic Oscillations (BAO) from BOSS DR12 and the Cosmic Microwave Background (CMB) of Planck 2015. After calculating the evolution of the equation of state as well as the deceleration parameters, we find that with a logarithmic form for the BD scalar field the phantom crossing can be achieved in the late time of cosmic evolution. Unlike the conventional theory of holographic dark energy in standard cosmology (\(\omega_{D}=0\)), our model results in a late time accelerated expansion. It is also shown that the cosmic coincidence problem may be resolved in the proposed model. We execute the statefinder and Om diagnostic tools and demonstrate that interaction term does not play a significant role. Based on the observational data sets used in this paper it seems that the best value with \(1\sigma \) and \(2\sigma \) confidence interval are \(\varOmega_{m}=0.268^{+0.008~+0.010}_{-0.007~-0.009}\), \(\alpha =3.361^{+0.332~+0.483} _{-0.401~-0.522}\), \(\beta =5.560^{+0.541~+0.780}_{-0.510~-0.729}\), \(c=0.777^{+0.023~+0.029}_{-0.017~-0.023}\) and \(b^{2} =0.045\), according to which we find that the proposed model in the presence of interaction is compatible with the recent observational data.     

Superhump and outburst activity of the cataclysmic variable RZ LMi in the $$- and optical passbands

An analysis of the new \(U,B,V,R_{C},I_{C}\)-photometry of the cataclysmic variable RZ LMi obtained in 2016–17 showed the largest \((U-B)\) colour excess in quiescence as well as during the decline of brightness, associated with the outbursts activity. The smallest \((U-B)\) colour excess was found during the brightness increase from the quiescence. In contrast to the \((U-B)\) colour index, the \((B-V),(V-R_{C}),(R_{C}-I_{C})\) colour indices exhibits the largest colour excesses near the maximum of the outburst and the smallest during the quiescence. The \((B-V)\) colour index showed also a large excess 1–2 days before a minimum. The detailed study of superhumps during the maximum of activity reveals the largest \((U-B)\) colour excess at the time of the minimum brightness of superhumps. The \((B-R_{C})\) colour index exhibits a similar behaviour, but with a phase shift of \(+0.1\mbox{--}{+}0.2\) period of superhumps. The tracks in two-colour and colour-magnitude diagrams during superoutbursts are compared with the data for other cataclysmic variables during their outbursts as well as with published theoretical calculations.     

Stability analysis for non-minimally coupled dark energy models in the Palatini formalism

In this paper, we use the method of global analysis to study the stability of de-Sitter solutions in an universe dominated by a scalar field dark energy, which couples non-minimally with the Ricci scalar defined in the Palatini formalism. Effective potential and phase-space diagrams are introduced to describe qualitatively the de-Sitter solutions and their stabilities. We find that for the simple power-law function \(V(\varphi)=V_{0}\varphi^{n}\) there are no stable de-Sitter solutions. While for some more complicated potentials, i.e. \(V(\varphi)=V_{0}\varphi^{n}+\varLambda\) and \(V(\varphi)=V_{0} (e ^{-\lambda\varphi}+e^{\lambda\varphi})^{2}\), stable de-Sitter solutions can exist.     

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

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.     

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.     

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.     

Suggestion for searh of ethylene oxide ($$-C2H4O) in a osmi objet

Ethylene oxide (\(c\)-C₂H₄O) and its isomer acetaldehyde (CH₃CHO) are important organic molecules because of their potential role in the formation of amino acids. The \(c\)-C₂H₄O molecule is a \(b\)-type asymmetric top molecule and owing to half-spin of each of the four hydrogen atoms, it has two distinct ortho (nuclear spin one) and para (nuclear spin zero and two) species. It has been detected in the Sgr B2N. Using the rotational and centrifugal distortion constants along with the electric dipole moment, we have calculated energies of 100 rotational levels of each of the ortho and para species of \(c\)-C₂H₄O molecule and the Einstein \(A\)-coefficients for radiative transitions between the levels. The values of Einstein \(A\)-coefficients along with the scaled values for the collisional rate coefficients are used for solving a set of statistical equilibrium equations coupled with the equations of radiative transfer.Brightness-temperatures of five rotational transitions of each of the ortho and para species of \(c\)-C₂H₄O molecule are investigated. Out of these ten transitions, three transitions are found to show the anomalous absorption and rest seven are found to show the emission feature. We have also investigated seven transitions observed unblended in the Sgr B2(N). We have found that the transitions \(3_{3 0} - 3_{2 1}\) (23.134 GHz), \(2_{2 0} - 2_{1 1}\) (15.603 GHz), \(3_{3 1} - 3_{2 2}\) (39.680 GHz) and \(1_{1 1} - 0_{0 0}\) (39.582 GHz) may play important role for the identification of ethylene oxide in a cosmic object.     

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.     

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.     

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.     

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.