Piezoelectric energy harvesting is a process in which energy in the form of kinetic movements can be harvested and converted into useful electrical energy using piezoelectric materials. Metal-organic frameworks (MOFs) have a huge potential for piezoelectric energy harvesting owing to their high flexibility, structural tunability, and very low dielectric constants due to their high porosity. The piezoelectric constant d relevant for piezoelectric energy harvesting depends on the piezoelectric constant e and the flexibility of the structure (i.e. mechanical properties). The mechanical properties of MOFs have previously been extensively studied but the piezoelectric constant e was never explored for MOFs. In this work, we generate a database of piezoelectric properties, specifically e for around ∼1608 previously synthesized non-centrosymmetric MOF structures. The calculations were performed using the density functional perturbation theory (DFPT) method. The highest piezoelectric constant e obtained in this work is approximately ∼2.76 C m⁻², which is significantly higher than that of the flexible organic piezoelectric polymer polyvinylidene fluoride (PVDF) and its copolymers. In this work, we analyze and identify structural factors that influence the values of the piezoelectric constant for high-performing MOFs. Based on that, a series of guidelines for the design of MOF structures that can lead to a high piezoelectric constant e are presented. One class of high-performing piezoelectric MOFs is based on polar patterns of O—(short)—Mo—(long)—O unequal bond length, reminiscent of ferroelectric inorganic oxides. This class could have potential for ferroelectricity, meaning that the bond length pattern could be reversed by external electrical field. We substantiate this by showing experimentally via SHG-microscopy that the O—(short)—Mo—(long)—O unequal bond lengths are indeed malleable by external conditions.

By tuning the steric environment and free pore space in metal-organic frameworks, a large variety of rotor dynamics of the organic linkers can appear. Nitrofunctionalized MIL-53 is a terephthalate-linker-based MOF that shows coupled rotor dynamics between the neighboring linkers along the pore direction. Here, we use classical molecular dynamics up to 6 × 2 × 2 supercells to investigate the range of the correlated linker dynamics. Interestingly, we observe an PNPNPNPN... conformational arrangement (P = nearly planar and N = nonplanar) for the conformations of the linkers in a row along the pore direction in the MOF. We identified correlated linker dynamics emerging among the direct and next nearest neighboring linkers along the pore. Due to 180° rotational flips of the planar linkers along the pore, (1) a change in the width of librations in their direct neighbors (PN) is observed; (2) intriguingly, their next nearest planar neighbors (PP) rotate between 0° and ±180° to reattain aligned (0°, 0°) or (±180°, ±180°) conformations. The presence of correlated dynamics in such linkers over long-length scales occurring at nanoseconds time scales is desirable for applications like ferroelectric switching or diffusion control via geared linker rotation, and this work provides insight into the design for such applications.

Metal-organic frameworks (MOFs) are a class of nanoporous crystalline materials with very high structural tunability. They possess a very low dielectric permittivity ϵr due to their porosity and hence are favorable for piezoelectric energy harvesting. Even though they have huge potential as piezoelectric materials, a detailed analysis and structure-property relationship of the piezoelectric properties in MOFs are lacking so far. This work focuses on a class of cubic non-centrosymmetric MOFs, namely, zeolitic imidazolate frameworks (ZIFs) to rationalize how the variation of different building blocks of the structure, that is, metal node and linker substituents affect the piezoelectric constants. The piezoelectric tensor for the ZIFs is computed from ab initio theoretical methods. From the calculations, we analyze the different contributions to the final piezoelectric constant d14, namely, the clamped ion (e140) and the internal strain (e14int) contributions and the mechanical properties. For the studied ZIFs, even though e14 (e140 + e14int) is similar for all ZIFs, the resultant piezoelectric coefficient d14 calculated from piezoelectric constant e14 and elastic compliance constant s44 varies significantly among the different structures. It is the largest for CdIF-1 (Cd2+ and -CH3 linker substituent). This is mainly due to the higher elasticity or flexibility of the framework. Interestingly, the magnitude of d14 for CdIF-1 is higher than II-VI inorganic piezoelectrics and of a similar magnitude as the quintessential piezoelectric polymer polyvinylidene fluoride.

The organic components in metal-organic frameworks (MOFs) are unique: they are embedded in a crystalline lattice, yet, as they are separated from each other by tunable free space, a large variety of dynamic behavior can emerge. These rotational dynamics of the organic linkers are especially important due to their influence over properties such as gas adsorption and kinetics of guest release. To fully exploit linker rotation, such as in the form of molecular machines, it is necessary to engineer correlated linker dynamics to achieve their cooperative functional motion. Here, we show that for MIL-53, a topology with closely spaced rotors, the phenylene functionalization allows researchers to tune the rotors' steric environment, shifting linker rotation from completely static to rapid motions at frequencies above 100 MHz. For steric interactions that start to inhibit independent rotor motion, we identify for the first time the emergence of coupled rotation modes in linker dynamics. These findings pave the way for function-specific engineering of gear-like cooperative motion in MOFs.