CN111241655A - 用于分子动力学计算机模拟的神经网络力场计算算法 - Google Patents
用于分子动力学计算机模拟的神经网络力场计算算法 Download PDFInfo
- Publication number
- CN111241655A CN111241655A CN201911182094.8A CN201911182094A CN111241655A CN 111241655 A CN111241655 A CN 111241655A CN 201911182094 A CN201911182094 A CN 201911182094A CN 111241655 A CN111241655 A CN 111241655A
- Authority
- CN
- China
- Prior art keywords
- features
- neural network
- rotationally
- covariant
- invariant
- 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
- 238000013528 artificial neural network Methods 0.000 title claims abstract description 58
- 238000004422 calculation algorithm Methods 0.000 title description 46
- 238000000329 molecular dynamics simulation Methods 0.000 title description 34
- 238000005094 computer simulation Methods 0.000 title description 3
- 238000012617 force field calculation Methods 0.000 title description 2
- 238000000034 method Methods 0.000 claims abstract description 40
- 239000013598 vector Substances 0.000 claims abstract description 37
- 238000000205 computational method Methods 0.000 claims abstract description 4
- 230000006870 function Effects 0.000 claims description 30
- 238000012549 training Methods 0.000 claims description 23
- 238000004364 calculation method Methods 0.000 claims description 16
- 239000000463 material Substances 0.000 claims description 16
- 238000004088 simulation Methods 0.000 claims description 13
- 239000000446 fuel Substances 0.000 claims description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 238000000576 coating method Methods 0.000 claims description 4
- 238000010612 desalination reaction Methods 0.000 claims description 3
- 238000006555 catalytic reaction Methods 0.000 claims description 2
- 230000004931 aggregating effect Effects 0.000 claims 1
- 125000004429 atom Chemical group 0.000 description 59
- 229920003171 Poly (ethylene oxide) Polymers 0.000 description 13
- 238000010586 diagram Methods 0.000 description 10
- 125000004432 carbon atom Chemical group C* 0.000 description 9
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 8
- 238000006243 chemical reaction Methods 0.000 description 7
- 125000004430 oxygen atom Chemical group O* 0.000 description 7
- 238000003775 Density Functional Theory Methods 0.000 description 6
- 239000000126 substance Substances 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 229920000642 polymer Polymers 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 230000006399 behavior Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000005610 quantum mechanics Effects 0.000 description 3
- 238000000367 ab initio method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 2
- 238000013434 data augmentation Methods 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 150000002500 ions Chemical class 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 239000008204 material by function Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 229920000052 poly(p-xylylene) Polymers 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 241000969130 Atthis Species 0.000 description 1
- 239000005711 Benzoic acid Substances 0.000 description 1
- 241001432959 Chernes Species 0.000 description 1
- 241000321453 Paranthias colonus Species 0.000 description 1
- 239000011362 coarse particle Substances 0.000 description 1
- 238000013034 coating degradation Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 239000003792 electrolyte Substances 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- UQSQSQZYBQSBJZ-UHFFFAOYSA-N fluorosulfonic acid Chemical compound OS(F)(=O)=O UQSQSQZYBQSBJZ-UHFFFAOYSA-N 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000037427 ion transport Effects 0.000 description 1
- 229910001386 lithium phosphate Inorganic materials 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 125000004433 nitrogen atom Chemical group N* 0.000 description 1
- 230000002688 persistence Effects 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 239000005365 phosphate glass Substances 0.000 description 1
- 238000004869 quantum mechanical method Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000013526 transfer learning Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- TWQULNDIKKJZPH-UHFFFAOYSA-K trilithium;phosphate Chemical compound [Li+].[Li+].[Li+].[O-]P([O-])([O-])=O TWQULNDIKKJZPH-UHFFFAOYSA-K 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N10/00—Quantum computing, i.e. information processing based on quantum-mechanical phenomena
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- 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
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computing Systems (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Software Systems (AREA)
- Data Mining & Analysis (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Mathematical Physics (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Molecular Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Databases & Information Systems (AREA)
- Medical Informatics (AREA)
- Chemical & Material Sciences (AREA)
- Crystallography & Structural Chemistry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种用于使用神经网络力场(NNFF)来模拟多元素系统内的元素的运动的计算方法。该方法包括:接收多元素系统的局部环境的多个旋转不变特征和多个旋转协变特征的组合;和基于多个旋转不变特征和多个旋转协变特征的组合以及NNFF来预测用于多元素系统内的每个元素的力向量,以获得多元素系统内的元素的模拟运动。
Description
技术领域
本发明大体上涉及用于诸如电化学和水过滤设备之类的材料系统中的分子动力学计算机模拟中的原子力的直接预测的神经网络力场(NNFF)计算算法。
背景技术
分子动力学是用于模拟真实操作压力和温度条件下的材料系统中的原子的运动的计算材料科学方法。存在用于计算在原子运动模拟中使用的基础原子力的方法。一种方法是从头算量子力学方法。这种方法非常准确,但是由于应用该方法所需的大量的计算资源,所以也非常昂贵。虽然存在消耗较少计算资源的其他方法,但是这些其他方法并没有提供同样高的准确度。
发明内容
在第一实施例中,公开了一种用于使用神经网络力场(NNFF)来模拟多元素系统内的元素的运动的计算方法。该方法包括接收多元素系统的局部环境的多个旋转不变特征和多个旋转协变特征的组合。该方法还包括基于多个旋转不变特征、多个旋转协变特征的组合以及NNFF来预测用于多元素系统内的每个元素的力向量,以获得多元素系统内的元素的模拟运动。
在第二实施例中,公开了有形地体现用于软件程序的计算机可读指令的非暂时性计算机可读介质。该软件程序可由计算设备的处理器执行以提供操作。该操作包括接收多元素系统的局部环境的多个旋转不变特征和多个旋转协变特征的组合。该操作还包括基于多个旋转不变特征和多个旋转协变特征的组合以及NNFF来预测用于多元素系统内的每个元素的力向量,以获得多元素系统内的元素的模拟运动。
在第三实施例中,公开了一种用于使用NNFF来模拟多元素系统内的元素的运动的计算机系统。该计算机系统包括模拟计算机,该模拟计算机具有用于执行计算机可读指令的处理器和用于维持计算机可执行指令的存储器。计算机可执行指令在由处理器执行时执行以下功能:接收多元素系统的局部环境的多个旋转不变特征和多个旋转协变特征的组合;和基于多个旋转不变特征和多个旋转协变特征的组合以及NNFF来预测用于多元素系统内的每个元素的力向量,以获得多元素系统内的元素的模拟运动。
附图说明
图1描绘了示出针对不同的现有技术分子动力学(MD)力计算方法的尺度(scale)和准确度之间的权衡的示意图形。
图2是现有技术Behler-Parrinello方法下的神经网络力场(NNFF)能量预测方案的神经网络节点示意图。
图3描绘了根据一个实施例的NNFF算法的神经网络节点示意。
图4是使用一个实施例的NNFF算法的预测力对比块聚环氧乙烷(PEO)的量子力学计算力的图形。
图5A示出了一个实施例的NNFF算法的指纹神经网络(NN)的图。
图5B、5C和5D描绘了示出使用一个实施例的NNFF算法的指纹NN来分别预测块PEO的H、C和O原子的结果的图形。
图6A示出了在一个实施例的NNFF算法中被组合以产生组合NN的指纹NN和逆投影NN的图。
图6B、6C和6D描绘了示出使用一个实施例的组合NNFF算法分别预测块PEO的H、C和O原子的结果的图形。
图7是在一个或多个实施例中可以用于实现NNFF算法(例如图6A的NNFF算法)的计算平台的示意图。
具体实施方式
根据需要,本文公开了本发明的详细实施例;然而,应当理解,所公开的实施例仅仅是本发明的示例,其可以以各种且替代形式来体现。附图不一定按比例绘制;一些特征可能被放大或最小化以示出特定组件的细节。因此,本文公开的具体结构和功能细节不应被解释为限制性的,而仅仅作为用于教导本领域技术人员以各种方式采用本发明的代表性基础。
术语“基本上”在本文中可以用于描述公开或要求保护的实施例。术语“基本上”可以修改在本公开中公开或要求保护的值或相对特性。在这种情况下,“基本上”可以表示其修改的值或相对特性在值或相对特性的±0%、0.1%、0.5%、1%、2%、3%、4%、5%或10%之内。
分子动力学(MD)方法对于研究物理现象是有益的,所述物理现象诸如但不限于离子迁移、化学反应、以及材料系统(诸如设备或功能材料)中的材料块和表面降解。这种材料系统的非限制性示例包括燃料电池、表面涂层、电池、水脱盐和水过滤。存在用于计算在原子运动模拟中使用的基础原子力的方法。从头算量子力学方法非常精确,但由于应用该方法所需的大量的计算资源,所以也非常昂贵。
神经网络已经被用来拟合和预测量子力学能量。这些方法被称为神经网络力场(NNFF)。使用量子力学能量来预测能量相对于原子位置和力的导数。然而,这些方法在计算上也很广泛。鉴于前述内容,需要的是一种提供适当水平的准确度同时消耗合理量的计算资源的用于计算原子力的计算方法。
分子动力学使用原子位置(以及可能为电荷、键或其他结构信息)来计算每个原子的原子间力,其因此用于修改模拟中的原子的速度。原子的结果所得轨迹用于描述物理现象,诸如但不限于电池和燃料电池中的离子迁移、块和表面材料降解期间的化学反应、以及固态材料相变。取决于用于计算原子力的基础方法,在模拟的准确度和尺寸(通过原子数和模拟动力学时间度量)之间存在权衡。如上所述,一种准确而昂贵的方法使用从头算量子力学方法,称为从头算分子动力学(AIMD)。
图1描绘了示出不同现有技术MD力计算方法(例如AIMD 102、力场MD 104和粗粒度MD 106)的尺度和准确度之间的权衡。y轴108是准确度,而x轴110是尺度(例如100个原子、1-10K个原子和100K个原子)。AIMD 102可以用于建模抗粘连涂层降解。通过力场拟合112,力场MD可以用于建模非晶态电解质和质子交换膜燃料电池(PEMFC)双极板(BPP)氧化物涂层。粗粒度MD 106可以用于建模PEMFC全氟磺酸聚合物。
AIMD针对许多相关应用在计算上要求太高。例如,根据一个估计,在具有两百(200)个核的高性能计算机上的AIMD模拟将花费六十五(65)年完成,以研究具有约五百(500)个原子的非晶态磷酸锂玻璃在五百(500)纳秒内的动力学行为。
AIMD的替代方案是经典分子动力学。使用该替代方案,原子间力被拟合成一些封闭的函数形式,以使得相比于AIMD加速几个数量级。这种替代方法已用于模拟大型系统,诸如聚合物和非晶态材料。根据一个估计,使用经典分子动力学方法的上述的相同模拟在不到两(2)天内完成。
使用预先限定的函数形式的经典分子动力学模型的缺点之一是,它们在某些条件下保持有效的假设内限制了原子间相互作用的行为。例如,用于原子之间的键的具有二次函数形式的经典力场假设了原子之间的化学键不会断裂。该假设与自然行为相反,在自然行为中分子中的化学键在化学反应期间形成和断裂。存在用于模拟化学反应的相对更复杂的预设函数形式。例如,反作用力场(诸如ReaxFF)是这样一种模型。这些模型经受难以参数化并且需要大量相对复杂的手工制作的封闭函数形式的缺点。
如上所述,存在需要计算量子力学能量的NNFF算法。NNFF的一种实现使用了Behler和Parrinello所发展的旋转不变特征(以其他方式称为Behler-Parrinello方法),如以下文章中所述:J. Behler和M. Parrinello, "Generalized Neural-NetworkRepresentation of High-Dimensional Potential-Energy Surfaces", Phys. Rev.Lett. 98, 146401 (2007);M. Gastegger, J. Behler, 和P. Marquetand, "MachineLearning Molecular Dynamics for Simulation of Infrared Spectra", Chern. Sci.2017, 8, 6924;以及J. Behler, "Atom-centered Symmetry Functions forConstructing High-Dimensional Neural Network Potentials", J. Chern. Phys.,134, 074106 (2011),全部通过引用以其整体并入本文。
在该实现中,在预定义半径Rc内围绕中心原子i的原子被用作向神经网络中的输入。不是将原子的原子笛卡尔坐标用作输入特征,而是将原子笛卡尔坐标转换为指纹特征向量。这些向量在不应改变局部原子结构的能量坐标的物理操作下是不变的。这些物理操作包括旋转、平移和置换(permutation)。旋转不变性是指相对于空间中的任意点(诸如中心原子i)以相同旋转来旋转局部原子环境的整个笛卡尔坐标时的能量不变性。平移不变性是指将局部原子环境的整个笛卡尔坐标移动相同位移时的能量不变性。置换不变性是指在改变局部原子环境中的原子的有序顺序并保持其位置固定时的能量不变性。除了满足这三(3)个条件之外,所选的指纹特征向量还可以满足以下准则(其被假设成与系统的物理性一致):来自远离中心原子i的原子的衰减贡献(向下到0,针对R>Rc);以及随着原子位置在局部原子环境中平滑地变化而平滑地变化。
其中,、、和是用户限定的G指纹函数的参数,如J. Behler, "Atom-centeredSymmetry Functions for Constructing High-Dimensional Neural NetworkPotentials", J. Chern. Phys., 134, 074106 (2011)中所述。通常选择这些参数的值,以适当覆盖考虑下的局部原子系统相关的长度尺度(通常在截止半径Rc<10 A之内)。截止函数fc取决于Rc并且可以被限定如下:
置换不变性仅对于相同元素的原子为真。例如,交换位置1处的C原子和位置2处的C原子描述了具有相同物理性的两个系统,因此,置换不变性成立。作为另一示例,交换位置1处的C原子和位置2处的O原子描述了具有不同物理性的两个系统,并且置换不变性不成立。因此,在研究多元素系统时,以上G函数的求和应仅针对具有相同元素(对于两体G)或相同元素对(对于三体G)的原子进行。给定围绕中心原子i的感兴趣的MD系统中的m数量的元素,则针对每个参数化的两体函数计算m个不同的G,并且针对每个参数化的三体函数计算m(m+1)/2个不同的G。由于这些函数的构造(包括旋转不变性),这些函数提供了用于学习量子力学的旋转不变性质的有益表示。然而,使用这些不变特征的机器学习系统限于用来预测旋转不变输出,诸如能量。
基于量子力学的许多方法(例如,密度泛函理论)被用于计算系统的能量E。然而,这些方法不能方便地将电子能分为来自单独原子的贡献。因此,在一个或多个实施例中,计算用于系统中的每个原子的特征向量,然后将其用于训练或预测单个标量能量量E,如图2所示。图2描绘了现有Behler-Parrinello方法下的NNFF能量预测方案的神经网络节点示意图。原子坐标环境R1至RN 202分别用作NNFF算法200的输入。原子坐标环境R1至RN 202用于计算对称函数G1至GN 204。每个原子坐标环境R1至RN 202与每个对称函数G1至GN 204相互关连。使用对称函数G1至GN 204训练原子NN 206以获得原子能208,其用于获得总能量E210。
用于AIMD应用的NNFF方法经受一个或多个限制。例如,在训练周期期间,每个原子NN(通常,每个元素一个NN)被共同训练为单个能量量E,其关联于感兴趣的系统。因此,NN是相互依赖的,因为每个需要共同用于预测单个能量值E,并且仅针对预测训练集中使用的确切MD系统而优化。没有独立的样本,其可以被包括以用于训练包含类似于的局部原子环境的不同MD系统。其次,在MD运行时间期间(强制预测周期),执行昂贵的计算以计算原子i上的力向量。是总系统能量E相对于单个原子的位移的导数。然而,预测E需要知道针对MD系统中的所有原子的G。因此,计算需要计算指纹向量3N次,如以下等式所示。
其中,指示笛卡尔x、y和z轴,N是MD系统中的原子的总数,并且是描述原子i的指纹函数的总数。因此,系统中的所有N个原子的计算力以增长。在某些实现中,某些算法修改可以将指纹向量计算的次数减少到3P次,其中P是原子的局部原子环境中的原子的平均数量,从而使计算成本尺度减少到。然而,此限制要求针对一个预测的计算(取决于所选的Rc,P的范围可能从100到400个原子),而不是直接预测中的计算。因此,具有直接预测的NNFF算法可以将基于NNFF的MD计算加速两(2)到三(3)个数量级。
在一个或多个实施例中,公开了一种计算机系统,该计算机系统被配置成实现一种算法,该算法使用神经网络直接从分子结构预测量子力学原子力而不使用能量导数。在某些实施例中,与需要量子力学能量的计算的已知NNFF算法相比,该算法的计算上需求可以少两(2)至(3)个数量级(例如,高性能计算(HPC)系统的更少计算时间)。
HPC系统执行某些任务的能力可以使用基准进行量化,诸如来自计算代码的应用基准,诸如标准性能评估公司的SPEC HPC2002和SPEC MPI2007套件(www.spec.org/hpg)和美国国家科学基金会(NSF)HPC系统获取过程(其由此通过引用以其整体并入本文)中所表示的那些。
在一个或多个实施例中,经典分子动力学模型的更简单的函数形式被替换为神经网络。可以根据提供的[输入,输出]数据集容易地调整这样的网络。围绕原子i的原子的局部环境可以用于预测原子i所经历的力。在通常的NNFF模型中,中心原子i被位移,然后计算原子i和i附近的所有其他原子的能量贡献的变化。总系统能量对原子i位移的数值导数是原子i所经历的原子力。该方法很耗时,因为需要计算来自围绕原子i的所有原子的能量贡献的变化以计算力。在一个或多个实施例中,本文公开的系统直接预测原子i上的力。该计算系统的速度可以是整个系统方法的两(2)到三(3)个量级那么快,因为无需计算来自围绕原子i的其他原子的附加贡献,同时维持必要的准确度。
在一个实施例中,公开了一种计算算法,其被配置成通过分析围绕中心原子的局部原子环境来直接预测原子上的力向量。可以在局部原子环境的旋转不变特征上训练神经网络。经训练的神经网络被配置成预测MD模拟内的原子力。在一个或多个实施例中,旋转不变特征可以适于预测相对少量的旋转不变中间状态。随后,可以将多个旋转协变特征添加到这些旋转不变中间状态,其然后可以用在神经网络的随后部分上以预测力向量值。力向量值表示旋转协变输出。
在一个实施例中,因为该计算算法使得能够对每个训练样本进行直接预测(不同于需要知道模拟帧中的所有数据样本的E预测),所以来自不同系统或材料集的样本可以被组合用于力场训练和/或发展。例如,一个或多个实施例的计算算法可以用于发展用于材料A和B之间的界面的大型原子系统的力场。使用该算法,可以通过组合材料A、材料B以及材料A和B之间的界面的单独较小结构来创建训练集。不能使用需要将E预测作为例程的一部分的算法来创建此训练集,因为从计算观点出发,从量子力学计算组合的大型原子系统的E太昂贵。
可以将一个或多个实施例的计算算法集成到MD模拟引擎中,以执行MD模拟,其准确度接近明显昂贵得多的AIMD方法。例如,一个或多个实施例的MD模拟可以有用于对与燃料电池、水脱盐、催化、涂层和电池中使用的功能材料相关的原子/离子动力学和化学反应进行建模。
在一个或多个实施例中,一个或多个实施例的原子力预测算法被集成到MD软件中。合适的MD软件的非限制性示例是大规模原子/分子大规模并行模拟器(LAMMPS)。LAMMPS可从位于新墨西哥阿尔伯克基的Sandia国家实验室获得。LAMMPS使用消息传递接口(MPI)来并行通信。LAMMPS是开源软件,根据GNU通用公共许可证的条款分布。
该算法能够进行直接力预测的主要原理是为每个原子坐标环境限定一组k个内轴,其随局部原子环境结构(旋转协变向量)旋转,然后将旋转协变力向量投影到上以获得内力向量。因为随后被投影回笛卡尔坐标,所以k是固定数且大于或等于3。即使和两者都是旋转协变的,这些向量也是相对于局部原子结构本身的刚性向量(两个向量都随结构旋转)。因此,原子位置的旋转变换对应地旋转且。因此,内部投影向量是旋转不变的。在一个或多个实施例中,旋转不变特征向量被用于预测旋转不变中间状态。旋转不变和旋转协变被用于执行逆投影并获得旋转协变输出笛卡尔中心原子力。下表1示出了不变和协变向量符号的概要。
限定了一组内轴,其中k指示围绕中心原子i的一组至少一些原子。原子j是否被视为组k的一部分的准则是用户限定的,并且可以取决于所考虑的MD系统进一步优化。一旦确定此原子分组准则,则应用如下约束,即:应是该组中的原子相对于中心原子i的所有相对位置向量的加权平均。该约束可以数学地写成:
其中,是如本文所述的轴函数。在一个或多个实施例中,设置与针对计算的Rc相同的针对计算的Rc。在其他实施例中,设置相对于局部原子坐标的变化更稳定的向量,因为仅用于向量投影参考轴。针对聚环氧乙烷(PEO)聚合物(C2H4O)n使用以下参数化函数形式来代替:
在用于PEO的NNFF发展的示例中,下表2中概述了分组准则:
一旦限定了从计算的方式,则将投影(使用量子力学)到上以找到一组内力投影。因为和两者都是旋转不变的,所以可以训练第一神经网络(NN1)来直接执行 预测。作为下一步,执行从内力回到笛卡尔原子力的逆投影。在一个实施例中,因为关系成立,所以使用矩阵伪逆关系执行逆投影以找到最小平方差解:。实际中,这可能不给出合适的结果,因为该组内轴可能不被约束为针对整个3D-xyz空间的正交(或接近于正交)基。因此,如果基向量重叠,则少量的噪声(无论是由于误差还是由于原子坐标的扰动)都可能导致数值不稳定性。虽然如果在预测中存在零误差则这不是问题,但NN1的预测中的任何误差都可能由于计算期间的伪逆方法而放大几个数量级。在一个或多个实施例中,这使得矩阵最小平方或伪逆方法不适合进行力向量的逆投影。
可以训练第二较小的神经网络(NN2)来解决逆投影问题。和用作针对NN2的输入,并且NN2将预测为输出。在一个或多个实施例中,NN2通过执行数据扩充来考虑某种程度的旋转对称性。在一个或多个实施例中,不仅,而且也被用作训练样本的一部分,其中是随机3D旋转矩阵(许多不同的可以用于NN2数据扩充)。
作为下一步,可以通过将NN1和NN2的权重和偏差转移到具有与NN1和NN2相同的节点维度的组合NN架构中(转移学习)来将NN1和NN2组合成单个神经网络。可以通过使用(训练数据集在组合网络上运行较少数量的训练周期来减少和/或最小化组合NN的误差。
一个或多个实施例中的上述的NNFF算法的修改使得力预测能够简化,其中旋转不变和旋转协变特征用于使得能够在仅给出局部原子环境坐标和围绕原子i的元素的情况下直接预测原子i上的原子力。图3描绘了根据一个实施例的NNFF算法300的神经网络节点示意。原子坐标环境R1到RN分别用作NNFF算法300的输入302a到302n。原子坐标环境R1到RN用于计算对称函数G1到GN(称为图3的对称函数304a到304n)和轴函数A1至AN(称为图3的对称函数306a至306n)。使用对称函数G1到GN训练第一神经网络NN1(在图3称为308)以获得内力向量至(称为图3的310a至310n)。轴函数A1至AN和内力向量到用于训练第二神经网络NN2(称为图3的312)。NNFF算法300的输出是笛卡尔坐标系中的原子力F1至FN(称为图3的314a至314n)。在一个或多个实施例中,图3的NNFF算法300没有训练数据样本之间的相互依赖性。因此,可以执行直接预测,这显著提高了力预测的计算速度。
图4是使用一个实施例的NNFF算法预测的402对比使用量子力学密度泛函理论(DFT)模拟所计算的404的图形400。图形400示出了使用两种不同方法的PEO的H原子上的力。NNFF算法和量子力学DFT模拟之间的误差可以被表达为。使用此表达式,平均误差为0.01(eV/A)2,其在可接受的范围内。在此示例中,在训练周期期间没有看到来自测试集的样本。
图5A示出了一个实施例的NNFF算法的指纹神经网络(NN)502的图500。指纹NN的旋转不变分量504对应于内力投影中间状态506的预测。图5B、5C和5D分别描绘了图形508、510和512,其示出了使用一个实施例的NNFF算法来分别预测从块PEO取得的原子快照的测试集中的H、C和O原子的结果。使用上述的误差方程,NNFF算法和量子力学DFT模拟之间的误差对于H、C和O原子的预测分别为0.01(eV/A)2、0.04(eV/A)2和0.16(eV/A)2,这在可接受的范围内。在这些示例中,在训练周期期间没有看到来自测试集的样本。
图6A示出了在一个实施例的NNFF算法中被组合以产生组合NN 606的指纹NN 602和逆投影NN 604的图600。根据一个或多个实施例,NN 602和604的权重被转移并且组合的NN 606被重新训练。指纹NN 602是NNFF算法的旋转不变分量,而逆投影NN 604是NNFF算法的旋转协变分量。组合的NN 606预测笛卡尔原子力。图6B、6C和6D分别描绘了图形608、610和612,其示出了使用一个实施例的NNFF算法来分别预测从块PEO取得的原子快照的测试集中的H、C和O原子的结果。表示到笛卡尔x、y和z轴上的投影。使用上述的误差方程,NNFF算法和量子力学DFT模拟之间的误差对于H、C和O原子的预测分别为0.03(eV/A)2、0.58(eV/A)2和0.42(eV/A)2,这在可接受的范围内。在这些示例中,在训练周期期间没有看到来自测试集的样本。
一个或多个实施例的NNFF算法和/或方法使用诸如图7所示的计算平台700之类的计算平台来实现。计算平台700可以包括处理器702、存储器704和非易失性存储706。处理器702可以包括选自高性能计算(HPC)系统的一个或多个设备,包括高性能核、微处理器、微控制器、数字信号处理器、微计算机、中心处理单元、现场可编程门阵列、可编程逻辑器件、状态机、逻辑电路、模拟电路、数字电路、或基于驻留在存储器704中的计算机可执行指令来操纵信号(模拟或数字)的任何其他设备。存储器704可以包括单个存储器设备或多个存储器设备,包括但不限于随机存取存储器(RAM)、易失性存储器、非易失性存储器、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、闪速存储器、高速缓存存储器、或能够存储信息的任何其他设备。非易失性存储器706可以包括一个或多个持久性数据存储设备,诸如硬盘驱动器、光盘驱动器、磁带驱动器、非易失性固态设备、云存储、或能够持久地存储信息的任何其他设备。
处理器702可以被配置成读入存储器704并执行驻留在非易失性存储706的NNFF软件模块708中并且体现一个或多个实施例的NNFF算法和/或方法的计算机可执行指令。处理器702可以进一步被配置成读入存储器704并执行驻留在非易失性存储706的MD软件模块710(诸如LAMMPS)中并且体现MD算法和/或方法的计算机可执行指令。软件模块708和710可以包括操作系统和应用。软件模块708和710可以从使用多种编程语言和/或技术所创建的计算机程序来编译或解释,包括但不限于并且单独或组合Java、C、C++、C#、Objective C、Fortran、Pascal、Java Script、Python、Perl和PL/SQL。
在由处理器702执行时,NNFF软件模块708和MD软件模块710的计算机可执行指令可以使计算平台700实现本文分别公开的NNFF算法和/或方法以及MD算法和/或方法中的一种或多种。非易失性存储器706还可以包括支持本文描述的一个或多个实施例的功能、特征和过程的NNFF数据712和MD数据714。
体现本文描述的算法和/或方法的程序代码能够作为程序产品以各种不同的形式单独或共同分布。可以使用其上具有用于使处理器执行一个或多个实施例的各方面的计算机可读程序指令的计算机可读存储介质来分布程序代码。本质上是非临时性的计算机可读存储介质可以包括易失性和非易失性以及可移动和不可移动的有形介质,其采用用于存储诸如计算机可读指令、数据结构、程序模块或其他数据之类的信息的任何方法或技术来实现。计算机可读存储介质还可以包括RAM、ROM、可擦可编程只读存储器(EPROM)、电可擦可编程只读存储器(EEPROM)、闪速存储器或其他固态存储技术、便携式压缩盘只读存储器(CD-ROM)或其他光学存储、盒式磁带、磁带、磁盘存储或其他磁存储设备、或可以用于存储所需信息并可以由计算机读取的任何其他介质。可以将计算机可读程序指令从计算机可读存储介质下载到计算机、另一种类型的可编程数据处理装置、或另一设备,或者经由网络下载到外部计算机或外部存储设备。
存储在计算机可读介质中的计算机可读程序指令可以用于指导计算机、其他类型的可编程数据处理装置、或其他设备以特定方式起作用,使得存储在计算机可读介质中的指令产生制造品,包括实现流程图或图中指定的功能、动作和/或操作的指令。在某些替代实施例中,流程图和图中指定的功能、动作和/或操作可以按照一个或多个实施例被重新排序、串行处理、和/或同时处理。此外,任何流程图和/或图可以包括比按照一个或多个实施例图示的节点或框更多或更少的节点或框。
虽然已经通过各种实施例的描述说明了本发明的全部,并且虽然已经对这些实施例进行了相当详细的描述,但是申请人的意图不是将所附权利要求的范围限制或以任何方式限于这样的细节。附加的优点和修改对于本领域技术人员将是显而易见的。因此,本发明在其更广泛的方面不限于所示出和描述的具体细节、代表性装置和方法、以及说明性示例。因此,在不脱离总的发明构思的精神或范围的情况下,可以与这样的细节有所不同。
Claims (20)
1.一种用于使用神经网络力场(NNFF)来模拟多元素系统内的元素的运动的计算方法,所述方法包括:
接收多元素系统的局部环境的多个旋转不变特征和多个旋转协变特征的组合;和
基于多个旋转不变特征和多个旋转协变特征的组合以及NNFF来预测用于多元素系统内的每个元素的力向量,以获得多元素系统内的元素的模拟运动。
2.根据权利要求1所述的计算方法,还包括基于多个旋转不变特征来训练第一神经网络(NN1)以预测多个旋转不变中间状态。
3.根据权利要求2所述的计算方法,还包括基于多个旋转不变中间状态和多个旋转协变特征来训练第二神经网络(NN2)以预测多个旋转协变输出。
4.根据权利要求3所述的计算方法,还包括将NN1和NN2集中成第三神经网络(NN3),其用于获得多元素系统内的元素的模拟运动。
5.根据权利要求4所述的计算方法,其中,所述多元素系统是材料系统,所述元素是所述材料系统内的原子,并且所述局部环境是原子局部环境。
9.根据权利要求1所述的计算方法,其中,所述多个旋转不变特征包括一个或多个Behler-Parrinello型指纹特征向量。
10.根据权利要求1所述的计算方法,其中,所述多元素系统是以下系统之一的一部分:燃料电池、水脱盐、催化、涂层和电池。
11.一种非暂时性计算机可读介质,其有形地体现用于软件程序的计算机可读指令,所述软件程序可由计算设备的处理器执行以提供包括以下内容的操作:
接收具有元素的多元素系统的局部环境的多个旋转不变特征和多个旋转协变特征的组合;和
基于多个旋转不变特征和多个旋转协变特征的组合以及神经网络力场(NNFF)来预测用于多元素系统内的每个元素的力向量,以获得多元素系统内的元素的模拟运动。
12.根据权利要求11所述的非暂时性计算机可读介质,其中,所述软件程序可由所述计算设备的处理器执行以提供基于多个旋转不变特征来训练第一神经网络(NN1)以预测多个旋转不变中间状态的进一步操作。
13.根据权利要求12所述的非暂时性计算机可读介质,其中,所述软件程序可由所述计算设备的处理器执行以提供基于多个旋转不变中间状态和多个旋转协变特征来训练第二神经网络(NN2)以预测多个旋转协变输出的进一步操作。
14.根据权利要求13所述的非暂时性计算机可读介质,其中,所述软件程序可由所述计算设备的处理器执行以提供将NN1和NN2集中成用于获得多元素系统内的元素的模拟运动的第三神经网络(NN3)的进一步操作。
15.根据权利要求11所述的非暂时性计算机可读介质,其中,所述多个旋转不变特征包括一个或多个Behler-Parrinello型指纹特征向量。
16.一种用于使用神经网络力场(NNFF)来模拟多元素系统内的元素的运动的计算机系统,其包括具有用于执行计算机可读指令的处理器和用于维持计算机可执行指令的存储器的模拟计算机,所述计算机可执行指令在由所述处理器执行时执行以下功能:
接收具有元素的多元素系统的局部环境的多个旋转不变特征和多个旋转协变特征的组合;和
基于多个旋转不变特征和多个旋转协变特征的组合以及NNFF来预测用于多元素系统内的每个元素的力向量,以获得多元素系统内的元素的模拟运动。
17.根据权利要求16所述的计算机系统,其中,所述计算机可执行指令在由所述处理器执行时执行基于多个旋转不变特征来训练第一神经网络(NN1)以预测多个旋转不变中间状态的进一步功能。
18.根据权利要求17所述的计算机系统,其中,所述计算机可执行指令在由所述处理器执行时执行基于多个旋转不变中间状态和多个旋转协变特征来训练第二神经网络(NN2)以预测多个旋转协变输出的进一步功能。
19.根据权利要求18所述的计算机系统,其中,所述计算机可执行指令在由所述处理器执行时执行将NN1和NN2集中成用于获得多元素系统内的元素的模拟运动的第三神经网络(NN3)的进一步功能。
20.根据权利要求16所述的计算机系统,其中,所述多个旋转不变特征包括一个或多个Behler-Parrinello型指纹特征向量。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US16/202,890 US11455439B2 (en) | 2018-11-28 | 2018-11-28 | Neural network force field computational algorithms for molecular dynamics computer simulations |
US16/202890 | 2018-11-28 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111241655A true CN111241655A (zh) | 2020-06-05 |
Family
ID=70546431
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911182094.8A Pending CN111241655A (zh) | 2018-11-28 | 2019-11-27 | 用于分子动力学计算机模拟的神经网络力场计算算法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US11455439B2 (zh) |
CN (1) | CN111241655A (zh) |
DE (1) | DE102019217222A1 (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102017127276A1 (de) * | 2017-08-30 | 2019-02-28 | Taiwan Semiconductor Manufacturing Co., Ltd. | Standardzellen und abwandlungen davon innerhalb einer standardzellenbibliothek |
US11734478B2 (en) * | 2019-03-18 | 2023-08-22 | Lawrence Livermore National Security, Llc | Spliced soft-core interaction potential for filling small-scale enclosures |
CN111875826B (zh) * | 2020-07-22 | 2022-11-25 | 上海大学 | 粗粒化分子动力学/蒙特卡洛模拟方法及质子交换膜材料 |
CN113515854B (zh) * | 2021-06-10 | 2022-09-27 | 太原理工大学 | 一种评价荷电活化水雾雾化降尘性能的分子动力学方法 |
CN113704961B (zh) * | 2021-07-14 | 2024-04-19 | 湖南大学 | 一种质子交换膜燃料电池瞬态运行性能的预测方法 |
US20230095631A1 (en) * | 2021-09-24 | 2023-03-30 | Robert Bosch Gmbh | Neural network force field computational training routines for molecular dynamics computer simulations |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20120085136A (ko) * | 2011-10-04 | 2012-07-31 | 주식회사 켐에쎈 | 순수한 유기화합물의 표준상태의 이상기체 절대 엔트로피를 예측하는 다중선형회귀-인공신경망 모형 |
KR20120085148A (ko) * | 2011-10-05 | 2012-07-31 | 주식회사 켐에쎈 | 순수한 유기화합물의 표준상태 생성에너지를 예측하는 다중선형회귀-인공신경망 모형 |
CN104715096A (zh) * | 2013-12-12 | 2015-06-17 | 中国科学院大连化学物理研究所 | Bp神经网络预测二肽模型多极展开属性计算方法 |
CN108140131A (zh) * | 2015-10-04 | 2018-06-08 | 艾腾怀斯股份有限公司 | 用于将卷积网络应用于空间数据的系统和方法 |
CN108804869A (zh) * | 2018-05-04 | 2018-11-13 | 深圳晶泰科技有限公司 | 基于神经网络的分子结构和化学反应能量函数构建方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0810413D0 (en) | 2008-06-06 | 2008-07-09 | Cambridge Entpr Ltd | Method and system |
-
2018
- 2018-11-28 US US16/202,890 patent/US11455439B2/en active Active
-
2019
- 2019-11-07 DE DE102019217222.9A patent/DE102019217222A1/de active Pending
- 2019-11-27 CN CN201911182094.8A patent/CN111241655A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20120085136A (ko) * | 2011-10-04 | 2012-07-31 | 주식회사 켐에쎈 | 순수한 유기화합물의 표준상태의 이상기체 절대 엔트로피를 예측하는 다중선형회귀-인공신경망 모형 |
KR20120085148A (ko) * | 2011-10-05 | 2012-07-31 | 주식회사 켐에쎈 | 순수한 유기화합물의 표준상태 생성에너지를 예측하는 다중선형회귀-인공신경망 모형 |
CN104715096A (zh) * | 2013-12-12 | 2015-06-17 | 中国科学院大连化学物理研究所 | Bp神经网络预测二肽模型多极展开属性计算方法 |
CN108140131A (zh) * | 2015-10-04 | 2018-06-08 | 艾腾怀斯股份有限公司 | 用于将卷积网络应用于空间数据的系统和方法 |
CN108804869A (zh) * | 2018-05-04 | 2018-11-13 | 深圳晶泰科技有限公司 | 基于神经网络的分子结构和化学反应能量函数构建方法 |
Non-Patent Citations (4)
Title |
---|
BEHLER, J, ET AL.: "Generalized Neural-Network Representation of High-Dimensional Potential-Energy Surfaces", 《PHYSICAL REVIEW LETTERS》, 2 April 2007 (2007-04-02), pages 1 - 4 * |
BEHLER, J: "Atom-centered symmetry functions for constructing high-dimensional neural network potentials", 《JOURNAL OF CHEMICAL PHYSICS》, 16 February 2011 (2011-02-16), pages 1 - 13 * |
LI, H, ET AL.: "Application of Artificial Neural Networks for Catalysis: A Review. Catalysts", 《CATALYSTS》, 18 October 2017 (2017-10-18), pages 1 - 19 * |
雷博;饶金理;覃武;何克波;唐玲;霍明;江建军;: "分子动力学和人工神经网络的算法耦合研究", 功能材料, no. 02, 20 February 2009 (2009-02-20), pages 284 - 290 * |
Also Published As
Publication number | Publication date |
---|---|
DE102019217222A1 (de) | 2020-05-28 |
US20200167439A1 (en) | 2020-05-28 |
US11455439B2 (en) | 2022-09-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111241655A (zh) | 用于分子动力学计算机模拟的神经网络力场计算算法 | |
AU2017414716B2 (en) | Operator averaging within quantum computing systems | |
Thomas et al. | Tensor field networks: Rotation-and translation-equivariant neural networks for 3d point clouds | |
Thompson et al. | Spectral neighbor analysis method for automated generation of quantum-accurate interatomic potentials | |
CN110023966B (zh) | 使用量子计算仿真材料 | |
CN111951899A (zh) | 用于分子动力学计算机模拟的图形神经网络力场计算算法 | |
US11170141B2 (en) | Graph transformer neural network force field for prediction of atomic forces and energies in molecular dynamic simulations | |
JP7305714B2 (ja) | 量子シミュレーションのための平面波双対基底 | |
Qiao et al. | Unite: Unitary n-body tensor equivariant network with applications to quantum chemistry | |
CN117747012A (zh) | 一种晶体性质预测方法、装置、计算机设备和储存介质 | |
Gubaev et al. | Performance of two complementary machine-learned potentials in modelling chemically complex systems | |
Ellis et al. | A non-critical string approach to black holes, time and quantum dynamics | |
Zeng et al. | A nearsighted force-training approach to systematically generate training data for the machine learning of large atomic structures | |
US11815945B2 (en) | Efficient scaling of neural-network interatomic potential prediction on CPU clusters | |
Hu et al. | Efficient Machine Learning Force Field for Large-Scale Molecular Simulations of Organic Systems | |
US20230359929A1 (en) | Orbital mixer machine learning method for predicting an electronic structure of an atomic system | |
Gubaev et al. | Exploring the limits of machine-learned potentials for chemically complex multicomponent systems | |
Pandey et al. | Evolutionary Analysis: Classical Bits to Quantum Qubits | |
Kashyap | Quantum Computers, Computing, and Machine Learning: A Review | |
Takahashi et al. | Influence of Recognition Performance on Recurrent Neural Network Using Phase-Change Memory as Synapses | |
Lou | Hybrid Path Integral Monte Carlo/Molecular Dynamics Approach for the Simulation of Rigid Rotating and Translating Molecules | |
Ahmed et al. | NeuSpin: Design of a Reliable Edge Neuromorphic System Based on Spintronics for Green AI | |
Smith | Fast Machine-Learned Simulations for Future Detectors | |
CN115790822A (zh) | 一种确定分子振动频率的方法及装置 | |
CN115829039A (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 |