- Research article
- Open Access
The effect of temperature and pressure on the crystal structure of piperidine
© Budd et al.; licensee Springer. 2015
- Received: 29 October 2014
- Accepted: 5 February 2015
- Published: 12 April 2015
The response of molecular crystal structures to changes in externally applied conditions such as temperature and pressure are the result of a complex balance between strong intramolecular bonding, medium strength intermolecular interactions such as hydrogen bonds, and weaker intermolecular van der Waals contacts. At high pressure the additional thermodynamic requirement to fill space efficiently becomes increasingly important.
The crystal structure of piperidine-d11 has been determined at 2 K and at room temperature at pressures between 0.22 and 1.09 GPa. Unit cell dimensions have been determined between 2 and 255 K, and at pressures up to 2.77 GPa at room temperature. All measurements were made using neutron powder diffraction. The crystal structure features chains of molecules formed by NH…N H-bonds with van der Waals interactions between the chains. Although the H-bonds are the strongest intermolecular contacts, the majority of the sublimation enthalpy may be ascribed to weaker but more numerous van der Waals interactions.
Analysis of the thermal expansion data in the light of phonon frequencies determined in periodic DFT calculations indicates that the expansion at very low temperature is governed by external lattice modes, but above 100 K the influence of intramolecular ring-flexing modes also becomes significant. The principal directions of thermal expansion are determined by the sensitivity of different van der Waals interactions to changes in distance. The principal values of the strain developed on application of pressure are similarly oriented to those determined in the variable-temperature study, but more isotropic because of the need to minimise volume by filling interstitial voids at elevated pressure.
- Unit Cell Volume
- Unit Cell Dimension
- Sublimation Enthalpy
In cyclohexane-II the molecules are distributed in layers which stack in the ABCABC… pattern characteristic of cubic close packing. In piperidine and piperazine the NH groups are capable of forming hydrogen bonds. In piperazine the arrangement of the molecular centroids is the same as in cyclohexane-II, but the molecules are rotated to enable H-bonds to form between the layers . The structure of piperidine can also be considered to consist of layers, but the distribution of the molecules with respect to inversion centres in the space group makes the topology hexagonally close-packed. H-bonds form only between alternate layers, and the structure can be considered to be a hybrid of cyclohexane and piperazine. The rotations necessary for H-bond formation optimise interactions between the NH groups, but at the expense of breaking the van der Waals interactions formed within the layers of cyclohexane, and piperidine has a somewhat lower packing coefficient  (0.66) than either cyclohexane or piperazine (0.71).
The low packing coefficient suggests that piperidine may be susceptible to formation of different phases. This suggestion is supported by differential scanning calorimetry data presented by Parkin , which show that piperidine appears to form one phase on cooling from the liquid, which then transforms to another phase at 239 K. However, the higher temperature phase is not recovered on warming. In addition, during preliminary crystal growth experiments Parkin et al. obtained a weakly-diffracting form of piperidine with a cell volume approximately half of that of the phase they eventually characterised. The cell dimensions of this smaller-volume phase were a = 7.033(3), b = 5.224(3), c = 7.852(4) Å and β = 108.03(3)°. However, despite numerous attempts in which the cooling rate and sample size were varied this phase was never observed again.
The primary aim of the present investigation was therefore to determine if new forms of piperidine (in its perdeuterated form) could be identified under different conditions of temperature and pressure. The data collected also enable a detailed analysis of the thermodynamic properties of piperidine which are analysed with the aid of periodic DFT simulations and PIXEL packing energy calculations.
Piperidine-d11 was obtained from Aldrich and used as received.
Variable temperature measurements
Variable temperature time-of-flight neutron powder diffraction data were recorded using the HRPD instrument at ISIS . As piperidine is liquid under ambient conditions, piperidine-d11 was poured into a liquid nitrogen chilled stainless steel mortar  and cold-ground before being loaded into a rectangular aluminium sample can fitted with a heater. The sample was then placed in a cryostat held at 100 K. After confirming the sample was in the known phase, the temperature was reduced to 2 K and data collected at this temperature. Data were then collected in 5 K steps from 5 K to 250 K; there was no indication of a phase transition in this range. No evidence of a pre-melting transition was found on increasing the temperature further to 255 K.
In order to investigate the possibility of a second phase close to the melting point, the sample was removed from the cryostat and allowed to melt. It was then transferred to a cylindrical vanadium can containing glass wool (to promote formation of a randomly-oriented powder) before being placed in a cryostat held at 245 K. This should have been well above the transition temperature (239 K) reported in ref. , however upon crystallisation the known low temperature phase was obtained.
Quenching a similar sample of piperidine in liquid nitrogen also resulted in the sample crystallising into the same phase. No evidence for a second phase of piperidine was observed.
High pressure measurements
Ambient temperature, high pressure neutron powder diffraction data were collected using the time-of-flight technique on the PEARL beamline high pressure facility at ISIS [7,8]. Piperidine-d11 was contained in a null-scattering Ti-Zr alloy capsule gasket  and loaded into a Paris-Edinburgh press [10-12]. The sample was loaded with powdered silica wool as previous work has shown it to aid formation of a well-oriented powder when crystallising liquids in situ; a small pellet of lead was included as a pressure marker. The pressure was calculated from the refined lead cell parameter using a Birch-Murnaghan equation of state [13,14] with V 0 = 30.3128 Å3, B = 41.92 GPa, B’ = 5.72. These parameters were derived by Fortes [15,16] by refitting data obtained in three earlier studies [17-19].
On increasing the pressure to 0.22 GPa the sample crystallised in the monoclinic phase observed at low temperature. Diffraction data were collected at 0.22, 0.49, 0.80, 1.09 and 1.36 GPa for ca 700 μA h, with shorter (ca 100 μA h) collections at 2.06 and 2.77 GPa. Peak broadening was pronounced above 1.09 GPa, because of the lack of a hydrostatic medium, and therefore the sample was decompressed to approximately 0.23 GPa. Over a period of four hours the pressure decreased slightly resulting in the sample melting. Recompressing the sample to 0.30 GPa resulted in the sample recrystallising in the same phase.
Crystal and refinement data
Cell setting, space group
a, b, c (Å)
8.6994(17), 5.2552(9), 11.9045(16)
8.5969(12), 5.2010(7), 11.7936(12)
8.5150(14), 5.1577(8), 11.6988(14)
8.4452(15), 5.1204(9), 11.6181(16)
8.59695(4), 5.21506(2), 11.93271(4)
α, β, γ (°)
90, 96.468(17), 90
90, 96.507(14), 90
90, 96.532(15), 90
90, 96.558(17), 90
90, 96.8790(4), 90
Dcalc (g cm−3)
Time of flight
Time of flight
Time of flight
Time of flight
Time of flight
Fitted range of d (Å)
0.75 – 4.17
0.75 – 4.17
0.75 – 4.17
0.75 – 4.12
0.83 – 2.51
9 term Chebychev polynomial
6 term Chebychev polynomial
Back-to-back exponential convoluted with Voigt function
Number of parameters
During refinement of the data collected at 2 K the piperidine molecules were treated using the Z-matrix formalism in TOPAS-Academic. Bond lengths and angles were refined, though all C-D distances were constrained to be equal. Anisotropic displacement parameters were modelled using the TLS formalism [23,24]. In spite of cold-grinding it was necessary to include a fourth order spherical harmonic correction for preferred orientation. The Rietveld refinement profile is shown in Figure 1c with crystal and refinement data given in Table 1. Pawley fits  were carried out on the rapid-scan data to obtain the lattice parameters.
DFT calculations were performed on piperidine-d11 using the plane-wave pseudopotential method in the CASTEP code  as incorporated into Materials Studio version 7 . The PBE exchange-correlation functional  was used with norm-conserving pseudopotentials and a basis set cut-off energy of 950 eV. Brillouin zone integrations were performed with a Monkhorst-Pack  k-point grid spacing of 0.07 Å−1, corresponding to a grid of 2 × 3 × 1 points. These parameters gave an energy convergence of < 0.1 meV/atom.
For analysis of the effect of pressure and temperature on intermolecular interactions geometry optimisations were carried out while holding the unit cell parameters fixed to those observed experimentally. The total energy convergence tolerance was 5×10−6 eV/atom, with a maximum force tolerance of 0.01 eV Å−1 and a maximum displacement of 5×10−4 Å. The SCF convergence criterion was 5×10−7 eV/atom, and the space group symmetry was retained during geometry optimisation.
For frequency calculations, the geometry optimisation was carried out in two stages. The coordinates and unit cell dimensions were optimised using the Tkatchenko-Scheffler correction for dispersion (DFT-D) [30,31] starting from the experimentally determined structure determined at 2 K. The convergence criteria at this stage were the same as those given above with the additional requirement of a maximum stress tolerance of 0.02 GPa. The optimised cell parameters and volume were: a = 8.537417 [expt at 2 K: 8.59695(4)] Å, b = 5.169532 [5.21506(2)] Å, c = 11.923783 [11.93271(4)] Å, β = 96.834764 [96.8790(4)]° and V = 522.509783 [531.14(1)] Å3, reproducing the cell dimensions to within 1% and the volume to 1.6%. Linear response  calculations are not yet possible in CASTEP with dispersion-corrected functionals, and so a second stage of optimisation was carried out with the pure PBE functional holding the cell dimensions fixed. At this stage the energy convergence tolerance was 5×10−7 eV/atom, with a maximum force tolerance of 0.002 eV Å−1 and a maximum displacement of 5×10−5 Å; the SCF convergence criterion was 5×10−10 eV/atom. The maximum force on any atom at convergence was 0.00012 eV Å−1. Phonon density-of-states and dispersion curves were calculated with Fourier interpolation . No imaginary frequencies were observed at any point in the Brillouin zone. Generation of structures distorted along phonon modes was accomplished with MODE_FOLLOW .
Electron densities were calculated using Gaussian09  at the MP2 level of theory with the 6-31G** basis set using molecular geometries derived from the crystal structures. The PIXEL method [36-40], as implemented in the program OPiX , was then used to calculate the intermolecular interaction energies.
Other programs used
Structures were visualised with XP , MERCURY  and DIAMOND . Strain tensor calculations were carried out using a locally written program , based on the discussion by Ohashi and Burnham [46,47] and employing the JACOBI subroutine Numerical Recipes . Equation-of-state calculations were performed using EOSfit . ORIGIN was used for fitting of the unit cell parameters to variable temperature unit cell data. Animated GIFs were created from images generated in MERCURY using the web-site http://www.gifmaker.me/#003.
The crystal structure of piperidine has been previously determined from X-ray data at 150 K. Variable temperature neutron powder diffraction revealed the structure exists between 2 K and 255 K. Attempts to access another phase by crystallisations via rapid cooling to 245 K and 77 K were not successful. Crystallisation via the application of pressure resulted in the same phase as obtained at low temperature. No phase transitions were evident upon compression to 2.77 GPa.
The crystal structure of piperidine
Intermolecular contact distances in piperidine-d 11 as determined by neutron powder diffraction, and molecule-molecule contact energies E estimated in PIXEL calculations
− E COUL / kJmol −1
− E POL / kJmol −1
− E DISP / kJmol −1
E REP / kJmol −1
− E TOT / kJmol −1
T or P
T = 2 K
P = 0.22 GPa
P = 1.09 GPa
Contacts formed within H-bonded chains
[−x + 1,y − ½, −z + ½]
[−x + 1,y + ½, −z + ½]
[x,y + 1,z]
[x,y − 1,z]
Contacts formed between H-bonded chains
[−x + 2,y − ½, −z + ½]
[−x + 2,y + ½, −z + ½]
[−x + 2, −y, −z + 1]
[−x + 1, −y, −z + 1]
[x, −y + ½,z + ½]
[x, −y + ½,z − ½]
[x, −y − ½,z + ½]
[x, −y − ½,z − ½]
The effect of temperature on the structure of piperidine
Fitting of volume versus temperature data using equation (5) yields values of θ, V(T 0) and X which can be interpreted in terms of vibrational, structural and heat capacity data.
Einstein model fitting parameters
X 0 /Å or Å 3
X 1 /Å or Å 3
X 2 /Å or Å 3
Θ 1 /K
Θ 2 /K
Structures of piperidine were optimised by periodic DFT using cell dimensions fixed to the experimentally-determined values at 2 K and in 25 K steps from 25 K to 250 K. The resulting structures were then combined into a movie, available in the Additional file 2: Movie S1, enabling the thermal expansion to be visualised. The H-bonded chains are arranged in rows which run parallel to c; these rows move apart as temperature is increased. Interestingly, the largest expansion component is the one approximately aligned with the axis (a) associated with the largest Einstein temperatures given in Table 3; similarly the smallest component lies approximately along c, which is associated with the smallest Einstein temperatures. In part, this may be the result of applying Equ. 5 to cell dimension data, when it strictly applies only to volume. However, the strain tensors components obtained from the variable temperature and pressure experiments adopt a similar alignment, suggesting that the anisotropy of thermal expansion is determined by the compressibility of the structure along different directions (see below).
The effect of pressure on the unit cell dimensions of piperidine
Comparison of eigenvalues and eigenvectors of the temperature and pressure induced strain tensors
T (255 K versus 2 K)
P (0.22 GPa versus 0.80 GPa)
[0.113 0.000 0.021]
[0.115 0.000 0.016]
[0.000 0.188 0.000]
[0.000 0.190 0.000]
[−0.016 0.000 0.081]
[−0.009 0.000 0.083]
The bulk modulus (B) refined for a third order Birch-Murnaghan equation of state  is 6.4(4) GPa, though the data-set used to calculate this quantity is admittedly rather limited. The values of V 0 and B’ refined to 557.0(15) Å3 and 8.2(8) respectively (χ2 = 0.74). Molecular solids typically have B < 30 GPa and the following B values are useful for comparison : Ru3(CO)12 6.6 GPa, L-alanine  and salicylaldoxime  13.3 GPa, NaCl 25 GPa, quartz 37 GPa, ceramics 50–300 GPa and diamond 440 GPa. Piperidine is thus quite soft compared to other hydrogen bonded materials.
The effect of pressure on the structure of piperidine
The neutron powder diffraction data obtained at elevated pressure were of medium resolution and modest statistical precision, and during structure refinement the piperidine molecules were treated as rigid bodies with bond distances and angles taken from the X-ray structure at 150 K. To test the validity of the rigid-body assumption, the experimental structures were optimised using periodic DFT, holding the unit cell dimensions fixed to experimental values. Such calculations have been shown to reproduce experimentally determined atomic positions with crystal packing similarity valuesa of < 0.1 Å . For piperidine, crystal packing similarities [43,66] of the experimental and optimised structures were between 0.021 and 0.037 Å, confirming that the rigid-body assumption is reasonable.
Over the course of the first one or two kbars of applied pressure, provided the symmetry remains unchanged and volume changes are small, the path of compression might be expected to follow the path of a low-energy lattice vibration. Specifically, the mode would be expected to lie at the Γ-point and belong to the totally symmetric representation, in the present case, where the lattice point symmetry is 2/m, A 1g. The lowest totally symmetric phonon, which has a calculated frequency of 36 cm−1 is shown as an animation in the Additional file 4: Movie S3a) and can be compared to a sequential animation of the structures obtained at 0.22 and 0.49 GPa (Additional file 4: Movie S3b). The similarity between the two animations is striking, suggesting that structure determinations at pressures of a few kbar can potentially be used to visualise low-frequency phonons.
Comparison of the effects of pressure and temperature
The strain induced by cooling or compression is characterised by similarly orientated strain tensors, with largest principal axis lying approximately along a, the middle lying exactly along b and the smallest approximately along c (Figure 3c). Though the principal axes are similarly orientated, the tensor describing the response to pressure is more isotropic.
In an attempt to gain some insight into the structural reasons for these similarities and differences the DFT-optimised structure at 0.22 GPa was reoptimised with first the a and then the c axis reduced in length by 5%. PIXEL calculations were then carried out on all three optimised structures, referred to as ‘opt’, ‘A’, and ‘C’ respectively. The lattice energies were −55.2, −53.8 and −53.1 kJ mol−1, respectively; though the differences between these figures are small, they are consistent with the greater compressibility along a. The H-bond energies (kJ mol−1) are −21.7 (opt), −21.4 (A) and −21.9 (C), indicating that compression along c actually stabilises the H-bonds, while compression along a destabilises them. The H-bond distances (Å) are opt 2.136, A 2.074 and C 2.111, and though in both A and C the electrostatic term is more negative as a result of the shorter N…H distance, it is outweighed by a more positive repulsion term in A.
We have explored the crystal structure of piperidine under a range of different conditions of temperature and pressure. Although a previous study had indicated another phase exists, no new phases have been identified here, and it is possible that the alternative phase identified by Parkin on the basis of cell dimensions was a metastable form. The crystal structure of piperidine is characterised by chains of molecules linked by NH…N hydrogen bonds which are disposed about a 21 axis running along b. The unit cell of the metastable phase has a b axis length very similar to that of the b-axis of the phase discussed in this paper, and both phases are likely to feature similar H-bond chains. Though our original crystal growth experiments and DSC measurements indicated formation of a metastable form, there is no mention of transitions in the report of heat capacity measurements in ref. , and it is possible that its formation is dependent on rather specific crystallisation conditions which remain unidentified. Deuteration is also known to influence formation of alternative phases .
Though the H-bonds have a characteristic and easily recognisable structural signature, the same is not true of van der Waals interactions , but their influence can be appreciated with the aid of packing energy calculations, such as the PIXEL method used here. In the case of piperidine, though H-bonds are the most energetic contacts, around 60% of the sublimation enthalpy is the result of van der Waals interactions in which individual interatomic distances teeter on the brink of insignificance when assessed in the conventional manner using van der Waals radii.
Periodic DFT calculations reproduce the experimental cell dimensions of piperidine at 2 K to within 1%. The phonon frequencies calculated using this optimised structure, though they are obtained on the basis of a harmonic model for the energies of small atomic displacements, also reproduce experimental heat capacity data well. These calculations enable the thermal expansion of piperidine to be understood in terms of a simple quantitative model for how energy in the form of heat is absorbed into lattice vibrations.
The response of piperidine to pressure is superficially similar to its response to cooling. Although the principal directions of strain are essentially the same under both sets of varying conditions, the tensor describing the strain developed under pressure is more isotropic. PIXEL calculations demonstrate that the energetic imbalance between the weak van der Waals interactions and the relatively strong H-bonds determines the anisotropy of the strain tensor on cooling. Under pressure the need to fill space efficiently competes with minimisation of free energy though intermolecular contacts, making the strain less anisotropic.
aThis is a parameter quantifying the similarity of two crystal structures, and it is defined as the root-mean square deviation calculated for a cluster of 15 molecules taken from each structure. Crystal packing similarities were calculated in MERCURY.
We thank STFC and ISIS through The Centre for Molecular Structure and Dynamics and EPSRC for studentship funding to LEB. We also thank the ISIS Facility for provision of neutron beam time on the HRPD and PEARL instruments. We also thank the referees for their careful reading of the manuscript.
- Kahn R, Fourme R, Andre D, Renaud M. Crystal structures of cyclohexane I and II. Acta Crystallogr Sec B: Struct Crystallogr Crystal Chem. 1973;29:131–8.View ArticleGoogle Scholar
- Funnell NP, Dove MT, Goodwin AL, Parsons S, Tucker MG. Local structure correlations in plastic cyclohexane-a reverse Monte Carlo study. J Phys Condens Matter. 2013;25:454204.View ArticleGoogle Scholar
- Parkin A, Oswald IDH, Parsons S. Structures of piperazine, piperidine and morpholine. Acta Crystallogr Sec B-Struct Sci. 2004;60:219–27.View ArticleGoogle Scholar
- Kitigorodski AI. Molecular crystals and molecules. New York: Academic Press; 1973.Google Scholar
- Ibberson RM. Design and performance of the new supermirror guide on HRPD at ISIS. Nucl Instrum Meth Phys Res Section A: Accelerators Spectrometers Detectors Assoc Equip. 2009;600:47–9.View ArticleGoogle Scholar
- Ibberson RM. A simple technique for preparing low-melting-point samples for neutron powder diffraction. J Appl Crystallogr. 1996;29:498–500.View ArticleGoogle Scholar
- ISIS. Dedicated facility for high pressure diffraction. In: ISIS 96 - ISIS facility annual report 1995–96, Rutherford Appleton laboratory 1996, RAL-TR-96-050. 1996. p. 61–2.Google Scholar
- ISIS. PEARL - pressure and engineering research line. In: ISIS 97 - ISIS facility annual report 1996–97, Rutherford Appleton laboratory 1997, RAL-TR-97-050. 1997. p. 28–9.Google Scholar
- Marshall WG, Francis DJ. Attainment of near-hydrostatic compression conditions using the Paris-Edinburgh cell. J Appl Crystallogr. 2002;35:122–5.View ArticleGoogle Scholar
- Besson JM, Nelmes RJ, Halem G, Loveday JS, Weill G, Hull S. Neutron powder diffraction above 10 GPa. Phys B Condens Matter. 1992;180–181:907–10.View ArticleGoogle Scholar
- Besson JM, Nelmes RJ. New developments in neutron-scattering methods under high pressure with the Paris-Edinburgh cells. Phys B Condens Matter. 1995;213&214:31–6.View ArticleGoogle Scholar
- Le Godec Y, Dove MT, Redfern SAT, Tucker MG, Marshall WG, Syfosse G, et al. Recent developments using the Paris-Edinburgh cell for neutron diffraction at high pressure and high temperature and some applications. High Pressure Res. 2003;23:281–7.View ArticleGoogle Scholar
- Angel RJ. Equations of state. Rev Mineral Geochem. 2002;41:35–59.View ArticleGoogle Scholar
- Birch F. Finite elastic strain of cubic crystals. Phys Rev. 1947;71:809–24.View ArticleGoogle Scholar
- Fortes AD, Wood IG, Alfredsson M, Vocadlo L, Knight KS, Marshall WG, et al. The high-pressure phase diagram of ammonia dihydrate. High Pressure Res. 2007;27:201–12.View ArticleGoogle Scholar
- Fortes AD, Wood IG, Alfredsson M, Vocadlo L, Knight KS, Marshall WG, et al. The high-pressure phase diagram of ammonia dihydrate (Corrigendum). High Pressure Res. 2012;32:337-337.Google Scholar
- Kuznetsov AZ, Dmitriev V, Dubrovinsky L, Prakapenka V, Weber HP. FCC-HCP phase boundary in lead. Solid State Commun. 2002;122:125–7.View ArticleGoogle Scholar
- Miller RA, Schuele DE. Pressure derivatives of the elastic constants of lead. J Phys Chem Solid. 1969;30:589–600.View ArticleGoogle Scholar
- Waldorf DL, Alers GA. Low-temperature elastic moduli of lead. J Appl Phys. 1962;33:3266–9.View ArticleGoogle Scholar
- Rietveld HM. Profile refinement method for nuclear and magnetic structures. J Appl Crystallogr. 1969;2:65–71.View ArticleGoogle Scholar
- Coelho A. TOPAS-academic: general profile and structure analysis software for powder diffraction data. Version 5. Brisbane, Australia; 2013.Google Scholar
- Jarvinen M. Application of symmetrized harmonics expansion to correction of the preferred orientation effect. J Appl Crystallogr. 1993;26:525–31.View ArticleGoogle Scholar
- Dunitz JD. X-ray analysis and the structure of organic molecules. Basel, Switzerland: VCH; 1995.View ArticleGoogle Scholar
- Schomaker V, Trueblood KN. Rigid-body motion of molecules in crystals. Acta Crystallogr Section B: Struct Crystallogr Crystal Chem. 1968;24:63–76.View ArticleGoogle Scholar
- Pawley GS. Unit-cell refinement from powder diffraction scans. J Appl Crystallogr. 1981;14:357–61.View ArticleGoogle Scholar
- Clark SJ, Segall MD, Pickard CJ, Hasnip PJ, Probert MJ, Refson K, et al. First principles methods using CASTEP. Z Krist. 2005;220:567–70.View ArticleGoogle Scholar
- Inc A. Materials studio version 7. Cambridge, UK: Materials Studio; 2013.Google Scholar
- Perdew JP, Burke K, Ernzerhof M. Generalized gradient approximation made simple. Phys Rev Lett. 1996;77:3865–8.View ArticleGoogle Scholar
- Monkhorst HJ, Pack JD. Special points for Brillouin-zone integrations. Phys Rev B. 1976;13:5188–92.View ArticleGoogle Scholar
- Tkatchenko A, Scheffler M. Accurate molecular van der Waals interactions from ground-state electron density and free-atom reference data. Phys Rev Lett. 2009;102:073005.View ArticleGoogle Scholar
- McNellis ER, Meyer J, Reuter K. Azobenzene at coinage metal surfaces: role of dispersive van der Waals interactions. Phys Rev B. 2009;80:205414.View ArticleGoogle Scholar
- Refson K, Tulip PR, Clark SJ. Variational density-functional perturbation theory for dielectrics and lattice dynamics. Phys Rev B. 2006;73:155114.View ArticleGoogle Scholar
- Gonze X, Charlier J-C, Teter MP. Interatomic force constants from first principles: the case of alpha-quartz. Phys Rev B. 1994;50:13035–8.View ArticleGoogle Scholar
- Refson K. MODE_FOLLOW. A routine for generating structures distorted along phonon modes calculated using CASTEP. Chilton, OXON, UK: ISIS Facility, Rutherford-Appleton Laboratory; 2014.Google Scholar
- Frisch GWT MJ, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Scalmani G, et al. Gaussian 09. Wallingford, USA: Gaussian Inc; 2009.Google Scholar
- Gavezzotti A. Calculation of intermolecular interaction energies by direct numerical integration over electron densities. I. Electrostatic and polarization energies in molecular crystals. J Phys Chem B. 2002;106:4145–54.View ArticleGoogle Scholar
- Gavezzotti A. Calculation of intermolecular interaction energies by direct numerical integration over electron densities. 2. An improved polarization model and the evaluation of dispersion and repulsion energies. J Phys Chem B. 2003;107:2344–53.View ArticleGoogle Scholar
- Gavezzotti A. Calculation of lattice energies of organic crystals: the PIXEL integration method in comparison with more traditional methods. Z Krist. 2005;220:499–510.View ArticleGoogle Scholar
- Gavezzotti A. Quantitative ranking of crystal packing modes by systematic calculations on potential energies and vibrational amplitudes of molecular dimers. J Chem Theory Comput. 2005;1:834–40.View ArticleGoogle Scholar
- Gavezzotti A. Molecular aggregation - structure analysis and molecular simulation of crystals and liquids. 1st ed. New York: Oxford University Press; 2007.Google Scholar
- Gavezzotti A. OPiX. Milan, Italy; 2003.Google Scholar
- Sheldrick GM. SHELXTL-XP version 6.01. Göttingen, Germany and Madison, Wisconsin, USA: Bruker AXS, Madison Wisc; 2001.Google Scholar
- Macrae CF, Bruno IJ, Chisholm JA, Edgington PR, McCabe P, Pidcock E, et al. Mercury CSD 2.0 - new features for the visualization and investigation of crystal structures. J Appl Cryst. 2008;41:466–70.View ArticleGoogle Scholar
- CrystalImpact. DIAMOND version 3.0, visual crystal structure information system. Bonn, Germany: CrystalImpact; 2004.Google Scholar
- Parsons S. STRAIN, program for calculation of linear strain tensors. 2003.Google Scholar
- Ohashi Y, Burnham CW. Clinopyroxene lattice deformations. Roles of chemical substitution and temperature. Am Mineral. 1973;58:843–9.Google Scholar
- Hazen RM, Finger LW. Comparative crystal chemistry: temperature, pressure, composition and the variation of crystal structure. Chichester, New York, USA: John Wiley and Sons; 1982. p. 81.Google Scholar
- Press WH, Teukolsky SA, Vetterling WT, Flannery BP. Numerical recipes in Fortran. Secondth ed. Cambridge, England: Cambridge University Press; 1992.Google Scholar
- Angel RJ, Gonzalez-Platas J, Alvaro M. EosFit7c and a Fortran module (library) for equation of state calculations. Z Krist. 2014;229:405–19.Google Scholar
- Bruno IJ, Cole JC, Edgington PR, Kessler M, Macrae CF, McCabe P, et al. New software for searching the Cambridge structural database and visualizing crystal structures. Acta Crystallogr Sec B: Struct Sci. 2002;B58:389–97.View ArticleGoogle Scholar
- Gould RO, Taylor P, Thorpe M. PUCKER, a program for conformational analysis of 5–8 membered rings. Edinburgh, UK: University of Edinburgh; 1995.Google Scholar
- Messerly JF, Todd SS, Finke HL, Good WD, Gammon BE. Condensed-phase heat-capacity studies and derived thermodynamic properties for six cyclic nitrogen compounds. J Chem Thermodyn. 1988;20:209–24.View ArticleGoogle Scholar
- Hossenlopp IA, Archer DG. Enthalpies of vaporization of piperidine and 1,2-dimethylbenzene; gas-phase isobaric heat capacities of piperidine. J Chem Thermodyn. 1988;20:1061–8.View ArticleGoogle Scholar
- Chickos JS, Acree Jr WE. Enthalpies of sublimation of organic and organometallic compounds. 1910–2001. J Phys Chem Ref Data. 2002;31:537–698.View ArticleGoogle Scholar
- Fortes AD, Wood IG, Vocadlo L, Brand HEA, Knight KS. Crystal structures and thermal expansion of α-MgSO4 and β-MgSO4 from 4.2 to 300 K by neutron powder diffraction. J Appl Crystallogr. 2007;40:761–70.View ArticleGoogle Scholar
- Knight KS. A neutron powder diffraction determination of the thermal expansion tensor of crocoite (PbCrO4) between 60 K and 290 K. Mineral Mag. 1996;60:963–72.View ArticleGoogle Scholar
- Schofield PF, Knight KS, Stretton IC. Thermal expansion of gypsum investigated by neutron powder diffraction. Am Mineral. 1996;81:847–51.Google Scholar
- Cemič L. Thermodynamics in mineral sciences. Berlin: Springer-Verlag; 2005.Google Scholar
- Wallace DC. Thermodynamics of crystals. Mineola, New York: Dover Publications; 1998.Google Scholar
- Schlenker JL, Gibbs GV, Boisen Jr MB. Thermal expansion coefficients for monoclinic crystals: a phenomenological approach. Am Mineral. 1975;60:828–33.Google Scholar
- Saunders GA, Yogurtcu YK, Macdonald JE, Pawley GS. The elastic behavior of orthorhombic sulfur under pressure. Proc R Soc London Series: Math Phys Eng Sci. 1986;407:325–42.View ArticleGoogle Scholar
- Slebodnick C, Zhao J, Angel R, Hanson BE, Song Y, Liu Z, et al. High pressure study of Ru3(CO)12 by X-ray diffraction, Raman, and infrared spectroscopy. Inorg Chem. 2004;43:5245–52.View ArticleGoogle Scholar
- Funnell NP, Dawson A, Francis D, Lennie AR, Marshall WG, Moggach SA, et al. The effect of pressure on the crystal structure of L-alanine. CrystEngComm. 2010;12:2573–83.View ArticleGoogle Scholar
- Wood PA, Forgan RS, Henderson D, Parsons S, Pidcock E, Tasker PA, et al. Effect of pressure on the crystal structure of salicylaldoxime-I, and the structure of salicylaldoxime-II at 5.93 GPa. Acta Crystallogr Sec B: Struct Sci. 2006;62:1099–111.View ArticleGoogle Scholar
- Binns J, Healy MR, Parsons S, Morrison CA. Assessing the performance of density functional theory in optimizing molecular crystal structure parameters. Acta Crystallogr B. 2014;70:259–67.View ArticleGoogle Scholar
- Childs SL, Wood PA, Rodriguez-Hornedo N, Reddy LS, Hardcastle KI. Analysis of 50 crystal structures containing carbamazepine using the materials module of mercury CSD. Cryst Growth Des. 2009;9:1869–88.View ArticleGoogle Scholar
- Crawford S, Kirchner MT, Blaser D, Boese R, David WIF, Dawson A, et al. Isotopic polymorphism in pyridine. Angewandte Chem Int Ed. 2009;48:755–7.View ArticleGoogle Scholar
- Gavezzotti A. The “sceptical chymist”: intermolecular doubts and paradoxes. CrystEngComm. 2013;15:4027–35.View ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.