Abstract
The paper illustrates the Activity Weighted Velocities (AWV) methodology to compute Vibrational Circular Dichroism (VCD) anharmonic spectra from Density Functional Theory (DFT) molecular dynamics. AWV calculates the spectra by the Fourier Transform of the time correlation functions of velocities, weighted by specific observables: the Atomic Polar Tensors (APTs) and the Atomic Axial Tensors (AATs). Indeed, AWV shows to correctly reproduce the experimental spectra for systems in the gas and liquid phases, both in the case of weakly and strongly interacting systems. The comparison with the experimental spectra is striking especially in the fingerprint region, as demonstrated by the three benchmark systems discussed: (1S)-Fenchone in the gas phase, (S)-(-)-Propylene oxide in the liquid phase, and (R)-(-)-2-butanol in the liquid phase. The time evolution of APTs and AATs can be adequately described by a linear combination of the tensors of a small set of appropriate reference structures, strongly reducing the computational cost without compromising accuracy. Additionally, AWV allows the partition of the spectral signal in its molecular components without any expensive postprocessing and any localization of the charge density or the wave function.