Adiabatic One- and Two-Photon Excited States in Phenylene-Based Conjugated Oligomers: A Quantum-Chemical Study

Controlling Intramolecular Singlet Fission Dynamics via Torsional Modulation of Through-Bond versus Through-Space Couplings

Covalent dimers, particularly pentacenes, are the dominant platform for developing a mechanistic understanding of intramolecular singlet fission (iSF). Numerous studies have demonstrated that a photoexcited singlet state in these structures can rapidly and efficiently undergo exciton multiplication to form a correlated pair of triplets within a single molecule, with potential applications from photovoltaics to quantum information science. One of the most significant barriers limiting such dimers is the fast recombination of the triplet pair, which prevents spatial separation and the formation of long-lived triplet states. There is an ever-growing need to develop general synthetic strategies to control the evolution of triplets following iSF and enhance their lifetime. Here, we rationally tune the dihedral angle and interchromophore separation between pairs of pentacenes in a systematic series of bridging units to facilitate triplet separation. Through a combination of transient optical and spin-resonance techniques, we demonstrate that torsion within the linker provides a simple synthetic handle to tune the fine balance between through-bond and through-space interchromophore couplings that steer iSF. We show that the full iSF pathway from femtosecond to microsecond timescales is tuned through the static coupling set by molecular design and structural fluctuations that can be biased through steric control. Our approach highlights a straightforward design principle to generate paramagnetic spin pair states with higher yields.

Adiabatic Excited-State Molecular Dynamics with an Explicit Solvent: NEXMD-SANDER Implementation

We present a method to link the Nonadiabatic EXcited-state Molecular Dynamics (NEXMD) package to the SANDER package supplied by AMBERTOOLS to provide excited-state adiabatic quantum mechanics/molecular mechanics (QM/MM) simulations. NEXMD is a computational package particularly developed to perform simulations of the photoexcitation and subsequent nonadiabatic electronic and vibrational energy relaxation in large multichromophoric conjugated molecules involving several coupled electronic excited states. The NEXMD-SANDER exchange has been optimized in order to achieve excited-state adiabatic dynamics simulations of large conjugated materials in a QM/MM environment, such as an explicit solvent. Dynamics of a substituted polyphenylene vinylene oligomer (PPV3-NO₂) in vacuum and different explicit solvents has been used as a test case by performing comparative analysis of changes in its optical spectrum, state-dependent conformational changes, and quantum bond orderings. The method has been tested and compared with respect to previous implicit solvent implementations. Also, the impact on the expansion of the QM region by including a variable number of solvent molecules has been analyzed. Altogether, these results encourage future implementations of NEXMD simulations using the same combination of methods.

Excited-state structural relaxation and exciton delocalization dynamics in linear and cyclic π-conjugated oligothiophenes

π-Conjugated oligothiophene is considered a chain segment of its polymeric counterpart with simper excited-state dynamics and spectral signatures.

The Role of Linkers in the Excited-State Dynamic Planarization Processes of Macrocyclic Oligothiophene 12-Mers

Linkers adjoining chromophores play an important role in modulating the structure of conjugated systems, which is bound up with their photophysical properties. However, to date, the focus of works dealing with linker effects was limited only to linear π-conjugated materials, and there have been no detailed studies on cyclic counterparts. Herein we report the linker effects on the dynamic planarization processes of π-conjugated macrocyclic oligothiophene 12-mers, where the different ratio between ethynylene and vinylene linkers was chosen to control the backbone rigidity. By analyzing transient fluorescence spectra, we demonstrate that the connecting linkers play a crucial role in the excited-state dynamics of cyclic conjugated systems. Faster dynamic planarization, longer exciton delocalization length, and higher degree of planarity were observed in vinylene inserted cyclic oligothiophenes. Molecular dynamics simulations and density functional theory calculations also stress the importance of the role of linkers in modulating the structure of cyclic oligothiophenes.

Nonadiabatic excited-state molecular dynamics: modeling photophysics in organic conjugated materials

To design functional photoactive materials for a variety of technological applications, researchers need to understand their electronic properties in detail and have ways to control their photoinduced pathways. When excited by photons of light, organic conjugated materials (OCMs) show dynamics that are often characterized by large nonadiabatic (NA) couplings between multiple excited states through a breakdown of the Born-Oppenheimer (BO) approximation. Following photoexcitation, various nonradiative intraband relaxation pathways can lead to a number of complex processes. Therefore, computational simulation of nonadiabatic molecular dynamics is an indispensable tool for understanding complex photoinduced processes such as internal conversion, energy transfer, charge separation, and spatial localization of excitons. Over the years, we have developed a nonadiabatic excited-state molecular dynamics (NA-ESMD) framework that efficiently and accurately describes photoinduced phenomena in extended conjugated molecular systems. We use the fewest-switches surface hopping (FSSH) algorithm to treat quantum transitions among multiple adiabatic excited state potential energy surfaces (PESs). Extended molecular systems often contain hundreds of atoms and involve large densities of excited states that participate in the photoinduced dynamics. We can achieve an accurate description of the multiple excited states using the configuration interaction single (CIS) formalism with a semiempirical model Hamiltonian. Analytical techniques allow the trajectory to be propagated "on the fly" using the complete set of NA coupling terms and remove computational bottlenecks in the evaluation of excited-state gradients and NA couplings. Furthermore, the use of state-specific gradients for propagation of nuclei on the native excited-state PES eliminates the need for simplifications such as the classical path approximation (CPA), which only uses ground-state gradients. Thus, the NA-ESMD methodology offers a computationally tractable route for simulating hundreds of atoms on ~10 ps time scales where multiple coupled excited states are involved. In this Account, we review recent developments in the NA-ESMD modeling of photoinduced dynamics in extended conjugated molecules involving multiple coupled electronic states. We have successfully applied the outlined NA-ESMD framework to study ultrafast conformational planarization in polyfluorenes where the rate of torsional relaxation can be controlled based on the initial excitation. With the addition of the state reassignment algorithm to identify instances of unavoided crossings between noninteracting PESs, NA-ESMD can now be used to study systems in which these so-called trivial unavoided crossings are expected to predominate. We employ this technique to analyze the energy transfer between poly(phenylene vinylene) (PPV) segments where conformational fluctuations give rise to numerous instances of unavoided crossings leading to multiple pathways and complex energy transfer dynamics that cannot be described using a simple Förster model. In addition, we have investigated the mechanism of ultrafast unidirectional energy transfer in dendrimers composed of poly(phenylene ethynylene) (PPE) chromophores and have demonstrated that differential nuclear motion favors downhill energy transfer in dendrimers. The use of native excited-state gradients allows us to observe this feature.

Experimental and theoretical study of the n-doped successive polyanions of oligocruciform molecular wires: up to five units of charge

The electronic structure of polyanions of sterically encumbered triisopropylsilyl-substituted linear and cyclic oligo(phenyleneethynylene)s (Monomer, Trimer, Pentamer, and Triangle) is investigated by electron paramagnetic resonance (EPR), electron nuclear double resonance (ENDOR), and UV/Vis-near-infrared (NIR) spectroscopies, cyclic voltammetry, and theoretical calculations (DFT). Increasing anion orders are generated sequentially in vacuo at room temperature by chemical reaction with potassium metal up to the pentaanion. The relevance of these compounds acting as electron reservoirs is thus demonstrated. Even-order anions are EPR silent, whereas the odd species exhibit different signatures, which are identified after comparison of the measured hyperfine couplings by ENDOR spectroscopy with those predicted by DFT calculations. With increasing size of the oligomers the electron spin density is first distributed over the backbone carbon atoms for the monoanions, and then further localized at the outer phenyl rings for the trianion species. Examination of the UV/Vis-NIR spectra indicates that the monoanions (T(.-) , P(.-) ) exhibit two transitions in the Vis-NIR region, whereas a strong absorption in the IR region is solely observed for higher reduced states. Electronic transitions of the neutral monoanions and trianions are redshifted with increasing oligomer size, whereas for a given oligomer a blueshift is observed upon increasing the charge, which suggests a localization of the spin density.

Conformational disorder and ultrafast exciton relaxation in PPV-family conjugated polymers

We report combined experimental and theoretical studies of excitation relaxation in poly[2-methoxy,5-(2'-ethyl-hexoxy)-1,4-phenylenevinylene] (MEH-PPV), oligophenylenevinylene (OPV) molecules of varying length, and model PPV chains. We build on the paradigm that the basic characteristics of conjugated polymers are decided by conformational subunits defined by conjugation breaks caused by torsional disorder along the chain. The calculations reported here indicate that for conjugated polymers like those in the PPV family, these conformational subunits electronically couple to neighboring subunits, forming subtly delocalized collective states of nanoscale excitons that determine the polymer optical properties. We find that relaxation among these exciton states can lead to a decay of anisotropy on ultrafast time scales. Unlike in Forster energy transfer, the exciton does not necessarily translate over a large distance. Nonetheless, the disorder in the polymer chain means that even small changes in the exciton size or location has a significant effect on the relaxation pathway and therefore the anisotropy decay.

Theoretical Study of Spectroscopic Properties of Dimethoxy-p-Phenylene-Ethynylene Oligomers: Planarization of the Conjugated Backbone

The optical spectra of the dimethoxy-p-phenylene-ethynylene oligomers (up to n = 10) are calculated by DFT and TD-DFT methods. It is found that the conformational rotations around the cylindrical triple-bonded carbon links impact significantly the optical spectrum. The effective conjugation length (ECL) of the oligomer is obtained by extrapolating the HOMO-LUMO gap to infinite chain length with an alternative exponential function. The spectral shift is mainly dependent on the high pi-conjugation segment of oligomers, resulting from the planarization of the backbone. Although the rotational barrier is very low, the calculated results further indicate that rotation about the cylindrical triple bond still interrupts the conjugation of rod-like oligomers to some extent, and leads to an angle-dependent HOMO-LUMO gap. The results are helpful to interpret the conformational-dependent spectroscopic phenomena of p-phenyleneethynylene oligomers and polymers (PPEs) observed in ensemble and single molecule spectroscopy.

Investigation of the Electronic Spectra and Excited-State Geometries of Poly(para-phenylene vinylene) (PPV) and Poly(para-phenylene) (PP) by the Symmetry-Adapted Cluster Configuration Interaction (SAC-CI) Method

The symmetry-adapted cluster-configuration interaction (SAC-CI) method has been used to investigate the optical and geometric properties of the oligomers of poly(para-phenylene vinylene) (PPV) and poly(para-phenylene) (PP). Vertical singlet and triplet absorption spectra and emission spectra have been calculated accurately; the mean average deviation from available experimental results lies within 0.2 eV. The chain length dependence of the transition energies has been improved in comparison to earlier TDDFT and MRSDCI calculations. The present analysis suggests that conventional TDDFT with the B3LYP functional should be used carefully, as it can provide inaccurate estimates of the chain length dependence of the excitation energies of these molecules with long pi conjugation. The T1 state was predicted to be at a lower energy, by 1.0-1.5 eV for PPV and by 0.9-1.7 eV for PP, than the S1 state, which indicates a localized T1 state with large exchange energy. By calculating the SAC-CI electron density difference between the ground and excited states, the geometry relaxations due to excitations can be analyzed in detail using electrostatic force theory. For trans-stilbene, the doubly excited 21Ag state was studied, and the calculated transition energy of 4.99 eV agrees very well with the experimental value of 4.84 eV. In contrast to previous ab initio calculations, we predict this doubly excited 21Ag state to lie above the 11Bu state.

Quantum Simulation of Solution Phase Intramolecular Electron Transfer Rates in Betaine-30

Mixed quantum-classical atomistic simulations have been carried out to investigate the mechanistic details of excited state intramolecular electron transfer in a betaine-30 molecule in acetonitrile. The key electronic degrees of freedom of the solute molecule are treated quantum mechanically using the semiempirical Pariser-Parr-Pople Hamiltonian, including the solvent influence on electronic structure. The intramolecular vibrational modes are also treated explicitly at a quantum level, with the remaining elements treated classically using empirical potentials. The electron-transfer rate, corresponding to S1 --> S0 relaxation, is evaluated via time-dependent perturbation theory with the explicit inclusion of the dynamics of solvation and intramolecular conformation. The calculations reveal that, while solvation dynamics is critical to the rate, the intramolecular torsional dynamics also plays an important role. The importance of the use of multiple high-frequency quantum modes is also discussed.

Modeling the Effects of Torsional Disorder on the Spectra of Poly- and Oligo-(p-phenyleneethynylenes)

The absorption spectra of phenyleneethynylene oligomers show an unusual change in shape with oligomer length. The unusual aspects of the spectra arise from rotation of the phenylene rings about the long axis of the oligomer. In the ground electronic state, the barrier to this rotation is low and the spectra in room temperature come from an ensemble of different structures. In the excited electronic state, the barrier to rotation is substantially higher, giving rise to strong nonlinear electron-phonon coupling. A multidimensional semiempirical model that includes these effects is developed for the photophysics of phenyleneethynylene oligomers. The ground-state energy is modeled with a molecular mechanics expression, and the excitation energy is modeled with an exciton model. Intermediate Neglect of Differential Overlap/Singles Configuration Interaction (INDO/SCI) calculations verify the exciton model and provide initial estimates of the model parameters. These parameters generate the qualitative features seen in experimental spectra. Inclusion of entropy effects from the multiple torsional coordinates is essential. Refinement of the parameters yields quantitative agreement with experiment. This agreement shows that coupling to torsional motion is a major factor in the spectroscopy and photophysics of these conjugated polymers.

Synthesis and Optical Properties of Bis(oligophenyleneethynylenes)

Bis(oligophenyleneethynylenes) 1-4 were prepared as representative members of a new class of potential nonlinear optical materials. The optical properties of 1-4 were examined for evidence of restricted rotation of the aryl rings when compared to their single-strand precursors, which could potentially increase their nonlinear response through more effective conjugation. The effect of altering the electron density of the terminating functional group of these compounds on their properties was also investigated.

Optical Excitations in Carbon Architectures Based on Dodecadehydrotribenzo[18]annulene

The origin of excitations in multi-chromophore carbon network substructures based on dodecadehydrotribenzo[18]annulene has been investigated by steady-state and photon echo spectroscopy, configuration interaction (CIS and CIS(D)), and time-dependent density functional theory (TD-DFT). 1,4-diphenylbutadiyne, the simplest structural subunit within the annulene, was used in modeling the spectroscopic studies to explain the origin of excitations in the macrocycles. The optical excitations in longer linear systems were found to be similar to its diphenylacetylene analogue. However, the results from dodecadehydrotribenzo[18]annulene and other multichromophore networks systems illustrate the possibility of strong intramolecular interactions and the formation of delocalized excited states. Calculations were carried out to explain the basic similarities and differences in excitations of the model compounds such as diphenylbutadiyne and the macrocycles. The fundamental excitation in these systems can be primarily described as a pi --> pi* transition. Two low-energy resonances were observed from experiment for the annulene systems, and possible explanations for these low-energy resonances in the macrocycles are explored. The significant difference found in the calculated oscillator strength of the two low-energy bands for the macrocycles as well as the dynamics of solvent interactions was further investigated by three-pulse photon echo measurements. A simple exciton model was developed to discuss the excitations in the larger macrocycles. The results from this model were found to be in good agreement with the TD-DFT calculations.

Time-dependent density-functional calculations of S[sub 0]–S[sub 1] transition energies of poly(p-phenylene vinylene)

We used time-dependent density-functional-theory (TDDFT) approaches to calculate absorption (S(0)-->S(1)) and emission (S(1)-->S(0)) transition energies of poly(p-phenylene vinylene) (PPV). The absorption and emission energies were estimated to be 2.44 and 2.16 eV, respectively, from the extrapolation of calculated results for oligomers. Comparisons with available experimental data demonstrated that TDDFT is a very reliable tool for investigating the electronic transitions of PPV.

Spin-dependent electron-hole capture kinetics in luminescent conjugated polymers

The recombination of electron-hole pairs injected in extended conjugated systems is modeled as a multistep interconversion relaxation in monoexcited electronic state space, mediated by electron-phonon coupling. The computed ratio of triplet-to-singlet exciton formation times r=tau(T)/tau(S) increases from 0.9 for a model dimer to 2.5 for a 32-unit chain, in good agreement with recent experiments. We rationalize the conjugation-length dependence of r in terms of spin-specific energetics and mutual vibronic coupling of the excited states.

Conformational dynamics of photoexcited conjugated molecules

Time-dependent photoexcitation and optical spectroscopy of pi-conjugated molecules is described using a new method for the simulation of excited state molecular dynamics in extended molecular systems with sizes up to hundreds of atoms. Applications are made to poly(p-phenylene vinylene) oligomers. Our analysis shows self-trapping of excitations on about six repeat units in the course of photoexcitation relaxation, identifies specific slow (torsion) and fast (bond-stretch) nuclear motions strongly coupled to the electronic degrees of freedom, and predicts spectroscopic signatures of molecular conformations.