Rapid anharmonic vibrational corrections derived from partial Hessian analysis

Vibrational analysis within a partial Hessian framework can successfully describe the vibrational properties of a variety of systems where the vibrational modes of interest are localized within a specific region of the system. We have developed a new approach to calculating anharmonic frequencies based on vibrational frequencies and normal modes obtained from a partial Hessian analysis using second-order vibrational perturbation theory and the transition optimized shifted Hermite method. This allows anharmonic frequencies for vibrational modes that are spatially localized to be determined at a significantly reduced computational cost. Several molecular systems are examined in order to demonstrate the effectiveness of this method including organic molecules adsorbed on the Si(100)-2×1 surface, model peptides in solution, and the C-H stretching region of polycyclic aromatic hydrocarbons. Overall, for a range of systems, anharmonic frequencies calculated using the partial Hessian approach are found to be in close agreement with the results obtained using full anharmonic calculations while providing a significant reduction in computational cost.     

Vibrational modes in partially optimized molecular systems

In this paper the authors develop a method to accurately calculate localized vibrational modes for partially optimized molecular structures or for structures containing link atoms. The method avoids artificially introduced imaginary frequencies and keeps track of the invariance under global translations and rotations. Only a subblock of the Hessian matrix has to be constructed and diagonalized, leading to a serious reduction of the computational time for the frequency analysis. The mobile block Hessian approach (MBH) proposed in this work can be regarded as an extension of the partial Hessian vibrational analysis approach proposed by Head [Int. J. Quantum Chem. 65, 827 (1997)]. Instead of giving the nonoptimized region of the system an infinite mass, it is allowed to move as a rigid body with respect to the optimized region of the system. The MBH approach is then extended to the case where several parts of the molecule can move as independent multiple rigid blocks in combination with single atoms. The merits of both models are extensively tested on ethanol and di-n-octyl-ether.     

Analytical hessian fitting schemes for efficient determination of force‐constant parameters in molecular mechanics

Building upon our recently developed partial Hessian fitting (PHF) method (Wang et al., J. Comput. Chem. 2016 , 37, 2349), we formulated and implemented two other rapid force‐field parameterization schemes called full Hessian fitting (FHF) and internal Hessian fitting (IHF), and comparisons were made among these three parameterization schemes to assess their performance. FHF minimizes deviation between the Hessian matrices in Cartesian coordinates computed by quantum mechanics (QM) and molecular mechanics (MM), to determine the best possible MM force‐constant parameters. While PHF requires step‐by‐step fittings of 3 × 3 partial Hessian matrices, FHF compares the lower triangular part of the QM and MM Hessian matrices, which allows simultaneous determination of all force‐constant parameters. In addition to this simple FHF scheme, IHF was developed such that it considers the Hessian matrices in redundant internal coordinates, where all possible internal coordinates that arise from the user‐defined interatomic connectivity are utilized. The results show that IHF performs best overall, followed by PHF and then FHF. Python‐based programing codes were developed to automate various tedious steps involved in the parameterization processes. © 2017 Wiley Periodicals, Inc.     

Partial Hessian vibrational analysis: the localization of the molecular vibrational energy and entropy

The partial Hessian vibrational analysis (PHVA), in which only a subblock of the Hesssian matrix is diagonalized to yield vibrational frequencies for partially optimized systems, is extended to the calculation of vibrational enthalpy and entropy changes for chemical reactions. The utility of this method is demonstrated for various deprotonation reactions by reproducing full HVA values to within 0.1–0.4 kcal/mol, depending on the number atoms included in the PHVA. When combined with the hybrid effective fragment potential method [Gordon MS, et al. (2001) J Phys Chem A 105:293–307], the PHVA method can provide (harmonic) free-energy changes for localized chemical reactions in very large systems. Received: 21 September 2001 / Accepted: 30 October 2001 / Published online: 22 March 2002     

Cartesian formulation of the mobile block Hessian approach to vibrational analysis in partially optimized systems

Partial optimization is a useful technique to reduce the computational load in simulations of extended systems. In such nonequilibrium structures, the accurate calculation of localized vibrational modes can be troublesome, since the standard normal mode analysis becomes inappropriate. In a previous paper [A. Ghysels et al., J. Chem. Phys. 126, 224102 (2007)], the mobile block Hessian (MBH) approach was presented to deal with the vibrational analysis in partially optimized systems. In the MBH model, the nonoptimized regions of the system are represented by one or several blocks, which can move as rigid bodies with respect to the atoms of the optimized region. In this way unphysical imaginary frequencies are avoided and the translational/rotational invariance of the potential energy surface is fully respected. In this paper we focus on issues concerning the practical numerical implementation of the MBH model. The MBH normal mode equations are worked out for several coordinate choices. The introduction of a consistent group-theoretical notation facilitates the treatment of both the case of a single block and the case of multiple blocks. Special attention is paid to the formulation in terms of Cartesian variables, in order to provide a link with the standard output of common molecular modeling programs.     

QM/MM vibrational mode tracking

Vibrational spectroscopy is a powerful tool to investigate the structure and dynamics of biomolecules. When small subsystems of large molecules such as active centers of enzymes are studied, quantum chemical calculations based on quantum mechanics/molecular mechanics (QM/MM) coupling schemes are a valuable means to interpret the spectra. The goal of this work is a methodological pilot study on how to selectively and thus efficiently extract certain vibrational information for extended molecular systems described by QM/MM methods. This is achieved by an extension of the mode tracking algorithm and a comparison with the partial Hessian diagonalization approach. After validating the methodology for the CO stretching vibration of 2-butanone and a delocalized CO stretch in acetylacetone, the stretching and bending modes of the CO ligand in CO myoglobin are tracked. Such systems represent an ideal application for mode tracking, because only a few strongly localized vibrations are sought for, while the large remainder of the molecule is of interest only as far as it affects these local vibrations. This influence is treated exactly by mode tracking.     

Acceleration of Correlation-corrected Vibrational Self-consistent Field Calculation Times for Large Polyatomic Molecules

Acceleration of the correlation-corrected Vibrational self-consistent field (CC-VSCF) method for anharmonic calculations of vibrational states of polyatomic molecules is described. The acceleration assumes pairwise additive interactions between different normal modes, and employs orthogonality of the single-mode vibrational wavefunctions. This greatly reduces the effort in computing correlation effects between different vibrational modes, which is treated by second order perturbation theory in CC-VSCF. The acceleration can improve the scaling of the overall computational effort from N ⁶ to N ⁴, where N is the number of vibrational modes. Sample calculation times, using semi-empirical potential surfaces (PM3), are given for a series of glycine peptides. Large computational acceleration, and significant reduction of the scaling of the effort with system size, is found and discussed.     

Vibrational Coupling between Organic and Inorganic Sublattices of Hybrid Perovskites

The interplay of electronic and nuclear degrees of freedom in semiconductor hybrid organic–inorganic perovskites determines many of their fundamental photophysical properties. For instance, charge carriers are dressed with phonons, that is, form polarons, and combination modes composed of strongly mixed localized vibrations and delocalized phonons can provide pathways for electronic energy relaxation and dissipation. Mixing of the different types of nuclear motion in vibrational combination modes requires their strong coupling. The direct measurement of coupling between the high‐frequency N?H stretch modes of the organic methylammonium and formamidinium ions and low‐frequency Pb?I phonon modes of the inorganic sub‐lattice in hybrid organic–inorganic perovskites is presented. The results reveal direct and substantial coupling between the non‐covalently interacting organic and inorganic sub‐lattices.     

Intensity‐Carrying Modes in Raman and Raman Optical Activity Spectroscopy

We describe a quantum‐chemical approach for the determination of modes with maximum Raman and Raman optical activity (ROA) intensity by maximizing the intensities with respect to the Raman and Raman optical activity intensity, respectively, which is shown to lead to eigenvalue equations. The intensity‐carrying modes are in general hypothetical modes and do not directly correspond to a certain normal mode in the spectrum. However, they provide information about those molecular distortions leading to intense bands in the spectrum. Modes with maximum Raman intensity are presented for propane‐1,3‐dione, propane‐1,3‐dionate, and Λ‐tris(propane‐1,3‐dionato)cobalt(III). Moreover, the mode with highest ROA intensity is examined for this chiral cobalt complex and also for the (chiral) amino acid L ‐tryptophan. The Raman and ROA high‐intensity modes are an optimal starting guess for intensity‐tracking calculations, in which selectively normal modes with high Raman or ROA intensity are converged. We present the first Raman and ROA intensity‐tracking calculations. These reveal a high potential for large molecules, for which the selective calculation of normal modes with high intensity is desirable in view of the large computational effort required for the calculation of Raman and ROA polarizability property tensors.     

Comparing normal modes across different models and scales: Hessian reduction versus coarse‐graining

Dimension reduction is often necessary when attempting to reach longer length and time scales in molecular simulations. It is realized by constraining degrees of freedom or by coarse‐graining the system. When evaluating the accuracy of a dimensional reduction, there is a practical challenge: the models yield vectors with different lengths, making a comparison by calculating their dot product impossible. This article investigates mapping procedures for normal mode analysis. We first review a horizontal mapping procedure for the reduced Hessian techniques, which projects out degrees of freedom. We then design a vertical mapping procedure for the “implosion” of the all‐atom (AA) Hessian to a coarse‐grained scale that is based upon vibrational subsystem analysis. This latter method derives both effective force constants and an effective kinetic tensor. Next, a series of metrics is presented for comparison across different scales, where special attention is given to proper mass‐weighting. The dimension‐dependent metrics, which require prior mapping for proper evaluation, are frequencies, overlap of normal mode vectors, probability similarity, Hessian similarity, collectivity of modes, and thermal fluctuations. The dimension‐independent metrics are shape derivatives, elastic modulus, vibrational free energy differences, heat capacity, and projection on a predefined basis set. The power of these metrics to distinguish between reasonable and unreasonable models is tested on a toy alpha helix system and a globular protein; both are represented at several scales: the AA scale, a Gō‐like model, a canonical elastic network model, and a network model with intentionally unphysical force constants. Published 2012 Wiley Periodicals, Inc.     

Analytical derivative techniques for TDDFT excited-state properties: Theory and application

We review our recent work on the methodology development of the excited-state properties for the molecules in vacuum and liquid solution.The general algorithms of analytical energy derivatives for the specific properties such as the first and second geometrical derivatives and IR/Raman intensities are demonstrated in the framework of the time-dependent density functional theory(TDDFT).The performance of the analytical approaches on the calculation of excited-state energy Hessian has also been shown.It is found that the analytical approaches are superior to the finite-difference method on the computational accuracy and efficiency.The computational cost for a TDDFT excited-state Hessian calculation is only 2–3 times as that for the DFT ground-state Hessian calculation.With the low computational complexity of the developed analytical approaches,it becomes feasible to realize the large-scale numerical calculations on the excited-state vibrational frequencies,vibrational spectroscopies and the electronic-structure parameters which enter the spectrum calculations of electronic absorption and emission,and resonance Raman spectroscopies for medium-to large-sized systems.     

Correlating the vibrational spectra of structurally related molecules: A spectroscopic measure of similarity

Using catastrophe theory and the concept of a mutation path, an algorithm is developed that leads to the direct correlation of the normal vibrational modes of two structurally related molecules. The mutation path is defined by weighted incremental changes in mass and geometry of the molecules in question, which are successively applied to mutate a molecule into a structurally related molecule and thus continuously converting their normal vibrational spectra from one into the other. Correlation diagrams are generated that accurately relate the normal vibrational modes to each other by utilizing mode‐mode overlap criteria and resolving allowed and avoided crossings of vibrational eigenstates. The limitations of normal mode correlation, however, foster the correlation of local vibrational modes, which offer a novel vibrational measure of similarity. It will be shown how this will open new avenues for chemical studies. © 2017 Wiley Periodicals, Inc.     

Vibrational Spectroscopy of the Dehydrogenated Uracil Radical by Autodetachment of Dipole‐Bound Excited States of Cold Anions

Molecules with large enough dipole moments can bind an electron by the dipole field, which has little effect on the molecular core. A molecular anion can be excited to a dipole‐bound state, which can autodetach by vibronic coupling. Autodetachment spectroscopy of a complex anion cooled in a cryogenic ion trap is reported. Vibrational spectroscopy of the dehydrogenated uracil radical is obtained by a dipole‐bound state with partial rotational resolution. Fundamental frequencies for 21 vibrational modes of the uracil radical are reported. The electron affinity of the uracil radical is measured accurately to be 3.4810±0.0006 eV and the binding energy of the dipole‐bound state is measured to be 146±5 cm^?1. The rotational temperature of the trapped uracil anion is evaluated to be 35 K.     

Chemical mechanism of surface‐enhanced raman scattering spectrum of pyridine adsorbed on Ag cluster: Ab initio molecular dynamics approach

The surface‐enhanced Raman scattering (SERS) spectrum of pyridine adsorbed on Ag₂₀ cluster (pyridine‐Ag₂₀) at room temperature is calculated by performing ab initio molecular dynamics simulations in connection with a Fourier transform of the polarizability autocorrelation function to investigate the static chemical enhancement behind the SERS spectrum. The five enhanced vibrational modes of pyridine, namely, υ_6a, υ₁, υ₁₂, υ_9a, and υ_8a, can be assigned and identified by using a new analytical scheme, namely, single‐frequency‐pass filter, which is based on a Fourier transform filtering technique. To understand the factors evoking the enhancement in the SERS spectrum, the dynamic properties of molecular structures and charges for both of the free pyridine and adsorbed pyridine are analyzed. The calculated results indicate that the vibrational amplitudes of adsorbed pyridine are enhanced due to both of the electron transfer from pyridine to Ag₂₀ cluster and the softening of pyridine bond. In addition, the N‐Ag stretching within pyridine‐Ag₂₀ will couple with these five vibrational modes of pyridine. Consequently, the electron transfer between pyridine and Ag₂₀ cluster induced by different molecular vibrational modes prompts the redistribution of electron density of pyridine. These factors collectively cause the noticeable change in polarizability during molecular vibrations and hence result in the enhancement of Raman peaks. © 2013 Wiley Periodicals, Inc.     

Computational spectroscopy of ubiquitin: comparison between theory and experiments

Using the constrained molecular dynamics simulation method in combination with quantum chemistry calculation, Hessian matrix reconstruction, and fragmentation approximation methods, the authors have established computational schemes for numerical simulations of amide I IR absorption, vibrational circular dichroism (VCD), and two-dimensional (2D) IR photon echo spectra of the protein ubiquitin in water. Vibrational characteristic features of these spectra in the amide I vibration region are discussed. From the semiempirical quantum chemistry calculation results on an isolated ubiquitin, amide I local mode frequencies and vibrational coupling constants were fully determined. It turns out that the amide I local mode frequencies of ubiquitin in both gas phase and aqueous solution are highly heterogeneous and site dependent. To directly test the quantitative validity of thus obtained spectroscopic properties, they compared the experimentally measured amide I IR, 2D IR, and electronic circular dichroism spectra with experiments, and found good agreements between theory and experiments. However, the simulated VCD spectrum is just qualitatively similar to the experimentally measured one. This indicates that, due to delicate cancellations between the positive and negative VCD contributions, the prediction of protein VCD spectrum is critically relied on quantitative accuracy of the theoretical model for predicting amide I local mode frequencies. On the basis of the present comparative investigations, they found that the site dependency of amide I local mode frequency, i.e., diagonal heterogeneity of the vibrational Hamiltonian matrix in the amide I local mode basis, is important. It is believed that the present computational methods for simulating various vibrational and electronic spectra of proteins will be of use in further refining classical force fields and in addressing the structure-spectra relationships of proteins in solution.     

Localization transition of instantaneous normal modes and liquid diffusion

We analyze the structure of instantaneous normal modes (INM) associated with the Hessian matrix of a liquid. Utilizing a scaling theory developed in the framework of Anderson localization studies, we unambiguously identify the crossover point in the INM spectrum between extended and localized modes. We establish a relation between the unstable, delocalized INMs and the liquid diffusion coefficient that appears valid over a large temperature range, covering both Arrhenius and non-Arrhenius regimes of temperature dependence. These results suggest a possible route to theoretically relate dynamics to thermodynamical properties of the liquid via the tomography of the INMs.     

Partial hessian fitting for determining force constant parameters in molecular mechanics

We present a new protocol for deriving force constant parameters that are used in molecular mechanics (MM) force fields to describe the bond‐stretching, angle‐bending, and dihedral terms. A 3 × 3 partial matrix is chosen from the MM Hessian matrix in Cartesian coordinates according to a simple rule and made as close as possible to the corresponding partial Hessian matrix computed using quantum mechanics (QM). This partial Hessian fitting (PHF) is done analytically and thus rapidly in a least‐squares sense, yielding force constant parameters as the output. We herein apply this approach to derive force constant parameters for the AMBER‐type energy expression. Test calculations on several different molecules show good performance of the PHF parameter sets in terms of how well they can reproduce QM‐calculated frequencies. When soft bonds are involved in the target molecule as in the case of secondary building units of metal‐organic frameworks, the MM‐optimized geometry sometimes deviates significantly from the QM‐optimized one. We show that this problem is rectified effectively by use of a simple procedure called Katachi that modifies the equilibrium bond distances and angles in bond‐stretching and angle‐bending terms. © 2016 Wiley Periodicals, Inc.     

Modeling the Electronic Absorption Spectra of the Indocarbocyanine Cy3

Accurate modeling of optical spectra requires careful treatment of the molecular structures and vibronic, environmental, and thermal contributions. The accuracy of the computational methods used to simulate absorption spectra is limited by their ability to account for all the factors that affect the spectral shapes and energetics. The ensemble-based approaches are widely used to model the absorption spectra of molecules in the condensed-phase, and their performance is system dependent. The Franck–Condon approach is suitable for simulating high resolution spectra of rigid systems, and its accuracy is limited mainly by the harmonic approximation. In this work, the absorption spectrum of the widely used cyanine Cy3 is simulated using the ensemble approach via classical and quantum sampling, as well as, the Franck–Condon approach. The factors limiting the ensemble approaches, including the sampling and force field effects, are tested, while the vertical and adiabatic harmonic approximations of the Franck–Condon approach are also systematically examined. Our results show that all the vertical methods, including the ensemble approach, are not suitable to model the absorption spectrum of Cy3, and recommend the adiabatic methods as suitable approaches for the modeling of spectra with strong vibronic contributions. We find that the thermal effects, the low frequency modes, and the simultaneous vibrational excitations have prominent contributions to the Cy3 spectrum. The inclusion of the solvent stabilizes the energetics significantly, while its negligible effect on the spectral shapes aligns well with the experimental observations.     

液态RbCl的态密度

给出了在分子动力学模拟基础上Fumi-Tosi势离子液体的正则模式分析方法,用Fumi-Tosi势(包括长程势)代替Lennard-Jones势,并且用等效Coulomb势处理长程Coulomb作用.讨论了Hessian矩阵元的计算方法和Hessian矩阵特征值的计算方法.计算实践表明，取用余误差函数形式的等效库仑势,可以合理地得到Hessian矩阵和态密度.液态RbCl中构型平均态密度的数值结果表明，液态RbCl的态密度表现出与Lennard-Jones液体的态密度相仿的特点.