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-(+)- orS-(−)-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 betweenR- andS-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.

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 (C₆H₅(CH₂)₂NH₃)₂PbI₄ 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.

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 H₂ 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.

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.

For the prototypical two-dimensional hybrid organic–inorganic perovskites (2D HOIPs) (AE4T)PbX₄ (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.

The electronic properties of hybrid organic–inorganic semiconductor interfaces depend strongly on the alignment of the electronic carrier levels in the organic/inorganic components. In the present work, we address this energy level alignment from first principles theory for two paradigmatic organic–inorganic semiconductor interfaces, the singlet fission materials tetracene and pentacene on H/Si(111), using all-electron density functional theory calculations with a hybrid exchange–correlation functional. For isolated tetracene on H/Si(111), a type I-like heterojunction (lowest-energy electron and hole states on Si) is found. For isolated pentacene, the molecular and semiconductor valence band edges are degenerate. For monolayer films, we show how to construct supercell geometries with up to 1192 atoms, which minimize the strain between the inorganic surface and an organic monolayer film. Based on these models, we predict the formation of type II heterojunctions (electron states on Si, hole-like states on the organic species) for both acenes, indicating that charge separation at the interface between the organic and inorganic components is favored. The paper discusses the steps needed to find appropriate low-energy interface geometries for weakly bonded organic molecules and films on inorganic substrates from first principles, a necessary prerequisite for any computational level alignment prediction.