A scheme for rapid evaluation of the intermolecular three-body polarization effect in water clusters

The ability to accurately and rapidly evaluate the intermolecular many-body polarization effect of the water system is very important for computer simulations of biomolecule in aqueous. In this paper, a scheme is proposed based on the polarizable dipole-dipole interaction model and used to rapidly estimate the intermolecular many-body polarization effect in water clusters. We use a bond-dipole-based polarization function to evaluate the polarization energy. We regard two OH bonds of a water molecule as two bond-dipoles and set the permanent OH bond-dipole moment of a water molecule to be 1.51 Debye. We estimate the induced OH bond-dipole moment via a simple formula in which only one correction factor is needed. This scheme is then applied to tens of water clusters to calculate the three- and four-body interaction energies. The three-body interaction energies of 93 water clusters produced by our scheme are compared with those produced by the counterpoise-corrected CCSD(T)/aug-cc-pVDZ, MP2/aug-cc-pVDZ, M06-2X/jul-cc-pVTZ methods, by the AMOEBApro13, iAMOEBA, AMOEBA+, AMOEBA+(CF) methods, and by the MB-pol method. The four-body interaction energies of 47 water clusters yielded by our scheme are compared with those yielded by the counterpoise-corrected MP2/aug-cc-pVDZ and M06-2X/ jul-cc-pVTZ methods, by the AMOEBApro13, AMOEBA+, AMOEBA+(CF) methods, and by the MB-pol method. The comparison results show that the scheme proposed in this paper can reproduce the counterpoise-corrected CCSD(T)/aug-cc-pVDZ three-body interaction energies and reproduce the counterpoise-corrected MP2/aug-cc-pVDZ four-body interaction energies both accurately and efficiently. We anticipate the scheme proposed here can be useful for computer simulations of liquid water and aqueous solutions.

Effects of Morphology and Solvent/Temperature on THz Spectra: Take Nucleosides as Example

Water molecules were easy to combine with organic molecules and embed into the lattice of solid molecules to form a hydrate. Compared with anhydrous compounds, a hydrate has completely different physical and chemical properties. In this paper, terahertz (THz) spectra of five nucleosides in the solid and liquid phases were studied experimentally by Fourier-transform infrared spectroscopy (FTIR) in the frequency of 0.5-9 THz. In addition, the lattice energy, geometric structure, and vibration spectrum of the molecular crystal of the nucleosides were analyzed theoretically by the generalized energy-based fragmentation approach under periodic boundary conditions (denoted as PBC-GEBF). Furthermore, different nucleoside molecular morphology (monomer, polymer, and crystal), solvent (implicit and explicit water), and temperature/theoretical model effect on the THz spectra were mainly investigated. It was found that in the low-frequency band, the vibrational modes were generally originated from the collective vibration of all molecules involved (more than 99% of them were vibration; only less than 1% of them were rotation and translation), which can reflect the molecular structure and spatial distribution of different substances. The Gibbs free energy of thymidine monomer, dimer, tetramer, and crystal was studied. It was found that the cell-stacking energy had the greatest influence on the spectrum, indicating that only the crystal structure constrained by the periodic boundary conditions could well describe the experimental results. In addition, hydrophobic forces dominated the formation of new chemical bonds and strong inter-molecular interactions; the free water had little contribution to the THz spectrum of nucleosides, while crystalline water had a great influence on the spectrum.

Terahertz spectra and weak intermolecular interactions of nucleosides or nucleoside drugs

In this paper, terahertz (THz) spectra of four DNA nucleosides (Adenosine, Thymidine, Cytidine and Guanosine) and two nucleoside derivatives (Ribavirin and Entecavir, first time reported) in the solid phase were studied experimentally by Fourier Transform Infrared Spectroscopy (FTIR) in the frequency of 1-10 THz. The lattice energy, geometric structure, vibration spectrum of them were analyzed theoretically by the generalized energy-based fragmentation approach under periodic boundary conditions (denoted as PBC-GEBF) and the density functional theory (DFT). The intra- and inter-molecular weak interactions corresponding to the vibrational modes of the crystal, polymer and monomer were obtained, with the help of the potential energy distribution (PED) and reduced density gradient (RDG) methods. It was found that the sum of electronic and thermal free energies increased from the monomer to polymer, and from the polymer to crystal. For example, the inter-molecular interaction energy from the monomer to dimer of adenosine increased 6.969 kcal/mol, and that from the dimer to crystal (the periodic boundary conditions were considered) increased 666.792 kcal/mol. Therefore, only the crystal structure constrained the periodic boundary conditions could well describe the experimental results, although the former scholars chose the monomer or polymer as the initial configuration due to the limitation of computing resources and methods. In THz band, the vibrational modes were generally originated from the collective vibration (more than 99% of them were vibration, only less than 1% of them were rotation and translation) of all molecules involved, which could reflect the molecular structure and spatial distribution of different substances. In order to accurately identify the spectra, we studied the location, type and contribution of all weak interactions, and found that the strong characteristic peaks corresponding to the strong hydrogen bonds came from inter-molecular, while the weak hydrogen bonds mainly originated from intra- and inter-molecular, the out-of-plane bending made the largest contribution, accounting for more than 90%. Furthermore, taking guanine, guanosine and two guanosine derivatives (Ribavirin and Entecavir) as examples, the differences of weak interaction among them caused by different molecular configuration, arrangement and substituent position were studied, and the fundamental reason of THz spectrum change was found. This research can lay a foundation for crystal engineering, supramolecular chemistry, molecular recognition and self-assembly, protein-ligand interaction, etc.

A unified and flexible formulation of molecular fragmentation schemes

We present a flexible formulation for energy-based molecular fragmentation schemes. This framework does not only incorporate the majority of existing fragmentation expansions but also allows for flexible formulation of novel schemes. We further illustrate its application in multi-level approaches and for electronic interaction energies. For the examples of small water clusters, a small protein, and protein-protein interaction energies, we show how this flexible setup can be exploited to generate a well-suited multi-level fragmentation expansion for the given case. With such a setup, we reproduce the electronic protein-protein interaction energy of ten different structures of a neurotensin and an extracellular loop of its receptor with a mean absolute deviation to the respective super-system calculations below 1 kJ/mol.

Fragment-based quantum mechanical approach to biomolecules, molecular clusters, molecular crystals and liquids

To study large molecular systems beyond the system size that the current state-of-the-art ab initio electronic structure methods could handle, fragment-based quantum mechanical (QM) approaches have been developed over the past years, and proved to be efficient in dealing with large molecular systems at various ab initio levels. According to the fragmentation approach, a large molecular system can be divided into subsystems (fragments), and subsequently the property of the whole system can be approximately obtained by taking a proper combination of the corresponding terms of individual fragments. Therefore, the standard QM calculation of a large system could be circumvented by carrying out a series of calculations on small fragments, which significantly promotes computational efficiency. The electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) method is one of the fragment-based QM approaches which has been developed by our research group in recent years. This Perspective presents the theoretical framework of this fragmentation method and its applications in biomolecules, molecular clusters, molecular crystals and liquids, including total energy calculation, protein-ligand/protein binding affinity prediction, geometry optimization, vibrational spectrum simulation, ab initio molecular dynamics simulation, and prediction of excited-state properties.

Fantasy versus reality in fragment-based quantum chemistry

Since the introduction of the fragment molecular orbital method 20 years ago, fragment-based approaches have occupied a small but growing niche in quantum chemistry. These methods decompose a large molecular system into subsystems small enough to be amenable to electronic structure calculations, following which the subsystem information is reassembled in order to approximate an otherwise intractable supersystem calculation. Fragmentation sidesteps the steep rise (with respect to system size) in the cost of ab initio calculations, replacing it with a distributed cost across numerous computer processors. Such methods are attractive, in part, because they are easily parallelizable and therefore readily amenable to exascale computing. As such, there has been hope that distributed computing might offer the proverbial "free lunch" in quantum chemistry, with the entrée being high-level calculations on very large systems. While fragment-based quantum chemistry can count many success stories, there also exists a seedy underbelly of rarely acknowledged problems. As these methods begin to mature, it is time to have a serious conversation about what they can and cannot be expected to accomplish in the near future. Both successes and challenges are highlighted in this Perspective.

Variational Formulation of the Generalized Many-Body Expansion with Self-Consistent Charge Embedding: Simple and Correct Analytic Energy Gradient for Fragment-Based ab Initio Molecular Dynamics

The many-body expansion (MBE) and its extension to overlapping fragments, the generalized (G)MBE, constitute the theoretical basis for most fragment-based approaches for large-scale quantum chemistry. We reformulate the GMBE for use with embedding charges determined self-consistently from the fragment wave functions, in a manner that preserves the variational nature of the underlying self-consistent field method. As a result, the analytic gradient retains the simple "sum of fragment gradients" form that is often assumed in practice, sometimes incorrectly. This obviates (without approximation) the need to solve coupled-perturbed equations, and we demonstrate stable, fragment-based ab initio molecular dynamics simulations using this technique. Energy conservation fails when charge-response contributions to the Fock matrix are neglected, even while geometry optimizations and vibrational frequency calculations may yet be accurate. Stable simulations can be recovered by means of straightforward modifications introduced here, providing a general paradigm for fragment-based ab initio molecular dynamics.

Structures, relative stability and binding energies of neutral water clusters, (H2O)2–30

We have revised the structures of neutral water clusters, (H₂O)_n=2–30, with the affordable M06-2X functional, presenting up to 25 isomers for each cluster size.

DFT study of water adsorption on lignite molecule surface

High moisture content is a main characteristic of low-rank coal, such as lignite. Numerous oxygen containing functional groups in lignite make it represent some special properties, and these functional groups affect the adsorption mechanisms of water molecules on lignite surface. This study reports some typical water · · · lignite conformations, along with a detailed analysis of the geometry, electrostatic potential distribution, reduced density gradient of interaction, and interaction energy decomposition. The results show that water molecules tend to aggregate around functional groups, and hydrogen bonds play a dominant role in the interaction. The adsorption energy of water cluster on lignite surface is larger than that of isolated water molecule, a good linear relationship between the interaction distance and adsorption energy of layers has been found. Since water is a polar molecule, the local minima and maxima of electrostatic potential in conformations increase along with more water adsorbing on lignite surface. Reduced density gradient analysis shows that H-bonds, van der Waals interaction, and a little steric make up the interaction between water cluster and lignite molecule. In these studied conformations which mainly are H-bond complexes, electrostatic and exchange repulsion play a dominant role, whereas polarization and dispersion make relatively small contribution to the interaction. Attractive and repulsive interaction both affect the stability of water · · · lignite conformations.