Calculations of structure, dynamics and thermodynamics in molecular science all rely on the underlying potential energy surface (PES). Recent advances allow us to visualise this high-dimensional object in a compact fashion, locate global minima efficiently, and sample multistep pathways to obtain rate constants. These methods have been applied to a wide variety of systems, including clusters, glasses and biomolecules, and enable us to treat dynamics on the experimental timescale and beyond.