METHOD FOR QUANTITATIVE ANALYSIS OF ATOMIC COMPONENTS OF MATERIALS BY LIBS SPECTROSCOPY MEASUREMENTS

Abstract: A method based on the LIBS technique is described, which allows to determime, without calibration of the measurement system, the concentration of atomic components in solid, liquid and aerial samples. The method comprises: (a) obtaining the plasma temperature (T) of at least one species of the emitted radiation; (b) determining the partition function U (T) of each species of the emitted radiation at the plasma temperature; (c) calculating the concentration value of each species once deducted an experimental factor (F); (d) calculating the concentration value of each non-individuated species by Saha equation; (e) calculating the concentration of atomic components as sum of the corresponding species concentrations; (f) measurement of the concentration of components by eliminating the experimental factor through normalization.

Description: - 1 - Method for quantitative analysis of atomic components of materials by LIBS spectroscopy measurements DESCRIPTION The present invention relates to a method for quantitative analysis of atomic components of materials by LIBS spectroscopy measurements without calibration. Quantitative and real time analysis of the elementary composition of materials is of great interest in many fields, from industrial production to environment diagnostics, from on-line control of product quality to monitoring of industrial exhausts. At present, the most widespread systems of quantitative analysis require sampling of the materials, with a consequent off-line characterisation. Generally, this procedure involves an increase of the times for the acquisition of the necessary data and, therefore, of costs and risks. LIBS (Laser Induced Breakdown Spectroscopy) or LIPS (Laser Induced Plasma Spectroscopy) technique is a technique for quick analysis of sample's elementary constituents, that can be applied in situ and gives results in real time. Current systems based on LIBS relay of calibration curves in order to determine the concentration of the various elements [D.A. Cremers, M.J. Ferris and M. Davies: "Transportable Laser Induced Breakdown Spectroscopy (LIBS) instrument for field-based soil analysis". SPIE Vol. 2835 (190-200) 1996]. The present invention proposes a method based on LIBS technology for determining the concentration of the various elements present in the sample, without use of calibration curves. The method of the present invention and embodiments thereof are defined in the claims. The application of the present method allows: i) to perform quantitative analysis by the LIBS technique, without need of reference samples, and without calibrations; PO7EP99/01842 ii) to reduce measurement costs and times, because the method is self- consistent (the experimental quantities that are necessary to determine the concentrations are absolute quantities that are obtained from the same measure). The invention is described herein, with illustrative and non-limiting purpose, for some applications that are considered as being particularly significant. a) Working principle of the LIBS technique LIBS is based on the analysis of the spectrum of the radiation emitted by the atoms of the sample, when they are excited and ionised by a laser pulse that is properly focused on the target. The beam of a pulsed laser is focused on a sample through an appropriate optical system. If the energy of the laser pulse is sufficiently high, the density of photons in the focal spot is such as to trigger multiphotonic ionisation (even if the phenomenon has a very low probability). The first free electrons are accelerated by the electric field of the radiation and they ionise other atoms by impact; a chain reaction which produces the breakdown, with its characteristic detonation noise and flash of light, is thus generated. The plasma (composed of neutral atoms, ions and free electrons), which is formed in the first instants of the laser pulse, absorbs energy from the laser pulse, thus reaching temperatures of the order of about ten eN (~ 10 5 K). The material contained in the focal spot (characteristic volume - 0.1 mm 3 ) is thus heated and ionised. Once the pulse and heating phase are over, the plasma cools down and the atomic species generate their specific spectrum. By resolving the evolution of plasma through time, it is possible to discriminate various radiation regimes: a) bremsstrahlung radiation; b) emission of the spectrum lines of the elements; c) emission of vibration bands of simple molecules, formed as a consequence of the impacts. - 3 - LEBS analysis is based on the measurement of the lines indicated at point b). To this purpose, the observation of the spectrum must start with a certain delay with respect to the laser pulse, in order to reduce the effects of the presence of bremsstrahlung radiation. The typical working regime thus individuated corresponds to plasma temperatures of the order of 1 eV. The duration of the observation is established on the bases of the characteristics of the apparatus (energy and duration of the laser pulse) in order to maximise the signal/noise ratio and to limit the temperature variations due to the dynamics of plasma. In the case in which during the time interval of observation the condition of LTE (Local Thermal Equilibrium) ^Principles of Laser Plasma " Ed. by G. Bekefi. J. Wiley & Sons, 1976] is verified, the plasma is characterised by its temperature and by the electronic density. The atomic elements that form the plasma are present in different ionisation states. It is called species herein an element having a certain ionisation state (e.g.: neutral magnesium [Mgl] and ionised magnesium [Mgll] are two species of the same element). In the case of LTE, the population of the atomic levels is described by Boltzmann distribution; the ratio between populations of two levels belonging to subsequent ionisation states of a same element is given by Saha-Boltzmann equation: N e N, (z) ^ c 0nx« m10 2"1 ≤ g Z 6ϊ,1-τr 6. ^, -G d - E 1 --^^ L (1) N„ (Z - I) ' n v / τ *y where: N e : electronic density (cm "3 ); N n (Z-l): level n population, of the species with (Z-l) charge; N,(Z): fundamental level population, of the species with (Z) charge; g, z : degeneration of the fundamental level of the species with (Z) charge; g n z" ': degeneration of the n level of the species with (Z-l) charge; E "{Α r l): ionisation energy of the (Z-l) charge species initially at level (ft, I); T e : electronic temperature of plasma. T e , which appears in equation 1, is the temperature of the distribution of the free electrons in the plasma; while T plasma temperature stands for temperature of the bonded electrons. In working condition (LTE), the difference between these temperatures is negligible. The intensity of each spectral line, characteristic of each species a present in the plasma, depends on the concentration of the same species N α and on the population of the starting level, that is on the plasma temperature T, according to the formula: I£ = Ν. g * exι V u (r A k ) ,τ,' A ki (2) where I ki is the number of photons emitted per volume unit and for time unit at the wavelength λ, characteristic of the transistion between the levels k and i of the species a, g k is the degeneration of the upper level, E k the energy of the upper level, k B the Boltzmann constant, A ki the probability of transition between the two levels k and i, U a (T) the partition function of the species a at the temperature T, according to the expression: U fl (T) = ∑g 4 exp^ (3) Tj b) Procedure for the analysis at the base of the proposed method The value of the quantity T is obtained through measurements of the intensity of the radiation emitted by electronic transitions. —λ The measured intensity values, i , are connected with the emitted intensity, I k x i , by a proportion factor F, depending on the experimental apparatus and on the measurement conditions (the spectral response of the —λ experimental apparatus is evaluated in the calculation of the quantity Ik,). Such factor F must be constant during a measurement run, that is the acquisition of all the portions of spectrum that are necessary to the measurement must be taken under the experimental conditions. Therefore the following relation is valid: l ~ i= Flt (4) By utilising the measured values of the intensity of the lines emitted by a same species and the relative spectroscopic data (E fo g k , A h , that are available in literature), it is possible to calculate the temperature of the plasma. By substituting in eq. 4 the expression I kl , from equation 2, and by taking the natural logarithm of both terms the result is: In- r h In EN, (5) SΛ, U a (T) k B T In the LTE condition, the temperature T is the same for all species and U a (T) is determined by the temperature. Equation 5 is therefore the equation of a straight line: y b = m + nx kl (6) where: x kl =E k Vi,= ln (7iN A) m = n (FN β / U a (T)) n = -1 / k B T For each species, the values (x h , y k which characterise each line are reported on the plane (E, ln( gA)) for each species the value of the slope n and of the intercept m is obtained by linear regression. The slope n thus obtained gives the value of the temperature of the plasma: T= -l / k B n. (7) With the value T calculated through equation 7 (possibly averaged on the different species), the value of the U a (T) is determined through equation 3. At this point, by using the intercept value , it is possible to calculate the value of the concentration FN a , that is the product of the concentration of the species a times the experimental factor F. In this way, one proceeds to the measurement of FN a for all the species that have been identified (Note: once T has been found, the measure of a single line will be sufficient to determine FN a ). In working conditions (T ~ 1 eV), only spectra that are emitted by the first two species of each element, that is by the neutral state and by the first ionisation one, are observed: indeed, the relative abundance of atoms in states of ionisation above the first is lower than 10 "3 . Once the plasma temperature and the species concentrations (times the F Factor) have been obtained, the concentration of the elements (hereinafter indicated by C,) can be calculated. The procedure is done according to the following scheme: A) Calculation of the concentrations of the elements times the experimental factor F. For each element the following cases can be encountered: i) the concentrations FN a for both species are available: then the concentration of the element is obtained through the sum of the concentrations of the two species. In addition, through eq. 1, it is possible to calculate N e (that in this case is the only unknown quantity). Then, the procedure continues according to point B) ii) the concentration for only one of the two species is available and the value N e is known (this case occurs when there is at least one element for which the concentrations of both species have been measured): the concentration of the element can be obtained by calculating the concentration of the other species, through eq. 1, and by summing the concentrations of the two species. Then, the procedure continues according to point B) iii) the concentrations for only one of the two species is available and - 7 - the value N e is not known (this case occurs when it is possible to measure the lines of a single species for each one of the identified element): in this case N e is calculated as value of convergence of an iterative algorithm based on Saha equation and on the concentrations of the species that have been individuated. That is, a reasonable starting value for N e must be set (in the typical working conditions N e ~ 10 18 electrons/cm 3 ) and then the procedure continues according to point ii) for the determination of the concentrations of the elements, times the F factor; then it is necessary to proceed with the calculation of the real concentrations according to point B). Once the F factor has been eliminated and the concentrations of the elements have been obtained, N e is recalculated as the sum of the concentrations N { of the ionised species (for the neutrality condition of plasma this corresponds to: Tot Λ^ = ∑N. (II), ι'=l where Tot is the number of individuated elements). This new value of N e is reintroduced at point iii) of the algorithm and the cycle is repeated until the value of N e converges. B) Determination of the experimental factor F and calculation of the real concentrations. It is possible to proceed in two ways: i) to utilise an internal reference. When the concentration C ϊ™" nl of an element of the sample is known (e.i. Nitrogen in atmosphere), it is possible to determine the factor F as the ratio between the concentration of the same element, as measured with the LIBS method according to point A), and the known concentration P< ( vFC e e l l e e m m m en t t) 'LIBS ro C "i kknnoowwnn ' / element and then to rescale the other concentrations with the value F, thus - 8 - determined; ii) to normalise each concentration with respect to the sum of all concentrations: ( -1 _ (J ^j ) UBS Q\ ; ~ Tot > y' j j where Tot indicates the number of individuated elements. The concentrations calculated according the methods described in point ii) are affected by a percentage error that is equal to at least the concentration of the elements that have not been individuated. The measurements taken according to the method described in point i) are affected only by a percentage error that is determined by the precision by which the reference is known, but such method is only applicable when an internal reference is available. The algorithm described in points A) and B) is summarised in the •flow chart shown in Fi gure 6. The following is a block scheme describing the procedure for the quantitative analysis herein presented. Determination of the working regime: a) energy of the laser pulse and repetition rate; b) focusing and receiving optics; c) delay and duration of the acquisition gate. d) spectrometer set up Spectra Acquisition at fixed working spot Spectra Analysis: a) individuation of the lines; b) measurement of the wavelength at mid line; c) measurement of the line integral. Qualitative analysis: identification of the species through correspondence between the measurement at mid line and proper database, according to the maximum probability criterion. Quantitative analysis: a) calculation of the temperature T on at least one species; b) calculation of the partition function of each species; c) calculation of the concentration of the species times the experimental factor F; d) calculation of the concentration of the species that have not been individuated by Saha equation; e) calculation of the concentration of the elements as a sum of the concentrations of the species; f) measurement of the concentration of the elements through elimination of the F factor by means of normalisation or internal reference. -10- Description of possible applications: There are described hereinafter different applications of the LIBS technique here described. In all applications the same spectrum analysis procedure that is described in the previous block scheme is utilised. a) Characterisation of metal alloys In the characterisation of metal alloys the interest is usually on the measurement of the concentration of the major components, that is of the components that are present in percentage that is at least 0.1%. In other cases it can otherwise be important to determine the presence of impurities in the alloys, in a quick way, giving priority to measurement sensitivity rather than accuracy. A suitable apparatus is shown in Fig. 1 and consists of a laser radiation source 1, a focusing optic system 2, a first mirror 3, a support 4 for a metal alloy sample 5, a second mirror 6, an optic system 7 for receiving the emitted signal, an optic fibre 8, a spectrometer 9, a timing system 10 for discriminating the spectrum lines from the continuous emission, a computer 11, an electronic control and interfacing system 12 and a software for analysis of data. For the analysis of metal alloys, typical working conditions are: - density of laser radiation flow on the target: ~ 10 8 Wcm "2 - repetition rate of 10 Hz - reception with quartz optics (transparent to UN) - use of a system for the movement of the sample to control the irradiated spot. b) Monitoring of trace elements in the atmosphere In the analysis of traces in the atmosphere the interest is on the measurement of pollutans or other elements that are present in the form of aerosol, submicrometric particulate or vapours, in quantities in the order of ppm. A suitable apparatus is shown in Figure 2 and consists of a laser -11- radiation source 21, a focusing optic system 22 for directing the laser radiation in the gas, an optic system 24, an optic fibre 25, a spectrometer 26, a timing system 27 for discriminating the spectrum lines from the continuous emission, a computer 28, an electronic control and interfacing system 29 and a software for analysis of data. For the monitoring of the elements present in traces in the atmosphere, typical working conditions are: - density of laser radiation flow at the focal spot: ~ 5.10 9 Wcm "2 - repetition rate of 10 Hz (or greater) - reception with quartz optics (transparent to UN). c) Analysis of soils A suitable apparatus for carrying out fast in situ or laboratory analysis of soils with the method according to the invention is shown in Fig. 3 and comprises a laser radiation source 31, a mirror 32 for directing the laser radiation towards on a soil sample 33 on a support 34, a focusing lens 35, a collecting optics 37, an otical fibre 38, an optical signal analysing system 39, a computer 40 and a software for analysis of data. d) Analysis of cultural heritage A suitable apparatus for carrying out the analysis of cultural heritage with the method according to the invention is shown in Fig. 4 and comprises a laser radiation source 41, a frequency multiplier 42, a group of four mirrors 43-46 for selectively directing the laser radiation towards a target 47 either directly (LIPS analysis) or through a frequency multiplier [for LIF (Laser Induced Frequency) analysis], two further mirrors 48-49 for directing the emitted signsl towards a sensor 50, a computer 51 and a software for analysis of data. e) Atomic emission spectroscopy A suitable apparatus for carrying out atomic emission spectroscopy with the method according to the invention is shown in Fig. 5 and comprises a sample 61 on whose surface a plasma 62 is generated by a laser source or -12- any other radiation source such as electric spark, microwaves and so on. The emitted signal reaches a dispersion reticule 63 inside a housing 64 provided with a plurality of slits 65, each faced to a respective photomultiplier 66. The photomultipliers collect the integral line intensity corresponding to specific emission lines of the materials, and these data are processed using the spectral analysis procedure here described.

Claims: -13- C L A I M S 1. Method for quantitative measurement of the concentration of atomic components of solid, liquid or gaseous materials, characterized by comprising: a) spectral analysis of the radiation emitted by a sample excited and ioinized by an electronic radiation source according to the LIBS technique; b) qualitative analysis of the emitted radiation; c) quantitative analysis without calibration of the emitted radiation: 2. Method according to claim 1, characterized in that said quantitative analysis comprises: a) obtaining the plasma temperature (T) of at least one species of the emitted radiation; b) determining the partition function U β (T) of each species of the emitted radiation at the plasma temperature (T); c) calculating the concentration value of each species once deducted an experimental factor (F); d) calculating the concentration value of each non-individuated species by Saha equation; e) calculating the concentration of the . atomic components as sum of the corresponding species concentrations; f) measurement of the concentration of components by eliminating the experimental factor (F) through normalization. 3. Method according to claim 1 or 2, characterized in that said spectral analysis comprises: a) individuation of the spectra lines; b) measurement of the wavelength at a mid line; c) measurement of the line integral. 4. Method according to any claims 1, 2 or 3, characterized in that said qualitative analysis comprises identification of the species through correspondence between the wavelength measurement at mid line and proper -14- database, according to the maximum probability criterion. 5. Method according to any claims 1 to 4, characterized by its application for characterisation of metal alloys. 6. Method according to any claims 1 to 4, characterized by its application for monitoring of trace elements in the atmosphere. 7. Method according to any claims 1 to 4, characterized by its application for in situ or laboratory analysis of soils. 8. Method according to any claims 1 to 4, characterized by its application for analysis of cultural heritage. 9. Method according to any claims 1 to 4, characterized by its application for atomic emission spectroscopy.

Inventors: CIUCCI ALESSANDRO, PALLESCHI VINCENZO, RASTELLI SIMONE, TOGNONI ELISABETTA, CIUCCI, ALESSANDRO, PALLESCHI, VINCENZO, RASTELLI, SIMONE, TOGNONI, ELISABETTA
Applicants: CIUCCI ALESSANDRO, CONSIGLIO NAZIONALE RICERCHE, PALLESCHI VINCENZO, RASTELLI SIMONE, TOGNONI ELISABETTA, CIUCCI, ALESSANDRO, CONSIGLIO NAZIONALE DELLE RICERCHE, PALLESCHI, VINCENZO, RASTELLI, SIMONE, TOGNONI, ELISABETTA
Application year: 1999
Main classification: G01N21
IPC codes: G01N21