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

Hamiltonian Flow Equations for a Dirac Particle in an External Potential

Alexander B. Bylev, Hans J. Pirner

Published Jan 1, 1998• 12 citations
OpenAlex

Abstract

We derive and solve the Hamiltonian flow equations for a Dirac particle in an external static potential. The method shows a general procedure for the set up of continuous unitary transformations to reduce the Hamiltonian to a quasidiagonal form. Supported by DFG Grant 436 RUS 113/205/1 on Russian German CooperationThe well known Foldy-Wouthuysen (FW) transformation [1] is a unitary transformation that separates big and small components in the Dirac equation for a particle interacting with an external electromagnetic field. If the external field is stationary then the FW transformation reduces to a unitary transformation of the Hamiltonian H of the Dirac particle H ′ = e iS (⃗α(⃗p − e ⃗ A(⃗x)) + βm + eA0(⃗x))e −iS (1) where S can be expanded in a sequence of hermitean operators in powers of 1/m. Recently a novel method to diagonalize Hamiltonians via continuous unitary transformations has been proposed [2], where the Hamiltonian H and the S are considered to be functions of a continuous variable, the flow parameter l. This new method has been applied to several topics like the n − orbital model [2], superconductivity and the electron phonon interaction [3, 4]. We want to apply it here to the Dirac particle in an external field in order to explore its possibilities with a well known problem. Thereby we can gain experience in solving problems like the elimination of negative energy states in near light cone QCD [5], where the solution has yet to be discovered. For infinitesimal steps exp(η(l)dl) a continuous unitary transformation of the Hamiltonian H leads to the so-called Hamiltonian flow equation dH(l) dl = [η(l), H(l)] (2) with initial condition H(0) = H. The transformation matrix η(l) is antihermitian and related to S as follows η(l) = ( d dl eiS(l))e −iS(l). The parameter l has dimension of energy −2 and characterizes the amount of diagonalization which has been achieved. In the energy representation the matrix elements of the Hamiltonian H will be grouped in a band of width

Related papers

A Flow Equation Approach Striving Towards an Energy-Separating Hamiltonian Unitary Equivalent to the Dirac Hamiltonian with Coupling to Electromagnetic Fields

Nils Schopohl, Nadine Suzan Cetin

The Dirac Hamiltonian $H^{\left(D\right)}$ for relativistic charged fermions minimally coupled to (possibly time-dependent) electromagnetic fields is transformed with a purpose-built flow equation method, so that the result of that transformation is unitary equivalent to $H^{\left(D\right)}$ and granted to strive towards a limiting value $H^{\left(NW\right)}$ commuting with the Dirac $β$-matrix. Upon expansion of $H^{\left(NW\right)}$ to order $\frac{v^2}{c^2}$ the nonrelativistic Hamiltonian $H^{\left(SP\right)}$ of Schrödinger-Pauli quantum mechanics emerges as the leading order term adding to the rest energy $mc^2$. All the relativistic corrections to $H^{\left(SP\right)}$ are explicitly taken into account in the guise of a Magnus type series expansion, the series coefficients generated to order $\left(\frac{v^{2}}{c^{2}}\right)^{n}$ for $n\geq2$ comprising partial sums of iterated commutators only. In the special case of static fields the equivalence of the flow equation method with the well known energy-separating unitary transformation of Eriksen is established on the basis of an exact solution of a reverse flow equation transforming the $β$-matrix into the energy-sign operator associated with $H^{\left(D\right)}$. That way the identity $H^{\left(NW\right)}=β\sqrt{H^{\left(NW\right)}H^{\left(NW\right)}}$ is established implying $H^{\left(NW\right)}$ being determined unambiguously.

arXiv:2207.12825•Jul 26, 2022•0 citations

Brownian path integral from Dirac equation: a probabilistic approach to the Foldy-Wouthuysen transformation

Gian Fabrizio De Angelis, Maurizio Serva

The spectrum of the Dirac Hamiltonian H D has a positive and a negative part, the first corresponds to the energy levels of the electron, the second to minus the energy levels of the positron. The Foldy-Wouthuysen transformation is a tool to obtain a decomposition of H D in a direct sum H D =H e- (+)-H e+ , where H e- is the Hamiltonian of the electron and H e+ is that of the positron. Unfortunately, this decomposition has been obtained in a closed form only when the external field is purely magnetic, while, in the presence of an electric field, only a perturbative expansion is available. The authors give a path integral representation of the semigroup phi (x, t)=exp(-(t/)H e- ) phi 0 (x) in an external electromagnetic field. Their formula is the relativistic version, for Dirac particles, of the well known Feynman-Kac-Ito formula.

Journal of Physics A Mathematical and General•Dec 21, 1992•2 citations

The Necessary and Sufficient Conditions for Transformation from Dirac Representation to Foldy-Wouthuysen Representation

В. П. Незнамов

The paper describes conditions for transformation from the Dirac representation to the Foldy-Wouthuysen representation. The necessary condition is the block-diagonal transformation of Hamiltonian relative to the upper and lower components of the wave function. The sufficient condition is the wave function transformation law described by relation (6). It has been demonstrated that the unitary transformations offered by the authors of the papers [14], [15], [16], [17] do not satisfy the sufficiency condition (6) and, hence, they are not the Foldy-Wouthuysen transformations. In applications, the matrix elements of any operator in the FW representation can be calculated, according to (6), using the normalized two-component wave functions in the Dirac representation known for the given problem.

arXiv:0804.0333•Apr 2, 2008•3 citations

Exact Foldy-Wouthuysen transformation of the Dirac-Pauli Hamiltonian by Kutzelnigg's method

Dah-Wei Chiou, Tsung-Wei Chen

We apply Kutzelnigg's method for the Foldy-Wouthuysen (FW) transformation of the Dirac-Pauli Hamiltonian. The exact FW transformations exist and agree with those obtained by Eriksen's method for two special cases. In the weak-field limit of static and homogeneous electromagnetic fields, the long-held speculation is rigorously proven, by mathematical induction on the orders of $1/c$ in the power series, that the FW transformed Dirac-Pauli Hamiltonian is in full agreement with the classical counterpart, which is the sum of the orbital Hamiltonian for the Lorentz force equation and the spin Hamiltonian for the Thomas-Bargmann-Michel-Telegdi equation.

•May 18, 2014•0 citations

Unitary transformation for new diagonalization of Dirac Hamiltonian

Hong-Yi Fan, Ze Wu

The Foldy–Wouthuysen transformation was originally developed to elucidate the nonrelativistic limit of the Dirac equation. In this paper, we introduce a novel unitary transformation for the Hamiltonian of the Dirac equation. The unitary transformation operator is given by [Formula: see text] which results in the diagonal Hamiltonians [Formula: see text] and [Formula: see text]. We also apply this new unitary transformation to wave function calculations, and elaborate on its significance.

Modern Physics Letters A•Oct 31, 2025•0 citations

Separation of Dirac's Hamiltonian by Van Vleck transformation

Flemming Jørgensen

The now classic Foldy–Wouthuysen transformation (FWT) was introduced as successive unitary transformations. This fundamental idea has become the standard in later developments such as the Douglas–Kroll transformation (DKT) – but it is not the only possibility. FWT can be seen as a simple special case of the general Van Vleck transformation (VVT) which besides the successive version has another, known as the canonical because of a series of nice mathematical properties discovered gradually over time. The aim of the present paper is to compare the two approaches – which give identical results in the lower orders, but not in the higher. After having recapitalised both, we apply them to Dirac's Hamiltonian for the electron in a constant electromagnetic field, written with so few assumptions about the operators that the mathematical techniques stand out separated from the terminology of relativistic quantum mechanics. FWT for a free particle is dealt with by a recent geometric approach to VVT. The original FWT is continued through the next non-zero orders. DKT is considered with special weight on equivalent formulations of the generalised and the optimised forms introduced by Wolf, Reiher and Hess.

Molecular Physics•Oct 12, 2016•2 citations
We derive and solve the Hamiltonian flow equations for a Dirac particle in an external static potential. The method shows a general procedure for the set up of continuous unitary transformations to reduce the Hamiltonian to a quasidiagonal form. Supported by DFG Grant 436 RUS 113/205/1 on Russian German CooperationThe well known Foldy-Wouthuysen (FW) transformation [1] is a unitary transformation that separates big and small components in the Dirac equation for a particle interacting with an external electromagnetic field. If the external field is stationary then the FW transformation reduces to a unitary transformation of the Hamiltonian H of the Dirac particle H ′ = e iS (⃗α(⃗p − e ⃗ A(⃗x)) + βm + eA0(⃗x))e −iS (1) where S can be expanded in a sequence of hermitean operators in powers of 1/m. Recently a novel method to diagonalize Hamiltonians via continuous unitary transformations has been proposed [2], where the Hamiltonian H and the S are considered to be functions of a continuous variable, the flow parameter l. This new method has been applied to several topics like the n − orbital model [2], superconductivity and the electron phonon interaction [3, 4]. We want to apply it here to the Dirac particle in an external field in order to explore its possibilities with a well known problem. Thereby we can gain experience in solving problems like the elimination of negative energy states in near light cone QCD [5], where the solution has yet to be discovered. For infinitesimal steps exp(η(l)dl) a continuous unitary transformation of the Hamiltonian H leads to the so-called Hamiltonian flow equation dH(l) dl = [η(l), H(l)] (2) with initial condition H(0) = H. The transformation matrix η(l) is antihermitian and related to S as follows η(l) = ( d dl eiS(l))e −iS(l). The parameter l has dimension of energy −2 and characterizes the amount of diagonalization which has been achieved. In the energy representation the matrix elements of the Hamiltonian H will be grouped in a band of width
The Dirac Hamiltonian H(D)H^{\left(D\right)}H(D) for relativistic charged fermions minimally coupled to (possibly time-dependent) electromagnetic fields is transformed with a purpose-built flow equation method, so that the result of that transformation is unitary equivalent to H(D)H^{\left(D\right)}H(D) and granted to strive towards a limiting value H(NW)H^{\left(NW\right)}H(NW) commuting with the Dirac βββ-matrix. Upon expansion of H(NW)H^{\left(NW\right)}H(NW) to order v2c2\frac{v^2}{c^2}c2v2​ the nonrelativistic Hamiltonian H(SP)H^{\left(SP\right)}H(SP) of Schrödinger-Pauli quantum mechanics emerges as the leading order term adding to the rest energy mc2mc^2mc2. All the relativistic corrections to H(SP)H^{\left(SP\right)}H(SP) are explicitly taken into account in the guise of a Magnus type series expansion, the series coefficients generated to order (v2c2)n\left(\frac{v^{2}}{c^{2}}\right)^{n}(c2v2​)n for n≥2n\geq2n≥2 comprising partial sums of iterated commutators only. In the special case of static fields the equivalence of the flow equation method with the well known energy-separating unitary transformation of Eriksen is established on the basis of an exact solution of a reverse flow equation transforming the βββ-matrix into the energy-sign operator associated with H(D)H^{\left(D\right)}H(D). That way the identity H(NW)=βH(NW)H(NW)H^{\left(NW\right)}=β\sqrt{H^{\left(NW\right)}H^{\left(NW\right)}}H(NW)=βH(NW)H(NW)​ is established implying H(NW)H^{\left(NW\right)}H(NW) being determined unambiguously.
The spectrum of the Dirac Hamiltonian H D has a positive and a negative part, the first corresponds to the energy levels of the electron, the second to minus the energy levels of the positron. The Foldy-Wouthuysen transformation is a tool to obtain a decomposition of H D in a direct sum H D =H e- (+)-H e+ , where H e- is the Hamiltonian of the electron and H e+ is that of the positron. Unfortunately, this decomposition has been obtained in a closed form only when the external field is purely magnetic, while, in the presence of an electric field, only a perturbative expansion is available. The authors give a path integral representation of the semigroup phi (x, t)=exp(-(t/)H e- ) phi 0 (x) in an external electromagnetic field. Their formula is the relativistic version, for Dirac particles, of the well known Feynman-Kac-Ito formula.
The paper describes conditions for transformation from the Dirac representation to the Foldy-Wouthuysen representation. The necessary condition is the block-diagonal transformation of Hamiltonian relative to the upper and lower components of the wave function. The sufficient condition is the wave function transformation law described by relation (6). It has been demonstrated that the unitary transformations offered by the authors of the papers [14], [15], [16], [17] do not satisfy the sufficiency condition (6) and, hence, they are not the Foldy-Wouthuysen transformations. In applications, the matrix elements of any operator in the FW representation can be calculated, according to (6), using the normalized two-component wave functions in the Dirac representation known for the given problem.
We apply Kutzelnigg's method for the Foldy-Wouthuysen (FW) transformation of the Dirac-Pauli Hamiltonian. The exact FW transformations exist and agree with those obtained by Eriksen's method for two special cases. In the weak-field limit of static and homogeneous electromagnetic fields, the long-held speculation is rigorously proven, by mathematical induction on the orders of 1/c1/c1/c in the power series, that the FW transformed Dirac-Pauli Hamiltonian is in full agreement with the classical counterpart, which is the sum of the orbital Hamiltonian for the Lorentz force equation and the spin Hamiltonian for the Thomas-Bargmann-Michel-Telegdi equation.
The now classic Foldy–Wouthuysen transformation (FWT) was introduced as successive unitary transformations. This fundamental idea has become the standard in later developments such as the Douglas–Kroll transformation (DKT) – but it is not the only possibility. FWT can be seen as a simple special case of the general Van Vleck transformation (VVT) which besides the successive version has another, known as the canonical because of a series of nice mathematical properties discovered gradually over time. The aim of the present paper is to compare the two approaches – which give identical results in the lower orders, but not in the higher. After having recapitalised both, we apply them to Dirac's Hamiltonian for the electron in a constant electromagnetic field, written with so few assumptions about the operators that the mathematical techniques stand out separated from the terminology of relativistic quantum mechanics. FWT for a free particle is dealt with by a recent geometric approach to VVT. The original FWT is continued through the next non-zero orders. DKT is considered with special weight on equivalent formulations of the generalised and the optimised forms introduced by Wolf, Reiher and Hess.