几类溶剂分子的表面张力和黏度的分子动力学研究.pdfVIP

摘要

中文摘要

很多化学反应和物理化学过程的发生都需要特定的溶剂环境。比如日用化学

品行业的增泡剂和增粘剂、清洗行业的清洁剂、化工行业的树脂和乳液、石油等

液体的输送、润滑剂的润滑过程等，都与溶剂体系的表面张力、黏度等物理化学

指标有着密切的关系。长期以来，人们需要通过实验手段测定某种溶剂分子或混

合溶剂的表面张力和黏度。尽管实验室测定方法在确保数据准确性和可靠性方面

具有显著优势，并且对于科学研究和产品质量控制等方面至关重要，但同时也需

要注意到，这种方法的实施可能涉及到较为复杂的过程，以及对专业技术和设备

的高要求。随着计算机计算能力的不断增加，现在人们可以通过分子动力学模拟

来得到某一种溶剂或者混合溶剂的表面张力、黏度等物理化学性质。但是，对于

这些物理化学指标，尚缺乏标准化的模拟方法。分子动力学模拟是一种强大的计

算工具，可用于研究物质的各种性质。本文着重探讨了利用分子动力学模拟计算

表面张力和黏度的方法。首先介绍了分子动力学模拟的基本原理，包括

GROMACS原理、力场的选择以及模拟过程中的积分算法等。接着详细阐述了

计算表面张力和黏度的方法。最后，讨论了分子动力学模拟在计算表面张力和黏

度方面的优势与局限性，并对未来的研究方向进行了展望。

在本研究中，我们使用GROMACS分子动力学软件包，应用CHARMM、

OPLS-AA、GAFF三种常见的小分子力场，计算了苯系物、醇类、烷烃类等多种

分子的表面张力和黏度。探索了适用于多种溶剂分子体系的分子动力学模拟计算

ZNPT

表面张力和黏度的一般流程。我们通过对初始结构能量最小化，延长轴，

模拟，NVT模拟后收集数据。首先就目前的计算结果，我们发现在表面张力的

计算中，三种力场都可以很好的用于计算链状烷烃，环己烷和正十六烷表现出相

对较大的差值，对于苯系物CHARMM力场表现出较好的一致性，对于醇类，

CHARMM力场模拟的数值与实验值的差距逐渐收敛，而二氯甲烷，氯仿和二甲

基亚砜结果却没有规律可循。在黏度的计算中，从结果看，模拟值与实验值基本

吻合。OPLS力场在计算二氯甲烷、氯仿的黏度值时比较具有优势，其他分子的

计算还需要进一步研究，CHARMM力场在计算醇类和二甲基亚砜的黏度值时比

I

摘要

较具有优势，其他分子比如烷烃卤代烃，OPLS力场可能更有优势。不同力场对

同一物质的黏度预测存在差异，这和力场特性和适用范围有关，所以力场参数的

优化对提高计算的准确性至关重要。不同力场对同一物质的表面张力预测也存在

差异，这和力场在分子间相互作用的描述有关。

综上所述，本研究为应用分子动力学手段模拟溶剂的表面张力和黏度提供了

可行的模拟方法，虽然还存在诸多的问题，我们也希望可以为今后通过分子动力

学手段模拟溶剂分子的物理化学性质提供系统性的工作基础。

关键词：

分子动力学模拟，表面张力，黏度，GROMACS模拟

II

Abstract

Abstract

Manychemicalreactionsandphysicochemicalprocessesrequireaspecific

solventenvironment.Forexample,foamingagentsandtackifiersinthedaily

chemicalsindustry,detergentsinthecleaningindustry,resinsandemulsionsinthe

chemicalindustry,transportationofliquidssuchaspetroleum,andlubrication

processesoflubricantsareallclosel