Electrostatically Designing Materials and Interfaces.

Collective electrostatic effects arise from the superposition of electrostatic potentials of periodically arranged (di)polar entities and are known to crucially impact the electronic structures of hybrid interfaces. Here, it is discussed, how they can be used outside the beaten paths of materials design for realizing systems with advanced and sometimes unprecedented properties. The versatility of the approach is demonstrated by applying electrostatic design not only to metal-organic interfaces and adsorbed (complex) monolayers, but also to inter-layer interfaces in van der Waals heterostructures, to polar metal-organic frameworks (MOFs), and to the cylindrical pores of covalent organic frameworks (COFs). The presented design ideas are straightforward to simulate and especially for metal-organic interfaces also their experimental implementation has been amply demonstrated. For van der Waals heterostructures, the needed building blocks are available, while the required assembly approaches are just being developed. Conversely, for MOFs the necessary growth techniques exist, but more work on advanced linker molecules is required. Finally, COF structures exist that contain pores decorated with polar groups, but the electrostatic impact of these groups has been largely ignored so far. All this suggest that the dawn of the age of electrostatic design is currently experienced with potential breakthroughs lying ahead.

Designing Accurate Moment Tensor Potentials for Phonon-Related Properties of Crystalline Polymers.

The phonon-related properties of crystalline polymers are highly relevant for various applications. Their simulation is, however, particularly challenging, as the systems that need to be modeled are often too extended to be treated by ab initio methods, while classical force fields are too inaccurate. Machine-learned potentials parametrized against material-specific ab initio data hold the promise of being extremely accurate and also highly efficient. Still, for their successful application, protocols for their parametrization need to be established to ensure an optimal performance, and the resulting potentials need to be thoroughly benchmarked. These tasks are tackled in the current manuscript, where we devise a protocol for parametrizing moment tensor potentials (MTPs) to describe the structural properties, phonon band structures, elastic constants, and forces in molecular dynamics simulations for three prototypical crystalline polymers: polyethylene (PE), polythiophene (PT), and poly-3-hexylthiophene (P3HT). For PE, the thermal conductivity and thermal expansion are also simulated and compared to experiments. A central element of the approach is to choose training data in view of the considered use case of the MTPs. This not only yields a massive speedup for complex calculations while essentially maintaining DFT accuracy, but also enables the reliable simulation of properties that, so far, have been entirely out of reach.

Impact of hydrophilic side chains on the thin film transistor performance of a benzothieno-benzothiophene derivative.

Side-chain engineering in molecular semiconductors provides a versatile toolbox for precisely manipulating the material's processability, crystallographic properties, as well as electronic and optoelectronic characteristics. This study explores the impact of integrating hydrophilic side chains, specifically oligoethylene glycol (OEG) units, into the molecular structure of the small molecule semiconductor, 2,7-bis(2(2-methoxy ethoxy)ethoxy) benzo[b]benzo[4,5] thieno[2,3-d] thiophene (OEG-BTBT). The investigation includes a comprehensive analysis of thin film morphology and crystallographic properties, along with the optimization of deposition parameters for improving the device performance. Despite the anticipated benefits, such as enhanced processability, our investigation into OEG-BTBT-based organic field-effect transistors (OFETs) reveals suboptimal performance marked by a low effective charge carrier mobility, a low on/off ratio, and a high threshold voltage. The study unveils bias stress effects and device degradation attributed to the high ionization energy of OEG-BTBT alongside the hydrophilic nature of the ethylene-glycol moieties, which lead to charge trapping at the dielectric interface. Our findings underscore the need for a meticulous balance between electronic properties and chemical functionalities in molecular semiconductors to achieve stable and efficient performance in organic electronic devices.

Determining Structures of Layer-by-Layer Spin-Coated Zinc Dicarboxylate-Based Metal-Organic Thin Films.

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.

Predicting Spin-Dependent Phonon Band Structures of HKUST-1 Using Density Functional Theory and Machine-Learned Interatomic Potentials.

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.

Anisotropic Phonon Bands in H-Bonded Molecular Crystals: The Instructive Case of α-Quinacridone.

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).

Electrostatic Design of the Nanoscale Internal Surfaces of Porous Covalent Organic Frameworks.

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.

Sweep-Character-Dependent Switching of the Conductance State in Ferrocene-Substituted Thiofluorene Self-Assembled Monolayers.

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.

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.

Concept of Embedded Dipoles as a Versatile Tool for Surface Engineering.

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.

Discovering structure-property relationships for the phonon band structures of hydrocarbon-based organic semiconductor crystals: the instructive case of acenes.

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.

Maximizing the Carrier Mobilities of Metal-Organic Frameworks Comprising Stacked Pentacene Units.

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.

Avoiding the Center-Symmetry Trap: Programmed Assembly of Dipolar Precursors into Porous, Crystalline Molecular Thin Films.

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.

Porous Honeycomb Self-Assembled Monolayers: Tripodal Adsorption and Hidden Chirality of Carboxylate Anchored Triptycenes on Ag.

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.

Understanding the origin of serrated stacking motifs in planar two-dimensional covalent organic frameworks.

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.

First-principles calculations of hybrid inorganic-organic interfaces: from state-of-the-art to best practice.

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.

Mechanism of mediated alkali peroxide oxidation and triplet versus singlet oxygen formation.

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.

Strategies for Controlling Through-Space Charge Transport in Metal-Organic Frameworks via Structural Modifications.

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.

Electrostatic Design of Polar Metal-Organic Framework Thin Films.

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.

Final-State Simulations of Core-Level Binding Energies at Metal-Organic Hybrid Interfaces: Artifacts Caused by Spurious Collective Electrostatic Effects.

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.