|Negative thermal expansion and metallophilicity in Cu[Co(CN)]|
|Adam F. Sapnik, Xiaofei Liu, Hanna L. B. Boström, Chloe S. Coates, Alistair R. Overy, Emily M. Reynolds, Alexandre Tkatchenko, and Andrew L. Goodwin|
|We report the synthesis and structural characterisation of the molecular framework copper(I) hexacyanocobaltate(III), Cu[Co(CN)], which we find to be isostructural to H[Co(CN)] and the colossal negative thermal expansion material Ag[Co(CN)]. Using synchrotron X-ray powder diffraction measurements, we find strong positive and negative thermal expansion behaviour respectively perpendicular and parallel to the trigonal crystal axis: = 25.4(5) MK and = 43.5(8) MK. These opposing effects collectively result in a volume expansivity = 7.4(11) MK that is remarkably small for an anisotropic molecular framework. This thermal response is discussed in the context of the behaviour of the analogous H- and Ag-containing systems. We make use of density-functional theory with many-body dispersion interactions (DFT+MBD) to demonstrate that CuCu metallophilic (‘cuprophilic’) interactions are significantly weaker in Cu[Co(CN)] than AgAg interactions in Ag[Co(CN)], but that this lowering of energy scale counterintuitively translates to a more moderate—rather than enhanced—degree of structural flexibility. The same conclusion is drawn from consideration of a simple lattice dynamical model, which we also present here. Our results demonstrate that strong interactions can actually be exploited in the design of ultra-responsive materials if those interactions are set up to act in tension.|
State Key Laboratory of Mechanics and Control of Mechanical Structures, Key Laboratory for Intelligent Nano Materials and Devices of Ministry of Education, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
Diamond Light Source, Chilton, Oxfordshire OX11 0DE, U.K.
Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg††footnotetext: † Electronic Supplementary Information (ESI) available: X-ray powder diffraction measurements and Rietveld fits, lattice parameter data and structural parameter data, ab initio van der Waals parameters and energies, GULP input files. See DOI: 10.1039/b000000x/
The development of responsive materials often exploits weak interactions as key design elements because lower interaction energies heighten the sensitivity of a material to external perturbations.1, 2, 3, 4 It is no accident, for example, that the weak inter-molecular forces in molecular crystals generally allow more extreme responses to changes in temperature5, 6 and pressure7, 8 than is possible in conventional inorganic ceramics, the structures of which are held together by strong ionic and covalent bonding networks. In this context, supramolecular interactions assume a particular importance, given that their energy scales are so much lower than those of electrostatic or covalent interactions. Hence the prevalence of hydrogen-bonding,9 halogen-bonding,10 –,11 van der Waals (vdW),5 host–guest,12, 13 and metallophilic14 interactions amongst many of the important materials in the field.
Thermal expansion behaviour is a straightforward measure of responsiveness: it quantifies the effect of temperature on the linear dimensions of a material.15 Compounds with large thermal expansion coefficients often show extreme and counterintuitive responses to pressure,16, 17 and may harbour various other anomalous elastic properties, such as negative Poisson’s ratios18 or thermosalient19, 20 effects. So it is perhaps unsurprising that some of the most extreme (‘colossal’) thermal expansion known has been observed in framework materials whose lattice dimensions are a function of weak metallophilic interactions.14, 17, 21 The canonical system of this type is Ag[Co(CN)], which adopts a lattice structure22 that can flex in such a way as to vary argentophilic AgAg separations without affecting covalent interactions within the lattice itself.23, 24 A geometric consequence of this flexing behaviour is that the positive thermal expansion (PTE) of argentophilic interactions (i.e. increase in separation with increasing temperature) is translated into a negative thermal expansion (NTE) effect in a perpendicular direction [Fig. 1]. The same mechanism operates under application of hydrostatic pressure, such that volume compression actually results in linear expansion for a particular set of directions17—so-called negative linear compressibility (NLC).25, 26, 27 NTE and NLC are valuable material properties, exploitable in the design of athermal composites used in optical devices and next-generation pressure sensors.
In seeking to design even more responsive analogues of Ag[Co(CN)], we considered the possibility of replacing Ag by Cu. Metallophilic interactions involving Cu ions are perhaps less well studied than argentophilic and aurophilic interactions, but are expected to be weaker given the reduced polarisability of the shell.28, 29 Hence, by the arguments discussed above, Cu[Co(CN)] has always been an obvious candidate for extreme thermomechanical response. To the best of our knowledge, this system has never previously been reported: the difficulty of preparing the phase is likely a consequence of the propensity for Cu to disproportionate under the aqueous reaction conditions used to prepare the family of materials related to Ag[Co(CN)].23 We have recently exploited the Cu reduction protocol developed in Ref. 30 to allow access to otherwise unrealisable Cu(I)-containing frameworks,31 suggesting that a similar synthetic approach may provide an alternative synthetic entry point to Cu[Co(CN)].
Here we validate such an approach, reporting the synthesis, crystal structure, and thermal expansion behaviour of Cu[Co(CN)]. Using a combination of high-resolution synchrotron X-ray diffraction measurements and Rietveld refinement, we show the system to be isostructural with Ag[Co(CN)] and H[Co(CN)].22, 32, 33, 34 Variable-temperature (100–598 K) diffraction measurements allow determination of the corresponding coefficients of thermal expansion , which we find to be substantially less extreme than those of Ag[Co(CN)] (even if they remain large in the context of the behaviour conventional inorganic solids35). In particular, our data give MK and MK; cf MK and MK for Ag[Co(CN)].23 In order to rationalise this more moderate thermomechanical response in terms of the relative strengths of CuCu and AgAg metallophilic interactions, we carry out a series of ab initio calculations. Our analysis suggests (i) that cuprophilic interactions are indeed weaker than argentophilic interactions in this family, and (ii) the more extreme thermomechanical response of the Ag-containing compound is a result of the balance of metallophilic and electrostatic interaction energies rather than a signature of particularly weak argentophilicity. Lattice dynamical calculations using a highly simplified interaction model relevant to the entire A[Co(CN)] structural family lead to the same conclusions. Our results suggest that competing interactions—rather than low-energy interactions per se—might be key in the design of ultra-responsive materials.
All reagents were obtained from commercial suppliers and used as received.
2.1 Copper(I) hexacyanocobaltate(III)
We prepared polycrystalline samples of copper(I) hexacyanocobaltate(III) following a modification of the reduction protocol reported in Refs. 30, 31. A saturated aqueous solution of copper(II) sulfate (Sigma Aldrich, 99%; 0.17705 g) was added dropwise to a concentrated aqueous solution of sodium bisulfite (Sigma Aldrich, 0.05771 g), present in stoichiometric excess, to afford a mint-green solution. The solution was stirred for 30 min, after which time an aqueous solution of potassium hexacyanocobaltate(III) (Sigma Aldrich, 97%, 0.12288 g; stoichiometric with respect to copper) was added dropwise to afford a pale blue precipitate. The solution was stirred for a further 2 h, and the pale-blue solid product formed was isolated by filtration, washing (HO) and drying under vacuum. The solid contained a mixture of copper(I) hexacyanocobaltate(III) and Prussian-blue-structured potassium copper(II) hexacyanocobaltate(III), a seemingly inescapable by-product of this synthetic strategy.
Copper(I) hexacyanocobaltate(III) could also be obtained in impure form using mechanochemical synthesis. Stoichiometric quantities of solid tetrakis(acetonitrilo)copper(I) hexafluorophosphate (Chem Cruz, 98%, 0.41346) and potassium hexacyanocobaltate (Sigma Aldrich, 97%, 0.12288g) were combined in an agate mortar, and intimately mixed via solid-state grinding for 30 min. An obvious colour change from white to pale blue occurred during this process. The resulting solid was washed (HO) and dried to afford a mixture of copper(I) hexacyanocobaltate(III), potassium copper(II) hexacyanocobaltate(III) and at least one further unidentified product. Given the reduced purity of this product, the solution-phase product described above was used for all diffraction measurements carried out in this study.
2.2 Powder X-ray diffraction
High-resolution synchrotron powder diffraction measurements were carried out using the I11 beamline at the Diamond Light Source. A finely-ground sample of copper(I) hexacyanocobaltate(III), prepared as above, was loaded into a borosilicate capillary (0.5 mm diameter) and mounted on the diffractometer. Diffraction patterns were recorded using the Mythen2 point sensitive detector over the angular range –, using an X-ray wavelength Å calibrated by refinement of a silicon NIST 640c standard. Each measurement consisted of two scans of 5 s exposure, offset relative to one another by . The sample temperature was controlled using an Oxford Cryostream (100–500 K) and a Cyberstar hot air blower (523–598 K). Diffraction patterns were measured at intervals of 25 K between 100 and 500 K and again between 523 and 598 K.
Both Pawley and Rietveld refinements were carried out using TOPAS Academic (version 4.1).36 We employed a modified Thompson–Cox–Hasting pseudo-Voigt (TCHZ) peak shape, combined with a simple axial divergence correction and a Stephens anisotropic peak broadening term.37 The potassium copper(II) hexacyancobaltate(III) impurity phase was modelled using Pawley refinement of the double-metal cyanide cell (10 Å).38 Rietveld refinement of the Cu[Co(CN)] phase made use of a starting model based on the known structure of Ag[Co(CN)].22 Refinement was stable for all temperature points, provided that Co–C/C–N bond distance restraints and a single isotropic displacement parameter for all atom types were used in the Rietveld model. Sequential (seed-batch) Rietveld refinements, where the starting structural parameters for each temperature point were those used at the preceding temperature, provided structural models with physically-sensible temperature dependencies for K. For the temperature regime K, we found that the positional coordinates of the C and N atoms and the value of showed strong covariance, and hence we have reduced confidence in the absolute values of these parameters. This regime corresponds to the temperature range over which decomposition of the KCu[Co(CN)] phase appears to set in.
2.3 Thermal expansivity determination
Thermal expansivities were calculated using the PASCal software.39 We employed estimated temperature uncertainties of 5 K and fitted the principal axis expansivities using linear functions. For internal consistency with the uniaxial expansivities, the volume expansivity was determined using the trace of the expansivity tensor40 rather than via the direct – fit given by PASCal.39
2.4 Ab initio calculations
Ab initio calculations were performed within the FHI-aims code,41 using the numeric atom-centred orbital tier 1 basis set for the wavefunction and a -point mesh for the Brillouin zone sampling. The Perdew-Burke-Ernzernhof (PBE) functional42 was used to model the semilocal exchange-correlation energy. To describe the non-local dispersion energies, we used both the interatomic pairwise Tkatchenko-Scheffler (TS) method,43 as well as the many-body dispersion (MBD) method, which includes many-body dipolar interatomic interactions to all orders in perturbation theory.44, 45 The lattice constants were obtained from unit cell relaxations with cell angles fixed to experimental values. Full a posteriori relaxation of the unit cell proved the reliability of this scheme.
2.5 Lattice dynamical calculations
Lattice dynamical calculations made use of the GULP software (version 4.4).46 Cell optimisations were carried out under constant pressure conditions and at , with strains constrained by symmetry. Dispersion interactions were modelled using a Buckingham potential with vanishingly small repulsive term, and the ‘c6’ flag was activated to employ Ewald-type summation. For all calculations, checks were carried out to ensure optimisation convergence and to verify the conservation of angle terms in the parameterisation.
3 Results and discussion
3.1 Crystal structure of Cu[Co(CN)]
The ambient-temperature X-ray powder diffraction pattern of our Cu[Co(CN)] sample is shown in Fig. 2, where it is compared to that of Ag[Co(CN)] as reported in Ref. 23. The structural similarity between the two phases is immediately evident, as is the presence of a substantial quantity of an impurity phase. We could account for the entire diffraction pattern using two components, one based on the Ag[Co(CN)] structure-type (space group symmetry ) and one with the cubic Prussian blue structure (space group symmetry ). This second phase would be consistent with the formation of KCu[Co(CN)] during synthesis, which is certainly feasible on chemical grounds.47, 48 A Pawley fit using this two-phase model confirms our assignment of space group symmetries and rules out the presence of any additional crystalline phases [Fig. 2(d)]. We note that there is good (if fortuitous) distinction between the diffraction profiles of the two phases present, which allows us to clearly distinguish the corresponding lattice parameters and their thermal expansion behaviour (see SI).
Having established the space group symmetry of Cu[Co(CN)] we proceeded to carry out a Rietveld refinement, employing a starting model based on the lattice parameters obtained during Pawley fitting and the published atom coordinates of Ag[Co(CN)].22 We continued to model the KCu[Co(CN)] phase using a Pawley fit—indeed this is the case for all subsequent refinements and is not discussed further. We found good stability in the refinement of our structural model of Cu[Co(CN)], obtaining a -value of ; the corresponding fit is shown in Fig. 2(e) and the relevant structural details are summarised in Table 1. A representation of the crystal structure itself is given in Fig. 3. All refined bond lengths are chemically sensible: we find a Co–C distance of 1.832(11) Å, which is similar to that in Ag[Co(CN)] (Co–C Å);22 likewise the Cu–N separation of 1.887(10) Å is comparable to that found in CuCN (Cu–C/N– Å).49
|300 K (experimental)||0 K (estimated)|
A property of the particular space group symmetry of Cu[Co(CN)] is that the CuCu separation is directly related to the lattice parameters:
Hence we find Å, which lies at the very upper bound of CuCu separations for which cuprophilic interactions are considered relevant.50 One crude measure of the strength of metallophilic interactions is to consider the ratio of the observed interatomic distance to the sum of the corresponding vdW radii.28 Using our room-temperature lattice parameters and the vdW radii given in Ref. 51 we obtain a ratio of 1.00 for Cu[Co(CN)], which is remarkably similar to the corresponding value for Ag[Co(CN)] (0.99).23 So, at face value, one might expect comparable metallophilic interaction strengths for the two systems.
3.2 Thermal expansion behaviour
Having collected a series of X-ray diffraction patterns over the temperature range 100–598 K, we carried out Rietveld refinements for each data set using the same approach described above. We obtained satisfactory refinements for all temperatures, albeit with some signs of increased uncertainties at the very highest temperatures—i.e. close to the onset of decomposition of the Prussian blue phase. The temperature dependence of the lattice parameters, illustrated in Fig. 4(a), was observed to be approximately linear over this entire temperature range. As in nearly all members of this structural family, Cu[Co(CN)] exhibits an NTE effect parallel to the crystal axis, and PTE effects in perpendicular directions (i.e., including the and crystal axes). Hence the basic thermomechanical response of this system can be understood in terms of the same ‘wine-rack’ mechanism illustrated in Fig. 1. The remaining structural parameters also show linear temperature dependencies [Fig. 4(b–d)]; taken together these values allow us to estimate a set of 0 K structural parameters that may prove useful for comparison against e.g. ab initio studies [Table 1].
Coefficients of thermal expansion were determined using linear fits to the lattice parameter data,39 and are given in Table 2. What is immediately apparent is that the magnitudes of both PTE and NTE effects in Cu[Co(CN)] are substantially smaller than those in the Ag-containing system. Consequently, Cu[Co(CN)] is not a colossal thermal expansion material, and its thermomechanical response shares more in common with other Cu-containing networks such as -Cu[C(CN)] (Ref. 31) and CuCN (Ref. 52) than with Ag[Co(CN)] and Ag[Fe(CN)].14 We will come to rationalise the differences in behaviour of the copper(I) and silver(I) hexacyanocobaltates below, but include first some additional analysis of the trends in lattice parameters we observe using our newly-measured data.
The ‘wine-rack’ mechanism that is thermally activated in this system can be interrogated directly using the so-called mechanical building unit (XBU) approach.2 We make use of the pair of transformations
which relate the unit cell dimensions to the framework strut length and framework angle [Fig. 3]. Using these same relationships, we can recast the lattice expansivities in terms of XBU expansivities, obtaining MK and MK. Hence the bulk of the thermal expansion response can be accounted for by changes in the framework geometry (); the lattice expansivities attributable to this flexing mechanism alone are MK and MK, where we use the prime notation to indicate calculation from . The observation indicates that the Co–CN–Cu–NC–Co ‘struts’ from which the framework structure of Cu[Co(CN)] is assembled actually contract with increasing temperature. This behaviour is likely due to thermal activation of transverse vibrational modes where lateral displacements of the chain (maximal at the Cu site) require shortening of the CoCo vector.53, 54 Such a mechanism is implicated in the uniaxial NTE of CuCN itself ( MK, Ref. 52, 49), and is presumably tempered here somewhat relative to that system by the increased strength of Co–C vs Cu–C bonds.47
One consequence of the negative value of is that the volume coefficient of thermal expansion of Cu[Co(CN)] is unusually small for systems in this particular family. Formally, this situation arises because of the fortuitous equivalence , which is the geometric requirement for .***Note that , and hence . Hence this material has the unusual property of (approximately) temperature-independent density despite its relatively large linear thermal expansivities. At face value, this property may be expected to result in unusually extreme uniaxial compressibilities under application of hydrostatic pressure, since small changes in volume would appear to be linked to large changes in lattice dimensions. However, we anticipate by analogy to related systems that the XBU compressibility is actually positive rather than negative, and so a small need not require a large bulk modulus.55, 2 Nevertheless we expect the particular uniaxial compressibility corresponding to the crystal axis to be negative, and so investigation of the NLC behaviour of Cu[Co(CN)] could prove a fruitful avenue of future research.
3.3 Ab initio calculations
The observation of more moderate thermal expansion behaviour in Cu[Co(CN)] relative to that in Ag[Co(CN)] poses a simple question: does this situation arise because cuprophilic interactions are actually stronger than argentophilic interactions, and hence less susceptible to changes in temperature?
In order to answer this question, we turn to ab initio calculations, which if carried out so as to include consideration of vdW interactions allow direct quantification of the metallophilic interactions in both compounds. We begin by reporting the 0 K structure for Cu[Co(CN)] obtained computationally and demonstrate that the inclusion of dispersive interactions is necessary to obtain good consistency with our experimental results. By mapping out the potential energy surface (PES) for all three A[Co(CN)] systems (A = H, Cu, Ag) across a variety of lattice strains and then taking into account the variation in vdW energies at each point, we extract the free-atom and in-solid (effective) coefficients. The value of these coefficients for each atom type A acts as a measure of the strength of metallophilic interactions in the corresponding A[Co(CN)] system.
|exp. (0 K)||PBE||TS||MBD|
The unit cell dimensions obtained in our DFT + vdW calculations are given in Table 3. The influence of dispersion energy on the lattice constants is large, just as is now known to be the case for Ag[Co(CN)].56 Our PBE calculation without vdW interactions overestimates and underestimates . Upon including dispersion interactions the lattice constants move closer to the experimental values. We note that the enhanced cohesive MBD energy for Cu[Co(CN)] arises from the collective effect of vdW interactions and the self-consistent polarisation in the unit cell.56 In Figure 5 we show a representative section of the PES for the three calculation regimes, and Figure 6 shows the TS and MBD vdW energies as a function of the individual and lattice constants. Our results make clear that the vdW energy depends more strongly on than it does on . Since the framework strut length is more rigid than the framework angle , then to lower the total energy the lattice simply contracts along (and ) while expanding along . Hence the same mechanism explains the qualitative change in lattice constants observed both as a result of using different vdW calculation methods and as a result of an increase in the polarisability of atom A. Indeed because the MBD energy depends almost linearly on the lattice constants it behaves as an effective pressure on the lattice, equivalent to 1.22 GPa along and 1.76 GPa along .
To compare the strength of cuprophilic interactions in Cu[Co(CN)] with that of argentophilic interactions in Ag[Co(CN)] we further analysed our DFT+vdW results. Our basic approach was to parameterise the vdW contribution to the TS-vdW energy in terms of dispersion coefficients and vdW radii for each atom type. In the PBE+TS calculations, the free-atom coefficient and vdW radii are used as the initial input parameters. The effect of the local chemical environment is taken into account by calculating the effective in-solid and as described in Ref. 43. Table 4 lists our results for the free-atom vdW parameters and the effective parameters for A[Co(CN)] (A = Ag, Cu, H) at the experimental lattice constants. We find that the argentophilic interactions are indeed stronger than cuprophilic interactions in these systems, as both the free-atom and effective values are larger by % for Ag relative to Cu. For completeness we note that the effect of the local chemical environment on the coefficients is to reduce the dispersion coefficients.
3.4 Lattice dynamical calculations
We supplement these high-level ab initio results with a series of extremely simple lattice-dynamical calculations that also allow us to estimate the relative strengths of metallophilic interactions in Cu[Co(CN)] and Ag[Co(CN)]. The approach we use is to develop the very simplest abstraction of all three A[Co(CN)] systems (A = H, Cu, Ag) that captures the key interactions responsible for their thermomechanical response. We parameterise this model with sufficiently few variables that six experimental observables (the two independent lattice parameters for each of the three systems) can be used to estimate the relative magnitudes of the metallophilic interactions in the A = Cu, Ag compounds.
The same structural model is used for all three systems: crystal symmetry, with a single anion (mass ) of charge at position and a cation () with charge at position [Fig. 7(a)]. These charges reflect the approximate Mulliken charges determined for H[Co(CN)] and Ag[Co(CN)] in Ref. 57 and are consistent with the Hirshfeld and Bader charges obtained in our own ab initio calculations (see SI). We refer to the anion using the symbol X (formally this corresponds to the [Co(CN)] ion), giving the unit cell composition AX. This structural model is then decorated with three interaction potentials: first, a harmonic bond potential between neighbouring A and X sites
second, a harmonic bond angle potential governing A–X–A triplets
and, third, (in the case of Cu and Ag systems) dispersive interactions between neighbouring A sites intended to reflect the empirical -dependence of metallophilic interactions58
In this way the lattice energy is given by the sum
In order to reduce the number of parameters involved in this model, we make the following assumptions. First, we take the effective charges at X and A sites to be system-independent. We justify this assumption by noting that the Mulliken charges reported for H[Co(CN)] and Ag[Co(CN)] vary more greatly by calculation method than they do between systems;57 the A = Cu case is intermediate to the A = H and A = Ag cases (see SI). Second, we take the flexing stiffness and equilbrium angle also to be system-independent. This is probably reasonable given that both terms will be governed by the chemistry of the [Co(CN)] anion, which is common to all three systems. Third, we take the (system-dependent) values of as the sum of bond lengths Co–CC–NN–A determined crystallographically: we use the values from Ref. 34 for A = H, from Ref. 17 for A = Ag, and from our present study for A = Cu. Fourth, we assume that the term is negligible for the A = H case, which is likely given the small electron density expected at the A site for this system (and is certainly supported by our own ab initio results as given in Table 4).
Our starting point is to determine a set of parameters that, when used to drive a lattice-dynamical geometry optimisation, result in the closest possible agreement between 0 K (derived from experiment) and relaxed cell parameters for A = H. We make the additional requirement that should be as close to as possible. The parameter set we obtain is listed in Table 5, together with a comparison of the experimental and simulated lattice parameters; the corresponding match in framework geometry is illustrated in Fig. 7(b).†††We found the quality of fit was relatively insensitive to changes in of up to ca 25% of its value. Variations in this parameter did affect the absolute values of the compressibilities determined subsequently; however the same trend in magnitudes of compressibilities shown in Fig. 8 was found in all cases. We note that we do not attach any particular physical meaning to the parameter values in our model, since (in particular) the charge distribution we use is a heavily simplified representation of reality. Nevertheless it is reassuring that even this simple model allows robust geometry optimisation to a physically-sensible state.
Having used the geometry of the A = H system to determine all of the system-independent parameter values, we proceeded to optimise the geometry of analogous models for A = Cu and Ag. In each case the value of was updated according to the experimental bond lengths, and only the value of was varied in order to obtain the closest match between calculated and experimental (0 K extrapolated) lattice parameters. The corresponding parameter values and optimised cell dimensions are again summarised in Table 5; we note that the level of agreement () is encouraging given the simplicity of the lattice dynamical model we have used. Also encouraging is that, for both compounds, the lattice parameters are overestimated in the absence of a metallophilic contribution to the lattice enthalpy. This mirrors the results of vastly higher-level ab initio geometry optimisations,57, 59 and indicates that the electrostatic contribution to the free energy (the single component of our model acting to increase ) operates in tension with the metallophilic interactions. While we do not attach any importance to the absolute values of the parameters that emerge from our calculations, what we do think is meaningful is the observation that is larger for A = Ag than for A = Cu. In other words, the experimental unit cell dimensions for Cu[Co(CN)] and Ag[Co(CN)] are consistent with stronger argentophilic interactions in the latter than cuprophilic interactions in the former. Moreover, the ratio of cuprophilic:argentophilic interaction strengths we deduce from our simple lattice-dynamical model is essentially the same as that obtained in our ab initio calculations: (Cu)/(Ag) = 61% vs 70%, respectively.
3.5 Flexibility from competing interactions
So our various calculations converge on the same scenario whereby cuprophilic interactions in Cu[Co(CN)] are weaker than argentophilic interactions in Ag[Co(CN)] by 30–40%. One obvious question remains: how is this observation consistent with the more moderate thermal expansion behaviour of the Cu-containing compound?
Here is the isothermal specific heat, the molar volume, the mean effective Grüneisen parameter and the uniaxial compressibilities. We estimate that the pre-factor varies by not more than between the A = Cu and A = Ag systems,‡‡‡Here we have made use of three relationships: first, that appears to be relatively system-independent;57 second, that the ratio of the values for A = Cu and Ag will be approximately equal to the ratio of the terms, since the low-energy phonon dispersion will be dominated by heavy-atom displacements; and third, we use the experimental molar volumes. such that a comparison of compressibilities for the two compounds provides a reasonable first-order approximation to the relative thermal expansivities. We concern ourselves with compressibilities rather than expansivities since the former are obtainable directly from the calculations (both ab initio and lattice-dynamical) described above. The relative compressibilities for all three compounds are illustrated graphically in Fig. 8. What is evident is that the Cu-containing compound exhibits intermediate behaviour to the H- and Ag-containing systems, despite its relatively weaker metallophilic interactions. The qualitative similarity to the relative thermal expansivities is striking, particularly given the (necessary) omission of anharmonic contributions from our calculations which likely contribute substantially to the experimental behaviour.61
4 Concluding remarks
We are led to the counterintuitive conclusion that stronger interactions can actually make a material more compliant: Ag[Co(CN)] exhibits colossal thermomechanical responses but Cu[Co(CN)] does not, despite the energy scale associated with metallophilic interactions being larger in the former than in the latter. Of course the key here is that metallophilic interactions are net attractive, and act in tension with the (repulsive) electrostatic component.59, 62 Any effective harmonic potential can be made increasingly shallow by the addition of attractive terms, as illustrated in Fig. 9. This is the nub of the physics at play in this family: in the absence of metallophilic interactions, the frameworks are not especially mechanically responsive but they do become so as metallophilicity is introduced.
Hence the conventional materials design rules are reversed, and we anticipate that the member of the A[Co(CN)] family likely to show the most extreme thermomechanical response is actually the as-yet-unrealised compound Au[Co(CN)]. It was shown in Ref. 57 that this system is likely to have a particularly compliant structure, although the degree of compliance will depend heavily on the strength of the aurophilic interaction contribution to the lattice enthalpy. Given the notorious difficulty of accessing aqueous Au(I) chemistry, it is not yet clear how Au[Co(CN)] might be accessed synthetically. A viable alternative is the (also unrealised) compound Fe[Au(CN)]—i.e., with Co(III) replaced by Fe(III) and the CN ion orientations reversed—which by analogy to Fe[Ag(CN)] should in principle be accessible via reaction of aqueous Fe-containing solutions with KAu(CN).63 The observation14 of qualitatively similar ‘colossal’ thermal expansion in Ag[Co(CN)] and Ag[Fe(CN)] suggests that chemical substitution at the trivalent metal site is unlikely to influence the degree of thermomechanical response observed.
From a computational perspective, one key implication of our study is the importance of obtaining accurate descriptions of vdW interactions in compliant framework materials. This importance is particularly acute for systems such as Cu[Co(CN)] and Ag[Co(CN)] where the PES is anomalously shallow as a result of competition between vdW and electrostatic contributions. A key challenge in this regard is the treatment of finite-temperature effects; i.e. anharmonicity. We anticipate that the discovery of anomalous mechanics in increasingly many systems based on vdW-type interactions64, 65 will motivate further research effort along precisely these lines.
The synchrotron diffraction measurements were carried out at the Diamond Light Source (I11 Beamline). We are extremely grateful for the award of a Block Allocation Grant, which made this work possible, and for the assistance in data collection provided by S. J. Cassidy (Oxford) and the I11 beamline staff. A.F.S., A.R.O. and A.L.G. gratefully acknowledge funding through the European Research Council (Grant 279705). C.S.C. and A.L.G. acknowledge funding through the Leverhulme Trust (Grant RPG-2015-292). A.R.O. thanks the Diamond Light Source for a Studentship. This project received funding from the European Union (EU) Horizon 2020 Research and Innovation Programme under Marie Sklodowska-Curie Grant Agreement 641887 (project acronym DEFNET).
- Coudert 2015 F.-X. Coudert, Chem. Mater., 2015, 27, 1905–1916.
- Ogborn et al. 2012 J. M. Ogborn, I. E. Collings, S. A. Moggach, A. L. Thompson and A. L. Goodwin, Chem. Sci., 2012, 3, 3011–3017.
- Horike et al. 2009 S. Horike, S. Shimomura and S. Kitagawa, Nature Chem., 2009, 1, 695–704.
- Krishna et al. 2016 G. R. Krishna, R. Devarapalli, G. Lal and C. M. Reddy, J. Am. Chem. Soc., 2016, 138, 13561–13567.
- Das et al. 2010 D. Das, T. Jacobs and L. J. Barbour, Nature Mater., 2010, 9, 36–39.
- Goodwin 2010 A. L. Goodwin, Nature Mater., 2010, 9, 7–8.
- Shepherd et al. 2012 H. J. Shepherd, T. Palamarciuc, P. Rosa, P. Guionneau, G. Molnár, J.-F. Létard and A. Bousseksou, Angew. Chem. Int. Ed., 2012, 51, 3910–3914.
- Duyker et al. 2016 S. G. Duyker, V. K. Peterson, G. J. Kearley, A. J. Studer and C. J. Kepert, Nature Chem., 2016, 8, 270–275.
- Li et al. 2014 W. Li, A. Thirumurugan, P. T. Barton, Z. Lin, S. Henke, H. H.-M. Yeung, M. T. Wharmby, E. G. Bithell, C. J. Howard and A. K. Cheetham, J. Am. Chem. Soc., 2014, 136, 7801–7804.
- Jones et al. 2014 R. H. Jones, K. S. Knight, W. G. Marshall, J. Clews, R. J. Darton, D. Pyatt, S. J. Coles and P. N. Horton, CrystEngComm, 2014, 16, 237–243.
- Yot et al. 2012 P. G. Yot, Q. Ma, J. Haines, Q. Yang, A. Ghoufi, T. Devic, C. Serre, V. Dmitriev, G. Férey, C. Zhong and G. Maurin, Chem. Sci., 2012, 3, 1100–1104.
- Engel et al. 2014 E. R. Engel, V. J. Smith, C. X. Bezuidenhout and L. J. Barbour, Chem. Commun., 2014, 50, 4238–4241.
- Salles et al. 2010 F. Salles, G. Maurin, C. Serre, P. L. Llewellyn, C. Knöfel, H. J. Choi, Y. Filinchuk, L. Oliviero, A. Vimont, J. R. Long and G. Férey, J. Am. Chem. Soc., 2010, 132, 13782–13788.
- Goodwin et al. 2008 A. L. Goodwin, D. A. Keen, M. G. Tucker, M. T. Dove, L. Peters and J. S. O. Evans, J. Am. Chem. Soc., 2008, 130, 9660–9661.
- Ho and Taylor 1998 C. Y. Ho and R. E. Taylor, Thermal Expansion of Solids, ASM International, Ohio, 1998.
- Krishnan 1979 Thermal Expansion of Crystals, ed. R. Krishnan, Pergamon Press, Oxford, 1979, vol. 22.
- Goodwin et al. 2008 A. L. Goodwin, D. A. Keen and M. G. Tucker, Proc. Natl. Acad. Sci., U.S.A., 2008, 105, 18708–18713.
- Greaves et al. 2011 G. N. Greaves, A. L. Greer, R. S. Lakes and T. Rouxel, Nature Mater., 2011, 10, 823–837.
- Panda et al. 2014 M. K. Panda, T. Runčevski, S. Chandra Sahoo, A. A. Belik, N. K. Nath, R. E. Dinnebier and P. Naumov, Nature Commun., 2014, 5, 4811.
- Panda et al. 2016 M. K. Panda, R. Centore, M. Causà, A. Tuzi, F. Borbone and P. Naumov, Sci. Rep., 2016, 6, 29610.
- Korčok et al. 2009 J. L. Korčok, M. J. Katz and D. B. Leznoff, J. Am. Chem. Soc., 2009, 131, 4866–4871.
- Ludi and Güdel 1968 A. Ludi and H. U. Güdel, Helv. Chim. Acta, 1968, 51, 1762–1765.
- Goodwin et al. 2008 A. L. Goodwin, M. Calleja, M. J. Conterio, M. T. Dove, J. S. O. Evans, D. A. Keen, L. Peters and M. G. Tucker, Science, 2008, 319, 794–797.
- Conterio et al. 2008 M. J. Conterio, A. L. Goodwin, M. G. Tucker, D. A. Keen, M. T. Dove, L. Peters and J. S. O. Evans, J. Phys.: Cond. Matt., 2008, 20, 255225.
- Baughman et al. 1998 R. H. Baughman, S. Stafström, C. Cui and S. O. Dantas, Science, 1998, 279, 1522–1524.
- Lakes and Wojciechowski 2008 R. Lakes and K. W. Wojciechowski, Phys. Stat. Sol. B, 2008, 245, 545–551.
- Cairns and Goodwin 2015 A. B. Cairns and A. L. Goodwin, Phys. Chem. Chem. Phys., 2015, 17, 20449–20465.
- Jansen 1987 M. Jansen, Angew. Chem. Int. Ed. Engl., 1987, 26, 1098–1110.
- O’Grady and Kaltsoyannis 2004 E. O’Grady and N. Kaltsoyannis, Phys. Chem. Chem. Phys., 2004, 6, 680–687.
- Bessler et al. 2001 K. E. Bessler, L. A. d. P. Calzavara, V. M. Deflon and E. Niquet, Acta Cryst. E, 2001, 57, m522–m523.
- Hunt et al. 2016 S. J. Hunt, M. J. Cliffe, J. A. Hill, A. B. Cairns, N. P. Funnell and A. L. Goodwin, CrystEngComm, 2016, 17, 361–369.
- Pauling and Pauling 1968 L. Pauling and P. Pauling, Proc. Natl. Acad. Sci., 1968, 60, 362–367.
- Haser et al. 1977 R. Haser, B. Bonnet and J. Roziere, J. Mol. Structure, 1977, 40, 177–189.
- Keen et al. 2010 D. A. Keen, M. T. Dove, J. S. O. Evans, A. L. Goodwin, L. Peters and M. G. Tucker, J. Phys.: Cond. Matt., 2010, 22, 404202.
- Lind 2012 C. Lind, Materials, 2012, 5, 1125–1154.
- 36 A. A. Coelho, TOPAS-Academic, version 4.1 (computer software), Coelho software technical report.
- Stephens 1999 P. W. Stephens, J. Appl. Cryst., 1999, 32, 281–289.
- Ratuszna and Małecki 2000 A. Ratuszna and G. Małecki, Mater. Sci. Forum, 2000, 321-324, 947–953.
- Cliffe and Goodwin 2012 M. J. Cliffe and A. L. Goodwin, J. Appl. Cryst., 2012, 45, 1321–1329.
- Nye 1957 J. F. Nye, Physical Properties of Crystals, Oxford University Press, Oxford, 1957.
- Blum et al. 2009 V. Blum, R. Gehrke, F. Hanke, P. Havu, V. Havu, X. Ren, K. Reuter and M. Scheffler, Comput. Phys. Commun., 2009, 180, 2175–2196.
- Perdew et al. 1996 J. P. Perdew, K. Burke and M. Ernzerhof, Phys. Rev. Lett., 1996, 77, 3865–3868.
- Tkatchenko and Scheffler 2009 A. Tkatchenko and M. Scheffler, Phys. Rev. Lett., 2009, 102, 073005.
- Tkatchenko et al. 2012 A. Tkatchenko, R. A. DiStasio Jr., R. Car and M. Scheffler, Phys. Rev. Lett., 2012, 108, 236402.
- Ambrosetti et al. 2014 A. Ambrosetti, A. M. Reilly, R. A. DiStasio Jr. and A. Tkatchenko, J. Chem. Phys., 2014, 140, 18A508.
- Gale 1997 J. D. Gale, J. Chem. Soc., Faraday Trans., 1997, 93, 629–637.
- Sharpe 1976 A. G. Sharpe, The Chemistry of Cyano Complexes of the Transition Metals, Academic Press, London, 1976.
- Widmann et al. 2005 A. Widmann, H. Kahlert, H. Wulff and F. Scholz, J. Sol. St. Electrochem., 2005, 9, 380–389.
- Hibble et al. 2004 S. J. Hibble, S. G. Eversfield, A. R. Cowley and A. M. Chippindale, Angew. Chem. Int. Ed., 2004, 43, 628–630.
- Katz et al. 2008 M. J. Katz, K. Sakai and D. B. Leznoff, Chem. Soc. Rev., 2008, 37, 1884–1895.
- Batsanov 2001 S. S. Batsanov, Inorg. Mater., 2001, 37, 871–885.
- Hibble et al. 2010 S. J. Hibble, G. B. Wood, E. J. Biblé, A. H. Pohl, M. G. Tucker, A. C. Hannon and A. M. Chippindale, Z. Krist., 2010, 225, 457–462.
- Mary et al. 1996 T. A. Mary, J. S. O. Evans, T. Vogt and A. W. Sleight, Science, 1996, 272, 90–92.
- Barerra et al. 2005 G. D. Barerra, J. A. O. Bruno, T. H. K. Barron and N. L. Allan, J. Phys.: Cond. Matt., 2005, 17, R217–R252.
- Adamson et al. 2015 J. Adamson, T. C. Lucas, A. B. Cairns, N. P. Funnell, M. G. Tucker, A. K. Kleppe, J. A. Hriljac and A. L. Goodwin, Physica B, 2015, 479, 35–40.
- Liu et al. 2016 X. Liu, J. Hermann and A. Tkatchenko, J. Chem. Phys., 2016, 145, 241101.
- Calleja et al. 2008 M. Calleja, A. L. Goodwin and M. T. Dove, J. Phys.: Cond. Matt., 2008, 20, 255226.
- Magnko et al. 2002 L. Magnko, M. Schweizer, G. Rauhut, M. Schutz, H. Stoll and H.-J. Werner, Phys. Chem. Chem. Phys., 2002, 4, 1006–1013.
- Fang et al. 2014 H. Fang, M. T. Dove and K. Refson, Phys. Rev. B, 2014, 90, 054302.
- Munn 1972 R. W. Munn, J. Phys. C.: Solid State Phys., 1972, 5, 535–542.
- Dove and Fang 2016 M. T. Dove and H. Fang, Rep. Prog. Phys., 2016, 79, 066503.
- Hermet et al. 2013 P. Hermet, J. Catafesta, J.-L. Bantignies, C. Levelut, D. Maurin, A. B. Cairns, A. L. Goodwin and J. Haines, J. Phys. Chem. C, 2013, 117, 12848–12857.
- Hill et al. 2015 J. A. Hill, A. B. Cairns, J. J. K. Lim, S. J. Cassidy, S. J. Clarke and A. L. Goodwin, CrystEngComm, 2015, 17, 2925–2928.
- Cairns et al. 2013 A. B. Cairns, J. Catafesta, C. Levelut, J. Rouquette, A. van der Lee, L. Peters, A. L. Thompson, V. Dmitriev, J. Haines and A. L. Goodwin, Nature Mater., 2013, 12, 212–216.
- Woodall et al. 2013 C. H. Woodall, C. M. Beavers, J. Christensen, L. E. Hatcher, M. Intissar, A. Parlett, S. J. Teat, C. Reber and P. R. Raithby, Angew. Chem. Int. Ed., 2013, 52, 9691–9694.