Measuring molecular force fields: Terahertz, inelastic neutron scattering, Raman, FTIR, DFT, and BOMD molecular dynamics of solid l-serine

Abstract

A vibrational analysis of polycrystalline l-serine is provided using experimental terahertz, FTIR, Raman and inelastic neutron scattering (INS) spectra, calculated INS spectra – and Born–Oppenheimer molecular dynamics (BOMD) simulations from which the power spectra for the electronegative elements are compared to the THz spectra. Corrections are made to density functional theory (DFT) calculations for van der Waals interactions. Assignments and potential energy distributions are included for all 3 N = 336 normal modes of an eight molecule supercell, including those for 48 non-bonded whole molecule translating and rotating vibrations, of which three are acoustic modes, usually not considered. Calculated and observed frequencies differ by an average 3 cm −1 ( s = 4). The INS spectrum of these modes below 100 cm −1, calculated from energy second derivatives, show a remarkable similarity to the experimental 10 K spectra. The calculated low frequency modes are insensitive to small changes in cell parameters and geometry. THz intensities are represented by power spectra and not calculated explicitly. Nevertheless, power spectra of 13 ps BOMD trajectories at classical temperatures of 20 K, 400 K, and 500 K are markedly similar to the experimental terahertz spectra at 77 K and 298 K. Calculations on a serine crystal supercell 2 × 2 × 2 molecules deep appear to include, in a crude but fortuitously accurate way, enough of the principle out of phase dispersion to yield a match with experimental frequencies and intensities.