Ab Initio Anharmonic Analysis of Vibrational Spectra of Uracil Using the Numerical-Analytic Implementation of Operator Van Vleck Perturbation Theory

Accuracy of Different Electronic Basis Set Families for Anharmonic Molecular Vibrations: A Comprehensive Benchmark Study

In this work, the accuracy and convergence of different electronic basis set families for the computation of anharmonic molecular vibrational spectroscopic calculations are benchmarked. A series of 39 different basis sets from different families following their hierarchy are assessed on VSCF and VSCF-PT2 algorithms with commonly used MP2 and DFT based B3LYP-D potentials for a set of molecular systems. Such an effort has been validated in a previous work ( J. Phys. Chem. A 2020, 124, 9203-9221) with split-valence basis sets for fundamentals and intensities. Here, fundamental transitions, vibrationally excited states, and intensities are compared with the experimental data to estimate the accuracy for a series of Jensen, Dunning, Calendar, Karlsruhe, and Sapporo basis set families. The convergence of basis sets are also compared with the large ANO basis set. Comprehensive statistical error analysis in terms of accuracy and precision was carried out to assess the performance of each basis set. It is observed that the improvement for the calculated harmonic and anharmonic values from the smaller basis sets to the medium (i.e., triple-ξ) is considerable. Beyond this, from medium to large basis sets, the convergence is slow and mostly posits nearly converged values. Basis sets with and without diffuse functions offer characteristically different accuracies and convergence patterns. Finally, recommendations are given on the choice of basis set chosen as black-box which can balance between accuracy and computational time, estimation of the errors, and their selections especially for large molecules.

Excited-State DMRG Made Simple with FEAST

We introduce DMRG[FEAST], a new method for optimizing excited-state many-body wave functions with the density matrix renormalization group (DMRG) algorithm. Our approach applies the FEAST algorithm, originally designed for large-scale diagonalization problems, to matrix product state wave functions. We show that DMRG[FEAST] enables the stable optimization of both low- and high-energy eigenstates, therefore overcoming the limitations of state-of-the-art excited-state DMRG algorithms. We demonstrate the reliability of DMRG[FEAST] by calculating anharmonic vibrational excitation energies of molecules with up to 30 fully coupled degrees of freedom.

Neural Network Potential Energy Surfaces for Small Molecules and Reactions

We review progress in neural network (NN)-based methods for the construction of interatomic potentials from discrete samples (such as ab initio energies) for applications in classical and quantum dynamics including reaction dynamics and computational spectroscopy. The main focus is on methods for building molecular potential energy surfaces (PES) in internal coordinates that explicitly include all many-body contributions, even though some of the methods we review limit the degree of coupling, due either to a desire to limit computational cost or to limited data. Explicit and direct treatment of all many-body contributions is only practical for sufficiently small molecules, which are therefore our primary focus. This includes small molecules on surfaces. We consider direct, single NN PES fitting as well as more complex methods that impose structure (such as a multibody representation) on the PES function, either through the architecture of one NN or by using multiple NNs. We show how NNs are effective in building representations with low-dimensional functions including dimensionality reduction. We consider NN-based approaches to build PESs in the sums-of-product form important for quantum dynamics, ways to treat symmetry, and issues related to sampling data distributions and the relation between PES errors and errors in observables. We highlight combinations of NNs with other ideas such as permutationally invariant polynomials or sums of environment-dependent atomic contributions, which have recently emerged as powerful tools for building highly accurate PESs for relatively large molecular and reactive systems.

Overtone spectroscopy of v(C=O) stretching vibration of hexafluoroacetone: Experimental and ab initio determination of peak positions, absolute intensities, and band shapes

Infrared overtone spectra of the ν(C = O) stretching vibration (ν₁) of gaseous hexafluoroacetone ((CF₃)₂C = O, HFA) were recorded in the spectral range of 7450-3300 cm^-1 with a resolution of 0.1 cm^-1. Experimental absolute IR intensities and vibrational band centers of the overtones 2ν₁, 3ν₁, 4ν₁ of HFA were measured and compared with their ab initio counterparts, calculated by the second-order canonical vibrational perturbation theory (CVPT2). A hybrid potential energy surface (PES) was evaluated using the quantum-mechanical models MP2/cc-pVQZ for harmonic and MP2/cc-pVTZ for anharmonic parts. Cubic surfaces of dipole moment components were determined using the MP2/cc-pVTZ model. The predicted IR intensities for the first and second overtones reproduced the experimental values with a discrepancy of 6% and 9%, respectively. A weak Fermi resonance between the first overtone 2ν₁ and combination tone ν₁ + ν₂ + ν₈ was predicted. The appropriate model was employed for simulating the bands studied using: (i) the asymmetric-top vibration-rotational structure of the ν₁ mode, (ii) the inhomogeneous band structure due to contributions of the vibrational states populated at the room temperature, and (iii) the homogeneous broadening of each vibration-rotation transition due to the intramolecular vibrational redistribution (IVR). The rotational and vibrational anharmonic constants were taken from the ab initio calculations, whereas the IVR data were obtained from the experimental data of Chekalin et al. (JPCA 2014, 118:955) with a time resolution of ≈100 fs. A good level of agreement of the predicted shapes of bands studied with experiment is achieved.

Plumbing Potentials for Molecules with Up To Tens of Atoms: How to Find Saddle Points and Fix Leaky Holes

Potential energy surfaces (PESs) play an indispensable role in molecular dynamics but are notoriously difficult to flesh out properly in large-dimensional spaces. In particular, the undetected presence of PES holes, i.e., unphysical saddle points beyond which the potential energy drops arbitrarily, can have devastating effects on both classical and quantum dynamics calculations. In this study, the Crystal algorithm is developed as a tool for efficiently and accurately finding PES holes, as well as legitimate saddle points, even in very large-dimensional configuration spaces. The approach is applied to three large-dimensional PESs for molecular systems of current interest: uracil, naphthalene, and formic acid dimer. Low-lying PES holes are discovered and located for the first two systems-including naphthalene, for which no holes were previously suspected, to the best of our knowledge. Likewise, the double-well, double-proton-transfer isomerization saddle point for formic acid dimer is also located.

Vibrations of the guanine-cytosine pair in chloroform: an anharmonic computational study

We compute at the anharmonic level the vibrational spectra of the Watson-Crick dimer formed by guanosine (G) and cytidine (C) in chloroform, together with those of G, C and the most populated GG dimer. The spectra for deuterated and partially deuterated GC are also computed. We use DFT calculations, with B3LYP and CAM-B3LYP as reference functionals. Solvent effects from chloroform are included via the Polarizable Continuum Model (PCM), and by performing tests on models including up two chloroform molecules. Both B3LYP and CAM-B3LYP calculations reproduce the shape of the experimental spectra well in the fingerprint region (1500-1700 cm-1) and in the N-H stretching region (2800-3600 cm-1), with B3LYP providing better quantitative agreement with experiments. According to our calculations, the N-H amido streching mode of G falls at ∼2900 cm-1, while the N-H amino of G and C falls at ∼3100 cm-1 when hydrogen-bonded, or ∼3500 cm-1 when free. Overtone and combination bands strongly contribute to the absorption band at ∼3300 cm-1. Inclusion of bulk solvent effects significantly increases the accuracy of the computed spectra, while solute-solvent interactions have a smaller, though still noticeable, effect. Some key aspects of the anharmonic treatment of strongly vibrationally coupled supermolecular systems and the related methodological issues are also discussed.

Simulated NIR spectra as sensitive markers of the structure and interactions in nucleobases

Near-infrared (near-IR; NIR) spectroscopy is continuously advancing in biophysical and biochemical fields of investigation. For instance, recent progresses in NIR hyperspectral imaging of biological systems may be noted. However, interpretation of NIR bands for biological samples is difficult and creates a considerable barrier in exploring the full potential of NIR spectroscopy in bioscience. For this reason, we carried out a systematic study of NIR spectra of adenine, cytosine, guanine, and thymine in polycrystalline state. Interpretation of NIR spectra of these nucleobases was supported by anharmonic vibrational analysis using Deperturbed Vibrational Second-Order Perturbation Theory (DVPT2). A number of molecular models of nucleobases was applied to study the effect of the inter-molecular interactions on the NIR spectra. The accuracy of simulated NIR spectra appears to depend on the intra-layer interactions; in contrast, the inter-layer interactions are less influential. The best results were achieved by combining the simulated spectra of monomers and dimers. It is of particular note that in-plane deformation bands are far more populated than out-of-plane ones and the importance of ring modes is relatively small. This trend is in contrast to that observed in mid-IR region. As shown, the local, short-range chemical neighborhood of nucleobase molecules influence their NIR spectra more considerably. This suggests that NIR spectra are more sensitive probe of the nucleobase pairing than mid-IR ones. The obtained results allow, for the first time, to construct a frequency correlation table for NIR spectra of purines and pyrimidines.

Vibrational investigation of nucleobases by means of divide and conquer semiclassical dynamics

In this work, we report a computational study of the vibrational features of four different nucleobases employing the divide-and-conquer semiclassical initial value representation molecular dynamics method. Calculations are performed on uracil, cytosine, thymine, and adenine. Results show that the overall accuracy with respect to experiments is within 20 wavenumbers, regardless of the dimensionality of the nucleobase. Vibrational estimates are accurate even in the complex case of cytosine, where two relevant conformers are taken into account. These results are promising in the perspective of future studies on more complex systems, such as nucleotides or nucleobase pairs.

Energetics and Vibrational Signatures of Nucleobase Argyrophilic Interactions

This study investigates the interactions of both purine (adenine and guanine) and pyrimidine (cytosine, thymine, and uracil) nucleobases with a pair of silver atoms (Ag2). Full geometry optimizations were performed on several structures of each nucleobase/Ag2 complex and the corresponding isolated monomers using the M06-2X density functional with a correlation consistent triple-ζ basis set augmented with diffuse functions on all atoms and a relativistic pseudopotential on Ag (aug-cc-pVTZ for H, C, N, and O and aug-cc-pVTZ-PP for Ag; denoted aVTZ). Harmonic vibrational frequency computations indicate that each optimized structure corresponds to a minimum on the M06-2X/aVTZ potential energy surface. Relative electronic energies for interactions between Ag2 and each nucleobase were compared to elucidate energetic differences between isomers. Further analysis of the changes in vibrational frequencies, infrared intensities, and Raman scattering activities reveals how different Ag2 binding sites might be differentiated spectroscopically. These results provide molecular-level insight into the interactions between nucleobases and silver, which may lead to better understanding and interpretation of surface-enhanced Raman spectroscopy experiments on nucleobases and related systems.

Application of classical simulations for the computation of vibrational properties of free molecules

In this study, we investigate the ability of classical molecular dynamics (MD) and Monte-Carlo (MC) simulations for modeling the intramolecular vibrational motion. These simulations were used to compute thermally-averaged geometrical structures and infrared vibrational intensities for a benchmark set previously studied by gas electron diffraction (GED): CS₂, benzene, chloromethylthiocyanate, pyrazinamide and 9,12-I₂-1,2-closo-C₂B₁₀H₁₀. The MD sampling of NVT ensembles was performed using chains of Nose-Hoover thermostats (NH) as well as the generalized Langevin equation thermostat (GLE). The performance of the theoretical models based on the classical MD and MC simulations was compared with the experimental data and also with the alternative computational techniques: a conventional approach based on the Taylor expansion of potential energy surface, path-integral MD and MD with quantum-thermal bath (QTB) based on the generalized Langevin equation (GLE). A straightforward application of the classical simulations resulted, as expected, in poor accuracy of the calculated observables due to the complete neglect of quantum effects. However, the introduction of a posteriori quantum corrections significantly improved the situation. The application of these corrections for MD simulations of the systems with large-amplitude motions was demonstrated for chloromethylthiocyanate. The comparison of the theoretical vibrational spectra has revealed that the GLE thermostat used in this work is not applicable for this purpose. On the other hand, the NH chains yielded reasonably good results.

Anharmonic vibrational analysis of s-trans and s-cis conformers of acryloyl fluoride using numerical-analytic Van Vleck operator perturbation theory

A new gas-phase infrared (IR) spectrum of acryloyl fluoride (ACRF, CH₂CHCFO) with a resolution of 0.1cm^-1 in the range 4000-450cm^-1 was measured. Theoretical ab initio molecular structures, full quartic potential energy surfaces (PES), and cubic surfaces of dipole moments and polarizability tensor components (electro-optical properties, EOP) of the s-trans and s-cis conformers of the ACRF were calculated by the second-order Møller-Plesset electronic perturbation theory with a correlation consistent Dunning triple-ζ basis set. The numerical-analytic implementation of the second-order operator canonical Van Vleck perturbation theory was employed for predicting anharmonic IR and Raman scattering (RS) spectra of ACRF. To improve the anharmonic predictions, harmonic frequencies were replaced by their counterparts evaluated with the higher-level CCSD(T)/cc-pVTZ model, to form a "hybrid" PES. The original operator representation of the Hamiltonian is analytically reduced to a quasi-diagonal form, integrated in the harmonic oscillator basis and diagonalized to account for strong resonance couplings. Double canonical transformations of EOP expansions enabled prediction of integral intensities of both fundamental and multi-quanta transitions in IR/RS spectra. Enhanced band shape analysis reinforced the assignments. A thorough interpretation of the new IR experimental spectra and existing matrix-isolation literature data for the mixture of two conformers of ACRF was accomplished, and a number of assignments clarified.

Anharmonic Effects on Vibrational Spectra Intensities: Infrared, Raman, Vibrational Circular Dichroism, and Raman Optical Activity

The aim of this paper is 2-fold. First, we want to report the extension of our virtual multifrequency spectrometer (VMS) to anharmonic intensities for Raman optical activity (ROA) with the full inclusion of first- and second-order resonances for both frequencies and intensities in the framework of the generalized second-order vibrational perturbation theory (GVPT2) for all kinds of vibrational spectroscopies. Then, from a more general point of view, we want to present and validate the performance of VMS for the parallel analysis of different vibrational spectra for medium-sized molecules (IR, Raman, VCD, ROA) including both mechanical and electric/magnetic anharmonicity. For the well-known methyloxirane benchmark, careful selection of density functional, basis set, and resonance thresholds permitted us to reach qualitative and quantitative agreement between experimental and computed band positions and shapes. Next, the whole series of halogenated azetidinones is analyzed, showing that it is now possible to interpret different spectra in terms of mass, electronegativity, polarizability, and hindrance variation between closely related substituents, chiral spectroscopies being particular effective in this connection.