Theoretical and computational chemistry
Potential-energy surfaces and free-energy calculations
Connect molecular geometry to potential-energy landscapes, then use statistical mechanics and enhanced sampling to estimate free-energy differences and barriers.
IntuitionLandscapes and accessible states
A potential-energy surface assigns an energy to each molecular geometry. Its slopes indicate forces and its valleys suggest stable structures, but temperature changes which regions are statistically populated: free energy, not potential energy alone, governs equilibrium probabilities.
SchoolReading an energy diagram
Definition: Potential-energy surface
For fixed electronic state and nuclear coordinates R, V(R) is the potential energy. A stationary point has zero gradient; a local minimum has positive curvature in all unconstrained directions, while a first-order saddle has one downhill direction.
Example: A barrier on a reaction coordinate
Along a simplified one-dimensional coordinate, a reactant minimum is 12 kJ mol⁻¹ below a transition-state maximum. What is the potential-energy barrier from reactant to transition state?
Solution
The barrier is V‡ − Vreactant = 12 kJ mol⁻¹. It is a potential-energy difference, not automatically the activation free energy at finite temperature.
UndergraduateFrom microscopic states to free energy
Definition: Potential of mean force
For a collective variable ξ, the free-energy profile W(ξ) is defined by the equilibrium probability density p(ξ): W(ξ)=−kBT ln p(ξ)+constant. It integrates out all coordinates not explicitly retained in ξ.
The equilibrium populations of two basins determine their free-energy difference, not just the energy at each minimum. Entropic contributions reflect how much configuration space each basin occupies.
| Method | Strategy | Caution |
|---|---|---|
| Umbrella sampling | Bias overlapping windows | Overlap and reweighting matter |
| Metadynamics | Accumulate bias along collective variables | Poor variables hide slow modes |
| Free-energy perturbation | Reweight samples between states | Needs phase-space overlap |
AdvancedBarriers, coordinates, and uncertainty
A reaction coordinate is a reduced description; a poor choice can hide orthogonal slow rearrangements and produce hysteresis. Free-energy methods must demonstrate overlap, reversibility or convergence across windows, and uncertainty estimates; a visually smooth profile alone is insufficient evidence.
ResearchFrontier: rare events and reliable free energies
References
- Understanding Molecular Simulation: From Algorithms to Applications · D. Frenkel, B. Smit, 2002
- Optimized Monte Carlo data analysis · A. M. Ferrenberg, R. H. Swendsen, 1989
- Calculating free energy differences from computer simulations using perturbation theory · R. W. Zwanzig, 1954