We report a single-crystal nanoindentation study of the negative compressibility material zinc(II) dicyanoaurate. The material exhibits a particularly strong elastic recovery, which we attribute to the existence of supramolecular helices that function as atomic-scale springs—storing mechanical energy during compressive stress and inhibiting plastic deformation. Our results are consistent with the relationship noted by Cheng and Cheng [Appl. Phys. Lett. 73, 614 (1998)] between elastic recovery and the ratio of material hardness to Young’s modulus. Drawing on comparisons with other framework materials containing helical motifs, we suggest helices as an attractive design element for imparting resistance to plastic deformation in functional materials.

Framework materials—solids with structures assembled from a combination of nodes and linkers—are well-known to harbour a variety of interesting mechanical properties.1 This is especially true for systems with long and flexible linkers, since flexibility generally amplifies mechanical response.2 One topical example is the broad family of metal–organic frameworks (MOFs), members of which are known to exhibit, e.g., anomalously low elastic moduli,3 extreme mechanical anisotropy,4 negative Poisson’s ratios,5 and a propensity for negative linear compressibility (NLC).6 Taken at face value, these properties suggest a range of attractive applications for open frameworks such as MOFs and MOF-like systems, including in stimuli-responsive materials (sensors), actuators, and shock absorbers.7,8

Fundamental studies of framework mechanical response often (rightly) focus on elastic behaviour, and yet the relevance of the elastic regime to practical applications relies on the extent to which plastic deformation can be avoided under application of stress. Often the origin of this stress is an external pressure (e.g., hydrostatic or uniaxial), but guest sorption is an alternative mechanism of particular relevance to MOFs.9 A key measure of resistance to irreversible deformation is the elastic recovery, We, which quantifies the proportion of work resulting in elastic deformation.10 Put simply, materials with high elastic recoveries will regain their structure and shape after stress cycling: the elastic recovery of rubber, for example, is nearly unity, whereas brittle ceramics have values of We close to zero. Consequently there is clear merit in developing an understanding of the elastic recovery of open frameworks so as to help predict the extent to which their anomalous mechanical responses might be exploited in practice.

Here we study the elastic recovery of zinc(II) dicyanoaurate, Zn[Au(CN)2]2, a MOF-like system known for its extreme NLC behaviour.11 Whereas the uniaxial compressibilities K=(/p)T of conventional materials are usually in the range of 5–10 TPa−1 (i.e., a linear contraction of 0.5%–1% for each 1 GPa of applied pressure), Zn[Au(CN)2]2 exhibits a remarkably large and negative compressibility along the hexagonal axis of its quartzlike structure, with Kc = −42(5) TPa−1.12 In other words, its structure expands on compression by roughly 4% for each GPa. This extraordinary elastic property places Zn[Au(CN)2]2 in a unique position for application in pressure sensors and shock amplifiers.13 Yet remarkably little is known regarding its mechanical behaviour beyond its (elastic) compressibility and thermal expansivity.14 As an experimental technique, nanoindentation methods have provided crucial insight into the mechanical properties of flexible MOFs.15–17 Consequently we sought to determine the experimental nanoindentation moduli, hardness, and elastic recovery of Zn[Au(CN)2]2. Our key result is the discovery of a large elastic recovery, which we attribute to the same supramolecular motifs thought to be responsible for NLC itself.

Before presenting our results we summarise briefly the structural chemistry of Zn[Au(CN)2]2. Its crystal structure was first reported in Ref. 18, where the relationship to the structure of SiO2-quartz was clearly noted [Fig. 1(a)]. Like that of β-quartz, the crystal symmetry of Zn[Au(CN)2]2 is hexagonal, and so from a mechanical perspective the important crystal directions are the [100] and [001] axes. The quartzlike nets of Zn[Au(CN)2]2 are assembled from tetrahedral Zn2+ centres, connected via almost-linear dicyanoaurate ([Au(CN)2]) linkers. In this structure, the ZnZn distance is so large (ca 1 nm) that six separate quartzlike nets interpenetrate.19 While each net is covalently distinct, neighbouring nets interact via aurophilic (AuAu) interactions, the existence of which is inferred from the relatively short AuAu contacts.20 These interactions connect to form a set of aurophilic helices lying along directions perpendicular to the hexagonal axis. Variable-pressure crystallographic measurements suggest that the extreme compressibility of this system arises because these aurophilic helices function as supramolecular “springs”.6 Just as a steel spring is more compressible than steel itself, so too are the aurophilic helices in Zn[Au(CN)2]2 more compressible than the AuAu interactions from which they are made. We will come to show that this same supramolecular motif may also be responsible for a strong elastic recovery of the material in these directions.

FIG. 1.

(a) Representation of the crystal structure of Zn[Au(CN)2]2. Zn atoms in gray, C in black, and N in blue, and Au atoms in gold. The top view is parallel to [001]; the bottom view is parallel to [11¯0]. (b) Aurophilic interactions between neighbouring frameworks form helices oriented perpendicular to the hexagonal c-axis, such as the one shown here in red. (c) The compound crystallises as elongated hexagonal bipyramids. The well-defined crystal facets make Zn[Au(CN)2]2 an ideal candidate for nanoindentation studies.

FIG. 1.

(a) Representation of the crystal structure of Zn[Au(CN)2]2. Zn atoms in gray, C in black, and N in blue, and Au atoms in gold. The top view is parallel to [001]; the bottom view is parallel to [11¯0]. (b) Aurophilic interactions between neighbouring frameworks form helices oriented perpendicular to the hexagonal c-axis, such as the one shown here in red. (c) The compound crystallises as elongated hexagonal bipyramids. The well-defined crystal facets make Zn[Au(CN)2]2 an ideal candidate for nanoindentation studies.

Close modal

For our study, we prepared single crystals of Zn[Au(CN)2]2 using the hydrothermal synthesis approach reported in Ref. 11 (see supplementary material for details). Nanoindentation measurements were carried out at ambient temperature using an MTS NanoIndenter©XP, equipped with a continuous stiffness measurement (CSM) module. The instrument was placed within an isolation cabinet that shielded against thermal instability and acoustic interference. A three-sided pyramidal Berkovich indenter with a sharp tip (end radius 50 nm) was used to measure the indentation modulus and hardness.21 Calibration was performed using a fused silica standard, with elastic modulus (E = 72 GPa) and hardness (H = 9 GPa). Thermal drifts were ensured to be consistently low. Indentation depth was 2000 nm for all measurements and the inter-indent spacing was typically 35–40 μm to avoid interference between neighbouring indentations. We probed two orientations of Zn[Au(CN)2]2 crystals, cleaved and polished so as to expose (100) and (001) crystal faces. Continuous stiffness measurement allowed us to collect depth-dependent mechanical data. The indentation modulus was determined using the Oliver and Pharr method;22 the hardness was determined by the ratio of the applied pressure and contact area at maximum indentation.

Typical load–displacement curves for Zn[Au(CN)2]2 are shown in Fig. 2(a). The particular shape of these curves reflects the pyramidal nature of the Berkovich tip used in our study, which results in an increasing contact area at higher loads—and hence decreasing variation in displacement depth on indentation. Our data are characteristic of the nanoindentation response of framework materials, showing contributions from both elastic and plastic deformation;22 the latter is evident in that loading and unloading curves are not coincident.

FIG. 2.

Typical load–displacement plots obtained from nanoindentation measurements of Zn[Au(CN)2]2 using the (100) and (001) facets. The elastic recovery We can be calculated from the ratio of the areas under the unloading and loading curves, as described in the text. The inset shows the effective elastic moduli extracted from these data, given as a function of indentation depth h. (b) Microscopy images of the (100) and (001) surfaces (left and right, respectively) following nanoindentation to a common maximum penetration depth hmax. Note that the degree of residual deformation is much greater for (001) than for (100); this is consistent with the smaller elastic recovery of the former relative to the latter.

FIG. 2.

Typical load–displacement plots obtained from nanoindentation measurements of Zn[Au(CN)2]2 using the (100) and (001) facets. The elastic recovery We can be calculated from the ratio of the areas under the unloading and loading curves, as described in the text. The inset shows the effective elastic moduli extracted from these data, given as a function of indentation depth h. (b) Microscopy images of the (100) and (001) surfaces (left and right, respectively) following nanoindentation to a common maximum penetration depth hmax. Note that the degree of residual deformation is much greater for (001) than for (100); this is consistent with the smaller elastic recovery of the former relative to the latter.

Close modal

Indentation moduli were determined using continuous stiffness measurements on loading, collected for at least twelve crack-free indents for both (100) and (001) crystal faces [inset to Fig. 2(a); see supplementary material]. Using a representative value of the effective isotropic Poisson’s ratio (ν=0.2), we obtain the values E100 = 14.5(3) GPa and E001 = 21.4(5) GPa; these values vary by not more than 30% when recalculated with alternate values of ν (see supplementary material).

Indentation moduli are not direct measurements of the Young’s moduli, but we treat the former as a useful approximation to the latter for the remainder of our analysis (e.g., as in Ref. 22). The moduli are low in comparison with corresponding values for many other coordination polymers23 but remain higher than the values reported for most porous MOFs.24 Zn[Au(CN)2]2 is not porous, and yet its number density (ρ=0.0512 atoms Å−3) shares more in common with open framework materials such as the zeolitic imidazolate framework (ZIF) family of MOFs (ρ<0.08 atoms Å−3) than it does with non-porous frameworks (ρ>0.08 atoms Å−3). Hence, we suggest that the mechanical behaviour of Zn[Au(CN)2]2 is MOF-like, despite the absence of any appreciable pore volume in its structure. This borderline behaviour is consistent with the observation that the isostructural material Zn[Ag(CN)2]2 is known to act as a host for the inclusion of inorganic guests.25 

The elastic response of Zn[Au(CN)2]2 shows clear evidence of anisotropy: the elastic modulus along [100] is 32% lower than along [001]. This anisotropy is consistent with the orientation within the crystal structure of aurophilic helices—supramolecular motifs that might be expected to be particularly compliant. Ab initio calculations of the Young’s moduli arrive at the same conclusion, with E100 lower than E001 by 29% (see supplementary material for details).

Analysis of load–displacement curves measured during nanoindentation also yields effective hardness values H in addition to the moduli E. We find H100 = 2.23(4) GPa and H001 = 2.19(7) GPa, suggesting that hardness is essentially isotropic for Zn[Au(CN)2]2. The magnitude of these values is again intermediate to those for dense frameworks, on the one hand, and ZIFs, on the other hand.26 

Our key finding relates to the elastic recovery, We, which was determined from our nanoindentation measurements by comparing the area under the loading and unloading curves—measures of the work done on indentation or release, respectively10,27:

We=WelasticWtotal=hfhmaxPdh0hmaxPdh,
(1)

where hmax is the maximum tip displacement, hf the final tip displacement, and P(h),P(h) represent the load applied during loading and unloading cycles at displacement h.

Figure 2(a) shows clearly that the elastic recovery for indentation in the [100] direction is greater than for the [001] direction. The numerical values we obtain using Eq. (1) are We = 62.7(3)% and 43.6(6)% for the (100) and (001) faces, respectively. These values are consistent with the recovery observed directly by inspection of the crystal post-indentation, as shown explicitly in Fig. 2(b). The (100) face shows a much higher resistance to irreversible deformation: the indents are clearly more shallow and show less cracking compared to those for the (001) face. A high elastic recovery leads to improved fracture toughness—i.e., resistance to cracking.28 Likewise, the degree of material pile-up around residual indents is much greater for the (001) face than for the (100) face, which is consistent with the amplification of material uplift as a result of the formation of cracks.

It is known that there is a linear empirical relationship between We and the ratio H/E of hardness to Young’s modulus.10,22,29 This relationship allows us to place the elastic recovery behaviour of Zn[Au(CN)2]2 in the context of that of other MOFs and coordination polymers, by exploiting published H, E data for the ZIF family30 as well as some dense hybrid frameworks: copper phosphonoacetate (polymorphs 1 and 2, referred to here as CuPA1 and CuPA2),28 cerium oxalate formate (CeOx),31 and zinc phosphate phosphonoacetate hydrate (ZnPPA).32 Our results are shown graphically in Fig. 3, where we have also included in this figure the original data of Ref. 10. What is clear is that the We value we measure for Zn[Au(CN)2]2 along [100] is amongst the largest for all these systems and is comparable to that of fused silica, which is widely used as a high-recovery standard in nanoindentation studies. The behaviour of Zn[Au(CN)2]2 along [001] for Zn[Au(CN)2]2 is typical of a dense hybrid framework, such that the system as a whole behaves like a MOF along a but like a dense framework along c.

FIG. 3.

Empirical relationship between elastic recovery We and H/E. Data are given for Zn[Au(CN)2]2 (green symbols), a variety of other porous (red symbols), and dense (blue symbols) framework materials. Also shown are the corresponding data for conventional materials as collated in Ref. 10 (black and white symbols). The data fall naturally into three clusters: systems with low elastic recoveries (brittle ceramics), those with moderate values of We (dense frameworks), and a small set of materials with high, polymer-like elastic recoveries. This last set includes many of the porous ZIFs. Zinc dicyanoaurate behaves as if a dense framework when stress is applied along the hexagonal c axis but as if an open MOF-type system when compressed along [100]. The insets show a representation of the structures of ZIF-zni as discussed in the text.

FIG. 3.

Empirical relationship between elastic recovery We and H/E. Data are given for Zn[Au(CN)2]2 (green symbols), a variety of other porous (red symbols), and dense (blue symbols) framework materials. Also shown are the corresponding data for conventional materials as collated in Ref. 10 (black and white symbols). The data fall naturally into three clusters: systems with low elastic recoveries (brittle ceramics), those with moderate values of We (dense frameworks), and a small set of materials with high, polymer-like elastic recoveries. This last set includes many of the porous ZIFs. Zinc dicyanoaurate behaves as if a dense framework when stress is applied along the hexagonal c axis but as if an open MOF-type system when compressed along [100]. The insets show a representation of the structures of ZIF-zni as discussed in the text.

Close modal

We suggest it is the supramolecular helix motif—instrumental in allowing NLC6—that is also responsible for the remarkably high elastic recovery of Zn[Au(CN)2]2 along [100]. As already discussed, the helix is an object with the particular mechanical property of being able to accommodate bulk compression without substantial compression of the material from which it is constructed. Hence compression of Zn[Au(CN)2]2 along the [100] direction can in principle be accommodated without substantial changes in any bonding interactions within the solid. In this way, the strain induced during nanoindentation is less likely to result in bond breaking and fracture than might otherwise be the case. Instead the supramolecular springs effectively store the energy transferred during indentation (i.e., work done) and release this energy reversibly upon removal of the indenter tip.

If this is indeed the case for Zn[Au(CN)2]2, then it is natural to question whether similar supramolecular motifs may be responsible for strong elastic recoveries in other materials. The data of Ref. 30 reveal that ZIF-zni shows an elastic recovery (We = 69%) that is comparable to that found in Zn[Au(CN)2]2. At face value, it is perhaps unexpected that this particular member of the ZIF family should show a particularly large value of We since it is both the densest and the least compliant ZIF: its Young’s moduli fall in the range E = 7.5–8.5 GPa.33,34 The material crystallises in the tetragonal space-group I41cd and within the network are chains of Zn2+ cations and imidazolate linkers that run along the c-axis, forming a triple-helical structure as shown in the inset to Fig. 3. By analogy with the arguments presented above, we suggest that this triple helix would also provide a high resistance to permanent deformation. Indeed this triple-helix motif is also found in collagen—crucial supporting components of cartilage, ligaments, tendons, bone, and skin—which itself has a high elastic recovery.35 

Of course, there are likely to be many other supramolecular motifs—such as hydrogen bonding, or low-dimensional components—that can allow efficient storage of mechanical energy during compression. Nevertheless an interesting future avenue of research would be to explore the elastic recovery properties of helical systems in a broader sense, including elemental solids such as selenium,36,37 liquid crystals,38 and supramolecular assemblies.39 Our key result is that the chemical motif responsible for interesting elastic response (NLC) in Zn[Au(CN)2]2 also confers resistance to plastic deformation, establishing helical structures as an attractive motif in the design of resilient materials.

See supplementary material for a detailed analysis of our nanoindentation data together with a summary of our ab initio density functional theory calculations.

C.S.C. and A.L.G. gratefully acknowledge financial support from the E.R.C. (No. 279705) and the Leverhulme Trust (Grant No. RPG-2015-292). J.-C.T. acknowledges research funding from the EPSRC (Nos. EP/K031503/1 and EP/N014960/1). M.R.R. thanks the EPSRC for a DTA postgraduate scholarship, the STFC for a CMSD Award (13-05), and the Rutherford Appleton Laboratory for access to the SCARF cluster.

1.
B. F.
Hoskins
and
R.
Robson
,
J. Am. Chem. Soc.
112
,
1546
(
1990
).
2.
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.
136
,
7801
(
2014
).
3.
J. C.
Tan
,
B.
Civalleri
,
C.-C.
Lin
,
L.
Valenzano
,
R.
Galvelis
,
P.-F.
Chen
,
T. D.
Bennett
,
C.
Mellot-Draznieks
,
C. M.
Zicovich-Wilson
, and
A. K.
Cheetham
,
Phys. Rev. Lett.
108
,
095502
(
2012
).
4.
A. U.
Ortiz
,
A.
Boutin
,
A. H.
Fuchs
, and
F.-X.
Coudert
,
Phys. Rev. Lett.
109
,
195502
(
2012
).
5.
G. N.
Greaves
,
A. L.
Greer
,
R. S.
Lakes
, and
T.
Rouxel
,
Nat. Mater.
10
,
823
(
2011
).
6.
R. H.
Baughman
,
S.
Stafström
,
C.
Cui
, and
S. O.
Dantas
,
Science
279
,
1522
(
1998
).
7.
A. E.
Aliev
,
J.
Oh
,
M. E.
Kozlov
,
A. A.
Kuznetsov
,
S.
Fang
,
A. F.
Fonseca
,
R.
Ovalle
,
M. D.
Lima
,
M. H.
Haque
,
Y. N.
Gartstein
,
M.
Zhang
,
A. A.
Zakhidov
, and
R. H.
Baughman
,
Science
323
,
1575
(
2009
).
8.
G. M.
Spinks
,
G. G.
Wallace
,
L. S.
Fifield
,
L. R.
Dalton
,
A.
Mazzoldi
,
D.
De Rossi
,
I. I.
Khayrullin
, and
R. H.
Baughman
,
Adv. Mater.
14
,
1728
(
2002
).
9.
S.
Krause
,
V.
Bon
,
I.
Senkovska
,
U.
Stoeck
,
D.
Wallacher
,
D. M.
Többens
,
S.
Zander
,
R. S.
Pillai
,
G.
Maurin
,
F.-X.
Coudert
, and
S.
Kaskel
,
Nature
532
,
348
(
2016
).
10.
Y.-T.
Cheng
and
C.-M.
Cheng
,
Appl. Phys. Lett.
73
,
614
(
1998
).
11.
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
,
Nat. Mater.
12
,
212
(
2013
).
12.
H. L.
Zhang
,
S.
Lu
,
M. P. J.
Punkkinen
,
Q.-M.
Hu
,
B.
Johansson
, and
L.
Vitos
,
Phys. Rev. B
82
,
132409
(
2010
).
13.
Z. G.
Nicolaou
and
A. E.
Motter
,
Nat. Mater.
11
,
608
(
2012
).
14.
A. L.
Goodwin
,
B. J.
Kennedy
, and
C. J.
Kepert
,
J. Am. Chem. Soc.
131
,
6334
(
2009
).
15.
J. P.
Best
,
J.
Michler
,
J.
Liu
,
Z.
Wang
,
M.
Tsotsalas
,
X.
Maeder
,
S.
Röse
,
V.
Oberst
,
J.
Liu
,
S.
Walheim
,
H.
Gliemann
,
P. G.
Weidler
,
E.
Redel
, and
C.
Wöll
,
Appl. Phys. Lett.
107
,
101902
(
2015
).
16.
X.-L.
Qu
,
D.
Gui
,
Z.-L.
Zheng
,
R.
Li
,
H.-L.
Han
,
X.
Li
, and
P.-Z.
Li
,
Dalton Trans.
45
,
6983
(
2016
).
17.
C. L.
Hobday
,
R. J.
Marshall
,
C. F.
Murphie
,
J.
Sotelo
,
T.
Richards
,
D. R.
Allan
,
T.
Düren
,
F.-X.
Coudert
,
R. S.
Forgan
,
C. A.
Morrison
,
S. A.
Moggach
, and
T. D.
Bennett
,
Angew. Chem., Int. Ed.
55
,
2401
(
2016
).
18.
B. F.
Hoskins
,
R.
Robson
, and
N. V. Y.
Scarlett
,
Angew. Chem., Int. Ed.
34
,
1203
(
1995
).
19.
S. R.
Batten
,
Curr. Opin. Solid State Mater. Sci.
5
,
107
(
2001
).
20.
M. J.
Katz
,
T.
Ramnial
,
H.-Z.
Yu
, and
D. B.
Leznoff
,
J. Am. Chem. Soc.
130
,
10662
(
2008
).
21.
E. S.
Berkovich
,
Ind. Diamond Rev.
11
,
129
(
1951
).
22.
W. C.
Oliver
and
G. M.
Pharr
,
J. Mater. Res.
19
,
3
(
2004
).
23.
Z.
Zhang
,
X.
Jiang
,
G.
Feng
,
Z.
Lin
,
B.
Hu
, and
W.
Li
,
J. Solid State Chem.
233
,
289
(
2016
).
24.
D. F.
Bahr
,
J. A.
Reid
,
W. M.
Mook
,
C. A.
Bauer
,
R.
Stumpf
,
A. J.
Skulan
,
N. R.
Moody
,
B. A.
Simmons
,
M. M.
Shindel
, and
M. D.
Allendorf
,
Phys. Rev. B
76
,
184106
(
2007
).
25.
J. A.
Hill
,
A. L.
Thompson
, and
A. L.
Goodwin
,
J. Am. Chem. Soc.
138
,
5886
(
2016
).
26.
J. C.
Tan
and
A. K.
Cheetham
,
Chem. Soc. Rev.
40
,
1059
(
2011
).
27.
T.
Miura
,
Y.
Benino
,
R.
Sato
, and
T.
Komatsu
,
J. Eur. Ceram. Soc.
23
,
409
(
2003
).
28.
J. C.
Tan
,
C. A.
Merrill
,
J. B.
Orton
, and
A. K.
Cheetham
,
Acta Mater.
57
,
3481
(
2009
).
29.
A.
Bolshakov
and
G. M.
Pharr
,
J. Mater. Res.
13
,
1049
(
1998
).
30.
J. C.
Tan
,
T. D.
Bennett
, and
A. K.
Cheetham
,
Proc. Natl. Acad. Sci. U. S. A.
107
,
9938
(
2010
).
31.
J. C.
Tan
,
J. D.
Furman
, and
A. K.
Cheetham
,
J. Am. Chem. Soc.
131
,
14252
(
2009
).
32.
M.
Kosa
,
J. C.
Tan
,
C. A.
Merrill
,
M.
Krack
,
A. K.
Cheetham
, and
M.
Parrinello
,
Chem. Phys. Chem.
11
,
2332
(
2010
).
33.
R.
Lehnert
and
F.
Seel
,
Z. Anorg. Allg. Chem.
464
,
187
(
1980
).
34.
T. D.
Bennett
,
A. L.
Goodwin
,
M. T.
Dove
,
D. A.
Keen
,
M. G.
Tucker
,
E. R.
Barney
,
A. K.
Soper
,
E. G.
Bithell
,
J. C.
Tan
, and
A. K.
Cheetham
,
Phys. Rev. Lett.
104
,
115503
(
2010
).
35.
K. E.
Aifantis
,
S.
Shrivastava
, and
G. M.
Odegard
,
J. Mater. Sci.: Mater. Med.
22
,
1375
(
2011
).
36.
D. R.
McCann
,
L.
Cartz
,
R. E.
Schmunk
, and
Y. D.
Harker
,
J. Appl. Phys.
43
,
1432
(
1972
).
37.
A. B.
Cairns
and
A. L.
Goodwin
,
Phys. Chem. Chem. Phys.
17
,
20449
(
2015
).
38.
Z.-G.
Zheng
,
Y.
Li
,
H. K.
Bisoyi
,
L.
Wang
,
T. J.
Bunning
, and
Q.
Li
,
Nature
531
,
352
(
2016
).
39.
E.
Yashima
,
N.
Ousaka
,
D.
Taura
,
K.
Shimomura
,
T.
Ikai
, and
K.
Maeda
,
Chem. Rev.
116
,
13752
(
2016
).

Supplementary Material