CN113345530A - 一种基于分子动力学的二元体系相互扩散系数模拟方法 - Google Patents
一种基于分子动力学的二元体系相互扩散系数模拟方法 Download PDFInfo
- 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
Links
- 229910002056 binary alloy Inorganic materials 0.000 title claims abstract description 66
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000004088 simulation Methods 0.000 title claims abstract description 26
- 238000000329 molecular dynamics simulation Methods 0.000 title claims abstract description 19
- 238000004364 calculation method Methods 0.000 claims abstract description 24
- 230000009471 action Effects 0.000 claims abstract description 17
- 230000008859 change Effects 0.000 claims abstract description 4
- 238000009792 diffusion process Methods 0.000 claims description 20
- 230000008569 process Effects 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 6
- 238000005452 bending Methods 0.000 claims description 4
- 230000009881 electrostatic interaction Effects 0.000 claims description 4
- 230000003993 interaction Effects 0.000 claims description 4
- 230000000737 periodic effect Effects 0.000 claims description 3
- 238000000465 moulding Methods 0.000 claims 1
- 238000012360 testing method Methods 0.000 abstract description 7
- OKKJLVBELUTLKV-UHFFFAOYSA-N Methanol Chemical compound OC OKKJLVBELUTLKV-UHFFFAOYSA-N 0.000 description 9
- NJMWOUFKYKNWDW-UHFFFAOYSA-N 1-ethyl-3-methylimidazolium Chemical compound CCN1C=C[N+](C)=C1 NJMWOUFKYKNWDW-UHFFFAOYSA-N 0.000 description 6
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Natural products C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 6
- 239000002608 ionic liquid Substances 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 150000001768 cations Chemical class 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 150000001450 anions Chemical class 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000005210 holographic interferometry Methods 0.000 description 1
- 150000002500 ions Chemical class 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001225 nuclear magnetic resonance method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/70—Machine 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:计算互相扩散系数,计算公式如下,
其中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中,互相扩散系数的计算公式经过最小二乘法处理后得到的计算公式如下:
其中,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,
用分子动力学软件Lammps进行后续过程模拟,为避免周期性边界条件对后续二元体系的单项扩散产生影响,在模型的下端添加一层挡板,而后采用系宗NPT使模型在目标温度及压力下进行能量、密度的平衡,压力控制方向设定为Z方向。
步骤S2:将平衡后的两个二元体系模型在系宗NVT的作用下,进行目标温度与压力下的相互扩散模拟,一定时间内对单个二元体系模型中的分子数目进行统计。
对单个二元体系模型中的分子数目进行统计的具体步骤如下:
将两个平衡后的二元体系模型的X、Y截面拼接在一起,使两个二元体系模型的原子不重叠,然后在系宗NVT的作用下,使二元体系模型在目标温度下进行分子动力学模拟,在一定时间内,输出原子轨迹,对特定分子的原子进行分类和统计,确定单个二元体系模型中特定分子数目随时间变化的规律。
步骤S3:计算互相扩散系数,计算公式如下,
其中C为混合体系的浓度,C1和C2分别为两个二元体系模型的浓度,Z代表的是任一二元体系模型中扩散方向的距离,τ代表的是相互扩散的时间,D12代表的是相互扩散系数。
互相扩散系数的计算公式经过最小二乘法处理后得到的计算公式如下:
其中,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所示。
利用如下公式进行扩散系数的计算:
用数学计算软件对上公式进行最小二乘法处理,可有计算公式:
其中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)
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的作用下,使二元体系模型在目标温度下进行分子动力学模拟,在一定时间内,输出原子轨迹,对特定分子的原子进行分类和统计,确定单个二元体系模型中特定分子数目随时间变化的规律。
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)
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力场计算结果数据处理的方法 |
-
2021
- 2021-07-02 CN CN202110749565.XA patent/CN113345530A/zh active Pending
Patent Citations (7)
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)
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 |