Elimination of translational and rotational motions in nuclear orbital plus molecular orbital theory

The nuclear orbital plus molecular orbital (NOMO) theory was developed in order to determine the nonadiabatic nuclear and electronic wave functions. This study presents a formulation to remove the contamination of rotational motion as well as translational motion in the NOMO theory. We have formulated the translation- and rotation-free (TRF)-NOMO theory by introducing the TRF Hamiltonian. The principal moment of inertia, which is the denominator in the rotational Hamiltonian, is expanded in a Taylor series. The zeroth-order of the Taylor expansion corresponds to a rigid-body rotator. The first-order terms contribute the coupling between the vibration and the rotation. Hartree-Fock equations have been derived in the framework of the TRF-NOMO theory. Numerical assessments, which were preformed for H2, D2, T2, mu2 (muon dimmer), and H2O, confirmed the importance of the TRF treatment.     

General orbital invariant MP2-F12 theory

A general form of orbital invariant explicitly correlated second-order closed-shell Moller-Plesset perturbation theory (MP2-F12) is derived, and compact working equations are presented. Many-electron integrals are avoided by resolution of the identity (RI) approximations using the complementary auxiliary basis set approach. A hierarchy of well defined levels of approximation is introduced, differing from the exact theory by the neglect of terms involving matrix elements over the Fock operator. The most accurate method is denoted as MP2-F12/3B. This assumes only that Fock matrix elements between occupied orbitals and orbitals outside the auxiliary basis set are negligible. For the chosen ansatz for the first-order wave function this is exact if the auxiliary basis is complete. In the next lower approximation it is assumed that the occupied orbital space is closed under action of the Fock operator [generalized Brillouin condition (GBC)]; this is equivalent to approximation 2B of Klopper and Samson [J. Chem. Phys. 116, 6397 (2002)]. Further approximations can be introduced by assuming the extended Brillouin condition (EBC) or by neglecting certain terms involving the exchange operator. A new approximation MP2-F12/3C, which is closely related to the MP2-R12/C method recently proposed by Kedzuch et al. [Int. J. Quantum Chem. 105, 929 (2005)] is described. In the limit of a complete RI basis this method is equivalent to MP2-F12/3B. The effect of the various approximations (GBC, EBC, and exchange) is tested by studying the convergence of the correlation energies with respect to the atomic orbital and auxiliary basis sets for 21 molecules. The accuracy of relative energies is demonstrated for 16 chemical reactions. Approximation 3C is found to perform equally well as the computationally more demanding approximation 3B. The reaction energies obtained with smaller basis sets are found to be most accurate if the orbital-variant diagonal Ansatz combined with localized orbitals is used for the first-order wave function. This unexpected result is attributed to geminal basis set superposition errors present in the formally more rigorous orbital invariant methods.     

Colle-Salvetti-type correction for electron-nucleus correlation in the nuclear orbital plus molecular orbital theory

A Colle-Salvetti (CS)-type electron-nucleus correction in the nuclear orbital plus molecular orbital theory is proposed. The CS-type correction is designed to satisfy the cusp condition for the electron-nucleus interaction. Since the CS-type correction is expressed in terms of the electron and nucleus densities, its evaluation is computationally feasible. Numerical assessment confirms that the CS-type correction performs well for the small G2 set.     

Molecular orbital studies of polyethylene deformation

Young's modulus E for polyethylene in the chain direction is calculated with molecular orbital theory applied to n-alkanes C₃H₈ through n-C₁₃H₂₈ and analyzed with the cluster-difference method. Semiempirical CNDO, MNDO, and AM1 models and ab initio HF/STO-3G, HF/6-31G, HF/6-31G*, and MP2/6-31G* models are used. Cluster-difference results, when extrapolated to infinite chain length, give E in good agreement with moduli evaluated with molecular cluster or crystal orbital methods, provided minimal basis sets are employed. E decreases from 495 GPa (CNDO) to 336 GPa (MP2/6-31G*) as the level of theory is improved, consistent with established behaviors of the various models. Our calculations do not reproduce earlier molecular cluster or crystal orbital results, which gave E < 330 GPa. The most rigorous MP2/6-31G* model is known to overestimate force constants by ∼ 11%; the scaled modulus E = 299 GPa is in good accord with E = 306 GPa from recent calculations based on experimental vibration frequencies. © 1996 John Wiley & Sons, Inc.     

Analytic second-order energy derivatives in natural orbital functional theory

The analytic energy gradients in the atomic orbital representation have recently been published (Mitxelena and Piris in J Chem Phys 146:014102, 2017) within the framework of the natural orbital functional theory (NOFT). We provide here an alternative expression for them in terms of natural orbitals, and use it to derive the analytic second-order energy derivatives with respect to nuclear displacements in the NOFT. The computational burden is shifted to the calculation of perturbed natural orbitals and occupancies, since a set of linear coupled-perturbed equations obtained from the variational Euler equations must be solved to attain the analytic Hessian at the perturbed geometry. The linear response of both natural orbitals and occupation numbers to nuclear geometry displacements need only specify the reconstruction of the second-order reduced density matrix in terms of occupation numbers.     

Analytical optimization of orbital exponents in Gaussian-type functions for molecular systems based on MCSCF and MP2 levels of fully variational molecular orbital method

We have analyzed the basis function series in molecular systems by optimization of orbital exponents in Gaussian-type functions (GTFs) including the electron correlation effects with multiconfiguration self-consistent field (MCSCF) and M?ller?CPlesset second-order perturbation (MP2) methods. First, we have derived and implemented the gradient formulas of MCSCF and MP2 energies with respect to GTF exponent, as well as GTF center and nuclear geometry, based on the fully variational molecular orbital (FVMO) method. Second, we have applied these electron-correlated FVMO methods to H₂, LiH, and hydrocarbon (CH₄, C₂H₆, C₂H₄, and C₂H₂) molecules. We have clearly demonstrated that the optimized exponent values with electron-correlated methods are different from those with simple Hartree?CFock method, since adequate basis functions for adequate virtual orbitals are indispensable to describe the accurate wave function and geometry for electron-correlated calculations.     

Møller-Plesset (MP2) perturbation theory for large molecules

Summary A novel formulation of MP2 theory is presented which starts from the Laplace transform MP2 ansatz, and subsequently moves from a molecular orbital (MO) representation to an atomic orbital (AO) representation. Consequently, the new formulation is denoted AO-MP2. As in traditional MP2 approaches electron repulsion integrals still need to be transformed. Strict bounds on the individual MP2 energy contribution of each intermediate four-index quantity allow to screen off numerically insignificant integrals with a single threshold parameter. Implicit in our formulation is a bound to two-particle density matrix elements. For small molecules the computational cost for AO-MP2 calculations is about a factor of 100 higher than for traditional MO-based approaches, but due to screening the computational effort in larger systems will only grow with the fourth power of the size of the system (or less) as is demonstrated both in theory and in application. MP2 calculations on (non-metallic) crystalline systems seem to be a feasible extension of the Laplace transform approach. In large molecules the AO-MP2 ansatz allows massively parallel MP2 calculations without input/output of four-index quantities provided that each processor has in-core memory for a limited number of two-index quantities. Energy gradient formulas for the AO-MP2 approach are derived.Dedicated to Prof. W. Kutzelnigg whose books on theoretical chemistry aroused my interest in this field     

Density functional theory treatment of electron correlation in the nuclear-electronic orbital approach

This paper presents the nuclear-electronic orbital density functional theory [NEO-DFT(ee)] method for including electron-electron correlation and nuclear quantum effects self-consistently in quantum chemical calculations. The NEO approach is designed to treat a relatively small number of nuclei quantum mechanically, while the remaining nuclei are treated classically. In the NEO-DFT(ee) approach, the correlated electron density is used to obtain the nuclear molecular orbitals, and the resulting nuclear density is used to obtain the correlated electron density during an iterative procedure that continues until convergence of both the nuclear and electronic densities. This approach includes feedback between the correlated electron density and the nuclear wavefunction. The application of this approach to bihalides and acetylene indicates that the nuclear quantum effects do not significantly impact the electron correlation energy, but the quantum nuclear energy is enhanced in the NEO-DFT(ee) B3LYP method. The excellent agreement of the NEO-DFT(ee)-optimized bihalide structures with the vibrationally averaged geometries from grid-based quantum dynamical methods provides validation for the NEO-DFT(ee) approach. Electron-proton correlation could be included by the development of an electron-nucleus correlation functional. Alternatively, explicit electron-proton correlation could be included directly into the NEO self-consistent-field framework with Gaussian-type geminal functions.     

Hybrid correlation models based on active-space partitioning: correcting second-order Moller-Plesset perturbation theory for bond-breaking reactions

Moller-Plesset second-order (MP2) perturbation theory breaks down at molecular geometries which are far away from equilibrium. We decompose the MP2 energy into contributions from different orbital subspaces and show that the divergent behavior of the MP2 energy comes from the excitations located within a small (or sometimes even the minimal) active space. The divergent behavior of the MP2 energy at large interfragment distances may be corrected by replacing a small number of terms by their more robust counterparts from coupled-cluster (CCSD) theory. We investigated several schemes of such a substitution, and we find that a coupling between the active-space CCSD and the remaining MP2 amplitudes is necessary to obtain the best results. This naturally leads us to an approach which has previously been examined in the context of cost-saving approximations to CCSD for equilibrium properties by Nooijen [J. Chem. Phys. 111, 10815 (1999)]. The hybrid MP2-CCSD approach, which has the same formal scaling as conventional MP2 theory, provides potential curves with a correct shape for bond-breaking reactions of BH, CH(4), and HF. The error of the MP2-CCSD method (measured against full configuration-interaction data) is smaller than that of MP2 at all interfragment separations and is qualitatively similar to that of full CCSD.     

Compared performances of the molecular orbital and density functional theories for fragments of biomolecules

As a preliminary to quantum-mechanical studies of nucleic acids, the possibilities of the molecular orbital (MO) method at the MP2 level and of the density functional theory (DFT) have been examined for typical fragments of nucleosides: methanol, tetrahydrofuran, 3-hydroxytetrahydrofuran and uracil. The results obtained from both methods using polarized split-valence basis sets are found to be consistent with each other, except for small but significant differences that are possibly related to the neglect of DFT self-interaction corrections.     

Gaussian-4 theory using reduced order perturbation theory

Two modifications of Gaussian-4 (G4) theory [L. A. Curtiss et al., J. Chem. Phys. 126, 084108 (2007)] are presented in which second- and third-order perturbation theories are used in place of fourth-order perturbation theory. These two new methods are referred to as G4(MP2) and G4(MP3), respectively. Both methods have been assessed on the G3/05 test set of accurate experimental data. The average absolute deviation from experiment for the 454 energies in this test set is 1.04 kcalmol for G4(MP2) theory and 1.03 kcalmol for G4(MP3) theory compared to 0.83 kcalmol for G4 theory. G4(MP2) is slightly more accurate for enthalpies of formation than G4(MP3) (0.99 versus 1.04 kcalmol), while G4(MP3) is more accurate for ionization potentials and electron affinities. Overall, the G4(MP2) method provides an accurate and economical method for thermochemical predictions. It has an overall accuracy for the G3/05 test set that is much better than G3(MP2) theory (1.04 versus 1.39 kcalmol) and even better than G3 theory (1.04 versus 1.13 kcalmol). In addition, G4(MP2) does better for challenging hypervalent systems such as H(2)SO(4) and for nonhydrogen species than G3(MP2) theory.     

Ab initio quantum chemical calculation of electron density,electrostatic potential,and electric field of biomolecule based on fragment molecular orbital method

Efficient quantum chemical calculations of electrostatic properties, namely, the electron density (EDN), electrostatic potential (ESP), and electric field (EFL), were performed using the fragment molecular orbital (FMO) method. The numerical errors associated with the FMO scheme were examined at the HF, MP2, and RI‐MP2 levels of theory using 4 small peptides. As a result, the FMO errors in the EDN, ESP, and EFL were significantly smaller than the magnitude of the electron correlation effects, which indicated that the FMO method provides sufficiently accurate values of electrostatic properties. In addition, an attempt to reduce the computational effort was proposed by combining the FMO scheme and a point charge approximation. The error due to this approximation was examined using 2 proteins, prion protein and human immunodeficiency virus type 1 protease. As illustrative examples, the ESP values at the molecular surface of these proteins were calculated at the MP2 level of theory.     

Improved correlation energy extrapolation schemes based on local pair natural orbital methods

It is well-known that the basis set limit is difficult to reach in correlated post Hartree-Fock ab initio calculations. One possible route forward is to employ basis set extrapolation schemes. In order to avoid prohibitively expensive calculations, the highest level calculation (typically based on the "gold standard" coupled cluster theory with single, double, and perturbative triple excitations, CCSD(T)) is only performed with the smallest basis set, and the remaining basis set incompleteness is estimated at a lower level of theory, typically second-order M?ller-Plesset perturbation theory (MP2). In this work, we provide a comprehensive investigation of alternative schemes where the MP2 extrapolation is replaced by the coupled-electron pair approximation, version 1 (CEPA/1) or the local pair natural orbital version of this method (LPNO-CEPA/1). It is shown that the MP2 method achieves apparent accuracy only due to error cancellation. Systematically more accurate results at small additional computational cost are obtained if the MP2 step is replaced by LPNO-CEPA/1. The errors of LPNO-CEPA/1 relative to canonical CEPA/1 are negligible. Owing to the highly systematic nature of the deviations between canonical and LPNO methods, basis set extrapolation reduces the LPNO errors in the total energies by 1 order of magnitude (~0.2 kcal/mol) and errors in energy differences to essentially zero. Using the CCSD(T)/LPNO-CEPA/1-based extrapolation scheme, new reference values are proposed for the recently published S66 set of interaction energies. The deviations between the new values and the original interactions energies are mostly very small but reach values up to 0.3 kcal/mol.     

The ellipsoidal gaussian basis in molecular orbital theory

New analytic integral formulas are presented for the potential energy integrals over ellipsoidal Gaussian basis functions [ exp (-x² - y² - z²)] that enter into solving the conventional expansion self-consistent field equations. Near minimal atomic orbital bases combined from large nuclear-centered primitive Gaussian sets are used in test calculations on the HF and CO molecules. The ellipsoidal exponential parameters for the valence atomic orbitals are fully optimized using a single scale factor for each atomic orbital and nuclear coordinate. The results are compared with those obtained using an unoptimized nuclear centered double-zeta spherical Gaussian basis.     

Semiempirical molecular orbital theory in carcinogenesis research

It is generally believed that most chemical carcinogens exert their effects through covalent modification of the nucleic acid base sites. Significant clarification of the underlying chemistry involved in these early and critical events is available through the techniques of computational chemistry. The role of semiempirical molecular orbital theory in this area is reviewed using examples focusing on

1. (i) the conversion of precarcinogens to reactive electrophiles,

2. (ii) the regiochemistry of nucleic acid adduct formation,

3. (iii) the chemical and physical consequences of nucleic acid base modification, and

4. (iv) empirical structure activity relationships.

The symbiotic relationship that can exist between the semiempirical and ab initio molecular orbital procedures is emphasized, as is the continuing niche for the semiempirical methodology in the study of systems that are too large, or in situations that are insufficiently cost effective, for the application of ab initio methods.     

Two-dimensional hückel molecular orbital theory

The conventional Hückel molecular orbital (HMO) theory is extended to include hydrocarbons with sp-hybridized C atoms in addition to the sp²-hybridization. MOs are constructed as linear combinations of 2p-orbitals in two dimensions perpendicular to each other. The resulting two parts of the π system (viz. π′ and π″) are assumed to be independent. The total π-bond order is assumed to be obtained by the addition of contributions from π′ and π″. The two-dimensional HMO theory covers hydrocarbons with acetylenic (CC) bonds and cumulated CC bonds. The theory is applied to vinylacetylene, diacetylene, allene, butatriene and divinylacetylene. CC bond lengths calculated from the theoretical π-bond orders are compared with experimental data. Agreements between calculated and observed bond lengths within ±0.03 Å are found.     

Density functional theory with fractional orbital occupations

In contrast to the original Kohn-Sham (KS) formalism, we propose a density functional theory (DFT) with fractional orbital occupations for the study of ground states of many-electron systems, wherein strong static correlation is shown to be described. Even at the simplest level represented by the local density approximation (LDA), our resulting DFT-LDA is shown to improve upon KS-LDA for multi-reference systems, such as dissociation of H(2) and N(2), and twisted ethylene, while performing similar to KS-LDA for single-reference systems, such as reaction energies and equilibrium geometries. Because of its computational efficiency (similar to KS-LDA), this DFT-LDA is applied to the study of the singlet-triplet energy gaps (ST gaps) of acenes, which are "challenging problems" for conventional electronic structure methods due to the presence of strong static correlation effects. Our calculated ST gaps are in good agreement with the existing experimental and high-level ab initio data. The ST gaps are shown to decrease monotonically with the increase of chain length, and become vanishingly small (within 0.1 kcal/mol) in the limit of an infinitely large polyacene. In addition, based on our calculated active orbital occupation numbers, the ground states for large acenes are shown to be polyradical singlets.     

Ab initio and density functional theory study of structures and energies for dimethoxymethane as a model for the anomeric effect

Ab initio molecular orbital theory and density functional theory calculations have been carried out on dimethoxymethane as a model for the anomeric effect. We optimized various conformations of dimethoxymethane using Gaussian 92 at the MP2/6-311 + + G**, MP2/DZP + Diffuse, MP2/6-31G**, and Becke3LYP/6-31G** levels of theory. These methods were evaluated based on their performance in reproducing structures and energies of dimethoxymethane when compared to experiment. This study also examined the structure and energy of dimethoxymethane as a function of dihedral angles for examining the anomeric effect at the MP2/6-31G** and Becke3LYP/6-31G** levels of theory. These calculations are qualitatively consistent with the anomeric effect observations in carbohydrates and with earlier calculations. Quantitative comparisons with earlier results reveal that dimethoxymethane has lower total energies, smaller rotational barriers, and shorter bond lengths than was previously determined. The Becke3LYP calculations were also compared to the MP2 results. The density functional theory findings show that the minimum energy structures correspond well with experimental and MP2 data. The total and relative energies from molecular orbital theory and density functional theory vary to some extent. Contour plots of the relative energies of dimethoxymethane were evaluated and compared to a relative energy contour plot determined by MM3. The contour plots were similar, showing slightly larger changes in energies for the MP2 results than for the Becke3LYP results, which in turn were slightly larger than the MM3 results. Density functional theory calculations are an excellent alternative method of calculation due to increased speed and reliable accuracy of the density functional calculations. These results will serve as a benchmark for modelling the anomeric effect in carbohydrates. © 1996 John Wiley & Sons, Inc.     

Decomposition of nuclear magnetic resonance spin-spin coupling constants into active and passive orbital contributions

The theory of the J-OC-PSP (decomposition of J into orbital contributions using orbital currents and partial spin polarization) method is derived to distinguish between the role of active, passive, and frozen orbitals on the nuclear magnetic resonance (NMR) spin-spin coupling mechanism. Application of J-OC-PSP to the NMR spin-spin coupling constants of ethylene, which are calculated using coupled perturbed density functional theory in connection with the B3LYP hybrid functional and a [7s,6p,2d/4s,2p] basis set, reveal that the well-known pi mechanism for Fermi contact (FC) spin coupling is based on passive pi orbital contributions. The pi orbitals contribute to the spin polarization of the sigma orbitals at the coupling nuclei by mediating spin information between sigma orbitals (spin-transport mechanism) or by increasing the spin information of a sigma orbital by an echo effect. The calculated FC(pi) value of the SSCC (1)J(CC) of ethylene is 4.5 Hz and by this clearly smaller than previously assumed.