CN108593695A - 从测量深度剖析数据直接获取薄膜元素成分分布的方法 - Google Patents
从测量深度剖析数据直接获取薄膜元素成分分布的方法 Download PDFInfo
- Publication number
- CN108593695A CN108593695A CN201810317643.7A CN201810317643A CN108593695A CN 108593695 A CN108593695 A CN 108593695A CN 201810317643 A CN201810317643 A CN 201810317643A CN 108593695 A CN108593695 A CN 108593695A
- Authority
- CN
- China
- Prior art keywords
- depth
- fathoms
- profile data
- film element
- directly acquiring
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
- G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明涉及一种从测量深度剖析数据直接获取薄膜元素成分分布的方法,主要包括:(1)利用深度剖析技术对薄膜样品中元素的分布进行测量,并将测量得到的强度‑溅射时间谱转换成浓度‑深度谱;(2)MRI参数的确定:包括原子混合长度、信息深度、样品粗糙度三个参数的确定;(3)采用正则化方法中的Tikhonov‑TV方法,并结合其牛顿阻尼法进行求解,获取原始成分浓度‑深度分布。本发明利用MRI模型和反卷积方法,来获得原始成分分布。针对反卷积过程中,低信噪比和实验数据不足等问题,采用正则化方法中的Tikhonov‑TV方法,并结合其牛顿阻尼法,可以对噪声有着较好的容纳能力,对于数据点较少或不足的情况也具有较优的模拟效果,可以准确地获得原始成分浓度‑深度分布。
Description
技术领域
本发明属于表面分析技术领域,尤其涉及一种从测量深度剖析数据直接获取薄膜元素成分分布的方法。
背景技术
深度剖析定量分析,一方面能更好的了解深度剖析技术的缺陷,提供定量标准与解释,另一方面,通过直接卷积或反卷积,能获得更为真实的成分深度分布。目前,卷积理论在深度剖析定量分析上占主导地位;而反卷积理论的优势则体现在提升测量深度剖面质量上。在深度剖析中,直接卷积的方法只有假设原始浓度分布,才能通过调整分辨率函数的相应参数(如原子混合度,粗糙度等)来达到最佳的拟合目的。而影响反卷积理论应用的原因是,直接反卷积方法获得样品的深度浓度分布的过程很复杂,同时由于数据精度不足,或是低信噪比,将对反卷积结果造成极大的误差。
目前Gautier等人通过SIMS获得了在硅衬底上硅硼掺杂薄膜的实验数据,并利用由Dowsett提出的响应函数模型[4]和迭代反卷积的方法对实验谱进行了定量分析,拟合得到的纳米层结构厚度结果与原始结构较为接近。但是由算法本身限制迭代反卷积方法对原始结构为平滑函数较为适合,当原始结构为方波时,该算法就不再适用了。
发明内容
针对现有技术中存在的问题,本发明提供一种从测量深度剖析数据直接获取薄膜元素成分分布的方法。
为了实现上述的目的,采用如下的技术方案:
一种从测量深度剖析数据直接获取薄膜元素成分分布的方法,主要包括以下步骤:
(1)利用深度剖析技术对薄膜样品进行深度剖析测量,并将测量得到的强度-溅射时间谱转换成浓度-深度谱;
(2)MRI参数的确定:包括原子混合长度、信息深度、样品粗糙度三个参数的确定;
(3)采用正则化方法中的Tikhonov-TV方法,并结合其牛顿阻尼法进行求解,获取原始成分深度浓度分布。
卷积与反卷积的对比:卷积是在已知深度分辨率函数的条件下,假定原始的层结构分布,通过迭代方法模拟计算的深度谱,并与测量的深度谱进行比较直至达到最佳拟合。而反卷积方法是在已知深度分辨率函数的条件下,通过反卷积方法直接获取原始的成分分布。
Tikhonov-TV方法与其他反卷积方法对比:其他反卷积方法,如迭代法等,仅在原始分布为平滑函数时较为适用,当原始成分分布为非平滑函数,如方波时,该方法便不适用了;其他反卷积方法,如傅里叶方法,对方波也适用,但当含有噪声时,便会对其结果产生巨大的影响。而采用正则化方法中的Tikhonov-TV方法却能够较好的解决低信噪比与数据不足的情况。
进一步的,步骤(1)所述深度剖析技术包括AES、XPS、SIMS、GDOES。
进一步的,步骤(2)中所述原子混合长度的确认方法为:利用SRIM软件及相应的溅射条件,对样品的膜层进行溅射模拟,取纵向岐离值作为MRI模型原子混合长度的值。
进一步的,步骤(2)中所述样品粗糙度的确定可以为利用AFM仪器进行测量。
进一步的,步骤(2)中所述信息深度的确认方法为:对于XPS、AES深度谱,利用NIST数据库查到相应光电子/俄歇电子的逃逸深度值,再与其出射角的余弦相乘;对于GDOES、SIMS深度谱,由于获取的信息源于样品的最表层,探测技术信息深度的影响可忽略不计,即λ=0。
与现有技术相比,本发明利用MRI模型和反卷积方法,来获得原始成分分布。针对反卷积过程中,低信噪比和实验数据不足等问题,采用正则化方法中的Tikhonov-TV方法,并结合其牛顿阻尼法,可以对噪声有着较好的容纳能力。对于数据点较少或不足的情况也具有较优的模拟效果,可以准确地获得原始成分深度浓度分布。
附图说明
图1是本发明从测量深度剖析数据直接获取薄膜元素成分分布的方法的流程图;
图2是本发明利用MRI模型与反卷积理论获取原始成分浓度分布的模拟结果,左图是方波,右图是平滑函数:
图3是本发明对含噪声的数据的模拟结果图;
图4是本发明对不连续数据的模拟结果图;
图5是实施例2测得的Ni/Cu多层膜的AES溅射谱;
图6是实施例2转换成的Ni/Cu多层膜的深度浓度谱;
图7是实施例2反卷后获得的Ni元素的原始成分分布;
图8是实施例2反卷后获得的Cu元素的原始成分分布。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步地详细描述。
实施例1
一种从测量深度剖析数据直接获取薄膜元素成分分布的方法如图1所示,主要包括以下步骤:
先利用AES、XPS、SIMS、GDOES深度剖析技术对薄膜样品进行深度剖析测量,并将测量得到的强度-溅射时间谱转换成浓度-深度谱。
再利用MRI(Mixing length-Roughness-Information depth)模型与反卷积方法来获得原始成分分布。MRI模型包含三个参数分别为原子混合长度(Mixing length),粗糙度(Roughness)以及测试仪器的信息深度(Information depth),这三个MRI参数有明确的物理意义,其中原子混合长度表示由于级联效应使得表层以下的原子与表层原子之间的混合,导致深度剖析谱的展宽;粗糙度包括样品的初始粗糙度和溅射诱导产生的粗糙度;信息深度表示由于仪器所获取的信号并非来自最外层的粒子单层收集到的,而是距表层一定深度范围内的信息。这三个参数既可以通过理论计算得到,也可以通过实验测量得到。其表达式分别如公式(1):
其中w为溅射导致的原子混合长度,σ为样品表面的粗糙度,λ为来自测量信号的信息深度。
MRI模型的解析表达式如如公式(2):
该解析式仅适用于MRI参数(w,σ,λ)为常数以及无择优溅射存在的情况。
在采用反卷积方法获取原始成分分布时,由于低信噪比与数据不足等,存在着不适定问题。为此,采用正则化方法解决该问题。正则化的核心是设法构造一个连续算子(正则算子)去逼近不连续算子,将不适定问题转化为适定问题,从而得到原问题的近似解。运用正则化方法,将反卷积转化为求解公式(3)(TV-Tikhonov function)的最小值。
其中α是正则化参数,A是一个与深度分辨率函数相关的矩阵,X表示成分深度浓度分布,b代表归一化信号强度。
TV(X)是一个全偏分(total variation)方程,表达如公式(4):
但是,公式(3)是非线性的,不利于数值方法的计算,比如最速下降法或牛顿法,为了解决这个缺陷,可用来近似表示欧几里得范数|X|,其中β是一个非常小的修正参数,所以,公式(4)可由公式(5)表示
根据公式(4)中微分的概念,可得
DiX=(Xi+1-Xi)/Δz i=1,....,n-1(6)
其中D是n*n-1维矩阵:
假设一个离散化的函数表达式:
其中由公式(3)可得:
利用牛顿法对公式(9)求解,首先,计算J(X)的梯度函数
grad J(X)=L(X)·X (10)
其中L(X)=Δz·DT·diag(ψ'(X))·D (11)
其中J(X)的海瑟矩阵为:
Hess J(X)=L(X)+L'(X)·X (12)
其中L(X)可由公式(11)获得,且
L'(X)·X=Δz·DT·diag(2(DX)2·ψ”(X))·D (13)
由公式(8)、公式(10)、公式(11),可得到公式(9)的梯度表达式:
gradT(X)=AT·(AX-b)+α·L(X)·X (14)
根据公式公式(12)、公式(13)、公式(14),可得到公式(9)的海瑟矩阵
Hess T(X)=AT·A+α·L(X)+α·L'(X)·X (15)
算法:Tikhonov-TV方法的牛顿阻尼解法
令j=0;
选择初始迭代向量X(0)(一般取零向量);
计算梯度:gi=AT·(AX-b)+α·L(X)·X;
计算海瑟矩阵:H=AT·A+α·L(X)+α·L'(X)·X;
计算下降方向:Sj=-H-1·gj;
算出最优步长:τj=argminτ>0T(X(j)+τSj);
计算第j+1次迭代解:X(j+1)=X(j)+τjSj;
更新当前迭代次数:j=j+1;
如果当前解符合要求,则停止;否则转到第(3)步继续迭代。
采用上述方法模拟结果如下图2所示,从图2中可以看出,该方法能有效获得原始成分分布。有很多算法在针对平滑函数时能获得较好的反卷结果,但对于方波确不行。但是本发明的方法在针对方波和平滑函数时都能获得较好的反卷结果。图3噪声以a*sin(z)的形式加入到信号中,由图3可知,本方法对噪声有着较好的容纳能力。图4中标注“数据间隔(5)”表示每隔五个数据点获取一个数据点作为新的原始数据,用以模拟数据点较少或不足的情况。模拟结果显示,该方法针对数据点不足有着良好的结果。
实施例2
本实施例利用MRI模型结合Tihonov-TV方法的牛顿阻尼解法,对Ni/Cu量多层膜的AES深度谱进行为例,说明具体从测量深度剖析数据直接获取薄膜元素成分分布的方法,具体包括以下步骤:
(一)利用俄歇电子能谱(AES)对Ni/Cu量多层膜进行深度剖析实验,归一化的测量数据(强度-溅射时间谱)如图5所示。
(二)将图5中的时间强度谱转换为浓度-深度谱(时间轴乘以溅射速率,溅射速率可由台阶仪获得的总的厚度推算出来),如图6所示。
(三)利用SRIM程序模拟MRI模型中的原子混合长度w,用AFM测量出样品表面粗糙度σ;由NIST数据库查得电子的逃逸深度,与电子发射角余弦值相乘作为其信息深度λ。(可得,w=1.6nm,λ=1.1nm,σ=1.6nm)。
(四)将相应的参数(如原子混合度等)代入Tikhonov-TV中的牛顿阻尼解法,直接获取了Ni与Cu原始层结构的成分分布。如图7和图8所示。
以上所揭露的仅为本发明的较佳实施例而已,当然不能以此来限定本发明之权利范围,因此依本发明权利要求所作的等同变化,仍属本发明所涵盖的范围。
Claims (5)
1.一种从测量深度剖析数据直接获取薄膜元素成分分布的方法,其特征在于,主要包括以下步骤:
(1)利用深度剖析技术对薄膜样品中元素的分布进行测量,并将测量得到的强度-溅射时间谱转换成浓度-深度谱;
(2)MRI参数的确定:包括原子混合长度、信息深度、样品粗糙度三个参数的确定;
(3)采用正则化方法中的Tikhonov-TV方法,并结合其牛顿阻尼法进行求解,获取原始成分深度浓度分布。
2.根据权利要求1所述从测量深度剖析数据直接获取薄膜元素成分分布的方法,其特征在于,步骤(1)所述深度剖析技术包括AES、XPS、SIMS、GDOES。
3.根据权利要求1所述从测量深度剖析数据直接获取薄膜元素成分分布的方法,其特征在于,步骤(2)中所述原子混合长度的确认方法为:利用SRIM软件及相应的溅射条件,对样品的膜层进行溅射模拟,取纵向岐离值作为MRI模型原子混合长度的值。
4.根据权利要求2所述从测量深度剖析数据直接获取薄膜元素成分分布的方法,其特征在于,步骤(2)中所述样品粗糙度的确定可以为利用AFM仪器进行测量。
5.根据权利要求2所述从测量深度剖析数据直接获取薄膜元素成分分布的方法,其特征在于,步骤(2)中所述信息深度的确认方法为:对于XPS、AES深度谱,利用NIST数据库查到相应光电子/俄歇电子的逃逸深度值,再与其出射角的余弦相乘;对于GDOES、SIMS深度谱,信息深度为0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810317643.7A CN108593695B (zh) | 2018-04-10 | 2018-04-10 | 从测量深度剖析数据直接获取薄膜元素成分分布的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810317643.7A CN108593695B (zh) | 2018-04-10 | 2018-04-10 | 从测量深度剖析数据直接获取薄膜元素成分分布的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108593695A true CN108593695A (zh) | 2018-09-28 |
CN108593695B CN108593695B (zh) | 2020-10-09 |
Family
ID=63621808
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810317643.7A Active CN108593695B (zh) | 2018-04-10 | 2018-04-10 | 从测量深度剖析数据直接获取薄膜元素成分分布的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108593695B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113281266A (zh) * | 2020-02-19 | 2021-08-20 | 神华(北京)光伏科技研发有限公司 | 膜层物质的光谱联用分析方法及其应用 |
CN115060755A (zh) * | 2022-08-18 | 2022-09-16 | 季华实验室 | 未知样品层结构深度剖析方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1060532A (zh) * | 1991-11-06 | 1992-04-22 | 复旦大学 | 无标样法sims定量分析 |
US20100224822A1 (en) * | 2009-03-05 | 2010-09-09 | Quebec Metal Powders, Ltd. | Insulated iron-base powder for soft magnetic applications |
JP2013036872A (ja) * | 2011-08-09 | 2013-02-21 | Ricoh Co Ltd | ラマン分光装置およびラマン分光測定法 |
CN106770620A (zh) * | 2016-11-21 | 2017-05-31 | 汕头大学 | 用溅射深度剖析技术确定薄膜中元素成分深度分布的方法 |
-
2018
- 2018-04-10 CN CN201810317643.7A patent/CN108593695B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1060532A (zh) * | 1991-11-06 | 1992-04-22 | 复旦大学 | 无标样法sims定量分析 |
US20100224822A1 (en) * | 2009-03-05 | 2010-09-09 | Quebec Metal Powders, Ltd. | Insulated iron-base powder for soft magnetic applications |
JP2013036872A (ja) * | 2011-08-09 | 2013-02-21 | Ricoh Co Ltd | ラマン分光装置およびラマン分光測定法 |
CN106770620A (zh) * | 2016-11-21 | 2017-05-31 | 汕头大学 | 用溅射深度剖析技术确定薄膜中元素成分深度分布的方法 |
Non-Patent Citations (2)
Title |
---|
F. BOULSINA ET AL.: "Deconvolution of SIMS depth profiles: Towards simple and faster techniques", 《APPLIED SURFACE SCIENCE》 * |
康红利 等: "溅射深度剖析定量分析及其应用研究进展", 《汕头大学学报(自然科学版)》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113281266A (zh) * | 2020-02-19 | 2021-08-20 | 神华(北京)光伏科技研发有限公司 | 膜层物质的光谱联用分析方法及其应用 |
CN113281266B (zh) * | 2020-02-19 | 2023-06-30 | 神华(北京)光伏科技研发有限公司 | 膜层物质的光谱联用分析方法及其应用 |
CN115060755A (zh) * | 2022-08-18 | 2022-09-16 | 季华实验室 | 未知样品层结构深度剖析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108593695B (zh) | 2020-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sobola et al. | Influence of scanning rate on quality of AFM image: Study of surface statistical metrics | |
Stoffelen | Toward the true near‐surface wind speed: Error modeling and calibration using triple collocation | |
CN103854305B (zh) | 一种基于多尺度建模的模型传递方法 | |
CN109358095B (zh) | 一种混合气体各组分气体浓度的定量测量方法和系统 | |
CN111879725B (zh) | 一种基于权重系数的光谱数据修正方法 | |
Pires et al. | Porosity distribution by computed tomography and its importance to characterize soil clod samples | |
Rab et al. | Comparison of Monte Carlo simulation, least square fitting and calibration factor methods for the evaluation of measurement uncertainty using direct pressure indicating devices | |
CN108593695A (zh) | 从测量深度剖析数据直接获取薄膜元素成分分布的方法 | |
Zhu et al. | Estimating soil water content from surface digital image gray level measurements under visible spectrum | |
Takeda et al. | Detection of gravitational redshift on the solar disk by using iodine-cell technique | |
Yang et al. | A relaxed matrix inversion method for retrieving water constituent concentrations in case II waters: The case of Lake Kasumigaura, Japan | |
CN106126879B (zh) | 一种基于稀疏表示技术的土壤近红外光谱分析预测方法 | |
KR100769566B1 (ko) | 신경망을 이용한 박막 두께 측정 방법, 장치 및 이를 위한기록매체 | |
Shard et al. | A simple approach to measuring thick organic films using the XPS inelastic background | |
Zhang et al. | An improved evaporation quantification method: Accounting for soil deformability and side evaporation effect | |
CN112818762A (zh) | 一种大尺寸复合材料及其夹层结构快速无损检测方法 | |
CN116932996A (zh) | 基于改进Hapke模型的土壤有机质含量的估测方法 | |
CN110132188A (zh) | 一种基于多元素x射线特征光谱综合分析的涂渗层厚度计算方法 | |
Brousaides | The radiosonde hygristor and low relative humidity measurements | |
CN109444055A (zh) | 一种盐碱地土壤盐分含量的高光谱计算方法 | |
CN108896511A (zh) | 一种智能化修复光谱分析仪谱图形变的方法 | |
Ren et al. | Hyper-Spectral Estimation Model of Soil Organic Matter Based on Generalized Greyness of Grey Number. | |
Vašát et al. | Predicting oxidizable carbon content via visible-and near-infrared diffuse reflectance spectroscopy in soils heavily affected by water erosion. | |
Chen et al. | Discussion on remote sensing estimation of soil nutrient contents | |
Morawiec | A method of texture inhomogeneity estimation from reflection pole figures |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230403 Address after: No. 419 Zhongjie Road, Qiandeng Town, Kunshan City, Suzhou City, Jiangsu Province, 215300 Patentee after: KUNSHAN SHUHOW INSTRUMENT TECHNOLOGY Co.,Ltd. Address before: 515000 No. 243, University Road, Shantou, Guangdong Patentee before: SHANTOU University |
|
TR01 | Transfer of patent right |