Simulations of the OKE Response in Simple Liquids Using a Polarizable and a Nonpolarizable Force Field

Recently polarizable force fields are becoming increasingly popular for molecular dynamics simulations. As the signal obtained in the optical Kerr effect (OKE) experiment is due to the polarizability dynamics of the investigated system, a study is conducted in order to compare the experimental results with those obtained with the polarizable AMOEBA force field. The comparison is made in the frequency domain; however, time domain data are also included. The selected molecular systems are the isotropic carbon tetrachloride molecule, the anisotropic chloroform, carbon disulfide and acetone molecules, and the hydrogen-bonded water and methanol molecules. Different dipole-induced-dipole (DID) method variants are used for calculation of the OKE response, showing the importance of use of the all-atom approach with preoptimized atomic polarizabilities. In order to obtain a good intermolecular to intramolecular components amplitude ratio, the isotropic polarizability in the Thole correction needs to be updated between iterations. The convergence of the spectra calculated with different DID variants is also considered, and the approach that appears to be the best gives a very good approximation after three iterations. The comparison of the experimental and simulated spectra shows a rather good agreement for the non-hydrogen-bonded molecules, although the contribution of the reorientation of anisotropic molecules is overestimated. In the case of the hydrogen-bonded molecules, the theoretical spectra are far from the experimental ones. The highly overestimated librational bands indicate excessive polarizability anisotropy introduced by the potential model. Finally, in order to verify the significance of different components of the AMOEBA model, it is gradually simplified and compared with a simple reference potential model. Removal of polarizability shows a tremendous change in the case of hydrogen-bonded liquids, whereas for the other molecules it is of minor importance. The non-hydrogen-bonded liquids are, however, more sensitive to the presence of atomic multipoles in the model.