Structure and Spectrum of Binary Classic Systems Confined in a Parabolic Trap

The static and dynamic properties of the two-dimensional classic system of two-species interacting charged particles in a parabolic trap are studied. The ground state energy and configuration for different kinds of binary systems are obtained by Monte Carlo simulation and Newton optimization. The spectrum and normal modes vectors can be gained by diagonalizing the dynamical matrix of the system. It is found that the total particle number, particle number and mass-to-charge ratio of each species are decisive factors for the system structure and spectrum. The three intrinsic normal modes of single species Coulomb clusters are inherent, concluded from our numerical simulations and analytical results.     

ELECTRONIC STRUCTURE AND CRYSTALLINE ELECTRIC FIELD IN Rco5(R= Y,La, Ce,Pr, Nd,Sm, Gd,and Tb)

Self consistent charge and spin polarized local spin-density approximation functional theory calculations based on the discrete variational method have been performed for RCo₅(R=Y, La, Ce, Pr, Nd, Sm, Gd, and Tb) compounds. The partial density of states of the Pr atom in the PrCo₆Co₁₂ cluster is established to be strikingly similar to that of the Ce atom in the CeCo₆Co₁₂ cluster, supporting the suggestion that the Pr atom is valence fluctuating. The radii ＜r_4f＞ and ＜r₄f²＞ of the 4f electrons of the R atom from La to Tb, except Ce, show the lanthanide contraction. The crystalline electric field (CEF) parameter A⁰₂ at the R site is calculated using a real charge distribution ρ(R) in the cluster, except for Pr and Nd, and is in agreement with that evalu ated based on the single-ion model. This result shows that the CEF parameter A⁰₂ is mainly determined by the near electronic structure. There exists a hybridization in a certain degree between the light rare-earth R-4f and Co-3d orbitals in some single-electron-molecular-orbitals, which are n ear the Fermi energy level and occupied by electrons. For light rare-earths the R-4f electrons in R Co₆Co₁₂(R=Y, La, Ce, Pr, Nd, and Sm) clu sters are not localized entirely and a small amount of the R-4f electrons have itinerant properties.     

氧化镁纳米多晶的微结构和磁性

实验发现,宏观晶体是非磁性的氧化镁时,其多晶样品有弱铁磁性.本文用第一性原理电子结构方法研究了氧化镁表面、纳米颗粒和晶界的磁性.计算结果表明:绝缘的氧化镁表面可以是导电的,并且有与之相关的铁磁性;磁性表面的共同特征是在表面上有氧原子富集,包括(111)表面的纯氧原子层,(114)表面的氧原子链;其他高晶面指数表面也会有氧原子富集区域;氧化镁纳米颗粒的磁性出现在高晶面指数表面以及不同晶面交界的棱及其顶角等有氧原子富集的区域,这种由氧原子富集而形成的磁性有巡游特征.氧化镁Σ7[111]和Σ5[001]晶界的计算结果表明:在没有氧原子富集的情况下,多晶样品中晶界的磁性很弱,而在有氧原子富集的情况下,晶界磁性比较强.因此可以推断多晶样品的磁性主要出现在多晶表面、晶粒包围孔隙、微裂纹界面、晶界和其他晶体缺陷等有氧原子富集的区域.这种残余磁性可以通过热处理等结构优化过程而削弱甚至消除.     

支撑材料基因工程的高通量材料集成计算平台

总结材料计算和数据管理面临的挑战,阐述高通量材料集成计算的内涵,介绍我们在国内率先研制出的高通量材料集成计算平台和软件框架MatCloud.平台支持与不同高性能计算集群的动态绑定,基于图形化界面的第一性原理计算流程设计、持久化、解析和执行;大批量计算作业的生成、提交、运行监控,及容错和纠错的自动流程实现;支持数据的提取、规范化加工、及自动存储等.MatCloud形成了一个高通量材料集成计算软件框架,支持个性化插件的快速研发.     

Cu2ZnSnS4晶界性质与光伏效应的第一性原理研究

本文应用第一性原理电子结构计算方法研究了锌黄锡矿Cu₂ZnSnS₄ (CZTS)晶界的性质: 包括微结构和电子结构及其对光伏效应的影响. 计算结果表明: 从p-n结区扩散过来空穴可以翻越一定势垒后被晶界俘获, 晶界进一步提供载流子扩散的快速通道, 使得这些空穴可以快速运动到阳极. 少数载流子电子在晶界中心区附近感受到很高的静电势垒, 但其高势垒两侧存在的势阱可以束缚少量电子. 对多数载流子空穴, 晶界中心则是势阱, 势阱两边有阻止空穴扩散到晶界中心的势垒. 由于CZTS晶体的易解理面是(112), 晶界面与(112)面平行的扭转晶界∑ 3^*[221]和∑ 6^*[221] 等不破坏原有晶体的基本结构, 它们的晶界能很小, 而且其电子结构与晶体内部基本相同, 因此尽管它们大量存在于CZTS材料中, 但是对材料性质仅有很小的影响. 通过比较晶体、晶界、空腔的表面和纳米棒的电子结构和光吸收系数, 我们可以看出: 这些微结构会在带隙内引入新的能级(复合中心), 同时高的孔隙率会降低(大于1.3 eV)光的吸收系数, 因此提高CZTS薄膜的致密度是提高CZTS太阳能电池效率关键.     

硅中硼离子注入的带电缺陷动力学模拟

为准确描述硼离子注入硅后缺陷/杂质的动力学物理过程,获得硼浓度空间分布及其演化行为,构建一个跨尺度带电缺陷动力学模型,考虑离子注入缺陷的产生及其演化的多种微观过程,包括缺陷电荷态和带电缺陷间的反应、硼-自间隙团簇（BICs）演化以及缺陷与载流子相互作用等物理过程。模拟得到与实验一致的硼浓度深度分布。结果表明：BICs对硼浓度的深度分布起主要作用,而间隙硼（B_I）导致硼浓度分布向深处扩展;计及缺陷的不同电荷态修正自间隙（I）和硼间隙（B_I）的扩散系数,从而更准确地描述硼浓度分布。模型揭示了硼离子注入硅发生的物理过程和微观机理,证明BICs和缺陷真实的电荷态是描述硼浓度分布的重要因素,为半导体器件制造与研发提供理论指导。     

三元铁基金属氮化物的高压复分解反应合成北大核心CSCD

以二元金属氧化物(氧化铁、氧化钴、氧化镍)和六方氮化硼为反应前驱体,在大腔体压机提供的高温高压条件下(5GPa、1 673K),通过发生复合高压固相复分解(HPSSM)反应合成组分可调控的三元铁基金属氮化物圆球状块体材料ε-Fe_(3-x)M_xN_(1+δ)(M=Co,Ni)。并利用X射线粉末衍射(XPD)、场发射扫描电子显微镜(FE-SEM)、高分辨率透射电子显微镜(HRTEM)等多种材料表征手段对高压合成的三元铁基金属氮化物进行结构表征,同时基于密度泛函理论(DFT)的第一性原理计算探究压力对HPSSM反应的影响。研究结果表明,高压密闭环境有利于制备高质量金属氮化物,HPSSM反应合成法是制备铁基金属氮化物块体材料的一种有效方法。     

ELECTRONIC STRUCTURE AND CRYSTALLINE ELECTRIC FIELD IN Rco5(R= Y, La, Ce, Pr, Nd, Sm, Gd, and Tb)

Self consistent charge and spin polarized local spin-density approximation functional theory calculations based on the discrete variational method have been performed for RCo₅(R=Y, La, Ce, Pr, Nd, Sm, Gd, and Tb) compounds. The partial density of states of the Pr atom in the PrCo₆Co₁₂ cluster is established to be strikingly similar to that of the Ce atom in the CeCo₆Co₁₂ cluster, supporting the suggestion that the Pr atom is valence fluctuating. The radii ＜r_4f＞ and ＜r₄f²＞ of the 4f electrons of the R atom from La to Tb, except Ce, show the lanthanide contraction. The crystalline electric field (CEF) parameter A⁰₂ at the R site is calculated using a real charge distribution ρ(R) in the cluster, except for Pr and Nd, and is in agreement with that evalu ated based on the single-ion model. This result shows that the CEF parameter A⁰₂ is mainly determined by the near electronic structure. There exists a hybridization in a certain degree between the light rare-earth R-4f and Co-3d orbitals in some single-electron-molecular-orbitals, which are n ear the Fermi energy level and occupied by electrons. For light rare-earths the R-4f electrons in R Co₆Co₁₂(R=Y, La, Ce, Pr, Nd, and Sm) clu sters are not localized entirely and a small amount of the R-4f electrons have itinerant properties.     

Anisotropic properties of TaS2

The anisotropic properties of 1T- and 2H-TaS2 axe investigated by the density functional theory within the framework of full-potential linearized augmented plane wave method. The band structures of 1T- and 2H-TaS2 exhibit anisotropic properties and the calculated electronic specific-heat coefficient γ of 2H-TaS2 accords well with the existing experimental value. The anisotropic frequency-dependent dielectric functions including the effect of the Drude term are analysed, where the ε^xx（ω） spectra corresponding to the electric field E perpendicular to the z axis show excellent agreement with the measured results except for the ε1^xx（ω） of 1T-TaS2 below the energy level of 2.6 eV which is due to the lack of the enough CDW information for reference in our calculation. Furthermore, based on the values of optical effective mass ratio P of 1T and 2H phases it is found that the anisotropy in 2H-TaS2 is stronger than that in 1T-TaS2.     

氮原子修饰点缺陷石墨烯结构和性质的理论研究

氮原子掺杂石墨烯对基于石墨烯的器件和催化研究具有重要的应用价值.本文采用基于密度泛函理论的计算方法,研究了氮原子修饰的C-Bridge(碳原子吸附在石墨烯碳碳键桥位)、C-Top(碳原子吸附在石墨烯一个碳原子上方)和C7557(碳原子对吸附在石墨烯碳环上方)三种不同点缺陷类型的石墨烯物理性质.讨论不同缺陷石墨烯结构在用氮原子进行修饰前后体系的稳定性、电子结构等;计算得到了缺陷处原子的分波态密度(PDOS)图,分析了原子间的相互作用;模拟出氮原子修饰后缺陷石墨烯恒流模式的STM图像,以便和实验上得出的图像进行对比.计算结果表明,对于所选取的三种不同缺陷,氮原子能够较稳定地吸附在缺陷表面.C-Bridge和C-Top缺陷结构本身具有磁矩,经氮原子修饰后结构磁矩消失.与之相反,C7557缺陷结构本身没有磁矩,经氮原子修饰后缺陷体系带有磁矩.另外,C-Bridge和CTop两种不同缺陷结构石墨烯经过氮原子修饰后,体系几何结构变得完全一样.