Curvature sign convention Created 2026-10-14 Updated 2026-10-25
Definitions of the Riemann curvature tensor can differ by an overall sign, independent of the metric signature. Reversing that definition reverses the Ricci tensor, Ricci scalar and Einstein tensor, as well as the curvature side of the geodesic deviation identity and the Palatini identity. A derivation must choose the contraction and derivative ordering explicitly before combining formulas from different conventions. Changing a metric signature alone is a separate operation.
Deviation vector 2026-10-07
For a geodesic variation , the deviation vector describes infinitesimal separation at equal parameter values. With , commuting parameter derivatives give . A torsion-free connection then gives . For an affinely parametrized variation, geodesic deviation states with the convention .
Einstein field equations Created 2026-09-24 Updated 2026-10-18
With signature and vanishing cosmological constant, the Einstein field equations areIn four spacetime dimensions, contraction with gives the useful trace relation .
In the usual commutator convention for the Riemann curvature tensor, : is the gravitational constant, whereas is the Einstein tensor. The Ricci tensor and scalar curvature encode contractions of curvature. The stress-energy tensor includes energy density, momentum flow and stresses, so pressure also gravitates. A cosmological constant adds to the left-hand side. See Tong's Einstein equations and matter coupling.
For a concrete consequence, set and . The trace gives , and substitution gives . Yet the exterior Schwarzschild spacetime still has the tidal forces computed under geodesic deviation: vanishing Ricci tensor does not imply vanishing Riemann curvature tensor. Geometry carries gravitational effects through empty space. Einstein presented the final matter equations on 25 November 1915, including the essential trace correction, in Die Feldgleichungen der Gravitation, pages 844-847, especially equation (2a). His symbols and sign conventions differ from modern ones; the displayed formula uses modern notation.
In a parallel-propagated orthonormal frame, orthogonal geodesic deviation in constant sectional curvature reduces to three scalar equations . For , : generic initial data grow exponentially, but the pure decaying mode is a counterexample to universal growth. Initially comoving neighbors have . For , , so the deviation is bounded and oscillatory. Simultaneous focusing requires compatible vector initial data and is not automatic for every oscillatory solution.
Flat-space anisotropic beam shear focusing 2026-10-14
In flat spacetime take a twist-free beam's two screen Jacobi factors to be and . They solve the zero-curvature geodesic deviation equation. Their product is the displayed area, and subtracting the trace from their optical matrix gives the displayed nonnegative shear. Direct differentiation yields , an explicit check of the focusing sign.
General relativity Created 2026-09-24 Updated 2026-10-13
General relativity describes gravity through the geometry of spacetime. Its metric tensor specifies distances, proper time measured by clocks, and the causal structure that limits communication. Matter and radiation affect this geometry through the Einstein field equations; freely falling test particles follow geodesics of the resulting geometry. A planet's orbit and a falling apple therefore belong to the same geometric story. The Newtonian limit recovers familiar Newtonian gravity when fields are weak and speeds are small compared with the speed of light. See David Tong's introduction to general relativity.
The most direct way to understand spacetime curvature is to compare neighboring freely falling bodies. Even when each body feels weightless, their separation can change: this is geodesic deviation, observed as a tidal force. Curvature concerns these relative measurements, rather than the shape of a drawing in an imagined surrounding space. A second accessible measurement is gravitational time dilation: identical clocks maintained at different heights accumulate different proper times. Together, tides and clocks turn an abstract geometry into laboratory questions.
At first order, geodesic deviation gives for a small detector in a plane gravitational wave in linearized gravity. For physical initial displacement zero and initial velocity , integration gives . A bounded periodic response selects . Prescribed positions alone do not select that response when the wave has no wave-free past.
Local tidal acceleration scaling Created 2026-10-14 Updated 2026-10-25
The geodesic deviation equation makes the relative gravitational acceleration proportional to the separation and a typical Riemann curvature tensor component . Thus tidal acceleration vanishes as at fixed bounded curvature. A uniform electric field can give different accelerations at the same point, independently of separation. Relative tidal changes accumulated over a fixed long time need not vanish as a fraction of the initial separation; a local spacetime experiment must control its duration as well as its spatial size.
Newtonian tidal curvature 2026-10-14
Using , the Newton–Cartan gravitational connection has only the displayed Hessian components and their antisymmetric partners in the last two indices. The geodesic deviation equation gives relative acceleration . The Poisson equation for Newtonian gravity is therefore . A spatially uniform gravitational acceleration has zero tidal curvature even if the connection coefficients do not vanish in the original coordinates.
Use the curvature sign conventionwhich has the derivative ordering printed in the source for . A Christoffel symbol is of order , so its derivative is of order , whereas a product of two connections is of order . Their ratio is . The metric is stationary, hence the temporal derivative term vanishes. ConsequentlyHere , and . The weak-field tidal curvature in a stationary metric is
Printed sign inconsistency. The source's positive half-Hessian and its final negative-prefactor expression are the opposites of this result. A concrete check at is , giving , while the printed final expression gives . Thus both printed equalities cannot hold with its stated derivative ordering. Reversing the definition of the Riemann curvature tensor would give the printed final sign, but would also require reversing the curvature term in the geodesic deviation equation. The convention used here and in the next part gives the physical radial stretching and transverse compression.
Geodesic deviation compares neighboring members of a family of freely falling trajectories. Subtracting their geodesic equations to first order in the connecting vector gives a covariant relative-acceleration equation governed by the Riemann curvature tensor. For a slowly moving satellite, dominates the spatial components. Choose simultaneous spatial separation in a locally transported frame, take the weak-field approximation, and neglect terms quadratic in the small separation and higher velocity corrections. The leading spatial equation becomesReplacing the transported covariant derivatives by ordinary derivatives is a leading local or Newtonian approximation; a rotating orbital frame would introduce additional inertial terms. In the convention established above, the instantaneous tidal acceleration stretches radial separations and compresses transverse ones.
The growth expected for nearby different-radius circular orbits is especially transparent from Kepler's third law. If the radii are and , their coordinate angular frequencies differ by . Initially aligned orbits therefore acquire an along-track separation . Over the interval in which the nearby-orbit approximation remains valid, this produces a growing separation without needing to solve the deviation equation.
This growth is not universal or indefinitely unbounded. Equal-radius co-orbiting geodesics with a fixed phase offset have constant separation, and slightly inclined equal-radius circular orbits have bounded, periodically varying separation. Even different-radius circular orbits have bounded exact positions; their relative angle eventually wraps around. The secular drift statement assumes distinct radii and describes the local differential-orbit regime.
Physical tidal stretching and the early differential drift of neighboring circular orbits
. Details. Take the real physical wave and write . The specified amplitude has zero trace, so its metric perturbation equals its trace-reversed metric perturbation. It is in transverse-traceless gauge:to first order. Assume the freely falling particles are initially at rest in these coordinates. Since and in coordinates ,The geodesic equations therefore admitUniqueness preserves these initial conditions. Thus the coordinate locations remain fixed; this conclusion does not apply to arbitrary initially moving particles.
Coordinate position is not measured separation. In the central observer's local orthonormal frame in spacetime, and in the long-wavelength limit appropriate to a small detector, the transverse physical separations areto first order. In particular the two principal semiaxes of the ring areFor a point at angular position the radial displacement isThe freely falling ring in transverse-traceless gauge alternates between ellipses elongated in the and directions. Its area is unchanged to first order. This is the plus gravitational wave polarization; the axes exchange their roles after half a period.
An observer at the center detects these changing distances, for example by exchanging light signals with the particles. Equivalently geodesic deviation gives the local tidal accelerationto the same order. There is physical relative motion and nonzero Riemann curvature tensor, even though the selected coordinate positions are constant. For a detector comparable to a wavelength, its light-travel-time response must be included rather than treating all distances as simultaneously measured local lengths.
The paper uses , the negative of a common torsion tensor convention. Keep that sign throughout. Also, denotes the covariant Hessian, so the second derivative acts on the lower index of the covector ; it is not just a composition of two directional derivatives on a scalar.
Since the covariant derivative of a scalar field is its ordinary derivative,Commuting the ordinary partial derivatives therefore givesFor a vector field, write . This is a tensor with one covariant and one contravariant index. HenceSubtracting the expression with exchanged, the first two terms combine toThe coefficients of first derivatives of in those two terms cancel. The last terms give . Identifying the curvature coefficient yieldsHere the bracket denotes antisymmetrization of the tensorial second covariant derivative, consistently with the preceding calculation.
To check tensoriality directly, antisymmetrize the supplied transformation law of the connection coefficients in the two lower indices. The inhomogeneous term vanishes because mixed coordinate second derivatives are symmetric. The remaining transformation law iswhich is precisely that of a tensor. Intrinsically, with the paper's sign,The derivative terms involving the components of cancel against their Lie bracket, leaving . This also shows that the expression is linear over smooth functions in each argument.
Finally take the torsion-free connection and the commuting parameter fields along the geodesic family. Vanishing torsion gives , while affine geodesic parametrization gives . The curvature operator in the specified convention isThereforeIn components this is the geodesic deviation equationThe connecting field represents infinitesimal separation within the family, so this is a linear equation for relative free-fall acceleration, rather than an exact coordinate difference for widely separated worldlines.
Displacement must be specified relative to an observer. In transverse-traceless gauge the free particles can have fixed coordinates, while their measured proper lengths oscillate. In a local, freely falling orthonormal frame, the geodesic deviation equation for an initial separation is, to first order in ,Its general solution with physical displacement and initial physical velocity isThe source has been running indefinitely, so there is no wave-free initial state. Its initial conditions for an eternal monochromatic gravitational wave need more than positions to fix . The usual bounded, periodic response selects ; if the particles instead start at physical rest at , set and retain the linear term. Both solve the same linear geodesic deviation problem. The distinction between fixed TT coordinates and measured separation is explained in [ref].
Choose so the wave phase at the detector is . For the bounded response the four initial placements then giveThese are displacements from the specified positions at ; positions are . Equivalently, relative to mean positions the usual oscillatory displacement is to first order. The response repeats after .
Both plus polarization and cross polarization occur with equal amplitude and a quarter-cycle phase difference. Their superposition is circular gravitational-wave polarization, with handedness reversed when the viewing direction or source rotation is reversed. The principal stretching axis turns through half the wave phase. There is no leading longitudinal motion.
One full circularly polarized gravitational-wave cycle, with the four particles at their specified initial positions at phase zero
. Details. Past exam of the general relativity II course of the University of Oxford 2017 1 c Solution 2026-10-25
The geodesic deviation equation compares neighboring geodesics by covariant derivatives. Replacing those derivatives by ordinary derivatives of coordinate components requires a parallel transport frame along the reference geodesic. There is also an affine parameter normalization to specify: the argument of is , and it becomes only after choosing .
For the intended local tidal calculation, take the reference timelike geodesicHere , so is proper time. Both and its first partial derivatives vanish on . All Christoffel symbols therefore vanish along this curve, and its coordinate basis is parallel transport along it. With the Riemann curvature tensor from part (b),The Jacobi field equation thus becomes, in this adapted frame,A purely transverse initial Jacobi field, with initially, remains transverse. This transverse tidal equation of a Brinkmann plane wave describes the relative tidal acceleration of the freely falling bodies.
In vacuum write the symmetric matrix asIts instantaneous eigenvalues are . The tidal acceleration stretches along one principal axis and squeezes along the perpendicular axis, with the axes allowed to vary with . This describes accelerations; it does not imply monotonically growing and shrinking separations for arbitrary initial velocities. For the simple constant profile , ,Thus the geodesic deviation can be defocusing in one direction and oscillatory in the other. For a varying profile, diagonalizing at each instant does not decouple the ordinary differential equations without also including derivatives of the rotating eigenvectors.
Qualification of the coordinate claim. The displayed equation for is not true for arbitrary coordinate displacements along arbitrary reference geodesics. For synchronized neighboring curves with and , the exact geodesic equations areTheir variation givesThe second expression need not vanish. A concrete vacuum counterexample is and the one-parameter family of timelike geodesicsEach has . Differentiating with respect to produces a Jacobi fieldThis even has , so it is a physical orthogonal separation. The covariant derivative form of geodesic deviation remains correct; the unrestricted assertion about ordinary coordinate derivatives needs the adapted reference and frame used above. If instead is a constant affine parameter normalization, the transverse tidal coefficient is .
Past exam of the mathematical physics course of the University of Oxford 2007 B 7 Solution 2026-10-31
The coordinate derivatives of the smooth two-parameter map commute. With a torsion-free Levi-Civita connection,because is symmetric in its lower indices. Each member of the family is a geodesic, so . ConsequentlyThe commutator of covariant derivatives is the action of the Riemann curvature tensor. With the source's index and sign convention this givesThis is geodesic deviation: free fall relative acceleration measures curvature, even where an individual observer can eliminate the connection at one point.
For Earth's exterior Newtonian field, and acceleration is . Nearby freely falling particles have relative acceleration . At the observation point align radially upward and horizontally. Since ,The weak-field geodesic-deviation equation therefore identifiesIn these units time is measured in seconds as in the stated Newtonian comparison; if the zeroth coordinate is instead, the curvature components acquire the corresponding factor. A small freely falling cloud is squeezed in the two horizontal directions and stretched radially. These are tidal effects, not the uniform acceleration removed by a freely falling frame.
Their characteristic inverse-square time is , so . One concrete thought experiment is an ideal satellite moving in a circular orbit just above a spherical Earth, ignoring atmosphere and terrain. Its orbital angular frequency is and its period is minutes. Equivalently, initially nearby free fall trajectories develop order-one tidal changes on this characteristic scale. This is the surface tidal time scale of a spherical gravitating body.
Past exam of the mathematics course of the University of Cambridge 2001 iii Paper 68 2 Solution Created 2026-10-03 Updated 2026-10-14
Use , background Minkowski spacetime metric tensor , and with small metric perturbation. All indices and derivatives in this linear approximation use . The linearized inverse metric is , and the linearized Levi-Civita connection isDefine the trace-reversed metric perturbation , with .
An infinitesimal coordinate change gives the linearized coordinate gauge transformationThe metric perturbation therefore has redundant components. The gauge invariance of the linearized Riemann tensor follows because its change contains third derivatives of cancelling by commutation of partial derivatives. Conversely, flat linearized metric perturbations are locally pure gauge on a contractible patch. Thus linearized curvature, rather than a particular coordinate value of , detects the physical disturbance. Gauge freedom can also be restricted by global topology or boundary conditions.
Keeping the paper's curvature sign convention, its linearized Ricci tensor and scalar areThese are the negatives of the alternative frequently used curvature convention. Put . The Linearized Einstein equations in vacuum areUnder a gauge change , so solving imposes the Lorenz gauge in linearized gravity. The field equations then become . The remaining residual gauge symmetry of linearized gravity has ; Lorenz gauge does not exhaust the coordinate freedom.
A Fourier mode obeys and . Hence the plane gravitational waves in linearized gravity propagate on the background light cones. Four transversality conditions and four residual gauge amplitudes leave physical degrees of freedom. The explicit plane-wave reduction to transverse-traceless gauge can set , spatial trace zero and . For a wave in the direction, write its physical spatial strain as , so the spatial metric tensor is :There are two gravitational wave polarizations, plus polarization and cross polarization. Rotating transverse axes through rotates their amplitude pair through . Nonconstant profiles with nonzero second derivative have nonzero linearized Riemann curvature tensor and cannot be removed as pure gauge.
The TT coordinate separation and measured separation distinction explains a detector's response. Freely falling particles initially at rest can keep fixed TT coordinates, since , while their proper separation changes. For a short arm along a unit vector ,This is the tidal response described by geodesic deviation; coordinate motion by itself is not the measured signal.
Finally, gravitational-wave energy is second order in the perturbation, so it is absent from a first-order vacuum equation. In a short-wavelength averaging regime the averaged stress-energy of transverse gravitational waves isFor a plane wave its energy flux is . The averaging scale and weak-field assumptions matter: this does not assign a coordinate-independent local gravitational energy density to an arbitrary first-order perturbation.
Past exam of the mathematics course of the University of Cambridge 2001 iii Paper 69 5 Solution Created 2026-10-03 Updated 2026-10-14
Interpret first as an infinitesimal screen area transported along generators of a future event horizon, not as the total area of an arbitrary trapped surface. Horizon generators form an affinely parametrized, hypersurface-orthogonal null geodesic congruence. The two-dimensional screen metric is positive definite. Its optical tensor has zero twist and decomposes asThe scalar here agrees with the supplied derivative definition: nullness and the affine geodesic equation remove the longitudinal terms in the contracted norm, leaving .
To derive the focusing equation, vary neighbouring geodesics and let be their screen Jacobi map. The geodesic deviation equation is , where . Hence satisfies . Taking its screen trace yields the Null Raychaudhuri equationSet . Since , the area-square-root optical focusing equation isThe PDF's plus sign before is inconsistent with its own positive shear definition. This is an actual source error, not an antisymmetrization or curvature-sign convention.
For an explicit check, take a flat-space twist-free beam with screen Jacobi factors and , for . ThenIts Ricci term is zero, so the result equals and has the opposite sign to the printed equation. This is flat-space anisotropic beam shear focusing.
Under the null convergence condition , the correct equation gives concavity of each local area square root. Equivalently . If a horizon generator had , integration gives a focal point within affine distance at most . Such a generator could not remain on an achronal boundary beyond the focal point. Under the usual global hypotheses that horizon generators remain regular and future complete, this contradicts their being generators of the event horizon. Hence and . Integrating over the horizon patches, and including newly joining generators which can add area, proves Hawking's area theorem. Global predictability/completeness assumptions are essential; the local differential inequality alone is not a theorem about arbitrary trapped-surface areas.
For the perfect fluid stress tensor, null contraction eliminates the pressure-metric term:A nonzero null vector cannot be orthogonal to a timelike fluid velocity. Thus the perfect-fluid null energy condition is preciselySeparate requirements or are not needed for this null focusing argument. With the Einstein equation, .
A cosmological constant contributes only a term proportional to , whose null contraction is zero:Therefore a nonzero cosmological constant does not change the fluid condition or the local area-focusing argument, though it can change the global asymptotics and which horizons and completeness hypotheses are appropriate.
Past exam of the mathematics course of the University of Cambridge 2002 iii Paper 72 2 Solution Created 2026-10-03 Updated 2026-10-14
The Pound-Rebka experiment measures a small gravitational redshift by comparing nuclear emission and resonant absorption at different heights in Earth's field. The transition of iron-57 supplies the gamma-ray frequency standard. The Mössbauer effect permits recoil-free emission and absorption by nuclei bound in solids, giving a sufficiently narrow resonance. A source and absorber separated vertically by about have a predicted fractional frequency difference of order . A calibrated Doppler effect from slow source motion measures the shift through its effect on resonant absorption. Interchanging the upper and lower arrangements reverses the gravitational contribution and helps separate it from intrinsic source-absorber offsets; temperature-dependent frequency shifts must also be controlled. The observed change agrees with the weak-field prediction.
The prediction follows either from an accelerating laboratory and the equivalence principle, or directly from a static metric. With , write . The time-translation Killing vector field gives a conserved photon energy , while a static observer has four-velocity and measures . ThusFor an upward ray over height , ; downward propagation has the opposite sign. The equivalent compensating Doppler speed has magnitude , about for this baseline. Identical clocks at rest have , so a clock at higher potential runs faster relative to the common static coordinate time. This does not mean that the nuclear transition changes its locally measured frequency with height: the comparison involves transporting a photon between different local proper-time standards.
The weak equivalence principle asserts universal free fall for test bodies. The Einstein equivalence principle adds local Lorentz invariance and local position invariance for nongravitational experiments in a freely falling laboratory. The redshift comparison tests the clock aspect of local position invariance and supports a metric description of gravity. The strong equivalence principle extends these requirements to self-gravitating bodies and local gravitational experiments, including the contribution of gravitational binding energy. Their outcomes should be independent of the external location and motion of a freely falling laboratory when tidal effects are negligible. A nuclear redshift experiment is therefore not a test of the whole strong equivalence principle, and its success does not uniquely establish the gravitational field equations.
The geometrical realization uses the Levi-Civita connection of a spacetime metric. At any event, normal coordinates make and the Christoffel symbols vanish. Freely falling test particles follow geodesics; local nongravitational equations take their special-relativistic form at that event. The Riemann curvature tensor need not vanish, so finite laboratories can measure tidal geodesic deviation. The strong principle concerns removing a uniform external gravitational field locally, not removing curvature throughout an extended region.
For a perfect fluid, local rest-frame isotropy and the absence of heat flux or viscosity give , where is rest-frame energy density. With the PDF's signature , and projects spatially. Its covariant form isThe pressure sign follows from the negative spatial metric; it gives positive diagonal stresses in the rest frame. Local stress-energy conservation, , produces the two projectionsThe first is the energy equation, including pressure work; the second is the relativistic Euler equation. Its rest-frame spatial pressure force is .
Further assumptions are needed to derive the Einstein field equations. Take a local generally covariant metric equation, linear in curvature and of second differential order, and demand compatibility with stress-energy conservation. The available rank-two tensors are , and a constant multiple of . Defining , the contracted Bianchi identity gives . Hence the combination is identically divergence-free only when . This determines the Einstein tensor ; a cosmological constant term remains allowed.
The coupling sign must be matched to the actual curvature convention. Here represents , and a static weak field has . With negligible pressure and cosmological term, a trial equation implies . Comparing with the Poisson equation fixes . Thus the Einstein equation in the sequential-derivative curvature convention isEquivalently, , with . Pressure therefore gravitates as well as supplying a material force. The cosmological constant is not fixed by the local Newtonian matching.
Setting for the action normalization, the same metric equation follows from the Einstein-Hilbert action , where , together with a matter action obeying . After the metric boundary term is removed or supplied separately, the bulk variation is proportional to . This also makes the sign convention explicit. General covariance, the field equation and the contracted Bianchi identity are mutually consistent with the fluid conservation laws. The experiment motivates the local metric description; the curvature and conservation assumptions, plus the Newtonian limit, supply the additional route to the field equation.
Past exam of the mathematics course of the University of Cambridge 2003 iii Paper 54 2 Solution Created 2026-10-03 Updated 2026-10-14
Work to first order in a metric perturbation , with , raising perturbation indices with and discarding quadratic terms. The Christoffel symbols areUsing the printed curvature convention gives the linearized Riemann curvature operatorFor definiteness contract . The corresponding first-order Ricci tensor isThis is the reversed sign convention specified in curvature-sign convention in Killing derivative identities; the vacuum equation is unaffected by an overall curvature sign.
An infinitesimal passive coordinate change changes the perturbation by . Substitution into gives zero change, since partial derivatives commute. Equivalently the change in curvature is the Lie derivative of the background curvature, which vanishes in Minkowski spacetime. Thus linearized curvature, and the linearized Weyl tensor in vacuum, are gauge-invariant tidal observables, whereas individual components of are generally gauge dependent.
The trace-reversed metric perturbation transforms by . Its divergence changes by . Solving therefore imposes Lorenz gauge in linearized gravity, . With the stated sign convention the Einstein tensor isso the linearized vacuum field equation becomes . The residual gauge symmetry of linearized gravity has . This makes a hyperbolic evolution problem after the gauge condition and its initial consistency constraints have been imposed.
To separate actual gauge invariants from a gauge choice, use gauge-invariant Minkowski perturbations. For spatial indices contracted with , decomposewhere , and . Spatial falloff or Fourier-mode conventions fix the decomposition. Write and with . Direct substitution in the gauge change givesHence the invariant combinations areA representative with has , and . The explicit vacuum constraints for gauge-invariant Minkowski perturbations follow fromFor a nonzero spatial Fourier mode, or regular decaying fields on all space, forces , then forces , and the spatial trace gives , hence . Only remains. This removal of the scalar and vector sectors does not exclude singular Coulomb perturbations outside a source, boundary-driven harmonic potentials or spatially homogeneous modes; those require separate boundary and gauge conventions.
For a propagating plane mode , the wave equation gives . A transverse-traceless spatial tensor perpendicular to a fixed propagation direction has two independent components. For propagation along , define the physical strain in this mostly-minus convention, so the spatial line element is . It has the formThese are the two gravitational wave polarizations, plus polarization and cross polarization. A rotation in the transverse plane rotates the polarization pair through twice the angle; cross is plus rotated through degrees. Lorenz constraints and residual gauge freedom give the same count: ten symmetric metric components minus four constraints and four gauge freedoms leave two radiative degrees of freedom.
In transverse-traceless gauge, particles initially at rest keep fixed spatial coordinates since , but their measured separation changes. For a short separation of length in unit direction , . The gauge-invariant curvature gives the corresponding geodesic deviation:Here the last equation is in a freely falling orthonormal frame, with the printed reversed curvature convention. Coordinate separation and physical separation must not be confused. The linearized vacuum equations describe propagation and tidal measurements; wave energy transfer and an effective wave stress-energy tensor arise at quadratic order, beyond the linear approximation.


