第一章：GROMACS入门——分子动力学到底在干什么？

大家好，欢迎来到这门课。我做了十年分子模拟，头一回接触GROMACS的时候，说实话也挺懵的。一堆命令、各种输入文件，完全不知道从哪里下手。但别担心，这章咱们就把地基打牢。

说白了，分子动力学（MD）就是让原子动起来，然后看它们怎么相互作用。你想想看，在纳米尺度下观察蛋白质怎么折叠、药物分子怎么钻进靶点，这就是MD干的事。

1.1 什么是分子动力学？

先回答最基础的问题：MD是怎么工作的？

• 每个原子被看作经典粒子（忽略量子效应）
• 原子间有相互作用力（力场描述）
• 牛顿第二定律：F = ma，算加速度，更新位置
• 循环往复，得到轨迹

时间步长通常设2飞秒（fs）。为什么是2fs？因为氢原子的振动周期大约10fs，步长太大体系会炸，太小又算不动。我踩过这个坑——曾把步长设为4fs试跑一个小体系，结果能量直接飞了。嗯，这里要注意。

1.2 为什么选GROMACS？

市面上MD软件不少：AMBER、NAMD、CHARMM、LAMMPS……为什么我偏爱GROMACS？

特性	GROMACS	其他软件
运行速度	极快（GPU优化好）	各有千秋
力场支持	广泛（OPLS、AMBER、CHARMM等）	通常专一
分析工具	内置上百个工具	常需额外脚本
学习曲线	中等	普遍偏陡

我在项目中处理过百万原子级的膜蛋白体系，GROMACS能在单张GPU上跑到每天数百纳秒。这种效率，其他软件很难比。

1.3 安装GROMACS——一个让新手头疼的环节

实话实说，我第一次装GROMACS花了整整一下午。现在有conda了，方便很多。但如果你要用GPU加速，还是老老实实源码编译吧。

1.3.1 Conda安装（最快方式）

conda create -n gromacs -c conda-forge gromacs
conda activate gromacs

一行搞定。适合只想跑跑小体系练手的朋友。

1.3.2 源码编译（推荐生产环境）

我个人习惯用CMake编译，控制力更强。

wget https://ftp.gromacs.org/pub/gromacs/gromacs-2024.1.tar.gz
tar xf gromacs-2024.1.tar.gz
cd gromacs-2024.1
mkdir build && cd build
cmake .. -DGMX_BUILD_OWN_FFTW=ON -DGMX_GPU=CUDA -DCUDA_TOOLKIT_ROOT_DIR=/usr/local/cuda
make -j8
make install

1.4 GROMACS的核心概念——文件类型与命令体系

GROMACS有一套自己的文件格式。做模拟，你得先搞懂这些。

• .gro：坐标文件，存原子位置和速度
• .top：拓扑文件，描述原子类型、键、角度、二面角
• .itp：包含拓扑文件（常用于分子片段）
• .mdp：参数文件，控制模拟条件（温度、压力、步长等）
• .tpr：合并后的可执行文件，GROMACS直接读它跑
• .xtc / .trr：轨迹文件，存运动路径（.xtc压缩版更常用）
• .edr：能量文件，存温度、压力、势能等

1.5 GROMACS命令：你最先接触的四个gmx命令

GROMACS 5.0之后的命令统一以gmx开头。这里列出最常用的四个：

1. gmx pdb2gmx：从PDB文件生成拓扑和坐标文件
2. gmx editconf：设置模拟盒子大小
3. gmx solvate：往盒子中加溶剂（水）
4. gmx grompp：把拓扑、坐标、参数打包成.tpr文件

1.6 第一个完整模拟流程

跑一个水中的溶菌酶（lysozyme）模拟，看看完整的命令链长啥样：

# 1. 生成拓扑
gmx pdb2gmx -f 1AKI.pdb -o processed.gro -water spce

# 2. 定义盒子
gmx editconf -f processed.gro -o box.gro -c -d 1.0 -bt cubic

# 3. 加溶剂
gmx solvate -cp box.gro -cs spc216.gro -o solv.gro -p topol.top

# 4. 加离子（电荷中和）
gmx grompp -f ions.mdp -c solv.gro -p topol.top -o ions.tpr
gmx genion -s ions.tpr -o solv_ions.gro -p topol.top -pname NA -nname CL -neutral

# 5. 能量最小化
gmx grompp -f em.mdp -c solv_ions.gro -p topol.top -o em.tpr
gmx mdrun -v -deffnm em

# 6. NVT平衡
gmx grompp -f nvt.mdp -c em.gro -r em.gro -p topol.top -o nvt.tpr
gmx mdrun -v -deffnm nvt

# 7. NPT平衡
gmx grompp -f npt.mdp -c nvt.gro -r nvt.gro -t nvt.cpt -p topol.top -o npt.tpr
gmx mdrun -v -deffnm npt

# 8. 成品模拟（MD）
gmx grompp -f md.mdp -c npt.gro -t npt.cpt -p topol.top -o md.tpr
gmx mdrun -v -deffnm md

代码有点长，但别慌。你只需要理解每一步在干什么就行。后面的章节会一个个拆开讲。

1.7 避坑指南：新手最常见的三个错误

我见过太多人倒在起跑线上了。这里说三个最常见的坑：

• 用了错误的力场：pdb2gmx时选错力场，后面全白搭。仔细看输出日志。
• 盒子太小：溶质离盒壁至少1.0 nm，不然周期性镜像会相互作用。我用过一个案例，蛋白质自身的镜像在模拟中挨在了一起——结果完全不能用。
• 忽略电荷中和：体系净电荷不为零，模拟会崩溃。gmx genion加离子时别忘了-neutral参数。

好了，这一章的内容就是这些。把基础概念搞明白，把GROMACS装上，然后找个结构文件试试手。别怕出错，出错了再回来翻这章内容就行。