Anisotropic fluid distribution in higher dimensional general theory of relativity

We have obtained static and spherically symmetric self-gravitating solution of the field equations for anisotropic distribution of matter in higher- dimensional in the context of Einstein’s general theory of relativity. This work is an extension of the previous work of Hector Rago (Astrophys. Space Sci. 183:333, 1991) for four dimensional space-time. The solutions are matched to the analytical solutions for spherically symmetric self gravitating distribution of anisotropic matter obtained by Hector Rago (1991) for n=2.     

Tolman's cosmological models

Oscillating Tolman universes have been identified as solutions of the Friedmann equation with the bulk viscosity dissipation. Their dynamics and some thermodynamical properties are briefly discussed. Problems such as the Eternal Return philosophy and its modern incarnation into the doctrine of the Universe as a global oscillator, Tipler's no-return theorem, etc., are touched upon.     

Steady state obliquity of a rigid body in the spin–orbit resonant problem: application to Mercury

In this paper, we consider the elliptic collinear solutions of the classical n-body problem, where the n bodies always stay on a straight line, and each of them moves on its own elliptic orbit with the same eccentricity. Such a motion is called an elliptic Euler–Moulton collinear solution. Here we prove that the corresponding linearized Hamiltonian system at such an elliptic Euler–Moulton collinear solution of n-bodies splits into \((n-1)\) independent linear Hamiltonian systems, the first one is the linearized Hamiltonian system of the Kepler 2-body problem at Kepler elliptic orbit, and each of the other \((n-2)\) systems is the essential part of the linearized Hamiltonian system at an elliptic Euler collinear solution of a 3-body problem whose mass parameter is modified. Then the linear stability of such a solution in the n-body problem is reduced to those of the corresponding elliptic Euler collinear solutions of the 3-body problems, which for example then can be further understood using numerical results of Martínez et al. on 3-body Euler solutions in 2004–2006. As an example, we carry out the detailed derivation of the linear stability for an elliptic Euler–Moulton solution of the 4-body problem with two small masses in the middle.     

Radiating Fluid Spheres in the Effective Variables Approximation

We study the evolution of spherically symmetric radiating fluid distributions using the effective variables method, implemented ab initio in Schwarzschild coordinates. To illustrate the procedure and to establish some comparison with the original method, we integrate numerically the set of equations at the surface for two different models. The first model is derived from the Schwarzschild interior solution. The second model is inspired in the Tolman VI solution.     

General vacuum-universe solutions in Brans-Dicke cosmology

By a rescalation of the scalar field ? of the Jordan-Brans and Dicke cosmology, the general solutions of the Friedmannian ‘vacuum’ Universe are obtained. Only the flat space solution was previously known. Each solution is caracterized by the sign of the second time derivative of the rescaled field ψ≡?R ³ (R being the scale factor of the Robertson-Walker line-element): \(\ddot \psi\) = 0 (flat space), \(\ddot \psi\) < 0 (closed space), and \(\ddot \psi\) > 0 (open space), so that the solutions are mutually exclusive. Of these, the open space one is damped-oscillatory andR attains its absolute minimum, equal to zero, in only one of the two ‘extreme’ cycles. Otherwise,R min remains positive. If the ?-field is dominant near the singularity, these solutions may have physical significance. Also obtained, by the method mentioned above, is the general flat space solution for a ‘dust’ Universe and from it a closed space ‘dust’ solution. Both were found before by different authors, each one using a different method and, therefore, seemed up to now unrelated.     

An exact solution to determination of an open orbit

We present an exact solution of the equations for orbit determination of a two body system in a hyperbolic or parabolic motion. In solving this problem, we extend the method employed by Asada, Akasaka and Kasai (AAK) for a binary system in an elliptic orbit. The solutions applicable to each of elliptic, hyperbolic and parabolic orbits are obtained by the new approach, and they are all expressed in an explicit form, remarkably, only in terms of elementary functions. We show also that the solutions for an open orbit are recovered by making a suitable transformation of the AAK solution for an elliptic case.     

Formation of voids in the Universe within the Lemaître–Tolman model

We develop models of void formation starting from a small initial fluctuation at recombination and growing to a realistic present-day density profile in agreement with observations of voids. The model construction is an extension of previously developed algorithms for finding a Lemaître–Tolman metric that evolves between two profiles of either density or velocity specified at two times. Of the four profiles of concern (those of density and velocity at recombination and at the present day), two can be specified and the other two follow from the derived model.
We find that, in order to reproduce the present-day void density profiles, the initial velocity profile is more important than the initial density profile.
Extrapolation of current cosmic microwave background (CMB) observations to the scales relevant to protovoids is very uncertain. Even so, we find that it is very difficult to make both the initial density and velocity fluctuation amplitudes small enough and still obtain a realistic void by today.     

A new approach to impulsive rendezvous near circular orbit

A new approach is presented for the problem of planar optimal impulsive rendezvous of a spacecraft in an inertial frame near a circular orbit in a Newtonian gravitational field. The total characteristic velocity to be minimized is replaced by a related characteristic-value function and this related optimization problem can be solved in closed form. The solution of this problem is shown to approach the solution of the original problem in the limit as the boundary conditions approach those of a circular orbit. Using a form of primer-vector theory the problem is formulated in a way that leads to relatively easy calculation of the optimal velocity increments. A certain vector that can easily be calculated from the boundary conditions determines the number of impulses required for solution of the optimization problem and also is useful in the computation of these velocity increments. Necessary and sufficient conditions for boundary conditions to require exactly three nonsingular non-degenerate impulses for solution of the related optimal rendezvous problem, and a means of calculating these velocity increments are presented. A simple example of a three-impulse rendezvous problem is solved and the resulting trajectory is depicted. Optimal non-degenerate nonsingular two-impulse rendezvous for the related problem is found to consist of four categories of solutions depending on the four ways the primer vector locus intersects the unit circle. Necessary and sufficient conditions for each category of solutions are presented. The region of the boundary values that admit each category of solutions of the related problem are found, and in each case a closed-form solution of the optimal velocity increments is presented. Similar results are presented for the simpler optimal rendezvous that require only one-impulse. For brevity degenerate and singular solutions are not discussed in detail, but should be presented in a following study. Although this approach is thought to provide simpler computations than existing methods, its main contribution may be in establishing a new approach to the more general problem.     

A comparative study of Newtonian,Kustaanheimo/Stiefel,and Sperling/Burdet optimal trajectories

An optimal trajectory problem is formulated in each of three sets of equations, and the resulting solutions are numerically compared. The three formulations are the classical Newtonian (N), the Kustaanheimo/Stiefel (K/S), and the Sperling/Burdet (S/B). The last two solutions are first regularized by the classical Sundman technique and the K/S solution is transformed before the optimization problem is posed. A novel technique is developed for generating initial control vectors for each solution. Numerically generated derivatives (central differences) are used by a type of gradient, Newton-Raphson iterator to converge the two-point boundary value problems. The results indicate that, although the K/S and S/B formulations are more difficult to express mathematically than the Newtonian formulation, the transformed solutions are significantly more numerically stable than the Newtonian solution when the perturbing acceleration is less than a minimum value (T/W _o=0.05 for the particular example problem treated).     

Symetries du Probleme desN Corps,Regularisation des Solutions Centrales Singulieres

The newtonian problem ofn mass points bodies is invariant by several changes of spatio-temporal variables. These symmetries correspond to arbitrary choices of the referential and they are related via Noether's theorem or by its generalization to conservative quantities of the motion. Forn=2 the author has defined two families of symmetriesS ₁ andS ₂ changing the eccentricity of a solution. The family of symmetries,S ₁, is associated to the arbitrary choice of thezero level of the potential and may related unbounded and bounded solutions. The family of symmetries,S ₂, is related to a possibleaffinity of the configurations space. Via a symmetry of theS ₂ family a zero angular momentum solution is equivalent to a non-zero angular momentum solution. Via a product of two symmetries of each family, denoted byS ₁.S ₂, any solution of the two-body problem is equivalent to a circular solution. In this paper it is shown that some of these transformations may be generalized to symmetries changing the quantityC ² H in then-body problem, whereC is the angular momentum andH is the energy. The extension is easily made to central solutions of then-body problem because involving several synchroneous two-body problems. We consider for exposition then=3 case. The principal results may be resumed by the following propositions:

1. The two families of symmetriesS ₁ andS ₂ are described by a spatial transformation product of an instantaneous homothethy and an instantaneous rotation completed by a change of temporal variable.
2. TheS ₁ family of symmetries may relate unbounded and bounded central solutions of the same type, i.e. unaligned or aligned.
3. TheS ₂ family of symmetries may regularize multiple collisions among central solutions of the same type.

Therefore any central solution, via a symmetryS ₁ orS ₂ orS ₁.S ₂, is equivalent to a central circular solution of the same type. That is a form of regularization.     

Free time minimizers for the three-body problem

Free time minimizers of the action (called “semi-static” solutions by Mañe in International congress on dynamical systems in Montevideo (a tribute to Ricardo Mañé), vol 362, pp 120–131, 1996) play a central role in the theory of weak KAM solutions to the Hamilton–Jacobi equation (Fathi in Weak KAM Theorem in Lagrangian Dynamics Preliminary Version Number 10, 2017). We prove that any solution to Newton’s three-body problem which is asymptotic to Lagrange’s parabolic homothetic solution is eventually a free time minimizer. Conversely, we prove that every free time minimizer tends to Lagrange’s solution, provided the mass ratios lie in a certain large open set of mass ratios. We were inspired by the work of Da Luz and Maderna (Math Proc Camb Philos Soc 156:209–227, 1980) which showed that every free time minimizer for the N-body problem is parabolic and therefore must be asymptotic to the set of central configurations. We exclude being asymptotic to Euler’s central configurations by a second variation argument. Central configurations correspond to rest points for the McGehee blown-up dynamics. The large open set of mass ratios are those for which the linearized dynamics at each Euler rest point has a complex eigenvalue.     

A steady-state solution for warped accretion discs

We consider a thin accretion disc warped due to the Bardeen–Petterson effect, presenting both analytical and numerical solutions for the situation in which the two viscosity coefficients vary with radius as a power law, with the two power-law indices not necessarily equal. The analytical solutions are compared with numerical ones, showing that our new analytical solution is more accurate than the previous one, which overestimated the inclination change in the outer disc. Our new analytical solution is appropriate for moderately warped discs, while for extremely misaligned discs only a numerical solution is appropriate.     

Determining parameters of Moon’s orbital and rotational motion from LLR observations using GRAIL and IERS-recommended models

The aim of this work is to combine the model of orbital and rotational motion of the Moon developed for DE430 with up-to-date astronomical, geodynamical, and geo- and selenophysical models. The parameters of the orbit and physical libration are determined in this work from lunar laser ranging (LLR) observations made at different observatories in 1970–2013. Parameters of other models are taken from solutions that were obtained independently from LLR. A new implementation of the DE430 lunar model, including the liquid core equations, was done within the EPM ephemeris. The postfit residuals of LLR observations make evident that the terrestrial models and solutions recommended by the IERS Conventions are compatible with the lunar theory. That includes: EGM2008 gravitational potential with conventional corrections and variations from solid and ocean tides; displacement of stations due to solid and ocean loading tides; and precession-nutation model. Usage of these models in the solution for LLR observations has allowed us to reduce the number of parameters to be fit. The fixed model of tidal variations of the geopotential has resulted in a lesser value of Moon’s extra eccentricity rate, as compared to the original DE430 model with two fit parameters. A mixed model of lunar gravitational potential was used, with some coefficients determined from LLR observations, and other taken from the GL660b solution obtained from the GRAIL spacecraft mission. Solutions obtain accurate positions for the ranging stations and the five retroreflectors. Station motion is derived for sites with long data spans. Dissipation is detected at the lunar fluid core-solid mantle boundary demonstrating that a fluid core is present. Tidal dissipation is strong at both Earth and Moon. Consequently, the lunar semimajor axis is expanding by 38.20 mm/yr, the tidal acceleration in mean longitude is \(-25.90 {{}^{\prime \prime }}/\mathrm{cy}^2\), and the eccentricity is increasing by \(1.48\times 10^{-11}\) each year.     

Global solutions of adiabatic accretion flows with isothermal shocks in Kerr black hole geometry

This paper presents global solutions of adiabatic accretion flows with isothermal shocks in Kerr black hole geometry. It is known that in the previously studied cases, where the flow including the shock is either entirely adiabatic or entirely isothermal, there can be no more than one stable shock solution, and the solution can only be of α –x type. However, the solution topology in the present case shows remarkable new characteristics: for the same flow parameters there can be two stable shock solutions satisfying physical boundary conditions, and the solution can be of three types, namely α– x , x –α and α–α type. In addition, shocks in the present case occur for a parameter region different from that for Rankine–Hugoniot shocks. These results greatly increase the possibilities of shock formation in astrophysical flows. It is also significant that the effects of frame-dragging of a rapid Kerr black hole on the shock formation are discovered. Finally, a brief comparison is made between shocked inviscid flows and two types of shock-free viscous flows, namely those of Shakura & Sunyaev and Narayan & Yi, and some comments are made about the fact that numerous authors who have studied transonic global solutions of accretion flows have found no shocks.     

Relativistic parametric embedding class I solutions of cold stars in Karmarkar space-time continuum

The aim of this paper is to explore a new parametric class of relativistic solutions to the Einstein field equations describing a spherically symmetric, static distribution of anisotropic fluid spheres to study the behavior of some of the cold stars in the setting of Karmarkar space-time continuum. We develop models of stellar objects for a range of parameter values of n and analyze their behavior through graphical representation. For each of these models, we have found that the metric potentials are well behaved inside the stellar interior and the physical parameters such as density, radial and tangential pressures, red-shift, radial speed, radial pressure density ratio and energy conditions display a continuous decrease from the center to surface of the stars whereas the mass, anisotropy, adiabatic indexes and compactification factor show a monotonous increase which imply that the proposed solution satisfy all the physical aspects of a realistic stellar objects. The stability of the solutions are verified by examining various stability aspects, viz., Zeldovich criteria, causality condition, Bondi condition, equilibrium condition (TOV-equation) and stable static criteria in connection to their cogency.     

RR trianguli australis: A new analysis of a semi-detached system of unusually small dimensions

RR TrA is an Algol-type eclipsing binary with a period of 0 _. ^d 7131. The components are believed to be of spectral types B9 V and F0 IV. The analysis presented is based on previously unreduced observations made by Somerville in 1961–63. Orbital elements are obtained using the Wood model, and these indicate that the system is semi-detached, and that the eclipses are marginally complete. The components appear to be of unusually small size for their spectral classes. The solution presented is compared with earlier ones for this system obtained by Kwee and by Giuricin and Mardirossian, and similarities and differences between these solutions are noted.     

Orbites periodiques de satellites

In a previous paper (Stellmacher, 1981, hereafter mentioned as Paper I), we have given an algorithm for the construction of periodic orbits in a rotating frame, for satellites around an oblate planet. In the present paper, we apply this theory to the Mimas-Tethys case; we obtain the following results:

1. Without resonance, it is possible to find a rotating system in which the solution is a periodic one. The angular velocity of this rotating frame is calculated as function of the masses of the two satellites.
2. Including the resonant terms and assuming an exact commensurability of the implied frequencies, we demonstrate that the condition for periodic solutions in the rotating system as defined in (a) is: the initial position of the satellites at conjunction lies on an axis defined by (Ω₁+Ω₂)/2 or (Ω₁+Ω₂)/2 + π/2;Ω₁ and Ω₂ are the longitudes of the ascending nodes of the satellite's orbits. The solution still is a periodic one, thus all the conjunction occur in either axis.
3. In the Mimas Tethys case there is only approximately commensurability between these frequencies. The two satellites are considered as oscillators whose amplitudes and phases are functions of time. The equation of the libration can be established; we find the usual form, but for each satellite the generating solution is a periodic solution (as defined in Paper I), but not a Keplerian one. It follows a determination of the masses which slightly differs from that given by Kozai (1957), when the same values of the observed quantities are used for calculations.
4. The equation of the libration is: $$\ddot z + n_1^2 h^2 \sin z + n_1 q\dot z\sin z = 0$$

     

General class of Bianchi cosmological models with Λ in creation-field cosmology

The solutions of Einstein’s equations with cosmological constant (Λ) in the presence of a creation field have been obtained for general class of anisotropic cosmological models. We have obtained the cosmological solutions for two different scenarios of average scale factor. In first case, we have discussed three different types of physically viable cosmological solutions of average scale factor for the general class of Bianchi cosmological models by using a special law for deceleration parameter which is linear in time with a negative slope. In second case, we have discussed another three different forms of cosmological solutions by using the average scale factor in three different scenarios like Intermediate scenario, Logamediate scenario and Emergent scenario. All physical parameters are calculated and discussed in each physical viable cosmological model. We examine the nature of creation field and cosmological constant is dominated the early Universe but they do not survive for long time and finally tends to zero for large cosmic time t. We have also discussed the all energy conditions in each cases.     

Exploding radiating viscous spheres in general relativity

The influence of viscosity on the gravitational collapse in radiating fluid spheres is investigated. The interior solution is matched with the Vaidya metric at the boundary of the fluid distribution. Prescribing an equation of state to take into account the degree of induced anisotropy by the viscosity and using the Herrera, Jiménez and Ruggeri method, we obtain an explicit Tolman VI-like exploding model. The sphere explodes more violently when the anisotropy due to the viscosity is smaller. The shear viscosity diminishes with the expansion of the distribution of matter.     

Dilaton stabilization in string cosmology. I

A new mechanism is proposed for stabilization of the scalar dilaton field within the framework of lowenergy string gravitation with loop corrections to the dilaton coupling functions. It is based on the assumption that the loop corrections generate a kinetic dilaton function, which is singular for some finite value of the dilaton field. For a nongravitational source of the barotropic type, the system of equations describing the evolution of homogeneous and isotropic cosmological models is represented in the form of a thirdorder, autonomous, dynamical system. The behavior of the general solution in the vicinity of singularities of the dilaton coupling function is investigated by methods of the qualitative theory of dynamical systems. It is shown that there is a class of solutions, different from solutions of the general theory of relativity, with a constant dilaton. The conditions under which these solutions are an attractor for a general solution with a variable dilaton are determined. Translated from Astrofizika, Vol. 43, No. 1, pp. 123-136, January–March, 2000.