where F QM,nx • T -1 • R m represents the contribution of the MM part, which indicates that the two perturbations both act on the QM region (i.e., F QM,nx ) in the current framework. It is incorporated into the dipole moment matrix and its derivative matrix in the last line.