Binding energies of benzene on coinage metal surfaces: Equal stability on different metals

Dispersion-Corrected DFT-D4 Study of the Adsorption of Halobenzenes and 1,3,5-Trihalobenzenes on the Cu(111) Surface─Effect of Sigma Hole Properties

Halogen bonds, characterized by directionality, tunability, hydrophobicity, and variable sizes, are ideal noncovalent interactions to design and control the formation of self-assembled nanostructures. The specific self-assembly cases formed by the halogen-bonding interaction have been well studied by scanning tunneling microscopy (STM) experiments and density functional theory (DFT) calculations. However, there is a lack of systematic theoretical adsorption studies on halogenated molecules. In this work, the adsorption of halobenzenes and 1,3,5-trihalobenzenes on the Cu(111) surface was examined by dispersion-corrected DFT methods. The adsorption geometries, noncovalent molecule-surface interactions, electronic densities, and electrostatic potential maps were examined for their most stable adsorption sites using the DFT-D4 method. Our calculations revealed that the iodo compounds favor a different adsorption geometry from aryl chlorides and bromides. Down the halogen group (Cl to I), the adsorption energy increases and the distance between the halogen atom and Cu surface decreases, which indicates stronger molecule-surface interactions. This is supported by the changes in the density of states upon adsorption. Noncovalent interaction analysis was also employed to further understand the nature and relative strength of the molecule-surface interactions. Electrostatic potential maps revealed that the positive character of the halogen sigma hole becomes stronger upon adsorption. Thus, surface adsorption of the halogenated molecule will enhance the formation of intermolecular halogen bonds. The present theoretical findings are expected to contribute toward a more comprehensive understanding of halogen bonding on the Cu(111) surface.

Temperature-programmed desorption of large molecules: influence of thin film structure and origin of intermolecular repulsion

Although the exact knowledge of the binding energy of organic adsorbates on solid surfaces is of vital importance for the realization of molecular nanostructures and the theoretical modelling of molecule-substrate interactions, an experimental determination is by no means trivial. Temperature-programmed desorption (TPD) is a widely used technique that can provide such information, but a quantitative analysis requires detailed knowledge of the pre-exponential factor of desorption and is therefore rarely performed on a quantitative level for larger molecules that often exhibit notable mutual intermolecular interactions. Here, we provide a thorough anlysis of TPD data of monolayers of pentacene and perfluoropentacene adsorbed on Au(111) that serve as a model system for polycyclic aromatic hydrocarbons adsorbed on noble metal surfaces. We show that the pre-exponential factor varies by several orders of magnitude with the surface coverage and evolves in a step-like fashion due to the sudden activation of a rotational degree of freedom during thermally controlled monolayer desorption. Using complementary coverage-dependent work function measurements, the interface dipole moments were determined. This allows to identify the origin and quantify the relative contributions of the lateral intermolecular interactions, which we modelled by force field calculations. This analysis clearly shows that the main cause for intermolecular repulsion are electrostatic interactions between the intramolecular charge distributions, while interface dipoles play only a minor role.

Crucial impact of exchange between layers on temperature programmed desorption

Desorption of molecules from surfaces constitutes an elementary process that is fundamental in both natural and application-oriented fields, including dewetting, weathering and catalysis. A powerful method to investigate desorption processes is temperature-programmed desorption (TPD), which offers the potential to provide mechanistic insights into the desorption kinetics. Using TPD, the desorption order, the energy barrier as well as the entropy change upon desorption can be accessed. In the past, several analysis methods have been developed for TPD data. These methods have in common that they rely on the Polanyi-Wigner equation, which requires proposing a desorption mechanism with a single (or at least dominating) desorption path. For real systems, however, several coupled desorption paths can be easily envisioned, which cannot be disentangled. Here, we analyse the influence of exchange between the first and the second adsorbate layer on the desorption process. We present a kinetic model, in which molecules can desorb directly from the first layer or change into the second layer and desorb from there. Interestingly, considering this additional desorption pathway alters the desorption spectrum considerably, even if the transient second-layer occupation remains as small as 4 × 10^-6 monolayers. We show that the impact of this layer exchange can be described by a modified Polanyi-Wigner equation. Our study demonstrates that layer exchange can crucially impact the TPD data.

Describing adsorption of benzene, thiophene, and xenon on coinage metals by using the Zaremba-Kohn theory-based model

Semilocal (SL) density functional approximations (DFAs) are widely applied but have limitations due to their inability to incorporate long-range van der Waals (vdW) interaction. Non-local functionals (vdW-DF, VV10, and rVV10) or empirical methods (DFT+D, DFT+vdW, and DFT+MBD) are used with SL-DFAs to account for such missing interaction. The physisorption of a molecule on the surface of the coinage metals (Cu, Ag, and Au) is a typical example of systems where vdW interaction is significant. However, it is difficult to find a general method that reasonably describes both adsorption energy and geometry of even the simple prototypes of cyclic and heterocyclic aromatic molecules such as benzene (C₆H₆) and thiophene (C₄H₄S), respectively, with reasonable accuracy. In this work, we present an alternative scheme based on Zaremba-Kohn theory, called DFT+vdW-dZK. We show that unlike other popular methods, DFT+vdW-dZK and particularly SCAN+vdW-dZK give an accurate description of the physisorption of a rare-gas atom (xenon) and two small albeit diverse prototype organic molecules on the (111) surfaces of the coinage metals.

r2SCAN-3c: A "Swiss army knife" composite electronic-structure method

The recently proposed r²SCAN meta-generalized-gradient approximation (mGGA) of Furness and co-workers is used to construct an efficient composite electronic-structure method termed r²SCAN-3c. To this end, the unaltered r²SCAN functional is combined with a tailor-made triple-ζ Gaussian atomic orbital basis set as well as with refitted D4 and geometrical counter-poise corrections for London-dispersion and basis set superposition error. The performance of the new method is evaluated for the GMTKN55 database covering large parts of chemical space with about 1500 data points, as well as additional benchmarks for non-covalent interactions, organometallic reactions, and lattice energies of organic molecules and ices, as well as for the adsorption on polar salt and non-polar coinage-metal surfaces. These comprehensive tests reveal a spectacular performance and robustness of r²SCAN-3c: It by far surpasses its predecessor B97-3c at only twice the cost and provides one of the best results of all semi-local density-functional theory (DFT)/QZ methods ever tested for the GMTKN55 database at one-tenth of the cost. Specifically, for reaction and conformational energies as well as non-covalent interactions, it outperforms prominent hybrid-DFT/QZ approaches at two to three orders of magnitude lower cost. Perhaps, the most relevant remaining issue of r²SCAN-3c is self-interaction error (SIE), owing to its mGGA nature. However, SIE is slightly reduced compared to other (m)GGAs, as is demonstrated in two examples. After all, this remarkably efficient and robust method is chosen as our new group default, replacing previous composite DFT and partially even expensive high-level methods in most standard applications for systems with up to several hundreds of atoms.

Engineering of TMDC-OSC hybrid interfaces: the thermodynamics of unitary and mixed acene monolayers on MoS2

Hybrid systems of two-dimensional (2D) materials such as transition metal dichalcogenides (TMDCs) and organic semiconductors (OSCs) have become subject of great interest for future device architectures. Although OSC-TMDC hybrid systems have been used in first device demonstrations, the precise preparation of ultra-thin OSC films on TMDCs has not been addressed. Due to the weak van der Waals interaction between TMDCs and OSCs, this requires precise knowledge of the thermodynamics at hand. Here, we use temperature-programmed desorption (TPD) and Monte Carlo (MC) simulations of TPD traces to characterize the desorption kinetics of pentacene (PEN) and perfluoropentacene (PFP) on MoS2 as a model system for OSCs on TMDCs. We show that the monolayers of PEN and PFP are thermally stabilized compared to their multilayers, which allows preparation of nominal monolayers by selective desorption of multilayers. This stabilization is, however, caused by entropy due to a high molecular mobility rather than an enhanced molecule-substrate bond. Consequently, the nominal monolayers are not densely packed films. Molecular mobility can be suppressed in mixed monolayers of PEN and PFP that, due to intermolecular attraction, form highly ordered films as shown by scanning tunneling microscopy. Although this reduces the entropic stabilization, the intermolecular attraction further stabilizes mixed films.

Next-Generation Nonlocal van der Waals Density Functional

The fundamental ideas for a nonlocal density functional theory-capable of reliably capturing van der Waals interactions-were already conceived in the 1990s. In 2004, a seminal paper introduced the first practical nonlocal exchange-correlation functional called vdW-DF, which has become widely successful and laid the foundation for much further research. However, since then, the functional form of vdW-DF has remained unchanged. Several successful modifications paired the original functional with different (local) exchange functionals to improve performance, and the successor vdW-DF2 also updated one internal parameter. Bringing together different insights from almost 2 decades of development and testing, we present the next-generation nonlocal correlation functional called vdW-DF3, in which we change the functional form while staying true to the original design philosophy. Although many popular functionals show good performance around the binding separation of van der Waals complexes, they often result in significant errors at larger separations. With vdW-DF3, we address this problem by taking advantage of a recently uncovered and largely unconstrained degree of freedom within the vdW-DF framework that can be constrained through empirical input, making our functional semiempirical. For two different parameterizations, we benchmark vdW-DF3 against a large set of well-studied test cases and compare our results with the most popular functionals, finding good performance in general for a wide array of systems and a significant improvement in accuracy at larger separations. Finally, we discuss the achievable performance within the current vdW-DF framework, the flexibility in functional design offered by vdW-DF3, as well as possible future directions for nonlocal van der Waals density functional theory.

Binding and electronic level alignment of π-conjugated systems on metals

We review the binding and energy level alignment of π-conjugated systems on metals, a field which during the last two decades has seen tremendous progress both in terms of experimental characterization as well as in the depth of theoretical understanding. Precise measurements of vertical adsorption distances and the electronic structure together with ab initio calculations have shown that most of the molecular systems have to be considered as intermediate cases between weak physisorption and strong chemisorption. In this regime, the subtle interplay of different effects such as covalent bonding, charge transfer, electrostatic and van der Waals interactions yields a complex situation with different adsorption mechanisms. In order to establish a better understanding of the binding and the electronic level alignment of π-conjugated molecules on metals, we provide an up-to-date overview of the literature, explain the fundamental concepts as well as the experimental techniques and discuss typical case studies. Thereby, we relate the geometric with the electronic structure in a consistent picture and cover the entire range from weak to strong coupling.

Stability and Exchange Processes in Ionic Liquid/Porphyrin Composite Films on Metal Surfaces

In light of increasing interest in the development of organic-organic multicomponent heterostructures on metals, this molecular-scale study investigates prototypical composite systems of ultrathin porphyrin and ionic liquid (IL) films on metallic supports under well-defined ultrahigh vacuum conditions. By means of angle-resolved X-ray photoelectron spectroscopy, we investigated the adsorption, stability, and thermal exchange of the resulting films after sequential physical vapor deposition of the free-base porphyrin 5,10,15,20-tetraphenylporphyrin, 2H-TPP, and the IL 1-methyl-3-octylimidazolium hexafluorophosphate, [C₈C₁Im][PF₆], on Ag(111) and Au(111). 2H-TPP shows two-dimensional growth of up to two closed molecular layers on Ag(111) and Au(111) and three-dimensional island growth for thicker films. IL films on top of a monolayer of 2H-TPP exhibit Stranski-Krastanov-like growth and are stable up to 385 K. The 2H-TPP layer leads to destabilization of the IL films, compared to the IL in direct contact with the bare metals, by inhibiting the specific adsorption of the ions on the metal surfaces. When the porphyrin is deposited on top of [C₈C₁Im][PF₆] at low temperature, the 2H-TPP molecules adsorb on top of the IL film at first but replace the IL at the IL/metal interfaces upon heating above 240 K. This exchange process is most likely driven by the higher adsorption energy of 2H-TPP on Ag(111) and Au(111) surfaces, as compared to the IL. The behavior observed on Ag(111) and Au(111) is identical. The results are highly relevant for the stability of porphyrin/IL-based thin film catalyst systems and molecular devices, and more generally, stacked organic multilayer architectures.

Nonadditivity of the Adsorption Energies of Linear Acenes on Au(111): Molecular Anisotropy and Many-Body Effects

Adsorption energies of chemisorbed molecules on inorganic solids usually scale linearly with molecular size and are well described by additive scaling laws. However, much less is known about scaling laws for physisorbed molecules. Our temperature-programmed desorption experiments demonstrate that the adsorption energy of acenes (benzene to pentacene) on the Au(111) surface in the limit of low coverage is highly nonadditive with respect to the molecular size. For pentacene, the deviation from an additive scaling of the adsorption energy amounts to as much as 0.7 eV. Our first-principles calculations explain the observed nonadditive behavior in terms of anisotropy of molecular polarization stemming from many-body electronic correlations. The observed nonadditivity of the adsorption energy has implications for surface-mediated intermolecular interactions and the ensuing on-surface self-assembly. Thus, future coverage-dependent studies should aim to gain insights into the impact of these complex interactions on the self-assembly of π-conjugated organic molecules on metal surfaces.

Understanding non-covalent interactions in larger molecular complexes from first principles

Non-covalent interactions pervade all matter and play a fundamental role in layered materials, biological systems, and large molecular complexes. Despite this, our accumulated understanding of non-covalent interactions to date has been mainly developed in the tens-of-atoms molecular regime. This falls considerably short of the scales at which we would like to understand energy trends, structural properties, and temperature dependencies in materials where non-covalent interactions have an appreciable role. However, as more reference information is obtained beyond moderately sized molecular systems, our understanding is improving and we stand to gain pertinent insights by tackling more complex systems, such as supramolecular complexes, molecular crystals, and other soft materials. In addition, accurate reference information is needed to provide the drive for extending the predictive power of more efficient workhorse methods, such as density functional approximations that also approximate van der Waals dispersion interactions. In this perspective, we discuss the first-principles approaches that have been used to obtain reference interaction energies for beyond modestly sized molecular complexes. The methods include quantum Monte Carlo, symmetry-adapted perturbation theory, non-canonical coupled cluster theory, and approaches based on the random-phase approximation. By considering the approximations that underpin each method, the most accurate theoretical references for supramolecular complexes and molecular crystals to date are ascertained. With these, we also assess a handful of widely used exchange-correlation functionals in density functional theory. The discussion culminates in a framework for putting into perspective the accuracy of high-level wavefunction-based methods and identifying future challenges.