Generalized Bond Polarizability model for more accurate atomistic modeling of Raman spectra

TL;DR

This paper introduces a generalized bond polarizability model (GBPM) that significantly improves the accuracy of Raman spectra predictions from molecular dynamics simulations by including non-traditional atom pairs, aligning well with ab initio results.

Contribution

The paper develops and validates a generalized bond polarizability model that enhances Raman spectra modeling accuracy for various molecules and complex systems.

Findings

GBPM accurately reproduces ab initio polarizability and Raman spectra.

Improved agreement with experimental spectra for liquid water.

Enhanced modeling of complex ferroelectric BaTiO3 systems.

Abstract

Raman spectroscopy is an important tool for studies of molecules, liquids and solids. While Raman spectra can be obtained theoretically from molecular dynamics (MD) simulations, this requires the calculation of the electronic polarizability along the simulation trajectory. First-principles calculations of electronic polarizability are computationally expensive, motivating the development of atomistic models for the evaluation of the changes in the electronic polarizability with the changes in the atomic coordinates of the system. The bond polarizability model (BPM) is one of the oldest and simplest such atomistic models, but cannot reproduce the effects of angular vibrations, leading to inaccurate modeling of Raman spectra. Here, we demonstrate that the generalization of BPM through inclusion of terms for atom pairs that are traditionally considered to be not involved in bonding dramatically improves the accuracy of polarizability modeling and Raman spectra calculations. The generalized BPM (GBPM) reproduces the ab initio polarizability and Raman spectra for a range of tested molecules (SO2, H2S, H2O, NH3, CH4, CH3OH and CH3CH2OH) with high accuracy and also shows significantly improved agreement with ab initio results for the more complex ferroelectric BaTiO3 systems. For liquid water, the anisotropic Raman spectrum derived from atomistic MD simulations using GBPM evaluation of polarizability shows significantly improved agreement with the experimental spectrum compared to the spectrum derived using BPM. Thus, GBPM can be used for the modeling of Raman spectra using large-scale molecular dynamics and provides a good basis for the further development of atomistic polarizability models.