Theoretical Study of Endohedral Fullerenes M@C60 (M = Li,Na, or K) in Periodic Boundary Conditions

Journal of Experimental and Theoretical Physics - Density functional calculations in periodic boundary conditions (PBCs) were performed to investigate the structural and electronic properties of...     

内嵌复合物M@C60（M=Li,Na,K,Rb,Cs）的分子力学研究

利用分子力学方法计算了碱金属内嵌复合物Ｍ＠Ｃ６０中Ｍ与Ｃ６０之间的相互作用，考察了Ｍ在Ｃ６０笼内的平衡位置。研究表明Ｌｉ和Ｎａ的平衡位置偏离Ｃ６０分子的中心，Ｋ、Ｒｂ和ＣＳ的平衡位置在Ｃ６０分子的中心。平衡位置的确定取决于色散作用和排斥作用的大小。最后，讨论了碱金属原子进入笼内的可能机制。     

Theoretical investigations on the orientational dependence of electron transport through porphyrin molecular wire

The effect of molecular orientation on the electron transport behavior of single porphyrin sandwiched between two gold (111) electrodes is investigated by density functional theory calculations combined with non-equilibrium Green’s function method. The results show that the porphyrin with parallel connection to gold (111) electrodes is more conductive than the porphyrin with diagonal connection to gold (111) electrodes. The mechanism of the difference of electron transport for these two molecular junctions is analyzed from the transmission spectra and the molecular projected self-consistent Hamiltonian states. It is found that the intrinsic nature of the molecule, such as the π-conjugated framework and the strength of molecule–electrode coupling, are the essential reason for generating this difference of electron transport for the two molecular systems.     

Structural and optical properties of the M@C59X cages (X=N,B and M=Li,Na)

Using B3LYP/6-31G* density functional level of theory, the structural and optical properties of the C₆₀ and M@C₅₉X cages have been investigated. Results indicate that the charge on C atoms and band gap of C₆₀ cage are changed dramatically with the substitution of one B or N atom at one of the C sites and the Li and Na atom encapsulations in the C₆₀ cage. The Mulliken analyses show that the charge is transferred completely between the alkali atoms and the C₅₉X cage. The substitutional and encapsulation doping (SED) reduce the optical gaps of the C₆₀ cage. Also, the oscillator strengths of the absorption peaks are dependent on dopant types.     

Nonlinear optical and structural properties of M@C N endohedrals (M = Li,Ca and Sc,N = 60 and 70)

We have investigated the static (ω = 0) and frequency-dependent nonlinear optical (NLO) properties of the M@C _N endohedrals (M = Li, Ca, Sc, N = 60 and 70) using the SSH (Su–Schrieffer–Heeger) approximation and sum-over-state (SOS) approach. Also, we study the effects of displacement and alkali, alkaline earth and lanthanide metal atoms and type of cage on the hyperpolarisabilities of the M@C _N endohedrals. The hyperpolarisability magnitudes and spectra are in agreement with experiment and the work of others using the SSH approximation. Our results indicate that the cage-type effect on the NLO spectra of M@C _N endohedrals is dramatic. Also, atom type has little effect on the highest peak value. These relationships between the atom and cage type and hyperpolarisability values may be beneficial to experimentalists when designing new NLO materials with large NLO responses.     

The compressive mechanical properties of Cn (n=20, 60, 80, 180) and endohedral M@C60 (M=Na,Al, Fe) fullerene molecules

The compressive mechanical properties of C_n (n = 20, 60, 80, 180) and endohedral M@C₆₀ (M = Na, Al, Fe) fullerene molecules are investigated using a quantum molecular dynamics (QMD) technique. Energy–strain curves, force–strain curves, endurance load, failure strain corresponding to the endurance load, and compressive stiffness of the fullerene molecules are obtained. The compressive mechanical properties of C₂₀, C₆₀, C₈₀, C₁₈₀ and M@C₆₀ (M = Na, Al, Fe) are discussed. The results show that the larger the magic number n of an empty fullerene, the higher its endurance load and compressive stiffness, but the lower its failure strain, and comparing to the empty C₆₀ fullerene, all the M@C60 molecules have greater endurance capability and failure strain.     

B与N掺杂对单层石墨纳米带自旋极化输运的影响

通过第一性原理计算研究了具有锯齿状边沿并且具有反铁磁构型的单层石墨纳米带的自旋极化输运.研究发现,在中心散射区同一位置掺入单个B和N原子,尽管对整个体系磁矩的影响完全相同,但对两个自旋分量电流的影响却完全相反.掺B时,自旋向上的电流显著大于自旋向下的电流;而掺N时,自旋向下的电流显著大于自旋向上的电流.这是由于不管掺B还是掺N都将打破自旋简并,使得导带和价带中自旋向上的能级比自旋向下的能级更高.掺B引入空穴,使完全占据的价带变为部分占据,从而自旋向上的能级正好处于费米能级,使得电子透射能力更强、电流更大,而自旋向下的能级则离费米能级较远使电子透射的能力较弱.掺N则引入电子,使得原来全空的导带变为部分占据,从而费米能级穿过导带中自旋向下的能级,使得自旋向下的电子比自旋向上的电子透射能力更强. 关键词： 自旋极化输运 单层石墨纳米带 第一性原理 非平衡格林函数     

密度泛函理论计算内掺富勒烯M@C60H60(M=Li、Na)的几何结构和电子性质

采用密度泛函理论中的广义梯度近似(generalized gradient approximation,简称GGA),对M@C_(60)H_(60) (M=Li、Na)几何结构和电子性质进行计算研究.发现M原子的平衡位置处在偏心位置处,并且稳定的存在于一个围绕中心的球体内;掺杂能计算表明:M@C_(60)H_(60)需要在一定的实验条件下才能被合成出来;电子性质分析表明:M原子掺入到C_(60)H_(60)中,对费米能级附近有一定贡献,并产生了1 μB的净磁矩.     

The spin-dependent electronic transport properties of M(dcdmp)2 (M = Cu,Au, Co,Ni) molecular devices based on zigzag graphene nanoribbon electrodes

The spin-dependent electronic transport properties of M(dcdmp)₂ (M = Cu, Au, Co, Ni; dcdmp = 2,3-dicyano-5,6-dimercaptopyrazyne) molecular devices based on zigzag graphene nanoribbon (ZGNR) electrodes were investigated by density functional theory combined nonequilibrium Green's function method (DFT-NEGF). Our results show that the spin-dependent transport properties of the M(dcdmp)₂ molecular devices can be controlled by the spin configurations of the ZGNR electrodes, and the central 3d-transition metal atom can introduce a larger magnetism than that of the nonferrous metal one. Moreover, the perfect spin filtering effect, negative differential resistance, rectifying effect and magnetic resistance phenomena can be observed in our proposed M(dcdmp)₂ molecular devices.     

Effect of the encapsulation of Li atom on the electronic transport properties of C20F20 cage

Carrying out theoretical calculations using a self-consistent ab initio approach that combines the non-equilibrium Green′s function formalism with density functional theory, we investigate the effect of the center encapsulation of Li atom on the electronic transport properties of C₂₀F₂₀ cage sandwiched between two bulk gold electrodes. The results show that the electrical conductivity of the endohedral complex Li@C₂₀F₂₀ becomes better than that of the empty C₂₀F₂₀ in the bias voltages ranging from 0 to 1.2 V. The novel negative differential resistance behavior in the I-V characteristic curves can be observed by inserting Li atom into C₂₀F₂₀ cage. The mechanism for the negative differential resistance behavior of Li@C₂₀F₂₀ is suggested.     

Contact position and width effect of graphene electrode on the electronic transport properties of dehydrobenzoannulenne molecule under bias

By applying nonequilibrium Green?s function formalism in combination with density functional theory, we have investigated the electronic transport properties of dehydrobenzoannulenne molecule attached to different positions of the zigzag graphene nanoribbons (ZGNRs) electrode. The different contact positions are found to drastically turn the transport properties of these systems. The negative differential resistance (NDR) effect can be found when the ZGNRs electrodes are mirror symmetry under the xz midplane, and the mechanism of NDR has been explained. Moreover, parity limitation tunneling effect can be found in a certain symmetry two-probe system and it can completely destroy electron tunneling process. The present findings might be useful for the application of ZGNRs-based molecular devices.     

M+AsF6-(M=Li,Na,K,Rb和Cs)离子对的结构和振动光谱

利用abinitio计算方法研究了AsF6-阴离子和M AsF6-(M =Li ,Na ,K ,Rb 和Cs )直接接触离子对的结构和光谱行为。结果表明,具有C3V结构三齿相互作用的M AsF6-最稳定,二齿配位的结构只有在特定条件下才可能稳定存在。当形成离子对后,阳离子与AsF6-的相互作用将改变AsF6-的结构,这其中Li 的影响最大。另外,AsF6-光谱的变化可用来指认溶液中离子的缔合和离子对的种类。     

Effect of boundary treatments on quantum transport current in the Green’s function and Wigner distribution methods for a nano-scale DG-MOSFET

In this paper, we conduct a study of quantum transport models for a two-dimensional nano-size double gate (DG) MOSFET using two approaches: non-equilibrium Green’s function (NEGF) and Wigner distribution. Both methods are implemented in the framework of the mode space methodology where the electron confinements below the gates are pre-calculated to produce subbands along the vertical direction of the device while the transport along the horizontal channel direction is described by either approach. Each approach handles the open quantum system along the transport direction in a different manner. The NEGF treats the open boundaries with boundary self-energy defined by a Dirichlet to Neumann mapping, which ensures non-reflection at the device boundaries for electron waves leaving the quantum device active region. On the other hand, the Wigner equation method imposes an inflow boundary treatment for the Wigner distribution, which in contrast ensures non-reflection at the boundaries for free electron waves entering the device active region. In both cases the space-charge effect is accounted for by a self-consistent coupling with a Poisson equation. Our goals are to study how the device boundaries are treated in both transport models affects the current calculations, and to investigate the performance of both approaches in modeling the DG-MOSFET. Numerical results show mostly consistent quantum transport characteristics of the DG-MOSFET using both methods, though with higher transport current for the Wigner equation method, and also provide the current–voltage (I–V) curve dependence on various physical parameters such as the gate voltage and the oxide thickness.     

Transport through a triple quantum dot system: Formation of resonance band and its application as a spin filter

A three-quantum-dot spin filter based on nonequilibrium Green?s function technique is proposed with external magnetic flux, Rashba spin orbit interaction, and intradot coulomb interaction taken into consideration. Numerical results indicate a spin filter can be made efficient by adjusting external magnetic flux and Rashba spin orbit interaction. Moreover, the formation of a resonance band is discussed through calculation. It is observed that the possibility of transition from one peak to other three peaks in the conductance spectrum increases with increasing interdot coupling strength.     

Elastic properties of MBH4 (M=Na, K, Rb, Cs)

The elastic properties of MBH₄ (M=Na, K, Rb, Cs) are studied by first-principles calculations using a projected augmented plane-wave method and finite strain method. KBH₄ is found to be the most isotropic material among MBH₄, which can be attributed to the contribution of C₄₄. The Debye temperatures are also calculated using theoretical elastic constants and the change along the series is explained.     

Electronic structures of efficient MBiO_-3(M= Li,Na, K,Ag) photocatalyst

In order to deepen the understanding of the relationship between fundamental properties(including: microstructure and composition) and photocatalytic performance, four bismuthate compounds, including: LiBiO_3, NaBiO_3, KBiO_3, and AgBiO_3, are regarded as research examples in the present work, because they have particular crystal structures and similar compositions. Using density functional theory calculations, their structural, electronic, and optical properties are investigated and compared systematically. First of all, the calculated results of crystal structures and optical properties are in agreement with available published experimental data. Based on the calculated results, it is found that the tunneled or layered micro-structural properties lead to the stronger interaction between bismuth and oxygen, and the weaker interaction between alkaline-earth metal and [BiO_6] octahedron, resulting in the feature of multi-band gaps in the cases of LiBiO_3,NaBiO_3, and KBiO_3. This conclusion is supported by the case of AgBiO_3, in which the feature of multi-band gaps disappears, due to the stronger interaction between the noble metal and [BiO_6] octahedron. These properties have significant advantages in the photocatalytic performance: absorbing low energy photons, rapidly transferring energy carriers. Furthermore, the features of electronic structures of bismuthate compounds are well reflected by the absorption spectra, which could be confirmed by experimental measurements in practice. Combined with the calculated results, it could be considered that the crystal structures and compositions of the photocatalyst determine the electronic structures and optical properties,and subsequently determine the corresponding photocatalytic performance. Thus, a novel Bi-based photocatalyst driven by visible-light could be designed by utilizing specific compositions to form favorable electronic structures or specific micro-structures to form a beneficial channel for energy carriers.     

内嵌富勒烯[M@C_(70)](M=Li,Na,K,Be,Mg,Ca)与[Mg(H_2O)_6]~(2+)相互作用的密度泛函理论研究

使用密度泛函理论方法,对内嵌多种碱金属及碱土金属原子(M=Li,Na,K,Be,Mg,Ca)的富勒烯C70与水合Mg(II)离子之间的相互作用进行了理论研究.首先对各原子(M=Li,Na,K,Be,Mg,Ca)嵌入C70后的形成能进行了讨论,之后计算了[M@C70]各体系与[Mg(H_2O)6]2~+的相互作用能,并采用自然键轨道理论(NBO)研究了电荷转移的情况,最后进行了电子密度拓扑分析.结果表明,内嵌原子半径越大,[M@C70]各体系的热力学稳定性就越高,转移至[Mg(H_2O)6]2~+的电荷也随之增加,二者之间的相互作用属于闭壳层相互作用及共价作用.     

(GaAs)n(n=1-4)原子链电子输运性质的理论计算

本文利用密度泛函理论结合非平衡格林函数的方法, 对 (GaAs)_n(n=1-4)直线原子链与Au(100)-3×3两半无限电极耦合构成Au-(GaAs)_n-Au纳米结点的电子输运性质进行了第一性原理计算. 在各结点拉伸过程中, 对其结构进行了优化, 得到各结点稳定平衡结构时Ga-As的平均键长分别为0.220, 0.224, 0.223, 0.223 nm, 平衡电导分别为2.328G₀, 1.167G₀, 0.639G₀, 1.237G₀; 通过对结点投影态密度的计算, 发现电子传输主要是通过Ga, As原子中p_x与p_y电子轨道相互作用形成的π键进行的. 在0-2 V的电压范围内, 对于(GaAs)_n(n=1-3)的原子链的电流随电压增大而增大, I-V曲线呈线性关系, 表现出类似金属导电行为; 对于(GaAs)₄原子链在0.6-0.7 V, 0.8-0.9 V的电压范围内却存在负微分电阻现象.