Shkrob Ilya A, Glover William J, Larsen Ross E, Schwartz Benjamin J
Chemistry Division, Argonne National Laboratory, 9700 South Cass Avenue, Argonne, Illinois 60439, USA.
J Phys Chem A. 2007 Jun 21;111(24):5232-43. doi: 10.1021/jp0682816. Epub 2007 May 27.
Adiabatic mixed quantum/classical (MQC) molecular dynamics (MD) simulations were used to generate snapshots of the hydrated electron in liquid water at 300 K. Water cluster anions that include two complete solvation shells centered on the hydrated electron were extracted from the MQC MD simulations and embedded in a roughly 18 Ax18 Ax18 A matrix of fractional point charges designed to represent the rest of the solvent. Density functional theory (DFT) with the Becke-Lee-Yang-Parr functional and single-excitation configuration interaction (CIS) methods were then applied to these embedded clusters. The salient feature of these hybrid DFT(CIS)/MQC MD calculations is significant transfer (approximately 18%) of the excess electron's charge density into the 2p orbitals of oxygen atoms in OH groups forming the solvation cavity. We used the results of these calculations to examine the structure of the singly occupied and the lower unoccupied molecular orbitals, the density of states, the absorption spectra in the visible and ultraviolet, the hyperfine coupling (hfcc) tensors, and the infrared (IR) and Raman spectra of these embedded water cluster anions. The calculated hfcc tensors were used to compute electron paramagnetic resonance (EPR) and electron spin echo envelope modulation (ESEEM) spectra for the hydrated electron that compared favorably to the experimental spectra of trapped electrons in alkaline ice. The calculated vibrational spectra of the hydrated electron are consistent with the red-shifted bending and stretching frequencies observed in resonance Raman experiments. In addition to reproducing the visible/near IR absorption spectrum, the hybrid DFT model also accounts for the hydrated electron's 190-nm absorption band in the ultraviolet. Thus, our study suggests that to explain several important experimentally observed properties of the hydrated electron, many-electron effects must be accounted for: one-electron models that do not allow for mixing of the excess electron density with the frontier orbitals of the first-shell solvent molecules cannot explain the observed magnetic, vibrational, and electronic properties of this species. Despite the need for multielectron effects to explain these important properties, the ensemble-averaged radial wavefunctions and energetics of the highest occupied and three lowest unoccupied orbitals of the hydrated electrons in our hybrid model are close to the s- and p-like states obtained in one-electron models. Thus, one-electron models can provide a remarkably good approximation to the multielectron picture of the hydrated electron for many applications; indeed, the two approaches appear to be complementary.
采用绝热混合量子/经典(MQC)分子动力学(MD)模拟来生成300 K液态水中水合电子的瞬时图像。从MQC MD模拟中提取出包含以水合电子为中心的两个完整溶剂化壳层的水簇阴离子,并将其嵌入到一个大致为18 Å×18 Å×18 Å的分数点电荷矩阵中,该矩阵旨在代表其余的溶剂。然后将采用Becke-Lee-Yang-Parr泛函的密度泛函理论(DFT)和单激发组态相互作用(CIS)方法应用于这些嵌入的簇。这些混合DFT(CIS)/MQC MD计算的显著特征是,过量电子的电荷密度有相当一部分(约18%)转移到形成溶剂化腔的OH基团中氧原子的2p轨道上。我们利用这些计算结果来研究这些嵌入的水簇阴离子的单占据和较低未占据分子轨道的结构、态密度、可见和紫外吸收光谱、超精细耦合(hfcc)张量以及红外(IR)和拉曼光谱。计算得到的hfcc张量用于计算水合电子的电子顺磁共振(EPR)和电子自旋回波包络调制(ESEEM)光谱,其与碱性冰中捕获电子的实验光谱吻合良好。计算得到的水合电子振动光谱与共振拉曼实验中观察到的红移弯曲和拉伸频率一致。除了再现可见/近红外吸收光谱外,混合DFT模型还解释了水合电子在紫外区的190 nm吸收带。因此,我们的研究表明,为了解释水合电子几个重要的实验观测性质,必须考虑多电子效应:不允许过量电子密度与第一壳层溶剂分子的前沿轨道混合的单电子模型无法解释该物种观测到的磁性、振动和电子性质。尽管需要多电子效应来解释这些重要性质,但在我们的混合模型中,水合电子最高占据和三个最低未占据轨道的系综平均径向波函数和能量与单电子模型中得到的类s和类p态接近。因此,对于许多应用而言,单电子模型可以为水合电子的多电子图像提供非常好的近似;实际上,这两种方法似乎是互补的。