QEq方法在ReaxFF中的应用与局限性 QEq方法简介 QEq(电荷平衡方法)是ReaxFF中最常用的电荷分布方法，其主要特点包括： 基本原理：在每个分子动力学步骤中，系统会寻找使静电能量最小化的电荷分布 约束条件：整个系统必须保持电中性（总电荷为零） 核心思想：电子可以在系统中流动，直到每个原子的"电化学势"（或电负性）相等 数学表达： min_q E_ele(q) = […]
lammps
14 posts
LAMMPS完成reaxff模拟后统计产物分子 输入脚本示例 # ..... units real atom_style charge read_data data.CHO pair_style […]
如何在LAMMPS运算中断后继续计算 使用restart文件继续计算 预防性设置 在运行模拟之前，建议在input文件中设置断点保存： restart 1000 restart_file.1 restart_file.2 这种设置的优点： 每1000步会轮流写入两个文件 […]
Roy Kid 中国科学技术大学 物理化学硕士在读 在做高分子体系中添加大直径填料的时候遇到了一个问题：相同的体系，加入填料后的模拟速度是不加入的五分之一。查看output的时候发现，neigh一项占去了总时间的一半以上。这就很诡异，因为正常情况应该是Pair占到七成左右比较合理。 想要解决问题，就需要从lammps计算原子对之间的势的算法说起。 首先，我要计算一个原子受到那些原子的作用，首先得去寻找周围的原子。这样，为了避免每个timestep都要去搜索，lammps先拉一个叫neighbor list 的表，这个表里储存了所有周遭原子的信息。这个表的范围是多大呢，(cutoff + […]
