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 is calculated from an integral over the (localized) charge distribution. If there is a shift, , in the origin, the dipole transforms as
If is time dependent, we obtain the centerofcharge world line by taking , i.e., from . 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 (also obtained by an integral over the current distribution) and write

By allowing the displacement to take complex values, , Equation (114), can be generalized to
so that the complex centerofcharge is given by orWe 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.
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 shearfree 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 shearfree NGCs comes from Maxwell theory in Minkowski space. We review the general theory of the behavior of asymptoticallyflat or vanishing Maxwell fields assuming throughout that there is a nonvanishing total charge, . As stated in Section 2, the Maxwell field is described in terms of its complex tetrad components, (). In a Bondi coordinate/tetrad system the asymptotic integration is relatively simple [35, 25] resulting in the radial behavior (the peeling theorem):
where the leading coefficients of , () satisfy the evolution equations:The formal integration procedure is to take as an arbitrary function of () (the free broadcasting data), then integrate the second, for , with a timeindependent spinweight function of integration and finally integrate the first, for . Using a slight modification of this, namely from the spherical harmonic expansion, we obtain,
with the harmonic coefficients related to each other by the evolution equations: The physical meaning of the coefficients are etc. For later use, the complex dipole is written as . Note that the 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 (114).Later in this section it will be shown that we can find a unique complex world line, , (the world line associated with a shearfree NGC), that is closely related to the . From this complex world line we can define the intrinsic complex dipole moment, .Update
However, we first discuss a particular Maxwell field, , where one of its eigenvectors is a tangent field to a shearfree 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.
The parametric form of the general NGC was given earlier by Equation (72),
The parameters labeled the individual members of the congruence while 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, (or the Bondi coordinates) and the geodesic coordinates . 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 (74). We restrict ourselves to the special case of the coordinates and tetrad associated with the from a shearfree NGC. Though we are dealing with a real shearfree twisting congruence, the congruence, as we saw, is generated by a complex analytic world line in the complexified Minkowski space, . The complex parameter, , must in the end be chosen so that the ‘’ of Equation (90) is real. The Minkowski metric and the spin coefficients associated with this geodesic system can be calculated [27] in the frame. Unfortunately, it must be stated parametrically, since the explicitly appears via the and can not be directly eliminated. (An alternate choice of these geodesic coordinates is to use the instead of the . 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, , (an eigenvector of the Maxwell tensor, ), is a tangent vector of a shearfree 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 spincoefficient version of the Maxwell equations, using the geodesic tetrad, the choice of as the principle null vector ‘congruence’ is just the statement that

This allows a very simple exact integration of the remaining Maxwell components [35].
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 (), Equation (125). As the detailed calculations are long [35] and take us too far afield, we only give an outline here. The integration of the radial Maxwell equations leads to the asymptotic behavior,
with The order expression is known in terms of (). The function is given byRemark: In this case of the complex Liénard–Wiechert Maxwell field, the determines the intrinsic centerofcharge world line, rather than the relative centerofcharge line.
The remaining unknowns, , are determined by the last of the Maxwell equations,
which have been obtained from Equations (59) and (60) via the null rotation between the Bondi and geodesic tetrads and the associated Maxwell field transformation, namely, withThese remaining equations depend only on , which, in turn, is determined by . In other words, the solution is driven by the complex line, . As they now stand, Equations (130) appear to be difficult to solve, partially due to the implicit description of the .
Actually they are easily solved when the independent variables are changed, via Equation (86), from to the complex . They become, after a bit of work,
with the solution 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 . 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 (131) – (136), back to a Bondi coordinate/tetrad system. This can only be done by approximations. After a lengthy calculation [35], we find the Bondi peeling behavior
with Next, treating the world line, as discussed earlier, as a small deviation from the straight line, , i.e., byAgain to first order, Equations (142), (143) and (144) yield
the known electromagnetic dipole field, with a Coulomb charge, . One then has the physical interpretation of is the complex dipole moment; (the electric plus ‘’ times magnetic dipole) and 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, , and the intrinsic complex center of charge, , is true only at linear order. The secondorder 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 (14) and (134) to the geodesic system, it would have resulted in

This example was intended to show how physical meaning could be attached to the complex world line associated with a shearfree NGC. In this case and later in the case of asymptoticallyflat 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 flatspace asymptoticallyvanishing Maxwell field (with nonvanishing total charge) and for a vacuum asymptoticallyflat 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.
We return now to the general asymptoticallyvanishing Maxwell field, Equations (117) and (120), and its transformation behavior under the null rotation around ,
with , being one of our shearfree angle fields defined by a world line, . The leading components of the Maxwell fields transform asThe ‘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 shearfree congruence can be thought of as having its origin on the complex line, .) From the physical identifications of charge, dipole moments, etc., of Equation (123), we can obtain the transformation law of these physical quantities. In particular, the harmonic of , or, equivalently, the complex dipole, transforms as
where the notation means extract only the harmonic from a Clebsch–Gordon expansion of . 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 could be chosen so that . However, though this cannot be done in general, the harmonics of can be made to vanish by the appropriate choice of the . 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 function on , say, , given in a specific Bondi coordinate system, , that has a harmonic expansion given, for constant , by

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

the harmonic expansion at constant would be different. The new coefficients are extracted by the twosphere integral taken at constant :
It is in this rather obvious sense that the expansions are not unique. From Equation (97), using the real values of the complex cuts, i.e., (in place of , the transformation, Equation (153), and harmonic extraction yields . This is proportional to the transformed complex dipole moment. This transformation, a functional of the form, is decidedly nontensorial: in fact it is very nonlocal and nonlinear.Though it is clear that extracting 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.
The complex center of charge is defined by the vanishing of the complex dipole moment ; in other words,
determines the (up to this point) arbitrary complex world line, . In practice we do this only up to second order with the use of only the harmonics. The approximation we are using is to consider the charge as zeroth order and the dipole moments and the spatial part of the complex world line as first order.From Equation (153),
with the identifications, Equation (123), for and , we have to first order (with ), This is exactly the same result as we obtained earlier in Equation (116), via the charge and current distributions in a fixed Lorentz frame.Carrying this calculation [35] to second order, we find the secondorder relationship between the intrinsic complex dipole and the intrinsic dipole,
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,

with this uniquely determined world line is referred to as the Maxwell UCF.
http://www.livingreviews.org/lrr20096 
This work is licensed under a Creative Commons License. Problems/comments to 