Coherence, signifying concurrent electron-vibrational dynamics in complex natural and man-made systems, is currently a subject of intense study. Understanding this phenomenon is important when designing carrier transport in optoelectronic materials. Here, excited state dynamics simulations reveal a ubiquitous pattern in the evolution of photoexcitations for a broad range of molecular systems. Symmetries of the wavefunctions define a specific form of the non-adiabatic coupling that drives quantum transitions between excited states, leading to a collective asymmetric vibrational excitation coupled to the electronic system. This promotes periodic oscillatory evolution of the wavefunctions, preserving specific phase and amplitude relations across the ensemble of trajectories. The simple model proposed here explains the appearance of coherent exciton-vibrational dynamics due to non-adiabatic transitions, which is universal across multiple molecular systems. The observed relationships between electronic wavefunctions and the resulting functionalities allows us to understand, and potentially manipulate, excited state dynamics and energy transfer in molecular materials.

Dendrimers are arrays of coupled chromophores, where the energy of each unit depends on its structure and conformation. The light harvesting and energy funneling properties are strongly dependent on their highly branched conjugated architecture. Herein, the photoexcitation and subsequent ultrafast electronic energy relaxation and redistribution of a first generation dendrimer (1) are analyzed combining theoretical and experimental studies. Dendrimer 1 consists of three linear phenylene-ethynylene (PE) units, or branches, attached in the meta position to a central group opening up the possibility of inter-branch energy transfer. Excited state dynamics are explored using both time-resolved spectroscopy and non-adiabatic excited state molecular dynamics simulations. Our results indicate a subpicosecond loss of anisotropy due to an initial excitation into several states with different spatial localizations, followed by exciton self-trapping on different units. This exciton hops between branches. The absence of an energy gradient leads to an ultrafast energy redistribution among isoenergetic chromophore units. At long times we observe similar probabilities for each branch to retain significant contributions of the transition density of the lowest electronic excited-state. The observed unpolarized emission is attributed to the contraction of the electronic wavefunction onto a single branch with frequent interbranch hops, and not to its delocalization over the whole dendrimer.

The non-adiabatic excited state molecular dynamics (NA-ESMD) approach is applied to investigate photoexcited dynamics and relaxation pathways in a spiro-linked conjugated polyfluorene at room (T = 300 K) and low (T = 10 K) temperatures. This dimeric aggregate consists of two perpendicularly oriented weakly interacting α-polyfluorene oligomers. The negligible coupling between the monomer chains results in an initial absorption band composed of equal contributions of the two lowest excited electronic states, each localized on one of the two chains. After photoexcitation, an efficient ultrafast localization of the entire electronic population to the lowest excited state is observed on the time scale of about 100 fs. Both internal conversion between excited electronic states and vibronic energy relaxation on a single electronic state contribute to this process. Thus, photoexcited dynamics of the polyfluorene dimer follows two distinct pathways with substantial temperature dependence on their efficiency. One relaxation channel involves resonance electronic energy transfer between the monomer chains, whereas the second pathway concerns the relaxation of the electronic energy on the same chain that has been initially excited due to electron-phonon coupling. Despite the slower vibrational relaxation, a more efficient ultrafast electronic relaxation is observed at low temperature. Our numerical simulations analyze the effects of molecular geometry distortion during the electronic energy redistribution and suggest spectroscopic signatures reflecting complex electron-vibrational dynamics.

High-level theoretical benchmark investigations of the UV-vis absorption spectra of paradigmatic polycyclic aromatic hydrocarbons as models for graphene quantum dots

We explore the diverse origins of unpolarized absorption and emission of molecular polygons, consisting of π-conjugated oligomer chains held in a bent geometry by strain controlled at the vertex units. For this purpose, we make use of atomistic nonadiabatic excited-state molecular dynamics simulations of a bichromophore molecular polygon (digon) with bent chromophore chains. Both structural and photoexcited dynamics were found to affect polarization features. Bending strain induces exciton localization on individual chromophore units of the conjugated chains. The latter display different transition dipole moment orientations, a feature not present in the linear oligomer counterparts. In addition, bending makes exciton localization very sensitive to molecular distortions induced by thermal fluctuations. The excited-state dynamics reveals an ultrafast intramolecular energy redistribution that spreads the exciton equally among spatially separated chromophore fragments within the molecular system. As a result, digons become virtually unpolarized absorbers and emitters, in agreement with recent experimental studies on the single-molecule level.

Photoexcitation of multichromophoric light harvesting molecules induces a number of intramolecular electronic energy relaxation and redistribution pathways that can ultimately lead to ultrafast exciton self-trapping on a single chromophore unit.

Simulation of electronic dynamics in realistically large molecular systems is a demanding task that has not yet achieved the same level of quantitative prediction already realized for its static counterpart. This is particularly true for processes occurring beyond the Born–Oppenheimer regime. Non-adiabatic molecular dynamics (NAMD) simulations suffer from two convoluted sources of error: numerical algorithms for dynamics and electronic structure calculations. While the former has gained increasing attention, particularly addressing the validity of ad hoc methodologies, the effect of the latter remains relatively unexplored. Indeed, the required accuracy for electronic structure calculations to reach quantitative agreement with experiment in dynamics may be even more strict than that required for static simulations. Here, we address this issue by modeling the electronic energy transfer in a donor–acceptor–donor (D–A–D) molecular light harvesting system using fewest switches surface hopping NAMD simulations. In the studied system, time-resolved experimental measurements deliver complete information on spectra and energy transfer rates. Subsequent modeling shows that the calculated electronic transition energies are “sufficiently good” to reproduce experimental spectra but produce over an order of magnitude error in simulated dynamical rates. We further perform simulations using artificially shifted energy gaps to investigate the complex relationship between transition energies and modeled dynamics to understand factors affecting non-radiative relaxation and energy transfer rates.

Cyclic and bent conjugated molecular systems have tunable optical, structural, and dynamical features that differentiate them from their linear counterparts. Examples of such systems are [n]­cycloparaphenylenes (CPPs), which consist of nanorings composed of n para-linked benzene units. Circular geometry and tunability of π-orbital overlaps and bending strains enrich them with unique physicochemical and electronic properties compared to those of the corresponding linear oligoparaphenylenes. Herein, we explore the changes of these properties on alkyl-tethered-p-heptaphenylenes by modifying the methylene tether lengths from 1 to 19 carbons, leading to a gradual linearization of the conjugated backbone conformation. For this purpose, the photoinduced internal conversion processes of different alkyl-tethered-p-heptaphenylenes are simulated using nonadiabatic excited-state molecular dynamics. We found that the greater the strain introduced on the conjugated system, the slower the electronic and vibrational energy relaxation process. All bent p-heptaphenylenes exhibit similar patterns of intramolecular energy redistribution that finally spatially localize the exciton on phenylene units in the middle of the conjugated chain. This behavior is opposite to the random exciton localization previously reported for [n]­CPPs. Moreover, the nonadiabatic S2 → S1 electronic transition activates specific collective asymmetric vibrational excitations that promote periodic oscillatory evolution of the excitonic wave function before an excessive energy dissipates into the bath degrees of freedom.