Oxidation of Ni-Cr and Ni-Cr-Mo was studied in operando with near ambient pressure x-ray photoelectron spectroscopy in the Cabrera-Mott regime. The oxidation temperature was 200°C—a severely diffusion-limited regime. The near-surface alloy is Cr-enriched after the reduction of native oxide in vacuum, and especially so for Ni-15Cr-6Mo. Mo-cations are integrated into the oxide and Mo(VI) dominates at the surface. The surface chemistry-driven promotion of chromia by Mo predicted by theory is negated by the limited surface diffusion of reactants. Preoxidation processing is proposed to control the oxide properties for the use of Ni-Cr superalloys at low temperatures.
INTRODUCTION
NiCr alloys are highly coveted for various technical applications as they exhibit a combination of superior corrosion resistance and high-temperature strength.1-4 These superalloys are used in a number of challenging environments including the nuclear fuel industry, the petrochemical industry, as well as for marine applications.5-7 The Ni base provides the mechanical strength for the system while Cr establishes the corrosion/oxidation resistance upon the formation of a passive Cr2O3 film.1,3-4,8 A critical concentration of 10 wt% to 15 wt% Cr is necessary for the Cr2O3 layer to protect the underlying alloy from further degradation. The properties of NiCr alloys are further enhanced through the addition of minor alloying elements (MAE) such as Mo and W.9-14 For example, alloying with Mo reduces catastrophic breakdown events, namely localized crevice corrosion, thereby improving overall corrosion resistance.8,10-11,14-16 In high-temperature alloys, molybdenum is added to enhance mechanical strength by solid solution hardening.17 However, Mo forms volatile oxides which can rupture the passive chromium oxide that is formed in these alloys.18-19 The spallation becomes relevant at temperatures greater than 600°C and severely degrades the operational lifetime of components.
The mechanisms which are responsible for the beneficial impact of MAEs on the quality of the oxide layer remain a subject of discussion in the literature.13 From an aqueous corrosion viewpoint, the solute vacancy interaction model suggests Mo has a direct effect on the transport of point defects and influences oxide film growth.14 The model argues Mo4+ and Mo6+ cations become incorporated into Cr2O3 films as substitutional defects and electrostatically interact with cation vacancies and reduce their mobility. The decreased mobility slows the formation of cation vacancy condensates which have been proposed to initiate pit formation.14,16 To date, most of the work on MAEs and their role in enhancing oxide layer stability13-15 has focused on processes in the oxide and mostly neglected the initial alloy-oxygen interaction, and their role in oxide growth. Recent studies give an atomic scale view into the oxidation process on the reaction of O2 with Ni-Cr alloys using operando transmission electron microscopy, and scanning tunneling microscopy and spectroscopy (STM/STS).20-27 In those references the competition between Ni and Cr oxidation is illustrated at temperatures, where the oxidation process is kinetically limited and below the volatilization temperature of Mo-oxides. Atomic scale imaging using STM and operando studies with near ambient pressure-x-ray photoelectron spectroscopy (NAP-XPS) have been extended to compositionally more complex alloys such as Ni-Cr-W and stainless steels.22-23,25,28 It was shown experimentally that the addition of W to Ni-Cr alloys enhances chromia formation for oxidation at T > 400°C. The experimental results are supported by density functional theory (DFT) calculations, which identify the favorable site for O-adsorption at W-Cr bridge sites.22 This agrees with the DFT study by Samin and Taylor29 who investigated oxygen adsorption for Ni-Cr-Mo alloys and reported an increased O adsorption energy also for Cr adjacent to Mo atoms. The oxygen adsorption energy was recently introduced as a descriptor to predict competition between oxidation processes on compositionally complex alloy surfaces.30
Discrepancies exist in the literature on the effect of Mo addition when alloys are oxidized at temperatures below 500°C. Greeff, et al.,31 observed no influence of Mo on the oxidation of Fe-9Cr-1Mo steel after 1,000 L oxidation at 157°C to 758°C. However, Mathieu and Landolt32 report at 384°C that the oxidation rate of Fe-24Cr-11Mo alloy is lower than the corresponding alloy with Mo. It was observed that Mo accumulates at the oxide/metal interface in hexavalent and metallic states and this barrier layer slows down the oxidation rate. Montemor, et al.,33 also observed the presence of Mo at the interface and in the internal chromium oxide layer on stainless steels oxidized at 250°C to 450°C, whereas the chromium oxide contribution increased with the addition of Mo to stainless steel. All literature reports for low T oxidation study thicker oxide scales, whereas the work presented here focuses on the Cabrera-Mott regime with oxide thicknesses in and below the 10 nm range.
The current work observes oxide evolution for different NiCr alloys including Ni-5Cr, Ni-15Cr, Ni-30Cr, and Ni-15Cr-6Mo by wt% using NAP-XPS. The oxidation temperature is 200°C and the diffusion of reactants from bulk to surface is therefore suppressed.34-36 Previous work on the nucleation and growth of NiO islands on Ni-5Cr(100) surfaces at 300°C combined STM experiments and modeling based on the Johnson-Mehl-Avrami-Kolmogorov equation as a framework.37 The capture radius for adsorbed oxygen atoms is only about 0.4 nm from the edge of an oxide island which confirms the locality of this process and the limits posed on the diffusion length. The oxide islands do not yet display the characteristic NiO structure which is clearly visible at higher oxidation temperatures. The island structure is closer to a diffusion-limited aggregation growth mode where reduced ad-atom mobility prevents reaching the thermodynamically favored crystalline structure.20-21,38
In our present study, the oxidation reaction was probed in real time—in operando—by using NAP-XPS. The oxide growth was performed at p(O2) = 0.0133 mbar and this positions our work in the Cabrera-Mott regime. The near-surface composition of the alloy is replicated in the oxide composition and Mo is integrated predominantly as Mo(IV) and Mo(VI) at the alloy/gas interface. Alloying with Mo increases the near-surface Cr inventory determined by the reduction of the native oxide by annealing prior to the oxidation experiment.
EXPERIMENTAL PROCEDURE
The NAP-XPS experiments were performed at Brookhaven National Laboratory (BNL) using the ambient pressure photoelectron spectroscopy (APPES) endstation at the in situ and operando soft x-ray spectroscopy (IOS/23-ID-2) beamline at the National Synchrotron Light Source II (NSLS-II). Details about the beamline and endstation are described elsewhere.39 The NiCr samples were cast into rectangular prisms, cold-worked to between 22% and 24% strain, and then recrystallized at 900°C for 1 h, with the exception of the NiCrMo alloy, which was recrystallized at 1,300°C. The samples were cast by John Perepezko, University of Wisconsin-Madison.22 Energy dispersive spectroscopy (EDS) was used to confirm sample homogeneity and composition. The samples were then polished to a mirror finish and cut into suitable dimensions for the XPS sample holder.
(a) Tabulated compositions for the bulk alloy measured with EDS, and the surface composition of the alloy after annealing and prior to the oxidation experiment measured with XPS. Due to the sample heterogeneity, the XPS composition cannot be corrected for mean free path variation between Ni and Cr core levels. Representative core level spectra of the clean alloy Ni-15Cr: (b) Ni2p3/2 and (c) Cr2p3/2.
(a) Tabulated compositions for the bulk alloy measured with EDS, and the surface composition of the alloy after annealing and prior to the oxidation experiment measured with XPS. Due to the sample heterogeneity, the XPS composition cannot be corrected for mean free path variation between Ni and Cr core levels. Representative core level spectra of the clean alloy Ni-15Cr: (b) Ni2p3/2 and (c) Cr2p3/2.
Composition for all samples as a function of oxidation time at T = 200°C, p(O2) = 0.01 Torr. Ni, Cr, Mo metal, and corresponding oxide contributions are obtained from a fit of the core levels. The sum of all contributions is set to 1, and the signal intensity reflects the change in the relative contributions of the constituents over time. Note that the time scales differ between oxidation experiments. (a) Ni-5Cr, (b) Ni-15Cr, (c) Ni-30Cr, and (d) Ni-15Cr-5Mo.
Composition for all samples as a function of oxidation time at T = 200°C, p(O2) = 0.01 Torr. Ni, Cr, Mo metal, and corresponding oxide contributions are obtained from a fit of the core levels. The sum of all contributions is set to 1, and the signal intensity reflects the change in the relative contributions of the constituents over time. Note that the time scales differ between oxidation experiments. (a) Ni-5Cr, (b) Ni-15Cr, (c) Ni-30Cr, and (d) Ni-15Cr-5Mo.
The spectra were analyzed with the fitting software, KolXPD† (XPS analysis software kolibri.net). All spectra were fit with a Shirley background, a Doniach-Sunjic function for the Ni(0), Cr(0), and Mo(0) metal peaks to account for their inherent asymmetry,41-46 and all oxide peaks were fit using a Voigt function. The multiplet splitting for Ni2p and Cr2p oxide peaks is described by Biesinger, et al.,41-45 and a constrained fit is used which follows the published reference spectra. The Ni2p1/2 core level was not included in our measurements, and preference was given to achieving shorter data acquisition time to more fully capture the oxidation process. A detailed description of the fit process is included in the Supplemental Materials S1. The Au4f7/2 peak of an Au foil mounted on the sample holder was positioned at a binding energy of 84.0 eV for energy calibration. The core levels were corrected by their respective cross sections47 for all calculations of composition.
RESULTS
The spectra and composition of the pristine alloys are summarized in Figure 1, which includes the bulk and surface composition of the alloys as measured by EDS and XPS at 800°C (600°C), and 200°C immediately prior to the start of the oxidation experiment. Figures 1(b) and (c) show the core level spectra of the pristine alloy surface for Ni-15Cr as a representative example. The Ni2p3/2 and Cr2p3/2 core levels are positioned at 853.8 eV and 574.4 eV, respectively, and the satellite accompanying Ni2p3/2 is readily apparent on the high binding energy side. The oxygen concentration on the surface is <1% and no oxide contributions are detected in the Ni and Cr core level spectra. In Figure 1(c), the XPS measurements reveal a Cr-enriched alloy surface for all alloys after oxide reduction, which is retained to a lesser degree after cooling down to 200°C with the exception of Ni-15Cr. This alloy shows a near-surface Cr depletion compared to the bulk composition, which is tentatively attributed to the presence of more rapid diffusion pathways defined by microstructure, or thermal history during and after casting. Future experiments using temperature-programmed reduction of native oxides will clarify the interplay between native oxide reduction and microstructure, which ultimately defines the near-surface composition.48-52
(a) Evolution of the oxide layer composition with Ni-oxide and Cr-oxide contributions for all alloys. The measurement of Ni-5Cr was interrupted for technical reasons. Examples of core level fits (b) Ni2p3/2 and (c) Cr2p3/2. Red symbols: measured spectra, black line: envelope of the fit, broken black line: Shirley background. (b) Ni 2p3/2: blue lines—NiO multiplet, green line—Ni3+ or defect peak, no alloy contribution. (c) Cr2p3/2: blue lines—Cr2O3 multiplet, green line—Cr(0), and dark line marking a low-intensity peak on the high binding energy—not assigned.
(a) Evolution of the oxide layer composition with Ni-oxide and Cr-oxide contributions for all alloys. The measurement of Ni-5Cr was interrupted for technical reasons. Examples of core level fits (b) Ni2p3/2 and (c) Cr2p3/2. Red symbols: measured spectra, black line: envelope of the fit, broken black line: Shirley background. (b) Ni 2p3/2: blue lines—NiO multiplet, green line—Ni3+ or defect peak, no alloy contribution. (c) Cr2p3/2: blue lines—Cr2O3 multiplet, green line—Cr(0), and dark line marking a low-intensity peak on the high binding energy—not assigned.
Figure 3(a) illustrates the interplay between NiO and chromia contributions, and Figures 3(b) and (c) are representative fits for Ni2p3/2 and Cr2p3/2 core levels close to the endpoint of the oxidation experiment. The Cr-oxide is identified as chromia by core level fitting, and detailed information about the fit process is given in the Supplemental Material S1. The chromia contributions at the first oxidation step are about 12% for Ni-5Cr and Ni-15Cr, and significantly larger for Ni-30Cr and Ni-15Cr-Mo with 53% and 34%, respectively. The error bar in the composition calculation is about 5%. The Mo-oxide contribution is not included in this graph and will be discussed separately. The chromia signal is reduced by 10% to 15% over the course of the oxidation time which can be linked to NiO overgrowing chromia islands formed early in the oxidation.
Oxidation of Cr is thermodynamically favorable,2-4,54 and the dominance of NiO at low oxidation temperatures is attributed to a lower nucleation barrier due to epitaxial lattice matching at the alloy/oxide interface and sluggish bulk diffusion of Cr, and consequently reduced reactant supply.34-36,55 Earlier work on NiCr oxidation with STM also illustrates the inherent morphological complexity which is superimposed on the chemical complexity seen in XPS: NiO tends to grow in a layer-by-layer growth mode driven by interfacial, epitaxial alignment between Ni and NiO lattices, while chromia prefers an island growth mode with small clusters or nodules with high aspect ratios.20-21
Analysis of the Ni2p3/2 spectra shown in Figure 3(b) revealed an additional component at a binding energy of 855.0 eV, which is tentatively assigned to Ni(III) ions. Fit examples with and without this component are included in the Supplemental Material S1, and it is evident that the associated spectra cannot be described by including spinel NiCr2O4 multiplets (Supplemental Material in Volders, et al.,22 for details). Figure 3(c) illustrates the excellent quality of the chromia phase for Ni-15Cr-6Mo, which is indistinguishable from the pristine chromia reference described in the literature. The high quality of chromia on the ternary alloy enhances the protective function of the oxide layer.
Ni(III) contribution for all alloys as a function of oxidation time. The sum for NiO+Ni(III) is set to 100%. The fit process is described in the Supplemental Material S1.
Ni(III) contribution for all alloys as a function of oxidation time. The sum for NiO+Ni(III) is set to 100%. The fit process is described in the Supplemental Material S1.
The concentration of Ni(III) ions is constant within the error bar throughout the oxidation process, except for two data points for Ni-15Cr-6Mo at the start of the oxidation time, and Ni5Cr. This could indicate a different defect population in this alloy, but a conclusive interpretation would require observing shorter timescales of oxidation not accessible in the present experiment. The overall trend in Ni(III) contribution implies that Ni(III) ions are either distributed in the NiO phases, or at the chromia-NiO grain boundaries and are not unique to the alloy/oxide interface. The associated vacancies might serve as precursors to the formation of Kirkendall voids.26
Attenuation of the metal signal for all alloys after the first oxidation step. (a) Metal signal as a function of oxidation time includes all M(0) contributions—Ni(0) for Ni-5Cr, Ni(0) and Cr(0) for Ni-15Cr and Ni-30Cr, and Ni(0) and Mo(0) for Ni-15Cr-6Mo. Elements not listed were not detected. The error bar reflects the precision of the fit process. (b) Tabulated values of the metal signal at the first oxidation step recorded after 360 s exposure. The initial M(0) concentration is 100%.
Attenuation of the metal signal for all alloys after the first oxidation step. (a) Metal signal as a function of oxidation time includes all M(0) contributions—Ni(0) for Ni-5Cr, Ni(0) and Cr(0) for Ni-15Cr and Ni-30Cr, and Ni(0) and Mo(0) for Ni-15Cr-6Mo. Elements not listed were not detected. The error bar reflects the precision of the fit process. (b) Tabulated values of the metal signal at the first oxidation step recorded after 360 s exposure. The initial M(0) concentration is 100%.
The total metal signal includes only Ni(0) for Ni-5Cr, Ni(0) and Cr(0) for Ni-15Cr, and Ni-30Cr, and Ni(0) and Mo(0) but no Cr(0) for Ni-15Cr-6Mo. The total metal signal as shown in Figure 5 starting at the second and all following oxidation steps, only includes the elements which are detected in their metallic state. This number excludes the oxide (cationic) contributions to the respective core level. For Ni-5Cr, no Cr(0) is detected anymore, and the same holds for Ni-15Cr-6Mo. For the two other alloys, we still can see a Cr(0) signal. For oxidation times >1,500 s, more than 80% of the respective alloy signal is either Cr(0) (Ni-15Cr, Ni-30Cr) or Mo(0) (Ni-15Cr-6Mo) with Ni(0) as a balance to 100%. This is in contrast to oxidation experiments at 500°C, where the alloys were completely depleted in Cr(0) and Mo(0) in the near-interface region.22 However, if the overlayer (oxide) thickness is close to the photoelectron information depth the element-specific variation in IMFP has to be considered in the interpretation. The imfps in the alloy for hν = 1,120 eV are as follows: Ni2p—0.6 nm, Cr2p—–1.0 nm, and Mo3d—1.4 nm. All values were calculated with Quases (Quases-IMFP-TPP2M Ver.3.0†) which follows the method described by Tanuma, et al.71 The IMFP values for oxides are about 10% to 20% larger than for metals.
Once the detection limit is approached for Ni(0), the element with the smallest IMFP, the alloy composition cannot be measured anymore. The detection limit is set here by the ability to separate Ni(0) from the NiO contribution to the core level and is about 1% to 1.5%. Once this limit is reached Cr(0) and Mo(0) with their larger IMFPs and consequently larger intensity will disproportionally contribute to the alloy signal. An experimental solution can be to measure each element with a different photon energy to adjust all IMFP’s to the same value as done in the work by Larsson, et al.,72 but this comes at the cost of acquisition rate. In conclusion, we can confirm that Cr(0) is present at the alloy/oxide interface for Ni-15Cr and Ni-30Cr, but not for Ni-5Cr and Ni-15Cr-Mo; and Mo(0) is present at the Ni-15Cr-6Mo interface. Quantitative information on the near-interface alloy composition cannot be obtained, and enrichment factors as defined by Castle and Asami cease to be informative.73 Cr(0) or Mo(0) are not fully depleted at the interface for select alloys which might be explained by a kinetic bottleneck, where transport through the alloy to the interface is more rapid than the oxidation rate.
Mo3d core level at select times during the oxidation process. Red symbols: measured spectra, black line: envelope of the fit, thin black line: Shirley background, green lines: Mo(0), blue lines: Mo(VI)—MoO3, and purple lines: Mo(IV)—MoO2. (a) pristine alloy, (b) and (c) M(0) and oxides, (d) photon energy of 850 eV, and (e) photon energy of 500 eV to probe the near-surface region.
Mo3d core level at select times during the oxidation process. Red symbols: measured spectra, black line: envelope of the fit, thin black line: Shirley background, green lines: Mo(0), blue lines: Mo(VI)—MoO3, and purple lines: Mo(IV)—MoO2. (a) pristine alloy, (b) and (c) M(0) and oxides, (d) photon energy of 850 eV, and (e) photon energy of 500 eV to probe the near-surface region.
(a) Contributions of different oxidation states to the Mo signal as a function of oxidation time. (b) Oxidation state contributions as a function of photon energy and oxidation time. The variation in photon energy yields a depth profile with 500 eV as the most surface-sensitive measurement. The oxidation time for each measurement is indicated by the numbers on the centerline.
(a) Contributions of different oxidation states to the Mo signal as a function of oxidation time. (b) Oxidation state contributions as a function of photon energy and oxidation time. The variation in photon energy yields a depth profile with 500 eV as the most surface-sensitive measurement. The oxidation time for each measurement is indicated by the numbers on the centerline.
DISCUSSION
Comparison of surface composition and oxide composition for Cr. The first oxidation step measured at 360 s and the one at 2,500 s are compared, and the data are also presented in Figures 1 and 4.
The pre-Cabrera-Mott regime which represents the early stages of oxidation, was studied in more detail in our recent publication22 on the oxidation of Ni-Cr and Ni-Cr-W alloys (T = 400°C and 500°C, p(O2) = 10−6 Torr), where the addition of W favored chromia formation by a large margin. This is explained by preferred O-adsorption at Cr-sites adjacent to W as determined by DFT calculations.29,76-79 The oxygen adsorption energy was identified as a suitable descriptor for the affinity for oxidation of a specific element and predictions for Ni-Cr-Mo mirror our work on Ni-Cr-W.29 The favored adsorption sites between (W, Mo) and Cr present a facile pathway to chromia formation.
The direct comparison of the Ni-Cr-W experiments with our present work has to be approached with caution: first, we intentionally target here an oxidation parameter space that represents very limited reactant mobility, second, the higher p(O2) accelerates oxide growth rates and opens the Cabrera-Mott regime to direct observation but limits our view of the early oxidation steps. The experimental trade-off is that we worked close to the limits of information depth for XPS, which will be remedied in future work by the use of higher photon energies. The gap in processing conditions, both with respect to p(O2) and oxidation temperature will be filled in future work.
Mo is incorporated into the oxide layer as Mo(VI), which is enriched at the oxide/gas interface and smaller contributions of Mo(IV) are present throughout the oxide. The rapid Cr-depletion of the alloy seen at higher oxidation temperatures is not replicated and Cr, and Mo, are still present at the alloy/oxide interface—a phenomenon that might be attributed to a diffusion rate mismatch at the interface. The oxide layer is nearly continuous after the first oxidation step as judged by the attenuation of the alloy signal, and we therefore observe the Cabrera-Mott oxidation regime.
In the low-temperature regime probed here, the limiting factor that controls oxide composition is not set anymore by alloy surface chemistry but by the availability of reactants. This is also evident in the comparison of Ni-15Cr, with its unusually low Cr inventory, and the other alloys: the oxide growth is nearly identical when Ni-15Cr and Ni-5Cr are compared and both alloys have nearly the same near-surface alloy composition prior to the oxidation experiment. The oxide composition can therefore be controlled by thermal cycling or other processing protocols that determine near-surface compositions: an oxide layer dominated by chromia can be achieved if the pretreatment of the alloy leads to near-surface Cr-enrichment. The Mo concentration is not affected by the thermal reduction step, and Mo is integrated in the oxide in concentrations similar to the bulk.37
The Supplemental Material includes STM images of oxide islands recorded after oxidation of a Ni-10Cr-6Mo(100) surface at 300°C which illustrate the relatively poor crystallinity of the oxide islands through the combined action of limited reactant diffusion and a short oxygen capture radius (Supplemental Material Figure S3).37 The significant concentration of defects, at least at the island perimeter, seems on first sight to contradict the low concentration of Ni3+ defects cited in Figure 4. However, the defect inventory is most certainly much more diverse and numerous types of defects are present which are not captured in XPS, and the Ni3+ concentration is measured for much thicker oxide layers. A direct quantitative comparison is therefore not valid.
We propose that elemental reaction steps in the surface chemistry are suppressed: while Oads depends only on the local adsorption sites on the alloy, the nucleation and sustained growth of chromia requires the diffusion of Cr atoms. At higher temperatures, chromia nucleation and growth proceed in all likelihood close to the primary oxygen adsorption sites: Cr-W-O or Cr-Mo-O. At 200°C, on the other hand, the severely limited surface diffusion of Cr (and Mo) removes the benefit of the Cr-Mo-O site to favor chromia over NiO. As a consequence, the oxide composition only depends on the local environment at the alloy surface and scales with near-surface alloy composition. The computational description of MAE effects therefore has to include an additional descriptor that relates to surface and near-surface diffusion coefficients.
CONCLUSIONS
In the regime of very limited diffusion the near-surface composition of the alloys is imprinted on the oxide. Alloy processing prior to oxidation can be used to modulate the near-surface composition and, in our work, vacuum annealing removed the native oxide resulting in a pristine alloy surface. Even for variable bulk composition shown in the comparison of Ni-5Cr and Ni-15Cr, a near identical near-surface composition can be reached and determines chromia content in the oxide as seen in Figure 8. Cr enrichment in the near-surface region triggered by processing (annealing) is particularly high for Ni-15Cr-6Mo; an observation also made in prior work for Ni-Cr-W alloys and attributed to accelerated diffusion and segregation of Cr at the surface of the ternary alloy. Experimental evidence from our work supports these ideas, but a quantitative study of diffusion coefficients and paths is still necessary for confirmation. Hence, the near-surface composition in the alloy is reflected in the oxide layer and surface chemistry becomes of minor importance in this regime. As a consequence, alloy processing that favors a near-surface Cr-enrichment can be used to seed a chromia-rich oxide layer, and enhance the protective function. On the other hand, we cannot rely on knowledge of the bulk composition to uniquely predict oxide composition and functionality.
Trade name.
ACKNOWLEDGMENTS
We gratefully acknowledge the support of Prof. John Perepezko, University of Wisconsin-Madison, who cast the alloys used in this work. This work was supported by the Office of Naval Research Multi-University Research Initiative “Understanding Atomic Scale Structure in Four Dimensions to Design and Control Corrosion Resistant Alloys” on Grant No. N00014-14-1-0675 under the direction of Dr. David Shifler. Reinke acknowledges the support by NSF DMR Metals and Metallic Nanostructures Award No. 2004326. This research used resources of the 23-ID-2 (IOS) beamline of the National Synchrotron Light Source II, a U.S. Department of Energy (DOE) Office of Science User Facility operated for the DOE Office of Science by Brookhaven National Laboratory under Contract No. DE-SC0012704. The support for V.A.A. to conduct this work by Prof. E. Opila (UVa) is gratefully acknowledged. The experimental data from this report will be made available if a reasonable request is made to the corresponding author.
AUTHOR CONTRIBUTIONS
G.R. and V.A.A. performed the experiments with the guidance of I.W. and A.H. The data analysis was performed by C.V., G.R., and P.R. The results were discussed by all authors. P.R. and C.V. wrote the manuscript which was proofread and approved for submission by all authors. P.R. conceived and led the study.
DECLARATION OF COMPETING INTEREST
The authors declare no competing financial interest.