Accurate calculation of anharmonic vibrational frequencies of medium sized molecules using local coupled cluster methods

Zero-Point-Energy Driven Isotopic Exchange of the [H3O]- anion Probed by Mid-Infrared Action Spectroscopy

We present the first observation of vibrational transitions in the [H3O]- anion, an intermediate in the anion-molecule reaction of water, H2O, and hydride, H-, using a laser-induced isotopic H/D exchange reaction action spectroscopy scheme applied to anions. The observed bands are assigned as the fundamental and first overtone of the H2O-H- vibrational stretching mode, based on anharmonic calculations within the vibrational perturbation theory and vibrational configuration interaction. Although the D2O·D- species has the lowest energy, our experiments confirm the D2O·H- isotope to be a sink of the H/D exchange reaction. Ab initio calculations corroborate that the formation of D2O·H- is favored, as the zero-point-energy difference is larger between D2 and H2 than between D2O·H- and D2O·D-.

Performance of Vibrational Self-Consistent Field Theory for Accurate Potential Energy Surfaces: Fundamentals, Excited States, and Intensities

The performance of vibrational structure calculations beyond harmonic approximation in the framework of the vibrational self-consistent field method with second-order perturbation corrections (VSCF-PT2) is investigated in conjunction with very accurate potential energy surfaces (PESs) given by various coupled-cluster electronic structure theories. The quality of anharmonic calculations depends on the accuracy of the underlying multidimensional PES obtained from its functional form, which is given by the level of electronic structure theory. Two such highest levels of typical coupled-cluster electronic structure methods, CCSD and the ″gold standard″ CCSD(T), along with their variants such as CCD, CR-CCL (completely renormalized CR-CC(2,3) approach), and CCSD(TQ) are tested for the construction of accurate anharmonic potentials without any fitting or ad hoc scaling and using cc-pVTZ basis sets. The accuracy of VSCF-PT2 theory in comparison to experimental values is tested for a series of 16 molecules with 135 fundamental bands, 64 overtones, and combination bands and also for 39 intensities. It is found that CCD and CCSD bind the potential tighter than CCSD(T) and the computed VSCF-PT2 transitions are more blue-shifted showing higher deviation from the experiment. In general, VSCF-PT2 results computed at the CCSD(T) potential offer a good cost/accuracy ratio, with the mean absolute deviation and the mean absolute percentage error with the experiment being ∼16 cm^-1 and 1.38, respectively, for fundamentals. Additionally, while the CR-CCL and CCSD(TQ) methods offer similar levels of accuracies as compared to CCSD(T), the former offers a better accuracy/cost ratio than the latter and is a suitable alternative to CCSD(T).

Vibrational spectrum and photochemistry of phosphaketene HPCO

The vibrational spectra of the simplest phosphaketene HPCO and its isotopologue DPCO in solid Ar-matrices at 12.0 K have been analyzed with the aid of the computations at the CCSD(T)-F12a/cc-pVTZ-F12 level using configuration-selective vibrational configuration interaction (VCI). In addition to the four IR fundamentals, four overtone and ten combination bands have been unambiguously identified. Furthermore, the photochemistry of HPCO in the matrix has been investigated for the first time. Upon UV-light irradiation (365 or 266 nm), CO-elimination occurs by forming the parent phosphinidene HP that can be trapped by ˙NO to yield the elusive phosphinimine-N-oxyl radical HPNO˙. In contrast, an excimer laser (193 nm) irradiation of HPCO causes additional decomposition to H˙ and ˙PCO with concomitant formation of the long-sought phosphaethyne HOCP.

DLPNO-MP2 second derivatives for the computation of polarizabilities and NMR shieldings

We present a derivation and efficient implementation of the formally complete analytic second derivatives for the domain-based local pair natural orbital second order Møller-Plesset perturbation theory (MP2) method, applicable to electric or magnetic field-response properties but not yet to harmonic frequencies. We also discuss the occurrence and avoidance of numerical instability issues related to singular linear equation systems and near linear dependences in the projected atomic orbital domains. A series of benchmark calculations on medium-sized systems is performed to assess the effect of the local approximation on calculated nuclear magnetic resonance shieldings and the static dipole polarizabilities. Relative deviations from the resolution of the identity-based MP2 (RI-MP2) reference for both properties are below 0.5% with the default truncation thresholds. For large systems, our implementation achieves quadratic effective scaling, is more efficient than RI-MP2 starting at 280 correlated electrons, and is never more than 5-20 times slower than the equivalent Hartree-Fock property calculation. The largest calculation performed here was on the vancomycin molecule with 176 atoms, 542 correlated electrons, and 4700 basis functions and took 3.3 days on 12 central processing unit cores.

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.

Characterizing structures, energetics, and spectra of species on the 1,3[H, C, As] potential energy surfaces: A high-level theoretical contribution

The ground and the low lying electronic states of structures on the ^1,3[H, C, As] potential energy surfaces were investigated with the highly correlated theoretical approaches CCSD(T), CCSD(T)-F12b, and CASSCF/MRCI along with the series of correlation consistent (aug-cc-pVnZ, n = D, T, Q, 5) basis sets. Energetic and spectroscopic parameters were obtained at the complete basis set limit, and the effect of core-valence correlation on these properties evaluated. Fundamental frequencies were also computed with the variational configuration interaction (VCI) approach. Heats of formation at 0 and 298.15 K were estimated for HCAs and CH, AsH, CAs, and HCAs, as well as the calculation of ionization potentials for HCAs. Comparisons of the present results with literature ones for the systems HCN/HNC, HCP/HPC highlight similarities and differences among these systems. Altogether, this investigation provides a very reliable characterization of the species on the surfaces and should guide future experimental studies on these systems.

Computing Bulk Phase Resonance Raman Spectra from ab Initio Molecular Dynamics and Real-Time TDDFT

We present our novel approach for computing resonance Raman (RR) spectra of periodic bulk phase systems from ab initio molecular dynamics, including solvent influence and some anharmonic effects. Based on real-time time-dependent density functional theory, we obtain the RR spectra for all laser wavelengths in one pass. We compute the RR spectrum of uracil in aqueous solution, which is in good agreement with experiment. This is the first simulation of a bulk phase RR spectrum.

Probing vibrational coupling via a grid-based quantum approach-an efficient strategy for accurate calculations of localized normal modes in solid-state systems

In this work an approach to investigate the properties of strongly localized vibrational modes of functional groups in bulk material and on solid-state surfaces is presented. The associated normal mode vectors are approximated solely on the basis of structural information and obtained via diagonalization of a reduced Hessian. The grid-based Numerov procedure in one and two dimensions is then applied to an adequate scan of the respective potential surface yielding the associated vibrational wave functions and energy eigenvalues. This not only provides a detailed description of anharmonic effects but also an accurate inclusion of the coupling between the investigated vibrational states on a quantum mechanical level. All results obtained for the constructed normal modes are benchmarked against their analytical counterparts obtained from the diagonalization of the total Hessian of the entire system. Three increasingly complex systems treated at quantum chemical level of theory have been considered, namely the symmetric and asymmetric stretch vibrations of an isolated water molecule, hydroxyl groups bound to the surface of GeO2 (001), α-quartz(001) and Rutil (001) as well as crystalline Li2 NH serving as an example for a bulk material. While the data obtained for the individual systems verify the applicability of the proposed methodology, comparison to experimental data demonstrates the accuracy of this methodology despite the restriction to limit this methodology to a few selected vibrational modes. The possibility to investigate vibrational phenomena of localized normal modes without the requirement of executing costly harmonic frequency calculations of the entire system enables the application of this method to cases in which the determination of normal modes is prohibitively expensive or not available for a particular level of theory. © 2018 Wiley Periodicals, Inc.

Efficient and Accurate Born-Oppenheimer Molecular Dynamics for Large Molecular Systems

An efficient scheme for the calculation of Born-Oppenheimer molecular dynamics (BOMD) simulations is introduced. It combines the corrected small basis set Hartree-Fock (HF-3c) method by Sure and Grimme [J. Comput. Chem. 2013, 43, 1672], extended Lagrangian BOMD (XL-BOMD) by Niklasson et al. [J. Chem. Phys. 2009, 130, 214109], and the calculation of the two electron integrals on graphics processing units (GPUs) [J. Chem. Phys. 2013, 138, 134114; J. Chem. Theory Comput. 2015, 11, 918]. To explore the parallel performance of our strong scaling implementation of the method, we present timings and extract, as its validation and first illustrative application, high-quality vibrational spectra from simulated trajectories of β-carotene, paclitaxel, and liquid water (up to 500 atoms). We conclude that the presented BOMD scheme may be used as a cost-efficient and reliable tool for computing vibrational spectra and thermodynamics of large molecular systems including explicit solvent molecules containing 500 atoms and more. Simulating 50 ps of maitotoxin (nearly 500 atoms) employing time steps of 0.5 fs requires ∼3 weeks on 12 CPUs (Intel Xeon E5 2620 v3) with 24 GPUs (AMD FirePro 3D W8100).

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.

Effects of Global and Local Anharmonicities on the Thermodynamic Properties of Sulfuric Acid Monohydrate

We use state-of-the-art electronic structure calculation methods and large basis sets to obtain reliable values for the thermodynamic properties of sulfuric acid monohydrate and study the effects of vibrational anharmonicity on these properties. We distinguish between two forms of vibrational anharmonicity: local anharmonicity, which refers to the anharmonicity of the vibrational modes of a given cluster conformer, and global anharmonicity, which originates from accounting for the presence of different conformers in the first place. In our most accurate approach, we solve the nuclear Schrödinger equation variationally for the intermolecular large-amplitude motions, thus quantum-mechanically accounting for the presence of higher-energy conformers for both reactants and products, while using the standard vibrational perturbational approach for the other vibrational modes. This results in a value of -11.0 kJ/mol for the reaction Gibbs free energy at 298.15 K. When standard vibrational perturbational approaches are employed, the effects of local anharmonicity depend heavily on the choice of the electronic structure calculation basis set. In fact, better results can often be achieved by combining a simple harmonic treatment for the vibrational partition function with a statistical mechanical accounting of global anharmonicity. Thus, we recommend that future studies that intend to include anharmonicity start by accounting for the presence of higher-energy conformers and only then consider whether local anharmonicity calculations are feasible and necessary.

Vibrational anharmonicity of small gold and silver clusters using the VSCF method

We study the vibrational spectra of small neutral gold (Au₂–Au₁₀) and silver (Ag₂–Au₅) clusters using the vibrational self-consistent field method (VSCF) in order to account for anharmonicity.

Pushing the limits in accurate vibrational structure calculations: anharmonic frequencies of lithium fluoride clusters (LiF)n, n = 2-10

The vibrational spectra of a series of small lithium fluoride clusters, i.e. (LiF)n, n = 2-10, were studied by vibrational configuration interaction (VCI) calculations relying on potential energy surfaces including three-mode coupling terms and being obtained from explicitly correlated local coupled cluster calculations. Due to the account for anharmonicity effects, the simulated spectra allow for a direct comparison with experimental data and may thus help to identify clusters in the experiments. Even structurally closely related clusters can clearly be distinguished by infrared spectroscopy. The largest system in this study required more than 1000 basis functions in the electronic structure calculations and more than 10(7) configurations in the vibrational structure calculations and became computationally feasible only due to a combination of different approximations and highly parallelized algorithms.

Bracketing subtle conformational energy differences between self-solvated and stretched trifluoropropanol

The interconversion landscape between stretched and folded trifluoropropanol conformations is characterized by Raman jet spectroscopy and theory.

Anharmonicities and coherent vibrational dynamics of phosphate ions in bulk H2O

2D IR spectroscopy reveals Fermi resonances and long lived quantum beats for phosphate ions in water.

Constrained molecular vibration-rotation Hamiltonians: contravariant metric tensor

Here, I present a practical recipe for obtaining contravariant vibration-rotation metric tensors, and thus the kinetic energy operators, when some degrees of freedom are constrained rigidly. An element of the contravariant metric tensor is obtained as a sum of dot products of contravariant measuring vectors, which are obtained from their unconstrained counterparts by adding a frozen mode correction. The present method applies in principle for any choice of shape coordinates and a body-frame for which the contravariant measuring vectors can be evaluated. In contrast to the existing methods, the present method does not involve evaluation of covariant metric tensors, matrix inversions, chain rules of derivation, or numerical differentiation. It is applied in the sequel paper [L. Partanen, J. Pesonen, E. Sjöholm, and L. Halonen, J. Chem. Phys. 139, 144311 (2013)] to study the effects of several different approximations to the kinetic energy operator, when the two large-amplitude OH-torsional motions in H2SO4 are of interest.