Electronic Circular Dichroism in Exciton-Coupled Dimers: Vibronic Spectra from a General All-Coordinates Quantum-Dynamical Approach

Quantum-Classical Protocol for Efficient Characterization of Absorption Lineshape and Fluorescence Quenching upon Aggregation: The Case of Zinc Phthalocyanine Dyes

A quantum-classical protocol that incorporates Jahn-Teller vibronic coupling effects and cluster analysis of molecular dynamics simulations is reported, providing a tool for simulations of absorption spectra and ultrafast nonadiabatic dynamics in large molecular photosystems undergoing aggregation in solution. Employing zinc phthalocyanine dyes as target systems, we demonstrated that the proposed protocol provided fundamental information on vibronic, electronic couplings and thermal dynamical effects that mostly contribute to the absorption spectra lineshape and the fluorescence quenching processes upon dye aggregation. Decomposing the various effects arising upon dimer formation, the structure-property relations associated with their optical responses have been deciphered at atomistic resolution.

How the Interplay among Conformational Disorder, Solvation, Local, and Charge-Transfer Excitations Affects the Absorption Spectrum and Photoinduced Dynamics of Perylene Diimide Dimers: A Molecular Dynamics/Quantum Vibronic Approach

In this contribution we present a mixed quantum-classical dynamical approach for the computation of vibronic absorption spectra of molecular aggregates and their nonadiabatic dynamics, taking into account the coupling between local excitations (LE) and charge-transfer (CT) states. The approach is based on an adiabatic (Ad) separation between the soft degrees of freedom (DoFs) of the system and the stiff vibrations, which are described by the quantum dynamics (QD) of wave packets (WPs) moving on the coupled potential energy surfaces (PESs) of the LE and CT states. These PESs are described with a linear vibronic coupling (LVC) Hamiltonian, parameterized by an overlap-based diabatization on the grounds of time-dependent density functional theory computations. The WPs time evolution is computed with the multiconfiguration time-dependent Hartree method, using effective modes defined through a hierarchical representation of the LVC Hamiltonian. The soft DoFs are sampled with classical molecular dynamics (MD), and the coupling between the slow and fast DoFs is included by recomputing the key parameters of the LVC Hamiltonians, specifically for each MD configuration. This method, named Ad-MD|gLVC, is applied to a perylene diimide (PDI) dimer in acetonitrile and water solutions, and it is shown to accurately reproduce the change in the vibronic features of the absorption spectrum upon aggregation. Moreover, the microscopic insight offered by the MD trajectories allows for a detailed understanding of the role played by the fluctuation of the aggregate structure on the shape of the vibronic spectrum and on the population of LE and CT states. The nonadiabatic QD predicts an extremely fast (∼50 fs) energy transfer between the two LEs. CT states have only a moderate effect on the absorption spectrum, despite the fact that after photoexcitation they are shown to acquire a fast and non-negligible population, highlighting their relevance in dictating the charge separation and transport in PDI-based optical devices.

Bright Frenkel Excitons in Molecular Crystals: A Survey

We computed the optical properties of a large set of molecular crystals (∼2200 structures) composed of molecules whose lowest excited states are strongly coupled and generate wide excitonic bands. Such bands are classified in terms of their dimensionality (1-, 2-, and 3-dimensional), the position of the optically allowed state in relation with the excitonic density of states, and the presence of Davydov splitting. The survey confirms that one-dimensional aggregates are rare in molecular crystals highlighting the need to go beyond the simple low-dimensional models. Furthermore, this large set of data is used to search for technologically interesting and less common properties. For instance, we considered the largest excitonic bandwidth that is achievable within known molecular crystals and identified materials with strong super-radiant states. Finally, we explored the possibility that strong excitonic coupling can be used to generate emissive states in the near-infrared region in materials formed by molecules with bright visible absorption and we could identify the maximum allowable red shift in this material class. These insights with the associated searchable database provide practical guidelines for designing materials with interesting optical properties.

A computational study of the vibrationally-resolved electronic circular dichroism spectra of single-chain transoid and cisoid oligothiophenes in chiral conformations

We simulate the vibronic profile of the electronic circular dichroism (ECD) spectra of oligothiophenes in cisoid and transoid chiral arrangements. We consider oligomers of different lengths, from two to fifteen units, and investigate extensively how the ECD spectral shapes depend on the inter-ring torsions. In general, the molecular structures we consider are not stationary points of the ground state potential energy surface. Therefore, in order to perform vibronic calculations, we present a new computational protocol able to define reduced-dimensionality models where the effect of the off-equilibrium modes is removed. This is done adopting a description of the vibrational motions in curvilinear internal coordinates, and vertical harmonic models coupled with an iterative application of projectors to define energy Hessians, and therefore effective normal modes, in the space complementary to the one of the off-equilibrium coordinates. Although we consider both Franck-Condon and Herzberg-Teller contributions, the results show that transoid twisted ribbons always give rise to monosignated ECD spectra, while bi-signated and multi-signated spectra are expected for cisoid helices. These findings are explained on the basis of the different transition strengths of the lowest electronic states imparted by the different spatial arrangement, that is almost linear for transoid structures and more globular for cisoid ones. We predicted the chiroptical response of a large number of possible molecular arrangements. These data are employed to critically discuss the experimental ECD of polythiophenes in different experimental conditions, forming either aggregates or host-guest complexes. The method here proposed to perform vibronic calculations in reduced-dimensionality models is of general applicability and its potential interest goes beyond the practical application presented here.

Expanded Theory of H- and J-Molecular Aggregates: The Effects of Vibronic Coupling and Intermolecular Charge Transfer

The electronic excited states of molecular aggregates and their photophysical signatures have long fascinated spectroscopists and theoreticians alike since the advent of Frenkel exciton theory almost 90 years ago. The influence of molecular packing on basic optical probes like absorption and photoluminescence was originally worked out by Kasha for aggregates dominated by Coulombic intermolecular interactions, eventually leading to the classification of J- and H-aggregates. This review outlines advances made in understanding the relationship between aggregate structure and photophysics when vibronic coupling and intermolecular charge transfer are incorporated. An assortment of packing geometries is considered from the humble molecular dimer to more exotic structures including linear and bent aggregates, two-dimensional herringbone and "HJ" aggregates, and chiral aggregates. The interplay between long-range Coulomb coupling and short-range charge-transfer-mediated coupling strongly depends on the aggregate architecture leading to a wide array of photophysical behaviors.

Vibrational Structure in Magnetic Circular Dichroism Spectra of Polycyclic Aromatic Hydrocarbons

Absorption and magnetic circular dichroism (MCD) spectroscopies are powerful and simple methods to discriminate among various compounds. Polycyclic aromatic hydrocarbons provide particularly strong signal, which, for example, facilitates their detection in the environment. However, interpretation of the spectra is often based on quantum-chemical simulations, providing a limited precision only. In the present work, we use time-dependent density functional theory and complete active space second-order perturbation theories to understand spectral features observed in a series of naphthalene, anthracene, phenanthrene, and three larger compounds. The electronic computations provided reasonable agreement with the experiment for the smaller molecules, while a large error persisted for the bigger ones. However, many discrepancies could be explained by vibrational splitting of the electronic transitions across the entire spectral range. Compared to plain absorption, MCD spectral bands and their vibrational splitting were more specific for each aromatic molecule. The computational tools allowing simulations of detailed vibrational features in the electronic spectra thus promise to open a qualitatively new chapter in the spectroscopy of aromatic compounds.

The shape of the electronic circular dichroism spectrum of (2,6-dimethylphenyl)(phenyl)methanol: interplay between conformational equilibria and vibronic effects

Analysis of the interplay between conformational equilibria, solvent effects and vibronic contributions in the ECD spectra.

Circularly Polarized Luminescence from Axially Chiral BODIPY DYEmers: An Experimental and Computational Study

With our new home-built circularly polarized luminescence (CPL) instrument, we measured fluorescence and CPL spectra of the enantiomeric pairs of two quasi-isomeric BODIPY DYEmers 1 and 2, endowed with axial chirality. The electronic circular dichroism (ECD) and CPL spectra of these atropisomeric dimers are dominated by the exciton coupling between the main π-π* transitions (550-560 nm) of the two BODIPY rings. Compound 1 has strong ECD and CPL spectra (g_lum =4×10^-3 ) well reproduced by TD-DFT and SCS-CC2 (spin-component scaled second-order approximate coupled-cluster) calculations using DFT-optimized ground- and excited-state structures. Compound 2 has weaker ECD and CPL spectra (g_lum =4×10^-4 ), partly due to the mutual cancellation of electric-electric and electric-magnetic exciton couplings, and partly to its conformational freedom. This compound is computationally very challenging. Starting from the optimized excited-state geometries, we predicted the wrong sign for the CPL band of 2 using TD-DFT with the most recommended hybrid and range-separated functionals, whereas SCS-CC2 or a DFT functional with full exact exchange provided the correct sign.

Resolving Electronic Transitions in Synthetic Fluorescent Protein Chromophores by Magnetic Circular Dichroism

The detailed electronic structures of fluorescent chromophores are important for their use in imaging of living cells. A series of green fluorescent protein chromophore derivatives is examined by magnetic circular dichroism (MCD) spectroscopy, which allows the resolution of more bands than plain absorption and fluorescence. Observed spectral patterns are rationalized with the aid of time-dependent density functional theory (TDDFT) computations and the sum-over-state (SOS) formalism, which also reveals a significant dependence of MCD intensities on chromophore conformation. The combination of organic and theoretical chemistry with spectroscopic techniques also appears useful in the rational design of fluorescence labels and understanding of the chromophore's properties. For example, the absorption threshold can be heavily affected by substitution on the phenyl ring but not much on the five-member ring, and methoxy groups can be used to further tune the electronic levels.

A simple dimeric model accounts for the vibronic ECD spectra of chiral polythiophenes in their aggregated states

Aggregates of chiral polythiophenes (PTs) show exciton-coupled electronic circular dichroism (ECD) spectra with a unique vibronic structure, which can be reproduced by quantum dynamical simulations on the coupled states of small dimeric models.

Circularly polarized luminescence spectroscopy reveals low-energy excited states and dynamic localization of vibronic transitions in CP43

Circularly polarized luminescence (CPL) spectroscopy is an established but relatively little-used technique that monitors the chirality of an emission. When applied to photosynthetic pigment assemblies, we find that CPL provides sensitive and detailed information on low-energy exciton states, reflecting the interactions, site energies and geometries of interacting pigments. CPL is the emission analog of circular dichroism (CD) and thus spectra explore the optical activity only of fluorescent states of the pigment-protein complex and consequently the nature of the lowest-energy excited states (trap states), whose study is a critical area of photosynthesis research. In this work, we develop the new approach of temperature-dependent CPL spectroscopy, over the 2-120 K temperature range, and apply it to the CP43 proximal antenna protein of photosystem II. Our results confirm strong excitonic interactions for at least one of the two well-established emitting states of CP43 named "A" and "B". Previous structure-based models of CP43 spectra are evaluated in the light of the new CPL data. Our analysis supports the assignments of Shibata et al. [Shibata et al. J. Am. Chem. Soc. 135 (2013) 6903-6914], particularly for the highly-delocalized B-state. This state dominates CPL spectra and is attributed predominantly to chlorophyll a's labeled Chl 634 and Chl 636 (alternatively labeled Chl 43 and 45 by Shibata et al.). The absence of any CPL intensity in intramolecular vibrational sidebands associated with the delocalized "B" excited state is attributed to the dynamic localization of intramolecular vibronic transitions.

Probing deactivation pathways of DNA nucleobases by two-dimensional electronic spectroscopy: first principles simulations

The SOS//QM/MM [Rivalta et al., Int. J. Quant. Chem., 2014, 114, 85] method consists of an arsenal of computational tools allowing accurate simulation of one-dimensional (1D) and bi-dimensional (2D) electronic spectra of monomeric and dimeric systems with unprecedented details and accuracy. Prominent features like doubly excited local and excimer states, accessible in multi-photon processes, as well as charge-transfer states arise naturally through the fully quantum-mechanical description of the aggregates. In this contribution the SOS//QM/MM approach is extended to simulate time-resolved 2D spectra that can be used to characterize ultrafast excited state relaxation dynamics with atomistic details. We demonstrate how critical structures on the excited state potential energy surface, obtained through state-of-the-art quantum chemical computations, can be used as snapshots of the excited state relaxation dynamics to generate spectral fingerprints for different de-excitation channels. The approach is based on high-level multi-configurational wavefunction methods combined with non-linear response theory and incorporates the effects of the solvent/environment through hybrid quantum mechanics/molecular mechanics (QM/MM) techniques. Specifically, the protocol makes use of the second-order Perturbation Theory (CASPT2) on top of Complete Active Space Self Consistent Field (CASSCF) strategy to compute the high-lying excited states that can be accessed in different 2D experimental setups. As an example, the photophysics of the stacked adenine-adenine dimer in a double-stranded DNA is modeled through 2D near-ultraviolet (NUV) spectroscopy.

Determination of the absolute configuration of perylene quinone-derived mycotoxins by measurement and calculation of electronic circular dichroism spectra and specific rotations

Altertoxins I-III, alterlosins I and II, alteichin (alterperylenol), stemphyltoxins I-IV, stemphyperylenol, stemphytriol, 7-epi-8-hydroxyaltertoxin I, and 6-epi-stemphytriol are mycotoxins derived from perylene quinone, for which the absolute configuration was not known. Electronic circular dichroism (ECD) spectra were calculated for these compounds and compared with measured spectra of altertoxins I-III, alteichin, and stemphyltoxin III and with reported Cotton effects. Specific rotations were calculated and compared with reported specific rotations. The absolute configuration of all the toxins, except for stemphyltoxin IV, could thus be determined. The validity of the assignment was high whenever reported ECD data were available for comparison, and the validity was lower when the assignment was based only on the comparison of calculated and reported specific rotations. ECD spectra are intrinsically different for toxins with a biphenyl substructure and for toxins derived from dihydroanthracene.