Christopher Nolan's science fiction movie Interstellar offers a variety of opportunities for students in elementary courses on general relativity theory. This paper describes such opportunities, including: (i) At the motivational level, the manner in which elementary relativity concepts underlie the wormhole visualizations seen in the movie; (ii) At the briefest computational level, instructive calculations with simple but intriguing wormhole metrics, including, e.g., constructing embedding diagrams for the threeparameter wormhole that was used by our visual effects team and Christopher Nolan in scoping out possible wormhole geometries for the movie; (iii) Combining the proper reference frame of a camera with solutions of the geodesic equation, to construct a lightraytracing map backward in time from a camera's local sky to a wormhole's two celestial spheres; (iv) Implementing this map, for example, in Mathematica, Maple or Matlab, and using that implementation to construct images of what a camera sees when near or inside a wormhole; (v) With the student's implementation, exploring how the wormhole's three parameters influence what the camera sees—which is precisely how Christopher Nolan, using our implementation, chose the parameters for Interstellar's wormhole; (vi) Using the student's implementation, exploring the wormhole's Einstein ring and particularly the peculiar motions of star images near the ring, and exploring what it looks like to travel through a wormhole.
I. INTRODUCTION
A. The context and purposes of this paper
In 1988, in connection with Carl Sagan's novel Contact,^{1} later made into a movie,^{2} one of the authors published an article in this journal about wormholes as a tool for teaching general relativity.^{3}
This article is a followup, a quarter century later, in the context of Christopher Nolan's movie Interstellar^{4} and Kip Thorne's associated book The Science of Interstellar.^{5} Like Contact, Interstellar has real science built into its fabric, thanks to a strong science commitment by the director, screenwriters, producers, and visual effects team, and thanks to Thorne's role as an executive producer.
Although wormholes were central to the theme of Contact and to many movies and TV shows since then, including Star Trek and Stargate, none of these have depicted correctly a wormhole as it would be seen by a nearby human. Interstellar is the first to do so. The authors of this paper, together with Christopher Nolan who made key decisions, were responsible for that depiction.
This paper has two purposes: (i) To explain how Interstellar's wormhole images were constructed and explain the decisions made on the way to their final form and (ii) to present this explanation in a way that may be useful to students and teachers in elementary courses on general relativity.
B. The status of wormholes in the real universe
Before embarking on these explanations, we briefly describe physicists' current understanding of wormholes based on much research done since 1988. For a thorough and readable, but nontechnical review, see the recent book Time Travel and Warp Drives by Allen Everett and Thomas Roman.^{6} For reviews that are more technical, see papers by Friedman and Higuchi^{7} and by Lobo.^{8}
In brief, physicists' current understanding is this:
There is no known mechanism for making wormholes, either naturally in our universe or artificially by a highly advanced civilization, but there are speculations; for example, that wormholes in hypothetical quantum foam on the Planck scale, $ G \u210f / c 3 \u223c 10 \u2212 35 \u2009 m$, might somehow be enlarged to macroscopic size.^{6,9}
Any creation of a wormhole where initially there is none would require a change in the topology of space, which would entail, in classical, nonquantum physics, both negative energy and closed timelike curves (the possibility of backward time travel)—according to theorems by Frank Tipler and Robert Geroch.^{7} It is likely the laws of physics forbid this. Likely, but not certain.
A wormhole will pinch off so quickly that nothing can travel through it, unless it has “exotic matter” at its throat—matter (or fields) that, at least in some reference frames, has negative energy density. Although such negative energy density is permitted by the laws of physics (e.g., in the Casimir effect, the electromagnetic field between two highly conducting plates), there are quantum inequalities that limit the amount of negative energy that can be collected in a small region of space and how long it can be there; and these appear to place severe limits on the sizes of traversable wormholes (wormholes through which things can travel at the speed of light or slower).^{6} The implications of these inequalities are not yet fully clear, but it seems likely that, after some strengthening, they will prevent macroscopic wormholes like the one in Interstellar from staying open long enough for a spaceship to travel through. Likely, but not certain.
The research leading to these conclusions has been performed ignoring the possibility that our universe, with its four spacetime dimensions, resides in a higher dimensional bulk with one or more large extra dimensions, the kind of bulk envisioned in Interstellar's “fifth dimension.” Only a little is known about how such a bulk might influence the existence of traversable wormholes, but one intriguing thing is clear: Properties of the bulk can, at least in principle, hold a wormhole open without any need for exotic matter in our four dimensional universe (our “brane”).^{8} But the words “in principle” just hide our great ignorance about our universe in higher dimensions.
In view of this current understanding, it seems very unlikely to us that traversable wormholes exist naturally in our universe, and the prospects for highly advanced civilizations to make them artificially are also pretty dim.
Nevertheless, the distances from our solar system to others are so huge that there is little hope, using rocket technology, for humans to travel to other stars in the next century or two;^{10} so wormholes, quite naturally, have become a staple of science fiction.
C. The genesis of our research on wormholes
This paper is a collaboration between Caltech physicist Kip Thorne, and computer graphics artists at Double Negative Visual Effects in London. We came together in May 2013, when Christopher Nolan asked us to collaborate on building, for Interstellar, realistic images of a wormhole, and also a fast spinning black hole and its accretion disk, with ultrahigh (IMAX) resolution and smoothness. We saw this not only as an opportunity to bring realistic wormholes and black holes into the Hollywood arena but also an opportunity to create images of wormholes and black holes for relativity and astrophysics research.
Elsewhere^{12} we describe the simulation code that we wrote for this project—DNGR for “Double Negative Gravitational Renderer”—and the blackhole and accretiondisk images we generated with it, and also some new insights into gravitational lensing by black holes that it has revealed. In this paper, we focus on wormholes, which are much easier to model mathematically than Interstellar's fast spinning black hole, and are far more easily incorporated into elementary courses on general relativity.
In our modelling of Interstellar's wormhole, we pretended we were engineers in some arbitrarily advanced civilization, and that the laws of physics place no constraints on the wormhole geometries our construction crews can build. (This is almost certainly false; the quantum inequalities mentioned above, or other physical laws, likely place strong constraints on wormhole geometries, if wormholes are allowed at all—but we know so little about those constraints that we chose to ignore them.) In this spirit, we wrote down the spacetime metrics for candidate wormholes for the movie, and then proceeded to visualize them.
D. Overview of this paper
We begin in Sec. II by presenting the spacetime metrics for several wormholes and visualizing them with embedding diagrams. These metrics include, most importantly, the threeparameter “Dneg wormhole” metric used in our work on the movie Interstellar. Then we discuss adding a Newtoniantype gravitational potential to our Dneg metric, to produce the gravitational pull that Christopher Nolan wanted, and the potential's unimportance for making wormhole images.
In Sec. III, we describe how light rays, traveling backward in time from a camera to the wormhole's two celestial spheres, generate a map that can be used to produce images of the wormhole and of objects seen through or around it; and we discuss our implementations of that map to make the images seen in Interstellar. In the Appendix, we present a fairly simple computational procedure by which students can generate their own map and thence their own images.
In Sec. IV, we use our own implementation of the map to describe the influence of the Dneg wormhole's three parameters on what the camera sees. Then, in Secs. V and VI, we discuss Christopher Nolan's use of these kinds of implementations to choose the parameter values for Interstellar's wormhole; we discuss the resulting wormhole images that appear in Interstellar, including that wormhole's Einstein ring, which can be explored by watching the movie or its trailers, or in students' own implementations of the raytracing map; and we discuss images made by a camera travelling through the wormhole, that do not appear in the movie. Finally in Sec. VII, we present brief conclusions.
Scattered throughout the paper are suggestions of calculations and projects for students in elementary courses on general relativity. And throughout, as is common in relativity, we use “geometrized units” in which Newton's gravitational constant G and the speed of light c are set equal to unity, so time is measured in length units: $ 1 \u2009 s = c \xd7 1 \u2009 s = 2.998 \xd7 10 8 \u2009 m$; and mass is expressed in length units: $ 1 \u2009 kg = ( G / c 2 ) \xd7 1 \u2009 kg = 0.742 \xd7 10 \u2212 27 \u2009 m$; and the mass of the Sun is 1.476 km.
II. SPACETIME METRICS FOR WORMHOLES AND EMBEDDING DIAGRAMS
In general relativity, the curvature of spacetime can be expressed mathematically in terms of a spacetime metric. In this section, we review a simple example of this: the metric for an Ellis wormhole; and then we discuss the metric for the Double Negative (Dneg) wormhole that we designed for Interstellar.
A. The Ellis wormhole
In 1973 Homer Ellis^{13} introduced the following metric for a hypothetical wormhole, which he called a “drainhole”^{14}
where r is a function of the coordinate ℓ given by
and ρ is a constant.
As always in general relativity, one does not need to be told anything about the coordinate system in order to figure out the spacetime geometry described by the metric; the metric by itself tells us everything. Deducing everything is a good exercise for students. Here is how we do so.
First, in −dt^{2} the minus sign tells us that t, at fixed ℓ, θ, $\varphi $, increases in a timelike direction; and the absence of any factor multiplying −dt^{2} tells us that t is, in fact, proper time (physical time) measured by somebody at rest in the spatial, $ { \u2113 , \theta , \varphi }$, coordinate system.
Second, the expression $ r 2 ( d \theta 2 + sin 2 \theta \u2009 d \varphi 2 )$ is the familiar metric for the surface of a sphere with circumference 2πr and surface area 4πr^{2}, written in spherical polar coordinates $ { \theta , \varphi }$, so the Ellis wormhole must be spherically symmetric. As we would in flat space, we shall use the name “radius” for the sphere's circumference divided by 2π, i.e., for r. For the Ellis wormhole, this radius is $ r = \rho 2 + \u2113 2$.
Third, from the plus sign in front of dℓ^{2} we infer that ℓ is a spatial coordinate; and since there are no cross terms dℓdθ or $ d \u2113 \u2009 d \varphi $, the coordinate lines of constant θ and $\varphi $, with increasing ℓ, must be radial lines; and since dℓ^{2} has no multiplying coefficient, ℓ must be the proper distance (physical distance) traveled in that radial direction.
Fourth, when ℓ is large and negative, the radii of spheres $ r = \rho 2 + \u2113 2$ is large and approximately equal to $  \u2113 $. When ℓ increases to zero, r decreases to its minimum value ρ. And when ℓ increases onward to a very large (positive) value, r increases once again, becoming approximately ℓ. This tells us that the metric represents a wormhole with throat radius ρ, connecting two asymptotically flat regions of space, ℓ → −∞ and ℓ → +∞.
In Hartle's textbook,^{11} a number of illustrative calculations are carried out using Ellis's wormhole metric as an example. The most interesting is a computation, in Sec. VII, of what the twodimensional equatorial surfaces (surfaces with constant t and θ = π/2) look like when embedded in a flat 3dimensional space, the embedding space. Hartle shows that equatorial surfaces have the form shown in Fig. 1—a form familiar from popular accounts of wormholes.
Figure 1 is called an “embedding diagram” for the wormhole. We discuss embedding diagrams further in Sec. II B 3, in the context of our Dneg wormhole.
Thomas Müller and colleagues^{15} have visualized an Ellis wormhole in various environments by methods similar to those that we lay out below.
B. The double negative threeparameter wormhole
The Ellis wormhole was not an appropriate starting point for our Interstellar work. Christopher Nolan, the movie's director, wanted to see how the wormhole's visual appearance depends on its shape, so the shape had to be adjustable, which is not the case for the Ellis wormhole.
So for Interstellar we designed a wormhole with three free shaping parameters and produced images of what a camera orbiting the wormhole would see for various values of the parameters. Christopher Nolan and Paul Franklin, the leader of our Dneg effort, then discussed the images; and based on them, Nolan chose the parameter values for the movie's wormhole.
In this section, we explain our threeparameter Double Negative (Dneg) wormhole in three steps: First, a variant with just two parameters (the length and radius of the wormhole's interior) and with sharp transitions from its interior to its exteriors; then a variant with a third parameter, called the lensing length, that smooths the transitions; and finally a variant in which we add a gravitational pull.
1. Wormhole with sharp transitions
Our wormhole with sharp transitions is a simple cylinder of length 2a, whose cross sections are spheres, all with the same radius ρ; this cylinder is joined at its ends onto flat threedimensional spaces with balls of radius ρ removed. This wormhole's embedding diagram is shown in Fig. 2. As always, the embedding diagram has one spatial dimension removed, so the wormhole's cross sections appear as circles rather than spheres.
Using the same kinds of spherical polar coordinates as for the Ellis wormhole above, the spacetime metric has the general wormhole form (1) with
2. Dneg wormhole without gravity
Our second step is to smooth the transitions between the wormhole interior $  \u2113  < a$ (the cylinder) and the two external universes $  \u2113  > a$. As we shall see, the smoothed transitions give rise to gravitational lensing (distortions) of the star field behind each wormhole mouth. Such gravitational lensing is a big deal in astrophysics and cosmology these days (see, e.g., the Gravitational Lensing Resource Letter^{16}) and, as we discuss in Sec. V C, it shows up in a rather weird way in Interstellar near the edges of the wormhole image.
Somewhat arbitrarily, we chose to make the transition have approximately the same form as that from the throat (horizon) of a nonspinning black hole to the external universe in which the hole lives. Such a hole's metric (the “Schwarzschild metric”) has a form that is most simply written using radius r as the outward coordinate rather than proper distance ℓ:
where $M$ is the black hole's mass. Comparing the spatial part of this metric (t = constant) with our general wormhole metric (1), we see that $ d \u2113 = \xb1 d r / 1 \u2212 2 M / r$, which can easily be integrated to obtain the proper distance traveled as a function of radius, ℓ(r). What we want, however, is r as a function of ℓ, and we want it in an analytic form that is easy to work with. So for our Dneg wormhole, we choose a fairly simple analytic function that is roughly the same as the Schwarzschild r(ℓ).
Specifically, outside the wormhole's cylindrical interior, we chose
(Students might want to compare this graphically with the inverse of the Schwarzschild $ \u2113 = \u222b d r / 1 \u2212 2 M / r$, plotting, e.g., r − ρ for our wormhole as a function of $  \u2113  \u2212 a$; and r − 2M of Schwarzschild as a function of distance from the Schwarzschild horizon r = 2M.) Within the wormhole's cylindrical interior, we chose, of course,
Equations (5) for r(ℓ), together with our general wormhole metric (1), describe the spacetime geometry of the Dneg wormhole without gravity.
For the Schwarzschild metric, the throat radius ρ is equal to twice the black hole's mass (in geometrized units), $ \rho = 2 M$. For our Dneg wormhole, we choose the two parameters ρ and $M$ to be independent: they represent the wormhole's radius and the gentleness of the transition from the wormhole's cylindrical interior to its asymptotically flat exterior.
We shall refer to the ends of the cylindrical interior, ℓ = ±a, as the wormhole's mouths. They are spheres with circumferences 2πρ.
3. Embedding diagrams for the Dneg wormhole
We construct embedding diagrams for the Dneg wormhole (and any other spherical wormhole) by comparing the spatial metric of the wormhole's twodimensional equatorial surface $ d s 2 = d \u2113 2 + r 2 ( \u2113 ) \u2009 d \varphi 2$ with the spatial metric of the embedding space. Doing so is a good exercise for students. For the embedding space, we choose cylindrical coordinates with the symmetry axis along the wormhole's center line. Then (as in Figs. 1 and 2), the embedding space and the wormhole share the same radial coordinate r and angular coordinate $\varphi $, so with z the embeddingspace height above the wormhole's midplane, the embeddingspace metric is $ d s 2 = d z 2 + d r 2 + r 2 \u2009 d \varphi 2$. Equating this to the wormhole metric, we see that^{17} $ d z 2 + d r 2 = d \u2113 2$, which gives us an equation for the height z of the wormhole surface as a function of distance ℓ through the wormhole:
By inserting the Dneg radius function (5) into this expression and performing the integral numerically, we obtain the wormhole shapes shown in Fig. 3 (and also Figs. 7 and 9).
The actual shape of this embedding diagram depends on two dimensionless ratios of the Dneg metric's three parameters: the wormhole's lengthtodiameter ratio 2a/2ρ = a/ρ, and its ratio $ M / \rho $. For chosen values of these ratios, the wormhole's size is then fixed by its interior radius ρ, which Christopher Nolan chose to be one kilometer in Interstellar, so with the technology of the movie's era the wormhole's gravitational lensing of our galaxy's star field can be seen from Earth, but barely so.^{18}
In the embedding diagram of Fig. 3, instead of depicting $M$ we depict the lateral distance $W$ in the embedding space over which the wormhole's surface changes from vertical to 45 degrees. This $W$ is related to $M$ by^{19}
We call this $W$ the wormhole's Lensing width, and we often use it in place of $M$ as the wormhole's third parameter.
4. Dneg wormhole with gravity
Christopher Nolan asked for the movie's spacecraft Endurance to travel along a trajectory that gives enough time for the audience to view the wormhole up close before Cooper, the pilot, initiates descent into the wormhole's mouth. Our Double Negative team designed such a trajectory, which required that the wormhole have a gravitational acceleration of order the Earth's, ∼10 m/s^{2}, or less. This is so weak that it can be described accurately by a Newtonian gravitational potential Φ of magnitude $  \Phi  \u226a c 2 = 1$ (see below), that shows up in the time part of the metric. More specifically, we modify the wormhole's metric (1) to read
The sign of Φ is negative (so the wormhole's gravity will be attractive), and spherical symmetry dictates that it be a function only of ℓ.
According to the equivalence principle, the gravitational acceleration experienced by a particle at rest outside or inside the wormhole (at fixed spatial coordinates $ { \u2113 , \theta , \varphi } = constant$) is the negative of that particle's 4acceleration. Since the 4acceleration is orthogonal to the particle's 4velocity, which points in the time direction, its gravitational acceleration is purely spatial in the coordinate system $ { t , \u2113 , \theta , \varphi }$. It is a nice exercise for students to compute the particle's 4acceleration and thence its gravitational acceleration. The result, aside from negligible fractional corrections of order $  \Phi $, is
where $ e \u2113 \u0302$ is the unit vector pointing in the radial direction. Students may have seen an equation analogous to Eq. (8) when space is nearly flat, and a calculation in that case which yields Eq. (9) for g (e.g., Sec. VI of Hartle^{11}). Although for the wormhole metric (8), with r given by Eqs. (5) or (2), space is far from flat, Eq. (9) is still true—a deep fact that students would do well to absorb and generalize.
It is reasonable to choose the gravitational acceleration $ g =  g  =  d \Phi / d \u2113 $ to fall off as ∼1/(distance)^{2} as we move away from the wormhole mouth; or at least faster than ∼1/(distance). Integrating $ g =  d \Phi / d \u2113 $ radially and using this rapid falloff, the student can deduce that the magnitude of Φ is of order g times the wormhole's radius ρ. With a gravitational acceleration $ g =  g  \u2272 10 \u2009 m / s 2$ and ρ = 1 km, this gives $  \Phi  \u223c  g  \rho \u2272 10 4 \u2009 ( m / s ) 2 \u223c 10 \u2212 12$. Here, we have divided by the speed of light squared to bring this into our geometrized units.
Such a tiny gravitational potential corresponds to a slowing of time near the wormhole by the same small amount, no more than a part in 10^{12} [cf. the time part of the metric (8)]. This is so small as to be utterly unimportant in the movie, and so small that, when computing the propagation of light rays through the wormhole, to ultrahigh accuracy we can ignore Φ and use the Dneg metric without gravity. We shall do so.
III. MAPPING A WORMHOLE'S TWO CELESTIAL SPHERES ONTO A CAMERA'S SKY
A. Foundations for the Map
A camera inside or near a wormhole receives light rays from light sources and uses them to create images. In this paper we shall assume, for simplicity, that all the light sources are far from the wormhole, so far that we can idealize them as lying on “celestial spheres” at ℓ → −∞ (lower celestial sphere; Saturn side of the wormhole in the movie Interstellar) and ℓ → +∞ (upper celestial sphere; Gargantua side in Interstellar); see Fig. 4. (Gargantua is a supermassive black hole in the movie that humans visit.) Some light rays carry light from the lower celestial sphere to the camera's local sky (e.g., Ray 1 in Fig. 4); others carry light from the upper celestial sphere to the camera's local sky (e.g., Ray 2). Each of these rays is a null geodesic through the wormhole's spacetime.
On each celestial sphere we set up spherical polar coordinates $ { \theta \u2032 , \varphi \u2032}$, which are the limits of the spherical polar coordinates $ { \theta , \varphi }$ as ℓ → ±∞. We draw these two celestial spheres in Fig. 5, a diagram of the threedimensional space around each wormhole mouth, with the curvature of space not shown. Notice that we choose to draw the north polar axes θ = 0 pointing away from each other and the south polar axes θ = π pointing toward each other. This is rather arbitrary, but it feels comfortable to us when we contemplate the embedding diagram of Fig. 4.
We assume the camera moves at speeds very low compared to light speed (as it does in Interstellar), so relativistic aberration and doppler shifts are unimportant, Therefore, when computing images the camera makes, we can treat the camera as at rest in the $ { \u2113 , \theta , \varphi }$ coordinate system.
We can think of the camera as having a local sky on which there are spherical polar coordinates $ { \theta cs , \varphi cs}$ (“cs” for camera sky; not to be confused with celestial sphere!); Fig. 5. In more technical language, $ { \theta cs , \varphi cs}$ are spherical polar coordinates for the tangent space at the camera's location.
A light ray that heads backward in time from the camera (e.g., Ray 1 or 2 in Fig. 4), traveling in the $ { \theta cs , \varphi cs}$ direction, ultimately winds up at location $ { \theta \u2032 , \varphi \u2032}$ on one of the wormhole's two celestial spheres. It brings to $ { \theta cs , \varphi cs}$ on the camera's sky an image of whatever was at $ { \theta \u2032 , \varphi \u2032}$ on the celestial sphere.
This means that the key to making images of what the camera sees is a rayinduced map from the camera's sky to the celestial spheres: $ { \theta \u2032 , \varphi \u2032 , s}$ as a function of $ { \theta cs , \varphi cs}$, where the parameter s tells us which celestial sphere the backward light ray reaches, the upper one (s = +) or the lower one (s = −).
In the Appendix, we sketch a rather simple computational procedure by which students can compute this map and then, using it, can construct images of wormholes and their surroundings. (We also describe a Mathematica implementation of this procedure by this paper's computationally challenged author Kip Thorne.)
B. Our DNGR mapping and image making
To produce the IMAX images needed for Interstellar, at Double Negative we developed a much more sophisticated implementation of the map within a computer code that we call DNGR^{12} (Double Negative Gravitational Renderer). In DNGR, we use ray bundles (light beams) to do the mapping rather than just light rays. We begin with a circular light beam about one pixel in size at the camera and trace it backward in time to its origin on a celestial sphere using the ray equations (A7), plus the general relativistic equation of geodesic deviation, which evolves the beam's size and shape. At the celestial sphere, the beam is an ellipse, often highly eccentric. We integrate up the image data within that ellipse to deduce the light traveling into the camera's circular pixel. We also do spatial filtering to smooth artifacts and time filtering to mimic the behavior of a movie camera (when the image is changing rapidly), and we sometimes add lens flare to mimic the effects of light scattering and diffraction in a movie camera's lens.
Elsewhere^{12} we give some details of these various “bells and whistles” for a camera orbiting a black hole rather than a wormhole. They are essentially the same for a wormhole.
However, fairly nice images can be produced without any of these bells and whistles, using the simple procedure described in the Appendix and thus are within easy reach of students in an elementary course on general relativity.
IV. THE INFLUENCE OF THE WORMHOLE'S PARAMETERS ON WHAT THE CAMERA SEES
For Christopher Nolan's perusal in choosing Interstellar's wormhole parameters, we used our map to make images of the galaxy in which the black hole Gargantua resides, as viewed from the Saturn side of the wormhole; see below. But for this paper, and the book^{5} that Thorne has written about the science of Interstellar, we find it more instructive, pedagogically, to show images of Saturn and its rings as seen through the wormhole from the Gargantua side. This section is a more quantitative version of a discussion of this in Chapter 15 of that book.^{5}
Figure 6 shows the simple Saturn image that we placed on the lower celestial sphere of Fig. 5, and a star field that we placed on the upper celestial sphere (the Gargantua side of the wormhole). Both images are mapped from the celestial sphere onto a flat rectangle with azimuthal angle $\varphi $ running horizontally and polar angle θ vertically. In computer graphics, this type of image is known as a longitudelatitude map.^{21}
A. Influence of the Wormhole's length
In Fig. 7, we explore the influence of the wormhole's length on the camerasky image produced by these two celestial spheres. Specifically, we hold the wormhole's lensing width fixed at a fairly small value, $ W = 0.05 \rho $, and we vary the wormhole's length from 2a = 0.01ρ (top picture), to 2a = ρ (middle picture), to 2a = 10ρ (bottom picture).
Because Saturn and its rings are white and the sky around it is black, while the star field on the Gargantua side of the wormhole is blue, we can easily identify the edge of the wormhole mouth as the transition from blackandwhite to blue. (The light's colors are preserved as the light travels near and through the wormhole because we have assumed the wormhole's gravity is weak, $  \Phi  \u226a 1$; there are no significant gravitational frequency shifts.)
Through a short wormhole (top), the camera sees a large distorted image of Saturn nearly filling the right half of the wormhole mouth. This is the primary image, carried by light rays that travel on the shortest possible paths through the wormhole from Saturn to camera, such as path 1 in Fig. 8. There is also a very thin, lenticular, secondary image of Saturn, barely discernable, near the left edge of the wormhole mouth. It is brought to the camera by light rays that travel around the left side of the wormhole (e.g., path 2 in Fig. 8)—a longer route than for the primary image. The lenticular structure at the lower right is blue, so it is a secondary gravitationally lensed image of the blue star field that resides on the camera's side of the wormhole.
As the wormhole is lengthened (middle of Fig. 7), the primary and secondary images move inward and shrink in size. A lenticular tertiary image emerges from the mouth's right edge, carried by rays like 3 in Fig. 8 that wrap around the wormhole once; and a fourth (faint) lenticular image emerges from the left side, carried by rays like 4 that wrap around the wormhole in the opposite direction, one and a half times.
As the wormhole is lengthened more and more (bottom of Fig. 7), the existing images shrink and move inward toward the mouth's center, and new images emerge, one after another, from the right then left then right… sides of the mouth.
For a short wormhole, all these images were already present, very near the wormhole's edge; but they were so thin as to be unresolvable. Lengthening the wormhole moved them inward and made them thick enough to see.
B. Influence of the Wormhole's lensing width
In Fig. 9, we explore the influence of the wormhole's lensing width on what the camera sees. We hold its length fixed and fairly small: equal to its radius, 2a = ρ.
For small lensing width $ W = 0.014 \rho $ (top), the transition from the wormhole's cylindrical interior to its asymptotically flat exterior is quite sharp; so, not surprisingly, the camera sees an exterior, blue star field that extends with little distortion right up to the edge of the wormhole mouth.
By contrast, when the lensing width is larger, $ W = 0.43 \rho $ (bottom), the external star field is greatly distorted by gravitational lensing. The dark cloud on the upper left side of the wormhole is enlarged and pushed out of the cropped picture, and we see a big secondary image of the cloud on the wormhole's lower right and a tertiary image on its upper left. We also see lensing of the wormhole mouth itself: it is enlarged; and lensing of the image that comes through the wormhole from the Saturn side. The lenticular secondary image of Saturn near the mouth's left edge is thickened, while the primary image is shrunken a bit and moved inward to make room for a new tertiary image on the right.
V. INTERSTELLAR'S WORMHOLE
After reviewing images analogous to Figs. 7 and 9, but with Saturn replaced by the stars and nebulae of Interstellar's distant galaxy (the galaxy on the Gargantua side of the wormhole), Christopher Nolan made his choice for the parameters of Interstellar's wormhole.
He chose a very short wormhole: length 2a = 0.01ρ as in the top panel of Fig. 7; for greater lengths, the multiple images would be confusing to a mass audience. And he chose a modest lensing width: $ W = 0.05 \rho $ also as in the top panel of Fig. 7 and in between the two lensing widths of Fig. 9. This gives enough gravitational lensing to be interesting (see below), but far less lensing than for a black hole, thereby enhancing the visual distinction between Interstellar's wormhole and its black hole Gargantua.
A. Interstellar's distant galaxy
For Interstellar, a team under the leadership of authors Paul Franklin and Eugénie von Tunzelmann constructed images of the distant galaxy through a multistep process.
The distant end of the wormhole was imagined to be in the distant galaxy and closer to its center than we are to the center of our Milky Way. Consequently, the view of the surrounding galaxy must be recognisably different from the view we have from Earth: larger and brighter nebulae, more dense dust, with brighter and more numerous visible stars. This view was created as an artistic task.
Nebulae were painted (by texture artist Zoe Lord), using a combination of space photography and imagination, covering a range of color palettes. These were combined with layers of painted bright space dust and dark, silhouetted dust channels, to create a view of the galaxy with as much visual depth and complexity as possible.
Star layout was achieved by taking real star data as seen from Earth and performing various actions to make the view different: the brightest stars were removed from the data set (to avoid recognisable constellations) and the brightnesses of all the other stars were increased and shuffled. The result was a believably naturallooking star layout which was unrecognisable compared to our familiar view of the night sky from Earth.
Figure 10 is one of our distantgalaxy images, showing nebulae, space dust, and stars.
B. View through Interstellar's Wormhole
When we place this distantgalaxy image on the upper celestial sphere of Fig. 5 and place a simple star field on the lower celestial sphere, within which the camera resides, then the moving camera sees the wormhole images shown in Interstellar and its trailers; for example, Fig. 11.
Students can create similar images, using their implementation of the map described in the Appendix, and putting Fig. 10 on the upper celestial sphere. They could be invited to explore how their images change as the camera moves farther from the wormhole, closer, and through it, and as the wormhole parameters are changed.
C. The Einstein ring
Students could be encouraged to examine closely the changing image of the wormhole in Interstellar or one of its trailers, on a computer screen where the student can move the image back and forth in slow motion. Just outside the wormhole's edge, at the location marked by a dotted circle in Fig. 11, the star motions (induced by camera movement) are quite peculiar. On one side of the dotted circle, stars move rightward; on the other, leftward. The closer a star is to the circle, the faster it moves; see Fig. 12.
The circle is called the wormhole's Einstein ring. This ring is actually the ring image, on the camera's local sky, of a tiny light source that is precisely behind the wormhole and on the same end of the wormhole as the camera. That location, on the celestial sphere and precisely opposite the camera, is actually a caustic (a singular, focal point) of the camera's past light cone. As the camera orbits the wormhole, causing this caustic to sweep very close to a star, the camera sees two images of the star, one just inside the Einstein ring and the other just outside it, move rapidly around the ring in opposite directions. This is the same behavior as occurs with the Einstein ring of a black hole (see, e.g., Fig. 2 of our paper on blackhole lensing^{12}) and any other spherical gravitational lens, and it is also responsible for long, lenticular images of distant galaxies gravitationally lensed by a more nearby galaxy.^{22}
Students, having explored the wormhole's Einstein ring in a DVD or trailer of the movie, could be encouraged to go learn about Einstein rings and/or figure out for themselves how these peculiar star motions are produced. They could then use their own implementation of our map to explore whether their explanation is correct.
VI. TRIP THROUGH THE WORMHOLE
Students who have implemented the map (described in the Appendix) from the camera's local sky to the celestial spheres could be encouraged to explore, with their implementation, what it looks like to travel through the Dneg wormhole for various parameter values.
We ourselves did so, together with Christopher Nolan, as a foundation for Interstellar's wormhole trip. Because the wormhole Nolan chose to visualize from the outside (upper left of Fig. 7; images in Figs. 10 and 12) is so short and its lensing width so modest, the trip was quick and not terribly interesting, visually—not at all what Nolan wanted for his movie. So we generated additional throughthewormhole clips for him, with the wormhole parameters changed. For a long wormhole, the trip was like traveling through a long tunnel, too much like things seen in previous movies. None of the clips, for any choice of parameters, had the compelling freshness that Nolan sought.
Moreover, none had the right feel. Figure 13 illustrates this problem. This figure shows stills from a trip through a moderately short wormhole with a/ρ = 0.5 (stills that students could replicate with their implementation). Although these images are interesting, the resulting animated sequence is hard for an audience to interpret. The view of the wormhole appears to scale up from its center, growing in size until it fills the frame, and until none of the starting galaxy is visible; at this point only the new galaxy can be seen, because we now are actually inside that new galaxy. This is hard to interpret visually. Because there is no parallax or other relative motion in the frame, to the audience it looks like the camera is zooming into the center of the wormhole using the camera's zoom lens. In the visual grammar of filmmaking, this tells the audience that we are zooming in for a closer look but we are still a distance from the wormhole; in reality we are travelling through it, but this is not how it feels.
It was important for the audience to understand that the wormhole allows the Endurance (the movie's space ship) to take a shortcut through the higher dimensional bulk. To foster that understanding, Nolan asked the visual effects team to convey a sense of travel through an exotic environment, one that was thematically linked to the exterior appearance of the wormhole but also incorporated elements of passing landscapes and the sense of a rapidly approaching destination. The visual effects artists at Double Negative combined existing DNGR visualisations of the wormhole's interior with layers of interpretive effects animation derived from aerial photography of dramatic landscapes, adding lensbased photographic effects to tie everything in with the rest of the sequence. The end result was a sequence of shots that told a story comprehensible by a general audience while resembling the wormhole's interior, as simulated with DNGR.
VII. CONCLUSION
As we wrote this paper, we became more and more enthusiastic about the educational opportunities provided by our Interstellar experience. The tools we used in building, scoping out, and exploring Interstellar's wormhole—at least those discussed in this paper—should be easily accessible to fourth year undergraduates studying relativity, as well as to graduate students. And the movie itself, and our own route to the final wormhole images in the movie, may be a strong motivator for students.
ACKNOWLEDGMENTS
For extensive advice on our wormhole visualizations, the authors thank Christopher Nolan. For contributions to DNGR and its wormhole applications, the authors thank members of the Double Negative R&D team Sylvan Dieckmann, Simon Pabst, Shane Christopher, PaulGeorge Roberts, and Damien Maupu; and also Double Negative artists Zoe Lord, Fabio Zangla, Iacopo di Luigi, Finella Fan, Tristan Myles, Stephen Tew, and Peter Howlett. The construction of our code DNGR was funded by Warner Bros. Entertainment Inc., for generating visual effects for the movie Interstellar. The authors thank Warner Bros. for authorizing this code's additional use for scientific research and physics education, and in particular the work reported in this paper.
APPENDIX: THE RAYINDUCED MAP FROM THE CAMERA'S LOCAL SKY TO THE TWO CELESTIAL SPHERES
In this Appendix, we describe our fairly simple procedure for generating the map from points $ { \theta cs , \varphi cs}$ on the camera's local sky to points $ { \theta \u2032 , \varphi \u2032 , s}$ on the wormhole's celestial sphere, with s = + for the upper celestial sphere and s = − for the lower.
1. The ray equations
As we discussed in Sec. III A, the map is generated by light rays that travel backward in time from the camera to the celestial spheres. In the language of general relativity, these light rays are null (lightlike) geodesics and so are solutions of the geodesic equation
Here, the $ \Gamma \alpha \u2009 \mu \nu $ are Christoffel symbols (also called connection coefficients) that are constructable from first derivatives of the metric coefficients, and ζ is the socalled affine parameter, which varies along the geodesic.
This form of the geodesic equation is fine for analytical work, but for numerical work, it is best rewritten in the language of Hamiltonian mechanics. Elsewhere^{23} one of us will discuss, pedagogically, the advantages and the underpinnings of this Hamiltonian rewrite.
There are several different Hamiltonian formulations of the geodesic equation. The one we advocate is sometimes called the “superHamiltonian” because of its beauty and power, but we will stick to the usual word “Hamiltonian.” The general formula for this Hamiltonian is^{23,24}
Here, g^{μν} are the contravariant components of the metric, x^{α} is the coordinate of a photon traveling along the ray, and p_{α} is the generalized momentum that is canonically conjugate to x^{α}, which turns out to be the same as the covariant component of the photon's 4momentum. Hamilton's equations, with the affine parameter ζ playing the role of time, take the standard form
In the first of Eqs. (A3), the metric raises the index on the covariant momentum, so it becomes $ p \alpha = d x \alpha / d \zeta $, an expression that may be familiar to students. The second expression may not be so familiar, but it can be given as an exercise for students to show that the second equation, together with $ p \alpha = d x \alpha / d \zeta $, is equivalent to the usual form (A1) of the geodesic equation.
For the general wormhole metric (1), the (super) Hamiltonian (A2) has the simple form
Because this (super) Hamiltonian is independent of the time coordinate t and of the azimuthal coordinate $\varphi $, p_{t} and $ p \varphi $ are conserved along a ray [cf. Eq. (A3b)]. Since $ p t = d t / d \zeta = \u2212 p t$, changing the numerical value of p_{t} merely renormalizes the affine parameter ζ; so without loss of generality, we set p_{t} = −1, which implies that ζ is equal to time t [Eq. (A6)]. Since photons travel at the speed of light, ζ is also distance travelled (in our geometrized units where the speed of light is one).
We use the notation b for the conserved quantity $ p \varphi $:
Students should easily be able to show that, because we set p_{t} = −1, this b is the ray's impact parameter relative to the (arbitrarily chosen^{25}) polar axis. Because the wormhole is spherical, there is a third conserved quantity for the rays, its total angular momentum, which (with p_{t} = −1) is the same as its impact parameter B relative to the hole's center
By evaluating Hamilton's equations for the wormhole Hamiltonian (A4) and inserting the conserved quantities on the righthand side, we obtain the following ray equations:
which reaffirms that ζ = t (up to an additive constant); and, replacing ζ by t:
These equations may seem like an overly complicated way to describe a ray. Complicated, maybe, but near ideal for simple numerical integrations. They are stable and in all respects well behaved everywhere except the poles θ = 0 and θ = π, and they are easily implemented in studentfriendly software such as Mathematica, Maple, and Matlab.
2. Procedure for generating the map
It is an instructive exercise for students to verify the following procedure for constructing the map from the camera's local sky to the two celestial spheres:

Choose a camera location $ ( \u2113 c , \theta c , \varphi c )$. It might best be on the equatorial plane, θ_{c} = π/2, so the coordinate singularities at θ = 0 and θ = π are as far from the camera as possible.
 Set up local Cartesian coordinates centered on the camera, with x along the direction of increasing ℓ (toward the wormhole on the Saturn side; away from the wormhole on the Gargantua side), y along the direction of increasing $\varphi $, and z along the direction of decreasing θ,Here $ e \u2113 \u0302 , \u2009 e \theta \u0302$, and $ e \varphi \u0302$ are unit vectors that point in the ℓ, θ, and $\varphi $ directions (the hats tell us their lengths are one). Figure 5 shows these camera basis vectors for the special case where the camera is in the equatorial plane. The minus sign in our choice $ e z = \u2212 e \theta \u0302$ makes the camera's e_{z} parallel to the wormhole's polar axis on the Gargantua side of the wormhole, where ℓ is positive.$ e x = e \u2113 \u0302 , \u2003 e y = e \varphi \u0302 , \u2003 e z = \u2212 e \theta \u0302 .$

Set up a local spherical polar coordinate system for the camera's local sky in the usual way, based on the camera's local Cartesian coordinates; cf. Eq. (A9a).
 Choose a direction $ ( \theta c s , \u2009 \varphi c s )$ on the camera's local sky. The unit vector N pointing in that direction has Cartesian componentsBecause of the relationship (A8) between bases, the direction n of propagation of the incoming ray that arrives from direction −N has components in the global spherical polar basis$ N x = sin \u2009 \theta c s \u2009 cos \u2009 \varphi c s , \u2003 N y = sin \u2009 \theta c s \u2009 sin \u2009 \varphi c s , N z = cos \u2009 \theta c s .$$ n \u2113 \u0302 = \u2212 N x , \u2003 n \varphi \u0302 = \u2212 N y , \u2003 n \theta \u0302 = + N z .$
 Compute the incoming light ray's canonical momenta from(it's a nice exercise for students to deduce these equations from the relationship between the covariant components of the photon 4momentum and the components on the unit basis vectors). Then compute the ray's constants of motion from$ p \u2113 = n \u2113 \u0302 , \u2003 p \theta = r n \theta \u0302 , \u2003 p \varphi = r \u2009 sin \u2009 \theta n \varphi \u0302$$ b = p \varphi = r \u2009 sin \u2009 \theta n \varphi \u0302 , B 2 = p \theta 2 + p \varphi 2 \u2009 sin 2 \theta = r 2 ( n \theta \u0302 2 + n \varphi \u0302 2 ) .$

Take as initial conditions for ray integration that at t = 0 the ray begins at the camera's location $ ( \u2113 , \theta , \varphi ) = ( \u2113 c , \theta c , \varphi c )$ with canonical momenta (A9c) and constants of motion (A9d). Numerically integrate the ray equations (A7), subject to these initial conditions, from t = 0 backward along the ray to time t_{i} = −∞ (or some extremely negative, finite initial time t_{i}). If ℓ(t_{i}) is negative, then the ray comes from location $ { \theta \u2032 , \varphi \u2032} = { \theta ( t i ) , \varphi ( t i )}$ on the Saturn side of the wormhole, s = −. If ℓ(t_{i}) is positive, then the ray comes from location $ { \theta \u2032 , \varphi \u2032} = { \theta ( t i ) , \varphi ( t i )}$ on the Gargantua side of the wormhole, s = +.
3. Implementing the map
Evaluating this map numerically should be a moderately easy task for students. Kip Thorne, the author among us who is least adept at numerical work, did it using Mathematica, and then used that map—a numerical table of $ { \theta \u2032 , \varphi \u2032 , s}$ as a function of $ { \theta cs , \varphi cs}$—to make camerasky images of whatever was placed on the two celestial spheres. For image processing, Thorne first built an interpolation of the map using the Mathematica command ListInterpolation; and he then used this interpolated map, together with Mathematica's command ImageTransformation, to produce the camerasky image from the images on the two celestial spheres.
References
See, e.g., chapter 13 of The Science of Interstellar^{4}.
The polar axis is arbitrary because the wormhole's geometry is spherically symmetric.