Malek Majid, Danaie Mohammad
Faculty of Electrical and Computer Engineering, Semnan University, Semnan, Iran.
J Mol Model. 2023 Oct 9;29(11):332. doi: 10.1007/s00894-023-05740-z.
In this paper, we simulate a single-molecule diode to calculate the effective coupling and investigate the conductivity, as well as the effect of the electric field on these two parameters. First, we obtain the molecule states and energies at 0 V. The next step is to calculate the electrode/molecule coupling using the obtained electrode and molecule Hamiltonian. The electrode/molecule coupling depends on distance. By increasing the distance from 5 to 5.5 angstroms, the coupling decreases from 0.004 to 0.0002 eV. After calculating the electrode/molecule coupling, which is the most significant parameter in electron transfer, the results can be used to obtain the current-voltage and conductivity curves of the device. Simulation results demonstrate that externally applied electric field to the benzene molecule (isolated molecule) can cause a reduction in the effective coupling between the Au electrode and benzene, leading to decreased current and conductivity. Additionally, the applied electric field narrows the gap between the HOMO and LUMO energy levels.
We conducted this computational work using Gaussian 09 software and a MATLAB code, both of which are based on the density functional theory (DFT) approach and the self-consistent field (SCF) method. For DFT calculations, we employed the three-parameter Beck hybrid exchange functional (B3), hybridized with the nonlocal correlation functional developed by Lee, Yang, and Parr (LYP). All optimizations were performed with triple-zeta polarized (TZP) split-valence 6-311G basis sets. The final step involved calculating the electrode/molecule coupling using the Huckel method and integrating the site-to-state transformation with Huckel parameters and the Fermi golden rule. After this calculation, we obtained the current-voltage and conductivity curves using MATLAB software.
在本文中,我们模拟了一个单分子二极管以计算有效耦合,并研究其电导率以及电场对这两个参数的影响。首先,我们获得了0 V时的分子状态和能量。下一步是使用所得到的电极和分子哈密顿量来计算电极/分子耦合。电极/分子耦合取决于距离。通过将距离从5埃增加到5.5埃,耦合从0.004 eV降低到0.0002 eV。在计算出电极/分子耦合(这是电子转移中最重要的参数)之后,所得结果可用于获得该器件的电流-电压曲线和电导率曲线。模拟结果表明,对苯分子(孤立分子)施加外部电场会导致金电极与苯之间的有效耦合降低,从而导致电流和电导率下降。此外,施加的电场会缩小最高占据分子轨道(HOMO)和最低未占据分子轨道(LUMO)能级之间的差距。
我们使用高斯09软件和一个MATLAB代码进行了这项计算工作,这两者均基于密度泛函理论(DFT)方法和自洽场(SCF)方法。对于DFT计算,我们采用了三参数贝克混合交换泛函(B3),并与由李、杨和帕尔构建的非局域相关泛函(LYP)进行了杂化。所有优化均使用三重ζ极化(TZP)分裂价6 - 311G基组进行。最后一步是使用休克尔方法计算电极/分子耦合,并将位点到状态的变换与休克尔参数和费米黄金规则进行整合。完成此计算后,我们使用MATLAB软件获得了电流-电压曲线和电导率曲线。