Linear scaling local correlation approach for solving the coupled cluster equations of large systems

A linear scaling local correlation approach is proposed for approximately solving the coupled cluster doubles (CCD) equations of large systems in a basis of orthogonal localized molecular orbitals (LMOs). By restricting double excitations from spatially close occupied LMOs into their associated virtual LMOs, the number of significant excitation amplitudes scales only linearly with molecular size in large molecules. Significant amplitudes are obtained to a very good approximation by solving the CCD equations of various subsystems, each of which is made up of a cluster associated with the orbital indices of a subset of significant amplitudes and the local environmental domain of the cluster. The combined effect of these two approximations leads to a linear scaling algorithm for large systems. By using typical thresholds, which are designed to target an energy accuracy, our numerical calculations for a wide range of molecules using the 6-31G or 6-31G* basis set demonstrate that the present local correlation approach recovers more than 98.5% of the conventional CCD correlation energy.     

A localized molecular-orbital assembler approach for Hartree-Fock calculations of large molecules

We describe an alternative fragment-based method, the localized molecular-orbital assembler method, for Hartree-Fock (HF) calculations of macromolecules. In this approach, a large molecule is divided into many small-size fragments, each of which is capped by its local surroundings. Then the conventional HF calculations are preformed on these capped fragments (or subsystems) and the canonical molecular orbitals of these systems are transferred into localized molecular orbitals (LMOs). By assembling the LMOs of these subsystems into a set of LMOs of the target molecule, the total density matrix of the target molecule is constructed and correspondingly the HF energy or other molecular properties can be approximately computed. This approach computationally achieves linear scaling even for medium-sized systems. Our test calculations with double-zeta and polarized double-zeta basis sets demonstrate that the present approach is able to reproduce the conventional HF energies within a few millihartrees for a broad range of molecules.     

A fragment energy assembler method for Hartree-Fock calculations of large molecules

We present a fragment energy assembler approach for approximate Hartree-Fock (HF) calculations of macromolecules. In this method, a macromolecule is divided into small fragments with appropriate size, and then each fragment is capped by its neighboring fragments to form a subsystem. The total energy of the target system is evaluated as the sum of the fragment energies of all fragments, which are available from conventional HF calculations on all subsystems. By applying the method to a broad range of molecules, we demonstrate that the present approach could yield satisfactory HF energies for all studied systems.     

A natural linear scaling coupled-cluster method

It is shown that using an appropriate localized molecular orbital (LMO) basis, one is able to calculate coupled-cluster singles and doubles (CCSD) wave functions and energies for very large systems by performing full CCSD calculations on small subunits only. This leads to a natural linear scaling coupled-cluster method (NLSCC), in which total correlation energies of extended systems are evaluated as the sum of correlation energy contributions from individual small subunits within that system. This is achieved by defining local occupied orbital correlation energies. These are quantities, which in the LMO basis become transferable between similar molecular fragments. Conventional small scale existing molecular CCSD codes are all that is needed, the local correlation effect being simply transmitted via the appropriate LMO basis. Linear scaling of electronic correlation energy calculations is thus naturally achieved using the NLSCC approach, which in principle can treat nonperiodic extended systems of infinite basis set size. Results are shown for alkanes and several polyglycine molecules and the latter compared to recent results obtained via an explicit large scale LCCSD calculation. (c) 2004 American Institute of Physics.     

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.     

Kohn-Sham Density Matrix and the Kernel Energy MethodSCI北大核心CSCD

The kernel energy method(KEM) has been shown to provide fast and accurate molecular energy calculations for molecules at their equilibrium geometries.KEM breaks a molecule into smaller subsets,called kernels,for the purposes of calculation.The results from the kernels are summed according to an expression characteristic of KEM to obtain the full molecule energy.A generalization of the kernel expansion to density matrices provides the full molecule density matrix and orbitals.In this study,the kernel expansion for the density matrix is examined in the context of density functional theory(DFT) Kohn-Sham(KS) calculations.A kernel expansion for the one-body density matrix analogous to the kernel expansion for energy is defined,and is then converted into a normalizedprojector by using the Clinton algorithm.Such normalized projectors are factorizable into linear combination of atomic orbitals(LCAO) matrices that deliver full-molecule Kohn-Sham molecular orbitals in the atomic orbital basis.Both straightforward KEM energies and energies from a normalized,idempotent density matrix obtained from a density matrix kernel expansion to which the Clinton algorithm has been applied are compared to reference energies obtained from calculations on the full system without any kernel expansion.Calculations were performed both for a simple proof-of-concept system consisting of three atoms in a linear configuration and for a water cluster consisting of twelve water molecules.In the case of the proof-of-concept system,calculations were performed using the STO-3 G and6-31 G(d,p) bases over a range of atomic separations,some very far from equilibrium.The water cluster was calculated in the 6-31 G(d,p) basis at an equilibrium geometry.The normalized projector density energies are more accurate than the straightforward KEM energy results in nearly all cases.In the case of the water cluster,the energy of the normalized projector is approximately four times more accurate than the straightforward KEM energy result.The KS density matrices of this study are applicable to quantum crystallography.     

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.     

The tautomers of uracil: A local correlation treatment

The relative energies of all six uracil tautomers have been determined at the MP 4(SDQ )/6-31G** level, using both conventional correlation theory and the Local Correlation method. Geometries were optimized at the SCF /6-31G* level with offset forces. Comparison of our energies with energies from structures optimized at the SCF level supports the conclusion that offset forces are an advantageous alternative to correlated geometry optimization. The Local Correlation method compares very well with conventional Møller–Plesset theory, recovering at least 98.4% of the conventional correlated energy in all cases. More importantly, the relative energies also show good agreement with the conventional results, even for these delocalized systems. CPU timings show a substantial computational savings for the Local Correlation method over the conventional method. The results of the local method using Boys localization are compared with those using Pipek–Mezey localization. The dioxo tautomer ( 1 ) is predicted to be the most stable. The ( 1 )–( 3 ) and ( 1 )–( 4 ) energy differences are found to be within the bounds estimated from experimental work. © 1993 John Wiley & Sons, Inc.     

Semiempirical PM5 molecular orbital study on chlorophylls and bacteriochlorophylls: comparison of semiempirical, ab initio, and density functional results

The semiempirical PM5 method has been used to calculate fully optimized structures of magnesium-bacteriochlorin, magnesium-chlorin, magnesium-porphin, mesochlorophyll a, chlorophylls a, b, c(1), c(2), c(3), and d, and bacteriochlorophylls a, b, c, d, e, f, g, and h with all homologous structures. Hartree-Fock/6-31G* ab initio and density functional B3LYP/6-31G* methods were used to optimize structures of methyl chlorophyllide a, chlorophyll c(1), and methyl bacteriochlorophyllides a and c for comparison. Spectroscopic transition energies of the chromophores and their 1:1 or 1:2 solvent complexes were calculated with the Zindo/S CIS method. The self-consistent reaction field model was used to estimate solvent shifts. The PM5 calculations predict planar structure of the porphyrin ring and central position of the four coordinated magnesium atoms in all pigments studied, in accord with the experimental, ab initio, and density functional results, a significant improvement as compared to the older semiempirical PM3 approach. Only small differences in PM5 and B3LYP/6-31G* or Hartree-Fock/6-31G* minimum energy geometries of the reference molecules were observed. Calculations show that in 1:1 solvent complexes, where the magnesium atom is five coordinated, the magnesium atom is shifted out of the plane of the porphyrin ring towards the solvent molecule, while the hexa coordinated 1:2 complexes are again planar. The PM5 method gives atomic charges that are comparable with those obtained from the Hartree-Fock/6-31G* and B3LYP/6-31G* calculations. The single point ZINDO/S CIS calculations with PM5 minimum energy structure gave excellent correlations between calculated and experimental transition energies of the chlorophylls and bacteriochlorophylls studied. Such correlations may be used for prediction of transition energies of the chromophores in protein binding sites. Calculations also predict existence of dark electronic states below the main Soret absorption band in all chromophores studied. The results suggest that the semiempirical PM5 method is a fairly reliable and computationally efficient method in predicting molecular parameters of porphyrin-like molecules.     

Generalized energy-based fragmentation approach for computing the ground-state energies and properties of large molecules

We present a generalized energy-based fragmentation (GEBF) approach for approximately predicting the ground-state energies and molecular properties of large molecules, especially those charged and polar molecules. In this approach, the total energy (or properties) of a large molecule can be approximately obtained from energy (or properties) calculations on various small subsystems, each of which is constructed to contain a certain fragment and its local surroundings within a given distance. In the quantum chemistry calculation of a given subsystem, those distant atoms (outside this subsystem) are modeled as background point charges at the corresponding nuclear centers. This treatment allows long-range electrostatic interaction and polarization effects between distant fragments to be taken into account approximately, which are very important for polar and charged molecules. We also propose a new fragmentation scheme for constructing subsystems. Our test calculations at the Hartree-Fock and second-order M?ller-Plesser perturbation theory levels demonstrate that the approach could yield satisfactory ground-state energies, the dipole moments, and static polarizabilities for polar and charged molecules such as water clusters and proteins.     

Effective computational modeling of constitutional isomerism and aggregation states of explicit solvates of lithiated phenylacetonitrile

We present the first calculations which accurately account for the position of metalation and aggregation state of lithiated nitriles. Solvation is found to be a key determinant of structure. Five known solvates of lithiated phenylacetonitrile were examined computationally to determine the minimum level of theory required to reproduce the observed X-ray and multinuclear NMR structures. In all cases Hartree-Fock 3-21G energies of explicit solvates calculated at PM3 geometries correctly predict the observed N-lithiated constitutional isomer. Selected density functional theory (B3LYP/6-31+G*//PM3) energy calculations reproduce this trend. We also show that 3-21G//PM3 calculations which do not include explicit solvent molecules, or which include water as a model for diethyl ether, may lead to incorrect predictions of the preferred constitutional isomer. 3-21G//PM3 energies also adequately account for observed aggregation states of the TMEDA, diethyl ether, and THF solvates. Finally, calculations of THF-solvated monomers up to the B3LYP/6-31+G*//B3LYP/6-31+G level indicate a significant (6.8 kcal/mol) preference for N-lithiation.     

Quantum mechanical investigation of the inner-sphere reorganization energy of cyclooctatetraene/cyclooctatetraene radical anion. Part I

The inner-sphere reorganization energy of the electron self-exchange of the couple cyclooctatetraene/cyclooctatetraene radical anion has been investigated by quantum mechanical calculations. The more stable Jahn Teller distorted B2g conformation of the radical anion has been used in this study. Two different theories have been applied in this first part. The harmonic approximation in the classical Marcus scheme has been modified by using projected force constants, which are obtained from the complete force constant matrix and the geometry changes of the molecule during the ET (introduced by Mikkelsen). A different approach (introduced by Nelsen) combines the different energies of the neutral and radical anion with and without relaxation corresponding to the vertical ionization potential and the vertical electron affinity. The electronic energies of the neutral molecule and the radical anion differ dramatically applying three different levels of quantum mechanical calculations (UAM1, UB3LYP, PMP2 with three different basis sets with and without diffuse functions). Nevertheless the Nelsen method gives almost consistent results for the inner-sphere reorganization energies: 120.1 kJ/mol for semiempirical UAM1 method, 159.3 kJ/mol, 156.4 kJ/mol and 158.3 kJ/mol for density functional UB3LYP/6-31G*, UB3LYP/6-31++G* and UB3LYP/AUG-cc-pVDZ calculations and 192.5 kJ/mol for ab-initio PMP2/6-31G* investigations, respectively. These values are in agreement with earlier experimental work supposing the total reorganization energy to be larger than 38 kcal/mol assuming an electron self-exchange rate of 10(4) M(-1) s(-1). The simple harmonic approximation of Marcus relation has not yet been applied for a molecule like cyclooctatetraene with large torsional geometry changes. Using the projected force constants after scaling, considerably different results for the inner-sphere reorganization energy have been calculated: 738.1 kJ/mol for the UB3LYP/6-31G*, 743.3 kJ/mol for UB3LYP/6-31++G* and 759.1 kJ/mol for UB3LYP/AUG-cc-pVDZ level of theory. Comparison with our concentration dependent EPR experiments are controversial to the earlier experimental results, but the latter supports the assumption that the electron self-exchange occurs in a time scale so that the molecules cannot complete their vibrational motions. Therefore the projected Marcus relation is not valid for cyclooctatetraene/cyclooctatetraene radical anion including a large torsional change during the electron transfer.     

Structures and energies of D-galactose and galabiose conformers as calculated by ab initio and semiempirical methods

Optimized geometries and total energies of some conformers of alpha- and beta-D-galactose have been calculated using the RHF/6-31G* ab initio method. Vibrational frequencies were computed at the 6-31G* level for the conformers that favor internal hydrogen bonding, in order to evaluate their enthalpies, entropies, Gibbs free energies, and then their structural stabilities. The semiempirical AM1, PM3, MNDO methods have also been performed on the conformers GG, GT, and TG of alpha- and beta-D-galactose. In order to test the reliability of each semiempirical method, the obtained structures and energies from the AM1, PM3, and MNDO methods have been compared to those achieved using the RHF/6-31G* ab initio method. The MNDO method has not been investigated further, because of the large deviation in the structural parameters compared with those obtained by the ab initio method for the galactose. The semiempirical method that has yielded the best results is AM1, and it has been chosen to perform structural and energy calculations on the galabiose molecule (the disaccharides constituted by two galactose units alpha 1,4 linked). The goal of such calculations is to draw the energy surface maps for this disaccharide. To realize each map, 144 different possible conformations resulting from the rotations of the two torsional angles psi and phi of the glycosidic linkage are considered. In each calculation, at each increment of psi and phi, using a step of 30 degrees from 0 to 330 degrees, the energy optimization is employed. In this article, we report also calculations concerning the galabiose molecule using different ab initio levels such as RHF/6-31G*, RHF/6-31G**, and B3Lyp/6-31G*.     

Probing furanose ring conformation by gas-phase computational methods: energy profile and structural parameters in methyl beta-D-arabinofuranoside as a function of ring conformation

The potential energy surface of methyl beta-D-arabinofuranoside (3) has been studied by ab initio molecular orbital (HF/6-31G) and density functional theory (B3LYP/6-31G) calculations via minimization of the 10 possible envelope conformers. The partial potential energy surface identified that the global minimum and lowest energy northern conformer was E(2). In the HF calculations, (2)E was the most stable southern conformer, while the density functional theory methods identified (4)E as the local minimum in this hemisphere. Additional calculations at higher levels of theory showed that the B3LYP-derived energies of many of the envelope conformers of 3 are dependent upon the basis set used. It has also been demonstrated that B3LYP/6-31+G//B3LYP/6-31G single point energies are essentially the same as those obtained from full geometry optimizations at the B3LYP/6-31+G level. The northern and southern minima of the B3LYP/6-31+G surface are, respectively, the E(2) and (2)E conformers. The B3LYP/6-31G geometries were used to study the relationship between ring conformation and various structural parameters including bond angles, dihedral angles, bond lengths, and interatomic distances.     

On the keto-enol tautomerization of malonaldehyde: an effective fragment potential study

Recent studies have mapped the keto-enol tautomerization of malonaldehyde through a general transition structure that leads exclusively to the Z isomer of the enol. However, it will be shown that a competing general transition structure exists that leads to both the E and Z isomers of the enol at the B3LYP/6-31G(d,p) and MP2/6-31G(d,p) levels of theory. Both the RHF- and DFT-based effective fragment potential methods have been used to model solvation effects, and the results are compared with full ab initio calculations. It is found that two bridging water molecules with two discrete DFT-based effective fragment potential solvent waters at the MP2/6-31G(d,p) level of ab initio theory provides the most computationally effective model for solvent effects in this system. It is shown that the relative energies for this QM/MM model differ from the full MP2/6-31G(d,p) energies by an average absolute relative difference of 2.2 kcal mol-1 across the reaction path when the zero-point vibrational energy correction is included.     

Structures of protonated arginine dimer and bradykinin investigated by density functional theory: further support for stable gas-phase salt bridges

The gas-phase structures and energetics of both protonated arginine dimer and protonated bradykinin were investigated using a combination of molecular mechanics with conformational searching to identify candidate low-energy structures, and density functional theory for subsequent minimization and energy calculations. For protonated arginine dimer, a good correlation (R = 0.88) was obtained between the molecular mechanics and EDF1 6-31+G* energies, indicating that mechanics with MMFF is suitable for finding low-energy conformers. For this ion, the salt-bridge or ion-zwitterion form was found to be 5.7 and 7.2 kcal/mol more stable than the simple protonated or ion-molecule form at the EDF1 6-31++G** and B3LYP 6-311++G** levels. For bradykinin, the correlation between the molecular mechanics and DFT energies was poor (R = 0.28), indicating that many low-energy structures are likely passed over in the mechanics conformational searching. This result suggests that structures of this larger peptide ion obtained using mechanics calculations alone are not necessarily reliable. The lowest energy structure of the salt-bridge form of bradykinin is 10.6 kcal/mol lower in energy (EDF1) than the lowest energy simple protonated form at the 6-311G* level. Similarly, the average energy of all salt-bridge structures investigated is 13.6 kcal/mol lower than the average of all the protonated forms investigated. To the extent that a sufficient number of structures are investigated, these results provide some additional support for the salt-bridge form of bradykinin in the gas phase.     

Theoretical determination of standard oxidation and reduction potentials of chlorophyll-a in acetonitrile

QM/MM calculations were performed on ethyl chlorophyllide-a and its radical cation and anion, by using the density functional (DF) B3LYP method to determine the molecular characteristics, and a molecular mechanics (MM) method to simulate the solvating medium. The presence of the solvent was accounted for during the optimization of the geometry of the 85-atom chlorophyll-a system by using an ONIOM methodology. A total of 24 solvent molecules were explicitly considered during the optimization process, and these were treated by the universal force field (UFF) method. Initially, the split-valence 3-21G basis set was used for optimizing the geometry of the 85-atom species, neutral, cation and anion. Electronic energies were then determined for the optimized species by making use of the polarized 6-31G(d) basis set. The ionization energy calculated (6.0 eV) is in very good agreement with the observed one (6.1 eV). The MM+ force field was used to investigate the dynamics of the acetonitrile molecules around the neutral species as well as the radical ions of chlorophyll. The required atomic charges on all the atoms were obtained from calculations on all involved molecules at the DFT/6-31G(d) level. Randomly sampled configurations were used to determine the first solvation layer contribution to the free energy of solvation of various species. A truncated 46-atom model of ethyl chlorophyllide-a was used to evaluate the thermal energies of neutral chlorophyll molecule relative to its two radical ions in the gas phase. Born energy, Onsager energy, and the Debye-Huckel energy of the chlorophyll-solvent aggregate were added as perturbative corrections to the free energy of solvation that was initially obtained through molecular dynamics method for the same complex. These calculations yield the oxidation potential as 0.75 +/- 0.32 V and the reduction potential -1.18 +/- 0.31 V at 298.15 K. The calculated values are in good agreement with the experimental midpoint potentials of +0.76 and -1.04 V, respectively.     

An ab initio study of intermolecular hydrogen bonding between small peptide fragments

Several small peptide fragments are investigated with ab initio (Hartree-Fock) calculations, using Gaussian basis sets. Complexation energies, net atomic charges, and optimum geometries are obtained. The geometries predicted by the STO -6G, and 6–31G* basis sets are quite similar, whereas the binding energies obtaiend by the 6–31G calculations are higher than those obtained with STO -6G and 6–31G* basis sets.