结合案例图分析，我把这个项目的完整过程、参数设置、结果解读和踩坑经验整理如下。
客户要做的是磷脂双分子层（lipid bilayer）的粗粒度模拟，目的是快速获得膜的密度分布和有序参数，用于理解膜结构和流动性。全原子模拟成本高，粗粒度（MARTINI）可以把时间尺度放大到微秒级。
我用的软件是 GROMACS + MARTINI 2.0 力场。体系是 DPPC 或 POPC 双分子层，包含约 200 个脂质分子和水珠子。先做能量最小化，然后 NPT 平衡 500 ns，最后跑 2 µs 生产轨迹。
参数设置：时间步长 20 fs（MARTINI 常用 20–25 fs），温度 323 K，压力 1 bar，半各向同性压力耦合；脂质头部和尾部用不同珠子类型表示；每 10 ps 保存一帧。
第一张案例图是双分子层沿膜法线方向（z）的归一化密度分布。不同颜色代表不同组分：水（青色）、磷酸头部（橙色）、甘油（绿色）、酰基链 CH₂（深蓝）、末端 CH₃（红色）。可以看到膜中心是疏水尾部密集区，两侧是水-头部界面。膜厚度（水-水距离减去两侧水渗透深度）约 3.8 nm，和实验值 3.9 nm 很接近，说明粗粒度参数化合理。
第二张案例图是酰基链的序参数 |S_CD| 随碳位置 i 的变化。|S_CD|≈0.20 出现在链中部（i≈6–10），接近界面处稍高（i≈2–4），链末端（i>12）下降。这个分布说明链中部有一定无序性（液晶态），而靠近甘油骨架的部分更有序；末端 CH₃ 完全无序。数值 0.20 与实验 NMR 测得的 DPPC 序参数一致。
第一，MARTINI 的温度需要校准。标准 MARTINI 2.0 对 DPPC 的玻璃化转变温度描述偏低，通常要把模拟温度提高 10–15 K 才能和实验对应。我用 323 K 对应实验 310 K 体系，结果合理。
第二，水珠子数量要足够。我最初按 1:10 水合比，膜发生塌缩；提到 1:20 后膜才稳定。
第三，粗粒度不能看氢键细节。如果客户要解释特定氢键，必须回全原子模型做回补映射（back-mapping）。
下次再做膜体系粗粒度模拟，我会把”密度分布+序参数+面积/脂质”作为标准输出。MARTINI 适合看大趋势和长时间尺度，但不适合定量氢键或精确扩散系数。我会在报告里明确说明 coarse-grained 的适用范围。
计算工程师团队：500+ ｜ 成功计算案例：5W+ ｜ 1W+ 助力顶刊发布经验 需要模拟代算服务可以联系我们，可免费评估项目的计算费用和免费出计算方案。
