Comparing a Gauge-Invariant Formulation and a “Conventional Complete Gauge-Fixing Approach” for l = 0, 1 Mode Perturbations on the Schwarzschild Background Spacetime

1. Introduction

Through the ground-based gravitational-wave detectors [1,2,3,4], many events of gravitational waves, mainly from black hole-black hole coalescences, have now been detected. We are now at the stage where there is no doubt about the existence of gravitational waves due to their direct observations. One of the future directions of gravitational-wave astronomy will be a precise science through the statistics of many events. Toward further development of gravitational-wave science, the projects of future ground-based gravitational-wave detectors [5,6] are also progressing to achieve more sensitive detectors and some projects of space gravitational-wave antenna are also progressing [7,8,9,10]. Although there are many targets of these detectors, the Extreme-Mass-Ratio-Inspiral (EMRI), which is a source of gravitational waves from the motion of a stellar mass object around a supermassive black hole, is a promising target of the Laser Interferometer Space Antenna [7]. Since the mass ratio of this EMRI is very small, we can describe the gravitational waves from EMRIs through black hole perturbations [11]. Furthermore, the sophistication of higher-order black hole perturbation theories is required to support gravitational-wave physics as a precise science. The motivation of our series of papers, Refs. [12,13,14,15,16] and this paper, is in this theoretical sophistication of black hole perturbation theories toward higher-order perturbations.

Although realistic black holes have their angular momentum and we have to consider the perturbation theory of a Kerr black hole for direct applications to the EMRI, we may say that further sophistications are possible even in perturbation theories on the Schwarzschild background spacetime. From the pioneering works by Regge and Wheeler [17] and Zerilli [18,19], there have been many studies on the perturbations in the Schwarzschild background spacetime [20,21,22,23,24,25,26,27,28,29,30,31,32,33]. In these works, perturbations are decomposed through the spherical harmonics $Y_{lm}$ because of the spherical symmetry of the background spacetime, and $l=0,1$ modes should be separately treated. These modes correspond to the monopole and dipole perturbations. Due to these separate treatments, “gauge-invariant” treatments for $l=0$ and $l=1$ modes were unclear.

Owing to this situation, in the previous papers [12,13,14,15,16], we proposed the strategy of the gauge-invariant treatments of these $l=0,1$ mode perturbations, which is declared as Proposal 2.1 in Sec. 2 of this paper below. One of the important premises of our gauge-invariant perturbations is the distinction between the first-kind gauge and the second-kind gauge. The first-kind gauge is essentially the choice of the coordinate system on the single manifold, and we often use this first-kind gauge when we predict or interpret the measurement results of experiments and observations. On the other hand, the second-kind gauge is the choice of the point identifications between the points on the physical spacetime ${\mathcal{M}}_{ϵ}$ and the background spacetime $\mathcal{M}$. This second-kind gauge has nothing to do with our physical spacetime ${\mathcal{M}}_{ϵ}$. Although this difference is extensively explained in the Part I paper [14], we also emphasize this difference in Sec. 2 of this paper. The proposal in the Part I paper [14] is a part of our developments of the general formulation of a higher-order gauge-invariant perturbation theory on a generic background spacetime toward unambiguous sophisticated nonlinear general-relativistic perturbation theories [34,35,36,37,38,39]. Although we have been applied this general framework to cosmological perturbations [40,41,42,43,44,45,46,47], we applied it to black hole perturbations in the series of papers, i.e., Refs. [12,13,14,15,16] and this paper. Even in cosmological perturbation theories, the same problem as the above $l=0,1$-mode problem exists as gauge-invariant treatments of homogeneous modes of perturbations. In this sense, we can expect that the proposal in the previous paper [14] will be a clue to the same problem in gauge-invariant perturbation theory on the generic background spacetime.

In the Part I paper [14], we also derived the linearized Einstein equations in a gauge-invariant manner following Proposal 2.1. Perturbations on the spherically symmetric background spacetime are classified into even- and odd-mode perturbations. In the same paper [14], we also gave the strategy to solve the odd-mode perturbations including $l=0,1$ modes. Furthermore, we also derived the formal solutions for the $l=0,1$ odd-mode perturbations to the linearized Einstein equations following Proposal 2.1. In the Part II paper [15], we gave the strategy to solve the even-mode perturbations including $l=0,1$ modes, and we also derive the formal solutions for the $l=0,1$ even-mode perturbations following Proposal 2.1. In the Part III paper [16], we find the fact that the derived solutions in the Part II paper [15] realize the linearized version of two exact solutions: one is the Lemaître-Tolman-Bondi (LTB) solution [48] and the other is the non-rotating C-metric [49,50]. Due to this fact, we conclude that the solutions for even-mode perturbations derived in the Part II paper [15] are physically reasonable. This series of papers is the full paper version of our short paper [12]. Furthermore, brief discussions on the extension to the higher-order perturbations are given in the short paper [13].

On the other hand, it is well-known the fact that we cannot construct gauge-invariant variables for $l=0,1$ modes in a similar manner to $l\ge 2$ modes if we decompose the metric perturbations through the spherical harmonics as the scalar harmonics from the starting point. For this reason, we usually use gauge-fixing approaches. Furthermore, it is often said that “gauge-invariant formulations in general-relativistic perturbations are equivalent to complete gauge-fixing approaches.” In this paper, we check this statement through the comparison of our proposed gauge-invariant formulation, in which we introduce singular harmonics at once and regularize them after the derivation of the mode-by-mode Einstein equation, and a “conventional complete gauge-fixing approach”, in which we use the spherical harmonic functions $Y_{lm}$ as the scalar harmonics $S_{\delta }$ from the starting point. As the result of this comparison, we conclude that our gauge-invariant formulation and the above “conventional complete gauge-fixing approach” are different, though these two formulations lead similar solutions for $l=0,1$-mode perturbations. More specifically, there is a case where the boundary conditions and initial conditions are restricted in a “conventional complete gauge-fixing approach” where we use the decomposition of the metric perturbation by the spherical harmonics $Y_{lm}$ from the starting point.

The organization of this paper is as follows. In Sec. 2, we briefly review the premise of our series of papers [12,13,14,15,16], which are necessary for the ingredients of this paper. We also emphasize the difference between the concepts of the first- and the second-kind gauges and summarize the linearized Einstein equations for $l=0,1$ modes and their formal solutions in Sec. 2. In Sec. 3, we specify the rule of our comparison between our gauge-invariant formulation and a conventional gauge-fixing approach and summarize the gauge transformation rules for the metric perturbations of $l=0,1$ modes, because the above statement “gauge-invariant formulations in general-relativistic perturbations are equivalent to complete gauge-fixing approaches” includes some ambiguous. In Sec. 4, we discuss the linearized Einstein equations for $l=1$ odd-modes and their solutions in the conventional gauge-fixing approach. In Sec. 5, we discuss the linearized Einstein equations for $l=1$ even-modes and their solutions in the conventional gauge-fixing approach. In Sec. 6, we derive the solution to the linearized Einstein equations for $l=0$ modes through the complete gauge-fixing and discuss the comparison with linearized LTB solution. Sec. 7 is devoted to the summary of this paper and discussions based on our results.

We use the notation used in the previous papers [12,13,14,15,16] and the unit $G=c=1$, where G is Newton’s constant of gravitation and c is the velocity of light.

2. Brief Review of a Gauge-Invariant Treatment of $l=0,1$ Modes

In this section, we briefly review the premises of our series of papers [12,13,14,15,16] which are necessary for the ingredients of this paper. In Sec. 2.1, we review our framework of the gauge-invariant perturbation theory [34,35]. This is an important premise of the series of our papers [12,13,14,15,16] and this paper. In Sec. 2.2, we review the gauge-invariant perturbation theory on spherically symmetric spacetimes which includes our proposal in Refs. [12,13,14,15,16]. In Sec. 2.3, we summarize the $l=1$ odd-mode linearized Einstein equations and their “formal” solutions. In Sec. 2.4, we summarize the $l=0,1$ even-mode linearized Einstein equations and their “formal” solutions. The equations and their solutions in Sec. 2.3 and Sec. 2.4 are derived based on our proposal in Refs. [12,13,14,15,16]. These are necessary for the arguments within this paper.

2.1. General Framework of Gauge-Invariant Perturbation Theory

In any perturbation theory, we always treat two spacetime manifolds. One is the physical spacetime $({\mathcal{M}}_{\mathrm{ph}},{\overline{g}}_{ab})$, which is identified with our nature itself, and we want to describe this spacetime $({\mathcal{M}}_{\mathrm{ph}},{\overline{g}}_{ab})$ by perturbations. The other is the background spacetime $(\mathcal{M},g_{ab})$, which is prepared as a reference by hand. Note that these two spacetimes are distinct. Furthermore, in any perturbation theory, we always write equations for the perturbation of the variable Q as follows:

$$Q(“p”)=Q_0\left(p\right)+\delta Q\left(p\right).$$

(1)

Equation (1) gives a relation between variables on different manifolds. Actually, $Q(“p”)$ in Equation (1) is a variable on the physical spacetime ${\mathcal{M}}_{ϵ}={\mathcal{M}}_{\mathrm{ph}}$, whereas $Q_0\left(p\right)$ and $\delta Q\left(p\right)$ are variables on the background spacetime $\mathcal{M}$. Because we regard Equation (1) as a field equation, Equation (1) includes an implicit assumption of the existence of a point identification map ${\mathcal{X}}_{ϵ}$: $\mathcal{M}\to {\mathcal{M}}_{ϵ}$: $p\in \mathcal{M}\mapsto “p”\in {\mathcal{M}}_{ϵ}$. This identification map is a gauge choice in general-relativistic perturbation theories. This is the notion of the second-kind gauge pointed out by Sachs [51,52,53,54]. Note that this second-kind gauge is a different notion from the degree of freedom of the coordinate transformation on a single manifold, which is called the first-kind gauge [14,40,41]. This distinction between the first- and the second-kind of gauges extensively explained in the Part I paper [14] and is also important to understand the ingredients of this paper.

To compare the variable Q on ${\mathcal{M}}_{ϵ}$ with its background value $Q_0$ on $\mathcal{M}$, we use the pull-back ${\mathcal{X}}_{ϵ}^{*}$ of the identification map ${\mathcal{X}}_{ϵ}$:$\mathcal{M}\to {\mathcal{M}}_{ϵ}$ and we evaluate the pulled-back variable ${\mathcal{X}}_{ϵ}^{*}Q$ on the background spacetime $\mathcal{M}$. Furthermore, in perturbation theories, we expand the pull-back operation ${\mathcal{X}}_{ϵ}^{*}$ to the variable Q with respect to the infinitesimal parameter $ϵ$ for the perturbation as

$$\begin{array}{c}\hfill {\mathcal{X}}_{ϵ}^{*}Q=Q_0+ϵ{}_{\mathcal{X}}{}^{\left(1\right)}Q+O\left({ϵ}^2\right).\end{array}$$

(2)

Equation (2) are evaluated on the background spacetime $\mathcal{M}$. When we have two different gauge choices ${\mathcal{X}}_{ϵ}$ and ${\mathcal{Y}}_{ϵ}$, we can consider the gauge transformation, which is the change of the point-identification ${\mathcal{X}}_{ϵ}\to {\mathcal{Y}}_{ϵ}$. This gauge transformation is given by the diffeomorphism ${\Phi }_{ϵ}$$:=$ ${\left({\mathcal{X}}_{ϵ}\right)}^{-1}\circ {\mathcal{Y}}_{ϵ}$ :$\mathcal{M}$→$\mathcal{M}$. Actually, the diffeomorphism ${\Phi }_{ϵ}$ induces a pull-back from the representation ${\mathcal{X}}_{ϵ}^{*}{Q}_{ϵ}$ to the representation ${\mathcal{Y}}_{ϵ}^{*}{Q}_{ϵ}$ as ${\mathcal{Y}}_{ϵ}^{*}{Q}_{ϵ}\left(q\right)={\Phi }_{ϵ}^{*}{\mathcal{X}}_{ϵ}^{*}{Q}_{ϵ}\left(q\right)$ at any point $q\in \mathcal{M}$. From general arguments of the Taylor expansion [55,56,57], the pull-back ${\Phi }_{ϵ}^{*}$ is expanded as

$$\begin{array}{ccc}\hfill {\mathcal{Y}}_{ϵ}^{*}{Q}_{ϵ}\left(q\right)& =& {\mathcal{X}}_{ϵ}^{*}{Q}_{ϵ}\left(q\right)+ϵ{\pounds }_{{\xi }_{\left(1\right)}}{\mathcal{X}}_{ϵ}^{*}{Q}_{ϵ}\left(q\right)+O\left({ϵ}^2\right),\hfill \end{array}$$

(3)

where ${\xi }_{\left(1\right)}^a$ is the generator of ${\Phi }_{ϵ}$. From Equations (2) and (3), the gauge transformation for the first-order perturbation ${}^{\left(1\right)}Q$ is given by

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}{}^{\left(1\right)}Q\left(q\right)-{}_{\mathcal{X}}{}^{\left(1\right)}Q\left(q\right)& =& {\pounds }_{{\xi }_{\left(1\right)}}{Q}_0\left(q\right)\hfill \end{array}$$

(4)

at any point $q\in \mathcal{M}$. We also employ the order by order gauge invariance as a concept of gauge invariance [34,35]. We call the kth-order perturbation ${}_{\mathcal{X}}{}^{\left(k\right)}Q$ as gauge invariant if and only if

$$\begin{array}{c}\hfill {}_{\mathcal{X}}{}^{\left(k\right)}Q\left(q\right)={}_{\mathcal{Y}}{}^{\left(k\right)}Q\left(q\right)\end{array}$$

(5)

for any gauge choice ${\mathcal{X}}_{ϵ}$ and ${\mathcal{Y}}_{ϵ}$ at any point of $q\in \mathcal{M}$.

Based on the above setup, we proposed a formulation to construct gauge-invariant variables of higher-order perturbations [34,35]. First, we expand the metric on the physical spacetime ${\mathcal{M}}_{ϵ}$, which was pulled back to the background spacetime $\mathcal{M}$ through a gauge choice ${\mathcal{X}}_{ϵ}$ as

$$\begin{array}{ccc}\hfill {\mathcal{X}}_{ϵ}^{*}{\overline{g}}_{ab}& =& g_{ab}+ϵ{}_{\mathcal{X}}{h}_{ab}+O\left({ϵ}^2\right).\hfill \end{array}$$

(6)

Although the expression (6) depends entirely on the gauge choice ${\mathcal{X}}_{ϵ}$, henceforth, we do not explicitly express the index of the gauge choice ${\mathcal{X}}_{ϵ}$ in the expression if there is no possibility of confusion. The important premise of our formulation of higher-order gauge-invariant perturbation theory was the following conjecture [34,35] for the linear metric perturbation $h_{ab}$:

Conjecture 2.1.

If the gauge transformation rule for a perturbative pulled-back tensor field $h_{ab}$ from the physical spacetime ${\mathcal{M}}_{ϵ}$ to the background spacetime $\mathcal{M}$ is given by ${}_{\mathcal{Y}}{h}_{ab}$ −${}_{\mathcal{X}}{h}_{ab}$=${\pounds }_{{\xi }_{\left(1\right)}}{g}_{ab}$ with the background metric $g_{ab}$, there then exist a tensor field ${\mathcal{F}}_{ab}$ and a vector field $Y^a$ such that $h_{ab}$ is decomposed as $h_{ab}$$=:$${\mathcal{F}}_{ab}$+${\pounds }_Y{g}_{ab}$, where ${\mathcal{F}}_{ab}$ and $Y^a$ are transformed as ${}_{\mathcal{Y}}{\mathcal{F}}_{ab}$−${}_{\mathcal{X}}{\mathcal{F}}_{ab}$= 0 and ${}_{\mathcal{Y}}{Y}^a$−${}_{\mathcal{X}}{Y}^a$=${\xi }_{\left(1\right)}^a$ under the gauge transformation, respectively.

We call ${\mathcal{F}}_{ab}$ and $Y^a$ as the gauge-invariant and gauge-dependent parts of $h_{ab}$, respectively.

The proof of Conjecture 2.1 is highly nontrivial [36,37,38], and it was found that the gauge-invariant variables are essentially non-local. Despite this non-triviality, once we accept Conjecture 2.1, we can decompose the linear perturbation of an arbitrary tensor field ${}_{\mathcal{X}}{}^{\left(1\right)}Q$, whose gauge transformation is given by Equation (4), through the gauge-dependent part $Y_a$ of the metric perturbation in Conjecture 2.1 as

$$\begin{array}{c}\hfill {}_{\mathcal{X}}{}^{\left(1\right)}Q={}^{\left(1\right)}\mathcal{Q}+{\pounds }_{{}_{\mathcal{X}}Y}{Q}_0,\end{array}$$

(7)

where ${}^{\left(1\right)}\mathcal{Q}$ is the gauge-invariant part of the perturbation ${}_{\mathcal{X}}{}^{\left(1\right)}Q$. As examples, the linearized Einstein tensor ${}_{\mathcal{X}}{}^{\left(1\right)}G_a^b$ and the linear perturbation of the energy-momentum tensor ${}_{\mathcal{X}}{}^{\left(1\right)}T_a^b$ are also decomposed as

$$\begin{array}{c}\hfill {}_{\mathcal{X}}{}^{\left(1\right)}{G}_a^b={}^{\left(1\right)}{\mathcal{G}}_a^b\left[\mathcal{F}\right]+{\pounds }_{{}_{\mathcal{X}}Y}{G}_a^b,{}_{\mathcal{X}}{}^{\left(1\right)}{T}_a^b={}^{\left(1\right)}{\mathcal{T}}_a^b+{\pounds }_{{}_{\mathcal{X}}Y}{T}_a^b,\end{array}$$

(8)

where $G_{ab}$ and $T_{ab}$ are the background values of the Einstein tensor and the energy-momentum tensor, respectively [35,46]. Using the background Einstein equation $G_a^b=8\pi T_a^b$, the linearized Einstein equation ${}_{\mathcal{X}}{}^{\left(1\right)}{G}_{ab}=8\pi {}_{\mathcal{X}}{}^{\left(1\right)}{T}_{ab}$ is automatically given in the gauge-invariant form

$$\begin{array}{c}\hfill {}^{\left(1\right)}{\mathcal{G}}_a^b\left[\mathcal{F}\right]=8\pi {}^{\left(1\right)}{\mathcal{T}}_a^b\left[\mathcal{F},\varphi \right]\end{array}$$

(9)

even if the background Einstein equation is nontrivial. Here, “$\varphi$” in Equation (9) symbolically represents the matter degree of freedom.

Finally, we comment on the coordinate transformation induced by the gauge transformation ${\Phi }_{ϵ}$ of the second-kind [12,40]. To see this, we introduce the coordinate system $\{O_{\alpha },{\psi }_{\alpha }\}$ on the background spacetime $\mathcal{M}$, where $O_{\alpha }$ are open sets on the background spacetime and ${\psi }_{\alpha }$ are diffeomorphisms from $O_{\alpha }$ to ${\mathbb{R}}^4$ ($4=dim\mathcal{M}$). The coordinate system $\{O_{\alpha },{\psi }_{\alpha }\}$ is the set of collections of the pair of open sets $O_{\alpha }$ and diffeomorphism ${\psi }_{\alpha }:O_{\alpha }\mapsto {\mathbb{R}}^4$. If we employ a gauge choice ${\mathcal{X}}_{ϵ}$ of the second kind, we have the correspondence of the physical spacetime ${\mathcal{M}}_{ϵ}$ and the background spacetime $\mathcal{M}$. Together with the coordinate system ${\psi }_{\alpha }$ on $\mathcal{M}$, this correspondence ${\mathcal{X}}_{ϵ}$ between ${\mathcal{M}}_{ϵ}$ and $\mathcal{M}$ induces the coordinate system on ${\mathcal{M}}_{ϵ}$. Actually, ${\mathcal{X}}_{ϵ}\left(O_{\alpha }\right)$ for each $\alpha$ is an open set of ${\mathcal{M}}_{ϵ}$. Then, ${\psi }_{\alpha }\circ {\mathcal{X}}_{ϵ}^{-1}$ becomes a diffeomorphism from an open set ${\mathcal{X}}_{ϵ}\left(O_{\alpha }\right)\subset {\mathcal{M}}_{ϵ}$ to ${\mathbb{R}}^4$. This diffeomorphism ${\psi }_{\alpha }\circ {\mathcal{X}}_{ϵ}^{-1}$ induces a coordinate system of an open set on ${\mathcal{M}}_{ϵ}$. When we have two different gauge choices ${\mathcal{X}}_{ϵ}$ and ${\mathcal{Y}}_{ϵ}$ of the second kind, ${\psi }_{\alpha }\circ {\mathcal{X}}_{ϵ}^{-1}$ : ${\mathcal{M}}_{ϵ}$ ↦ ${\mathbb{R}}^4$ ($\left\{x^{\mu }\right\}$) and ${\psi }_{\alpha }\circ {\mathcal{Y}}_{ϵ}^{-1}$ : ${\mathcal{M}}_{ϵ}$ ↦${\mathbb{R}}^4$ ($\left\{y^{\mu }\right\}$) become different coordinate systems on ${\mathcal{M}}_{ϵ}$. We can also consider the coordinate transformation from the coordinate system ${\psi }_{\alpha }\circ {\mathcal{X}}_{ϵ}^{-1}$ to another coordinate system ${\psi }_{\alpha }\circ {\mathcal{Y}}_{ϵ}^{-1}$. Because the gauge transformation ${\mathcal{X}}_{ϵ}\to {\mathcal{Y}}_{ϵ}$ is induced by the diffeomorphism ${\Phi }_{ϵ}$ $:=$ ${\left({\mathcal{X}}_{ϵ}\right)}^{-1}\circ {\mathcal{Y}}_{ϵ}$, this diffeomorphism ${\Phi }_{ϵ}$ induces the coordinate transformation as

$$\begin{array}{c}\hfill y^{\mu }\left(q\right):=x^{\mu }\left(p\right)=\left({\left({\Phi }_{ϵ}^{-1}\right)}^{*}{x}^{\mu }\right)\left(q\right)\end{array}$$

(10)

in the passive point of view [34,55], where p, q∈$\mathcal{M}$ are identified to the same point $``p"$∈${\mathcal{M}}_{ϵ}$ by the gauge choices ${\mathcal{X}}_{ϵ}$ and ${\mathcal{Y}}_{ϵ}$, respectively. If we represent this coordinate transformation in terms of the Taylor expansion (3), we have the coordinate transformation

$$\begin{array}{c}\hfill y^{\mu }\left(q\right)=x^{\mu }\left(q\right)-ϵ{\xi }_{\left(1\right)}^{\mu }\left(q\right)+O\left({ϵ}^2\right).\end{array}$$

(11)

We should emphasize that the coordinate transformation (11) is not the starting point of the gauge transformation but a result of the above framework.

On the other hand, we may consider the point replacement by the one-parameter diffeomorphism $s={\Psi }_{\lambda }\left(r\right)$, where $s,r\in {\mathcal{M}}_{ϵ}$ and $\lambda$ is an infinitesimal parameter satisfying ${\Psi }_{\lambda =0}\left(r\right)=r$. The pull-back ${\Psi }_{\lambda }^{*}$ of any tensor field Q on ${\mathcal{M}}_{ϵ}$ is given by

$$\begin{array}{ccc}\hfill Q\left(s\right)& =& Q\left({\Psi }_{\lambda }\left(r\right)\right)=\left({\Psi }_{\lambda }^{*}\right)Q\left(r\right)\hfill \\ & =& Q\left(r\right)+\lambda {\left.{\pounds }_{\zeta }Q\right|}_{\lambda =0}\left(r\right)+O\left({\lambda }^2\right),\hfill \end{array}$$

(12)

where ${\zeta }^a$ is the generator of the pull-back ${\Psi }_{\lambda }^{*}$. Equation (12) is just the formula of the Taylor expansion on the manifold ${\mathcal{M}}_{ϵ}$ without using any coordinate system. At the same time, Equation (12) is the definitions of the Lie derivative ${\pounds }_{\zeta }$ and its generator ${\zeta }^a$. In many literatures, this formula is derived from the coordinate transformation

$$\begin{array}{c}\hfill y^{\mu }\left(s\right):=x^{\mu }\left(r\right)+\lambda {\zeta }^{\mu }\left(r\right)+O\left({\lambda }^2\right)\end{array}$$

(13)

as explained in the Part I paper [14]. The formula (12) of the Taylor expansion defined without any coordinate system and the second term in the right-hand side of Equation (12) represents the actual difference of the tensor fields $Q\left(s\right)$ and $Q\left(r\right)$ in the different point $s\ne r$ and its physical meaning.

If we consider a formal metric decomposition ${\overline{g}}_{ab}={}^{\left(0\right)}{g}_{ab}+\lambda h_{ab}+O\left({\lambda }^2\right)$ within ${\mathcal{M}}_{ϵ}$ as an example of the tensor field Q in Equation (12), the formula (12) is given by

$$\begin{array}{c}\hfill {}^{\left(0\right)}g_{ab}\left(s\right)+\lambda h_{ab}\left(s\right)={}^{\left(0\right)}g_{ab}\left(r\right)+\lambda \left(h_{ab}\left(r\right)+{\left.{\pounds }_{\zeta }{}^{\left(0\right)}g_{ab}\right|}_r\right)+O\left({\lambda }^2\right).\end{array}$$

(14)

Since the formula of the Taylor expansion (12) is derived from the infinitesimal coordinate transformation rule (13) in many literatures, the term ${\left.{\pounds }_{\zeta }{}^{\left(0\right)}g_{ab}\right|}_r$ in the right-hand side of Equation (14) is often called “the degree of freedom of the coordinate transformation” and “unphysical degree of freedom”. However, since the formula (14) is just the re-expression of the Taylor expansion (12) defined without any coordinate system and the term ${\left.{\pounds }_{\zeta }{}^{\left(0\right)}g_{ab}\right|}_r$ in Equation (14) should have its physical meaning. If we insist that the term ${\left.{\pounds }_{\zeta }{}^{\left(0\right)}g_{ab}\right|}_r$ which appears due to the point replacement on the single manifold is “unphysical”, this directly leads the statement that “the famous arguments of the Killing vector and the symmetry of the spacetime are physically meaningless.” For this reason, we have to emphasize that we cannot regard the second term of the right-hand side of Equation (12) as an “unphysical degree of freedom.”

In our series of papers [12,13,14,15,16] and in this paper, we classify the term ${\left.{\pounds }_{\zeta }{}^{\left(0\right)}g_{ab}\right|}_r$ in Equation (14) as one of gauge-degree of freedom of the first kind, since this term can be eliminate the infinitesimal coordinate transformation (13) which is the coordinate transformation within the single manifold ${\mathcal{M}}_{ϵ}$. Furthermore, Equation (14) does not mean $h_{ab}\left(s\right)$=$h_{ab}\left(r\right)$+ ${\left.{\pounds }_{\zeta }{}^{\left(0\right)}g_{ab}\right|}_r$ but it just means ${}^{\left(0\right)}g_{ab}\left(s\right)$=${}^{\left(0\right)}g_{ab}\left(r\right)$+ $\lambda {\left.{\pounds }_{\zeta }{}^{\left(0\right)}g_{ab}\right|}_r$ +$O\left({\lambda }^2\right)$. Moreover, we also have to emphasize that the infinitesimal coordinate transformation (13) is essentially different from the infinitesimal coordinate transformation (11). Actually, the infinitesimal coordinate transformation (13) is the replacement of the point within the single manifold ${\mathcal{M}}_{ϵ}$. On the other hand, the coordinate transformation (11) is just the change of the coordinate label at the same point in the background spacetime $\mathcal{M}$.

2.2. Linear perturbations on spherically symmetric background spacetime

Here, we consider the 2+2 formulation of the perturbation of a spherically symmetric background spacetime, which originally proposed by Gerlach and Sengupta [25,26,27,28]. Spherically symmetric spacetimes are characterized by the direct product $\mathcal{M}={\mathcal{M}}_1\times S^2$ and their metric is

$$\begin{array}{ccc}\hfill g_{ab}& =& y_{ab}+r^2{\gamma }_{ab},\hfill \end{array}$$

(15)

$$\begin{array}{ccc}\hfill y_{ab}& =& y_{AB}{\left(d{x}^A\right)}_a{\left(d{x}^B\right)}_b,{\gamma }_{ab}={\gamma }_{pq}{\left(d{x}^p\right)}_a{\left(d{x}^q\right)}_b,\hfill \end{array}$$

(16)

where $x^A=(t,r)$, $x^p=(\theta ,\varphi )$, and ${\gamma }_{pq}$ is the metric on the unit sphere. In the Schwarzschild spacetime, the metric (15) is given by

$$\begin{array}{ccc}\hfill y_{ab}& =& -f{\left(dt\right)}_a{\left(dt\right)}_b+f^{-1}{\left(dr\right)}_a{\left(dr\right)}_b,f=1-{\displaystyle \frac{2M}{r}},\hfill \end{array}$$

(17)

$$\begin{array}{ccc}\hfill {\gamma }_{ab}& =& {\left(d\theta \right)}_a{\left(d\theta \right)}_b+{sin}^2\theta {\left(d\varphi \right)}_a{\left(d\varphi \right)}_b={\theta }_a{\theta }_b+{\varphi }_a{\varphi }_b,\hfill \end{array}$$

(18)

$$\begin{array}{ccc}\hfill {\theta }_a& =& {\left(d\theta \right)}_a,{\varphi }_a=sin\theta {\left(d\varphi \right)}_a.\hfill \end{array}$$

(19)

On this background spacetime $(\mathcal{M},g_{ab})$, the components of the metric perturbation is given by

$$\begin{array}{c}\hfill h_{ab}=h_{AB}{\left(d{x}^A\right)}_a{\left(d{x}^B\right)}_b+2h_{Ap}{\left(d{x}^A\right)}_{(a}{\left(d{x}^p\right)}_{b)}+h_{pq}{\left(d{x}^p\right)}_a{\left(d{x}^q\right)}_b.\end{array}$$

(20)

Here, we note that the components $h_{AB}$, $h_{Ap}$, and $h_{pq}$ are regarded as components of scalar, vector, and tensor on $S^2$, respectively. In the Part I paper [14], we showed the linear independence of the set of harmonic functions

$$\begin{array}{c}\hfill \left\{S_{\delta },{\widehat{D}}_p{S}_{\delta },{ϵ}_{pq}{\widehat{D}}^q{S}_{\delta },{\displaystyle \frac{1}{2}}{\gamma }_{pq}{S}_{\delta },\left({\widehat{D}}_p{\widehat{D}}_q-{\displaystyle \frac{1}{2}}{\gamma }_{pq}\right)S_{\delta },2{ϵ}_{r(p}{\widehat{D}}_{q)}{\widehat{D}}^r{S}_{\delta }\right\},\end{array}$$

(21)

where ${\widehat{D}}_p$ is the covariant derivative associated with the metric ${\gamma }_{pq}$ on $S^2$, ${\widehat{D}}^p={\gamma }^{pq}{\widehat{D}}_q$, ${ϵ}_{pq}={ϵ}_{\left[pq\right]}=2{\theta }_{[p}{\varphi }_{q]}$ is the totally antisymmetric tensor on $S^2$. In the set of harmonic function (21), the scalar harmonic function $S_{\delta }$ is given by

$$\begin{array}{c}\hfill S_{\delta }=\left\{\begin{array}{ccccc}{Y}_{lm}& & for& & l\ge 2;\\ k_{(\widehat{\Delta }+2)m}& & for& & l=1;\\ k_{\left(\widehat{\Delta }\right)}& & for& & l=0.\end{array}\right.\end{array}$$

(22)

Here, functions $k_{\left(\widehat{\Delta }\right)}$ and $k_{(\widehat{\Delta }+2)m}$ are the kernel modes of the derivative operator $\widehat{\Delta }$ and $[\widehat{\Delta }+2]$, respectively, and we employ the explicit form of these functions as

$$\begin{array}{ccc}\hfill k_{\left(\widehat{\Delta }\right)}& =& 1+\delta ln{\left({\displaystyle \frac{1-cos\theta }{1+cos\theta }}\right)}^{1/2},\delta \in \mathbb{R},\hfill \end{array}$$

(23)

$$\begin{array}{ccc}\hfill k_{(\widehat{\Delta }+2,m=0)}& =& cos\theta +\delta \left({\displaystyle \frac{1}{2}}cos\theta ln{\displaystyle \frac{1+cos\theta }{1-cos\theta }}-1\right),\delta \in \mathbb{R},\hfill \end{array}$$

(24)

$$\begin{array}{ccc}\hfill k_{(\widehat{\Delta }+2,m=\pm 1)}& =& \left[sin\theta +\delta \left(+{\displaystyle \frac{1}{2}}sin\theta ln{\displaystyle \frac{1+cos\theta }{1-cos\theta }}+cot\theta \right)\right]e^{\pm i\varphi }.\hfill \end{array}$$

(25)

Then, we consider the mode decomposition of the components $\{h_{AB},h_{Ap},h_{pq}\}$ as follows:

$$\begin{array}{ccc}\hfill h_{AB}& =& \sum _{l,m}{\tilde{h}}_{AB}{S}_{\delta },\hfill \end{array}$$

(26)

$$\begin{array}{ccc}\hfill h_{Ap}& =& r\sum _{l,m}\left[{\tilde{h}}_{\left(e1\right)A}{\widehat{D}}_p{S}_{\delta }+{\tilde{h}}_{\left(o1\right)A}{ϵ}_{pq}{\widehat{D}}^q{S}_{\delta }\right],\hfill \end{array}$$

(27)

$$\begin{array}{ccc}\hfill h_{pq}& =& r^2\sum _{l,m}\left[{\displaystyle \frac{1}{2}}{\gamma }_{pq}{\tilde{h}}_{\left(e0\right)}{S}_{\delta }+{\tilde{h}}_{\left(e2\right)}\left({\widehat{D}}_p{\widehat{D}}_q-{\displaystyle \frac{1}{2}}{\gamma }_{pq}{\widehat{D}}^r{\widehat{D}}_r\right)S_{\delta }+2{\tilde{h}}_{\left(o2\right)}{ϵ}_{r(p}{\widehat{D}}_{q)}{\widehat{D}}^r{S}_{\delta }\right].\hfill \end{array}$$

(28)

Since the linear independence of each element of the set of harmonic function (21) is guaranteed, the one-to-one correspondence between the components $\{h_{AB},$$h_{Ap},$ $h_{pq}\}$ and the mode coefficients $\{{\tilde{h}}_{AB},$ ${\tilde{h}}_{\left(e1\right)A},$${\tilde{h}}_{\left(o1\right)A},$${\tilde{h}}_{\left(e0\right)},$ ${\tilde{h}}_{\left(e2\right)},$${\tilde{h}}_{\left(o2\right)}\}$ with the decomposition formulae (26)-(28) is guaranteed including $l=0,1$ mode if $\delta \ne 0$. Then, the mode-by-mode analysis including $l=0,1$ is possible when $\delta \ne 0$. However, the mode functions (23)–(25) are singular if $\delta \ne 0$. When $\delta =0$, we have $k_{\left(\widehat{\Delta }\right)}\propto Y_{00}$ and $k_{(\widehat{\Delta }+2)m}\propto Y_{1m}$, but the linear dependence of the set of harmonics (21) is lost in this case. Because of this situation, we proposed the following strategy:

Proposal 2.1.

We decompose the metric perturbation $h_{ab}$ on the background spacetime with the metric (15)–(18) through Equations (26)–(28) with the harmonic function $S_{\delta }$ given by Equation (22). Then, Equations (26)–(28) become invertible including $l=0,1$ modes. After deriving the mode-by-mode field equations such as linearized Einstein equations by using the harmonic functions $S_{\delta }$, we choose $\delta =0$ as the regular boundary condition for solutions when we solve these field equations.

As shown in the Part I paper [14], once we accept Proposal 2.1, the Conjecture 2.1 becomes the following statement:

Theorem 1.

If the gauge-transformation rule for a perturbative pulled-back tensor field $h_{ab}$ to the background spacetime $\mathcal{M}$ is given by ${}_{\mathcal{Y}}{h}_{ab}$−${}_{\mathcal{X}}{h}_{ab}$= ${\pounds }_{{\xi }_{\left(1\right)}}{g}_{ab}$ with the background metric $g_{ab}$ with spherically symmetry, there then exist a tensor field ${\mathcal{F}}_{ab}$ and a vector field $Y^a$ such that $h_{ab}$ is decomposed as $h_{ab}$$=:$${\mathcal{F}}_{ab}$+ ${\pounds }_Y{g}_{ab}$, where ${\mathcal{F}}_{ab}$ and $Y^a$ are transformed into ${}_{\mathcal{Y}}{\mathcal{F}}_{ab}$− ${}_{\mathcal{X}}{\mathcal{F}}_{ab}$= 0 and ${}_{\mathcal{Y}}{Y}^a$ −${}_{\mathcal{X}}{Y}^a$=${\xi }_{\left(1\right)}^a$ under the gauge transformation, respectively.

Actually, the gauge-dependent variable $Y_a$ is given by

$$\begin{array}{ccc}\hfill Y_a& :=& \sum _{l,m}{\tilde{Y}}_A{S}_{\delta }{\left(d{x}^A\right)}_a+\sum _{l,m}\left({\tilde{Y}}_{\left(e1\right)}{\widehat{D}}_p{S}_{\delta }+{\tilde{Y}}_{\left(o1\right)}{ϵ}_{pq}{\widehat{D}}^q{S}_{\delta }\right){\left(d{x}^p\right)}_a,\hfill \end{array}$$

(29)

where

$$\begin{array}{ccc}\hfill {\tilde{Y}}_A& :=& r{\tilde{h}}_{\left(e1\right)A}-{\displaystyle \frac{r^2}{2}}{\overline{D}}_A{\tilde{h}}_{\left(e2\right)},\hfill \end{array}$$

(30)

$$\begin{array}{ccc}\hfill {\tilde{Y}}_{\left(e1\right)}& :=& {\displaystyle \frac{r^2}{2}}{\tilde{h}}_{\left(e2\right)},\hfill \end{array}$$

(31)

$$\begin{array}{ccc}\hfill {\tilde{Y}}_{\left(o1\right)}& :=& -r^2{\tilde{h}}_{\left(o2\right)}.\hfill \end{array}$$

(32)

Furthermore, including $l=0,1$ modes, the components of the gauge-invariant part ${\mathcal{F}}_{ab}$ of the metric perturbation $h_{ab}$ is given by

$$\begin{array}{ccc}\hfill {\mathcal{F}}_{AB}& =& \sum _{l,m}{\tilde{F}}_{AB}{S}_{\delta },\hfill \end{array}$$

(33)

$$\begin{array}{ccc}\hfill {\mathcal{F}}_{Ap}& =& r\sum _{l,m}{\tilde{F}}_A{ϵ}_{pq}{\widehat{D}}^q{S}_{\delta },{\widehat{D}}^p{\mathcal{F}}_{Ap}=0,\hfill \end{array}$$

(34)

$$\begin{array}{ccc}\hfill {\mathcal{F}}_{pq}& =& {\displaystyle \frac{1}{2}}{\gamma }_{pq}{r}^2\sum _{l,m}\tilde{F}{S}_{\delta },\hfill \end{array}$$

(35)

where ${\tilde{F}}_{AB}$, ${\tilde{F}}_A$, and $\tilde{F}$ are given by

$$\begin{array}{ccc}\hfill {\tilde{F}}_{AB}& :=& {\tilde{h}}_{AB}-2{\overline{D}}_{(A}{\tilde{Y}}_{B)},\hfill \end{array}$$

(36)

$$\begin{array}{ccc}\hfill {\tilde{F}}_A& :=& {\tilde{h}}_{\left(o1\right)A}+r{\overline{D}}_A{\tilde{h}}_{\left(o2\right)},\hfill \end{array}$$

(37)

$$\begin{array}{ccc}\hfill \tilde{F}& :=& {\tilde{h}}_{\left(e0\right)}-{\displaystyle \frac{4}{r}}{\tilde{Y}}_A{\overline{D}}^{A}r+{\tilde{h}}_{\left(e2\right)}l(l+1).\hfill \end{array}$$

(38)

Thus, we have constructed gauge-invariant metric perturbations on the Schwarzschild background spacetime including $l=0,1$ modes.

Furthermore, from Equations (33)–() for the gauge-invariant variables ${\mathcal{F}}_{AB}$, ${\mathcal{F}}_{Ap}$, and ${\mathcal{F}}_{pq}$ and gauge-dependent variables $Y_a$ defined by Equations (29)–(), the original components $\{h_{AB},h_{Ap},h_{pq}\}$ of the metric perturbation (26)–() are given by

$$\begin{array}{ccc}\hfill h_{ab}& =& {\mathcal{F}}_{ab}+{\pounds }_Y{g}_{ab}.\hfill \end{array}$$

(39)

This is the assertion of Theorem 1.

The aim of our proposal is to add a greater degree of freedom of the metric perturbation so that the decomposition (39) is guaranteed even in the case of $l=0,1$ modes. Therefore, it is impossible to reach the final expression (39) if we treat the metric perturbation by using $S_{\delta =0}=Y_{lm}$ in the decomposition formulae (26)–(28) from the starting point.

To see this more specifically, we show an explicit example of the $l=1$ odd mode perturbation. If we apply our proposal to the decomposition formulae (26)–(28), $l=1$ odd-mode perturbation is given by

$$\begin{array}{c}\hfill h_{ab}=2r{\tilde{h}}_{\left(o1\right)A}^{(l=1)}{ϵ}_{pq}{\widehat{D}}^q{S}_{\delta }{\left(d{x}^A\right)}_{(a}{\left(d{x}^p\right)}_{b)}+2r^2{\tilde{h}}_{\left(o2\right)}^{(l=1)}{ϵ}_{r(p}{\widehat{D}}_{q)}{\widehat{D}}^r{S}_{\delta }{\left(d{x}^p\right)}_{(a}{\left(d{x}^q\right)}_{b)}.\end{array}$$

(40)

If we use $S_{\delta =0}=Y_{lm}$ from the starting point, $2{ϵ}_{r(p}{\widehat{D}}_{q)}{\widehat{D}}^r{S}_{\delta }=0$. So, the mode coefficient ${\tilde{h}}_{\left(o2\right)}^{(l=1)}$ does not appear. When $\delta \ne 0$, the gauge-transformation rule of the variables ${\tilde{h}}_{\left(o1\right)A}^{(l=1)}$ and ${\tilde{h}}_{\left(o2\right)}^{(l=1)}$ are given by

$$\begin{array}{c}\hfill {}_{\mathcal{Y}}{\tilde{h}}_{\left(o1\right)A}^{(l=1)}-{}_{\mathcal{X}}{\tilde{h}}_{\left(o1\right)A}^{(l=1)}=r{\overline{D}}_A\left({\displaystyle \frac{1}{r}}{\zeta }_{\left(o\right)}\right),{}_{\mathcal{Y}}{\tilde{h}}_{\left(o2\right)}^{(l=1)}-{}_{\mathcal{X}}{\tilde{h}}_{\left(o2\right)}^{(l=1)}=-{\displaystyle \frac{1}{r}}{\zeta }_{\left(o\right)},\end{array}$$

(41)

where the generator of the gauge transformation is given by ${\xi }_A=0$ and ${\xi }_p=r{\zeta }_{\left(o\right)}{ϵ}_{pq}{\widehat{D}}^q{S}_{\delta }$. If we use $S_{\delta =0}=Y_{lm}$ from the starting point, the gauge-transformation rule for the variable ${\tilde{h}}_{\left(o2\right)}^{(l=1)}$ does not appear, either. However, in our case, this gauge transformation appears due to the fact ${ϵ}_{r(p}{\widehat{D}}_{q)}{\widehat{D}}^r{S}_{\delta }\ne 0$. Inspecting these gauge-transformation rules, we rewrite Equation (40) as

$$\begin{array}{ccc}\hfill h_{ab}& =& 2r\left[\left({\tilde{h}}_{\left(o1\right)A}^{(l=1)}+r{\overline{D}}_A{\tilde{h}}_{\left(o2\right)}\right)-r{\overline{D}}_A{\tilde{h}}_{\left(o2\right)}\right]{ϵ}_{pq}{\widehat{D}}^q{S}_{\delta }{\left(d{x}^A\right)}_{(a}{\left(d{x}^p\right)}_{b)}\hfill \\ & & +2r^2{\tilde{h}}_{\left(o2\right)}^{(l=1)}{ϵ}_{r(p}{\widehat{D}}_{q)}{\widehat{D}}^r{S}_{\delta }{\left(d{x}^p\right)}_{(a}{\left(d{x}^q\right)}_{b)}\hfill \\ & =& 2r\left[{\tilde{F}}_A-r{\overline{D}}_A\left({\displaystyle \frac{1}{r^2}}{\tilde{Y}}_{\left(o1\right)}\right)\right]{ϵ}_{pq}{\widehat{D}}^q{S}_{\delta }{\left(d{x}^A\right)}_{(a}{\left(d{x}^p\right)}_{b)}\hfill \\ & & +2r^2\left({\displaystyle \frac{1}{r^2}}{\tilde{Y}}_{\left(o1\right)}\right){ϵ}_{r(p}{\widehat{D}}_{q)}{\widehat{D}}^r{S}_{\delta }{\left(d{x}^p\right)}_{(a}{\left(d{x}^q\right)}_{b)}.\hfill \end{array}$$

(42)

In Equation (42), ${\tilde{F}}_A$ is the gauge-invariant variable for the metric perturbation defined by Equation (37) and ${\tilde{Y}}_{\left(o1\right)}$ is the gauge-dependent part of the metric perturbation defined by Equation (32). The terms of ${\tilde{Y}}_{\left(o1\right)}$ in Equation (42) can be written in the form of the Lie derivative of the background metric as in Equation (39). When $\delta =0$, Equation (42) is given by

$$\begin{array}{ccc}\hfill h_{ab}& =& 2r\left[{\tilde{F}}_A-r{\overline{D}}_A\left({\displaystyle \frac{1}{r^2}}{\tilde{Y}}_{\left(o1\right)}\right)\right]{ϵ}_{pq}{\widehat{D}}^q{S}_{\delta }{\left(d{x}^A\right)}_{(a}{\left(d{x}^p\right)}_{b)}.\hfill \end{array}$$

(43)

This is also given in the form (39).

When we construct the gauge-invariant variable ${\tilde{F}}_A$ in Equations (42) and (43), the existence of the coefficient ${\tilde{h}}_{\left(o2\right)}^{(l=1)}={\tilde{Y}}_{\left(o1\right)}/r^2$ is essential which does not appear in the treatment using the spherical harmonics $Y_{lm}$ from the starting point. Actually, we cannot define the gauge-invariant variable for $l=1$ odd-mode perturbation in the same manner as $l\ge 2$-mode case. Our proposal 2.1, make the degree of freedom of the metric perturbations increase. We have to emphasize that we can reach the decomposition (39) owing to this additional degree of freedom.

To discuss the linearized Einstein equation (9) and the linear perturbation of the continuity equation

$$\begin{array}{c}\hfill {\nabla }_a{}^{\left(1\right)}{\mathcal{T}}_b^a=0\end{array}$$

(44)

of the gauge-invariant energy-momentum tensor ${}^{\left(1\right)}{\mathcal{T}}_b^a:=g^{ac}{}^{\left(1\right)}{\mathcal{T}}_{bc}$ on a vacuum background spacetime, we consider the mode-decomposition of the gauge-invariant part ${}^{\left(1\right)}{\mathcal{T}}_{bc}$ of the linear perturbation of the energy-momentum tensor through the set (21) of the harmonics as follows:

$$\begin{array}{ccc}\hfill {}^{\left(1\right)}{\mathcal{T}}_{ab}& =& \sum _{l,m}{\tilde{T}}_{AB}{S}_{\delta }{\left(d{x}^A\right)}_a{\left(d{x}^B\right)}_b+r\sum _{l,m}\left\{{\tilde{T}}_{\left(e1\right)A}{\widehat{D}}_p{S}_{\delta }+{\tilde{T}}_{\left(o1\right)A}{ϵ}_{pr}{\widehat{D}}^r{S}_{\delta }\right\}2{\left(d{x}^A\right)}_{(a}{\left(d{x}^p\right)}_{b)}\hfill \\ & & +r^2\sum _{l,m}\left\{{\tilde{T}}_{\left(e0\right)}{\displaystyle \frac{1}{2}}{\gamma }_{pq}{S}_{\delta }+{\tilde{T}}_{\left(e2\right)}\left({\widehat{D}}_p{\widehat{D}}_q{S}_{\delta }-{\displaystyle \frac{1}{2}}{\gamma }_{pq}{\widehat{D}}_r{\widehat{D}}^r{S}_{\delta }\right)\right.\hfill \\ & & \left.+{\tilde{T}}_{\left(o2\right)}2{ϵ}_{r(p}{\widehat{D}}_{q)}{\widehat{D}}^r{S}_{\delta }\right\}{\left(d{x}^p\right)}_a{\left(d{x}^q\right)}_b.\hfill \end{array}$$

(45)

In the Part I paper [14], we derived the linearized Einstein equations, discussed the odd-mode perturbation ${\tilde{F}}_{Ap}$ in Equation (34), and derived the $l=1$ odd-mode solutions to these equations. The Einstein equation for even mode ${\tilde{F}}_{AB}$ and $\tilde{F}$ in Equations (33) and (35) also derived in the Part I paper [14], and $l=0,1$ even-mode solutions are derived in the Part II paper [15]. Since these solutions include the Kerr parameter perturbation and the Schwarzschild mass parameter perturbation of the linear order in the vacuum case, these are physically reasonable. Then, we conclude that our proposal is also physically reasonable. Furthermore, we also checked that our derived solutions include the linearized LTB solution and non-rotating C-metric with the Schwarzschild background in the Part III paper [16].

The purpose of this paper is to compare our proposed gauge-invariant treatments for $l=0,1$-mode perturbations on the Schwarzschild background spacetime with conventional “complete gauge-fixing treatments” in which we use the spherical harmonics $Y_{lm}$ from the starting point. For this purpose, the linearized Einstein equations for the $l=0,1$-modes and their solutions based on our proposal are necessary. Therefore, we review them below.

2.3. $l=1$ Odd-Mode Linearized Einstein Equations and Solutions

As derived in the Part I paper [14], the $l=1$ odd-mode part in the linearized Einstein equations are simplified as the constraint equation

$$\begin{array}{c}\hfill {\overline{D}}_D\left(r{\tilde{F}}^D\right)=0,\end{array}$$

(46)

and the evolution equation

$$\begin{array}{ccc}& & -\left[{\overline{D}}^D{\overline{D}}_D-{\displaystyle \frac{2}{r^2}}\right]\left(r{\tilde{F}}_A\right)-{\displaystyle \frac{2}{r^2}}\left({\overline{D}}^{D}r\right)\left({\overline{D}}_{A}r\right)\left(r{\tilde{F}}_D\right)+{\displaystyle \frac{2}{r}}\left({\overline{D}}^{D}r\right){\overline{D}}_A\left(r{\tilde{F}}_D\right)\hfill \\ & =& 16\pi r{\tilde{T}}_{\left(o1\right)A},\hfill \end{array}$$

(47)

where we choose ${\tilde{T}}_{\left(o2\right)}=0$ by hand. Furthermore, we have the continuity equation

$$\begin{array}{c}\hfill {\overline{D}}^C{\tilde{T}}_{\left(o1\right)C}+{\displaystyle \frac{3}{r}}\left({\overline{D}}^{D}r\right){\tilde{T}}_{\left(o1\right)D}=0\end{array}$$

(48)

for the $l=1$ odd-mode matter perturbation which is derived from the divergence of the first-order perturbation of the energy-momentum tensor.

We also note that the constraint (46) comes from the traceless part of the $(q,p)$-component of the Einstein equation (9), which is the coefficient of the mode function $2{ϵ}_{r(p}{\widehat{D}}_{q)}{\widehat{D}}^p{S}_{\delta }$. If $\delta =0$, i.e., the scalar harmonics $S_{\delta }$ is the spherical harmonics $Y_{lm}$, $2{ϵ}_{r(p}{\widehat{D}}_{q)}{\widehat{D}}^p{S}_{\delta }=0$ for $l=1$. Therefore, the constraint (46) does not appear when we use the spherical harmonics $Y_{lm}$ from the starting point.

The explicit strategy to solve these odd-mode perturbations and $l=0,1$ mode solutions was discussed in the Part I paper [14]. As the result, we obtained the $l=1$ odd-mode “formal” solution as follows:

$$\begin{array}{c}2{\mathcal{F}}_{Ap}{\left(d{x}^A\right)}_{(a}{\left(d{x}^p\right)}_{b)}=6M{r}^2\left[\int dr{\displaystyle \frac{a_1(t,r)}{r^4}}\right]{sin}^2\theta {\left(dt\right)}_{(a}{\left(d\varphi \right)}_{b)}+{\pounds }_V{g}_{ab},\hfill \end{array}$$

(49)

$$\begin{array}{c}{V}_a=\left({\beta }_{1}t+{\beta }_0+W_{\left(o\right)}(t,r)\right)r^2{sin}^2\theta {\left(d\varphi \right)}_a,\hfill \end{array}$$

(50)

where

$$\begin{array}{ccc}\hfill a_1(t,r)& =& -{\displaystyle \frac{16\pi }{3M}}{r}^{3}f\int dt{\tilde{T}}_{\left(o1\right)r}+a_{10}\hfill \\ & =& -{\displaystyle \frac{16\pi }{3M}}\int dr{r}^3{\displaystyle \frac{1}{f}}{\tilde{T}}_{\left(o1\right)t}+a_{10}.\hfill \end{array}$$

(51)

The constant $a_{10}$ in Equation (51) corresponds to the Kerr parameter. Furthermore, the function $W_{\left(o\right)}(t,r)$ in Equation (50) is related to the solution to the $l=1$-mode Regge-Wheeler equation

$$\begin{array}{c}\hfill {\partial }_t^2{Z}_{\left(o\right)}-f{\partial }_r\left(f{\partial }_r{Z}_{\left(o\right)}\right)+{\displaystyle \frac{1}{r^2}}f\left[l(l+1)-3(1-f)\right]Z_{\left(o\right)}=16\pi f^2{\tilde{T}}_{\left(o1\right)r},\end{array}$$

(52)

where $Z_{\left(o\right)}=rf{\partial }_r{W}_{\left(o\right)}$. We have to solve Equation (52) to obtain the function $W_{\left(o\right)}(t,r)$. In this sense, the solution () should be regarded as the “formal” one.

2.4. $l=0,1$ even-mode linearized Einstein equations and solutions

The $l=0,1$ even-mode part of the linearized Einstein equation (9) is summarized as follows:

$$\begin{array}{ccc}\hfill {\tilde{F}}_D^D& =& 0,\hfill \end{array}$$

(53)

$$\begin{array}{ccc}\hfill {\overline{D}}^D{\tilde{\mathbb{F}}}_{AD}-{\displaystyle \frac{1}{2}}{\overline{D}}_A\tilde{F}& =& 16\pi r{\tilde{T}}_{\left(e1\right)A},\hfill \end{array}$$

(54)

where the variable ${\tilde{\mathbb{F}}}_{AB}$ is the traceless part of the variable ${\tilde{F}}_{AB}$ defined by

$$\begin{array}{c}\hfill {\tilde{\mathbb{F}}}_{AB}:={\tilde{F}}_{AB}-{\displaystyle \frac{1}{2}}{y}_{AB}{\tilde{F}}_C^C\end{array}$$

(55)

and we choose the component of the energy-momentum tensor so that ${\tilde{T}}_{\left(e2\right)}=0$ by hand. Owing to the same choice ${\tilde{T}}_{\left(e2\right)}=0$, we also have the evolution equations

$$\begin{array}{ccc}& & \left({\overline{D}}_D{\overline{D}}^D+{\displaystyle \frac{2}{r}}\left({\overline{D}}^{D}r\right){\overline{D}}_D-{\displaystyle \frac{(l-1)(l+2)}{r^2}}\right)\tilde{F}-{\displaystyle \frac{4}{r^2}}\left({\overline{D}}_{C}r\right)\left({\overline{D}}_{D}r\right){\tilde{\mathbb{F}}}^{CD}=16\pi S_{\left(F\right)},\hfill \end{array}$$

(56)

$$\begin{array}{ccc}& & S_{\left(F\right)}:={\tilde{T}}_C^C+4\left({\overline{D}}_{D}r\right){\tilde{T}}_{\left(e1\right)}^D,\hfill \\ & & \left[-{\overline{D}}_D{\overline{D}}^D-{\displaystyle \frac{2}{r}}\left({\overline{D}}_{D}r\right){\overline{D}}^D+{\displaystyle \frac{4}{r}}\left({\overline{D}}^D{\overline{D}}_{D}r\right)+{\displaystyle \frac{l(l+1)}{r^2}}\right]{\tilde{\mathbb{F}}}_{AB}\hfill \\ & & +{\displaystyle \frac{4}{r}}\left({\overline{D}}^{D}r\right){\overline{D}}_{(A}{\tilde{\mathbb{F}}}_{B)D}-{\displaystyle \frac{2}{r}}\left({\overline{D}}_{(A}r\right){\overline{D}}_{B)}\tilde{F}\hfill \end{array}$$

(57)

$$\begin{array}{ccc}& & =16\pi S_{\left(\mathbb{F}\right)AB},\hfill \end{array}$$

(58)

$$\begin{array}{ccc}\hfill S_{\left(\mathbb{F}\right)AB}& :=& {\tilde{T}}_{AB}-{\displaystyle \frac{1}{2}}{y}_{AB}{\tilde{T}}_C^C-2\left({\overline{D}}_{(A}\left(r{\tilde{T}}_{\left(e1\right)B)}\right)-{\displaystyle \frac{1}{2}}{y}_{AB}{\overline{D}}^D\left(r{\tilde{T}}_{\left(e1\right)D}\right)\right)\hfill \\ & & +2y_{AB}\left({\overline{D}}^{C}r\right){\tilde{T}}_{\left(e1\right)C},\hfill \end{array}$$

(59)

for the variable $\tilde{F}$ and the traceless variable ${\tilde{\mathbb{F}}}_{AB}$ with $l=0,1$. We also have to take into account the even-mode part of the continuity equation as follows:

$$\begin{array}{ccc}& & {\overline{D}}^C{\tilde{T}}_C^B+{\displaystyle \frac{2}{r}}\left({\overline{D}}^{D}r\right){\tilde{T}}_D^B-{\displaystyle \frac{1}{r}}l(l+1){\tilde{T}}_{\left(e1\right)}^B-{\displaystyle \frac{1}{r}}\left({\overline{D}}^{B}r\right){\tilde{T}}_{\left(e0\right)}=0,\hfill \end{array}$$

(60)

$$\begin{array}{ccc}& & {\overline{D}}^C{\tilde{T}}_{\left(e1\right)C}+{\displaystyle \frac{3}{r}}\left({\overline{D}}^{C}r\right){\tilde{T}}_{\left(e1\right)C}+{\displaystyle \frac{1}{2r}}{\tilde{T}}_{\left(e0\right)}=0,\hfill \end{array}$$

(61)

for $l=0,1$.

We also note that the constraint (53) comes from the traceless part of the $(q,p)$-component of the Einstein tensor (9), which is the coefficient of the mode function $\left({\widehat{D}}_p{\widehat{D}}_q-{\displaystyle \frac{1}{2}}{\gamma }_{pq}{\widehat{D}}^r{\widehat{D}}_r\right)S_{\delta }$. If $\delta =0$, i.e., the scalar harmonics $S_{\delta }$ is the spherical harmonics $Y_{lm}$, $\left({\widehat{D}}_p{\widehat{D}}_q-{\displaystyle \frac{1}{2}}{\gamma }_{pq}{\widehat{D}}^r{\widehat{D}}_r\right)S_{\delta }=0$ for $l=1$. Therefore, the constraint (53) does not appear when we use the spherical harmonics $Y_{lm}$ from the starting point.

2.4.1. $l=1$ Even-Mode Solution

In the Part II paper [15], we derived the $l=1$ solution to Equations (53), (), (56), (), (60), and () with $l=1$. For $m=0$ mode, in the Part II paper [15], we derived the following “formal” solution to the linearized Einstein equation

$$\begin{array}{ccc}\hfill {\mathcal{F}}_{ab}& =& {\pounds }_V{g}_{ab}-{\displaystyle \frac{16\pi r^2}{3(1-f)}}\left[f^2\left\{{\displaystyle \frac{1+f}{2}}{\tilde{T}}_{rr}+rf{\partial }_r{\tilde{T}}_{rr}-{\tilde{T}}_{\left(e0\right)}-4{\tilde{T}}_{\left(e1\right)r}\right\}{\left(dt\right)}_a{\left(dt\right)}_b\right.\hfill \\ & & +{\displaystyle \frac{2r}{f}}\left\{{\partial }_t{\tilde{T}}_{tt}-{\displaystyle \frac{3f(1-f)}{2r}}{\tilde{T}}_{tr}\right\}{\left(dt\right)}_{(a}{\left(dr\right)}_{b)}\hfill \\ & & +{\displaystyle \frac{r}{f}}\left\{{\partial }_r{\tilde{T}}_{tt}-{\displaystyle \frac{3(1-3f)}{2rf}}{\tilde{T}}_{tt}\right\}{\left(dr\right)}_a{\left(dr\right)}_b\hfill \\ & & \left.+r^2{\tilde{T}}_{tt}{\gamma }_{ab}\right]cos\theta ,\hfill \end{array}$$

(62)

where the vector field $V_a$ is given by

$$\begin{array}{ccc}\hfill V_a& :=& -r{\partial }_t{\Phi }_{\left(e\right)}cos\theta {\left(dt\right)}_a+\left({\Phi }_{\left(e\right)}-r{\partial }_r{\Phi }_{\left(e\right)}\right)cos\theta {\left(dr\right)}_a-r{\Phi }_{\left(e\right)}sin\theta {\left(d\theta \right)}_a.\hfill \end{array}$$

(63)

Here, the variable ${\Phi }_{\left(e\right)}$ is a solution to the $l=1$ Zerilli equation

$$\begin{array}{ccc}& & -{\partial }_t^2{\Phi }_{\left(e\right)}+f{\partial }_r\left[f{\partial }_r{\Phi }_{\left(e\right)}\right]-{\displaystyle \frac{f(1-f)}{r^2}}{\Phi }_{\left(e\right)}\hfill \\ & =& {\displaystyle \frac{4\pi r}{3(1-f)}}\left[3(1-3f){\tilde{T}}_{tt}+(1+f)f^2{\tilde{T}}_{rr}-2rf{\partial }_r{\tilde{T}}_{tt}\right.\hfill \\ & & \left.+2r{f}^3{\partial }_r{\tilde{T}}_{rr}-2f^2{\tilde{T}}_{\left(e0\right)}-8f^2{\tilde{T}}_{\left(e1\right)r}\right].\hfill \end{array}$$

(64)

If we obtain the solution to the $l=1$ Zerilli equation (64), we can write the explicit form of the solution (62) through the generator $V_a$ defined by Equation (63). In this sense, we have to regard the solution (62) is a “formal” one.

2.4.2. $l=0$ Even-Mode Solution

On the other hand, for the $l=0$-mode, we may choose ${\tilde{T}}_{\left(e1\right)A}=0$ in Equations (53), (54), (56)–(61) with $l=0$. Then, we derived the $l=0$ mode solution

$$\begin{array}{ccc}\hfill {\mathcal{F}}_{ab}& =& {\displaystyle \frac{2}{r}}\left(M_1+4\pi \int dr{\displaystyle \frac{r^2}{f}}{T}_{tt}\right)\left({\left(dt\right)}_a{\left(dt\right)}_a+{\displaystyle \frac{1}{f^2}}{\left(dr\right)}_a{\left(dr\right)}_a\right)\hfill \\ & & +2\left[4\pi r\int dt\left({\displaystyle \frac{1}{f}}{\tilde{T}}_{tt}+f{\tilde{T}}_{rr}\right)\right]{\left(dt\right)}_{(a}{\left(dr\right)}_{b)}+{\pounds }_V{g}_{ab},\hfill \end{array}$$

(65)

where

$$\begin{array}{c}\hfill V_a=\left({\displaystyle \frac{f}{4}}{\rm Y}+{\displaystyle \frac{rf}{4}}{\partial }_r{\rm Y}-{\displaystyle \frac{r\Xi \left(r\right)}{(1-3f)}}+f\int dr{\displaystyle \frac{2\Xi \left(r\right)}{f{(1-3f)}^2}}\right){\left(dt\right)}_a+{\displaystyle \frac{1}{4f}}r{\partial }_t{\rm Y}{\left(dr\right)}_a.\end{array}$$

(66)

Here, the variable $\tilde{F}=:{\partial }_t{\rm Y}$ must satisfy the equation

$$\begin{array}{ccc}& & -{\displaystyle \frac{1}{f}}{\partial }_t^2{\rm Y}+{\partial }_r\left(f{\partial }_r{\rm Y}\right)+{\displaystyle \frac{3(1-f)}{r^2}}{\rm Y}-{\displaystyle \frac{8}{r^3}}\int dt{m}_1(t,r)-{\displaystyle \frac{4}{1-3f}}{\partial }_r\Xi \left(r\right)\hfill \\ & =& 16\pi \int dt\left(-{\displaystyle \frac{1}{f}}{\tilde{T}}_{tt}+f{\tilde{T}}_{rr}\right),\hfill \end{array}$$

(67)

where $m_1(t,r)$ is given by

$$\begin{array}{c}\hfill m_1(t,r)=4\pi \int dr{\displaystyle \frac{r^2}{f}}{\tilde{T}}_{tt}+M_1=4\pi \int dt{r}^{2}f{\tilde{T}}_{tr}+M_1,\end{array}$$

(68)

$M_1$ is the constant corresponding to the Schwarzschild mass perturbation, and $\Xi \left(r\right)$ is an arbitrary function of r.

3. Rule of Comparison and Gauge-Transformation Rules in a Conventional Gauge-Fixing

The main purpose of this paper is to check the following statement. “Gauge invariant formulations of perturbations are equivalent to complete gauge fixing approaches.” However, this statement is too ambiguous to check if we have the background knowledge explained in Sec. 2. For example, there is no explanation of the terminology “gauge” in this statement. Therefore, we have to clarify this statement as follows:

Our rule for comparison

First of all, we assume the terminology “gauge” in the statement “Gauge invariant formulations of perturbations are equivalent to complete gauge fixing approaches” is the second-kind gauge which is explained in Sec. 2.1. Therefore, the gauge-transformation rule for this statement is given by (4). Furthermore, the degree of freedom that changes under this gauge-transformation rule is regarded as “unphysical.” The “gauge fixing” in the above statement is a specification of some perturbative variables through the degree of freedom of the generator ${\xi }_{\left(1\right)}^a$. Furthermore, “complete gauge fixing” in the above statement is a specification of some perturbative variables through the “entire” degree of freedom of the generator ${\xi }_{\left(1\right)}^a$.

These are all rules of our game to compare our gauge-invariant formulation and a “conventional complete gauge-fixing approach” in which we use the spherical harmonics $Y_{lm}$ from the starting point.

3.1. Metric Perturbations

Based on the above conceptual premise, we consider the metric perturbation $h_{ab}$ on the background Schwarzschild spacetime $(\mathcal{M},g_{ab})$ whose metric $g_{ab}$ is given by Equations (15)–(19). We consider the components of the metric perturbation $h_{ab}$ as Equation (20) and the decomposition of these components $\{h_{AB},h_{Ap},h_{pq}\}$ as Equations (26)–(28) but we concentrate on the case $S_{\delta }=Y_{lm}=:S$, where $Y_{lm}$ is the conventional spherical harmonics.

Since we choose $S=Y_{lm}$, we have ${\widehat{D}}_{p}S$= ${ϵ}_{pq}{\widehat{D}}^{q}S$= 0, $\left({\widehat{D}}_p{\widehat{D}}_q-{\displaystyle \frac{1}{2}}{\gamma }_{pq}{\widehat{D}}^r{\widehat{D}}_r\right)S$ =$2{ϵ}_{r(p}{\widehat{D}}_{q)}{\widehat{D}}^{r}S$= 0 for $l=0$ modes. For $l=1$ modes, $\left({\widehat{D}}_p{\widehat{D}}_q-{\displaystyle \frac{1}{2}}{\gamma }_{pq}{\widehat{D}}^r{\widehat{D}}_r\right)S$ =$2{ϵ}_{r(p}{\widehat{D}}_{q)}{\widehat{D}}^{r}S$= 0. Due to these facts, we cannot construct gauge-invariant variables for $l=0,1$ modes in a similar manner to the derivation of Equations (36)–() for $l\ge 2$ modes. Then, we cannot use the gauge-invariant formulation reviewed in Sec. 2. Instead, we have to fix the second-kind gauge degree of freedom to exclude the “unphysical” degree of freedom.

3.2. Conventional gauge-transformation rules for $l=0,1$ metric perturbations

Here, we consider the gauge-transformation rule

$$\begin{array}{c}\hfill {}_{\mathcal{Y}}{h}_{ab}-{}_{\mathcal{X}}{h}_{ab}={\pounds }_{\xi }{g}_{ab}=2{\nabla }_{(a}{\xi }_{b)}.\end{array}$$

(69)

From this gauge-transformation rule, we can derive the gauge-transformation rule for the components $\{h_{AB},h_{Ap},h_{qp}\}$ defined by Equation (20) is given by

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}{h}_{AB}-{}_{\mathcal{X}}{h}_{AB}& =& {\pounds }_{\xi }{g}_{AB}={\nabla }_A{\xi }_B+{\nabla }_B{\xi }_A={\overline{D}}_A{\xi }_B+{\overline{D}}_B{\xi }_A,\hfill \end{array}$$

(70)

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}{h}_{Ap}-{}_{\mathcal{X}}{h}_{Ap}& =& {\pounds }_{\xi }{g}_{Ap}={\nabla }_A{\xi }_p+{\nabla }_p{\xi }_A={\overline{D}}_A{\xi }_p+{\widehat{D}}_p{\xi }_A-{\displaystyle \frac{2}{r}}\left({\overline{D}}_{A}r\right){\xi }_p,\hfill \end{array}$$

(71)

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}{h}_{pq}-{}_{\mathcal{X}}{h}_{pq}& =& {\pounds }_{\xi }{g}_{pq}={\nabla }_p{\xi }_q+{\nabla }_q{\xi }_p={\widehat{D}}_p{\xi }_q+{\widehat{D}}_q{\xi }_p+2r\left({\overline{D}}^{A}r\right){\gamma }_{pq}{\xi }_A\hfill \end{array}$$

(72)

from the formulae summarized in Appendix B of the Part I paper [14]. Here, we consider the Fourier transformation of ${\xi }_a=:{\xi }_A{\left(d{x}^A\right)}_a+{\xi }_p{\left(d{x}^p\right)}_a$ as

$$\begin{array}{c}\hfill {\xi }_A=\sum _{l,m}{\zeta }_{A}S,{\xi }_p=r\sum _{l,m}\left({\zeta }_{\left(e\right)}{\widehat{D}}_{p}S+{\zeta }_{\left(o\right)}{ϵ}_{pq}{\widehat{D}}^{p}S\right).\end{array}$$

(73)

From these Fourier transformation rules, Equations (70)–(72) are given by

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}{h}_{AB}-{}_{\mathcal{X}}{h}_{AB}& =& {\overline{D}}_A{\xi }_B+{\overline{D}}_B{\xi }_A=\sum _{l,m}2{\overline{D}}_{(A}{\zeta }_{B)}S,\hfill \\ \hfill {}_{\mathcal{Y}}{h}_{Ap}-{}_{\mathcal{X}}{h}_{Ap}& =& {\overline{D}}_A{\xi }_p+{\widehat{D}}_p{\xi }_A-{\displaystyle \frac{2}{r}}\left({\overline{D}}_{A}r\right){\xi }_p\hfill \\ & =& r\sum _{l,m}\left[\left({\displaystyle \frac{1}{r}}{\zeta }_A+{\overline{D}}_A{\zeta }_{\left(e\right)}-{\displaystyle \frac{1}{r}}\left({\overline{D}}_{A}r\right){\zeta }_{\left(e\right)}\right){\widehat{D}}_{p}S\right.\hfill \end{array}$$

(74)

$$\begin{array}{ccc}& & \left.+\left({\overline{D}}_A{\zeta }_{\left(o\right)}-{\displaystyle \frac{1}{r}}\left({\overline{D}}_{A}r\right){\zeta }_{\left(o\right)}\right){ϵ}_{pq}{\widehat{D}}^{p}S\right],\hfill \\ \hfill {}_{\mathcal{Y}}{h}_{pq}-{}_{\mathcal{X}}{h}_{pq}& =& {\widehat{D}}_p{\xi }_q+{\widehat{D}}_q{\xi }_p+2r\left({\overline{D}}^{A}r\right){\gamma }_{pq}{\xi }_A\hfill \\ & =& r^2\sum _{l,m}\left[{\displaystyle \frac{4}{r}}\left(-{\displaystyle \frac{l(l+1)}{2}}{\zeta }_{\left(e\right)}+\left({\overline{D}}^{A}r\right){\zeta }_A\right){\displaystyle \frac{1}{2}}{\gamma }_{pq}S\right.\hfill \\ & & +{\displaystyle \frac{2}{r}}{\zeta }_{\left(e\right)}\left({\widehat{D}}_p{\widehat{D}}_{q}S-{\displaystyle \frac{1}{2}}{\gamma }_{pq}{\widehat{D}}^r{\widehat{D}}_{r}S\right)\hfill \end{array}$$

(75)

$$\begin{array}{ccc}& & \left.-{\displaystyle \frac{1}{r}}{\zeta }_{\left(o\right)}2{ϵ}_{r(q}{\widehat{D}}_{p)}{\widehat{D}}^{r}S\right].\hfill \end{array}$$

(76)

Here, we used the property of the mode function S is an eigen function of the Laplacian ${\widehat{D}}_r{\widehat{D}}^r:={\gamma }^{rs}{\widehat{D}}_r{\widehat{D}}_s$ as

$$\begin{array}{c}\hfill {\widehat{D}}_r{\widehat{D}}^{r}S=-l(l+1)S,S=Y_{lm}.\end{array}$$

(77)

3.2.1. $l\ge 2$ Perturbations

For $l\ge 2$ modes, the set of mode functions

$$\begin{array}{c}\hfill \left\{S,{\widehat{D}}_{p}S,{ϵ}_{pq}{\widehat{D}}^{q}S,{\displaystyle \frac{1}{2}}{\gamma }_{pq}S,\left({\widehat{D}}_p{\widehat{D}}_q-{\displaystyle \frac{1}{2}}{\gamma }_{pq}{\widehat{D}}_r{\widehat{D}}^r\right)S,2{ϵ}_{r(p}{\widehat{D}}^r{\widehat{D}}_{q)}S\right\}\end{array}$$

(78)

is a linear-independent set of the tensor field of the rank 0, 1, and 2, even if $S=Y_{lm}$. Then, we may compare Equations (26)–() and Equations (74)–(). As the result of this comparison, we obtain the gauge-transformation rules for the mode functions $\{{\tilde{h}}_{AB},{\tilde{h}}_{\left(e1\right)A},{\tilde{h}}_{\left(o1\right)A},{\tilde{h}}_{\left(e0\right)},{\tilde{h}}_{\left(e2\right)},{\tilde{h}}_{\left(o2\right)}\}$ as

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}\tilde{h}_{AB}-{}_{\mathcal{X}}\tilde{h}_{AB}& =& 2{\overline{D}}_{(A}{\zeta }_{B)},\hfill \end{array}$$

(79)

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}\tilde{h}_{\left(e1\right)A}-{}_{\mathcal{X}}\tilde{h}_{\left(e1\right)A}& =& {\displaystyle \frac{1}{r}}{\zeta }_A+{\overline{D}}_A{\zeta }_{\left(e\right)}-{\displaystyle \frac{1}{r}}\left({\overline{D}}_{A}r\right){\zeta }_{\left(e\right)},\hfill \end{array}$$

(80)

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}\tilde{h}_{\left(o1\right)A}-{}_{\mathcal{X}}\tilde{h}_{\left(o1\right)A}& =& {\overline{D}}_A{\zeta }_{\left(o\right)}-{\displaystyle \frac{1}{r}}\left({\overline{D}}_{A}r\right){\zeta }_{\left(o\right)},\hfill \end{array}$$

(81)

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}\tilde{h}_{\left(e0\right)}-{}_{\mathcal{X}}\tilde{h}_{\left(e0\right)}& =& {\displaystyle \frac{4}{r}}\left(-{\displaystyle \frac{l(l+1)}{2}}{\zeta }_{\left(e\right)}+\left({\overline{D}}^{A}r\right){\zeta }_A\right),\hfill \end{array}$$

(82)

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}\tilde{h}_{\left(e2\right)}-{}_{\mathcal{X}}\tilde{h}_{\left(e2\right)}& =& {\displaystyle \frac{2}{r}}{\zeta }_{\left(e\right)},\hfill \end{array}$$

(83)

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}\tilde{h}_{\left(o2\right)}-{}_{\mathcal{X}}\tilde{h}_{\left(o2\right)}& =& -{\displaystyle \frac{1}{r}}{\zeta }_{\left(o\right)}.\hfill \end{array}$$

(84)

From these gauge-transformation rules, it is well-known that we can construct gauge-invariant variables and the gauge-dependent variables and Conjecture 2.1 is valid for $l\ge 2$ modes.

3.2.2. $l=1$ Perturbations

As shown in the Appendix A in Ref. [14], for $l=1$ modes,

$$\begin{array}{c}\hfill {\widehat{D}}_p{Y}_{1m}\ne 0\ne {ϵ}_{pq}{\widehat{D}}^q{Y}_{1m}\end{array}$$

(85)

but

$$\begin{array}{c}\hfill \left({\widehat{D}}_p{\widehat{D}}_q-{\displaystyle \frac{1}{2}}{\widehat{D}}_r{\widehat{D}}^r\right)Y_{1m}=0=2{ϵ}_{r(p}{\widehat{D}}_{q)}{\widehat{D}}^r{Y}_{1m}.\end{array}$$

(86)

In this case, the Fourier transformation (26)–() of the metric perturbation with $S_{\delta }=S=Y_{lm}$ is given by

$$\begin{array}{ccc}\hfill h_{AB}& =& \sum _m{\tilde{h}}_{AB}S,\hfill \end{array}$$

(87)

$$\begin{array}{ccc}\hfill h_{Ap}& =& \sum _{m}r{\tilde{h}}_{\left(e1\right)A}{\widehat{D}}_{p}S+\sum _{m}r{\tilde{h}}_{\left(o1\right)A}{ϵ}_{pq}{\widehat{D}}^{q}S,\hfill \end{array}$$

(88)

$$\begin{array}{ccc}\hfill h_{pq}& =& \sum _m{\displaystyle \frac{r^2}{2}}{\gamma }_{pq}{\tilde{h}}_{\left(e0\right)}S.\hfill \end{array}$$

(89)

The gauge-transformation rules (74)–(76) for each m mode are given by

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}{\tilde{h}}_{AB}-{}_{\mathcal{X}}{\tilde{h}}_{AB}& =& 2{\overline{D}}_{(A}{\zeta }_{B)},\hfill \end{array}$$

(90)

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}{\tilde{h}}_{\left(e1\right)A}-{}_{\mathcal{X}}{\tilde{h}}_{\left(e1\right)A}& =& {\displaystyle \frac{1}{r}}{\zeta }_A+{\overline{D}}_A{\zeta }_{\left(e\right)}-{\displaystyle \frac{1}{r}}\left({\overline{D}}_{A}r\right){\zeta }_{\left(e\right)},\hfill \end{array}$$

(91)

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}{\tilde{h}}_{\left(o1\right)A}-{}_{\mathcal{X}}{\tilde{h}}_{\left(o1\right)A}& =& {\overline{D}}_A{\zeta }_{\left(o\right)}-{\displaystyle \frac{1}{r}}\left({\overline{D}}_{A}r\right){\zeta }_{\left(o\right)},\hfill \end{array}$$

(92)

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}{\tilde{h}}_{\left(e0\right)}-{}_{\mathcal{X}}{\tilde{h}}_{\left(e0\right)}& =& -{\displaystyle \frac{4}{r}}{\zeta }_{\left(e\right)}+{\displaystyle \frac{4}{r}}\left({\overline{D}}^{A}r\right){\zeta }_A.\hfill \end{array}$$

(93)

Comparing the gauge-transformation rules (79)–(84) for $l\ge 2$ with the gauge-transformation rules (90)–(93), it is easy to find the gauge-transformation rule (83) and (84) for $l\ge 2$ do not appear in the $l=1$ mode gauge transformation. This is due to the fact Equations (86), and the reason why we cannot construct gauge-invariant variables for $l=1$ mode perturbations in a similar method to the $l\ge 2$ case.

3.2.3. $l=0$ Perturbations

The spherical harmonic function $S=Y_{lm}$ is constant when $l=0$. In this case, the Fourier transformation (26)–(28) of the metric perturbation with $S_{\delta }=S=Y_{lm}$ is given by

$$\begin{array}{ccc}\hfill h_{AB}& =& {\tilde{h}}_{AB}S,\hfill \end{array}$$

(94)

$$\begin{array}{ccc}\hfill h_{Ap}& =& 0,\hfill \end{array}$$

(95)

$$\begin{array}{ccc}\hfill h_{pq}& =& {\displaystyle \frac{r^2}{2}}{\gamma }_{pq}{\tilde{h}}_{\left(e0\right)}S.\hfill \end{array}$$

(96)

The gauge-transformation rules (74)–() are given by

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}{\tilde{h}}_{AB}-{}_{\mathcal{X}}{\tilde{h}}_{AB}& =& 2{\overline{D}}_{(A}{\zeta }_{B)}={\overline{D}}_A{\zeta }_B+{\overline{D}}_B{\zeta }_A,\hfill \end{array}$$

(97)

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}{\tilde{h}}_{\left(e0\right)}-{}_{\mathcal{X}}{\tilde{h}}_{\left(e0\right)}& =& {\displaystyle \frac{4}{r}}\left({\overline{D}}^{A}r\right){\zeta }_A.\hfill \end{array}$$

(98)

It is easy to see that many gauge-transformation rules in (79)–(84) are missing. This is due to the fact

$$\begin{array}{c}\hfill {\widehat{D}}_{p}S={ϵ}_{pq}{\widehat{D}}^{q}S=0,\left({\widehat{D}}_q{\widehat{D}}_p-{\displaystyle \frac{1}{2}}{\gamma }_{pq}{\widehat{D}}_r{\widehat{D}}^r\right)S=2{ϵ}_{r(p}{\widehat{D}}_{q)}{\widehat{D}}^{r}S=0\end{array}$$

(99)

This is the reason why we cannot construct gauge-invariant variables for $l=0$ mode perturbations in a similar method to the $l\ge 2$ case.

4. l=1 Odd-Mode Perturbation in the Conventional Approach

As discussed in Sec. 3.2.2, $l=1$ mode perturbations are given by Equations (87)–(89). In particular, among these expressions of the $l=1$ mode perturbation, odd-mode perturbation is given by

$$\begin{array}{ccc}\hfill h_{AB}& =& 0,h_{Ap}=r{\tilde{h}}_{\left(o1\right)A}{ϵ}_{pq}{\widehat{D}}^{q}S,h_{pq}=0.\hfill \end{array}$$

(100)

As the property of the $l=1$ spherical harmonics, we obtain the condition (86). For odd-mode perturbation we have

$$\begin{array}{c}\hfill 2{ϵ}_{r(p}{\widehat{D}}_{q)}{\widehat{D}}^r{Y}_{1m}=0.\end{array}$$

(101)

Here, we apply the notation which is introduced by Equations (A14)–(A15) and we obtain

$$\begin{array}{c}\hfill {\overline{h}}_{AB}=0,{\overline{h}}_A^D=0,{\overline{h}}^{pD}:={\gamma }^{pq}{y}^{DE}r{\tilde{h}}_{\left(o1\right)E}{ϵ}_{qr}{\widehat{D}}^{r}S,{\overline{h}}_p^q=0,{\overline{h}}^{pq}=0,\end{array}$$

(102)

for $l=1$ odd mode perturbations. Through this notation of the metric perturbation, the linear perturbations (A16)–(A18) of the Einstein tensor are given by

$$\begin{array}{ccc}\hfill {}^{\left(1\right)}{G}_A^B& =& 0,\hfill \end{array}$$

(103)

$$\begin{array}{ccc}\hfill {}^{\left(1\right)}{G}_A^q& =& {\displaystyle \frac{1}{2r}}\left[{\overline{D}}_A{\overline{D}}^C{\tilde{h}}_{\left(o1\right)C}-{\overline{D}}_C{\overline{D}}^C{\tilde{h}}_{\left(o1\right)A}-{\displaystyle \frac{2}{r}}\left({\overline{D}}^{C}r\right){\overline{D}}_C{\tilde{h}}_{\left(o1\right)A}\right.\hfill \\ & & +{\displaystyle \frac{3}{r}}\left({\overline{D}}^{C}r\right){\overline{D}}_A{\tilde{h}}_{\left(o1\right)C}-{\displaystyle \frac{1}{r}}\left({\overline{D}}_{A}r\right){\overline{D}}^C{\tilde{h}}_{\left(o1\right)C}+{\displaystyle \frac{1}{2r^2}}\left({\overline{D}}_{C}r\right)\left({\overline{D}}^{C}r\right){\tilde{h}}_{\left(o1\right)A}\hfill \\ & & \left.-{\displaystyle \frac{2}{r^2}}\left({\overline{D}}_{A}r\right)\left({\overline{D}}^{C}r\right){\tilde{h}}_{\left(o1\right)C}+{\displaystyle \frac{3}{2r^2}}{\tilde{h}}_{\left(o1\right)A}\right]{ϵ}^{qr}{\widehat{D}}_{r}S,\hfill \end{array}$$

(104)

$$\begin{array}{ccc}\hfill {}^{\left(1\right)}{G}_p^q& =& -{\displaystyle \frac{1}{2r^2}}{\overline{D}}^C\left(r{\tilde{h}}_{\left(o1\right)C}\right){\gamma }^{qr}2{ϵ}_{s(p}{\widehat{D}}_{r)}{\widehat{D}}^{s}S=0.\hfill \end{array}$$

(105)

The final equality in Equation (105) is due to the property of the spherical harmonics $Y_{lm}$ (101). Furthermore, the expression (104) is gauge invariant under the gauge transformation rule (). On the other hand, the ${\overline{D}}^C\left(r{\tilde{h}}_{\left(o1\right)C}\right)$ in Equation (105) is not gauge invariant as shown in below. However, the gauge-invariance of the linear-order Einstein tensor ${}^{\left(1\right)}{G}_a^b$ is guaranteed by the identity $2{ϵ}_{s(p}{\widehat{D}}_{r)}{\widehat{D}}^{s}S=0$ for the $l=1$ odd-mode perturbations.

For the $l=1$ odd-mode perturbations, the components of the linearized energy-momentum tensor are summarized as

$$\begin{array}{c}\hfill {}^{\left(1\right)}{\mathcal{T}}_A^B=0,{}^{\left(1\right)}{\mathcal{T}}_A^q={\displaystyle \frac{1}{r}}{\tilde{T}}_{\left(o1\right)A}{ϵ}^{qr}{\widehat{D}}_{r}S,{}^{\left(1\right)}{\mathcal{T}}_p^q=0.\end{array}$$

(106)

Then, the linearized Einstein equations for the odd-mode perturbations are given by

$$\begin{array}{ccc}& & {\overline{D}}_A{\overline{D}}^C{\tilde{h}}_{\left(o1\right)C}-{\overline{D}}_C{\overline{D}}^C{\tilde{h}}_{\left(o1\right)A}\hfill \\ & & -{\displaystyle \frac{2}{r}}\left({\overline{D}}^{C}r\right){\overline{D}}_C{\tilde{h}}_{\left(o1\right)A}+{\displaystyle \frac{3}{r}}\left({\overline{D}}^{C}r\right){\overline{D}}_A{\tilde{h}}_{\left(o1\right)C}-{\displaystyle \frac{1}{r}}\left({\overline{D}}_{A}r\right){\overline{D}}^C{\tilde{h}}_{\left(o1\right)C}\hfill \\ & & +{\displaystyle \frac{1}{2r^2}}\left({\overline{D}}_{C}r\right)\left({\overline{D}}^{C}r\right){\tilde{h}}_{\left(o1\right)A}-{\displaystyle \frac{2}{r^2}}\left({\overline{D}}_{A}r\right)\left({\overline{D}}^{C}r\right){\tilde{h}}_{\left(o1\right)C}+{\displaystyle \frac{3}{2r^2}}{\tilde{h}}_{\left(o1\right)A}\hfill \\ & =& 16\pi {\tilde{T}}_{\left(o1\right)A},\hfill \end{array}$$

(107)

Although the equation

$$\begin{array}{c}\hfill {\overline{D}}^C\left(r{\tilde{h}}_{\left(o1\right)C}\right)=0\end{array}$$

(108)

does not appear from Equation (105) due to the fact that $2{ϵ}_{s(p}{\widehat{D}}_{r)}{\widehat{D}}^s{Y}_{1m}=0$, we may use Equation (108) as a gauge condition. As noted in Sec. 3.2.2, the gauge-transformation rule for the variable ${\tilde{h}}_{\left(o1\right)C}$ is given by Equation (92). From the gauge transformation rule (92), we consider the gauge-transformation rule for the left-hand side of Equation (108) as

$$\begin{array}{ccc}\hfill {\overline{D}}^A\left(r{}_{\mathcal{Y}}{\tilde{h}}_{\left(o1\right)A}\right)-{\overline{D}}^A\left(r{}_{\mathcal{X}}{\tilde{h}}_{\left(o1\right)A}\right)& =& r{\overline{D}}^A{\overline{D}}_A{\zeta }_{\left(o\right)}-\left({\overline{D}}^A{\overline{D}}_{A}r\right){\zeta }_{\left(o\right)}.\hfill \end{array}$$

(109)

Using the background Einstein equation (Equation (B67) of Appendix B in Ref. [14]), we consider the equation

$$\begin{array}{c}\hfill r{\overline{D}}^C{\overline{D}}_C{\zeta }_{\left(o\right)}-{\displaystyle \frac{1}{r}}\left(1-\left({\overline{D}}^{C}r\right)\left({\overline{D}}_{C}r\right)\right){\zeta }_{\left(o\right)}=-{\overline{D}}^C\left(r{}_{\mathcal{X}}{\tilde{h}}_{\left(o1\right)C}\right)\end{array}$$

(110)

More explicitly, Equation (110)

$$\begin{array}{c}\hfill -{\partial }_t^2{\zeta }_{\left(o\right)}+f{\partial }_r\left(f{\partial }_r{\zeta }_{\left(o\right)}\right)-{\displaystyle \frac{f}{r^2}}\left(1-f\right){\zeta }_{\left(o\right)}=-{\displaystyle \frac{f}{r}}{\overline{D}}^C\left(r{}_{\mathcal{X}}{\tilde{h}}_{\left(o1\right)C}\right)\end{array}$$

(111)

If we choose ${\zeta }_{\left(o\right)}$ as a special solution to Equation (111) or equivalently Equation (110), Equation (108) is regarded as a gauge condition in the $\mathcal{Y}$-gauge, i.e.,

$$\begin{array}{ccc}\hfill {\overline{D}}^C\left(r{}_{\mathcal{Y}}{\tilde{h}}_{\left(o1\right)C}\right)& =& 0.\hfill \end{array}$$

(112)

We have to emphasize that this is not a complete gauge fixing, since there is a room of the choice of ${\zeta }_{\left(o\right)}$ which satisfy the homogeneous equation of Equation (110):

$$\begin{array}{c}\hfill r{\overline{D}}^C{\overline{D}}_C{\zeta }_{\left(o\right)}-{\displaystyle \frac{1}{r}}\left(1-\left({\overline{D}}^{C}r\right)\left({\overline{D}}_{C}r\right)\right){\zeta }_{\left(o\right)}=0,\end{array}$$

(113)

i.e.,

$$\begin{array}{c}\hfill -{\partial }_t^2{\zeta }_{\left(o\right)}+f{\partial }_r\left(f{\partial }_r{\zeta }_{\left(o\right)}\right)-{\displaystyle \frac{f}{r^2}}\left(1-f\right){\zeta }_{\left(o\right)}=0.\end{array}$$

(114)

Under the gauge choice (112), Equation (107) is given by

$$\begin{array}{ccc}& & \left[{\overline{D}}_C{\overline{D}}^C-{\displaystyle \frac{2}{r^2}}\right]\left(r{\tilde{h}}_{\left(o1\right)A}\right)-{\displaystyle \frac{2}{r^2}}\left({\overline{D}}_{A}r\right)\left({\overline{D}}^{C}r\right)\left(r{\tilde{h}}_{\left(o1\right)C}\right)+{\displaystyle \frac{2}{r}}\left({\overline{D}}^{C}r\right){\overline{D}}_A\left(r{\tilde{h}}_{\left(o1\right)C}\right)\hfill \\ & =& 16\pi r{\tilde{T}}_{\left(o1\right)A}.\hfill \end{array}$$

(115)

Although the variable $r{\tilde{h}}_{\left(o1\right)C}$ is a gauge-dependent variable which is different from the gauge-invariant variable $r{\tilde{F}}_A$, we can obtain the equation (115) if we replace the variable $r{\tilde{F}}_A$ in Equation (47) with the variable $r{\tilde{h}}_{\left(o1\right)C}$.

We also note that the gauge condition (112) coincides with Equation (46), though we have to replace the gauge-invariant variable $r{\tilde{F}}^D$ with the gauge-dependent variable $r{\tilde{h}}_{\left(o1\right)C}$ in Equation (46) to confirm this coincidence. Furthermore, we also take into account the continuity equation (48) of the $l=1$ odd-mode perturbation of the energy-momentum tensor. Thus, the equations to be solved for the $l=1$ odd-mode perturbation are the gauge condition (108), the evolution equation (115), and the continuity equation (48). These equations coincide with Equations (46)–(48) except for the fact that the variable to be obtained is not the gauge-invariant $r{\tilde{F}}^D$ variable but the gauge-dependent variable $r{\tilde{h}}_{\left(o1\right)C}$. Then, through the same logic in Ref. [14], we obtain the solution to the gauge condition (112), the evolution equation (115), and the continuity equation (48) as follows:

$$\begin{array}{ccc}& & 2r{\tilde{h}}_{\left(o1\right)C}{sin}^2\theta {\left(d{x}^C\right)}_{(a}{\left(d\varphi \right)}_{b)}=6M{r}^2\left[\int dr{\displaystyle \frac{a_1(t,r)}{r^4}}\right]{sin}^2\theta {\left(dt\right)}_{(a}{\left(d\varphi \right)}_{b)}\hfill \end{array}$$

$$\begin{array}{ccc}& & +{\pounds }_V{g}_{ab},,\hfill \end{array}$$

(116)

$$\begin{array}{ccc}& & V_a=({\beta }_{1}t+{\beta }_0+W_{\left(o\right)}(t,r))r^2{sin}^2\theta {\left(d\varphi \right)}_a,\hfill \end{array}$$

(117)

Here, we concentrate only on the $m=0$ solution, ${\beta }_1$ and ${\beta }_0$ are constant, and $Z_{\left(o\right)}=rf{\partial }_r{W}_{\left(o\right)}(t,r)$ is an arbitrary function which satisfies the equation

$$\begin{array}{c}\hfill {\partial }_t^2{Z}_{\left(o\right)}-f{\partial }_r\left(f{\partial }_r{Z}_{\left(o\right)}\right)+{\displaystyle \frac{1}{r^2}}f\left[2-3(1-f)\right]Z_{\left(o\right)}=16\pi f^2{\tilde{T}}_{\left(o1\right)r}.\end{array}$$

(118)

$a_1(t,r)$ in Equation (116) is given by Equation (6.44) in Ref. [14], i.e.,

$$\begin{array}{ccc}\hfill a_1(t,r)& =& -{\displaystyle \frac{16\pi }{3M}}{r}^{3}f\int dt{\tilde{T}}_{\left(o1\right)r}+a_{10}\hfill \\ & =& -{\displaystyle \frac{16\pi }{3M}}\int dr{r}^3{\displaystyle \frac{1}{f}}{\tilde{T}}_{\left(o1\right)t}+a_{10},\hfill \end{array}$$

(119)

where $a_{10}$ is the perturbative Kerr parameter.

Note that the formal solution described by Equations (116)–(119) has the same form as the formal solution described by Equations (49)–(52). However, we have to emphasize that the variable ${\mathcal{F}}_{Ap}$ in Equation (49) is gauge invariant as noted by Equation (43) and the vector field $V_a$ in Equation (49) is also gauge invariant. On the other hand, the variable ${\tilde{h}}_{\left(o1\right)C}$ is still gauge-dependent. Actually, there are remaining gauge degrees of freedom of the generator ${\zeta }_{\left(o\right)}$ which satisfy Equation (114) as noted above. This is clear from the fact that the gauge transformation rule (92) with Equation (114) gives a still non-trivial transformation rule. This remaining gauge degree of freedom is so-called “residual gauge”. For this reason, there is a possibility that $V_a$ in Equation (116) includes this residual gauge.

To clarify whether $V_a$ in Equation (116) includes the “residual gauge”, we have to confirm Equation (114). Within our rules to compare our gauge-invariant formulation with a conventional gauge-fixed approach, we regard the terms in $V_a$ in Equation (116) which satisfies Equation (114) is the second-kind gauge degree of freedom and we regard these degrees of freedom as “unphysical degree of freedom.” To clarify this “unphysical degree of freedom,” we introduce the indicator function ${\mathfrak{R}}_{\left(o\right)}[*]$ as

$$\begin{array}{c}\hfill {\mathfrak{R}}_{\left(o\right)}\left[{\zeta }_{\left(o\right)}\right]:=-{\partial }_t^2{\zeta }_{\left(o\right)}+f{\partial }_r\left(f{\partial }_r{\zeta }_{\left(o\right)}\right)-{\displaystyle \frac{f}{r^2}}\left(1-f\right){\zeta }_{\left(o\right)}.\end{array}$$

(120)

If ${\mathfrak{R}}_{\left(o\right)}\left[{\zeta }_{\left(o\right)}\right]=0$, we should regard ${\zeta }_{\left(o\right)}$ is the second-kind gauge and we regard ${\zeta }_{\left(o\right)}$ as “unphysical degree of freedom.”

We can easily confirm that

$$\begin{array}{c}\hfill {\mathfrak{R}}_{\left(o\right)}\left[({\beta }_{1}t+{\beta }_0)r\right]=0.\end{array}$$

(121)

Then, according to our above rule of comparison, we have to conclude that the rigid rotating term $({\beta }_{1}t+{\beta }_0)r$ in $V_a$ in Equation (117) should be regarded as the second-kind gauge degree of freedom which is “unphysical degree of freedom.”

Next, we consider the term $W_{\left(o\right)}(t,r)r$ in Equation (117). In this case, the direct calculation of the indicator ${\mathfrak{R}}_{\left(o\right)}\left[W_{\left(o\right)}(t,r)r\right]$ is useless. However, it is useful to consider the r-derivative of ${\mathfrak{R}}_{\left(o\right)}\left[W_{\left(o\right)}(t,r)r\right]$ and we can show that

$$\begin{array}{ccc}\hfill rf{\partial }_r\left({\displaystyle \frac{1}{r}}{\mathfrak{R}}_{\left(o\right)}\left[W_{\left(o\right)}(t,r)r\right]\right)& =& -{\partial }_t^2{Z}_{\left(o\right)}+f{\partial }_r\left(f{\partial }_r{Z}_{\left(o\right)}\right)-{\displaystyle \frac{f}{r^2}}\left(2-3(1-f)\right)Z_{\left(o\right)},\hfill \end{array}$$

(122)

where we used $Z_{\left(o\right)}=rf{\partial }_r{W}_{\left(o\right)}(t,r)$. Since the left-hand side of Equation (122) coincides with the right-hand side of Equation (118). Through Equation (118), we obtain

$$\begin{array}{ccc}\hfill rf{\partial }_r\left({\displaystyle \frac{1}{r}}{\mathfrak{R}}_{\left(o\right)}\left[W_{\left(o\right)}(t,r)r\right]\right)& =& -16\pi f^2{\tilde{T}}_{\left(o1\right)r}.\hfill \end{array}$$

(123)

or, equivalently,

$$\begin{array}{ccc}\hfill {\mathfrak{R}}_{\left(o\right)}\left[W_{\left(o\right)}(t,r)r\right]& =& -16\pi r\int dr{\displaystyle \frac{f}{r}}{\tilde{T}}_{\left(o1\right)r}.\hfill \end{array}$$

(124)

Then, if we consider the case where ${\tilde{T}}_{\left(o1\right)r}\ne 0$, we have nonvanishing ${\mathfrak{R}}_{\left(o\right)}\left[W_{\left(o\right)}(t,r)r\right]$. This means that the $W_{\left(o\right)}(t,r)$ term does not belong to the second-kind gauge, and we have to regard the degree of freedom $W_{\left(o\right)}(t,r)$ in Equation (117) is a “physical one”.

On the other hand, even in the case where ${\tilde{T}}_{\left(o1\right)r}=0$, Equation (122) is valid. If $W_{\left(o\right)}(t,r)$ belongs to the second-kind gauge, i.e., ${\mathfrak{R}}_{\left(o\right)}\left[W_{\left(o\right)}(t,r)r\right]=0$, Equation (122) yields

$$\begin{array}{c}\hfill -{\partial }_t^2{Z}_{\left(o\right)}+f{\partial }_r\left(f{\partial }_r{Z}_{\left(o\right)}\right)-{\displaystyle \frac{f}{r^2}}\left(2-3(1-f)\right)Z_{\left(o\right)}=0.\end{array}$$

(125)

However, even if $Z_{\left(o\right)}$ is a solution to Equation (125), we cannot directly yield ${\mathfrak{R}}_{\left(o\right)}\left[W_{\left(o\right)}(t,r)r\right]=0$. Therefore, considering the following two sets of function $W_{\left(o\right)}(t,r)r$ as

$$\begin{array}{ccc}\hfill {\mathfrak{G}}_{\left(o\right)}& :=& \left\{\left.r{W}_{\left(o\right)}(t,r)\right|{\mathfrak{R}}_{\left(o\right)}\left[W_{\left(o\right)}(t,r)r\right]=0\right\},\hfill \\ \hfill {\mathfrak{H}}_{\left(o\right)}& :=& \left\{\left.r{W}_{\left(o\right)}(t,r)\right|Z_{\left(o\right)}=rf{\partial }_r{W}_{\left(o\right)}(t,r),\right.\hfill \end{array}$$

(126)

$$\begin{array}{ccc}& & \left.-{\partial }_t^2{Z}_{\left(o\right)}+f{\partial }_r\left(f{\partial }_r{Z}_{\left(o\right)}\right)-{\displaystyle \frac{f}{r^2}}\left(2-3(1-f)\right)Z_{\left(o\right)}=0.\right\},\hfill \end{array}$$

(127)

we obtain the relation

$$\begin{array}{c}\hfill {\mathfrak{G}}_{\left(o\right)}\subset {\mathfrak{H}}_{\left(o\right)}.\end{array}$$

(128)

This indicates that a part of solutions to Equation (125) should be regarded as a gauge degree of freedom of the second kind which is the “unphysical degree of freedom.”

Furthermore, to obtain the explicit solution $W_{\left(o\right)}(t,r)$, we have to solve Equation (118) with appropriate boundary conditions. Equation (118) is an inhomogeneous second-order linear differential equation for $Z_{\left(o\right)}$ and its boundary conditions are adjusted by the homogeneous solutions to Equation (118), i.e., the element of the set of function ${\mathfrak{H}}_{\left(o\right)}$. However, a part of this homogeneous solution to Equation (118) should be regarded as an unphysical degree of freedom in the “complete gauge fixing approach” as mentioned above. Therefore, in the conventional “complete gauge fixing approach”, the boundary conditions for Equation (118) is restricted. In this sense, a conventional “complete gauge-fixing approach” includes a stronger restriction than our proposed gauge-invariant formulation.

5. l=1 Even-Mode Perturbation in the Conventional Approach

As discussed in Sec. 3.2.2, $l=1$ mode perturbations are given by Equations (87)–(89). In particular, among these expressions of the $l=1$ mode perturbations, even-mode perturbation is given by

$$\begin{array}{c}\hfill h_{AB}={\tilde{h}}_{AB}S,h_{Ap}=r{\tilde{h}}_{\left(e1\right)A}{\widehat{D}}_{p}S,h_{pq}={\displaystyle \frac{r^2}{2}}{\gamma }_{pq}{\tilde{h}}_{\left(e0\right)}S.\end{array}$$

(129)

As the property of the $l=1$ spherical harmonics, we obtain the condition (86). For even-mode perturbation, we have

$$\begin{array}{c}\hfill \left({\widehat{D}}_p{\widehat{D}}_q-{\displaystyle \frac{1}{2}}{\gamma }_{pq}{\widehat{D}}_r{\widehat{D}}^r\right)Y_{1m}=0.\end{array}$$

(130)

The gauge-transformation rules for the $l=1$ even-mode perturbations are given by Equations (90), (91), and (93). Inspecting the gauge-transformation rules (91) and (93), we define the variable $\tilde{H}$ by

$$\begin{array}{ccc}\hfill \tilde{H}& :=& {\tilde{h}}_{\left(e0\right)}-4\left({\overline{D}}^{A}r\right){\tilde{h}}_{\left(e1\right)A}.\hfill \end{array}$$

(131)

The gauge-transformation rule of $\tilde{H}$ is given by

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}\tilde{H}-{}_{\mathcal{X}}\tilde{H}& =& -4\left({\overline{D}}^{A}r\right){\overline{D}}_A{\zeta }_{\left(e\right)}-{\displaystyle \frac{4}{r}}\left(1-\left({\overline{D}}^{A}r\right)\left({\overline{D}}_{A}r\right)\right){\zeta }_{\left(e\right)}.\hfill \end{array}$$

(132)

Furthermore, inspecting gauge-transformation rules (90) and () we also define the variable ${\tilde{H}}_{AB}$ by

$$\begin{array}{ccc}\hfill {\tilde{H}}_{AB}& :=& {\tilde{h}}_{AB}-{\overline{D}}_A\left(r{\tilde{h}}_{\left(e1\right)B}\right)-{\overline{D}}_B\left(r{\tilde{h}}_{\left(e1\right)A}\right).\hfill \end{array}$$

(133)

The gauge-transformation of ${\tilde{H}}_{AB}$ is given by

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}{\tilde{H}}_{AB}-{}_{\mathcal{X}}{\tilde{H}}_{AB}& =& -2r{\overline{D}}_A{\overline{D}}_B{\zeta }_{\left(e\right)}+2\left({\overline{D}}_A{\overline{D}}_{B}r\right){\zeta }_{\left(e\right)}.\hfill \end{array}$$

(134)

The definitions (131) and (133) of the variables $\tilde{H}$ and ${\tilde{H}}_{AB}$, respectively, are analogous to the gauge-invariant variable $\tilde{F}$ and ${\tilde{F}}_{AB}$ define by Equations (36) and (38) for $l\ge 2$ modes, respectively. However, the variables $\tilde{H}$ and ${\tilde{H}}_{AB}$ are not gauge invariant. The employment of the variables $\tilde{H}$ and ${\tilde{H}}_{AB}$ corresponds to the gauge fixing of the gauge degree of freedom of the generator ${\zeta }_A$ and the specification of the component ${\tilde{h}}_{\left(e1\right)A}$ of the metric perturbation. On the other hand, the gauge-transformation rules (90), (91), and (93) includes the gauge degree of freedom of the generator ${\zeta }_{\left(e\right)}$. This gauge degree of freedom appears in the gauge-transformation rules (132) and (134).

The linearized Einstein tensor for $l=1$ even modes in terms of ${\tilde{H}}_{AB}$ and $\tilde{H}$ are given as

$$\begin{array}{ccc}\hfill y_{CB}{}^{\left(1\right)}{G}_A^C/S& =& \left[-{\displaystyle \frac{1}{2}}{\overline{D}}_C{\overline{D}}^C+{\displaystyle \frac{3}{r^2}}-{\displaystyle \frac{1}{r}}\left({\overline{D}}^{C}r\right){\overline{D}}_C\right]{\tilde{H}}_{AB}+{\overline{D}}_{(A}{\overline{D}}^C{\tilde{H}}_{B)C}-{\displaystyle \frac{1}{2}}{\overline{D}}_A{\overline{D}}_B{\tilde{H}}_C^C\hfill \\ & & +{\displaystyle \frac{2}{r}}\left({\overline{D}}^{C}r\right)\left[{\overline{D}}_{(A}{\tilde{H}}_{B)C}-{\displaystyle \frac{1}{r}}\left({\overline{D}}_{C}r\right){\tilde{H}}_{AB}\right]-{\displaystyle \frac{1}{2}}{\overline{D}}_A{\overline{D}}_B\tilde{H}-{\displaystyle \frac{1}{r}}\left({\overline{D}}_{(A}r\right){\overline{D}}_{B)}\tilde{H}\hfill \\ & & +{\displaystyle \frac{1}{2}}{y}_{AB}\left[\left[{\overline{D}}_C{\overline{D}}^C+{\displaystyle \frac{3}{r}}\left({\overline{D}}^{C}r\right){\overline{D}}_C\right]\tilde{H}+\left\{{\overline{D}}_D{\overline{D}}^D-{\displaystyle \frac{5}{r^2}}\right\}{\tilde{H}}_C^C\right.\hfill \\ & & +{\displaystyle \frac{1}{r}}\left({\overline{D}}^{C}r\right)\left\{2{\overline{D}}_C{\tilde{H}}_D^D+{\displaystyle \frac{3}{r}}\left({\overline{D}}_{C}r\right){\tilde{H}}_D^D\right\}-{\overline{D}}^C{\overline{D}}^D{\tilde{H}}_{CD}\hfill \\ & & \left.-{\displaystyle \frac{2}{r}}\left({\overline{D}}^{D}r\right)\left\{2{\overline{D}}^C{\tilde{H}}_{CD}+{\displaystyle \frac{1}{r}}\left({\overline{D}}^{C}r\right){\tilde{H}}_{CD}\right\}\right],\hfill \end{array}$$

(135)

$$\begin{array}{ccc}\hfill {}^{\left(1\right)}{G}_A^q& =& \left[{\displaystyle \frac{1}{2r^2}}{\overline{D}}^C{\tilde{H}}_{AC}-{\displaystyle \frac{1}{2r^2}}{\overline{D}}_A{\tilde{H}}_C^C+{\displaystyle \frac{1}{2r^3}}\left({\overline{D}}_{A}r\right){\tilde{H}}_C^C-{\displaystyle \frac{1}{4r^2}}{\overline{D}}_A\tilde{H}\right]{\widehat{D}}^{q}S,\hfill \end{array}$$

(136)

$$\begin{array}{ccc}\hfill {\gamma }_{qr}{}^{\left(1\right)}{G}_p^r& =& \left[+{\displaystyle \frac{1}{4}}{\overline{D}}_C{\overline{D}}^C\tilde{H}+{\displaystyle \frac{1}{2r}}\left({\overline{D}}^{C}r\right){\overline{D}}_C\tilde{H}+{\displaystyle \frac{1}{2}}\left\{+{\overline{D}}_C{\overline{D}}^C+{\displaystyle \frac{1}{r}}\left({\overline{D}}^{C}r\right){\overline{D}}_C-{\displaystyle \frac{1}{r^2}}\right\}{\tilde{H}}_D^D\right.\hfill \\ & & \left.-{\displaystyle \frac{1}{2}}\left\{{\overline{D}}_C+{\displaystyle \frac{2}{r}}\left({\overline{D}}_{C}r\right)\right\}{\overline{D}}_D{\tilde{H}}^{CD}\right]{\gamma }_{pq}S\hfill \\ & & -{\displaystyle \frac{1}{2r^2}}{\tilde{H}}_C^C\left({\widehat{D}}_p{\widehat{D}}_{q}S-{\displaystyle \frac{1}{2}}{\gamma }_{pq}{\widehat{D}}^r{\widehat{D}}_{r}S\right).\hfill \end{array}$$

(137)

Since there is the identity (130) for $l=1$ mode perturbations, the last term of Equation (137) does not appear for $l=1$ even modes. Furthermore, the expressions (135) and (136) are gauge invariant under the gauge-transformation rules (132) and (134). Even in Equation (137) the component ${\gamma }_{qr}{}^{\left(1\right)}{G}_p^r$ is gauge invariant except for the last term in Equation (137). Actually, the variable ${\tilde{H}}_C^C$ is gauge dependent as shown below. However, the gauge-invariance ${\gamma }_{qr}{}^{\left(1\right)}{G}_p^r$ is guaranteed by the identity (130).

On the other hand, the $l=1$ even components of them are given by

$$\begin{array}{ccc}\hfill {}^{\left(1\right)}{\mathcal{T}}_A^B& =& S{\tilde{T}}_A^B,\hfill \end{array}$$

(138)

$$\begin{array}{ccc}\hfill {\gamma }_{qr}{}^{\left(1\right)}{\mathcal{T}}_A^r& =& {\displaystyle \frac{1}{r}}{\tilde{T}}_{\left(e1\right)A}{\widehat{D}}_{q}S,\hfill \end{array}$$

(139)

$$\begin{array}{ccc}\hfill {\gamma }_{qr}{}^{\left(1\right)}{\mathcal{T}}_p^r& =& {\tilde{T}}_{\left(e0\right)}{\displaystyle \frac{1}{2}}{\gamma }_{pq}S+{\tilde{T}}_{\left(e2\right)}\left({\widehat{D}}_p{\widehat{D}}_{q}S-{\displaystyle \frac{1}{2}}{\gamma }_{pq}{\widehat{D}}_r{\widehat{D}}^{r}S\right)\hfill \end{array}$$

(140)

Here again, we note that Equation (130) satisfies as a mathematical identity for $l=1$ mode perturbations. Therefore, the last term of Equation (140) does not appear for $l=1$ even modes.

If equation (130) is not satisfied, we have the equation

$$\begin{array}{c}\hfill \left[-{\displaystyle \frac{1}{2r^2}}{\tilde{H}}_C^C\right]\left({\widehat{D}}_p{\widehat{D}}_{q}S-{\displaystyle \frac{1}{2}}{\gamma }_{pq}{\widehat{D}}^r{\widehat{D}}_{r}S\right)=8\pi {\tilde{T}}_{\left(e2\right)}\left({\widehat{D}}_p{\widehat{D}}_{q}S-{\displaystyle \frac{1}{2}}{\gamma }_{pq}{\widehat{D}}_r{\widehat{D}}^{r}S\right)\end{array}$$

(141)

as one of the components of the linearized Einstein equation. However, Equation (130) implies that this equation is identically satisfied. Therefore, this equation does not restrict ${\tilde{H}}_C^C$ nor ${\tilde{T}}_{\left(e2\right)}$ and there is no other restriction of them. Here, we note that the tensor field ${\tilde{H}}_{AB}$ is not gauge-invariant as shown in Equation (134). Since we may freely choose ${\tilde{H}}_C^C$ and ${\tilde{T}}_{\left(e2\right)}$, we choose ${\tilde{T}}_{\left(e2\right)}=0$ by hand and choose

$$\begin{array}{c}\hfill H_C^C=0.\end{array}$$

(142)

as a gauge condition.

Actually, from the gauge-transformation rule (134), we obtain

$$\begin{array}{ccc}\hfill y^{AB}{}_{\mathcal{Y}}{\tilde{H}}_{AB}-y^{AB}{}_{\mathcal{X}}{\tilde{H}}_{AB}& =& -2r\left[-{\displaystyle \frac{1}{f}}{\partial }_t^2{\zeta }_{\left(e\right)}+{\partial }_r\left(f{\partial }_r{\zeta }_{\left(e\right)}\right)-{\displaystyle \frac{1-f}{r^2}}{\zeta }_{\left(e\right)}\right].\hfill \end{array}$$

(143)

Then, if we choose ${\zeta }_{\left(e\right)}$ as a special solution to the equation

$$\begin{array}{ccc}\hfill {\displaystyle \frac{1}{2r}}{y}^{AB}{}_{\mathcal{X}}{\tilde{H}}_{AB}& =& -{\displaystyle \frac{1}{f}}{\partial }_t^2{\zeta }_{\left(e\right)}+{\partial }_r\left(f{\partial }_r{\zeta }_{\left(e\right)}\right)-{\displaystyle \frac{1-f}{r^2}}{\zeta }_{\left(e\right)},\hfill \end{array}$$

(144)

we may regard that ${}_{\mathcal{Y}}{\tilde{H}}_C^C=0$ in the $\mathcal{Y}$-gauge. We have to note that this gauge fixing (143) is not complete gauge fixing. There is a remaining degree of freedom in the choice of ${\zeta }_{\left(e\right)}$ as the homogeneous solution ${\zeta }_{\left(e\right)h}$ to the wave equation

$$\begin{array}{c}\hfill -{\displaystyle \frac{1}{f}}{\partial }_t^2{\zeta }_{\left(e\right)h}+{\partial }_r\left(f{\partial }_r{\zeta }_{\left(e\right)h}\right)-{\displaystyle \frac{1-f}{r^2}}{\zeta }_{\left(e\right)h}=0.\end{array}$$

(145)

Due to the gauge condition (142), the components (135)–(137) of the linearized Einstein tensor and components (138)–(140) of the linear-order energy-momentum tensor yield the linearized Einstein equations. The $(A,p)$-components (equivalently, $(p,B)$-components) of the linearized Einstein equations are given by

$$\begin{array}{c}\hfill {\overline{D}}^C{\tilde{H}}_{AC}-{\displaystyle \frac{1}{2}}{\overline{D}}_A\tilde{H}=16\pi r{\tilde{T}}_{\left(e1\right)A}.\end{array}$$

(146)

Since the tensor ${\tilde{H}}_{AC}$ is traceless due to the gauge condition (142), equation (146) coincides with Equation (54).

From the trace part of the $(p,q)$-component of the linearized Einstein equation and Equation (146) yields the component (61) of the continuity equation of the linearized energy-momentum tensor.

The trace part of $(A,B)$-component of the linearized Einstein equation with Equation (146) gives

$$\begin{array}{c}\hfill {\overline{D}}_C{\overline{D}}^C\tilde{H}+{\displaystyle \frac{2}{r}}\left({\overline{D}}_{C}r\right){\overline{D}}^C\tilde{H}-{\displaystyle \frac{4}{r^2}}\left({\overline{D}}^{C}r\right)\left({\overline{D}}^{D}r\right){\tilde{H}}_{CD}=16\pi \left({\tilde{T}}_C^C+4\left({\overline{D}}^{D}r\right){\tilde{T}}_{\left(e1\right)D}\right).\end{array}$$

(147)

Since ${\tilde{H}}_{CD}$ is a traceless tensor, Equation (147) coincides with the $l=1$ version of Equation (56) with the source term (57). However, we have to emphasize that $\tilde{H}$ nor ${\tilde{H}}_{CD}$ are not gauge invariant as shown above.

Finally, the traceless part of $(A,B)$-component of the linearized Einstein equation with Equations (146) and (147) yields

$$\begin{array}{ccc}& & \left[-{\overline{D}}_C{\overline{D}}^C-{\displaystyle \frac{2}{r}}\left({\overline{D}}^{C}r\right){\overline{D}}_C+{\displaystyle \frac{4}{r}}\left({\overline{D}}^C{\overline{D}}_{C}r\right)+{\displaystyle \frac{2}{r^2}}\right]{\tilde{H}}_{AB}+{\displaystyle \frac{4}{r}}\left({\overline{D}}^{C}r\right){\overline{D}}_{(A}{\tilde{H}}_{B)C}\hfill \\ & & -{\displaystyle \frac{2}{r}}\left({\overline{D}}_{(A}r\right){\overline{D}}_{B)}\tilde{H}\hfill \\ & =& 16\pi \left[{\tilde{T}}_{AB}-{\displaystyle \frac{1}{2}}{y}_{AB}{\tilde{T}}_C^C-2\left({\overline{D}}_{(A}\left(r{\tilde{T}}_{\left(e1\right)B)}\right)-{\displaystyle \frac{1}{2}}{y}_{AB}{\overline{D}}^C\left(r{\tilde{T}}_{\left(e1\right)C}\right)\right)\right.\hfill \\ & & \left.+2y_{AB}\left({\overline{D}}^{D}r\right){\tilde{T}}_{\left(e1\right)D}\right],\hfill \end{array}$$

(148)

Since the variable $H_{AB}$ is traceless due to the gauge condition (142), Equation (148) coincides with the $l=1$ version of Equation (58) with the source term (59).

In addition to these linearized Einstein equations, the following linearized perturbations of the continuity equation should be satisfied

$$\begin{array}{ccc}& & {\overline{D}}^C{\tilde{T}}_C^B+{\displaystyle \frac{2}{r}}\left({\overline{D}}^{D}r\right){\tilde{T}}_D^B-{\displaystyle \frac{2}{r}}{\tilde{T}}_{\left(e1\right)}^B-{\displaystyle \frac{1}{r}}\left({\overline{D}}^{B}r\right){\tilde{T}}_{\left(e0\right)}=0,\hfill \end{array}$$

(149)

$$\begin{array}{ccc}& & {\overline{D}}^C{\tilde{T}}_{\left(e1\right)C}+{\displaystyle \frac{3}{r}}\left({\overline{D}}^{C}r\right){\tilde{T}}_{\left(e1\right)C}+{\displaystyle \frac{1}{2r}}{\tilde{T}}_{\left(e0\right)}=0.\hfill \end{array}$$

(150)

Together with the gauge condition (142), Equations (146)–() coincide with Equations (53)–(61) with $l=1$ in our gauge-invariant formulation developed in Refs. [12,13,14,15,16]. In Ref. [15], we derived the formal solution to Equations (53)–(61) with $l=1$ as Equation (62) with Equations (63), (64), and $m=0$. Due to the coincidence of the set of equations, this formal solution should be the $m=0$ formal solution to Equations (142), and Equations (146)–(150). Then, as the $m=0$ solution, we obtain

$$\begin{array}{ccc}\hfill {\mathcal{H}}_{ab}& :=& H_{AB}cos\theta {\left(d{x}^A\right)}_{(a}{\left(d{x}^B\right)}_{b)}+{\displaystyle \frac{r^2}{2}}Hcos\theta {\gamma }_{ab}\hfill \\ & =& {\pounds }_V{g}_{ab}-{\displaystyle \frac{16\pi r^2}{3(1-f)}}\left[f^2\left\{{\displaystyle \frac{1+f}{2}}{\tilde{T}}_{rr}+rf{\partial }_r{\tilde{T}}_{rr}-{\tilde{T}}_{\left(e0\right)}-4{\tilde{T}}_{\left(e1\right)r}\right\}{\left(dt\right)}_a{\left(dt\right)}_b\right.\hfill \\ & & +{\displaystyle \frac{2r}{f}}\left\{{\partial }_t{\tilde{T}}_{tt}-{\displaystyle \frac{3f(1-f)}{2r}}{\tilde{T}}_{tr}\right\}{\left(dt\right)}_{(a}{\left(dr\right)}_{b)}\hfill \\ & & +{\displaystyle \frac{r}{f}}\left\{{\partial }_r{\tilde{T}}_{tt}-{\displaystyle \frac{3(1-3f)}{2rf}}{\tilde{T}}_{tt}\right\}{\left(dr\right)}_a{\left(dr\right)}_b)\hfill \\ & & \left.+r^2{\tilde{T}}_{tt}{\gamma }_{ab}\right]cos\theta ,\hfill \end{array}$$

(151)

where ${\gamma }_{ab}={\gamma }_{pq}{\left(d{x}^p\right)}_a{\left(d{x}^q\right)}_b$, $V_a=g_{ab}{V}^b$ is the vector field defined by

$$\begin{array}{c}\hfill V_a:=-r{\partial }_t{\Phi }_{\left(e\right)}cos\theta {\left(dt\right)}_a+({\Phi }_{\left(e\right)}-r{\partial }_r{\Phi }_{\left(e\right)})cos{\left(dr\right)}_a-r{\Phi }_{\left(e\right)}sin\theta {\left(d\theta \right)}_a\end{array}$$

(152)

and ${\Phi }_{\left(e\right)}$ is a solution to the equation

$$\begin{array}{c}-{\displaystyle \frac{1}{f}}{\partial }_t^2{\Phi }_{\left(e\right)}+{\partial }_r\left[f{\partial }_r{\Phi }_{\left(e\right)}\right]-{\displaystyle \frac{1-f}{r^2}}{\Phi }_{\left(e\right)}={\displaystyle \frac{16\pi r}{3(1-f)}}{S}_{\left({\Phi }_{\left(e\right)}\right)},\hfill \end{array}$$

(153)

$$\begin{array}{ccc}& & S_{\left({\Phi }_{\left(e\right)}\right)}={\displaystyle \frac{1}{2}}r{\partial }_t{T}_{tr}-{\displaystyle \frac{1}{2}}r{\partial }_r{\tilde{T}}_{tt}+{\displaystyle \frac{1-4f}{2f}}{\tilde{T}}_{tt}-{\displaystyle \frac{f}{2}}{\tilde{T}}_{rr}-f{\tilde{T}}_{\left(e1\right)r}.\hfill \end{array}$$

(154)

As in the case of the $l=1$ odd-mode perturbation, the tensor field ${\mathcal{F}}_{ab}$ in Equation (62) is gauge invariant. Then, $V_a$ in Equation (62), which defined by Equation (63) is also gauge invariant. Therefore, we regarded $V_a$ in Equation (62) as the first kind gauge. On the other hand, the tensor field ${\mathcal{H}}_{ab}$ in Equation (151) is gauge-dependent in the sense of the second kind due to the gauge-transformation rules (132) and (134) for the components $\tilde{H}$ and ${\tilde{H}}_{AB}$, respectively. Actually, there is a remaining gauge degree of freedom of the generator ${\zeta }_{\left(e\right)}$ which satisfies Equation (145) as noted above. This remaining gauge degree of freedom is so-called “residual gauge”. For this reason, there is a possibility that $V_a$ in Equation (151) includes this “residual gauge”, while $V_a$ in Equation (62) is gauge invariant.

Since we already fixed the gauge-degree of freedom ${\zeta }_{\left(e\right)}$ so that Equation (142) is satisfied, the remaining gauge degree of freedom ${\zeta }_{\left(e\right)h}$ must satisfy the equation (145). Therefore, to clarify whether $V_a$ in Equation (151) includes the “residual gauge,” we have to confirm Equation (145). Within our rules to compare our gauge-invariant formulation with a conventional gauge-fixed approach, we regard the term in $V_a$ in Equation (151) satisfies Equation (145) is the second-kind gauge degree of freedom and we regard this degree of freedom as an “unphysical degree of freedom.” As in the case of the $l=1$ odd-mode perturbations, we introduce the indicator function ${\mathfrak{R}}_{\left(e\right)}[*]$ as

$$\begin{array}{c}\hfill {\mathfrak{R}}_{\left(e\right)}\left[{\zeta }_{\left(e\right)h}\right]:=-{\displaystyle \frac{1}{f}}{\partial }_t^2{\zeta }_{\left(e\right)h}+{\partial }_r\left(f{\partial }_r{\zeta }_{\left(e\right)h}\right)-{\displaystyle \frac{1-f}{r^2}}{\zeta }_{\left(e\right)h}.\end{array}$$

(155)

If ${\mathfrak{R}}_{\left(e\right)}\left[{\zeta }_{\left(e\right)}\right]=0$, we should regard ${\zeta }_{\left(e\right)}$ is the second-kind gauge degree of freedom and we regard ${\zeta }_{\left(e\right)}$ as an “unphysical degree of freedom.”

To clarify the second-kind gauge degree of freedom, we consider the gauge-transformation rule of the tensor field ${\mathcal{H}}_{ab}$ through the gauge-transformation rules (132) and (134) as

$$\begin{array}{ccc}& & {}_{\mathcal{Y}}\mathcal{H}_{ab}-{}_{\mathcal{X}}\mathcal{H}_{ab}\hfill \\ & :=& \left({}_{\mathcal{Y}}{H}_{AB}-{}_{\mathcal{X}}{H}_{AB}\right)cos\theta {\left(d{x}^A\right)}_{(a}{\left(d{x}^B\right)}_{b)}+{\displaystyle \frac{r^2}{2}}\left({}_{\mathcal{Y}}H-{}_{\mathcal{X}}H\right)cos\theta {\gamma }_{ab}\hfill \end{array}$$

$$\begin{array}{ccc}& =& -2r\left[{\partial }_t^2{\zeta }_{\left(e\right)h}-{\displaystyle \frac{f(1-f)}{2r}}{\partial }_r{\zeta }_{\left(e\right)h}+{\displaystyle \frac{f(1-f)}{2r^2}}{\zeta }_{\left(e\right)h}\right]cos\theta {\left(dt\right)}_a{\left(dt\right)}_b\hfill \\ & & -4r\left[{\partial }_t{\partial }_r{\zeta }_{\left(e\right)h}-{\displaystyle \frac{1-f}{2rf}}{\partial }_t{\zeta }_{\left(e\right)h}\right]cos\theta {\left(dt\right)}_{(a}{\left(dr\right)}_{b)}\hfill \\ & & -2r\left[{\partial }_r^2{\zeta }_{\left(e\right)h}+{\displaystyle \frac{1-f}{2rf}}{\partial }_r{\zeta }_{\left(e\right)h}-{\displaystyle \frac{1-f}{2r^{2}f}}{\zeta }_{\left(e\right)h}\right]cos\theta {\left(dr\right)}_a{\left(dr\right)}_b\hfill \end{array}$$

$$\begin{array}{ccc}& & -2r^2\left(f{\partial }_r{\zeta }_{\left(e\right)h}+{\displaystyle \frac{1}{r}}(1-f){\zeta }_{\left(e\right)h}\right)cos\theta {\gamma }_{ab}\hfill \end{array}$$

(156)

$$\begin{array}{ccc}& =& {\pounds }_W{g}_{ab},\hfill \end{array}$$

(157)

where

$$\begin{array}{ccc}\hfill W_a& =& -r{\partial }_t{\zeta }_{\left(e\right)h}cos\theta {\left(dt\right)}_a+\left({\zeta }_{\left(e\right)h}-r{\partial }_r{\zeta }_{\left(e\right)h}\right)cos\theta {\left(dr\right)}_a-r{\zeta }_{\left(e\right)h}sin\theta {\left(d\theta \right)}_a.\hfill \end{array}$$

(158)

Comparing Equation (157) with the generator (158) and Equation (151) with the generator (152), there is the possibility that the Lie derivative term ${\pounds }_V{g}_{ab}$ in the solution (151) is “residual gauge degree of freedom” with the identification

$$\begin{array}{c}\hfill {\Phi }_{\left(e\right)}={\zeta }_{\left(e\right)h}.\end{array}$$

(159)

For this reason, we check the indicator (155). From Equation (153), we obtain

$$\begin{array}{c}\hfill {\mathfrak{R}}_{\left(e\right)}\left[{\Phi }_{\left(e\right)}\right]:={\displaystyle \frac{16\pi r}{3(1-f)}}{S}_{\left({\Phi }_{\left(e\right)}\right)}.\end{array}$$

(160)

Therefore, the variable ${\Phi }_{\left(e\right)}$ should be regarded as the second-kind gauge degree of freedom and is the “unphysical degree of freedom” if $S_{\left({\Phi }_{\left(e\right)}\right)}=0$, while in the non-vacuum case $S_{{\Phi }_{\left(e\right)}}\ne 0$, the indicator (160) yields ${\mathfrak{R}}_{\left(e\right)}\left[{\Phi }_{\left(e\right)}\right]\ne 0$. This means that ${\Phi }_{\left(e\right)}$ is not the second-kind gauge degree of freedom but is the “physical degree of freedom” in the non-vacuum case $S_{{\Phi }_{\left(e\right)}}\ne 0$. Due to this existence of the source term, the identification (160) is impossible in the case of the non-vacuum situation. Rather, this term is gauge-invariant in the sense of the second kind, and we regard this term as a gauge-degree of freedom of the first kind in the non-vacuum case.

Although we have ${\mathfrak{R}}_{\left(e\right)}\left[{\Phi }_{\left(e\right)}\right]=0$ in the case where $S_{{\Phi }_{\left(e\right)}}=0$ and we should regard ${\Phi }_{\left(e\right)}$ is a “gauge degree of freedom of the second kind”, Equation (161) indicates that

$$\begin{array}{c}\hfill {\mathfrak{R}}_{\left(e\right)}\left[{\Phi }_{\left(e\right)}\right]=-{\displaystyle \frac{1}{f}}{\partial }_t^2{\Phi }_{\left(e\right)}+{\partial }_r\left(f{\partial }_r{\Phi }_{\left(e\right)}\right)-{\displaystyle \frac{1-f}{r^2}}{\Phi }_{\left(e\right)}=0.\end{array}$$

(161)

This coincides with the left-hand side of Equation (153). Therefore, as in the case of the $l=1$ odd-mode perturbations, we consider the following two sets of function ${\Phi }_{\left(e\right)}$ as

$$\begin{array}{ccc}\hfill {\mathfrak{G}}_{\left(e\right)}& :=& \left\{\left.{\Phi }_{\left(e\right)}\right|{\mathfrak{R}}_{\left(o\right)}\left[{\Phi }_{\left(e\right)}\right]=0\right\},\hfill \end{array}$$

(162)

$$\begin{array}{ccc}\hfill {\mathfrak{H}}_{\left(e\right)}& :=& \left\{\left.{\Phi }_{\left(e\right)}\right|-{\displaystyle \frac{1}{f}}{\partial }_t^2{\Phi }_{\left(e\right)}+{\partial }_r\left(f{\partial }_r{\Phi }_{\left(e\right)}\right)-{\displaystyle \frac{1-f}{r^2}}{\Phi }_{\left(e\right)}=0\right\}.\hfill \end{array}$$

(163)

From Equation (161), we obtain the relation

$$\begin{array}{c}\hfill {\mathfrak{G}}_{\left(e\right)}={\mathfrak{H}}_{\left(e\right)}.\end{array}$$

(164)

This indicates that any homogeneous solution (without source term) to Equation (153) should be regarded as a gauge degree of freedom of the second kind which is the “unphysical degree of freedom.”

On the other hand, to obtain the explicit solution ${\Phi }_{\left(e\right)}$ in the case $S_{{\Phi }_{\left(e\right)}}\ne 0$, which is regarded as a “physical degree of freedom,” we have to solve Equation (153) with appropriate boundary conditions. Equation (153) is an inhomogeneous second-order linear differential equation for ${\Phi }_{\left(e\right)}$ and its boundary conditions are adjusted by the homogeneous solutions to Equation (153), i.e., the element of the set of function ${\mathfrak{H}}_{\left(e\right)}$. However, any homogeneous solution to Equation (153) should be regarded as an “unphysical degree of freedom” in the “complete gauge fixing approach” as mentioned above. Therefore, in the conventional “complete gauge fixing approach”, we have to impose the boundary conditions for Equation (118) using a homogeneous solution, which is regarded as an “unphysical degree of freedom,” to obtain a “physical solution” to Equation (118) with nonvanishing source term $S_{{\Phi }_{\left(e\right)}}\ne 0$. This situation is a dilemma. In this sense, as in the case of $l=1$ odd-mode perturbations, a conventional “complete gauge-fixing approach” includes a stronger restriction than our proposed gauge-invariant formulation.

6. l=0 Mode Perturbation in the Conventional Approach

6.1. Einstein Equations for $l=0$ Mode Perturbations

Here, we consider the $l=0$ mode perturbations which are described by Equations (94)–(96). Substitution of Equations (94)–(96) into Equations (A16)–(A18), we obtain the $l=0$ mode perturbations of the linearized Einstein tensor. We use the background Einstein equations (Equation (B67) and (B68) of Appendix B in the Part I paper [14]). Then, we obtain the non-trivial components of the $l=0$ mode Einstein equation ${}^{\left(1\right)}{G}_a^b=8\pi {}^{\left(1\right)}{T}_a^b$ are given by

$$\begin{array}{ccc}& & -{\displaystyle \frac{1}{2}}{\overline{D}}_C{\overline{D}}^C{\tilde{h}}_A^B+{\displaystyle \frac{1}{2}}{\overline{D}}^B{\overline{D}}^C{\tilde{h}}_{AC}+{\displaystyle \frac{1}{2}}{\overline{D}}_A{\overline{D}}_C{\tilde{h}}^{BC}-{\displaystyle \frac{1}{2}}{\overline{D}}_A{\overline{D}}^B{\tilde{h}}_C^C\hfill \\ & & -{\displaystyle \frac{1}{r}}\left({\overline{D}}^{C}r\right){\overline{D}}_C{\tilde{h}}_A^B+{\displaystyle \frac{1}{r}}\left({\overline{D}}^{C}r\right){\overline{D}}^B{\tilde{h}}_{AC}+{\displaystyle \frac{1}{r}}\left({\overline{D}}_{C}r\right){\overline{D}}_A{\tilde{h}}^{BC}-{\displaystyle \frac{2}{r^2}}\left({\overline{D}}^{C}r\right)\left({\overline{D}}_{C}r\right){\tilde{h}}_A^B\hfill \\ & & -{\displaystyle \frac{1}{2r}}\left({\overline{D}}^{B}r\right){\overline{D}}_A{\tilde{h}}_{\left(e0\right)}-{\displaystyle \frac{1}{2r}}\left({\overline{D}}_{A}r\right){\overline{D}}^B{\tilde{h}}_{\left(e0\right)}-{\displaystyle \frac{1}{2}}{\overline{D}}_A{\overline{D}}^B{\tilde{h}}_{\left(e0\right)}+{\displaystyle \frac{2}{r^2}}{\tilde{h}}_A^B\hfill \\ & & +y_A^B\left(-{\displaystyle \frac{1}{2}}{\overline{D}}_C{\overline{D}}_D{\tilde{h}}^{CD}+{\displaystyle \frac{1}{2}}{\overline{D}}_C{\overline{D}}^C{\tilde{h}}_D^D-{\displaystyle \frac{2}{r}}\left({\overline{D}}_{D}r\right){\overline{D}}_C{\tilde{h}}^{CD}+{\displaystyle \frac{1}{r}}\left({\overline{D}}^{C}r\right){\overline{D}}_C{\tilde{h}}_D^D\right.\hfill \\ & & -{\displaystyle \frac{1}{r^2}}\left({\overline{D}}_{C}r\right)\left({\overline{D}}_{D}r\right){\tilde{h}}^{CD}+{\displaystyle \frac{3}{2r^2}}\left({\overline{D}}^{D}r\right)\left({\overline{D}}_{D}r\right){\tilde{h}}_C^C-{\displaystyle \frac{3}{2r^2}}{\tilde{h}}_C^C\hfill \\ & & \left.+{\displaystyle \frac{1}{2r^2}}{\tilde{h}}_{\left(e0\right)}+{\displaystyle \frac{3}{2r}}\left({\overline{D}}^{C}r\right){\overline{D}}_C{\tilde{h}}_{\left(e0\right)}+{\displaystyle \frac{1}{2}}{\overline{D}}_C{\overline{D}}^C{\tilde{h}}_{\left(e0\right)}\right)\hfill \end{array}$$

(165)

$$\begin{array}{ccc}& =& 8\pi {\tilde{T}}_A^B,\hfill \\ & & {\overline{D}}_C{\overline{D}}^C{\tilde{h}}_D^D-{\overline{D}}_C{\overline{D}}_D{\tilde{h}}^{CD}-{\displaystyle \frac{2}{r}}\left({\overline{D}}_{C}r\right){\overline{D}}_D{\tilde{h}}^{CD}+{\displaystyle \frac{1}{r}}\left({\overline{D}}^{C}r\right){\overline{D}}_C{\tilde{h}}_D^D\hfill \\ & & +{\displaystyle \frac{1}{2}}{\overline{D}}_C{\overline{D}}^C{\tilde{h}}_{\left(e0\right)}+{\displaystyle \frac{1}{r}}\left({\overline{D}}^{C}r\right){\overline{D}}_C{\tilde{h}}_{\left(e0\right)}\hfill \end{array}$$

$$\begin{array}{ccc}& =& 8\pi {\tilde{T}}_{\left(e0\right)}.\hfill \end{array}$$

(166)

Here, we decomposition of the component ${\tilde{h}}_{AB}$ as

$$\begin{array}{c}\hfill {\tilde{h}}_{AB}=:{\tilde{\mathbb{H}}}_{AB}+{\displaystyle \frac{1}{2}}{y}_{AB}{\tilde{h}}_C^C.\end{array}$$

(167)

In terms of the variables defined by Equation (169), the linearization of the Einstein equations for $l=0$ mode are summarized as follows. The trace of ${}^{\left(1\right)}{G}_A^B=8\pi {}^{\left(1\right)}{T}_A^B$ for $l=0$ mode is given by

$$\begin{array}{ccc}& & -{\displaystyle \frac{1}{r^2}}{\tilde{h}}_C^C-{\displaystyle \frac{2}{r}}\left({\overline{D}}_{C}r\right)\left(+{\overline{D}}_D{\tilde{\mathbb{H}}}^{DC}+{\displaystyle \frac{1}{r}}\left({\overline{D}}_{D}r\right){\tilde{\mathbb{H}}}^{CD}\right)\hfill \\ & & +{\displaystyle \frac{1}{2}}{\overline{D}}_C{\overline{D}}^C{\tilde{h}}_{\left(e0\right)}+{\displaystyle \frac{2}{r}}\left({\overline{D}}_{C}r\right){\overline{D}}^C{\tilde{h}}_{\left(e0\right)}+{\displaystyle \frac{1}{r^2}}{\tilde{h}}_{\left(e0\right)}=8\pi {\tilde{T}}_C^C.\hfill \end{array}$$

(168)

The traceless part of ${}^{\left(1\right)}{G}_A^B=8\pi {}^{\left(1\right)}{T}_A^B$ for $l=0$ mode is given by

$$\begin{array}{ccc}& & -{\displaystyle \frac{1}{2}}\left[{\overline{D}}_C{\overline{D}}^C+{\displaystyle \frac{2}{r}}\left({\overline{D}}^{C}r\right){\overline{D}}_C+{\displaystyle \frac{4}{r^2}}\left({\overline{D}}^{C}r\right)\left({\overline{D}}_{C}r\right)-{\displaystyle \frac{4}{r^2}}\right]{\tilde{\mathbb{H}}}_{AB}\hfill \\ & & +{\overline{D}}_{(A}{\overline{D}}^C{\tilde{\mathbb{H}}}_{B)C}-{\displaystyle \frac{1}{2}}{y}_{AB}{\overline{D}}_C{\overline{D}}_D{\tilde{\mathbb{H}}}^{CD}+{\displaystyle \frac{2}{r}}\left({\overline{D}}^{C}r\right)\left({\overline{D}}_{(A}{\tilde{\mathbb{H}}}_{B)C}-{\displaystyle \frac{1}{2}}{y}_{AB}{\overline{D}}^D{\tilde{\mathbb{H}}}_{DC}\right)\hfill \\ & & +{\displaystyle \frac{1}{r}}\left(\left({\overline{D}}_{(A}r\right){\overline{D}}_{B)}{\tilde{h}}_E^E-{\displaystyle \frac{1}{2}}{y}_{AB}\left({\overline{D}}^{C}r\right){\overline{D}}_C{\tilde{h}}_D^D\right)-{\displaystyle \frac{1}{2}}\left[{\overline{D}}_A{\overline{D}}_B{\tilde{h}}_{\left(e0\right)}-{\displaystyle \frac{1}{2}}{y}_{AB}{\overline{D}}_C{\overline{D}}^C{\tilde{h}}_{\left(e0\right)}\right]\hfill \\ & & -{\displaystyle \frac{1}{r}}\left[\left({\overline{D}}_{(A}r\right){\overline{D}}_{B)}{\tilde{h}}_{\left(e0\right)}-{\displaystyle \frac{1}{2}}{y}_{AB}\left({\overline{D}}^{C}r\right){\overline{D}}_C{\tilde{h}}_{\left(e0\right)}\right]=8\pi \left[{\tilde{T}}_{AB}-{\displaystyle \frac{1}{2}}{y}_{AB}{\tilde{T}}_C^C\right].\hfill \end{array}$$

(169)

Furthermore, ${}^{\left(1\right)}{G}_p^q=8\pi {}^{\left(1\right)}{T}_p^q$ for $l=0$ mode is given by

$$\begin{array}{ccc}& & {\displaystyle \frac{1}{2}}{\overline{D}}_C{\overline{D}}^C{\tilde{h}}_{\left(e0\right)}+{\displaystyle \frac{1}{r}}\left({\overline{D}}^{C}r\right){\overline{D}}_C{\tilde{h}}_{\left(e0\right)}-{\overline{D}}_C{\overline{D}}_D{\tilde{\mathbb{H}}}^{CD}-{\displaystyle \frac{2}{r}}\left({\overline{D}}_{C}r\right){\overline{D}}_D{\tilde{\mathbb{H}}}^{CD}+{\displaystyle \frac{1}{2}}{\overline{D}}_C{\overline{D}}^C{\tilde{h}}_D^D\hfill \\ & =& 8\pi {\tilde{T}}_{\left(e0\right)}.\hfill \end{array}$$

(170)

6.2. Gauge-Fixing for $l=0$ Mode Perturbations

We consider the gauge-transformation rules (97) and (98) of the mode coefficient ${\tilde{h}}_{AB}$ and ${\tilde{h}}_{\left(e0\right)}$. If these gauge-transformation rules are those of the second kind, we should exclude these gauge degrees of freedom through some gauge-fixing procedure, because the degree of freedom of the second-kind gauge is the unphysical degree of freedom.

In the static chart in which the metric $y_{AB}$ is given by Equation (17). Through the static chart (17), the gauge-transformation rule (98) is given by

$$\begin{array}{c}\hfill {}_{\mathcal{Y}}{\tilde{h}}_{\left(e0\right)}-{}_{\mathcal{X}}{\tilde{h}}_{\left(e0\right)}={\displaystyle \frac{4}{r}}f{\zeta }_r.\end{array}$$

(171)

Then we may choose ${\zeta }_r$ so that

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}{\tilde{h}}_{\left(e0\right)}& =& {}_{\mathcal{X}}{\tilde{h}}_{\left(e0\right)}+{\displaystyle \frac{4}{r}}f{\zeta }_r=0,\hfill \end{array}$$

(172)

i.e.,

$$\begin{array}{c}\hfill {\zeta }_r=-{\displaystyle \frac{r}{4f}}{}_{\mathcal{X}}{\tilde{h}}_{\left(e0\right)}.\end{array}$$

(173)

Through this gauge-fixing, we may regard ${\tilde{h}}_{\left(e0\right)}=0$ in the $\mathcal{Y}$-gauge.

As the gauge-fixing for ${\tilde{h}}_{AB}$, we fix the gauge so that ${\tilde{h}}_{AB}$ is traceless. This gauge-fixing so that ${\tilde{h}}_{AB}$ is traceless makes easy to compare with our gauge-invariant expression in Refs. [12,13,14,15,16]. The trace of the gauge transformation rue Equation (97), we obtain

$$\begin{array}{ccc}\hfill y^{AB}{}_{\mathcal{Y}}{\tilde{h}}_{AB}-y^{AB}{}_{\mathcal{X}}{\tilde{h}}_{AB}& =& {}_{\mathcal{Y}}{\tilde{h}}_C^C-{}_{\mathcal{X}}{\tilde{h}}_C^C=2{\overline{D}}^C{\zeta }_C.\hfill \end{array}$$

(174)

In the static chart (17) of the Schwarzschild spacetime, the gauge transformation rule (174) is given by

$$\begin{array}{ccc}\hfill {}_{\mathcal{Y}}{\tilde{h}}_C^C-{}_{\mathcal{X}}{\tilde{h}}_C^C& =& -2f^{-1}{\partial }_t{\zeta }_t+2{\partial }_r\left(f{\zeta }_r\right).\hfill \end{array}$$

(175)

From this gauge-transformation rule and the gauge fixing (173), we may fix the gauge degree of freedom ${\zeta }_t$ so that

$$\begin{array}{ccc}\hfill 0& =& {}_{\mathcal{Y}}{\tilde{h}}_C^C\hfill \\ & =& {}_{\mathcal{X}}{\tilde{h}}_C^C-2f^{-1}{\partial }_t{\zeta }_t+2{\partial }_r\left(-{\displaystyle \frac{r}{4}}{}_{\mathcal{X}}{\tilde{h}}_{\left(e0\right)}\right).\hfill \end{array}$$

(176)

Through this gauge-fixing, we may regard ${\tilde{h}}_C^C=0$. However, we have to emphasize that the degree of freedom of the choice of the generator ${\zeta }_t$ remains the degree of freedom of an arbitrary function of r.

Thus, through the gauge-fixing (173) and (176), we may regard

$$\begin{array}{c}\hfill {\tilde{h}}_{\left(e0\right)}={\tilde{h}}_C^C=0\end{array}$$

(177)

Under the gauge-fixing condition (177) linearization of the Einstein equations for $l=0$ mode are summarized as

$$\begin{array}{c}\hfill -\left({\overline{D}}_{C}r\right){\overline{D}}_D{\tilde{\mathbb{H}}}^{DC}-{\displaystyle \frac{1}{r}}\left({\overline{D}}_{C}r\right)\left({\overline{D}}_{D}r\right){\tilde{\mathbb{H}}}^{CD}=4\pi r{\tilde{T}}_C^C,\end{array}$$

(178)

$$\begin{array}{ccc}& & -{\displaystyle \frac{1}{2}}\left[{\overline{D}}_C{\overline{D}}^C+{\displaystyle \frac{2}{r}}\left({\overline{D}}^{C}r\right){\overline{D}}_C+{\displaystyle \frac{4}{r^2}}\left({\overline{D}}^{C}r\right)\left({\overline{D}}_{C}r\right)-{\displaystyle \frac{4}{r^2}}\right]{\tilde{\mathbb{H}}}_{AB}+{\overline{D}}_{(A}{\overline{D}}^C{\tilde{\mathbb{H}}}_{B)C}\hfill \\ & & -{\displaystyle \frac{1}{2}}{y}_{AB}{\overline{D}}_C{\overline{D}}_D{\tilde{\mathbb{H}}}^{CD}+{\displaystyle \frac{2}{r}}\left({\overline{D}}^{C}r\right)\left(+{\overline{D}}_{(A}{\tilde{\mathbb{H}}}_{B)C}-{\displaystyle \frac{1}{2}}{y}_{AB}{\overline{D}}^D{\tilde{\mathbb{H}}}_{DC}\right)\hfill \\ & =& 8\pi \left[{\tilde{T}}_{AB}-{\displaystyle \frac{1}{2}}{y}_{AB}{\tilde{T}}_C^C\right],\hfill \end{array}$$

(179)

and

$$\begin{array}{c}\hfill -{\overline{D}}_C{\overline{D}}_D{\tilde{\mathbb{H}}}^{CD}-{\displaystyle \frac{2}{r}}\left({\overline{D}}_{C}r\right){\overline{D}}_D{\tilde{\mathbb{H}}}^{CD}=8\pi {\tilde{T}}_{\left(e0\right)}.\end{array}$$

(180)

6.3. Component Expression of the Linearized $l=0$ Gauge-Fixed field equations

Here, we consider the component representations of Equations (178)–(180). To do this, we consider the components of the traceless tensor ${\tilde{\mathbb{H}}}_{AB}$ as

$$\begin{array}{c}\hfill {\tilde{\mathbb{H}}}_{AB}=:X_{\left(e\right)}\left[{\left(dt\right)}_A{\left(dt\right)}_B+f^{-2}{\left(dr\right)}_A{\left(dr\right)}_B\right]+2Y_{\left(e\right)}{\left(dt\right)}_{(A}{\left(dr\right)}_{B)}.\end{array}$$

(181)

Here, we use the background Einstein equation (B65) of Appendix B in the Part I paper [14], i.e.,

$$\begin{array}{c}\hfill {\partial }_{r}f={\displaystyle \frac{1-f}{r}}.\end{array}$$

(182)

Then, we obtain

$$\begin{array}{c}\hfill {\partial }_r^{2}f={\partial }_r\left({\displaystyle \frac{1-f}{r}}\right)=-{\displaystyle \frac{2(1-f)}{r^2}}\end{array}$$

(183)

Through Equations (182) and (183), Equation (178) is given by

$$\begin{array}{c}\hfill rf{\partial }_t{Y}_{\left(e\right)}-rf{\partial }_r{X}_{\left(e\right)}-f{X}_{\left(e\right)}=-4\pi r^2\left({\tilde{T}}_{tt}-f^2{\tilde{T}}_{rr}\right).\end{array}$$

(184)

The $(A,B)=(t,t)$ and $(A,B)=(r,r)$ components of Equation (179) yield the same equation as

$$\begin{array}{c}\hfill rf{\partial }_t{Y}_{\left(e\right)}=4\pi r^2\left({\tilde{T}}_{tt}+f^2{\tilde{T}}_{rr}\right).\end{array}$$

(185)

The $(A,B)=(t,r)$ component of Equation (179) is given by

$$\begin{array}{c}\hfill {\partial }_t{X}_{\left(e\right)}=8\pi rf{\tilde{T}}_{tr}.\end{array}$$

(186)

In terms of the components (181), Equation (180) is given by

$$\begin{array}{c}\hfill -{\partial }_t^2{X}_{\left(e\right)}-f^2{\partial }_r^2{X}_{\left(e\right)}-{\displaystyle \frac{2}{r}}{f}^2{\partial }_r{X}_{\left(e\right)}+2f^2{\partial }_t{\partial }_r{Y}_{\left(e\right)}+{\displaystyle \frac{1}{r}}f(1+f){\partial }_t{Y}_{\left(e\right)}=8\pi f^2{\tilde{T}}_{\left(e0\right)}.\end{array}$$

(187)

On the other hand, the even-mode perturbation of the divergence of the energy-momentum tensor is given by Equations (60) and (60). However, Equation (60) does not appear due to ${\widehat{D}}_{p}S={ϵ}_{pq}{\widehat{D}}^{q}S=0$. Then, the non-trivial $l=0$ components of the divergence of the energy-momentum tensor are given by

$$\begin{array}{c}\hfill {\overline{D}}^C{\tilde{T}}_C^B+{\displaystyle \frac{2}{r}}\left({\overline{D}}^{D}r\right){\tilde{T}}_D^B-{\displaystyle \frac{1}{r}}\left({\overline{D}}^{B}r\right){\tilde{T}}_{\left(e0\right)}=0.\end{array}$$

(188)

The $B=t$ component of Equation (188) is given by

$$\begin{array}{c}\hfill {\partial }_t{\tilde{T}}_{tt}-f^2{\partial }_r{\tilde{T}}_{rt}-{\displaystyle \frac{1}{r}}f(1+f){\tilde{T}}_{rt}=0.\end{array}$$

(189)

and $B=r$ component of Equation (188) is given by

$$\begin{array}{c}\hfill -f{\partial }_t{\tilde{T}}_{tr}+f^3{\partial }_r{\tilde{T}}_{rr}+{\displaystyle \frac{1}{2r}}(1-f){\tilde{T}}_{tt}+{\displaystyle \frac{1}{2r}}{f}^2(3+f){\tilde{T}}_{rr}-{\displaystyle \frac{1}{r}}{f}^2{\tilde{T}}_{\left(e0\right)}=0.\end{array}$$

(190)

Substituting Equation (185) into Equation (184), we obtain

$$\begin{array}{c}\hfill {\partial }_r\left(r{X}_{\left(e\right)}\right)=r{\partial }_r{X}_{\left(e\right)}+X_{\left(e\right)}=8\pi {\displaystyle \frac{r^2}{f}}{\tilde{T}}_{tt}.\end{array}$$

(191)

Substituting Equations (185), (186), and (191) into Equation (187), we obtain

$$\begin{array}{ccc}\hfill 0& =& -f{\partial }_t{\tilde{T}}_{tr}+f^3{\partial }_r{\tilde{T}}_{rr}+{\displaystyle \frac{1}{2r}}(1-f){\tilde{T}}_{tt}+{\displaystyle \frac{1}{2r}}{f}^2(3+f){\tilde{T}}_{rr}-{\displaystyle \frac{1}{r}}{f}^2{\tilde{T}}_{\left(e0\right)}.\hfill \end{array}$$

This coincides with Equation (190) which indicates that Equation (187) does not give any new information other than Equation (190).

Here, we consider the integrability of Equation (186) and (191) as follows:

$$\begin{array}{ccc}\hfill {\partial }_r\left({\partial }_t\left(r{X}_{\left(e\right)}\right)\right)-{\partial }_t\left({\partial }_r\left(r{X}_{\left(e\right)}\right)\right)& =& {\partial }_r\left(8\pi r^{2}f{\tilde{T}}_{tr}\right)-{\partial }_t\left(8\pi {\displaystyle \frac{r^2}{f}}{\tilde{T}}_{tt}\right)\hfill \\ & =& -8\pi r^2\left[+{\partial }_t{\tilde{T}}_{tt}-f^2{\partial }_r{\tilde{T}}_{tr}-{\displaystyle \frac{1}{r}}f(1+f){\tilde{T}}_{tr}\right]\hfill \\ & =& 0.\hfill \end{array}$$

(192)

Here, we used Equation (189) in the last equality. This means that the t-component (189) of the continuity equation guarantees the integrability of Equations (186) and (191). Then, from Equation (191), we may write the solution to Equations (186) and (191) under the integrability condition (189) as

$$\begin{array}{c}\hfill X_{\left(e\right)}={\displaystyle \frac{1}{r}}\left[2M_1+8\pi \int dr{\displaystyle \frac{r^2}{f}}{\tilde{T}}_{tt}\right],\end{array}$$

(193)

where $M_1$ is the constant of integration. This $M_1$ corresponds to the perturbation of the Schwarzschild mass parameter.

On the other hand, in the linearized Einstein equation, there is no equation for ${\partial }_r{Y}_{\left(e\right)}$ which guarantees the integrability condition for Equation (185). However, we may write the solution to the equation (185) as

$$\begin{array}{c}\hfill Y_{\left(e\right)}={\displaystyle \frac{4\pi r}{f}}\int dt\left({\tilde{T}}_{tt}+f^2{\tilde{T}}_{rr}\right)+Y_{\left(e\right)0}\left(r\right).\end{array}$$

(194)

where $Y_{\left(e\right)0}\left(r\right)$ is an arbitrary function of r. There is no equation that determines the arbitrary function $Y_{\left(e\right)0}\left(r\right)$ of r within the linearized Einstein equations.

From the definition (181) of the metric perturbation ${\tilde{\mathbb{H}}}_{AB}$ and the solutions (193) and (194), we obtain

$$\begin{array}{ccc}\hfill {\tilde{\mathbb{H}}}_{AB}& =& {\displaystyle \frac{2}{r}}\left(M_1+4\pi \int dr{\displaystyle \frac{r^2}{f}}{\tilde{T}}_{tt}\right)\left({\left(dt\right)}_A{\left(dt\right)}_B+f^{-2}{\left(dr\right)}_A{\left(dr\right)}_B\right)\hfill \\ & & 2\left[4\pi r\int dt\left({\displaystyle \frac{1}{f}}{\tilde{T}}_{tt}+f{\tilde{T}}_{rr}\right)+Y_{\left(e\right)0}\left(r\right)\right]{\left(dt\right)}_{(A}{\left(dr\right)}_{B)}.\hfill \end{array}$$

(195)

Here, the components of the energy-momentum tensor satisfy the continuity equations (189) and (190).

To interpret the term of $Y_{\left(e\right)0}\left(r\right)$ in Equation (195), we consider the term ${\pounds }_V{g}_{ab}$ with the generator $V_a$ whose components are given by

$$\begin{array}{c}\hfill V_a=V_t\left(r\right){\left(dt\right)}_a.\end{array}$$

(196)

Then, the nonvanishing components of ${\pounds }_V{g}_{ab}$ are given by

$$\begin{array}{ccc}\hfill {\pounds }_V{g}_{tr}& =& f{\partial }_r\left({\displaystyle \frac{1}{f}}{V}_t\right).\hfill \end{array}$$

(197)

Choosing $V_t$ so that

$$\begin{array}{c}\hfill f{\partial }_r\left({\displaystyle \frac{1}{f}}{V}_t\right)=Y_{\left(e\right)0}\left(r\right),V_t=f\int dr{\displaystyle \frac{1}{f}}{Y}_{\left(e\right)0}\left(r\right),\end{array}$$

(198)

we obtain

$$\begin{array}{ccc}\hfill {\pounds }_V{g}_{tr}& =& Y_{\left(e\right)0}\left(r\right).\hfill \end{array}$$

(199)

Then, the solution (195) is given by

$$\begin{array}{ccc}\hfill {\tilde{\mathbb{H}}}_{ab}& =& {\displaystyle \frac{2}{r}}\left(M_1+4\pi \int dr{\displaystyle \frac{r^2}{f}}{\tilde{T}}_{tt}\right)\left({\left(dt\right)}_a{\left(dt\right)}_b+f^{-2}{\left(dr\right)}_a{\left(dr\right)}_b\right)\hfill \\ & & +8\pi r\int dt\left({\displaystyle \frac{1}{f}}{\tilde{T}}_{tt}+f{\tilde{T}}_{rr}\right){\left(dt\right)}_{(a}{\left(dr\right)}_{b)}+{\pounds }_V{g}_{ab},\hfill \end{array}$$

(200)

where

$$\begin{array}{c}\hfill V_a=\left(f\int dr{\displaystyle \frac{1}{f}}{Y}_{\left(e\right)0}\left(r\right)\right){\left(dt\right)}_a.\end{array}$$

(201)

As noted just after Equation (176), there is still the remaining degree of freedom of gauge whose generator is given by

$$\begin{array}{c}\hfill {\zeta }_a={\zeta }_t\left(r\right)d{t}_a,\end{array}$$

(202)

where ${\zeta }_t\left(r\right)$ is an arbitrary function of r. Therefore, we may regard the degree of freedom of $V_a$ given in Equation (201) can be eliminated as the second-kind gauge degree of freedom (202). According to our rule of the comparison between our gauge-invariant formulation and a “conventional complete gauge-fixing approach,” we have to regard that the degree of freedom of $V_a$ given in Equation (201) is “unphysical.”

Thus, as the “physical solution” for $l=0$-mode linearized Einstein equation based of a “conventional complete gauge-fixing approach,” we obtain

$$\begin{array}{ccc}\hfill h_{ab}& =& {\displaystyle \frac{2}{r}}\left(M_1+4\pi \int dr{\displaystyle \frac{r^2}{f}}{\tilde{T}}_{tt}\right)\left({\left(dt\right)}_a{\left(dt\right)}_b+f^{-2}{\left(dr\right)}_a{\left(dr\right)}_b\right)\hfill \\ & & +8\pi r\int dt\left({\displaystyle \frac{1}{f}}{\tilde{T}}_{tt}+f{\tilde{T}}_{rr}\right){\left(dt\right)}_{(a}{\left(dr\right)}_{b)}.\hfill \end{array}$$

(203)

This coincides with the solution obtained in Refs. [12,15] except for the terms of the Lie derivative of the background metric $g_{ab}$. Therefore, we conclude that the $l=0$ solution except for the terms of the Lie derivative of the background metric $g_{ab}$ obtained in Refs. [12,15] can be also obtained as (203) through a complete gauge-fixing approach.

The difference between our solution in Refs. [12,15] is in the term of the Lie derivative of the background metric. As shown above, according to our rule of comparison between our gauge-invariant formulation and a “conventional complete gauge-fixing approach,” all terms of the Lie derivative of the background metric should be regarded as the second-kind gauge degree of freedom, and these are “unphysical degree of freedom” in the above solution (203). On the other hand, in our gauge-invariant formulation developed in Refs. [12,13,14,15,16], $l=0,1$-mode metric perturbations are given in a gauge-invariant form and the terms of the Lie derivative of the background metric is included in these gauge-invariant variables. Since the second-kind gauge degree of freedom is completely excluded in our gauge-invariant formulation, we cannot regard these terms of the Lie derivative of the background metric as an “unphysical degree of freedom.” Therefore, we regard the terms of the Lie derivative of the background metric in these gauge-invariant variables as first-kind gauges that have some physical meaning. This point is the essential difference of the solutions in Refs. [12,13,14,15,16] and the above solutions based on a “conventional complete gauge-fixing approach.” This difference leads to the confusion of the interpretation of the perturbative Tolman Bondi solution, as shown in the next subsection.

6.4. Comparing with Lemaître-Olman-Bondi Solution

6.4.1. Perturbative expression of the LTB solution on Schwarzschild background spacetime

Here, we consider the Lemaître-Tolman-Bondi (LTB) solution [48] which is an exact solution to the Einstein equation with the matter field

$$\begin{array}{c}\hfill T_{ab}=\rho u_a{u}_b,u_a=-{\left(d\tau \right)}_a,\end{array}$$

(204)

and the metric

$$\begin{array}{ccc}\hfill g_{ab}& =& -{\left(d\tau \right)}_a{\left(d\tau \right)}_b+{\displaystyle \frac{{\left({\partial }_{R}r\right)}^2}{1+f\left(R\right)}}{\left(dR\right)}_a{\left(dR\right)}_b+r^2{\gamma }_{ab},r=r(\tau ,R).\hfill \end{array}$$

(205)

This solution is a spherically symmetric solution to the Einstein equation. The function $r=r(\tau ,R)$ satisfies the differential equation

$$\begin{array}{c}\hfill {\left({\partial }_{\tau }r\right)}^2={\displaystyle \frac{F\left(R\right)}{r}}+f\left(R\right).\end{array}$$

(206)

Here, we note that $F\left(R\right)$ is an arbitrary function of R representing the dust matter’s initial distribution. $f\left(R\right)$ is also an arbitrary function of R that represents the initial distribution of the energy of the dust field in the Newtonian sense. The solution to Equation (206) is given in the three cases

(i) 

$f\left(R\right)>0$ :

$$\begin{array}{c}\hfill r={\displaystyle \frac{F\left(R\right)}{2f\left(R\right)}}(cosh\eta -1),{\tau }_0\left(R\right)-\tau ={\displaystyle \frac{F\left(R\right)}{2f{\left(R\right)}^{3/2}}}(sinh\eta -\eta ),\end{array}$$

(207)

(ii) 

$f\left(R\right)<0$ :

$$\begin{array}{c}\hfill r={\displaystyle \frac{F\left(R\right)}{-2f\left(R\right)}}(1-cos\eta ),{\tau }_0\left(R\right)-\tau ={\displaystyle \frac{F\left(R\right)}{2{(-f\left(R\right))}^{3/2}}}(\eta -sin\eta ),\end{array}$$

(208)

(iii) 

$f\left(R\right)=0$ :

$$\begin{array}{c}\hfill r={\left({\displaystyle \frac{9F\left(R\right)}{4}}\right)}^{1/3}{\left[{\tau }_0\left(R\right)-\tau \right]}^{2/3}.\end{array}$$

(209)

The energy density $\rho$ is given by

$$\begin{array}{c}\hfill 8\pi \rho ={\displaystyle \frac{{\partial }_{R}F}{\left({\partial }_{R}r\right)r^2}}.\end{array}$$

(210)

The LTB solution includes the three arbitrary functions $f\left(R\right)$, $F\left(R\right)$, and ${\tau }_0\left(R\right)$.

Here, we consider the vacuum case $\rho =0$. In this case, from Equation (210), we have

$$\begin{array}{c}\hfill {\partial }_{R}F=0,F=2M,\end{array}$$

(211)

where M is the Schwarzschild metric with the mass parameter M. Furthermore, we consider the case $f\left(R\right)=0$. Here, we chose ${\tau }_0=R$, i.e., ${\partial }_R{\tau }_0=1$. In this case, Equation (209) yields

$$\begin{array}{ccc}\hfill {\left(dR\right)}_a& =& {\left(d\tau \right)}_a+{\left({\displaystyle \frac{2M}{r}}\right)}^{-1/2}{\left(dr\right)}_a.\hfill \end{array}$$

(212)

Then, the metric (205) is given by

$$\begin{array}{ccc}\hfill g_{ab}& =& -{\left(d\tau \right)}_a{\left(d\tau \right)}_b+{\left({\partial }_{R}r\right)}^2{\left(dR\right)}_a{\left(dR\right)}_b+r^2{\gamma }_{ab}\hfill \\ & =& -f{\left(dt\right)}_a{\left(dt\right)}_b+f^{-1}{\left(dr\right)}_a{\left(dr\right)}_b+r^2{\gamma }_{ab},f=1-{\displaystyle \frac{2M}{r}},\hfill \end{array}$$

(213)

where ${\left(dt\right)}_a$ is defined by

$$\begin{array}{c}\hfill {\left(dt\right)}_a:={\left(d\tau \right)}_a-f^{-1}{\left(1-f\right)}^{1/2}{\left(dr\right)}_a.\end{array}$$

(214)

We also note that the degree of freedom of the choice of the coordinate R is completely fixed through the choice ${\tau }_0=R$, though there remains a degree of freedom of the choice of $R=R\left(\tilde{R}\right)$ in the exact solution (205).

Now, we consider the perturbation of the Schwarzschild spacetime which is derived by the exact LTB solution (205) so that

$$\begin{array}{ccc}\hfill F\left(R\right)& =& 2\left[M+ϵm_1\left(R\right)\right]+O\left({ϵ}^2\right),\hfill \end{array}$$

(215)

$$\begin{array}{ccc}\hfill f\left(R\right)& =& 0+ϵf_1\left(R\right)+O\left({ϵ}^2\right),\hfill \end{array}$$

(216)

$$\begin{array}{ccc}\hfill {\tau }_0\left(R\right)& =& R+ϵ{\tau }_1\left(R\right)+O\left({ϵ}^2\right).\hfill \end{array}$$

(217)

Through these perturbations (215)–(217), we consider the perturbative expansion of the function r which is determined by Equation (206):

$$\begin{array}{c}\hfill r(\tau ,R)=r_s(\tau ,R)+ϵr_1(\tau ,R)+O\left({ϵ}^2\right).\end{array}$$

(218)

Here, the function $r_s(\tau ,R)$ is given by Equation (209), i.e.,

$$\begin{array}{c}\hfill r_s(\tau ,R)=r(\tau ,R)={\left({\displaystyle \frac{9M}{2}}\right)}^{1/3}{\left[R-\tau \right]}^{2/3}.\end{array}$$

(219)

In Equations (217) and (219), we chose the background value of the function ${\tau }_0\left(R\right)$ to be R.

Through this perturbative expansion, we evaluate $O\left({ϵ}^1\right)$ perturbation of Equation (206) through Equation (212) as

$$\begin{array}{c}\hfill {(1-f)}^{1/2}\left({\partial }_{\tau }{r}_1\right)+{\displaystyle \frac{m_1\left(R\right)}{r}}-{\displaystyle \frac{M}{r^2}}{r}_1+{\displaystyle \frac{1}{2}}{f}_1\left(R\right)=0,\end{array}$$

(220)

where we used Equation (212) and the replacement $r_s\to r$. The solution to Equation (222) is given by

$$\begin{array}{ccc}\hfill r_1& =& {\left({\displaystyle \frac{M}{6}}\right)}^{1/3}{\displaystyle \frac{m_1\left(R\right)}{M}}{\left[R-\tau \right]}^{2/3}-{\displaystyle \frac{3}{20}}{\left({\displaystyle \frac{6}{M}}\right)}^{1/3}{f}_1\left(R\right){\left[R-\tau \right]}^{+4/3}\hfill \\ & & +B\left(R\right){\left[R-\tau \right]}^{-1/3}.\hfill \end{array}$$

(221)

From the comparison with Equation (219), $B\left(R\right)$ is the perturbation of the ${\tau }_1\left(R\right)$ as ${\tau }_0\left(R\right)=R+{\tau }_1\left(R\right)$ in the exact solution (207)–(209). Furthermore, the solution (221) can be also derived from the exact solution (207)–(209). From Equation (210), the perturbative dust energy density given by

$$\begin{array}{c}\hfill 8\pi \rho ={\displaystyle \frac{2{\partial }_R{m}_1\left(R\right)}{\left({\partial }_{R}r\right)r^2}}.\end{array}$$

(222)

Through the perturbative solution (221), the metric (205) is given by

$$\begin{array}{ccc}\hfill g_{ab}& =& -{\left(d\tau \right)}_a{\left(d\tau \right)}_b+{\left({\partial }_{R}r\right)}^2{\left(dR\right)}_a{\left(dR\right)}_b+r^2{\gamma }_{ab}\hfill \\ & & +ϵ\left[\left(2\left({\partial }_R{r}_1\right)-f_1\left({\partial }_{R}r\right)\right)\left({\partial }_{R}r\right){\left(dR\right)}_a{\left(dR\right)}_b+2r{r}_1{\gamma }_{ab}\right]+O\left({ϵ}^2\right)\hfill \\ & =:& g_{ab}^{\left(0\right)}+ϵ{}_{\mathcal{X}}{h}_{ab}+O\left({ϵ}^2\right).\hfill \end{array}$$

(223)

As shown in Equation (213), the background metric $g_{ab}^{\left(0\right)}$ is given by the Schwarzschild metric in the static chart. On the other hand, the linear order perturbation ${}_{\mathcal{X}}{h}_{ab}$ (in the gauge ${\mathcal{X}}_{ϵ}$) is given by

$$\begin{array}{c}\hfill {}_{\mathcal{X}}{h}_{ab}:=\left(2\left({\partial }_R{r}_1\right)-f_1\left(R\right)\left({\partial }_{R}r\right)\right)\left({\partial }_{R}r\right){\left(dR\right)}_a{\left(dR\right)}_b+2r{r}_1{\gamma }_{ab}.\end{array}$$

(224)

Here, we fixed the second-kind gauge so that

$$\begin{array}{c}\hfill {\mathcal{X}}_{ϵ}:(\tau ,R,\theta ,\varphi )\in {\mathcal{M}}_{ph}\mapsto (\tau ,R,\theta ,\varphi )\in \mathcal{M}.\end{array}$$

(225)

Through this second-kind gauge choice and the choice of the background radial coordinate R in Equation (213), the degree of freedom of the choice of the radial coordinate $R=R\left(\tilde{R}\right)$ is also completely fixed even in the linearized version of the exact solution (224), which though there remains a degree of freedom of the choice of $R=R\left(\tilde{R}\right)$ in the exact solution (205).

Of course, if we employ the different gauge choice ${\mathcal{Y}}_{ϵ}$ from the above gauge-choice ${\mathcal{X}}_{ϵ}$, we obtain the different expression of the metric perturbation ${}_{\mathcal{Y}}{h}_{ab}$. Actually, we may choose ${\mathcal{Y}}_{ϵ}$ as the identification of

$$\begin{array}{c}\hfill {\mathcal{Y}}_{ϵ}:(\tau +ϵ{\xi }^{\tau }(\tau ,R),R+ϵ{\xi }^R(\tau ,R),\theta ,\varphi )\in {\mathcal{M}}_{ph}\mapsto (\tau ,R,\theta ,\varphi )\in \mathcal{M}.\end{array}$$

(226)

In this identification, the metric on ${\mathcal{M}}_{ph}$ pulled back to $\mathcal{M}$ is given by

$$\begin{array}{ccc}\hfill {}_{{\mathcal{Y}}_{ϵ}}g_{ab}& =& {}_{{\mathcal{X}}_{ϵ}}g_{ab}+ϵ{\pounds }_{\xi }{g}_{ab}^{\left(0\right)}+O\left({ϵ}^2\right)\hfill \\ & =& g_{ab}^{\left(0\right)}+ϵ\left({}_{\mathcal{X}}{h}_{ab}+{\pounds }_{\xi }{g}_{ab}^{\left(0\right)}\right)+O\left({ϵ}^2\right),\hfill \end{array}$$

(227)

where ${\xi }^a={\xi }^{\tau }{\left({\partial }_{\tau }\right)}^a+{\xi }^R{\left({\partial }_R\right)}^a$ is the generator of second-kind gauge transformation ${\mathcal{X}}_{ϵ}\to {\mathcal{Y}}_{ϵ}$.

6.4.2. Expression of the Perturbative LTB Solution in Static Chart

Here, we consider the expression of the linear perturbation ${}_{\mathcal{X}}{h}_{ab}$ given by Equation (224). From Equations (212) and (214) with $F=2M$, we obtain

$$\begin{array}{ccc}\hfill {\left(dR\right)}_a& =& {\left(dt\right)}_a+f^{-1}{(1-f)}^{-1/2}{\left(dr\right)}_a,f=1-{\displaystyle \frac{2M}{r}},\hfill \end{array}$$

(238)

$$\begin{array}{ccc}\hfill {\left(d\tau \right)}_a& =& {\left(dt\right)}_a+f^{-1}{(1-f)}^{1/2}{\left(dr\right)}_a.\hfill \end{array}$$

(229)

First, we consider the perturbation of the energy-momentum tensor of the matter field. In the case of the LTB solution, the matter field is characterized by the dust field whose energy-momentum tensor (204) is given by

$$\begin{array}{c}\hfill T_{ab}=\rho u_a{u}_b,u_a=-{\left(d\tau \right)}_a,u^a={\left({\partial }_{\tau }\right)}^a.\end{array}$$

(230)

In our case, the linearized Einstein equation gives Equation (222), i.e.,

$$\begin{array}{c}\hfill 8\pi \rho ={\displaystyle \frac{2{\partial }_R{m}_1\left(R\right)}{\left({\partial }_{R}r\right)r^2}}={\displaystyle \frac{{\partial }_R{m}_1\left(R\right)}{4\pi r^2}}{(1-f)}^{-1/2}.\end{array}$$

(231)

On the other hand, substituting Equation (229) into Equation (230), we obtain

$$\begin{array}{ccc}\hfill T_{ab}& =& \rho {\left(d\tau \right)}_a{\left(d\tau \right)}_b\hfill \\ & =& \rho \left({\left(dt\right)}_a+f^{-1}{(1-f)}^{1/2}{\left(dr\right)}_a\right)\left({\left(dt\right)}_b+f^{-1}{(1-f)}^{1/2}{\left(dr\right)}_b\right)\hfill \\ & =& \rho {\left(dt\right)}_a{\left(dt\right)}_b+\rho {\displaystyle \frac{{(1-f)}^{1/2}}{f}}2{\left(dt\right)}_{(a}{\left(dr\right)}_{b)}+\rho {\displaystyle \frac{1-f}{f^2}}{\left(dr\right)}_a{\left(dr\right)}_b.\hfill \end{array}$$

(232)

Then, we obtain the components of the energy-momentum tensor for the static coordinate $(t,r)$ as

$$\begin{array}{c}\hfill {\tilde{T}}_{tt}=\rho ,{\tilde{T}}_{tr}={\displaystyle \frac{{(1-f)}^{1/2}}{f}}\rho ,{\tilde{T}}_{rr}={\displaystyle \frac{1-f}{f^2}}\rho ,{\tilde{T}}_{\left(e0\right)}=0.\end{array}$$

(233)

From this component, we can confirm the continuity equations (189) and (190).

As derived in Ref. [16], the linearized Tolman-Bondi solution (224) with the Schwarzschild background is given by

$$\begin{array}{ccc}\hfill {}_{\mathcal{X}}{h}_{ab}& =& {\displaystyle \frac{2m_1\left(R\right)}{r}}\left[{\left(dt\right)}_a{\left(dt\right)}_b+{\displaystyle \frac{1}{f^2}}{\left(dr\right)}_a{\left(dr\right)}_b\right]+{\displaystyle \frac{2-f}{f{(1-f)}^{1/2}}}{\displaystyle \frac{m_1\left(R\right)}{r}}2{\left(dt\right)}_{(a}{\left(dr\right)}_{b)}\hfill \\ & & +{\pounds }_{V_{\left(LTB\right)}}{g}_{ab},\hfill \end{array}$$

(234)

where $V_{\left(LTB\right)a}$ are given by

$$\begin{array}{c}\hfill V_{\left(LTB\right)a}:=\left[{(1-f)}^{1/2}{r}_1+{\displaystyle \frac{1}{2}}f\int dt{f}_1\left(R\right)\right]{\left(dt\right)}_a+{\displaystyle \frac{r_1}{f}}{\left(dr\right)}_a.\end{array}$$

(235)

As shown in Ref. [16], the $l=0$ solution (203) with the components (230) of the linearized energy-momentum tensor realize the first line of Equation (234). In this sense, the solutions (203) of the $l=0$ mode perturbations are justified by the LTB solutions.

However, we emphasize that the linearized LTB solution (234) does have the term ${\pounds }_{V_{LTB}}{g}_{ab}$. As noted by Equations (228) and (229), we can always eliminate the terms of the Lie derivative of the background metric as the second-kind gauge-degree of freedom. Therefore, according to our rule of the comparison, the term ${\pounds }_{V_{\left(LTB\right)}}{g}_{ab}$ should be regarded as the second-kind gauge-degree of freedom in the “conventional gauge-fixing approach” discussed in this paper. In this case, we have to regard that the perturbation $f_1\left(R\right)$ of the initial distribution $f\left(R\right)$ of the energy of dust field in Equation (206) is an “unphysical degree of freedom” in spite that the behavior of the LTB solution crucially depends on the signature of the function $f\left(R\right)$ as shown in Equations (207)–(209).

7. Summary and Discussions

In this paper, we have discussed comparison of our gauge-invariant formulation for $l=0,1$ perturbations on the Schwarzschild background spacetime proposed in the series of papers [12,13,14,15,16] and a “conventional complete gauge-fixing approach.” It is well-known that we cannot construct gauge-invariant variables for $l=0,1$ mode perturbations through a similar manner to the construction of the gauge-invariant variables for $l\ge 2$ mode perturbations if we use the decomposition formula (26)–(28) with the spherical harmonic functions $Y_{lm}$ as the scalar harmonics $S_{\delta }$. In our gauge-invariant formulation for $l=0,1$ perturbations on the Schwarzschild background spacetime, we proposed the introduction of the singular harmonic function at once. Due to this, we can construct gauge-invariant variables for $l=0,1$-mode perturbations in a similar manner to the $l\ge 2$ modes of perturbations. After deriving the mode-by-mode perturbative Einstein equations in terms of the gauge-invariant variables, we impose the regularity on the introduced singular harmonics when we solve the derived Einstein equations. Based on this proposal, we derived formal solutions to the $l=0,1$-mode linearized Einstein equations without the specification of the components of the linear perturbation of the energy-momentum tensor [12,14,15]. Our proposal enables us to develop higher-order perturbations of the Schwarzschild spacetime [13]. Furthermore, we check that our derived solutions realized the linearized version of the LTB solution and non-rotating C-metric [15]. In this sense, we conclude that our proposal is physically reasonable. On the other hand, it is often said that “gauge-invariant formulations in general-relativistic perturbations are equivalent to complete gauge-fixing approaches.” For this reason, we check this statement through the comparison of our gauge-invariant formulation and a “conventional complete gauge-fixing approach” in which we use the spherical harmonic functions $Y_{lm}$ as the scalar harmonics $S_{\delta }$ from the starting point.

After reviewing the concept of “gauges” in general relativistic perturbation theories and our proposed gauge-invariant formulation for the $l=0,1$-mode perturbations and our derived $l=0,1$-mode solutions, we considered $l=1$ odd-mode perturbations, $l=1$ even-mode perturbations, and $l=0$ even-mode perturbations, separately. As a result, it is shown that we can actually derive similar solutions even in the treatment of the “conventional complete gauge-fixing approach.” However, we have to emphasize that the derived solutions are slightly different from those derived based on our gauge-invariant formulation from a conceptual point of view.

In the case of $l=1$ odd-mode perturbations, we obtained the formal solution to the linearized Einstein equation through our proposed gauge-invariant formulation. This formal solution includes the term of the Lie derivative of the background metric. In our gauge-invariant formulation, we describe the solutions only through gauge-invariant variables. Therefore, we should regard the term of the Lie derivative of the background metric is gauge invariant. On the other hand, in the conventional gauge-fixing approach where we use the spherical harmonics $Y_{lm}$ from the starting point, we cannot construct gauge-invariant variable for $l=1$ odd-mode perturbations in a similar manner to those for $l\ge 2$ mode perturbations. For this reason, we have to treat gauge-dependent variables for perturbations. However, the linearized Einstein equations and the continuity equations of the linearized energy-momentum tensor for $l=1$ odd-mode perturbations in terms of these gauge-dependent variables have the completely same form that is derived through our gauge-invariant formulation. Therefore, the same formal solutions derived through our gauge-invariant formulation should be the formal solutions to these linearized Einstein equations and the continuity equation of the linearized energy-momentum tensor in terms of gauge-dependent variables. As noted above, we treat gauge-dependent variables in the conventional gauge-fixing approach. Therefore, the above Lie derivative terms of the background metric may include the gauge degree of freedom of the second kind, which should be regarded as an “unphysical degree of freedom.” Checking the residual gauge degree of freedom, we conclude that the above Lie derivative terms of the background metric include the second-kind gauge degree of freedom. However, in our formal solution, there is a variable that should be obtained by solving the $l=1$ Regge-Wheeler equation. We conclude that the solution to this $l=1$ Regge-Wheeler equation is not the gauge degree of freedom of the second kind but a physical degree of freedom. Furthermore, we have to impose appropriate boundary conditions to solve this $l=1$ Regge-Wheeler equation. Since the $l=1$ Regge-Wheeler equation is an inhomogeneous linear second-order partial differential equation, the boundary conditions for an inhomogeneous linear second-order partial differential equation are adjusted by the homogeneous solutions to this linear second-order partial differential equation. According to the check of the residual gauge degree of freedom, we have to conclude that a part of homogeneous solutions to this equation is the gauge degree of freedom of the second kind. We have to eliminate this part from our consideration because this gauge degree of freedom is “unphysical.” This is the restriction of the boundary conditions of the linearized Einstein equations.

In the case of $l=1$ even-mode perturbations, the situation is worse than the $l=1$ odd-mode case. As in the $l=1$ odd-mode case, we obtained the formal solution to the linearized Einstein equation through our proposed gauge-invariant formulation. This formal solution includes the term of the Lie derivative of the background metric, which is gauge-invariant within our proposed gauge-invariant formulation. On the other hand, in the conventional gauge-fixing approach where we use the spherical harmonics $Y_{lm}$ from the starting point, we cannot construct gauge-invariant variable for $l=1$ even-mode perturbations in a similar manner to those for $l\ge 2$ mode perturbations. For this reason, we have to treat gauge-dependent variables for perturbations. As in the $l=1$ odd-mode case, the linearized Einstein equations and the continuity equations of the linearized energy-momentum tensor for $l=1$ even-mode perturbations in terms of these gauge-dependent variables have the completely same form that is derived through our gauge-invariant formulation. Therefore, the same formal solutions derived through our gauge-invariant formulation should be the formal solutions to these linearized Einstein equations and the continuity equation of the linearized energy-momentum tensor in terms of gauge-dependent variables. As noted above, we treat gauge-dependent variables in the conventional gauge-fixing approach as in the $l=1$ odd-mode case. Therefore, the above Lie derivative terms of the background metric may include the second-kind gauge degree of freedom, which should be regarded as an “unphysical degree of freedom.” However, in our formal solution, there is a variable that should be obtained by solving the $l=1$ Zerilli equation. We conclude that the solution to this $l=1$ Zerilli equation is not the gauge degree of freedom of the second kind but the physical degree of freedom in the non-vacuum case. Furthermore, we have to impose appropriate boundary conditions to solve this $l=1$ Zerilli equation in the non-vacuum case. Since the $l=1$ Zerilli equation is an inhomogeneous linear second-order partial differential equation. The boundary conditions for an inhomogeneous linear second-order partial differential equation are adjusted by the homogeneous solutions to this linear second-order partial differential equation. According to the check of the residual gauge degree of freedom, we have to conclude that all homogeneous solutions to this equation are the gauge degree of freedom of the second kind, i.e., all homogeneous solution to the $l=1$ Zerilli equation is “unphysical.” We have to eliminate these homogeneous solutions from our consideration because these are regarded as “unphysical.” Due to this situation, we have to say that we have to impose the boundary conditions using these “unphysical degrees of freedom” to obtain the “physical solution” through $l=1$ Zerilli equation. This is a dilemma.

In the case of $l=0$ even-mode perturbations, we obtain the complete gauge fixed solution. This solution does not include any term of the Lie derivative of the background metric. All terms of the Lie derivative of the background metric are regarded as the gauge degree of freedom of the second kind. On the other hand, the solution that is derived through our proposed gauge-invariant formulation includes the terms of the Lie derivatives of the background metric which is regarded as the first-kind gauge in our gauge-invariant formulation, i.e., these are “physical.” This difference leads a problem when we compare the derived solution with the linearized LTB solution. In the linearized LTB solution, the initial energy distribution $f_1\left(R\right)$ of the dust field in the Newtonian sense is included in the terms of the Lie derivative of the background metric. When we realize this linearized exact solution by the solution obtained by our gauge-invariant formulation, this initial energy distribution is regarded as a physical degree of freedom. On the other hand, when we realize this linearized exact solution through the solution obtained by a conventional complete gauge-fixing approach, we have to regard this initial energy distribution of the dust field as the gauge degree of freedom of the second kind and have to regard it “unphysical.” Since the behavior of the exact LTB solution crucially depends on this initial energy distribution of the dust field, it is not natural to regard that this initial degree of freedom is unphysical.

In summary, we have to conclude that there is a case where the boundary conditions and initial conditions are restricted in the conventional complete gauge-fixing approach where we use the decomposition of the metric perturbation by the spherical harmonics $Y_{lm}$ from the starting point. On the other hand, such a situation does not occur in our proposed gauge-invariant formulation. As a theory of physics, this point should be regarded as the incompleteness of the conventional complete gauge-fixing approach where we use the decomposition of the metric perturbation by the spherical harmonics $Y_{lm}$ from the starting point. This is the main result of this paper.

Let us discuss why such a difference occurs between our proposed gauge-invariant formulation and a conventional complete gauge-fixed approach where we use the decomposition of the metric perturbation by the spherical harmonics $Y_{lm}$ from the starting point. The aim of the introduction of the singular harmonic functions in our proposed gauge-invariant formulation is to increase the degree of freedom to clarify the distinction between the gauge degree of freedom of the second kind and the physical degree of freedom. As emphasized in Sec. 3, if we choose $S=Y_{lm}$, we have ${\widehat{D}}_{p}S$=${ϵ}_{pq}{\widehat{D}}^{q}S$= 0, $\left({\widehat{D}}_p{\widehat{D}}_q-{\displaystyle \frac{1}{2}}{\gamma }_{pq}{\widehat{D}}^r{\widehat{D}}_r\right)S$ =$2{ϵ}_{r(p}{\widehat{D}}_{q)}{\widehat{D}}^{r}S$= 0 for $l=0$ modes and $\left({\widehat{D}}_p{\widehat{D}}_q-{\displaystyle \frac{1}{2}}{\gamma }_{pq}{\widehat{D}}^r{\widehat{D}}_r\right)S$ =$2{ϵ}_{r(p}{\widehat{D}}_{q)}{\widehat{D}}^{r}S$= 0 for $l=1$ modes. This is the essential reason for the fact that we cannot construct gauge-invariant variables for $l=0,1$ mode metric perturbation in a similar manner to the derivation of gauge-invariant variables for $l\ge 2$ modes. Due to these vanishing vector- or tensor-harmonics, the associated mode coefficients do not appear, and we cannot construct gauge-invariant variables for $l=0,1$-modes in a similar manner to the case of $l\ge 2$ modes. In our series of papers [12,13,14,15,16], we regarded that this is due to the lack of the degree of freedom. For this reason, in our proposed gauge-invariant formulation, we introduced singular harmonic functions $k_{\widehat{\Delta }}$ and $k_{\widehat{\Delta +2}m}$ for $l=0$ and $l=1$ modes, respectively. Owing to the introduction of these singular harmonic functions, we have ${\widehat{D}}_{p}S$≠ 0 ≠${ϵ}_{pq}{\widehat{D}}^{q}S$ and $\left({\widehat{D}}_p{\widehat{D}}_q-{\displaystyle \frac{1}{2}}{\gamma }_{pq}{\widehat{D}}^r{\widehat{D}}_r\right)S$ ≠ 0 ≠$2{ϵ}_{r(p}{\widehat{D}}_{q)}{\widehat{D}}^{r}S$ for $l=0$ modes and $\left({\widehat{D}}_p{\widehat{D}}_q-{\displaystyle \frac{1}{2}}{\gamma }_{pq}{\widehat{D}}^r{\widehat{D}}_r\right)S$ ≠ 0 ≠$2{ϵ}_{r(p}{\widehat{D}}_{q)}{\widehat{D}}^{r}S$ for $l=1$ modes. Thanks to these non-vanishing vector- or tensor-harmonics, the associated mode coefficients do appear, and we can construct gauge-invariant variables for $l=0,1$-modes in a similar manner to the case of $l\ge 2$ modes. This leads to the development of the gauge-invariant perturbation theory for all modes [12,14,15,16] and the development of the higher-order gauge-invariant perturbations [13].

As noted in the Part I paper [14], the decomposition using the spherical harmonics $Y_{lm}$ from the starting point corresponds to the imposition of the regular boundary condition on $S^2$ on the metric perturbations from the starting point. In this sense, our introduction of the singular harmonic functions corresponds to the change of the boundary conditions on $S^2$. As shown in the Part I paper [14], if we introduce this change of boundary conditions on $S^2$, we can clearly distinguish the gauge degree of freedom of the second kind and the physical degree of freedom, i.e., we can easily construct gauge-invariant variables for $l=0,1$-mode perturbations. This indicates that the imposition of the boundary conditions on $S^2$ and the construction of the gauge-invariant variables does not commute as the procedure of calculations. This is the appearance of the non-locality of $l=0,1$-mode perturbations as pointed out in Ref. [38]. Due to these reasons, we reached the conceptual difference between the $l=0,1$-mode solutions discussed in this paper. Namely, in the conventional complete gauge-fixing approach where we use the decomposition by the harmonic function $Y_{lm}$ from the starting point, the degree of freedom of the metric perturbations is not sufficient so that we can distinguish the gauge degree of freedom of the second kind, i.e., “unphysical modes” and “physical modes.” In spite of this situation of the lack of degree of freedom, if we dare to carry out the “complete gauge-fixing” as the “elimination of unphysical modes,” this “complete gauge-fixing” leads to the restriction of the boundary conditions and initial conditions as shown in this paper.

On the other hand, in our proposed gauge-invariant formulation, there are no conceptual difficulties, such as the restriction of the boundary conditions and initial conditions pointed out in this paper, due to the sufficient degree of freedom of the metric perturbations. Incidentally, due to this sufficient degree of freedom of the metric perturbations through the introduction of the singular harmonic functions, our proposed gauge-invariant variables are equivalent to variables of the complete gauge-fixing within our proposed formulation in which the degree of freedom of the metric perturbations is sufficiently extended. Furthermore, we can develop higher-order perturbation theory without gauge-ambiguities if we apply our proposal and there are wide applications of the higher-order gauge-invariant perturbations on the Schwarzschild background spacetime as briefly discussed in Ref. [13]. We leave these further developments of the application of our formulation to concrete problems of perturbations on the Schwarzschild background spacetime as future works.