Dynamical effects in ab initio NMR calculations: Classical force fields fitted to quantum forces

Abstract

NMR chemical shifts for an L-alanine molecular crystal are calculated using ab initio plane wave density functional theory. Dynamical effects including anharmonicity may be included by averaging chemical shifts over an ensemble of structural configurations generated using molecular dynamics (MD). The time scales required mean that ab initio MD is prohibitively expensive. Yet the sensitivity of chemical shifts to structural details requires that the methodologies for performing MD and calculating NMR shifts be consistent. This work resolves these previously competing requirements by fitting classical force fields to reproduce ab initio forces. This methodology is first validated by reproducing the averaged chemical shifts found using ab initio molecular dynamics. Study of a supercell of L-alanine demonstrates that finite size effects can be significant when accounting for dynamics.