Quantum Mechanical Fragment Methods Based on Partitioning Atoms or Partitioning Coordinates

Abstract

Conspectus The development of more efficient and more accurate ways to represent reactive potential energy surfaces is a requirement for extending the simulation of large systems to more complex systems, longer-time dynamical processes, and more complete statistical mechanical sampling. One way to treat large systems is by direct dynamics fragment methods. Another way is by fitting system-specific analytic potential energy functions with methods adapted to large systems. Here we consider both approaches. First we consider three fragment methods that allow a given monomer to appear in more than one fragment. The first two approaches are the electrostatically embedded many-body (EE-MB) expansion and the electrostatically embedded many-body expansion of the correlation energy (EE-MB-CE), which we have shown to yield quite accurate results even when one restricts the calculations to include only electrostatically embedded dimers. The third fragment method is the electrostatically embedded molecular tailoring approach (EE-MTA), which is more flexible than EE-MB and EE-MB-CE. We show that electrostatic embedding greatly improves the accuracy of these approaches compared with the original unembedded approaches. Quantum mechanical fragment methods share with combined quantum mechanical/molecular mechanical (QM/MM) methods the need to treat a quantum mechanical fragment in the presence of the rest of the system, which is especially challenging for those parts of the rest of the system that are close to the boundary of the quantum mechanical fragment. This is a delicate matter even for fragments that are not covalently bonded to the rest of the system, but it becomes even more difficult when the boundary of the quantum mechanical fragment cuts a bond. We have developed a suite of methods for more realistically treating interactions across such boundaries. These methods include redistributing and balancing the external partial atomic charges and the use of tuned fluorine atoms for capping dangling bonds, and we have shown that they can greatly improve the accuracy. Finally we present a new approach that goes beyond QM/MM by combining the convenience of molecular mechanics with the accuracy of fitting a potential function to electronic structure calculations on a specific system. To make the latter practical for systems with a large number of degrees of freedom, we developed a method to interpolate between local internal-coordinate fits to the potential energy. A key issue for the application to large systems is that rather than assigning the atoms or monomers to fragments, we assign the internal coordinates to reaction, secondary, and tertiary sets. Thus, we make a partition in coordinate space rather than atom space. Fits to the local dependence of the potential energy on tertiary coordinates are arrayed along a preselected reaction coordinate at a sequence of geometries called anchor points; the potential energy function is called an anchor points reactive potential. Electrostatically embedded fragment methods and the anc...

Explicit solvent simulations of the aqueous oxidation potential and reorganization energy for neutral molecules: gas phase, linear solvent response, and non-linear response contributions

“…Recent advances in fragment quantum mechanical models suggest that improved electronic descriptions of extended systems, at decreased computational cost, may soon be within reach. [124][125][126]…”

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

“…Based on this chemical intuition, the system is divided into many individual subsystems (fragments) and subsequently the properties of the whole system can be obtained by taking a linear combination of the properties of these fragments. Over the past decade, many fragmentation QM methods have been proposed 46 , 47 , including the fragment molecular orbital (FMO) method 48 , the systematic fragmentation method (SFM) 49 , the molecular tailoring approach (MTA) 50 , the molecular fractionation with conjugate caps (MFCC) method 21 , 51 , the adjustable density matrix assembler (ADMA) method 52 , the electrostatically embedded many-body (EE-MB) expansion approach 53 , the explicit polarizatioin (X-Pol) potential 54 .…”

Efficient Calculation of Absorption Spectra in Solution: Approaches for Selecting Representative Solvent Configurations and for Reducing the Number of Explicit Solvent Molecules

“…However, the large system size required to reproduce a real solution and large number of configurations necessary to cover the entire phase space often make full QM methods implausible. Recent advances in full QM methods try to partition the large solution system into molecular fragments, and the computational cost is reduced to some extent [11][12][13][14] . However, the problem of configurational sampling still exists.…”