Electronic polarization effect on low-frequency infrared and Raman spectra of aprotic solvent: Molecular dynamics simulation study with charge response kernel by second order Møller–Plesset perturbation method

Ambiguities in Decomposing Molecular Polarizability into Atomic Charge Flow and Induced Dipole Contributions

The molecular dipole polarizability can be decomposed into components corresponding to the charge flow between atoms and changes in atomic dipole moments. Such decompositions are recognized to depend on how atoms are defined within a molecule, as, for example, by Hirshfeld, iterative Stockholder, or quantum topology partitioning of the electron density. For some of these, however, there are significant differences between the numerical results obtained by analytical response methods and finite field calculations. We show that this difference is due to analytical response methods accounting for (only) the change in electron density by a perturbation, while finite field methods may also include a component corresponding to a perturbation-dependent change in the definition of an atom within a molecule. For some atom-in-molecule definitions, such as the iterative Hirshfeld, iterative Stockholder, and quantum topology methods, the latter effect significantly increases the charge flow component. The decomposition of molecular polarizability into atomic charge flow and induced dipole components thus depends on whether the atom-in-molecule definition is taken to be perturbation-dependent.

Dynamic construction of refractive index-dependent vibrations using surface plasmon-phonon polaritons

One of the fundamental hurdles in infrared spectroscopy is the failure of molecular identification when their infrared vibrational fingerprints overlap. Refractive index (RI) is another intrinsic property of molecules associated with electronic polarizability, but with limited contribution to molecular identification in mixed environments currently. Here, we investigate the coupling mode of localized surface plasmon and surface phonon polaritons for vibrational de-overlapping. The coupling mode is sensitive to the molecular refractive index, attributed to the RI-induced vibrational variations of surface phonon polaritons (SPhP) within the Reststrahlen band, referred to as RI-dependent SPhP vibrations. The RI-dependent SPhP vibrations are linked to molecular RI features. According to the deep-learning-augmented demonstration of bond-breaking-bond-making dynamic profiling in biological reaction, we substantiate that the RI-dependent SPhP vibrations effectively disentangle overlapping vibrational modes, achieving a 92% identification accuracy even for the strongly overlapping vibrational modes in the reaction. Our findings offer insights into the realm of light-matter interaction and provide a valuable toolkit for biomedicine applications.

Unifying Charge-Flow Polarization Models

We show that several models where electric polarization in molecular systems is modeled by charge-flow between atoms can all be considered as different manifestations of a general underlying mathematical structure. The models can be classified according to whether they employ atomic or bond parameters and whether they employ atom/bond hardness or softness. We show that an ab initio calculated charge response kernel can be considered as the inverse screened Coulombic matrix projected onto the zero-charge subspace, and this may provide a method for deriving charge screening functions to be used in force fields. The analysis suggests that some models contain redundancies, and we argue that a parameterization of charge-flow models in terms of bond softness is preferable as it depends on local quantities and decay to zero upon bond dissociation, while bond hardness depends on global quantities and increases toward infinity upon bond dissociation.

Using atomic charges to model molecular polarization

We review different models for introducing electric polarization in force fields, with special focus on methods where polarization is modelled at the atomic charge level. While electric polarization has been included in several force fields, the common approach has been to focus on atomic dipole polarizability. Several approaches allow modelling electric polarization by using charge-flow between charge sites instead, but this has been less exploited, despite that atomic charges and charge-flow is expected to be more important than atomic dipoles and dipole polarizability. A number of challenges are required to be solved for charge-flow models to be incorporated into polarizable force fields, for example how to parameterize the models and how to make them computational efficient.

A wave-function based approach for polarizable charge model: Systematic comparison of polarization effects on protic, aprotic, and ionic liquids

We first describe a wave-function based formalism of polarizable charge model by starting from the Hartree product ansatz for the total wave function and making the second-order expansion of individual molecular energies with the use of partial charge operators. The resulting model is shown to be formally equivalent to the charge response kernel model that starts from the linear-response approximation to partial charges, and also closely related to a family of fluctuating charge models that are based on the electronegativity equalization principle. We then apply the above model to a systematic comparison of polarization effects on qualitatively different liquids, namely, protic solvents (water and methanol), an aprotic polar solvent (acetonitrile), and imidazolium-based ionic liquids. Electronic polarization is known to decelerate molecular motions in conventional solvents while it accelerates them in ionic liquids. To obtain more insights into these phenomena, we consider an effective decomposition of total polarization energy into molecular contributions, and show that their statistical distribution is well-correlated with the acceleration/deceleration of molecular motions. In addition, we perform effective nonpolarizable simulations based on mean polarized charges, and compare them with fully polarizable simulations. The result shows that the former can reproduce structural properties of conventional solvents rather accurately, while they fail qualitatively to reproduce acceleration of molecular motions in ionic liquids.

Polarizable Force Field for Protein with Charge Response Kernel

We present a molecular mechanical force field for polypeptides and proteins involving the electronic polarization effect described with the charge response kernel. All of the electrostatic parameters for 20 amino acids are obtained by ab initio electronic structure calculations and combined with the AMBER99 force field. The refittings of dihedral angle parameters in the torsional potentials are performed so as to reproduce the ab initio optimized geometries and relative energies for the conformers of dipeptides. The present force field is applied to molecular dynamics simulation calculations of the extended alanine tetra and cyclic pentapeptides in aqueous solution. The infrared spectra are calculated in order to analyze the charge polarization effect on the spectral profiles.