Hydrogen Absorption in Metal Thin Films and Heterostructures Investigated in Situ with Neutron and X-ray Scattering

Due to hydrogen possessing a relatively large neutron scattering length, hydrogen absorption and desorption behaviors in metal thin films can straightforwardly be investigated by neutron reflectometry. However, to further elucidate the chemical structure of the hydrogen absorbing materials, complementary techniques such as high resolution X-ray reflectometry and diffraction remain important too. Examples of work on such systems include Nband Pd-based multilayers, where Nb and Pd both have strong affinity to hydrogen. W/Nb and Fe/Nb multilayers were measured in situ with unpolarized and polarized neutron reflectometry under hydrogen gas charging conditions. The gas-pressure/hydrogen-concentration dependence, the hydrogen-induced macroscopic film swelling as well as the increase in crystal lattice plane distances of the films were determined. Ferromagnetic-Co/Pd multilayers were studied with polarized neutron reflectometry and in situ ferromagnetic resonance measurements to understand the effect of hydrogen absorption on the magnetic properties of the system. This electronic effect enables a novel approach for hydrogen sensing using a magnetic readout scheme.


Introduction: Hydrogen Absorption in Metals-Bulk versus Thin Films
Thin films and layers of metals with chemical affinity to hydrogen such as Nb, Pd, Pt, V, Mg, Ti, the rare earth metals and related alloys play an important role in many technological and engineering materials and devices.A prime area of research is chemical energy storage in the vast field of metal hydrides [1][2][3].Here, the thin-film approach is of particular interest as it allows studying hydrogen-metal interactions on the nanoscale, which may show distinct novel features as well as non-equilibrium states of hydride materials that would not exist in bulk form [4][5][6].
Generally, the absorbed hydrogen occupies interstitial sites in the host crystal lattice.As a result, the host atoms are displaced from their regular lattice positions and the associated strain/stress fields will give rise to altered physical properties.Metal hydrides can show disordered and ordered arrangements of hydrogen sublattices, depending on the hydrogen concentration and temperature.
Metals 2016, 6, 125 3 of 20 emphasized that the adhesive forces at the metal-ceramic interface must be enormous as the in-plane stress resulting from hydrogen loading exceeded the yield stress of bulk Nb by at least one order of magnitude [25].It turns out that the yield stress for thin film/substrate systems may significantly differ from the corresponding bulk yield stress of that particular film material.For example, it was recently demonstrated that nanometer thick ultra-thin epitaxial Nb films can absorb very high amounts of hydrogen, c = 1 (H/Nb), reversibly without any noticeable plastic deformation or delamination from the substrate [26].The measured value of the hydrogen-induced stress was on the order of ´10 GPa, i.e., over an order of magnitude higher than the Nb bulk yield stress.The effective yield stress of a particular thin film/substrate system will depend on a number of factors, including the specific material combinations, layer thicknesses, crystalline quality and intrinsic stress levels of the films, temperature, etc.
Metals 2016, 6, 125 3 of 20 in-plane stress resulting from hydrogen loading exceeded the yield stress of bulk Nb by at least one order of magnitude [25].It turns out that the yield stress for thin film/substrate systems may significantly differ from the corresponding bulk yield stress of that particular film material.For example, it was recently demonstrated that nanometer thick ultra-thin epitaxial Nb films can absorb very high amounts of hydrogen, c = 1 (H/Nb), reversibly without any noticeable plastic deformation or delamination from the substrate [26].The measured value of the hydrogen-induced stress was on the order of −10 GPa, i.e., over an order of magnitude higher than the Nb bulk yield stress.The effective yield stress of a particular thin film/substrate system will depend on a number of factors, including the specific material combinations, layer thicknesses, crystalline quality and intrinsic stress levels of the films, temperature, etc. Charging bulk Nb with hydrogen, for example, usually results in an isotropic, three-dimensional lattice expansion of the Nb.In the α and α′ phases, the experimentally determined relation between relative expansion Δd/d0 of lattice parameters and the hydrogen concentration c is Δd/d0 = 0.058•c (H/Nb), and it is linear up to c (H/Nb) = 1 [21].For a hydrogen absorbing thin film being firmly clamped to a rigid substrate, one would find an anisotropic, one-dimensional out-of-plane lattice expansion.Due to the clamping to the substrate, no layer expansion within the layer plane is possible, and the proportionality factor will differ from the case of hydrogen charged bulk material.Theoretically, one would expect to find Δd/d0 = 0.136•c (H/Nb) [14].
The second important effect is the spatial proximity between a hydrogen absorbing layer and any neighboring material.Close to the interface, within the first few monolayers, drastic deviations from bulk behavior may occur due to charge-transfer processes altering the absorption potential for hydrogen [27].Thus, one can expect a film-thickness dependence of the average hydrogen solubility.Figure 2 shows this effect schematically.Charging bulk Nb with hydrogen, for example, usually results in an isotropic, three-dimensional lattice expansion of the Nb.In the α and α 1 phases, the experimentally determined relation between relative expansion ∆d/d 0 of lattice parameters and the hydrogen concentration c is ∆d/d 0 = 0.058¨c (H/Nb), and it is linear up to c (H/Nb) = 1 [21].For a hydrogen absorbing thin film being firmly clamped to a rigid substrate, one would find an anisotropic, one-dimensional out-of-plane lattice expansion.Due to the clamping to the substrate, no layer expansion within the layer plane is possible, and the proportionality factor will differ from the case of hydrogen charged bulk material.Theoretically, one would expect to find ∆d/d 0 = 0.136¨c (H/Nb) [14].
The second important effect is the spatial proximity between a hydrogen absorbing layer and any neighboring material.Close to the interface, within the first few monolayers, drastic deviations from bulk behavior may occur due to charge-transfer processes altering the absorption potential for hydrogen [27].Thus, one can expect a film-thickness dependence of the average hydrogen solubility.Figure 2 shows this effect schematically.It should be noted that the absorption of hydrogen by Nb lattice vacancies is energetically strongly preferred to absorption into interstitial sites.First-principle calculations of Nb hydride formation showed that a single vacancy can accommodate six hydrogen atoms in the symmetrically equivalent lowest-energy sites and additional hydrogen in the nearby interstitial sites affected by the strain field.This indicates that a vacancy can serve as a nucleation center for hydride phase formation [19].The trapping in vacancies was experimentally confirmed by positron annihilation experiments [28] where four hydrogen atoms were detected in one vacancy.The issue of hydrogen absorption and trapping in thin films (<100 nm) due to strain field effects is extensively discussed in the literature, see references 17-25 in [19], and in [29].
In view of the important role that the surface and surface-near layers play in the hydrogen absorption and diffusion process, it is imperative that the structural effects of hydrogen within these systems must be well understood.However, it is experimentally difficult to study hydrogen-induced structural changes in nanoscale layers.In this regard, neutrons provide a unique opportunity to directly probe hydrogen absorbed in metallic systems by scattering methods [8,[10][11][12][30][31][32].Often, X-rays are used to assess the structure of hydrogen absorbing materials and they are capable of observing effects such as lattice and layer expansion of the host crystal with increased hydrogen concentration [33][34][35][36][37][38].However, X-rays are nearly completely insensitive to the hydrogen atoms themselves since the X-ray scattering length is determined by the atomic number.In contrast, the scattering length for neutrons varies randomly from element to element, and, as a consequence, hydrogen in a given material usually scatters as strongly, or sometimes stronger, than the absorbing medium.The fact that both hydrogen and deuterium atoms strongly scatter neutrons opens up the opportunity for isotope contrast variation [39].These combined properties make neutron scattering a unique tool to identify the location and concentration of hydrogen in hydride materials.
Specifically for thin films, neutron reflectometry has proven to be an important technique when studying systems where the layered geometry is relevant to the behavior of the material or to their eventual applications.Neutron reflectometry is capable of determining the chemical profile of hydrogen and the host compound across thin films and multilayer samples up to roughly 200 nm in thickness.The depth sensitivity of this technique can elucidate properties such as hydrogen concentration gradients or the qualities of buried interfaces.Besides hydrogen sensitivity, neutrons have an additional advantage over X-rays in that they are strongly sensitive to magnetic moments, and, as such, polarized neutron reflectometry can be used to obtain magnetic depth profiles, which can be correlated with structural information.It should be noted that the absorption of hydrogen by Nb lattice vacancies is energetically strongly preferred to absorption into interstitial sites.First-principle calculations of Nb hydride formation showed that a single vacancy can accommodate six hydrogen atoms in the symmetrically equivalent lowest-energy sites and additional hydrogen in the nearby interstitial sites affected by the strain field.This indicates that a vacancy can serve as a nucleation center for hydride phase formation [19].The trapping in vacancies was experimentally confirmed by positron annihilation experiments [28] where four hydrogen atoms were detected in one vacancy.The issue of hydrogen absorption and trapping in thin films (<100 nm) due to strain field effects is extensively discussed in the literature, see references 17-25 in [19], and in [29].
In view of the important role that the surface and surface-near layers play in the hydrogen absorption and diffusion process, it is imperative that the structural effects of hydrogen within these systems must be well understood.However, it is experimentally difficult to study hydrogen-induced structural changes in nanoscale layers.In this regard, neutrons provide a unique opportunity to directly probe hydrogen absorbed in metallic systems by scattering methods [8,[10][11][12][30][31][32].Often, X-rays are used to assess the structure of hydrogen absorbing materials and they are capable of observing effects such as lattice and layer expansion of the host crystal with increased hydrogen concentration [33][34][35][36][37][38].However, X-rays are nearly completely insensitive to the hydrogen atoms themselves since the X-ray scattering length is determined by the atomic number.In contrast, the scattering length for neutrons varies randomly from element to element, and, as a consequence, hydrogen in a given material usually scatters as strongly, or sometimes stronger, than the absorbing medium.The fact that both hydrogen and deuterium atoms strongly scatter neutrons opens up the opportunity for isotope contrast variation [39].These combined properties make neutron scattering a unique tool to identify the location and concentration of hydrogen in hydride materials.
Specifically for thin films, neutron reflectometry has proven to be an important technique when studying systems where the layered geometry is relevant to the behavior of the material or to their eventual applications.Neutron reflectometry is capable of determining the chemical profile of hydrogen and the host compound across thin films and multilayer samples up to roughly 200 nm in thickness.The depth sensitivity of this technique can elucidate properties such as hydrogen concentration gradients or the qualities of buried interfaces.Besides hydrogen sensitivity, neutrons have an additional advantage over X-rays in that they are strongly sensitive to magnetic moments, and, as such, polarized neutron reflectometry can be used to obtain magnetic depth profiles, which can be correlated with structural information.This script focuses on utilizing both neutron and X-ray scattering techniques to study thin films and multilayers containing hydrogen absorbing metal components.In the first part of this review, we discuss hydrogen behavior in polycrystalline Fe/Nb multilayers.The results are compared with the case of epitaxially grown single-crystalline W/Nb [001] superlattices.The latter show distinctively different hydrogen absorption/desorption behavior due to a much stronger impact of the initial adherence to the substrate.The second portion of the work presented here reports on Pd/Co multilayers as a hybrid metallic-spintronic system for sensing hydrogen using a magnetic readout based on the ferromagnetic resonance technique.In this case, polarized neutron reflectometry proved especially useful as it allowed the observation of the magnetism of the Co layer simultaneous with hydrogen absorption studies.

Materials and Applications
2.1.1.Hydrogen in Fe/Nb and W/Nb Films Due to its high corrosion resistance, Nb is widely used across many types of applications, including hypoallergenic medical device and jewelry coatings, and linings in chemical plants.Despite the pure Nb surface being quite reactive, under most conditions, Nb will readily form a very thin oxide layer (Nb 2 O 5 being the most stable oxide), which is responsible for its chemical stability.However, a cold-worked piece of Nb metal manufactured for a real-world application will still possess numerous vacancies, dislocations and grain boundaries, and these can very effectively transport hydrogen from an external source, e.g., from processing materials, to the Nb bulk interior.Even relatively low hydrogen concentrations at levels of 100 ppm may substantially deteriorate particular properties, e.g., the ductility of Nb alloy construction materials [40].
For the work presented here, polycrystalline thin Nb layers of variable thickness were prepared by growing Fe/Nb multilayers on oxidized Si substrates [41].These samples were prepared by ion-beam sputtering at room temperature in an ultra-high vacuum chamber with a base pressure <5 ˆ10 ´9 mbar, where the Nb layers grew strongly (110) textured.Epitaxial Nb layers were prepared by growing W/Nb(001) superlattices at elevated temperatures of 200 ˝C on etched MgO(001) substrates by magnetron sputtering [42,43].Both Fe/Nb and W/Nb multilayers had a 5 nm thick Pd capping layer to prevent oxidation and to facilitate hydrogen absorption.Both systems were prepared in the form of multilayers to optimize the intensity for the scattering experiments discussed below.
The effect of hydrogen on Nb thin films was initially studied by us due to observed hydrogen-induced changes of the magnetic coupling and the related giant magnetoresistivity effect in Fe/Nb multilayers [10].The hydrogen-induced structural changes in the polycrystalline Fe/Nb system [41] were later compared with the behavior of epitaxially grown W/Nb superlattices [43].
For the Fe/Nb system, the Fe layer thickness, d Fe , was kept constant at 2.6 nm, while the Nb thickness, d Nb , varied from 2 to 100 nm.In the W/Nb system, the W layer thickness, d W , was 2.6 nm, and d Nb was 3 or 10 nm.The number of repeats of each Fe/Nb bilayer was chosen to provide a strong reflectivity signal allowing for the characterization of the Nb layers.Typical interface roughnesses, as determined by X-ray reflectivity, were 0.5 nm for Fe/Nb and 0.3 nm for W/Nb.A summary of the samples can be found in Table 1.Pd is one of the most utilized materials in hydrogen gas sensing due to its ability to reversibly absorb hydrogen [44].The absorption process is accompanied by several changes to the physical properties of Pd, including the aforementioned lattice expansion of up to several percent and an increase in resistivity as high as 80% for PdH 0.7 at 25 ˝C [45].Pd is also quite selective regarding hydrogen absorption and exhibits low sensitivity to other gases, such as CO, Cl 2 , SO 2 , H 2 S, NO x , and hydrocarbons, making it a suitable material in hydrogen sensors [46].The challenge to reach the application stage lies in how to most effectively sense these hydrogen-induced changes to Pd.Current devices have several drawbacks: poor sensitivity, complex detection systems, slow response time, high power consumption, or potential flammability [18].
A novel approach to this problem lies in exploiting one of the lesser known aspects of Pd: magnetic properties.Although Pd itself is paramagnetic, when it neighbors ferromagnetic metals, such as Co, Fe, and NiFe, a perpendicular magnetic anisotropy (PMA) develops at the interface [47].For ultra-thin ferromagnetic films (~0.5-1 nm), PMA aligns the magnetization of the ferromagnetic layer in the out-of-plane direction.For increased film thickness of the ferromagnet, the magnetization direction lies in the film plane, but the presence of PMA at the interface can affect a number of magnetic properties, including ferromagnetic resonance (FMR).Monitoring the magnetic response of these systems may then be a way to approach the design of novel hydrogen gas sensors.
Here, work is presented on a (10 nm) Co/(10 nm) Pd bilayer deposited via sputtering on a Si substrate [15], with the Pd layer on top allowing more Pd surface area to be exposed to external hydrogen gas.

Neutron Reflectometry
Reflectometry uses small-angle specular reflection (θ i = θ f ) to probe the chemical depth profile across a surface, thin film, or multilayer as a function of depth.The geometry of this technique is shown in Figure 3 where the scattering vector, Q, is given by Q = k f ´ki .Q is related to the wavelength of the incoming neutrons, λ, and the angle of reflection, θ, by Metals 2016, 6, 125 6 of 20 Pd is one of the most utilized materials in hydrogen gas sensing due to its ability to reversibly absorb hydrogen [44].The absorption process is accompanied by several changes to the physical properties of Pd, including the aforementioned lattice expansion of up to several percent and an increase in resistivity as high as 80% for PdH0.7 at 25 °C [45].Pd is also quite selective regarding hydrogen absorption and exhibits low sensitivity to other gases, such as CO, Cl2, SO2, H2S, NOx, and hydrocarbons, making it a suitable material in hydrogen sensors [46].The challenge to reach the application stage lies in how to most effectively sense these hydrogen-induced changes to Pd.Current devices have several drawbacks: poor sensitivity, complex detection systems, slow response time, high power consumption, or potential flammability [18].
A novel approach to this problem lies in exploiting one of the lesser known aspects of Pd: magnetic properties.Although Pd itself is paramagnetic, when it neighbors ferromagnetic metals, such as Co, Fe, and NiFe, a perpendicular magnetic anisotropy (PMA) develops at the interface [47].For ultra-thin ferromagnetic films (~0.5-1 nm), PMA aligns the magnetization of the ferromagnetic layer in the out-of-plane direction.For increased film thickness of the ferromagnet, the magnetization direction lies in the film plane, but the presence of PMA at the interface can affect a number of magnetic properties, including ferromagnetic resonance (FMR).Monitoring the magnetic response of these systems may then be a way to approach the design of novel hydrogen gas sensors.
Here, work is presented on a (10 nm) Co/(10 nm) Pd bilayer deposited via sputtering on a Si substrate [15], with the Pd layer on top allowing more Pd surface area to be exposed to external hydrogen gas.

Neutron Reflectometry
Reflectometry uses small-angle specular reflection (θi = θf) to probe the chemical depth profile across a surface, thin film, or multilayer as a function of depth.The geometry of this technique is shown in Figure 3 where the scattering vector, Q, is given by Q = kf − ki.Q is related to the wavelength of the incoming neutrons, λ, and the angle of reflection, θ, by  At each surface and interface, some fraction of the incoming beam is reflected and the other fraction is transmitted and refracted.As neutrons pass through each layer, they are exposed to an effective nuclear scattering potential given as where b is the scattering length of the element, N is its atomic density and m is the neutron mass.The term bN is the scattering length density (SLD) and is related to the chemical makeup of a material as well as its physical density of the contained elements (for compounds: b total ¨Ntotal = b 1 ¨N1 + b 2 ¨N2 + . . .+ b i ¨Ni ) [48].Reflectometry allows one to model the SLD as a function of sample depth, which can be used to link features like hydrogen concentration to the change in SLD from a material's expected bulk value.
As previously mentioned, neutrons as scattering probes are particularly well suited to investigate hydrogen-sensitive materials because the magnitude of the hydrogen scattering length is of the same order as for many other metals.However, hydrogen has a negative scattering length, which can lead to its easy detection by comparing reflections from hydrogenated and non-hydrogenated materials.An example of this is demonstrated in Figure 4, where neutron reflectivity curves are simulated for a bulk Pd surface (red solid line), which would have an SLD of 402 µm ´2, as compared to the α 1 phase of PdH (blue dashed line), where the SLD changes to 126 µm ´2.At each surface and interface, some fraction of the incoming beam is reflected and the other fraction is transmitted and refracted.As neutrons pass through each layer, they are exposed to an effective nuclear scattering potential given as where b is the scattering length of the element, N is its atomic density and m is the neutron mass.The term bN is the scattering length density (SLD) and is related to the chemical makeup of a material as well as its physical density of the contained elements (for compounds: btotal [48].Reflectometry allows one to model the SLD as a function of sample depth, which can be used to link features like hydrogen concentration to the change in SLD from a material's expected bulk value. As previously mentioned, neutrons as scattering probes are particularly well suited to investigate hydrogen-sensitive materials because the magnitude of the hydrogen scattering length is of the same order as for many other metals.However, hydrogen has a negative scattering length, which can lead to its easy detection by comparing reflections from hydrogenated and non-hydrogenated materials.An example of this is demonstrated in Figure 4, where neutron reflectivity curves are simulated for a bulk Pd surface (red solid line), which would have an SLD of 402 μm −2 , as compared to the α′ phase of PdH (blue dashed line), where the SLD changes to 126 μm −2 .Neutron reflectivity measurements on the Fe/Nb and W/Nb samples were performed at the V6 reflectometer [49,50] at the Hahn-Meitner-Institut Berlin, Germany, and the POSY2 reflectometer [51] at the Intense Pulsed Neutron Source at Argonne National Laboratory, Argonne, IL, USA.

Polarized Neutron Reflectometry (PNR)
Polarized neutron reflectometry (PNR) exploits the spin of the neutron as an additional tool for observing magnetic depth profiles across samples.In addition to the nuclear scattering potential in Equation (2), neutrons travelling through a magnetic layer will be exposed to a magnetic scattering potential given by where B is the magnetic induction vector, s is the neutron spin, and bmag is the magnetic scattering length (1 μB/atom results in bmag = 2.695 fm).Depending on the direction of the incoming neutron Neutron reflectivity measurements on the Fe/Nb and W/Nb samples were performed at the V6 reflectometer [49,50] at the Hahn-Meitner-Institut Berlin, Germany, and the POSY2 reflectometer [51] at the Intense Pulsed Neutron Source at Argonne National Laboratory, Argonne, IL, USA.

Polarized Neutron Reflectometry (PNR)
Polarized neutron reflectometry (PNR) exploits the spin of the neutron as an additional tool for observing magnetic depth profiles across samples.In addition to the nuclear scattering potential in Equation (2), neutrons travelling through a magnetic layer will be exposed to a magnetic scattering potential given by where B is the magnetic induction vector, ŝ is the neutron spin, and b mag is the magnetic scattering length (1 µ B /atom results in b mag = 2.695 fm).Depending on the direction of the incoming neutron spin, the total scattering potential (combining Equations ( 2) and ( 3)) is then V tot = V n ˘Vm .By performing reflectivity measurements with two incoming, antiparallel neutron polarization states (+ and ´), the difference between the R+ and R´channels essentially excludes the nuclear scattering and can be used to ascertain a profile of the magnetization across thin layers.Ferromagnetic Resonance (FMR) is the effect of spin precession of the macroscopic magnetization vector in an external magnetic field and is driven experimentally by microwave power [52].The FMR signal for a magnetic sample is observed as a sharp increase in the absorption of microwave power by the sample in a magnetic resonance condition.FMR has been applied to Pd/Co bilayers and multilayers in order to study how magnetic properties associated with perpendicular magnetic anisotropy (PMA) change with the hydrogen concentration in the Pd layer in order to investigate the feasibility of this system to act as a metallic spintronic hydrogen sensor [15].In Ref. [15], it was reported that hydrogen absorption resulted in narrowing and shifting of the FMR absorption peak for the Co layer, indicating a coupling between these two properties.
In situ FMR and PNR measurements were carried out as a function of external hydrogen gas pressure at the PLATYPUS beamline [53][54][55] at the OPAL reactor at ANSTO, Lucas Heights, Australia.Figure 5 schematically shows the custom-made scattering chamber with FMR capability.The samples were placed face-up across the microwave coplanar transmission line, and aluminum windows mounted at the front and the rear of the chamber allowed the neutron beam to enter and exit.The gas atmosphere within the chamber was controlled by a mass flow controller system capable of delivering pure N 2 gas or N 2 /H 2 gas mixtures of ď3.5% hydrogen gas at slightly above normal atmospheric pressure and at an ambient temperature of 23 ˝C.The external magnetic field was generated using a 1 T electromagnet.spin, the total scattering potential (combining Equations ( 2) and ( 3)) is then Vtot = Vn ± Vm.By performing reflectivity measurements with two incoming, antiparallel neutron polarization states (+ and −), the difference between the R+ and R− channels essentially excludes the nuclear scattering and can be used to ascertain a profile of the magnetization across thin layers.Ferromagnetic Resonance (FMR) is the effect of spin precession of the macroscopic magnetization vector in an external magnetic field and is driven experimentally by microwave power [52].The FMR signal for a magnetic sample is observed as a sharp increase in the absorption of microwave power by the sample in a magnetic resonance condition.FMR has been applied to Pd/Co bilayers and multilayers in order to study how magnetic properties associated with perpendicular magnetic anisotropy (PMA) change with the hydrogen concentration in the Pd layer in order to investigate the feasibility of this system to act as a metallic spintronic hydrogen sensor [15].In Ref. [15], it was reported that hydrogen absorption resulted in narrowing and shifting of the FMR absorption peak for the Co layer, indicating a coupling between these two properties.
In situ FMR and PNR measurements were carried out as a function of external hydrogen gas pressure at the PLATYPUS beamline [53][54][55] at the OPAL reactor at ANSTO, Lucas Heights, Australia.Figure 5 schematically shows the custom-made scattering chamber with FMR capability.The samples were placed face-up across the microwave coplanar transmission line, and aluminum windows mounted at the front and the rear of the chamber allowed the neutron beam to enter and exit.The gas atmosphere within the chamber was controlled by a mass flow controller system capable of delivering pure N2 gas or N2/H2 gas mixtures of 3.5% hydrogen gas at slightly above normal atmospheric pressure and at an ambient temperature of 23 °C.The external magnetic field was generated using a 1 T electromagnet.

Results and Discussion
Nb Thin Films and Multilayers

Fe/Nb Multilayers
In order to study the correlation between hydrogen absorption and resulting changes of the structural properties of Nb films, Fe/Nb heterostructures were measured with neutron reflectometry  (NR) and X-ray diffraction/reflectometry [41][42][43].The multilayer geometry is advantageous as the scattering signals are strongly amplified.Using Fe as a partnering element is an apt choice, as it does not absorb hydrogen but is sufficiently transparent for hydrogen diffusion, which was experimentally confirmed.Multilayer samples were charged in situ with hydrogen at various pressures from 10 ´4 to 900 mbar.Measurements were made at a sample temperature of 185 ˝C to allow for relatively fast hydrogen absorption and to avoid phase separation that can occur at lower temperatures.For the Fe/Nb system (and for W/Nb discussed below), we present experimental data collected in the first cycle of hydrogen absorption.
Sample NR curves for polycrystalline samples are shown in Figure 6, comparing the systems in the virgin state and at the largest hydrogen pressure.The strong peaks in the top three sections are Bragg peaks resulting from the superlattice geometry.There is a noticeable shift in the peak position to a lower Q value with increased hydrogen content, indicating an expansion of the Nb layers within the heterostructures.In addition, we observe a significant increase in the Bragg peak intensity with increasing hydrogen absorption due to the increase in scattering contrast between Fe and Nb + H layers. (NR) and X-ray diffraction/reflectometry [41][42][43].The multilayer geometry is advantageous as the scattering signals are strongly amplified.Using Fe as a partnering element is an apt choice, as it does not absorb hydrogen but is sufficiently transparent for hydrogen diffusion, which was experimentally confirmed.Multilayer samples were charged in situ with hydrogen at various pressures from 10 −4 to 900 mbar.Measurements were made at a sample temperature of 185 °C to allow for relatively fast hydrogen absorption and to avoid phase separation that can occur at lower temperatures.For the Fe/Nb system (and for W/Nb discussed below), we present experimental data collected in the first cycle of hydrogen absorption.Sample NR curves for polycrystalline samples are shown in Figure 6, comparing the systems in the virgin state and at the largest hydrogen pressure.The strong peaks in the top three sections are Bragg peaks resulting from the superlattice geometry.There is a noticeable shift in the peak position to a lower Q value with increased hydrogen content, indicating an expansion of the Nb layers within the heterostructures.In addition, we observe a significant increase in the Bragg peak intensity with increasing hydrogen absorption due to the increase in scattering contrast between Fe and Nb + H layers. PNR was also performed due to the ferromagnetism of the Fe layers.While the Fe magnetism was found to be almost entirely unaffected by the Nb hydrogen concentration, the ability to measure two data channels for each pressure allowed for more certainty in the determination of the SLD of each material.
The data were fitted to ascertain the SLDs and thicknesses of the Nb and Fe layers.From the SLDs, the hydrogen concentrations can be calculated from [41]: PNR was also performed due to the ferromagnetism of the Fe layers.While the Fe magnetism was found to be almost entirely unaffected by the Nb hydrogen concentration, the ability to measure two data channels for each pressure allowed for more certainty in the determination of the SLD of each material.
The data were fitted to ascertain the SLDs and thicknesses of the Nb and Fe layers.From the SLDs, the hydrogen concentrations can be calculated from [41]: Metals where the subscripts represent Nb before and after the introduction of hydrogen, d is the thickness of the Nb layer, and c H is the number of hydrogen atoms per Nb metal atom.Data sets could be fitted well assuming a homogenous distribution of hydrogen throughout the Nb thickness.A sample of the fitted PNR data is shown in Figure 7.The resulting SLDs and layer thicknesses are shown in Table 2.
where the subscripts represent Nb before and after the introduction of hydrogen, d is the thickness of the Nb layer, and cH is the number of hydrogen atoms per Nb metal atom.Data sets could be fitted well assuming a homogenous distribution of hydrogen throughout the Nb thickness.A sample of the fitted PNR data is shown in Figure 7.The resulting SLDs and layer thicknesses are shown in Table 2.The most striking differences between the high and low hydrogen pressure data occur in the Nb layer.First, there is a decrease in SLD with hydrogen, from 380 to 170 μm −2 , which is accounted for by the negative hydrogen scattering length.There is also an increase in the Nb thickness by nearly 10 nm or 10%.This data can be used with Equation (4) to calculate the hydrogen concentration, and the example data indicated a hydrogen concentration of 95% (H/Nb) at 900 mbar external hydrogen pressure, i.e., nearly one hydrogen atom for each Nb atom.This analysis was performed for several samples over a range of hydrogen pressures, leading to the development of a Nb/hydrogen concentration/pressure phase diagram shown in Figure 8 [56].Here, it is demonstrated that in  The most striking differences between the high and low hydrogen pressure data occur in the Nb layer.First, there is a decrease in SLD with hydrogen, from 380 to 170 µm ´2, which is accounted for by the negative hydrogen scattering length.There is also an increase in the Nb thickness by nearly 10 nm or 10%.This data can be used with Equation (4) to calculate the hydrogen concentration, and the example data indicated a hydrogen concentration of 95% (H/Nb) at 900 mbar external hydrogen pressure, i.e., nearly one hydrogen atom for each Nb atom.This analysis was performed for several samples over a range of hydrogen pressures, leading to the development of a Nb/hydrogen concentration/pressure phase diagram shown in Figure 8 [56].Here, it is demonstrated that in thin-film form, the phase diagram is drastically different than in bulk, particularly for samples with Nb thicknesses below 5 nm, where the peak hydrogen concentration reached is much less than the bulk value, and the pressures for obtaining a particular hydrogen concentration are orders of magnitude higher.As the W/Nb system shows similar behavior, we will discuss the causes for thickness dependence of the hydrogen absorption of both systems together in the Section 3.1.2.
thin-film form, the phase diagram is drastically different than in bulk, particularly for samples with Nb thicknesses below 5 nm, where the peak hydrogen concentration reached is much less than the bulk value, and the pressures for obtaining a particular hydrogen concentration are orders of magnitude higher.As the W/Nb system shows similar behavior (see Figure 10), we will discuss the causes for thickness dependence of the hydrogen absorption of both systems together in the final Section of 3.1.2.

W/Nb Multilayers
This section describes hydrogen absorption of Nb layers embedded in single-crystal (epitaxial) W/Nb multilayers and demonstrates how the degree of structural perfection impacts the hydrogen solubility in comparison to the previously presented polycrystalline Fe/Nb system.Similar to Fe, W does not accumulate hydrogen in any significant way due to its positive solution enthalpy for hydrogen.However, at T = 185 °C, W as well as Fe is sufficiently transparent to allow a homogeneous charging of all the Nb layers in the layer stack.
In contrast to the Fe/Nb system, the W/Nb samples were charged with deuterium gas.The reason is as follows: In neutron reflectivity measurements of the Fe/Nb system, the scattering contrast between Fe (high SLD) and Nb + H (lower SLD) continuously increases with increasing hydrogen absorption due to the negative scattering length of hydrogen (the averaged bH of hydrogen gas with naturally occurring isotope fractions is −3.739 fm).In contrast, hydrogen absorption in Nb layers of W/Nb would decrease the scattering contrast between W and Nb which would have made analysis of the data more difficult.Using deuterium with its positive scattering length of bD = 6.671 fm avoids this problem.It has previously been shown on bulk [57,58] and on thin layers [59] that the solubility of deuterium in Nb deviates insignificantly from that of hydrogen.Note that, in order to facilitate the readability of the text, we continue to use the term hydrogen.
Results of in situ neutron reflectivity measurements on W/Nb multilayers are shown in Figure 9.The measurement on the uncharged sample is represented by red circles while the blue triangles represent the spectrum measured after hydrogen charging at   [56] with permission from Elsevier.

W/Nb Multilayers
This section describes hydrogen absorption of Nb layers embedded in single-crystal (epitaxial) W/Nb multilayers and demonstrates how the degree of structural perfection impacts the hydrogen solubility in comparison to the previously presented polycrystalline Fe/Nb system.Similar to Fe, W does not accumulate hydrogen in any significant way due to its positive solution enthalpy for hydrogen.However, at T = 185 ˝C, W as well as Fe is sufficiently transparent to allow a homogeneous charging of all the Nb layers in the layer stack.
In contrast to the Fe/Nb system, the W/Nb samples were charged with deuterium gas.The reason is as follows: In neutron reflectivity measurements of the Fe/Nb system, the scattering contrast between Fe (high SLD) and Nb + H (lower SLD) continuously increases with increasing hydrogen absorption due to the negative scattering length of hydrogen (the averaged b H of hydrogen gas with naturally occurring isotope fractions is ´3.739 fm).In contrast, hydrogen absorption in Nb layers of W/Nb would decrease the scattering contrast between W and Nb which would have made analysis of the data more difficult.Using deuterium with its positive scattering length of b D = 6.671 fm avoids this problem.It has previously been shown on bulk [57,58] and on thin layers [59] that the solubility of deuterium in Nb deviates insignificantly from that of hydrogen.Note that, in order to facilitate the readability of the text, we continue to use the term hydrogen.
Results of in situ neutron reflectivity measurements on W/Nb multilayers are shown in Figure 9.The measurement on the uncharged sample is represented by red circles while the blue triangles represent the spectrum measured after hydrogen charging at P H 2 = 900 mbar.The influence of the hydrogen is clearly reflected in both the shift of the Bragg reflections which originate from the W/Nb bilayer period and their apparent increase in intensity due to the change in scattering contrast.Due to the lower double-layer thickness of the sample with 3 nm Nb, however, only one Bragg reflection can be found in the experimentally accessible Q region.
The neutron reflectivity measurement of the sample with [2.6 nm W/3.0 nm Nb]60 features a specific characteristic: in the Q region of about 1.1 nm −1 , the structural Bragg reflection is significantly asymmetrical.Such an effect can be observed if the sample possesses a slight layer thickness gradient in the lateral direction, which is present here and needed to be taken into account in the data analysis process.
Figure 10 shows the experimentally determined hydrogen solubility curves of the two W/Nb samples with 3 nm and 10 nm Nb, respectively, together with the theoretical curve of Nb bulk material.
For the sample with 3 nm Nb, a maximal hydrogen concentration of about c (H/Nb) = 40 at.% can be found, whereas this value amounts to about 80 at.% for the sample with 10 nm Nb (both at 2 H P = 900 mbar).The latter approaches the solubility of bulk samples at high pressures.As for the Fe/Nb multilayer system, the plateau-like regions-where most of the hydrogen is absorbed-depend on the Nb layer thickness.For 10 nm Nb, such a region is one order of magnitude, and, for 3 nm Nb, such a region is two orders of magnitude higher in pressure than for the Nb bulk material.The concentration/pressure gradient within this region increases with thinner Nb layer thickness.This indicates that the critical temperature for phase separation is significantly reduced for thin Nb layers compared to thicker Nb layers.
In order to elucidate the impact of hydrogen loading on the Nb crystal structure, in situ X-ray diffraction measurements were performed on the samples [2.6 nm W/10 nm Nb]60 and [2.6 nm W/3 nm Nb]120, see Figure 11, during the hydrogen charging at increasing hydrogen pressures.The data were taken with the scattering vector directing out-of-plane.The influence of the hydrogen is clearly reflected in both the shift of the Bragg reflections which originate from the W/Nb bilayer period and their apparent increase in intensity due to the change in scattering contrast.Due to the lower double-layer thickness of the sample with 3 nm Nb, however, only one Bragg reflection can be found in the experimentally accessible Q region.
The neutron reflectivity measurement of the sample with [2.6 nm W/3.0 nm Nb] 60 features a specific characteristic: in the Q region of about 1.1 nm ´1, the structural Bragg reflection is significantly asymmetrical.Such an effect can be observed if the sample possesses a slight layer thickness gradient in the lateral direction, which is present here and needed to be taken into account in the data analysis process.
Figure 10 shows the experimentally determined hydrogen solubility curves of the two W/Nb samples with 3 nm and 10 nm Nb, respectively, together with the theoretical curve of Nb bulk material.
For the sample with 3 nm Nb, a maximal hydrogen concentration of about c (H/Nb) = 40 at.% can be found, whereas this value amounts to about 80 at.% for the sample with 10 nm Nb (both at P H 2 = 900 mbar).The latter approaches the solubility of bulk samples at high pressures.As for the Fe/Nb multilayer system, the plateau-like regions-where most of the hydrogen is absorbed-depend on the Nb layer thickness.For 10 nm Nb, such a region is one order of magnitude, and, for 3 nm Nb, such a region is two orders of magnitude higher in pressure than for the Nb bulk material.The concentration/pressure gradient within this region increases with thinner Nb layer thickness.This indicates that the critical temperature for phase separation is significantly reduced for thin Nb layers compared to thicker Nb layers.
In order to elucidate the impact of hydrogen loading on the Nb crystal structure, in situ X-ray diffraction measurements were performed on the samples [2.6 nm W/10 nm Nb] 60 and [2.6 nm W/3 nm Nb] 120 , see Figure 11, during the hydrogen charging at increasing hydrogen pressures.The data were taken with the scattering vector directing out-of-plane.A reflection characteristic of the thin-film Nb crystal structure is located at a Bragg angle of Θ « 27.3 ˝for both samples, a value that is shifted to smaller reflection angles when compared to bulk samples indicating larger Nb/Nb lattice plane distances in the out-of-plane direction.This is most likely due to a close matching of the W and Nb in-plane lattice parameters due to the epitaxy.We call the peak at this scattering angle Phase I reflection.Charging the sample with hydrogen changes the spectrum as expected.Most significantly the intensity of the Phase I reflection is highly reduced at higher pressures, its full width at half maximum increased, and its position shifted to smaller Θ values.At higher gas pressures, a new reflection can be identified at Θ « 25.9 ˝.We assume that a transition occurs between these two phases: first a slight one-dimensional lattice expansion occurs (the Phase I peak shifts towards lower angles), followed by a transition of Nb volume into a vertically strongly expanded phase with significantly higher hydrogen content (Phase II).Note that Phase II grows at the expense of Phase I (see bottom panels of Figure 11) which shows the time dependence of the Nb crystal transformation upon hydrogenation.Phase II is most pronounced at P H 2 = 900 mbar for both samples.
Our interpretation of this transformation is that above a Nb thickness dependent critical hydrogen pressure, the Nb lattice cannot keep its structural coherence to the W and the substrate lattice.It seems that the initially only slightly out-of-plane expanded Nb (marked as Phase I) undergoes a non-coherent transition to a largely out-of-plane expanded Phase II at hydrogen pressures above 1 mbar (see below for details).We note that only a very small fraction of the Nb is still in Phase I at 900 mbar, but the sum of the integrated intensities of the Phases I and II peaks stays constant for all hydrogen pressures.
Figure 12 compares relative changes of the Nb interplanar spacing and the Nb layer thickness for the two in situ hydrogen-charged W/Nb samples.As evidenced by the in situ X-ray diffraction measurements, the data sets reveal a second, unexpected hydrogen absorption mechanism which seems to be of equal importance to the hydrogen absorption on interstitial sites.A reflection characteristic of the thin-film Nb crystal structure is located at a Bragg angle of Θ ≈ 27.3° for both samples, a value that is shifted to smaller reflection angles when compared to bulk samples indicating larger Nb/Nb lattice plane distances in the out-of-plane direction.This is most likely due to a close matching of the W and Nb in-plane lattice parameters due to the epitaxy.We call the peak at this scattering angle Phase I reflection.Charging the sample with hydrogen changes the spectrum as expected.Most significantly the intensity of the Phase I reflection is highly reduced at higher pressures, its full width at half maximum increased, and its position shifted to smaller Θ values.At higher gas pressures, a new reflection can be identified at Θ ≈ 25.9°.We assume that a transition occurs between these two phases: first a slight one-dimensional lattice expansion occurs (the Phase I peak shifts towards lower angles), followed by a transition of Nb volume into a vertically strongly expanded phase with significantly higher hydrogen content (Phase II).Note that Phase II grows at the expense of Phase I (see bottom panels of Figure 11) which shows the time dependence of the Nb crystal transformation upon hydrogenation.Phase II is most pronounced at 2 H P = 900 mbar for both samples.
Our interpretation of this transformation is that above a Nb thickness dependent critical hydrogen pressure, the Nb lattice cannot keep its structural coherence to the W and the substrate lattice.It seems that the initially only slightly out-of-plane expanded Nb (marked as Phase I) undergoes a non-coherent transition to a largely out-of-plane expanded Phase II at hydrogen pressures above 1 mbar (see below for details).We note that only a very small fraction of the Nb is still in Phase I at 900 mbar, but the sum of the integrated intensities of the Phases I and II peaks stays constant for all hydrogen pressures.
Figure 12 compares relative changes of the Nb interplanar spacing and the Nb layer thickness for the two in situ hydrogen-charged W/Nb samples.As evidenced by the in situ X-ray diffraction measurements, the data sets reveal a second, unexpected hydrogen absorption mechanism which seems to be of equal importance to the hydrogen absorption on interstitial sites.The data show that the relative increase in the Nb layer thickness of both samples is significantly larger than the expansion of the Nb interplanar spacing.This result is consistent with what was found in the Fe/Nb system except that during the hydrogen absorption at no time a distinct transformation from a low out-of-plane state to a high out-of-plane state could be observed in the Fe/Nb system.Overall, the effects of changes in Nb layer thickness and Nb interplanar spacing through hydrogen charging during a first charging cycle in Fe/Nb and W/Nb (Phase II) are comparable.Also numerically, the same values can approximately be found.In both systems, the layer expansion in the direction of growth is considerably larger than the corresponding lattice expansion.This clearly demonstrates that, besides the well-known process of absorption at interstitial sites that entails expansion of the atomic lattice planes, there is an additional hydrogen absorption mechanism at work, which leads to a large macroscopic swelling of the Nb film.
An explanation for the large macroscopic swelling effect and for the distinct transformation between Phase I and Phase II in W/Nb was found by performing extended X-ray absorption fine structure (EXAFS) studies on the same series of samples [43].It turned out that Phase II is related to a state in which, due to the large hydrogen concentrations, the in-plane Nb lattice parameter is significantly expanded indicating loss of coherency with the W lattice.In contrast, the Nb in-plane lattice parameter related to the low hydrogen concentration Phase I is smaller, still epitaxially matching the W in-plane lattice parameter.For the (001) orientated epitaxial [2.6 nm W/10 nm Nb] 60 multilayer, EXAFS revealed a small 2.5% in-plane expansion of the Nb lattice under hydrogen absorption.The in-plane expansion of 2.5% and the corresponding out-of-plane lattice expansion of 6.2% result in a total increase in volume of 11.5%.This amount coincides with the relative increase of the Nb layer thickness expansion of 11.1% as determined from reflectivity experiments.Hence, one only has to take into account the boundary condition of a fixed in-plane area (i.e., a one-dimensional degree of freedom to expand macroscopically in the out-of-plane direction) in order to explain quantitatively the results in Figure 12.This model implies a three-dimensional rearrangement of Nb atoms caused by the massive mechanical lateral strain and its relaxation with increasing hydrogen charging within the Nb layers, especially near the interfaces to the substrate and the unloaded Fe or W layers.As a result, individual Nb atoms are squeezed out of existing lattice planes and start to form additional planes, which finally causes the anomalous large expansion of the layer thickness.Note that the newly created partial lattice planes imply a large amount of additional dislocations.The latter are evidenced by the fact that the out-of-plane X-ray diffraction peaks broaden significantly upon hydrogen absorption.For [2.6 nm W/10 nm Nb] 60 , for example, we find a decrease in coherence length in growth direction from 9.0 nm after preparation, to 5.5 nm after hydrogen absorption at 900 mbar.Since lattice imperfections effectively trap hydrogen atoms, the new effect likely also explains the extraordinarily high hydrogen concentrations found in our Nb films and the partial non-reversibility of the hydrogen absorption process itself.

Co/Pd Bilayers
A combination of in situ PNR and FMR measurements was performed on Co/Pd bilayer films in order to study the magnetic effects of hydrogen absorption in this ferromagnet/metallic system and to correlate FMR changes with structural and/or magnetic changes in each individual layer.
Figure 13 shows a comparison of the experimental R+ reflectivity data and the theoretical model fits of a Co/Pd bilayer before (solid red line) and after (solid blue line) absorption of hydrogen from an ambient pressure N 2 /H 2 gas mixture containing 3.5% hydrogen partial pressure (this corresponds to an absolute hydrogen partial pressure of «35.5 mbar within the loading chamber).The experimental data set prior to hydrogen absorption was fitted using a model consisting of d Pd = 13.50 nm and d Co = 8.36 nm.However, a best fit to the data was obtained by considering an interfacial region of thickness d interface = 1.43 nm and SLD interface = 162 µm ´2 in the virgin state.This interface layer had a rather large roughness of 1.0 nm at its transition to the Pd layer.After exposure to a 3.5% hydrogen gas partial pressure the interfacial layer was found to undergo a 25% reduction in SLD accompanied by an increase in the in-plane magnetic moment by 66%, and a sharpening of its interfacial roughness to the Pd layer from 1.0 nm to 0.2 nm.Due to the negative scattering length of hydrogen, these observed changes to structure and magnetism highlight the presence of hydrogen at the Co/Pd interface, which have additionally resulted in modifications to the PMA strength within the Pd lattice.Determining the Pd SLDs and exact layer thicknesses before and after hydrogen exposure and using Equation ( 4) established the number density ratio of hydrogen atoms to Pd atoms to be 30% in the 3.5% hydrogen partial pressure state.Similar to the Fe/Nb (see Figure 2 in Ref. [41]) and W/Nb (see Figure 5 in Ref. [43]) systems, we found that, after the initial hydrogen absorption and desorption cycle was performed on the Co/Pd sample, the Pd layer did not return to its original thickness or density, demonstrating that the presence of hydrogen within the Pd led to the creation of irreversible deformations.It is significant, however, to note that a reproducible and reversible Pd structure was established between the second and third hydrogen absorption and desorption cycling experiments.This "training" effect needs to be taken into account for technical applications.Reversibility and training effects of the PdHx system have been extensively studied in References [60,61].
The in situ FMR signal responses of the 8.36 nm Co/13.50 nm Pd bilayer obtained on PLATYPUS during two hydrogen absorption and desorption cycles are shown in Figure 14.During these measurements, the microwave frequency was kept constant and equal to 3.4 GHz, and the applied field was swept across the resonance to produce the plots shown in the figure.Ensuring efficient neutron reflection from the film requires taking FMR measurements in the "film facing up" geometry (see Section 2.2.3).For this film orientation, the FMR signal is significantly reduced because of film separation from the microwave transmission line by the 0.5-mm-thick silicon substrate of the film.In order to compensate for this signal reduction, the applied-field modulation technique combined with lock-in signal detection was utilized [62].For this reason, the registered FMR traces have the shape of the first derivative of a Lorentzian instead of the Lorentzian itself.
The characteristic resonance changes its location with respect to the applied magnetic field depending on the hydrogen concentration within the Pd layer, shifting to lower applied fields upon hydrogen absorption, and shifting upwards upon hydrogen desorption.The downshift of the resonance peak implies that the PMA strength is reduced when hydrogen atoms are incorporated into the Pd crystal lattice.Reproducible responses analogous to those observed in the PNR data between the second and third hydrogen absorption and desorption cycling stages were also evident in the FMR data.Determining the Pd SLDs and exact layer thicknesses before and after hydrogen exposure and using Equation ( 4) established the number density ratio of hydrogen atoms to Pd atoms to be 30% in the 3.5% hydrogen partial pressure state.Similar to the Fe/Nb (see Figure 2 in Ref. [41]) and W/Nb (see Figure 5 in Ref. [43]) systems, we found that, after the initial hydrogen absorption and desorption cycle was performed on the Co/Pd sample, the Pd layer did not return to its original thickness or density, demonstrating that the presence of hydrogen within the Pd led to the creation of irreversible deformations.It is significant, however, to note that a reproducible and reversible Pd structure was established between the second and third hydrogen absorption and desorption cycling experiments.This "training" effect needs to be taken into account for technical applications.Reversibility and training effects of the PdH x system have been extensively studied in References [60,61].
The in situ FMR signal responses of the 8.36 nm Co/13.50 nm Pd bilayer obtained on PLATYPUS during two hydrogen absorption and desorption cycles are shown in Figure 14.During these measurements, the microwave frequency was kept constant and equal to 3.4 GHz, and the applied field was swept across the resonance to produce the plots shown in the figure.Ensuring efficient neutron reflection from the film requires taking FMR measurements in the "film facing up" geometry (see Section 2.2.3).For this film orientation, the FMR signal is significantly reduced because of film separation from the microwave transmission line by the 0.5-mm-thick silicon substrate of the film.In order to compensate for this signal reduction, the applied-field modulation technique combined with lock-in signal detection was utilized [62].For this reason, the registered FMR traces have the shape of the first derivative of a Lorentzian instead of the Lorentzian itself.
The characteristic resonance changes its location with respect to the applied magnetic field depending on the hydrogen concentration within the Pd layer, shifting to lower applied fields upon hydrogen absorption, and shifting upwards upon hydrogen desorption.The downshift of the resonance peak implies that the PMA strength is reduced when hydrogen atoms are incorporated into the Pd crystal lattice.Reproducible responses analogous to those observed in the PNR data between the second and third hydrogen absorption and desorption cycling stages were also evident in the FMR data.

Conclusions
In summary, the combination of the scattering methods neutron reflectometry, X-ray reflectometry, high-angle X-ray diffraction and EXAFS, all performed in situ during various hydrogen charging and discharging cycles, revealed substantial new insights in understanding the detailed physics related to hydrogen absorption and desorption in thin metal films.
Our case studies of the individual metal systems Fe/Nb, W/Nb and Co/Pd revealed that such structures absorb hydrogen in large quantities.Focusing on structural property changes upon hydrogenation, for Fe/Nb and W/Nb, we found that the relative out-of-lane expansion of the Nb layers is considerably larger than the relative increase of the Nb interplanar spacing, which indicates two distinctly different mechanisms of hydrogen absorption.In polycrystalline Fe/Nb multilayers, hydrogen expands the Nb interplanar spacing in a continuous way as a function of the external hydrogen pressure.In contrast, in epitaxial W/Nb multilayers-due to a much stronger impact of the initial adherence to the substrate-the Nb lattice expansion is discontinuous and can be regarded as a structural phase transition from exclusively out-of-plane to a three-dimensionally expanded state at low and high hydrogen pressures, respectively.
The hydrogen-induced structural changes, such as layer and lattice plane expansions, that we present in this manuscript are most pronounced during the first hydrogen absorption process.The initial absorption process results in large plastic deformations as the original state of the layer thickness and lattice constants cannot be recovered during the following hydrogen desorption process.The second process and following hydrogen absorption/desorption cycles, however, show, to a large extent, reversible behaviors of layer thickness and lattice spacings.
Our studies of the effect of hydrogen absorption on the magnetic properties of our systems revealed that the magnetic exchange coupling in Fe/Nb multilayers can reversibly be switched by hydrogen charging and discharging.In the Co/Pd system, the magnetic anisotropy at the interfaces can reversibly be changed by absorption and desorption of hydrogen.These two observed electronic effects enable a novel approach for hydrogen sensing using a magnetic readout scheme.

Conclusions
In summary, the combination of the scattering methods neutron reflectometry, X-ray reflectometry, high-angle X-ray diffraction and EXAFS, all performed in situ during various hydrogen charging and discharging cycles, revealed substantial new insights in understanding the detailed physics related to hydrogen absorption and desorption in thin metal films.
Our case studies of the individual metal systems Fe/Nb, W/Nb and Co/Pd revealed that such structures absorb hydrogen in large quantities.Focusing on structural property changes upon hydrogenation, for Fe/Nb and W/Nb, we found that the relative out-of-lane expansion of the Nb layers is considerably larger than the relative increase of the Nb interplanar spacing, which indicates two distinctly different mechanisms of hydrogen absorption.In polycrystalline Fe/Nb multilayers, hydrogen expands the Nb interplanar spacing in a continuous way as a function of the external hydrogen pressure.In contrast, in epitaxial W/Nb multilayers-due to a much stronger impact of the initial adherence to the substrate-the Nb lattice expansion is discontinuous and can be regarded as a structural phase transition from exclusively out-of-plane to a three-dimensionally expanded state at low and high hydrogen pressures, respectively.
The hydrogen-induced structural changes, such as layer and lattice plane expansions, that we present in this manuscript are most pronounced during the first hydrogen absorption process.The initial absorption process results in large plastic deformations as the original state of the layer thickness and lattice constants cannot be recovered during the following hydrogen desorption process.The second process and following hydrogen absorption/desorption cycles, however, show, to a large extent, reversible behaviors of layer thickness and lattice spacings.
Our studies of the effect of hydrogen absorption on the magnetic properties of our systems revealed that the magnetic exchange coupling in Fe/Nb multilayers can reversibly be switched by hydrogen charging and discharging.In the Co/Pd system, the magnetic anisotropy at the interfaces can reversibly be changed by absorption and desorption of hydrogen.These two observed electronic effects enable a novel approach for hydrogen sensing using a magnetic readout scheme.

Figure 1 .
Figure 1.Hydrogen-induced three-dimensional isotropic expansion of the host lattice (bulk) versus one-dimensional lattice expansion (thin films).

Figure 1 .
Figure 1.Hydrogen-induced three-dimensional isotropic expansion of the host lattice (bulk) versus one-dimensional lattice expansion (thin films).

Figure 2 .
Figure 2. Schematic characteristic of the potential energy (chemical potential) of hydrogen atoms at a metal/metal-hydride/metal interface.Example: Fe (small solubility for hydrogen) and Nb (large solubility for hydrogen).The red dots indicate preferred hydrogen absorption sites.(a) Fe and Nb layers separated; (b) Fe and Nb layers in contact.

Figure 2 .
Figure 2. Schematic characteristic of the potential energy (chemical potential) of hydrogen atoms at a metal/metal-hydride/metal interface.Example: Fe (small solubility for hydrogen) and Nb (large solubility for hydrogen).The red dots indicate preferred hydrogen absorption sites.(a) Fe and Nb layers separated; (b) Fe and Nb layers in contact.

Figure 4 .
Figure 4. Simulated neutron reflectivity curves of a bulk Pd surface (red solid line) and a Pd-hydride bulk surface (blue dashed line).

Figure 4 .
Figure 4. Simulated neutron reflectivity curves of a bulk Pd surface (red solid line) and a Pd-hydride bulk surface (blue dashed line).

3 .
Results and Discussion 3.1.Nb Thin Films and Multilayers 3.1.1.Fe/Nb Multilayers In order to study the correlation between hydrogen absorption and resulting changes of the structural properties of Nb films, Fe/Nb heterostructures were measured with neutron reflectometry Metals 2016, 6, 125 9 of 20

Figure 6 .
Figure 6.Unpolarized neutron reflectivity data for several polycrystalline Fe/Nb heterostructures.The red curves represent data taken in the virgin (pre-hydrogen) state, while the blue curves represent data taken at the highest hydrogen pressure of 900 mbar.

Figure 6 .
Figure 6.Unpolarized neutron reflectivity data for several polycrystalline Fe/Nb heterostructures.The red curves represent data taken in the virgin (pre-hydrogen) state, while the blue curves represent data taken at the highest hydrogen pressure of 900 mbar.

Figure 7 .
Figure 7. PNR data (circles) and fits (black lines) for a polycrystalline Fe/Nb heterostructure.Figure adapted with permission from Ref. [41].Copyrighted by the American Physical Society.

Figure 7 .
Figure 7. PNR data (circles) and fits (black lines) for a polycrystalline Fe/Nb heterostructure.Figure adapted with permission from Ref. [41].Copyrighted by the American Physical Society.

Figure 8 .
Figure 8. Phase diagram for hydrogen absorbing Nb thin films.Reprinted from [56] with permission from Elsevier.

Figure 8 .
Figure 8. Phase diagram for hydrogen absorbing Nb thin films.Reprinted from [56] with permission from Elsevier.

Figure 9 .
Figure 9. Neutron reflectivity curves of W/Nb multilayers with different Nb layer thicknesses as indicated, measured in the as-prepared state (red circles), and at a hydrogen pressure of 900 mbar (blue triangles).

Figure 9 .
Figure 9. Neutron reflectivity curves of W/Nb multilayers with different Nb layer thicknesses as indicated, measured in the as-prepared state (red circles), and at a hydrogen pressure of 900 mbar (blue triangles).

Figure 10 .
Figure 10.Hydrogen solubility curves of epitaxial [2.6 nm W/dNb]n multilayers and Nb bulk material as indicated.

Figure 11 .
Figure 11.Cu Kα X-ray diffraction on the multilayers [2.6 nm W/10 nm Nb]60 (left); and [2.6 nm W/3 nm Nb]120 (right) as a function of increasing external hydrogen pressure, and (bottom) time-resolved measurement at 1 mbar and 6.5 mbar, respectively.The dashed lines indicate the positions of the Nb(002) and W(002) bulk reflections.

Figure 11 .
Figure 11.Cu Kα X-ray diffraction on the multilayers [2.6 nm W/10 nm Nb]60 (left); and [2.6 nm W/3 nm Nb]120 (right) as a function of increasing external hydrogen pressure, and (bottom) time-resolved measurement at 1 mbar and 6.5 mbar, respectively.The dashed lines indicate the positions of the Nb(002) and W(002) bulk reflections.

Figure 11 .
Figure 11.Cu K α X-ray diffraction on the multilayers [2.6 nm W/10 nm Nb] 60 (left); and [2.6 nm W/3 nm Nb] 120 (right) as a function of increasing external hydrogen pressure, and (bottom) time-resolved measurement at 1 mbar and 6.5 mbar, respectively.The dashed lines indicate the positions of the Nb(002) and W(002) bulk reflections.

Figure 12 .
Figure 12.Comparison between hydrogen-induced relative increase of the out-of-plane Nb(110) interplanar spacing (top); and increase of the Nb layer thickness (bottom) in W/Nb multilayers for various hydrogen pressures.

Figure 12 .
Figure 12.Comparison between hydrogen-induced relative increase of the out-of-plane Nb(110) interplanar spacing (top); and increase of the Nb layer thickness (bottom) in W/Nb multilayers for various hydrogen pressures.
the Pd layer from 1.0 nm to 0.2 nm.Due to the negative scattering length of hydrogen, these observed changes to structure and magnetism highlight the presence of hydrogen at the Co/Pd interface, which have additionally resulted in modifications to the PMA strength within the Pd lattice.

Figure 13 .
Figure 13.Comparison of R+ polarized neutron reflectivity data (circles) and fits (solid lines) for a Co/Pd bilayer before (red) and after (blue) exposure to 3.5% hydrogen gas partial pressure.

Figure 13 .
Figure 13.Comparison of R+ polarized neutron reflectivity data (circles) and fits (solid lines) for a Co/Pd bilayer before (red) and after (blue) exposure to 3.5% hydrogen gas partial pressure.

Figure 14 .
Figure 14.In situ FMR response of a 8.36 nm Co/13.50 nm Pd bilayer at various stages during the hydrogen absorption and desorption cycling PNR experiment.

Figure 14 .
Figure 14.In situ FMR response of a 8.36 nm Co/13.50 nm Pd bilayer at various stages during the hydrogen absorption and desorption cycling PNR experiment.

Table 1 .
List of samples for the Fe/Nb and W/Nb measurements.The description in the brackets is the thickness of each layer in a single bilayer.The number in subscript represents how many times each bilayer was repeated.

Table 1 .
List of samples for the Fe/Nb and W/Nb measurements.The description in the brackets is the thickness of each layer in a single bilayer.The number in subscript represents how many times each bilayer was repeated.

Table 2 .
SLD and layer thicknesses for the sample data shown in Figure7.

Table 2 .
SLD and layer thicknesses for the sample data shown in Figure7.