Dierential geometry for Physicists

Similar documents
Solutions of M3-4A16 Assessed Problems # 3 [#1] Exercises in exterior calculus operations

BACKGROUND IN SYMPLECTIC GEOMETRY

Survey on exterior algebra and differential forms

Math 225B: Differential Geometry, Final

Gravitation: Tensor Calculus

GEOMETRIC QUANTIZATION

EXERCISES IN POISSON GEOMETRY

Symplectic geometry of deformation spaces

Compatible Discretization Schemes

CALCULUS ON MANIFOLDS. 1. Riemannian manifolds Recall that for any smooth manifold M, dim M = n, the union T M =

MATH 4030 Differential Geometry Lecture Notes Part 4 last revised on December 4, Elementary tensor calculus

A DIFFERENTIAL GEOMETRIC APPROACH TO FLUID MECHANICS

Modern Geometric Structures and Fields

[#1] Exercises in exterior calculus operations

The Geometry of Euler s equation. Introduction

Seminar Geometrical aspects of theoretical mechanics

Deformations of coisotropic submanifolds in symplectic geometry

Exercises in Geometry II University of Bonn, Summer semester 2015 Professor: Prof. Christian Blohmann Assistant: Saskia Voss Sheet 1

Caltech Ph106 Fall 2001

Symplectic and Poisson Manifolds

Symplectic Geometry versus Riemannian Geometry

Elements of differential geometry

Lecture I: Constrained Hamiltonian systems

Chap. 1. Some Differential Geometric Tools

Math 868 Final Exam. Part 1. Complete 5 of the following 7 sentences to make a precise definition (5 points each). Y (φ t ) Y lim

Hamiltonian flows, cotangent lifts, and momentum maps

Gauge Fixing and Constrained Dynamics in Numerical Relativity

An Introduction to Symplectic Geometry

Solutions to the Hamilton-Jacobi equation as Lagrangian submanifolds

Curves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations,

fy (X(g)) Y (f)x(g) gy (X(f)) Y (g)x(f)) = fx(y (g)) + gx(y (f)) fy (X(g)) gy (X(f))

Contact and Symplectic Geometry of Monge-Ampère Equations: Introduction and Examples

Section 2. Basic formulas and identities in Riemannian geometry

PHY411 Lecture notes Part 2

Reminder on basic differential geometry

LECTURE 1: LINEAR SYMPLECTIC GEOMETRY

Eva Miranda. UPC-Barcelona and BGSMath. XXV International Fall Workshop on Geometry and Physics Madrid

Let V, W be two finite-dimensional vector spaces over R. We are going to define a new vector space V W with two properties:

NOTES ON DIFFERENTIAL FORMS. PART 3: TENSORS

CS 468. Differential Geometry for Computer Science. Lecture 13 Tensors and Exterior Calculus

M4P52 Manifolds, 2016 Problem Sheet 1

THE LAGRANGIAN AND HAMILTONIAN MECHANICAL SYSTEMS

[#1] R 3 bracket for the spherical pendulum

Differential Forms, Integration on Manifolds, and Stokes Theorem

GEOMETRICAL TOOLS FOR PDES ON A SURFACE WITH ROTATING SHALLOW WATER EQUATIONS ON A SPHERE

NOTES ON DIFFERENTIAL FORMS. PART 1: FORMS ON R n

LIE ALGEBROIDS AND POISSON GEOMETRY, OLIVETTI SEMINAR NOTES

Hamiltonian Systems of Negative Curvature are Hyperbolic

1. Geometry of the unit tangent bundle

Non-associative Deformations of Geometry in Double Field Theory

DIFFERENTIAL FORMS AND COHOMOLOGY

CS 468 (Spring 2013) Discrete Differential Geometry

Lorentzian elasticity arxiv:

William P. Thurston. The Geometry and Topology of Three-Manifolds

THE NEWLANDER-NIRENBERG THEOREM. GL. The frame bundle F GL is given by x M Fx

Homogeneous para-kähler Einstein manifolds. Dmitri V. Alekseevsky

HOMOGENEOUS AND INHOMOGENEOUS MAXWELL S EQUATIONS IN TERMS OF HODGE STAR OPERATOR

INTRODUCTION TO ALGEBRAIC GEOMETRY

Poincaré (non-holonomic Lagrange) Equations

Math 550 / David Dumas / Fall Problems

SOME HIGHLIGHTS OF SYMPLECTIC GEOMETRY, SUMMER SCHOOL UTRECHT, GEOMETRY, AUGUST 16, 2017

Relativistic Mechanics

The BRST complex for a group action

Towards Discrete Exterior Calculus and Discrete Mechanics for Numerical Relativity

COTANGENT MODELS FOR INTEGRABLE SYSTEMS

Dynamical systems on Leibniz algebroids

Introduction to finite element exterior calculus

Twisted Poisson manifolds and their almost symplectically complete isotropic realizations

LECTURE 5: COMPLEX AND KÄHLER MANIFOLDS

The geometry of Euler s equations. lecture 2

2.20 Fall 2018 Math Review

Poisson geometry of b-manifolds. Eva Miranda

Lie algebra cohomology

Patrick Iglesias-Zemmour

Vortex Equations on Riemannian Surfaces

Hamilton-Jacobi theory

Hamiltonian Dynamics

arxiv:solv-int/ v1 24 Mar 1997

HONGYU HE. p = mẋ; q = x. be the Hamiltonian. It represents the energy function. Then H q = kq,

Smooth Dynamics 2. Problem Set Nr. 1. Instructor: Submitted by: Prof. Wilkinson Clark Butler. University of Chicago Winter 2013

PREQUANTIZATION OF SYMPLECTIC SUPERMANIFOLDS

Complex manifolds, Kahler metrics, differential and harmonic forms

LECTURE 28: VECTOR BUNDLES AND FIBER BUNDLES

Lagrangian submanifolds and generating functions

Lecture: Lorentz Invariant Dynamics

Lecture Notes a posteriori for Math 201

Applications of Differential Geometry to Physics

Vector fields and differential forms. William G. Faris

Many of the exercises are taken from the books referred at the end of the document.

Invariant Lagrangian Systems on Lie Groups

SYMPLECTIC MANIFOLDS, GEOMETRIC QUANTIZATION, AND UNITARY REPRESENTATIONS OF LIE GROUPS. 1. Introduction

Synchro-Betatron Motion in Circular Accelerators

Invariant Nonholonomic Riemannian Structures on Three-Dimensional Lie Groups

Differential Geometry MTG 6257 Spring 2018 Problem Set 4 Due-date: Wednesday, 4/25/18

5 Constructions of connections

Dirac structures. Henrique Bursztyn, IMPA. Geometry, mechanics and dynamics: the legacy of J. Marsden Fields Institute, July 2012

COMPUTING THE POISSON COHOMOLOGY OF A B-POISSON MANIFOLD

5.3 Surface Theory with Differential Forms

LECTURE 9: MOVING FRAMES IN THE NONHOMOGENOUS CASE: FRAME BUNDLES. 1. Introduction

The prototypes of smooth manifolds

Transcription:

Dierential geometry for Physicists (What we discussed in the course of lectures) Marián Fecko Comenius University, Bratislava Syllabus of lectures delivered at University of Regensburg in June and July 2012 Contents 1 Smooth manifolds 2 1.1 Manifolds................................ 2 1.2 Vectors and vector elds on manifolds................ 2 1.3 Tensors in linear algebra........................ 3 1.4 Tensors and tensor elds on manifolds................ 4 1.5 Mapping of tensors induced by mapping of manifolds........ 4 1.6 Lie derivative, isometries, Killing vectors............... 5 2 Dierential forms 6 2.1 Volumes of parallelepipeds and forms in linear algebra....... 6 2.2 Dierential calculus of forms on manifolds.............. 6 2.3 Closed and exact forms......................... 7 2.4 Forms on manifolds - integral calculus................ 7 3 Hamiltonian mechanics and symplectic manifolds 8 3.1 How Poisson tensor emerges...................... 8 3.2 How symplectic form emerges..................... 8 3.3 Hamiltonian elds and Poisson brackets............... 9 3.4 Symmetries and conserved quantities (no action integral)...... 9 3.5 Canonical transformations....................... 10 3.6 Poincaré integral invariants...................... 10 3.7 Volume form and (the corresponding) divergence of a vector eld. 10 3.8 Algebra of observables of classical mechanics............. 11 3.9 Cotangent bundle T M as a phase space............... 11 3.10 Time-dependent Hamiltonians - what changes are needed..... 11 e-mail: fecko@fmph.uniba.sk 1

3.11 Action integral, Hamilton's principle................. 12 3.12 Symmetries and conserved quantities (based on action integral).. 12 3.13 Poincaré-Cartan integral invariants.................. 12 4 Field theory in the language of forms 13 4.1 The Hodge star operator........................ 13 4.2 Forms in E 3 and how vector analysis drops out........... 14 4.3 Forms in Minkowski space E 1,3.................... 14 4.4 Maxwell equations in terms of forms................. 15 4.5 Some immediate consequences..................... 15 4.6 The action formulation of Maxwell equations............ 16 1 Smooth manifolds 1.1 Manifolds Local (as opposed to global) coordinates x i, change of coordinates x i (x) on overlaps of domains, charts, atlas, conguration spaces in mechanics as manifolds (double plane and spherical pendulum), (smooth) mapping of manifolds, injective, surjective and bijective mappings, coordinate presentation y a (x) of a mapping of manifolds, a curve γ on a manifold M (as a mapping γ : R M of manifolds), its coordinate presentation x i (t), a function ψ on a manifold M (as a mapping ψ : M R of manifolds), its coordinate presentation ψ(x i ) 1.2 Vectors and vector elds on manifolds Curves tangent at a point, tangency as an equivalence relation γ 1 γ 2, equivalence class [γ] (of curves; γ is a representative), linear combination of such equivalence classes, directional derivative, algebra F(M) of functions on M, four (equivalent) denitions of a vector in m M (vector as: equivalence class of curves, derivation of algebra of functions plus Leibniz rule, expression a i i x, set of numbers a i ), the tangent space T m M at m M - the (linear) space of vectors in m M (in any of the four versions), vector eld V as a rst order dierential operator V i (x) i, 2

transformation law for its components V i (x) (from V i i = V i i ), integral curve of a vector eld, equations ẋ i = V i for nding integral curves, ow Φ t : M M of a vector eld, m γ(0) γ(t) =: Φ t (m) 1.3 Tensors in linear algebra Finite dimensional linear space L, the dual space L, a basis e a L, the dual basis e a L, dened by e a, e b = δ a b, tensor of type (p, q) in L (as a multi-linear map), identication of well known examples, components t a...b c...d, space T q p (L) of tensors of type (p, q) in L, various roles played by the same tensor, tensor product and its properties, the basis e a e b e c e d in T p q (L) induced by a basis e a in L, the unit tensor metric tensor the inverse metric tensor (co-metric) lowering of indices raising of indices contraction operation as a mapping e a e a δ a b, g = g ab e a e b, g ab e a e b, : L L (v a v a := g ab v b ), : L L (α a α a := g ab α b ), C : T p q (L) T p 1 q 1 (L) 3

1.4 Tensors and tensor elds on manifolds A tensor in m on M: take T m M as L from linear algebra, the space T p q (M) of tensor elds of type (p, q) on M, the gradient dψ of a function ψ as a covector eld on M, dx i as the dual basis to i, the (coordinate) basis dx i dx j k l in T p q (M) (on a coordinate patch) induced by a (coordinate) basis i for vector elds, metric tensor g = g ij dx i dx j on a manifold, Riemannian manifold, transformation of components under change of coordinates from and where g = g ij dx i dx j = g i j dxi dx j dx i = J i k dxk, J i k (x) xi / x j is the (always square) Jacobi matrix of the transformation of coordinates, the Euclidean metric tensor in the plane g ij = δ ij, i.e. g = dx dx + dy dy = dr dr + r 2 dφ dφ, the length functional γ t2 t 1 γ dt t2 t 1 g( γ, γ)dt t2 t 1 g ij ẋ i ẋ j dt, the gradient as a vector eld i.e. ψ := dψ ( ψ) i = g ij j ψ 1.5 Mapping of tensors induced by mapping of manifolds Push-forward v f v of a vector in x M to a vector in f(x) N (induced by f : M N), explicitly v i i (J a i v i ) a, 4

or, on the coordinate basis, where f : i f i = J a i a, J a i (x) y a (x)/ x j is the Jacobi matrix (in general non-square) of the mapping f, pull-back of a function ψ on N, the result being a function f ψ = ψ f on M (in coordinates ψ(y) ψ(y(x))), pull-back α f α of a covector in f(x) N to a covector in x M (induced by f : M N), given as f α, v := α, f v, on the coordinate basis where dy a f dy a = dy a (x) = J a i dx i, J a i (x) y a (x)/ x j is the Jacobi matrix of the mapping f (the same as for push-forward of vectors), an alternative (equivalent) formula is f (α a (y)dy a ) = α a (y(x))dy a (x) = α a (y(x))j a i (x)dx i, pull-back of a metric tensor - the induced metric tensor, its example in Lagrangian mechanics, "curved" and "at" torus, general properties of pull-back - behavior on tensor product and linear combination, commutation with taking of the gradient (f d = df on functions) 1.6 Lie derivative, isometries, Killing vectors Pull-back of a metric tensor w.r.t. the innitesimal ow generated by a vector eld V, comparison of a tensor in x with that pulled-back from Φ ϵ (x), denition of Lie derivative, component computation of Lie derivative, general properties of Lie derivative, isometry = such f : M M which preserves length of any curve, one-parameter group (ow) of isometries, Killing equations L V g = 0, 5

commutator of vector elds [V, W ], Lie algebra of solutions of Killing equations, explicit solution for the common plane, explicit solution for (pseudo)-euclidean spaces, Lorentz algebra (rotations and boosts = hyperbolic rotations = pseudo-rotations), Poincaré algebra (rotations, boosts and translations), strain tensor (small deformations of continuous media) 2 Dierential forms 2.1 Volumes of parallelepipeds and forms in linear algebra Parallelepiped associated with several vectors, degenerate parallelepiped, a p-form in L, anti-symmetrization operation the wedge product t i...j t [i...j], α β (bilinear, associative, graded commutative), basis forms expression e a e b, α = (1/p!)α a...b e a e b of a general p-form, how it helps in practical wedge multiplication, interior product its nice properties i v α, 2.2 Dierential calculus of forms on manifolds Comma operation does not work as tensor operation, it works as an operation on forms (when followed by square brackets), a good idea t i...j t [i...j,k], exterior derivative it is nilpotent d : Ω p Ω p+1, dd = 0 6

and it obeys graded Leibniz rule d(α β) = dα β + ( 1) p α dβ for α = a p-form, Cartan's (magic) formula L V = di V + i V d, Cartan's formulas for dα(u, V,... ), Lie derivative commutes with the exterior one [L V, d] = 0, formula [L V, i W ] = i [V,W ] 2.3 Closed and exact forms Closed form - such α that exact form - such α that dα = 0, α = dβ for some β (if it exists, β is called potential), the freedom β β + d(... ), exact closed always holds because of dd = 0, closed exact not always holds, but it does hold locally, e.g. in a coordinate patch (more detailed information about the relation between closed and exact forms on M is given by "derham cohomology" theory) 2.4 Forms on manifolds - integral calculus An intuitive picture why exactly dierential forms are integrated, Stokes formula α = dα, D some special cases (Newton-Leibniz formula, area under the graph of a function, Green's formula, integration by parts), the formula α = f α f(d) D D 7

3 Hamiltonian mechanics and symplectic manifolds 3.1 How Poisson tensor emerges Hamiltonian equations ẋ a = H/ p a ṗ a = H/ x a as linear relations between dots of (all) coordinates and (all) components od dh, most easily visible when (x a, p a ) coordinates are relabeled to z A, A = 1,... 2n, we get ż A = (dh) B P BA with P AB ( ) 0 1 1 0 it denes (skew-symmetric) Poisson tensor P of type ( 2 0), which raises index on the gradient dh. Poisson bracket is then and Hamilton equations read {f, g} = P(df, dg) γ = V H V H = P(dH,. ) 3.2 How symplectic form emerges Poisson tensor P P AB of type ( 2 0) denes a (skew-symmetric) tensor ω ωab of type ( 0 2) as its "inverse" via, so we obtain a 2-form P AC ω CB = δb A, ( ) 0 1 ω AB, 1 0 ω = (1/2)ω AB dz A dz B = dp a dx a which is closed (dω = 0) and non-degenerate (det ω AB 0); any closed and nondegenerate 2-form is called symplectic form; Hamilton equations read γ = V H i VH ω = dh 8

3.3 Hamiltonian elds and Poisson brackets Hamiltonian eld V f generated by a function f is dened either as V f = P(df,. ) (i.e. via raising of index on df with the help of P), or, in terms of the symplectic form, as V f obeying i Vf ω = df (denitions are, for non-degenerate P, equivalent); Poisson brackets may be expressed in several ways, too, as {f, g} = P(df, dg) = ω(v f, V g ) = V f g = V g f ; closure of Hamiltonian elds w.r.t. the commutator: [V f, V g ] = V {f,g} ; Jacobi identity for Poisson bracket is equivalent to dω = 0 (i.e. to the fact that ω is closed), invariance of ω w.r.t. (any) Hamiltonian ow L Vf ω = 0, comparison with isometries and Killing vectors 3.4 Symmetries and conserved quantities (no action integral) What is a symmetry of a Hamiltonian triple (M, ω, H), Cartan symmetries, exact Cartan symmetries, i.e. Hamiltonian elds V f, whose generators f obey and the bijection {H, f} = 0 V f f onto conserved quantities f, new solutions γ s (t) Φ f s (γ(t)) from an old solution γ(t) and a symmetry V f Φ f s 9

3.5 Canonical transformations Darboux theorem (canonical "appearance" of a closed 2-form), its manifestation in symplectic case, transformations of coordinates which preserve canonical "appearance" of ω, dp a dq a = dp a dq a dp a (q, p) dq a (q, p) their manifestation on the appearance of Hamilton equations, two ways of explicit work with canonical transformations (generators and generating functions) 3.6 Poincaré integral invariants A form invariant w.r.t. to a vector eld, consequences for their integrals; on a phase space, ω ω ω ω ω ω... are forms invariant w.r.t. any Hamiltonian eld, Poincaré (absolute) integral invariants ω ω ω ω ω ω..., D 2 D 4 D 6 Liouville theorem as a particular case, canonical Liouville volume form Poincaré (relative) integral invariants θ θ ω c 1 c 3 where Ω ω ω, c 5 θ ω ω..., θ = p a dx a dθ = ω c 1 = c 3 = = 0 3.7 Volume form and (the corresponding) divergence of a vector eld Volume form - a non-zero n-form on n-dimensional manifold, coordinate expression divergence dened as Ω = f(x)dx 1 dx 1 dx n, L V Ω =: (div V )Ω, its geometrical meaning (it measures how volumes are changing due to the ow of V ), divergence-less elds (the ow preserves volumes), examples - metric divergence, symplectic divergence (Hamiltonian elds are divergence-less) 10

3.8 Algebra of observables of classical mechanics Pure states (m (M, ω)) and observables (f F(M)) in Hamiltonian mechanics, prediction of the result of measurement of an observable in a state (f(m) R), time development of pure states m Φ t (m), (Schrödinger-like picture), time development of observables f f t := Φ t f, (Heisenberg-like picture), equivalence of the two pictures f(φ t (m)) = (Φ t f)(m), two products in the algebra F(M) (fg and {f, g}), Hamiltonian ows preserve whole structure of the algebra, more general states (probabilistic distributions ρ on M), their time evolution ρ ρ t := Φ tρ, equivalence of the two pictures (Φ tρ)fω = equations of motion M M ρ(φ t f)ω, t f t = {H, f t } t ρ t = {H, ρ t } 3.9 Cotangent bundle T M as a phase space Construction of T M from M (in our context, construction of the phase space T M associated with a conguration space M), canonical coordinates (x a, p a ), canonical (exact) symplectic form on T M, ber bundle, base space, total space, ber over x, projection 3.10 Time-dependent Hamiltonians - what changes are needed Extended phase space, Hamilton equations in the form Hamilton equations in the form α a = 0 β a = 0, i γ dσ = 0, where dσ = α a β a = d(p a dx a Hdt) 11

3.11 Action integral, Hamilton's principle How action integral for Hamiltonian mechanics can be constructed: S[γ] = σ, how one can easily see that the action is stationary on solutions of Hamilton equations, the dierence between endpoints xation in Lagrangian and Hamiltonian cases 3.12 Symmetries and conserved quantities (based on action integral) If the innitesimal ow Φ ϵ of a vector eld ξ on extended phase space does not change the action, i.e. S[Φ ϵ γ] = S[γ], we speak of a symmetry of the action; then it turns out that the function i ξ σ is the corresponding conserved quantity; e.g. for we get that it is a symmetry if and that γ ξ = t t H = 0 i ξ σ = H (so, the time translation is a symmetry if t H = 0 and the corresponding conserved quantity is H itself, the energy) 3.13 Poincaré-Cartan integral invariants On extended phase space, integrals of σ σ dσ σ dσ dσ..., over any two closed surfaces (i.e. c = 0) of (corresponding) odd dimension encircling the same tube of trajectories (solutions of Hamilton equations) are equal (relative Poincaré-Cartan integral invariants, integrals σ c 1 σ dσ c 3 σ dσ dσ c 5..., 12

where σ = p a dx a Hdt ; if the surfaces lie in hyper-planes of constant time, we return to relative Poincaré integral invariants); similarly, integrals of dσ dσ dσ dσ dσ dσ... over any two domains of (corresponding) even dimension connected by trajectories are equal, absolute Poincaré-Cartan integral invariants = integrals dσ D 2 dσ dσ D 4 dσ dσ dσ D 6... 4 Field theory in the language of forms 4.1 The Hodge star operator Metric volume form (xation of f(x) in Ω = f(x)dx 1 dx 2 dx n ), as a linear isomorphism from p-forms to (n p)-forms, it holds = ±ˆ1 Important scalar product on forms α, β := From one gets where D α β d(α β) = dα β + ˆηα d β dα, β = α, δβ + α β D δ := 1 d ˆη If the surface integral vanishes, we obtain d + = δ 13

4.2 Forms in E 3 and how vector analysis drops out In E 3, the most general 0,1,2 and 3-forms read f A.dr B.dS = (B.dr) hdv = h, so all of them are parametrized either by scalar (0- and 3-forms) or by vector (1- and 2-forms) elds, therefore the exterior derivative Ω 0 d Ω 1 d Ω 2 d Ω 3 produces eectively three dierential operators, of type scalar vector, vector vector and vector scalar, respectively; they are just well-known grad, rot=curl and div operations from vector analysis; then (the general form-version of the) Stokes theorem gives standard integral identities relating integrals of neighboring dimensions 4.3 Forms in Minkowski space E 1,3 In E 1,3, each form may be decomposed as α = dt ŝ + ˆr where ŝ and ˆr are "spatial", i.e. they do not contain dt; that's why spatial forms may be expressed as forms in E 3, but with t present in components, e.g. a 2-form reads α = dt a.dr + b.ds with a(t, r) b(t, r) ; similarly, operations on forms can be expressed through "spatial" operations; in particular dα d(dt ŝ + ˆr) = dt ( tˆr ˆdŝ) + ˆdˆr and e.g. on 2-forms we get α (dt ŝ + ˆr) = dt (ˆ ˆr) + ˆ ˆηŝ ; d(dt a.dr + b.ds) = dt ( t b curl a).ds + (div b)dv and (dt a.dr + b.ds) = dt b.dr a.ds 14

4.4 Maxwell equations in terms of forms Introducing F = dt E.dr B.dS we get df = dt ( t B curl E).dS + ( div B)dV so that the homogeneous half of Maxwell equations simply reads df = 0 ; now F = (dt E.dr B.dS) = dt ( B).dr E.dS and so d F = dt ( t E + curl B).dS (div E)dV Therefore the inhomogeneous half of Maxwell equations reads d F = J where J := dt ( j.ds) + ρdv is the source 3-form. Equivalently, δf = j where j := J 4.5 Some immediate consequences Existence of 4-potential: df = 0 F = da for some 1-form A. Gauge freedom A A A + dψ (Then F := da = da = F.) Consistency of d F = J needs dj = 0 (= the continuity equation, local conservation of the charge) 15

4.6 The action formulation of Maxwell equations From we get S[A] := 1 da, da A, j 2 S[A + ϵa] := S[A] + ϵ( a, δda + j ) So the variation principle gives the equation When denoting δda = j F = da we get (both halves of) Maxwell equations df = 0 δf = j References [1] M.Fecko: Dierential geometry and Lie groups for physicists, Cambridge University Press 2006 (paperback 2011) [2] M.Fecko: Dierential Geometry in Physics, lecture notes from Regensburg 2007, http://sophia.dtp.fmph.uniba.sk/fecko/regensburg_2007.html 16