An official website of the United States government
An approach to generate anharmonic potential energy surfaces for both linear and bent XY 2 -type molecules from their equilibrium geometries, Hessians, and total atomization energies alone is presented. Two key features of the potential energy surfaces are that (a) they reproduce the harmonic behavior around the equilibrium geometries exactly and (b) they have the correct limiting behavior with respect to total bond dissociation. The potentials are constructed from two diatomic potentials, for which both the Morse or Varshni potentials are tested, and a triatomic potential, for which modified forms of the Anderson- n potential are tested. Potential energy surfaces for several linear and bent molecules are constructed from ab initio data, and the third-order derivatives of these surfaces at their equilibrium geometries are compared to the results of finite difference computations. For bent molecules, the vibrational spectra predicted by vibrational configuration interaction calculations on these surfaces are compared to experiment. A modified version of the Anderson- n potential, in combination with the Varshni potential, is demonstrated to predict vibrational frequencies associated with bond angle bending an average of 20 cm −1 below the harmonic oscillator approximation and with a fourfold reduction in the root-mean-square deviation from experiment compared to the harmonic oscillator approximation.
more »
« less
Potential energy surfaces and dynamic properties via ab initio composite and density functional approaches
Abstract Vibrational spectroscopy enables critical insight into the structural and dynamic properties of molecules. Presently, the majority of theoretical approaches to spectroscopy employ wavefunction‐basedab initioor density functional methods that rely on the harmonic approximation. This approximation breaks down for large molecules with strongly anharmonic bonds or for molecules with large internuclear separations. An alternative to these methods involves generating molecular anharmonic potential energy surfaces (potentials) and using them to extrapolate the vibrational frequencies. This study examines the efficacy of density functional theory (DFT) and the correlation consistent Composite Approach (ccCA) in generating anharmonic frequencies from potentials of small main group molecules. Vibrational self‐consistent field Theory (VSCF) and post‐VSCF methods were used to calculate the fundamental frequencies of these molecules from their potentials. Functional choice, basis set selection, and mode‐coupling are also examined as factors in influencing accuracy. The absolute deviations for the calculated frequencies using potentials at the ccCA level of theory were lower than the potentials at the DFT level. With DFT resulting in bending modes that are better described than those of ccCA, a multilevel DFT:ccCA approach where DFT potentials are used for single vibrational mode potentials and ccCA is used for vibrational mode‐mode couplings can be utilized for larger polyatomic systems. The frequencies obtained with this multilevel approach using VCIPSI‐PT2 were closer to experimental frequencies than the scaled harmonic frequencies, indicating the success of utilizing post‐VSCF methods to generate more accurate representations of computed infrared spectra.
more »
« less
- Award ID(s):
- 2154526
- PAR ID:
- 10554470
- Publisher / Repository:
- Wiley
- Date Published:
- Journal Name:
- Journal of Computational Chemistry
- Volume:
- 45
- Issue:
- 16
- ISSN:
- 0192-8651
- Page Range / eLocation ID:
- 1352 to 1363
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Two-color, two-photon laser-induced fluorescence experiments were performed to probe the intermolecular interactions within the Ar + I2(E, vE = 0–3) potential energy surfaces. Spectra were recorded using the lowest-energy T-shaped level and an excited intermolecular vibrational level with bending excitation within the Ar + I2(B, vB = 23) potential as intermediate levels to guide the spectral assignments. Progressions of intermolecular stretching and bending levels bound within the Ar + I2(E, vE) potentials were identified, and their vibrational frequencies were determined. The harmonic frequency and anharmonic constant for the bending vibrational mode were determined to be ωe(b) ∼ 34.8 cm−1 and ωeχe(b) ∼ 0.3 cm−1. The frequency and anharmonic constant for the stretching mode were found to be the same as reported previously [V.V. Baturo, et al. Chem. Phys. Lett. 647 (2016) 161], ωe(s) = 37.2(1.1) cm−1 and ωeχe(s) = 1.8(2) cm−1.more » « less
-
Ion receptors are molecular hosts that bind ionic guests, often with great selectivity. The interplay of solvation and ion binding in anion host-guest complexes in solution governs the binding efficiency and selectivity of such ion receptors. To gain molecular-level insight into the intrinsic binding properties of octamethyl calix[4]pyrrole (omC4P) host molecules with halide guest ions, we performed cryogenic ion vibrational spectroscopy (CIVS) of omC4P in complexes with fluoride, chloride, and bromide ions. We interpret the spectra using density functional theory, describing the infrared spectra of these complexes with both harmonic and anharmonic second-order vibrational perturbation theory (VPT2) calculations. The NH stretching modes of the pyrrole moieties serve as sensitive probes of the ion binding properties, as their frequencies encode the ion-receptor interactions. While scaled harmonic spectra reproduce the experimental NH stretching modes of the chloride and bromide complexes in broad strokes, the high proton affinity of fluoride introduces strong anharmonic effects. As a result, the spectrum of F−·omC4P is not even qualitatively captured by harmonic calculations, but it is recovered very well by VPT2 calculations. In addition, the VPT2 calculations recover the intricate coupling of the NH stretching modes with overtones and combination bands of CH stretching and NH bending modes and with low-frequency vibrations of the omC4P macrocycle, which are apparent for all halide ion complexes investigated here. A comparison of the CIVS spectra with infrared spectra of solutions of the same ion-receptor complexes in d3-acetonitrile and d6-acetone shows how ion solvation changes the ion-receptor interactions for the different halide ions.more » « less
-
The terminal alkyne C≡C stretch has a large Raman scattering cross section in the “silent” region for biomolecules. This has led to many Raman tag and probe studies using this moiety to study biomolecular systems. A computational investigation of these systems is vital to aid in the interpretation of these results. In this work, we develop a method for computing terminal alkyne vibrational frequencies and isotropic transition polarizabilities that can easily and accurately be applied to any terminal alkyne molecule. We apply the discrete variable representation method to a localized version of the C≡C stretch normal mode. The errors of (1) vibrational localization to the terminal alkyne moiety, (2) anharmonic normal mode isolation, and (3) discretization of the Born–Oppenheimer potential energy surface are quantified and found to be generally small and cancel each other. This results in a method with low error compared to other anharmonic vibrational methods like second-order vibrational perturbation theory and to experiments. Several density functionals are tested using the method, and TPSS-D3, an inexpensive nonempirical density functional with dispersion corrections, is found to perform surprisingly well. Diffuse basis functions are found to be important for the accuracy of computed frequencies. Finally, the computation of vibrational properties like isotropic transition polarizabilities and the universality of the localized normal mode for terminal alkynes are demonstrated.more » « less
-
Terahertz (THz) spectroscopy has been put forth as a non-contact, analytical probe to characterize the intermolecular interactions of biologically active molecules, specifically as a way to understand, better develop, and use active pharmaceutical ingredients. An obstacle towards fully utilizing this technique as a probe is the need to couple features in the THz regions to specific vibrational modes and interactions. One solution is to use density functional theory (DFT) methods to assign specific vibrational modes to signals in the THz region, coupling atomistic insights to spectral features. Here, we use open source planewave DFT packages that employ ultrasoft pseudopotentials to assess the infrared (IR) response of organic compounds and complex co-crystal formulations in the solid state, with and without dispersion corrections. We compare our DFT computed lattice parameters and vibrational modes to experiment and comment on how to improve the agreement between theory and modeling to allow for THz spectroscopy to be used as an analytical probe in complex biologically relevant systems.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1002/jcc.27333
