Thin films of crystalline solids with substantial free volume built from organic chromophores and metal secondary building units (SBUs) are promising for engineering new optoelectronic properties through control of interchromophore coupling. Zn-based SBUs are especially relevant in this case because they avoid quenching the chromophore's luminescence. We find that layer-by-layer spin-coating using Zn acetate dihydrate and benzene-1,4-dicarboxylic acid (H2BDC) and biphenyl-4,4'-dicarboxylic acid (H2BPDC) linkers readily produces crystalline thin films. However, analysis of the grazing-incidence wide-angle X-ray scattering (GIWAXS) data reveals the structures of these films vary significantly with the linker, and with the metal-to-linker molar ratio used for fabrication. Under equimolar conditions, H2BPDC creates a type of structure like that proposed for SURMOF-2, whereas H2BDC generates a different metal-hydroxide-organic framework. Large excess of Zn2+ ions causes the growth of layered zinc hydroxides, irrespective of the linker used. Density functional theory (DFT) calculations provide structural models with minimum total energy that are consistent with the experimentally observed diffractograms. In the broader sense, this work illustrates the importance in this field of careful structure determination, e. g., by utilizing GIWAXS and DFT simulations to determine the structure of the obtained crystalline metal-organic thin films, such that properties can be rationally engineered and explained.
The present study focuses on the spin-dependent vibrational properties of HKUST-1, a metal-organic framework with potential applications in gas storage and separation. Employing density functional theory (DFT), we explore the consequences of spin couplings in the copper paddle wheels (as the secondary building units of HKUST-1) on the material's vibrational properties. By systematically screening the impact of the spin state on the phonon bands and densities of states in the various frequency regions, we identify asymmetric -COO- stretching vibrations as being most affected by different types of magnetic couplings. Notably, we also show that the DFT-derived insights can be quantitatively reproduced employing suitably parametrized, state-of-the-art machine-learned classical potentials with root-mean-square deviations from the DFT results between 3 cm-1 and 7 cm-1. This demonstrates the potential of machine-learned classical force fields for predicting the spin-dependent properties of complex materials, even when explicitly considering spins only for the generation of the reference data used in the force-field parametrization process.
Metal-Organic Frameworks , Density Functional Theory , Phonons , Algorithms
Phonons play a crucial role in the thermodynamic and transport properties of solid materials. Nevertheless, rather little is known about phonons in organic semiconductors. Thus, we employ highly reliable quantum mechanical calculations for studying the phonons in the α-polymorph of quinacridone. This material is particularly interesting, as it has highly anisotropic properties with distinctly different bonding types (H-bonding, π-stacking, and dispersion interactions) in different spatial directions. By calculating the overlaps of modes in molecular quinacridone and the α-polymorph, we associate Γ-point phonons with molecular vibrations to get a first impression of the impact of the crystalline environment. The situation becomes considerably more complex when analyzing phonons in the entire 1st Brillouin zone, where, due to the low symmetry of α-quinacridone, a multitude of avoided band crossings occur. At these, the character of the phonon modes typically switches, as can be inferred from mode participation ratios and mode longitudinalities. Notably, avoided crossings are observed not only as a function of the length but also as a function of the direction of the phonon wave vector. Analyzing these avoided crossings reveals how it is possible that the highest frequency acoustic band is always the one with the largest longitudinality, although longitudinal phonons in different crystalline directions are characterized by fundamentally different molecular displacements. The multiple avoided crossings also give rise to a particularly complex angular dependence of the group velocities, but combining the insights from the various studied quantities still allows drawing general conclusions, e.g., on the relative energetics of longitudinal vs transverse deformations (i.e., compressions and expansions vs slips of neighboring molecules). They also reveal how phonon transport in α-quinacridone is impacted by the reinforcing H-bonds and by π-stacking interactions (resulting from a complex superposition of van der Waals, charge penetration, and exchange repulsion).
It is well established that the collective action of assemblies of dipoles determines the electronic structure of surfaces and interfaces. This raises the question, to what extent the controlled arrangement of polar units can be used to also tune the electronic properties of the inner surfaces of materials with nanoscale pores. In the present contribution, state-of-the-art density-functional theory calculations are used to show for the prototypical case of covalent organic frameworks (COFs) that this is indeed possible. Decorating pore walls with assemblies of polar entities bonded to the building blocks of the COF layers triggers a massive change of the electrostatic energy within the pores. This, inevitably, also changes the relative alignment between electronic states in the framework and in guest molecules and is expected to have significant impacts on charge separation in COF heterojunctions, on redox reactions in COFs-based electrodes, and on (photo)catalysis.
Self-assembled monolayers (SAMs) of ferrocene-substituted thiofluorene on Au(111) exhibit two distinct conductance states (CSs) in two-terminal junctions featuring a sharp tip of eutectic GaIn as the top electrode. The occurrence of these states and the resulting effective rectification by the SAM depend on the way the bias voltage is swept; when the junction is only negatively biased, the original, high CS is preserved, whereas the junction is switched to a low CS when applying only positive biases. This results in an exceptionally high effective rectification ratio (RR) of â¼2100 already at voltages as low as 0.1 V. In contrast, when sweeping the junction alternatingly to the maximum positive and negative bias voltages (as usually performed in the literature), fully symmetric J-V curves are observed. That is, for the present SAM, rectification disappears, and the effective RR is ≈1. It is noteworthy that whether the junction in these symmetric sweeps is in the high or low CS depends on the polarity of the first sweep. We attribute the occurrence of the two CSs to a (quasi) non-reversible oxidation of the ferrocenes in combination with structural changes in the monolayer geometry. The observed sweeping dependence of the conductivity switching is an additional parameter that needs to be considered when interpreting experimental J-V curves, especially when dealing with redox-active systems.
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.
Controlling the physical and chemical properties of surfaces and interfaces is of fundamental relevance in various areas of physical chemistry and a key issue of modern nanotechnology. A highly promising strategy for achieving that control is the use of self-assembled monolayers (SAMs), which are ordered arrays of rodlike molecules bound to the substrate by a suitable anchoring group and carrying a functional tail group at the other end of the molecular backbone. Besides various other applications, SAMs are frequently used in organic electronics for the electrostatic engineering of interfaces by controlling the interfacial level alignment. This is usually achieved by introducing a dipolar tail group at the SAM-semiconductor interface. Such an approach, however, also changes the chemical character of that interface, for example, affecting the growth of subsequent layers. A strategy for avoiding this complication is to embed polar groups into the backbones of the SAM-forming molecules. This allows disentangling electronic interface engineering and the nucleation of further layers, such that both can be optimized independently. This novel concept was successfully demonstrated for both aliphatic and aromatic SAMs on different application-relevant substrates, such as gold, silver, and indium tin oxide. Embedding, for example, ester and pyrimidine groups in different orientations into the backbones of the SAM-forming molecules results in significant work-function changes. These can then be fine-tuned over a wide energy range by growing mixed monolayers consisting of molecules with oppositely oriented polar groups. In such systems, the variation of the work function is accompanied by pronounced shifts of the peaks in X-ray photoelectron spectra, which demonstrates that electrostatically triggered core-level shifts can be as important as the well-established chemical shifts. This illustrates the potential of X-ray photoelectron spectroscopy (XPS) as a tool for probing the local electrostatic energy within monolayers and, in systems like the ones studied here, makes XPS a powerful tool for studying the composition and morphology of binary SAMs. All these experimental observations can be rationalized through simulations, which show that the assemblies of embedded dipolar groups introduce a potential discontinuity within the monolayer, shifting the energy levels above and below the dipoles relative to each other. In molecular and monolayer electronics, embedded-dipole SAMs can be used to control transition voltages and current rectification. In devices based on organic and 2D semiconductors, such as MoS2, they can reduce contact resistances by several orders of magnitude without adversely affecting film growth even on flexible substrates. By varying the orientation of the embedded dipolar moieties, it is also possible to build p- and n-type organic transistors using the same electrode materials (Au). The extensions of the embedded-dipole concept from hybrid interfaces to systems such as metal-organic frameworks is currently underway, which further underlines the high potential of this approach.
Nanotechnology , Semiconductors , Electrodes , Electronics , Gold/chemistry , Nanotechnology/methods
By studying the low-frequency phonon bands of a series of crystalline acenes, this article lays the foundation for the development of structure-property relationships for phonons in organic semiconductors. Combining state-of-the art quantum-mechanical simulations with simple classical models, we explain how and why phonon frequencies and group velocities do or do not change when varying the molecular and crystal structures of the materials.
Liquid-phase, quasi-epitaxial growth is used to stack asymmetric, dipolar organic compounds on inorganic substrates, permitting porous, crystalline molecular materials that lack inversion symmetry. This allows material fabrication with built-in electric fields. A new programmed assembly strategy based on metal-organic frameworks (MOFs) is described that facilitates crystalline, noncentrosymmetric space groups for achiral compounds. Electric fields are integrated into crystalline, porous thin films with an orientation normal to the substrate. Changes in electrostatic potential are detected via core-level shifts of marker atoms on the MOF thin films and agree with theoretical results. The integration of built-in electric fields into organic, crystalline, and porous materials creates possibilities for band structure engineering to control the alignment of electronic levels in organic molecules. Built-in electric fields may also be used to tune the transfer of charges from donors loaded via programmed assembly into MOF pores. Applications include organic electronics, photonics, and nonlinear optics, since the absence of inversion symmetry results in a clear second-harmonic generation signal.
Charge transport properties of metal-organic frameworks (MOFs) are of distinct interest for (opto)electronic applications. In contrast to the situation in molecular crystals, MOFs allow an extrinsic control of the relative arrangement of π-conjugated entities through the framework architecture. This suggests that MOFs should enable materials with particularly high through-space charge carrier mobilities. Such materials, however, do not yet exist, despite the synthesis of MOFs with, for example, seemingly ideally packed stacks of pentacene-bearing linkers. Their rather low mobilities have been attributed to dynamic disorder effects. Using dispersion-corrected density functional theory calculations, we show that this is only part of the problem and that targeted network design involving comparably easy-to-implement structural modifications have the potential to massively boost charge transport. For the pentacene stacks, this is related to the a priori counterintuitive observation that the electronic coupling between neighboring units can be strongly increased by increasing the stacking distance.
Molecules with tripodal anchoring to substrates represent a versatile platform for the fabrication of robust self-assembled monolayers (SAMs), complementing the conventional monopodal approach. In this context, we studied the adsorption of 1,8,13-tricarboxytriptycene (Trip-CA) on Ag(111), mimicked by a bilayer of silver atoms underpotentially deposited on Au. While tripodal SAMs frequently suffer from poor structural quality and inhomogeneous bonding configurations, the triptycene scaffold featuring three carboxylic acid anchoring groups yields highly crystalline SAM structures. A pronounced polymorphism is observed, with the formation of distinctly different structures depending on preparation conditions. Besides hexagonal molecular arrangements, the occurrence of a honeycomb structure is particularly intriguing as such an open structure is unusual for SAMs consisting of upright-standing molecules. Advanced spectroscopic tools reveal an equivalent bonding of all carboxylic acid anchoring groups. Notably, density functional theory calculations predict a chiral arrangement of the molecules in the honeycomb network, which, surprisingly, is not apparent in experimental scanning tunneling microscopy (STM) images. This seeming discrepancy between theory and experiment can be resolved by considering the details of the actual electronic structure of the adsorbate layer. The presented results represent an exemplary showcase for the intricacy of interpreting STM images of complex molecular films. They are also further evidence for the potential of triptycenes as basic building blocks for generating well-defined layers with unusual structural motifs.
Covalent organic frameworks (COFs) have attracted significant attention due to their chemical versatility combined with a significant number of potential applications. Of particular interest are two-dimensional COFs, where the organic building units are linked by covalent bonds within a plane. Most properties of these COFs are determined by the relative arrangement of neighboring layers. These are typically found to be laterally displaced, which, for example, reduces the electronic coupling between the layers. In the present contribution we use dispersion-corrected density-functional theory to elucidate the origin of that displacement, showing that the common notion that the displacement is a consequence of electrostatic repulsions of polar building blocks can be misleading. For the representative case of COF-1 we find that electrostatic and van der Waals interactions would, actually, favor a cofacial arrangement of the layers and that Pauli repulsion is the crucial factor causing the serrated AA-stacking. A more in-depth analysis of the electrostatic contribution reveals that the "classical" Coulomb repulsion between the boroxine building blocks of COF-1 suggested by chemical intuition does exist, but is overcompensated by attractive effects due to charge-penetration in the phenylene units. The situation becomes more involved, when additionally allowing the interlayer distance to relax for each displacement, as then the different distance-dependences of the various types of interactions come into play. The overall behavior calculated for COF-1 is recovered for several additional COFs with differently sized π-systems and topologies, implying that the presented results are of more general relevance.
The computational characterization of inorganic-organic hybrid interfaces is arguably one of the technically most challenging applications of density functional theory. Due to the fundamentally different electronic properties of the inorganic and the organic components of a hybrid interface, the proper choice of the electronic structure method, of the algorithms to solve these methods, and of the parameters that enter these algorithms is highly non-trivial. In fact, computational choices that work well for one of the components often perform poorly for the other. As a consequence, default settings for one materials class are typically inadequate for the hybrid system, which makes calculations employing such settings inefficient and sometimes even prone to erroneous results. To address this issue, we discuss how to choose appropriate atomistic representations for the system under investigation, we highlight the role of the exchange-correlation functional and the van der Waals correction employed in the calculation and we provide tips and tricks how to efficiently converge the self-consistent field cycle and to obtain accurate geometries. We particularly focus on potentially unexpected pitfalls and the errors they incur. As a summary, we provide a list of best practice rules for interface simulations that should especially serve as a useful starting point for less experienced users and newcomers to the field.
Aprotic alkali metal-O2 batteries face two major obstacles to their chemistry occurring efficiently, the insulating nature of the formed alkali superoxides/peroxides and parasitic reactions that are caused by the highly reactive singlet oxygen (1O2). Redox mediators are recognized to be key for improving rechargeability. However, it is unclear how they affect 1O2 formation, which hinders strategies for their improvement. Here we clarify the mechanism of mediated peroxide and superoxide oxidation and thus explain how redox mediators either enhance or suppress 1O2 formation. We show that charging commences with peroxide oxidation to a superoxide intermediate and that redox potentials above ~3.5 V versus Li/Li+ drive 1O2 evolution from superoxide oxidation, while disproportionation always generates some 1O2. We find that 1O2 suppression requires oxidation to be faster than the generation of 1O2 from disproportionation. Oxidation rates decrease with growing driving force following Marcus inverted-region behaviour, establishing a region of maximum rate.
In recent years, optical and electronic properties of metal-organic frameworks (MOFs) have increasingly shifted into the focus of interest of the scientific community. Here, we discuss a strategy for conveniently tuning these properties through electrostatic design. More specifically, based on quantum-mechanical simulations, we suggest an approach for creating a gradient of the electrostatic potential within a MOF thin film, exploiting collective electrostatic effects. With a suitable orientation of polar apical linkers, the resulting non-centrosymmetric packing results in an energy staircase of the frontier electronic states reminiscent of the situation in a pin-photodiode. The observed one dimensional gradient of the electrostatic potential causes a closure of the global energy gap and also shifts core-level energies by an amount equaling the size of the original band gap. The realization of such assemblies could be based on so-called pillared layer MOFs fabricated in an oriented fashion on a solid substrate employing layer by layer growth techniques. In this context, the simulations provide guidelines regarding the design of the polar apical linker molecules that would allow the realization of MOF thin films with the (vast majority of the) molecular dipole moments pointing in the same direction.
In recent years, charge transport in metal-organic frameworks (MOFs) has shifted into the focus of scientific research. In this context, systems with efficient through-space charge transport pathways resulting from π-stacked conjugated linkers are of particular interest. In the current manuscript, we use density functional theory-based simulations to provide a detailed understanding of such MOFs, which, in the present case, are derived from the prototypical Zn2(TTFTB) system (with TTFTB4- corresponding to tetrathiafulvalene tetrabenzoate). In particular, we show that factors such as the relative arrangement of neighboring linkers and the details of the structural conformations of the individual building blocks have a profound impact on bandwidths and charge transfer. Considering the helical stacking of individual tetrathiafulvalene (TTF) molecules around a screw axis as the dominant symmetry element in Zn2(TTFTB)-derived materials, the focus, here, is primarily on the impact of the relative rotation of neighboring molecules. Not unexpectedly, changing the stacking distance in the helix also plays a distinct role, especially for structures which display large electronic couplings to start with. The presented results provide guidelines for achieving structures with improved electronic couplings. It is, however, also shown that structural defects (especially missing linkers) provide major obstacles to charge transport in the studied, essentially one-dimensional systems. This suggests that especially the sample quality is a decisive factor for ensuring efficient through-space charge transport in MOFs comprising stacked π-systems.
Core-level energies are frequently calculated to explain the X-ray photoelectron spectra of metal-organic hybrid interfaces. The current paper describes how such simulations can be flawed when modeling interfaces between physisorbed organic molecules and metals. The problem occurs when applying periodic boundary conditions to correctly describe extended interfaces and simultaneously considering core hole excitations in the framework of a final-state approach to account for screening effects. Since the core hole is generated in every unit cell, an artificial dipole layer is formed. In this work, we study methane on an Al(100) surface as a deliberately chosen model system for hybrid interfaces to evaluate the impact of this computational artifact. We show that changing the supercell size leads to artificial shifts in the calculated core-level energies that can be well beyond 1 eV for small cells. The same applies to atoms at comparably large distances from the substrate, encountered, for example, in extended, upright-standing adsorbate molecules. We also argue that the calculated work function change due to a core-level excitation can serve as an indication for the occurrence of such an artifact and discuss possible remedies for the problem.
In recent years, the photophysical properties of crystalline metal-organic frameworks (MOFs) have become increasingly relevant for their potential application in light-emitting devices, photovoltaics, nonlinear optics and sensing. The availability of high-quality experimental data for such systems makes them ideally suited for a validation of quantum mechanical simulations, aiming at an in-depth atomistic understanding of photophysical phenomena. Here we present a computational DFT study of the absorption and emission characteristics of a Zn-based surface-anchored metal-organic framework (Zn-SURMOF-2) containing anthracenedibenzoic acid (ADB) as linker. Combining band-structure and cluster-based simulations on ADB chromophores in various conformations and aggregation states, we are able to provide a detailed explanation of the experimentally observed photophysical properties of Zn-ADB SURMOF-2: The unexpected (weak) red-shift of the absorption maxima upon incorporating ADB chromophores into SURMOF-2 can be explained by a combination of excitonic coupling effects with conformational changes of the chromophores already in their ground state. As far as the unusually large red-shift of the emission of Zn-ADB SURMOF-2 is concerned, based on our simulations, we attribute it to a modification of the exciton coupling compared to conventional H-aggregates, which results from a relative slip of the centers of neighboring chromophores upon incorporation in Zn-ADB SURMOF-2.
Metal-Organic Frameworks/chemistry , Anthracenes/chemistry , Density Functional Theory , Dimerization , Molecular Conformation , Solvents/chemistry , Zinc/chemistry
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.
The molecular organization and electronic properties of dithiocarbamate (DTC) anchored self-assembled monolayers (SAMs) linked to Au(111) substrates are studied by a combination of X-ray photoelectron spectroscopy (XPS), near-edge X-ray absorption fine structure (NEXAFS) spectroscopy, and state-of-the-art density functional theory calculations. For that, several piperidine/piperazine precursors with different architecture and substitution patterns are selected. The presented data show that the DTC anchor provides a useful building block for monomolecular self-assembly on coinage metals with both sulfur atoms bonded to the substrate in a way similar to what is usually observed for the more commonly applied thiolate docking group. The combination of the DTC group with the quite flexible piperidine/piperazine cyclic linkers results in a dense molecular packing with an upright orientation of the terminal moieties. The latter comprise phenyl rings bearing various substituents, which enables tuning the interfacial dipole over a wide range. Simulations on two prototypical DTC-docked SAMs help to better understand the experimental observations and provide insight into the local origin of the SAM-induced shifts in the electrostatic energy. In particular, a comparison of measured and simulated XP spectra reveals the significant contribution of the DTC group to the interfacial dipole.
