Qingxin
Zeng
*,
Wenbo
Qiu
,
Chengxi
Li
,
Yan
Sun
,
Jian
Hao
and
Yinwei
Li
*
Laboratory of Quantum Functional Materials Design and Application, School of Physics and Electronic Engineering, Jiangsu Normal University, Xuzhou 221116, China. E-mail: zengqx@jsnu.edu.cn; yinwei_li@jsnu.edu.cn
First published on 8th May 2024
A material with the “hidden” negative linear compressibility (NLC) will expand along a specific crystal direction upon uniformly compression to a critical pressure; such materials are thought to be promising candidates for non-linear actuators, switches and sensors. Herein, we use density functional theory (DFT) calculations to uncover the hidden NLC in a V-shaped molecular crystal, bis(5-amino-1,2,4-triazol-3-yl)methane (BATZM). The calculations indicate that the crystal is normally compressed over the pressure range of 0–3 GPa while it expands along the b-axis when the external hydrostatic pressure exceeds 3 GPa. The compressive behavior of the BATZM crystal is modulated by inter-molecular hydrogen bonds, which act as highly compressible springs at low pressures but robust struts at high pressures. Hence, the crystal prefers to compress the hydrogen bonds coupled with PLC at first and flatten the molecules, coupled with later NLC to resist the increasing external pressure. The compressive behavior of BATZM provides a strategy to design more hidden NLC materials via the rational use of the hydrogen bonds.
To date, known mechanisms for continuous NLC effects include (1) polyhedral rotation, (2) helices, (3) framework hinging and (4) layer sliding.1,11 Polyhedral rotation usually occurs in inorganic systems that contain polyhedral atom groups, which rotate under high pressure to give rise to elongation along one crystal direction, such as the NLC effect in YFe(CN)6, which arises from the rotation of YN6 units.12 The helix mechanism has been proposed in NLC effects in Se and Te.1,13 Metal–organic frameworks or molecular frameworks exhibiting NLC mostly follow the framework hinging mechanism, in which geometrical deformation compensates for the bond shrinkage in their crystal structure under high pressure. Wine-rack and deformed honey-comb motifs are two typical crystal structure models that give rise to NLC properties.14,15 The layer sliding mechanism is proposed in the case of Co(SCN)2(pyrazine)2, in which the pressure-driven sliding of 2D infinite-Co-pyrazine-layers leads to a moderate NLC effect.11
Most reported NLC materials exhibit NLC behavior at the beginning of compression at ambient conditions because they possess rigid struts in their crystal structures. However, a few “normal” materials also exhibit NLC phenomena only when compressed to a critical pressure without a phase transition; this phenomenon is known as “hidden” NLC.16–21 For example, the crystal of ammonium oxalate monohydrate shrinks along all directions at pressures below 5.1 GPa; thereafter, the c-axis begins to expand, which is explained by the pressure-induced strengthening of hydrogen bonds.21 Materials with hidden NLC are thought to be applicable as non-linear actuators, switches and sensors, as the NLC response is “off” at low pressures but can be switched “on” above a certain critical pressure.16 However, hidden NLC materials are very rare. Therefore, it is a vital issue to seek out more materials with these properties.
Our strategy to acquire hidden NLC materials is to focus on materials with weak interactions (such as van der Waals interactions and hydrogen bonds), which will be easily compressed at low pressure but could potentially be strengthened to a hard-to-compress level by pressure. The V-shaped organic molecule bis(5-amino-1,2,4-triazol-3-yl)methane (BATZM, C5H8N8) aggregates to form a crystal via intermolecular N–H⋯N hydrogen bonds. The V-shaped molecule itself is naturally anisotropic and would exhibit an abnormal compressibility in one direction by reason of the closing/opening of the “V” angle driven by pressure. However, as the connector between two molecules, the original weak hydrogen bond acts as a buffer to help release the stress applied on the crystal, which delays the potential abnormal compressive behavior of the BATZM crystal to exhibit hidden NLC behavior. In this work, we have investigated the elastic constants and compressive behavior of the BATZM crystal under high pressure using density functional theory (DFT) calculations. The results show that BATZM exhibits a hidden NLC effect along the b-axis modulated by the hydrogen bond at pressures higher than 3 GPa. The findings demonstrate the feasibility of our strategy to identify materials with so-called hidden NLC, which should shed light on the discovery of other materials with such a fantastic property.
The experimental crystal structure data were used as the beginning model for full structural relaxation to obtain the ground state for BATZM at 1 atm. To take into account the effect of the hydrogen bonds on the compressive behavior of this organic material, different van der Waals interaction correction methods were used for lattice relaxions at 0 GPa, as shown in Table 1. The DFT dispersion correction with Becke–Johnson damping (DFT-D3(BJ)) van der Waals scheme gives the best description of the lattice parameters (especially the b-axis) and intermolecular N4⋯N4 distance in the N4–H4B⋯N4 hydrogen bond. DFT-D3(BJ) is selected for further high-pressure calculations. The calculated lattice parameters a = 18.3938, b = 19.9132 and c = 4.1998 Å at 0 GPa are in excellent agreement with experimentally reported values, with discrepancies of −1.28%, 0.10% and 2.54%, respectively.
a (Å) | b (Å) | c (Å) | N4⋯N4 distance (Å) | |
---|---|---|---|---|
Expt. | 18.6320 | 19.9330 | 4.3095 | 3.204 |
DFT-GGA | 18.4504 | 20.2671 | 4.4751 | 3.323 |
DFT-D2 | 18.3641 | 20.5132 | 4.3157 | 3.277 |
DFT-D3 | 18.3640 | 20.5264 | 4.3203 | 3.280 |
DFT-D3(BJ) | 18.3938 | 19.9132 | 4.1998 | 3.120 |
DFT-D4 | 18.3316 | 20.3767 | 4.4461 | 3.310 |
DFT-TS | 18.3659 | 20.5228 | 4.3122 | 3.277 |
DFT-TS/HI | 18.3606 | 20.5192 | 4.3184 | 3.279 |
DFT-dDsC | 18.3534 | 20.5052 | 4.3162 | 3.277 |
optB86-vdW | 18.3722 | 20.5277 | 4.3151 | 3.279 |
optPBE-vdW | 18.3947 | 20.5415 | 4.3217 | 3.282 |
optB88-vdW | 18.3700 | 20.5292 | 4.3158 | 3.279 |
To investigate the mechanical response of BATZM to pressure, we applied hydrostatic pressures from 0 to 6 GPa in full lattice relaxations of BATZM with the DFT-D3(BJ) dispersion correction. No phase transitions occurred over the whole pressure range. The pressure-dependent cell volume data were fitted to a very smooth V(P) curve utilizing the third-order Birch–Murnaghan equation of state, calculated using PASCal,29 as depicted in Fig. 2a. The fitting yields the zero-pressure bulk modulus B0 with a value of 23.07 GPa, which is comparable to those of other hydrogen-bonded materials such as 3-methyl-4-nitropyridine N-oxide (B0 = 23 GPa for phase-I).30Fig. 2b shows the evolution of the relative change rates of the lattice parameters as a function of pressure. It is clearly seen that the BATZM crystal exhibits highly anisotropic compression, with the c-axis undergoing a remarkable reduction of about 12% up to 6 GPa, while the a-axis decreased by about 2%. Fantastically, below 3 GPa, the length of the b-axis decreased normally with increasing pressure, exhibiting a greater compressibility than that of the a-axis, but began to increase abnormally when the pressure was raised above 3 GPa, giving rise to an anomalous NLC behavior.
The online tool PASCal calculates linear compressibility based on the l–p curve fitted according to an empirical expression with the form l = l0 + λ(p − pc)ν,29 which is only applicable to those cases in which the length of the lattice parameter l changes monotonically as a function of pressure over the whole fitting pressure range. As the b-axis exhibited an increase after the initial decrease above 3 GPa, PASCal is inapplicable to calculate the linear compressibility of the b-axis, and hence we fitted the l–p curve to the fifth-order polynomial function l = a0 + a1p + a2p2 + ⋯ to obtain a highly smooth l–p curve,31,32 as shown in Fig. S1 (ESI†), for linear compressibility calculations on the basis of −dl/(ldp). As depicted in Fig. 2c, the linear compressibility of the lattice parameters indicates a large mechanical anisotropy in the structure. The a-axis has a nearly constant compressibility with the value varying around 4 TPa−1 over the whole compression process. The small Ka indicates that the a-axis is relatively hard to compress. The compressibility of the c-axis at 0 GPa is initially large with a value of about 26 TPa−1, in accordance with the rapid shortening of inter-molecular distances. It is worth stressing that the compressibility of the b-axis is positive at the beginning of the compression process; however, it becomes negative after compression to a simulated pressure of about 3.1 GPa, as shown in Fig. 2c; hence, the b-axis increases with pressure in the pressure range of 3 to 6 GPa. The mean compressibility along the b-axis calculated using the formula Kl = −Δl/(l0Δp) over the pressure range of 3–6 GPa is −2.0 TPa−1, which is comparable to that of many known NLC materials, such as methanol monohydrate (Ka = −3.1 TPa−1),32 silver(I) tricyanomethanide (Ka = −3.5 TPa−1, Kc = −4.0 TPa−1),33 silver(I) 2-methylimidazolate (K∼c = −4.32 TPa−1),34 ammonium oxalate monohydrate (Kb = −2.3 TPa−1)21 and perovskite MOFs [NH4][Zn(HCOO)3] (Kc = −1.8 TPa−1)35 and [NH2NH3][Co(HCOO)3] (Kb = −2.3 TPa−1).36 It is worth emphasizing that the PLC-to-NLC axis in other materials has the minimum change rate in length (namely the hardest axis) at low pressure.16,18–21,36,37 Taking ammonium oxalate monohydrate as an example, the PLC-to-NLC axis (b-axis) slightly contracted by about 0.1%, which was far less than the a-axis (∼2%) and c-axis (∼7%).19 In contrast, the b-axis in BATZM is not the hardest and contracts faster than the a-axis, as shown in Fig. 2b.
We then turned to the structural basic building unit (BBU)34 parameters to elucidate the origin of the NLC along the b-axis of BATZM from 3 to 6 GPa. As depicted in the inset of Fig. 3a, two fundamental parameters r and θ are assigned to the length of two arms (i.e., C3⋯N4 distance) of the symmetrical V-shaped motif and hinge angle (i.e. C3⋯N4⋯C3) between two edges, respectively. The changes in r and θ cooperatively determine the deformations of two dependent parameters l and d, where l represents the opening distance of V-shape characterized by the intra-molecular N4⋯N4 distance and d (parallel to the c-axis) describes the height of the V-shape. It is clearly seen in Fig. 3a that r and θ remain nearly unchanged in the low-pressure (LP) region, especially from 0 to approximately 1 GPa and begin to change faster after 1 GPa; however, in the high-pressure (HP) region, r remains nearly unchanged while θ undergoes a tremendous rise. Therefore, l and d are invariable in the LP region, while in the HP region, l and d undergo a colossal increase and decrease, respectively, following geometrical relationships l = 2rsin(θ/2) and d = rcos(θ/2). As d is parallel to the c-axis, and l is approximately parallel to the b-axis, with discrepancies of 4.96° and 7.08° at 0 GPa and 6 GPa, respectively, the changes in the lengths of the b- and c-axes can be ascribed to the changes in l and d, respectively. The shortening of d leads to the continuous shrinkage of the c-axis during the entire compression process, while the elongation of l in HP dominates the abnormal expansion of the b-axis.
The intermolecular hydrogen bonds were also monitored to comprehensively understand the variation in the lattice parameters. The hydrogen bonds in BATZM, including the intra- and inter-molecular ones, are all of the N–H⋯N type, as depicted in Fig. 1. Previous reports in the literature indicate that N–H⋯N bond lengths (N⋯N distances) vary from 2.526 to 3.793 Å and are distributed intensively around 2.98 Å.38–41 The intermolecular N4–H4B⋯N4 interaction within the bc-plane (shown as a green dashed line in Fig. 1b) with an initial large D⋯A distance of 3.1198 Å is quite weak and highly compressible. As shown in Fig. 3b, N4–H4B⋯N4 is compressed dramatically in the LP region and later becomes stronger and harder to compress in the HP region due to the pressure-driven shortening of the D⋯A distance. The persistence of l and θ indicates the absence of geometrical deformation in the BATZM molecule in the LP region; therefore, the shrinkage of the b-axis in the LP region is mainly caused by N4–H4B⋯N4 shortening. In the HP region, N4–H4B⋯N4 is too robust to shrink drastically; hence, the increase in strain to balance the elevated external pressure is mainly generated from the opening effect of the V-shaped molecule. Thus, the BBU parameter θ begins to increase and NLC along the b-axis occurs, as shown in Fig. 3a. N4–H4A⋯N2 and N1–H1⋯N3 (Fig. 1c) are arranged mainly along the a-axis and hence determine the evolution of the length of this axis under pressure. At ambient pressure, the N⋯N distance in N1–H1⋯N3 is 2.7071 Å, which makes the hydrogen bond robust enough to avoid suffering tremendous contractions. The nearly incompressible N1–H1⋯N3 finally gives rise to a modest contraction of the a-axis.
It can be seen from Fig. 3c that all of N–H bonds in the three types of N–H⋯N hydrogen bonds simultaneously elongate under pressure as a result of the increased electrostatic interaction between the acceptor N and H, which demonstrates the enhancement of the hydrogen bonds.42 The elongation of D–H segments in enhanced hydrogen bonds is very common in hydrogen-bonded organics such as biurea,43 glycinium oxalate44 and sulfamide.45 The Hirshfeld surface produced using CrystalExplorer was analysed to visualize intermolecular interactions, especially the hydrogen bonds.46 As shown in Fig. 3d, red spots indicate contacts shorter than van der Waals radii, while white and blue regions represent distances equal to and longer than the van der Waals radii, respectively. At high pressure, H⋯N contacts are apparently intensified, as revealed by expanded red spots. The most prominent change occurs in N4–H4B⋯N4 (the two ends of the molecule), which is nearly invisible at ambient pressure (Fig. 3d, top) but becomes a large red spot in the left/right side of the Hirshfeld surface at 6 GPa (Fig. 3d, bottom). The 2D fingerprint of the N–H distance shown in Fig. S2 (ESI†) exhibits a prominent shift toward the origin, which also implies the shortening and consequent strengthening of N⋯H nonbonds.
A schematic of the mechanism of the abnormal mechanical behavior is depicted in Fig. 4. To simplify the V-shaped molecule and while preserving its main characteristics, two solid lines are used to represent rigid planar triazol groups linking end N atoms from amino groups and hinge C atom from the substituted methane group together. The original N4–H4B⋯N4 hydrogen bond is described as a relaxed spring. In the LP region, the initial relaxed spring is relatively weaker than the rigid BATZM molecule; hence, the former is more sensitive to external pressure stimuli. Under compression, the molecule tends to exhibit drastic shrinkage in the length of spring but hold its molecular shape at the same time. The drastic compression of the spring strengthens its stiffness from ambient pressure up to the critical pressure; the stiff spring changes subtly, and the flattening of the V-shaped BATZM molecule starts to be activated by the external pressure stimulus under higher pressures. The flattening of the molecular shape exceeds the sum of decreases in bond length and intermolecular distance in amplitude, resulting in NLC behavior along the b-axis.
Fig. 4 Schematic of the mechanism of the abnormal mechanical behavior along the b-axis of the BATZM crystal. |
To investigate whether the hidden NLC is dispersion-correction-independent, we also performed full lattice relaxion without dispersion corrections over the pressure range of 0–10 GPa. The relative changes in the lattice parameters are shown in Fig. S3a (ESI†). The length of the b-axis first decreases and later starts to increase above 6 GPa, confirming the existence of the hidden NLC along the b-axis. However, the critical pressure for the PLC–NLC conversion is somewhat higher than obtained in the dispersion corrected results using DFT-D3(BJ). The change in the BBU parameter θ also exhibits such a delay: it changes very slightly over the pressure range of 0–3 GPa (Fig. S3b, ESI†) compared to the change in the range of 0–1 GPa in the DFT-D3(BJ) result. The uncorrected (by DFT-GGA) initial intermolecular N4⋯N4 distance at 0 GPa (3.323 Å) is longer than the DFT-D3(BJ) corrected one (3.120 Å), which plays as a key role in the conversion of PLC–NLC; the longer initial N4⋯N4 distance requires a higher pressure to compress it to an incompressible trend; hence, the uncorrected result shows a higher critical pressure (6 GPa) than the corrected result (3 GPa, DFT-D3(BJ)). Taking into account the fact that the experimental initial N4⋯N4 distance (3.204 Å at ambient condition) is intermediate between those of DFT-GGA and DFT-D3(BJ), the experimentally observed critical pressure for the conversion of PLC–NLC may presumably fall in between 3 and 6 GPa.
Footnote |
† Electronic supplementary information (ESI) available: Lattice parameter–pressure curve fitting, 2D fingerprint of N–H distance. See DOI: https://doi.org/10.1039/d4cp01102c |
This journal is © the Owner Societies 2024 |