進(jìn)行l(wèi)ammps模擬，如果初始模型不合理、原子間距離過(guò)近，容易導(dǎo)致體系能量過(guò)高，原子速度過(guò)快致使原子飛出box。
前文介紹過(guò)“delete_atoms”命令可刪除重疊或距離過(guò)近原子，但在大多數(shù)模型中，并不需要?jiǎng)h除原子，有些模型也不允許刪除原子，這時(shí)候就需要用到minimize命令。
1. minimize命令
minimize命令通過(guò)力的計(jì)算排開(kāi)距離過(guò)近的原子，使體系的能量降到一個(gè)較低的水平，便于后續(xù)的弛豫操作。
minimize語(yǔ)法為：

minimize etol ftol maxiter maxeval

后面四個(gè)參數(shù)為minimize運(yùn)行停止條件，其中etol、ftol分別為能量和力的容差，maxiter、maxeval為迭代最大步數(shù)。
在最小化過(guò)程中，達(dá)到四個(gè)條件中的任意一個(gè)時(shí)，最小化即停止運(yùn)行。
力或者能量已經(jīng)降低到一個(gè)較低水平，或者力或能量還沒(méi)有降低到一個(gè)低值，但設(shè)定的迭代步數(shù)已經(jīng)用完，都會(huì)停止最小化。
2. min_style命令
minimize最小化采用的算法由min_style命令確定，min_style語(yǔ)法為：

min_style style style = cg or hftn or sd or quickmin or fire or fire/old or spin or spin/cg or spin/lbfgs

最常用的為共軛梯度法（cg）、最速下降法（sd）。
3. 使用舉例

min_style cg minimize 1.0e-4 1.0e-6 100 1000

如果已經(jīng)進(jìn)行了最小化，但是后續(xù)馳豫仍然會(huì)報(bào)錯(cuò)，可查看log文件，排查minimize運(yùn)行停止的條件，適當(dāng)調(diào)小能量和力的容差或者加大迭代步數(shù)。
如以上語(yǔ)句可改為：

min_style cg minimize 1.0e-15 1.0e-15 10000 10000

mimimize可有效降低出錯(cuò)概率，誰(shuí)用誰(shuí)知道，用了都說(shuō)好的一個(gè)命令。

更多l(xiāng)ammps教程，關(guān)注微信公眾號(hào)：lammps加油站

總結(jié)

以上是生活随笔為你收集整理的lammps教程：minimize命令使用详解的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。