3.4 Asymptotically-vanishing Maxwell fields

3.4.1 A prelude

The basic starting idea in this work is simple. It is in the generalizations and implementations where difficulties arise.

Starting in Minkowski space in a fixed given Lorentzian frame with spatial origin, the electric dipole moment −→ D E is calculated from an integral over the (localized) charge distribution. If there is a shift, −→R, in the origin, the dipole transforms as

−→ −→ −→ D ∗E = D E − q R . (114 )
If −→ D E is time dependent, we obtain the center-of-charge world line by taking −→ D ∗E = 0, i.e., from −→ −→ −1 R = D Eq. It is this idea that we want to generalize and extend to gravitational fields.

The first generalization is formal and somewhat artificial: shortly it will become quite natural. We introduce, in addition to the electric dipole moment, the magnetic dipole moment −→ D M (also obtained by an integral over the current distribution) and write

−→ −→ −→ D ℂ = D E + iD M .

By allowing the displacement −→R to take complex values, −→R ℂ, Equation (114View Equation), can be generalized to

−D→∗ = −→D − q−→R , (115 ) ℂ ℂ ℂ
so that the complex center-of-charge is given by −→ D ∗ℂ = 0 or
−→ −→ R ℂ = D ℂq−1. (116 )

We emphasize that this is done in a fixed Lorentz frame and only the origin is moved. In different Lorentz frames there will be different complex centers of charge.

Later, directly from the general asymptotic Maxwell field itself (satisfying the Maxwell equations), we define the asymptotic complex dipole moment and give its transformation law, including transformations between Lorentz frames. This yields a unique complex center of charge independent of the Lorentz frame.

3.4.2 Asymptotically-vanishing Maxwell fields: General properties

In this section, we describe how a complex center of charge for asymptotically vanishing Maxwell fields in flat spacetime can be found by using the shear-free NGCs congruences, constructed from solutions of the homogeneous good cut equation, to transform certain Maxwell field components to zero. Although this serves as a good example for our later methods in asymptotically flat spacetimes, the reader may wish to skip ahead to Section 4, where we go directly to gravitational fields in a setting of greater generality.

Our first set of applications of shear-free NGCs comes from Maxwell theory in Minkowski space. We review the general theory of the behavior of asymptotically-flat or vanishing Maxwell fields assuming throughout that there is a nonvanishing total charge, q. As stated in Section 2, the Maxwell field is described in terms of its complex tetrad components, (ϕ0,ϕ1,ϕ2). In a Bondi coordinate/tetrad system the asymptotic integration is relatively simple [35Jump To The Next Citation Point25] resulting in the radial behavior (the peeling theorem):

ϕ00 −4 ϕ0 = -3-+ O (r ), (117 ) r0 ϕ1 = ϕ1-+ O (r−3), r2 ϕ02 −2 ϕ2 = ---+ O (r ), r
where the leading coefficients of r, (0 0 0 ϕ0,ϕ1,ϕ2) satisfy the evolution equations:
˙0 0 ϕ0 + ∂ϕ1 = 0, (118 ) ˙ϕ01 + ∂ϕ02 = 0. (119 )

The formal integration procedure is to take ϕ02 as an arbitrary function of (uB,ζ, ¯ζ) (the free broadcasting data), then integrate the second, for ϕ01, with a time-independent spin-weight s = 0 function of integration and finally integrate the first, for ϕ0 0. Using a slight modification of this, namely from the spherical harmonic expansion, we obtain,

ϕ00 = ϕ00iY11i + ϕ00ijY21ij + ..., (120 ) 0 0 0 0 0 ϕ1 = q + ϕ1iY1i + ϕ 1ijY2ij + ..., (121 ) ϕ02 = ϕ02iY1−i1 + ϕ02ijY −2i1j + ..., (122 )
with the harmonic coefficients related to each other by the evolution equations:
ϕ00 = 2qηi(uret)Y11i + Qiℂj′Y21ij + ..., (123 ) √ -- ϕ0 = q + √2q-ηi′(u )Y 0 + --2Qij′′Y 0 + ..., 1 ret 1i 6 ℂ 2ij 0 i′′ −1 1 ij′′′ − 1 ϕ2 = − 2qη (uret)Y1i − 3Q ℂ Y2ij + ...
The physical meaning of the coefficients are
q = total electric charge, (124 ) i i i i qη = D ℂ = complex (electric & magnetic ) dipole moment = D E + iD M, Qijℂ = complex (electric & magnetic ) quadrupole moment,
etc. For later use, the complex dipole is written as i i D ℂ(uret) = qη (uret). Note that the i D ℂ is defined relative to a given Bondi system. This is the analogue of a given origin for the calculations of the dipole moments of Equation (114View Equation).

Later in this section it will be shown that we can find a unique complex world line, ξa(τ) = (ξ0,ξi), (the world line associated with a shear-free NGC), that is closely related to the i η (uret). From this complex world line we can define the intrinsic complex dipole moment, i i D ℐℂ = qξ (s ).UpdateJump To The Next Update Information

However, we first discuss a particular Maxwell field, F ab, where one of its eigenvectors is a tangent field to a shear-free NGC. This solution, referred to as the complex Liénard–Wiechert field is the direct generalization of the ordinary Liénard–Wiechert field. Though it is a real solution in Minkowski space, it can be thought of as arising from a complex world line in complex Minkowski space.

3.4.3 A coordinate and tetrad system attached to a shear-free NGC

The parametric form of the general NGC was given earlier by Equation (72View Equation),

xa = u (ˆla + ˆna) − Lmˆa − ¯L ˆma + (r∗ − r )ˆla. (125 ) B 0
The parameters (uB, ζ, ¯ζ) labeled the individual members of the congruence while r∗ was the affine parameter along the geodesics. An alternative interpretation of the same equation is to consider it as the coordinate transformation between the coordinates, xa (or the Bondi coordinates) and the geodesic coordinates ∗ ¯ (uB,r ,ζ, ζ). Note that these coordinates are not Bondi coordinates, though, in the limit, at + ℑ, they are. The associated (geodesic) tetrad is given as a function of these geodesic coordinates, but with Minkowskian components by Equations (74View Equation). We restrict ourselves to the special case of the coordinates and tetrad associated with the L from a shear-free NGC. Though we are dealing with a real shear-free twisting congruence, the congruence, as we saw, is generated by a complex analytic world line in the complexified Minkowski space, a a z = ξ (τ ). The complex parameter, τ, must in the end be chosen so that the ‘uB’ of Equation (90View Equation) is real. The Minkowski metric and the spin coefficients associated with this geodesic system can be calculated [27Jump To The Next Citation Point] in the (uB,r∗,ζ,ζ¯) frame. Unfortunately, it must be stated parametrically, since the τ explicitly appears via the a ξ (τ) and can not be directly eliminated. (An alternate choice of these geodesic coordinates is to use the τ instead of the uB. Unfortunately, this leads to an analytic flat metric on the complexified Minkowski space, where the real spacetime is hard to find.)

The use and insight given by this coordinate/tetrad system is illustrated by its application to a special class of Maxwell fields. We consider, as mentioned earlier, the Maxwell field where one of its principle null vectors, l∗a, (an eigenvector of the Maxwell tensor, Fabl∗a = λl∗b), is a tangent vector of a shear-free NGC. Thus, it depends on the choice of the complex world line and is therefore referred to as the complex Liénard–Wiechert field. (If the world line was real it would lead to the ordinary Liénard–Wiechert field.) We emphasize that though the source can formally be thought of as a charge moving on the complex world line, the Maxwell field is a real field on real Minkowski space. It will have a real (distributional) source at the caustics of the congruence. Physically, its behavior is very similar to real Liénard–Wiechert fields, the essential difference is that the electric dipole is now replaced by the combined electric and magnetic dipoles. The imaginary part of the world line determines the magnetic dipole moment.

In the spin-coefficient version of the Maxwell equations, using the geodesic tetrad, the choice of l∗a as the principle null vector ‘congruence’ is just the statement that

∗ ∗a ∗b ϕ0 = Fabl m = 0.

This allows a very simple exact integration of the remaining Maxwell components [35Jump To The Next Citation Point].

3.4.4 Complex Liénard–Wiechert Maxwell field

The complex Liénard–Wiechert fields (which we again emphasize are real Maxwell fields) are formally given by the (geodesic) tetrad components of the Maxwell tensor in the null geodesic coordinate system (∗ uB, r ,ζ, ¯ζ), Equation (125View Equation). As the detailed calculations are long [35Jump To The Next Citation Point] and take us too far afield, we only give an outline here. The integration of the radial Maxwell equations leads to the asymptotic behavior,

∗ ϕ 0 = 0, (126 ) ϕ ∗= ρ2 ϕ∗0, (127 ) 1∗ ∗10 2 ϕ 2 = ρ ϕ2 + O (ρ ), (128 )
ρ = − (r∗ + iΣ)−1, (129 ) -- 2iΣ = ∂¯L + L ˙¯L − ∂L − ¯LL˙.
The order expression is known in terms of (ϕ ∗0,ϕ ∗0 1 2). The function L (u ,ζ, ¯ζ) B is given by
L(uB, ζ, ¯ζ) = ∂(τ)G (τ,ζ,ζ¯), u = G (τ,ζ, ¯ζ) = ξa(τ)ˆl(ζ,ζ¯), B a
with ξa(τ) an arbitrary complex world line that determines the shear-free congruence whose tangent vectors are the Maxwell field eigenvectors.

Remark: In this case of the complex Liénard–Wiechert Maxwell field, the ξa determines the intrinsic center-of-charge world line, rather than the relative center-of-charge line.

The remaining unknowns, ϕ∗10,ϕ∗20, are determined by the last of the Maxwell equations,

∂ϕ ∗0+ L ˙ϕ∗0+ 2 ˙Lϕ∗0 = 0, (130 ) 1∗0 1∗0 1∗0 ∗0 ∂ϕ 2 + L ˙ϕ2 + ˙Lϕ2 = ϕ˙1 ,
which have been obtained from Equations (59View Equation) and (60View Equation) via the null rotation between the Bondi and geodesic tetrads and the associated Maxwell field transformation, namely,
a ∗a a ¯L- a L- a ∗−2 l → l = l − r m − r ¯m + O (r ), (131 ) L ma → m ∗a = ma − -na, (132 ) a ∗a a r n → n = n , (133 )
∗0 0 0 2 0 ϕ 0 = 0 = ϕ0 − 2Lϕ 1 + L ϕ2, (134 ) ϕ ∗0= ϕ0 − L ϕ0, (135 ) 1∗0 10 2 ϕ 2 = ϕ 2. (136 )

These remaining equations depend only on ¯ L(uB, ζ,ζ), which, in turn, is determined by a ξ (τ). In other words, the solution is driven by the complex line, a ξ (τ). As they now stand, Equations (130View Equation) appear to be difficult to solve, partially due to the implicit description of the L (uB,ζ,ζ¯).

Actually they are easily solved when the independent variables are changed, via Equation (86View Equation), from -- (uB, ζ,ζ) to the complex -- (τ,ζ,ζ). They become, after a bit of work,

2 0 ∂ (τ)(V ϕ1) = 0, (137 ) ∂ (V ϕ0) = ϕ0′, (138 ) (τ) 2 1 V = ξa′(τ)ˆla(ζ, ¯ζ), (139 )
with the solution
∗0 q −2 ϕ 1 = --V , (140 ) 2q -- ϕ ∗02 = --V −1∂(τ)[V −1∂τV ]. 2
q being the Coulomb charge.

Though we now have the exact solution, unfortunately it is in complex coordinates where virtually every term depends on the complex variable τ, via ξa(τ). This is a severe impediment to a full description and understanding of the solution in the real Minkowski space.

In order to understand its asymptotic behavior and physical content, one must transform it, via Equations (131View Equation) – (136View Equation), back to a Bondi coordinate/tetrad system. This can only be done by approximations. After a lengthy calculation [35Jump To The Next Citation Point], we find the Bondi peeling behavior

ϕ = r− 3ϕ0+ O (r−4), (141 ) 0 0 ϕ1 = r− 2ϕ01 + O (r−3), ϕ = r− 1ϕ0+ O (r−2), 2 2
( ) 0 −2 1 2 −1-- −1 ′ ϕ 0 = q LV + 2L V ∂(τ)[V V ] , (142 ) q ( -- ) ϕ01 = ---- 1 + LV ∂ (τ)[V− 1V′] , (143 ) 2V 2 -- ϕ0 = − qV −1∂(τ)[V −1V ′], (144 ) 2 2 V = ξa′ˆla(ζ, ¯ζ). (145 )
Next, treating the world line, as discussed earlier, as a small deviation from the straight line, a a ξ (τ) = τδ0, i.e., by
a ( 0 i ) ξ (τ) = τ + δξ (τ),ξ (τ ) , ξi(τ) ≪ 1, δξ0(τ) = second order.
The GCF and its inverse (see Section 6) are given, to first order, by
√ -- √ -- √ -- -- uret = 2uB = 2G = τ − --2ξi(τ)Y01i(ζ,ζ ), (146 ) √ -- 2 2 i 0 -- τ = uret +----ξ(uret)Y1i(ζ,ζ). (147 ) 2

Again to first order, Equations (142View Equation), (143View Equation) and (144View Equation) yield

0 i 1 ϕ0 = 2qξ (√uret)Y1i, (148 ) ϕ01 = q + 2qξi′(uret)Y10i, 0 i′′ − 1 ϕ2 = − 2qξ (uret)Y 1i ,
the known electromagnetic dipole field, with a Coulomb charge, q. One then has the physical interpretation of ξa(u ) : qξa(u ) ret ret is the complex dipole moment; (the electric plus ‘i’ times magnetic dipole) and a ξ (uret) is the complex center of charge, the real part being the ordinary center of charge, while the imaginary part is the ‘imaginary’ magnetic center of charge. This simple relationship between the Bondi form of the complex dipole moment, qξi(uret), and the intrinsic complex center of charge, ξa(τ), is true only at linear order. The second-order relationship is given later.

Reversing the issue, if we had instead started with an exact complex Liénard–Wiechert field but now given in a Bondi coordinate/tetrad system and performed on it the transformations, Equations (14View Equation) and (134View Equation) to the geodesic system, it would have resulted in

ϕ∗0 = 0.

This example was intended to show how physical meaning could be attached to the complex world line associated with a shear-free NGC. In this case and later in the case of asymptotically-flat spacetimes, when the GCF is singled out by either the Maxwell field or the gravitational field, it will be referred to it as a UCF. For either of the two cases, a flat-space asymptotically-vanishing Maxwell field (with nonvanishing total charge) and for a vacuum asymptotically-flat spacetime, there will be a unique UCF. In the case of the Einstein–Maxwell fields there, in general, will be two UCFs, one for each field.

3.4.5 Asymptotically vanishing Maxwell fields & shear-free NGCs

We return now to the general asymptotically-vanishing Maxwell field, Equations (117View Equation) and (120View Equation), and its transformation behavior under the null rotation around a n,

a ∗a a ¯L- a L- a −2 l → l = l − r m − rm¯ + 0(r ), (149 ) L ma → m ∗a = ma − --na, a ∗a a r n → n = n ,
with -- L (uB,ζ,ζ ) = ξa(τ )mˆa, being one of our shear-free angle fields defined by a world line, za = ξa(τ). The leading components of the Maxwell fields transform as
∗0 0 0 2 0 ϕ 0 = ϕ0 − 2Lϕ1 + L ϕ2, (150 ) ϕ ∗10= ϕ01 − Lϕ02, (151 ) ∗0 0 ϕ 2 = ϕ2. (152 )

The ‘picture’ to adopt is that the new ϕ ∗s are now given in a tetrad defined by the complex light cone (the generalized light cone) with origin on the complex world line. (This is obviously formal and physically nonsense, but mathematically quite sound, as the shear-free congruence can be thought of as having its origin on the complex line, ξa(τ).) From the physical identifications of charge, dipole moments, etc., of Equation (123View Equation), we can obtain the transformation law of these physical quantities. In particular, the l = 1 harmonic of ϕ0 0, or, equivalently, the complex dipole, transforms as

0∗ 0 0 2 0 ϕ0i = ϕ0i − 2(L ϕ1)|i + (L ϕ2)|i, (153 )
where the notation W | i means extract only the l = 1 harmonic from a Clebsch–Gordon expansion of W. A subtlety and difficulty of this extraction process is explained/clarified below.

If by some accident the Maxwell field was a complex Liénard–Wiechert field, a world line ξa(τ ) could be chosen so that ϕ∗0 = 0 0. However, though this cannot be done in general, the l = 1 harmonics of ϕ ∗0 0 can be made to vanish by the appropriate choice of the ξa(τ). This is the means by which a unique world line is chosen.

A Clarifying Comment:

An important observation, obvious but easily overlooked, concerning the spherical harmonic expansions, is that, in a certain sense, they lack uniqueness. As this issue is significant, its clarification is important.

Assume that we have a particular spin-s function on ℑ+, say, -- η(s)(uB, ζ,ζ), given in a specific Bondi coordinate system, -- (uB,ζ,ζ ), that has a harmonic expansion given, for constant uB, by

-- l,(ijk...) (s) η(s)(uB, ζ,ζ) = l,(Σijk...)η(s) (uB)Yl,(ijk...)

If exactly the same function was given on different cuts or slices, say,

-- uB = G (s, ζ,ζ), (154 )
∗ -- ( -- -) η(s)(s,ζ,ζ) = η(s) G (s, ζ,ζ),ζ,ζ ,

the harmonic expansion at constant s would be different. The new coefficients are extracted by the two-sphere integral taken at constant s:

∫ ----(s) η(∗ls,()ijk...)(s) = η∗(s)(s,ζ,ζ)Y l,(ijk...)d2S. (155 ) S2
It is in this rather obvious sense that the expansions are not unique. From Equation (97View Equation), using the real values of the complex cuts, i.e.,
( -- -) uB = GR s,Λ (s,ζ,ζ),ζ,ζ , (156 )
(in place of -- uB = G(s,ζ,ζ )), the transformation, Equation (153View Equation), and harmonic extraction yields ϕ0∗(s) 0i. This is proportional to the transformed complex dipole moment. This transformation, a functional of the form,
0∗ 0 0 0 a ϕ 0i(s) = Γ i(ϕ0,ϕ 1,ϕ2,ξ ), (157 )
is decidedly nontensorial: in fact it is very nonlocal and nonlinear.

Though it is clear that extracting 0∗ ϕ0i(s) with this relationship is available in principle, in practice it is impossible to do it exactly and all examples are done with approximations: essentially using slow motion for the complex world line.

3.4.6 The complex center of charge

The complex center of charge is defined by the vanishing of the complex dipole moment 0∗ ϕ0i(s); in other words,

Γ i(ϕ00,ϕ01,ϕ02,ξa) = 0 (158 )
determines the (up to this point) arbitrary complex world line, a ξ (τ ). In practice we do this only up to second order with the use of only the (l = 0,1,2) harmonics. The approximation we are using is to consider the charge q as zeroth order and the dipole moments and the spatial part of the complex world line as first order.

From Equation (153View Equation),

0∗ 0 0 ϕ0i = Γ i ≈ ϕ0i − 2L ϕ1|i = 0 (159 )
with the identifications, Equation (123View Equation), for q and D ℂ, we have to first order (with s ≈ uret ≈ τ),
i i i D ℂ(uret) = qη (uret) = qξ (uret). (160 )
This is exactly the same result as we obtained earlier in Equation (116View Equation), via the charge and current distributions in a fixed Lorentz frame.

Carrying this calculation [35] to second order, we find the second-order relationship between the intrinsic complex dipole and the intrinsic dipole,

Diℐ:ℂ = qξi(s), Diℂ = qηi(s), (161 ) i ξk(s) = ηk(s) + iη0I(s)ηk′(s) − -𝜖ijkηi(s)ηj′(s). (162 ) 2

In Section 5, these ideas are applied to GR, with the complex electric and magnetic dipoles being replaced by the complex combination of the mass dipole and the angular momentum.

The GCF,

uB = ξa(τ)ˆla(ζ, &tidle;ζ),

with this uniquely determined world line is referred to as the Maxwell UCF.

  Go to previous page Go up Go to next page