Bryce G.
Mullens
a,
Frederick P.
Marlton
*b,
Matilde
Saura-Múzquiz
c,
Michelle
Everett
d,
Cheng
Li
d,
Alicia M.
Manjon-Sanz
d,
Matthew G.
Tucker
d,
Frederic
Poineau
e,
James
Louis-Jean
e,
Supratik
Mukherjee
f,
Subrata
Mondal
f,
Ganapathy
Vaitheeswaran
*g and
Brendan J.
Kennedy
a
aSchool of Chemistry, The University of Sydney, Sydney, New South Wales 2006, Australia. E-mail: brendan.kennedy@sydney.edu.au
bCentre for Clean Energy Technology, School of Mathematical and Physical Sciences, Faculty of Science, University of Technology Sydney, Sydney, New South Wales 2007, Australia. E-mail: frederick.marlton@uts.edu.au
cDepartamento de Física de Materiales, Facultad de Ciencias Físicas, Universidad Complutense de Madrid, 28040, Madrid, Spain
dNeutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
eDepartment of Chemistry, University of Nevada Las Vegas, Las Vegas, Nevada 89154, USA
fAdvanced Centre of Research in High Energy Materials (ACRHEM), University of Hyderabad, Prof. C. R. Rao Road, Gachibowli, 500046 Hyderabad, Telangana, India
gSchool of Physics, University of Hyderabad, Prof. C. R. Rao Road, Gachibowli, 500046 Hyderabad, Telangana, India. E-mail: vaithee@uohyd.ac.in
First published on 25th November 2024
A detailed variable temperature neutron total scattering study of the potential nuclear waste matrix TlTcO4 was conducted. The long-range average structure of TlTcO4 undergoes an orthorhombic Pnma to tetragonal I41/amd phase transition below 600 K, consistent with previous synchrotron X-ray diffraction studies. However, several anomalies were observed in the Rietveld refinements to the neutron powder diffraction data, such as large atomic displacement parameters at low temperature and a shortening of the Tc–O bond distance upon heating. Modelling the short-range local structure of both the low- and high-temperature data required a lowering of symmetry to the monoclinic P21/c model due to the stereochemical activity of the Tl+ 6s2 lone pairs. Density functional theory calculations also verified this model to have a lower ground state energy than the corresponding long-range average structure. It is concluded that at low temperatures, the Tl+ 6s2 lone pairs are ‘frozen’ into the structure. Upon heating, the rigid TcO4 tetrahedra begin to rotate, as governed by the Γ3+ and M4+ modes. However, there is a disconnect between the two length scales, with the 6s2 lone pair electrons remaining stereochemically active on the local scale, as observed in the neutron pair distribution function fits. The orthorhombic Pnma to tetragonal I41/amd phase transition is seemingly the result of a change in the correlation length of the Tl+ 6s2 lone pairs, leading to a larger unit cell volume due to their uncorrelated displacements.
In developing technologies for the safe storage of Tc-containing waste, it is necessary to consider the chemistry of the commonly encountered oxidation states; namely Tc4+ and Tc7+. The Tc4+ oxidation state is only soluble in acidic conditions (below pH ∼3), meaning it often precipitates out as solid TcO2·nH2O under environmental conditions.2,3 Tc7+ forms the water-soluble and highly stable pertechnetate anion TcO4− which is much more mobile and can lead to the corrosion of Fe- and Zr-based storage containers.4–6 As a consequence, several materials have been investigated for storing Tc-containing waste to prevent its leaching into the environment. These storage mechanisms include metal oxides,7 glasses,8,9 and metal–organic frameworks (MOFs).10,11
Since Tc has no stable isotopes, it is common for studies to use Re as a nonradioactive surrogate.12–16 This is due to the predicted similarities in the chemical behaviour of Tc7+ and Re7+, which are attributed to their identical oxidation states and similar ionic radii (Tc7+ = 0.37 vs. Re7+ = 0.38 Å in tetrahedra geometry).17 However, a growing body of literature shows that Tc7+ and Re7+ do not always act identically. Lukens et al. reported that Tc7+ is significantly more volatile than Re7+, leading to a decrease in the amount stabilised in waste glasses.18 Gan et al. described the different behaviours of Tc and Re under various redox conditions, reflecting the lower standard reduction potential of Re7+ (ReO4−/ReO2 = 0.510 V) compared to Tc7+ (TcO4−/TcO2 = 0.738 V).19 These subtle differences are also evident in other oxides. For example, the perovskites CaTcO3 and SrTcO3 have exceptional magnetic properties, but their AReO3 analogues have yet to be synthesised.20–22 This suggests that differences in the electronic structures between the two cations may be important. Re has more radially-extended 5d valence orbitals than the 4d valence orbitals in Tc, potentially impacting its reactivity and redox stability.23 Similar differences have also been observed in other 4d and 5d metal oxides, where the more radially-extended 5d orbitals of Ta5+ lead to more covalent character in Ta–O bonds compared to Nb–O,24 and the lighter MoO4 tetrahedra undergo a greater amount of rotational disorder compared to the heavier WO4 analogues.25
Pertechnetate salts of the type ATcO4 (A+ = NH4, K, Cs) were among the first Tc oxides isolated and studied.26–28 Subsequently, several powder diffraction and single crystal studies have reported the structures and phase transitions in ATcO4 and AReO4 salts.26–31 Recent synchrotron X-ray diffraction (SXRD) studies of some ATcO4 and AReO4 oxides by Chay et al., Kennedy et al., and Weaver et al. showed that most are isostructural at room temperature.30,32,33 For both Tc and Re, the structures of the four salts (A+ = Na, K, Rb, Ag) were refined to the CaWO4, tetragonal scheelite, structure (space group I41/a, #88). In this structure, the AO8 dodecahedra contain two distinct A–O distances, whereas the BO4 tetrahedra possess a single B–O distance (see Fig. 1).34 The BO4 tetrahedra are isolated from each other and connected to the AO8 dodecahedra via bridging oxygen anions. Bastide showed that the structures of various ABO4 compounds can be predicted from the rA/rO and rB/rO ratios, where rX is the ionic radii of the A- and B-site cations and oxygen anion, respectively.35 This explains why the ATcO4 and AReO4 salts mostly exhibit the same room temperature structure (rTc/rO ∼ rRe/rO), as well as explaining why RbTcO4 and RbReO4 (rRb/rO ∼ 1.13) exhibit the tetragonal I41/a scheelite-type structure whilst the CsTcO4 and CsReO4 (rCs/rO ∼ 1.23) analogues exhibit an orthorhombic pseudo-scheelite structure (space group Pnma, #62. See Fig. 1 and Fig. S1, ESI†).28,30,32–34,36,37 Exceptions to these ionic radii rules have been observed, as reflected in the rich structural chemistry of the LnTaO4 (Ln3+ = lanthanoid La-Lu) oxides that is not replicated in the corresponding LnNbO4 oxides.38–43
Although Re7+ and Tc7+ are similar in size, previous studies have described differences between the ATcO4 and AReO4 salts. Reynolds et al. highlighted differences in the thermal behaviour between NH4TcO4 and NH4ReO4 attributed to subtle differences in the B–O (B7+ = Tc, Re) bonds within the BO4 tetrahedra.23 More dramatic differences are evident in the two thallium salts TlTcO4 and TlReO4. The former adopts the orthorhombic Pnma pseudo-scheelite structure at room temperature, whereas TlReO4 exhibits a monoclinic P21/c (space group #14) structure (Fig. 1).32,33,44 The temperature dependence of the two structures is also different, with TlReO4 exhibiting an unusual re-entrant phase transition.45–47 A neutron total scattering study by Saura-Múzquiz et al. revealed that the local structure of TlReO4 was always best fit to the monoclinic P21/c space group, irrespective of the long-range average structure, with a change in the correlation length of the Tl+ 6s2 lone pairs used to explain the re-entrant phase transition.48 This behaviour appears to be unique to TlReO4 and has not been reported in the closely related oxide TlTcO4 that undergoes an orthorhombic-to-tetragonal phase transition upon heating.33
In this work, the temperature dependence of the structure of TlTcO4 was investigated using neutron total scattering and pair distribution function (PDF) analysis. The thermal behaviour of TlTcO4 is analysed, and in-depth modelling of its long-range average and short-range local structure is presented. These results are then compared to TlReO4 and other 6s2-containing ABO4 materials to determine the origin of the differences between these materials.
*** Caution! *** Thallium salts are potentially fatal if swallowed. Appropriate personal protective equipment and safe work procedures were employed during sample preparation and sample manipulation.
Thallium pertechnetate was prepared by cation metathesis from KTcO4 and TlF in deionised water. A thallium(I) fluoride (1.223 g, 5.48 mmol) solution was added dropwise to a solution of KTcO4 (1.095 g, 5.42 mmol) and TlTcO4 precipitated immediately. The mixture was mechanically stirred at room temperature for 1 hour. The solid was allowed to settle, the supernatant was discarded, and the resulting white solid was washed with cold deionised water (3 × 5 mL), isopropanol (3 × 5 mL), and diethyl ether (3 × 5 mL) before being dried at 120 °C for 12 hours. After cooling, TlTcO4 was removed from the vial, weighed (1.477 g, 4.02 mmol, yield 74.2%) and finally sealed in a 6 mm vanadium can and shipped to Oak Ridge National Laboratory (ORNL) for neutron scattering measurements.
Neutron scattering data suitable for neutron PDF (NPDF) analysis were collected with a centre wavelength of 0.8 Å using the POWGEN diffractometer at the Spallation Neutron Source (SNS) located at ORNL.49 The data were corrected for instrument background, the incident neutron spectrum, absorption, and multiple scattering events before normalisation. The software PDFgetN3 was used to calculate the NPDF with a Qmax of 30.0 Å−1. A Savitzky–Golay filter was applied to the S(Q) for Q > 9.0 Å−1 to reduce the noise in the broad regions of the reciprocal space data. Additional details are provided in the ESI (Fig. S2–S13).†
Structural refinements using the Rietveld method were carried out using the program TOPAS 6.50 The background of each pattern was estimated using a 12th-order Chebyshev polynomial. The scale factor, lattice parameters, atomic positions, and atomic displacement parameters (ADPs) were refined simultaneously with the peak profile parameters. Partial NPDFs were calculated using RMCProfile,51 and the crystal structures were drawn using VESTA.52
The theoretical calculations for TlTcO4 were performed using the projector augmented wave (PAW) method,53 as implemented in the Vienna Ab initio Simulation Package (VASP) version 5.4.4.54 The Perdew–Burke–Ernzerhof (PBE) functional revised for solids (PBEsol) and the Armiento-Mattsson 2005 (AM05) exchange–correlation functionals were employed to study the ground state structure and electronic interactions.55–57 The basis sets considered for each atom were as follows: Tl: 6s26p1, Tc: 4d55s2, O: 2s22p4. The plane-wave cutoff energy was set to 700 eV. A k-point mesh of 11 × 11 × 7 was used in the irreducible Brillouin zone, utilising Monkhorst Pack (MP) grids.58 The force minimisation and convergence criteria were set to 10−9 eV for electronic convergence and 10−3 eV for ionic relaxation, respectively. To compute the electronic structure of TlTcO4, the Tran–Blaha modified Becke–Johnson (TB-mBJ) potential was employed in conjunction with the traditional generalised gradient approximation (GGA) method incorporated in VASP.59–62 The TB-mBJ potential is known to effectively address the band gap of materials, particularly for insulators and semiconductors. The impact of spin–orbit coupling (SOC) was also considered. Our results indicate that the inclusion of SOC has almost no impact on the partial density of states (PDOS) for the atoms. Isosurface maps illustrating bonding interactions through the electron localisation function (ELF) were generated for the structures of TlTcO4 determined from the NPDF fits at 50 K. The visualisation of the ELF was subsequently performed using VESTA.52
A previous SXRD study concluded that TlTcO4 underwent a phase transition from orthorhombic Pnma to tetragonal I41/amd at ∼500 K.33 This was based on a lack of evidence of the (114) reflection, which is allowed in I41/a but forbidden in I41/amd by the 2h + l = 4n reflection condition. The orthorhombic Pnma to tetragonal I41/amd phase transition is allowed to be continuous (by condensation of the M4+ mode, see Fig. S1, ESI†). However, other orthorhombic Pnma pseudo-scheelite structures (such as CsBO4, B7+ = Tc, Ru, Re, Os; and RbRuO4) display a first-order phase transition from orthorhombic Pnma to tetragonal I41/a.32,33,36,63 These two tetragonal structures differ by the rotation of the BO4 tetrahedra that is allowed in I41/a but not in I41/amd. Attempts to fit both models to the high temperature 600 K dataset returned similar qualities of fit (Rwp = 2.00% for I41/a vs. 2.02% for I41/amd; see Fig. S15, ESI†). Therefore, the tetragonal I41/amd model was adopted at 600 K in line with the suggestion of Kennedy et al.33 The quality of the fit to the 600 K data demonstrates the quality of the sample, further highlighting the anomalous nature of the 50 K dataset. The Tl+ cation is 12-coordinate in the tetragonal I41/amd structure with four Tl–O bonds at 2.942(7) and a further eight at 3.352(4). The bond distortion index of the TlO12 polyhedra is noticeably smaller than seen at 50 K (0.0567 vs. 0.08205).
The temperature dependence of the unit cell parameters and Tc–O distances is displayed in Fig. 3. Due to the sparsity of data points in the current NPD study, values from the previous SXRD study are also presented,33 and these are in good agreement with the current study. The lattice parameters display conventional thermal expansion, with the a and b lattice parameters becoming equal above 500 K, corresponding to the orthorhombic-to-tetragonal phase transition. The temperature dependence of the ADPs shows the values for the Tc7+ and O2− atoms increasing between 50 and 300 K before dropping at 450 K (Fig. S16, ESI†). This is thought to be a result of the rare isosymmetric orthorhombic Pnma to orthorhombic Pnma phase transition at ∼400 K.33
![]() | ||
Fig. 3 Temperature dependence of the (a) and (b) lattice parameters, (c) unit cell volume, and (d) the Tc–O bond distances for TlTcO4 derived from the Rietveld refinements of the neutron powder diffraction data. Where not shown, the error bars are smaller than the symbols. The dashed lines in (a)–(c) are the results from the previous synchrotron X-ray diffraction study.33 The dotted line in (c) drawn to guide the eyes. |
As evident from Fig. 3, the unit cell volume of the tetragonal I41/amd phase is significantly larger than that of the extrapolated orthorhombic Pnma phase. Although only one temperature point is taken above the phase transition in the current NPD study, the previous SXRD refinements show the same trend.33 Similar behaviour was recently described in TlReO4, and Saura-Múzquiz et al. proposed that ordering of the Tl+ 6s2 lone pairs in the lower symmetry phase caused a more efficient packing of the structure.48 By analogy, we postulate that the 6s2 lone pairs in TlTcO4 are ordered in the orthorhombic phase and disordered in the tetragonal phase (Fig. S17, ESI†). This hypothesis is supported through a comparison of the monoclinic I2/b and tetragonal I41/a structures of NdNbO4 (Nd3+: 6s0) and BiVO4 (Bi3+: 6s2).24,64,65 The tetragonal I41/a unit cell volume in NdNbO4 is smaller than the unit cell volume extrapolated from the monoclinic I2/b phase. In contrast, it is larger in BiVO4 due to the disordering of the 6s2 lone pairs (Fig. S18, ESI†).
An additional advantage of NPD is that neutrons are scattered isotropically due to the short range of their interaction. This makes them more sensitive to the small displacements of atoms from their equilibrium positions, leading to greater precision in the refined ADPs compared to SXRD studies.66 Although oxygen is a relatively weak X-ray scatterer (Z = 8) compared to Tl+ (Z = 81) and Tc7+ (Z = 43), the coherent neutron scattering length of the three atoms are comparable (bO = 5.8, bTc = 6.8 fm, bTl = 8.8), leading to greater sensitivity to the Tc–O bond distances. As observed in Fig. 3, the Tc–O bond lengths decrease upon heating. This has been observed in several other ABO4 structures, including CaMoO4 and CsReO4,36,67 and has been correlated with the onset of rotational disorder in the rigid BO4 tetrahedra units. It is possible that the orthorhombic-to-tetragonal phase transition in TlTcO4 may be related to the onset of rotational disorder in the rigid TcO4 units. However, observation of the short-range local structure is required to verify this possibility.
![]() | ||
Fig. 4 Variable temperature neutron pair distribution function data of TlTcO4 at 50 (blue), 300 (green), 450 (orange), and 600 K (red). |
Stereochemically active ns2 lone pairs often cause cations to disorder through well-defined displacements from their centrosymmetric positions (giving a static contribution to the ADPs) as well as increased thermal motion (giving a dynamic contribution).69–71 Data were measured at 50 K to minimise the thermal disorder, resulting in sharper peaks at high r. Fits were conducted using the long-range average orthorhombic Pnma model over a range of 1.0–30 Å (Fig. 5). This model returned a poor fit especially at low r, with a large residual observed across the entire fitting range. The failure to accurately fit the high r region to the orthorhombic model suggests a subtle symmetry lowering. The feature at ∼2.84 Å was particularly poorly fitted, with partial PDFs revealing this region corresponds to the Tl–O distances within the TlO10 polyhedra and the O–O distances across the TlO10 polyhedra and TcO4 tetrahedra (Fig. S19, ESI†). Fits were attempted with both rigid (single Tc–O distance and O–Tc–O angles constrained to be 109.5°) and flexible TcO4 tetrahedra. However, neither approach was able to reproduce the sharp feature at ∼2.84 Å. Fits were also performed to the higher symmetry tetragonal I41/a and I41/amd scheelite-type structures, which also returned poor fits (Fig. S20, ESI†).
Fits were then undertaken using the monoclinic P21/c structural model seen in other 6s2-containing ABO4 materials, such as the long-range average and short-range local structure of TlReO4 at room temperature and the short-range local structure of PbWO4.48,72 This monoclinic model involves more complex distortion of the TlO10 polyhedra associated with short-range disorder of the Tl+ 6s2 lone pair electrons and provides the best fit to the data over a range of 1.0–30 Å, as shown in Fig. 5. This implies that the local structures of TlTcO4 and TlReO4 are the same despite their different long-range average structures. A key feature of the monoclinic P21/c TlReO4 structure is the displacement of the Tl+ cations from the centre of the Tl+ polyhedra, resulting in more significant variability in the Tl–O bond lengths and rotation of the ReO4 tetrahedra. Evidently, similar displacements occur in TlTcO4, with the same features observed in the 100 K dataset taken on the NOMAD diffractometer (Fig. S21 and S22, ESI†), as well as the 300 and 450 K NPDF datasets (Fig. S23–S25, ESI†). As detailed below, density functional theory (DFT) calculations show the ground state energy of the monoclinic P21/c structure to be lower than that of the orthorhombic Pnma structure (Tables S2 and S3, ESI†).
Fits were also performed across different distance ‘windows’, with better fits achieved by using 10 Å windows across the data (Fig. S26, ESI†). This also resulted in a more distorted structure at low r and differences in the TcO4 volume. This suggests nanodomains may form at lower temperatures related to the ‘freezing in’ of Tl+ 6s2 lone pair disorder and the isosymmetric orthorhombic Pnma to orthorhombic Pnma phase transition that are difficult to model, also potentially explaining the suboptimal fitting to the NPD at 50 K. In fitting the NPDF in 10 Å windows, the less-local fits provide an ‘averaging’ of these more-local features to better fit the data. Further confirmation using additional techniques, such as transmission electron microscopy or solid-state nuclear magnetic resonance spectroscopy, would be required to better model these contributions.
There is an apparent disconnect between the long-range average and short-range local structures of TlTcO4. At 50 K, the long-range NPD Bragg data were best fit using the orthorhombic Pnma model consisting of TlO10 polyhedra with several different Tl–O bond lengths and distorted TcO4 tetrahedra. The short-range NPDF data revealed that the TcO4 tetrahedra are rigid and a greater distribution of Tl–O bond distances within the Tl+ polyhedra (Fig. S27, ESI†), where the long and flexible Tl–O bonds allow for the rigid TcO4 tetrahedra to rotate freely. This degree of rotational freedom in the ABO4 scheelite-derived structures is due to the A–O–B connectivity that acts as a ‘hinge’.67 Unlike in perovskite-derived structures where the BX6 octahedra are corner sharing, the TcO4 tetrahedra are not directly connected, relaxing the requirement of coherent tetrahedra rotations. It is this lack of connectivity that leads to incoherent TcO4 tetrahedra rotations and translations that Rietveld refinements may model as a combination of increased ADPs and Tc–O bond shortening, as also seen in RbReO4 and CsReO4.34,36 Furthermore, the monoclinic P21/c structural model can be derived from the tetragonal I41/amd aristotype by rotation of the TcO4 tetrahedra about the a-axis (similar to those observed in the orthorhombic Pnma structure) governed by a M4+ mode, and about the b-axis (analogous to that observed in the tetragonal I41/a structure) governed by a Γ3+ mode (Fig. S1, ESI†). It is postulated that locally both the M4+ and Γ3+ modes are present, and the lack of direct connectivity reduces the correlation of the rotations. This explains why the long-range average structure is best fitted to the orthorhombic Pnma (M4+) structure, as the TcO4 tetrahedra rotations governed by the Γ3+ mode are uncorrelated and separated by ∼18 Å.
The fitting process was repeated for the 600 K dataset (Fig. S28, ESI†), recalling that 600 K is above the orthorhombic-to-tetragonal phase transition.33 All four models (orthorhombic Pnma, tetragonal I41/a and I41/amd, and monoclinic P21/c) returned a similar goodness of fit, most likely due to the greater amount of thermal motion and broader NPDF peaks at higher temperatures. Despite these broadened features, the Tc–O peak at ∼1.72 Å remains relatively unchanged between the 50 and 600 K datasets, while the prominent peak at ∼2.84 Å shows significantly less thermal broadening compared to the other peaks within the NPDF datasets reflecting the rigid TcO4 tetrahedra (Fig. 4). Although the local features at ∼2.84 Å are best fit using the monoclinic P21/c model, the less local (10–30 Å) range is fit equally well between the four models (Fig. S29, ESI†). This suggests that the Tl+ 6s2 lone pairs are still sterically active at high temperatures, but with a short (>10 Å) correlation length.
Fig. 6 illustrates the band structure and PDOS for TlTcO4 in the monoclinic P21/c structure. Essentially identical results were obtained for calculations assuming the orthorhombic Pnma structure (Fig. S30, ESI†). This shows that TlTcO4 is a direct band gap semiconductor with a band gap of 2.53 eV, as estimated using the TB-mBJ potential. The valence band is predominantly compromised of the O p and Tl s states, while the conduction band is mainly derived from Tc d states. The Tl 6s states are present at the top of the valence band and show significant hybridisation with the O 2p states, with minimal involvement from the Tc states. This is a characteristic of stereochemical active 6s2 electrons.74,75 This is demonstrated through ELF calculations, illustrated for the 〈001〉 plane in Fig. 7, that indicate localisation of the Tl 6s states and presumed stereochemical activity. Conversely, the O 2p states exhibit moderate localisation, while the Tc states show very little.
![]() | ||
Fig. 8 The A–O, B–O bond lengths, and O–B–O bond angles from (a)–(c) TlTcO4 from this work, (d)–(f) TlReO4 taken from Saura-Múzquiz et al.,48 (g)–(i) PbWO4 taken from Mullens et al.,72 and (j)–(l) BiVO4 taken from Sleight et al.76 All data is taken from the neutron pair distribution function fits at 300 K to the monoclinic P21/c model, except for BiVO4 (crossed histogram bars), which is the long-range average structure in the monoclinic I2/b space group. |
The A–O–B connectivity makes it reasonable to assume that distortions of the AO8 polyhedra and BO4 tetrahedra are related. Across the 5p block elements (A = Tl+, Pb2+, Bi3+), the oxidation state increases resulting in shorter A–O bonds and smaller AO8 polyhedra (Table 1). As seen in Fig. 8, the larger oxidation state A-site cations result in less flexible AO8 polyhedra and a smaller distribution of the A–O bond lengths. As the AO8 polyhedra get smaller, the flexibility of the BO4 tetrahedra becomes more important. While the TcO4 and ReO4 tetrahedra remain rigid, the WO4 tetrahedra distort through relaxation of the O–W–O bond angles and the VO4 tetrahedra distort through variations in both the O–V–O bond angles and the V–O bond lengths. Likewise, the choice of nd0 element on the B-site affects the degree of polarizability, with the VO4 tetrahedra in BiVO4 (V: 3d0) exhibiting a greater amount of distortion than the TcO4 tetrahedra in TlTcO4 (Tc: 4d0) and the ReO4 tetrahedra in TlReO4 (Re: 5d0).
The PDOS for TlTcO4 is very similar to that of TlReO4 and PbWO4 (Fig. S31, ESI†), despite these materials having different long-range average structures at 300 K.48,72 The PDOS is, however, significantly different from that reported for BiVO4 where the Bi 6s states are significantly further below the Fermi level (∼9 eV vs. 1 eV for TlTcO4).73 Presumably strong hybridisation of the 6s states of Pb and Tl with O p or metal d states pushes these levels closer to the conduction band edge, whereas the higher oxidation state of Bi reduces the hybridisation of the Bi 6s states with the oxygen states. Indeed, near the Fermi level BiVO4 features only O p states in the valence bands. The implication is that the 6s electrons do not participate significantly in bonding, rather they remain as lone pairs and exhibit an obvious stereochemical influence. It can be speculated that this difference in the PDOS may be correlated with the photocatalytic properties of these oxides.
The question remains as to why TlTcO4 and TlReO4 have the same monoclinic P21/c local structure but different long-range average structures. Previous studies have shown that all other ATcO4 and AReO4 salts (A+ = Na, K, Rb, Ag, Cs) are isostructural.30,32,33 Bastide's analysis of the rA/rO and rB/rO ratios in the ABO4 materials suggests that TlReO4 should adopt a tetragonal I41/a structure similar to RbReO4,35 as the two A-site cations are approximately the same size (Tl+ = 1.59 vs. Rb+ = 1.61 Å).17 TlReO4 exhibits the tetragonal I41/a structure at low (<150 K) and high temperatures (>400 K), with the room temperature monoclinic P21/c structure the result of long-range ordering of the Tl+ 6s2 lone pairs.48 Since Tc7+ is slightly smaller than Re7+ (0.37 vs. 0.38 Å),17 as evident from the sharp B–O peaks at 1.72 and 1.74 Å in the NPDF data of TlTcO4 and TlReO4 respectively, it is possible that this small change in the B-site cation ionic radii is sufficient to favour the formation of the long-range orthorhombic Pnma, rather than the tetragonal I41/a structure (rTc/rO ∼ 0.27 vs. rRe/rO ∼ 0.28).
The long-range average monoclinic P21/c structure of TlReO4 at room temperature is derived from the tetragonal I41/amd aristotype by rotations of the ReO4 tetrahedra associated with the irreps Γ3+ and M4+.32 In this structure, TlReO4 has two distinct ReO4 tetrahedra sites rotated about the b-axis (associated with the irrep Γ3+, analogous to the tetragonal I41/a to tetragonal I41/amd phase transition), and one distinct ReO4 tetrahedra site rotated about the a-axis (associated with the irrep M4+, analogous to the orthorhombic Pnma to tetragonal I41/amd phase transition). Although the same irreps are present in the local structure of TlTcO4, the opposite trend is observed. There are two distinct TcO4 tetrahedra sites rotated about the a-axis, and one rotated about the b-axis (Fig. S32, ESI†). Evidently, the ‘averaging’ of these leads to the observation of different long-range average structures.
A further difference between TlTcO4 and TlReO4 is the lack of a re-entrant phase transition observed in TlTcO4. Phase transitions can be driven by the ability of the 6s2 lone pairs to macroscopically align, leading to a lowering of symmetry and a decrease in unit cell volume due to the more compact structure. This is seen in TlReO4 where the room temperature structure is monoclinic P21/c due to the ability of the 6s2 to align macroscopically, before an increase in temperature leads to greater thermal motion and a breakdown of this long-range ordering.48 No such macroscopic ordering is observed in TlTcO4 due to the larger thermal motion of the TcO4 tetrahedra compared to ReO4 (Fig. S33, ESI†). This was similarly seen in work by Amarasinghe et al., where greater rotational disorder was observed in the lighter MoO4 tetrahedra of NaRE(BO4)2 (RE3+ = Ln, Y; B6+ = Mo, W) compared to the heavier WO4 tetrahedra (95 vs. 184 amu).25 As Tc7+ is almost half of the atomic mass of Re7+ (98 vs. 186 amu), a similar sensitivity to temperature and larger thermal displacements is expected. The greater susceptibility of the TcO4 tetrahedra to rotate may also explain the differences in the thermal expansion behaviour of the two ammonium salts NH4TcO4 and NH4ReO4 at different temperatures.23
Several anomalies were observed in the Rietveld refinements to the NPD data, such as large ADPs at low temperature and a shortening of the Tc–O bond distance upon heating. Using NPDF, these structural anomalies were resolved by devising a model with longer, more flexible Tl–O bonds that allow for the rotation of rigid TcO4 tetrahedra units. This model was then supported by DFT calculations. Considering both the long-range average and short-range local structures of TlTcO4, a more holistic structural model was developed. At low temperatures, the stereochemical activity of the Tl+ 6s2 lone pairs is ‘frozen’ into the structure, resulting in static rotations of the TcO4 tetrahedra. As the material is heated, the TcO4 tetrahedra begin to rotate as governed by the Γ3+ and M4+ modes. The Γ3+ rotations appear to be less correlated than the M4+ modes, resulting in the observation of the orthorhombic Pnma structure in the long range. At ∼500 K, the thermal motion is large enough such that the rotations shift from coherent to incoherent, resulting in an ‘averaging out’ of the rotations and the observation of the higher symmetry tetragonal I41/amd structure. However, the 6s2 lone pairs remain stereochemically active on the local scale, as observed in the NPDF fits. The orthorhombic Pnma to tetragonal I41/amd phase transition is seemingly the result of a change in the correlation length of the Tl+ 6s2 lone pairs, leading to a larger unit cell volume due to their uncorrelated displacements.
This study of TlTcO4 was put into the larger context of other ABO4 structures containing 6s2 lone pairs, including TlReO4, PbWO4, and BiVO4. Each of these oxides has a different long-range average structure at 300 K. However, common trends in the local scale distortions necessary to accommodate the stereochemically active 6s2 lone pair were identified. Although TlTcO4 and TlReO4 disorder in a similar way, a combination of the size of the A-site polyhedra and the choice of nd0 B-site cation leads to different disordering phenomena. These include large ranges of A–O bonds, flexibility in the O–B–O tetrahedra bond angles, and distortion of the B–O bond lengths. This gives a rich tapestry of mechanisms that can be further implemented when fine-tuning the local scale interactions of atoms for the next generation of functional energy materials. The role of polarizability of the B-site cations is worthy of further exploration, especially with the large distortions of the VO4 tetrahedra observed and the lack of NPDF data available for BiVO4. Unfortunately, no such BiBO4 analogous structures are currently known, with BiBO4 [B5+ = Nb (4d0), Sb (4d10), Ta (5d0)] crystallising in different structures consisting of BO6 octahedra.77–79
Optimisation of the structures using the commonly employed PBE, PBEsol, and AM05 exchange–correlation functionals within VASP indicated that the monoclinic P21/c structure was more stable than the orthorhombic Pnma model that describes the long-range average structure of TlTcO4 below 500 K. However, for each of these functionals the tetragonal I41/a structure was predicted to be the energetic ground state, demonstrating the challenges of accurately describing the 6s2 lone pair electrons and polarizability of the Tl+ cation. While it may be possible to identify a hybrid functional that correctly predicts the relative energy of these phases, it is evident that there is no unique universal exchange–correlation functional appropriate for these types of oxides, as seen with other ABO4 metal oxides.38,73
Finally, this work adds to the suite of studies demonstrating the necessity of considering both the long-range average and short-range local structure of ABO4 scheelite-related materials, which have found application as photocatalysts, selective oxidation catalysts, and potential nuclear waste forms.33,76,80 Relationships such as these are of great significance, and additional experimental studies of the local structures, partnered with DFT calculations, are required to identify additional ways of fine-tuning these local disordering features.
Footnote |
† Electronic supplementary information (ESI) available: Variable temperature Rietveld refinements of TlTcO4, space group analysis at 600 K, atomic displacement parameters, depiction of structures, partial PDF fits, density functional theory calculations, and variable r range NPDF fits. See DOI: https://doi.org/10.1039/d4cp03707c |
This journal is © the Owner Societies 2024 |