An efficient fragment-based approach for predicting the ground-state energies and structures of large molecules

An efficient fragment-based approach for predicting the ground-state energies and structures of large molecules at the Hartree-Fock (HF) and post-HF levels is described. The physical foundation of this approach is attributed to the "quantum locality" of the electron correlation energy and the HF total energy, which is revealed by a new energy decomposition analysis of the HF total energy proposed in this work. This approach is based on the molecular fractionation with conjugated caps (MFCC) scheme (Zhang, D. W.; Zhang, J. Z. H. J. Chem. Phys. 2003, 119, 3599), by which a macromolecule is partitioned into various capped fragments and conjugated caps formed by two adjacent caps. We find that the MFCC scheme, if corrected by the interaction between non-neighboring fragments, can be used to predict the total energy of large molecules only from energy calculations on a series of small subsystems. The approach, named as energy-corrected MFCC (EC-MFCC), computationally achieves linear scaling with the molecular size. Our test calculations on a broad range of medium- and large molecules demonstrate that this approach is able to reproduce the conventional HF and second-order Moller-Plesset perturbation theory (MP2) energies within a few millihartree in most cases. With the EC-MFCC optimization algorithm described in this work, we have obtained the optimized structures of long oligomers of trans-polyacetylene and BN nanotubes with up to about 400 atoms, which are beyond the reach of traditional computational methods. In addition, the EC-MFCC approach is also applied to estimate the heats of formation for a series of organic compounds. This approach provides an appealing approach alternative to the traditional additivity rules based on either bond or group contributions for the estimation of thermochemical properties.     

A statistical approach to resonance energies of large molecules

For large conjugated molecules, resonance energies of SCF MO quality are not available. Here we outline a method of determining molecular resonance energies by combining a graph theoretical approach to aromaticity with a statistical analysis of random Kekule valence structures. The approach involves construction of random Kekule valence forms and subsequent enumeration of conjugated circuits within each such structure.     

Systematic fragmentation of large molecules by annihilation

A new version of systematic molecular fragmentation is presented which provides a hierarchy of estimates for the energy, and other properties, of large molecules with a computation time that scales linearly with the size of the molecule. This method is combined with an algorithm which ensures that the evaluation of the fragment compositions is efficient for very large molecules. The method is illustrated using protein structures derived from NMR spectroscopy.     

An approach to computing electrostatic charges for molecules

We present an approach for deriving net atomic charges from ab initio quantum mechanical calculations using a least squares fit of the quantum mechanically calculated electrostatic potential to that of the partial charge model. Our computational approach is similar to those presented by Momany [J. Phys. Chem., 82 , 592 (1978)], Smit, Derissen, and van Duijneveldt [Mol. Phys., 37 , 521 (1979)], and Cox and Williams [J. Comput. Chem., 2 , 304 (1981)], but differs in the approach to choosing the positions for evaluating the potential. In this article, we present applications to the molecules H₂O, CH₃OH, (CH₃)₂O, H₂CO, NH₃, (CH₃O)₂PO, deoxyribose, ribose, adenine, 9-CH₃ adenine, thymine, 1-CH₃ thymine, guanine, 9-CH₃ guanine, cytosine, 1-CH₃ cytosine, uracil, and 1-CH₃ uracil. We also address the question of inclusion of “lone pairs,” their location and charge.     

Reduction of the virtual space for coupled-cluster excitation energies of large molecules and embedded systems

We investigate how the reduction of the virtual space affects coupled-cluster excitation energies at the approximate singles and doubles coupled-cluster level (CC2). In this reduced-virtual-space (RVS) approach, all virtual orbitals above a certain energy threshold are omitted in the correlation calculation. The effects of the RVS approach are assessed by calculations on the two lowest excitation energies of 11 biochromophores using different sizes of the virtual space. Our set of biochromophores consists of common model systems for the chromophores of the photoactive yellow protein, the green fluorescent protein, and rhodopsin. The RVS calculations show that most of the high-lying virtual orbitals can be neglected without significantly affecting the accuracy of the obtained excitation energies. Omitting all virtual orbitals above 50 eV in the correlation calculation introduces errors in the excitation energies that are smaller than 0.1 eV. By using a RVS energy threshold of 50 eV, the CC2 calculations using triple-ζ basis sets (TZVP) on protonated Schiff base retinal are accelerated by a factor of 6. We demonstrate the applicability of the RVS approach by performing CC2/TZVP calculations on the lowest singlet excitation energy of a rhodopsin model consisting of 165 atoms using RVS thresholds between 20 eV and 120 eV. The calculations on the rhodopsin model show that the RVS errors determined in the gas-phase are a very good approximation to the RVS errors in the protein environment. The RVS approach thus renders purely quantum mechanical treatments of chromophores in protein environments feasible and offers an ab initio alternative to quantum mechanics/molecular mechanics separation schemes.     

Decomposition of electronic properties and calculation of solvation energies for surface-active molecules

For ions and polar molecules located near a surface, we derived and analyzed the general formulas to calculate their solvation energies. St. Petersburg State University. Translated fromZhurnal Strukturnoi Khimii, Vol. 35, No. 6, pp. 3–12, November-December, 1994. Translated by A. Arbuznikov     

Localized electron pair theory for the calculation of ground state energies of large molecules

The quantum mechanical energy is examined in which groups of one, two, three, and four localized electron pairs found within a molecule are separately computed. From these results, the interaction energies of the electron pairs taken one, two, three, and four at a time form the terms of a convergent molecular mechanics like expansion of the molecular ground state energy. This procedure can be used with any size consistent quantum mechanical method. The computational time for large molecules depends chiefly upon the order needed in the energy expansion to obtain sufficient convergence and not on the particular quantum mechanical method used. Preliminary results within the framework of a semiempirical CNDO/2 model Hamiltonian show at the Hartree–Fock and Møller–Plesset perturbation levels that relative energies converge to within a few tenths of a kcal/mol of the exact values at the four body level for molecules that have little delocalization. In strained ring and aromatic systems, convergence is however not nearly as rapid. Results can be improved somewhat by using larger interacting fragments containing two or more electron pairs over three or more atomic centers. © 1992 by John Wiley & Sons, Inc.     

ESCAPE: A novel approach for a fast estimation of dynamic correlation energies: Application to large organic molecules

Advanced wave function-based quantum chemical ab initio methods, such as CCSD(T), are able to calculate the energies of small- to medium-sized molecules with chemical accuracy. Unfortunately, these methods scale quite unfavorably with the size of the system and are getting too time consuming—and too expensive—for larger molecules. In order to be able to treat larger organic molecules, we propose a novel scheme for a quick and reliable estimate of molecular correlation energies, which we call ESCAPE (ES timation of C orrelA tion energies by P air E nergies). It is based on the pair correlation energies for localized molecular orbitals that have been generated by CCSD[T] and fitted to suitable functional forms. All fit parameters are stored in a large parameter file. Aiming at chemical accuracy (±1 kcal/mol), we have first limited our approach to aliphatic hydrocarbons. The total molecular CCSD[T] correlation energies of a training set of 41 aliphatic hydrocarbons could be reproduced with a mean absolute error (MAE) of 0.56 kcal/mol or 0.11%. A similar accuracy could be obtained for a test set of 11 additional hydrocarbons with up to eight carbon atoms (MAE of 0.65 kcal/mol or 0.09%). In a more critical test, we checked the small energy differences for a set of 13 isomerization reactions. The comparison with experimental data showed that we could reach chemical accuracy as well. Our estimate (MAE of 0.55 kcal/mol) is slightly inferior to the CCSD[T] result (MAE of 0.17 kcal/mol), but superior to SCF, DFT/B3LYP, and DFT/B3LYP + D3. Moreover, in all cases, we obtained the correct sign, that is, the correct equilibrium structure. A similar accuracy could be reached in an application to the three lowest isomers of the C₆₀ molecule. Using the example of a set of eight alcohols, we were able to proof the method's ability for molecules including heteroatoms. Three fast steps are necessary for the application to any aliphatic hydrocarbon or alcohol: (1) An SCF calculation at the selected molecular geometry; it can be fast since a medium size basis set is generally sufficient. (2) The localization of the occupied molecular orbitals and determination of their properties (center of charge and spatial extent). (3) Estimate of the correlation energy using the existing parameter file. © 2019 Wiley Periodicals, Inc.     

Local CC2 electronic excitation energies for large molecules with density fitting

A new local method for the computation of electronic excitation energies of singlet states in extended molecular systems is presented. It is based on the CC2 model and local approximations to the wave functions. In the proposed method the singles excitations are treated nonlocally and local restrictions are imposed on doubles amplitudes only. The accuracy of the new method was tested by calculating several lowest excited states for 14 molecules and comparing them with canonical CC2 values. Deviations of the local excitation energies from the canonical reference values do not exceed 0.05 eV for all test molecules and all states in the lower energy range investigated in this work. The method uses the density-fitting approximation for all two-electron integrals, which considerably simplifies the computational complexity of the individual diagrams. A combination of the local approximations and the powerful density-fitting technique leads to a low-scaling method, capable to treat molecular systems comprised of 100 atoms and more in a basis of a polarized double zeta quality. A test calculation for a system consisting of 127 atoms and 370 active electrons without symmetry is presented to show the efficiency of the new method.     

Molecular spectroscopy for ground-state transfer of ultracold RbCs molecules

We perform one- and two-photon high resolution spectroscopy on ultracold samples of RbCs Feshbach molecules with the aim to identify a suitable route for efficient ground-state transfer in the quantum-gas regime to produce quantum gases of dipolar RbCs ground-state molecules. One-photon loss spectroscopy allows us to probe deeply bound rovibrational levels of the mixed excited (A(1)Σ(+)-b(3)Π)0(+) molecular states. Two-photon dark state spectroscopy connects the initial Feshbach state to the rovibronic ground state. We determine the binding energy of the lowest rovibrational level |v' = 0, J' = 0> of the X(1)Σ(+) ground state to be D = 3811.5755(16) cm(-1), a 300-fold improvement in accuracy with respect to previous data. We are now in the position to perform stimulated two-photon Raman transfer to the rovibronic ground state.     

Convergence radii of the perturbation expansions for the ground-state energies of finite Hubbard models

The ground state energies of finite Hubbard molecules are calculated by numerically solving the Lieb–Wu equations for a complex Hubbard repulsion parameter U. From the positions of the singular points located in the complex plane, the radii of convergence of the perturbation expansions for the ground state energies are determined.     

An efficient approach for the calculation of Franck-Condon integrals of large molecules

A general and efficient approach for the calculation of Franck-Condon integrals (FCIs) of large molecules is presented. In a first step, by exploiting the diagonally dominant and sparse structure of the Duschinsky matrix, a model system is constructed for which the Duschinsky matrix takes a block-diagonal form. For each of these blocks separately, the FCIs are calculated discarding all below a certain threshold. From those integrals retained the FCIs of the model system are obtained by simple multiplication. These serve as an estimate for the FCIs of the exact system which are calculated for those integrals which lie above a certain threshold. By systematically decreasing the threshold, the simulation can be reliably converged to the exact result with an arbitrary accuracy. Using this scheme, a considerable reduction of the number of FCIs which have to be calculated is achieved which leads to an improved scaling behavior of the computational effort with system size. The approach has been tested thoroughly for a set of molecules including difficult cases. For the larger systems a speedup of up to three orders of magnitude compared to an exact calculation is observed while the errors can be kept negligible. With this approach accurate calculations of FCIs are feasible also for large molecules encountered in "real-life" chemistry, especially biochemistry and material science.     

A new fragment-based approach for calculating electronic excitation energies of large systems

We present a new fragment-based scheme to calculate the excited states of large systems without necessity of a Hartree-Fock (HF) solution of the whole system. This method is based on the implementation of the renormalized excitonic method [M. A. Hajj et al., Phys. Rev. B 72, 224412 (2005)] at ab initio level, which assumes that the excitation of the whole system can be expressed by a linear combination of various local excitations. We decomposed the whole system into several blocks and then constructed the effective Hamiltonians for the intra- and inter-block interactions with block canonical molecular orbitals instead of widely used localized molecular orbitals. Accordingly, we avoided the prerequisite HF solution and the localization procedure of the molecular orbitals in the popular local correlation methods. Test calculations were implemented for hydrogen molecule chains at the full configuration interaction, symmetry adapted cluster/symmetry adapted cluster configuration interaction, HF/configuration interaction singles (CIS) levels and more realistic polyene systems at the HF/CIS level. The calculated vertical excitation energies for lowest excited states are in reasonable accordance with those determined by the calculations of the whole systems with traditional methods, showing that our new fragment-based method can give good estimates for low-lying energy spectra of both weak and moderate interaction systems with economic computational costs.     

A refined cluster-in-molecule local correlation approach for predicting the relative energies of large systems

A refined cluster-in-molecule (CIM) method for local correlation calculations of large molecules is presented. In the present work, two new strategies are introduced to further improve the CIM approach: (1) Some medium-range electron correlation energies, which are neglected in the previous CIM approach, are taken into account. (2) A much simpler procedure using only a distance threshold is used to construct various clusters. To cover the medium-range correlation effect as much as possible, some two-atom-centered clusters are built, in addition to one-atom-centered clusters. Our test calculations at the second order perturbation theory (MP2) level show that the refined CIM method can recover about 99.9% of the conventional MP2 correlation energy using an appropriate distance threshold. The accuracy of the present CIM method is capable of providing reliable relative energies of medium-sized systems such as polyalanines with 10 residues, and water molecules with 50 water molecules. For polyalanines with up to 30 residues, we have demonstrated that the computational cost of the CIM-MP2 calculation increases linearly with the molecular size, but the required memory and disc-space do not need to increase for large systems. The improved CIM method has been used to compute the relative energy of ice-like (H(2)O)(96) clusters (with 2400 basis functions) and to predict the dimerization energy of a double-helical foldamer (with 2330 basis functions). The present CIM method is expected to be a practical local correlation method for describing the relative energies of large systems.     

Long nanotube dielectric properties used as sensors of large molecules: a semicontinuum approach

A semicontinuum approach on the basis of an effective polarizability tensor per length and radius units is used to describe the dielectric response of a long single wall nanotube to the adsorption of an extended molecule. Changes in the permittivity ratio of the nanotube+molecule over the nanotube alone, which are directly connected to frequency shifts of the nanotube in a resonator configuration due to the presence of the molecule, provide a test of sensitivity of the system. The behavior of this ratio is analyzed for linear and circular geometries of the molecule, as a function of the tube characteristics (length and radius) and of the molecular size and polarizability distribution. Extension to three dimensional systems with a large set of polarizable centers is discussed in terms of self-polarization of the centers and morphology of the surface of the sensed system.     

Ab initio quality one-electron properties of large molecules: development and testing of molecular tailoring approach

The development of a linear-scaling method, viz. "molecular tailoring approach" with an emphasis on accurate computation of one-electron properties of large molecules is reported. This method is based on fragmenting the reference macromolecule into a number of small, overlapping molecules of similar size. The density matrix (DM) of the parent molecule is synthesized from the individual fragment DMs, computed separately at the Hartree-Fock (HF) level, and is used for property evaluation. In effect, this method reduces the O(N(3)) scaling order within HF theory to an n.O(N'(3)) one, where n is the number of fragments and N', the average number of basis functions in the fragment molecules. An algorithm and a program in FORTRAN 90 have been developed for an automated fragmentation of large molecular systems. One-electron properties such as the molecular electrostatic potential, molecular electron density along with their topography, as well as the dipole moment are computed using this approach for medium and large test chemical systems of varying nature (tocopherol, a model polypeptide and a silicious zeolite). The results are compared qualitatively and quantitatively with the corresponding actual ones for some cases. This method is also extended to obtain MP2 level DMs and electronic properties of large systems and found to be equally successful.