破乳剂界面微观特性分子模拟方法

界面微观特性采用全原子分子动力学模拟方法进行研究，所有模拟计算均通过 Gromacs2020.6 程序[120]完成，并以十一烷作为油相。十一烷、聚碳酸醚酯类破乳剂和腰 果酚基两性表面活性剂（CNBS）的分子通过 Gaussian 2016 的 B3LYP/def2SVP 水平和 xtb [115]的 GFN2-xTB 方法相结合的方法进行结构优化。Restrained Electrostatic Potential （RESP）电荷计算由 Gaussian 2016 和 Multiwfn [116]软件在 B3LYP/def2TZVP 水平完成， General Amber force field （GAFF）力场参数采用北京科音自然科学研究中心开发的 Sobtop 软件[121]获得。水分子、Na +、Cl -使用 Amber99sb-ildn 力场参数。分子模拟体系 使用 Packmol 软件[122]构建，盒子尺寸为 10 nm×10 nm×15 nm，包括 2 层油层（10 nm×10 nm×2.5 nm）、2 层表面活性剂层（10 nm×10 nm×2.5 nm）和 1 层水层（10 nm×10 nm×5 nm）。每层油层填充 400 个油分子，每层水层填充 15896 个水分子及 52 个 NaCl 分子， 表面活性剂层填充量根据实验设计设置。整个分子模拟过程的温度保持在 45℃，先进 行能量最小化优化，然后进行 0.1 MPa 压力下进行等温等压（NPT）系综模拟，总模拟 时长为 1 ns，时间步长为 2 fs。