Vibronic coupling in molecular crystals: A Franck-Condon Herzberg-Teller model of H-aggregate fluorescence based on quantum chemical cluster calculations.

Here, we present a general approach to treating vibronic coupling in molecular crystals based on atomistic simulations of large clusters. Such clusters comprise model aggregates treated at the quantum chemical level embedded within a realistic environment treated at the molecular mechanics level. As we calculate ground and excited state equilibrium geometries and vibrational modes of model aggregates, our approach is able to capture effects arising from coupling to intermolecular degrees of freedom, absent from existing models relying on geometries and normal modes of single molecules. Using the geometries and vibrational modes of clusters, we are able to simulate the fluorescence spectra of aggregates for which the lowest excited state bears negligible oscillator strength (as is the case, e.g., ideal H-aggregates) by including both Franck-Condon (FC) and Herzberg-Teller (HT) vibronic transitions. The latter terms allow the adiabatic excited state of the cluster to couple with vibrations in a perturbative fashion via derivatives of the transition dipole moment along nuclear coordinates. While vibronic coupling simulations employing FC and HT terms are well established for single-molecules, to our knowledge this is the first time they are applied to molecular aggregates. Here, we apply this approach to the simulation of the low-temperature fluorescence spectrum of para-distyrylbenzene single-crystal H-aggregates and draw comparisons with coarse-grained Frenkel-Holstein approaches previously extensively applied to such systems.

Dynamics of guest molecules in PHTP inclusion compounds as probed by solid-state NMR and fluorescence spectroscopy.

Partially deuterated 1,4-distyrylbenzene () is included into the pseudohexagonal nanochannels of perhydrotriphenylene (PHTP). The overall and intramolecular mobility of is investigated over a wide temperature range by (13)C, (2)H NMR as well as fluorescence spectroscopy. Simulations of the (2)H NMR spectral shapes reveal an overall wobble motion of in the channels with an amplitude of about 4 degrees at T = 220 K and 10 degrees at T = 410 K. Above T = 320 K the wobble motion is superimposed by localized 180 degrees flips of the terminal phenyl rings with a frequency of 10(6) Hz at T = 340 K. The activation energies of both types of motions are around 40 kJ mol(-1) which imply a strong sterical hindrance by the surrounding PHTP channels. The experimental vibrational structure of the fluorescence excitation spectra of is analyzed in terms of small amplitude ring torsional motions, which provide information about the spatial constraints on by the surrounding PHTP host matrix. Combining the results from NMR and fluorescence spectroscopy as well as of time-dependent density functional calculations yields the complete potential surfaces of the phenyl ring torsions. These results, which suggest that intramolecular mobility of is only reduced but not completely suppressed by the matrix, are corroborated by MD simulations. Unrealistically high potential barriers for phenyl ring flips are obtained from MD simulations using rigid PHTP matrices which demonstrate the importance of large amplitude motions of the PHTP host lattice for the mobility of the guest molecules.