\section{Introduction}

During the last eighty years researches in different areas of physics have found a vast number of applications of the two-photon (TP) states. In astrophysics, a first step was taken by Breit and Teller \cite{Br-Tel} who applied the calculations of Goeppert-Mayer \cite{GM} of the TP emission for the 2S-1S atomic transition showing that this process is the principal cause of the decay of interstellar hydrogen atoms and the continuum radiation from gaseous nebulae. Generalization of these calculations for highly excited levels of hydrogen atoms have been carried out \cite{Sun}. In addition, the processes of TP decay has been considered recently as a possible mechanism for the temperature and polarization anisotropies of the cosmic microwave background radiation \cite{Peeb}.

Two-photon emission is one of the principal tool in experiments testing violations of the Bell inequalities in atomic cascade and nuclei \cite {Asp,Klei,Gis}. These investigations test the fundamental question of quantum entanglement. Additionally, TP transitions have been proposed as a tool to measure weak interaction properties, viz., the problem of parity violations in atomic transitions \cite{Sch,Frat}. The previously indicated applications of the TP states concern interaction of the electromagnetic field in creation processes of TP by atomic systems, as well as, the properties of entanglement of many-particle states.

Landau \cite{Land} and Yang \cite{Yang} have analyzed symmetry properties of TP under spatial rotations and inversions in the center of mass reference frame of the system. Landau and Yang established selection rules for the TP decay of positronium and mesons. It is shown \cite{Land} that the TP annihilation of orthopositronium is forbidden since there is no state of the TP with angular momentum $J=1$ and the decay is primarily into three gamma photons. In addition, Yang \cite{Yang} also points out the kind of correlation that exists between the planes of polarizations of the TP emitted in the decay of the $\pi ^{0}$ meson.

Landau \cite{Land} introduces two-photon spherical states (TPSS) that are eigenfunctions of energy, parity $P$, angular momentum $J$, and its projection $M$ on the fixed quantization axis. Landau analyzes in detail the general transformation properties of TPSS wave functions in momentum space, in the center of mass reference frame. Table I indicates the number $N_{JM}^{(P )}$ of possible TPSS with given angular momentum $J$, projection $M$, and parity $P$.
\begin{table}[th]
\caption{All possible two-photon spherical states $|JMP\rangle$ ($k= 1, 2,\cdots)$}
\label{table:two
\centering
\begin{tabular}{ccc}
\hline\hline
$J$ &  $N_{JM}^{(+)}$ & $N_{JM}^{( -)} $\\[0.5ex] \hline
0 & $1$ & $1$ \\
1 & $-$ & $-$ \\
2$k$&2&1\\
2$k$+1 & 1 & $-$\\[1ex] \hline
\end{tabular}
\end{table}

The explicit expressions for the TPSS appears in the Wick-Jacob review paper \cite{Jacob} as an application of their helicity formalism. They first introduce TP spherical helicity states $\left\vert JM;\lambda \lambda ^{\prime }\right\rangle $ (see Chap. 2) as eigenfunctions of angular momentum $J$, and its projection $M$ and the total helicity of photons. They then show that each Landau TPSS can be represented as a sum or difference of two TP spherical helicity states with opposite total helicities \cite{Jacob, LLIV} (see Table II below).

The one-particle spherical states together with the theory of angular polarization correlations have been developed both in the classical \cite{Jack} and in the quantum theories \cite{Akh}.  The present article is devoted to analyzing the angular distribution and the polarization correlations in the TPSS.

The structure of this article is as follows. In Sec. II, after a short, preliminary consideration of the explicit form of each two-photon Landau spherical states, we analyze their angular distribution in momentum space relative to the quantization axes. In Sec. III, two-photon polarization matrices are derived for each Landau state. The correlations between the polarizations of the photons is considered in Sec. IV. Finally, the conclusions are given in Sec. V.

\section{Angular distributions}

In the center of mass reference frame, the wave function in momentum space of TP depends only on the direction of the relative momentum $\mathbf{n}\left(\theta ,\varphi \right)$. The Landau TPSS states \cite{Jacob,LLIV} in Table II are given in terms  of the two-photon spherical helicity states  $\left\vert JM;\lambda \lambda ^{\prime }\right\rangle $ defined by
\begin{equation}
\left\vert JM;\lambda \lambda ^{\prime }\right\rangle =\sqrt{\frac{2J+1}
4\pi }}D_{\lambda -\lambda ^{\prime },M}^{\left( J\right) }\left( \mathbf{n
\right) \left\vert \lambda ,\lambda ^{\prime }\right\rangle,  \tag{1}
\end{equation
which are eigenfunctions of the angular momentum, the projection of the angular momentum, and the total helicity of the photon pair. In (1), $D_{\lambda M}^{( J)}$ is the Wigner function and $|\lambda,\lambda^{\prime}\rangle \equiv | \textbf{n};\lambda^{\prime} \otimes |-\textbf{n} ;\lambda^{\prime}\rangle$; $(\lambda,\lambda ^{\prime}=\pm 1)$ are eigenfunctions of the total helicity with eigenvalue $\lambda -\lambda ^{\prime }$. In succeeding calculations, we use the expression
\begin{equation}
D_{\lambda M}^{\left( J\right) }\left( \mathbf{n}\right) =e^{iM\varphi
}d_{\lambda M}^{\left( J\right) }\left( \theta \right)  \tag{2}
\end{equation
for the Wigner function \cite{LLIV}. In (1), $J=0, 2, 3\cdots $ and the absence of $J=1$ is a direct consequence of the transversality of the two-photon system in the center of mass reference frame \cite{LLIV}. Each TPSS in Table II is either the sum or the difference of two-photon spherical helicity states given in (1) with opposite total helicity.

\begin{table*}[ht]
\caption{All possible two-photon spherical states classified by angular momentum $J$ and parity $P$}
\centering
\begin{tabular}{ccc}
\hline\hline
$J$  &\ $P = +1$ &\ $P=-1$  \\ [0.5ex]
\hline
0 & $ \frac{1}{\sqrt{2}}(|00;11\rangle + |00;-1-1\rangle)$ & $\frac{1}{\sqrt{2}}(|00;11\rangle - |00;-1-1\rangle)$  \\
1 & $-$ & $-$ \\
 & a. \hspace{.01in}$\frac{1}{\sqrt{2}}(|JM;11\rangle + |JM;-1-1\rangle)$ &   \\
 2$k$& &$\frac{1}{\sqrt{2}}(|JM;11\rangle - |JM;-1-1\rangle)$  \\
 & b. \hspace{.01in} $\frac{1}{\sqrt{2}}(|JM;1-1\rangle + |JM;-11\rangle)$ &  \\
2$k$+1 &  $\frac{1}{\sqrt{2}}(|JM;1-1\rangle - |JM;-11\rangle)$ & $-$  \\ [1ex]
\hline
\end{tabular}
\label{table:Table I}
\end{table*}

There are two axes in these problems, viz., the quantization axis $Z$ that is defined by the source of radiation and $\mathbf{n}$. Yang \cite{Yang} uses plane waves propagating along and opposite the direction of $\mathbf{n}$. If the quantization axis is coincident with $\mathbf{n}$, the
\begin{equation}
d_{\lambda M}^{\left( J\right) }\left( 0\right) =\delta _{\lambda M}  \tag{3}
\end{equation}
and our Table II corresponds to Yang's table.

The angular distribution of the TP radiation is actually the density of registration of photon pairs at given direction $\mathbf{n}$, which is defined by the square of the modulus of the wave functions in Table II. Using (1) for the two-photon spherical helicity states and (2) for the Wigner function, we obtain the angular distribution for Landau's TPSS
\begin{equation}
w_{\left\vert \Lambda \right\vert }^{JM}=\frac{2J+1}{8\pi }\left[ d_{\Lambda
M}^{\left( J\right) 2}+d_{-\Lambda M}^{\left( J\right) 2}\right],   \tag{4}
\end{equation
which depends on the modulus of the total helicity $\Lambda =\lambda -\lambda^{\prime }$, viz., $\left\vert \Lambda \right\vert =0, 2$. Using the Clebsch-Gordan series representation of two Wigner functions product \cite{Varsh}, we obtain the final expression for the angular distribution (4)
\begin{equation}
w_{m}^{JM}\left( \theta \right) =A_{JM}\sum\limits_{n=0}^{J}\left(4n+1\right) \Big{(}
\begin{array}{ccc}
J & J & 2n \\
m & -m &
\end{array
\Big{)}
\Big{(}
\begin{array}{ccc}
J & J & 2n \\
M & -M &
\end{array
\Big{)} P_{2n}\left( \cos \theta \right), \hspace{0.2in} (m=0, 2)  \tag{5}
\end{equation
where $A_{JM}=\left( -1\right) ^{M}(2J+1)/(4\pi)$.

Table III indicates the angular distributions of the TPSS in accordance with Landau's Table I.
\begin{table}[th]
\caption{Angular distribution of the two-photon spherical states}
\label{table:two
\centering
\begin{tabular}{ccc}
\hline\hline
$J$ & \ $P = +1$ & \ $P=-1$ \\[0.5ex] \hline
0 & $\frac{1}{4\pi}$ & $\frac{1}{4\pi}$ \\
1 & $-$ & $-$ \\
 & a.\hspace{.04in}$w_{0}^{JM}$ &  \\
 2$k$ &  &$w_{0}^{JM}$\\
  & b.\hspace{.04in}$w_{2}^{JM}$ &  \\
 2$k$+1  & $w_{2}^{JM}$ & $-$ \\[1ex] \hline
\end{tabular}
\end{table}

\section{Polarization matrices}

We introduce normalized two-photon polarization matrix $\hat{\rho}$ \cite{F1,vem} for each TPSS of Table II via the pure spherical state density operator
\begin{equation}
\left\vert PJM\right\rangle \left\langle PJM\right\vert =w_{\left\vert \Lambda \right\vert }^{JM}\hat{\rho},  \tag{6}
\end{equation}
where $w_{\left\vert \Lambda \right\vert }^{JM}$ is the angular distribution of the radiation that is coincident with the trace of the density operator. Substituting the expression of every state in Table II into (6) and using the factorization via a Kronecker product of the Pauli matrices $\hat{\sigma}_{i}\otimes \hat{\sigma}_{j}^{\prime }$ \cite{Blum},  we obtain the entries in Table IV for the normalized two-photon polarization matrices
\begin{table}[th]
\caption{Two-photon polarization matrices of Landau's states}
\label{table:two
\centering
\begin{tabular}{ccc}
\hline\hline
$J$ & \ $P = +1$ & \ $P=-1$ \\[0.5ex] \hline
0 & $\hat{\rho}^{+}$ & $\hat{\rho}^{-}$ \\
1 & $-$ & $-$ \\
 & a. $\hspace{.04in}$$\hat{\rho}^{+}$ & {--} \\
 2$k$ &  &$\hat{\rho}^{-}$\\
  & b. \hspace{.04in}$\hat{\rho}^{JM}_{e}$ &  \\
 2$k$+1  & $\hat{\rho}_{o}^{JM}$ & $-$ \\[1ex] \hline
\end{tabular}
\end{table}

In Table IV, the two-photon polarization density matrices $\hat{\rho}^{\pm }$
\begin{equation}
\hat{\rho}^{\pm }=\frac{1}{4}\left\{ \hat{I}\otimes I^{\prime }+\hat{\sigma
_{3}\otimes \hat{\sigma}_{3}^{\prime }\pm \left( \hat{\sigma}_{1}\otimes
\hat{\sigma}_{1}^{\prime }-\hat{\sigma}_{2}\otimes \hat{\sigma}_{2}^{\prime
}\right) \right\}   \tag{7}
\end{equation
do not depend on $J$ and $M$ and are associated with even and odd parities of TPSS. They describe fully entangled states of polarization of the photon pair, i.e., the photons are not polarized individually but the two-photon system is polarized completely \cite{vem,F2}). Now $\hat{\rho}_{e,o}^{JM}$

\begin{equation}
\hat{\rho}_{e,o}^{JM}  = \frac{1}{4} \Big{\{} \hat{I}\otimes I^{\prime }+\xi^{JM}\left( \hat{\sigma}_{3}\otimes \hat{I}^{\prime }-\hat{I}\otimes \hat{\sigma}_{3}^{\prime }\right)\\
 -\hat{\sigma}_{3}\otimes \hat{\sigma}_{3}^{\prime }
 \pm \zeta ^{JM}(\hat{\sigma}_{1}\otimes \hat{\sigma}_{1}^{\prime}+\hat{\sigma}_{2}\otimes \hat{\sigma}_{2}^{\prime} \Big{\}} \tag{8}
\end{equation}
depends on $J$ and $M$ and are associated with TPSS with even ($e$) and odd ($o$) angular momentum $J$ $\left( J\geq 2\right)$. In (8), the polarization parameters $\xi ^{JM}$ and $\zeta ^{JM}$ are defined by \begin{equation}
\xi ^{JM}=\frac{2J+1}{8\pi w_{2}}\left[ d_{2M}^{\left( J\right)
2}-d_{-2M}^{\left( J\right) 2}\right] \text{ }  \tag{9.a}
\end{equation
\begin{equation}
\zeta ^{JM}=\frac{2J+1}{4\pi w_{2}}d_{2M}^{\left( J\right) }d_{-2M}^{\left(
J\right) }  \tag{9.b}
\end{equation
that depend on the polar angle $\theta $ formed by the quantization axis $Z$ and $\mathbf{n}$. The matrix (8) describes correlated states of the individually, partially circularly polarized photons wherein
the system of two photons is polarized completely \cite{F1, vem}.

In the case of photons flying parallel to the quantization axes ($\theta =0$), the matrices $\hat{\rho}_{e,o}^{JM}$ describe uncorrelated states of the two fully, circularly polarized photons. Actually, using identity (3), we find $\xi ^{JM}=\delta _{M,2}-\delta _{M,-2}$ and $\zeta ^{JM}=0$ for the polarization parameters (9). Accordingly, the matrices (8) become
\begin{equation*}
\hat{\rho}_{e,o}^{JM}=\left\{
\begin{array}{c}
\left(\hat{I}+\hat{\sigma}_{3}\right) \otimes \left( \hat{I}^{\prime }-\hat{\sigma}_{3}^{\prime }\right) \text{ } \hspace{.15in} (M=2) \\
\left(\hat{I}-\hat{\sigma}_{3}\right) \otimes \left( \hat{I}^{\prime }+\hat{\sigma}_{3}^{\prime }\right). \text{ } \hspace{.15in}(M=-2
\end{array
\right.
\end{equation*}

The next simple case for expression (8) is the case for photons flying perpendicular to the quantization axes ($\theta =\pi/2$). Using the symmetry property of the $d_{M^{\prime }M}^{\left( J\right) }$ function
\begin{equation*}
d_{M^{\prime }M}^{\left( J\right) }\left( \pi -\theta \right) =\left(
-1\right) ^{J-M}d_{-M^{\prime }M}^{\left( J\right) }\left( \theta \right),
\end{equation*
 we find $d_{2M}^{\left( J\right) }\left( \pi/2\right) =\left( -1\right) ^{J-M}d_{-2M}^{\left( J\right) }\left( \pi/2\right) $ and, with the aid of (4) and (9), we obtian
\begin{equation*}
\xi ^{JM}=0,\text{ }\zeta ^{JM}=\left( -1\right) ^{J-M}
\end{equation*
then the polarization matrices (8) reduces to

\begin{equation}
\hat{\rho}_{e}^{JM} = \hat{\rho}_{o}^{JM} = \frac{1}{4}\big{\{}\hat{I}\otimes I^{\prime}-\hat{\sigma}_{3}\otimes \hat{\sigma}_{3}^{\prime}
 +(-1)^{M}(\hat{\sigma}_{1}\otimes \hat{\sigma}_{1}^{\prime}+\hat{\sigma}_{2}\otimes \hat{\sigma}_{2}^{\prime})\big{\}},\tag{10}
\end{equation}
i.e., at the angle $\theta =\pi /2$ the polarization matrix of the two-photon system depends only on the projection $M$ of the angular momentum  and again we come to a fully entangled polarization state of the individual photons when the photon pair is completely polarized.

Finally, using the Clebsch-Gordan series representation of two Wigner functions product \cite{Varsh}, we obtain the expressions for the polarization parameters $\xi^{JM}$ and $\zeta^{JM}$ of (9) in the form
\begin{equation}
\xi ^{JM}\left( \theta \right) =\frac{A_{JM}}{w_{2}}\sum\limits_{n=1}^{J}\left( 4n-1\right)
\Big{(}\begin{array}{ccc}
J & J & 2n-1 \\
2 & -2 &
\end{array}\Big{)
\Big{(}\begin{array}{ccc}
J & J & 2n-1 \\
M & -M &
\end{array}\Big{)
 P_{2n-1}\left( \cos \theta \right).   \tag{11.a}
\end{equation

\begin{equation}
\zeta ^{JM}\left( \theta \right) =\frac{A_{JM}}{w_{2}}\sum\limits_{n=2}^{J
\left( 4n+1\right) \left[ \frac{\left( 2n-4\right) !}{\left( 2n+4\right) !
\right] ^{1/2}\Big{(}
\begin{array}{ccc}
J & J & 2n \\
2 & 2 & -
\end{array
\Big{)}
\Big{(}
\begin{array}{ccc}
J & J & 2n \\
M & -M &
\end{array
\Big{)} P_{2n}\left( \cos \theta \right).  \tag{11.b}
\end{equation}

\section{Polarization correlations}

The angular distributions (5) of the two-photon radiation contains information on the quantum numbers $J$ and $M$. Information related to the parity $P$ of TPSS is determined by linear polarization correlation measurements \cite{Yang}.

Photon polarization are measured by means of polarization analyzers whose efficiency depends on the value of the helicity. The efficiency of the polarization analyzer is determined by the efficiency matrix with
orientation vector $\boldsymbol{\eta }$ ($\eta =1$) having the form of a Stokes matrix, viz.,
\begin{equation}
\hat{\epsilon}=\frac{1}{2}(\hat{I}+\text{\textbf{$\eta $}}\cdot \mathbf{\hat
\sigma}}),  \tag{12}
\end{equation
which acts in the helicity Hilbert subspace reducing the mixed quantum state of polarization into a pure quantum state \cite{RB}.

The two-photon polarization is measured by two polarization analyzers $\hat{\epsilon}$ and $\hat{\epsilon}^{\prime}$ measuring the helicity states of the individual photons propagating along $\mathbf{n}$ and $\mathbf{n}^{\prime}$, respectively. The result of the measurement is defined as the probability of the reduction of the two-photon polarization matrix $\hat{\rho}$ into the pure quantum state of two uncorrelated photons described by $\hat{\epsilon}\otimes \hat{\epsilon}^{\prime }$, i.e.,

\begin{equation}
W(\mathbf{\eta ,\eta }^{\prime })=\textup{Tr}(\hat{\rho}\hat{\epsilon}\otimes \hat
\epsilon}^{\prime }).  \tag{13}
\end{equation}

In other words, $W(\mathbf{\eta ,\eta }^{\prime })$ is the probability of registering a photon in the definite helicity state $\left\vert \mathbf{n;}\lambda \right\rangle $ in coincidence with the other photon in the helicity state $\left\vert \mathbf{n}^{\prime }\mathbf{;}\lambda ^{\prime}\right\rangle $. Measurement of the probability $W(\mathbf{\eta,\eta }^{\prime })$ involves definite orientations and the kind of analyzers used. These measurements characterize the quantum state of the photon pair in an unambiguous way \cite{vem, ma}.

The efficiency matrix $\hat{\epsilon}$ of an analyzer measuring the linear polarization $\mathbf{\eta }$ having oriented angle $\psi $ relative to axes $x$ for the photon propagation direction $\mathbf{n}$ is given by \cite{Blum}
\begin{equation}
\hat{\epsilon}=\frac{1}{2}\left( \hat{I}-\cos(2\psi) \hat{\sigma}_{1}-\sin(2\psi) \hat{\sigma}_{2}\right).  \tag{14.a}
\end{equation
The efficiency matrix $\hat{\epsilon}^{\prime }$ $\left( \mathbf{\eta }^{\prime }\right)$ for the photon propagating along $\mathbf{-n}$ is similarly
\begin{equation}
\hat{\epsilon}^{\prime }=\frac{1}{2}\left( \hat{I}^{\prime }-\cos(2\psi^{\prime}) \hat{\sigma}_{1}^{\prime }+\sin(2\psi^{\prime }) \hat{\sigma}_{2}^{\prime }\right).  \tag{14.b}
\end{equation
On substituting (14) in (13), we find for the two-photon polarization density matrices $\hat{\rho}^{\pm }$ (7)
\begin{equation}
W^{\pm }(\psi ,\psi ^{\prime })=\textup{Tr}(\hat{\rho}^{\pm }\hat{\epsilon}\otimes
\hat{\epsilon}^{\prime })=\frac{1}{2}\left\{
\begin{array}{c}
\cos ^{2}\left( \psi -\psi ^{\prime }\right)  \\
\sin ^{2}\left( \psi -\psi ^{\prime }\right).
\end{array
\right . \tag{15}
\end{equation

Probability (15) indicates that in the state with positive parity, the photons are polarized predominately with parallel linear polarizations. On the other hand, with negative parity, the photons have linear polarization that are predominately perpendicular to each other. In addition, the photon polarization correlations do not have angular dependence, viz., photons are invariant relative to $Z$ to spatial rotations.

These type of correlations were mentioned by Yang \cite{Yang}. In addition, $W^{+}(\psi ,\psi ^{\prime })$ is the well-known $\left( \mathbf{ee}^{\prime }\right)^{2}$ polarization correlation for the Goeppert-Mayer's two-photon emission \cite{GM, Br-Tel}. Owing to momentum and parity conservation for the 2S-1S atomic transition, we have the emission of the two photons in the $\left\vert 00+1\right\rangle $\ Landau state.

For the polarization density matrices $\hat{\rho}_{e,o}^{JM}$ (8) we find, in an analogous way,
\begin{equation}
W_{e,o}(\psi ,\psi ^{\prime })=\textup{Tr}(\hat{\rho}_{e,o}\hat{\epsilon}\otimes \hat
\epsilon}^{\prime })=\frac{1}{4}\left\{
\begin{array}{c}
1+\zeta ^{JM}\cos[2\left(\psi +\psi ^{\prime}\right)]  \\
1-\zeta ^{JM}\cos[ 2\left( \psi +\psi ^{\prime }\right)],
\end{array
\right.   \tag{16}
\end{equation}
which clearly depends on $\theta$.

It should be remarked that in (16), the angles $\psi$ and $\psi^{\prime}$ in (14) are defined in terms of the $x$ axis. In addition, the direction of \textbf{n} is arbitrary relative to the quantization axis $Z$ and so one can take the direction of $x$ as being normal to the plane (\textbf{n},\textbf{Z)}. Consider the choice $\psi^{\prime} = 0$, that is, the analyzer $\hat{\epsilon}^{\prime}$ is orientated along the $x$-axis, then (16) becomes
\begin{equation}
W_{e,o}(\psi ,\psi ^{\prime })=\textup{Tr}(\hat{\rho}_{e,o}\hat{\epsilon}\otimes \hat
\epsilon}^{\prime })=\frac{1}{4}\left\{
\begin{array}{c}
1+\zeta ^{JM}\cos(2\psi)  \\
1-\zeta ^{JM}\cos(2\psi).
\end{array
\right.   \tag{17}
\end{equation}
If \textbf{n} and $Z$ are parallel, that is, $\theta=0$, then the photon correlation disappears since $\zeta^{JM}=0$.

Expression (17) for the correlation is further simplified for the case $\theta= \pi/2$ (see (10)), we then have
\begin{equation}
W_{e}(\psi ,\psi ^{\prime })=W_{o}(\psi ,\psi ^{\prime }) =\frac{1}{2}\left\{
\begin{array}{c}
\cos^2\psi  \hspace{.1in} (M=even) \\
\sin^2\psi  \hspace{.1in} (M=odd).
\end{array
\right.   \tag{18}
\end{equation}
Accordingly, if one of the linear analyzers is oriented normal to the (\textbf{n},$Z$) plane, then for M even the photons are polarized predominantly parallel to the plane while for M odd, the polarization is predominantly perpendicular to the plane.

In closing, we note that consideration of the correlation between circularly polarized photons does not provide any information about the parity of the two-photon state.

\section{Conclusion}

We have analyzed in detail the angular-polarization properties of the Landau two-photon spherical states in momentum space in the center of mass reference frame of the two-photon system. The angular distributions do not depend on the parity and, for fixed values of $J$ and $M$, are defined by two different functions of the polar angle between the relative momentum and the quantization axes. The two-photon polarization density matrices are derived for each values of $J$, $M$, and $P$. The linear polarization correlations of individual photons are analyzed in details. We find, besides the ordinary correlation laws given in terms of $sin$ and $cos$ of the angle between the orientation of the analyzers for $J\geq 2$, also the correlations in terms of $sin$ and $cos$ of the sum of the orientation angles of analyzers.\\


\section{Introduction}

Black holes are shown to carry soft hairs \cite{Hawking:2016msc,Compere:2016hzt,Afshar:2016wfy,Mao:2016pwq,Grumiller:2016kcp,Hawking:2016sgy,Chu:2018tzu,Choi:2018oel,Choi:2019fuq,Takeuchi:2021rrq}. The soft hairs in three dimensions are purely microscopic \cite{Afshar:2016uax,Sheikh-Jabbari:2016npa}. In four dimensions, the same feature has been verified by checking the Hawking radiation of soft-hairy black holes \cite{Chu:2018tzu,Javadinezhad:2018urv,Compere:2019rof,Lin:2020gva,Takeuchi:2021ibg}. More precisely, the Hawking fluxes from Schwarzschild black hole is unaffected by their linearized supertranslation hair \cite{Chu:2018tzu}. The spectrum of Hawking radiation of the Schwarzschild black hole is unchanged after including the dressing of asymptotic states with clouds of soft photons and soft gravitons \cite{Javadinezhad:2018urv}. Unruh effect is invariant under supertranslations \cite{Compere:2019rof}. The Hawking flux from Schwarzschild with non-linear supertranslation hair is still the same as the Schwarzschild's one \cite{Lin:2020gva} (see also \cite{Takeuchi:2021ibg}). However counter examples are presented in \cite{Chu:2018tzu,Chiang:2020lem,Wen:2021ahw}. In particular, the Hawking temperature of dynamical black holes is found to be modified by soft hair \cite{Chu:2018tzu,Chiang:2020lem}. Complementary to these investigations where only scalar particles are involved, we study in this paper the Hawking radiation of Dirac particles from soft-hairy black holes using the tunneling method.

The tunneling method proposed by Parikh and Wilczek \cite{Parikh:1999mf} is a very effective way of studying Hawking radiation. In general, the tunneling method involves calculating the imaginary part of the action for the classically forbidden path to evaluate the radiation spectrum. Applying the tunneling method, we examine the Hawking radiation of Dirac particles from Schwarzschild and Vaidya black holes with linear supertranslation hair. We find that the radiation spectrums and the Hawking temperatures in both cases are unchanged after soft hair implant on black holes. It may not be surprising that the linear supertranslation hair of Schwarzschild black hole does not affect the Hawking radiation of Dirac particles as the scalar case \cite{Chu:2018tzu}. For the Vaidya case, our result is different from the one of scalar particles in \cite{Chu:2018tzu}. When applying the  tunneling method for Dirac particles, one needs to properly take care of the angular-dependent part of the action \cite{Kerner:2007rr,Li:2008zra,DiCriscienzo:2008dm,Chen:2008vi,Kerner:2008qv,Chen:2008ge,Jiang:2008gq,Chen:2014xgj}. The simplest assumption that solving only the radial equation at a fixed angle may not be a solution of the Dirac equation. We find that the linear supertranslation hair dependence in the Hawking spectrum of the Vaidya black hole is precisely canceled by taking account the contribution from the angular-dependent part. Our results furnish an important confirmation of the microscopic nature of soft hairs in four dimensions.

\section{Schwarzschild black hole with supertranslation hair}

The line element of the supertranslated Schwarzschild black hole in the advanced Bondi coordinates reads \cite{Hawking:2016sgy}
\begin{multline}\label{sch-super}
\td s^2=-\left(V-\frac{M D^2f}{r^2}\right) \td v^2 + 2 \td v \td r - D_A \left(2 V f + D^2 f \right) \td v \td \Theta^A \\
+\left(r^2 {\bar\gamma}_{AB} + 2 r D_A D_B f - r {\bar\gamma}_{AB} D^2 f \right) \td \Theta^A \td \Theta^B \ , \quad V=1-\frac{2M}{r} \ ,
\end{multline}
where ${\bar\gamma}_{AB}$ is the metric of the unit two sphere on the boundary and $D_A$ is the covariant derivative with respect to ${\bar\gamma}_{AB}$, $\Theta^A=(z, \bar z)$ are the angular coordinates, $f$ is an arbitrary function of the angular variables that characterizes the supertanslation. Up to order $f$ term, $g_{vv}=-\frac{1}{r^2}(r-r_h)(r+r_f)$ where $r_f=\frac12 D^2 f$ and  $r_h=2M + \frac12 D^2 f$ is the location of event horizon. Since the supertranslation is infinitesimal, we drop the terms with order higher than $\cO(f)$. Physically, the solution can be considered as throwing in an asymmetric null shock into the Schwarzschild black hole.

The crucial step of computing Hawking radiation by tunneling method is the application of WKB approximation to the covariant Dirac equation \cite{Kerner:2007rr}. We consider massive spinor field $\Psi$ with mass $\mu$. The covariant Dirac equation reads
\begin{equation}\label{Dirac}
i \hbar \gamma^a e_a^\mu \nabla_\mu \Psi - \mu \Psi = 0 \ ,
\end{equation}
where $\nabla_\mu$ is the spinor-covariant derivative defined by $\nabla_\mu=\partial_\mu + \frac14 {\omega^{ab}}_\mu \gamma_{[a}\gamma_{b]}$ and ${\omega^{ab}}_\mu$ is the spin connection. The cotetrad fields $e^a_\mu$ according to the metric \eqref{sch-super} are chosen as
\begin{align}
&e^0_\mu=\left[ -\frac12 g_{vv}, -1 , - g_{vA}\right] \ ,\\
&e^1_\mu=\left[1,0,0,0\right] \ ,\\
&e^2_\mu=\left[0,0,f_1(x^\mu),f_2(x^\mu)\right] \ ,\\
&e^3_\mu=\left[0,0,0,f_3(x^\mu)\right] \ ,
\end{align}
where $f_1=\sqrt{g_{zz}},\,f_2=\frac{g_{z\bar z}}{\sqrt{g_{zz}}},\,f_3=\sqrt{g_{z\bar z} - \frac{g_{z\bar z}^2}{g_{zz}}}$ and the fundamental matrix is of the form
\begin{equation}
\eta_{ab}=\eta^{ab}=\begin{pmatrix}
0 & -1 & 0 & 0  \\ -1 & 0 & 0 & 0  \\ 0 & 0 & 1 & 0  \\ 0 & 0 & 0 & 1
\end{pmatrix}\ .
\end{equation}
We will pick for the $\gamma$ matrices
\begin{align}
&\gamma^0=\frac{1}{\sqrt{2}}\begin{pmatrix}
i & 0 \\ 0 & - i
\end{pmatrix}  + \frac{1}{\sqrt{2}} \begin{pmatrix}
0 & \sigma^3 \\ \sigma^3 & 0
\end{pmatrix}\ ,\quad \gamma^2=\begin{pmatrix}
0 & \sigma^1 \\ \sigma^1 & 0
\end{pmatrix}\ ,\\
&\gamma^1=\frac{1}{\sqrt{2}}\begin{pmatrix}
i & 0 \\ 0 & - i
\end{pmatrix}  - \frac{1}{\sqrt{2}} \begin{pmatrix}
0 & \sigma^3 \\ \sigma^3 & 0
\end{pmatrix} \ ,
\quad
\gamma^3=\begin{pmatrix}
0 & \sigma^2 \\ \sigma^2 & 0
\end{pmatrix}\ ,
\end{align}
where the matrices $\sigma^k(k = 1,2,3)$ are the Pauli matrices. The $\gamma^a$ matrices satisfy the anti-commutation relations $\{\gamma^a,\gamma^b\}=2\eta^{ab}\times \mathbb{1}$.
The spin-up and spin-down spinor fields are assumed to have the form \cite{Kerner:2007rr}
\begin{equation}
\Psi_{\uparrow}=\begin{pmatrix}
A(x^\mu)  \\  0  \\ B(x^\mu)   \\ 0
\end{pmatrix}\text{exp}\left[\frac{i}{\hbar}I_{\uparrow}(x^\mu)\right] \ ,\quad
\Psi_{\downarrow}=\begin{pmatrix}
0 \\ C(x^\mu)  \\  0  \\ D(x^\mu)
\end{pmatrix}\text{exp}\left[\frac{i}{\hbar}I_{\downarrow}(x^\mu)\right] \ .
\end{equation}
We will only demonstrate the spin-up case explicitly while the spin-down case is fully in analogue.

Inserting the spin-up ansatz into the covariant Dirac equation \eqref{Dirac} and applying the WKB approximation, one can arrive at the following four equations at the leading order in the expansion of $\hbar$:
\begin{align}
&\sqrt{2} B \left[\partial_r I (g_{vv} - 2) - 2 \partial_v I \right] + A \left[4 \mu + \sqrt2 i (2 + g_{vv}) \partial_r I - 2\sqrt2 i \partial_v I \right] = 0 \ ,\\
&B \left[ (f_2 + i f_3 ) (\partial_z I - g_{vz} \partial_r I) - f_1 (\partial_{\bar z} I - g_{v\bar z} \partial_r I ) \right]= 0 \ ,\label{2}\\
&\sqrt{2} A \left[\partial_r I (g_{vv} - 2) - 2 \partial_v I \right] + B \left[4 \mu - \sqrt2 i (2 + g_{vv}) \partial_r I + 2\sqrt2 i \partial_v I \right] = 0 \ ,\\
&A \left[ (f_2 + i f_3 ) (\partial_z I - g_{vz} \partial_r I) - f_1 (\partial_{\bar z} I - g_{v\bar z} \partial_r I ) \right]= 0 \ .\label{4}
\end{align}
For the first and third equations, the non-trivial solution for A and B exists if and only if the determinant of the coefficient matrix vanishes. Then we obtain
\begin{equation}\label{r}
g_{vv} (\partial_r I)^2 - 2 \partial_v I \partial_r I + \mu^2= 0 \ .
\end{equation}
Equations \eqref{2} and \eqref{4} yield the same equation regardless of $A$ or $B$, i.e.,
\begin{equation}\label{angle}
(f_2 + i f_3 ) (\partial_z I - g_{vz} \partial_r I) - f_1 (\partial_{\bar z} I - g_{v\bar z} \partial_r I )=0 \ .
\end{equation}
Because $\frac{\partial}{\partial_v}$ is a Killing vector in the supertranslated spacetime, we can separate the variables for $I(v,r,z,\bar z)$ as follows:
\begin{equation}
I =- \omega v + R(r,z,\bar z) + K \ ,
\end{equation}
where $\omega$ is the Dirac particle's energy and $K$ is a complex constant. Then \eqref{r} can be solved out as
\begin{align}
&R_+=\int \td r \, \frac{\omega + \sqrt{\omega^2 - \mu^2 g_{vv}}}{-g_{vv}} + \Phi_+(z,\bar z) \ ,\label{+} \\
&R_-=\int \td r \, \frac{\omega - \sqrt{\omega^2 - \mu^2 g_{vv}}}{-g_{vv}} + \Phi_-(z,\bar z) \ ,\label{-}
\end{align}
where $R_+$ denotes the outgoing solution and $R_-$ denotes the incoming solution. The angular-dependent part of the action for outgoing solution is usually the same as the incoming solution for rotating black holes, see, e.g., \cite{Li:2008zra,DiCriscienzo:2008dm,Chen:2008vi,Kerner:2008qv,Chen:2008ge,Jiang:2008gq,Chen:2014xgj}. Though it is not necessary to be a real function, it will be canceled out when computing the probability of tunneling. However in the present case, equation \eqref{angle} involves the $r$-derivative of the action. A priori, it is necessary to consider them as different functions. Nevertheless we will argue that $\Phi_+(z,\bar z)$ and $\Phi_-(z,\bar z)$ can have real solutions. We will consider the case that the two expressions $\partial_z I - g_{vz} \partial_r I$ and $\partial_{\bar z} I - g_{v\bar z} \partial_r I$ in \eqref{angle} are complex conjugate to each other. The complex conjugate of the first expression for the incoming solution is
\begin{equation}
\partial_{\bar z} \left[\int \td r \, \frac{\omega - \sqrt{\omega^2 - \mu^2 g_{vv}}}{-g_{vv}} + \xbar\Phi_-(z,\bar z)\right] - g_{v\bar z} \frac{\omega - \sqrt{\omega^2 - \mu^2 g_{vv}}}{-g_{vv}} \ ,
\end{equation}
where we have used the fact that $\int \td r \, \frac{\omega - \sqrt{\omega^2 - \mu^2 g_{vv}}}{-g_{vv}}$ is real and $\xbar\Phi_-(z,\bar z)$ is the complex conjugate function of $\Phi_-(z,\bar z)$. Comparing this equation to the second expression
\begin{equation}
\partial_{\bar z} \left[\int \td r \, \frac{\omega - \sqrt{\omega^2 - \mu^2 g_{vv}}}{-g_{vv}} + \Phi_-(z,\bar z)\right] - g_{v\bar z} \frac{\omega - \sqrt{\omega^2 - \mu^2 g_{vv}}}{-g_{vv}} \ ,
\end{equation}
one obtains that $\xbar\Phi_-(z,\bar z)=\Phi_-(z,\bar z)$. Hence $\Phi_-(z,\bar z)$ is real.

For the outgoing solution, the expression $\int \td r \, \frac{\omega + \sqrt{\omega^2 - \mu^2 g_{vv}}}{-g_{vv}}$ is not real. But the imaginary part of this expression is $4\pi \omega M$ \cite{Chu:2018tzu}, i.e., a constant that will be killed by a $z$ or $\bar z$ derivative. Hence, one can show that $\Phi_+(z,\bar z)$ is real following the same argument as the incoming case.

The probabilities of crossing the outer horizon each way are respectively given by
\begin{align}
&P_{\text{out}}=\text{exp}\left[-2 \text{Im}I_+\right]= \text{exp}\left[- 8\pi \omega M  - 2 \text{Im}K \right]\ ,\\
&P_{\text{in}}=\text{exp}\left[-2 \text{Im}I_-\right]=\text{exp}\left[-2 \text{Im}K\right] \ ,
\end{align}
where we set $\hbar$ to unity. Any incoming classical particles should have a $100\%$ chance of entering the black hole horizon. So we set $K=0$. This implies that the probability of a particle tunneling from inside to outside of the horizon is given by
\begin{equation}
\Gamma= \text{exp}\left[- 8\pi \omega M\right] \ ,
\end{equation}
which recovers the result of the scalar particles tunneling from the supertranslated Schwarzschild black hole \cite{Chu:2018tzu}. The Hawking radiation is independent of the supertranslation hair of the black hole. Identifying $\Gamma$ with the Boltzmann factor, exp$(-\frac{\omega}{T} )$, we can read off the Hawking
temperature as $\frac{1}{8\pi M}$.

\section{Vaidya black hole with supertranslation hair}

The line element of the supertranslated Vaidya black hole in the advanced Bondi coordinates is \cite{Chu:2018tzu}
\begin{multline}
\td s^2=-\left(V - \frac{2fM'(v)}{r} - \frac{M D^2f}{r^2}\right) \td v^2 + 2 \td v \td r -  D_A \left(2 V f + D^2 f \right) \td v \td \Theta^A \\
+\left(r^2 \gamma_{AB} + 2 r D_A D_B f - r \gamma_{AB} D^2 f \right) \td \Theta^A \td \Theta^B \ ,
\end{multline}
where $V=1-\frac{2M(v)}{r}$ and $M'(v)=\frac{\td M(v)}{\td v}$. Up to order $f$ term, $g_{vv}=-\frac{1}{r^2}(r-r_h)(r+r_f)$ where $r_f=\frac12 D^2 f$ and $r_h=2M + 2f M' + \frac12 D^2 f$ is the location of trapping horizon. Physically, the solution can be considered as collapse of a certain energy momentum tensor with soft charges.

The line element of the supertranslated Vaidya black hole has the same form as the supertranslated Schwarzschild one. Hence the covariant Dirac equation with the same choice of $\gamma$ matrices and the same ansatz for spin-up spinor fields, after applying the WKB approximation, will lead to the same formulas as the supertranslated Schwarzschild case, i.e., \eqref{r} and \eqref{angle}. Following closely \cite{Chu:2018tzu}, we define the invariant energy of particle $\omega$ in favor of the Kodama-like vector which yields $\partial_v I=-\omega$. Substituting this back to \eqref{r}, one obtains
\begin{align}
&I_+=-\omega v + \int \td r \, \frac{\omega + \sqrt{\omega^2 - \mu^2 g_{vv}}}{-g_{vv}} + \Phi_+(z,\bar z) + K \ ,\label{plus} \\
&I_-=-\omega v + \int \td r \, \frac{\omega - \sqrt{\omega^2 - \mu^2 g_{vv}}}{-g_{vv}} + \Phi_-(z,\bar z) + K \ .\label{minus}
\end{align}
Since function $\int \td r \, \frac{\omega - \sqrt{\omega^2 - \mu^2 g_{vv}}}{-g_{vv}}$ is real \cite{Chu:2018tzu}, $\Phi_-(z,\bar z)$ is real applying the same argument as previous section. We set $K=0$ due to the fact that the probability of any incoming classical particles crossing the horizon is unity. A direct computation shows that \cite{Chu:2018tzu}
\begin{equation}\label{Im}
\text{Im}\int \td r \, \frac{\omega + \sqrt{\omega^2 - \mu^2 g_{vv}}}{-g_{vv}}=4\pi \omega M + 4\pi \omega f M' \ ,
\end{equation}
where higher orders in $f$ has been neglected. Applying the same argument as previous section, \eqref{Im} will lead to
\begin{equation}
\text{Im}\Phi_+(z,\bar z)=-4\pi \omega  f M' \ .
\end{equation}
Hence
\begin{equation}
\text{Im}I_+=4\pi \omega M \ .
\end{equation}
In the end, the probability of a particle tunneling from inside to outside of the horizon is given by
\begin{equation}
\Gamma= \text{exp}\left[- 8\pi \omega M\right] \ ,
\end{equation}
which is different from the result of scalar particles tunneling from the supertranslated Vaidya black hole \cite{Chu:2018tzu}. Note that the Hawking radiation is again independent of the supertranslation hair when taking into account the angular contribution. By identifying $\Gamma$ with the Boltzmann factor, exp$(-\frac{\omega}{T} )$, the Hawking
temperature is $\frac{1}{8\pi M}$.

\section{Conclusion}

In this work, we examine the semi-classical property of soft hair implants on black holes. Applying the tunneling method, we compute the Hawking radiation of Dirac particles from Schwarzschild and Vaidya black holes with linear supertranslation hair. Our calculation provides the first example that Hawking radiation of dynamical black hole is not modified by supertranslation hair. This is another evidence to support the microscopic nature of soft hairs.

\section*{Acknowledgments}

This work is supported in part by the National Natural Science Foundation of China under Grant No. 11905156, No. 11975164 and No. 11935009.

\providecommand{\href}[2]{#2}\begingroup\raggedright