Tractor beams, pressor beams, and stressor beams in general relativity

The metrics of general relativity generally fall into two categories: Those which are solutions of the Einstein equations for a given source energy-momentum tensor, and the"reverse engineered"metrics -- metrics bespoke for a certain purpose. Their energy-momentum tensors are then calculated by inserting these into the Einstein equations. This latter approach has found frequent use when confronted with creative input from fiction, wormholes and warp drives being the most famous examples. In this paper, we shall again take inspiration from fiction, and see what general relativity can tell us about the possibility of a gravitationally induced tractor beam. We will base our construction on warp drives and show how versatile this ansatz alone proves to be. Not only can we easily find tractor beams (attracting objects); repulsor/pressor beams are just as attainable, and a generalization to"stressor"beams is seen to present itself quite naturally. We show that all of these metrics would violate various energy conditions. This will provide an opportunity to ruminate on the meaning of energy conditions as such, and what we can learn about whether an arbitrarily advanced civilization might have access to such beams.

To the best of our knowledge, no really focussed work has been carried out on putting tractor/pressor beams into a coherent general relativistic context. (Acoustic tractor beams [24][25][26][27], matter wave tractor beams [28], or optical tweezers [29], seem to be the closest one gets in the current scientific literature.) Herein we shall analyze tractor/pressor/stressor beams from a general relativistic perspective. The basic idea is to significantly modify and adapt the "warp drive" spacetimes [13][14][15][16][17][18] in a suitable manner, giving them a "beam like" profile, and analysing the induced stresses and forces. Instead of a spaceship riding inside a warp bubble, we will assume that the warp field is in the form of a "beam" generated to pull/repel a target. The mechanisms by which this field is generated is beyond the scope of this article. We will assume that some arbitrarily advanced civilisation [30,31] might have developed the appropriate beam generation technology.
Specifically, we shall assume for convenience that the modified warp drive space-times are oriented in the z direction and give them a uniform transverse profile in the x and y directions, typically of the form f (x 2 + y 2 ). Doing so, one obtains a "beam" rather than a "warp bubble". Note that in this work we will let the (t, z) dependence remain arbitrary.
As always, when working in this area of speculative physics, including wormholes, and warp drives, and now tractor/pressor/stressor beams, a major justification for undertaking this exercise is to push general relativity to the breaking point; in the hope that the resulting wreckage will tell us something interesting -possibly even about quantum gravity [5,16].
After first analysing Natário's generic warp drive case [14], we will consider three special cases: 1. We modify the Alcubierre fixed-flow-direction warp field. 2. We modify the Natário zero-expansion warp field. 3. We modify the zero-vorticity warp field.
We shall also illustrate each of these three cases with some specific examples based on beams with a Gaussian profile.
A recurring theme in the analysis will be the use of the classical point-wise energy conditions (null, weak, strong, and dominant; abbreviated NEC, WEC, SEC, and DEC, respectively) [32][33][34][35]. They can be considered as an attempt to remain as agnostic as possible about underlying equations of state. While the energy conditions do not seem to be fundamental physics, they are at the very least a very good sanity check on just how weird the physics is getting [33,36,37]. We already know of examples of violations at microscopic scales (e.g., Hawking radiation) and mesoscopic scales (e.g., Casimir effect). No macroscopic violations of the energy conditions are known up to this point, except at truly cosmological scales -and they violate only some of the energy conditions (the accelerated expansion of the universe violates the strong and dominant energy conditions, but not the null and weak energy conditions [38][39][40][41]). Therefore, besides the violation of the energy conditions not being an absolute prohibition, it is an indication that one should look very carefully at the underlying physics [36,37]. For more background on the energy conditions see .
For the sake of full transparency, we should also mention that our interest in these topics was rekindled and inspired by three recent papers [70][71][72]. Unfortunately, significant parts of those three papers are incorrect, misguided, and/or misleading. See reference [18] for details.

When things need to be moved
One of Wheeler's adages that became standard general relativity folklore is the famous saying that "space-time tells matter how to move; matter tells space-time how to curve". From many a practical point of view, questions regarding objects' movement are less about the how and more about the ought -things are wanted elsewhere from where they are now. It is this logistical perspective that we shall address in the following: How can we ensure that general relativity does the job of moving an object (like a cow, [preferably a spherical cow, in vacuum], or a Corellian CR90 corvette) for us?
The key ingredient will be to limit ourselves to test field cases, where we neglect the mass of the objects we want to move, how they interact with space-time and with the matter we put in space-time to move them. This reduces the core physics question to one of forces: We want to use the pressures encoded in the stress-energy tensor of a beam-like field to move target test masses.
-3 - The field is assumed to be sourced by someone on the left at negative z, the target-a flat cow in the tractor field space-time-on the right at positive z. Choosing the source and target provides for a distinction between tractor and pressor (or repulsor) fields. Details concerning this particular beam configuration can be found in section 6.2.1. The parameters of equation (6.37) that we have chosen are: A = 0.5, B = C = 1.0. The purple line in the density plot for the zero-vorticity beam indicates the location where the energy density is zero.
The primary force-related calculation we shall undertake is this: If the beam is pointed in the z direction, then one calculates the stress-energy component T zz (t, x, y, z), and integrates it over the entire transverse x-y plane to find the net force: Here the + sign corresponds to a beam impinging on the target from the left, whereas the − sign corresponds to a beam impinging on the target from the right. There is an approximation being made here, that the beam is narrow with respect to the target, so that it is a good approximation to integrate over the entire transverse x-y plane. If the beam is instead wide compared to the size of the target then one should instead use the approximation Here T zz (t, 0, 0, z) is the on-axis stress, and A is the cross sectional area of the target. For a beam of intermediate widths, (comparable to the size of the target), one would in principle need to calculate but this is unnecessarily complicated for the primary issues we wish to address.
The quantity F (t, z) is the net force the beam exerts on some target located at some position z at time t. For convenience we shall henceforth assume that the field is generated by someone positioned on the left, and that the target be positioned to the right of the generator, (see Figure 1), thus allowing us to restrict attention to the plus sign in equations (2.1) -(2.3). We shall furthermore assume that the target will move under the influence of the field, while the "generator" will not, and-as mentioned above-both behave as test fields. This setup provides for a simple characterization of the effect of the field: If F (t, z) < 0, corresponding to attraction, we call this a tractor beam. If F (t, z) > 0, corresponding to repulsion, we call this a pressor beam 1 . On the other hand, the definition of a stressor beam can be a little trickier. The reason for this being that, independent of the overall sign of F (t, z), one can quite generally define a beam which has significantly varying pressure across the cross-sectional area of the target. In this way, there might be a certain ambiguity about when a specific beam would be considered to be a tractor/pressor or a stressor beam, since this would depend on the properties of the target material -such as its elasticity and ultimate yield strength, and so on.
However, for most "applications", we expect the T zz component for a tractor (pressor) beam to not vary too greatly over its region of influence on the target. A quick measure of when a beam would behave as a stressor beam is given by: Here A is the cross-sectional area of the target exposed to the beam, and σ material is the ultimate yield stress of the material making up the target.
While equations (2.1)-(2.3) are universally valid, both for standard general relativity, and for modified theories of gravity, we will focus mainly on standard general relativity. Therefore, using the Einstein equations, we have for a narrow beam in terms of the Einstein tensor: while for a wide beam These are the key equations we will be using in the following sections.
As usual, we are using geometrodynamic units, where G Newton → 1 and c → 1. If one wishes to reinstate SI units, then in terms of the Stoney force and It is worthwhile mentioning that the magnitude of the Stoney force is truly enormous -some 1.2 × 10 44 Newtons. Accordingly, relatively small spacetime curvatures (weakfield gravity) can still lead to significant human-scale forces and stresses. It is beyond the scope of the present article to consider just how weak the weak fields can be before the test field approximations for the target mass break down.

Kinematics
Our tractor/pressor/stressor beams will be based on modifications of Natário style generic warp drives [14][15][16][17][18]. The generic form of the space-time metric line element is Note that the lapse is unity, N → 1, the spatial slices are flat, g ij → δ ij , and the "flow" vector v i (t, x, y, z) is the negative of what is (in the ADM decomposition) usually called the "shift" vector [78][79][80][81].
A kinematically useful quantity is the vorticity of the flow field, and its square, ω · ω.
The constant-t spatial slices have covariant normal n a = ∂ a t = (−1, 0, 0, 0) a , whose contravariant components are the future-pointing 4-velocity n a = (1, v i ). Observers that "go with the flow", moving with 4-velocity n a , are geodesics, and are often called Eulerian.
-7 - In the current context this can be recast as [18]: • The Gauss-Mainardi equations yield the Eulerian energy flux: • The 3 × 3 stress tensor is somewhat messier, and can be expressed in terms of the extrinsic curvature and its Lie derivatives [18]: (3.6) For the various explicit examples we consider below, we shall instead often use ab initio calculations instead of this general (but relatively intractable) result.
• In contrast, the trace of the 3 × 3 stress tensor is somewhat easier to deal with.
For the average pressurep we have [18]: ∇ a (Kn a ). (3.8) These are the key stress-energy components we need for the current task. For further discussion on these and related issues see references [18,[78][79][80].
An immediate consequence of these general results is that once appropriate fall-off conditions are imposed at spatial infinity one has This implies that violations of the WEC and NEC are unavoidable [18], and we shall see similar results repeatedly recurring in the subsequent discussion.

Beam profile
In this section, we will discuss the kinematics and general properties of the stress-energy tensor of such beams, including the forces key to our interpretation of them.

Beam kinematics
For our purposes we shall choose a factorized "beam" profile for the flow vector, one that respects axial symmetry around the z-axis: We shall refer to f (x 2 + y 2 ) and h(x 2 + y 2 ) as profile functions, whereas v(t, z) and k(t, z) will be referred to as envelope functions.
Note the explicit presence of x and y in the flow components v x and v y , precisely to maintain axial symmetry. Furthermore, Useful definitions of the average transverse width of the beam are to consider and/or Both of these characterizations of average width depend only on the profile functions, not on the envelope functions.
Far away from the beam axis, as x 2 + y 2 → ∞, we will demand that both profile functions tend to zero: f (x 2 + y 2 ) → 0 and h(x 2 + y 2 ) → 0, in order that the beam asymptotically reduces to flat Minkowski space. All of the t and z dependence is encoded in the two functions v(t, z) and k(t, z). Since one wants the beam to be of finite length, and not stretch all the way across the universe, one should demand both lim z→±∞ v(t, z) → 0 and lim z→±∞ k(t, z) → 0, again ensuring an asymptotic approach to Minkowski space.
More precisely, we shall demand sufficiently rapid fall-off at spatial infinity, which will then also allow integration by parts unrestricted by boundary terms. We shall also enforce smooth on-axis behaviour by demanding that the profile functions and their derivatives be finite on the beam axis. These structural assumptions for the flow vector is basically our definition of what we mean by a "beam" directed along the z-axis.
The previously introduced vorticity (3.2) for our beam geometry reduces to: The square of the vorticity, will show up quite often in subsequent calculations.

Stress-energy basics
If we now additionally impose the factorization conditions (4.1)-(4.2)-(4.3) appropriate to a beam geometry, then the axial symmetry imposes additional constraints on the stress-energy tensor. Specifically: and This implies in particular that and Similarly, for the x-directed and y-directed fluxes, we have: x, y, z) are specific scalar functions that can be explicitly calculated when required. However, the F i (t, x, y, z) are not the most interesting quantities for our purposes. We shall instead be more focussed on the comoving energy density ρ(t, x, y, z), the stress-energy component T zz (t, x, y, z), the flux component f z (t, x, y, z) directed along the beam axis, and the average stressp(t, x, y, z).
We now continue our calculations using the generic beam-like flow (4.1)-(4.2)-(4.3) . As yet, we impose no extra restriction on the four functions v(t, z), k(t, z), f (x 2 + y 2 ), and h(x 2 + y 2 ), apart from the previously mentioned asymptotic conditions. Namely that f (x 2 + y 2 ) → 0 and h(x 2 + y 2 ) → 0 away from the beam axis and both lim

Force
In order to calculate the force (2.1), let us now investigate T zz (t, x, y, z) for this factorized flow, and integrate this over the x-y plane. For T zz (t, x, y, z) we find: Here, using the shorthand u = x 2 + u 2 , we have: Without detailed calculation we can immediately deduce: and Using this, we find that in the narrow beam approximation This is a sum of negative definite and positive definite terms, thus allowing the generic beam to potentially be fine-tuned as either a tractor or a pressor (or even a stressor).
In contrast, in the wide beam approximation we need to evaluate T zz (t, 0, 0, z). Note Consequently (4.30) So in the wide-beam approximation the force exerted on the target is This is of indefinite sign, depending delicately on the envelope functions, potentially allowing either tractor/pressor behaviour.

Flux
The flux in the z-direction, as defined in equation (3.5), is given by: Thence, For the x-direction Similarly, for the y-direction and again by appealing to anti-symmetry, Consequently, for the general tractor/pressor/stressor beam we always have the net flux integrating to zero: Thence, at least in the narrow-beam approximation, we never need to worry about the net fluxes impinging on the target, they always quietly cancel. However, even if the net fluxes seen by Eulerian observers cancel, there might be significant fluctuations around zero over the cross-sectional area of the target. For instance, on axis we have It is now the envelope functions v(t, z) and k(t, z) that primarily drive the localized on-axis fluxes in the wide-beam approximation.

Off-diagonal stress components
Similar steps can be applied to equation (4.14) concerning the T xz and T yz components: implying (using anti-symmetry under x ←→ −x and y ←→ −y respectively) Finally, from equation (4.13) we get: implying (now using either anti-symmetry under x ←→ −x, or anti-symmetry under y ←→ −y) Combining all the above, the integral R 2 Tâb dxdy is purely diagonal, all off-diagonal elements vanish: This really is just a consequence of the assumed axial symmetry of our beam. These observations have the effect of focussing our attention on the diagonal components of the (integrated) stress-energy.

Eulerian energy density
For the Eulerian comoving energy density in this generic beam we find: Then, after an integration by parts, Now we also integrate over z and apply appropriate boundary conditions at z = ±∞ (where the beam has to switch off by definition) to discard the first term, which is a total derivative. Then,

Weak energy condition
This puts us into a good position to have a first look at an energy condition, this time the WEC. Let us do another integration by parts, again invoking suitable boundary conditions, to replace +∞ -15 -But this is now actually a perfect square: (4.55) The integrand appearing above is just 1 4 of the square of the vorticity ( ω · ω), see equation (4.10), so that this is equivalent to This should not come as a surprise, given it is just equation (3.9).
Accordingly, in this generic tractor/pressor/stressor beam configuration, if the Eulerian comoving energy density is positive anywhere, then it must be negative somewhere else -so the WEC is certainly violated.

Null energy condition
Now consider the NEC: Take equation (3.8) and integrate over all space. Note and so (We have already seen in the previous subsection that this last quantity is nonpositive.) Accordingly, in this generic tractor/pressor/stressor beam configuration, if the quantity (ρ +p) is positive anywhere, then it must be negative somewhere else -so the NEC is certainly violated.
-16 -Now, given that the NEC is the weakest of all standard, classical, point-wise energy conditions, we have that all the other energy conditions will also be violated. This has to hold for all tractor/pressor/stressor configurations based on modifications of the generic Natário warp drive. Furthermore, this is completely in accord with what we saw happen for generic warp drive space-times [18].

Special Cases
We now consider three special cases that link our tractor/pressor/stressor discussion back to various previous warp drive analyses [13][14][15][16][17][18]. The connections between the envelope and profile functions of the generic Natário case described by equations (4.1)-(4.3) and those appearing in these special cases is summarised in table 1. Table 1: A summary of the connection between the generic Natário metric, its envelope functions k and v, and its profile functions h and f on the one hand, and the various functions appearing in the special cases considered section 5. .

Modified Alcubierre warp flow
For this particular special case we will assume the field to be oriented along a fixed direction, for convenience taken to be the z direction. This corresponds to taking the flow field to be: For this modified Alcubierre flow field the vorticity is and hence -17 -Now, using the result that for the Alcubierre warp field T zz = 3ρ, obtained in [17], a standard computation yields [13,[16][17][18]: This is already enough to guarantee that both the weak energy condition (WEC) and null energy condition (NEC) are violated in this space-time [17,18]. Calculating the net force we obtain for a narrow beam: But, given our factorization assumption, the stress reduces to Under this assumption the force factorizes to That is, using u = x 2 + y 2 , This is always a tractor beam. The x-y integral is just some positive dimensionless number characterizing the shape of the beam. (Recall that our convention was to always put the target to the right of the generator. If we flip target and generator, so that the target is now on the left and the beam impinges on the target from the right, then there is a sign flip for the force F (t, z), and with F (t, z) > 0 the target is still attracted to the generator.) If we instead assume a wide beam, one can immediately deduce that in this case equation (4.31) will always reduce to zero, as either k or h is zero.

Zero-expansion beam
Now consider a zero-expansion flow field subject to ∂ i v i = 0. Starting with the generic flow field appropriate to an axisymmetric beam, we have: v x (t, x, y, z) = k(t, z) x h(x 2 + y 2 ) (5.11) v y (t, x, y, z) = k(t, z) y h(x 2 + y 2 ) (5.12) v z (t, x, y, z) = v(t, z) f (x 2 + y 2 ). (5.13) -18 -Then, in order to ensure zero expansion, we must enforce Separating variables, one finds for some separation constant C.
Then, without loss of generality, we can enforce: Therefore, the zero-expansion flow field can be rewritten in terms of only two free functions v(t, z) and h( This flow field automatically satisfies axial symmetry, a beam-like profile, and zero expansion. So this is indeed suitable for describing a zero expansion "beam". The vorticity for this beam is easily evaluated as (5.21)

Force
Let us now calculate T zz (t, x, y, z) for this flow field, and then integrate over the x-y plane in order to obtain the net force. For T zz (t, x, y, z) we find: (5.22) -19 -Again using the shorthand u = x 2 + y 2 , we can explicitly calculate: Now consider the integrals over the x-y plane. But first note that Because you want the beam to die off far away from the beam axis, you want h(x 2 +y 2 ) = h(u) → 0 as x 2 + y 2 = u → ∞. So we can already extract some limited information regarding the integrals: Overall, for the zero-expansion narrow beam we now have The first term is indefinite (even though the coefficient is positive), the second term is positive semi-definite, and the third term is negative semi-definite. So the zeroexpansion narrow beam can be tuned to be either a tractor or a pressor, (or even a stressor). One cannot say more about the force F (t, z) without making a specific choice for the profile h(u), and the envelope function v(t, z).
If we now consider a wide beam, then we should look on axis and evaluate T zz (t, 0, 0, z).
-20 -Specifically, we see: (5.35) So in the wide beam limit of a zero-expansion beam we have

Energy conditions
For this zero expansion space-time we have K = tr(K ij ) = 0, and so from equations (3.3) and (3.7) it is immediate that This is enough to guarantee that both the WEC and NEC are violated [18], but for the sake of completeness we perform an explicit calculation.

WEC:
For the Eulerian energy density we note: (5.38) Using the shorthand u = x 2 + y 2 , we can explicitly calculate: This is almost a (negative) perfect square: By performing an integration over the x-y plane, this can then be fully written as the sum of negative perfect squares: This is more than sufficient to guarantee WEC violation somewhere on each x-y plane. Now let us also integrate over dz.
Now, note that So, after an integration by parts, But this now implies that we have a (negative) perfect square: ( ω · ω) dudz ≤ 0, (5.49) as expected. Again, this is more than sufficient to guarantee WEC violation somewhere on each spatial slice, apart from also verifying internal consistency of the formalism.

NEC:
To prove the violation of the NEC, we must now look at the quantity [ρ + T zz ]: Consequently, integrating over the x-y plane we have: This is now a (negative) sum of squares, thereby guaranteeing violation of the NEC. This is again a useful consistency check on the formalism.

Zero-vorticity beam
Let us now consider a zero-vorticity beam described by the flow field: The stress component T zz (t, x, y, z) will again be somewhat complicated. However once one integrates over the x-y plane we shall soon see that That is, there is no net force once you integrate over the entire 2-plane. We shall soon see, however, that there are regions of both repulsion/attraction at various points on the 2-plane. This is best interpreted as a stressor beam.

Force
Explicitly calculating the stress component T zz (t, x, y, z) we find: Now, using again u = x 2 + y 2 for compactness, explicit computation yields Noting again that and observing that each of the P i is a pure derivative, one has So there is no net force.
Note, however, that: which will in general not equal zero. In this way we see that, while the force integrated all over the x−y plane sums up to zero, this does not imply an identically zero force. On the contrary, different parts of the target will be pulled, while others pushed, creating a perfect example of a stressor beam. Furthermore, this means that-ignoring matters of material properties-a wide zero-vorticity beam could potentially still act as a tractor or pressor beam, whereas a narrow zero-vorticity beam would not.

Energy density and null energy condition
Calculating the Eulerian (comoving) energy density we find: Note that R 2 R 0 dxdy = 0, whereas after an integration by parts: Hence, That is, Now, given the fall-off conditions on Φ(t, z), namely that v(t, z) = ∂ z Φ(t, z) z→±∞ −→ 0, we have that: Therefore, if the zero vorticity stressor beam has positive density anywhere, then it must have negative energy density somewhere else. Thence, this zero-vorticity configuration violates the WEC. This is fully in agreement with the general warp-drive analysis presented in [18].
Furthermore, since we have already seen R 2 T zz dx dy = 0, it automatically follows that R 3 T zz dxdydz = 0, and thence we have R 3 (ρ+T zz ) dx dy dz = 0. So, just like before, if the zero vorticity stressor beam has (ρ+T zz ) positive anywhere, then this quantity must be negative somewhere else. Therefore this zero-vorticity configuration also violates the NEC. Again, this is fully in agreement with the general warp-drive analysis presented in [18], and is a useful consistency check on the fact that zero-vorticity flow fields do indeed violate the NEC.
- 25 -In lieu of direct knowledge how one would actually build a tractor beam, one is left with two extremes: General considerations or modelling of specific possibilities. Our discussion in sections 3 and 4 are based on the generic Natário warp drive, in a sense a compromise between the two. While fixing, for example, a certain (3 + 1) split and flat spatial slices in this split, it still retains a large amount of freedom. Section 5 then considered more constrained choices found in the literature, while still retaining some freedom to choose certain functions appearing therein. In this section, we shall illustrate the results of the preceding sections for specific profile functions f and h and specific envelope functions k, v. As a first step, we shall start by imposing a Gaussian profile by fixing the functions f (x 2 + y 2 ) and h(x 2 + y 2 ) to be Gaussian functions. In a second step, we will then employ envelope functions that contribute to the stressenergy tensor between the positions of the generator (at z generator 0) and the target (at z target > z generator ), while vanishing exactly outside of some certain region (−b, b) on the z-axis. More specifically, we will adopt a kind of smooth "bump function". Naturally, these are by far not the only choices, and they contain a certain amount of arbitrariness. Nevertheless, this should give a good idea of what can be done if an arbitrarily advanced civilization could impose stress-energy sources in such a targeted way.

Gaussian beam profiles
As Gaussian beam profiles are very popular toy models in optics and acoustics, they are an obvious starting point for investigating our tractor beams. Let us then provide a few specific examples based on Gaussian beam profiles in the following discussion.

Generic Gaussian beam
Let us consider a generic Gaussian beam, where we set the two profile functions to be identical Gaussians with width parameters a: Then, for the net force exerted by this Gaussian beam (in the narrow field limit), we find a particularly simple factorized form: z) . (6.12) Note that the behaviour switches from pressor to tractor when the beam satisfies the two critical conditions: So, adjusting the two envelope functions is the determining factor in choosing tractor/pressor/stressor behaviour.

Alcubierre-based Gaussian beam
Let us now consider a Gaussian beam based on the modified Alcubierre flow field. Take f (x, y) = exp(−[x 2 + y 2 ]/a 2 ), then from (5.6) and (5.7) we ultimately see Note this Gaussian profile implies that T zz (t, x, y, z) is zero on the z axis, rises to a maximum for (x 2 + y 2 ) ∼ a 2 , and then very rapidly decays as you move further off axis. For the total net force on the x-y plane this Gaussian beam gives: Putting back all the appropriate dimensions, we obtain in SI units Here F * is again the Stoney force. Note that, as expected, this is always a tractor beam.

Zero-expansion Gaussian beam
Looking now at a zero-expansion Gaussian beam, we set h(x 2 +y 2 ) = exp(−(x 2 + y 2 )/a 2 ). Then, using (5.17)-(5.18)-(5.19) and (5.23)-(5.27), we have Thence, for the relevant integrals So, for the Gaussian zero expansion beam, we see that This can be either be a pressor or a tractor beam, depending on the choice of the envelope function. Now consider the wide beam limit. For a Gaussian zero-expansion beam equation (5.35) for T zz reads: In SI units, Again, this can be either be a pressor or a tractor beam, depending on the choice of the envelope function.

Zero vorticity Gaussian beam
If we now take a specific Gaussian profile f (x 2 + y 2 ) = exp(−[x 2 + y 2 ]/a 2 ) then, for a zero vorticity beam we find: We can explicitly check that The sign of the P i (x, y) and consequently the sign of T zz (t, x, y, z) can and will change near x 2 + y 2 ∼ a 2 , so the spatially target will be alternately pushed and pulled -which is why we classify this case as a stressor beam. The calculation on-axis (x = y = 0) gives us: which results in: So, in the wide-beam limit, As we can see, this is another "tunable" case, which can behave either as a pressor or a tractor beam, depending on the choice of the envelope function.

Envelope functions
In order to be able to visualize some of the properties of tractor/pressor/stressor beams we shall now impose two different possibilities for the envelope functions v(t, z) and k(t, z). This will allow us to plot the force field generated by these functions and the energy density distribution necessary to create them. All of the calculations done in the previous sections are still completely valid here.

Illustrating Gaussian beams
In figure 1, used in the Introduction to describe where target and generator are with respect to the tractor field, we also plotted the energy densities and forces of the nontrivial beam configurations described above. To produce those plots, we imposed a Gaussian envelope together with a Gaussian profile for the defining functions: where, for the plotting, we used A = 0.5, B = C = 1.0, and we evaluated the energy density and forces at t = 1. Note that F (t, z) for both the narrow Alcubierre and the wide zero-vorticity beams are always negative for this specific setup, implying a tractor beam behaviour, while the other beam configurations allow for a tractor/pressor behaviour depending on the positioning of the target. It is also nice to notice how nontrivial is the cancelation of the energy density along the spatial 3-slices for the zero vorticity case, given by equation (5.71) and represented in figure 1-(b).

Bump functions
Another, much more brutal way of enforcing that fall-off conditions be fulfilled is by Use this to then define .
In a last step, define for real numbers a and b As we are interested in functions satisfying appropriate fall-off conditions at infinity, this example fulfils this by construction in the most trivial way possible: It vanishes for sufficiently large positive or negative values of x. Furthermore, as we are specifying the metric by hand, the Einstein equation will tell us the required sources; just as in all the calculations of this paper. Neither the Gaussian beams nor beams based on such smooth bump functions differ in this regard from each other, and the general analysis of the previous sections will still hold. Nevertheless, using such smooth bump functions for the envelope functions v or k is an intriguing way to model a tractor beam that only contributes to the stress energy on the z-axis between "generating device" and "target".
The algebra becomes arbitrarily involved in this case; for this reason we opt to only show our results and the functions we chose. The bump function used is which, depending on the specific (special) case plotted, was used for v, k, or Φ. The profile functions were again chosen to be the Gaussians, as described in section 6.1.1, which also allows an easier comparison with the plots shown in figure 1. In figure 3, the parameters are t = −1, a = 2, b = 10, and D = 1. In figure 4, the parameters are t = −1, a = 2, b = 4, and D = 1. Just this minor variation produced noticeable changes in the forces and energy density. The choice of t can also produce significant differences, but this is not shown here, as it adds little to the discussion. Again, note how non-trivial is the distribution of the energy density for the zero vorticity beams, which sums up to zero when integrated over any 3-spatial slice. It is also interesting to notice the different behaviour for distinct types of beam, varying from constant pull forces (e.g. the Alcubierre case) up to elaborate push/pull behaviours (e.g. wide zero-expansion). This reveals the great diversity of mechanisms one can create by varying the envelope functions only. Behaviour for different types of profile functions might possibly create yet other interesting scenarios, which we will leave for the enthusiastic reader.
-33 - Setting aside the issue of the magnitude of the Stoney force (which can be taken care of by an appropriately small pre-factor in our functions), we in particular like to draw attention to the force of the zero expansion beam in figure 3a: A target positioned to the right at z ≈ 10 would be accelerated to the left, then travel for a while at near constant velocity, before it is decelerated again. Sufficient fine-tuning thus allows for safe docking or boarding.
- 34 -In this article we have seen how to analyze tractor/pressor/stressor beams within the framework of standard general relativity. The analysis was made based on modified warp drive spacetimes, by creating a "beam like" profile. A general case based on Natário's warp field was analyzed, followed by specific cases and examples.
As expected, we have seen that in this case, just like with warp drives and traversable wormholes, the violation of the NEC, and so of all the classical point-wise energy conditions, is unavoidable. A closely related statement is still true even if one moves beyond Einstein gravity. The key point is that it is ultimately the focussing properties of the tractor/pressor/stressor beams, warp fields and traversable wormholes that translate into convergence conditions [99][100][101][102][103][104][105], and thence into [effective] energy conditions. Whenever you can rearrange the equations of motion in the form then the effective energy-momentum tensor [T effective ] ab will consequently violate the NEC and so violate all the classical point-wise energy conditions.
However, a significant question remains open: Are energy conditions truly fundamental physics? Probably not, (indeed, almost certainly not). But the energy conditions are certainly good diagnostics for unusual physics -and, as we have seen, the physics of these tractor/pressor/stressor beams is certainly extremely unusual -comparable in weirdness to that of traversable wormholes and warp drives. This is not an absolute prohibition on tractor/pressor/stressor beams, but it is an invitation to think very carefully about the underlying physics.