膜脂有序参数（Order Parameter S_cd）绘图工具
膜脂有序参数（Order Parameter S_cd）绘图工具order_param · Trajectory Dynamics
概述
有序参数 Scd（deuterium order parameter）刻画脂酰链每个碳位上 C–H 键相对膜法向的有序程度：Scd = |⟨3cos²θ−1⟩/2|，θ 为 C–D/C–H 键与膜法向夹角。数值 0 = 完全无序（各向同性），~0.35 = 高度有序（全反式伸直链）。它是膜模拟中与²H NMR 实验直接可比的核心观测量。
本工具解析 gmx order 输出，绘制 Scd 随碳位（C2–C16 等）分布曲线。sn-1 与 sn-2 两条链可叠画对比；小分子/蛋白插入或温度变化前后的 Scd 曲线对比，即“扰动如何沿链传播”的标准证据图。
输入数据格式
| 格式 | 说明 |
|---|---|
| GROMACS .xvg | gmx order 输出：碳位序号 + S_cd（sn-1/sn-2 各一份数据集） |
| 多体系对比 | 纯膜 vs 含药物/蛋白膜各一份 .xvg 叠画 |
| 兼容格式 | .xvg / .txt / .csv / .zip |
结果解读
典型 POPC 链的 Scd 剖面：C2–C8 平台区（0.18–0.22，链有序段），向链尾单调下降至 <0.05（末端甲基自由旋转）。小分子嵌入通常压低头部区平台（增加无序）；胆固醇则显著抬升整链（有序化）——方向相反的两类效应在曲线上方向明确。sn-1/sn-2 的系统性差异属链化学结构差异，正常。
对比实验需注意模拟口径：力场（Berger vs Lipid17/CHARMM36）的系统差可达 ±0.03；C2 位与甘油骨架耦合常呈异常跳点，文献惯例从 C3 起比较。叠加²H NMR 数据点对照时纵轴一致（同 Scd 定义）。
使用步骤
- 为脂酰链生成碳位索引（gmx order 文档的 make_ndx 流程）。
- 运行
gmx order -od deuter分别对 sn-1/sn-2。 - 拖入 .xvg，两条链自动分色。
- 导出 600 dpi PNG / 矢量 PDF。
所有解析与渲染均在浏览器本地完成：上传入口就是文件选择器本身，数据不经过任何服务器， 断网环境同样可用。可以先用本页工作台的“一键载入示例数据”体验全流程，再换自己的模拟输出。
常见问题
S_cd 怎么和 NMR 实验对比？
²H NMR 的四极裂分 Δν_Q 与 S_cd 直接相关（Δν_Q = (3/4)·(e²qQ/h)·S_cd），文献常把实测 S_cd 直接标在图上对比。碳位对齐与链归属（sn-1/sn-2）先核对，混合位置的比较是常见错误。
为什么 C2 位的 S_cd 总是异常？
C2 与甘油骨架的二面角耦合，其 C–H 键取向统计偏离简单的链内假设，多数力场在此位置给出系统偏差。惯例从 C3 开始与实验/NMR 对比。
小分子使膜“变软”还是“变硬”怎么看？
看曲线方向：S_cd 全链压低 = 流动化（软化，如短链醇、麻醉剂类）；抬升 = 有序化（胆固醇及其类似物）。只压低头部区而链尾不变的局部效应，提示分子嵌在头部-界面区。
