On Gravitational Chirality as the Genesis of Astrophysical Jets

It has been suggested that single and double jets observed emanating from certain astrophysical objects may have a purely gravitational origin. We discuss new classes of plane-fronted and pulsed gravitational wave solutions to the equation for perturbations of Ricci-flat spacetimes around Minkowski metrics, as models for the genesis of such phenomena. These solutions are classified in terms of their chirality and generate a family of non-stationary spacetime metrics. Particular members of these families are used as backgrounds in analysing time-like solutions to the geodesic equation for test particles. They are found numerically to exhibit both single and double jet-like features with dimensionless aspect ratios suggesting that it may be profitable to include such backgrounds in simulations of astrophysical jet dynamics from rotating accretion discs involving electromagnetic fields.


Introduction
Many astrophysical phenomena find an adequate explanation in the context of Newtonian gravitation and Einstein's description of gravitation is routinely used (together with Maxwell's theory of electromagnetism and the use of time-like spacetime geodesics to model the histories of massive point test particles) to analyse a vast range of phenomena where non-Newtonian effects are manifest. However, there remain a number of intriguing astrophysical phenomena suggesting that our current understanding is incomplete. These include the large scale dynamics of the observed Universe and a detailed dynamics of certain compact stellar objects interacting with their environment. In this note we address the question of the dynamical origin of the extensive "cosmic jets" that have been observed emanating from a number of compact rotating sources. Such jets often contain radiating plasmas and are apparently the result of matter accreting on such sources in the presence of magnetic fields. One of the earliest models to explain these processes suggested that the gravitational fields of rotating black holes surrounded by a magnetised "accretion disc" could provide a viable mechanism [1]. More recently, the significance of magneto-hydrodynamic processes in transferring angular momentum and energy into collimated jet structures has been recognised [2,3,4]. Many of these models implicitly assume the existence of a magnetosphere in a stationary gravitational field and employ "force-free electrodynamics" in their development. To our knowledge, a dynamical model that fully accounts for all the observed aspects of astrophysical jets does not exist.
However in recent years there has been mounting evidence, both theoretical and numerical, suggesting that the genesis of jets may have a purely gravitational origin. By the genesis of such phenomena we mean a mechanism that initiates the plasma collimation process whereby electrically charged matter arises from initial distributions of neutral matter in a background gravitational field. In [5,6], the authors carefully analyse the properties of a class of Ricci-flat cylindrically symmetric spacetimes possessing time-like and null geodesics that approach attractors confining massive particles to cylindrical spacetime structures. Additional studies [7,8,9] of the asymptotic behaviour of test particles on time-like geodesics with large Newtonian speeds relative to a class of co-moving observers have given rise to the notion of cosmic jets associated with different types of gravitational collapse scenarios satisfying certain Einstein-Maxwell field systems. There has also been a recent approach based on certain approximations within a linearised gravitational framework involving "gravito-magnetic fields" generated by non-relativistic matter currents [10]. All these investigations auger well for the construction of models for astrophysical jets that include non-Newtonian gravitational fields as well as electromagnetically induced plasma interactions.
Although astrophysical jets involve both gravitational and electromagnetic interactions with matter it is natural to explore the structure of electrically neutral test particle geodesics in non-stationary, anisotropic background metric spacetimes as a first approximation to what is undoubtedly a complex dynamical process. In this paper we explore geodesics that relax any assumption of a cylindrically Killing-symmetric background spacetime metric. Furthermore, to this end we construct particular exact solutions to the linearised Einstein vacuum equations which are then used to numerically calculate time-like geodesics in non-stationary backgrounds. The use of the linearised Einstein vacuum equations facilitates the construction of families of complex eigensolutions with definite chirality that are used to construct real spacetime metrics exhibiting families of timelike geodesics possessing particular jet-like characteristics on space-like hyper-surfaces. Test particles on such time-like geodesics exhibit, in general, a well defined sense of "handed-ness" in space that we argue may offer a mechanism that initiates a flow of matter into directed jet-like structures. In particular, we construct families of plane-fronted gravitational wave metrics and new non-stationary metrics having propagating pulselike characteristics with bounded components in three-dimensional spatial domains. These are analogous to particular exact solutions of the vacuum Maxwell equations which we have recently shown can be used to model single-cycle electromagnetic laser pulses [11,12].
In section 1, our notation is established in the context of Einstein's vacuum equations on a spacetime manifold M ( g) with a Lorentzian metric g and their linearisation about a flat Minkowski metric η on a region U ⊂ M ( g) . Section 2 describes a family of complex plane-fronted gravitational wave exact solutions to the linearised field equations that are eigen-tensors of the complex axial-symmetry operator about any particular direction in space. We call the eigenvalues of this operator "chirality" and show how the associated "helicity" of the real-part solutions with respect to their direction of propagation is correlated with certain properties of time-like geodesics associated with the real linearised spacetime metric. The simplest harmonic plane-fronted gravitational wave modes have helicity ±2 and, in general, induce a (modulated) non-circular helicoidal motion of massive testparticles that are initially arranged (with non-zero speeds) in a ring that lies in a plane orthogonal to the direction of propagation of such waves. The set composed of the individual particle paths form a uni-directional bi-jet-like array in space emanating from such a ring in a background harmonic wave of definite helicity.
By contrast in section 3, we give a construction of exact solutions to the linearised equations with pulse-like perturbations that are bounded in all spatial dimensions. The simplest gravitational pulse-like solutions have chirality zero and are shown to generate a real non-singular Ricci curvature scalar field on U with well defined loci in spacetime emanating from a core containing a maximum at a particular event. Future-pointing timelike geodesics that emanate from events on rings containing massive test particles centred on this core can give rise to single and oppositely directed jet-like arrays in space, transverse to the plane of the rings (bipolar outflows).
It is assumed that exact analytic solutions to the linearised Einstein vacuum equations and their associated time-like geodesics will exhibit features that persist to some degree beyond the linearisation regime and, in particular, offer an approach to a better understanding of the genesis of observable cosmic jets in models that include charged matter with plasma interactions. Based on a numerical exploration of particular time-like geodesics associated with background metrics constructed from complex eigen-solutions with definite chirality, we conclude that it may be profitable to include these non-Newtonian gravitational backgrounds in simulations of cosmic jet dynamics from rotating accretion discs involving electromagnetic fields.

The Linearised Einstein Vacuum System
In any matter free domain of spacetime U ⊂ M ( g) , an Einsteinian gravitational field is described by a symmetric covariant rank-two tensor g with Lorentzian signature that satisfies the vacuum Einstein equation Ein ( g) = 0 where and Ric ( g) is the Levi-Civita Ricci tensor associated with the torsion-free, metric-compatible Levi-Civita connection ∇ ( g) . A coordinate independent linearisation of (1) can be found in [13,14]. In particular a linearisation about a flat Minkowski spacetime metric η on U determines the linearised metric 1 g = η + h and to first order one writes g = g + O(κ 2 ). The variable κ is a parameter in h used to keep track of the expansion order and in any η−orthonormal coframe 2 on U. Since we explore the source-free Einstein equation (relevant to the motion of test-matter far from any sources) the physical scale associated with any linearised solutions must be fixed by the solutions themselves rather than any coupling to self-gravitating matter. Furthermore since only dimensionless relative scales have any significance we define the tensor h to be a perturbation of η on U relative to any local η−orthonormal coframe {e a } provided where e a ∈ T * U, X b ∈ T U, e a (X b ) = δ a b . It should be noted that the perturbation order of the η−covariant derivative of a tensor and its η−trace relative to such a coframe is not necessarily of the same order as that assigned to the tensor. Thus perturbation order is not synonymous with "scale" in this context. We use the conditions (3) to define perturbative Lorentzian spacetime domains P U ⊂ U to be regions where max 0≤a,b≤3 |h(X a , X b )| < 1.
The real tensor h ≡ µ (η) ψ with ψ ≡ Re(ψ) may be constructed from any complex covariant symmetric rank two tensor ψ satisfying [14]: Here and below, ∇ (η) denotes the operator of Levi-Civita covariant differentiation associated with η, Xa X a and for all covariant symmetric rank two tensors T on U: 1 Throughout this paper we assign physical dimensions of length 2 to the tensors g and η. The Ricci curvature scalar associated with g then has the dimensions of length −2 . In a g-orthonormal coframe the components of g are {−1, 1, 1, 1} and in an η−orthonormal coframe the components of η are {−1, 1, 1, 1}. This does not imply that components of the tensor field η are necessarily constant in an arbitrary coframe on U .
2 An arbitrary local coframe is a set of 1−forms {e a } on U satisfying e 0 ∧ e 1 ∧ e 2 ∧ e 3 = 0. If ψ = ψ ab e a ⊗ e b in any such coframe, Tr (η) ψ = ψ ab η ab with η ab η bc = δ a c in terms of the Kronecker symbols δ a c .
If any worldline C has components C µ (τ ) in any local chart on U with coordinates {x µ } andĊ µ ≡ ∂ τ C µ then where Γ µ αβ denotes a Christoffel symbol associated with ∇ (g) .
In the following sections only solutions to (5) that lie in the perturbative domains P U are displayed. The worldline of an idealised observer in U is modelled by the integral curve C V of a future-pointing time-like unit vector field V , (i.e. g(V, V ) = −1). At any event in U the g−orthogonal decomposition ofĊ with respect to an observer C V : Relative to C V , the observed "Newtonian speed" of the proper-time parameterised worldline C at any event V V = 0 the observer is said to be geodesic otherwise it will be accelerating. If there exists a local coordinate system {t, ξ 1 , ξ 2 , ξ 3 } on U with ∂ t time-like, in which C V can be parameterised monotonically with λ as t = λ, ξ 1 = ξ 1 (0), ξ 2 = ξ 2 (0), ξ 3 = ξ 3 (0) such an observer is said to be at rest in this coordinate system. Although any particular time-like worldline defines a local "rest observer" in some chart, only the existence of a family of rest observers in a particular chart Φ U on U provides a way to interpret the Newtonian velocity of any event on a time-like worldline that is not necessarily a rest-observer in U. In units 3 with c 0 = 1, a point particle of rest-mass m 0 with a worldline C, when observed by C V , has energy and 3-momentum values at any event on C given by

Gravitational Waves
The recent experimental detection of gravitational waves lends credence to the existence of non-stationary metrics on spacetime domains in the vicinity of colliding black holes. There are good reasons to believe that such domains where c 0 is a fundamental constant with the physical dimensions of speed. In SI units c 0 = 3 × 10 8 m/sec. The dimensionless amplitudes in (8) define amplitudes L 0 q 1 and L 0 q 2 with physical dimensions of length. In this scheme with g(Ċ,Ċ) = −1, the parameter τ has length dimensions and its conversion to a parameter τ with dimensions of a clock time is given by τ = τ /c 0 = τ L 0 /c 0 where τ is a dimensionless parameter. may exist in the vicinity of other non-stationary astrophysical processes. In this section we explore the properties of particular proper time parameterised plane-fronted gravitational wave spacetimes in the perturbative domains P U defined above.
In a local chart Φ U possessing dimensionless coordinates 3 {T, R, θ, Z} with T ≥ 0, R > 0, θ ∈ [0, 2π) and |Z| ≥ 0 on a spacetime domain U ⊂ M ( g) , a local coframe C adapted to these co-ordinates is {e 0 = dT, e 1 = dR, e 2 = R dθ, e 3 = dZ} in terms of the exterior derivative operator d. With η given by (2), this coframe is η−orthonormal but not in general g−orthonormal. Such a chart facilitates the coordination of a series of massive test particles initially arranged in a series of concentric rings with different values of R lying in spatial planes with different values of Z at T = 0. Furthermore, we define on U: For σ = ±1, m 1 , m 2 ∈ Z, particular non-harmonic complex plane-fronted gravitational wave solutions to (4) take the form where W is an arbitrary complex-valued function that is at least twice differentiable on U. Waves with opposite values of σ propagate in opposite directions along the axis Z. These tensors satisfy 4 and will be called complex eigen-tensors of definite chirality m 1 + m 2 . Such tensors give rise to real tensors (ψ σ m1+m2 ) ≡ Re(ψ σ m1+m2 ) with "helicity" (m 1 + m 2 )σ. If |m 1 | ≥ 1, |m 2 | ≥ 1 and m 1 m 2 > 0, it may be shown that Lap (η) ψ σ m1+m2 = 0, Div (η) ψ σ m1+m2 = 0 and Tr (η) ψ σ m1+m2 = 0. Since Lap (η) and Div (η) are linear operators, the superposition σ,n C σ n ψ σ n with arbitrary complex constant coefficients C σ n is also a solution of (4). From the simplest trace-free solution with m 1 = 1 and m 2 = 1 one may construct from (ψ σ 2 ) the real helicity 2σ symmetric covariant metric tensor g σ 2 on U: Pure harmonic plane-fronted gravitational waves arise when the function W in (7) is a bounded trigonometric function of its arguments. The higher complex chirality solutions are related to ψ σ 2 by the formula: The real tensors (ψ σ 2 ) remain trace-free, divergence-free solutions to Lap (η) (ψ σ 2 ) = 0 under the finite real rotation map θ → θ + θ 0 with arbitrary real constant θ 0 .
We have explored numerically the structure of time-like geodesics on U associated with g σ 2 for pure harmonic plane-fronted gravitational waves where: with dimensionless angular frequency Ω and real dimensionless amplitudes q 1 , q 2 . In these gravitational wave spacetimes the geodesic vector field ∂ T ∈ T U satisfies g σ 2 (∂ T , ∂ T ) = −1. In U, the metric tensor g 1 2 has non-zero 4 For any vector field ξ, the operator L ξ denotes Lie differentiation with respect to ξ.
components in the coframe C: The integral curves of the vector field ∂ T provide, in Φ U , a family of geodesic rest observers on U. In this chart we write the proper-time parameterised geodesic curves as Only three of the equations in (6) are independent. In the L 0 T = τ gauge the equations with µ = 1, 2, 3 are . These initial conditions must be chosen consistently with the condition g σ 2 (Ċ,Ċ)| τ =0 = −1 which is then preserved for all τ . Furthermore the parameters q 1 , q 2 , Ω are chosen so that throughout the numerical integration of (5) all components of (ψ σ 2 ) in the η−orthogonal coframe evaluated on C for τ ≥ 0 remain less than unity in absolute value. This ensures that all computed geodesics lie in the perturbative domains P U . The solutions to (5) wheṅ test particles on these geodesics remain at rest relative to rest observers. Solutions R(τ ), θ(τ ), Z(τ ) are displayed individually or as single space curves in Cartesian domains labelled with dimensionless coordinates X, Y, Z where X = R cos(θ) and Y = R sin(θ). For each solution in P U one can calculate the square of the Newtonian speed v at any τ relative to a rest observer: To display the general features of the geodesic solutions to (4) that exhibit jet-like features for specified initial conditions we select particular values of the parameters Ω, q 1 , q 2 , L 0 and helicity 2σ where σ = ±1. Having the notion of a thick accretion disc in mind, we choose initial (τ = 0) positions of individual geodesics to be equally spaced in θ on rings with fixed R in planes with fixed Z. On the left in figure 1 a set of 8 such time-like geodesics is displayed in an X-Y projection, each for an evolution proper time from τ = 0 to τ = 4 × 10 4 . Each geodesic exhibits a segment of a helicoidal structure in a background gravitational wave of definite helicity 2. The sense of the helicoidal motion is reversed for waves of opposite helicity as demonstrated in the figure on the right. Although the average motion of the helicoids along the Z−axis is determined by the direction of the wave, the motions R(τ ) and Z(τ ) of a test particle with an initial non-zero radial speed are not immediately monotonic in τ , as evident in figure 2. The corresponding parametric behaviours of the solutions R(τ ), θ(τ ), Z(τ ) are displayed in figure 3.
The left hand side plot in figure 4 displays a short segment of a typical helicoid on a scale that makes visible a cycloidal-like oscillatory modulation. The spatial periodicity of this modulation is determined by the value of Ω and decreases as this parameter increases. The right hand side of figure 4 displays on a lower resolution (thereby suppressing the cycloidal modulation) a segment of a typical helicoid that exhibits helicoidal drift and a decrease of helicoidal radius as τ increases. Relative to a rest observer, the dimensionless Newtonian speed v(τ ) of a test particle on each helicoid is oscillatory but in general tends to unity for large values of τ .  displays the evolution of 8 geodesics initially arranged on two different rings, each of radius 0.06. Four of them emanate from a ring lying in the plane Z = 0.01 with initial θ values of 0, π/2, π, 3π/ and four of them from a ring with the same radius and same θ values lying in the plane Z = −0.01. The gravitational wave with helicity 2 is propagating in the negative Z−direction. The figure on the right has τ for all geodesics varying from zero to 2 × 10 3 while that on the right shows the evolution from zero to 5 × 10 3 . Both figures capture the cycloidal-like modulations and the general uni-directional emanations induced by geodesics with definite helicity but different angular locations on identical radius rings above and below the plane Z = 0. It is also evident that geodesics emanating from the values θ = 0 and θ = π initially execute different motions from those starting at θ = π/2  and θ = 3π/2. Figure 6 displays the evolution of 16 geodesics initially arranged on two different rings, each with radius 1. Eight of them emanate from a ring lying in the plane Z = 0.245 with initial θ values of nπ/4 for n = 0, .., 7 and eight of them from a ring with the same radius and same θ values lying in the plane Z = −0.245. The gravitational wave with helicity 2 is propagating in the negative Z−direction. The proper time τ for all geodesics varies from zero to 5.8 × 10 5 . With this graphical resolution no cycloidal-like modulations are visible. It is evident that uni-directional emanations from different initial angular locations are being induced by harmonic waves with definite helicity irrespective of their origins on rings above or below the plane Z = 0. It also suggests that for a fixed evolution time a relativistic uni-directional bi-jet structure is dominant. The detailed properties of time-like geodesics in the perturbative domains of these gravitational wave spacetimes depend on the values assigned to a number of their parameters but the emergent general features have been displayed in figures (1)- (6). It appears that, in background definite helicity harmonic waves, all time-like geodesics with non-zero initial Newtonian speeds relative to rest-observers evolve into uni-directional modulated helicoids. Being plane-fronted these characteristics are independent of the initial values of R occupied by test matter.
Although plane-fronted gravitational waves (like plane-fronted electromagnetic waves) are physical idealisations they offer a possible mechanism for the initiation of uni-directional jet-like patterns when incident on local distributions of matter. In the next section we explore a class of non-stationary chiral solutions to the field equations (4) that are neither plane-fronted nor harmonic gravitational waves.

Compact Gravitational Pulses
In [11,12], following pioneering work by Synge [15], Brittingham [16], Stewart [13] and Ziolkowski [17], we have shown how to construct complex analytic solutions to the Maxwell vacuum equations describing propagating pulse-like electromagnetic fields. These have found direct application for analysing the behaviour of singlecycle laser pulses. One notable feature of such solutions is that, although they do not have a unique direction of propagation (and hence helicity), they are derived from solutions with definite chirality. To construct real spacetime metrics analytically with analogous pulse-like structures, we now outline how the techniques used to construct antisymmetric Maxwell tensors F satisfying 5 dF = 0, δF = 0 can be generalised to construct symmetric second rank covariant tensors satisfying (4). Unlike the gravitational wave solutions above, the metric perturbations are compact in all spatial directions and it will be shown describe non-stationary spacetimes possessing time-like geodesics with jet-like families of space-curves.
If α is any complex (four times differentiable) scalar field on U then Hence if Lap (η) α = 0, then Lap (η) ψ 0 = 0 and Div (η) ψ 0 = 0 are satisfied with 6 ψ 0 = ∇ (η) dα. Furthermore since ∇ (η) is the torsion-free Levi-Civita connection, ψ 0 is symmetric and trace-free with respect to η. Hence, such complex α in general give rise to real metrics: The linear divergence operator δ acting on p−forms in spacetime is defined by d in terms of of the Hodge map associated with a metric g. 6 For any scalar α and metric tensor g, ∇ (g) α = dα is independent of g. with curvature. The particular complex solution α of relevance here is given in the (T, R, θ, Z) chart Φ U above as with κ, Q 1 , Q 2 strictly positive real dimensionless constants. The scalar α(T, R, Z) is then singularity-free in T , R and Z and clearly axially-symmetric with respect to rotations about the Z-axis. It also gives rise to an axially-symmetric complex (zero chirality) tensor ψ 0 satisfying 7 L ∂ θ ψ 0 = 0. In U, the real axially-symmetric metric tensor g 0 has non-zero components in the coframe C: where α ≡ Re α and satisfies L ∂ θ g 0 = 0. A rest observer field in this spacetime metric is 1/ 1 − ∂ 2 T T ψ 0 ∂ T ∈ 7 Since ∇ (η) is a flat connection, if K is an η−Killing vector then the operator ∇ (η) L K = L K ∇ (η) on all tensors.

T U.
Complex symmetric tensors ψ m with integer chirality m > 0 satisfying Lap (η) ψ m = 0, Div (η) ψ m = 0, Tr (η) ψ m = 0 and 1 i L ∂ θ ψ m = mψ m may be generated from ψ 0 by repeated covariant differentiation with respect to a particular η−null and η−Killing complex vector field 7 S: Solutions with negative integer chirality can be obtained by complex conjugation of the positive chirality complex eigen-solutions. Each ψ m defines a spacetime metric g m = η + ψ m = η + Re ψ m on U ∈ M (gm) which, for m = 0, is not axially symmetric: L ∂ θ g m = 0. An indication of the nature of the spacetime geometry determined by g m on M (gm) is given by analysing the structure of the associated Ricci curvature scalar R gm (T, R, Z). Unlike in the gravitational wave spacetimes, this scalar is not identically zero. For m = 0 it is axially symmetric and in the chart Φ U its independence of θ means that for values of a fixed radius its structure can be displayed for a range of T and Z values given a choice of parameters (Q 1 , Q 2 , κ, L 0 ). Regions where R g0 (T, 1, Z) change sign are clearly visible in the right side of figure 8 where a 2-dimensional density plot shows a pair of prominent loci that separately approach the future (T ≥ 0) light-cone of the event at {R = 1, T = 0, Z = 0}. A more detailed graphical description of R g0 (T, 1, Z) is given in the left hand 3-dimensional plot in figure 8 where an initial pulse-like maximum around T 0 evolves into a pair of enhanced loci with localised peaks at values of Z with opposite signs when T ≥ 1. In this presentation the maximum pulse height has been normalised to unity. This characteristic behaviour is similar to that possessed by α (T, 1, Z). It suggests that "tidal forces" (responsible for the geodesic deviation of neighbouring geodesics [18,19,20]) are concentrated in spacetime regions where components of the Riemann tensor of g 0 have pulse-like behaviour in domains similar to those possessed by R g0 (T, R, Z).
Explicit formulae for R g0 (T, R, Z) and H(T, R, θ, Z) are not particularly illuminating 8 . However, for fixed values of the parameters (Q 1 , Q 2 , κ, L 0 ), their values can be plotted numerically in order to gain some insight into their relative magnitudes in any perturbative domain P U . With Z fixed at zero, figure 7 displays such plots as functions of R and a set of T values. It is clear that in perturbative domains the curvature scalar may exceed unity. Since in general: where Q is a non-singular rational function of its arguments and the tensor h is, by definition, of order κ, figure 7 demonstrates that relative tensor κ−orders are not, in general, indicators of their corresponding relative magnitudes.
The time-like geodesics of g 0 satisfying (5) on P U are found numerically following the procedures given in the previous section. To emulate matter on a thick accretion disc test particles are placed initially on rings with fixed R lying in planes with fixed Z with specified initial conditions in perturbative domains. Since g 0 is axially symmetric with respect to Z, all geodesics will inherit this symmetry. The figures below exhibit time-like geodesics starting at events with L 0 T = 0, R = R(0), Z = Z(0) and evolving to L 0 T = τ max , R = R(τ max ), Z = Z(τ max ) for various values of θ(0). We assign to a collection of evolved time-like geodesics the value of the dimensionless aspect ratio A(τ max ) defined by: On the left of figure 9 six geodesics associated with pulse parameters Q 1 = Q 2 = 1 are shown emanating from six locations with θ values 0, π/3, 2π/3, π, 4π/3, 5π/3 on a ring with radius 10 −4 in the plane Z = 0.735 and six from similarly arranged points on rings of the same radius at Z = 0.245, Z = −0.245 and Z = −0.735. The initial locations are not resolved in these figures. The 24 geodesics each evolve from τ = 0 to τ = 10 4 and clearly display an axially symmetric bi-directional jet structure from the rings in conformity with the expectations based on the spacetime structure of R g0 (T, 1, Z) in figure 8. The figure on the right resolves the structure of the directional bi-jet array for 0 ≤ τ ≤ 100. A single uni-directional jet-array arises when only one ring is populated with matter. This figure demonstrates that the jets from the sources at Z = ±0.245 have a dimensionless aspect ratio A(10 4 ) = 64.7 much greater than those produced from the sources at Z = ±0.735 where A(10 4 ) = 3.22. The equality of the aspect ratios for a pair of jets produced from sources placed symmetrically around Z = 0 is a result of the choice Q 1 = Q 2 . When Q 1 = Q 2 , it is possible to produce, from a pair of sources on rings arranged symmetrically with respect to Z = 0, a pair of uni-directional jets with members of each jet having different aspect ratios. This is demonstrated in figure 10 where the jet sources on the left have larger separations in Z than those in the centre, all other initial conditions and pulse parameters being the same. If all such sources are in evidence, one produces two pairs of oppositely directed jets, as shown in the right-hand figure, but having different aspect ratios:

Summary and Concluding Remarks
In this article, we have exploited a linearisation of the Einstein vacuum equations about a Minkowski spacetime to construct families of spacetime metrics describing non-stationary perturbed geometries. These metrics are constructed from complex, symmetric, covariant tensors ψ with definite chirality satisfying the tensor equations: Lap (η) ψ = 0, Div (η) ψ = 0 and Tr (η) ψ = 0.
The family of solutions in section 2 describe plane-fronted gravitational waves of definite helicity while those in section 3 are gravitational pulses of definite chirality. We have explored the nature of the time-like geodesics in perturbative spacetime domains associated with the lowest allowed chirality solutions in each family.
Using suitably arranged massive test particles to emulate a thick accretion disc, together with a particular family of fiducial observers, we have displayed a number of characteristic features of these geodesics in each On the left, six geodesics are shown emanating from six locations with θ values 0, π/3, 2π/3, π, 4π/3, 5π/3 on a ring with radius 10 −4 in the plane Z = 0.735 and six from similarly arranged points on rings of the same radius at Z = −0.735. The initial locations are not resolved in these figures. The 12 geodesics each evolve from τ = 0 to τ = 10 4 and clearly display an axially symmetric uni-directional jet structure from the rings. All geodesics are generated with the additional initial conditionsṘ(0) = 0,Ż(0) = 0, θ(0) = 0.4 and the background perturbation pulse has parameters L 0 = 1, Q 1 = 1, Q 2 = 3, κ = 1/6. The figure in the centre shows an oppositely directed jet evolving from similar initial conditions, but with initial Z = 0.245 and Z = −0.245. The figure on the right displays the "asymmetric" jet structure obtained by merging both pairs of sources with dominant component belonging to the jet in the left-hand figure having aspect ratio A(10 4 ) = 178.8. background metric. Within the context of a non-dimensional scheme, solution parameters can be chosen that result in characteristic spatial jet-like patterns in three-dimensions. These have specific dimensionless aspect ratios relative to well-defined directions of a background harmonic gravitational wave or background gravitational pulse and the corresponding orthogonal subspace. For a helicity ±2 plane-fronted harmonic wave incident at T = 0 on a bounded region of matter in the vicinity of the spatial plane Z = 0 in three-dimensions, parameters (q 1 = q 2 ) can be chosen to yield uni-directional bi-jets: i.e. a pair of time-like geodesic families with Newtonian speeds approaching the speed of light in the Minkowski vacuum, relative to geodesic observers, for large proper times. At such times, each pair lies in the domain Z > 0 (Z < 0) if the wave propagates in the direction of increasing (decreasing) Z. For the chirality zero gravitational pulse incident at T = 0 on a similar disposition of matter, one finds that (with Q 1 = Q 2 ) a pair of oppositely directed jet-like structures arise: i.e. a pair of time-like geodesic families with Newtonian speeds approaching terminal values less than the speed of light for both Z > 0 and Z < 0. For a pulse with Q 1 = Q 2 , we have demonstrated the existence of a pair of uni-directional jet-like structures from particular initial conditions. In all these cases, the structures have well-defined aspect ratios that can be calculated numerically. The propagation characteristics for T > 0 of the pulse responsible for these jet structures in space is discernible from features of the non-zero Ricci scalar curvature associated with the perturbed spacetime domains.
We have also stressed that by linearising only the gravitational field equations and analysing the exact geodesic equations of motion in perturbative spacetime domains, one can capture the full effects of "tidal accelerations" on matter produced by the curvature tensor (and its contractions) associated with the metric perturbations. This opens up the possibility of a gravito-ionisation process whereby extended electrically neutral micro-matter can be split into electrically charged components by purely gravitational forces, leading to modifications of matter worldlines by the presence of Lorentz forces.
We conclude that both background spacetime families discussed, separately or in superposition with higher chirality solutions, may offer a non-Newtonian gravitational mechanism for the initialisation of a dynamic process leading to astrophysical jet structures, particularly since it is unlikely that such a phenomenon originates from a unique set of initial conditions. The nature of possible physical sources for the gravitational backgrounds discussed in this article and the influence of electromagnetic interactions needed to set our results into a particular astrophysical context will be discussed elsewhere.