Excitonic Hamiltonians for Calculating Optical Absorption Spectra and Optoelectronic Properties of Molecular Aggregates and Solids.

Rational design of disordered molecular aggregates and solids for optoelectronic applications relies on our ability to predict the properties of such materials using theoretical and computational methods. However, large molecular systems where disorder is too significant to be considered in the perturbative limit cannot be described using either first principles quantum chemistry or band theory. Multiscale modeling is a promising approach to understanding and optimizing the optoelectronic properties of such systems. It uses first-principles quantum chemical methods to calculate the properties of individual molecules, then constructs model Hamiltonians of molecular aggregates or bulk materials based on these calculations. In this paper, we present a protocol for constructing a tight-binding Hamiltonian that represents the excited states of a molecular material in the basis of Frenckel excitons: electron-hole pairs that are localized on individual molecules that make up the material. The Hamiltonian parametrization proposed here accounts for excitonic couplings between molecules, as well as for electrostatic polarization of the electron density on a molecule by the charge distribution on surrounding molecules. Such model Hamiltonians can be used to calculate optical absorption spectra and other optoelectronic properties of molecular aggregates and solids.

Molecular Mechanics Simulations and Improved Tight-Binding Hamiltonians for Artificial Light Harvesting Systems: Predicting Geometric Distributions, Disorder, and Spectroscopy of Chromophores in a Protein Environment.

We present molecular mechanics and spectroscopic calculations on prototype artificial light harvesting systems consisting of chromophores attached to a tobacco mosaic virus (TMV) protein scaffold. These systems have been synthesized and characterized spectroscopically, but information about the microscopic configurations and geometry of these TMV-templated chromophore assemblies is largely unknown. We use a Monte Carlo conformational search algorithm to determine the preferred positions and orientations of two chromophores, Coumarin 343 together with its linker and Oregon Green 488, when these are attached at two different sites (104 and 123) on the TMV protein. The resulting geometric information shows that the extent of disorder and aggregation properties and therefore the optical properties of the TMV-templated chromophore assembly are highly dependent on both the choice of chromophores and the protein site to which they are bound. We use the results of the conformational search as geometric parameters together with an improved tight-binding Hamiltonian to simulate the linear absorption spectra and compare with experimental spectral measurements. The ideal dipole approximation to the Hamiltonian is not valid because the distance between chromophores can be very small. We found that using the geometries from the conformational search is necessary to reproduce the features of the experimental spectral peaks.

Absorption Spectra for Disordered Aggregates of Chromophores Using the Exciton Model.

Optimizing the optical properties of large chromophore aggregates and molecular solids for applications in photovoltaics and nonlinear optics is an outstanding challenge. It requires efficient and reliable computational models that must be validated against accurate theoretical methods. We show that linear absorption spectra calculated using the molecular exciton model agree well with spectra calculated using time-dependent density functional theory and configuration interaction singles for aggregates of strongly polar chromophores. Similar agreement is obtained for a hybrid functional (B3LYP), a long-range corrected hybrid functional (ωB97X), and configuration interaction singles. Accounting for the electrostatic environment of individual chromophores in the parametrization of the exciton model with the inclusion of atomic point charges significantly improves the agreement of the resulting spectra with those calculated using all-electron methods; different charge definitions (Mulliken and ChelpG) yield similar results. We find that there is a size-dependent error in the exciton model compared with all-electron methods, but for aggregates with more than six chromophores, the errors change slowly with the number of chromophores in the aggregate. Our results validate the use of the molecular exciton model for predicting the absorption spectra of bulk molecular solids; its formalism also allows straightforward extension to calculations of nonlinear optical response.

Single molecule charge transport: from a quantum mechanical to a classical description.

This paper explores charge transport at the single molecule level. The conductive properties of both small organic molecules and conjugated polymers (molecular wires) are considered. In particular, the reasons for the transition from fully coherent to incoherent charge transport and the approaches that can be taken to describe this transition are addressed in some detail. The effects of molecular orbital symmetry, quantum interference, static disorder and molecular vibrations on charge transport are discussed. All of these effects must be taken into account (and may be used in a functional way) in the design of molecular electronic devices. An overview of the theoretical models employed when studying charge transport in small organic molecules and molecular wires is presented.

Mechanism of charge transport along zinc porphyrin-based molecular wires.

In this study charge transport along zinc porphyrin-based molecular wires is simulated, considering both bandlike and hopping mechanisms. It is shown that bandlike transport simulations yield significantly overestimated hole mobility values. On the basis of kinetic and thermodynamic considerations, it is inferred that charge transport along zinc porphyrin-based molecular wires occurs by small polaron hopping. Hole mobility values on the order of 0.1 cm(2) V(-1) s(-1) are found from small polaron hopping simulations, which agree well with previously reported experimental results. It is suggested that the experimentally observed increase of the charge carrier mobility on formation of supramolecular ladderlike structures is determined by two factors. One of these is an increase of charge transfer integrals between monomer units due to molecular wire planarization. A more important factor is the reduction of the amount of energetic disorder along the molecular wire and in its environment. General guidelines for determining the mechanism of charge transport along molecular wires are discussed.