221A Miscellaneous Notes: Continuity Equation 1 The Continuity Equation
221A Miscellaneous Notes: Continuity Equation 1 The Continuity Equation
Continuity Equation
where the surface element dS ~ is defined outwards and integrated over the
surface ∂V of the volume V . In our case, because we contained the system
completely in the volume, there is nothing at the surface and this term van-
ishes. For the first term, interchanging the order of the volume integral and
the time derivative, it becomes
∂ Z
dV ρ (4)
∂t V
Therefore, the equation is
∂ Z
dV ρ = 0, (5)
∂t V
and hence Z
dV ρ (6)
V
1
is conserved .
In fact, even if
R
the volume does not contain the full system, the amount
of the change in dV ρ is given by the amount of flow out from the volume,
∂ Z Z
~ · ~,
dV ρ = − dS (7)
∂t
which gives the interpretation that the ~(~x, t) is the current of the conserved
quantity, while ρ(~x, t) is the density of the conserved quantity.
2 Examples
This equation shows up in many many different contexts in physics.
2.3 Diffusion
The heat conduction equation is just one example of diffusion equations. For
example, the diffusion of gas is described by
∂n
= D∆n, (10)
∂t
2
where n is the number density, D the diffusion constant. The conserved
quantity is the total number of gas molecules
Z
N= dV n, (11)
~
and the current is ~ = D∇n.
2.4 Electromagnetism
Another example is in the electromagnetism. Look at the Maxwell’s equa-
tions,
~ ·E
~ = 1
∇ ρ, (12)
0
~ ×B
∇ ~ = 1E ~˙ + µ0~, (13)
c2
~ ×E
∇ ~ = −B,~˙ (14)
~ ·B
∇ ~ = 0. (15)
The electric charge density ρ and the electric current density ~ appear. By
taking the time derivative of the first equation, we find
∇ ~˙ = 1 ρ̇,
~ ·E (16)
0
while the divergence of the second equation is
~ · (∇
∇ ~ = 1∇
~ × B) ~˙ + µ0 ∇
~ ·E ~ · ~. (17)
c 2
~ · (∇
Note that ∇ ~ × B)
~ = 0, and hence the second equation becomes
∇ ~˙ = −c2 µ0 ∇
~ ·E ~ · ~. (18)
Comparing Eqs. (16,18), we find that the Maxwell’s equations are consistent
only if
~ · ~.
ρ̇ = −c2 µ0 0 ∇ (19)
Of course, c2 = 1/(µ0 0 ) and this is nothing but the continuity equation. The
conserved quantity is the total electric charge,
Z
Q= dV ρ. (20)
V
3
The continuity equation can be written in a manifestly Lorentz-invariant
fashion. The charge density and the current form a four-vector
~ the continuity
and together with xµ = (ct, ~x) and hence ∂µ = ( 1c ∂t , ∇),
equation becomes
∂µ j µ = 0. (22)
∇µ T µν = 0, (24)
is 100%. If this is the case at a given moment, it should stay that way for
the rest of the evolution; it is supposed to be somewhere in the whole space
anytime. What it means is that the total probability is conserved, and hence
∂P/∂t = 0. Another equivalent statement is that, once the wave function is
4
normalized, it stays normalized forever. This is a consequence of the unitarity
of the time-evolution operator,
hψ, t|ψ, ti = hψ, 0|e+iHt/h̄ e−iHt/h̄ |ψ, 0i = hψ, 0|ψ, 0i. (26)
It sounds funny to say that the probability is conserved, but it has to be.
Correspondingly, there must be a probability current. What we are looking
for is the current ~ that satisfies the continuity equation Eq. (1) with ρ = ψ ∗ ψ.
In the case of a free particle, together with the Schrödinger equation,
∂ h̄2
ih̄ ψ=− ∆ψ, (27)
∂t 2m
we find !
∂ ∗ ih̄ ih̄ ~ · ~.
ψ ψ = ψ∗ ∆ψ − ∆ψ ∗ ψ = −∇ (28)
∂t 2m 2m
By inspection, we find
!
∗h̄ ~ h̄ ~ ∗
~ = ψ ∇ψ − ∇ψ ψ. (29)
2mi 2mi
1 p~ p~
~ = 3
=ρ . (30)
(2πh̄) m m
One important caveat is that this form of the probability current may
not be true depending on the form of the Hamiltonian, e.g., the presence of
the vector potential, if the wave function has more than one particles, etc.
The explicit form needs to be determined knowing what the system is.