where c is the center of mass of the system and K a user-defined variable with units kcal/ Å2 , which scales the pressure. The resulting force field corresponds to putting the nuclei into a parabolic potential centered at c. The compression force field F ext is added to the nuclear gradient during a geometry optimization, converging at a stationary point on the pressure-modified PES.