- Research article
- Open Access
Characterisation of temperature-dependent phase transitions in 2,2-trimethylenedioxy- 4,4,6,6-tetrachlorocyclotriphosphazene, N3P3Cl4[O(CH2)3O]
© Coles et al 2007
- Received: 02 April 2007
- Accepted: 18 July 2007
- Published: 18 July 2007
The crystal structure of 2,2-trimethylenedioxy-4,4,6,6-tetrachlorocyclo triphosphazene has been determined at 120, 274 and 293 K. The result at 293 K confirms the room temperature Cmc21 structure, but at the lower temperatures the space group is Pna21. Nevertheless the basic structure remains the same, with only small displacements of the atoms, amounting to an average of 25 pm between 120 and 293 K.
X-ray diffraction and DSC results indicate that the phase transition takes place in two steps between 274 – 293 K and provides an understanding of previous NQR results. In the intermediate temperature range the molecules are displaced from their room temperature positions in such a way as to give an average structure with Cmc21 symmetry.
The overall phase transition is consistent with the occurrence of a soft lattice mode at room temperature in which a large displacement of the molecule in the x-direction is coupled with a flexing motion about an axis defined by the nitrogen atoms in the N1 and N3 positions.
- Differential Scanning Calorimetry
- Wyckoff Position
- Room Temperature Structure
- Variable Temperature Experiment
In order to maintain consistency with previous work on compounds 1 – 3 , we have used the original 2,4,6 numbering of the phosphorus atoms for the cyclophosphazene ring rather than change to the 1,3,5 numbering using IUPAC convention.
2.1. Synthesis and crystallisation
Compound 2 was prepared as in the literature , characterised by 1H and 31P NMR  and was crystallised from dichloromethane. The melting temperature of compound 2 in the literature  is given as 429–430 K but is measured as 427 K by Differential Scanning Calorimetry (DSC) in this work.
2.2 Differential Scanning Calorimetry
Thermal measurements were made using a Shimadzu Differential Scanning Calorimeter (DSC60). Samples were examined under a polarising Olympus microscope at 100× and good quality crystals were selected for the measurements. The selected crystals were then placed in custom made Al micro-crucibles and weighed on a Sartorius electronic microbalance. Since attention was paid to crystal quality the sample size used was small (0.326 mg). A second sample (1.982 mg) was also tested. Cooling to 323 K was performed using liquid nitrogen and samples were heated to room temperature at 5 K/min. A total number of 5 cycles were performed of cooling and heating. After the 5th cycle samples were melted by heating at 5 K/min to 438 K and then cooled to liquid nitrogen temperatures and then re-heated to room temperature. Nitrogen was used as purge gas with a flow rate of 60 cm3/min. The reference crucible was an empty Al micro-crucible.
2.3. Structure determinations
Intensity data were recorded on a Nonius KappaCCD diffractometer situated at the window of Bruker Nonius FR591 rotating anode generator equipped with a molybdenum target (λ Mo-kα = 0.71073 Å) and driven by COLLECT  and DENZO  software. Structures were determined using the direct methods procedure in SHELXS-97 and refined by full-matrix least squares on F2 using SHELXL-97 . Data were corrected for absorption effects by means of comparison of equivalent reflections using the program SADABS . Non-hydrogen atoms were refined anisotropically, whilst hydrogen atoms were fixed in idealised positions (0.99 Å) with their displacement parameters riding on the values of their parent atoms (1.2 U eq ). The temperature of the sample was controlled using an Oxford Cryosystems 700 series Cryostream Plus cryostat . Pertinent data collection and refinement parameters are openly available for download for the 120 K  (for CIF data see Additional File 1), 274 K (for CIF data see Additional File 2) and 293 K  structures (for CIF data see Additional File 3). CIF data has also been deposited with the Cambridge Crystallographic Data Centre (deposition numbers CCDC641683 – CCDC641685).
Although a number of diffractometer manufacturers provide rudimentary control of low temperature devices via their data collection software, it is generally not flexible enough to carry out complex variable temperature experiments in an automated fashion. To address the problems posed by the chemistry under investigation in this study, it was necessary to write software that could automate the collection of multiple data sets under various temperature-change regimes. The application is written in Python and utilizes the Nonius Python libraries of the COLLECT  software, along with the associated programmer's manual. The libraries include interfaces to both the KappaCCD and the Oxford Cryosystems cryostream enabling a synergistic approach to the experiment. The software allows the following variables to be controlled; start temperature, final temperature, increment (determines the number of steps), ramp rate and dwell time (to allow a crystal to equilibrate at its new temperature). The type of diffraction experiment can be set as either full data collection or unit cell determination, and the usual choices of experimental parameters are available. In order to avoid time consuming processing of multiple data sets, data reduction is carried out automatically and therefore, once the experiment has been designed, it can be left to run with no further human input.
Phase transition in compound 2
Although there is a very simple relationship between the low and room temperature structures of this compound, with one structure realisable from the other with only very small displacements of the atoms, the phase transition does not appear to follow the simple second-order model of Landau but is more complex, taking place in at least two stages over a temperature range of about fifteen degrees. The phenomenon of a two-stage transition occurring for such a simple symmetry relationship has been recognised previously . In such cases it is believed that there is an intermediate phase, in which there is a long range order corresponding to the high-temperature structure, but the local symmetry conforms to the low-temperature phase, with fluctuations occurring.
In order to compare the low temperature structure to the room temperature structure in Cmc21, it must first be transformed to a non-standard setting of space group Pna21, namely Pbn21, by interchanging the a- and b-axes. Then, using International Tables  we find these two space groups have a simple maximal index-2 group-subgroup relationship, which is class-equivalent rather than translation-equivalent, and so the primitive unit cell must double in volume in going from room temperature to low temperature. (However, since the room temperature space group is C-centred, the conventional cell remains the same.) This symmetry relationship is borne out by the actual structures, since the low-temperature structure is realisable from the other with only small displacements of the atoms, during which the asymmetric unit doubles and the molecular mirror plane disappears. In terms of the Wyckoff positions, atoms occupying the special positions (4a) in Cmc21 with fixed x-coordinate occupy general positions (4a) in Pbn21 with no such restriction, while atoms in the general positions (8b) of Cmc21 split into two separate sets of general positions (4a) in Pbn21.
If we apply the phenomenological theory of such transitions , we can recognise two order parameters which might contribute to the Gibbs energy of the system and lead to this behaviour. Firstly we note that the displacements leading to the low-temperature structure arise in an ordered manner, such that molecules separated by the vector (1/2, 1/2, 0), related by a lattice translation at room temperature, are displaced in opposite directions in the low temperature structure, where they are related only by a glide operation. Such an effect can arise in a "soft mode" transition, so-called  since it can be explained in terms of a temperature dependence of the frequency of a normal mode of crystal vibration. At some temperature there ceases to be a restoring force associated with that particular set of atomic displacements, since the potential well changes from having a symmetric minimum corresponding to the atomic positions at room temperature to having a symmetric double minimum. When the resulting instability becomes sufficiently pronounced, the crystal transforms co-operatively to one of the two possible lower symmetry states with equivalent structures.
Since the group-subgroup relationship is not to a translation-equivalent subgroup, we can locate where the propagation vector of this mode lies in reciprocal space away from the Brillioun zone origin, and so correlate it with the actual atomic displacements. Specifically, since the additional Bragg reflections of the low-temperature phase occur for (h + k) odd, the displacements of pairs of atoms (1/2, 1/2, 0) apart in the unit cell must be 180 degrees out-of-phase in the normal mode concerned. Given the loss of the mirror plane present in the high-temperature structure, we expect them to lie perpendicular to this plane, i.e. along the a-axis.
An attempt was made to describe the displacements in the low temperature phase in terms of the symmetry mode analysis of Aroyo and Perez-Mato . In order to apply this analysis to a molecule, we could consider it as a rigid body, in which case it will give rise to just six lattice modes. However this makes the assumption that any internal deformations are so stiff as to have insignificant amplitude as compared to the lattice modes. This is not necessarily the same as saying that the free molecule can adopt only one or several conformations; for, as we shall see in this case, a moiety normally considered rigid may show significant deformation, or one capable of adopting several conformations remain unaffected through the transition.
Now three of these lattice modes involve displacements of the centre of mass, which correlate with the translations of the free molecule. These can be treated in an identical way as for an individual atom by the Aroyo-Perez-Mato method. However the other three are librational modes, which correspond to rotations of the free molecule. Taking as an example a libration L x about the x-axis, we may treat these by applying the method to two sets of Wyckoff positions, the one which contains the point (x, y + v, z) and the one containing (x, y, z + w), where v and w are small displacements. In order to carry out this analysis, since the three structure determinations chosen for comparison had each been refined independently in the standard setting of the appropriate space group, the cell axes and coordinates of the low temperature structures were transformed to bring each of the reference molecules into correspondence with each other and with the reference molecule at room temperature, fixing as identical the z-coordinates of the centre of mass in all three cases.
The centre of mass was found to be displaced along x by 25 pm at 120 K and 17 pm at 274 K. Since the displacements of the atoms from their room temperature positions at 274 K proved to be uniformly two-thirds of those at 120 K, the subsequent analysis concentrated on comparing the structures at 120 K and 293 K. No satisfactory fit could be found of the low temperature molecule to the room temperature molecule by simple rotational displacements about the centre of mass. This difficulty was traced to an internal deformation of the molecule, when it became clear from the observed torsion angles that the molecule at 120 K was folded along the N1 – N3 axis. Using this model to fit the internal non-bonded distances indicated a folding angle of 9.1°, which was apportioned in transforming the coordinates in the ratio 0.35 : 0.65, the inverse ratio of the relevant moments of inertia about N1 – N3. After readjusting to the centre of mass, the positions of the non-hydrogen atoms in the two structures could be brought into good agreement (standard error 5 pm) by very small rotational adjustments.
The other phenomenon that should be considered is macroscopic strain, i.e. there may be a ferroelastic aspect to the phase transition. Elongation or compression of the crystal in the appropriate direction should alter the form of the potential function for the normal mode in question. Since the displacements associated with the soft mode lie along a, we would predict the direction of maximum thermal expansion to lie in that direction, and as well as any actual isothermal shear. The lattice parameters were compared  over the range 200 to 320 K. The thermal expansion remains fairly uniform in the range 200 to 290 K with principal components of the expansivity tensor αa = 99(3), αb = 45(3) and αc = 44(3) × 10-6 K-1, with αa being by far the largest as would be expected. However there also appears to be a small isothermal shear, affecting mainly the a-axis, at about 295 K (full data and plots of unit cell length as a function of temperature are provided in Figure S3, see Additional File 5).
The NQR  clearly shows the relative symmetries of the room and low temperature phases in the number of lines observed, and the way the corresponding lines appear to converge supports a model where there are no abrupt changes across the transition. However the lines broaden to the point of disappearance over a temperature range of about 20 K, from a little below the lower transition to near the upper transition temperature. This suggests that the fluctuations present in this temperature range must be such as to produce a very large range of individual chlorine atom environments.
The DSC measurements indicate that the transition occurs as a two stage process with an intermediate crystallisation event. The transition is observed to be reproducible and become sharper as the sample is annealed on cycling. A sharp melting point is observed and the behaviour of the bulk material is shown to mimic that of the small amount of pure material.
The various observations may be interpreted in the following way in terms of cooling the crystal. Above room temperature the molecule lies in a crystallographic mirror plane in space group Cmc21, but there are large vibrational displacements in the x-direction, which couple to a flexing motion along the long molecular axis. At approximately 290 K there is a weak first order isosymmetric transition to an intermediate state in which the molecules occupy two symmetry related positions in close proximity to each other, i.e. the Wyckoff positions of Cmc21 break down to those of Pbn21, but these are occupied in a disordered fashion so as to maintain on average Cmc21 symmetry.
As the temperature is lowered further through this intermediate region, the two molecular sites move steadily apart, and their correlation in terms of the displacement directions of neighbouring molecules becomes progressively longer in range, until about 274 K the system undergoes a second weak first order transition, where the crystal assumes the ordered low temperature structure in Pbn21.
The authors would like to thank the EPSRC for funding the National Crystallographic Service (Southampton, UK). We thank the Shin Nisso Kako Co Ltd for gifts of N3P3Cl6, and Gebze Institute of Technology (GIT) Research Fund for partial support (AU and AK) and the financial support provided (grant no: TBAG-1817) of the Scientific and Technical Research Council of Turkey.
- Contractor SR, Hursthouse MB, Shaw (née Gözen) LS, Shaw RA, Yilmaz H: Structural Investigations of Phosphorus-Nitrogen Compounds. 1. The Structures of Three Monospiro Compounds: N3P3[O(CH2)2O]Cl4, N3P3[O(CH2)3O]Cl4 and N3P3[O(CH2)4O]Cl4. The Relationship of OPO Bond Angles in PO2N2 Tetrahedra with 31P Chemical Shifts. Acta Crystallogr. 1985, B41: 122-131.View ArticleGoogle Scholar
- Hursthouse MB, Porte AL, Shaw (née Gözen) LS, Shaw RA: X-Ray crystallographic and C1-35 NQR spectroscopic studies on chlorocyclophosphazenes. Phosphorus Sulfur. 1986, 28: 213-220. 10.1080/03086648608072811.View ArticleGoogle Scholar
- Connelly A, Dalgleish VH, Harkins P, Keat R, Porte AL, Raitt I, Shaw RA: Chlorine-35 NQR studies of some chlorocyclotriphosphazatrienes. J Magnetic Resonance. 1978, 30: 439-450.Google Scholar
- Contractor SR, Hursthouse MB, Parkes HG, Shaw (née Gözen) LS, Shaw RA, Yilmaz H: Preparation, spectroscopic and structural investigations of spiro derivatives of hexachlorocyclotriphosphazatriene. Phosphorus and Sulfur. 1983, 18: 149-152. 10.1080/03086648308075989.View ArticleGoogle Scholar
- Alkubaisi AH, Parkes HG, Shaw RA: Phosphorus-Nitrogen Compounds. Part 58. The Reactions of Hexachlorocyclotriphosphazatriene with Ethane-, 1,3-Propane- and 1,4-Butane-diols. Spiro, Ansa, Bridged and Dangling Derivatives and their 31P and 1H Nuclear Magnetic Resonance Spectra. Heterocycles. 1989, 28: 347-358.View ArticleGoogle Scholar
- Hooft RWW: COLLECT: Data Collection Software. Edited by: Nonius BV. 1998, Delft, The NetherlandsGoogle Scholar
- Otwinowski Z, Minor W: Processing of X-ray diffraction data collected in oscillation mode. Methods in Enzymology. 1997, 276: 307-View ArticleGoogle Scholar
- Sheldrick GM: SHELX97: Programs for structure solution and refinement. 1997, University of Göttingen, GermanyGoogle Scholar
- Sheldrick GM: SADABS. Version 2.10. 2003, Bruker AXS Inc., Madison, Wisconsin, USAGoogle Scholar
- Cosier J, Glazer AM: A nitrogen-gas-stream cryostat for general X-ray diffraction studies. J Appl Crystallogr. 1986, 19: 105-107. 10.1107/S0021889886089835.View ArticleGoogle Scholar
- Data collection and refinement parameters for the 120 K structure. [http://ecrystals.chem.soton.ac.uk/299]
- Data collection and refinement parameters for the 274 K structure. [http://ecrystals.chem.soton.ac.uk/300]
- Data collection and refinement parameters for the 293 K structure. [http://ecrystals.chem.soton.ac.uk/301]
- Toledano J-C, Toledano P: The Landau Theory of Phase Transitions: Application to Structural, Incommensurate, Magnetic, and Liquid Crystal Systems. 1987, World Scientific, SingaporeView ArticleGoogle Scholar
- Wondratschek H, Müller U: International Tables for Crystallography. Volume A1: Symmetry relations between space groups. 2004, Springer, BerlinGoogle Scholar
- Cochran W: The Dynamics of Atoms in Crystals. 1979, Arnold, LondonGoogle Scholar
- Aroyo MI, Perez-Mato JM: Symmetry-Mode Analysis of Displacive Phase Transitions using International Tables for Crystallography. Acta Crystallogr. 1998, A54: 19-30.View ArticleGoogle Scholar
- Burnett MN, Johnson CK: ORTEP-III: Oak Ridge Thermal Ellipsoid Plot Program for Crystal Structure Illustrations. 1996, Oak Ridge National Laboratory Report ORNL-6895View ArticleGoogle Scholar
- Rutherford JS: On comparing lattice parameters among isostructural molecular crystals. Acta Chim Hung – Models in Chemistry. 1997, 134: 395-405.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.