JP2023127779A - 混合物の回折パターン分析装置、方法、プログラム及び情報記憶媒体 - Google Patents
混合物の回折パターン分析装置、方法、プログラム及び情報記憶媒体 Download PDFInfo
- Publication number
- JP2023127779A JP2023127779A JP2022031672A JP2022031672A JP2023127779A JP 2023127779 A JP2023127779 A JP 2023127779A JP 2022031672 A JP2022031672 A JP 2022031672A JP 2022031672 A JP2022031672 A JP 2022031672A JP 2023127779 A JP2023127779 A JP 2023127779A
- Authority
- JP
- Japan
- Prior art keywords
- pattern
- fitting
- intensity ratio
- unknown
- diffraction
- 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
- 239000000203 mixture Substances 0.000 title claims abstract description 57
- 238000004458 analytical method Methods 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title description 28
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 238000002441 X-ray diffraction Methods 0.000 claims description 21
- 230000008859 change Effects 0.000 claims description 11
- 238000004445 quantitative analysis Methods 0.000 claims description 6
- 239000000126 substance Substances 0.000 description 22
- 238000011156 evaluation Methods 0.000 description 17
- 230000006870 function Effects 0.000 description 14
- 238000010586 diagram Methods 0.000 description 12
- 238000012937 correction Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 9
- 239000000523 sample Substances 0.000 description 7
- 229910018072 Al 2 O 3 Inorganic materials 0.000 description 4
- 238000005457 optimization Methods 0.000 description 3
- 230000000704 physical effect Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000012921 fluorescence analysis Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000012916 structural analysis Methods 0.000 description 2
- 238000003991 Rietveld refinement Methods 0.000 description 1
- OMOVVBIIQSXZSZ-UHFFFAOYSA-N [6-(4-acetyloxy-5,9a-dimethyl-2,7-dioxo-4,5a,6,9-tetrahydro-3h-pyrano[3,4-b]oxepin-5-yl)-5-formyloxy-3-(furan-3-yl)-3a-methyl-7-methylidene-1a,2,3,4,5,6-hexahydroindeno[1,7a-b]oxiren-4-yl] 2-hydroxy-3-methylpentanoate Chemical compound CC12C(OC(=O)C(O)C(C)CC)C(OC=O)C(C3(C)C(CC(=O)OC4(C)COC(=O)CC43)OC(C)=O)C(=C)C32OC3CC1C=1C=COC=1 OMOVVBIIQSXZSZ-UHFFFAOYSA-N 0.000 description 1
- PNEYBMLMFCGWSK-UHFFFAOYSA-N aluminium oxide Inorganic materials [O-2].[O-2].[O-2].[Al+3].[Al+3] PNEYBMLMFCGWSK-UHFFFAOYSA-N 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 229910052593 corundum Inorganic materials 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000000634 powder X-ray diffraction Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000013076 target substance Substances 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 229910001845 yogo sapphire Inorganic materials 0.000 description 1
Images
Classifications
-
- 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
- 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/2055—Analysing diffraction patterns
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/05—Investigating materials by wave or particle radiation by diffraction, scatter or reflection
- G01N2223/056—Investigating materials by wave or particle radiation by diffraction, scatter or reflection diffraction
- G01N2223/0566—Investigating materials by wave or particle radiation by diffraction, scatter or reflection diffraction analysing diffraction pattern
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/10—Different kinds of radiation or particles
- G01N2223/101—Different kinds of radiation or particles electromagnetic radiation
- G01N2223/1016—X-ray
Landscapes
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Crystallography & Structural Chemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
【課題】一部の成分についてだけ回折パターンが既知である場合にも、それらの成分の回折パターンの強度比を算出することができる混合物の回折パターン分析装置を提供する。【解決手段】混合物の回折パターン分析装置は、標的成分を示すとともに形状パラメータにより形状変更が可能な既知の標的パターン係る項と、残差群を示す未知パターンに係る項と、を含むフィッティングパターンを用いる。前記形状パラメータに所与の値を付与し、且つ前記未知パターンを初期パターンに設定した状態で、前記フィッティングパターンを観測パターンにフィッティングさせる。その後、前記未知パターンを変更することにより、前記フィッティングパターンを前記観測パターンにフィッティングさせる。複数の前記形状パラメータを用いて上記フィッティングを実行し、いずれかの形状パラメータに係る計算結果を選択する。【選択図】図2
Description
本発明は混合物の回折パターン分析装置、方法、プログラム及び情報記憶媒体に関し、X線回折の観測パターンに含まれる1又は複数の既知回折パターンの強度比を算出する技術に関する。
X線回折法を用いて混合物の定量分析を行うことができる。実際に観測される混合物の回折パターンには、各成分に由来する既知回折パターンが重なって含まれている。定量分析を実施する場合には、実際に観測される回折パターンにおける、各成分に由来する既知回折パターンの強度比を計算する。このような強度比が分かれば、例えば下記特許文献1~3に記載された方法により、各成分の重量分率を算出することができる。
上記従来の方法によれば、観測される回折パターンを複数の既知回折パターンに分解する際、すべての可能性のある成分について既知回折パターンを用意する必要がある。しかしながら、すべての成分について既知回折パターンを用意するのは、実際には困難な場合が多い。また、既知回折パターンの強度比が正確に分かれば、その反射として未知成分の回折パターンが分かることになる。これにより、その成分の構造解析等も実施することができるようになる。
本発明は上記課題に鑑みてなされたものであって、その目的は、一部の成分についてだけ回折パターンが既知である場合にも、それらの成分の回折パターンの強度比を精度よく算出することができる混合物の回折パターン分析方法、装置、プログラム及びコンピュータ可読情報記憶媒体を提供することにある。
上記課題を解決するために、混合物の回折パターン分析装置は、X線回折の観測パターンを取得する観測パターン取得手段と、標的成分を示すとともに形状パラメータにより形状変更が可能な既知の標的パターンに第1強度比を乗じた項と、1以上の残差成分からなる残差群を示す未知パターンに第2強度比を乗じた項と、を含み、前記第1強度比、前記第2強度比及び前記未知パターンがフィッティングパラメータであるフィッティングパターンを用い、前記形状パラメータに所与の値を付与し、且つ前記未知パターンを初期パターンに設定した状態で、前記第1強度比及び前記第2強度比を変更することにより、前記フィッティングパターンを前記観測パターンにフィッティングさせる第1フィッティング手段と、前記第1フィッティング手段によるフィッティングの後、前記第1強度比及び前記第2強度比の変更を制限しつつ、前記未知パターンを変更することにより、前記フィッティングパターンを前記観測パターンにフィッティングさせる第2フィッティング手段と、複数の前記形状パラメータを生成するパラメータ生成手段と、複数の前記形状パラメータのそれぞれを用いて、前記第1フィッティング手段及び前記第2フィッティング手段を動作させる制御手段と、前記未知パターンに基づいて、複数の前記形状パラメータのうちいずれかに係る前記第1強度比を選択する選択手段と、を含む。
ここで、前記形状パラメータは、前記標的成分の格子定数を含んでよい。また、前記形状パラメータは、前記既知の標的パターンに含まれるピークの幅に関する情報を含んでよい。
また、前記計算結果選択手段は、前記未知パターンを示す値に対して正値に変換する関数を適用してよい。また該関数の値の総和を算出してよい。そして、該総和に基づいて、複数の前記形状パラメータのいずれかに係る前記第1強度比を選択してよい。
また、前記計算結果選択手段は、前記未知パターンを示す値の変化量に対して正値に変換する関数を適用してよい。また該関数の値の総和を算出してよい。そして、該総和に基づいて、複数の前記形状パラメータのいずれかに係る前記第1強度比を選択してよい。
また、複数の前記形状パラメータのそれぞれを用いて、前記第1フィッティング手段と前記第2フィッティング手段を複数回ずつ動作させてよい。
また、前記フィッティングパターンは、さらにフィッティングパラメータである仮パターンの項を含んでよい。前記第1フィッティング手段は、前記第1強度比及び前記第2強度比に加えて、前記仮パターンを変更することにより、前記フィッティングパターンを前記観測パターンにフィッティングさせてよい。前記第2フィッティング手段は、前記仮パターンの項の少なくとも一部を前記残差群に係る項が吸収するように前記未知パターンを変更してよい。
また、前記第1強度比及び前記第2強度比に基づいて前記標的成分の定量分析を行ってよい。
また、本発明に係る混合物の回折パターン分析方法は、X線回折の観測パターンを取得する観測パターン取得ステップと、標的成分を示すとともに形状パラメータにより形状変更が可能な既知の標的パターンに第1強度比を乗じた項と、1以上の残差成分からなる残差群を示す未知パターンに第2強度比を乗じた項と、を含み、前記第1強度比、前記第2強度比及び前記未知パターンがフィッティングパラメータであるフィッティングパターンを用いる第1フィッティングステップであって、前記形状パラメータに所与の値を付与し、且つ前記未知パターンを初期パターンに設定した状態で、前記第1強度比及び前記第2強度比を変更することにより、前記フィッティングパターンを前記観測パターンにフィッティングさせる第1フィッティングステップと、前記第1フィッティングステップによるフィッティングの後、前記第1強度比及び前記第2強度比の変更を制限しつつ、前記未知パターンを変更することにより、前記フィッティングパターンを前記観測パターンにフィッティングさせる第2フィッティングステップと、前記形状パラメータを複数生成するパラメータ生成ステップと、を含み、複数の前記形状パラメータのそれぞれを用いて、前記第1フィッティングステップ及び前記第2フィッティングステップを実行させ、前記未知パターンに基づいて、複数の前記形状パラメータのいずれかに係る前記第1強度比を選択する。
また、本発明に係るプログラムは、X線回折の観測パターンを取得する観測パターン取得手段、標的成分を示すとともに形状パラメータにより形状変更が可能な既知の標的パターンに第1強度比を乗じた項と、1以上の残差成分からなる残差群を示す未知パターンに第2強度比を乗じた項と、を含み、前記第1強度比、前記第2強度比及び前記未知パターンがフィッティングパラメータであるフィッティングパターンを用いる第1フィッティング手段であって、前記形状パラメータに所与の値を付与し、且つ前記未知パターンを初期パターンに設定した状態で、前記第1強度比及び前記第2強度比を変更することにより、前記フィッティングパターンを前記観測パターンにフィッティングさせる第1フィッティング手段、前記第1フィッティング手段によるフィッティングの後、前記第1強度比及び前記第2強度比の変更を制限しつつ、前記未知パターンを変更することにより、前記フィッティングパターンを前記観測パターンにフィッティングさせる第2フィッティング手段と、前記形状パラメータを複数生成するパラメータ生成手段、複数の前記形状パラメータのそれぞれを用いて、前記第1フィッティング手段及び前記第2フィッティング手段を動作させる制御手段、及び前記未知パターンに基づいて、複数の前記形状パラメータのいずれかに係る前記第1強度比を選択する選択手段としてコンピュータを機能させるためのプログラムである。
また、本発明に係る情報記憶媒体は、上記プログラムを格納したコンピュータ可読情報記憶媒体である。
以下、本発明の一実施形態について図面に基づき詳細に説明する。
(システム構成)
図1は、本発明の一実施形態に係る分析システムの構成を示す図である。同図に示すように本分析システム10は、X線回折装置12、分析装置14、記憶部16及び表示部18を含んでいる。
図1は、本発明の一実施形態に係る分析システムの構成を示す図である。同図に示すように本分析システム10は、X線回折装置12、分析装置14、記憶部16及び表示部18を含んでいる。
X線回折装置12は、粉末X線回折測定を行う。具体的にはX線回折装置12は、試料物質に既知波長のX線を入射し、回折X線の強度を測定する。回折角度2θの値ごとのX線強度のデータは観測パターンとしてX線回折装置12から分析装置14に出力される。
ここで、本システムにより分析される試料物質は混合物であり、この混合物は1以上の標的成分(標的物質)と残差群を含んでいる。標的成分(Target component)とは、定量測定の対象となる成分である。残差群は1以上の残差成分(標的成分以外の物質)からなる。残差群の一例としては、どのような化学組成の成分がどのくらいの混合比で存在しているかが分かっているものの、その回折パターンについては未知であるものが挙げられる。また、他の例としては、混合物全体の化学組成は蛍光分析などにより分かっているものの、残差群については化学組成や混合比が分かっていないものが挙げられる。
分析装置14は、例えば公知のコンピュータシステムにより構成されており、演算装置やメモリを含んでいる。分析装置14にはSSD(Solid State Disk)やHDD(Hard Drive Disk)などのコンピュータ可読情報記憶媒体により構成された記憶部16が接続されている。記憶部16には、本発明の一実施形態に係る分析プログラムが記憶されており、この分析プログラムを分析装置14が実行することにより、本発明の一実施形態に係る分析方法が具現化される。
記憶部16にはさらに、各標的成分単体のX線回折パターンが標的パターンとして事前に記憶されている。これらの標的パターンは、コンピュータシステムによって所定のアルゴリズムを実行することにより生成される。標的パターンは、標的成分に関する種々の物性データから計算されるが、形状パラメータの変更によりその形状を調整できるようになっている。ここでは形状パラメータとして、回折線の位置を決める格子定数、及び回折線の幅を決める半値幅が用いられる。
記憶部16にはさらに、残差群のX線回折パターンである未知パターンの初期パターンが記憶されている。本実施形態では、この初期パターンを変更しながら、未知パターンを真のパターンに近づけていくものである。初期パターンは、残差群に含まれる成分のうち主なもののX線回折パターンであってもよい。或いは初期パターンは、残差群に含まれる複数成分の各X線回折パターンの線形結合であってもよい。しかし、初期パターンはこれに限らない。後述のように計算過程において初期パターンは適切なパターンに変更されるから、初期パターンとして、残差群に実際には含まれていない物質のX線回折パターンを用いることもできる。
記憶部16には、さらに各標的成分の化学組成情報(含有する原子の種類や原子量)や物性データを記憶している。なお、記憶部16は、試料である混合物全体の化学組成情報を記憶してもよい。
表示部18は、分析装置14による分析結果を表示する表示デバイスである。例えば表示部18は、標的パターンと未知パターンの強度比や、各標的成分の重量分率や残差群全体の重量分率を表示する。
(理論的背景)
ここで、分析装置14により行われるX線回折パターン分析の理論的背景について説明する。
ここで、分析装置14により行われるX線回折パターン分析の理論的背景について説明する。
試料である混合物がK個の成分を持つとき、k番目の重量分率wkは次式(1)で表される。
SkはLp補正を施したk番目の成分の強度の総和であり観測強度に相当する。またakは、記憶部16に記憶されているk番目の成分の化学組成情報から計算されるパラメータで、単位重量当たりの散乱強度の逆数に相当する。akは次式(2)で表される。
ここで、Mkはk番目の成分の化学式量である。nki'はk番目の成分を構成するi番目の原子に含まれている電子の数である。Σはk番目の成分を構成する全ての原子についての和を意味する。
次に、標的成分の個数をKT個(k=1~KT)、残差成分をK-KT個(k=KT+1~K)とすると、標的成分であるk番目の成分の重量分率wk(k=1~KT)は次式(3)で表される。
ここで、SRはLp補正を施した残差群の観測強度の総和である。すなわち、SRは次式(4)を意味する。
但し本実施形態では、残差成分の観測強度の総和SKT+1、SKT+2、…、SKをそれぞれ算出するのではなく、残差群全体に関するSRを算出する点が特徴の一つである。
また、aR_avは、残差群全体の化学組成情報から計算されるものである。aR_avの計算方法については後述する。
残差群全体の重量分率wRは次式(5)で表される。
標的成分についてak(k=1~KT)は既知であり、観測パターンに含まれる標的パターンの強度比からSk(k=1~KT)も算出可能である。後述するように、aR_avも複数の計算方法が存在する。さらに、以下に説明するように残差群に係る未知パターンも計算可能であり、その強度比からSRも計算可能である。このため、式(3)及び(5)により、各標的成分についての重量分率wk、及び残差群についての重量分率wRを求めることができる。
(フィッティングパターン)
フィッティングパターンは次式(6)で表される。
フィッティングパターンは次式(6)で表される。
ここでiは各回折角度を示す(i=1~N)。Yi
calcは、フィッティングパターンにおけるi番目の回折角度での強度を示す。Sck
Tは、k番目の標的成分の強度比を示す。Yki
Tは、k番目の標的成分の回折パターンにおけるi番目の回折角度での強度を示す。上述のようにYki
Tは形状パラメータによりパターン形状の調整が可能である。ここでは仮に格子定数及び半値幅の2つの形状パラメータにそれぞれ所与の値が設定されているものとする。ScRは、残差群の強度比を示す。Sci×Yi
Rは、残差群の回折パターンである未知パターンを示す。このうち、Yi
Rは、未知パターンの初期パターンを示す。具体的には、初期パターンにおけるi番目の回折角度での強度を示す。Sciは初期パターンに乗算される、i番目の回折角度での強度に対する補正因子である。Sciはフィッティング開始時にはすべて1に設定される。未知パターンの積分強度を一定にするため、Sciは次式(7)の拘束条件を有する。
すなわち、式(6)により示されるフィッティングパターンYi
calcは、標的成分を示す標的パターンYki
Tにその強度比を乗じた項と、残差群を示す未知パターンにその強度比を乗じた項とを含む。そして、2つの強度比と未知パターンがフィッティングパラメータである。
また、式(6)において、Yi
TMPは、仮パターンにおけるi番目の回折角度での強度を示す。Yi
TMPは、例えば各項の係数をフィッティングパラメータとする多項式やフーリエ級数を採用することができる。なお、Yi
TMPは、未知パターンの収束を良好なものとする役割の暫定的な項であり、計算終了時には零、又は極めて零に近い値となる。
(第1フィッティングステップ)
フィッティングの際には、まず未知パターンを初期パターンであるYi Rにした状態で、第1強度比であるSck T、第2強度比であるScR、及びYi TMPを変更し、フィッティングパターンYi calcを観測パターンYi obsにフィッティングさせる。具体的には、未知パターンを初期パターンであるYi Rにするため、Sciをすべて1に設定する。例えば最小二乗法等を用いて、式(6)に示されるYi calcとX線回折装置12から得られた観測パターンYi obsとの差が最小となるよう、Sck T、ScR及びYi TMPを決定する。
フィッティングの際には、まず未知パターンを初期パターンであるYi Rにした状態で、第1強度比であるSck T、第2強度比であるScR、及びYi TMPを変更し、フィッティングパターンYi calcを観測パターンYi obsにフィッティングさせる。具体的には、未知パターンを初期パターンであるYi Rにするため、Sciをすべて1に設定する。例えば最小二乗法等を用いて、式(6)に示されるYi calcとX線回折装置12から得られた観測パターンYi obsとの差が最小となるよう、Sck T、ScR及びYi TMPを決定する。
(第2フィッティングステップ)
次に、Sck T、ScR及びYi TMPを第1フィッティングステップで決定した値に固定した状態で、未知パターンSci×Yi Rを変更し、フィッティングパターンYi calcを観測パターンYi obsにフィッティングさせる。ここでは、Sciを変更することにより、未知パターンSci×Yi Rを変更する。
次に、Sck T、ScR及びYi TMPを第1フィッティングステップで決定した値に固定した状態で、未知パターンSci×Yi Rを変更し、フィッティングパターンYi calcを観測パターンYi obsにフィッティングさせる。ここでは、Sciを変更することにより、未知パターンSci×Yi Rを変更する。
具体的には、次式(8)によりSciを計算する。
式(8)に示されるSciは、式(7)の要件を満足しない。そこで次式(9)によりSciを正規化する。
ここで、Sci
newは正規化後のSciを示し、Sci
oldは式(8)の左辺を示す。また、SA及びSBkは次式(10)及び(11)で表される。
その後、正規化後のSciを用いて、再び第1フィッティングステップを実行する。すなわち、Yi
obsとYi
calcの誤差が収束するまで、第1フィッティングステップ及び第2フィッティングステップを複数回繰り返して実行する。なお、式(8)は、式(6)の右辺第1項及び第2項の和が観測パターンYi
obsに等しいと仮定した場合のSciの値を示している。これにより、仮パターンYi
TMPの値は、残差群に係る第2項により吸収されるようになる。このため計算終了時には、仮パターンYi
TMPの値は零、又は極めて零に近い値に収束する。これによりフィッティングパターンYi
calcは観測パターンYi
obsにフィッティングされる。
以上の説明では標的パターンYki
Tの2つの形状パラメータ(格子定数及び半値幅)にそれぞれ所与の値が設定されているものとした。本実施形態では、それらの値を変更しながら、複数回にわたって上述のパラメータフィッティングを実行する。これにより、標的パターンYki
Tの2つの形状パラメータの値を最適化する。後述のように標的パターンYki
Tに誤差が含まれていると、当該誤差に対応する成分(対応誤差成分)が残差群に係る未知パターンにも含まれることになる。これにより、未知パターンが不正確なものとなり、最終的な計算結果であるSck
T、ScRにも誤差が含まれることになる。そこで、標的パターンYki
Tの2つの形状パラメータの値を最適化することで、このような対応誤差成分を最小化することができる。どのような形状パラメータの値が最適であるかは、後述するように未知パターンSci×Yi
Rに基づいて判断することができる。
図2は、分析装置14の動作を示すフロー図である。分析装置14により実行される分析プログラムは同図において左側に示される形状パラメータ直接探索プログラム(S201~S204)と、右側に示されるパターンフィッティングプログラム(S101~S103)と、両プログラムの橋渡しを担うプログラム(S301~S303)を含んでいる。
分析装置14では、まず標的パターンYki
T、及びその可変パラメータである形状パラメータXpを指定する(S101)。ここでは形状パラメータXpとして標的成分の格子定数及びピークの半値幅の2つが用いられる。形状パラメータXpのセット{Xp}は形状パラメータ直接探索プログラムに通知される(S301)。
形状パラメータ直接探索プログラムは、最適化アルゴリズムを実行するものであり、ここでは最適化アルゴリズムの一例として直接探索法であるネルダー-ミード法(シンプレックス法/アメーバ法)が用いられる。形状パラメータ直接探索プログラムは、まず形状パラメータXpのセット{Xp}に対応する2次元シンプレックスの初期状態を設定する。具体的には2次元シンプレックスの3つの頂点のそれぞれに初期座標を設定する(S201)。各初期座標は、それぞれ予め記憶部16に記憶された座標が設定されてもよいし、図示しない入力デバイスにより手入力で設定されてもよい。各頂点に対応する形状パラメータのセット{Xp}はパターンフィッティングプログラムに転送される(S302)。
パターンフィッティングプログラムは、形状パラメータの各セット{Xp}に基づいてシンプレックスの各頂点に対応する標的パターンYki
Tを生成する(S102)。この際、記憶部16に事前に記憶されている標的成分に関する各種物性データも参照される。そして、各頂点に対応する標的パターンYki
Tを順に用いて後述のパターンフィッティングを実行する(S103)。
すべての頂点についてパターンフィッティングを終えると、各頂点に対応するフィッティングパターンYi
calcのScR、Sci、Yi
Rを形状パラメータ直接探索プログラムに転送する(S303)。
形状パラメータ直接探索プログラムは、各頂点に対応するフィッティングパターンYi
calcのScR、Sci、Yi
Rに基づいて、各頂点に対応する評価値Snを計算する(S203)。評価値Snは未知パターンSci×Yi
Rに基づくものである。一例として評価値Snの次式(12)に示されるものであってよい。ここでnは任意の自然数である(1以上の正数であってもよい)。
式(12)に示される評価値Snは、未知パラメータSci×Yi
Rの各回折角度における値の変化量に対して正値に変換する関数(ここでは絶対値のべき乗)を適用し、該関数の値の総和を演算するものである。これにより、未知パラメータSci×Yi
Rの変化値の総量を評価することができる。未知パラメータSci×Yi
Rに上述の対応誤差成分が含まれる場合には、未知パラメータSci×Yi
Rの変化量の総量が大きくなる傾向にある。そこで、式(12)に示される評価値Snを最小化することで、標的パターンYki
Tに対して適切な形状パターンのセット{Xp}を選択するとともに、未知パラメータSci×Yi
Rに含まれる対応誤差成分を小さくすることができる。なお、ここでは上記関数が適用される値を未知パターンそのものの変化量としたが、未知パターンにScRを乗じた値の変化量に上記関数を適用してもよい。或いはYi
Rが固定値の場合にはSciは未知パターンと同視できるから、Sciの変化量に上記関数を適用してもよい。
評価値Snが十分に収束していない場合には、2次元シンプレックスの3頂点の座標を変更し(S202)、再び、S302以降の処理を実行する。座標の変更には、ネルダー-ミード法の4つの操作(鏡像、拡大、収縮、縮小)が用いられる。
評価値Snの収束は、例えばシンプレックスの3頂点に対応する3つの評価値Snが所定条件(相互に十分に近い値であるかどうか)を満足するかにより判断してよい。或いは、直前のシンプレックスに対応する評価値Snの最小値と、今回のシンプレックスに対応する評価値Snの最小値が所定条件(十分に近い値であるかどうか)を満足するかにより判断してもよい。
形状パラメータ直接探索プログラムは、S204において評価値Snが十分に収束していると判断すれば処理を終了する。このとき、評価値Snの最小値に対応するSck
T、ScRが計算結果として出力される。すなわち、未知パターンSci×Yi
Rに基づいて評価値Snが計算され、該評価値Snに基づいてSck
T、ScRが選択される。
なお、評価値Snの変形例として次式(13)を採用してもよい。ここでもnは任意の自然数である(1以上の正数であってもよい)。
式(13)に示される評価値Snは、未知パターンSci×Yi
Rの各回折角度における値に対して正値に変換する関数(ここでは絶対値のべき乗)を適用し、その関数の値の総和を演算するものである。これにより、未知パターンSci×Yi
Rの振幅の総量を評価することができる。未知パターンSci×Yi
Rに上述の対応誤差成分が含まれる場合には、未知パターンSci×Yi
Rの振幅の総量が大きくなる傾向にある。そこで、式(13)に示される評価値Snを最小化することで、標的パターンYki
Tに対して適切な形状パターンのセット{Xp}を選択するとともに、未知パターンSci×Yi
Rに含まれる対応誤差成分を小さくすることができる。なお、ここでは上記関数が適用される値を未知パターンそのものとしたが、未知パターンにScRを乗じた値に上記関数を適用してもよい。或いはYi
Rが固定値の場合にはSciは未知パターンと同視できるから、Sciに上記関数を適用してもよい。
図3は、分析装置14における上述のパターンフィッティング(S103)の処理を示すフロー図である。
分析装置14は、まずX線回折装置12から観測パターンYi
obsを取得する(S1031)。さらに、図2のS102において生成された複数の標的パターンYki
Tのうち1つを読み出す(S1032)。その後、式(10)及び(11)によりSA及びSBkを算出する(S1033)。さらに、初期パターンYi
Rを記憶部16から読み出す(S1034)。
その後、補正パターンSciの値をすべて1に初期化し(S1035)、上述の第1フィッティングステップを実行する(S1036)。すなわち、式(6)に示されるYi
calcと観測パターンYi
obsとの誤差が最小化されるようにして、フィッティングパラメータであるSck
T、ScR及びYi
TMPを決定する。S1036においては、式(6)に示されるYi
calcが取得され、S1032、S1034及びS1035で得られる値が代入される。
次に、式(8)により補正パターンSci(正規化前)を算出し(S1037)、式(9)によりそれを正規化する(S1038)。
以上のS1035~S1038の処理を、Yi
calcと観測パターンYi
obsとの誤差が収束条件を満足するまで繰り返す(S1039)。この処理をすべての標的パターンYki
Tについて実行し、フィッティングパターンYi
calcの第2項に含まれるScR、Sci、Yi
Rを出力する。
以上のようにして形状パラメータを最適化しつつ、フィッティングパターンYi
calcを観測パターンYi
obsに一致させると、その結果を用いて定量分析を行うことができる。また、未知パターンSci×Yi
Rを用い、リートベルト法等により未知成分の構造解析を実施することができる。
定量分析を実施する場合、Sck
T及びScRを用いてSk及びSRの値を算出する。例えば、Yki
TとYi
Rとが予め規格化されていれば、SkはSck
Tと等しく、SRはScRと等しい。そして、それらの値を式(3)に代入し、重量分率wkを算出する。また重量分率wkの値を式(5)に代入し、残差群に係る重量分率wRを算出する。
(aR_avの算出方法(1))
ここで、aR_avの算出方法について説明する。
ここで、aR_avの算出方法について説明する。
残差群に対してどのような化学組成の成分がどれくらいの混合比で存在しているか分かっている場合、それらの情報から直接aR_avを求めることができる。
すなわち、残差群が物質A(WAg)及び物質B(WBg)からなるとき、残差群は次式(14)に示される散乱強度を与える。
この散乱強度を残差群の総重量で割れば、単位重量当たりの散乱強度、即ちaR_avが求まる。すなわち、aR_avは次式(15)で与えられる。
式(15)をK-T個の成分から成る残差群に一般化し、重量分率wkを用いてaR_avを表すと、次式(16)のようになる。
ここで、ak'については化学組成情報から式(2)を用いて算出することができる。したがって、残差群にどのような化学組成の成分がどのくらいの混合比で存在しているかが分かっている場合には、同式(16)によりaR_avを算出することができる。
(aR_avの算出方法(2))
次に、混合物全体(バッチ組成)の化学組成情報が分かっているが、残差群に関して化学組成情報が分かっていない場合について説明する。バッチ組成の化学組成情報は、例えばバッチ組成に対して蛍光分析等を適用することで求めることができる。或いは、揮発性成分がないと仮定できる場合には、混合物の合成に用いた原料の化学組成情報をそのまま用いることができる。
次に、混合物全体(バッチ組成)の化学組成情報が分かっているが、残差群に関して化学組成情報が分かっていない場合について説明する。バッチ組成の化学組成情報は、例えばバッチ組成に対して蛍光分析等を適用することで求めることができる。或いは、揮発性成分がないと仮定できる場合には、混合物の合成に用いた原料の化学組成情報をそのまま用いることができる。
このような場合、バッチ組成の化学組成情報を式(2)に代入し、混合物試料全体に対するakを算出する。この値をaBと表記する。
式(16)と同様に、バッチ組成のaBは次式(17)で表される。
式(17)を変形し、残差群に係るaR_avは次式(18)により表される。
同式(18)において、aB及びak'は既知であるが、重量分率wR及びwk'は未知である。そこで、例えばaR_avの初期値をaBと仮定して、式(3)(5)より重量分率wR及びwk'を計算し、これを再び式(18)に代入し、aR_avを再度算出する。これを繰り返すことにより、真値に近いaR_avを計算することができる。
図4は、標的成分の重量分率wk、及び残差群の重量分率wRを計算する処理を示すフロー図である。同図に示すように、分析装置14は標的成分の化学組成情報を記憶部16から読み出し、式(2)によりakを計算する(S401)。さらに、分析装置14は残差群に係るaR_avを算出する(S402)。例えば、残差群にどのような化学組成の成分がどれくらいの混合比で存在しているか分かっている場合、式(16)によりaR_avを算出する。
その後、式(3)により標的成分の重量分率wkを算出する(S403)。さらに、式(5)により残差群の重量分率wRを算出する(S404)。そして、重量分率wk及びwRを表示部18に表示させる(S405)。
図5は、変数aR_avの計算手順を示すフロー図である。同図に示す処理は、図4に示されるS402の処理の一例である。分析装置14は、記憶部16からバッチ組成の化学組成情報を読み出し、バッチ組成に係るaB
-1を式(17)により計算する(S4021)。次に、aR_av
-1の初期値としてaB
-1を設定し(S4022)、式(3)及び式(5)により重量分率wk及びwRを計算する(S4023)。それらの値を式(18)に代入し、aR_av
-1を計算する(S4024)。aR_av
-1が収束条件を満足するまでS4023及びS4024の処理を繰り返し、収束条件を満足すればaR_av
-1を出力する(S4026)。この値は、S403の処理で用いられてよい。
図6及び図7は、分析システム10による混合物試料のX線回折パターンの分析例を示す図である。試料はα-Al2O3とγ-Al2O3の混合物(混合比75:25)である。標的成分はα- Al2O3である。残差群は単成分でありγ- Al2O3である。残差群の回折パターンは未知である。但し、図6では、標的パターンYki
Tの形状パラメータのセット{Xp}を最適化していない。一方、図7では、形状パラメータ直接探索プログラムにより形状パラメータのセット{Xp}を最適化している。
図6において、符号20には観測パターンYi
obsが示されている。符号22には補正パターンSciが示されている。補正パターンSciには標的成分のピークの位置に鋭い立ち上がりや立下りが表れている(一例として符号24参照)。これは標的パターンYki
Tの誤差に対応してするものであり、上述の対応誤差成分である。
一方、図7では、符号26に示される補正パターンSciにおいて、標的成分のピークの位置に生じる立上がりや立下りが小さくなっている(一例として符号28参照)。このように、形状パラメータ直接探索プログラムにより形状パラメータのセット{Xp}を最適化すると、補正パターンSciにおいて、標的成分のピークの位置に表れる対応誤差成分を小さくすることができる。この結果、標的成分の重量分率wk、及び残差群の重量分率wRを精度よく算出することが可能となる。なお、本実施形態の方法によりγ-Al2O3の重量分率を計算したところ、真値に対して誤差は1%程度であった。
以上説明した本実施形態によれば、一部の成分についてだけ回折パターンが既知である場合にも、回折パターンが既知の成分のみならず、回折パターンが未知の残差群についても、強度比を算出することができ、またそれら成分の重量分率を精度よく算出することができる。また、標的パターンYki
Tを形状パラメータのセット{Xp}により形状変更可能とし、形状パラメータのセット{Xp}の各値を最適化アルゴリズムにより最適化したので、標的成分及び残差群について重量分率をさらに精度よく算出することができる。また、未知パターンSci×Yi
R精度よく生成できるので、未知成分の構造解析や成分同定を精度よく行うことができる。
なお、本発明は上記実施形態に限定されるものではなく、種々の変形実施が可能である。例えば、各パターンの強度比は上記説明とは異なる方法により算出してもよい。
10 分析システム、12 X線回折装置、14 分析装置、16 記憶部、18 表示部、20 観測パターン、22 未知パターン、24,28 対応誤差成分。
Claims (11)
- X線回折の観測パターンを取得する観測パターン取得手段と、
標的成分を示すとともに形状パラメータにより形状変更が可能な既知の標的パターンに第1強度比を乗じた項と、1以上の残差成分からなる残差群を示す未知パターンに第2強度比を乗じた項と、を含み、前記第1強度比、前記第2強度比及び前記未知パターンがフィッティングパラメータであるフィッティングパターンを用い、前記形状パラメータに所与の値を付与し、且つ前記未知パターンを初期パターンに設定した状態で、前記第1強度比及び前記第2強度比を変更することにより、前記フィッティングパターンを前記観測パターンにフィッティングさせる第1フィッティング手段と、
前記第1フィッティング手段によるフィッティングの後、前記第1強度比及び前記第2強度比の変更を制限しつつ、前記未知パターンを変更することにより、前記フィッティングパターンを前記観測パターンにフィッティングさせる第2フィッティング手段と、
複数の前記形状パラメータを生成するパラメータ生成手段と、
複数の前記形状パラメータのそれぞれを用いて、前記第1フィッティング手段及び前記第2フィッティング手段を動作させる制御手段と、
前記未知パターンに基づいて、複数の前記形状パラメータのいずれかに係る前記第1強度比を選択する選択手段と、
を含む、混合物の回折パターン分析装置。 - 請求項1に記載の混合物の回折パターン分析装置において、
前記形状パラメータは、前記標的成分の格子定数を含む、混合物の回折パターン分析装置。 - 請求項1又は2に記載の混合物の回折パターン分析装置において、
前記形状パラメータは、前記既知の標的パターンに含まれるピークの幅に関する情報を含む、混合物の回折パターン分析装置。 - 請求項1乃至3のいずれかに記載の混合物の回折パターン分析装置において、
前記計算結果選択手段は、前記未知パターンを示す値に対して正値に変換する関数を適用するとともに、該関数の値の総和を算出し、該総和に基づいて、複数の前記形状パラメータのいずれかに係る前記第1強度比を選択する、混合物の回折パターン分析装置。 - 請求項1乃至3のいずれかに記載の混合物の回折パターン分析装置において、
前記計算結果選択手段は、前記未知パターンを示す値の変化量に対して正値に変換する関数を適用するとともに、該関数の値の総和を算出し、該総和に基づいて、複数の前記形状パラメータのいずれかに係る前記第1強度比を選択する、混合物の回折パターン分析装置。 - 請求項1乃至5のいずれかに記載の回折パターン分析装置において、
複数の前記形状パラメータのそれぞれを用いて、前記第1フィッティング手段と前記第2フィッティング手段を複数回ずつ動作させることを特徴とする混合物の回折パターン分析装置。 - 請求項1乃至6のいずれかに記載の回折パターン分析装置において、
前記フィッティングパターンは、さらにフィッティングパラメータである仮パターンの項を含み、
前記第1フィッティング手段は、前記第1強度比及び前記第2強度比に加えて、前記仮パターンを変更することにより、前記フィッティングパターンを前記観測パターンにフィッティングさせ、
前記第2フィッティング手段は、前記仮パターンの項の少なくとも一部を前記残差群に係る項が吸収するように前記未知パターンを変更する、
ことを特徴とする混合物の回折パターン分析装置。 - 請求項1乃至7のいずれかに記載の回折パターン分析装置において、
前記第1強度比及び前記第2強度比に基づいて前記標的成分の定量分析を行うことを特徴とする混合物の回折パターン分析装置。 - X線回折の観測パターンを取得する観測パターン取得ステップと、
標的成分を示すとともに形状パラメータにより形状変更が可能な既知の標的パターンに第1強度比を乗じた項と、1以上の残差成分からなる残差群を示す未知パターンに第2強度比を乗じた項と、を含み、前記第1強度比、前記第2強度比及び前記未知パターンがフィッティングパラメータであるフィッティングパターンを用いる第1フィッティングステップであって、前記形状パラメータに所与の値を付与し、且つ前記未知パターンを初期パターンに設定した状態で、前記第1強度比及び前記第2強度比を変更することにより、前記フィッティングパターンを前記観測パターンにフィッティングさせる第1フィッティングステップと、
前記第1フィッティングステップによるフィッティングの後、前記第1強度比及び前記第2強度比の変更を制限しつつ、前記未知パターンを変更することにより、前記フィッティングパターンを前記観測パターンにフィッティングさせる第2フィッティングステップと、
前記形状パラメータを複数生成するパラメータ生成ステップと、を含み、
複数の前記形状パラメータのそれぞれを用いて、前記第1フィッティングステップ及び前記第2フィッティングステップを実行させ、
前記未知パターンに基づいて、複数の前記形状パラメータのいずれかに係る前記第1強度比を選択する、
混合物の回折パターン分析方法。 - X線回折の観測パターンを取得する観測パターン取得手段、
標的成分を示すとともに形状パラメータにより形状変更が可能な既知の標的パターンに第1強度比を乗じた項と、1以上の残差成分からなる残差群を示す未知パターンに第2強度比を乗じた項と、を含み、前記第1強度比、前記第2強度比及び前記未知パターンがフィッティングパラメータであるフィッティングパターンを用いる第1フィッティング手段であって、前記形状パラメータに所与の値を付与し、且つ前記未知パターンを初期パターンに設定した状態で、前記第1強度比及び前記第2強度比を変更することにより、前記フィッティングパターンを前記観測パターンにフィッティングさせる第1フィッティング手段、
前記第1フィッティング手段によるフィッティングの後、前記第1強度比及び前記第2強度比の変更を制限しつつ、前記未知パターンを変更することにより、前記フィッティングパターンを前記観測パターンにフィッティングさせる第2フィッティング手段と、
前記形状パラメータを複数生成するパラメータ生成手段、
複数の前記形状パラメータのそれぞれを用いて、前記第1フィッティング手段及び前記第2フィッティング手段を動作させる制御手段、及び
前記未知パターンに基づいて、複数の前記形状パラメータのいずれかに係る前記第1強度比を選択する選択手段
としてコンピュータを機能させるためのプログラム。 - 請求項10に記載のプログラムを格納したコンピュータ可読情報記憶媒体。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2022031672A JP2023127779A (ja) | 2022-03-02 | 2022-03-02 | 混合物の回折パターン分析装置、方法、プログラム及び情報記憶媒体 |
DE102023104071.5A DE102023104071A1 (de) | 2022-03-02 | 2023-02-20 | Vorrichtung und verfahren zum analysieren eines beugungsmusters einer mischung, und informationsspeichermedium |
US18/174,673 US20230280290A1 (en) | 2022-03-02 | 2023-02-27 | Device and method for analyzing diffraction pattern of mixture, and information storage medium |
CN202310190251.XA CN116698892A (zh) | 2022-03-02 | 2023-03-02 | 混合物的衍射图案分析装置、方法以及信息存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2022031672A JP2023127779A (ja) | 2022-03-02 | 2022-03-02 | 混合物の回折パターン分析装置、方法、プログラム及び情報記憶媒体 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2023127779A true JP2023127779A (ja) | 2023-09-14 |
Family
ID=87572202
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2022031672A Pending JP2023127779A (ja) | 2022-03-02 | 2022-03-02 | 混合物の回折パターン分析装置、方法、プログラム及び情報記憶媒体 |
Country Status (4)
Country | Link |
---|---|
US (1) | US20230280290A1 (ja) |
JP (1) | JP2023127779A (ja) |
CN (1) | CN116698892A (ja) |
DE (1) | DE102023104071A1 (ja) |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6658089B2 (ja) | 2016-02-26 | 2020-03-04 | 日本ポリプロ株式会社 | ポリオレフィンの製造方法 |
JP7052249B2 (ja) | 2017-08-08 | 2022-04-12 | 株式会社リコー | 画像形成装置、液体を吐出する装置 |
JP6930737B2 (ja) | 2018-04-02 | 2021-09-01 | 株式会社リガク | 非晶質相の定量分析装置、非晶質相の定量分析方法、及び非晶質相の定量分析プログラム |
-
2022
- 2022-03-02 JP JP2022031672A patent/JP2023127779A/ja active Pending
-
2023
- 2023-02-20 DE DE102023104071.5A patent/DE102023104071A1/de active Pending
- 2023-02-27 US US18/174,673 patent/US20230280290A1/en active Pending
- 2023-03-02 CN CN202310190251.XA patent/CN116698892A/zh active Pending
Also Published As
Publication number | Publication date |
---|---|
US20230280290A1 (en) | 2023-09-07 |
CN116698892A (zh) | 2023-09-05 |
DE102023104071A1 (de) | 2023-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
De Nolf et al. | XRDUA: crystalline phase distribution maps by two-dimensional scanning and tomographic (micro) X-ray powder diffraction | |
JP3699723B2 (ja) | 材料サンプルの分析 | |
EP2024734B1 (en) | A method of determining the feasibility of a proposed structure analysis process | |
US7516031B2 (en) | Apparatus and method for correcting for aberrations | |
Tomii et al. | Renormalization of domain-wall bilinear operators with short-distance current correlators | |
Mendenhall et al. | An implementation of the fundamental parameters approach for analysis of X-ray powder diffraction line profiles | |
Schuetzke et al. | Enhancing deep-learning training for phase identification in powder X-ray diffractograms | |
JPH09297112A (ja) | 構造パラメータ解析装置及び解析方法 | |
Zendler et al. | VITESS 3–virtual instrumentation tool for the European Spallation Source | |
Zhang et al. | New capabilities for enhancement of RMCProfile: instrumental profiles with arbitrary peak shapes for structural refinements using the reverse Monte Carlo method | |
Bicudo et al. | Study of I= 0 bottomonium bound states and resonances in S, P, D, and F waves with lattice QCD static-static-light-light potentials | |
JP6930737B2 (ja) | 非晶質相の定量分析装置、非晶質相の定量分析方法、及び非晶質相の定量分析プログラム | |
JP2023127779A (ja) | 混合物の回折パターン分析装置、方法、プログラム及び情報記憶媒体 | |
Toby et al. | Determination of standard uncertainties in fits to pair distribution functions | |
JP7414280B2 (ja) | 混合物の回折パターン分析方法、装置、プログラム及び情報記憶媒体 | |
US10416101B2 (en) | Model independent grazing incidence X-ray reflectivity | |
JP4115542B2 (ja) | ディフラクトグラムによる無標準相分析法 | |
Müller et al. | Analysis of the random measurement error of areal 3D coordinate measurements exclusively based on measurement repetitions | |
Spinolo et al. | Cohen's method revisited | |
RU2240543C2 (ru) | Способ рентгенофлуоресцентного анализа элементного состава вещества | |
JP3954936B2 (ja) | ブラッグ反射条件シミュレーション装置およびブラッグ反射測定システム | |
Iro | Further development of the software jGiXa for grazing incidence X-ray fluoresence analysis | |
Charnock et al. | Planck confronts large scale structure: methods to quantify discordance | |
BARCELLOS LINS | Development of a 2D-XRF scanner for cultural heritage applications | |
Hubert | Structural refinement of single crystals using digital-large angle convergent beam electron diffraction: PhD Thesis |