LAMMPS（Large-scale Atomic/Molecular Massively Parallel Simulator）是一款专为材料模拟设计的高性能分子动力学软件，支持在并行计算环境中高效运行。其核心目标是实现大规模原子系统的动态行为分析，广泛应用于固体、液体及复杂界面体系的研究。该工具由美国能源部资助，最初由桑迪亚国家实验室开发，现已成为开源社区的重要项目，遵循 GNU GPLv2 许可协议发布。

本示例通过三种不同的指令机制——velocity、fix move 和 fix setforce，演示如何精确控制单个原子的运动轨迹。每种方法对应不同的物理原理与应用场景，适用于不同需求的仿真设计。

一、输入文件配置解析

以下为核心输入脚本 in.txt 的结构化说明：

# 模拟系统设置
units           lj
atom_style      atomic
boundary        f f f
region          box block 0 150 0 300 0 10
create_box      1 box
create_atoms    1 single 30.0 10.0 5.0
create_atoms    1 single 70.0 10.0 5.0
create_atoms    1 single 110.0 10.0 5.0
mass            1 1.0
pair_style      none

• 单位系统：采用无量纲的 Lennard-Jones 单位制，便于跨体系比较。
• 原子样式：使用 atomic 类型，仅记录位置与类型信息。
• 边界条件：所有方向设为固定边界（f），防止原子逃逸或周期性重复。
• 空间区域：定义一个长方体模拟域，尺寸为 (150×300×10)，用于限制原子活动范围。
• 原子创建：在指定坐标处生成三个原子，初始位于 Y=10.0 平面，横向分布以方便观察。
• 质量设定：原子质量为 1.0，在 setforce 中直接影响加速度。
• 相互作用：关闭粒子间势能（none），避免干扰运动控制逻辑。

二、原子运动控制策略

1. 通过 velocity 设置初速度（匀速运动）

group           vel_group id 1
velocity        vel_group set 0.0 50.0 0.0

该方式通过直接赋予初始速度，使原子在后续积分中依据牛顿定律持续运动。由于无外力作用，速度保持恒定，表现为匀速直线运动。

2. 使用 fix move 强制路径控制（强制匀速）

group           move_group id 2
velocity        move_group set 0.0 0.0 0.0
fix             2 move_group move linear 0.0 30.0 0.0

fix move 是一种主动干预机制，它绕过标准动力学积分过程，按指定速率直接更新原子位置。无论系统状态如何，原子均以设定速度沿指定方向移动，适用于理想化路径追踪。

3. 利用 fix setforce 实现加速运动

group           force_group id 3
velocity        force_group set 0.0 0.0 0.0
fix             3 force_group setforce 0.0 20.0 0.0
fix             1 all nve

通过施加恒定外力，结合牛顿第二定律 F = ma，产生恒定加速度。在此例中，加速度为 a = 20.0 / 1.0 = 20.0（LJ 单位）。因此原子将做匀加速运动，位移随时间平方增长。

三、输出与运行参数

dump            1 all atom 1 trajectory.xyz
dump_modify     1 sort id
timestep        0.01
thermo          10
thermo_style    custom step temp ke pe etotal
run             500

• dump 命令每步保存一次原子坐标，生成 trajectory.xyz 文件，可用于可视化。
• dump_modify sort id 确保输出按原子编号排序，提升可读性。
• timestep 设为 0.01，总模拟时间为 5.0 单位时间。
• thermo 输出热力学数据，便于监控能量守恒情况。
• run 500 表示执行 500 步模拟。

四、预期运动结果对比

三种方法分别体现了经典力学中的不同运动规律，适用于从理想运动建模到真实力场响应的多种场景。

五、可视化与验证

模拟结束后，下载 trajectory.xyz 文件，使用 VMD、Ovito 等工具加载，即可直观观察三个原子在时间轴上的运动路径差异。该过程有助于理解各类 fix 命令的本质区别及其在实际仿真实验中的适用性。