QUANTITATIVE PHASE ANALYSIS DEVICE, QUANTITATIVE PHASE ANALYSIS METHOD, AND NON-TRANSITORY COMPUTER-READABLE STORAGE MEDIUM STORING QUANTITATIVE PHASE ANALYSIS PROGRAM
20200173938 ยท 2020-06-04
Inventors
Cpc classification
International classification
Abstract
A quantitative phase analysis device includes: a unit for acquiring a powder diffraction pattern of the sample; a unit for acquiring information on a plurality of crystalline phases; a unit for acquiring a fitting function for each of the plurality of crystalline phases; a unit for executing whole-powder pattern fitting for the powder diffraction pattern by using the acquired fitting functions, to thereby acquire a fitting result; and a unit for calculating a weight ratio of the plurality of crystalline phases based on the fitting result. Each fitting function is selected from the group consisting of a first fitting function using an integrated intensity obtained by whole-powder pattern decomposition, a second fitting function using an integrated intensity obtained by observation or calculation, and a third fitting function using a profile intensity obtained by observation or calculation.
Claims
1. A quantitative phase analysis device, which is configured to perform quantitative phase analysis of crystalline phases contained in a sample based on a powder diffraction pattern of the sample, the quantitative phase analysis device comprising: at least one processor; and at least one memory device that stores a plurality of instructions, which when executed by the at least one processor, cause the at least one processor to: acquire the powder diffraction pattern of the sample; acquire information on a plurality of crystalline phases contained in the sample; acquire a fitting function for each of the plurality of crystalline phases; execute whole-powder pattern fitting for the powder diffraction pattern of the sample by using the fitting function for each of the plurality of crystalline phases, to thereby acquire a fitting result; and calculate a weight ratio of the plurality of crystalline phases based on the fitting result, wherein the fitting function for each of the plurality of crystalline phases is one fitting function selected from the group consisting of a first fitting function that uses an integrated intensity obtained by whole-powder pattern decomposition, a second fitting function that uses an integrated intensity obtained by observation or calculation, and a third fitting function that uses a profile intensity obtained by observation or calculation.
2. The quantitative phase analysis device according to claim 1, wherein calculating the weight ratio includes calculating a weight fraction by using an intensity-composition (IC) formula.
3. The quantitative phase analysis device according to claim 1, wherein the quantitative phase analysis device is configured to select, for the plurality of crystalline phases, two or more types of fitting functions among the first fitting function to the third fitting function.
4. A quantitative phase analysis method for performing quantitative phase analysis of crystalline phases contained in a sample based on a powder diffraction pattern of the sample, the quantitative phase analysis method comprising: acquiring the powder diffraction pattern of the sample; acquiring information on a plurality of crystalline phases contained in the sample; acquiring a fitting function for each of the plurality of crystalline phases; executing whole-powder pattern fitting for the powder diffraction pattern of the sample by using the fitting function for each of the plurality of crystalline phases, to thereby acquire a fitting result; and calculating a weight ratio of the plurality of crystalline phases based on the fitting result, wherein the fitting function for each of the plurality of crystalline phases is one fitting function selected from the group consisting of a first fitting function that uses an integrated intensity obtained by whole-powder pattern decomposition, a second fitting function that uses an integrated intensity obtained by observation or calculation, and a third fitting function that uses a profile intensity obtained by observation or calculation.
5. A non-transitory computer-readable storage medium storing a quantitative phase analysis program for performing quantitative phase analysis of crystalline phases contained in a sample based on a powder diffraction pattern of the sample, the quantitative phase analysis program causing a computer to execute processing comprising: acquiring the powder diffraction pattern of the sample; acquiring information on a plurality of crystalline phases contained in the sample; acquiring a fitting function for each of the plurality of crystalline phases; executing whole-powder pattern fitting for the powder diffraction pattern of the sample by using the fitting function for each of the plurality of crystalline phases, to thereby acquire a fitting result; and calculating a weight ratio of the plurality of crystalline phases based on the fitting result, wherein the fitting function for each of the plurality of crystalline phases is one fitting function selected from the group consisting of a first fitting function that uses an integrated intensity obtained by whole-powder pattern decomposition, a second fitting function that uses an integrated intensity obtained by observation or calculation, and a third fitting function that uses a profile intensity obtained by observation or calculation.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
[0039]
[0040]
[0041]
[0042]
[0043]
[0044]
DETAILED DESCRIPTION
[0045] Now, an embodiment of the present invention is described with reference to the drawings. For clearer illustration, some sizes, shapes, and the like are schematically illustrated in the drawings in comparison to actual ones. However, the sizes, the shapes, and the like are merely examples, and do not limit understanding of the present invention. Further, like elements as those described relating to the drawings already referred to are denoted by like reference symbols herein and in each of the drawings, and detailed description thereof is sometimes omitted as appropriate.
[0046]
[0047] The quantitative phase analysis device 1 according to this embodiment includes an analysis unit 2, an information input unit 3, an information output unit 4, and a storage unit 5. The quantitative phase analysis device 1 is achieved by a computer used in general, and further includes a read only memory (ROM) (not shown) and a random access memory (RAM) (not shown). The ROM and the RAM form internal memories of the computer. The storage unit 5 is a recording medium, and may be formed of a semiconductor memory, a hard disk drive, or other such arbitrary recording medium. In this case, the storage unit 5 is installed inside the computer, but may be installed outside the computer. The storage unit 5 may be a single recording medium, or may be formed of a plurality of recording mediums. The quantitative phase analysis device 1 is connected to an X-ray diffractometer 11 and an input device 13. The X-ray diffractometer 11 is configured to subject a sample having a powder shape to X-ray diffraction measurement to measure X-ray diffraction data on the sample and output the measured X-ray diffraction data to the information input unit 3 of the quantitative phase analysis device 1. The input device 13 is realized with a keyboard, a mouse, a touch panel, or the like. The information input unit 3 is an interface or the like to be connected to the X-ray diffractometer 11 and the input device 13. The analysis unit 2 is configured to acquire the X-ray diffraction data from the information input unit 3 and subject the X-ray diffraction data to preprocessing to generate a powder diffraction pattern of the sample. In this case, the preprocessing refers to processing, for example, smoothing of data, or removal of a K2 component. The powder diffraction pattern generated by the analysis unit 2 is input and stored in the storage unit 5. The X-ray diffractometer 11 may include an analysis unit (data processing unit), and the analysis unit of the X-ray diffractometer 11 may subject X-ray diffraction data to be measured to preprocessing to generate a powder diffraction pattern of the sample and output the powder diffraction pattern of the sample to the information input unit 3 of the quantitative phase analysis device 1. The analysis unit 2 is configured to acquire the powder diffraction pattern of the sample from the storage unit 5 (or the information input unit 3) and perform quantitative phase analysis of crystalline phases contained in the sample based on the powder diffraction pattern to output a weight ratio of the crystalline phases subjected to quantitative phase analysis to the information output unit 4 as an analysis result. The information output unit 4 is an interface or the like to be connected to a display device 12 and is configured to output the weight ratio of the crystalline phases to the display device 12. The display device 12 is configured to display the analysis result of the quantitative phase analysis.
[0048]
[0049] [Step S1: Powder Diffraction Pattern Acquisition Step]
[0050] A powder diffraction pattern of a sample is acquired (S1: powder diffraction pattern acquisition step). The powder diffraction pattern of the sample is stored in the storage unit 5. Alternatively, as described above, the X-ray diffractometer 11 may include an analysis unit (data processing unit) and subject X-ray diffraction data on a sample to be measured to preprocessing to generate a powder diffraction pattern of the sample. Then, the X-ray diffractometer 11 may output the powder diffraction pattern of the sample to the information input unit 3 of the quantitative phase analysis device 1. The analysis unit 2 of the quantitative phase analysis device 1 acquires the powder diffraction pattern of the sample from the storage unit 5 (or the information input unit 3). In the powder diffraction pattern, a horizontal axis represents a diffraction angle 2 indicating a peak position, and a vertical axis represents a spectrum indicating an intensity of a diffraction X-ray. In this case, the diffraction angle 2 is an angle formed by an incident X-ray direction and a diffraction X-ray direction. The X-ray diffraction data on the sample measured by the X-ray diffractometer 11 may be input to the information input unit 3 or stored in the storage unit 5. In this case, the analysis unit 2 acquires the X-ray diffraction data on the sample from the information input unit 3 or the storage unit 5 and subjects the X-ray diffraction data on the sample to preprocessing to generate the powder diffraction pattern of the sample.
[0051] [Step S2: Qualitative Phase Analysis Result Acquisition Step]
[0052] Information on a plurality of crystalline phases contained in the sample is acquired (S2: qualitative phase analysis result acquisition step). The analysis unit 2 identifies crystalline phases based on the positions and the intensities of a diffraction lines (peaks) of the powder diffraction pattern of the sample acquired in the step S1. That is, the analysis unit 2 acquires information on a plurality of crystalline phases contained in the sample by qualitative phase analysis. In this case, the information on the crystalline phases contains chemical compositions thereof, information on polymorphism when the crystalline phases have polymorphism having different crystal structures, and a plurality of peak positions of powder diffraction patterns of the crystalline phases. The information may further contain intensities at the plurality of peak positions of the powder diffraction patterns of the crystalline phases.
[0053] Based on the peak positions and the peak intensities of the powder diffraction pattern of the sample acquired in the step S1, the analysis unit 2 acquires information on the plurality of crystalline phases contained in the sample by subjecting the sample to qualitative phase analysis. However, the present invention is not limited thereto, and the information input unit 3 may acquire the information on the plurality of crystalline phases contained in the sample, which is a result of the qualitative phase analysis of the sample, from the input device 13.
[0054] [Step S3: Fitting Function Acquisition Step]
[0055] A fitting function for each of a plurality of crystalline phases contained in the sample is acquired (Step S3: fitting function acquisition step). A user determines to execute fitting by using one fitting function, which is selected from the group consisting of first to third fitting functions, for the powder diffraction pattern of each of the plurality of crystalline phases based on information on the powder diffraction pattern of the sample acquired in Step S1 and information on the plurality of crystalline phases acquired in Step S2. The user uses the input device 13 to input a fitting function to be used for each of the plurality of crystalline phases. The analysis unit 2 acquires a fitting function for each of the plurality of crystalline phases, which is input to the input device 13 by the information input unit 3.
[0056] Now, a description is given of the first to third fitting functions. When the powder diffraction pattern y(2) of the sample can be regarded as superimposition of a background intensity y(2).sub.back and the powder diffraction pattern y(2).sub.k of each of the K crystalline phases, the powder diffraction pattern y(2) of the sample is represented by the following Expression 10.
y(2)=y(2).sub.back+.sub.k=1.sup.Ky(2).sub.k(10)
[0057] The powder diffraction pattern y(2).sub.k of each crystalline phase can be represented by various forms, and serves as the fitting function. An integrated intensity acquired through whole-powder pattern decomposition based on a Pawley method is used as the first fitting function, and is represented by the following Expression 11.
y(2).sub.k=.sub.jI.sub.jkP(2).sub.jk(11)
[0058] In this expression, P(2).sub.jk represents a normalized profile function describing the profile shape. A function defined on an interval [, +], for example, a pseudo-Voigt function, is used as P(2). However, in actuality, the function may be considered to have a value only around the peak position of each diffraction line.
[0059] An integrated intensity obtained through observation or calculation, which is input from the outside, is used as the second fitting function, and is represented by the following Expression 12.
y(2)=.sub.jSc.sub.kI.sub.jkP(2).sub.jk(12)
[0060] In this expression, Sc.sub.k represents a scale factor, and is defined by I.sub.jk=Sc.sub.kI.sub.jk. {I.sub.jk}, which is a set of integral intensities, may be a set of integral intensities separately measured (or calculated) for a single-phase sample of the k-th crystalline phase, or may be a function of the crystal structure parameter. In fitting, the set {I.sub.jk} of integral intensities is fixed, and the scale factor Sc.sub.k is refined instead.
[0061] A profile intensity obtained through observation or calculation, which is input from the outside, is used as the third fitting function, and is represented by the following Expression 13.
y(2).sub.k=Sc.sub.ky(2).sub.k(13)
[0062] In this expression, Sc.sub.k represents a scale factor similarly to the second fitting function. y(2).sub.k may represent a profile intensity separately measured (or calculated) for a single-phase sample of the k-th crystalline phase, and may be calculated at the time of fitting based on the crystal structure parameter. In fitting, the profile intensity y(2).sub.k is fixed, and the scale factor Sc.sub.k is refined instead.
[0063] The user determines which one of the first to third fitting functions is to be used to perform fitting for the powder diffraction pattern of each of the plurality of crystalline phases based on the information on the powder diffraction pattern of the sample and the information on the plurality of crystalline phases. The first fitting function is desired to be selected when the crystalline phase has a high crystallinity and the symmetry of the crystal is relatively high. In contrast, the second or third fitting function is desired to be used for a crystalline phase having a low symmetry and exhibiting a complex diffraction pattern including a large number of peaks. In particular, when the crystallinity is low and the peak profile has collapsed, the third fitting function may be used.
[0064] [Step S4: Whole-Powder Pattern Fitting Step]
[0065] The fitting function for each of the plurality of crystalline phases acquired in Step S3 is used to execute whole-powder pattern fitting for the powder diffraction pattern of the sample to acquire a result thereof (Step S4: whole-powder pattern fitting step). In this case, the fitting function to be used for whole-powder pattern fitting is Expression 10, and the fitting function y(2).sub.k of the k-th crystalline phase described in Expression 10 is any one of the first to third fitting functions.
[0066] When the first fitting function is used, the parameter of a model required for calculating the shape of the profile is: (a) a parameter for determining a full width half maximum (FWHM), (b) a parameter for determining the shape of the profile, and (c) a lattice constant of the k-th crystalline phase. When the first fitting function is used, the initial value of the integrated intensity I.sub.jk is not required.
[0067] When the second fitting function is used, the parameter of a model required for calculating the shape of the profile includes the above-mentioned parameters (a) to (c) similarly to the case of using the first fitting function, and further includes the scale factor and the parameter of the integrated intensity determined in advance. As described above, the parameter of the integrated intensity is fixed in fitting.
[0068] When the third fitting function is used, the parameter of a model required for calculating the shape of the profile includes the scale factor and data on the measured or calculated profile function y(2) from which the background intensity is removed. As described above, the profile function y(2) is fixed in fitting.
[0069] Even when any one of the first to third fitting functions is used, the parameter of the background intensity y(2).sub.back is required as described in Expression 10. Those parameters are optimized by fitting to acquire a result thereof.
[0070] [Step S5: Weight Ratio Calculation Step]
[0071] The weight ratio of the plurality of crystalline phases is calculated based on the fitting result acquired in Step S4 (weight ratio calculation step).
[0072] When the first fitting function is used, Expression 7 is used to calculate the parameter S.sub.k of the k-th crystalline phase based on the integrated intensity I.sub.jk.
[0073] When the second fitting function is used, the following Expression 14 is used to calculate the parameter S.sub.k of the k-th crystalline phase based on the scale factor Sc.sub.k and the integrated intensity I.sub.jk.
S.sub.k=Sc.sub.k.sub.j=1.sup.N.sup.
[0074] When the third fitting function is used, the following Expression 15 is used to calculate the parameter S.sub.k of the k-th crystalline phase.
S.sub.k=Y.sub.k=Sc.sub.k.sub.2.sub.
[0075] Now, a description is given of derivation of Expression 15. G.sub.jk, which is defined in Expression 2, can be regarded as a continuous function G(2) with respect to the diffraction angle 2. When G(2) is multiplied by both sides of Expression 11 and integrated over the finite 2 range of [2.sub.L, 2.sub.H], the following Expression 16 is obtained. The integrated value obtained through this integration is represented by Y.sub.k. It suffices that the 2 range contains a sufficient number of diffraction lines required for performing quantitative phase analysis as described above.
Y.sub.k=.sub.2.sub.
[0076] As described above, the profile function P(2).sub.jk may be considered to have a value only around the peak position of each diffraction line, and G(2) to be multiplied by the profile function P(2).sub.jk may be considered to have a constant value within that range without producing a significant difference. Further, the profile function P(2).sub.jk is normalized so as to satisfy P(2).sub.jkd(2)=1. Thus, Y.sub.k may be considered to be equal to the parameter S.sub.k given in Expression 7. Thus, Expression 15 is derived based on Expression 13 and Expression 16.
[0077] As described above, the substance parameter a.sub.k is acquired when the chemical composition of the crystalline phase is determined by qualitative phase analysis. Further, even when the sample contains an uncertain crystalline phase (uncertain crystalline phase) whose chemical composition is uncertain, the substance parameter a.sub.k of such a substance may be estimated. Thus, the weight factor W.sub.k of the k-th crystalline phase is calculated by using the substance parameter a.sub.k and the parameter S.sub.k acquired from the fitting result acquired in Step S4.
[0078] Thus, it is possible to calculate the weight ratio of the plurality of crystalline phases contained in the sample by using the weight factor W.sub.k. Further, it is possible to calculate the weight fraction w.sub.k of the k-th crystalline phase by using Expression 8 or the IC formula (Expression 9).
[0079] Now, a description is given of an example of performing quantitative phase analysis for a mixture sample by using a quantitative phase analysis method according to the embodiment.
[0080]
[0081] In Step S1, the observed powder diffraction pattern of the sample is acquired.
[0082] In Step S2, information on the plurality of crystalline phases contained in the sample is acquired. As shown in
[0083] In Step S3, the fitting function for each of the plurality of crystalline phases contained in the sample is acquired. In this case, as shown in
[0084] Quartz being the crystalline phase A has a relatively high symmetry of the crystal structure. Thus, even when there are other components at the time of whole-powder pattern decomposition, fitting is easy to perform, and hence the user selects the first fitting function as the fitting function to be used for the crystalline phase A. In this case, the integrated intensity parameters {I.sub.jA} are all optimized at the time of fitting.
[0085] Albite being the crystalline phase B is one type of feldspar, which is a triclinic crystal system having a low symmetry of the crystal structure. For example, 810 diffraction lines are known to exist within a range in which 2 takes a value of from 5 to 80. The user selects the second fitting function as the fitting function to be used for the crystalline phase B. The set {I.sub.jB} of observed integral intensities, which is obtained by whole-powder pattern decomposition already performed for the sample containing Albite as a single component, is acquired. In this case, the set of observed integral intensities acquired in advance are set as the integrated intensity parameters {I.sub.jB}, and are multiplied by the scale factor for use in calculation of the profile intensity.
[0086]
[0087] Kaolinite being the crystalline phase C is one type of clay mineral, and is formed by weathering and degeneration of feldspar, for example. Kaolinite is known to have a low crystallinity. Thus, among a plurality of peaks included in the observed powder diffraction pattern of Kaolinite, there are a large number of adjacent peaks superimposed onto each other because those peaks become wider, with the result that the shape of the diffraction profile is not clear. This causes a difficulty in applying the whole-powder pattern decomposition method, which is based on the clear shape of the diffraction profile. Thus, the user selects the third fitting function as the fitting function to be used for the crystalline phase C. The background intensity is removed from the observed powder diffraction pattern of a sample containing Kaolinite as a single component to be acquired as the observed profile intensity y(2).sub.c. In this case, the observed profile intensity acquired in advance is set as the observed profile intensity y(2).sub.c, and the observed profile intensity y(2).sub.c is multiplied by the scale factor to be used for calculating the profile intensity.
[0088]
[0089] In Step S4, the fitting function for each of the plurality of crystalline phases is used to execute whole-powder pattern fitting for the powder diffraction pattern of the sample and acquire a result thereof. As shown in
[0090] In Step S5, the weight ratio of the plurality of crystalline phases is calculated based on the fitting result acquired in Step S4. In this case, the weight fraction w.sub.k is acquired, and the result is shown in
[0091] In the above, a description has been given of the quantitative phase analysis method according to the embodiment. With the quantitative phase analysis method according to the embodiment, at the time of fitting using the first to third fitting functions, it is possible to execute fitting as long as the input parameter can be determined, and it is possible to implement quantitative phase analysis as long as the substance parameters a.sub.k of the plurality of crystalline phases contained in the sample can be determined or estimated by qualitative phase analysis.
[0092] In the quantitative phase analysis method according to the embodiment, when the sample contains the plurality of crystalline phases, the fitting function for each of the plurality of crystalline phases may be any one of the first to third fitting functions. Even when the fitting function for the powder diffraction pattern of the sample includes any one of one type of, two types of, and three types of fitting functions among the first to third fitting functions, it is possible to execute the whole-powder pattern fitting for the powder diffraction pattern of the sample, and to calculate the weight ratio of the plurality of crystalline phases.
[0093] In particular, in the related-art quantitative phase analysis method, it is extremely difficult to execute whole-powder pattern fitting for the powder diffraction pattern of the sample by using two or more types of fitting functions at the same time among the first to third fitting functions. In contrast, in the quantitative phase analysis method according to the embodiment, two or more types of fitting functions may be selected among the first to third fitting functions. The two or three types of fitting functions can be used at the same time to execute whole-powder pattern fitting, to thereby achieve a significant effect.
[0094] In the quantitative phase analysis method according to the embodiment, the user determines which one of the first to third fitting functions is appropriate for the fitting function corresponding to each of the plurality of crystalline phases contained in the sample, based on the diffraction pattern of the sample. However, the configuration is not limited thereto. The analysis unit 2 may automatically perform determination by itself, and determine the fitting function for each of the plurality of crystalline phases to acquire a result thereof.
[0095] In the embodiment, the powder diffraction pattern of the sample to be subjected to whole-powder pattern fitting includes the background intensity, and thus the fitting function also includes the background intensity. However, the configuration is not limited thereto. The background may be removed in the preprocessing, and the powder diffraction pattern of the sample to be subjected to whole-powder pattern fitting may not include the background intensity. In this case, the fitting function does not include the background intensity.
[0096] In the above, the quantitative phase analysis device, the quantitative phase analysis method, and the quantitative phase analysis program according to the embodiment of the disclosure have been described. The disclosure is not limited to the above-mentioned embodiment and can be widely applied. For example, the powder diffraction pattern in the above-mentioned embodiment is obtained by X-ray diffraction measurement. However, the powder diffraction pattern is not limited thereto and may be obtained by another measurement, for example, neutron diffraction measurement. Further, various approximations, such as the determination of diffraction lines included in the powder diffraction pattern and the distribution of intensities of diffraction lines that are superimposed onto each other or close to each other, are considered, as required. In the quantitative phase analysis method in the above-mentioned embodiment, the weight ratio of the plurality of crystalline phases is calculated, but another quantitative ratio, for example, a molar ratio, may be calculated based on such weight ratio.