LAMMPS – 晨阳博客 https://zhangchenyang.cn 欢迎来到张晨阳的博客 Mon, 12 Apr 2021 08:58:44 +0000 zh-Hans hourly 1 https://wordpress.org/?v=6.9.4 [LAMMPS命令合集]units命令 https://zhangchenyang.cn/151.html Mon, 12 Apr 2021 08:55:51 +0000 http://zhangchenyang.cn/?p=151 继续阅读 [LAMMPS命令合集]units命令]]> units是LAMMPS中最常用的命令之一，其用途为设置模拟过程的单位

语法

units stytle

style = lj or real or metal or si or cgs or electron or micro or nano

例子

units metal
units lj

用法

style为real时

mass = grams/mole 质量=克/摩尔
distance = Angstroms 距离=埃米（10-9米）
time = femtoseconds 时间=飞秒（10-15秒）
energy = Kcal/mole 能量=千卡/摩尔
velocity = Angstroms/femtosecond 速度=埃米/飞秒
force = Kcal/mole-Angstrom 力=能量/距离
torque = Kcal/mole 力矩=千卡/摩尔
temperature = Kelvin 温度=开尔文单位
pressure = atmospheres 压力=大气压
dynamic viscosity = Poise 动力学黏度=帕
charge = multiple of electron charge (1.0 is a proton) 电荷=电子电荷倍数（1.0就是一个质子）
dipole = charge*Angstroms 偶极子=电荷*埃米
electric field = volts/Angstrom 电场=伏/埃米
density = gram/cm^dim 密度=克/平方厘米
]]> LAMMPS模拟命令顺序 https://zhangchenyang.cn/144.html Mon, 12 Apr 2021 08:00:39 +0000 http://zhangchenyang.cn/?p=144 继续阅读 LAMMPS模拟命令顺序]]> LAMMPS模拟过程涉及较多的命令，一般来说主要分为以下几个流程

/////////////////

//初始参数

/////////////////

 

//模型初始化

units

boundary

atom_style

 

//设置键势类型

bond_style

//设置键角类型

angle_style

//设置二面角势类型

dihederal_style

//设置对势类型

pair_style

//读取结构数据

read_data

 

//设置邻近原子参数

neighbor

neigh_modify

 

/////////////////

//计算区域

/////////////////

//设置多久输出一次数据

//弛豫

nve

//温度设置

 

//设置系综

nve nvt npt

 

//热力学输出（显示在屏幕上）

thermo_style

thermo

 

//设置步长

timestep

 

//*将步数的计数器设置为指定值，设为0即清空

reset_timestep 0

 

//开始计算

run 1000（计算多少步，模拟时长为 步数*时间步长）
]]> [转]region命令 https://zhangchenyang.cn/139.html Mon, 12 Apr 2021 07:08:57 +0000 http://zhangchenyang.cn/?p=139 继续阅读 [转]region命令]]> region命令用于定义一个空间几何区域。

使用语法:

region ID style args keyword arg ...
  • ID = 待定义区域的ID
  • style = delete / block / cone / cylinder / plane / prism / sphere / union / intersect
delete = no args
  block args = xlo xhi ylo yhi zlo zhi
    xlo,xhi,ylo,yhi,zlo,zhi = 各个维度上的范围值
  cone args = dim c1 c2 radlo radhi lo hi
    dim = x or y or z = 圆锥台轴线方向
    c1,c2 = 圆锥台轴线在另外两个维度上的坐标值
    radlo,radhi = 圆锥台上下底面的半径
    lo,hi = 圆锥台轴线方向的范围值
  cylinder args = dim c1 c2 radius lo hi
    dim = x or y or z = 圆柱体轴线方向
    c1,c2 = 圆柱体轴线在两位两个维度上的坐标值
    radius = 圆柱体半径，可以使用变量
    lo,hi = 圆柱体轴线方向的范围值
  plane args = px py pz nx ny nz
    px,py,pz = 平面上一个点的坐标
    nx,ny,nz = 平面法矢坐标
  prism args = xlo xhi ylo yhi zlo zhi xy xz yz
    xlo,xhi,ylo,yhi,zlo,zhi = 平行六面体在各个维度上的范围值
    xy = y方向在x方向的倾斜值
    xz = z方向在x方向的倾斜值
    yz = z方向在y方向的倾斜值
  sphere args = x y z radius
    x,y,z = 球心坐标
    radius = 球体半径，可以使用变量
  union args = N reg-ID1 reg-ID2 ...
    N = 要进行并集操作的区域的数目，必须大于等于2
    reg-ID1,reg-ID2, ... = 要进行并集操作的区域的ID
  intersect args = N reg-ID1 reg-ID2 ...
    N = 要进行交集操作的区域的数目，必须大于等于2
    reg-ID1,reg-ID2, ... = 要进行交集操作的区域的ID
  • keyword = side or units or move or rotate
side value = in or out
    in = 指定几何体内侧作为区域
    out = 指定几何体外侧作为区域
  units value = lattice or box
    lattice = 以晶格距离作为几何距离单位
    box = 以模拟盒子作为几何距离单位
  move args = v_x v_y v_z
    v_x,v_y,v_z = x, y, z 方向的位移量，equal类型的变量
  rotate args = v_theta Px Py Pz Rx Ry Rz
    v_theta = 旋转角，equal类型的变量，弧度单位
    Px,Py,Pz = 旋转操作的中心点
    Rx,Ry,Rz = 旋转轴矢量

使用举例

region 1 block -3.0 5.0 INF 10.0 INF INF
region 2 sphere 0.0 0.0 0.0 5 side out
region void cylinder y 2 3 5 -5.0 EDGE units box
region 1 prism 0 10 0 10 0 10 2 0 0
region outside union 4 side1 side2 side3 side4
region 2 sphere 0.0 0.0 0.0 5 side out move v_left v_up NULL
]]>