Nonadiabatic charge dynamics in novel solar cell materials

An Exhaustive Quantum-Classical Study of C6F6+ Using the Newly Formulated Parallel TDDVR Method

We recently implemented our parallelized quantum-classical dynamical approach, known as the Time-Dependent Discrete Variable Representation (TDDVR) method, which is applied to the spectroscopically important hexafluorobenzene (HFBz) radical cation, where several conical intersections exist in their seven lowest excited electronic states (S11B2u, S21E1g, S31B1u, S41E1u, and S51A2u) considering degeneracy among potential energy surfaces (PESs), to demonstrate their various dynamical aspects. This new parallel version shows almost linear scalability with increasing number of computing processors. To get photoelectron (PE) spectra, Mass-Analyzed Threshold Ionization (MATI) spectra, population dynamics, and many other dynamical observables, the first-principles dynamics is applied at the state-of-the-art level to the corresponding Hamiltonian, where the Jahn-Teller (JT) and pseudo-Jahn-Teller (PJT) type interactions are involved in those coupled seven electronic states. The quantum-classical method is used to thoroughly analyze the effects of these couplings on the nuclear dynamics of the involved electronic states, and the findings are compared with those observables obtained from experiments. Intrinsic dynamical properties are explained using the reduced densities of the wave packet (WP) in a coupled electronic manifold. The PE and MATI spectra of HFBz computed using TDDVR are found to be in good agreement with earlier experimental data and other theoretically simulated spectra.

Compression of Organic Molecules Coupled with Hydrogen Bonding Extends the Charge Carrier Lifetime in BA2SnI4

Two-dimensional (2D) metal halide perovskites, such as BA2SnI4 (BA═CH3(CH2)3NH3), exhibit an enhanced charge carrier lifetime in experiments under strain. Experiments suggest that significant compression of the BA molecule, rather than of the inorganic lattice, contributes to this enhancement. To elucidate the underlying physical mechanism, we apply a moderate compressive strain to the entire system and subsequently introduce significant compression to the BA molecules. We then perform ab initio nonadiabatic molecular dynamics simulations of nonradiative electron-hole recombination. We observe that the overall lattice compression reduces atomic motions and decreases nonadiabatic coupling, thereby delaying electron-hole recombination. Additionally, compression of the BA molecules enhances hydrogen bonding between the BA molecules and iodine atoms, which lengthens the Sn-I bonds, distorts the [SnI6]4- octahedra, and suppresses atomic motions further, thus reducing nonadiabatic coupling. Also, the elongated Sn-I bonds and weakened antibonding interactions increase the band gap. Altogether, the compression delays the nonradiative electron-hole recombination by more than a factor of 3. Our simulations provide new and valuable physical insights into how compressive strain, accommodated primarily by the organic ligands, positively influences the optoelectronic properties of 2D layered halide perovskites, offering a promising pathway for further performance improvements.

Water Solvent Reorganization upon Ultrafast Resonant Stimulated X-ray Raman Excitation of a Metalloporphyrin Dimer

We propose an X-ray Raman pump-X-ray diffraction probe scheme to follow solvation dynamics upon charge migration in a solute molecule. The X-ray Raman pump selectively prepares a valence electronic wavepacket in the solute, while the probe provides information about the entire molecular ensemble. A combination of molecular dynamics and ab initio quantum chemistry simulations is applied to a Zn-Ni porphyrin dimer in water. Using time-resolved X-ray diffraction and pair distribution functions, we extracted solvation shell dynamics.

Nonadiabatic Field with Triangle Window Functions on Quantum Phase Space

Recent progress on the constraint coordinate-momentum phase space (CPS) formulation of finite-state quantum systems has revealed that the triangle window function approach is an isomorphic representation of the exact population-population correlation function of the two-state system. We use the triangle window (TW) function and the CPS mapping kernel element to formulate a novel useful representation of discrete electronic degrees of freedom (DOFs). When it is employed with nonadiabatic field (NaF) dynamics, a new variant of the NaF approach (i.e., NaF-TW) is proposed. The NaF-TW expression of the population of any adiabatic state is always positive semidefinite. Extensive benchmark tests of model systems in both the condensed phase and gas phase demonstrate that the NaF-TW approach is able to faithfully capture the dynamical interplay between electronic and nuclear DOFs in a broad region, including where the states remain coupled all the time, as well as where the bifurcation characteristic of nuclear motion is important.

Nonadiabatic Field on Quantum Phase Space: A Century after Ehrenfest

Nonadiabatic transition dynamics lies at the core of many electron/hole transfer, photoactivated, and vacuum field-coupled processes. About a century after Ehrenfest proposed "Phasenraum" and the Ehrenfest theorem, we report a conceptually novel trajectory-based nonadiabatic dynamics approach, nonadiabatic field (NAF), based on a generalized exact coordinate-momentum phase space formulation of quantum mechanics. It does not employ the conventional Born-Oppenheimer or Ehrenfest trajectory in the nonadiabatic coupling region. Instead, in NAF the equations of motion of the independent trajectory involve a nonadiabatic nuclear force term in addition to an adiabatic nuclear force term of a single electronic state. A few benchmark tests for gas phase and condensed phase systems indicate that NAF offers a practical tool to capture the correct correlation of electronic and nuclear dynamics for processes where the states remain coupled all the time as well as for the asymptotic region where the coupling of electronic states vanishes.

Nonadiabatic Dynamics Simulations for Photoinduced Processes in Molecules and Semiconductors: Methodologies and Applications

Nonadiabatic dynamics (NAMD) simulations have become powerful tools for elucidating complicated photoinduced processes in various systems from molecules to semiconductor materials. In this review, we present an overview of our recent research on photophysics of molecular systems and periodic semiconductor materials with the aid of ab initio NAMD simulation methods implemented in the generalized trajectory surface-hopping (GTSH) package. Both theoretical backgrounds and applications of the developed NAMD methods are presented in detail. For molecular systems, the linear-response time-dependent density functional theory (LR-TDDFT) method is primarily used to model electronic structures in NAMD simulations owing to its balanced efficiency and accuracy. Moreover, the efficient algorithms for calculating nonadiabatic coupling terms (NACTs) and spin-orbit couplings (SOCs) have been coded into the package to increase the simulation efficiency. In combination with various analysis techniques, we can explore the mechanistic details of the photoinduced dynamics of a range of molecular systems, including charge separation and energy transfer processes in organic donor-acceptor structures, ultrafast intersystem crossing (ISC) processes in transition metal complexes (TMCs), and exciton dynamics in molecular aggregates. For semiconductor materials, we developed the NAMD methods for simulating the photoinduced carrier dynamics within the framework of the Kohn-Sham density functional theory (KS-DFT), in which SOC effects are explicitly accounted for using the two-component, noncollinear DFT method. Using this method, we have investigated the photoinduced carrier dynamics at the interface of a variety of van der Waals (vdW) heterojunctions, such as two-dimensional transition metal dichalcogenides (TMDs), carbon nanotubes (CNTs), and perovskites-related systems. Recently, we extended the LR-TDDFT-based NAMD method for semiconductor materials, allowing us to study the excitonic effects in the photoinduced energy transfer process. These results demonstrate that the NAMD simulations are powerful tools for exploring the photodynamics of molecular systems and semiconductor materials. In future studies, the NAMD simulation methods can be employed to elucidate experimental phenomena and reveal microscopic details as well as rationally design novel photofunctional materials with desired properties.

What is Missing in the Mean Field Description of Spatial Distribution of Population? Important Role of Auxiliary Wave Packets in Trajectory Branching

When the traditional Ehrenfest mean field approach is employed to simulate nonadiabatic dynamics, an effective wave packet (WP) on the average potential energy surface (PES) is utilized to describe the nuclear motion. In the fully quantum picture, however, the WP components on different adiabatic PESs gradually separate in space because they evolve under different velocities and forces. Due to trajectory branching of the WP components, proper decoherence needs to be taken into account, and the spatial distribution of population cannot be described by a single effective WP. Here, we propose an auxiliary branching corrected mean field (A-BCMF) method, where trajectories of auxiliary WPs on adiabatic PESs are introduced. As benchmarked in the three standard Tully models, A-BCMF not only gives correct channel populations but also captures an accurate time-dependent spatial distribution of population. Thereby, we reveal the important role of auxiliary WPs in solving intrinsic problems of the widely used mean field description of nonadiabatic dynamics.

Fantastical excited state optimized structures and where to find them

The quantum chemistry community has developed analytic forces for approximate electronic excited states to enable walking on excited state potential energy surfaces (PES). One can thereby computationally characterize excited state minima and saddle points. Always implicit in using this machinery is the fact that an excited state PES only exists within the realm of the Born-Oppenheimer approximation, where the nuclear and electronic degrees of freedom separate. This work demonstrates through ab initio calculations and simple nonadiabatic dynamics that some excited state minimum structures are fantastical: they appear to exist as stable configurations only as a consequence of the PES construct, rather than being physically observable. Each fantastical structure exhibits an unphysically high predicted harmonic frequency and associated force constant. This fact can serve as a valuable diagnostic of when an optimized excited state structure is non-observable. The origin of this phenomenon can be attributed to the coupling between different electronic states. As PESs approach one another, the upper surface can form a minimum that is very close to a near-touching point. The force constant, evaluated at this minimum, relates to the strength of the electronic coupling rather than to any characteristic excited state vibration. Nonadiabatic dynamics results using a Landau-Zener model illustrate that fantastical excited state structures have extremely short lifetimes on the order of a few femtoseconds. Their appearance in a calculation signals the presence of a nearby conical intersection through which the system will rapidly cross to a lower surface.

Charge Recombination Dynamics in a Metal Halide Perovskite Simulated by Nonadiabatic Molecular Dynamics Combined with Machine Learning

Nonadiabatic coupling (NAC) plays a central role in driving nonadiabatic dynamics in various photophysical and photochemical processes. However, the high computational cost of NAC limits the time scale and system size of quantum dynamics simulation. By developing a machine learning (ML) framework and applying it to a traditional CH₃N₃PbI₃ perovskite, we demonstrate that the various ML algorithms (XGBoost, LightGBM, and random forest) combined with three descriptors (sine matrix, MBTR, and SOAP) can predict accurate NACs that all agree well with the direct calculations, particularly for the combination of LightGBM and sine matrix descriptor showing the best performance with a high correlation coefficient of ≤0.87. The simulated nonradiative electron-hole recombination time scales agree well with each other between the NACs obtained from direct calculations and ML prediction. The study shows the advantage in accelerating quantum dynamics simulations using ML algorithms.

Great Influence of Organic Cation Motion on Charge Carrier Dynamics in Metal Halide Perovskite Unraveled by Unsupervised Machine Learning

Unsupervised machine learning combined with time-dependent density functional theory reveals the significant influence of organic cation on the charge carrier lifetime of FAPbI₃ (FA = HC(NH₂)₂⁺) by analyzing their mutual information (MI) between the geometric features and the nonadiabatic coupling (NAC) and bandgap. Analysis of MI values demonstrates that the NAC and bandgap are dominated by the orientation and shape of the inorganic octahedron because iodine and lead atoms are composed of the band edge states. Counterintuitively, the correlated motion promotes the contribution of the FA cation to the NAC; in particular, one type of FA rotation even supersedes the influence of the velocities of the lead and iodine atoms due to the enhanced hydrogen bond interaction. Our study demonstrates the importance of the correlated motion on the excited-state lifetimes of FAPbI₃, which provides a guidance for optimizing the optoelectronic properties of metal halide perovskites.

Transition Path Flight Times and Nonadiabatic Electronic Transitions

Transition path flight times are studied for scattering on two electronic surfaces with a single crossing. These flight times reveal nontrivial quantum effects such as resonance lifetimes and nonclassical passage times and reveal that nonadiabatic effects often increase flight times. The flight times are computed using numerically exact time propagation and compared with results obtained from the Fewest Switches Surface Hopping (FSSH) method. Comparison of the two methods shows that the FSSH method is reliable for transition path times only when the scattering is classically allowed on the relevant adiabatic surfaces. However, where quantum effects such as tunneling and resonances dominate, the FSSH method is not adequate to accurately predict the correct times and transition probabilities. These results highlight limitations in methods which do not account for quantum interference effects, and suggest that measuring flight times is important for obtaining insights from the time-domain into quantum effects in nonadiabatic scattering.

Electric Field Effects on Photoelectrochemical Water Splitting: Perspectives and Outlook

The grand challenges in renewable energy lie in our ability to comprehend efficient energy conversion systems, together with dealing with the problem of intermittency via scalable energy storage systems. Relatively little progress has been made on this at grid scale and two overriding challenges still need to be addressed: (i) limiting damage to the environment and (ii) the question of environmentally friendly energy conversion. The present review focuses on a novel route for producing hydrogen, the ultimate clean fuel, from the Sun, and renewable energy source. Hydrogen can be produced by light-driven photoelectrochemical (PEC) water splitting, but it is very inefficient; rather, we focus here on how electric fields can be applied to metal oxide/water systems in tailoring the interplay with their intrinsic electric fields, and in how this can alter and boost PEC activity, drawing both on experiment and non-equilibrium molecular simulation.

Extending the Time Scales of Nonadiabatic Molecular Dynamics via Machine Learning in the Time Domain

A novel methodology for direct modeling of long-time scale nonadiabatic dynamics in extended nanoscale and solid-state systems is developed. The presented approach enables forecasting the vibronic Hamiltonians as a direct function of time via machine-learning models trained directly in the time domain. The use of periodic and aperiodic functions that transform time into effective input modes of the artificial neural network is demonstrated to be essential for such an approach to work for both abstract and atomistic models. The best strategies and possible limitations pertaining to the new methodology are explored and discussed. An exemplary direct simulation of unprecedentedly long 20 picosecond trajectories is conducted for a divacancy-containing monolayer black phosphorus system, and the importance of conducting such extended simulations is demonstrated. New insights into the excited states photophysics in this system are presented, including the role of decoherence and model definition.

Significance of the Chemical Environment of an Element in Nonadiabatic Molecular Dynamics: Feature Selection and Dimensionality Reduction with Machine Learning

Using supervised and unsupervised machine learning (ML) on features generated from nonadiabatic (NA) molecular dynamics (MD) trajectories under the classical path approximation, we demonstrate that mutual information with the NA Hamiltonian can be used for feature selection and model simplification. Focusing on CsPbI₃, a popular metal halide perovskite, we observe that the chemical environment of a single element is sufficient for predicting the NA Hamiltonian. The conclusion applies even to Cs, although Cs does not contribute to the relevant wave functions. Interatomic distances between Cs and I or Pb and the octahedral tilt angle are the most important features. We reduce a typical 360-parameter ML force-field model to just a 12-parameter NA Hamiltonian model, while maintaining a high NA-MD simulation quality. Because NA-MD is a valuable tool for studying excited state processes, overcoming its high computational cost through simple ML models will streamline NA-MD simulations and expand the ranges of accessible system size and simulation time.

Increasing Efficiency of Nonadiabatic Molecular Dynamics by Hamiltonian Interpolation with Kernel Ridge Regression

Nonadiabatic (NA) molecular dynamics (MD) goes beyond the adiabatic Born-Oppenheimer approximation to account for transitions between electronic states. Such processes are common in molecules and materials used in solar energy, optoelectronics, sensing, and many other fields. NA-MD simulations are much more expensive compared to adiabatic MD due to the need to compute excited state properties and NA couplings (NACs). Similarly, application of machine learning (ML) to NA-MD is more challenging compared with adiabatic MD. We develop an NA-MD simulation strategy in which an adiabatic MD trajectory, which can be generated with a ML force field, is used to sample excitation energies and NACs for a small fraction of geometries, while the properties for the remaining geometries are interpolated with kernel ridge regression (KRR). This ML strategy allows for one to perform NA-MD under the classical path approximation, increasing the computational efficiency by over an order of magnitude. Compared to neural networks, KRR requires little parameter tuning, saving efforts on model building. The developed strategy is demonstrated with two metal halide perovskites that exhibit complicated MD and are actively studied for various applications.

Reduced-dimensional surface hopping with offline-online computations

Molecular dynamics simulations often classically evolve the nuclear geometry on adiabatic potential energy surfaces (PESs), punctuated by random hops between energy levels in regions of strong coupling, in an algorithm known as surface hopping. However, the computational expense of integrating the geometry on a full-dimensional PES and computing the required couplings can quickly become prohibitive as the number of atoms increases. In this work, we describe a method for surface hopping that uses only important reaction coordinates, performs all expensive evaluations of the true PESs and couplings only once before simulating dynamics (offline), and then queries the stored values during the surface hopping simulation (online). Our Python codes are freely available on GitHub. Using photodissociation of azomethane as a test case, this method is able to reproduce experimental results that have thus far eluded ab initio surface hopping studies.

Dependence between Structural and Electronic Properties of CsPbI3: Unsupervised Machine Learning of Nonadiabatic Molecular Dynamics

Using unsupervised machine learning on the trajectories from a nonadiabatic molecular dynamics simulation with time-dependent Kohn-Sham density functional theory, we elucidated the structural parameters with the largest influence on nonradiative recombination of charge carriers in CsPbI₃, which forms the basis for solar energy and optoelectronic applications. The I-I-I angles between PbI₆ octahedra, followed by the Cs-I distance, have the strongest impact on the bandgap and the nonadiabatic coupling. The importance of the Cs-I distance is unexpected, because Cs does not contribute to electron and hole wave functions. The nonadiabatic coupling is most influenced by static properties, which is also surprising, given its explicit dependence on atomic velocities.

Commutator Matrix in Phase Space Mapping Models for Nonadiabatic Quantum Dynamics

We show that a novel, general phase space mapping Hamiltonian for nonadiabatic systems, which is reminiscent of the renowned Meyer-Miller mapping Hamiltonian, involves a commutator variable matrix rather than the conventional zero-point-energy parameter. In the exact mapping formulation on constraint space for phase space approaches for nonadiabatic dynamics, the general mapping Hamiltonian with commutator variables can be employed to generate approximate trajectory-based dynamics. Various benchmark model tests, which range from gas phase to condensed phase systems, suggest that the overall performance of the general mapping Hamiltonian is better than that of the conventional Meyer-Miller Hamiltonian.

Branching and phase corrected surface hopping: A benchmark of nonadiabatic dynamics in multilevel systems

Since the seminal work of Tully [J. Chem. Phys. 93, 1061 (1990)], two-level scattering models have been extensively adopted as the standard benchmark systems to assess the performance of different trajectory surface hopping methods for nonadiabatic dynamics simulations. Here, we extend the branching and phase corrections to multilevel systems and combine them with both the traditional fewest switches surface hopping (FSSH) and its variant global flux surface hopping (GFSH) algorithms. To get a comprehensive evaluation of the proposed methods, we construct a series of more challenging and diverse three-level and four-level scattering models and use exact quantum solutions as references. Encouragingly, both FSSH and GFSH with the branching and phase corrections produce excellent and nearly identical results in all investigated systems, indicating that the new surface hopping methods are robust to describe multilevel problems and the reliability is insensitive to the definition of self-consistent hopping probabilities in the adiabatic representation. Furthermore, the branching correction is found to be especially important when dealing with strongly repulsive potential energy surfaces, which are common in realistic systems, thus promising for general applications.

Materials Design and Optimization for Next-Generation Solar Cell and Light-Emitting Technologies

We review some of the most potent directions in the design of materials for next-generation solar cell and light-emitting technologies that go beyond traditional solid-state inorganic semiconductor-based devices, from both the experimental and computational standpoints. We focus on selected recent conceptual advances in tackling issues which are expected to significantly impact applied literature in the coming years. Specifically, we consider solution processability, design of dopant-free charge transport materials, two-dimensional conjugated polymeric semiconductors, and colloidal quantum dot assemblies in the fields of experimental synthesis, characterization, and device fabrication. Key modeling issues that we consider are calculations of optical properties and of effects of aggregation, including recent advances in methods beyond linear-response time-dependent density functional theory and recent insights into the effects of correlation when going beyond the single-particle ansatz as well as in the context of modeling of thermally activated fluorescence.

Concentric Approximation for Fast and Accurate Numerical Evaluation of Nonadiabatic Coupling with Projector Augmented-Wave Pseudopotentials

We develop an efficient and accurate method for numerical evaluation of nonadiabatic (NA) coupling in the Kohn-Sham representation with projector augmented-wave (PAW) pseudopotentials that are commonly used in electronic structure calculations on nanoscale, condensed matter, and molecular systems. Without additional cost, the method provides an order of magnitude improvement in accuracy compared to the current technique, while it is 3-4 orders of magnitude faster than the exact evaluation. Atomic displacements over typical time steps in molecular dynamics (MD) simulations are much smaller than the size of the PAW core region, and therefore, evaluation of the NA in the core is simplified. The accuracy is demonstrated with three condensed matter systems. The method is robust to variation in the MD time step. The accurate NA coupling evaluation also helps in maintaining phase-consistency of the NA coupling and identifying trivial crossings of adiabatic states. The approach stimulates NAMD applications to modeling of modern materials and processes.

Negative Zero-Point-Energy Parameter in the Meyer-Miller Mapping Model for Nonadiabatic Dynamics

The celebrated Meyer-Miller mapping model has been a useful approach for generating practical trajectory-based nonadiabatic dynamics methods. It is generally assumed that the zero-point-energy (ZPE) parameter is positive. The constraint implied in the conventional Meyer-Miller mapping Hamiltonian for an F-electronic-state system actually requires γ∈(-1/F, ∞) for the ZPE parameter for each electronic degree of freedom. Both negative and positive values are possible for such a parameter. We first establish a rigorous formulation to construct exact mapping models in the Cartesian phase space when the constraint is applied. When nuclear dynamics is approximated by the linearized semiclassical initial value representation, a negative ZPE parameter could lead to reasonably good performance in describing dynamic behaviors in typical spin-boson models for condensed-phase two-state systems, even at challenging zero temperature.

Phonon-Mediated Interlayer Charge Separation and Recombination in a MoSe2/WSe2 Heterostructure

Monolayer transition metal dichalcogenides bear great potential for photodetection and light harvesting due to high absorption coefficients. However, these applications require dissociation of strongly bound photogenerated excitons. The dissociation can be achieved by vertically stacking different monolayers to realize band alignment that favors interlayer charge transfer. In such heterostructures, the reported recombination times vary strongly, and the charge separation and recombination mechanisms remain elusive. We use two color pump-probe microscopy to demonstrate that the charge separation in a MoSe₂/WSe₂ heterostructure is ultrafast (∼200 fs) and virtually temperature independent, whereas the recombination accelerates strongly with temperature. Ab initio quantum dynamics simulations rationalize the experiments, indicating that the charge separation is temperature-independent because it is barrierless, involves dense acceptor states, and is promoted by higher-frequency out-of-plane vibrations. The strong temperature dependence of the recombination, on the other hand, arises from a transient indirect-to-direct bandgap modulation by low-frequency shear and layer breathing motions.

Accurate Computation of Nonadiabatic Coupling with Projector Augmented-Wave Pseudopotentials

Synergy of nonadiabatic molecular dynamics with real-time time-dependent density functional theory has led to significant progress in modeling excited-state dynamics in nanoscale and condensed matter systems over the past decade. Nonadiabatic coupling (NAC) is the central quantity in such simulations, and its accurate and efficient evaluation is an enduring challenge in time-dependent Kohn-Sham theory, particularly in conjunction with planewave basis sets and projector augmented-wave (PAW) pseudopotentials because of the complexity of the PAW "all-electron" wave function. We report a method for rigorous evaluation of the NAC with PAW wave functions and demonstrate an efficient approximation to the rigorous NAC that gives comparable accuracy. As a validation, we intensely examine the NAC matrix elements calculated using both pseudo- and all-electron wave functions under the PAW formalism in six representative systems. The approximate NAC obtained with pseudowave functions is close to the exact all-electron NAC, with the largest deviations observed when subshell d-electrons are involved in the transitions. The developed approach provides a rigorous and convenient methodology for the numerical computation of NAC in the Kohn-Sham theory framework.

Branching Corrected Mean Field Method for Nonadiabatic Dynamics

When nonadiabatic dynamics are described on the basis of trajectories, severe trajectory branching occurs when the nuclear wave packets on some potential energy surfaces are reflected while those on the remaining surfaces are not. As a result, the traditional Ehrenfest mean field (EMF) approximation breaks down. In this study, two versions of the branching corrected mean field (BCMF) method are proposed. Namely, when trajectory branching is identified, BCMF stochastically selects either the reflected or the nonreflected group to build the new mean field trajectory or splits the mean field trajectory into two new trajectories with the corresponding weights. As benchmarked in six standard model systems and an extensive model base with two hundred diverse scattering models, BCMF significantly improves the accuracy while retaining the high efficiency of the traditional EMF. In fact, BCMF closely reproduces the exact quantum dynamics in all investigated systems, thus highlighting the essential role of branching correction in nonadiabatic dynamics simulations of general systems.

Long-Lived Hot Electron in a Metallic Particle for Plasmonics and Catalysis: Ab Initio Nonadiabatic Molecular Dynamics with Machine Learning

Multiple experiments provide evidence for photovoltaic, catalytic, optoelectronic, and plasmonic processes involving hot, i.e., high energy, electrons in nanoscale materials. However, the mechanisms of such processes remain elusive, because electrons rapidly lose energy by relaxation through dense manifolds of states. We demonstrate a long-lived hot electron state in a Pt nanocluster adsorbed on the MoS₂ substrate. For this purpose, we develop a simulation technique, combining classical molecular dynamics based on machine learning potentials with ab initio nonadiabatic molecular dynamics and real-time time-dependent density functional theory. Choosing Pt₂₀/MoS₂ as a prototypical system, we find frequent shifting of a top atom in the Pt particle occurring on a 50 ps time scale. The distortion breaks particle symmetry and creates unsaturated chemical bonds. The lifetime of the localized state associated with the broken bonds is enhanced by a factor of 3. Hot electrons aggregate near the shifted atom and form a catalytic reaction center. Our findings prove that distortion of even a single atom can have important implications for nanoscale catalysis and plasmonics and provide insights for utilizing machine learning potentials to accelerate ab initio investigations of excited state dynamics in condensed matter systems.

Optoelectronic Properties of Two-Dimensional Bromide Perovskites: Influences of Spacer Cations

Two-dimensional (2D) halide perovskites have displayed unique emission properties, making them potential candidates for next-generation light-emitting devices. Here, we combine nonadiabatic molecular dynamics and time-domain density functional theory to investigate the fundamental mechanisms of carrier recombination processes. Considering monolayer bromide perovskites with dissimilar organic spacer molecules, n-butylammonium (BA) and phenylethylammonium (PEA) cations, we find a strong correlation between temperature-induced structural fluctuations and nonradiative carrier recombination rates in these materials. The more flexible geometry of (BA)₂PbBr₄ compared to that of (PEA)₂PbBr₄, results in faster electron-hole recombination and shorter carrier lifetime, diminishing the photoluminescence quantum yield for softer 2D perovskites. Reduced structural fluctuations in relatively rigid (PEA)₂PbBr₄ not only indicate of a longer carrier lifetime but also suggest a narrower emission line width, implying a higher purity of the emitted light. Our ab initio modeling of excited state properties in 2D perovskites conveys material designing strategies to fine-tune perovskite emissions for solid-state lighting applications.

Control of Charge Carrier Dynamics in Plasmonic Au Films by TiOx Substrate Stoichiometry

Plasmonic excitations in noble metals have many fascinating properties and give rise to a broad range of applications. We demonstrate, using nonadiabatic molecular dynamics combined with time-domain density functional theory, that the chemical composition and stoichiometry of substrates can have a strong influence on charge dynamics. By changing oxygen content in TiO₂, including stoichiometric, oxygen rich, and oxygen poor phases, and Ti metal, one can alter lifetimes of charge carriers in Au by a factor of 5 and control the ratio of electron-to-hole relaxation rates by a factor of 10. Remarkably, a thin TiO_x substrate greatly alters charge carrier properties in much thicker Au films. Such large variations stem from the fact that the Ti and O atoms are much lighter than Au, and their vibrations are much faster at dissipating the energy. The control over a particular charge carrier and an energy range depends on the Au and TiO_x level alignment, and the interfacial interaction strength. These factors are easily influenced by the TiO_x stoichiometry. In particular, oxygen rich and poor TiO₂ can be used to control holes and electrons, respectively, while metallic Ti affects both charge carriers. The detailed atomistic analysis of the interfacial and electron-vibrational interactions generates the fundamental understanding of the properties of plasmonic materials needed to design photovoltaic, photocatalytic, optoelectronic, sensing, nanomedical, and other devices.

Strain Controls Charge Carrier Lifetimes in Monolayer WSe2: Ab Initio Time Domain Analysis

Mono- and few-layer transition metal dichalcogenides (TMDs) are among the most appealing candidates for electronic and optoelectronic devices. During synthesis, TMDs actively interact with substrates, which induce notable strain and influence significantly charge carriers in TMDs. By performing time-domain ab initio simulations on monolayer WSe₂, we demonstrate that charge carrier lifetimes vary by a factor of 3 within a typical 1% strain range, the bandgap changes by 0.2 eV, and electron-phonon interactions vary by 60%. Fortuitously, the most common tensile strain extends the lifetimes. The changes arise because of modifications in interatomic interactions. Further, compared to the optimized structure, at ambient temperature the bandgap drops by 0.1 eV and fluctuates by 0.1 eV. WSe₂ obeys linear response within 1% strain; however, further strain leads to nonlinear qualitative changes in WSe₂ electronic properties. The conduction band is affected more strongly than the valence band. Charges couple to phonons within a 100-400 cm^-1 frequency range, with the strongest coupling to in-plane and out-of-plane modes at 250 cm^-1. The reported findings agree with the available experiments and should be generic to other 2D materials. The strain effects need to be considered during TMD synthesis and provide means to control and tune TMD properties for 2D device applications.

Lattice Expansion in Hybrid Perovskites: Effect on Optoelectronic Properties and Charge Carrier Dynamics

Hybrid halide perovskites frequently undergo structural expansion due to various stimuli, significantly affecting their electronic properties and in particular their charge carrier dynamics. It is essential to atomistically model how geometric changes modify electronic characteristics that are important for applications such as light harvesting and lighting. Using ab initio simulations, here we investigate the structural dynamics and optoelectronic properties of FAPbI₃ under tensile strain. The applied strain leads to elongation of the Pb-I bonds and a decrease in the level of PbI₆ octahedral tilting, which manifests as blue-shifts in band gaps. Nonadiabatic molecular dynamics simulations further reveal that charge carrier recombination rates moderately decrease in these expanded lattices. The complex influence of lattice dynamics on electron-phonon scattering results in a longer carrier lifetime, which is advantageous for efficient solar cells. By providing detailed information about the structure-property relationships, this work emphasizes the role of controlled lattice expansion in enhancing the electronic functionalities of hybrid perovskites.

Influence of Defects on Excited-State Dynamics in Lead Halide Perovskites: Time-Domain ab Initio Studies

This Perspective summarizes recent research into the excited-state dynamics in lead halide perovskites that are of paramount importance for photovoltaic and photocatalytic applications. Nonadiabatic molecular dynamics combined with time-domain ab initio density functional theory allows one to mimic time-resolved spectroscopy experiments at the atomistic level of detail. The focus is placed on realistic aspects of perovskite materials, including point defects, surfaces, grain boundaries, mixed stoichiometries, dopants, and interfaces. The atomistic description of the quantum dynamics of electron and hole trapping and recombination, provided by the time-domain ab initio simulations, generates important insights into the mechanisms of charge and energy losses and guides the development of high-performance perovskite solar cell devices.

Understanding structures and properties of phosphorene/perovskite heterojunction toward perovskite solar cell applications

Two-dimensional black phosphorus (phosphorene) has drawn much attention in recent years due to its excellent electronic and optical properties. In this manuscript, we employ ab initio calculations to investigate the structural origin of the phosphorene/perovskite heterostructure. The calculations suggest that the chemical stability and the mechanical stability depend on the surface terminations, and the mechanical stability of the phosphorene/perovskite heterojunction should be further improved. The weak interactions between the P atoms in the phosphorene and the under-coordinated Pb atoms at the perovskite surfaces, as well as the weak interfacial charge transfer characters, are proposed to be mainly responsible for the moderate heterostructure stability. Suggestions to improve the stability of the heterojunction are provided. This study helps the fundamental understanding of the interaction between the phosphorene and the halide perovskite materials, and could provide a foundation for the better understanding of the low-dimensional materials in perovskite-based optoelectronic devices.

Subspace Surface Hopping with Size-Independent Dynamics

We present a subspace surface hopping strategy to deal with complex surface crossings in nonadiabatic dynamics. By focusing on only important adiabatic states, we make subspace crossing correction (SCC) in the framework of the standard fewest switches surface hopping (FSSH) and the global flux surface hopping (GFSH). The resulting SCC-FSSH and SCC-GFSH approaches show much better performance than the counterparts using all adiabatic states for surface hopping. As demonstrated in a series of Holstein models with up to over 1000 molecular sites, both SCC-FSSH and SCC-GFSH show excellent size independence with a large time step size of 1 fs. Especially, SCC-GFSH does not refer to nonadiabatic couplings at all and gives a more proper description of superexchange, and thus, it is promising for realistic applications with complex potential energy surfaces.

Grain Boundary Facilitates Photocatalytic Reaction in Rutile TiO2 Despite Fast Charge Recombination: A Time-Domain ab Initio Analysis

TiO₂ is an excellent photocatalytic and photovoltaic material but suffers low efficiency because of deep trap states giving rise to fast charge and energy losses. Using a combination of time-domain density functional theory and nonadiabatic molecular dynamics, we demonstrate that grain boundaries (GBs), which are common in polycrystalline TiO₂, accelerate nonradiative electron-hole recombination by a factor of 3. Despite GBs increase the band gap without creating deep trap states, and accelerate coherence loss, they enhance nonadiabatic electron-phonon coupling, and facilitate the relaxation. Importantly, electrons accumulated at the boundaries together with the relatively long-lived excite state favor photocatalytic reaction. Our study rationalizes the experimental observations and provides valuable perspectives for improving the device performance by defect engineering.

Crossing Classified and Corrected Fewest Switches Surface Hopping

In the traditional fewest switches surface hopping (FSSH), trivial crossings between uncoupled or weakly coupled states have highly peaked nonadiabatic couplings and thus are difficult to deal with in the preferred, adiabatic representation. Here, we classify surface crossings into four general types and propose a parameter-free crossing corrected FSSH (CC-FSSH) algorithm, which could treat multiple trivial crossings within a time interval. As examples, Holstein Hamiltonians with different parameters are adopted to mimic electron dynamics in tens to hundreds of molecules, which suffer from severe trivial crossing problems. Using existed surface hopping approaches as references, we show that CC-FSSH exhibits significantly fast time interval convergence and weak system size dependence. In all cases, a reliable description is achieved with a large time interval of 1 fs. With a simple formalism and the ability to describe complex surface crossings, CC-FSSH could potentially simulate general nonadiabatic dynamics in nanoscale materials with a high efficiency.

Unraveling the Mechanism of Photoinduced Charge-Transfer Process in Bilayer Heterojunction

Charge transfer is a fundamental process that determines the performance of solar cell devices. Although great efforts have been made, the detailed mechanism of charge-transfer process across the two-dimensional van der Waals (vdW) heterostructure remains elusive. Here, on the basis of the ab initio nonadiabatic molecular dynamics simulation, we model the photoinduced charge-transfer dynamics at the InSe/InTe vdW heterostructures. Our results show that carriers can follow either the R-scheme or Z-scheme transfer path, depending on the coupling between the interlayer states at the band-edge positions. In addition, the charge-transfer dynamics can be effectively controlled by the external parameters, such as strains and interlayer stacking configurations. The predicated electron-hole recombination lifetime in the R-scheme transfer path is up to 1.4 ns, whereas it is shortened to 1.2 ps in the Z-scheme transfer path. The proposed R-scheme and Z-scheme are further verified by the quantum transport simulations on the basis of the density functional theory (DFT) method combined with nonequilibrium Green's functions (NEGF-DFT). The analysis reveals that the system dominated by the Z-scheme shows better performance, which can be attributed to the built-in electric field that facilitates the charge transfer. Our work may pave the way for the designing of next-generation devices for light detecting and harvesting.