计算化学公社

标题: 如何求出每个原子在各方向上的力 [打印本页]

作者
Author: 1154975925 时间: 2024-9-21 17:27
标题: 如何求出每个原子在各方向上的力
为了计算原子受到的电场力，模拟两次，一次带电场一次不带电场，设置nstfout通过gmxdump -f all.trr |grep "f. 977" > force.txt 导出每个原子的受力，将两次force.txt里的原子受力相减是不是就是每个原子受到的电场力。gmx dump可以导出原子每一帧的X Y Z方向上的力，两次相减就是每个原子在各方向上受到的电场力，这样子可以吗

作者
Author: sobereva 时间: 2024-9-21 23:16
根据F=Eq，直接根据电场矢量E和原子电荷q就知道电场力了，根本不用这么麻烦

作者
Author: 1154975925 时间: 2024-9-22 15:01

sobereva 发表于 2024-9-21 23:16
根据F=Eq，直接根据电场矢量E和原子电荷q就知道电场力了，根本不用这么麻烦

是的我搞错了，因为电场力是单方向的

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