A Direct, Quantitative Connection between Molecular Dynamics Simulations and Vibrational Probe Line Shapes

A quantitative connection between molecular dynamics simulations and vibrational spectroscopy of probe-labeled systems would enable direct translation of experimental data into structural and dynamical information. To constitute this connection, all-atom molecular dynamics (MD) simulations were performed for two SCN probe sites (solvent-exposed and buried) in a calmodulin-target peptide complex. Two frequency calculation approaches with substantial nonelectrostatic components, a quantum mechanics/molecular mechanics (QM/MM)-based technique and a solvatochromic fragment potential (SolEFP) approach, were used to simulate the infrared probe line shapes. While QM/MM results disagreed with experiment, SolEFP results matched experimental frequencies and line shapes and revealed the physical and dynamic bases for the observed spectroscopic behavior. The main determinant of the CN probe frequency is the exchange repulsion between the probe and its local structural neighbors, and there is a clear dynamic explanation for the relatively broad probe line shape observed at the “buried” probe site. This methodology should be widely applicable to vibrational probes in many environments

Overcoming the Failure of Correlation for Out-of-Plane Motions in a Simple Aromatic: Rovibrational Quantum Chemical Analysis of <i>c</i>‑C₃H₂

Truncated, correlated, wave function methods either produce imaginary frequencies (in the extreme case) or nonphysically low frequencies in out-of-plane motions for carbon and adjacent atoms when the carbon atoms engage in π bonding. Cyclopropenylidene is viewed as the simplest aromatic hydrocarbon, and the present as well as previous theoretical studies have shown that this simple molecule exhibits this behavior in the two out-of-plane bends (OPBs). This nonphysical behavior has been treated by removing nearly linear dependent basis functions according to eigenvalues of the overlap matrix, by employing basis sets where the <i>spd</i> space saturatation is balanced with higher angular momentum functions, by including basis set superposition/incompleteness error (BSSE/BSIE) corrections, or by combining standard correlation methods with explicitly correlated methods to produce hybrid potential surfaces. However, this work supports the recently described hypothesis that the OPB problem is both a method and a basis set effect. The correlated wave function’s largest higher-order substitution term comes from a π → π* excitation where constructive interference of both orbitals artificially stabilizes the OPB. By employing schema to overcome this issue, the symmetric OPB ν₉ is the predicted to be the second-brightest transition, and it will be observed very close to 775 cm^–1. However, more work from the community is required to formulate better how carbon atoms interact with their adjacent atoms in π-bonded systems. Such bonds are ubiquitous in all of chemistry and beyond