Anharmonic Vibrational Spectroscopy Calculations for Proton-Bound Amino Acid Dimers

Dual Basis Approach for Ab Initio Anharmonic Calculations of Vibrational Spectroscopy: Application to Microsolvated Biomolecules

A dual electronic basis set approach is introduced for more efficient but accurate calculations of the anharmonic vibrational spectra in the framework of the vibrational self-consistent field (VSCF) theory. In this approach, an accurate basis set is used to compute the vibrational spectra at the harmonic level. The results are used to scale the potential surface from a more modest but much more efficient basis set. The scaling is such that at the harmonic level the new, scaled potential agrees with one of the accurate basis sets. The approach is tested in the application of the microsolvated, protected amino acid Ac-Phe-OMe, using the scaled anharmonic hybrid potential in the VSCF and VSCF-PT2 algorithms. The hybrid potential method yields results that are in good accord with the experiment and very close to those obtained in calculations with the high-level, very costly potential from the large basis set. At the same time, the hybrid potential calculations are considerably less expensive. The results of the hybrid calculations are much more accurate than those computed from the potential surface corresponding to the modest basis set. The results are very encouraging for using the hybrid potential method for inexpensive yet sufficiently accurate anharmonic calculations for the spectra of large biomolecules.

Structure of Proton-Bound Methionine and Tryptophan Dimers in the Gas Phase Investigated with IRMPD Spectroscopy and Quantum Chemical Calculations

The structures of three proton-bound dimers (Met₂H⁺, MetTrpH⁺, and Trp₂H⁺) are investigated in the gas phase with infrared multiple photon disassociation (IRMPD) spectroscopy in combination with quantum chemical calculations. Their IRMPD spectra in the range of 600–1850 cm^–1 are obtained experimentally using an FT-ICR mass spectrometer and the CLIO free electron laser as an IR light source. The most abundant conformers are elucidated by comparing the IRMPD spectra with harmonic frequencies obtained at the B3LYP-GD3BJ/6-311++G** level of theory. Discrepancies between the experimental and theoretical data in the region of 1500–1700 cm^–1 are attributed to the anharmonicity of the amino bending modes. We confirm the result of a previous IRMPD study that the structure of gas-phase Trp₂H⁺ is charge-solvated but find that there are more stable structures than originally reported (Feng, R.; Yin, H.; Kong, X. Rapid Commun. Mass Spectrom.2016, 30, 24–28). In addition, gas-phase Met₂H⁺ and MetTrpH⁺ have been revealed to have charge-solvated structures. For all three dimers, the most stable conformer is found to be of type A. The spectrum of Met₂H⁺, however, cannot be explained without some abundance of type B charge-solvated conformers as well as salt-bridged structures.

Protonated glycine supramolecular systems: the need for quantum dynamics

Quantum mechanical simulations unequivocally explain experimental IR spectra of protonated supramolecular systems.

IR spectroscopy is one of the most commonly employed techniques to study molecular vibrations and interactions. However, characterization of experimental IR spectra is not always straightforward. This is the case of protonated glycine supramolecular systems like Gly₂H⁺ and (GlyH + nH₂), whose IR spectra raise questions which have still to find definitive answers even after theoretical spectroscopy investigations. Specifically, the assignment of the conformer responsible for the spectrum of the protonated glycine dimer (Gly₂H⁺) has led to much controversy and it is still debated, while structural hypotheses formulated to explain the main experimental spectral features of (GlyH + nH₂) systems have not been theoretically confirmed. We demonstrate that simulations must account for quantum dynamical effects in order to resolve these open issues. This is achieved by means of our divide-and-conquer semiclassical initial value representation technique, which approximates the quantum dynamics of high dimensional systems with remarkable accuracy and outperforms not only the commonly employed but unfit scaled-harmonic approaches, but also pure classical dynamics simulations. Besides the specific insights concerning the two particular cases presented here, the general conclusion is that, due to the widespread presence of protonated systems in chemistry, quantum dynamics may play a prominent role and should not be totally overlooked even when dealing with large systems including biological structures.

Anharmonic vibrational spectra and mode-mode couplings analysis of 2-aminopyridine

Vibrational spectra of 2-aminopyridine (2AP) have been analyzed using the vibrational self-consistence field theory (VSCF), correlated corrected vibrational self-consistence field theory (CC-VSCF) and vibrational perturbation theory (VPT2) at B3LYP/6-311G(d,p) framework. The mode-mode couplings affect the vibrational frequencies and intensities. The coupling integrals between pairs of normal modes have been obtained on the basis of quartic force field (2MR-QFF) approximation. The overtone and combination bands are also assigned in the FTIR spectrum with the help of anharmonic calculation at VPT2 method. A statistical analysis of deviations shows that estimated anharmonic frequencies are closer to the experiment over harmonic approximation. Furthermore, the anharmonic correction has also been carried out for the dimeric structure of 2AP. The fundamental vibration bands have been assigned on the basis of potential energy distribution (PED) and visual look over the animated modes. Other important molecular properties such as frontier molecular orbitals and molecular electrostatics potential mapping have also been analyzed.

A Decapeptide Hydrated by Two Waters: Conformers Determined by Theory and Validated by Cold Ion Spectroscopy

The intrinsic structures of biomolecules in the gas phase may not reflect their native solution geometries. Microsolvation of the molecules bridges the two environments, enabling a tracking of molecular structural changes upon hydration at the atomistic level. We employ density functional calculations to compute a large pool of structures and vibrational spectra for a gas-phase complex, in which a doubly protonated decapeptide, gramicidin S, is solvated by two water molecules. Though most vibrations of this large complex are treated in a harmonic approximation, the water molecules and the vibrations of the host ion coupled to them are locally described by a quantum mechanical vibrational self-consistent field theory with second-order perturbation correction (VSCF-PT2). Guided and validated by the available cold ion spectroscopy data, the computational analysis identifies structures of the three experimentally observed conformers of the complex. They, mainly, differ by the hydration sites, of which the one at the Orn side chain is the most important for reshaping the peptide toward its native structure. The study demonstrates the ability of a quantum chemistry approach that intelligently combines the semiempirical and ab initio computations to disentangle a complex interplay of intra- and intermolecular hydrogen bonds in large molecular systems.

Interpretation of semiclassical transition moments through wave function expansion of dipole moment functions with applications to the OH stretching spectra of simple acids and alcohols

Semiclassical description of molecular vibrations has provided us with various computational approximations and enhanced our conceptual understanding of quantum mechanics. In this study, the transition moments of the OH stretching fundamental and overtone intensities (Δv = 1-6) of some alcohols and acids are calculated by three kinds of semiclassical methods, correspondence-principle (CP) approximation, quasiclassical approximation, and uniform WKB approximation, and their respective transition moments are compared to those by the quantum theory. On the basis of the local mode picture, the one-dimensional potential energy curves and the dipole moment functions (DMFs) were obtained by density functional theory calculations and then fitted to Morse functions and sixth-order polynomials, respectively. It was shown that both the transition energies and the absorption intensities derived in the semiclassical methods reproduced their respective quantum values. In particular, the CP approximation reproduces the quantum transition moments if the formula given by Naccache is used for the action integral value. On the basis of these semiclassical results, we present a picture to understand the small variance in the overtone intensities of these acids and alcohols. Another important result is the ratios of semiclassical-to-quantum transition moment are almost independent of the applied molecules even with a great molecular variance of the DMFs, and they depend only on the nature of the semiclassical approximations and the quantum number. The difference between the semiclassical and quantum transition moments was analyzed in terms of a hitherto unrecognized concept that the Fourier expansion of the time dependent DMF in the CP treatment is a kind of the wave function expansion method using trigonometric functions as the quotient functions. For a Morse oscillator, we derive the analytic and approximate expressions of the quotient functions in terms of the bond displace coordinate in both the CP and the quantum mechanical frameworks and discuss the methodological dependence of the calculated transition moments. As a byproduct, we have found a simple derivation of the DMF expression first derived by Timm and Mecke long time ago.