Homogeneous Projective Coordinates for the Bondi-Metzner-Sachs Group

Giampiero Esposito $\,{}^{1,2,}$, Giuseppe Filiberto Vitale$\,{}^{1}$

$\,{}^{1}$Dipartimento di Fisica “Ettore Pancini”, Complesso Universitario di Monte S. Angelo, Via Cintia Edificio 6, 80126 Napoli, Italy

$\,{}^{2}$INFN Sezione di Napoli, Complesso Universitario di Monte S. Angelo, Via Cintia Edificio 6, 80126 Napoli, Italy

ABSTRACT This paper studies the Bondi-Metzner-Sachs group in homogeneous projective coordinates, because it is then possible to write all transformations of such a group in a manifestly linear way. The 2-sphere metric, Bondi-Metzner-Sachs metric, asymptotic Killing vectors, generators of supertranslations, as well as boosts and rotations of Minkowski spacetime, are all re-expressed in homogeneous projective coordinates. Last, the integral curves of vector fields which generate supertranslations are evaluated in detail. This work prepares the ground for more advanced applications of the differential geometry of asymptotically flat spacetimes in projective coordinates.

gesposito@na.infn.it, giuseppef.vitale@gmail.com

1 Introduction

The Bondi-Metzner-Sachs [1, 2, 3] asymptotic symmetry group of asymptotically flat spacetime has received again much attention over the last decade by virtue of its relevance for black-hole physics [4, 5, 6], the group-theoretical structure of general relativity [7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20] and the infrared structure of fundamental interactions [21, 22, 23, 24]. The appropriate geometric framework can be summarized as follows. In spacetime models for which null infinity can be defined, the cuts of null infinity are spacelike two-surfaces orthogonal to the generators of null infinity [25]. On using the familiar stereographic coordinate

$$\zeta=e^{i\varphi}\cot{\theta\over 2},$$ (1.1)

the first half of Bondi-Metzner-Sachs transformations read as

$$\zeta^{\prime}=f(\zeta)={(a\zeta+b)\over(c\zeta+d)}=f_{\Lambda}(\zeta),$$ (1.2)

where the matrix $\Lambda=\left(\begin{matrix}a&b\cr c&d\end{matrix}\right)$ has unit determinant $(ad-bc)=1$ and belongs therefore to the group ${\rm SL}(2,{\mathbb{C}})$. The resulting projective version of the special linear group can be defined as the space of pairs

$${\rm PSL}(2,{\mathbb{C}})=\left\{(f,\Lambda)|\;f:\zeta\in{\mathbb{C}}\rightarrow f_{\Lambda}(\zeta),\Lambda\in{\rm SL}(2,{\mathbb{C}})\right\},$$ (1.3)

i.e., the group of fractional linear maps $f_{\Lambda}$ according to Eq. (1.2) with the associated matrix $\Lambda$. Since

$$f_{\Lambda}(\zeta)={(a\zeta+b)\over(c\zeta+d)}={(-a\zeta-b)\over(-c\zeta-d)}=f_{-\Lambda}(\zeta),$$ (1.4)

one can write that ${\rm PSL}(2,{\mathbb{C}})$ is the quotient space ${\rm SL}(2,{\mathbb{C}})/\delta$, where $\delta$ is the homeomorphism defined by

$$\delta(a,b,c,d)=(-a,-b,-c,-d).$$ (1.5)

The fractional linear maps (1.2) can be defined for all values of $\zeta$ upon requiring that

$$f_{\Lambda}(\infty)={a\over c},\;f_{\Lambda}\left(-{d\over c}\right)=\infty.$$ (1.6)

Moreover, under fractional linear maps, lengths along the generators of null infinity scale according to

$$du^{\prime}=K_{\Lambda}(\zeta)du,$$ (1.7)

where the conformal factor is given by [25, 19]

$$K_{\Lambda}(\zeta)={{1+|\zeta|^{2}}\over{|a\zeta+b|^{2}+|c\zeta+d|^{2}}}.$$ (1.8)

By integration, Eq. (1.7) yields the second half of Bondi-Metzner-Sachs transformations:

$$u^{\prime}=K_{\Lambda}(\zeta)\Bigr{[}u+\alpha(\zeta,{\overline{\zeta}})\Bigr{]}.$$ (1.9)

As was pointed out in Ref. [19], the complex homogeneous coordinates associated to the Bondi-Metzner-Sachs transformation (1.2) have modulus $\leq 1$, which is the equation of a unit circle, and are

$$z_{0}=e^{i\varphi\over 2}\cos{\theta\over 2},\;z_{1}=e^{-i{\varphi\over 2}}\sin{\theta\over 2}.$$ (1.10)

In other words, upon remarking that

$$\zeta={z_{0}\over z_{1}},$$ (1.11)

Eq. (1.2) is equivalent to the linear transformation law

$$\left(\begin{matrix}z_{0}^{\prime}\cr z_{1}^{\prime}\end{matrix}\right)=\left(\begin{matrix}a&b\cr c&d\end{matrix}\right)\left(\begin{matrix}z_{0}\cr z_{1}\end{matrix}\right).$$ (1.12)

The next step of the program initiated in Ref. [19] consists in realizing that, much in the same way as the affine transformations in the Euclidean plane

$$x^{\prime}=x+a,\;y^{\prime}=y+b,$$ (1.13)

can be re-expressed with the help of a $3\times 3$ matrix in the form

$$\left(\begin{matrix}1&0&a\cr 0&1&b\cr 0&0&1\end{matrix}\right)\left(\begin{matrix}x\cr y\cr z\end{matrix}\right)=\left(\begin{matrix}x+a\cr y+b\cr 1\end{matrix}\right),$$ (1.14)

one can further re-express Eq. (1.12) with the help of a $3\times 3$ matrix in the form

$$\left(\begin{matrix}w_{0}^{\prime}\cr w_{1}^{\prime}\cr w_{2}^{\prime}\end{matrix}\right)=\left(\begin{matrix}1&0&0\cr 0&a&b\cr 0&c&d\end{matrix}\right)\left(\begin{matrix}w_{0}\cr w_{1}\cr w_{2}\end{matrix}\right),$$ (1.15)

with the understanding that Eq. (1.12) is the restriction to the unit circle $\Gamma$ of the map (1.15), upon defining

$$\left.w_{0}\right|_{\Gamma}=1,\;\left.w_{1}\right|_{\Gamma}=z_{0},\;\left.w_{2}\right|_{\Gamma}=z_{1}.$$ (1.16)

The work in Ref. [19] has outlined the resulting geometric picture, where $(w_{0},w_{1},w_{2})$ are viewed as homogeneous coordinates in a complex projective plane. In our paper we have instead a less abstract and more concrete task: since the Bondi-Metzner-Sachs transformation (1.2) becomes linear when expressed in terms of $z_{0}$ and $z_{1}$, we are aiming to develop the Bondi-Metzner-Sachs formalism with the associated Killing vector fields by using the pair of variables $(z_{0},z_{1})$ instead of $(\zeta,{\overline{\zeta}})$. For this purpose, the homogeneous projective coordinates for the 2-sphere are studied in Sect. 2, while the Bondi-Sachs metric in homogeneous coordinates is considered in Sect. 3. Asymptotic Killing fields for supertranslations are evaluated in Sect. 4, while their flow is investigated in Sect. 5. Concluding remarks and open problems are presented in Sect. 6, while technical details are provided in the Appendices.

2 Homogeneous coordinates on the $2$-sphere

It is useful, as an instrument to develop the BMS formalism in homogeneous coordinates, to re-write the 2-sphere metric in the desired coordinates. By using the definition (1.10), we get

$$z_{0}z_{1}=\sin(\theta/2)\cos(\theta/2)=\dfrac{\sin(\theta)}{2}\Rightarrow\theta=\sin^{-1}(2z_{0}z_{1})\text{,}$$ (2.1)

while for $\varphi$ we obtain

$$\dfrac{z_{0}}{z_{1}}=e^{i\varphi}\cot(\theta/2)\Rightarrow\varphi=-i\,\log\left(\tan(\theta/2)\dfrac{z_{0}}{z_{1}}\right)\text{.}$$ (2.2)

By virtue of the identity

$$\tan(\theta/2)==\dfrac{2\sin(\theta/2)\cos(\theta/2)}{2\cos^{2}(\theta/2)}=\dfrac{\sin(\theta)}{1+\sqrt{1-\sin^{2}(\theta)}},$$ (2.3)

we obtain for $\varphi$ the more convenient expression

$$\varphi=-i\,\log\left(\frac{2z^{2}_{0}}{1+\sqrt{1-4z^{2}_{0}z^{2}_{1}}}\right).$$

In order to re-express the $2$-sphere metric, let us evaluate

$$d\theta^{2} = \dfrac{d\theta}{dz_{0}}\dfrac{d\theta}{dz_{0}}dz_{0}^{2}+\dfrac{d\theta}{dz_{1}}\dfrac{d\theta}{dz_{1}}dz_{1}^{2}+2\dfrac{d\theta}{dz_{0}}\dfrac{d\theta}{dz_{1}}dz_{0}dz_{1}$$ (2.4)

$$= \frac{4z^{2}_{1}}{1-4z^{2}_{0}z^{2}_{1}}dz_{0}^{2}+\frac{4z^{2}_{0}}{1-4z^{2}_{0}z^{2}_{1}}dz_{1}^{2}+\frac{8z_{0}z_{1}}{1-4z^{2}_{0}z^{2}_{1}}dz_{0}dz_{1}\text{,}$$

while

$$\sin^{2}(\theta)\,d\varphi^{2} = 4z^{2}_{0}z^{2}_{1}\,d\varphi^{2}=4z^{2}_{0}z^{2}_{1}\left\{\dfrac{d\varphi}{dz_{0}}\dfrac{d\varphi}{dz_{0}}dz_{0}^{2}+\dfrac{d\varphi}{dz_{1}}\dfrac{d\varphi}{dz_{1}}dz_{1}^{2}+2\dfrac{d\varphi}{dz_{0}}\dfrac{d\varphi}{dz_{1}}dz_{0}dz_{1}\right\}$$ (2.5)

$$= \dfrac{-16z^{2}_{1}\left(1-2z^{2}_{0}z^{2}_{1}+\sqrt{1-4z^{2}_{0}z^{2}_{1}}\right)^{2}}{\left(1-4z^{2}_{0}z^{2}_{1}+\sqrt{1-4z^{2}_{0}z^{2}_{1}}\right)^{2}}dz_{0}^{2}-\dfrac{64z^{6}_{0}z^{4}_{1}}{\left(1-4z^{2}_{0}z^{2}_{1}+\sqrt{1-4z^{2}_{0}z^{2}_{1}}\right)^{2}}dz_{1}^{2}$$

$$- \dfrac{32z^{3}_{0}z^{3}_{1}}{1-4z^{2}_{0}z^{2}_{1}}dz_{0}dz_{1}.$$

Eventually, we obtain the metric for the 2-sphere in homogeneous coordinates

$$\Omega_{2} = d\theta^{2}+\sin^{2}(\theta)d\varphi^{2}=\sum_{\mu,\nu=0}^{1}g_{\mu\nu}dz^{\mu}dz^{\nu}$$ (2.6)

$$- 4z^{2}_{1}\left(\frac{1-4z^{2}_{0}z^{2}_{1}+2\sqrt{1-4z^{2}_{0}z^{2}_{1}}}{1-4z^{2}_{0}z^{2}_{1}}\right)dz_{0}^{2}+8z_{0}z_{1}dz_{0}dz_{1}$$

$$- 4z^{2}_{0}\left(\frac{1-4z^{2}_{0}z^{2}_{1}-2\sqrt{1-4z^{2}_{0}z^{2}_{1}}}{1-4z^{2}_{0}z^{2}_{1}}\right)dz_{1}^{2}.$$

At this stage, upon defining the real-valued function

$$\gamma(z_{0},z_{1})=\frac{2}{\sqrt{1-4z^{2}_{0}z^{2}_{1}}}=\frac{2}{\cos\theta},$$ (2.7)

we can write the matrix of metric components in the form

$$\gamma_{AB}=\begin{pmatrix}-4z^{2}_{1}\left(1+\gamma\right)&4z_{0}z_{1}\\ &\\ 4z_{0}z_{1}&-4z^{2}_{0}\left(1-\gamma\right)\end{pmatrix},$$ (2.8)

with non-vanishing determinant $-16z^{2}_{0}z^{2}_{1}\gamma^{2}$ and inverse matrix

$$\gamma^{AB}=\begin{pmatrix}\dfrac{1-\gamma}{4z^{2}_{1}\gamma^{2}}&\dfrac{1}{4z_{0}z_{1}\gamma^{2}}\\ &\\ \dfrac{1}{4z_{0}z_{1}\gamma^{2}}&\dfrac{1+\gamma}{4z^{2}_{0}\gamma^{2}}\end{pmatrix}.$$ (2.9)

We can see from (2.1) that the terms

$$2z_{0}z_{1}=\sin(\theta)\rightarrow 4z_{0}^{2}z^{2}_{1}=\sin^{2}(\theta)\rightarrow 1-4z^{2}_{0}z^{2}_{1}=\cos^{2}(\theta)\rightarrow 4z_{0}z_{1}=2\sin(\theta),$$

are real-valued, whereas

$$z^{2}_{0}=e^{i\varphi}\cos^{2}(\theta/2),\;z^{2}_{1}=e^{-i\varphi}\sin^{2}(\theta/2)\hskip 28.45274pt$$

are complex.

3 Bondi-Sachs metric in homogeneous coordinates

We can now write the retarded Bondi-Sachs (hereafter BS) metric in homogeneous coordinates with the help of the previous formulae. For this purpose, let us first write the general BS metric in the form

$$ds^{2}=-Udu^{2}-2e^{2\beta}dudr+h_{AB}\left(dx^{A}+\frac{1}{2}U^{A}du\right)\left(dx^{B}+\frac{1}{2}U^{B}du\right).$$ (3.1)

On passing from $(\theta,\varphi)$ to $(z_{0},z_{1})$ coordinates, we find the metric components of (3.1) expressed as follows:

$$g_{uu}=-U+{1\over 4}h_{z_{0}z_{0}}(U^{z_{0}})^{2}+{1\over 4}h_{z_{1}z_{1}}(U^{z_{1}})^{2}+{1\over 2}h_{z_{0}z_{1}}U^{z_{0}}U^{z_{1}},$$ (3.2)

$$g_{ur}=-e^{2\beta},$$ (3.3)

$$g_{uz_{0}}={1\over 2}(h_{z_{0}z_{0}}U^{z_{0}}+h_{z_{0}z_{1}}U^{z_{1}}),$$ (3.4)

$$g_{uz_{1}}={1\over 2}(h_{z_{0}z_{1}}U^{z_{0}}+h_{z_{1}z_{1}}U^{z_{1}}),$$ (3.5)

$$g_{z_{0}z_{0}}=h_{z_{0}z_{0}},\;g_{z_{0}z_{1}}=h_{z_{0}z_{1}},\;g_{z_{1}z_{1}}=h_{z_{1}z_{1}}.$$ (3.6)

The Bondi gauge $\partial_{r}\det\left(r^{-2}g_{AB}\right)=0$ implies that [26] $\gamma^{AB}C_{AB}=0$, where $\gamma^{AB}$ is given in Eq. (2.9). With our coordinates, this relation reads as

$$\gamma^{AB}C_{AB}=0\Leftrightarrow g^{z_{0}z_{0}}C_{z_{0}z_{0}}+g^{z_{1}z_{1}}C_{z_{1}z_{1}}+2g^{z_{0}z_{1}}C_{z_{0}z_{1}}=0.$$

We no longer have the simple result $C_{z\overline{z}}=0$ for the mixed component as in the stereographic coordinates, because in homogeneous coordinates we obtain

$$\dfrac{1-\gamma}{4z^{2}_{1}\gamma^{2}}C_{z_{0}z_{0}}+\dfrac{1+\gamma}{4z^{2}_{0}\gamma^{2}}C_{z_{1}z_{1}}+\dfrac{1}{2z_{0}z_{1}\gamma^{2}}C_{z_{0}z_{1}}=0,$$ (3.7)

which implies that

$$C_{z_{0}z_{1}}=-\dfrac{1}{2}\left(1-\gamma\right)\dfrac{z_{0}}{z_{1}}C_{z_{0}z_{0}}-\dfrac{1}{2}\left(1+\gamma\right)\dfrac{z_{1}}{z_{0}}C_{z_{1}z_{1}}.$$ (3.8)

The angular components of the metric are

$$g_{z_{0}z_{0}}=r^{2}\gamma_{z_{0}z_{0}}+rC_{z_{0}z_{0}}+{\mathcal{O}}(r),\;g_{z_{1}z_{1}}=r^{2}\gamma_{z_{1}z_{1}}+rC_{z_{1}z_{1}}+{\mathcal{O}}(r),$$

$$g_{z_{0}z_{1}}=r^{2}\gamma_{z_{0}z_{1}}+rC_{z_{0}z_{1}}+{\mathcal{O}}(r)$$

$$=r^{2}\gamma_{z_{0}z_{1}}-r\left(\dfrac{z_{0}}{z_{1}}\dfrac{\left(1-\gamma\right)}{2}C_{z_{0}z_{0}}+\dfrac{z_{1}}{z_{0}}\dfrac{\left(1+\gamma\right)}{2}C_{z_{1}z_{1}}\right)+{\mathcal{O}}(r),$$

where, of course, $\gamma_{AB}$ is given in Eq. (2.8). These formulae, jointly with the falloff conditions

$$U(u,r,x^{A})=1-\dfrac{2m(u,r,x^{A})}{r}+\dfrac{U_{2}(u,x^{A})}{r^{2}}+\mathcal{O}(r^{-3})$$

$$\beta(u,r,x^{A})=\dfrac{\beta_{1}(u,x^{A})}{r}+\dfrac{\beta_{2}(u,x^{A})}{r^{2}}+\mathcal{O}(r^{-3})$$

$$U^{A}(u,r,x^{B})=\dfrac{U^{A}_{2}(u,x^{B})}{r^{2}}+\dfrac{U^{A}_{3}(u,x^{B})}{r^{3}}+\mathcal{O}(r^{-4})$$

$$g_{AB}(u,r,x^{A})=r^{2}\gamma_{AB}(x^{A})+rC_{AB}(u,x^{A})+D_{AB}(u,x^{A})+\mathcal{O}(r^{-1}),$$ (3.9)

help to rewrite

$$g_{uu}=-\left(1-\dfrac{2m}{r}\right)+\mathcal{O}(r^{-2}).$$ (3.10)

Upon assuming that $\beta_{1}/r\ll 1$, we get

$$g_{ur}=-{\rm exp}\left(\dfrac{2\beta_{1}}{r}+\mathcal{O}(r^{-2})\right)=-1-\dfrac{2\beta_{1}}{r}+\mathcal{O}(r^{-2}),$$ (3.11)

while for $g_{uz_{0}}$ and $g_{uz_{1}}$ we find

$$g_{uz_{0}} =\dfrac{1}{2}\left(r^{2}\gamma_{z_{0}z_{0}}+rC_{z_{0}z_{0}}\right)\left(\dfrac{U^{z_{0}}_{2}}{r^{2}}+\dfrac{U^{z_{0}}_{3}}{r^{3}}\right)+\dfrac{1}{2}\left\{r^{2}\gamma_{z_{0}z_{1}}-r\left[\dfrac{z_{0}}{z_{1}}\dfrac{\left(1-\gamma\right)}{2}C_{z_{0}z_{0}}\right.\right.$$

$$\left.\left.+\dfrac{z_{1}}{z_{0}}\dfrac{\left(1+\gamma\right)}{2}C_{z_{1}z_{1}}\right]\right\}\left(\dfrac{U^{z_{1}}_{2}}{r^{2}}+\dfrac{U^{z_{1}}_{3}}{r^{3}}\right)$$

$$=\dfrac{\gamma_{z_{0}z_{0}}}{2}U^{z_{0}}_{2}+\dfrac{\gamma_{z_{0}z_{1}}}{2}U^{z_{1}}_{2}+\dfrac{1}{r}\left[\dfrac{\gamma_{z_{0}z_{0}}}{2}U^{z_{0}}_{3}+\dfrac{C_{z_{0}z_{0}}}{2}U^{z_{0}}_{2}+\dfrac{\gamma_{z_{0}z_{1}}}{2}U^{z_{1}}_{3}\right.$$

$$\left.-\dfrac{z_{0}}{z_{1}}\dfrac{\left(1-\gamma\right)}{4}C_{z_{0}z_{0}}U^{z_{1}}_{2}-\dfrac{z_{1}}{z_{0}}\dfrac{\left(1+\gamma\right)}{4}C_{z_{1}z_{1}}U^{z_{1}}_{2}\right]+\mathcal{O}(r^{-2})$$ (3.12)

and

$$g_{uz_{1}} =\dfrac{1}{2}\left(r^{2}\gamma_{z_{1}z_{1}}+rC_{z_{1}z_{1}}\right)\left(\dfrac{U^{z_{1}}_{2}}{r^{2}}+\dfrac{U^{z_{1}}_{3}}{r^{3}}\right)+\dfrac{1}{2}\left\{r^{2}\gamma_{z_{0}z_{1}}-r\left[\dfrac{z_{0}}{z_{1}}\dfrac{\left(1-\gamma\right)}{2}C_{z_{0}z_{0}}\right.\right.$$

$$\left.\left.+\dfrac{z_{1}}{z_{0}}\dfrac{\left(1+\gamma\right)}{2}C_{z_{1}z_{1}}\right]\right\}\left(\dfrac{U^{z_{0}}_{2}}{r^{2}}+\dfrac{U^{z_{0}}_{3}}{r^{3}}\right)$$

$$=\dfrac{\gamma_{z_{1}z_{1}}}{2}U^{z_{1}}_{2}+\dfrac{\gamma_{z_{0}z_{1}}}{2}U^{z_{0}}_{2}+\dfrac{1}{r}\left[\dfrac{\gamma_{z_{1}z_{1}}}{2}U^{z_{1}}_{3}+\dfrac{C_{z_{1}z_{1}}}{2}U^{z_{1}}_{2}+\dfrac{\gamma_{z_{0}z_{1}}}{2}U^{z_{0}}_{3}\right.$$

$$\left.-\dfrac{z_{0}}{z_{1}}\dfrac{\left(1-\gamma\right)}{4}C_{z_{0}z_{0}}U^{z_{0}}_{2}-\dfrac{z_{1}}{z_{0}}\dfrac{\left(1+\gamma\right)}{4}C_{z_{1}z_{1}}U^{z_{0}}_{2}\right]+\mathcal{O}(r^{-2}),$$ (3.13)

where use has been made of (3.8). Eventually, we get the matrix of Bondi metric components

$$g_{\mu\nu}=\begin{pmatrix}-\left(1-\dfrac{2m}{r}\right)&-1-\dfrac{2\beta_{1}}{r}&g_{uz_{0}}&g_{uz_{1}}\\ \\ -1-\dfrac{2\beta_{1}}{r}&0&0&0\\ \\ g_{uz_{0}}&0&r^{2}\gamma_{z_{0}z_{0}}+rC_{z_{0}z_{0}}&r^{2}\gamma_{z_{0}z_{1}}+rC_{z_{0}z_{1}}\\ \\ g_{uz_{1}}&0&r^{2}\gamma_{z_{0}z_{1}}+rC_{z_{0}z_{1}}&r^{2}\gamma_{z_{1}z_{1}}+rC_{z_{1}z_{1}}\end{pmatrix}\hskip 5.69046pt+\mathcal{O}(r^{-2}).$$ (3.14)

The gauge condition ${\rm det}\left(g_{AB}/r^{2}\right)=0$, instead of giving a solution for $D_{AB}$ such as in stereographic coordinates, gives us a condition for $C_{AB}$

$${\rm det}\left(g_{AB}\right)={\rm det}\begin{pmatrix}r^{2}\gamma_{z_{0}z_{0}}+rC_{z_{0}z_{0}}+D_{z_{0}z_{0}}&r^{2}\gamma_{z_{0}z_{1}}+rC_{z_{0}z_{1}}+D_{z_{0}z_{1}}\\ \\ r^{2}\gamma_{z_{0}z_{1}}+rC_{z_{0}z_{1}}+D_{z_{0}z_{1}}&r^{2}\gamma_{z_{1}z_{1}}+rC_{z_{1}z_{1}}+D_{z_{1}z_{1}}\end{pmatrix}$$

$$=r^{4}\left(\gamma_{z_{0}z_{0}}\gamma_{z_{1}z_{1}}-\gamma^{2}_{z_{0}z_{1}}\right)+r^{3}\left(\gamma_{z_{0}z_{0}}C_{z_{1}z_{1}}+\gamma_{z_{1}z_{1}}C_{z_{0}z_{0}}-2r^{3}\gamma_{z_{0}z_{1}}C_{z_{0}z_{1}}\right)$$

$$+r^{2}\left(\gamma_{z_{0}z_{0}}D_{z_{1}z_{1}}+C_{z_{0}z_{0}}C_{z_{1}z_{1}}+\gamma_{z_{1}z_{1}}D_{z_{0}z_{0}}-C^{2}_{z_{0}z_{1}}-2\gamma_{z_{0}z_{1}}D_{z_{0}z_{1}}\right)+\mathcal{O}(r),$$

$${\rm det}\left(\dfrac{g_{AB}}{r^{2}}\right)\Rightarrow\, \gamma_{z_{0}z_{0}}D_{z_{1}z_{1}}+C_{z_{0}z_{0}}C_{z_{1}z_{1}}+\gamma_{z_{1}z_{1}}D_{z_{0}z_{0}}-C^{2}_{z_{0}z_{1}}-2\gamma_{z_{0}z_{1}}D_{z_{0}z_{1}}\overset{!}{=}0$$

$$\Longrightarrow D_{z_{0}z_{1}}=\dfrac{\gamma_{z_{0}z_{0}}D_{z_{1}z_{1}}+C_{z_{0}z_{0}}C_{z_{1}z_{1}}+\gamma_{z_{1}z_{1}}D_{z_{0}z_{0}}-C^{2}_{z_{0}z_{1}}}{2\gamma_{z_{0}z_{1}}}$$

$$C^{2}_{z_{0}z_{1}}=C_{z_{0}z_{0}}C_{z_{1}z_{1}}.$$

In order to determine the various coefficients in the falloff conditions, we require that the Bondi metric should satisfy the Einstein equations

$$G_{\mu\nu}\equiv R_{\mu\nu}-\dfrac{1}{2}Rg_{\mu\nu}=8\pi GT_{\mu\nu}.$$

Upon restricting to the vacuum case $T=0$, in the limit as $r$ approaches $\infty$ in the Einstein tensor, first looking at $G_{rr}$, and neglecting the terms of order $\mathcal{O}(r^{-4})$, we get

$$G_{rr}=-\dfrac{4\beta_{1}}{r^{3}}+\mathcal{O}(r^{-4})\overset{!}{=}0\Rightarrow\beta_{1}\equiv 0.$$

Upon looking at $G_{rz_{0}}$ and $G_{rz_{1}}$ respectively, we get lengthy relations for $U^{z_{1}}_{2}$ and $U^{z_{0}}_{2}$, compared to the stereographic coordinates case, which depend on other coefficients. However, we still manage to solve directly for $U^{z_{0}}_{2}$ and $U^{z_{1}}_{2}$. On studying $G_{rA}=0$ we find

$$U^{z_{0}}_{2}=\dfrac{2z_{0}z_{1}\left(C_{z_{1}z_{1}}U^{z_{1}}_{2}+\gamma_{z_{0}z_{1}}U^{z_{0}}_{3}+\gamma_{z_{1}z_{1}}U^{z_{1}}_{3}\right)}{z^{2}_{1}\left(1+\gamma\right)C_{z_{1}z_{1}}+z^{2}_{0}\left(1-\gamma\right)C_{z_{0}z_{0}}}=-\dfrac{C_{z_{1}z_{1}}U^{z_{1}}_{2}+2\gamma_{z_{0}z_{1}}U^{z_{0}}_{3}}{C_{z_{0}z_{1}}},$$ (3.15)

and

$$U^{z_{1}}_{2}=\dfrac{2z_{0}z_{1}\left(C_{z_{0}z_{0}}U^{z_{0}}_{2}+\gamma_{z_{0}z_{0}}U^{z_{0}}_{3}+\gamma_{z_{0}z_{1}}U^{z_{1}}_{3}\right)}{z^{2}_{1}\left(1+\gamma\right)C_{z_{1}z_{1}}+z^{2}_{0}\left(1-\gamma\right)C_{z_{0}z_{0}}}=-\dfrac{C_{z_{0}z_{0}}U^{z_{0}}_{2}+2\gamma_{z_{0}z_{1}}U^{z_{1}}_{3}}{C_{z_{0}z_{1}}},$$ (3.16)

where we recall that $C_{z_{0}z_{1}}$ is given in Eq. (3.8). By virtue of Eqs. (3.12) and (3.13) we find eventually the metric in the form

$$ds^{2}=-du^{2}-2dudr+2\left(r^{2}\gamma_{z_{0}z_{1}}+rC_{z_{0}z_{1}}\right)dz_{0}dz_{1}+\dfrac{2m}{r}du^{2}+\left(r^{2}\gamma_{z_{0}z_{0}}+rC_{z_{0}z_{0}}\right)dz^{2}_{0}$$ (3.17)

$$+\left(r^{2}\gamma_{z_{1}z_{1}}+rC_{z_{1}z_{1}}\right)dz^{2}_{1}$$

$$+\biggr{[}{\gamma_{z_{0}z_{0}}\over 2}U_{2}^{z_{0}}+{\gamma_{z_{0}z_{1}}\over 2}U_{2}^{z_{1}}+{1\over r}\biggr{(}{C_{z_{0}z_{0}}\over 2}U_{2}^{z_{0}}+{C_{z_{0}z_{1}}\over 2}U_{2}^{z_{1}}+{\gamma_{z_{0}z_{0}}\over 2}U_{3}^{z_{0}}+{\gamma_{z_{0}z_{1}}\over 2}U_{3}^{z_{1}}\biggr{)}\biggr{]}dudz_{0}$$

$$+\biggr{[}{\gamma_{z_{0}z_{1}}\over 2}U_{2}^{z_{0}}+{\gamma_{z_{1}z_{1}}\over 2}U_{2}^{z_{1}}+{1\over r}\biggr{(}{C_{z_{0}z_{1}}\over 2}U_{2}^{z_{0}}+{C_{z_{1}z_{1}}\over 2}U_{2}^{z_{1}}+{\gamma_{z_{0}z_{1}}\over 2}U_{3}^{z_{0}}+{\gamma_{z_{1}z_{1}}\over 2}U_{3}^{z_{1}}\biggr{)}\biggr{]}dudz_{1}$$

$$+\mathcal{O}(r^{-2}).$$

Now we are ready to evaluate the BMS generators in homogeneous coordinates in order to determine the supertranslations.

4 Asymptotic Killing fields

After finding the most general Bondi metric in homogeneous coordinates satisfying the asymptotically flat spacetime falloffs, our aim is to find the most general vector fields $\xi$ satisfying the Bondi gauge condition and the asymptotically flat spacetime falloffs. As is well known, the Killing vectors solve by definition the equations

$$\left(\mathcal{L}_{\xi}g\right)_{\mu\nu}=\xi^{\rho}\partial_{\rho}g_{\mu\nu}+g_{\mu\rho}\partial_{\nu}\xi^{\rho}+g_{\nu\rho}\partial_{\mu}\xi^{\rho}=0.$$

Moreover, the preservation of the Bondi gauge condition yields [26]

$$\left(\mathcal{L}_{\xi}g\right)_{rr}=0\text{,}\hskip 5.69046pt\left(\mathcal{L}_{\xi}g\right)_{rA}=0\hskip 5.69046pt\text{and}\hskip 5.69046ptg^{AB}\left(\mathcal{L}_{\xi}g\right)_{AB}=0.$$ (4.1)

From these relations one can calculate the four components of $\xi^{\mu}$. At this stage, instead of repeating the detailed calculations already available, for example, in Ref. [26], we can compute the asymptotic Killing fields in homogeneous coordinates by using the familiar transformation law of vector fields. In other words, the work in Ref. [26] has defined the stereographic variable (we write $\psi$ rather than $z$ used in Ref. [26], in order to avoid confusion with our $\zeta$ in Eq. (1.1))

$$\psi=e^{i\varphi}\tan{\theta\over 2}={1\over{\bar{\zeta}}},$$ (4.2)

and has found, in Bondi coordinates $u,r,\theta,\varphi$, the asymptotic Killing fields $\xi_{T}^{+}$ where the components depend on a function $f$ and on the Bondi coordinates. On denoting as usual by $Y_{l}^{m}$ the spherical harmonics on the $2$-sphere, one finds [26]

$$\xi_{T}^{+}\biggr{|}_{f=Y_{0}^{0}}={\partial\over\partial u},$$ (4.3)

$$\xi_{T}^{+}\biggr{|}_{f=Y_{1}^{0}}={(1-\psi{\bar{\psi}})\over(1+\psi{\bar{\psi}})}\left({\partial\over\partial u}-{\partial\over\partial r}\right)+{\psi\over r}{\partial\over\partial\psi}+{{\bar{\psi}}\over r}{\partial\over\partial{\bar{\psi}}},$$ (4.4)

$$\xi_{T}^{+}\biggr{|}_{f=Y_{1}^{1}}={\psi\over(1+\psi{\bar{\psi}})}\left({\partial\over\partial u}-{\partial\over\partial r}\right)+{\psi^{2}\over 2r}{\partial\over\partial\psi}-{1\over 2r}{\partial\over\partial{\bar{\psi}}},$$ (4.5)

$$\xi_{T}^{+}|_{f=Y_{1}^{-1}}={{\bar{\psi}}\over(1+\psi{\bar{\psi}})}\left({\partial\over\partial u}-{\partial\over\partial r}\right)-{1\over 2r}{\partial\over\partial\psi}+{{\bar{\psi}}^{2}\over 2r}{\partial\over\partial{\bar{\psi}}}.$$ (4.6)

Now by virtue of the basic identities

$${\partial\over\partial\psi}={\partial z_{0}\over\partial\psi}{\partial\over\partial z_{0}}+{\partial z_{1}\over\partial\psi}{\partial\over\partial z_{1}},$$ (4.7)

$${\partial\over\partial{\bar{\psi}}}={\partial z_{0}\over\partial{\bar{\psi}}}{\partial\over\partial z_{0}}+{\partial z_{1}\over\partial{\bar{\psi}}}{\partial\over\partial z_{1}},$$ (4.8)

and upon exploiting the formulae (A7)-(A10) in the Appendix, we find

$$\xi_{T}^{+}\biggr{|}_{f=Y_{1}^{0}}={2\over\gamma}\left({\partial\over\partial u}-{\partial\over\partial r}\right)+{z_{0}\over 2r}{(2-\gamma)\over\gamma}{\partial\over\partial z_{0}}+{z_{1}\over 2r}{(2+\gamma)\over\gamma}{\partial\over\partial z_{1}},$$ (4.9)

$$\xi_{T}^{+}\biggr{|}_{f=Y_{1}^{1}} = {z_{0}\over 2z_{1}}{(\gamma-2)\over\gamma}\left({\partial\over\partial u}-{\partial\over\partial r}\right)+{1\over r}{(z_{0})^{2}\over z_{1}}\left({1\over 4}-{1\over\gamma(\gamma+2)}\right){\partial\over\partial z_{0}}$$ (4.10)

$$+ {z_{0}\over 2r}\left({1\over(\gamma+2)}-{(\gamma+2)\over 2\gamma}\right){\partial\over\partial z_{1}},$$

$$\xi_{T}^{+}\biggr{|}_{f=Y_{1}^{-1}} = {z_{1}\over 2z_{0}}{(\gamma+2)\over\gamma}\left({\partial\over\partial u}-{\partial\over\partial r}\right)-{z_{1}\over 2r}\left({1\over(\gamma-2)}+{(\gamma-2)\over 2\gamma}\right){\partial\over\partial z_{0}}$$ (4.11)

$$+ {1\over r}{(z_{1})^{2}\over z_{0}}\left({1\over 4}-{1\over\gamma(\gamma-2)}\right){\partial\over\partial z_{1}}.$$

Now we denote by $\xi_{0},\xi_{1},\xi_{2},\xi_{3}$ the vector fields (4.3), (4.9), (4.10) and (4.11), respectively. Nontrivial Lie brackets among them involve $\xi_{1},\xi_{2},\xi_{3}$ only. With our notation, we can re-write Eqs. (4.9)-(4.11) in the form

$$\xi_{1}=A_{11}\left({\partial\over\partial u}-{\partial\over\partial r}\right)+A_{12}{\partial\over\partial z_{0}}+A_{13}{\partial\over\partial z_{1}},$$ (4.12)

$$\xi_{2}=A_{21}\left({\partial\over\partial u}-{\partial\over\partial r}\right)+A_{22}{\partial\over\partial z_{0}}+A_{23}{\partial\over\partial z_{1}},$$ (4.13)

$$\xi_{3}=A_{31}\left({\partial\over\partial u}-{\partial\over\partial r}\right)+A_{32}{\partial\over\partial z_{0}}+A_{33}{\partial\over\partial z_{1}},$$ (4.14)

where the values taken by the $A_{ij}$ functions can be read off from (4.9)-(4.11). At this stage, a patient evaluation proves that such vector fields have vanishing Lie brackets:

$$[\xi_{1},\xi_{2}]=[\xi_{2},\xi_{3}]=[\xi_{3},\xi_{1}]=0.$$ (4.15)

The result is simple, but the actual proof requires several details, for which we refer the reader to Appendix B.

Figure 1: Numerical evaluation of the integral curve for the supertranslation vector field (4.9). The initial conditions (5.14) are taken to be $u=0,r=1,z_{0}=e^{i{\pi\over 8}}\cos{\pi\over 8},z_{1}=e^{-i{\pi\over 8}}\sin{\pi\over 8}$.

Figure 2: Numerical evaluation of the integral curve for the supertranslation vector field (4.10). The initial conditions (5.14) are taken to be $u=0,r=1,z_{0}=e^{i{\pi\over 8}}\cos{\pi\over 8},z_{1}=e^{-i{\pi\over 8}}\sin{\pi\over 8}$.

Figure 3: Numerical evaluation of the integral curve for the supertranslation vector field (4.11). The initial conditions (5.14) are taken to be $u=0,r=1,z_{0}=e^{i{\pi\over 8}}\cos{\pi\over 8},z_{1}=e^{-i{\pi\over 8}}\sin{\pi\over 8}$.

Figure 4: Numerical evaluation of the integral curve for the supertranslation vector field (4.9). The initial conditions (5.14) are taken to be $u=1,r=1,z_{0}=e^{i{\pi\over 8}}\cos{\pi\over 4},z_{1}=e^{-i{\pi\over 8}}\sin{\pi\over 4}$.

Figure 5: Numerical evaluation of the integral curve for the supertranslation vector field (4.10). The initial conditions (5.14) are taken to be $u=1,r=1,z_{0}=e^{i{\pi\over 8}}\cos{\pi\over 4},z_{1}=e^{-i{\pi\over 8}}\sin{\pi\over 4}$.

Figure 6: Numerical evaluation of the integral curve for the supertranslation vector field (4.11). The initial conditions (5.14) are taken to be $u=1,r=1,z_{0}=e^{i{\pi\over 8}}\cos{\pi\over 4},z_{1}=e^{-i{\pi\over 8}}\sin{\pi\over 4}$.

Figure 7: Numerical evaluation of the integral curve for the supertranslation vector field (4.9). The initial conditions (5.14) are taken to be $u=0,r=1,z_{0}=e^{i{\pi\over 4}}\cos{\pi\over 12},z_{1}=e^{-i{\pi\over 4}}\sin{\pi\over 12}$. In this particular case, the real parts meet at a single point.

Figure 8: Numerical evaluation of the integral curve for the supertranslation vector field (4.10). The initial conditions (5.14) are taken to be $u=0,r=1,z_{0}=e^{i{\pi\over 4}}\cos{\pi\over 12},z_{1}=e^{-i{\pi\over 4}}\sin{\pi\over 12}$.

Figure 9: Numerical evaluation of the integral curve for the supertranslation vector field (4.11). The initial conditions (5.14) are taken to be $u=0,r=1,z_{0}=e^{i{\pi\over 4}}\cos{\pi\over 12},z_{1}=e^{-i{\pi\over 4}}\sin{\pi\over 12}$.

5 Flow of supertranslation vector fields

In order to appreciate that the familiar geometric constructions are feasible also in projective coordinates, we now consider the flow of supertranslation vector fields (4.9)-(4.11). For example, by virtue of (2.7), and defining $p=(u,r,z_{0},z_{1})$, the task of finding the flow of the supertranslation vector fields (4.9), (4.10) and (4.11) consists of solving a system of nonlinear and coupled differential equations. For this purpose, we denote by $\sigma,\Sigma,\chi$, respectively, the appropriate flow, and define

$$\delta(W;\tau,p)=\sqrt{1-4\Bigr{(}W^{3}(\tau,p)W^{4}(\tau,p)\Bigr{)}^{2}},$$ (5.1)

where $W=\sigma,\Sigma,\chi$, respectively, with components $W^{1},W^{2},W^{3},W^{4}$. Hence we study the following coupled systems of nonlinear differential equations:

$${d\sigma^{1}\over d\tau}=\delta(\sigma;\tau,p),$$ (5.2)

$${d\sigma^{2}\over d\tau}=-\delta(\sigma;\tau,p),$$ (5.3)

$${d\sigma^{3}\over d\tau}={\sigma^{3}(\tau,p)\over 2\sigma^{2}(\tau,p)}(\delta(\sigma;\tau,p)-1),$$ (5.4)

$${d\sigma^{4}\over d\tau}={\sigma^{4}(\tau,p)\over 2\sigma^{2}(\tau,p)}(\delta(\sigma;\tau,p)+1),$$ (5.5)

$${d\Sigma^{1}\over d\tau}={\Sigma^{3}(\tau,p)\over 2\Sigma^{4}(\tau,p)}(1-\delta(\Sigma;\tau,p)),$$ (5.6)

$${d\Sigma^{2}\over d\tau}=-{\Sigma^{3}(\tau,p)\over 2\Sigma^{4}(\tau,p)}(1-\delta(\Sigma;\tau,p)),$$ (5.7)

$${d\Sigma^{3}\over d\tau}={(\Sigma^{3}(\tau,p))^{2}\over 4\Sigma^{2}(\tau,p)\Sigma^{4}(\tau,p)}\left[1-\delta(\Sigma;\tau,p)+{\delta(\Sigma;\tau,p)\over(1+\delta(\Sigma;\tau,p))}\right],$$ (5.8)

$${d\Sigma^{4}\over d\tau}={\Sigma^{3}(\tau,p)\over 4\Sigma^{2}(\tau,p)}\left[{\delta(\Sigma;\tau,p)\over(1+\delta(\Sigma;\tau,p))}-1-\delta(\Sigma;\tau,p)\right],$$ (5.9)

$${d\chi^{1}\over d\tau}={\chi^{4}(\tau,p)\over 2\chi^{3}(\tau,p)}(1+\delta(\chi;\tau,p)),$$ (5.10)

$${d\chi^{2}\over d\tau}=-{\chi^{4}(\tau,p)\over 2\chi^{3}(\tau,p)}(1+\delta(\chi;\tau,p)),$$ (5.11)

$${d\chi^{3}\over d\tau}=-{\chi^{4}(\tau,p)\over 4\chi^{2}(\tau,p)}\left[{\delta(\chi;\tau,p)\over(1-\delta(\chi;\tau,p)}+1-\delta(\chi;\tau,p)\right],$$ (5.12)

$${d\chi^{4}\over d\tau}={(\chi^{4}(\tau,p))^{2}\over 4\chi^{2}(\tau,p)\chi^{3}(\tau,p)}\left[1+\delta(\chi;\tau,p)-{\delta(\chi;\tau,p)\over(1-\delta(\chi;\tau,p))}\right],$$ (5.13)

with the initial conditions

$$W^{1}(0,p)=u,\;W^{2}(0,p)=r,\;W^{3}(0,p)=z_{0},\;W^{4}(0,p)=z_{1}.$$ (5.14)

The resulting equations can only be solved numerically, to the best of our knowledge, and such solutions are displayed in Figures from 1 to 9. Since the desired solutions are complex-valued, we have displayed both real and imaginary parts, with three choices of initial conditions.

6 Concluding remarks and open problems

As far as we can see, the interest of our investigation lies in having shown that homogeneous projective coordinates lead to a fully computational scheme for all applications of the BMS group. This might pay off when more advanced properties will be studied. In particular, we have in mind the concept of superrotations [21, 22, 26] on the one hand, and the physical applications of the Segre manifold advocated in Ref. [19]. In other words, since our Eq. (1.15) contains Eq. (1.12), which in turn is just a re-expression of the BMS transformation (1.2), one might aim at embedding the study of BMS symmetries into the richer mathematical framework of complex analysis [27] and algebraic geometry. Our paper has tried to prepare the ground for such a synthesis.