【中南大学仲奇凡教授课题组】改质煤沥青的分子结构建模及其原子尺度的焦化机理研究(Fuel 2023)

Posted · Add Comment

相关背景 煤沥青是煤焦油蒸馏工艺的副产品,由于其具有含碳量高、杂元素含量低、廉价且来源广泛等优势。对煤沥青进行改质可以获得更高的结焦值和软化点,改质煤沥青是石墨电极、碳纤维、碳纳米管、C/C复合材料等高附加值碳制品的重要优质原料。改质煤沥青的焦化是制备高附加值碳素产品的必经工序,它决定了产品的真密度、机械强度、电导率、石墨化度等性能。探索改质煤沥青焦化机理对优化高附加值碳材料的制备工艺具有重要意义。然而,现有科学文献对MCTP的微观结构组成及焦化过程的深度演化机制研究较少,本文精确构建了改质煤沥青的分子结构模型,并基于ReaxFF分子动力学方法对改质煤沥青的焦化机理在原子尺度上进行了深度研究,为实现改质煤沥青更高值化的利用,以及高附加值碳材料制备工艺的改进提供理论了依据。 研究亮点 通过一系列的检测分析(XRD, 13C NMR, FT-IR, XPS, MALDI-TOF-MS等)获得了MCTP的微观结构特征和组成。MCTP的碳骨架主要由芳香结构组成,芳香结构单元主要为苯和萘。脂肪族结构主要以甲基和亚甲基的形式存在。MCTP的分子量分布集中在1500 Da范围内。MCTP中含O、N和S的官能团分别以醚、季氮和无机硫的形式存在。在此基础上,建立了能够表征改质煤沥青微观结构特征和组成的分子模型(C93H59NO)。利用AMS2019软件中ADF模块计算改质煤沥青分子模型的量化性质,得到的13C NMR和FTIR图谱与实际检测结果吻合,验证了模型的合理性。 运用所建立的分子模型、ReaxFF MD模拟(AMS2019软件ReaxFF模块)、TG-MS、XRD和SAXS分析,揭示了改质煤沥青焦化过程中的挥发份去除及结焦机理。改质煤沥青焦化过程中挥发物的去除是由改质煤沥青分子边缘结构的破坏和活性自由基的产生引起的。主要挥发产物为H2、H2O、CO、CH4和C2H4。挥发分主要在430 ~ 900 K的温度范围内析出。对于焦化过程,XRD和SAXS分析表明,沥青焦样品的结晶和石墨化程度在挥发分去除阶段(低于973 K)被破坏,在高温阶段(>高于973 K)得到改善,Lc、Nc、La、Rg和孔隙率在1573 K时分别达到3.81 nm、12.03 nm、1.68 nm、18.54 Å和9.9%。在ReaxFF MD模拟过程中,沥青焦核的真密度(最终为2.2 g/cm3)、RDF临界峰强度、sp2杂化键比例(最终为41.4%)和六元环比例(最终为63.3%)均呈现先下降后上升的趋势。实际检测分析(XRD、SAXS)和ReaxFF MD模拟得到的结论可相互映证。沥青焦核的形成可分为两个阶段:焦核的形成阶段和有序化阶段。在形核阶段,改质煤沥青分子中的开环反应和脂肪链的交联导致了初级焦核的形成,主要通过三条反应路径,在路径一中,芳香烃通过芳香甲基直接连接;在路径二中,由边缘六元环断裂形成脂肪链,然后通过脂肪链结合实现芳烃的聚合;在路径三中,两个芳烃边缘的C原子相互结合形成五元或七元环结构。在有序化阶段,在初级焦核内会发生缩聚和芳构化反应,最终形成沥青焦微晶核。 图1 文章摘要图 图2 研究流程示意图 图3 改质煤沥青分子结构模型及验证 图4 改质煤沥青结焦机理图 参考文献 Zihan You, Jin Xiao, Gang Wang, Zhen Yao, Ye Wan, Qifan Zhong, Molecular representation and atomic-level coking evolution investigation of modified coal tar pitch […]

氟化氢刻蚀二氧化硅的ReaxFF分子动力学模拟

Posted · Add Comment

首尔大学化工学院化学与生物工程学院、韩国国民大学高级材料工程学院,最近开发了一个用于Si/O/H/F系统,以模拟HF刻蚀剂刻蚀SiO2的ReaxFF力场。其中,利用DFT计算得到的训练集,包括反应物/产物的结构、键离解能、价角畸变、SiO2团簇与SiO2板与HF气体的反应等,对ReaxFF参数进行了优化。使用ReaxFF计算的结构和能量与QM训练集很好地匹配。 作者利用该力场,对活性HF分子对SiO2衬底的刻蚀过程进行了分子动力学模拟。研究了不同入射能量下,HF刻蚀剂的刻蚀率和反应产物数。这些模拟表明,该力场对SiO2刻蚀过程中原子表面反应,能够很好的描述。AMS的Molecule Gun功能能够模拟刻蚀模型。 AMS中包含拟合、优化ReaxFF的功能,帮助用户自建新力场,但该功能要求用户对反应、数值拟合有一定的了解。 参考文献: Dong Hyun Kim, Seung Jae Kwak, Jae Hun Jeong, Suyoung Yoo, Sang Ki Nam, YongJoo Kim*, and Won Bo Lee*, Molecular Dynamics Simulation of Silicon Dioxide Etching by Hydrogen Fluoride Using the Reactive Force Field, ACS Omega 2021

分子动力学网络课程
“Atomistic Simulation of Thermal Transport Across Interfaces”

Posted · Add Comment

QuantumWise 将于近期举办“跨界面的热输运模拟”的网络课程。 网络课程时间:2017 年 2月 14 日 4-6 PM 或 2 月 15 日 3-5 AM(均为北京时间) 课程安排如下: 入门讲座(半小时):介绍热输运基本概念 操作演示(一小时):特别关注非平衡分子动力学模拟和基于非平衡格林函数的声子输运,举例介绍如何模拟透过硅晶界的热传导。 答疑(半小时) 本次网络课程以英文进行。 点此注册 成功注册后,您会收到一封确认邮件,提示您如何加入本次网络课程。 有疑问请联系我们:sales@fermitech.com.cn 或 info@quantumwise.com。 上期“使用 VNL 和 ATK-Classical 进行分子动力学模拟”的网络课程讲稿下载。 参加者准备: 没有任何计算模拟经验的人也可以从中获益,但建议参加者应该有基本的原子级别建模和模拟的经验。 参加者可以提前在自己的电脑上安装 VNL-ATK,以便在操作演示中进行练习。所有用户都可以申请 30 天免费试用 license(点此申请),学术机构的研究者可以申请 ATK-Classical 的永久免费 license(点此申请)。 深入了解VNL-ATK中的经验力场计算引擎ATK-Classical,请参考: 简介:http://www.fermitech.com.cn/vnl-atk/atk-classical/ 文章: http://arxiv.org/abs/1701.02495 与ATK-Classical相关的中文实例教程有: 分子动力学基础 模拟气相沉积薄膜生长过程 模拟离子轰击单层石墨烯 缺陷碳纳米管的杨氏模量 界面热导的模拟 VNL-ATK是先进的材料与器件模拟平台,除经验力场之外,还支持DFT、DFTB、半经验量子力学计算,可以模拟材料和电子器件的各种性质(了解 VNL-ATK 的更多功能)。

分子动力学网络课程
“How Atoms Move: Introduction to Molecular Dynamics Simulations with VNL and ATK-Classical”

Posted · Add Comment

QuantumWise 将于近期举办“使用 VNL 和 ATK-Classical 进行分子动力学模拟”的网络课程。 网络课程时间:2017 年 1 月 24 日 4-6 PM 或 1 月 26 日 3-5 AM(均为北京时间) 课程安排如下: 入门讲座(半小时):介绍分子动力学的基本原理、方法和应用 操作演示(一小时) 答疑(半小时) 本次网络课程以英文进行。 点此注册 成功注册后,您会收到一封确认邮件,提示您如何加入本次网络课程。 有疑问请联系我们:sales@fermitech.com.cn 或 info@quantumwise.com。 参加者准备: 没有任何计算模拟经验的人也可以从中获益,但建议参加者应该有基本的原子级别建模和模拟的经验。 参加者可以提前在自己的电脑上安装 VNL-ATK,以便在操作演示中进行练习。所有用户都可以申请 30 天免费试用 license(点此申请),学术机构的研究者可以申请 ATK-Classical 的永久免费 license(点此申请)。 深入了解VNL-ATK中的经验力场计算引擎ATK-Classical,请参考: 简介:http://www.fermitech.com.cn/vnl-atk/atk-classical/ 文章: http://arxiv.org/abs/1701.02495 与ATK-Classical相关的中文实例教程有: 分子动力学基础 模拟气相沉积薄膜生长过程 模拟离子轰击单层石墨烯 缺陷碳纳米管的杨氏模量 界面热导的模拟 VNL-ATK是先进的材料与器件模拟平台,除经验力场之外,还支持DFT、DFTB、半经验量子力学计算,可以模拟材料和电子器件的各种性质(了解 VNL-ATK 的更多功能)。

VNL作为LAMMPS的用户界面

Posted · Add Comment

概述 LAMMPS 是一个应用广泛的经验力场分子动力学程序,使用VNL可以: 快速构建多种模型并导出LAMMPS的文件; 导入LAMMPS 分子动力学轨迹,并使用VNL中的分析功能进行分析。 现在申请【获取永久免费学术版的VNL图形界面】,即可获得这些功能! 使用VNL快速构建结构模型 导入数据库中的结构 导出LAMMPS的数据格式 构建块体、表面、界面结构 构建分子模型 构建随机合金、多晶体系 构建无定形体系 使用packmol将分子充填到晶胞 更多详细介绍,请参见:【VNL的建模功能】。 使用 VNL 分析 LAMMPS 的结果 VNL支持超大体系的显示和分析 原子数可达百万级 VNL可以直接由分子动力学轨迹得到体系的各种分布函数和性质: 径向分布函数 速度/动能分布 速度自相关函数 局域质量密度分布 局域应力分析 从MD轨迹计算声子DOS 配位数分析 角分布函数 均方位移 Evolution of average nearest neighbor number with tine 中子散射结构因子 弹性常数(应力-应变曲线) 多种函数重叠作图 支持柱状图、线状图 显示MD模拟过程的原子速度 显示MD动画并支持自定义显示、按帧导出结构 导出动画文件 相关实例教程 使用VNL产生无定形结构 英文教程:Generating Amorphous Structures VNL和LAMMPS一起使用进行分子动力学研究 教程:在VNL中导入LAMMPS轨迹文件。 […]

QuantumATK:更强大、更灵活的材料动力学模拟

Posted · Add Comment

概述 动力学模拟是一种重要的原子级模拟方法,通过求解原子运动的经典力学牛顿方程对相空间进行采样,不仅可以研究体系在相空间的演化过程,还可以通过产生的系列结构(系综)通过统计方法得到体系在非零温度下的各种性质。 动力学过程中的原子间相互作用力则可以通过多种方法求得,可以是密度泛函理论,也可以是经验力场。 使用QuantumATK进行材料动力学模拟 可以用多种能量-力计算方法密度泛函理论(DFT-LCAO 和 DFT-PlaneWave):支持几千原子级别超大体系的计算半经验量子力学模型(SemiEmpirical):支持密度泛函紧束缚近似(DFTB)经验力场(Forcefield):支持几百种经验势参数(含ReaxFF)机器学习力场:高效率平衡DFT的精度和力场的速度支持多种系综和理论方法NVE velocity verletNPT with stress maskNPT/NVT(Berendsen)NPT MelchionnaNPT with stress maskLangevin多种方法初始化速度灵活的结构限制固定部分原子分别固定 x、y、z 坐标MD 过程中固定质心限定布拉维格子(可同时设定目标应力)刚性限定原子相对位置关系支持多种动力学方法平衡态分子动力学非平衡态分子动力学(RNEMD)计算热导time-stamped force-bias Monte Carlo长时域的动力学方法Metadyamics(PLUMED):更快的对能量(自由能)面进行采样,获得大范围的结构-能量信息自适应动力学蒙特卡罗方法(adaptive kinetic Monte Carlo):研究结构变化与机理可控制局域温度、设定升温速率所有恒温器、恒压器支持线性升温或降温计算过程中分析Python 脚本支持计算过程中分析或添加其他限制条件部分电荷分析可视化原子速度轨迹或单个结构分析工具(同时支持导入 VASP、LAMMPS 结果分析)径向分布、角分布函数速度自相关函数局域质量密度分布配位数分析均方位移最近邻数中子散射结构因子速度/动能分布局域结构(Voronoi)中心对称性从轨迹计算声子态密度使用脚本可以方便的对超大体系的部分原子进行以上各种分析力学性质力和应力(解析 Hellmann-Feynman)弹性常数(应力-应变曲线)局域应力 QuantumATK:高效的动力学引擎 基于QuantumATK高效的DFT、SE和ForceField计算引擎和MD代码的优化,分子动力学计算速度有明显优势。其中DFT、SE、ForceField均支持MPI大规模并行,并获得极大的速度提升。 基于DFT的MD计算速度测试结果 7056原子的分子动力学(水分子):在64MPI并行时,第一步MD耗时4分钟。 3220原子的分子动力学(硅固体):144MPI进程并行,1步耗时约10分钟。 基于ForceField的MD并行加速测试(1百万SiO2原子) DFT+ForceField混合方法动力学 QuantumATK提供更灵活的分子动力学模拟框架,可以在一个动力学模拟过程中混合使用DFT和ForceField。 图:LiFePO4中Li在外电场下的扩散动力学模拟。考虑电场的同时,DFT的引入可以反映原子电荷在动力学过程中的涨落。 time-stamped force-bias Monte Carlo方法 使用time-stamped force-bias Monte Carlo方法(基于DFT、DFTB、力场)代替普通的分子动力学,可以研究更长时间的平衡、沉积、无定型化、扩散、快速熔化–退火、对罕见事件采样。模拟过程可以恒定温度,也可以是线性升温降温。 与MD类似,但是原子位移使用Monte Carlo方法采样将结构在给定温度下平衡并采样(图:HfO2的fbMC与普通MD的结果比较) 模拟时长可以是传统MD的100倍(图:SiO2 的 800 ps fbMC动力学轨迹分析) 与普通MD混合使用(图:fbMC恒温—fbMC降温—NPT MD恒温) […]

 
  • 标签

  •