At the same time, high physical resolution is required to describe molecular and atomic polarization, and chemical reactivity. And yet changes in molecular electronic structure and reactivity are too expensive to describe more than a handful of individual molecules or surface atoms. Classcial MD simulations can be extended to represent bulk systems, but they rely on force fields that are parameterized for the bulk and thus may fail to accurately represent the molecular structure of the interface.