Editing Slater determinant

{{short description|Function that can be used to build the wave function of a multi-fermionic system}}
In [[quantum mechanics]], a '''Slater determinant'''  is an expression that describes the [[wave function]] of a multi-[[fermion]]ic system. It satisfies [[Skew-symmetric matrix|anti-symmetry]] requirements, and consequently the [[Pauli principle]], by changing [[Plus and minus signs|sign]] upon exchange of two fermions.<ref name="Atkins">Molecular Quantum Mechanics Parts I and II: An Introduction to QUANTUM CHEMISTRY (Volume 1), P. W. Atkins, Oxford University Press, 1977, {{ISBN|0-19-855129-0}}.</ref> Only a small subset of all possible many-body fermionic wave functions can be written as a single Slater determinant, but those form an important and useful subset because of their simplicity.

The Slater determinant arises from the consideration of a wave function for a collection of electrons, each with a wave function known as the [[spin-orbital]] <math>\chi(\mathbf{x})</math>, where <math>\mathbf{x}</math> denotes the position and spin of a single electron. A Slater determinant containing two electrons with the same spin orbital would correspond to a wave function that is zero everywhere.

The Slater determinant is named for [[John C. Slater]], who introduced the determinant in 1929 as a means of ensuring the antisymmetry of a many-electron wave function,<ref>{{cite journal |last1=Slater |first1=J. |title=The Theory of Complex Spectra |journal=Physical Review |volume=34 |issue=2 |pages=1293–1322 |year=1929  |doi=10.1103/PhysRev.34.1293 |bibcode = 1929PhRv...34.1293S }}</ref> although the wave function in the determinant form first appeared independently in Heisenberg's<ref>{{cite journal |last1 = Heisenberg |first1 = W. |title = Mehrkörperproblem und Resonanz in der Quantenmechanik |journal = Zeitschrift für Physik |year = 1926 |volume = 38 |issue = 6–7 |pages = 411–426 |doi= 10.1007/BF01397160 |bibcode = 1926ZPhy...38..411H |s2cid = 186238286 }}</ref> and Dirac's<ref>{{cite journal |last1 = Dirac |first1 = P. A. M. |title = On the Theory of Quantum Mechanics |journal = Proceedings of the Royal Society A |year = 1926 |volume = 112 |issue = 762 |pages = 661–677 |doi= 10.1098/rspa.1926.0133 |bibcode = 1926RSPSA.112..661D |doi-access = free }}</ref><ref>{{Cite web |title=Slater determinant in nLab |url=https://ncatlab.org/nlab/show/Slater+determinant#references |access-date=2023-11-08 |website=ncatlab.org}}</ref> articles three years earlier.

== Definition ==

=== Two-particle case ===

The simplest way to approximate the wave function of a many-particle system is to take the product of properly chosen [[Orthogonality (quantum mechanics)|orthogonal]] wave functions of the individual particles. For the two-particle case with coordinates <math>\mathbf{x}_1</math> and <math>\mathbf{x}_2</math>, we have

:<math>
\Psi(\mathbf{x}_1, \mathbf{x}_2) = \chi_1(\mathbf{x}_1) \chi_2(\mathbf{x}_2).
</math>

This expression is used in the [[Hartree–Fock method#Brief history|Hartree method]] as an [[ansatz]] for the many-particle wave function and is known as a [[Hartree product]]. However, it is not satisfactory for [[fermions]] because the wave function above is not antisymmetric under exchange of any two of the fermions, as it must be according to the [[Pauli exclusion principle]]. An antisymmetric wave function can be mathematically described as follows:

:<math>
\Psi(\mathbf{x}_1, \mathbf{x}_2) = -\Psi(\mathbf{x}_2, \mathbf{x}_1).
</math>

This does not hold for the Hartree product, which therefore does not satisfy the Pauli principle. This problem can be overcome by taking a [[linear combination]] of both Hartree products:

:<math>
\begin{aligned}
 \Psi(\mathbf{x}_1, \mathbf{x}_2) &= \frac{1}{\sqrt{2}} \{\chi_1(\mathbf{x}_1) \chi_2(\mathbf{x}_2) - \chi_1(\mathbf{x}_2) \chi_2(\mathbf{x}_1)\} \\
  &= \frac{1}{\sqrt2}\begin{vmatrix} \chi_1(\mathbf{x}_1) & \chi_2(\mathbf{x}_1) \\ \chi_1(\mathbf{x}_2) & \chi_2(\mathbf{x}_2) \end{vmatrix},
\end{aligned}
</math>

where the coefficient is the [[normalization factor]]. This wave function is now antisymmetric and no longer distinguishes between fermions (that is, one cannot indicate an ordinal number to a specific particle, and the indices given are interchangeable). Moreover, it also goes to zero if any two spin orbitals of two fermions are the same. This is equivalent to satisfying the Pauli exclusion principle.

=== Multi-particle case ===

The expression can be generalised to any number of fermions by writing it as a [[determinant]]. For an ''N''-electron system, the Slater determinant is defined as<ref name="Atkins" /><ref name="Szabo">{{cite book
  | last1 = Szabo
  | first1 = A.
  | last2= Ostlund| first2= N. S.
  | title = Modern Quantum Chemistry
  | publisher = Dover Publishing | year = 1996 | location = Mineola, New York
  | isbn = 0-486-69186-1}}</ref>

:<math>
\begin{aligned}
 \Psi(\mathbf{x}_1, \mathbf{x}_2, \ldots, \mathbf{x}_N) &=
  \frac{1}{\sqrt{N!}}
  \begin{vmatrix} \chi_1(\mathbf{x}_1) & \chi_2(\mathbf{x}_1) & \cdots & \chi_N(\mathbf{x}_1) \\
                      \chi_1(\mathbf{x}_2) & \chi_2(\mathbf{x}_2) & \cdots & \chi_N(\mathbf{x}_2) \\
                      \vdots & \vdots & \ddots & \vdots \\
                      \chi_1(\mathbf{x}_N) & \chi_2(\mathbf{x}_N) & \cdots & \chi_N(\mathbf{x}_N)
  \end{vmatrix} \\
 &\equiv | \chi _1, \chi _2, \cdots, \chi _N \rangle \\
 &\equiv | 1, 2, \dots, N \rangle,
\end{aligned}
</math>

where the last two expressions use a shorthand for Slater determinants: The normalization constant is implied by noting the number N, and only the one-particle wavefunctions (first shorthand) or the indices for the fermion coordinates (second shorthand) are written down. All skipped labels are implied to behave in ascending sequence. The linear combination of Hartree products for the two-particle case is identical with the Slater determinant for ''N'' = 2. The use of Slater determinants ensures an antisymmetrized function at the outset. In the same way, the use of Slater determinants ensures conformity to the [[Pauli principle]]. Indeed, the Slater determinant vanishes if the set  <math>\{\chi_i\}</math> is [[linearly dependent]]. In particular, this is the case when two (or more) spin orbitals are the same. In chemistry one expresses this fact by stating that no two electrons with the same spin can occupy the same spatial orbital.

== Example: Matrix elements in a many electron problem ==
{{See also|Slater-Condon rules}}
Many properties of the Slater determinant come to life with an example in a non-relativistic many electron problem.<ref name="ReferenceA">Solid State Physics - Grosso Parravicini - 2nd edition pp.140-143</ref>
* ''The one particle terms of the Hamiltonian will contribute in the same manner as for the simple Hartree product, namely the energy is summed and the states are independent''
* ''The multi-particle terms of the Hamiltonian will introduce exchange term to lower of the energy for the anti-symmetrized wave function''

Starting from a [[molecular Hamiltonian]]:
<math display="block">\hat{H}_\text{tot} = \sum_i \frac{\mathbf{p}^2_i}{2 m} + \sum_I \frac{\mathbf{P}^2_I}{2 M_I} + \sum_i V_\text{nucl}(\mathbf{r_i}) + \frac{1}{2}\sum_{i \ne j} \frac{e^2}{|\mathbf{r}_i-\mathbf{r}_j|} + \frac{1}{2}\sum_{I \ne J} \frac{Z_I Z_J e^2}{|\mathbf{R}_I-\mathbf{R}_J|}</math>
where <math>\mathbf{r}_i</math> are the electrons and <math>\mathbf{R}_I</math>are the nuclei and
: <math>V_\text{nucl}(\mathbf{r})= - \sum_I \frac{Z_I e^2}{|\mathbf{r}-\mathbf{R}_I|}</math>

For simplicity we freeze the nuclei at equilibrium in one position and we remain with a simplified Hamiltonian
: <math>\hat{H}_e = \sum^N_i \hat{h}(\mathbf{r}_i) + \frac{1}{2}\sum^N_{i \ne j} \frac{e^2}{r_{ij}} </math>
where
: <math>\hat{h}(\mathbf{r}) = \frac{\hat{\mathbf{p}}^2}{2m} + V_\text{nucl}(\mathbf{r})</math>

and where we will distinguish in the Hamiltonian between the first set of terms as <math>\hat{G}_1</math> (the "1" particle terms) 
and the last term <math>\hat{G}_2</math> (the "2" particle term) which contains exchange term for a Slater determinant.

: <math>\hat{G}_1 =\sum^N_i \hat{h}(\mathbf{r}_i) </math>
: <math>\hat{G}_2 =\frac{1}{2} \sum^N_{i \ne j} \frac{e^2}{r_{ij}}</math>

The two parts will behave differently when they have to interact with a Slater determinant wave function. We start to compute the expectation values of one-particle terms
: <math>\langle\Psi_0 |G_1 | \Psi_0\rangle = \frac{1}{N!}\langle \det\{\psi_1 ... \psi_N\}|G_1|\det\{\psi_1 ... \psi_N\}\rangle </math>
In the above expression, we can just select the identical permutation in the determinant in the left part, since all the other N! − 1 permutations would give the same result as the selected one. We can thus cancel N! at the denominator
: <math>\langle\Psi_0 |G_1 | \Psi_0\rangle = \langle\psi_1 ... \psi_N|G_1|\det\{\psi_1 ... \psi_N\}\rangle</math>
Because of the orthonormality of spin-orbitals it is also evident that only the identical permutation survives in the determinant on the right part of the above matrix element
: <math>\langle\Psi_0 |G_1 | \Psi_0\rangle = \langle\psi_1 ... \psi_N|G_1|\psi_1 ... \psi_N\rangle</math>
This result shows that the anti-symmetrization of the product does not have any effect for the one particle terms and it behaves as it would do in the case of the simple Hartree product.

And finally we remain with the trace over the one-particle Hamiltonians
: <math>\langle\Psi_0 |G_1 | \Psi_0\rangle = \sum_i \langle\psi_i|h|\psi_i\rangle</math>
Which tells us that to the extent of the one-particle terms the wave functions of the electrons are independent of each other and the expectation value of total system is given by the sum of expectation value of the single particles.

For the two-particle terms instead
: <math>\langle\Psi_0 |G_2 | \Psi_0\rangle = \frac{1}{N!}\langle\det\{\psi_1 ... \psi_N\}|G_2|\det\{\psi_1 ... \psi_N\}\rangle = \langle\psi_1 ... \psi_N|G_2|\det\{\psi_1 ... \psi_N\}\rangle</math>

If we focus on the action of one term of <math>G_2</math>, it will produce only the two terms
: <math> \langle\psi_1(r_1,\sigma_1) ... \psi_N(r_N, \sigma_N) |\frac{e^2}{r_{12}}|\mathrm{det}\{\psi_1(r_1,\sigma_1) ... \psi_N(r_N,\sigma_N)\}\rangle= \langle\psi_1\psi_2|\frac{e^2}{r_{12}}|\psi_1\psi_2\rangle - \langle\psi_1\psi_2|\frac{e^2}{r_{12}}|\psi_2\psi_1\rangle</math>

And finally
<math display="block">\langle\Psi_0 |G_2 | \Psi_0\rangle = \frac{1}{2}\sum_{i\ne j}\left[ \langle\psi_i \psi_j | \frac{e^2}{r_{ij}} |
\psi_i \psi_j \rangle - \langle\psi_i \psi_j | \frac{e^2}{r_{ij}} |
\psi_j \psi_i \rangle \right] </math>

which instead is a mixing term. The first contribution is called the "coulomb" term or "coulomb" integral and the second is the "exchange" term or exchange integral. Sometimes different range of index in the summation is used <math display="inline">\sum_{ij}</math> since the Coulomb and exchange contributions exactly cancel each other for <math>i = j</math>.

It is important to notice explicitly that the exchange term, which is always positive for local spin-orbitals,<ref>See appendix I in {{ cite journal|doi=10.1103/RevModPhys.23.69 | title=New Developments in Molecular Orbital Theory |first1= C. C. J. |last1= Roothaan | page=69 | volume=23 |year=1951 |journal= Reviews of Modern Physics | issue=69| url=http://elib.bsu.by/handle/123456789/154388 }}</ref> is absent in the simple Hartree product. Hence the electron-electron repulsive energy <math>\langle\Psi_0 |G_2 |\Psi_0\rangle</math> on the antisymmetrized product of spin-orbitals is always lower than the electron-electron repulsive energy on the simple Hartree product of the same spin-orbitals. Since exchange bielectronic integrals are different from zero only for spin-orbitals with parallel spins, we link the decrease in energy with the physical fact that electrons with parallel spin are kept apart in real space in Slater determinant states.

== As an approximation ==

Most fermionic wavefunctions cannot be represented as a Slater determinant.  The best Slater approximation to a given fermionic wave function can be defined to be the one that maximizes the [[orbital overlap|overlap]] between the Slater determinant and the target wave function.<ref>{{ cite journal|doi=10.1103/PhysRevA.89.012504 | title=Optimal multiconfiguration approximation of an ''N''-fermion wave function |first1=J. M. |last1=Zhang | first2 = Marcus | last2 = Kollar |page=012504 | volume=89 |year=2014 |journal= Physical Review A | issue=1 |arxiv = 1309.1848 |bibcode = 2014PhRvA..89a2504Z| s2cid=17241999 }}</ref> The maximal overlap is a geometric measure of [[quantum entanglement|entanglement]] between the fermions.

A single Slater determinant is used as an approximation to the electronic wavefunction in [[Hartree–Fock method|Hartree–Fock theory]]. In more accurate theories (such as [[configuration interaction]] and [[MCSCF]]), a linear combination of Slater determinants is needed.

== Discussion ==

The word "'''detor'''" was proposed by [[S. Francis Boys|S. F. Boys]] to refer to a Slater determinant of orthonormal orbitals,<ref>{{ cite journal| title=Electronic wave functions I. A general method of calculation for the stationary states of any molecular system|author-link=S. Francis Boys|first=S. F. |last=Boys|page=542| volume=A200 |year=1950|journal= Proceedings of the Royal Society|issue=1063|doi=10.1098/rspa.1950.0036|bibcode=1950RSPSA.200..542B|s2cid=122709395|url=http://elib.bsu.by/handle/123456789/154387}}</ref> but this term is rarely used.

Unlike [[fermions]] that are subject to the Pauli exclusion principle, two or more [[bosons]] can occupy the same single-particle quantum state. Wavefunctions describing systems of identical [[bosons]] are symmetric under the exchange of particles and can be expanded in terms of [[Permanent (mathematics)|permanent]]s.

==See also ==
* [[Antisymmetrizer]]
* [[Atomic orbital|Electron orbital]]
* [[Fock space]]
* [[Quantum electrodynamics]]
* [[Quantum mechanics]]
* [[Physical chemistry]]
* [[Hund's rule]]
* [[Hartree–Fock method]]

==References==
{{reflist|2}}

==External links==
* [http://www.cond-mat.de/events/correl13/manuscripts/koch.pdf Many-Electron States] in E.  Pavarini, E. Koch, and U. Schollwöck: Emergent Phenomena in Correlated Matter, Jülich 2013, {{ISBN|978-3-89336-884-6}}

[[Category:Quantum mechanics]]
[[Category:Quantum chemistry]]
[[Category:Theoretical chemistry]]
[[Category:Computational chemistry]]
[[Category:Determinants]]
[[Category:Pauli exclusion principle]]