abstract: The band structure of iron-based superconductors gives rise to yet another scenario for the appearance of Dirac fermions. A viewpoint on "Observation of Dirac cone electronic dispersion in BaFe$_2$As$_2$": (Richard et.al., PRL 104, 137001 (2010)).
author: M. Zahid Hasan; B. Andrei Bernevig
date: 2024-10-02
title: Dirac cone in iron-based superconductors

Superconductivity above 30 K was realized in the mid 80s with the discovery of cuprate high-temperature superconductors . The monopoly of copper was broken in 2008 with the discovery of iron-based superconductors - of which there are now many families - which have a maximum transition temperature (55 K) comparable to that of single-layer cuprates . Their complex phase diagrams are poorly understood, but we do know that most have a magnetic phase that changes to a superconducting one as the concentration of electrons or holes in the bulk increases.

The parent magnetic state of iron-based superconductors determines their superconducting transition temperature. In their paper in *Physical Review Letters* , Pierre Richard and collaborators from China, Japan, and the US report band-structure measurements of the parent magnetic state of iron-based superconductors that show linearly dispersed bands in the shape of a cone, which implies that the carriers follow dispersion relations for massless relativistic fermions that obey the Dirac equation. These results confirm an existing theoretical prediction and earlier experimental results on the Dirac band structure and the orbital nature of the magnetic order in the parent superconductor . Richard *et al.* present angle-resolved photoemission (ARPES) studies of high-quality single crystals of BaFe$_2$As$_2$, which has a magnetic transition temperature of 138 K when undoped. They map the Fermi surfaces and find slightly anisotropic band crossings arranged in the shape of a Dirac cone. Dirac fermions are also known to exist in other condensed-matter systems such as in the bulk and surface of topological insulators , in graphene , and cuprate superconductors .

In iron-arsenide compounds, the parent state is a collinear antiferromagnet, with a commensurate ordering wave vector ($\pi$, 0) or (0, $\pi$). This wave vector also appears in the band structure of the normal state of iron-based superconductors as a momentum-space vector (nesting vector) connecting the Fermi surface's hole pockets that surround the $\Gamma$ point and the electron pockets around the $M$ point in the unfolded Brillouin zone (see Fig. 1). This band structure is a consequence of the multiorbital nature of the system; all of iron's five 3$d$ orbitals are necessary for a full electronic description of the system. Of these, the most important for understanding the present experimental results are the *d$_{xz}$, d$_{xy}$, and d$_{yz}$* orbitals.

Plausible explanations for the low-temperature magnetic ordering vector involve either a spin-density wave (SDW) formation tendency consistent with the nesting wave vector or an antiferromagnetic state formation tendency due to next nearest neighbor spin-spin interactions . In the magnetic phase, the Fermi surfaces reorganize, the Brillouin zone folds, and gaps due to electron-hole scattering are seen in the electronic spectrum. SDW formation usually leads to the opening of an electronic gap by destroying the Fermi surface if the Bloch vectors along the nesting vectors can span the full Fermi surface. Strong interactions can also gap the SDW state. In the iron-based compounds, as the present experiment shows, the situation is different: due to the topology of the Fermi surface a coupling between the Bloch wave vectors of the electron and hole pockets sets the SDW-induced electronic gap.

In accounting for the nature of the SDW state, a full theoretical analysis \[6\] based on a five-orbital model shows that the Bloch band of the electron Fermi surface pocket is odd under reflection, while the outer hole pocket of the Fermi surface is even along the SDW ordering axis. The SDW matrix element between them vanishes and gives rise to a nodal (not fully gapped) SDW, which can at most be semimetallic. This gapless SDW magnetic state survives even for strong interactions . In experiments, the authors observe a high-intensity point in the ARPES data on the $\Gamma$?M symmetry line for the folded band structure of magnetically ordered BaFe$_2$As$_2$. This feature, occurring roughly where the electron band hybridizes with a hole band, disappears above the SDW transition temperature and has a dispersion characteristic of a relativistic Dirac cone. The observed Dirac cone is thus related to the nature of the SDW state in these materials. The apex of the Dirac cone is situated slightly above the Fermi energy, which would give rise to small hole pockets in the ARPES data. Its dispersion has a small degree of anisotropy. Away from the cone the band is gapped and the authors' data provide the variation of the gap function. A fourfold symmetric electron density pattern is observed around the $\Gamma$ and M points. This property could be intrinsic, such as in a weak-SDW mode, or arise from the superposition of twin domains expected to form under the structural distortion and SDW transitions. The latter scenario was originally proposed in Ref. .

The occurrence of Dirac fermions in the system implies that the SDW magnetic state is always metallic. The experimental observations, however, cannot distinguish between a weak-coupling SDW and an antiferromagnetic picture in which some more localized bands are also present. The discovery of the Dirac cone shows that theoretical models that neglect orbital symmetries of iron-based superconductors may not explain the important aspects of the physics. More experiments are needed to address several puzzles. The 122 materials are known to have rather strong $z$-axis dispersion and the evolution of the Dirac cone with the $z$-axis momentum needs to be studied. In several materials, the SDW state coexists with the superconducting state. One possibility for the superconducting order parameter is that of the sign changing s$^\pm$-wave type, which is nodeless. If the superconductor is indeed of the nodeless s$^\pm$ type, does the gapless Dirac fermion still survive in the phase where SDW and superconductivity co-exist? If it does not, then could it be that the superconducting gap opening of the Dirac cone gives rise to a superconductor with topological properties ? The observation of a Dirac fermion in iron-based materials introduces a new system with protected cones, after cuprate superconductors , graphene , and topological insulators . It is likely to inspire new research and unravel the connection between nodal SDW magnetism and the nature of superconductivity.

Iron-based superconductors are the latest example of the Dirac physics that has recently energized condensed matter physics. In cuprate superconductors, Dirac nodes (points where the two cones meet) appear in the structure of the superconducting order parameter, caused by the electron-electron interactions in the parent magnetic Mott insulator. In graphene, degenerate Dirac fermions appear in the noninteracting band structure due to the ${C_3}$ symmetry of the lattice, along with inversion and time-reversal symmetry. In topological insulators, spin-polarized Dirac fermions appear in the strong spin-orbit coupled band structure at the edge of the insulating bulk and enjoy extended protection due to time-reversal symmetry . In iron-based superconductors, Dirac fermions appear in the magnetic SDW state band structure as a result of the point-group symmetry of the orbitals making up the Fermi surfaces involved in the nesting process. Even though they enjoy different degrees of protection against disorder and interactions, the Dirac fermions present in all these materials are fascinating examples of how different systems can lead to similar profound low-energy electron behavior.


abstract: We present a method to semiclassically compute the pair creation rate of bosons and fermions in de Sitter spacetime. The results in the bosonic case agree with the ones in the literature. We find that for the constant electric field the fermionic and bosonic pair creation rate are the same. This analogy of bosons and fermions in the semiclassical limit is known from several flat spacetime examples.
author: Clément Stahl; Strobel Eckhard
date: 1 June 2015
title: Semiclassical fermion pair creation in de Sitter spacetime

# Introduction

Quantum field theory (QFT) in curved spacetime is one way of merging Einstein's theory of gravitation and QFT in the usual Minkowski spacetime within a self consistent framework . One of its central results is the discovery of particle creation in a time-dependent gravitational field . This mechanism is believed to generate the primordial cosmic inhomogeneities that serve as seeds for the observed large-scale structure of the universe.
A lot of studies of QFT in curved spacetime are investigating de Sitter (dS) spacetime. It has constant scalar curvature and is maximally symmetric for any given dimension, in the same way as Minkowski spacetime. In cosmology, dS space is used as a model for both the early stage of inflation (for reviews see ) and the late stage of acceleration of the expansion (an introductory review is given in ). In addition to its relevance for cosmology, dS space may give hints to understand the quantum nature of spacetime . These facts motivate the investigation of physical processes in this spacetime.
Originally, particle creation was studied under a strong field background originating from a time dependent vector potential . Since then, it has been explored in a more general set-up, where the effects of the gravitational field and the electric field are both contributing to the creation of pairs. The Schwinger mechanism in curved spacetime, more specifically, in dS space, has recently started to be studied in depth .
Better understanding Quantum electrodynamics (QED) in dS spacetime could be insightful for a couple of problems. It is an interesting framework for the study of false vacuum decay or bubble nucleation . Constrains on magnetogenesis scenarios were put in via the backreaction of the created pairs and its induced current. In some models of preheating, it might give clues on the open problem of baryogenesis. Via the AdS-CFT correspondence, it has also been used as a playing field to test the ER=EPR conjecture . Additionally it might help to better understand renormalization schemes in curved space-time and the relations of these schemes to each other.
The purpose of this paper is to investigate spin 1/2 pair production under the influence of an external electric field, in four-dimensional dS spacetime ($\text{dS}_4$). This problem has already been studied for bosons in two-dimensional dS spacetime ($\text{dS}_2$) and in $\text{dS}_4$ . In spin 1/2 pair production in $\text{dS}_2$ was investigated but the 4D analog case is much less studied. We propose a derivation of the number of pairs created in the semiclassical limit following the ideas of . We find that there is no difference between bosons and fermions in this limit.
The paper is organized as follows: in section we investigate the Schwinger effect in $\text{dS}_4$. After describing the basic equations, we show how the Dirac equation becomes a coupled second order differential equation in the presence of both electric and gravitational terms. We review the derivation of the bosonic pair creation rate and calculate the fermionic one with the help of a semiclassical expansion for general electric fields in dS spacetime in section . We use this result in section to compute the bosonic and fermionic pair creation rate of a constant electric field in dS spacetime. Finally we draw some conclusions in section .

# Dirac field in dS Space-time

We consider QED coupled to a Dirac fermion in $\text{dS}_4$. In order to study Schwinger pair production, we assume that the gravitational and electrical field are background fields and that the fermionic field is dynamic. The action is given by $$\begin{aligned}
\label{action}
S=\int \text{d}^4 \text{x} \sqrt{-g} \mathcal{L}=\int \text{d}^4 \text{x} \sqrt{-g} \left[  -\frac{1}{\kappa}R +\frac{i}{2} (\bar{\psi} \gamma^{\mu} \nabla_{\mu} \psi -\nabla_{\mu} \bar{\psi}\gamma^{\mu}\psi) -m \bar{\psi} \psi  -\frac{1}{4} F_{\mu \nu}F^{\mu \nu} \right],
\end{aligned}$$ where the field strength is defined in the usual way $F_{\mu \nu}\equiv \partial_{\mu} A_{\nu}-\partial_{\nu} A_{\mu}$.
The $\text{dS}_4$ spacetime we want to study is described by the metric $$\text{ds}^2=a(\eta)^2 (\text{d}\eta^2-\text{d}\textbf{x}^2),$$ with signature (+,-,-,-). Here $\eta$ is the conformal time which is parametrized by the Hubble factor in the following way $$\eta= -\frac{1}{a(\eta)H},\hspace{1cm} a(\eta)^2 H \equiv \frac{\text{d}a(\eta)}{\text{d} \eta }, \hspace{1cm} (-\infty<\eta<0).$$ For the description of spinors the tetrad field is used. It is related to the metric through the relation $$\begin{aligned}
g_{\mu \nu}=e^a_{\mu} e^b_{\nu} \eta_{ab},
\end{aligned}$$ where $\eta_{ab}$ is the usual Minkowski metric. Throughout this paper, we use Greek indices for spacetime indices ($\mu,\nu=\eta,x,y,z$) and Latin indices for tetrad ones ($a,b=\eta,x,y,z$). Applying the tetrad formalism to $\text{dS}_4$, one gets $$\label{tetrad} e_{\mu}^a=
\begin{pmatrix}
  a(\eta) & 0 & 0 &0  \\
  0 & a(\eta) & 0 & 0 \\
  0 & 0 & a(\eta) & 0 \\
  0 & 0 & 0 & a(\eta)
 \end{pmatrix}.$$ The covariant derivative for fermion fields is defined as $$\begin{aligned}
  \nabla_{\mu} \equiv \hbar\partial_{\mu} +ie A_{\mu}(x) -\frac{{\text{i}}}{4}\omega^{ab}_\mu\sigma_{ab},
\end{aligned}$$ with the commutator of the gamma matrices $$\begin{aligned}
 \sigma_{ab}&\equiv\frac{{\text{i}}}{2}[\gamma_a,\gamma_b]
\end{aligned}$$ and the spin connection are defined as $$\begin{aligned}
\begin{split}
\omega_{\mu}^{ab}&\equiv \frac{1}{4}\left[e^{b\alpha}(x) \partial_{\mu}e^a_{\alpha}(x)-e^{a\alpha}(x) \partial_{\mu}e^b_{\alpha}(x)+ e^{a \alpha}(x) \partial_{\alpha}e^b_{\mu}(x)-e^{b\alpha}(x) \partial_{\alpha} e^{a \mu}(x)\right.\\ & \hspace{2cm}\left.+e^{b\nu}(x)e^{a \lambda}(x)e_{c \mu}(x)\partial_{\lambda} e^c_{\nu}(x)-e^{a\nu}(x)e^{b \lambda}(x)e_{c \mu}(x)\partial_{\lambda} e^c_{\nu}(x) \right]. \label{eq:spinconnection}
\end{split}
\end{aligned}$$ The non-zero components of the spin connection () can be shown to be $$\begin{aligned}
 \label{connection}
 \omega_{1}^{01}=\omega_{2}^{02}=\omega_{3}^{30}=-\omega_{1}^{10}=-\omega_{2}^{20}=\omega_{3}^{30}=\frac{a'(\eta)}{2a(\eta)},

\end{aligned}$$ where prime denotes derivative with respect to conformal time $\eta$.
The gamma matrices $\gamma^{\mu}$ are related to the gamma matrices in the tangent flat space $\underline{\gamma^a}$ viz. $$\begin{aligned}
\underline{\gamma^a} \equiv \gamma_{\mu} e^a_{\mu}.
\end{aligned}$$ We will work with the Dirac representation of the gamma matrices, *i.e.* $$\begin{aligned}
 \gamma^{j}=\begin{pmatrix}
             0 &\sigma^j\\
             -\sigma^j & 0
            \end{pmatrix},
&&
  \gamma^{0}=\begin{pmatrix}
             I_2 &0\\
             0 & -I_2
            \end{pmatrix},
&& \text{ where }&&
   \sigma^x=\begin{pmatrix}
             0 & 1\\
             1 & 0
            \end{pmatrix},
&&
   \sigma^y=\begin{pmatrix}
             0 & -{\text{i}}\\
             {\text{i}}& 0
            \end{pmatrix},
            &&
   \sigma^z=\begin{pmatrix}
             1 & 0\\
             0 & -1
            \end{pmatrix}.
\end{aligned}$$ Varying the action with respect to the spinor field gives the Dirac equation $$\left(i \gamma^{\mu} \nabla_{\mu} -m \right) \psi(\textbf{x}, \eta) =0.$$ Using Eqs. () and (), this equation becomes $$\label{eq:diracavantcahnge}
\left\{i \left(\gamma^{\mu} \hbar\partial_{\mu}+\frac{3}{2} a H \underline{\gamma^0} + i e A_{\mu} \gamma^{\mu} \right) -m \right\} \psi(\textbf{x}, \eta)=0.$$ One now considers the auxiliary field $\Psi(\textbf{x}, \eta) = a^{3/2}(\eta) \psi(\textbf{x}, \eta)$ which can be thought of as the equivalent of the Mukhanov-Sasaki variable in inflation models . With this substitution the Dirac equation takes the form $$\left\{\gamma^{\mu}  ( i \hbar\partial_{\mu} - e A_{\mu})  - m \right\} \Psi(\textbf{x}, \eta)=0. \label{eq:Dirac}$$ We will also decompose this field in momentum modes according to $$\begin{aligned}
 \Psi(\textbf{x},\eta)\sim{\text{e}}^{\frac{{\text{i}}}{\hbar}\textbf{k}.\textbf{x}}\psi_{\textbf{k}}.\label{TF}
\end{aligned}$$ To solve the Dirac equation, for the purpose of the calculation of the pair creation rate, it is often useful to use the squared version of the Dirac equation because of its similarities to the Klein-Gordon equation, see e.g. . The squared Dirac equation can be found using $$\begin{aligned}
  \Psi(\textbf{x},\eta)&=\gamma^{\mu}\left[(i \partial_{\mu}-e A_{\mu}(\eta)) +m a(\eta)\right]\phi(\textbf{x},\eta) \label{eq:squared}

\end{aligned}$$ with $$\begin{aligned}
 \phi_{\textbf{k}}(\eta)&=\begin{pmatrix}\phi_1(\eta)\\  \phi_2(\eta) \end{pmatrix},&&  \phi_i(\eta)=\begin{pmatrix} \phi_i^+(\eta)\\  \phi_i^-(\eta)\end{pmatrix}.

\end{aligned}$$ in the Dirac equation. In the previous equation, $\phi_{\textbf{k}}(\eta)$ is the Fourier transform (in the sense of ()) of $\phi(\textbf{x},\eta)$. Here we consider a background vector potential for the electromagnetic sector such that $$A_{\mu}\equiv A(\eta)\delta^{z}_{\mu}.$$ For such fields the squared Dirac equation takes the form $$\begin{aligned}
&\left(\hbar^2 \partial_{\eta}^2+\omega_{\textbf{k}}(\eta)^2-i\hbar ma'(\eta) \right)\phi_1^{\pm} \pm i\hbar e A'(\eta) \phi_2^{\pm}(\eta)=0, \\
&\left(\hbar^2 \partial_{\eta}^2+\omega_{\textbf{k}}(\eta)^2+i\hbar ma'(\eta) \right)\phi_2^{\pm} \pm i\hbar e A'(\eta) \phi_1^{\pm}(\eta)=0,   \label{to solve}
\end{aligned}$$ where the effective pulsation and the kinetical momentum are defined as $$\begin{aligned}
& \omega_{\textbf{k}}(\eta)^2\equiv p_z(\eta)^2+k_{\perp}^2+m^2 a(\eta)^2, \label{eq:omega}
& p_z(\eta)\equiv k_z+e A(\eta),
&& k_\perp^2\equiv k_x^2+k_y^2.
\end{aligned}$$ This equation can be compared to the equivalent bosonic problem. The equation of motion derived from the Klein-Gordon equation is (see Eq (2.13) of ) $$\label{boson to solve}
\left(\hbar^2 \partial_{\eta}^2+\omega_{\textbf{k},\text{boson}}(\eta)^2 \right) q_{\textbf{k}} =0,$$ with $$\begin{aligned}
\omega_{\textbf{k}. \text{boson}}(\eta)^2=\omega_{\textbf{k}}(\eta)^2-\frac{2}{\eta^2}. \label{eq:omegaboson}
\end{aligned}$$ The equation of the bosonic problem () can be understood as a harmonic oscillator with a time dependent pulsation. The two other terms in Eq. () are new and due to the fermionic nature of the particles considered. On the one hand, the mass term was already derived *e.g*. in where no background electric field was considered. On the other hand, the electric term is present for instance in in flat spacetime. Such that this equation is a generalization of the Dirac equation in curved spacetime with background electric and gravitational field.
The squared Dirac equation is also analogous to the Dirac equation for two-component fields in flat spacetime. In a method was used to semiclassically compute the pair creation rate for these fields. It is possible to use the same method for the case studied here. Instead of looking for a solution of the squared Dirac-equation () we will however use the ansatz $$\begin{aligned}
\label{newansatz}
 \psi_{\vec{k},\uparrow}(\eta)=\begin{pmatrix}
                      -{\text{i}}k_\perp \, \psi_2^+(\eta)\\
                       (k_x+{\text{i}}k_y)\, \psi_2^+(\eta)\\
                      -{\text{i}}k_\perp \, \psi_1^+(\eta)\\
                      -(k_x+{\text{i}}k_y)\, \psi_1^+(\eta)\\
                     \end{pmatrix}, &&
   \psi_{\vec{k},\downarrow}(\eta)=\begin{pmatrix}
                      (k_x-{\text{i}}k_y)\, \psi_2^-(\eta)\\
                      {\text{i}}k_\perp \, \psi_2^-(\eta)\\
                      -(k_x-{\text{i}}k_y) \, \psi_1^-(\eta)\\
                      {\text{i}}k_\perp\, \psi_1^-(\eta)\\
                     \end{pmatrix} .
\end{aligned}$$ This ansatz can be derived by finding the solution of the squared equation analogous to and then use () to construct a solution for the Dirac equation (). Observe that $\psi_{\vec{k},\uparrow}(\eta)$ and $\psi_{\vec{k},\downarrow}(\eta)$ are independent since $$\begin{aligned}
 \psi_{\textbf{k},\uparrow}(\eta)^\dagger\cdot \psi_{\textbf{k},\downarrow}(\eta)=0.
\end{aligned}$$ Putting () in the Dirac equation () leads to the equations $$\begin{aligned}
 {\text{i}}\hbar\, \psi_1^{'\pm}(\eta)+ m a(\eta)\, \psi_1^\pm(\eta)\pm (p_z(\eta)+{\text{i}}k_\perp)\,\psi_2^\pm(\eta)=0,\label{to solve 1}\\
 {\text{i}}\hbar\, \psi_2^{'\pm}(\eta)- m a(\eta)\, \psi_2^\pm(\eta)\pm (p_z(\eta)-{\text{i}}k_\perp)\,\psi_1^\pm(\eta)=0,\label{to solve 2}
\end{aligned}$$ which we will solve in section .

# Semiclassical number of pairs in dS spacetime

The semiclassical expansion is considering cases where a vacuum state for the produced particles exits in the asymptotic future. That is true if the background fields are evolving slowly. To quantify the slow varying background more precisely, we introduce a dimensionless slowness parameter $T$ by replacing the scale factor $a(\eta)$ by a family of functions $a_T(\eta) \equiv a(\eta/T)$. Doing that, in the limit of infinitely slow varying backgrounds, $T \rightarrow \infty$, the derivatives of $a(\eta)$ will tend to zero. Orders of $T$ are usually called adiabatic orders but it can be noticed that the only place where $T$ is involved is in the derivative $\partial_{\eta}$ of the Dirac equation so it is possible to "formally" pose $T=1/\hbar$ and expand in power of $\hbar$.
In this section, we will first review the calculation of the number of particles in the bosonic case, which can be carried out with flat spacetime techniques and then present the fermionic case. The strategy will be the following:

-   Reformulate the equation of motion in terms of an equation for the mode functions $\alpha(\eta)$, $\beta(\eta)$.

-   Perform a multiple integral iteration to compute $|\beta|^2$.

-   Calculate the integrals with a semiclassical saddle point approximation to derive the number of created pairs for each momentum mode $\mathbf{k}$.

## Equations for the mode functions

In this section we derive the equations for the mode functions from the Klein-Gordon and Dirac equation respectively. The aim is to construct equations in which $\alpha'(\eta)$ depends only on $\beta(\eta)$ and vice versa, in order to perform the multiple integral iteration of the next section.

### Bosonic case

To compute the semiclassical pair creation rate we start from (). The form of these equations is the same as in flat spacetime, the only difference being the specific time dependence of the fields. We will shortly review the well know techniques of flat spacetime (see e.g. ) applied to our field configurations. One can use the ansatz (which is inspired by a WKB expansion) $$\begin{aligned}
  q_{\textbf{k}}(\eta)= &\frac{ \alpha(\eta)}{\sqrt{\omega_{\textbf{k}}(\eta)}}e^{-\frac{{\text{i}}}{2} K_0(\eta) }+\frac{\beta(\eta)}{\sqrt{\omega_{\textbf{k}}(\eta)}}e^{\frac{{\text{i}}}{2} K_0(\eta)},\label{eq:WKB-ansatzboson}\\
 q'_{\textbf{k}}(\eta)=&-\frac{{\text{i}}\omega_{\textbf{k}}'(\eta)}{\hbar}\left[\frac{ \alpha(\eta)}{\sqrt{\omega_{\textbf{k}}(\eta)}}{\text{e}}^{-\frac{{\text{i}}}{2}K_0(\eta)}-\frac{\beta(\eta)}{\sqrt{\omega_{\textbf{k}}(\eta)}}{\text{e}}^{\frac{{\text{i}}}{2}K_0(\eta)}\right]\label{eq:WKB-ansatz2boson},
\end{aligned}$$ where $$\begin{aligned}
K_0(\eta)=\frac{2}{\hbar} \int_{-\infty}^{\eta} \omega_{\textbf{k}}(\tau) d \tau, \label{eq:K}
\end{aligned}$$ where $\alpha(\eta)$, $\beta(\eta)$ are the mode functions.
The momentum spectrum of the pair creation rate is defined as $$\begin{aligned}
n_{\textbf{k}} \equiv \lim_{\eta\rightarrow 0} \left|\beta(\eta)\right|^2, \label{eq:trans}
\end{aligned}$$ with the boundary conditions $$\begin{aligned}
 \beta(-\infty)=0, && \alpha(-\infty)=1. \label{eq:boundarycond}
\end{aligned}$$ Using this, together with the ansatz ()-() in the Klein-Gordon equation () we find that the mode functions are connected through coupled differential equations $$\begin{aligned}
 \alpha'(\eta)&=\frac{\omega_{\textbf{k}}'(\eta)}{2\omega_{\textbf{k}}(\eta)}{\text{e}}^{ i K_0(\eta)}\beta(\eta),\label{eq:alphaboson}\\
 \beta'(\eta)&=\frac{\omega_{\textbf{k}}'(\eta)}{2\omega_{\textbf{k}}(\eta)}{\text{e}}^{-i K_0(\eta)}\alpha(\eta)\label{eq:betaboson}.
\end{aligned}$$

### Fermionic case

To derive analogous equations in the fermionic case we start from the equations () and (). The initial idea of this work comes from the similarities between this equations and Eq. (11-12) of . For the semiclassical treatment one makes the following ansatz $$\begin{aligned}
\psi^{\pm}_{1}(\eta)&=\frac{C_\pm}{\sqrt{\omega_{\textbf{k}}(\eta)}}\frac{\sqrt{p(\eta)}}{\sqrt{p_z(\eta)-{\text{i}}k_\perp}}\left(\alpha^\pm(\eta)[\omega_{\textbf{k}}(\eta)-m a(\eta)]{{\text{e}}^{-\frac{{\text{i}}}{2}K(\eta)}} +\beta^\pm(\eta)[\omega_{\textbf{k}}(\eta)+m a(\eta)]{{\text{e}}^{\frac{{\text{i}}}{2}K(\eta)}}\right),\label{eq:ansatz1} \\
 \psi^{\pm}_{2}(\eta)&=\frac{\mp C_\pm}{\sqrt{\omega_{\textbf{k}}(\eta)}}\frac{\sqrt{p(\eta)}}{\sqrt{p_z(\eta)+{\text{i}}k_\perp}}\left(\alpha^\pm(\eta)[\omega_{\textbf{k}}(\eta)+m a(\eta)]{{\text{e}}^{-\frac{{\text{i}}}{2}K(\eta)}} +\beta^\pm(\eta)[\omega_{\textbf{k}}(\eta)-m a(\eta)]{{\text{e}}^{\frac{{\text{i}}}{2}K(\eta)}}\right), \label{eq:ansatz2}
\end{aligned}$$ with the integrals $$\begin{aligned}
 K(\eta)& \equiv K_0(\eta)+K_1(\eta),\label{eq:K_s}\\
%K_0(\eta)& \equiv \frac{2}{\hbar}\int^{\eta} \omega_{\textbf{k}}(\tau) d\tau, \label{eq:Kint}\\
 K_1(\eta)& \equiv k_\perp\int_{-\infty}^{\eta} \frac{m a(\tau)p_z'(\tau)}{\omega_{\textbf{k}}(\tau)p(\tau)^2}d\tau \label{eq:K_xy}.

\end{aligned}$$ Using the ansatz ()-() in () and (), we find that the mode functions are connected through coupled differential equations $$\begin{aligned}
 \alpha'^\pm(\eta)&=-\frac{\omega_{\textbf{k}}'(\eta)}{2 \omega_{\textbf{k}}(\eta)}G_{\alpha}(\eta) {\text{e}}^{{\text{i}}K(\eta)}\beta^\pm(\eta),\label{eq:alphaferm}\\
 \beta'^\pm(\eta)&=\frac{\omega_{\textbf{k}}'(\eta)}{2 \omega_{\textbf{k}}(\eta)}G_{\beta}(\eta) {\text{e}}^{-{\text{i}}K(\eta)}\alpha^\pm(\eta).\label{eq:betaferm}
\end{aligned}$$ with $$\begin{aligned}
& G_{\alpha}=\frac{ma(\eta)}{p(\eta)}-\frac{\omega_{\textbf{k}}(\eta)ma'(\eta)-{\text{i}}k_\perp p_z'(\eta)}{p(\eta)\omega_{\textbf{k}}'(\eta)}, \\
&G_{\beta}=\frac{ma(\eta)}{p(\eta)}-\frac{\omega_{\textbf{k}}(\eta)ma'(\eta)+{\text{i}}k_\perp p_z'(\eta)}{p(\eta)\omega_{\textbf{k}}'(\eta)},
\end{aligned}$$ representing the fermionic corrections to the analog bosonic case.

## Multiple integral iteration

In this section we will perform the multiple integral iteration for the fermionic case. Because of the similar form of the equations ()-() and ()-() the results for the scalar case can be derived from the ones obtained below by setting $G_\alpha(\eta)=G_\beta(\eta)=1$ and $K(\eta)\rightarrow K_0(\eta)$.
By iteratively using Eqs. () and () and the boundary conditions () one finds $$\begin{split}
& \beta^{\pm}(0)=\sum_{m=0}^\infty \int_{-\infty}^\infty d\eta_0\frac{\omega_k'(\eta_0)}{2\omega_k(\eta_0)}G_{\beta}(\eta_0){\text{e}}^{-iK(\eta_0)} \\
& \times \prod_{n=1}^m \int_{-\infty}^{\eta_{n-1}} d\tau_n\frac{\omega_k'(\tau_n)}{2\omega_k(\tau_n)}G_{\alpha}(\tau_n){\text{e}}^{iK(\tau_n)} \int_{-\infty}^{\tau_n}d\eta_n\frac{\omega_k'(\eta_n)}{2\omega_k(\eta_n)} G_{\beta}(\eta_n){\text{e}}^{-iK(\eta_n)}. \\\label{eq:multiint}
 \end{split}$$ As described in these integrals are dominated by the classical turning points given by $$\begin{aligned}
\label{TP}
 \omega_k(\eta_p^\pm)=0.
\end{aligned}$$ It can be shown (see [^1]) that for one pair of simple turning points the momentum spectrum of the pair creation rate () in a semiclassical saddlepoint approximation is equal to $$\begin{aligned}
 n^{\text{fermion}}_{\textbf{k}}=\left|{\text{e}}^{-iK(\eta_p^{-})}\right|^2 \label{eq:MomentumSpectrumConst}.
\end{aligned}$$ For the bosonic case we can find $$\begin{aligned}
 n_{\textbf{k}}^{\text{scalar}}=\left|{\text{e}}^{-i K_0(\eta_p^{-})}\right|^2 .
\end{aligned}$$ by using the substitution discussed above. We thus find that analogous to the case of two-component fields in flat space-time the difference between fermions and bosons is a factor of the form $\exp(K_1(\eta_p^-))$.

# Pair creation rate of a constant electric field in dS spacetime

In this section we derive the pair creation rate for a constant electric field, which is described in $\text{dS}_4$ by $$A(\eta)=\frac{E}{H^2 \eta}.$$ This is due to the fact that a co-moving observer, with a four-velocity $u^{\mu}$, would measure an electric field of $$E_{\mu} = u^{\nu} F_{\mu \nu}= a E\delta^{z}_{\mu},$$ which leads to a constant field strength since $E_{\mu}E^{\mu}=E^2$.
For convenience we introduce $$\mu^2 = \frac{e^2 E ^2}{H^4}+ \frac{m^2}{H^2}= \lambda^2 + \gamma^2,$$ where $\lambda= \frac{eE}{H^2}$ represent the electric field divided by the Hubble rate as usually in cosmological spacetime and $\gamma= \frac{m}{H}$ is the mass term divided by the Hubble rate. These two quantities represent the electric and gravitational contribution respectively. We can now check in which limit the semiclassical approximation holds. It requires the rate of change of the background to be small in the asymptotic future, that is $$\begin{aligned}
\left(\frac{\omega'_{\textbf{k}}(\eta)}{\omega_{\textbf{k}}^2(\eta)}\right)^2 \underset{\eta \rightarrow 0}{\sim} \mu^{-2} &&\text{ and }&& \left(\frac{\omega''_{\textbf{k}}(\eta)}{\omega_{\textbf{k}}^3(\eta)}\right) \underset{\eta \rightarrow 0}{\sim} 2\mu^{-2}
\end{aligned}$$ being small, we see that this is the case when $\mu \gg 1$. Comparing to the bosonic case we see that in the limit $\mu \gg 1$ the "-2" in () is negligible, so that the fermionic and bosonic pulsation are the same. To compute the pair creation rate we first have to compute the integrals $K_0(\eta)$ and $K_1(\eta)$ defined in () and () respectively. The value of $\eta$ at the turning point () is found to be $$\eta_p^{-}=\frac{-\lambda\frac{k_z}{k}-{\text{i}}\sqrt{\gamma^2+\lambda^2\left(1-\left(\frac{k_z}{k}\right)^2\right)}}{k}.\label{eq:TP}$$ The imaginary parts of $K_0(\eta_p^-)$ and $K_1(\eta_p^-)$ are the only ones contributing to (). One can show that $$\begin{aligned}
&{\text{Im}}[K_0(\eta_p^-)]=-\pi\left(\mu-\frac{k_z}{k} \lambda\right)\theta(k_z \lambda),\\
& {\text{Im}}[K_1(\eta_p^-)]=0,
\end{aligned}$$ Where $\theta(x)$ is the Heaviside step function. It was introduced since the real part of the turning point () has to be negative for the turning point to be inside of the closed contour which is needed in the approximation of (). We thus find that only pairs with a momentum in $z$-direction which has the same sign as $\lambda$ are produced in the semiclassical limit. This is what was called pair production in "screening" direction or "downward" tunneling in . The case of opposite sign, called "upward" tunneling is suppressed in the semiclassical limit.
$K_1(\eta_p^-)$ was the only difference between bosons and fermions and is not contributing to the number of pairs created. Thus we find that the number of pairs in the semiclassical limit for both bosons and fermions is given by $$\label{main}
%\boxed{
n_{\textbf{k}} = \exp\left[-2\pi \left(\mu-\frac{k_z}{k} \lambda\right)\right] \theta(k_z \lambda).
%}$$ We will close this section by performing the flat spacetime limit and making the relation with the bosonic case of explicit. The definition of the pair production rate is $$\Gamma \equiv \frac{1}{(2 \pi)^3 V} \int d^3 \textbf{k} \,n_{\textbf{k}},$$
where $V= a(\eta)^4 d\eta$ is the unit four volume of the spacetime. As in , an estimate for the moment when most of the particles are created can be found by analyzing when the adiabaticity is violated. This gives $$\label{estimate}
\eta \sim -\frac{\mu}{k} .$$ Using this the $k$-integral can be changed into a time integral. Going to spherical coordinates the $k_z$ integral can be performed. Putting everything together, one finds $$\Gamma= \frac{H^4}{(2\pi)^3} \frac{\mu^3}{|\lambda|} \left({\text{e}}^{2 \pi |\lambda|}-1\right) {\text{e}}^{-2 \pi \mu}.$$ We can compare this with Eq. (2.33) of , that we reproduce here $$\Gamma = \frac{H^4}{(2 \pi)^3} \frac{(|\mu|^2+1/4)^{3/2}}{\sinh(2 |\mu| \pi)} \left\{\frac{H^2}{eE} \sinh\left(\frac{2 \pi eE}{H^2} \right) + 2 \pi e^{-2 |\mu| \pi} \right\}.$$ We see that in the semiclassical limit $\mu \gg 1$ and $\lambda \gg 1$ the expressions are equal. As in the bosonic case the physical number density $n$ of produced pairs at the time $\eta$ is given by $$n=\frac{1}{a(\eta)^3} \int_{-\infty}^{\eta} d\tau a(\tau)^4 \Gamma= \frac{\Gamma}{3H}.$$ The fact that it is constant shows that the dilution from the expansion of the universe is exactly compensated by the particles created from Schwinger and gravitational particle creation. Hence in the semiclassical limit, the population of fermions is always dominated by the particles created within a Hubble time. The vacuum decay rate is defined for fermions as $$\Upsilon=\log(1-|\beta_{\textbf{k}}|^2).$$ The limit in which the Hubble parameter is negligible compared to the gravitational and electrical strength corresponds to the limit to flat spacetime. After some calculations which are analogous to the ones in one can find the Minkowski limit by taking $H\rightarrow0$, which gives $$\label{M4}
\begin{split}
& \lim_{H\rightarrow 0} \Gamma = \frac{(eE)^2}{(2 \pi)^3} \exp \left(-\frac{\pi m^2}{|eE|} \right), \\
& \lim_{H\rightarrow 0} \Upsilon =\sum_{i=1}^{\infty} \frac{1}{i^2} \frac{(eE)^2}{(2 \pi)^3} \exp \left(-\frac{i\pi m^2}{|eE|} \right).
\end{split}$$ These are the familiar results for Schwinger pair production in Minkowski spacetime .

# Conclusions

In this work, we investigated the fermionic pair creation rate by the combination of an electric and a gravitational field in $\text{dS}_4$. We first presented the basic equations of our setup and derived the corresponding Dirac equation in section . In section , we proposed a semiclassical approximation to compute the number of pairs produced. This approximation is based on a saddle point approximation of the integrals in (). To make the comparison to the analog bosonic case easier, we first reviewed its computation within our formalism. Then we presented the computation of the number of fermions created for a constant electric field. It is the main result of this paper shown in (). The pair creation rate, in this limit, is the same in the fermionic case as in the bosonic case. This equivalence between fermions and bosons in the semiclassical limit occurs also for one-component fields in flat spacetime when there is only one pair of turning points (see e.g. ). The limit to flat Minkowski spacetime is presented in () and agrees with the usual expression for the Schwinger effect.
Going beyond the semiclassical limit is possible but is out of the scope of this paper. To do this one can compute a more general quantity: the fermionic induced current. The current is proportional to the number of pairs created in the semiclassical limit but allows to explore a regime where the notion of adiabatic vacuum and of particles does not necessarily exist. It is a more precise way of describing the Schwinger effect because its definition does not depend on the time of creation of the pairs and the rough estimate () can be avoided. A calculation of the current induced by an electric field in $\text{dS}_2$ is performed in .
Throughout this paper the electric and gravitational fields have been assumed to be external sources. One of the next steps could be to consider backreaction effects of the newly created particles to the external electric and gravitational fields. To do so, for the electromagnetic side, the electric current needs to be computed and be plugged into the equivalent of the Maxwell equation in curved spacetime. For the gravitational side, the density of particle needs to be plugged in the Friedman equation. Another direction could be to look for cosmological application *e.g*. magnetogenesis, baryogenesis and the relation between dark energy and dark matter in the context of neutrino physics.

# Acknowledgements

The authors thank Fernanda Gomes De Oliveira, Hendrik Ludwig, Pereira Jonas, Carlos Argüelles and She Sheng Xue for fruitful discussions. They are also grateful to all the people who were involved with the organization of the 2nd César Lattes Meeting in Rio de Janeiro. CS and ES are supported by the Erasmus Mundus Joint Doctorate Program by Grant Number 2012-1710 and 2013-1471 from the EACEA of the European Commission respectively.

[^1]: The detailed intermediate steps can be found in Eqs. (32)-(38) of . Observe that in the current paper the integration contour is closed in the lower imaginary half plane because of opposite convention for the phases in ().