Analysis of Hydrogen in Inorganic Materials and Coatings: A Critical Review

: The currently used bulk analysis and depth proﬁling methods for hydrogen in inorganic materials and inorganic coatings are reviewed. Bulk analysis of hydrogen is based on fusion of macroscopic samples in an inert gas and the detection of the thereby released gaseous H 2 using inert gas fusion (IGF) and thermal desorption spectroscopy (TDS). They offer excellent accuracy and sensitivity. Depth proﬁling methods involve glow discharge optical emission spectroscopy and mass spectrometry (GDOES and GDMS), laser-induced breakdown spectroscopy (LIBS), secondary ion mass spectrometry (SIMS), nuclear reaction analysis (NRA), and elastic recoil detection analysis (ERDA). The principles of all these methods are explained in terms of the methodology, calibration procedures, analytical performance, and major application areas. The synergies and the complementarity of various methods of hydrogen analysis are described. The existing literature about these methods is critically evaluated, and major papers concerning each method are listed.


Introduction
The presence of hydrogen in various materials and the ways in which it affects their properties has been the subject of attention of the scientific community for decades. A major area developed in this context since the first half of 20th century was hydrogen in metals [1,2]. Other application areas gradually followed, reflecting the progress in materials technology ever since. The analysis of hydrogen is special in many respects, and it makes sense to treat it separately from other elements. In this context, it is worth mentioning the electron structure of the hydrogen atom; hydrogen has only one electron and does not possess deep inner shells. Therefore, many common methods of material characterization, such as X-ray fluorescence analysis, electron microanalysis (EDX/WDX systems), electron spectroscopies (Auger electron spectroscopy, X-ray photoelectron spectroscopy), and some other cannot be used for hydrogen, for fundamental reasons. Methods for hydrogen analysis are diverse, in every respect, and information about them is scattered through the literature, ranging from analytical chemistry to metallurgy, nuclear science, and other application areas. The aim of this review is to categorize these methods and describe their principles and analytical performance, so that the reader may understand how they work, the kind of information they can provide, and their figures of merit, as well as get an idea of how to address the specific analytical tasks related to hydrogen arising in various areas of science and technology.
Analytical tasks involving hydrogen analysis can be divided into two large groups: (1) bulk analysis concerns the volume of the sample, i.e., where the concentration of hydrogen in the sample as a whole is of interest, and (2) depth profile analysis (depth profiling) in which the distribution of hydrogen is measured as function of the depth beneath the sample surface. It concerns various thin film systems, materials with a modified composition close to the surface, e.g., materials modified by various chemical treatments, ion implantation, etc. This very basic division of the methods of hydrogen analysis is reflected in the structure of this review.

Bulk Analysis of Hydrogen
The earliest analyses of hydrogen concerned metallurgical applications. This is also the area in which the most robust methods of hydrogen analysis emerged. In ferrous metallurgy, hydrogen in steels is an important issue, in connection with unwanted effects such as hydrogen embrittlement [3][4][5], hydrogen flaking [6], and problems caused by hydrogen in welding [7][8][9]. Another area worth mentioning is electrochemical charging of metals by hydrogen and cathodic protection against it [10]. In nonferrous metallurgy, issues related to hydrogen are important in nickel alloys [11], aluminum [12,13], titanium [14], and magnesium [15], as reviewed, e.g., in [16]. Another metal extensively studied in connection with hydrogen is palladium [17] and some other precious metals, e.g., in the role of hydrogen-storage materials. In the production and processing of engineering metals and alloys, it is usually desirable to minimize the amount of hydrogen present in the product.
Typical ranges of hydrogen concentration to be analyzed in metallurgical applications can be assessed on the basis of the plots in Figure 1. The solubility of hydrogen in various metals depends on temperature and crystallographic structure; for example, in austenite (γ-Fe, fcc lattice), it is higher than in ferrite (α-Fe or δ-Fe, bcc lattice). In ferrous metallurgy, welding, casting, etc., it is common to distinguish between "diffusible" and "residual" hydrogen [8,9], whereby the former can be released by diffusion at room or only slightly elevated temperature, unlike the latter. Residual hydrogen concentration depends, in addition to the solubility in the melt, on the conditions of solidification and/or the subsequent heat treatment of the material under study. The content of hydrogen in a material is expressed in various units: atomic or weight fraction (percentage) or the volume of gaseous molecular hydrogen in certain amount of the material, assuming atmospheric pressure and room temperature. These are the ordinate units in Figure 1. This unit, y = cm 3 H 2 at 0.1 MPa per 100 g of material, is linked to the weight percentage of hydrogen in the material, c H , by Avogadro's law.
It is worth noting that weight concentrations of hydrogen are often lower than atomic concentrations by 1-2 orders of magnitude in number, because the atomic mass of hydrogen, M r (H) = 1, is low compared to other elements (e.g., M r (Fe) ≈ 56). This needs to be considered when comparing the sensitivity of analytical methods for hydrogen and other elements, as sensitivity is usually expressed in weight units.
Another issue, relevant to the analysis of hydrogen, is transport properties; hydrogen exhibits a high diffusion coefficient in some materials. Therefore, sample preparation and handling may be critical. In many cases, samples need to be kept at a low temperature before analysis (liquid nitrogen), to prevent losses by diffusion to the sample surface. The diffusion coefficient of hydrogen interstitials in bcc iron metal (α-Fe) is D ≈ 7 × 10 −5 cm 2 ·s −1 at room temperature (25 • C), which yields a characteristic diffusion length d of~0.2 mm (d 2 = 6 Dt) within t = 1 s [3]. On the other hand, in austenite (γ-Fe, fcc lattice), diffusion of hydrogen is slower by orders of magnitude, and austenitic stainless steels are even used as stable reference materials for hydrogen analysis, with an indefinite shelf-life. As an example, the grade 1.4546.9 steel with 0.97 ± 0.05 ppm H (Stahl-H1 by BAM Berlin) can be permanently stored at 4 • C [18].

Figure 1.
Hydrogen solubility in various solids and melts. Note that temperature is expressed in Kelvin. Adapted according to [19]. Copyright (2004) The Japan Society of Applied Physics.

Inert Gas Fusion
Inert gas fusion, (IGF) [20], is a classical method of hydrogen analysis, supported by commercial instruments from several manufacturers. It has become the industry standard in this area and a benchmark for other methods. Sometimes, it is also called gas fusion analysis (GFA); both names are used. Experimental setup is relatively simple; it does not require high vacuum or expensive mass spectrometric detection systems. The sample is placed in a flow of inert gas, typically argon, at atmospheric pressure, and hydrogen present in the sample is released in the form of H2 gas by heating the sample, typically high above the melting point. The amount of hydrogen in the (Ar + H2) gas flow after the furnace is measured by a thermal conductivity (TC) detector [21,22], as the presence of hydrogen in argon substantially enhances the thermal conductivity of the gas. Alternatively, hydrogen in the (Ar + H2) mixture after the furnace can be oxidized to water vapor, H2O, and detected by an infrared absorption detector (IR) [23]; it is based upon the absorption of infrared radiation as it passes through a volume of gas. The advantage of IR detectors is that, by placing a suitable filter into the beam of originally broadband radiation, it is possible to achieve absorption in a narrow region around a specific wavelength. Therefore, the detector can be made selective to a certain molecular gas, such as H2O, without sensing other molecules. Other gases of interest in the IGF analysis also have strong absorption bands in the IR region, such as carbon monoxide, CO, and carbon dioxide, CO2, and can, thus, be detected by IR detectors. Unlike IR detectors, TC detectors are not selective; every component of the gas mixture contributes to its thermal conductivity. Consequently, it is necessary to modify the potentially complex gas mixture originating from the sample so that the TC detector is reached by a simple binary mixture, i.e., the carrier gas + the gas to be measured. In addition to hydrogen, there are other elements in the sample causing the release of gases, typically oxygen and nitrogen. In the IGF method, the sample is melted in a graphite crucible, heated by a high current passing through it, at a power of several kilowatts. Oxygen from the sample reacts with the graphite to form CO and, in some cases, a mixture of CO and CO2. Nitrogen from the sample is released as N2 gas. To facilitate fusion of refractory metals, a flux is usually added to the sample, either pure nickel or pure tin.
The characteristics of the detectors described above and the nature of the gases potentially released from the sample imply different possible configurations of an IGF ana- Hydrogen solubility in various solids and melts. Note that temperature is expressed in Kelvin. Adapted according to [19]. Copyright (2004) The Japan Society of Applied Physics.

Inert Gas Fusion
Inert gas fusion, (IGF) [20], is a classical method of hydrogen analysis, supported by commercial instruments from several manufacturers. It has become the industry standard in this area and a benchmark for other methods. Sometimes, it is also called gas fusion analysis (GFA); both names are used. Experimental setup is relatively simple; it does not require high vacuum or expensive mass spectrometric detection systems. The sample is placed in a flow of inert gas, typically argon, at atmospheric pressure, and hydrogen present in the sample is released in the form of H 2 gas by heating the sample, typically high above the melting point. The amount of hydrogen in the (Ar + H 2 ) gas flow after the furnace is measured by a thermal conductivity (TC) detector [21,22], as the presence of hydrogen in argon substantially enhances the thermal conductivity of the gas. Alternatively, hydrogen in the (Ar + H 2 ) mixture after the furnace can be oxidized to water vapor, H 2 O, and detected by an infrared absorption detector (IR) [23]; it is based upon the absorption of infrared radiation as it passes through a volume of gas. The advantage of IR detectors is that, by placing a suitable filter into the beam of originally broadband radiation, it is possible to achieve absorption in a narrow region around a specific wavelength. Therefore, the detector can be made selective to a certain molecular gas, such as H 2 O, without sensing other molecules. Other gases of interest in the IGF analysis also have strong absorption bands in the IR region, such as carbon monoxide, CO, and carbon dioxide, CO 2 , and can, thus, be detected by IR detectors. Unlike IR detectors, TC detectors are not selective; every component of the gas mixture contributes to its thermal conductivity. Consequently, it is necessary to modify the potentially complex gas mixture originating from the sample so that the TC detector is reached by a simple binary mixture, i.e., the carrier gas + the gas to be measured. In addition to hydrogen, there are other elements in the sample causing the release of gases, typically oxygen and nitrogen. In the IGF method, the sample is melted in a graphite crucible, heated by a high current passing through it, at a power of several kilowatts. Oxygen from the sample reacts with the graphite to form CO and, in some cases, a mixture of CO and CO 2 . Nitrogen from the sample is released as N 2 gas. To facilitate fusion of refractory metals, a flux is usually added to the sample, either pure nickel or pure tin.
The characteristics of the detectors described above and the nature of the gases potentially released from the sample imply different possible configurations of an IGF analyzer, with the potential to analyze H, O, and N simultaneously or in various combinations. Two such configurations, relevant to the analysis of hydrogen, are presented in Figure 2. The reactions utilized to modify the composition of the gas on its path from the furnace to the detector and the corresponding reagents are in Table 1. In the analysis with TC detection (Figure 2a), the gas mixture entering the detector may contain nitrogen in addition to hydrogen. To enable the analysis of this mixture based on thermal conductivity, a column with a molecular sieve is put before the detector. Hydrogen passes through the column much more quickly than nitrogen (similarly to gas chromatography), and it is, thus, possible to separate H 2 from N 2 using a suitable time window in which the TC signal is collected.

Thermal Desorption Analysis
Thermal desorption analysis (TDA) or thermal desorption spectrometry (TDS) [20,[24][25][26] is closely related to IGF and also uses thermal decomposition of the sample. However, instead of a rapid heating above the melting point, the sample is heated slowly, and the rate at which hydrogen is released is registered as a function of temperature. Typical temperature ramp rates are 1-10 °C/min, while the temperature range usually does not exceed 1000 °C. Quartz tube furnaces are used most frequently, with nitrogen as the carrier gas and TC detection. This experimental setup is also called hot extraction. The total amount of hydrogen released throughout the analysis or in a selected temperature region is obtained by integrating the desorption curve. To increase sensitivity, a quadrupole mass spectrometer may be used as a detector (thermal desorption mass spectrometry, TDMS). Typical detection limits based on 1 g sample mass are, thus, in the sub-ppb range (~0.1 ng H) [27].
Applications involve the analysis of "trapped" and "diffusible" hydrogen in steels [24,[27][28][29], hydrogen trapping [25] (see Figure 3), hydrogen-induced second phases, and hydrogen embrittlement mechanisms in Ti alloys [14,24] and hydrogen-storage materials [30]. TDS analysis, together with thermogravimetry, was also used to characterize Ti-H reference materials for the analysis of higher concentrations of hydrogen [31]. For completeness, it is also worth mentioning here the more traditional low-temperature volumetric method of the determination of diffusible hydrogen by its collection at room or slightly elevated temperature under mercury [8,9] or glycerine.  Detection limits are in the sub-ppm range, based on 1 g sample mass, i.e., in the sub-µg range of the absolute amount of hydrogen. The analysis is relative, based on calibration with appropriate reference materials. Alternatively, calibration by known additions of H 2 gas into the carrier gas is also possible. This makes IGF a primary reference method for hydrogen analysis. The analyses are metrologically traceable to fundamental measurement units, in a very wide range of concentrations: from sub-ppm up to the concentrations in hydrides, typically several weight percent of H. IGF analyses are typically performed on dedicated instruments and are considered routine. The typical duration of an IGF analysis is a few minutes.

Thermal Desorption Analysis
Thermal desorption analysis (TDA) or thermal desorption spectrometry (TDS) [20,[24][25][26] is closely related to IGF and also uses thermal decomposition of the sample. However, instead of a rapid heating above the melting point, the sample is heated slowly, and the rate at which hydrogen is released is registered as a function of temperature. Typical temperature ramp rates are 1-10 • C/min, while the temperature range usually does not exceed 1000 • C. Quartz tube furnaces are used most frequently, with nitrogen as the carrier gas and TC detection. This experimental setup is also called hot extraction. The total amount of hydrogen released throughout the analysis or in a selected temperature region is obtained by integrating the desorption curve. To increase sensitivity, a quadrupole mass spectrometer may be used as a detector (thermal desorption mass spectrometry, TDMS). Typical detection limits based on 1 g sample mass are, thus, in the sub-ppb range (~0.1 ng H) [27].
Applications involve the analysis of "trapped" and "diffusible" hydrogen in steels [24,[27][28][29], hydrogen trapping [25] (see Figure 3), hydrogen-induced second phases, and hydrogen embrittlement mechanisms in Ti alloys [14,24] and hydrogen-storage materials [30]. TDS analysis, together with thermogravimetry, was also used to characterize Ti-H reference materials for the analysis of higher concentrations of hydrogen [31]. For completeness, it is also worth mentioning here the more traditional low-temperature volumetric method of the determination of diffusible hydrogen by its collection at room or slightly elevated temperature under mercury [8,9] or glycerine.

Analysis of Hydrogen in Coatings and Thin Films
The analysis of hydrogen in coatings can be divided into two categories: (1) analysis of total amount of hydrogen in a coating and (2) a depth-resolved analysis (depth profiling). The latter case is relevant especially if the depth distribution of hydrogen is uneven. The total amount of hydrogen can then be established by the integration of the depth profile (see Section 3.1.1., Equations (4) and (5)). In principle, for the former task, the methods described in Section 2 can be used, provided that the amount of hydrogen in the coating is high enough and if the contribution of hydrogen present in the substrate is negligible or can be subtracted, on the basis of another analysis, of the plain substrate, without the coating (a "blank" measurement). As an illustration, a 1 μm thick TiH2 layer on a 0.2 mm thick substrate (e.g., stainless steel) can be considered. The hydrogen concentration in the coating is ~4% (in mass units) and the density is ρ(TiH2) = 3.76 g·cm −3 [32], which corre-

Analysis of Hydrogen in Coatings and Thin Films
The analysis of hydrogen in coatings can be divided into two categories: (1) analysis of total amount of hydrogen in a coating and (2) a depth-resolved analysis (depth profiling). The latter case is relevant especially if the depth distribution of hydrogen is uneven. The total amount of hydrogen can then be established by the integration of the depth profile (see Section 3.1.1, Equations (4) and (5)). In principle, for the former task, the methods described in Section 2 can be used, provided that the amount of hydrogen in the coating is high enough and if the contribution of hydrogen present in the substrate is negligible or can be subtracted, on the basis of another analysis, of the plain substrate, without the coating (a "blank" measurement). As an illustration, a 1 µm thick TiH 2 layer on a 0.2 mm thick substrate (e.g., stainless steel) can be considered. The hydrogen concentration in the coating is~4% (in mass units) and the density is ρ(TiH 2 ) = 3.76 g·cm −3 [32], which corresponds to an amount of hydrogen in the coating of 15 µg·cm −2 . For the IGF or TDS analysis, it is feasible to consider a relatively large sample, up to say 10 cm 2 , which, for a one-side coverage by the coating, means 150 µg of hydrogen. This is more than by two orders of magnitude more than the detection limit of IGF or TDS; hence, it is high enough for accurate analysis (see Sections 2.1 and 2.2). To keep the amount of hydrogen in the substrate comparable (same) with that in the coating, the hydrogen concentration in the substrate should be (4/200) × ρ(TiH 2 )/ρ(steel) ≈ 0.0095% = 95 ppm, while ρ(steel) = 7.9 g·cm −3 . Hydrogen concentrations in steel around 100 ppm are much higher than what can be expected, provided that the process of making the coating does not itself add a large amount of hydrogen into the substrate. Hence, analyses like this are feasible. If the sensitivity for hydrogen is enhanced by mass spectrometric detection (TDMS), this approach should also work for coatings with a much lower hydrogen concentration than that in TiH 2 . TDMS analyses of titanium deuteride films on a quartz substrate were described in [33]. The size of the sample may also play a role; for example, a 10 cm 2 sample may seem too large to be sacrificed for a single analysis. An advantage of this is, however, that the IGF and TDA methods are accurate and "absolute", in the sense mentioned in Section 2. This does not necessarily apply to some other methods listed below.
Another fact worth mentioning is that the most common methods of elemental analysis in laboratories involved in materials science are those implemented in electron microscopes: electron microanalysis (EDS/WDS) and electron energy loss spectroscopy (EELS). As neither of these can be used to analyze hydrogen, and because, in many labs, no direct methods of hydrogen analysis are available, the effects of hydrogen are sometimes treated indirectly, without actual analysis. Such an approach was used, e.g., for hydrogen in boron-doped diamond-like carbon layers, in which hydrogen in the coating was sensed by a peak at 188 eV in the X-ray photoelectron spectrum (XPS) of the 1s level of boron, indicating the presence of B-H bonds [34]. Although information about the B-H bonds in the spectrum exists, this can hardly be regarded as a true analysis of hydrogen.

Emission Spectroscopy and Mass Spectrometry
The two emission spectroscopic methods listed below, GDOES (glow discharge optical emission spectroscopy) and LIBS (laser-induced breakdown spectroscopy), are well established for general analytical purposes and are quite common in laboratories involved in solid-state analytical chemistry. The measurements themselves are relatively simple; however, the subsequent interpretation of the data, to establish the concentration of hydrogen, may not be easy, particularly in complex matrices. The emission spectrum of atomic hydrogen has been known for more than a century and was the central point in the development of the quantum theory, during the first third of the 20th century. The following are the lines of choice for the analysis of hydrogen: λ = 121.57 nm (Lyman α) and three lines of the Balmer series, λ = 656.27 nm (H α ), λ = 486.13 nm (H β ), and λ = 434.05 nm (H γ ). The line at 121.57 nm is the most sensitive one but requires a spectrometer designed for the vacuum UV region, as this radiation is strongly absorbed by atmospheric oxygen. Such instruments are common in GDOES, not so much in LIBS. The optical system has to be either evacuated or purged with an inert gas. The second most sensitive hydrogen line is H α , λ = 656.27 nm, in the visible region (red), and it is used in both LIBS and GDOES.
Related to GDOES is glow discharge mass spectrometry (GDMS). It also uses a glow discharge atomization and ionization source, attached to a mass spectrometer, and the analysis is based on the signal of selected ions (see below). Another mass-spectrometric method is secondary ion mass spectrometry (SIMS), in which the ions of the analyzed material are not created in a glow discharge plasma but by bombarding the surface under study using a focused ion beam, under high-vacuum conditions. GDMS and SIMS are more expensive techniques than GDOES and are also less common.

Glow Discharge Optical Emission Spectroscopy (GDOES) and Glow Discharge Mass Spectrometry (GDMS)
In GDOES and GDMS [35,36], the sample acts as cathode in a glow discharge (GD) in argon and is atomized by cathodic sputtering. The sputtered atoms are excited and ionized in the GD plasma and emit characteristic radiation that is analyzed by an optical spectrometer (GDOES). In GDMS, the ions are extracted from the plasma and analyzed by a mass spectrometer. The sample is atomized continuously, layer-by-layer, as the erosion crater proceeds into the depth, and this makes both methods suitable for depth profiling, provided that the signal is registered as function of time. Total depths that can be analyzed range from tens of nanometers to >100 µm. Typical sputtering rates for metallic matrices arẽ 10-100 nm/s and the total duration of the analysis depends on the thickness of the layer under study. Nonmetallic matrices such as nitrides, oxides, and carbides have typically lower sputtering rates than metals. Depth resolution is depth-dependent, typically at~10% of the depth below the original sample surface. GDOES and GDMS depth profiles, not quantified and involving hydrogen, are shown in Figure 4a [37] and Figure 5. A sufficiently high sampling rate, necessary for GDMS depth profiling of several elements in very thin layers, is possible by using a time-of-flight (TOF) mass spectrometer [38].
A crucial point in the analysis by GDOES and GDMS is quantification, i.e., the conversion of raw data (qualitative depth profiles of signal intensity versus time of sputtering) into the composition-versus-depth (quantitative) profiles. Both GDOES and GDMS are used as relative methods; quantification is based on empirical a priori information acquired by analyzing suitable reference samples. Various calibration models and quantification schemes have been developed for both GDOES [35,39] and GDMS [40,41].
In GDOES, the approximation of matrix-independent emission yields (MIEY) is widely used, under which the intensity I E , λ of an emission line λ of an element E is proportional to its concentration c E,M in the matrix M analyzed, as well as to the sputtering rate q M of that matrix.
where R E,λ is a constant called the emission yield, expressing the relative sensitivity of the line λ in the specific instrument and under the specific discharge conditions used. Emission yields are treated as constants, independent of the matrix analyzed. Equation (2) shows that the quantities controlling the spectral response are interlinked, as the sputtering rate q M depends on the matrix (see the plot in Figure 6). Hence, in a qualitative profile, a higher intensity of the line is not necessarily caused by a higher concentration of the respective element. The abscissa of the profile may also be distorted; the time needed to pass different layers is not only proportional to their thickness, but also inversely proportional to their sputtering rates. The effects of a variable sputtering rate may be very substantial; for example, the sputtering rate of diamond-like carbon is~30 times lower than that of iron [42]. This calibration model, Equation (2), allows for accurate quantitative multi-matrix analyses, provided that suitable reference materials with a known composition and known sputter rates are available. Another prerequisite is that virtually all elements present at significant concentrations are measured. The achievable accuracy is well comparable with best methods routinely used for bulk compositional analysis (e.g., X-ray fluorescence, XRF). The presence of hydrogen, however, adds to the complexity, as hydrogen affects GD excitation and ionization of other elements. Such effects were first described in chromium [43] and have been extensively studied in both GDOES [44,45] and GDMS [46]. Various correction schemes have been proposed to handle such situations, e.g., in the analysis of N-implanted Ti [47]. Accurate quantitative analysis of hydrogen by GDOES usually requires the calibration and quantification to be matrix-matched and supported by suitable reference materials. They include bulk standards, such as a hot-pressed powder of TiH 2 mixed with Cu [48] and hydrogenated Zr-Nb alloy [49], or various coatings, e.g., electrodeposited Zn(H), magnetron-sputtered WC(H), or amorphous hydrogenated silicon, a-Si:H, produced by chemical vapor deposition [37,50].     If there is some hydrogen in the coating but not in the substrate, the GDOES analysis can be verified and linked to an "absolute" method, such as IGF or TDS, as follows: quantification of an unknown sample, based on Equation (2), results in the first step in the calculated concentrations of the elements present, cE,M(t), and calculated sputtering rate, If there is some hydrogen in the coating but not in the substrate, the GDOES analysis can be verified and linked to an "absolute" method, such as IGF or TDS, as follows: quantification of an unknown sample, based on Equation (2), results in the first step in the calculated concentrations of the elements present, c E,M (t), and calculated sputtering rate, q M (t), as functions of time, t. The conversion of the time coordinate t into depth x, x = x(t), is based on the definition of the sputtering rate (the sputtered mass-per-time).
where S is the area of the erosion crater, and ρ M (x) is the density at the depth x. The depth coordinate x can then be obtained by integration of Equation (3), provided that the density ρ M is known as function of the composition. There is, however, no universal formula expressing the density of a material as function of its composition; hence, the resulting profile as function of x is only approximate. The total mass of hydrogen in the sputtered area of the coating, m H , is By substituting from Equation (3), this can be rewritten as i.e., the integration over depth x can be replaced by integration over time, while the density is implicitly included in the sputtering rate q M (t). Equation (5) gives, therefore, an accurate correspondence between the measured depth profile and the total amount of hydrogen m H . m H can then be compared with the amount of hydrogen as obtained by IGF or TDS. GDOES is a rapid and affordable method, suitable for obtaining basic semiquantitative information about hydrogen in a coating or for comparison of hydrogen depth profiles in a series of similar samples that differ, e.g., by the conditions of preparation. A flat sample is needed, with an area exceeding the dimensions of the sealing o-ring of the GD source (typical diameter of several mm). The results of GDOES can eventually be further refined by another method or by comparison with GDOES data from matrix-specific standards, created specifically for the application under study. GDOES analyses of hydrogen described in the literature include, in addition to the already mentioned applications, hydrogen in various thin-film and multilayer systems [51], e.g., in diamond-like carbon films [42,52]. The former paper is also a comparison with ERDA (elastic recoil detection analysis; see Section 3.2.1). Other applications involve Si/C anodes of Li-ion batteries [53], analysis of deuterium in W [54] and Ti, and research into plasma-facing materials for thermonuclear fusion [55]. In that study, integrated GDOES depth profiles were compared with analyses by TDS, as mentioned above. Compared to GDOES, much fewer analyses of hydrogen have been reported using GDMS. They include, e.g., analysis of H in gold layers [56] or a-Si:H solar cells [38].

Laser-Induced Breakdown Spectrometry (LIBS)
In laser-induced breakdown spectrometry (LIBS) [57], the sample is irradiated by a laser shot which ablates a small amount of sample material, generating a plasma plume. Its light emission is collected and analyzed by an optical spectrometer. Unlike GDOES, the small diameter of the focused laser beam enables spatially resolved analysis [58]. However, the performance of LIBS in accurate multi-element and multi-matrix analyses is inferior to that of GDOES. A simple but still accurate semi-empirical calibration model, similar to the MIEY approximation in GDOES, does not exist in LIBS, because of the much more complex physics involved in the atomization (ablation) and excitation processes. There have been attempts to interpret the LIBS spectra on the basis of the concepts of fundamental plasma physics, using plasma parameters such as excitation temperatures (calibration-free (CF) LIBS [59][60][61]). Their analytical outcome, however, has so far been rather semi-quantitative, with few exceptions, such as the analysis of glass [62]. Best results can be achieved with matrix-specific calibrations of hydrogen in simple matrices (e.g., a homogeneous matrix consisting of two or few elements only). This is how most of the hydrogen analyses mentioned below were done. A rudimentary depth profiling by LIBS is possible through successive analyses at the same spot [63]; however, the depth range is limited compared to GDOES at both ends, both for thin and thick coatings. Furthermore, only a coarse depth sampling is possible, as the depth step is equal to the depth of the ablation crater resulting from a single laser shot.
Applications of the analysis of hydrogen by LIBS involve hydrogen and deuterium analysis in plasma-facing materials for thermonuclear fusion devices [60,61,63,64], materials for nuclear fission reactors, such as Zircalloy [65], hydrogen in weldments [66], hydrogen-storage materials (MgH 2 ) [67], isotopic analysis of hydrogen in titanium [68,69], and even hydrogen analysis in H-bearing minerals on Mars by the rover Curiosity [70]. Isotopic analysis is traditionally a domain of mass spectrometry; however, the isotopic shift in the hydrogen spectrum also allows distinguishing between hydrogen and deuterium by optical emission. The wavelength difference between H α and D α lines is 0.179 nm (see Figure 7), and this is sufficient for the analysis if a spectrometer with adequate resolution is used [71,72]. In the analysis of H and D in beryllium layers, intensities of only partially resolved D α and H α lines measured with the spectral resolution of~0.1 nm were established by deconvolution, assuming a Gaussian apparatus function and a Lorentzian component caused by Stark broadening [64]. It is worth mentioning that the LIBS signal of hydrogen is enhanced if the analysis is done in helium, typically at a pressure of several torr [72,73]. The likely mechanism of this effect is Penning ionization of hydrogen by the He I, 2s 3 S 1 , and 2s 1 S 0 metastables (19.82, 20.62 eV), followed by recombination of the H + ions in the plasma, thereby populating the excited atomic levels from which the analytical hydrogen lines originate [73].

021, 2,
shift in the hydrogen spectrum also allows distinguishing bet ium by optical emission. The wavelength difference between (see Figure 7), and this is sufficient for the analysis if a spectr lution is used [71,72]. In the analysis of H and D in berylliu partially resolved Dα and Hα lines measured with the spectra established by deconvolution, assuming a Gaussian apparatu component caused by Stark broadening [64]. It is worth men of hydrogen is enhanced if the analysis is done in helium, typi torr [72,73]. The likely mechanism of this effect is Penning io He I, 2s 3 S1, and 2s 1 S0 metastables (19.82, 20.62 eV), followed ions in the plasma, thereby populating the excited atomic leve hydrogen lines originate [73].

Secondary Ion Mass Spectrometry (SIMS)
Secondary ion mass spectrometry (SIMS) [74,75] provides elemental depth profiles, elemental mapping, and isotopic analysis. The sample surface is sputtered in vacuum with a beam of primary ions (usually O 2 + or Cs + ), at an energy in the keV range; a small fraction of the sputtered sample material is ionized, and these secondary ions are extracted and analyzed using a mass spectrometer: quadrupole, magnetic sector, or time-of-flight (TOF). Compared to GDOES, SIMS is a much more sensitive technique; SIMS depth profiles may span over several orders of magnitude of the concentrations of some elements. Furthermore, unlike GDOES, elemental mapping is possible, either by direct imaging using ion optics in a magnetic sector instrument or by scanning a finely focused primary beam. Experimentally, SIMS is a more demanding technique than GDOES, as it requires ultra-high-vacuum conditions, especially when analyzing hydrogen, and it is also much more expensive. Unlike GDOES, SIMS is not used for a general bulk analysis of medium and high concentrations of different elements, because of poor accuracy. Typical sputtering rates in SIMS analyses are much lower than in GDOES, typically by a factor of 10 2 , which also affects the duration of analysis. Depth profiling by SIMS is typically by two orders of magnitude slower than GDOES, for the same layer thickness. The same applies to typical analyzed depths. Unlike GDOES, SIMS can be operated in a "static" mode, providing information about topmost atomic layers on the surface. On the other hand, with GDOES, it is easy to analyze layers as thick as~100 µm, whilst typical depths easily accessible by SIMS are in the range of few hundreds of nanometers only.
The methodology of hydrogen analysis by SIMS was described, e.g., in reviews [76,77]. Cs + sputtering and detection by H − secondary ions is recommended, while a detection limit of H in silicon of~2 × 10 18 atoms/cm 3 was reported in [77], which corresponds tõ 40 ppm (atomic) or~1.4 ppm (weight). It is interesting that this performance was achieved 40 years ago [76] and has not substantially improved ever since. It should also be noted that there is a tradeoff among sensitivity, depth resolution, and lateral resolution [75], and the detection limit mentioned above can hardly be achieved in compositional mapping. SIMS measurements only provide qualitative information, as relative secondary ion intensities are different for different elements and depend on the matrix analyzed and instrumental factors. For quantification, therefore, so-called relative sensitivity factors (RSFs) are required to convert the measured secondary ion intensities into atomic concentrations in the investigated matrix. The RSF is defined as where i indicates an impurity isotope, M indicates the matrix isotope chosen for reference, n i is the impurity atom density (atoms/cm 3 ) in the matrix, I is the secondary ion intensity (counts/s), and A represents the corresponding isotope abundance. A great effort was devoted in the past to describing the physics of the secondary ion emission. Various physical models were proposed, such as bond-breaking, electron tunneling, and local thermal equilibrium (LTE) models [78]. However, they have not become the basis for a useful quantification scheme for complex matrices or a calibration-free analysis of different elements in the same matrix, a would-be analogy to CF-LIBS. The RSFs of different elements in silicon differ by several orders of magnitude [79]. Only matrix-matched calibrations with dedicated reference samples provide an acceptable analytical performance. In some well-defined matrices of technical importance, such as pure Si or ZnO, RSFs are established using ion implanted standards with known fluences [77]. An example of this approach is the depth profile analysis of hydrogenated amorphous silicon, a-Si:H, using implanted deuterium as internal standard [80]. Concerning the depth scale, virtually no attempts have been reported in the literature to correct for the variable sputtering rate of different layers. Hence, it is very likely that the depth accuracy in the analysis of multilayer systems by SIMS is worse than in GDOES.
Analytical applications involving the analysis of hydrogen by SIMS are diverse, and it is not the intention here to give a comprehensive listing. In the area of semiconductors, in addition to the already mentioned a-Si:H [80], SIMS has been used in several studies to monitor the distribution and migration of hydrogen in heterojunction solar cells, where an intrinsic hydrogen-rich amorphous layer is used to passivate interfacial defects. On annealing, the initially amorphous film of a-SiC x (p):H crystallizes into nanocrystalline nc-SiC x (p), and hydrogen effuses from that layer. Re-hydrogenation can be performed via hydrogen diffusion from a sacrificial layer of SiN x :H [81]. This is an example of a sequence of several technological steps in which hydrogen plays a major role.
In the area of metals, an example of a sophisticated analysis involving 2D SIMS mapping is the study of the redistribution of deuterium in the AISI 304L austenitic steel subjected to external mechanical load [82]. Another example is the analysis of hydrogen diffusion through the oxide layer formed on nickel base alloys in water at a high temperature and high pressure, employing a deuterium tracer [83]. In this study, the SIMS analyses were verified by comparison with elastic recoil detection analysis (ERDA). Depth profiles of 1 H and 2 D implanted into niobium oxide were reported in [84] (see Figure 8). Analyses by SIMS of hydrogen in materials for nuclear applications involve, e.g., the diffusion of deuterium through a tungsten coating into ferritic steel [85] or through oxide layers formed on Zircalloy-4 and other zirconium alloys [86].
Hydrogen 2021, 2, analysis can be made using either bulk analysis methods, typically IGF or TDS with I mass spectrometric detection (see Section 2), or SIMS, especially if microanalysis or tope analysis is requested, eventually in combination with FTIR (Fourier-transform in red spectroscopy). An example is the analysis of volcanic glasses [87,88]. A high-res tion magnetic sector instrument was needed to resolve 18 OH ions from 19 F. In [87], ma specific SIMS calibrations of hydrogen, based on the 1 H − / 30 Si − ratios, were shown for ba andesite, and rhyolite. Similar calibrations for olivine, pyroxene, and garnet were sented in [88,89]. A review about the applications of SIMS in geochemistry was prov in [90]. Hydrogen analysis in quartz and quartz glasses was described in [91].

Nuclear Methods
The methods described in this section employ high-energy ion beams, in the M region, and they require instruments and facilities common in nuclear physics, suc particle accelerators, high-energy ion beam optics, and the corresponding detection tems. Such experimental setups are expensive and appear to make these techniques tively less commonly used for material characterization, but they are nevertheless c mercially available. They are frequently used to verify, complement, and help qua analyses made by other, more easily accessible methods (see, e.g., [37,61,83,84,92]).

Elastic Recoil Detection Analysis (ERDA)
In elastic recoil detection analysis (ERDA) [93,94], the sample is irradiated at a glan angle by MeV ions and the energy of the recoiling target atoms is measured, in the inci Another application area worth mentioning is the analysis of hydrogen in minerals and rocks. Hydrogen-bearing species of interest in such materials are H 2 O and OH − . Their analysis can be made using either bulk analysis methods, typically IGF or TDS with IR or mass spectrometric detection (see Section 2), or SIMS, especially if microanalysis or isotope analysis is requested, eventually in combination with FTIR (Fourier-transform infrared spectroscopy). An example is the analysis of volcanic glasses [87,88]. A highresolution magnetic sector instrument was needed to resolve 18 OH ions from 19 F. In [87], matrix-specific SIMS calibrations of hydrogen, based on the 1 H − / 30 Si − ratios, were shown for basalt, andesite, and rhyolite. Similar calibrations for olivine, pyroxene, and garnet were presented in [88,89]. A review about the applications of SIMS in geochemistry was provided in [90]. Hydrogen analysis in quartz and quartz glasses was described in [91].

Nuclear Methods
The methods described in this section employ high-energy ion beams, in the MeV region, and they require instruments and facilities common in nuclear physics, such as particle accelerators, high-energy ion beam optics, and the corresponding detection systems. Such experimental setups are expensive and appear to make these techniques relatively less commonly used for material characterization, but they are nevertheless commercially available. They are frequently used to verify, complement, and help quantify analyses made by other, more easily accessible methods (see, e.g., [37,52,61,83,84,92]).

Elastic Recoil Detection Analysis (ERDA)
In elastic recoil detection analysis (ERDA) [93,94], the sample is irradiated at a glancing angle by MeV ions and the energy of the recoiling target atoms is measured, in the incident hemisphere (forward recoil). Target atoms that are lighter than the incident ions recoil with higher energy than that of forward scattered projectile ions and are, therefore, easier to detect. ERDA is, thus, a suitable method for determination and depth profiling of light isotopes in the presence of heavier elements. This is complementary to a related technique called the Rutherford backscattering spectroscopy (RBS) [93], in which it is the energy spectrum of the scattered projectile ions that carries the analytical information. RBS is suitable for analysis of heavier isotopes in a matrix consisting of lighter atoms. The opposite situation would lead to a high background signal and the consequential deterioration of sensitivity (see Figure 9). In ERDA, ion beams from He to Au with energies in the range 0.3-2 MeV/amu are typically used. The yield of recoiled hydrogen atoms is proportional to the concentration of H and can be established on the basis of the intensity of the incident beam and the geometry of the experiment, using the Rutherford recoil cross-section for H [93]. This makes ERDA a quantitative technique, with a typical accuracy of~10% relative. If the recoil originates from the surface, the energy E r of the recoil atom is determined only by kinematics (conservation of mass and energy) and can be expressed as where E p is the energy of the incident ions, and k = 4m p m r m p + m r 2 cos 2 ϕ where m p and m r represent the mass of the projectile and the recoil atom, respectively, and ϕ is the angle of detection with respect to the incident beam.
Heavier atoms than H have higher energy E r as per Equation (7) and appear as separate groups in the energy spectrum. If the recoil originates deeper in the sample, the detected energy of the recoil atoms is lower due to the energy loss of the incoming ion and the recoil atom on their paths inside the sample. Therefore, the energy spectrum of recoil atoms also contains information about their depth distribution, which can be used for depth profiling. Systems with an additional time-of-flight detector (TOF-ERDA) have been introduced to support some depth profiling applications [95]. A typical detection limit of ERDA for hydrogen is about 0.1 atomic percent [94]. An example of an ERDA spectrum of an amorphous carbon film with surface contamination by H and O on silicon substrate is shown in Figure 9b.
Hydrogen in carbon films is a frequent application of ERDA. This involves amorphous carbon [96,97] (see Figure 9) and diamond-like films [52,98,99]. H-implanted a-Si [95] and a-Si, a-Ge [100] are examples where hydrogen analysis by ERDA was used in semiconductors. Analysis of hydrogen and deuterium by ERDA in oxide layers on nickel alloys was described in [83]. Analyses of hydrogen in metallic layers by ERDA have also involved, e.g., Al/Cu/Ag/Au multilayers [101], Mg/Ni films [102,103], and Fe/Ti films [104]. Figure 9. RBS (a) and ERDA (b) spectra of amorphous carbon film on silicon. Impurity levels are H ≈ 1 at.%, Ar ≈ 1 at.%, and O ≈ 2 at.%. In the RBS spectrum, hydrogen is not visible and the peak of carbon (Mr = 12) is superposed on a background caused by silicon (Mr = 28). Reprinted from [97] with the permission of AIP Publishing. 1 × 10 3

Nuclear Reaction Analysis (NRA)
In nuclear reaction analysis (NRA) [105], the sample is bombarded by a beam of accelerated and energy-monochromatized ions that may induce a nuclear reaction with the nuclide in the sample to be analyzed, while one or several particles are emitted, with characteristic energies, and can be used for detection. The most common reaction for hydrogen analysis is the 1 H( 15 N, αγ) 12 C reaction; the sample is bombarded with 15 N ions, which, upon collisions with protons ( 1 H) in the sample, change into carbon 12 C, while α and γ particles are emitted, whereby the latter are used for detection. The reaction yield, YB, determines the sensitivity to detect an areal density of hydrogen nuclei (nH (cm −2 )) in the sample. It is proportional to the flux (Ni (cm −2 ·s −1 )) of incident projectiles and to the energydependent cross-section of the nuclear reaction, σ(Ε) (cm 2 ).
= (9) where subscript H denotes hydrogen ( 1 H) and i denotes the incident ion, i.e., 15 N. This reaction is resonant; its cross-section σ(Ε) has a narrow peak at ER = 6.385 MeV (the next resonance of the 1 H( 15 N, αγ) 12 C reaction occurs at 13.35 eV), with a width of only 1.8 keV.
The ratio of the on-resonance to the nearby off-resonance cross-section is >10 4 . This makes this reaction very suitable for depth profiling; if the incident ions have a higher energy than the resonance energy, they are decelerated on their path through the sample material . RBS (a) and ERDA (b) spectra of amorphous carbon film on silicon. Impurity levels are H ≈ 1 at.%, Ar ≈ 1 at.%, and O ≈ 2 at.%. In the RBS spectrum, hydrogen is not visible and the peak of carbon (M r = 12) is superposed on a background caused by silicon (M r = 28). Reprinted from [97] with the permission of AIP Publishing.

Nuclear Reaction Analysis (NRA)
In nuclear reaction analysis (NRA) [105], the sample is bombarded by a beam of accelerated and energy-monochromatized ions that may induce a nuclear reaction with the nuclide in the sample to be analyzed, while one or several particles are emitted, with characteristic energies, and can be used for detection. The most common reaction for hydrogen analysis is the 1 H( 15 N, αγ) 12 C reaction; the sample is bombarded with 15 N ions, which, upon collisions with protons ( 1 H) in the sample, change into carbon 12 C, while α and γ particles are emitted, whereby the latter are used for detection. The reaction yield, Y B , determines the sensitivity to detect an areal density of hydrogen nuclei (n H (cm −2 )) in the sample. It is proportional to the flux (N i (cm −2 ·s −1 )) of incident projectiles and to the energy-dependent cross-section of the nuclear reaction, σ(E) (cm 2 ).
where subscript H denotes hydrogen ( 1 H) and i denotes the incident ion, i.e., 15 N. This reaction is resonant; its cross-section σ(E) has a narrow peak at E R = 6.385 MeV (the next resonance of the 1 H( 15 N, αγ) 12 C reaction occurs at 13.35 eV), with a width of only 1.8 keV. The ratio of the on-resonance to the nearby off-resonance cross-section is >10 4 .
This makes this reaction very suitable for depth profiling; if the incident ions have a higher energy than the resonance energy, they are decelerated on their path through the sample material and reach the resonance energy at a certain depth. The reaction then proceeds only in a very narrow layer around that particular depth, see Figure 4b. Typical stopping powers S = (dE/dz) for 15 N near the 6.4 MeV range are on the order of~1-4 keV/nm for elemental targets; thus, the typical near-surface depth resolution is between~1 and a few nanometers. At larger depths, depth resolution worsens, which is a phenomenon referred to as straggling [105]. The stopping power for a given matrix can be calculated on the basis of the areal number densities (cm −2 ) of the elements present in the matrix and their stopping cross-sections ε (for 15 N at 6.385 MeV). Those are known and tabulated, e.g., in [105]. The stopping power S is required not only to convert the energy abscissa of the NRA spectrum into the probing depth but also to derive absolute H concentrations from γ-yield measurements at bulk targets. An extensive collection of experimental (dE/dz) values is, e.g., the database in [106]. For analytical interpretation of NRA measurements, dedicated reference materials have also been proposed, e.g., hydrogenated tantalum [107]. The 1 H( 15 N, αγ) 12 C reaction provides good sensitivity, with detection limits around~100 atomic ppm H or even lower, as well as a practical probing range (2-4 µm). There are few other nuclear reactions for the analysis of 1 H [105]; however, the 1 H( 15 N, αγ) 12 C reaction provides generally the best performance and is most widely used. It should also be noted that two NRA systems for depth-resolved two-dimensional hydrogen mapping have been developed in Japan [108], providing lateral resolution of 30 and 150 µm.
A comparison of NRA (the 1 H( 15 N, αγ) 12 C nuclear reaction) and ERDA was presented in [105]. NRA seems to outperform ERDA in terms of depth resolution and sensitivity (on the order of 10 weight ppm for NRA vs. 50-100 for ERDA). Furthermore, the processing of NRA data requires the use of less complex models (particle-matter interaction) than ERDA to extract distribution maps or H concentration profiles. However, in terms of lateral resolution, because of the size of the primary beam, ERDA outperforms the micro-beam NRA.
In addition to the 1 H isotope, a common NRA technique exists for deuterium, namely, the D( 3 He, p) 4 He reaction [109]. It produces protons of~13 MeV and α particles with 2-4 MeV energy that are easily identified. The D( 3 He, p) 4 He cross-section has a broad (~350 keV) peak at 650 keV. Since the stopping powers for the product α particles and protons are small, the reaction offers a very poor depth resolution (~500 nm near the surface). The large resonance width allows quantifying the D content in the entire nearsurface region [109].
NRA is a mature method and has been used in many application areas. It is often combined with ERDA to get more complete information. Hydrogen analysis by NRA in metals has involved, e.g., Al/Cu/Ag/Au multilayers [101], Mg/Ni films [102], tantalum [107], Zr/Cu/Al glassy alloys [108], and deuterium in W and WO 3 layers [61], in niobium [84,92], and beryllium [110]. Concerning semiconductors, NRA has been used, e.g., to analyze hydrogen in a-Si:H [37,111], silicon carbide [112], and silicon nitride [113]. Nonmetals are represented, e.g., by the analysis of deuterium implanted into carbon [114] and hydrogen analysis in diamond-like carbon coatings [115]. A number of other applications of NRA to hydrogen analysis were listed in [105].

Conclusions and Perspectives
Hydrogen is an element that affects many areas of science and technology, ranging from unwanted effects in classical metallurgy to electrochemical systems, corrosion of metals, the development of hydrogen-storage materials, semiconductors, various coatings, and thin-film systems. Hydrogen is also of interest in geochemistry, nuclear technologies, and many other areas. This yields an ever-increasing demand for the analysis of hydrogen in a very diverse range of materials, processes, and applications. The currently used bulk analysis and depth profiling methods for hydrogen in inorganic materials and inorganic coatings were reviewed in this paper. The bulk analysis methods of hydrogen are based on fusion or controlled heating of macroscopic samples and the detection of hydrogen, H 2 , in the released gases. They offer excellent accuracy and sensitivity. Depth profiling methods involve optical emission spectroscopy (GDOES and LIBS), mass spectrometry (GDMS and SIMS), and nuclear methods (ERDA and NRA). The principles of all these methods were explained in terms of their methodology, calibration procedures, analytical performance, and major application areas. The synergies and the complementarity of various methods of hydrogen analysis were described. The existing literature about all these methods was categorized and critically evaluated, while major papers concerning each method were selected and listed, in the context relevant to every aspect mentioned in the review.
Lastly, a few general remarks about hydrogen analysis should be made. Virtually all the methods listed here have been around for decades. Some have been developed to such a level that they can be considered as mature, with only little potential for further major improvements. This applies, e.g., to bulk analysis methods (Section 2). Concerning emission spectroscopies (Sections 3.1.1 and 3.1.2), a major issue is not the experiment itself but rather the interpretation of the data. In GDOES and GDMS, a number of processes involving hydrogen occur in the glow discharge plasma that are far from being completely described and understood. However, they significantly affect the measured intensities [43][44][45][46][47]. Likewise, in LIBS, analytical interpretation of the spectra is a very complex task, especially in the calibration-free approach. Further work on the fundamentals of these methods is desirable. In nuclear methods, significant instrumental developments have occurred relatively recently, such as hydrogen mapping by NRA.
Interesting is also the point of view of applications. Whereas the bulk analysis of hydrogen is considered rather routine and the relevant instrument manufacturers (e.g., LECO, Bruker, Eltra, Horiba, and some others) keep extensive libraries of application notes for a wide range of materials, for the methods listed in Section 3, the methodology for a particular analytical task is usually nontrivial and is, therefore, more frequently described in the scientific literature. In the papers focusing primarily on the coating under study, as well as its preparation, properties, etc., analytical aspects are sometimes treated incompletely, the results are often interpreted without considering the uncertainty of the analyses, etc. More attention to such aspects is usually paid in journals devoted to analytical science or instrumentation. A very fruitful approach to the analysis of coatings, particularly of hydrogen, is when the analytical methodology is developed hand in hand with the process itself, e.g., by preparing hydrogen-containing sample(s) that can be characterized by an independent method and subsequently used as a reference in more routine analyses, when optimizing the deposition technology itself.