Glover William J, Larsen Ross E, Schwartz Benjamin J
Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, USA.
J Chem Phys. 2008 Oct 28;129(16):164505. doi: 10.1063/1.2996350.
The charge-transfer-to-solvent (CTTS) reactions of solvated atomic anions serve as ideal models for studying the dynamics of electron transfer: The fact that atomic anions have no internal degrees of freedom provides one of the most direct routes to understanding how the motions of solvent molecules influence charge transfer, and the relative simplicity of atomic electronic structure allows for direct contact between theory and experiment. To date, molecular dynamics simulations of the CTTS process have relied on a single-electron description of the atomic anion-only the electron involved in the charge transfer has been treated quantum mechanically, and the electronic structure of the atomic solute has been treated via pseudopotentials. In this paper, we examine the severity of approximating the electronic structure of CTTS anions with a one-electron model and address the role of electronic exchange and correlation in both CTTS electronic structure and dynamics. To do this, we perform many-electron mixed quantum/classical molecular dynamics simulations of the ground- and excited-state properties of the aqueous sodium anion (sodide). We treat both of the sodide valence electrons quantum mechanically and solve the Schrodinger equation using configuration interaction with singles and doubles (CISD), which provides an exact solution for two electrons. We find that our multielectron simulations give excellent general agreement with experimental results on the CTTS spectroscopy and dynamics of sodide in related solvents. We also compare the results of our multielectron simulations to those from one-electron simulations on the same system [C. J. Smallwood et al., J. Chem. Phys. 119, 11263 (2003)] and find substantial differences in the equilibrium CTTS properties and the nonadiabatic relaxation dynamics of one- and two-electron aqueous sodide. For example, the one-electron model substantially underpredicts the size of sodide, which in turn results in a dramatically different solvation structure around the ion. The one-electron model also misses the existence of an entire manifold of bound CTTS excited states and predicts an absorption spectrum that is blueshifted from that in the two-electron model by over 2 eV. Even the use of a quantum mechanics/molecular mechanics (QM/MM)-like approach, where we calculated the electronic structure with our CISD method using solvent configurations generated from the one-electron simulations, still produced an absorption spectrum that was shifted approximately 1 eV to the blue. In addition, we find that the two-electron model sodide anion is very polarizable: The instantaneous dipole induced by local fluctuating electric fields in the solvent reaches values over 14 D. This large polarizability is driven by an unusual solvation motif in which the solvent pushes the valence electron density far enough to expose the sodium cation core, a situation that cannot be captured by one-electron models that employ a neutral atomic core. Following excitation to one of the bound CTTS excited states, we find that one of the two sodide valence electrons is detached, forming a sodium atom:solvated electron contact pair. Surprisingly, the CTTS relaxation dynamics are qualitatively similar in both the one- and two-electron simulations, a result we attribute to the fact that the one-electron model does correctly describe the symmetry of the important CTTS excited states. The excited-state lifetime of the one-electron model, however, is over three times longer than that in the two-electron model, and the detachment dynamics in the two-electron model is correlated with the presence of solvent molecules that directly solvate the cationic atomic core. Thus, our results make it clear that a proper treatment of anion electron structure that accounts for electronic exchange and correlation is crucial to understanding CTTS electronic structure and dynamics.
溶剂化原子阴离子的电荷转移到溶剂(CTTS)反应是研究电子转移动力学的理想模型:原子阴离子没有内部自由度这一事实为理解溶剂分子的运动如何影响电荷转移提供了最直接的途径之一,并且原子电子结构的相对简单性使得理论与实验能够直接联系起来。迄今为止,CTTS过程的分子动力学模拟依赖于对原子阴离子的单电子描述——只有参与电荷转移的电子被量子力学处理,而原子溶质的电子结构则通过赝势来处理。在本文中,我们研究了用单电子模型近似CTTS阴离子电子结构的严重程度,并探讨了电子交换和关联在CTTS电子结构及动力学中的作用。为此,我们对水合钠离子(钠化物)的基态和激发态性质进行了多电子混合量子/经典分子动力学模拟。我们对钠化物的两个价电子都进行量子力学处理,并使用单双激发组态相互作用(CISD)求解薛定谔方程,该方法为两个电子提供了精确解。我们发现,我们的多电子模拟在CTTS光谱和钠化物在相关溶剂中的动力学方面与实验结果总体上非常吻合。我们还将多电子模拟的结果与同一系统的单电子模拟结果[C. J. Smallwood等人,《化学物理杂志》119, 11263 (2003)]进行了比较,发现单电子和双电子水合钠化物在平衡CTTS性质和非绝热弛豫动力学方面存在显著差异。例如,单电子模型大大低估了钠化物的尺寸,这反过来又导致离子周围的溶剂化结构有很大不同。单电子模型还忽略了一整套束缚CTTS激发态的存在,并预测的吸收光谱比双电子模型中的吸收光谱蓝移超过2 eV。即使使用类似量子力学/分子力学(QM/MM)的方法,即我们使用CISD方法根据单电子模拟生成的溶剂构型来计算电子结构,仍然产生了大约蓝移1 eV的吸收光谱。此外,我们发现双电子模型的钠化物阴离子非常容易极化:溶剂中局部波动电场诱导的瞬时偶极矩达到超过14 D的值。这种大的极化率是由一种不寻常的溶剂化模式驱动的,在这种模式中,溶剂将价电子密度推得足够远以暴露出钠离子核心,而采用中性原子核的单电子模型无法捕捉到这种情况。在激发到束缚CTTS激发态之一后,我们发现钠化物的两个价电子之一被分离,形成钠原子:溶剂化电子接触对。令人惊讶的是,单电子和双电子模拟中的CTTS弛豫动力学在定性上是相似的,我们将这一结果归因于单电子模型确实正确地描述了重要CTTS激发态的对称性。然而,单电子模型的激发态寿命比双电子模型长三倍多,并且双电子模型中的分离动力学与直接溶剂化阳离子核心的溶剂分子的存在相关。因此,我们的结果清楚地表明,正确处理考虑电子交换和关联的阴离子电子结构对于理解CTTS电子结构和动力学至关重要。