In quantum mechanics and quantum field theory, a Schrödinger field, named after Erwin Schrödinger, is a quantum field which obeys the Schrödinger equation.[1] While any situation described by a Schrödinger field can also be described by a many-body Schrödinger equation for identical particles, the field theory is more suitable for situations where the particle number changes.
A Schrödinger field is also the classical limit of a quantum Schrödinger field, a classical wave which satisfies the Schrödinger equation. Unlike the quantum mechanical wavefunction, if there are interactions between the particles the equation will be nonlinear. These nonlinear equations describe the classical wave limit of a system of interacting identical particles.
The path integral of a Schrödinger field is also known as a coherent state path integral, because the field itself is an annihilation operator whose eigenstates can be thought of as coherent states of the harmonic oscillations of the field modes.
Schrödinger fields are useful for describing Bose–Einstein condensation, the Bogolyubov - de Gennes equation of superconductivity, superfluidity, and many-body theory in general. They are also a useful alternative formalism for nonrelativistic quantum mechanics.
A Schrödinger field is the nonrelativistic limit of a Klein–Gordon field.
A Schrödinger field is a quantum field whose quanta obey the Schrödinger equation. In the classical limit, it can be understood as the quantized wave equation of a Bose Einstein condensate or a superfluid.
A Schrödinger field has the free field Lagrangian
L=\psi\dagger\left(i{\partial\over\partialt}+{\nabla2\over2m}\right)\psi.
When
\psi
\psi
If the particles interact with an external potential
V(x)
S=\intxt\psi\dagger\left(i{\partial\over\partialt}+{\nabla2\over2m}\right)\psi-\psi\dagger(x)\psi(x)V(x).
l{L}=i\psi\dagger\partialo\psi-
1 | |
2m |
(\partiali\psi\dagger\partiali\psi)-V\psi\dagger\psi
i\partialo\psi(x\mu)=\left(
-\Delta | |
2m |
+V(\vec{x})\right)\psi(x\mu)
-i\partialo\psi\dagger(x\mu)=\left(
-\Delta | |
2m |
+V(\vec{x})\right)\psi\dagger(x\mu)
\phii(x)
Ei
\psi(x)=\sumi\psii\phii(x).
S=\intt\sumi
\dagger\left( | |
\psi | |
i |
i{\partial\over\partialt}-Ei\right)\psii
To see the equivalence, note that decomposed into real and imaginary parts the action is:
S=\intt\sumi2\psir{d\psii\overdt}-Ei(\psi
2 | |
r |
+
2) | |
\psi | |
i |
\psir
S=\intt\sumi{1\overEi}\left(
d\psii | |
dt |
\right)2-Ei
2 | |
\psi | |
i |
Ei
V(x1,x2)
S=\intxt\psi\dagger\left(i
\partial | |
\partialt |
+{\nabla2\over2m}\right)\psi-\intxy\psi\dagger(y)\psi\dagger(x)V(x,y)\psi(x)\psi(y).
A pair-potential is the non-relativistic limit of a relativistic field coupled to electrodynamics. Ignoring the propagating degrees of freedom, the interaction between nonrelativistic electrons is the Coulomb repulsion. In 2+1 dimensions, this is:
V(x,y)={j2\over|y-x|}.
V(x1,x2)=λ\delta(x1-x2)
S=\intx\psi\dagger\left(i{\partial\over\partialt}+{\nabla2\over2m}\right)\psi+λ\intx\psi\dagger\psi\dagger\psi\psi
λ
The potentials can include many-body contributions. The interacting Lagrangian is then:
Li=\intx
\dagger(x | |
\psi | |
n) |
V(x1,x2,...,xn)\psi(x1)\psi(x2) … \psi(xn).
These types of potentials are important in some effective descriptions of close-packed atoms. Higher order interactions are less and less important.
The canonical momentum association with the field
\psi
\Pi(x)=i\psi\dagger.
The canonical commutation relations are like an independent harmonic oscillator at each point:
[\psi(x),\psi\dagger(y)]=\delta(x-y).
The field Hamiltonian is
H=S-\int\Pi(x){d\overdt}\psi=\int{|\nabla\psi|2\over2m}+\intxy\psi\dagger(x)\psi\dagger(y)V(x,y)\psi(x)\psi(y)
and the field equation for any interaction is a nonlinear and nonlocal version of the Schrödinger equation. For pairwise interactions:
i{\partial\over\partialt}\psi=-{\nabla2\over2m}\psi+\left(\intyV(x,y)\psi\dagger(y)\psi(y)\right)\psi(x).
The expansion in Feynman diagrams is called many-body perturbation theory. The propagator is
G(k)={1\overi\omega-{k2\over2m}}.
The interaction vertex is the Fourier transform of the pair-potential. In all the interactions, the number of incoming and outgoing lines is equal.
The many body Schrödinger equation for identical particles describes the time evolution of the many-body wavefunction ψ(x1, x2...xN) which is the probability amplitude for N particles to have the listed positions. The Schrödinger equation for ψ is:
i | \partial |
\partialt |
\psi=\left(
| |||||||
2m |
+
| |||||||
2m |
+ … +
| |||||||
2m |
+V(x1,x2,...,xN)\right)\psi
H=
| |||||||
2m |
+
| |||||||
2m |
+ … +
| |||||||
2m |
+V(x1,...,xN).
\psi(x1,x2,...)=\psi(x2,x1,...) forbosons
\psi(x1,x2,...)=-\psi(x2,x1,...) forfermions
Since the particles are indistinguishable, the potential V must be unchanged under permutations.If
V(x1,...,xN)=V1(x1)+V2(x2)+ … +VN(xN)
then it must be the case that
V1=V2= … =VN
V(x1...,xN)=V1,2(x1,x2)+V1,3(x2,x3)+V2,3(x1,x2)
then
V1,2=V1,3=V2,3
In the Schrödinger equation formalism, the restrictions on the potential are ad-hoc, and the classical wave limit is hard to reach. It also has limited usefulness if a system is open to the environment, because particles might coherently enter and leave.
A Schrödinger field is defined by extending the Hilbert space of states to include configurations with arbitrary particle number. A nearly complete basis for this set of states is the collection:
|N;x1,\ldots,xN\rangle
\psi0|0\rangle+\intx\psi1(x)|1;x\rangle+
\int | |
x1x2 |
\psi2(x1,x2)|2;x1x2\rangle+\ldots
\psi0
In order to reproduce the Schrödinger description, the inner product on the basis states should be
\langle1;x1|1;y1\rangle=\delta(x1-y1)
\langle2;x1x2|2;y1y2\rangle=\delta(x1-y1)\delta(x2-y2)\pm\delta(x1-y2)\delta(x2-y1)
There are natural operators in this Hilbert space. One operator, called , is the operator which introduces an extra particle at x. It is defined on each basis state:
\psi\dagger(x)\left|N;x1,...,xn\right\rangle=\left|N+1;x1,...,xn,x\right\rangle
Another operator removes a particle at x, and is called
\psi
\psi\dagger
\psi
\psi(x)\left|N;x1,...,xN\right\rangle=\delta(x-x1)\left|N-1;x2,...,xN\right\rangle+\delta(x-x2)\left|N-1;x1,x3,...,xN\right\rangle+ …
The position basis is an inconvenient way to understand coincident particles because states with a particle localized at one point have infinite energy, so intuition is difficult. In order to see what happens when two particles are at exactly the same point, it is mathematically simplest either to make space into a discrete lattice, or to Fourier transform the field in a finite volume.
The operator
\psi\dagger(k)=\intxe-ikx\psi\dagger(x)
\psi(k)=\intxeikx\psi(x)
If the potential energy for interaction of infinitely distant particles vanishes, the Fourier transformed operators in infinite volume create states which are noninteracting. The states are infinitely spread out, and the chance that the particles are nearby is zero.
The matrix elements for the operators between non-coincident points reconstructs the matrix elements of the Fourier transform between all modes:
\psi\dagger(k)\psi\dagger(k')-\psi\dagger(k')\psi\dagger(k)=0
\psi(k)\psi(k')-\psi(k')\psi(k)=0
\psi(k)\psi\dagger(k')-\psi(k')\psi\dagger(k)=\delta(k-k')
The commutation relations now determine the operators completely, and when the spatial volume is finite, there are no conceptual hurdle to understand coinciding momenta because momenta are discrete. In a discrete momentum basis, the basis states are:
|n1,n2,...nk\rangle
\dagger(k)|...,n | |
\psi | |
k,\ldots\rangle |
=\sqrt{nk+1}|...,nk+1,\ldots\rangle
\psi(k)\left|...,nk,\ldots\right\rangle=\sqrt{nk}\left|...,nk-1,\ldots\right\rangle
So that the operator
\sumk\psi\dagger(k)\psi(k)=\intx\psi\dagger(x)\psi(x)
Now it is easy to see that the matrix elements of and have harmonic oscillator commutation relations too.
[\psi(x),\psi(y)]=[\psi\dagger(x),\psi\dagger(y)]=0
[\psi(x),\psi\dagger(y)]=\delta(x-y)
The operator which removes and replaces a particle, acts as a sensor to detect if a particle is present at x. The operator acts to multiply the state by the gradient of the many body wavefunction. The operator
H=-\intx\psi\dagger(x){\nabla2\over2m}\psi(x)
\psi\daggeri{d\overdt}\psi=\psi\dagger{-\nabla2\over2m}\psi
i{\partial\over\partialt}\psi={-\nabla2\over2m}\psi
To add interactions, add nonlinear terms in the field equations. The field form automatically ensures that the potentials obey the restrictions from symmetry.
The field Hamiltonian which reproduces the equations of motion is
H={\nabla\psi\dagger\nabla\psi\over2m}
The Heisenberg equations of motion for this operator reproduces the equation of motion for the field.
To find the classical field Lagrangian, apply a Legendre transform to the classical limit of the Hamiltonian.
L=\psi\dagger\left(i{\partial\over\partialt}+{\nabla2\over2m}\right)\psi
Although this is correct classically, the quantum mechanical transformation is not completely conceptually straightforward because the path integral is over eigenvalues of operators ψ which are not hermitian and whose eigenvectors are not orthogonal. The path integral over field states therefore seems naively to be overcounting. This is not the case, because the time derivative term in L includes the overlap between the different field states.
The non-relativistic limit as
c\toinfty
\hat{a}p,\hat{b}p
\hat{a}(x)=\intd\Omegap\hat{a}
-ip ⋅ x | |
pe |
, \hat{b}(x)=\intd\Omegap\hat{b}
-ip ⋅ x | |
pe |
such that
\hat\phi(x)=\hata(x)+\hatb\dagger(x)
\hat{A}(x)
\hat{B}(x)
\hat{a}(x)= |
| ||||
\sqrt{2mc2 |
which factor out a rapidly oscillating phase due to the rest mass plus a vestige of the relativistic measure, the Lagrangian density
L=(\hbar
\mu\phi | |
c) | |
\mu\phi\partial |
\dagger-(mc2)2\phi\phi\dagger
\begin{align} L &=\left(\hbarc\right)2
\mu\hat{a} | |
\left(\partial | |
\mu\hat{a}\partial |
\dagger+
\mu\hat{b} | |
\partial | |
\mu\hat{b}\partial |
\dagger+ … \right)-\left(mc2\right)2\left(\hat{a}\hat{a}\dagger+\hat{b}\hat{b}\dagger+ … \right)\\ &=
1 | |
2mc2 |
\left[\left(\hbar
| ||||
c\right) |
\hat{A}+\partial0\hat{A}\right)\left(
imc | |
\hbar |
\hat{A}\dagger+\partial0\hat{A}\dagger\right)-\left(\hbar
x\hat{A} | |
c\right) | |
x\hat{A}\partial |
\dagger+(A ⇒ B)+ … -\left(mc2\right)2\left(\hat{A}\hat{A}\dagger+\hat{B}\hat{B}\dagger+ … \right)\right]\\ &=
\hbar2 | |
2m |
\left[
imc | |
\hbar |
\dagger-\hat{A}\partial | |
\left(\partial | |
0\hat{A}\hat{A} |
0\hat{A}\dagger\right)+ \partial\mu\hat{A}\partial\mu\hat{A}\dagger+(A ⇒ B)+ … \right] \end{align}
where terms proportional to
\pm2imc2t/\hbar | |
e |
\begin{align} LA&=i\hbar\hat{A}\dagger\hat{A}'+
\hbar2 | \left[ | |
2m |
1 | |
c2 |
\hat{A}'{\hat{A}'}\dagger-
x\hat{A} | |
\partial | |
x\hat{A}\partial |
\dagger\right]\\ &=i\hbar\hat{A}\dagger\hat{A}'+
\hbar2 | |
2m |
\left[
x\hat{A} | |
-\left(\partial | |
x\left(\hat{A}\partial |
\dagger\right)-
x\hat{A} | |
\hat{A}\partial | |
x\partial |
\dagger\right)\right]\\ &=i\hbar\hat{A}\dagger\hat{A}'+
\hbar2 | |
2m |
x\hat{A} | |
\hat{A}\partial | |
x\partial |
\dagger. \end{align}
The final Lagrangian takes the form[3]
L=
1 | |
2 |
\left[ \hat{A}\dagger\left(i\hbar
\partial | |
\partialt |
+
\hbar2\nabla2 | |
2m |
\right)\hat{A} +\hat{B}\dagger\left(i\hbar
\partial | |
\partialt |
+
\hbar2\nabla2 | |
2m |
\right)\hat{B} +h.c. \right].
\pm2imc2t/\hbar | |
e |
Compare for example to
F(x)=\sin(x)+\sin(10x)/10