Oxford Lectures 2025
From Werner KRAUTH
My 2025 Public Lectures at the University of Oxford (UK), entitled Algorithms and computations in theoretical physics, run from 21 January 2025 through 11 March 2025.
Contents |
Lecture 1: 21 January 2025
We start our parallel exploration of physics and of computing with the concept of sampling, the process of producing examples (“samples”) of a probability distribution. In week 1, we consider “direct” sampling (the examples are obtained directly) and, among the many connections to physics, will come across the Maxwell distribution. In 1859, it marked the beginning of the field of statistical physics.
Here are the notes for the first lecture (21 January 2025).
Lecture 2: 28 January 2025
We start with a Special Topic A, where we supplement the first Lecture with two all-important special topics that explore the meaning of convergence in statistics, and the fundamental usefulness of statistical reasoning. We can discuss them here, in the direct-sampling framework, but they are more generally relevant. The strong law of large numbers, for example, will turn into the famous ergodic theorem for Markov chains.
Here are the notes for Special Topic A (28 January 2025).
Then, in this second lecture proper, we consider Markov-chain sampling, from the adults' game on the Monte Carlo heliport to modern ideas on non-reversibility. We discuss a lot of theory, but also six ten-line pseudo-code algorithms, none of them approximate, and all of them as intricate as they are short.
Here are the notes for Lecture 2.
Lecture 3: 04 February 2025
In this third lecture, we consider Markov-chain sampling in an abstract setting in one dimension. We discuss some theory, but also seven ten-line pseudo-code algorithms, none of them approximate, and all of them as intricate as they are short. At the end, we discuss the foundations of statistical mechanics, as seen in a one- dimensional example.
Here are the notes for Lecture 3.
- Here is the Metropolis algorithm for a single particle in an anharmonic potential
- Here is the factorized Metropolis algorithm for a single particle in an anharmonic potential
- Here is the Zig-zag algorithm for a single particle in an anharmonic potential
Lecture 4: 17 February 2025
In this fourth lecture, we treat classical many-particle systems that interact with hard-sphere potentials, a case of importance for the foundations of statistical mechanics but also for its applications. We move from Newtonian mechanics to Boltzmann mechanics and from classical mechanics to statistical mechanics, in a way that is again entirely example-based. We then treat the case of one-dimensional hard spheres and discover the Asakura--Oosawa interaction, the fifth force in nature. We will end with a discussion of the double nature of pressure: kinematic and thermodynamic.
Here are the notes for Lecture 4.
- Here is the event-driven algorithm for four hard disks in a square box, in a few dozen lines of code.
Lecture 5: 18 February 2025
In this fourth lecture, we treat classical many-particle systems that interact with hard-sphere potentials, a case of importance for the foundations of statistical mechanics but also for its applications. We move from Newtonian mechanics to Boltzmann mechanics and from classical mechanics to statistical mechanics, in a way that is again entirely example-based. We then treat the case of one-dimensional hard spheres and discover the Asakura--Oosawa interaction, the fifth force in nature. We will end with a discussion of the double nature of pressure: kinematic and thermodynamic.
Here are the notes for Lecture 4.
- Here is the Lévy-construction algorithm for the harmonic chain, in nine lines of code. It directly samples configurations.
- Here is the Metropolis algorithm for the harmonic chain, in .. lines of code. It samples configurations using the algorithm developed by Metropolis et al. in (1953).
- Here is the Hamiltonian Monte Carlo algorithm for the harmonic chain, in .. lines of code. It implements the algorithm developed by Duane et al. (1988).
