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.     

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.     

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

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.     

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.     

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.     

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

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.     

Subsurface Zonal and Meridional Flow During Cycles 23 and 24

We study the solar-cycle variation of subsurface flows from the surface to a depth of 16 Mm. We have used ring-diagram analysis to analyze Dopplergrams obtained with the Michelson Doppler Imager (MDI) Dynamics Program, the Global Oscillation Network Group (GONG), and the Helioseismic and Magnetic Imager (HMI) instrument. We combined the zonal and meridional flows from the three data sources and scaled the flows derived from MDI and GONG to match those from HMI observations. In this way, we derived their temporal variation in a consistent manner for Solar Cycles 23 and 24. We have corrected the measured flows for systematic effects that vary with disk positions. Using time-depth slices of the corrected subsurface flows, we derived the amplitudes and times of the extrema of the fast and slow zonal and meridional flows during Cycles 23 and 24 at every depth and latitude. We find an average difference between maximum and minimum amplitudes of \(8.6 \pm0.4~\mbox{m}\,\mbox{s}^{-1}\) for the zonal flows and \(7.9 \pm0.3~\mbox{m}\,\mbox{s}^{-1}\) for the meridional flows associated with Cycle 24 averaged over a depth range from 2 to 12 Mm. The corresponding values derived from GONG data alone are \(10.5 \pm0.3~\mbox{m}\,\mbox{s}^{-1}\) for the zonal and \(10.8 \pm0.3~\mbox{m}\,\mbox{s}^{-1}\) for the meridional flow. For Cycle 24, the flow patterns are precursors of the magnetic activity. The timing difference between the occurrence of the flow pattern and the magnetic one increases almost linearly with increasing latitude. For example, the fast zonal and meridional flow appear \(2.1 \pm 0.6\) years and \(2.5\pm 0.6\) years, respectively, before the magnetic pattern at \(30^{\circ}\) latitude in the northern hemisphere, while in the southern hemisphere, the differences are \(3.2 \pm 1.2\) years and \(2.6 \pm 0.6\) years. The flow patterns of Cycle 25 are present and have reached \(30^{\circ}\) latitude. The amplitude differences of Cycle 25 are about 22% smaller than those of Cycle 24, but are comparable to those of Cycle 23. Moreover, polynomial fits of meridional flows suggest that equatorward meridional flows (counter-cells) might exist at about \(80^{\circ}\) latitude except during the declining phase of the solar cycle.     

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

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.     

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.     

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.     

Discovery of GeV gamma-ray emission from the LMC B0443-6657 with the <Emphasis Type="Italic">Fermi</Emphasis> Large Area Telescope

We report the discovery of gamma-ray detection from the Large Magellanic Cloud (LMC) B0443-6657 using the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. LMC B0443-6657 is a flat-spectrum radio source, possibly associated with a supernova remnant in the Large Magellanic Cloud (LMC N4). Employing the LAT data of 8 years, our results show a significant excess (\(>9.4\sigma \)) of gamma rays in the range of 0.2–100 GeV above the gamma-ray background. A power-law function is found to adequately describe the 0.2–\(100\mbox{ GeV}\)\(\gamma \)-ray spectrum, which yields a photon flux of \(3.27\pm 0.53\ \text{photon}\,\mbox{cm}^{2}\,\mbox{s}^{-1}\) with a photon index of \(2.35\pm 0.11\), corresponding to an isotropic gamma-ray luminosity of \(5.3\times 10^{40}~\mbox{erg}\,\mbox{s}^{-1}\). The hadronic model predicts a low X-ray and TeV flux while the leptonic model predicts an observable flux in these two energy bands. The follow-up observations of the LMC B0443-6657 in X-ray or TeV band would distinguish the radiation models of gamma rays from this region.     

Improved Determination of the Location of the Temperature Maximum in the Corona

The most used method to calculate the coronal electron temperature [\(T_{\mathrm{e}} (r)\)] from a coronal density distribution [\(n_{\mathrm{e}} (r)\)] is the scale-height method (SHM). We introduce a novel method that is a generalization of a method introduced by Alfvén (Ark. Mat. Astron. Fys. 27, 1, 1941) to calculate \(T_{\mathrm{e}}(r)\) for a corona in hydrostatic equilibrium: the “HST” method. All of the methods discussed here require given electron-density distributions [\(n_{\mathrm{e}} (r)\)] which can be derived from white-light (WL) eclipse observations. The new “DYN” method determines the unique solution of \(T_{\mathrm{e}}(r)\) for which \(T_{\mathrm{e}}(r \rightarrow \infty) \rightarrow 0\) when the solar corona expands radially as realized in hydrodynamical solar-wind models. The applications of the SHM method and DYN method give comparable distributions for \(T_{\mathrm{e}}(r)\). Both have a maximum [\(T_{\max}\)] whose value ranges between 1?–?3 MK. However, the peak of temperature is located at a different altitude in both cases. Close to the Sun where the expansion velocity is subsonic (\(r < 1.3\,\mathrm{R}_{\odot}\)) the DYN method gives the same results as the HST method. The effects of the other free parameters on the DYN temperature distribution are presented in the last part of this study. Our DYN method is a new tool to evaluate the range of altitudes where the heating rate is maximum in the solar corona when the electron-density distribution is obtained from WL coronal observations.     

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.     

Density variation effect on multi-ions with kinetic Alfven wave around cusp region—a kinetic approach

The kinetic Alfven waves in the presence of homogeneous magnetic field plasma with multi-ions effect are investigated. The dispersion relation and normalised damping rate are derived for low-\(\beta\) plasma using kinetic theory. The effect of density variation of \(\text{H}^{+}\), \(\text{He}^{+}\) and \(\text{O}^{+}\) ions is observed on frequency and damping rate of the wave. The variation of frequency (\(\omega\)) and normalised damping rate (\(\gamma / \varOmega_{H^{ +}} \)) of the wave are studied with respect to \(k_{ \bot} \rho_{j}\), where \(k_{ \bot} \) is the perpendicular wave number, \(\rho_{j}\) is the ion gyroradius and \(j \) denotes \(\text{H}^{+}\), \(\text{He}^{+}\) and \(\text{O}^{+}\) ions. The variation with \(k_{ \bot} \rho_{j}\) is considered over wide range. The parameters appropriate to cusp region are used for the explanation of results. It is found that with hydrogen and helium ions gyration, the frequency of wave is influenced by the density variation of \(\text{H}^{+}\) and \(\text{He}^{+}\) ions but remains insensitive to the change in density of \(\text{O}^{+}\) ions. For oxygen ion gyration, the frequency of wave varies over a short range only for \(\text{O}^{+}\) ion density variation. The wave shows damping at lower altitude due to variation in density of lighter \(\text{H}^{+}\) and \(\text{He}^{+}\) ions whereas at higher altitude only heavy \(\text{O}^{+}\) ions contribute in wave damping. The damping of wave may be due to landau damping or energy transfer from wave to particles. The present study signifies that the both lighter and heavier ions dominate differently to change the characteristics of kinetic Alfven wave and density variation is also an important parameter to understand wave phenomena in cusp region.     

Analysis of the exoplanet containing system Kepler-13

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