CN113345530A - 一种基于分子动力学的二元体系相互扩散系数模拟方法 - Google Patents

一种基于分子动力学的二元体系相互扩散系数模拟方法 Download PDF

Info

Publication number
CN113345530A
CN113345530A CN202110749565.XA CN202110749565A CN113345530A CN 113345530 A CN113345530 A CN 113345530A CN 202110749565 A CN202110749565 A CN 202110749565A CN 113345530 A CN113345530 A CN 113345530A
Authority
CN
China
Prior art keywords
binary system
model
molecular dynamics
interdiffusion
under
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202110749565.XA
Other languages
English (en)
Inventor
陈伟
许晨滨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qingdao University of Science and Technology
Original Assignee
Qingdao University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Qingdao University of Science and Technology filed Critical Qingdao University of Science and Technology
Priority to CN202110749565.XA priority Critical patent/CN113345530A/zh
Publication of CN113345530A publication Critical patent/CN113345530A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/70Machine learning, data mining or chemometrics

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computing Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Health & Medical Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Chemical & Material Sciences (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Medical Informatics (AREA)
  • Evolutionary Computation (AREA)
  • Data Mining & Analysis (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于分子动力学的二元体系相互扩散系数模拟方法,具体包括以下步骤:构建两个浓度具有差异的二元体系模型,在系宗NPT的作用下,达到目标温度与压力下的能量和密度平衡;将平衡后的两个二元体系模型在系宗NVT的作用下,进行目标温度与压力下的相互扩散模拟,一定时间内对单个二元体系模型中的分子数目进行统计;计算互相扩散系数。采用上述的一种基于分子动力学的二元体系相互扩散系数模拟方法,通过构建二元体系模型并对其进行单向互扩散模拟,统计中特定区域中分子数目变化情况,计算二元体系的相互扩散系数,避免外界干扰,提高计算精度,同时提高测试效率与应用范围并极大地降低成本。

Description

一种基于分子动力学的二元体系相互扩散系数模拟方法
技术领域
本发明涉及流体热物理性质计算技术领域,尤其涉及一种基于分子动力学的二元体系相互扩散系数模拟方法。
背景技术
扩散是势差作用下的分子热运动现象,广泛地存在于实际生活中。相互扩散是异类原子由浓度等差异引起的相互迁移、相互渗透的过程,表征这类过程传递行为快慢的参数称为相互扩散系数,是科研、工程领域重要的热物理性质。
目前,测量二元体系相互扩散系数通常采用实验的方法,例如隔膜电池法、同位素法、核磁共振法以及全息干涉法等。但采用实验的方法有着实验设备昂贵、条件要求高、结果处理过程复杂等缺陷。
发明内容
本发明的目的是提供一种基于分子动力学的二元体系相互扩散系数模拟方法,具有操作简易、花费较少和更为直观的优点,同时减少了实验测试过程中的不确定因素,可以弥补实验测试的诸多不足,大大提高了测试效率与准确性。
为实现上述目的,本发明提供了一种基于分子动力学的二元体系相互扩散系数模拟方法,具体包括以下步骤:
步骤S1:构建两个浓度具有差异的二元体系模型,在系宗NPT的作用下,达到目标温度与压力下的能量和密度平衡;
步骤S2:将平衡后的两个二元体系模型在系宗NVT的作用下,进行目标温度与压力下的相互扩散模拟,一定时间内对单个二元体系模型中的分子数目进行统计;
步骤S3:计算互相扩散系数,计算公式如下,
Figure BDA0003145529980000021
其中C为混合体系的浓度,C1和C2分别为两个二元体系模型的浓度,Z代表的是任一二元体系模型中扩散方向的距离,τ代表的是相互扩散的时间,D12代表的是相互扩散系数。
进一步的,在步骤S1中,在系宗NPT的作用下,达到目标温度与压力下的能量和密度平衡的具体步骤如下:
用分子动力学软件Packmol及Moltemplate构建二元体系模型,力场形式采用OPLS-AA:共价键伸缩势和键角弯曲势采用harmonic势函数,二面角扭曲势只包括Fourier展开式的前三项,van der Waals相互作用采用Lennard-Jones 12-6势函数,静电相互作用采用库伦势函数;
用分子动力学软件Lammps进行后续过程模拟,为避免周期性边界条件对后续二元体系的单项扩散产生影响,在模型的下端添加一层挡板,而后采用系宗NPT使模型在目标温度及压力下进行能量、密度的平衡,压力控制方向设定为Z方向。
进一步的,在步骤S2中,对单个二元体系模型中的分子数目进行统计的具体步骤如下:
将两个平衡后的二元体系模型的X、Y截面拼接在一起,使两个二元体系模型的原子不重叠,然后在系宗NVT的作用下,使二元体系模型在目标温度下进行分子动力学模拟,在一定时间内,输出原子轨迹,对特定分子的原子进行分类和统计,确定单个二元体系模型中特定分子数目随时间变化的规律。
进一步的,在步骤S3中,互相扩散系数的计算公式经过最小二乘法处理后得到的计算公式如下:
Figure BDA0003145529980000031
其中,C(z,τ)的计算方法为较高浓度或较低浓度体系模型中随时间变化的特定分子的个数与该模型区域体积的比值;C1的计算方法为单独驰豫后较高浓度体系模型中特定分子的个数与该模型区域体积的比值;C2的计算方法为单独驰豫后较低浓度体系模型中特定分子的个数与该模型区域体积的比值。
因此,本发明采用上述一种基于分子动力学的二元体系相互扩散系数模拟方法,具有以下有益效果:
(1)、本发明采用分子动力学模拟,通过Packmol及Moltemplate建立二元体系模型,通过大规模原子并行模拟器Lammps输出原子运动轨迹,达到快速计算的目的,通过统计一定区域内特定分子的数目,继而根据所提出的公式计算得到目标温度、压力、浓度条件下二元体系的相互扩散系数,时间短、计算方便、效率高,避免了繁琐的实验过程,且规避了外界干扰,结果准确;
(2)计算机可以连续工作,通过增加计算核数可以大大提高计算效率,可以预测不同温度、压力等其它情况下二元体系的相互扩散系数,为后续科学研究及工程实践提供理论依据。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1是本发明中所述方法的流程图;
图2是本发明中通过Packmol及Moltemplate建立的二元体系模型示意图;
图3是本发明中两个浓度具有差异的二元体系模型单项扩散示意图;
图4是本发明中高浓度体系模型中特定分子数量随时间变化示意图;
图5是本发明中计算出的以离子液体质量分数为20%的[emim]BF4/CH3OH二元体系为例在常温常压下的互扩散系数。
图中:阳离子1、阴离子2、甲醇分子3。
具体实施方式
实施例
如图1所示,一种基于分子动力学的二元体系相互扩散系数模拟方法,具体包括以下步骤:
步骤S1:构建两个浓度具有差异的二元体系模型,两个浓度具有差异很微小。在系宗NPT的作用下,达到目标温度与压力下的能量和密度平衡。
在系宗NPT的作用下,达到目标温度与压力下的能量和密度平衡的具体步骤如下:
用分子动力学软件Packmol及Moltemplate构建二元体系模型,力场形式采用OPLS-AA:共价键伸缩势和键角弯曲势采用harmonic势函数,二面角扭曲势只包括Fourier展开式的前三项,van der Waals相互作用采用Lennard-Jones 12-6势函数,静电相互作用采用库伦势函数,具体势函数作用形式如下:
E(rN)=Ebonds+Eangles+Edihedrals+Enonbonded
Figure BDA0003145529980000041
Figure BDA0003145529980000042
Figure BDA0003145529980000051
Figure BDA0003145529980000052
用分子动力学软件Lammps进行后续过程模拟,为避免周期性边界条件对后续二元体系的单项扩散产生影响,在模型的下端添加一层挡板,而后采用系宗NPT使模型在目标温度及压力下进行能量、密度的平衡,压力控制方向设定为Z方向。
步骤S2:将平衡后的两个二元体系模型在系宗NVT的作用下,进行目标温度与压力下的相互扩散模拟,一定时间内对单个二元体系模型中的分子数目进行统计。
对单个二元体系模型中的分子数目进行统计的具体步骤如下:
将两个平衡后的二元体系模型的X、Y截面拼接在一起,使两个二元体系模型的原子不重叠,然后在系宗NVT的作用下,使二元体系模型在目标温度下进行分子动力学模拟,在一定时间内,输出原子轨迹,对特定分子的原子进行分类和统计,确定单个二元体系模型中特定分子数目随时间变化的规律。
步骤S3:计算互相扩散系数,计算公式如下,
Figure BDA0003145529980000053
其中C为混合体系的浓度,C1和C2分别为两个二元体系模型的浓度,Z代表的是任一二元体系模型中扩散方向的距离,τ代表的是相互扩散的时间,D12代表的是相互扩散系数。
互相扩散系数的计算公式经过最小二乘法处理后得到的计算公式如下:
Figure BDA0003145529980000061
其中,C(z,τ)的计算方法为较高浓度或较低浓度体系模型中随时间变化的特定分子的个数与该模型区域体积的比值;C1的计算方法为单独驰豫后较高浓度体系模型中特定分子的个数与该模型区域体积的比值;C2的计算方法为单独驰豫后较低浓度体系模型中特定分子的个数与该模型区域体积的比值。
根据上述方法计算离子液体质量分数为19.5%的[emim]BF4/CH3OH二元体系在常温常压下的相互扩散系数,
如图2所示,左侧为盒子体系大小97A*97A*104A的离子液体质量分数为16.1%的[emim]BF4/CH3OH体系模型,含有的总原子数为92994个;右侧为盒子体系大小97A*97A*104A离子液体质量分数为23.6%的[emim]BF4/CH3OH体系模型含有的总原子数为99714个。力场形式采用OPLS-AA:共价键伸缩势和键角弯曲势采用harmonic势函数,二面角扭曲势只包括Fourier展开式的前三项,van der Waals相互作用采用Lennard-Jones 12-6势函数,静电相互作用采用库伦势函数。保持挡板不受力的情况,设置控压方向为Z方向,保证X、Y截面大小不变,时间步长为1fs,在系宗NPT的作用下进行常温、常压驰豫,模拟时长为1ns。
如图3所示,将驰豫后的模拟盒子以无挡板的一面为拼接面,形成可单向扩散的高、低浓度的[emim]BF4/CH3OH体系模型,时间步长为1fs,在系宗NVT的作用下进行常温、常压的模拟,运行70万步,时长为700ps,每500步记录一组原子轨迹数据,共计1400组数据。通过对原子轨迹数据进行处理,统计原较高浓度的模拟盒子中特定分子(实施例中统计的为阳离子)数目,将特定分子汇总到一个图上,并把图中的数量点与模拟运行时间一一对应,如图4所示。
利用如下公式进行扩散系数的计算:
Figure BDA0003145529980000071
用数学计算软件对上公式进行最小二乘法处理,可有计算公式:
Figure BDA0003145529980000072
其中C代表的是混合体系的浓度,Z代表的是任一二元体系模型中扩散方向的距离,τ代表的是相互扩散的时间,D12代表的是相互扩散系数。C(z,τ)的计算方法为较高浓度或较低浓度体系模型中随时间变化的某特定分子的个数与该模型区域体积的比值。C1的计算方法为单独驰豫后较高浓度体系模型中某特定分子的个数与该模型区域体积的比值。C2的计算方法为单独驰豫后较低浓度体系模型中某特定分子的个数与该模型区域体积的比值。在本实施例中,C(z,τ)的计算结果如图5右侧曲线所示,C1的计算结果为0.000584个/A3,C2的计算结果为0.000365个/A3。将数值代入上述公式中,可以通过计算得出离子液体质量分数为19.5%的[emim]BF4/CH3OH二元体系在常温常压下的相互扩散系数,如图5所示。
模拟结果说明:本发明可以运用到分子动力学模拟实例计算的过程中,可以真实模拟计算不同温度、压力条件下某浓度二元体系的相互扩散系数。
因此,本发明采用上述一种基于分子动力学的二元体系相互扩散系数模拟方法,具有操作简易、花费较少和更为直观的优点,同时减少了实验测试过程中的不确定因素,可以弥补实验测试的诸多不足,大大提高了测试效率与准确性。
最后应说明的是:以上实施例仅用以说明本发明的技术方案而非对其进行限制,尽管参照较佳实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对本发明的技术方案进行修改或者等同替换,而这些修改或者等同替换亦不能使修改后的技术方案脱离本发明技术方案的精神和范围。

Claims (4)

1.一种基于分子动力学的二元体系相互扩散系数模拟方法,其特征在于,具体包括以下步骤:
步骤S1:构建两个浓度具有差异的二元体系模型,在系宗NPT的作用下,达到目标温度与压力下的能量和密度平衡;
步骤S2:将平衡后的两个二元体系模型在系宗NVT的作用下,进行目标温度与压力下的相互扩散模拟,一定时间内对单个二元体系模型中的分子数目进行统计;
步骤S3:计算互相扩散系数,计算公式如下,
Figure FDA0003145529970000011
其中C为混合体系的浓度,C1和C2分别为两个二元体系模型的浓度,Z代表的是任一二元体系模型中扩散方向的距离,τ代表的是相互扩散的时间,D12代表的是相互扩散系数。
2.根据权利要求1所述的一种基于分子动力学的二元体系相互扩散系数模拟方法,其特征在于:在步骤S1中,在系宗NPT的作用下,达到目标温度与压力下的能量和密度平衡的具体步骤如下:
用分子动力学软件Packmol及Moltemplate构建二元体系模型,力场形式采用OPLS-AA:共价键伸缩势和键角弯曲势采用harmonic势函数,二面角扭曲势只包括Fourier展开式的前三项,van der Waals相互作用采用Lennard-Jones 12-6势函数,静电相互作用采用库伦势函数;
用分子动力学软件Lammps进行后续过程模拟,为避免周期性边界条件对后续二元体系的单项扩散产生影响,在模型的下端添加一层挡板,而后采用系宗NPT使模型在目标温度及压力下进行能量、密度的平衡,压力控制方向设定为Z方向。
3.根据权利要求1所述的一种基于分子动力学的二元体系相互扩散系数模拟方法,其特征在于:在步骤S2中,对单个二元体系模型中的分子数目进行统计的具体步骤如下:
将两个平衡后的二元体系模型的X、Y截面拼接在一起,使两个二元体系模型的原子不重叠,然后在系宗NVT的作用下,使二元体系模型在目标温度下进行分子动力学模拟,在一定时间内,输出原子轨迹,对特定分子的原子进行分类和统计,确定单个二元体系模型中特定分子数目随时间变化的规律。
4.根据权利要求1所述的一种基于分子动力学的二元体系相互扩散系数模拟方法,其特征在于:在步骤S3中,互相扩散系数的计算公式经过最小二乘法处理后得到的计算公式如下:
Figure FDA0003145529970000021
其中,C(z,τ)的计算方法为较高浓度或较低浓度体系模型中随时间变化的特定分子的个数与该模型区域体积的比值;C1的计算方法为单独驰豫后较高浓度体系模型中特定分子的个数与该模型区域体积的比值;C2的计算方法为单独驰豫后较低浓度体系模型中特定分子的个数与该模型区域体积的比值。
CN202110749565.XA 2021-07-02 2021-07-02 一种基于分子动力学的二元体系相互扩散系数模拟方法 Pending CN113345530A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110749565.XA CN113345530A (zh) 2021-07-02 2021-07-02 一种基于分子动力学的二元体系相互扩散系数模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110749565.XA CN113345530A (zh) 2021-07-02 2021-07-02 一种基于分子动力学的二元体系相互扩散系数模拟方法

Publications (1)

Publication Number Publication Date
CN113345530A true CN113345530A (zh) 2021-09-03

Family

ID=77482310

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110749565.XA Pending CN113345530A (zh) 2021-07-02 2021-07-02 一种基于分子动力学的二元体系相互扩散系数模拟方法

Country Status (1)

Country Link
CN (1) CN113345530A (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104063561A (zh) * 2014-07-11 2014-09-24 哈尔滨工业大学 一种利用计算机模拟纳米材料与环境本底纳米污染物在水环境中相互作用的方法
CN109885895A (zh) * 2019-01-25 2019-06-14 南京航空航天大学 一种基于分子动力学的材料表面结冰形核过程的监测方法
CN110021380A (zh) * 2019-04-25 2019-07-16 济南大学 一种基于分子动力学模拟探究玻璃体系中各原子扩散性质的方法
CN110097927A (zh) * 2019-05-10 2019-08-06 青岛理工大学 基于分子动力学测试电场作用下离子扩散系数的方法
JP2019186133A (ja) * 2018-04-13 2019-10-24 パナソニックIpマネジメント株式会社 燃料電池の発電性能を予測するシミュレーション方法
JP2019186200A (ja) * 2018-04-12 2019-10-24 パナソニックIpマネジメント株式会社 パラメータ決定方法および細孔内のガスまたはイオンの輸送性を求めるシミュレーション方法
CN110767267A (zh) * 2019-09-30 2020-02-07 华中科技大学 一种基于Python对ReaxFF力场计算结果数据处理的方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104063561A (zh) * 2014-07-11 2014-09-24 哈尔滨工业大学 一种利用计算机模拟纳米材料与环境本底纳米污染物在水环境中相互作用的方法
JP2019186200A (ja) * 2018-04-12 2019-10-24 パナソニックIpマネジメント株式会社 パラメータ決定方法および細孔内のガスまたはイオンの輸送性を求めるシミュレーション方法
JP2019186133A (ja) * 2018-04-13 2019-10-24 パナソニックIpマネジメント株式会社 燃料電池の発電性能を予測するシミュレーション方法
CN109885895A (zh) * 2019-01-25 2019-06-14 南京航空航天大学 一种基于分子动力学的材料表面结冰形核过程的监测方法
CN110021380A (zh) * 2019-04-25 2019-07-16 济南大学 一种基于分子动力学模拟探究玻璃体系中各原子扩散性质的方法
CN110097927A (zh) * 2019-05-10 2019-08-06 青岛理工大学 基于分子动力学测试电场作用下离子扩散系数的方法
CN110767267A (zh) * 2019-09-30 2020-02-07 华中科技大学 一种基于Python对ReaxFF力场计算结果数据处理的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
张颖等: "液相互扩散系数测量充注方法的改进", 《热科学与技术》 *
郑燕升等: "室温离子液体的分子动力学模拟", 《化学进展》 *
颜彩繁: "《固体物理专题实验》", 31 August 2020, 南开大学出版社 *

Similar Documents

Publication Publication Date Title
Bonnet et al. Airfrans: High fidelity computational fluid dynamics dataset for approximating reynolds-averaged navier–stokes solutions
CN109470553B (zh) 应用Laplace变换实现沥青混合料材料参数快速转换的方法
CN106596087A (zh) 一种新型可实现全组件摩擦辨识的装置和方法
CN112185473A (zh) 一种分子动力学仿真测试系统
CN115544675B (zh) 高超声速飞行器防热材料表面催化特性多尺度预测方法
CN105824987A (zh) 一种基于遗传算法的风场特征统计分布模型建立方法
Fu et al. Nonlinear vibration analysis of a rotor system with parallel and angular misalignments under uncertainty via a Legendre collocation approach
CN110648724B (zh) 一种电渗析传质过程的数值模拟方法
CN113345530A (zh) 一种基于分子动力学的二元体系相互扩散系数模拟方法
CN104807606A (zh) 一种matlab-stm32混合动力试验系统及其试验方法
Zhu et al. Compilation of program-loading spectrum for milling of a motorized spindle based on cutting force model
CN106326574A (zh) 一种计及焊接变形和残余应力的极限强度数值计算方法
Kalinitchev New model of multicomponent masstransfer kinetics in bifunctional matrix of nanocomposites and results of simulation
Srinivasan et al. Latency tolerance through parallelization of time in scientific applications
CN112414708A (zh) 一种滚珠丝杠进给系统的轴向动态特性参数识别方法
Tam et al. Characterizing mechanical properties of polymeric material: A bottom-up approach
Fekete et al. Investigation on geometric imperfections of tensile test specimens using optical full-field measurements and digital twin-based simulations
Rao et al. Research on prediction method on RUL of motor of CNC machine based on deep learning
Hu et al. Investigation of the dynamic plastic hardening of metal thin-walled tube under liquid impact forming
te Braake Improving the PX4 software-in-the-loop multirotor UAV simulator accuracy
Abbitt et al. Implementing the digital design process for the development of a centrifugal fan impeller in the undergraduate engineering curriculum
CN117077336A (zh) 一种天然气管网动态气源占比跟踪方法
Yin et al. Evaluation of task specific measurement uncertainty for gear measuring instrument using VGMI
Guo et al. Application of ant colony algorithm to volumetric thermal error modeling and compensation of a CNC machine tool
CN112257312A (zh) 一种透平叶片材料的微细片状颗粒群冲蚀模型参数化建模方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination