Computing derivatives of a gravity potential by using automatic differentiation

A new method, based on automatic differentiation technique, has been proposed in this paper to compute the derivatives of the gravity potential. Using this method we can obtain derivatives up to any order. Instead of explicit expressions of the derivatives we use an iterative scheme to simultaneously compute the value of all the desired derivatives. The algorithm here presented can be easily parallelized by using OpenMP with the consequent improvement in CPU-time efficiency.     

Optimizing expansions of the Earth’s gravity potential and its derivatives over ellipsoidal harmonics

A set of spherical harmonics is the most widely used representation of the Earth’s gravity potential. This series converges outside and on the surface of a reference sphere enveloping the Earth. However, the Earth’s surface is better approximated by the reference ellipsoid—a compressed ellipsoid of revolution that covers the entire Earth. The gravity potential can be expanded in a series over ellipsoidal harmonics on the surface of the reference ellipsoid and on the surface of other external confocal ellipsoids of revolution. In contrast to spherical harmonics, depending on the associated Legendre functions of the first kind, ellipsoidal harmonics depend also on the associated Legendre functions of the second kind. The latter contain the very slowly converging hypergeometric Gauss series. The number of series increases with increasing the order of their derivatives. In this work, we derived new series for the gravitational potential of the Earth and its derivatives over ellipsoidal harmonics. Starting from the first order derivative, all the series corresponding to higher order derivatives depend on the same two hypergeometric Gauss series. The latter converges considerably faster than that for the hypergeometric series previously used when computing the gravity potential and its derivatives.     

Deep versus shallow origin of gravity anomalies, topography and volcanism on Earth, Venus and Mars

The relation between gravity anomalies, topography and volcanism can yield important insights about the internal dynamics of planets. From the power spectra of gravity and topography on Earth, Venus and Mars we infer that gravity anomalies have likely predominantly sources below the lithosphere up to about spherical harmonic degree l=30 for Earth, 40 for Venus and 5 for Mars. To interpret the low-degree part of the gravity spectrum in terms of possible sublithospheric density anomalies we derive radial mantle viscosity profiles consistent with mineral physics. For these viscosity profiles we then compute gravity and topography kernels, which indicate how much gravity anomaly and how much topography is caused by a density anomaly at a given depth. With these kernels, we firstly compute an expected gravity-topography ratio. Good agreement with the observed ratio indicates that for Venus, in contrast to Earth and Mars, long-wavelength topography is largely dynamically supported from the sublithospheric mantle. Secondly, we combine an empirical power spectrum of density anomalies inferred from seismic tomography in Earth’s mantle with gravity kernels to model the gravity power spectrum. We find a good match between modeled and observed gravity power spectrum for all three planets, except for 2?l?4 on Venus. Density anomalies in the Venusian mantle for these low degrees thus appear to be very small. We combine gravity kernels and the gravity field to derive radially averaged density anomaly models for the Martian and Venusian mantles. Gravity kernels for l?5 are very small on Venus below ≈800 km depth. Thus our inferences on Venusian mantle density are basically restricted to the upper 800 km. On Mars, gravity anomalies for 2?l?5 may originate from density anomalies anywhere within its mantle. For Mars as for Earth, inferred density anomalies are dominated by l=2 structure, but we cannot infer whether there are features in the lowermost mantle of Mars that correspond to Earth’s Large Low Shear Velocity Provinces (LLSVPs). We find that volcanism on Mars tends to occur primarily in regions above inferred low mantle density, but our model cannot distinguish whether or not there is a Martian analog for the finding that Earth’s Large Igneous Provinces mainly originate above the margins of LLSVPs.     

Temporal Variations of Perturbed Elliptic Elements: A Semi-Analytical Approach

The classical treatment of implied differences on the orbital ellipticelements from the errors involved at an initial epoch is not possible toapply if we consider a long interval of integration, because there is atemporal variation for all the partial derivatives of the elements withrespect to all the variations in the initial ones. We propose asemi-analytical method to compute these partial derivatives by solving a setof initial value problems which are obtained from the planetary Lagrangeequations and their partial derivatives with respect to all the variationsin the initial elements.     

Mean rates of the orbital elements of a satellite perturbed by a lens shaped mass concentration

Long arc gravity analysis of lunar orbiter tracking data in the past has been carried out with the help of averaged equations of motion, in which short period effects have been suppressed. This procedure has required that the harmonic terms in the gravity potential be averaged over an orbital period. In the present paper, we extend this technique to mass points and mass discs in the gravity field. This required the evaluation of expressions for the mean rates of the orbit elements for a satellite perturbed by a lens shaped mass concentration. Corresponding expressions for the perturbations due to a mass point are obtained in the limit as the lens radius goes to zero. The derived equations have been programmed on the UNIVAC 1108 computer, and the results checked by numerical differencing.     

Small body surface gravity fields via spherical harmonic expansions

Conventional gravity field expressions are derived from Laplace’s equation, the result being the spherical harmonic gravity field. This gravity field is said to be the exterior spherical harmonic gravity field, as its convergence region is outside the Brillouin (i.e., circumscribing) sphere of the body. In contrast, there exists its counterpart called the interior spherical harmonic gravity field for which the convergence region lies within the interior Brillouin sphere that is not the same as the exterior Brillouin sphere. Thus, the exterior spherical harmonic gravity field cannot model the gravitation within the exterior Brillouin sphere except in some special cases, and the interior spherical harmonic gravity field cannot model the gravitation outside the interior Brillouin sphere. In this paper, we will discuss two types of other spherical harmonic gravity fields that bridge the null space of the exterior/interior gravity field expressions by solving Poisson’s equation. These two gravity fields are obtained by assuming the form of Helmholtz’s equation to Poisson’s equation. This method renders the gravitational potentials as functions of spherical Bessel functions and spherical harmonic coefficients. We refer to these gravity fields as the interior/exterior spherical Bessel gravity fields and study their characteristics. The interior spherical Bessel gravity field is investigated in detail for proximity operation purposes around small primitive bodies. Particularly, we apply the theory to asteroids Bennu (formerly 1999 RQ36) and Castalia to quantify its performance around both nearly spheroidal and contact-binary asteroids, respectively. Furthermore, comparisons between the exterior gravity field, interior gravity field, interior spherical Bessel gravity field, and polyhedral gravity field are made and recommendations are given in order to aid planning of proximity operations for future small body missions.     

Axisymmetric smoothed particle hydrodynamics with self-gravity

The axisymmetric form of the hydrodynamic equations within the smoothed particle hydrodynamics (SPH) formalism is presented and checked using idealized scenarios taken from astrophysics (free fall collapse, implosion and further pulsation of a Sun-like star), gas dynamics (wall heating problem, collision of two streams of gas) and inertial confinement fusion (ablative implosion of a small capsule). New material concerning the standard SPH formalism is given. That includes the numerical handling of those mass points which move close to the singularity axis, more accurate expressions for the artificial viscosity and the heat conduction term and an easy way to incorporate self-gravity in the simulations. The algorithm developed to compute gravity does not rely in any sort of grid, leading to a numerical scheme totally compatible with the Lagrangian nature of the SPH equations.     

Emergent universe in Hor?ava gravity

Recently Hor̆ava has proposed a non-relativistic renormalisable gravity theory with higher spatial derivatives in four dimensions which reduces to Einstein’s gravity at large distances with a non-vanishing cosmological constant but with improved UV behaviour. In this paper, we have considered the Friedman-Lema?tre-Robertson-Walker cosmological model in Hor̆ava gravity and the emergent scenario for all values of the spatial curvature k (=0,±1) has been studied. As a result, there are constraints on the parameters involved.     

Numerical computation of incomplete elliptic integrals of a general form

We present an algorithm to compute the incomplete elliptic integral of a general form. The algorithm efficiently evaluates some linear combinations of incomplete elliptic integrals of all kinds to a high precision. Some numerical examples are given as illustrations. This enables us to numerically calculate the values and the partial derivatives of incomplete elliptic integrals of all kinds, which are essential when dealing with many problems in celestial mechanics, including the analytic solution of the torque-free rotational motion of a rigid body around its barycenter.     

Exterior gravitation of a polyhedron derived and compared with harmonic and mascon gravitation representations of asteroid 4769 Castalia

The exterior gravitation of a constant-density polyhedron is derived analytically in closed form. Expressions for potential, attraction, and gravity gradient matrix involve one logarithm term per edge and one arctangent term per face, The Laplacian can be used to determine whether a field point is inside or outside the polyhedron, This polyhedral method is well suited to evaluating the gravitational field of an irregularly shaped body such as an asteroid or comet, Conventional harmonic and mascon potential and attraction expressions suffer large errors when evaluated close to a polyhedral model of asteroid 4769 Castalia.     

Nonsingular expansions of the gravity potential and its derivatives at satellite altitudes in the ellipsoidal coordinate system

The series in ellipsoidal harmonics for derivatives of the Earth’s gravity potential are used only on the reference ellipsoid enveloping the Earth due to their very complex mathematical structure. In the current study, the series in ellipsoidal harmonics are constructed for first- and second-order derivatives of the potential at satellite altitudes; their structure is similar to the series on the reference ellipsoid. The point P is chosen at a random satellite altitude; then, the ellipsoid of revolution is described, which passes through this point and is confocal to the reference ellipsoid. An object-centered coordinate system with the origin at the point P is considered. Using a sequence of transformations, the nonsingular series in ellipsoidal harmonics is constructed for first and second derivatives of the potential in the object-centered coordinate system. These series can be applied to develop a model of the Earth’s potential, based on combined use of surface gravitational force measurements, data on the satellite orbital position, its acceleration, or measurements of the gravitational force gradients of the first and second order. The technique is applicable to any other planet of the Solar System.     

Partial derivatives of eclipse functions αoc and αtr

Expressions are given for partial derivatives of eclipse functions with respect to geometrical depth,p, and the ratio of radii,k. The derivatives are evaluated for critical combinations ofp andk at which indeterminacies occur and the resulting expressions are listed. All expressions are given in a form suitable for numerical evaluation. Notation employed is that of Merrill.     

The roche coordinates and their use in hydrodynamics or celestial mechanics

The aim of the present paper will be to introduce a new system of curvilinear coordinateshereafter referred to as Roche coordinates-in which spheres of constant radius are replaced by equipotential surfaces of a rotating gravitational dipole (which consists of two discrete points of finite mass, revolving around their common center of gravity); while the remaining coordinates are orthogonal to the equipotentials. It will be shown that the use of such coordinates offers a new method of approach to the solution of certain problems of particle dynamics (such as, for instance, the construction of certain types of trajectories in the restricted problem of three bodies); as well as of the hydrodynamics of gas streams in close binary systems, in which the equipotential surfaces of their components distorted by axial rotation and mutual tidal interaction constitute essential boundary conditions.Following a general outline of the problem in Section 1, the Roche coordinates associated with the equipotentials of a rotating gravitational dipole will be constructed in the plane case (Section 2), and their geometrical properties discussed. In Section 3, we shall transform the fundamental equations of hydrodynamics to their forms appropriate in the curvilinear Roche coordinates. The metric coefficients of this transformation will be formulated in a closed form in Section 4 in terms of the respective partial derivatives of the potential; while in Section 5 analytic expressions for the Roche coordinates will be given in the orbital plane of the dipole, which are exact as far as the distortion of the equipotential curves from circular form can be described by the second, third and, fourth harmonics.The concluding Section 6 will be devoted to a formulation of the equations of a mass-point in the restricted problem of three bodies in the Roche coordinates. Three special cases will be considered: (a) motion in the neighborhood of the equipotential curves; (b) motion in the direction normal to such curves; and (c) motion in the neighbourhood of the Lagrangian points. It will be shown that motion in one coordinate is possible only in limiting cases which will be enumerated; but twodimensional motions in which one velocity component is very much smaller than the other invite further study.A generalization of the plane Roche coordinates to three dimensions, with application to additional classes of problems, is being postponed for a subsequent paper.     

Applying an Expansion of the Second Derivatives of the Earth’s Gravitational Potential in Modified Spherical Harmonics for Constructing and Estimating Its Models

Second-order derivatives of the Earth’s potential in a local north-oriented coordinate system are expanded in series of modified spherical harmonics. Linear relations are derived between the spectral coefficients of these series and the spectrum of the geopotential. Based on these relations, recurrent procedures are developed for estimating the geopotential coefficients from the spectrum of each derivative and, conversely, for simulating the spectrum from a known geopotential model. The very simple structure of the expressions for the derivatives is convenient for estimating the coefficients of the geopotential by the least squares method at a certain step of processing satellite gradiometry data. Since the new series are orthogonal, the method with a quadrature formula can be applied, which helps avoid aliasing errors caused by the truncation of the series. The spectral coefficients of the derivatives are estimated using the derived relations for different models on an average orbital sphere of the GOCE satellite and at other altitudes above the Earth’s surface.

     

Representation of light pressure resultant force and moment as a tensor series

In this article, we address the problem of the determination of light pressure upon space structures with a complex geometric shape. For each surface element, we enforce a condition that it can interact with light only from its front side, a condition represented in the form of series of Chebyshev polynomials of the first kind. This Chebyshev expansion enables the use of a series of tensors of increasing rank for determination of the force and moment acting on the sail. We obtain expressions for the determination of light pressure on space structures of complex geometry, taking into account self-shadowing and reflections within the structure. We also give the expressions for tensor parametrization using the specularity coefficient in case of specular -diffuse reflection. For these expressions, we calculated the principal moment and force upon two-sided flat solar sail, spherical and cylindrical bodies, and approximated light pressure upon the proposed space-based observatory Millimetron. The proposed expressions can be used in the ballistic analysis of solar sails and other space objects significantly affected by radiation pressure. Also, these results can be used to analyze the dynamics of large-scale space structures around their center of gravity under light pressure.     

Tunneling analysis under the influences of Einstein–Gauss–Bonnet black holes gravity theory

We consider an equation of motion for Glashow–Weinberg–Salam model and apply the semiclassical Hamilton–Jacobi process and WKB approximation in order to compute the tunneling probability of W-bosons in the background of electromagnetic field to analyze the quantum gravity effects of charged black hole(BH) in Einstein–Gauss–Bonnet gravity theory. After this, we examine the quantum gravity influences on the generalized Lagrangian field equation. We make clear that quantum gravity effects leave the remnants on the tunneling radiation becomes non-thermal. Moreover, we analyze the graphical behavior of quantum gravity influences on corrected Hawking temperature for spin-1 particles for charged BHs.     

Effects of a deformation of a star on the gravitational lensing

We study analytically a gravitational lens due to a deformed star, which is modelled by using a monopole and a quadrupole moment. Positions of the images are discussed for a source on the principal axis. We present explicit expressions for the lens equation for this gravitational lens as a single real 10th-order algebraic equation. Furthermore, we compute an expression for the caustics as a discriminant for the polynomial. Another simple parametric representation of the caustics is also presented in a more tractable form. A simple expression for the critical curves is obtained to clarify a topological feature of the critical curves; the curves are simply connected if and only if the distortion is sufficiently large.     

Gluon polarization at finite temperature and density

We calculate the gluon self-mass and the QCD coupling constant at finite temperature in the real-time formalism up to the first loop level. The expressions are derived in a form that is valid for all temperature ranges in QCD. Using these results the dynamically generated mass of gluons and the plasma screening length can be determined from their effective potential.     

Effects of inflationary bubbles on the polarization and temperature anisotropies of the cosmic microwave background

The broad X-ray iron line, detected in many active galactic nuclei, is likely to be produced by fluorescence from the X-ray-illuminated central parts of an accretion disc close to a supermassive black hole. The time-averaged shape of the line can be explained most naturally by a combination of special and general relativistic effects. Such line profiles contain information about the black hole spin and the accretion disc, as well as the geometry of the emitting region, and may help to test general relativity in the strong gravity regime. In this paper we embark on the computation of the temporal response of the line to the illuminating flux. Previous studies concentrated on the calculation of reverberation signatures from static sources illuminating the disc. In this paper we focus on the more physically justified case of flares located above the accretion disc and corotating with it. We compute the time-dependent iron line, taking into account all general relativistic effects, and show that its shape is of a very complex nature, and we also present light curves accompanying the iron line variability. We suggest that present and future X-ray satellites like XMM or Constellation-X may be capable of detecting features present in the computed reverberation maps.     

Gravitational collapse in generalized Vaidya space-time for Lovelock gravity theory

In this work, we have assumed the generalized Vaidya solution in Lovelock theory of gravity in (n+2)-dimensions. It has been shown that Gauss-Bonnet gravity, dimensionally continued Lovelock gravity and pure Lovelock gravity can be constructed by suitable choice of parameters. We have investigated the occurrence of singularities formed by the gravitational collapse in above three particular forms of Lovelock theory of gravity. The dependence of the nature of singularity on the existence of radial null geodesic for Vaidya space-time has been specially considered. In all the three models, we have shown that the nature of singularities (naked singularity or black hole) completely depend on the parameters. Choices of various parameters are shown in tabular form. In Gauss-Bonnet gravity theory, it can be concluded that the possibility of naked singularity increases with increase in dimensions. In dimensionally continued Lovelock gravity, the naked singularity is possible for odd dimensions for several values of parameters. In pure Lovelock gravity, only black hole forms due to the gravitational collapse for any values of parameters. It has been shown that when accretion is taking place on a collapsing object, it is highly unlikely to get a black hole. Finally on considering the phantom era in the expanding universe it is observed that there is no possibility of formation of a black hole if we are in the Gauss-Bonnet gravity considering the accreting procedure upon a collapsing object.