Hydration structure and dynamics, ultraviolet–visible and fluorescence spectra of caffeine in ambient liquid water. A combined classical molecular dynamics and quantum chemical study

Abstract

The hydration structure and related dynamics of caffeine diluted in ambient liquid water have been extensively studied by performing classical molecular dynamics simulations, using our previously developed potential model of caffeine and the TIP4P/2005 water model. The results obtained have revealed that the first hydration shell of caffeine contains on average 56 water molecules. Hydrated caffeine forms in total 3 hydrogen bonds with its neighbor water molecules, with the Ocaffeine … Hwater hydrogen bonds exhibiting similar lifetimes with the ones corresponding to Owater … Hwater hydrogen bonds in liquid water. The self-diffusion coefficient of caffeine has been found to be four times lower than the corresponding value for water, being also in agreement with recent experimental measurements. The presence of water molecules inside the solvation shell of caffeine changes significantly their low-frequency intermolecular vibrations, as reflected on the calculated atomic velocity time correlation functions and corresponding spectral densities. Using the estimated average intermolecular structure of the first hydration shell of caffeine, the molecular cluster caffeine@W56 was optimized via quantum chemical calculations and subsequently the time-dependent density functional theory was used in order to predict the ultraviolet–visible and fluorescence spectra of hydrated caffeine. The results obtained are in agreement with recent experimental studies, which have proposed that such spectroscopic measurements can be used for the direct determination of alkaloids in aqueous extracts of natural products. In this framework, multi-scale molecular modelling providing accurate predictions of experimental data could also be a very useful tool, linking theoretical physical chemistry with analytical chemistry applications.