Sharma Deepika, Das Banshi, Chandra Amalendu
Department of Chemistry, Indian Institute of Technology, Kanpur, Uttar Pradesh 208016, India.
J Phys Chem B. 2023 Aug 3;127(30):6714-6725. doi: 10.1021/acs.jpcb.3c03241. Epub 2023 Jul 21.
The terahertz (THz) absorption spectrum of water is calculated at varying temperatures from 260 to 340 K from molecular dynamics simulations using a polarizable potential model of water. The current calculations produce the known experimental features of the THz spectrum of water such as the hydrogen-bond stretch mode at ∼200 cm and librational mode at ∼600 cm. The peak positions generally show a red shift with an increase of temperature due to the weakening of the hydrogen bonds at higher temperatures. Overall, the changes of the spectrum with temperature are found to be in good agreement with experimental results. The total THz spectrum at a given temperature is dissected into self- and cross-correlation contributions and also into contributions from permanent dipoles, induced dipoles, and permanent-induced dipole correlations. It is shown that while the peak at ∼200 cm due to hydrogen-bond stretching primarily comes from fluctuations of induced dipoles, the librational peak at around 600 cm originates mainly from fluctuations of the permanent dipoles. Also, through calculations of self- and cross-correlations, it is shown that the broad librational peak arises from the superposition of several components like the antisymmetric libration, symmetric libration, and also self-dipole correlations. The length-scale-resolved calculations of cross-correlations reveal the contributions from different solvation shells to the total cross-component of the THz spectrum and how such length-scale-resolved components change with temperature. Results are also presented for the dielectric relaxation of water over different length scales and temperatures.
利用水的极化势模型,通过分子动力学模拟计算了260至340K不同温度下水的太赫兹(THz)吸收光谱。当前的计算结果展现出了水的太赫兹光谱已知的实验特征,比如在200cm处的氢键拉伸模式以及在600cm处的摆动模式。由于高温下氢键减弱,峰值位置通常随温度升高而出现红移。总体而言,发现光谱随温度的变化与实验结果高度吻合。给定温度下的总太赫兹光谱被分解为自相关和交叉相关贡献,也被分解为永久偶极子、诱导偶极子以及永久 - 诱导偶极子相关性的贡献。结果表明,虽然~200cm处由氢键拉伸引起的峰值主要来自诱导偶极子的涨落,但600cm左右的摆动峰值主要源于永久偶极子的涨落。此外,通过自相关和交叉相关计算表明,宽摆动峰源于反对称摆动、对称摆动以及自偶极子相关性等多个分量的叠加。交叉相关的长度尺度分辨计算揭示了不同溶剂化壳层对太赫兹光谱总交叉分量的贡献,以及这种长度尺度分辨分量如何随温度变化。还给出了不同长度尺度和温度下水的介电弛豫结果。