基于分子干涉函数的光子-原子相干散射截面计算方法研究.pdf

第58 卷第3期2024年3月原子能科学技术Atomic Energy Science and TechnologyVol.58,No.3Mar.2024基于分子干涉函数的光子-原子相干散射截面计算方法研究徐宁，祖铁军*，曹良志，吴宏春（西安交通大学核科学与技术学院，陕西西安7 1 0 0 49)摘要：X射线衍射在物质结构分析和材料无损检测领域有着广泛的应用，其基本物理原理为光子与物质发生的相干散射。传统的相干散射截面计算方法基于独立原子形状因子近似方法，忽略了光子动量转移较小时与原子发生相互作用时的分子干涉效应，影响相干散射截面的计算精度。因此，为了获得光子动量转移较小时精确的相干散射截面，本文在核数据处理程序NECP-Atlas中对基于分子干涉函数的光子-原子相干散射截面计算方法进行研究，利用分子动力学模拟方法计算分子干涉函数，对蒙特卡罗程序使用的ACE格式数据库中的原子形状因子进行修正，并给出了模拟得到的水分子和乙醇分子的分子干涉函数，对基于独立原子形状因子近似方法和考虑分子干涉效应计算得到的水和乙醇的散射成像结果进行了对比分析。数值结果显示：基于分子动力学模拟得到的分子干涉函数计算得到的水的散射成像结果与文献结果吻合较好；同时，当光子动量转移较小时，分子干涉效应对相干散射的次级光子角度分布有着显著影响。本文建立的光子-原子相干散射截面计算方法可显著提高光子动量转移较小时的相干散射次级光子角度分布计算精度，可为X射线衍射模拟提供数据基础。关键词：核数据处理；光子；相干散射；原子形状因子；分子干涉效应；NECP-Atlas中图分类号：TL325doi:10.7538/yzk.2023.youxian.0472文献标志码：A文章编号：1 0 0 0-6 9 31(2 0 2 4)0 3-0 57 3-0 8Research on Calculation Method of Molecular Interference Effectof Photon-atomic Coherent ScatteringXU Ning,ZU Tiejun*,CAO Liangzhi,WU Hongchun(School of Nuclear Science and Technology,Xian Jiaotong University,Xian 710049,China)Abstract:Coherent scattering is the theoretical basis of X-ray diffraction,which iswidely used in the field of materials.A functional module to calculate the nondestructivetesting of photon-atomic and photon-molecular coherent scattering cross sections isdeveloped in the nuclear data processing code NECP-Atlas.The photon-atomic coherentscattering cross section is processed in NECP-Atlas with the same method as that ofNJOY2016.In practice,photon interacts with materials which usually consist of numer-ous polyatomic molecules.Therefore,the photon-molecular coherent scattering crosssections should be calculated instead of photon-atomic coherent scattering cross sections收稿日期：2 0 2 3-0 6-30；修回日期：2 0 2 3-0 7-2 8基金项目：国家自然科学基金（1 2 0 7 51 8 3）；国家重点研发计划（2 0 2 2 YFB1902600）*通信作者：祖铁军574to acquire more accurate simulated results.First,the independent atom model(IAM)was implemented to calculate the photon-molecular coherent scattering cross sectionswith the assumption that the atom was isolated.This approximation works well whenthe momentum transfer of incident photons is large.At small momentum transfer,theinterference effects of various atoms cannot be ignored.Then,the molecular interfer-ence model was introduced to consider the influence on angular distribution of secondaryphotons due to the interference effects.The molecular interference functions of waterand ethanol were simulated using molecular dynamics simulation,and then applied inthe ACE library.The atomic form factors of the atoms in water and ethanol moleculeswere modified by the molecular interference function.To quantify the influence ofmolecular interference effects on angular distribution of secondary photons,an imagingsystem was simulated.Numerical results show that the molecular interference effectshave a significantly influence on the angular distribution of photons at small momentumtransfer,while at large momentum transfer,the molecular interference effects can beignored.Key words:nuclear data processing;photon;coherent scattering;atomic form factor;molecular interference effect;NECP-Atlas光子在核能系统中有着重要的应用。在核反应堆释热计算 1 1 和辐射屏蔽计算 2 中需要考虑光子影响，同时，在材料无损检测和医学散射成像领域，X射线衍射也有着重要应用 3。X射线衍射的基本物理原理为光子与物质发生相干散射，为了实现对于X射线衍射过程的数值模拟，需要提供精确的光子-原子相干散射截面。光子与原子发生相互作用时，原子核外的电子在电磁场的作用下产生受迫振动，电子的受迫振动会产生一个与入射光子频率相同的交变电磁场，该交变电磁场会向周围辐射与入射光子频率相同的次级光子。由于次级光子频率相同，且具有固定的相位差，彼此之间会产生相互干涉，因此光子与原子发生的上述相互作用过程称为相干散射，也称为瑞利散射。相干散射产生的次级光子与人射光子具有相同的频率，反应前后光子的能量不会发生改变，仅光子的方向改变，因此相干散射也是弹性散射过程。传统的相干散射微分截面计算采用独立原子形状因子近似方法(4修正Thomson散射截面，该方法可精确计算光子动量转移较大时的相干散射微分截面 5。当光子动量转移较小时，采用上述假设会导致相干散射微分截面的原子能科学技术第58 卷计算存在较大偏差，无法准确描述相干散射次级光子的角度分布。因此，对于上述情况光子相干散射微分截面的计算，需要采用更加准确的计算方法。本文在传统的独立原子形状因子近似的基础上，通过考虑分子干涉效应将传统的相干散射微分截面计算方法拓展至光子动量转移较小的场景，可实现对于全能量段光子相干散射微分截面的准确计算。同时，将修正后的原子形状因子应用于蒙特卡罗程序，对比分析采用独立原子形状因子近似方法和采用分子干涉函数修正后的原子形状因子对相干散射次级光子角度分布的影响。1相干散射截面计算方法1.1独独立原子形状因子近似光子与1 个自由电子的弹性散射截面可由Thomson散射公式表示：dom=3gm(1+2)du8式中：:d 为 Thomson 微分散射截面;on 为duThomson散射截面;为平均散射角余弦。不同于Thomson散射公式描述的光子与1个自由电子的弹性散射过程，原子核外的电(1)第3期徐宁等：基于分子干涉函数的光子-原子相干散射截面计算方法研究子会受到原子核对核外电子的束缚，不同轨道上的电子能量存在差异，即原子核外的电子存在电子密度分布；同时对于除了氢元素以外的其他元素，其中性原子的原子核外都存在不止1个电子 6 。基于以上两种原因，无法直接采用Thomson散射公式对相干散射截面进行计算，需要对其进行修正以考虑上述两个因素的影响。目前，比较常用的修正方法为独立原子形状因子近似方法。对于1 个孤立原子,其原子核外电子分布不受其他原子核外电子分布影响，假设其原子核外电子密度分布为球形对称分布，此时相干散射微分截面可表示为：docob=da I(2,Z)l2du-du:dc为相干散射微分截面;f(a,z)为独式中du立原子形状因子；为光子动量转移；Z为元素质子数。动量转移可表示为：m.c2sin(0/2)2E入2hcm.csin220.607 44X 2k2式中：0 为散射角度；me为电子静止质量；h为普朗克常数；c为光速；k为以电子静止能量为单位的人射光子能量；入为人射光子波长；E为入射光子能量。采用独立原子形状因子近似修正Thom-son散射截面仅可考虑孤立原子内部的电子密度分布对相干散射微分截面的影响，忽略了分子内干涉效应和分子间干涉效应对相干散射微分截面的影响。基于独立原子形状因子近似方法计算分子的相干散射微分截面时只需要对分子中各原子的微分散射截面按照核子密度加权即可。dZf(a,Z.)dudu式中om(（E)为基于独立原子形状因子近似计算的分子相干散射截面；f（，Z.)为元素i的独立原子形状因子。该近似方法仅适用于当光子动量转移较大时的相干散射微分截面计算，当光子动量转移较小时，忽略分子内和分子间干涉效应将会对575相干散射次级光子角度分布的计算精度产生影响。1.2分子干涉函数计算分子干涉效应可分为分子内部干涉效应和分子间干涉效应。对于分子内部干涉效应可采用Debye公式计算。考虑分子内部干涉效应的原子形状因子为：Fpebye()=Zf.(,Z,)f,(,Z,)式中：FDebye（）为采用Debye公式计算的原子形状因子;f,（,Z,)为元素j的独立原子形状因子；r;为分子内元素i和i之间的距离。对于分子间干涉效应一般采用分子干涉函(2)数进行修正，分子干涉函数S()81可表示为：Zf:(a)f;(a)s;()S()=三(Zf(,z.)式中，si（)为局部分子干涉函数。在已知分子结构的情况下，局部分子干涉函数可通过分子动力学模拟进行计算 9 。s(a)=4xo.-(g(n)-1)sin(4xdr(3)(7)式中：po为原子数密度，即单位体积内的原子数;g(r)为径向分布函数，表示在1 个参考原子i周围距离为r处，原子类型为i的原子在dr厚度的球壳内的局部密度相对于参考原子数密度的比值。考虑分子内部干涉效应和分子间干涉效应的原子形状因子F()可表示为：F()=4元r采用式（8)计算得到的原子形状因子修正(4)Thomson散射截面可以得到考虑分子干涉效应的光子-原子相干散射微分截面：doThF()dudu1.3分子干涉效应在蒙特卡罗程序中的实现在蒙特卡罗程序中，相干散射微分截面 1 0可表示为：)sin(4元ari)4元r(5)（6)4元cr(8)(9)576原子能科学技术第58 卷p(u)du=元(1+)/f(,Z)/2duOCoh(E)式中：p（）为相干散射次级光子角度分布；r为电子半径。为了获得合适的概率密度分布，,将式(1 0)中对的微分转化为的微分：p(a)da2=p()d22元A(amx,Z)1+f(,Z)d2(ak）o c o h(E)2A(cmax,Z)式中：=mec2；max为光子动量转移的最大值；V2hcA,Z)为积分原子形状因子。由式(1 1)可知，在蒙特卡罗程序中进行相干散射模拟时，需要光子-原子ACE格式数据库中提供的独立原子形状因子f(Z)和积分原子形状因子A（，Z)。基于独立原子形状因子近似计算的积分原子形状因子可表示为：2f(,Z)dz?A(?,Z)=J。Z2当考虑分子干涉效应时，需要利用分子干涉函数对原子形状因子进行修正，考虑分子干涉函数的原子形状因子可表示为：F2(,Z)=S()f(,Z)将考虑分子干涉效应的原子形状因子代入式(1 0)和式(1 1)可得到考虑分子干涉效应的原子形状因子FMI(,Z)和积分原子形状因子AMI(?,Z):FMI(,Z)=VS()f(,Z)3.5a3.0F2.5人2.0F1.51.0F0.50.050S()f(,Z)dr?AMI(,Z)=JoZ(10)基于式（1 4)和式（1 5)对ACE格式数据库中需要的原子形状因子和积分原子形状因子进行计算，得到考虑分子干涉效应的光子-原子ACE格式数据库。d=2楼数值结果与验证2.1分子干涉函数计算对于结构参数已知的分子，其分子干涉函(11)数可通过分子动力学模拟程序进行计算。本文分别针对水分子和乙醇分子进行分子动力学模拟，获得各原子之间的径向分布函数,其中，水分子采用TIP4P/2005模型 1 1，乙醇分子采用LorenWilliam模型 1 2 。对于水分子和乙醇分子在模拟过程中均采用 OPLS-AA力场，利用GROMACS分子动力学程序L131分别对水分子和乙醇分子进行模拟，获得水分子和乙醇分子各原子间的径向分布函数，如图1 所示。基于图1 计算得到的水分子和乙醇分子的(12)径向分布函数，利用式(7)可计算得到水分子和乙醇分子的局部分子干涉函数，如图2 所示。最后，对局部分子干涉函数进行加权计算，可得到水分子和乙醇分子的分子干涉函数，如图3(13)所示。其中，将水分子的分子干涉函数模拟结果与实验测量结果 1 4-1 8 进行对比，可看出基于分子动力学模拟获得的水分子的分子干涉函数与实验测量结果吻合较好，可以用于修正独立原子形状因子，获得考虑分子干涉效应的原子(14)形状因子。6bH-H-O-H0-0124r/A图1 水分子(a)和乙醇分子(b)中各原子间的径向分布函数Fig.1Radial distribution function of water(a)and ethanol(b)(15)He-He-Oe-He5.Oe-Oe46632106811012012345678910r/A第3期徐宁等：基于分子干涉函数的光子-原子相干散射截面计算方法研究a1.41.21.0(as0.80.60.40.20.00.0Fig.2Partial molecular interference function of water(a)and ethanol(b)1.41.21.00.8(x)S0.60.40.20.00.0将模拟计算得到的水分子和乙醇分子的分子干涉函数分别修正，基于独立原子形状因子近似方法给出的水分子和乙醇分子中各元素的原子形状因子，可得到上述两种分子中各元素考虑分子干涉效应的原子形状因子,如图4、5所示。5772.0b-H-H-O-H.0-00.40.8X/A-1图2 水分子(a)和乙醇分子(b)的局部分子干涉函数aSPring-8SSRL-1SSRL-2APS本文计算结果0.40.8X/A-1图3水分子(a)和乙醇分子(b)的分子干涉函数Fig.3 Molecular interference function of water(a)and ethanol(b)He-He-Oe-He1.6Oe-Oe1.20.80.40.01.21.61.21.62.02.00.01.41.21.00.8F(x)S0.6F0.4F0.2F0.040.02.2分子干涉效应对相干散射角度分布的影响分析为了验证分子干涉效应对相干散射角度分布的影响，本文分别制作了采用独立原子形状因子近似和考虑分子干涉效应的光子-原子ACE格式数据库，利用蒙特卡罗程序对图6 所0.40.40.8x/A-10.81.2x/A-11.21.6b1.62.02.0a1.00.80.60.40.20.00.0图4水分子中H元素（a)和O元素(b)基于独立原子形状因子近似和考虑分子干涉效应的原子形状因子Fig.4 Independent atomic form factor approximation and atomic form factor modifiedby molecular interference effect of H(a)and O(b)elements of waterb一独立原子形状因子近似8一分子干涉效应6420.40.8x/A-1一独立原子形状因子近似一分子干涉效应1.21.62.00.00.40.8X/A-11.21.62.0578654320.00.40.81.21.6 2.0 x/A-1图5乙醇分子中C元素(a)、H 元素(b)和O元素(c)基于独立原子形状因子近似和考虑分子干涉效应的原子形状因子Fig.5 Independent atomic form factor approximation and atomic form factor modifiedby molecular interference effect of C(a),H(b)and O(c)elements of ethanol示装置进行了模拟，并对水和乙醇的成像效果进行了对比分析。图6 装置中的光子源S采用管电压为1 0 0 kV的X射线管，其光子源能谱如图7 所示，光子源距离样品中心的距离SAD为1 0 0 cm，光子源距离探测板W的距离SDD为1 53.6 cm，样品为0.5cmX0.5cmX0.5cm的立方体,探测板宽度为41 cml10。图8 示出采用独立原子形状因子近似和考虑分子干涉效应的水分子散射成像模拟结果与样品SSAD图6 成像装置示意图 1 9 Fig.6Ilustration of imaging geometryl1o原子能科学技术第58 卷a1.0一独立原子形状因子近似一分子干涉效应SDDb一独立原子形状因子近似0.8一分子干涉效应0.60.40.20.00.0 0.40.81.21.62.0X/A-1Poludniowski等 1 9 计算结果的对比。由图8可看出，基于本文计算得到的分子干涉函数计算得到的水分子成像结果与Poludniowski等的计算结果吻合较好。通过图8 的对比分析可发现，考虑分子干涉效应模拟得到的水分子散射成像在光子动量转移较小时与采用独立原子0.070.060.050.040.030.02W0.010.000.00Fig.7Spectrum of photon sourceC80一独立原子形状因子近似一分子干涉效应64200.00.40.8 1.21.62.0 x/A-10.020.04光子能量/MeV图 7 光子源能谱0.060.080.101086420-20图：采用独立原子形状因子近似(a)和考虑分子干涉效应（b)的水分子响应信号对比Fig.8 Comparison of scattering image of water using independentatom form factor approximation(a)and molecular interference effect(b)-Poludniowski等a本文模拟计算结果1-100位置/cm10864201020b20-10位置/cm-Poludniowski等-本文模拟计算结果101020第3期徐宁等：基于分子干涉函数的光子-原子相干散射截面计算方法研究形状因子近似模拟得到的结果存在明显的差异。对于光子与水分子散射成像结果可看出，考虑分子干涉效应会导致相干散射的前向散射显著减少。图9 示出本文模拟得到的乙醇分子579的散射成像模拟结果，由图9 的结果对比可得出与水分子散射成像模拟相同的结论，即考虑分子干涉效应时，相干散射的前向散射明显减少。a32b320-20图9 采用独立原子形状因子近似(a)和考虑分子干涉效应(b)的乙醇分子响应信号对比Fig.9Comparison of scattering image of ethanol using independentatom form factor approximation(a)and molecular interference effect(b)3结论本文介绍了基于独立原子形状因子近似的传统相干散射截面计算方法，同时为了实现对于X射线衍射的数值模拟，本文介绍了考虑分子干涉效应的相干散射截面计算方法。基于分子动力学模拟获得了水分子和乙醇分子的分子干涉函数,对光子-原子ACE格式数据库中相应数据进行了修正，并利用蒙特卡罗程序对水分子和乙醇分子的散射成像效果进行了数值模拟。数值结果显示，对于水分子和乙醇分子，考虑分子干涉效应会导致光子的前向散射显著减少。参考文献：1 祖铁军，徐宁，尹文，等核数据处理软件NECP-Atlas中的光子相关数据计算方法研究 J原子能科学技术，2 0 2 2,56（5）：9 6 9-9 7 7.ZU Tiejun,XU Ning,YIN Wen,et al.Researchon generation method for photon-related data innuclear data processing code NECP-AtlasJJ.Atomic Energy Science and Technology,2022,56(5):969-977(in Chinese).2祖铁军，徐宁，曹良志，等.NECP-Atlas中屏蔽数据库制作模块的开发与验证J核动力工程，2022,43(3):15-20.ZU Tiejun,XU Ning,CAO Liangzhi,et al.1-100位置/cm010201-20-10位置/cmDevelopment and verification of shielding data-base generation module in NECP-AtlasJJ.Nu-clear Power Engineering，2 0 2 2,43(3):1 5-2 0(inChinese).3HUBBELL J H.Review of photon interactioncross section data in the medical and biologicalcontextJI.Physics in Medicine and Biology,1999,44:1-22.4DERMOTT E C,HUBBELL J H,KISSEL L.EPDL97:The evaluated photon data libraryR.America:Lawrence Livermore NationalLaboratory,1997.5SOPER A K.Joint structure refinement of X-rayand neutron diffraction data on disordered materi-als:Application to liquid waterJJ.Journal ofPhysics:Condensed Matter,2007,19:335206.6SKINNER L B,HUANG C,SCHLESINGERD,et al.Benchmark oxygen-oxygen pair-distri-bution function of ambient water from X-ray dif-fraction measurements with a wide Q-rangeJ.The Journal of Chemical Physics,2013,138(7):074506.7BOKE A.The effect of molecular interference oncoherent scatteringJJ.Journal of Balikesir Uni-versity Institute of Science and Technology,2017,19(2)：1 2 3-1 36.8WANG Jiahu,AWADH N,VEDENE H.Chemical binding and electron correlation effects01020580in X-ray and high energy electron scatteringJ.The Journal of Chemical Physics,1994,101(6):4 842-4 854.9SORENSON J M,HURA G,GLAESER R M,et al.What can X-ray scattering tell us about theradial distribution functions of waterJ.TheJournal of Chemical Physics,2000,113（2 0):9149-9161.10 ROMANOP K,HORELIK N E,HERMAN BR,et al.OpenMC:A state-of-the-art MonteCarlo code for research and developmentJ.Annals of Nuclear Energy,2015,82:90-97.11 ABASCAL J L,VEGA C.A general purposemodel for the condensed phases of water:TIP4P/2005J.The Journal of Chemical Phys-ics,2005,123(23)：2 3450 5.12 ZHANG Cuijuan,YANG Xiaoning.Moleculardynamics simulation of ethanol/water mixturesfor structure and diffusion propertiesJ.FluidPhaseEquilibria,2005,231(1):1-10.13 van der SPOEL D,LINDAHL E,HESS B,etal.GROMACS:Fast,flexible,and freeJ.Journal of Computational Chemistry,2005,26(16):1 701-1 718.14 SKINNER L B,BENMORE C J,PARISE J B.Molecular arrangement in water:Random but not原子能科学技术第58 卷quiteJ.Journal of Physics:Condensed Matter,2012,24338001.15 SKINNERLB,BENMORECJ,SHYAMB,etal.Structure of the floating water bridge andwater in an electric fieldJ.Proceedings of theNational Academy of Sciences，2 0 1 2，1 0 9:16463.16 HART R T,BENMORE C J,NEUEFEIND JC,et al.Temperature dependence of isotopicquantum effects in waterJ.Physical ReviewLetters,2005,94:047801.17 HUANG C,WIKFELDT K T,NORDLUND D,et al.Wide-angle X-ray diffraction and moleculardynamics study of medium-range order in ambientand hot waterJJ.Physical Chemistry ChemistryPhysical,2011,13:19997-20 007.18 FU L,BIENENSTOCK A,BRENNAN S.X-ray study of the structure of liquid waterJ.The Journal of Chemical Physics,2009,131:234702.19 POLUDNIOWSKI G,EVANS P M,WEBB S.Rayleigh scatter in kilovoltage X-ray imaging:Isthe independent atom approximation good enoughJJ.Physics in Medicine&Biology,2009,54:6931.