Introduction

A comprehensive review of the crystal chemistry of hydroxide perovskites having stoichiometries AB(OH)3 and A2BB′(OH)6 is given by Mitchell et al. (2017). The structure of hydroxide perovskite minerals and synthetic analogues (Fig. 1) consists of a framework of corner-linked B(OH)6 octahedra and a vacant A site. All oxygens in hydroxide perovskites are both H donors and acceptors. The octahedra can be occupied by homovalent cations giving the stoichiometry B(OH)3, where B = Fe3+, Ga3+, In3+, or heterovalent with BB′(OH)6 where B = Na, Ca, Mg, Zn, Cu2+, Fe2+, Mn2+, Nb5+, and B′ = Sn4+, Si, Ge, and Sb. Heterovalent varieties have an ordered arrangement of alternating B(OH)6 and B′(OH)6 octahedra and hydroxide perovskites with B2+B4+ and B1+B5+ combinations having been observed in nature.

Fig. 1
figure 1

Crystal structure of a tetragonal P42/n and b cubic Pn\(\bar{3}\) hydroxide perovskites relevant to the present study. Successive octahedra when viewed along the axes of the cubic aristotype have either the reverse or the same sense of rotation. a A reverse rotation along [001] for the tetragonal a+a+c tilt system of P42/n, and b the same sense for the cubic a+a+a+ tilt system. The two different tilt systems have significantly different hydrogen-bonding topologies. Examples of hydrogen-bonded O–H···O bridges are shown as dashed red lines. A single B–O–B′ angle is shown in green in a. c The acute angle formed by [001]-Mn–O(3) in tetrawickmanite MnSn(OH)6, which corresponds to [001]-Cu–O(apical) in mushistonite. The orange triangle is used to calculate the likely magnitude of the component of the Cu–O(apical) bond that is parallel to the c axis

The absence of an A cation leads to a high degree of rotation of octahedra as the framework contracts around the vacant A site. Furthermore, the tilting angles of octahedra, as referred to the pseudo-cubic axes of the untilted Pm\(\bar{3}\)m aristotype (Glazer 1972), have a very limited range of values of 14°–17° and are almost independent of composition and space group. For example, the Pn\(\bar{3}\) structures of CaSn(OH)6 (burtite), MnSn(OH)6 (wickmanite), and MgSn(OH)6 (schoenfliesite) have tilts of 14°–17°, and the P42/n structures of stottite (Welch and Wunder 2012) and tetrawickmanite (Lafuente et al. 2015) both have tilts of 17° about [001] and 15° about [100]/[010]. MgSi(OH)6 with monoclinic space group P21/n (Welch and Wunder, 2012) has tilts of 16° [110], 14° [1–10], and 16° [001]. The very limited range of tilt angles shown by diverse hydroxide perovskites strongly points to a major role for hydrogen bonding in these structures, i.e., the principal control on tilting is the strength of hydrogen bonds.

Hydrogen bonding is a key feature of these structures, and the connectivity of hydrogen-bonded networks correlates with tilt system (Mitchell et al. 2017). For example, hydroxide perovskites with cubic space groups Pn\(\bar{3}\) and Im\(\bar{3}\) (both having tilt system a+a+a+) have isolated squares of O–H···O linkages without extended connectivity. Tetragonal hydroxide perovskites having space group P42/n and the corresponding tilt system a+a+c, such as stottite FeGe(OH)6, have a combination of isolated squares and crankshafts (Welch and Wunder 2012). In tetragonal structures, all H atom sites of crankshafts are half-occupied, whereas the H site of the square is fully occupied (e.g., Lafuente et al. 2015; Welch and Kleppe 2016). All oxygens are both hydrogen-bond donors and acceptors. In the monoclinic synthetic high-pressure perovskite MgSi(OH)6 with tilt system aac+ (Wunder et al. 2011; Welch and Wunder 2012), the hydrogen-bonded topology comprises crankshafts in which the sites of all four non-equivalent H atoms are half-occupied, and there are interleaved isolated zigzag chains in which the remaining H atom site is fully occupied. Correlations between space group and tilt system for hydroxide perovskites are summarised by Mitchell et al. (2017).

Hydroxide perovskites have potential technological importance and a range of properties that includes proton conduction (Jena et al. 2004), flame retardancy (Zhang et al. 2007), and photocatalytic activity (Fu et al. 2009). Natural occurrences are very rare, due to their low temperatures of formation and their reaction to other Sn-bearing phases (Mitchell 2002; Nefedov et al. 1977). Consequently, experimental studies have made extensive use of synthetic analogues, which are easily prepared following the method of Strunz and Contag (1960), although the products are typically very fine-grained and preclude study by single-crystal XRD.

The absence of the A cation makes it possible to study compositional and compressional effects of an octahedral framework related to perovskite. Studies of the compressibility of hydroxide perovskites (Ross et al. 2002; Welch and Crichton 2002; Welch et al. 2005) have shown that the absence of A site cations results in much higher compressibilities compared with perovskites. Bulk moduli of hydroxide perovskites are typically in the range 70–80 GPa, compared with 150–200 GPa for perovskites, e.g., GdFeO3K0 = 182 GPa (Ross et al. 2004) and CaSnO3K0 = 163 GPa (Kung et al. 2001). Burtite–CaSn(OH)6 has an exceptionally low bulk modulus of 47 GPa (Welch and Crichton 2002).

Although several equation-of-state studies of hydroxide perovskites have been made, the mechanisms of compression are unknown in detail. The highly tilted frameworks of ambient hydroxide perovskites suggest that a change in compression mechanism could occur at modest pressures. This behaviour is unexplored in hydroxide perovskites at the atomic scale.

The objective of this investigation is to determine the effects of composition and pressure upon the structural state of the mushistonite–vismirnovite solid solution series, CuxZn1−xSn(OH)6. Structure studies indicate that vismirnovite-ZnSn(OH)6 is cubic with SG Pn\(\bar{3}\) (Cohen-Addad 1968), and mushistonite–CuSn(OH)6, is tetragonal P42/nnm (Morgernstern-Badarau 1976). The Pn\(\bar{3}\) structure has tilt system a+a+a+ (Glazer 1972) and is consistent with other cubic members of the schoenfliesite group, such as burtite, wickmanite, and schoenfliesite (Basciano and Peterson 1998). However, space group P42/nnm was shown to be very implausible for mushistonite by Welch and Kleppe (2016), as it involves a tilt system (a0b+b+) with a zero tilt and a mirror plane bisecting the octahedra. Both features are improbable, because tilting of octahedra is expected for hydroxide perovskites, due to contraction around the empty A site and strong hydrogen bonding leads to further octahedral rotation. Furthermore, the topology of P42/nnm forces octahedra to become extremely distorted due to the bisecting mirror plane. Recent studies demonstrate that tetragonal hydroxide perovskites crystallize in space group P42/n (Ross et al. 1988; Lafuente et al. 2015; Welch and Kleppe 2016) with tilt system a+a+c.

Unfortunately, crystals of natural or synthetic mushistonite of a quality suitable for structure determination by single-crystal XRD have not been available. Consequently, powder-diffraction data for synthetic analogues have to be tested against various space groups. Peak-splitting associated with lowering of metric symmetry from cubic to tetragonal is usually easily detected, although the correct choice of tetragonal space group from the powder-diffraction data is more challenging. However, the only tetragonal space group known for double hydroxide perovskites is P42/n, which is, therefore, the clear candidate for mushistonite.

Mushistonite is unique among tetragonal hydroxide perovskites in having a “prolate” cell metric with c > a. Tetrawickmanite, stottite, and Ga(OH)3 (soehngeite) have an “oblate” cell metric (a > c). Furthermore, the “eccentricity” of the metric (c:a) is 1.07 for mushistonite, whereas the value for the three other tetragonal phases is only 0.99. The strongly prolate nature of the mushistonite unit cell points to a significant contribution from Jahn–Teller distortion of the Cu(OH)6 octahedron. Such a distortion in mushistonite is expected as it is a common feature of compounds with octahedrally coordinated transition metals having an odd number of d-electrons, such as d9 Cu2+. Consideration of the tetragonal (P42/n) structures of MnSn(OH)6, FeGe(OH)6, and Ga(OH)3 indicates that the acute angle between [001] and the B−O(3) bond is 20°–22° for these structures (Fig. 1c). A simple trigonometric calculation shows that the component of the B−O(3) bond that is parallel to [001] is 92% of the bond length for these structures. The prolate cell metric of CuSn(OH)6 suggests that the long Cu−O(apical) bond corresponds to that forming the [001]−B−O(3) angle, and so is also likely to have a component parallel to [001] of at least 90%. Hence, the behaviour of the c parameter of CuSn(OH)6 is likely to be highly correlated with that of the Cu−O(apical) bond of the Cu(OH)6 octahedron.

Methods

Sample preparation

The samples of CuxZn1−xSn(OH)6 with x = 0, 0.04, 0.1, 0.2, …, 0.9 and 1 were synthesised using the method of Strunz and Contag (1960). The reagents were 5 M aqueous solutions of K2SnO3.3H2O, CuCl2, and ZnCl2 (Sigma-Aldrich). The synthesis reaction involved is as follows:

$$\begin{aligned} &{\text{K}}_{ 2} {\text{SnO}}_{ 3} \cdot 3 {\text{H}}_{ 2} {\text{O}}_{{({\text{aq}})}} + x{\text{CuCl}}_{{ 2({\text{aq}})}} + \, ( 1- x){\text{ZnCl}}_{{ 2({\text{aq}})}} \\ &\quad\to {\text{Cu}}_{x} {\text{Zn}}_{ 1- x} {\text{Sn}}\left( {\text{OH}} \right)_{{ 6({\text{s}})}} + {\text{ 2 KCl}}_{{({\text{s}})}} . \end{aligned}$$

Solutions were mixed with magnetic stirrer at room temperature. Hydroxide perovskite was precipitated immediately. The well-crystallized precipitate was then washed repeatedly until all crystallized KCl was gone, as confirmed by powder-XRD.

Powder X-ray diffraction at ambient conditions

The powdered samples were loaded into circular well mount holders. Powder-diffraction data were collected in Bragg–Brentano geometry using a Panalytical X’Pert Pro MPD equipped with primary Ge(111) monochromator, X’celerator position sensitive detector (opening angle of 2.122°), incident and diffracted-beam Soller slits (0.02 rad), and ¼° divergence and ½° anti-scatter slits. X-ray tube operating conditions were 45 kV and 40 mA. Measurements were performed with Cu Kα1 radiation in the range from 5 to 125° 2θ, at a step size 0.02°, and an effective counting time of 100 s per step.

Synchrotron powder X-ray diffraction at high-pressure conditions

High-pressure measurements on cubic Cu0.4Zn0.6Sn(OH)6 and tetragonal CuSn(OH)6 were obtained using a diamond-anvil cell (DAC) fitted with rhenium gasket pre-indented to a thickness of 40 μm and drilled with a 150 μm sample hole. Liquid neon was used as the pressure-transmitting medium (Hemley et al. 1989), and the maximum pressure used was 17 GPa, which is its quasi-hydrostatic limit. Powdered samples of two key compositions, Cu0.4Zn0.6Sn(OH)6 (cubic, Pn\(\bar{3}\)) and CuSn(OH)6 (tetragonal, P42/n), were loaded into the DAC with a ruby sphere for pressure calibration using the ruby-fluorescence method (Mao et al. 1986). Diffraction experiments were collected on beamline I15 at the Diamond Light Source (UK) from 10−4 to 30 GPa at 298 K. Experiments were performed with a monochromatic beam (λ = 0.30996 Å) collimated to 30 μm in diameter. A MAR345 image plate detector (MarResearch) was used with a sample-to-detector distance of 450 mm and exposure times of 15–25 min. The sample-to-detector distance was calibrated with silicon powder and LaB6. The powder-diffraction images have been integrated with Fit2d.

Rietveld refinement

All XRD powder data were refined using GSAS (Larson and Von Dreele 2004), and refined parameters were background, peak profiles, unit cell, oxygen position, and isotropic displacement factors (Uiso). The B and B′ sites lie at 000 and ½ ½ ½, respectively, and so only their Uiso values were refined. Starting models for Rietveld refinements were the cubic Pn\(\bar{3}\) structure of vismirnovite (Cohen-Addad 1968), and the tetragonal P42/n structure of tetrawickmanite (Lafuente et al. 2015). The latter structure model was preferred over that reported for mushistonite by Morgernstern-Badarau (1976) for which the Sn octahedron shows a very unreasonable distortion due to an incorrect choice of space group (P42/nnm). Low intense diffraction peaks from the rhenium gasket and the neon pressure transmission medium were observed in some experiments. Here, structures for rhenium (Swanson and Fuyat 1953; Duffy et al. 1999) and neon (Hemley et al. 1989) were included in the refinements.

Results

XRD at ambient conditions

The solid products of all experiments were hydroxide perovskite and KCl. The latter was completely removed by washing as described above. Figure 2 shows XRD powder patterns for the Cu–Zn series. One-phase fields of cubic and tetragonal CuxZn1−xSn(OH)6 occur within the compositional ranges XCu = 0−0.4 and 0.9−1, respectively. Bulk compositions 0.5 ≤ XCu ≤ 0.8 produced coexisting cubic and tetragonal phases.

Fig. 2
figure 2

Powder-XRD patterns of CuxZn1-xSn(OH)6. Synthesis products were cubic within the compositional range XCu = 0−0.4 and tetragonal in the range 0.9−1. Bulk compositions 0.5 ≤ XCu ≤ 0.8 produced coexisting cubic and tetragonal phases

Variations of lattice parameters with composition are shown in Fig. 3. A very minor change can be observed for the cubic phase with increasing Cu content (0 ≤ XCu ≤ 0.4). However, for XCu > 0.4, there is a significant increase in metric anisotropy for the tetragonal phase with cell parameter atetr decreasing and ctetr increasing. Bond length and bond angles of the cubic phase are shown in Fig. 4 and listed in Table 1. The most significant changes from XCu = 0 to 0.4 are an increase in Cu/Zn–O bond lengths and a decrease in bond angles B–O–B′ (B = Cu, Zn, B′ = Sn). No bond-length data are presented for the tetragonal phase due to low precision of refined oxygen positions; Uiso values had to be refined as an overall value.

Fig. 3
figure 3

Lattice parameters of CuxZn1-xSn(OH)6 solid solutions

Fig. 4
figure 4

Bond length and bond angles of cubic CuxZn1-xSn(OH)6

Table 1 Structural parameters of cubic CuxZn1−xSn(OH)6 refined in Pn\(\bar{3}\)

High-pressure behaviour of Cu0.4Zn0.6Sn(OH)6 and CuSn(OH)6

Crystal structures of Cu0.4Zn0.6Sn(OH)6 and CuSn(OH)6 from laboratory XRD refinements were taken as starting models for Rietveld refinement using the GSAS program. Peak profiles, unit cell parameters, oxygen coordinates, and isotropic displacement factors were refined. A fixed Uiso of 0.004 Å2 for the Cu/Zn site (the refined value for the ambient out-of-cell data set) had to be used for the cubic structure due to a persistent non-positive-definite value being obtained otherwise. Furthermore, an overall Uiso value for oxygen had to be refined for the CuSn(OH)6 due to Uiso for two of the three oxygens becoming non-positive-definite. The constraints on oxygen displacement parameters lead to some uncertainty in the refined positions of these atoms, and so no bond-length data are presented for CuSn(OH)6.

Representative high-pressure powder-diffraction patterns for the two phases are shown in Fig. 5, and typical results of Rietveld refinement of high-pressure patterns for both phases are shown in Fig. 6. Unit cell parameters for the high-pressure experiments are given in Table 2. Normalised volume and axial data as functions of pressure are shown in Figs. 7 and 8, respectively. Structural data are given in Tables 3 and 4.

Fig. 5
figure 5

Representative high-pressure PXRD patterns of a cubic Cu0.4Zn0.6Sn(OH)6 and b tetragonal CuSn(OH)6. ↓ sign indicates decompression data

Fig. 6
figure 6

Rietveld refinement of a cubic Cu0.4Zn0.6Sn(OH)6 at 13.8 GPa and b tetragonal CuSn(OH)6 at 14.1 GPa

Table 2 Unit cell parameters of Cu0.4Zn0.6Sn(OH)6 and CuSn(OH)6 with pressure
Fig. 7
figure 7

Normalised plots of the variation of the unit cell volume with pressure for Cu0.4Zn0.6Sn(OH)6 and CuSn(OH)6. Filled circles are data collected on compression and open symbols for data collected on decompression. Lines are EOSBM2 fits to compression data. CuSn(OH)6 shows a change in compression behaviour above 7 GPa. The dashed black line is the PV curve of burtite–CaSn(OH)6 (Welch and Crichton 2002) and the dotted black line is the PV curve for vismirnovite-ZnSn(OH)6 from Welch and Crichton (unpublished data)

Fig. 8
figure 8

Normalised plots of the variation of unit cell axes with pressure. Compression data have filled symbols and decompression data have open symbols. Lines are EOSBM2 fits to the data for CuSn(OH)6

Table 3 Cu0.4Zn0.6Sn(OH)6 with pressure: structural and statistical parameters of the Rietveld refinement in Pn\(\bar{3}\)
Table 4 Bond lengths and angles with pressure of Cu0.4Zn0.6Sn(OH)6

On compression, the PV data for the cubic Cu0.4Zn0.6Sn(OH)6 decrease smoothly to 17 GPa, with no indication of an inflexion or discontinuity. A fit to a second-order Birch–Murnaghan equation-of-state (BMEOS) for the compression data gives K0 = 75.8(4) GPa. A very similar fit was obtained to third-order, which gives K0 = 75(1) GPa and K0′ = 4.0(2), indicating that a second-order model describes the PV data sufficiently. The decompression data for Cu0.4Zn0.6Sn(OH)6 lie on the compression trend.

Variations with pressure of bond angles and bond lengths of the cubic phase are shown in Figs. 9 and 10, respectively. The (Zn,Cu)–O–Sn bond angle decreases markedly on compression to 7 GPa and remains approximately constant at 131° between 8 and 17 GPa. Within the small standard errors on their values, the lengths of (Zn,Cu)–O and Sn–O bonds decrease linearly to 17 GPa, although there is a suggestion of a steepening of the trend for the (Zn,Cu)–O bond at ~ 8 GPa.

Fig. 9
figure 9

Variations of (Zn,Cu)–O and Sn–O bond lengths of Cu0.4Zn0.6Sn(OH)6 with pressure

Fig. 10
figure 10

Variation of the (Zn,Cu)–O–Sn angle and volumes of octahedra of Cu0.4Zn0.6Sn(OH)6 with pressure

Above 7 GPa, the compression curve of CuSn(OH)6 cannot be fitted by a single equation-of-state. The non-conformity between the trends below and above 7 GPa is seen as a slight steepening of the compression curve from 10 to 17 GPa. The axial compression curves also register this discontinuity. Furthermore, the discontinuity in the a parameter appears to show a change between two linear trends in which there is decrease in the compressibility at 10–17 GPa (Fig. 8). In contrast, the trend of the c parameter at 10–17 GPa shows a steepening that implies an increased compressibility.

Normalised PV data to 7 GPa for CuSn(OH)6 were fitted to a second-order BMEOS (Fig. 7), which gives K0 = 59.7(9) GPa. A third-order fit gives a similar value of K0 = 59(2) GPa with K′ = 4.6(3). Second-order fits of the axial data to 7 GPa (Fig. 8) give Ka0 = 79(2) GPa and Kc0 = 38.0(3) GPa. However, the fit to the a parameter data lies slightly below the data points, which vary almost linearly over this pressure range. Figures 7, 8 also show that the decompression curves of CuSn(OH)6 are very different from that obtained on compression, and that this effect is primarily due to the behaviour of the c parameter. This clear difference between compression and decompression curves implies that there is hysteresis in CuSn(OH)6 that is not seen in the cubic Cu0.4Zn0.6Sn(OH)6 phase. The decompression data for the c parameter fall well below the corresponding compression curve. Furthermore, the value of the ambient c parameter recovered on decompression is 7.989(6) Å compared with the initial value of 8.1010(7) Å.

Discussion

Structural trends of CuxZn1−xSn(OH)6 solid solutions at ambient conditions

Incorporation of Cu into CuxZn1−xSn(OH)6 hydroxide perovskite in the range 0 ≤ XCu ≤ 0.4 results in only a minor change of the cubic cell parameter. There is a decrease in B–O–B′ angle and an increase in B–O bond lengths (B = Cu, Zn, B′ = Sn) (Fig. 4; Table 1). At the local level, increasing distortion of the octahedral framework is expected with increasing Cu content due to Jahn–Teller distortion associated with Cu(OH)6 octahedra. For compositions 0 ≤ XCu ≤ 0.4, an increasing proportion of Cu(OH)6 octahedra results in an increase of the average B–O bond length, while long-range cubic symmetry is retained. Structure expansion by lengthening of B–O bonds is compensated by increasing tilting of neighbouring octahedra (reflected in a decrease of the B–O–B′ angle), leading to only a minor change in cell dimensions.

Higher Cu contents lead to the formation of CuSnCuSn quartets being more probable. It is conceivable that strain associated with aggregated Cu-rich domains (propagated Jahn–Teller distortion) results in critical behaviour and a change in overall symmetry Pn\(\bar{3}\) → P42/n for compositions with XCu > 0.4. Although precise bond-length data are unavailable for the tetragonal phase, it is reasonable to assume that the increase of unit cell parameter c is most likely correlated with increasing distortion of the long apical B–O(3) bond, which is the major connector for Cu(OH)6 octahedra along the [001] direction.

Structural trends of Cu0.4Zn0.6Sn(OH)6 and CuSn(OH)6 as functions of pressure

Normalised PV curves of cubic Cu0.4Zn0.6Sn(OH)6 and tetragonal CuSn(OH)6 are clearly separated (Fig. 7), with the latter being significantly more compressible. The compression curve of Cu0.4Zn0.6Sn(OH)6 decreases smoothly without inflexion or discontinuity, and compares closely with those curves of γ-Fe(OH)3 (Welch et al. 2005), FeGe(OH)6 (Ross et al. 2002), In(OH)3, MgSn(OH)6, and ZnSn(OH)6 (Welch, unpublished data) given in Table 5. The similarity is exemplified in Fig. 7 by the compression curve for synthetic ZnSn(OH)6 [K0 = 73(1) GPa, K0′ = 4]. Only burtite–CaSn(OH)6 has a significantly different bulk modulus [K0 = 47.4(4) GPa, K0′ = 4, Welch and Crichton 2002]. With the single exception of burtite, the similar bulk moduli (K0 = 72–80 GPa) for a wide range of hydroxide perovskites point to a compression behaviour that is largely independent of composition, symmetry, or tilt system.

Table 5 Bulk moduli of hydroxide perovskites

The bulk modulus of the tetragonal CuSn(OH)6 is much lower than that of cubic Cu0.4Zn0.6Sn(OH)6, and there is a strong compressional axial anisotropy: Kc0 = 38.0(3) GPa and Ka0 = 79(2) GPa. Furthermore, the axial modulus Kc0 of CuSn(OH)6 is even lower than K0 of burtite. At the atomic level, it would appear reasonable to attribute this strong compressional axial anisotropy to the presence of a long Cu–O bond that has a significant component (> 90%) parallel to [001], and arises from Jahn–Teller distortion of the Cu(OH)6 octahedron.

We interpret the trends in Cu/Zn–O and Sn–O bond lengths (Fig. 9), and Cu/Zn–O–Sn angle and octahedral volumes (Fig. 10) for Cu0.4Zn0.6Sn(OH)6 in terms of octahedral tilting as being the major mechanism to 7 GPa, but becoming less significant at higher pressures, whereas compression of octahedra continues almost linearly to 17 GPa. Furthermore, it is clear from the ambient and high-pressure studies presented here that Jahn–Teller distortion associated with Cu has a significant role in the crystal chemistry and physical behaviour of hydroxide perovskites.