Introduction
Compounds that belong to the oxyhydroxide group (M 3+OOH, M = Al, Sc, V, Cr, Mn, Fe, Ga, Rh and In) have been investigated widely under non ambient conditions, such as high pressure (HP) and/or variable temperature (Li et al., Reference Li, Snoussi, Li, Wang, Yang and Gao2010; Sano-Furukawa et al., Reference Sano-Furukawa, Yagi, Okada, Gotou and Kikegawa2012; Suzuki, Reference Suzuki2013, Reference Suzuki2018; Zhang et al., Reference Zhang, Cui, Yang, Zhang, Tang, Wu and Cui2017; Shito et al., Reference Shito, Okamoto, Sato, Watanabe, Ohashi, Fuchizaki, Kuribayashi and Suzuki2019; Tang et al., Reference Tang, Huang, Zhang and Cui2020; Ito et al., Reference Ito, Ikeda, Sakamaki, Kuribayashi and Suzuki2022). A common phase transition sequence observed in this group is from the diaspore-type structure (α-AlOOH, Pnma) at ambient pressure to a high pressure InOOH-type structure (Pmn21) (see Figs 1 and 2), and then further to the YOOH-type structure (P21/m) (Bolotina et al., Reference Bolotina, Molchanov, Dyuzheva, Lityagina and Bendeliani2008; Ito et al., Reference Ito, Ikeda, Ban, Kubota, Sakamaki, Kuribayashi and Suzuki2021). The exact values of the transition pressures between the α-AlOOH-type and InOOH-type and then between the InOOH-type and YOOH-type differ significantly depending on the composition. For example, the transitions between the α-AlOOH-type and the InOOH-type for GaOOH, AlOOH and ScOOH occurs at 8.5 GPa, 18.0 GPa, and 4.4 GPa, respectively but the phase transition from the InOOH-type to the YOOH-type for InOOH and ScOOH occurs at 50.8 GPa and 8.1 GPa, respectively (Ito et al., Reference Ito, Ikeda, Ban, Kubota, Sakamaki, Kuribayashi and Suzuki2021).
From a mineralogical perspective, (M 3+OOH) oxyhydroxide minerals are interesting as potential carriers of water into the Earth's lower mantle (Verma et al., Reference Verma, Modak and Stixrude2018). The deep Earth mantle was originally thought to be devoid of water due to the observed decomposition of hydrous minerals at high temperatures and pressures at depths of ~150 km (5 GPa). Recent discoveries, however, have revealed hydrous minerals stable at much greater pressures, some retaining water even in the lower mantle. Researchers are now actively investigating the high-pressure behaviours of dense hydrous magnesium silicates, including phases A, superhydrous phase B, phase D and phase H, to understand their role in the deep Earth's water cycle (Xu et al., Reference Xu, Inoue, Kakizawa, Noda and Gao2021; Li et al., Reference Li, Speziale, Koch-Mueller, Husband and Liermann2022).
Diaspore is also a mineral of practical importance. Cryptocrystalline or earthy diaspore is an important component of bauxite, the most important aluminium ore. It forms as the final product of intense weathering of various types of rocks. Its prismatic or tabular crystals are found in some alkali pegmatites and products of hydrothermal alteration of aluminosilicate minerals. Transparent, well-formed diaspore crystals, that sometimes change colour under different lighting, are cut as gemstones (‘zultanite’).
The first crystallographic investigations of diaspore were performed almost a century ago (Ewing, Reference Ewing1935; Busing and Levy, Reference Busing and Levy1958). α-AlOOH diaspore is the thermodynamically stable phase at low pressure. Under equilibrium conditions (high temperature) it transforms to δ-AlOOH above 18 GPa (Suzuki et al., Reference Suzuki, Ohtani and Kamada2000; Komatsu et al., Reference Komatsu, Kuribayashi, Sano, Ohtani and Kudoh2006; Simonova et al., Reference Simonova, Bykova, Bykov, Kawazoe, Simonov, Dubrovinskaia and Dubrovinsky2020), but when compressed at ambient temperature, it metastably retains the original structure to pressures as high as 50 GPa (Friedrich et al., Reference Friedrich, Haussuehl, Boehler, Morgenroth, Juarez-Arellano and Winkler2007b). Diaspore crystallises in space group 62. In the ICSD [Inorganic Crystal Structure Database, https://icsd.products.fiz-karlsruhe.de/] all deposited diaspore structures are described in Pbnm, however the alternative, equivalent Pnma setting better complies with the current IUCr edition of International Tables for Crystallography (Hahn, Reference Hahn2002). The stability of this structure within various pressure ranges and potential phase transitions has already been investigated widely by various groups (Xu et al., Reference Xu, Hu, Ming, Huang and Xie1994; Pawley et al., Reference Pawley, Redfern and Holland1996; Grevel et al., Reference Grevel, Burchard, Fasshauer and Peun2000; Friedrich et al., Reference Friedrich, Wilson, Haussuehl, Winkler, Morgenroth, Refson and Milman2007a; Reference Friedrich, Haussuehl, Boehler, Morgenroth, Juarez-Arellano and Winkler2007b; Cedillo et al., Reference Cedillo, Torrent and Cortona2016; Sugiura et al., Reference Sugiura, Arima, Nagai and Sugiyama2018; Huang et al., Reference Huang, Xu, Liu, Li, Ye, Chen, Kuang, Chi, Fan, Ma and Zhou2021).
δ-AlOOH has been investigated systematically and its equilibrium p–T phase boundary with diaspore has been established (Sano-Furukawa et al., Reference Sano-Furukawa, Kagi, Nagai, Nakano, Fukura, Ushijima, Iizuka, Ohtani and Yagi2009). The differences in compressibility between δ-AlOOH and the deuterated form δ-AlOOD have also been studied (Suzuki, Reference Suzuki2022). It is known that δ-AlOOH undergoes further subtle phase transitions at pressures from 6.1 to 8.2 GPa (Kuribayashi et al., Reference Kuribayashi, Sano-Furukawa and Nagase2014). Dense hydrous magnesium silicates within subducting slabs undergo a phase transition leading to formation of a mixture of H2O and anhydrous minerals. δ-AlOOH, however, has demonstrable stability within the lower mantle's pressure–temperature conditions, though questions remain about its behaviour in the lowermost mantle, where temperature increases rapidly with depth. Its thermodynamic stability extends to pressures and temperatures up to 134 GPa and 2300 K, with theoretical predictions suggesting a transition to a pyrite-type structure above 170 GPa, confirmed experimentally at ~190 GPa (Ohira et al., Reference Ohira, Ohtani, Sakai, Miyahara, Hirao, Ohishi and Nishijima2014). However, the formation of pyrite-type AlOOH remains unverified, highlighting the need for further high-pressure, high-temperature experiments. Recent high-pressure experiments propose that δ-AlOOH, which is structurally similar to phase H, could remain stable in the lower mantle (Thompson et al., Reference Thompson, Campbell and Tsuchiya2021). It has also been shown that Al can separate from bridgmanite and CaSiO3 perovskite to form δ-AlOOH in the presence of H2O. Thus, δ-AlOOH is a promising candidate as a primary carrier of water deep within the Earth.
The first study describing electron density distribution in diaspore (α-AlOOH) was published in 1979 (Hill, Reference Hill1979), in which, on the basis of 597 reflections, the valence expansion/contraction coefficients, as well as atomic charges were estimated. In this paper, we compare electron density distribution in the α-AlOOH and δ-AlOOH phases. In this study we performed a new high-pressure, high-resolution synchrotron XRD experiment at 2.5 GPa, and analysed these experimental data using multipole refinement of electron density according to the Hansen–Coppens model (Hansen and Coppens, Reference Hansen and Coppens1978; Coppens et al., Reference Coppens, Gururow, Leung, Stevens, Becker and Yang1979). We have also optimised the structure of α-AlOOH by DFT calculations using CRYSTAL17 software (Dovesi et al., Reference Dovesi, Orlando, Civalleri, Roetti, Saunders and Zicovich-Wilson2005, Reference Dovesi, Erba, Orlando, Zicovich-Wilson, Civalleri, Maschio, Rérat, Casassa, Baima, Salustro and Kirtman2018). Calculated on the basis of this optimisation, dynamical structure factors ($\tilde{F}_{hkl}$) (Erba et al., Reference Erba, Ferrabone, Orlando and Dovesi2013) have allowed us to obtain theoretical electron density distributions, which can be used as a benchmark for experimental results. Furthermore, to test the reusability of older X-ray data sets and their usefulness in more advanced modern quantum crystallographic approaches available today, we made an attempt to perform similar analysis using the ambient pressure dataset published by Hill (Hill, Reference Hill1979). Those data were used as an input for a full multipole model refinement, although the quality of those results turned out to be quite limited. As a comparison benchmark, we also conducted X-ray measurements on diaspore at ambient pressure with our laboratory diffractometer.
We also performed a multipole model refinement for δ-AlOOH, using a published set of experimental XRD data by Komatsu et al. (Reference Komatsu, Kuribayashi, Sano, Ohtani and Kudoh2006). The data from the crystallographic information file deposited with that paper was used as a starting point to conduct multipole model refinement in the same way as it was conducted for our experimental data for the α-AlOOH phase. In addition optimisation of the δ-AlOOH structure with CRYSTAL17 and calculations of dynamical structure factors ($\tilde{F}_{hkl}$) allowed us to obtain a theoretical density distribution. As a result, for both α and δ phases, we were able to obtain charge-density distributions on the basis of experimental data and theoretical calculations. Using these results, we attempted to compare the α and δ phases from the point of view of electron density distribution.
Describing a 3D structure with the use of polyhedra is just a convenient simplification. When one considers differences in electron density distribution, it is more beneficial to compare the shapes of atomic basins rather than the geometry of coordination polyhedra. Hence, in this study, we are using atomic basins rather than classic polyhedra to describe both investigated phases of AlOOH.
Experimental
Experimental quantitative charge-density studies
Aspherical modelling allows up to 32 additional parameters per atom, compared to routine refinements (Hansen and Coppens, Reference Hansen and Coppens1978; Coppens et al., Reference Coppens, Gururow, Leung, Stevens, Becker and Yang1979; Coppens and Coppens, Reference Coppens and Coppens1997). The independent atom model (IAM) refines only 9 parameters, i.e. positions and atomic displacement parameters (ADPs), related to atomic thermal vibrations. However, a high resolution is necessary to secure a sufficient observation/parameter ratio and to properly deconvolute the valence electron density from the thermal motions of atoms. High-resolution reflections are strongly associated with nuclear positions and with ADPs. The valence electron density contributes mostly to the low-angle diffraction reflections. The data set should be sufficiently complete to avoid systematic effects in the refinement – ideally 100% complete to a resolution higher that sinθ/λ = 1. Both accurate and precise measurements of the reflection intensities are needed to minimise systematic effects such as sample absorption, extinction and, very importantly, absorption by the apparatus itself, namely the diamonds and the metal gasket (Gajda et al., Reference Gajda, Stachowicz, Makal, Sutula, Parafiniuk, Fertey and Wozniak2020, Reference Gajda, Zhang, Parafiniuk, Dera and Wozniak2022; Stachowicz et al., Reference Stachowicz, Gajda, Huć, Parafiniuk, Makal, Sutuła, Fertey and Woźniak2023). Furthermore, the contribution of valence electrons to the total reflection intensity never exceeds a few percent, which reinforces the need for very accurate intensity measurements. Apart from accurate data correction, modern area-detector technologies offer the possibility to improve measurement precision by repeated collection of the same intensities (redundancy) (Sanjuan-Szklarz et al., Reference Sanjuan-Szklarz, Woinska, Domagala, Dominiak, Grabowsky, Jayatilaka, Gutmann and Wozniak2020).
The most common, aspherical quantitative model of experimental charge density is based on a finite spherical harmonic expansion of the electronic part of the charge distribution about each atomic centre, called a pseudo-atom. In the formalism of Hansen and Coppens (Hansen and Coppens, Reference Hansen and Coppens1978; Coppens et al., Reference Coppens, Gururow, Leung, Stevens, Becker and Yang1979) the pseudo-atom electron density is defined as:
where ρc(r) and ρv(r) are spherical core and valence densities, respectively. The third term contains the sum of the angular functions d lm±(θ,ϕ) to take into account the aspherical deformations. The angular functions d lm±(θ,ϕ) are real spherical harmonic functions. Coefficients P v and P lm± are populations for the valence and deformation density multipoles, respectively. κ and κ' are scaling parameters introduced to make the valence and deformation densities expand or contract. In the Hansen–Coppens formalism, P v, P lm±, κ and κ' are refinable parameters together with the atomic coordinates and thermal coefficients. Least-squares refinements were performed against the measured intensities of reflections (i.e. F 2(hkl)) obtained by single-crystal X-ray diffraction. This requires a data resolution of 0.50 Å and 100% data completeness. Starting atomic coordinates and anisotropic displacement parameters are taken from the ordinary spherical refinement stage and freely refined. Each atom was assigned core and spherical-valence scattering factors derived from the Su and Coppens (Reference Su and Coppens1997) wavefunctions. A single-ζ Slater-type radial function multiplied by density-normalised spherical harmonics was used for describing the valence deformation terms. The multipole expansion was truncated at the hexadecapole level for Al and O atoms.
Once an aspherical atomic electron density ρ(r) was defined, then it could be used to obtain aspherical atomic form factors and aspherical structural factors for a given crystal.
In AIM theory (Bader, Reference Bader1994; Popelier, Reference Popelier1996), a many electron system is separated into subsystems (ionic basins) by zero-flux surfaces (ZFSs). Any point on this surface satisfies the equation n•∇ρ(r) = 0, where ∇ρ(r) is the gradient vector field of the electron density, r is a point on the zero-flux surface that separates two fragments, and n is the vector normal to the surface at that point. Further analysis of the gradient vector field of electron density results in localisation of the extremes of the electron density by finding critical points (CP) at which ∇ρ(rCP) = 0. Particularly useful are bond critical points – the weakest points in bonds which define their properties. Integrating properties over ionic basins is one of the cornerstones of AIM theory because it yields valuable information such as the integrated charges, the volumes of atoms/ions, their energies, and electronic populations as well as higher multiple moments and polarisabilities (Angyan et al., Reference Angyan, Jansen, Loos, Hattig and Hess1994).
Space-group setting
On the basis of powder-diffraction results, the δ-AlOOH phase structure of diaspore was first proposed as space group P21nm (Suzuki et al., Reference Suzuki, Ohtani and Kamada2000). Three space groups were proposed as stable structures from first principle calculations Pnnm, P21nm and Pn21m (Tsuchiya et al., Reference Tsuchiya, Tsuchiya, Tsuneyuki and Yamanaka2002). Finally, a single-crystal X-ray data collection led to the P21nm space group (Komatsu et al., Reference Komatsu, Kuribayashi, Sano, Ohtani and Kudoh2006). Again, this is a non-standard setting of the space group. However in this paper, we describe α-ALOOH with the use of the standard setting Pnma. To be consistent, we used the standard setting of Pmn21 for the lower symmetry phase of diaspore. This setting is convenient to compare this structure with the Pnma structure. As described by Komatsu and co-workers (Komatsu et al., Reference Komatsu, Kuribayashi, Sano, Ohtani and Kudoh2006), a phase transition from the α to δ phase was possible at quite severe conditions, 18 GPa and 1273 K for 10 hrs. Generally, it is expected that HP causes shrinking of the unit cell and transformation to a crystal of higher symmetry. In our case the Z number changed from 4 in the alpha phase, to 2 in the delta phase. The equivalent volume/density difference between the two phases at ambient pressure is ca. 4.6%. The delta phase indeed has a lower volume and higher density.
Data quality
The primary goal of this paper is to compare electron density distribution in the crystal structure of the α and δ phase of AlOOH. As current experimental and literature data originated from different sources (different radiation sources) and were collected under different conditions (ambient or non-ambient pressure) the differences in the quality of these data sets are notable. This fact should be taken into consideration when comparing these results. The most basic parameters of particular data sets such as resolution, the number of measured reflections, total completeness, wavelength and pressure conditions are presented in Table 1.
The oldest literature data set, which was in fact the inspiration for the current investigation, was an attempt to refine the electron density distribution in diaspore conducted by Hill (Reference Hill1979). Although the resolution of this data set seems to be quite high, in fact, the completeness of data is rather poor. The completeness within the range typical for routine structural studies is close to 100% but varies between 50–70% for the high-resolution range. We have made an attempt to use this data set (the full list of reflections is available with the original publication) to refine the electron density according to the Hansen–Coppens multipole model (Hansen and Coppens, Reference Hansen and Coppens1978; Coppens et al., Reference Coppens, Gururow, Leung, Stevens, Becker and Yang1979). However, the results of refinement against the full hkl data set (791 reflections) are rather poor because the highest peaks and deepest holes in the residual density map are +1.78 and –1.95, respectively.
Our own experimental data set, which is the basis for the current considerations, was measured at the CRISTAL beamline of the SOLEIL synchrotron facility and was collected under high pressure. The experiment was conducted under pressure using a diamond-anvil cell, because we wanted to check if datasets collected in this way are suitable for obtaining the electron density distribution. If successful, further, more comprehensive, tests will be carried out. Ideally, the set of pressure points showing step by step how electron density changes as a function of pressure before and after a pressure triggered phase transition.
Although a single crystal was placed in the diamond-anvil cell, which causes a restriction of access to reciprocal space, due to the use of short-wavelength synchrotron X-ray radiation it was possible to collect X-ray diffraction data of a resolution comparable to the literature data (Hill, Reference Hill1979) with even better completeness. This data set will be referred to as the SOLEIL experimental data.
As a benchmark we have also collected a dataset of diaspore at ambient pressure on our laboratory diffractometer (λ = 0.56087). This data will be referred to as the Lab. Diffract. Experiment. These data are our only experimental data which are, in practise, 100% complete.
The 4th experimental data set considered in this paper is literature data for the δ-AlOOH phase (Komatsu et al., Reference Komatsu, Kuribayashi, Sano, Ohtani and Kudoh2006). The crystallographic information file attached to this manuscript contains a list of all the measured intensities of reflections, allowing an attempt to refine the electron density for this structure according to the Hansen–Coppens multipole model (Hansen and Coppens, Reference Hansen and Coppens1978; Coppens et al., Reference Coppens, Gururow, Leung, Stevens, Becker and Yang1979). This data set was collected at the synchrotron facility (in Table 1 named ‘δ-AlOOH experiment’). However, the beam wavelength was relatively long λ = 0.7007(1) Å. Only 262 reflections were measured which gave a completeness of only 62%. Such a poor resolution and completeness will affect the results of any refinement, especially of the multipole model of electron density refinement. Our analysis of the given list of reflections revealed that only a half of the potentially accessible reciprocal space was measured. Therefore to improve the results, we decided to support the experimental data with the theoretical.
As the first step of preparing the input data set for the multipole refinement based on theoretical calculations, the crystal structure of δ-AlOOH (reflections from Komatsu et al., Reference Komatsu, Kuribayashi, Sano, Ohtani and Kudoh2006) was optimised in CRYSTAL17 (Dovesi et al., Reference Dovesi, Orlando, Civalleri, Roetti, Saunders and Zicovich-Wilson2005, Reference Dovesi, Erba, Orlando, Zicovich-Wilson, Civalleri, Maschio, Rérat, Casassa, Baima, Salustro and Kirtman2018). We used the B3LYP (Lee et al., Reference Lee, Yang and Parr1988; Becke, Reference Becke1993) exchange–correlation function corrected for dispersion by Grimme's D3 (Grimme et al., Reference Grimme, Antony, Ehrlich and Krieg2010) correction in conjunction with the pob-TZVP-rev2 basis sets (Vilela Oliveira et al., Reference Vilela Oliveira, Laun, Peintinger and Bredow2019). Calculations were conducted for ambient pressure. Optimisation of the atomic position and cell volume was allowed. Convergence criterion on the root mean square of the gradient (TOLDEG = 0.00085), as well as on the displacement (TOLDEX = 0.0009), was also employed. The calculation grid had 75 radial points and maximum number of 974 angular points in the regions relevant for chemical bonding (XLGRID). The truncation criteria for bi-electronic integrals were adopted (TOLINTEG = 7 7 7 9 30). A Pack-Monkhorst/Gilat shrinking factor was also used (8 – shrinking factor in reciprocal space, 16 – shrinking factor for Gilat net). The convergence accelerator (BROYDEN) was also used. Fock matrices and KS matrices were mixed in a 50:50 ratio.
The dynamical structure factors ($\tilde{F}_{hkl}$) were calculated (Erba et al., Reference Erba, Ferrabone, Orlando and Dovesi2013). The list of structure factors was built in such a way that the resolution given in the experimental data was preserved and, additionally, not measured but potentially accessible reflections were added to the list. As a result, theoretical data which correspond to the experimental δ-AlOOH data set had a comparable resolution but completeness equal to 100% (‘δ-AlOOH CRYSTAL17’ in Table 1).
The theoretical data set corresponding to α-AlOOH (based on the SOLEIL experiment) was prepared similarly. In the case of the α-AlOOH data set, the theoretical list of reflections was cut slightly to obtain resolution and completeness comparable to the δ-AlOOH experimental data.
These prepared structure factors were used as the input data sets for multipole refinement which was carried out in the same way as the refinements of the experimental data sets. There are two major differences between the experimental and theoretical data sets which should be underlined here. Firstly, the experimental hkl file contains the experimentally measured reflection intensities (I(hkl) ≈ F(hkl)2) with their sample standard deviations. Secondly, whereas theoretical hkl files contain calculated dynamical structure factors ($\tilde{F}_{hkl}$) and because they are calculated, they have no sample standard deviations although they were used in the optimisation in the same form as the experimental values.
Data reduction
Data reduction for all the frames collected was performed using CrysAlisPRO software (CrysAlis Pro, 2014). Next, the structures were solved and refined with ShelXS (Sheldrick, Reference Sheldrick2008) and ShelXL (Sheldrick, Reference Sheldrick2015), respectively, within the Olex2 suite (Dolomanov et al., Reference Dolomanov, Bourhis, Gildea, Howard and Puschmann2009). Then, the intensities for each of the measurements were merged using Sortav (Blessing, Reference Blessing1995) implemented in the WinGX program suite (Farrugia, Reference Farrugia2012). The merged reflection intensity data sets were used subsequently as an input for the XD2016 (Volkov et al., Reference Volkov, Macchi, Farrugia, Gatti, Mallinson, Richter and Koritsanszky2016) program.
Parameters describing results of IAM refinement of only the experimental data (including literature data) are presented in Table 2. Table 3 presents results of multipole refinement conducted on the basis of experimental data as well as theoretical calculations. For all the datasets in Table 3 the same weighting scheme was used: w2 = 1/[s 2(F o2)]. Literature references indicate which original data were used as the starting point for refinement. The number of parameters determined in each refinement presented in Table 3 was slightly different. For α-AlOOH (theoretical data), it was possible to refine the κ and κ’ parameters for both Al and O atoms. For α-AlOOH (SOLEIL experiment), κ’ was not refined. For α-AlOOH (based on Hill's data), only κ for oxygen atoms was refined and the dipole parameter for hydrogen atom was not refined. In the case of both δ-AlOOH, the number of parameters was lower because for Al only the monopole was refined. The value of S (goodness of fit) for purely computational and for computationally augmented data was unreasonably high due to the fact that no uncertainty estimates for reflection intensities were available from the calculations, and standard deviations were arbitrarily set to unit values of sigma to fulfil the software requirements.
The CSD entries with deposition numbers 2298866–2298867 and 2307563–2307564 contain the supplementary crystallographic data for the α-AlOOH system investigated in this paper. This data can be obtained freely via http://www.ccdc.cam.ac.uk/data_request/cif, or by contacting [email protected] or the Cambridge Crystallographic Data Centre directly. The crystallographic information files have also been deposited with the Principal Editor of Mineralogical Magazine and are available as Supplementary material (see below)
Results and discussion
Phase comparison
The crystal structure of diaspore α-AlOOH is relatively simple. The crystal system is orthorhombic with four molecules within the unit cell. Two different mirror planes pass through each unit cell and all of the atoms are placed on these planes. As a result each atom is at a special position with the occupation factor equal to ½. In the asymmetric part of the unit cell, there are: one aluminium, one hydrogen atom, and two types of oxygen atoms (one of them is bonded with hydrogen and the other is not) (see Fig. 1). Each Al atom is surrounded by six oxygen atoms (three bonded with hydrogens and three not bonded). In fact, the most common structural element of the crystal structure is an AlO6 octahedron. Octahedra are arranged in a specific way (see Fig. 2).
In the α-AlOOH structure, the edge-sharing AlO6 octahedra are organised in chains along the y-direction. Within such a chain, the octahedra have the same orientation. Each chain of octahedra is connected with three other chains, with one of them sharing edges, and two others sharing vertices. Tunnels in the crystal structure along the z-axis are formed between the chains connected in this manner.
To simplify this description, we can distinguish four repeating types of chains within the diaspore structure (marked as a to d and depicted with different colours in Fig. 2. Some voids between polyhedra are clearly visible along the [010] direction and are filled by the O–H‧‧‧O hydrogen bonds, omitted for clarity in Fig. 2. At ambient temperature, the pressure of 18 GPa is not high enough to trigger the phase transition by itself, due to the high-activation barrier. The sample must also be heated for many hours at a temperature close to 1000°C, which illustrates that a significant amount of energy must be delivered to the system to break the existing connections and rearrange the AlO6 blocks.
In the δ phase, the atomic arrangement after phase transition seems to be very similar to that observed in the α-AlOOH phase. Four atoms in the asymmetric unit are located on a mirror plane. The Al cation at the centre of the octahedral site is coordinated by two types of oxygen atoms. In the δ-phase, the oxygen atoms are placed at the intersection of mirror and glide planes. As the oxygen atoms in the δ-phase have fewer degrees of freedom, it is expected that the deformation of octahedra could be smaller. When the structure is depicted in the form of AlO6 octahedra, it is easy to understand how it has changed in comparison with the α-AlOOH phase (see Fig. 2). The structure of the δ-phase is simplified because instead of four repeating chains of AlO6 octahedra in α-AlOOH there are only two of them.
In Fig. 3, one can clearly see that a schematic concept of the polyhedra does not necessarily correspond directly to the shapes and orientations of the atomic basins. First, in the polyhedral representation, ligand atoms are placed at the vertices of the polyhedra, whereas in the basins representation of the ligands are located inside the atomic basins. For that reason, all the atomic basins which belong to the central polyhedron in Fig. 3 are coloured orange. Second, the volumes of polyhedra have nothing in common with the volumes of corresponding atomic basins. In the case of atomic basins, the electron density which belongs to a particular atom is inside the 3D boundaries. For polyhedra, their boundaries have, in fact, no physical meaning.
Atomic basins
Atomic basins are a good representation of the distribution and partition of electron density. They can be compared within groups of corresponding to typical structural elements, for example, the AlO6 blocks in the structure of AlOOH. Examples of such aggregates for both phases are depicted in Fig. 4. To distinguish the existence of two types of oxygen atoms, the first one with bonds to hydrogens and the second one with the non-bonding interactions, particular basins are coloured in purple and orange, respectively. The same colouring scheme is used also in the case of the δ-AlOOH phase.
Arrangement of particular atomic basins is worth consideration in the context of their closest surroundings. Atomic basins describing the AlO6 octahedron in the crystal structure of the α-AlOOH phase are presented in Fig. 5. The α-phase is relatively densely packed. Atomic basins belonging to the oxygen atoms cover the Al basin completely from each side. The basins should fill in the space of the unit cell completely (taking hydrogens into account). The phase transition causes rearrangement of the building blocks (see Fig. 6). The higher pressure δ-AlOOH phase results in a lower symmetry structure. The oxygen atomic basins have much more irregular shape than those in the α-phase. The basins of the type 1 oxygen atoms (bonding with hydrogen atoms, purple shapes) are smaller and their shape is slimmer in comparison with bulky type 2 oxygen basins (orange shapes). Atomic basins describing the AlO6 group in the crystal structure of the δ-AlOOH phase are presented in Fig. 6.
As a consequence of rearrangement of electron density distribution (EDD), the atomic positions are affected and phase transition can occur. Such rearrangement of EDD can be visualised by maps of total electron density (see Fig. 7, second row). One can see that redistribution of EDD is significant and results in unquestionable phase transition. In δ-AlOOH, all atoms are placed at special positions on the mirror plane (the same situation as with α-AlOOH). That is why it is convenient to choose such a plane to compare molecular arrangements in both space groups.
At the first glance, it seems that transition from Pnma to Pnm21 results only in a twofold screw axes along the z-direction instead of the a glide plane perpendicular to the z-axis. In fact, we are losing more symmetry elements. As a result of the phase transition from the α-AlOOH to δ-AlOOH phase, the a glides perpendicular to the z-axis disappear as well as the twofold screw axes along the x and y axes.
However, such two-dimensional maps show only some small fragments of the structure and the 3D visualisations (i.e. atomic basins) could be much more convenient in some aspects).
Information about the fractional atomic coordinates and U eq values for both types of refined structures, α-AlOOH and δ-AlOOH are presented in Table 4.
Integrated atomic values
Integrated atomic charges and volumes are included in the Tables 5 and 6, respectively. The tables only contain results for α-AlOOH because basin integration for δ-AlOOH was not successful. One of the reasons could be the relatively poor data for δ-AlOOH.
* volume of the experimental unit cell; **volume of the unit cell calculated theoretically
In general, integrated atomic charges and volumes are important indicators of the quality of the multipole refinement. The refined molecule (or pair of anion and cation) should be neutral and the sum of all the atomic basins (within unit cell) together should be equal to the volume of the unit cell. Only minimal discrepancies can be allowed.
For charges, the Hill's data leads to relatively significant underestimation of the charge on Al(1) (ca. +1.47 e–) compared to the modern ambient experimental (2.29 e–) and theoretical data (2.57 e–), as well as the 2.5 GPa synchrotron results (2.33 e–). In addition, the literature charge value for O(2) is over two times smaller than that calculated for datasets collected in this work. However, literature data indicates a neutral molecule, whereas our calculations show small deviations from 0. Additionally, due to the specific calculation procedure our data are not accompanied by defined systematic error. Neither the experimental nor theoretical total ionic charges take formal values.
For the atomic volumes (Table 6), there was quite good agreement between the experimental atomic volumes calculated on the basis of the experimental data from SOLEIL and the results obtained by applying CRYSTAL17. The only exception was the atomic volume of H(1) which was far smaller when the experimental data from SOLEIL were used (1.12 Å3) compared to 3.19 Å3 for the theoretical data from CRYSTAL17. In each case the sum of all atomic basins (within unit cell) is very close to the volume of the whole unit cell.
Bond critical points
Selected parameters at bond critical points (BCP) such as electron density and Laplacian values are given in Table 7. In the case of bond-critical points, the spread of results was bigger than that of the integrated atomic values. It was not clear why the BCP was not found for the O(1)–H(1) bond in α-AlOOH (CRYSTAL17) or δ-AlOOH (experimental data). In our opinion it might have been due to the limited quality of the data.
Especially interesting was comparison of the Laplacian values. According to the theory of atoms in molecules (AIM) there are shared and closed-shell interactions (Bader, Reference Bader1994, Reference Bader2006; Zhang et al., Reference Zhang, Ying, Wu, Hiberty and Shaik2009). The shared interaction has a large value for charge density and a negative Laplacian value at the BCP, whereas the closed-shell interaction has a small charge-density value and a positive Laplacian value at the BCP. In this way, the charge density and Laplacian values at the BCP can characterise the nature of the bonding.
Positive Laplacian values usually indicates ionic bonding whereas negative Laplacian values indicate covalent bonding. In the case of positive Laplacian values there is another possible explanation: a charge-shift bond where the covalent term is repulsive rather than attractive. In Table 7 we see that O–H bonds are clearly shared interactions (covalent bonds), because they have relatively high charge-density values and negative Laplacian values, as expected. However, these results are no longer so clear in the case of the Al–O interactions. For both phases α and δ results correspond with synchrotron experiments and their theoretical benchmarks show positive Laplacian values (closed-shell interaction). The results of our experiments conducted on our laboratory diffractometer shows that the Laplacian values of the Al–O bonds are no longer clearly positive but rather gently negative. Taking into account that the dataset from our laboratory diffractometer is the only one which is 100% complete it raises the question how much such a sensitive parameter such as the Laplacian (second derivative of charge density) could be affected by low data completeness. We think that this might be such a case.
Al–O bonds and distortion of AlO6 octahedra.
Bond length values and polyhedral distortion parameters are in Tables 8 and 9, respectively. For the bond lengths, there was a good agreement between the final results from our laboratory experiments (diffractometer with Ag X-ray source) with Hill's results. The Al(1)–O(1) bond lengths in δ-AlOOH were slightly shorter than in α-AlOOH and the reverse was true for the Al(1)–O(2) bond length (see Table 8). In general, the bond lengths agreed very well. The largest differences were for the valence angles, for example, when Hill's data were compared with the SOLEIL data set, the differences were as large as 8.19(6)° for the O(1)–Al(1)–O(2) or O(1)–Al(1)–O(2) valence angles.
α: (i) –x, ½+y, 1–z; (ii) ½–x, 1–y, –½+z; δ: (i) 1–x, ½–y, ½+z; (ii) –x, ½–y, –½+z
Note – the symbols are described in the text.
The parameter d mean, used in Table 9, refers to the average Al–O distance in the AlO6 octahedron. Other parameters such as ζ and Σ describe the stretching and angular distortions, respectively. Parameters ζ, Δ and Σ are defined as follows:
Parameter ζ is the average of the sum of the deviation of 6 unique Al–O bond lengths around the central Al atom (d i) from the average value (d mean). Parameter Δ is the octahedral distortion parameter (Lufaso and Woodward, Reference Lufaso and Woodward2004). Parameter Σ is the sum of the deviation of 12 unique cis O–Al–O angles (ϕi) from 90°. All the mentioned parameters were calculated with use of OctaDist (Ketkaew et al., Reference Ketkaew, Tantirungrotechai, Harding, Chastanet, Guionneau, Marchivie and Harding2021).
Hydrogen bonds
Several scientific papers considered the issue of symmetrisation of hydrogen bonds in the structure of δ-AlOOH. (Sano-Furukawa et al., Reference Sano-Furukawa, Komatsu, Vanpeteghem and Ohtani2008, Reference Sano-Furukawa, Yagi, Okada, Gotou and Kikegawa2009, Reference Sano-Furukawa, Hattori, Komatsu, Kagi, Nagai, Molaison, dos Santos and Tulk2018; Cortona, Reference Cortona2017; Kang et al., Reference Kang, Feng, Yuan, Ye, Zhu, Huo, Li and Wu2017; Pillai et al., Reference Pillai, Jha, Padmalal, Maurya and Chamyal2018; Trybel et al., Reference Trybel, Meier, Wang and Steinle-Neumann2021; Luo et al., Reference Luo, Umemoto and Wentzcovitch2022) At ambient pressure, the hydrogen atom in the O(1)–H···O(2) hydrogen bond is located much closer to the O(1) donor atom than to the O(2) acceptor atom, and the potential energy surface for the hydrogen has a double well shape. However, pressure significantly changes this configuration. The energy barrier separating the two minima first becomes shallow enough that the hydrogen atom becomes disordered between two alternative sites, and eventually assumes a symmetrical position midway between O(1) and O(2), resulting in two subtle phase transitions. The structures investigated in our paper were measured (calculated) either at ambient conditions or at relatively low pressure, which is why such phenomena were not observed. The bond length between hydrogen and heavy atoms determined on the basis of X-ray measurements were usually shortened in comparison to neutron data. This is because of the shift of electron density in this bond towards the heavy atom and because of thermal motion of the H-atom which shortened the X–H bond length. The experimental data for α-AlOOH (Hill, Reference Hill1979) reported a 0.9886(8) Å O–H distance, whereas experimental data for δ-AlOOH (Komatsu et al., Reference Komatsu, Kuribayashi, Sano, Ohtani and Kudoh2006) reported an O–H distance equal to 0.81(4) Å. In the case of our multipolar refinement, when the data were incomplete, we decided to fix O–H bond length and not refine it. In the case of 100% complete data from a diffractometer (α-AlOOH Lab. Diffr. Experiment), the O–H bond was fully refined and determined as equal to 0.88(3) Å. This value was of course shorter than those obtained on the basis of neutron measurements but is typical for X-ray measurements.
Conclusions
Comparison of data sets which were collected under different circumstances and have different qualities is quite challenging, especially when a particular set of data has lower-than-expected completeness. However, even though some details seemed not to be determined precisely, it was still possible to focus on general trends. The α–δ phase transitions in AlOOH were triggered by high pressure and temperature, and caused significant reorganisation of the distribution of electron density. As a consequence, the crystal symmetry was lowered, significant reduction of the unit cell took place, and the basic AlO6 building blocks of the structure were rearranged.
This study demonstrates that when we compare currently collected datasets with those collected thirty or forty years ago we face problems. Issues caused by factors such as differences in hardware (detectors) and software (for example different treatment of weak reflections) quality. That is why straightforward comparison of older and current datasets is not always possible.
However, comparing the multiple data sets examined in this paper together (including historical data), we found common features, and showed that augmenting experimental data with computational results can help even incomplete data to reveal important information. Although the high-pressure high-resolution data of α-AlOOH are not 100% complete, they still clearly correspond with previously collected data and allow one to make suppositions about how phase transitions affect crystal structures.
Acknowledgements
Financial support for this work provided by the National Science Centre, Poland (OPUS grant No. UMO-2019/33/B/ST10/02671 awarded to KW) is gratefully acknowledged. This work was partly carried out at the Biological and Chemical Research Centre, University of Warsaw, established within the project co-financed by the European Union from the European Regional Development Fund under the Operational Programme Innovative Economy, 2007–2013. The work was accomplished at the TEAM TECH Core Facility for crystallographic and biophysical research to support the development of medicinal products sponsored by the Foundation for Polish Science (FNP). The synchrotron radiation experiments were performed at the SOLEIL (beamline CRISTAL) proposal No. 20210474.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1180/mgm.2024.22.
Competing interests
The authors declare none.