Next Article in Journal
Energy, Economic, and Environmental Performance of a Single-Family House in Chile Built to Passivhaus Standard
Next Article in Special Issue
Challenges in the Long-Term Behaviour of Highly Radioactive Materials
Previous Article in Journal
Sustainable Alternate Materials for Concrete Production from Renewable Source and Waste
Previous Article in Special Issue
Long-Term Aging of Chernobyl Fuel Debris: Corium and “Lava”
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Plutonium-Doped Monazite and Other Orthophosphates—Thermodynamics and Experimental Data on Long-Term Behavior

1
Geology Department, Moscow State University, 119991 Moscow, Russia
2
V.G. Khlopin Radium Institute, 194021 St. Petersburg, Russia
3
A.F. Ioffe Institute, 194021 St. Petersburg, Russia
4
A.N. Frumkin Institute of Physical Chemistry and Electrochemistry RAS, 119071 Moscow, Russia
*
Authors to whom correspondence should be addressed.
Sustainability 2021, 13(3), 1203; https://doi.org/10.3390/su13031203
Submission received: 15 December 2020 / Revised: 16 January 2021 / Accepted: 20 January 2021 / Published: 24 January 2021
(This article belongs to the Special Issue Challenges in the Long-Term Behaviour of Highly Radioactive Materials)

Abstract

:
The paper consists of two main parts: a microscopic and spectroscopic investigation of the single crystal of 17-year-old 238Pu-doped Eu-monazite, and a theoretical calculation of the properties of several structural types of orthophosphates. It is shown that actinide-doped monazite is prone to the formation of mechanically weak, poorly crystalline crust, presumably consisting of rhabdophane. Its formation is likely promoted by the formation of peroxides and, potentially, acidic compounds, due to the radiolysis of atmospheric moisture. The calculations of mixing the enthalpies and Gibbs energies of binary solid solutions of Pu and rare earth element (REE) phosphates that were performed for the principal structural types—monazite, xenotime, rhabdophane—show that, in the case of light REEs, the plutonium admixture is preferentially redistributed into the rhabdophane. This process strongly affects the behavior of actinides, leached from a monazite-based waste form. The applications of these results for the development of actinide waste forms are discussed. The current data on the behavior of real actinide-doped monazite suggest that this type of ceramic waste form is not very resistant, even in relatively short time periods.

1. Introduction

The modeling of the long-term behavior of crystalline forms for actinides immobilization greatly benefits from investigation of the natural minerals comprising radioactive elements, such as uranium and thorium. Direct extrapolation of the results for the natural minerals to synthetic ceramics is complicated by markedly different concentrations of radionuclides and dose rate effects; nevertheless, some hints can be obtained. Monazite—ideally, a rare earth element phosphate (REEPO4) with monoclinic structure—occupies a special place among actinide-containing minerals with remarkable chemical stability, as natural monazite is almost never present in metamict, that is, an amorphized state. Indeed, only very few works report diffuse X-ray diffraction patterns of natural Th-rich monazite [1]. Moreover, the samples studied in [1] were likely huttonite–monazite solid solutions, which are less radiation-resistant than pure monazite. The radiation damage in natural monazite is largely confined to small domains, eventually enriched in impurities [2,3]. A low critical temperature of recovery from radiation damage in monazite [4] leads to gradual annealing of natural grains and explains the lack of metamict samples (for review, see [5]), although the exact mechanism remains unclear [6]. Together with the remarkable flexibility of the monazite crystalline lattice and its chemical stability [7], these factors led to proposals of using monazite-based materials for the long-term immobilization of excess actinides [8,9,10].
However, to be industrially attractive, real ceramic waste forms should possess a reasonably high loading capacity to actinides, and realistic dose rates are several orders of magnitude higher than what is encountered in natural minerals containing admixtures of U and Th. Studies of radiation damage induced by ion beams or due to self-irradiation from isotopes with high specific activity, such as 238Pu or 244Cm, show that, in these conditions, monazite can readily be turned amorphous [4,5,11,12,13]. Moreover, although monazite is chemically rather stable, it is nevertheless subject to dissolution and other types of alterations. Cases of significant alteration of natural monazite with dramatic losses of Th, U, and radiogenic Pb are not uncommon (e.g., [14,15]).
An unusual type of alteration of 238Pu-doped single crystal Eu-monazite and polycrystalline La-monazite is observed on timescales of just several years of storage in ambient conditions: the formation of “blisters”, or crust, which consists of Pu-containing rhabdophane (hydrated REE phosphate) with variable composition [16,17,18]. Its formation is almost certainly caused by intense damage and ionization by α-particles. The newly formed rhabdophane consists of submicron grains and is only poorly retained on the monazite surface. Whereas rhabdophane formation during a monazite-fluid interaction is a well-known phenomenon (e.g., [19]), its appearance, due to contact with atmospheric moisture, is unexpected, and certainly detrimental for eventual waste forms. Even although rhabdophane is poorly soluble [20], its nanocrystallinity may lead to the remobilization of actinides due to the destructive influence of high energy ions on nanoparticles [21]. Actinide-bearing nanoparticles are very mobile in geological environments and their migration is a subject of serious concern related to the safety of radioactive waste disposal. An investigation into the nanoparticles of LaPO4 doped with α-emitting 225Ac showed that the retention of daughter products, including radioactive isotopes, by such nanoparticles is not very high; losses up to 50% were observed on a timescale of one month, which corresponds to 3T1/2 of the parent isotope (note the existence of several energetic decays in the 225Ac chain) [22]. The aim of the latter work differed from the examination of the structural changes of the nanoparticles and no relevant microscopic observations were performed. Some of the daughters may have been lost due to recoils; some might have been damaged, or they may have even exploded the nanoparticles following the mechanism described in [23]. In any case, the paper [22] demonstrates that even on moderate a timescale, monazite-based nanoparticles are not that resistant against loss of radionuclides.
This contribution concerns the behavior of Pu-containing REE orthophosphates and contains two independent but related parts. First, we present the experimental data behavior of aged (17-year-old) 238Pu-doped Eu-monazite single crystal. This part of the paper reveals the influence of the storage medium on the development of secondary mineralization on the surfaces of the actinide-doped monazite. In the second part, we describe the atomistic calculations of the thermodynamic functions, which reveal peculiarities of structural control on the behavior of the solid solution in Pu-REEPO4 phases. The results of both experimental and modeling studies are discussed in the context of the application of monazite-based waste forms for the immobilization of actinides.

2. Materials and Methods

2.1. Samples and Experimental Details

Monazite single crystals with the composition of Eu0.937Pu0.063PO4 doped with 238Pu (4.9 wt% 238Pu + 1.1 wt% of other Pu isotopes) were grown using the flux method in December 2003 at the V.G. Khlopin Radium institute [24]. All the samples were stored at 20–25 °C in ambient conditions. We note that, due to a rather short T1/2 (87.7 y), the use of 238Pu allowed for the investigation of long-term effects from Pu decay in reasonably short periods of time, although the dose rate is much higher than in the case of, for example, 239Pu, and thus, direct comparison should be done with caution. The choice of europium as a matrix component was driven by the intention of synthesis of self-luminescing light sources and by the possibility of detecting minor changes in the matrix structure using Eu photoluminescence (PL) spectra. In addition, the crystallochemical behavior of Eu is similar to that of Am, another actinide of high importance for the immobilization of fractionated waste.
In [16,17] we presented the results of the characterization of these samples of various ages (spanning in range from 11 to 15 years) using spectroscopy, X-ray, electron diffraction and absorption, and optical and electron microscopy. The samples preserved crystallinity, as shown by single crystal XRD and Raman scattering, but developed a whitish shell made of Pu-containing rhabdophane. There is strong evidence that the extent of the rhabdophane formation was influenced by the morphology of the monazite substrate [17].
In contrast to the samples studied in [16,17], the single crystal from the current work was hermetically sealed in a glass ampoule with ambient air several months after the synthesis. After 17 years of storage at 20–25 °C, the sample was studied using optical microscopy (Olympus BX51), and Raman and photoluminescence spectroscopy using a Renishaw inVia Reflex spectrometer. All studies were performed in situ in the ampoule; thus, interferences from the glass wall were present in Raman spectra. Lasers with 405 and 633 nm and long working distance objectives of 50× and 100× magnification were employed. Laser power was kept to a reasonable minimum and no changes in the sample under the probing beam were observed either visually or spectroscopically.

2.2. Atomistic Simulation

The semiempirical method of structural simulation described in detail in [25] was used for the theoretical estimation of the incorporation energies of plutonium impurities into the monazite (space group P21/n), xenotime (I41/amd), and anhydrous rhabdophane (C2) TRPO4 structures (TR = Ln, Y). We note that, due to crystallochemical reasons, not all compounds mentioned above exist for every single rare earth element (REE), some of them are metastable or absent (e.g., [26]). Following [25], “virtual” phosphates were developed for these cases. The following effective charges were used: q (REE/Pu) = 1.6 e0, q (P) = 1.2 e0, and q (O) = −0.7 e0. We employed Morse potentials developed in [25], which provided an excellent description of monazite and xenotime crystal structures and reproduced available experimental S0(298 K) values with good accuracy:
(r) = DM∙[exp (−2α (rr0)) − 2 exp (−α (rr0))]
The three variable parameters DM (eV), α (Å−1), and r0 (Å) denote, respectively, the energy of dissociation, the parameter of “softness” of the chemical bond, and the optimal interatomic distance in a contact pair. The parameters DM, α, r0, and Rmax for light rare earth elements were optimized using the experimental structural data for the corresponding monazites from [27]. The parameters for heavy rare earth elements (Tb–Lu) were additionally optimized using the xenotime structural data from [28]. The Pu–O interaction parameters were optimized using the PuPO4 monazite structural data from [27]. All the calculations were performed using GULP 4.1 and 4.5 software [29].

3. Results

3.1. Optical Microscopy and Spectroscopy

The gas inside the ampoule possessed a reddish-brown color. The Raman spectra of the gas phase contained two strong peaks centered at 750 and 1319 cm−1, due to NO2 [30]. The formation of nitrogen dioxide was also observed by us in a sealed ampoule with 238Pu-doped zircon of comparable age. The formation of NO2 was due to the radiolysis of the sealed gases ([31] and refs. therein).
Optical photographs (Figure 1) of the studied monazite crystal were made in a hermetically sealed ampoule (see Section 2.1 for detail). The monazite crystal was rather uniformly covered by a thick white crust, which appeared to be considerably more developed in comparison with other samples from the same synthesis batch, which were stored in conditions of nonrestricted airflow [16,17]. The degree of the crust development and propensity to form bubble-like formations appeared to be determined by the morphology of the monazite crystal. The crust started cracking in some places, indicating brittle behavior. One apex of the specimen broke away, revealing a dark brown color of the monazite core; however, light absorption in the crust may have altered the color perception.
The Raman and photoluminescence spectra of the shell and of the exposed core are shown in Figure 2, together with the representative spectra from our previous studies. The Raman spectra of the bulk monazite and crust (Figure 2b, Curves 1 and 2) were consistent with the monazite–rhabdophane system. However, in the crust, the intensity of the asymmetric P–O stretch in the phosphate groups (1045 cm−1) was markedly enhanced and dominated the spectrum; the symmetric stretch was very weak. At present, we do not have an unambiguous explanation of this observation, but a strong distortion of the PO43− tetrahedra, and, possibly, protonation, are likely the causes.
In this work, the photoluminescence spectra of the crust were recorded using 405 nm excitation, which may have influenced the comparison with previous results. As expected, the PL spectra of the crust were dominated by the hypersensitive Eu 5D07F3 transition, peaking at 617 nm. The transition was rather narrow and consisted of a single peak. The 5D07F4 transition was very weak. The absence of the splitting of the 5D0 -> 7F0 peak (576 nm) indicates a single site for Eu3+, but its exact symmetry remains uncertain due to the ambiguities in the analysis of the 5D0 -> 7F1 transition. In any case, the local environment of the europium ion in the crust formed on the specimen in a sealed ampoule differed from that in other studied samples.
The principal results of the experimental study are confirmation of the rather rapid formation of secondary phases on the surfaces of the actinide-doped monazite samples observed in our previous works. This secondary phase is likely dominated by rhabdophane-like material, is mechanically unstable, as shown by the presence of cracks, and consists of submicron grains. An important difference of the current work is that it shows that the rhabdophane formation proceeds even in a closed system and is likely even accelerated in such an environment. The detrimental effect of the rhabdophane formation is emphasized by the preferential partitioning of Pu into this phase, as shown by the modeling results presented in the next section.

3.2. Thermodynamic Properties

3.2.1. Estimation of Interaction Energy and Enthalpy of Mixing

The reliability of the employed set of interatomic potentials and of the calculation approach as a whole is supported by a very good correspondence of the calculated and experimentally obtained crystal structures and the thermodynamic data of the monazites and xenotimes [25]. Using this set of potentials, the structural, elastic, and thermodynamic properties of the dehydrated rhabdophanes were calculated and summarized in Table 1.
The mixing properties of the binary systems of rare-earth elements (TR = Ln, Y) with plutonium were calculated (Table 2). The Gibbs free energy of mixing was calculated using the following equation:
∆Gmix = ∆Hmix − T∆Smix = x1x2 (x2Q1 + x1Q2) − T∆Smix
where ∆Hmix is the formation enthalpy at T = 0 K of a binary solid solution comprising molar fractions xi of the components. The interaction energies Q1 and Q2 can be calculated with the following equations:
Q1 = Edef (TR2 in TR1PO4) + Estr (TR1PO4) − Estr (TR2PO4)
Q2 = Edef (TR1 in TR2PO4) + Estr (TR2PO4) − Estr (TR1PO4)
where Edef is the energy of an isolated defect consisting of TR1 isomorphous impurity in a phosphate of TR2; Estr is the phosphate structural energy per formula unit.
Table 3 shows the enthalpy of mixing of the equimolar compositions of the solid solutions calculated using Equations (2) and (3) under the assumption that the interaction energy depends linearly on the composition of a solid solution. It is worth noting that the estimated values ∆Hmix strongly depend on ∆R.

3.2.2. Calculation of the Properties of Mixing for Solid Solutions in Systems LaPO4–PuPO4, EuPO4–PuPO4, GdPO4–PuPO4 and TbPO4–PuPO4 by the Supercell Method

The Gibbs free energies in the systems LaPO4–PuPO4, EuPO4–PuPO4, GdPO4–PuPO4, and TbPO4–PuPO4 were calculated using the large supercell approach [33]. Three types of structure—monazite, xenotime, and rhabdophane—were considered. A 4 × 4 × 4 supercell contains, depending on the structure, 192 to 256 isomorphic positions in the cationic sublattice (Figure 3). The choice of the elements is explained by the existence of our own experimental data for La and Pu [12,16,17,18]; Gd and Tb were selected as representatives of light (La–Gd) and heavy (Tb–Lu) lanthanides. The calculated Gibbs energies and mixing enthalpies are shown in Table 4.
The experimental [18,34] and theoretical [35,36] data suggest that the optimum amount of Pu in monazite-based waste forms should not exceed ~15 at.%. Based on this observation, the TRPO4–PuPO4 binary solid solutions containing from 5% to 15% plutonium were also studied. Supercells with four different atomic configurations representing a disordered solid solution as close as possible were used. The dependencies of the enthalpy of mixing on the composition of a solid solution for three polymorphic modifications of rare earth phosphates are shown in Figure 4. For all systems except TbPO4–PuPO4 with a xenotime structure, the values of mixing enthalpy obtained by the supercell method are systematically smaller for a small percentage than the values obtained in the approximation of infinite dilution [25]. Figure 5 shows the dependence of the obtained free energy values for solid solutions of Pu in REE phosphates with different structures and in a broad temperature range.

4. Discussion

The results of the present work provide important information about the potential use of monazite-based ceramics or compounds for the long-term immobilization of actinides. Whereas the investigation of natural samples and extensive laboratory studies of the broad range of non-radioactive and low-radioactive synthetic REE-monazite ceramics and crystals gave promising results, the experiments with industry-relevant levels of doping with actinides show that chemical and radiation stability of highly radioactive material is not sufficiently high.
Comparison of the present work with our previous studies suggests that rhabdophane formation on the surfaces of actinide-doped monazite crystals is stimulated not only (or solely) by the interaction with atmospheric moisture, but is markedly enhanced by radiolysis effects. It is reasonable to suggest that the intense ionization in the vicinity of a Pu-doped monazite surface, either external or freshly exposed by mechanical cracking due to swelling, may produce, not only NO2 observed by us in significant amounts, but other oxidizers, such as hydrogen peroxide or even some acidic compounds as well. We note here that the treatment of monazite with hot concentrated acids is a common approach to release REEs, Th, and U (monazite “cracking”, e.g., [37]). The investigation into the synthetic La-monazite leaching in acidic media [38] indicated the formation of a thin rhabdophane layer, which controls subsequent processes.
The marked dependency of the extent of the secondary phase formation on crystallographic orientation is not yet fully understood. It may be caused by preferential incorporation of a Pu impurity in different growth sectors (see discussion in [17]), or reflect a kind of epitaxial relationship. This issue presents certain applied interest, since the mechanical properties of La-monazite ceramics notably depend on the morphology of comprising grains [39]. Moreover, some important thermodynamic properties of monazite (e.g., thermal capacity) also demonstrated marked morphological effects [40]. Consequently, an understanding of the crystallographic control of the monazite alteration may help to reduce its impact on the waste form stability.
According to the results presented in Section 3.2, the formation of the rhabdophane phase on the monazite-based actinide-loaded waste form is not that benign. The enthalpy of mixing and the Gibbs free energy of Pu-REEPO4 solid solutions with various structures indicate that, for compounds of light REE, the plutonium tends to redistribute into the rhabdophane phase; for heavy REE, orthophosphates monazite is the main Pu host. Rhabdophane is generally less stable against chemical alteration, and thus, its propensity to concentrate actinides is clearly detrimental. Even recalling that rhabdophane may transform back into monazite upon heating, the thermodynamically-favored concentration of Pu in the former phase may lead to the formation of free Pu-phase, for example, an oxide, readily forming mobile colloids. Together with the obvious nanocrystallinity of the rhabdophane, the eventual Pu-phase may easily lead to the remobilization of an actinide from the waste form and their near- and far-field migration.

5. Conclusions

We present new experimental and computational results relevant for the understanding of the long-term behavior of monazite-based waste forms loaded with several atomic percentages of plutonium. It has been shown that conclusions about the high performance of monazite waste forms based on the studies of simulant materials, that is, those with rare earth elements only, are too optimistic. It is usually assumed that a waste form will experience degradation only in the relatively distant future, when engineering barriers and containers will (partly) fail and fluids will be in contact with the form. We have shown that the marked degradation of Pu-doped monazite ceramics and even single crystals readily occurs even during storage in air. Presumably, the radiolysis of atmospheric moisture and the surrounding air by α-particles dramatically accelerates monazite alteration. This process might be slowed down by the storage of the waste forms in inert, dry atmospheres, but such an approach will make the repository even more costly and will set new criteria for hermeticity of the canisters on very long timescales. Thermodynamic calculations suggest that, in the case of orthophosphates of light, REE plutonium favors less stable rhabdophane structures. Recalling that rhabdophane almost always consists of mechanically fragile and, importantly, micron- and submicron crystallites, the rapid formation of rhabdophane and/or other secondary phases demonstrated in our studies raises major concerns about the eventual large-scale use of monazite-based forms in real life.

Author Contributions

P.M., N.E.—thermodynamic calculations and writing; B.B.—synthesis and storage of samples; A.A.—spectroscopic measurements; A.S.—conceptualization, writing. All authors have read and agreed to the published version of the manuscript.

Funding

Experimental part of this research was funded by the Ministry of Science and Higher Education (grant agreement no. 075–15–2020–782).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Karkhanavala, M.D.; Shankar, J. An X-ray study of natural monazite: I. Proc. Indian Acad. Sci. 1954, 40, 67–71. [Google Scholar] [CrossRef]
  2. Meldrum, A.; Boatner, L.A.; Weber, W.J.; Ewing, R.C. Radiation damage in zircon and monazite. Geochim. Cosmochim. Acta 1998, 62, 2509–2520. [Google Scholar] [CrossRef]
  3. Seydoux-Guillaume, A.M.; Wirth, R.; Nasdala, L.; Gottschalk, M.; Montel, J.M.; Heinrich, W. An XRD, TEM and Raman study of experimentally annealed natural monazite. Phys. Chem. Miner. 2002, 29, 240–253. [Google Scholar] [CrossRef]
  4. Meldrum, A.; Boatner, L.A.; Ewing, R.C. Displacive irradiation effects in the monazite- and zircon-structure orthophosphates. Phys. Rev. B 1997, 56, 13805–13814. [Google Scholar] [CrossRef]
  5. Ewing, R.C.; Meldrum, A.; Wang, L.; Wang, S.X. Radiation-Induced Amorphization. Rev. Miner. Geochem. 2000, 39, 319–361. [Google Scholar] [CrossRef]
  6. Nasdala, L.; Akhmadaliev, S.; Burakov, B.E.; Chanmuang, N.C.; Škoda, R. The absence of metamictisation in natural monazite. Sci. Rep. 2020, 10, 14676. [Google Scholar] [CrossRef]
  7. Clavier, N.; Podor, R.; Dacheux, N. Crystal chemistry of the monazite structure. J. Eur. Ceram. Soc. 2011, 31, 941–976. [Google Scholar]
  8. McCarthy, G.L.; White, W.B.; Pfoertsch, D.E. Synthesis of nuclear waste monazites, ideal actinide hosts for geological disposal. Mater. Res. Bull. 1978, 13, 1239–1245. [Google Scholar] [CrossRef]
  9. Boatner, L.A.; Beall, G.W.; Abraham, M.M.; Finch, C.B.; Huray, P.C.; Rapaz, M. Monazite and other lanthanide orthophosphates as alternative actinide waste forms. In Scientific Basis for Nuclear Waste Management; Northrup, C.J.M., Jr., Ed.; Plenum Press: New York, NY, USA, 1980; Volume 2, pp. 289–296. [Google Scholar]
  10. Dacheux, N.; Clavier, N.; Podor, R. Monazite as a promising long-term radioactive waste matrix: Benefits of high-structural flexibility and chemical durability. Am. Miner. 2013, 98, 833–847. [Google Scholar] [CrossRef]
  11. Karioris, F.G.; Gowda, K.; Cartz, L. Heavy ion bombardment on monoclinic ThSiO4, ThO2, and monazite. Radiat. Eff. Lett. 1981, 58, 1–3. [Google Scholar] [CrossRef]
  12. Burakov, B.E.; Yagovkina, M.A.; Garbuzov, V.M.; Kitsay, A.A.; Zirlin, V.A. Self-irradiation of monazite ceramics: Contrasting behavior of PuPO4 and (La,Pu)PO4 doped with Pu-238. In Mater. Res. Symp. Proc. Scientific Basis for Nuclear Waste Management XXVIII; Hanchar, J.M., Stroes-Gascoyne, S., Browning, L., Eds.; Materials Research Society: Warrendale, PA, USA, 2004; Volume 824, p. CC4.1. [Google Scholar] [CrossRef]
  13. Nasdala, L.; Grötzschel, R.; Probst, S.; Bleisteiner, B. Irradiation damage in monazite (CePO4): An example to establish the limits of Raman confocality and depth resolution. Can. Mineral. 2010, 48, 351–359. [Google Scholar] [CrossRef]
  14. Poitrasson, F.; Chenery, S.; Shepherd, T.J. Electron microprobe and LA-ICP-MS study of monazite hydrothermal alteration: Implications for U-Th-Pb geochronology and nuclear ceramics. Geochim. Cosmochim. Acta 2000, 64, 3283–3297. [Google Scholar] [CrossRef]
  15. Williams, M.L.; Jercinovic, M.J.; Harlov, D.E.; Budzyń, B.; Hetherington, C.J. Resetting monazite ages during fluid-related alteration. Chem. Geol. 2011, 283, 218–225. [Google Scholar]
  16. Shiryaev, A.A.; Nickolsky, M.S.; Averin, A.A.; Grigoriev, M.S.; Zubavichus, Y.V.; Vlasova, I.E.; Petrov, V.G.; Burakov, B.E. Structural peculiarities of aged 238Pu-doped monazite. MRS Adv. 2016, 1, 4275. [Google Scholar] [CrossRef]
  17. Shiryaev, A.A.; Burakov, B.E.; Nickolsky, M.S.; Yapaskurt, V.O.; Pavlushin, A.D.; Grigoriev, M.S.; Vlasova, I.E. Surface features on aged 238Pu-doped Eu-monazite. Radiochim. Acta 2020, 108, 353–360. [Google Scholar] [CrossRef]
  18. Shiryaev, A.A.; Burakov, B.E.; Yapaskurt, V.O.; Egorov, A.V.; Vlasova, I.E. Microstructure of Aged 238Pu-doped La-monazite Ceramic and Peculiarities of its X-ray Emission Spectra. MRS Adv. 2020, 5, 1–7. [Google Scholar] [CrossRef]
  19. Du Fou de Kerdaniel, E.; Clavier, N.; Dacheux, N.; Terra, O.; Podor, R. Actinide solubility-controlling phases during the dissolution of phosphate ceramics. J. Nucl. Mater. 2007, 362, 451–458. [Google Scholar] [CrossRef]
  20. Shelyug, A.; Mesbah, A.; Szenknect, S.; Clavier, N.; Dacheux, N.; Navrotsky, A. Thermodynamics and stability of rhabdophanes, hydrated rare earth phosphates REPO4 n H2O. Front. Chem. 2018, 6, 604. [Google Scholar] [CrossRef]
  21. Lu, F.; Shen, Y.; Sun, X.; Dong, Z.; Ewing, R.C.; Lian, J. Size dependence of radiation-induced amorphization and recrystallization of synthetic nanostructured CePO4 monazite. Acta Mater. 2013, 61, 2984–2992. [Google Scholar] [CrossRef]
  22. Woodward, J.; Kennel, S.J.; Stuckey, A.; Osborne, D.; Wall, J.; Rondinone, A.J.; Standaert, R.F.; Mirzadeh, S. LaPO4 Nanoparticles Doped with Actinium-225 that Partially Sequester Daughter Radionuclides. Bioconjugate Chem. 2011, 22, 766–776. [Google Scholar] [CrossRef]
  23. Shiryaev, A.A.; Hinks, J.A.; Marks, N.A.; Greaves, G.; Valencia, F.J.; Donnelly, S.E.; González, R.I.; Kiwi, M.; Trigub, A.L.; Bringa, E.M.; et al. Ion implantation in nanodiamonds: Size effect and energy dependence. Sci. Rep. 2018, 8, 5099. [Google Scholar] [CrossRef] [PubMed]
  24. Burakov, B.E.; Ojovan, M.I.; Lee, W.E. Crystalline Materials for Actinide Immobilisation, Materials for Engineering; Imperial College Press: London, UK, 2010; Volume 1, p. 197. [Google Scholar]
  25. Eremin, N.N.; Marchenko, E.I.; Petrov, V.G.; Mitrofanov, A.A.; Ulanova, A.S. Solid solutions of monazites and xenotimes of lanthanides and plutonium: Atomistic model of crystal structures, point defects and mixing properties. Comput. Mater. Sci. 2019, 157, 43–50. [Google Scholar] [CrossRef]
  26. Boatner, L.A. Synthesis, Structure, and Properties of Monazite, Pretulite, and Xenotime. Rev. Miner. Geochem. 2002, 48, 87–121. [Google Scholar] [CrossRef]
  27. Jardin, R.; Pavel, C.C.; Raison, P.E.; Bouëxière, D.; Santa-Cruz, H.; Konings, R.J.; Popa, K. The high-temperature behaviour of PuPO4 monazite and some other related compounds. J. Nucl. Mat. 2008, 378, 167–171. [Google Scholar] [CrossRef]
  28. Ni, Y.; Hughes, J.M.; Mariano, A.M. Crystal chemistry of the monazite and xenotime structures. Amer. Miner. 2005, 80, 21–26. [Google Scholar] [CrossRef]
  29. Gale, J.D.; Rohl, A.L. The general utility lattice program (GULP). Mol. Simul. 2003, 29, 291–341. [Google Scholar] [CrossRef]
  30. Hardwick, J.L.; Brand, J.C.D. Anharmonic potential constants and the large amplitude bending vibration in nitrogen dioxide. Can. J. Phys. 1976, 54, 80–91. [Google Scholar] [CrossRef]
  31. Reed, D.; Van Konynenburg, R. Effect of Ionizing Radiation on Moist Air Systems. MRS Proc. 1987, 112, 393–404. [Google Scholar] [CrossRef] [Green Version]
  32. Mesbah, A.; Clavier, N.; Elkaim, E.; Gausse, C.; Kacem, I.B.; Szenknect, S.; Dacheux, N. Monoclinic Form of the Rhabdophane Compounds: REEPO4·0.667H2O. Cryst. Growth Des. 2014, 14, 5090–5098. [Google Scholar] [CrossRef]
  33. Eremin, N.N.; Deyanov, R.Z.; Urusov, V.S. Choice of the supercell with the optimum atomic configuration in simulation of disordered solid solutions. Glass Phys. Chem. 2008, 34, 9–18. [Google Scholar] [CrossRef]
  34. Arinicheva, Y.; Popa, K.; Scheinost, A.C.; Rossberg, A.; Dieste-Blanco, O.; Raison, P.; Cambriani, A.; Somers, J.; Bosbach, D.; Neumeier, S. Structural investigations of (La,Pu)PO4 monazite solid solutions: XRD and XAFS study. J. Nucl. Mat. 2017, 493, 404–411. [Google Scholar] [CrossRef]
  35. Kowalski, P.M.; Beridze, G.; Vinograd, V.; Bosbach, D. Heat capacities of lanthanide and actinide monazite-type ceramics. J. Nucl. Mater. 2015, 464, 147–154. [Google Scholar] [CrossRef]
  36. Ji, Y.; Kowalski, P.M.; Kegler, P.; Huittinen, N.; Marks, N.A.; Vinograd, V.L.; Arinicheva, Y.; Neumeier, S.; Bosbach, D. Rare-Earth Orthophosphates from Atomistic Simulations. Front. Chem. 2019, 7, 197. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Kumari, A.; Panda, R.; Jha, M.K.; Kumar, J.R.; Lee, J.Y. Process development to recover rare earth metals from monazite mineral: A review. Miner. Eng. 2015, 79, 102–115. [Google Scholar] [CrossRef]
  38. Arinicheva, Y.; Gausse, C.; Neumeier, S.; Brandt, F.; Rozov, K.; Szenknect, S.; Dacheux, N.; Bosbach, D.; Deissmann, G. Influence of temperature on the dissolution kinetics of synthetic LaPO4-monazite in acidic media between 50 and 130 °C. J. Nucl. Mater. 2018, 509, 488–495. [Google Scholar] [CrossRef]
  39. Neumeier, S.; Arinicheva, Y.; Clavier, N.; Podor, R.; Bukaemskiy, A.; Modolo, G.; Dacheux, N.; Bosbach, D. The effect of the synthesis route of monazite precursors on the microstructure of sintered pellets. Progr. Nucl. Energy 2016, 92, 298–305. [Google Scholar] [CrossRef]
  40. Bryukhanova, K.I.; Nikiforova, G.E.; Tyurin, A.V.; Kondrat’eva, O.N.; Gavrichev, K.S. Effect of the Particle Habit on the Heat Capacity and Thermodynamic Functions of EuPO4 in the Range 7–1600 K. Russ. J. Inorg. Chem. 2020, 65, 681–687. [Google Scholar] [CrossRef]
Figure 1. Optical photographs of the 238Pu-doped Eu-monazite crystal. (A,B): As-grown crystals and the same samples 14 months later [24], before sealing in the ampoule. Arrows in (B) show tiny particles spontaneously detached from the crystals. Scale bar is the same for (A) and (B), and is 1 mm. (CE): Crystal in a sealed ampoule. Variations in the degree of development of the secondary phases on different crystallographic faces were observed. The brightness of the photographs (CE) was increased to partly offset the brownish background due to NO2. Scale bar is the same for (CE) and is 0.5 mm.
Figure 1. Optical photographs of the 238Pu-doped Eu-monazite crystal. (A,B): As-grown crystals and the same samples 14 months later [24], before sealing in the ampoule. Arrows in (B) show tiny particles spontaneously detached from the crystals. Scale bar is the same for (A) and (B), and is 1 mm. (CE): Crystal in a sealed ampoule. Variations in the degree of development of the secondary phases on different crystallographic faces were observed. The brightness of the photographs (CE) was increased to partly offset the brownish background due to NO2. Scale bar is the same for (CE) and is 0.5 mm.
Sustainability 13 01203 g001
Figure 2. Spectra of the 238Pu-doped Eu-monazite sample in comparison with some data from [16,17]. (A) photoluminescence spectra; excitation wavelength and age of the sample is indicated for every curve. (B) Raman spectra. Curve 1: core of the crystal from the present work (17 years old); 2: shell of the present crystal; 3: bulk of 15-year-old sample; 4: shell of 15-year-old crystal; 5: bulk of 11-year-old crystal. Curves 1, and 2 were recorded using 633 nm excitation; 3–5, with 532 nm excitation.
Figure 2. Spectra of the 238Pu-doped Eu-monazite sample in comparison with some data from [16,17]. (A) photoluminescence spectra; excitation wavelength and age of the sample is indicated for every curve. (B) Raman spectra. Curve 1: core of the crystal from the present work (17 years old); 2: shell of the present crystal; 3: bulk of 15-year-old sample; 4: shell of 15-year-old crystal; 5: bulk of 11-year-old crystal. Curves 1, and 2 were recorded using 633 nm excitation; 3–5, with 532 nm excitation.
Sustainability 13 01203 g002
Figure 3. Rhabdophane 4 × 4 × 4 supercell of Eu182Pu10P192O768 composition. Left: [100]; right: [001] plane.
Figure 3. Rhabdophane 4 × 4 × 4 supercell of Eu182Pu10P192O768 composition. Left: [100]; right: [001] plane.
Sustainability 13 01203 g003
Figure 4. Enthalpy of mixing of the systems LaPO4–PuPO4, EuPO4–PuPO4, LaPO4–PuPO4, GdPO4–PuPO4, and TbPO4–PuPO4 as a function of Pu content. Symbols: calculation with the 4 × 4 × 4 supercell model for different atomic configurations; mean square deviations are shown. Solid lines: calculation of the infinite dilution approximation.
Figure 4. Enthalpy of mixing of the systems LaPO4–PuPO4, EuPO4–PuPO4, LaPO4–PuPO4, GdPO4–PuPO4, and TbPO4–PuPO4 as a function of Pu content. Symbols: calculation with the 4 × 4 × 4 supercell model for different atomic configurations; mean square deviations are shown. Solid lines: calculation of the infinite dilution approximation.
Sustainability 13 01203 g004
Figure 5. Free energies of solid solution for LaPO4–PuPO4, EuPO4–PuPO4, GdPO4–PuPO4, and TbPO4–PuPO4 systems as a function of Pu content. The calculations were performed for monazite, xenotime, and rhabdophane structures at 200, 400, and 600 K. The mean square deviations of values for different atomic configurations are shown. The lines are drawn to guide the eye.
Figure 5. Free energies of solid solution for LaPO4–PuPO4, EuPO4–PuPO4, GdPO4–PuPO4, and TbPO4–PuPO4 systems as a function of Pu content. The calculations were performed for monazite, xenotime, and rhabdophane structures at 200, 400, and 600 K. The mean square deviations of values for different atomic configurations are shown. The lines are drawn to guide the eye.
Sustainability 13 01203 g005
Table 1. Structural, elastic, and thermodynamic characteristics of dehydrated rhabdophanes.
Table 1. Structural, elastic, and thermodynamic characteristics of dehydrated rhabdophanes.
ParameterLaPO4CePO4PuPO4PrPO4NdPO4SmPO4SmPO4 (Exp., [32])EuPO4GdPO4
Volume, Å3594.05580.70573.78571.41562.57545.70540.97 (1)539.80534.56
a, Å12.645412.542812.499812.471012.416912.286212.1443 (1)12.227512.1754
b, Å7.30087.24167.21677.20017.16897.09357.0178 (1)7.05967.0295
c, Å6.43466.39336.36076.36376.31996.26146.3476 (1)6.25346.2459
β,°90.0090.0090.0090.0090.0090.0090.02 (1)90.0090.00
Structural energy, eV−34.314−34.582−34.687−34.802−35.220−35.506-−35.568−35.648
K, GPa94.5893.0139.2639.8646.1752.69-50.9850.01
G, GPa41.9443.2945.3844.0547.5050.82-49.8349.18
Q1(TRPO4-PuPO4), kJ/mol0.9640.097-0.0390.2641.767-2.6873.624
Q2(PuPO4-TRPO4), kJ/mol0.9870.101-0.0390.2631.745 2.7423.856
ΔH0,5 (TRPO4-PuPO4), kJ/mol0.24390.0248-0.00970.06580.4390-0.67860.9350
Table 2. Calculated energies of TR–Pu binary solid solutions of several considered structural types (monazite [25], xenotime, rhabdophane), calculated by the method of point defects in the limit of infinite dilution. The radii of the spheres used were 8.5 and 18.5 Å, respectively.
Table 2. Calculated energies of TR–Pu binary solid solutions of several considered structural types (monazite [25], xenotime, rhabdophane), calculated by the method of point defects in the limit of infinite dilution. The radii of the spheres used were 8.5 and 18.5 Å, respectively.
StructureParameterLaPO4CePO4PrPO4NdPO4SmPO4EuPO4GdPO4TbPO4DyPO4YPO4HoPO4ErPO4TmPO4YbPO4LuPO4
MonaziteQ11.060.10.040.462.213.284.3810.4217.5519.5524.0131.1137.7644.1751.71
Q21.060.120.050.462.153.134.1610.0816.0017.7121.3627.1832.5837.7143.66
η1.000.910.891.001.031.051.051.031.101.101.121.161.161.171.18
XenotimeQ11.470.110.060.312.514.416.5524.7432.1534.4839.7247.9055.3563.0472.13
Q21.440.110.060.312.464.196.0418.7923.9825.7629.1934.7839.8144.9050.82
η1.021.001.041.001.021.051.081.321.341.341.361.381.391.401.42
RhabdophaneQ10.9640.0970.0390.2641.7672.6873.62411.75316.65117.89721.50226.93332.14937.36143.553
Q20.9870.1010.0390.2631.7452.7423.85612.65516.28817.60920.29924.83329.09233.34338.344
η0.980.961.001.001.010.980.940.931.021.021.061.081.111.121.14
Table 3. The enthalpy of the mixing of equimolar compositions of binary solid solutions of rare earth and plutonium phosphates. Virtual minerals are marked with an asterisk.
Table 3. The enthalpy of the mixing of equimolar compositions of binary solid solutions of rare earth and plutonium phosphates. Virtual minerals are marked with an asterisk.
Monazite Solid SolutionΔH0,5, kJ/molXenotime Solid SolutionΔH0,5, kJ/molRhabdophane Solid SolutionΔH0,5, kJ/mol
LaPO4–PuPO40.26LaPO4 *–PuPO4 *0.36LaPO4–PuPO4 *0.24
CePO4–PuPO40.03CePO4 *–PuPO4 *0.03CePO4–PuPO4 *0.02
PuPO4–PrPO40.01PuPO4 *–PrPO4 *0.01PuPO4 *–PrPO40.01
PuPO4–NdPO40.12PuPO4 *–NdPO4 *0.08PuPO4 *–NdPO40.07
PuPO4–SmPO40.54PuPO4 *–SmPO4 *0.62PuPO4 *–SmPO40.44
PuPO4–EuPO40.80PuPO4 *–EuPO4 *1.08PuPO4 *–EuPO40.68
PuPO4–GdPO41.07PuPO4 *–GdPO4 *1.57PuPO4 *–GdPO40.93
PuPO4–TbPO4 *2.56PuPO4 *–TbPO45.44PuPO4 *–TbPO4 *3.05
PuPO4–DyPO4 *4.19PuPO4 *–DyPO47.02PuPO4 *–DyPO4 *4.12
PuPO4–YPO4 *4.66PuPO4 *–YPO47.53PuPO4 *–YPO4 *4.44
PuPO4–HoPO4 *5.67PuPO4 *–HoPO48.61PuPO4 *–HoPO4 *5.23
PuPO4–ErPO4 *7.29PuPO4 *–ErPO410.34PuPO4 *–ErPO4 *6.47
PuPO4–TmPO4 *8.79PuPO4 *–TmPO411.90PuPO4 *–TmPO4 *7.66
PuPO4–YbPO4 *10.23PuPO4 *–YbPO413.49PuPO4 *–YbPO4 *8.84
PuPO4–LuPO4 *11.92PuPO4 *–LuPO415.37PuPO4 *–LuPO4 *10.24
Table 4. Mixing enthalpies and Gibbs free energies in systems LaPO4–PuPO4, EuPO4–PuPO4, GdPO4–PuPO4, and TbPO4–PuPO4 at a Pu fraction of 10%.
Table 4. Mixing enthalpies and Gibbs free energies in systems LaPO4–PuPO4, EuPO4–PuPO4, GdPO4–PuPO4, and TbPO4–PuPO4 at a Pu fraction of 10%.
ΔH0,1
kJ/mol
ΔG0,1 (200 K)
kJ/mol
ΔG0,1 (400 K)
kJ/mol
ΔG0,1 (600 K)
kJ/mol
LaPO4–PuPO4 system
Monazite0.085−0.4588−1.0027−1.5473
Xenotime0.119−0.4271−0.9743−1.5397
Rhabdophane0.067−0.4722−1.0113−1.5505
EuPO4–PuPO4 system
Monazite0.268−0.2579−0.7845−1.3117
Xenotime0.367−0.1820−0.7335−1.2850
Rhabdophane0.194−0.3457−0.8862−1.4269
GdPO4–PuPO4 system
Monazite0.376−0.1529−0.6819−1.2111
Xenotime0.533−0.0235−0.5833−1.5726
Rhabdophane0.269−0.2764−0.8224−1.3684
TbPO4–PuPO4 system
Monazite0.8670.3088−0.2351−0.7787
Xenotime2.0681.44010.8101−1.7786
Rhabdophane0.8970.82570.2557−0.3034
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mikhailova, P.; Burakov, B.; Eremin, N.; Averin, A.; Shiryaev, A. Plutonium-Doped Monazite and Other Orthophosphates—Thermodynamics and Experimental Data on Long-Term Behavior. Sustainability 2021, 13, 1203. https://doi.org/10.3390/su13031203

AMA Style

Mikhailova P, Burakov B, Eremin N, Averin A, Shiryaev A. Plutonium-Doped Monazite and Other Orthophosphates—Thermodynamics and Experimental Data on Long-Term Behavior. Sustainability. 2021; 13(3):1203. https://doi.org/10.3390/su13031203

Chicago/Turabian Style

Mikhailova, Polina, Boris Burakov, Nikolai Eremin, Alexei Averin, and Andrey Shiryaev. 2021. "Plutonium-Doped Monazite and Other Orthophosphates—Thermodynamics and Experimental Data on Long-Term Behavior" Sustainability 13, no. 3: 1203. https://doi.org/10.3390/su13031203

APA Style

Mikhailova, P., Burakov, B., Eremin, N., Averin, A., & Shiryaev, A. (2021). Plutonium-Doped Monazite and Other Orthophosphates—Thermodynamics and Experimental Data on Long-Term Behavior. Sustainability, 13(3), 1203. https://doi.org/10.3390/su13031203

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop