基于不同势函数的正构烷烃物理性质分子动力学模拟.pdf

燃烧科学与技术 Journal of Combustion Science and Technology 2023，29(4)：)：451-459DOI 10.11715/rskxjs.R202305023 收稿日期：2022-03-16 基金项目：国家自然科学基金资助项目(52106154)；江苏省自然科学基金资助项目(BK20190855).作者简介：陈飞龙(1998 )，男，硕士研究生，fl_.通信作者：张延志，男，博士，副教授，.基于不同势函数的正构烷烃物理性质分子动力学模拟 陈飞龙1，张延志1，何志霞2(1.大连理工大学能源与动力学院，大连 116024；2.江苏大学能源研究院，镇江 212013)摘 要：采用分子动力学模拟详细比较了 4 种常用势函数对正构烷烃物理性质预测精度的影响，提出了一种新的修正力场并对其进行广泛的验证，验证主要包括饱和两相密度、饱和蒸气压和自扩散系数等热物理和输运性质之后，将改进力场应用于二元组分燃料热物性研究中，探究了不同混合比例条件下通过每种组分加权平均的方法估算混合物饱和热物性的准确度研究发现，改进力场对正庚烷和正十二烷物理性质的预测与实验值较为吻合；在温度较低、正庚烷摩尔分数较高或者温度较高、正庚烷摩尔分数较低的工况下，通过加权平均得到的饱和物性值较为准确 关键词：分子动力学模拟；物理性质；势函数；二元组分 中图分类号：TK401 文献标志码：A 文章编号：1006-8740(2023)04-0451-09 Molecular Dynamics Simulation of Physical Properties of n-Alkane Based on Different Potential Models Chen Feilong1，Zhang Yanzhi1，He Zhixia2(1.School of Energy and Power，Dalian University of Technology，Dalian 116024，China；2.Institute for Energy Research，Jiangsu University，Zhenjiang 212013，China)Abstract：In this paper，molecular dynamics simulation was used to compare in detail the effect of four com-monly used potential models on the prediction accuracy of the physical properties of n-alkane，and a new modified model was proposed and extensively validated，which mainly includes the thermophysical and transport properties such as saturated two-phase density，saturated vapor pressure and self-diffusion coefficient.After that，the modi-fied potential model was applied to the study of the thermophysical properties of binary-component fuels，and the accuracy of estimating the saturated physical properties of the mixture by the weighted average method of each component at different mixing ratios was explored.It was found that the predictions from the modified model show good agreement with the measurements for n-heptane and n-dodecane.When the temperature is lower and the mole fraction of n-heptane is higher，or when the temperature is higher and the mole fraction of n-heptane is lower，the saturated physical property value obtained by the weighted average is more accurate.Keywords：molecular dynamics simulation；physical property；potential model；binary-component 燃烧科学与技术 第 29 卷 第 4 期 452 步入全球现代化的 21 世纪以来，化石燃料具有潜在的能源短缺的危机，燃烧化石燃料(主要是碳氢燃料)而产生的二氧化碳是加快全球变暖的温室气体的主要来源之一 面对能源短缺和环境恶化的双重危机，节约能源和清洁低碳化发展已成为全球共识 除了发展可再生能源和探索替代能源之外，提高碳氢燃料的燃烧利用效率至关重要 实现这一目标的基础和前提是对碳氢化合物各项物理性质的定量估计，而对于燃料在一些极端工况如高温高压的环境条件下或者是不同组分混合物热力学以及输运性质的实验数据比较匮乏，且实验成本较高，因此用理论预测的手段来补充很有必要 状态方程是预测流体热物性的有效手段，但是先决条件是需要进行各种假设、需要提供经验参数以及混合规则等，存在对于不同环境不同组分的适用性问题，同时计算效率和精度等问题也同样需要考虑 为了避免传统数值模拟和实验所遇到的问题，分子动力学(molecular dynamics，MD)模拟技术近年来被越来越地多应用于流体性质的研究 MD 模拟的优势是不需要引入对所研究物理过程的任何假设，只要有能正确描述物质粒子间作用力的势函数，适用于对各项热力学以及输运性质的预测 MD 计算中，势函数作为描述粒子间作用力近似数学表达方法的选取极为关键 碳氢化合物的常用分子模型主要分为全原子模型和联合原子模型，前者视每个原子为一个位点，后者将每个甲基 CH3和亚甲基 CH2基团分别视为一个联合原子来简化计算 最早出现的有 Jorgensen 等1-4提出的 OPLS-AA和 OPLS-UA 力场，广泛应用于蛋白质、烃类和醇类等复杂的大分子物质 之后 Siepmann 等5-6提出的SKS 力场，更适用于模拟链长较短的烷烃 除此之外，又相继出现了 Martin 等7通过拟合烷烃的临界温度和饱和液相密度的实验数据提出的 TraPPE-UA 力场以及 Nath 等8提出的对短长链烷烃都适用的NERD 力场 邓磊9比较了 OPLS-AA 力场以及一种联合原子力场对正庚烷液滴的蒸发特性的影响，结果表明，在采用截断半径法的计算下所选用的联合原子力场比 OPLS-AA 力场对正庚烷粒子间作用力的描述更优，而且全原子模型计算成本较高，尤其对于长链正构烷烃 然而，到目前为止，详细比较不同常用联合原子势函数对正构烷烃物理性质预测精度的研究比较缺乏 本文选取了内燃机燃料中典型组分正庚烷和正十二烷为研究对象，比较了不同势函数对正构烷烃物理性质(包括饱和及输运物性)的影响，并提出了一种新的修正力场并对其进行了烷烃饱和两相密度、饱和蒸气压和自扩散系数等热物理性质的验证 另外，将改进力场应用于双组分燃料(正庚烷和正十二烷)热物性的研究中，主要探究了不同混合比例条件下通过将每种组分饱和物性按摩尔比例加权平均的方法来估算混合物饱和物性的准确度 1 模拟方法 1.1 势函数模型 分子动力学模拟中，势函数或者力场模型是作为粒子间作用力的近似描述，力场的选取对于整个模拟的过程可谓重中之重，因此在进行模拟前需要选择准确合适的力场 每种物质有各自不同适合于本身的力场参数，同一种物质也有适合于表述各种性质而对应的不同的力场参数 为了节约计算成本，本文采用了常用的 4 种适用于碳氢化合物的联合原子力场模型 OPLS-UA、SKS、TraPPE-UA 以及 NERD 此外，基于 SKS 非键势能、键伸缩势能以及键角弯曲势能和另一种五参数约束的二面角扭曲势能10，本文提出一种新的修正力场(命名为 SKS-New)联 合 原 子 力 场 的 非 键 势 能 可 以 用 截 断 的Lennard-jones12-6 势来描述：126126cc()4=LJijijijijijijijijurrrrr(1)式中：ij和ij分别是尺寸参数和相互作用势阱深度；ijr是 L-J 势为 0 处原子间的距离；cr为截断半径 通过增大截断半径可使 MD 模拟结果更精确，但也会增加计算成本 综合考虑计算精度与成本，本文截断半径取c1.5=rnm 采用 Lorentz-Berthelot 混合规则来确定不同种类原子间相互作用的 L-J 参数：()12ijijijij=+=(2)在最初的 SKS 和 TraPPE-UA 的力场模型中，将所有键视为固定长度的刚性键，然而 LAMMPS 平台不提供正构烷烃等大分子的键长约束，因此本研究中统一使用谐波势来代替刚性键拉伸描述键伸缩势能：eq2br1()2iiuk rr=(3)式中：rk为键伸缩系数；eqir为平衡键长，本研究中除OPLS-UA 力场之外其他力场的键伸缩势能参数统一取自 NERD 模型8 以往的研究表明，对于烷烃的蒸陈飞龙等：基于不同势函数的正构烷烃物理性质分子动力学模拟 燃烧科学与技术 453 发特性模拟，这种改动没有对结果带来差异11 同样用一种谐波势来描述键角弯曲势能：eq21/2()=iiuk(4)式中：k为键弯曲系数；eqi为平衡角 用 OPLS 二面角扭曲势4描述前 4 种力场的二面角扭曲势能：()1211(1 cos2)1cos22=+uVV 31(1cos3)2+V(5)改进力场的二面角扭曲势能则由余弦多项式表示的多谐波势10来描述：511cos=nnnuA(6)此处为相连 4 个联合原子之间的二面角，nV和nA皆为二面角扭曲系数 各联合原子模型力场的 4 项势能参数列于表 1 1.2 热物性统计方法及计算公式 本文将模拟盒子沿 x 方向(即垂直于界面方向)平均划分为sN个切片，由每层切片中的原子速度、原子数、以及原子受力的微观性质，来统计最终气液平衡中均匀两相的温度、饱和密度、饱和压力以及界面层表面张力的空间分布等宏观性质 第k块切片内的分子数密度()k：s()()n kkV=(7)ss/xyzVL L LN=式中：()n k为第k块切片中的分子数；sV为每块切片的体积；xL、yL、zL分别为模拟盒 x、y、z 方向上的 长度 表 1 不同力场的势参数 Tab.1 Potential parameters for different force fields 参数 OPLS-UA SKS TraPPE-UA NERD SKS-New 31CH/(kJ mol)0.732 6 0.947 9 0.814 9 0.864 8 0.947 9 21CH/(kJ mol)0.493 9 0.390 8 0.382 5 0.380 8 0.390 8 3CH/nm 0.390 5 0.393 0.375 0.391 0.393 2CH/nm 0.390 5 0.393 0.395 0.393 0.393 kb/(kJmol-1nm-2)21.771 5 8.027 2 8.027 2 8.027 2 8.027 2 eqir/nm 0.152 6 0.154 0.154 0.154 0.154 k/(kJmol-1rad-2)519.688 519.688 519.688 519.688 519.688 eq/i()112 114 114 114 114 Vn,An/(kJmol-1)V15.906 V2-1.134 V313.165 V15.906 V2-1.134 V313.165 V15.906 V2-1.134 V313.165 V15.906 V2-1.134 V313.165 A120.204 A20.679 A3-26.041 A4-1.360 A5-2.104 第k块切片的温度()T k：()21B1()3()n kiiiT kmvn k k=(8)式中：Bk为玻尔兹曼常数；im、iv分别表示第i个粒子的质量和速度 第k块切片的压力()p k、法向应力N()pk、切向应力T()pk以及局部表面张力()k12：Bs1()()()=ijijijr u rp kkk TV(9)2NBs1()()()ijijijijxu rpkkk TrV=(10)Ts1()()2=Bpkkk TV 22()2ijijijijijyzu rr+(11)NT()()()=kpkpkx(12)式中：ijx、ijy和ijz分别为粒子i和j的总距离ijr在x、y、z方向上的分量，x为每块切片的厚度 纯液相模型中，时间足够长的情况下，粒子的运动可以近似看作随机游走，均方位移(MSD)与时间成线性关系，根据爱因斯坦公式，自扩散系数D由均方位移和时间的关系计算所得：2()(0)lim6=tr trDt(13)式中()r t为 t 时刻分子的位置，对所有分子取平均值.1.3 模拟设置 本文比较了 5 种力场对从低温到接近临界温度(正庚烷临界温度为 540K、正十二烷临界温度为658K)下气液平衡模拟纯组分正构烷烃饱和物性和纯液相烷烃的自扩散特性的影响以及研究了不同摩燃烧科学与技术 第 29 卷 第 4 期 454 尔比(正庚烷：正十二烷)的双组分烷烃在正庚烷临界温度以下的饱和物性 系统的初始构型如图 1 所示，初始构型气液气构型用于饱和物性的计算；纯液相构型用于自扩散系数的计算 因考虑实际计算效率和对结果的影响，不同的模拟过程选择不同的适合的分子数以及模拟盒尺寸，详细模拟细节在表2中总结.（a）计算饱和物性的气液气构型 （b）计算自扩散系数的纯液相构型 图 1 初始构型 Fig.1 Initial configuration 表 2 不同模拟过程的初始设置 Tab.2 Initial settings for different simulation processes 模型 模拟盒尺寸/nm3 分子数 温度范围/K正庚烷 3088 2 048 293533 正十二烷 407.687.68 1 536 450650 液相正庚烷 4.78353 500 273368 液相正十二烷 5.835.933 576 308364 11 混合物 408.458.45 2 100 413533 12 混合物 408.458.45 2 100 413533 21 混合物 408.458.45 2 100 413533 气液气构型中模拟盒子为长方体，初始时刻液相分子位于盒子中心，两端为气相区 为消除模拟体系的规模限制带来的边界效应，在模拟盒 3 个方向上均采用周期性边界条件，整个模拟过程均使用正则(NVT)系综对系统温度进行恒定采用 velocity-verlet 算法对原子运动方程进行积分，原子初始速度按高斯分布随机分配 取时间步长为 210-6ns，模拟总时长为 3ns，其中前 2ns 为弛豫时间，后 1ns 为统计时间 当整个系统的温度、压力、总能和势能 4个指标仅在固定值附近波动时，则认定达到平衡状态 纯液相模型为立方模拟盒，选择合适的模拟盒尺寸，使初始构型的液体密度在不同温度下接近真实值，本次模拟中初始密度取自 NIST 数据库13 先用NPT 系综运行 1ns 使整个系统压力稳定在常压附近，然后再使用 NVT 系综运行 2ns，同样取最后 1ns作为统计时间 2 结果与讨论 2.1 纯正庚烷的气液平衡 MD模拟 气相与液相的相交处为气液界面，体系平衡状态下即可获得稳定界面 MD 模拟很好地再现了气液两相平衡特征，各力场下气液相密度分布曲线如图 2 所示，在整个体系达到平衡之后，根据这些分布，可以识别出密度较高的液相区域、密度较低的气相区域和气液界面 5 种联合原子力场在低于 473K 以下时都有较好的模拟表现，能够完整地模拟正庚烷从 293473K 的饱和过程，有明显的密度梯度存在，温度的升高导致液相的密度下降，而导致气相的密度增加；而在高于 473K 时，因接近于临界温度(540K)，密度梯度缩小而界面区慢慢消失，5 种力场中只有 OPLS-UA 和改进力场 SKS-New 能完整模拟所有给定温度下的气液平衡过程，在高于 473K 时也能有明显的气液密度梯度，SKS 力场在 533K 时密度梯度消失，而TraPPE-UA 和 NERD 力场在 513K 和 533K 时模拟效果较差，已全部变为气相 2.2 对纯组分正庚烷和正十二烷的热物性验证 2.2.1 饱和物性 由于界面轮廓随时间波动，气液界面的位置在不同的模拟中有所不同 尽管有些波动，气液界面的统计位置依然稳定 平均统计了 5 种力场、两种纯组分在各温度下两边气相区密度和中间液相区密度，来与NIST 数据库13的两相密度进行对比 所得气液共存曲线如图 3 所示 正庚烷气液共存曲线 5 种力场的模拟结果总体上大都与 NIST 值相符，其中 TraPPE-UA 和 NERD 力场在 413K 以下的低温区气液两相区吻合较好，在高于 413K 时误差增大，气相区密度与 NIST 值相比偏大而液相区密度偏小，由于二者在513K 以上之后就全部变为气相，因此无 513K 以上的气液共存密度，SKS 力场在 533K 时亦是如此 OPLS-UA 力场对所有温度下的液相密度模拟结果都略微偏大，而 SKS 及提出的改进力场的模拟结果与 NIST 值最为吻合，但由于新力场在 533K 时还能保持良好的气液两相平衡，说明改进力场对临界温度的预测效果最好 因为选取精确的势函数是表征模拟准确性的前提，力场描述的粒子间吸引力越小，在给定温度下从中心液核弛豫到气相区的分子越多，使得液相区密度偏小而气相区密度偏大，从模拟结果来看，OPLS-UA 力场表述的正构烷烃分子吸引力偏大，TraPPE-UA 和 NERD 力场则偏小，而 SKS 和改陈飞龙等：基于不同势函数的正构烷烃物理性质分子动力学模拟 燃烧科学与技术 455 进力场相对而言最为准确 （a）OPLS-UA （b）SKS （c）TraPPE-UA （d）NERD （e）SKS-New 图 2 不同力场下正庚烷的密度分布曲线 Fig.2 Density distribution curves of n-heptane under different force fields （a）正庚烷 （b）正十二烷 图 3 不同力场下的气液共存曲线 Fig.3 Vapor-liquid coexistence curves under different force fields 不同力场下正十二烷从 450650K 的模拟结果与正庚烷大致趋势相同，同样是 SKS 和改进力场的模拟结果与 NIST 值较为符合，同样 SKS、TraPPE-UA 和 NERD 这 3 种力场在温度升高时气相密度会偏大、液相密度会偏小，而 OPLS-UA 力场则相反.从气液共存密度的模拟结果来看，改进力场 SKS-燃烧科学与技术 第 29 卷 第 4 期 456 New 对正十二烷的描述优于其他 4 种力场，改进力场在每个温度下的模拟结果能相对比较准确地模拟出气液相密度，且在 650K 即接近临界温度(658K)时也能有明显气液界面，在 5 种力场中甚至比 SKS力场对高温度下的描述更为准确，预测临界点的效果更好 这可能是由于改进力场势能参数所用的二面角扭转势能参数比原 SKS 力场更为准确，同时在链更长的烷烃里所包含的二面角更多，因此对二面角扭曲势能更为准确的描述会使得模拟结果更为精确 燃料的饱和蒸气压对于发动机内混合气形成以及后续燃烧过程有着重要影响，因此为验证此改进力场参数是否能准确表述烷烃其他物性，用气液气模型计算了烷烃的不同温度下的饱和蒸气压 将气液平衡状态下两端饱和气相区的压力统计平均值即可得到饱和蒸气压，MD 模拟正庚烷和正十二烷的计算结果与实验值14-15和 NIST 数据库13的对比如图 4 所示 两种纯组分不同力场下的饱和蒸气压都随温度升高而增大，模拟结果与前面两相密度的效果大致相同，5 种力场中还是 SKS 和改进力场 SKS-New 的模拟结果相对较为准确，其他 3 种力场与实验值或数据库数据偏差较大，尤其在对正十二烷的模拟中更加明显 SKS 和改进力场模拟的饱和蒸气压都随温度升高而误差增大，而对于正十二烷改进力场的饱和蒸气压在高温区也能较好地吻合且比 SKS 力场表现更好，这也印证了上面对于改进力场势能参数更适用于 （a）正庚烷 （b）正十二烷 图 4 不同力场下的饱和蒸气压 Fig.4 Saturated vapor pressure under different force fields 表述中长链烷烃的推断 2.2.2 纯液相的自扩散系数 除了对气液平衡体系下正构烷烃进行不同力场的 MD 模拟验证外，还通过建立纯组分纯液相的立方初始构型来模拟计算液相的自扩散系数，用不同力场模拟的正庚烷和正十二烷结果与实验值16-17对比分别如图 5(a)和(b)所示 自扩散系数的总体趋势随着温度升高而增大，与前面气液模型计算的两相密度和饱和蒸气压效果不同，在种力场中最为相符的是OPLS-UA 力场，其次为 SKS-New 改进力场和 SKS力场，而另外两种力场还是与实验值相比偏大 结合上一小节的研究来看，这可能是因为 OPLS-UA 力场描述的粒子间吸引力在 5 种力场中最大，因此模拟过程中分子在自由运动中的均方位移最小，通过均方位移计算的自扩散系数也就最小，而其他力场对应的描述的粒子间吸引力越小，则计算的自扩散系数越大 （a）正庚烷 （b）正十二烷 图 5 不同力场下的液相自扩散系数 Fig.5Self-diffusion coefficient of liquid under different force fields 2.3 正庚烷/正十二烷二元组分的饱和物性 鉴于上小节改进力场对两种正构烷烃良好的模拟表现，用改进力场模拟了不同摩尔比例为 12、11、21 的正庚烷与正十二烷混合物的饱和过程，为了探究正构烷烃混合物与单组分的饱和密度以及饱和蒸气压之间的联系，取温度范围 413533K(正庚烷临界点以下)对于混合物，采用将纯组分模拟值按混合物摩尔比例加权平均的方法计算饱和物性陈飞龙等：基于不同势函数的正构烷烃物理性质分子动力学模拟 燃烧科学与技术 457 值来与混合物实际模拟值simY作对比，探究加权平均的方法是否能符合实际混合物的饱和性质 混合物物性的加权值mixY计算公式如下：mix1=niiiYx y(14)式中：ix为纯组分的摩尔分数，iy为纯组分物性的模拟值 结果如图 6 所示，混合物物性并不是单组分物性的简单加权平均，由于正十二烷分子对正庚烷分子的吸附作用，本该在接近临界温度时全部变为气相的正庚烷分子比预计更少地从液核弛豫到气相区，因此液相密度的加权值始终比实际模拟值偏小，相应地，气相密度和饱和蒸气压始终偏大，且总体趋势似乎是随着温度增大而偏差越多 （a）液相密度 （b）气相密度 （c）饱和蒸气压 图 6 混合物物性的加权值与实际模拟值 Fig.6 Weighted and actual simulated values of mixture properties 为了更直观地比较混合物物性加权值与实际模拟值的差异，统计了两者的相对误差随不同温度和混合比的变化关系 相对误差的计算公式如下：mixsimsim100%YYY=(15)图 7 展示了混合物物性加权平均值与实际模拟值的相对误差随温度、混合比(用混合物中正庚烷的摩尔分数1x来表示)的变化关系，可以看出混合物的饱和物性不是随温度或者混合比的单一变化，而是受到两者的双重影响，在温度较低的时候，物性的相对误差随正庚烷摩尔分数的减小而增大，而在温度较高的时候，相对误差又是随着正庚烷摩尔分数的减小而减小，而整体随温度的变化趋势是先增大、后减小、（a）液相密度 （b）气相密度 （c）饱和蒸气压 图 7 相对误差随温度、混合比的变化关系 Fig.7Variation of relative error with temperature and mixing ratio 燃烧科学与技术 第 29 卷 第 4 期 458 再增大 结合图 6 可以理解为，在温度较低时，混合比占影响因素的主导地位，而正庚烷组分越多，可以平衡双组分之间的性质差异，因此更符合按加权平均的估算；在温度较高时，温度占影响因素的主导地位，由温度增大带来的对相对误差的影响大于混合比，而在接近正庚烷临界温度时，混合物的性质更接近于正十二烷，因此正庚烷组分越少相对误差越小；在两个因素的共同影响作用下，相对误差在温度和混合比的中间区域会出现一个小的峰值，但随着温度升高，由于正庚烷弛豫到为气相区的能力削弱，气相密度的加权值偏大得越来越多，因此相对误差整体会随着温度不断增大 结果表明，采用将单组分正构烷烃饱和物性按混合物摩尔比例加权平均的方法去估算混合物的饱和物性时，在温度较低、正庚烷比重较高的条件下，加权平均得到的饱和物性值较为准确；在温度较高、正庚烷比重较低的条件下，加权平均得到的饱和物性值较为准确 3 结 论 采用 MD 模拟比较了不同势函数对不同温度下纯组分正庚烷和正十二烷燃料物理性质的影响，并提出了一种新的修正力场并对其进行了烷烃饱和两相密度、饱和蒸气压和自扩散系数等热物性的验证 还将改进力场应用于二元组分烷烃的模拟中，主要研究了不同摩尔比(12、11、21)的双组分烷烃在正庚烷临界温度下的饱和过程 本文的研究结果如下.(1)建立了单组分烷烃气液平衡的 MD 模型 比较了不同势函数对于预测两种正构烷烃两相饱和物性的精度差异，分析了温度对于饱和状态下两相密度以及饱和蒸气压的影响 随温度升高，液相密度减小，气相密度与饱和蒸气压增大，5 种力场皆能良好地模拟烷烃在近临界温度下的饱和过程 (2)提出了一种改进的适用于正构烷烃的修正力场，比较了与其他 4 种常用力场对预测饱和两相密度、饱和蒸气压以及自扩散系数的精度 其中 SKS力场对于正庚烷的饱和性质模拟最为准确，但对于临界点的预测不如改进力场 SKS-New；改进力场对于正十二烷的饱和性质模拟相较于 SKS 力场要更优，且对两种正构烷烃的临界温度的预测更准确；OPLS-UA 力场对于自扩散系数的模拟最为准确，其次为改进力场 综合来看，改进力场对正构烷烃的粒子间作用力描述最优且更适用于中长链烷烃 (3)研究了双组分烷烃饱和过程 探究了不同混合比和温度条件下通过将纯组分饱和物性按不同摩尔比加权平均的方法估算混合物饱和物性的准确度 在估算饱和物性时，选择温度较低、正庚烷比重较高或者是温度较高、正庚烷比重较低的工况下加权平均更为准确 参考文献：1 Jorgensen W L，Tirado-Rives J.The OPLSoptimized potentials for liquid simulations potential functions for proteins，energy minimizations for crystals of cyclic peptides and crambinJ.Journal of the American Chemical Society，1988，110(6)：1657-1666.2 Jorgensen W L，Maxwell D S，Tirado-Rives J.Devel-opment and testing of the OPLS all-atom force field on conformational energetics and properties of organic liq-uidsJ.Journal of the American Chemical Society，1996，118(45)：11225-11236.3 Jorgensen W L.Optimized intermolecular potential func-tions for liquid alcohols J.The Journal of Physical Chemistry，1986，90(7)：1276-1284.4 Jorgensen W L，Madura J D，Swenson C J.Optimized intermolecular potential functions for liquid hydrocar-bonsJ.Journal of the American Chemical Society，1984，106(22)：6638-6646.5 Siepmann J I，Karaborni S，Smit B.Vapor-liquid equilibria of model alkanesJ.Journal of the American Chemical Society，1993，115(14)：6454-6455.6 Siepmann J I，Karaborni S，Smit B.Simulating the critical behaviour of complex fluidsJ.Nature，1993，365：330-332.7 Martin M G，Siepmann J I.Transferable potentials for phase equilibria(1)：United-atom description of n-alkanesJ.The Journal of Physical Chemistry B，1998，102(14)：2569-2577.8 Nath S K，Escobedo F A，de Pablo J J.On the simula-tion of vapor-liquid equilibria for alkanesJ.The Jour-nal of Chemical Physics，1998，108(23)：9905-9911.9 邓 磊.亚/超临界环境下气液界面特性的分子动力学模拟D.大连：大连理工大学能源与动力学院，2015.Deng Lei.Molecular Dynamics Simulation of Liquid-Gas Interface in Sub/Supercritical SurroundingsD.Da-陈飞龙等：基于不同势函数的正构烷烃物理性质分子动力学模拟 燃烧科学与技术 459 lian：School of Energy and Power Engineering，Dalian University of Technology，2015(in Chinese).10 Kalyanasundaram V，Spearot D E，Malshe A P.Mo-lecular dynamics simulation of nanoconfinement induced organization of n-decaneJ.Langmuir：the ACS Jour-nal of Surfaces and Colloids，2009，25(13)：7553-7560.11 Mo G，Qiao L.A molecular dynamics investigation of n-alkanes vaporizing into nitrogen：Transition from sub-critical to supercriticalJ.Combustion and Flame，2017，176：60-71.12 Lee S-H.Pressure analyses at the planar surface of liquid-vapor argon by a test-area molecular dynamics simula-tionJ.Bulletin of the Korean Chemical Society，2012，33(9)：3039-3042.13 NIST.NIST Chemistry WebBookEB/OL.https:/webbook.nist.gov/chemistry/，2022-03-16.14 Weber L A.Vapor pressure of heptane from the triple point to the critical pointJ.Journal of Chemical&En-gineering Data，2000，45(2)：173-176.15 Morgan D L，Kobayashi R.Direct vapor pressure meas-urements of ten n-alkanes in the C10-C28 rangeJ.Fluid Phase Equilibria，1994，97：211-242.16 Fishman E.Self-diffusion in liquid n-pentane and n-heptaneJ.The Journal of Physical Chemistry，1955，59(5)：469-472.17 Harris K R，Ganbold B，Price W S.Viscous calibration liquids for self-diffusion measurementsJ.Journal of Chemical&Engineering Data，2015，60(12)：3506-3517.(责任编辑：隋韶颖)