Pseudo-Relativistic Hartree-Fock Method

• Pseudo-relativistic Hartree-Fock is a method that retains relativistic one-body dynamics while enforcing antisymmetry through explicit exchange terms.
• It provides accurate control of binding energies, ionization limits, and density distributions across atomic, nuclear, and molecular systems.
• The approach unifies various formulations—from pseudo-differential kinetic operators in atoms to covariant Dirac equations in nuclei—into a coherent computational framework.

Searching arXiv for recent and foundational papers on pseudo-relativistic Hartree-Fock and closely related relativistic Hartree-Fock formulations. Pseudo-relativistic Hartree-Fock denotes a class of mean-field formalisms in which relativistic kinematics or Dirac structure are retained at the one-particle level while antisymmetry is enforced through Hartree-Fock exchange. In the atomic mathematical formulation, the defining ingredient is the pseudo-relativistic kinetic operator $\alpha^{-1}(\sqrt{-\alpha^2\Delta+1}-1)$ coupled to Coulomb attraction, direct Coulomb repulsion, and exchange (Dall'Acqua et al., 2010). In nuclear and molecular implementations, the same general Hartree-Fock logic is realized through covariant Dirac-based equations, nonlocal Fock terms, and practical reorganizations such as relativistic Hartree-Fock with PKA1, density-dependent RMF mappings, and unified 4C/Q4C/X2C Hamiltonians (Li et al., 2019, Gmuca et al., 2021, Liu, 2023).

1. Conceptual scope and shared structure

Across the literature, the method appears in several technically distinct but structurally related forms. In each case, the one-body dynamics are relativistic or pseudo-relativistic, and the many-body approximation is Hartree-Fock rather than pure Hartree: direct mean fields are supplemented by exchange terms, either kept explicitly or absorbed into effective couplings.

Context	One-body dynamics	Exchange treatment
Atomic pseudo-relativistic HF	$\alpha^{-1}(\sqrt{-\alpha^2\Delta+1}-1)$ with Coulomb attraction	Explicit direct and exchange Coulomb energies in the density-matrix functional
Nuclear RHF/DHF	Dirac spinors with scalar, vector, tensor self-energies	Explicit Fock terms, including $\pi$-meson and $\rho$-tensor couplings through exchange
Molecular 4C/Q4C/X2C HF	Four-component, quasi-four-component, or exact-two-component relativistic HF equations	Exchange retained in a unified Hamiltonian structure after pmf/MDM and 1CSC approximations

In the atomic model, the basic unknown is a one-particle density matrix $\gamma$ satisfying $0\le \gamma\le I$ and $\operatorname{Tr}[T(p)\gamma]<\infty$. In nuclear RHF, the basic object is the Dirac spinor with upper and lower radial components, and pseudo-spin structure is analyzed through the lower component. In molecular relativistic HF, the basic object is the relativistic Fock matrix in 4C, Q4C, or X2C representation. This suggests that “pseudo-relativistic Hartree-Fock” is best understood as a methodological family rather than a single equation.

2. Atomic pseudo-relativistic Hartree-Fock functional

For an atom with nuclear charge $Z$ and $N$ electrons, in units where $\hbar=m=e=1$, the many-body Hamiltonian is given by

$\alpha^{-1}(\sqrt{-\alpha^2\Delta+1}-1)$0

The pseudo-relativistic aspect lies in the kinetic term: the Laplacian is replaced by a square-root pseudo-differential operator. In Fourier variables,

$\alpha^{-1}(\sqrt{-\alpha^2\Delta+1}-1)$1

The Hartree-Fock description is expressed through the density-matrix functional

$\alpha^{-1}(\sqrt{-\alpha^2\Delta+1}-1)$2

with

$\alpha^{-1}(\sqrt{-\alpha^2\Delta+1}-1)$3

$\alpha^{-1}(\sqrt{-\alpha^2\Delta+1}-1)$4

and

$\alpha^{-1}(\sqrt{-\alpha^2\Delta+1}-1)$5

For a Slater determinant $\alpha^{-1}(\sqrt{-\alpha^2\Delta+1}-1)$6, the corresponding $\alpha^{-1}(\sqrt{-\alpha^2\Delta+1}-1)$7 is the orthogonal projection onto $\alpha^{-1}(\sqrt{-\alpha^2\Delta+1}-1)$8, and the ground-state energy is

$\alpha^{-1}(\sqrt{-\alpha^2\Delta+1}-1)$9

The direct term $\pi$0 is the Hartree contribution; the exchange term $\pi$1 is the defining Hartree-Fock correction (Dall'Acqua et al., 2010).

A central structural restriction is the subcritical condition

$\pi$2

The estimate

$\pi$3

implies that this regime guarantees stability and a well-defined quadratic form. The existence theorem used in the cited work states that if

$\pi$4

then a minimizer exists (Dall'Acqua et al., 2010).

3. Rigorous large-$\pi$5 properties and analytical machinery

The atomic pseudo-relativistic Hartree-Fock model supports a relativistic version of the ionization conjecture. The main theorem states that if $\pi$6, then there exists a constant $\pi$7, depending only on $\pi$8, such that whenever a Hartree-Fock minimizer exists,

$\pi$9

Hence the maximal negative ionization $\rho$0 is bounded uniformly in $\rho$1 and $\rho$2 as long as $\rho$3 stays fixed below $\rho$4. The same work also proves that the ionization energy of a neutral atom,

$\rho$5

is bounded by a universal constant independent of $\rho$6 and $\rho$7 (Dall'Acqua et al., 2010).

The proof strategy compares Hartree-Fock to Thomas-Fermi theory. The neutral Thomas-Fermi potential is

$\rho$8

and the corresponding TF equation is

$\rho$9

Sommerfeld asymptotics give

$\gamma$0

for large $\gamma$1, providing the screening scale against which the HF density is measured.

The pseudo-relativistic analysis differs sharply from the nonrelativistic case because only $\gamma$2 is available, rather than $\gamma$3. The proof therefore uses Daubechies’ inequality,

$\gamma$4

together with a Daubechies-Lieb-Yau type inequality controlling $\gamma$5 for potentials with Coulomb singularities. A further technical ingredient is a pseudo-relativistic IMS-type formula,

$\gamma$6

whose nonlocal error operator $\gamma$7 has a kernel involving the Bessel function $\gamma$8. Space is then decomposed into near-nucleus, intermediate, and outer regions; the screened potential

$\gamma$9

is controlled iteratively, leading to the excess-charge bound. The same analysis yields a uniform bound on the HF radius of the outer $0\le \gamma\le I$0 electrons with the same $0\le \gamma\le I$1-type behavior as in the nonrelativistic case, and an explicit decay estimate for the difference between the HF and TF mean-field potentials (Dall'Acqua et al., 2010).

4. Covariant relativistic Hartree-Fock in nuclei

In nuclear structure, relativistic Hartree-Fock is formulated within covariant density functional theory. The cited study uses the effective Lagrangian PKA1, which includes not only the usual Hartree mean fields but also Fock exchange terms, and incorporates $0\le \gamma\le I$2-meson and $0\le \gamma\le I$3-tensor couplings through the exchange channel. In spherical symmetry, the nucleon Dirac spinor is written as

$0\le \gamma\le I$4

where $0\le \gamma\le I$5 and $0\le \gamma\le I$6 are the upper and lower radial components (Li et al., 2019).

Variation of the RHF functional yields the radial Dirac equations

$0\le \gamma\le I$7

$0\le \gamma\le I$8

with vector, scalar, and tensor self-energies $0\le \gamma\le I$9, and nonlocal Fock terms $\operatorname{Tr}[T(p)\gamma]<\infty$0. To handle nonlocal exchange, the paper localizes these contributions through equivalent local potentials $\operatorname{Tr}[T(p)\gamma]<\infty$1.

Pseudo-spin symmetry is treated in the relativistic picture as the near-degeneracy of

$\operatorname{Tr}[T(p)\gamma]<\infty$2

which are recast as pseudo-spin partners with

$\operatorname{Tr}[T(p)\gamma]<\infty$3

The proton doublet of interest is

$\operatorname{Tr}[T(p)\gamma]<\infty$4

Exact pseudo-spin symmetry is associated with

$\operatorname{Tr}[T(p)\gamma]<\infty$5

so the lower component of the Dirac spinor carries the relevant pseudo-spin structure (Li et al., 2019).

To analyze the symmetry directly, the lower component $\operatorname{Tr}[T(p)\gamma]<\infty$6 is shown to satisfy a Schrödinger-like equation,

$\operatorname{Tr}[T(p)\gamma]<\infty$7

where

$\operatorname{Tr}[T(p)\gamma]<\infty$8

is the pseudo-centrifugal barrier. The Hartree and Fock pieces are decomposed as

$\operatorname{Tr}[T(p)\gamma]<\infty$9

$Z$0

with pseudo-spin-orbit potentials

$Z$1

and

$Z$2

The derivative of the direct scalar-plus-vector field is therefore the basic symmetry-breaking ingredient.

For even-even $Z$3 and $Z$4 isotones from $Z$5 (Si) to $Z$6 (Ni), with pairing treated by BCS and the finite-range Gogny D1S force, the central observable is the pseudo-spin orbital splitting

$Z$7

The systematic result is strong pseudo-spin-symmetry violation in sulfur, followed by progressive restoration from sulfur to nickel; by nickel, the $Z$8-shell orbit ordering can even invert. The tensor-force contributions introduced naturally by the Fock terms do affect the splitting, become weaker near spin saturation at $Z$9, and are enhanced again in spin-unsaturated systems, but they only partially explain the observed systematics. The overall trend is dominated by the Hartree terms. The decisive mechanism is the evolution of the proton central density profile: bubble-like in silicon, central-bumped in sulfur, and progressively central-flat toward nickel. This self-consistent density rearrangement changes $N$0, strengthens the pseudo-spin-orbit potential in sulfur, and restores pseudo-spin symmetry toward heavier isotones (Li et al., 2019).

5. Exchange, localization, and effective local reformulations

A recurrent theme in pseudo-relativistic Hartree-Fock is the tension between explicit Fock nonlocality and the desire for local or quasi-local effective descriptions. In symmetric nuclear matter, the exchange part of the linear Dirac-Hartree-Fock model has been evaluated in parameter-free closed form and reorganized as a density functional. Starting from

$N$1

the total potential energy density can be rewritten as

$N$2

where the effective couplings $N$3 and $N$4 are density-dependent and encode Fock correlations (Gmuca et al., 2021).

The exchange energy is built from analytically evaluated shape functions $N$5, $N$6, and $N$7. In the linear case, the resulting RMF couplings are generated entirely by exchange effects; in the nonlinear extension, meson self-interactions introduce additional medium dependence through functions such as

$N$8

A central implication is that Fock exchange need not be discarded to obtain an RMF-like local functional: it can be absorbed into density-dependent effective vertices. The paper further reports that the Fock-induced density dependence is around the $N$9 level in the couplings, that $\hbar=m=e=1$0-exchange is the dominant exchange contribution, and that the pion contributes only through exchange and is more visible at larger Fermi momentum (Gmuca et al., 2021).

In molecular relativistic Hartree-Fock, a related but computationally distinct reorganization is developed for 4C, Q4C, and X2C theories. The key assumption is the atomic nature of the small component. In an RKB basis,

$\hbar=m=e=1$1

Because small components are strongly localized around nuclei, the model density matrix approximation replaces the molecular small-component density by a superposition of atomic ones,

$\hbar=m=e=1$2

leading to a static pre-molecular mean field

$\hbar=m=e=1$3

The Fock matrix is then partitioned so that, after the pmf correction is folded into the one-electron term, only the nonrelativistic-like two-electron term built from the large-component density is iterated in the SCF cycle (Liu, 2023).

A complementary simplification is the one-center small-component approximation. It neglects interatomic overlaps of small-component basis functions and reduces relativistic integral classes according to relations such as

$\hbar=m=e=1$4

$\hbar=m=e=1$5

$\hbar=m=e=1$6

Under the MDM/pmf and 1CSC approximations, the 4C, Q4C, and X2C mean-field and many-electron Hamiltonians are stated to share precisely the same structure and accuracy, and the same normal-ordered post-HF Hamiltonian form is retained: $\hbar=m=e=1$7 The paper explicitly states that the associated approximation errors are orders of magnitude smaller than truncation errors in the one- and many-particle bases and uncertainties of experimental measurements (Liu, 2023).

6. Related Hartree reductions, limitations, and common misconceptions

A common source of confusion is the distinction between pseudo-relativistic Hartree and pseudo-relativistic Hartree-Fock. The pseudo-relativistic Hartree equation

$\hbar=m=e=1$8

is directly relevant to pseudo-relativistic mean-field theory, but it does not include exchange or spin effects. For $\hbar=m=e=1$9, $\alpha^{-1}(\sqrt{-\alpha^2\Delta+1}-1)$00, the kernel is Coulomb and the model reduces to the standard pseudo-relativistic Hartree equation. The analysis proceeds through a local realization in the half-space,

$\alpha^{-1}(\sqrt{-\alpha^2\Delta+1}-1)$01

a variational functional in $\alpha^{-1}(\sqrt{-\alpha^2\Delta+1}-1)$02, a penalization scheme, minimax construction, and concentration-compactness. The resulting solutions exhibit single-spike concentration near the minimum set

$\alpha^{-1}(\sqrt{-\alpha^2\Delta+1}-1)$03

with exponential decay and convergence, after translation and rescaling, to a least-energy solution of the limiting autonomous problem (Cingolani et al., 2015).

This Hartree theory captures mean-field self-interaction, a scalar effective potential, and nonlocal interaction via convolution, but Hartree-Fock would add exchange terms, antisymmetry of fermionic states, spin degrees of freedom, and matrix-valued or projector-valued structure. The Hartree equation is therefore a reduced model rather than a substitute for pseudo-relativistic Hartree-Fock (Cingolani et al., 2015).

Two further misconceptions are corrected by the relativistic HF literature itself. First, in the nuclear PKA1 analysis, pseudo-spin-symmetry restoration is not primarily a tensor-force story: tensor-force contributions influence the splitting but cannot account for the full isotonic trend, whereas Hartree terms dominate the overall systematics (Li et al., 2019). Second, in molecular relativistic HF, 4C, Q4C, and X2C are not presented as unrelated theories in spirit; under a common treatment of small components and relativistic integrals, they can be written in the same algebraic form, differing mainly in representation and decoupling strategy (Liu, 2023).

Taken together, these works define pseudo-relativistic Hartree-Fock as a broad methodology for relativistic mean-field many-body theory. Its atomic form provides rigorous control of binding, ionization, and screening in the large-$\alpha^{-1}(\sqrt{-\alpha^2\Delta+1}-1)$04 regime; its nuclear form resolves shell evolution and pseudo-spin structure through explicit Hartree and Fock channels; and its molecular form shows how four-component and two-component relativistic Hamiltonians can be unified without abandoning Hartree-Fock exchange [(Dall'Acqua et al., 2010); (Li et al., 2019); (Liu, 2023)].