Abstract: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.The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jctc.1c00060.Computational details including all-atom model of BR, settings of MD simulations, QM/MM calculations, geometry optimization, and vibrational calculations, protocol of equilibration, modified CHARMM force field parameters for retinal, harmonic frequencies obtained from partial Hessian matrices, effect of the size of a group, VQDPT2 frequencies obtained by full and reduced dimensional calculations for each group, harmonic frequencies and IR intensities in terms of normal and local coordinates, vibrational frequencies and IR intensities, harmonic frequencies obtained by QM/MM calculations using QM1–QM4, IR spectra obtained by multiresolution PES and 2MR-GridPES and 3MR-QFF, IR spectra obtained by the harmonic approximation at the level of RI-MP2 and B3LYP-D3, IR difference spectrum (K-BR) calculated by VQDPT2 at the level of B3LYP-D3/aug-cc-pVDZ using an isolated system, snapshot structure of a membrane system, plots of rmsd, P–P distance as a function of simulation time, VQDPT2 spectra obtained from snapshot structures taken every 2 ns between 6 and 20 ns, VQDPT2 spectra obtained from X = X(t) and X = Xeq, and VQDPT2 spectra obtained from eight snapshot structures that satisfy r1 > 5.1 Å and r2 > 1.2 Å (PDF)This article has not yet been cited by other publications.

TD-DFT analysis of excitation of the closed-form spironaphthopyran in methanol solution: The contribution of vibronic transitions and intermolecular hydrogen bonds

Exploring non-covalent interactions in excited states: beyond aromatic excimer models

TD‐DFT analysis of the excitation of H‐dimers of cationic dyes in an aqueous solution using functionals without additional dispersion correction

Optical Excitations of Chlorophyll $a$ and $b$ Monomers and Dimers

Intramolecular excimers of open forms of 2H-benzopyran, 2H- and 3H-naphthopyrans in solution: TD-DFT/DFT analysis

Experimental and Theoretical Electron Momentum Spectroscopic Study of the Valence Electronic Structure of Tetrahydrofuran under Pseudorotation

Theoretical investigations on the excited-state intramolecular proton transfer in the solvated 2-hydroxy-1-naphthaldehyde carbohydrazone

Photoisomerization of Spiropyran in Solution: A Surface Hopping Investigation

Solvent conditions effect on the excited state intramolecular proton transfer mechanism and photophysical property of 1′-hydroxy-2′-acetonaphthone: A DFT/TD-DFT analysis

A novel non-fluorescent excited state intramolecular proton transfer phenomenon induced by intramolecular hydrogen bonds: an experimental and theoretical investigation

Excited State Electronic Structure of Dimethyl Disulfide Involved in Photodissociation at ∼200 nm

Excited-state Proton Transfer of 4-Hydroxyl-1, 8-Naphthalimide Derivatives: A Combined Experimental and Theoretical Investigation

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

Intramolecular proton or hydrogen-atom transfer in the ground- and excited-states of 2-hydroxybenzophenone: A theoretical study

Intramolecular O-H⋯S hydrogen bonding in threefold symmetry: Line broadening dynamics from ultrafast 2DIR-spectroscopy and ab initio calculations

How do spin‐scaled double hybrids designed for excitation energies perform for noncovalent excited‐state interactions? An investigation on aromatic excimer models

The interplay of excitonic delocalization and vibrational localization in optical lineshapes: A variational polaron approach

Photoexcitation Dynamics of 4-Aminopthalimide in Solution Investigated Using Femtosecond Time-Resolved Infrared Spectroscopy

How do spin-scaled double hybrids designed for excitation energies perform for noncovalent excited-state interactions? An investigation on aromatic excimer models

DFT/TDDFT based study to decipher the excited state intramolecular proton transfer process for 3-benzothiazol-2-yl-2-hydroxy-5-methyl-benzaldehyde-oxime