基于蒙特卡罗模拟的透射式XRD背景的定量分析.pdf

第5 8 卷第2 期2024年2 月原子能科学技术Atomic Energy Science and TechnologyVol.58,No.2Feb.2024基于蒙特卡罗模拟的透射式XRD背景的定量分析袁靖茜，黄宁*，何泽，彭博，王鹏（四川大学原子核科学技术研究所辐射物理及技术教育部重点实验室，四川成都6 10 0 6 4）摘要：XRD背景主要来源于空气、装置和样品中光子的散射，XRD分析需准确扣除背景，常用的扣背景方法具有一定的局限。蒙特卡罗模拟可实现光子的甄别与统计，因此可将其用于XRD背景来源的定量分析。相干散射的模拟中需考虑物质的Form Factor,这可由Baro等的拟合式、独立原子模型或德拜公式给出；在2 0 keV能量及透射几何模式下，通过对骨、尼龙、水等物质的XRD模拟，对探测器中光子来源进行了区分和统计。结果表明，德拜公式计算的空气背景大于独立原子模型；空气散射在小角度时背景占比大，随散射角的增加从10 0%减小至40%以下，而样品中康普顿散射的计数占比则从0 增大至60%以上。该方法可对不同样品组成以及复杂装置结构的XRD背景进行定量分析。关键词：蒙特卡罗模拟；XRD背景；FormFactor;Geant4中图分类号：TL81doi:10.7538/yzk.2023.youxian.0270Quantitative Analysis of Transmission XRD BackgroundYUAN Jingxi,HUANG Ning,HE Ze,PENG Bo,WANG Peng(Key Laboratory of Radiation Physics and Technology of Ministry of Education,Institute of Nuclear Science and Technology,Sichuan University,Chengdu 610064,China)Abstract:The first prerequisite for obtaining correct information from an X-ray diffrac-tion pattern is to accurately remove the background.The most commonly used methodin experiments is direct subtraction,which ignores the absorption of the sample andother interactions within the sample such as the photoelectric effect and the counting ofphotons produced by Compton scattering.The“rolling bal algorithm and the use ofsoftware such as Jade deduction background as well as other deduction methods alsohave their own limitations.Monte Carlo simulation can identify and count the particles,so it is very important to use it for quantitative analysis of XRD background.In thisstudy,the atomic Form Factor of each element in the air was calculated by the Barofitting formula and the molecular Form Factor was given by the Debye formula and文献标志码：ABased on Monte Carlo Simulation文章编号：10 0 0-6 9 31(2 0 2 4)0 2-0 441-10收稿日期：2 0 2 3-0 4-2 8；修回日期：2 0 2 3-0 6-2 8基金项目：四川省重大科技专项（2 0 2 0 ZDZX0004）*通信作者：黄宁442compared with the calculation results of the independent atomic model.Geant4 was usedto simulate XRD for crystalline substances such as bone and nylon and amorphoussubstances such as water,and the background and real signals were effectively separatedin the results,then the XRD background was quantitatively analyzed on this basis.Theresults show that there is a huge difference in counting between the independent atomicmodel and Debyes formula,and the Debye method produces more scattered particles.XRD background comes from the air scattering and other types of interactions betweenphotons and matter in the sample,the former contributing the most at small angles,theproportion decreases from 100%to less than 40%with the increase of the scatteringangle,and the proportion of the latter increases from O to more than 60%with the scat-tering angle.The results indirectly illustrate the inaccuracy of the background subtrac-tion methods in the experiment.Such as the direct subtraction method not only ignoresthe absorption of X-ray by the sample,but also ignores the background count caused byother interactions of X-ray photons in the sample.To accurately remove the back-ground,all the counts of background mentioned above must be considered.The MonteCarlo method was used to simulate the XRD background,and the simulation was closelyconnected with the practice.The same XRD background pattern can be obtained byconstructing the same geometric model as the experiment,setting the same physicalparameters as the experiment,and using the same sample Form Factor as the experi-ment.Applying the simulation results to actual experiments may become a moreaccurate method of background deduction.Key words:Monte Carlo simulation;XRD background;Form Factor;Geant4XRD(X-ray diffraction)技术是一种不可或缺的结构测试手段，广泛应用于物相分析、应力测量及结晶度的测定等应用中1-3。在进行XRD分析时，正确扣除背景是得到准确XRD图谱的必要前提4。实际中使用最广泛的扣背景方法是直接减除法，它在没有样品的条件下进行XRD测量，获得的图谱即为背景图谱。在透射测量几何中，由于没有样品，直接减除法忽略了样品的吸收以及样品中物质与光子的非相干散射所产生的背景，因此不适用于较大或较厚样品的测量。此外，Sternberg提出的“滚球算法”适合晶体材料，其原理是保留衍射图谱中较大的梯度变化而扣除较小的梯度变化，而当样品没有尖锐的衍射峰时衍射信号会被当做背景扣除，因此不适用于非晶材料的背景处理5 。Zheng等6 结合实验提出了一种新颖的方法，通过改变样品的位置，得到了较直接减除法更准确的背景，但该方法需实验测量。此外，还有背景扣除方法如小波分析和SNIP算法，但这两种方法均以纯粹信号的角度进行背景扣原子能科学技术第5 8 卷除,未考虑物理过程7-8 。软X射线光子穿过物质时，主要发生光电效应、康普顿散射和相干散射。相干散射是指人射光子与原子核外电子发生弹性碰撞，入射方向变化而能量不变的散射。具有相同能量的相干散射光子之间可能发生相互干涉，这种相互干涉可发生在原子内不同电子之间、分子内不同原子之间、不同分子与分子之间。而具有特定结构的物质如晶体，它们的空间结构长程有序排列，由其散射的光子之间具有固定的相位差从而会产生大规模的相互干涉，这种大规模的相干散射在宏观上表现出衍射现象 。蒙特卡罗模拟可对物理过程进行仿真，并实现光子的甄别与统计，因此，本文使用该方法对XRD进行模拟，以实现背景中各计数成分及占比的定量分析。一般,XRD图谱中的背景来源于空气、实验装置对X射线的散射和样品中物质与光子的非相干散射，以及X射线荧光、致辐射等，其中空气散射对背景的贡献较大，样品中物质与光子的非相干散射以康普顿第2 期散射为主5 。XRD图谱的衍射信号来源于物质中不同密度分布的电子对X射线的大规模的相干散射L9,该相干散射通常用Form Factor加以描述，是原子或分子相干散射截面与单电子相干散射截面的比值。使用精准的Form Factor值与相干散射物理过程，是实现XRD模拟的关键。基于蒙特卡罗方法的粒子输运模拟软件有Geant4、EG S4等。Geant4是由欧洲核子研究中心（CERN)开发的一款模拟软件，具有广泛的应用领域，包括辐射防护、核反应堆设计、热核聚变研究、宇宙射线探测器设计等EGS4是一种电子和光子输运的蒙特卡罗模拟软件，由斯坦福大学等开发，主要应用于医学物理、核医学和辐射防护等领域10 11。多年来，众多学者基于 Geant4 与 EGS4 进行了诸多的相干散射模拟研究。19 9 6 年，Leliveld等9 1实现了液体水、尼龙等样品的模拟，并详细阐述了蒙特卡罗模拟相干散射的原理与要点。与此同时,Peplow、Tartari等12-13 针对生物组织（包括脂肪、骨等）的 Form Factor 进行了一系列研究,基于 EGS4实现了相干散射模拟，并更新了生物组织的Form Factor。2 0 2 0 年,Paterno 等14 基于 Geant4实现了相干散射模拟，对之前关于生物组织模拟研究进行了总结，并更新了Geant4相干散射的数据库。本文通过Baro等的拟合式与德拜公式计算空气Form Factor，代替Geant4中以独立原子模型方法计算空气FormFactor。模拟中调用Paterno的Form Factor数据库，以得到较准确的XRD图谱15 。通过模拟具有的可轻松实现和调整几何模型、可考虑和区分光子与物质的各种相互作用、可对探测器上的光子来源进行区分和记录的优势，实现模拟场景下XRD背景的定量分析。1理论与方法1.1相干散射一个人射光子与原子的相干散射截面可表示为：0)-/F(g,2 dgado oh(0)d5do Th(0)2=r:1+cos*0d52袁靖茜等：基于蒙特卡罗模拟的透射式XRD背景的定量分析443其中:doTh（)/d2为非偏振光的汤姆逊散射截面，表示1个光子与1个电子发生相干散射的概率；re为电子经典半径；F（q,Z)为该原子的Form Factor,是动量转移q和原子序数Z 的函数。q可表示为：2E.qsin(0/2)C其中：E为人射光子能量；c为真空中光速；0 为人射光子的散射角。在很多文献中，也用=q/2h表示FormFactor,单位为A-1,h为普朗克常量。因此1个原子的相干散射截面可表示为：Ccoh(E)=ra(1+cos?0)|F(q,Z)sin 0do(4)1）原子 Form Factor 的计算根据式(5)可知,计算原子相干散射截面的关键在于如何获取准确的原子FormFactor。为此,Hubbell 等16-17 做了大量的研究，通过不同的原子模型，结合实验获取了原子序数11oo的原子FormFactor，并将这些数据制为表格。截至目前，Hubbell 的原子Form Factor数据仍使用广泛。为更高效使用这些数据，许多人对其进行插值拟合，其中，Bar6等18 的拟合式与原数据拟合结果较好，将各元素FormFactor的计算公式拟合如下：F(q,Z)=r(z,Z)=zlta+a2+a2(1+a4+asat)2max(f(,Z),Fk(q,Z)if Z 10 and f(,Z)2其中：a1va2va3va4vas为各元素的拟合参数，可通过查表获得；Fk（q Z)可表示为：Fk(q,Z)=sin(2yaretan Q)YQ(1+Q)其中,Q=q/2meca,a=(Z-5/16),=1/137,=V1-,q=（2 E/c)s in（0/2)=mec/20.6074,=20.6074q/(mec)，m e 与c分别为电子静止质量与光速。利用Bar等的拟合式（式（5）、（6)），可在保证准确性的前提下，较快捷获取各元素(1)d5(2)(3)(5)(（6)Form Factor，同时，拟合式的计算使得Hub-bell的表格数据变得连续，更方便FormFactor的使用与研究。4442）分子FormFactor的计算对于分子FormFactor的计算，实验上通常的做法是先利用独立原子模型（IAM)计算FiAM：将分子内各原子FormFactor按质量权重累加而得到分子FormFactor。再用FiAM乘以干涉函数，以修正独立原子模型计算过程对干涉效应的忽略，具体的计算过程可表示为：(7)(FM(q)=FiAM(q)s(q)其中:s(q)为干涉函数;FM（q)为经干涉函数修正后的分子FormFactor。在这种表示方法下，干涉函数是表征干涉效应的关键信息，它通常是一个震荡函数，并且一般从实验测量的衍射图谱中提取获得。也可用德拜公式对分子的Form Factor进行计算，如Narten就使用该公式计算了水分子的Form Factor，且与实验结果对比吻合较好19 。相较于独立原子模型对各原子FormFactor的质量权重累加，德拜公式考虑了分子内原子之间的相互干涉，即分子内各原子FormFactor应是一种“振幅叠加关系”而不是“强度叠加关系”。德拜公式可表示为：IN=22FF.NN,sin2元srk2元srk其中：F,、Fk 为原子Form Factor;s为原子间的散射矢量；r丞为两原子间的原子间距。1.2空气FormFactor的计算基于1.1节中阐述的计算方法可对空气的FormFactor进行计算。结合式（5）、（6）可得到空气中各元素的原子FormFactor，图1示出了空气中部分元素的计算结果，其中横坐标=q/2h=sin(0/2)/入，用以替代q表示物质的Form Factor。空气是气体，各分子间具有较大的距离而不会发生分子与分子之间的相互干涉。若用独立原子模型十干涉函数修正的方法表示空气的FormFactor,则干涉函数s(q)=19。显然，此时的干涉函数已失去了修正的作用，因此使用独立原子模型单独表示空气FormFactor不准确。独立原子模型的不足也被很多研究者提及9 141。相较而言，使用考量了分子内原子间的相互干涉的德拜公式计算空气Form Factor原子能科学技术第5 8 卷则更准确。本文考虑的空气成分列于表1，此外，表格里还包含了计算FormFactor所需要的原子间距与质量占比2 0 1815129630E0.0图1空气中部分原子FormFactorFig.1Partial atomic Form Factor in air表1空气成分Table1Aircomposition原子成分间距/pmN2TN-N=145.002ro-0=120.9CO2rc-0=116.0,ro-0=232.0(8)杂质（气态水为代表）稀有气体（Ar为代表）需注意，表1中考虑的空气杂质和稀有气体以质量占比最高的气态水和Ar为代表，其中Ar为原子气体，故没有分子内原子间距这一参数。因此，利用表中参数与德拜公式，可计算出空气内各分子Form Factor，另外，将计算结果与独立原子模型进行比较，结果如图2 所示。从图2 可看出，忽略了分子内原子之间相互干涉的独立原子模型下各分子Form Factor与德拜公式的计算结果差距很大，这在较小时尤为明显，也与客观事实相符：干涉效应主要发生在较小时17 。如上所言，空气FormFactor不需考虑分子与分子之间的相互干涉，因此可把图2 中各分子FormFactor根据质量权重累加，得到空气FormFactor，如图3所示。通过图3可明显CH0.30.6X/A-1TH-0=97.6,rHH=155.9一0.91.2体积质量含量/%含量/%78.090.75529820.900.231 0550.040.0006080.030.0001860.940.012.8511.5第2 期看出德拜公式与独立原子模型的计算结果差异,图中结果显示，aFiAM，是由于德拜公式中分子内原子之间的干涉相加;0.18 A-1a0.4 A-1时,FialyeFiAM，是分子内原子之间的干涉相消所致，这说明在不需考虑分子与分子之间相互干涉的情况下,德拜公式计算的空气 Form Factor考虑了分子内的所有干涉,拥有准确的物理意义。5005Fgebye(CO2)Figebye(O)400Fgebye(N2)300Figebye(H,0)一-FAM(CO2)200-FiAM(O2)-FiAM(N2)-FiAM(H,O)10000.0FDebye由德拜公式计算的分子FormFactor；FIAM由独立原子模型计算的分子FormFactor图2 由德拜公式与独立原子模型计算的空气中各分子FormFactorFig.2Form Factor of each molecule in air calculatedby Debye formula and independent atomic model250200A150F10050F0.0图3由德拜公式与独立原子模型计算的空气Form FactorFig.3AirForm Factor calculated by Debyeformula and independent atomic model1.3模拟方法Paterno等14 更新了Geant4 中相干散射多种物质的Form Factor数据库，参考其部分代码，本研究实现了空气Form Factor的计算与多种样品的XRD背景计数探究。本次模拟袁靖茜等：基于蒙特卡罗模拟的透射式XRD背景的定量分析源Z图4Geant4模拟的设备及几何Fig.4Equipment and geometry in Geant4 simulation0.30.6x/A-10.30.6X/A-1445的设备及其几何如图4所示。模拟中整个空间填充材料为空气；源为单能2 0 keV的点源，光子发射方向为之轴，发散角为0.7；样品(本文也称为Phantom)中心在（O,O，d i)处，它是底面半径为2 cm、厚度为1 mm的圆柱体，这里的填充材料可被任意指定；探测器中心在（0，0，di十d)处是一半径10 cm的理想探测器,厚度为0.2 cm。diPhantom0.9模拟分为实验1和实验2。实验1中Phantom填充材料为空气，相当于直接减除法对空气背景的直接测量。模拟实验1的目的是对比独立原子模型与德拜公式计算的FormFactor对空气背景计数的影响。在Geant4中，允许用户以输入文件的方式指定样品材料的FormFactor，如果用户不提供这个文件，Geant4会按独立原子模型计算样品材料的Form Factor。因此，当没有传入 Form FactorFFRAMFDebye0.9d数据文件时，可得到独立原子模型计算的Form Factor的模拟结果。模拟实验2 将Phantom换成了液体水、尼龙等不同的物质，这些物质被许多研究者用于相干散射模拟，它们的Form Factor来自实验测量，数据较准确，可获得较好模拟结果，因此更方便于背景的定量分析。在模拟中，背景信号与样品衍射信号将被区分，这是蒙特卡罗模拟的优势所在：可对计数光子进行甄别，即根据每个光子在输运过程中发生的不同种相互作用类型或发生作用的位置对光子加以区分，因此可实现对来自空气与样品的计数光子分别统计。2结果与讨论2.1林模拟结果的对比验证在进行模拟实验1、2 的结果展示前，以液探测器446体水为例，将关于XRD的模拟结果与Paterno等的工作、Chaparian、K o s a n e t z k y 等2 1-2 3 的实验结果进行了对比验证，与Paterno程序的对比如图5 所示。图5 a是相同模拟条件下本团队的工作与Paterno程序模拟结果的对比，即有无空气MI效应的总计数对比，图5 b来源于图5 a中样品衍射信号的提取。本团队的工作是在Paterno等工作的基础上增加了对空气MI效应的描述。从图5 a可发现，空气的MI效应会对低角a4000300062000100000图6 中不同颜色的曲线分别代表本团队的模拟工作与Kosanetzky、Ch a p a r i a n 等的实验结果。3条曲线整体存在一定差异，受实验中不同几何、仪器等多种因素影响。本团队的模拟与Kosanetzky的实验结果相比,两曲线在小角度时基本吻合，大角度出现较小差异，而与Chaparian的结果相比，小角度时存在差异，大角度基本吻1.00.8F0.60.20.00.0图6 Geant4模拟与实验结果的对比Fig.6Comparison of Geant4with experimental results原子能科学技术第5 8 卷度的总计数造成显著影响，但在角度增大时，MI效应的影响持续降低，这是由于空气的MI效应主要发生在小角度，随着角度的增大、MI截面的减少,Geant4中独立原子模型的计算结果越趋近于实际的相干散射截面。图5 b中，两曲线基本重合略有差异，曲线的微小差异说明了不同的空气相干散射截面会造成样品衍射信号的不同，但影响不大。总计数在低角度上的差异，衍射信号的基本一致，均符合预期，说明了空气MI效应会影响低角度的背景计数。一本工作b1.0二Paterno程序结果0.80.6F0.4H0.20.051510散射角/)a有无空气MI效应的总计数对比;b归一化的衍射信号对比图5 XRD与Paterno程序的对比结果Fig.5 Comparison of XRD with Paterno procedureGeant4模拟-Kosanetzky的实验测量Chaparian的实验测量工0.10.2x/A-1一本工作Paterno程序结果15200.30.4250合。因此整体看来，模拟结果介于两者之间。2.2XRD背景的定量分析验证模拟方法后，进行了模拟实验1与模拟实验2。模拟实验1的空气背景计数I的结果如图7 所示。可看到，德拜公式的计数FormFactor明显多于独立原子模型FormFactor计数，在低角度0 5 时，差异在1倍左右；同时，德拜公式计算的空气背景曲线与实验更相符：即低角度时，空气背景产生大量计数，随着角度的增大，背景计数降低。因此使用德拜公式计算的Form Factor进行空气背景的探究，使得空气散射产生的背景拥有更准确的计数。模拟实验2 是对不同样品的XRD背景探究。首先，各样品的XRD图谱模拟结果如图8所示。同组数据中的I、I b、I。均来自同一模拟实验,其中，I为总计数，I为背景计数，I。为样品衍射信号计数，且满足I.=I十I。此外，对每次模拟的样品衍射信号计数I，取相对强度，可看到这些计数曲线均具有较明显的衍射特征峰；从I的曲线规律中可看到，背景在小510散射角/）152025第2 期104103102100图7 不同FormFactor计算方法下空气背景计数对比Fig.7Comparison of air background countsin different Form Factor calculation methods袁靖茜等：基于蒙特卡罗模拟的透射式XRD背景的定量分析510散射角/)447角度时几乎与总计数相等，因此，可认为总计数-lpebye-IAM1520中小角度的大部分计数来自于背景，但随着角度的增加，背景计数减小。在获取了各样品的XRD图谱的前提下，对模拟2 中不同样品的背景曲线进行了对比，即I对比图，结果如图9 所示。结果显示，除样品为骨外，各实验的背景曲线在散射角小于2510时大致重合，这是由于不同实验中样品均是底面积2 cm、厚度1mm的圆柱体，而除骨外的所有样品基本上均是轻元素，对X射线光子的吸收都很小，因此在散射角小于10 时也与400030002.000100040003000数20001000W0500040003000志2.00010000540003.00020001000011液体水的计数曲线510散射角/)骨的计数曲线510散射角/）NH,NO,的计数曲线510散射角/)尼龙的计数曲线510散射角/)图8各样品的XRD模拟Fig.8XRD simulation of each sample1.00.80.6F0.40.20.0152015201520Is1520液体水的衍射信号工2501.00.80.60.40.20.052501.00.80.6F0.40.2250.0%1.0F0.80.6F0.40.20.0525055NH,NO,的衍射信号510散射角/)尼龙的衍射信号510散射角/）10散射角/(）骨的衍射信号1015散射角/()1520152015202025252525448“无样品的计数曲线基本重合。而样品为骨的计数曲线明显低于其他计数曲线，原因是骨中含有Ca，对X射线光子的吸收明显大于其他样品，这说明了含有重元素的样品的吸收对背景会产生较大影响。而在散射角大于10 后，各样品的背景曲线展现了较大差异，且明显多于“无样品”的计数曲线，这是由于在这个散射角区域，X射线光子与样品的非相干散射相互作用(以康普顿散射为主，其他可忽略不计)会产生大量的计数光子。对总计数中各类型光子计数进行对比探究，结果如图10 所示，其中实验1041031021010图9 不同样品的背景计数Fig.9Background counting of different samples原子能科学技术第5 8 卷样品为液体水和骨。从图10 a、c 可看到，总计数I由3种光子计数构成：样品衍射信号计数I空气背景计数Iair和样品中光子与物质的非相干散射相互作用所引起的计数Iother。图10 b、d 的结果表明,Iai在整个散射角范围内占比很大，特别是散射角小于5 时，占比接近1,虽然其随角度的增大而剧烈降低，但散射角最大时，依旧贡献了30%40%的占比；Iother与样品的体积、厚度和各反应截面有关，它在散射角较小时贡献了较低的计数占比，但在散射角持续增大时，占比超过空气背景计数的占比，最终达到6 0%以上，因此在本模拟的XRD背景构成中，样品中光NHANO;样品尼龙样品PMMA样品水样品骨样品无样品M15子与物质的非相干散射相互作用所引起的计数是不可忽略的。图9、10 是对背景的定量分析结果，结果展示了空气散射、样品的吸收和光子与物质的非相干散射相互作用对背景的影响，这也间接说明了实验中背景扣除方法的不准确来源，如直1015散射角/)2025接减除法不止忽略了样品对X射线的吸收，同时也忽略了X射线光子在样品中发生的非相干散射相互作用而引起的背景计数。要精确去除背景，就必须同时考虑3个因素。使用蒙特200001500010000500000200001500010000500000alairIother液体水510散射角/)lairIother骨510散射角/）图10 总计数中不同类型计数对比Fig.1o Comparisons of different counts in total countb1.0-lotherlair0.80.6F0.4F0.2F0.01520C1520液体水12501.00.8F0.6F0.4F0.20.052505510散射角/（）-lotherlair骨1015散射角/）15202025d25第2 期卡罗的方法模拟XRD背景，可将模拟与实际紧密联系，构建与实验相同的几何模型、设置与实验相同的物理参数、使用与实验相同的样品FormFactor就可得到与实验结果相近的XRD背景图谱，而将模拟结果用于实际实验，可能会成为一种更准确的背景扣除方法。3结论使用Geant4模拟了独立原子模型和德拜公式计算FormFactor时的空气背景和几种样品的XRD图谱，以及XRD背景的来源及占比。结果表明，德拜公式考虑了分子内的相互干涉，得到的空气散射背景更准确；液体水的模拟结果与Paterno、K o s a n e t z k y、C h a p a r ia n 等工作的对比有效验证了本文模拟方法的准确性；分离的样品衍射信号均具有较明显的衍射特征峰；样品的吸收会影响整个背景曲线的计数；空气背景在小角度时贡献巨大，随着角度的增加，占比逐渐减小，而此时来自样品中非相干散射的贡献不可忽略。本研究使用蒙特卡罗方法实现了对XRD背景的定量分析，可考虑各种物理过程的贡献，得到背景的来源与占比。针对不同的样品组成以及复杂的装置结构，可采用同样的方法进行模拟，是一种方便、高效、精确度高的XRD背景分析方法。参考文献：1 BISH D,BLAKE D,VANIMAN D,et al.Thefirst X-ray diffraction measurements on MarsJJ.IUCRJ,2014,1(6):514-522.2BORTOLOTTI M,LUTTEROTTI L,BOROVINE,et al.Combined XRD-XRF cluster analysisfor automatic chemical and crystallographic sur-face mappingsJJ.Powder Diffraction,2019,34(S1):36-41.3 DALEY T,OPUNI K B,RAJ E,et al.Monito-ring the process of formation of ZnO from ZnO2using in situ combined XRD/XAS techniqueJ.Journal of Physics Condensed Matter,2021,doi:10.1088/1361-648X/abfb91.4KRIEGER M,CHAMBERSJ L,CHRISTOP GG,et al.Data-collection in protein crystallogra-phy-capllary effects and background corrections袁靖茜等：基于蒙特卡罗模拟的透射式XRD背景的定量分析2009.6ZHENG Y,SPELLER R,GRIFFITHS J.Anovel method to remove the background fromX-ray diffraction signalJI.Physics in Medicine&.Biology,2018,doi:10.1088/1361-6560/aaac9e.7方勇，曾立波，雷俊锋，等。一种新的X射线能谱背景扣除方法J.分析测试学报，2 0 0 1（3）：2 3-27.FANG Yong,ZENG Libo,LEI Junfeng,et al.A new method of background subtraction forenergy dispersive X-ray spectraJ.Journal ofInstrumental Analysis,2001(3):23-27(in Chi-nese).8文刘远.基于MATLAB藏药X射线衍射图谱的数据处理D.青海：青海师范大学，2 0 15.9LELIVELD C J,MAAS J G,BOM V R,et al.Monte Carlo modeling of coherent scattering:Influence of interferenceJ.IEEE TransactionsonNuclearScience,1996,43(6):3315-3321(inChinese).10邱睿，李君利，武祯，等四种蒙特卡罗程序的比较计算J原子能科学技术，2 0 0 8，42（12）：1 149-1 152.QIU Rui,LI Junli,WU Zhen,et al.Comparisoncalculation of four Monte-Carlo codes J.AtomicEnergy Science and Technology,2008,42(12):1 149-1 152(in Chinese).11 GUATELLI S,CUTAJAR D,OBORN B,et al.Introduction to the Geant4 Simulation toolkitJ.Concepts and Trends in Medical RadiationDosimetry,2011:303-322.12PEPLOW D E，V ERG H ESE K.M e a s u r e dmolecular coherent scattering form factors ofanimal tissues,plastics and human breast tissueJ.Physics inMedicine&Biology,1998,43(9):2431-2452.13 TARTARI A,TAIBI A,BONIFAZZI C,et al.Updating of form factor tabulations for coherentscattering of photons in tissuesJ.Physics inMedicine&.Biology,2002,47(1):163-175.14 PATERNO G,CARDARELLI P,GAMBACCINIM,et al.Comprehensive data set to includeint