Emergent $\alpha$-like fermionic vacuum structure and entanglement in the hyperbolic de Sitter spacetime

TL;DR

This paper uncovers a new vacuum structure in hyperbolic de Sitter spacetime for Dirac fields, revealing a one-parameter family of vacua affecting quantum correlations and entanglement between causally disconnected regions.

Contribution

It introduces a novel orthonormalisation approach for global modes, leading to a continuous parameter-dependent family of de Sitter vacua with implications for quantum entanglement.

Findings

Emergence of a one-parameter family of vacua characterized by $ heta_{RL}$.

Dependence of entanglement and Rénnyi entropies on the vacuum parameter.

Identification of parameter values corresponding to thermal particle spectra.

Abstract

We report a non-trivial feature of the vacuum structure of free massive or massless Dirac fields in the hyperbolic de Sitter spacetime. Here we have two causally disconnected regions, say $R$ and $L$ separated by another region, $C$. We are interested in the field theory in $R\cup L$ to understand the long range quantum correlations between $R$ and $L$. There are local modes of the Dirac field having supports individually either in $R$ or $L$, as well as global modes found via analytically continuing the $R$ modes to $L$ and vice versa. However, we show that unlike the case of a scalar field, the analytic continuation does not preserve the orthogonality of the resulting global modes. Accordingly, we need to orthonormalise them following the Gram-Schmidt prescription, prior to the field quantisation in order to preserve the canonical anti-commutation relations. We observe that this prescription naturally incorporates a spacetime independent continuous parameter, $\theta_{\rm RL}$, into the picture. Thus interestingly, we obtain a naturally emerging one-parameter family of $\alpha$-like de Sitter vacua. The values of $\theta_{\rm RL}$ yielding the usual thermal spectra of massless created particles are pointed out. Next, using these vacua, we investigate both entanglement and R\'enyi entropies of either of the regions and demonstrate their dependence on $\theta_{\rm RL}$.

Equations

$$t_R=i\left(\tau -\frac{\pi }{2}\right),r_R=i\rho ,\left(t_R\ge 0,r_R\ge 0\right);t_L=i\left(-\tau -\frac{\pi }{2}\right),r_L=i\rho ,\left(t_L\ge 0,r_L\ge 0\right)$$

$$t_C=\tau ,r_C=i\left(\rho -\frac{\pi }{2}\right)\left(-\frac{\pi }{2}\le t_C\le \frac{\pi }{2},0\le r_C<\infty \right).$$

$$d{s}_{R,L}^2=H^{-2}\left(-d{t}_{R,L}^2+{\sinh }^2{t}_{R,L}(d{r}_{R,L}^2+{\sinh }^2{r}_{R,L}d{\Omega }^2)\right)$$

$$d{s}_C^2=H^{-2}\left(d{t}_C^2+{\cos }^2{t}_C(-d{r}_C^2+{\cosh }^2{r}_{C}d{\Omega }^2)\right)$$

$$\displaystyle\Psi^{+R}_{(-)}=\left(\begin{array}[]{l}u_{p}(z_{R})\\ v_{p}(z_{R})\end{array}\right)\chi^{(-)}_{pjm}\qquad\qquad\Psi^{-R}_{(-)}=\left(\begin{array}[]{l}-v_{p}^{\star}(z_{R})\\ \ \ u_{p}^{\star}(z_{R})\end{array}\right)\chi^{(-)}_{pjm}$$

$$\displaystyle\Psi^{-R}_{(+)}=\left(\begin{array}[]{l}\ \ u_{p}^{\star}(z_{R})\\ -v_{p}^{\star}(z_{R})\end{array}\right)\chi^{(+)}_{pjm}\qquad\qquad\Psi^{+R}_{(+)}=\left(\begin{array}[]{l}v_{p}(z_{R})\\ u_{p}(z_{R})\end{array}\right)\chi^{(+)}_{pjm}.$$

$$\widetilde{/\nabla }{\chi }_{pjm}^{(\pm )}=\pm ip{\chi }_{pjm}^{(\pm )},$$

$$u_p(z_R)$$

$$v_p(z_R)$$

$$\displaystyle\Psi^{+L}_{(-)}(z_{L})=\left(\begin{array}[]{r}v_{p}(z_{L})\\ -u_{p}(z_{L})\end{array}\right)\chi^{(-)}_{pjm}\qquad\Psi^{-L}_{(-)}(z_{L})=\left(\begin{array}[]{r}u^{\star}_{p}(z_{L})\\ \ \ v^{\star}_{p}(z_{L})\end{array}\right)\chi^{(-)}_{pjm}$$

$$\displaystyle\Psi^{-L}_{(+)}(z_{L})=\left(\begin{array}[]{r}\ \ v^{\star}_{p}(z_{L})\\ u^{\star}_{p}(z_{L})\end{array}\right)\chi^{(+)}_{pjm}\qquad\Psi^{+L}_{(+)}(z_{L})=\left(\begin{array}[]{r}-u_{p}(z_{L})\\ v_{p}(z_{L})\end{array}\right)\chi^{(+)}_{pjm}$$

$$\left({\Psi }_{(a)},{\Psi }_{(b)}\right):=\int {\sinh }^{3}t\sqrt{h}drd\theta d\varphi {\Psi }_{(a)}^{†}{\Psi }_{(b)}$$

$$\left({\Psi }_{(s)}^{+R}(p,j,m),{\Psi }_{(s^{\prime })}^{+R}(p^{\prime },j^{\prime },m^{\prime })\right)=\delta (p-p^{\prime }){\delta }_{j{j}^{\prime }}{\delta }_{m{m}^{\prime }}{\delta }_{s{s}^{\prime }}\left({\Psi }_{(s)}^{+R}(p,j,m),{\Psi }_{(s^{\prime })}^{-R}(p^{\prime },j^{\prime },m^{\prime })\right)=0etc.$$

$$\left({\Psi }_{(s)}^{+L}(p,j,m),{\Psi }_{(s^{\prime })}^{+L}(p^{\prime },j^{\prime },m^{\prime })\right)=\delta (p-p^{\prime }){\delta }_{j{j}^{\prime }}{\delta }_{m{m}^{\prime }}{\delta }_{s{s}^{\prime }}\left({\Psi }_{(s)}^{+L}(p,j,m),{\Psi }_{(s^{\prime })}^{-L}(p^{\prime },j^{\prime },m^{\prime })\right)=0etc.$$

$$\displaystyle\Psi^{+R}_{(-)}(z_{R})=\left(\begin{array}[]{r}\ \ u_{p}(z_{R})\\ v_{p}(z_{R})\end{array}\right)\chi^{(-)}_{pjm},\qquad\Psi^{+R}_{(-)}(z_{L})=\left(\begin{array}[]{r}\lambda_{1}\,v_{p}(z_{L})+\lambda_{2}\,u_{p}^{\star}(z_{L})\\ -\lambda_{1}\,u_{p}(z_{L})+\lambda_{2}\,v_{p}^{\star}(z_{L})\end{array}\right)\chi^{(-)}_{pjm}$$

$$\displaystyle\Psi^{+R}_{(+)}(z_{R})=\left(\begin{array}[]{r}\ \ v_{p}(z_{R})\\ u_{p}(z_{R})\end{array}\right)\chi^{(+)}_{pjm},\qquad\Psi^{+R}_{(+)}(z_{L})=\left(\begin{array}[]{r}-\lambda_{1}\,u_{p}(z_{L})+\lambda_{2}\,v_{p}^{\star}(z_{L})\\ \lambda_{1}\,v_{p}(z_{L})+\lambda_{2}\,u_{p}^{\star}(z_{L})\end{array}\right)\chi^{(+)}_{pjm}$$

$$\displaystyle\Psi^{-R}_{(-)}(z_{R})=\left(\begin{array}[]{r}-v^{\star}_{p}(z_{R})\\ u^{\star}_{p}(z_{R})\end{array}\right)\chi^{(-)}_{pjm},\qquad\Psi^{-R}_{(-)}(z_{L})=\left(\begin{array}[]{r}\ \ \lambda_{1}u^{\star}_{p}(z_{L})-\lambda_{2}^{\star}v_{p}(z_{L})\\ \lambda_{1}v^{\star}_{p}(z_{L})+\lambda_{2}^{\star}u_{p}(z_{L})\end{array}\right)\chi^{(-)}_{pjm}$$

$$\displaystyle\Psi^{-R}_{(+)}(z_{R})=\left(\begin{array}[]{r}u^{\star}_{p}(z_{R})\\ -v^{\star}_{p}(z_{R})\end{array}\right)\chi^{(+)}_{pjm},\qquad\Psi^{-R}_{(+)}(z_{L})=\left(\begin{array}[]{r}\ \ \lambda_{1}v^{\star}_{p}(z_{L})+\lambda_{2}^{\star}u_{p}(z_{L})\\ \lambda_{1}u^{\star}_{p}(z_{L})-\lambda_{2}^{\star}v_{p}(z_{L})\end{array}\right)\chi^{(+)}_{pjm}$$

$${\lambda }_1=\frac{\sinh \frac{m\pi }{H}}{\cosh \pi p},{\lambda }_2=\frac{e^{-\pi p}{\left(\Gamma \left(\frac{1}{2}-ip\right)\right)}^2}{\Gamma \left(\frac{1}{2}-ip-\frac{im}{H}\right)\Gamma \left(\frac{1}{2}-ip+\frac{im}{H}\right)}.$$

$$\displaystyle\Psi^{+L}_{(-)}(z_{L})=\left(\begin{array}[]{r}v_{p}(z_{L})\\ -u_{p}(z_{L})\end{array}\right)\chi^{(-)}_{pjm},\qquad\Psi^{+L}_{(-)}(z_{R})=\left(\begin{array}[]{r}-\lambda_{1}\,u_{p}(z_{R})+\lambda_{2}\,v_{p}^{\star}(z_{R})\\ -\lambda_{1}\,v_{p}(z_{R})-\lambda_{2}\,u_{p}^{\star}(z_{R})\end{array}\right)\chi^{(-)}_{pjm}$$

$$\displaystyle\Psi^{+L}_{(+)}(z_{L})=\left(\begin{array}[]{r}-u_{p}(z_{L})\\ v_{p}(z_{L})\end{array}\right)\chi^{(+)}_{pjm},\qquad\Psi^{+L}_{(+)}(z_{R})=\left(\begin{array}[]{r}-\lambda_{1}\,v_{p}(z_{R})-\lambda_{2}\,u_{p}^{\star}(z_{R})\\ -\lambda_{1}\,u_{p}(z_{R})+\lambda_{2}\,v_{p}^{\star}(z_{R})\end{array}\right)\chi^{(+)}_{pjm}$$

$$\displaystyle\Psi^{-L}_{(-)}(z_{L})=\left(\begin{array}[]{r}u^{\star}_{p}(z_{L})\\ \ \ v^{\star}_{p}(z_{L})\end{array}\right)\chi^{(-)}_{pjm},\qquad\Psi^{-L}_{(-)}(z_{R})=\left(\begin{array}[]{r}\lambda_{1}v^{\star}_{p}(z_{R})+\lambda_{2}^{\star}u_{p}(z_{R})\\ -\lambda_{1}u^{\star}_{p}(z_{R})+\lambda_{2}^{\star}v_{p}(z_{R})\end{array}\right)\chi^{(-)}_{pjm}$$

$$\displaystyle\Psi^{-L}_{(+)}(z_{L})=\left(\begin{array}[]{r}\ \ v^{\star}_{p}(z_{L})\\ u^{\star}_{p}(z_{L})\end{array}\right)\chi^{(+)}_{pjm},\qquad\Psi^{-L}_{(+)}(z_{R})=\left(\begin{array}[]{r}\ -\lambda_{1}u^{\star}_{p}(z_{R})+\lambda_{2}^{\star}v_{p}(z_{R})\\ \lambda_{1}v^{\star}_{p}(z_{R})+\lambda_{2}^{\star}u_{p}(z_{R})\end{array}\right)\chi^{(+)}_{pjm}$$

$${\left({\Psi }_{(-)}^{+L},{\Psi }_{(-)}^{-R}\right)}_{\mathrm{G}}={\left({\Psi }_{(-)}^{+L},{\Psi }_{(-)}^{-R}\right)}_{z=z_L=1}+{\left({\Psi }_{(-)}^{+L},{\Psi }_{(-)}^{-R}\right)}_{z=z_R=1}$$

$${\left({\Psi }_{(-)}^{+L},{\Psi }_{(-)}^{-R}\right)}_{\mathrm{G}}=-2{\lambda }_2^{\star }$$

$${\left({\Psi }_{(-)}^{+L},{\Psi }_{(-)}^{-R}\right)}_{\mathrm{G}}={\left({\Psi }_{(+)}^{+L},{\Psi }_{(+)}^{-R}\right)}_{\mathrm{G}}=-{\left({\Psi }_{(-)}^{+R},{\Psi }_{(-)}^{-L}\right)}_{\mathrm{G}}=-{\left({\Psi }_{(+)}^{+R},{\Psi }_{(+)}^{-L}\right)}_{\mathrm{G}}=-2{\lambda }_2^{\star }\ne 0$$

$${\tilde{\Psi }}_{{\theta }_{\mathrm{RL}}}:={\Psi }_{(-)}^{-R}+\frac{2{\lambda }_2\Delta {\theta }_1}{N_b^2}{\Psi }_{(-)}^{+L}{\Psi }_{{\theta }_{\mathrm{RL}}}^{\prime }:={\Psi }_{(-)}^{+L}+\frac{2{\lambda }_2^{\star }\Delta {\theta }_2}{N_b^2}{\Psi }_{(-)}^{-R}$$

$$\Delta {\theta }_1=\frac{{\cos }^2{\theta }_{\text{RL}}}{1-\frac{2|{\lambda }_2|}{N_b^2}{\sin }^2{\theta }_{\text{RL}}}\Delta {\theta }_2=\frac{{\sin }^2{\theta }_{\text{RL}}}{1+\frac{2|{\lambda }_2|}{N_b^2}{\cos }^2{\theta }_{\text{RL}}}\left(0\le {\theta }_{\mathrm{RL}}\le \frac{\pi }{2}\right)$$

$$N_b^2=\left(1+{\lambda }_1^2+|{\lambda }_2{|}^2\right)$$

$${\Psi }_1=\frac{1}{N_1{N}_b}\left({\Psi }_{(-)}^{+R}-\frac{2{\lambda }_2\Delta {\theta }_1}{N_b^2}{\Psi }_{(-)}^{-L}\right){\Psi }_2=\frac{1}{N_2{N}_b}\left({\Psi }_{(+)}^{+R}-\frac{2{\lambda }_2\Delta {\theta }_2}{N_b^2}{\Psi }_{(+)}^{-L}\right)$$

$${\Psi }_3=\frac{1}{N_1{N}_b}\left({\Psi }_{(-)}^{+L}+\frac{2{\lambda }_2\Delta {\theta }_1}{N_b^2}{\Psi }_{(-)}^{-R}\right){\Psi }_4=\frac{1}{N_2{N}_b}\left({\Psi }_{(+)}^{+L}+\frac{2{\lambda }_2\Delta {\theta }_2}{N_b^2}{\Psi }_{(+)}^{-R}\right)$$

$${\Psi }_5=\frac{1}{N_1{N}_b}\left({\Psi }_{(+)}^{-R}+\frac{2{\lambda }_2^{\star }\Delta {\theta }_1}{N_b^2}{\Psi }_{(+)}^{+L}\right){\Psi }_6=\frac{1}{N_2{N}_b}\left({\Psi }_{(-)}^{-R}+\frac{2{\lambda }_2^{\star }\Delta {\theta }_2}{N_b^2}{\Psi }_{(-)}^{+L}\right)$$

$${\Psi }_7=\frac{1}{N_1{N}_b}\left({\Psi }_{(+)}^{-L}-\frac{2{\lambda }_2^{\star }\Delta {\theta }_1}{N_b^2}{\Psi }_{(+)}^{+R}\right){\Psi }_8=\frac{1}{N_2{N}_b}\left({\Psi }_{(-)}^{-L}-\frac{2{\lambda }_2^{\star }\Delta {\theta }_2}{N_b^2}{\Psi }_{(-)}^{+R}\right)$$

$$N_1^2=1+\frac{4|{\lambda }_2{|}^2\Delta {\theta }_1^2}{N_b^4}-\frac{8|{\lambda }_2{|}^2\Delta {\theta }_1}{N_b^4}{N}_2^2=1+\frac{4|{\lambda }_2{|}^2\Delta {\theta }_2^2}{N_b^4}-\frac{8|{\lambda }_2{|}^2\Delta {\theta }_2}{N_b^4}$$

$$\Psi =\int dp\sum _{jms}\left(c_{(s)pjm}^R{\Psi }_{(s)pjm}^{+R}+d_{(s)pjm}^{R†}{\Psi }_{(s)pjm}^{-R}+c_{(s)pjm}^L{\Psi }_{(s)pjm}^{+L}+d_{(s)pjm}^{L†}{\Psi }_{(s)pjm}^{-L}\right),$$

$${\left[c_{(s)pjm}^R,c_{(s^{\prime })p^{\prime }{j}^{\prime }{m}^{\prime }}^{R†}\right]}_{+}={\left[d_{(s)pjm}^R,d_{(s^{\prime })p^{\prime }{j}^{\prime }{m}^{\prime }}^{R†}\right]}_{+}=\delta \left(p-p^{\prime }\right){\delta }_{s{s}^{\prime }}{\delta }_{j{j}^{\prime }}{\delta }_{m{m}^{\prime }}$$

$${\left[c_{(s)pjm}^L,c_{(s^{\prime })p^{\prime }{j}^{\prime }{m}^{\prime }}^{L†}\right]}_{+}={\left[d_{(s)pjm}^L,d_{(s^{\prime })p^{\prime }{j}^{\prime }{m}^{\prime }}^{L†}\right]}_{+}=\delta \left(p-p^{\prime }\right){\delta }_{s{s}^{\prime }}{\delta }_{j{j}^{\prime }}{\delta }_{m{m}^{\prime }}$$

$$c_{(s)pjm}^R|0{⟩}_R=d_{(s)pjm}^R|0{⟩}_R=0,c_{(s)pjm}^L|0{⟩}_L=d_{(s)pjm}^L|0{⟩}_L=0.$$

Full text

Keywords : Hyperbolic de Sitter, fermionic vacua, quantum entanglement

1 Introduction

The de Sitter (dS) spacetime is a maximally symmetric manifold endowed with a positive cosmological constant. It is physically interesting in many ways. First, owing to the recent observed phase of the accelerated cosmic expansion, there is a strong possibility that our current universe is dominated by a small but positive cosmological constant at large scales. Second, the high degree of homogeneity and isotropy of the current universe at large scales indicate that our early universe went through an inflationary phase described by a quasi de Sitter spacetime [1]. Being accelerated with expansion and endowed with a cosmological event horizon, the dS offers interesting thermal and other field theoretic properties, we refer our reader to e.g., [2, 3, 4, 5, 6] and also references therein.

It is interesting to investigate the long range quantum correlations between two observers in the dS space, not only causally separated but so by a large distance, say of the order of the superhorizon size. This issue was first addressed in [7] for a scalar field theory using the coordinatisation of [8, 9], known as the hyperbolic or open chart describing two casually disconnected expanding regions in dS, as denoted by regions $R$ and $L$ in 1. Since $R$ and $L$ are separated by an entire causally disconnected region $C$, the framework described by 1 offers a very natural stage to investigate such long range non-local quantum correlations. Being motivated by this, we wish to investigate in the following the entanglement properties of the Dirac fermionic vacua in the hyperbolic dS.

Let us first briefly review the case of a real scalar field [7]. One first defines orthonormal local basis mode functions having supports either in $R$ or in $L$, with definite positive or negative frequency behaviour in the asymptotic past. One makes a field expansion using them and defines the local vacuum as a direct product between the vacua of $R$ and $L$. However, if there is any correlation between these local states, clearly there must exist some mode functions having supports in *both * $R$ and $L$. Such global modes are obtained by analytically continuing the local modes from one region into the another along a complex path going through $C$ [8, 9]. One then makes a field quantisation using the global modes and defines a suitable global vacuum. The field quantisations in terms of the local and global modes give a Bogoliubov relation, yielding in turn an expansion of the global vacuum in terms of the local states. It follows that the states belonging to $R$ and $L$ are entangled. The entanglement entropy density is computed using the reduced density operator, found by tracing out the states belonging to either of the regions. Being originated from the long range correlations, the entanglement entropy thus found will not be proportional to any area. This procedure will be more explicit in the due course of the discussion.

A lot of effort has been given to explore quantum entanglement in dS so far, in various coordinatisations, e.g. [10]-[28]. They not only involve the computations of the entanglement or the Rényi entropy (e.g. [29] and references therein), but also studies of other measures like Bell’s inequality, entanglement negativity and discord, in the Bunch-Davies or more general $\alpha$-vacua (e.g. [30, 31, 32] and references therein). We further refer our reader to e.g. [33, 34, 35, 36, 37, 38, 39, 40, 41] for discussions on various aspects of quantum correlation in dS, including their possible observational consequences.

Most of the references cited above focus on the scalar field theory and discussions on other spin fields seem sparse. In [42], the entanglement properties of a Dirac fermion in the hyperbolic dS was discussed and certain qualitative differences with a scalar field were pointed out. See [39, 43, 44, 45] for interesting aspects of fermionic entanglement in cosmological spacetimes. See also [46, 47] for discussions on fermionic entanglement in the Rindler spacetime.

In this work we wish to point out a further qualitative difference of the fermionic field theory in the hyperbolic dS from that of the scalar [7], which seems to have been missed in the earlier literature, as follows. After reviewing the construction of the local orthonormal modes and the global ones in 2, we show in 3 that those global modes are not orthogonal, as opposed to the scalar field theory. It follows from a simple and generic result of the canonical quantum field theory (see [48, 49] and references therein), that if one attempts to do field quantisation with such non-orthogonal global modes, the resulting Bogoliubov structure would *not * preserve the desired anti-commutation relations for the creation and annihilation operators corresponding to the field quantisation of the global modes. Thus we need to orthonormalise those global modes before making any sensible field quantisation. We argue in 3 that such orthonormalisation is never unique a priori in the present scenario, giving rise to a continuous, one parameter family of global vacua. In other words, we wish to demonstrate a natural emergence of de Sitter $\alpha$-like vacua for Dirac fermions in the hyperbolic de Sitter spacetime. This is the main result of this paper. Using such vacua, we investigate next the fermionic entanglement properties in 4.

We shall work with the mostly positive signature of the metric in $(3+1)$-dimensions and will set $c=G=\hbar=\kappa_{B}=1$ throughout.

2 The Dirac modes in hyperbolic dS

In this section we shall review the basic geometry of the hyperbolic dS, the local Dirac modes and their analytic continuation to form the global modes. The detail of the following can be seen in [8, 9, 42, 50, 51].

The dS with hyperbolic spatial slicing is obtained by analytic continuation of the Euclidean sphere $S^{4}$ onto the Lorentzian sector [8, 9]. Such analytic continuation gives rise to three distinct and causally disconnected regions of the global dS spacetime, say $R$, $L$ and $C$, as shown in 1. The relationship between the Lorentzian $(t,r)$ and Euclidean coordinates $(\tau,\rho)$ of these three regions is given by,

[TABLE]

The metrics in the three regions read,

[TABLE]

[TABLE]

where $d\Omega^{2}$ is the metric on ${\cal S}^{2}$ and $r$ and $t$ are dimensionless in the units of $H$.

We shall briefly review now the solutions of the free Dirac equation in our regions of interest $R$ and $L$ [42, 50, 51]. In the region $R$, the four Dirac modes are given by,

[TABLE]

where $z_{R}=\cosh t_{R}$ and a ‘star’ denotes complex conjugation. $\chi^{(\pm)}(r,\theta,\phi)$ are two orthonormal spatial 2-spinors defined on the 3-hyperboloid, satisfying the eigenvalue equation

[TABLE]

The ‘tilde’ denotes differentiation over the 3-hyperboloid and since it is an open manifold, the eigenvalue $p$ is continuos and positive. The temporal parts $u_{p}$, $v_{p}$ appearing in 8 are given by the hypergeometric functions

[TABLE]

The four modes in the region $L$, consistent with 2 are given by,

[TABLE]

where as earlier $u_{p}(z_{L})$ and $v_{p}(z_{L})$ are given by 15 with $z_{L}=\cosh t_{L}$.

It is easy to see that in the asymptotic past, $t_{R}\to 0$ or $t_{L}\to 0$, each of the sets 8 and 20 splits into two positive frequency (representing particles) and two negative frequency (representing anti-particles) modes.

The definition of the Dirac inner product between any two modes $\Psi_{(a)}$ and $\Psi_{(b)}$, is given by

[TABLE]

Since the local mode functions have supports only in their respective regions, the sets given in 8 and 20 are trivially mutually orthogonal. Also, since the inner product is independent of time, one can take, in order to simplify the algebra, the Cauchy surface of the integration to be either the $t_{R}=0$ or the $t_{L}=0$ hypersurface in the relevant region in 1. Using appropriate simplifications of the hypergeometric functions in this limit [52] and also by using the orthonormality of $\chi^{(\pm)}$, it is easy to check that the modes in $R$ or $L$ are orthonormal, i.e.,

[TABLE]

with $s,s^{\prime}=\pm$ and $(\Psi_{(s)}^{R},\,\Psi_{(s^{\prime})}^{L})=0$, trivially.

Note also that since we are interested in constructing field theory in $R\cup L$ and the region $C$ is causally disconnected from them, we shall not be concerned about the modes in $C$ in this work. We further refer our reader to [42, 53] on the normalisability of the mode functions for Dirac spinors and massive vectors in that region.

In order to understand the quantum correlations of fields located in $R$ and $L$, we also require, apart from the above local modes, the notion of global modes which have support in $R\cup L$. First noting that the functions in 15 have branch points at $z=\pm 1$ and $\infty$, using 2 and some identities of the hypergeometric function [52], such global modes are achieved by analytically continuing the local modes of $R$ to $L$ or of $L$ to $R$ along a complex path going through $C$ [42]. We review this procedure in A. The resulting global mode functions originating from the $R\to L$ analytic continuation are given by,

[TABLE]

where

[TABLE]

Likewise we have for the $L\to R$ analytic continuation,

[TABLE]

Since each of the above mode functions have supports in both the regions, when normalised, they are supposed to be the global versions of the local modes appearing in 8, 20. However, we shall see below that the modes of 32 and 53 do not form an orthogonal set under the global inner product, in contrast to the scalar field theory [8]. As we have discussed in 1, we cannot simply treat such non-orthogonal modes as our global basis modes in the field quantisation [48, 49]. Hence we first need to orthonormalise the modes of 32, 53.

3 Constructing the global orthonormal modes

The normalisation integration for the global modes looks formally similar to that of the local ones 26, with the difference that the integration hypersurface now must exist in $R\cup L$. Following [8], we choose it to be the $(t_{L}=0)\cup(t_{R}=0)$ Cauchy surface for convenience. Then e.g., for the pair $\{\Psi^{+L}_{(-)}$, $\Psi^{-R}_{(-)}\}$ we have,

[TABLE]

The suffix ‘G’ stands for global, whereas the inner products on the right hand side are local. We can express $\Psi^{-R}_{(-)}$ in terms of $z_{L}$ (for the first term on the right hand side) and $\Psi^{+L}_{(-)}$ in terms of $z_{R}$ (for the second term on the right hand side) respectively via 32 and 53. Then using 27, we find (after suppressing the various $\delta$-functions for the sake of brevity)

[TABLE]

where $\lambda_{2}$ is given by 48. It can further be checked in the similar manner that the eight global fermionic modes of 32, 53 can be grouped into four pairs such that the members of any pair are not orthogonal (although inter-pair orthogonality is satisfied) with respect to the global inner product, given by

[TABLE]

Thus evidently we cannot use these modes as our global basis modes, for they would lead to a non-preservation of the canonical anti-commutation relations when used as basis of field expansion, e.g. [48, 49]. Hence we need to make suitable linear combinations between the members of each pair of 69, to orthogonalise them via standard Gram-Schmidt procedure.

Now, we note from 32, 53 that in doing so, we are basically superposing positive and negative frequency modes. Thus in order to accommodate sufficient generality a priori in our orthogonalisation scheme, we must treat both the solutions simultaneously in an equal footing. This can be achieved by introducing a continuous parametrisation to obtain two orthogonal global modes,

[TABLE]

where $\Delta\theta_{1}$ and $\Delta\theta_{2}$ are one-parameter functions given by

[TABLE]

and

[TABLE]

The ‘angle’ $\theta_{\rm RL}$ does not depend upon any spacetime coordinate. It is easy to check that the two mode functions defined in 70 are indeed orthogonal under the global inner product.

In the same spirit, by choosing appropriate linear combinations for the three other pairs in 69 one can generate orthogonal pairs of global mode functions. The full set of eight orthonormal global modes is given by,

[TABLE]

where we have written for the normalisations,

[TABLE]

Using 69, the orthonormality of these modes under the global inner product can explicitly be checked at once. To the best of our knowledge, the above issue of orthonormalisation for the fermionic global modes in hyperbolic dS was not addressed in the earlier literature.

Note that the quantisation of the Dirac field with the modes of 3 will effectively thus give de Sitter $\alpha$-vacua like structure (see e.g. [30, 31, 32, 54]). However, we emphasise that unlike the usual cases of such vacua, introducing the parametrisation $\theta_{\rm RL}$ was a priori necessary in our current scenario, in order to maintain sufficient generality in the orthogonalisation procedure. This is the main result of this work. We shall see later that $\theta_{\rm RL}=0,\pi/2$ values correspond to the usual thermal distribution of created massless particles in $R$ or $L$.

4 Computation of the entanglement and the Rényi entropies

4.1 The field quantisation and the Bogoliubov coefficients

Let us first make a field quantisation in terms of the local modes of 8, 20,

[TABLE]

where $s=\pm$. The operators are postulated to satisfy the anti-commutation relations

[TABLE]

and all other anti-commutators vanish. We define the local vacua $|0\rangle_{R},|0\rangle_{L}$,

[TABLE]

Likewise, we can also expand the Dirac field in terms of the orthonormal global modes, 3. How do we identify the creation and annihilation operators here? Recalling scalar field’s case [8], we note from 48 that in the limit $p\to\infty$, both $\lambda_{1}$ and $\lambda_{2}$ are vanishing, showing in this limit we do not have any analytically continued modes in 32 and 53 and accordingly, 3 reduces to the local mode functions, having well defined positive or negative energy characteristic in the asymptotic past. Since in that limiting scenario we do not have any trouble with identifying the creation and annihilation operators, we may make the following expansion of the field $\Psi$ in terms of the global modes,

[TABLE]

where we interpret $a_{1},\dots,a_{4}$ and $b_{1},\dots,b_{4}$ as the annihilation operators related to the global modes. The global vacuum is defined as

[TABLE]

We now equate 74 and 4.1, and successively take eight inner products with both sides with respect to the eight global basis modes, 3. While doing so, we need to use 32 or 53 in order to express the global modes in terms of the local ones. For the global mode $\Psi_{1pjm}$ for example, we obtain

[TABLE]

where we have written

[TABLE]

and have suppressed the eigenvalues $p,j,m,$ for the sake of brevity. Similarly, by taking the inner products with the seven other modes $\Psi_{2},\,\Psi_{3}\dots,\Psi_{8}$, we obtain

[TABLE]

where we have written

[TABLE]

Using now 4.1 it can be checked that the global operators satisfy the desired anti-commutation relations,

[TABLE]

where $\sigma=1,\,2,\,3,\,4$ with all other anti-commutations vanishing.

We once again emphasise here that had we not properly orthonormalised our global modes, we would not have retained the above canonical anti-commutation structure, essential to preserve the spin-statistics of the field theory.

Now, by observing the right hand sides of 79, 80, it becomes evident that the set of global operators ($a_{1}$, $a_{3}$, $b_{2}$, $b_{4}$) and ($a_{2}$, $a_{4}$, $b_{1}$, $b_{3}$) form two disjoint sectors, for their operator contents on the right hand side are different. This implies that the global vacuum defined in 78 can also be split into two subspaces,

[TABLE]

where $|0\rangle^{1}$ is defined by $a_{2}|0\rangle^{(1)}=a_{4}|0\rangle^{(1)}=b_{1}|0\rangle^{(1)}=b_{3}|0\rangle^{(1)}=0$, whereas $|0\rangle^{(2)}$ is such that $a_{1}|0\rangle^{(2)}=a_{3}|0\rangle^{(2)}=b_{2}|0\rangle^{(2)}=b_{4}|0\rangle^{(2)}=0$. For the rest of this paper, we shall work with $|0\rangle^{(1)}$ only. As long as we are only concerned with the computation of the entanglement and Rényi entropies, the other subspace, $|0\rangle^{(2)}$, will produce identical results.

From 79, 80, it is clear that $|0\rangle^{(1)}$ can be constructed as a squeezed state over the local vacua defined in 76,

[TABLE]

where $\xi_{ij}$’s are four complex numbers and $N$ is the normalisation. Also, since the operators $c$’s and $d$’s anti-commute (4.1), we may further decompose the vacua defined in 76 as,

[TABLE]

where $|0_{c}\rangle_{R},~{}|0_{d}\rangle_{R},~{}|0_{c}\rangle_{L},~{}|0_{d}\rangle_{L}$ are annihilated respectively by $c_{(+)}^{R},~{}d_{(+)}^{R},~{}c_{(+)}^{L},~{}d_{(+)}^{L}$.

We now expand the exponential in 83 keeping in mind the various anti-commutations and then annihilate $|0\rangle^{(1)}$ by $a_{2}$, $a_{4}$, $b_{1}$, $b_{3}$, to obtain,

[TABLE]

Also, the normalisation in 83 reads,

[TABLE]

Putting these all in together, we may now explicitly write down the global vacuum in terms of the local states,

[TABLE]

The above expression shows that there will be pair creation in $R$ and $L$ with respect to the global vacuum. Since the states belonging to $R$ and $L$ cannot be factored out, those pairs will be entangled. Also, since $|0\rangle^{(1)}$ depends upon $\theta_{\rm RL}$ through $\xi_{1}$ and $\xi_{2}$, the particle creation and the entanglement will also depend upon it.

Let us then first compute, as a check of consistency of the entire framework, the expectation value of the local number operator with respect to the global vacuum, which will give us the number of created particles at a given mode. We shall do it in the massless limit only, for which from 84 and 86 we have,

[TABLE]

For $\theta_{\text{RL}}\to 0,\pi/2$, we reproduce the familiar Fermi-Dirac distribution,

[TABLE]

Away from these values of $\theta_{\rm RL}$, the spectrum is not thermal. For $\alpha$-vacua in the dS spacetime with flat spatial slicing, such non-thermality was also noted earlier in [32].

4.2 The entanglement entropy

We start with the full density operator $\rho_{\rm global}=|0\rangle^{(1)}\,{}^{(1)}\langle 0|$ and trace over the states of, say $L$ region, inaccessible to an observer in $R$. We thus obtain the reduced density operator in $R$,

[TABLE]

Tracing over the $R$ states will give identical results. We find a matrix representation of $\rho_{R}$ in terms of the $R$-state vectors,

[TABLE]

The entanglement entropy for a single mode (with a given $p$ and $\theta_{\rm RL}$ value) is given by

[TABLE]

where $\lambda_{R}^{(i)}$’s are the eigenvalues of 93. The entanglement entropy per unit volume between $R$ and $L$ is obtained by integrating over all $p$ modes. The final expression of the entanglement entropy is obtained by further integrating the result over the purely spatial section of 2. Since $S(p,m;\theta_{\rm RL})$ has no spatial dependence, the integration, being done over a non-compact space, would diverge. Accordingly, one needs to put a cut off at some ‘large’ radial distance. The resultant regularised volume integral equals $2\pi$, see [7] for details. The regularised entanglement entropy then equals

[TABLE]

where ${\cal D}(p)=(1/4+p^{2})/(2\pi^{2})$ is the apprpriate measure in the momentum space corresponding to the spatial section of 2 for the spin-1/2 field [51, 55]. The integral in 95 is convergent and can be calculated numerically.

We have plotted various characteristics of the entanglement entropy in 2 and 3 as a function of the parameter $\nu^{2}=9/4-m^{2}/H^{2}$ subject to different values of $\theta_{\rm RL}$. We may chiefly note the following features.

a) In 2, the curves corresponding to $\theta_{\rm RL}=0,\,\pi/2$, are exactly coincident and they show maximal $R-L$ entanglement for all values of $\nu^{2}$. The coincidence corresponds to the fact that the coefficients $\xi_{1,2}(\theta_{\rm RL}=0)=\pm\xi_{1,2}(\theta_{\rm RL}=\pi/2)$, in 84. b) While most of the curves in 2 are monotonic, the curve corresponding to $\theta_{\rm RL}=\pi/3$ shows extrema. c) For any given value of $\theta_{\rm RL}$, the entanglement entropy has its maximum value in the massless limit, $\nu^{2}=9/4$. This might correspond to the fact that in this limit the creation of particle-antiparticle pairs is energetically most favourable. Since such pairs are entangled, 86, it is perhaps reasonable to expect that the entanglement entropy also gets its maximum value in the massless limit.

4.3 The Rényi entropy

Before we end, we wish to briefly discuss the Rényi entropy, a one-parameter generalisation of entanglement entropy, e.g. [29]

[TABLE]

For $q\to 1$, 96 reduces to the entanglement entropy, as can be easily seen by using the L’Hopital’s rule to the above equation.

The final Rényi entropy, akin to the expression for the final regularised entanglement entropy, is given by,

[TABLE]

We have plotted $S_{q}(m,\theta_{\text{RL}})$ in 4(a), 4(b) as earlier with respect to the parameter $\nu^{2}=9/4-m^{2}/H^{2}$, with different values of $\theta_{\rm RL}$. We note chiefly that as a whole, the qualitative nature of the Rényi entropy for different $q$ values remains the same as that of the entanglement entropy. In particular, a) the values $\theta_{\rm RL}=0,\pi/2$ gives maximum Rényi entropy for all values of $\nu^{2}$ and b) the extrema for $\theta=\pi/3$ is still present.

5 Discussions

In this paper we have investigated the entanglement and the Rényi entropies between the $R$ and $L$ states of the Dirac field in the hyperbolic dS spacetime, 2, 3, as a measure of the long range non-local quantum correlations between these two regions.

The chief results of these paper could be summarised as follows. First, the natural emergence of the continuous, one parameter family of global modes (cf. 3) and vacua, 86. Such vacua have structures similar to the de Sitter $\alpha$-vacua, though they have originated here from the mere necessity of an a priori general orthonormalisation scheme for the global modes. Such orthonormalisation of the modes is necessary to preserve the canonical anti-commutation structure of the field theory [48, 49]. Second, we have seen in 4.2, 4.3 that a) $\theta_{\rm RL}=0,\pi/2$ reproduces the thermal spectrum for the created massless particles. The dependence of the entanglement and the Rényi entropies on the parametrisation $\theta_{\rm RL}$ was depicted in 2, 3, 4(a), 4(b).

We recall that instead of taking $\alpha$ as a usual independent parameter, one can also take it to be momentum dependent, see [30] for a discussion on the scalar field theory. Note that even though we have not taken any momentum dependence in $\theta_{\rm RL}$ here, the coefficients of the linear combinations in the global modes, 3, are indeed momentum dependent. This effectively makes our construction qualitatively similar to that of [30]. Nevertheless, it will be interesting on top of this to further allow explicit momentum dependence in $\theta_{\rm RL}$. In this case we need to make a suitable ansatz for it, such that the mode by mode normalisability of the states is achieved and also the various momentum integrals we encounter converge.

It seems interesting to investigate the effects of $\theta_{\rm RL}$ into the other measures of quantum correlations e.g., the entanglement negativity, the violation of Bell’s inequality and the quantum discord etc, in order to quantify it further. We hope to address these issues in our future work.

Acknowledgements

We would like to acknowledge R. Basu, S. Chakraborty, R. Gupta, A. Lahiri, S. Mukherjee, B. Sathiapalan and I. S. Tyagi for fruitful discussions. SB’s research is partially supported by the ISIRD grant 9-289/2017/IITRPR/704. SC is partially supported by ISIRD grant 9-252/2016/IITRPR/708.

Appendix A The analytic continuation and the global modes

Here we review the construction of the global modes found via the analytic continuation of the local mode functions of 8, 20, [42]. Since the hypergeometric function has branch points at $1$ and $\infty$, the functions $u_{p}(z_{R})$ and $v_{p}(z_{R})$ in 15 have branch points at $\pm 1$ and $\infty$. We join $\pm 1$ with a cut on the real line and further join either of them to $\infty$ in any possible way we wish, ensuring however, that we do not hit the cut to $\infty$ during the process described below. The spatial functions $\chi^{(\pm)}$ however, do not undergo any formal changes under this procedure.

We divide the real $z$ line into three regions : a) $z\geq 1$ ($z_{R}$ or the $R$ region) b) $-1\leq z\leq 1$ (the $C$ region) and c) $z<-1$. Note that both $z_{R}$ and $z_{L}$ appearing in the mode functions are greater than unity. Thus when we analytically continue an $R$ mode to the $L$ region via $C$, we eventually reach at negative $z$ values. Hence while continuing analytically, we shall take $z_{R}=-z_{L}$ with $z_{L}>1$.

We now extend $u_{p}(z_{R})$ from $R$ to $C$ by taking $z_{R}-1=e^{-\pi i}(1-z_{R})$ and change the variable by $z_{R}=-z_{L}$, as mentioned above. Since $z_{L}>1$, we perform $1-z_{L}=e^{\pi i}(z_{L}-1)$, to obtain the continued form of $u_{p}(z_{R})$,

[TABLE]

In order to cast the above expression into the form of our initial mode functions, we recall the identity, e.g. [52]

[TABLE]

Substituting this into 98 with $z\equiv(1+z_{L})/2$ and also using

[TABLE]

we finally obtain

[TABLE]

where $\lambda_{1}$, $\lambda_{2}$ are given by 48.

Similarly we find,

[TABLE]

The opposite procedure, i.e. the analytic continuation $L\to R$ yields formally exactly similar results.