Combining Lindblad Master Equation and Surface Hopping to Evolve Distributions of Quantum Particles

Modeling Intermolecular Coulombic Decay with Non-Hermitian Real-Time Time-Dependent Density Functional Theory

In this work, we investigate the capability of using real-time time-dependent density functional theory (RT-TDDFT) in conjunction with a complex absorbing potential (CAP) to simulate the intermolecular Coulombic decay (ICD) processes following the ionization of an inner-valence electron. We examine the ICD dynamics in a series of noncovalent bonded dimer systems, including hydrogen-bonded and purely van der Waals (VdW)-bonded systems. In comparison to previous work, we show that RT-TDDFT simulations with a CAP correctly capture the ICD phenomenon in systems exhibiting a stronger binding energy. The calculated time scales for ICD of the studied systems are in the range of 5-50 fs, in agreement with previous studies. However, there is a breakdown in the accuracy of the methodology for the pure VdW-bonded systems. Overall, the presented RT-TDDFT/CAP methodology provides a powerful tool for differentiating between competing electronic relaxation pathways following inner-valence or core ionization without necessitating any a priori assumptions.

Stochastic Schrödinger Equation Approach to Real-Time Dynamics of Anderson-Holstein Impurities: An Open Quantum System Perspective

We develop a stochastic Schrödinger equation (SSE) framework to simulate the real-time dynamics of Anderson-Holstein (AH) impurities coupled to a continuous Fermionic bath. The bath degrees of freedom are incorporated through fluctuating terms determined by exact system-bath correlations, which is derived in an ab initio manner. We show that such an SSE treatment provides a middle ground between numerically expansive microscopic simulations and macroscopic master equations. Computationally, the SSE model enables efficient numerical methods for propagating stochastic trajectories. We demonstrate that this approach not only naturally provides microscopically detailed information unavailable from reduced models but also captures effects beyond master equations, thus serving as a promising tool to study open quantum dynamics emerging in physics and chemistry.

Electron- versus Spin-Phonon Coupling Governs the Temperature-Dependent Carrier Dynamics in the Topological Insulator Bi2Te3

Ultrafast charge and spin dynamics have immense effects on the applications of topological insulators (TIs). By performing spin-adiabatic nonadiabatic molecular dynamics simulations in the presence of electron-phonon (e-ph) and spin-phonon couplings, we investigate temperature-dependent intra- and interband charge and spin relaxation dynamics via the bulk and surface paths in the three-dimensional TI Bi2Te3. The e-ph coupling dominates charge relaxation in the bulk path, and the relaxation rate is positively correlated with temperature due to the large energy gaps and weak spin polarization. Conversely, the relaxation dynamics exhibits an opposite temperature dependence in the surface path because of electron re-excitation and spin mismatching induced by spin-phonon coupling, which arises from small energy gaps and strong spin polarization. The two mechanisms rationalize the charge carriers being long-lived in the bulk and surface phases at low and room temperature, respectively. Additionally, strong thermal fluctuations of the topological states' magnetic moments destroy the spin-momentum locking and trigger backscattering at room temperature.

Top-ranked cycle flux network analysis of molecular photocells

We introduce a top-ranked cycle flux ranking scheme of network analysis to assess the performance of molecular junction solar cells. By mapping the Lindblad master equation to the quantum-transition network, we propose a microscopic Hamiltonian description underpinning the rate equations commonly used to characterize molecular photocells. Our approach elucidates the paramount significance of edge flux and unveils two pertinent electron transfer pathways that play equally important roles in robust photocurrent generation. Furthermore, we demonstrate that nonradiative loss processes impede the maximum power efficiency of photocells, which may otherwise be above the Curzon-Ahlborn limit. These findings shed light on the intricate functionalities that govern molecular photovoltaics and offer a comprehensive approach to address them in a systematic way.

Ultrafast and Nanoscale Energy Transduction Mechanisms and Coupled Thermal Transport across Interfaces

The coupled interactions among the fundamental carriers of charge, heat, and electromagnetic fields at interfaces and boundaries give rise to energetic processes that enable a wide array of technologies. The energy transduction among these coupled carriers results in thermal dissipation at these surfaces, often quantified by the thermal boundary resistance, thus driving the functionalities of the modern nanotechnologies that are continuing to provide transformational benefits in computing, communication, health care, clean energy, power recycling, sensing, and manufacturing, to name a few. It is the purpose of this Review to summarize recent works that have been reported on ultrafast and nanoscale energy transduction and heat transfer mechanisms across interfaces when different thermal carriers couple near or across interfaces. We review coupled heat transfer mechanisms at interfaces of solids, liquids, gasses, and plasmas that drive the resulting interfacial heat transfer and temperature gradients due to energy and momentum coupling among various combinations of electrons, vibrons, photons, polaritons (plasmon polaritons and phonon polaritons), and molecules. These interfacial thermal transport processes with coupled energy carriers involve relatively recent research, and thus, several opportunities exist to further develop these nascent fields, which we comment on throughout the course of this Review.

Ab Initio Molecular Cavity Quantum Electrodynamics Simulations Using Machine Learning Models

We present a mixed quantum-classical simulation of polariton dynamics for molecule-cavity hybrid systems. In particular, we treat the coupled electronic-photonic degrees of freedom (DOFs) as the quantum subsystem and the nuclear DOFs as the classical subsystem and use the trajectory surface hopping approach to simulate non-adiabatic dynamics among the polariton states due to the coupled motion of nuclei. We use the accurate nuclear gradient expression derived from the Pauli-Fierz quantum electrodynamics Hamiltonian without making further approximations. The energies, gradients, and derivative couplings of the molecular systems are obtained from the on-the-fly simulations at the level of complete active space self-consistent field (CASSCF), which are used to compute the polariton energies and nuclear gradients. The derivatives of dipoles are also necessary ingredients in the polariton nuclear gradient expression but are often not readily available in electronic structure methods. To address this challenge, we use a machine learning model with the Kernel ridge regression method to construct the dipoles and further obtain their derivatives, at the same level as the CASSCF theory. The cavity loss process is modeled with the Lindblad jump superoperator on the reduced density of the electronic-photonic quantum subsystem. We investigate the azomethane molecule and its photoinduced isomerization dynamics inside the cavity. Our results show the accuracy of the machine-learned dipoles and their usage in simulating polariton dynamics. Our polariton dynamics results also demonstrate the isomerization reaction of azomethane can be effectively tuned by coupling to an optical cavity and by changing the light-matter coupling strength and the cavity loss rate.

Incorporating Lindblad Decay Dynamics into Mixed Quantum-Classical Simulations

We derive the $\mathcal{L}$-MFE method to incorporate Lindblad jump operator dynamics into the mean-field Ehrenfest (MFE) approach. We map the density matrix evolution of Lindblad dynamics onto pure state coefficients using trajectory averages. We use simple assumptions to construct the $\mathcal{L}$-MFE method that satisfies this exact mapping. This establishes a method that uses independent trajectories which exactly reproduces Lindblad decay dynamics using a wavefunction description, with deterministic changes of the magnitudes of the quantum expansion coefficients, while only adding on a stochastic phase. We further demonstrate that when including nuclei in the Ehrenfest dynamics, the $\mathcal{L}$-MFE method gives semi-quantitatively accurate results, with the accuracy limited by the accuracy of the approximations present in the semiclassical MFE approach. This work provides a general framework to incorporate Lindblad dynamics into semiclassical or mixed quantum-classical simulations.

Modeling Non-adiabatic Dynamics in Nanoscale and Condensed Matter Systems

Rapid, far-from-equilibrium processes involving excitation of electronic, vibrational, spin, photon, topological, and other degrees of freedom form the basis of modern technologies, including electronics and optoelectronics, solar energy harvesting and conversion to electrical and chemical energy, quantum information processing, spin- and valleytronics, chemical detection, and medical therapies. Such processes are studied experimentally with various time-resolved spectroscopies that allow scientists to track system's evolution on ultrafast time scales and at close to atomistic level of detail. The availability of various forms of lasing has made such measurements easily accessible to many experimental groups worldwide, to study atoms and small molecules, nanoscale and condensed matter systems, proteins, cells, and mesoscopic materials. The experimental work necessitates parallel theoretical efforts needed to interpret the experiments and to provide insights that cannot be gained through measurements due to experimental limitations.Non-adiabatic (NA) molecular dynamics (MD) allows one to study processes at the atomistic level and in the time domain most directly mimicking the time-resolved experiments. Atomistic modeling takes full advantage of chemical intuition and principles that guide design and fabrication of molecules and materials. It provides atomistic origins of quasi-particles, such as holes, excitons, trions, plasmons, phonons, polarons, polaritons, spin-waves, momentum-resolved and topological states, electrically and magnetically polarized structures, and other abstract concepts. An atomistic description enables one to study realistic aspects of materials, which necessarily contain defects, dopants, surfaces, interfaces, passivating ligands, and solvent layers. Often, such realistic features govern material properties and are hard to account for phenomenologically. NA-MD requires few approximations and assumptions. It does not need to assume that atomic motions are harmonic, that electrons are Drude oscillators, that coupling between different degrees of freedom is weak, that dynamics is Markovian or has short memory, or that evolution occurs by exponential kinetics of transitions between few states. The classical or semiclassical treatment of atomic motions constitutes the main approximation of NA-MD and is used because atoms are 3-5 orders of magnitude heavier than electrons. NA-MD is limited by system size, typically hundreds or thousands of atoms, and time scale, picoseconds to nanoseconds. The quality of NA-MD simulations depends on the electronic structure method used to obtain excited state energies and NA couplings.NA-MD has been largely popularized and advanced in the chemistry community that focuses on molecules. Modeling far-from-equilibrium dynamics in nanoscale and condensed matter systems often has to account for other types of physics. At the same time, condensed phase NA-MD allows for approximations that may not work in molecules. Focusing on the recent NA-MD developments aimed at studying excited state processes in nanoscale and condensed phases, this Account considers how the phenomena important on the nanoscale can be incorporated into NA-MD and what approximations can be made to increase its efficiency with complex systems and processes.

Approaching Disordered Quantum Dot Systems by Complex Networks with Spatial and Physical-Based Constraints

This paper focuses on modeling a disordered system of quantum dots (QDs) by using complex networks with spatial and physical-based constraints. The first constraint is that, although QDs (=nodes) are randomly distributed in a metric space, they have to fulfill the condition that there is a minimum inter-dot distance that cannot be violated (to minimize electron localization). The second constraint arises from our process of weighted link formation, which is consistent with the laws of quantum physics and statistics: it not only takes into account the overlap integrals but also Boltzmann factors to include the fact that an electron can hop from one QD to another with a different energy level. Boltzmann factors and coherence naturally arise from the Lindblad master equation. The weighted adjacency matrix leads to a Laplacian matrix and a time evolution operator that allows the computation of the electron probability distribution and quantum transport efficiency. The results suggest that there is an optimal inter-dot distance that helps reduce electron localization in QD clusters and make the wave function better extended. As a potential application, we provide recommendations for improving QD intermediate-band solar cells.

The 2021 ultrafast spectroscopic probes of condensed matter roadmap

In the 60 years since the invention of the laser, the scientific community has developed numerous fields of research based on these bright, coherent light sources, including the areas of imaging, spectroscopy, materials processing and communications. Ultrafast spectroscopy and imaging techniques are at the forefront of research into the light-matter interaction at the shortest times accessible to experiments, ranging from a few attoseconds to nanoseconds. Light pulses provide a crucial probe of the dynamical motion of charges, spins, and atoms on picosecond, femtosecond, and down to attosecond timescales, none of which are accessible even with the fastest electronic devices. Furthermore, strong light pulses can drive materials into unusual phases, with exotic properties. In this roadmap we describe the current state-of-the-art in experimental and theoretical studies of condensed matter using ultrafast probes. In each contribution, the authors also use their extensive knowledge to highlight challenges and predict future trends.

Ab initio nonadiabatic molecular dynamics of charge carriers in metal halide perovskites

Photoinduced nonequilibrium processes in nanoscale materials play key roles in photovoltaic and photocatalytic applications. This review summarizes recent theoretical investigations of excited state dynamics in metal halide perovskites (MHPs), carried out using a state-of-the-art methodology combining nonadiabatic molecular dynamics with real-time time-dependent density functional theory. The simulations allow one to study evolution of charge carriers at the ab initio level and in the time-domain, in direct connection with time-resolved spectroscopy experiments. Eliminating the need for the common approximations, such as harmonic phonons, a choice of the reaction coordinate, weak electron-phonon coupling, a particular kinetic mechanism, and perturbative calculation of rate constants, we model full-dimensional quantum dynamics of electrons coupled to semiclassical vibrations. We study realistic aspects of material composition and structure and their influence on various nonequilibrium processes, including nonradiative trapping and relaxation of charge carriers, hot carrier cooling and luminescence, Auger-type charge-charge scattering, multiple excitons generation and recombination, charge and energy transfer between donor and acceptor materials, and charge recombination inside individual materials and across donor/acceptor interfaces. These phenomena are illustrated with representative materials and interfaces. Focus is placed on response to external perturbations, formation of point defects and their passivation, mixed stoichiometries, dopants, grain boundaries, and interfaces of MHPs with charge transport layers, and quantum confinement. In addition to bulk materials, perovskite quantum dots and 2D perovskites with different layer and spacer cation structures, edge passivation, and dielectric screening are discussed. The atomistic insights into excited state dynamics under realistic conditions provide the fundamental understanding needed for design of advanced solar energy and optoelectronic devices.