
Nuclear, particle & quantum-information sequence · Interactive lesson
Electric charge, U(1) phase and superselection
Charge conservation is the one ledger line nobody argues about; this is where you earn that confidence. Follow one phase angle from a symmetry of the field to a continuity equation, to an operator commuting with H, and on to a rule that forbids superposing a proton with a neutron — then see what an electron-decay search pays to test it.
UPV-19.04 · Lesson 4 of 9 in this unit. Position is not a completion score.Nuclear, particle & quantum-information sequence. These units continue into nuclear and particle physics and quantum information. Follow your university’s prerequisites, mathematical depth and assessment; this independent pathway is not a universal programme specification.
Build the model
One angle for every charged field, and a ledger no observable can cheat.
The key idea
Electric charge is conserved because the Hamiltonian is unchanged when every charged field is rotated by one common phase; the Hermitian generator of that rotation is the charge operator Q, and [H, Q] = 0 keeps its value fixed.
Example. In neutron decay, n → p + e⁻ + ν̄ₑ, the charge line reads 0 → (+1) + (−1) + 0 in units of e. No measured process has ever failed it, and the search for e⁻ → ν + γ puts the electron's lifetime above 6.6 × 10²⁸ years.
The complete physical model and assumptions
Rotate every charged field by one common phase, ψq → exp(iqα/e) ψq, and nothing measurable changes: every observable is a product in which the charge created equals the charge destroyed, so the angle cancels. That invisibility is a continuous symmetry, and Noether's theorem turns it into a current. For the Schrödinger field ρ = q|ψ|² and j = (qℏ/m) Im(ψ*∇ψ), the equation of motion makes ∂ρ/∂t + ∇⋅j = 0 hold identically, and Q = ∫ρ d³x is constant. As operators it is one line: the rotation is U(α) = exp(iαQ/e) with Hermitian generator Q, invariance of H is U H U† = H, and that is [H, Q] = 0 — so d⟨Q⟩/dt = (i/ℏ)⟨[H, Q]⟩ vanishes and charge eigenstates stay so. The global symmetry buys exactly that: it does not fix e or explain why every charge is a multiple of it, and a tiny charge-changing term is still a violation. Promoting α to a function of position costs a photon, Dμ = ∂μ + (iq/ℏ)Aμ with Aμ shifting to compensate; gauge invariance then forbids a photon mass and makes ∂μ jμ = 0 an identity of Maxwell's equations. Gauss's law then puts Q on a surface at infinity, so Q commutes with every local observable and no experiment sees the phase between components of different charge: a|Q₁⟩ + b|Q₂⟩ is, to every probe, a mixture. That is superselection. The test is an absence — no 255.5 keV photon from e⁻ → νγ, and an electron lifetime above 6.6 × 10²⁸ years.
1A phase you cannot measure is a symmetry
Take every field that carries charge q and multiply it by the same phase, ψq → exp(iqα/e) ψq, with one angle α for the whole system. A neutral field is untouched; an electron field turns by exp(−iα), a proton field by exp(+iα), a Δ⁺⁺ field by exp(+2iα). Every observable is built from products such as ψq* ψq, or the neutron-decay term (ψₚ* ψₙ)(ψₑ* ψν), in which the charge created equals the charge destroyed, so the angle cancels and nothing measurable changes: the rotation is a symmetry. Because α is continuous the rotations form a one-parameter unitary group, U(α) = exp(iαQ/e), and Stone's theorem hands you a Hermitian generator — the charge operator Q, whose eigenvalue on a state holding n elementary charges is ne. Two features matter. The angle lives on a circle, α ≡ α + 2π, so the group is the compact U(1), which is why the eigenvalues of Q/e can all be integers. And this is not the overall phase of the state vector, which is unobservable for any system and generates nothing; it is a relative rotation between components of different charge.
2Noether turns the phase into a current
Noether's theorem promises a conserved current for every continuous symmetry, and for the Schrödinger field you can check it by hand. With iℏ ∂ψ/∂t = −(ℏ²/2m)∇²ψ + Vψ and V real, differentiate |ψ|²: ∂(ψ*ψ)/∂t = (iℏ/2m)(ψ*∇²ψ − ψ∇²ψ*), the potential terms cancelling exactly because V is real — the condition that makes H Hermitian. A complex, absorptive V is just as blind to the phase rotation yet drains |ψ|², so the symmetry alone is not enough: Noether's theorem also needs a real action. The bracket is a divergence, ∇⋅(ψ*∇ψ − ψ∇ψ*), so with ρ = q|ψ|² and j = (qℏ/m) Im(ψ*∇ψ) = (qℏ/2mi)(ψ*∇ψ − ψ∇ψ*) you have ∂ρ/∂t + ∇⋅j = 0 identically, whatever the state. Test it on a plane wave ψ = A exp(i(kx − ωt)): ρ = q|A|² and j = ρ(ℏk/m) x̂, charge density times velocity as a current should be, and both are constant, so the equation reads 0 + 0. Multiply ψ by exp(iα) and neither ρ nor j notices; that blindness is the symmetry, and the current is its bookkeeping. Relativistically the same content is the one equation ∂μ jμ = 0 with jμ = (cρ, j), and for a Dirac field jμ = qc ψ̄ γμ ψ.
3Invariance of H is the commutator [H, Q] = 0
In operator language the argument is one line. The symmetry is U(α) = exp(iαQ/e); H is invariant when U H U† = H for every α; expand to first order in α and that is [H, Q] = 0. The Heisenberg equation then gives d⟨Q⟩/dt = (i/ℏ)⟨[H, Q]⟩ = 0 in every state, and because H and Q commute they share an eigenbasis, so a state of definite charge stays one under time evolution. The S-matrix is built from H, so it commutes with Q too: every reaction balances the charge column. Now look at what a violation would be. In the two-state space of |ν⟩ (Q = 0) and |e⁻⟩ (Q = −e), Q = (e/2)(σz − 1), and a mixing term H′ = εσₓ gives [H′, Q] = −iεe σy, so d⟨Q⟩/dt = (εe/ℏ)⟨σy⟩ — for ε = 1 μeV that is 1.5 × 10⁹ e per second in a state with ⟨σy⟩ = 1. A conservation law is not approximately true when the breaking term is small; it is false, at a rate the term sets, and a rate is what experiments bound. The commutator does not say why the eigenvalues of Q are multiples of e, or why the proton and positron charges agree to a part in 1021. Those facts need more than a global U(1).
4Gauge the phase and a massless photon appears
Let the angle vary from point to point, ψ → exp(iqχ(x)/ℏ) ψ. The derivative now picks up an unwanted piece, ∂μ ψ → exp(iqχ/ℏ)(∂μ + (iq/ℏ) ∂μ χ)ψ, and the kinetic term is no longer invariant. Repair it with a vector field that shifts to cancel the intruder: with Dμ = ∂μ + (iq/ℏ)Aμ and Aμ → Aμ − ∂μ χ, the combination Dμ ψ rotates exactly as ψ does. The field strength Fμν = ∂μ Aν − ∂ν Aμ is invariant, so a Maxwell term is allowed; a mass term m²Aμ Aμ is not, so the photon is massless — the present limit is mγ < 10⁻¹⁸ eV. The field equation ∂μ Fμν = μ₀ jν is where the payoff sits: F is antisymmetric, so ∂ν ∂μ Fμν = 0 identically, and ∂ν jν = 0 is forced by the equations themselves. Charge conservation is no longer one symmetry among many; it is a consistency condition of electrodynamics, and a non-conserved source could not be coupled to the photon at all. The same structure gives Gauss's law, Q = ε₀ ∮E⋅dA over any surface enclosing the system, which means the total charge can be read from arbitrarily far away.
5Why charge sectors never interfere
Gauss's law makes Q measurable on a sphere at spatial infinity, spacelike to any apparatus, so no local operation can change it: every physical observable A commutes with Q. Sandwich [Q, A] = 0 between charge eigenstates: (Q₁ − Q₂)⟨Q₁|A|Q₂⟩ = 0, so ⟨Q₁|A|Q₂⟩ = 0 whenever Q₁ ≠ Q₂. Now take |ψ⟩ = a|Q₁⟩ + b exp(iφ)|Q₂⟩ with |a|² = 0.6 and |b|² = 0.4, and a probe with ⟨A⟩₁ = 2 and ⟨A⟩₂ = −1. The cross term vanishes, so ⟨A⟩ = 0.6(2) + 0.4(−1) = 0.8 for every φ: the relative phase is invisible, and the state cannot be told from the mixture 0.6|Q₁⟩⟨Q₁| + 0.4|Q₂⟩⟨Q₂|. Compare spin, where Sz is also conserved by a rotationally symmetric H but Sₓ connects the sectors: the same weights give ⟨σₓ⟩ = 2√(0.6 × 0.4) cos φ = 0.98 cos φ, and the phase is plainly there. Conservation does not imply superselection; the extra ingredient is that nothing connects the sectors. The rule concerns the total charge of a closed system. A subsystem may sit in a superposition of charges, as a BCS condensate carries an indefinite number of Cooper pairs, with the environment holding the complement.
6What an electron-decay search actually tests
Because the rule is exact, its test is a search for an absence. The cleanest charge-violating process is e⁻ → ν + γ: it balances every other ledger line and fails only charge. Two-body kinematics fixes the signature — an electron at rest shares mₑ c² = 511.0 keV equally between two massless products, so the photon carries 255.5 keV. Borexino watched the 2.5 × 10³¹ electrons of a 75.5 t fiducial volume for 408 days, saw no excess there, and set τ > 6.6 × 10²⁸ years at 90% confidence, some 5 × 10¹⁸ times the age of the universe. The number does not prove conservation; it bounds the coupling of any charge-changing term, and background at 255 keV, not the electron count, sets it — a background-free search of the same exposure would reach 10³¹ years. Theory expects zero, not small: Okun and Zeldovich argued in 1978 that a charge-changing vertex with a massless photon drags in an unbounded cloud of soft photons, so a small violation is not even consistent. It does not test quantisation: the equality of proton and electron charges to a part in 10²¹ is a separate question the global U(1) leaves open.
Key equations
Read the units, assumptions and limits before substituting.
Global U(1) rotation and its generator
ψq → exp(iqα/e) ψq, U(α) = exp(iαQ/e)
One common angle for every charged field. If H is unchanged for every α, the symmetry is exact and Q, the generator, is what it conserves.
α a dimensionless angle, α ≡ α + 2π; Q Hermitian with eigenvalues ne; e = 1.602 × 10⁻¹⁹ CNoether current of the phase symmetry
ρ = q|ψ|², j = (qℏ/m) Im(ψ* ∇ψ)
The conserved current is the symmetry's own bookkeeping: with V real, the Schrödinger equation makes this pair obey ∂ρ/∂t + ∇⋅j = 0 identically.
Schrödinger field; ρ in C m⁻³, j in A m⁻², ψ in m⁻³⁄², m the particle mass in kgContinuity, then the conserved charge
∂μ jμ = 0 ⇔ ∂ρ/∂t + ∇⋅j = 0 ⇒ dQ/dt = 0, Q = ∫ρ d³x
Local bookkeeping becomes a global law: a region can lose charge only through its surface, so all of space loses none.
jμ = (cρ, j); Q in coulombs; needs j → 0 faster than 1/r² on the surface at infinityConservation as a commutator
d⟨Q⟩/dt = (i/ℏ)⟨[H, Q]⟩ = 0 ⇔ [H, Q] = 0
H and Q share an eigenbasis, so a charge eigenstate stays one forever; any charge-changing term in H shows up as a non-zero commutator.
H in J, Q in C, so [H, Q]/ℏ is in C s⁻¹; holds in every state when the commutator vanishes as an operatorGauging the phase
ψ → exp(iqχ(x)/ℏ) ψ, Aμ → Aμ − ∂μ χ, Dμ = ∂μ + (iq/ℏ)Aμ
Makes the local phase unobservable at the price of a photon, and ∂μ Fμν = μ₀ jν then gives ∂ν jν = 0 as an identity.
χ(x) a gauge function in V s, so qχ/ℏ is the local angle; Aμ = (φ/c, −A) in V s m⁻¹; a term m²Aμ Aμ is not invariant, so mγ = 0Superselection of charge
[A, Q] = 0 for every observable A ⇒ ⟨Q₁|A|Q₂⟩ = 0 for Q₁ ≠ Q₂
Charge sectors behave as classical alternatives: a state spanning two charges is, to every probe, the mixture of them.
so ⟨A⟩ = |a|²⟨A⟩₁ + |b|²⟨A⟩₂ in a|Q₁⟩ + b|Q₂⟩; the relative phase of a and b is unobservableInteractive lab
Predict first. Explore the model. Compare the result.
Explore the model
Try this. Drag g to zero, which is what charge superselection asserts about every observable in nature: the fringe collapses onto the dashed mixture line and moving φ changes nothing the probe can see. Then set g = 1 for a spin-like observable and watch p = 0.5 give the deepest fringes.
Live results
- MIXTURE VALUE 2p − 1
- 0.00
- FRINGE AMPLITUDE 2g√(p(1−p))
- 0.60
- ⟨A⟩ AT THIS φ
- 0.30
- ‖[A, Q]‖
- 0.60 e
Results follow this lesson’s model and assumptions. Displayed values are rounded.
Predict → change one variable → check
Plan a change before running it. Every other input stays at the saved baseline.
Saved baseline. Weight p of the Q = 0 sector: 0.50; Charge-changing element g of the probe: 0.60; Relative phase φ: 60 °
This experiment is kept only while this page is open. It does not change saved course progress.
Catch the common trap
The tempting shortcut is that the overall phase of a wavefunction is unobservable, so charge conservation is just that fact dressed up. It is not. Multiplying the entire state vector by exp(iα) is invisible for any system; it is a property of rays in Hilbert space and constrains nothing. The U(1) that conserves charge acts differently on different fields: an electron field picks up exp(−iα), a proton field exp(+iα), a neutrino field nothing at all, and only the relative phases between components of different charge matter. A Hamiltonian is invariant under that rotation only if every term carries equal charge in and out; a term turning an electron into a neutrino is Hermitian, well defined, and not invariant, so [H, Q] ≠ 0 and charge leaks. If the symmetry were the trivial global phase, [H, Q] = 0 would hold for every H and electron-decay searches would have nothing to test.
Practice & worked examples
3 authored examples. Make an attempt, then reveal the reasoning one step at a time.
EasyExample 1
An electron beam is described by ψ = A exp(i(kx − ωt)) with |A|² = 1.0 × 10¹² m⁻³ and k = 8.5 × 10⁹ m⁻¹. Find the charge density and current density carried by the Noether current, check the continuity equation, and show that a global phase leaves both unchanged.
Sketch your answer or reasoning, then compare it with the authored solution. Your writing is not automatically graded or saved.
0 of 6 steps revealed.
- Charge density: ρ = q|ψ|² = (−1.602 × 10⁻¹⁹ C)(1.0 × 10¹² m⁻³) = −1.60 × 10⁻⁷ C m⁻³.
- Current: ψ*∇ψ = |A|²(ik) x̂, so Im(ψ*∇ψ) = k|A|² x̂ and j = (qℏ/m) k|A|² x̂ = ρ (ℏk/m) x̂ — charge density times the velocity ℏk/m.
- Velocity: v = ℏk/m = (1.055 × 10⁻³⁴ J s)(8.5 × 10⁹ m⁻¹)/(9.109 × 10⁻³¹ kg) = 9.84 × 10⁵ m s⁻¹.
- So j = ρv = (−1.60 × 10⁻⁷ C m⁻³)(9.84 × 10⁵ m s⁻¹) = −0.158 A m⁻²: 0.158 A m⁻² of conventional current flowing against the beam, because q < 0.
- Continuity: ρ is uniform and constant and j is uniform, so ∂ρ/∂t = 0 and ∇⋅j = 0 — the equation reads 0 + 0 = 0, and no charge accumulates anywhere along the beam.
- Phase test: replace ψ by exp(iα)ψ. Then |ψ|² and ψ*∇ψ are unchanged, so ρ and j are blind to α; that invisibility is the U(1) symmetry, and j is its Noether current.
Answer. ρ = −1.60 × 10⁻⁷ C m⁻³, j = −0.158 A m⁻² along x (0.158 A m⁻² of conventional current against the beam), ∂ρ/∂t + ∇⋅j = 0 + 0, and both are unchanged by exp(iα).
MediumExample 2
In the two-state space spanned by |ν⟩ (Q = 0) and |e⁻⟩ (Q = −e), write Q as a matrix. A diagonal H₀ is then joined by a hypothetical charge-mixing term H′ = εσₓ with ε = 1.0 μeV. Compute [H, Q], then the initial rate d⟨Q⟩/dt for the state (|ν⟩ + i|e⁻⟩)/√2.
Sketch your answer or reasoning, then compare it with the authored solution. Your writing is not automatically graded or saved.
0 of 5 steps revealed.
- Basis |ν⟩ = (1, 0)ᵀ, |e⁻⟩ = (0, 1)ᵀ, so Q = diag(0, −e) = (e/2)(σz − 1): check, (e/2)diag(1 − 1, −1 − 1) = diag(0, −e).
- [H₀, Q] = 0 because both are diagonal — the charge-conserving part passes. Since [σₓ, σz] = −2iσy, [H′, Q] = ε(e/2)[σₓ, σz] = −iεe σy ≠ 0: the mixing term is the violation and the commutator is its size.
- Heisenberg: d⟨Q⟩/dt = (i/ℏ)⟨[H, Q]⟩ = (i/ℏ)(−iεe)⟨σy⟩ = (εe/ℏ)⟨σy⟩.
- For (1, i)ᵀ/√2: σy (1, i)ᵀ = (−i⋅i, i⋅1)ᵀ = (1, i)ᵀ, so ⟨σy⟩ = +1, while ⟨Q⟩ = −e/2 at t = 0.
- ε/ℏ = (1.0 × 10⁻⁶ eV)/(6.582 × 10⁻¹⁶ eV s) = 1.52 × 10⁹ s⁻¹, so d⟨Q⟩/dt = +1.5 × 10⁹ e per second: at that initial rate the missing half-unit of charge would be restored in about 0.3 ns.
Answer. [H, Q] = −iεe σy and d⟨Q⟩/dt = (εe/ℏ)⟨σy⟩ = +1.5 × 10⁹ e s⁻¹ for that state. A micro-electronvolt of mixing is not a small violation: the lifetime limit on e⁻ → νγ is a bound on ε.
HardExample 3
A liquid-scintillator detector holds 100 t of pseudocumene (C₉H₁₂) in its fiducial volume and runs 408 days searching for e⁻ → ν + γ. Find the photon energy that tags the decay, the number of electrons watched, the lifetime limit if the analysis can exclude at most 560 signal events above background at 90% confidence, and the limit a background-free search of the same exposure would have set.
Sketch your answer or reasoning, then compare it with the authored solution. Your writing is not automatically graded or saved.
0 of 5 steps revealed.
- Two-body decay of an electron at rest into two massless products: each carries mₑ c²/2 = 511.0/2 = 255.5 keV. The tag is a 255.5 keV γ line, sitting among ordinary β and γ backgrounds, which is why the search is background-limited.
- Electrons per molecule of C₉H₁₂: 9 × 6 + 12 × 1 = 66. Molar mass 9(12.011) + 12(1.008) = 120.20 g mol⁻¹, so 100 t = 1.000 × 10⁸ g holds 1.000 × 10⁸/120.20 = 8.32 × 10⁵ mol, or 8.32 × 10⁵ × 6.022 × 10²³ = 5.01 × 10²⁹ molecules.
- Nₑ = 66 × 5.01 × 10²⁹ = 3.31 × 10³¹ electrons. Exposure 408 d = 1.117 yr, so Nₑ t = 3.31 × 10³¹ × 1.117 = 3.70 × 10³¹ electron-years.
- With Γ = 1/τ ≪ 1/t the expected number of decays is Nₑ t Γ. Requiring Nₑ t Γ < 560 gives Γ < 560/(3.70 × 10³¹ yr) = 1.51 × 10⁻²⁹ yr⁻¹, so τ > 6.6 × 10²⁸ yr — about 4.8 × 10¹⁸ times the 1.38 × 10¹⁰ yr age of the universe.
- A background-free search that sees zero events excludes 2.30 expected events at 90%: τ > 3.70 × 10³¹/2.30 = 1.6 × 10³¹ yr. The factor of about 240 between the two limits is the background: what bounds the lifetime is what the detector cannot see through, not how many electrons it holds.
Answer. Eγ = 255.5 keV; Nₑ ≈ 3.3 × 10³¹; τ > 6.6 × 10²⁸ yr at 90% confidence, against about 1.6 × 10³¹ yr if the search had no background.
Quick check
Select an answer and explain it before checking.
Review & continue
Electric charge is conserved because the Hamiltonian is unchanged when every charged field is rotated by one common phase; the Hermitian generator of that rotation is the charge operator Q, and [H, Q] = 0 keeps its value fixed.
- Explain the relationship without looking at the equation.
- Check a sign, a unit and a model assumption in your working.
- Return to the lab or worked example if a step is still unclear.