Vibrational Fermi Resonance in Atomically Thin Black Phosphorus

Fermi resonance is a phenomenon involving the hybridization of two coincidentally quasi-degenerate states that is observed in the vibrational or electronic spectra of molecules. Despite numerous examples in molecular systems, vibrational Fermi resonances in dispersive semiconducting systems remain largely unexplored due to the rarity of occurrence. Here we report a vibrational Fermi resonance in atomically thin black phosphorus. The Fermi resonance arises via anharmonic mixing of a fundamental Raman mode and a Davydov component of an infrared mode, leading to a doublet with mixed character. The extent of Fermi coupling can be modulated by the application of external biaxial strain. The consequences of Fermi hybridization are revealed by electronic resonance effects in the thickness-dependent and excitation-wavelength-dependent Raman spectrum, which is predicted by ab initio hybrid functional simulations including excitonic interactions. This work reveals new insight into electron-phonon coupling in black phosphorus and demonstrates a novel method for modulating Fermi resonances in 2D semiconductors.

Vibrational Entanglement through the Lens of Quantum Information Measures

We introduce a quantum information analysis of vibrational wave functions to understand complex vibrational spectra of molecules with strong anharmonic couplings and vibrational resonances. For this purpose, we define one- and two-modal entropies to guide the identification of strongly coupled vibrational modes and to characterize correlations within modal basis sets. We evaluate these descriptors for multiconfigurational vibrational wave functions which we calculate with the n-mode vibrational density matrix renormalization group algorithm. Based on the quantum information measures, we present a vibrational entanglement analysis of the vibrational ground and excited states of CO2, which display strong anharmonic effects due to the symmetry-induced and accidental (near-) degeneracies. We investigate the entanglement signature of the Fermi resonance and discuss the maximally entangled state arising from the two degenerate bending modes.

Similarity scores of vibrational spectra reveal the atomistic structure of pentapeptides in multiple basins

Vibrational spectroscopy combined with theoretical calculations is a powerful tool for analyzing the interaction and conformation of peptides at the atomistic level. Nonetheless, identifying the structure becomes increasingly difficult as the peptide size grows large. One example is acetyl-SIVSF-N-methylamide, a capped pentapeptide, whose atomistic structure has remained unknown since its first observation [T. Sekiguchi, M. Tamura, H. Oba, P. Çarçarbal, R. R. Lozada-Garcia, A. Zehnacker-Rentien, G. Grégoire, S. Ishiuchi and M. Fujii, Angew. Chem., Int. Ed., 2018, 57, 5626-5629]. Here, we propose a novel conformational search method, which exploits the structure-spectrum correlation using a similarity score that measures the agreement of theoretical and experimental spectra. Surprisingly, the two conformers have distinctly different energy and geometry. The second conformer is 25 kJ mol-1 higher in energy than the other, lowest-energy conformer. The result implies that there are multiple pathways in the early stage of the folding process: one to the global minimum and the other to a different basin. Once such a structure is established, the second conformer is unlikely to overcome the barrier to produce the most stable structure due to a vastly different hydrogen bond network of the backbone. Our proposed method can characterize the lowest-energy conformer and kinetically trapped, high-energy conformers of complex biomolecules.

Vibrational heat-bath configuration interaction with semistochastic perturbation theory using harmonic oscillator or VSCF modals

Vibrational heat-bath configuration interaction (VHCI)-a selected configuration interaction technique for vibrational structure theory-has recently been developed in two independent works [J. H. Fetherolf and T. C. Berkelbach, J. Chem. Phys. 154, 074104 (2021); A. U. Bhatty and K. R. Brorsen, Mol. Phys. 119, e1936250 (2021)], where it was shown to provide accuracy on par with the most accurate vibrational structure methods with a low computational cost. Here, we eliminate the memory bottleneck of the second-order perturbation theory correction using the same (semi)stochastic approach developed previously for electronic structure theory. This allows us to treat, in an unbiased manner, much larger perturbative spaces, which are necessary for high accuracy in large systems. Stochastic errors are easily controlled to be less than 1 cm-1. We also report two other developments: (i) we propose a new heat-bath criterion and an associated exact implicit sorting algorithm for potential energy surfaces expressible as a sum of products of one-dimensional potentials; (ii) we formulate VHCI to use a vibrational self-consistent field (VSCF) reference, as opposed to the harmonic oscillator reference configuration used in previous reports. Our tests are done with quartic and sextic force fields, for which we find that with VSCF, the minor improvements to accuracy are outweighed by the higher computational cost associated the matrix element evaluations. We expect VSCF-based VHCI to be important for more general potential representations, for which the harmonic oscillator basis function integrals are no longer analytic.

An approximation to the vibrational coupled-cluster method for CH-stretching of large molecules: application to naphthalene and anthracene

We propose an approximation to the vibrational coupled-cluster method (VCCM) to describe the CH-stretching region of the vibrational spectrum of large molecules. The vibrational modes of a molecule are divided into two sets: the target set and the bath set. The target set includes the CH stretches and the modes that are strongly coupled with the CH stretches and/or involve strong Fermi resonances with a CH stretch fundamental. The rest of the modes are in the bath set. First, the effective harmonic oscillator (EHO) approximation is invoked for the whole system to obtain the zeroth-order frequencies and modified potentials. The effects of interaction between the bath set and the target sets are included in the modified potential from the EHO calculation. The VCCM equations are constructed with the modified potential from the EHO calculations and for the target set only. The transition energies and intensities are calculated using such a truncated VCCM approximation. The proposed method is applied to calculate the IR spectra of naphthalene and anthracene. The results with three different criteria for selecting the modes in the target set are compared with the experimental IR spectra.

Spectroscopic and Theoretical Studies on Conformational Stability of Benzyl Methyl Ether

Conformer-selected electronic and vibrational spectra of benzyl methyl ether and its terminal methyl group-substituted derivative in a supersonic jet have been measured using ultraviolet (UV)-UV hole burning and fluorescence-detected infrared spectroscopy to investigate the conformational stability of flexible molecules. Various quantum chemical calculations as well as experimental observations reveal the coexistence of three conformers with different CCOC dihedral angles and side-chain orientations relative to the benzene ring plane. Vibrational analysis in the excited state with time-dependent density functional theory and IR simulations containing anharmonic coupling sufficiently reproduce the experimental results, suggesting that these three conformers can be distinguished into one gauche-conformer and two trans-ones with respect to the CCOC dihedral angle. We also observe that the gauche conformer exhibits higher-frequency CH₂ modes. The natural bond orbital analysis indicates that this phenomenon is attributed to the electron delocalization from the non-bonding orbitals and the C-O orbitals associated with the neighboring oxygen atom, which leads to a conformer dependence of the methylene C-H bond strength.

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).

Retinal Vibrations in Bacteriorhodopsin are Mechanically Harmonic but Electrically Anharmonic: Evidence From Overtone and Combination Bands

Fundamental vibrations of the chromophore in the membrane protein bacteriorhodopsin (BR), a protonated Schiff base retinal, have been studied for decades, both by resonance Raman and by infrared (IR) difference spectroscopy. Such studies started comparing vibrational changes between the initial BR state (all-trans retinal) and the K intermediate (13-cis retinal), being later extended to the rest of intermediates. They contributed to our understanding of the proton-pumping mechanism of BR by exploiting the sensitivity of fundamental vibrational transitions of the retinal to its conformation. Here, we report on new bands in the 2,500 to 1,800 cm⁻¹ region of the K-BR difference FT-IR spectrum. We show that the bands between 2,500 and 2,300 cm⁻¹ originate from overtone and combination transitions from C-C stretches of the retinal. We assigned bands below 2,300 cm⁻¹ to the combination of retinal C-C stretches with methyl rocks and with hydrogen-out-of-plane vibrations. Remarkably, experimental C-C overtone bands appeared at roughly twice the wavenumber of their fundamentals, with anharmonic mechanical constants ≤3.5 cm⁻¹, and in some cases of ∼1 cm⁻¹. Comparison of combination and fundamental bands indicates that most of the mechanical coupling constants are also very small. Despite the mechanical quasi-harmonicity of the C-C stretches, the area of their overtone bands was only ∼50 to ∼100 times smaller than of their fundamental bands. We concluded that electrical anharmonicity, the second mechanism giving intensity to overtone bands, must be particularly high for the retinal C-C stretches. We corroborated the assignments of negative bands in the K-BR difference FT-IR spectrum by ab initio anharmonic vibrational calculations of all-trans retinal in BR using a quantum-mechanics/molecular mechanics approach, reproducing reasonably well the small experimental anharmonic and coupling mechanical constants. Yet, and in spite accounting for both mechanical and electrical anharmonicities, the intensity of overtone C-C transitions was underestimated by a factor of 4–20, indicating room for improvement in state-of-the-art anharmonic vibrational calculations. The relatively intense overtone and combination bands of the retinal might open the possibility to detect retinal conformational changes too subtle to significantly affect fundamental transitions but leaving a footprint in overtone and combination transitions.

Sum-frequency vibrational spectroscopy of methanol at interfaces due to Fermi resonance

We present a theoretical method of studying sum-frequency vibrational spectroscopy for the CH3 group of methanol at interfaces due to Fermi resonance, which provides a novel and untraditional point of view with respect to traditional approaches.

Towards complete assignment of the infrared spectrum of the protonated water cluster H+(H2O)21

The spectroscopic features of protonated water species in dilute acid solutions have been long sought after for understanding the microscopic behavior of the proton in water with gas-phase water clusters H⁺(H₂O)_n extensively studied as bottom-up model systems. We present a new protocol for the calculation of the infrared (IR) spectra of complex systems, which combines the fragment-based Coupled Cluster method and anharmonic vibrational quasi-degenerate perturbation theory, and demonstrate its accuracy towards the complete and accurate assignment of the IR spectrum of the H⁺(H₂O)₂₁ cluster. The site-specific IR spectral signatures reveal two distinct structures for the internal and surface four-coordinated water molecules, which are ice-like and liquid-like, respectively. The effect of inter-molecular interaction between water molecules is addressed, and the vibrational resonance is found between the O-H stretching fundamental and the bending overtone of the nearest neighboring water molecule. The revelation of the spectral signature of the excess proton offers deeper insight into the nature of charge accommodation in the extended hydrogen-bonding network underpinning this aqueous cluster.

Protonated water species have been the subject of numerous experimental and computational studies. Here the authors provide a nearly complete assignment of the experimental IR spectrum of the H⁺(H₂O)₂₁ water cluster based on high-level wavefunction theory and anharmonic vibrational quasi-degenerate perturbation theory.

Advances, challenges and perspectives of quantum chemical approaches in molecular spectroscopy of the condensed phase

The purpose of this review is to demonstrate advances, challenges and perspectives of quantum chemical approaches in molecular spectroscopy of the condensed phase. Molecular spectroscopy, particularly vibrational spectroscopy and electronic spectroscopy, has been used extensively for a wide range of areas of chemical sciences and materials science as well as nano- and biosciences because it provides valuable information about structure, functions, and reactions of molecules. In the meantime, quantum chemical approaches play crucial roles in the spectral analysis. They also yield important knowledge about molecular and electronic structures as well as electronic transitions. The combination of spectroscopic approaches and quantum chemical calculations is a powerful tool for science, in general. Thus, our article, which treats various spectroscopy and quantum chemical approaches, should have strong implications in the wider scientific community. This review covers a wide area of molecular spectroscopy from far-ultraviolet (FUV, 120-200 nm) to far-infrared (FIR, 400-10 cm-1)/terahertz and Raman spectroscopy. As quantum chemical approaches, we introduce several anharmonic approaches such as vibrational self-consistent field (VSCF) and the combination of periodic harmonic calculations with anharmonic corrections based on finite models, grid-based techniques like the Numerov approach, the Cartesian coordinate tensor transfer (CCT) method, Symmetry-Adapted Cluster Configuration-Interaction (SAC-CI), and the ZINDO (Semi-empirical calculations at Zerner's Intermediate Neglect of Differential Overlap). One can use anharmonic approaches and grid-based approaches for both infrared (IR) and near-infrared (NIR) spectroscopy, while CCT methods are employed for Raman, Raman optical activity (ROA), FIR/terahertz and low-frequency Raman spectroscopy. Therefore, this review overviews cross relations between molecular spectroscopy and quantum chemical approaches, and provides various kinds of close-reality advanced spectral simulation for condensed phases.

Anharmonic Vibrational Calculations Based on Group-Localized Coordinates: Applications to Internal Water Molecules in Bacteriorhodopsin

An efficient anharmonic vibrational method is developed exploiting the locality of molecular vibration. Vibrational coordinates localized to a group of atoms are employed to divide the potential energy surface (PES) of a system into intra- and inter-group contributions. Then, the vibrational Schrödinger equation is solved based on a PES, in which the inter-group coupling is truncated at the harmonic level while accounting for the intra-group anharmonicity. The method is applied to a pentagonal hydrogen bond network (HBN) composed of internal water molecules and charged residues in a membrane protein, bacteriorhodopsin. The PES is calculated by the quantum mechanics/molecular mechanics (QM/MM) calculation at the level of B3LYP-D3/aug-cc-pVDZ. The infrared (IR) spectrum is computed using a set of coordinates localized to each water molecule and amino acid residue by second-order vibrational quasi-degenerate perturbation theory (VQDPT2). Benchmark calculations show that the proposed method yields the N-D/O-D stretching frequencies with an error of 7 cm-1 at the cost reduced by more than five times. In contrast, the harmonic approximation results in a severe error of 150 cm-1. Furthermore, the size of QM regions is carefully assessed to find that the QM regions should include not only the pentagonal HBN itself but also its HB partners. VQDPT2 calculations starting from transient structures obtained by molecular dynamics simulations have shown that the structural sampling has a significant impact on the calculated IR spectrum. The incorporation of anharmonicity, sufficiently large QM regions, and structural samplings are of essential importance to reproduce the experimental IR spectrum. The computational spectrum paves the way for decoding the IR signal of strong HBNs and helps elucidate their functional roles in biomolecules.

Reduced Two-Electron Interactions in Anharmonic Molecular Vibrational Calculations Involving Localized Normal Coordinates

Spatially localized vibrational normal mode coordinates are shown to reduce the importance of calculating the full set of two-electron terms in the molecular electronic Schrödinger equation. Electron correlation and dispersion interactions become less significant in (E,E)-1,3,5,7-octatetraene vibrational self-consistent field calculations when displacing remote atoms along multiple coordinates. Electron correlation interactions between spatially remote modes are also found to be less important compared to their corresponding uncorrelated interaction terms. Attenuation of the Coulomb operator indicates that the two-electron terms between remote electrons become less important for accurately describing the strongly contributing mode-coupling terms between sets of localized vibrational modes.

Exploring the anharmonic vibrational structure of carbon dioxide trimers

Our previously developed mbCO2 potential [O. Sode and J. N. Cherry, J. Comput. Chem. 38, 2763 (2017)] is used to describe the vibrational structure of the intermolecular motions of the CO₂ trimers: barrel-shaped and cyclic trimers. Anharmonic corrections are accounted for using the vibrational self-consistent field theory, vibrational second-order Møller-Plesset perturbation (VMP2) theory, and vibrational configuration interaction (VCI) methods and compared with experimental observations. For the cyclic structure, we revise the assignments of two previously observed experimental peaks based on our VCI and VMP2 results. We note that the experimental band observed near 13 cm^-1 is the out-of-phase out-of-plane degenerate motion with E″ symmetry, while the peak observed at 18 cm^-1 likely corresponds to the symmetric out-of-plane torsion A″ vibration. Since the VCI treatment of the vibrational motions accounts for vibrational mixing and delocalization, overtones and combination bands were also observed and quantified in the intermolecular regions of the two trimer isomers.

Anharmonic Phonon Dispersion in Polyethylene

The second-order Green's function method for anharmonic crystals has been applied to an infinite, periodic chain of polyethylene taking into account up to quartic force constants. The frequency-independent approximation to the Dyson self-energy gives rise to numerous divergent resonances, which are fortuitous. Instead, solving the Dyson equation self-consistently with a frequency-dependent self-energy resists divergences from resonances or zero-frequency acoustic vibrations. The calculated anharmonic phonon dispersion, which nonetheless displays many true resonances, and anharmonic phonon density of states furnish hitherto unknown details that explain smaller features of observed vibrational spectra.

Chiral Inversion of Amino Acids in Antiparallel β-Sheets at Interfaces Probed by Vibrational Sum Frequency Generation Spectroscopy

A parallel study of protein variants with all (l-), all (d-), or mixed (l-)/(d-) amino acids can be used to assess how backbone architecture versus side chain identity determines protein structure. Here, we investigate the secondary structure and side chain orientation dynamics of the antiparallel β-sheet peptide LK7β (Ac-Leu-Lys-Leu-Lys-Leu-Lys-Leu-NH2) composed of all (l-), all (d-), or alternating (l-Leu)/(d-Lys) amino acids. Using interface-selective vibrational sum frequency generation spectroscopy (VSFG), we observe that the alternating (l-)/(d-) peptide lacks a resonant C-H stretching mode compared to the (l-) and (d-) variants and does not form antiparallel β-sheets. We rationalize our observations on the basis of density functional theory calculations and molecular dynamics (MD) simulations of LK7β at the air-water interface. Irrespective of the handedness of the amino acids, leucine side chains prefer to orient toward the hydrophobic air phase while lysine side chains prefer the hydrophilic water phase. These preferences dictate the backbone configuration of LK7β and thereby the folding of the peptide. Our MD simulations show that the preferred side chain orientations can force the backbone of a single strand of (l-) LK7β at the air-water interface to adopt β-sheet Ramachandran angles. However, denaturation of the β-sheets at pH = 2 results in a negligible chiral VSFG amide I response. The combined computational and experimental results lend critical support to the theory that a chiral VSFG response requires macroscopic chirality, such as in β-sheets. Our results can guide expectations about the VSFG optical responses of proteins and should improve understanding of how amino acid chirality modulates the structure and function of natural and de novo proteins at biological interfaces.

Efficiently Calculating Anharmonic Frequencies of Molecular Vibration by Molecular Dynamics Trajectory Analysis

Two efficient methods, the Eckart frame algorithm and the multiorder derivative algorithm, for vibrational frequency calculation directly based on the raw data of atomic trajectory from the state-of-the-art first-principles molecular dynamics simulation are presented. The Eckart frame approach is robust to retrieve the full set of anharmonic fundamental frequencies of any molecule from the atomic trajectory for a sufficiently long molecular dynamics simulation at a temperature close to 0 K. In addition to the fundamental vibrational frequencies, the multiorder derivative approach is universal for the calculations of vibrational frequencies based on the molecular dynamics result in a wide range of temperatures. The accuracy, efficiency, and applicability of these two methods are demonstrated through several successful examples in calculating the anharmonic fundamental vibrational frequencies of methane, ethylene, water, and cyclobutadiene.

Understanding the anharmonic vibrational structure of the carbon dioxide dimer

Understanding the vibrational structure of the CO₂ system is important to confirm the potential energy surface and interactions in such van der Waals complexes. In this work, we use our previously developed mbCO2 potential function to explore the vibrational structure of the CO₂ monomer and dimer. The potential function has been trained to reproduce the potential energies at the CCSD(T)-F12b/aug-cc-pVTZ level of electronic structure theory. The harmonic approximation, as well as anharmonic corrections using vibrational structure theories such as vibrational self-consistent field, vibrational second-order Møller-Plesset perturbation, and vibrational configuration interaction (VCI), is applied to address the vibrational motions. We compare the vibrational results using the mbCO2 potential function with traditional electronic structure theory results and to experimental frequencies. The anharmonic results for the monomer most closely match the experimental data to within 3 cm^-1, including the Fermi dyad frequencies. The intermolecular and intramolecular dimer frequencies were treated separately and show good agreement with the most recent theoretical and experimental results from the literature. The VCI treatment of the dimer vibrational motions accounts for vibrational mixing and delocalization, such that we observe the dimer Fermi resonance phenomena, both in the intramolecular and intermolecular regions.

Anharmonic Vibrational Analysis of Biomolecules and Solvated Molecules Using Hybrid QM/MM Computations

Quantum mechanics/molecular mechanics (QM/MM) calculations are applied for anharmonic vibrational analyses of biomolecules and solvated molecules. The QM/MM method is implemented into a molecular dynamics (MD) program, GENESIS, by interfacing with external electronic structure programs. Following the geometry optimization and the harmonic normal-mode analysis based on a partial Hessian, the anharmonic potential energy surface (PES) is generated from QM/MM energies and gradients calculated at grid points. The PES is used for vibrational self-consistent field (VSCF) and post-VSCF calculations to compute the vibrational spectrum. The method is first applied to a phosphate ion in solution. With both the ion and neighboring water molecules taken as a QM region, IR spectra of representative hydration structures are calculated by the second-order vibrational quasi-degenerate perturbation theory (VQDPT2) at the level of B3LYP/cc-pVTZ and TIP3P force field. A weight-average of IR spectra over the structures reproduces the experimental spectrum with a mean absolute deviation of 16 cm^–1. Then, the method is applied to an enzyme, P450 nitric oxide reductase (P450nor), with the NO molecule bound to a ferric (Fe^III) heme. Starting from snapshot structures obtained from MD simulations of P450nor in solution, QM/MM calculations have been carried out at the level of B3LYP-D3/def2-SVP(D). The spin state of Fe^III(NO) is likely a closed-shell singlet state based on a ratio of N–O and Fe–NO stretching frequencies (ν_N–O and ν_Fe–NO) calculated for closed- and open-shell singlet states. The calculated ν_N–O and ν_Fe–NO overestimate the experimental ones by 120 and 75 cm^–1, respectively. The electronic structure and solvation of Fe^III(NO) affect the structure around the heme of P450nor leading to an increase in ν_N–O and ν_Fe–NO.

Reduced Basis Set Dependence in Anharmonic Frequency Calculations Involving Localized Coordinates

Localized normal coordinates are known to be effective in speeding up anharmonic frequency calculations by reducing the complexity of the nuclear Hamiltonian and wave function. Displacing atoms in localized coordinates can also cause relatively small changes in the electronic structure, which can be exploited for further computational efficiency improvements during ab initio electronic structure calculations of the potential energy surface by reducing the electronic basis set dependence. Three different schemes for reducing the basis set dependence have been investigated in this work. These include combining localized coordinate schemes with general mixed basis sets, distance based force-field reductions, and using coordinate specific basis sets. The importance of accurately describing electronic interactions is found to diminish both for multicoordinate terms involving the displacement of remote atoms and when describing the interactions between more remote atoms within specific coordinates.

Hydrogen bond induced enhancement of Fermi resonances in N–H⋯N hydrogen bonded complexes of anilines

Enhancement of Fermi resonance intensities due to the formation of N–H⋯N hydrogen bonding of anilines with alkyl amines is analyzed using a two-state deperturbation model.

Intermediate vibrational coordinate localization with harmonic coupling constraints

Optimized normal coordinates can significantly improve the speed and accuracy of vibrational frequency calculations. However, over-localization can occur when using unconstrained spatial localization techniques. The unintuitive mixtures of stretching and bending coordinates that result can make interpreting spectra more difficult and also cause artificial increases in mode-coupling during anharmonic calculations. Combining spatial localization with a constraint on the coupling between modes can be used to generate coordinates with properties in-between the normal and fully localized schemes. These modes preserve the diagonal nature of the mass-weighted Hessian matrix to within a specified tolerance and are found to prevent contamination between the stretching and bending vibrations of the molecules studied without a priori classification of the different types of vibration present. Relaxing the constraint can also be used to identify which normal modes form specific groups of localized modes. The new coordinates are found to center on more spatially delocalized functional groups than their fully localized counterparts and can be used to tune the degree of vibrational correlation energy during anharmonic calculations.

Trimethylamine-N-oxide: its hydration structure, surface activity, and biological function, viewed by vibrational spectroscopy and molecular dynamics simulations

Vibrational spectroscopy and molecular simulations revealed the hydrophilicity and hydrophobicity of TMAO in aqueous solution.

Examining the impact of harmonic correlation on vibrational frequencies calculated in localized coordinates

Carefully choosing a set of optimized coordinates for performing vibrational frequency calculations can significantly reduce the anharmonic correlation energy from the self-consistent field treatment of molecular vibrations. However, moving away from normal coordinates also introduces an additional source of correlation energy arising from mode-coupling at the harmonic level. The impact of this new component of the vibrational energy is examined for a range of molecules, and a method is proposed for correcting the resulting self-consistent field frequencies by adding the full coupling energy from connected pairs of harmonic and pseudoharmonic modes, termed vibrational self-consistent field (harmonic correlation). This approach is found to lift the vibrational degeneracies arising from coordinate optimization and provides better agreement with experimental and benchmark frequencies than uncorrected vibrational self-consistent field theory without relying on traditional correlated methods.

Efficient anharmonic vibrational spectroscopy for large molecules using local-mode coordinates

This article presents a general computational approach for efficient simulations of anharmonic vibrational spectra in chemical systems. An automated local-mode vibrational approach is presented, which borrows techniques from localized molecular orbitals in electronic structure theory. This approach generates spatially localized vibrational modes, in contrast to the delocalization exhibited by canonical normal modes. The method is rigorously tested across a series of chemical systems, ranging from small molecules to large water clusters and a protonated dipeptide. It is interfaced with exact, grid-based approaches, as well as vibrational self-consistent field methods. Most significantly, this new set of reference coordinates exhibits a well-behaved spatial decay of mode couplings, which allows for a systematic, a priori truncation of mode couplings and increased computational efficiency. Convergence can typically be reached by including modes within only about 4 Å. The local nature of this truncation suggests particular promise for the ab initio simulation of anharmonic vibrational motion in large systems, where connection to experimental spectra is currently most challenging.

Approximate first-principles anharmonic calculations of polyatomic spectra using MP2 and B3LYP potentials: comparisons with experiment

Anharmonic vibrational spectroscopy calculations using MP2 and B3LYP computed potential surfaces are carried out for a series of molecules, and frequencies and intensities are compared with those from experiment. The vibrational self-consistent field with second-order perturbation correction (VSCF-PT2) is used in computing the spectra. The test calculations have been performed for the molecules HNO3, C2H4, C2H4O, H2SO4, CH3COOH, glycine, and alanine. Both MP2 and B3LYP give results in good accord with experimental frequencies, though, on the whole, MP2 gives very slightly better agreement. A statistical analysis of deviations in frequencies from experiment is carried out that gives interesting insights. The most probable percentage deviation from experimental frequencies is about -2% (to the red of the experiment) for B3LYP and +2% (to the blue of the experiment) for MP2. There is a higher probability for relatively large percentage deviations when B3LYP is used. The calculated intensities are also found to be in good accord with experiment, but the percentage deviations are much larger than those for frequencies. The results show that both MP2 and B3LYP potentials, used in VSCF-PT2 calculations, account well for anharmonic effects in the spectroscopy of molecules of the types considered.