ISRN Mathematical PhysicsVolume 2012 (2012), Article ID 260951, 15 pagesdoi:10.5402/2012/260951

Research Article

Mimicking General Relativity with Newtonian Dynamics

E. Goulart

Instituto de Cosmologia Relatividade Astrofisica (ICRA-Brasil/CBPF) Rua Dr. Xavier Sigaud, 150, 22290-180 Rio de Janeiro, RJ, Brazil

Received 4 August 2011; Accepted 21 September 2011

Academic Editors: I. Radinschi, W. Sabra, and D. Singleton

Copyright © 2012 E. Goulart. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

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.

1. 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 . 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.

2. 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 “” in a static potential . We are interested in the subset of particle trajectories with a fixed energy according to Newton’s mechanics. It is instructive to begin with an arbitrary spatial coordinate system () and three-dimensional euclidean space with metric . Also, instead of using the usual time coordinate “” as a parameter I use the arclength “” along the curve, where Defining the normalized tangent vector 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 “” is automatically known if is given, we just need as initial conditions the position and the direction of propagation . If we specify only the initial position, the solutions of (2.2) characterize all possible trajectories passing through with fixed momentum .

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 is given by the extrema of the optical length where “” is the velocity of light in vacuum. By varying the action and again parametrizing the paths with the arclength “” 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 to express it in terms of an adimensional quantity . Because the index of refraction is also an adimensional quantity, one may define the map where “+” stands if and “−” if . 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 “,” it is possible to assign a static potential “” 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.

2.1. 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 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 in (2.8). We obtain, after a straightforward calculation and definition (2.9), the equation 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 by 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 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.

3. 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 is a normalized timelike vector with respect to Minkowski’s background metric, that is, and represents the motion of the dielectric medium. In the simplest case where the medium is at rest in laboratory’s frame, we have 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 . Our next aim is to write this equation as a quadratic form in a four-dimensional riemannian manifold. We thus define an effective spacetime with coordinates 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 in terms of , that is, . Defining the tangent four-vector as we see that (3.3) is equivalent to the expression with the effective metric 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 “” is real, and thus the potential is such that . 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 and “+” if . The four-dimensional geodesics of this geometry is given by We first note that, because the potential depends only on position by hypothesis we have We thus obtain that, because , the “0” component of the equation is identically satisfied. Also, because we have , 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 in the geodesic equation to obtain Newton’s equations). We thus arrive at the following conclusion.

The trajectories of particles with given energy 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 Trajectories are simply given by the condition . 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 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 () 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.”

4. 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.

4.1. The Geometry of Harmonic Oscillations

We start with a one-dimensional oscillator with potential energy . The () dimensional effective line element (3.10) is of the form where “” is the amplitude of the oscillator. Note that this metric is hyperbolic only if . Thus, this effective spacetime is such that there exists two geometrical boundaries given by the conditions . This is an immediate consequence of the fact that we are considering only trajectories such that the total energy is fixed, and thus the motion is not defined anywhere.

The null trajectories may be obtained by imposing the condition , that is, 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 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 is not well defined in the limit . In this limit the effective volume 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 .

Figure 1: Harmonic oscillator space-time diagram.

At a given point of the diagram there exist only two admissible trajectories with fixed energy . These two curves determine at each point 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 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.

4.2. The Geometry of Newtonian Gravitation

We now turn to the motion of point test particles with mass “” in the Newtonian gravitational potential. As in the previous section, we are interested in the motions associated to a given total energy . 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 . We will consider only the case of bounded orbits, that is, . The other situations will be considered elsewhere.

4.2.1. Effective Metrics

In this case we have a maximum admissible radii 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 . This is because outside this region, that is, for , the movement is not allowed for the considered energies .

This fact has an interesting geometrical interpretation in our scheme. First, let us note that the square root of the determinant is given by This relation implies that the element of effective volume vanishes for the values and . 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 One immediately obtains for the Ricci scalar the expression This object diverges for the values and (note that ). Thus, the effective () 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 . As we mentioned before, all possible trajectories with this energy cannot scape from the region . 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 ) where As usual, we assume that the orbit lies in the plane . Thus, (4.10) ensures that, if and 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 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 We finally turn to (4.9). As expected, this is not a simple equation. Nevertheless, we can bypass this difficulty by considering the constraint (4.8). We obtain This is the equation of motion of the radial part written in terms of the parameter . By introducing the auxiliary variable we obtain Deriving with respect to we arrive at the familiar Newtonian equation for the orbit which admits the solution where because always. Thus, we succeeded to find the orbits around the field of a sphericall object with mass “” in terms of the null geodesics of the effective metric (4.4).

4.3. 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 . The effective metric is given by the expression where 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.

4.3.1. The Case of the Electric Dipole

One interesting situation appears in the geometrization of the static electric dipole. If the dipole moment is oriented in the direction, we obtain, for large distances, For the sake of simplicity lets us consider and positive energies . The effective line element reads where . Note that this effective spacetime only makes sense if This is an immediate consequence of the fact that the energy 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 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 , 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].

5. 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 () effective spacetimes in a similar way as it appears in the context of analog gravity models. The mechanical effective metric is given explicitly in terms of the potential energy and the total energy 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.

Acknowledgments

The author would like to thank M. Novello, F. T. Falciano, S. E. P. Bergliaffa, and M. Borba for usefull comments. This work is supported by Faperj, Brazil.

References

M. Novello and E. Goulart, “Beyond analog gravity: the case of exceptional dynamics,” Classical and Quantum Gravity, vol. 28, no. 14, Article ID 145022, 2011.

G. Krein, G. Menezes, and N. F. Svaiter, “Analog model for quantum gravity effects: phonons in random fluids,” Physical Review Letters, vol. 105, Article ID 131301, 2010.

M. Visser and C. Molina-Paris, “Acoustic geometry for general relativistic barotropic irrotational fluid flow,” New Journal of Physics, vol. 12, Article ID 095014, 2010.

C. Barcelo, S. Liberati, and M. Visser, “Refringence, field theory and normal modes,” Classical and Quantum Gravity, vol. 19, no. 11, pp. 2961–2982, 2002.

M. Visser, C. Barcelo, and S. Liberati, “Bi-refringence versus bi-metricity,” http://arxiv.org/abs/gr-qc/0204017v1

T. A. Jacobson and G. E. Volovik, “Event horizons and ergoregions in 3He,” Physical Review D, vol. 58, 7 pages, 1998.

R. Schutzhold, “Emergent horizons in the laboratory,” Classical and Quantum Gravity, vol. 25, Article ID 114011, 2008.

M. Novello, “Effective geometry in nonlinear electrodynamics,” International Journal of Modern Physics A, vol. 17, no. 29, pp. 4187–4196, 2002.

M. Novello and J. M. Salin, “Nonlinear electrodynamics can generate a closed spacelike path for photons,” Physical Review D, vol. 63, Article ID 103516, 5 pages, 2001.

M. Novello and J. M. Salim, “Effective electromagnetic geometry,” Physical Review D, vol. 63, Article ID 083511, 4 pages, 2001.

M. Novello, V. A. De Lorenci, J. M. Salim, and R. Klippert, “Geometrical aspects of light propagation in nonlinear electrodynamics,” Physical Review D, vol. 61, no. 4, 10 pages, 2000.

C. Mayoral, A. Recati, A. Fabbri, R. Parentani, R. Balbinot, and I. Carusotto, “Acoustic white holes in flowing atomic Bose-Einstein condensates,” New Journal of Physics, vol. 13, Article ID 025007, 2011.

U. Leonhardt and T. G. Philbin, “Chapter 2 transformation optics and the geometry of light,” Progress in Optics, vol. 53, pp. 69–152, 2009.

S. Weinfurtner, E. W. Tedford, M. C. J. Penrice, W. G. Unruh, and G. A. Lawrence, “Measurement of stimulated hawking emission in an analogue system,” Physical Review Letters, vol. 106, 4 pages, 2011.

C. Barcelo, S. Liberati, and M. Visser, “Analogue gravity,” Living Reviews in Relativity, vol. 8, no. 12, 2005.

M. Novello, M. Visser, and G. Volovik, Eds., Artificial Black Holes, World Scientific, 2002.

W. Gordon, “Zur lichtfortp anzung nach der relativitat-stheorie,” Annalen der Physik, vol. 72, pp. 421–456, 1923.

W. G. Unruh, “Experimental black-hole evaporation?” Physical Review Letters, vol. 46, no. 21, pp. 1351–1353, 1981.

R. Dugas, A History of Mechanics, Dover Publications, 1988.

J. Lutzen, Mechanistic Images in Geometric Form, Oxford University Press, 2005.

J. L. Synge, “On the geometry of dynamics,” Transactions of the Royal Society, 1926.

J. L. Synge, “Mechanical Models of Spaces With Positive-Definite Line-Element,” Annals of Mathematics, vol. 36, no. 3, 1935.

C. Lanczos, The Variational Principles of Mechanics, Dover Publications, 4th edition, 1986.

A. Lichnerowicz, “Sur la transformation des equations de la dynamique,” Comptes rendus de l'Académie des Science de Paris, vol. 223, pp. 649–651, 1946.

J. Evans and M. Rosenquist, American Journal of Physics, vol. 54, no. 10, p. 876, 1986.

J. Evans, K. K. Nandi, and A. Islam, “The optical-mechanical analogy in general relativity: new methods for the paths of light and of the planets,” American Journal of Physics, vol. 64, no. 11, pp. 1404–1415, 1996.

J. Evans, K. K. Nandi, and A. Islam, “The optical-mechanical analogy in general relativity: exact Newtonian forms for the equations of motion of particles and photons,” General Relativity and Gravitation, vol. 28, no. 4, pp. 413–439, 1996.

J. L. Synge, Geometrical Optics: An Introduction to Hamilton's Method, Cambridge University Press, 1937.

A. Lichnerowicz, Ed., Elements de Calcul Tensoriel, Jacques Gabay, 1987.

D. A. Genov, S. Zhang, and X. Zhang, “Mimicking celestial mechanics in metamaterials,” Nature, vol. 5, no. 9, pp. 687–692, 2009.

U. Leonhardt and T. G. Philbin, “General relativity in electrical engineering,” New Journal of Physics, vol. 8, article 247, 2006.

U. Leonhardt and T. Tyc, “Broadband invisibility by non-euclidean cloaking,” Science, vol. 323, no. 5910, pp. 110–112, 2009.

V. M. Shalaev, “Physics: transforming light,” Science, vol. 322, no. 5900, pp. 384–386, 2008.

D. Schurig, J. J. Mock, B. J. Justice et al., “Metamaterial electromagnetic cloak at microwave frequencies,” Science, vol. 314, no. 5801, pp. 977–980, 2006.

R. A. Shelby, D. R. Smith, and S. Schultz, “Experimental verification of a negative index of refraction,” Science, vol. 292, no. 5514, pp. 77–79, 2001.

O. Hess, “Optics: farewell to flatland,” Nature, vol. 455, no. 7211, pp. 299–300, 2008.