MDAnalysis同一Universe多Reader导出轨迹坐标异常解决方法

MDAnalysis 2.2.0版本的ChainReader（多轨迹串联加载时的默认Reader实现）存在坐标缓存同步逻辑缺陷：当你直接遍历u.trajectory.readers获取子Reader、再迭代reader.trajectory时，除第一个默认激活的Reader外，其余子Reader的逐帧迭代不会主动更新主Universe绑定的原子位置缓存，protein.atoms持有的坐标始终停留在该子Reader被访问时的第一帧数值，最终导出的轨迹所有帧坐标完全重复。原始轨迹本身无损坏，问题出在API调用方式没有触发正确的坐标同步。

不需要重新加载拓扑或轨迹文件，完全基于已加载全部3条轨迹的Universe对象操作，通过ChainReader的帧切片机制自动处理跨轨迹的坐标同步，运行效率最高：

import MDAnalysis as mda
u = mda.Universe("11159_dyn_117.pdb", "11156_trj_117.xtc", "11157_trj_117.xtc", "11158_trj_117.xtc")
protein = u.select_atoms("protein")
protein.write("protein.pdb")

current_frame = 0
for traj_idx, sub_reader in enumerate(u.trajectory.readers, 1):
    traj_frame_num = len(sub_reader)
    with mda.Writer(f"{traj_idx}.xtc", protein.n_atoms) as xtc_writer:
        # 按子轨迹的帧范围切片总串联轨迹，自动触发坐标同步
        for ts in u.trajectory[current_frame : current_frame + traj_frame_num]:
            xtc_writer.write(protein.atoms)
    current_frame += traj_frame_num

• 全程复用初始加载的Universe对象，没有重复IO开销，内存占用最低
• 走ChainReader原生的帧切换逻辑，避免直接操作子Reader带来的缓存不同步问题，导出的坐标和原始轨迹完全一致
• 代码兼容MDAnalysis 2.0~2.7全系列稳定版本，不需要升级或降级包版本

• 将导出的单条xtc和之前保存的protein.pdb加载为新Universe，逐帧比对原子位置数组和原始轨迹对应帧的蛋白原子坐标，无偏差
• 用PyMOL等可视化工具加载导出轨迹，可观察到蛋白构象逐帧正常变化，无全帧静止的异常

内容的提问来源于stack exchange，提问作者Francho Nerín Fonz