Source author record

Reinhard J. Maurer

Reinhard J. Maurer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2024arXiv

First-principles Nonadiabatic Dynamics of Molecules at Metal Surfaces with Vibrationally Coupled Electron Transfer

Accurate description of nonadiabatic dynamics of molecules at metal surfaces involving electron transfer has been a longstanding challenge for theory. Here, we tackle this problem by first constructing high-dimensional neural network diabatic potentials including state crossings determined by constrained density functional theory, then applying mixed quantum-classical surface hopping simulations to evolve coupled electron-nuclear motion. Our approach accurately describes the nonadiabatic effects in CO scattering from Au(111) without empirical parameters and yields results agreeing well with experiments under various conditions for this benchmark system. We find that both adiabatic and nonadiabatic energy loss channels have important contributions to the vibrational relaxation of highly vibrationally excited CO(vi = 17), whereas relaxation of low vibrationally excited states of CO(vi = 2) is weak and dominated by nonadiabatic energy loss. The presented approach paves the way for accurate first-principles simulations of electron transfer mediated nonadiabatic dynamics at metal surfaces.

preprint2022arXiv

Direct experimental evidence for substrate adatom incorporation into a molecular overlayer

While the phenomenon of metal substrate adatom incorporation into molecular overlayers is generally believed to occur in several systems, the experimental evidence for this relies on the interpretation of scanning tunnelling microscopy (STM) images, which can be ambiguous and provides no quantitative structural information. We show that surface X-ray diffraction (SXRD) uniquely provides unambiguous identification of these metal adatoms. We present the results of a detailed structural study of the Au(111)-F4TCNQ system, combining surface characterisation by STM, low energy electron diffraction and soft X-ray photoelectron spectroscopy with quantitative experimental structural information from normal incidence X-ray standing waves (NIXSW) and SXRD, together with dispersion corrected density functional theory (DFT) calculations. Excellent agreement is found between the NIXSW data and the DFT calculations regarding the height and conformation of the adsorbed molecule, which has a twisted geometry rather than the previously supposed inverted bowl shape. SXRD measurements provide unequivocal evidence for the presence and location of Au adatoms, while the DFT calculations show this reconstruction to be strongly energetically favoured.

preprint2022arXiv

Long-range dispersion-inclusive machine learning potentials for structure search and optimization of hybrid organic-inorganic interfaces

The computational prediction of the structure and stability of hybrid organic-inorganic interfaces provides important insights into the measurable properties of electronic thin film devices, coatings, and catalyst surfaces and plays an important role in their rational design. However, the rich diversity of molecular configurations and the important role of long-range interactions in such systems make it difficult to use machine learning (ML) potentials to facilitate structure exploration that otherwise require computationally expensive electronic structure calculations. We present an ML approach that enables fast, yet accurate, structure optimizations by combining two different types of deep neural networks trained on high-level electronic structure data. The first model is a short-ranged interatomic ML potential trained on local energies and forces, while the second is an ML model of effective atomic volumes derived from atoms-in-molecules partitioning. The latter can be used to connect short-range potentials to well-established density-dependent long-range dispersion correction methods. For two systems, specifically gold nanoclusters on diamond (110) surfaces and organic $π$-conjugated molecules on silver (111) surfaces, we train models on sparse structure relaxation data from density functional theory and show the ability of the models to deliver highly efficient structure optimizations and semi-quantitative energy predictions of adsorption structures.

preprint2022arXiv

NQCDynamics.jl: A Julia Package for Nonadiabatic Quantum Classical Molecular Dynamics in the Condensed Phase

Accurate and efficient methods to simulate nonadiabatic and quantum nuclear effects in high-dimensional and dissipative systems are crucial for the prediction of chemical dynamics in condensed phase. To facilitate effective development, code sharing and uptake of newly developed dynamics methods, it is important that software implementations can be easily accessed and built upon. Using the Julia programming language, we have developed the NQCDynamics.jl package which provides a framework for established and emerging methods for performing semiclassical and mixed quantum-classical dynamics in condensed phase. The code provides several interfaces to existing atomistic simulation frameworks, electronic structure codes, and machine learning representations. In addition to the existing methods, the package provides infrastructure for developing and deploying new dynamics methods which we hope will benefit reproducibility and code sharing in the field of condensed phase quantum dynamics. Herein, we present our code design choices and the specific Julia programming features from which they benefit. We further demonstrate the capabilities of the package on two examples of chemical dynamics in condensed phase: the population dynamics of the spin-boson model as described by a wide variety of semi-classical and mixed quantum-classical nonadiabatic methods and the reactive scattering of H2 on Ag(111) using the Molecular Dynamics with Electronic Friction method. Together, they exemplify the broad scope of the package to study effective model Hamiltonians and realistic atomistic systems.

preprint2022arXiv

Thermodynamic Driving Forces for Substrate Atom Extraction by Adsorption of Strong Electron Acceptor Molecules

A quantitative structural investigation is reported, aimed at resolving the issue of whether substrate adatoms are incorporated into the monolayers formed by strong molecular electron acceptors deposited onto metallic electrodes. A combination of normal-incidence X-ray standing waves, low energy electron diffraction, scanning tunnelling microscopy and X-ray photoelectron spectroscopy measurements demonstrate that the systems TCNQ and F4TCNQ on Ag(100) lie at the boundary between these two possibilities and thus represent ideal model systems with which to study this effect. A room-temperature commensurate phase of adsorbed TCNQ is found not to involve Ag adatoms, but to adopt an inverted bowl configuration, long predicted but not previously identified experimentally. By contrast, a similar phase of adsorbed F4TCNQ does lead to Ag adatom incorporation in the overlayer, the cyano endgroups of the molecule being twisted relative to the planar quinoid ring. Density functional theory (DFT) calculations show that this behaviour is consistent with the adsorption energetics. Annealing of the commensurate TCNQ overlayer phase leads to an incommensurate phase that does appear to incorporate Ag adatoms. Our results indicate that the inclusion (or exclusion) of metal atoms into the organic monolayers is the result of both thermodynamic and kinetic factors.

preprint2016arXiv

Ab-initio tensorial electronic friction for molecules on metal surfaces: nonadiabatic vibrational relaxation

Molecular adsorbates on metal surfaces exchange energy with substrate phonons and low-lying electron-hole pair excitations. In the limit of weak coupling, electron-hole pair excitations can be seen as exerting frictional forces on adsorbates that enhance energy transfer and facilitate vibrational relaxation or hot-electron mediated chemistry. We have recently reported on the relevance of tensorial properties of electronic friction [Phys. Rev. Lett. 116, 217601 (2016)] in dynamics at surfaces. Here we present the underlying implementation of tensorial electronic friction based on Kohn-Sham Density Functional Theory for condensed phase and cluster systems. Using local atomic-orbital basis sets, we calculate nonadiabatic coupling matrix elements and evaluate the full electronic friction tensor in the classical limit. Our approach is numerically stable and robust as shown by a detailed convergence analysis. We furthermore benchmark the accuracy of our approach by calculation of vibrational relaxation rates and lifetimes for a number of diatomic molecules at metal surfaces. We find friction-induced mode-coupling between neighboring CO adsorbates on Cu(100) in a c(2x2) overlayer to be important to understand experimental findings.

preprint2016arXiv

Adsorption structures and energetics of molecules on metal surfaces: Bridging experiment and theory

Adsorption geometry and stability of organic molecules on surfaces are key parameters that determine the observable properties and functions of hybrid inorganic/organic systems (HIOSs). Despite many recent advances in precise experimental characterization and improvements in first-principles electronic structure methods, reliable databases of structures and energetics for large adsorbed molecules are largely amiss. In this review, we present such a database for a range of molecules adsorbed on metal single-crystal surfaces. The systems we analyze include noble-gas atoms, conjugated aromatic molecules, carbon nanostructures, and heteroaromatic compounds adsorbed on five different metal surfaces. The overall objective is to establish a diverse benchmark dataset that enables an assessment of current and future electronic structure methods, and motivates further experimental studies that provide ever more reliable data. Specifically, the benchmark structures and energetics from experiment are here compared with the recently developed van der Waals (vdW) inclusive density-functional theory (DFT) method, DFT+vdW$^{\mathrm{surf}}$. In comparison to 23 adsorption heights and 17 adsorption energies from experiment we find a mean average deviation of 0.06 Å and 0.16 eV, respectively. This confirms the DFT+vdW$^{\mathrm{surf}}$ method as an accurate and efficient approach to treat HIOSs. A detailed discussion identifies remaining challenges to be addressed in future development of electronic structure methods, for which the here presented benchmark database may serve as an important reference.

preprint2016arXiv

Communication: Charge-Population Based Dispersion Interactions for Molecules and Materials

We introduce a system-independent method to derive effective atomic C$_6$ coefficients and polarizabilities in molecules and materials purely from charge population analysis. This enables the use of dispersion-correction schemes in electronic structure calculations without recourse to electron-density partitioning schemes and expands their applicability to semi-empirical methods and tight-binding Hamiltonians. We show that the accuracy of our method is en par with established electron-density partitioning based approaches in describing intermolecular C$_6$ coefficients as well as dispersion energies of weakly bound molecular dimers, organic crystals, and supramolecular complexes. We showcase the utility of our approach by incorporation of the recently developed many-body dispersion (MBD) method [Tkatchenko et al., Phys. Rev. Lett. 108, 236402 (2012)] into the semi-empirical Density Functional Tight-Binding (DFTB) method and propose the latter as a viable technique to study hybrid organic-inorganic interfaces.

preprint2016arXiv

Global Structure Search for Molecules on Surfaces: Efficient Sampling with Curvilinear Coordinates

Efficient structure search is a major challenge in computational materials science. We present a modification of the basin hopping global geometry optimization approach that uses a curvilinear coordinate system to describe global trial moves. This approach has recently been shown to be efficient in structure determination of clusters [Nano Letters 15, 8044-8048 (2015)] and is here extended for its application to covalent, complex molecules and large adsorbates on surfaces. The employed, automatically constructed delocalized internal coordinates are similar to molecular vibrations, which enhances the generation of chemically meaningful trial structures. By introducing flexible constraints and local translation and rotation of independent geometrical subunits we enable the use of this method for molecules adsorbed on surfaces and interfaces. For two test systems, trans-$β$-ionylideneacetic acid adsorbed on a Au(111) surface and methane adsorbed on a Ag(111) surface, we obtain superior performance of the method compared to standard optimization moves based on Cartesian coordinates.

preprint2016arXiv

Thermal and electronic fluctuations of flexible adsorbed molecules: Azobenzene on Ag(111)

We investigate the thermal and electronic collective fluctuations that contribute to the finite-temperature adsorption properties of flexible adsorbates on surfaces on the example of the molecular switch azobenzene C$_{12}$H$_{10}$N$_{2}$ on the Ag(111) surface. Using first-principles molecular dynamics simulations we obtain the free energy of adsorption that accurately accounts for entropic contributions, whereas the inclusion of many-body dispersion interactions accounts for the electronic correlations that govern the adsorbate binding. We find the adsorbate properties to be strongly entropy-driven, as can be judged by a kinetic molecular desorption prefactor of 10$^{24}$ s$^{-1}$ that largely exceeds previously reported estimates. We relate this effect to sizable fluctuations across structural and electronic observables. Comparison of our calculations to temperature-programmed desorption measurements demonstrates that finite-temperature effects play a dominant role for flexible molecules in contact with polarizable surfaces, and that recently developed first-principles methods offer an optimal tool to reveal novel collective behavior in such complex systems.

preprint2015arXiv

Interfacial charge rearrangement and intermolecular interactions: Density-functional theory study of free-base porphine adsorbed on Ag(111) and Cu(111)

We employ dispersion-corrected density-functional theory to study the adsorption of tetrapyrrole 2H-porphine (2H-P) at Cu(111) and Ag(111). Various contributions to adsorbate-substrate and adsorbate-adsorbate interactions are systematically extracted to analyze the self-assembly behavior of this basic building block to porphyrin-based metal-organic nanostructures. This analysis reveals a surprising importance of substrate-mediated van der Waals interactions between 2H-P molecules, in contrast to negligible direct dispersive interactions. The resulting net repulsive interactions rationalize the experimentally observed tendency for single molecule adsorption.

preprint2015arXiv

Many-body dispersion effects in the binding of adsorbates on metal surfaces

A correct description of electronic exchange and correlation effects for molecules in contact with extended (metal) surfaces is a challenging task for first-principles modeling. In this work we demonstrate the importance of collective van der Waals dispersion effects beyond the pairwise approximation for organic--inorganic systems on the example of atoms, molecules, and nanostructures adsorbed on metals. We use the recently developed many-body dispersion (MBD) approach in the context of density-functional theory [Phys. Rev. Lett. 108, 236402 (2012); J. Chem. Phys. 140, 18A508 (2014)] and assess its ability to correctly describe the binding of adsorbates on metal surfaces. We briefly review the MBD method and highlight its similarities to quantum-chemical approaches to electron correlation in a quasiparticle picture. In particular, we study the binding properties of xenon, 3,4,9,10-perylene-tetracarboxylic acid (PTCDA), and a graphene sheet adsorbed on the Ag(111) surface. Accounting for MBD effects we are able to describe changes in the anisotropic polarizability tensor, improve the description of adsorbate vibrations, and correctly capture the adsorbate--surface interaction screening. Comparison to other methods and experiment reveals that inclusion of MBD effects improves adsorption energies and geometries, by reducing the overbinding typically found in pairwise additive dispersion-correction approaches.

preprint2015arXiv

Spin Manipulation by Creation of Single-Molecule Radical Cations

All-trans-retinoic acid (ReA), a closed-shell organic molecule comprising only C, H, and O atoms, is investigated on a Au(111) substrate using scanning tunneling microscopy and spectroscopy. In dense arrays single ReA molecules are switched to a number of states, three of which carry a localized spin as evidenced by conductance spectroscopy in high magnetic fields. The spin of a single molecule may be reversibly switched on and off without affecting its neighbors. We suggest that ReA on Au is readily converted to a radical by the abstraction of an electron.

preprint2014arXiv

Temperature-dependent templated growth of porphine thin films on the (111) facets of copper and silver

The templated growth of the basic porphyrin unit, free-base porphine (2H-P), is characterized by means of X-ray photoelectron spectroscopy (XPS) and near-edge X-ray absorption fine-structure (NEXAFS) spectroscopy measurements and density functional theory (DFT). The DFT simulations allow the deconvolution of the complex XPS and NEXAFS signatures into contributions originating from five inequivalent carbon atoms, which can be grouped into C-N and C-C bonded species. Polarization-dependent NEXAFS measurements reveal an intriguing organizational behavior: On both Cu(111) and Ag(111), for coverages up to one monolayer, the molecules adsorb undeformed and parallel to the respective metal surface. Upon increasing the coverage, however, the orientation of the molecules in the thin films depends on the growth conditions. Multilayers deposited at low temperatures (LT) exhibit a similar average tilting angle (30 degree relative to the surface plane) on both substrates. Conversely, for multilayers grown at room temperature a markedly different scenario exists. On Cu(111) the film thickness is self-limited to a coverage of approximately two layers, while on Ag(111) multilayers can be grown easily and, in contrast to the bulk 2H-P crystal, the molecules are oriented perpendicular to the surface. This difference in molecular orientation results in a modified line-shape of the C 1s XPS signatures that is dependent on the incident photon energy, which is explained by comparison with depth-resolved DFT calculations. Simulations of ionization energies for differently stacked molecules show no indication for a packing-induced modification of the multilayer XP spectra, thus indicating that the comparison of single molecule calculations to multilayer data is justified.

preprint2013arXiv

Excited-state potential-energy surfaces of metal-adsorbed organic molecules from Linear Expansion Δ-Self-Consistent Field Density-Functional Theory (ΔSCF-DFT)

Accurate and efficient simulation of excited state properties is an important and much aspired cornerstone in the study of adsorbate dynamics on metal surfaces. To this end, the recently proposed linear expansion ΔSelf-Consistent Field (leΔSCF) method by Gavnholt et al. [Phys. Rev. B 78, 075441 (2008)] presents an efficient alternative to time consuming quasi-particle calculations. In this method the standard Kohn-Sham equations of Density-Functional Theory are solved with the constraint of a non-equilibrium occupation in a region of Hilbert-space resembling gas-phase orbitals of the adsorbate. In this work we discuss the applicability of this method for the excited-state dynamics of metal-surface mounted organic adsorbates, specifically in the context of molecular switching. We present necessary advancements to allow for a consistent quality description of excited-state potential-energy surfaces (PESs), and illustrate the concept with the application to Azobenzene adsorbed on Ag(111) and Au(111) surfaces. We find that the explicit inclusion of substrate electronic states modifies the topologies of intra-molecular excited-state PESs of the molecule due to image charge and hybridization effects. While the molecule in gas phase shows a clear energetic separation of resonances that induce isomerization and backreaction, the surface-adsorbed molecule does not. The concomitant possibly simultaneous induction of both processes would lead to a significantly reduced switching efficiency of such a mechanism.

preprint2012arXiv

Assessing computationally efficient isomerization dynamics: Delta-SCF density-functional theory study of azobenzene molecular switching

We present a detailed comparison of the S0, S1 (n -> π*) and S2 (π-> π*) potential energy surfaces (PESs) of the prototypical molecular switch azobenzene as obtained by Delta-self-consistent-field (Delta-SCF) Density-Functional Theory (DFT), time-dependent DFT (TD-DFT) and approximate Coupled Cluster Singles and Doubles (RI-CC2). All three methods unanimously agree in terms of the PES topologies, which are furthermore fully consistent with existing experimental data concerning the photo-isomerization mechanism. In particular, sum-method corrected Delta-SCF and TD-DFT yield very similar results for S1 and S2, when based on the same ground-state exchange-correlation (xc) functional. While these techniques yield the correct PES topology already on the level of semi-local xc functionals, reliable absolute excitation energies as compared to RI-CC2 or experiment require an xc treatment on the level of long-range corrected hybrids. Nevertheless, particularly the robustness of Delta-SCF with respect to state crossings as well as its numerical efficiency suggest this approach as a promising route to dynamical studies of larger azobenzene-containing systems.

preprint2012arXiv

Molecular simulation of surface reorganization and wetting in crystalline cellulose I and II

Cellulose is one of the most versatile substances in the world. Its immense variety of applications was in recent years complemented by interesting nanotechnological applications. The fabrication of complex cellulose-based materials asks for thorough understanding of the surface structure and its interactions with adsorbates. In this study we investigate several surface model systems of nanotechnological interest, which are obtained by reorganization of the cellulose-vacuum or cellulose-water interfaces of slabs of crystalline cellulose. To do this, we constructed supercells from crystallographic data, and then optimized the interface structures. From the bulk and surface systems we calculated structural properties such as unit cell parameters, dihedral conformation distributions, density profiles and hydrogen bonding. The results suggest that no overall geometrical restructuring occurs at the interface. However, the hydrogen bond network is strongly reconstructed, as is inferred from the dihedral conformations and hydrogen bond occurrences, although only within the first few layers. This holds for low index close packed structures as well as for high index loosely packed surfaces. Replacing the vacuum by ambient pressure water molecules we find less rearrangements of the cellulose surface, because the water allows formation of hydrogen bonds similar to those in the bulk phase. The water near the cellulose surface shows, however, strong structural changes. We observe reduced mobility of the water molecules, which corresponds to a cooling of water by about 30 degrees, in a slab that is about 10 Å thick. Although structuring and adsorption is observed on all surfaces, no actual penetration of water into the cellulose structure could be observed. This suggests that pure water is not sufficient to produce cellulose swelling at mesoscopic timescales.