Jezierska Aneta, Panek Jarosław J, Koll Aleksander
University of Wrocław, Faculty of Chemistry, F. Joliot-Curie 14, 50-383 Wrocław, Poland.
Chemphyschem. 2008 Apr 21;9(6):839-46. doi: 10.1002/cphc.200700769.
Car-Parrinello molecular dynamics simulations in vacuum and in the solid state are performed on a strongly anharmonic system, namely, 2-(N-diethylamino-N-oxymethyl)-4,6-dichlorophenol, to investigate its molecular and spectroscopic properties. The investigated compound contains two slightly different molecules in the crystal cell with very short intramolecular hydrogen bonds (of 2.400 and 2.423 A), as determined previously by neutron diffraction. The vibrational properties of the compound are studied on the basis of standard approaches, that is, Fourier transformation of the autocorrelation functions of the atomic velocities and dipole moments. Then, the trajectory obtained from ab initio molecular dynamics is sampled and the obtained snapshots are used to solve the vibrational Schrödinger equations and to calculate the O--H stretching envelope as a superposition of the 0-->1 transition. Using an envelope method, the a posteriori quantum effects are included in the O--H stretching. In addition, NMR spectra are calculated also using the obtained snapshots. One- and two-dimensional potentials of mean force (1D and 2D pmf) are derived to explain the details of the proton dynamics. The computational results are supported by NMR experimental data. In addition, the calculated results are compared with previously published X-ray, neutron diffraction, and spectroscopic descriptions. A detailed analysis of the bridged proton's dynamics is thus obtained. The application of 1D and 2D pmf in a system with a strong bridged-proton delocalization is also demonstrated.
在真空和固态条件下,对强非谐体系2-(N-二乙氨基-N-氧甲基)-4,6-二氯苯酚进行了Car-Parrinello分子动力学模拟,以研究其分子和光谱性质。如先前通过中子衍射所确定的,该研究化合物在晶胞中包含两个略有不同的分子,具有非常短的分子内氢键(键长分别为2.400 Å和2.423 Å)。基于标准方法,即原子速度和偶极矩自相关函数的傅里叶变换,研究了该化合物的振动性质。然后,对从头算分子动力学得到的轨迹进行采样,并使用得到的快照来求解振动薛定谔方程,并将O-H伸缩包络计算为0→1跃迁的叠加。使用包络方法,将后验量子效应纳入O-H伸缩中。此外,还使用得到的快照计算了核磁共振谱。推导了一维和二维平均力势(1D和2D pmf)以解释质子动力学的细节。计算结果得到了核磁共振实验数据的支持。此外,将计算结果与先前发表的X射线、中子衍射和光谱描述进行了比较。由此获得了对桥连质子动力学的详细分析。还展示了1D和2D pmf在具有强桥连质子离域的体系中的应用。