Christopher J.
Ridley
*a,
Craig L.
Bull
ab and
Nicholas P.
Funnell
a
aISIS Neutron and Muon Source, Rutherford Appleton Laboratory, Chilton, Didcot OX11 0QX, UK. E-mail: christopher.ridley@stfc.ac.uk
bSchool of Chemistry, University of Edinburgh, David Brewster Road, Edinburgh EH9 3FJ, UK
First published on 22nd December 2021
We report the first full-structural analysis of two ambient temperature phases (II, III) of LiOD as a function of pressure, and present a direct-structural confirmation of a hydrogen-bonded network in the high-pressure phase-III. LiOD has been measured using neutron powder diffraction up to a maximum pressure of 4.1 GPa. The sample is observed to transform from tetragonal phase-II (P4/nmm) to phase-III at a pressure of 1.16(3) GPa. The previously suggested monoclinic structure of phase-III, isostructural with NaOD, is shown to strongly misfit the measured data, while a new tetragonal structure (I41/acd), suggested in the literature from a first-principles calculation, is found to fit the pattern extremely well. Neutron diffraction data from this new tetragonal structure agree with earlier spectroscopic evidence for a hydrogen-bonded network within the system. This tetragonal structure has not been observed in any other alkali-metal hydroxide before now. Details of the hydrogen bond lengths, lattice strains through the transition, and the compressibility of each phase are presented. This revised structure may have implications as to how LiOD is used for high-pressure/high-temperature materials synthesis.
LiOH is unusual as an alkali-hydroxide, in that the structure of phase-II doesn’t support, or perhaps more accurately, only has an extremely weakly formed hydrogen-bonded network, as the hydroxyl group arrangement leads to a large H⋯O separation (approximately 3.5 Å), with nearest neighbour hydroxyl groups being antialigned.6 In contrast, the heavier-alkali hydroxides (K, Rb, Cs) are more extensively hydrogen bonded,7 while NaOH is anomalous; the hydrogenous form of which shows no hydrogen bonding, while the deuterated form has a low-temperature phase with a hydrogen-bonded network. While all alkali-metal hydroxides have long O⋯O separations, NaOH is at the limit where the effects of hydrogen bonding on the structure are negligible, while the deuterated counterpart has a low-temperature phase transition to an extensively hydrogen-bonded state. When the structure is substituted with D, there is sufficient contraction in the hydrogen bond length to form a stronger hydrogen-bonded network at low temperatures, which is one contributing factor to the favourability of a new structure.8 This effect can similarly be induced in the hydrogenous form with a small amount of applied pressure at low temperatures.9 Unlike heavier alkali-hydroxides, which show low-temperature phases with zig-zag hydrogen-bonded chains,10 LiOH-II is the stable phase down to 10 K.11 Above 648 K at ambient pressure, a first-order transition to an unknown structure, phase-I, is observed close to the melting temperature, where the formation of Li2O is suppressed.11,12
Structurally, LiOH is unique among the other alkali-metal hydroxides. Under ambient conditions, it crystallises with tetragonal symmetry (P4/nmm),13 designated as phase-II. In this phase, each –OH group forms the cap of a square pyramid formed of Li atoms, which are each coordinated to four O atoms. Neighbouring –OH groups are anti-aligned in the ab-plane, and are aligned along the c-axis, analogous to the α-PbO structure14 (see Fig. 1). The closest related alkali-metal hydroxide, NaOH, instead crystallises with orthorhombic (Cmcm) symmetry with a similar pyramidal motif, but with the Na coordinated to five O atoms.15 KOH and RbOH both crystallise with monoclinic (P21/m) symmetry, but the –OH libration is such that exact H positions are difficult to determine, though the pyramidal motif is assumed in the average structure, with the metal ion being essentially coordinated to six O atoms, forming a distorted octahedral NaCl-type structure.16
Fig. 1 Three proposed structures of LiOH. (left) The known ambient pressure/temperature tetragonal phase-II (P4/nmm). The square pyramidal units are stacked anti-aligned in the ab-plane, and aligned along the c-axis. (middle) The proposed monoclinic structure (P21/c) of phase-III (Adams et al.18). The motif is changed, forming a more distorted Li substructure. Zig-zag hydrogen bonding chains (not shown) are formed along the b-axis. (right) proposed high-pressure phase-III tetragonal structure (I41/acd) (Hermann et al.21). For simplicity, the structure is shown only in the ac-plane, though the –OH capped, square pyramidal units from the ambient structure still feature, albeit flattened, and in a more complex arrangement. The hydrogen bonding in this tetragonal structure (not shown) forms linear chains along the a- and b-axes, in alternating stacked layers along c. Both high-pressure structures have considerably shorter O⋯H distances. The two structures proposed for phase-III are listed in Table 1. Note: the Li⋯O ‘bonds’ are illustrated differently to indicate that the Li+ is coordinated to four OH− through ‘lithium-bonding’.22 |
At ambient temperature, LiOH has been shown via IR and Raman spectroscopy to undergo a pressure-induced phase transition above approximately 0.7 GPa.17 Neutron diffraction was used to measure LiOD, and confirmed the phase transition to phase-III at 1.7 GPa, showing a large isotope effect. This was indexed with a monoclinic (P21/c) cell,18 assuming a structure similar to that of the low-temperature phase of NaOD where the Li remains coordinated to four O atoms (see Fig. 1, and Table 1), but the motif changes to an irregular form.19 This was inferred from relative peak intensities, as the data were of insufficient quality for a Rietveld refinement.18 Importantly, both the –OH and –OD Raman stretching bands were observed to broaden and soften slightly with pressure, indicating a slight weakening of the O–H bond, and the formation of hydrogen bonding in phase-III. Ab initio molecular dynamics confirm the monoclinic symmetry,20 and reproduce the observations from Raman and IR of the –OH mode, while first-principles density function theory calculations have suggested an alternative, energetically more favourable structure for phase-III.21 In this case the symmetry remains tetragonal (I41/acd), which instead suggests a rearrangement of the hydroxide ions, where the local environment remains very similar to that of phase-II, but with a more complex rotated arrangement of the pyramidal units. The pyramidal units themselves are also flattened, with the Li and O atoms essentially planar to each other.
The only previous structural study of this material was limited in scope due to the low data quality, and resolution, and therefore did not report any information on the structural evolution of LiOD-II with pressure, or the structure of the new phase-III.18 The present study has collected neutron powder-diffraction data at high pressure with improved quality and reports the refined structure of LiOD-III, and the evolution of LiOD-II leading up to it. This is therefore the first direct-structural evidence of the tetragonal phase-III, and the formation of a hydrogen-bonded network in LiOD. This may have important implications for the chemistry of Li-ion cathode materials synthesised from LiOH under high pressure.
High-pressure time-of-flight neutron powder diffraction was performed on the PEARL instrument at the ISIS Neutron and Muon Source.24 For the high-pressure data, the sample was loaded into a null-scattering encapsulated TiZr gasket, with a mixture (1:1 by volume) of perdeuterated iso-/n-pentanes added as a pressure transmitting medium. A V3 Paris-Edinburgh press was used to apply load to the sample, using single toroidal ZrO2/Al2O3 anvils. The sample pressure was determined from the equation of state of Li2CO3.25 All data were reduced and corrected for attenuation using Mantid,26 and Rietveld refinements were performed using Fullprof.27
Present study | Mair29 | Dachs13 | |
---|---|---|---|
NPD | NPD | NSXD | |
a-axis (Å) | 3.5458(2) | 3.549 | 3.55 |
c-axis (Å) | 4.3458(4) | 4.334 | 4.33 |
Volume (Å3) | 54.640(6) | 54.588 | 54.57 |
Oz | 0.1933(8) | 0.1938(4) | 0.1951(5) |
Dz | 0.3928(9) | 0.410(1) | 0.4069(14) |
B Li11 (Å2) | 0.49(6) | 0.38(2) | 0.85(12) |
B Li33 (Å2) | — | 0.93(5) | 2.4(3) |
B O11 (Å2) | 0.49(3) | 0.34(3) | 1.28(4) |
B O33 (Å2) | — | 0.68(2) | 2.01(11) |
B D11 (Å2) | 5.14(10) | 1.28(4) | 5.7(2) |
B D33 (Å2) | — | 0.85(4) | 2.8(2) |
Li–O (Å) | 1.959(2) | 1.963 | 1.96 |
O–D (Å) | 0.867(9) | 0.937(7) | 0.917(8) |
Ka,II = −2.62(12) × 10−2 Å GPa−1 |
Kc,II = −2.52(12) × 10−1 Å GPa−1 |
The D⋯O separation was determined to reduce from 3.431(11) to 3.298(12) Å, and the O–D bond reduced from approximately 0.891(11) to 0.835(12) Å over this pressure range. Full structural parameters and the EOS fit are provided in the ESI.†
The sample was observed to begin transforming from phase-II to -III at 1.16(3) GPa, where phase coexistence between the two indicated a first-order transition. By 1.25 GPa the sample was pure phase-III. The data at this pressure were fitted against both proposed models in Table 1. From a LeBail fit, the lattice parameters for the monoclinic structure were refined to be similar to those reported by Adams et al.18 However, when the NaOD-type structure was fitted to the data, there were clear intensity mismatches around principal reflections. Fully relaxing the structure led to little improvement in the overall fit and an unphysical model. In contrast, the tetragonal model of phase-III fitted much better, providing a closer intensity match to the peaks at approximately 2.6 Å. To verify the robustness of the model, the deuterium atom was moved away from the special position, and placed onto a general (x, y, z)-coordinate with reduced occupancy, and allowed to refine. Within error, it returned to the original position, fully occupying a 16e site, supporting the antiferroelectric-type arrangement of the –OH groups.
Fig. 3 shows a comparison of the fits against the two models and also reports the goodness-of-fit parameters. The value of the global χ2 drops from 4.98 for the monoclinic structure, to 1.13 when the tetragonal structure is fitted to the data, metrically demonstrating a significant improvement in fit. The phase transition results in a large volume discontinuity (approximately 10%), and significant reduction in compressibility of the sample, with V0 = 376.3(4) Å3 (Z = 16), K0 = 39.1(8) GPa. In phase-III the a/b-axes and c-axis compress at a comparable rate:
Ka,III = −4.91 (11) × 10−2 Å GPa−1 |
Kc,III = −4.93 (12) × 10−2 Å GPa−1 |
At 1.41(5) GPa the O–D bond length in phase-III is 1.122(12) Å which stays essentially unchanged with pressure, and the D⋯O separation is 1.869(12) Å, reducing only slightly to 1.807(11) Å by 4.1 GPa. This is significantly reduced from that in phase-II, and clearly indicative of hydrogen bonding in the system (this value is approximately 1.9 Å for ice VIII). The lengthening of the O–D bond, and shortening of the D⋯O separation, is clear evidence for the emergence of more extensive hydrogen bonding in this material in phase-III, consistent with spectroscopic observations of the O–H stretching frequency.17 No additional features were observed in the diffraction spectra up to the maximum pressure of 4.1 GPa. The refined atomic coordinates for this phase are listed in the ESI.†
The bulk modulus of phase-II, K0 = 11.4(7) GPa, is very low compared to a number of other group-I and group-II hydroxides, being typically 50% and 70% lower.37 with NaOH phase-III (also known as α-NaOH) being the one exception. NaOH phase-III has a reported K0 = 25.0(4) GPa, though with a negative pressure derivative K′ = −16.0(7).38 If the negative curvature in the low-pressure region is ignored, and the data from Beck and Lederer38 are refitted with a second-order Birch–Murnaghan EOS, then K0 ≈ 12 GPa, closer to the value determined for LiOD in the present study, though this completely fails to fit the data below 0.2 GPa for NaOH phase-III. As the ambient forms of KOH, RbOH, and CsOH have similar compressibilities with a hydrogen-bonded network, suggesting that electrostatics don't dominate this property, the same might be expected for LiOD and NaOH without a hydrogen-bonded network. Their high-pressure forms contrast this, with LiOD phase-III being significantly less compressible (K0 = 39.1(8) GPa) than NaOH phase-V (K0 = 29.7(4) GPa,37), despite having similar transition pressures, and LiOD phase-III having approximately 20% lower density at ∼ 1 GPa. While one might expect the high-density form of LiOD-III to be more compressible than NaOH-V based purely on electrostatics, this is not found to be the case, indicating that the hydrogen- and ionic-bonding play a much more balanced role in the behaviour of this material.
Relating the difference in compressibilities to the hydrogen bonding in LiOD is non-trivial. Both NaOH/D phase-V have been shown to demonstrate no softening of the O–H stretching mode with pressure, thought to be due to bent hydrogen bonds.9 This was later demonstrated to be the case in NaOD by Loveday et al.,34 who also posited that if both nearest and next-nearest O atoms are part of the hydrogen-bonded network, then the latter is weakened in favour of the former with further applied pressure. In LiOD phase-III, the hydrogen bonding is purely linear, and this seemingly contributes to a lower overall compressibility compared with the canted structure of NaOH phase-V.
The length of the hydrogen bond in phase-III is shorter than that found in ice VIII, and so might conventionally be called a ‘strong’ hydrogen bond, however, this fails to consider the effect of strain on the system caused by the repulsion between the terminal O2− ions.39 The symmetry of the phase-III structure constrains the hydrogen bond to be linear (∠O–D⋯O ≡ 180o), such that each O atom has a nearest neighbour parallel and perpendicular to it. The O⋯O separation in LiOD phase-III was determined to be 2.992(14) Å parallel-to, and 2.9933(6) Å perpendicular-to the hydrogen bonds, this results in a slightly more symmetric bond than would be energetically optimal, and is subsequently less ‘strong’ than in ice.39 This reasoning lies in good agreement with the measured O–H stretching frequency from Raman and infrared,17 and the calculated weighted bond valence sum at 1.41(5) GPa is 0.864 v.u., consistent with the trend observed by Lutz et al.25
The transition between phase-II and -III of LiOD is reconstructive, whereby the system restructures through a rotation of the tetrahedral units with the O–H bond moving into the ab-plane, and a change in translational symmetry. As such, there is no direct group-subgroup relation between the two phases. It is therefore expected that the transition is driven by homogeneous strain on the lattice. As detailed by Carpenter et al.,40 the variation of lattice parameters at a phase transition can be interpreted in terms of macroscopic strain, separate from the expected changes due to thermal expansion or length contraction due to applied stress, caused by the microscopic rearrangement of the structure. While the strains in LiOD cannot be directly allocated to an active symmetry mode, due to the reconstructive nature of the transition, the determined strains are still informative, as it is assumed that the progression between phase-II and -III occurs through an intermediate subgroup common to both. As both phases are tetragonal, it can be assumed that the unstrained sample would have cubic-symmetry, allowing the calculation of lattice strains relative to a pseudo-cubic cell (axis length, ao, see ESI†) in each phase. The calculated strain (see ESI†) is analogous to the ‘tetragonal strain’ of the system, though is corrected for the expected lattice contraction due to pressure. The absolute values between the two phases cannot be directly compared, as the rearrangement of the system does not allow for any meaningful correction for different values of Z. The volumetric strain can be corrected for Z, and shows an overall reduction in strain at the transition, and qualitatively εs reduces with applied pressure in phase-II, changes discontinuously at the point of transition, and then increases approximately linearly with further applied pressure in phase-III (Fig. 4). This shows that decreasing the c-axis length in phase-II, and therefore D⋯O separation, reduces the strain in the system. When reduced sufficiently the system then rearranges into a more complex hydrogen-bonded network where the cell is more constrained, resulting in an increase in strain with further applied pressure. This positive slope of εs for phase-III suggests that the system tends towards a lower symmetry form at higher pressure still; this would be consistent with a transition to LiOH phase-IV at higher pressures, though we have insufficient data to expand on this further.21
The onset of phase-III is found to occur at approximately 0.6 GPa lower pressure than previously reported for LiOD.18 The reason for this discrepancy is not clear but may be related to differences in the pressure medium used between the two studies. Fluorinert is known to be non-hydrostatic at much lower pressures than a mixture of pentanes. Alternatively, it may be related to the rate of compression differing between the experiments.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1ma01024g |
This journal is © The Royal Society of Chemistry 2022 |