Physics Formulary - WEVERS 1ed
Physics Formulary - WEVERS 1ed
Dear reader,
This document contains a 108 page LATEX file which contains a lot equations in physics. It is written
at advanced undergraduate/postgraduate level. It is intended to be a short reference for anyone who
works with physics and often needs to look up equations.
This, and a Dutch version of this file, can be obtained from the author, Johan Wevers
([email protected]).
It can also be obtained on the WWW. See http://www.xs4all.nl/~johanw/index.html, where also
Postscript and PDF versions are available.
If you find any errors or have any comments, please let me know. I am always open for suggestions
and possible corrections to the physics formulary.
The Physics Formulary is made with MiKTEX.
If you prefer the notation in which vectors are typefaced in boldface, uncomment the redefinition of
the \vec command in the TEX file and recompile the file.
This work is licenced under the Creative Commons Attribution 3.0 Netherlands License. To view a
copy of this licence, visit http://creativecommons.org/licenses/by/3.0/nl/ or send a letter to
Creative Commons, 171 Second Street, Suite 300, San Francisco, California 94105, USA.
Johan Wevers
Contents
Contents I
Physical Constants 1
1 Mechanics 2
1.1 Point-kinetics in a fixed coordinate system . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.1.1 Definitions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.1.2 Polar coordinates . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2 Relative motion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.3 Point-dynamics in a fixed coordinate system . . . . . . . . . . . . . . . . . . . . . . . . 2
1.3.1 Force, (angular)momentum and energy . . . . . . . . . . . . . . . . . . . . . . . 2
1.3.2 Conservative force fields . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.3.3 Gravitation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.3.4 Orbital equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.3.5 The virial theorem . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.4 Point dynamics in a moving coordinate system . . . . . . . . . . . . . . . . . . . . . . 4
1.4.1 Apparent forces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.4.2 Tensor notation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.5 Dynamics of masspoint collections . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.5.1 The centre of mass . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.5.2 Collisions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.6 Dynamics of rigid bodies . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.6.1 Moment of Inertia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.6.2 Principal axes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.6.3 Time dependence . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.7 Variational Calculus, Hamilton and Lagrange mechanics . . . . . . . . . . . . . . . . . 7
1.7.1 Variational Calculus . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.7.2 Hamilton mechanics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.7.3 Motion around an equilibrium, linearization . . . . . . . . . . . . . . . . . . . . 7
1.7.4 Phase space, Liouville’s equation . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.7.5 Generating functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
3 Relativity 13
3.1 Special relativity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
3.1.1 The Lorentz transformation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
3.1.2 Red and blue shift . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
3.1.3 The stress-energy tensor and the field tensor . . . . . . . . . . . . . . . . . . . 14
3.2 General relativity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
3.2.1 Riemannian geometry, the Einstein tensor . . . . . . . . . . . . . . . . . . . . . 14
3.2.2 The line element . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
3.2.3 Planetary orbits and the perihelium shift . . . . . . . . . . . . . . . . . . . . . 16
3.2.4 The trajectory of a photon . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
3.2.5 Gravitational waves . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
3.2.6 Cosmology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
4 Oscillations 18
4.1 Harmonic oscillations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
4.2 Mechanic oscillations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
4.3 Electric oscillations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
4.4 Waves in long conductors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
4.5 Coupled conductors and transformers . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
4.6 Pendulums . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
5 Waves 20
5.1 The wave equation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
5.2 Solutions of the wave equation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
5.2.1 Plane waves . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
5.2.2 Spherical waves . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
5.2.3 Cylindrical waves . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
5.2.4 The general solution in one dimension . . . . . . . . . . . . . . . . . . . . . . . 21
5.3 The stationary phase method . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
5.4 Green functions for the initial-value problem . . . . . . . . . . . . . . . . . . . . . . . 22
5.5 Waveguides and resonating cavities . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
5.6 Non-linear wave equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
6 Optics 24
6.1 The bending of light . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
6.2 Paraxial geometrical optics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
6.2.1 Lenses . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
6.2.2 Mirrors . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
6.2.3 Principal planes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
6.2.4 Magnification . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
6.3 Matrix methods . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
6.4 Aberrations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
6.5 Reflection and transmission . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
6.6 Polarization . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
6.7 Prisms and dispersion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
6.8 Diffraction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
6.9 Special optical effects . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
6.10 The Fabry-Perot interferometer . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
7 Statistical physics 30
7.1 Degrees of freedom . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
7.2 The energy distribution function . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
7.3 Pressure on a wall . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
7.4 The equation of state . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
7.5 Collisions between molecules . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
Physics Formulary by ir. J.C.A. Wevers III
8 Thermodynamics 33
8.1 Mathematical introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
8.2 Definitions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
8.3 Thermal heat capacity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
8.4 The laws of thermodynamics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
8.5 State functions and Maxwell relations . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
8.6 Processes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
8.7 Maximal work . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
8.8 Phase transitions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
8.9 Thermodynamic potential . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
8.10 Ideal mixtures . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
8.11 Conditions for equilibrium . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
8.12 Statistical basis for thermodynamics . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
8.13 Application to other systems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
9 Transport phenomena 39
9.1 Mathematical introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
9.2 Conservation laws . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
9.3 Bernoulli’s equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41
9.4 Characterising of flows by dimensionless numbers . . . . . . . . . . . . . . . . . . . . . 41
9.5 Tube flows . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
9.6 Potential theory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
9.7 Boundary layers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
9.7.1 Flow boundary layers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
9.7.2 Temperature boundary layers . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
9.8 Heat conductance . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
9.9 Turbulence . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
9.10 Self organization . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
10 Quantum physics 45
10.1 Introduction to quantum physics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
10.1.1 Black body radiation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
10.1.2 The Compton effect . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
10.1.3 Electron diffraction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
10.2 Wave functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
10.3 Operators in quantum physics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
10.4 The uncertainty principle . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
10.5 The Schrödinger equation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
10.6 Parity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
10.7 The tunnel effect . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
10.8 The harmonic oscillator . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
10.9 Angular momentum . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
10.10 Spin . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
10.11 The Dirac formalism . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
10.12 Atomic physics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
10.12.1 Solutions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
10.12.2 Eigenvalue equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
10.12.3 Spin-orbit interaction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
10.12.4 Selection rules . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
10.13 Interaction with electromagnetic fields . . . . . . . . . . . . . . . . . . . . . . . . . . 50
10.14 Perturbation theory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
10.14.1 Time-independent perturbation theory . . . . . . . . . . . . . . . . . . . . . . 50
10.14.2 Time-dependent perturbation theory . . . . . . . . . . . . . . . . . . . . . . . 51
IV Physics Formulary by ir. J.C.A. Wevers
11 Plasma physics 54
11.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
11.2 Transport . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
11.3 Elastic collisions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
11.3.1 General . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
11.3.2 The Coulomb interaction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
11.3.3 The induced dipole interaction . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
11.3.4 The centre of mass system . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
11.3.5 Scattering of light . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
11.4 Thermodynamic equilibrium and reversibility . . . . . . . . . . . . . . . . . . . . . . . 57
11.5 Inelastic collisions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
11.5.1 Types of collisions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
11.5.2 Cross sections . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
11.6 Radiation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
11.7 The Boltzmann transport equation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
11.8 Collision-radiative models . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
11.9 Waves in plasma’s . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
13 Theory of groups 71
13.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
13.1.1 Definition of a group . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
13.1.2 The Cayley table . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
13.1.3 Conjugated elements, subgroups and classes . . . . . . . . . . . . . . . . . . . . 71
13.1.4 Isomorfism and homomorfism; representations . . . . . . . . . . . . . . . . . . . 72
13.1.5 Reducible and irreducible representations . . . . . . . . . . . . . . . . . . . . . 72
13.2 The fundamental orthogonality theorem . . . . . . . . . . . . . . . . . . . . . . . . . . 72
13.2.1 Schur’s lemma . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
13.2.2 The fundamental orthogonality theorem . . . . . . . . . . . . . . . . . . . . . . 72
13.2.3 Character . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
13.3 The relation with quantum mechanics . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
13.3.1 Representations, energy levels and degeneracy . . . . . . . . . . . . . . . . . . . 73
13.3.2 Breaking of degeneracy by a perturbation . . . . . . . . . . . . . . . . . . . . . 73
13.3.3 The construction of a base function . . . . . . . . . . . . . . . . . . . . . . . . 73
13.3.4 The direct product of representations . . . . . . . . . . . . . . . . . . . . . . . 74
13.3.5 Clebsch-Gordan coefficients . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
13.3.6 Symmetric transformations of operators, irreducible tensor operators . . . . . . 74
13.3.7 The Wigner-Eckart theorem . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
13.4 Continuous groups . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
13.4.1 The 3-dimensional translation group . . . . . . . . . . . . . . . . . . . . . . . . 75
13.4.2 The 3-dimensional rotation group . . . . . . . . . . . . . . . . . . . . . . . . . . 76
13.4.3 Properties of continuous groups . . . . . . . . . . . . . . . . . . . . . . . . . . . 76
13.5 The group SO(3) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
13.6 Applications to quantum mechanics . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78
13.6.1 Vectormodel for the addition of angular momentum . . . . . . . . . . . . . . . 78
13.6.2 Irreducible tensor operators, matrixelements and selection rules . . . . . . . . . 78
13.7 Applications to particle physics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
14 Nuclear physics 81
14.1 Nuclear forces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81
14.2 The shape of the nucleus . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
14.3 Radioactive decay . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82
14.4 Scattering and nuclear reactions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
14.4.1 Kinetic model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
14.4.2 Quantum mechanical model for n-p scattering . . . . . . . . . . . . . . . . . . . 83
14.4.3 Conservation of energy and momentum in nuclear reactions . . . . . . . . . . . 84
14.5 Radiation dosimetry . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 84
15.13.3 Quantumchromodynamics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 94
15.14 Path integrals . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
15.15 Unification and quantum gravity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
16 Astrophysics 96
16.1 Determination of distances . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
16.2 Brightness and magnitudes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96
16.3 Radiation and stellar atmospheres . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
16.4 Composition and evolution of stars . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
16.5 Energy production in stars . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
The ∇-operator 99
Mechanics
When the acceleration is constant this gives: v(t) = v0 + at and s(t) = s0 + v0 t + 21 at2 .
For the unit vectors in a direction ⊥ to the orbit ⃗et and parallel to it ⃗en holds:
⃗v d⃗r v ⃗e˙t
⃗et = = ⃗e˙t = ⃗en ; ⃗en =
|⃗v | ds ρ |⃗e˙t |
d⃗
p d(m⃗v ) d⃗v dm m=const
F⃗ (⃗r, ⃗v , t) = = =m + ⃗v = m⃗a
dt dt dt dt
Chapter 1: Mechanics 3
Z2 Z2
The work A, delivered by a force, is A = F⃗ · d⃗s = F cos(α)ds
1 1
∂U
τ =−
∂θ
P⃗ P
Hence, the conditions for a mechanical equilibrium are: Fi = 0 and ⃗τi = 0.
The force of friction is usually proportional to the force perpendicular to the surface, except when
the motion starts, when a threshold has to be overcome: Ffric = f · Fnorm · ⃗et .
So the work delivered by a conservative force field depends not on the trajectory covered but only on
the starting and ending points of the motion.
1.3.3 Gravitation
The Newtonian law of gravitation is (in GRT one also uses κ instead of G):
m1 m2
F⃗g = −G 2 ⃗er
r
The gravitational potential is then given by V = −Gm/r. From Gauss law it then follows: ∇2 V =
4πGϱ.
ℓ
r(θ) = , or: x2 + y 2 = (ℓ − εx)2
1 + ε cos(θ − θ0 )
with
L2 2W L2 ℓ ℓ k
ℓ= ; ε2 = 1 + 2 =1− ; a= =
Gµ2 Mtot G2 µ3 Mtot a 1 − ε2 2W
a is half the length of the
√ long axis of the elliptical orbit in case the orbit is closed. Half the length
of the short axis is b = aℓ. ε is the excentricity of the orbit. Orbits with an equal ε are of equal
shape. Now, 5 types of orbits are possible:
5. k > 0 and ε > 1: a hyperbole, curved away from the centre of force.
Other combinations are not possible: the total energy in a repulsive force field is always positive so
ε > 1.
If the surface between the orbit covered between t1 and t2 and the focus C around which the planet
moves is A(t1 , t2 ), Kepler’s 2nd law is
LC
A(t1 , t2 ) = (t2 − t1 )
2m
Kepler’s 3rd law is, with T the period and Mtot the total mass of the system:
T2 4π 2
3
=
a GMtot
d2 xα
m = Fα
dt2
will be transformed into:
d2 xα dxβ dxγ
m 2
+ Γα
βγ = Fα
dt dt dt
dxβ dxγ
The apparent forces are taken from he origin to the effect side in the way Γα
βγ .
dt dt
⃗ = m1⃗r1 + m2⃗r2
R
m1 + m2
With ⃗r = ⃗r1 − ⃗r2 , the kinetic energy becomes: T = 12 Mtot Ṙ2 + 12 µṙ2 , with the reduced mass µ given
1 1 1
by: = +
µ m1 m2
The motion within and outside the centre of mass can be separated:
⃗˙ outside = ⃗τoutside ;
L ⃗˙ inside = ⃗τinside
L
1.5.2 Collisions
With collisions, where B are the coordinates of the collision and C an arbitrary other position, holds:
p⃗ = m⃗vm is constant, and T = 21 m⃗vm2 is constant. The changes in the relative velocities can be derived
from: S⃗ = ∆⃗ ⃗ C = CB
p = µ(⃗vaft − ⃗vbefore ). Further holds ∆L ⃗ × S,
⃗ p⃗ ∥ S⃗ =constant and L ⃗ w.r.t. B is
constant.
where I is the moment of inertia with respect to a central axis, which is given by:
X
I= mi⃗ri 2 ; T ′ = Wrot = 12 ωIij ⃗ei⃗ej = 12 Iω 2
i
Object I Object I
Zb
du d
δ L(q, q̇, t)dt = 0 with δ(a) = δ(b) = 0 and δ = (δu)
dx dx
a
The Hamiltonian of a Harmonic oscillator is given by H(x, p) = p2 /2m + 21 mω 2 x2 . With new coordi-
p √
nates (θ, I), obtained by the canonical transformation x = 2I/mω cos(θ) and p = − 2Imω sin(θ),
with inverse θ = arctan(−p/mωx) and I = p2 /2mω + 12 mωx2 it follows: H(θ, I) = ωI.
The Hamiltonian of a charged particle with charge q in an external electromagnetic field is given by:
1 2
H= ⃗ + qV
p⃗ − q A
2m
This Hamiltonian can be derived from the Hamiltonian of a free particle H = p2 /2m with the
⃗ and H → H − qV . This is elegant from a relativistic point of view: this is
transformations p⃗ → p⃗ − q A
equivalent to the transformation of the momentum 4-vector pα → pα − qAα . A gauge transformation
on the potentials Aα corresponds with a canonical transformation, which make the Hamilton equations
the equations of motion for the system.
With T = 21 (Mik q̇i q̇k ) one receives the set of equations M q̈ + V q = 0. If qi (t) = ai exp(iωt) is
substituted, this set of equations has solutions if det(V −ω 2 M ) = 0. This leads to the eigenfrequencies
aT V ak
of the problem: ωk2 = Tk . If the equilibrium is stable holds: ∀k that ωk2 > 0. The general solution
ak M ak
is a superposition if eigenvibrations.
8 Physics Formulary by ir. J.C.A. Wevers
∂ϱ
{ϱ, H} + =0
∂t
For an arbitrary quantity A holds:
dA ∂A
= {A, H} +
dt ∂t
Liouville’s theorem can than be written as:
Z
dϱ
= 0 ; or: pdq = constant
dt
one can derive the following Hamilton equations with the new Hamiltonian K:
dQi ∂K dPi ∂K
= ; =−
dt ∂Pi dt ∂Qi
Now, a distinction between 4 cases can be made:
dF1 (qi , Qi , t)
1. If pi q̇i − H = Pi Qi − K(Pi , Qi , t) − , the coordinates follow from:
dt
∂F1 ∂F1 ∂F1
pi = ; Pi = − ; K=H+
∂qi ∂Qi ∂t
dF2 (qi , Pi , t)
2. If pi q̇i − H = −Ṗi Qi − K(Pi , Qi , t) + , the coordinates follow from:
dt
∂F2 ∂F2 ∂F2
pi = ; Qi = ; K=H+
∂qi ∂Pi ∂t
dF3 (pi , Qi , t)
3. If −ṗi qi − H = Pi Q̇i − K(Pi , Qi , t) + , the coordinates follow from:
dt
∂F3 ∂F3 ∂F3
qi = − ; Pi = − ; K=H+
∂pi ∂Qi ∂t
dF4 (pi , Pi , t)
4. If −ṗi qi − H = −Pi Qi − K(Pi , Qi , t) + , the coordinates follow from:
dt
∂F4 ∂F4 ∂F4
qi = − ; Qi = ; K=H+
∂pi ∂Pi ∂t
⃗ = µ0 (H
⃗ +M
⃗ ) = µ0 µr H,
⃗ ⃗ = P m/Vol, µ0 nm20
B M ⃗ µr = 1 + χm , with χm =
3kT
Q1 Q2 ⃗
F⃗12 = ⃗er ; ⃗ = F
E
4πε0 εr r2 Q
The Lorentzforce is the force which is felt by a charged particle that moves through a magnetic field.
The origin of this force is a relativistic transformation of the Coulomb force: F⃗L = Q(⃗v ×B
⃗ ) = l(I×
⃗ B⃗ ).
The magnetic field in point P which results from an electric current is given by the law of Biot-Savart,
also known as the law of Laplace. In here, d⃗l ∥ I⃗ and ⃗r points from d⃗l to P :
⃗P = µ0 I ⃗
dB dl × ⃗er
4πr2
If the current is time-dependent one has to take retardation into account: the substitution I(t) →
I(t − r/c) has to be applied.
Z2
The potentials are given by: V12 = − ⃗ · d⃗s and A
E ⃗ = 1B⃗ × ⃗r.
2
1
10 Physics Formulary by ir. J.C.A. Wevers
Here, the freedom remains to apply a gauge transformation. The fields can be derived from the
potentials as follows:
⃗
E⃗ = −∇V − ∂ A , B
⃗ =∇×A ⃗
∂t
⃗ = ⃗v × E.
Further holds the relation: c2 B ⃗
The energy density of the electromagnetic wave of a vibrating dipole at a large distance is:
p20 sin2 (θ)ω 4 p20 sin2 (θ)ω 4 ck 4 |⃗
p |2
w = ε0 E 2 = sin2 (kr − ωt) , ⟨w⟩t = , P =
16π 2 ε0 r2 c4 2
32π ε0 r c 2 4 12πε0
Chapter 2: Electricity & Magnetism 11
∂2 ∂2
2 µ ∂ ⃗ 2 µ ∂ ⃗
∇ − εµ 2 − E=0, ∇ − εµ 2 − B=0
∂t ρ ∂t ∂t ρ ∂t
This results in a damped wave: E ⃗ = E exp(−k ′′⃗n · ⃗r ) exp(i(k ′⃗n · ⃗r − ωt)). If the material is a good
r
µω
conductor, the wave vanishes after approximately one wavelength, k = (1 + i) .
2ρ
2.6 Multipoles
∞ ′ l
1 1X r Q X kn
Because = Pl (cos θ) the potential can be written as: V =
|⃗r − ⃗r ′ | r 0 r 4πε n rn
For the lowest-order terms this results in:
Monopole: l = 0, k0 = ρdV
R
Dipole: l = 1, k1 = r cos(θ)ρdV
R
1. The electric dipole: dipole moment: p⃗ = Ql⃗e , where ⃗e goes from ⊕ to ⊖, and F⃗ = (⃗ ⃗ ext ,
p · ∇)E
and W = −⃗ p·E⃗ out .
⃗ Q 3⃗p · ⃗r ⃗ out
Electric field: E ≈ − p⃗ . The torque is: ⃗τ = p⃗ × E
4πεr3 r2
√
2. The magnetic dipole: dipole moment: if r ≫ A: µ ⃗ = I⃗ × (A⃗e⊥ ), F⃗ = (⃗ ⃗ out
µ · ∇)B
2
mv⊥ ⃗ out
|µ| = , W = −⃗ µ×B
2B
⃗ −µ 3µ · ⃗r ⃗ out
Magnetic field: B = −µ ⃗ . The moment is: ⃗τ = µ ⃗ ×B
4πr3 r2
⃗
E⃗ dep = E ⃗0 = − N P
⃗ mat − E
ε0
⃗ dep = H
H ⃗ mat − H
⃗ 0 = −N M⃗
N is a constant depending only on the shape of the object placed in the field, with 0 ≤ N ≤ 1. For
a few limiting cases of an ellipsoid holds: a thin plane: N = 1, a long, thin bar: N = 0, a sphere:
N = 31 .
Relativity
where
1
γ=r
2
1 − v2
c
The velocity difference ⃗v ′ between two observers transforms according to:
−1
⃗v1 · ⃗v2 ⃗v1 · ⃗v2
⃗v ′ = γ 1− ⃗v2 + (γ − 1) ⃗
v 1 − γ⃗
v 1
c2 v12
⃗ = Q (1 − β 2 )⃗er
E
4πε0 r2 (1 − β 2 sin2 (θ))3/2
Length, mass and time transform according to: ∆tr = γ∆t0 , mr = γm0 , lr = l0 /γ, with 0 the
quantities in a co-moving reference frame and r the quantities in a frame moving with velocity v
w.r.t. it. The proper time τ is defined as: dτ 2 = ds2 /c2 , so ∆τ = ∆t/γ. For energy and momentum
14 Physics Formulary by ir. J.C.A. Wevers
The equations of motion for a charged particle in an EM field become with the field tensor:
dpα
= qFαβ uβ
dτ
1. The geodesic postulate: free falling particles move along geodesics of space-time with the proper
time τ or arc length s as parameter. For particles with zero rest mass R (photons), the use of a
free parameter is required because for them holds ds = 0. From δ ds = 0 the equations of
motion can be derived:
d2 xα β
α dx dx
γ
+ Γ βγ =0
ds2 ds ds
Chapter 3: Relativity 15
µ
The Ricci tensor is a contraction of the Riemann tensor: Rαβ := Rαµβ , which is symmetric: Rαβ =
Rβα . The Bianchi identities are: ∇λ Rαβµν + ∇ν Rαβλµ + ∇µ Rαβνλ = 0.
The Einstein tensor is given by: Gαβ := Rαβ − 21 g αβ R, where R := Rα
α
is the Ricci scalar, for which
p
holds: ∇β Gαβ = 0. With the variational principle δ (L(gµν ) − Rc2 /16πκ) |g|d4 x = 0 for variations
R
8πκ 8πκ
Gαβ = Tαβ , which can also be written as Rαβ = (Tαβ − 21 gαβ Tµµ )
c2 c2
For empty space this is equivalent to Rαβ = 0. The equation Rαβµν = 0 has as only solution a flat
space.
The Einstein equations are 10 independent equations, which are of second order in gµν . From this,
the Laplace equation from Newtonian gravitation can be derived by stating: gµν = ηµν + hµν , where
|h| ≪ 1. In the stationary case, this results in ∇2 h00 = 8πκϱ/c2 .
8πκ
The most general form of the field equations is: Rαβ − 21 gαβ R + Λgαβ = Tαβ
c2
where Λ is the cosmological constant. This constant plays a role in inflatory models of the universe.
In general holds: ds = gµν dx dx . In special relativity this becomes ds2 = −c2 dt2 + dx2 + dy 2 + dz 2 .
2 µ ν
r > 2m: r
r r t
u = − 1 exp cosh
2m 4m 4m
r
r r t
v = − 1 exp sinh
2m 4m 4m
r < 2m: r
r r t
u = 1− exp sinh
2m 4m 4m
r
r r t
v = 1− exp cosh
2m 4m 4m
r = 2m: here, the Kruskal coordinates are singular, which is necessary to eliminate the coordi-
nate singularity there.
The line element in these coordinates is given by:
32m3 −r/2m 2
ds2 = − e (dv − du2 ) + r2 dΩ2
r
The line r = 2m corresponds to u = v = 0, the limit x0 → ∞ with u = v and x0 → −∞ with u = −v.
The Kruskal coordinates are only singular on the hyperbole v 2 − u2 = 1, this corresponds with r = 0.
On the line dv = ±du holds dθ = dφ = ds = 0.
For the metric outside a rotating, charged spherical mass the Newman metric applies:
2mr − e2 r2 + a2 cos2 θ
2 2 2
ds = 1− 2 c dt − dr2 − (r2 + a2 cos2 θ)dθ2 −
r + a2 cos2 θ r2 − 2mr + a2 − e2
(2mr − e2 )a2 sin2 θ 2a(2mr − e2 )
2
r2 + a2 + sin θdφ2
+ sin2 θ(dφ)(cdt)
r2 + a2 cos2 θ r2 + a2 cos2 θ
where m = κM/c2 , a = L/M c and e = κQ/ε0 c2 . √
A rotating charged black hole has an event horizon with RS = m + m2 − a2 − e2 .
Near rotating black holes frame dragging occurs√because gtφ ̸= 0. For the Kerr metric (e = 0, a ̸= 0)
then follows that within the surface RE = m + m2 − a2 cos2 θ (de ergosphere) no particle can be at
rest.
du d2 u
+ u − 3mu = 0
dφ dφ2
3.2.6 Cosmology
If for the universe as a whole is assumed:
1. There exists a global time coordinate which acts as x0 of a Gaussian coordinate system,
2. The 3-dimensional spaces are isotrope for a certain value of x0 ,
3. Each point is equivalent to each other point for a fixed x0 .
then the Robertson-Walker metric can be derived for the line element:
R2 (t)
ds2 = −c2 dt2 + 2 2 2
2 (dr + r dΩ )
kr
r02 1 − 2
4r0
Oscillations
where Ψ̂ is the amplitude. A superposition of several harmonic oscillations with the same frequency
results in another harmonic oscillation:
X
Ψ̂i cos(αi ± ωt) = Φ̂ cos(β ± ωt)
i
with:
P
Ψ̂i sin(αi ) X XX
i
tan(β) = P and Φ̂2 = Ψ̂2i + 2 Ψ̂i Ψ̂j cos(αi − αj )
Ψ̂i cos(αi ) i j>i i
i
dn x(t)
Z
x(t)
For harmonic oscillations holds: x(t)dt = and = (iω)n x(t).
iω dtn
F F
x= , and for the velocity holds: ẋ = √
m(ω02 − ω 2 ) + ikω i Cmδ + k
ω ω0
where δ = − . The quantity Z = F/ẋ is called the impedance of the system. The quality of the
ω0 ω √
Cm
system is given by Q = .
k
The frequency with minimal
√ |Z| is called velocity resonance frequency. This is equal to ω0 . In the
resonance curve |Z|/ Cm is plotted
√ against ω/ω0 . The width of this curve is characterized by the
points where |Z(ω)| = |Z(ω0 )| 2. In these points holds: R = X and δ = ±Q−1 , and the width is
2∆ωB = ω0 /Q.
The damping frequency r ωD is a measure for the time in which an oscillating system comes to rest. It
1
is given by ωD = ω0 1 − . A weak damped oscillation (k 2 < 4mC) dies out after TD = 2π/ωD .
4Q2
For a critical damped oscillation (k 2 = 4mC) holds ωD = 0. A strong damped oscillation (k 2 > 4mC)
drops like (if k 2 ≫ 4mC) x(t) ≈ x0 exp(−t/τ ).
Chapter 4: Oscillations 19
4.6 Pendulums
The oscillation time T = 1/f , and for different types of pendulums is given by:
Oscillating spring: T = 2π m/C if the spring force is given by F = C · ∆l.
p
Physical pendulum: T = 2π I/τ with τ the moment of force and I the moment of inertia.
p
2lm
Torsion pendulum: T = 2π
p
I/κ with κ = the constant of torsion and I the moment
πr4 ∆φ
of inertia.
Mathematical pendulum: T = 2π
p
l/g with g the acceleration of gravity and l the length of
the pendulum.
Chapter 5
Waves
The phase velocity is given by vph = ω/k. The group velocity is given by:
dω dvph k dn
vg = = vph + k = vph 1 −
dk dk n dk
where n is the refractive index of the medium. If vph does not depend on ω holds: vph = vg . In a
dispersive medium it is possible that vg > vph or vg < vph , and vg · vf = c2 . If one wants to transfer
information with a wave, e.g. by modulation of an EM wave, the information travels with the velocity
at with a change in the electromagnetic field propagates. This velocity is often almost equal to the
group velocity.
For some media, the propagation velocity follows from:
s
gλ 2πγ 2πh
Surface waves on a liquid: v = + tanh
2π ϱλ λ
√
where h is the depth of the liquid and γ the surface tension. If h ≪ λ holds: v ≈ gh.
The equation for a harmonic traveling plane wave is: u(⃗x, t) = û cos(⃗k · ⃗x ± ωt + φ)
If waves reflect at the end of a spring this will result in a change in phase. A fixed end gives a phase
change of π/2 to the reflected wave, with boundary condition u(l) = 0. A lose end gives no change
in the phase of the reflected wave, with boundary condition (∂u/∂x)l = 0.
If an observer is moving w.r.t. the wave with a velocity vobs , he will observe a change in frequency:
f vf − vobs
the Doppler effect. This is given by: = .
f0 vf
1 ∂ 2 (ru) ∂ 2 (ru)
− =0
v 2 ∂t2 ∂r2
with general solution:
f (r − vt) g(r + vt)
u(r, t) = C1 + C2
r r
1 ∂2u 1 ∂
∂u
− r =0
v 2 ∂t2 r ∂r ∂r
This is a Bessel equation, with solutions which can be written as Hankel functions. For sufficient
large values of r these are approximated by:
û
u(r, t) = √ cos(k(r ± vt))
r
where bm ∈ IR. Substituting u(x, t) = Aei(kx−ωt) gives two solutions ωj = ωj (k) as dispersion
relations. The general solution is given by:
Z∞
u(x, t) = a(k)ei(kx−ω1 (k)t) + b(k)ei(kx−ω2 (k)t) dk
−∞
Because in general the frequencies ωj are non-linear in k there is dispersion and the solution cannot
be written any more as a sum of functions depending only on x ± vt: the wave front transforms.
d
significantly to the integral are areas with a stationary phase, determined by (kx − ω(k)t) = 0.
dk
Now the following approximation is possible:
Z∞ v
N u
i(kx−ω(k)t)
X u 2π
t d2 ω(k ) exp −i 41 π + i(ki x − ω(ki )t)
a(k)e dk ≈
i
−∞ i=1 dki2
p
This gives a discrete solution ψℓ with eigenvalue γℓ2 : k = εµω 2 − γℓ2 . For ω < ωℓ , k is
imaginary and the wave is damped. Therefore, ωℓ is called the cut-off frequency. In rectangular
conductors the following expression can be found for the cut-off frequency for modes TEm,n of
TMm,n :
2
λℓ = p
(m/a)2 + (n/b)2
3. Ez and Bz are zero everywhere: the Transversal electromagnetic mode (TEM). Than holds:
√
k = ±ω εµ and vf = vg , just as if here were no waveguide. Further k ∈ IR, so there exists no
cut-off frequency.
In a rectangular, 3 dimensional resonating cavity with edges a, b and c the possible wave numbers
n1 π n2 π n3 π
are given by: kx = , ky = , kz = This results in the possible frequencies f = vk/2π
a b c
in the cavity: r
v n2x n2y n2z
f= + +
2 a2 b2 c2
For a cubic cavity, with a = b = c, the possible number of oscillating modes NL for longitudinal waves
is given by:
4πa3 f 3
NL =
3v 3
Because transversal waves have two possible polarizations holds for them: NT = 2NL .
This is an energy equation. Energy is conserved if the left-hand side is 0. If x2 > 1/β, the right-hand
side changes sign and an increase in energy changes into a decrease of energy. This mechanism limits
the growth of oscillations.
The Korteweg-De Vries equation is given by:
∂u ∂u ∂u ∂3u
+ − au + b2 3 = 0
∂t ∂x | {z∂x} ∂x }
| {z
non−lin dispersive
This equation is for example a model for ion-acoustic waves in a plasma. For this equation, soliton
solutions of the following form exist:
−d
u(x − ct) = 2
cosh (e(x − ct))
Optics
In these equations the following signs are being used for refraction at a spherical surface, as is seen
by an incoming light ray:
Quantity + −
R Concave surface Convex surface
f Converging lens Diverging lens
v Real object Virtual object
b Virtual image Real image
6.2.2 Mirrors
For images of mirrors holds:
2
h2 1
1 1 1 2 1
= + = + −
f v b R 2 R v
where h is the perpendicular distance from the point the light ray hits the mirror to the optical
axis. Spherical aberration can be reduced by not using spherical mirrors. A parabolical mirror has
no spherical aberration for light rays parallel with the optical axis and is therefore often used for
telescopes. The used signs are:
Quantity + −
R Concave mirror Convex mirror
f Concave mirror Convex mirror
v Real object Virtual object
b Real image Virtual image
6.2.4 Magnification
b
The linear magnification is defined by: N = −
v
αsyst
The angular magnification is defined by: Nα = −
αnone
where αsys is the size of the retinal image with the optical system and αnone the size of the retinal
image without the system. Further holds: N · Nα = 1. For a telescope holds: N = fobjective /focular .
The f-number is defined by f /Dobjective .
1 0
1. Transfer along length l: MR =
l/n 1
1 −D
2. Refraction at a surface with dioptric power D: MT =
0 1
6.4 Aberrations
Lenses usually do not give a perfect image. Some causes are:
1. Chromatic aberration is caused by the fact that n = n(λ). This can be partially corrected
with a lens which is composed of more lenses with different functions ni (λ). Using N lenses
makes it possible to obtain the same f for N wavelengths.
2. Spherical aberration is caused by second-order effects which are usually ignored; a spherical
surface does not make a perfect lens. Incomming rays far from the optical axis will more bent.
3. Coma is caused by the fact that the principal planes of a lens are only flat near the principal
axis. Further away of the optical axis they are curved. This curvature can be both positive or
negative.
4. Astigmatism: from each point of an object not on the optical axis the image is an ellipse
because the thickness of the lens is not the same everywhere.
6. Distorsion gives abberations near the edges of the image. This can be corrected with a com-
bination of positive and negative lenses.
where E0r is the reflected amplitude and E0t the transmitted amplitude. Then the Fresnel equations
are:
tan(θi − θt ) sin(θt − θi )
r∥ = , r⊥ =
tan(θi + θt ) sin(θt + θi )
2 sin(θt ) cos(θi ) 2 sin(θt ) cos(θi )
t∥ = , t⊥ =
sin(θt + θi ) cos(θt − θi ) sin(θt + θi )
The following holds: t⊥ − r⊥ = 1 and t∥ + r∥ = 1. If the coefficient of reflection R and transmission
T are defined as (with θi = θr ):
Ir It cos(θt )
R≡ and T ≡
Ii Ii cos(θi )
6.6 Polarization
Ip Imax − Imin
The polarization is defined as: P = =
Ip + Iu Imax + Imin
where the intensity of the polarized light is given by Ip and the intensity of the unpolarized light
is given by Iu . Imax and Imin are the maximum and minimum intensities when the light passes a
polarizer. If polarized light passes through a polarizer Malus law applies: I(θ) = I(0) cos2 (θ) where
θ is the angle of the polarizer.
The state of a light ray can be described by the Stokes-parameters: start with 4 filters which each
transmits half the intensity. The first is independent of the polarization, the second and third are
linear polarizers with the transmission axes horizontal and at +45◦ , while the fourth is a circular
polarizer which is opaque for L-states. Then holds S1 = 2I1 , S2 = 2I2 − 2I1 , S3 = 2I3 − 2I1 and
S4 = 2I4 − 2I1 .
The state of a polarized light ray can also be described by the Jones vector:
E0x eiφx
⃗
E=
E0y eiφy
6.8 Diffraction
Fraunhofer diffraction occurs far away from the source(s). The Fraunhofer diffraction of light passing
through multiple slits is described by:
2 2
I(θ) sin(u) sin(N v)
= ·
I0 u sin(v)
where u = πb sin(θ)/λ, v = πd sin(θ)/λ. N is the number of slits, b the width of a slit and d the
distance between the slits. The maxima in intensity are given by d sin(θ) = kλ.
The diffraction through a spherical aperture with radius a is described by:
2
I(θ) J1 (ka sin(θ))
=
I0 ka sin(θ)
The diffraction pattern of a rectangular aperture at distance R with length a in the x-direction and
b in the y-direction is described by:
2 2
sin(α′ ) sin(β ′ )
I(x, y)
=
I0 α′ β′
In case n2 = n3 ̸= n1 , which happens e.g. at trigonal, hexagonal and tetragonal crystals there
is one optical axis in the direction of n1 . Incident light rays can now be split up in two parts:
the ordinary wave is linear polarized ⊥ the plane through the transmission direction and the
optical axis. The extraordinary wave is linear polarized in the plane through the transmission
direction and the optical axis. Dichroism is caused by a different absorption of the ordinary
and extraordinary wave in some materials. Double images occur when the incident ray makes
an angle with the optical axis: the extraordinary wave will refract, the ordinary will not.
Retarders: waveplates and compensators. Incident light will have a phase shift of ∆φ =
2πd(|n0 − ne |)/λ0 if an uniaxial crystal is cut in such a way that the optical axis is parallel with
the front and back plane. Here, λ0 is the wavelength in vacuum and n0 and ne the refractive
indices for the ordinary and extraordinary wave. For a quarter-wave plate holds: ∆φ = π/2.
The Kerr-effect: isotropic, transparent materials can become birefringent when placed in an
⃗ The difference in refractive index in
electric field. In that case, the optical axis is parallel to E.
2
the two directions is given by: ∆n = λ0 KE , where K is the Kerr constant of the material. If
the electrodes have an effective length ℓ and are separated by a distance d, the retardation is
given by: ∆φ = 2πKℓV 2 /d2 , where V is the applied voltage.
The Pockels or linear electro-optical effect can occur in 20 (from a total of 32) crystal symmetry
classes, namely those without a centre of symmetry. These crystals are also piezoelectric: their
polarization changes when a pressure is applied and vice versa: P⃗ = pd+ε0 χE. ⃗ The retardation
3
in a Pockels cell is ∆φ = 2πn0 r63 V /λ0 where r63 is the 6-3 element of the electro-optic tensor.
The Faraday effect: the polarization of light passing through material with length d and to
which a magnetic field is applied in the propagation direction is rotated by an angle β = VBd
where V is the Verdet constant.
C̆erenkov radiation arises when a charged particle with vq > vf arrives. The radiation is
emitted within a cone with an apex angle α with sin(α) = c/cmedium = c/nvq .
Statistical physics
Because vibrational degrees of freedom account for both kinetic and potential energy they count
double. So, for linear molecules this results in a total of s = 6n − 5. For non-linear molecules this
gives s = 6n − 6. The average energy of a molecule in thermodynamic equilibrium is ⟨Etot ⟩ = 21 skT .
Each degree of freedom of a molecule has in principle the same energy: the principle of equipartition.
ℏ2
Wrot = l(l + 1) = Bl(l + 1) , Wvib = (v + 12 )ℏω0
2I
The vibrational levels are excited if kT ≈ ℏω, the rotational levels of a hetronuclear molecule are
excited if kT ≈ 2B. For homonuclear molecules additional selection rules apply so the rotational
levels are well coupled if kT ≈ 6B.
1
1. Even s: s = 2l: c(s) =
(l − 1)!
2l
2. Odd s: s = 2l + 1: c(s) = √
π(2l − 1)!!
Chapter 7: Statistical physics 31
ZZZ Z∞ Zπ Z2π
3
d N= nAvτ cos(θ)P (v, θ, φ)dvdθdφ
0 0 0
From this follows for the particle flux on the wall: Φ = 14 n ⟨v⟩. For the pressure on the wall then
follows:
2mv cos(θ)d3 N 2
d3 p = , so p = n ⟨E⟩
Aτ 3
1
N m v2
pV = ns RT =
3
Here, ns is the number of moles particles and N is the total number of particles within volume V . If
the own volume and the intermolecular forces cannot be neglected the Van der Waals equation can
be derived:
an2s
p + 2 (V − bns ) = ns RT
V
There is an isotherme with a horizontal point of inflection. In the Van der Waals equation this
corresponds with the critical temperature, pressure and volume of the gas. This is the upper limit of
the area of coexistence between liquid and vapor. From dp/dV = 0 and d2 p/dV 2 = 0 follows:
8a a
Tcr = , pcr = , Vcr = 3bns
27bR 27b2
For the critical point holds: pcr Vm,cr /RTcr = 38 , which differs from the value of 1 which follows from
the general gas law.
Scaled on the critical quantities, with p∗ := p/pcr , T ∗ = T /Tcr and Vm∗ = Vm /Vm,cr with Vm := V /ns
holds:
3
p∗ + Vm∗ − 13 = 38 T ∗
∗
(Vm ) 2
Gases behave the same for equal values of the reduced quantities: the law of the corresponding states.
A virial expansion is used for even more accurate views:
1 B(T ) C(T )
p(T, Vm ) = RT + + + ···
Vm Vm2 Vm3
The Boyle temperature TB is the temperature for which the 2nd virial coefficient is 0. In a Van der
Waals gas, this happens at TB = a/Rb. The inversion temperature Ti = 2TB .
with a minimum ϵ at r = rm . The following holds: D ≈ 0.89rm . For the Van der Waals coefficients
a and b and the critical quantities holds: a = 5.275NA2 D3 ϵ, b = 1.3NA D3 , kTkr = 1.2ϵ and Vm,kr =
3.9NA D3 .
A more simple model for intermolecular forces assumes a potential U (r) = ∞ for r < D, U (r) = ULJ
for D ≤ r ≤ 3D and U (r) = 0 for r ≥ 3D. This gives for the potential energy of one molecule:
Z 3D
Epot = U (r)F (r)dr.
D
with F (r) the spatial distribution function in spherical coordinates, which for a homogeneous distri-
bution is given by: F (r)dr = 4nπr2 dr.
Some useful mathematical relations are:
Z∞ Z∞ √ Z∞
n −x 2n −x2 (2n)! π 2
x e dx = n! , x e dx = 2n+1
, x2n+1 e−x dx = 12 n!
n!2
0 0 0
Chapter 8
Thermodynamics
8.2 Definitions
1 ∂p
The isochoric pressure coefficient: βV =
p ∂T V
1 ∂V
The isothermal compressibility: κT = −
V ∂p T
1 ∂V
The isobaric volume coefficient: γp =
V ∂T p
1 ∂V
The adiabatic compressibility: κS = −
V ∂p S
For an ideal gas holds: Cmp − CmV = R. Further, if the temperature is high enough to thermalize
all internal rotational and vibrational degrees of freedom, holds: CV = 12 sR. Hence Cp = 12 (s + 2)R.
For their ratio now follows γ = (2 + s)/s. For a lower T one needs only to consider the thermalized
degrees of freedom. For a Van der Waals gas holds: CmV = 12 sR + ap/RT 2 .
In general holds:
2
∂p ∂V ∂V ∂p
Cp − CV = T · = −T ≥0
∂T V ∂T p ∂T p ∂V T
Because (∂p/∂V )T is always < 0, the following is always valid: Cp ≥ CV . If the coefficient of
expansion is 0, Cp = CV , and also at T = 0K.
dQ
dS ≥
T
If the only processes occurring are reversible holds: dS = d Qrev /T . So, the entropy difference after
a reversible process is:
Z2
d Qrev
S2 − S1 =
T
1
I
d Qrev
So, for a reversible cycle holds: = 0.
T
I
d Qirr
For an irreversible cycle holds: < 0.
T
The third law of thermodynamics is (Nernst):
∂S
lim =0
T →0 ∂X T
From this it can be concluded that the thermal heat capacity → 0 if T → 0, so absolute zero
temperature cannot be reached by cooling through a finite number of steps.
From the total differential and the definitions of CV and Cp it can be derived that:
∂p ∂V
T dS = CV dT + T dV and T dS = Cp dT − T dp
∂T V ∂T p
for an enlarged surface holds: d Wrev = −γdA, with γ the surface tension. From this follows:
∂U ∂F
γ= =
∂A S ∂A T
8.6 Processes
Work done
The efficiency η of a process is given by: η =
Heat added
Cold delivered
The Cold factor ξ of a cooling down process is given by: ξ =
Work added
Reversible adiabatic processes
For adiabatic processes holds: W = U1 − U2 . For reversible adiabatic processes holds Poisson’s
equation: with γ = Cp /CV one gets that pV γ =constant. Also holds: T V γ−1 =constant and
T γ p1−γ =constant. Adiabatics exhibit a greater steepness p-V diagram than isothermics because
γ > 1.
Isobaric processes
R2
Here holds: H2 − H1 = 1
Cp dT . For a reversible isobaric process holds: H2 − H1 = Qrev .
The throttle process
This is also called the Joule-Kelvin effect and is an adiabatic expansion of a gas through a porous
material or a small opening. Here H is a conserved quantity, and dS > 0. In general this is
accompanied with a change in temperature. The quantity which is important here is the throttle
coefficient: " #
∂T 1 ∂V
αH = = T −V
∂p H Cp ∂T p
The inversion temperature is the temperature where an adiabatically expanding gas keeps the same
temperature. If T > Ti the gas heats up, if T < Ti the gas cools down. Ti = 2TB , with for TB :
[∂(pV )/∂p]T = 0. The throttle process is e.g. applied in refridgerators.
The Carnotprocess
The system undergoes a reversible cycle with 2 isothemics and 2 adiabatics:
so G has a twist in the transition point. In a two phase system Clapeyron’s equation is valid:
dp S α − Sm
β
rβα
= m β
=
dT α
Vm − Vm (Vm − Vmβ )T
α
For an ideal gas one finds for the vapor line at some distance from the critical point:
p = p0 e−rβα/RT
There exist also phase transitions with rβα = 0. For those there will occur only a discontinuity in the
second derivates of Gm . These second-order transitions appear at organization phenomena.
A phase-change of the 3rd order, so with e.g. [∂ 3 Gm /∂T 3 ]p non continuous arises e.g. when ferro-
magnetic iron changes to the paramagnetic state.
Chapter 8: Thermodynamics 37
where xi = ni /n is the molar fraction of component i. The molar volume of a mixture of two
components can be a concave line in a V -x2 diagram: the mixing contracts the volume.
The thermodynamic
P potentials are not independent in a multiple-phase
P system. It can be derived
that ni dµi = −SdT + V dp, this gives at constant p and T : xi dµi = 0 (Gibbs-Duhmen).
i i
Each component has as much µ’s as there are phases. The number of free parameters in a system
with c components and p different phases is given by f = c + 2 − p.
For the thermodynamic potentials holds: µi = µ0i + RT ln(xi ) < µ0i . A mixture of two liquids is
rarely ideal: this is usually only the case for chemically related components or isotopes. In spite of
this holds Raoult’s law for the vapour pressure holds for many binary mixtures: pi = xi p0i = yi p.
Here is xi the fraction of the ith component in liquid phase and yi the fraction of the ith component
in gas phase.
A solution of one component in another gives rise to an increase in the boiling point ∆Tk and a
decrease of the freezing point ∆Ts . For x2 ≪ 1 holds:
RTk2 RT 2
∆Tk = x2 , ∆Ts = − s x2
rβα rγβ
with rβα the evaporation heat and rγβ < 0 the melting heat. For the osmotic pressure Π of a solution
0
holds: ΠVm1 = x2 RT .
The most probable distribution, that with the maximum value for P , is the equilibrium state. When
Stirling’s equation, ln(n!) ≈ n ln(n)−n is used, one finds for a discrete system the Maxwell-Boltzmann
distribution. The occupation numbers in equilibrium are then given by:
N Wi
ni = gi exp −
Z kT
P
The state sum Z is a normalization constant, given by: Z = gi exp(−Wi /kT ). For an ideal gas
i
holds:
V (2πmkT )3/2
Z=
h3
The entropy can then be defined as: S = k ln(P ) . For a system in thermodynamic equilibrium this
becomes: N
U Z U Z
S= + kN ln + kN ≈ + k ln
T N T N!
V (2πmkT )3/2
For an ideal gas, with U = 32 kT then holds: S = 25 kN + kN ln
N h3
Transport phenomena
dX ∂X
This results in general to: = + (⃗v · ∇)X .
dt ∂t
ZZZ ZZZ ZZ
d ∂
From this follows that also holds: Xd3 V = Xd3 V + ⃝ X(⃗v · ⃗n )d2 A
dt ∂t
where the volume V is surrounded by surface A. Some properties of the ∇ operator are:
Here, ⃗v is an arbitrary vector field and ϕ an arbitrary scalar field. Some important integral theorems
are: ZZ ZZZ
Gauss: ⃝ (⃗v · ⃗n )d2 A = (div⃗v )d3 V
I ZZ
Stokes for a scalar field: (ϕ · ⃗et )ds = (⃗n × gradϕ)d2 A
I ZZ
Stokes for a vector field: (⃗v · ⃗et )ds = (rot⃗v · ⃗n )d2 A
ZZ
This results in: ⃝ (rot⃗v · ⃗n )d2 A = 0
ZZ ZZZ
2
Ostrogradsky: ⃝ (⃗n × ⃗v )d A = (rot⃗v )d3 A
ZZ ZZZ
2
⃝ (ϕ⃗n )d A = (gradϕ)d3 V
T can be split in a part pI representing the normal tensions and a part T′ representing the shear
stresses: T = T′ + pI, where I is the unit tensor. When viscous aspects can be ignored holds:
divT= −gradp.
When the flow velocity is ⃗v at position ⃗r holds on position ⃗r + d⃗r:
The quantity L:=grad⃗v can be split in a symmetric part D and an antisymmetric part W. L = D + W
with
1 ∂vi ∂vj 1 ∂vi ∂vj
Dij := + , Wij := −
2 ∂xj ∂xi 2 ∂xj ∂xi
When the rotation or vorticity ω⃗ = rot⃗v is introduced holds: Wij = 12 εijk ωk . ω
⃗ represents the local
⃗ 1 ⃗
rotation velocity: dr · W = 2 ω × dr.
For a Newtonian liquid holds: T′ = 2ηD. Here, η is the dynamical viscosity. This is related to the
shear stress τ by:
∂vi
τij = η
∂xj
For compressible media can be stated: T′ = (η ′ div⃗v )I + 2ηD. From equating the thermodynamical
and mechanical pressure it follows: 3η ′ + 2η = 0. If the viscosity is constant holds: div(2D) =
∇2⃗v + grad div⃗v .
The conservation laws for mass, momentum and energy for continuous media can be written in both
integral and differential form. They are:
Integral notation:
ZZZ ZZ
∂
1. Conservation of mass: ϱd3 V + ⃝ ϱ(⃗v · ⃗n )d2 A = 0
∂t
ZZZ ZZ ZZZ ZZ
∂ 3 2
2. Conservation of momentum: ϱ⃗v d V + ⃝ ϱ⃗v (⃗v · ⃗n )d A = f0 d V + ⃝ ⃗n · T d2 A
3
∂t
ZZZ ZZ
∂
3. Conservation of energy: ( 2 v + e)ϱd V + ⃝ ( 12 v 2 + e)ϱ(⃗v · ⃗n )d2 A =
1 2 3
∂t
ZZ ZZZ ZZ
− ⃝ (⃗q · ⃗n )d2 A + (⃗v · f⃗0 )d3 V + ⃝ (⃗v · ⃗n T)d2 A
Differential notation:
∂ϱ
1. Conservation of mass: + div · (ϱ⃗v ) = 0
∂t
∂⃗v
2. Conservation of momentum: ϱ + (ϱ⃗v · ∇)⃗v = f⃗0 + divT = f⃗0 − gradp + divT′
∂t
ds de p dϱ
3. Conservation of energy: ϱT =ϱ − = −div⃗q + T′ : D
dt dt ϱ dt
Here, e is the internal energy per unit of mass E/m and s is the entropy per unit of mass S/m.
⃗ is the heat flow. Further holds:
⃗q = −κ∇T
∂E ∂e ∂E ∂e
p=− =− , T = =
∂V ∂1/ϱ ∂S ∂s
so
∂e ∂h
CV = and Cp =
∂T V ∂T p
Chapter 9: Transport phenomena 41
For compressible flows holds: 12 v 2 + gh + p/ϱ =constant along aRline of flow. If also holds rot⃗v = 0
and the entropy is equal on each streamline holds 12 v 2 + gh + dp/ϱ =constant everywhere. For
incompressible flows this becomes: 12 v 2 + gh + p/ϱ =constant everywhere. For ideal gases with
constant Cp and CV holds, with γ = Cp /CV :
1 2 γ p c2
2v + = 12 v 2 + = constant
γ−1ϱ γ−1
With a velocity potential defined by ⃗v = gradϕ holds for instationary flows:
Z
∂ϕ 1 2 dp
+ 2 v + gh + = constant everywhere
∂t ϱ
Ma: (velocity)/(speed of sound): objects moving faster than approximately Ma = 0,8 produce
shockwaves which propagate with an angle θ with the velocity of the object. For this angle
holds Ma= 1/ arctan(θ).
Now, the dimensionless Navier-Stokes equation becomes, with x′ = x/L, ⃗v ′ = ⃗v /V , grad′ = Lgrad,
∇′2 = L2 ∇2 and t′ = tω:
∂⃗v ′ ′ ′ ′ ′ ⃗g ∇′2⃗v ′
Sr + (⃗
v · ∇ )⃗
v = −grad p + +
∂t′ Fr Re
1 dp 2
v(r) = − (R − r2 )
4η dx
ZR
π dp 4
For the volume flow holds: ΦV = v(r)2πrdr = − R
8η dx
0
For non viscous media, if p = p(ϱ) and all forces are conservative, Kelvin’s theorem can be derived:
dΓ
=0
dt
Chapter 9: Transport phenomena 43
For rotationless flows a velocity potential ⃗v = gradϕ can be introduced. In the incompressible case
follows from conservation of mass ∇2 ϕ = 0. For a 2-dimensional flow a flow function ψ(x, y) can be
defined: with ΦAB the amount of liquid flowing through a curve s between the points A and B:
ZB ZB
ΦAB = (⃗v · ⃗n )ds = (vx dy − vy dx)
A A
and the definitions vx = ∂ψ/∂y, vy = −∂ψ/∂x holds: ΦAB = ψ(B) − ψ(A). In general holds:
∂2ψ ∂2ψ
+ = −ωz
∂x2 ∂y 2
1 ∂ψ ∂ϕ ∂ψ 1 ∂ϕ
vr = = , vθ = − =
r ∂θ ∂r ∂r r ∂θ
Q
For source flows with power Q in (x, y) = (0, 0) holds: ϕ = ln(r) so that vr = Q/2πr, vθ = 0.
2π
For a dipole of strength Q in x = a and strength −Q in x = −a follows from superposition: ϕ =
−Qax/2πr2 where Qa is the dipole strength. For a vortex holds: ϕ = Γθ/2π.
If an object is surrounded by an uniform main flow with ⃗v = v⃗ex and such a large Re that viscous
effects are limited to the boundary layer holds: Fx = 0 and Fy = −ϱΓv. The statement that Fx = 0 is
d’Alembert’s paradox and originates from the neglection of viscous effects. The lift Fy is also created
by η because Γ ̸= 0 due to viscous effects. Henxe rotating bodies also create a force perpendicular to
their direction of motion: the Magnus effect.
Z∞ Z∞
2 ∗ ∂vx
ϑv = (v − vx )vx dy , δ v = (v − vx )dy and τ0 = −η
∂y y=0
0 0
∂vx dp
The boundary layer is released from the surface if = 0. This is equivalent with =
∂y y=0 dx
12ηv∞
.
δ2
∂T κ ∂2T
= +Φ
∂t ϱc ∂x2
where Φ is a source term. If Φ = 0 the solutions for harmonic oscillations at x = 0 are:
T − T∞ x x
= exp − cos ωt −
Tmax − T∞ D D
p
with D = 2κ/ωϱc. At x = πD the temperature variation is in anti-phase with the surface. The
one-dimensional solution at Φ = 0 is
x2
1
T (x, t) = √ exp −
2 πat 4at
9.9 Turbulence
√
The time scale of turbulent velocity variations τt is of the order of: τt = τ Re/Ma2 with τ the
molecular time scale. For the velocity of the particles holds: v(t) = ⟨v⟩ + v ′ (t) with ⟨v ′ (t)⟩ = 0. The
Navier-Stokes equation now becomes:
where SRij = −ϱ ⟨vi vj ⟩ is the turbulent stress tensor. Boussinesq’s assumption is: τij = −ϱ vi′ vj′ .
It is stated that, analogous to Newtonian media: SR = 2ϱνt ⟨D⟩. Near a boundary holds: νt = 0, far
away of a boundary holds: νt ≈ νRe.
From this follows that in a two-dimensional flow the energy flux goes towards large values of k: larger
structures become larger at the expanse of smaller ones. In three-dimensional flows the situation is
just the opposite.
Chapter 10
Quantum physics
This is also written as ⟨f (t)⟩ = ⟨Φ|f |Φ⟩. The normalizing condition for wavefunctions follows from
this: ⟨Φ|Φ⟩ = ⟨Ψ|Ψ⟩ = 1.
∆A · ∆B ≥ 12 | ⟨ψ|[A, B]|ψ⟩ |
From this follows: ∆E · ∆t ≥ 12 ℏ, and because [x, px ] = iℏ holds: ∆px · ∆x ≥ 21 ℏ, and ∆Lx · ∆Ly ≥
1
2 ℏLz .
ℏ2 2 ∂ψ
− ∇ ψ + U ψ = Eψ = iℏ
2m ∂t
The linear combination of the solutions of this equation give the general solution. In one dimension
it is: X Z
iEt
ψ(x, t) = + dE c(E)uE (x) exp −
ℏ
ℏ
The current density J is given by: J = (ψ ∗ ∇ψ − ψ∇ψ ∗ )
2im
∂P (x, t)
The following conservation law holds: = −∇J(x, t)
∂t
10.6 Parity
The parity operator in one dimension is given by Pψ(x) = ψ(−x). If the wavefunction is split in even
and odd functions, it can be expanded into eigenfunctions of P:
[P, H] = 0. The functions ψ + = 21 (1+P)ψ(x, t) and ψ − = 12 (1−P)ψ(x, t) both satisfy the Schrödinger
equation. Hence, parity is a conserved quantity.
10.10 Spin
For the spin operators are defined by their commutation relations: [Sx , Sy ] = iℏSz . Because the
spin operators do not act in the physical space (x, y, z) the uniqueness of the wavefunction is not a
criterium here: also half odd-integer values are allowed for the spin. Because [L, S] = 0 spin and
angular momentum operators do not have a common set of eigenfunctions. The spin operators are
⃗⃗ ⃗σ , with
given by S = 21 ℏ⃗
⃗ 0 1 ⃗ 0 −i ⃗ 1 0
⃗σ x = , ⃗σ y = , ⃗σ z =
1 0 i 0 0 −1
The eigenstates of Sz are called spinors: χ = α+ χ+ + α− χ− , where χ+ = (1, 0) represents the state
with spin up (Sz = 21 ℏ) and χ− = (0, 1) represents the state with spin down (Sz = − 12 ℏ). Then the
probability to find spin up after a measurement is given by |α+ |2 and the chance to find spin down
is given by |α− |2 . Of course holds |α+ |2 + |α− |2 = 1.
The electron will have an intrinsic magnetic dipole moment M ⃗ due to its spin, given by M⃗ =
⃗
−egS S/2m, with gS = 2(1 + α/2π + · · ·) the gyromagnetic ratio. In the presence of an external
magnetic field this gives a potential energy U = −M ⃗ · B.⃗ The Schrödinger equation then becomes
(because ∂χ/∂xi ≡ 0):
∂χ(t) egS ℏ ⃗
iℏ = ⃗σ · Bχ(t)
∂t 4m
with ⃗σ = (⃗⃗σ x , ⃗
⃗σ y , ⃗ ⃗ = B⃗ez there are two eigenvalues for this problem: χ± for E =
⃗σ z ). If B
±egS ℏB/4m = ±ℏω. So the general solution is given by χ = (ae−iωt , beiωt ). From this can be
derived: ⟨Sx ⟩ = 21 ℏ cos(2ωt) and ⟨Sy ⟩ = 21 ℏ sin(2ωt). Thus the spin precesses about the z-axis with
frequency 2ω. This causes the normal Zeeman splitting of spectral lines.
The potential operator for two particles with spin ± 21 ℏ is given by:
1 ⃗ ⃗
V (r) = V1 (r) + (S1 · S2 )V2 (r) = V1 (r) + 21 V2 (r)[S(S + 1) − 23 ]
ℏ2
This makes it possible for two states to exist: S = 1 (triplet) or S = 0 (Singlet).
1 ∂2 m20 c2
2
∇ − 2 2− 2 ψ(⃗x, t) = 0
c ∂t ℏ
Clm
Ylm = √ Plm (cos θ)eimφ
2π
For an atom or ion with one electron holds: Rlm (ρ) = Clm e−ρ/2 ρl L2l+1
n−l−1 (ρ)
with ρ = 2rZ/na0 with a0 = ε0 h2 /πme e2 . The Lji are the associated Laguere functions and the Plm
are the associated Legendre polynomials:
The fine structure in atomic spectra arises from this. With gS = 2 follows for the average magnetic
moment: M ⃗ av = −(e/2me )gℏJ,
⃗ where g is the Landé-factor:
⃗ · J⃗
S j(j + 1) + s(s + 1) − l(l + 1)
g =1+ =1+
J2 2j(j + 1)
For atoms with more than one electron the following limiting situations occur:
50 Physics Formulary by ir. J.C.A. Wevers
1. L − S coupling: for small atoms the electrostatic interaction is dominant and the state can be
characterized by L, S, J, mJ . J ∈ {|L − S|, ..., L + S − 1, L + S} and mJ ∈ {−J, ..., J − 1, J}. The
spectroscopic notation for this interaction is: 2S+1 LJ . 2S + 1 is the multiplicity of a multiplet.
2. j −j coupling: for larger atoms the electrostatic interaction is smaller than the Li ·si interaction
of an electron. The state is characterized by ji ...jn , J, mJ where only the ji of the not completely
filled subshells are to be taken into account.
The energy difference for larger atoms when placed in a magnetic field is: ∆E = gµB mJ B where
g is the Landé factor. For a transition between two singlet states the line splits in 3 parts, for
∆mJ = −1, 0 + 1. This results in the normal Zeeman effect. At higher S the line splits up in more
parts: the anomalous Zeeman effect.
Interaction with the spin of the nucleus gives the hyperfine structure.
where µ is the reduced mass of the system. The term ∼ A2 can usually be neglected, except for very
⃗ = B⃗ez it is given by e2 B 2 (x2 + y 2 )/8µ.
strong fields or macroscopic motions. For B
When a gauge transformation A ⃗′ = A ⃗ − ∇f , V ′ = V + ∂f /∂t is applied to the potentials the
wavefunction is also transformed according to ψ ′ = ψeiqef /ℏ with qe the charge of the particle.
Because f = f (x, t), this is called a local gauge transformation, in contrast with a global gauge
transformation which can always be applied.
(1) (2)
When cnk and En are being expanded into λ: cnk = λcnk + λ2 cnk + · · ·
(1) (2)
En = En0 + λEn + λ2 En + · · ·
Chapter 10: Quantum physics 51
(1)
and this is put into the Schrödinger equation the result is: En = ⟨ϕn |H1 |ϕn ⟩ and
⟨ϕm |H1 |ϕn ⟩
c(1)
nm = if m ̸= n. The second-order correction of the energy is then given by:
En0 − Em0
X | ⟨ϕk |H1 |ϕn ⟩ |2 X ⟨ϕk |λH1 |ϕn ⟩
En(2) = 0 . So to first order holds: ψn = ϕn + ϕk .
0
En − Ek En0 − Ek0
k̸=n k̸=n
In case the levels are degenerated the above does not hold. In that case an orthonormal set eigen-
functions ϕni is chosen for each level n, so that ⟨ϕmi |ϕnj ⟩ = δmn δij . Now ψ is expanded as:
X X (1) X
ψn = N (λ) αi ϕni + λ cnk βi ϕki + · · ·
i k̸=n i
0 (1) 0
Eni = Eni +λEni is approximated by Eni := En0 . Substitution in the Schrödinger equation and taking
(1)
|αi |2 = 1.
P P
dot product with ϕni gives: αi ⟨ϕnj |H1 |ϕni ⟩ = En αj . Normalization requires that
i i
Zt
i(En0 − Ek0 )t′
λ
follows: c(1)
n (t) = ⟨ϕn |V (t′ )|ϕk ⟩ exp dt′
iℏ ℏ
0
1. Particles with a half-odd integer spin (Fermions): ψtotal must be antisymmetric w.r.t. inter-
change of the coordinates (spatial and spin) of each pair of particles. The Pauli principle results
from this: two Fermions cannot exist in an identical state because then ψtotal = 0.
2. Particles with an integer spin (Bosons): ψtotal must be symmetric w.r.t. interchange of the
coordinates (spatial and spin) of each pair of particles.
For a system of two electrons there are 2 possibilities for the spatial wavefunction. When a and b are
the quantum numbers of electron 1 and 2 holds:
ψS (1, 2) = ψa (1)ψb (2) + ψa (2)ψb (1) , ψA (1, 2) = ψa (1)ψb (2) − ψa (2)ψb (1)
Because the particles do not approach each other closely the repulsion energy at ψA in this state is
smaller. The following spin wavefunctions are possible:
√
χA = 21 2[χ+ (1)χ− (2) − χ+ (2)χ− (1)] ms = 0
χ+√(1)χ+ (2) ms = +1
1
χS = 2[χ + (1)χ− (2) + χ+ (2)χ− (1)] ms = 0
2
χ− (1)χ− (2) ms = −1
Because the total wavefunction must be antisymmetric it follows: ψtotal = ψS χA or ψtotal = ψA χS .
52 Physics Formulary by ir. J.C.A. Wevers
10.15.2 Molecules
The wavefunctions of atom a and b are ϕa and ϕb . If the 2 atoms approach each other there are
two possibilities:
√ the total wavefunction approaches the bonding function with lower √ total energy
ψB = 12 2(ϕa +ϕb ) or approaches the anti-bonding function with higher energy ψAB = 12 2(ϕa −ϕb ).
If a molecular-orbital is symmetric w.r.t. the connecting axis, like a combination of two s-orbitals it
is called a σ-orbital, otherwise a π-orbital, like the combination of two p-orbitals along two axes.
⟨ψ|H|ψ⟩
The energy of a system is: E = .
⟨ψ|ψ⟩
The energy calculated with this method is always higher than the real energy if ψ is only an ap-
proximation for the solutions of Hψ = Eψ. Also, if there are more functions to be chosen, the
function
P which gives the lowest energy is the best approximation. Applying this to the function
ψ = ci ϕi one finds: (Hij − ESij )ci = 0. This equation has only solutions if the secular determinant
|Hij − ESij | = 0. Here, Hij = ⟨ϕi |H|ϕj ⟩ and Sij = ⟨ϕi |ϕj ⟩. αi := Hii is the Coulomb integral and
βij := Hij the exchange integral. Sii = 1 and Sij is the overlap integral.
The first approximation in the molecular-orbital theory is to place both electrons of a chemical bond in
the bonding orbital: ψ(1, 2) = ψB (1)ψB (2). This results in a large electron density between the nuclei
and therefore a repulsion. A better approximation is: ψ(1, 2) = C1 ψB (1)ψB (2) + C2 ψAB (1)ψAB (2),
with C1 = 1 and C2 ≈ 0.6.
In some atoms, such as C, it is energetical more suitable to form orbitals which are a linear combination
of the s, p and d states. There are three ways of hybridization in C:
√
1. SP-hybridization: ψsp = 21 2(ψ2s ± ψ2pz ). There are 2 hybrid orbitals which are placed on one
line under 180◦ . Further the 2px and 2py orbitals remain.
√ p √ √
2. SP2 hybridization: ψsp2 = ψ2s / 3+c1 ψ2pz +c2 ψ2py , where (c1 , c2 ) ∈ {( 2/3, 0), (−1/ 6, 1/ 2)
√ √
, (−1/ 6, −1/ 2)}. The 3 SP2 orbitals lay in one plane, with symmetry axes which are at an
angle of 120◦ .
3. SP3 hybridization: ψsp3 = 21 (ψ2s ± ψ2pz ± ψ2py ± ψ2px ). The 4 SP3 orbitals form a tetraheder
with the symmetry axes at an angle of 109◦ 28′ .
where ρlk = c∗k cl . ρ is hermitian, with Tr(ρ) = 1. Further holds ρ = ri |ψi ⟩⟨ψi |. The probability to
P
find eigenvalue an when measuring A is given by ρnn if one uses a basis of eigenvectors of A for {ϕk }.
For the time-dependence holds (in the Schrödinger image operators are not explicitly time-dependent):
dρ
iℏ = [H, ρ]
dt
Chapter 10: Quantum physics 53
For a macroscopic system in equilibrium holds [H, ρ] = 0. If all quantumstates with the same energy
are equally probable: Pi = P (Ei ), one can obtain the distribution:
e−En /kT X
Pn (E) = ρnn = with the state sum Z = e−En /kT
Z n
The thermodynamic quantities are related to these definitions as follows: F = −kT ln(Z), U = ⟨H⟩ =
P ∂ P
pn En = − ln(Z), S = −k Pn ln(Pn ). For a mixed state of M orthonormal quantum states
n ∂kT n
with probability 1/M follows: S = k ln(M ).
The distribution function for the internal states for a system in thermal equilibrium is the most
probable function. This function can be found by taking the maximum of the function P which gives
the number of states with Stirling’s equation: ln(n!) ≈ n ln(n) − n, and the conditions nk = N and
P k
nk Wk = W . For identical, indistinguishable particles which obey the Pauli exclusion principle the
k
possible number of states is given by:
Y gk !
P =
nk !(gk − nk )!
k
This results in the Fermi-Dirac statistics. For indistinguishable particles which do not obey the
exclusion principle the possible number of states is given by:
Y g nk
k
P = N!
nk !
k
This results in the Bose-Einstein statistics. So the distribution functions which explain how particles
are distributed over the different one-particle states k which are each gk -fold degenerate depend on
the spin of the particles. They are given by:
N gk
1. Fermi-Dirac statistics: integer spin. nk ∈ {0, 1}, nk =
P Z g exp((E k − µ)/kT ) + 1
with ln(Zg ) = gk ln[1 + exp((Ei − µ)/kT )].
N gk
2. Bose-Einstein statistics: half odd-integer spin. nk ∈ IN , nk =
P Z g exp((E k − µ)/kT ) − 1
with ln(Zg ) = − gk ln[1 − exp((Ei − µ)/kT )].
Here,
P Zg is the large-canonical state sum and µ the chemical potential. It is found by demanding
nk = N , and for it holds: lim µ = EF , the Fermi-energy. N is the total number of particles. The
T →0
Maxwell-Boltzmann distribution can be derived from this in the limit Ek − µ ≫ kT :
N Ek X Ek
nk = exp − with Z = gk exp −
Z kT kT
k
With the Fermi-energy, the Fermi-Dirac and Bose-Einstein statistics can be written as:
gk
1. Fermi-Dirac statistics: nk = .
exp((Ek − EF )/kT ) + 1
gk
2. Bose-Einstein statistics: nk = .
exp((Ek − EF )/kT ) − 1
Chapter 11
Plasma physics
11.1 Introduction
ne
The degree of ionization α of a plasma is defined by: α =
ne + n0
where ne is the electron density and n0 the density of the neutrals. If a plasma contains also negative
charged ions α is not well defined.
The probability that a test particle collides with another is given by dP = nσdx where σ is the cross
section. The collision frequency νc = 1/τc = nσv. The mean free path is given by λv = 1/nσ. The
rate coefficient K is defined by K = ⟨σv⟩. The number of collisions per unit of time and volume
between particles of kind 1 and 2 is given by n1 n2 ⟨σv⟩ = Kn1 n2 .
The potential of an electron is given by:
s r
−e r ε0 kTe Ti ε0 kTe
V (r) = exp − with λD = ≈
4πε0 r λD e2 (ne Ti + ni Te ) ne e2
because charge is shielded in a plasma. Here, λD is the Debye length. For distances < λD the
plasma cannot be assumed to be quasi-neutral. Deviations of charge neutrality by thermic motion
are compensated by oscillations with frequency
s
ne e2
ωpe =
me ε0
The distance of closest approximation when two equal charged particles collide for a deviation of π/2
−1/3
is 2b0 = e2 /(4πε0 21 mv 2 ). A “neat” plasma is defined as a plasma for which holds: b0 < ne ≪
λD ≪ Lp . Here Lp := |ne /∇ne | is the gradient length of the plasma.
11.2 Transport
1
Relaxation times are defined as τ = 1/νc . Starting with σm = 4πb20 ln(ΛC ) and with 2 mv
2
= kT it
can be found that: √ √
4πε2 m2 v 3 8 2πε20 m(kT )3/2
τm = 4 0 =
ne ln(ΛC ) ne4 ln(ΛC )
For momentum transfer between electrons and ions holds for a Maxwellian velocity distribution:
√ √ √ √
6π 3ε20 me (kTe )3/2 6π 3ε20 mi (kTi )3/2
τee = ≈ τei , τii =
ne e4 ln(ΛC ) ni e4 ln(ΛC )
The energy relaxation times for identical particles are equal to the momentum relaxation times.
Because for e-i collisions the energy p transfer is only ∼ 2me /mi this is a slow process. Approximately
holds: τee : τei : τie : τieE = 1 : 1 : mi /me : mi /me .
The relaxation for e-o interaction is much more complicated. For T > 10 eV holds approximately:
−2/5
σeo = 10−17 ve , for lower energies this can be a factor 10 lower.
Chapter 11: Plasma physics 55
The diffusion coefficient D is defined by means of the flux Γ by ⃗Γ = n⃗vdiff = −D∇n. The equation of
continuity is ∂t n + ∇(nvdiff ) = 0 ⇒ ∂t n = D∇2 n. One finds that D = 31 λv v. A rough estimate gives
τD = Lp /D = L2p τc /λ2v . For magnetized plasma’s λv must be replaced with the cyclotron radius. In
electrical fields also holds J⃗ = neµE
⃗ = e(ne µe + ni µi )E
⃗ with µ = e/mνc the mobility of the particles.
The Einstein ratio is:
D kT
=
µ e
Because a plasma is electrically neutral electrons and ions are strongly coupled and they don’t diffuse
independent. The coefficient of ambipolar diffusion Damb is defined by ⃗Γ = ⃗Γi = ⃗Γe = −Damb ∇ne,i .
From this follows that
kTe /e − kTi /e kTe µi
Damb = ≈
1/µe − 1/µi e
In an external magnetic field B0 particles will move in spiral orbits with cyclotron radius ρ = mv/eB0
and with cyclotron frequency Ω = B0 e/m. The helical orbit is perturbed by collisions. A plasma is
called magnetized if λv > ρe,i . So the electrons are magnetized if
√
ρe me e3 ne ln(ΛC )
= √ 2 <1
λee 6π 3ε0 (kTe )3/2 B0
Magnetization of only the electrons is sufficient to confine the plasma reasonable because they are
coupled to the ions by charge neutrality. In case of magnetic confinement holds: ∇p = J⃗ × B. ⃗
Combined with the two stationary Maxwell equations for the B-field these form the ideal magneto-
hydrodynamic equations. For a uniform B-field holds: p = nkT = B 2 /2µ0 .
If both magnetic and electric fields are present electrons and ions will move in the same direction.
⃗ = Er ⃗er + Ez ⃗ez and B
If E ⃗ = Bz ⃗ez the E ⃗ ×B ⃗ drift results in a velocity ⃗u = (E
⃗ ×B
⃗ )/B 2 and the
˙
velocity in the r, φ plane is ṙ(r, φ, t) = ⃗u + ρ
⃗(t).
Because the influence of a particle vanishes at r = λD holds: σ = π(λ2D − b20 ). Because dp = d(mv) =
mv0 (1 − cos χ) a cross section related to momentum transfer σm is given by:
ln(v 4 )
Z
1 λD
σm = (1 − cos χ)I(Ω)dΩ = 4πb20 ln 1 = 4πb20 ln := 4πb20 ln(ΛC ) ∼
sin( 2 χmin ) b0 v4
where ln(ΛC ) is the Coulomb-logarithm. For this quantity holds: ΛC = λD /b0 = 9n(λD ).
m2 sin(χ)
tan(θ) =
m1 + m2 cos(χ)
∆f 2v
= sin( 21 χ)
f c
This gives for the scattered energy Escat ∼ nλ40 /(λ2 − λ20 )2 with n the density. If λ ≫ λ0 it is
called Rayleigh scattering. Thomson sccattering is a limit of Compton scattering, which is given by
λ′ − λ = λC (1 − cos χ) with λC = h/mc and cannot be used any more if relativistic effects become
important.
Chapter 11: Plasma physics 57
8πhν 3 2πn √
1 E
ρ(ν, T )dν = dν , N (E, T )dE = E exp − dE
c3 exp(hν/kT ) − 1 (πkT )3/2 kT
“Detailed balancing” means that the number of reactions in one direction equals the number of
reactions in the opposite direction because both processes have equal probability if one corrects for
the used phase space. For the reaction
→
X X
Xforward ← Xback
forward back
nx h3
η̂x = e−Ekin /kT
gx (2πmx kT )3/2
where g is the statistical weight of the state and n/g := η. For electrons holds g = 2, for excited
states usually holds g = 2j + 1 = 2n2 .
With this one finds for the Boltzmann balance, Xp + e− ← → X1 + e− + (E1p ):
nB
p gp Ep − E1
= exp
n1 g1 kTe
→ X+ + 2e− :
And for the Saha balance, Xp + e− + (Epi ) ← 1
nSp n+ h3
1 ne Epi
= + exp
gp g1 ge (2πme kTe )3/2 kTe
Because the number of particles on the left-hand side and right-hand side of the equation is different,
a factor g/Ve remains. This factor causes the Saha-jump.
From microscopic reversibility one can derive that for the rate coefficients K(p, q, T ) := ⟨σv⟩pq holds:
gp ∆Epq
K(q, p, T ) = K(p, q, T ) exp
gq kT
m1 m2 (v1 − v2 )2
E=
2(m1 + m2 )
Some types of inelastic collisions important for plasma physics are:
→ Aq + e −
1. Excitation: Ap + e− ←
→ Ap + hf
2. Decay: Aq ←
58 Physics Formulary by ir. J.C.A. Wevers
→ A+ + 2e−
3. Ionisation and 3-particles recombination: Ap + e− ←
→ Ap + hf
4. radiative recombination: A+ + e− ←
5. Stimulated emission: Aq + hf → Ap + 2hf
→ AB+ + e−
6. Associative ionisation: A∗∗ + B ←
→ Ar+ + Ne + e−
7. Penning ionisation: b.v. Ne∗ + Ar ←
→ A + B+
8. Charge transfer: A+ + B ←
→ A + A+
9. Resonant charge transfer: A+ + A ←
πZ 2 e4 ∆Eq,q+1
Then follows for the transition p → q: σpq (E) =
(4πε0 )2 E(∆E)2pq
1 1 1.25βE
For ionization from state p holds to a good approximation: σp = 4πa20 Ry − ln
Ep E Ep
A[1 − B ln(E)]2
For resonant charge transfer holds: σex =
1 + CE 3.3
11.6 Radiation
In equilibrium holds for radiation processes:
Here, Apq is the matrix element of the transition p → q, and is given by:
8π 2 e2 ν 3 |rpq |2
Apq = with rpq = ⟨ψp |⃗r |ψq ⟩
3ℏε0 c3
c3 Apq Bpq gq
Bpq = and =
8πhν 3 Bqp gp
The natural line width is usually ≪ than the broadening due to the following two mechanisms:
Chapter 11: Plasma physics 59
The natural broadening and the Stark broadening result in a Lorentz profile of a spectral line: kν =
1 1 2 2
2 k0 ∆νL /[( 2 ∆νL ) + (ν − ν0 ) ]. The total line shape is a convolution of the Gauss- and Lorentz profile
and is called a Voigt profile.
The number of transitions p → q is given by np Bpq ρ and by np nhf ⟨σa c⟩ = np (ρdν/hν)σa c where dν
is the line width. Then follows for the cross section of absorption processes: σa = Bpq hν/cdν.
The background radiation in a plasma originates from two processes:
1. Free-Bound radiation, originating from radiative recombination. The emission is given by:
C1 zi ni ne hc
εf b = 2 √ 1 − exp − ξf b (λ, Te )
λ kTe λkTe
with C1 = 1.63 · 10−43 Wm4 K1/2 sr−1 and ξ the Biberman factor.
2. Free-free radiation, originating from the acceleration of particles in the EM-field of other parti-
cles:
C1 zi ni ne hc
εf f = 2 √ exp − ξf f (λ, Te )
λ kTe λkTe
F (⃗r, ⃗v , t) = Fr (⃗r, t) · Fv (⃗v , t) = F1 (x, t)F2 (y, t)F3 (z, t)F4 (vx , t)F5 (vy , t)F6 (vz , t)
dF ∂F ∂F
Then the BTE is: = + ∇r · (F ⃗v ) + ∇v · (F⃗a ) =
dt ∂t ∂t coll−rad
Assuming that v does not depend on r and ai does not depend on vi , holds ∇r · (F ⃗v ) = ⃗v · ∇F
and ∇v · (F⃗a ) = ⃗a · ∇v F . This is also true in magnetic fields because ∂ai /∂xi = 0. The velocity
R is
separated
R in a thermal velocity ⃗
v t and a drift velocity w.
⃗ The total density is given by n = F d⃗v
and ⃗v F d⃗v = nw.
⃗
The balance equations can be derived by means of the moment method:
Z
∂n ∂n
1. Mass balance: (BTE)d⃗v ⇒ + ∇ · (nw)
⃗ =
∂t ∂t cr
Z
dw⃗
2. Momentum balance: (BTE)m⃗v d⃗v ⇒ mn + ∇T′ + ∇p = mn ⟨⃗a ⟩ + R
⃗
dt
Z
3 dp 5
3. Energy balance: (BTE)mv 2 d⃗v ⇒ + p∇ · w⃗ + ∇ · ⃗q = Q
2 dt 2
60 Physics Formulary by ir. J.C.A. Wevers
⃗ i(kz−ωt)
⃗ = B0⃗ez +B̂e
For disturbances in the z-direction in a cold, homogeneous, magnetized plasma: B
and n = n0 + n̂ei(kz−ωt) (external E fields are screened) follows, with the definitions α = ωp /ω and
β = Ω/ω and ωp2 = ωpi 2 2
+ ωpe :
1 −iβs
0
X 1 − βs2 1 − βs2
J⃗ = ⃗ ⃗ , with ⃗
⃗σ E ⃗σ = iε0 ω αs2 iβs 1
0
1 − βs2 1 − βs2
s
0 0 1
where the sum is taken over particle species s. The dielectric tensor E, with property:
⃗k · (⃗E⃗ · E)
⃗ =0
⃗ ⃗
is given by E⃗ = I⃗ − ⃗
⃗σ /iε0 ω.
X αs2 X α 2 βs
s
X
With the definitions S = 1 − , D = , P = 1 − αs2
s
1 − βs2 s
1 − βs2 s
follows:
S −iD 0
⃗⃗
E= iD S 0
0 0 P
The eigenvalues
√ of this √ hermitian matrix are R = S + D, L = S − D, λ3 = P , with eigenvectors
⃗er = 21 2(1, i, 0), ⃗el = 21 2(1, −i, 0) and ⃗e3 = (0, 0, 1). ⃗er is connected with a right rotating field for
which iEx /Ey = 1 and ⃗el is connected with a left rotating field for which iEx /Ey = −1. When k
makes an angle θ with B ⃗ one finds:
P (n2 − R)(n2 − L)
tan2 (θ) =
S(n2 − RL/S)(n2 − P )
where n is the refractive index. From this the following solutions can be obtained:
A. θ = 0: transmission in the z-direction.
1. P = 0: Ex = Ey = 0. This describes a longitudinal linear polarized wave.
2. n2 = L: a left, circular polarized wave.
3. n2 = R: a right, circular polarized wave.
B. θ = π/2: transmission ⊥ the B-field.
1. n2 = P : the ordinary mode: Ex = Ey = 0. This is a transversal linear polarized wave.
2. n2 = RL/S: the extraordinary mode: iEx /Ey = −D/S, an elliptical polarized wave.
Resonance frequencies are frequencies for which n2 → ∞, so vf = 0. For these holds: tan(θ) = −P/S.
For R → ∞ this gives the electron cyclotron resonance frequency ω = Ωe , for L → ∞ the ion cyclotron
resonance frequency ω = Ωi and for S = 0 holds for the extraordinary mode:
mi Ω2i m2i Ω2i Ω2i
2
α 1− = 1− 2 2 1− 2
me ω 2 me ω ω
Cut-off frequencies are frequencies for which n2 = 0, so vf → ∞. For these holds: P = 0 or R = 0 or
L = 0.
In the case that β 2 ≫ 1 one finds Alfvén waves propagating parallel to the field lines. With the
Alfvén velocity
Ωe Ωi 2
vA = 2 2 c
ωpe + ωpi
p
follows: n = 1 + c/vA , and in case vA ≪ c: ω = kvA .
Chapter 12
ne (⃗r + T⃗ ) = ne (⃗r )
This periodicity is suitable to use Fourier analysis: n(⃗r ) is expanded as:
X
n(⃗r ) = nG exp(iG ⃗ · ⃗r )
G
with ZZZ
1 ⃗ · ⃗r )dV
nG = n(⃗r ) exp(−iG
Vcell
cell
⃗ is the reciprocal lattice vector. If G
G ⃗ is written as G ⃗ = v1⃗b1 + v2⃗b2 + v3⃗b3 with vi ∈ IN , it follows for
⃗
the vectors bi , cyclically:
⃗bi = 2π ⃗ai+1 × ⃗ai+2
⃗ai · (⃗ai+1 × ⃗ai+2 )
⃗
The set of G-vectors determines the Röntgen diffractions: a maximum in the reflected radiation occurs
if: ∆k = G with ∆⃗k = ⃗k − ⃗k ′ . So: 2⃗k · G
⃗ ⃗ ⃗ = G2 . From this follows for parallel lattice planes (Bragg
reflection) that for the maxima holds: 2d sin(θ) = nλ.
The Brillouin zone is defined as a Wigner-Seitz cell in the reciprocal lattice.
d2 us
Fs = M = C(us+1 − us ) + C(us−1 − us )
dt2
Assuming that all solutions have the same time-dependence exp(−iωt) this results in:
This gives: us = exp(iKsa). Substituting the later two equations in the fist results in a system of
linear equations, which has only a solution if their determinant is 0. This gives:
4C
ω2 = sin2 ( 12 Ka)
M
Only vibrations with a wavelength within the first Brillouin Zone have a physical significance. This
requires that −π < Ka ≤ π.
12.3.3 Phonons
The quantum mechanical excitation of a crystal vibration with an energy ℏω is called a phonon.
Phonons can be viewed as quasi-particles: with collisions, they behave as particles with momentum
ℏK. Their total momentum is 0. When they collide, their momentum need not be conserved: for a
normal process holds: K1 + K2 = K3 , for an umklapp process holds: K1 + K2 = K3 + G. Because
phonons have no spin they behave like bosons.
64 Physics Formulary by ir. J.C.A. Wevers
In the Einstein model for the thermal heat capacity one considers only phonons at one frequency, an
approximation for optical phonons.
Chapter 12: Solid state physics 65
0 hhhh - H
hhhh diamagnetism
hhh h
12.4.1 Dielectrics
The quantum mechanical origin of diamagnetism is the Larmorprecession of the spin of the electron.
Starting with a circular electron orbit in an atom with two electrons, there is a Coulomb force Fc
and a magnetic force on each electron. If the magnetic part of the force is not strong enough to
significantly deform the orbit holds:
s 2
2 Fc (r) eB 2 eB eB eB
ω = ± ω = ω0 ± (ω0 + δ) ⇒ ω = ω0 ± + · · · ≈ ω0 ± = ω0 ± ωL
mr m m 2m 2m
Here, ωL is the Larmor frequency. One electron is accelerated, the other decelerated. Hence there
is a net circular current which results
in
a magnetic
moment µ ⃗ . The circular current is given by
I = −ZeωL /2π, and ⟨µ⟩ = IA = Iπ ρ2 = 32 Iπ r2 . If N is the number of atoms in the crystal it
⃗ =µ
follows for the susceptibility, with M ⃗N:
µ0 M µ0 N Ze2
2
χ= =− r
B 6m
12.4.2 Paramagnetism
Starting with the splitting of energy levels in a weak magnetic field: ∆Um − µ⃗ ·B ⃗ = mJ gµB B,
P withPa distribution fm ∼ exp(−∆Um /kT ),Pone finds P
and for the average magnetic moment ⟨µ⟩ =
J = 2J + 1 and m2J = 32 J(J + 1)(J +
P
fm µ/ fm . After linearization and because mJ = 0,
1
2 ) it follows that:
µ0 M µ0 N ⟨µ⟩ µ0 J(J + 1)g 2 µ2B N
χp = = =
B B 3kT
This is the Curie law, χp ∼ 1/T .
12.4.3 Ferromagnetism
A ferromagnet behaves like a paramagnet above a critical temperature Tc . To describe ferromagnetism
⃗ E = λµ0 M
a field BE parallel with M is postulated: B ⃗ . From there the treatment is analogous to the
paramagnetic case:
C
µ0 M = χp (Ba + BE ) = χp (Ba + λµ0 M ) = µ0 1 − λ M
T
µ0 M C
From this follows for a ferromagnet: χF = = which is Weiss-Curie’s law.
Ba T − Tc
If BE is estimated this way it results in values of about 1000 T. This is clearly unrealistic and suggests
another mechanism. A quantum mechanical approach from Heisenberg postulates an interaction
between two neighbouring atoms: U = −2J S ⃗i · S
⃗j ≡ −⃗µ·B⃗ E . J is an overlap integral given by:
J = 3kTc /2zS(S + 1), with z the number of neighbours. A distinction between 2 cases can now be
made:
66 Physics Formulary by ir. J.C.A. Wevers
⃗p · (S
U = −2J S ⃗p−1 + S
⃗p+1 ) ≈ µ ⃗p
⃗p · B ⃗ p = −2J (S
with B ⃗p−1 + S
⃗p+1 )
gµB
⃗
dS 2J ⃗
The equation of motion for the magnons becomes: = ⃗p−1 + S
Sp × (S ⃗p+1 )
dt ℏ
⃗p =
From here the treatment is analogous to phonons: postulate traveling waves of the type S
⃗u exp(i(pka − ωt)). This results in a system of linear equations with solution:
ℏω = 4JS(1 − cos(ka))
The probability density |ψ(+)|2 is high near the atoms of the lattice and low in between. The
probability density |ψ(−)|2 is low near the atoms of the lattice and high in between. Hence the
energy of ψ(+) is also lower than the energy of ψ)(−). Suppose that U (x) = U cos(2πx/a), than the
bandgap is given by:
Z1
Egap = U (x) |ψ(+)|2 − |ψ(−)|2 dx = U
12.7 Semiconductors
The band structures and the transitions between them of direct and indirect semiconductors are
shown in the figures below. Here it is assumed that the momentum of the absorbed photon can be
neglected. For an indirect semiconductor a transition from the valence- to the conduction band is
also possible if the energy of the absorbed photon is smaller than the band gap: then, also a phonon
is absorbed.
E conduction E
6 band 6
•
•
6
6
ωg Ω ω
◦ ◦
absorption absorption
6 6
.
.....
..
....
. .
....
-E
.. - E
ℏωg Eg + ℏΩ
Direct semiconductor Indirect semiconductor
So indirect semiconductors, like Si and Ge, cannot emit any light and are therefore not usable to
fabricate lasers. When light is absorbed holds: ⃗kh = −⃗ke , Eh (⃗kh ) = −Ee (⃗ke ), ⃗vh = ⃗ve and mh = −m∗e
if the conduction band and the valence band have the same structure.
Instead of the normal electron mass one has to use the effective mass within a lattice. It is defined
by:
2 −1
∗ F dp/dt dK 2 d E
m = = =ℏ =ℏ
a dvg /dt dvg dk 2
with E = ℏω and vg = dω/dk and p = ℏk.
With the distribution function fe (E) ≈ exp((µ − E)/kT ) for the electrons and fh (E) = 1 − fe (E) for
the holes the density of states is given by:
∗ 3/2 p
1 2m
D(E) = 2
E − Ec
2π ℏ2
with Ec the energy at the edge of the conductance band. From this follows for the concentrations of
the holes p and the electrons n:
Z∞ 3/2
m∗ kT
µ − Ec
n= De (E)fe (E)dE = 2 exp
2πℏ2 kT
Ec
3
kT p Eg
For the product np follows: np = 4 m∗e mh exp −
2πℏ2 kT
For an intrinsic (no impurities) semiconductor holds: ni = pi , for a n − type holds: n > p and in a
p − type holds: n < p.
An exciton is a bound electron-hole pair, rotating on each other as in positronium. The excitation
energy of an exciton is smaller than the bandgap because the energy of an exciton is lower than the
energy of a free electron and a free hole. This causes a peak in the absorption just under Eg .
12.8 Superconductivity
12.8.1 Description
A superconductor is characterized by a zero resistivity if certain quantities are smaller than some
critical values: T < Tc , I < Ic and H < Hc . The BCS-model predicts for the transition temperature
Tc :
−1
Tc = 1.14ΘD exp
U D(EF )
while experiments find for Hc approximately:
T2
Hc (T ) ≈ Hc (Tc ) 1 − 2 .
Tc
Chapter 12: Solid state physics 69
··
··
··
··
··
··
··
-H · - H
Hc Hc1 Hc2
Type I Type II
The transition to a superconducting state is a second order thermodynamic state transition. This
means that there is a twist in the T − S diagram and a discontinuity in the CX − T diagram.
∂ψ1 ∂ψ2
iℏ = ℏT ψ2 , iℏ = ℏT ψ1
∂t ∂t
ℏT is the effect of the coupling of the electrons, or the transfer interaction through the insulator. The
√ √
electron wavefunctions are written as √ ψ1 = n1 exp(iθ1 ) and ψ2 = n2 exp(iθ2 ). Because a Cooper
pair exist of two electrons holds: ψ ∼ n. From this follows, if n1 ≈ n2 :
∂ψ1 ∂ψ2
iℏ = ℏT ψ2 − eV ψ1 and iℏ = ℏT ψ1 + eV ψ2
∂t ∂t
2eV t
This gives: J = J0 sin θ2 − θ1 − .
ℏ
Hence there is an oscillation with ω = 2eV /ℏ.
−A⃗ −B⃗
J⃗ = or rotJ⃗ =
µ0 λ2L µ0 λ2L
⃗ = B/λ
⃗ 2.
p
where λL = ε0 mc2 /nq 2 . From this follows: ∇2 B L
⃗
The Meissner effect is the solution of this equation: B(x) = B0 exp(−x/λL ). Magnetic fields within
a superconductor drop exponentially.
∞ Z∞ Z∞
X
n 1 1
And, when (−1) = 2 follows: sin(px)dx = cos(px)dx = .
n=0
p
0 0
Chapter 13
Theory of groups
13.1 Introduction
13.1.1 Definition of a group
G is a group for the operation • if:
1. ∀A,B∈G ⇒ A • B ∈ G: G is closed.
If also holds:
5. ∀A,B∈G ⇒ A • B = B • A the group is called Abelian or commutative.
A conjugacy class is the maximum collection of conjugated elements. Each group can be split up in
conjugacy classes. Some theorems:
E is a class itself: for each other element in this class would hold: A = XEX −1 = E.
E is the only class which is also a subgroup because all other classes have no unit element.
The physical interpretation of classes: elements of a group are usually symmetry operations which
map a symmetrical object into itself. Elements of one class are then the same kind of operations.
The opposite need not to be true.
72 Physics Formulary by ir. J.C.A. Wevers
13.2.3 Character
The character of a representation is given by the trace of the matrix and is therefore invariant for
base transformations: χ(j) (R) = Tr(Γ(j) (R))
X
Also holds, with Nk the number of elements in a conjugacy class: χ(i)∗ (Ck )χ(j) (Ck )Nk = hδij
k
n
X
Theorem: ℓ2i = h
i=1
Chapter 13: Theory of groups 73
where Γ(n) is an irreducible, unitary representation of the symmetry group G of the system. Each
n corresponds with another energy level. One can purely mathematical derive irreducible represen-
tations of a symmetry group and label the energy levels with a quantum number this way. A fixed
(n)
choice of Γ(n) (R) defines the base functions ψν . This way one can also label each separate base
function with a quantum number.
Particle in a periodical potential: the symmetry operation is a cyclic group: note the operator
describing one translation over one unit as A. Then: G = {A, A2 , A3 , . . . , Ah = E}.
The group is Abelian so all irreducible representations are one-dimensional. For 0 ≤ p ≤ h−1 follows:
ℓn1 = dim(Γ(n1 ) )
ℓn ℓn2 = dim(Γ(n2 ) )
ℓn3 = dim(Γ(n3 ) )
Spectrum H0 Spectrum H
(n)
Theorem: The set of ℓn degenerated eigenfunctions ψν with energy En is a basis for an ℓn -
dimensional irreducible representation Γ(n) of the symmetry group.
In general the space D(1) ⊗ D(2) can be split up in a number of invariant subspaces:
X
Γ(1) ⊗ Γ(2) = ni Γ(i)
i
A useful tool for this reduction is that for the characters hold:
X
χ(1) (R)χ(2) (R) = ni χ(i) (R)
i
(j)
X
and the inverse transformation by: φ(i)
κ φλ = φ(aκ)
µ (akµ|iκjλ)
akµ
(j)
If Γ(j) is irreducible they are called irreducible tensor operators A(j) with components Aκ .
Chapter 13: Theory of groups 75
P (j)
An operator can also be decomposed into symmetry types: A = ak , with:
jk
!
ℓj X (j)∗
a(j)
κ = Γκκ (R) (PR APR−1 )
h
R∈G
Theorem: Matrix elements Hij of the operator H which is invariant under ∀A∈G , are 0 between
states which transform according to non-equivalent irreducible unitary representations or according
(i) (i)
to different rows of such a representation. Further ⟨φκ |H|ψκ ⟩ is independent of κ. For H = 1 this
becomes the previous theorem.
This is applied in quantum mechanics in perturbation theory and variational calculus. Here one tries
(i)
to diagonalize H. Solutions can be found within each category of functions φκ with common i and
κ: H is already diagonal in categories as a whole.
Perturbation calculus can be applied independent within each category. With variational calculus
the try function can be chosen within a separate category because the exact eigenfunctions transform
according to a row of an irreducible representation.
(i)
⟨φλ |A(j) (k) (i) (j)
κ |ψµ ⟩ = (iλ|jκkµ)⟨φ ∥A ∥ψ
(k)
⟩
This theorem can be used to determine selection rules: the probability of a dipole transition is given
by (⃗ϵ is the direction of polarization of the radiation):
8π 2 e2 f 3 |r12 |2
PD = with r12 = ⟨l2 m2 |⃗ϵ · ⃗r |l1 m1 ⟩
3ℏε0 c3
Further it can be used to determine intensity ratios: if there is only one value of a the ratio of
the matrix elements are the Clebsch-Gordan coefficients. For more a-values relations between the
intensity ratios can be stated. However, the intensity ratios are also dependent on the occupation of
the atomic energy levels.
If α and θ are very small and are expanded to second order, and the corresponding terms are put
equal with ⃗n · J⃗ = Jx cos θ + Jz sin θ, it follows from the αθ term: Jx Jy − Jy Jx = iℏJz .
For the translation group holds: g(⃗a) = constant = g(⃗0 ) because g(a⃗n )d⃗a′ = g(⃗0 )d⃗a and d⃗a′ = d⃗a.
This leads to the fundamental orthogonality theorem:
Z Z
(j) 1
Γ(i)∗
µν (R)Γαβ (R)dR = δ ij δ µα δνβ dR
ℓi
G G
In the performed derivation α is the rotational angle around the z-axis. This expression is valid for
all rotations over an angle α because the classes of SO(3) are rotations around the same angle around
an axis with an arbitrary orientation.
Via the fundamental orthogonality theorem for characters one obtains the following expression for
the density function (which is normalized so that g(0) = 1):
sin2 ( 21 α)
g(α) =
( 12 α)2
With this result one can see that the given representations of SO(3) are the only ones: the character
of another representation χ′ would have to be ⊥ to the already found ones, so χ′ (α) sin2 ( 21 α) = 0∀α
⇒ χ′ (α) = 0∀α. This is contradictory because the dimension of the representation is given by χ′ (0).
Because fermions have an half-odd integer spin the states ψsms with s = 21 and ms = ± 12 constitute
a 2-dim. space which is invariant under rotations. A problem arises for rotations over 2π:
ψ 12 ms → e−2πiSz /ℏ ψ 12 ms = e−2πims ψ 12 ms = −ψ 21 ms
78 Physics Formulary by ir. J.C.A. Wevers
However, in SO(3) holds: Rz,2π = E. So here holds E → ±II. Because observable quantities can
always be written as ⟨ϕ|ψ⟩ or ⟨ϕ|A|ψ⟩, and are bilinear in the states, they do not change sign if the
states do. If only one state changes sign the observable quantities do change.
The existence of these half-odd integer representations is connected with the topological properties
of SO(3): the group is two-fold coherent through the identification R0 = R2π = E.
D(j1 ) ⊗ D(j2 ) = D(j1 +j2 ) ⊕ D(j1 +j2 −1) ⊕ ... ⊕ D(|j1 −j2 |)
The states can be characterized by quantum numbers in two ways: with j1 , m1 , j2 , m2 and with
j1 , j2 , J, M . The Clebsch-Gordan
P coefficients, for SO(3) called the Wigner coefficients, can be chosen
real, so: ψj1 j2 JM = ψj1 m1 j2 m2 (j1 m1 j2 m2 |JM )
m
P1 m2
ψj1 m1 j2 m2 = ψj1 j2 JM (j1 m1 j2 m2 |JM )
JM
characters this means one can choose base operators which transform like Y1m (θ, φ). These turn
out to be the spherical components:
(1) 1 (1) (1) 1
A+1 = − √ (Ax + iAy ), A 0 = Az , A−1 = √ (Ax − iAy )
2 2
3. A cartesian tensor of rank 2: Tij is a quantity which transforms under rotations like Ui Vj , where
⃗ and V
U ⃗ are vectors. So Tij transforms like PR Tij P −1 = P Tkl Dki (R)Dlj (R), so like D(1) ⊗
R
kl
D(1) = D(2) ⊕ D(1) ⊕ D(0) . The 9 components can be split in 3 invariant subspaces with
dimension 1 (D(0) ), 3 (D(1) ) and 5 (D(2) ). The new base operators are:
⃗ ·V
I. Tr(T ) = Txx + Tyy + Tzz . This transforms as the scalar U ⃗ , so as D(0) .
II. The 3 antisymmetric components Az = 21 (Txy − Tyx ), etc. These transform as the vector
⃗ ×V
U ⃗ , so as D(1) .
III. The 5 independent components of the traceless, symmetric tensor S:
Sij = 12 (Tij + Tji ) − 13 δij Tr(T ). These transform as D(2) .
Chapter 13: Theory of groups 79
⟨αjm|A⃗ · J⃗ |αjm⟩
⃗
⟨αjm′ |A|αjm⟩ = ⟨αjm′ |J⃗ |αjm⟩
j(j + 1)ℏ2
as a linear combination of the 3 Pauli-matrices σi . So these matrices are a choice for the operators
⃗
of SU(2). One can write: SU(2)={exp(− 12 i⃗σ · Θ)}.
In abstraction, one can consider an isomorphic group where only the commutation rules are considered
to be known regarding the operators Ti : [T1 , T2 ] = iT3 , etc.
In elementary particle physics the Ti can be interpreted e.g. as the isospin operators. Elementary
particles can be classified in isospin-multiplets, these are the irreducible representations of SU(2).
The classification is:
⃗ ⃗
1. The isospin-singlet ≡ the identical representation: e−iT ·Θ = 1 ⇒ Ti = 0
2. The isospin-doublet ≡ the faithful representation of SU(2) on 2 × 2 matrices.
The group SU(3) has 8 free parameters. (The group SU(N ) has N 2 − 1 free parameters). The
Hermitian, traceless operators are 3 SU(2)-subgroups in the ⃗e1⃗e2 , ⃗e1⃗e3 and the ⃗e2⃗e3 plane. This gives
9 matrices, which are not all 9 linear independent. By taking a linear combination one gets 8 matrices.
8
∂ D ∂ X
In the Lagrange density for the colour force one has to substitute → := − Ti Aix
∂x Dx ∂x i=1
The terms of 3rd and 4th power in A show that the colour field interacts with itself.
Chapter 14
Nuclear physics
and A = Z + N , in the droplet or collective model of the nucleus the binding energy Ebind is given by:
energy. This is just the opposite for electrons, which is an indication that the L − S interaction
is not electromagnetical. The energy of a 3-dimensional harmonic oscillator is E = (N + 23 )ℏω.
N = nx + ny + nz = 2(n − 1) + l where n ≥ 1 is the main oscillator number. Because −l ≤ m ≤ l and
82 Physics Formulary by ir. J.C.A. Wevers
ms = ± 21 ℏ there are 2(2l + 1) substates which exist independently for protons and neutrons. This
gives rise to the so called magical numbers: nuclei where each state in the outermost level are filled
are particulary stable. This is the case if N or Z ∈ {2, 8, 20, 28, 50, 82, 126}.
l = 0 gives rise to monopole vibrations, density vibrations, which can be applied to the theory of
neutron
√ stars. l = 1 gives dipole vibrations, l = 2 quadrupole, with a2,0 = β cos γ and a2,±2 =
1
2 2β sin γ where β is the deformation factor and γ the shape parameter. The multipole moment
is given by µl = Zerl Ylm (θ, φ). The parity of the electric moment is ΠE = (−1)l , of the magnetic
moment ΠM = (−1)l+1 .
⃗L = e ⃗ ⃗ S = gS e S.
⃗
There are 2 contributions to the magnetic moment: M L and M
2mp 2mp
where gS is the spin-gyromagnetic ratio. For protons holds gS = 5.5855 and for neutrons gS =
−3.8263. The z-components of the magnetic moment are given by ML,z = µN ml and MS,z = gS µN mS .
⃗ = gI (e/2mp )I.
The resulting magnetic moment is related to the nuclear spin I according to M ⃗ The
z-component is then Mz = µN gI mI .
1. α-decay: the nucleus emits a He2+ nucleus. Because nucleons tend to order themselves in groups
of 2p+2n this can be considered as a tunneling of a He2+ nucleus through a potential barrier.
The tunnel probability P is
s Z
incoming amplitude 1
P = = e−2G with G = 2m [V (r) − E]dr
outgoing amplitude ℏ
3. Electron capture: here, a proton in the nucleus captures an electron (usually from the K-shell).
5. γ-decay: here the nucleus emits a high-energetic photon. The decay constant is given by
2l
P (l) Eγ Eγ R
λ= ∼ ∼ 10−4l
ℏω (ℏc)2 ℏc
where l is the quantum number for the angular momentum and P the radiated power. Usually
the decay constant of electric multipole moments is larger than the one of magnetic multipole
moments. The energy of the photon is Eγ = Ei − Ef − TR , with TR = Eγ2 /2mc2 the recoil
l i f
p is Π = Π · Π .
energy, which can usually be neglected. The parity of the emitted radiation
With I the quantum number of angular momentum of the nucleus, L = ℏ I(I + 1), holds the
following selection rule: |I⃗i − I⃗f | ≤ ∆l ≤ |I⃗i + I⃗f |.
sin(kr + δ0 )
If the potential is approximately rectangular holds: ψ0′ = C
kr
sin2 (δ0 ) 4π sin2 (δ0 )
Z
The cross section is then σ(θ) = so σ = σ(θ)dΩ =
k2 k2
ℏ2 k 2 /2m
At very low energies holds: sin2 (δ0 ) =
W0 + W
4π X 2
with W0 the depth of the potential well. At higher energies holds: σ = sin (δl )
k2
l
84 Physics Formulary by ir. J.C.A. Wevers
mk )c2 .
P
the total energy Q gained or required is given by Q = (m1 + m2 −
k>2
The minimal required kinetic energy T of P1 in the laboratory system to initialize the reaction is
P
m1 + m2 + mk
T = −Q
2m2
If Q < 0 there is a threshold energy.
The fluention is given by Φ = dN/dA. The flux is given by ϕ = dΦ/dt. The energy loss is defined
by Ψ = dW/dA, and the energy flux density ψ = dΨ/dt. The absorption coefficient is given by
µ = (dN/N )/dx. The mass absorption coefficient is given by µ/ϱ.
The radiation dose X is the amount of charge produced by the radiation per unit of mass, with unit
C/kg. An old unit is the Röntgen: 1Ro= 2.58 · 10−4 C/kg. With the energy-absorption coefficient
µE follows:
dQ eµE
X= = Ψ
dm Wϱ
where W is the energy required to disjoin an elementary charge.
The absorbed dose D is given by D = dEabs /dm, with unit Gy=J/kg. An old unit is the rad: 1
rad=0.01 Gy. The dose tempo is defined as Ḋ. It can be derived that
µE
D= Ψ
ϱ
The Kerma K is the amount of kinetic energy of secundary produced particles which is produced per
mass unit of the radiated object.
The equivalent dose H is a weight average of the absorbed dose per type of radiation, where for each
type radiation the effects on biological material is used for the weight factor. These weight factors
are called the quality factors. Their unit is Sv. H = QD.
P If the absorption is not equally distributed
also weight factors w per organ need to be used: H = wk Hk . For some types of radiation holds:
Radiation type Q
Röntgen, gamma radiation 1
β, electrons, mesons 1
Thermic neutrons 3 to 5
Fast neutrons 10 to 20
protons 10
α, fission products 20
Chapter 15
The coefficients c can be interpreted as follows: |cn1 n2 ··· |2 is the probability to find n1 particles
with momentum ⃗k1 , n2 particles with momentum ⃗k2 , etc., and ⟨Ψ(t)|Ψ(t)⟩ = |cni (t)|2 = 1. The
P
expansion of the states in time is described by the Schrödinger equation
d
i |Ψ(t)⟩ = H|Ψ(t)⟩
dt
where H = H0 + Hint . H0 is the Hamiltonian for free particles and keeps |cni (t)|2 constant, Hint is
the interaction Hamiltonian and can increase or decrease a c2 at the cost of others.
All operators which can change occupation numbers can be expanded in the a and a† operators. a is
the annihilation operator and a† the creation operator, and:
√
a(⃗ki )|n1 n2 · · · ni · · ·⟩ = ni |n1 n2 · · · ni − 1 · · ·⟩
† ⃗
√
a (ki )|n1 n2 · · · ni · · ·⟩ = ni + 1 |n1 n2 · · · ni + 1 · · ·⟩
Because the states are normalized holds a|0⟩ = 0 and a(⃗ki )a† (⃗ki )|ni ⟩ = ni |ni ⟩. So aa† is an occupation
number operator. The following commutation rules can be derived:
[a(⃗ki ), a(⃗kj )] = 0 ,
[a† (⃗ki ), a† (⃗kj )] = 0 , [a(⃗ki ), a† (⃗kj )] = δij
The conjugated field is, analogous to momentum in classical mechanics, defined as:
∂L
Πα (x) =
∂ Φ̇α
With this, the Hamilton density becomes H(x) = Πα Φ̇α − L(x).
Quantization of a classical field is analogous to quantization in point mass mechanics: the field
functions are considered as operators obeying certain commutation rules:
[Φα (⃗x), Φβ (⃗x ′ )] = 0 , [Πα (⃗x), Πβ (⃗x ′ )] = 0 , [Φα (⃗x), Πβ (⃗x ′ )] = iδαβ (⃗x − ⃗x ′ )
The field operators contain a volume V , which is used as normalization factor. Usually one can take
the limit V → ∞.
In general it holds that the term with e−ikx , the positive frequency part, is the creation part, and the
negative frequency part is the annihilation part.
the coefficients have to be each others hermitian conjugate because Φ is hermitian. Because Φ has
only one component this can be interpreted as a field describing a particle with spin zero. From this
follows that the commutation rules are given by [Φ(x), Φ(x′ )] = i∆(x − x′ ) with
Z
1 sin(ky) 3
∆(y) = d k
(2π)3 ωk
∆(y) is an odd function which is invariant for proper Lorentz transformations (no mirroring). This
is consistent with the previously found result [Φ(⃗x, t, Φ(⃗x ′ , t)] = 0. In general holds that ∆(y) = 0
outside the light cone. So the equations obey the locality postulate.
The Lagrange density is given by: L(Φ, ∂ν Φ) = − 21 (∂ν Φ∂ν Φ + m2 Φ2 ). The energy operator is given
by: Z X
H = H(x)d3 x = ℏωk a† (⃗k )a(⃗k )
⃗
k
Chapter 15: Quantum field theory & Particle physics 87
This is a continuity equation ⇒ conservation law. Which quantity is conserved depends on the
symmetry. The above Lagrange density is invariant for a change in phase Φ → Φeiθ : a global gauge
transformation. The conserved quantity is the current density Jµ (x) = −ie(Φ∂µ Φ∗ −Φ∗ ∂µ Φ). Because
this quantity is 0 for real fields a complex field is needed to describe charged particles. When this
field is quantized the field operators are given by
1 X 1 ⃗ ikx 1 X 1 † ⃗ ikx
Φ(x) = √ √ a(k )e + b† (⃗k )e−ikx , Φ† (x) = √ √ a (k )e + b(⃗k )e−ikx
V ⃗ 2ωk V ⃗ 2ωk
k k
From this follows that a† a := N+ (⃗k ) is an occupation number operator for particles with a positive
charge and b† b := N− (⃗k ) is an occupation number operator for particles with a negative charge.
p )e−i(⃗p·⃗x±Et)
ψ(x) = ur± (⃗
Here, r is an indication for the direction of the spin, and ± is the sign of the energy. With the
notation v r (⃗
p ) = ur− (−⃗
p ) and ur (⃗
p ) = ur+ (⃗
p ) one can write for the dot products of these spinors:
′ E ′ E ′
ur+ (⃗
p )ur+ (⃗
p) = δrr′ , ur− (⃗
p )ur− (⃗
p) = δrr′ , ur+ (⃗
p )ur− (⃗
p) = 0
M M
88 Physics Formulary by ir. J.C.A. Wevers
Because of the factor E/M this is not relativistic invariant. A Lorentz-invariant dot product is defined
by ab := a† γ4 b, where a := a† γ4 is a row spinor. From this follows:
′ ′ ′
p )ur (⃗
ur (⃗ p )v r (⃗
p ) = δrr′ , v r (⃗ p )v r (⃗
p ) = −δrr′ , ur (⃗ p) = 0
Combinations of the type aa give a 4 × 4 matrix:
2 2
X −iγλ pλ + M X −iγλ pλ − M
ur (⃗
p )ur (⃗
p) = , v r (⃗
p )v r (⃗
p) =
r=1
2M r=1
2M
The Lagrange density which results in the Dirac equation and having the correct energy normalization
is:
∂
L(x) = −ψ(x) γµ + M ψ(x)
∂xµ
and the current density is Jµ (x) = −ieψγµ ψ.
and r
MX 1 X †
p )e−ipx + dr (⃗ p )eipx
ψ(x) = √ p )ur (⃗
cr (⃗ p )v r (⃗
V E r
p
⃗
Here, c† and c are the creation respectively annihilation operators for an electron and d† and d the
creation respectively annihilation operators for a positron. The energy operator is given by
X 2
X
c†r (⃗ p )d†r (⃗
H= Ep⃗ p ) − dr (⃗
p )cr (⃗ p)
p
⃗ r=1
To prevent that the energy of positrons is negative the operators must obey anti commutation rules
in stead of commutation rules:
p ), c†r′ (⃗
[cr (⃗ p ), d†r′ (⃗
p )]+ = [dr (⃗ p )]+ = δrr′ δpp′ , all other anti commutators are 0.
The field operators obey
[ψα (x), ψβ (x′ )] = 0 , [ψα (x), ψβ (x′ )] = 0 , [ψα (x), ψβ (x′ )]+ = −iSαβ (x − x′ )
∂
with S(x) = γλ − M ∆(x)
∂xλ
The anti commutation rules give besides the positive-definite energy also the Pauli exclusion principle
and the Fermi-Dirac statistics: because c†r (⃗ p )c†r (⃗
p ) = −c†r (⃗
p )c†r (⃗
p ) holds: {c†r (p)}2 = 0. It appears
to be impossible to create two electrons with the same momentum and spin. This is the exclusion
principle. Another way to see this is the fact that {Nr+ (⃗ p )}2 = Nr+ (⃗ p ): the occupation operators
have only eigenvalues 0 and 1.
To avoid infinite vacuum contributions to the energy and charge the normal product is introduced.
The expression for the current density now becomes Jµ = −ieN (ψγµ ψ). This product is obtained by:
Expand all fields into creation and annihilation operators,
Keep all terms which have no annihilation operators, or in which they are on the right of the
creation operators,
In all other terms interchange the factors so that the annihilation operators go to the right.
By an interchange of two fermion operators add a minus sign, by interchange of two boson
operators not. Assume hereby that all commutators are zero.
Chapter 15: Quantum field theory & Particle physics 89
The operators obey [am (⃗k ), a†m′ (⃗k )] = δmm′ δkk′ . All other commutators are 0. m gives the polariza-
tion direction of the photon: m = 1, 2 gives transversal polarized, m = 3 longitudinal polarized and
m = 4 timelike polarized photons. Further holds:
In spite of the fact that A4 = iV is imaginary in the classical case, A4 is still defined to be hermi-
tian because otherwise the sign of the energy becomes incorrect. By changing the definition of the
inner product in configuration space the expectation values for A1,2,3 (x) ∈ IR and for A4 (x) become
imaginary.
If the potentials satisfy the Lorentz gauge condition ∂µ Aµ = 0 the E and B operators derived
from these potentials will satisfy the Maxwell equations. However, this gives problems with the
commutation rules. It is now demanded that only those states are permitted for which holds
∂A+µ
|Φ⟩ = 0
∂xµ
∂Aµ
This results in: = 0.
∂xµ
From this follows that (a3 (⃗k ) − a4 (⃗k ))|Φ⟩ = 0. With a local gauge transformation one obtains
N3 (⃗k ) = 0 and N4 (⃗k ) = 0. However, this only applies to free EM-fields: in intermediary states in
interactions there can exist longitudinal and timelike photons. These photons are also responsible for
the stationary Coulomb potential.
Here, the T -operator means a time-ordered product: the terms in such a product must be ordered in
increasing time order from the right to the left so that the earliest terms act first. The S-matrix is
then given by: Sij = ⟨Φi |S|Φj ⟩ = ⟨Φi |Φ(∞)⟩.
The interaction Hamilton density for the interaction between the electromagnetic and the electron-
positron field is: Hint (x) = −Jµ (x)Aµ (x) = ieN (ψγµ ψAµ )
When this is expanded as: Hint = ieN (ψ + + ψ − )γµ (ψ + + ψ − )(A+ −
µ + Aµ )
90 Physics Formulary by ir. J.C.A. Wevers
eight terms appear. Each term corresponds with a possible process. The term ieψ + γµ ψ + A− µ acting
− + +
on |Φ⟩ gives transitions where Aµ creates a photon, ψ annihilates an electron and ψ annihilates a
positron. Only terms with the correct number of particles in the initial and final state contribute to a
matrix element ⟨Φi |S|Φj ⟩. Further the factors in Hint can create and thereafter annihilate particles:
the virtual particles.
The expressions for S (n) contain time-ordered products of normal products. This can be written
as a sum of normal products. The appearing operators describe the minimal changes necessary to
change the initial state into the final state. The effects of the virtual particles are described by the
(anti)commutator functions. Some time-ordened products are:
T {Φ(x)Φ(y)} = N {Φ(x)Φ(y)} + 21 ∆F (x − y)
n o n o
F
T ψα (x)ψβ (y) = N ψα (x)ψβ (y) − 21 Sαβ (x − y)
F
T {Aµ (x)Aν (y)} = N {Aµ (x)Aν (y)} + 21 δµν Dµν (x − y)
The term 21 ∆F (x − y) is called the contraction of Φ(x) and Φ(y), and is the expectation value of the
time-ordered product in the vacuum state. Wick’s theorem gives an expression for the time-ordened
product of an arbitrary number of field operators. The graphical representation of these processes
are called Feynman diagrams. In the x-representation each diagram describes a number of processes.
The contraction functions can also be written as:
−2i eikx −2i iγµ pµ − M 4
Z Z
∆F (x) = lim d 4
k and S F
(x) = lim eipx 2 d p
ϵ→0 (2π)4 k 2 + m2 − iϵ ϵ→0 (2π)4 p + M 2 − iϵ
In the expressions for S (2) this gives rise to terms δ(p + k − p′ − k ′ ). This means that energy and
momentum is conserved. However, virtual particles do not obey the relation between energy and
momentum.
Geometrical conservation laws are invariant under Lorentz transformations and the CPT-operation.
These are:
1. Mass/energy because the laws of nature are invariant for translations in time.
2. Momentum because the laws of nature are invariant for translations in space.
3. Angular momentum because the laws of nature are invariant for rotations.
Dynamical conservation laws are invariant under the CPT-operation. These are:
1. Electrical charge because the Maxwell equations are invariant under gauge transformations.
2. Colour charge is conserved.
3. Isospin because QCD is invariant for rotations in T-space.
4. Baryon number and lepton number are conserved but not under a possible SU(5) symmetry of
the laws of nature.
5. Quarks type is only conserved under the colour interaction.
6. Parity is conserved except for weak interactions.
Quarks exist in three colours but because they are confined these colours cannot be seen directly. The
color force does not decrease with distance. The potential energy will become high enough to create
a quark-antiquark pair when it is tried to disjoin an (anti)quark from a hadron. This will result in
two hadrons and not in free quarks.
3. There is no connection with the neutrino: the decay of the Λ particle through: Λ → p+ + π −
and Λ → n0 + π 0 has also these properties.
The CP-symmetry was found to be violated by the decay of neutral Kaons. These are the lowest
possible states with a s-quark so they can decay only weakly. The following holds: C|K0 ⟩ = η|K0 ⟩
where η is a phase factor. Further holds P|K0 ⟩ = −|K0 ⟩ because K0 and K0 have an intrinsic parity
of −1. From this follows that K0 and K0 are not eigenvalues of CP: CP|K0 ⟩ = |K0 ⟩. The linear
combinations √ √
|K01 ⟩ := 21 2(|K0 ⟩ + |K0 ⟩) and |K02 ⟩ := 21 2(|K0 ⟩ − |K0 ⟩)
Chapter 15: Quantum field theory & Particle physics 93
are eigenstates of CP: CP|K01 ⟩ = +|K01 ⟩ and CP|K02 ⟩ = −|K02 ⟩. A base of K01 and K02 is practical while
describing weak interactions. For colour interactions a base of K0 and K0 is practical because then
the number u−number u is constant. The expansion postulate must be used for weak decays:
to find a final state with CP= −1 is 21 | K02 |K0 |2 , the probability of CP=+1 decay
The probability
is 21 | K01 |K0 |2 .
The relation between the mass eigenvalues of the quarks (unaccented) and the fields arising in the
weak currents (accented) is (u′ , c′ , t′ ) = (u, c, t), and:
′
d 1 0 0 1 0 0 cos θ1 sin θ1 0
s′ = 0 cos θ2 sin θ2 0 1 0 − sin θ1 cos θ1 0
′
b 0 − sin θ2 cos θ2 0 0 eiδ 0 0 1
1 0 0 d
0 cos θ3 sin θ3 s
0 − sin θ3 cos θ3 b
Here the Lk are the generators of the gauge group (the “charges”) and the Akµ are the gauge fields.
g is the matching coupling constant. The Lagrange density for a scalar field becomes:
L = − 12 (Dµ Φ∗ Dµ Φ + M 2 Φ∗ Φ) − 14 Fµν
a
Faµν
a
and the field tensors are given by: Fµν = ∂µ Aaν − ∂ν Aaµ + gcalm Alµ Am
ν .
0 1 0 0
the left- and right- handed solutions of the Dirac equation for neutrino’s are given by:
ψL = 21 (1 + γ5 )ψ and ψR = 12 (1 − γ5 )ψ
It appears that neutrino’s are always left-handed while antineutrino’s are always right-handed. The
hypercharge Y , for quarks given by Y = B + S + C + B∗ + T′ , is defined by:
Q = 21 Y + T3
so [Y, Tk ] = 0. The group U(1)Y ⊗SU(2)T is taken as symmetry group for the electroweak interaction
because the generators of this group commute. The multiplets are classified as follows:
94 Physics Formulary by ir. J.C.A. Wevers
e−
R νeL e−
L uL d′L uR dR
1 1
T 0 2 2 0 0
1 1 1 1
T3 0 2 − 2 2 − 2 0 0
1 4
Y −2 −1 3 3 − 32
Now, 1 field Bµ (x) is connected with gauge group U(1) and 3 gauge fields A ⃗ µ (x) are connected with
SU(2). The total Lagrange density (minus the fieldterms) for the electron-fermion field now becomes:
′
g⃗ 1 g ψνe,L
L0,EW = −(ψνe,L , ψeL )γ µ ∂µ − i A µ · ( 1
2 ⃗
σ ) − i
2 ℏ µ B · (−1) −
ℏ ψeL
g′
ψeR γ µ ∂µ − 12 i (−2)Bµ ψeR
ℏ
Φ+ iω + √
0√
Φ= = and ⟨0|Φ|0⟩ =
Φ0 (v + ϕ − iz)/ 2 v/ 2
Because this expectation value ̸= 0 the SU(2) symmetry is broken but the U(1) symmetry is not.
When the gauge fields in the resulting Lagrange density are separated one obtains:
√ √
Wµ− = 12 2(A1µ + iA2µ ) , Wµ+ = 21 2(A1µ − iA2µ )
gA3µ − g ′ Bµ
Zµ = p ≡ A3µ cos(θW ) − Bµ sin(θW )
g 2 + g ′2
g ′ A3µ + gBµ
Aµ = p ≡ A3µ sin(θW ) + Bµ cos(θW )
g 2 + g ′2
where θW is called the Weinberg angle. For this angle holds: sin2 (θW ) = 0.255 ± 0.010. Relations
for the masses of the field quanta can be obtained from the remaining terms: MW = 12 vg and
gg ′
= g ′ cos(θW ) = g sin(θW )
p
MZ = 12 v g 2 + g ′2 , and for the elementary charge holds: e = p
g 2 + g ′2
Experimentally it is found that MW = 80.022 ± 0.26 GeV/c2 and MZ = 91.187 ± 0.007 GeV/c2 .
According to the weak theory this should be: MW = 83.0 ± 0.24 GeV/c2 and MZ = 93.8 ± 2.0
GeV/c2 .
15.13.3 Quantumchromodynamics
Coloured particles interact because the Lagrange density is invariant for the transformations of the
group SU(3) of the colour interaction. A distinction can be made between two types of particles:
2. “Coloured” particles: the generators T⃗ are 8 3 × 3 matrices. There exist three colours and three
anticolours.
The Lagrange density for coloured particles is given by
X X
LQCD = i Ψk γ µ Dµ Ψk + a
Ψk Mkl Ψl − 41 Fµν Faµν
k k,l
The gluons remain massless because this Lagrange density does not contain spinless particles. Because
left- and right- handed quarks now belong to the same multiplet a mass term can be introduced. This
term can be brought in the form Mkl = mk δkl .
Astrophysics
vr = V cos(θ) , vt = V sin(θ) = ωR
-5
where θ is the angle between the star and the point of convergence and -4
R̂ the distance in pc. This results, with vt = vr tan(θ), in: Type 1
⟨M ⟩ -3
vr tan(θ) 1′′ -2
R= ⇒ R̂ = -1 Type 2
ω p
0
where p is the parallax in arc seconds. The parallax is then given by RR-Lyrae
1
0,1 0,3 1 3 10 30 100
4.74µ
p= P (days) →
vr tan(θ)
with µ de proper motion of the star in ′′ /yr. A method to determine the distance of objects which
are somewhat further away, like galaxies and star clusters, uses the period-Brightness relation for
Cepheids. This relation is shown in the above figure for different types of stars.
where πFν is the monochromatic radiation flux. At the position of an observer this is πfν , with
fν = (R/r)2 Fν if absorption is ignored. If Aν is the fraction of the flux which reaches Earth’s surface,
the transmission factor is given by Rν and the surface of the detector is given by πa2 , then the
apparent brightness b is given by:
Z∞
2
b = πa fν Aν Rν dν
0
b1 1
= (100) 5 (m2 −m1 ) = (2.512)m2 −m1
b2
Chapter 16: Astrophysics 97
because the human eye perceives lightintensities logaritmical. From this follows that m2 − m1 =
2.5 ·10 log(b1 /b2 ), or: m = −2.5 ·10 log(b) + C. The apparent brightness of a star if this star would be
at a distance of 10 pc is called the absolute brightness B: B/b = (r̂/10)2 . The absolute magnitude is
then given by M = −2.5 ·10 log(B) + C, or: M = 5 + m − 5 ·10 log(r̂). When an interstellar absorption
of 10−4 /pc is taken into account one finds:
If a detector detects all radiation emitted by a source one would measure the absolute bolometric
magnitude. If the bolometric correction BC is given by
R
10 Energy flux received 10 fν dν
BC = 2.5 · log = 2.5 · log R
Energy flux detected fν Aν Rν dν
holds: Mb = MV − BC where MV is the visual magnitude. Further holds
10 L
Mb = −2.5 · log + 4.72
L⊙
where X is the mass fraction of H, Y the mass fraction of He and Z the mass fraction of the other
elements. Further holds:
Convection will occur when the star meets the Schwartzschild criterium:
dT dT
<
dr conv dr rad
Otherwise the energy transfer takes place by radiation. For stars in quasi-hydrostatic equilibrium
hold the approximations r = 21 R, M (r) = 12 M , dM/dr = M/R, κ ∼ ϱ and ε ∼ ϱT µ (this last
assumption is only valid for stars on the main sequence). For pp-chains holds µ ≈ 5 and for the CNO
chains holds µ = 12 tot 18. It can be derived that L ∼ M 3 : the mass-brightness relation. Further
holds: L ∼ R4 ∼ Teff
8
. This results in the equation for the main sequence in the Hertzsprung-Russel
diagram:
10
log(L) = 8 ·10 log(Teff ) + constant
The ∇-operator
In cartesian coordinates (x, y, z) holds:
⃗ ∂ 1 ∂ 1 ∂
∇ = ⃗er + ⃗eθ + ⃗eφ
∂r r ∂θ r sin θ ∂φ
∂f 1 ∂f 1 ∂f
gradf = ⃗er + ⃗eθ + ⃗eφ
∂r r ∂θ r sin θ ∂φ
∂ar 2ar 1 ∂aθ aθ 1 ∂aφ
div ⃗a = + + + +
∂r r r ∂θ r tan θ r sin θ ∂φ
1 ∂aφ aθ 1 ∂aθ 1 ∂ar ∂aφ aφ
rot ⃗a = + − ⃗er + − − ⃗eθ +
r ∂θ r tan θ r sin θ ∂φ r sin θ ∂φ ∂r r
∂aθ aθ 1 ∂ar
+ − ⃗eφ
∂r r r ∂θ
∂2f 2 ∂f 1 ∂2f 1 ∂f 1 ∂2f
∇2 f = 2
+ + 2 2 + 2 + 2 2
∂r r ∂r r ∂θ r tan θ ∂θ r sin θ ∂φ2
General orthonormal curvelinear coordinates (u, v, w) can be obtained from cartesian coordinates by
the transformation ⃗x = ⃗x(u, v, w). The unit vectors are then given by:
1 ∂⃗x 1 ∂⃗x 1 ∂⃗x
⃗eu = , ⃗ev = , ⃗ew =
h1 ∂u h2 ∂v h3 ∂w
where the factors hi set the norm to 1. Then holds:
1 ∂f 1 ∂f 1 ∂f
gradf = ⃗eu + ⃗ev + ⃗ew
h1 ∂u h2 ∂v h3 ∂w
1 ∂ ∂ ∂
div ⃗a = (h2 h3 au ) + (h3 h1 av ) + (h1 h2 aw )
h1 h2 h3 ∂u ∂v ∂w
1 ∂(h3 aw ) ∂(h2 av ) 1 ∂(h1 au ) ∂(h3 aw )
rot ⃗a = − ⃗eu + − ⃗ev +
h2 h3 ∂v ∂w h3 h1 ∂w ∂u
1 ∂(h2 av ) ∂(h1 au )
− ⃗ew
h1 h2 ∂u ∂v
1 ∂ h2 h3 ∂f ∂ h3 h1 ∂f ∂ h1 h2 ∂f
∇2 f = + +
h1 h2 h3 ∂u h1 ∂u ∂v h2 ∂v ∂w h3 ∂w
100 The SI units
The SI units
Basic units Derived units with special names
Quantity Unit Sym. Quantity Unit Sym. Derivation
Length metre m Frequency hertz Hz s−1
Mass kilogram kg Force newton N kg · m · s−2
Time second s Pressure pascal Pa N · m−2
Therm. temp. kelvin K Energy joule J N·m
Electr. current ampere A Power watt W J · s−1
Luminous intens. candela cd Charge coulomb C A·s
Amount of subst. mol mol El. Potential volt V W · A−1
El. Capacitance farad F C · V−1
Extra units El. Resistance ohm Ω V · A−1
Plane angle radian rad El. Conductance siemens S A · V−1
solid angle sterradian sr Mag. flux weber Wb V·s
Mag. flux density tesla T Wb · m−2
Inductance henry H Wb · A−1
Luminous flux lumen lm cd · sr
Illuminance lux lx lm · m−2
Activity bequerel Bq s−1
Absorbed dose gray Gy J · kg−1
Dose equivalent sievert Sv J · kg−1
Prefixes
quetta Q 1030 deci d 10−1
ronna R 1027 centi c 10−2
yotta Y 1024 milli m 10−3
zetta Z 1021 micro µ 10−6
exa E 1018 nano n 10−9
peta P 1015 pico p 10−12
tera T 1012 femto f 10−15
giga G 109 atto a 10−18
mega M 106 zepto z 10−21
kilo k 103 yocto y 10−24
hecto h 102 ronto r 10−27
deca da 10 quecto q 10−30