We present an analysis of the classic Alcubierre metric based on conformal gravity, rather than standard general relativity. The main characteristics of the resulting warp drive remain the same as in the original study by Alcubierre, that is, effective superluminal motion is a viable outcome of the metric. We show that for particular choices of the shaping function, the Alcubierre metric in the context of conformal gravity does not violate the weak energy condition, as was the case of the original solution. In particular, the resulting warp drive does not require the use of exotic matter. Therefore, if conformal gravity is a correct extension of general relativity, superluminal motion via an Alcubierre metric might be a realistic solution, thus allowing faster-than-light interstellar travel.

1. Introduction

In 1994, Alcubierre introduced the so-called Warp Drive Metric (WDM), within the framework of general relativity (GR), which allows in principle for super-luminal motion, that is, faster-than-light travel [1]. This superluminal propulsion is achieved, respectively, by expanding and contracting the space time behind and in front of a spaceship, while the spacecraft is left inside a locally flat region of space time, within the so-called warp bubble.

In this way, the spaceship can travel at arbitrarily high speeds, without violating the laws of special and general relativity, or other known physical laws. Furthermore, the spacecraft and its occupants would also be at rest in flat space time, thus immune from high accelerations and unaffected by special relativistic effects, such as time dilation. Enormous tidal forces would only be present near the edge of the warp bubble, which can be made large enough to accommodate the volume occupied by the ship.

However, Alcubierre [1] was also the first to point out that this hypothetical solution of Einstein’s equations of GR would violate all three standard energy conditions (weak, dominant, and strong; see [24] for definitions). In particular, the violation of the weak energy condition (WEC) implies that negative energy density is required to establish the Alcubierre WDM, thus making it practically impossible to achieve this type of super-luminal motion, unless large quantities of exotic matter (i.e., with negative energy density) can be created. Since our current knowledge of this type of exotic matter is limited to some special effects in quantum field theory (such as the Casimir effect), it is unlikely that the Alcubierre WDM can be practically established within the framework of general relativity.

Following Alcubierre’s seminal paper, many other studies appeared in the literature, either proposing alternatives to the original warp drive mechanism ([5, 6]) or refining and analyzing in more detail the original idea ([719]). However, all these studies were conducted using standard GR and could not avoid the violation of the WEC, meaning that some exotic matter would always be required for faster-than-light travel. Similar issues also exist in other well-known GR solutions for super-luminal motion, such as space-time wormholes [20].

Einstein’s general relativity and the related “Standard Model” of cosmology have been highly successful in describing our universe, from the solar system up to the largest cosmological scales, but recently these theories have also led to a profound crisis in our understanding of its ultimate composition. From the original discovery of the expansion of the universe, which resulted in standard big bang cosmology, scientists have progressed a long way towards our current picture, in which the contents of the universe are today described in terms of two main components, dark matter (DM) and dark energy (DE), accounting for most of the observed universe, with ordinary matter just playing a minor role.

Since there is no evidence available yet as to the real nature of dark matter and dark energy, alternative gravitational and cosmological theories are being developed, in addition to standard explanations of dark matter/dark energy invoking the existence of exotic new particles also yet to be discovered. In line with these possible new theoretical ideas, conformal gravity (CG) has emerged as a nonstandard extension of Einstein’s GR, based on a possible symmetry of the universe: the conformal symmetry, that is, the invariability of the space-time fabric under local “stretching” of the metric (for reviews see [21, 22]). This alternative theory has been reintroduced in recent years (following the original work by Weyl [2325]), leading to cosmological models which do not require the existence of DM and DE ([2632]).

At the quantum level conformal gravity, as well as other theories with higher derivatives, was thought to be affected by the presence of “ghosts,” leading to possible instabilities of the quantum version of the theory. However, recent studies [26] have shown that CG as a quantum theory is both renormalizable and unitary, thus providing a solution to the ghost problem.

In view of a possible extension of Einstein’s general relativity into conformal gravity, in this paper we have reconsidered the Alcubierre WDM, basing it on CG rather than standard GR. In Section 2, we review the fundamental principles of CG and the calculation of the stress-energy tensor in this gravitational theory. In Section 3, we consider the Alcubierre metric in CG and compute the energy density for different shaping functions of the metric.

In particular, we will show that for certain shaping functions, the Alcubierre metric in the context of conformal gravity does not violate the weak energy condition, as was the case of the original solution. This analysis continues in Section 4, where we study other energy conditions and estimate the total energy required for this CG warp drive. Finally, in Section 5, we conclude that if CG is a correct extension of GR, super-luminal motion via an Alcubierre metric might be a realistic possibility, thus enabling faster-than-light interstellar travel without requiring exotic matter.

2. Conformal Gravity and the Stress-Energy Tensor

Weyl in 1918 ([2325]) developed the “conformal” generalization of Einstein’s relativity by introducing the conformal (or Weyl) tensor, a special combination of the Riemann tensor , the Ricci tensor , and the curvature (or Ricci) scalar [33] as follows: where is invariant under the local transformation of the following metric: The factor determines the amount of local “stretching” of the geometry, hence the name “conformal” for a theory invariant under all local stretchings of the space-time (see [22] and references therein for more details).

This conformally invariant generalization of GR was found to be a fourth-order theory, as opposed to the standard second-order general relativity, since the field equations originating from a conformally invariant Lagrangian contain derivatives up to the fourth order of the metric, with respect to the space-time coordinates. Following the works done by Bach [34], Lanczos [35], and others, CG was ultimately based on the following Weyl or conformal action (In this paper we adopt a metric signature and we follow the sign conventions of Weinberg [33]. In this section we will leave fundamental constants, such as and , in all equations, but later we will use geometrized units (, ), or c.g.s. units when needed.): or on the following equivalent expression, differing from the previous one only by a topological invariant: where and is the gravitational coupling constant of conformal gravity (see [21, 3638].) (In these cited papers, is considered a dimensionless constant by using natural units. Working with c.g.s. units, we can assign dimensions of an action to the constant so that the dimensionality of (5) will be correct.) Under the conformal transformation in (2), the Weyl tensor transforms as , while the conformal action is locally conformally invariant, the only general coordinate scalar action with such properties.

Bach [34] introduced the gravitational field equations in the presence of a stress-energy tensor (We follow here the convention [21] of introducing the stress-energy tensor so that the quantity has the dimensions of an energy density.) as as opposed to Einstein’s standard equations as where the “Bach tensor” [34] is the equivalent in CG of the Einstein curvature tensor on the left-hand side of (6).

has a very complex structure and can be defined in a compact way as [39] or in an expanded form as ([37, 40]) involving derivatives up to the fourth order of the metric with respect to space-time coordinates.

Therefore, in conformal gravity, the stress-energy tensor is computed by combining together (5) and (8) as This form of the tensor will be used in the following sections, in connection to the Alcubierre metric, to compute the energy density and other relevant quantities.

For this purpose, we have developed a special Mathematica program which enables us to compute all the tensor quantities of both GR and CG, for any given metric. In particular, this program can compute the conformal tensor in (1), the Bach tensor in (8), or the stress-energy tensor in (9) by performing all the necessary covariant derivatives. Given the complexity of these types of computations, we have tested extensively our program against the results for and computed by Kazanas and Mannheim in [38] for different metrics, obtaining a perfect agreement.

3. Alcubierre Metric, Shaping Functions, and the Weak Energy Condition

The original Alcubierre metric [1] considered a spaceship traveling along the -axis, with motion described by a function and spaceship velocity . Using the 3+1 formalism of GR, the metric was written in Cartesian coordinates as () where is the distance from the spaceship position as and is a “form function” or “shaping function” which needs to have values and , respectively, inside and outside the warp bubble, while it can have an arbitrary shape in the transition region of the warp bubble itself.

The original shaping function used by Alcubierre was where basically indicates the radius of the spherical warp bubble, while relates to the bubble thickness, which decreases with increasing values of . In the following, we will refer to the function in (12) as the “Alcubierre shaping function” (ASF).

We will show that the particular form of the shaping function can play an important role in the energy conditions for the WDM. In our analysis we tested several different functions obeying the general requirements for outlined above. In addition to the Alcubierre function above, in this paper we will also use the following: where is a positive integer. Since this particular function for is used by Hartle to illustrate the warp drive in his textbook [20], we will refer to the function in (13) as the “Hartle shaping function” (HSF).

The top panels in Figure 1 illustrate the differences between the Alcubierre shaping function (top left panel, for as used originally by Alcubierre) and the Hartle shaping function (top right panel, for ). All functions in this figure are computed for a fixed value of and at time , when the spaceship is located at the origin. All quantities shown in the different panels are plotted as a function of the coordinate of the spaceship motion and of the transverse cylindrical coordinate . Similar coordinates will also be used in the other figures. (The cylindrical coordinate should be considered as non-negative and all quantities in the figures plotted only for . However, for illustrative purposes and also to follow similar figures in the literature (such as those in [1, 11, 12]), we decided to let run on negative values in all figures, except in the last one, where we restrict for a correct energy calculation.)

The expansion/contraction function of the volume elements behind/in front of the spaceship was also computed by Alcubierre as [1] and is illustrated for (in geometrized units, i.e., in traditional units) in Figure 1(b), for the two different shaping functions. Again, the choice of the parameters , , and is the same as in the top panels in the figure. The expansion for the ASF is the same as Figure  1 in [1], while the corresponding for the HSF is slightly different but still shows expansion of the normal volume elements behind the spaceship and contraction in front of it.

The weak energy condition ([24]) requires that for all timelike vectors . Alcubierre has also shown that for the Eulerian observers in the warp drive metric, and for their 4-velocity , the following relation holds [1]: which implies that the energy density is negative everywhere for any choice of the shaping function and, therefore, the WEC is violated (also the dominant energy condition, DEC, and the strong energy condition, SEC, are violated in the analysis based on GR [1]).

This violation of the WEC in GR is illustrated in Figure 1(c), where is calculated using (15) for both shaping functions. Although the results in the two panels are slightly different, they obviously show negative energy densities and therefore a complete violation of the WEC.

The situation is different if we compute the energy density in the framework of CG, following (9), setting for simplicity, and using the completely contravariant form of the stress-energy tensor, instead of the covariant one. As seen in the bottom row of Figure 1, the energy density in CG is completely different from the one calculated within the framework of GR. In the bottom left panel is computed with the ASF and the resulting function is in part positive and in part negative, thus still violating the WEC.

However, the bottom right panel shows computed with the HSF and in this case the energy density is completely non-negative, showing that the WEC is verified and no exotic matter is needed to establish the warp drive. (In this panel of Figure 1 (and also in the other figures) we show a step discontinuity at the edge of the warp bubble, common to all our CG solutions. These discontinuities can be mathematically replaced by appropriate delta functions, as described by the equations used in the Appendix to model the Hartle shaping function. However, these delta functions do not need to be included in the physical solutions for in CG shown in the figures. Therefore, these energy densities are completely non-negative and they fully determine the expansion of the volume elements required for the warp drive effect.) This non-negative energy density plot in the bottom right panel of Figure 1 is the main result of our paper as it shows that—if CG is the correct extension of GR—it might be possible to establish a warp drive without having to use negative energy (mass), thus overcoming the main difficulty of the warp drive mechanism.

The explicit expression of in CG, computed with our Mathematica program, is rather cumbersome and is reproduced in (A.1) of the Appendix. Here we present just the graphical computation of in Figure 1, or in the other figures in this paper. We have tested the validity of these results by running the program in several different ways, including computing the stress-energy tensor for a simplified three-dimensional Alcubierre metric (coordinates , , and only), always obtaining consistent results. The differences in the two plots at the bottom of the figure can be attributed to the different shaping functions and their derivatives up to the fourth order. All these derivatives enter the complex expression of in CG, as in the master equation (9), and their interplay ultimately determines the shape of , or of the other components, in a way which is hard to predict before the actual computation is performed.

In Figure 2, we analyze the dependence of the energy density on the spaceship velocity . In this case we consider only the Hartle shaping function for and , and we compute in CG () for speeds ranging from the subluminal to the super-luminal (for we obtain the same function as in the bottom right panel of Figure 1). The shape of the energy density function is about the same for speeds up to , although the function values increase with speed. For higher velocities, the function develops two “downward lobes” which eventually become negative for speeds . This implies that the WEC is verified for speeds up to , while at higher velocities exotic matter would be required to sustain the warp drive.

This apparent “speed limit” at about might be raised or overcome completely by adopting a different shaping function, instead of the HSF used here, but this analysis would go beyond the scope of this work. We do not ascribe any physical significance to this new speed limit, which is probably just an artificial product of the shaping functions used in this analysis.

In any case, the results reported in Figure 2 show that a warp drive in CG with positive energy density is possible for a wide range of spaceship velocities; therefore, if CG is the correct extension of GR, the Alcubierre warp drive might be a viable mechanism for super-luminal travel.

In Figure 3, we present the other components of the stress-energy tensor. These were computed with the same Mathematica program, following (9) with , leading to even more complex expressions than the one for (we will omit to report these expressions for brevity). To simplify the computation, we used the covariant components and adopted cylindrical coordinates around the -axis, , instead of the Cartesian coordinates of the original Alcubierre metric. The results in the different panels are labeled accordingly. We recall that the stress-energy tensor is symmetric, , and only the nonzero components are illustrated in this figure, under similar conditions used before (CG with HSF and parameters , , ).

The shapes of the other components of are more complex than the one of , but they can all be determined analytically, either in covariant or contravariant form, using our program. If this exact form of the stress-energy tensor could be established in the region surrounding the spacecraft, warp drive motion would be feasible within the framework of conformal gravity.

We also want to point out that we set the spaceship motion in the positive direction of the -axis (setting ), and this is reflected in the symmetry, or lack thereof, of the components of . While some components, such as (or ), , , , , appear to be symmetric under the exchange , the other components, and , are not symmetric under this exchange and, therefore, these components must contain information about the spaceship direction of motion. This argument overcomes the objection addressed in [12] that since is symmetric about the plane, there is uncertainty in where the space-time is expanded/contracted, thus making it impossible for the spaceship to know in which direction of the -axis, positive or negative, to move. Rather, Figure 3 shows that a “bias” towards one of the two possible directions is induced by some components of .

Figure 4 illustrates one last dependence of the energy density on the parameters used. In this case, we set , and consider the Hartle shaping function as in (13), while varying the integer parameter . In addition to our standard value, , we have also tried values from to , as shown in the figure.

In general, increasing the value increases the internal volume of the warp bubble, where space time is flat, and therefore decreases the thickness of the bubble wall where the space-time distortion takes place, in a way similar to that of the parameter in the original Alcubierre warp drive. However, increasing also increases the energy required to establish the warp drive (as we checked by integrating the functions in Figure 4). Therefore, it appears to be convenient to use a low value for this parameter. As shown in the different panels of the figure, the first value, , does not work since it does not leave a flat space-time volume inside the bubble, while seems to create a very small volume inside the bubble. Therefore, intermediate values such as = 4–6 would appear to be more adequate to establish the warp drive. Meanwhile, the solution for would require much more energy and would not give any advantage, except reducing the thickness of the warp bubble. We also checked that changing the value of does not have a strong effect on the “speed limit” of , reported above for the case . Thus, this value of the parameter seems to be the most adequate for this type of solutions.

4. Other Energy Conditions and Warp Drive Energy Estimate

In the previous section, we have discussed at length the weak energy condition (WEC) for the conformal gravity Alcubierre warp drive. We have seen that, if the Hartle shaping function is used, this condition is not violated for a wide range of spaceship velocities, including super-luminal speeds. In this section, we will briefly analyze the other main energy conditions and estimate the energy necessary to establish the warp drive in CG.

The dominant energy condition (DEC) is reported in the literature ([24]) as for any , , or equivalently as assuming the WEC plus the additional condition that is a nonspacelike vector, that is, . It is easy to see that using as a vector the 4-velocity of the Eulerian observers [1], the previous condition for the DEC becomes .

Figure 5 illustrates the violation of the DEC for our standard solution (AWD with HSF and , , ). The plotted function is not negative everywhere, as required by the DEC, but shows a violation for the central portion of the warp bubble. Even if this energy condition appears to be violated, this does not notably affect the feasibility of our CG warp drive. We recall that the DEC is usually related to the standard perfect fluid stress-energy tensor, , where here and are the fluid density and pressure, while is the fluid 4-velocity. In this context, the DEC requires , but this condition is not required in general by all classical forms of matter [3]; therefore, its violation in our case is not particularly significant.

On the contrary, our standard solution also verifies the strong energy condition (SEC) which states that for all timelike vectors . Again, using the Eulerian 4-velocity vector in place of , the previous condition is equivalent to . Since the scalar is identically zero for all our solutions, as checked using our Mathematica program, the SEC is equivalent to , which is the WEC already verified in Section 3.

Finally, we want to estimate the energy necessary to establish our CG warp drive, under reasonable conditions. For this purpose, in Figure 6 we computed once again the energy density for our AWD with the Hartle shaping function (, , ), but this time for , a reasonable radius for a warp bubble enclosing our spaceship.

Figure 6 illustrates this solution, plotted only for , as this is the correct interval for the transverse coordinate . The cylindrical symmetry of this solution can also be better appreciated in this type of plot. We then followed the procedure outlined in [8, 11], to integrate the local energy density over the proper volume, in cylindrical coordinates at time over all space, obtaining the total energy as where is the determinant of the spatial metric on the constant time hypersurface. Since we assume that the spaceship is traveling at constant velocity, , the total energy is also constant with time. In the last equation, we reinstated a factor of to obtain the correct dimensions (see Section 2) and also inserted an overall multiplicative factor , which corresponds to the conformal gravity coupling constant in (9). This factor is necessary since our computation of in Figure 6 was done assuming .

Therefore, we need to know the CG value for in order to complete our energy estimation. Unfortunately, the value of this coupling constant is not well determined yet. The only value in the literature is reported by Mannheim [30] as since this coupling constant has the dimensions of action. Inserting this value for in (16), we obtain the energy estimate as where we also converted this energy into equivalent mass and compared our result with the solar mass .

The estimate in (18) would imply that an enormous amount of (standard) mass or energy is needed to establish our warp drive at a velocity equal to the speed of light, with a reasonable size for the warp bubble. However, the CG value of is not well established, since the number in (17) represents only an estimate of the macroscopic value of this coupling constant. This does not need to be the same as the microscopic value associated with the fundamental theory, which could be reduced by a factor of , where could be the number of occupied baryonic states in a galaxy (), or possibly the number of baryons in the universe () [41].

Therefore, our estimate could be reduced by many orders of magnitude. Moreover, the energy necessary to establish the warp drive might also be decreased by using a more efficient shaping function, an analysis which we leave for a future study on the subject.

5. Conclusions

In this paper, we have analyzed in detail the Alcubierre warp drive mechanism within the framework of conformal gravity. We have seen that a particular choice of the shaping function (Hartle shaping function, instead of the original Alcubierre one) can overcome the main limitation of the AWD in standard general relativity, namely, the violation of the weak energy condition.

In fact, we have shown that for a wide range of spaceship velocities, the CG solutions do not violate the WEC, and, therefore, the AWD mechanism might be viable, if CG is the correct extension of the current gravitational theories. All the components of the stress-energy tensor can be analytically calculated, using a Mathematica program based on conformal gravity. Thus, a warp drive can, at least in principle, be fully established following our computations.

We have also checked two other main energy conditions: the SEC is always verified, while the DEC is violated, at least in the case we considered. Finally, we estimated the energy needed to establish a reasonable warp drive at the speed of light. This energy depends critically on the value of , the CG coupling constant, which is not well known. Therefore, this estimate will need to be refined in future studies.


Energy Density Expression in Conformal Gravity

We present here the expression for the energy density in conformal gravity, computed using our Mathematica program. This is the general form of for any shaping function and its derivatives, up to the fourth order. The energy density is a function of coordinates , , , and, therefore, it has a cylindrical symmetry around the -axis. The distance is defined in (11), is the spaceship velocity, and is the conformal gravity coupling constant:

In the main part of our work, we used the Hartle shaping function in (13), or more explicitly as where indicates the positive part of the function. The derivatives of the HSF, up to the fourth order, were computed in terms of the Heaviside step function and the Dirac delta function , also using the relations for the derivatives of these special functions ; ; as follows:

By inserting these derivatives into (A.1), we obtained the expression used to compute in Figure 1 (bottom right panel) and Figures 2, 4, and 6.


This work was supported by a grant from the Frank R. Seaver College of Science and Engineering, Loyola Marymount University. The authors would like to acknowledge suggestions and clarifications by Dr. P. Mannheim and also thank the anonymous reviewers for the useful comments received.