在LAMMPS中,单位制定义了模拟中使用的物理量的单位,包括距离、时间、质量、能量等。LAMMPS支持多种单位制,用户可以根据研究的物理系统选择合适的单位制。
以下是LAMMPS中常用的单位制:
1. 单位制概览 (units
命令)
在LAMMPS脚本中,可以通过以下命令选择单位制:
units <style>
其中 <style>
是以下几种单位制之一:
单位制名称 | 主要应用领域 | 特点 |
---|---|---|
lj | Lennard-Jones势模型、无量纲模拟 | 无量纲单位,适用于分子动力学简单模型 |
real | 化学、软物质、分子模拟 | 使用常见的实验单位,如Å、fs、kcal/mol |
metal | 金属、合金、固体材料模拟 | 使用金属领域的单位,如Å、ps、eV |
si | 工程学仿真、连续介质力学 | 国际单位制(SI单位),如m、kg、J |
cgs | 历史上使用的厘米-克-秒单位制 | 厘米-克-秒(CGS)系统 |
electron | 电子学、量子点、等离子体 | 用于电子物理模拟,基于原子单位制 |
micro | 微米级系统、微流控、MEMS | 微米级长度、毫秒时间单位 |
nano | 纳米级系统、纳米材料 | 纳米级长度、纳秒时间单位 |
2. 各单位制详细信息
以下是所有单位制的物理量对应单位的具体说明:
(1) Lennard-Jones 单位制 (lj
)
- 特点:无量纲,适用于简化的模型研究,如Lennard-Jones势。
- 单位定义:
- 距离: σ \sigma σ(Lennard-Jones参数)
- 能量:ϵ(Lennard-Jones参数)
- 时间: ( m σ 2 / ϵ ) 1 / 2 (m\sigma^2/\epsilon)^{1/2} (mσ2/ϵ)1/2
- 质量:m(粒子质量)
- 温度: ϵ / k B \epsilon/k_B ϵ/kB
- 压力: ϵ / σ 3 \epsilon/\sigma^3 ϵ/σ3
- 动量: m ⋅ σ / t m \cdot \sigma / t m⋅σ/t
(2) 实验单位制 (real
)
- 特点:适用于化学、软物质、分子模拟等实验相关的研究。
- 单位定义:
- 距离:Å(埃,1Å = 1 0 − 10 10^{-10} 10−10 m)
- 时间:fs(飞秒, 1 0 − 15 10^{-15} 10−15 秒)
- 质量:amu(原子质量单位)
- 能量:kcal/mol
- 温度:K(开尔文)
- 压力:atm(大气压,1 atm = 101325 Pa)
(3) 金属单位制 (metal
)
- 特点:用于金属材料模拟,特别是涉及电子伏特(eV)的计算。
- 单位定义:
- 距离:Å
- 时间:ps(皮秒, 1 0 − 12 10^{-12} 10−12秒)
- 质量:amu
- 能量:eV(电子伏特, 1 e V ≈ ( 1.602 × 1 0 − 19 ) J 1 eV ≈ (1.602 \times 10^{-19}) J 1eV≈(1.602×10−19)J)
- 温度:K
- 压力:bars( 1 b a r = 1 0 5 P a 1 \ bar = 10^5 \ Pa 1 bar=105 Pa)
(4) 国际单位制 (si
)
- 特点:工程仿真中广泛使用的单位,基于米-千克-秒。
- 单位定义:
- 距离:m(米)
- 时间:s(秒)
- 质量:kg(千克)
- 能量:J(焦耳)
- 温度:K
- 压力:Pa(帕斯卡)
(5) CGS 单位制 (cgs
)
- 特点:基于厘米-克-秒,历史上常用于流体力学研究。
- 单位定义:
- 距离:cm(厘米)
- 时间:s(秒)
- 质量:g(克)
- 能量:erg(尔格, 1 e r g = 1 0 − 7 J 1 \ \ erg = 10^{-7}J 1 erg=10−7J)
- 温度:K
- 压力:dyne/cm²
(6) 电子单位制 (electron
)
- 特点:用于电子学、等离子体领域的模拟,基于原子单位制。
- 单位定义:
- 距离:Å
- 时间:fs
- 质量:电子质量 m e = 9.109 × 1 0 − 31 k g m_e = 9.109 \times 10^{-31}kg me=9.109×10−31kg
- 能量:eV
- 温度:K
- 压力:Pa
(7) 微米单位制 (micro
)
- 特点:适用于微米尺度下的流体、微纳机电系统(MEMS)仿真。
- 单位定义:
- 距离:μm(微米, 1 0 − 6 10^{-6} 10−6米)
- 时间:ms(毫秒, 1 0 − 3 10^{-3} 10−3秒)
- 质量:pg(皮克克, 1 0 − 12 10^{-12} 10−12克)
- 能量:pJ(皮焦耳)
- 温度:K
- 压力:Pa
(8) 纳米单位制 (nano
)
- 特点:用于纳米级系统模拟。
- 单位定义:
- 距离:nm(纳米, 1 0 − 9 10^{-9} 10−9米)
- 时间:ns(纳秒, 1 0 − 9 10^{-9} 10−9秒)
- 质量:attogram( 1 0 − 18 10^{-18} 10−18克)
- 能量:eV
- 温度:K
- 压力:Pa
3. 选择单位制的注意事项
- 与物理模型匹配:选择单位制时,应确保其与研究系统的物理模型和参数吻合。例如:
- 金属材料模拟:推荐使用
metal
单位。 - 化学系统或生物分子:推荐使用
real
单位。 - 无量纲模拟:使用
lj
单位。
- 金属材料模拟:推荐使用
- 单位转换:如果实验数据或文献参数的单位与LAMMPS的单位不匹配,需要在脚本中进行单位转换。
- 一致性:所有输入参数(如力场参数、温度、盒子尺寸等)必须与选择的单位制一致。例如,在
real
单位下,能量单位是 kcal/mol,不能直接使用 eV。
4. 示例
(1) 使用 real
单位
units real
lattice fcc 4.0
pair_style lj/cut 10.0
pair_coeff 1 1 0.1 3.5
- 单位:Å、fs、kcal/mol、atm。
- Lennard-Jones势参数 ϵ = 0.1 k c a l / m o l , σ = 3.5 A ˚ \epsilon = 0.1 \ kcal/mol,\sigma = 3.5\ Å ϵ=0.1 kcal/mol,σ=3.5 A˚)。
(2) 使用 metal
单位
units metal
lattice bcc 2.87
pair_style eam
pair_coeff * * Fe.eam.alloy Fe
- 单位:Å、ps、eV、bar。
- 晶格常数 a = 2.87 A ˚ a = 2.87 \ Å a=2.87 A˚。
5. 总结
LAMMPS支持多种单位制,以满足不同研究领域的需求。选择正确的单位制对于构建准确的模拟模型至关重要。常见单位制包括:
lj
(无量纲)用于简化模型。real
用于分子模拟。metal
用于固体材料模拟。si
和cgs
用于工程学应用。