Past exam of the quantum condensed matter course of the University of Cambridge 2014 1 a Solution Created 2026-10-18 Updated 2026-10-20
Use a representative conduction-electron density , roughly one mobile electron per atom in a dense metal, and the bare electron mass for this estimate. The plasma frequency is thenThis is an order-of-magnitude answer because no particular metal or carrier density is specified. It corresponds to an ordinary frequency , excitation energy , and free-space wavelength near .
The physical oscillation is collective motion of the electron gas against the positive ionic background. Displacing the electrons relative to the ions produces a restoring electric field. In the simplest uniform-displacement picture, : the frequency increases with the density of mobile electric charge and decreases with carrier inertia. Quantizing this longitudinal collective oscillation gives a plasmon.
The same scale controls the ideal collisionless Drude model optical response. Its cold collisionless plasma dielectric response is . Below , the negative relative permittivity makes a transverse electromagnetic wave evanescent, producing strong reflection; above it, the ideal model admits propagation. Real-metal collisions, interband photon absorption and background polarization modify this simple optical edge. Thus is a collective charge-response scale, not a single-electron scattering rate or a sharp transparency threshold for every real metal.
Past exam of the quantum condensed matter course of the University of Cambridge 2014 4 Solution Created 2026-10-18 Updated 2026-10-20
Nondegenerate carrier product. Assume three-dimensional parabolic energy bands near their extrema, spin degeneracy two, and no additional uncounted valley multiplicity. Near the conduction edge, a Boltzmann approximation to the Fermi-Dirac distribution givesAngular integration and the Gaussian radial integral yieldFor valence electron holes, the empty-state probability is . Measuring electron hole kinetic energy downwards from similarly givesThe factors two count spin once in each band. Multiplying cancels the chemical potential and proves the semiconductor mass-action law:In an intrinsic neutral semiconductor, . In a doped equilibrium semiconductor the same product holds under the stated nondegenerate assumptions, even when . Extra valleys or several valence branches must be included in the effective density-of-states factors rather than counted inconsistently.
Equilibrium junction. Put the p region on the left and the n region on the right. Before contact their majority carriers are electron holes and electrons, respectively. Diffusion initially sends electrons from n to p and electron holes from p to n. This leaves negative ionized acceptors on the p side and positive ionized donors on the n side, creating a depletion region and an electric field directed from n to p. At equilibrium this field balances diffusion, so the net electric current vanishes.
For a homojunction the gap remains constant while the conduction band and valence band edges bend downwards from p to n by . The equilibrium chemical potential is flat. Since , the electrostatic potential rises from p to n. The carrier profiles satisfyWith fully ionized nondegenerate dopants, the neutral p side has , , while the neutral n side has , . Thus the electron concentration rises across the junction and the electron hole concentration falls, with both small relative to the far-side majority densities in much of the depleted region. The ideal depletion approximation neglects these mobile populations when calculating the fixed space electric charge, but the actual equilibrium profiles are smooth.
Equilibrium band bending, a flat chemical potential and complementary electron/hole profiles across a homojunction
. Details. The built-in potential follows from these Boltzmann factors:There is no applied terminal bias; the internal electrostatic potential variation does not contradict a constant electrochemical potential.
Opposing carrier fluxes under reverse bias. Use for the magnitude of reverse bias and for reverse electric current, as in the PDF. The junction barrier increases by . The Shockley diode equation in that orientation can be writtenThe reverse-directed contribution is supplied by minority electrons on the p side and minority electron holes on the n side. They reach the depletion edge by diffusion and are swept across by the built-in field. The opposing contribution is associated with majority carriers thermally crossing the barrier in the forward direction; reverse bias suppresses it by the Boltzmann factor. At zero bias detailed balance makes the two equal, while at large enough reverse bias the opposing flux becomes negligible.
In the ideal neutral-region diffusion model,Here are minority-carrier diffusion coefficients, their diffusion lengths, and is junction area. The two terms inside this positive prefactor are the electron and electron hole contributions; both flow in the same reverse-current direction. This separates the two carrier-species contributions from the two opposing fluxes in the Shockley diode equation. Calling the saturation term a drift current refers to its passage across the depleted region, and does not eliminate the minority-carrier diffusion needed to supply it.
Reading the measured curves. At , the ideal equation gives , initially proportional to reverse electric potential difference. At , it approaches a plateau . The thermal electric potential differences at the three temperatures are about –, consistent with a rise at small reverse biases and near saturation at larger ones. Higher temperature produces a much larger plateau because thermal electron-hole excitation raises exponentially.
The lower-temperature curve has a noticeable residual electric potential difference dependence rather than a perfectly flat plateau. Depletion-region generation, surface leakage or other departures from the ideal diffusion model can contribute. The sharp upturn at the largest reverse electric potential differences is reverse breakdown, outside the unextended Shockley diode equation; the plot alone does not identify a unique microscopic breakdown mechanism. The original horizontal axis duplicates the final decade label, so a precise breakdown electric potential difference cannot be read reliably. Neither issue prevents estimating the gap from the intermediate-bias saturation region.
Band-gap estimate. Use plateau electric currents at a common intermediate reverse bias, away from both the small-bias rise and the breakdown upturn. Approximate readings are , –, and . Kelvin temperatures are required in the activation law.
For fixed effective masses, . If the remaining transport prefactor is approximately constant, extracting a semiconductor band gap from diode saturation current therefore usesUsing the widest temperature separation and the rounded readings givesNeglecting all algebraic prefactors and fitting instead gives about . These differ by roughly , a material correction on the graph's precision scale. Thus an honest graph-based conclusion isThe nonflat low-temperature plateau, rough electric current readings, unprovided temperature dependence of , and possible gap variation prevent a more precise extraction. If depletion-region generation dominated rather than the specified mechanism, its activation slope could involve a different energy scale; the displayed estimate explicitly uses the diffusion-dominated saturation model requested here.
For an electron wavepacket in one isolated band under a uniform electric field,Because , acceleration through successive Brillouin zones makes the velocity periodic rather than increasing indefinitely. For ,for a positive field along the chosen axis. Integration givesThe excursion amplitude of a Bloch-oscillating wavepacket is about its mean and peak-to-peak extent . The turning of velocity corresponds semiclassically to Bragg reflection near a zone boundary.
Coherent Bloch oscillations require scattering/dephasing times long compared with and a field weak enough to avoid substantial interband tunnelling. Ordinary metals usually scatter too rapidly; semiconductor superlattices and cold atoms in periodic optical potentials provide cleaner realizations. A completely filled band has canceling velocities, so its net current does not oscillate simply because individual wavepackets can.[ref]
Constant force sweeps crystal momentum linearly, but a periodic band produces oscillating velocity and bounded displacement
. Source. Details. Past exam of the quantum condensed matter course of the University of Cambridge 2016 4 Solution Created 2026-10-18 Updated 2026-10-20
Relaxation-time model. The Drude model treats mobile charge carriers as independent particles accelerated by the Lorentz force between scattering events. Collisions randomize the drift momentum toward local equilibrium and are represented by a uniform, field-independent momentum relaxation time. For the present calculation the carrier density, isotropic mass and relaxation time are fixed, the fields are uniform and stationary, and quantum orbital quantization and intercarrier correlations are neglected. The steady drift velocity obeysHere is the number of carriers of that species per unit volume, their signed charge, their inertial or isotropic band mass, the momentum relaxation time, and their single-species electrical conductivity at zero field.
One-carrier tensor. With and , solving for the electric field givesSince and ,Thus the resistivity tensor, defined by , isIts off-diagonal signs follow from the stated axes and signed . For zero transverse electric current, is the Hall electrical resistivity; reversing charge reverses it. The single-species longitudinal electrical resistivity is independent of in this ideal isotropic model.
Adding electrons and electron holes. Let both species have density and charges , with . The stipulated equal zero-field electrical conductivity means equal positive charge carrier mobility , even if their masses and momentum relaxation times separately differ. Put and . Inverting each electrical resistivity matrix givesBoth species experience the same electric field, so their electric current densities add. Therefore it is their conductivity tensors, rather than their resistivity tensors, that must be summed. The Hall parts cancel andThis is compensated two-carrier magnetoresistance. Although the net Hall field vanishes, each carrier species is bent in its own direction; maintaining the total longitudinal electric current requires an increasing longitudinal field. Equal carrier numbers alone do not ensure zero Hall response if charge carrier mobilities differ; both equalities given here are used.
Large-field limit and quantum oscillations. Expanding the exact diagonal entry,For ,Here the last equality applies to either species with its own mass and transport time; the stipulated charge carrier mobilities give the same threshold. This is the classical strong-field condition : carriers turn substantially before their drift momentum relaxes.
Observable quantum oscillations additionally require sufficiently sharp Landau levels and a sufficiently cold Fermi-Dirac distribution, approximatelyWhen , the disorder field scale is comparable to the onset scale for large classical magnetoresistance. The thermal condition can demand a much larger field, and when even the disorder criterion is stronger. The quantitative damping factors in Q2(c) show that these are crossover scales rather than sharp visibility boundaries. A large smooth magnetoresistance therefore does not itself demonstrate Landau level; it can exist while thermal smearing suppresses oscillations. Conversely, the asymptotic condition is stronger than merely reaching unity, so there is no universal ordering of the two experimental onset fields without temperatures and lifetime information.
Scaling function. The material's actual zero-field longitudinal electrical resistivity is . Dividing the exact answer by this value givesThe assumption that stays constant is essential to eliminating further temperature-dependent parameters. is a dimensionful scaling variable; is its dimensionless version. The corresponding fractional increase is .
What the supplied data test. In the original PDF's left panel, the field dependence is strongly temperature dependent. In the right panel, plotting the fractional magnetoresistance against the squared rescaled field makes the displayed 125, 150, 175, 200 and 225 K curves collapse approximately onto one straight line, over the field intervals where they overlap. This supports the proposed one-parameter scaling and, on that range, the quadratic field law. The 100 K curve does not collapse: it bends away and has a different slope in the rescaled plot. The indicated scaling is therefore supported for the displayed 125–225 K data, not for the entire displayed temperature range. The graph alone does not locate a sharp transition between 100 and 125 K or establish behavior above 225 K.
This empirical collapse is Kohler's rule. It is consistent with a common temperature-dependent charge carrier mobility or relaxation scale and a fixed carrier population; it is not a unique proof of exact compensation or of one microscopic scattering mechanism. A change in relative scattering times, carrier population or Fermi-surface structure could spoil the collapse, and this figure alone does not select the cause of the 100 K departure.
There is a factor-of-two normalization convention in the printed text. The introduced for each species gives , whereas the final paragraph calls “.” If denotes the experimentally measured zero-field electrical resistivity, the consistent slope isIf one instead keeps the literal single-species definition , the same model has . The choice rescales and offsets the axes by temperature-independent constants, so it does not change which temperature curves collapse, but it must be resolved before inferring a carrier density from the slope.
Past exam of the quantum condensed matter course of the University of Cambridge 2018 Section B B2 Solution Created 2026-10-18 Updated 2026-10-20
Equilibrium junction. Put the p-type side on the left and n-type side on the right. When they are joined, electrons diffuse from n to p and electron holes diffuse from p to n. This leaves negative ionized acceptors on the p-side and positive ionized donors on the n-side of the depletion region. The resulting built-in electric field points from n to p, opposing further diffusion. At equilibrium there is one constant Fermi level, and drift and diffusion currents cancel for each carrier species.
In the nondegenerate, fully ionized limit, the neutral-region carrier densities areThe conduction band and valence band lie higher on the p-side and bend down towards the n-side. The gap between them stays fixed for a homojunction. If increases towards the n-side, both edges shift by , whileThus rises and falls continuously across the junction; neither mobile population is an exact step function. Their depletion exposes fixed dopant charge density, rather than an absence of every kind of charge. The built-in voltage is
Light emission. Forward bias lowers the barrier and injects minority carriers across the p-n junction. In a light-emitting diode, electron-hole recombination converts some of their energy into photons. A direct band gap permits efficient radiative transitions with nearly conserved crystal momentum; nonradiative channels compete. Near-edge emission has and . Material composition and confinement select the emission colour. Reverse-bias saturation is not the ordinary LED operating mode.
Reverse current. Take and to mean reverse voltage and reverse current magnitudes. The Shockley diode equation in this convention isThe bias-independent term is the supply of thermally excited minority carriers from the neutral regions, swept across the junction by its electric field. The opposing injection current is reduced by the extra barrier energy and is . Their equality gives zero net current at zero bias; for , the reverse current approaches until leakage or breakdown invalidates this ideal law. The collecting field need not create the minority carriers; their thermal supply limits the diffusion-dominated current.
For an ideal long junction of area , steady diffusion with recombination gives the diode minority-carrier diffusion currentIt depends on doping, geometry, carrier diffusion coefficients and lifetimes, and especially the intrinsic carrier concentration. Depletion-region generation current and surface leakage can add nonideal contributions to the measured reverse current.
Band-gap estimate from the data. Read a nearly flat region well above the small-bias knee and below the high-bias upturn, rather than estimating in the breakdown region. For a diffusion-dominated junction with nearly temperature-independent gap and effective masses,so , where includes transport and doping factors. Extracting a semiconductor band gap from diode saturation current amounts to measuring the slope of against : it is . If the prefactor is taken approximately constant over a short interval,Approximate readings from the original PDF plateaus areUsing the two warmer curves and givesKeeping the factor but taking the residual transport factor constant changes this estimate to approximately . The justified graphical estimate is thereforeThe colder curve is a genuine limitation of the supplied data: the to pair gives only about with a constant prefactor, and the endpoint pair gives about . Thus all three cannot support one precise diffusion-only activation slope. A depletion-region generation current, approximately proportional to rather than , can have a roughly half-gap activation energy; other leakage and temperature-dependent transport can also alter the slope. The PDF alone does not identify their relative contributions. A reliable precision gap measurement would use more temperatures and verify the transport regime instead of forcing one fit through these three readings.
Past exam of the quantum condensed matter course of the University of Cambridge 2021 Section B 6 b Solution Created 2026-10-18 Updated 2026-10-20
At equilibrium, majority electrons undergo diffusion from n to p, while majority electron holes diffuse from p to n. Their departure exposes positive donors on the n side and negative acceptors on the p side. The resulting electric field points from n to p and opposes further diffusion. Equilibrium is reached when the chemical potential is constant across the junction. The built-in potential is this electrostatic adjustment, not a discontinuity of the material's band gap.
Since electronic band energies shift by , both the conduction band and valence band undergo band bending down from p to n by at zero applied bias. Their separation remains the same gap. In the nondegenerate thermal equilibrium limit the built-in voltage can also be expressed using neutral-edge densities as .
Take the n-contact potential and its electrochemical potential as reference. Applying to the p side changes the total depletion-layer drop toThus is reverse bias: it increases the barrier and widens the depletion region. The bands have horizontal tangents at the neutral edges in this approximation and a smooth cubic bend within the linearly graded layer.
A biased junction is not in global thermal equilibrium. Its contact electron electrochemical potentials obeyUnder reverse bias the p-contact level is higher by . A single horizontal equilibrium chemical potential cannot describe this state; electron and hole quasi-Fermi levels describe the nonequilibrium populations. The sketch marks the contact levels. Detailed spatial quasi-Fermi profiles require transport and recombination information beyond the electrostatic doping profile.
At equilibrium the diffusion currents and drift currents cancel separately for each carrier species. Electron diffusion is from n to p, while the electric field drives electrons from p to n; hole diffusion is from p to n, while hole drift is from n to p. The net current is zero, although these opposing microscopic fluxes are present. In conventional-current signs, both majority diffusion contributions point p to n and the balancing drift contributions point n to p.
Forward bias, , reduces the barrier and enhances majority-carrier diffusion and minority carrier injection. The ideal low-injection Shockley diode equation giveswhere positive current is p to n. For reverse bias, majority diffusion is suppressed and the field sweeps thermally supplied minority carriers across the junction, giving approximately . This saturation is limited by carrier supply rather than a current proportional to the field. Depletion-layer generation, recombination and eventual diode reverse breakdown modify the ideal result.
Let the neutral edges be and and take the exterior electric field to vanish. Integrating between the edges gives , so . This is also the condition of equal and opposite integrated depletion charges.
Past exam of the quantum condensed matter course of the University of Cambridge 2022 Section A 1 Solution Created 2026-10-18 Updated 2026-10-20
The drift velocity is the small ensemble-average velocity acquired by charge carriers in an applied electric field. It is distinct from their much larger random thermal or Fermi velocity. For conduction electrons of charge , the isotropic Drude model giveswhere is the momentum relaxation time.
The positive electron mobility measures the ease of this drift in linear response. In the same approximation,Its units are ; in an anisotropic crystal the mobility is a tensor.
The band effective mass describes the acceleration implied by the curvature of the dispersion relation. The group velocity is , while . ThusFor an isotropic parabolic minimum this reduces to and . It incorporates the periodic crystal potential and need not equal the bare electron mass. Negative electronic curvature near a band maximum is conveniently described by positive-mass electron holes.
Use the shock frame, with upstream and downstream fluid velocities and pointing in the same direction. These are the Strong-shock Rankine-Hugoniot conditions for a nonrelativistic monatomic gas with adiabatic index . The two scattering media approach one another at relative speed , to first order in . An energetic electron repeatedly crosses the shock wave, scatters from magnetic turbulence carried by the local plasma, and returns. Each scattering approximately conserves its energy in the local fluid rest frame while changing its direction; comparing successive fluid frames produces a systematic energy and momentum gain. The energy comes from the converging bulk flow, through the moving magnetic field irregularities and their associated electric fields. This is diffusive shock acceleration, a first-order Fermi acceleration process.
Let be the positive cosine of the crossing angle with the shock wave's normal. The Lorentz transformation between fluid frames gives, for either crossing direction with the corresponding sign of the normal reversed,For an almost isotropic population, the number crossing a unit area is weighted by its normal velocity. Thus the normalized crossing distribution is , rather than a uniform distribution in , andFor an ultrarelativistic electron, and its speed , soMore generally gives , while a nonrelativistic kinetic energy has . A complete upstream–downstream–upstream cycle contains two crossings and has ultrarelativistic fractional energy gain .
The initially cold electrons need an injection stage: shock wave heating and small-scale plasma processes must first produce particles with speed much larger than that can return across the front. The result is a description of that accelerated population; a cold electron cannot simply be assigned at its first encounter.
Thermal bremsstrahlung is radiation from the acceleration of electrons in the Coulomb force of ions. At fixed temperature, the encounter rate per electron is proportional to ion number density , and the number of radiating electrons per unit volume is . Therefore in a dilute, optically thin plasma. This is a binary-encounter rate, not a one-particle emission process.
For a fixed encounter geometry, the ionic electric field and hence the electron acceleration are proportional to ionic charge . The Larmor formula is quadratic in acceleration, giving at fixed species number densities. The dimensionless Gaunt factor corrects the simple classical result for the detailed quantum scattering and radiation calculation, after averaging over the thermal electron distribution. It depends on frequency and temperature and, in a more complete notation, ionic charge.
Past exam of the relativity course of the University of Cambridge 2012 Paper 1 2 c Solution 2026-10-20
With Cartesian Minkowski spacetime coordinates , introduce the four-current and electromagnetic four-potential . The antisymmetric electromagnetic field tensor has components and . Maxwell equations becomeThe time and space components reproduce Gauss's law and the Ampère-Maxwell equation; the dual equation gives absence of magnetic monopoles and Faraday's law of induction. Antisymmetry gives charge conservation . The relativistic Lorentz force is . Every displayed equation transforms tensorially under a Lorentz transformation; the electric field and magnetic field are different components of the same object.
Two independent algebraic electromagnetic field invariants areThe second is a pseudoscalar: it is invariant under proper Lorentz transformations but changes sign under spatial reflection. For perpendicular fields, the sign of the first distinguishes a frame with purely magnetic field from one with purely electric field. Both invariants vanish for a plane electromagnetic wave with and , although the field itself is nonzero.
Past exam of the relativity course of the University of Cambridge 2014 Paper 1 2 b Solution 2026-10-20
Use coordinates , signature , and the convention , . The antisymmetric electromagnetic field tensor combines six field components:It transforms as . Lowering indices changes the time–space signs: , while . The covariant tensor instead uses the inverse coordinate Jacobian on each lower index.
For a frame moving with velocity relative to the original one,Electric and magnetic fields therefore mix under a boost; they are frame-dependent parts of one tensor. For example, a purely transverse electric field acquires a magnetic field in a boosted frame. The electromagnetic field invariants and restrict which fields can be eliminated by a frame change. The covariant force law guarantees consistent particle dynamics in all inertial frames.
Past exam of the relativity course of the University of Cambridge 2015 Paper 1 3 Solution 2026-10-20
Use Cartesian coordinates with and . Spatial dot and cross products below use the Euclidean three-metric. The alternating tensor is given the orientation .
Vector and covector transformation. The vector is independent of coordinates, and the chain rule gives . Likewise the one-form is independent of coordinates. Comparing coefficients provesTheir contraction is unchanged because the forward and inverse Jacobians cancel. Here lowering an index uses the metric, not an identification of the numerical upper and lower components.
Why antisymmetrization matters. Write . Differentiating the covector law givesThe last term obstructs the general tensor transformation law. It is symmetric in , so it cancels from . ThereforeEquivalently is the exterior derivative of the electromagnetic four-potential one-form. The statement that is not a tensor needs a qualification: under affine Cartesian Lorentz transformations it does transform tensorially, because the second derivatives vanish. The failure concerns arbitrary nonlinear coordinate changes, including curvilinear coordinates in flat spacetime. This distinction is precisely the difference between general covariance and covariance under the restricted inertial-frame transformations.[ref]
Lecture Notes on General Relativity by Sean M. Carroll. Section 2, Manifolds, equations 2.26–2.27: the extra second-derivative term and its cancellation under antisymmetrization
Electric and magnetic components. Choose , so , and . ThenThus the lower-index field tensor isThis is a derived component representation of the field, with the source's sign convention.
Raised indices and the first invariant. Raising both indices changes the time-space signs and leaves the spatial block unchanged:In the scalar contraction, each independent antisymmetric pair is counted twice, givingUnder a Lorentz transformation each inverse transformation from a lower index cancels its matching transformation from an upper index, leaving the same scalar. In particular a plane electromagnetic wave has in these units and makes this invariant vanish, while a pure electric field makes it negative.
The printed dual convention. The PDF explicitly defines , with no . Each independent pair therefore contributes twice: for example . The complete result isThe time-space pairs contribute and the spatial pairs contribute another . The conventional Hodge dual includes and would instead give the contraction ; it must not silently replace this paper's definition. With fixed orientation the latter contraction is a pseudoscalar: invariant under proper Lorentz transformations but odd under orientation reversal if compared using a fixed alternating symbol.
Energy density and flux. Substitution into the given mixed stress-energy tensor yieldsFor the first formula, and . For the second, . The minus sign on arises from its lowered spatial index; it does not reverse the physical energy flow. Raising that index gives , and symmetry gives , with the Poynting vector in units. Thus is energy density and represents energy flux. In SI coordinates , restore and , with . This keeps the mixed-index signs and the units distinct.
Past exam of the relativity course of the University of Cambridge 2018 Section A 1 b Solution 2026-10-18
The required invariant is obtained by tensor contraction, rather than by an ordinary sum of squares of the entries of the electromagnetic field tensor. With , raising both indices givesThe electromagnetic field tensor is an antisymmetric second-rank tensor, so every independent off-diagonal entry occurs twice in its contraction. The three electric contributions each acquire a minus sign, whereas the magnetic contributions retain their sign. ConsequentlyAll indices are contracted, so this quantity is a Lorentz scalar. Since the speed of light is invariant,is Lorentz invariant. Both terms have the units of squared electric field, in SI units.
Past exam of the relativity course of the University of Cambridge 2019 Section B B4 Solution 2026-10-18
Use signature and local inertial coordinates . Write and , where is rest mass density. The perfect fluid in general relativity has componentsHere is energy density, is energy flux, is momentum density, and is the flux of the th momentum component through a surface normal to . In the local rest frame these reduce to energy density , zero energy and momentum flow, and isotropic pressure . The equation expresses local stress-energy conservation. Its temporal and spatial components give energy and momentum balance; in flat coordinates, for example, .
For relativistic dust, . Contracting its divergence with , and using metric compatibility together with , givesConsequently . Substituting this continuity equation back into the divergence leaves . Wherever matter is present,so the individual particles follow geodesics.
The electromagnetic four-force law must preserve the particle's rest mass. Thus , and requiring this for arbitrary timelike four-velocities makes the electromagnetic field tensor antisymmetric. Its six independent components comprise the electric field and magnetic field. With the stated signature and , one consistent convention isIt gives the spatial force and energy transfer . Under a Lorentz transformation, transforms as a rank-two contravariant tensor; its electric and magnetic components consequently mix. A four-potential supplies , unchanged by . The Maxwell equations areThese combine the source equations and the homogeneous equations into covariant tensor statements; is the four-current.
For charged relativistic dust, the matter tensor exchanges four-momentum with the electromagnetic field: , while their total stress-energy tensor is conserved. Identical particles of rest mass and electric charge have . Contraction with again gives mass continuity because . Removing that term givesUsing zero divergence for the matter tensor alone in this last step would omit the electromagnetic Lorentz force.





