This article is to be merged with Perturbation theory (quantum mechanics).
Quantum Mechanical Theory[edit]
Hamiltonian and Force Operator[edit]
From the differential geometric point of view, a parameterized Hamiltonian is considered as a function defined on the parameter manifold that maps each particular set of parameters
to an Hermitian operator
that acts on the Hilbert space. The parameters here can be external field, interaction strength, or driving parameters in the quantum phase transition. Let
and
be the nth eigenenergy and eigenstate of
respectively. In the language of deferential geometry, the states
form a vector bundle over the parameter manifold, on which derivatives of these states can be defined. The perturbation theory is to answer the following question: given
and
at a reference point
, how to estimate the
and
at
close to that reference point.
Without loss of generality, the coordinate system can be shifted, such that the reference point
is set to be the origin. The following linearized model Hamiltonian is frequently used
![{\displaystyle H(x^{\mu })=H(0)+x^{\mu }F_{\mu }.}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/1dcb9512aae63a8989e095069370860ae62df314)
If the parameters
are considered as generalized coordinates, then
should be identified as the generalized force operators related to those coordinates. Different indices
label the different forces along different directions in the parameter manifold. For example, if
denotes the external magnetic field in the
-direction, then
should be the magnetization in the same direction.
Perturbation as Power Series Expansion[edit]
The validity of the perturbation theory lies on the adiabatic assumption, which assumes the eigenenergies and eigenstates of the Hamiltonian are smooth functions of parameters such that their values in the vicinity region can be calculated in power series (like Taylor expansion) of the parameters:
![{\displaystyle E_{n}(x^{\mu })=E_{n}+x^{\mu }\partial _{\mu }E_{n}+{\frac {1}{2!}}x^{\mu }x^{\nu }\partial _{\mu }\partial _{\nu }E_{n}+\cdots ,}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/64ebfb04af68ae53643dff8d324e6a8e5bde8767)
![{\displaystyle |n(x^{\mu })\rangle =|n\rangle +x^{\mu }|\partial _{\mu }n\rangle +{\frac {1}{2!}}x^{\mu }x^{\nu }|\partial _{\mu }\partial _{\nu }n\rangle +\cdots .}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/7a0be9d94139bc7c68a459c98b4f90ca14aa2c09)
Here
denotes the derivative with respect to
. When applying to the state
, it should be understood as the Lie derivative if the vector bundle is equipped with non-vanishing connection. All the terms on the right-hand-side of the series are evaluated at
, e.g.
and
. This convention will be adopted in though out this section, that all functions without the parameter dependence explicitly stated are assumed to be evaluated at the origin. The power series may converge slowly or even not converging when the energy levels are close to each other. The adiabatic assumption breaks down when there is energy level degeneracy, and hence the perturbation theory is not applicable in that case.
Hellmann-Feynman Theorems[edit]
The above power series expansion can be readily evaluated if there is a systematic approach to calculate the derivates to any order. Using the chain rule, the derivatives can be broken down to the single derivative on either the energy or the sate. The Hellmann-Feynman theorems are used to calculated these single derivatives. The first Hellmann-Feynman theorem gives the derivative of the energy,
![{\displaystyle \partial _{\mu }E_{n}=\langle n|\partial _{\mu }H|n\rangle }](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/8a18caca6250298e16f827e7bb07633319e7ee7e)
The second Hellmann-Feynman theorem gives the derivative of the state (resolved by the complete basis with
),
![{\displaystyle \langle m|\partial _{\mu }n\rangle ={\frac {\langle m|\partial _{\mu }H|n\rangle }{E_{n}-E_{m}}}}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/96e4f40daf63564f616b57fe31a5e965e58417e7)
![{\displaystyle \langle \partial _{\mu }m|n\rangle ={\frac {\langle m|\partial _{\mu }H|n\rangle }{E_{m}-E_{n}}}}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/1ecac03d5c1a3064ce3b1fae21e6e5678fbf147a)
The theorems can be simply derived by applying the differential operator
to both sides of the Schrödinger equation
, which reads
![{\displaystyle \partial _{\mu }H|n\rangle +H|\partial _{\mu }n\rangle =\partial _{\mu }E_{n}|n\rangle +E_{n}|\partial _{\mu }n\rangle }](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/0a2663b6ab13f083c8bb23484285163fc6dd90e1)
Then overlap with the state
from left and make use of the Schrödinger equation again
,
![{\displaystyle \langle m|\partial _{\mu }H|n\rangle +E_{m}\langle m|\partial _{\mu }n\rangle =\partial _{\mu }E_{n}\langle m|n\rangle +E_{n}\langle m|\partial _{\mu }n\rangle }](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/73cbc230f3732eb716514a79ab0bdbc61ee9ec71)
Given that the states from a set of orthonormal basis
, by rearranging the terms in the above equation, it is straight forward to prove both the Hellmann-Feynman theorems. With the differential rules given by the Hellmann-Feynman theorems, the perturbative correction to the energies and states can be calculated systematically.
Correction of Energy[edit]
To the second order
![{\displaystyle E_{n}(x^{\mu })=\langle n|H|n\rangle +\langle n|\partial _{\mu }H|n\rangle x^{\mu }+\sum _{m\neq n}{\frac {\langle n|\partial _{\nu }H|m\rangle \langle m|\partial _{\mu }H|n\rangle }{E_{n}-E_{m}}}x^{\mu }x^{\nu }+\cdots }](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/940ae94e101b3ae149238b7ce6eed6856eb395e4)
Correction of State[edit]
To the second order
$$
|n(x^\mu)\rangle =|n\rangle +\sum _{m\neq n} \frac{\langle m|\partial_\mu H|n\rangle }{E_n-E_m}|m\rangle x^\mu$$
![{\displaystyle +\left(\sum _{m\neq n}\sum _{l\neq n}{\frac {\langle m|\partial _{\mu }H|l\rangle \langle l|\partial _{\nu }H|n\rangle }{(E_{n}-E_{m})(E_{n}-E_{l})}}|m\rangle -\sum _{m\neq n}{\frac {\langle m|\partial _{\mu }H|n\rangle \langle n|\partial _{\nu }H|n\rangle }{(E_{n}-E_{m})^{2}}}|m\rangle -{\frac {1}{2}}\sum _{m\neq n}{\frac {\langle n|\partial _{\mu }H|m\rangle \langle m|\partial _{\nu }H|n\rangle }{(E_{n}-E_{m})^{2}}}|n\rangle \right)x^{\mu }x^{\nu }+\cdots }](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/51581d1c6d7cc59daca562fe7b06445d9233b2cb)
Effective Hamiltonian[edit]
Let
be the Hamiltonian completely restricted either in the low-energy subspace
or in the high-energy subspace
, such that there is no matrix element in
connecting the low- and the high-energy subspaces, i.e.
if
. Let
be the coupling terms connecting the subspaces. Then when the high energy degrees of freedoms are integrated out, the effective Hamiltonian in the low energy subspace reads
![{\displaystyle H_{mn}^{\text{eff}}\left(x^{\mu }\right)=\langle m|H|n\rangle +\langle m|\partial _{\mu }H|n\rangle x^{\mu }+{\frac {1}{2!}}\sum _{l\in {\mathcal {H}}_{H}}\left({\frac {\langle m|\partial _{\mu }H|l\rangle \langle l|\partial _{\nu }H|n\rangle }{E_{m}-E_{l}}}+{\frac {\langle m|\partial _{\nu }H|l\rangle \langle l|\partial _{\mu }H|n\rangle }{E_{n}-E_{l}}}\right)x^{\mu }x^{\nu }+\cdots }](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/0fdbcb859e462db4c8b66a424bd5e71efb939642)
Here
are restricted in the low energy subspace. The above result can be derived by power series expansion of
.
Quantum Field Theory[edit]
For finite temperature case, we calculate the correction to the free energy instead of the ground state energy.
Kernel, Free Energy and Propagator[edit]
In quantum field theory, instead of Hamiltonian, an action is used to describe a model. In general, the action of a free field
takes the bilinear form,
![{\displaystyle S=\phi ^{\dagger }\cdot K\cdot \phi }](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/f602a6be73604727c200d9c8ec26b267d08c15bd)
which involves the convolution with the action kernel
. The action kernel can be
for diffusive dynamics, or
for wave dynamics. In any case, the free energy is given by
![{\displaystyle F=\eta \mathrm {Tr} \ln K}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/08f82b21a250a4e50f62f95a295574fe2e061fe5)
where
for bosonic field and
for fermionic field. The perturbation expansion of
can be calculated in Taylor series term by term as long as we know how to take partial derivatives of the action kernel. Because the free energy involves the logarithm of
, so after first order derivative,
will appear in the formula. Therefore the propagator
![{\displaystyle G=-K^{-1}}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/cb25fb6bf6fdcbd9e7c63416b5505856c874bcc3)
is introduced to denote the inverse of the action kernel.
The free energy and the propagator in the quantum field theory take the role of the energy and the state in quantum mechanics respectively. So the goal of perturbation theory at the field theory level is to calculate the free energy and the propagator in power series, given the action kernel of the model.
The Derivative Formula[edit]
Each term in the Taylor expansion involves evaluation of derivatives. The derivative formula for the propagator reads
![{\displaystyle \partial _{\mu }G=G\cdot \partial _{\mu }K\cdot G}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/c4f9a89342bf1bea04f4ac77288703681e84488f)
This can be derived by starting from the identity
, and take partial derivative on both sides
![{\displaystyle \partial _{\mu }G\cdot K+G\cdot \partial _{\mu }K=0}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/e6b60cced5f6515c00978a519bae06b34cbbc6a4)
then the derivative formula is easily obtained by rearranging the terms.
The derivative of free energy to any order can be calculated by recursively using the derivative formula for the propagator. For example,
![{\displaystyle \partial _{\mu }F=-\eta \mathrm {Tr} \,G\cdot \partial _{\mu }K}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/85045788d00c4e318655fa972016df0951f33a4c)
![{\displaystyle \partial _{\mu }\partial _{\nu }F=-\eta \mathrm {Tr} (G\cdot \partial _{\mu }\partial _{\nu }K+G\cdot \partial _{\mu }K\cdot G\cdot \partial _{\nu }K)}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/3572e0f44b573b969f864999cb42a359ca9d1be7)
The expression grows complicated with the order quickly. For single variable, if the derivative of the Kernel is cut-off at the second order, i.e.
, then the derivative of free energy has the following simple form
![{\displaystyle \partial ^{n}F=-\eta \,(n-1)!\,\mathrm {Tr} (G\cdot \partial K)^{n}}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/8b51b188d63aa4cbb3bc737f529593239492ee5d)
Free Energy Correction[edit]
Then the correction of free energy for the single parameter perturbation can be obtained by plugging the derivative formula into the Taylor series,
![{\displaystyle F=F_{0}-\eta \sum _{n=1}^{\infty }{\frac {1}{n}}\mathrm {Tr} (G_{0}\cdot \partial K)^{n}x^{n}}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/57067eaae3ab2334461d902cc3ed6963f7e0f573)
and
are respectively the free energy and the propagator evaluated at the origin
of the perturbation parameter. The series may be formally expressed as a logarithm
![{\displaystyle F=F_{0}+\eta \mathrm {Tr} \ln(1-(G_{0}\cdot \partial K)x)}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/20d2b5e7a949f18587c741e210ca95bea37de65a)
Dynson's Equation and Self-Energy Correction[edit]
Notice the relation between the free energy and the propagator
. The above free energy correction implies
![{\displaystyle \eta \mathrm {Tr} \ln(-G^{-1})=\eta \mathrm {Tr} \ln(-G_{0}^{-1})+\eta \mathrm {Tr} \ln(1-(G_{0}\cdot \partial K)x)=\eta \mathrm {Tr} \ln(-G_{0}^{-1}+x\partial K))}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/c7069111dc2120f1181cd85a04147d90bc448d93)
which is the Dynson's equation of the propagator,
![{\displaystyle G^{-1}=G_{0}^{-1}-x\partial K}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/a23872ddc7104e09f6cb042dc969fc4ffeef5c6b)
One immediately identifies
to be the bare self-energy, which becomes obvious from the fact that the Dynson's equation can be also obtained from the series expansion of
as
, with
and
, while
carries the energy (action) dimension.
The Dynson's equation may be formally solved as
![{\displaystyle G=(G_{0}^{-1}-\Sigma _{0})^{-1}=(1-G_{0}\cdot \Sigma _{0})^{-1}G_{0}}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/b1d731168852beff70da6116aae014fe7c120e1b)
however the difficulty lies in the evaluation of the matrix inverse. A way to circumvent the inversion is to use perturbative expansion, such that
![{\displaystyle G=G_{0}+G_{0}\cdot \Sigma _{0}\cdot G_{0}+G_{0}\cdot \Sigma _{0}\cdot G_{0}\cdot \Sigma _{0}\cdot G_{0}+\cdots =G_{0}\cdot \sum _{n=1}^{\infty }(\Sigma _{0}\cdot G_{0})^{n}}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/022879afe01cf75634a994e859580d3f5efb0ec0)
This leads to the definition of the self-energy, which follows form
![{\displaystyle G=G_{0}+G_{0}\cdot \Sigma \cdot G_{0}}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/e5a4c44c375202ed4555183f0679427e46400e19)
meaning that the normal propagator is dressed from the bare propagator by the self-energy correction.
Effective Hamiltonian under Perturbation[edit]
To obtain an effective Hamiltonian under perturbation, we can calculate the effective propagator by series expansion, then restore the Hamiltonian by
.....
![{\displaystyle \Sigma =\Sigma _{0}+\Sigma _{0}\cdot G_{0}\cdot \Sigma _{0}+\Sigma _{0}\cdot G_{0}\cdot \Sigma _{0}\cdot G_{0}\cdot \Sigma _{0}+\ldots =\Sigma _{0}\cdot \sum _{n=1}^{\infty }(G_{0}\cdot \Sigma _{0})^{n}.}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/bb08fed21e659072cff2c4934a4efa3c7179fab6)
The bare self-energy
stands for the perturbation (connecting low and high energy sectors), and
contains the dynamics of both energy sectors independently. The resulting self-energy can then be projected to the low-energy subspace to produce the effective Hamiltonian.
The normal self-energy obtained by the above formula is a function of frequency
. To obtain the effective correction near energy
, the trick is to evaluate the following residue
![{\displaystyle H_{\text{eff}}(\omega _{0})=\mathrm {Res} _{\omega \rightarrow \omega _{0}}{\frac {\Sigma (\omega )}{\omega -\omega _{0}}}.}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/bf8ee41db526f4a8ea325d6a792bae9720a1ecbc)
This is a way to regularize
around
by taking the 0th order term in its Laurent series.
However, when the low-energy Hamiltonian is not trivially zero, we need to set
to the low-energy Hamiltonian
itself, and the effective Hamiltonian reads
![{\displaystyle H_{\text{eff}}=H_{0}+\Sigma (H_{0}).}](https://wikimedia.riteme.site/api/rest_v1/media/math/render/svg/e8987a26ef8522190de04cf2f8d85aee5f54e203)