ReaxFF教程：使用ReaxFF中巨正则蒙特卡洛（GCMC）方法模拟锂化放电过程的电压特性

Posted · Add Comment

英文教程 点击此处下载英文教程 英文教程阅读说明 准备知识 该功能版本号≥2017.211，下载地址：https://www.scm.com/support/downloads/development-snapshots/，svnRevision大于r62841的都可以 熟悉ReaxFF建模、保存任务、提交任务等基本操作 熟悉巨正则系综的含义更好 知道如何在命令行里面运行命令（本教程需要在命令行运行一次分析命令，得到电压特性数据） 知道如何将数据从文本格式转移到Excel中，并生成曲线，参考：如何从图谱中导出数据到Excel表格 阅读说明 该教程非常详细，实际上初学者也可以照着说明一步一步做出来，当然结果会有差异，这是正常的。本教程主要包含如下几个小节： 从外部导入Alpha硫晶体的*.cif文件，分别采用NPT系综低温弛豫（ADFinput > ReaxFF > Main > Task: Molecular Dynamics）、能量最小化（ADFinput > ReaxFF > Main > Task: Geometry Optimization）来得到其平衡结构。需要说明的是，教程中能量最小化得到的能量是-8535.99 kcal/mol，实际我们得到的能量可能更低。 计算Li原子的化学势：计算体心立方锂晶体的8*8*8超胞（共512个Li原子）的能量，除以512，即得到Li原子的化学势-37.70kcal/mol 对Alpha硫晶体进行巨正则系综蒙特卡洛模拟 从上一步GCMC模拟生成的*.rxkf文件中，抽取电压特性等数据：将LiVoltageProfile.py文件和*.rxkf放到同一个文件夹内，在ADF环境变量对当前命令行生效的前提下，在命令行执行：startpython LiVoltageProfile.py gcmc_test.rxkf。将生成一个文件，名为voltage_profile.out。该文件实际上是一个文本文件，给出了4列数据，其中第二列是LixS中的x值，第三列是电压，第四列是体积变化。用户可以自行用这两列数据生成教程中的曲线 上面这几步的计算模拟文件（请点击此处下载），在阅读困难的情况下，可以辅助您阅读教程，但由于GCMC这一步很耗时间（对应文件gcmc_test.adf），可能需要一天左右才能完成，因此设置的参数数值放宽了条件，没有教程里面那么严格，因此这一步仅供参考。

Read More

ADF中ReaxFF的新功能：局域温度分布、原子电荷分布的可视化

Posted · Add Comment

ADF最新开发版r63428(最新开发版下载地址：https://www.scm.com/support/downloads/development-snapshots/)，能够计算原子的温度以及在空间上的均值，并自动保存在轨迹文件中。该功能通过在*.run文件中control区域添加如下两行内容开启： 1 localt Request local temperature averaging between atoms 3.0 localr Set the awarness radius to 3.0 Å (Default 5.0) 其中1 localt表示开启该功能；3.0 localr表示设置局域范围半径设置为3.0埃，默认值为5.0埃。 通过ADFMovie窗口中，为每个原子着色，来表征原子的温度分布，不同颜色表示不同温度： 除了原子温度以外，如上图所示，还增加了加速度、压力、速度、原子电荷分布等数据的图形化显示。 例如下面的视频演示了温度梯度、热的波动情况： 该模型是用ReaxFF的分子枪功能，入射一个石墨烯球形分子到石墨烯片层上，导致的温度传递。分子枪功能参考：水分子在石墨烯表面形成水簇的模拟（分子入射、沉积）或者：https://www.scm.com/doc/Tutorials/ReaxFF/Bouncing_a_buckyball.html

Read More

ADF Highlight：神奇的钍-砷多键（ Nature Comm., 2017）

Posted · Add Comment

参考文献： E. P. Wildman, G. Balázs, A. J. Wooles, M. Scheer, and S. T. Liddle, Triamidoamine thorium-arsenic complexes with parent arsenide, arsinidiide and arsenido structural motifs, Nature Comm. 8, 14769 (2017). 2017年Nature Communication报道了钍与体积较为庞大的triamidoamine配体形成钍-砷多重键，从而得到稳定的结构。文章使用ADF的相对论密度泛函理论，研究了ThAs(H)K、ThAs(H)Th、ThAsTh)的键合作用，并展示了从5f、6d、7s轨道生成的极化多重共价键。 这项研究进一步提高了对锕系配合物多重键的基础认知，最终有助于核废料处理过程中的配体设计。另外可以参考：novel uranium complexes。 以及其它相关文献： Thorium–phosphorus triamidoamine complexes containing Th–P single- and multiple-bond interactions, Nat. Comm. 7, 12884 (2016). Isolation of Elusive HAsAsH in a Crystalline Diuranium(IV) […]

Read More

ADF Highlight：次锕系元素中的共价键合（JACS, 2017）

Posted · Add Comment

参考文献： M. P. Kelley, J. Su, M. Urban, M. Luckey, E, R. Batista, P. Yang, and J. C. Shafer, On the Origin of Covalent Bonding in Heavy Actinides, J. Am. Chem. Soc., 2017, 139, 9901–9908 (2017). 在核燃料循环中，将高放射性次锕系元素从镧系和其他裂变产物中分离出来，对于使用过的核燃料回收再利用，以及核废物储存，都是至关重要的。通常通过锕系元素与配体之间的相互作用差异来进行分离。上世纪四十年代，学界已经热烈争论过锕系-配体之间的“共价”概念，假定次锕系元素形成了部分共价键，导致其5f轨道与配体轨道混合。 著名的洛斯阿拉莫斯国家实验室使用ADF中的相对论密度泛函理论，研究了Am、Cm、Bk、Cf离子与DPA（dipicolinate）配体之间的成键。通过实验热力学测量，观察到Cm和Bk之间络合能的差异；并发现这个差异与锕系5f轨道能量降低到配体的配位氧原子2p轨道能级附近，以及5f轨道的离域性有关。分离过程常常有赖于络合能的这种微小差异，而能量的（近）简并导致共价的机制，可能为锕系分离方法设计带来新的策略。  

Read More

ReaxFF教程：ReaxFF计算锂离子扩散系数

Posted · Add Comment

最新的ADF2017在ReaxFF中增加了“通过轨迹直接计算”、“通过高温外推”两种方法，计算（致密）材料中离子扩散系数。 体系范例、流程参考文献： ReaxFF molecular dynamics simulations on lithiated sulfur cathode materials, Md Mahbubul Islam, Alireza Ostadhossein,  Oleg Borodin,  A. Todd Yeates,  William W. Tipton,  Richard G. Hennig,  Nitin Kumar  and  Adri C. T. van Duin* ,Phys. Chem. Chem. Phys., 2015,17, 3383-3393 本教程以LixS为例，计算锂离子的扩散系数。 技能要求 熟悉ReaxFF的基本操作（参考费米维基ADF知识库） 一点点命令行操作常识（需要在命令行执行2个命令处理结果，从而得到扩散系数。其中Windows系统参考：Windows如何进行ADF的命令行操作；Linux、Mac可以直接打开命令行，但前提是ADF的环境变量在当前命令行是生效的，否则找不到所使用的这些ADF内部命令） 将文本数据导入Excel表格（参考：如何从图谱中导出数据到Excel表格） 英文教程下载 请点击此处。 教程说明 第一节：从网上下载cif文件，导入ADFinput中，并Edit → Crystal → Map atoms to (0..1)适调晶格坐标，并使用NPT系综，进行低温“弛豫”的方法，优化结构 第二节：生成两种不同x值的LixS无定形结构 第三节：通过NVT系综进行退火：升温、迅速降温；之后用NPT系综使体系达到平衡（因为退火时，结构可能发生收缩，但由于NVT系综的晶格常数是不变的，因此这样的结构并不稳定，需要用NPT系综使其达到稳态）。教程中温度上升、下降是退火过程的结果；密度增大，是NPT模拟过程的结果 第四节：这里开始，正式计算扩散系数：基于前面得到的结果，再次进行NVT计算：1600K、110000 steps、10000 non-reactive steps。注意需要加一个关键词：1 itrout，也就是itrout这个变量赋值为1。MD完成之后，进入命令行，分别执行两个命令，从而最终在生成的*.log文件中得到扩散系数 […]

Read More

QuantumATK独有的新功能：非平衡态格林函数方法研究半无限表面模型

Posted · Add Comment

概述 QuantumATK从最新版本开始引入了一个全新的、独一无二的半无限表面模型。与非平衡态格林函数方法配合，半无限表面模型可以比表面slab模型更完善的模拟表面体系。描述此方法的文章已经在Phys. Rev. B 96,195309 上发表（预印文章参见：https://arxiv.org/abs/1707.02141）。 模型 片层（Slab）模型 与其他的周期性模型程序类似，QuantumATK也可以用传统的Slab模型来描述表面体系，但Slab模型有很大的缺陷和局限： Slab最大的不足是无法模拟实际表面下方通常是无限大的块体材料; 由于厚度有限，Slab中的电子容易体现出量子限制的效应； 两个表面之间可能相互影响； 很难正确的在表面方向模拟外加电场； 经常需要表面钝化、偶极校正等额外补救措施。 单电极表面（One-probe surface）模型（或半无限表面模型） 为此，QuantumATK 基于 DFT 和格林函数方法方法开发了真正可以模拟半无限表面体系的模型，即将一个表面 Slab 模型耦合于半无限的块体结构上（见下图）。 这种模型有以下几个独特的优势： 算法复杂度降低，特别适合大体系计算； 表面性质对表面层数的依赖显著降低； 只需很少的层数就可以再现块体的电子态； 可以正确的施加垂直表面方向的电场，模拟电场对表面体系的影响。 应用 文章报道了半无限表面模型的原理和应用，这些应用实例展示了半无限表面模型和格林函数方法的精确性，也证实了这个模型在表面体系研究中比传统模型具有明显的优势。 计算过渡金属的功函数 计算贵金属和拓扑绝缘体的表面态 Ge(001)|Si薄膜半导体异质结构的能带对齐 电场对碘在Pt(111)上吸附的影响   NanoLab高级图形用户界面：专注于研究，更快获得结果 NanoLab 图形用户界面丰富易用的功能可以让用户专注于研究项目的科学问题，专心思考科学问题，更快的发现新材料、创建新结构，避免在数据的导入、导出、处理、作图等琐碎的问题上浪费时间。NanoLab 可以： 最强大的材料表面结构建模工具 直观的选择表面方向和表面超胞 最合理的表面结构优化方法 快速构建各种结构模型 内嵌晶体结构数据库 搜索在线晶体结构数据 应用领域 与QuantumATK中的结构优化、能量计算、CI-NEB过渡态搜索等功能结合，这种模型还在表面化学、催化等领域有广泛的潜在应用。 表面功函数 表面态 表面反应过渡态与催化 拓扑绝缘体 相关的实例教程 表面结构与吸附 表面分子吸附体系建模：中文教程、英文教程 硅表面重构研究：中文教程、英文教程 CO在Pd(100)表面的吸附：英文教程 […]

Read More

VNL-ATK 2017新版发布

Posted · Add Comment

QuantumWise已于7月1日正式发布了VNL-ATK 2017版，新版引入了多种新的计算方法和大量的性能优化。   更新概要 大幅改进电声耦合计算性能 新的“special thermal displacement”近似方法，高效描述声子辅助隧穿 新增两种计算精确带隙的方法 改进默认参数设置 创新性的Wigner-Seitz方法加速大超胞的动态矩阵计算 输出文件改为HDF5格式 新增多种酷炫的2D、3D作图选项 变晶胞NEB过渡态搜索（模拟相变） 连接Materials Projects和COD数据库（或自建数据库） 完整的投影能带和DOS License信息 已经获得2017之后版本授权的用户，直接下载安装即可使用； 最后版本为2016与之前的版本的用户，欢迎联系我们（sales@fermitech.com.cn）付费升级事宜； 欢迎各位新老用户联系我们（sales@fermitech.com.cn）获取最新版的试用许可。 本次更新详情 电声耦合计算 显著降低内存消耗（稀疏矩阵存储） 典型测试显示内存需求从800GB降为1.3GB 计算迁移率的新方法：常数弛豫时间方法（与BoltzTraP代码类似） 弛豫时间可以由实验确定，也可以从动态矩阵计算并在能量/k点取平均 对于电极材料重复形成的器件（长的纳米线、纳米管、二维片层）体系，非弹性电流计算速度提高100000倍，因为仅计算电极的动态矩阵和哈密顿量导数 将所有声子模式在能量范围内分组求和（近似），提速10-100倍 突破性进展 ATK2017引入了一种全新的特殊热位移（special thermal displacement；STD）近似，来考虑声子散射对IV曲线的影响。 STD-Landauer方法详见（PRB文章已发表），这种方法将全电子-声子耦合计算减少为中间区域的动态矩阵（LOE/XLOE计算所必须的）和每个偏压（温度）下的一个器件计算。基本思想是根据全部声子的正则系综平均生成单个的原子位移组合，即可包含所有的温度效应。 文章演示了这种方法是如何高效的研究声子对硅p-n结电流、硅块体和纳米线的迁移率的影响随温度变化的效应（图）。 总体性能改进 费米面以上能带数现自动设置 提速约2倍；过去版本可以手动设置，但多数用户都未注意到这个选项 新的自能存储选项 NoStorage（大大减少大体系的内存占用，但是每次重新计算） SaveInMemory（快速，但耗内存） StoreOnDisk（快速，不消耗内存，但消耗硬盘空间） 中间区域不同，但是电极相同时，重复使用自能计算结果 SparseGreenFunction改进 使用可以减少大洁面体系的内存消耗 尤其对于有限偏压时，每个能量（k）点使用多进程可以大大减少每个MPI进程的内存 No performance overhead anymore for the distributed Pulay mixer […]

Read More

ADF Highlight：六氯化镅中的5f共价电子（JACS, 2017）

Posted · Add Comment

参考文献： J. N. Cross, J. Su, E, R. Batista, S. K. Cary, W. J. Evans, S. A. Kozimor, V. Mocko, B. L. Scott, B. W. Stein, C. J. Windorff, P. Yang, Covalency in Americium(III) Hexachloride,  J. Am. Chem. Soc.,2017,139, 8667–8677 (2017). 本文结合Spin-Orbit Coupling DFT与跃迁偶极矩方法，计算分析了六氯化镅中氯K边XAS，观察到了Am-Cl键中5f与3d轨道混合的直接证据。 结合实验上氯K边XAS的数据，对比AmCl63- 和等电子的EuCl63-，确认了Seaborg在1954年提出的假设：AmIII 的5f轨道共价比EuIII的4f轨道混合更牢固。 ADF软件中，激发态计算orbital-difference DFT的设置方式：ADF input > Properties > Excitations(UV/Vis), CD > Method > […]

Read More

ReaxFF教程：优化ReaxFF力场

Posted · Add Comment

费米维基发布了综合使用ADF模块和ReaxFF模块，对ReaxFF力场进行优化的中文基础教程。详细教程请见：http://www.fermitech.com.cn/wiki/doku.php?id=adf:optimizeff。ADF软件集成了使用蒙特卡洛方法优化力场参数的脚本，目前已经可以较为方便的使用。  

Read More

【QuantumATK亮点文章】非对称的金属-半导体-金属MoS2异相结

Posted · Add Comment

科学家使用QuantumATK研究了非对称的金属-半导体-金属MoS2异相结。对不同掺杂浓度的透射和β-和β*-界面处的能垒研究表明，异相结的非对称能垒高度与掺杂浓度无关。这些结果显示了MoS2这个热门材料在电路元件设计中的一些有趣的性质。 使用QuantumATK对金属-半导体、半导体金属界面的结构进行优化得到的结果与实验一致。密度泛函理论和非平衡态格林函数方法则用来研究电子态和电子输运性质。 VNL中提供了非常方便的界面建模工具，用来构造并优化了金属-半导体、半导体-金属界面器件。两个界面的有效电子能垒有0.2V的差别，与掺杂浓度无关。 VNL中集成的Projected Local Density of States工具可以方便快捷的研究器件的输运性质和能带对齐，并直接作图，免去研究者编程和数据处理的痛苦。相关的教程参见：这里。 参考文献 [1] D. Saha, S. Mahapatra, “Asymmetric Junctions in Metallic-Semiconducting-Metallic Heterophase MoS₂”, IEEE Transactions on Electron Devices, vol.PP, no.99, pp.1 (2017)

Read More
 
  • 原位聚合限域实现高效蓝光钙钛矿 LED(Nature 2026)本文使用元宝 Hy3 preview 模型总结归纳 。 摘要 金属卤化物钙钛矿因其优异的光致发光特性，已成为颇具前景的发光二极管半导体材料。然而，其性能仍受限于在衬底上原位合成钙钛矿纳米晶过程中“高结晶度”与“小尺寸”之间的固有矛盾。本文报道了通过原位聚合驱动的纳米晶限域策略，合成了由高质量纳米晶组成的钙钛矿薄膜，从而实现了高效的蓝色钙钛矿发光二极管。原位形成的聚合物网络在钙钛矿纳米晶生长过程中施加了纳米尺度的空间限域作用，使纳米晶兼具小尺寸和 83% 的高光致发光量子产率。此外，具有足够配位位点的可聚合单体延长了钙钛矿团簇的晶格重排时间，从而提升了纳米晶的结晶度。所合成的钙钛矿纳米晶用于制备发光二极管后，在 491 [...]
  • 面向亚 5 nm MOSFET 的 Si 兼容高性能二维纳电子器件研究背景 随着集成电路器件尺寸持续向亚 5 nm 尺度推进，传统硅基 MOSFET 在短沟道效应、量子隧穿及功耗控制等方面面临严重瓶颈，难以满足未来高性能与低功耗芯片的发展需求。如何在极限尺度下实现高驱动电流、低泄漏及良好的栅控能力，并能够与现有硅基工艺相兼容，有效利用现有器件制造技术，成为下一代晶体管研究和应用亟须解决的关键问题。 研究内容 该研究针对超短沟道 MOSFET 性能提升这一核心问题，基于一类具有硅基工艺高度兼容性的二维单层材料 MNH2（M N = C Si Ge），在保证高载流子输运能力的同时实现与现有工艺体系的良好兼容，为器件实际应用提供基础。该工作系统地开展了器件物理与量子输运研究，并构建适用于亚 5 nm 尺度的双栅 MOSFET 模型，深入分析其输运特性与性能极限。 [...]
  • 基于统计形状模型的成人二尖瓣形态变异概述 二尖瓣反流（MR）影响着全球超过 2400 万人，是发病和死亡的主要原因之一。若不及时治疗，可能会进展为心房颤动、肺动脉高压，最终导致心力衰竭。经导管二尖瓣置换术（TMVR）因二尖瓣解剖结构的复杂性和高度变异性而面临独特挑战。为确保手术成功，TMVR 器械必须实现安全锚定、有效封堵和左心室流出通畅。然而，目前的临床前设计范式因对人群水平的二尖瓣解剖变异性认识不足而受到阻碍。 统计形状模型（SSM）为量化解剖结构变异性和构建具有代表性的 3D 解剖虚拟患者队列提供了一个稳健的框架，可作为生物力学计算机模拟临床试验的基础。本研究利用 72 例成人的前瞻性心电门控对比增强冠状动脉 [...]
  • 铜/银立方八面体团簇对卤素阴离子的封装机制及光电性质调控研究背景 在纳米材料科学的前沿领域，原子级精确的贵金属分子团簇因其卓越的稳定性和明确的结构，被视为构建新型功能材料的理想基石。通过引入阴离子作为模板剂，引导金属纳米团簇形成独特的笼状空间构型，并以此为平台深入探究“结构-性能关系”。然而，相比于研究较深入的金纳米团簇，铜与银团簇，尤其在阴离子封装特性上的对比研究仍显不足。近日，知名学者 Alvaro Muñoz-Castro 在无机化学领域权威期刊 Inorganic Chemistry 上发表了题为 “Copper and [...]

  • 标签

    3D打印 ABAQUS ADF AMS ATK Workshop BAND COSMO-RS DES DFT EDA EDA-NOCV NanoLab NMR NOCV OLED QATK Highlight QuantumATK reaxff Simpleware SOC SOCME TADF TDDFT ZORA 二维材料 催化 分子动力学 半导体 反应分子动力学 太阳能电池 密度泛函 工业 微电子 数字岩石 深共晶溶剂 溶解度 热解 燃烧 电声耦合 能量分解 自旋轨道耦合 迁移率 钙钛矿 镧系元素 骨科

  • 关于费米科技

    费米科技以促进工业级模拟与仿真的应用为宗旨，致力于推广基于原子级别模拟技术和基于图像模型的仿真技术，为学术和工业研究机构提供研发咨询、软件部署、技术攻关等全方位的服务。费米科技提供的模拟方案具有面向应用、模型新颖、功能丰富、计算高效、简单易用的特点，已经服务于众多的学术和工业用户。

    欢迎加入我们！（点击链接）

  • 最近更新

  • 联系方式

    • 留言板：点击留言
    • 邮箱：sales_at_fermitech.com.cn
    • 电话：010-80393990
    • QQ： 1732167264

  • 订阅费米科技新闻

    • 邮件订阅：
      您可以使用常用的邮件地址接收费米科技定期发送的产品更新和新闻。
      点击这里马上订阅！
    • 微信订阅：
      微信扫描右侧二维码。
    •