lammps教程：如何計算單個原子的體積
發(fā)布時間：
大家好，我是小馬老師。本文介紹如何使用lammps命令計算單原子的體積。lammps提供了compute voronoi/atom命令計算單原子體積。基本
大家好，我是小馬老師。
本文介紹如何使用lammps命令計算單原子的體積。
lammps提供了compute voronoi/atom命令計算單原子體積。
基本原理可參考voronoi算法，根據(jù)voronoi算法，把單個原子所占據(jù)的空間劃分為一個多邊形，也稱為泰森多邊形，多邊形的體積即為該原子的體積。
命令例句：
compute 1 all voronoi/atom
該例句共輸出兩個計算結(jié)果，第一個結(jié)果為單原子的體積，第二個結(jié)果為該原子多面體的面數(shù)，也可以理解為相鄰原子的數(shù)目。
因此，引用單原子體積時，直接使用c_1[1]即可。
原子體積的單位為距離單位的立方，如A3。
配合compute stress命令就可以計算并輸出單原子應力：
compute 1 all stress/atom NULLncompute 2 all voronoi/atom nvariable stressx atom c_1[1]/c_2[1]ndump 1 all custom 1000 dump.xyz id type x y z v_stressx
公眾號：lammps加油站
