research communications
Utilizing Hirshfeld surface calculations, non-covalent interaction (NCI) plots and the calculation of interaction energies in the analysis of molecular packing
aResearch Centre for Crystalline Materials, School of Science and Technology, Sunway University, 47500 Bandar Sunway, Selangor Darul Ehsan, Malaysia, and bDepartment of Physics, Bhavan's Sheth R. A. College of Science, Ahmedabad, Gujarat 380001, India
*Correspondence e-mail: edwardt@sunway.edu.my
The analysis of atom-to-atom and/or residue-to-residue contacts remains a favoured mode of analysing the molecular packing in crystals. In this contribution, additional tools are highlighted as methods for analysis in order to complement the `crystallographer's tool', PLATON [Spek (2009). Acta Cryst. D65, 148–155]. Thus, a brief outline of the procedures and what can be learned by using Crystal Explorer [Spackman & Jayatilaka (2009). CrystEngComm 11, 19–23] is presented. Attention is then directed towards evaluating the nature, i.e. attractive/weakly attractive/repulsive, of specific contacts employing NCIPLOT [Johnson et al. (2010). J. Am. Chem. Soc. 132, 6498–6506]. This is complemented by a discussion of the calculation of energy frameworks utilizing the latest version of Crystal Explorer. All the mentioned programs are free of charge and straightforward to use. More importantly, they complement each other to give a more complete picture of how molecules assemble in molecular crystals.
Keywords: Hirshfeld surface calculations; non-covalent interaction plots; interaction energies; molecular packing.
1. Introduction
A widely employed approach to describe the packing of molecular compounds in their crystals is based on describing specific atom-to-atom contacts, such as in conventional A—H⋯D hydrogen bonding. This analysis is often extended into the highly popular supramolecular synthon approach (Desiraju, 1995), whereby residue-to-residue contacts are evaluated as exemplified in the familiar eight-membered carboxylic acid synthon, i.e. {⋯HOCO}2. Often sharing the directionality, robustness and utility in molecular packing that characterizes hydrogen bonding is the very well documented phenomenon of halogen bonding (Cavallo et al., 2016). The electrostatic attraction between ostensibly two (partially) negatively charged entities in halogen bonding is ascribed to an anisotropic distribution of electron density around the halogen atom (X) in that at the tip of the C—X bond, there is an electron-deficient region, a so-called polar cap or σ-hole (Brinck et al., 1992; Murray et al., 2007); π-hole interactions rely on a similar concept (Bauzá et al., 2015). Such σ-hole considerations are now employed to rationalize (Kolář & Hobza, 2016) the very long-documented secondary bonding interactions (Alcock, 1972; Haiduc, 1997), more recently repackaged in terms of the participating atoms, e.g. tetrel bonding for interactions involving Group 14 elements (Bauzá et al., 2013), pnictogen (Group 15; Scheiner, 2013), chalcogen (Group 16; Wang et al., 2009) and even aerogen bonding, i.e. interactions involving noble gases (Bauzá & Frontera, 2015), but see Edwards et al. (2017). In molecular crystals of organic molecules, π-systems are well known to participate in identifiable points of contact between molecules by π–π stacking interactions (Janiak, 2000) and C—H⋯π contacts (Nishio, 2004). In the realm of coordination chemistry, chelate rings can also have π-character (Masui, 2001) and participate in analogous π–π interactions, where one or both of the interacting rings is a chelate ring (Malenov et al., 2017), and C—H⋯π(chelate) interactions (Sredojević et al., 2007; Tiekink & Zukerman-Schpector, 2011). Metals themselves can associate in the solid state as most famously demonstrated by gold which, owing to relativistic effects, has a significant propensity to form Au⋯Au (aurophilic) interactions that provide comparable energies of stabilization to their crystals as do conventional hydrogen-bonding interactions and, indeed, can be competitive with these (Schmidbaur, 2000; Schmidbaur & Schier, 2008; Tiekink, 2014). Gold can also form Au⋯π(arene) interactions in their crystals (Caracelli et al., 2013), as do many elements, including main-group elements in low oxidation states (Caracelli et al., 2016). In the case of the latter, it is the lone-pair of electrons on the heavy element that interacts with the π-system, a phenomenon that arises owing to the σ-hole present at the tip of the lone-pair of electrons, analogous to that discussed above. Far from being mere curiosities, interactions involving metals, chelate rings and secondary bonding interactions impart stabilization energies to molecular packing akin to conventional hydrogen bonding (Tiekink, 2017).
Confronted by a myriad of different types of identifiable points of contact between ions and molecules in their crystals – the list of interactions cited above is not exhaustive and further types are sure to be appreciated and documented in the coming years – describing molecular packing in detail can be a challenge. Ton Spek's program PLATON (Spek, 2009), downloadable free of charge for academic users from https://www.platonsoft.nl/spek/xraysoft/, is an excellent starting point for such an analysis as the geometric parameters characterizing close contacts, including for non-conventional contacts such as element(lone pair)⋯π(arene) and C—H⋯π(chelate ring) contacts, are evident in the output from PLATON. All that is required is the final, validated (CIF). An obvious limitation here is the application of distance criteria to determine the presence of a contact: interactions can and do extend beyond sums of van der Waals radii (Boese et al., 2001; Dance, 2003). Other tools are also freely available to gain further insight into the way molecules pack in their crystals. Two in particular form the focus of this contribution, namely Hirshfeld surface analysis and non-covalent interaction plots. As indicated below, these are easy-to-use programs and can provide complementary information useful for the study of molecular packing. The purpose of the present contribution is to highlight the use of Crystal Explorer (Spackman & Jayatilaka, 2009) and NCIPLOT (Johnson et al., 2010) in the analysis of crystals and to provide pointers to get the most out of these programs. The interested reader is referred to the original cited papers for more detailed information of the theory behind the different approaches described herein.
A caveat: the popularity and importance of identifying atom-to-atom/residue-to-residue contacts notwithstanding, the perennial question facing those trying to understand how and why molecular crystals form is nicely summarized by the `egg causality dilemma' – what came first, the chicken or the egg? In the present context, are the identified intermolecular interactions responsible for directing the way molecules assemble in crystals or are the identified intermolecular interactions the result of the formation of crystals? Developing this last point further, in a global molecular packing approach, molecules assemble to minimize free space so that protrusions (`bumps') in a molecule are accommodated by impressions (`craters') of symmetry-related molecules. It is also noted in this context of close packing considerations that around 83% of molecular compounds crystallize in one of six close-packing space groups (Allen, 2002). In salient comments underscoring the above are the observations by Dunitz & Gavezzotti (2009) when writing specifically about weak interactions involving hydrogen atoms, in effect, that these atoms have to be accommodated somewhere and are unlikely to adopt repulsive configurations; the same ideas apply equally to other intermolecular interactions. Whatever the origin of the intermolecular interactions revealed in crystals, their identification and analysis, especially in a systematic and thorough manner, is surely a worthwhile enterprise.
2. Hirshfeld surface analysis and two-dimensional fingerprint plots
2.1. Preamble
The analysis of calculated Hirshfeld surfaces has become an invaluable tool for crystallographers and crystal engineers alike as this provides additional insight into weak intermolecular interactions influential in the packing of molecules in crystals. A Hirshfeld surface is defined by the density weight function of the specific molecule of interest (i.e. the pro-molecule) over the same sum of density of its nearest neighbour (i.e. the pro-crystal), thereby resulting in a 0.5 arbitrary units isosurface, which is similar to that of a van der Waals surface but, unlike the latter, takes into consideration neighbouring molecules and hence provides information about intermolecular interactions (McKinnon et al., 2007; Spackman & Jayatilaka, 2009). The Hirshfeld surfaces can be mapped with different properties namely, dnorm, electrostatic potential, shape-index and curvedness. These are useful to accumulate additional information on weak intermolecular interactions. Crystal Explorer 17 (Turner et al., 2017) may be downloaded from https://crystalexplorer.scb.uwa.edu.au/downloads.html.
The Hirshfeld surfaces mapped over dnorm utilize the function of normalized distances de and di, where de and di are the distances from a given point on the surface to the nearest atom outside and inside, respectively. The blue, white and red colour conventions used for the dnorm-mapped Hirshfeld surfaces recognize the interatomic contacts as longer, at van der Waals separations and short interatomic contacts, respectively. The views of Hirshfeld surfaces mapped over the electrostatic potential obtained using the computational chemistry package Tonto (Jayatilaka & Grimwood, 2003), integrated into the Crystal Explorer 17 program, also enables the visualization of the donors and acceptors of intermolecular interactions through blue and red regions around the participating atoms corresponding to positive and negative electrostatic potential on the surface, respectively. Being a powerful quantum chemistry package for wave-function calculation and surface generation, Tonto can be used as an alternative to popular quantum chemistry packages, e.g. Gaussian16 (Frisch et al., 2106), and is available in Crystal Explorer 17. The package uses Hartree–Fock/DFT theory wave-function calculations based on the input CIF.
The useful measures of curvature, namely curvedness and shape-index, introduced by Koendrink (Koenderink, 1990; Koenderink & Doorn, 1992), provide further chemical insight into molecular packing. A surface with low curvedness designates a flat region and may be indicative of π–π stacking in the crystal. On the other hand, a Hirshfeld surface with high curvedness is highlighted as dark-blue edges, which is indicative of an absence of π–π stacking. The shape-index is a qualitative measure of shape and is sensitive to subtle changes in surface shape, particularly in a flat region. Two shape indices differing by sign represent complementary `bumps and hollows'. The blue bump-shape and shape-index > 1 belongs to the donor, and that representing a red hollow with index < 1 corresponds to the acceptor of an intermolecular interaction.
The two-dimensional fingerprint plot derived from a Hirshfeld surface (Spackman & McKinnon, 2002; McKinnon et al., 2004) provides a convenient visual summary of the frequency of each combination of de and di across the surface of a molecule. It is a highly useful method to summarize complex information contained in a crystal. The colour of each point corresponding to the relative area of a (de, di) pair is recognized as the contribution from different interatomic contacts: blue, green and red correspond to small, moderate and greatest contributions whereas an uncoloured region indicates no contribution to the Hirshfeld surface. A fingerprint plot delineated into specific interatomic contacts contains information related to specific intermolecular interactions.
To conduct the above calculations, one should employ the final validated Crystal Explorer 17; by default, the program will adjust the X—H bond lengths to their neutron-derived values. Specific examples of how each of the above can be applied in the analysis of molecular compounds follows.
as the input to2.2. Illustrative examples
The first example concerns the generation and interpretation of Hirshfeld surfaces calculated over dnorm. Fig. 1(a) shows the chemical structure of (4-nitrophenyl)methyl 2,3-dihydro-1H-pyrrole-1-carboxylate (C12H10N4O), (I), which was reported recently (Zukerman-Schpector, Soto-Monsalve et al., 2018). To view the characteristic red spots indicating specific points of contact in the crystal, the Hirshfeld surface mapped over dnorm was calculated with the default setting of arbitrary units range; rotation of the generated plot enables the identification regions of interest, e.g. Fig. 1(b) and (c). The red spots can be classified as bright, diminutive and faint to correlate (qualitatively) with the strength of intermolecular contact, i.e. as potential hydrogen bonds, weak interactions or short interatomic contacts.
In Fig. 1(b), the bright-red spots near the hydrogen (indicated with `1') and oxygen (`2') atoms indicate donors and acceptors of a potential C—H⋯O interaction. The diminutive-red spot near the nitro-oxygen atom (`3') represents its participation as an acceptor in a comparatively weak C—H⋯O contact (with a pyrrole-hydrogen atom). Additional faint-red spots arising from short interatomic contacts can be viewed by reducing the range of arbitrary units in the calculation by modifying the value of the negative arbitrary unit, as is apparent from Fig. 1(c) where additional red regions are highlighted as 4–6 (Zukerman-Schpector, Soto-Monsalve et al., 2018).
The Hirshfeld surface mapped over the calculated electrostatic potential for (I) is shown as the image in Fig. 2. Here, the blue and red regions around the different atoms correspond to positive and negative electrostatic potentials, respectively, e.g. red regions are apparent around the oxygen atoms participating in the C—H⋯O contacts mentioned above.
The Hirshfeld surfaces mapped with other properties like shape-index and curvedness can be employed to describe the effect of weak intermolecular interactions in a crystal, e.g. in the crystal of 3-[(1Z)-{2-[bis({[(2-methylphenyl)methyl]sulfanyl})methylidene]hydrazin-1-ylidene}methyl]benzene-1,2-diol (C23H22N2O2S2) (II), Fig. 3(a) (Yusof et al., 2018). As an example, the donor and the acceptors of intermolecular C—H⋯π contacts can be recognized as blue and red regions around the participating atoms on the Hirshfeld surfaces mapped over shape-index properties corresponding to C—H⋯π/π⋯H—C (often abbreviated as C⋯H/H⋯C) contacts, as shown in Fig. 3(b) for (II).
Thus far, the focus has been upon conventional interactions such as C—H⋯O and C—H⋯π. It is noted that Hirshfeld surface analysis is equally useful for identifying other, less-common points of contact between molecules. This point is exemplified in the crystal of the halide-rich salt 2-{[2,8-bis(trifluoromethyl)quinolin-4-yl](hydroxy)methyl}piperidin-1-ium trichloroacetate (C19H17N2OF6+·C2Cl3O2−) (III), Fig. 4(a) (Wardell et al., 2018). In (III), C—X⋯π contacts involving X = F and Cl are evident, and these are readily apparent when the Hirshfeld surface is mapped with the shape-index property as shown in Fig. 4(b) and (c).
To examine the influence of π–π stacking on the molecular packing, an analysis of the Hirshfeld surface mapped over the shape-index and curvedness properties can be instructive. The structure of (I), Fig. 1(a), is used an example as the crystal features π(pyrrole)–π(nitrobenzene) stacking interactions (Zukerman-Schpector, Soto-Monsalve et al., 2018). Two views of the Hirshfeld surface mapped over the shape-index property are shown in Fig. 5(a) and (b). From these, the π–π stacking between the rings is indicated by the appearance of small blue regions surrounding bright-red spots within the respective five- and six-membered rings. The presence of π–π stacking is also evident as the flat regions around the pyrrole and benzene rings on the Hirshfeld surface mapped over curvedness, Fig. 5(c).
Two-dimensional fingerprint plots can also be used to analyse the calculated Hirshfeld surface of a molecule. Typically, the overall fingerprint plot is calculated, encompassing all intermolecular contacts, as well as the delineated (or decomposed) fingerprint plots, which focus on specific interactions. This is illustrated for the structure of S-benzyl 3-[1-(6-methylpyridin-2-yl)ethylidene]dithiocarbazate (C16H17N3S2) (IV), Fig. 6(a) (Omar et al., 2018). The overall fingerprint plot for (IV) is shown in Fig. 6(b) and those delineated into H⋯H, C⋯H/H⋯C, S⋯H/H⋯S and N⋯H/H⋯N interactions are shown in Fig. 6(c)–(f), respectively. While it is likely there are other identifiable points of contact that can be highlighted in the crystal, these may be of limited significance and do not require detailed discussion nor illustration. In the present case of (IV), the relative percentage contributions to the overall Hirshfeld surface are presented in Table 1. Ideally, in the absence of rounding-up errors, the relative percentage contributions should sum to 100%. Fingerprint plots would normally be presented for the more significant contributions to the surface unless a special feature of the molecular packing deserves highlighting. As seen from Fig. 6(b), the overall two-dimensional fingerprint plot is the sum of the delineated plots, having features drawn from the plots shown in Fig. 6(c)–(f). It is usually the case that the main contribution to the overall surface arises from H⋯H contacts. Also noteworthy is that while often forming the focus of discussion, conventional hydrogen bonding often makes relatively small percentage contributions to the overall surface.
|
When evaluating fingerprint plots, peaks/tips/features occurring at values less than the sum of the van der Waals radii need to be looked for. For example, in the case of H⋯H contacts, Fig. 6(c), the tip occurs at de + di < 2.40 Å, i.e. less than 2 × the van der Waals radius of hydrogen, suggestive of some of sort of contact, whether it be attractive or repulsive. The same procedure is followed for all other contacts. In (IV), the forceps-like tips in Fig. 6(d) and (e) correspond to interactions less than the sum of the respective van der Waals radii but, not so in Fig. 6(f).
Hirshfeld surface analyses are equally useful for assessing multi-component crystals, including solvates, salts and structures with Z′ > 1. In these situations, not only should the overall fingerprint plots be plotted but also those for the individual components. In a recent study where four cations and four anions comprised the crystallographic unit, distinctive features were evident in the fingerprint plots and in the relative percentage contributions of different interactions to the Hirshfeld surfaces for each individual component of the structure, which enabled the confirmation of the (Jotani et al., 2019). The calculation of Hirshfeld surfaces over the electrostatic potential will indicate interacting regions of the constituent molecules and can often be a useful starting point for analysis. Less confidence in interpretation will be likely in structures featuring disorder.
3. Non-covalent interaction plots
It is a fair assumption that under ambient conditions molecules, by and large, assemble into crystals optimizing attractive interactions while at the same time minimizing repulsive interactions. Given the nature and broad range of different intermolecular interactions now widely discussed in the crystallographic literature, it is salient to confirm whether such interactions are indeed attractive and therefore, stabilizing. In their landmark paper entitled `Revealing Noncovalent Interactions', Yang and co-workers (Johnson et al., 2010; Contreras-García et al., 2011) put forward a convenient, rapid and user-friendly approach to enable the discrimination between attractive and repulsive interactions. The method relies solely on the three-dimensional atomic coordinates and is equally applicable to macromolecular systems. The program NCIPLOT may be downloaded, again without charge, from https://www.lct.jussieu.fr/pagesperso/contrera/nciplot.html.
In short, the reduced density gradient is plotted as a function of the density (mapped as isosurfaces) over the molecule of interest. The sign of the second Hessian eigenvalue times the electron density [i.e. sign(λ2)ρ in atomic units] enables the identification of attractive/stabilizing (favourable) or repulsive (unfavourable) interactions. The derived results are readily visualized employing the VMD (visual molecular dynamics) molecular graphics viewer (Humphrey et al., 1996), which is freely available from https://www.ks.uiuc.edu/Research/vmd/.
The nature of the specific interactions is highlighted through a red–blue–green colour scheme on the calculated isosurface. A strong attractive interaction is indicated in blue whereas red indicates a strong repulsive interaction. Weak interactions are highlighted by a green isosurface.
A recent example employing this approach is illustrated for the structure of 4-(4-acetyl-5-methyl-1H-1,2,3-triazol-1-yl)benzonitrile (V), Fig. 7(a) (Zukerman-Schpector, Dias et al., 2018). One identified contact between centrosymmetrically related molecules is a carbonyl-C=O⋯π(triazolyl) interaction, where the interacting species are approximately parallel. As seen in the images of Fig. 7(b), a green isosurface is evident between the participating residues suggesting the interaction is weakly attractive. Fig. 7(c) shows an overall plot of the reduced density gradient versus the electron density times the sign of the second Hessian eigenvalue. Fig. 7(d) is an expanded version of Fig. 7(c) highlighting the weakly attractive nature of carbonyl-C=O⋯π(triazolyl) interaction in the negative region of the plot.
4. Interaction energies and energy frameworks
A new feature has been recently incorporated into Crystal Explorer 17 to enable the calculation of pair-wise interaction energies within a crystal by summing up four energy components comprising electrostatic (Eele), polarization (Epol), dispersion (Edis) and exchange-repulsion (Erep) (Turner et al., 2015). Users may apply two energy models available in the software to perform the calculation, i.e. CE-B3LYP/6-31G(d,p) and CE-HF/3-21G, which have been appropriately scaled to reproduce the B3LYP-D2/6-31G(d,p) counterpoise-corrected energies with a small mean absolute deviation of 2.4 and 4.7 kJ mol−1 for the respective models based on a set of crystal structures covering neutral organic molecules, organic salts, solvates, coordination compounds and radicals (Turner et al., 2015; Mackenzie et al., 2017).
The calculation of interaction energies is generally straightforward for crystal structures with Z′ = 1, whereby users simply need to generate a cluster of molecules within a radius of 3.8 Å (i.e. the default value for molecules comprising light atoms) for a selected reference molecule and subsequently subject it to energy calculation upon setting the relevant parameters such molecular charge, multiplicity and energy model. For multi-component crystals or crystals with Z′ > 1, the wave-functions for each unique molecule need to be calculated prior to obtaining the interaction energies for a cluster of molecules, as this is to ensure that the terrain of energy will encompass all pair-wise energies between the unique molecules, be they hetero-molecules (A⋯B) or within homo-molecules (A⋯A′ or B⋯B′), etc. An example of such a calculation was demonstrated in a recent study on the comprising two molecules of 2,2′-thiodibenzoic acid (S1 and S2) and four molecules of triphenylphosphane oxide (P1, P2, P3 and P4) in the (VI), Fig. 8(a) (Tan & Tiekink, 2018). Here, the pair-wise energy was first obtained for the respective pairs of interacting molecules (i.e. S1⋯P1, S1⋯P4, S2⋯P3 and S2⋯P4) prior to the calculation of interaction energies within 3.8 Å for the S1 and S2 clusters.
Useful information can be obtained upon the successful calculation of interaction energies. For instance, the calculation results in a colour-coded molecular cluster related to the specific interaction energy, Fig. 8(b). The individual energy components (Eele, Epol, Edis and Erep) as well as the sum of energy components (Etot) for the interactions relative to the reference molecule (based on the colour scheme) are provided in the accompanying table under the information dialogue; the individual energy components are not scaled but the Etot is scaled according to the relevant energy model (Mackenzie et al., 2017). Apart from these energy data, other information can be obtained from the generated table, such as the existence of rotational symmetry operations with respect to the reference molecule (Symop), the centroid-to-centroid distance between the reference molecule and interacting molecules (R), as well as the number of pair(s) of interacting molecules with respect to the reference molecule (N), which is useful in calculating a lattice energy of a crystal.
As mentioned in the Hirshfeld surface analysis and two-dimensional fingerprint plots section, Hirshfeld surface analysis is used to identify any close contacts present in a crystal through mapping of dnorm on the pro-molecule surface, and the strength of the close contacts may be estimated qualitatively through the intensity of the red spots observed on the surface or via the di + de contact distance as determined from a delineated fingerprint plot. With the availability of an immensely useful feature in the newly released Crystal Explorer 17, users may now quantify the strength of contacts by calculating the interaction energies and correlate this information with the results of the Hirshfeld surface analysis. This feature is especially useful in crystal engineering, for which it can be applied to compare and subsequently fine-tune the strength of interactions for any closely related analogues in designing structures with specific interactions for desirable applications. This idea is illustrated in a recent study of the 2:1 formed between 2,2′-dithiodibenzoic acid and benzoic acid (C14H10O2S2·C7H6O2), (VII), Fig. 9(a) (Tan & Tiekink, 2019). The interactions between the carboxylic acid residues via {⋯HOC=O}2 synthons in 2,2′-dithiodibenzoic acid (Humphrey & Wood, 2003), Fig. 9(b), are the same as in the structure of this conformer in the 1:2 with benzoic acid, and about 10 × greater than a benzene-C—H⋯O(hydroxyl) interaction, Fig. 9(c).
An option also exists in the new version of the Crystal Explorer 17 software to simulate energy frameworks, i.e. a graphical representation of the individual energy components depicted as cylinders joining the centroids of interacting molecular pairs, in which Eele, Edis and Etot are, respectively, colour-coded in red, green and blue, and with the radius of the corresponding cylinders proportional to the magnitude of interaction energy (Turner et al., 2015).
The simulation of the energy framework is an extended feature established based on the calculation of interaction energies. To simulate a framework, users first need to obtain the wave-functions for all unique pairs of interacting molecules as described earlier. Subsequently, a cluster of molecules within an appropriate number of unit cells needs to be generated depending on the completeness of the framework, e.g. a cluster of molecules within 2 × 2 × 2 unit cells may be a good start. Upon the completion of the energy calculations for the molecular cluster within the unit cells, the frameworks can be obtained through manifestation of the corresponding cylinder rods; these may need to be adjusted by a scale factor for direct comparison. An appropriate can be set to omit any weak interactions for purposes of clarity. An illustrative example is given in Fig. 10 for the structure of (VII) (Tan & Tiekink, 2019).
The calculation of energy frameworks was developed to better understand the topology of the overall interaction energies between the constituents of a crystal. For example, such an approach has found application in rationalization of the mechanical behaviour of drugs with relation to their tabletability (the ease of forming a tablet from a powder) (Turner et al., 2015). The importance of this functionality can be clearly seen when it is applied to polymorphs, as it allows users to directly compare the topological differences of the energy components between the structures, and potentially enable the correlation of energy frameworks with the physicochemical properties or packing behaviour of the polymorph of interest. As an example, the calculated energy frameworks for two conformational polymorphs of 4-(2H-1,3-benzodioxol-5-yl)-1-(4-methylphenyl)-1H-pyrazol-5-amine (C15H13N3O2) (VIII), Fig. 11(a) (Gajera et al., 2013; Jotani et al., 2015) is described. One polymorph is triclinic with Z′ = 2 (Gajera et al., 2013) while the other is monoclinic with Z′ = 1 (Jotani et al., 2015). The main difference between the two polymorphs is that one of the independent molecules in the triclinic form adopts a syn disposition for the dioxolyl fused-ring system with respect to the amino substituent connected to the central pyrazolyl ring but the other adopts an anti-arrangement (Gajera et al., 2013). In the monoclinic form, the molecules appear entirely in the syn form (Jotani et al., 2015). Through a powder X-ray diffraction study, it was found that the syn- and anti-orientations exist in 3:1 ratio (Jotani et al., 2015). This result is affirmed by a study of the energy frameworks for the polymorphs in which the monoclinic form exhibits a more compact framework in comparison to the triclinic form, as evidenced by the relatively thicker cylindrical radius at the same scale factor, which gives an indication that greater stabilization energies exist in the monoclinic system, Fig. 11(b)–(d).
5. Conclusion
The ready availability and ease of use of Crystal Explorer 17, including the calculation of energy frameworks, and NCIPLOT suggests these should be routinely employed tools in describing the molecular packing, as they complement the geometric analysis provided by the indispensable tool, PLATON. In short, utilizing these additional tools will ensure that the practitioner will get the most out of their experiments.
Acknowledgements
Sunway University Sdn Bhd is thanked for supporting the crystallographic laboratory in the Research Centre for Crystalline Materials at Sunway University.
References
Alcock, N. W. (1972). Adv. Inorg. Chem. Radiochem. 15, 1–58. CrossRef CAS Google Scholar
Allen, F. H. (2002). Acta Cryst. B58, 380–388. Web of Science CrossRef CAS IUCr Journals Google Scholar
Bauzá, A. & Frontera, A. (2015). Angew. Chem. Int. Ed. 54, 7340–7343. Google Scholar
Bauzá, A., Mooibroek, T. J. & Frontera, A. (2013). Angew. Chem. Int. Ed. 52, 12317–12321. Google Scholar
Bauzá, A., Mooibroek, T. J. & Frontera, A. (2015). ChemPhysChem, 16, 2496–2517. Web of Science PubMed Google Scholar
Boese, R., Kirchner, M. T., Dunitz, J., Filippini, G. & Gavezzotti, A. (2001). Helv. Chim. Acta, 84, 1561–1577. CrossRef CAS Google Scholar
Brinck, T., Murray, J. S. & Politzer, P. (1992). Int. J. Quantum Chem. 44, 57–64. CrossRef Google Scholar
Caracelli, I., Zukerman-Schpector, J., Haiduc, I. & Tiekink, E. R. T. (2016). CrystEngComm, 18, 6960–6978. CrossRef CAS Google Scholar
Caracelli, I., Zukerman-Schpector, J. & Tiekink, E. R. T. (2013). Gold Bull. 46, 81–89. Web of Science CrossRef CAS Google Scholar
Cavallo, G., Metrangolo, P., Milani, R., Pilati, T., Priimagi, A., Resnati, G. & Terraneo, G. (2016). Chem. Rev. 116, 2478–2601. Web of Science CrossRef CAS PubMed Google Scholar
Contreras-García, J., Johnson, E. R., Keinan, S., Chaudret, R., Piquemal, J.-P., Beratan, D. N. & Yang, W. (2011). J. Chem. Theory Comput. 7, 625–632. Web of Science PubMed Google Scholar
Dance, I. (2003). New J. Chem. 27, 22–27. Web of Science CrossRef CAS Google Scholar
Desiraju, G. R. (1995). Angew. Chem. Int. Ed. Engl. 34, 2311–2327. CrossRef CAS Web of Science Google Scholar
Dunitz, J. D. & Gavezzotti, A. (2009). Chem. Soc. Rev. 38, 2622–2633. Web of Science CrossRef PubMed CAS Google Scholar
Edwards, A. J., Mackenzie, C. F., Spackman, P. R., Jayatilaka, D. & Spackman, M. A. (2017). Faraday Discuss. 203, 93–112. CrossRef CAS Google Scholar
Frisch, M. J., Trucks, G. W., Schlegel, H. B., Scuseria, G. E., Robb, M. A., Cheeseman, J. R., Scalmani, G., Barone, V., Mennucci, B., Petersson, G. A., Nakatsuji, H., Caricato, M., Li, X., Hratchian, H. P., Izmaylov, A. F., Bloino, J., Zheng, G., Sonnenberg, J. L., Hada, M., Ehara, M., Toyota, K., Fukuda, R., Hasegawa, J., Ishida, M., Nakajima, T., Honda, Y., Kitao, O., Nakai, H., Vreven, T., Montgomery, J. A. Jr, Peralta, J. E., Ogliaro, F., Bearpark, M., Heyd, J. J., Brothers, E., Kudin, K. N., Staroverov, V. N., Kobayashi, R., Normand, J., Raghavachari, K., Rendell, A., Burant, J. C., Iyengar, S. S., Tomasi, J., Cossi, M., Rega, N., Millam, J. M., Klene, M., Knox, J. E., Cross, J. B., Bakken, V., Adamo, C., Jaramillo, J., Gomperts, R., Stratmann, R. E., Yazyev, O., Austin, A. J., Cammi, R., Pomelli, C., Ochterski, J. W., Martin, R. L., Morokuma, K., Zakrzewski, V. G., Voth, G. A., Salvador, P., Dannenberg, J. J., Dapprich, S., Daniels, A. D., Farkas, Ö., Foresman, J. B., Ortiz, J. V., Cioslowski, J. & Fox, D. J. (2016). Gaussian16. Gaussian, Inc., Wallingford, Connecticut, USA. Google Scholar
Gajera, N. N., Patel, M. C., Jotani, M. M. & Tiekink, E. R. T. (2013). Acta Cryst. E69, o736–o737. CrossRef IUCr Journals Google Scholar
Haiduc, I. (1997). Coord. Chem. Rev. 158, 325–358. CrossRef CAS Google Scholar
Humphrey, W., Dalke, A. & Schulten, K. (1996). J. Mol. Graph. 14, 33–38. Web of Science CrossRef CAS PubMed Google Scholar
Humphrey, S. M. & Wood, P. T. (2003). Acta Cryst. E59, o1364–o1366. CrossRef IUCr Journals Google Scholar
Janiak, C. (2000). J. Chem. Soc. Dalton Trans. pp. 3885–3896. Web of Science CrossRef Google Scholar
Jayatilaka, D. & Grimwood, D. J. (2003). Tonto: A Fortran Based Object-Oriented System for Quantum Chemistry and Crystallography. In: Computational Science – ICCS 2003. ICCS 2003. Lecture Notes in Computer Science, Vol. 2660, edited by P. M. A. Sloot, D. Abramson, A. V. Bogdanov, Y. E. Gorbachev, J. J. Dongarra & A. Y. Zomaya. Berlin, Heidelberg: Springer. Google Scholar
Johnson, E. R., Keinan, S., Mori-Sánchez, P., Contreras-García, J., Cohen, A. J. & Yang, W. (2010). J. Am. Chem. Soc. 132, 6498–6506. Web of Science CrossRef CAS PubMed Google Scholar
Jotani, M. M., Gajera, N. N., Patel, M. C., Sung, H. H. Y. & Tiekink, E. R. T. (2015). Acta Cryst. E71, 1121–1124. CrossRef IUCr Journals Google Scholar
Jotani, M. M., Wardell, J. L. & Tiekink, E. R. T. (2019). Z. Kristallogr. Cryst. Mat. 234, 43–57. CrossRef CAS Google Scholar
Koenderink, J. J. (1990). Solid Shape. Cambridge MA: MIT Press. Google Scholar
Koenderink, J. J. & van Doorn, A. J. (1992). Image Vis. Comput. 10, 557–564. CrossRef Google Scholar
Kolář, M. H. & Hobza, P. (2016). Chem. Rev. 116, 5155–5187. Web of Science PubMed Google Scholar
Mackenzie, C. F., Spackman, P. R., Jayatilaka, D. & Spackman, M. A. (2017). IUCrJ, 4, 575–587. Web of Science CrossRef CAS PubMed IUCr Journals Google Scholar
Malenov, D. P., Janjić, G. V., Medaković, V. B., Hall, M. B. & Zarić, S. D. (2017). Coord. Chem. Rev. 345, 318–341. Web of Science CrossRef CAS Google Scholar
Masui, H. (2001). Coord. Chem. Rev. 219–221, 957–992. Web of Science CrossRef CAS Google Scholar
McKinnon, J. J., Jayatilaka, D. & Spackman, M. A. (2007). Chem Commun. pp. 3814–3816. Web of Science CrossRef Google Scholar
McKinnon, J. J., Spackman, M. A. & Mitchell, A. S. (2004). Acta Cryst. B60, 627–668. Web of Science CrossRef CAS IUCr Journals Google Scholar
Murray, J. S., Lane, P., Clark, T. & Politzer, P. (2007). J. Mol. Model. 13, 1033–1038. CrossRef CAS Google Scholar
Nishio, M. (2004). CrystEngComm, 6, 130–158. Web of Science CrossRef CAS Google Scholar
Omar, S. A., Chah, C. K., Ravoof, T. B. S. A., Jotani, M. M. & Tiekink, E. R. T. (2018). Acta Cryst. E74, 261–266. CrossRef IUCr Journals Google Scholar
Scheiner, S. (2013). Acc. Chem. Res. 46, 280–288. Web of Science CrossRef CAS PubMed Google Scholar
Schmidbaur, H. (2000). Gold Bull. 33, 3–10. CrossRef CAS Google Scholar
Schmidbaur, H. & Schier, A. (2008). Chem. Soc. Rev. 37, 1931–1951. Web of Science CrossRef PubMed CAS Google Scholar
Spackman, M. A. & Jayatilaka, D. (2009). CrystEngComm, 11, 19–32. Web of Science CrossRef CAS Google Scholar
Spackman, M. A. & McKinnon, J. J. (2002). CrystEngComm, 4, 378–392. Web of Science CrossRef CAS Google Scholar
Spek, A. L. (2009). Acta Cryst. D65, 148–155. Web of Science CrossRef CAS IUCr Journals Google Scholar
Sredojević, D., Bogdanović, G. A., Tomić, Z. D. & Zarić, S. D. (2007). CrystEngComm, 9, 793–798. Google Scholar
Tan, S. L. & Tiekink, E. R. T. (2018). Acta Cryst. E74, 1764–1771. CrossRef IUCr Journals Google Scholar
Tan, S. L. & Tiekink, E. R. T. (2019). Acta Cryst. E75, 1–7. CrossRef IUCr Journals Google Scholar
Tiekink, E. R. T. (2014). Coord. Chem. Rev. 275, 130–153. Web of Science CrossRef CAS Google Scholar
Tiekink, E. R. T. (2017). Coord. Chem. Rev. 345, 209–228. Web of Science CrossRef CAS Google Scholar
Tiekink, E. R. T. & Zukerman-Schpector, J. (2011). Chem. Commun. 47, 6623–6625. Web of Science CrossRef CAS Google Scholar
Turner, M. J., Mckinnon, J. J., Wolff, S. K., Grimwood, D. J., Spackman, P. R., Jayatilaka, D. & Spackman, M. A. (2017). Crystal Explorer 17. The University of Western Australia. Google Scholar
Turner, M. J., Thomas, S. P., Shi, M. W., Jayatilaka, D. & Spackman, M. A. (2015). Chem. Commun. 51, 3735–3738. Web of Science CrossRef CAS Google Scholar
Wang, W., Ji, B. & Zhang, Y. (2009). J. Phys. Chem. A, 113, 8132–8135. CrossRef Google Scholar
Wardell, J. L., Jotani, M. M. & Tiekink, E. R. T. (2018). Acta Cryst. E74, 1851–1856. CrossRef IUCr Journals Google Scholar
Yusof, E. N. Md., Tiekink, E. R. T., Jotani, M. M., Simone, M. I., Page, A. J. & Ravoof, T. B. S. A. (2018). J. Mol. Struct. 1171, 650–657. Google Scholar
Zukerman-Schpector, J., Soto-Monsalve, M., De Almeida Santos, R. H., Garcia, A. L. L., Correia, C. R. D., Jotani, M. M. & Tiekink, E. R. T. (2018). Acta Cryst. E74, 371–375. CrossRef IUCr Journals Google Scholar
Zukerman-Schpector, J., Dias, C. da S., Schwab, R. S., Jotani, M. M. & Tiekink, E. R. T. (2018). Acta Cryst. E74, 1195–1200. 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.