CN113936746A - 多原子体系射线衍射的快速分析方法、系统、终端及介质 - Google Patents

多原子体系射线衍射的快速分析方法、系统、终端及介质 Download PDF

Info

Publication number
CN113936746A
CN113936746A CN202111535638.1A CN202111535638A CN113936746A CN 113936746 A CN113936746 A CN 113936746A CN 202111535638 A CN202111535638 A CN 202111535638A CN 113936746 A CN113936746 A CN 113936746A
Authority
CN
China
Prior art keywords
dimension
space
lattice point
summation
calculation
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.)
Granted
Application number
CN202111535638.1A
Other languages
English (en)
Other versions
CN113936746B (zh
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.)
Institute of Fluid Physics of CAEP
Original Assignee
Institute of Fluid Physics of CAEP
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 Institute of Fluid Physics of CAEP filed Critical Institute of Fluid Physics of CAEP
Priority to CN202111535638.1A priority Critical patent/CN113936746B/zh
Publication of CN113936746A publication Critical patent/CN113936746A/zh
Application granted granted Critical
Publication of CN113936746B publication Critical patent/CN113936746B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/20—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by using diffraction of the radiation by the materials, e.g. for investigating crystal structure; by using scattering of the radiation by the materials, e.g. for investigating non-crystalline materials; by using reflection of the radiation by the materials
    • G01N23/207—Diffractometry using detectors, e.g. using a probe in a central position and one or more displaceable detectors in circumferential positions
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/141—Discrete Fourier transforms
    • G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/17—Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F9/00—Arrangements for program control, e.g. control units
    • G06F9/06—Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
    • G06F9/46—Multiprogramming arrangements
    • G06F9/50—Allocation of resources, e.g. of the central processing unit [CPU]
    • G06F9/5005—Allocation of resources, e.g. of the central processing unit [CPU] to service a request
    • G06F9/5027—Allocation of resources, e.g. of the central processing unit [CPU] to service a request the resource being a machine, e.g. CPUs, Servers, Terminals
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B15/00—ICT specially adapted for analysing two-dimensional or three-dimensional molecular structures, e.g. structural or functional relations or structure alignment
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B5/00—ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B50/00—ICT programming tools or database systems specially adapted for bioinformatics
    • G16B50/30—Data warehousing; Computing architectures
    • 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/20—Identification of molecular entities, parts thereof or of chemical compositions
    • 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/90—Programming languages; Computing architectures; Database systems; Data warehousing
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Health & Medical Sciences (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Chemical & Material Sciences (AREA)
  • Software Systems (AREA)
  • General Health & Medical Sciences (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Databases & Information Systems (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Evolutionary Biology (AREA)
  • Biotechnology (AREA)
  • Medical Informatics (AREA)
  • Biophysics (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Bioethics (AREA)
  • Discrete Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Molecular Biology (AREA)
  • Physiology (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明公开了多原子体系射线衍射的快速分析方法、系统、终端及介质，涉及射线衍射技术领域，其技术方案要点是：建立分维度编号的二维或三维的实空间的格点阵列；确定原子与距离最近的实空间格点的关联；建立与衍射实验范围和精度对应的倒空间的格点阵列；采用泰勒展开对传统算法作近似处理，分别对各维度的变量求和计算出倒空间格点的振幅；以倒空间格点的振幅作为输入，并通过插值计算得到任意转移动量上的振幅。本发明通过采用泰勒展开的方式，使得基于原子位置的射线衍射计算可以按空间维度分离变量求和，每次求和只涉及一个维度，大大减小计算量的同时得到精确度按展开阶数指数上升的近似结果。

Description

多原子体系射线衍射的快速分析方法、系统、终端及介质
技术领域
本发明涉及射线衍射技术领域，更具体地说，它涉及多原子体系射线衍射的快速分析方法、系统、终端及介质。
背景技术
射线衍射是探测物质微介观结构的重要手段，溶液X射线衍射分析等技术已被广泛应用于DNA结构分析等领域，只在散射角一个维度上开展的射线衍射实验和数值模拟技术已经发展得非常成熟。随着相干X射线光源的发展，利用多维衍射实验分析物质微介观结构已逐渐成为可能，而与之相应的数值模拟技术则还停留在实验室阶段。多维衍射能够提供更多的材料微介观信息，其复杂性远高于一维衍射，现有数值模拟技术并不能直接移植。
现有的多维衍射技术中，只有相干衍射成像的相关数值模拟技术发展得较为成熟，相干成像关心的客体在每一个单元内部被视作是均匀的，而原子则更像一个点而不像一个均匀的单元，因此相关数值模拟技术难以移植。而现有的基于原子位置的射线衍射计算方法，如L. Soulard, et al., Observation of phase transitions in shocked tinby molecular dynamics, 2020；孙伟，等，分子动力学模拟纳米晶体铜的结构与性能，1998，对每个转移动量、对每个原子求和，或者甚至对原子编号二次求和，计算代价随体系变大增长很快，应用到微米尺度的微介观物质所需的计算代价难以接受。
因此，如何研究设计一种能够克服上述缺陷的多原子体系射线衍射的快速分析方法、系统、终端及介质是我们目前急需解决的问题。
发明内容
为解决现有技术中的不足，本发明的目的是提供多原子体系射线衍射的快速分析方法、系统、终端及介质，通过采用泰勒展开的方式，使得基于原子位置的射线衍射计算可以按空间维度分离变量求和，每次求和只涉及一个维度，大大减小计算量的同时得到精确度按展开阶数指数上升的近似结果，可以在现代超级计算机上实现尺度可达微米以上的生物、医药等样品的精确到原子层面的射线衍射数值模拟，现代超级计算机如具有GPU-CPU异构的计算机。
本发明的上述技术目的是通过以下技术方案得以实现的：
第一方面，提供了多原子体系射线衍射的快速分析方法，包括以下步骤：
建立分维度编号的二维或三维的实空间的格点阵列；
确定原子与距离最近的实空间格点的关联；
建立与衍射实验范围和精度对应的倒空间的格点阵列；
采用泰勒展开对传统算法作近似处理，分别对各维度的变量求和计算出倒空间格点的振幅；
以倒空间格点的振幅作为输入，并通过插值计算得到任意转移动量上的振幅。
进一步的，所述泰勒展开的计算公式具体为：
Figure 775715DEST_PATH_IMAGE001
其中，
Figure 709167DEST_PATH_IMAGE002
表示自然对数的底数；
Figure 488904DEST_PATH_IMAGE003
表示倒空间中格点上的转移动量在x方向上的分量；
Figure 946430DEST_PATH_IMAGE004
是任一原子位置与实空间关联格点之间的位置差在x方向上的分量；
Figure 59005DEST_PATH_IMAGE005
是单位虚数；
Figure 209364DEST_PATH_IMAGE006
是泰勒展开的阶数，其值越大精确程度越高，但计算代价越大。
进一步的，所述各维度的变量求和包括第一个维度求和、第二个维度求和以及第三个维度求和；
第一个维度求和的计算公式为：
Figure 273135DEST_PATH_IMAGE007
其中，
Figure 268773DEST_PATH_IMAGE008
、
Figure 298171DEST_PATH_IMAGE009
分别表示另两个维度上的格点编号；
Figure 353851DEST_PATH_IMAGE010
表示与格点关联的不同原子的编号；
Figure 373760DEST_PATH_IMAGE011
表示倒空间的格点在第一个维度上的编号；
Figure 235406DEST_PATH_IMAGE012
表示倒空间格点上的转移动量在第一个维度上的分量；
Figure 322573DEST_PATH_IMAGE013
表示原子的散射强度；
Figure 283575DEST_PATH_IMAGE014
表示原子位置坐标在第一个维度上的分量；
Figure 587518DEST_PATH_IMAGE015
、
Figure 190537DEST_PATH_IMAGE016
分别表示原子与关联格点的位置差在第二、第三两个维度上的分量；幂指数
Figure 302850DEST_PATH_IMAGE017
和
Figure 467377DEST_PATH_IMAGE018
从零一直取到泰勒展开的阶数；
Figure 727457DEST_PATH_IMAGE019
表示第一个维度求和的计算结果；
第二个维度求和的计算公式为：
Figure 134168DEST_PATH_IMAGE020
其中，
Figure 694462DEST_PATH_IMAGE021
表示倒空间的格点在第二个维度上的编号；
Figure 537612DEST_PATH_IMAGE022
表示倒空间格点上的转移动量在第二个维度上的分量；
Figure 816146DEST_PATH_IMAGE023
表示在第二维第三维上编号分别为
Figure 495389DEST_PATH_IMAGE024
和
Figure 910190DEST_PATH_IMAGE025
的实空间格点坐标在第二个维度上的分量；
Figure 321580DEST_PATH_IMAGE026
表示第二个维度求和的计算结果；
第三个维度求和的计算公式为：
Figure 57717DEST_PATH_IMAGE027
其中，
Figure 806230DEST_PATH_IMAGE028
表示倒空间的格点在第三个维度上的编号；
Figure 75538DEST_PATH_IMAGE029
表示倒空间格点上的转移动量在第三个维度上的分量；
Figure 48042DEST_PATH_IMAGE030
表示在第三维上编号为
Figure 271475DEST_PATH_IMAGE031
的实空间格点坐标在第三个维度上的分量；
Figure 558100DEST_PATH_IMAGE032
表示第三个维度求和的计算结果。
进一步的，若所述实空间的格点阵列为二维，则依据泰勒展开直接进行变量求和；若所述实空间的格点阵列为三维，则可依据泰勒展开直接进行变量求和，或在第一个维度上建立的格点基础上进行变量求和。
进一步的，若所述实空间的格点阵列和/或倒空间的格点阵列是二维，则每个格点有两个编号；
若所述格点阵列是三维，则每个格点有三个编号；
任意一个编号不同，表示不同的格点；
格点阵列覆盖的空间范围和稀密程度由涉及的转移动量范围和计算要求的精确度决定。
进一步的，所述原子与距离最近的实空间格点的关联确定过程具体为：
为每个格点建立一个原子基团列表；
原子基团列表中的原子基团中存在至少一个原子在第二、三维度上的投影到相应格点的距离最近，而非离另一个格点更近；
构建一个调用函数，通过调用函数对每个格点遍历关联的原子基团中与该格点关联的原子。
进一步的，所述各维度的变量求和计算采用分布式计算和存储开展并行计算。
第二方面，提供了多原子体系射线衍射的快速分析系统，包括：
第一阵列模块，用于建立分维度编号的二维或三维的实空间的格点阵列；
关联确定模块，用于确定原子与距离最近的实空间格点的关联；
第二阵列模块，用于建立与衍射实验范围和精度对应的倒空间的格点阵列；
分维计算模块，用于采用泰勒展开对传统算法作近似处理，分别对各维度的变量求和计算出倒空间格点的振幅；
插值计算模块，用于以倒空间格点的振幅作为输入，并通过插值计算得到任意转移动量上的振幅。
第三方面，提供了一种计算机终端，包含存储器、处理器及存储在存储器并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如第一方面中任意一项所述的多原子体系射线衍射的快速分析方法。
第四方面，提供了一种计算机可读介质，其上存储有计算机程序，其特征是，所述计算机程序被处理器执行可实现如第一方面中任意一项所述的多原子体系射线衍射的快速分析方法。
与现有技术相比，本发明具有以下有益效果：
1、本发明提出的多原子体系射线衍射的快速分析方法，通过采用泰勒展开的方式，使得基于原子位置的射线衍射计算可以按空间维度分离变量求和，每次求和只涉及一个维度，大大减小计算量的同时得到精确度按展开阶数指数上升的近似结果，可以在现代超级计算机上实现尺度可达微米以上的生物、医药等样品的精确到原子层面的射线衍射数值模拟，现代超级计算机如具有GPU-CPU异构的计算机；尤其是在多角度衍射、复色光衍射、大体系衍射等多种情形下，所需的计算代价显著小于传统的直接对每个原子和转移动量求和的计算方法；且原子数量越多，转移动量点越多，计算分析方法在减小计算代价上的益处越显著；
2、本发明对于不同的转移动量，每个原子的位置也只需要访问一次，非常适用于原子数量特别大、不适合将所有位置存储起来、只能临时算出来的情形；如果采用传统方法，原子位置需要计算多次，而本发明中的方法只需要计算一次。
3、本发明提供的分析方法易实现并行化和分布计算。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解，构成本申请的一部分，并不构成对本发明实施例的限定。在附图中：
图1是本发明实施例中的流程图；
图2是本发明实施例中的系统框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。
实施例1：多原子体系射线衍射的快速分析方法，以染色体为例，如图1所示，具体由以下步骤实现。
步骤一、建立分维度编号的二维或三维的实空间的格点阵列。
若格点阵列是二维，则每个格点有两个编号；若格点阵列是三维，则每个格点有三个编号；任意一个编号不同，表示不同的格点；格点阵列覆盖的空间范围和稀密程度由涉及的转移动量范围和计算要求的精确度决定。
选择转移动量分布范围最小的方向为第一个维度；以垂直于第一个维度的方向为第二个维度和第三个维度。由于在第二、三两维上拟采用快速傅里叶变换，在第二、三维度上建立大小为1.024μm x 1.024μm 二维实空间格点阵列，格点间距为1nm。
步骤二、确定原子与距离最近的实空间格点的关联。
由于原子坐标是按照碱基对存储的，为了建立原子与实空间格点的关联，为每个格点建立一个碱基对列表，每个列表中包含且仅包含如下的碱基对：碱基对中存在至少一个原子在第二、三维度上的投影到该格点的距离最近而不是离另一个格点更近；构建一个调用函数，是其可以对每个格点遍历关联的碱基对中与该格点关联的原子。
步骤三、建立与衍射实验范围和精度对应的倒空间的格点阵列。
倒空间格点在第一维上的边界以转移动量在该维度上的分布范围为边界，间距取1/μm；在第二、三两个维度上的边界取为0-1/nm，该两个纬度每个维度上的数量为1024个。
步骤四、采用泰勒展开对传统算法作近似处理，分别对各维度的变量求和计算出倒空间格点的振幅。
泰勒展开的计算公式具体为：
Figure 104750DEST_PATH_IMAGE001
其中，
Figure 248155DEST_PATH_IMAGE002
表示自然对数的底数；
Figure 457420DEST_PATH_IMAGE003
表示倒空间中格点上的转移动量在x方向上的分量；
Figure 252463DEST_PATH_IMAGE004
是任一原子位置与实空间关联格点之间的位置差在x方向上的分量；
Figure 293100DEST_PATH_IMAGE005
是单位虚数；
Figure 810669DEST_PATH_IMAGE006
是泰勒展开的阶数，其值越大精确程度越高，但计算代价越大。在本实施例中，泰勒级数的阶数取5-10阶。
需要说明的是，若实空间的格点阵列为二维，则依据泰勒展开直接进行变量求和；若实空间的格点阵列为三维，则可依据泰勒展开直接进行变量求和，或在第一个维度上建立的格点基础上进行变量求和。
各维度的变量求和包括第一个维度求和、第二个维度求和以及第三个维度求和。
（1）第一个维度求和的计算公式为：
Figure 507229DEST_PATH_IMAGE007
其中，
Figure 637122DEST_PATH_IMAGE008
、
Figure 407632DEST_PATH_IMAGE009
分别表示另两个维度上的格点编号；
Figure 96102DEST_PATH_IMAGE010
表示与格点关联的不同原子的编号；
Figure 14379DEST_PATH_IMAGE011
表示倒空间的格点在第一个维度上的编号；
Figure 712077DEST_PATH_IMAGE012
表示倒空间格点上的转移动量在第一个维度上的分量；
Figure 166454DEST_PATH_IMAGE013
表示原子的散射强度；
Figure 822564DEST_PATH_IMAGE014
表示原子位置坐标在第一个维度上的分量；
Figure 696979DEST_PATH_IMAGE015
、
Figure 463946DEST_PATH_IMAGE016
分别表示原子与关联格点的位置差在第二、第三两个维度上的分量；幂指数
Figure 38409DEST_PATH_IMAGE017
和
Figure 803103DEST_PATH_IMAGE018
从零一直取到泰勒展开的阶数；
Figure 633656DEST_PATH_IMAGE019
表示第一个维度求和的计算结果。
需要说明的是，第一个维度求和的计算量最大的，在这一步中，倒空间格点在第一个维度上编号不同的涉及的计算完全不相关，采用分布式计算和存储开展并行计算，分布式存储和计算在第E步之前无需交互；在计算中需要遍历与格点关联的原子时，通过调用步骤二中的调用函数得到原子信息和位置；由于各格点关联的原子互不相同，这一步采用CPU或GPU并行加速；式3中m和n不同的，也采用CPU或GPU加速，并且优先采用GPU加速。
（2）第二个维度求和的计算公式为：
Figure 735473DEST_PATH_IMAGE020
其中，
Figure 430022DEST_PATH_IMAGE021
表示倒空间的格点在第二个维度上的编号；
Figure 631196DEST_PATH_IMAGE022
表示倒空间格点上的转移动量在第二个维度上的分量；
Figure 11362DEST_PATH_IMAGE023
表示在第二维第三维上编号分别为
Figure 588974DEST_PATH_IMAGE024
和
Figure 138029DEST_PATH_IMAGE025
的实空间格点坐标在第二个维度上的分量；
Figure 978946DEST_PATH_IMAGE026
表示第二个维度求和的计算结果。
需要说明的是，各b指标都需要对m进行求和，这一步也优先采用GPU并行加速；不同的A采用CPU并行计算，而指标b到B的转换则采用快速傅立叶变换技术。
（3）第三个维度求和的计算公式为：
Figure 846408DEST_PATH_IMAGE027
其中，
Figure 962131DEST_PATH_IMAGE028
表示倒空间的格点在第三个维度上的编号；
Figure 162431DEST_PATH_IMAGE029
表示倒空间格点上的转移动量在第三个维度上的分量；
Figure 439828DEST_PATH_IMAGE030
表示在第三维上编号为
Figure 60165DEST_PATH_IMAGE031
的实空间格点坐标在第三个维度上的分量；
Figure 714001DEST_PATH_IMAGE032
表示第三个维度求和的计算结果。
需要说明的是，不同的指标A和B采用CPU并行计算，从c到C的计算采用快速傅里叶变换。
步骤五、此步骤计算量很小，先将步骤四中分布存储和计算的数据归集起来之后，再开展计算。以倒空间格点的振幅作为输入，并通过插值计算得到任意转移动量上的振幅。
实施例2：多原子体系射线衍射的快速分析系统，如图2所示，包括第一阵列模块、关联确定模块、第二阵列模块、分维计算模块和插值计算模块。
其中，第一阵列模块，用于建立分维度编号的二维或三维的实空间的格点阵列；关联确定模块，用于确定原子与距离最近的实空间格点的关联；第二阵列模块，用于建立与衍射实验范围和精度对应的倒空间的格点阵列；分维计算模块，用于采用泰勒展开对传统算法作近似处理，分别对各维度的变量求和计算出倒空间格点的振幅；插值计算模块，用于以倒空间格点的振幅作为输入，并通过插值计算得到任意转移动量上的振幅。
工作原理：本发明通过采用泰勒展开的方式，使得基于原子位置的射线衍射计算可以按空间维度分离变量求和，每次求和只涉及一个维度，大大减小计算量的同时得到精确度按展开阶数指数上升的近似结果，可以在现代超级计算机上实现尺度可达微米以上的生物、医药等样品的精确到原子层面的射线衍射数值模拟，现代超级计算机如具有GPU-CPU异构的计算机；尤其是在多角度衍射、复色光衍射、大体系衍射等多种情形下，所需的计算代价显著小于传统的直接对每个原子和转移动量求和的计算方法；且原子数量越多，转移动量点越多，计算分析方法在减小计算代价上的益处越显著。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质（包括但不限于磁盘存储器、CD-ROM、光学存储器等）上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备（系统）、和计算机程序产品的流程图和／或方框图来描述的。应理解可由计算机程序指令实现流程图和／或方框图中的每一流程和／或方框、以及流程图和／或方框图中的流程和／或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能的步骤。
以上的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.多原子体系射线衍射的快速分析方法，其特征是，包括以下步骤：
建立分维度编号的二维或三维的实空间的格点阵列；
确定原子与距离最近的实空间格点的关联；
建立与衍射实验范围和精度对应的倒空间的格点阵列；
采用泰勒展开对传统算法作近似处理，分别对各维度的变量求和计算出倒空间格点的振幅；
以倒空间格点的振幅作为输入，并通过插值计算得到任意转移动量上的振幅。
2.根据权利要求1所述的多原子体系射线衍射的快速分析方法，其特征是，所述泰勒展开的计算公式具体为：
Figure 704159DEST_PATH_IMAGE001
其中，
Figure 981819DEST_PATH_IMAGE002
表示自然对数的底数；
Figure 682927DEST_PATH_IMAGE003
表示倒空间中格点上的转移动量在x方向上的分量；
Figure 281399DEST_PATH_IMAGE004
是任一原子位置与实空间关联格点之间的位置差在x方向上的分量；
Figure 158088DEST_PATH_IMAGE005
是单位虚数；
Figure 75491DEST_PATH_IMAGE006
是泰勒展开的阶数。
3.根据权利要求1所述的多原子体系射线衍射的快速分析方法，其特征是，所述各维度的变量求和包括第一个维度求和、第二个维度求和以及第三个维度求和；
第一个维度求和的计算公式为：
Figure 670420DEST_PATH_IMAGE007
其中，
Figure 72583DEST_PATH_IMAGE008
、
Figure 334937DEST_PATH_IMAGE009
分别表示另两个维度上的格点编号；
Figure 656197DEST_PATH_IMAGE010
表示与格点关联的不同原子的编号；
Figure 613788DEST_PATH_IMAGE011
表示倒空间的格点在第一个维度上的编号；
Figure 180161DEST_PATH_IMAGE012
表示倒空间格点上的转移动量在第一个维度上的分量；
Figure 297022DEST_PATH_IMAGE013
表示原子的散射强度；
Figure 930128DEST_PATH_IMAGE014
表示原子位置坐标在第一个维度上的分量；
Figure 765229DEST_PATH_IMAGE015
、
Figure 368249DEST_PATH_IMAGE016
分别表示原子与关联格点的位置差在第二、第三两个维度上的分量；幂指数
Figure 480561DEST_PATH_IMAGE017
和
Figure 645089DEST_PATH_IMAGE018
从零一直取到泰勒展开的阶数；
Figure 701907DEST_PATH_IMAGE019
表示第一个维度求和的计算结果；
第二个维度求和的计算公式为：
Figure 983983DEST_PATH_IMAGE020
其中，
Figure 75436DEST_PATH_IMAGE021
表示倒空间的格点在第二个维度上的编号；
Figure 315925DEST_PATH_IMAGE022
表示倒空间格点上的转移动量在第二个维度上的分量；
Figure 594459DEST_PATH_IMAGE023
表示在第二维第三维上编号分别为
Figure 318044DEST_PATH_IMAGE024
和
Figure 139370DEST_PATH_IMAGE025
的实空间格点坐标在第二个维度上的分量；
Figure 409814DEST_PATH_IMAGE026
表示第二个维度求和的计算结果；
第三个维度求和的计算公式为：
Figure 175645DEST_PATH_IMAGE027
其中，
Figure 924158DEST_PATH_IMAGE028
表示倒空间的格点在第三个维度上的编号；
Figure 960509DEST_PATH_IMAGE029
表示倒空间格点上的转移动量在第三个维度上的分量；
Figure 277221DEST_PATH_IMAGE030
表示在第三维上编号为
Figure 530348DEST_PATH_IMAGE031
的实空间格点坐标在第三个维度上的分量；
Figure 348131DEST_PATH_IMAGE032
表示第三个维度求和的计算结果。
4.根据权利要求1所述的多原子体系射线衍射的快速分析方法，其特征是，若所述实空间的格点阵列为二维，则依据泰勒展开直接进行变量求和；若所述实空间的格点阵列为三维，则可依据泰勒展开直接进行变量求和，或在第一个维度上建立的格点基础上进行变量求和。
5.根据权利要求1所述的多原子体系射线衍射的快速分析方法，其特征是，若所述实空间的格点阵列和/或倒空间的格点阵列是二维，则每个格点有两个编号；
若所述格点阵列是三维，则每个格点有三个编号；
任意一个编号不同，表示不同的格点；
格点阵列覆盖的空间范围和稀密程度由涉及的转移动量范围和计算要求的精确度决定。
6.根据权利要求1所述的多原子体系射线衍射的快速分析方法，其特征是，所述原子与距离最近的实空间格点的关联确定过程具体为：
为每个格点建立一个原子基团列表；
原子基团列表中的原子基团中存在至少一个原子在第二、三维度上的投影到相应格点的距离最近，而非离另一个格点更近；
构建一个调用函数，通过调用函数对每个格点遍历关联的原子基团中与该格点关联的原子。
7.根据权利要求1所述的多原子体系射线衍射的快速分析方法，其特征是，所述各维度的变量求和计算采用分布式计算和存储开展并行计算。
8.多原子体系射线衍射的快速分析系统，其特征是，包括：
第一阵列模块，用于建立分维度编号的二维或三维的实空间的格点阵列；
关联确定模块，用于确定原子与距离最近的实空间格点的关联；
第二阵列模块，用于建立与衍射实验范围和精度对应的倒空间的格点阵列；
分维计算模块，用于采用泰勒展开对传统算法作近似处理，分别对各维度的变量求和计算出倒空间格点的振幅；
插值计算模块，用于以倒空间格点的振幅作为输入，并通过插值计算得到任意转移动量上的振幅。
9.一种计算机终端，包含存储器、处理器及存储在存储器并可在处理器上运行的计算机程序，其特征是，所述处理器执行所述程序时实现如权利要求1-7中任意一项所述的多原子体系射线衍射的快速分析方法。
10.一种计算机可读介质，其上存储有计算机程序，其特征是，所述计算机程序被处理器执行可实现如权利要求1-7中任意一项所述的多原子体系射线衍射的快速分析方法。
CN202111535638.1A 2021-12-16 2021-12-16 多原子体系射线衍射的快速分析方法、系统、终端及介质 Active CN113936746B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111535638.1A CN113936746B (zh) 2021-12-16 2021-12-16 多原子体系射线衍射的快速分析方法、系统、终端及介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111535638.1A CN113936746B (zh) 2021-12-16 2021-12-16 多原子体系射线衍射的快速分析方法、系统、终端及介质

Publications (2)

Publication Number Publication Date
CN113936746A true CN113936746A (zh) 2022-01-14
CN113936746B CN113936746B (zh) 2022-02-22

Family

ID=79288972

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111535638.1A Active CN113936746B (zh) 2021-12-16 2021-12-16 多原子体系射线衍射的快速分析方法、系统、终端及介质

Country Status (1)

Country Link
CN (1) CN113936746B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114969029A (zh) * 2022-04-22 2022-08-30 南京航空航天大学 一种基于泰勒展开的多维插值方案及其维间点取值方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE1958234A1 (de) * 1968-12-16 1970-07-09 Ibm Mikroskopie-Geraet zur Untersuchung kristallinen Materials
US20020035455A1 (en) * 2000-09-15 2002-03-21 Xinhui Niu Generation of a library of periodic grating diffraction signals
US20100082306A1 (en) * 2008-01-02 2010-04-01 Mills Randell L System and method of computing the nature of atoms and molecules using classical physical laws
US20170138851A1 (en) * 2014-07-24 2017-05-18 Nxgen Partners Ip, Llc System and method for multi-parameter spectroscopy
WO2020084272A1 (fr) * 2018-10-26 2020-04-30 Ixblue Système et procédé d'interferomètre atomique multi-axes
WO2021108715A1 (en) * 2019-11-26 2021-06-03 Sail Sv Llc X-ray imaging system
CN113160376A (zh) * 2021-05-13 2021-07-23 中国工程物理研究院流体物理研究所 一种多色光照明三维粒子成像系统及方法
CN113447507A (zh) * 2021-06-28 2021-09-28 中国科学技术大学 一种x射线衍射信号的采集方法、系统、设备及存储介质

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE1958234A1 (de) * 1968-12-16 1970-07-09 Ibm Mikroskopie-Geraet zur Untersuchung kristallinen Materials
US20020035455A1 (en) * 2000-09-15 2002-03-21 Xinhui Niu Generation of a library of periodic grating diffraction signals
US20100082306A1 (en) * 2008-01-02 2010-04-01 Mills Randell L System and method of computing the nature of atoms and molecules using classical physical laws
US20170138851A1 (en) * 2014-07-24 2017-05-18 Nxgen Partners Ip, Llc System and method for multi-parameter spectroscopy
WO2020084272A1 (fr) * 2018-10-26 2020-04-30 Ixblue Système et procédé d'interferomètre atomique multi-axes
WO2021108715A1 (en) * 2019-11-26 2021-06-03 Sail Sv Llc X-ray imaging system
CN113160376A (zh) * 2021-05-13 2021-07-23 中国工程物理研究院流体物理研究所 一种多色光照明三维粒子成像系统及方法
CN113447507A (zh) * 2021-06-28 2021-09-28 中国科学技术大学 一种x射线衍射信号的采集方法、系统、设备及存储介质

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
BD BUTLER等: "Interpretation of displacement-caused diffuse scattering using the Taylor expansion", 《ACTA CRYSTALLOGRAPHICA SECTION A》 *
VALENTINA SHUMAKOVA等: "Comparative Study of Harmonic Generation in Air and Argon in Light Filaments Driven by Circularly Polarized Mid-IR Pulses", 《2019 CONFERENCE ON LASERS AND ELECTRO-OPTICS EUROPE & EUROPEAN QUANTUM ELECTRONICS CONFERENCE (CLEO/EUROPE-EQEC)》 *
姬洪等: "应用X射线倒易空间衍射图研究CdZnTe晶片", 《第八届全国X射线衍射学术会议论文集》 *
孙肇玲: "电子衍射斑点与倒格点关系的探讨", 《河北工学院学报》 *
方渡飞等: "较小动量转移下γ射线的相干散射", 《物理学报》 *
朱红波等: "量子微腔中原子束偏转与衍射的横向动量效应", 《中国科学(A辑)》 *
陈小辉等: "基于单晶原位X射线衍射直接测量任意冲击方向下的晶格应变", 《第十四届全国物理力学学术会议缩编文集》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114969029A (zh) * 2022-04-22 2022-08-30 南京航空航天大学 一种基于泰勒展开的多维插值方案及其维间点取值方法

Also Published As

Publication number Publication date
CN113936746B (zh) 2022-02-22

Similar Documents

Publication Publication Date Title
Zhou et al. Mmrotate: A rotated object detection benchmark using pytorch
Hand et al. Automated tracking of migrating cells in phase‐contrast video microscopy sequences using image registration
CN113936746B (zh) 多原子体系射线衍射的快速分析方法、系统、终端及介质
Chen et al. Two-Loop integrals for CP-even heavy quarkonium production and decays: Elliptic Sectors
CN109858158B (zh) 一种计算流体力学模拟的参数配置方法及系统
Huang et al. A systematic approach for online minimizing volume difference of multiple chambers in machining processes based on high-definition metrology
CN105405100B (zh) 一种稀疏驱动sar图像重建正则化参数自动选择方法
Beier et al. GiST scan acceleration using coprocessors
Deng et al. Prediction of transonic flow over supercritical airfoils using geometric-encoding and deep-learning strategies
Li et al. Research on three-dimensional reconstruction technology of line laser scanning scene based on Otsu method
Ranno et al. Neural representations of cryo-EM maps and a graph-based interpretation
Sane et al. Scalable in situ computation of Lagrangian representations via local flow maps
CN108629127B (zh) 一种周向非对称缺陷的涡流检测磁场的半解析计算方法
Huang et al. Parallel Performance and Optimization of the Lattice Boltzmann Method Software Palabos Using CUDA
Song et al. Phase demodulation of single frame projection fringe pattern based on deep learning
Patelli et al. Computational optimization strategies for the simulation of random media and components
Lobera et al. High performance computing for a 3-D optical diffraction tomographic application in fluid velocimetry
Moreno et al. HPC enables efficient 3D membrane segmentation in electron tomography
Taylor et al. Cavity volume and free energy in many-body systems
Huang et al. 3D reconstruction based on model registration using RANSAC-ICP algorithm
Peng et al. An Improved Monocular-Vision-Based Method for the Pose Measurement of the Disc Cutter Holder of Shield Machine
Dvoynishnikov et al. Phase images processing using parallel programming
Iryanto et al. Irregular grid interpolation using radial basis function for large cylindrical volume
CN108804859B (zh) 一种对称逐次局部枚举拉丁超立方试验设计方法
Qiao et al. An Efficient Spatio-Textual Skyline Query Processing Algorithm Based on Spark

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
GR01 Patent grant
GR01 Patent grant