Abstract
Inelastic neutron scattering (INS) is uniquely sensitive to hydrogen due to its comparatively large thermal neutron scattering cross-section (82 b). Consequently, the inclusion of water in real samples presents significant challenges to INS data analysis due directly to the scattering strength of hydrogen. Here, we investigate uranyl fluoride (UO2F2) with inelastic neutron scattering. UO2F2 is the hydrolysis product of uranium hexafluoride (UF6), and is a hygroscopic, uranyl-ion containing particulate. Raman spectral signatures are commonly used for inferential understanding of the chemical environment for the uranyl ion in UO2F2, but no direct measurement of the influence of absorbed water molecules on the overall lattice dynamics has been performed until now. To deconvolute the influence of waters on the observed INS spectra, we use density functional theory with full spectral modeling to separate lattice motion from water coupling. In particular, we present a careful and novel analysis of the Q-dependent Debye–Waller factor, allowing us to separate spectral contributions by mass, which reveals preferential water coupling to the uranyl stretching vibrations. Coupled with the detailed partial phonon densities of states calculated via DFT, we infer the probable adsorption locations of interlayer waters. We explain that a common spectral feature in Raman spectra of uranyl fluoride originates from the interaction of water molecules with the uranyl ion based on this analysis. The Debye–Waller analysis is applicable to all INS spectra and could be used to identify light element contributions in other systems.
Similar content being viewed by others
Introduction
Fluorinated uranium compounds are ubiquitous in the nuclear fuel cycle but remain a largely understudied component of actinide chemistry. One such example where chemical understanding is still incomplete is uranyl fluoride (UO2F2), a byproduct of the exposure of uranium hexafluoride (UF6) to water. UF6 itself is a highly corrosive and reactive gas, and is the primary chemical form used in enrichment operations for nuclear fuel production. When exposed to atmospheric water vapor UF6 converts into uranyl fluoride (UO2F2), releasing hydrogen fluoride (HF). Recently, we showed that the hygroscopic UO2F2 adsorbs environmental water molecules in a crystal hydrate form1. These waters have a significant effect on spectroscopic signatures, in particular the Raman spectra, but accessing this information quantitatively proves difficult for a number of experimental techniques due to various intrinsic limitations based on cross sections or selection rules for scattering.
UO2F2 is comprised of layers of uranyl ions hexagonally coordinated by fluoride ligands. At ambient conditions anhydrous uranyl fluoride is a hygroscopic powder that can undergo a hydration reaction resulting in uranyl fluoride hydrate2. The transition to the hydrated state can be reversed through annealing under vacuum conditions at 115 °C1. The exact locations of retained water molecules in the anhydrous structure were heretofore unknown, but the appearance of stacking disorder in x-ray and neutron diffraction patterns offered tantalizing hints of interlayer waters2,3.
The uranyl ion is highly Raman active, with a symmetric stretching vibration around 113 meV in the anhydrous form and around 108 meV in the various hydrate structures4. Because it the most intense Raman active mode in UO2F2 and in most uranyl systems, the symmetric stretching vibration of the uranyl ion alone has been used to identify the underlying chemistry5,6. In minerals, the frequency and shape of the symmetric stretch mode is sensitive to the uranyl ligand coordination environment, including to the presence of dispersion bonded water molecules7,8,9,10,11. Nevertheless, the precise location and the influence of these remnant water molecules on the broader lattice dynamics has not been well understood.
Water dynamics within non-mineral forms of uranium solids remains at the forefront of our efforts to more fully-understand solid-phase uranium chemistry of the nuclear fuel cycle. These efforts are motivated by the need to improve analytical spectroscopic methods for the identification and characterization of non-ideal, and otherwise degraded, fuel cycle solids. The goal of this investigation is to identify the existence and understand the influence of nonremovable water on the lattice dynamics in uranyl fluoride. To achieve this, we combine inelastic neutron scattering (INS) on the in-situ vacuum-annealed anhydrous structure with density functional theory (DFT) calculations. We expect the total water content to be small in the anhydrous structure such that the observed neutron scattering spectra contain motion of both light and heavy nuclei. Further, unlike in Raman or infrared scattering, there are no selection rules governing scattering of neutrons from the phonons of the UO2F2 host lattice, allowing collection of a more complete picture of the dynamics. Nevertheless, H atoms coupled to the lattice dynamics complicate analysis, because there is no intrinsic mechanism for separating contributions in INS by mass. To tackle this problem, we offer a novel mechanism to analyze neutron scattering spectra based on application of the equipartition theorem to directly identify motion of adsorbed water in the UO2F2 structure. We further employ a self-consistent procedure for extracting the phonon density of states from the measured data enabling direct comparison to the phonon density of states from DFT12,13,14. These combined approaches allow us for the first time, to directly ascertain the dynamic affects of the adhered water in the anhydrous uranyl fluoride. The new analysis technique presented here is directly applicable to understanding the water dynamics in a wide variety of materials where residual water is retained in the crystalline structure.
Materials and Methods
DFT calculations on the anhydrous UO2F2 structure were previously performed by the authors using the VASP software package15,16,17,18,19. Anhydrous UO2F2 is spacegroup \(R\bar{3}m\) (structure shown in Fig. 1); in the hexagonal description, the linear uranyl ion is equatorially coordinated in the hexagonal a/b plane by six F ligands and stacked in a three-layer repeat motif in the c direction. Of particular importance, it was necessary to employ a van der Waals functional to appropriately model the interlayer structure, which is stabilized only by dispersion interactions20,21,22,23. An on-site Coulomb interaction was also necessary to correctly determine the bond distances and lattice parameters24,25. Phonons were calculated using density functional perturbation theory on a 4 × 4 × 4 rhombohedral supercell. A kinetic energy cutoff of 685 eV and a 4 × 4 × 4, \({\rm{\Gamma }}\)-centered k-point mesh were used. A Blöchl correction was applied to partial orbital occupations and forces were converged to 0.5 meV/Å26. Details regarding these initial calculations can be found in ref.19. Here, we are extending those results with calculations of the neutron spectroscopic response. To do this, we used the O’Climax software package27. (For details about the production of \(S(Q,\omega )\) from calculated phonon eigenvectors, see ref.27 and references therein).
Calculations were also performed with a 2 × 2 × 1 hexagonal unit cell containing an adsorbed water molecule. All computational parameters were the same as in the anhydrous structure except that a 2 × 2 × 4, \({\rm{\Gamma }}\)-centered k-point mesh was used. Phonon frequencies and eigenvectors were calculated with density functional perturbation theory after geometric relaxation. Figure 1, panel (b), shows the position of a relaxed water molecule, forming H bonds with uranyl oxygens and perturbing the linear O–U–O unit.
INS measurements were collected on the Wide-Angular Range Chopper Spectrometer (ARCS) at SNS28. Measurements were collected at 15 K with 20, 80, 175, and 700 meV incident energy neutrons. Neutrons were collected at each incident energy for approximately 4 h for a consistent beam current. Anhydrous uranyl fluoride was prepared by hydrolysis of UF6 (UF6 + 2H2O −> UO2F2 + 4HF) in a reaction chamber containing approximately 20% relative humidity (23 °C). Uranyl fluoride forms small (<1 μm) particulates from this reaction, which were subsequently collected on a series of silicon plates. The material was baked under a steady flow of 10 mL/s N2 at 150 °C for 18 h, a procedure which has been shown to produce anhydrous uranyl fluoride without decomposing the crystal structure2,29. The sample was further exposed to vacuum and baked at 150 °C in the beamline for several hours prior to data collection on ARCS to drive off additional surface-adsorbed water. Measurement of the empty, 8 mm diameter Al sample can at the same temperature and incident energies was subtracted from presented data.
Results and Discussion
In Fig. 2, we plot the background-subtracted neutron dynamic structure factor, \(S(Q,\omega )\), as collected on ARCS (top panels) for three incident neutron energies (20, 80, and 175 meV) and a comparison to the same quantity derived from the DFT calculation (bottom panels). The measured phonon spectra contains a low-energy (<50 meV) region containing crystalline phonon modes and an intermediate region near 110 meV with uranyl ion stretching vibrations.
Simulated \(S(Q,\omega )\), as produced with the O’Climax code for the same set of incident neutron energies are shown in Fig. 2 panels (d)–(f). O’Climax calculates the instruments’ energy resolution profile given input on the instrument geometry, moderator emission time, and chopper rotation frequencies. It then applies a Gaussian broadening to the points and bins them in 0.08, 0.1, and 0.2 meV bins, which were chosen to match the energy binning in the data sets. The actual energy resolution on ARCS depends on both incident and outgoing neutron energy, but for 20, 80, and 175 meV incident energy neutrons it is approximately 0.5, 2, and 4 meV at the elastic line respectively. Agreement is observed between the calculated and experimental spectra in terms of reproduction of the key peak frequencies. The symmetric and asymmetric uranyl stretching vibrations are observed at 110 and 118 meV. The frequencies of the same vibrations are calculated via DFT to be 110 and 135 meV. The peak frequency at 135 meV disagrees with the data for computational reasons described in Kirkegaard et al.19. The assignment of these peaks as the uranyl stretching vibrations is unambiguous because the DFT calculation gives the phonon eigenvectors exactly. Observation of coherent phonons is noted with peaks at Q = 3.8 and 4.9 Å−1 in panel (a) and near 5, 7, and 9 Å−1 in panel (b), which is also observed in panels (d) and (e).
A fundamental quantity derivable from \(S(Q,\omega )\) is the generalized phonon density of states (gdos), which can also be determined from the DFT calculation. Because the data is collected at several incident energies, individual data sets must be knit together using a self-consistent procedure to take advantage of both the better energy resolution of lower energy neutrons and the dynamic range of higher energy neutrons. Details of this procedure are outlined in refs12,13,14. In addition to knitting together multiple incident neutron energies, the procedure also corrects for multiple scattering effects, such that the data points in Fig. 3 represent the one-phonon density of states. Generation of this quantity is important because it can be directly compared to the DFT-calculated one-phonon density of states, given by:
where σi is the (total) neutron cross-section of each isotope (with mass Mi), \({e}_{i}(j,Q)\) is the phonon eigenvector describing the displacement of each atom along phonon branch j at wavevector Q, with energy \(\omega (j,Q)\)30. DFT calculations provide both \(\omega (j,Q)\) and \({e}_{i}(j,Q)\). Calculation of the full dispersion is done via phonopy with a 48 × 48 × 48 k-point mesh sampling of the Brillouin zone31.
Normalization of both curves is done assuming
with N representing the number of atoms in the primitive unit cell (5). We select \({E}_{{\rm{\min }}}=0\,{\rm{meV}}\) and \({E}_{{\rm{\max }}}=145\,{\rm{meV}}\), above which there is only molecular scattering. Due to the appearance of additional H atom coupling, this normalization is not correct for the measured data, but nevertheless allows a qualitative comparison to the calculated \({g}_{1}(\omega )\).
We note several key features in Fig. 3, panel (a). First, the position of the main peaks are at 15, 32, 48, 55, 112, and 122 meV. The frequency of the asymmetric uranyl stretching vibration (observed at 122 meV) is overestimated by the calculation, but its assignment is unambiguous. Second, broadening of the experimental gdos is seen near 60 meV, extending up to 80 meV, which is not modeled by the calculation. Third, the modeled data significantly underestimates the gdos in the uranyl stretching region (from 90 to 130 meV). We speculate that this additional scattering intensity in the observed gdos is due to preferential water coupling to the phonon modes involving motion of the uranyl ions, and this is explored in Section 4. Last, the measured gdos is also weaker than the calculated values at the 48 meV peak (Eu uranyl rocking mode, which includes compensatory motion of F ligands) and the 15 meV peak (A1g describing equatorial F ligand motion along the hexagonal c axis).
Panel (b) of Fig. 3 shows the partial gdos extracted from DFT calculations. Although the low-energy lattice modes (<70 meV) show a mix of O and F motion, the uranyl stretching region is almost entirely represented by O motion. Motion of U is attenuated by its mass (238 amu) such that its contribution to the measured and calculated gdos is negligible at all energies. Further projection of O and F motion onto hexagonal a/b and c axes in Fig. 4 indicate that the O motions are polarized along c at 112 and 122 meV (uranyl stretching region, Fig. 4a inset) and polarized along a/b below 60 meV (representing primarily uranyl rocking modes). Fluorine motion is predominantly polarized along a/b for the 42 and 55 meV peaks and along c for the 18 meV peak.
Panel (c) of Fig. 3 shows the partial gdos for the structure in Fig. 1(b) with an adsorbed interlayer water molecule. Due to the large neutron cross-section and low mass of H, only one significant peak representing H motion is observed in this region, near 89 meV and is shown in the inset. This phonon mode is a librational mode with H motion polarized along the a/b axis (Rz). We therefore assign the 97.3 meV peak in the data to this mode. Because the observed intensity of the 97 meV mode is low, and it originates from H motion, we further confirm that the quantity of water in this sample is small. Although H is the dominant scatterer in the adsorbed water DFT structure, it is also noteworthy that the O motion near 47 meV is broader and weaker than in the anhydrous case.
Figure 5 shows a decomposition of the uranyl stretching region and the molecular vibrational region (300–600 meV) into constituent peaks. The uranyl stretching region shown in Fig. 5a decomposes into a set of four Gaussian peaks centered at 97.3, 109.0, 112.0, and 122.0 meV, as enumerated in Table 1. Previous Raman measurements showed that the (Raman-active) symmetric stretching vibration centered at 112 meV manifests a smaller, redshifted shoulder peak19. The IR-active asymmetric uranyl stretching peak at 122 meV is in agreement with experiment32. A smaller blueshifted shoulder is observed near 129 meV, but because of its weak intensity, we do not assign it per se.
In the molecular region (Fig. 5, panel [b]), peaks can be assigned to O–H stretching modes at 370 and 430 meV as listed in Table 133. The apparent increased background is due to the presence of a hydrogen recoil peak, which adds an incoherent background at higher incident neutron energies but no vibrational information. The existence of higher energy molecular vibrations and a hydrogen recoil peak confirm the presence of water in the experimental sample. As surface-adsorbed water was removed prior to experiment, the remaining water is further evidence that the layered UO2F2 crystal supports the localization of water from the environment, either during formation or immediately afterwards1,3. In Fig. 1(b), H atoms form H bonds with uranyl oxygens at two distinct distances, O–H–O = 2.639 and 2.820 Å. As tabulated by Libowitzky, these distances would correspond to O–H stretching frequencies of about 366 and 425 meV34. Thus, both the 370 and 430 meV could correspond to O–H stretches, related to interactions with different uranyl ions. The appearance of multiple O–H stretching frequencies originating from a single water molecule has previously been observed in mineral hydrates such as natrolite and scolecite35.
A vibrational peak is observed at 534 meV, which is too high to be assignable to fundamental water vibrations. This peak is likely a combination band (O–H plus the uranyl vibrations at 112 meV). A water bending mode near 220 meV was observed but difficult to assign directly due to relatively weak signal to background. However, a similar O–H bend plus uranyl vibration is observed in Fig. 5 at 329 meV, which allows us to indirectly assign \({\nu }_{2}-{{\rm{H}}}_{2}{\rm{O}}\) to 212 meV.
Finally, we explored two alternative models with hydrogenated species. One, we considered uranyl fluoride with an interlayer OH− ion, and two, a model in which a F ligand were replaced by an OH− ion. Both models failed to reproduce key features of the data; in particular, the OH− ligand model resulted in a complete layer collapse and significant red-shifting of the uranyl stretching frequencies due to formation of interlayer O–H–F hydrogen bonds, whereas the free OH− ion model suggested the appearance of several strong peaks around 65 meV that are not observed.
Water Coupling
The observed water O–H stretching frequencies are both lower than for liquid water and fall within the range of “strong” hydrogen bonding34. In solids, the \({\nu }_{2}\) frequency is shifted to higher energy due to H bonds, which is also observed here despite the larger uncertainty36. Observed water fundamental frequencies suggest strong H bonding.
In Fig. 3 the major disagreement between theory and experiment is the mismatch in intensity between the higher energy uranyl stretching region (90–130 meV) compared to the generally good agreement below 90 meV (both frequencies and intensities). The high energy molecular region (Fig. 5 panel [b]) indicates the existence of O–H stretching frequencies, as well as the presence of a hydrogen recoil peak. Thus, despite having been vacuum annealed, anhydrous UO2F2 contains residual water. Our previous neutron diffraction and quasielastic neutron spectroscopy measurements also indicated the existence of a strongly bound set of water molecules that could not be removed without substantial (>250 °C) annealing that can result in oxidation and damage1,3. Those water molecules are an intrinsic part of the structure.
To understand the interaction of water molecules with the lattice modes, we examine the Q-dependence of the observed data. The Q-dependence of the intensity of phonon bands measured in INS experiments can be described by
where bi is the scattering length for isotope i and nj is the thermal population factor of mode j, with frequency \({\omega }_{j}\), \({e}_{i}^{j}\) is the eigenvector describing the displacement of atom i for mode j, Mi are the masses of species atom i, Ri is the position of atom i, and N is the number of atoms37. The \({W}_{i}(Q)\) are the Debye–Waller factors, which are
The mean-square displacement of atoms from their equilibrium positions is \(\langle {u}^{2}\rangle \). In a powder average, incoherent approximation, Equation 3 can be simplified:
where \(A(\omega )\) contains all information related to the intensity as a function of frequency (the scattering lengths, polarization average, thermal factor, and frequency denominator). The \(S(Q,\omega )\) experimental data are binned in 1.5 meV bins and fit to Equation 5 with two free parameters (A and u2). The critical observation in this system is the clear difference in Q-dependence of different peaks as shown in Fig. 2. The mode at 112 meV (uranyl stretching vibration) is peaked near 9.8 Å−1, whereas the 32 meV mode is centered at 16.2 Å−1. Two examples of this Q-dependence, as well as their fit to Equation 5, are shown in Fig. 6, which clearly shows a sharper attenuation in Q for the higher energy uranyl stretching mode region. For a given \(\omega \), only certain Q values are accessible based on instrument geometry and neutron kinematics, but data from all accessible Q values for each energy bin is included in each fit. The mean-square displacement, \(\langle {u}^{2}\rangle \), is inversely proportional to mass as a requirement of the equipartition theorem38. As a consequence, modes involving significant contributions from H atoms will demonstrate stronger Q attenuation.
In Fig. 7 we fit \(S(Q,\omega )\) from Fig. 2 to Equation 5, by binning \(S(Q,\omega )\) into 1.5 meV slices and extracting \(\langle {u}^{2}\rangle \) for each slice. For reference, we include the previously reported neutron-weighted density of states from Fig. 3. The strongest Q attenuation (i.e., highest Debye–Waller factor) is observed in the uranyl stretching region (90–130 meV), indicating that the observed data represent constituents with larger mean-square displacements (H). Cross-referencing this region with Fig. 4, panel (b), we note that the mode decomposition for the uranyl stretching modes is nearly entirely the motion of the O atoms along the crystallographic c axis.
Water couples most strongly to the uranyl vibrations in UO2F2 based on the following pieces of evidence: One, the calculated gdos (which includes only U, O, and F motions) is much lower than the observed gdos in the uranyl stretching region, indicating a strong scatterer at the same frequency (H) is present; two, significant Q attenuation is noted as per Fig. 7 selectively for the 90–130 meV uranyl stretching region; three, the computed partial phonon gdos in Fig. 3, panel (b) reveals that the uranyl oxygens dominate atomic motion at those frequencies, disambiguating the assignment of those modes to uranyl stretches and four, there is a significant density of vibrational states of O–H stretching plus O–U–O stretching combination bands in the observed gdos.
Because of the layered nature of the anhydrous UO2F2 structure, it is logical that the adsorbed water molecules will form H bonds with the uranyl oxygens. The above evidence indicates that this intuition is correct. A possible position of adsorbed water molecules is shown in Fig. 1, panel (b), as relaxed with DFT. This interlayer water position accommodates between two and four H-bonds with neighboring uranyl ions and fills the space between interlayer equatorial F ligands (separated by 5.6 Å). With respect to the attenuation of the measured gdos compared to the calculated gdos at 48 and 15 meV (Fig. 3(a)), recall the resolution of the partial density of states into components of motion (Fig. 4) indicate that the motion at 15 meV is entirely F atoms moving along the hexagonal c axis, and the motion at 48 meV is primarily motion of O atoms moving along the a/b axes. We observe that both of these atomic motions collide with the location of adsorbed water molecules, attenuating the gdos by a mechanism of steric hindrance39,40,41. In particular, the effective oscillator mass is increased by the coupling of O/F atoms to the H motion, reducing the expected scattering power. Previous INS measurements also identified signatures of this effect in the mean-square displacement42,43. Whereas, the atomic motions of the uranyl oxygens in the uranyl stretching region (strong H coupling) are polarized along c and not in collision with adsorbed waters. In short, a vector drawn between uranyl oxygens or ligand F atoms to the proposed water position in Fig. 1 panel (b) is collinear with the phonon polarization at the same frequencies where we observe additional attenuation in the gdos. The effect of the additional water molecule on the gdos clearly corroborates this interpretation in Fig. 3, panel (c), because the 15 and 48 meV modes are preferentially attenuated compared to the same modes in the anhydrous structure, and significant water H motion is present at those frequencies.
According to the inset in Fig. 3, panel (c), the primary observable in the gdos for adsorbed water molecules is the appearance of the Rz librational peak. The observation of this mode in the data is strong evidence for the presence of interlayer water, especially because the anhydrous UO2F2 structure has no vibrational modes to explain this mode. The model structure also directly shows coupling between adsorbed water molecules and the uranyl vibrations as evidenced by the overlap of gdos between uranyl O and water H atoms.
In Table 1 we identify two peaks at 109 and 112 meV associated with the symmetric stretching mode of the uranyl ion. Raman spectra of the same vibrational mode consistently show a similar two peak decomposition, always with a redshift shoulder1,4,44. For the symmetric mode, the vibrational frequency is proportional to the square root of the oscillator mass (16 amu for O). The existence of transient Oyl–H bonds increases the effective mass of the oscillator by 1 amu. Considering most uranyl ions are not near water molecules (as required by previous thermal gravimetric measurements), the fundamental oscillation frequency at 112 meV represents unperturbed uranyl ions2. Adding a nearby H-bonded water molecule shifts the frequency by
compared to the measured value of 109.0 meV. This phenomenon has been observed in a variety of other water-containing uranyl species, such as uranyl hydroxide, ammonium diuranate, ammonium uranyl carbonate, sodium diuranate, and uranyl acetate45,46,47. This mechanism of H-bonding damping may be responsible for the line shape of uranyl stretching vibrations. Redshift in the uranyl stretching vibrations in Fig. 3, panel (c), due to the presence of water corroborates this hypothesis. More investigation, such as careful isotopic substitution, could further improve the interpretive value of this commonly used vibrational signature. Previously, we attempted to equilibrate samples of hydrated UO2F2 with deuterated saturated salts, but INS measurements indicated only a fraction of water molecules were exchangeable1. Entrainment of D2O into UO2F2 may thus require formation in a D2O environment.
Conclusions
Using a combination of INS and DFT modeling, the vibrational spectra of anhydrous UO2F2 was explored. An enhancement of scattering intensity in the 90–130 meV uranyl stretching region compared to the modeled intensity suggested the presence of a strong scatterer, which was confirmed by the appearance of O–H stretching modes in the molecular vibrational region (>350 meV) and the appearance of a water librational mode at 97 meV. Observation of H atoms in anhydrous UO2F2 corroborated previous investigations that postulated the existence of strongly bound water1.
Detailed analysis of the Q-dependence of the intensity of \(S(Q,\omega )\) indicates preferential attenuation of intensity in the uranyl stretching region. Correlation of this attenuation with oxygen motion (indicated by the partial phonon gdos), as well as the observed scattering intensity enhancement in that region, indicates a strong and preferential coupling between adsorbed, nonremovable water molecules and the uranyl stretching vibrations (particularly the symmetric stretching mode). Steric hindrance between adsorbed water molecules and the backbone structure, as indicated by the attenuation of certain vibrational modes whose atomic motions collide with the locations of those adsorbed waters, independently substantiate this argument.
Data Availability
Data sets and DFT input files can be made available upon reasonable request to the corresponding author.
References
Miskowiec, A. et al. Structural phase transitions and water dynamics in uranyl fluoride hydrates. J. Phys. Chem. A 119, 11900–10, https://doi.org/10.1021/acs.jpca.5b09296 (2015).
Morato, F., Fulconis, J. M., Rouquérol, F. & Fourcade, R. Study of the dehydration process of uranyl difluoride hydrates stable under usual conditions of temperature, pressure and atmospheric moisture. J. Fluorine Chem. 91, 69–73, https://doi.org/10.1016/S0022-1139(98)00214-0 (1998).
Miskowiec, A. et al. Time-dependent water dynamics in hydrated uranyl fluoride. Mol. Phys. 114, 61–71, https://doi.org/10.1080/00268976.2015.1084056 (2016).
Armstrong, D., Jarabek, R. J. & Fletcher, W. H. Micro-Raman spectroscopy of selected solid UxOyFz compounds. Appl. Spectrosc. 43, 461–468, https://doi.org/10.1366/0003702894203002 (1989).
McGlynn, S. P., Smith, J. K. & Neely, W. C. Electronic structure, spectra, and magnetic properties of oxycations. III. Ligation effects on the infrared spectrum of the uranyl ion. J. Chem. Phys. 35, 105–116, https://doi.org/10.1063/1.1731876 (1961).
Bullock, J. I. Raman and infrared spectroscopic studies of the uranyl ion: the symmetric stretching frequency, force constants, and bond lengths. J. Chem. Soc. A 0, 781–784, https://doi.org/10.1039/J19690000781 (1969).
Frost, R. L. & Čejka, J. A Raman spectroscopic study of the uranyl carbonate rutherfordine. J. Raman Spectrosc. 38, 1488–1493, https://doi.org/10.1002/jrs.1796 (2007).
Frost, R. L., Čejka, J., Ayoko, G. A. & Dickfos, M. J. Raman spectroscopic study of the multi-anion uranyl mineral schroeckingerite. J. Raman Spectrosc. 38, 1609–1614, https://doi.org/10.1002/jrs.1827 (2007).
Frost, R. L., Carmody, O., Erickson, K. L., Weier, M. L. & Čejka, J. Molecular structure of the uranyl mineral andersonite—a Raman spectroscopic study. J. Mol. Struct. 703, 47–54, https://doi.org/10.1016/j.molstruc.2004.06.001 (2004).
Bonales, L. J., Colmenero, F., Cobos, J. & Timon, V. Spectroscopic Raman characterization of rutherfordine: a combined DFT and experimental study. Phys. Chem. Chem. Phys. 18, 16575–84, https://doi.org/10.1039/c6cp01510g (2016).
Stefaniak, E. A. et al. Combined SEM/EDX and micro-Raman spectroscopy analysis of uranium minerals from a former uranium mine. J. Hazard. Mater. 168, 416–423, https://doi.org/10.1016/j.jhazmat.2009.02.057 (2009).
Arnold, O. et al. Mantid—Data analysis and visualization package for neutron scattering and μSR experiments. Nucl. Instrum. Meth. A 764, 156–166, https://doi.org/10.1016/j.nima.2014.07.029 (2014).
Kresch, M., Delaire, O., Stevens, R., Lin, J. Y. Y. & Fultz, B. Neutron scattering measurements of phonons in nickel at elevated temperatures. Phys. Rev. B 75, 104301, https://doi.org/10.1103/PhysRevB.75.104301 (2007).
Lin, J. Y. Y., Islam, F. & Kresch, M. Multiphonon: Phonon density of states tools for inelastic neutron scattering powder data. J. Open Source Software 3, 440, https://doi.org/10.21105/joss.00440 (2018).
Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 47, 558–561, https://doi.org/10.1103/PhysRevB.47.558 (1993).
Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251–14269, https://doi.org/10.1103/PhysRevB.49.14251 (1994).
Kresse, G. & FurthmÜller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169–11186, https://doi.org/10.1103/PhysRevB.54.11169 (1996).
Kresse, G. & FurthmÜller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comp. Mat. Sci. 6, 15–50, https://doi.org/10.1016/0927-0256(96)00008-0 (1996).
Kirkegaard, M. C., Langford, J., Steill, J., Anderson, B. & Miskowiec, A. Vibrational properties of anhydrous and partially hydrated uranyl fluoride. J. Chem. Phys. 146, 024502, https://doi.org/10.1063/1.4973430 (2017).
Klimeš, J., Bowler, D. R. & Michaelides, A. Van der Waals density functionals applied to solids. Phys. Rev. B 83, 195131, https://doi.org/10.1103/PhysRevB.83.195131 (2011).
Román-Pérez, G. & Soler, J. M. Efficient implementation of a van der Waals density functional: application to double-wall carbon nanotubes. Phys. Rev. Lett. 103, 096102, https://doi.org/10.1103/PhysRevLett.103.096102 (2009).
Klimeš, J., Bowler, D. R. & Michaelides, A. Chemical accuracy for the van der Waals density functional. J. Phys. Condens. Mat. 22, 022201, https://doi.org/10.1088/0953-8984/22/2/022201 (2010).
Dion, M., Rydberg, H., Schöder, E., Langreth, D. C. & Lundqvist, B. I. Van der Waals density functional for general geometries. Phys. Rev. Lett. 92, 246401, https://doi.org/10.1103/PhysRevLett.92.246401 (2004).
Anisimov, V. I., Zaanen, J. & Andersen, O. K. Band theory and Mott insulators: Hubbard U instead of stoner. I. Phys. Rev. B 44, 943–954, https://doi.org/10.1103/PhysRevB.44.943 (1991).
Dudarev, S. L., Botton, G. A., Savrasov, S. Y., Humphreys, C. J. & Sutton, A. P. Electron-energy-loss spectra and the structural stability of nickel oxide: An LSDA+U study. Phys. Rev. B 57, 1505–1509, https://doi.org/10.1103/PhysRevB.57.1505 (1998).
Blöchl, P. E., Jepsen, O. & Andersen, O. K. Improved tetrahedron method for Brillouin-zone integrations. Phys. Rev. B 49, 16223–16233, https://doi.org/10.1103/PhysRevB.49.16223 (1994).
Cheng, Y. Q., Daemen, L. L., Kolesnikov, A. I. & Ramirez-Cuesta, A. J. Simulation of inelastic neutron scattering spectra using OCLIMAX. J. Chem. Theory Comput. 15, 1974–1982, https://doi.org/10.1021/acs.jctc.8b01250 (2019).
Abernathy, D. L. et al. Design and operation of the wide angular-range chopper spectrometer ARCS at the Spallation Neutron Source. Rev. Sci. Instrum. 83, 015114, https://doi.org/10.1063/1.3680104 (2012).
Ferris, L. M. & Baird, F. G. Decomposition of uranyl fluoride between 700 °C and 950 °C. J. Electrochem. Soc. 107, 305–308, https://doi.org/10.1149/1.2427685 (1960).
Fultz, B. Vibrational thermodynamics of materials. Prog. Mater. Sci. 55, 247–352, https://doi.org/10.1016/j.pmatsci.2009.05.002 (2010).
Togo, A. & Tanaka, I. First principles phonon calculations in materials science. Scr. Mater. 108, 1–5, https://doi.org/10.1016/j.scriptamat.2015.07.021 (2015).
Hoekstra, H. R. Uranium-oxygen bond lengths in uranyl salts: Uranyl fluoride and uranyl carbonate. Inorg. Chem. 2, 492–496, https://doi.org/10.1021/ic50007a016 (1962).
Crupi, V., Longo, F., Majolino, D. & Venuti, V. Vibrational properties of water molecules adsorbed in different zeolitic frameworks. J. Phys.-Condens. Mat 18, 3563–3580, https://doi.org/10.1088/0953-8984/18/15/004 (2006).
Libowitzky, E. Correlation of O–H stretching frequencies and O–H–O hydrogen bond lengths in minerals. Monatshefte fur Chemie 130, 1047–1059, https://doi.org/10.1007/BF03354882 (1999).
Kolesov, B. A. & Geiger, C. A. Behavior of H2O molecules in the channels of natrolite and scolecite: A Raman and IR spectroscopic investigation of hydrous microporous silicates. Am. Mineral. 91, 1039–1048, https://doi.org/10.2138/am.2006.2049 (2006).
Falk, M. The frequency of the H–O–H bending fundamental in solids and liquids. Spectrochim. Acta 40A, 43–48, https://doi.org/10.1016/0584-8539(84)80027-6 (1984).
Taraskin, S. N. & Elliott, S. R. Connection between the true vibrational density of states and that derived from inelastic neutron scattering. Phys. Rev. B 55, 117–123, https://doi.org/10.1103/PhysRevB.55.117 (1997).
Colella, R., Zhang, Y., Sutter, J. P., Ehrlich, S. N. & Kycia, S. W. Debye-Waller factors in a quasicrystal. Phys. Rev. B 63, 014202, https://doi.org/10.1103/PhysRevB.63.014202 (2000).
Li, C. W. et al. Orbitally driven giant phonon anharmonicity in SnSe. Nature Physics 11, 1063–1069, https://doi.org/10.1038/nphys3492 (2015).
Skoug, E. J. & Morelli, D. T. Role of lone-pair electrons in producing minimum thermal conductivity in nitrogen-group chalcogenide compounds. Phys. Rev. Lett. 107, 235901, https://doi.org/10.1103/PhysRevLett.107.235901 (2011).
Keppens, V. et al. Localized vibrational modes in metallic solids. Nature 395, 876, https://doi.org/10.1038/27625 (1998).
Frick, B. & Richter, D. Change of the vibrational dynamics near the glass transition in polyisobutylene: Inelastic neutron scattering on a nonfragile polymer. Phys. Rev. B 47, 14795–14804, https://doi.org/10.1103/PhysRevB.47.14795 (1993).
Line, C. M. B. & Kearley, G. J. An inelastic incoherent neutron scattering study of water in small-pored zeolites and other water-bearing minerals. J. Chem. Phys. 112, 9058–9067, https://doi.org/10.1063/1.481518 (2000).
Kirkegaard, M. C., Miskowiec, A., Ambrogio, M. W. & Anderson, B. B. Evidence of a nonphotochemical mechanism for the solid-state formation of uranyl peroxide. Inorg. Chem. 57, 5711–5715, https://doi.org/10.1021/acs.inorgchem.8b00512 (2018).
Ho, D. M. et al. Raman spectroscopy of uranium compounds and the use of multivariate analysis for visualization and classification. Forensic Sci. Int. 251, 61–68, https://doi.org/10.1016/j.forsciint.2015.03.002 (2015).
Berlizov, A., Ho, D. M. L., Nicholl, A., Fanghänel, T. & Mayer, K. Assessing hand-held Raman spectrometer FirstDefender RM for nuclear safeguards applications. J. Radioanal. Nucl. Ch. 307, 285–295, https://doi.org/10.1007/s10967-015-4160-z (2015).
Palacios, M. L. & Taylor, S. H. Characterization of uranium oxides using in situ micro-Raman spectroscopy. Appl. Spectrosc. 54, 1372, https://doi.org/10.1366/0003702001951057 (2000).
Acknowledgements
A portion of this research used resources at the Spallation Neutron Source, a Department of Energy, Office of Science User Facility operated by the Oak Ridge National Laboratory. This research also used resources of the Compute and Data Environment for Science (CADES) at Oak Ridge National Laboratory, which is supported by the U.S. Department of Energy, Office of Science under contract number DE-AC05-00OR22725. A.E.S. would like to thank the Department of Homeland Security Domestic Nuclear Detection Office National Technical Nuclear Forensics Center for the fellowship funding. The authors would like to thank an anonymous reviewer for significantly improving the spectral assignment based on their comments.
Author information
Authors and Affiliations
Contributions
A.M. performed the experiments, analyzed the INS spectra, prepared all figures and the text of the manuscript. J.N. calculated the phonon densities of state from the measured data and significantly contributed to the analysis approach. M.C.K. performed the DFT calculations. A.E.S. performed the full spectral modeling. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Miskowiec, A., Niedziela, J.L., Kirkegaard, M.C. et al. Analysis of Water Coupling in Inelastic Neutron Spectra of Uranyl Fluoride. Sci Rep 9, 10476 (2019). https://doi.org/10.1038/s41598-019-46675-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-019-46675-x
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.