Mimicking General Relativity with Newtonian Dynamics

The aim of this paper is twofolded. 1 Showing that Newtonian mechanics of point particles in static potentials admits an alternative description in terms of effective riemannian spacetimes. 2 Using the above geometrization scheme to investigate aspects of the gravitational field as it appears in the Einstein theory. It is shown that the mechanical 3 1 effective metrics are quite similar to Gordon’s metric, as it is suggested by the well-known optical-mechanical analogy. Some special potentials are worked out.


Introduction
The last years have witnessed a great interest in the study of analog gravity 1-14 .The underlying idea is to investigate kinematical aspects of Einstein's theory using physical situations where a kind of geometrization procedure can also be applied.The main motivation to the relativist is the fact that it is possible to assign an effective spacetime to describe propagation features of different physical systems.The generality behind these effective metric techniques is being investigated in diverse contexts such as hydrodynamics, optics inside media, nonlinear field theory, superfluids, and quantum condensates see 15, 16 for a complete review .
Although the first analog models appeared in the early days of general relativity 17 it was only in the eighties that they called the attention of the community as a whole.The main reason for this is Unruh's seminal paper on "experimental black hole evaporation" 18 .Unruh conjectured that analogue models could help us to understand not only classical aspects of GR but also some features of quantum field theory in curved spacetimes.In particular, Unruh showed that an analogue black hole may present some properties of gravitational ones as far as the quantum thermal radiation is concerned.Since then, it became a common practice to use analog gravity as a tool to mimic different aspects of the gravitational theory: artificial black holes, Hawking radiation, emergent gravity, quantum field theory in curved backgrounds, and others.
In this paper, I adopt a different perspective concerning "artificial" gravity.I analyze if effective spacetime techniques are suitable to describe some aspects of Newtonian mechanics of point particles and vice versa.As it is known, the relation between curved geometries and particle dynamics was studied at least since the time of Jacobi 19 .After him, this geometrization scheme was to be the concern of Liouville, Lipschitz, Thomson, Tait, and Hertz see 19, 20 for a detailed discussion .With the advent of tensor calculus, it became clear that there existed a map between the trajectories of certain mechanical systems in configuration space and the geodesics of a curved manifold.This fact was explored to solve mechanical problems of holonomic constraints by Ricci, Levi-Civita and in the relativistic era by Synge 21,22 , Lanczos 23 , Lichnerowicz 24 , and others.Nevertheless, despite the intense use of geometrical techniques in the context of dynamics, it seems that the relation between mechanics and geometry was not clearly appreciated in the literature of analog gravity see, nevertheless, 25-27 .This is, perhaps, because many of the above maps did not take the time coordinate into account in the geometrization scheme.
It will be shown that Newton's mechanics of point particles in static potentials may provide a very simple analog model of gravitation.At first this statement seems to be suspicious once we will obtain an alternative description of Newtonian trajectories in terms of pseudo-Riemannian spacetimes without any reference to relativity theory.Nevertheless, using the well-known optical-mechanical analogy as a starting point we will see that, indeed, it is possible to give an effective spacetime description of the motion instead of using the traditional description in terms of forces.We will see that the resulting analog model works as a counterpart of Gordon's metric in the context of optics, while the trajectory of the particle is mapped into null geodesics of an effective spacetime with metric g μν .Finally, we explore some specific potentials and explicitly calculate the effective metrics.
This new geometrized scenario of mechanics may give us some interesting hints to the study of future models because of three main reasons. 1 It introduces curved spacetimes from a very simple and well understood physics. 2 It may extend our laboratory perspectives to measure effective gravitation using, for instance, electronic optics.3 It may provide an interesting and smooth transition to the issue of quantization since we are working with the mechanics of particles instead of fields.

The Optical-Mechanical Analogy Revisited
The well-known optical-mechanical analogy OMA has been discussed often and from many different points of view 25, 27 .The analogy may be understood as a formal tool which maps mechanical systems into optical ones and vice versa.In this sense it implies that the experience and insights developed in one area may be extrapolated to solve problems in the other.Although the OMA is quite familiar, I briefly discuss its main points to fix notation and to make the paper self-contained.
Let us start by considering a particle with mass "m" in a static potential V x k .We are interested in the subset of particle trajectories with a fixed energy E 0 according to Newton's mechanics.It is instructive to begin with an arbitrary spatial coordinate system x k k 1, 2, 3 and three-dimensional euclidean space with metric g ik x .Also, instead of using the usual time coordinate "t" as a parameter I use the arclength "l" along the curve, where dl 2 g ik ẋi ẋk dt 2 .

2.1
Defining the normalized tangent vector u i to the curve it is straightforward to show that Newton's second law can be cast in the compact form where ";" is defined as the three-dimensional covariant derivative, and the modulus of the momentum is related to the potential by the usual energy relation It is important to note here that, because the modulus "p" is automatically known if x k is given, we just need as initial conditions the position x k 0 and the direction of propagation u k 0 .If we specify only the initial position, the solutions of 2.2 characterize all possible trajectories passing through x k 0 with fixed momentum p x k 0 .A similar situation appears in the context of geometrical optics.According to Fermat's principle see, for instance 28 , the ray equation inside an heterogeneous isotropic material with index of refraction n x k is given by the extrema of the optical length λ ≡ ct λ ndl n g ik ẋi ẋk dt, 2.4 where "c" is the velocity of light in vacuum.By varying the action and again parametrizing the paths with the arclength "l" it follows In the same way, because the velocity of light is automatically given in terms of the index of refraction, one only needs to specify position and direction to obtain a given trajectory.By comparing 2.5 with 2.2 one immediately recognizes the similarity between the equations.Thus, there exists a clear analogy between Newtonian trajectories for point particles in static conservative fields and the light rays inside isotropic heterogeneous media.
To complete the analogy we divide 2.2 by p 0 ≡ 2m|E 0 | to express it in terms of an adimensional quantity p ≡ p/p 0 .Because the index of refraction is also an adimensional quantity, one may define the map where " " stands if E 0 > 0 and "−" if E 0 < 0. In this way, the content of the OMA can be formulated as follows.
i To a given optical trajectory inside an isotropic material with index of refraction "n," it is possible to assign a static potential "V " given by relation 2.6 such that it mimics this trajectory in the context of Newtonian mechanics or vice versa.
It is important to note that the parallelism extends only to the geometrical form of the extremal curve.As was noted by Lanczos 23 , the evolution of the systems in time is not equivalent.We will explore this fact later.

The Relation between Paths and Three-Dimensional Effective Geodesics
It is a well-known fact in optics that Fermat's principle may be mapped in the problem of finding geodesics of a three-dimensional curved manifold 23, 29 .If one interprets the optical length dλ ≡ ndl as an infinitesimal element of arc, it is immediate to see that the extremal of the action 2.4 is completely equivalent to the geodesics of the riemannian geometry given by the metric In other words, the light path is such that, if it is parametrized by the optical length λ it satisfies where the effective Christoffel symbol is such that Although this is a purely formal property of the action principle, one can convince himself/herself by explicitly substituting the identity λ ndl in 2.8 .We obtain, after a straightforward calculation and definition 2.9 , the equation

2.10
By calculating explicitly the first term one finally obtains which is no more than the ray equation 2.5 .The identification trajectories/geodesics appears often in analog models of gravitation based on optics and is currently being explored in diverse contexts such as in the architecture of metamaterials, cloaking devices, negative refraction structures, and perfect magnifying lenses see 30-36 and references therein .
Less discussed in the literature of analogue gravity is the geometrical version of the mechanical counterpart of Fermat's principle.Nevertheless, the similarity between 2.2 and 2.5 is not fortuitous.Newton's equation as well is given in terms of an action principle of the form 2.4 which is Jacobi's principle.One only needs to replace n by p in 2.4 to obtain 2.2 .Thus, there exists also in mechanics a geometrical interpretation of the motion as described in its configuration space.The trajectories 2.2 may be alternatively obtained from the geodesics of the effective riemannian space given by the metric 2.12 Note however that, although the four-dimensional generalization of 2.7 made by Gordon 17 became the paradigm of analog gravity, the same did not happen with its mechanical version.The aim of the next section is to explore this last possibility.

Gordon's Metric and Newtonian Trajectories
Gordon was the first to develop effective metric techniques in the context of analog models.He was interested in trying to describe dielectric media by an effective metric while at the same time using the gravitational field in Einstein ansatz to mimic a dielectric medium.Gordon showed that the trajectories of light inside a dielectric medium was such that they could be mapped in null geodesics of a four-dimensional pseudo-Riemannian metric given by where v μ is a normalized timelike vector with respect to Minkowski's background metric, that is, η μν v μ v ν 1 and represents the motion of the dielectric medium.In the simplest case where the medium is at rest in laboratory's frame, we have v μ δ μ 0 and Gordon's metric acquire the simpler form Now, we are going to show that there exists an entirely equivalent situation in Newtonian mechanics of point paticles.We will see that Gordon's metric provides the apparatus to describe Newtonian motions by means of the optical-mechanical analogy.Although it is possible to give a formal demonstration of our proposition we are going to adopt a simpler route by showing that Newton's second law in the form 2.2 is completely equivalent to a null geodesic on the four-dimensional curved spacetime given by the metric 3.2 .
We first note that the energy relation 2.3 can be put in the suggestive form where the first term has the index of refraction form given by 2.6 and p i ≡ dx i /dλ.Our next aim is to write this equation as a quadratic form in a four-dimensional riemannian manifold.We thus define an effective spacetime M with coordinates x μ ≡ λ, x 1 , x 2 , x 3 with λ given by 2.4 and 2.6 note that the coordinate λ has the dimensions of length .
We are going to parametrize the trajectory of the particle in M in terms of λ, that is, x μ λ .Defining the tangent four-vector p μ as we see that 3.3 is equivalent to the expression g μν p μ p ν 0, 3.5 with the effective metric g μν given by the Gordon metric 3.2 and the index of refraction given by the optical-mechanical relation 2.6 .Relation 3.5 is the mechanical analogue of the light cone in relativity.Note that, for the effective metric to make sense it has to carry a hyperbolic Lorentzian signature.This is in complete agreement with the fact that the modulus of the three-dimensional momentum "p" is real, and thus the potential is such that V x ≤ E 0 .We thus succeeded to write the energy relation of Newtonian particles with fixed energy in terms of a "dispersion relation" in a pseudo-Riemannian geometry.Nevertheless, contrarily to our relativistic intuition, the trajectories of massive particles are not mapped into timelike geodesics, but are tangent to effective null curves.
To conclude our four-dimensional geometrization we are going to show that, indeed, the null geodesics of Gordon's metric coincide with Newton's equation in the form 2.2 .To make our calculations simpler, we note that null geodesics are insensitive to conformal transformations.We thus make the conformal map by the substitution to obtain a new manifold with metric given by where "−" stands if E 0 > 0 and " " if E 0 < 0. The four-dimensional geodesics of this geometry is given by We first note that, because the potential V x k depends only on position by hypothesis we have 3.9

ISRN Mathematical Physics 7
We thus obtain that, because x 0 ≡ λ, the "0" component of the equation is identically satisfied.Also, because we have Γ k μ0 0, the spatial components are immediately reduced to an equation of the form 2.8 , which is identically satisfied as a consequence of Jacobi principle one can convince herself/himself of this fact by explicitly substituting dλ pvdt in the geodesic equation to obtain Newton's equations .We thus arrive at the following conclusion.
The trajectories of particles with given energy E 0 in Newton's mechanics of static potentials are such that they can be derived from the null geodesics of a curved four-dimensional spacetime given by Gordon's metric up to a conformal transformation.
We thus define the effective four-dimensional element of distance as

3.10
Trajectories are simply given by the condition d s 2 0. Note that the effective geometry is dependent on the total energy of the system.This means that we have to deal with different geometries to describe motions with different energies.Note also that, in general, the mechanical Riemann tensor R αβμν associated to the effective metric will be nonzero.At this point it is appropriate to use some words of Barcelo et al. in 15 : "It is quite remarkable that even though the underlying" particle "dynamics is Newtonian, nonrelativistic, and takes place in flat space plus time", the trajectories "are governed by a curved 3 1 dimensional Lorentzian pseudo-Riemannian spacetime geometry.For practitioners of general relativity this observation describes a very simple and concrete physical model for certain classes of Lorentzian spacetimes, including black holes.It is also potentially of interest to practitioners of particle "mechanics in that it provides a simple concrete introduction to Lorentzian differential geometric techniques."

Effective Geometries from Newtonian Potentials
We now turn to some explicit situations where our geometrization procedure may be applied.Our aim is to show that some mechanical properties of the systems may as well be described in terms of an effective geometrical language.

The Geometry of Harmonic Oscillations
We start with a one-dimensional oscillator with potential energy V kx 2 /2.The 1 1 dimensional effective line element 3.10 is of the form where "a" is the amplitude of the oscillator.Note that this metric is hyperbolic only if |x| < a.Thus, this effective spacetime is such that there exists two geometrical boundaries given by the conditions x ±a.This is an immediate consequence of the fact that we are considering only trajectories such that the total energy E 0 is fixed, and thus the motion is not defined anywhere.
The null trajectories may be obtained by imposing the condition d s 0, that is,

4.2
In this simple case, it is possible to integrate this equation to obtain the implicit relation A diagram for the trajectories is given in Figure 1.The set of all possible null curves x λ given by the implicit relation 4.3 characterizes the space-time diagram for the oscillator.Note however that, although the curves exhibit a periodic behavior in terms of λ the derivative dx/dλ is not well defined in the limit x → ±a.In this limit the effective volume − g dx dλ 1 − x/a dx dλ also vanishes, and one can interpret this as a consequence of the fact that the effective spacetime only makes sense inside the boundaries given by x ±a.
At a given point of the diagram there exist only two admissible trajectories with fixed energy E 0 .These two curves determine at each point P a conoid structure in the same way as the propagation of light characterizes the light cone in relativity.We see from the figure that the intersection between the curves gives the mechanical analogue of the light cone at each spacetime point.The domain of causality of the point P is represented by the shaded regions.The curves 4.3 are null geodesics with respect to the effective metric.Nevertheless, a careful inspection in 4.1 shows that this metric is flat.Thus, metric 4.1 is simply Minkowski metric in a curvilinear coordinate system.We are now going to investigate a more complicated situation where the Riemann tensor does not vanish.

The Geometry of Newtonian Gravitation
We now turn to the motion of point test particles with mass "m" in the Newtonian gravitational potential.As in the previous section, we are interested in the motions associated to a given total energy E 0 .According to our previous discussion, the null geodesics of the respective effective metric 3.10 are such that they reproduce ordinary Newtonian trajectories in the gravitational potential.
To be specific, let us consider the spherically symmetric potential V r −GMm/r.We will consider only the case of bounded orbits, that is, E 0 < 0. The other situations will be considered elsewhere.

E 0 < 0 Effective Metrics
In this case we have a maximum admissible radii r H ≡ −GMm/E 0 for the trajectories.One obtains for the effective metric, in spherical coordinates, the expression Analogously to the harmonic oscillator case, we obtain an effective metric that admits a hyperbolic Lorentzian signature only in a particular domain.In the present case this domain is given by the condition 0 < r < r H .This is because outside this region, that is, for r > r H , the movement is not allowed for the considered energies E 0 −GMm/r H .This fact has an interesting geometrical interpretation in our scheme.First, let us note that the square root of the determinant is given by

4.5
This relation implies that the element of effective volume dV eff ≡ − g d 4 x vanishes for the values r 0 and r r H .This is a hint that the confinement of the motion between these values of the radial coordinate has a deep geometrical meaning.
In fact, the components of the Riemann tensor are rr H sin 2 θ 3r H − 4r r H − r .

ISRN Mathematical Physics
One immediately obtains for the Ricci scalar the expression This object diverges for the values r 0 and r r H note that R < 0 .Thus, the effective 3 1 spacetime is delimited by two singularities in the effective geometry.The first one is quite similar to Schwarzchild singularity and admits a similar interpretation due to the strength of the gravitational field.It appears because the Newtonian gravitational potential also diverges at the origin.The second singularity has, nevertheless, a different and unexpected origin.
It appears because we confined our attention only to the movement of particles with a fixed energy E 0 .As we mentioned before, all possible trajectories with this energy cannot scape from the region r < r H .This fact appears in the geometrical description as a kind of geometrical barrier that do not allow the null geodesics to escape from the allowed region.
To develop some additional feeling on our geometrical description of dynamics, it is instructive to solve the well-known problem of finding the trajectories in the Newtonian gravitational potential in terms of the effective null geodesics 3.8 .We recall that the following quantity is a constant of the motion: From the geodesic equation we obtain for the sake of conciseness I use dx/dλ ≡ ẋ r A r ṙ2 B r θ2 sin 2 θ φ2 0, 4.9 θ C r ṙ θ − sin θ cos θ φ2 0, 4.10 φ C r φ ṙ 2cotgθ φ θ 0, 4.11 where

4.12
As usual, we assume that the orbit lies in the plane θ π/2.Thus, 4.10 ensures that, if θ π/2 and θ 0 initially, the orbit remains always in this plane.Equation 4.11 may be integrated directly.We obtain that the following quantity is also conserved where 0 ≤ < 1 because r ≤ r H always. Thus, we succeeded to find the orbits around the field of a sphericall object with mass "M" in terms of the null geodesics of the effective metric 4.4 .

The Geometry of Electrostatic Potentials
The conservative potentials that appear in electromagnetism are, perhaps, the most interesting from the point of view of analog gravity in our context.This is not only because electromagnetic fields are simple to handle in laboratory sizes, but because many interesting configurations of static potentials occur in the advanced field of electronic optics.I will concentrate here on simple electrostatic configurations.The potential is obtained by the familiar formula while the potential energy is V x qφ.The effective metric is given by the expression where q is the charge of the test charge.Our next task is to calculate this metric for a specific case and to use some of its geometrical properties to infer about mechanical properties of the system.

The Case of the Electric Dipole
One interesting situation appears in the geometrization of the static electric dipole.If the dipole moment p is oriented in the z direction, we obtain, for large distances, φ r, θ p 4π 0 cos θ r 2 .

4.22
For the sake of simplicity lets us consider pq > 0 and positive energies E 0 > 0. The effective line element reads where κ 2 ≡ pq/4π 0 E 0 .Note that this effective spacetime only makes sense if This is an immediate consequence of the fact that the energy E 0 is fixed from the beginning and there are regions that are not allowed for the particle motion.In our scheme, this fact has an interesting geometrical interpretation.Note that the determinant of the metric vanishes for some values of the radii.At first we could suspect that the coordinate system we are using is not well defined in this region.Nevertheless, we obtain for the scalar of curvature the expression which diverges in the abovementioned region.Thus, the description of charged particle motions such as electrons or ions in terms of null geodesics of an effective spacetime forces us to deal with a singular geometry in the case of the electric dipole.Nevertheless, the singularity is not so simple as it appears in Schwarzschild geometry, but it is determined by a compact bidimensional region given by the condition r 2 κ 2 cos θ .

4.27
The null geodesics are not defined in the interior region.Note, however, that the entire region outside the singularity is well defined.Nevertheless, the geodesics are not so simple to be integrated as in the case of the solar system.Also, it is immediate to see that this geometry is asymptotically flat.Thus, for very large values of r, where the field of the dipole tends to vanish, the trajectories are described by the null geodesics of Minkowski spacetime.As a final remark, I would like to point out that the geometrical framework may introduce techniques of general relativity and riemannian geometry in the realm of electron optics.We suspect that it is possible to envisage situations where electronic dispositives may be projected using these new tools.This last statement is strongly suggested by the recent achievements of transformation optics 31-35 .

Conclusions
The relation between geometry and dynamics is an interesting problem that has been investigated from many different perspectives along the years.In this paper it is shown that Newtonian mechanics of point particles in static potentials admits a geometrical interpretation in terms of 3 1 effective spacetimes in a similar way as it appears in the context of analog gravity models.The mechanical effective metric g μν is given explicitly in terms of the potential energy V x and the total energy E 0 and is similar to Gordon's metric in optics.From this last similarity, we inferred that mechanical systems may provide a very simple arena for the study of analogue gravity.This new arena may bring interesting theoretical and experimental challenges.I investigated some explicit examples, where the use of geometrical techniques were used to infer about some mechanical properties.Thus, it may be natural to speak about effective horizons, singularities and other structures that appear in general relativity in mechanical problems also.In particular it was considered the motion of particles in the gravitational field of a central mass and in the potential of an electric dipole.The null geodesics of the above geometries coincides with the usual trajectories as they're described in the Newtonian ansatz.It may be interesting to investigate if these geometrical picture can be sustained after first quantization.
r H r φ ξ,4.13where ξ is a constant.It is straightforward to show that this equation implies the usual conservation of angular momentum.In fact, from 4.13 , it follows, after a reparametrization that