Toward more efficient CCSD(T) calculations of intermolecular interactions in model hydrogen-bonded and stacked dimers

Interaction energies of the model H-bonded complexes, the formamide and formamidine dimers, as well as the stacked formaldehyde and ethylene dimers are calculated by the coupled cluster CCSD(T) method. These systems serve as a model for H-bonded and stacking interactions, typical in molecules participating in biological systems. We use the optimized virtual orbital space (OVOS) technique, by which the dimension of the space of virtual orbitals in coupled cluster CCSD(T) calculations can be significantly reduced. We demonstrate that when the space of virtual orbitals is reduced to 50% of the full space, which means reducing computational demands by 1 order of magnitude, the interaction energies for both H-bonded and stacked dimers are affected by no more than 0.1 kcal/mol. This error is much smaller than the error when interaction energies are calculated using limited basis sets.     

Are atoms intrinsic to molecular electronic wavefunctions? I. The FORS model

The model of the full optimized reaction space describes the electronic structure of a molecule in terms of the best wave-function that can be obtained as a superposition of all those configurations which are generate possible occupancies and couplings from a “formal minimal basis” of valence, orbitals on the constituent atoms. These configurations span a “full reaction space”, and MC SCF optimization of the orbitals in terms of an extended set of quantitative basis orbitals determines the full optimized reaction space (FORS). Basic justifications, methodological specifics and sample applications are discussed.     

Quantum molecular mechanics—a noniterative procedure for the fast Ab Initio calculation of closed shell systems

In this article, we advance the foundations of a strategy to develop a molecular mechanics method based not on classical mechanics and force fields but entirely on quantum mechanics and localized electron‐pair orbitals, which we call quantum molecular mechanics (QMM). Accordingly, we introduce a new manner of calculating Hartree–Fock ab initio wavefunctions of closed shell systems based on variationally preoptimized nonorthogonal electron pair orbitals constructed by linear combinations of basis functions centered on the atoms. QMM is noniterative and requires only one extremely fast inversion of a single sparse matrix to arrive to the one‐particle density matrix, to the electron density, and consequently, to the ab initio electrostatic potential around the molecular system, or cluster of molecules. Although QMM neglects the smaller polarization effects due to intermolecular interactions, it fully takes into consideration polarization effects due to the much stronger intramolecular geometry distortions. For the case of methane, we show that QMM was able to reproduce satisfactorily the energetics and polarization effects of all distortions of the molecule along the nine normal modes of vibration, well beyond the harmonic region. We present the first practical applications of the QMM method by examining, in detail, the cases of clusters of helium atoms, hydrogen molecules, methane molecules, as well as one molecule of HeH⁺ surrounded by several methane molecules. We finally advance and discuss the potentialities of an exact formula to compute the QMM total energy, in which only two center integrals are involved, provided that the fully optimized electron‐pair orbitals are known. © 2012 Wiley Periodicals, Inc.     

Are atoms sic to molecular electronic wavefunctions? II. Analysis of fors orbitals

Electronic wavefunctions that describe molecules in the full optimized reaction space (FORS) are multiconfigurational wavefunctions which are invariant under non-singular linear transformations of the occupied molecular orbitals. They offer therefore a considerably wider scope for orbital interpretations than the single-configuration Hartree-Fock approximation. For example they can be analyzed in terms of natural MOs and in terms of localized MOs. The latter turn out to be remarkably atomic in character and a new localization procedure can be formulated which yields atom-adapted molecular orbitals. These have the character of minimal-basis-set AOs that are optimally adapted to the molecular environment and furnish an unambigious atomic population analysis. On the other hand, chemically adapted molecular orbitals can be defined by an appropriate compromise between natural orbitals and localized orbitals. The freedom to use, as configuration-generating molecular orbitals, atom-adapted FORS MOs as well as chemically adapted FORS MOs makes FORS wavefunctions particularly suitable for chemical interpretations. The ensuing analysis establishes the minimal basis set (in molecule-adapted form) as a theoretically sound concept for the understanding of accurate molecular wavefunctions. An illustrative example is discussed.     

Orbital-optimized opposite-spin scaled second-order correlation: an economical method to improve the description of open-shell molecules

Coupled-cluster methods based on Brueckner orbitals are well known to resolve the problems of symmetry breaking and spin contamination that are often associated with Hartree-Fock orbitals. However, their computational cost is large enough to prevent application to large molecules. Here the authors present a simple approximation where the orbitals are optimized with the mean-field energy plus a correlation energy taken as the opposite-spin component of the second-order many-body correlation energy, scaled by an empirically chosen parameter (recommended as 1.2 for general applications). This "optimized second-order opposite-spin" (abbreviated as O2) method requires fourth-order computation on each orbital iteration. O2 is shown to yield predictions of structure and frequencies for closed-shell molecules that are very similar to scaled second-order Moller-Plesset methods. However, it yields substantial improvements for open-shell molecules, where problems with spin contamination and symmetry breaking are shown to be greatly reduced.     

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.     

Using the local elevation method to construct optimized umbrella sampling potentials: Calculation of the relative free energies and interconversion barriers of glucopyranose ring conformers in water

A method is proposed to combine the local elevation (LE) conformational searching and the umbrella sampling (US) conformational sampling approaches into a single local elevation umbrella sampling (LEUS) scheme for (explicit‐solvent) molecular dynamics (MD) simulations. In this approach, an initial (relatively short) LE build‐up (searching) phase is used to construct an optimized biasing potential within a subspace of conformationally relevant degrees of freedom, that is then used in a (comparatively longer) US sampling phase. This scheme dramatically enhances (in comparison with plain MD) the sampling power of MD simulations, taking advantage of the fact that the preoptimized biasing potential represents a reasonable approximation to the negative of the free energy surface in the considered conformational subspace. The method is applied to the calculation of the relative free energies of β‐D ‐glucopyranose ring conformers in water (within the GROMOS 45A4 force field). Different schemes to assign sampled conformational regions to distinct states are also compared. This approach, which bears some analogies with adaptive umbrella sampling and metadynamics (but within a very distinct implementation), is shown to be: (i) efficient (nearly all the computational effort is invested in the actual sampling phase rather than in searching and equilibration); (ii) robust (the method is only weakly sensitive to the details of the build‐up protocol, even for relatively short build‐up times); (iii) versatile (a LEUS biasing potential database could easily be preoptimized for small molecules and assembled on a fragment basis for larger ones). © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Calculation of transition matrix elements by nonsingular orbital transformations

A general strategy is described for the evaluation of transition matrix elements between pairs of full class CI wave functions built up from mutually nonorthogonal molecular orbitals. A new method is proposed for the counter‐transformation of the linear expansion coefficients of a full CI wave function under a nonsingular transformation of the molecular‐orbital basis. The method, which consists in a straightforward application of the Cauchy–Binet formula to the definition of a Slater determinant, is shown to be simple and suitable for efficient implementation on current high‐performance computers. The new method appears mainly beneficial to the calculation of miscellaneous transition matrix elements among individually optimized CASSCF states and to the re‐evaluation of the CASCI expansion coefficients in Slater‐determinant bases formed from arbitrarily rotated (e.g., localized or, conversely, delocalized) active molecular orbitals. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

On the stability of Be3: a benchmark complete active space self-consistent field + averaged quadratic coupled cluster study

The optimized geometries and binding energies for the linear and triangular isomers of the beryllium trimer have been obtained through benchmark multireference averaged quadratic coupled cluster (AQCC) calculations using very large complete active space SCF (CASSCF) references (12 active electrons in 13 and 14 orbitals). Geometries were optimized with the cc-pV5Z basis, while the binding energies (including counterpoise correction) were obtained with the significantly larger aug-cc-pV5Z basis set. The binding energies (27.3 and 16.3 kcal/mol for the equilateral and linear isomers, respectively) are larger than the previous full CI benchmark values, while the corresponding Be-Be equilibrium distances of 4.101 and 4.088 a.u. are smaller. In view of the near-size consistency character of the CASSCF + AQCC method, the fact that all 12 electrons are fully correlated, the active reference space includes 14 orbitals, and the very large basis set used here, we propose to consider these results as reference data for Be(3). Using the electron pair localization function obtained at the CASSCF(12,15) level, it is clearly illustrated that the 2p orbitals lying in the molecular plane play a dominant role in the bonding pattern for the equilateral isomer.     

Symmetry and equivalence restrictions in electronic structure calculations

A simple method for obtaining MCSCF orbitals and CI natural orbitals adapted to degenerate point groups, with full symmetry and equivalence restrictions, is described. Among several advantages accruing from this method are the ability to perform atomic SCF calculations on states for which the SCF energy expression cannot be written in terms of Coulomb and exchange integrals over real orbitals, and the generation of symmetry-adapted atomic natural orbitals for use in a recently proposed method for basis set contraction.     

Parallel, linear-scaling building-block and embedding method based on localized orbitals and orbital-specific basis sets

We present a linear scaling method for the energy minimization step of semiempirical and first-principles Hartree-Fock and Kohn-Sham calculations. It is based on the self-consistent calculation of the optimum localized orbitals of any localization method of choice and on the use of orbital-specific basis sets. The full set of localized orbitals of a large molecule is seen as an orbital mosaic where each tessera is made of only a few of them. The orbital tesserae are computed out of a set of embedded cluster pseudoeigenvalue coupled equations which are solved in a building-block self-consistent fashion. In each iteration, the embedded cluster equations are solved independently of each other and, as a result, the method is parallel at a high level of the calculation. In addition to full system calculations, the method enables to perform simpler, much less demanding embedded cluster calculations, where only a fraction of the localized molecular orbitals are variational while the rest is frozen, taking advantage of the transferability of the localized orbitals of a given localization method between similar molecules. Monitoring single point energy calculations of large poly(ethylene oxide) molecules and three dimensional carbon monoxide clusters using an extended Huckel Hamiltonian are presented.     

Optimization of auxiliary basis sets for the LEDO expansion and a projection technique for LEDO-DFT

We present a systematic procedure for the optimization of the expansion basis for the limited expansion of diatomic overlap density functional theory (LEDO-DFT) and report on optimized auxiliary orbitals for the Ahlrichs split valence plus polarization basis set (SVP) for the elements H, Li--F, and Na--Cl. A new method to deal with near-linear dependences in the LEDO expansion basis is introduced, which greatly reduces the computational effort of LEDO-DFT calculations. Numerical results for a test set of small molecules demonstrate the accuracy of electronic energies, structural parameters, dipole moments, and harmonic frequencies. For larger molecular systems the numerical errors introduced by the LEDO approximation can lead to an uncontrollable behavior of the self-consistent field (SCF) process. A projection technique suggested by L?wdin is presented in the framework of LEDO-DFT, which guarantees for SCF convergence. Numerical results on some critical test molecules suggest the general applicability of the auxiliary orbitals presented in combination with this projection technique. Timing results indicate that LEDO-DFT is competitive with conventional density fitting methods.     

Coupled-cluster with active space selected higher amplitudes: performance of seminatural orbitals for ground and excited state calculations

The active space approach for coupled-cluster models is generalized using the general active space concept and implemented in a string-based general coupled-cluster code. Particular attention is devoted to the choice of orbitals on which the subspace division is based. Seminatural orbitals are proposed for that purpose. These orbitals are obtained by diagonalizing only the hole-hole and particle-particle block of the one-electron density of a lower-order method. The seminatural orbitals are shown to be a good replacement for complete active space self-consistent field orbitals and avoid the ambiguities with respect to the reference determinant introduced by the latter orbitals. The seminatural orbitals also perform well in excited state calculations, including excited states with strong double excitation contributions, which usually are difficult to describe with standard coupled-cluster methods. A set of vertical excitation energies is obtained and benchmarked against full configuration interaction calculations, and alternative hierarchies of active space coupled-cluster models are proposed. As a simple application the spectroscopic constants of the C(2) B (1)Delta(g) and B(') (1)Sigma(g) (+) states are calculated using active space coupled-cluster methods and basis sets up to quadruple-zeta quality in connection with extrapolation and additivity schemes.     

Benchmarking NMR indirect nuclear spin-spin coupling constants: SOPPA, SOPPA(CC2), and SOPPA(CCSD) versus CCSD

Accurate calculations of NMR indirect nuclear spin-spin coupling constants require especially optimized basis sets and correlated wave function methods such as CCSD or SOPPA(CCSD). Both methods scale as N(6), where N is the number of orbitals, which prevents routine applications to molecules with more than 10-15 nonhydrogen atoms. We have therefore developed a modification of the SOPPA(CCSD) method in which the CCSD singles and doubles amplitudes are replaced by CC2 singles and doubles amplitudes. This new method, called SOPPA(CC2), scales only as N(5), like the original SOPPA-method. The performance of the SOPPA(CC2) method for the calculation of indirect nuclear spin-spin coupling constants is compared to SOPPA and SOPPA(CCSD) employing a set of benchmark molecules. We also investigate the basis set dependence by employing three different basis sets optimized for spin-spin coupling constants, namely the HuzIV-su4, ccJ-pVTZ, and ccJ-pVQZ basis sets. The results of the corresponding CCSD calculations are used as a theoretical reference.     

Atomic orbitals in molecules: general electronegativity and improvement of Mulliken population analysis

An approach of atomic orbitals in molecules (AOIM) has been developed to study the atomic properties in molecules, in which the molecular orbitals are expressed in terms of the optimized minimal atomic orbitals. The atomic electronegativities are calculated using Pauling's electronegativity of free atom and are employed to find the electronegativity equilibrium in molecules and to describe the amphoteric properties of the transition metals from the groups 4 to 10. AOIM can also improve the numerical stability and accuracy of the original Mulliken population analysis.     

Optimized linear combinations of simple exponentials for atomic systems

Hartree-Fock wave functions for the He and Be isoelectronic sequences of ions are calculated using orbitals which are linear combinations of simple exponential functions. By a full optimization of the exponents and coefficients close approximations to the HartreeFock energies were obtained. To the same order of accuracy the resulting Hartree–Fock orbitals require fewer basis functions than used previously. A number of difficulties which arise in the numerical procedures as the size of the basis set is increased are analysed in detail. Similar results are obtained for the Li sequence using the Unrestricted HartreeFock method with and without projection.     

A hybrid ab initio/free electron computational model for conjugated dye molecules: Simple cyanines and oxonols

Justifications developed for the application the free electron model to the π‐orbitals of conjugated molecules suggest that the optical properties of these molecules would be well described by a one‐dimensional free electron model with a potential chosen to reproduce the energy level spacing of the ground state occupied π‐orbitals. Such a hybrid ab initio/free electron modeling approach, where the free electron potential parameters are optimized on a molecule‐by‐molecule basis, is developed, and applied to a series of simple cyanine and oxonol dyes. The ensuing predictions for λ_max, oscillator strengths, and redox properties compare well to available experimental information. Two important strengths of this approach are that no explicit calculations of the excited electronic state are required, and that the ab initio determination of the occupied π‐orbital level spacing considers all the electrons (π and σ) of the entire molecule in a specified geometry, environment, etc. This second characteristic gives the ability to efficiently model modifications of the optical properties of conjugated molecules resulting from chemical and/or physical modifications occuring within and remote to the conjugated region of the molecule. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 943–953, 2000     

Superior performance of Mukherjee’s state-specific multi-reference coupled-cluster theory at the singles and doubles truncation scheme with localized active orbitals

The state-specific multi-reference coupled-cluster (SS-MRCC) theory of Mukherjee et al., in its singles and doubles truncation scheme (SS-MRCCSD), misses important couplings between the virtual functions reached by single and double excitations from different model functions. Since the SS-MRCC theory is not invariant with respect to the transformations among the active orbitals, the results are dependent on the active orbitals chosen. We demonstrate in this paper with results for potential energy curves for several example molecules involving single and multiple bond dissociation that the performance of SS-MRCCSD is significantly improved if localized active orbitals are used. The improvement is remarkable both in terms of the non-parallelity error and the magnitude of correlation energy recovered vis-a-vis the full configuration interaction results with the same basis set. The results bolster our claim that SS-MRCCSD with localized orbitals is an accurate general theory for potential energy surfaces.     

Complete active-space configuration interaction with optimized orbitals: Application to Li2

We performed CAS –CI calculations on Li₂ using a set of molecular orbitals (MO ) optimized with a procedure that, in the case of highly symmetric molecules, permits extraction of a small set of MO out of a large set of atomic orbitals (AO ). The dimension of the CAS –CI space was of about 12 million symmetry-adapted determinants. We determined some spectroscopic constants of Li₂ with three different atomic basis sets of increasing quality. The values obtained with the largest atomic basis set are very close to the experimental results. © 1995 John Wiley & Sons, Inc.     

CASSCF optimization problem for a group of excited states

A second-order version of the CASSCF approach to the optimization problem for a single (ground or excited) state and a group of excited states (involving, if necessary, also the ground state) is proposed. In contrast to the already existing methods, in the frameworks of our approach, there arises no need in completing the set of states to be optimized to the full basis set of configuration function space. Generation of secondary orbitals in the course of orbital optimization is also not required. All necessary integral transformations are performed only with active orbitals. These certainly attractive features of our approach appear due to employing the Gauss parametrization of average electronic energy domain, which is nonstandard in quantum chemistry. © 1992 John Wiley & Sons, Inc.