lammps中class2条件下的pair_coeff参数转换成cut下的pair_coeff参数需要怎么转换
时间: 2024-09-20 18:03:43 浏览: 36
Fe_mm.eam_cvff_lammps的cvff力场参数文件_分子动力学_
5星 · 资源好评率100%
在LAMMPS中,从`lj/cut/coul/class2`(通常用于分子模拟)切换到`lj/cut/coul/cut`(主要用于金属或其他硬核物质),`pair_coeff`参数的转换主要涉及两部分:
1. ** LJ (Lennard-Jones) 项**:
- `lj/cut/coul/class2`通常会有一个单一的`sigma`(直径)和`epsilon`(势能参数)。对于`lj/cut/coul/cut`,这两个参数保持不变,因为它们描述的是分子间的相互作用强度。
2. ** Coulomb (库仑) 项**:
- 在`lj/cut/coul/class2`中,库仑力可能包括短程排斥和长程屏蔽(通过`kspace_style`控制)。如果使用的是`lj/cut/coul/cut`,则需要明确指定每个粒子对的库仑截断半径(`cutoff`)。
- 如果`lj/cut/coul/class2`中没有明确指定长程库仑截断半径 (`cutoff2`),它会默认等于`cutoff`。
- 要从`class2`转换到`cut`,只需确保提供每个粒子对的单独截断半径,比如:
```python
# 原始 class2 类型
pair_coeff 1 1 100.0 3.5 # 这里假设只有一个截断距离
# 转换为 cut 类型
pair_style lj/cut/coul/cut 10.0 # 设置单个截断半径
pair_coeff 1 1 100.0 3.5 10.0 # 现在需要为1-1配对指定单独的cutoff
```
如果你的`lj/cut/coul/class2`配置中有长程库仑作用并且使用了`kspace_style`(如PPP或Ewald),那么在切换到`lj/cut/coul/cut`时,你还需决定如何处理长程效应,因为`lj/cut/coul/cut`不支持长程屏蔽。
阅读全文