lammps计算偶极矩遇到的相关问题

yjw1204338757

各位老师好，最近想用lammps软件计算一块板（总共十三个层原子层）的偶极矩，我查手册发现有一个compute dipole/chunk 命令，按照手册里面说的命令运行之后，发现是计算了每个原子的偶极矩（图一）。compute cc1 ban chunk/atom molecule
compute myChunk ban dipole/chunk cc1
fix 8 all ave/time 1000 1 1000 c_myChunk[*] file tmp.out mode vector

在网上又搜到一个命令是计算所有原子的偶极矩。对比发现在z方向两种结果算出的偶极矩相差很大（法一：将所有单个原子的z方向的偶极矩相加，法二直接看z方向偶极矩）。


variable ban atom 1
compute mychunk ban chunk/atom v_ban nchunk once ids once
compute dipole ban dipole/chunk mychunk
variable dipole_x equal c_dipole[1][1]
variable dipole_y equal c_dipole[1][2]
variable dipole_z equal c_dipole[1][3]
fix 6 ban ave/correlate 5 1000 1000 v_dipole_x v_dipole_y v_dipole_z type auto/upper ave running
fix 7 ban ave/time 1000 1 1000 v_dipole_x v_dipole_y v_dipole_z file dipole.data

xuteng

你好，这个问题解决了吗，我是计算一个水滴的偶极矩，没有思路呢。

xkui

你好，请问这个问题解决了吗？我计算的时候也出现这个问题。而且也出现了和你数据类似的情况，就是某一方向偶极矩值全是正（或负）。感觉是有问题的，这种现象就像是提供了一个统一方向的外加电场。导致方向一致性了。。但是模拟的时候并没有加电场呀。我感觉数据应该是每个方向都是在0上下震荡的，有正有负的。。

yanghaifei

我也遇到了这个问题 请问怎么解决的