Accuracy of spectroscopic constants of diatomic molecules from ab initio calculations

Complex potential energy surfaces with projected CAP technique: Vibrational excitation of N2

The projected complex absorbing potential (CAP) technique is one of the methods that allow one to extend the bound state methods for computing resonances' energies and widths. Here, we explore the accuracy of the potential energy curves generated with different electronic structure theory methods in combination with the projected CAP technique by considering resonant vibrational excitation (RVE) of N2 by electron impact as a model process. We report RVE cross sections computed using the boomerang model with potential energy curves obtained with CAP-based extended multistate complete active space perturbation theory (XMS-CASPT2) and equation of motion coupled-cluster method for electron attachment with single and double substitution (EOM-EA-CCSD) methods. We also compare potential energy curves computed with several electronic structure methods, including XMS-CASPT2, EOM-EA-CCSD, multireference configuration interaction with singles (MR-CIS) and singles and doubles (MR-CISD). A good agreement is observed between the experiment and simulated RVE cross sections obtained with the potential energy curves generated with XMS-CASPT2 and EOM-EA-CCSD methods, thus highlighting the potential of the projected CAP technique combined with accurate electronic structure methods for dynamical simulations of the processes that proceed through metastable electronic states.

On the performance of composite schemes in determining equilibrium molecular structures

Determination of equilibrium molecular structures is an essential ingredient in predicting spectroscopic parameters that help in identifying molecular carriers of microwave transitions. Here, the performance of two different ab initio composite approaches for obtaining equilibrium structures, "energy scheme" and "geometry scheme," is explored and compared to semi-experimental equilibrium structures. This study is performed for a set of 11 molecules which includes diatomics, linear triatomics, and a few non-linear molecules. The ab initio calculations were performed using three tiers of composite chemical recipes. The current results show that as the overall rigor of calculation is increased, the semi-experimental and the ab initio numbers agree to within 0.0003 Å for all molecules in the test set. The composite approach based on correcting the potential energy surface (energy scheme) and the one based on correcting the geometry directly (geometry scheme) show excellent agreement with each other. This work represents a step toward development of efficient and highly accurate procedures for computing ab initio equilibrium structures.

Accurate Property Prediction by Second Order Perturbation Theory: The REMP and OO-REMP Hybrids

The prediction of molecular properties such as equilibrium structures or vibrationalwavenumbers is a routine task in computational chemistry. If very high accuracy is required, however, the use of computationally demanding ab initio wavefunction methods is mandatory. We present property calculations utilizing the REMP and OO-REMP hybrid perturbation theories showing that with the latter approach, very accurate results are obtained at second order in perturbation theory. Specifically, equilibrium structures and harmonic vibrational wavenumbers as well as dipole moments of closed and open shell molecules were calculated and compared to the best available experimental results or very accurate calculations.OO-REMP is capable of predicting bond lengths of small closed and open shell molecules with an accuracy of 0.2 pm and 0.5 pm, respectively, often within the range of experimental uncertainty. Equilibrium harmonic vibrational wavenumbers are predicted with an accuracy better than 20 cm−1 . Dipole moments of small closed and open shell molecules are reproduced with a relative error of less than 3 %. Across all investigated properties it turns out that a 20 %:80 % MP:RE mixing ratio consistently provides the best results. This is in line with our previous findings featuring closed and open shell reaction energies.

Atomic Mean-Field Approach within Exact Two-Component Theory Based on the Dirac-Coulomb-Breit Hamiltonian

An extension of the exact two-component theory with atomic mean-field integrals (the X2CAMF scheme) to the treatment of the Breit term together with efficient implementation using an atomic Dirac-Coulomb-Breit Hartree-Fock program is reported. The accuracy of the X2CAMF scheme for treating the contributions from the Breit term to the molecular properties is demonstrated using benchmark calculations of equilibrium bond lengths, harmonic frequencies, and dipole moments for molecules containing elements across the periodic table. Calculations of the properties for molecules containing period four elements aiming at high accuracy as well as for Th- and U-containing molecules are also presented and compared with experimental results to demonstrate the usefulness of the X2CAMF scheme in combination with accurate treatments of electron correlation by the coupled-cluster (CC) methods. The combination of CC methods and the X2CAMF scheme shows potential to extend the accuracy of CC calculations to heavy elements, e.g., to computational heavy-element thermochemistry.

Quantitative Comparison of Experimental and Computed IR-Spectra Extracted from Ab Initio Molecular Dynamics

Experimentally measured infrared spectra are often compared to their computed equivalents. However, the accordance is typically characterized by visual inspection, which is prone to subjective judgment. The primary challenge for a similarity-based analysis is that the artifacts introduced by each approach are very different and, therefore, may require preprocessing steps to determine and correct impeding irregularities. To allow for automated objective assessment, we propose a practical and comprehensive workflow involving scaling factors, a novel baseline correction scheme, and peak smoothing. The resulting spectra can then easily be compared quantitatively using similarity measures, for which we found the Pearson correlation coefficient to be the most suitable. The proposed procedure is then applied to compare the agreement of the experimental infrared spectra from the NIST Chemistry Web book with the calculated spectra using standard harmonic frequency analysis and spectra extracted from ab initio molecular dynamics simulations at different levels of theory. We conclude that the direct, quantitative comparison of calculated and measured IR spectra might become a novel, sophisticated approach to benchmark quantum-chemical methods. In the present benchmark, simulated spectra based on ab initio molecular dynamics show in general better agreement with the experiment than static calculations.

Convergence of nuclear magnetic shieldings and one-bond 1 J(11 B1 H) indirect spin-spin coupling constants in small boron molecules

Self-consistent field Hartree-Fock, density functional theory, and coupled-cluster calculations of the nuclear magnetic shielding constants of BH and BH₃ molecules have been conducted to characterize the convergence of individual results obtained with correlation-consistent and polarization-consistent basis sets. The individual ¹¹ B and ¹ H NMR parameters were estimated in the complete basis set limit and compared with benchmark results. Only the KT3 density functional accurately reproduced ¹¹ B shielding in BH molecule.

Theoretical studies of atmospheric molecular complexes interacting with NIR to UV light

Theoretical structural and spectroscopic data for weakly bonded atmospheric complexes of formaldehyde interacting with Earth’s electromagnetic spectrum.

Computing Bulk Phase Raman Optical Activity Spectra from ab initio Molecular Dynamics Simulations

We present our novel methodology for computing Raman optical activity (ROA) spectra of liquid systems from ab initio molecular dynamics (AIMD) simulations. The method is built upon the recent developments to obtain magnetic dipole moments from AIMD and to integrate molecular properties by using radical Voronoi tessellation. These techniques are used to calculate optical activity tensors for large and complex periodic bulk phase systems. Only AIMD simulations are required as input, and no time-consuming perturbation theory is involved. The approach relies only on the total electron density in each time step and can readily be combined with a wide range of electronic structure methods. To the best of our knowledge, these are the first computed ROA spectra for a periodic bulk phase system. As an example, the experimental ROA spectrum of liquid (R)-propylene oxide is reproduced very well.

Analytic calculations of anharmonic infrared and Raman vibrational spectra

Using a recently developed recursive scheme for the calculation of high-order geometric derivatives of frequency-dependent molecular properties [Ringholm et al., J. Comp. Chem., 2014, 35, 622], we present the first analytic calculations of anharmonic infrared (IR) and Raman spectra including anharmonicity both in the vibrational frequencies and in the IR and Raman intensities. In the case of anharmonic corrections to the Raman intensities, this involves the calculation of fifth-order energy derivatives-that is, the third-order geometric derivatives of the frequency-dependent polarizability. The approach is applicable to both Hartree-Fock and Kohn-Sham density functional theory. Using generalized vibrational perturbation theory to second order, we have calculated the anharmonic infrared and Raman spectra of the non- and partially deuterated isotopomers of nitromethane, where the inclusion of anharmonic effects introduces combination and overtone bands that are observed in the experimental spectra. For the major features of the spectra, the inclusion of anharmonicities in the calculation of the vibrational frequencies is more important than anharmonic effects in the calculated infrared and Raman intensities. Using methanimine as a trial system, we demonstrate that the analytic approach avoids errors in the calculated spectra that may arise if numerical differentiation schemes are used.

Estimation of isotropic nuclear magnetic shieldings in the CCSD(T) and MP2 complete basis set limit using affordable correlation calculations

A linear correlation between isotropic nuclear magnetic shielding constants for seven model molecules (CH2 O, H2 O, HF, F2 , HCN, SiH4 and H2 S) calculated with 37 methods (34 density functionals, RHF, MP2 and CCSD(T)), with affordable pcS-2 basis set and corresponding complete basis set results, estimated from calculations with the family of polarization-consistent pcS-n basis sets is reported. This dependence was also supported by inspection of profiles of deviation between CBS estimated nuclear shieldings and shieldings obtained with the significantly smaller basis sets pcS-2 and aug-cc-pVTZ-J for the selected set of 37 calculation methods. It was possible to formulate a practical approach of estimating the values of isotropic nuclear magnetic shielding constants at the CCSD(T)/CBS and MP2/CBS levels from affordable CCSD(T)/pcS-2, MP2/pcS-2 and DFT/CBS calculations with pcS-n basis sets. The proposed method leads to a fairly accurate estimation of nuclear magnetic shieldings and considerable saving of computational efforts.

Benchmarking density-functional theory calculations of NMR shielding constants and spin-rotation constants using accurate coupled-cluster calculations

Accurate sets of benchmark nuclear-magnetic-resonance shielding constants and spin-rotation constants are calculated using coupled-cluster singles-doubles (CCSD) theory and coupled-cluster singles-doubles-perturbative-triples [CCSD(T)] theory, in a variety of basis sets consisting of (rotational) London atomic orbitals. The accuracy of the calculated coupled-cluster constants is established by a careful comparison with experimental data, taking into account zero-point vibrational corrections. Coupled-cluster basis-set convergence is analyzed and extrapolation techniques are employed to estimate basis-set-limit quantities, thereby establishing an accurate benchmark data set. Together with the set provided for rotational g-tensors and magnetizabilities in our previous work [O. B. Lutnæs, A. M. Teale, T. Helgaker, D. J. Tozer, K. Ruud, and J. Gauss, J. Chem. Phys. 131, 144104 (2009)], it provides a substantial source of consistently calculated high-accuracy data on second-order magnetic response properties. The utility of this benchmark data set is demonstrated by examining a wide variety of Kohn-Sham exchange-correlation functionals for the calculation of these properties. None of the existing approximate functionals provide an accuracy competitive with that provided by CCSD or CCSD(T) theory. The need for a careful consideration of vibrational effects is clearly illustrated. Finally, the pure coupled-cluster results are compared with the results of Kohn-Sham calculations constrained to give the same electronic density. Routes to future improvements are discussed in light of this comparison.

Accurate structure, thermodynamic and spectroscopic parameters from CC and CC/DFT schemes: the challenge of the conformational equilibrium in glycine

The structures, relative stabilities, and infrared spectra of the six low-energy conformers of glycine have been characterized using a state-of-the-art quantum-mechanical approach allowing the bond distances, conformational enthalpies and vibrational frequencies to be determined well within the chemical accuracy. Transition state structures governing interconversion among the different energy minima have also been characterized. In detail, the gas-phase thermodynamic properties (at 15 K and 410 K) of the glycine conformers considered have been obtained with a 1 kJ mol(-1) accuracy, and it has been shown that the employment of DFT geometries usually reduces such accuracy by at most 0.1 kJ mol(-1). Regarding molecular structures, the use of two different composite schemes allowed us to further confirm the suitability of a rather cost-effective approach and provide geometrical parameters with an overall accuracy better than 0.002 Å for distances and 1 degree for angles. Thanks to a hybrid CC/DFT approach, the infrared spectra of all conformers considered and of several deuterated isotopologues have been reproduced (when experimental data were available) or predicted with an accuracy of 10 cm(-1). Finally, the joint thermodynamic and spectroscopic investigation allowed us to shed some light on the possible observation of elusive conformers. On the whole, the high accuracy of the computational results allows us to draw a fully consistent interpretation of the available experimental data and to obtain a more complete characterization of the potential energy surface of glycine.

Implementation of the incremental scheme for one-electron first-order properties in coupled-cluster theory

A fully automated parallelized implementation of the incremental scheme for coupled-cluster singles-and-doubles (CCSD) energies has been extended to treat molecular (unrelaxed) first-order one-electron properties such as the electric dipole and quadrupole moments. The convergence and accuracy of the incremental approach for the dipole and quadrupole moments have been studied for a variety of chemically interesting systems. It is found that the electric dipole moment can be obtained to within 5% and 0.5% accuracy with respect to the exact CCSD value at the third and fourth orders of the expansion, respectively. Furthermore, we find that the incremental expansion of the quadrupole moment converges to the exact result with increasing order of the expansion: the convergence of nonaromatic compounds is fast with errors less than 16 mau and less than 1 mau at third and fourth orders, respectively (1 mau=10(-3)ea(0)(2)); the aromatic compounds converge slowly with maximum absolute deviations of 174 and 72 mau at third and fourth orders, respectively.

Extrapolation to the complete basis set limit of structural parameters: comparison of different approaches

For the first time, different procedures for extrapolating structural parameters to the complete basis set limit are compared. More precisely, a well-tested scheme, based on extrapolated gradients, is employed for validating the extrapolation procedure directly applied to geometrical parameters by means of different extrapolation functions. For this purpose, a set of small systems (mostly triatomics), also containing second-row elements, has been chosen.

H2O, H2, HF, F2 and F2O nuclear magnetic shielding constants and indirect nuclear spin-spin coupling constants (SSCCs) in the BHandH/pcJ-n and BHandH/XZP Kohn-Sham limits

Good performance of segmented contracted basis sets XZP, where X = D, T, Q and 5, for obtaining H(2)O, H(2), HF, F(2) and F(2)O nuclear isotropic shielding constants in the BHandH Kohn-Sham basis set limit was shown. The results of two- and three-parameter complete basis set limit extrapolation schemes were compared with experimental results, earlier literature data and benchmark ab initio results. Similar convergence patterns of shieldings obtained from calculations using general purpose XZP basis sets and from polarization-consistent basis sets pcS-n and pcJ-n, where n = 0, 1, 2, 3 and 4, designed to accurately predict magnetic properties were observed. On the contrary, the SSCCs were more sensitive to the XZP basis set size and generally less accurate than those estimated using pcJ-n basis set family. The BHandH density functional markedly outperforms B3LYP method in predicting heavy atom shieldings and SSCCs values in the studied systems.

Perturbative treatment of scalar-relativistic effects in coupled-cluster calculations of equilibrium geometries and harmonic vibrational frequencies using analytic second-derivative techniques

An analytic scheme for the computation of scalar-relativistic corrections to nuclear forces is presented. Relativistic corrections are included via a perturbative treatment involving the mass-velocity and the one-electron and two-electron Darwin terms. Such a scheme requires mixed second derivatives of the nonrelativistic energy with respect to the relativistic perturbation and the nuclear coordinates and can be implemented using available second-derivative techniques. Our implementation for Hartree-Fock self-consistent field, second-order Moller-Plesset perturbation theory, as well as the coupled-cluster level is used to investigate the relativistic effects on the geometrical parameters and harmonic vibrational frequencies for a set of molecules containing light elements (HX, X=F, Cl, Br; H2X, X=O, S; HXY, X=O, S and Y=F, Cl, Br). The focus of our calculations is the basis-set dependence of the corresponding relativistic effects, additivity of electron correlation and relativistic effects, and the importance of core correlation on relativistic effects.

The 1,2,3-Tridehydrobenzene Triradical: 2B or Not 2B? The Answer is 2A!

The molecular and electronic structure of 1,2,3-tridehydrobenzene was investigated by a variety of computational methods. The two lowest electronic states of the triradical are the (2)B(2) and (2)A(1) doublet states characterized by different interactions of the unpaired electrons. Vertically, the two states are well separated in energy-by 4.9 and 1.4 eV, respectively. However, due to different bonding patterns, their equilibrium structures are very different and, adiabatically, the two states are nearly degenerate. The adiabatic energy gap between the (2)B(2) and (2)A(1) states is estimated to be 0.7-2.1 kcal/mol, in favor of the (2)A(1) state. Harmonic vibrational frequencies and anharmonic corrections were calculated for both states. Comparison with the three experimentally observed IR transitions supports the assignment of the (2)A(1) ground state for the triradical with a weakly bonding distance of 1.67-1.69 A between the meta radical centers.

Basis Set Limit CCSD(T) Harmonic Vibrational Frequencies

Benchmark, frozen-core CCSD(T) equilibrium harmonic vibrational frequencies of 12 closed-shell and five open-shell molecules are computed to within 1 cm-1 of the basis set limit using the explicitly correlated CCSD(T)-R12 method. The convergence of the standard CCSD(T) method with the one-particle basis sets of Dunning and co-workers is examined and found to be slow, with mean and maximum absolute errors of 1.3 and 3.5 cm-1 remaining at the cc-pV6Z level. Finite basis set effects do not appear to introduce systematic errors in equilibrium harmonic frequencies, and mean absolute errors reduce by a factor of 2 for each basis set cardinal number increment. The convergence of individual equilibrium harmonic frequencies is not guaranteed to be monotonic due to the associated shift in the equilibrium structure. The inclusion of computed scalar relativistic effects and previously available corrections for core-valence correlation and higher-order excitations in the cluster operator results in an agreement with experimentally derived harmonic frequencies of 0.1, 0.3, and -0.4 cm-1 for HF, N2, and CO, respectively. F2 continues to present a challenge to computational chemistry with an error of 3.2 cm-1, primarily resulting from the high basis set dependence of the quadruples contribution.

Hydrogen bonded OH-stretching vibration in the water dimer

We have calculated the frequencies and intensities of the hydrogen-bonded OH-stretching transitions in the water dimer complex. The potential-energy curve and dipole-moment function are calculated ab initio at the coupled cluster with singles, doubles, and perturbative triples level of theory with correlation-consistent Dunning basis sets. The vibrational frequencies and wavefunctions are found from a numerical solution to a one-dimensional Schrödinger equation. The corresponding transition intensities are found from numerical integration of these vibrational wavefunctions with the ab initio calculated dipole moment function. We investigate the effect of counterpoise correcting both the potential-energy surface and dipole-moment function. We find that the effect of using a numeric potential is significant for higher overtones and that inclusion of a counterpoise correction for basis set superposition error is important.

Efficient correlation-corrected vibrational self-consistent field computation of OH-stretch frequencies using a low-scaling algorithm

The authors present a new computational scheme to perform accurate and fast direct correlation-corrected vibrational self-consistent field (CC-VSCF) computations for a selected number of vibrational modes, which is aimed at predicting a few vibrations in large molecular systems. The method is based on a systematic selection of vibrational mode-mode coupling terms, leading to the direct ab initio construction of a sparse potential energy surface. The computational scaling of the CC-VSCF computation on the generated surface is then further reduced by using a screening procedure for the correlation-correction contributions. The proposed method is applied to the computation of the OH-stretch frequency of five aliphatic alcohols. The authors investigate the influence of different pseudopotential and all-electron basis sets on the quality of the correlated potential energy surfaces computed and on the OH-stretch frequencies calculated for each surface. With the help of these test systems, the authors show that their method offers a computational scaling that is two orders of magnitude lower than a standard CC-VSCF method and that it is of equal accuracy.

Addition by subtraction in coupled-cluster theory: A reconsideration of the CC and CI interface and the nCC hierarchy

The nCC hierarchy of coupled-cluster approximations, where n guarantees exactness for n electrons and all products of n electrons are derived and applied to several illustrative problems. The condition of exactness for n=2 defines nCCSD=2CC, with nCCSDT=3CC and nCCSDTQ=4CC being exact for three and four electrons. To achieve this, the minimum number of diagrams is evaluated, which is less than in the corresponding CC model. For all practical purposes, nCC is also the proper definition of a size-extensive CI. 2CC is also an orbitally invariant coupled electron pair approximation. The numerical results of nCC are close to those for the full CC variant, and in some cases are closer to the full CI reference result. As 2CC is exact for separated electron pairs, it is the natural zeroth-order approximation for the correlation problem in molecules with other effects introduced as these units start to interact. The nCC hierarchy of approximations has all the attractive features of CC including its size extensivity, orbital invariance, and orbital insensitivity, but in a conceptually appealing form suited to bond breaking, while being computationally less demanding. Excited states from the equation of motion (EOM-2CC) are also reported, which show results frequently approaching those of EOM-CCSDT.

Intrinsic Dinitrogen Activation at Bare Metal Atoms

There is ongoing interest in metal complexes which bind dinitrogen and facilitate either its reduction or oxidation under mild conditions. In nature, the enzyme nitrogenase catalyzes this process, and dinitrogen fixation occurs under mild and ambient conditions at a metal-sulfur cluster in the center of the MoFe protein, but the mechanism of this process remains largely unknown. In the last few years, new important discoveries have been made in this field. In this review are discussed recent findings on the interaction of N(2) with metal atoms and metal-atom dimers from all groups of the periodic table as provided by gas-phase as well as matrix-isolation experiments. Intrinsic dinitrogen activation at such bare metal atoms is then related to corresponding processes at complexes, clusters, and surfaces.

The photoelectron spectrum of elusive cyclic-N3 and characterization of the potential energy surface and vibrational states of the ion

A potential energy surface is constructed for the ground X (1)A(1) electronic state of cyclic-N(3) (+) based on three-dimensional spline interpolation of ab initio points. The vibrational states of this molecular ion are calculated in the range up to 14 500 cm(-1) using hyperspherical coordinates and the coupled-channel (sector-adiabatic) approach. All the vibrational states are analyzed and assigned. The Franck-Condon overlaps of these states with the vibrational states of the neutral are calculated to predict the photoelectron spectrum of cyclic-N(3). Peak intensities are governed by the nodal structure of the vibrational wave functions and reflect the large geometric phase effect predicted for cyclic-N(3). Experimental validation may shed light on the existence of this elusive molecule and confirm the magnitude of the geometric phase effect.

Dual-basis second-order Møller-Plesset perturbation theory: A reduced-cost reference for correlation calculations

The resolution-of-the-identity (RI) approximation has placed the onus of the cost of a second-order Moller-Plesset (MP2) calculation on the underlying self-consistent field (SCF) calculation for many moderately sized molecules. A dual-basis approach to the SCF calculation, based on previous methods demonstrated for density functional theory, is combined with RI-MP2 calculations, and small basis subsets for cc-pVTZ, cc-pVQZ, and 6-311++G(3df,3pd) are presented. These subsets provide time savings of greater than 90%, with negligible errors in absolute and relative energies, compared to the associated full-basis counterpart. The method is tested with a series of rotational barriers, relative conformational energies of alanine tetrapeptides, as well as the full G3/99 molecular set. RI-MP2 calculations on alanine octapeptides (40 heavy atoms, 3460 basis functions), using cc-pVQZ, are presented. Results improve upon previous methods that diagonalize the virtual space separately.

High-Accuracy Quantum Mechanical Studies of π−π Interactions in Benzene Dimers

Although supramolecular chemistry and noncovalent interactions are playing an increasingly important role in modern chemical research, a detailed understanding of prototype noncovalent interactions remains lacking. In particular, pi-pi interactions, which are ubiquitous in biological systems, are not fully understood in terms of their strength, geometrical dependence, substituent effects, or fundamental physical nature. However, state-of-the-art quantum chemical methods are beginning to provide answers to these questions. Coupled-cluster theory through perturbative triple excitations in conjunction with large basis sets and extrapolations to the complete basis set limit have provided definitive results for the binding energy of several configurations of the benzene dimer, and benchmark-quality ab initio potential curves are being used to calibrate new density functional and force-field models for pi-pi interactions. Studies of substituted benzene dimers indicate flaws in the conventional wisdom about substituent effects in pi-pi interactions. Three-body and four-body interactions in benzene clusters have also been examined.