A comparative analysis of Hartree-Fock and Kohn-Sham orbital energies

Hartree-Fock and Kohn-Sham orbital energies, the latter computed with several different exchange/correlation functionals, are compared and analyzed for 12 molecules. The Kohn-Sham energies differ significantly from experimental ionization energies, but by amounts that are, for a given molecule and exchange/correlation functional, roughly the same for all of the valence orbitals. With the exchange/correlation functionals used, the energy of the highest occupied Kohn-Sham orbital does not approximate the corresponding ionization potential any better than do the other orbital energies. Received: 24 October 1997 / Accepted 31 October 1997     

DFT‐Aided Vibrational Circular Dichroism Spectroscopy Study of (−)‐S‐cotinine

The implementation of a strategy comprising the use of vibrational circular dichroism spectroscopy and DFT calculations allows determination of the solution‐state conformational distribution in (?)‐S‐cotinine, giving further proof of the extra conformer‐discriminating potential of this experimental technique, which may offer unique molecular fingerprints of subtly dissimilar molecular conformers of chiral samples. Natural bond orbital electronic structure calculations of the rotational barrier height between the two main conformers of the species indicate that hyperconjugative effects are the key force governing the conformational equilibrium. The negligible effect of the solvent’s polarity over both structure and conformational energy profile supports this result.     

An ab initio study of heterodienophiles addition to 2,3-diaza-1,3-butadiene: An example of endo-lone-pair effect on the reaction energy barrier

Transition states for the Diels-Alder reactions of 2,3-diaza-1,3-butadiene with ethylene, formaldehyde, formaldimine, cis- and trans- diazene, and nitrosyl hydride were located by ab initio molecular orbital calculations. The bond orders of the new forming bonds have been used to determine the asynchronicity of the reactions. Ab initio calculations show that the energy barrier for the hetero-Diels-Alder reactions is relatively high. The highest energy barrier of 34.76 kcal/mol calculated at the MP4/6-31G*//MP2/6-31G* level was found for the exo-cis-diazene addition to 2,3-diaza-1,3-butadiene. In all cases, when two diastereomeric transition structures are possible, the one with the endo hydrogen, exo lone pair was predicted to have a lower activation barrier. This behavior can be explained by the n-π and n-n loan pair repulsion interaction between the dienophile and diene heteroatoms in the corresponding transition state. The barrier is higher for those reactions which in the transition state have more lone electron pairs. Also, the barrier is higher when the lone pairs are endo oriented than when they are exo oriented in the transition state. © 1996 by John Wiley & Sons, Inc.     

(Z)-1,2-二苯基-1,3-丁二烯与单取代的乙烯衍生物区位选择性环加成反应

(Z)-1,2-二苯基-1,3-丁二烯分别与甲基乙烯基酮、丙烯酰胺、丙烯腈以及丙烯醛发生环加成反应,均生成顺式和反式的连位取代环己烯衍生物。以IR、¹HNMR和MS鉴定了产物的结构。利用CNDO/2法计算了分子轨道能量和系数,由前线分子轨道理论解释了这种环加成反应的区位选择性。     

Theoretical studies of cyclopropylsilylenes: The structures and stability of cyclopropylsilylene C3H5SiH

Ab initio molecular orbital calculations at the G2(MP2) level have been carried out on cyclopropylsilylene C₃H₅SiH. Four equilibrium structures were located. Like H₂Si, the ground state of C₃H₅SiH is singlet and the triplet is the low‐lying excited state. The singlet–triplet separation energy is 127.9 kJ/mol. The cis‐trans isomerization path of singlet cyclopropylsilylene was investigated by intrinsic reaction coordinate (IRC) calculations. The calculations show that no gauche conformers exist along the potential energy curve of the cis‐trans isomerization and the isomerization happens with a barrier of 30.1 kJ/mol. Changes (ΔH and ΔG) in thermodynamic functions, equilibrium constant K(T), and A factor and reaction rate constant k(T) in Eyring transition state theory of the cis‐trans isomerization were also calculated. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Π delocalisation in 1,3-butadiene: a nonempirical molecular orbital study

A value of −0.33 eV or −7.6 kcal mol⁻¹ has been obtained for the vertical delocalisation energy of trans-1,3-butadiene from a nonempirical molecular orbital calculation on the π system. The result agrees well enough with ab initio calculations to suggest that a simplified approach need not be semiempirical. In a basis of orthogonalised atomic orbitals the central bond order is found to be 0.295 (Hückel value 0.447) for the delocalised structure and 0.125 for the localised (Hückel value zero). Core resonance integrals between neighbouring atoms, the analogues of Hückel's β, have theoretical values of −3.9 and −3.2 eV compared with −3.6 eV in benzene. Received: 11 May 1999 / Accepted: 22 July 1999 / Published online: 2 November 1999     

Ab initio Mechanism Study on the Reaction of Chlorine Atom with Formic Acid

The potential energy surface(PES) for the reaction of Cl atom with HCOOH is predicted using ab initio molecular orbital calculation methods at UQCIDS(T,full)6-311 G(3df,2p)//UMP2(full)/6-311 G(d,P) level of theory with zero-point vibrational energy (ZPVE) correction.The calculated results show that the reaction mechanism of Cl atom with formic acid is a C-site hydrogen abstraction reaction from cis-HOC(H)O molecule by Cl atom with a 3.73kJ/mol reaction barrier height,leading to the formation of cis-HOCO radical which will reacts with Cl atom or other molecules in such a reaction system.Because the reaction barrier height of O-site hydrogen abstraction reaction from cis-HOC(H)O molecule by Cl atom which leads to the formation of HCO2 radical is 67.95kJ/mol,it is a secondary reaction channel in experiment,This is in good agreement with the prediction based on the previous experiments.     

General method for symmetry orbitals and tensors in electronic structure calculations

The symmetry orbital tensor (SOT) method, which makes full use of symmetries in all point groups and can be applied to the self-consistent field (SCF) and post-SCF calculations, is introduced. The principal feature of this method is the definition of the symmetry orbitals (SOs). Any element in a molecular point group will transform one SO to another equivalent SO or simply to itself, and no mixture among SOs exists. Thus, although the SOs for non-Abelian point groups may adapt to reducible representations, their transformation properties are much simpler than in conventional treatments. This article also presents a general scheme to generate SOs for all point groups. The direct products of N SOs form an Nth-rank SOT group, and each matrix element between SOTs is the product of a physical factor and a geometric factor. Compared with the canonical molecular orbitals, the use of SOs can noticeably reduce the computation efforts by decreasing the number of integrals needed in the SCF calculations or the number of configurations needed in the configuration interaction (CI) calculations. The SOT-SCF and SOT-CI approaches are formulated and a preliminary SOT-SCF program is written. Pilot calculations demonstrate the value of the SOT approach, at least at the closed-shell Hartree–Fock level. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 305–321, 1999     

Molecular states of atoms. I. Orbital energy parameters

Atomic valence state energies are analyzed to obtain values of orbital energy parameters that may be used in semiempirical molecular orbital calculations. Difficulty in defining the interaction between orbitals with non-integer electron populations is systematically avoided by distinguishing between a valence state and a molecular state of an atom, only the latter state having non-integer spin paired orbital occupancy. Application of the virial theorem to the molecular state enables a value for the orbital kinetic energy to be obtained from the valence state orbital energy parameters once an arbitrary configuration is defined as reference. The orbitals then are eigenfunctions of the atomic Fock operator for that reference molecular state and, with their energy parameters, may be employed as a fixed basis set for molecular orbital calculations.     

Ab initio through‐space/bond interaction analysis of the long C–C bonds in Bi(anthracene‐9,10‐dimethylene) photoisomers

The bi(anthracene‐9,10‐dimethylene) photoisomer has remarkably long C–C single bonds. To examine the lengthening of the C–C bond, we propose a novel procedure for quantitatively analyzing orbital interactions in a molecule at the level of the ab initio molecular orbital method. In this procedure, we can cut off the specific through‐space/bond interactions in a molecule by artificially increasing the absolute magnitude of the exponents in a Gaussian function. Then, the spatial orbital interactions were perfectly cut off, and, each term that makes up the total energy, that is, the nuclear–electron attractions, the electron–electron repulsions, and the nuclear–nuclear repulsions cancel each other. Several model molecules of the photoisomer were analyzed by this procedure. It was found that the orbital interaction between the p orbital on the benzene ring and the σ* orbital on the C–C bond in question, σ→σ* electron transfer through π orbital, weakens the C–C bond efficiently when these orbitals were located in the “periplanar” conformation. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Density-functional theory with effective potential expressed as a direct mapping of the external potential: applications to atomization energies and ionization potentials

In this paper the authors further develop and apply the direct-mapping density functional theory to calculations of the atomization energies and ionization potentials. Single-particle orbitals are determined by solving the Kohn-Sham [Phys. Rev. A. 140, 1133 (1965)] equations with a local effective potential expressed in terms of the external potential. A two-parametric form of the effective potential for molecules is proposed and equations for optimization of the parameters are derived using the exchange-only approximation. Orbital-dependent correlation functional is derived from the second-order perturbation theory in its Moller-Plesset-type zeroth-order approximation based on the Kohn-Sham orbitals and orbital energies. The total atomization energies and ionization potentials computed with the second-order perturbation theory were found to be in agreement with experimental values and benchmark results obtained with ab initio wave mechanics methods.     

How Ionization Catalyzes Diels-Alder Reactions

The catalytic effect of ionization on the Diels-Alder reaction between 1,3-butadiene and acrylaldehyde has been studied using relativistic density functional theory (DFT). Removal of an electron from the dienophile, acrylaldehyde, significantly accelerates the Diels-Alder reaction and shifts the reaction mechanism from concerted asynchronous for the neutral Diels-Alder reaction to stepwise for the radical-cation Diels-Alder reaction. Our detailed activation strain and Kohn-Sham molecular orbital analyses reveal how ionization of the dienophile enhances the Diels-Alder reactivity via two mechanisms: (i) by amplifying the asymmetry in the dienophile's occupied π-orbitals to such an extent that the reaction goes from concerted asynchronous to stepwise and thus with substantially less steric (Pauli) repulsion per reaction step; (ii) by enhancing the stabilizing orbital interactions that result from the ability of the singly occupied molecular orbital of the radical-cation dienophile to engage in an additional three-electron bonding interaction with the highest occupied molecular orbital of the diene.     

Indo-Mo Model with s·p Separation

A newly developed self-consistent-field molecular orbital theory is described and tested. The intermediate neglect of differential overlap (INDO) approximation is used, and all interaction integrals are differentiated according to their dependence upon 2s and 2p AO's. The bonding parameter β_l is reformulated so that the model is calibrated to only one specific molecular property, namely, the ionization energy of Hi. We expect that this model will not be biased toward any special families of molecules.     

Model study of the impact of orbital choice on the accuracy of coupled-cluster energies. II. Valence-universal coupled-cluster method

The impact of the choice of molecular orbital sets on the results of the valence-universal coupled cluster method involving up to three-body amplitudes (VU-CCSDT) was studied for the H4 model. This model offers a straightforward way of representing all possible symmetry-adapted orbitals. Moreover, the degree of quasi-degeneracy of its lowest ¹A₁ states can be varied over a wide range by changing its geometry. Calculations were performed both for 13 sets of standard quantum chemical orbitals and for a vast variety of nonstandard orbital sets defined by nodes of a two-dimensional orbital grid. The performance of various standard orbital sets in VU-CCSDT calculations is compared. It is also documented that for every quasi-degeneracy region there exist nonstandard orbital sets which allow one to obtain more accurate VU-CCSDT energies than the standard orbital sets. In an attempt to provide a general interpretation for some of the alternative orbital sets, we defined a set of orbitals which maximize the proximity of the model and target spaces—maximum proximity orbitals (MPO). It is demonstrated that outside the strong quasi-degeneracy region the energies obtained for the VU-CCSDT approach based on the MPOs are more accurate than for the standard orbital sets. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 67: 221–237, 1998     

Criteria for pericyclic and pseudopericyclic character of electrocyclization of (Z)-1,2,4,6-heptatetraene and (2Z)-2,4,5-hexatriene-1-imine

The electrocyclic reaction mechanisms of (Z)-1,2,4,6-heptatetraene and (2Z)-2,4,5-hexatriene-1-imine were studied by ab initio MO methods. The activation energy barrier height of the electrocyclic reaction of (Z)-1,2,4,6-heptatetraene is extremely a low energy barrier of 8.58 kcal/mol by a MRMP method. The activation energy barrier height of the electrocyclic ring closure of the trans-type of (2Z)-2,4,5-hexatriene-1-imine is lower by 3.18 kcal/mol than that of (Z)-1,2,4,6-heptatetraene. These low energy barriers come from some orbital interactions relating to allene group. For the reaction of (Z)-1,2,4,6-heptatetraene, the interactions of the vertical and side π orbitals of the allene group with another terminal π orbital are important at the transition state. The interaction of the vertical π orbital of allene group with a lone pair orbital of N atom is dominant at the transition state of the reaction of the trans-type of (2Z)-2,4,5- hexatriene-1-imine. The electrocyclic mechanism of the cis-type of (2Z)-2,4,5-hexatriene-1-imine was also discussed. Contribution of the Mark S. Gordon 65th Birthday Festschrift issue.     

Quantum chemistry beyond Born–Oppenheimer approximation on a quantum computer: A simulated phase estimation study

We present an efficient quantum algorithm for beyond‐Born–Oppenheimer molecular energy computations. Our approach combines the quantum full configuration interaction method with the nuclear orbital plus molecular orbital method. We give the details of the algorithm and demonstrate its performance by classical simulations. Two isotopomers of the hydrogen molecule (H₂, HT) were chosen as representative examples and calculations of the lowest rotationless vibrational transition energies were simulated. © 2016 Wiley Periodicals, Inc.     

On the use of isovalued surfaces to determine molecule shape and reaction pathways

Summary Novel insights into local molecule structure and reactivity can be gained from viewing isovalued surfaces of the molecular electron density, electrostatic potential and molecular orbitals rendered as colored, 3-D objects. For example, drawing positive and negative electrostatic isopotential surfaces partitions the molecule into regions subject to nucleophilic or electrophilic attack. Similarly, coloring isodensity surfaces to indicate the magnitude of the gradient of the electron density maps the molecule surface into regions of high and low electronegativity.A basic understanding of reaction mechanisms can also come from viewing and manipulating isovalued surfaces. A theory of molecular interactions, based upon second-order perturbation theory, provides for the decomposition of the intermolecular interaction energy into steric, electrostatic and orbital interactions. Color figures illustrate the docking of reactant molecular densities, electrostatic potentials and orbitals on low-energy pathways. The figures are used to visualize the steric, electrostatic and orbital contributions to molecular interaction energy. The visualization not only identifies low-energy reaction pathways, but it frequently reveals local interactions which determine the magnitude of the total interaction energy. Similar insight is not easily obtained by simple evaluation of the total interaction energy. Approximate transition states, built from structures along low-energy approach pathways, are excellent starting points for transition state searches.CAChe^TM Technical Report No. 1, Tektronix Inc., Beaverton, OR.     

Extremely localized molecular orbitals (ELMO): a non-orthogonal Hartree-Fock method

A new optimization method for extremely localized molecular orbitals (ELMO) is derived in a non-orthogonal formalism. The method is based on a quasi Newton-Raphson algorithm in which an approximate diagonal-blocked Hessian matrix is calculated through the Fock matrix. The Hessian matrix inverse is updated at each iteration by a variable metric updating procedure to account for the intrinsically small coupling between the orbitals. The updated orbitals are obtained with approximately n ² operations. No n ³ processes such as matrix diagonalization, matrix multiplication or orbital orthogonalization are employed. The use of localized orbitals allows for the creation of high-quality initial “guess” orbitals from optimized molecular orbitals of small systems and thus reduces the number of iterations to converge. The delocalization effects are included by a Jacobi correction (JC) which allows the accurate calculation of the total energy with a limited number of operations. This extension, referred to as ELMO(JC), is a variational method that reproduces the Hartree-Fock (HF) energy with an error of less than 2 kcal/mol for a reduced total cost compared to standard HF methods. The small number of variables, even for a very large system, and the limited number of operations potentially makes ELMO a method of choice to study large systems. Received: 30 December 1996 / Accepted: 5 June 1997