The thermal-expansion properties [1] of substances are very important in materials design; for example, cracks form when joined materials expand or contract by different amounts upon heating. The most famous example of a substance that contracts when heated is ice: it transforms into water, which has a higher density than ice. Negative thermal expansion (NTE) in solids is relatively rare, although examples have been found in zeolites.[2] The underlying physics of NTE remains poorly understood. Herein, we show, on the basis of molecular simulations, that the recently synthesized isoreticular metal-organic frameworks (IRMOFs) consistently have negative thermal-expansion coefficients and are by far the most contracting materials known. Our simulations point to two competing effects: a local effect, where all bond lengths increase with temperature, and a second long-range effect, where the thermal movement of the linker molecules leads to a shorter average distance between corners upon heating.MOFs are a new class of nanoporous materials that have good stability, large void volumes, and well-defined tailorable cavities of uniform size. Their potential appears great, because these are precisely the properties needed for catalysis, separations, and storage/release applications. [3] MOFs generally consist of metal or metal-oxygen vertices interconnected by rigid or semirigid organic molecules. A large variety of MOFs, featuring different linker molecules and different types of bonding between the vertices with the linkers, have been produced by various research groups. The specific examples shown in Figure 1 are IRMOFs developed by Yaghi and co-workers. [4][5][6][7][8] In general, the IRMOFs consist of zinc-oxygen complexes connected by carboxylate-terminated linkers, forming a three-dimensional lattice of cubic cavities.Molecular simulations of adsorption in MOFs have shown very good agreement with experiment, [9][10][11][12][13] and it is interesting to note that simulations of diffusion in MOFs preceded experiments by almost two years. [9,14] In addition to predicting macroscopic observables, simulations can also provide useful molecular-level insights. To systematically investigate the thermal properties of MOFs, we herein simulate the (cubic) structures of several IRMOFs of varying linker length. We obtain information about the unit-cell length (L) as a function of temperature and about adsorbate loadings (q) as a function of pressure. We show that the experimental data scattered in the literature (for different adsorbate loadings and temperatures) are, in fact, consistent, and we elucidate the different effects of temperature and loading on L at the microscopic level.Various models for MOF flexibility have recently appeared. [15][16][17] Our flexible framework model for IRMOF-1, IRMOF-10, and IRMOF-16 is described in the Supporting Information. It is similar in spirit to the model of Greathouse and Allendorf, [15] but differs in the treatment of the carboxylate group and has the advantage of being calibrated to ex...