Efficient implementation and performance analysis of the independent electron surface hopping method for dynamics at metal surfaces.

Independent electron surface hopping (IESH) is a computational algorithm for simulating the mixed quantum-classical molecular dynamics of adsorbate atoms and molecules interacting with metal surfaces. It is capable of modeling the nonadiabatic effects of electron-hole pair excitations on molecular dynamics. Here, we present a transparent, reliable, and efficient implementation of IESH, demonstrating its ability to predict scattering and desorption probabilities across a variety of systems, ranging from model Hamiltonians to full dimensional atomistic systems. We further show how the algorithm can be modified to account for the application of an external bias potential, comparing its accuracy to results obtained using the hierarchical quantum master equation. Our results show that IESH is a practical method for modeling coupled electron-nuclear dynamics at metal surfaces, especially for highly energetic scattering events.

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 the 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 the 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 the condensed phase: the population dynamics of the spin-boson model as described by a wide variety of semiclassical 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.

On the optical anisotropy in 2D metal-halide perovskites.

Two-dimensional metal-halide perovskites (MHPs) are versatile solution-processed organic/inorganic quantum wells where the structural anisotropy creates profound anisotropy in their electronic and excitonic properties and associated optical constants. We here employ a wholistic framework, based on semiempirical modeling (k·p/effective mass theory calculations) informed by hybrid density functional theory (DFT) and multimodal spectroscopic ellipsometry on (C6H5(CH2)2NH3)2PbI4 films and crystals, that allows us to link the observed optical properties and anisotropy precisely to the underlying physical parameters that shape the electronic structure of a layered MHP. We find substantial frequency-dependent anisotropy in the optical constants and close correspondence between experiment and theory, demonstrating a high degree of in-plane alignment of the two-dimensional planes in both spin-coated thin films and cleaved single crystals made in this study. Hybrid DFT results elucidate the degree to which organic and inorganic frontier orbitals contribute to optical transitions polarized along a particular axis. The combined experimental and theoretical approach enables us to estimate the fundamental electronic bandgap of 2.65-2.68 eV in this prototypical 2D perovskite and to determine the spin-orbit coupling (ΔSO = 1.20 eV) and effective crystal field (δ = -1.36 eV) which break the degeneracy of the frontier conduction band states and determine the exciton fine structure. The methods and results described here afford a better understanding of the connection between structure and induced optical anisotropy in quantum-confined MHPs, an important structure-property relationship for optoelectronic applications and devices.

Organic-to-inorganic structural chirality transfer in a 2D hybrid perovskite and impact on Rashba-Dresselhaus spin-orbit coupling.

Translation of chirality and asymmetry across structural motifs and length scales plays a fundamental role in nature, enabling unique functionalities in contexts ranging from biological systems to synthetic materials. Here, we introduce a structural chirality transfer across the organic-inorganic interface in two-dimensional hybrid perovskites using appropriate chiral organic cations. The preferred molecular configuration of the chiral spacer cations, R-(+)- or S-(-)-1-(1-naphthyl)ethylammonium and their asymmetric hydrogen-bonding interactions with lead bromide-based layers cause symmetry-breaking helical distortions in the inorganic layers, otherwise absent when employing a racemic mixture of organic spacers. First-principles modeling predicts a substantial bulk Rashba-Dresselhaus spin-splitting in the inorganic-derived conduction band with opposite spin textures between R- and S-hybrids due to the broken inversion symmetry and strong spin-orbit coupling. The ability to break symmetry using chirality transfer from one structural unit to another provides a synthetic design paradigm for emergent properties, including Rashba-Dresselhaus spin-polarization for hybrid perovskite spintronics and related applications.

Frenkel-Holstein Hamiltonian applied to absorption spectra of quaterthiophene-based 2D hybrid organic-inorganic perovskites.

For the prototypical two-dimensional hybrid organic-inorganic perovskites (2D HOIPs) (AE4T)PbX4 (X = Cl, Br, and I), we demonstrate that the Frenkel-Holstein Hamiltonian (FHH) can be applied to describe the absorption spectrum arising from the organic component. We first model the spectra using only the four nearest neighbor couplings between translationally inequivalent molecules in the organic herringbone lattice as fitting parameters in the FHH. We next use linear-response time-dependent density functional theory (LR-TDDFT) to calculate molecular transition densities, from which extended excitonic couplings are evaluated based on the atomic positions within the 2D HOIPs. We find that both approaches reproduce the experimentally observed spectra, including changes in their shape and peak positions. The spectral changes are correlated with a decrease in excitonic coupling from X = Cl to X = I. Importantly, the LR-TDDFT-based approach with extended excitonic couplings not only gives better agreement with the experimental absorption line shape than the approach using a restricted set of fitted parameters but also allows us to relate the changes in excitonic coupling to the underlying geometry. We accordingly find that the decrease in excitonic coupling from X = Cl to Br to I is due to an increase in molecular separation, which in turn can be related to the increasing Pb-X bond length from Cl to I. Our research opens up a potential pathway to predicting optoelectronic properties of new 2D HOIPs from ab initio calculations and to gain insight into structural relations from 2D HOIP absorption spectra.

Molecular engineering of organic-inorganic hybrid perovskites quantum wells.

Semiconductor quantum-well structures and superlattices are key building blocks in modern optoelectronics, but it is difficult to simultaneously realize defect-free epitaxial growth while fine tuning the chemical composition, layer thickness and band structure of each layer to achieve the desired performance. Here we demonstrate the modulation of the electronic structure-and consequently the optical properties-of organic semiconducting building blocks that are incorporated between the layers of perovskites through a facile solution processing step. Self-aggregation of the conjugated organic molecules is suppressed by functionalization with sterically demanding groups and single crystalline organic-perovskite hybrid quantum wells (down to one-unit-cell thick) are obtained. The energy and charge transfers between adjacent organic and inorganic layers are shown to be fast and efficient, owing to the atomically flat interface and ultrasmall interlayer distance of the perovskite materials. The resulting two-dimensional hybrid perovskites are very stable due to protection given by the bulky hydrophobic organic groups.

Direct-Bandgap 2D Silver-Bismuth Iodide Double Perovskite: The Structure-Directing Influence of an Oligothiophene Spacer Cation.

Three-dimensional (3D) hybrid organic-inorganic lead halide perovskites (HOIPs) feature remarkable optoelectronic properties for solar energy conversion but suffer from long-standing issues of environmental stability and lead toxicity. Associated two-dimensional (2D) analogues are garnering increasing interest due to superior chemical stability, structural diversity, and broader property tunability. Toward lead-free 2D HOIPs, double perovskites (DPs) with mixed-valent dual metals are attractive. Translation of mixed-metal DPs to iodides, with their prospectively lower bandgaps, represents an important target for semiconducting halide perovskites, but has so far proven inaccessible using traditional spacer cations due to either intrinsic instability or formation of competing non-perovskite phases. Here, we demonstrate the first example of a 2D Ag-Bi iodide DP with a direct bandgap of 2.00(2) eV, templated by a layer of bifunctionalized oligothiophene cations, i.e., (bis-aminoethyl)bithiophene, through a collective influence of aromatic interactions, hydrogen bonding, bidentate tethering, and structural rigidity. Hybrid density functional theory calculations for the new material reveal a direct bandgap, consistent with the experimental value, and relatively flat band edges derived principally from Ag-d/I-p (valence band) and Bi-p/I-p (conduction band) states. This work opens up new avenues for exploring specifically designed organic cations to stabilize otherwise inaccessible 2D HOIPs with potential applications for optoelectronics.

Unified description of H-atom-induced chemicurrents and inelastic scattering.

The Born-Oppenheimer approximation (BOA) provides the foundation for virtually all computational studies of chemical binding and reactivity, and it is the justification for the widely used "balls and springs" picture of molecules. The BOA assumes that nuclei effectively stand still on the timescale of electronic motion, due to their large masses relative to electrons. This implies electrons never change their energy quantum state. When molecules react, atoms must move, meaning that electrons may become excited in violation of the BOA. Such electronic excitation is clearly seen for: (i) Schottky diodes where H adsorption at Ag surfaces produces electrical "chemicurrent;" (ii) Au-based metal-insulator-metal (MIM) devices, where chemicurrents arise from H-H surface recombination; and (iii) Inelastic energy transfer, where H collisions with Au surfaces show H-atom translation excites the metal's electrons. As part of this work, we report isotopically selective hydrogen/deuterium (H/D) translational inelasticity measurements in collisions with Ag and Au. Together, these experiments provide an opportunity to test new theories that simultaneously describe both nuclear and electronic motion, a standing challenge to the field. Here, we show results of a recently developed first-principles theory that quantitatively explains both inelastic scattering experiments that probe nuclear motion and chemicurrent experiments that probe electronic excitation. The theory explains the magnitude of chemicurrents on Ag Schottky diodes and resolves an apparent paradox--chemicurrents exhibit a much larger isotope effect than does H/D inelastic scattering. It also explains why, unlike Ag-based Schottky diodes, Au-based MIM devices are insensitive to H adsorption.

Hydrogen collisions with transition metal surfaces: Universal electronically nonadiabatic adsorption.

Inelastic scattering of H and D atoms from the (111) surfaces of six fcc transition metals (Au, Pt, Ag, Pd, Cu, and Ni) was investigated, and in each case, excitation of electron-hole pairs dominates the inelasticity. The results are very similar for all six metals. Differences in the average kinetic energy losses between metals can mainly be attributed to different efficiencies in the coupling to phonons due to the different masses of the metal atoms. The experimental observations can be reproduced by molecular dynamics simulations based on full-dimensional potential energy surfaces and including electronic excitations by using electronic friction in the local density friction approximation. The determining factors for the energy loss are the electron density at the surface, which is similar for all six metals, and the mass ratio between the impinging atoms and the surface atoms. Details of the electronic structure of the metal do not play a significant role. The experimentally validated simulations are used to explore sticking over a wide range of incidence conditions. We find that the sticking probability increases for H and D collisions near normal incidence-consistent with a previously reported penetration-resurfacing mechanism. The sticking probability for H or D on any of these metals may be represented as a simple function of the incidence energy, Ein, metal atom mass, M, and incidence angle, 𝜗in. S=(S0+a⋅Ein+b⋅M)*(1-h(𝜗in-c)(1-cos(𝜗in-c)d⋅h(Ein-e)(Ein-e))), where h is the Heaviside step function and for H, S0 = 1.081, a = -0.125 eV-1, b=-8.40⋅10-4 u-1, c = 28.88°, d = 1.166 eV-1, and e = 0.442 eV; whereas for D, S0 = 1.120, a = -0.124 eV-1, b=-1.20⋅10-3 u-1, c = 28.62°, d = 1.196 eV-1, and e = 0.474 eV.

Electron-hole pair excitation determines the mechanism of hydrogen atom adsorption.

How much translational energy atoms and molecules lose in collisions at surfaces determines whether they adsorb or scatter. The fact that hydrogen (H) atoms stick to metal surfaces poses a basic question. Momentum and energy conservation demands that the light H atom cannot efficiently transfer its energy to the heavier atoms of the solid in a binary collision. How then do H atoms efficiently stick to metal surfaces? We show through experiments that H-atom collisions at an insulating surface (an adsorbed xenon layer on a gold single-crystal surface) are indeed nearly elastic, following the predictions of energy and momentum conservation. In contrast, H-atom collisions with the bare gold surface exhibit a large loss of translational energy that can be reproduced by an atomic-level simulation describing electron-hole pair excitation.

An accurate full-dimensional potential energy surface for H-Au(111): Importance of nonadiabatic electronic excitation in energy transfer and adsorption.

We have constructed a potential energy surface (PES) for H-atoms interacting with fcc Au(111) based on fitting the analytic form of the energy from Effective Medium Theory (EMT) to ab initio energy values calculated with density functional theory. The fit used input from configurations of the H-Au system with Au atoms at their lattice positions as well as configurations with the Au atoms displaced from their lattice positions. It reproduces the energy, in full dimension, not only for the configurations used as input but also for a large number of additional configurations derived from ab initio molecular dynamics (AIMD) trajectories at finite temperature. Adiabatic molecular dynamics simulations on this PES reproduce the energy loss behavior of AIMD. EMT also provides expressions for the embedding electron density, which enabled us to develop a self-consistent approach to simulate nonadiabatic electron-hole pair excitation and their effect on the motion of the incident H-atoms. For H atoms with an energy of 2.7 eV colliding with Au, electron-hole pair excitation is by far the most important energy loss pathway, giving an average energy loss ≈3 times that of the adiabatic case. This increased energy loss enhances the probability of the H-atom remaining on or in the Au slab by a factor of 2. The most likely outcome for H-atoms that are not scattered also depends prodigiously on the energy transfer mechanism; for the nonadiabatic case, more than 50% of the H-atoms which do not scatter are adsorbed on the surface, while for the adiabatic case more than 50% pass entirely through the 4 layer simulation slab.