Research Article | Open Access

Zhi-Gang Li, Kai Li, Li-Yuan Dong, Tian-Meng Guo, Muhammad Azeem, Wei Li, Xian-He Bu, "Acoustic Properties of Metal-Organic Frameworks", *Research*, vol. 2021, Article ID 9850151, 11 pages, 2021. https://doi.org/10.34133/2021/9850151

# Acoustic Properties of Metal-Organic Frameworks

#### Abstract

Metal-organic frameworks (MOFs) have attracted significant attention in the past two decades due to their diverse physical properties and associated functionalities. Although numerous advances have been made, the acoustic properties of MOFs have attracted very little attention. Here, we systematically investigate the acoustic velocities and impedances of 19 prototypical MOFs *via* first-principle calculations. Our results demonstrate that these MOFs exhibit a wider range of acoustic velocities, higher anisotropy, and lower acoustic impedances than their inorganic counterparts, which are ascribed to their structural diversity and anisotropy, as well as low densities. In addition, the piezoelectric properties, which are intimately related to the acoustic properties, were calculated for 3 MOFs *via* density functional perturbation theory, which reveals that MOFs can exhibit significant piezoelectricity due to the ionic contribution. Our work provides a comprehensive study of the fundamental acoustic properties of MOFs, which could stimulate further interest in this new exciting field.

#### 1. Introduction

Metal-organic frameworks (MOFs), constructed by metal ions/clusters and organic ligands *via* directional coordination bonds, have attracted wide attention during the past two decades since the homogeneous hybrid nature gives rise to the integration of properties from both the inorganic and organic components [1]. On the one hand, porous MOFs have shown promising application potential in gas storage and separation, catalysis, and drug delivery [2]; on the other hand, this emergent class of materials can exhibit diverse physical properties, including fluorescence, ferromagnetism, ferroelectricity, and multiferroicity [3–7]. Although there are numerous studies on the diverse properties and functionalities of MOFs [8], little attention has been paid to exploring an important area, namely, the acoustic properties [9]. Very recently, the application of acoustics and energy-transfer mechanisms in MOFs has been experimentally validated [10, 11], which prompts the necessity to systematically study the acoustic properties of MOFs. The most fundamental aspect of acoustics is to study the propagation of sound in materials, which includes the understanding of the acoustic velocity and acoustic impedance. Both of these parameters are closely related to the elastic properties of a given material and are determined by the stiffness constant and density [12]. Compared with purely inorganic (i.e., oxide) and organic (i.e., polymer) materials, MOFs have intermediate stiffness and density, which could enable them to show very different acoustic properties [13–15]. In this context, it is highly interesting to explore the fundamental acoustic velocity and acoustic impedance of MOFs and compare them with those of their purely inorganic and organic counterparts. In addition, it is well known that MOFs exhibit significant structural diversity and chemical variability; therefore, they are expected to show diverse acoustic properties, similar to their other properties [16]. Furthermore, the comprehensive study of the basic acoustic properties of MOFs is also of great importance for exploring their application potential in acoustics, such as acoustic switching devices and sensors.

In this work, we present a systematic study of the fundamental acoustic properties of 19 prototypical MOFs *via* first-principle calculations and compare them with those of traditional inorganic materials having similar structures. At the same time, potential acoustic applications of MOFs are proposed based on these properties. Moreover, the piezoelectric properties of several MOFs are calculated using density functional perturbation theory, and the advantageous properties of these MOFs over their inorganic counterparts are comprehensively discussed from the viewpoint of potential applications.

#### 2. Results

##### 2.1. Crystal Structure

To study the acoustic properties of MOFs, we calculated the acoustic quantities of 22 materials, which included 19 prototypical MOFs and 3 inorganic materials for comparison. According to the chemical compositions and structural features, these 19 MOFs can be classified into 7 categories, including isoreticular MOF-5 [Zn_{4}O(BDC)_{3}, ] [17]; ZIF-8 [Zn(mim)_{2}, ] and MAF-7 [Zn(mtz)_{2}, ] [18, 19]; DMOF-1sq and DMOF-1loz [Cu_{2}(BDC)_{2}(DABCO), ] [20]; MIL-47 [V^{IV}(O)(BDC)], MIL-53(Al)lp [Al(OH)(BDC)], and MIL-53(Ga)lp [Ga(OH)(BDC)] [21, 22]; QMOF-1 [Zn(ISN)_{2}, ] and QMOF-2 [InH(BDC)_{2}] [23]; metal formate frameworks [AB(HCOO)_{3}, where , MA (CH_{3}NH_{3}^{+}), Gua (C(NH_{2})_{3}^{+}), DMA ((CH_{3})_{2}NH_{2}^{+}), Hyz (NH_{2}NH_{3}^{+}), NH_{4}; , Co, Zn, Mg] [24–30]; and a unique dense organic-inorganic framework, DABCOH_{2}K(ClO_{4})_{3} [31]. In addition, the three traditional inorganic materials, quartz, BaTiO_{3}, and sodalite [Na_{4}Al_{3}(SiO_{4})_{3}], are selected for comparison with the structurally similar QMOFs, perovskite formates, and ZIF-8 and MAF-7, respectively. The structures of these materials are depicted in Figure 1.

##### 2.2. Acoustic Velocity

Generally, there are two types of acoustic waves: pressure waves (p-waves, longitudinal waves) and shear waves (s-waves, transverse waves). Here, we divide the velocity of shear waves into the velocities of pure-s-waves () and quasi-s-waves (). The quasi-s-waves have a small longitudinal wave component, which means that they are not purely transverse waves. The velocities of p-waves () and s-waves are determined by the Young’s moduli and shear moduli of the materials, respectively. As the Young’s moduli of most materials are greater than their shear moduli, the magnitudes of these wave velocities have the following trend: .

The calculated independent elastic constants of 15 materials by us and 7 materials by others all conform to the elastic stability criterion and are tabulated in Table S1 along with some experimentally measured densities [32–34]. Meanwhile, the acoustic velocities were obtained by solving the Christoffel equations based on the elastic constants and densities. The 3D surfaces of the acoustic velocities for these 22 materials are presented in Figs. S1-S21.

Figure 2 shows the maximum and minimum values of , , and in these 22 materials, which demonstrate that acoustic waves propagate in these MOF crystals with a wide range of velocities. For the maximum values of all three velocities, there are no obvious boundaries between the 19 MOFs and 3 inorganic materials. These results are reasonable considering the medium maximum moduli and low densities of MOFs. A similar trend is also found between the porous and dense MOFs. Notably, the maximum velocity values of some MOFs are even higher than those of quartz, BaTiO_{3},_{,} and sodalite. For example, the value of MIL-47 is 9.94 km/s, which is 138% larger than that of BaTiO_{3}. In terms of the minimum values of the three velocities, and are comparable with those of inorganic materials. However, the values of dense MOFs are comparable with those of inorganic materials and generally higher than those of porous MOFs. This is because the minimum shear moduli of some porous MOFs studied here are significantly lower than those of the dense systems and inorganic materials. For instance, MIL-53(Ga)lp has the lowest of 0.219 km/s among all 19 MOFs.

The ratios of the maximum and minimum velocities of each type of wave for these 22 materials were calculated as indicators to evaluate their anisotropy (), and the results are shown in Figure 3. Notably, the values of DMOF-1loz, DMOF-1sq, MIL-47, MIL-53(Al)lp, and MIL-53(Ga)lp are up to an order of magnitude higher than those of the 3 inorganic materials, which could give rise to some potential applications unavailable in conventional materials, as discussed below. In particular, MIL-53(Ga)lp has the highest acoustic velocity anisotropy, with anisotropy indices of the longitudinal wave velocity () and pure shear wave velocity () of 2.35 and 15.17, respectively.

The 3D surface contours and projected 2D plots of the acoustic velocities for MOF-5 are displayed in Fig. S1, which demonstrates that its , , and are in the ranges of 1.67-3.13, 1.67-3.67, and 5.50-6.83 km/s, respectively. The and of MOF-5 are approximately half of those of BaTiO_{3}, which also has a ReO_{3}-related structure, while the and are close to those of BaTiO_{3}. In addition, both the and values of MOF-5 approximate those of BaTiO_{3}. These findings disclose that acoustic waves can propagate similarly in porous MOFs and dense oxides. Furthermore, the anisotropy of the 3 types of wave velocities in MOF-5 is also distinct, and (1.19) is significantly smaller than (1.88) and (2.19). As the Young’s modulus of MOF-5 is determined by the connectivity of organic ligands, the isotropic periodic extension of the organic linkers in the structure leads to an isotropic Young’s modulus corresponding to a small . However, unlike the p-wave velocity, the s-wave velocity in MOF-5 is affected by the shear modulus. For the square-shaped building block in the structure, it is much more easily sheared along the organic linker direction (edge direction) than along the diagonal direction, hence resulting in dramatic anisotropy in and [32].

The 3D surface contours and 2D plots of the acoustic velocities for ZIF-8 and MAF-7 are displayed in Figs. S2 and S3. Although they are isomorphic with the sodalite structure, the , , and of ZIF-8 are all slightly smaller than those of MAF-7. According to formulas (1)–(3), the difference is attributed to the enhanced Young’s modulus of MAF-7 due to the electron-donating effect and its lower density (ZIF-8, 1.14 g/cm^{3}; MAF-7, 1.08 g/cm^{3}). In addition, the and values of ZIF-8 and MAF-7 are only approximately 26-31% of those of sodalite, while the values are approximately 53-62% of the value of sodalite. These phenomena arise from the more obvious difference in the Young’s moduli between ZIF-8/MAF-7 and sodalite than that in the shear moduli. Moreover, both ZIF-8 and MAF-7 exhibit very small anisotropy in sound velocities, which is reminiscent of that of sodalite.

DMOFs have unusual guest-dependent dynamic behaviour: the framework shrinks to a lozenge (DMOF-1loz) when the guest is present and expands to a square (DMOF-1sq) when the guest is released. Along with the change in unit cell volume from 1147 (DMOF-1sq) to 1114 Å^{3} (DMOF-1loz), the density increases from 0.826 to 0.850 g/cm^{3}. This drastic structural change also leads to a significant difference in framework stiffness and corresponding acoustic velocities (Figs. S4 and S5). The , , , and values of DMOF-1lsq are all larger than those of DMOF-1loz, while the and exhibit an inverse trend. Notably, the anisotropies in and of both the lozenge and square phases are very high, which are 7.55 and 11.43 for the former and 5.31 and 8.69 for the latter.

The members of the MIL family of MOFs have similar structures but different densities (MIL-47, , , 1.57 g/cm^{3}), which leads to acoustic velocities in the same sequence as the density (Figures 4(a) and 4(b), S6, S7). For the Al and Ga analogues of MIL-53, the strength difference in the shorter Al-O and longer Ga-O coordination bonds is responsible for their different moduli and corresponding acoustic velocities. In addition, the anisotropy in for all three frameworks is significantly high and is 7.69, 8.65, and 15.17 for MIL-47, MIL-53(Al)lp, and MIL-53(Ga)lp, respectively. In particular, the anisotropy in of MIL-53Ga is the highest among all 22 materials, and such a large anisotropy is even larger than those of many 2D crystals. In terms of potential applications, MIL-53(Ga)lp could be utilized as a new type of acoustic switch through which acoustic waves can travel approximately 14 times faster in one direction than in the other, as shown in Figure 4(c).

QMOF-1 and QMOF-2 are constructed by ZnO and InO chains with isonicotinate and terephthalate linkers, respectively. The bond length of Zn-O is shorter than that of In-O, so the elastic moduli of QMOF-1 are larger than those of QMOF-2. In addition, the density of QMOF-1 (1.52 g/cm^{3}) is less than that of QMOF-2 (1.76 g/cm^{3}). Both aspects result in 53-65%, 101-111%, and 39-55% higher , , and of QMOF-1 than those of QMOF-2 (Figs. S8, S9). Furthermore, the anisotropies in acoustic velocities of both MOFs are small, with the highest value less than 1.8. The , , and values of QMOF-2 are approximately 47-58%, 48-57%, and 57-77% of those of quartz. However, the , , and values are approximately 15%, 2%, and 8% higher than those of quartz, although the , , and values are approximately 4%, 28%, and 20% lower.

As shown in Figures 2 and S10-17, the values of all formate frameworks are generally larger than those of porous MOFs and are broadly similar to that of the most structurally related BaTiO_{3}. Taking [DMA][Mg(HCOO)_{3}] as an example, its , , and are approximately 101-110%, 72-92%, and 93-146% of those of BaTiO_{3}. In addition, these formates generally have small anisotropy indices below 1.8 apart from [MA][Mn(HCOO)_{3}] (, ). The organic-inorganic framework DABCOH_{2}K(ClO_{4})_{3} has a density of 3.15 g/cm^{3}, which is much higher than those of all aforementioned MOFs. Considering the similar elastic moduli of DABCOH_{2}K(ClO_{4})_{3} and perovskite formates, its acoustic velocities are generally smaller (Fig. S18). Compared with BaTiO_{3}, its , , and are approximately 57-83%, 48-60%, and 52% lower, although its is almost the same.

##### 2.3. Acoustic Impedance

Acoustic impedance () describes the ratio of acoustic pressure to sound flow propagating in a material, which can be defined as the product of the acoustic velocity and the density of the material: , where and are the acoustic velocity and density, respectively. The acoustic impedance can also be divided into transverse acoustic impedance and longitudinal acoustic impedance according to the different acoustic velocities. The longitudinal acoustic impedance is widely used in practice, so we calculated the maximum and minimum longitudinal acoustic impedances of these 22 materials.

As seen from Figure 5, the longitudinal acoustic impedances of MOF-5, ZIF-8, and MAF-7 have very similar values and are the smallest among all 22 materials. The values of ZIF-8 (, MPa∙s∙m^{-1}) and MAF-7 (, MPa∙s∙m^{-1}) are only approximately 30% of those of their inorganic counterpart, sodalite (, MPa∙s∙m^{-1}), which is attributed to their much lower elastic moduli and density. For the MIL family, their density and increase in the following sequence: , which leads to the same trend in their acoustic impedances: MIL-47 ( MPa∙s∙m^{-1}) < MIL-53(Al)lp ( MPa∙s∙m^{-1}) < MIL-53(Ga)lp ( MPa∙s∙m^{-1}). Both QMOF-1 and QMOF-2 have quartz-like structures, so the study of their acoustic impedances is of great significance for their application in piezoelectric sensors. The and of QMOF-1 and QMOF-2 are approximately 46% and 68% and approximately 38% and 51% of those of quartz, respectively. For the metal formates, their and values are in the range of 3.72-5.26 and 4.29-9.17 MPa∙s∙m^{-1}. Taking the DMAMg(HCOO)_{3} perovskite as an example, its and are approximately 45% and 40% of those of the perovskite oxide BaTiO_{3}. By comparison, it can be seen that the acoustic impedances of MOFs are generally lower than those of traditional inorganic materials.

##### 2.4. Piezoelectric Properties

The acoustic and piezoelectric properties of materials are closely related since acoustic waves can be detected by piezoelectric materials or vice versa. According to the type of energy conversion, the piezoelectric effects can be divided into direct and indirect effects. The direct piezoelectric effect occurs when a piezoelectric material generates a dipole moment under external stress, while the indirect piezoelectric effect refers to the strain generated in piezoelectric materials when they are placed in an electric field [35].

Among the 19 crystals, QMOF-1, QMOF-2, and DMAMg(HCOO)_{3} have piezoelectric properties owing to their 3, 622, and *m* noncentrosymmetric point groups, respectively. The voltage signal harvested by our experiment also confirmed the piezoelectric property of QMOF-1 (Fig. S22). Because of the structural similarity, the origin of the piezoelectric effect of QMOF-1 is reminiscent of that of quartz. As shown in Figure 6, Zn-O polyhedra and isonicotinate linkers are considered negatively and positively charged components, respectively. In the absence of stress, the positive and negative charge centres overlap, which leads to no macroscopic spontaneous polarization. However, negative and positive charges are generated on the top and bottom surfaces of the crystal, respectively, when the *a*-axis is stressed.

To further analyse the piezoelectric properties of QMOF-1, QMOF-2, and DMAMg(HCOO)_{3}, the density functional perturbation theory (DFPT) calculation method is used to obtain their piezoelectric tensors. It is well known that the piezoelectric stress tensor , elastic compliance matrix , and piezoelectric strain tensor can be related *via* the following formula: [36, 37]. The full piezoelectric stress and strain tensors of all three MOFs are listed in Tables S2, S4, and S5. In all 3 MOFs, the ionic contribution is much larger than the electronic contribution to the piezoelectric stress tensors. For QMOF-1, the piezoelectric strain tensor (23.64 pC/N) is 31.3 times larger than that of quartz ( pC/N), and its is 4.60 pC/N. QMOF-2 only has a tensor (2.52 pC/N) due to its 622 symmetry, which is approximately 3.5 times higher than that of quartz. The and values of DMAMg(HCOO)_{3} are 3.19 and -10.5 pC/N, respectively, which are only approximately 2% and 4% of those of the commercial piezoelectric perovskite oxide BaTiO_{3}. In addition, the value of DMAMg(HCOO)_{3} is 2-3 orders of magnitude lower than those of high-performance hybrid perovskite piezoelectrics, such as [TMCM][MnCl_{3}] ( ammonium, pC/N) and (TMFM)_{x}(TMCM)_{1–x}CdCl_{3} ( ammonium, pC/N) [38, 39].

The piezoelectric strain tensor in these 3 MOFs is mainly affected by two factors: the magnitude of the atomic displacement caused by the strain and the charge change caused by the atomic displacement. The Born effective charge (BEC) is defined as the change in electron polarization caused by ion displacement and is the second derivative of the energy with respect to the displacement and electric field. In this regard, the BEC () can directly reflect the contribution of the second factor to the piezoelectric strain tensor. The BEC can be expressed as: where , , and represent the atomic displacement, electric field, and strain, respectively. denotes the total energy of the material in the ground state, and is the unit cell volume of the material.

To explore the direction-dependent piezoelectric effects, the BEC of the zinc ions along different directions in QMOF-1 was calculated, and the results are listed in Table S3. There are three nonequivalent zinc ions in the unit cell of QMOF-1, and the sum of their BECs along the -axis is approximately half of those along the /-axis. As the zinc ion normally shows a positive divalence, it is prone to exhibit charge transfer if its BEC is larger than 2. The BECs of the three zinc elements are 1.59, 1.59, and 1.58 along the -axis, while they are 2.11, 3.79, and 2.33 along the -axis and 3.38, 1.69, and 3.16 along the -axis. Such a significant difference indicates that transfer of their charges is more difficult along the -axis than along the /-axis, hence giving rise to smaller but larger and .

#### 3. Discussion

Our above analyses demonstrate the tight structural-acoustic property relationships of MOFs and suggest some tuning approaches *via* crystal engineering. First, variation of the organic ligands can result in prominent differences in the elastic moduli of MOFs [14] and the corresponding acoustic velocities and impedance for the same structural topology and similar densities. As this strategy has been well demonstrated by the difference in acoustic properties between the prototypical ZIF-8 and MAF-7 by simply varying the ligand, it could be widely applied in other MOF systems. Second, metal nodes play an important role in determining the framework rigidity and therefore substantially influence the corresponding acoustic properties. For the MIL-53 family, the Al analogue is stiffer than the Ga counterpart due to the shorter coordination bonds [32], which leads to its higher sound velocities but lower acoustic anisotropy. Similarly, the Jahn-Teller effect and ligand splitting stabilization energy could be employed to modulate the elasticity and corresponding acoustic properties of MOFs [40]. Furthermore, the coordination number of metal nodes would be an alternative way to tune the acoustic properties of MOFs. For instance, MOFs with highly coordinated metal nodes show higher resistance towards shear stress than those with low-coordinated structures and therefore exhibit higher transverse wave velocities. Third, the acoustic properties of MOFs could also be tuned and engineered by varying host-guest interactions. For instance, the much stronger hydrogen bonding in [GUA][Mn(HCOO)_{3}] leads to twice the Young’s moduli compared to those of the analogous [AZE][Mn(HCOO)_{3}] (), which naturally results in its much larger acoustic velocities and impedance [41]. On the other hand, gas molecule-filled MOFs exhibit lower acoustic velocities owing to unaltered elastic moduli but significantly increased density compared to their parent phases [42].

In terms of the piezoelectricity of MOFs, very few studies have been reported, and more research efforts need to be devoted to this field, especially considering the availability of vast numbers of noncentrosymmetric MOFs with diverse structures and topologies. When sufficient piezoelectric MOFs are identified and their structure-property relationship is corroborated, the next goal would be to optimize their performance to catch up with those of their conventional ceramic counterparts. In this regard, screening more noncentrosymmetric MOFs using DFPT calculations with the aid of artificial intelligence would be a promising method if combined with complementary experimental confirmation.

#### 4. Conclusions

In conclusion, the acoustic properties of 19 MOFs were systematically investigated *via* density functional theory (DFT) calculations and by solving the Christoffel equations. Our results show that the MOFs have a very wide range of acoustic velocities due to their structural diversity and abundant chemical compositions. Meanwhile, some MOFs can exhibit very high anisotropy in acoustic velocities, which could be applicable in acoustic switches. In addition, the low densities and moduli of MOFs are responsible for their small acoustic impedances compared with their inorganic counterparts. Furthermore, the piezoelectric properties, which are closely related to the acoustic properties, were calculated for two quartz-like and one perovskite-like MOFs. Our results reveal that these MOFs can exhibit significant piezoelectricity, which mainly originates from the ionic contribution. Finally, we propose some strategies for engineering the acoustic properties of MOFs and discuss future directions for studying the piezoelectricity of MOFs. This work offers an original study of the fundamental acoustic properties of MOFs, which we believe could inspire further attention in this largely unexplored area.

#### 5. Methods

Structural optimizations and calculations of elastic constants based on DFT were performed using a plane-wave basis set as implemented in the Vienna ab initio simulation package (VASP) [43–45]. Projector augmented wave (PAW) pseudopotentials were employed to describe ion-electron interactions [46]. The considered PAW pseudopotentials and numbers of valence electrons were H (1), C (4), N (5), O (6), Na_pv (7), Mg (2), Al (3), Si (4), Cl (7), K_sv (9), Ca_sv (10), Ti_sv (12), V_sv (13), Mn_pv (13), Fe (8), Co (9), Cu (11), Zn (12), Ga_d (13), and Ba_sv (10). The GGA (PBEsol) functional was used to further accurately describe the interactions within each MOF structure. Considering the abundant van der Waals interactions existed in MOF systems, dispersion corrected calculations are known to be crucial in the correct optimization of their crystal structures. Specifically, GGA (PBEsol) underestimates the binding energy, hence resulting in larger lattice parameters for these molecular systems. By applying the dispersion corrections, more reasonable binding energies could be obtained, which gives rise to more reliable lattice parameters [47, 48]. As elastic constants are obtained by the second derivative of the total energy with respective to specific strain and strongly correlated to the relative change of the crystal structure, obtaining accurate cell parameters are of vital importance in the calculations. In this regard, the van der Waals corrected zero damping DFT-D3 method was chosen in our calculations. Considering that the calculation of elastic constants requires a high standard for the structures of materials, the atomic positions and cell parameters were completely relaxed when the structure was optimized. The elastic constants were calculated by setting six finite distortions of the lattice and the step size to 0.015 Å. The total energy converged to within 10^{-7} eV, and the residual forces on each atom were less than 0.005 eV/Å. The energy cutoff for the plane-wave basis and the smallest allowed spacing between -points of these compounds were set to 500 eV and 0.5 Å^{-1}, respectively.

The acoustic velocity was obtained by solving the Christoffel equation with the material density and the elastic constants calculated by DFT [49]. The Christoffel equation can be expressed as: where is the elastic displacement field and , , , and are the Christoffel matrix, acoustic velocity in the crystal, density, and direction vector, respectively.

The second derivative of the total energy reflects many properties of the calculated structure, such as the Born effective charge (BEC) and piezoelectric and dielectric tensors. The density functional perturbation theory proposed by Baroni et al. in 1987 can effectively calculate the derivative of the total energy with respect to strain, atomic displacement, and electric field [50–53]. The expression for the total energy (, , ) under perturbation is as follows: where , , and represent the atomic displacement, electric field, and strain, respectively. (, , ) represents the total energy of the material in the ground state, is the unit cell volume of the material and is the unit cell volume of the material under perturbation, and is the electric field polarization.

The piezoelectric effect reflects the deformation degree of a material under an electric field, which can be evaluated by the piezoelectric stress tensor (). The piezoelectric tensors calculated by DFPT are divided into the electronic contribution () and ionic contribution (). The former is an ion-clamped piezoelectric tensor that ignores the effects of atomic relaxation; in contrast, the latter takes into account the effect of the displacement of atoms under strain on the piezoelectric tensor. Therefore, the piezoelectric tensor can be expressed as:

In this part of the DFPT calculation, we used the previously highly optimized crystal structure as the calculation model and improved the self-consistent field convergence standard to 10^{-8} eV.

#### Data Availability

All data needed in the paper are present in the paper and in the supplementary section. Additional data related to this paper may be requested from the authors.

#### Conflicts of Interest

The authors declare no conflict of interest.

#### Authors’ Contributions

Z.-G.L. and K.L. contributed equally to this work. Z.-G.L., K.L., W.L., and X.-H.B. designed and organized the research. Theoretical calculations were performed by Z.-G.L. and L.-Y.D. Data analysis was performed by Z.-G.L., K.L., and W.L. Z.-G.L., K.L., and W.L. wrote the full manuscript, and all authors discussed the results and commented on the manuscript.

#### Acknowledgments

We gratefully acknowledge the financial support from the National Natural Science Foundation of China (Nos. 91856124, 22035003, 21975132, and 21991143) and the Fundamental Research Funds for the Central Universities (No. 63196006).

#### Supplementary Materials

Table S1: summary of the density (g/cm^{3}) and elastic constants (GPa) of the MOFs and inorganic materials obtained from the DFT calculations. Table S2: summary of the transverse wave velocities and longitudinal wave velocities (m/s) of the MOFs and inorganic materials. Table S3: clamped-ion (C/m^{2}), relaxed-ion (C/m^{2}), and piezoelectric strain tensor (pC/N) for QMOF-1. Table S4: summary of the BEC of QMOF-1 obtained from the DFPT calculation. Table S5: clamped-ion (C/m^{2}), relaxed-ion (C/m^{2}), and piezoelectric strain tensor (pC/N) for QMOF-2. Table S6: clamped-ion (C/m^{2}), relaxed-ion ê (C/m^{2}), and piezoelectric strain tensor (pC/N) for DMAMg(HCOO)_{3}. Figure S1: (a) 3D surfaces of the acoustic velocity for MOF-5. (b–d) 2D plots of the acoustic velocity for MOF-5 projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S2: (a) 3D surfaces of the acoustic velocity for ZIF-8. (b–d) 2D plots of the acoustic velocity for ZIF-8 projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S3: (a) 3D surfaces of the acoustic velocity for MAF-7. (b–d) 2D plots of the acoustic velocity for MAF-7 projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S4: (a) 3D surfaces of the acoustic velocity for DMOF-1loz. (b–d) 2D plots of the acoustic velocity for DMOF-1loz projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S5: (a) 3D surfaces of the acoustic velocity for DMOF-1sq. (b–d) 2D plots of the acoustic velocity for DMOF-1sq projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S6: (a) 3D surfaces of the acoustic velocity for MIL-47. (b–d) 2D plots of the acoustic velocity for MIL-47 projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S7: (a) 3D surfaces of the acoustic velocity for MIL-53(Al)lp. (b–d) 2D plots of the acoustic velocity for MIL-53(Al)lp projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S8: (a) 3D surfaces of the acoustic velocity for QMOF-1. (b–d) 2D plots of the acoustic velocity for QMOF-1 projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S9: (a) 3D surfaces of the acoustic velocity for QMOF-2. (b–d) 2D plots of the acoustic velocity for QMOF-2 projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S10: (a) 3D surfaces of the acoustic velocity for NaMn(HCOO)_{3}. (b–d) 2D plots of the acoustic velocity for NaMn(HCOO)_{3} projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S11: (a) 3D surfaces of the acoustic velocity for NH_{4}Co(HCOO)_{3}. (b–d) 2D plots of the acoustic velocity for NH_{4}Co(HCOO)_{3} projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S12: (a) 3D surfaces of the acoustic velocity for NH_{4}Zn(HCOO)_{3}. (b–d) 2D plots of the acoustic velocity for NH_{4}Zn(HCOO)_{3} projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S13: (a) 3D surfaces of the acoustic velocity for GuaZn(HCOO)_{3}. (b–d) 2D plots of the acoustic velocity for GuaZn(HCOO)_{3} projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S14: (a) 3D surfaces of the acoustic velocity for HyzZn(HCOO)_{3}. (b–d) 2D plots of the acoustic velocity for HyzZn(HCOO)_{3} projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S15: (a) 3D surfaces of the acoustic velocity for MAMn(HCOO)_{3}. (b–d) 2D plots of the acoustic velocity for MAMn(HCOO)_{3}projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S16: (a) 3D surfaces of the acoustic velocity for MAZn(HCOO)_{3}. (b–d) 2D plots of the acoustic velocity for MAZn(HCOO)_{3} projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S17: (a) 3D surfaces of the acoustic velocity for DMAMg(HCOO)_{3}. (b–d) 2D plots of the acoustic velocity for DMAMg(HCOO)_{3} projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S18: (a) 3D surfaces of the acoustic velocity for DABCOH_{2}K(ClO_{4})_{3}. (b–d) 2D plots of the acoustic velocity for DABCOH_{2}K(ClO_{4})_{3} projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S19: (a) 3D surfaces of the acoustic velocity for quartz. (b–d) 2D plots of the acoustic velocity for quartz projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S20: (a) 3D surfaces of the acoustic velocity for BaTiO_{3}. (b–d) 2D plots of the acoustic velocity for BaTiO_{3} projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S21: (a) 3D surfaces of the acoustic velocity for Na_{4}Al_{3}(SiO_{4})_{3}. (b–d) 2D plots of the acoustic velocity for Na_{4}Al_{3}(SiO_{4})_{3} projected normal to the (100), (010), and (001) planes, respectively. The green, red, and blue lines represent , , and , respectively. Figure S22: (a) Schematic diagram of piezoelectric energy harvester based on polycrystalline pressed samples of QMOF-1. (b) Generated output voltages for QMOF-1 under certain mechanical press recorded by digital oscilloscope. The experiments were conducted based on a sandwich-structured Cu electrode-polycrystalline samples-Cu electrode device, under a periodic impact created by a mechanical contactor. The output voltage was measured using an oscilloscope (Keysight DSOX3024T).* (Supplementary Materials)*

#### References

- A. K. Cheetham, C. N. R. Rao, and R. K. Feller, “Structural diversity and chemical trends in hybrid inorganic–organic framework materials,”
*Chemical Communications*, vol. 38, no. 46, pp. 4780–4795, 2006. View at: Publisher Site | Google Scholar - C. N. R. Rao, A. K. Cheetham, and A. Thirumurugan, “Hybrid inorganic–organic materials: a new family in condensed matter physics,”
*Journal of Physics-Condensed Matter*, vol. 20, no. 8, article 083202, 2008. View at: Publisher Site | Google Scholar - A. K. Cheetham and C. N. R. Rao, “There’s room in the middle,”
*Science*, vol. 318, no. 5847, pp. 58-59, 2007. View at: Publisher Site | Google Scholar - Z. C. Hu, B. J. Deibert, and J. Li, “Luminescent metal-organic frameworks for chemical sensing and explosive detection,”
*Chemical Society Reviews*, vol. 43, no. 16, pp. 5815–5840, 2014. View at: Publisher Site | Google Scholar - E. Coronado, “Molecular magnetism: from chemical design to spin control in molecules, materials and devices,”
*Nature Reviews Materials*, vol. 5, no. 2, pp. 87–104, 2020. View at: Publisher Site | Google Scholar - P. Jain, A. Stroppa, D. Nabok et al., “Switchable electric polarization and ferroelectric domains in a metal- organic-framework,”
*Npj Quantum Materials*, vol. 1, no. 1, article 16012, 2016. View at: Publisher Site | Google Scholar - A. Stroppa, P. Barone, P. Jain, J. M. Perez-Mato, and S. Picozzi, “Hybrid improper ferroelectricity in a multiferroic and magnetoelectric metal-organic framework,”
*Advanced Materials*, vol. 25, no. 16, pp. 2284–2290, 2013. View at: Publisher Site | Google Scholar - B. Li, H. M. Wen, Y. J. Cui, W. Zhou, G. Qian, and B. Chen, “Emerging multifunctional metal-organic framework materials,”
*Advanced Materials*, vol. 28, no. 40, pp. 8819–8860, 2016. View at: Publisher Site | Google Scholar - Q. R. S. Miller, S. K. Nune, H. T. Schaef et al., “Microporous and flexible framework acoustic metamaterials for sound attenuation and contrast agent applications,”
*ACS Applied Materials Interfaces*, vol. 10, no. 51, pp. 44226–44230, 2018. View at: Publisher Site | Google Scholar - H. Ahmed, A. R. Rezk, J. J. Richardson et al., “Acoustomicrofluidic assembly of oriented and simultaneously activated metal- organic frameworks,”
*Nature Communications*, vol. 10, no. 1, article 2282, 2019. View at: Publisher Site | Google Scholar - K. Leong, M. E. Foster, B. M. Wong et al., “Energy and charge transfer by donor-acceptor pairs confined in a metal-organic framework: a spectroscopic and computational investigation,”
*Journal of Materials Chemistry A*, vol. 2, no. 10, pp. 3389–3398, 2014. View at: Publisher Site | Google Scholar - S. Haussühl,
*Physical properties of crystals: an introduction*, WILEY-VCH Verlag GmbH & Co. KgaA, 2007. - J. C. Tan and A. K. Cheetham, “Mechanical properties of hybrid inorganic–organic framework materials: establishing fundamental structure–property relationships,”
*Chemical Society Reviews*, vol. 40, no. 2, pp. 1059–1080, 2011. View at: Publisher Site | Google Scholar - W. Li, S. Henke, and A. K. Cheetham, “Research update: mechanical properties of metal-organic frameworks – influence of structure and chemical bonding,”
*APL Materials*, vol. 2, no. 12, article 123902, 2014. View at: Publisher Site | Google Scholar - J. Heinen, N. C. Burtch, K. S. Walton, and D. Dubbeldam, “Flexible force field parameterization through fitting on the ab initio-derived elastic tensor,”
*Journal of Chemical Theory and Computation*, vol. 13, no. 8, pp. 3722–3730, 2017. View at: Publisher Site | Google Scholar - W. Xu, B. Tu, Q. Liu et al., “Anisotropic reticular chemistry,”
*Nature Reviews Materials*, vol. 5, no. 10, pp. 764–779, 2020. View at: Publisher Site | Google Scholar - J. L. C. Rowsell, A. R. Millward, K. S. Park, and O. M. Yaghi, “Hydrogen sorption in functionalized metal−organic frameworks,”
*Journal of the American Chemical Society*, vol. 126, no. 18, pp. 5666-5667, 2004. View at: Publisher Site | Google Scholar - K. S. Park, Z. Ni, A. P. Cote et al., “Exceptional chemical and thermal stability of zeolitic imidazolate frameworks,”
*Proceedings of the National Academy of Sciences of the United States of America*, vol. 103, no. 27, pp. 10186–10191, 2006. View at: Publisher Site | Google Scholar - J. P. Zhang, A. X. Zhu, R. B. Lin, X. L. Qi, and X. M. Chen, “Pore surface tailored SOD-type metal-organic zeolites,”
*Advanced Materials*, vol. 23, no. 10, pp. 1268–1271, 2011. View at: Publisher Site | Google Scholar - D. N. Dybtsev, H. Chun, and K. Kim, “Rigid and flexible: a highly porous metal–organic framework with unusual guest-dependent dynamic behavior,”
*Angewandte Chemie-International Edition*, vol. 43, no. 38, pp. 5033–5036, 2004. View at: Publisher Site | Google Scholar - P. G. Yot, Q. Ma, J. Haines et al., “Large breathing of the MOF MIL-47(VIV) under mechanical pressure: a joint experimental–modelling exploration,”
*Chemical Science*, vol. 3, no. 4, pp. 1100–1104, 2012. View at: Publisher Site | Google Scholar - Y. Liu, J. H. Her, A. Dailly, A. J. Ramirez-Cuesta, D. A. Neumann, and C. M. Brown, “Reversible structural transition in MIL-53 with large temperature hysteresis,”
*Journal of the American Chemical Society*, vol. 130, no. 35, pp. 11813–11818, 2008. View at: Publisher Site | Google Scholar - J. Sun, L. Weng, Y. Zhou et al., “QMOF-1 and QMOF-2: three-dimensional metal–organic open frameworks with a quartzlike topology,”
*Angewandte Chemie-International Edition*, vol. 41, no. 23, pp. 4471–4473, 2002. View at: Publisher Site | Google Scholar - E. Eikeland, N. Lock, M. Filsø et al., “Alkali metal ion templated transition metal formate framework materials: synthesis, crystal structures, ion migration, and magnetism,”
*Inorganic Chemistry*, vol. 53, no. 19, pp. 10178–10188, 2014. View at: Publisher Site | Google Scholar - G. C. Xu, X. M. Ma, L. Zhang, Z. M. Wang, and S. Gao, “Disorder−order ferroelectric transition in the metal formate framework of [NH
_{4}][Zn(HCOO)_{3}],”*Journal of the American Chemical Society*, vol. 132, no. 28, pp. 9588–9590, 2010. View at: Publisher Site | Google Scholar - G. C. Xu, W. Zhang, X. M. Ma et al., “Coexistence of magnetic and electric orderings in the metal–formate frameworks of [NH
_{4}][M(HCOO)_{3}],”*Journal of the American Chemical Society*, vol. 133, no. 38, pp. 14948–14951, 2011. View at: Publisher Site | Google Scholar - Z. Wang, B. Zhang, T. Otsuka, K. Inoue, H. Kobayashi, and M. Kurmoo, “Anionic NaCl-type frameworks of [Mn
^{II}(HCOO)^{3−}], templated by alkylammonium, exhibit weak ferromagnetism,”*Dalton Transactions*, vol. 4, no. 15, pp. 2209–2216, 2004. View at: Publisher Site | Google Scholar - K. L. Hu, M. Kurmoo, Z. M. Wang, and S. Gao, “Metal-organic perovskites: synthesis, structures, and magnetic properties of [C(NH
_{2})_{3}][MII(HCOO)_{3}] (M=Mn, Fe, Co, Ni, Cu, and Zn; C(NH2)3= Guanidinium),”*Chemistry - A European Journal*, vol. 15, no. 44, pp. 12050–12064, 2009. View at: Publisher Site | Google Scholar - Z. Wang, X. Zhang, S. R. Batten, M. Kurmoo, and S. Gao, “[CH
_{3}NH_{2}(CH_{2})_{2}NH_{2}CH_{3}][M_{2}(HCOO)_{6}] (M = Mn^{II}and Co^{II}): weak ferromagnetic metal formate frameworks of unique binodal 6-connected (4^{12}·6^{3})(4^{9}·6^{6}) topology, templated by a diammonium cation,”*Inorganic Chemistry*, vol. 46, no. 21, pp. 8439–8441, 2007. View at: Publisher Site | Google Scholar - S. Chen, R. Shang, K. L. Hu, Z. M. Wang, and S. Gao, “[NH
_{2}NH_{3}][M(HCOO)_{3}] (M = Mn^{2+}, Zn^{2+}, Co^{2+}and Mg^{2+}): structural phase transitions, prominent dielectric anomalies and negative thermal expansion, and magnetic ordering,”*Inorganic Chemistry Frontiers*, vol. 1, no. 1, pp. 83–98, 2014. View at: Publisher Site | Google Scholar - Z. M. Jin, Y. J. Pan, X. F. Li, M. L. Hu, and L. Shen, “Diazabicyclo[2.2.2]octane-1,4-diium occluded in cubic anionic coordinated framework: the role of trifurcated hydrogen bonds of N-H⋯O and C-H⋯O,”
*Journal of Molecular Structure*, vol. 660, no. 1-3, pp. 67–72, 2003. View at: Publisher Site | Google Scholar - A. U. Ortiz, A. Boutin, A. H. Fuchs, and F. X. Coudert, “Anisotropic elastic properties of flexible metal-organic frameworks: how soft are soft porous crystals?”
*Physical Review Letters*, vol. 109, no. 19, article 195502, 2012. View at: Publisher Site | Google Scholar - J. C. Tan, B. Civalleri, C. C. Lin et al., “Exceptionally low shear modulus in a prototypical imidazole-based metal-organic framework,”
*Physical Review Letters*, vol. 108, no. 9, article 095502, 2012. View at: Publisher Site | Google Scholar - D. F. Bahr, J. A. Reid, W. M. Mook et al., “Mechanical properties of cubic zinc carboxylate IRMOF-1 metal-organic framework crystals,”
*Physical Review B*, vol. 76, no. 18, article 184106, 2007. View at: Publisher Site | Google Scholar - W. G. Nelson,
*Piezoelectric materials: structure, properties and applications*, Nova Science Pub Inc, 2010. - Z. X. Wang, H. Zhang, F. Wang et al., “Superior transverse piezoelectricity in a halide perovskite molecular ferroelectric thin film,”
*Journal of the American Chemical Society*, vol. 142, no. 29, pp. 12857–12864, 2020. View at: Publisher Site | Google Scholar - H. Wang, H. H. Liu, Z. Y. Zhang et al., “Large piezoelectric response in a family of metal-free perovskite ferroelectric compounds from first-principles calculations,”
*Npj Computational Materials*, vol. 5, no. 1, article 17, 2019. View at: Publisher Site | Google Scholar - Y. M. You, W. Q. Liao, D. W. Zhao et al., “An organic-inorganic perovskite ferroelectric with large piezoelectric response,”
*Science*, vol. 357, no. 6348, pp. 306–309, 2017. View at: Publisher Site | Google Scholar - W. Q. Liao, D. W. Zhao, Y. Y. Tang et al., “A molecular perovskite solid solution with piezoelectricity stronger than lead zirconate titanate,”
*Science*, vol. 363, no. 6432, pp. 1206–1210, 2019. View at: Publisher Site | Google Scholar - L. J. Ji, S. J. Sun, Y. Qin, K. Li, and W. Li, “Mechanical properties of hybrid organic-inorganic perovskites,”
*Coordination Chemistry Reviews*, vol. 391, pp. 15–29, 2019. View at: Publisher Site | Google Scholar - W. Li, A. Thirumurugan, P. T. Barton et al., “Mechanical tunability via hydrogen bonding in metal-organic frameworks with the perovskite architecture,”
*Journal of the American Chemical Society*, vol. 136, no. 22, pp. 7801–7804, 2014. View at: Publisher Site | Google Scholar - A. Ba, A. Kovalenko, C. Aristégui, O. Mondain-Monval, and T. Brunet, “Soft porous silicone rubbers with ultra-low sound speeds in acoustic metamaterials,”
*Scientific Reports*, vol. 7, no. 1, article 40106, 2017. View at: Publisher Site | Google Scholar - G. Kresse and J. Furthmüller, “Efficient iterative schemes forab initiototal-energy calculations using a plane-wave basis set,”
*Physical Review B*, vol. 54, no. 16, pp. 11169–11186, 1996. View at: Publisher Site | Google Scholar - G. Kresse and J. Furthmüller, “Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set,”
*Computational Materials Science*, vol. 6, no. 1, pp. 15–50, 1996. View at: Publisher Site | Google Scholar - G. Kresse and J. Hafner, “Ab initiomolecular dynamics for liquid metals,”
*Physical Review B*, vol. 47, no. 1, pp. 558–561, 1993. View at: Publisher Site | Google Scholar - P. E. Blochl, “Projector augmented-wave method,”
*Physical Review B*, vol. 50, no. 24, pp. 17953–17979, 1994. View at: Publisher Site | Google Scholar - S. Grimme, J. Antony, S. Ehrlich, and H. Krieg, “A consistent and accurateab initioparametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu,”
*Journal of Chemical Physics*, vol. 132, no. 15, article 154104, 2010. View at: Publisher Site | Google Scholar - B. Chan and G. E. Ball, “A benchmark ab initio and DFT study of the structure and binding of methane in the
*σ*–alkane complex_{CpRe(CO)2(CH4)},”*Journal of Chemical Theory and Computation*, vol. 9, no. 5, pp. 2199–2208, 2013. View at: Publisher Site | Google Scholar - B. A. Auld,
*Acoustic Fields and Waves in Solid*, John Wiley, 1973. - D. R. Hamann, X. Wu, K. M. Rabe, and D. Vanderbilt, “Erratum: Metric tensor formulation of strain in density-functional perturbation theory,”
*Physical Review B*, vol. 72, no. 7, article 079901, 2005. View at: Publisher Site | Google Scholar - A. P. Sakhya, A. Dutta, S. Shannigrahi, and T. P. Sinha, “Electronic structure, optical dielectric constant and born effective charge of EuAlO
_{3},”*Journal of Physics and Chemistry of Solids*, vol. 88, pp. 1–7, 2016. View at: Publisher Site | Google Scholar - A. Filippetti and N. A. Spaldin, “Strong-correlation effects in Born effective charges,”
*Physical Review B*, vol. 68, no. 4, article 045111, 2003. View at: Publisher Site | Google Scholar - S. Baroni, “Green’s-function approach to linear response in solids,”
*Physical Review Letters*, vol. 58, no. 18, pp. 1861–1864, 1987. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2021 Zhi-Gang Li et al. Exclusive Licensee Science and Technology Review Publishing House. Distributed under a Creative Commons Attribution License (CC BY 4.0).