1. The many-electron problem and QMC routes
Position in the course: Lesson 1 of 10. Complete the preceding derivation and use the explained exercises to check understanding.
Prerequisites: Wavefunctions; variational principle; Monte Carlo
Learning goal: explain the mathematical steps, reproduce the analytical examples, and state the conditions under which the conclusions hold.
1. What is being integrated?
In the Born–Oppenheimer approximation, nuclei are fixed parameters and electronic coordinates form \(R=(\mathbf r_1,\ldots,\mathbf r_N)\), a point in \(3N\) dimensions. Atomic units set \(\hbar=m_e=e=4\pi\epsilon_0=1\). For nonrelativistic electrons and local Coulomb interactions,
\(E_{NN}\) is a constant for fixed nuclei. The wavefunction includes spin as well as coordinates; a fixed spin assignment is often used for a spin-independent Hamiltonian, with antisymmetry enforced within same-spin sectors. Energy differences require a consistent nuclear term. QMC samples electronic configurations rather than single-electron orbitals; random sampling does not remove electron correlation or antisymmetry from the model.
2. Three distinct targets
VMC evaluates an expectation in a chosen trial state, using \(|\Psi_T|^2\) as sampling density. DMC applies imaginary-time projection to improve an amplitude, with a sign or nodal constraint for fermions. Finite-temperature path-integral Monte Carlo instead samples a density matrix trace \(\operatorname{Tr}e^{-\beta H}\). These approaches share stochastic integration but solve different mathematical problems. Increasing VMC samples improves an integral, not the trial wavefunction; increasing DMC projection time does not by itself remove fixed-node bias.
3. Derive the variational bound before sampling
Expand a normalized admissible trial state in eigenstates: \(\Psi_T=\sum_nc_n\Psi_n\), with \(\sum_n|c_n|^2=1\). Then
This requires a self-adjoint bounded-below Hamiltonian and a trial state in the appropriate quadratic-form domain and symmetry sector. The result is an exact expectation bound, not a guarantee that every noisy sample mean lies above \(E_0\). A restricted symmetry sector may have a different lowest energy. Nonlocal pseudopotential treatments and fixed-phase methods require separately stated qualifications.
4. A worked two-component state
Suppose \(E_0=-1\) and \(E_1=-0.5\) hartree, with excited weight \(|c_1|^2=0.2\). The trial energy is \(0.8(-1)+0.2(-0.5)=-0.9\) hartree. Reducing the excited weight to \(0.1\) improves energy to \(-0.95\) hartree. Averaging the original state more precisely cannot produce that improvement in its exact expectation. The construction illustrates a distinction between systematic wavefunction error and finite sampling error, not an actual calculation for a material.
Dimensionality reflects information content: a one-electron density discards much of a correlated many-body amplitude. QMC often retains an explicit many-electron trial function even when orbitals come from a one-particle method. The calculation is neither automatically exact nor simply DFT with random numbers. Identify the trial ansatz and Hamiltonian independently, and match boundary conditions, nuclear terms, potentials, and units before comparing energies.
5. Exercises and explained solutions
Exercise: Does replacing a Coulomb Hamiltonian by a pseudopotential preserve the same exact reference energy?
Solution
No. It changes the modeled Hamiltonian. A variational or convergence claim must identify which Hamiltonian is being bounded or solved; comparisons to all-electron values require a controlled approximation discussion.
Exercise: Is a QMC walker a physical electron trajectory?
Solution
Usually no. A walker is a full many-electron configuration used to represent a distribution or amplitude. Its sampling evolution is not real-time quantum dynamics.
Original teaching schematic of the mathematics or algorithm; it is not simulation or experimental data.
6. Sources and connections
Related theory: Classical Monte Carlo · Molecular methods · Density functional theory
Software connection: CP2K · Quantum ESPRESSO · Gaussian
These software courses provide related background on energies, orbitals, or convergence management; they do not imply that the Monte Carlo or QMC examples on this page were executed there.
7. Related theory and practice
Quantum mechanics · Molecular methods · Monte Carlo · Molecular dynamics · QMCPACK