\documentclass[aps,prb,reprint,noshowkeys,superscriptaddress]{revtex4-1}
\usepackage{graphicx,dcolumn,bm,xcolor,microtype,multirow,amscd,amsmath,amssymb,amsfonts,physics,mhchem,pifont,wrapfig}

\usepackage[utf8]{inputenc}
\usepackage[T1]{fontenc}
\usepackage{txfonts}

\usepackage[normalem]{ulem}
\definecolor{darkgreen}{RGB}{0, 180, 0}
\newcommand{\titou}[1]{\textcolor{red}{#1}}
\newcommand{\trashPFL}[1]{\textcolor{red}{\sout{#1}}}
\newcommand{\PFL}[1]{\titou{(\underline{\bf PFL}: #1)}}

\usepackage[
	colorlinks=true,
    citecolor=blue,
    breaklinks=true
	]{hyperref}
\urlstyle{same}

\newcommand{\ie}{\textit{i.e.}}
\newcommand{\eg}{\textit{e.g.}}
\newcommand{\ctab}[1]{\multicolumn{1}{c}{#1}}
\newcommand{\mc}{\multicolumn}
\newcommand{\fnm}{\footnotemark}
\newcommand{\fnt}{\footnotetext}
\newcommand{\mr}{\multirow}

% addresses
\newcommand{\LCPQ}{Laboratoire de Chimie et Physique Quantiques (UMR 5626), Universit\'e de Toulouse, CNRS, UPS, France}
\newcommand{\VRIJE}{Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modelling, FEW, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands}

\begin{document}	

\title{Transient Uniform Electron Gases}

\author{Pierre-Fran\c{c}ois \surname{Loos}}
\email{loos@irsamc.ups-tlse.fr}
\affiliation{\LCPQ}
\author{Paola \surname{Gori-Giorgi}}
\affiliation{\VRIJE}
\author{Michael \surname{Seidl}}
\affiliation{\VRIJE}

\begin{abstract}
%\begin{wrapfigure}[11]{o}[-1.2cm]{0.4\linewidth}
%	\centering
%	\includegraphics[width=\linewidth]{TOC}
%\end{wrapfigure}
The uniform electron gas (UEG), a hypothetical system with finite homogenous electron density composed by an infinite number of electrons in a box of infinite volume, is the practical pillar of density-functional theory (DFT) and the foundation of the most acclaimed approximation of DFT, the local-density approximation (LDA). 
In the last thirty years, the knowledge of analytical parametrizations of the infinite UEG (IUEG) exchange-correlation energy has allowed researchers to perform millions of approximate electronic structure calculations for atoms, molecules, and solids.
Recently, it has been shown that the traditional concept of IUEG is not the unique example of UEGs, and systems, in their lowest-energy state, consisting of electrons that are confined to the surface of a sphere provide a beautiful new family of UEGs with more customizable properties.
Here, we show that, some of the excited states associated with these systems can be classified as transient UEGs (TUEGs) as their electron density is only homogenous for very specific values of the radius of the sphere.
\end{abstract}

\maketitle

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\textbf{Uniform electron gases.}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Alongside the two Hohenberg-Kohn theorems \cite{Hohenberg_1964} which put density-functional theory (DFT) on firm mathematical grounds and the Kohn-Sham (KS) formalism \cite{Kohn_1965} that makes DFT practically feasible, the uniform electron gas (UEG) \cite{Loos_2016} is one of the many pieces of the puzzle that have made DFT \cite{ParrBook} so successful in the past thirty years.
Indeed, apart from very few exceptions, most density-functional approximations are based, at some level at least, on the UEG via the so-called local-density approximation (LDA) \cite{Thomas_1927,Fermi_1927,Dirac_1930,Slater_1951,Ceperley_1980} which assumes that the electron density $\rho$ of an atom, a molecule, or a solid is locally uniform and has identical ``properties'' to the UEG with the same electron density.

Thanks to the construction of exchange-correlation LDA functionals \cite{Slater_1951,Vosko_1980,Perdew_1981,Perdew_1992,Chachiyo_2016} which can be loosely seen as a one-to-one mapping between a given value of the electron density and the exchange-correlation energy of the UEG, one can then straightforwardly compute, within KS-DFT, the electronic ground-state energy and properties of any molecules or materials with, nonetheless, a certain degree of approximation inherently associated with the approximate nature of the exchange-correlation LDA functional.
One can also access excited states via the time-dependent version of DFT. \cite{Runge_1984,Casida_1995,Petersilka_1996,UllrichBook}

As commonly done, the LDA can be refined by adding up new ingredients, such as the gradient of the density $\nabla \rho$ [which defines the generalized gradient approximation (GGA)], \cite{Perdew_1986,Becke_1988,Perdew_1996} the kinetic energy density $\tau$ (meta-GGA), \cite{Becke_1988b,Sun_2015} exact Hartree-Fock exchange (yielding the so-called hybrid functionals), \cite{Becke_1993,Adamo_1999} and others. 
Each of these quantities defines a new rung of the well-known Jacob ladder of DFT \cite{Perdew_2001} that is supposed to bring electronic structure theory calculations from the evil Hartree world to the chemical accuracy heaven.

The UEG, also known as jellium in some context, \cite{Loos_2016} is a hypothetical infinite substance where an infinite number of electrons ``bathe'' in a (uniform) positively charged jelly of infinite volume. \cite{ParrBook,Loos_2016}
It is can be ``created`` via a \textit{gedanken} experience by pouring electrons in an expandable box while keeping the ratio $\rho = n/V$ of the number of electrons $n$ and the volume of the box $V$ constant.
In the so-called thermodynamic limit where both $n$ and $V$ goes to infinity but $\rho$ remains finite, the electron density eventually becomes homogeneous.
In the following, this paradigm is named the infinite UEG (IUEG).
 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\textbf{Finite uniform electron gases.}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Recently, it has been shown that one can create finite UEGs (FUEGs) by placing a finite number of electrons confined to the surface of a sphere of radius $R$. \cite{Tempere_2002,Tempere_2007,Seidl_2007,Loos_2009a,Loos_2009c,Loos_2010e,Loos_2011b,Gill_2012,Loos_2018b} 
Of course, FUEGs only appear for well-defined electron numbers and electronic states. \cite{Rogers_2016,Rogers_2017}
In particular, the spin-unpolarized ground state of $n$ electrons on a sphere has a homogeneous density for $n = 2(L+1)^2$ (where $L \in \mathbb{N}$) for any $R$ values, and this holds also within the Hartree-Fock approximation. \cite{Loos_2011b}
This property comes from the addition theorem of the spherical harmonics \cite{NISTbook} $Y_{\ell m}(\bm{\Omega})$ (which are the spatial orbitals of the system in this particular case):
\begin{equation}
	\sum_{\ell=0}^L \sum_{m=-\ell}^{+\ell} Y_{\ell m}^*(\bm{\Omega}) Y_{\ell m}(\bm{\Omega}) = \frac{(L+1)^2}{2\pi^2}
\end{equation}
where $\bm{\Omega}=(\theta,\phi)$ gathers the polar and azimuthal angles, respectively.
Thanks to this amazing property, these FUEGs have be employed to construct alternative LDA functionals for both KS-DFT \cite{Loos_2014b,Loos_2017a} and ensemble DFT. \cite{Loos_2020g,Marut_2020}
Besides, hints of the equivalence of the FUEG and IUEG models have been found in the thermodynamic limit, \ie, when $L \to \infty$. \cite{Bowick_2002,Loos_2011b}

The two-electron case (\ie, $L = 0$) is of particular interest \cite{Seidl_2007,Loos_2009a,Loos_2018b} as it has been shown to be extremely useful for testing electronic structure methods \cite{Seidl_2007,Loos_2009a,Pedersen_2010,Loos_2012c,Schindlmayr_2013,Sun_2016,Loos_2018b} and is, furthermore, exactly solvable for a countably infinite set of $R$ values. \cite{Loos_2009c,Loos_2010e,Loos_2012}
In this case, the many-body Hamiltonian is simply
\begin{equation}
	\label{eq:H}
	\Hat{H} = - \frac{\nabla_1^2 + \nabla_2^2}{2} + \frac{1}{r_{12}}
\end{equation}
where $r_{12} = \abs{\bm{r}_1 - \bm{r}_2}$ is the interelectronic distance, \ie, the electrons interact Coulombically through the sphere.
Following Breit, \cite{Breit_1930} one can write the total electronic wave function as
\begin{equation}
\label{eq:Phi}
	\Phi (\bm{x}_1,\bm{x}_2)
	= \Xi(s_1,s_2) \chi (\bm{r}_1,\bm{r}_2) \Psi(r_{12})
\end{equation}
where $\Xi$, $\chi$ and $\Psi$ are the spin, angular and interelectronic wave functions, respectively, and $\bm{x}_i = (s_i,\bm{r}_i)$ is a composite coordinate gathering the spin coordinate $s_i$ and the spatial coordinate $\bm{r}_i$ associated with the $i$th electron.
The singlet and triplet wave functions read \cite{BetheBook}
\begin{subequations}
\begin{align}
	^1\Xi(s_1,s_2) & 
	= \frac{1}{\sqrt{2}} \qty[ \alpha(s_1) \beta(s_2) - \beta(s_1) \alpha(s_2) ]		\\
	^3\Xi(s_1,s_2) & = 
	\begin{cases}
		\alpha(s_1) \alpha(s_2),								\\
		\frac{1}{\sqrt{2}} \qty[ \alpha(s_1) \beta(s_2) + \beta(s_1) \alpha(s_2) ]	\\
		\beta(s_1) \beta(s_2).
	\end{cases}
\end{align}
\end{subequations}

After integration over the spin coordinates, the total electronic density is, by definition,  \cite{DavidsonBook} 
\begin{equation}
	\label{eq:rho}
	\rho(\bm{r}_1) = 2 \int \chi (\bm{r}_1,\bm{r}_2)^2 \Psi(r_{12})^2 d\bm{r}_2 
\end{equation}
and, following Seidl, \cite{Seidl_2007} one can decompose the square of the (totally symmetric) interelectronic wave function over the basis set of the Legendre polynomials 
\begin{equation}
	\Psi(r_{12})^2 = \sum_{\ell=0}^{\infty} c_{\ell} P_{\ell}(\cos \gamma) 
\end{equation}
where
\begin{equation}
	c_{\ell} = \frac{2\ell+1}{2} \int_0^\pi P_{\ell}(\cos \gamma) \Psi(r_{12})^2 \sin \gamma d\gamma
\end{equation}
and the interelectronic angle $\gamma$ is defined as $r_{12} = R \sqrt{2 - 2 \cos \gamma}$ or, equivalently, $\cos \gamma = \cos \theta_1 \cos \theta_2 + \sin \theta_1 \sin \theta_2 \cos(\phi_1 - \phi_2)$.
For the singlet ground state, we have $\chi_{^1S} (\bm{r}_1,\bm{r}_2) = 1$, and taking advantage of the addition theorem of the Legendre polynomials \cite{NISTbook}
\begin{equation}
	P_{\ell}(\cos \gamma) = \frac{4\pi}{2\ell+1} \sum_{m=-\ell}^{+\ell} Y_{\ell m}^*(\bm{\Omega}_1) Y_{\ell m}(\bm{\Omega}_2)
\end{equation}
and the orthonormality of the spherical harmonics, one finds that $c_\ell = 0$ for $\ell > 0$, demonstrating that the ground state of two electrons on a surface of a sphere has a uniform density $\rho_{^1S}(\bm{r}) = 2/(4\pi R^2)$ for any $R$ values.
This result is also true for any excited states with $^1S$ symmetry. 
For the other electronic states corresponding to higher total angular momentum, such as the lowest singlet and triplet $P$ states, \cite{Loos_2010e} the electron density is, in general, inhomogeneous, except in the very unlikely conditions discussed below. \cite{Seidl_2007}

%%% FIG 1 %%%
\begin{figure}
	\includegraphics[width=\linewidth]{3P}
	\caption{
	$c_0 - 2 c_1/3 + c_2/5$ as a function of the radius of the sphere $R$ for various states of $^3P$ symmetry.
	The root associated to each state is located by a marker.
	\label{fig:3P}}
\end{figure}
%%% %%% %%% %%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\textbf{Transient uniform electron gases.}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
As evidenced by Eq.~\eqref{eq:rho}, both the angular and interelectronic wave functions contribute to the electron density and a subtle interplay between these two quantities may result in a uniform density.
To illustrate this, let us consider explicit examples.

For the $^3P$ and $^1P$ two-electron states, the angular wave functions are \cite{Breit_1930,Seidl_2007,Loos_2009a,Loos_2009c,Loos_2010e}
\begin{subequations}
\begin{align}
	\label{eq:spin_3P}
	\chi_{^3P} (\bm{r}_1,\bm{r}_2) & = \cos \theta_1 - \cos \theta_2
	\\
	\label{eq:spin_1P}
	\chi_{^1P} (\bm{r}_1,\bm{r}_2) & = \cos \theta_1 + \cos \theta_2
\end{align}
\end{subequations}
After a lengthy derivation, the square of $\chi$ can be expanded in terms of spherical harmonics, and, after integration over the coordinates of the second electron, it yields
\begin{subequations}
\begin{align}
	\label{eq:3P}
	\rho_{^3P}(\bm{r}) & = 
	\qty(c_0 - \frac{c_2}{5}) Y_{00}(\bm{\Omega})^2 
	+ \qty(c_0 - \frac{2 c_1}{3} + \frac{c_2}{5}) Y_{10}(\bm{\Omega})^2
	\\
	\label{eq:1P}
	\rho_{^1P}(\bm{r}) & = 
	\qty(c_0 - \frac{c_2}{5}) Y_{00}(\bm{\Omega})^2 
	+ \qty(c_0 + \frac{2 c_1}{3} + \frac{c_2}{5}) Y_{10}(\bm{\Omega})^2
\end{align}
\end{subequations}
As readily seen from these expressions and remembering that $Y_{00}(\bm{\Omega})^2 = 1/(4\pi)$, the electron density is uniform if and only if the component associated with $Y_{10}(\bm{\Omega})^2$ vanishes, \ie,
\begin{equation}
	\label{eq:condition}
	c_0 \mp \frac{2 c_1}{3} + \frac{c_2}{5} = 0
\end{equation}

For each (ground or excited) $^3P$ states, there exists one and only one value of the radius, $R_\text{UEG}$, for which Eq.~\eqref{eq:condition} is fulfilled, and it can be computed numerically with great precision thanks to explicitly correlated calculations. \cite{Loos_2009a}
Figure \ref{fig:3P} shows the behavior of $c_0 - 2 c_1/3 + c_2/5$ as a function of $R$ for the $^3P$ ground state and its first and second excited states.
As one can see, $R_\text{UEG}$ is negative (which corresponds to attractive ``electron'' pair \cite{Seidl_2007,Seidl_2010} or exciton \cite{Pedersen_2010,Loos_2012c}) and gets larger (in absolute value) for excited states.
Because this model has an inhomogeneous electron density except for a unique $R$ value, we name these transient UEGs (TUEGs).
Note that this feature was first discovered by Seidl in Appendix A of Ref.~\onlinecite{Seidl_2007}. 
He also provided an estimate $R_\text{UEG} \approx -5.3$ for the $^3P$ ground state, a rather good estimate that we refine here to $R_\text{UEG} \approx -5.32527$.

For the $^1P$ states, the condition $c_0 + 2 c_1/3 + c_2/5 = 0$ [see Eq.~\eqref{eq:condition}] cannot be fulfilled and, hence, these states never exhibit uniform densities.
This further highlights the subtle balance that must be accomplished between the spin and spatial parts of the wave function [see Eq.~\eqref{eq:rho}] and can help us rationalizing why the $^1P$ states are TUEGs.
The $^3P$ spin wave function, $\chi_{^3P}$, defined in Eq.~\eqref{eq:spin_3P} has the natural tendency to pull apart the same-spin electron pair in accordance with the Pauli exclusion principle. 
The same physical effect can be obtained by increasing the value of $R$ (\ie, $R \gg 0$).
In such a case, the two electrons localize (or ``crystallize'') on opposite side of the sphere to form a Wigner crystal. \cite{Wigner_1934}
Oppositely, when $R \ll 0$, the two electrons are attracted to each other to form a pair of tightly bound electrons that freely move on the sphere. \cite{Seidl_2007,Seidl_2010}
For certain $R$ values, the attractive effect stemming from the spatial part of the wave function exactly compensates the Pauli exclusion principle originating from the spin part of the wave function to make the total electron density uniform, hence producing TUEGs.
In higher-energy excited states, the same-spin electrons are further away as compared to the ground state due to the larger number of nodes in the excited-state wave functions.
Therefore, the magnitude of the attractive effect has to be larger to compensate it, which corresponds to more negative values of $R$.


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\textbf{Concluding remarks.}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Here, we have introduced the concept of transient UEGs (TUEGs), a novel family of electron gases that exhibit, in very particular conditions, homogenous densities.
Using the electrons-on-a-sphere model, we have presented an example of such TUEGs created thanks to the competing effects of the Pauli exclusion principle and the creation of an attractive electron pair. 
TUEGs with larger number of electrons certainly exists and we hope to investigate these in the future.
The present concept might be useful in the future development of exchange-correlation functionals within DFT.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\textbf{Acknowledgements.}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (Grant agreement No.~863481).
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\bibliography{TrUEGs}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\end{document}