Intermolecular conical intersections in molecular aggregates

Machine learning photodynamics decode multiple singlet fission channels in pentacene crystal

Crystalline pentacene is a model solid-state light-harvesting material because its quantum efficiencies exceed 100% via ultrafast singlet fission. The singlet fission mechanism in pentacene crystals is disputed due to insufficient electronic information in time-resolved experiments and intractable quantum mechanical calculations for simulating realistic crystal dynamics. Here we combine a multiscale multiconfigurational approach and machine learning photodynamics to understand competing singlet fission mechanisms in crystalline pentacene. Our simulations reveal coexisting charge-transfer-mediated and coherent mechanisms via the competing channels in the herringbone and parallel dimers. The predicted singlet fission time constants (61 and 33 fs) are in excellent agreement with experiments (78 and 35 fs). The trajectories highlight the essential role of intermolecular stretching between monomers in generating the multi-exciton state and explain the anisotropic phenomenon. The machine-learning-photodynamics resolved the elusive interplay between electronic structure and vibrational relations, enabling fully atomistic excited-state dynamics with multiconfigurational quantum mechanical quality for crystalline pentacene.

Ultrafast Time-Domain Spectroscopy Reveals Coherent Vibronic Couplings upon Electronic Excitation in Crystalline Organic Thin Films

The coherent coupling between electronic excitations and vibrational modes of molecules largely affects the optical and charge transport properties of organic semiconductors and molecular solids. To analyze these couplings by means of ultrafast spectroscopy, highly ordered crystalline films with large domains are particularly suitable because the domains can be addressed individually, hence allowing azimuthal polarization-resolved measurements. Impressive examples of this are highly ordered crystalline thin films of perfluoropentacene (PFP) molecules, which adopt different molecular orientations on different alkali halide substrates. Here, we report polarization-resolved time-domain vibrational spectroscopy with 10 fs time resolution and Raman spectroscopy of crystalline PFP thin films grown on NaF(100) and KCl(100) substrates. Coherent oscillations in the time-resolved spectra reveal vibronic coupling to a high-frequency, 25 fs, in-plane deformation mode that is insensitive to the optical polarization, while the coupling to a lower-frequency, 85 fs, out-of-plane ring bending mode depends significantly on the crystalline and molecular orientation. Comparison with calculated Raman spectra of isolated PFP molecules in vacuo supports this interpretation and indicates a dominant role of solid-state effects in the vibronic properties of these materials. Our results represent a first step toward uncovering the role of anisotropic vibronic couplings for singlet fission processes in crystalline molecular thin films.

Quantum Computation of Conical Intersections on a Programmable Superconducting Quantum Processor

Conical intersections (CIs) are pivotal in many photochemical processes. Traditional quantum chemistry methods, such as the state-average multiconfigurational methods, face computational hurdles in solving the electronic Schrödinger equation within the active space on classical computers. While quantum computing offers a potential solution, its feasibility in studying CIs, particularly on real quantum hardware, remains largely unexplored. Here, we present the first successful realization of a hybrid quantum-classical state-average complete active space self-consistent field method based on the variational quantum eigensolver (VQE-SA-CASSCF) on a superconducting quantum processor. This approach is applied to investigate CIs in two prototypical systems─ethylene (C2H4) and triatomic hydrogen (H3). We illustrate that VQE-SA-CASSCF, coupled with ongoing hardware and algorithmic enhancements, can lead to a correct description of CIs on existing quantum devices. These results lay the groundwork for exploring the potential of quantum computing to study CIs in more complex systems in the future.

Unveiling Multiquantum Excitonic Correlations in Push-Pull Polymer Semiconductors

Bound and unbound Frenkel-exciton pairs are essential transient precursors for a variety of photophysical and biochemical processes. In this work, we identify bound and unbound Frenkel-exciton complexes in an electron push-pull polymer semiconductor using coherent two-dimensional spectroscopy. We find that the dominant A0-1 peak of the absorption vibronic progression is accompanied by a subpeak, each dressed by distinct vibrational modes. By considering the Liouville pathways within a two-exciton model, the imbalanced cross-peaks in one-quantum rephasing and nonrephasing spectra can be accounted for by the presence of pure biexcitons. The two-quantum nonrephasing spectra provide direct evidence for unbound exciton pairs and biexcitons with dominantly attractive force. In addition, the spectral features of unbound exciton pairs show mixed absorptive and dispersive character, implying many-body interactions within the correlated Frenkel-exciton pairs. Our work offers novel perspectives on the Frenkel-exciton complexes in semiconductor polymers.

Probing Conical Intersection in the Multipathway Isomerization of CH3Cl Using Coulomb Explosion

Ubiquitous ultrafast isomerization is paramount in photoexcited molecules, in which non-adiabatic coupling among multiple electronic states can occur. We use the pump-probe Coulomb explosion imaging method to study the isomerization of CH3Cl molecules. We find that the isomerization under our strong field pump-probe scheme proceeds along multiple pathways, which are encoded in several distinct branches of the time-resolved kinetic energy release spectra for the CH2++HCl+ Coulomb explosion channel. Apart from the isomerized dissociative pathway in neutral and cationic excited states, the pump laser can also induce coherent vibrational dynamics in two coupled intermediate states and set up the initial conditions for the two concurrently proceeding isomerization pathways. The isomerization of CH3Cl provides an intriguing example of a chemical reaction consisting of multiple pathways and non-adiabatic dynamics.

Full visible range two-dimensional electronic spectroscopy with high time resolution

Two-dimensional electronic spectroscopy (2DES) is a powerful method to study coherent and incoherent interactions and dynamics in complex quantum systems by correlating excitation and detection energies in a nonlinear spectroscopy experiment. Such dynamics can be probed with a time resolution limited only by the duration of the employed laser pulses and in a spectral range defined by the pulse spectrum. In the blue spectral range (<500 nm), the generation of sufficiently broadband ultrashort pulses with pulse durations of 10 fs or less has been challenging so far. Here, we present a 2DES setup based on a hollow-core fiber supercontinuum covering the full visible range (400-700 nm). Pulse compression via custom-made chirped mirrors yields a time resolution of <10 fs. The broad spectral coverage, in particular the extension of the pulse spectra into the blue spectral range, unlocks new possibilities for coherent investigations of blue-light absorbing and multichromophoric compounds, as demonstrated by a 2DES measurement of chlorophyll a.

Plasmon mediated coherent population oscillations in molecular aggregates

The strong coherent coupling of quantum emitters to vacuum fluctuations of the light field offers opportunities for manipulating the optical and transport properties of nanomaterials, with potential applications ranging from ultrasensitive all-optical switching to creating polariton condensates. Often, ubiquitous decoherence processes at ambient conditions limit these couplings to such short time scales that the quantum dynamics of the interacting system remains elusive. Prominent examples are strongly coupled exciton-plasmon systems, which, so far, have mostly been investigated by linear optical spectroscopy. Here, we use ultrafast two-dimensional electronic spectroscopy to probe the quantum dynamics of J-aggregate excitons collectively coupled to the spatially structured plasmonic fields of a gold nanoslit array. We observe rich coherent Rabi oscillation dynamics reflecting a plasmon-driven coherent exciton population transfer over mesoscopic distances at room temperature. This opens up new opportunities to manipulate the coherent transport of matter excitations by coupling to vacuum fields.

Time-Linear Quantum Transport Simulations with Correlated Nonequilibrium Green's Functions

We present a time-linear scaling method to simulate open and correlated quantum systems out of equilibrium. The method inherits from many-body perturbation theory the possibility to choose selectively the most relevant scattering processes in the dynamics, thereby paving the way to the real-time characterization of correlated ultrafast phenomena in quantum transport. The open system dynamics is described in terms of an "embedding correlator" from which the time-dependent current can be calculated using the Meir-Wingreen formula. We show how to efficiently implement our approach through a simple grafting into recently proposed time-linear Green's function methods for closed systems. Electron-electron and electron-phonon interactions can be treated on equal footing while preserving all fundamental conservation laws.

Coherent Vibronic Wavepackets Show Structure-Directed Charge Flow in Host-Guest Donor-Acceptor Complexes

Designing and controlling charge transfer (CT) pathways in organic semiconductors are important for solar energy applications. To be useful, a photogenerated, Coulombically bound CT exciton must further separate into free charge carriers; direct observations of the detailed CT relaxation pathways, however, are lacking. Here, photoinduced CT and relaxation dynamics in three host-guest complexes, where a perylene (Per) electron donor guest is incorporated into two symmetric and one asymmetric extended viologen cyclophane acceptor hosts, are presented. The central ring in the extended viologen is either p-phenylene (ExV²⁺) or electron-rich 2,5-dimethoxy-p-phenylene (ExMeOV²⁺), resulting in two symmetric cyclophanes with unsubstituted or methoxy-substituted central rings, ExBox⁴⁺ and ExMeOBox⁴⁺, respectively, and an asymmetric cyclophane with one of the central viologen rings being methoxylated ExMeOVBox⁴⁺. Upon photoexcitation, the asymmetric host-guest ExMeOVBox⁴⁺ ⊃ Per complex exhibits directional CT toward the energetically unfavorable methoxylated side due to structural restrictions that facilitate strong interactions between the Per donor and the ExMeOV²⁺ side. The CT state relaxation pathways are probed using ultrafast optical spectroscopy by focusing on coherent vibronic wavepackets, which are used to identify CT relaxations along charge localization and vibronic decoherence coordinates. Specific low- and high-frequency nuclear motions are direct indicators of a delocalized CT state and the degree of CT character. Our results show that the CT pathway can be controlled by subtle chemical modifications of the acceptor host in addition to illustrating how coherent vibronic wavepackets can be used to probe the nature and time evolution of the CT states.

Intermolecular-Type Conical Intersections in Benzene Dimer

The equilibrium and conical intersection geometries of the benzene dimer were computed in the framework of the conventional, linear-response time-dependent and spin-flipped time-dependent density functional theories (known as DFT, TDDFT and SF-TDDFT) as well as using the multiconfigurational complete active space self-consistent field (CASSCF) method considering the minimally augmented def2-TZVPP and the 6-31G(d,p) basis sets. It was found that the stacking distance between the benzene monomers decreases by about 0.5 Å in the first electronic excited state, due to the stronger intermolecular interaction energy, bringing the two monomers closer together. Intermolecular-type conical intersection (CI) geometries can be formed between the two benzene molecules, when (i) both monomer rings show planar deformation and (ii) weaker (approximately 1.6-1.8 Å long) C-C bonds are formed between the two monomers, with parallel and antiparallel orientation with respect to the monomer. These intermolecular-type CIs look energetically more favorable than dimeric CIs containing only one deformed monomer. The validity of the dimer-type CI geometries obtained by SF-TDDFT was confirmed by the CASSCF method. The nudged elastic band method used for finding the optimal relaxation path has confirmed both the accessibility of these intermolecular-type CIs and the possibility of the radiationless deactivation of the electronic excited states through these CI geometries. Although not as energetically favorable as the previous two CI geometries, there are other CI geometries characterized by the relative rotation of monomers at different angles around a vertical C-C axis.

Mixing of Excitons in Nanostructures Based on a Perylene Dye with CdTe Quantum Dots

Semiconductor quantum dots of the A₂B₆ group and organic semiconductors have been widely studied and applied in optoelectronics. This study aims to combine CdTe quantum dots and perylene-based dye molecules into advanced nanostructure system targeting to improve their functional properties. In such systems, new electronic states, a mixture of Wannier-Mott excitons with charge-transfer excitons, have appeared at the interface of CdTe quantum dots and the perylene dye. The nature of such new states has been analyzed by absorption and photoluminescence spectroscopy with picosecond time resolution. Furthermore, aggregation of perylene dye on the CdTe has been elucidated, and contribution of Förster resonant energy transfer has been observed between aggregated forms of the dye and CdTe quantum dots in the hybrid CdTe-perylene nanostructures. The studied nanostructures have strongly quenched emission of quantum dots enabling potential application of such systems in dissociative sensing.

From a one-mode to a multi-mode understanding of conical intersection mediated ultrafast organic photochemical reactions

This review discusses how ultrafast organic photochemical reactions are controlled by conical intersections, highlighting that decay to the ground-state at multiple points of the intersection space results in their multi-mode character.

Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals

Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.

Charge Delocalization and Vibronic Couplings in Quadrupolar Squaraine Dyes

Squaraines are prototypical quadrupolar charge-transfer chromophores that have recently attracted much attention as building blocks for solution-processed photovoltaics, fluorescent probes with large two-photon absorption cross sections, and aggregates with large circular dichroism. Their optical properties are often rationalized in terms of phenomenological essential state models, considering the coupling of two zwitterionic excited states to a neutral ground state. As a result, optical transitions to the lowest S1 excited state are one-photon allowed, whereas the next higher S2 state can only be accessed by two-photon transitions. A further implication of these models is a substantial reduction of vibronic coupling to the ubiquitous high-frequency vinyl-stretching modes of organic materials. Here, we combine time-resolved vibrational spectroscopy, two-dimensional electronic spectroscopy, and quantum-chemical simulations to test and rationalize these predictions for nonaggregated molecules. We find small Huang-Rhys factors below 0.01 for the high-frequency, 1500 cm^-1 modes in particular, as well as a noticeable reduction for those of lower frequency modes in general for the electronic S0 → S1 transition. The two-photon allowed state S2 is well separated energetically from S1 and has weak vibronic signatures as well. Thus, the resulting pronounced concentration of the oscillator strength in a narrow region relevant to the lowest electronic transition makes squaraines and their aggregates exceptionally interesting for strong and ultrastrong coupling of excitons to localized light modes in external resonators with chiral properties that can largely be controlled by the molecular architecture.

Plasmon-Enhanced Exciton Delocalization in Squaraine-Type Molecular Aggregates

Enlarging exciton coherence lengths in molecular aggregates is critical for enhancing the collective optical and transport properties of molecular thin film nanostructures or devices. We demonstrate that the exciton coherence length of squaraine aggregates can be increased from 10 to 24 molecular units at room temperature when preparing the aggregated thin film on a metallic rather than a dielectric substrate. Two-dimensional electronic spectroscopy measurements reveal a much lower degree of inhomogeneous line broadening for aggregates on a gold film, pointing to a reduced disorder. The result is corroborated by simulations based on a Frenkel exciton model including exciton-plasmon coupling effects. The simulation shows that localized, energetically nearly resonant excitons on spatially well separated segments can be radiatively coupled via delocalized surface plasmon polariton modes at a planar molecule-gold interface. Such plasmon-enhanced delocalization of the exciton wave function is of high importance for improving the coherent transport properties of molecular aggregates on the nanoscale. Additionally, it may help tailor the collective optical response of organic materials for quantum optical applications.

Nonadiabatic Dynamics Algorithms with Only Potential Energies and Gradients: Curvature-Driven Coherent Switching with Decay of Mixing and Curvature-Driven Trajectory Surface Hopping

Direct dynamics by mixed quantum-classical nonadiabatic methods is an important tool for understanding processes involving multiple electronic states. Very often, the computational bottleneck of such direct simulation comes from electronic structure theory. For example, at every time step of a trajectory, nonadiabatic dynamics requires potential energy surfaces, their gradients, and the matrix elements coupling the surfaces. The need for the couplings can be alleviated by employing the time derivatives of the wave functions, which can be evaluated from overlaps of electronic wave functions at successive time steps. However, evaluation of overlap integrals is still expensive for large systems. In addition, for electronic structure methods for which the wave functions or the coupling matrix elements are not available, nonadiabatic dynamics algorithms become inapplicable. In this work, building on recent work by Baeck and An, we propose new nonadiabatic dynamics algorithms that only require adiabatic potential energies and their gradients. The new methods are named curvature-driven coherent switching with decay of mixing (κCSDM) and curvature-driven trajectory surface hopping (κTSH). We show how powerful these new methods are in terms of computation time and accuracy as compared to previous mixed quantum-classical nonadiabatic dynamics algorithms. The lowering of the computational cost will allow longer nonadiabatic trajectories and greater ensemble averaging to be affordable, and the ability to calculate the dynamics without electronic structure coupling matrix elements extends the dynamics capability to new classes of electronic structure methods.

Distinguishing between coherent and incoherent signals in excitation-emission spectroscopy

The separation of incoherent emission signals from coherent light scattering often poses a challenge in (time-resolved) microscopy or excitation-emission spectroscopy. While in spectro-microscopy with narrowband excitation this is commonly overcome using spectral filtering, it is less straightforward when using broadband Fourier-transform techniques that are now becoming commonplace in, e.g., single molecule or ultrafast nonlinear spectroscopy. Here we show that such a separation is readily achieved using highly stable common-path interferometers for both excitation and detection. The approach is demonstrated for suppression of scattering from flavin adenine dinucleotide (FAD) and weakly emissive cryptochrome 4 (Cry4) protein samples. We expect that the approach will be beneficial, e.g., for fluorescence lifetime or Raman-based imaging and spectroscopy of various samples, including single quantum emitters.

Multimode two-dimensional vibronic spectroscopy. II. Simulating and extracting vibronic coupling parameters from polarization-selective spectra

Experimental demonstrations of polarization-selection two-dimensional Vibrational-Electronic (2D VE) and 2D Electronic-Vibrational (2D EV) spectroscopies aim to map the magnitudes and spatial orientations of coupled electronic and vibrational coordinates in complex systems. The realization of that goal depends on our ability to connect spectroscopic observables with molecular structural parameters. In this paper, we use a model Hamiltonian consisting of two anharmonically coupled vibrational modes in electronic ground and excited states with linear and bilinear vibronic coupling terms to simulate polarization-selective 2D EV and 2D VE spectra. We discuss the relationships between the linear vibronic coupling and two-dimensional Huang-Rhys parameters and between the bilinear vibronic coupling term and Duschinsky mixing. We develop a description of the vibronic transition dipoles and explore how the Hamiltonian parameters and non-Condon effects impact their amplitudes and orientations. Using simulated polarization-selective 2D EV and 2D VE spectra, we show how 2D peak positions, amplitudes, and anisotropy can be used to measure parameters of the vibronic Hamiltonian and non-Condon effects. This paper, along with the first in the series, provides the reader with a detailed description of reading, simulating, and analyzing multimode, polarization-selective 2D EV and 2D VE spectra with an emphasis on extracting vibronic coupling parameters from complex spectra.

Effect of varying the TD-lc-DFTB range-separation parameter on charge and energy transfer in a model pentacene/buckminsterfullerene heterojunction

Atomistic modeling of energy and charge transfer at the heterojunction of organic solar cells is an active field with many remaining outstanding questions owing, in part, to the difficulties in performing reliable photodynamics calculations on very large systems. One approach to being able to overcome these difficulties is to design and apply an appropriate simplified method. Density-functional tight binding (DFTB) has become a popular form of approximate density-functional theory based on a minimal valence basis set and neglect of all but two center integrals. We report the results of our tests of a recent long-range correction (lc) [A. Humeniuk and R. Mitrić, J. Chem. Phys. 143, 134120 (2015)] for time-dependent (TD) lc-DFTB by carrying out TD-lc-DFTB fewest switches surface hopping calculations of energy and charge transfer times using the relatively new DFTBABY [A. Humeniuk and R. Mitrić, Comput. Phys. Commun. 221, 174 (2017)] program. An advantage of this method is the ability to run enough trajectories to get meaningful ensemble averages. Our interest in the present work is less in determining exact energy and charge transfer rates than in understanding how the results of these calculations vary with the value of the range-separation parameter (R_lc = 1/μ) for a model organic solar cell heterojunction consisting of a gas-phase van der Waals complex P/F made up of a single pentacene (P) molecule together with a single buckminsterfullerene (F) molecule. The default value of R_lc = 3.03 a₀ is found to be much too small as neither energy nor charge transfer is observed until R_lc ≈ 10 a₀. Tests at a single geometry show that the best agreement with high-quality ab initio spectra is obtained in the limit of no lc (i.e., very large R_lc). A plot of energy and charge transfer rates as a function of R_lc is provided, which suggests that a value of R_lc ≈ 15 a₀ yields the typical literature (condensed-phase) charge transfer time of about 100 fs. However, energy and charge transfer times become as high as ∼300 fs for R_lc ≈ 25 a₀. A closer examination of the charge transfer process P^*/F → P⁺/F^- shows that the initial electron transfer is accompanied by a partial delocalization of the P hole onto F, which then relocalizes back onto P, consistent with a polaron-like picture in which the nuclei relax to stabilize the resultant redistribution of charges.