Bingjun Wu(吳秉駿), Jingkai Xia(夏經鎧), Shuo Zhang(張碩), Qiang Fu(傅強), Hui Zhang(章輝),Xiaoming Xie(謝曉明),4,5, and Zhi Liu(劉志),,§
1State Key Laboratory of Functional Materials for Informatics,Shanghai Institute of Microsystem and Information Technology,Chinese Academy of Sciences,Shanghai 200050,China
2Center for Transformative Science,ShanghaiTech University,Shanghai 201210,China
3School of Physical Science and Technology,ShanghaiTech University,Shanghai 201210,China
4Center for Excellence in Superconducting Electronics,Chinese Academy of Sciences,Shanghai 200050,China
5University of Chinese Academy of Sciences,Beijing 100049,China
Keywords: x-ray emission spectra and fluorescence, superconducting transition-edge sensor, rare earth elements,chemical composition analysis
The accurate analysis of the composition of each element in a compound is an essential and indispensable aspect of material analysis research.The chemical properties of different compounds,in the particular case of rare earth oxides(REOs),are determined by the specific composition and bonding behavior of different rare earth elements (REEs), highlighting the importance of precise determination of elemental composition in these materials.[1,2]And further in the practical research and application of the rare earth functional material,the trace level doped REEs with slightly different chemical compositions will lead to changes in the properties,which can be mainly attributed to the particularity of the electronic configuration of REEs.[1,3]In order to accurately determine the chemical composition of trace REEs in materials under laboratory conditions,it is necessary to adopt advanced analytical techniques like the methodology of x-ray fluorescence(XRF),which is crucial for gaining insight into more potential applications of REOs.[4–6]In XRF, the electron within the sample atom is excited when irradiated by a primary x-ray beam.This produces a hole in the initial electron orbital and results in an unstable atomic state.Then an electron from a higher energy orbital fills the hole and emits characteristic x-ray photons with quantized energy.[5]By collecting the emitted x-ray photons and analyzing the whole fluorescence spectrum,XRF can supply a large amount of information,including the peak energies, line shapes, electron binding energies, energy level widths,fluorescence yields,etc.[7,8]
In the XRF analysis, as the x-ray line energies are distributed in a certain range with a broadened extent between shells of differentZ-value elements, there would possibly be numerous peak overlaps, i.e., weaker spectral lines would be hidden beneath the extension range of spectral lines with higher intensity.[9–11]Such overlapping spectral feature is usually presented in the form of an intrinsic spectral peak broadening overlap or the existence of background spectral shape.[6,12]A better understanding of peak energies and spectral profiles would aid elemental composition extraction through peak-fitting analysis.In particular, for XRF research and analyses involving REEs,L-edge emission lines of high-ZREEs are preferred over the more intenseK-edge emission lines, which necessitate excitation energies beyond the capabilities of most laboratory-based x-ray sources, despite constraints caused by adjacent emission line overlap.[3,13]In addition, due to the limitation of the spectrometer resolution,the measured emission spectra of REEs may also overlap with theK-edge emission lines of 3d transition metals,such as Fe,Co, Ni, etc., which are usually with much higher abundance in natural materials.[14]The energy resolution required for those measurements has exceeded the capabilities provided by traditional elemental analysis instruments or spectroscopies,with the energy dispersive spectrometer (EDS) as one of the representatives.[11,15]Thus, in order to unravel the complicated XRF spectra with overlapping lines and to accurately measure the line energies and profiles for REEs,it is essential not only to fit and analyze the spectral data but also to explore,at a more fundamental level, the development of instruments with higher energy resolution.[16]
We present a newly developed instrument for XRF spectroscopy measurement here,which meets the experimental requirements including high detection efficiency,precise linearity and broad-band response.[16]Our scientific instrument is capable of measuring the absolute energies and line profiles of XRF emission lines,with a calibration process to linearize the energy response of each photon captured and to establish an XRF database of all other, unknown spectral features.[17–20]Our spectrometer used for this work consists of an array of 16-pixel(4×4)superconducting device called transition-edge sensor (TES), and a readout based on superconducting quantum interference device (SQUID) amplifiers.[20–22]Our detector can directly measure the energy of each incident xray photon, resulting in much better spectral resolution compared to a commercial silicon drift detector.The advantages of TES detectors over wavelength-dispersive grating spectrometers include a large collection area, high quantum efficiency and wide energy range,making them well-suited for studying a variety of materials and phenomena.[23,24]These attributes also enable obtaining high-resolution spectra from dilute and damage-sensitive samples, with dilutions as low as 0.001%(10 ppm), and a highP/Bratio for distinguishing weak signals from background fluctuations, which is essential for the analysis of trace elements.[20,25]Especially in the field of rare earth sciences, TES detectors have been proposed for a variety of applications, including high energy-resolution x-ray spectroscopy,material diagnostics,and mass spectrometry for neutral molecules.[11]
We have studied theL-edge emission lines of three adjacent and widely-studied REEs: terbium(Tb,Z=65),dysprosium(Dy,Z=66)and holmium(Ho,Z=67)for digging out their intrinsic information, with the energy calibration work carried out by theK-edge emission lines of the 3d transition metal elements.All theL-edge emission lines of the selected elements,from their lowest energy levelsL3M1(Lλ)up to the highestL2N4(Lγ1),lie in a range from 5.5 keV to 9.5 keV(see Fig.1), which also happens to be within the well-calibrated range of our instrument.[26,27]Thus,the preciseL-edge emission line energies and their full line profiles of each element can be well measured and estimated.

Fig.1.Most of the L-edge emission lines of all 17 REEs are shown,with the certain selected three REEs emphasized.The energies of terbium(Tb,Z=65),dysprosium(Dy,Z=66)and holmium(Ho,Z=67)lie in a range from 5.5 keV to 9.0 keV,with the lowest energy levels L3M1(Lλ)up to the highest L2N4(Lγ1).Data are rebuilt from NIST database.[26,27]
As the core component of the spectrometer, TES plays a key role in converting the energy signal into an electrical signal that can be calibrated and analyzed for statistical processing to obtain the spectrum.[28,29]The whole array of TESs is fabricated on a single silicon wafer.All TESs are operated at a superconducting transition temperature of about 90 mK,where the devices can absorb the energy of x-ray photons and be heated up to a non-superconducting state of higher temperature,and then transfer the energy to the surrounding bath by a weak thermal connection,thereby cooling back to the initial normal state.[16,30,31]
Thus,each TES consists of two components in close thermal contact: a heat absorber and a thermometer.[32]The absorbing layer of bismuth(Bi)with a thickness of about 2.4 μm and an area of about 500 μm×500 μm can convert x-ray photons into thermal energy with an absorption efficiency of 70%at 5.9 keV.The total detection efficiency contributed by each process of x-ray photons with energies from 0 keV–10 keV is shown in Fig.2.[7,20,33]The molybdenum-copper doublelayer film serves as a sensitive resistive thermometer, which can be held in its narrow superconducting transition by Joule self-heating.The composition and size are selected to achieve the required transition temperature and electrical resistance.Examples of the TES signal pulse diagram at various energies are shown in Fig.3, where the amplitude of the pulses indicates the amount of energy deposited.[30]

Fig.2.The total detection efficiency can be calculated by the product of all efficiency contributed to the whole process.Within the energy range of 5 keV–10 keV,the transmission process behaves stronger as the energy increases,while the absorption process of x-ray photons becomes weaker,accomplished by the heat absorber of TESs, which are made of bismuth(Bi)with a thickness of 2.4 μm.The total detection efficiency curve can be used for further calculation of the chemical composition of measured samples.Data are generated with the online tools.[7,33]

Fig.3.The shapes of seven different current pulses of a single TES behave similarly, and the pulse amplitude indicates the amount of energy deposited.By convention, each pulse signal with a current decrease is reversely plotted as a positive relative change of TES current.The absorption of different x-ray photon energies which fully span the calibrated energy range, causes a sharp rise of temperature and resistance, accompanied by a sudden drop in the bias current, followed by a recovering procedure to the initial superconducting state within a few milliseconds.
By the connection of the molybdenum leads,these signals of the TES chips are then transmitted to the cryogenic signal amplifier, which mainly includes SQUID array operating at 4 K,room temperature feedback amplifier circuit,control circuit and cryogenic cables.The room temperature feedback amplifier circuit possesses a bandwidth of 1 MHz.The TES array and the SQUID array are both located inside an adiabatic demagnetization refrigerator(ADR).The vacuum window and a series of thin aluminum filters are adopted for blocking infrared and visible radiation.[34]A schematic diagram of the spectrometer system is shown in Fig.4 below.The signal acquisition and analysis system mainly include analog-to-digital conversion (ADC) board, waveform analysis module and energy spectrum analysis software.[17]After the signal is transmitted to the analysis module, the pulse height is calibrated and fitted by offline analysis software, then the x-ray energy spectrum is finally obtained.After preliminary tests, the detector has an energy resolution of 5.2 eV@5.9 keV.

Fig.4.(a) The schematic diagram of the TES x-ray detector.(b) TES sensors are located in the 0.05 K stage of the ADR, at the forefront of the cryogenic probe and parallel to the window.(c)The entire spectrometer system is integrated within an ADR that contains vacuum stages at various temperatures.
The tested samples are compressed from different REO powders provided by Macklin?Reagent Company,with a purity of 99.99% according to the supplier’s certificate of analysis.These powdered oxides include terbium (III,IV) oxide(Tb4O7), dysprosium (III) oxide (Dy2O3) and holmium (III)oxide (Ho2O3).A precision electronic balance is used to weigh 100 mg of each oxide sample, and a set of tableting die and tableting machine are used to perform tableting samples on them under the same conditions.After preparation,the final compressed samples are all uniform discs with a diameter of 8 mm and a thickness of about 1 mm.Moreover,an evenly mixed sample piece of the three oxide components,fully grinded afterward and compressed under the same experimental conditions,is used for assessing the feasibility of our testing methods.Through accurate calculation of the molecular weight of each oxide,it is ensured that the concentration of each component remains consistent and comparable.
After accomplishing the compressing operation of REO samples,these compressed samples need to be fixed on a specific sample holder in order to control the stability of experimental parameters such as x-ray irradiation position,intensity and reflection traveling path, etc., meanwhile as well ensure the repeatability of the experiment.[12]Therefore, we adopt model designing and machining to manufacture a modularassemblable, sample-replaceable and reusable sample holder using polymethyl methacrylate(PMMA),a type of plastic material with high transparency and strong formability,as shown below in Fig.5.The sample holder can not only meet the requirements of rapid and convenient replacement of the sliced samples,but also realize both reflection and transmission measurement modes.As the compressed samples are steadily placed into the pre-designed groove with a diameter of 8 mm and a depth of about 1 mm inside the sample holder, and the incident angle of x-ray onto the sample is precisely fixed at 45°for maximized collection,the x-ray emitted by various elements of the sample will accurately reach our TES sensors through the designed path.

Fig.5.Both reflection (left, XES) and transmission (right, XAS) measurement modes can be realized by adopting the modular-assemblable,sample-replaceable and reusable sample holder made of PMMA.The sliced samples ought to be placed onto or into the pre-designed groove,so as to ensure the reflected or transmitted x-ray microbeam would accurately reach the TES sensors and finally can be collected and analyzed.
A commercial miniature x-ray tube,Mini-X2 with Au anode from AMPTEK?, is adopted to generate the initial incident x-ray.The electrons are accelerated from a cathode tube with the voltage of up to 15.0 kV towards the Au anode target,ensuring obtaining both bremsstrahlung and characteristic fluorescence emission lines.The current of the tube is set to 20.0 μA,i.e.,the power condition is controlled at 300 mW to stabilize the x-ray photon flux.The bias current of TES is set to 0.95 mA.The total test time per sample is about 16 h for the sufficient valid counts up to about 106.As is briefly explained in the introduction part previously, the incident x-ray hits the surface of the sample,thus stimulating the emitted xray with theL-edge emission lines of our elements of interest included, strikes the TES array, after traveling a distance of about 5 cm.Therefore,we ensure that only photons reflected,scattered, or absorbed and re-emitted by the target reach the sensors.When such XRF photons irradiate directly on the TESs,the energy contained in each photon is converted into a current pulse signal,and then recorded and saved.Through the accumulation of fluorescence signal data for several hours,all the pulse signals are counted according to the frequency distribution of pulse height, and theL-edge fluorescence spectrum of the sample in a wide energy range can be obtained.Afterwards, the pulse heights are accurately converted into energy scale by energy calibration.Then after peak fitting,the center energies,the peak shapes and the area ratio of different peaks can be accurately measured.[24,35,36]

Fig.6.The calibration curve for the XRF spectrum of one of the TES channels.Triangle anchors are obtained from the fitting results of K-edge emission lines of various 3d transition metals of standard samples, from Ti Kα1(4511 eV)to Cu Kβ1(8905 eV).The curve is well fitted according to the quadratic function,then can be utilized for the calibration process of the REOs spectrum afterward.
As the electrical signals collected by the spectrometer needed to be accurately converted into energy-scale counts,it is a crucial step for the energy calibration process of the instrument.For calibration, we adopt samples of various 3d transition metals (Z=22-30) as our standard samples,including titanium (22Ti), vanadium (23V), chromium (24Cr),manganese (25Mn), iron (26Fe), cobalt (27Co), nickel (28Ni),copper(29Cu),and zinc(30Zn).Each of theK-edge emission lines of these 3d transition metals contributes a vital point to the total calibration curve,which happens to meet up with the working energy range of our spectrometer,as shown in Fig.6 below.A quadratic function is selected as the fitting function to reduce the mean square error(MSE)of the calibration curve.
By the calibration process mentioned above, the XRF electrical signals ofL-edge emission lines of all the samples are converted into the absolute energy scale,and thus the spectrum can be obtained.To compare with the energy resolution of SDD detector,the experiment of the mixed sample is taken under similar conditions for TES and SDD, as the results are shown in Fig.7 below.The line shapes of fluorescence spectrum are usually theoretically described as the sum of multiple Lorentz(also known as Cauchy)distributions.Meanwhile,the intrinsic broadening effect of the instrument for each energy peak can be approximated to a Gaussian distribution.Therefore, the final energy spectrum is the convolution of Lorentz distribution and Gaussian distribution, which is also widely accepted as Voigt function.[37]And according to this,the peak fitting process of each calibrated energy spectrum is carried out respectively,and the relevant codes are developed and processed offline,which can achieve a large number of spectrum data batch processing.The peak fitting results of each region of the spectrum are shown in Fig.8 below, as fitted lines of different regions are presented in Figs.8(a)–8(d),respectively.
Since the respective weights of the three REOs of the evenly mixed sample are well calculated according to the ratio of moles of REEs at 1:1:1,their theoretical contribution to the mixing energy spectrum is consistent in probability.And at the same time,due to the difference in cross section and other factors of different element atoms,there may be some slight differences in x-ray absorption and reflection coefficients.Therefore,there will be some differences in the intensity of the multipleL-edge peaks across elements in the energy spectrum of the mixed sample.Further, the total detection efficiency in Fig.2 is also taken into consideration.Then,based on the intensity of each peak, we can establish a coefficient database of intensity contribution of each specific element,as shown in Fig.9 for instance.In addition, as for analyzing the energy spectrum of samples with unknown composition for specific elements,the actual concentration of each element can be obtained by intensity of the fitted peaks of the calibrated energy spectrum.Compared with the accuracy of the element composition analysis methodology based on SDD,our methodology would contribute an advance by one or two orders of magnitude.

Fig.7.The calibrated L-edge emission spectrum of the mixed sample of all REOs.The spectrum obtained by our spectrometer of TES (purple) is compared with the spectrum obtained under similar conditions by SDD detector (blue).The TES spectrum is the average result of all valid channels,with the standard deviation plotted as its margin.Vertical lines with elemental peak labels at the bottom of the figure indicate the energies of the L-edge emission lines of REEs.Data are collected from NIST database.[26,27]

Fig.8.The respective peak fitting results of each region of the average spectrum of all REOs.Original data are purple lines with the standard deviation margin,and fitted lines of different regions are plotted as magenta lines.The Lorentz function is chosen as the fitting function,while the energy resolution is under certain restrictions.(a) The fitting result of the measured spectrum in the energy range of Lλ regions.(b) The fitting result of the measured spectrum in the energy range of Lα and Lη regions,as Lα1 and Lα2 peaks are separated accordingly.(c)The fitting result of the measured spectrum in the energy range of Lβ regions,as each of the peaks of the overlapping parts is also separated accordingly.(d)The fitting result of the measured spectrum in the energy range of Lγ regions,with some of the L-edge peaks contributed by Au anode target ignored.

Fig.9.The intensity comparison result for regions of Tb:Ho and Dy:Ho across part of the TES channels.The intensity ratio can be calculated from the proportional relationship between the fitted areas of different regions.The total detection efficiency is taken into consideration as well.It shows small statistical fluctuation between different channels, while the fluctuation seems to be less obvious in larger areas of the region in Fig.8,i.e.,a more stable proportional coefficient can be obtained.
We have described the XRF spectroscopy measurement of REEs, by the newly developed spectrometer consisting of an array of 16-pixel(4×4)superconducting TESs with an energy resolution of 5.2 eV@5.9 keV.Further, we have developed an analytical methodology for material with REEs to figure out the element composition of each REE and its proportion according to the energy spectrum with higher accuracy.In this paper,we have explained the structures,principles and components of the spectrometer, and have elaborated the detailed process of our experiment.The XRF spectrum of evenly mixed samples of oxide components of terbium, dysprosium and holmium are obtained and well calculated out, with energy peaks fitted respectively and analyzed afterwards.According to the result of chemical component analysis,an XRF coefficient database of fitted spectrum peaks of REEs can be established.These measurements and analyses pose an update for the conventional XRF methodology, which could be extremely valuable for digging out the unrevealed energy lines and processes.In the future, through the measured data of XRF spectroscopic analysis of REOs,we can obtain information of various REEs.This methodology can be implemented into practical analyzing experiments for rare earth functional material, which contains trace level REEs and is sensitive to the slight changes in the doping proportion.We believe it can provide some powerful assistance in some cases of as yet undefined compositional analyses.Additionally, it has the potential to provide references for further applications, such as characterizing rare earth materials, analyzing their compositions, and aiding in resource development.Furthermore, the measurement is conducted using an x-ray tube with extremely low power,which implies that the TES-based spectrometer exhibits a high energy resolution and sensitivity,allowing some experiments that conventionally require synchrotron radiation sources to be conducted in small-scale laboratories, thus reducing the barriers for measurement.
Acknowledgments
Project supported by the National Major Scientific Research Instrument Development Project (Grant No.11927805), the National Key Research and Development Program of China (Grant No.2022YFF0608303),the NSFC Young Scientists Fund (Grant No.12005134),the Shanghai-XFEL Beamline Project (SBP) (Grant No.31011505505885920161A2101001), and the Shanghai Municipal Science and Technology Major Project (Grant No.2017SHZDZX02).