Definition (non-relativistic 3-current)
edit
Free spin-0 particle
edit
In non-relativistic quantum mechanics, the probability current j of the wave function Ψ of a particle of mass m in one dimension is defined as[2]
j
=
ℏ
2
m
i
(
Ψ
∗
∂
Ψ
∂
x
−
Ψ
∂
Ψ
∗
∂
x
)
=
ℏ
m
ℜ
{
Ψ
∗
1
i
∂
Ψ
∂
x
}
=
ℏ
m
ℑ
{
Ψ
∗
∂
Ψ
∂
x
}
,
{\displaystyle j={\frac {\hbar }{2mi}}\left(\Psi ^{*}{\frac {\partial \Psi }{\partial x}}-\Psi {\frac {\partial \Psi ^{*}}{\partial x}}\right)={\frac {\hbar }{m}}\Re \left\{\Psi ^{*}{\frac {1}{i}}{\frac {\partial \Psi }{\partial x}}\right\}={\frac {\hbar }{m}}\Im \left\{\Psi ^{*}{\frac {\partial \Psi }{\partial x}}\right\},}
where
Note that the probability current is proportional to a Wronskian
W
(
Ψ
,
Ψ
∗
)
.
{\displaystyle W(\Psi ,\Psi ^{*}).}
In three dimensions, this generalizes to
j
=
ℏ
2
m
i
(
Ψ
∗
∇
Ψ
−
Ψ
∇
Ψ
∗
)
=
ℏ
m
ℜ
{
Ψ
∗
∇
i
Ψ
}
=
ℏ
m
ℑ
{
Ψ
∗
∇
Ψ
}
,
{\displaystyle \mathbf {j} ={\frac {\hbar }{2mi}}\left(\Psi ^{*}\mathbf {\nabla } \Psi -\Psi \mathbf {\nabla } \Psi ^{*}\right)={\frac {\hbar }{m}}\Re \left\{\Psi ^{*}{\frac {\nabla }{i}}\Psi \right\}={\frac {\hbar }{m}}\Im \left\{\Psi ^{*}\nabla \Psi \right\}\,,}
where
∇
{\displaystyle \nabla }
denotes the del or gradient operator . This can be simplified in terms of the kinetic momentum operator ,
p
^
=
−
i
ℏ
∇
{\displaystyle \mathbf {\hat {p}} =-i\hbar \nabla }
to obtain
j
=
1
2
m
(
Ψ
∗
p
^
Ψ
−
Ψ
p
^
Ψ
∗
)
.
{\displaystyle \mathbf {j} ={\frac {1}{2m}}\left(\Psi ^{*}\mathbf {\hat {p}} \Psi -\Psi \mathbf {\hat {p}} \Psi ^{*}\right)\,.}
These definitions use the position basis (i.e. for a wavefunction in position space ), but momentum space is possible.
Spin-0 particle in an electromagnetic field
edit
The above definition should be modified for a system in an external electromagnetic field . In SI units , a charged particle of mass m and electric charge q includes a term due to the interaction with the electromagnetic field;[3]
j
=
1
2
m
[
(
Ψ
∗
p
^
Ψ
−
Ψ
p
^
Ψ
∗
)
−
2
q
A
|
Ψ
|
2
]
{\displaystyle \mathbf {j} ={\frac {1}{2m}}\left[\left(\Psi ^{*}\mathbf {\hat {p}} \Psi -\Psi \mathbf {\hat {p}} \Psi ^{*}\right)-2q\mathbf {A} |\Psi |^{2}\right]}
where A = A (r , t ) is the magnetic vector potential . The term q A has dimensions of momentum. Note that
p
^
=
−
i
ℏ
∇
{\displaystyle \mathbf {\hat {p}} =-i\hbar \nabla }
used here is the canonical momentum and is not gauge invariant , unlike the kinetic momentum operator
P
^
=
−
i
ℏ
∇
−
q
A
{\displaystyle \mathbf {\hat {P}} =-i\hbar \nabla -q\mathbf {A} }
.
In Gaussian units :
j
=
1
2
m
[
(
Ψ
∗
p
^
Ψ
−
Ψ
p
^
Ψ
∗
)
−
2
q
c
A
|
Ψ
|
2
]
{\displaystyle \mathbf {j} ={\frac {1}{2m}}\left[\left(\Psi ^{*}\mathbf {\hat {p}} \Psi -\Psi \mathbf {\hat {p}} \Psi ^{*}\right)-2{\frac {q}{c}}\mathbf {A} |\Psi |^{2}\right]}
where c is the speed of light .
Spin-s particle in an electromagnetic field
edit
If the particle has spin , it has a corresponding magnetic moment , so an extra term needs to be added incorporating the spin interaction with the electromagnetic field.
According to Landau-Lifschitz's Course of Theoretical Physics the electric current density is in Gaussian units:[4]
j
e
=
q
2
m
[
(
Ψ
∗
p
^
Ψ
−
Ψ
p
^
Ψ
∗
)
−
2
q
c
A
|
Ψ
|
2
]
+
μ
S
c
s
ℏ
∇
×
(
Ψ
∗
S
Ψ
)
{\displaystyle \mathbf {j} _{e}={\frac {q}{2m}}\left[\left(\Psi ^{*}\mathbf {\hat {p}} \Psi -\Psi \mathbf {\hat {p}} \Psi ^{*}\right)-{\frac {2q}{c}}\mathbf {A} |\Psi |^{2}\right]+{\frac {\mu _{S}c}{s\hbar }}\nabla \times (\Psi ^{*}\mathbf {S} \Psi )}
And in SI units:
j
e
=
q
2
m
[
(
Ψ
∗
p
^
Ψ
−
Ψ
p
^
Ψ
∗
)
−
2
q
A
|
Ψ
|
2
]
+
μ
S
s
ℏ
∇
×
(
Ψ
∗
S
Ψ
)
{\displaystyle \mathbf {j} _{e}={\frac {q}{2m}}\left[\left(\Psi ^{*}\mathbf {\hat {p}} \Psi -\Psi \mathbf {\hat {p}} \Psi ^{*}\right)-2q\mathbf {A} |\Psi |^{2}\right]+{\frac {\mu _{S}}{s\hbar }}\nabla \times (\Psi ^{*}\mathbf {S} \Psi )}
Hence the probability current (density) is in SI units:
j
=
j
e
/
q
=
1
2
m
[
(
Ψ
∗
p
^
Ψ
−
Ψ
p
^
Ψ
∗
)
−
2
q
A
|
Ψ
|
2
]
+
μ
S
q
s
ℏ
∇
×
(
Ψ
∗
S
Ψ
)
{\displaystyle \mathbf {j} =\mathbf {j} _{e}/q={\frac {1}{2m}}\left[\left(\Psi ^{*}\mathbf {\hat {p}} \Psi -\Psi \mathbf {\hat {p}} \Psi ^{*}\right)-2q\mathbf {A} |\Psi |^{2}\right]+{\frac {\mu _{S}}{qs\hbar }}\nabla \times (\Psi ^{*}\mathbf {S} \Psi )}
where S is the spin vector of the particle with corresponding spin magnetic moment μS and spin quantum number s .
It is doubtful if this formula is vaild for particles with an interior structure.[citation needed ] The neutron has zero charge but non-zero magnetic moment, so
μ
S
q
s
ℏ
{\displaystyle {\frac {\mu _{S}}{qs\hbar }}}
would be impossible (except
∇
×
(
Ψ
∗
S
Ψ
)
{\displaystyle \nabla \times (\Psi ^{*}\mathbf {S} \Psi )}
would also be zero in this case). For composite particles with a non-zero charge – like the proton which has spin quantum number s=1/2 and μS = 2.7927·μN or the deuteron (H-2 nucleus) which has s=1 and μS =0.8574·μN [5] – it is mathematically possible but doubtful.
Connection with classical mechanics
edit
The wave function can also be written in the complex exponential (polar ) form:[6]
Ψ
=
R
e
i
S
/
ℏ
{\displaystyle \Psi =Re^{iS/\hbar }}
where R, S are real functions of r and t .
Written this way, the probability density is
ρ
=
Ψ
∗
Ψ
=
R
2
{\displaystyle \rho =\Psi ^{*}\Psi =R^{2}}
and the probability current is:
j
=
ℏ
2
m
i
(
Ψ
∗
∇
Ψ
−
Ψ
∇
Ψ
∗
)
=
ℏ
2
m
i
(
R
e
−
i
S
/
ℏ
∇
R
e
i
S
/
ℏ
−
R
e
i
S
/
ℏ
∇
R
e
−
i
S
/
ℏ
)
=
ℏ
2
m
i
[
R
e
−
i
S
/
ℏ
(
e
i
S
/
ℏ
∇
R
+
i
ℏ
R
e
i
S
/
ℏ
∇
S
)
−
R
e
i
S
/
ℏ
(
e
−
i
S
/
ℏ
∇
R
−
i
ℏ
R
e
−
i
S
/
ℏ
∇
S
)
]
.
{\displaystyle {\begin{aligned}\mathbf {j} &={\frac {\hbar }{2mi}}\left(\Psi ^{*}\mathbf {\nabla } \Psi -\Psi \mathbf {\nabla } \Psi ^{*}\right)\\[5pt]&={\frac {\hbar }{2mi}}\left(Re^{-iS/\hbar }\mathbf {\nabla } Re^{iS/\hbar }-Re^{iS/\hbar }\mathbf {\nabla } Re^{-iS/\hbar }\right)\\[5pt]&={\frac {\hbar }{2mi}}\left[Re^{-iS/\hbar }\left(e^{iS/\hbar }\mathbf {\nabla } R+{\frac {i}{\hbar }}Re^{iS/\hbar }\mathbf {\nabla } S\right)-Re^{iS/\hbar }\left(e^{-iS/\hbar }\mathbf {\nabla } R-{\frac {i}{\hbar }}Re^{-iS/\hbar }\mathbf {\nabla } S\right)\right].\end{aligned}}}
The exponentials and R ∇R terms cancel:
j
=
ℏ
2
m
i
[
i
ℏ
R
2
∇
S
+
i
ℏ
R
2
∇
S
]
.
{\displaystyle \mathbf {j} ={\frac {\hbar }{2mi}}\left[{\frac {i}{\hbar }}R^{2}\mathbf {\nabla } S+{\frac {i}{\hbar }}R^{2}\mathbf {\nabla } S\right].}
Finally, combining and cancelling the constants, and replacing R 2 with ρ ,
j
=
ρ
∇
S
m
.
{\displaystyle \mathbf {j} =\rho {\frac {\mathbf {\nabla } S}{m}}.}
Hence, the spatial variation of the phase of a wavefunction is said to characterize the probability flux of the wavefunction. If we take the familiar formula for the mass flux in hydrodynamics:
j
=
ρ
v
,
{\displaystyle \mathbf {j} =\rho \mathbf {v} ,}
where
ρ
{\displaystyle \rho }
is the mass density of the fluid and v is its velocity (also the group velocity of the wave). In the classical limit, we can associate the velocity with
∇
S
m
,
{\displaystyle {\tfrac {\nabla S}{m}},}
which is the same as equating ∇S with the classical momentum p = m v however, it does not represent a physical velocity or momentum at a point since simultaneous measurement of position and velocity violates uncertainty principle . This interpretation fits with Hamilton–Jacobi theory , in which
p
=
∇
S
{\displaystyle \mathbf {p} =\nabla S}
in Cartesian coordinates is given by ∇S , where S is Hamilton's principal function .
The de Broglie-Bohm theory equates the velocity with
∇
S
m
{\displaystyle {\tfrac {\nabla S}{m}}}
in general (not only in the classical limit) so it is always well defined. It is an interpretation of quantum mechanics.
Motivation
edit
Continuity equation for quantum mechanics
edit
The definition of probability current and Schrödinger's equation can be used to derive the continuity equation , which has exactly the same forms as those for hydrodynamics and electromagnetism .[7]
For some wave function Ψ , let:
ρ
(
r
,
t
)
=
|
Ψ
|
2
=
Ψ
∗
(
r
,
t
)
Ψ
(
r
,
t
)
.
{\displaystyle \rho (\mathbf {r} ,t)=|\Psi |^{2}=\Psi ^{*}(\mathbf {r} ,t)\Psi (\mathbf {r} ,t).}
be the probability density (probability per unit volume, * denotes complex conjugate ). Then,
d
d
t
∫
V
d
V
ρ
=
∫
V
d
V
(
ψ
′
∗
ψ
+
ψ
∗
ψ
′
)
=
∫
V
d
V
(
−
i
ℏ
(
−
ℏ
2
2
m
∇
2
ψ
+
V
ψ
)
ψ
∗
+
i
ℏ
(
−
ℏ
2
2
m
∇
2
ψ
∗
+
V
ψ
∗
)
ψ
)
=
∫
V
d
V
i
ℏ
2
m
(
∇
2
ψ
ψ
∗
−
ψ
∇
2
ψ
∗
)
=
∫
V
d
V
∇
⋅
(
i
ℏ
2
m
(
ψ
∗
∇
ψ
−
ψ
∇
ψ
∗
)
)
=
∫
S
d
a
⋅
(
i
ℏ
2
m
(
ψ
∗
∇
ψ
−
ψ
∇
ψ
∗
)
)
{\displaystyle {\begin{aligned}{\frac {d}{dt}}\int _{\mathcal {V}}dV\,\rho &=\int _{\mathcal {V}}dV\,(\psi '^{*}\psi +\psi ^{*}\psi ')\\&=\int _{\mathcal {V}}dV\,\left(-{\frac {i}{\hbar }}\left(-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}\psi +V\psi \right)\psi ^{*}+{\frac {i}{\hbar }}\left(-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}\psi ^{*}+V\psi ^{*}\right)\psi \right)\\&=\int _{\mathcal {V}}dV\,{\frac {i\hbar }{2m}}(\nabla ^{2}\psi \psi ^{*}-\psi \nabla ^{2}\psi ^{*})\\&=\int _{\mathcal {V}}dV\,\nabla \cdot \left({\frac {i\hbar }{2m}}(\psi ^{*}\nabla \psi -\psi \nabla \psi ^{*})\right)\\&=\int _{\mathcal {S}}d\mathbf {a} \cdot \left({\frac {i\hbar }{2m}}(\psi ^{*}\nabla \psi -\psi \nabla \psi ^{*})\right)\\\end{aligned}}}
where V is any volume and S is the boundary of V .
This is the conservation law for probability in quantum mechanics. The integral form is stated as:
∫
V
(
∂
|
Ψ
|
2
∂
t
)
d
V
+
∫
V
(
∇
⋅
j
)
d
V
=
0
{\displaystyle \int _{V}\left({\frac {\partial |\Psi |^{2}}{\partial t}}\right)\mathrm {d} V+\int _{V}\left(\mathbf {\nabla } \cdot \mathbf {j} \right)\mathrm {d} V=0}
where
j
=
1
2
m
(
Ψ
∗
p
^
Ψ
−
Ψ
p
^
Ψ
∗
)
=
−
i
ℏ
2
m
(
ψ
∗
∇
ψ
−
ψ
∇
ψ
∗
)
=
ℏ
m
Im
(
ψ
∗
∇
ψ
)
{\displaystyle \mathbf {j} ={\frac {1}{2m}}\left(\Psi ^{*}{\hat {\mathbf {p} }}\Psi -\Psi {\hat {\mathbf {p} }}\Psi ^{*}\right)=-{\frac {i\hbar }{2m}}(\psi ^{*}\nabla \psi -\psi \nabla \psi ^{*})={\frac {\hbar }{m}}\operatorname {Im} (\psi ^{*}\nabla \psi )}
is the probability current or probability flux (flow per unit area).
Here, equating the terms inside the integral gives the continuity equation for probability:
∂
∂
t
ρ
(
r
,
t
)
+
∇
⋅
j
=
0
,
{\displaystyle {\frac {\partial }{\partial t}}\rho \left(\mathbf {r} ,t\right)+\nabla \cdot \mathbf {j} =0,}
and the integral equation can also be restated using the divergence theorem as:
∂
∂
t
∫
V
|
Ψ
|
2
d
V
+
{\displaystyle {\frac {\partial }{\partial t}}\int _{V}|\Psi |^{2}\mathrm {d} V+}
S
{\displaystyle \scriptstyle S}
j
⋅
d
S
=
0
{\displaystyle \mathbf {j} \cdot \mathrm {d} \mathbf {S} =0}
.
In particular, if Ψ is a wavefunction describing a single particle, the integral in the first term of the preceding equation, sans time derivative, is the probability of obtaining a value within V when the position of the particle is measured. The second term is then the rate at which probability is flowing out of the volume V . Altogether the equation states that the time derivative of the probability of the particle being measured in V is equal to the rate at which probability flows into V .
By taking the limit of volume integral to include all regions of space, a well-behaved wavefunction that goes to zero at infinities in the surface integral term implies that the time derivative of total probability is zero ie. the normalization condition is conserved.[8] This result is in agreement with the unitary nature of time evolution operators which preserve length of the vector by definition.
Transmission and reflection through potentials
edit
In regions where a step potential or potential barrier occurs, the probability current is related to the transmission and reflection coefficients, respectively T and R ; they measure the extent the particles reflect from the potential barrier or are transmitted through it. Both satisfy:
T
+
R
=
1
,
{\displaystyle T+R=1\,,}
where T and R can be defined by:
T
=
|
j
t
r
a
n
s
|
|
j
i
n
c
|
,
R
=
|
j
r
e
f
|
|
j
i
n
c
|
,
{\displaystyle T={\frac {|\mathbf {j} _{\mathrm {trans} }|}{|\mathbf {j} _{\mathrm {inc} }|}}\,,\quad R={\frac {|\mathbf {j} _{\mathrm {ref} }|}{|\mathbf {j} _{\mathrm {inc} }|}}\,,}
where j inc , j ref , j trans are the incident, reflected and transmitted probability currents respectively, and the vertical bars indicate the magnitudes of the current vectors. The relation between T and R can be obtained from probability conservation:
j
t
r
a
n
s
+
j
r
e
f
=
j
i
n
c
.
{\displaystyle \mathbf {j} _{\mathrm {trans} }+\mathbf {j} _{\mathrm {ref} }=\mathbf {j} _{\mathrm {inc} }\,.}
In terms of a unit vector n normal to the barrier, these are equivalently:
T
=
|
j
t
r
a
n
s
⋅
n
j
i
n
c
⋅
n
|
,
R
=
|
j
r
e
f
⋅
n
j
i
n
c
⋅
n
|
,
{\displaystyle T=\left|{\frac {\mathbf {j} _{\mathrm {trans} }\cdot \mathbf {n} }{\mathbf {j} _{\mathrm {inc} }\cdot \mathbf {n} }}\right|\,,\qquad R=\left|{\frac {\mathbf {j} _{\mathrm {ref} }\cdot \mathbf {n} }{\mathbf {j} _{\mathrm {inc} }\cdot \mathbf {n} }}\right|\,,}
where the absolute values are required to prevent T and R being negative.
Examples
edit
Plane wave
edit
For a plane wave propagating in space:
Ψ
(
r
,
t
)
=
A
e
i
(
k
⋅
r
−
ω
t
)
{\displaystyle \Psi (\mathbf {r} ,t)=\,Ae^{i(\mathbf {k} \cdot {\mathbf {r} }-\omega t)}}
the probability density is constant everywhere;
ρ
(
r
,
t
)
=
|
A
|
2
→
∂
|
Ψ
|
2
∂
t
=
0
{\displaystyle \rho (\mathbf {r} ,t)=|A|^{2}\rightarrow {\frac {\partial |\Psi |^{2}}{\partial t}}=0}
(that is, plane waves are stationary states ) but the probability current is nonzero – the square of the absolute amplitude of the wave times the particle's speed;
j
(
r
,
t
)
=
|
A
|
2
ℏ
k
m
=
ρ
p
m
=
ρ
v
{\displaystyle \mathbf {j} \left(\mathbf {r} ,t\right)=\left|A\right|^{2}{\hbar \mathbf {k} \over m}=\rho {\frac {\mathbf {p} }{m}}=\rho \mathbf {v} }
illustrating that the particle may be in motion even if its spatial probability density has no explicit time dependence.
Particle in a box
edit
For a particle in a box , in one spatial dimension and of length L , confined to the region
0
<
x
<
L
{\displaystyle 0<x<L}
, the energy eigenstates are
Ψ
n
=
2
L
sin
(
n
π
L
x
)
{\displaystyle \Psi _{n}={\sqrt {\frac {2}{L}}}\sin \left({\frac {n\pi }{L}}x\right)}
and zero elsewhere. The associated probability currents are
j
n
=
i
ℏ
2
m
(
Ψ
n
∗
∂
Ψ
n
∂
x
−
Ψ
n
∂
Ψ
n
∗
∂
x
)
=
0
{\displaystyle j_{n}={\frac {i\hbar }{2m}}\left(\Psi _{n}^{*}{\frac {\partial \Psi _{n}}{\partial x}}-\Psi _{n}{\frac {\partial \Psi _{n}^{*}}{\partial x}}\right)=0}
since
Ψ
n
=
Ψ
n
∗
{\displaystyle \Psi _{n}=\Psi _{n}^{*}}
Discrete definition
edit
For a particle in one dimension on
ℓ
2
(
Z
)
,
{\displaystyle \ell ^{2}(\mathbb {Z} ),}
we have the Hamiltonian
H
=
−
Δ
+
V
{\displaystyle H=-\Delta +V}
where
−
Δ
≡
2
I
−
S
−
S
∗
{\displaystyle -\Delta \equiv 2I-S-S^{\ast }}
is the discrete Laplacian, with S being the right shift operator on
ℓ
2
(
Z
)
.
{\displaystyle \ell ^{2}(\mathbb {Z} ).}
Then the probability current is defined as
j
≡
2
ℑ
{
Ψ
¯
i
v
Ψ
}
,
{\displaystyle j\equiv 2\Im \left\{{\bar {\Psi }}iv\Psi \right\},}
with v the velocity operator, equal to
v
≡
−
i
[
X
,
H
]
{\displaystyle v\equiv -i[X,\,H]}
and X is the position operator on
ℓ
2
(
Z
)
.
{\displaystyle \ell ^{2}\left(\mathbb {Z} \right).}
Since V is usually a multiplication operator on
ℓ
2
(
Z
)
,
{\displaystyle \ell ^{2}(\mathbb {Z} ),}
we get to safely write
−
i
[
X
,
H
]
=
−
i
[
X
,
−
Δ
]
=
−
i
[
X
,
−
S
−
S
∗
]
=
i
S
−
i
S
∗
.
{\displaystyle -i[X,\,H]=-i[X,\,-\Delta ]=-i\left[X,\,-S-S^{\ast }\right]=iS-iS^{\ast }.}
As a result, we find:
j
(
x
)
≡
2
ℑ
{
Ψ
¯
(
x
)
i
v
Ψ
(
x
)
}
=
2
ℑ
{
Ψ
¯
(
x
)
(
(
−
S
Ψ
)
(
x
)
+
(
S
∗
Ψ
)
(
x
)
)
}
=
2
ℑ
{
Ψ
¯
(
x
)
(
−
Ψ
(
x
−
1
)
+
Ψ
(
x
+
1
)
)
}
{\displaystyle {\begin{aligned}j\left(x\right)\equiv 2\Im \left\{{\bar {\Psi }}(x)iv\Psi (x)\right\}&=2\Im \left\{{\bar {\Psi }}(x)\left((-S\Psi )(x)+\left(S^{\ast }\Psi \right)(x)\right)\right\}\\&=2\Im \left\{{\bar {\Psi }}(x)\left(-\Psi (x-1)+\Psi (x+1)\right)\right\}\end{aligned}}}
References
edit
^ Paul, Wolfgang; Baschnagel, Jörg (1999). Stochastic Processes : From Physics to Finance . Berlin: Springer. p. 84. ISBN 3-540-66560-9 .
^ Quantum Field Theory, D. McMahon, Mc Graw Hill (USA), 2008, ISBN 978-0-07-154382-8
^ Quantum mechanics, Ballentine, Leslie E, Vol. 280, Englewood Cliffs: Prentice Hall, 1990.
^ see page 473, equation 115.4, L.D. Landau, E.M. Lifschitz. "COURSE OF THEORETICAL PHYSICS Vol. 3 – Quantum Mechanics" (PDF) . ia803206.us.archive.org (3rd ed.). Retrieved 29 April 2023 .
^ "Spin Properties of Nuclei" . www2.chemistry.msu.edu . Retrieved 29 April 2023 .
^ Analytical Mechanics , L.N. Hand, J.D. Finch, Cambridge University Press, 2008, ISBN 978-0-521-57572-0
^ Quantum Mechanics, E. Abers, Pearson Ed., Addison Wesley, Prentice Hall Inc, 2004, ISBN 978-0-13-146100-0
^ Sakurai, Jun John; Napolitano, Jim (2021). Modern quantum mechanics (3rd ed.). Cambridge: Cambridge University Press. ISBN 978-1-108-47322-4 .
Further reading
edit
Resnick, R.; Eisberg, R. (1985). Quantum Physics of Atoms, Molecules, Solids, Nuclei and Particles (2nd ed.). John Wiley & Sons. ISBN 0-471-87373-X .