In a seminal paper, Page found the exact formula for the average entanglement entropy for a pure random state. We consider the analogous problem for the ensemble of pure fermionic Gaussian states, which plays a crucial role in the context of random free Hamiltonians. Using recent results from random matrix theory, we show that the average entanglement entropy of pure random fermionic Gaussian states in a subsystem of $N_A$ out of $N$ degrees of freedom is given by ${\langle S_A\rangle }_{\mathrm{G}}=(N-\frac{1}{2})\mathrm{\Psi }\left(2N\right)+(\frac{1}{4}-N_A)\mathrm{\Psi }\left(N\right)+(\frac{1}{2}+N_A-N)\mathrm{\Psi }(2N-2N_A)-\frac{1}{4}\mathrm{\Psi }(N-N_A)-N_A$, where $\mathrm{\Psi }$ is the digamma function. Its asymptotic behavior in the thermodynamic limit is given by ${\langle S_A\rangle }_{\mathrm{G}}=N(log2-1)f+N(f-1)log(1-f)+\frac{1}{2}f+\frac{1}{4}log(1-f)+O(1/N)$, where $f=N_A/N\le 1/2$. Remarkably, its leading order agrees with the average over eigenstates of random quadratic Hamiltonians with number conservation, as found by Łydżba, Rigol, and Vidmar. Finally, we compute the variance in the thermodynamic limit, given by the constant ${lim}_{N\to \infty }{\left(\mathrm{\Delta }{S}_A\right)}_{\mathrm{G}}^2=\frac{1}{2}[f+f^2+log(1-f)]$.

• Received 19 March 2021
• Accepted 15 June 2021

DOI:https://doi.org/10.1103/PhysRevB.103.L241118