A new class of bulk viscous universe with time dependent deceleration parameter and Λ-term

A new class of exact solutions of Einstein’s field equations with a bulk viscous fluid for an LRS Bianchi type-Ia obtained by using a time dependent deceleration parameter and cosmological term Λ. The coefficient of bulk viscosity is assumed to be a power function of mass density (ξ=ξ ₀ ρ ⁿ ). We have obtained a general solution of the field equations from which six models of the universe are derived: exponential, polynomial and sinusoidal form respectively. The behaviour of these models of the universe are also discussed in the frame of reference of recent supernovae Ia observations.      

LRS Bianchi-V cosmology with decaying vacuum density and heat flow

In this paper we consider a locally-rotationally-symmetric (LRS) Bianchi type-V perfect fluid model with variable cosmological ‘constant’ representing the energy density of vacuum. The field equations are solved with and without heat conduction by using a variation law for the mean Hubble parameter, which is related to the average scale factor of the metric and yields a constant value of the deceleration parameter. A constant value of deceleration parameter generates power-law form of average scale factor which is used to find the exact solutions with and without heat conduction with decaying vacuum density. The solutions presented here satisfy all the necessary conditions for the physically acceptability. The thermodynamical relations in decaying vacuum fluid model are also studied in detail.     

Anisotropic Bianchi type-I models with constant deceleration parameter in general relativity

A special law of variation for Hubble’s parameter is presented in a spatially homogeneous and anisotropic Bianchi type-I space-time that yields a constant value of deceleration parameter. Using the law of variation for Hubble’s parameter, exact solutions of Einstein’s field equations are obtained for Bianchi-I space-time filled with perfect fluid in two different cases where the universe exhibits power-law and exponential expansion. It is found that the solutions are consistent with the recent observations of type Ia supernovae. A detailed study of physical and kinematical properties of the models is carried out.     

Bianchi type-V cosmological model with perfect fluid using negative constant deceleration parameter in a scalar tensor theory based on Lyra Manifold

An exact Bianchi type-V perfect fluid cosmological model is obtained in a scalar tensor theory proposed by Sen (Z. Phys. 149:311, 1957) based on Lyra Manifold in case of β is a constant and it is shown that this cosmological model exists only in the case of Radiation Universe (ρ=3p) if β is a function of ‘t’ using negative constant deceleration parameter. Some physical and geometrical properties of these models are discussed.     

Early Cosmological Models with Bulk Viscosity in Brans-Dicke Theory

The effect of time dependent bulk viscosity on the evolution of Friedmann models with zero curvature in Brans-Dicke theory is studied. The solutions of the field equations with ‘gamma-law’ equation of state p = (γ-1) ρ, where γ varies continuously as the Universe expands, are obtained by using the power-law relation φ = bR ⁿ , which lead to models with constant deceleration parameter. We obtain solutions for the inflationary period and radiation dominated era of the universe. The physical properties of cosmological solutions are also discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Posterior probability of the deceleration parameter q0 from quasars provided with a luminosity indicator

We have worked out a ’statistical algorithm’ for obtaining the posterior probability density of the deceleration parameter q₀ from quasars where there is a luminosity indicator available. We point out that the role of the luminosity indicator is to provide asecond estimate of individual luminosities after a first estimate has been obtained from measured brightness and redshift together with an assumed q₀. Discrimination of q₀ is to be sought in the statistical properties of the set of differences between the two estimates (the residuals). We show that the variance of the residuals and their correlation with redshifts (further refined to luminosity distances) are two independent test-statistics for q₀, whose known distributions then lead to the probability density sought. We have applied the above algorithm to a sample of flat-spectrum radio quasars with measured CIV, MgII and Ly α lines. A combined Baldwin’s relation was used for all 3 lines. Our result is that log q₀ is normally distributed with a mean value of + 0.270± 0.135 (s.d.), or, q₀ = + 1.86 ± 0.135 dex. This result, we believe, is the sharpest result so far published on q₀.     

Vertical stability of periodic solutions around the triangular equilibrium points

The vertical stability character of the families of short and long period solutions around the triangular equilibrium points of the restricted three-body problem is examined. For three values of the mass parameter less than equal to the critical value of Routh (μ _R ) i.e. for μ = 0.000953875 (Sun-Jupiter), μ = 0.01215 (Earth-Moon) and μ = μ _R = 0.038521, it is found that all such solutions are vertically stable. For μ > (μ _R ) vertical stability is studied for a number of ‘limiting’ orbits extended to μ = 0.45. The last limiting orbit computed by Deprit for μ = 0.044 is continued to a family of periodic orbits into which the well known families of long and short period solutions merge. The stability characteristics of this family are also studied.     

Analytic properties of Hansen coefficients

Hansen’s coefficients in the theory of elliptic motion with eccentricity e are studied as functions of the parameter η = (1 − e ²)^1/2. Their analytic behavior in the complex η plane is described and some symmetry relations are derived. In particular, for every Hansen coefficient, multiplication by suitable powers of e and η results in an entire analytic function of η. Consequently, Hansen’s coefficients can be in principle computed by means of rapidly convergent series in powers of η. A representation of Hansen’s coefficients in terms of two entire functions of e ² follows.      

Unified structure of analytical solutions for radiative transfer problems in plane-parallel,homogeneous media

The basic concepts for developing a system of analytic solutions for the standard problems of radiative transfer theory are discussed. These solutions, which are found using Ambartsumyan’s layer addition method in Sobolev’s probabilistic interpretation for radiative diffusion problems, are maximally compact and easily used in numerical computations. New expressions are obtained for the resolvents and the resolvent functions, as well as a unified structure for the form of an integral representation for solving different radiative transfer problems in semi-infinite media and in finite layers. Block diagrams of the sequence of stages for solving these problems are provided, where the Ambartsumyan function φ(η) (more precisely, 1/φ(η)) plays a fundamental role in the case of semi-infinite media while the functions a(η, τ₀ ) and b(η, τ₀) play an analogous role for finite layers.     

Accelerating universe with time variation of G and Λ

We study a gravitational model in which scale transformations play the key role in obtaining dynamical G and Λ. We take a non-scale invariant gravitational action with a cosmological constant and a gravitational coupling constant. Then, by a scale transformation, through a dilaton field, we obtain a new action containing cosmological and gravitational coupling terms which are dynamically dependent on the dilaton field with Higgs type potential. The vacuum expectation value of this dilaton field, through spontaneous symmetry breaking on the basis of anthropic principle, determines the time variations of G and Λ. The relevance of these time variations to the current acceleration of the universe, coincidence problem, Mach’s cosmological coincidence and those problems of standard cosmology addressed by inflationary models, are discussed. The current acceleration of the universe is shown to be a result of phase transition from radiation toward matter dominated eras. No real coincidence problem between matter and vacuum energy densities exists in this model and this apparent coincidence together with Mach’s cosmological coincidence are shown to be simple consequences of a new kind of scale factor dependence of the energy momentum density as ρ∼a ⁻⁴. This model also provides the possibility for a super fast expansion of the scale factor at very early universe by introducing exotic type matter like cosmic strings.     

Cosmological vacuum solutions in Brans and Dicke's scalar-tensor theory

The exact cosmological vacuum solutions of Brans and Dicke's scalas-tensor theory are derived when a power law is valid between the gravitational constant κ and the radius of curvatureR of the universe. There exist even in the case of the closed 3-dimensional space of positive curvature solutions with increasingR and κ with respect to the age of the universe, whereby the freely available parameter ω of the scalar-tensor theory can take all values greater than −3/2. Such solutions are contrary to Dirac's hypothesis as well as to Einstein-Mach's principle.     

Families of Asymmetric Periodic Orbits in Hill’s Problem of Three Bodies

We present five families of periodic solutions of Hill’s problem which are asymmetric with respect to the horizontal ξ axis. In one of these families, the orbits are symmetric with respect to the vertical η axis; in the four others, the orbits are without any symmetry. Each family consists of two branches, which are mirror images of each other with respect to the ξ axis. These two branches are joined at a maximum of Γ, where the family of asymmetric periodic solutions intersects a family of symmetric (with respect to the ξ axis) periodic solutions. Both branches can be continued into second species families for Γ → − ∞.     

Evolution of the General Solution of the Restricted Problem Covering Symmetric and Escape Solutions

The work presented in paper I (Papadakis, K.E., Goudas, C.L.: Astrophys. Space Sci. (2006)) is expanded here to cover the evolution of the approximate general solution of the restricted problem covering symmetric and escape solutions for values of μ in the interval [0, 0.5]. The work is purely numerical, although the available rich theoretical background permits the assertions that most of the theoretical issues related to the numerical treatment of the problem are known. The prime objective of this work is to apply the ‘Last Geometric Theorem of Poincaré’ (Birkhoff, G.D.: Trans. Amer. Math. Soc. 14, 14 (1913); Poincaré, H.: Rend. Cir. Mat. Palermo 33, 375 (1912)) and compute dense sets of axisymmetric periodic family curves covering the initial conditions space of bounded motions for a discrete set of values of the basic parameter μ spread along the entire interval of permissible values. The results obtained for each value of μ, tested for completeness, constitute an approximation of the general solution of the problem related to symmetric motions. The approximate general solution of the same problem related to asymmetric solutions, also computable by application of the same theorem (Poincaré-Birkhoff) is left for a future paper. A secondary objective is identification-computation of the compact space of escape motions of the problem also for selected values of the mass parameter μ. We first present the approximate general solution for the integrable case μ = 0 and then the approximate solution for the nonintegrable case μ = 10⁻³. We then proceed to presenting the approximate general solutions for the cases μ = 0.1, 0.2, 0.3, 0.4, and 0.5, in all cases building them in four phases, namely, presenting for each value of μ, first all family curves of symmetric periodic solutions that re-enter after 1 oscillation, then adding to it successively, the family curves that re-enter after 2 to 10 oscillations, after 11 to 30 oscillations, after 31 to 50 oscillations and, finally, after 51 to 100 oscillations. We identify in these solutions, considered as functions of the mass parameter μ, and at μ = 0 two failures of continuity, namely: 1. Integrals of motion, exempting the energy one, cease to exist for any infinitesimal positive value of μ. 2. Appearance of a split into two separate sub-domains in the originally (for μ = 0) unique space of bounded motions. The computed approximations of the general solution for all values of μ appear to fulfill the ‘completeness’ criterion inside properly selected sub-domains of the domain of bounded motions in the (x, C) plane, which means that these sub-domains are filled countably densely by periodic family curves, which form a laminar flow-line pattern. The family curves in this pattern may, or may not, be intersected by a ‘basic’ family curve segment of order from 1 up to 3. The isolated points generating asymptotic solutions resemble ‘sink’ points toward which dense sets of periodic family curves spiral. The points in the compact domain in the (x, C) plane resting outside the domain of bounded motions (μ = 0), including the gap between the two large sub-domains (μ > 0) created by the aforementioned split, generate escape motions. The gap between the two large sub-domains of bounded motions grows wider for growing μ. Also, a number of compact gaps that generate escape motions exist within the body of the two sub-domains of bounded motions. The approximate general solutions computed include symmetric, heteroclinic, asymptotic, collision and escape solutions, thus constituting one component of the full approximate general solution of the problem, the second and final component being that of asymmetric solutions.     

Some new exact solutions with finite central parameters and?uniform radial motion of sound

We present three new categories of exact and spherically symmetric Solutions with finite central parameters of the general relativistic field equations. Two well behaved solutions in curvature coordinates first category are being studied extensively. These solutions describe perfect fluid balls with positively finite central pressure, positively finite central density; their ratio is less than one and causality condition is obeyed at the centre. The outmarch of pressure, density, pressure-density ratio and the adiabatic speed of sound is monotonically decreasing for these solutions. Keeping in view of well behaved nature of these solutions, one of the solution (I₁) is studied extensively. The solution (I₁) gives us wide range of Schwarzschild parameter u (0.138≤u≤0.263), for which the solution is well behaved hence, suitable for modeling of Neutron star. For this solution the mass of Neutron star is maximized with all degree of suitability and by assuming the surface density ρ _b =2×10¹⁴ g/cm³. Corresponding to u=0.263, the maximum mass of Neutron star comes out to be 3.369 M _Θ with linear dimension 37.77 km and central and surface redshifts are 4.858 and 0.4524 respectively. We also study some well known regular solutions (T-4, D-1, D-2, H, A, P) of Einstein’s field equations in curvature coordinates with the feature of constant adiabatic sound speed. We have chosen those values of Schwarzschild parameter u for which, these solutions describe perfect fluid balls realistic equations of state. However, except (P) solution, all these solutions have monotonically non-decreasing feature of adiabatic sound speed. Hence (P) solution is having a well behaved model for uniform radial motion of sound. Keeping in view of well behaved nature of the solution for this feature and assuming the surface density; ρ _b =2×10¹⁴ g/cm³, the maximum mass of Neutron star comes out to be 1.34 M _Θ with linear dimension 28.74 km. Corresponding central and surface redshifts are 1.002 and 0.1752 respectively.     

The general solution of the HENON–HEILES problem

The general solution of the Henon–Heiles system is approximated inside a domain of the (x, C) of initial conditions (C is the energy constant). The method applied is that described by Poincaré as ‘the only “crack” permitting penetration into the non-integrable problems’ and involves calculation of a dense set of families of periodic solutions that covers the solution space of the problem. In the case of the Henon–Heiles potential we calculated the families of periodic solutions that re-enter after 1–108 oscillations. The density of the set of such families is defined by a pre-assigned parameter ε (Poincaré parameter), which ascertains that at least one periodic solution is computed and available within a distance ε from any point of the domain (x, C) for which the approximate general solution computed. The approximate general solution presented here corresponds to ε = 0.07. The same solution is further improved by “zooming” into four square sub-domain of (x, C), i.e. by computing sufficient number of families that reduce the density parameter to ε = 0.003. Further zooming to reduce the density parameter, say to ε = 10⁻⁶, or even smaller, although easily performable in both areas occupied by stable as well as unstable solutions, was found unnecessary. The stability of all members of each and all families computed was calculated and presented in this paper for both the large solution domain and for the sub-domains. The correspondence between areas of the approximate general solution occupied by stable periodic solutions and Poincaré sections with well-aligned section points and also correspondence between areas occupied by unstable solutions and Poincaré sections with randomly scattered section points is shown by calculating such sections. All calculations were performed using the Runge-Kutta (R-K) 8th order direct integration method and the large output received, consisting of many thousands of families is saved as “Atlas of the General Solution of the Henon–Heiles Problem,” including their stability and is available at request. It is concluded that approximation of the general solution of this system is straightforward and that the chaotic character of its Poincaré sections imposes no limitations or difficulties.     

A new class of LRS Bianchi type-II dark energy models with variable EoS parameter

A new class of dark energy models in a Locally Rotationally Symmetric Bianchi type-II (LRS B-II) space-time with variable equation of state (EoS) parameter and constant deceleration parameter have been investigated in the present paper. The Einstein’s field equations have been solved by applying a variation law for generalized Hubble’s parameter given by Berman: Nuovo Cimento 74:182 (1983) which generates two types of solutions for the average scale factor, one is of power-law type and other is of the exponential-law form. Using these two forms, Einstein’s field equations are solved separately that correspond to expanding singular and non-singular models of the universe respectively. The dark energy EoS parameter ω is found to be time dependent and its existing range for both models is in good agreement with the three recent observations of (i) SNe Ia data (Knop et al.: Astrophys. J. 598:102 (2003)), (ii) SNe Ia data collaborated with CMBR anisotropy and galaxy clustering statistics (Tegmark et al.: Astrophys. J. 606:702 (2004)) and latest (iii) a combination of cosmological datasets coming from CMB anisotropies, luminosity distances of high redshift type Ia supernovae and galaxy clustering (Hinshaw et al.: Astrophys. J. Suppl. 180:225 (2009); Komatsu et al. Astrophys. J. Suppl. 180:330 (2009)). The cosmological constant Λ is found to be a positive decreasing function of time and it approaches a small positive value at late time (i.e. the present epoch) which is corroborated by results from recent supernovae Ia observations. The physical and geometric behaviour of the universe have also been discussed in detail.     

Quantum vacuum fluctuations and dark energy

It is shown that the curvature of space-time induced by vacuum fluctuations of quantum fields should be proportional to the square of Newton’s constant G. This offers a possible explanation for the success of the formula ρ∼ Gm ⁶ c ² h ⁻⁴, ρ being the dark energy density and m a typical mass of elementary particles.     

The case for the Universe to be a quantum black hole

We present a necessary and sufficient condition for an object of any mass m to be a quantum black hole (q.b.h.): “The product of the cosmological constant Λ and the Planck’s constant ℏ, Λ and ℏ corresponding to the scale defined by this q.b.h., must be of order one in a certain universal system of units”. In this system the numerical values known for Λ are of order one in cosmology and about 10¹²² for Planck’s scale. Proving that in this system the value of the cosmological ℏ _c is of order one, while the value of ℏ for the Planck’s scale is about 10⁻¹²², both scales satisfy the condition to be a q.b.h., i.e. Λℏ≈1. In this sense the Universe is a q.b.h. We suggest that these objects, being q.b.h.’s, give us the linkage between thermodynamics, quantum mechanics, electromagnetism and general relativity, at least for the scale of a closed Universe and for the Planck’s scale. A mathematical transformation may refer these scales as corresponding to infinity (our universe) and zero (Planck’s universe), in a “scale relativity” sense.     

Effect of 3rd-degree gravity harmonics and Earth perturbations on lunar artificial satellite orbits

In a previous work we studied the effects of (I) the J ₂ and C ₂₂ terms of the lunar potential and (II) the rotation of the primary on the critical inclination orbits of artificial satellites. Here, we show that, when 3rd-degree gravity harmonics are taken into account, the long-term orbital behavior and stability are strongly affected, especially for a non-rotating central body, where chaotic or collision orbits dominate the phase space. In the rotating case these phenomena are strongly weakened and the motion is mostly regular. When the averaged effect of the Earth’s perturbation is added, chaotic regions appear again for some inclination ranges. These are more important for higher values of semi-major axes. We compute the main families of periodic orbits, which are shown to emanate from the ‘frozen eccentricity’ and ‘critical inclination’ solutions of the axisymmetric problem (‘J ₂ + J ₃’). Although the geometrical properties of the orbits are not preserved, we find that the variations in e, I and g can be quite small, so that they can be of practical importance to mission planning.     

Restricted Three-Body Problem: An Approximation of its General Solution – Part One – the Manifold of Symmetric Periodic Solutions

This paper gives the results of a programme attempting to exploit ‘la seule bréche’ (Poincaré, 1892, p. 82) of non-integrable systems, namely to develop an approximate general solution for the three out of its four component-solutions of the planar restricted three-body problem. This is accomplished by computing a large number of families of ‘solutions précieuses’ (periodic solutions) covering densely the space of initial conditions of this problem. More specifically, we calculated numerically and only for μ = 0.4, all families of symmetric periodic solutions (1st component of the general solution) existing in the domain D:(x ₀ ∊ [−2,2],C ∊ [−2,5]) of the (x ₀, C) space and consisting of symmetric solutions re-entering after 1 up to 50 revolutions (see graph in Fig. 4). Then we tested the parts of the domain D that is void of such families and established that they belong to the category of escape motions (2nd component of the general solution). The approximation of the 3rd component (asymmetric solutions) we shall present in a future publication. The 4th component of the general solution of the problem, namely the one consisting of the bounded non-periodic solutions, is considered as approximated by those of the 1st or the 2nd component on account of the `Last Geometric Theorem of Poincaré' (Birkhoff, 1913). The results obtained provoked interest to repeat the same work inside the larger closed domain D:(x ₀ ∊ [−6,2], C ∊ [−5,5]) and the results are presented in Fig. 15. A test run of the programme developed led to reproduction of the results presented by Hénon (1965) with better accuracy and many additional families not included in the sited paper. Pointer directions construed from the main body of results led to the definition of useful concepts of the basic family of order n, n = 1, 2,… and the completeness criterion of the solution inside a compact sub-domain of the (x ₀, C) space. The same results inspired the ‘partition theorem’, which conjectures the possibility of partitioning an initial conditions domain D into a finite set of sub-domains D i that fulfill the completeness criterion and allow complete approximation of the general solution of this problem by computing a relatively small number of family curves. The numerical results of this project include a large number of families that were computed in detail covering their natural termination, the morphology, and stability of their member solutions. Zooming into sub-domains of D permitted clear presentation of the families of symmetric solutions contained in them. Such zooming was made for various values of the parameter N, which defines the re-entrance revolutions number, which was selected to be from 50 to 500. The areas generating escape solutions have being investigated. In Appendix A we present families of symmetric solutions terminating at asymptotic solutions, and in Appendix B the morphology of large period symmetric solutions though examples of orbits that re-enter after from 8 to 500 revolutions. The paper concludes that approximations of the general solution of the planar restricted problem is possible and presents such approximations, only for some sub-domains that fulfill the completeness criterion, on the basis of sufficiently large number of families.