Gravitational memory effects and Bondi-Metzner-Sachs symmetry in scalar-tensor theories

The relation between gravitational memory effects and Bondi-Metzner-Sachs symmetry of the asymptotically flat spacetimes is studied in the scalar-tensor theory. For this purpose, the solutions to the equations of motion near the future null infinity are obtained in the generalized Bondi-Sachs coordinates with a suitable determinant condition. It turns out that the Bondi-Metzner-Sachs group is also a semi-direct product of an infinitesimal dimensional supertranslation group and the Lorentz group as in general relativity. There are also degenerate vacua in both the tensor and the scalar sectors in the scalar-tensor theory. The supertranslation relates the vacua in the tensor sector, while in the scalar sector, it is the Lorentz transformation that transforms the vacua to each other. So there are the tensor memory effect similar to the one in general relativity, and the scalar memory effect, which is new. The evolution equations for the Bondi mass and angular momentum aspects suggest that the null energy fluxes and the angular momentum fluxes across the null infinity induce the transition among the vacua in the tensor and the scalar sectors, respectively.


Introduction
Gravitational memory effect is a permanent displacement of test particles after the passage of gravitational waves (GWs) discovered in general relativity (GR) a long time ago [1][2][3]. It is caused by the energy fluxes traveling through the (future) null infinity I + in an asymptotically flat spacetime [4], which possesses an asymptotic symmetry group, called Bondi-Metzner-Sachs (BMS) group [5][6][7]. This group is a semi-direct product of an infinite dimensional supertranslation group and the Lorentz group. Under the action of supertranslations, the classical vacuum is degenerate. The energy fluxes induce the transition among the degenerate vacua and result in the memory effect. Memory effect gains an increasing interest due to the study of the infrared structure of gravity [4,[8][9][10] and the detection of GWs [11][12][13][14][15][16][17][18]. The former reveals the infrared triangle -the equivalence among memory effect, soft theorem [19] and asymptotic symmetry. It also stimulated the works [20][21][22][23][24][25] on the asymptotic symmetry on the black hole horizon and black hole information paradox [26]. The later opens up the possibility to examine the existence of this effect experimentally in the near future [27][28][29][30][31][32][33].
Alternative metric theories of gravity should also have memory effect, since they have two degrees of freedom like the two in GR, at least. In this work, we would like to study the memory effect in the scalar-tensor theory, the simplest alternative to GR. Unfortunately, there are essentially infinitely many scalar-tensor theories, so we would like to concentrate on the simplest one, i.e., the Brans-Dicke theory (BD) [34]. In spite of its simplicity, BD should be a good representative. This is because the more complicated scalar-tensor theories, such as Horndeski gravity [35], Degenerate Higher-Order Scalar-Tensor theories [36][37][38] and spatially covariant gravity [39,40], contain more general interactions for the metric field and the scalar field ϕ. None of the interactions should have impact on the asymptotic symmetry or the degenerate vacua. Therefore, BD should be a good example to study the memory effect in scalar-tensor theories.
As discussed in several previous works, there are three degrees of freedom in the scalartensor theory. Two of them excite the familiar + and × polarizations as in GR, and the third one, ϕ, excites the transverse breathing polarization, if it is massless [41][42][43]. Thus, one may call the GW corresponding to the + and × polarizations the tensor wave, and the one to the breathing polarization the scalar wave. There are thus two sectors in the scalar-tensor theory: the tensor and the scalar sectors. One expects that there should be memory effects in both sectors.
In this work, we would like to find out gravitational memory effects in these two sectors by solving the equations of motion of BD using a generalized Bondi-Sachs coordinate system {u, r, θ, φ} [44], with r proportional to the luminosity radius asymptotically (r → +∞). The resultant metric components and the scalar field ϕ will be expanded in powers of 1/r. The BMS symmetry is then determined by looking for the coordinate transformations that preserve the asymptotic behaviors of the metric components and ϕ. It turns out that the BMS transformations in BD are also parameterized by the arbitrary functions α on a unit 2-sphere, and the smooth conformal vector fields Y A of the round 2-metric with A = θ, φ. Therefore, analogously to GR, the BMS transformations in BD constitute a group, which is also a semi-direct product of the supertranslations (parameterized by α) and the Lorentz group (parameterized by Y A ). Then, once the vacua in the tensor and the scalar sectors are properly defined, one discovers that the vacua in both sectors are degenerate. Different vacua in the tensor sector are related to each other by supertranslations, as in GR, while those in the scalar sector by Lorentz transformations. Therefore, similarly to GR, the transition among the tensor vacua is induced by the null energy fluxes through I + , which include the energy flux of the scalar GW if it exists. This results in the tensor memory effect like the one in GR. However, the transition among the scalar vacua is induced by the angular momentum fluxes through I + . As a consequence, the scalar memory effect caused by this transition is absent in GR.
The memory effect mentioned above is also specifically called the displacement memory effect, as opposed to the spin memory, which is relatively a new effect in GR discovered in ref. [9]. This effect causes the time delay between two counter-orbiting massless particles in a circular orbit, and might be detected by LISA. From the results in section 2.2, one knows that there is also the spin memory effect in BD, and it is caused by the tensor GW. The scalar GW does not have this effect. So in this work, we will not elaborate on spin memory effect.
In the past, several works investigated the memory effect in scalar-tensor theories. Lang calculated the waveform of the GW generated by a compact binary system [45,46]. The tensor GW was calculated up to the second post-Newtonian (PN) order beyond the quadrupole approximation, while the scalar GW waveform is accurate up to 1.5 PN order. The memory effect occurs both in the tensor and the scalar sectors. Du and Nishizawa proposed the use of the memory effect as a probe to the nature of gravity, taking the scalartensor theory as an example [47]. The gravitational memory effect in massive gravity, higher derivative and infinite derivative gravity theories was also studied in refs. [48][49][50].
However, none of these works deciphered the relation between the memory effect and the BMS symmetry.
This work is organized in the following way. First, the equations of motion are solved in the generalized Bondi-Sachs coordinates in section 2. For that end, the boundary conditions for the metric components and ϕ are determined in section 2.1. Then, the solution is solved for in section 2.2, and this solution is used to calculate the geodesic deviation equation to reveal the GW polarization content of BD, establishing the connection to the previous works [41][42][43], in section 2.3. Second, based on the solution, the BMS transformations are determined in section 3, and the transformation properties of various quantities are also discussed there. The next section 4 focuses on the relation between the memory effect and the BMS symmetry. After that, there is a brief conclusion in section 5. In this work, the abstract index notation is used [51], and the speed of light in vacuum is c = 1. Most of the calculation was done with the help of xAct [52].

Asymptotically flat spacetimes in Brans-Dicke theory
In this section, the asymptotically flat spacetime, i.e., an isolated system, in Brans-Dicke theory is studied. The action of BD is given by [34], where ω is a constant, G 0 is the bare gravitational constant, and S m is the matter action. The variational principle leads to the following equations of motion, δg ab is the matter stress-energy tensor, T = g ab T ab is its trace, and T ab is the effective stress-energy tensor for ϕ, given by In fact, eq. (2.2c) follows from eqs. (2.2a) and (2.2b).
In a suitable set of coordinate system, e.g., the generalized Bondi-Sachs coordinates (u, r, x 2 = θ, x 3 = φ), the line element is [53] where β, V, U A and h AB (A, B = 2, 3) are six metric functions. After solving the equations of motion, these metric functions and the scalar field ϕ can be expanded in powers of 1/r. For that, one also needs know the boundary conditions, i.e., the asymptotic behavior of g ab and ϕ.

Boundary conditions
In GR, it is well-known that [5,6,53] where γ AB is the round metric on a unit 2 sphere, In addition, the determinant of h AB is By this, r is the luminosity radius [5,6]. In terms of the metric functions, one knows that However, in BD, the asymptotic behavior of g ab might be different. In fact, one can assume that ϕ = ϕ 0 + O r −1 with ϕ 0 a constant [54][55][56]. Then, the terms in the second round brackets in eq. (2.2a) begin at the order of O r −1 , while in GR, the matter stress-energy tensor T ab begins at one order higher, as T ab is usually quadratic in fields. Therefore, the asymptotic behavior of g ab in BD should be analyzed. One can perform a conformal transformationg ab = ϕg ab and define ϕ = eφ such that whereT =g ab T ab . Equation (2.9a) takes the similar form as in GR with a canonical scalar field √ 2ω + 3φ. Ifφ =φ 0 + O r −1 , again with a constantφ 0 , one expects thatg ab would have the same behavior as the one in GR. Therefore, one postulates that in BD, So r is no longer the luminosity radius, but it approaches a quantity proportional to the luminosity radius near I + . Comparing eqs. (2.5a) with (2.10a), one finds out that the condition on g ur is relaxed, so in BD, β = O r −1 . The remaining metric functions have exactly the same asymptotic behavior as those in GR, except that the determinant conditions (2.7) and (2.11) are different.

Asymptotic solutions
In order to solve the equations of motion near the null infinity, one can first assume Here, all of the expansion coefficients are functions of (u, x A ), and D A is the compatible covariant derivative of γ AB . The indices of c AB and d AB are lowered and raised by γ AB and its inverse γ AB . One can check that, the determinant condition (2.11) implies that where γ ABĉ AB = γ ABd AB = 0. Soĉ AB is transverse and traceless, and contains 2 degrees of freedom. This is nothing but the transverse-traceless part of the metric perturbation above the flat spacetime background [41][42][43]. The remaining degree of freedom is encoded in ϕ 1 . In the end, the components of the matter stress-energy tensor are assumed to have the following forms, as suggested in Ref. [57], based on the behavior of the radiative canonical scalar field in the flat spacetime. Here, again, all of the expansion coefficients are functions of (u, x A ). In the following, it is assumed that the matter stress-energy tensor is given. Following the method presented in Ref. [53], one defines the following quantities, k AB is the extrinsic curvature tensor of the null hypersurface u = const.. In terms of these quantities and the metric functions, the Christoffel symbols can be expressed, but are very complicated, which can be found in ref. [53].
The equations of motion can be solved in a special order [5,6,53,58]. First, one starts with the Einstein's equation (2.2a). The simplest one is the rr-component, which is From this, one can solve for ∂ r β, and further integrate the resulting equation to give rise to β, which has the following form, The rA-component of eq. (2.2a) is also very simple, which can be rearranged to have the following form where D A is the compatible covariant derivative for h AB . Integrating this equation leads to the series expansion for n A , 20) where N A is an integration function and can be named the Bondi angular momentum aspect, and By the definition (2.31), one can obtain the series expansion for U A , i.e., in a suitable way to obtain a fairly manageable expression, where R is the Ricci scalar for h AB , and the series expansions for β, φ, U A , n A , and h AB have been submitted to obtain the last line. With this, V can be solved for, that is, where m is an integration function, called the Bondi mass aspect [58], and the subscript r of the round brackets means to evaluate the integrand at r . The form of this equation suggests that one can replace V in the integration by the left-hand side, and repeat for enough times to obtain a series expansion for V truncated at a desired order. Thus, one finds out that Next, one considers the AB-component of eq. (2.2a), which is even more complicated. After some tedious manipulation, one arrives at where R A B is the Ricci tensor for h AB with one index raised, and The both sides of eq. (2.26) can be rewritten using the series expansions of the metric functions and ϕ, i.e., (2.28) It turns out that the terms in the brackets on the right-hand side add up to zero, so One can also try to integrate eq. (2.26), 30) with N A B a new integration function, and replaceL C D in the brackets with the right-hand side of the equation several times in order to obtain the series expansion forL A B up to an appropriate order. It turns out that, this treatment leads to an simple result, that is, which is called the news tensor [59]. This is also a symmetric and traceless tensor, i.e., N AB = N BA and γ AB N AB = 0. So it has two degrees of freedom, likeĉ AB . The nonvanishing of it signals the presence of the tensor GW.
Up to now, the series expansions for all the metric functions have been determined, i.e.,  32) where N =φ 1 , and dot means the partial u derivative. Second, the uA-component of eq. (2.2a) gives rise to the evolution of the Bondi angular momentum aspect, (2.33) Third, the equation of motion (2.2b) for ϕ leads tȯ 34) Similarly to N AB , a nonvanishing N implies that there is the scalar GW at I + . Finally, eq. (2.2c) results in The equations giving rise to the above results are extremely complicated and not very illuminating, so none of them is displayed.
Reflecting upon the computation, one realizes that the news tensor N AB and the scalar N =φ 1 play special roles. As long as N AB (u, x C ) and N (u, x A ) are prescribed on I + , one can determineĉ AB (u, x C ) and ϕ 1 (u, x A ), given their initial valuesĉ AB (u 0 , x C ) and ϕ 1 (u 0 , x A ) at a certain retarded time u 0 . Then by eqs. (2.32) and (2.33), the evolutions of the mass and the angular momentum aspects, m(u, x A ) and N A (u, x B ), are known with the knowledge of their initial values m(u 0 , x A ) and N A (u 0 , x B ), respectively. Similarly, eq. (2.29) also gives the evolution ofd AB (u, x C ), again once its initial valued AB (u 0 , x C ) is given. In the end, eq. (2.34) can now be used to fix ϕ 2 (u, x A ) with ϕ 2 (u 0 , x A ) provided. Indeed, by their definitions, N AB and N are directly associated with the degrees of freedom in BD.
Therefore, the metric components are completely determined up to certain orders in 1/r, (2.36d) Of course, ϕ has been expanded in eq. (2.12a) with the expansion coefficients satisfying eq. (2.34). Note the presence of the term with ln r in g uA . This leads to the existence of singularities on the unit 2-sphere [53]. To avoid that, we would like to set U A = 0, so When the scalar profile is trivial, i.e., ϕ n = 0 for n ≥ 1, the metric components (2.36) reduce to the GR's results [53], as one would expect. Now, a close inspection reveals that the coefficient of the 1/r term of g uu is m = m+ϕ 1 /2ϕ 0 . Here, we do not call m the Bondi mass aspect as in GR, because the evolution equation (2.32) suggests that m is consistent the original meaning of the Bondi mass aspect. Indeed, it is well-known that the outgoing GWs carry energy away from an isolated system, so a certain "mass" of the system should decrease. Now, define the following masses where the integrations are over the unit 2-sphere. Then, M decays according tȯ in which the second term in eq. (2.32) drops out because it is a total divergence. The first term in the square brackets above is the energy flux of the matter fields, such as the electromagnetic field. The second one represents the energy flux of the tensor GW, and the third is for the scalar GW [54,56]. However, the rate of the decrease in M (u) iṡ so the integrand in the second term is linear in N , which is not the energy flux of any sort. Therefore, this justifies that m should be called the Bondi mass aspect. In fact, M is the ADM mass [60], which can be understood by performing the coordinate transformation u = t − r + O r −1 .

Gravitational wave polarizations
In the previous sections, the asymptotic solution of an isolated system in BD has been calculated. It was claimed that when N AB and N are nonvanishing, there exist GWs at I + . Here, one can verify this statement. Previously, several works discussed GWs and the polarizations in scalar-tensor theories [41-43, 61, 62]. It is claimed that in a scalar-tensor theory with a massless scalar field, there are three polarizations. Two of them are the familiar + and × polarizations, and the third one is the transverse breathing polarization. Now, let us examine what they are in terms of N AB and N , orĉ AB and ϕ 1 . In fact, by simply counting the degrees of freedom, one immediately realizes thatĉ AB excites the + and × polarizations, while ϕ 1 the transverse breathing one. The GW is detected based on the geodesic deviation equation [51], where T a is the 4-velocity of a freely falling test particle, and S a is the deviation vector between adjacent test particles. To detect the GW of an isolated system, one places test particles far away from the isolated system. There are many ways to choose test particles. One of them is to place them at a fixed radius r 0 and the fixed direction x A 0 , and these particles are called BMS detectors [4]. In general, these particles will be accelerated, but as long as r 0 is very large, they approximately freely fall. In terms of the coordinates, one has [4] X u,r BMS (s) = X u,r geo (s) where the subscript BMS means the coordinates are for BMS detectors, and geo for the freely falling particles. So the 4-velocity of the BMS detector is approximately T a = (∂ u ) a at a far distance r 0 . Define (eÂ) a = r −1 (∂/∂x A ) a , and one can check that they are unit spatial vectors, normal to each other and to T a , at I + . Then, let S a = SÂ(eÂ) a , and the geodesic deviation equation is approximately, where the electric part of the Riemann tensor R abc d is Here, γÂB = δÂB, in fact. It is instructive to compare eq. (2.44) with the matrix (29) in ref. [42]. To make a sensible comparison, one must understand that eq. (2.44) gives only the transverse part of the matrix (29), i.e., the top-left 2 × 2 submatrix. The remaining part of the matrix (29) is related to R uÂuẑ and R uẑuẑ withẑ direction parallel to −(∂ u ) a +(∂ r ) a near I + . These Riemann tensor components are of higher orders in 1/r, so one can ignore them for the purpose of comparison. So, one clearly recognizes that it isĉ AB that is responsible for the presence of the + and × polarizations, and ϕ 1 for the transverse breathing one. When ∂ uĉAB = 0 andφ 1 = 0, that is, N AB = 0 and N = 0, there is no change in the deviation, so there are no GWs. Now, integrating eq. (2.43) twice results in 45) in which SB 0 is the initial deviation vector at the retarded time u 0 when there were no GWs, i.e., N AB (u 0 , x C ) = 0 and N (u 0 , x A ) = 0. A radiating isolated system will eventually settle down to a state in which no GWs can be emitted, and so N AB and N vanish again. But the deviation vector SÂ may not return to its initial value, that is, This means that there might exist a permanent change in the relative distances between test particles. This is the GW memory effect. In GR, memory effect is closely related to the asymptotic symmetry that preserves the boundary conditions (2.5) and the determinant condition (2.7). In the following, the relation between the memory effect and the asymptotic symmetry in BD is discussed. One begins with the asymptotic symmetry in BD.

Asymptotic symmetry
Like GR, the scalar-tensor theory is diffeomorphic. Among all the coordinate transformations, there exists a set of transformations that preserve the boundary conditions (2.10) and (2.12a). Such a set forms a group, and is called the Bondi-Metzner-Sachs (BMS) group [5][6][7]. Let ξ µ be a vector field that generates an infinitesimal BMS transformation, then one requires that [53] L ξ g rr = L ξ g rA = 0, (3.1a) In GR, one usually demands that g AB L ξ g AB = 0, but here, we impose a different condition given by eq. (3.1b). This is suggested by the determinant condition (2.11), and would reproduce GR's condition when ϕ is trivial. By some tedious calculation, one finds out that eqs. (3.1a) and (3.1b) determine the components of ξ µ , i.e., where f and Y A are arbitrary integration functions. Next, eq. (3.1c) implies the following properties of the integration functions, Finally, eqs. (3.1d) and (3.1e) lead to no new conditions on ξ µ . By eq. (3.3a), one obtains that with α(x A ) an integration function of x A . Equation (3.3b) shows that Y A depends only on x A , and it is actually a conformal Killing vector field for γ AB according to eq. (3.3c). Formally, eqs. (3.2a) and (3.2b) coincide with those in GR [53], which is because, in both GR and BD, it is demanded that g rr = g rA = 0 in any Bondi coordinates. Although the two theories have distinct determinant conditions, ξ r [refer to (3.2c)] still takes the same form as that in GR due to eq. (3.1b). The properties (3.3) and (3.4) satisfied by f and Y A are also the same as those in GR. Using the asymptotic expansions of the metric functions and the scalar field, one knows that Despite the fact that eq. (3.5) contains no expansion coefficients ϕ n of ϕ, and thus takes the same form as the one in GR, ϕ n show up at higher order terms in 1/r. These higher order terms are very complicated and not very illuminating, so none is displayed.
As shown by eq. (3.5), together with eq. (3.2a), an infinitesimal BMS transformation is generated by the arbitrary functions f and Y A , or α and Y A , referring to eq. (3.4). In particular, the transformation generated by α = 0 and Y A = 0 is called a supertranslation. The supertranslations form a subgroup of the BMS group, denoted by S . Among the supertranslations, those generated by with a n (n = 0, 1, 2, 3) constant are translations, and constitute the translation group T . The transformations generated by α = 0 and Y A = 0 are special. Since Y A is a conformal Killing vector field on a unit 2 sphere, if it is required to be smooth and finite, the set of such vector fields is isomorphic to the Lorentz algebra. So in this case, the transformations generated by such conformal Killing vector fields are called the Lorentz transformations, and the corresponding group is labeled by L . One may drop the requirement that Y A be smooth and finite, and the resulting transformations are called the super-rotations or superboosts [57,[63][64][65][66][67]. In this work, we will not consider the super-rotations/superboosts.
Under an infinitesimal BMS transformation ξ µ , the scalar field changes by δ ξ ϕ = L ξ ϕ, and so At the same time, the metric transforms according to δ ξ g µν = L ξ g µν , from which, we find out that is the Lie derivative, and a similar expression for L YdAB applies. Thus, the news tensor N AB and the scalar N change according to where L Y N AB has a similar expression to L YĉAB . Now, it is ready to decode the relation of memory effect, BMS symmetry and the degenerate vacua.

Memory effect and the degenerate vacua
As briefly introduced in section 2.3, memory effect is the permanent change in the deviation vector between the initial and the final states without any GWs. These states are called nonradiative [57]. In a nonradiative state, N AB = 0 and N = 0. Because of eq. (3.9), these states are invariant under any BMS transformation. The vanishing N AB and N does not implyĉ AB = 0 or ϕ 1 = 0. Ifĉ AB and ϕ 1 were zero, an arbitrary infinitesimal BMS transformation would generally result in nonvanishingĉ AB and ϕ 1 , referring to eqs. (3.7a) and (3.8a). We are more interested in the special case where the initial and the final states are both vacua. The definition of the vacuum in the scalar sector is simple, i.e., N = 0. But, for the tensor sector, it is more involved.
To define a vacuum state in the tensor sector, one sets up a Newman-Penrose (NP) tetrad {l a , n a , m a ,m a } at I + [68], andm a is the complex conjugate of m a . One can check that at I + , they satisfy the following normalization conditions, l a n a = m am a = −1, and other contractions are zero, and so these vector are all null. In addition, l a satisfies the geodesic equation. Now, one can calculate some NP variables related to the outgoing tensor GW, given by [69] Ψ 4 = C abcd n amb n cmd = − r 2 ∂ 2 uĉ ABm AmB + · · · , (4.3a) and where dots represent higher order terms, and is to take the imaginary part. The vacuum state is the one in which the leading order terms in the above expressions vanish, together with N AB = 0, in GR [59]. Here, in BD, we also define the vacuum state in the same way. So one has to solve forĉ AB such that the right-hand side of the last equation vanishes at the leading order. For this purpose, it is easier to use the standard complex coordinate ζ = e iφ cot θ 2 on the unit 2-sphere. In this coordinates, one has Therefore, the vanishing of the leading order term of the right-hand side of (4.3c) is equivalent to The general solution to this equation is [4,8] c ζζ = −2D 2 ζ C(ζ,ζ), (4.6) for some function C(ζ,ζ) on the unit 2-sphere. Thus, C(ζ,ζ) completely characterizesĉ ζζ in vacuum state.
This result has a profound implication. Let us perform an infinitesimal supertranslation with f = α(x A ) and Y A = 0 on the vacuum state, thenĉ AB transforms according to (4.7a) or, in the complex coordinate system, Therefore, the transformedĉ ζζ = −2D 2 ζ C with C = C + α is also a function on the unit 2-sphere. This means that the state described byĉ ζζ is still vacuum, so there is no a unique vacuum in the tensor sector. The degenerate vacua are related to each other via a supertranslation. However, an infinitesimal Lorentz transformation with α = 0 and Y A = 0 changes the vacuum state, because which does necessarily lead toĉ ζζ = −2D 2 ζ C for some C (ζ,ζ). So the resultant state is no longer a vacuum. Since we are now discussing the memory effect in BD, we should also study the transformation of ϕ 1 . By eq. (3.7a), an infinitesimal supertranslation does not change ϕ 1 in a vacuum state, but an infinitesimal Lorentz transformation generated by Y A does, So this means that in the scalar sector, there are also degenerate vacua with respect to the Lorentz transformation, instead of the supertranslation. Now, it is ready to discuss how memory effect is related to supertranslations and Lorentz transformations in BD. Consider an isolated system which is initially in a vacuum state before the retarded time u i , withĉ ζζ and ϕ 1 = ϕ i constant in u Then, it starts to radiate GWs, and eventually settles down to a new vacuum state withĉ ζζ and ϕ 1 = ϕ f after the time u f . Now, begin with the memory effect for the tensor GW. The above discussion shows that there exist two functions C(ζ,ζ) and C (ζ,ζ) on the unit 2-sphere such that c ζζ = −2D 2 ζ C andĉ ζζ = −2D 2 ζ C . The difference α = C − C generates an infinitesimal supertranslation, which relatesĉ ζζ toĉ ζζ in the following way, Therefore, the deviation vector changes permanently according to due to the tensor GW. This observation suggests the memory effect of the tensor GW is due to the supertranslation that induces the transition among the degenerate vacua in the tensor sector. However, the supertranslation does not change ϕ 1 as discussed above, so the memory effect caused by the scalar GW is not generated by a supertranslation. Rather, the Lorentz transformation does the work. Suppose after the radiation, the initial scalar field ϕ 1 changes to ϕ 1 . Then one may not be able to find a vector Y A such that δ Y ϕ 1 = ϕ f − ϕ i , but a finite Lorentz transformation should relate ϕ i to ϕ f . Therefore, it is the Lorentz transformation that generates the transition among the vacua in the scalar sector and thus the memory effect of the scalar GW.
With eqs. (2.32) and (2.33), one can find the changes inĉ AB and ϕ 1 . For the tensor memory effect, one integrates eq. (2.32) over u from u i to u f , and then, moves the terms to get Here, eq. (4.5) has been used. With eq. (4.10), one determines [4] α(ζ,ζ) = − dζ dζ γ ζ ζ G(ζ,ζ; ζ ,ζ ) ∆m + where G(ζ,ζ; ζ ,ζ ) is the Green's function, given by such that D 2 ζ D 2 ζ G(ζ,ζ; ζ ,ζ ) = −γ ζζ δ 2 (ζ − ζ ) + · · · . By eq. (4.13), one concludes that the tensor memory effect is due to the transition among the tensor vacua caused by the null energy fluxes passing through I + , exactly the same as in GR [4]. The null energy fluxes include the matter, the tensor GW and the scalar GW fluxes, as displayed in the square brackets in sequence. In terms of the terminology of ref. [70], the first term in the curly brackets of eq. (4.13) is the ordinary memory, and the remaining ones represent the null memory. Of course, the null memory caused by the scalar GW is new, which is called T memory inf ref. [47].
For the scalar memory effect, contract both sides of eq. (2.33) with D A , integrate over u from u i to u f , and act D −2 on both sides. Then, one obtains that where eq. (2.35) has been used, and G(ζ,ζ; ζ ,ζ ) = ln sin 2 Θ 2 is another Green's function satisfying ∂ ζ ∂ζG(ζ,ζ; ζ ,ζ ) = 2πδ 2 (ζ −ζ )−γ ζζ /2 [9], and if the terms in the square brackets were written in the ζ coordinate system, the expression would be very clumsy. So one can consider an experimentally relevant process whose initial state is some vacuum state with an initial ϕ i (andĉ ζζ ), and the final state contains a single stationary back hole. Due to the no-hair theorem, the final ϕ f = 0 [71,72]. Therefore, the left-hand side of eq. (4.14) is proportional to ϕ 2 i . The terms in the round brackets in the second line are the angular momentum fluxes of the matter, the tensor GW and the scalar GW, respectively [56]. This means that the scalar memory effect is related to the transition among the scalar vacua induced by the angular momentum fluxes through I + . The scalar memory effect is absent in GR, and is called S memory in ref. [47]. Again, in the terminology of ref. [70], one might call the effect corresponding to the round brackets the null scalar memory, and the one to the remaining terms the ordinary scalar memory.
As a final remark, we would like to brief discuss the spin memory effect in BD. By the argument in ref. [9], this effect is due to the leading order term of g uA component. According to eq. (2.36c), the leading order term of g uA contains no contributions from the scalar field ϕ. So there is also the similar spin memory effect in the tensor sector in BD to the one in GR, but in scalar sector, no new spin memory exists. Of course, this does not mean that ϕ plays no role in the spin memory. In fact, its angular momentum flux also generates the spin memory, just like other angular momentum fluxes [9]. Note that the spin memory effect is induced by the curl part of the angular momentum fluxes as clearly demonstrated by eq. (4.10) in ref. [57], while the scalar displacement memory effect is due to the divergence part of the angular momentum fluxes, referring to eq. (4.14).

Conclusion
From the above discussion, one understands that the memory effect of an isolated system in BD shares some similarities with that in GR and also has its own specialties. In both theories, there are memory effects for the tensor GW, which are induced by the passage of the null energy fluxes through I + . The supertranslations transform the degenerate vacua in the tensor sector, which makes the tensor memory effect possible. However, in BD, there exist the scalar GW. It not only contributes to the tensor memory effect by providing a new energy flux, but also has its own memory effect. This new memory effect is due to the angular momentum fluxes penetrating I + . There are also degenerate vacua in the scalar sector, which are transformed into each other via the Lorentz transformations.