---
abstract: 'We study the Hartree-Fock equation universally used in many-electron problems. We prove that for any negative constant there are only finite number of critical values of Hartree-Fock’s functional associated with eigenvalues less than the constant. We also prove that a negative accumulation point of negative eigenvalues is an eigenvalue.'
author:
- Sohei Ashida
title: 'Finiteness of Hartree-Fock energies associated with eigenvalues of the Hartree-Fock equation less than a negative constant'
---

Introduction {#firstsec}
============

In this paper we study the Hartree-Fock equation. We consider the space $L^2(\mathbb R^3;\mathbb C^2)$ of $L^2$ functions with a spin variable. We label the two components of $\varphi\in L^2(\mathbb R^3;\mathbb C^2)$ by $\frac{1}{2}$ and $-\frac{1}{2}$. Thus the two components are denoted by $\varphi(\cdot,\frac{1}{2})$ and $\varphi(\cdot,-\frac{1}{2})$. For $x\in \mathbb R^3$ and $\sigma\in \{\frac{1}{2},-\frac{1}{2}\}$ we use the bold-faced variable as $\mathbf x=(x,\sigma)$. We write $\int f(\mathbf y)d\mathbf y:=\int f(y,\frac{1}{2})dy+\int f(y,-\frac{1}{2})dy$, $\langle f,g\rangle:=\int f^*(\mathbf y)g(\mathbf y)d\mathbf y$, and $\lVert f\rVert:=\langle f,f\rangle^{1/2}$. Fix the number of electrons $n\in \mathbb N$, number of nuclei $N\in \mathbb N$, nuclear charges $Z_k\in \mathbb N,\ k=1,\dots,n$, and the positions of the nuclei $x_k\in \mathbb R^3,\ k=1,\dots,n$. After Hartree [@Ha] introduced the Hartree equation ignoring the symmetry with respect to exchange of variables, considering the symmetry the Hartree-Fock equation was introduced by Fock [@Fo] and Slater [@Sl] to obtain critical values and the critical points of the functional $\langle \Phi,H\Phi\rangle$, where $$\label{myeq1.1}
H:=-\frac{1}{2}\sum_{j=1}^N\Delta_{y_j}-\sum_{\substack{1\leq j \leq N\\ 1\leq k\leq n}}\frac{Z_k}{\lvert y_j-x_k\rvert}+\sum_{1\leq j <k\leq N}\frac{1}{\lvert y_j-y_k\rvert},$$ is an electronic Hamiltonian and $$\label{myeq2.2}
\Phi(\mathbf x_1,\dots,\mathbf x_N):=(N!)^{-1/2}\sum_{\tau\in\mathbf S_N}(\mathrm{sgn}\, \tau)\varphi_1(\mathbf x_{\tau(1)})\cdots\varphi_N(\mathbf x_{\tau(N)}),$$ is the Slater determinant. Here $\mathbf S_N$ is the symmetric group, $\mathrm{sgn}\, \tau$ is the signature of $\tau$, and $\varphi_i\in H^2(\mathbb R^3;\mathbb C^2):=\{\varphi\in L^2(\mathbb R^3;\mathbb C^2): \lVert \Delta\varphi\rVert<\infty\}$ with constraints $\langle \varphi_i,\varphi_j\rangle=\delta_{ij}$. The Hartree-Fock equation is the corresponding Euler-Lagrange equation: $$\label{myeq2.1}
\begin{split}
&-(\Delta\varphi_i)(x,\sigma)-\sum_{1\leq k\leq n}\frac{Z_k}{\lvert x-x_k\rvert}\varphi_i(x,\sigma)\\
&\quad+\sum_{j=1}^N\left[\int\frac{1}{\lvert x-y\rvert}\lvert \varphi_j(\mathbf y)\rvert^2d\mathbf y\varphi_i(x,\sigma)-\int\frac{1}{\lvert x-y\rvert}\varphi_j^*(\mathbf y)\varphi_i(\mathbf y)d\mathbf y\varphi_j(x,\sigma)\right]\\
&=\epsilon_i\varphi_i(x,\sigma),\ i=1,\dots,N,
\end{split}$$ with constraints $\langle \varphi_i,\varphi_j\rangle=\delta_{ij}$, where $\epsilon_i\in \mathbb R,\ i=1,\dots,N$ are Lagrange multipliers. Since the Slater determinant is a sum of products of functions, the Hartree-Fock equation obtained by the method of Lagrange multiplier is a system of nonlinear partial differential equations with unknown constants $\epsilon_i$ which are Lagrange multipliers. Since an eigenvalue of a linear partial differential operator such as an electronic Hamiltonian can be regarded as a Lagrange multiplier of a variational problem of a functional, we also call $(\epsilon_1,\dots,\epsilon_N)$ an eigenvalue if has a solution $(\varphi_1,\dots,\varphi_N)$, and corresponding $(\varphi_1,\dots,\varphi_N)$ an eigenfunction associated with $(\epsilon_1,\dots,\epsilon_N)$. One of the difficulty in the analysis of the eigenvalue problem is the nonlinearity of the equation. In addition, also the number of constraints of the corresponding minimization problem is a substantial difference. In the case of the linear eigenvalue problem of a Hamiltonian, only the norm of a function is constrained. Thus the eigenvalue (Lagrange multiplier) is equal to a critical value of the functional. However, in the case of the Hartree-Fock equation many Lagrange multipliers concerned with many constraints appear in the equation, and there is no clear relation between each multiplier and the critical value of the functional. Because of these reasons we can not use the methods for linear eigenvalue problems to study the critical values of Hartree-Fock’s functional and the eigenvalues of the Hartree-Fock equation. For example, we can not see immediately if the critical values or eigenvalues are countable or continuously distributed in some interval. This problem would be important in the study of convergence of the approximation methods such as the so-called self consistent field method.

Our main result is that for any negative constant there exist at most a finite number of critical values of Hartree-Fock’s functional with each component of eigenvalues smaller than the constant, which is a progress in the problems as above. To prove this result we first show that an accumulation point of such eigenvalues is also an eigenvalue. Next, we investigate the Frechet derivative of the Hartree-Fock equation and use the theorem by Fučik-Nečas-Souček-Souček [@FNSS] which shows the set of critical values near the one corresponding to the eigenfunction is one-point set, using implicit function theorem and the property of real-analytic functions in Souček-Souček [@SS].

Another fundamental problem is the existence of the eigenvalues and eigenfunctions which we do not deal with in this paper. Lieb-Simon [@LS] proved that if the number of electrons is smaller than or equal to the total charge of the nuclei, there exists a solution to the Hartree-Fock equation which minimizes the functional. Lions [@Li] proved that, under the same assumption on the total charge and the number of electrons there exists a sequence of eigenfunctions of the Hartree-Fock equation with nonpositive eigenvalues and corresponding values of the functional are converging to $0$. Lewin [@Le] constructed infinitely many solutions corresponding to critical values of Hartree-Fock’s functional less than the first energy threshold. If we ignore the symmetry and consider the expectation value of the electronic Hamiltonian with respect to a product of functions, the critical points satisfy the Hartree equation instead of the Hartree-Fock equation. For the existence of the solutions to the Hartree equation see [@Re; @Wo; @St; @LS; @Li].

In Section 3 the main result is stated and several lemmas needed in the proof of the main result are introduced. In Section 4 we prove the main result.

Main result and some preliminaries
==================================

Main result
-----------

Let $(\epsilon_1,\dots,\epsilon_N)$ be an eigenvalue and $(\varphi_1,\dots,\varphi_N)$ be an eigenfunction. Then let us call the value of $\langle \Phi,H\Phi\rangle$ the Hartree-Fock energy associated with $(\epsilon_1,\dots,\epsilon_N)$, where $H$ and $\Phi$ is defined by and respectively. Our main result is the following.

\[Finthm\] For any $\epsilon<0$ there exist at most finite number of Hartree-Fock energies associated with eigenvalues $(\epsilon_1,\dots,\epsilon_N)\in (-\infty,\epsilon)^N$ of the Hartree-Fock equation .

As mentioned in [@LS], the naive Euler-Lagrange equation for $\langle \Phi,H\Phi\rangle$ is more complicated than the Hartee-Fock equation , but after a unitary change $\varphi_i^{\mathrm{New}}=\sum a_{ij}\varphi_j^{\mathrm{Old}},$ with a unitary $N\times N$ matrix $a_{ij}$, is satisfied by $$(\varphi_1^{\mathrm{New}},\dots,\varphi_N^{\mathrm{New}}),$$ and some $(\epsilon_1,\dots,\epsilon_N)$. By Theorem \[Finthm\] we can see that if there exist infinitely many different critical values of $\langle \Phi,H\Phi\rangle$, there exist infinitely many nonnegative $\epsilon_i$ corresponding to the critical points or $\epsilon_i$ accumulates at $0$ for some $i$. Theorem \[Finthm\] does not prohibit existence of infinitely many critical values of $\langle \Phi,H\Phi\rangle$ less than a negative constant. Actually, it is proved that there exist such critical values in [@Le].

Decay and convergence of eigenfunctions
---------------------------------------

As in [@LS] when $(\varphi_1,\dots,\varphi_N)$ is an eigenfunction of , then the vector valued function $\Psi:={}^t(\varphi_1,\dots,\varphi_N)\in \bigoplus_{j=1}^NH^2(\mathbb R^3;\mathbb C^2)$ satisfies $A\Psi=0$ where $$A_{ij}:=\delta_{ij}(-\Delta+V(x)+R(x)-\epsilon_i)+Q_{ij}(x),$$ with $$V(x):=\sum_{1\leq k<n}\frac{Z_k}{\lvert x-x_k\rvert},$$ $$\label{myeq3.0.0}
R(x):=\sum_{j=1}^N\int\lvert x-y\rvert^{-1}\varphi_j^*(\mathbf y)\varphi_j(\mathbf y)d\mathbf y,$$ and $$Q_{ij}(x):=-\int\lvert x-y\rvert^{-1}\varphi_j^*(\mathbf y)\varphi_i(\mathbf y)d\mathbf y.$$ Hereafter, we identify $(\varphi_1,\dots,\varphi_N)$ with $\Psi$, and also call $\Psi$ an eigenfunction if $(\varphi_1,\dots,\varphi_N)$ is an eigenfunction. Since $R(x)$ and $Q_{ij}(x)$ are decaying potentials, as in the scalar valued case we have the exponential decay of eigenfuncitons basically by the method of Agmon [@Ag]. Because we need uniform exponential decay for a sequence of eigenfunctions, we estimate the decay of $R(x)$ and $Q_{ij}(x)$ uniformly under assumptions on eigenfunctions. For $r\in \mathbb R$ set $B_r:=\{x\in \mathbb R^3:\lvert x\rvert<r\}$ and denote by $\lVert \varphi\rVert_{L^2(\mathbb R^3\setminus B_r)}$ the $L^2$ norm on $\mathbb R^3\setminus B_r$. It is known that $\lvert x\rvert^{-1}$ is $\Delta$-bounded (see e.g. [@Ka Chapter V. §5]). Hence, there exists a constant $C>0$ such that the following inequality holds: $$\lVert \lvert x\rvert^{-1}u\rVert\leq C(\lVert\Delta u\rVert+\lVert u\rVert).$$

\[expdecay\] Let $\epsilon<\tilde \epsilon<0$, $c,\tilde r>0$. Then there exists $\tilde c>0$ such that for any eigenfunction $(\varphi_1,\dots,\varphi_N)$ of the Hartree-Fock equation associated with $(\epsilon_1,\dots,\epsilon_N)\in(-c,\epsilon)^N$ and satisfying $\lVert\Delta\varphi_i\rVert<c$ and $\lVert\varphi_i\rVert_{L^2(\mathbb R^3\setminus B_r)}<\frac{\tilde \epsilon-\epsilon}{4NC(c+1)},\ i=1,\dots,N$ for $r>\tilde r$ the following estimate holds. $$\lVert\exp(\lvert\tilde \epsilon\rvert^{1/2}\lvert x\rvert)\varphi_i(x)\rVert\leq \tilde c,\ i=1,\dots,N.$$

Let us first estimate the decay of $Q_{ij}(x)$. Since $\lvert x\rvert^{-1}$ is $\Delta$-bounded and the center of the Coulomb potential is irrelevant to the relative bound, for any $x\in \mathbb R^3$ the following holds $$\lVert\lvert x-y\rvert^{-1}\varphi_j\rVert_{L^2(\mathbb R^3_y)}\leq C(\lVert \Delta \varphi_j\rVert+\lVert \varphi_j\rVert)\leq C(c+1),$$\[myeq3.0\] where we used $\lVert \varphi_j\rVert=1$. We divide $Q_{ij}(x)$ into two parts: $$\begin{aligned}
Q_{ij}(x)=&-\int_{\lvert y\rvert<\lvert x\rvert/2}\lvert x-y\rvert^{-1}\varphi_j^*(\mathbf y)\varphi_i(\mathbf y)d\mathbf y\\
&\quad-\int_{\lvert y\rvert\geq\lvert x\rvert/2}\lvert x-y\rvert^{-1}\varphi_j^*(\mathbf y)\varphi_i(\mathbf y)d\mathbf y.\end{aligned}$$ The first term is estimated as $$\begin{aligned}
\left\lvert\int_{\lvert y\rvert<\lvert x\rvert/2}\lvert x-y\rvert^{-1}\varphi_j^*(\mathbf y)\varphi_i(\mathbf y)d\mathbf y\right\rvert&\leq2\lvert x\rvert^{-1}\int_{\lvert y\rvert<\lvert x\rvert/2}\lvert\varphi_j^*(\mathbf y)\varphi_i(\mathbf y)\rvert d\mathbf y\\
&\leq2\lvert x\rvert^{-1}\lVert \varphi_j\rVert\lVert \varphi_i\rVert=2\lvert x\rvert^{-1}.\end{aligned}$$ The second term is estimated as $$\begin{aligned}
\left\lvert\int_{\lvert y\rvert\geq\lvert x\rvert/2}\lvert x-y\rvert^{-1}\varphi_j^*(\mathbf y)\varphi_i(\mathbf y)d\mathbf y\right\rvert&\leq\lVert\lvert x-y\rvert^{-1}\varphi_j^*(\mathbf y)\rVert_{L^2(\mathbb R^3_y)}\lVert \varphi_i\rVert_{L^2(\mathbb R^3\setminus B_{\lvert x\rvert/2})}\\
&\leq C(c+1)\lVert \varphi_i\rVert_{L^2(\mathbb R^3\setminus B_{\lvert x\rvert/2})}.\end{aligned}$$ where we used . Thus by the assumption of the lemma for sufficiently large $r$ independent of $(\varphi_1,\dots,\varphi_N)$ we obtain $$\label{myeq3.0.1}
\begin{split}
\lvert Q_{ij}(x)\rvert&\leq2\lvert x\rvert^{-1}+C(c+1)\frac{\tilde \epsilon-\epsilon}{4NC(c+1)}\\
&< C(c+1)\frac{\tilde \epsilon-\epsilon}{4NC(c+1)}=\frac{\tilde \epsilon-\epsilon}{4N},\ \lvert x\rvert>r.
\end{split}$$ In the same way, for sufficiently large $r>0$ independent of $(\varphi_1,\dots,\varphi_N)$ we can show $$\label{myeq3.0.2}
\lvert R(x)\rvert< \frac{\tilde \epsilon-\epsilon}{4},\ \lvert x\rvert>r.$$

Let $\eta(r)\in C_0^{\infty}(\mathbb R)$ be a function such that $\eta(r)=r$ for $r<1$, $\lvert\eta'(r)\rvert\leq 1$ and $\lvert\eta(r/k)\rvert\geq\lvert\eta(r)\rvert$ for any $r>0$ and $k\in \mathbb N$. Set $\rho_k(x):=\lvert \tilde\epsilon\rvert^{1/2}\eta(\langle x\rangle/k)$ and $\chi_k(x):=e^{\rho_k(x)}$, where $\langle x\rangle:=\sqrt{1+\lvert x\rvert^2}$. By a direct calculation we have $$\begin{aligned}
\mathrm{Re}\, \langle(-\Delta\varphi_i),\chi_k^2\varphi_i\rangle&=\rVert\nabla(\chi_k\varphi_i)\lVert^2-\lVert (\nabla \chi_k)\varphi_i\rVert^2\\
&=\rVert\nabla(\chi_k\varphi_i)\lVert^2-\lVert (\nabla \rho_k)\chi_k\varphi_i\rVert^2\end{aligned}$$ By $A\Psi=0$ as in the beginning of this section we have $$\label{myeq3.0.2.1}
\begin{split}
0&=\sum_{i=1}^N\bigg\{\rVert\nabla(\chi_k\varphi_i)\lVert^2-\lVert (\nabla \rho_k)\chi_k\varphi_i\rVert^2+\langle(V(x)+R(x)-\epsilon_i)\chi_k\varphi_i,\chi_k\varphi_i\rangle\\
&\qquad+\sum_{j=1}^N\langle Q_{ij}(x)\chi_k\varphi_j,\chi_k\varphi_i\rangle\bigg\}\\
&\geq\sum_{i=1}^N\bigg\{\langle(V(x)+R(x)-\epsilon_i-\lvert\nabla\rho_k\rvert^2\\
&\qquad-\sum_{j=1}^N(\lvert Q_{ij}(x)\rvert+\lvert Q_{ji}(x)\rvert)/2)\chi_k\varphi_i,\chi_k\varphi_i\rangle\bigg\}\\
&=\sum_{i=1}^N\langle(S_i(x)-\epsilon_i-\lvert\nabla\rho_k\rvert^2)\chi_k\varphi_i,\chi_k\varphi_i\rangle,
\end{split}$$ where $S_i(x):=V(x)+R(x)-\sum_{j=1}^N(\lvert Q_{ij}(x)\rvert+\lvert Q_{ji}(x)\rvert)/2$ and we used $$\lvert\langle Q_{ij}(x)\chi_k\varphi_j,\chi_k\varphi_i\rangle\rvert\leq \frac{1}{2}\{\langle\lvert Q_{ij}(x)\rvert\chi_k\varphi_i,\chi_k\varphi_i\rangle+\langle\lvert Q_{ij}(x)\rvert\chi_k\varphi_j,\chi_k\varphi_j\rangle\}.$$ Since $\lvert V(x)\rvert\to0$ as $\lvert x\rvert\to\infty$, by and there exists $r_0>0$ independent of $(\varphi_1,\dots,\varphi_N)$ such that $\lvert S_i(x)\rvert<(\tilde\epsilon-\epsilon)/2$ for $\lvert x\rvert >r_0$. Thus noting $\lvert\nabla\rho_k\rvert^2<-\tilde \epsilon$ we have $$\begin{aligned}
&\sum_{i=1}^N\int_{\lvert x\rvert>r_0}(-\epsilon_i+\tilde \epsilon-(\tilde\epsilon-\epsilon)/2)\lvert\chi_k\varphi_i(\mathbf x)\rvert^2d\mathbf x\\
&\quad\leq\sum_{i=1}^N\int_{\lvert x\rvert>r_0}(-\epsilon_i-\lvert\nabla\rho_k\rvert^2-(\tilde\epsilon-\epsilon)/2)\lvert\chi_k\varphi_i(\mathbf x)\rvert^2d\mathbf x\\
&\quad\leq\sum_{i=1}^N\int_{\lvert x\rvert\leq r_0}(S_i(x)-\epsilon_i-\lvert\nabla\rho_k\rvert^2)\lvert\chi_k\varphi_i(\mathbf x)\rvert^2d\mathbf x\\
&\quad\leq\sum_{i=1}^N\sup_{\lvert x\rvert\leq r_0}\lvert\chi_k\rvert\int_{\lvert x\rvert\leq r_0}(S_i(x)-\epsilon_i-\lvert\nabla\rho_k\rvert^2)\lvert\varphi_i(\mathbf x)\rvert^2d\mathbf x,\end{aligned}$$ where we used in the second inequality. We can estimate as $\sup_{\lvert x\rvert\leq r_0}\lvert\chi_k\rvert\leq e^{\lvert\tilde\epsilon\rvert^{1/2}\sup\eta}$. Since $V$ and $\lvert x-y\rvert^{-1}$ are $\Delta$-bounded, $\lVert\varphi_i\rVert=1$ and $\lVert\Delta\varphi_i\rVert<c$, there exists a constant $C_{r_0}>0$ independent of $(\varphi_1,\dots,\varphi_N)$ such that $$\int_{\lvert x\rvert\leq r_0}(S_i(x)-\epsilon_i-\tilde \epsilon)\lvert\varphi_i(\mathbf x)\rvert^2d\mathbf x
\leq C_{r_0}.$$ Thus there exists $\tilde C_{r_0}>0$ such that $$\begin{aligned}
\tilde C_{r_0}&\geq\sum_{i=1}^N\int_{\lvert x\rvert>r_0}(-\epsilon_i+\tilde \epsilon-(\tilde\epsilon-\epsilon)/2)\lvert\chi_k\varphi_i(\mathbf x)\rvert^2d\mathbf x\\
&\geq(\tilde\epsilon-\epsilon)/2\sum_{i=1}^N\int_{\lvert x\rvert>r_0}\lvert\chi_k\varphi_i(\mathbf x)\rvert^2d\mathbf x\\
&\to(\tilde\epsilon-\epsilon)/2\sum_{i=1}^N\int_{\lvert x\rvert>r_0}\lvert e^{\lvert\tilde \epsilon\rvert^{1/2}\langle x\rangle}\varphi_i(\mathbf x)\rvert^2d\mathbf x,\end{aligned}$$ where we used $-\epsilon_i>-\epsilon$ and the monotone convergence theorem. The result of the lemma follows immediately from this inequality.

Using this uniform exponential decay we have the following lemma.

\[convseq\] Let $\mathcal E^m:=(\epsilon_1^m,\dots,\epsilon_N^m)\in(-\infty,0)^N,\ m=1,2,\dots$ be a sequence of eigenvalues converging to $\mathcal E:=(\epsilon_1,\dots,\epsilon_N)$ and $(\varphi_1^m,\dots,\varphi_N^m)$ be the associated eigenfunctions. Then $\mathcal E$ is an eigenvalue and there exists a subsequence of $(\varphi_1^m,\dots,\varphi_N^m)$ converging in $\bigoplus_{i=1}^NH^2(\mathbb R^3;\mathbb C^2)$ to an eigenfunction associated to $\mathcal E$. In particular, the set of eigenvalues in $(-\infty,0)^N$ is a closed set with respect to the topology induced by the one of $\mathbb R^N$.

Writing $$\label{myeq3.0.3}
R^m(x):=\sum_{j=1}^N\int\lvert x-y\rvert^{-1}(\varphi_j^m)^*(\mathbf y)\varphi_j^m(\mathbf y)d\mathbf y,$$ $$\label{myeq3.0.4}
(K^mw)(x,\sigma):=\sum_{j=1}^N\int\lvert x-y\rvert^{-1}(\varphi_j^m)^*(\mathbf y)w(\mathbf y)d\mathbf y\varphi_j^m(x,\sigma),$$ where $(\varphi_1^m,\dots,\varphi_N^m)$ are eigenfunctions associated with $\mathcal E^m$, the Hartree-Fock equation for $(\varphi_1^m,\dots,\varphi_N^m)$ is written as $$\label{myeq3.1}
-\Delta\varphi_i^m(x,\sigma)=(\epsilon_i^m-V(x)-R^m(x))\varphi_i^m(x,\sigma)+(K^m\varphi_i^m)(x,\sigma),\ i=1,\dots,N.$$ Here we note $V$ is $\Delta$-bounded with relative bound $0$. Thus for any $0<\gamma<1$ there exists $C>0$ such that $$\lVert Vu\rVert\leq \gamma\lVert \Delta u\rVert+C\lVert u\rVert.$$ Since the center of the Coulomb potential is irrelevant to the relative bound, for any $0<\tilde \gamma<1$ there exists $\tilde C>0$ such that for any $x\in \mathbb R^3$ $$\lVert \lvert x-y\rvert^{-1}u(y)\rVert_{L^2(\mathbb R^3_y)}\leq\tilde \gamma\lVert \Delta u\rVert_{L^2(\mathbb R^3_y)}+\tilde C\lVert  u\rVert_{L^2(\mathbb R^3_y)}.$$ Noting that $(\varphi_1^m,\dots,\varphi_N^m)$ satisfies $\lVert\varphi_j^m\rVert=1,\ j=1,\dots,N,$ we obtain $$\lVert R^m(x)\varphi_i^m(x)\rVert\leq N\tilde \gamma\lVert \Delta\varphi_i^m\rVert+N\tilde C\lVert \varphi_i^m\rVert,$$ $$\lVert K^m\varphi_i^m\rVert\leq N\tilde \gamma\lVert \Delta\varphi_i^m\rVert+N\tilde C\lVert \varphi_i^m\rVert.$$ Since $\mathcal E^m$ is convergent, there exists $\hat C>0$ such that $\lvert\epsilon_i^m\rvert<\hat C$ for any $i$ and $m$. Therefore, by we have $$\lVert\Delta\varphi_i^m\rVert\leq (\gamma+2N\tilde \gamma)\lVert\Delta \varphi_i^m\rVert+(C+2N\tilde C+\hat C)\lVert\varphi_i^m\rVert.$$ Since we can choose arbitrarily small $\gamma$ and $\tilde \gamma$, we can suppose $\gamma+2N\tilde \gamma<1$. Thus we obtain $$\label{myeq3.2}
\begin{split}
\lVert\Delta\varphi_i^m\rVert&\leq (1-\gamma-2N\tilde \gamma)^{-1}(C+2N\tilde C+\hat C)\lVert\varphi_i^m\rVert\\
&=(1-\gamma-2N\tilde \gamma)^{-1}(C+2N\tilde C+\hat C).
\end{split}$$ By and the Rellich selection theorem, for any $p\in \mathbb N$ there exists a subsequence $\varphi_i^{m_l}$ of $\varphi_i^m$ such that $$\lVert\varphi_i^{m_{l_1}}-\varphi_i^{m_{l_2}}\rVert_{L^2(B_{p})}+\lVert\nabla(\varphi_i^{m_{l_1}}-\varphi_i^{m_{l_2}})\rVert_{L^2(B_{p})}\to0,$$ as $l_1,l_2\to\infty$. Thus there exists a subsequence $\varphi_i^{m_k}$ of $\varphi_i^m$ such that $$\lVert\varphi_i^{m_{k_1}}-\varphi_i^{m_{k_2}}\rVert_{L^2(B_{k_0})}+\lVert\nabla(\varphi_i^{m_{k_1}}-\varphi_i^{m_{k_2}})\rVert_{L^2(B_{k_0})}<k_0^{-1},$$ where $k_0:=\min\{k_1,k_2\}$. Therefore, noting $\lVert\varphi_i^{m_k}\rVert=1$ for any $\epsilon>0$ there exists $\tilde k\in\mathbb N$ such that $\lVert\varphi_i^{m_k}\rVert_{L^2(\mathbb R^3\setminus B_{\tilde k})}=(1-\lVert\varphi_i^{m_k}\rVert_{L^2(B_{\tilde k})}^2)^{1/2}<\epsilon,\ \forall k>\tilde k$. Accordingly, by Lemma \[expdecay\] there exists $\tilde c>0$ such that $\lVert \langle x\rangle^2\varphi_i^{m_k}\rVert\leq\tilde c$ for any $i$ and sufficiently large $k$. Thus noting there exists $C_0,C_1>0$ such that for sufficiently large $k$, $\lVert \langle x\rangle\nabla \varphi_i^{m_k}\rVert\leq C_0(\lVert \langle x\rangle^2\varphi_i^{m_k}\rVert+\lVert \Delta \varphi_i^{m_k}\rVert)\leq C_1$, we obtain $$\begin{aligned}
&\lVert\varphi_i^{m_{k_1}}-\varphi_i^{m_{k_2}}\rVert+\lVert\nabla(\varphi_i^{m_{k_1}}-\varphi_i^{m_{k_2}})\rVert\\
&\leq \lVert\varphi_i^{m_{k_1}}-\varphi_i^{m_{k_2}}\rVert_{L^2(B_{k_0})}+\lVert\nabla(\varphi_i^{m_{k_1}}-\varphi_i^{m_{k_2}})\rVert_{L^2(B_{k_0})}\\
&\quad+\lVert\varphi_i^{m_{k_1}}-\varphi_i^{m_{k_2}}\rVert_{L^2(\mathbb R^3\setminus B_{k_0})}+\lVert\nabla(\varphi_i^{m_{k_1}}-\varphi_i^{m_{k_2}})\rVert_{L^2(\mathbb R^3\setminus B_{k_0})}\\
&\leq k_0^{-1}+k_0^{-1}\lVert\langle x\rangle(\varphi_i^{m_{k_1}}-\varphi_i^{m_{k_2}})\rVert_{L^2(\mathbb R^3\setminus B_{k_0})}\\
&\quad+k_0^{-1}\lVert\langle x\rangle\nabla(\varphi_i^{m_{k_1}}-\varphi_i^{m_{k_2}})\rVert_{L^2(\mathbb R^3\setminus B_{k_0})}\\
&\leq k_0^{-1}(1+2\tilde c+2C_1),\end{aligned}$$ for sufficiently large $k_0$. Thus $\varphi_i^{m_{k}}$ is a Cauchy sequence in $H^1(\mathbb R^3;\mathbb C^2)$.

Since $\lvert x-y\rvert^{-1}$ is $\Delta$-bounded, $\lvert x-y\rvert^{-1/2}$ is $\langle \nabla\rangle$-bounded. Using the Hartree-Fock equation we can see that there exists $C_2>0$ such that $$\begin{aligned}
&\lVert h(\varphi_i^{m_{k_1}}-\varphi_i^{m_{k_2}})\rVert\\
&\quad=\lVert (\epsilon_i^{m_{k_1}}-R^{m_{k_1}}(x))\varphi_i^{m_{k_1}}+K^{m_{k_1}}\varphi_i^{m_{k_1}}-(\epsilon_i^{m_{k_2}}-R^{m_{k_2}}(x))\varphi_i^{m_{k_2}}\\
&\qquad-K^{m_{k_2}}\varphi_i^{m_{k_2}}\rVert\\
&\quad\leq \sum_{i=1}^NC_2\lVert \varphi_i^{m_{k_1}}-\varphi_i^{m_{k_2}}\rVert_{H^1(\mathbb R^3)},\end{aligned}$$ where $h:=-\Delta+V$. Because $V$ is $\Delta$-bounded with relative bound smaller than $1$, $\Delta$ is $h$-bounded. Therefore, we can see that $\varphi_i^{m_{k}}$ is a Cauchy sequence in $H^2(\mathbb R^3;\mathbb C^2)$. Let $\varphi_i\in H^2(\mathbb R^3;\mathbb C^2)$ be the limit. Then the both sides of the Hartree-Fock equation converge in $L^2(\mathbb R^3;\mathbb C^2)$ and $$(h+R)\varphi_i+K\varphi_i=\epsilon_i\varphi_i,$$ where $R$ and $K$ are defined by and $$\label{myeq4.2}
(Kw)(x,\sigma):=\sum_{j=1}^N\int\lvert x-y\rvert^{-1}\varphi_j^*(\mathbf y)w(\mathbf y)d\mathbf y\varphi_j(x,\sigma),$$ respectively. Thus $\mathcal E_i$ is an eigenvalue with an eigenfunction $(\varphi_1,\dots,\varphi_N)$.

Real-analytic operators in Banach space
---------------------------------------

In this section following [@FNSS] we introduce the real-analytic operators and their property. Let $X$ and $Y$ be real Banach spaces. We denote the norm of $X$ by $\lVert \cdot\rVert$.

Let $D$ be an open subset of $X$. The mapping $F:D\to Y$ is said to be real-analytic on $D$ if the following conditions are satisfied:

-   For each $x\in D$ there exist Frechet derivatives of arbitrary orders $d^mF(x,\dots)$.

-   For each $x\in D$ there exists $\delta>0$ such that for any $h\in X$ satisfying $\lVert h\rVert<\delta$ one has $$F(x+h)=\sum_{m=0}^{\infty}\frac{1}{m!}d^mF(x,h,\dots,h).$$

\[onepoint\] Let $f$ be a real-analytic functional on a Banach space $Y_1$ and let $Y_2$ be another Banach space. Suppose that there exists a bilinear form $\langle\langle\cdot,\cdot\rangle\rangle$ on $Y_1\times Y_2$ such that for fixed $y\in Y_1$, $\langle\langle y,\cdot\rangle\rangle$ is continuous on $Y_2$ and $\langle\langle y,x\rangle\rangle=0$ for all $y\in Y_1$ implies $x=0$. For each $y\in Y_1$ suppose there exists $F(y)$ such that

-   $df(y,h)=\langle\langle h,F(y)\rangle\rangle$ for each $h\in Y_1$.

Let the operator

-   $F: Y_1\to Y_2$ is real-analytic.

Denote $\mathcal B_f:=\{y\in Y_1: f'(y)=0\}$ and let $y_0\in \mathcal B_f$. Suppose that

-   $$F'(y_0)=J+M,$$ where $J$ is an isomorphism of $Y_1$ onto $Y_2$ and $M$ is a compact operator.

Then there exists a neighborhood $U(y_0)$ in $Y_1$ of a point $y_0$ such that $f(\mathcal B_f\cap U(y_0))$ is a one-point set.

Proof of Theorem \[Finthm\]
===========================

In this section we prove Theorem \[Finthm\]. The method is based on analysis of Frechet derivative of Hartree-Fock equation at the limit point of eigenfunctions.

We assume there are infinitely many Hartree-Fock energies $\langle \Phi^m,H\Phi^m\rangle$ associated with eigenvalues $\mathcal E^m=(\epsilon^m_1,\dots,\epsilon^m_N)\in(-\infty,\epsilon)^N$ and eigenfunctions $\Psi^m={}^t(\varphi^m_1,\dots,\varphi^m_N)$ and satisfying $\langle \Phi^{m_1},H\Phi^{m_1}\rangle\neq \langle \Phi^{m_2},H\Phi^{m_2}\rangle,\ m_1\neq m_2$, where $\Phi^m$ is the Slater determinant of $\Psi^m$. We shall show this assumption leads to a contradiction.

*Step 1*. For $\Psi={}^t(\varphi_1,\dots,\varphi_N)\in \bigoplus_{i=1}^NH^2(\mathbb R^3;\mathbb C^2)$, define $R$ and $K$ as in and .

We shall show $$\label{myeq4.3}
\langle\varphi_i,R\varphi_i\rangle-\langle\varphi_i,K\varphi_i\rangle\geq0,\ i=1,\dots,N,$$ when $\langle\varphi_i,\varphi_j\rangle=\delta_{ij}$. By the symmetry we can assume $i=1$. Let us introduce the notation $$[ ij\vert kl]:=\int d\mathbf xd\mathbf y\lvert x-y\rvert^{-1}\varphi_i^*(\mathbf x)\varphi_j(\mathbf x)\varphi_k^*(\mathbf y)\varphi_l(\mathbf y).$$ Then $$\langle\varphi_1,R\varphi_1\rangle-\langle\varphi_1,K\varphi_1\rangle=\sum_{j=1}^N\{[11\vert jj]-[1j\vert j1]\}=\sum_{j\neq 1}\{[11\vert jj]-[1j\vert j1]\}.$$ Thus we only need to prove $[11\vert 22]-[12\vert 21]\geq 0$. Set $$\tilde \Phi_{12}(\mathbf x,\mathbf y):=\varphi_1(\mathbf x)\varphi_2(\mathbf y)-\varphi_1(\mathbf y)\varphi_2(\mathbf x).$$ Then we can see that $$[11\vert 22]-[12\vert 21]=\frac{1}{2}\int d\mathbf xd\mathbf y\lvert x-y\rvert^{-1}\lvert\tilde \Phi_{12}(\mathbf x,\mathbf y)\rvert^2\geq 0.$$

Set $h:=-\Delta+V$. Then the Hartree-Fock equation is written as $$h\varphi_i+R\varphi_i-K\varphi_i=\epsilon_i\varphi_i.$$ Multiplying $\varphi_i^*$ and integrating the both sides we obtain by $$\epsilon_i=\langle\varphi_i,h\varphi_i\rangle+\langle\varphi_i,R\varphi_i\rangle-\langle\varphi_i,K\varphi_i\rangle\geq \inf\sigma(h)>-\infty,$$ where $\sigma(h)$ is the spectra of $h$. Thus, if there are infinitely many eigenvalues $\mathcal E^m=(\epsilon_1^m,\dots,\epsilon_N^m)\in (-\infty,\epsilon)^N,\ m=1,2,\dots$, we have $\mathcal E^m\in (\inf\sigma(h),\epsilon)^N$ and there exists a subsequence $\mathcal E^{m_k}$ such that $\epsilon_i^{m_k}$ converges to $\epsilon_i^0\leq \epsilon$. By Lemma \[convseq\] $\mathcal E^0:=(\epsilon_1^0,\dots,\epsilon_N^0)$ is an eigenvalue and we may assume $\mathcal E^m$ converges to $\mathcal E^0$ and the associated eigenfunction $\Psi^m:={}^t(\varphi_1^m,\dots,\varphi_N^m)$ converges in $\bigoplus_{i=1}^NH^2(\mathbb R^3;\mathbb C^2)$ to the eigenfunction $\Psi^0:={}^t(\varphi_1^0,\dots,\varphi_N^0)$ associated with $\mathcal E^0$.

*Step 2*. Let us define $$R_i(x):=\sum_{j\neq i}\int\lvert x-y\rvert^{-1}\varphi_j^*(\mathbf y)\varphi_j(\mathbf y)d\mathbf y,$$ $$(K_iw)(x,\sigma):=\sum_{j\neq i}\int\lvert x-y\rvert^{-1}\varphi_j^*(\mathbf y)w(\mathbf y)d\mathbf y\varphi_j(x,\sigma),$$ Then the Hartree-Fock equation is written as $$\label{myeq4.4}
h\varphi_i+R_i\varphi_i-K_i\varphi_i-\epsilon_i\varphi_i=0,$$ Let $\Phi$ be the Slater determinant of $(\varphi_1,\dots,\varphi_N)$ and $$f(\varphi_1,\dots,\varphi_N):=\langle \Phi,H\Phi\rangle.$$ Denote by $Y_1:=(\bigoplus _{i=1}^NH^2(\mathbb R^3;\mathbb C^2))\times \mathbb R^N$ and $Y_2:=(\bigoplus _{i=1}^NL^2(\mathbb R^3;\mathbb C^2))\times \mathbb R^N$ Banach spaces with the norm of product Banach spaces regarding $H^2(\mathbb R^3;\mathbb C^2)$ and $L^2(\mathbb R^3;\mathbb C^2)$ as real Banach spaces with respect to multiplication by real numbers. Denoting $\Psi^j:={}^t(\varphi_1^j,\dots,\varphi_N^j)$ and $\mathcal E^j:=(\epsilon_1^j,\dots,\epsilon_N^j),\ j=1,2$ define a functional $\tilde f:Y_1\to \mathbb R$ by $$\tilde f(\Psi,\mathcal E):=f(\varphi_1,\dots,\varphi_N)-\sum_{i=1}^N\epsilon_i(\lVert \varphi_i\rVert^2-1),$$ and a bilinear form $\langle\langle\cdot,\cdot\rangle\rangle$ on $Y_1$ and $Y_2$ by $$\langle\langle[\Psi^1,\mathcal E^1],[\Psi^2,\mathcal E^2]\rangle\rangle:=\sum_{i=1}^N2\mathrm{Re}\, \langle \varphi_i^1,\varphi_i^2\rangle+\sum_{i=1}^N\epsilon_i^1\epsilon_i^2.$$ We also define a mapping $\tilde F: Y_1\to Y_2$ by $$\tilde F(\Psi,\mathcal E):=(\tilde F_1(\Psi,\mathcal E),\dots,\tilde F_N(\Psi,\mathcal E),\lVert \varphi_1\rVert^2-1,\dots,\lVert \varphi_N\rVert^2-1),$$ $$\tilde F_i(\Psi,\mathcal E):=h\varphi_i+R_i\varphi_i-K_i\varphi_i-\epsilon_i\varphi_i.$$ Then if $\Psi$ is an eigenfunction of the Hartree-Fock equation associated with $\mathcal E$, $\tilde F(\Psi,\mathcal E)=0$ holds. Moreover, we have $$\tilde f'(\Psi,\mathcal E)[\Psi^1,\mathcal E^1]=\langle\langle[\Psi^1,\mathcal E^1],\tilde F(\Psi,\mathcal E)\rangle\rangle,$$ where $\tilde f'(\Psi,\mathcal E)$ is the Frechet derivative of $\tilde f$. We can also see that if $\tilde F(\Psi,\mathcal E)=0$, then $\tilde f(\Psi,\mathcal E)=f(\Psi)$. Thus the eigenfunctions of the Hartree-Fock equation are critical points of $\tilde f$ and the corresponding critical values are Hartree-Fock energies.

*Step 3*. Now let us consider the Frechet derivative of $F: \bigoplus_{i=1}^NH^2(\mathbb R^3;\mathbb C^2)\to\bigoplus_{i=1}^NL^2(\mathbb R^3;\mathbb C^2)$ defined by $$F(\varphi_1,\dots,\varphi_N):=(\tilde F_1(\varphi_1,\dots,\varphi_N,\mathcal E^0),\dots,\tilde F_N(\varphi_1,\dots,\varphi_N,\mathcal E^0)).$$ For a mapping $G: \bigoplus_{i=1}^NH^2(\mathbb R^3;\mathbb C^2)\to L^2(\mathbb R^3;\mathbb C^2)$ we set $$G_{,j}w_j:=\lim_{t\to0}[G(\varphi_1^0,\dots,\varphi_i^0+w_i,\dots,\varphi_N^0)-G(\varphi_1^0,\dots,\varphi_i^0,\dots,\varphi_N^0)]/t.$$ Then $$\begin{aligned}
&[R_i\varphi_i]_{,i}w_i=R_iw_i,\\
&[R_i\varphi_i]_{,j}w_j=L_{ij}w_j+\bar L_{ij}w_j,\ j\neq i,\end{aligned}$$ where $$\begin{aligned}
(L_{ij}w)(x,\sigma)&:=\int\lvert x-y\rvert^{-1}(\varphi_j^0)^*(\mathbf y)w(\mathbf y)d\mathbf y\varphi_i^0(x,\sigma),\\
(\bar L_{ij}w)(x,\sigma)&:=\int\lvert x-y\rvert^{-1}\varphi_j^0(\mathbf y)w^*(\mathbf y)d\mathbf y\varphi_i^0(x,\sigma),\end{aligned}$$ and $$\begin{aligned}
&[K_i\varphi_i]_{,i}w_i=K_iw_i,\\
&[K_i\varphi_i]_{,j}w_j=T_{ij}w_j+\bar L_{ji}w_j,\ j\neq i,\end{aligned}$$ where $$T_{ij}(x):=\int\lvert x-y\rvert^{-1}(\varphi_j^0)^*(\mathbf y)\varphi_i^0(\mathbf y)d\mathbf y.$$ Set $W:={}^t(w_1,\dots,w_N)$. We define mappings $$\mathcal R,\mathcal T: \bigoplus_{i=1}^NH^2(\mathbb R^3;\mathbb C^2)\to\bigoplus_{i=1}^NL^2(\mathbb R^3;\mathbb C^2),$$ as $$\begin{aligned}
(\mathcal RW)_i&:=R_iw_i,\\
(\mathcal TW)_i&:=\sum_{j\neq i}T_{ij}w_j.\end{aligned}$$ We shall show that $\mathcal R-\mathcal T$ is a positive operator as an operator on the Hilbert space $\bigoplus_{i=1}^NL^2(\mathbb R^3;\mathbb C^2)$. For that purpose we introduce a function $$\begin{aligned}
&\tilde\Phi(\mathbf x_1,\dots,\mathbf x_{N})\\
&\quad:=2^{-1/2}\sum_{i=1}^N\sum_{j\neq i}(w_i(\mathbf x_1)\varphi_j(\mathbf x_2)-w_j(\mathbf x_1)\varphi_i(\mathbf x_2))\check \Phi_{ij}(\mathbf x_3,\dots,\mathbf x_N),\end{aligned}$$ where $$\begin{aligned}
&\check \Phi_{ij}(\mathbf x_3,\dots,\mathbf x_N):=\varphi_1(\mathbf x_{\kappa_{ij}(1)})\cdots\varphi_{i-1}(\mathbf x_{\kappa_{ij}(i-1)})\varphi_{i+1}(\mathbf x_{\kappa_{ij}(i+1)})\cdot\\
&\hspace{50pt}\cdots\varphi_{j-1}(\mathbf x_{\kappa_{ij}(j-1)})\varphi_{j+1}(\mathbf x_{\kappa_{ij}(j+1)})\cdots\varphi_N(\mathbf x_{\kappa_{ij}(N)}).\end{aligned}$$ Here $\kappa_{ij}$ is a map from $\{1,\dots, \check i,\dots,\check j,\dots,N\}$ onto $\{3,\dots,N\}$ where $\check i$ means that $i$ is excluded. We use the notation as $[\tilde i\tilde j\vert kl]$ where indices such as $\tilde i$ indicate that $\varphi_i$ is replaced by $w_i$, e.g. $$[\tilde i\tilde j\vert  k l]:=\int d\mathbf xd\mathbf y\lvert x-y\rvert^{-1}w_i^*(\mathbf x)w_j(\mathbf x)\varphi_k^*(\mathbf y)\varphi_l(\mathbf y).$$ Then we can calculate as $$\langle W,(\mathcal R-\mathcal T)W\rangle=\sum_{i=1}^N\sum_{j\neq i}^N\{[\tilde i\tilde i\vert jj]-[\tilde i\tilde j\vert ji]\}.$$ On the other hand we can calculate as $$\begin{aligned}
&\int d\mathbf x_1\cdots d\mathbf x_N\lvert x_1-x_2\rvert^{-1}\lvert \tilde \Phi(\mathbf x_1,\dots,\mathbf x_{N})\rvert^2\\
&\quad=2^{-1}\sum_{i=1}^N\sum_{j\neq i}\int d\mathbf x_1d\mathbf x_2\lvert x_1-x_2\rvert^{-1}\lvert w_i(\mathbf x_1)\varphi_j(\mathbf x_2)-w_j(\mathbf x_1)\varphi_i(\mathbf x_2)\rvert^2\\
&\quad=2^{-1}\sum_{i=1}^N\sum_{j\neq i}\{[\tilde i\tilde i\vert jj]+[\tilde j\tilde j\vert ii]-[\tilde i\tilde j\vert ji]-[\tilde j\tilde i\vert ij]\}\\
&\quad=\sum_{i=1}^N\sum_{j\neq i}\{[\tilde i\tilde i\vert jj]-[\tilde i\tilde j\vert ji]\}\\
&\quad=\langle W,(\mathcal R-\mathcal T)W\rangle.\end{aligned}$$ Since the left-hand side is positive, we can see that $\langle W,(\mathcal R-\mathcal T)W\rangle\geq0$.

Next we consider $h-\epsilon_i$. Denote the resolution of identity of $h$ by $E(\lambda)$. Then we can decompose $h$ as $$h=hE(\epsilon/2)+h(1-E(\epsilon/2)).$$ Since $\inf \sigma_{ess}(h)=0$, $hE(\epsilon/2)$ is a compact operator, where $\sigma_{ess}(h)$ is the essential spectra of $h$. Moreover, we have an inequality $h(1-E(\epsilon/2))\geq \epsilon/2$ of operators. Since $\epsilon_i\leq\epsilon$, the operator $h(1-E(\epsilon/2))-\epsilon_i$ is positive definite. Therefore, the operator $\mathcal H: \bigoplus_{i=1}^NH^2(\mathbb R^3;\mathbb C^2)\to\bigoplus_{i=1}^NL^2(\mathbb R^3;\mathbb C^2)$ defined by $$(\mathcal HW)_i:=(h-\epsilon_i)w_i,$$ is decomposed as a sum of positive definite operator $\mathcal H_1$ and a compact operator $\mathcal H_2$: $\mathcal H=\mathcal H_1+\mathcal H_2$. If we also define $\mathcal L$ and $\bar{\mathcal L}$ by $$\begin{aligned}
(\mathcal LW)_i=\sum_{j\neq 1}L_{ij}w_j,
(\bar{\mathcal L}W)_i=\sum_{j\neq 1}\bar L_{ij}w_j,\end{aligned}$$ the Frechet derivative of $F$ at $\Psi^0$ is $$F'(\Psi^0)=\mathcal H_1+\mathcal H_2+\mathcal R-\mathcal T+\mathcal L+\bar{\mathcal L}+{}^t\bar{\mathcal L}=\mathcal J+\mathcal M,$$ where $\mathcal J:=\mathcal H_1+\mathcal R-\mathcal T$ and $\mathcal M:=\mathcal H_2+\mathcal L+\bar{\mathcal L}+{}^t\bar{\mathcal L}$. Since $\mathcal R-\mathcal T$ is positive, $\mathcal J$ is positive definite and invertible. Moreover, since $L_{ij}$ and $\bar L_{ij}$ are compact operators, $\mathcal L+\bar{\mathcal L}+{}^t\bar{\mathcal L}$ is compact and $\mathcal M$ is also compact.

*Step 4*. We can easily see that $\tilde F$ is real-analytic as an operator between real Banach spaces $Y_1$ and $Y_2$. We can also see that $$\tilde F'(\Psi^0,\mathcal E^0)=J+M,$$ where $$\begin{aligned}
J[\Psi,\mathcal E]&:=[\mathcal J\Psi,\mathcal E],\\
M[\Psi,\mathcal E]&:=[\mathcal M\Psi-\mathcal E\Psi^0,2\mathrm{Re}\, \langle\varphi_1,\varphi_1^0\rangle-\epsilon_1,\dots,2\mathrm{Re}\, \langle\varphi_N,\varphi_N^0\rangle-\epsilon_N],\end{aligned}$$ where $\mathcal E\Psi^0={}^t(\epsilon_1\varphi_1^0,\dots,\epsilon_N\varphi_N^0)$. Since $M$ is a compact operator and $J$ is an isomorphism, we can apply Lemma \[onepoint\] for $\tilde f$ at $[\Psi^0,\mathcal E^0]$. Thus it is proved that there exists a neighborhood $U^0$ of $[\Psi^0,\mathcal E^0]$ such that $\tilde f(\mathcal B_{\tilde f}\cap U^0)$ is one-point set. Now remembering it has been assumed that the Hartree-Fock energy $f(\Psi^m)=\tilde f(\Psi^m)$ satisfy $f(\Psi^{m_1})\neq f(\Psi^{m_2})$ for $m_1\neq m_2$ from the beginning of the proof, and noting $f(\Psi^m)\to f(\Psi^0)$ by $\Psi^m\to\Psi^0$ in $\bigoplus_{i=1}^NH^2(\mathbb R^3,\mathbb C^2)$ we have a contradiction.

[99]{} S. Agmon, Lectures on exponential decay of solutions of second-order elliptic equations: Bounds on eigenfunctions of $N$-body Schrödinger operations, Princeton University Press, (1982) V. Fock, Näherungsmethode zur L" osing des quantenmechanischen Mehrkörperproblems, [*Z. phys.*]{} [**61**]{} (1930), 126–148. S. Fučik, J. Nečas, J. Souček, and V. Souček, Upper bound for the number of critical levels for nonlinear operators in Banach spaces of the type of second order nonlinear partial differential operators, [*J. Funct. Anal.*]{} [**11**]{} (1972), 314–333. D. Hartree, The wave mechanics of an atom with a non-coulomb central field. Part I. Theory and methods, [*Proc. Comb. Phil. Soc.*]{} [**24**]{} (1928), 89–132. T. Kato, Perturbation theory for linear operators, Berlin-Heidelberg-New York, Springer, (1976) M. Lewin, Existence of Hatree-Fock excited states for atoms and molecules, [*Lett. Math. Phys.*]{} [**108**]{} (2018), 985–1006. E. H. Lieb and B. Simon, The Hartree-Fock theory for Coulomb systems, [*Commun. math. Phys.*]{} [**53**]{} (1977), 185–194. P. L. Lions, Solutions of Hartree-Fock equations for Coulomb systems, [*Commun. math. Phys.*]{} [**109**]{} (1987), 33–97. M. Reeken, General theorem on bifurcation and its applications to the Hartree equation of Helium atom, [*J. Math. Phys.*]{} [**11**]{} (1970), 2505–2512. J. C. Slater, A note on Hartree’s method, [*Phys, Rev.*]{} [**35**]{} (1930), 210–211. J. Souček, and V. Souček, Morse-Sard theorem for real-analytic functions, [*Comment. Math. Univ. Carolinae*]{} [**13**]{} (1972), 45–51. C. A. Stuart, Existence theory for the Hartree equation, [*Arch. Rational Mech. Anal.*]{} [**51**]{} (1973), 60–69. J. H. Wolkowisky, Existence of solutions of the Hartree equations for N electrons an application of the Schauder-Tychonoff theorem, [*Indiana Univ. Math. J.*]{} [**22**]{} (1972), 551–568.
