Proton transfer in bulk water using the full adaptive QM/MM method: integration of solute- and solvent-adaptive approaches

The quantum mechanical/molecular mechanical (QM/MM) method is a hybrid molecular simulation technique that increases the accessibility of local electronic structures of large systems. The technique combines the benefit of accuracy found in the QM method and that of cost efficiency found in the MM method. However, it is difficult to directly apply the QM/MM method to the dynamics of solution systems, particularly for proton transfer. As explained in the Grotthuss mechanism, proton transfer is a structural interconversion between hydronium ions and solvent water molecules. Hence, when the QM/MM method is applied, an adaptive treatment, namely on-the-fly revisions on molecular definitions, is required for both the solute and solvent. Although several solvent-adaptive methods have been proposed, a full adaptive framework, which is an approach that also considers adaptation for solutes, remains untapped. In this paper, we propose a new numerical expression for the coordinates of the excess proton and its control algorithm. Furthermore, we confirm that this method can stably and accurately simulate proton transfer dynamics in bulk water.

On-the-fly determination of active region centers in adaptive-partitioning QM/MM

Quantum mechanics/molecular mechanics (QM/MM) methods are widely used in molecular dynamics (MD) simulations of large systems. By partitioning the system into active and environmental regions and treating them with different levels of theory, QM/MM methods achieve accuracy and efficiency at the same time. Adaptive-partitioning (AP) QM/MM allows the partition of the system to change during the MD simulation, making it possible to simulate processes in which the active and environmental regions exchange atoms or molecules, such as processes in solutions or solids. AP-QM/MM methods usually partition the system according to distances to centers of active regions. For energy-conserving AP-QM/MM methods, these centers are chosen beforehand and remain fixed during the MD simulation, making it difficult to simulate processes in which active regions may occur or vanish. In this paper, I develop an adaptive-center (AC) method that allows on-the-fly determination of the centers of active regions according to any geometrical criterion or any criterion dependent on the potential energy. The AC method is compatible with all existing energy-conserving AP-QM/MM methods, and the resulting potential energy surface is smooth. The application of the AC method is demonstrated with two examples in solid systems.

Extending scaled-interaction adaptive-partitioning QM/MM to covalently bonded systems

Quantum mechanics/molecular mechanics (QM/MM) is the method of choice for atomistic simulations of large systems that can be partitioned into active and environmental regions. Adaptive-partitioning (AP) methods extend the applicability of QM/MM, allowing active regions to change during the simulation. AP methods achieve continuous potential energy surface (PES) by introducing buffer regions in which atoms have both QM and MM characters. Most of the existing AP-QM/MM methods require multiple QM calculations per time step, which can be expensive for systems with many atoms in buffer regions. Although one can lower the computational cost by grouping atoms into fragments, this may not be possible for all systems, especially for applications in covalent solids. The SISPA method [Field, J. Chem. Theory Comput., 2017, 13, 2342] differs from other AP-QM/MM methods by only requiring one QM calculation per time step, but it has the flaw that the QM charge density and wavefunction near the buffer/MM boundary tend to those of isolated atoms/fragments. Besides, regular QM/MM methods for treating covalent bonds cut by the QM/MM boundary are incompatible with SISPA. Due to these flaws, SISPA in its original form cannot treat covalently bonded systems properly. In this work, I show that a simple modification to the SISPA method improves the treatment of covalently bonded systems. I also study the effect of correcting the charge density in SISPA by developing a density-corrected pre-scaled algorithm. I demonstrate the methods with simple molecules and bulk solids.

Analysis of Density Functional Tight Binding with Natural Bonding Orbitals

The description of chemical bonding by the density functional tight binding (DFTB) model is analyzed using natural bonding orbitals (NBOs) and compared to results from density functional theory (B3LYP/aug-cc-pVTZ) calculations. Several molecular systems have been chosen to represent fairly diverse bonding scenarios that include standard covalent bonds, hypervalent interactions, multicenter bonds, metal-ligand interactions (with and without the pseudo-Jahn-Teller effect), and through-space donor-acceptor interactions. Overall, the results suggest that DFTB3/3OB provides physically sound descriptions for the different bonding scenarios analyzed here, as reflected by the general agreement between DFTB3 and B3LYP NBO properties, such as the nature of the NBOs, the magnitudes of natural charges and bond orders, and the dominant donor-acceptor interactions. The degree of ligand-to-metal charge transfer and the ionic nature of pentavalent phosphate are overestimated, likely reflecting the minimal-basis nature of DFTB3/3OB. Moreover, certain orbital interactions, such as geminal interactions, are observed to be grossly overestimated by DFTB3 for hypervalent phosphate and several transition metal compounds that involve copper and nickel. The study indicates that results from NBO analysis can be instructive for identifying electronic structure descriptions at the approximate quantum-mechanical level that require improvement and thus for guiding the systematic improvement of these methods.

Improvement of Performance, Stability and Continuity by Modified Size-Consistent Multipartitioning Quantum Mechanical/Molecular Mechanical Method

For condensed systems, the incorporation of quantum chemical solvent effects into molecular dynamics simulations has been a major concern. To this end, quantum mechanical/molecular mechanical (QM/MM) techniques are popular and powerful options to treat gigantic systems. However, they cannot be directly applied because of temporal and spatial discontinuity problems. To overcome these problems, in a previous study, we proposed a corrective QM/MM method, size-consistent multipartitioning (SCMP) QM/MM and successfully demonstrated that, using SCMP, it is possible to perform stable molecular dynamics simulations by effectively taking into account solvent quantum chemical effects. The SCMP method is characterized by two original features: size-consistency of a QM region among all QM/MM partitioning and partitioning update. However, in our previous study, the performance was not fully elicited compared to the theoretical upper bound and the optimal partitioning update protocol and parameters were not fully verified. To elicit the potential performance, in the present study, we simplified the theoretical framework and modified the partitioning protocol.