Kerisit Sebastien, Rosso Kevin M
Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, Washington 99352, USA.
J Chem Phys. 2005 Dec 8;123(22):224712. doi: 10.1063/1.2137319.
Molecular-dynamics simulations and ab initio electronic structure calculations were carried out to determine the rate of charge transfer in stoichiometric wustite (FeO). The charge transfer of interest occurs by II/III valence interchange between nearest-neighbor Fe atoms, with the Fe(III) constituting a "hole" electronic defect. There are two possible nearest-neighbor charge transfers in the FeO lattice, which occur between edge-sharing or corner-sharing FeO(6) octahedra. Molecular-dynamics simulations predict charge-transfer rates of 3.7 x 10(11) and 1.9 x 10(9) s(-1) for the edge and corner transfers, respectively, in good agreement with those calculated using an ab initio cluster approach (1.6 x 10(11) and 8.0 x 10(8) s(-1), respectively). The calculated rates are also similar to those along the basal and c-axis directions in hematite (alpha-Fe(2)O(3)) determined previously. Therefore, as is the case for hematite, wustite is predicted to show anisotropic electrical conductivity. Our findings indicate that a rigid-ion model does not give acceptable results, thus showing the need to account for the change in polarizability of the system upon charge transfer. Our model achieves this by using a simple mechanical shell model. By calculating the electronic coupling matrix elements for many transition state configurations obtained from the molecular-dynamics simulations, we found evidence that the position of the bridging oxygen atoms can greatly affect the amount of electronic coupling between the donor and acceptor states. Finally, we address the effect of oxygen vacancies on the charge transfer. It was found that an oxygen vacancy not only creates a driving force for holes to transport away from the vacancy (or equivalently for electrons to diffuse toward the vacancy) but also lowers the free-energy barriers for charge transfer. In addition, the reorganization energy significantly differed from the nondefective case in a small radius around the defect.
进行了分子动力学模拟和从头算电子结构计算,以确定化学计量比的方铁矿(FeO)中的电荷转移速率。所关注的电荷转移是通过最近邻Fe原子之间的II/III价互换发生的,其中Fe(III)构成一个“空穴”电子缺陷。在FeO晶格中有两种可能的最近邻电荷转移,它们发生在边共享或角共享的FeO(6)八面体之间。分子动力学模拟预测边转移和角转移的电荷转移速率分别为3.7×10¹¹和1.9×10⁹ s⁻¹,这与使用从头算团簇方法计算的结果(分别为1.6×10¹¹和8.0×10⁸ s⁻¹)非常吻合。计算得到的速率也与之前确定的赤铁矿(α-Fe₂O₃)中沿基面和c轴方向的速率相似。因此,与赤铁矿的情况一样,预计方铁矿会表现出各向异性的电导率。我们的研究结果表明,刚性离子模型无法给出可接受的结果,从而表明需要考虑电荷转移时系统极化率的变化。我们的模型通过使用简单的机械壳模型实现了这一点。通过计算从分子动力学模拟获得的许多过渡态构型的电子耦合矩阵元,我们发现有证据表明桥连氧原子的位置会极大地影响供体和受体态之间的电子耦合量。最后,我们讨论了氧空位对电荷转移的影响。发现氧空位不仅为载流子提供了从空位迁移走的驱动力(或者等效地为电子向空位扩散提供驱动力),而且还降低了电荷转移的自由能垒。此外,在缺陷周围的小半径范围内,重组能与无缺陷情况有显著差异。