Exploring the Impact of the Linker Length on Heat Transport in Metal-Organic Frameworks.

Metal-organic frameworks (MOFs) are a highly versatile group of porous materials suitable for a broad range of applications, which often crucially depend on the MOFs' heat transport properties. Nevertheless, detailed relationships between the chemical structure of MOFs and their thermal conductivities are still largely missing. To lay the foundations for developing such relationships, we performed non-equilibrium molecular dynamics simulations to analyze heat transport in a selected set of materials. In particular, we focus on the impact of organic linkers, the inorganic nodes and the interfaces between them. To obtain reliable data, great care was taken to generate and thoroughly benchmark system-specific force fields building on ab-initio-based reference data. To systematically separate the different factors arising from the complex structures of MOF, we also studied a series of suitably designed model systems. Notably, besides the expected trend that longer linkers lead to a reduction in thermal conductivity due to an increase in porosity, they also cause an increase in the interface resistance between the different building blocks of the MOFs. This is relevant insofar as the interface resistance dominates the total thermal resistance of the MOF. Employing suitably designed model systems, it can be shown that this dominance of the interface resistance is not the consequence of the specific, potentially weak, chemical interactions between nodes and linkers. Rather, it is inherent to the framework structures of the MOFs. These findings improve our understanding of heat transport in MOFs and will help in tailoring the thermal conductivities of MOFs for specific applications.

Evaluating Computational Shortcuts in Supercell-Based Phonon Calculations of Molecular Crystals: The Instructive Case of Naphthalene.

Phonons crucially impact a variety of properties of organic semiconductor materials. For instance, charge- and heat transport depend on low-frequency phonons, while for other properties, such as the free energy, especially high-frequency phonons count. For all these quantities one needs to know the entire phonon band structure, whose simulation becomes exceedingly expensive for more complex systems when using methods like dispersion-corrected density functional theory (DFT). Therefore, in the present contribution we evaluate the performance of more approximate methodologies, including density functional tight binding (DFTB) and a pool of force fields (FF) of varying complexity and sophistication. Beyond merely comparing phonon band structures, we also critically evaluate to what extent derived quantities, like temperature-dependent heat capacities, mean squared thermal displacements, and temperature-dependent free energies are impacted by shortcomings in the description of the phonon bands. As a benchmark system, we choose (deuterated) naphthalene, as the only organic semiconductor material for which to date experimental phonon band structures are available in the literature. Overall, the best performance among the approximate methodologies is observed for a system-specifically parametrized second-generation force field. Interestingly, in the low-frequency regime also force fields with a rather simplistic model for the bonding interactions (like the General Amber Force Field) perform rather well. As far as the tested DFTB parametrization is concerned, we obtain a significant underestimation of the unit-cell volume resulting in a pronounced overestimation of the phonon energies in the low-frequency region. This cannot be mended by relying on the DFT-calculated unit cell, since with this unit cell the DFTB phonon frequencies significantly underestimate the experiments.

Surface Induced Phenytoin Polymorph. 1. Full Structure Solution by Combining Grazing Incidence X-ray Diffraction and Crystal Structure Prediction.

Understanding the behavior and properties of molecules assembled in thin layers requires knowledge of their crystalline packing. The drug phenytoin (5,5-diphenylhydantoin) is one of the compounds that can be grown as a surface induced polymorph. By using grazing incidence X-ray diffraction, the monoclinic unit cell of the new form II can be determined, but, due to crystal size and the low amount of data, a full solution using conventional structure solving strategies fails. In this work, the full solution has been obtained by combining computational structure generation and experimental results. The comparison between the bulk and the new surface induced phase reveals significant packing differences of the hydrogen-bonding network, which might be the reason for the faster dissolution of form II with respect to form I. The results are very satisfactory, and the method might be adapted for other systems, where, due to the limited amount of experimental data, one must rely on additional approaches to gain access to more detailed information to understand the solid-state behavior.

Surface Induced Phenytoin Polymorph. 2. Structure Validation by Comparing Experimental and Density Functional Theory Raman Spectra.

A method for structure solution in thin films that combines grazing incidence X-ray diffraction data analysis and crystal structure prediction was presented in a recent work (Braun et al. Cryst. Growth Des.2019, DOI: 10.1021/acs.cgd.9b00857). Applied to phenytoin form II, which is only detected in films, the approach gave a very reasonable, but not fully confirmed, candidate structure with Z = 4 and Z' = 2. In the present work, we demonstrate how, by calculating and measuring the crystal Raman spectrum in the low wavenumber energy region with the aim of validating the candidate structure, this can be further refined. In fact, we find it to correspond to a saddle point of the energy landscape of the system, from which a minimum of lower symmetry may be reached. With the new structure, with Z = 4 and Z' = 2, we finally obtain an excellent agreement between experimental and calculated Raman spectra.

Indexing of grazing-incidence X-ray diffraction patterns: the case of fibre-textured thin films.

Crystal structure solutions from thin films are often performed by grazing-incidence X-ray diffraction (GIXD) experiments. In particular, on isotropic substrates the thin film crystallites grow in a fibre texture showing a well defined crystallographic plane oriented parallel to the substrate surface with random in-plane order of the microcrystallites forming the film. In the present work, analytical mathematical expressions are derived for indexing experimental diffraction patterns, a highly challenging task which hitherto mainly relied on trial-and-error approaches. The six lattice constants a, b, c, α, ß and γ of the crystallographic unit cell are thereby determined, as well as the rotation parameters due to the unknown preferred orientation of the crystals with respect to the substrate surface. The mathematical analysis exploits a combination of GIXD data and information acquired by the specular X-ray diffraction. The presence of a sole specular diffraction peak series reveals fibre-textured growth with a crystallographic plane parallel to the substrate, which allows establishment of the Miller indices u, v and w as the rotation parameters. Mathematical expressions are derived which reduce the system of unknown parameters from the three- to the two-dimensional space. Thus, in the first part of the indexing routine, the integers u and v as well as the Laue indices h and k of the experimentally observed diffraction peaks are assigned by systematically varying the integer variables, and by calculating the three lattice parameters a, b and γ. Because of the symmetry of the derived equations, determining the missing parameters then becomes feasible: (i) w of the surface parallel plane, (ii) the Laue indices l of the diffraction peak and (iii) analogously the lattice constants c, α and ß. In a subsequent step, the reduced unit-cell geometry can be identified. Finally, the methodology is demonstrated by application to an example, indexing the diffraction pattern of a thin film of the organic semiconductor pentacenequinone grown on the (0001) surface of highly oriented pyrolytic graphite. The preferred orientation of the crystallites, the lattice constants of the triclinic unit cell and finally, by molecular modelling, the full crystal structure solution of the as-yet-unknown polymorph of pentacenequinone are determined.

Toward a Reliable Description of the Lattice Vibrations in Organic Molecular Crystals: The Impact of van der Waals Interactions.

This work assesses the reliability of different van der Waals (vdW) methods to describe lattice vibrations of molecular crystals in the framework of density functional theory (DFT). To accomplish this task, calculated and experimental lattice phonon Raman spectra of a pool of organic molecular crystals are compared. We show that the many-body dispersion (MBD@rsSCS) van der Waals method of Ambrosetti et al. and the pairwise method of Grimme et al. (D3-BJ) outperform the other tested approaches (i.e., the D2 method of Grimme, the TS method of Tkatchenko and Scheffler, and the nonlocal functional vdW-DF-optPBE of Klimes et al.). For the worse-performing approaches the results could not even be fixed by the introduction of scaling parameters, as commonly used for high-energy intramolecular vibrations. Interestingly, when using the experimentally determined unit cell parameters, DFT calculations using the PBE functional without corrections for long-range vdW interactions provide spectra of similar accuracy as the MBD@rsSCS and D3-BJ simulations.

DFT-Assisted Polymorph Identification from Lattice Raman Fingerprinting.

A combined experimental and theoretical approach, consisting of lattice phonon Raman spectroscopy and density functional theory (DFT) calculations, is proposed as a tool for lattice dynamics characterization and polymorph phase identification. To illustrate the reliability of the method, the lattice phonon Raman spectra of two polymorphs of the molecule 2,7-dioctyloxy[1]benzothieno[3,2-b]benzothiophene are investigated. We show that DFT calculations of the lattice vibrations based on the known crystal structures, including many-body dispersion van der Waals (MBD-vdW) corrections, predict experimental data within an accuracy of ≪5 cm-1 (≪0.6 meV). Due to the high accuracy of the simulations, they can be used to unambiguously identify different polymorphs and to characterize the nature of the lattice vibrations and their relationship to the structural properties. More generally, this work implies that DFT-MBD-vdW is a promising method to describe also other physical properties that depend on lattice dynamics like charge transport.