Multi-wavelength correlation studies of diffuse emission from the Aquila Rift

We present a multi-wavelength correlation study of diffuse ultraviolet radiation using GALEX observations towards the Aquila Rift. Apart from airglow and zodiacal emissions, we find a diffuse background of \(1300\mbox{--} 3700~\mbox{ph}\,\mbox{cm}^{-2}\,\mbox{s}^{-1}\,\mbox{sr}^{-1}\,\mathring{\mathrm{A}}^{-1}\) in the far-ultraviolet (FUV, 1350–1750 Å) band and \(1300\mbox{--}2800~\mbox{ph}\,\mbox{cm}^{-2}\,\mbox{s}^{-1}\,\mbox{sr}^{-1}\,\mathring{\mathrm{A}}^{-1}\) in the near-ultraviolet (NUV, 1750–2850 Å) band. The observed diffuse UV emissions are saturated with total as well as neutral hydrogen column density in the region due to high optical depth in UV (\(\tau \), 0.91–23.38). Higher values of FUV/NUV ratio in the region, greater than the threshold value of 0.6, along with the positive correlation between the ratio and FUV intensity are due to excess emission in the FUV band which is absent in the NUV band. We estimated the excess emission to be in the range \(\sim 400\mbox{--} 2700~\mbox{ph}\,\mbox{cm}^{-2}\,\mbox{s}^{-1}\,\mbox{sr}^{-1}\,\mathring{\mathrm{A}}^{-1}\), plausibly due to H₂ fluorescence, ion line emissions and two-photon continuum emissions from the region in the FUV band, which also shows saturation in optically thick regions with N(H₂) as well as \(\mbox{H}\alpha \) emissions. Since N(H₂) and \(\mbox{H}\alpha \) emissions spread all over the region, the excess emission from the field is composite in nature and a detailed spectroscopic analysis is needed to disentangle the contribution from individual components.     

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

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

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.

     

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

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.     

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

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.     

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.     

Kelvin–Helmholtz Instability in a Cool Solar Jet in the Framework of Hall Magnetohydrodynamics: A Case Study

We investigate the conditions under which the magnetohydrodynamic (MHD) modes in a cylindrical magnetic flux tube moving along its axis become unstable against the Kelvin–Helmholtz (KH) instability. We use the dispersion relations of MHD modes obtained from the linearized Hall MHD equations for cool (zero beta) plasma by assuming real wave numbers and complex angular wave frequencies/complex wave phase velocities. The dispersion equations are solved numerically at fixed input parameters and varying values of the ratio \(l_{\mathrm{Hall}}/a\), where \(l_{\mathrm{Hall}} = c/\omega_{\mathrm{pi}}\) (\(c\) being the speed of light, and \(\omega_{\mathrm{pi}}\) the ion plasma frequency) and \(a\) is the flux tube radius. It is shown that the stability of the MHD modes depends upon four parameters: the density contrast between the flux tube and its environment, the ratio of external and internal magnetic fields, the ratio \(l_{\mathrm{Hall}}/a\), and the value of the Alfvén Mach number defined as the ratio of the tube axial velocity to Alfvén speed inside the flux tube. It is found that at high density contrasts, for small values of \(l_{\mathrm{Hall}}/a\), the kink (\(m = 1\)) mode can become unstable against KH instability at some critical Alfvén Mach number (or equivalently at critical flow speed), but a threshold \(l_{\mathrm{Hall}}/a\) can suppress the onset of the KH instability. At small density contrasts, however, the magnitude of \(l_{\mathrm{Hall}}/a\) does not affect noticeably the condition for instability occurrence – even though it can reduce the critical Alfvén Mach number. It is established that the sausage mode (\(m = 0\)) is not subject to the KH instability.     

Shannon Entropy-Based Prediction of Solar Cycle 25

A new model is proposed to forecast the peak sunspot activity of the upcoming solar cycle (SC) using Shannon entropy estimates related to the declining phase of the preceding SC. Daily and monthly smoothed international sunspot numbers are used in the present study. The Shannon entropy is the measure of inherent randomness in the SC and is found to vary with the phase of an SC as it progresses. In this model each SC with length \(T_{\mathrm{cy}}\) is divided into five equal parts of duration \(T_{\mathrm{cy}}/5\). Each part is considered as one phase, and they are sequentially termed P1, P2, P3, P4, and P5. The Shannon entropy estimates for each of these five phases are obtained for the \(n\)th SC starting from \(n=10\,\mbox{--}\,23\). We find that the Shannon entropy during the ending phase (P5) of the \(n\)th SC can be efficiently used to predict the peak smoothed sunspot number of the \((n+1)\)th SC, i.e. \(S_{\mathrm{max}}^{n+1}\). The prediction equation derived in this study has a good correlation coefficient of 0.94. A noticeable decrease in entropy from 4.66 to 3.89 is encountered during P5 of SCs 22 to 23. The entropy value for P5 of the present SC 24 is not available as it has not yet ceased. However, if we assume that the fall in entropy continues for SC 24 at the same rate as that for SC 23, then we predict the peak smoothed sunspot number of 63±11.3 for SC 25. It is suggested that the upcoming SC 25 will be significantly weaker and comparable to the solar activity observed during the Dalton minimum in the past.     

A New Solar Imaging System for Observing High-Speed Eruptions: <Emphasis Type="Italic">Solar Dynamics Doppler Imager</Emphasis> (SDDI)

A new solar imaging system was installed at Hida Observatory to observe the dynamics of flares and filament eruptions. The system (Solar Dynamics Doppler Imager; SDDI) takes full-disk solar images with a field of view of \(2520~\mbox{arcsec} \times 2520~\mbox{arcsec}\) at multiple wavelengths around the \(\mathrm{H}\alpha\) line at 6562 Å. Regular operation was started in May 2016, in which images at 73 wavelength positions spanning from \(\mathrm{H}\alpha -9~\mathring{\mathrm{A}}\) to \(\mathrm{H}\alpha +9~\mathring{\mathrm{A}}\) are obtained every 15 seconds. The large dynamic range of the line-of-sight velocity measurements (\({\pm}\,400~\mbox{km}\,\mbox{s}^{-1}\)) allows us to determine the real motions of erupting filaments in 3D space. It is expected that SDDI provides unprecedented datasets to study the relation between the kinematics of filament eruptions and coronal mass ejections (CME), and to contribute to the real-time prediction of the occurrence of CMEs that cause a significant impact on the space environment of the Earth.     

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.     

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.     

A Check on the Validity of Magnetic Field Reconstructions

We investigate a method to test whether a numerically computed model coronal magnetic field \({\boldsymbol {B}}\) departs from the divergence-free condition (also known as the solenoidality condition). The test requires a potential field \({\boldsymbol {B}}_{0}\) to be calculated, subject to Neumann boundary conditions, given by the normal components of the model field \({\boldsymbol {B}}\) at the boundaries. The free energy of the model field may be calculated using \(\frac{1}{2\mu _{0}}\int ({\boldsymbol {B}}-{\boldsymbol {B}}_{0})^{2}\mathrm{d}V\), where the integral is over the computational volume of the model field. A second estimate of the free energy is provided by calculating \(\frac{1}{2\mu _{0}}\int {\boldsymbol {B}}^{2}\,\mathrm{d}V-\frac{1}{2\mu _{0}}\int {\boldsymbol {B}}_{0}^{2}\,\mathrm{d}V\). If \({\boldsymbol {B}}\) is divergence free, the two estimates of the free energy should be the same. A difference between the two estimates indicates a departure from \(\nabla \cdot {\boldsymbol {B}}=0\) in the volume. The test is an implementation of a procedure proposed by Moraitis et al. (Solar Phys.289, 4453, 2014) and is a simpler version of the Helmholtz decomposition procedure presented by Valori et al. (Astron. Astrophys.553, A38, 2013). We demonstrate the test in application to previously published nonlinear force-free model fields, and also investigate the influence on the results of the test of a departure from flux balance over the boundaries of the model field. Our results underline the fact that, to make meaningful statements about magnetic free energy in the corona, it is necessary to have model magnetic fields that satisfy the divergence-free condition to a good approximation.     

Supernova neutrinos: Flavor conversion independent of their mass

In extremely dense neutrino environments like in supernova core, the neutrino-neutrino refraction may give rise to self-induced flavor conversion. These neutrino flavor oscillations are well understood from the idea of the exponentially growing modes of the interacting oscillators in the flavor space. Until recently, the growth rates of these modes were found to be of the order of the vacuum oscillation frequency \(\Delta m^2/2E\) [\(\mathcal {O}(1~\mathrm{km}^{-1})\)] and were considered slow growing. However, in the last couple of years it was found that if the system was allowed to have different zenith-angle distributions for the emitted \(\nu _e\) and \(\bar{\nu }_e\) beams then the fastest growing modes of the interacting oscillators grew at the order of \(\mu =\sqrt{2} G_\mathrm{F}n_{\nu }\), a typical \(\nu \)–\(\nu \) interaction energy [\(\mathcal {O}(10^5~\mathrm{km}^{-1})\)]. Thus the growth rates are very large in comparison to the so-called ‘slow oscillations’ and can result in neutrino flavor conversion on a much faster scale. In fact, the point that the growth rates are no longer dependent on the vacuum oscillation frequency \(\Delta m^2/2E\), makes these ‘fast flavor conversions’ independent of \(\Delta m^2\) (thus mass) and energy. This is a surprising result as neutrino flavor conversions are considered to be the ultimate proof of massive neutrinos. However, the importance of this effect in the realistic astrophysical scenarios still remains to be understood.     

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.     

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.     

The extreme initial kinetic energy allowed by a collapsing turbulent core

We present high-resolution hydrodynamical simulations aimed at following the gravitational collapse of a gas core, in which a turbulent spectrum of velocity is implemented only initially. We determine the maximal value of the ratio of kinetic energy to gravitational energy, denoted here by \((\frac{E_{\mathrm{kin}} }{E_{\mathrm{grav}}} )_{\max}\), so that the core (i) will collapse around one free-fall time of time evolution or (ii) will expand unboundedly, because it has a value of \(\frac{E_{\rm kin}}{E_{\mathrm{grav}}}\) larger than \(( \frac{E_{\mathrm{kin}}}{E_{\mathrm{grav}}} )_{\mathrm{max}}\). We consider core models with a uniform or centrally condensed density profile and with velocity spectra composed of a linear combination of one-half divergence-free turbulence type and the other half of a curl-free turbulence type. We show that the outcome of the core collapse are protostars forming either (i) a multiple system obtained from the fragmentation of filaments and (ii) a single primary system within a long filament. In addition, some properties of these protostars are also determined and compared with those obtained elsewhere.     

Resonance tongues in the linear Sitnikov equation

In this paper, we deal with a Hill’s equation, depending on two parameters \(e\in [0,1)\) and \(\varLambda >0\), that has applications to some problems in Celestial Mechanics of the Sitnikov type. Due to the nonlinearity of the eccentricity parameter e and the coexistence problem, the stability diagram in the \((e,\varLambda )\)-plane presents unusual resonance tongues emerging from points \((0,(n/2)^2),\ n=1,2,\ldots \) The tongues bounded by curves of eigenvalues corresponding to \(2\pi \)-periodic solutions collapse into a single curve of coexistence (for which there exist two independent \(2\pi \)-periodic eigenfunctions), whereas the remaining tongues have no pockets and are very thin. Unlike most of the literature related to resonance tongues and Sitnikov-type problems, the study of the tongues is made from a global point of view in the whole range of \(e\in [0,1)\). Indeed, an interesting behavior of the tongues is found: almost all of them concentrate in a small \(\varLambda \)-interval [1, 9 / 8] as \(e\rightarrow 1^-\). We apply the stability diagram of our equation to determine the regions for which the equilibrium of a Sitnikov \((N+1)\)-body problem is stable in the sense of Lyapunov and the regions having symmetric periodic solutions with a given number of zeros. We also study the Lyapunov stability of the equilibrium in the center of mass of a curved Sitnikov problem.     

On the Pulse Shape of Ground-Level Enhancements

We study the temporal intensity profile, or pulse shape, of cosmic ray ground-level enhancements (GLEs) by calculating the rise \(( \tau_{\mathrm{r}})\) and decay \((\tau_{\mathrm{d}})\) times for a small subset of all available events. Although these quantities show very large inter-event variability, a linear dependence of \(\tau_{ \mathrm{d}} \approx 3.5 \tau_{\mathrm{r}}\) is found. We interpret these observational findings in terms of an interplanetary transport model, thereby including the effects of scattering (in pitch-angle) as these particles propagate from (near) the Sun to Earth. It is shown that such a model can account for the observed trends in the pulse shape, illustrating that interplanetary transport must be taken into account when studying GLE events, especially their temporal profiles. Furthermore, depending on the model parameters, the pulse shape of GLEs may be determined entirely by interplanetary scattering, obscuring all information regarding the initial acceleration process, and hence making a classification between impulsive and gradual events, as is traditionally done, superfluous.