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.

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.

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.

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.

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.

Exciton Coupling and Conformational Changes Impacting the Excited State Properties of Metal Organic Frameworks.

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.

Triptycene Tripods for the Formation of Highly Uniform and Densely Packed Self-Assembled Monolayers with Controlled Molecular Orientation.

When employing self-assembled monolayers (SAMs) for tuning surface and interface properties, organic molecules that enable strong binding to the substrate, large-area structural uniformity, precise alignment of functional groups, and control of their density are highly desirable. To achieve these goals, tripod systems bearing multiple bonding sites have been developed as an alternative to conventional monodentate systems. Bonding of all three sites has, however, hardly been achieved, with the consequence that structural uniformity and orientational order in tripodal SAMs are usually quite poor. To overcome that problem, we designed 1,8,13-trimercaptomethyltriptycene (T1) and 1,8,13-trimercaptotriptycene (T2) as potential tripodal SAM precursors and investigated their adsorption behavior on Au(111) combining several advanced experimental techniques and state-of-the-art theoretical simulations. Both SAMs adopt dense, nested hexagonal structures but differ in their adsorption configurations and structural uniformity. While the T2-based SAM exhibits a low degree of order and noticeable deviation from the desired tripodal anchoring, all three anchoring groups of T1 are equally bonded to the surface as thiolates, resulting in an almost upright orientation of the benzene rings and large-area structural uniformity. These superior properties are attributed to the effect of conformationally flexible methylene linkers at the anchoring groups, absent in the case of T2. Both SAMs display interesting electronic properties, and, bearing in mind that the triptycene framework can be functionalized by tail groups in various positions and with high degree of alignment, especially T1 appears as an ideal docking platform for complex and highly functional molecular films.

A dithiocarbamate anchoring group as a flexible platform for interface engineering.

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.

Tunneling Probability Increases with Distance in Junctions Comprising Self-Assembled Monolayers of Oligothiophenes.

Molecular tunneling junctions should enable the tailoring of charge-transport at the quantum level through synthetic chemistry but are hindered by the dominance of the electrodes. We show that the frontier orbitals of molecules can be decoupled from the electrodes, preserving their relative energies in self-assembled monolayers even when a top-contact is applied. This decoupling leads to the remarkable observation of tunneling probabilities that increase with distance in a series of oligothiophenes, which we explain using a two-barrier tunneling model. This model is generalizable to any conjugated oligomers for which the frontier orbital gap can be determined and predicts that the molecular orbitals that dominate tunneling charge-transport can be positioned via molecular design rather than by domination of Fermi-level pinning arising from strong hybridization. The ability to preserve the electronic structure of molecules in tunneling junctions facilitates the application of well-established synthetic design rules to tailor the properties of molecular-electronic devices.

Unconventional Current Scaling and Edge Effects for Charge Transport through Molecular Clusters.

Metal-molecule-metal junctions are the key components of molecular electronics circuits. Gaining a microscopic understanding of their conducting properties is central to advancing the field. In the present contribution, we highlight the fundamental differences between single-molecule and ensemble junctions focusing on the fundamentals of transport through molecular clusters. In this way, we elucidate the collective behavior of parallel molecular wires, bridging the gap between single molecule and large-area monolayer electronics, where even in the latter case transport is usually dominated by finite-size islands. On the basis of first-principles charge-transport simulations, we explain why the scaling of the conductivity of a junction has to be distinctly nonlinear in the number of molecules it contains. Moreover, transport through molecular clusters is found to be highly inhomogeneous with pronounced edge effects determined by molecules in locally different electrostatic environments. These effects are most pronounced for comparably small clusters, but electrostatic considerations show that they prevail also for more extended systems.

Exploring the driving forces behind the structural assembly of biphenylthiolates on Au(111).

In this contribution, we use dispersion-corrected density functional theory to study inter- and intramolecular interactions in a prototypical self-assembled monolayer (SAM) consisting of biphenylthiolates bonded to Au(111) via thiolate groups. The goal is to identify the nature of the interactions that drive the monolayer into a specific conformation. Particular focus is laid on sampling realistic structures rather than high symmetry model configurations. This is achieved by studying conceptually different local minimum structures of the SAM that are obtained via exploring the potential energy surface from systematically varied starting geometries. The six obtained packing motifs differ in the relative arrangement of the two molecules in the unit cell (co-planar versus herringbone) and in the intramolecular configuration (twisted versus planar rings). We find that van der Waals interactions within the organic adsorbate and between the adsorbate and substrate are the main reason that these molecular assemblies can form stable structures at all. The van der Waals interactions are, however, very similar for all observed motifs; by analyzing various types of interactions in the course of three notional SAM-formation steps, we find that the main driving force stabilizing the actual global minimum structure originates from electrostatic interactions between the molecules.

Understanding the adsorption of CuPc and ZnPc on noble metal surfaces by combining quantum-mechanical modelling and photoelectron spectroscopy.

Phthalocyanines are an important class of organic semiconductors and, thus, their interfaces with metals are both of fundamental and practical relevance. In the present contribution we provide a combined theoretical and experimental study, in which we show that state-of-the-art quantum-mechanical simulations are nowadays capable of treating most properties of such interfaces in a quantitatively reliable manner. This is shown for Cu-phthalocyanine (CuPc) and Zn-phthalocyanine (ZnPc) on Au(111) and Ag(111) surfaces. Using a recently developed approach for efficiently treating van der Waals (vdW) interactions at metal/organic interfaces, we calculate adsorption geometries in excellent agreement with experiments. With these geometries available, we are then able to accurately describe the interfacial electronic structure arising from molecular adsorption. We find that bonding is dominated by vdW forces for all studied interfaces. Concomitantly, charge rearrangements on Au(111) are exclusively due to Pauli pushback. On Ag(111), we additionally observe charge transfer from the metal to one of the spin-channels associated with the lowest unoccupied π-states of the molecules. Comparing the interfacial density of states with our ultraviolet photoelectron spectroscopy (UPS) experiments, we find that the use of a hybrid functionals is necessary to obtain the correct order of the electronic states.

X-ray based tools for the investigation of buried interfaces in organic electronic devices.

X-ray reflectivity combined with grazing incidence diffraction is a valuable tool for investigating organic multilayer structures that can be used in devices. We focus on a bilayer stack consisting of two materials (poly-(3-hexylthiophene)) (P3HT) and poly-(4-styrenesulfonic acid) (PSSA) spin cast from orthogonal solvents (water in the case of PSSA and chloroform or toluene for P3HT). X-ray reflectivity is used to determine the thickness of all layers as well as the roughness of the organic-organic hetero-interface and the P3HT surface. The surface roughness is found to be consistent with the results of atomic force microscopy measurements. For the roughness of P3HT/PSSA interface, we observe a strong dependence on the solvent used for P3HT deposition. The solvent also strongly impacts the texturing of the P3HT crystallites as revealed by grazing incidence diffraction. When applying the various PSSA/P3HT multilayers in organic thin-film transistors, we find an excellent correlation between the determined interface morphology, structure and the device performance.

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

Density-functional theory with screened van der Waals interactions for the modeling of hybrid inorganic-organic systems.

The electronic properties and the function of hybrid inorganic-organic systems (HIOS) are intimately linked to their interface geometry. Here we show that the inclusion of the many-body collective response of the substrate electrons inside the inorganic bulk enables us to reliably predict the HIOS geometries and energies. This is achieved by the combination of dispersion-corrected density-functional theory (the DFT+ van der Waals approach) [Phys. Rev. Lett. 102, 073005 (2009)], with the Lifshitz-Zaremba-Kohn theory for the nonlocal Coulomb screening within the bulk. Our method yields geometries in remarkable agreement (≈0.1 Å) with normal incidence x-ray standing wave measurements for the 3, 4, 9, 10-perylene-tetracarboxylic acid dianhydride (C(24)O(6)H(8), PTCDA) molecule on Cu(111), Ag(111), and Au(111) surfaces. Similarly accurate results are obtained for xenon and benzene adsorbed on metal surfaces.

Dimensionality effects in the electronic structure of organic semiconductors consisting of polar repeat units.

In conjugated organic molecules, excitation gaps typically decrease reciprocally with increasing the number of repeat units, n. This usually holds for individual molecules as well as for the corresponding bulk materials. Here, we show using density-functional theory calculations that a qualitatively different evolution is found for layers built from molecules consisting of polar repeat units. Whereas a 1/n-dependence is still observed in the case of isolated polar molecules, the global gap decreases essentially linearly with n in the corresponding 2D-periodic systems and vanishes beyond a certain molecular length, with the frontier states being localized at opposite ends of the layer. The latter is accompanied by a saturation of the dipole moment per molecule, an effect not observed in the isolated polar molecules. Interestingly, in both cases the limit of the gap for long (but finite) molecules differs qualitatively from that of infinite length obtained in 1D-periodic and 3D-periodic calculations, the latter serving as models for polymers and the bulk. We rationalize these dimensionality effects as a consequence of the potential gradient within the finite-length layers. They arise from the collective action of intra-molecular dipoles in the 2D periodic layers and can be traced back to surface effects.

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.

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.

Collectively induced quantum-confined Stark effect in monolayers of molecules consisting of polar repeating units.

The electronic structure of terpyrimidinethiols is investigated by means of density-functional theory calculations for isolated molecules and monolayers. In the transition from molecule to self-assembled monolayer (SAM), we observe that the band gap is substantially reduced, frontier states increasingly localize on opposite sides of the SAM, and this polarization in several instances is in the direction opposite to the polarization of the overall charge density. This behavior can be analyzed by analogy to inorganic semiconductor quantum-wells, which, as the SAMs studied here, can be regarded as semiperiodic systems. There, similar observations are made under the influence of a, typically external, electric field and are known as the quantum-confined Stark effect. Without any external perturbation, in oligopyrimidine SAMs one encounters an energy gradient that is generated by the dipole moments of the pyrimidine repeat units. It is particularly strong, reaching values of about 1.6 eV/nm, which corresponds to a substantial electric field of 1.6 × 10(7) V/cm. Close-lying σ- and π-states turn out to be a particular complication for a reliable description of the present systems, as their order is influenced not only by the docking groups and bonding to the metal, but also by the chosen computational approach. In the latter context we demonstrate that deliberately picking a hybrid functional allows avoiding pitfalls due to the infamous self-interaction error. Our results show that when aiming to build a monolayer with a specific electronic structure one can not only resort to the traditional technique of modifying the molecular structure of the constituents, but also try to exploit collective electronic effects.