Library of UV-Visible Absorption Spectra of Rare Earth Orthophosphates, LnPO4 (Ln = La-Lu, except Pm)

In recent times, rare earth orthophosphates ( L n PO 4 ) have shown great potential as efficient optical materials. They possess either m o n a z i t e or x e n o t i m e –type structures. These light or heavy rare earth bearing orthophosphates also exhibit an extraordinary stability over geological time scale in nature, ∼10 9 years. In the present contribution, we measure, collect, and present a library of absorption spectra of all the L n PO 4 hosts ( L n = La–Lu, except Pm) using their single crystal samples, to conclude that the observed spectral features for wavelengths longer than 200 nm were attributable to either Ln- or defect related centers, which corroborate the fact that they have a bandgap higher than 8.0 eV. The absorption band around wavelength, 275 nm, corresponds to defect absorption related to PO 3 centers and/or oxygen vacancies. The hosts can potentially be used to study and interpret unperturbed rare earth emissions due to absence of host related absorption above 300 nm. The information presented herein is expected to serve as a library of absorption spectra for geologists, physicists, material scientists, and chemists working in the field of rare earths.


Introduction
The lanthanide (Ln) or rare earth element (REE) are f-block elements in the periodic table which can be subdivided into two categories, namely light Ln (La-Gd) and heavy Ln (Tb-Lu). Ln concentration in nature between 10-300 ppm can be considered to be a primary source for their extraction, and almost all the Ln ions are extracted from monazite, xenotime, and bastnäsites [1,2]. Light Ln being the most abundant, are found at concentrations among the highest compared to other commodities in many primary deposits. Conversely, heavy Ln are much less abundant, similar to rare metals (Sn, W, Ta) and are therefore much rarer in natural deposits. Ln are frequently incorporated into the crystalline network of carbonates, oxides, silicates, or phosphates together with other more common chemical elements but the mineral industry can only valorize a limited number of Ln minerals: bastnäsite (fluorocarbonate), lamonazite and xenotime (phosphates), and loparite (oxide). Ln have also been mined from apatite in Russia and South Africa, but this production has been extremely limited. Except for the xenotime which may be enriched in heavy Ln, all these minerals are essentially characterized by an enrichment in light Ln. Bastnäsite is the main Ln ore and is mainly extracted from Chinese mines as well as from Mountain Pass in the United States until recently. The extraction of Ln from bastnäsite is primarily carried out by acid dissolution followed by solvent extraction and ion exchange methods. Another major source of Ln ions is monazite with nearly similar concentration of lighter Ln ions but with variable heavy Ln concentrations. Monazite occurs in a wide span of geological environment and known deposits are found worldwide [2]. Since the 1970s, ionic clays and lateritic clays, have also become an important source of heavy Ln [3]. This production, only located in the south of China is allowed by a low labor cost, and relatively simple extraction processes such as by in situ leaching with neutral or acid solutions. Despite the great variety of deposits, only 5 types are exploited: carbonatites (48% of world production), deposits of alkaline magmatism (2%), ion clay deposits (36%) laterite deposits (12% ) and placers (2%). In 2018, world production was estimated at 170,000 t LnO. Apart from these primary sources, minerals such as apatite, calcium fluorophosphates, silico-phosphates, vanado-phosphates, and arseno-phosphates also contain Ln ions in very low concentration and can be used to extract them but large scale extraction procedures are still in development.
One of the important rare earth sources, Ln orthophosphates (LnPO 4 ), are ceramic materials with interest in various application domains such as for luminescence-based temperature sensing [4], solid state lighting devices [5,6], geochronology [7], geothermometry [8,9], sensors and proton conductors [10], scintillators for X-ray and γ-ray detection-based medical imaging [11], catalysts [12], thermoluminescent phosphors [13], and laser hosts [14]. They also possess high durability and radiation damage resistance properties and thus are considered in nuclear waste management as potential immobilization matrices for minor actinides and plutonium [15,16]. Their natural analogues can contain significant amounts of actinides up to 50 wt. % of Th and U without harming their crystalline structure over geological time scales [17][18][19][20]. The LnPO 4 can be prepared through various synthesis routes which includes, but not limited to, the precipitation reaction (homogeneous mixing of solutions), high temperature solid state reaction, and solvo-thermal reactions, due to simplicity and better controllability as well as rapidness of the product formation [2]. However, to use them for practical application such as to validate them as model systems, one needs to investigate their single crystals. The advantage with single crystal samples is the absence of defects associated with grain boundaries and absence of entropic effects which otherwise favour imperfection in the microstructure of the monocrystals, such as impurities, inhomogeneous strain and crystallographic defects such as dislocations [21].
In this context, the aim of the present study is to collect and present results on the UV-Visible absorption spectra of LnPO 4 single crystal samples. The samples are obtained from the Smithsonian standard library, US. The samples are generally used for electron microprobe analysis (EMPA). The Smithsonian standard library contains absorption spectra of these samples at wavelengths >350 nm. Hence, the present results on one hand will provide additional information for the Smithsonian library, and on the other hand, will act as a library of reference spectra for the end users. We study the UV-Visible absorption spectra in the wavelength region, 190-920 nm, for these single crystal samples to overlap, and extend the existing database of absorption spectra for all the LnPO 4 samples.

Single Crystal Growth
The single crystal samples of LnPO 4 (Ln = La-Lu, expect Pm) were prepared by flux technique at Oak Ridge National Laboratory, USA [22], in accordance with the synthesis procedure discussed in Ref. [23]. In a typical procedure, Ln-oxides are placed in tightly covered platinum crucibles containing lead (II) hydrogen phosphate, PbHPO 4 , to form the flux matrix at the high temperature. This is followed by an annealing step at 1633 K for 16 h. The crucibles were quenched suddenly to room temperature after being cooled to 1173 K at a rate of 1 K/h. The crystals are later dissolved in HNO 3 to extract respective single crystal samples [22]. The samples are also known to possess some Pb impurity as discussed in Ref. [24], based on the the quantitative analyses, using Kα X-ray line for Sc, Lα lines for Y and the other REE elements, and the Mα line for Pb. However, the homogeneity of the samples has been tested using the homogeneity index and the samples have been found to be consistent in rare earth content [25].

Instrumentation
The high resolution two-dimensional (2D) images of the single crystals were obtained on a confocal laser microscope model LEXT OLS4000 from (Olympus, tokyo, Japan). The samples were illuminated through the objective and the same way back to the detector (bright field). Image mapping technique was used to combine several images laterally. The confocal microscope allow a sharp image over a wide z-range, which is not possible with a normal microscope. The room temperature absorption spectra were measured in the wavelength range, 190-920 nm, using a single beam photometer with diode array detector model number SPECORD S-600 UV/Vis spectrophotometer, at a resolution of 1 nm, procured from Analytik Jena AG, Germany. The unpolished samples were measured in the transmission mode using either 1.0 mm, 0.50 mm, or 0.25 mm sample holders prepared in-house. The lamp spectrum was calibrated in each case separately before measurements of the single crystal samples. The absorption data was converted to the absorption coefficient (α) by considering width of the samples obtained from the laser profile of the 2D images.

Surface Investigation
A confocal scanning laser microscope is used to obtain 2D images of different LnPO 4 (Ln = La-Lu, except Pm) crystal surfaces using white light source for surface illumination, Figure 1. As the Ln ions are divided into light and heavy Ln for order number 1-7 and 8-14; LnPO 4 can be categorized as light Ln bearing orthophosphates (La-Gd)PO 4 , and heavy Ln bearing orthophosphates (Tb-Lu)PO 4 . The profile from all the sample surfaces evince that the crystal growth is uniform, and the crystals possess rigid top surfaces, c f . Figure 1. The sample-type and corresponding dimensions for each of the sample are summarized in Table 1 (column 1 and 2). The surface for all the samples was hard and most of the crystals were rectangular in shape suggesting that the growth was favored in a particular direction under the influence of high temperature synthesis conditions.

Structural Evaluation
The variation of lattice parameters (a, b, c) and volume (V) of a unit cell for all the LnPO 4 samples, where the atomic number increases from La → Lu, are combined in Table 1 (column 3-6). The information is collected from Refs. [26][27][28][29][30][31][32]. It is to be remembered here that we have used same sample set in the present work.
With a general electronic configuration of [Xe]4f n 6s 2 (n = 1-14) for Ln ions, the electron enters the 4 f orbital, and except for La or Lu, all the elements have partially filled f -orbitals. Also, they are shielded from outside by s-orbitals leading to similar chemical properties for almost all the Ln ions. However, their nuclear charge is poorly shielded from the outermost electron leading to lanthanide contraction effect, i.e., there is a decrease in their atomic and ionic radii with an increasing atomic number. The decrease in ionic radii for different light and heavy Ln bearing orthophosphates are compared in Figure 2 using red colour and blue colour, respectively. The Ln bearing orthophosphate samples possess axial ratios around 1.125 (for a/b) and 1.050 (for a/c) for light Ln bearing samples. Similarly, for heavy Ln bearing orthophosphate samples the values are consistent at 1.000 and 0.875 for a/b and a/c, respectively.
The structural information for all the samples compiled in Table 1 (columns 7 and 8), show that there exists two different structures for light and heavy Ln bearing orthophosphate samples. The light Ln bearing samples possess P2 1 /n space group with monoclinic monazite-structure (a = b = c, On the other hand, the space group and structure for heavy Ln bearing orthophosphate samples is P4 1 /amd and tetragonal xenotime (a = b = c, α = β = γ ∼ 90 • ), respectively. The monazite orthophosphates possesses high thermal stability due to their compact structure which is made up of nine-coordinated LnO 9 polyhedra of Ln ions, and PO 4 tetrahedra, as shown in Figure 3(top). The monazite and xenotime orthophosphates are differentiated by crystallochemical partitioning, where the latter is formed by eight coordinated LnO 8 polyhedra of Ln ion and PO 4 tetrahedra, c f . Figure 3(bottom). Due to the greater ionic radius of La as compared to Lu, and a decrease in the ionic radii with an increasing atomic number for monazite orthophosphates, one corner of adjacent PO 4 tetrahedra rotate towards the Ln ion transforming LnO 8 to LnO 9 polyhedra.

Absorption Spectra
The valence band of orthophosphates is mainly formed by 2s, 2p states of oxygen (O) atoms with an admixture of phosphorous (P) 3s, 3p orbitals, mirroring covalent bonding effects of P and O atoms. The room temperature absorption spectra for LaPO 4 and GdPO 4 are shown in Figure 4A-a.
The spectrum for LaPO 4 shows absorption around 205 and 274 nm. The absorption peak around 205 nm can be assigned to the PO 3− 4 anions in accordance with an earlier report [35], while the absorption peak around 275 nm can be assigned to the defect related absorption in the host lattice. The nature of defect can be assigned to the charge transfer from an oxygen belonging to the PO 3 center and/or an oxygen vacancy [35,36].
The absorption spectrum for GdPO 4 overlaps exactly with the one for LaPO 4 , c f . Figure 4A-a. The ionic radii for La 3+ and Gd 3+ in 9-fold coordination are 1.216 and 1.107 Å, respectively [33]. For a change of ∼9% in the ionic radius of the two Ln ions, no change in the absorption features could be noticed suggesting that the origin of absorption bands in the wavelength region, 190-300 nm, is same for the two cases. The absorption spectrum for YPO 4 :Gd 3+ and GdPO 4 have also been reported in Ref. [37]. One can expect narrow absorption lines between wavelength, 170-280 nm, due to Gd 3+ 4f-4f transitions. We also noticed similar transitions at 275 and 308 nm, which corresponds to 8 S 7/2 → 6 I J and 6 P 7/2 , respectively. The same have been marked in Figure 4B, where the energy level diagram for LREEs (in vacuum) is shown for better comparison and understanding of the absorption transitions. One also expects an absorption at 201 nm due to the absorption from the ground state 8 S 7/2 to an emitting level 4 G 7/2 . However, such transitions could not be resolved in our results probably due to low resolution of the instrument. The results for the absorption of the host lattice are presented in Refs. [38,39] based on the vacuum ultra-violet (VUV) excitation spectra. The peak host sensitization bands of LaPO 4 and GdPO 4 are very close to each other at 159 and 160 nm, suggesting that they both have almost similar band gap [39].

CePO 4
The UV-visible absorption spectrum for CePO 4 shows large absorption between wavelength, 190-500 nm, as shown in Figure 4A-b. With a melting point of ∼2000 • C, CePO 4 is one of the most stable orthophosphate which can exist in two different forms, low temperature synthesized CePO 4 · xH 2 O with a hexagonal structure [40], and high temperature synthesized CePO 4 with a monazite structure [41]. As reported in Ref. [41,42], the absorption bands for high temperature synthesized sample (similar to our sample) are observed at 273, 257, 240, 215, and 202 nm due to 4f → 5d (1−5) transitions, typical for Ce 3+ containing/doped materials. For an unperturbed Ce 3+ ion (in vacuum), the 5d 1 level is expected around 52,000 cm −1 above the lowest lying Ce 3+ 4f levels. However, due to its presence in phosphate host, the red-shift (D) of the lowest 5d level of Ce 3+ , i.e., 5d 1 is observed at ∼12,700 cm −1 above the ground state, c f . Figure 4B. The red-shift is calculated using relation, D = Energy of free ion − (1240/5d 1 energy). The spectrum for our samples also shows large absorption >300 nm. In another report [43], the absorption spectrum for rod-like microstructures synthesized through hydrothermal method showed similar absorption features as is observed for our case. The authors assigned such a shift to a change in the bandgap. However, a shift of band gap from ∼8.0 eV from Ref. [41] to ∼3.2 eV from Ref. [43] does not appear realistic for such a large absorption above wavelength, 300 nm (marked "unknown" in Figure 4A-b).

PrPO 4
The absorption spectrum for PrPO 4 single crystal sample is shown in Figure 4A-c. PrPO 4 shows a strong absorption in the visible region as evinced by the green body colour of the sample. The results for a single crystal sample and PrPO 4 nanorods [44], are compared in our recent publication [6]. The absorption spectrum shows two set of transitions: (a) Pr 3+ 4f-4f transitions in the wavelength region, 440-500 nm (Pr 3+ 3 H 4 → 3 P 2 transition), and 580-610 nm (Pr 3+ 3 H 4 → 1 D 2 transition) as marked in Figure 4B for Pr 3+ in vacuum; and (b) two broad bands at 230 (O 2− → Pr 3+ charge-transfer band) and 330 nm (Pr 3+ 4f-5d band). The position of 4f-4f transitions in single crystal and nanorod samples are observed in the same wavelength region. The reader is referred to Ref. [6] for further information.

NdPO 4
The result of the room temperature absorption spectrum for NdPO 4 sample is shown in Figure 4A-d. The sample possesses a purple-red body colour, c f . Figure 1. It absorbs strongly in the green and yellow region of the electromagnetic spectrum. The sharp peaks are observed at 355, 455, 525, 580, 680, 748, and 865 nm matching closely with the absorption bands for Nd 3+ ion in vacuum as shown in Figure 4B and also in Ref. [45]. These sharp absorption bands belong to the Nd 3+ transitions in the 4f electronic shell from ground level 4 I 9/2 to higher 4f-levels, c f . Figure 4B. As with Pr 3+ , a broad band overlaps Nd 3+ 4f-4f transitions around 365 nm. The lowest of Pr 3+ 4f state lies within bandgap in LaPO 4 , while the lowest of Pr 3+ is within conduction band when compared in the vacuum referred binding diagram (VRBE) for the two cases [39]. Furthermore, the position of O 2− → Nd 3+ charge transfer band around 250 nm matches well with the expected position in the present host as per VRBE diagram.

SmPO 4
The room temperature absorption spectra for SmPO 4 sample is shown in Figure 4A-e. SmPO 4 absorbs weakly in the blue and blue-green region. The spectrum exhibit two different transition types: (a) broad band for wavelength <300 nm; and (b) sharp features between 300-550 nm (Sm 3+ 4f-4f 6 H 5/2 → higher multiplet transitions as marked in the energy level diagram, Figure 4B). The UV-Vis absorption spectrum for SmPO 4 single crystal and nano powders from Ref. [46], are compared in our recent work [4]. In another report [47], no absorption band was observed for SmPO 4 around 300 nm due to weak interactions between Sm and P ions. However, based on our investigations in Ref. [4], we infer that it is a culmination of multiple bands with the O 2− → Sm 3+ charge transfer absorption band around ∼200 nm, and at least one other defect related absorption band around ∼250 nm due to charge-transfer from an oxygen belonging to a PO − 3 center, and/or an oxygen vacancy [35,39,48]. The complexity in identification of absorption features of Sm 3+ arises as the free-ion configuration of Sm 3+ has 198 spin angular momentum (S), orbital angular momentum (L), and total angular momentum (J) states involving sextet, quartet, and doublet states. When such states experience crystal field, for example in C1 symmetry, each SLJ state splits into J ± 1/2 crystal sublevels to produce 2002 crystal states. To add to it, not all of these may participate in the absorption spectrum. The reader is referred to Ref. [4] for assignment of the absorption bands.

EuPO 4
The absorption spectrum for EuPO 4 is shown in Figure 4A-f. The samples possesses two broad bands at 223 and 278 nm, along with other sharp lines in the near-UV and blue region. In the published work of Ref. [35], the peak absorption for PO 3− 4 is observed at ∼160-200 nm, and the absorption longer than 200 nm (235 and 267 nm), could only be attributed to the rare earth/defect absorption. For our case, the sharp absorption lines at 316, 396, and 522 nm as also marked in the energy level diagram in Figure 4B, can be attributed to the intraconfigurational Eu 3+ ff transitions, while the other bands at 223 and 278 nm could be assigned to the O 2− → Eu 3+ charge-transfer and defect absorption, respectively. The reader is referred to Ref. [5] for further discussion on the absorption bands. The absorption spectrum at room temperature for TbPO 4 sample is shown in Figure 5A-a. The sample possess an off-white body colour. The absorption features between wavelength, 190-240 nm and 250-300 nm, can be assigned to the different spin-allowed and spin-forbidden 4f 8 → 4f 7 5d transitions of Tb 3+ , respectively [6,49,50]. The other transitions in the region, 300-500 nm, are assigned to Tb 3+ 4f-4f transitions ( 7 F 6 → higher multiplets). For example, the transitions at 352 and 490 could be marked for Tb 3+ 7 F 6 → 5 D 1 and 5 D 4 respectively, c f . Figure 5B. The absorption features of TbPO 4 in the VUV wavelength region, 140-160 nm, has been discussed and associated with the PO 4 3− anions in the recent publications [35,51], similar to the one observed for LaPO 4 and GdPO 4 .
The reader is referred to Ref. [6] for a complete assignment of absorption bands. The absorption spectrum for DyPO 4 sample is shown in Figure 5A-b. The sample possesses an ivory body color probably due to the lack of strong absorption in the visible region. The sharp line features are ascribed to the excitation from Dy 3+ 6 H 15/2 → higher electronic configuration (see Figure 5B), whereas the O 2− → Dy 3+ charge-transfer band is expected around 170 nm on theoretical calculations using an expression by Jørgensen [52]. The range is beyond the detection range of our instrument and could not be observed in the absorption spectrum. The prominent transitions are marked in Figure 5B for free Dy 3+ ion (in vacuum), for better understanding and visualization of absorption transitions. The reader is also referred to Ref. [6] for further information.

HoPO 4
The absorption spectrum for HoPO 4 sample is shown in Figure 5A-c. The sample colour for HoPO 4 is beige. The sample possesses numerous lines in the visible and UV-region. The absorption spectrum shows lines typical of Ho 3+ which can be assigned to intraconfigurational 4f-4f absorption lines typically originating from the lowest level, 5 I 8 . More interestingly, no broad band is observed in the wavelength region, 250-300 nm, similar to the one observed for light rare earth bearing orthophosphates. This can be due to two reasons: (a) the light and heavy rare earth bearing samples possess different crystal structures leading to different crystal field effects; and (b) there is a increase in the bandgap among light and heavy rare earth bearing orthophosphates, for example, a change of 0.55 eV from LaPO 4 to LuPO 4 [39], which is expected to influence the position of 4f-5d bands in the absorption spectra. A similar observation is also made for other heavy rare earth bearing orthophosphates (please see next sections). The prominent Ho 3+ transitions at 289,485, and 656 nm which belong to 5 I 8 → 3 G 3 , 5 F 3 and 5 F 5 , can be compared in Figure 5A-c and Figure 5B, respectively. Also, the complete assignment of the absorption bands is published in Ref. [53] for Ho 3+ doped YPO 4 . The sample also shows excellent photochromic properties as discussed in Ref. [54].

ErPO 4
The absorption spectrum for ErPO 4 sample is shown in Figure 5A-d along with energy level diagram for free Er 3+ ion (in vacuum) in Figure 5B. The absorption spectrum can be attributed to internal 4f-4f transitions of the configuration 4f 11 of the Er 3+ ion, from the ground state 4 I 15/2 to higher excited states. For example, the transitions at 380, 665 and 800 nm could be assigned to 4 I 15/2 → 4 G 11/2 , 4 F 9/2 , and 4 I 9/2 transitions respectively (see Figure 5B). The ErPO 4 sample possesses a pink colour probably due to strong absorption in the green region. By fitting the absorption spectral data using Judd-Ofelt model [36,55], Saltmarsh et al. predicted various radiative spectral properties such as radiative transition probability, fluorescence branching ratio and radiative decay time in Er 3+ doped LaPO 4 [56]. The spectrum for our sample (ErPO 4 ) is found to be similar to the ones reported therein and can be compared to the one presented in Figure 5B.
The absorption spectrum for TmPO 4 sample is shown in Figure 5A-e. The absorption spectrum shows sharp bands at 264, 361, 467, 693, and 792 nm. The bands can be attributed to the 3 H 6 → 1 P 2 , 1 D 2 , 1 G 4 , 7 F (2,3) , and 3 H 4 electronic transitions, as marked and compared for prominent cases in Figure 5B. The transitions are similar to one reported for Tm 3+ doped LaPO 4 [57]. Tm 3+ doped samples have found application in lasers, telecommunication and solar cells and is also known to produce upconversion luminescence in the visible region similar to Dy 3+ and Er 3+ ions. The energy level diagram and other details for the upconversion emission in Tm 3+ doped materials were reported in Ref. [58].

YbPO 4 and LuPO 4
The absorption spectra for YbPO 4 and LuPO 4 samples are shown in Figure 5A-f along with energy level diagram for free Yb 3+ ion (in vacuum) in Figure 5B. Both the samples are transparent in visible region with no absorption in the wavelength region, 250-750 nm. The spectrum is similar to the one observed for LaPO 4 and GdPO 4 . For Yb 3+ ions, one can expect absorption band around 975 nm (zero-phonon line) as the ion has only one excited state, 2 F 5/2 , located at about 10 4 cm −1 above the ground state 2 F 7/2 through the Stark splitting of the sub-levels due to the crystal field in a crystalline environment of Yb ions [59]. However, no such transition was observed in our case due to instrumental limitation, and hence marked with a symbol X in upward arrow, c f . Figure 5B. For LuPO 4 , the absorption spectrum shows a dip in the spectrum around 200 nm, which can be assigned to the PO 3− 4 anions in accordance with Ref. [35].

Conclusions
In the present work, we present a comprehensive library of UV-Vis absorption transitions for technologically important rare earth orthophosphate samples. The host absorption for all the phosphates is expected to be in the VUV-region which lies outside the detection range of our experimental setup. However, the information for the bandgaps of LaPO 4 , GdPO 4 , LuPO 4 and YPO 4 can be obtained from the literature. The host exciton creation energy (E ex ) for these hosts has been calculated in the literature using thermoluminescence investigations suggesting that the values are 8.00 eV (LaPO 4 ), 8.05 eV (GdPO 4 ), 8.55 eV (YPO 4 ), 8.60 eV (LuPO 4 ), respectively. Similar values for E ex can be obtained for other orthophosphate, CePO 4 and Eu/TbPO 4 . By combining this information, we can infer following major observations from the present results: 1. All the orthophosphates can be categorized with either monazite or xenotime structure for light and heavy bearing LnPO 4 hosts, respectively. 2. All the orthophosphates have a bandgap greater than 8.00 eV, with a change of ∼0.60 eV for LuPO 4 when compared to LaPO 4 . 3. The PO 3− 4 does not absorb wavelengths >200 nm with a peak absorption ∼180 nm. The absorption peak around 275 nm as observed mostly for light rare earth bearing orthophosphates, can be assigned to the defect absorption in the host lattice (charge transfer from an oxygen belonging to PO 3 center and/or oxygen vacancy). 4. The orthophosphate system can be considered to be an ideal system to study and compare absorption spectra of rare earth ions as the samples possess no absorption features >300 nm due to any other factors except that of Ln ions, c f . Figures 4A and 5A. 5. Some of the orthophosphates were already explored for numerous applications and the host possesses a fantastic opportunity for researchers working in the field. 6. For the LnPO 4 hosts, where the oxygen atoms are covalently bonded within the phosphate group, the luminescence efficiency is expected to be much lower. This is due to inherently higher saturation of charge within the PO 4 groups when compared to a single oxygen atom whose 2p shells remain unfilled in an ionic matrix. Funding: This research was funded by EIT RawMaterials grant number 16304.

Acknowledgments:
The authors thank the Smithsonian library, US, for providing samples with catalog # NMNH 168484 to NMNH168499 (16 samples), and Margret Fuchs for providing them to the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.