Perturbation expansions of stochastic wavefunctions for open quantum systems

Phonon-mediated ultrafast energy- and momentum-resolved hole dynamics in monolayer black phosphorus

The electron-phonon scattering plays a crucial role in determining the electronic, transport, optical, and thermal properties of materials. Here, we employ a non-Markovian stochastic Schrödinger equation (NMSSE) in momentum space, together with ab initio calculations for energy bands and electron-phonon interactions, to reveal the phonon-mediated ultrafast hole relaxation dynamics in the valence bands of monolayer black phosphorus. Our numerical simulations show that the hole can initially remain in the high-energy valence bands for more than 100 fs due to the weak interband scatterings, and its energy relaxation follows single-exponential decay toward the valence band maximum after scattering into low-energy valence bands. The total relaxation time of holes is much longer than that of electrons in the conduction band. This suggests that harnessing the excess energy of holes may be more effective than that of electrons. Compared to the semiclassical Boltzmann equation based on a hopping model, the NMSSE highlights the persistence of quantum coherence for a long time, which significantly impacts the relaxation dynamics. These findings complement the understanding of hot carrier relaxation dynamics in two-dimensional materials and may offer novel insights into harnessing hole energy in photocatalysis.

Phonon-Mediated Ultrafast Electron Relaxation Dynamics in Monolayer Black Phosphorus: Instantaneous Coherent Delocalization

Understanding carrier relaxation processes in semiconductors is crucial for designing high-performance optoelectronic and photocatalytic devices. Recent transient spectroscopic experiments on two-dimensional materials have revealed ultrafast optical responses within several tens of femtoseconds, which are usually ascribed to electron-electron scattering. Here, by conducting quantum dynamics simulations for monolayer black phosphorus, we show that electron-phonon scattering also profoundly influences the early stage of carrier dynamics. The photogenerated electron generally undergoes phonon-mediated instantaneous coherent delocalization in reciprocal space, accompanied by an entropy-driven sharp change in electronic energy. The distribution of the density of states controls the energy exchange between the electron and lattice vibrations. The phonon-induced quantum coherence significantly suspends the energy relaxation time, which is very beneficial for harvesting electron excess energy. These findings offer novel insights into the ultrafast carrier dynamics and energy flow in two-dimensional materials and may prompt new opportunities for regulation of carrier dynamic behaviors.

Substituent Effect on Vibrationally Resolved Absorption Spectra and Exciton Dynamics of Dipyrrolonaphthyridinedione Aggregates

Dipyrrolonaphthyridinedione (DPND) thin films exhibit interesting photophysical properties and singlet fission (SF) processes. A recent experimental work found that the alkyl substitution in the DPND skeleton has the remarkable influence on the characteristics of electronic absorption spectra and SF rates. Here, we theoretically elucidate the microscopic mechanism of the substituent effect on the optical properties and exciton dynamics of materials by combining the electronic structure calculations and the quantum dynamics simulations. The results show that the alkyl substituent has a minor effect on the single molecular properties but dramatically changes those of DPND aggregates via varying the intermolecular interactions. The aggregates of DPND with and without alkyl side chains exhibit the more likely characters of H-type aggregations. In the former (DPND6), the weak degree of mixing of intramolecular localized excited (LE) states and intermolecular charge transfer (CT) states makes the low-energy absorption band possess the predominant optical absorption, while in the latter (DPND), the CT and LE states are close in energy, together with their strong interaction, resulting in the substantial state-mixing, so that its two low-energy absorption bands have nearly equal oscillator strengths and a wide energy spacing of more than 0.5 eV. The simulation of exciton dynamics elucidates that the photoinitiated states in both aggregates cannot generate the free charge carrier because of the lack of enough driving forces. However, the population exchanges between LE and CT states in DPND aggregates are much faster than in DPND6 aggregates, indicating the different SF behaviors, consistent with the experimental observation.

Vibrationally resolved absorption spectra and ultrafast exciton dynamics in α-phase and β-phase zinc phthalocyanine aggregates

The vibrationally resolved absorption spectra and ultrafast exciton dynamics in α-phase and β-phase zinc phthalocyanine (ZnPc) aggregates are theoretically investigated using a non-Markovian stochastic Schrödinger equation combined with first-principles calculations. It is found that although similar double-peak structures arise in the Q-band region of the absorption spectra in both phases, these peaks are different in nature and exhibit distinct types of behavior with respect to the aggregation length. The analysis on the basis of an effective two-state model indicates that the two absorption peaks in the α phase are from mixing between the charge-transfer (CT) state and the bright Frenkel exciton (FE) state. By contrast, in the β-phase, the low-energy peak is solely contributed by a low-lying bright FE state, whereas the high-energy peak originates from the interplay between the CT state and another high-lying bright FE state. For the relaxation processes right after photoexcitation from the Q-band region, it is found that within the first dozens of femtoseconds the ZnPc aggregates of both phases tend to temporarily fall into some intermediate states where the population distribution and average electronic energy do not obviously evolve. In addition, it is found that the optical transition of the low-lying bright FE state in the β phase is not favorable for the formation of bound CT states due to the absence of enough driving forces.

Vibrationally Resolved Absorption Spectra and Exciton Dynamics in Zinc Phthalocyanine Aggregates: Effects of Aggregation Lengths and Remote Exciton Transfer

The vibrationally resolved absorption spectra and exciton dynamics in the α-zinc phthalocyanine aggregates are theoretically investigated by using a non-Markovian stochastic Schrödinger equation. The model Hamiltonian adopted for spectral and dynamic simulations explicitly includes the couplings for both nearest-neighbor and remote exciton transfer, and it is parametrized from first-principles calculations. The results indicate that aggregation lengths and remote exciton transfer significantly influence the relative energy alignment between delocalized Frenkel exciton (FE) and charge transfer (CT) states, which in turn strongly affects the relative intensities of the two absorption peaks in the Q-band region. Analytical formulas are derived to establish quantitative structure-spectra relationships in aggregates, and they offer simple patterns to extract electronic-state properties directly from absorption spectra. The dynamics simulations reveal that the light absorption can directly generate mixed states with both FE and CT features, but it is hard for the photoexcitation from the Q-band region to generate free carriers due to the high energies of charge-separated states.

A unified view of hierarchy approach and formula of differentiation

The stochastic differential equation is a powerful tool for describing the dynamics of a dissipative system in which noise characterizes the influence of the environment. For the Ornstein-Uhlenbeck noise, both the formula of differentiation and the hierarchy approach provide efficient numerical simulations, with the stochastic differential equation transformed into a set of coupled, linear ordinary differential equations. We show that while these two deterministic schemes result in different sets of equations, they can be regarded as two representations of an underlying linear-dynamics. Moreover, by manipulating the involved Ornstein-Uhlenbeck noise, we propose a unified algorithm that may reduce to the hierarchy approach or the formula of differentiation in different limits. We further analyze the numerical performance of this algorithm and find that the hierarchy approach appears to be more efficient for our numerical model studies.