Kundt geometries and memory effects in the Brans-Dicke theory of gravity

Memory effects are studied in the simplest scalar-tensor theory, the Brans--Dicke (BD) theory. To this end, we introduce, in BD theory, novel Kundt spacetimes (without and with gyratonic terms), which serve as backgrounds for the ensuing analysis on memory. The BD parameter $\omega$ and the scalar field ($\phi$) profile, expectedly, distinguishes between different solutions. Choosing specific localised forms for the free metric functions $H'(u)$ (related to the wave profile) and $J(u)$ (the gyraton) we obtain displacement memory effects using both geodesics and geodesic deviation. An interesting and easy-to-understand exactly solvable case arises when $\omega=-2$ (with $J(u)$ absent) which we discuss in detail. For other $\omega$ (in the presence of $J$ or without), numerically obtained geodesics lead to results on displacement memory which appear to match qualitatively with those found from a deviation analysis. Thus, the issue of how memory effects in BD theory may arise and also differ from their GR counterparts, is now partially addressed, at least theoretically, within the context of this new class of Kundt geometries.


I. INTRODUCTION
The detection of gravitational waves in binary mergers has opened up new prospects for testing theories of gravity in the strong field regime [1,2]. Gravitational wave memory is one such as-yet-unobserved strong field effect that can be used to test diverse theories of gravity. The gravitational wave memory effect is the residual permanent shift in the position (or velocity) caused due to the passage of a gravitational wave pulse [3].
The study of memory effects began in the work of Zel'dovich and Polnarev [4] who studied gravitational radiation emitted due to the motion of flybys/collapse of stars in a globular cluster. A few years later, Braginsky and Grishchuk [5] looked at the deviation of test particles in weak field, linearized gravity and coined the term memory effect to denote the change in the metric perturbation at early and late times. Geodesic deviation of test particles due to low frequency gravitational radiation at null infinity was investigated further by Ludvigsen [6]. Christodoulou, using full nonlinear GR, attributed the presence of memory to the transport of energy and momentum of gravitational waves to null infinity [7]. Further, this effect, related to non-linearity, was ascribed to gravitons produced by the radiation itself [8]. Memory effects are also possible in electrodynamics [9] and Yang-Mills theories [10,11].
Interesting theoretical links to memory effects have been conjectured, of late, in the context of soft theorems and BMS symmetries [12]. It has been noted that the nonlinear memory effect can also be understood as a BMS transformation relating two inequivalent Minkowski vacua at future null infinity caused by the passage of gravitational waves (see the review [13] and the references cited therein).
Apart from spacetime boundaries, memory effects can also be realized in the interior regions of a spacetime [14,15]. Such effects lead to permanent changes due to presence of gravitational waves and have been termed as persistent observables in [16]. Our work mainly focuses on one such observable named displacement memory observable and tries to calculate the memory in non-flat Kundt geometries.
Memory effects in non-flat backgrounds in GR have been studied in both dS [17,18] (motivations from cosmology) and AdS spacetimes [19]. In [19], the authors have showed how to isolate the gravitational wave contribution from the background spacetime by resorting to Fermi normal coordinates and solving the geodesic deviation equation. They treated the wave as a perturbation over AdS spacetime caused due to scattering of massive/massless particles. In our work, we adopt the same method for studying memory effects in Kundt spacetimes. However, in our case, the setting is non-perturbative, since we deal with exact spacetimes representing gravitational waves.
Kundt spacetimes are exact radiative geometries consisting of non-expanding, non-shearing and non-twisting null geodesic congruences (NGC) [20][21][22]. They admit various wave solutions (pp waves, Siklos waves [23,24]) related to the presence of NGC whose tangent vector is usually not covariantly constant. In general, the wave surfaces may not be Cartesian planes. This non-planarity can signal the presence of matter or a cosmological constant [25]. Gyratons (spinning relativistic sources) are solutions obtained as a sub-class of Kundt geometries [26][27][28]. Presence of gyratonic matter in a Kundt geometry imparts an angular momentum due to its intrinsic spin. Till date, most of the research around Kundt geometries have largely been focused within the realm of GR [29][30][31][32][33][34], though there is some recent work in Gauss-Bonnet [35] and quadratic gravity [36]. As far as we know, not much has been done on such geometries within the ambit of scalar tensor theories. Our article is one such attempt towards (a) contructing new Kundt-type solutions in the simplest scalar-tensor theory, BD gravity and, more importantly, (b) verifying/contrasting the presence/absence of memory effects w.r.t similar ones in GR.
There does exist previous work on memory effects in BD theory which are different from what we wish to pursue here. As is well known, in linearised gravity, the BD scalar field produces a breathing mode along with the two additional polarizations (+, ×) found in GR [37]. Lang computed GW waveforms for scalar and tensor modes separately in the PN approximation [38,39]. Du and Nishizawa proposed a test of gravity for scalar tensor theories [40]. They found two distinct sets of memory contributions: T-memory (tensor) and S-memory (scalar). Scalar memory is unique in such theories and was used as a tool to understand the Vainshtein screening mechanism in BD gravity [41]. Asymptotically flat spacetimes in BD theory and memory have been recently studied in [42,43]. The BMS group [44] is retained for the tensorial case. There are degenerate vacua for the scalar sector related via Lorentz transformations. The BMS charge algebra is computed in [45].
Studying memory effects for such Kundt wave spacetimes in GR was initiated by two of us in [46] through analysis of geodesic motion. Similar to exact plane wave spacetimes, one can construct sandwich waves here by choosing appropriate limiting profiles [14,15,22]. This serves as a qualitative toy model of a gravitational wave burst. Interesting distinctions occur between negative and positive constant curvature solutions, particularly for the latter, where we found a new frequency memory effect. In our analysis here, we consider memory in Kundt geometries without and with gyratonic terms. First, we construct explicit Kundt solutions for arbitrary ω. The case ω = −2 is special because it yields a spacetime with constant negative scalar curvature. For other ω we have variable positive or negative curvature. After constructing the solutions, we analyse geodesics with the intent of studying displacement memory. Thereafter, we move on to geodesic deviation. We do this by isolating the background, wave and gyratonic contributions to the deviation vector and their evolution equations. The coupled system of equations are then solved to obtain the behaviour of the deviation which helps us in analysing the presence of memory. We will also see how displacement memory is related to the BD scalar field and memory obtained via deviation shows the contributions of the background, gravitational wave and/or the gyratonic parts in the total deviation and hence, memory.
In our approach towards solving the deviation equation, our calculations are done in Fermi normal coordinates [19]. Here the coordinate system is locally Minkowskian and hence the notion of displacement and velocity memory effect is qualitatively similar to exact plane wave spacetimes [15,47,48]. In such Fermi coordinates, we construct parallely transported tetrads along a given timelike geodesic. The deviation vector is obtained w.r.t. the tetrad directions and then transformed back to the coordinate basis. We demonstrate the calculations in several specific cases (including an exactly solvable example) in the relevant section below.
The organization of our paper is as follows. In Sec. II we lay out the basic framework and the tools used in the paper. Section III deals with Kundt wave solutions without gyratonic terms and memory effects. Section IV covers the Kundt metric with gyratonic terms and aspects of memory. We summarise our work in Sec V with comments on possible extensions.
Relevant mathematical formulae used in the paper are listed in an Appendix.

A. Brans-Dicke gravity
Brans and Dicke, seeking motivation from Mach's principle, proposed their alternative theory [49], where the Newtonian gravitational constant (G) is related to the reciprocal of a scalar field. This link is based on the idea of variability of inertial mass at different points in spacetime. The action for the BD theory in the Jordan frame is given below.
Here, φ denotes the ambient scalar field, ω is the BD parameter and the L m denotes the matter Lagrangian. ω is a dimensionless parameter which is highly constrained from Solar System observations [50]. Different values of ω correspond to different theories.
Since its arrival on the scene in the early sixties, BD theory has been ruled out at times but has reappeared in different avatars serving diverse needs. For example, its ω = −1 limit is dilaton gravity, which emerges from string theory [51]. Similarly, extensions such as replacing a constant ω with ω(φ) leads to a broader class of theories which are actively pursued today in cosmological and astrophysical contexts [52]. Much of the relevance of BD, as well as scalar-tensor theories, rests on providing templates for comparative studies with GR, with the hope of discovering the new physics embodied in the theoretical constructs or ruling them out. Our work here, is also an attempt in a similar direction vis-a-vis the memory effect.
As mentioned before, we study Kundt geometries [20][21][22] in BD theory. Given our primary motivation with regard to distinguishing memory effects arising in GR and BD theory we first construct the solutions (Kundt-type). In our theoretical setup, we do not restrict the value of ω. Instead we solve the field equations for novel Kundt geometries in BD theory and choose the value of ω to study specific cases. Interestingly, we find that in both the The field equations are obtained by variation of g µν and φ. After performing a little algebra, we can write them in the standard form as shown below.
The box operator is constructed using the Kundt spacetime metric.

B. Kundt geometries and the geodesic analysis of memory
Let us first introduce the class known as as Kundt geometries. We also point out how geodesic analysis leads to our understanding of memory effects for such geometries. Later in Sec. III and IV we provide explicit calculations of memory for various solutions in BD theory which are special cases of the spacetimes mentioned below.

Kundt wave metric
The line element for a Kundt wave geometry is given as, The waves (denoted via the term H(u, x, y)) are viewed as propagating in the curved background spacetime [25,30,46]. The background curvature is dependent on P (u, x, y).

Kundt metric with gyraton/gyraton-like terms
The line element for a generalisation of the metric above is given as, The vector field k = ∂ v gives the NGC. The tangent to the spatial surfaces (P ∂ x , P ∂ y ) and k are orthogonal to each other. In this paper, we specifically work with gyratonic spacetimes Thereafter, we study geodesics and geodesic deviation to infer about memory effects.

Geodesic analysis of memory
Zhang et al. [14,47] have recently studied memory effects by analysing the evolution of geodesics in exact plane wave spacetimes. By choosing a Gaussian pulse for the polarization (radiative) term in the line element, they solved the geodesic equations numerically. The change in separation and velocity caused due to the passage of such a pulse was termed displacement and velocity memory effect. We further extended their analysis on predicting memory from geodesic analysis, for Kundt wave spacetimes in GR [46]. Expectedly, it was found that there exists a link between memory effects and the wave profile/curvature of the background spacetime. In order to test this hypothesis in an alternative theory of gravity, we construct and investigate memory for solutions in BD gravity with constant (negative) and varying (positive) curvature scalar. This indicates a direct comparison with results on memory in GR within the Kundt class of geometries.

C. Memory effects and geodesic deviation equation
Apart from a geodesic analysis, one can understand memory from geodesic deviation. In one of the seminal works on memory effect, Braginsky and Grishchuk [5] studied geodesic deviation between two test particles moving in a weak gravitational wave. The set of equations following from geodesic deviation was recast into an equation of a forced system. The forcing was shown to be caused due to the passage of the gravitational wave. Integrating the forced equation one could obtain the change in separation between the two particles. This separation was caused only due to the gravitational wave and hence was a residual change (memory) imprinted on the spacetime.
In general, the deviation between two neighbouring geodesics in curved spacetimes gets contributions from both gravitational wave part and the non-flat background. Thus, recovering the deviation solely due to the wave (to find the memory effect) is nontrivial. A step along this direction of calculating memory effect in AdS spacetime was initiated in [19].
They considered linearized gravitational wave perturbations over AdS background spacetime. Constructing a parallelly propagated tetrad along a timelike geodesic, they employed Fermi normal coordinates to separate the gravitational wave part from the background. We use the same method as given in [19] for our analysis. The only difference here is that in our case the spacetime itself is a radiative geometry while in the former case the perturbations act over a nonradiative background (i.e. AdS). Fermi coordinates were also used in [48] to find the memory for gravitational shock waves and gyratonic pp-waves. Physically it corresponds to a set of timelike inertial observers traversing along the central geodesic whose spatial distance scales are very small compared to the curvature length scale.
After obtaining the geodesic deviation in the tetrad basis, we transform them back to the coordinate basis. Thus, we determine separately, deviation arising due to different sources (e.g. background or the wave) in the coordinate basis. The total deviation will, of course, be qualitatively similar to results obtained from the geodesic analysis.
Let us now explicitly define Fermi normal coordinates and parallelly transported tetrads, which we will use to arrive at memory effects. Consider the geodesic deviation equation Here, ξ µ is the deviation vector between neighbouring timelike geodesics. U µ is the tangent vector along one of the geodesics. R µ νρσ being the Riemann curvature tensor. Along a chosen geodesic one can set up a coordinate system {t, Z i } such that the Christoffel connections are always zero along that curve. The spacetime curvature manifests itself through the Riemann curvature. Such basis sets are known as Fermi bases [53]. Along the geodesic, t denotes the proper time and Z i = 0. A parallely transported tetrad is denoted as e µ a . This satisfies Construction of such a tetrad depends on the spacetime metric. In the tetrad e µ a , the Greek indices are the spacetime coordinates while Latin indices are for the Fermi coordinates. We have e µ 0 = U µ (tangent vector) due to the geodesic equation. In the vicinity of the central geodesic, a neighbouring geodesic having separation ξ µ is re-written in the Fermi basis as, Hence, the geodesic deviation equation in Fermi coordinates (t = X 0 , X i ) becomes: t and λ both being affine parameters are related via affine transformations. Eq.(9) can be obtained by substituting Eq.(8) (and also using Eq. (7)) in Eq.(6). The spatial indices associated with the frame are denoted by i, j (R i 0j0 = R µ νρσ e i µ e ν 0 e ρ j e σ 0 ). The tetrads and metric are related via 1 η ij = e α i e β j g αβ . Since the deviation has contributions both from background and gravitational radiation, we assume that the total deviation vector is decomposed in the form: where the suffixes B, W are for background and wave respectively. As already mentioned, a similar analysis has been carried out in [19] to separate the radiation from background curvature.
The splitting of the Riemann tensor is done by noting the terms which are proportional . Such terms denote the gravitational wave contribution while the other terms are due to background curvature.
Thus, Eq.(9) separates into the two equations shown below.
Eq. (10)  Eq. (11) gives the memory effect in the tetrad frame. Once both background and wave deviation (i.e.Z i B and Z i W ) are known, we revert back to the coordinate basis using Eq.(8) to obtain ξ µ B , ξ µ W and ξ µ (= ξ µ B + ξ µ W ). In the case of gyratons, we split the deviation vector as: Apart from the familiar terms Z i B and Z i W , we also have deviation due to gyratons denoted by Z i G . The Riemann tensor corresponding to this deviation comprises of terms linked to W i (u, x, y) given in Eq. (5). The deviation equations corresponding to the three distinct effects are given by the Eqs.(64), (65) and (66) in Sec. IV B. We then go on to calculate ξ µ B , ξ µ G and ξ µ W in coordinate basis. The terms ξ µ G and ξ µ W correspond to coordinate memory effect for the gyraton and the wave respectively.
We have already emphasized the advantage of using Fermi-normal coordinates for our analysis. An important feature is that Eqs.(9), (10) and (11) take their respective forms only when the constructed tetrads are parallely transported.
It is also important to note that the results on memory from geodesic deviation are expected to match only qualitatively with those obtained from a geodesic analysis. This is because the deviation equation is perturbative by construction. In contrast there is no such restriction when we consider pairs of geodesics and differences in their separation caused by a pulse.

III. KUNDT WAVE METRIC
We now focus on finding solutions in BD theory which represent Kundt wave spacetimes generically given as in Eq.(4). The BD scalar field is assumed to be independent of v and hence, φ ≡ φ(u, x, y). The components of Eq.(2) which are relevant for solving the field equations are listed below. 2 We decompose the scalar field and the metric functions as: Adding equations (12) and (13) and using separation of variables from equation (18) results in ψ, xx +ψ, yy = 0 ( we know that G xx = G yy = 0 from the metric). The solution is, From the metric, G uv = P 2 ∆ log P (where, ∆ = (∂ xx + ∂ yy )). Using this in equation (15) gives,P Eq. (20) shows that background curvature is explicitly dependent on ω. The equations for the 'xu' and 'yu' components, as in (16), (17) and given the metric (4) we end up with Both the above equations reduce to the same equation after using the separation of variables.
We have, The Ricci scalar curvature is The component of G uu from the metric is given below.
Using equations (14) and (18) we get We set U = 1. Hence, from Eq.(23) we get α(u) as a constant . From equations (23), (25) and (19) we find that H ′ (u) is unconstrained. The xy dependent part of H(u, x, y) becomes This polarization term h(x, y) is different from GR. We will point out the consequences of this difference on the nature of the memory effect, contrasting it with GR. Thus from our generic analysis we find that only H ′ (u) and ω is unconstrained.
We now perform a coordinate transformation x = e X−a/2 cos Y, y = e X−a/2 sin Y . The metric in the new coordinates (u, v, X, Y ) becomes The reverse transformation shows that X = 1 2 (a + log[(x 2 + y 2 )]) ≡ 1 2 ψ[X]. This relation shows the imprint of the scalar field in the solution of the metric, via coordinate X. Thus, φ[X(u)] evaluated along the geodesic X(u) gives a measure of the memory effect for the scalar field due to the gravitational wave pulse.

A. Displacement memory using geodesics
The geodesic equations of coordinates X, Y for the metric line element in Eq.(27) are given below.
Geodesic equations for coordinate v is trivial (ü = 0). Hence, u acts as an affine parameter.
We try to solve Eqs. (28) and (29)  X which is given as, We use the transformation X = 1 2 q 1 ω+2 for ω = 2. Eq. (30) is thus transformed to the equation The above equation resembles that of a forced Levinson-Smith system which, generically, has an equation of the form [54] [55], x + a(x,ẋ)ẋ + g(x) = f (t) (32) where a(x,ẋ), g(x) and f (t) need to be specified. A comparison with the equation for q(u) given just above, shows the correspondence. Explicit solutions are not quite available, especially for the equation in our case, though a dynamical systems analysis exists [55].
The nature of the forcing term H ′ (u) (which encodes the effect of the gravitational wave pulse) dictates the behaviour of any solution. Therefore, it is likely that the analysis of such forced equations with nonlinearities may play a role in understanding the memory effect. with a known, nonlinear forced system, as noted above, seems worth investigating further, in future.
In principle, the geodesic equation given in Eq. (30) can be solved numerically for any value of ω for a given choice of the pulse H ′ (u). We have chosen H ′ (u) = 1 2 sech 2 u since it qualitatively resembles a gravitational wave pulse. However, as mentioned earlier, we are interested in observing memory effects for different choices of ω corresponding to different scalar curvature scenarios. The value ω = −2 is special since it is the only case where the Ricci scalar is constant and negative. For other ω one may have variable positive or negative R. We have chosen to illlustrate our analysis for the ω = −2 case with a ω value which yields a positive but variable Ricci scalar. Results for variable negative Ricci scalar are not very different and can be worked out easily too. We discuss each case briefly with corresponding plots.
• ω = −2 We consider a scenario where the scalar curvature is constant. Hence, ω = −2 and the scalar curvature is R = −8. The geodesic equation (30) becomes One can solve Eq.(33) analytically. SettingẊ/X = p, one finds that The constant is fixed by setting p = 0 at u → −∞ as initiallyẊ also vanishes. Solving X from the analytical form of p(u) yields Here, A denotes the initial position of the particle. One can check that as u → −∞, X = A whereas as u → +∞, X = 0. Thus, two geodesics starting with different initial coordinate values, eventually have a zero X value, after the passage of the pulse.
Let us now consider two different geodesics having initial X coordinate value as A 1 and A 2 . We find, as u → −∞, The final separation at u → +∞ is zero . Thus. we have displacement memory for the coordinate X.
A careful inspection reveals that Eq.(33) is invariant under u → −u. Hence, the analytical solution, X(u) = A 1+e −2u is also possible. Here, two geodesics both starting from X = 0 settle to two different final values (depending on A 1 − A 2 ). We illustrate both the analytical results below using plots.  We observe permanent displacement along X-direction (see Fig.(1)). This is qualitatively similar to GR ( see [46]) where we also observed constant separation after the passage of the gravitational wave pulse. Thus, for negative scalar curvature solutions, there is a qualitative agreement in the nature of the memory effect as found in GR and in BD (ω = −2) theory.
• ω = −2 The earlier geodesic analysis reveals that constant negative curvature case of ω = −2 is integrable. We have not been able to find analytical solutions for other values of ω. Resorting to numerics, we study the behaviour of the coordinate X for ω = +1.
The motivation for choosing ω = +1 is to compare between the results on memory for positive curvature solutions as obtained here in BD theory with those in GR. The Ricci scalar for ω = +1 becomes 4/(8X 3 ). Thus, there is a possibility of having negative scalar curvature in ω = +1 by choosing coordinate ranges where X < 0. However, in our entire analysis of memory effects, we have restricted our coordinate range beyond the singular region so that the solution conforms to the positive sign of scalar curvature.
We have thus avoided the negative scalar curvature region. Nevertheless, we have observed that the singularity does seem to influence pairs of geodesics through their gradual convergence towards it.
Eq.(30) is solved numerically in Mathematica 10 and the evolution of coordinate X is shown below. In Fig.(2), we find increasing separation between the geodesics after the departure of the pulse. This is in sharp contrast to the profiles obtained in GR. In the latter theory we found from geodesic analysis that positive curvature scenarios give rise to a frequency memory effect [46]. This is related to the different metric functions in the Kundt wave line element for the two theories. In BD theory, h(X, Y ) = log(2X) (obtained by solving the field equations) whereas in GR we took it as h(X, Y ) = 1 2 (X 2 −Y 2 ) (usual expression found in + polarization). We also find that the geodesics do not cross beyond X = 0. This is due to the singular nature of the metric as mentioned just above.

B. Geodesic deviation analysis of memory
We now turn towards discussing geodesic deviation following the method outlined in Section II C. First, we construct an orthonormal tetrad for the metric line element (27).
A similar construction was carried out in [56] for a different coordinate system. e 0 µ gives the tangent to the geodesic. The parallel transport condition (7) is only satisfied by e 0 µ (obeys the geodesic equations) and e 3 µ . Both e 1 µ , e 2 µ are not parallely transported. Hence, these two tetrads are rotated by an angleθ p = ωẎ /(2X). Since we haveẎ = 0 for all ω from the geodesic analysis, θ p is a constant. We take θ p = 0 so that the two tetrads e 1 µ , e 2 µ also satisfy Eq. (7).
The non-zero Riemann tensor components in the tetrad basis are shown below. Background Wave Substituting the expressions in Eq. (37) and Eq.(38) (for the Riemann tensor components in the tetrad basis) in Eqs. (10) and (11), we solve for the background and gravitational wave contributions to the geodesic deviation.
As pointed out earlier, we ultimately go over to the coordinate basis using Eq. (8). We find out ξ µ B , ξ µ W only along X, Y directions. The total deviation ξ µ is then obtained and compared with the results obtained from geodesics. This decomposition of deviation vectors into background and wave parts can be done because the geodesic deviation equation is linear. As the geodesic equations itself are non-linear, this method of decomposition is not possible. One may also directly integrate the geodesic deviation equation in the coordinate basis and come to similar conclusions as ours. However, using the tetrads, the equations simplify enough, as can be seen by comparing Eqs.(6) and (9). Interestingly, for ω = −2 we have an exact solution.
We now show that by solving the geodesic deviation equation we may obtain the memory effect. The change in the part of the deviation vector related to the background arises due to the P (u, x, y) term in Eq.(4). The change caused by the pulse H(u, x, y) is conventionally related to memory, largely because the pulse is viewed as the 'cause'. Our choice of Fermi coordinates simplifies the calculations to some extent, though it is surely possible to do everything in the coordinate basis as well.
Below, we use the previously chosen ω values (as in the geodesic analysis) to carry out the deviation analysis. Since u is an affine parameter we can replace the proper time t mentioned earlier, with u. ω = −2: As we have seen, the geodesic equations are analytically solvable in this case. We will use the solution X = A 1 + e −2u ,Ẏ = 0 for further investigation. In this scenario, two geodesics starting from zero initial value have two different final separations (see Fig.(1b)).
We have,B The analytical solutions for Eqs. (39) and (41) are The solutions for B(u) has hypergeometric functions denoted by 2 F 1 . C 1 , C 2 , C 3 are the constants of integration. We setK = 0 as u → −∞ to get rid of the other constant. W (u) can be easily obtained by subtracting B(u) from K(u). Since we start with zero initial velocities, we find Z 1 B and Z 1 W are constants. Eq. (43) shows that the total deviation K(u) is inverse of the geodesic solution X(u).
Reverting back to the coordinate basis, the deviation along coordinates X and Y becomes Eq. (44) shows how both the background and the wave deviation are simply proportional to X. This also follows from the Fig.(3a). Thus, the nature of the three deviations (background, wave and total) are all similar. The background and the gravitational wave both sum up to enhance the amplitude of the total deviation. Along the Y direction in Fig.(3b), we observe that the wave and background deviations cancel each other. The total deviation is a constant. Thus, the geodesic deviation analysis also gives the same qualitative result on separation, as obtained from the geodesics. This confirms our previous assertion that the geodesic analysis can only 'see' a total deviation (or separation). It is unable to retrieve the gravitational wave contribution from the total. In both plots (Fig.(3)), the blue line shows the gravitational wave memory effect.
Similar to the geodesic analysis, we perform the deviation analysis for ω = +1, which corresponds to a spacetime with a positive (but varying) Ricci scalar. The deviation equation is solved numerically. We follow the same steps as we did for ω = −2. From Eq.(37), we get that Z 1 B is a constant. The evolution equations of Z 2 B , Z 1 W and Z 2 W are solved numerically in Mathematica 10. Finally, we go over to the coordinate basis to state our results on memory. The total deviation plot in Fig.(4b) shows it to be constant, a feature also obtained from the geodesics. In Fig.(4a) for the deviation vector component along the X direction, we find that the deviation for the gravitational wave part is more than that for the background.
The net deviation grows monotonically with the onset of the pulse. This result is consistent with the geodesic separation behaviour shown in Fig.(2). The two analyses are not precisely equivalent due to the perturbative nature of geodesic deviation. Also, we note that displacement memory is seen along both directions while frequency memory is not found.

IV. KUNDT METRIC WITH GYRATON TERMS
Having discussed the nature of memory effect for the Kundt wave geometry in Brans-Dicke theory, in this section we will incorporate the contribution from the gyratonic terms in the generalized Kundt metric. As, discussed in section II B, such a metric is generically given as: where, P ≡ P (u, x, y), H ≡ H(u, x, y) and the gyratonic terms W i ≡ W i (u, x, y), Gyratonic spacetimes are ascribed to gravitational fields of spinning light beams [26,27]. The terms W i are responsible for angular momentum in the spacetime. Gauge transformations can locally set W i to zero [57]. But, globally it is not possible to remove W i and, thus, the spacetime retains its rotational nature.
As before, we first solve the field equations to obtain the metric functions and the scalar field and thereafter, using them, we study the memory effect via the geodesic equations and also, geodesic deviation.
From the G u v component of Eq.
(2) we have ω φ 2 ,v +φ φ ,vv = 0. We have considered a possible solution with φ, v = 0, i.e. the scalar field is independent of the null coordinate v. Adding the G y x and the G y y components of the Einstein tensor, we get: which is the same as Eq.(3) in presence of a traceless matter field. In our case here, since we are considering a vacuum solution, this equation is consistent.
From the G u u component we get the following equation: Similar to the previous analysis, in order to solve the above equations we decompose the metric functions and the scalar field as follows: φ(u, x, y) = α(u)ψ(x, y), Using this decomposition in Eq.(46) we get: Substituting the new functional forms in Eq.(47) we obtain: The equation that follows from the G x u component is: To solve the above equation we decompose the cross terms as follows: where J(u) is some function of u. We will see, while solving the geodesic equations that J(u) may be associated with angular momentum.
From Eq.(51) we note that Using the above result in Eq.(50) we find: Thus, α(u) and U(u) satisfy: This is the same as the Eq.(23) in the Kundt wave metric case. The last equation that will constrain the metric function is from the G uu component. It leads to, 1 8P 2 [4P 4 (H ,xx + H ,yy ) + 8P 3 H(P ,xx + P ,yy ) − 8P 2 H(P 2 ,x + P 2 ,y ) + 16P ,uu P − 32P 2 ,u ] = From both the Eqs. (48) and (51) one can see that: As in Kundt metric case we solve for U(u)=constant. This reduces Eq.(54) to: The solution is the same as for the Kundt wave case.
Thus, on solving the field equations we have obtained all the metric functions except H ′ (u) and J(u). These functions are unconstrained and are the source of the gravitational wave and the gyratonic contribution respectively. We will choose these functions judiciously while solving for the geodesics and the geodesic deviation.
Following our approach for Kundt waves, we rewrite the gyratonic metric in new coordinates defined via the coordinate transformations: x = e X−a/2 cos Y and y = e X−a/2 sin Y . As earlier, we find that the scalar ψ (realted to the φ) is twice of X.
The metric line element in coordinates (u,v,X,Y) is given as: We can see that on comparing with the Kundt wave metric, the only difference is the presence of the J(u)dudY term in the metric. We will see that this term will have its contribution to the overall memory effect.

A. Displacement memory from geodesic analysis
The geodesic equations corresponding to X and Y , for the metric line element given by Eq.(57) are as follows.
We find thatü = 0 from the geodesic equation of v. Therefore u may be chosen as as an affine parameter. From Eq.(59), it is clear that the geodesic along Y coordinate is dependent on J(u). This is different from the Kundt wave scenario where along Y there could be no evolution. Hence, Eq.(58) is dependent on both J(u) and H ′ (u). It is evident that the presence of cross terms involving the gyratonic contribution will have a significant impact on the evolution of the geodesics.
The displacement memory effect for coordinate X obtained below is equivalent to the memory due to the scalar field, as justified earlier.
Proceeding further, we note that the first integral of Eq.(59) gives, Here, C is a constant of integration which can be set to zero by noting thatẎ vanishes where J(u) = 0. Since X is functionally related to ψ[X(u)], we infer from Eq.(60) that the displacement memory along coordinate Y is both due to scalar field and the gyraton. As in the case without the gyraton term, we can, for ω = −2, convert the equation for X to that of a generalised Levinson-Smith system with the added feature that function g(x) in (32), is now g(x, t). We obtain, Asymptotic behaviour of Eq. (61) has been studied in [58].  The geodesic evolution show a constant shift along u as well as a displacement memory effect along both X and Y directions. Apart from the shift along u (which is present in the evolution of both X and Y ), Fig.(5b), i.e. the evolution of Y , also shows the role the gyratonic terms in the metric. The plot along the X direction in Fig.(5a) is similar (but shifted along u) to the Kundt wave case.
• ω = −2 Since the Ricci scalar is the same as for the Kundt wave metric, we find that for ω = 2 there is a singularity at X = 0. We have carefully restricted our geodesic analysis in the domain X > 0. The plots for the geodesics are shown below. From Fig.(6) we see that the separation slowly builds up along the X direction. We note from Eq.(60) thatẎ = J(u)/(4X). Since J(u) is a sech-squared pulse, we find thatẎ is zero at u → ±∞. The only significant contribution of the gyratonic pulse is centered at u = 0. Thus, for the gyratonic Kundt metric, both the X and Y coordinates have memory effect around u = 0 in contrast to memory only along the X direction for Kundt waves.
From the plots of the geodesics it is clear that the most visible contribution of the gyratonic term is in the evolution along for Y . As for the evolution along X, the plots are shifted because of the additional term dependent onẎ 2 ∝ J[u] 2 in the equation (58).

B. Geodesic deviation analysis of the memory effect
The orthonormal tetrads for the gyratonic Kundt metric (57) are written below.
The gyratonic term J(u) is only present in e 2 µ . We now check whether the constructed tetrads are parallely transported, i.e. (U µ = e µ 0 ).
Therefore, if we rotate e µ 1 and e µ 2 by an angle θ p , such that it satisfiesθ p = ωẎ 2X , we will obtain a parallely propagating tetrad. Thus, θ p is dependent on the gyratonic term J(u), following from Eq.(60). In the new tetrad basis, we calculate the nonvanishing Riemann components and separate them into background, gyratonic and the gravitational wave parts.
The deviation equation for the three parts are The deviation equations for the gyratonic Kundt metric are straightforward generalizations of Eqs. (10) and (11) where there was no contribution coming from the term (R i 0j0 ) G . The  Let us now summarise pointwise, the main results obtained in this article.
• In vacuum Brans-Dicke theory we have constructed solutions representing a Kundt-type line element without and with gyratonic terms. The solutions have the expected feature of two unspecified functions H ′ (u) and J(u) which represent the wave-profile and the gyraton term respectively. Both these functions can be specified while writing down explicit solutions.
We choose them to be localised fucntions (eg. proportional to sech 2 u) with the aim of studying memory effects.
• For both line elements with chosen profiles for the H ′ (u) and J(u) we first study geodesics and then geodesic deviation in order to obtain memory effects. Logical similarities between results using geodesics and geodesic deviation are visible in our results. Displacement memory via geodesics and memory found using solutions of the deviation equation explicitly demonstrate the dependencies on the presence and profiles of the functions H ′ and J. The role of the gyraton term is clearly visible in the memory effects.
• Displacement memory for ω = −2 and without the gyraton term is exactly solvable and easy to understand analytically-a feature we show while analysing this specific case with constant negative Ricci scalar.
• In our analysis of deviation we have used a parallely transported Fermi basis where the equations simplify. We have split the deviation vector into background, gravitational wave and gyratonic parts and have shown how each part influences the behaviour of the total deviation vector, thereby ensuring the existence of memory.
• The identification of the main geodesic equation for general ω with a known dynamical system, the Levinson-Smith system, is an interesting observation in our work. The lack of exact solutions of such systems was a hindrance in using this correspondence for understanding memory. As mentioned earlier, it may be possible to extend and enrich our work along these lines using inputs from mathematics and the theory of dynamical systems.
The memory effects shown here for the two different Kundt metrics are characteristic features of BD theory itself. This is apparent through the dependencies of memory on ω and/or the BD scalar field. Thus, there is an intrinsic difference with similar scenarios in GR at a qualitative as well as a quantitative level.
The geodesic deviation analysis carried out here (using the Fermi basis etc.) can readily be applied to calculate memory effects for any spacetime, in particular those where radiative behaviour is present. Further, one may also employ the B-matrix formalism as introduced in [59] to study the behaviour of the kinematic variables of timelike geodesic congruences.
One may also search for memories in impulsive gravitational wave spacetimes [60,61] and compare with the results obtained here, by setting appropriate limits.
Finally, it may be worthwhile identifying Kundt geometries in other alternative theories so that a comparison can be made on the nature and distinguishing features of the memory effects manifest in each such theory. Even within BD theory, there exists scope of finding newer Kundt spacetimes for which memory effects can always be investigated using the methods outlined here. We hope to return to these issues in future work.

APPENDIX Einstein tensor components
The Einstein tensor components of the generalized Kundt metric with gyratons are given below.
Note that by setting W 1 = W 2 = 0 we get the Einstein tensors for the Kundt wave metric given in Eq.(4).

Riemann tensor in tetrad frame for Kundt spacetimes with gyraton terms
The Riemann tensor in the parallely propagated tetrad frame for the gyratonic Kundt metric are as follows.