计算化学公社

标题: lammps做电场模拟时所加电场不生效（求助fix efield命令） [打印本页]

作者
Author: 事已至此先睡觉 时间: 2023-11-8 19:56
标题: lammps做电场模拟时所加电场不生效（求助fix efield命令）
本帖最后由事已至此先睡觉于 2023-11-9 14:26 编辑

在lammps里按手册格式使用fix efield命令对结构沿Y方向施加电场，模拟结果一点变化都没有，每一个group的温度都和不加电场的一样。
而且 fix 4 all efield 0.0 1.0 0.0 fix 4 all efield 0.0 10.0 0.0 variable E equal 10 / fix 4 all efield 0.0 v_E 0.0 的结果也都一样。
想求助各位大佬这种情况是命令设置不正确导致电场没生效，还是所用势函数没有涉及库仑力部分（我用的是Zhou合金EAM/alloy势）

2023/11/08
-------------------------------------------------------------------------------------------------------
使用set type ** charge ** 命令为每种原子施加电荷后，电场生效了，结果不同了
应该是用的势函数和data文件里没有电荷信息所以导致“施加的电场未生效”的结果

2023/11/09

作者
Author: Kin_Sing 时间: 2023-11-24 12:30
你好，加电场的问题解决了吗？因为我也遇到同样的问题，想请教一下您

作者
Author: 事已至此先睡觉 时间: 2023-11-24 16:33

Kin_Sing 发表于 2023-11-24 12:30
你好，加电场的问题解决了吗？因为我也遇到同样的问题，想请教一下您

可以的，你说说看

作者
Author: Kin_Sing 时间: 2023-11-24 17:18

事已至此先睡觉发表于 2023-11-24 16:33
可以的，你说说看

我在lammps的力场动力学模拟中，单独添加fix edield的交变电场，执行命令出现这样的问题：ERROR: Must not have electric field component in direction of periodic boundary when using charge equilibration with ReaxFF. (src/REAXFF/fix_qeq_reaxff.cpp:414) ，我判断是能量最小化的问题，然后我删掉能量最小化的命令去执行依然是这样的错误报告。然后我又删掉了随时间变化的电场和电场的能量（fix 20 all efield v_E 0 0
fix_modify 20 energy yes）去执行，结果发现命令是可以执行的，但是我删掉了随时间变化的电场的话去执行就相当于没有加家变电场吧，我是在X方向添加的交变电场。
方便加一下您微信吗？我的微信号是壹捌柒零零叁拾壹捌壹零

作者
Author: 事已至此先睡觉 时间: 2023-11-30 16:28

Kin_Sing 发表于 2023-11-24 17:18
我在lammps的力场动力学模拟中，单独添加fix edield的交变电场，执行命令出现这样的问题：ERROR: Must n ...

不好意思，前几天出差了，才看到，你的微信号有些奇怪我没加上，根据你的报错看，你的X方向是周期性边界，无法加电场，在别的帖子上也有人提到reaxff周期性边界无法加电场

作者
Author: Kin_Sing 时间: 2023-12-2 17:45

事已至此先睡觉发表于 2023-11-30 16:28
不好意思，前几天出差了，才看到，你的微信号有些奇怪我没加上，根据你的报错看，你的X方向是周期性边界 ...

好的，谢谢您。打扰您了，18700361810我的

作者
Author: nyt 时间: 2024-5-25 12:36
您好，使用fix efield命令前，自己对于原子set type **charge **这样不准吗？我尝试了您说的用ms中forcite模块施加电荷，但是力场信息缺少元素，请问您还有别的方法吗？

作者
Author: 事已至此先睡觉 时间: 2024-9-25 09:38

nyt 发表于 2024-5-25 12:36
您好，使用fix efield命令前，自己对于原子set type **charge **这样不准吗？我尝试了您说的用ms中forcite ...

我用了fix qeq 算的，自己加电荷了

作者
Author: chenym 时间: 2025-1-10 18:45
你好，我用set charge 设置电荷之后，因为挖了空位所以有电势差，想通过外加电场来平衡，因此使用fix 8 all efield -1.0e7 0.0 0.0，但是添加后仍显示 ERROR on proc 0: Non-numeric atom coords - simulation unstable (src/domain.cpp:551)其中我挖了十个+2价的Fe和23个+3价的Fe，不知道问题具体出现在哪里，还请指教，多谢！

作者
Author: Wolf123 时间: 2025-11-21 10:36
我在周期性体系中加电场时使用了fix ID group-ID qeq/rel/reaxff Nevery cutlo cuthi tolerance params gfile args，官网给出的gfile来自THEORY AND APPLICATIONS OF FLUCTUATING–CHARGE MODELS，我目前直接将其中的参数用到了计算中

欢迎光临计算化学公社 (http://bbs.keinsci.com/)