A model of the late universe with viscous Zel’ldovich fluid and decaying vacuum

Many have speculated about the presence of a stiff fluid in very early stage of the universe. Such a stiff fluid was first introduced by Zel’dovich. Recently the late acceleration of the universe was studied by taking bulk viscous stiff fluid as the dominant cosmic component, but the age predicted by such a model is less than the observed value. We consider a flat universe with viscous stiff fluid and decaying vacuum energy as the cosmic components and found that the model predicts a reasonable background evolution of the universe with de Sitter epoch as end phase of expansion. More over, the model also predicts a reasonable value for the age of the present universe. We also performed a dynamical system analysis of the model and found that the end de Sitter phase predicted by the model is stable.     

Dynamical system analysis of a five-dimensional cosmological model

A five-dimensional cosmological model including a single perfect fluid is studied in the framework of dynamical system analysis. All the critical points of the system are listed with their stability properties and some representative phase diagrams are explicitly shown. It is found that the stabilization of extra dimension is possible and the observed flatness of the three-dimensional space is provided for certain ranges of the equation of state parameter of the fluid. The model suggested here can be considered as a simplified model for examining the possible effects of the extra dimensions in the early universe.     

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.     

Element sets for high-order Poincaré mapping of perturbed Keplerian motion

The propagation and Poincaré mapping of perturbed Keplerian motion is a key topic in Celestial Mechanics and Astrodynamics, e.g., to study the stability of orbits or design bounded relative trajectories. The high-order transfer map (HOTM) method enables efficient mapping of perturbed Keplerian orbits using the high-order Taylor expansion of a Poincaré or stroboscopic map. The HOTM is only accurate close to the expansion point and therefore the number of revolutions for which the map is accurate tends to be limited. The proper selection of coordinates is of key importance for improving the performance of the HOTM method. In this paper, we investigate the use of different element sets for expressing the high-order map in order to find the coordinates that perform best in terms of accuracy. A new set of elements is introduced that enables extremely accurate mapping of the state, even for high eccentricities and higher-order zonal perturbations. Finally, the high-order map is shown to be very useful for the determination and study of fixed points and center manifolds of Poincaré maps.     

Fuel optimization for low-thrust Earth–Moon transfer via indirect optimal control

The problem of designing low-energy transfers between the Earth and the Moon has attracted recently a major interest from the scientific community. In this paper, an indirect optimal control approach is used to determine minimum-fuel low-thrust transfers between a low Earth orbit and a Lunar orbit in the Sun–Earth–Moon Bicircular Restricted Four-Body Problem. First, the optimal control problem is formulated and its necessary optimality conditions are derived from Pontryagin’s Maximum Principle. Then, two different solution methods are proposed to overcome the numerical difficulties arising from the huge sensitivity of the problem’s state and costate equations. The first one consists in the use of continuation techniques. The second one is based on a massive exploration of the set of unknown variables appearing in the optimality conditions. The dimension of the search space is reduced by considering adapted variables leading to a reduction of the computational time. The trajectories found are classified in several families according to their shape, transfer duration and fuel expenditure. Finally, an analysis based on the dynamical structure provided by the invariant manifolds of the two underlying Circular Restricted Three-Body Problems, Earth–Moon and Sun–Earth is presented leading to a physical interpretation of the different families of trajectories.     

Logarithmic spiral trajectories generated by Solar sails

Analytic solutions to continuous thrust-propelled trajectories are available in a few cases only. An interesting case is offered by the logarithmic spiral, that is, a trajectory characterized by a constant flight path angle and a fixed thrust vector direction in an orbital reference frame. The logarithmic spiral is important from a practical point of view, because it may be passively maintained by a Solar sail-based spacecraft. The aim of this paper is to provide a systematic study concerning the possibility of inserting a Solar sail-based spacecraft into a heliocentric logarithmic spiral trajectory without using any impulsive maneuver. The required conditions to be met by the sail in terms of attitude angle, propulsive performance, parking orbit characteristics, and initial position are thoroughly investigated. The closed-form variations of the osculating orbital parameters are analyzed, and the obtained analytical results are used for investigating the phasing maneuver of a Solar sail along an elliptic heliocentric orbit. In this mission scenario, the phasing orbit is composed of two symmetric logarithmic spiral trajectories connected with a coasting arc.     

Variational and symplectic integrators for satellite relative orbit propagation including drag

Orbit propagation algorithms for satellite relative motion relying on Runge–Kutta integrators are non-symplectic—a situation that leads to incorrect global behavior and degraded accuracy. Thus, attempts have been made to apply symplectic methods to integrate satellite relative motion. However, so far all these symplectic propagation schemes have not taken into account the effect of atmospheric drag. In this paper, drag-generalized symplectic and variational algorithms for satellite relative orbit propagation are developed in different reference frames, and numerical simulations with and without the effect of atmospheric drag are presented. It is also shown that high-order versions of the newly-developed variational and symplectic propagators are more accurate and are significantly faster than Runge–Kutta-based integrators, even in the presence of atmospheric drag.     

Capture into first-order resonances and long-term stability of pairs of equal-mass planets

Massive planets form within the lifetime of protoplanetary disks, and therefore, they are subject to orbital migration due to planet–disk interactions. When the first planet reaches the inner edge of the disk, its migration stops and consequently the second planet ends up locked in resonance with the first one. We detail how the resonant trapping works comparing semi-analytical formulae and numerical simulations. We restrict to the case of two equal-mass coplanar planets trapped in first-order resonances, but the method can be easily generalized. We first describe the family of resonant stable equilibrium points (zero-amplitude libration orbits) using series expansions up to different orders in eccentricity as well as a non-expanded Hamiltonian. Then we show that during convergent migration the planets evolve along these families of equilibrium points. Eccentricity damping from the disk leads to a final equilibrium configuration that we predict precisely analytically. The fact that observed multi-exoplanetary systems are rarely seen in resonances suggests that in most cases the resonant configurations achieved by migration become unstable after the removal of the protoplanetary disk. Here we probe the stability of the resonances as a function of planetary mass. For this purpose, we fictitiously increase the masses of resonant planets, adiabatically maintaining the low-amplitude libration regime until instability occurs. We discuss two hypotheses for the instability, that of a low-order secondary resonance of the libration frequency with a fast synodic frequency of the system, and that of minimal approach distance between planets. We show that secondary resonances do not seem to impact resonant systems at low amplitude of libration. Resonant systems are more stable than non-resonant ones for a given minimal distance at close encounters, but we show that the latter nevertheless play the decisive role in the destabilization of resonant pairs. We show evidence that as the planetary mass increases and the minimal distance between planets gets smaller in terms of mutual Hill radius, the region of stability around the resonance center shrinks, until the equilibrium point itself becomes unstable.     

Development of a Lunar Scintillometer as part of the national large optical telescope site survey

Ground layer turbulence is a very important site characterization parameter used to assess the quality of an astronomical site. The Lunar Scintillometer is a simple and effective site-testing device for measuring the ground layer turbulence. It consists of a linear array of photodiodes which are sensitive to the slight variations in the moon’s brightness due to scintillation by the lower layers of the Earth’s atmosphere. The covariance of intensity values between the non-redundant photodiode baselines can be used to measure the turbulence profile from the ground up to a height determined by the furthest pair of detectors. The six channel lunar scintillometer that has been developed at the Indian Institute of Astrophysics is based closely on an instrument built by the team led by Andrei Tokovinin of Cerro Tololo Inter-American Observatory (CTIO), Chile (Tokovinin et al., Mon. Not. R. Astron. Soc. 404(3), 1186–1196 2010). We have fabricated the instrument based on the existing electronic design, and have worked on the noise analysis, an EMI (Electromagnetic Induction) resistant PCB design and the software pipeline for analyzing the data from the same. The results from the instrument’s multi-year campaign at Mount Saraswati, Hanle is also presented.