Open Access Article
Craig L.
Bull
*ab,
Christopher J.
Ridley
a,
Nicholas P.
Funnell
a,
Sumit
Konar‡
b and
James
Cumby
b
aISIS Neutron and Muon Source, Rutherford Appleton Laboratory, Didcot, Oxon OX11 7XN, UK. E-mail: craig.bull@stfc.ac.uk
bSchool of Chemistry, University of Edinburgh, David Brewster Road, Edinburgh EH9 3FJ, Scotland, UK
First published on 13th March 2024
The pressure and temperature phase diagram of the halide-based perovskite, CsCaCl3 is investigated using neutron diffraction. At ambient pressure, it undergoes a cubic to tetragonal phase transition at approximately 95 K. The structural evolution is characterised by changes in Cs and Ca coordination, bond lengths, and polyhedral volumes. Heat capacity measurements reveal a discontinuity at the phase transition, where we use a two-term Debye model to show the presence of two distinct contributions. Under pressure, CsCaCl3 exhibits phase transitions from cubic to tetragonal to rhombohedral symmetry, which can be understood when comparing the relative compressibilites of the CaCl6 and CsCln (where n = 12 for the cubic and 8 for the tetragonal and rhombohedral phases). Overall, this study provides experimentally derived thermodynamic properties and a tentative phase diagram of CsCaCl3 as a function of pressure and temperature.
m) with corner shared BX6 octahedra surrounding a AX12 polyhedron (see Fig. 1). Of particular interest to this manuscript is the case where X represents a halide (F, Cl, Br, I), rather than oxygen. Fernández et al.1 previously demonstrated that materials with the general formula ABCl3 crystallise in this perovskite structure, with a tolerance factor between 0.81 and 0.9.2,3 Halide-based perovskites, as with the oxide analogues, often demonstrate temperature and pressure-driven polymorphism, with crystallographic symmetry deviating from or converging towards the cubic aristotype via octahedral rotation/distortion, or A-site cation displacement. For example, cubic CsSrCl3 transforms to tetragonal then orthorhombic then monoclinic symmetry at 385 K, 375 K and 362 K respectively.1,4 CsPbCl3 is also highly polymorphic, with multiple phase transitions over a very narrow temperature range.5 The application of hydrostatic pressure as a tuning parameter is an efficient way of generating materials with new crystalline structures with intrinsic physical properties.6 Pressure tends to increase the electron density therefore enhancing the repulsive forces within a material, often resulting in unexpected behaviour. The halide based perovskites are good analogues for perovskite oxides under high-pressure conditions, with the advantage of showing greater compressibilities compared to their oxide counterparts.7 For example, the pressure driven isostructural phase transition of CsPbCl3 at room temperature.8
At 293 K, CsCaCl3 is cubic and has been shown to transform to tetragonal symmetry upon cooling below 95 K.9 It is an effective scintillator material, and its doped derivatives find applications in fields such as medical imaging and high-energy physics experiments.10–12 More generally, inorganic chloride perovskites also find applications in light-emitting diodes, photodetectors and field effect transistors.13–16 Understanding the structural and thermodynamic/physical property relationships of such materials is key to the design of novel materials with tailored physical properties which can prepared by chemical tuning. The application of pressure allows a clean and rapid determination of structure–property relationships which provides information for materials design suitable for devices with specific physical properties. In this study we have used neutron diffraction to study the structure of CsCaCl3 as a function of temperature and pressure. Isothermal compression of CsCaCl3 at 290 and 190 K reveals two structural phase transitions allowing the determination of a tentative phase diagram.
:
ethanol (4
:
1 by volume) as a pressure-transmitting medium.19 The gasket was sealed between two single-toroidal anvils, machined from zirconia toughened alumina, within a V3 Paris–Edinburgh press.20 Diffraction patterns were obtained in the transverse scattering geometry in approximately 5 tonne steps up to a maximum of 60 tonnes applied load, collected for ∼2 hours per step. The data were focused and corrected for attenuation using Mantid.21 The crystallographic structures were Rietveld-refined against the resulting diffraction patterns using GSAS.22 Fits to the unit-cell data and subsequent determination of equations-of-state were determined using PASCal.23
super-lattice reflection of the pseudo-cubic cell (corresponding to the R-point of the Brillouin zone). The R-point reflection is also accompanied by a splitting of the cubic 400 reflection at ∼2.69 Å which index as the 220 and 004 reflections in the new tetragonal cell (space group I4/mcm, a0a0c+ tilt system24). Our data collections are coarse in temperature steps however, we estimate that the transition from the cubic to tetragonal phase occurs at ∼95 K which is in good agreement with previous experimental results.9Fig. 3 shows the thermal evolution of the pseudo-cubic lattice parameters for the cubic (a = apc) and tetragonal (
, cpc = ctet/2) phases.
In the tetragonal phase the Cs coordination reduces from 12 to 8 (with an associated point group change from m
m to
2m), the CaCl6 polyhedra remain unchanged in coordination number, but there is a change in point group symmetry from m
m to 4/m. In both the cubic and tetragonal phases the derived structural parameters (polyhedral volumes and bond lengths) can be obtained following the parameterisation of Knight et al.,25 a full tabulation of which is provided in the ESI.† These are derived from the lattice parameter (a) in the cubic phase, and the (a, c) lattice parameters, and the order parameter (ϕ) in the tetragonal phase. ϕ is the rotation of the CaCl6 octahedra about [001]. This can be derived from the atomic coordinate of the Cl atom, u (where tan(ϕ) = 4u).26 The refined structural parameters are given in Table 1 for the two data sets shown in Fig. 2, see the ESI,† for further temperatures.
m; Wyckoff positions: Cs 1b, Ca 1a, Cl 3d. For the tetragonal structure the space-group is I4/mcm, Cs 4b, Ca 4c, Cl(1) 4a, Cl(2)
| Parameter | 5.0(1) K | 260.0(1) K |
|---|---|---|
| Symmetry | Tetragonal | Cubic |
| Space group | I4/mcm |
Pm m |
| a (Å) | 7.5880(3) | 5.3991(3) |
| b (Å) | 7.5880(3) | 5.3991(3) |
| c (Å) | 10.7953(5) | 5.3991(3) |
| Unit-cell volume (Å3) | 621.57(7) | 157.38(3) |
| Cl(2) u | 0.0233 | — |
| wRp, Rp (%) | 1.3, 1.5 | 1.46, 1.68 |
| χ 2 | 1.02 | 0.6 |
The temperature evolution of the tetragonal lattice parameters have been fitted to a model related to an Einstein internal energy function;27 the fits are shown as solid lines in Fig. 3. The fit function, and resulting fit parameters are provided in the ESI.† The symmetry-adapted tetragonal strain (detailed in the ESI,† and ref. 25), etz, tends towards 0 upon approaching the phase transition, as the cubic phase is by necessity free of tetragonal-strain. Following the work of Knight et al.,25 it is possible to model the temperature dependence of etz with the following equation:
![]() | (1) |
![]() | (2) |
The two distinct Ca–Cl bonds both decrease in length with increasing temperature, while in the cubic phase the single bond-length increases linearly as the temperature is increased, see Fig. 4. This is mirrored in the temperature dependence of the CaCl6 polyhedral volume. This behaviour is analogous to that seen in the structurally-similar system RbCaF3.26 In the tetragonal phase the axial Cs–Cl(1) bond shows only a small variation with temperature, which is in contrast to the shorter bonded Cs–Cl(2) bond which rapidly increases in length with increasing temperature and this behaviour is mirrored by the longer non-bonding Cs–Cl(2) bond which reduces in length. In the cubic phase these three distinct bonds become equivalent and increase in length with increasing temperature. Overall, this results in a non-linear increase in the CsCl8 polyhedral volume with temperature, becoming a linear increase for the CsCl12 polyhedra in the cubic phase. The octahedral tilting (ϕ) is also shown, tending to zero as the phase transition to the cubic phase is approached.
The heat capacity and volumetric data were fitted in this way in the tetragonal phase, where the bulk modulus has been assumed to be temperature invariant, fixed at the value determined at ambient temperature for the tetragonal phase (see Table 3). It was found that a single-term Debye model fitted the volumetric data well (giving θD = 267(8) K), but it could not account for the measured temperature dependence in the heat capacity, see Fig. 5. Instead, a two-term Debye model was used, representing separated cation and anion behaviour, significantly improving the fit. Despite the limited number of data points, the fits converged stably; the determined thermodynamic parameters from these fits are summarised in Table 2. The weighted average value for the Debye temperature is θwD = 304 K, which is higher than that predicted recently from ab initio methods (233.7 K,30 225.67 K31). However both of these studies assumed cubic symmetry for the material, which may explain the discrepancy.
| Parameter | |
|---|---|
| V 0 (Å3) | 155.4(3) |
| γ 1 | 6(2) |
| θ D1 (K) | 461(12) |
| γ 2 | 2.0(5) |
| θ D2 (K) | 118(3) |
| Mixing parameter, z | 0.543(14) |
At the next pressure point measured (4.8(1) GPa), the relative intensity of the 220/004 reflections change and the pseudo-cubic 111 reflection is seen to broaden (suggesting a peak splitting§), indicative of a transition to lower symmetry due to further tilting of the CaCl6 octahedral units. This is caused by the freezing of phonons at the R-point (anti-phase tilt) or M-point (in-phase tilt) in reciprocal space of the cubic Pm
m structure. We do not observe an M-point reflection in the data set at high-pressure, though if it is very weak it could be masked by the background of the pressure–cell. The high sensitivity of neutron diffraction to chlorine would therefore indicate that, if present, an M-point distortion must be extremely small. We proceed assuming only an R-point distortion is present.
The pattern can be fitted with a rhombohedral phase, space group R
c (a−a−a−) for which only the R-point reflections are present, with a small fraction of the tetragonal phase coexisting. An attempt to fit the data against a lower-symmetry model, for example with orthorhombic Imma (a−a−c0) symmetry, did not yield a satisfactory fit to the data. The pattern remains mixed phase up to the maximum pressure measured of 6.4(1) GPa. Upon recovery to ambient conditions the R-point reflection is lost, as is the splitting of the pseudo-cubic reflections, showing full reversibility.
Fig. 7 shows the variation in the unit-cell volume with pressure in each of the three phases described above. Also shown are the derived Birch–Murnaghan equations of state (EoS) for each phase; the fitted EoS parameters are listed in Table 3. The value of B0 for the cubic phase of 25.3(9) GPa is significantly lower than would be expected for oxygen-based inorganic perovskites for example CaTiO3 the bulk modulus is 177 GPa,32 and is not unexpected for halide based inorganic perovskites, for example, CsPbCl3 has a B0 of 45.6 GPa.8 Theoretical prediction of B0 for CsCaCl3 suggest values between 23–32 dependent on the model used.33,34 The value of B0 is also comparable to polymeric halide compounds for example AuI which has a bulk modulus of 18 GPa.35 The bulk modulus determined for the tetragonal phase is greater than that of the cubic phase, as would be expected from a denser phase. However, unusually, the higher-pressure rhombohedral phase has a lower bulk modulus than the cubic phase. The values of B0 determined for the tetragonal and rhombohedral phases are unreliable (and only shown for completeness) as a result of the limited number of data points used for each fit. However, Fig. 7 clearly shows that the transition from cubic to tetragonal phase is close to second order in nature, given the almost continuous behaviour of the unit-cell volume across the phase-boundary and the first order like behaviour of the tetragonal to rhombohedral phase where a clear discontinuity in unit-cell volume is observed in addition to phase co-existence.
| Parameter | Cubic | Tetragonala | Rhombohedrala |
|---|---|---|---|
| a The number of data points available significantly limits the validity of the values determined in the two high pressure phases, and are shown only for completeness. | |||
| Axial compressibility | |||
| a (TPa−1) | 8.98(11) | 6.5 | 9.3(2) |
| b (TPa−1) | — | 6.5 | 8.5(4) |
| c (TPa−1) | — | 3.9 | 8.4(2) |
| EoS order | 3rd | 2nd | 2nd |
| V 0 (Å3) | 157.14(14) | 617.1(19) | 328.4(2) |
| B 0 (GPa) | 25.3(9) | 33.7(10) | 17(2) |
| B′ | 5.9(6) | — | — |
For the determination of tetragonal strain (etz) in the tetragonal high pressure phase it is possible to use the equations described previously. However a0, can be obtained by extrapolation of the cubic-phase EoS. Fig. 8 shows the variation in etz and apc with pressure for the tetragonal phase.
The structure of a rhombohedral symmetry perovskite with the R
c space group, has the atoms sited on Wyckoff positions: Cs 2a, Ca 2b, Cl 6e
. Here, w is the displacement of the Cl atom from the ideal cubic position and gives rise to tilting of the octahedra about the rhombohedral [111] axis. A value of w = 0 is the idealised cubic structure.36 The magnitude of the octahedral tilt in the rhombohedral phase (ω) is given by
.37 In the rhombohedral phase, the error in lattice parameter determination makes meaningful interpretation of the strains difficult, and hence only the strain tensors e11, e22 and e33 have been calculated. As described by Carpenter et al.38 the definition is valid for the strain as the determined values of the pseudo-cubic cell angle
is less then 90.5° and hence, can be defined as:
![]() | (3) |
. The variation of apc and that of e11 with pressure is shown in Fig. 8. This shows that in both phases, the strains are deviating away from zero, and are becoming more strained/distorted with increasing pressure, however the magnitude of the strain is small.
Fig. 9 shows the variation in the Cs–Cl and Ca–Cl bond lengths in all three phases with increasing pressure. In the cubic phase both the Cs–Cl and Ca–Cl bond lengths are decreasing as per the decrease in unit-cell dimensions. In both the tetragonal and rhombohedral phase the Ca–Cl bonds are decreasing overall in length with increasing pressure. In the tetragonal phase it is clear that the octahedra are becoming more distorted as the two distinct Ca–Cl bonds deviate from each other in length with increasing pressure. The same behaviour is seen for the behaviour of the bonds forming the CsCln polyhedra in both phases with increasing pressure. The variation of ω and ϕ for the rhombohedral and tetragonal phases respectively with increasing pressure at 290 K are also shown in Fig. 9. In both the tetragonal and rhombohedral phases the magnitude of octahedra tilt is increasing with increasing pressure. In the rhombohedral phase the octahedra are undistorted, however can homogeneously flatten or elongate along their triad axis of the octahedra and has been previously described by the factor ζ (the octahedral strain, defined following Megaw in the ESI†). The variation of ζ of the CaCl6 octahedra for the rhombohedral phase is shown in Fig. 9 and shows that the octahedra strain increase in value with increasing pressure.
It has been shown that the relative compressibilities of the polyhedral cation sites in the perovskite structure give an indication as to how the octahedra rotate upon compression, either by increase in rotation, tending towards lower symmetry, or a reduction in rotation and an increase in symmetry. A model based upon the bond valence concept for each polyhedra has been developed that predicts the compressibility of the polyhedra.39 The value Mi is the total estimated variation of bond valence in a polyhedral site due to the change of average bond distance.40 The relative values of the values of Mi for the two polyhedra predicts the behaviour upon compression (see ESI†). For CsCaCl3 based upon our current data and reference values elsewhere41 we derive a MCs/MCa ratio of ∼0.4 and predict that the CsCl12 polyhedra will be more compressible than the CaCl6 octahedra and as a result the structure will become more distorted with increasing pressure as is observed experimentally in the current study. The value of the ratio is significantly lower than that observed for oxide perovskites but may well reflect the difference in compressibilties of the oxide and halide polyhedra. Such a low value of the Mi ratio has also been observed for other 1+
:
2+ halide based materials such as KCaF3MK/MCa ratio of ∼0.5 which predicts the structural deformation of KCaF3 should increase with pressure which is also observed experimentally.42
Footnotes |
| † Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4tc00592a |
| ‡ Present address: Joseph Banks Laboratories, School of Chemistry, University of Lincoln, Lincoln LN6 7DL, UK. |
§ This follows the assumption that the sample is still under hydrostatic pressure conditions, as evidenced by the fact that the other reflections in the pattern remain sharp, while the absolute pressure remains below the hydrostatic limit for methanol:ethanol (4 : 1 by volume). |
| This journal is © The Royal Society of Chemistry 2024 |