文摘
The two-dimensional infrared spectrum of NMAH and NMAD in H2O and D2O is computed on the basis of force field parametrizations ranging from standard point charge (PC) to more elaborate multipolar (MTP) representations of the electrostatics. For the latter, the nonbonded parameters (MTP and van der Waals) were optimized to reproduce thermodynamic data. The frequency trajectory and frequency鈥揻requency correlation function (FFCF) are determined from explicit frequency calculations on 106 snapshots without using a more traditional 鈥渕apping鈥?approach. This allows us to both sample configurations and compute observables in a consistent fashion. In agreement with experiment, the FFCF shows one very rapid time scale (in the 50 fs range) followed by one or two longer time scales. In the case of three time scales, the intermediate one is 鈮?.5 ps or shorter, whereas the longest time scale can extend up to 2 or 3 ps. All interaction models lead to three time scales in the FFCF when fitted to an empirical parametrized form. When two time scales are assumed鈥攁s is usually done in the analysis of experimental data鈥攁nd the short time scale is fixed to the 蟿1 = 50鈥?00 fs range, the correlation time 蟿c from the simulations ranges from 0.7 to 1 ps, which agrees quite well with experimentally determined values. The major difference between MTP and PC models is the observation that the later decay times in the FFCF are longer for simulations with MTPs. Also, the amplitude of the FFCF is reduced when simulations are carried out with MTPs. Overall, however, PC-based models perform well compared to those based on MTPs for NMAD in D2O and can be recommended for such investigations in the context of peptide and protein simulations.