Origin and control of superlinear polarizability scaling in chemical potential equalization methods

Abstract

Many common chemical potential equalization (muEq) methods are known to suffer from a superlinear scaling of the polarizability with increasing molecular size that interferes with model transferability and prevents the straightforward application of these methods to large, biochemically relevant molecules. In the present work, we systematically investigate the origins of this scaling and the mechanisms whereby some existing methods successfully temper the scaling. We demonstrate several types of topological charge constraints distinct from the usual single molecular charge constraint that can successfully achieve linear polarizability scaling in atomic charge based equilibration models. We find the use of recently employed charge conservation constraints tied to small molecular units to be an effective and practical approach for modulating the polarizability scaling in atomic muEq schemes. We also analyze the scaling behavior of several muEq schemes in the bond representation and derive closed-form expressions for the polarizability scaling in a linear atomic chain model; for a single molecular charge constraint these expressions demonstrate a cubic dependence of the polarizability on molecular size compared with linear scaling obtainable in the case of the atom-atom charge transfer (AACT) and split-charge equilibration (SQE) schemes. Application of our results to the trans N-alkane series reveals that in certain situations, the AACT and SQE schemes can become unstable due to an indefinite Hessian matrix. Consequently, we discuss sufficient criteria for ensuring stability within these schemes.