摘要 分子动力学 (MD) 模拟对于预测不同分子体系的物理和化学性质至关重要。虽然全原子 (AA) MD 提供了高精度,但其计算成本高昂,这促使了粗粒度 MD (CGMD) 的发展。CGMD 将分子结构简化为具有代表性的微珠,以降低成本,但会牺牲精度。像 Martini3 这样的 CGMD 方法,经过实验数据校准后,在各个分子类别中具有良好的泛化能力,但往往无法满足特定领域应用的精度要求。本研究引入了一种基于贝叶斯优化的方法来优化 Martini3 拓扑结构,使其能够适应特定应用,从而确保精度和效率。优化后的 CG 势能适用于任何聚合度,提供与 AA 模拟相当的精度,同时保持与 CGMD 相当的计算速度。通过弥合效率和精度之间的差距,该方法推动了多尺度分子模拟的发展,使各个科学技术领域能够以经济高效的方式发现分子。 1. 引言粗粒度分子动力学 (CGMD) 1,2 已成为材料开发的重要工具,为了解聚合物 3 、蛋白质 4 和膜 5 等复杂分子系统提供了关键信息。CGMD 的主要优势在于它能够在更大长度尺度和更长时间范围内探索分子现象,超越了传统全原子分子动力学 (AAMD) 6–8 模拟的能力,后者通常提供更高的分辨率,因此特别擅长捕捉详细的界面相互作用 9 。具体而言,CGMD 通过将原子团有效地表示为珠子 10–15 来实现这种加速,从而将模拟能力在时间上从皮秒扩展到微秒,在空间上从纳米扩展到微米。因此,粗粒度技术为传统 AAMD 无法获得的复杂分子现象提供了前所未有的洞察,从而能够研究聚合物自组装行为等复杂现象 16 。新兴的CGMD建模工具集依赖于两个关键组件来学习潜在的分子间关系:珠子映射方案和珠子间相互作用的参数化。这些组件的开发主要采用两种方法:自上而下10–12和自下而上13–
主要关键词