REAlgo: Rapid and efficient algorithm for estimating MP2/CCSD energy gradients for large molecular clusters

This work reports the development of an algorithm for rapid and efficient evaluation of energy gradients for large molecular clusters employing correlated methods viz. second-order Møller-Plesset perturbation theory (MP2) theory and couple cluster singles and doubles (CCSD). The procedure segregates the estimation of Hartree-Fock (HF) and correlation components. The HF energy and gradients are obtained by performing a full calculation. The correlation energy is approximated as the corresponding two-body interaction energy. Correlation gradients for each monomer are approximated from the respective monomer-centric fragments comprising its immediate neighbours. The programmed algorithm is explored for the geometry optimization of large molecular clusters using the BERNY optimizer as implemented in the Gaussian suite of software. The accuracy and efficacy of the method are critically probed for a variety of large molecular clusters containing up to 3000 basis functions, in particular large water clusters. The CCSD level geometry optimization of molecular clusters containing ∼800 basis functions employing a modest hardware is also reported.

Hydrogen bond networks of dimethylsulfoxide (DMSO) pentamer

Understanding of clusters of dimethylsulfoxide (DMSO) is important in several applications in Chemistry. Despite its importance, very few studies of DMSO clusters, (DMSO)_n, have been reported in comparison to systems such as water clusters or methanol clusters. In order to provide further understanding of DMSO clusters, we investigated the structures and non-covalent interactions of the (DMSO)_n, n=5. Therefore, the potential energy surface (PES) of the DMSO pentamer has been examined using classical molecular dynamics. The structures generated using classical molecular dynamics are further optimized at the PW6B95D3/aug-cc-pVDZ level of theory. To comprehend the non-covalent bondings in the DMSO pentamer, we carried out a quantum theory of atoms in molecule (QTAIM) analysis. In addition, the effects of temperature on the structural stability is investigated between 20 and 500K. It comes out that seven different kind of non-covalent bondings can be found in DMSO pentamers.

Large-Sized Ammonia Clusters and Solvation Energies of the Proton in Ammonia

The absolute solvation energies (free energies and enthalpies) of the proton in ammonia are used to compute the pK_a of species embedded in ammonia. They are also used to compute the solvation energies of other ions in ammonia. Despite their importance, it is not possible to determine experimentally the solvation energies of the proton in a given solvent. We propose in this work a direct approach to compute the solvation energies of the proton in ammonia from large-sized neutral and protonated ammonia clusters. To undertake this investigation, we performed a geometry optimization of neutral and protonated ammonia 30-mer, 40-mer, and 50 mer to locate stable structures. These structures have been fully optimized at both APFD/6-31++g(d,p) and M06-2X/6-31++g(d,p) levels of theory. An infrared spectroscopic study of these structures has been provided to assess the reliability of our investigation. Using these structures, we have computed the absolute solvation free energy and the absolute solvation enthalpy of the proton in ammonia. It comes out that the absolute solvation free energy of the proton in ammonia is calculated to be -1192 kJ mol^-1 , whereas the absolute solvation enthalpy is evaluated to be -1214 kJ mol^-1 . © 2019 Wiley Periodicals, Inc.

Atlas of putative minima and low-lying energy networks of water clusters n = 3-25

We report a database consisting of the putative minima and ∼3.2 × 10⁶ local minima lying within 5 kcal/mol from the putative minima for water clusters of sizes n = 3-25 using an improved version of the Monte Carlo temperature basin paving (MCTBP) global optimization procedure in conjunction with the ab initio based, flexible, polarizable Thole-Type Model (TTM2.1-F, version 2.1) interaction potential for water. Several of the low-lying structures, as well as low-lying penta-coordinated water networks obtained with the TTM2.1-F potential, were further refined at the Møller-Plesset second order perturbation (MP2)/aug-cc-pVTZ level of theory. In total, we have identified 3 138 303 networks corresponding to local minima of the clusters n = 3-25, whose Cartesian coordinates and relative energies can be obtained from the webpage https://sites.uw.edu/wdbase/. Networks containing penta-coordinated water molecules start to appear at n = 11 and, quite surprisingly, are energetically close (within 1-3 kcal/mol) to the putative minima, a fact that has been confirmed from the MP2 calculations. This large database of water cluster minima spanning quite dissimilar hydrogen bonding networks is expected to influence the development and assessment of the accuracy of interaction potentials for water as well as lower scaling electronic structure methods (such as different density functionals). Furthermore, it can also be used in conjunction with data science approaches (including but not limited to neural networks and machine and deep learning) to understand the properties of water, nature's most important substance.

Pragmatic Many-Body Approach for Economic MP2 Energy Estimation of Molecular Clusters

We propose a procedure, within the many-body analysis (MBA) framework, for an economic yet accurate estimation of the correlated method-based energies of large molecular clusters employing Dunning's augmented basis sets. The basis of the procedure is to segregate the Hartree-Fock ( E_HF) and correlation energy ( E_C) estimations. E_HF is found to differ by tens of millihartrees (mH) from its full-calculation (FC) counterpart on truncating the MBA expansion at the two-body (MBA-2) level. On the contrary, E_C is estimated with smaller error on modest hardware with limited computation time at the (MBA-2) level. In view of this, we adopt a pragmatic method wherein the E_HF (accurate to five decimal places) is taken from the FC, whereas E_C is estimated at the MBA-2 level. This method is applied to a variety of medium to large molecular clusters at the MP2 level. Preliminary results at the CCSD(T) level for (H₂O)₁₆ and (H₂O)₁₇ are also reported with tremendous savings in wall-clock time and resources. The typical errors in MP2 and CCSD(T) energies per monomer are up to 0.1 and 0.2 mH, respectively. Thus the present method, balancing accuracy and computational economy, opens a way for estimating energies of large molecular clusters using correlated theories with large basis sets.

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.