Next Article in Journal
Equation of State of Strongly Magnetized Matter with Hyperons and Δ-Resonances
Previous Article in Journal
Squeezing the Parameter Space for Lorentz Violation in the Neutrino Sector with Additional Decay Channels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Curvature Invariants for the Accelerating Natário Warp Drive

1
Early Universe Cosmology and Strings (EUCOS) Group, Center for Astrophysics, Space Physics and Engineering Research (CASPER), Baylor University, Waco, TX 76798, USA
2
Department of Physics, Baylor University, Waco, TX 76798, USA
3
Institute for Advanced Studies at Austin, 11855 Research Blvd., Austin, TX 78759, USA
*
Author to whom correspondence should be addressed.
Submission received: 6 August 2020 / Revised: 3 September 2020 / Accepted: 10 September 2020 / Published: 16 September 2020

Abstract

:
A process for using curvature invariants is applied to evaluate the accelerating Natário warp drive. Curvature invariants are independent of coordinate bases and plotting the invariants is free of coordinate mapping distortions. While previous works focus mainly on the mathematical description of the warp bubble, plotting curvature invariants provides a novel pathway to investigate the Natário spacetime and its characteristics. For warp drive spacetimes, there are four independent curvature invariants the Ricci scalar, r 1 , r 2 , and w 2 . The invariant plots demonstrate how each curvature invariant evolves over the parameters of time, acceleration, skin depth and radius of the warp bubble. They show that the Ricci scalar has the greatest impact of the invariants on the surrounding spacetime. They also reveal key features of the Natário warp bubble such as a flat harbor in the center of it, a dynamic wake, and the internal structures of the warp bubble.
PACS:
04.20.-q; 04.20.Cv; 02.40.-k

1. Introduction

In Newtonian mechanics and special relativity, the velocity of any particle is fundamentally constrained by the speed of light, c. No particle can move through spacetime at a velocity greater than c in any reference frame. However, general relativity allows a global superluminal velocity while retaining a subluminal local velocity. Alcubierre demonstrated a solution for the Einstein field equations that allowed a spaceship to make a trip to a distant star in an arbitrarily short proper time [1]. He proposed a warped spacetime that pairs a local contraction of spacetime in front of the spaceship with a local expansion of spacetime behind the ship. While the spaceship remains within its own light cone and never exceeds c, globally the relative velocity, which is defined as proper spatial distance divided by proper time, may be much greater than c. Based on this principle, he named the faster-than-light (FTL) propulsion mechanism a “warp drive.”
FTL travel requires eight general conditions for a spaceship to carry a passenger safely [2]. First, the rocket equation is not required for travel by way of the warp. Second, the travel time through the FTL space warp to a distant star may be reduced to less than one year as seen both by the passengers in the warp and by stationary observers outside the warp. Third, the proper time as measured by the passengers should not be dilated by any relativistic effects. Fourth, any tidal-gravity accelerations acting on passengers inside of the spaceship may be reduced to be less than the acceleration of gravity near the Earth’s surface, g . Fifth, the local speed of the spaceship must be less than c. Sixth, the matter of the spaceship must not couple with any exotic material needed to generate the FTL space warp. Seventh, the FTL space warp should not generate an event horizon. Finally, the passengers should not encounter a singularity inside or outside of the FTL warp.
Traversable wormholes and warp drives are two known example spacetimes that satisfy these eight requirements [1,3,4,5,6,7,8,9,10,11]. While they are mathematical solutions to Einstein’s equations, building such devices are not achievable in the foreseeable future due to severe engineering constraints. In [12], the authors used the method of calculating and plotting curvature invariants to analyze several types of wormholes. In the present work, the authors will adapt this methodology to analyze the accelerating Natário Warp Drive.
Since Alcubierre, there has been considerable research into FTL warp drives. Krasnikov considered a non-tachyonic FTL warp bubble and showed it to be possible mathematically [6]. Van Den Broeck modified Alcubierre’s warp drive to have a microscopic surface area and a macroscopic volume inside. He showed that the modification allowed a warp bubble with more reasonable negative energy requirements of a few solar masses and that it has a more lenient violation of the null-energy-conditions (NEC). Later, Natário improved upon Alcubierre’s work by presenting a warp drive metric such that zero spacetime expansion occurs [3]. Instead of riding a contraction and expansion of spacetime, the warp drive may be observed to be “sliding” through the exterior spacetime at a constant global velocity. Loup expanded Natário’s work to encompass a global acceleration [4,5]. Finally, recent research has investigated the Einstein equations for the Alcubierre spacetime with a dust matter distribution as its source [13].
While the mathematics of a warp drive is well developed, mapping the spacetime around the warp drive remains unexplored until recently. Considering that a ship inside of a warp bubble is causally disconnected from the exterior, computer simulations of the surrounding spacetime are critical for the ship to map its journey and steer the warp bubble [3]. Alcubierre used the York time, which is defined as Θ = v s c x x s r s d f d r s , to map the volume expansion of a warp drive [1]. He plotted the York time to show how spacetime warped behind and in front of the spaceship. While the York time is appropriate when the 3-geometry of the hypersurfaces is flat, it will not contain all information about the curvature of spacetime in non-flat 3-geometries such as the accelerating Natátio warp drive spacetime. Alternatively, plotting the curvature invariants for a warp drive will display the spacetime’s curvature.
Christoffel proved that scalars constructed from the metric and its derivatives must be functions of the metric itself, the Riemann tensor, and its covariant derivatives [14]. In particular, curvature invariants are scalar products of Riemann, Ricci or Weyl tensors, or their covariant derivatives. Fourteen curvature invariants in ( 3 + 1 ) dimensions have been defined in the literature, but the total rises to seventeen when certain non-degenerate cases are taken into account [15]. Carminati and McLenaghan (CM) demonstrated a set of invariants that had several attractive properties. Their invariant set maintains general independence, requires each invariant to be of lowest possible degree, and contains a minimal independent set for any Petrov type and choice of Ricci tensor [16]. An invariant is considered to be independent if it cannot be written in terms of other members of the set that are of equal or lower degree. A set of invariants is independent if each element of the set is an independent invariant. The polynomial relationships, called syzgies, between the individual elements in a set of invariants may be used to determine independence. The syzgies of a set of invariants may be determined algorithmically. For the case of Class B spacetimes, the syzgies between the invariants further reduces the independent set of invariants to be only the four specific ones: R, r 1 , r 2 and w 2 [17]. Curvature invariants are the same regardless of your choice of coordinate. Consequently, a complicated spacetime, such as that of a warp drive, may be displayed without distortion when its curvature invariants are plotted.
Recent research by Overduin et al. [18] computed and plotted a set of independent curvature invariants for the hidden interiors of Kerr-Newman black holes. It produced visually stunning 3D plots of the surprisingly complex nature of spacetime curvature in Kerr-Newman black hole interiors. In addition, the calculation of curvature invariants in black holes is a rich topic of recent study [19,20,21,22,23]. This work motivated the present authors to undertake a similar study for the Natário warp drive. While the individual invariant functions are colossal in size and require a lengthy time to calculate, their plots can be quickly scanned and understood.

2. Method to Compute the Invariants

The CM curvature invariants can be calculated from any given line element. From the line element, the metric g α β is identified with the indices { α , β , γ , δ } ranging from { 0 , n 1 } , where n is the number of spacetime dimensions. The accelerating Natário warp drive was derived using n = 4 for its dimensions [4]. In spacetime, a complex null tetrad may be identified from a given line element. In this paper, a null tetrad ( l α , k α , m α , m ¯ α ) is found for the accelerating Natário metric. It is emphasized that the scalar invariants are independent of both the choice of coordinate parameterization and the choice of the null tetrad. A different choice of coordinates for each spacetime’s tetrad will result in invariants related by algebriac expressions to the ones plotted in this paper. From the metric, the affine connection Γ α β γ , Riemann tensor R α β γ δ , Ricci tensor R α β , Ricci scalar R, trace free Ricci tensor S α β and Weyl tensor C α β γ δ are calculated. The Newman-Penrose (NP) curvature components may be computed from the null tetrad, the trace-free Ricci tensor, and the Weyl tensor [24]. The thirteen different CM invariants are defined in [16]. Only four of these invariants are required by the syzgies for Class B spacetimes: the Ricci scalar (1), the first two Ricci invariants (2) (3), and the real component of the Weyl invariant J (4) [17]. In terms of the NP curvature coordinates, the CM invariants are:
R = g α β R α β ,
r 1 = 1 4 S α β S β α = 2 Φ 20 Φ 02 + 2 Φ 22 Φ 00 4 Φ 12 Φ 10 4 Φ 21 Φ 01 + 4 Φ 11 2 ,
r 2 = 1 8 S α β S γ α S β γ = 6 Φ 02 Φ 21 Φ 10 6 Φ 11 Φ 02 Φ 20 + 6 Φ 01 Φ 12 Φ 20 6 Φ 12 Φ 00 Φ 21 6 Φ 22 Φ 01 Φ 10 + 6 Φ 22 Φ 11 Φ 00 ,
w 2 = 1 8 C ¯ α β γ δ C ¯ α β ϵ η C ¯ γ δ ϵ η = 6 Ψ 4 Ψ 0 Ψ 2 6 Ψ 2 3 6 Ψ 1 2 Ψ 4 6 Ψ 3 2 Ψ 0 + 12 Ψ 2 Ψ 1 Ψ 3 .
The tetrad components of the traceless Ricci Tensor are Φ 00 through Φ 22 [24]. The complex tetrad components Ψ 0 to Ψ 5 are the six complex coefficients of the Weyl Tensor due to its tracelessness.

3. Warp Drive Spacetimes

Alcubierre and Natário developed warp drive theory using ( 3 + 1 ) dimensional ADM formalism. Spacetime is decomposed into space-like hypersurfaces parametrized by the value of an arbitrary time coordinate d x 0 [25,26]. Two nearby hypersurfaces, x 0 = constant and x 0 + d x 0 = constant, are separated by a proper time d τ = N ( x i , x 0 ) d x 0 . The ADM four-metric is
g α β = N 2 N i N j g i j N j N i g i j ,
where N is the lapse function, N i is the shift vector between hypersurfaces, and g i j is the 3-metric of the hypersurfaces with the indices { i j } ranging from { 1 , 3 } . A warp drive is defined as a globally hyperbolic spacetime ( M , g ) , where M = R 4 and g is given by the line element
d s 2 = d t 2 + i = 1 3 ( d x i X i d t ) 2 ,
for three unspecified bounded smooth functions ( X i ) = ( X , Y , Z ) in Cartesian coordinates [1,3]. The functions form a vector field given by X = X i x i = X x + Y y + Z z . X is a time-dependent vector field in Euclidean 3-space. Class B1 spacetimes include all spherical, planar and hyperbolic spacetimes [17]. Since the general warp drive spacetime is globally hyperbolic, any specific choice of X i will be a class B1 spacetime, and the complete set of curvature invariants will be Equation (1) through Equation (4).
Natário considered a specific choice of N and N i such that a net expansion/contraction of the surrounding spacetime is not necessary at all [3]. Instead, a volume-preserving warp drive “slides” the warp bubble region through space where the space in front of it is expanded and balanced by an opposite contraction of the space behind it. Originally, Natário only considered the warp bubble to slide at a constant superluminal velocity. Later, six line elements for the Natário spacetime metric with a constant acceleration were derived [4]. The specific equation for the Natário warp drive line element in the parallel covariant ( 3 + 1 ) dimensional ADM is
d s 2 = 1 2 X t + ( X t ) 2 ( X r s ) 2 ( X θ ) 2 d t 2 + 2 X r s d r s + X θ r s d θ d t d r s 2 r s 2 d θ 2 r s 2 sin 2 θ d ϕ 2 .
which is in the spherical coordinates, for which 0 r s < ; 0 θ π ; 0 φ 2 π , < t < , and a is the constant acceleration. The covariant shift vector components are given by
X t = 2 n ( r s ) a r s cos θ ,
X r s = 2 [ 2 n ( r s ) 2 + r s n ( r s ) ] a t cos θ ,
X θ = 2 n ( r s ) a t [ 2 n ( r s ) + r s n ( r s ) ] r s 2 sin θ .
The Natário warp drive continuous shape function is
n ( r s ) = 1 2 1 1 2 1 tanh [ σ ( r s ρ ) ] .
where σ is the skin depth of the warp bubble, and ρ is the radius of the warp bubble. Appendix A derives the comoving null tetrad for Equation (7). It is
l α = 1 2 1 X t + X r s 1 0 0 , k α = 1 2 1 X t X r s 1 0 0 , m α = 1 2 X θ 0 r i r sin θ , m ¯ α = 1 2 X θ 0 r i r sin θ .
The comoving null tetrad describes light rays traveling parallel to the warp bubble. The four CM invariants in Equation (1) through Equation (4) may be derived from Equations (7) and (12).

4. Invariants for the Accelerating Natário Warp Drive

The computations and plots for the four CM invariants in Equations (1)–(4) from the accelerating Natário warp drive are found by Mathematica® and presented herein. As a consequence of the choice of null tetrad (12), the plots are centered parallel to a spaceship’s flight. They represent the view from the spaceship’s bridge. The free parameters in the accelerating Natário warp drive are a, ρ , σ , and t in Equation (8) through Equation (11). To see the effect of the free parameters on the surrounding curvature invariants, each one was varied individually while maintaining all other parameters at constant values. Through this procedure, 684 plots were rendered for each invariant. For the remainder of the paper, a sample of these plots will be shared and inspected to see the effect of each parameter on the curvature invariants. Each subsection will include a sample of plots demonstrating the effect of each parameter and the remainder of the plots may be found in Appendix B.

4.1. Invariant Plots While Varying the Time

Figure 1 below shows how the Ricci scalar evolves from t = 0 s to t = 100 s, while setting the other free parameters to ρ = 100 m, σ = 50000 m 1 , and a = 1.0 m s 2 . In addition, the plots of the invariants r 1 , r 2 , and w 2 are included in Appendix B in Figure A1, Figure A2 and Figure A3. The figures provide rich details of the features in and around the warp bubble. Each plot has a safe harbor within ρ 100 m where the invariant’s magnitude is zero and the curvature is flat. This observation is consistent with a spaceship traveling at a velocity less than c in the interior of the harbor and experiencing only flat space throughout the entire time evolution. A symmetrical wake lies on both sides of the harbor and is characterized by large positive or negative magnitudes of the invariant. The warp bubble travels perpendicular to the wake from left to right in the plots. Its motion is driven by the constant curvature outside the harbor in the Ricci scalar. Initially, the curvature invariant has a small positive value in front of the harbor and a small negative value behind it as can be seen in Figure 1a. Afterwards, a constant negative curvature lies in front and behind of the harbor and wake. While the form of the Ricci scalar’s plots quickly reaches a constant shape, the Ricci scalar’s magnitude increases as time advances.
Choosing the same values for ρ , σ , and a as the Ricci scalar and varying the time, Figure A1 in Appendix B shows the time evolution for the r 1 invariant. It has many similar features to the Ricci scalar. It contains the safe harbor, a wake running perpendicular to the direction of motion, and its magnitude increases as time advances. The first difference is its positive magnitude and lack of internal structure in the wakes. The wake increases subtly in angular size as time increases. Finally, the space has no curvature in front and behind of the harbor and wake.
The invariant r 2 shares the same basic properties of the Ricci scalar and r 1 : the safe harbor, a wake, is asymptotically flat away from the bubble, and a linear increase with time. It is similar in shape to r 1 , but it increases in magnitude more drastically and has the same internal structure as the Ricci scalar.
The invariant w 2 has the same features as that of the Ricci scalar, r 1 , and r 2 : the safe harbor, a wake, and an increase in magnitude as time advances. Outside of the harbor and wake, it is asymptotically flat like the invariants r 1 and r 2 . The invariant’s wake contains a small amount of internal structure at lower time values, but as time reaches 100 s , the internal structures begins to form crenulations. As time continues to evolve, the crenulations travel out parallel to the length of the wake from the center of warp bubble. It can be speculated that this would cause an erratic flight path of the bubble as the crenulations’ ripple outwards. It can also be speculated that the interior structure of the Ricci scalar and invariant r 2 ’s wake would exhibit similar behavior at higher time values.
The shape of each invariant agrees with Natário’s description of the warp bubble sliding through spacetime. The positive and negative curvature from the wake will transport a ship in the harbor. The Ricci scalar has the greatest impact on the spacetime outside of the harbor and wake due to the constant curvature in this area. Each invariant’s magnitude experiences an increase as time advances. The invariant w 2 has the greatest impact on the wake as its magnitude increases the most as time advances. The wake’s shape also shows that there is internal structure to the warp bubble. Theoretically, the crenulations exhibited in the invariant w 2 would make navigation challenging. But, they can be avoided by accelerating the warp drive for only short periods. The increase suggests an engineering constraint on a maximum achievable global velocity.

4.2. Invariant Plots While Varying the Acceleration

Varying the acceleration of the invariants speeds up or slows down the time evolution in the previous section. Setting ρ = 100 m, σ = 50000 m 1 , and t = 1.0 s, Figure 2 below shows the variation of the acceleration of the Ricci scalar. The first plot of a = 0 m s 2 is consistent with the lapse function being zero for any given time period. The distance between hypersurfaces will be constant. The space will be flat and no warp bubble will form as is shown. The second plot corresponds to a time slice between Figure 1a,b. Similarly, the third plot is identical to Figure 1c, and the fourth plot corresponds to a time slice after Figure 1d. It can be concluded that modifying the acceleration parameter corresponds with modifying the rate of change of the hypersurfaces. Additionally, our analysis of the invariant’s shape and properties in Section 4.1 holds for the acceleration plots.
The plots of the invariants r 1 , r 2 and w 2 follow a similar process as the Ricci scalar. They are plotted in Figure A4, Figure A5 and Figure A6 in Appendix B. When a = 0 , their plots are identical to Figure 2a. Then, they can be seen as additional time slices between the plots shown in Figure A1a,b for r 1 and Figure A2c,d for r 2 . Some additional features are present in the plots. The invariants do warp themselves much more significantly and non-symmetricly than their counterparts for the Ricci scalar, as can be seen in Figure A5c.

4.3. Invariant Plots While Varying the Skin Depth of the Warp Bubble

Varying the skin depth of the warp bubble, σ , does not noticeably affect the invariant plots. Figure 3 below presents the plots for each of the four invariants while doubling the skin depth and setting ρ = 100 m, t = 1.0 s, and a = 1.0 m s 2 . Figure A7 in Appendix B presents the remaining invariants r 1 , r 2 , and w 2 , which have similar features to the Ricci Scalar. The shape of the invariant does not change in the figures after doubling. In the Natário line element, the skin depth occurs in tanh σ ( r s ρ ) of (11). σ moderates the rate at which the top hat function converges to be 0 outside the warp radius and 1 along the depth of warp’s edge near ρ . For the relatively high values of σ plotted, the top hat function must approach 1, and the overall depth of the warp bubble must be very small. If σ was plotted on the order of 1 m 1 , a greater effect on the invariants would be observed than what is shown with the attached plots.

4.4. Invariant Plots While Varying the Radius of the Warp Bubble

Varying the radius of the warp bubble, ρ , increases the size of the safe harbor inside each of the invariants. Figure 4 above presents the plots for each of the four invariants while doubling the skin depth and setting σ = 50,000 m 1 , a = 1.0 m s 2 , and t = 1 s. Figure A7 in Appendix B present the remaining invariants r 1 , r 2 , and w 2 , which have similar features to the Ricci Scalar. In the figures for each invariant, the radial coordinate doubles in radial size without affecting the shape of the plots. The safe harbor of ρ 100 m in the left hand column also doubles in size to ρ 200 m. The only other pertinent feature is in the internal structure of w 2 . The structures are reduced implying that they cluster near the center.

5. Conclusions

This paper demonstrates how computing and plotting curvature invariants for various parameters of warp drive spacetimes can reveal their underlying features. While the individual functions require pages to express, their plots can be easily read and understood. The plots give the magnitude of curvature at each point around the ship. Where the plots’ magnitudes are large, space is greatly warped and vice versa. Additionally, the rate at which spacetime is being folded may be approximated by observing the changes in slopes on the plots. This information can be used to map the spacetime around the ship and potentially aid in navigation.
In this research, the accelerating Natário warp drive metric was plotted for different choices of the parameters t , a , ρ , σ . In all plots of the listed parameters, the constant features include a safe harbor and a wake running perpendicular to the warp bubble’s direction of motion. The wake will carry the spaceship along as Natário predicted. The Ricci scalar revealed constant curvature outside of the harbor and wake. The other invariants had asymptotically flat space outside of the harbor and wake. As a consequence, the Ricci Scalar will have the greatest impact of the invariants on the spacetime outside of the harbor and wake. By varying time, the plot of each invariant experiences a sudden jump from positive curvature in the direction of motion to negative curvature. As time progresses, the shape of the R, r 1 , and r 2 invariants remains constant, but the magnitude of the invariants increases. Also, the angular arc widens slightly for the invariants r 1 and r 2 as time increases. After 100 s, the invariant w 2 begins to exhibit crenulations in the interior of the wake. Because the greatest invariant magnitudes occur for the invariant w 2 , it must have the greatest impact on the wake and the internal structures. By varying the acceleration, the invariant plots skip through the time slices and the internal structures become more prominent. Changing the skin depth did not change either the shape or magnitude of the invariant plots. As expected, doubling the radius also doubled the size of the safe harbor in the invariants plots without affecting the shape of the warp bubble. The invariant plots give a rich and detailed understanding of the warp bubble’s curvature.
Computing and plotting the invariant functions have significant advantages for the inspection of warp drives and their potential navigation. As mentioned previously, plotting the invariants has the advantage that they are free from coordinate mapping distortions, divergences, discontinuities or other artifacts of the chosen coordinates. Once the plots of the invariant functions reveal the location of any artifacts, their position can be related mathematically to the standard tensors, and their effect on an objects motion can then be analyzed. The invariant plots properly illustrate the entire underlying spacetime independent of a chosen coordinate system. A second advantage is the relative ease with which the invariants can be plotted. Software packages exist or can be developed to calculate the standard tensors. The aforementioned tensors lead to a chosen basis of invariants. While the CM invariants were chosen in this paper, other sets of invariants exist, such as the Cartan invariants and the Witten and Petrov invariants [23,27]. It is an open problem to inspect the curvature of the warp drive spacetimes in these invariant sets. It is expected that the main features identified in this paper will also hold in these different bases.
In addition to inspecting different invariant bases, further work can be done in mapping warp drive spacetimes. The work in this paper can be further expanded to greater time slices for each invariant. Potentially, crenulations like the ones observed for w 2 also exist within the internal structures observed in the wakes of the other invariants. As the crenulations ripple, navigating the warp drive will be increasingly difficult. To minimize the effect of the crenulations, the warp drive should be periodically turned on and off as it accelerates. It can be speculated that the increase in the magnitude of the curvature as time advances establishes a relationship between the magnitude of the acceleration and the amount of energy required to accelerate the warp bubble to arbitrarily high velocities. Since the acceleration skips through the warp drive’s time slices, greater values of the acceleration will reach high magnitudes of each invariant and curvature more quickly. This relationship suggests that a realistic warp drive would only be able to accelerate to some finite velocity that is potentially greater than c, and it is evidence for a less restrictive superluminal censorship theorem than previously considered in [28]. Instead, it is in accordance with the constraint that the net total energy stored in the warp bubble must be less than the total rest energy of the spaceship itself [29]. In addition, the crenulations imply that high frequency gravitational waves would be produced by an accelerating warp drive. Potentially, these waves could be detected by an extremely sensitive detector. A detector like the one proposed by [30] would be suitable. Another potential branch of the research presented in this paper is to calculate the geodesic equation for the accelerating Natário warp drive. The Christoffel symbols Γ i j k were calculated in the process to find the curvature invariants. It remains to plug them into the geodesic equation and solve given the initial conditions of being near Earth and at rest. Finally, the technique of plotting the invariants can be applied to other warp drive spacetimes such as Alcubierre, Krasnikov, Van Den Broeck, and Natário at a constant velocity.

Author Contributions

The individual contributions of the authors of this paper are: conceptualization, E.W.D. and G.B.C.; methodology, E.W.D. and B.M.; software, B.M.; validation, B.M., C.K.W., M.D.A., A.B., and C.E.; formal analysis, B.M. and E.W.D.; investigation, B.M.; resources, B.S., E.W.D. and G.B.C.; data curation, B.M.; writing—original draft preparation, B.M.; writing—review and editing, B.M., A.K., W.J., M.G., C.K.W., J.S.L., E.W.D., and G.B.C.; visualization, B.M.; supervision, E.W.D. and G.B.C.; project administration, E.D and G.B.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank F. Loup and D.D. McNutt for beneficial discussions. The authors would like to thank the reviewers at MDPI particles.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Null Vectors of the Natário Metric

A null tetrad contains two real null vectors, k and l , and two complex conjugate null vectors, m and m ¯ that satisfy the following algebraic relationships [24]:
e α = ( l , k , m , m ¯ ) ,
g α β = 2 m ( α m ¯ β ) 2 k ( α l β ) = 0 1 0 0 1 0 0 0 0 0 0 1 0 0 1 0 .
Given an orthonormal tetrad, E ff , it can be related to a complex null tetrad Equation (A1) by:
l α = 1 2 ( E 1 + E 2 ) , k α = 1 2 ( E 1 E 2 ) , m α = 1 2 ( E 3 + i E 4 ) , m ¯ α = 1 2 ( E 3 i E 4 ) .
The line element in Equation (7) has an orthonormal tetrad:
E 1 = 1 X t 0 0 0 , E 2 = X r s 1 0 0 , E 3 = X θ 0 r 0 , E 4 = 0 0 0 r sin θ .
Using Mathematica®, (A2) and (A4) satisfy the relationship
g α β = E α · E β .
By applying Equation (A3) to (A4), the null vectors in Equation (12) result.

Appendix B. Invariant Plots

Figure A1. The time, t, of the invariant r1 is plotted in the figures above. The other parameters were set to be σ = 50,000 m−1, ρ = 100 m, and a = 1.0 m s−2. The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the plot corresponds to a distance of approximately 33 m. The constant portion to the right of the harbor is to the front of the warp bubble, and the constant portion to the left of the harbor is the back of the warp bubble.
Figure A1. The time, t, of the invariant r1 is plotted in the figures above. The other parameters were set to be σ = 50,000 m−1, ρ = 100 m, and a = 1.0 m s−2. The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the plot corresponds to a distance of approximately 33 m. The constant portion to the right of the harbor is to the front of the warp bubble, and the constant portion to the left of the harbor is the back of the warp bubble.
Particles 03 00042 g0a1
Figure A2. The time, t, of the invariant r2 is plotted in the figures above. The other parameters were set to be σ = 50,000 m−1, ρ = 100 m, and a = 1.0 m s−2. The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the plot corresponds to a distance of approximately 33 m. The constant portion to the right of the harbor is to the front of the warp bubble, and the constant portion to the left of the harbor is the back of the warp bubble.
Figure A2. The time, t, of the invariant r2 is plotted in the figures above. The other parameters were set to be σ = 50,000 m−1, ρ = 100 m, and a = 1.0 m s−2. The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the plot corresponds to a distance of approximately 33 m. The constant portion to the right of the harbor is to the front of the warp bubble, and the constant portion to the left of the harbor is the back of the warp bubble.
Particles 03 00042 g0a2aParticles 03 00042 g0a2b
Figure A3. The time, t, of the invariant w2 is plotted in the figures above. The other parameters were set to be σ = 50,000 m−1, ρ = 100 m, and a = 1.0 m s−2. The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the plot corresponds to a distance of approximately 33 m. The constant portion to the right of the harbor is to the front of the warp bubble, and the constant portion to the left of the harbor is the back of the warp bubble.
Figure A3. The time, t, of the invariant w2 is plotted in the figures above. The other parameters were set to be σ = 50,000 m−1, ρ = 100 m, and a = 1.0 m s−2. The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the plot corresponds to a distance of approximately 33 m. The constant portion to the right of the harbor is to the front of the warp bubble, and the constant portion to the left of the harbor is the back of the warp bubble.
Particles 03 00042 g0a3aParticles 03 00042 g0a3b
Figure A4. The acceleration, a, of the invariant r1 is plotted in the figures above. The other parameters were set to be σ = 50,000 m−1, ρ = 100 m, and t = 1.0 s. The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the plot corresponds to a distance of approximately 33 m. The constant portion to the right of the harbor is to the front of the warp bubble, and the constant portion to the left of the harbor is the back of the warp bubble.
Figure A4. The acceleration, a, of the invariant r1 is plotted in the figures above. The other parameters were set to be σ = 50,000 m−1, ρ = 100 m, and t = 1.0 s. The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the plot corresponds to a distance of approximately 33 m. The constant portion to the right of the harbor is to the front of the warp bubble, and the constant portion to the left of the harbor is the back of the warp bubble.
Particles 03 00042 g0a4
Figure A5. The acceleration, a, of the invariant r2 is plotted in the figures above. The other parameters were set to be σ = 50,000 m−1, ρ = 100 m, and t = 1.0 s. The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the plot corresponds to a distance of approximately 33 m. The constant portion to the right of the harbor is to the front of the warp bubble, and the constant portion to the left of the harbor is the back of the warp bubble.
Figure A5. The acceleration, a, of the invariant r2 is plotted in the figures above. The other parameters were set to be σ = 50,000 m−1, ρ = 100 m, and t = 1.0 s. The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the plot corresponds to a distance of approximately 33 m. The constant portion to the right of the harbor is to the front of the warp bubble, and the constant portion to the left of the harbor is the back of the warp bubble.
Particles 03 00042 g0a5
Figure A6. The acceleration, a, of the invariant w2 is plotted in the figures above. The other parameters were set to be σ = 50,000 m−1, ρ = 100 m, and t = 1.0 s. The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the plot corresponds to a distance of approximately 33 m. The constant portion to the right of the harbor is to the front of the warp bubble, and the constant portion to the left of the harbor is the back of the warp bubble.
Figure A6. The acceleration, a, of the invariant w2 is plotted in the figures above. The other parameters were set to be σ = 50,000 m−1, ρ = 100 m, and t = 1.0 s. The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the plot corresponds to a distance of approximately 33 m. The constant portion to the right of the harbor is to the front of the warp bubble, and the constant portion to the left of the harbor is the back of the warp bubble.
Particles 03 00042 g0a6aParticles 03 00042 g0a6b
Figure A7. The warp bubble skin depth, σ, of the invariants r1, r2 and w2 are plotted in the figures above. The other parameters were set to be ρ = 100 m, and a = 1.0 m s™2, and t = 1 s. The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the plot corresponds to a distance of approximately 33 m. The constant portion to the right of the harbor is to the front of the warp bubble, and the constant portion to the left of the harbor is the back of the warp bubble.
Figure A7. The warp bubble skin depth, σ, of the invariants r1, r2 and w2 are plotted in the figures above. The other parameters were set to be ρ = 100 m, and a = 1.0 m s™2, and t = 1 s. The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the plot corresponds to a distance of approximately 33 m. The constant portion to the right of the harbor is to the front of the warp bubble, and the constant portion to the left of the harbor is the back of the warp bubble.
Particles 03 00042 g0a7aParticles 03 00042 g0a7b
Figure A8. The warp bubble radius ρ, of the invariants r1, r2 and w2 are plotted in the figures above. The other parameters were set to be σ = 50,000 m−1, a = 1.0 m s−2, and t = 1 s. The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the plot corresponds to a distance of approximately 33 m. The constant portion to the right of the harbor is to the front of the warp bubble, and the constant portion to the left of the harbor is the back of the warp bubble.
Figure A8. The warp bubble radius ρ, of the invariants r1, r2 and w2 are plotted in the figures above. The other parameters were set to be σ = 50,000 m−1, a = 1.0 m s−2, and t = 1 s. The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the plot corresponds to a distance of approximately 33 m. The constant portion to the right of the harbor is to the front of the warp bubble, and the constant portion to the left of the harbor is the back of the warp bubble.
Particles 03 00042 g0a8aParticles 03 00042 g0a8b

References

  1. Alcubierre, M. The Warp drive: Hyperfast travel within general relativity. Class. Quantum Gravity 1994, 11, L73. [Google Scholar] [CrossRef]
  2. Davis, E.W. Chapter 15: Faster-Than-Light Approaches in General Relativity. In Frontiers of Propulsion Science; (2nd Printing with Corrections, 2012); Progress in Astronautics & Aeronautics Series); American Institution of Aeronautics & Astronautics Press: Reston, VA, USA, 2009; Volume 227, pp. 473–509. [Google Scholar]
  3. Natário, J. Warp drive with zero expansion. Class. Quantum Gravity 2002, 19, 1157. [Google Scholar] [CrossRef]
  4. Loup, F. An Extended Version of the Natário Warp Drive Equation Based in the Original 3 + 1 ADM Formalism Which Encompasses Accelerations and Variable Velocities. Ph.D. Thesis, Residencia de Estudantes Universitas, Lisboa, Portugal, 2017. [Google Scholar]
  5. Loup, F. Six Different Natario Warp Drive Spacetime Metric Equations; Research Report; Residencia de Estudantes Universitas: Lisboa, Portugal, 2018. [Google Scholar]
  6. Krasnikov, S.V. Hyperfast travel in general relativity. Phys. Rev. D. 1998, 57, 4760. [Google Scholar] [CrossRef] [Green Version]
  7. Van Den Broeck, C. A ’Warp drive’ with reasonable total energy requirements. Class. Quantum Gravity 1999, 16, 3973. [Google Scholar] [CrossRef] [Green Version]
  8. Morris, M.S.; Thorne, K.S. Wormholes in space-time and their use for interstellar travel: A tool for teaching general relativity. Am. J. Phys. 1988, 56, 395. [Google Scholar] [CrossRef] [Green Version]
  9. Morris, M.S.; Thorne, K.S.; Yurtsever, U. Wormholes, time machines, and the weak energy conditions. Phys. Rev. Lett. 1988, 61, 1446–1449. [Google Scholar] [CrossRef] [Green Version]
  10. Visser, M. Lorentzian Wormholes: From Einstein to Hawking; AIP Press: New York, NY, USA, 1995. [Google Scholar]
  11. Lobo, F.S.N. Wormhole Basics. In Wormholes, Warp Drives and Energy Conditions; Lobo, F.S.N., Ed.; Springer: Cham, Switzerland, 2017; pp. 11–33. [Google Scholar] [CrossRef]
  12. Mattingly, B.; Kar, A.; Julius, W.; Gorban, M.; Watson, C.; Ali, M.D.; Baas, A.; Elmore, C.; Shakerin, B.; Davis, E.W.; et al. Curvature Invariants for Lorentzian Traversable Wormholes. Universe 2020, 6, 11. [Google Scholar] [CrossRef] [Green Version]
  13. Santos-Pereira, O.L.; Abreu, E.M.C.; Ribeiro, M.B. Dust content solutions for the Alcubierre warp drive spacetime. Eur. Phys. J. C 2020, 80, 786. [Google Scholar] [CrossRef]
  14. Christoffel, E.B. Ueber die Transformation der homogenen Differentialausdrücke zweiten Grades. J. Die Reine Angew. Math. 1869, 70, 46–70. [Google Scholar]
  15. Zakhary, E.; McIntosh, C.B.G. A Complete Set of Riemann Invariants. Gen. Relat. Gravit. 1997, 29, 539–581. [Google Scholar] [CrossRef]
  16. Carminati, J.; McLenaghan, R.G. Algebraic invariants of the Riemann tensor in a four-dimensional Lorentzian space. J. Math. Phys. 1991, 32, 3135–3140. [Google Scholar] [CrossRef] [Green Version]
  17. Santosuosso, K.; Pollney, D.; Pelavas, N.; Musgrave, P.; Lake, K. Invariants of the Riemann tensor for class B warped product spacetimes. Comput. Phys. Commun. 1998, 115, 381. [Google Scholar] [CrossRef] [Green Version]
  18. Overduin, J.; Coplan, M.; Wilcomb, K.; Henry, R.C. Curvature Invariants for Charged and Rotating Black Holes. Universe 2020, 6, 22. [Google Scholar] [CrossRef] [Green Version]
  19. Henry, R.C. Kretshmann Scalar for a Kerr-Newman Black Hole. Astrophys. J. 2000, 535, 350–353. [Google Scholar] [CrossRef] [Green Version]
  20. Baker, J.G.; Campanelli, M. Making use of geometrical invariants in black hole collisions. Phys. Rev. D 2000, 62, 127501. [Google Scholar] [CrossRef] [Green Version]
  21. Abdelqader, M.; Lake, K. Invariant characterization of the Kerr spacetime: Locating the horizon and measuring the mass and spin of rotating black holes using curvature invariants. Phys. Rev. D 2015, 91, 084017. [Google Scholar] [CrossRef] [Green Version]
  22. MacCallum, M.A.H. Spacetime Invariants and Their Uses. arXiv 2015, arXiv:1504.06857v1. [Google Scholar]
  23. Brooks, D.; MacCallum, M.A.H.; Gregoris, D.; Forget, A.; Coley, A.A.; Chavy-Waddy, P.C.; McNutt, D.D. Cartan Invariants and Event Horizon Detection, Extended Version. Gen. Relat. Gravit. 2018, 50, 37. [Google Scholar] [CrossRef]
  24. Stephani, H.; Kramer, D.; MacCallum, M.A.H.; Hoenselaers, C.; Herlt, E. Exact Solutions of Einstein’s Field Equations; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar] [CrossRef]
  25. Arnowitt, R.L.; Deser, S.; Misner, C.W. Dynamical Structure and Definition of Energy in General Relativity. Phys. Rev. 1959, 116, 1322. [Google Scholar] [CrossRef]
  26. Marquet, P. The Generalized Warp Drive Concept in the EGR Theory. Abraham Zelmanov J. 2009, 2, 261–287. [Google Scholar]
  27. Bronnikov, K.A. 1973 Scalar-tensor theory and scalar charge. Acta Phys. Pol. 2009, B4, 251–266. [Google Scholar]
  28. Barcelo, C.; Visser, M. Twilight for the energy conditions? Int. J. Mod. Phys. D 2002, 11, 1553. [Google Scholar] [CrossRef] [Green Version]
  29. Lobo, F.S.N.; Visser, M. Fundamental limitations on ’warp drive’ spacetimes. Class. Quantum Gravity 2004, 21, 5871–5892. [Google Scholar] [CrossRef] [Green Version]
  30. Woods, R.C.; Baker, R.M.L.; Li, F.; Stephenson, G.V.; Davis, E.W.; Beckwith, A.W. A new theoretical technique for the measurement of high frequency relic gravitational waves. J. Mod. Phys. 2011, 2, 498–518. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The time, t, of the Ricci Scalar R is plotted in the figures above. The other parameters were set to be σ = 50,000 m 1 , ρ = 100 m, and a = 1.0 m s 2 . The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the plot corresponds to a distance of approximately 33 m. The portion of small negative curvature to the right of the harbor is to the front of the warp bubble, and the portion to the left of the harbor is the back of the warp bubble.
Figure 1. The time, t, of the Ricci Scalar R is plotted in the figures above. The other parameters were set to be σ = 50,000 m 1 , ρ = 100 m, and a = 1.0 m s 2 . The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the plot corresponds to a distance of approximately 33 m. The portion of small negative curvature to the right of the harbor is to the front of the warp bubble, and the portion to the left of the harbor is the back of the warp bubble.
Particles 03 00042 g001
Figure 2. The acceleration, a, of the Ricci Scalar R is plotted in the figures above. The other parameters were set to be σ = 50,000 m 1 , ρ = 100 m, and t = 1 s. The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the plot corresponds to a distance of approximately 33 m. The portion of small negative curvature to the right of the harbor is to the front of the warp bubble, and the portion to the left of the harbor is the back of the warp bubble.
Figure 2. The acceleration, a, of the Ricci Scalar R is plotted in the figures above. The other parameters were set to be σ = 50,000 m 1 , ρ = 100 m, and t = 1 s. The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the plot corresponds to a distance of approximately 33 m. The portion of small negative curvature to the right of the harbor is to the front of the warp bubble, and the portion to the left of the harbor is the back of the warp bubble.
Particles 03 00042 g002
Figure 3. The warp bubble skin depth, σ, of the Ricci Scalar R is plotted in the figures above. The other parameters were set to be ρ = 100 m, a = 1.0 m s−2 and t = 1 s. The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the Z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The X and Y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line left hand column corresponds to a distance of approximately 33 m. The constant portion to the right of the harbor is to the front of the warp bubble, and the constant portion to the left of the harbor is the back of the warp bubble.
Figure 3. The warp bubble skin depth, σ, of the Ricci Scalar R is plotted in the figures above. The other parameters were set to be ρ = 100 m, a = 1.0 m s−2 and t = 1 s. The plots were rendered using Mathematica’s® RevolutionPlot3D. It plots a function by rotating the function around the Z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The X and Y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line left hand column corresponds to a distance of approximately 33 m. The constant portion to the right of the harbor is to the front of the warp bubble, and the constant portion to the left of the harbor is the back of the warp bubble.
Particles 03 00042 g003
Figure 4. The warp bubble radius, ρ, of the Ricci Scalar R is plotted in the figures above. The other parameters were set to be σ = 50,000 m 1 , a = 1.0 m s™2, and t = 1 s. The plots were rendered using Mathematica∙s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the left hand column corresponds to a distance of approximately 33 m and each radial line in the right hand column corresponds to a distance of approximately 67 m. The constant portion to the right of the harbor is to the front of the warp bubble, and the constant portion to the left of the harbor is the back of the warp bubble.
Figure 4. The warp bubble radius, ρ, of the Ricci Scalar R is plotted in the figures above. The other parameters were set to be σ = 50,000 m 1 , a = 1.0 m s™2, and t = 1 s. The plots were rendered using Mathematica∙s® RevolutionPlot3D. It plots a function by rotating the function around the z-axis. The amount that the space diverges from being flat is represented by the invariant function’s magnitude, and its amount is labeled on the vertical axis. The x and y-axis are displayed on the plots, and they display the distance from a spaceship in the flat portion at the center of each figure. Each radial line in the left hand column corresponds to a distance of approximately 33 m and each radial line in the right hand column corresponds to a distance of approximately 67 m. The constant portion to the right of the harbor is to the front of the warp bubble, and the constant portion to the left of the harbor is the back of the warp bubble.
Particles 03 00042 g004

Share and Cite

MDPI and ACS Style

Mattingly, B.; Kar, A.; Gorban, M.; Julius, W.; Watson, C.K.; Ali, M.D.; Baas, A.; Elmore, C.; Lee, J.S.; Shakerin, B.; et al. Curvature Invariants for the Accelerating Natário Warp Drive. Particles 2020, 3, 642-659. https://0-doi-org.brum.beds.ac.uk/10.3390/particles3030042

AMA Style

Mattingly B, Kar A, Gorban M, Julius W, Watson CK, Ali MD, Baas A, Elmore C, Lee JS, Shakerin B, et al. Curvature Invariants for the Accelerating Natário Warp Drive. Particles. 2020; 3(3):642-659. https://0-doi-org.brum.beds.ac.uk/10.3390/particles3030042

Chicago/Turabian Style

Mattingly, Brandon, Abinash Kar, Matthew Gorban, William Julius, Cooper K. Watson, M.D. Ali, Andrew Baas, Caleb Elmore, Jeffrey S. Lee, Bahram Shakerin, and et al. 2020. "Curvature Invariants for the Accelerating Natário Warp Drive" Particles 3, no. 3: 642-659. https://0-doi-org.brum.beds.ac.uk/10.3390/particles3030042

Article Metrics

Back to TopTop