# Tidal tensor

Last updated

In Newton's theory of gravitation and in various relativistic classical theories of gravitation, such as general relativity, the tidal tensor represents

## Contents

1. tidal accelerations of a cloud of (electrically neutral, nonspinning) test particles,
2. tidal stresses in a small object immersed in an ambient gravitational field.

The tidal tensor represents the relative acceleration due to gravity of two test masses separated by an infinitesimal distance. The component ${\displaystyle \Phi _{{a}{b}}}$ represents the relative acceleration in the ${\displaystyle {\hat {a}}}$ direction produced a displacement in the ${\displaystyle {\hat {b}}}$ direction.

## Tidal tensor for a spherical body

The most common example of tides is the tidal force around a spherical body (e.g., a planet or a moon). Here we compute the tidal tensor for the gravitational field outside an isolated spherically symmetric massive object. According to Newton's gravitational law, the acceleration a at a distance r from a central mass m is

${\displaystyle a=-Gm/r^{2}}$

(to simplify the math, in the following derivations we use the convention of setting the gravitational constant G to one. To calculate the differential accelerations, the results are to be multiplied by G.)

Let us adopt the frame in polar coordinates for our three-dimensional Euclidean space, and consider infinitesimal displacements in the radial and azimuthal directions, ${\displaystyle \partial _{r},\partial _{\theta },}$ and ${\displaystyle \partial _{\phi }}$, which are given the subscripts 1, 2, and 3 respectively.

${\displaystyle {\vec {\epsilon }}_{1}=\partial _{r},\;{\vec {\epsilon }}_{2}={\frac {1}{r}}\,\partial _{\theta },\;{\vec {\epsilon }}_{3}={\frac {1}{r\sin \theta }}\,\partial _{\phi }}$

We will directly compute each component of the tidal tensor, expressed in this frame. First, compare the gravitational forces on two nearby objects lying on the same radial line at distances from the central body differing by a distance h:

${\displaystyle m/(r+h)^{2}-m/r^{2}=-2mh/r^{3}+3mh^{2}/r^{4}+O(h^{3})}$

Because in discussing tensors we are dealing with multilinear algebra, we retain only first order terms, so ${\displaystyle \Phi _{11}=-2m/r^{3}}$. Since there is no acceleration in the ${\displaystyle \theta }$ or ${\displaystyle \phi }$ direction due to a displacement in the radial direction, the other radial terms are zero: ${\displaystyle \Phi _{12}=\Phi _{13}=0}$.

Similarly, we can compare the gravitational force on two nearby observers lying at the same radius ${\displaystyle r=r_{0}}$ but displaced by an (infinitesimal) distance h in the ${\displaystyle \theta }$ or ${\displaystyle \phi }$ direction. Using some elementary trigonometry and the small angle approximation, we find that the force vectors differ by a vector tangent to the sphere which has magnitude

${\displaystyle {\frac {m}{r_{0}^{2}}}\,\sin(\theta )\approx {\frac {m}{r_{0}^{2}}}\,{\frac {h}{r_{0}}}={\frac {m}{r_{0}^{3}}}\,h}$

By using the small angle approximation, we have ignored all terms of order ${\displaystyle O(h^{2})}$, so the tangential components are ${\displaystyle \Phi _{22}=\Phi _{33}=m/r^{3}}$. Again, since there is no acceleration in the radial direction due to displacements in either of the azimuthal directions, the other azimuthal terms are zero: ${\displaystyle \Phi _{21}=\Phi _{31}=0}$.

Combining this information, we find that the tidal tensor is diagonal with frame components ${\displaystyle \Phi _{{\hat {a}}{\hat {b}}}={\frac {m}{r^{3}}}\operatorname {diag} (-2,1,1)}$ This is the Coulomb form characteristic of spherically symmetric central force fields in Newtonian physics.

## Hessian formulation

In the more general case where the mass is not a single spherically symmetric central object, the tidal tensor can be derived from the gravitational potential ${\displaystyle U}$, which obeys the Poisson equation:

${\displaystyle \Delta U=4\pi \,\mu }$

where ${\displaystyle \mu }$ is the mass density of any matter present, and where ${\displaystyle \Delta }$ is the Laplace operator. Note that this equation implies that in a vacuum solution, the potential is simply a harmonic function.

The tidal tensor is given by the traceless part [1]

${\displaystyle \Phi _{ab}=J_{ab}-{\frac {1}{3}}\,{J^{m}}_{m}\,\eta _{ab}}$

of the Hessian

${\displaystyle J_{ab}={\frac {\partial ^{2}U}{\partial x^{a}\,\partial x^{b}}}}$

where we are using the standard Cartesian chart for E3, with the Euclidean metric tensor

${\displaystyle ds^{2}=dx^{2}+dy^{2}+dz^{2},\;-\infty

Using standard results in vector calculus, this is readily converted to expressions valid in other coordinate charts, such as the polar spherical chart

${\displaystyle ds^{2}=d\rho ^{2}+\rho ^{2}\,\left(d\theta ^{2}+\sin(\theta )^{2}\,d\phi ^{2}\right)}$
${\displaystyle 0<\rho <\infty ,\;0<\theta <\pi ,\;-\pi <\phi <\pi }$

### Spherically symmetric field

As an example, we can calculate the tidal tensor for a spherical body using the Hessian. Next, let us plug the gravitational potential ${\displaystyle U=-m/\rho }$ into the Hessian. We can convert the expression above to one valid in polar spherical coordinates, or we can convert the potential to Cartesian coordinates before plugging in. Adopting the second course, we have ${\displaystyle U=-m/{\sqrt {x^{2}+y^{2}+z^{2}}}}$, which gives

${\displaystyle \Phi _{ab}={\frac {m}{(x^{2}+y^{2}+z^{2})^{5/2}}}\,\left[{\begin{matrix}y^{2}+z^{2}-2x^{2}&-3xy&-3xz\\-3xy&x^{2}+z^{2}-2y^{2}&-3yz\\-3xz&-3yz&x^{2}+y^{2}-2z^{2}\end{matrix}}\right]}$

After a rotation of our frame, which is adapted to the polar spherical coordinates, this expression agrees with our previous result. The easiest way to see this is to set ${\displaystyle y,z}$ to zero so that the off-diagonal terms vanish and ${\displaystyle \rho =x}$, and then invoke the spherical symmetry.

## In General Relativity

In general relativity, the tidal tensor is generalized by the Riemann curvature tensor. In the weak field limit, the tidal tensor is given by the components ${\displaystyle \Phi _{ij}={R}_{i0j0}}$ of the curvature tensor.

## Related Research Articles

In mathematics, a spherical coordinate system is a coordinate system for three-dimensional space where the position of a point is specified by three numbers: the radial distance of that point from a fixed origin, its polar angle measured from a fixed zenith direction, and the azimuthal angle of its orthogonal projection on a reference plane that passes through the origin and is orthogonal to the zenith, measured from a fixed reference direction on that plane. It can be seen as the three-dimensional version of the polar coordinate system.

In mathematics and physics, Laplace's equation is a second-order partial differential equation named after Pierre-Simon Laplace, who first studied its properties. This is often written as

In physics, the Navier–Stokes equations are certain partial differential equations which describe the motion of viscous fluid substances, named after French engineer and physicist Claude-Louis Navier and Anglo-Irish physicist and mathematician George Gabriel Stokes. They were developed over several decades of progressively building the theories, from 1822 (Navier) to 1842–1850 (Stokes).

In mathematics, the Laplace operator or Laplacian is a differential operator given by the divergence of the gradient of a scalar function on Euclidean space. It is usually denoted by the symbols , , or . In a Cartesian coordinate system, the Laplacian is given by the sum of second partial derivatives of the function with respect to each independent variable. In other coordinate systems, such as cylindrical and spherical coordinates, the Laplacian also has a useful form. Informally, the Laplacian Δf (p) of a function f at a point p measures by how much the average value of f over small spheres or balls centered at p deviates from f (p).

Linear elasticity is a mathematical model of how solid objects deform and become internally stressed due to prescribed loading conditions. It is a simplification of the more general nonlinear theory of elasticity and a branch of continuum mechanics.

Stellar dynamics is the branch of astrophysics which describes in a statistical way the collective motions of stars subject to their mutual gravity. The essential difference from celestial mechanics is that the number of body

This is a list of some vector calculus formulae for working with common curvilinear coordinate systems.

In mathematics, a Killing vector field, named after Wilhelm Killing, is a vector field on a Riemannian manifold that preserves the metric. Killing fields are the infinitesimal generators of isometries; that is, flows generated by Killing fields are continuous isometries of the manifold. More simply, the flow generates a symmetry, in the sense that moving each point of an object the same distance in the direction of the Killing vector will not distort distances on the object.

The Kerr–Newman metric is the most general asymptotically flat, stationary solution of the Einstein–Maxwell equations in general relativity that describes the spacetime geometry in the region surrounding an electrically charged, rotating mass. It generalizes the Kerr metric by taking into account the field energy of an electromagnetic field, in addition to describing rotation. It is one of a large number of various different electrovacuum solutions, that is, of solutions to the Einstein–Maxwell equations which account for the field energy of an electromagnetic field. Such solutions do not include any electric charges other than that associated with the gravitational field, and are thus termed vacuum solutions.

The Schwarzschild solution describes spacetime under the influence of a massive, non-rotating, spherically symmetric object. It is considered by some to be one of the simplest and most useful solutions to the Einstein field equations.

In general relativity, a frame field is a set of four pointwise-orthonormal vector fields, one timelike and three spacelike, defined on a Lorentzian manifold that is physically interpreted as a model of spacetime. The timelike unit vector field is often denoted by and the three spacelike unit vector fields by . All tensorial quantities defined on the manifold can be expressed using the frame field and its dual coframe field.

In theoretical physics, Nordström's theory of gravitation was a predecessor of general relativity. Strictly speaking, there were actually two distinct theories proposed by the Finnish theoretical physicist Gunnar Nordström, in 1912 and 1913 respectively. The first was quickly dismissed, but the second became the first known example of a metric theory of gravitation, in which the effects of gravitation are treated entirely in terms of the geometry of a curved spacetime.

In the mathematical description of general relativity, the Boyer–Lindquist coordinates are a generalization of the coordinates used for the metric of a Schwarzschild black hole that can be used to express the metric of a Kerr black hole.

The Newman–Penrose (NP) formalism is a set of notation developed by Ezra T. Newman and Roger Penrose for general relativity (GR). Their notation is an effort to treat general relativity in terms of spinor notation, which introduces complex forms of the usual variables used in GR. The NP formalism is itself a special case of the tetrad formalism, where the tensors of the theory are projected onto a complete vector basis at each point in spacetime. Usually this vector basis is chosen to reflect some symmetry of the spacetime, leading to simplified expressions for physical observables. In the case of the NP formalism, the vector basis chosen is a null tetrad: a set of four null vectors—two real, and a complex-conjugate pair. The two real members asymptotically point radially inward and radially outward, and the formalism is well adapted to treatment of the propagation of radiation in curved spacetime. The Weyl scalars, derived from the Weyl tensor, are often used. In particular, it can be shown that one of these scalars— in the appropriate frame—encodes the outgoing gravitational radiation of an asymptotically flat system.

The Cauchy momentum equation is a vector partial differential equation put forth by Cauchy that describes the non-relativistic momentum transport in any continuum.

In mathematics, vector spherical harmonics (VSH) are an extension of the scalar spherical harmonics for use with vector fields. The components of the VSH are complex-valued functions expressed in the spherical coordinate basis vectors.

In general relativity, a point mass deflects a light ray with impact parameter by an angle approximately equal to

The Jeans equations are a set of partial differential equations that describe the motion of a collection of stars in a gravitational field. The Jeans equations relate the second-order velocity moments to the density and potential of a stellar system for systems without collision. They are analogous to the Euler equations for fluid flow and may be derived from the collisionless Boltzmann equation. The Jeans equations can come in a variety of different forms, depending on the structure of what is being modelled. Most utilization of these equations has been found in simulations with large number of gravitationally bound objects.

In the theory of Lorentzian manifolds, spherically symmetric spacetimes admit a family of nested round spheres. In such a spacetime, a particularly important kind of coordinate chart is the Schwarzschild chart, a kind of polar spherical coordinate chart on a static and spherically symmetric spacetime, which is adapted to these nested round spheres. The defining characteristic of Schwarzschild chart is that the radial coordinate possesses a natural geometric interpretation in terms of the surface area and Gaussian curvature of each sphere. However, radial distances and angles are not accurately represented.

In general relativity, the Weyl metrics are a class of static and axisymmetric solutions to Einstein's field equation. Three members in the renowned Kerr–Newman family solutions, namely the Schwarzschild, nonextremal Reissner–Nordström and extremal Reissner–Nordström metrics, can be identified as Weyl-type metrics.

## References

1. Baldauf, Tobias; Seljak, Uros; Desjacques, Vincent; McDonald, Patrick (13 January 2018). "Evidence for Quadratic Tidal Tensor Bias from the Halo Bispectrum". Physical Review D. 86 (8): 083540. arXiv:. Bibcode:2012PhRvD..86h3540B. doi:10.1103/PhysRevD.86.083540. S2CID   21681130.