Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Progress of Theoretical Physics

The Hartree-Fock State and Fermion Path-Integral

Jun Ishida

Published Feb 1, 1981• 2 citations
Publisher

Abstract

Following the path-integral formalism by Ohnuki and Kashiwa, we calculate the transition amplitude exactly for a model which has four-fermi type interaction. We also investigate the time development of the expectation value of field operator to see to what extent the classical correspondence works in quantum mechanics of fermions.
DOI
10.1143/ptp.65.717

Related papers

Path integral approach for many-fermion systems and an extension of the time-dependent Hartree-Fock method

Kazunari Kaneko

The quantization of an extended time-dependent Hartree-Fock equation is performed by using the path integral formalism. The path integral form of time evolution operator between initial and final states is obtained with the use of the completeness relation of the fermion coherent states. It is shown that the extended time-dependent Hartree-Fock equation is naturally derived as a classical limit of the path integral.

Physical Review C•Aug 1, 1984•6 citations

The Feynman Path Integral

N. Cabibbo, Luciano Maiani, Omar Benhar A

This Chapter introduces the Feynman Path Integral description of quantum transition amplitudes, considering the simplest case of a quantum system described by one quantum variable q and its conjugate momentum p . The transition amplitude from time https://www.w3.org/1998/Math/MathML" display="inline"> t 1 to time https://www.w3.org/1998/Math/MathML" display="inline"> t 2 = t 1 + T is obtained by discretising the time axis in N time steps of amplitude https://www.w3.org/1998/Math/MathML" display="inline"> ϵ , implying T=N https://www.w3.org/1998/Math/MathML" display="inline"> ϵ . For https://www.w3.org/1998/Math/MathML" display="inline"> ϵ → 0 , the Feynman Path Integral formula is derived: https://www.w3.org/1998/Math/MathML" display="inline"> < q N | e - i T H / ℏ | q 0 > = ∫ d [ q ( t ) ] e i S q ( t ) , q ' ( t ) / ℏ , where integration is over all trajectories https://www.w3.org/1998/Math/MathML" display="inline"> q ( t ) which satisfy the boundary conditions at https://www.w3.org/1998/Math/MathML" display="inline"> t 1 and t 2 . https://www.w3.org/1998/Math/MathML" display="inline"> S ( q , q ' ) is the classical Action, a function of https://www.w3.org/1998/Math/MathML" display="inline"> q and of its time derivative. The classical limit, corresponding to https://www.w3.org/1998/Math/MathML" display="inline"> ℏ → 0 , is considered, whereby the solution of the Feynman Path Integral reduces to the classical trajectory, https://www.w3.org/1998/Math/MathML" display="inline"> q c ( t ) , which minimises the classical Action. The lattice approximation , i.e. small, non vanishing https://www.w3.org/1998/Math/MathML" display="inline"> ϵ , with evaluation of the Feynman integral by numerical integration, provides a useful approximation in the case of the fundamental theory of the strong interactions (Quantum Chromodynamics, or QCD), to be considered later; Chapters 17 to19. This chapter is completed by consideration of the extension of gthe formalism to imaginary time, where the Feynman Path Integral becomes the Partition Function of Statistical Mechanics, and by the definition of the Green's functions, ground-state expectation values of products of quantum operators q(t) taken at different times, https://www.w3.org/1998/Math/MathML" display="inline"> t 1 , t 2 … t N .

CRC Press eBooks•Oct 22, 2025•0 citations

Grassmann Variables

I. J. R. Aitchison, Anthony J.G. Hey

In the path integral representation of quantum amplitudes ( chapter 16 ) the fields are regarded as classical functions. Matrix elements of time-ordered products of bosonic operators could be satisfactorily represented (see the discussion following ( 16.79 )). But something new is needed to represent, for example, the time-ordered product of two fermionic operators: there must be a sign difference between the two orderings, since the fermionic operators anticommute . Thus it seems that to represent amplitudes involving fermionic operators by path integrals we must think in terms of ‘classical’ anticommuting variables.

CRC Press eBooks•Apr 19, 2024•0 citations

Fermionic mean-field dynamics for spin systems beyond free fermions

Rishab Dutta, Marc Illa, Niranjan Govind +1 more

We introduce the fermionized time-dependent Hartree-Fock (fTDHF), a real-time quantum dynamics method for spin-1/2 Hamiltonians following their mapping to fermions via the Jordan-Wigner transformation. fTDHF is formally equivalent to exact dynamics in the case of free fermions and can efficiently handle non-local string operators arising from long-range interactions via transition matrix elements between non-orthogonal Slater determinants. We show that the fTDHF method can be implemented on a classical computer with a cost that scales polynomially with system size, and linearly with the time steps. We benchmark fTDHF against exact dynamics on three separate spin-1/2 models, representing adiabatic preparation of states with long-range correlations, disorder-driven observation of many-body localization, and particle production in the Schwinger model. For each of these systems, fTDHF is shown to reproduce the qualitative dynamics generated by the exact evolutions, while maintaining a simple physical picture due to its mean-field nature.

cond-mat.str-elCross-listedarXiv:2604.02584•The Journal of Chemical Physics•Apr 6, 2026•1 citation

Feynman path integral for fermions

Robert E. Pugh

With Lagrangian quantum field theory as a starting point, the general quantum transition amplitude is derived as a sum of action exponentials over all paths connecting the initial and final fermion states. For bookkeeping convenience in intermediate stages anticommuting parameters are introduced. All integrals are determined by completeness and unitarity. Grassmann calculus is not postulated.

Physical review. D. Particles, fields, gravitation, and cosmology/Physical review. D. Particles and fields•May 15, 1988•3 citations

Many-body computations by stochastic sampling in Hartree-Fock-Bogoliubov space

Hao Shi, Shiwei Zhang

We describe the computational ingredients for an approach to treat interacting fermion systems with pairing fields, based on path integrals in the space of Hartree-Fock-Bogoliubov (HFB) wave functions. The path integrals can be evaluated by Monte Carlo, via random walks of HFB wave functions whose orbitals evolve stochastically. The approach combines the advantage of HFB theory in paired fermion systems and many-body quantum Monte Carlo techniques. The properties of HFB states, written in the form of either product states or Thouless states, are discussed. The states preserve forms when propagated by generalized one-body operators. They can be stabilized for numerical iteration. Overlaps and one-body Green's functions between two such states can be computed. A constrained-path or phaseless approximation can be applied to the random walks of the HFB states if a sign problem or phase problem is present. The method is illustrated with an exact numerical projection in the Kitaev model, and in the Hubbard model with attractive interaction under an external pairing field.

arXiv:1611.02689•Physical review. B./Physical review. B•Jan 27, 2017•28 citations
This Chapter introduces the Feynman Path Integral description of quantum transition amplitudes, considering the simplest case of a quantum system described by one quantum variable q and its conjugate momentum p . The transition amplitude from time https://www.w3.org/1998/Math/MathML" display="inline"> t 1 to time https://www.w3.org/1998/Math/MathML" display="inline"> t 2 = t 1 + T is obtained by discretising the time axis in N time steps of amplitude https://www.w3.org/1998/Math/MathML" display="inline"> ϵ , implying T=N https://www.w3.org/1998/Math/MathML" display="inline"> ϵ . For https://www.w3.org/1998/Math/MathML" display="inline"> ϵ → 0 , the Feynman Path Integral formula is derived: https://www.w3.org/1998/Math/MathML" display="inline"> < q N | e - i T H / ℏ | q 0 > = ∫ d [ q ( t ) ] e i S q ( t ) , q ' ( t ) / ℏ , where integration is over all trajectories https://www.w3.org/1998/Math/MathML" display="inline"> q ( t ) which satisfy the boundary conditions at https://www.w3.org/1998/Math/MathML" display="inline"> t 1 and t 2 . https://www.w3.org/1998/Math/MathML" display="inline"> S ( q , q ' ) is the classical Action, a function of https://www.w3.org/1998/Math/MathML" display="inline"> q and of its time derivative. The classical limit, corresponding to https://www.w3.org/1998/Math/MathML" display="inline"> ℏ → 0 , is considered, whereby the solution of the Feynman Path Integral reduces to the classical trajectory, https://www.w3.org/1998/Math/MathML" display="inline"> q c ( t ) , which minimises the classical Action. The lattice approximation , i.e. small, non vanishing https://www.w3.org/1998/Math/MathML" display="inline"> ϵ , with evaluation of the Feynman integral by numerical integration, provides a useful approximation in the case of the fundamental theory of the strong interactions (Quantum Chromodynamics, or QCD), to be considered later; Chapters 17 to19. This chapter is completed by consideration of the extension of gthe formalism to imaginary time, where the Feynman Path Integral becomes the Partition Function of Statistical Mechanics, and by the definition of the Green's functions, ground-state expectation values of products of quantum operators q(t) taken at different times, https://www.w3.org/1998/Math/MathML" display="inline"> t 1 , t 2 … t N .
In the path integral representation of quantum amplitudes ( chapter 16 ) the fields are regarded as classical functions. Matrix elements of time-ordered products of bosonic operators could be satisfactorily represented (see the discussion following ( 16.79 )). But something new is needed to represent, for example, the time-ordered product of two fermionic operators: there must be a sign difference between the two orderings, since the fermionic operators anticommute . Thus it seems that to represent amplitudes involving fermionic operators by path integrals we must think in terms of ‘classical’ anticommuting variables.
We introduce the fermionized time-dependent Hartree-Fock (fTDHF), a real-time quantum dynamics method for spin-1/2 Hamiltonians following their mapping to fermions via the Jordan-Wigner transformation. fTDHF is formally equivalent to exact dynamics in the case of free fermions and can efficiently handle non-local string operators arising from long-range interactions via transition matrix elements between non-orthogonal Slater determinants. We show that the fTDHF method can be implemented on a classical computer with a cost that scales polynomially with system size, and linearly with the time steps. We benchmark fTDHF against exact dynamics on three separate spin-1/2 models, representing adiabatic preparation of states with long-range correlations, disorder-driven observation of many-body localization, and particle production in the Schwinger model. For each of these systems, fTDHF is shown to reproduce the qualitative dynamics generated by the exact evolutions, while maintaining a simple physical picture due to its mean-field nature.
We describe the computational ingredients for an approach to treat interacting fermion systems with pairing fields, based on path integrals in the space of Hartree-Fock-Bogoliubov (HFB) wave functions. The path integrals can be evaluated by Monte Carlo, via random walks of HFB wave functions whose orbitals evolve stochastically. The approach combines the advantage of HFB theory in paired fermion systems and many-body quantum Monte Carlo techniques. The properties of HFB states, written in the form of either product states or Thouless states, are discussed. The states preserve forms when propagated by generalized one-body operators. They can be stabilized for numerical iteration. Overlaps and one-body Green's functions between two such states can be computed. A constrained-path or phaseless approximation can be applied to the random walks of the HFB states if a sign problem or phase problem is present. The method is illustrated with an exact numerical projection in the Kitaev model, and in the Hubbard model with attractive interaction under an external pairing field.