Highly accurate equilibrium structure of the C2h symmetric N1‐to‐O2 hydrogen‐bonded uracil‐dimer

The highly accurate ab initio equilibrium geometry of the hydrogen‐bonded uracil dimer is derived using a composite geometry extrapolation scheme based on all‐electron, complete basis set extrapolated Møller–Plesset perturbation theory using the jun‐pwCV[T,Q]Z basis sets combined with a valence CCSD(T)/cc‐pVTZ high‐level correction. Geometrical changes on dimerization are discussed and the performance of the several density functional approximations (among others SCAN, ωB97M‐V, DSD‐PBEP86‐D3(BJ), and DSD‐PBEP86‐NL) is evaluated. Orbital‐optimized MP2.5 is discussed as a reduced‐cost alternative to the CCSD(T) gradient in the composite scheme. A new reference interaction energy is calculated with explicitly correlated F12‐CCSD(T).     

A Comparative Study of Hydrogen Bonding Using Density Functional Theory

The hydrogen‐bond energies of water dimer and water‐formaldehyde complexes have been studied using density functional theory (DFT). Basis sets up to aug‐cc‐pVXZ (X=D, T, Q) were used. It was found that counterpoise corrected binding energies using the aug‐cc‐pVDZ basis set are very close to those predicted with the aug‐cc‐pVQZ set. Comparative studies using various DFT functionals on these two systems show that results from B3LYP, mPW1PW91 and PW91PW91 functionals are in better agreements with those predicted using high‐level ab initio methods. These functionals were applied to the study of hydrogen bonding between guanine (G) and cytosine (C), and between adenine (A) and thy mine (T) base pairs. With the aug‐cc‐pVDZ basis set, the predicted binding energies of base pairs are in good agreement with the most elaborate ab initio results.     

Highly accurate benchmark calculations of the interaction energies in the complexes C6H6···C6X6 (X = F,Cl, Br,and I)

Different from the case of the benzene dimer, the differences between the interaction energies are always less than 0.50 kcal/mol for face‐to‐face eclipsed, face‐to‐face staggered, and parallel‐displaced configurations of all investigated complexes C₆H₆···C₆X₆ (X = F, Cl, Br, and I). Hence, it is a great challenge for quantum chemists to accurately calculate the interaction energies for the three configurations of the complexes C₆H₆···C₆X₆. This work demonstrates that results obtained with the PBE0 density functional combined with the D3 dispersion correction (PBE0‐D3) and the basis set def2‐TZVPP are in excellent agreement with the estimates of the coupled‐cluster singles, doubles, and perturbative triples [CCSD(T)] complete basis set (CBS) limit. The other finding in this study is that, in comparison with the gold‐standard CCSD(T)/CBS benchmark, the spin‐component scaled (SCS) zeroth‐order symmetry‐adapted perturbation theory (SAPT0), when paired with the basis set aug‐cc‐pVDZ, performs also very well, and its performance is even better than that of the PBE0‐D3/def2‐TZVPP method or the conventional SAPT/aug‐cc‐pVQZ method. The findings of this study are very significant because both PBE0‐D3/def2‐TZVPP and SCS‐SAPT0/aug‐cc‐pVDZ can deal with the systems with more than 200 atoms.     

Model studies of the optical rotation,and theoretical determination of its sign for β‐pinene and trans‐pinane

We have carried out extensive studies on the basis set dependence of the calculated specific optical rotation (OR) in molecules at the level of the time–dependent Hartree–Fock and density functional approximations. To reach the limits of the basis set saturation, we have devised an artificial model, the asymmetrically deformed (chiral) methane (CM) molecule. This small system permits to use basis sets which are prohibitively large for real chiral molecules and yet shows all the important features of the basis set dependence of the OR values. The convergence of the OR has been studied with n‐aug‐cc‐pVXZ basis sets of Dunning up to the 6–ζ. In a parallel series of calculations, we have used the recently developed large polarized (LPolX) basis sets. The relatively small LPolX sets have been shown to be competitive to very large n‐aug‐cc‐pVXZ basis sets. The conclusions reached in calculations of OR in CM concerning the usefulness of LPolX basis sets have been further tested on (S)‐methyloxirane and (S)‐fluoro‐oxirane. The smallest set of the LPolX family (LPol–ds) has been found to yield OR values of similar quality as those obtained with much larger Dunning's aug‐cc‐pVQZ basis set. These results have encouraged us to carry out the OR calculations with LPol–ds basis sets for systems as large as β‐pinene and trans‐pinane. In both cases, our calculations have lead to the correct sign of the OR value in these molecules. This makes the relatively small LPol–ds basis sets likely to be useful in OR calculations for large molecules. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

PW86–PW91 density functional calculation of vertical ionization potentials: Some implications for present‐day functionals

A total of 181 vertical ionization potentials (VIPs) of 41 molecules were calculated by density functional theory (DFT) employing the Perdew–Wang 1986 (PW86) exchange and Perdew–Wang 1991 (PW91) correlation functionals and using the aug‐cc‐pV5Z basis and experimental ground‐state geometries. The overall average absolute deviation (AAD) from experiment was found to be 0.55 eV and only 0.31 eV for linear molecules but 0.86 eV for nonplanar molecules. A number of VIPs were in error by over 2 eV. In particular, DFT performed most poorly when ionization was from an orbital with highly varying density gradients (which arise from the orbital's shape or compactness or through its density being distributed over a number of atoms). Indications are that many or all present‐day functionals suffer from the same failings. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 34–52, 2001     

Comparative study on formamide–water complex

The hydrogen bonding of 1:1 complexes formed between formamide and water molecule have been investigated systematically using Hartree–Fock (HF), hybrid density functional theory (B3LYP), and post‐Hartree–Fock (MP2 and CCSD(T)) methods with range of basis sets 6‐31G(d), cc‐pVXZ (X = D, T, Q) and aug‐cc‐pVYZ (Y = D, T). Three stable structures are considered on the potential energy surface of formamide and water system. The optimized geometric parameters and interaction energies for various isomers at different levels are estimated. The IR frequencies, intensities, and frequency shifts are reported. This study shows that B3LYP/aug‐cc‐pVDZ method gives better performance for formamide‐water complexes. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010.     

Many‐body interactions of carbon monoxide cyclic oligomers: A computational study

Structural properties and energetics for the optimized carbon monoxide cyclic oligomers are analyzed at the correlated ab initio second‐order Møller–Plesset (MP2) and density functional methods (B3LYP and mPW1PW), using Dunning's cc‐pVXZ (X = D, T, Q) basis set, augmented with diffuse functions. Many‐body interactions of the stable carbon monoxide cyclic oligomers, (CO)₄ and (CO)₅ are computed at the MP2/aug‐cc‐pVTZ level. Contributions of two‐ to five‐body terms to each of these oligomers for their interaction energies, including corrections for basis set superposition error (BSSE) are investigated by using function counterpoise and its generalized version. It has been found that three‐body terms are attractive in nature and essential in order to describe the cooperative effects in the stable cyclic CO oligomers. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

On the performance of long‐range‐corrected density functional theory and reduced‐size polarized LPol‐n basis sets in computations of electric dipole (hyper)polarizabilities of π‐conjugated molecules

Static longitudinal electric dipole (hyper)polarizabilities are calculated for six medium‐sized π‐conjugated organic molecules using recently developed LPol‐n basis set family to assess their performance. Dunning's correlation‐consistent basis sets of triple‐ζ quality combined with MP2 method and supported by CCSD(T)/aug‐cc‐pVDZ results are used to obtain the reference values of analyzed properties. The same reference is used to analyze (hyper)polarizabilities predicted by selected exchange‐correlation functionals, particularly those asymptotically corrected. © 2012 Wiley Periodicals, Inc.     

An analytical approach for computing Franck‐Condon integrals of harmonic oscillators with arbitrary dimensions

We have developed an analytical approach for computing Franck‐Condon integrals (FCIs) of harmonic oscillators (HOs) with arbitrary dimensions in which the mode‐mixing Duschinsky effect is taken into account. A general formula of FCIs of HOs was obtained and was applied to study the photoelectron spectroscopy of vinyl alcohol and ovalene (C₃₂H₁₄). The equilibrium geometries, harmonic vibrational frequencies and normal modes of vinyl alcohol, ovalene, and their cations were computed at the B3LYP/aug‐cc‐pVTZ or the B3LYP/6‐31G(d) level, from which Franck‐Condon factors were calculated and photoelectron spectra were simulated. The adiabatic ionization energies of vinyl alcohol were also computed by extrapolating the CCSD(T) energies to the complete basis set limit with aug‐cc‐pVXZ (X = D, T, Q, 5). The simulated photoelectron spectra of both vinyl alcohol and ovalene are in agreement with the experiments. The computed adiabatic ionization energies of syn‐ and anti‐vinyl alcohol are in consistent with the experiment within 0.008 eV and 0.014 eV, respectively. We show, for the first time, that the analytical approach of computing FCIs is also efficient and promising for the studies of vibronic spectra of macrosystems. © 2012 Wiley Periodicals, Inc.     

Computational NMR coupling constants: Shifting and scaling factors for evaluating 1JCH

Optimized shifting and/or scaling factors for calculating one‐bond carbon–hydrogen spin–spin coupling constants have been determined for 35 combinations of representative functionals (PBE, B3LYP, B3P86, B97‐2 and M06‐L) and basis sets (TZVP, HIII‐su3, EPR‐III, aug‐cc‐pVTZ‐J, ccJ‐pVDZ, ccJ‐pVTZ, ccJ‐pVQZ, pcJ‐2 and pcJ‐3) using 68 organic molecular systems with 88 ¹J_CH couplings including different types of hybridized carbon atoms. Density functional theory assessment for the determination of ¹J_CH coupling constants is examined, comparing the computed and experimental values. The use of shifting constants for obtaining the calculated coupling improves substantially the results, and most models become qualitatively similar. Thus, for the whole set of couplings and for all approaches excluding those using the M06 functional, the root‐mean‐square deviations lie between 4.7 and 16.4 Hz and are reduced to 4–6.5 Hz when shifting constants are considered. Alternatively, when a specific rovibrational contribution of 5 Hz is subtracted from the experimental values, good results are obtained with PBE, B3P86 and B97‐2 functionals in combination with HIII‐su3, aug‐cc‐pVTZ‐J and pcJ‐2 basis sets. Copyright © 2013 John Wiley & Sons, Ltd.     

Accurate spectroscopic constants of the lowest three electronic states in halonitrenes with multireference configuration interaction

Internally contracted multireference configuration interaction (icMRCI) calculations of the ground state (X³Σ⁻), the first excited state (a¹Δ) as well as the second excited state (b¹Σ⁺) have been performed for a series of halogenated nitrenes NXs (X = Cl, Br, and I). Accurate spectroscopic constants of these lowest three electronic states of each NX were obtained in this work using MRCI methods with aug‐cc‐pVXZ (X = T, Q, 5) basis sets and complete basis set (CBS) limit. In addition, various corrections, including the Davidson correction, scalar relativistic effect, core‐valence correlation, and spin‐orbit coupling effect, have been studied in calculating spectroscopic constants, especially for heavy‐atom nitrenes. Comparisons have been made with previous computational and experimental results where available. The icMRCI + Q calculations presented in this work provide a comprehensive series of results at a consistent high level of theory for all of the halogenated nitrenes.     

IR Low‐Temperature Matrix and ab Initio Study on β‐Alanine Conformers

For the first time the argon‐matrix low‐temperature FTIR spectra of β‐alanine are recorded. They reveal a quite complicated spectral pattern which suggests the presence of several β‐alanine conformers in the matrix. To interpret the spectra, the eighteen β‐alanine conformers, stable in the gas phase, are estimated at the B3LYP and MP2 levels combined with the aug‐cc‐pVDZ. Ten low‐energy structures are reoptimized at the QCISD/aug‐cc‐pVDZ and B3LYP and MP2 levels by using the aug‐cc‐pVTZ basis sets. Assignment of the experimental spectra is undertaken on the basis of the calculated B3LYP/aug‐cc‐pVDZ anharmonic IR frequencies as well as careful estimation of the conformer population. The presence of at least three β‐alanine conformers is demonstrated. The detailed analysis of IR spectra points to the possible presence of five additional β‐alanine conformers.     

Characterization of the potential energy surfaces of two small but challenging noncovalent dimers: (P2)2 and (PCCP)2

This work characterizes eight stationary points of the P₂ dimer and six stationary points of the PCCP dimer, including a newly identified minimum on both potential energy surfaces. Full geometry optimizations and corresponding harmonic vibrational frequencies were computed with the second‐order Møller–Plesset (MP2) electronic structure method and six different basis sets: aug‐cc‐pVXZ, aug‐cc‐pV(X+d)Z, and aug‐cc‐pCVXZ where X = T, Q. A new L‐shaped structure with C₂ symmetry is the only minimum for the P₂ dimer at the MP2 level of theory with these basis sets. The previously reported parallel‐slipped structure with C_2h symmetry and a newly identified cross configuration with D₂ symmetry are the only minima for the PCCP dimer. Single point energies were also computed using the canonical MP2 and CCSD(T) methods as well as the explicitly correlated MP2‐F12 and CCSD(T)‐F12 methods and the aug‐cc‐pVXZ (X = D, T, Q, 5) basis sets. The energetics obtained with the explicitly correlated methods were very similar to the canonical results for the larger basis sets. Extrapolations were performed to estimate the complete basis set (CBS) limit MP2 and CCSD(T) binding energies. MP2 and MP2‐F12 significantly overbind the P₂ and PCCP dimers relative to the CCSD(T) and CCSD(T)‐F12 binding energies by as much as 1.5 kcal mol^?1 for the former and 5.0 kcal mol^?1 for the latter at the CBS limit. The dominant attractive component of the interaction energy for each dimer configuration was dispersion according to several symmetry‐adapted perturbation theory analyses. © 2014 Wiley Periodicals, Inc.     

New global potential energy surface of the MgH2 system and dynamics studies of the reaction H + MgH → Mg + H2

A new global potential energy surface for the ground state of MgH₂ was constructed using the permutation invariant polynomial neural network method. About 70 000 ab initio energy points were calculated via the multi‐reference configuration interaction method method with aug‐cc‐pVTZ and aug‐cc‐pVQZ basis sets, and these points were used to construct the potential energy surface (PES). To avoid basis set superposition error, the basis set was extrapolated to the complete basis set limit using the two point energy extrapolation formula. The root mean square error of the present PES is only 8.85 meV. Initial state (v = 0, j = 0) dynamics studies were performed using the time‐dependent wave packet method with a second‐order split operator for the total angular momentum J up to a value of 50. Furthermore, the reaction probability, integral cross section, and thermal rate constant are reported and compared with available theoretical studies.     

Highly Accurate CCSD(T) and DFT–SAPT Stabilization Energies of H‐Bonded and Stacked Structures of the Uracil Dimer

The CCSD(T) interaction energies for the H‐bonded and stacked structures of the uracil dimer are determined at the aug‐cc‐pVDZ and aug‐cc‐pVTZ levels. On the basis of these calculations we can construct the CCSD(T) interaction energies at the complete basis set (CBS) limit. The most accurate energies, based either on direct extrapolation of the CCSD(T) correlation energies obtained with the aug‐cc‐pVDZ and aug‐cc‐pVTZ basis sets or on the sum of extrapolated MP2 interaction energies (from aug‐cc‐pVTZ and aug‐cc‐pVQZ basis sets) and extrapolated ΔCCSD(T) correction terms [difference between CCSD(T) and MP2 interaction energies] differ only slightly, which demonstrates the reliability and robustness of both techniques. The latter values, which represent new standards for the H‐bonding and stacking structures of the uracil dimer, differ from the previously published data for the S22 set by a small amount. This suggests that interaction energies of the S22 set are generated with chemical accuracy. The most accurate CCSD(T)/CBS interaction energies are compared with interaction energies obtained from various computational procedures, namely the SCS–MP2 (SCS: spin‐component‐scaled), SCS(MI)–MP2 (MI: molecular interaction), MP3, dispersion‐augmented DFT (DFT–D), M06–2X, and DFT–SAPT (SAPT: symmetry‐adapted perturbation theory) methods. Among these techniques, the best results are obtained with the SCS(MI)–MP2 method. Remarkably good binding energies are also obtained with the DFT–SAPT method. Both DFT techniques tested yield similarly good interaction energies. The large magnitude of the stacking energy for the uracil dimer, compared to that of the benzene dimer, is explained by attractive electrostatic interactions present in the stacked uracil dimer. These interactions force both subsystems to approach each other and the dispersion energy benefits from a shorter intersystem separation.     

Assessment of dispersion‐improved exchange‐correlation functionals for the simulation of CO2 binding by alcoholamines

In this study, 12 bound complexes were selected to construct a database for testing 15 dispersion‐improved exchange‐correlation (XC) functionals, including hybrid generalized gradient approximation (GGA), modified using the Grimme's pairwise strategy, and double hybrid XC functionals, for specifically characterizing the CO₂ binding by alcoholamines. Bound complexes were selected based on the characteristics of their hydrogen bonds, dispersion, and electrostatic (particularly between the positive charge of CO₂ and the lone pair of N of alcoholamines) interactions. The extrapolated binding energy from the aug‐cc‐pVTZ (ATZ) to aug‐cc‐pVQZ (AQZ) basis set at the CCSD(T)/CBS(MP2+DZ) level was used as the reference for the XC functional comparison. M06‐2X produced the optimal agreement if the optimized geometries at MP2/ATZ level were chosen for all the test bound complexes. However, M06‐L, ωB97X, and ωB97, and were preferred if the corresponding density functional theory (DFT) optimized geometries were adapted for the benchmark. Simple bimolecular reaction between CO₂ and monoethanolamine simulated using polarizable continuum solvation model confirmed that ωB97, ωB97X, and ωB97XD qualitatively reproduced the energetics of MP2 level. The inconsistent performance of the tested XC functionals, observed when using MP2 or DFT optimized geometries, raised concerns regarding using the single‐point ab initio correction combined with DFT optimized geometry, particularly for determining the nucleophilic attack by alcoholamines to CO₂. © 2014 Wiley Periodicals, Inc.     

Theoretical study of spectroscopic and molecular properties of several low‐lying electronic states of CO molecule

The potential energy curves (PECs) of eight low‐lying electronic states (X¹Σ⁺, a³Π, a′³Σ⁺, d³Δ, e³Σ^?, A¹Π, I¹Σ^?, and D¹Δ) of the carbon monoxide molecule have been studied by an ab initio quantum chemical method. The calculations have been performed using the complete active space self‐consistent field method, which is followed by the valence internally contracted multireference configuration interaction (MRCI) approach in combination with the correlation‐consistent aug‐cc‐pV5Z basis set. The effects on the PECs by the core‐valence correlation and relativistic corrections are included. The way to consider the relativistic corrections is to use the third‐order Douglas–Kroll Hamiltonian approximation at the level of a cc‐pV5Z basis set. Core‐valence correlation corrections are performed using the cc‐pCVQZ basis set. To obtain more reliable results, the PECs determined by the MRCI calculations are corrected for size‐extensivity errors by means of the Davidson modification (MRCI+Q). The spectroscopic parameters (D_e, T_e, R_e, ω_e, ω_ex_e, ω_ey_e, B_e, α_e, and γ_e) of these electronic states are calculated using these PECs. The spectroscopic parameters are compared with those reported in the literature. Using the Breit–Pauli operator, the spin–orbit coupling effect on the spectroscopic parameters is discussed for the a³Π electronic state. With the PECs obtained by the MRCI+Q/aug‐cc‐pV5Z+CV+DK calculations, the complete vibrational states of each electronic state have been determined. The vibrational manifolds have been calculated for each vibrational state of each electronic state. The vibrational level G(ν), inertial rotation constant B_ν, and centrifugal distortion constant D_ν of the first 20 vibrational states when the rotational quantum number J equals zero are reported and compared with the experimental data. Comparison with the measurements demonstrates that the present spectroscopic parameters and molecular constants determined by the MRCI+Q/aug‐cc‐pV5Z+CV+DK calculations are both reliable and accurate. © 2012 Wiley Periodicals, Inc.     

Structure and NMR spectra of cyclophanes with unsaturated bridges (cyclophenes)

The calculated structures of several known and hypothetical cyclophanes with ethylene bridges (cyclophenes) are reported together with experimental and calculated values of their NMR parameters. Of the exchange‐correlation functionals and basis sets used in this work, only the ωB97X‐D/6‐311++G(2d,2p) and ωB97X‐D/cc‐pVQZ yielded values of the C_sp3–C_sp3 bond length close to the experimental data, although significant differences still remain. As far as the NMR parameters are concerned, except for close‐lying signals, chemical shifts and coupling constants calculated at the ωB97X‐D/cc‐pVQZ level reproduce in most cases the experimental trends. Contrary to the calculations of geometries, an agreement between the values of the NMR parameters obtained at ωB97X‐D/cc‐pVQZ level and the experimental ones is the poorest compared with that of the ωB97X‐D/6‐311++G(2d,2p) one. Taking into account that the results of the different calculations show the same qualitative trends in most cases, we believe that they correctly describe the structure and properties of the hypothetical molecules studied here. Copyright © 2012 John Wiley & Sons, Ltd.     

Ab initio calculations of the Ar–ethane intermolecular potential energy surface using bond function basis sets

The intermolecular potential energy surface (PES) of argon with ethane has been studied by ab initio calculations at the levels of second‐order Møller–Plesset perturbation (MP2) theory and coupled‐cluster theory with single, double, and noniterative triple configurations (CCSD(T)) using a series of augmented correlation‐consistent basis sets. Two sets of bond functions, bf1 (3s3p2d) and bf2 (6s6p4d2f), have been added to the basis sets to show a dramatic and systematic improvement in the convergence of the entire PES. The PES of Ar–ethane is characterized by a global minimum at a near T‐shaped configuration with a well depth of 0.611 kcal mol^?1, a second minimum at a collinear configuration with a well depth of 0.456 kcal mol^?1, and a saddle point connecting the two minima. It is shown that an augmented correlation‐consistent basis set with a set of bond functions, either bf1 or bf2, can effectively produce results equivalent to the next larger augmented correlation‐consistent basis set, that is, aug‐cc‐pVDZ‐bf1 ≈ aug‐cc‐pVTZ, aug‐cc‐pVTZ‐bf1 ≈ aug‐cc‐pVQZ. Very importantly, the use of bond functions improves the PES globally, resulting accurate potential anisotropy. Finally, MP2 method is inadequate for accurate calculations, because it gives a potentially overestimated well depth and, more seriously, a poor potential anisotropy. © 2012 Wiley Periodicals, Inc.     

Spin‐component‐scaled double hybrids: An extensive search for the best fifth‐rung functionals blending DFT and perturbation theory

Following up on an earlier preliminary communication (Kozuch and Martin, Phys. Chem. Chem. Phys. 2011, 13 , 20104), we report here in detail on an extensive search for the most accurate spin‐component‐scaled double hybrid functionals [of which conventional double hybrids (DHs) are a special case]. Such fifth‐rung functionals approach the performance of composite ab initio methods such as G3 theory at a fraction of their computational cost, and with analytical derivatives available. In this article, we provide a critical analysis of the variables and components that maximize the accuracy of DHs. These include the selection of the exchange and correlation functionals, the coefficients of each component [density functional theory (DFT), exact exchange, and perturbative correlation in both the same spin and opposite spin terms], and the addition of an ad‐hoc dispersion correction; we have termed these parametrizations “DSD‐DFT” (Dispersion corrected, Spin‐component scaled, Double‐hybrid DFT). Somewhat surprisingly, the quality of DSD‐DFT is only mildly dependent on the underlying DFT exchange and correlation components, with even DSD‐LDA yielding respectable performance. Simple, nonempirical GGAs appear to work best, whereas meta‐GGAs offer no advantage (with the notable exception of B95c). The best correlation components appear to be, in that order, B95c, P86, and PBEc, while essentially any good GGA exchange yields nearly identical results. On further validation with a wider variety of thermochemical, weak interaction, kinetic, and spectroscopic benchmarks, we find that the best functionals are, roughly in that order, DSD‐PBEhB95, DSD‐PBEP86, DSD‐PBEPW91, and DSD‐PBEPBE. In addition, DSD‐PBEP86 and DSD‐PBEPBE can be used without source code modifications in a wider variety of electronic structure codes. Sample job decks for several commonly used such codes are supplied as electronic Supporting Information. Copyright © 2013 Wiley Periodicals, Inc.