Thermal energy management in metal-organic frameworks (MOFs) is an important, yet often neglected, challenge for many adsorption-based applications such as gas storage and separations. Despite its importance, there is insufficient understanding of the structure-property relationships governing thermal transport in MOFs. To provide a data-driven perspective into these relationships, here we perform large-scale computational screening of thermal conductivity k in MOFs, leveraging classical molecular dynamics simulations and 10,194 hypothetical MOFs created using the ToBaCCo 3.0 code. We found that high thermal conductivity in MOFs is favored by high densities (> 1.0 g cm−3), small pores (< 10 Å), and four-connected metal nodes. We also found that 36 MOFs exhibit ultra-low thermal conductivity (< 0.02 W m−1 K−1), which is primarily due to having extremely large pores (~65 Å). Furthermore, we discovered six hypothetical MOFs with very high thermal conductivity (> 10 W m−1 K−1), the structures of which we describe in additional detail.

Thermal transport in metal-organic frameworks (MOFs) is an essential but frequently overlooked property. Among the small number of existing studies on thermal transport in MOFs, even fewer have considered explicitly the influence of defects. However, defects naturally exist in MOF crystals and are known to influence many of their material properties. In this work, we investigate the influence of both randomly and symmetrically distributed defects on the thermal conductivity of the MOF UiO-66. Two types of defects were examined: missing linker and missing cluster defects. For symmetrically distributed (i.e., spatially correlated) defects, we considered three experimentally resolved defect nanodomains of UiO-66 with underlying topologies of bcu, reo, and scu. We observed that both randomly distributed missing linker and missing cluster defects typically decrease thermal conductivity, as expected. However, we found that the spatial arrangement of defects can significantly impact thermal conductivity. In particular, the spatially correlated missing linker defect nanodomain (bcu topology) displayed an intriguing anisotropy, with the thermal conductivity along a particular direction being higher than that of the defect-free UiO-66. We attribute this unusual defect-induced increase in thermal conductivity to the removal of the linkers perpendicular to the primary direction of heat transport. These perpendicular linkers act as phonon scattering sources such that removing them increases thermal conductivity in that direction. Moreover, we also observed an increase in phonon group velocity, which might also contribute to the unusual increase. Overall, we show that structural defects could be an additional lever to tune the thermal conductivity of MOFs.

The ability to control thermal transport is critical for the design of thermal rectifiers, logic gates, and transistors, although it remains a challenge to design materials that exhibit large changes in thermal conductivity with switching ratios suitable for practical applications. Here, we propose the use of flexible metal-organic frameworks, which can undergo significant structural changes in response to various stimuli, to achieve tunable switchable thermal conductivity. In particular, we use molecular dynamics simulations to show that the thermal conductivity of the flexible framework Fe(bdp) (bdp2- = 1,4-benzenedipyrazolate) becomes highly anisotropic upon transitioning from the expanded to the collapsed phase, with the conductivity decreasing by nearly an order of magnitude along the direction of compression. Our results add to a small but growing number of studies investigating metal-organic frameworks for thermal transport.

We experimentally and theoretically investigate the thermal conductivity and mechanical properties of polycrystalline HKUST-1 metal-organic frameworks (MOFs) infiltrated with three guest molecules: tetracyanoquinodimethane (TCNQ), 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane (F4-TCNQ), and (cyclohexane-1,4-diylidene)dimalononitrile (H4-TCNQ). This allows for modification of the interaction strength between the guest and host, presenting an opportunity to study the fundamental atomic scale mechanisms of how guest molecules impact the thermal conductivity of large unit cell porous crystals. The thermal conductivities of the guest@MOF systems decrease significantly, by on average a factor of 4, for all infiltrated samples as compared to the uninfiltrated, pristine HKUST-1. This reduction in thermal conductivity goes in tandem with an increase in density of 38% and corresponding increase in heat capacity of ∼48%, defying conventional effective medium scaling of thermal properties of porous materials. We explore the origin of this reduction by experimentally investigating the guest molecules' effects on the mechanical properties of the MOF and performing atomistic simulations to elucidate the roles of the mass and bonding environments on thermal conductivity. The reduction in thermal conductivity can be ascribed to an increase in vibrational scattering introduced by extrinsic guest-MOF collisions as well as guest molecule-induced modifications to the intrinsic vibrational structure of the MOF in the form of hybridization of low frequency modes that is concomitant with an enhanced population of localized modes. The concentration of localized modes and resulting reduction in thermal conductivity do not seem to be significantly affected by the mass or bonding strength of the guest species.

Mixtures of paraffin and carbon nanofillers have promising potential for thermal storage, as paraffin (the matrix) possesses high latent heat and the nanofillers compensate for the low thermal conductivity (TC) of paraffin. Understanding thermal transport in these materials is essential for practical applications, as weak thermal transport hinders fast charge/discharge of thermal energy. Here, we use non-equilibrium molecular dynamics (NEMD) simulations to study the interfacial thermal conductance (ITC) between graphene sheets and octadecane (C18H38) matrix under the limiting conditions of the sheets being parallel or perpendicular to the direction of the imposed heat flux. The results show that the systems containing thin graphene layers exhibit higher values of ITC. This study captures the asymptotic saturation of thermal conductance for the liquid phase of the perpendicular structure. Besides, given the greater number of structured layers of paraffin upon phase change, the ITC for the solid paraffin-graphene system is higher than the conductance of the liquid paraffin-graphene interface. We use the Pearson correlation coefficients of the vibrational power spectrum (VPS) of interfacing materials to explain the orders of magnitude variations of the observed ITC.

Diamine-appended variants of the metal-organic framework M2(dobpdc) (M = Mg, Mn, Fe, Co, Zn; dobpdc4- = 4,4'-dioxidobiphenyl-3,3'-dicarboxylate) exhibit exceptional CO2 capture properties owing to a unique cooperative adsorption mechanism, and thus hold promise for use in the development of energy- and cost-efficient CO2 separations. Understanding the nature of thermal transport in these materials is essential for such practical applications, however, as temperature rises resulting from exothermic CO2 uptake could potentially offset the energy savings offered by such cooperative adsorbents. Here, molecular dynamics (MD) simulations are employed in investigating thermal transport in bare and e-2-appended Zn2(dobpdc) (e-2 = N-ethylethylenediamine), both with and without CO2 as a guest. In the absence of CO2, the appended diamines function to enhance thermal conductivity in the ab-plane of e-2-Zn2(dobpdc) relative to the bare framework, as a result of noncovalent interactions between adjacent diamines that provide additional heat transfer pathways across the pore channel. Upon introduction of CO2, the thermal conductivity along the pore channel (the c-axis) increases due to the cooperative formation of metal-bound ammonium carbamates, which serve to create additional heat transfer pathways. In contrast, the thermal conductivity of the bare framework remains unchanged in the presence of zinc-bound CO2 but decreases in the presence of additional adsorbed CO2.

MOFs are of great interest for applications in gas storage, separation, catalysis and more recently, thermoelectrics (TE) due to their extremely high porosity (surface area), large chemical and mechanical tailorability defined by the choice of the metal node, the organic ligand, and infiltrated guest molecules, and low thermal conductivities. In particular, new emergent properties, such as electronic conductivity and energy transfer, have been achieved by infiltrating MOF pores with specialized ‘guest’ molecules. The thermal transport in infiltrated MOFs is critically important for understanding the efficiency of both gas adsorption and thermoelectric applications. Understanding the impacts on thermal transport upon MOF infiltration with an electrically conductive guest molecule is critical for realizing an efficient TE material. This work experimentally investigates the impacts on the thermal properties of HKUST-1 MOFs infiltrated with charge accepting TCNQ and F4-TCNQ guest molecules. We provide experimental evidence paired with Molecular Dynamics (MD) results to describe the thermal transport mechanism as a loading-dependent structural change within the MOF. Since HKUST-1 exhibits negative thermal expansion (NTE), we note the unusual phenomena that the material softens as it becomes denser with increasing temperatures. We show this is achieved by which the ligands take on a contorted configuration as they absorb more energy at higher temperatures. Upon infiltration, the MOF softens, however the softening (and therefore reduction in phonon group velocity) cannot be completely responsible for the large reduction in thermal conductivity that is experimentally observed (~71% at room temperature). Spectral Energy Density (SED) and MD results indicate that localized phonon modes from the adsorbates act as additional scattering mechanisms to reduce the thermal conductivity. Further, experimental evidence of the existence of low frequency ZA-like flexural modes (usually found only in 2D materials) within the highly porous 3D MOF structure likely contributes to a large portion of the thermal conductivity, and is suppressed upon infiltration with a guest molecule, causing the large reduction in thermal conductivity.

Whether the presence of adsorbates increases or decreases thermal conductivity in metal-organic frameworks (MOFs) has been an open question. Here we report observations of thermal transport in the metal-organic framework HKUST-1 in the presence of various liquid adsorbates: water, methanol, and ethanol. Experimental thermoreflectance measurements were performed on single crystals and thin films, and theoretical predictions were made using molecular dynamics simulations. We find that the thermal conductivity of HKUST-1 decreases by 40 – 80% depending on the adsorbate, a result that cannot be explained by effective medium approximations. Our findings demonstrate that adsorbates introduce additional phonon scattering in HKUST-1, which particularly shortens the lifetimes of low-frequency phonon modes. As a result, the system thermal conductivity is lowered to a greater extent than the increase expected by the creation of additional heat transfer channels. Finally, we show that thermal diffusivity is even more greatly reduced than thermal conductivity by adsorption.

Although metal-organic frameworks (MOFs) are promising materials for gas storage and separation applications, the heat released during the exothermic adsorption process can potentially negatively impact their practical utility. Thermal transport in MOFs has not been widely studied, and among the relatively few reports on the topic, MOFs have either been assumed to be defect free or the presence of defects was not discussed. However, defects naturally exist in MOFs and can also be introduced intentionally. Here, we investigate the effect of missing linker defects on the thermal conductivity of HKUST-1 using molecular dynamics (MDs) simulation and the Green-Kubo method. We found that missing linker defects, even at low concentrations, substantially reduce the thermal conductivity of HKUST-1. If not taken into account, the presence of defects could lead to significant discrepancies between experimentally measured and computationally predicted thermal conductivities.

We discuss the dependence of the propagon contribution to thermal conductivity on the medium range order (MRO) in amorphous silicon. Three different amorphous structures with the same size of 3.28 nm were studied. Among these three structures, two structures were constructed with experimentally observed MRO [M. M. J. Treacy and K. B. Borisenko, Science 335, 950 (2012)], and the other structure is from continuous random network (CRN), which lacks MRO and thus represents a randomized amorphous structure [G. Barkema and N. Mousseau, Phys. Rev. B 62, 4985 (2000)]. Using the simulated fluctuation electron microscopy and dihedral angle distribution, we confirm that the first two structures contain MRO in the length scale of 10–20 Å while the CRN structure does not. The transport of propagons in the MRO and CRN structures is compared using the dynamic structural factor calculation and normal mode decomposition of the molecular dynamics simulation data, showing noticeably longer lifetime of propagons in MRO structures than in the CRN structure. The propagon thermal conductivity in MRO structures is estimated to be 50% larger than that in the CRN structure.

With the aim of understanding heat transfer during structural changes in breathing porous crystals upon gas loading, we study the effect of pore expansion on the thermal conductivity of a series of idealized materials containing adsorbed gas using molecular dynamics simulations. We calculated the thermal conductivity in three main axes of the crystal lattice starting from a tilt angle of 40 degrees to represent the closed form of the crystal up to a tilt angle of 90 degrees to represent the open form. With no gas present, the thermal conductivity increases in the tilt direction with pore expansion whereas thermal conductivity in other directions remains unchanged. In the presence of adsorbed gas, porous crystals at all states of expansion experience reduced thermal conductivity due to phonon scattering introduced by gas-crystal interactions. Similarly, the thermal conductivity in the tilt direction increases as the pore expands; however, the increase is less pronounced compared to the case with no gas present in the pores. We also show that the diffusivity of gas increases during pore expansion, facilitating mass transport.

Graphene-C60 heterostructures assembled by van der Waals (vdW) interactions between graphene and C60 have shown exciting potential for multifunctional devices. Understanding thermal transport in graphene-C60 heterostructures is the key to guiding the design of vdW heterostructures with desired thermal transport properties. In this work, we report the first study of thermal transport in a graphene-C60 heterostructure and elucidate the importance of vdW interactions to heat conduction using molecular dynamics simulations. We find that the in-plane thermal conductivity of the graphene-C60 heterostructure is as high as about 234 W/(mK) at room temperature, exceeding those of most pure metals. As the vdW interaction parameter, chi, varies from 0.1 to 2, the in-plane thermal conductivity first increases then decreases. On the other hand, as vdW interactions increases, the interfacial thermal conductance between graphene and C60 is enhanced. Our study demonstrates that graphene-C60 heterostructures have high in-plane thermal conductivity and their interfacial thermal conductance is comparable to that of graphene-hexagonal boron-nitride (hBN) heterostructure. Graphene-C60 heterostructures are promising candidates for multifunctional devices with inherent heat dissipation capability. (C) 2019 Elsevier Ltd. All rights reserved.

The virial stress tensor-based instantaneous heat flux, which is used by LAMMPS, is only valid for the small subset of simulations that contain only pairwise interactions. For systems that contain many-body interactions using 3- or 4-body potentials, a more complete derivation is required. We have created a software patch to LAMMPS that implements the correct heat flux calculation approach for 3- and 4-body potentials, based on the derivation by Torii et al. ( J. Chem. Phys. 2008 , 128 , 044504 ) Using two example systems, the error in the uncorrected code for many-body potential heat flux is shown to be significant and reaches nearly 100% of the many-body potential heat flux for the systems we studied; hence, the error of the total heat flux calculation is proportional to the fraction of the total heat flux transferred through the many-body potentials. This error may have consequences for thermal conductivities calculated using the Green-Kubo method or any NEMD method that uses the instantaneous heat flux. We recommend that all researchers using LAMMPS for heat flux calculations where significant heat is transferred via the many-body potentials adopt the corrected code.

We report that the single interatomic potential, developed using Gaussian regression of data from density functional theory calculations, has high accuracy and flexibility to describe phonon transport with ab initio accuracy in two different atomistic configurations: perfect crystalline Si and crystalline Si with vacancies. The high accuracies of second- and third-order force constants from the Gaussian approximation potential (GAP) are demonstrated with phonon dispersion, Gruneisen parameter, three-phonon scattering rate, phonon-vacancy scattering rate, and thermal conductivity, all of which are very close to the results from density functional theory calculations. We also show that the widely used empirical potentials (Stillinger-Weber and Tersoff) produce much larger errors compared to the GAP. The computational cost of GAP is higher than the two empirical potentials, but five orders of magnitude lower than density functional theory calculations. Our work shows that GAP can provide a new opportunity for studying phonon transport in partially disordered crystalline phases with the high predictive power of ab initio calculation but at a feasible computational cost.

In practice, the usefulness of metal-organic frameworks (MOFs) for many gas storage applications depends on their ability to rapidly dissipate the heat generated during the exothermic adsorption process. MOFs can be precisely designed to have a wide variety of architectures which can allow tuning their thermal conductivity. In this work, we use molecular dynamics simulations to investigate the effect of interpenetration on the thermal conductivity of MOFs. We find that the addition of a parallel thermal transport pathway yields a thermal conductivity nearly the sum of the two constituent frameworks. This relationship holds for a variety of interpenetrating MOFs with different atomic masses and a wide range of interframework interaction parameters. We show that both the strength and range of interactions between constituent frameworks play a significant role in framework mobility as well as framework coupling which can result in deviation from this relationship. We propose a simple model to predict thermal conductivity of the interpenetrated framework based on thermal conductivities of individual frameworks and an interframework interaction parameter which can account for this deviation.

In practice, the usefulness of metal–organic frameworks (MOFs) for many gas storage applications depends on their ability to rapidly dissipate the heat generated during the exothermic adsorption process. MOFs can be precisely designed to have a wide variety of architectures which can allow tuning their thermal conductivity. In this work, we use molecular dynamics simulations to investigate the effect of pore geometry, gas loading, and interpenetration on the thermal conductivity of MOFs. We find that larger pores and increased gas loading reduce thermal conductivity. We also find that the addition of a parallel thermal transport pathway yields a thermal conductivity nearly the sum of the two constituent frameworks. This relationship holds for a variety of interpenetrating MOFs with different atomic masses and a wide range of interframework interaction parameters. We show that both the strength and range of interactions between constituent frameworks play a significant role in framework mobility as well as framework coupling which can result in deviation from this relationship. We propose a simple model to predict thermal conductivity of the interpenetrated framework based on thermal conductivities of individual frameworks and an interframework interaction parameter which can account for this deviation.

Methane adsorption into the metal-organic framework (MOF) HKUST-1 and the resulting heat generation and dissipation are investigated using molecular dynamics simulations. Transient simulations reveal that thermal transport in the MOF occurs two orders of magnitude faster than gas diffusion. A large thermal resistance at the MOF-gas interface (equivalent to 127 nm of bulk HKUST-1), however, prevents fast release of the generated heat. The mass transport resistance at the MOF-gas interface is equivalent to 1 nm of bulk HKUST-1 and does not present a bottleneck in the adsorption process. These results provide important insights into the application of MOFs for gas storage applications.

We have studied the mechanisms of heat transfer in a porous crystal-gas mixture system, motivated by the not insignificant challenge of quickly dissipating heat generated in metal-organic frameworks (MOFs) due to gas adsorption. Our study reveals that the thermal conductance of the system (crystal and gas) is dominated by lattice thermal conductivity in the crystal, and that conductance is reduced as the concentration of gas in the pores increases. This mechanism was observed from classical molecular simulations of a monatomic gas in an idealized porous crystal structure. We show that the decreased conductivity associated with increased gas concentration is due to phonon scattering in the crystal due to interactions with gas molecules. Calculations of scattering rates for two phonon modes reveal that scattering of the lowest frequency mode scales linearly with gas density. This result suggests that the probability of a phonon-gas collision is simply proportional to the number of gas molecules in the pore.