Machine learning based interatomic potential energy surface (PES) models are revolutionizing the field of molecular modeling. However, although much faster than electronic structure schemes, these models suffer from a lower efficiency as compared to
typical empirical force fields due to more sophisticated computations involved. Herein, we report a model compression scheme for boosting the performance of the Deep Potential (DP) model, a deep learning based PES model. This scheme, we call DP Compress, is an efficient post-processing step after the training of DP models (DP Train). DP Compress combines several DP-specific compression techniques, which typically speed up DP- based molecular dynamics simulations by an order of magnitude faster, and consume an order of magnitude less memory. We demonstrate that DP Compress is sufficiently accurate by testing a variety of physical properties of Cu, H2O, and Al-Cu-Mg systems. DP Compress applies to both CPU and GPU machines and is publicly available at https://github.com/deepmodeling/deepmd-kit.
During compression of a water dimer calculated with high-precision first-principles methods, the trends of H-bond and O-H bond lengths show quantum effect of the electronic structure. We found that the H-bond length keeps decreasing, while the O-H bo
nd length increases up to the stable point and decreases beyond it when the water dimer is further compressed. The remarkable properties are different from those observed in most previous researches which can be understood and extrapolated through classical simulation. The observations can be explained by the decrease in orbital overlap and change in the exchange repulsion interaction between water monomers. The dominant interaction between water monomers changes from electrostatic interaction to exchange repulsion at the turning point of the O-H bond length when the O...O distance is decreased. These findings highlight the quantum effect on the hydrogen bond in extreme conditions and play an important role in the recognition of the hydrogen bond structure and mechanism.
Since roaming was found as a new but common reaction path of isomerization, many of its properties, especially those of roaming transition state (TSR), have been studied recently on many systems. However, the mechanism of roaming is still not clear a
t the atomic level. In this work, we used first-principles calculations to illustrate the detailed structure of TSR in an internal isomerization process of nitrobenzene. The calculations distinctively show its nature of antiferromagnetic coupling between two roaming fragments. Moreover, the effect of dispersion is also revealed as an important issue for the stability of the TSR. Our work provides a new insight from the view of electronic structure towards the TSR and contributes to the basic understanding of the roaming systems.
Hydrogen bond (H-bond) covalency has recently been observed in ice and liquid water, while the penetrating molecular orbitals (MOs) in the H-bond region of most typical water dimer system, (H2O)2, have also been discovered. However, obtaining the qua
ntitative contribution of these MOs to the H-bond interaction is still problematic. In this work, we introduced the orbital-resolved electron density projected integral (EDPI) along the H-bond to approach this problem. The calculations show that, surprisingly, the electronic occupied orbital (HOMO-4) of (H2O)2 accounts for about 40% of the electron density at the bond critical point. Moreover, the charge transfer analysis visualizes the electron accumulating effect of the orbital interaction within the H-bond between water molecules, supporting its covalent-like character. Our work expands the classical understanding of H-bond with specific contributions from certain MOs, and will also advance further research into such covalency and offer quantitative electronic structure insights into intermolecular systems.
The self-consistent charge density functional tight-binding (DFTB) theory is a useful tool for realizing the electronic structures of large molecular complex systems. In this study, we analyze the electronic structure of C61, formed by fullerene C60
with a carbon adatom, using the fully localized limit and pseudo self-interaction correction methods of DFTB to adjust the Hubbard U parameter (DFTB+U). The results show that both the methods used to adjust U can significantly reduce the molecular orbital energy of occupied states localized on the defect carbon atom and improve the gap between highest occupied molecular orbital(HOMO) and lowest unoccupied molecular orbital(LUMO) of C61. This work will provide a methodological reference point for future DFTB calculations of the electronic structures of carbon materials.
The covalent-like characteristics of hydrogen bonds offer a new perspective on intermolecular interactions. Here, using density functional theory and post-Hartree-Fock methods, we reveal that there are two bonding molecular orbitals (MOs) crossing th
e O and H atoms of the hydrogen-bond in water dimer. Energy decomposition analysis also shows a non-negligible contribution of the induction term. These results illustrate the covalent-like character of the hydrogen bond between water molecules, which contributes to the essential understanding of ice, liquid water, related materials, and life sciences.
