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.     

Star-formation rate in compact star-forming galaxies

We use the data for the \(\text{H}\beta\) emission-line, far-ultraviolet (FUV) and mid-infrared 22 μm continuum luminosities to estimate star formation rates \(\langle \mbox{SFR} \rangle \) averaged over the galaxy lifetime for a sample of about 14000 bursting compact star-forming galaxies (CSFGs) selected from the Data Release 12 (DR12) of the Sloan Digital Sky Survey (SDSS). The average coefficient linking \(\langle \mbox{SFR} \rangle \) and the star formation rate \(\mbox{SFR}_{0}\) derived from the \(\text{H}\beta\) luminosity at zero starburst age is found to be 0.04. We compare \(\langle \mbox{SFR} \rangle \mbox{s}\) with some commonly used SFRs which are derived adopting a continuous star formation during a period of \({\sim}\,100~\mbox{Myr}\), and find that the latter ones are 2–3 times higher. It is shown that the relations between SFRs derived using a geometric mean of two star-formation indicators in the UV and IR ranges and reduced to zero starburst age have considerably lower dispersion compared to those with single star-formation indicators. We suggest that our relations for \(\langle \mbox{SFR} \rangle \) determination are more appropriate for CSFGs because they take into account a proper temporal evolution of their luminosities. On the other hand, we show that commonly used SFR relations can be applied for approximate estimation within a factor of \({\sim}\,2\) of the \(\langle \mbox{SFR} \rangle \) averaged over the lifetime of the bursting compact galaxy.     

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.     

Secular resonances between bodies on close orbits: a case study of the Himalia prograde group of jovian irregular satellites

The gravitational interaction between two objects on similar orbits can effect noticeable changes in the orbital evolution even if the ratio of their masses to that of the central body is vanishingly small. Christou (Icarus 174:215–229, 2005) observed an occasional resonant lock in the differential node \(\varDelta \varOmega \) between two members in the Himalia irregular satellite group of Jupiter in the N-body simulations (corresponding mass ratio \(\sim 10^{-9}\)). Using a semianalytical approach, we have reproduced this phenomenon. We also demonstrate the existence of two additional types of resonance, involving angle differences \(\varDelta \omega \) and \(\varDelta (\varOmega +\varpi )\) between two group members. These resonances cause secular oscillations in eccentricity and/or inclination on timescales \(\sim \)1 Myr. We locate these resonances in (a, e, i) space and analyse their topological structure. In subsequent N-body simulations, we confirm these three resonances and find a fourth one involving \(\varDelta \varpi \). In addition, we study the occurrence rates and the stability of the four resonances from a statistical perspective by integrating 1000 test particles for 100 Myr. We find \(\sim \)10 to 30 librators for each of the resonances. Particularly, the nodal resonance found by Christou is the most stable: 2 particles are observed to stay in libration for the entire integration.     

Constraints on dissipation in the deep interiors of Ganymede and Europa from tidal phase-lags

Jupiter’s satellites are subject to strong tidal forces which result in variations of the gravitational potential and deformations of the satellites’ surfaces on the diurnal tidal cycle. Such variations are described by the Love numbers \(k_2\) and \(h_2\) for the tide-induced potential variation due to internal mass redistribution and the radial surface displacement, respectively. The phase-lags \( \phi _{k_2}\) and \( \phi _{h_2}\) of these complex numbers contain information about the rheological and dissipative states of the satellites. Starting from interior structure models and assuming a Maxwell rheology to compute the tidal deformation, we calculate the phase-lags in application to Ganymede and Europa. For both satellites we assume a decoupling of the outer ice-shell from the deep interior by a liquid subsurface water ocean. We show that, in this case, the phase-lag difference \(\varDelta \phi = \phi _{k_2}- \phi _{h_2}\) can provide information on the rheological and thermal state of the deep interiors if the viscosities of the deeper layers are small. In case of Ganymede, phase-lag differences can reach values of a few degrees for high-pressure ice viscosities \({<}10^{14}\) Pa s and would indicate a highly dissipative state of the deep interior. In this case \(\varDelta \phi \) is dominated by dissipation in the high-pressure ice layer rather than dissipation within the ice-I shell. These phase lags would be detectable from spacecraft in orbit around the satellite. For Europa \(\varDelta \phi \) could reach values exceeding \(20^\circ \) and phase-lag measurements could help distinguish between (1) a hot dissipative silicate mantle which would in thermal equilibrium correspond to a very thin outer ice-I shell and (2) a cold deep interior implying that dissipation would mainly occur in a thick (several tens of km) outer ice-I shell. These measurements are highly relevant for ESA’s Jupiter Icy Moons Explorer (JUICE) and NASA’s Europa Multiple Flyby Mission, both targeted for the Jupiter system.     

The Poynting-Robertson effect in the Newtonian potential with a Yukawa correction

We consider a Yukawa-type gravitational potential combined with the Poynting-Robertson effect. Dust particles originating within the asteroid belt and moving on circular and elliptic trajectories are studied and expressions for the time rate of change of their orbital radii and semimajor axes, respectively, are obtained. These expressions are written in terms of basic particle parameters, namely their density and diameter. Then, they are applied to produce expressions for the time required by the dust particles to reach the orbit of Earth. For the Yukawa gravitational potential, dust particles of diameter \(10^{ - 3}\) m in circular orbits require times of the order of \(8.557 \times 10^{6}\) yr and for elliptic orbits of eccentricities \(e =0.1, 0.5\) require times of \(9.396 \times 10^{6}\) and \(2.129 \times 10^{6}\) yr respectively to reach Earth’s orbit. Finally, various cases of the Yukawa potential are studied and the corresponding particle times to reach Earth’s are derived per case along with numerical results for circular and various elliptical orbits.     

The Most Intense Electron-Scale Current Sheets in the Solar Wind

Previous analysis of magnetohydrodynamic-scale currents in high-speed solar wind near 1 AU suggests that the most intense current-carrying structures occur at electron scales and are characterized by average current densities on the order of \(1~\mbox{pA}/\mbox{cm}^{2}\). Here, this prediction is verified by examining the effects of the measurement bandwidth and/or measurement resolution on the analysis of synthetic solar wind signals. Assuming Taylor’s hypothesis holds for the energetically dominant fluctuations at kinetic scales, the results show that when \(\nu_{c}\gg \nu_{b}\), where \(\nu_{c}\) is the measurement bandwidth and \(\nu_{b} \approx 1/3~\mbox{Hz}\) is the break frequency, the average scale of the most intense fluctuations in the current density proxy is approximately \(1/\nu_{c}\), and the average peak current density is a weakly increasing function that scales approximately like \(\nu_{c}^{0.1}\).     

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.

     

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.     

Geodesics around oscillatons made of exponential scalar field potential

Some of the spherically symmetric solutions to the Einstein–Klein–Gordon (EKG) equations can describe the astronomical soliton objects made of a real time-dependent scalar fields. The solutions are known as oscillatons which are non-singular satisfying the flatness conditions asymptotically with periodic (separated) time-dependency. In this paper, we investigate the geodesic motion around an oscillaton. The Spherically Symmetric Geometry allows the bound orbits in the plan \(\theta=\frac{\pi}{2}\) under a given initial conditions. The potential for the scalar field \(\varPhi=\varPhi(r,t)\), is an exponential function of the form \(V(\varPhi)=V_{0}\exp(\lambda\sqrt{k_{0}}\varPhi)\).     

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.     

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.     

The global modulation of Galactic and Jovian electrons in the heliosphere

A full three-dimensional, numerical model is used to study the modulation of Jovian and Galactic electrons from 1 MeV to 50 GeV, and from the Earth into the heliosheath. For this purpose the very local interstellar spectrum and the Jovian electron source spectrum are revisited. It is possible to compute the former with confidence at kinetic energies \(E < 50~\mbox{MeV}\) since Voyager 1 crossed the heliopause in 2012 at \(\sim 122~\mbox{AU}\), measuring Galactic electrons at these energies. Modeling results are compared with Voyager 1 observations in the outer heliosphere, including the heliosheath, as well as observations at or near the Earth from the ISSE3 mission, and in particular the solar minimum spectrum from the PAMELA space mission for 2009, also including data from Ulysses for 1991 and 1992, and observations above 1 MeV from SOHO/EPHIN. Making use of the observations at or near the Earth and the two newly derived input functions for the Jovian and Galactic electrons respectively, the energy range over which the Jovian electrons dominate the Galactic electrons is determined so that the intensity of Galactic electrons at Earth below 100 MeV is calculated. The differential intensity for the Galactic electrons at Earth for \(E = 1~\mbox{MeV}\) is \(\sim 4\) electrons \(\mbox{m}^{-2}\,\mbox{s}^{-1}\,\mbox{sr}^{-1}\,\mbox{MeV}^{-1}\), whereas for Jovian electrons it is \(\sim 350\) electrons \(\mbox{m}^{-2}\,\mbox{s}^{-1}\,\mbox{sr}^{-1}\,\mbox{MeV}^{-1}\). At \(E = 30~\mbox{MeV}\) the two intensities are the same; above this energy the Jovian electron intensity quickly subsides so that the Galactic intensity completely dominates. At 6 MeV, in the equatorial plane the Jovian electrons dominate but beyond \(\sim 15~\mbox{AU}\) the Galactic intensity begins to exceed the Jovian intensity significantly.     

Investigation of the Geoeffectiveness of Disk-Centre Full-Halo Coronal Mass Ejections

We studied the occurrence and characteristics of geomagnetic storms associated with disk-centre full-halo coronal mass ejections (DC-FH-CMEs). Such coronal mass ejections (CMEs) can be considered as the most plausible cause of geomagnetic storms. We selected front-side full-halo coronal mass ejections detected by the Large Angle and Spectrometric Coronagraph onboard the Solar and Heliospheric Observatory (SOHO/LASCO) from the beginning of 1996 till the end of 2015 with source locations between solar longitudes E10 and W10 and latitudes N20 and S20. The number of selected CMEs was 66 of which 33 (50%) were deduced to be the cause of 30 geomagnetic storms with \(\mathrm{Dst} \leq- 50~\mbox{nT}\). Of the 30 geomagnetic storms, 26 were associated with single disk-centre full-halo CMEs, while four storms were associated, in addition to at least one disk-centre full-halo CME, also with other halo or wide CMEs from the same active region. Thirteen of the 66 CMEs (20%) were associated with 13 storms with \(-100~\mbox{nT} < \mbox{Dst} \leq- 50~\mbox{nT}\), and 20 (30%) were associated with 17 storms with \(\mbox{Dst}\leq- 100~\mbox{nT}\). We investigated the distributions and average values of parameters describing the DC-FH-CMEs and their interplanetary counterparts encountering Earth. These parameters included the CME sky-plane speed and direction parameter, associated solar soft X-ray flux, interplanetary magnetic field strength, \(B_{t}\), southward component of the interplanetary magnetic field, \(B_{s}\), solar wind speed, \(V_{sw}\), and the \(y\)-component of the solar wind electric field, \(E_{y}\). We found only a weak correlation between the Dst of the geomagnetic storms associated with DC-FH-CMEs and the CME sky-plane speed and the CME direction parameter, while the correlation was strong between the Dst and all the solar wind parameters (\(B_{t}\), \(B_{s}\), \(V_{sw}\), \(E_{y}\)) measured at 1 AU. We investigated the dependences of the properties of DC-FH-CMEs and the associated geomagnetic storms on different phases of solar cycles and the differences between Solar Cycles 23 and 24. In the rise phase of Solar Cycle 23 (SC23), five out of eight DC-FH-CMEs were geoeffective (\(\mbox{Dst} \leq- 50~\mbox{nT}\)). In the corresponding phase of SC24, only four DC-FH-CMEs were observed, three of which were nongeoeffective (\(\mbox{Dst} > - 50~\mbox{nT}\)). The largest number of DC-FH-CMEs occurred at the maximum phases of the cycles (21 and 17, respectively). Most of the storms with \(\mbox{Dst}\leq- 100~\mbox{nT}\) occurred at or close to the maximum phases of the cycles. When comparing the storms during epochs of corresponding lengths in Solar Cycles 23 and 24, we found that during the first 85 months of Cycle 23 the geoeffectiveness rate of the disk-centre full-halo CMEs was 58% with an average minimum value of the Dst index of \(- 146~\mbox{nT}\). During the corresponding epoch of Cycle 24, only 35% of the disk-centre full-halo CMEs were geoeffective with an average value of Dst of \(- 97~\mbox{nT}\).     

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.     

Direct and indirect capture of near-Earth asteroids in the Earth–Moon system

Near-Earth asteroids have attracted attention for both scientific and commercial mission applications. Due to the fact that the Earth–Moon \(\hbox {L}_{1}\) and \(\hbox {L}_{2}\) points are candidates for gateway stations for lunar exploration, and an ideal location for space science, capturing asteroids and inserting them into periodic orbits around these points is of significant interest for the future. In this paper, we define a new type of lunar asteroid capture, termed direct capture. In this capture strategy, the candidate asteroid leaves its heliocentric orbit after an initial impulse, with its dynamics modeled using the Sun–Earth–Moon restricted four-body problem until its insertion, with a second impulse, onto the \(\hbox {L}_{2}\) stable manifold in the Earth–Moon circular restricted three-body problem. A Lambert arc in the Sun-asteroid two-body problem is used as an initial guess and a differential corrector used to generate the transfer trajectory from the asteroid’s initial obit to the stable manifold associated with Earth–Moon \(\hbox {L}_{2}\) point. Results show that the direct asteroid capture strategy needs a shorter flight time compared to an indirect asteroid capture, which couples capture in the Sun–Earth circular restricted three-body problem and subsequent transfer to the Earth–Moon circular restricted three-body problem. Finally, the direct and indirect asteroid capture strategies are also applied to consider capture of asteroids at the triangular libration points in the Earth–Moon system.     

New holographic dark energy model with constant bulk viscosity in modified $$ gravity theory

The aim of this paper is to study new holographic dark energy (HDE) model in modified \(f(R,T)\) gravity theory within the framework of a flat Friedmann-Robertson-Walker model with bulk viscous matter content. It is thought that the negative pressure caused by the bulk viscosity can play the role of dark energy component, and drive the accelerating expansion of the universe. This is the motive of this paper to observe such phenomena with bulk viscosity. In the specific model \(f(R,T)=R+\lambda T\), where \(R\) is the Ricci scalar, \(T\) the trace of the energy-momentum tensor and \(\lambda \) is a constant, we find the solution for non-viscous and viscous new HDE models. We analyze new HDE model with constant bulk viscosity, \(\zeta =\zeta _{0}= \text{const.}\) to explain the present accelerated expansion of the universe. We classify all possible scenarios (deceleration, acceleration and their transition) with possible positive and negative ranges of \(\lambda \) over the constraint on \(\zeta _{0}\) to analyze the evolution of the universe. We obtain the solutions of scale factor and deceleration parameter, and discuss the evolution of the universe. We observe the future finite-time singularities of type I and III at a finite time under certain constraints on \(\lambda \). We also investigate the statefinder and \(\mathit{Om}\) diagnostics of the viscous new HDE model to discriminate with other existing dark energy models. In late time the viscous new HDE model approaches to \(\varLambda \mathit{CDM}\) model. We also discuss the thermodynamics and entropy of the model and find that it satisfies the second law of thermodynamics.     

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.     

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

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