research letters
ζ-Glycine: insight into the mechanism of a polymorphic phase transition
aISIS Facility, Rutherford–Appleton Laboratory, Harwell Science and Innovation Campus, Didcot, Oxfordshire OX11 0QX, UK, bSchool of Chemistry and Centre for Science at Extreme Conditions, The University of Edinburgh, King's Buildings, W. Mains Road, Edinburgh EH9 3FJ, UK, cScuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, Trieste 34136, Italy, dDepartment of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA 92093, USA, and eOak Ridge National Laboratory, PO Box 2008, Oak Ridge, TN 37831, USA
*Correspondence e-mail: degironc@sissa.it, kucukben@gmail.com, s.parsons@ed.ac.uk
Glycine is the simplest and most polymorphic amino acid, with five phases having been structurally characterized at atmospheric or high pressure. A sixth form, the elusive ζ phase, was discovered over a decade ago as a short-lived intermediate which formed as the high-pressure ɛ phase transformed to the γ form on decompression. However, its structure has remained unsolved. We now report the structure of the ζ phase, which was trapped at 100 K enabling neutron powder diffraction data to be obtained. The structure was solved using the results of a prediction procedure based on fully ab initio energy calculations combined with a for searching phase space. We show that the fate of ζ-glycine depends on its thermal history: although at room temperature it transforms back to the γ phase, warming the sample from 100 K to room temperature yielded β-glycine, the least stable of the known ambient-pressure polymorphs.
1. Introduction
Glycine, the simplest amino acid, exhibits greater phase diversity than any of the other naturally occurring amino acids. All phases consist of the zwitterionic form +H3N—CH2—CO2− and three polymorphs are known under ambient conditions. The α and β forms are both monoclinic (space groups P21/n and P21, respectively), with crystal structures composed of hydrogen-bonded layers. The γ form is trigonal (P31) and displays a three-dimensional hydrogen-bonded network featuring chains disposed about 31 screw axes. Of the three polymorphs, γ-glycine is the most thermodynamically stable at room temperature (Perlovich et al., 2001).
The effect of pressure on the polymorphs of glycine has been extensively investigated. Although α-glycine persists to 23 GPa (Murli et al., 2003), β-glycine undergoes a at only 0.8 GPa to the δ phase (P21/n; Dawson et al., 2005; Goryainov et al., 2005). This transition is displacive and fully reversible as a result of the topological similarity of the β and δ phases, consisting of a concerted inversion of the O—C—C—N torsion angles in half the molecules which otherwise retain their positions and orientations. Application of pressure to γ-glycine leads to ɛ-glycine (Pn; Boldyreva et al., 2005; Dawson et al., 2005; Moggach et al., 2015). The structure changes substantially, the threefold helices being replaced by a layered structure. The transition is sluggish and reconstructive, but is complete at between 4 and 5 GPa. Note that some authors refer to the δ and ɛ phases as β′ and δ, respectively.
The γ- to ɛ-glycine transition shows considerable hysteresis and the ɛ phase persists on decompression to 0.28 GPa (Moggach et al., 2015). Complete release of applied load yields yet another phase, ζ-glycine, which exists for as little as 30 min at room temperature before transforming back to the γ phase (Bordallo et al., 2008; Goryainov et al., 2006; Moggach et al., 2015). The ζ phase, which has been characterized by X-ray powder diffraction and Raman spectroscopy, is only observed upon decompression of the ɛ form; application of pressure to the γ phase yields the ɛ phase directly.
Although the ζ phase was first identified over a decade ago, its structure is unknown. On the basis of Raman data it has been proposed to have a layered structure similar to all the other phases except for γ-glycine (Bordallo et al., 2008). In this paper we describe the solution of the of ζ-glycine using high-quality neutron powder diffraction measurements in combination with prediction based on fully first-principles total energy calculations and an improved genetic search algorithm for searching the phase space.
2. Experimental
2.1. Neutron powder diffraction
A sample of γ-glycine-d5 (CDN Isotopes, Canada) was obtained by non-photochemical laser-induced nucleation of a supersaturated solution in D2O (Zaccaro et al., 2001). Colourless needle-like crystals were allowed to grow for 2 h at room temperature and then lightly crushed to form a polycrystalline sample. Pressure was applied to the sample using a Paris–Edinburgh press modified for use at low temperature, utilizing zirconia-toughened alumina anvils. The sample was enclosed in a null-scattering TiZr encapsulated gasket. A mixture of deuterated methanol and ethanol (4:1 v/v) was used as a pressure transmitting medium with a lead pellet as pressure marker. Diffraction data were obtained on the PEARL instrument at ISIS (Bull et al., 2016) and the resulting diffaction corrected for anvil attenuation. The sample was compressed to 5.56 GPa to yield a phase-pure sample of ɛ-glycine. The temperature was then rapidly lowered to 200 K and the load reduced to 0 tonnes. Then the temperature was further reduced to 100 K. Neutron powder diffraction data were collected over the following 5 h. Following data collection, the sample was warmed to room temperature while monitoring the diffraction pattern. Further experimental details are available in the supporting information.
2.2. prediction
The evolutionary algorithm as implemented in the USPEX package (Lyakhov et al., 2013) was used to search for the low-energy structures of glycine with Z = 2, 3 or 4 molecules in the In the first generation 30 structures were created randomly. The energetically least favourable 20% of the population was discarded. A fingerprint analysis of the remaining structures was performed (Oganov & Valle, 2009). The structures whose fingerprint was within an adimensional cosine distance of 0.01 from any lower energy structure were also discarded. The remaining structures were then eligible as parents and allowed to `procreate'. Thirty new structures of the next generation were created from the parents through the following operations: heredity (crossover of two structures) (40%), soft mutation (translation and rotation based on an estimate of soft vibrational modes) (20%), rotation of the molecule (20%) and random structure generation (20%). In addition, the three best parents were cloned directly to the next generation. In all simulations, the maximum number of generations was 20.
2.3. Ab initio calculations
For every structure generated by USPEX, the geometry and cell relaxation were calculated using a van der Waals density functional (vdW-DF) as implemented in the QUANTUM ESPRESSO package (Giannozzi et al., 2009). A cutoff of 80 Ryd and a charge density cutoff of 560 Ryd were used. The sampling resolution was gradually increased over three steps during relaxation, from 2π × 0.12 Å−1 to 2π × 0.10 Å−1 and finally 2π × 0.08 Å−1. Energies and geometries from the last step with the densest k-point grid were used throughout the study. Projector augmented wave (PAW) pseudopotentials were taken from the PSLibrary v0.1 (Küçükbenli et al., 2014). By using this setup all structures were fully relaxed within a convergence of less than 0.1 mRyd for the absolute total energy per molecule, 0.5 mRyd a.u.−1 for the forces on the atoms and less than 0.005 GPa for the pressure obtained from the stress tensor.
3. Results
3.1. prediction
Conventional methods for indexing and solving the structure of ζ-glycine failed to yield models which reproduced the experimental data. However, independently of the diffraction experiments, theoretical work had been carried out to assess the predictive power of a newly developed prediction (CSP) workflow using glycine as the test case. Although a number of CSP surveys of glycine have been described (Zhu et al., 2012; Lund et al., 2015; Chisholm et al., 2005), none has so far identified candidates which fit the experimental data for ζ-glycine, though a combined CSP–powder diffraction study has very recently yielded the of glycine dihydrate (Xu et al., 2017).
A critical component of CSP is the method of evaluating the lattice energies of candidate phases. This determines both energy ranking and the direction of the phase-space search. The highest quality energies are obtained using first-principles calculations based on periodic density functional theory (DFT) (Hohenberg & Kohn, 1964; Reilly et al., 2016), but despite some spectacular successes (Woodley & Catlow, 2008) molecular crystals are challenging because of the difficulties in accounting for dispersion interactions (Reilly et al., 2016).
In this study energies were calculated (QUANTUM ESPRESSO; Giannozzi et al., 2009) with a recently developed functional (vdW-DF; Berland et al., 2015) which takes dispersion into account via a non-local functional of the overall charge density. The method yields reliable energy differences between phases of molecular materials (Sabatini et al., 2012), including the ambient-pressure forms of glycine, without suffering from the transferability problems of empirical methods. The combination of this method for energy evaluation with the for searching phase space (USPEX; Lyakhov et al., 2013) assuming two, three or four molecules per yielded a very rich series of hypothetical structures for glycine (Fig. 1). Fuller details are available in the supporting information.
In addition to all the experimentally observed phases of glycine, the CSP analysis predicted several new phases within 2 kJ mol−1 of the most stable (γ) phase. In particular, the survey had identified a potential new polymorph in P1 with one molecule per (Z = 1; a = 5.0168, b = 4.7491 and c = 4.0593 Å; α = 95.8383, β = 105.6522 and γ = 64.8726°) that was very close in energy to the known β polymorph. This was proposed as ζ-glycine based on the agreement between previous mixed-phase X-ray diffraction experiments (Boldyreva et al., 2005) and theoretical assessment (see supporting information). The structural parameters of the new phase, which had been predicted completely ab initio, were found to form a suitable starting model for of the structure against the neutron diffraction data.
3.2. Neutron powder diffraction
Neutron powder diffraction data suitable for structure solution were collected for ζ-glycine by trapping the phase at 100 K, extending its lifetime for long enough to collect a high-quality diffraction pattern.
The powder diffraction data were indexed using the unit-cell dimensions obtained in the structure prediction multiplied by a factor of 0.98. This accounts for the tendency of the method used for CSP to overestimate slightly the unit-cell dimensions of other glycine polymorphs (Table S1). Pawley fitting yielded the following cell dimensions for ζ-glycine: a = 4.9307 (3), b = 4.54798 (4) and c = 3.9191 (3) Å, α = 95.550 (5), β = 105.250 (5) and γ = 64.938 (6)°, and V = 76.797 (11) Å3. These values can be transformed to a body-centred triclinic setting containing two molecules per which clarifies the structural relationship of ζ-glycine with the other phases. These transformed cell dimensions are compared with those of the other Z = 2 phases in Table 1. The predicted of ζ-glycine was refined using Rietveld methods (TOPAS-Academic; Coelho, 2015). The final fit is shown in Fig. 2(a).
|
When the sample of ζ-glycine prepared in this study was warmed from 100 K, the peaks due to the residual ɛ phase were first seen to disappear at 250 K. At 290 K the peaks due to ζ-glycine began to diminish, to yield not the expected γ phase but rather the β phase (Fig. 2b). Persistence of the ζ polymorph after the disappearance of the ɛ polymorph shows that the transition to β-glycine is from ζ-glycine rather than ɛ-glycine; the initial disappearance of ɛ-glycine is the continuation of the original ɛ to ζ transition halted by the temperature decrease.
4. Discussion
4.1. Structural relationships and phase stability
The data in Table 1 show that the unit-cell dimensions of the β, ɛ and ζ polymorphs are very closely related to each other. The structures, along with that of γ-glycine, are shown in Figs. S13–S16. The β, ɛ and ζ phases are layered structures in which the molecules in the reside near positions equivalent to and (Table S4). A crystallographic fractional coordinate of can be transformed to either by adding or by inverting it and adding 1. The β, ɛ and ζ polymorphs can simply be viewed as the result of various combinations of these operations along the x, y and z directions of the In the ζ phase the relationship is x + , y + , z + (I-centring); in the ɛ phase it is x + , −y + 1, z + (an n-glide plane), while in the β phase it is −x + 1, y + , −z + 1 (a 21 screw axis).
In the case of the β and ζ phases the symmetry operations preserve conformation, and this is why the cell dimensions of these two phases are so similar. In ɛ-glycine the two zwitterionic molecules are related by an n-glide and have opposite conformations, the N1—C2—C1—O2 torsion angles being ±16.0 (5)°. In ζ-glycine the torsion angle increases to 35.1 (4)° and all the molecules have the same conformation. The ɛ- to ζ-glycine transition can thus be viewed as a switch of conformation of half the molecules in the structure, and in this regard it is similar to the β- to δ-glycine transition at 0.8 GPa. Although γ-glycine shares the polarity of the ɛ and ζ phases along c, and since its is chiral all three molecules in the share the same conformation [τ(N1—C2—C1—O2) = 15.6 (1)°; Kvick et al., 1980], its structure consists of a three-dimensional hydrogen-bonded network rather than layers.
Symmetry-adapted perturbation theory (PSI4 software; Turney et al., 2012) and PIXEL calculations (Gavezzotti, 2005) (see supporting information) both indicate that the strongest intermolecular interaction in ζ-glycine, with an energy in excess of −100 kJ mol−1, is a `head-to-tail' hydrogen bond formed between ammonium and carboxylate groups. Repetition of this interaction builds a chain which runs along the c direction of the This same chain is present in all six known polymorphs, and the phases differ in the way these chains interact.
In ζ-glycine the chains are arranged in the ac plane, related by lattice repeats along a, forming a layer (Fig. 3a). This layer motif confirms Boldyreva's conclusion regarding the layered nature of the ζ phase made on the basis of vibrational spectroscopy (Bordallo et al., 2008). The layers are very similar, both in terms of geometry and in the molecule–molecule energies, to those seen for the ɛ phase. The intermolecular contact energy across the N1—H5⋯O1 hydrogen bonds (−26 kJ mol−1) is only about a quarter of that of the head-to-tail linkages, a result of the repulsive influence of neighbouring carboxylate groups.
In the ζ phase all layers are connected through equivalent N1—H4⋯O contacts (Fig. 3b). The layer separation is b/2 = 3.14 Å. Although the intermolecular contact energy for the pairs of molecules linked by N1—H4⋯O1 contacts (−59.4 kJ mol−1) is equivalent to a strong hydrogen bond, the geometry of the contact deviates markedly from linearity [∠N1—H4⋯O1 = 121.1 (10)°] with a relatively long H⋯O distance [2.274 (5) Å], and the interaction is more plausibly classified as a simple electrostatic interaction involving ammonium and carboxylate groups.
The relative energies (in kJ mol−1) of the six phases of glycine, calculated using DFT and including the approximation for van der Waals interactions described above, are: γ 0, α 0.084, β 1.055, ζ 1.070, ɛ 1.832 and δ 2.019. The energies refer to geometry-optimized structures at 0 K and 0 Pa. The ordering of the α, β and γ phases agrees with previous work (Marom et al., 2013; Perlovich et al., 2001; Sabatini et al., 2012), while the total range, which spans 2 kJ mol−1, agrees with Hunter's recent results on typical polymorph energy differences (Hunter & Prohens, 2017). The ζ phase is marginally less stable than the previously known ambient-pressure phases, but more stable than either of the high-pressure phases.
4.2. The fate of ζ-glycine on warming
Previous work (Boldyreva et al., 2005; Bordallo et al., 2008; Goryainov et al., 2006; Moggach et al., 2015) using X-ray and neutron powder diffraction and vibrational spectroscopy has shown that ζ-glycine transforms spontaneously to the γ phase at room temperature. By contrast, in the present investigation, warming the sample from 100 K to room temperature yielded β-glycine rather than the expected γ phase. The thermal history of the samples used in this and previous work were different: in previous studies all manipulations were carried out at room temperature, whereas here the sample had been cooled. Nevertheless the difference in behaviour is surprising, and this is the first time, to our knowledge, that a transition to the metastable β phase of glycine has been observed. Similar sensitivity of phase formation to thermal history has been observed, for example in paracetamol (Qi et al., 2008; Rossi et al., 2003), while in glycine itself the temperature of the thermal γ-to-α can be increased by ca 10 K by annealing the sample (Perlovich et al., 2001).
Over the course of the ɛ–ζ–γ transition (see the movie in the supporting information), the layers of ɛ-glycine first slide over one another with small molecular rotations to give ζ-glycine. Larger reorientations that disrupt the layer structure are needed as the system transforms to the γ phase.
As discussed above, the β and ζ phases are closely related. Both consist of chains composed of head-to-tail N1—H3⋯O2 hydrogen bonds which are connected into a layer via N1—H5⋯O1 hydrogen bonds. The difference between the phases is that the c axis is a polar direction in the ζ phase but not in the β phase. Therefore, in the ζ–β transition (see the movie in the supporting information) the layer structure is retained but a rotation of layers with respect to one another is needed. The lack of a strong directional preference in the electrostatic contacts between the layers (see above) may explain why such phase transitions involving rearrangements of layer stacking can occur.
The structure of ζ-glycine is significant because the phase is akin to a supramolecular reactive intermediate, providing insight into the mechanism of solid-state phase transitions. The results presented here, along with those of previous studies (Boldyreva et al., 2005; Bordallo et al., 2008; Goryainov et al., 2006; Moggach et al., 2015), indicate that two transformation pathways are available to ζ-glycine, one leading to the β phase and the other to the γ phase. One possibility is that the activation barrier to the β phase is lower, and it is possible that seeds of this phase begin to form at low temperature. In this interpretation, at ambient temperature there is enough thermal energy to form the thermodynamically more stable γ phase. Alternatively, the topological similarity of the ɛ and γ phases, in which all the head-to-tail chains have the same polarity, may enable a residual amount of the high-energy ɛ phase to seed the formation of γ-glycine. Further theoretical modelling and experimental investigations would provide valuable insight into the physical source of the difference between the low- and ambient-temperature behaviour.
5. Related literature
For additional literature relating to the supporting information, see Brandenburg (2004), Destro et al. (2000), Fortes et al. (2007, 2012), Gavezzotti (2011), Hohenstein & Sherrill (2010a,b, 2012), Jeziorski et al. (1994), Frisch et al. (2009), Macrae et al. (2008), Sheldrick (2001), Spek (2003) and Stone (2013).
Supporting information
https://doi.org/10.1107/S205225251701096X/lq5008sup1.cif
and TOPAS data. DOI:Further details relating to the https://doi.org/10.1107/S205225251701096X/lq5008sup2.pdf
prediction and crystallographic structure refinements. Figures showing the structures of all the phases discussed are given, along with an energy analysis of intermolecular contacts in the epsilon and zeta phases. DOI:Movie (animated gif, ezb.gif). DOI: https://doi.org/10.1107/S205225251701096X/lq5008sup3.gif
C2D5NO2 | γ = 104.136 (5)° |
Mr = 80.10 | V = 153.54 (2) Å3 |
Triclinic, I1 | Z = 2 |
Hall symbol: I1 | Cell parameters included in refinement |
a = 5.1000 (4) Å | Dx = 1.732 Mg m−3 |
b = 6.2850 (3) Å | Neutron radiation |
c = 5.4295 (3) Å | T = 100 K |
α = 85.815 (5)° | Particle morphology: Block |
β = 114.456 (5)° | Colourless |
PEARL, ISIS diffractometer | Scan method: Time-of-flight |
Radiation source: Pulsed neutron source | Absorption correction: empirical (using intensity measurements) Calculated using custom ISIS software. |
Specimen mounting: TiZr Capsule | |
Data collection mode: transmission |
Least-squares matrix: selected elements only | 0 restraints |
Rp = 2.866 | 22 constraints |
Rwp = 3.088 | Weighting scheme based on measured s.u.'s |
Excluded region(s): excluded short TOF region above | (Δ/σ)max = <0.001 |
Profile function: Voigt function convoluted with a back-to-back double exponential. | Background function: Chebychev polynomial |
65 parameters | Preferred orientation correction: Modelling of | was found to be unnecessary.
Geometry. Distances and angles were all constrained in the refinement so the sus are 0 [this won't be the case if the values are relcalculated from the using the coordinate data block above - getting the right values depends on having the fill covariance matrix available]. |
Refinement. The data in this cif come from a Rietveld refinement against neutron powder data. Details of refinement are given in the paper, and the Topas instructuion file is under _iucr_refine_instructions_details. Checkcif output: 091_ALERT_1_A No Wavelength found in CIF - 0.71073 Ang Assumed Please Check The data were collected by the time of flight method 155_ALERT_4_C The Triclinic Unitcell is NOT Reduced ·········. Please Do ! 128_ALERT_4_G Alternate Setting for Input Space Group I1 P1 Note The setting is chosen to show the relationship to the epsilon phase. The matrix relating the prediction and refinement cell settings is [-1 1 0] [ 0 -1 1] [ 1 0 1] 161_ALERT_4_C Missing or Zero su (esd) on x-coordinate for ··· N1 162_ALERT_4_C Missing or Zero su (esd) on y-coordinate for ··· N1 163_ALERT_4_C Missing or Zero su (esd) on z-coordinate for ··· N1 This atom was used to fix the origin. 340_ALERT_3_C Low Bond Precision on C-C Bonds ··············· 0.0080 Ang. 201_ALERT_2_G Isotropic non-H Atoms in Main Residue(s) ······. 5 Why ? 210_ALERT_3_G No Anisotropic ADP's Found in CIF ············.. Please Check The structure was refined against medium resolution neutron powder data. 741_ALERT_1_C Bond Calc 1.25 (2), Rep 1.25000 ······ Missing su O1 -C1 1.555 1.555 # 1 741_ALERT_1_C Bond Calc 1.250 (11), Rep 1.25000 ······ Missing su O2 -C1 1.555 1.555 # 2 ···plus numerous similar alerts The s.u.s on the diatnces and angles are all exactly 0 as they were constrained during refinement. 017_ALERT_1_G Check Consistency of Scattering Type D for H1 017_ALERT_1_G Check Consistency of Scattering Type D for H2 017_ALERT_1_G Check Consistency of Scattering Type D for H3 017_ALERT_1_G Check Consistency of Scattering Type D for H4 017_ALERT_1_G Check Consistency of Scattering Type D for H5 The sample was deuterated. 154_ALERT_1_G The su's on the Cell Angles are Equal ·········. 0.00500 Degree 720_ALERT_4_G Number of Unusual/Non-Standard Labels ·········. 5 Note 802_ALERT_4_G CIF Input Record(s) with more than 80 Characters ! Info No action taken |
x | y | z | Biso*/Beq | ||
N1 | 0.71144 | 0.68132 | 0.77926 | 0.58 (3) | |
C1 | 0.540 (2) | 0.777 (2) | 1.1166 (9) | 0.58 (3) | |
C2 | 0.518 (2) | 0.7995 (16) | 0.8265 (8) | 0.58 (3) | |
O1 | 0.308 (3) | 0.766 (4) | 1.1491 (13) | 0.58 (3) | |
O2 | 0.788 (2) | 0.7714 (17) | 1.2992 (6) | 0.58 (3) | |
H1 | 0.2893 (15) | 0.732 (3) | 0.6865 (9) | 0.69 (4) | |
H2 | 0.584 (4) | 0.9728 (15) | 0.7886 (14) | 0.69 (4) | |
H3 | 0.7178 (14) | 0.7244 (11) | 0.5909 (7) | 0.87 (5) | |
H4 | 0.625 (2) | 0.5116 (2) | 0.7800 (17) | 0.87 (5) | |
H5 | 0.9255 (5) | 0.7260 (15) | 0.9341 (9) | 0.87 (5) |
O1—C1 | 1.2500 | N1—H4 | 1.0500 |
O2—C1 | 1.2500 | N1—H5 | 1.0500 |
N1—C2 | 1.4800 | C2—H1 | 1.0900 |
C1—C2 | 1.5300 | C2—H2 | 1.0900 |
N1—H3 | 1.0500 | ||
O1—C1—O2 | 126.00 | H3—N1—H5 | 110.00 |
O1—C1—C2 | 117.00 | N1—C2—C1 | 112.00 |
O2—C1—C2 | 117.00 | N1—C2—H1 | 108.00 |
C2—N1—H5 | 109.00 | N1—C2—H2 | 108.00 |
C2—N1—H3 | 109.00 | C1—C2—H1 | 109.00 |
C2—N1—H4 | 109.00 | C1—C2—H2 | 109.00 |
H4—N1—H5 | 110.00 | H1—C2—H2 | 110.00 |
H3—N1—H4 | 110.00 | ||
O1—C1—C2—N1 | −144.9 (4) | O2—C1—C2—N1 | 35.1 (4) |
C2D5NO2 | V = 150.12 (2) Å3 |
Mr = 80.10 | Z = 2 |
Monoclinic, Pn | Cell parameters included in refinement |
Hall symbol: P -2yac | Dx = 1.772 Mg m−3 |
a = 5.0230 (4) Å | Neutron radiation |
b = 5.9846 (4) Å | T = 100 K |
c = 5.4946 (5) Å | Particle morphology: Block |
β = 114.654 (8)° | Colourless |
PEARL, ISIS diffractometer | Scan method: Time-of-flight |
Radiation source: Pulsed neutron source | Absorption correction: empirical (using intensity measurements) Calculated using custom ISIS software. |
Specimen mounting: TiZr Capsule | |
Data collection mode: transmission |
Least-squares matrix: selected elements only | 0 restraints |
Rp = 2.866 | 21 constraints |
Rwp = 3.088 | Weighting scheme based on measured s.u.'s |
Excluded region(s): excluded short TOF region above | (Δ/σ)max = <0.001 |
Profile function: Voigt function convoluted with a back-to-back double exponential. | Background function: Chebychev polynomial |
65 parameters | Preferred orientation correction: Modelling of | was found to be unnecessary.
Geometry. Distances and angles were all constrained in the refinement so the sus are 0 [this won't be the case if the values are relcalculated from the using the coordinate data block above - getting the right values depends on having the fill covariance matrix available]. |
Refinement. The data in this cif come from a Rietveld refinement against neutron powder data. Details of refinement are given in the paper. Checkcif output: 091_ALERT_1_A No Wavelength found in CIF - 0.71073 Ang Assumed Please Check Time of flight data. 340_ALERT_3_B Low Bond Precision on C-C Bonds ··············· 0.0110 Ang. 201_ALERT_2_G Isotropic non-H Atoms in Main Residue(s) ······. 5 Why ? 210_ALERT_3_G No Anisotropic ADP's Found in CIF ············.. Please Check The structure was refined against medium resolution neutron powder data. 161_ALERT_4_C Missing or Zero su (esd) on x-coordinate for ··· N1 163_ALERT_4_C Missing or Zero su (esd) on z-coordinate for ··· N1 These coordinates were used to fix the origin. 741_ALERT_1_C Bond Calc 1.25 (2), Rep 1.25000 ······ Missing su O1 -C1 1.555 1.555 # 1 741_ALERT_1_C Bond Calc 1.250 (14), Rep 1.25000 ······ Missing su O2 -C1 1.555 1.555 # 2 ···Plus similar alerts The s.u.s on the diatnces and angles are all exactly 0 as they were constrained during refinement. 017_ALERT_1_G Check Consistency of Scattering Type D for H1 017_ALERT_1_G Check Consistency of Scattering Type D for H2 017_ALERT_1_G Check Consistency of Scattering Type D for H3 017_ALERT_1_G Check Consistency of Scattering Type D for H4 017_ALERT_1_G Check Consistency of Scattering Type D for H5 The sample was deuterated. 720_ALERT_4_G Number of Unusual/Non-Standard Labels ·········. 5 Note 802_ALERT_4_G CIF Input Record(s) with more than 80 Characters ! Info No action taken |
x | y | z | Biso*/Beq | ||
O1 | 0.292 (3) | 0.795 (5) | 1.147 (3) | 0.42 (4) | |
O2 | 0.747 (2) | 0.667 (2) | 1.2731 (10) | 0.42 (4) | |
N1 | 0.734 | 0.6927 (12) | 0.78 | 0.42 (4) | |
C1 | 0.514 (2) | 0.751 (3) | 1.1073 (15) | 0.42 (4) | |
C2 | 0.500 (2) | 0.806 (2) | 0.8301 (13) | 0.42 (4) | |
H1 | 0.2887 (10) | 0.752 (4) | 0.6782 (16) | 0.50 (5) | |
H2 | 0.525 (4) | 0.986 (2) | 0.816 (2) | 0.50 (5) | |
H3 | 0.7163 (17) | 0.738 (2) | 0.5891 (10) | 0.63 (6) | |
H4 | 0.709 (3) | 0.5191 (12) | 0.790 (2) | 0.63 (6) | |
H5 | 0.9385 (9) | 0.743 (2) | 0.9265 (15) | 0.63 (6) |
O1—C1 | 1.2500 | N1—H4 | 1.0500 |
O2—C1 | 1.2500 | N1—H5 | 1.0500 |
N1—C2 | 1.4800 | C2—H1 | 1.0900 |
C1—C2 | 1.5300 | C2—H2 | 1.0900 |
N1—H3 | 1.0500 | ||
O1—C1—O2 | 126.00 | H3—N1—H5 | 110.00 |
O1—C1—C2 | 117.00 | N1—C2—C1 | 112.00 |
O2—C1—C2 | 117.00 | N1—C2—H1 | 108.00 |
C2—N1—H5 | 109.00 | N1—C2—H2 | 108.00 |
C2—N1—H3 | 109.00 | C1—C2—H1 | 109.00 |
C2—N1—H4 | 109.00 | C1—C2—H2 | 109.00 |
H4—N1—H5 | 110.00 | H1—C2—H2 | 110.00 |
H3—N1—H4 | 110.00 | ||
O1—C1—C2—N1 | −164.0 (5) | O2—C1—C2—N1 | 16.0 (5) |
C2D5NO2 | V = 157.71 (1) Å3 |
Mr = 80.10 | Z = 2 |
Monoclinic, P21 | Cell parameters included in refinement |
Hall symbol: P 2yb | Dx = 1.686 Mg m−3 |
a = 5.0907 (2) Å | Neutron radiation |
b = 6.25954 (16) Å | T = 290 K |
c = 5.38710 (19) Å | Particle morphology: Block |
β = 113.261 (4)° | Colourless |
PEARL, ISIS diffractometer | Scan method: Time-of-flight |
Radiation source: Pulsed neutron source | Absorption correction: empirical (using intensity measurements) Calculated using custom ISIS software. |
Specimen mounting: TiZr Capsule | |
Data collection mode: transmission |
Least-squares matrix: selected elements only | 0 restraints |
Rp = 3.220 | 21 constraints |
Rwp = 3.550 | Weighting scheme based on measured s.u.'s |
Excluded region(s): excluded short TOF region above | (Δ/σ)max = <0.001 |
Profile function: Voigt function convoluted with a back-to-back double exponential. | Background function: Chebychev polynomial |
62 parameters | Preferred orientation correction: Modelling of | was found to be unnecessary.
Geometry. Distances and angles were all constrained in the refinement so the sus are 0 [this won't be the case if the values are relcalculated from the using the coordinate data block above - getting the right values depends on having the fill covariance matrix available]. |
Refinement. The data in this cif come from a Rietveld refinement against neutron powder data. Details of refinement are given in the paper, and the Topas instructuion file is under _iucr_refine_instructions_details. Checkcif output: 091_ALERT_1_A No Wavelength found in CIF - 0.71073 Ang Assumed Please Check The data were collected by the time of flight method 162_ALERT_4_C Missing or Zero su (esd) on y-coordinate for ··· N1 This atom was used to fix the origin. 340_ALERT_3_C Low Bond Precision on C-C Bonds ··············· 0.0070 Ang. 201_ALERT_2_G Isotropic non-H Atoms in Main Residue(s) ······. 5 Why ? 210_ALERT_3_G No Anisotropic ADP's Found in CIF ············.. Please Check The structure was refined against medium resolution neutron powder data. 741_ALERT_1_C Bond Calc 1.252 (14), Rep 1.25000 ······ Missing su O1 -C1 1.555 1.555 # 1 741_ALERT_1_C Bond Calc 1.250 (9), Rep 1.25000 ······ Missing su O2 -C1 1.555 1.555 # 2 741_ALERT_1_C Bond Calc 1.480 (8), Rep 1.48000 ······ Missing su N1 -C2 1.555 1.555 # 3 ···plus other alerts like this The s.u.s on the diatnces and angles are all exactly 0 as they were constrained during refinement. 017_ALERT_1_G Check Consistency of Scattering Type D for H1 017_ALERT_1_G Check Consistency of Scattering Type D for H2 017_ALERT_1_G Check Consistency of Scattering Type D for H3 017_ALERT_1_G Check Consistency of Scattering Type D for H4 017_ALERT_1_G Check Consistency of Scattering Type D for H5 The sample was deuterated 720_ALERT_4_G Number of Unusual/Non-Standard Labels ·········. 5 802_ALERT_4_G CIF Input Record(s) with more than 80 Characters ! Info No action taken. |
x | y | z | Biso*/Beq | ||
O1 | 0.410 (2) | 0.378 (3) | 0.5987 (13) | 1.33 (4) | |
O2 | 0.8795 (16) | 0.3217 (13) | 0.7432 (8) | 1.33 (4) | |
N1 | 0.8522 (7) | 0.25 | 0.2370 (7) | 1.33 (4) | |
C1 | 0.6367 (16) | 0.3537 (16) | 0.5629 (10) | 1.33 (4) | |
C2 | 0.6146 (14) | 0.3646 (14) | 0.2716 (9) | 1.33 (4) | |
H1 | 0.4136 (8) | 0.292 (2) | 0.1385 (10) | 1.60 (5) | |
H2 | 0.621 (3) | 0.5315 (14) | 0.2161 (12) | 1.60 (5) | |
H3 | 0.8504 (12) | 0.2867 (10) | 0.0461 (9) | 2.00 (6) | |
H4 | 0.824 (2) | 0.08499 (18) | 0.2520 (17) | 2.00 (6) | |
H5 | 1.0465 (9) | 0.2989 (13) | 0.3890 (11) | 2.00 (6) |
O1—C1 | 1.2500 | N1—H4 | 1.0500 |
O2—C1 | 1.2500 | N1—H5 | 1.0500 |
N1—C2 | 1.4800 | C2—H1 | 1.0900 |
C1—C2 | 1.5300 | C2—H2 | 1.0900 |
N1—H3 | 1.0500 | ||
O1—C1—O2 | 126.00 | H3—N1—H5 | 110.00 |
O1—C1—C2 | 117.00 | N1—C2—C1 | 112.00 |
O2—C1—C2 | 117.00 | N1—C2—H1 | 108.00 |
C2—N1—H5 | 109.00 | N1—C2—H2 | 108.00 |
C2—N1—H3 | 109.00 | C1—C2—H1 | 109.00 |
C2—N1—H4 | 109.00 | C1—C2—H2 | 109.00 |
H4—N1—H5 | 110.00 | H1—C2—H2 | 110.00 |
H3—N1—H4 | 110.00 | ||
O1—C1—C2—N1 | −156.3 (3) | O2—C1—C2—N1 | 23.7 (3) |
C2D5NO2 | γ = 103.55 (3)° |
Mr = 80.10 | V = 155.85 (9) Å3 |
Triclinic, I1 | Z = 2 |
Hall symbol: I1 | Cell parameters included in refinement |
a = 5.1029 (16) Å | Dx = 1.707 Mg m−3 |
b = 6.3450 (12) Å | Neutron radiation |
c = 5.4331 (18) Å | T = 290 K |
α = 85.91 (3)° | Particle morphology: Block |
β = 114.26 (3)° | Colourless |
PEARL, ISIS diffractometer | Scan method: Time-of-flight |
Radiation source: Pulsed neutron source | Absorption correction: empirical (using intensity measurements) Calculated using custom ISIS software. |
Specimen mounting: TiZr Capsule | |
Data collection mode: transmission |
Least-squares matrix: selected elements only | 0 restraints |
Rp = 3.220 | 22 constraints |
Rwp = 3.550 | Weighting scheme based on measured s.u.'s |
Excluded region(s): excluded short TOF region above | (Δ/σ)max = <0.001 |
Profile function: Voigt function convoluted with a back-to-back double exponential. | Background function: Chebychev polynomial |
62 parameters | Preferred orientation correction: Modelling of | was found to be unnecessary.
Geometry. Distances and angles were all constrained in the refinement so the sus are 0 [this won't be the case if the values are relcalculated from the using the coordinate data block above - getting the right values depends on having the fill covariance matrix available]. |
Refinement. The data in this cif come from a Rietveld refinement against neutron powder data. Details of refinement are given in the paper. Checkcif output: 091_ALERT_1_A No Wavelength found in CIF - 0.71073 Ang Assumed Please Check The data were collected by the time of flight method 340_ALERT_3_B Low Bond Precision on C-C Bonds ··············· 0.0400 Ang. 201_ALERT_2_G Isotropic non-H Atoms in Main Residue(s) ······. 5 Why ? 210_ALERT_3_G No Anisotropic ADP's Found in CIF ············.. Please Check The structure was refined against medium resolution neutron powder data. 155_ALERT_4_C The Triclinic Unitcell is NOT Reduced ·········. Please Do ! 128_ALERT_4_G Suggested Setting for Input Space-Group I1 P1 Note The setting is chosen to show the relationship to the epsilon phase. The matrix relating the prediction and refinement cell settings is [-1 1 0] [ 0 -1 1] [ 1 0 1] 161_ALERT_4_C Missing or Zero su (esd) on x-coordinate for ··· N1 162_ALERT_4_C Missing or Zero su (esd) on y-coordinate for ··· N1 163_ALERT_4_C Missing or Zero su (esd) on z-coordinate for ··· N1 This atom was used to fix the origin. 723_ALERT_1_C Torsion Calc -147.21, Rep -149 (2) Dev··· 1.79 Sigma O1 -C1 -C2 -N1 1.555 1.555 1.555 1.555 1 741_ALERT_1_C Bond Calc 1.25 (11), Rep 1.25000 ······ Missing su O1 -C1 1.555 1.555 # 1 741_ALERT_1_C Bond Calc 1.25 (7), Rep 1.25000 ······ Missing su O2 -C1 1.555 1.555 # 2 ···plus numerous similar alerts The s.u.s on the diatnces and angles are all exactly 0 as they were constrained during refinement. 017_ALERT_1_G Check Consistency of Scattering Type D for H1 017_ALERT_1_G Check Consistency of Scattering Type D for H2 017_ALERT_1_G Check Consistency of Scattering Type D for H3 017_ALERT_1_G Check Consistency of Scattering Type D for H4 017_ALERT_1_G Check Consistency of Scattering Type D for H5 The sample was deuterated. 154_ALERT_1_G The su's on the Cell Angles are Equal ·········. 0.03000 Deg. 720_ALERT_4_G Number of Unusual/Non-Standard Labels ·········. 5 802_ALERT_4_G CIF Input Record(s) with more than 80 Characters ! Info No action taken |
x | y | z | Biso*/Beq | ||
N1 | 0.71144 | 0.68132 | 0.77926 | 0.7 (2) | |
C1 | 0.527 (13) | 0.767 (13) | 1.112 (5) | 0.7 (2) | |
C2 | 0.500 (11) | 0.784 (11) | 0.821 (5) | 0.7 (2) | |
O1 | 0.297 (16) | 0.75 (2) | 1.148 (7) | 0.7 (2) | |
O2 | 0.778 (11) | 0.767 (11) | 1.293 (3) | 0.7 (2) | |
H1 | 0.277 (7) | 0.703 (18) | 0.686 (5) | 0.8 (3) | |
H2 | 0.54 (2) | 0.955 (11) | 0.772 (8) | 0.8 (3) | |
H3 | 0.731 (8) | 0.737 (7) | 0.599 (5) | 1.0 (3) | |
H4 | 0.633 (12) | 0.5123 (11) | 0.766 (10) | 1.0 (3) | |
H5 | 0.917 (3) | 0.724 (9) | 0.943 (5) | 1.0 (3) |
O1—C1 | 1.2500 | N1—H4 | 1.0500 |
O2—C1 | 1.2500 | N1—H5 | 1.0500 |
N1—C2 | 1.4800 | C2—H1 | 1.0900 |
C1—C2 | 1.5300 | C2—H2 | 1.0900 |
N1—H3 | 1.0500 | ||
O1—C1—O2 | 126.00 | H3—N1—H5 | 110.00 |
O1—C1—C2 | 117.00 | N1—C2—C1 | 112.00 |
O2—C1—C2 | 117.00 | N1—C2—H1 | 108.00 |
C2—N1—H5 | 109.00 | N1—C2—H2 | 108.00 |
C2—N1—H3 | 109.00 | C1—C2—H1 | 109.00 |
C2—N1—H4 | 109.00 | C1—C2—H2 | 109.00 |
H4—N1—H5 | 110.00 | H1—C2—H2 | 110.00 |
H3—N1—H4 | 110.00 | ||
O1—C1—C2—N1 | −149 (2) | O2—C1—C2—N1 | 31 (2) |
Acknowledgements
We thank Professor Qiang Zhu of the University of Nevada Las Vegas who put our two groups in contact when he learned of our shared interest in ζ-glycine, and Drs L. Martin-Samos, L. Paulatto and F. Pellegrini for suggestions on improving the manuscript. We also thank Dr Andrew Alexander for assistance with crystal growth of γ-glycine-d5.
Funding information
The following funding is acknowledged: Science and Technology Facilities Council/ISIS (award No. RB1520301); Engineering and Physical Sciences Research Council (award No. DTA); SISSA, CINECA and PRACE (award No. 2011050736); H2020 European Research Council (award No. 676531).
References
Berland, K., Cooper, V. R., Lee, K., Schröder, E., Thonhauser, T., Hyldgaard, P. & Lundqvist, B. I. (2015). Rep. Prog. Phys. 78, 066501. Web of Science CrossRef PubMed Google Scholar
Boldyreva, E. V., Ivashevskaya, S. N., Sowa, H., Ahsbahs, H. & Weber, H.-P. (2005). Z. Kristallogr. Cryst. Mater. 220, 50–57. CAS Google Scholar
Bordallo, H. N., Boldyreva, E. V., Buchsteiner, A., Koza, M. M. & Landsgesell, S. (2008). J. Phys. Chem. B, 112, 8748–8759. Web of Science CrossRef PubMed CAS Google Scholar
Brandenburg, K. (2004). DIAMOND. Crystal Impact GbR, Bonn, Germany. Google Scholar
Bull, C. L., Funnell, N. P., Tucker, M. G., Hull, S., Francis, D. J. & Marshall, W. G. (2016). High. Press. Res. 36, 493–511. Web of Science CrossRef Google Scholar
Chisholm, J. A., Motherwell, S., Tulip, P. R., Parsons, S. & Clark, S. J. (2005). Cryst. Growth Des. 5, 1437–1442. Web of Science CrossRef CAS Google Scholar
Coelho, A. (2015). TOPAS-Academic. Version 5. Coelho Software, Brisbane, Australia. Google Scholar
Dawson, A., Allan, D. R., Belmonte, S. A., Clark, S. J., David, W. I., McGregor, P. A., Parsons, S., Pulham, C. R. & Sawyer, L. (2005). Cryst. Growth Des. 5, 1415–1427. Web of Science CSD CrossRef CAS Google Scholar
Destro, R., Roversi, P., Barzaghi, M. & Marsh, R. E. (2000). J. Phys. Chem. A, 104, 1047–1054. Web of Science CSD CrossRef CAS Google Scholar
Fortes, A. D., Wood, I. G., Alfredsson, M., Vocadlo, L., Knight, K. S., Marshall, W. G., Tucker, M. G. & Fernandez-Alonso, F. (2007). High Pressure Res. 27, 201–212. Web of Science CrossRef CAS Google Scholar
Fortes, A. D., Wood, I. G., Alfredsson, M., Vocadlo, L., Knight, K. S., Marshall, W. G., Tucker, M. G. & Fernandez-Alonso, F. (2012). High Pressure Res. 32, 337–337. Google Scholar
Frisch, M. J. et al. (2009). GAUSSIAN09. Gaussian Inc., Wallingford, Connecticut, USA. Google Scholar
Gavezzotti, A. (2005). Z. Kristallogr. 220, 499–510. Web of Science CrossRef CAS Google Scholar
Gavezzotti, A. (2011). New J. Chem. 35, 1360–1368. Web of Science CrossRef CAS Google Scholar
Giannozzi, P. et al. (2009). J. Phys. Condens. Matter, 21, 395502. Web of Science CrossRef PubMed Google Scholar
Goryainov, S. V., Boldyreva, E. V. & Kolesnik, E. N. (2006). Chem. Phys. Lett. 419, 496–500. Web of Science CrossRef CAS Google Scholar
Goryainov, S. V., Kolesnik, E. N. & Boldyreva, E. V. (2005). Phys. B Condens. Matter, 357, 340–347. Web of Science CrossRef CAS Google Scholar
Hohenberg, P. & Kohn, W. (1964). Phys. Rev. 136, B864–B871. CrossRef Web of Science Google Scholar
Hohenstein, E. G. & Sherrill, C. D. (2010a). J. Chem. Phys. 133, 014101. Web of Science CrossRef PubMed Google Scholar
Hohenstein, E. G. & Sherrill, C. D. (2010b). J. Chem. Phys. 133, 104107. Web of Science CrossRef PubMed Google Scholar
Hohenstein, E. G. & Sherrill, C. D. (2012). WIREs Comput. Mol. Sci. 2, 304–326. Web of Science CrossRef CAS Google Scholar
Hunter, C. A. & Prohens, R. (2017). CrystEngComm, 19, 23–26. Web of Science CrossRef CAS Google Scholar
Jeziorski, B., Moszynski, R. & Szalewicz, K. (1994). Chem. Rev. 94, 1887–1930. CrossRef CAS Web of Science Google Scholar
Küçükbenli, E., Monni, M., Adetunji, B. I., Ge, X., Adebayo, G. A., Marzari, N., de Gironcoli, S. & Dal Corso, A. (2014). arXiv, 1404.3015. Google Scholar
Kvick, Å., Canning, W. M., Koetzle, T. F. & Williams, G. J. B. (1980). Acta Cryst. B36, 115–120. CSD CrossRef CAS IUCr Journals Web of Science Google Scholar
Lund, A. M., Pagola, G. I., Orendt, A. M., Ferraro, M. B. & Facelli, J. C. (2015). Chem. Phys. Lett. 626, 20–24. Web of Science CrossRef CAS PubMed Google Scholar
Lyakhov, A. O., Oganov, A. R., Stokes, H. T. & Zhu, Q. (2013). Comput. Phys. Commun. 184, 1172–1182. Web of Science CrossRef CAS Google Scholar
Macrae, C. F., Bruno, I. J., Chisholm, J. A., Edgington, P. R., McCabe, P., Pidcock, E., Rodriguez-Monge, L., Taylor, R., van de Streek, J. & Wood, P. A. (2008). J. Appl. Cryst. 41, 466–470. Web of Science CrossRef CAS IUCr Journals Google Scholar
Marom, N., DiStasio, R. A., Atalla, V., Levchenko, S., Reilly, A. M., Chelikowsky, J. R., Leiserowitz, L. & Tkatchenko, A. (2013). Angew. Chem. Int. Ed. 52, 6629–6632. Web of Science CrossRef CAS Google Scholar
Moggach, S. A., Marshall, W. G., Rogers, D. M. & Parsons, S. (2015). CrystEngComm, 17, 5315–5328. Web of Science CSD CrossRef CAS Google Scholar
Murli, C., Sharma, S. M., Karmakar, S. & Sikka, S. K. (2003). Physica B, 339, 23–30. Web of Science CrossRef CAS Google Scholar
Oganov, A. R. & Valle, M. (2009). J. Chem. Phys. 130, 104504–104512. Web of Science CrossRef PubMed Google Scholar
Perlovich, G., Hansen, L. & Bauer-Brandl, A. (2001). J. Therm. Anal. Calorim. 66, 699–715. Web of Science CrossRef CAS Google Scholar
Qi, S., Avalle, P., Saklatvala, R. & Craig, D. Q. M. (2008). Eur. J. Pharm. Biopharm. 69, 364–371. Web of Science CrossRef PubMed CAS Google Scholar
Reilly, A. M. et al. (2016). Acta Cryst. B72, 439–459. Web of Science CrossRef IUCr Journals Google Scholar
Rossi, A., Savioli, A., Bini, M., Capsoni, D., Massarotti, V., Bettini, R., Gazzaniga, A., Sangalli, M. E. & Giordano, F. (2003). Thermochim. Acta, 406, 55–67. Web of Science CrossRef CAS Google Scholar
Sabatini, R., Küçükbenli, E., Kolb, B., Thonhauser, T. & de Gironcoli, S. (2012). J. Phys. Condens. Matter, 24, 424209. Web of Science CrossRef PubMed Google Scholar
Sheldrick, G. M. (2001). SHELXTL-XP. Version 6.01. University of Göttingen, Germany. Google Scholar
Spek, A. L. (2003). J. Appl. Cryst. 36, 7–13. Web of Science CrossRef CAS IUCr Journals Google Scholar
Stone, A. J. (2013). The Theory of Intermolecular Forces. Oxford University Press. Google Scholar
Turney, J. M. et al. (2012). WIREs Comput. Mol. Sci. 2, 556–565. Web of Science CrossRef CAS Google Scholar
Woodley, S. M. & Catlow, R. (2008). Nat. Mater. 7, 937–946. Web of Science CrossRef PubMed CAS Google Scholar
Xu, W., Zhu, Q. & Hu, C. T. (2017). Angew. Chem. 129, 2062–2066. CrossRef Google Scholar
Zaccaro, J., Matic, J., Myerson, A. S. & Garetz, B. A. (2001). Cryst. Growth Des. 1, 5–8. Web of Science CrossRef CAS Google Scholar
Zhu, Q., Oganov, A. R., Glass, C. W. & Stokes, H. T. (2012). Acta Cryst. B68, 215–226. Web of Science CrossRef IUCr Journals Google Scholar
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.