Calculation of the vibrational wave function of polyatomic molecules

Vibrational Corrections to NMR Spin-Spin Coupling Constants from Relativistic Four-Component DFT Calculations

Zero-point vibrational (ZPV) corrections to the nuclear spin–spin coupling constants have been calculated using four-component Dirac–Kohn–Sham DFT for H₂X (where X = O, S, Se, Te, Po), XH₃ (where X = N, P, As, Sb, Bi), and XH₄ (where X = C, Si, Ge, Sn, and Pb) molecules and for HC≡CPbH₃. The main goal was to study the influence of relativistic effects on the ZPV corrections and thus results calculated at relativistic and nonrelativistic approaches have been compared. The effects of relativity become notable for the ZPV corrections to the spin–spin coupling constants for compounds with lighter elements (selenium and germanium) than for the spin–spin coupling constants themselves. In the case of molecules containing heavier atoms, for instance BiH₃ and PbH₄, relativistic effects play a crucial role on the results and approximating ZPV corrections by the nonrelativistic results may lead to larger errors than omitting ZPV corrections altogether.

Quantum Chemical Approaches to the Calculation of NMR Parameters: From Fundamentals to Recent Advances

Quantum chemical methods for the calculation of indirect NMR spin–spin coupling constants and chemical shifts are always in progress. They never stay the same due to permanently developing computational facilities, which open new perspectives and create new challenges every now and then. This review starts from the fundamentals of the nonrelativistic and relativistic theory of nuclear magnetic resonance parameters, and gradually moves towards the discussion of the most popular common and newly developed methodologies for quantum chemical modeling of NMR spectra.

A new double-reference correction scheme for accurate and efficient computation of NMR chemical shieldings

Our novel correction procedure yields high-accuracy DFT predictions of absolute NMR shieldings and enables outliers due to relativistic effects or manifestly inadequate modelling of electron correlation to be easily and unambiguously identified.

Performance of the DLPNO-CCSD and recent DFT methods in the calculation of isotropic and dipolar contributions to 14N hyperfine coupling constants of nitroxide radicals

In the present study, the performance of a set of density functionals: BP86, PBE, OLYP, BEEF, PBEpow, TPSS, SCAN, PBEGXPBE, M06L, MN15L, B3LYP, PBE0, mPW1PW, B97, BHandHLYP, mPW1PW, B98, TPSS0, PBE1KCIS, SCAN0, M06, M06-2X, MN15, CAM-B3LYP, ωB97x, B2PLYP, and the B3LYP/N07D and PBE/N07D schemes in the calculation of the ¹⁴N anisotropic hyperfine coupling (HFC) constants of a set of 23 nitroxide radicals is evaluated. The results are compared with those obtained with the DLPNO-CCSD method and experimental HFC values. Harmonic contribution to the ¹⁴N HFC vibrational correction was calculated at the revPBE0/def2-TZVPP level and included in the evaluation. With the vibrational correction, the DLPNO-CCSD method yielded HFC values in good agreement with the experiment (mean absolute deviation (MAD) = 0.3 G for the dipole-dipole contribution and MAD = 0.8 G for the contact coupling contribution). The best DFT results are obtained using the M06 functional with MAD = 0.2 G for the dipole-dipole contribution and MAD = 0.7 G for the contact coupling contribution. In general, vibrational correction significantly improved most DFT functionals' performance but did not change its overall ranking.

Representing molecular ground and excited vibrational eigenstates with nuclear densities obtained from semiclassical initial value representation molecular dynamics

We present in detail and validate an effective Monte Carlo approach for the calculation of the nuclear vibrational densities via integration of molecular eigenfunctions that we have preliminary employed to calculate the densities of the ground and the excited OH stretch vibrational states in the protonated glycine molecule [Aieta et al., Nat Commun 11, 4348 (2020)]. Here, we first validate and discuss in detail the features of the method on a benchmark water molecule. Then, we apply it to calculate on-the-fly the ab initio anharmonic nuclear densities in the correspondence of the fundamental transitions of NH and CH stretches in protonated glycine. We show how we can gain both qualitative and quantitative physical insight by inspection of different one-nucleus densities and assign a character to spectroscopic absorption peaks using the expansion of vibrational states in terms of harmonic basis functions. The visualization of the nuclear vibrations in a purely quantum picture allows us to observe and quantify the effects of anharmonicity on the molecular structure, also to exploit the effect of IR excitations on specific bonds or functional groups, beyond the harmonic approximation. We also calculate the quantum probability distribution of bond lengths, angles, and dihedrals of the molecule. Notably, we observe how in the case of one type of fundamental NH stretching, the typical harmonic nodal pattern is absent in the anharmonic distribution.

Anharmonic quantum nuclear densities from full dimensional vibrational eigenfunctions with application to protonated glycine

The interpretation of molecular vibrational spectroscopic signals in terms of atomic motion is essential to understand molecular mechanisms and for chemical characterization. The signals are usually assigned after harmonic normal mode analysis, even if molecular vibrations are known to be anharmonic. Here we obtain the quantum anharmonic vibrational eigenfunctions of the 11-atom protonated glycine molecule and we calculate the density distribution of its nuclei and its geometry parameters, for both the ground and the O-H stretch excited states, using our semiclassical method based on ab initio molecular dynamics trajectories. Our quantum mechanical results describe a molecule elongated and more flexible with respect to what previously thought. More importantly, our method is able to assign each spectral peak in vibrational spectroscopy by showing quantitatively how normal modes involving different functional groups cooperate to originate that spectroscopic signal. The method will possibly allow for a better rationalization of experimental spectroscopy.

The Structure of the "Vibration Hole" around an Isotopic Substitution-Implications for the Calculation of Nuclear Magnetic Resonance (NMR) Isotopic Shifts

Calculations of nuclear magnetic resonance (NMR) isotopic shifts often rest on the unverified assumption that the “vibration hole”, that is, the change of the vibration motif upon an isotopic substitution, is strongly localized around the substitution site. Using our recently developed difference-dedicated (DD) second-order vibrational perturbation theory (VPT2) method, we test this assumption for a variety of molecules. The vibration hole turns out to be well localized in many cases but not in the interesting case where the H/D substitution site is involved in an intra-molecular hydrogen bond. For a series of salicylaldehyde derivatives recently studied by Hansen and co-workers (Molecules2019, 24, 4533), the vibrational hole was found to stretch over the whole hydrogen-bond moiety, including the bonds to the neighbouring C atoms, and to be sensitive to substituent effects. We discuss consequences of this finding for the accurate calculation of NMR isotopic shifts and point out directions for the further improvement of our DD-VPT2 method.

Efficient calculation of NMR isotopic shifts: Difference-dedicated vibrational perturbation theory

We present difference-dedicated second-order vibrational perturbation theory (VPT2) as an efficient method for the computation of nuclear magnetic resonance (NMR) isotopic shifts, which reflect the geometry dependence of the NMR property in combination with different vibration patterns of two isotopologues. Conventional calculations of isotopic shifts, e.g., by standard VPT2, require scanning the geometry dependence over the whole molecule, which becomes expensive rapidly as the molecule size increases. In DD-VPT2, this scan can be restricted to a small region around the substitution site. At the heart of DD-VPT2 is a set of localized vibration modes common to the two isotopologues and designed such that the difference between the vibration patterns is caught by a small subset of them (usually fewer than 10). We tested the DD-VPT2 method for a series of molecules with increasing size and found that this method provides results with the same quality as VPT2 and in good agreement with the experiment, with computational savings up to 95% and less numerical instabilities. The method is easy to automatize and straightforward to generalize to other molecular properties.

Towards reliable references for electron paramagnetic resonance parameters based on quantum chemistry: the case of verdazyl radicals

We present an efficient and accurate computational procedure to calculate properties measurable by EPR spectroscopy. We simulate a molecular dynamics (MD) trajectory by employing the quantum mechanically derived force field (QMDFF) [S. Grimme, J. Chem. Theory Comput., 2014, 10, 4497] and sample the trajectory at different time steps. For each snapshot EPR properties are calculated with a hybrid density functional theory (DFT) method. EPR spectra are simulated based on the averaged results. We applied the strategy to a number of previously published and novel verdazyl radicals, for which we recorded EPR spectra. The resulting simulated spectra are compatible with experiment already before employing an additional fitting step, in contrast to those from single point electronic-structure calculations. After the refinement, the experimental data are excellently reproduced, and the fitted EPR parameters do not deviate much from the calculated ones. This provides confidence in ascribing a direct physical meaning to the refined data in terms of experimental EPR parameters rather than merely considering them as mathematical fit parameters. We also find that couplings to hydrogen nuclei have a significant influence on the spectra of verdazyl radicals.

Optical Rotation Calculations for a Set of Pyrrole Compounds

Optical rotation of 14 molecules containing the pyrrole group is calculated by employing both time-dependent density functional theory (TDDFT) with the CAM-B3LYP functional and the second-order approximate coupled-cluster singles and doubles (CC2) method. All optical rotations have been provided using the aug-cc-pVDZ basis set at λ = 589 nm. The two methods predict similar results for both sign and magnitude for the optical rotation of all molecules. The obtained signs are consistent with experiments as well, although several conformers for four molecules needed to be studied to reproduce the experimental sign. We have also calculated excitation energies and rotatory strengths for the six lowest lying electronic transitions for several conformers of the two smallest molecules and found that each rotatory strength has various contributions for each conformer which can cause different optical rotations for different conformers of a molecule. Our results illustrate that both methods are able to reproduce the experimental optical rotations, and that the CAM-B3LYP functional, the least computationally expensive method used here, is an applicable and reliable method to predict the optical rotation for these molecules in line with previous studies.

Rotational constants and structure of para-difluorobenzene determined by femtosecond Raman coherence spectroscopy: A new transient type

Femtosecond Raman rotational coherence spectroscopy (RCS) detected by degenerate four-wave mixing is a background-free method that allows to determine accurate gas-phase rotational constants of non-polar molecules. Raman RCS has so far mostly been applied to the regular coherence patterns of symmetric-top molecules, while its application to nonpolar asymmetric tops has been hampered by the large number of RCS transient types, the resulting variability of the RCS patterns, and the 10(3)-10(4) times larger computational effort to simulate and fit rotational Raman RCS transients. We present the rotational Raman RCS spectra of the nonpolar asymmetric top 1,4-difluorobenzene (para-difluorobenzene, p-DFB) measured in a pulsed Ar supersonic jet and in a gas cell over delay times up to ∼2.5 ns. p-DFB exhibits rotational Raman transitions with ΔJ = 0, 1, 2 and ΔK = 0, 2, leading to the observation of J -, K -, A -, and C-type transients, as well as a novel transient (S-type) that has not been characterized so far. The jet and gas cell RCS measurements were fully analyzed and yield the ground-state (v = 0) rotational constants A0 = 5637.68(20) MHz, B0 = 1428.23(37) MHz, and C0 = 1138.90(48) MHz (1σ uncertainties). Combining the A0, B0, and C0 constants with coupled-cluster with single-, double- and perturbatively corrected triple-excitation calculations using large basis sets allows to determine the semi-experimental equilibrium bond lengths re(C1-C2) = 1.3849(4) Å, re(C2-C3) = 1.3917(4) Å, re(C-F) = 1.3422(3) Å, and re(C2-H2) = 1.0791(5) Å.

Vibrationally averaged dipole moments of methane and benzene isotopologues

DFT-B3LYP post-Born-Oppenheimer (finite-nuclear-mass-correction (FNMC)) calculations of vibrationally averaged isotopic dipole moments of methane and benzene, which compare well with experimental values, are reported. For methane, in addition to the principal vibrational contribution to the molecular asymmetry, FNMC accounts for the surprisingly large Born-Oppenheimer error of about 34% to the dipole moments. This unexpected result is explained in terms of concurrent electronic and vibrational contributions. The calculated dipole moment of C6H3D3 is about twice as large as the measured dipole moment of C6H5D. Computational progress is advanced concerning applications to larger systems and the choice of appropriate basis sets. The simpler procedure of performing vibrational averaging on the Born-Oppenheimer level and then adding the FNMC contribution evaluated at the equilibrium distance is shown to be appropriate. Also, the basis set choice is made by heuristic analysis of the physical behavior of the systems, instead of by comparison with experiments.

Comparing Vibrationally Averaged Nuclear Shielding Constants by Quantum Diffusion Monte Carlo and Second-Order Perturbation Theory

Using the method of modified Shepard's interpolation to construct potential energy surfaces of the H2O, O3, and HCOOH molecules, we compute vibrationally averaged isotropic nuclear shielding constants ⟨σ⟩ of the three molecules via quantum diffusion Monte Carlo (QDMC). The QDMC results are compared to that of second-order perturbation theory (PT), to see if second-order PT is adequate for obtaining accurate values of nuclear shielding constants of molecules with large amplitude motions. ⟨σ⟩ computed by the two approaches differ for the hydrogens and carbonyl oxygen of HCOOH, suggesting that for certain molecules such as HCOOH where big displacements away from equilibrium happen (internal OH rotation), ⟨σ⟩ of experimental quality may only be obtainable with the use of more sophisticated and accurate methods, such as quantum diffusion Monte Carlo. The approach of modified Shepard's interpolation is also extended to construct shielding constants σ surfaces of the three molecules. By using a σ surface with the equilibrium geometry as a single data point to compute isotropic nuclear shielding constants for each descendant in the QDMC ensemble representing the ground state wave function, we reproduce the results obtained through ab initio computed σ to within statistical noise. Development of such an approach could thereby alleviate the need for any future costly ab initio σ calculations.

Rovibrational and Temperature Effects in Theoretical Studies of NMR Parameters

The demand for high precision calculations of NMR shieldings (or their related values, chemical shifts δ) and spin–spin coupling constants facilitating and supporting detailed interpretations of NMR spectra increases hand in hand with the development of computational techniques and hardware resources. Highly sophisticated calculations including even relativistic effects are nowadays possible for these properties. However, NMR parameters depend not only on molecular structure and environment but also on molecular flexibility and temperature and the apparent success of theoretical predictions for molecular equilibrium geometries creates a demand for zero-point vibrational and temperature corrections. In this chapter we describe briefly the theory behind rovibrational corrections and review then some important contributions to this field.

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.

Ro-vibrational averaging of the isotropic hyperfine coupling constant for the methyl radical

We present the first variational calculation of the isotropic hyperfine coupling constant of the carbon-13 atom in the CH3 radical for temperatures T = 0, 96, and 300 K. It is based on a newly calculated high level ab initio potential energy surface and hyperfine coupling constant surface of CH3 in the ground electronic state. The ro-vibrational energy levels, expectation values for the coupling constant, and its temperature dependence were calculated variationally by using the methods implemented in the computer program TROVE. Vibrational energies and vibrational and temperature effects for coupling constant are found to be in very good agreement with the available experimental data. We found, in agreement with previous studies, that the vibrational effects constitute about 44% of the constant's equilibrium value, originating mainly from the large amplitude out-of-plane bending motion and that the temperature effects play a minor role.

New Experimental NMR Shielding Scales Mapped Relativistically from NSR: Theory and Application

The recently proposed relativistic mapping between nuclear magnetic resonance (NMR) shielding and nuclear spin-rotation (NSR) coupling tensors [J. Chem. Phys. 2013, 138, 134104] is employed to establish new experimental (more precisely, experimentally derived) absolute shielding constants for H and X in HX (X = F, Cl, Br, and I). The results are much more accurate than the old "experimental" values that were based on the well-known nonrelativistic mapping. The relativistic mapping is very robust in the sense that it is rather insensitive to the quality of one-particle basis sets and the treatment of electron correlation. Relativistic effects in the NSR coupling constants are also elucidated in depth.

The Dalton quantum chemistry program system

Dalton is a powerful general-purpose program system for the study of molecular electronic structure at the Hartree-Fock, Kohn-Sham, multiconfigurational self-consistent-field, Møller-Plesset, configuration-interaction, and coupled-cluster levels of theory. Apart from the total energy, a wide variety of molecular properties may be calculated using these electronic-structure models. Molecular gradients and Hessians are available for geometry optimizations, molecular dynamics, and vibrational studies, whereas magnetic resonance and optical activity can be studied in a gauge-origin-invariant manner. Frequency-dependent molecular properties can be calculated using linear, quadratic, and cubic response theory. A large number of singlet and triplet perturbation operators are available for the study of one-, two-, and three-photon processes. Environmental effects may be included using various dielectric-medium and quantum-mechanics/molecular-mechanics models. Large molecules may be studied using linear-scaling and massively parallel algorithms. Dalton is distributed at no cost from http://www.daltonprogram.org for a number of UNIX platforms.

Role of zero-point vibrational corrections to carbon hyperfine coupling constants in organic π radicals

By analyzing a set of organic π radicals, we demonstrate that zero-point vibrational corrections give significant contributions to carbon hyperfine coupling constants, in one case even inducing a sign reversal for the coupling constant. We discuss the implications of these findings for the computational analysis of electron paramagnetic spectra based on hyperfine coupling constants evaluated at the equilibrium geometry of radicals. In particular, we note that a dynamical description that involves the nuclear motion is in many cases necessary in order to achieve a semi-quantitatively predictive theory for carbon hyperfine coupling constants. In addition, we discuss the implications of the strong dependence of the carbon hyperfine coupling constants on the zero-point vibrational corrections for the selection of exchange-correlation functionals in density functional theory studies of these constants.

A second-order perturbation theory route to vibrational averages and transition properties of molecules: general formulation and application to infrared and vibrational circular dichroism spectroscopies

A general formulation to compute anharmonic vibrational averages and transition properties at the second-order of perturbation theory is derived from the Rayleigh-Schrödinger development. This approach is intended to be applicable to any property expanded as a Taylor series up to the third order with respect to normal coordinates or their associated momenta. The equations are straightforward to implement and can be easily adapted to various properties, as illustrated for the case of electric and magnetic dipole moments. From those, infrared and vibrational circular dichroism spectra can be readily obtained. This fully automatic procedure has been applied to several chiral molecules of small-to-medium sizes and compared to the standard double harmonic approximation and to experimental data.

Zero-point vibrational corrections to isotropic hyperfine coupling constants in polyatomic molecules

The present work addresses isotropic hyperfine coupling constants in polyatomic systems with a particular emphasis on a largely neglected, but a posteriori significant, effect, namely zero-point vibrational corrections. Using the density functional restricted-unrestricted approach, the zero-point vibrational corrections are evaluated for the allyl radical and four of its derivatives. In addition for establishing the numerical size of the zero-point vibrational corrections to the isotropic hyperfine coupling constants, we present simple guidelines useful for identifying hydrogens for which such corrections are significant. Based on our findings, we critically re-examine the computational procedures used for the determination of hyperfine coupling constants in general as well as the practice of using experimental hyperfine coupling constants as reference data when benchmarking and optimizing exchange-correlation functionals and basis sets for such calculations.

Calculation of NMR chemical shifts in organic solids: accounting for motional effects

NMR chemical shifts were calculated from first principles for well defined crystalline organic solids. These density functional theory calculations were carried out within the plane-wave pseudopotential framework, in which truly extended systems are implicitly considered. The influence of motional effects was assessed by averaging over vibrational modes or over snapshots taken from ab initio molecular dynamics simulations. It is observed that the zero-point correction to chemical shifts can be significant, and that thermal effects are particularly noticeable for shielding anisotropies and for a temperature-dependent chemical shift. This study provides insight into the development of highly accurate first principles calculations of chemical shifts in solids, highlighting the role of motional effects on well defined systems.

Dinuclear fluoro-peroxovanadium(v) complexes with symmetric and asymmetric peroxo bridges: syntheses, structures and DFT studies

Two new dinuclear fluoro peroxovanadium(v) complexes, Cs3[V2O2(O2)4F] x H2O (1) and Cs3[V2O2(O2)3F3] x 2HF x H2O (2), were prepared and characterized by elemental analysis, IR spectroscopy, thermal analysis and X-ray crystallography. While the anion in possesses an asymmetric structure with a micro-eta1:eta2 bridging peroxo group, the [V2O2(O2)3F3]3- ion in exhibits a symmetrical structure with a unique mu-fluoro and micro-eta2:eta2 peroxo double bridge. The X-ray structure data were compared with equilibrium and vibrationally-averaged (effective) DFT calculated geometries. The decomposition reactions of and in aqueous solution were studied by 51V NMR spectroscopy. The calculations of vibrationally averaged NMR chemical shifts (DFT-GIAO) were used to support the empirical assignment of NMR signals and afforded excellent agreement with experimental values for the studied peroxovanadium species. The ESI mass spectra of the prepared compounds are in accordance with the assignment of NMR spectra and with DFT study.

Molecular structure of the octatetranyl anion, C8H-: a computational study

The equilibrium molecular structure of the octatetranyl anion, C8H(-), which has been recently detected in two astronomical environments, is investigated with the aid of both ab initio post-Hartree-Fock and density functional theory (DFT) calculations. The model chemistry adopted in this study was selected after a series of benchmark calculations performed on molecular acetylene for which accurate gas-phase structural data are available. Geometry optimizations performed at the CCSD/6-311+G(2d,p), QCISD/6-311+G(2d,p), and MP4(SDQ)/6-311+G(2d,p) levels of theory yield for C8H(-) an interesting polyyne-type structure that defies the chemical formula displaying a simple alternation of triple and single carbon-carbon bonds, [:C[triple bond]C-C[triple bond]C-C[triple bond]C-C[triple bond]CH](1-). In the optimized geometry of C8H(-), as one proceeds from the naked carbon atom on one side of the chain to the CH unit on the opposite side of the chain, the short (formally triple) carbon-carbon bonds decrease in length from 1.255 to 1.213 A whereas the long (formally single) carbon-carbon bonds increase (albeit only slightly) in length from 1.362 to 1.378 A (CCSD results). In striking contrast, both MP2 and DFT (B3LYP and PBE0) calculations fail in reproducing the pattern of the carbon-carbon bond lengths obtained with the CCSD, QCISD, and MP4 methods. The structures of three shorter n-even chains, C(n)H(-) (n = 2, 4, and 6), along with those of four n-odd compounds (n = 3, 5, 7, and 9) are also investigated at the CCSD/6-311+G(2d,p) level of theory.