Intermolecular anharmonicity in molecular crystals: Interplay between experimental low-frequency dynamics and quantum quasi-harmonic simulations of solid purine

The intermolecular anharmonic potential of crystalline purine is probed by means of temperature-dependent terahertz time-domain spectroscopy, low-frequency Raman scattering, X-ray diffraction, and ab initio quasi-harmonic quantum-chemical simulations. As temperature increases, anharmonicity in the intermolecular interactions results in strongly anisotropic thermal expansion - with a negative thermal expansion along the b crystallographic axis - yielding corresponding bulk structural modiﬁcations. The observed thermally-induced shifts of most vibrational bands in the terahertz region of the spectra are shown to arise from volume-dependent thermal changes of the hydrogen-bond pattern along the a and b crystallographic axes.

2][3] Crystalline materials are often described from two distinct perspectives, either in their static structural or in their dynamic vibrational aspects. 4,57][8] From a fundamental standpoint both phenomena are a result of the shape of the potential energy surface (PES), with the PES minima and curvature governing equilibrium structures and dynamics, respectively.Experimentally, it is possible to characterise the PES-dependent effects using a variety of techniques, for example X-ray diffraction for structures and vibrational spectroscopy for dynamics. 9n particular, when molecular crystals are considered, the characterisation of their relatively flat PES proves particularly delicate and yet is key to the understanding of their structural, mechanical, and thermodynamic properties, which govern their rich polymorphism, for example. 10In this respect, lowfrequency vibrational methods such as terahertz time-domain spectroscopy (THz-TDS) and Raman scattering would constitute powerful tools for probing both the intra-and inter-molecular PES of such systems, since the relevant vibrations mostly involve motions of entire molecules, rather than only individual bonds. 11,12However, the interpretation of the low-frequency spectra is far more complicated than it is for their mid-infrared counterparts, since the corresponding vibrations strongly depend on the intermolecular interactions, thus implying that each material has a unique spectral fingerprint at THz frequencies. 13,14][20] Further complicating the analysis of molecular crystals is the presence of anharmonicity in the PES, a common phenomenon that strongly influences the bulk properties of materials, especially those that depend on temperature. 21,22While this effect is almost always present in solids it becomes especially pronounced in the case of weakly-bonded molecular crystals. 23This leads to the properties of molecular crystals being significantly influenced even by modest temperature changes close to ambient conditions; understanding their thermal properties is thus of critical importance for most potential applications. 246][27] An effective way of accounting for some anharmonic effects of solids (such as thermal expansion, temperature-dependent elastic response, etc.) within standard DFT is represented by the quasi-harmonic approximation (QHA), which requires the evaluation of lattice dynamics at different volumes and thus enables the description of volume-dependent

View Article Online
View Journal | View Issue thermal features. 28,29Only very recently has this methodology been used to investigate the thermal properties of molecular crystals, [30][31][32][33] and is here applied for the first time to interpret the temperature-dependent low-frequency spectral properties of a molecular solid.
Recently, it was shown that crystalline purine, an important component of DNA bases, exhibits significant changes in its THz-TDS spectrum as a function of temperature, strongly indicating that there exists significant intermolecular anharmonicity. 34,35he origin of such strong anharmonic features is explored here, by performing variable temperature THz-TDS measurements, low-frequency Raman scattering experiments, and quantummechanical quasi-harmonic simulations.Samples of crystalline purine were used as received (Sigma-Aldrich, 498%) and were either mixed with high-density polyethylene (10% w/w) for THz-TDS or simply ground for the Raman measurements.The THz-TDS data were acquired on a commercial TeraView Terapulse 4000 spectrometer from 20-90 cm À1 , and the lowfrequency Raman data utilized a Horiba LabRAM HR Evolution spectrometer operating at 532 nm coupled with an ultra-low frequency accessory from 25-150 cm À1 , each having a spectral resolution of 1 cm À1 .The results of these experiments, shown in Fig. 1, clearly demonstrate that, with the exception of the feature occurring at 46.5 cm À1 in both spectra, all of the peaks sharpen and shift to higher frequencies upon cooling to great effect.
Solid-state DFT simulations were performed in order to assign the spectra and determine the nature of the vibrational features.A developmental version of the Crystal14 software package was used for all simulations. 36The atoms were modeled using the all-electron def2-SVP basis set 37 and the Perdew-Burke-Ernzerhof (PBE) density functional, 38 and corrected for dispersion forces by means of the D3 correction by Grimme. 39e crystal structure of purine was fully-optimised with no constraints other than those given by the orthorhombic space group it belongs to (Pna2 1 ) prior to the simulation of vibrational properties, and the IR and Raman intensities were calculated analytically via a Coupled-Perturbed Kohn-Sham (CPKS) method, [41][42][43] both to a convergence of DE o 10 À10 Hartree (additional details available in the ESI †).The results of the vibrational simulation are shown in Fig. 1.As previously reported, the simulation slightly overestimates the low-frequency vibrational features, requiring a scaling of 0.85 in order to come into agreement with the experiment, and this factor was used for all predicted vibrational frequencies discussed here. 355][46][47][48] Vibrational frequencies are computed at four volumes (from a À2% compression to a +4% expansion with respect to the static equilibrium volume) after full volumeconstrained structural relaxations, which provide an accurate description of anisotropic thermally-induced structural changes.The results of the QHA simulations constitute a means of correlating structural modifications with temperature and thus of predicting thermally-induced shifts of spectral peaks as well.The computed frequencies and shifts (consistently scaled by 0.85) are reported in Table 1, along with the corresponding experimental IR and Raman ones.The absolute magnitude of the measured shifts between 300 K and 100 K (0.6, 1.6, 1.7, 2.2, 2.5, 5.1, 9.0 and 9.1 cm À1 ) is such to represent a formidable challenge to the DFT simulations in order to be reliably reproduced.Indeed, let us stress that a shift by 1 cm À1 corresponds to an energy difference of 0.003 kcal mol À1 (more than 300 times smaller than the very much sought for ''chemical accuracy'' in the theoretical chemistry field). 39,40Very tight numerical parameters, combined with the robustness of the adopted algorithms, have been used (energy convergence threshold for each atomic displacement set to 10 À10 Hartree and a reciprocal space sampling over 216 k Fig. 1 Variable temperature THz-TDS and low-frequency Raman of crystalline purine from 100 K to 300 K in 50 K intervals.The bottom panels show the simulated vibrational spectra (scaled by 0.85) as determined by equilibrium solid-state DFT calculations compared to the 100 K experimental spectra.
Table 1 Experimental (IR and Raman) and theoretical peak positions of low-frequency spectral features of crystalline purine as measured and simulated, through the QHA, at two temperatures (100 K and 300 K).The corresponding thermal shifts are also reported (100-300 K) points in the irreducible Brillouin zone) to attempt the first such characterisation.All of the vibrational frequencies are blueshifted upon cooling, which is reproduced by the QHA simulations.The predicted shifts match the experimentally-measured ones to a large extent for most peaks, particularly so in terms of relative amplitudes, which is a remarkable achievement if the required numerical accuracy is considered, as discussed above.A partial exception to this agreement are the IR modes occurring at 46.5 and 75.9 cm À1 in the 100 K spectrum, which have their shifts underestimated and overestimated (as well as their predicted intensities), respectively.Both modes are indeed believed to exhibit some intrinsic vibrational anharmonicity (not accounted for at the QHA level of theory); previous studies have reported that both features continue to sharpen and shift significantly upon cooling to 4 K, unlike the other absorptions that show little change below 100 K. 34,35 Further accentuating the intrinsic vibrational anharmonicity of the 75.9 cm À1 feature is the asymmetry of the vibrational band compared to the others, which is a strong indicator of this particular phenomenon. 49he QHA frequency shifts can be correlated with modifications of particular structural groups, by examination of the corresponding normal mode displacement vectors.Low-frequency vibrational features of solid purine are dominated by motion of the intermolecular hydrogen bonds, lying in the (ab) plane.Indeed, there is a clear correlation between the relative amount of frequency shifting and hydrogen bond perturbation, particularly so in the N-HÁ Á ÁH bond angle rather than bond length (Fig. 2).Both Table 1 and Fig. 2 confirm that a very large fraction of the observed temperature-dependent shifts of the low-frequency intermolecular spectral features of solid purine are due to thermally-induced structural changes rather than to intrinsic vibrational anharmonicity.A deeper investigation of the crystal structure of purine as a function of temperature is then needed, both experimentally and theoretically.
Quasi-harmonic lattice dynamics represents an effective way to determine the temperature dependence of structural and elastic properties of solids.Volumetric thermal expansion and temperature-dependent bulk modulus of solid purine are reported in the two bottom-right panels of Fig. 3.The sole inclusion of the zero-point vibrational effect results in a volume increase of 0.8% and a bulk modulus decrease of 10%, with respect to the static 0 K equilibrium DFT values.Thermal effects at 300 K further increase the volume by 1.2% and decrease the bulk modulus by 12% (2% and 20% with respect to the static 0 K DFT values, respectively).From these figures, it clearly emerges that even slight structural changes are dramatically affecting the mechanical properties of solid purine.In order to better highlight how the structure is affected by temperature, the thermal expansion of purine is reported in panel A of Fig. 3, as determined from the QHA simulation, which is found to be rather anisotropic with a large expansion along the c axis (where molecules are mainly held together by weak dispersive interactions), a smaller expansion along the a axis, and a small negative thermal expansion along the b axis (where stronger hydrogen-bond interactions are taking place).
From an experimental point of view, single crystals of purine were grown via sublimation of the as-received sample, and SC-XRD data obtained at three temperatures (100 K, 200 K, and 300 K) 50, 51 (additional details provided in the ESI †).The data confirm the overall QHA picture, with small structural changes to a and b (with the b axis showing a slight contraction), and with the c axis exhibiting the greatest amount of expansion.As regards absolute values, theoretical lattice parameters at 100 K turn out to be slightly underestimated by À0.27%, À2.15%, and À0.63% in the a, b, and c axes, respectively.The overall computed contraction of the unit cell may partially explain the need for the applied scaling factor to the vibrational frequencies.
The SC-XRD and QHA data enable a precise atomic-level analysis of the temperature-dependent structure of solid purine.It is clear that small thermally-induced structural changes largely affect the strength of the intermolecular interactions, which is reflected in large mechanical softening and spectral peak shifting.Intermolecular interactions along the a axis produce no relevant structural modifications other than an expansion, which linearly follows the corresponding lattice parameter.The analysis of the relative orientation of the molecules with respect to one another (for example the angle between the hydrogenbonded molecular sheets), rather than simply of the distances between them, provides insight.The angle formed by the molecules and the (ab) plane, (001) Miller plane, opens by 1.771 upon cooling, while the inter-chain distance simultaneously contracts (Fig. 4).Here, the contraction of the hydrogen bond is negated by the change in the packing angle, leading to the observed negative thermal expansion behaviour.
This work shows how some anharmonic thermal features of molecular crystals can effectively be understood at an atomistic level by coupling quasi-harmonic lattice-dynamical simulations with experimental THz-TDS, Raman scattering, and single crystal X-ray diffraction.The success of the present investigation suggests that quasi-harmonic simulations could provide a reliable description of the thermo-mechanical and thermodynamic properties of this class of weakly-bound solids.M. T. R. and J. A. Z. thank the UK Engineering and Physical Sciences Research Council for funding (EP/N022769/1).M. T. R. also thanks the European Molecular Biology Organization for travel funding.We also thank Sara Dampf and Adam Zaczek for their assistance with crystallisation.

Fig. 2
Fig.2The change in amine hydrogen bond angle (top) and bond distance (bottom) for each low-frequency normal mode determined from the equilibrium DFT simulation (solid line).The amount of QHA thermallyinduced shift is also shown (dashed line).Experimentally active vibrational bands are marked with blue stars and red diamonds for IR and Raman, respectively.

Fig. 3
Fig. 3 Results of the QHA simulations, showing (A) the thermal expansion coefficient of each crystal axis, (B, top) the predicted lattice constants as a function of temperature (with experimental values presented as points), (B, middle) the predicted volume as a function of temperature, (B, bottom) the bulk modulus as a function of temperature, and (C) the optimised 0 K crystal structure.