Multiscale modeling of molecular structure and optical properties of complex supramolecular aggregates

Supramolecular aggregates of synthetic dye molecules offer great perspectives to prepare biomimetic functional materials for light-harvesting and energy transport. The design is complicated by the fact that structure–property relationships are hard to establish, because the molecular packing results from a delicate balance of interactions and the excitonic properties that dictate the optics and excited state dynamics, in turn sensitively depend on this packing. Here we show how an iterative multiscale approach combining molecular dynamics and quantum mechanical exciton modeling can be used to obtain accurate insight into the packing of thousands of cyanine dye molecules in a complex double-walled tubular aggregate in close interaction with its solvent environment. Our approach allows us to answer open questions not only on the structure of these prototypical aggregates, but also about their molecular-scale structural and energetic heterogeneity, as well as on the microscopic origin of their photophysical properties. This opens the route to accurate predictions of energy transport and other functional properties.

Multiscale modeling resolves the molecular structure of a synthetic light-harvesting complex, unraveling the microscopic origin of its photophysical properties.

Supramolecular structures may self-assemble from a variety of building blocks, resulting in a wide range of advanced materials with attractive biomimetic, sensing, catalytic, optoelectronic and photonic functionalities.^1–10 The close-packed nanoscale organization of the individual molecules within a supramolecular system, held together via noncovalent interactions, gives rise to the aggregate''s (collective) properties. Assemblies consisting of dye molecules often exhibit unique collective optical properties and are of interest for opto-electronic applications as well as artificial light-harvesting complexes that mimic natural antenna systems of photosynthetic bacteria and plants.^11–13 For example, chlorosomal antenna complexes of photosynthetic green sulfur bacteria are self-assembled into multilayer tubular structures having bacteriochlorophyll pigments as building blocks.^14–16 The structure of these antenna complexes and the underlying molecular arrangement ensures that the process of light-harvesting and excitation energy transport is very efficient, even under extremely low light conditions.^17,18 The quest to recreate such efficiency under laboratory conditions has sparked numerous studies of synthetic self-assembled systems mimicking natural chlorosomes, e.g. using porphyrins,¹⁹ zinc chlorin,²⁰ and cyanine dyes.²¹ Of particular interest are the tubular aggregates of 3,3′-bis(2-sulfopropyl)-5,5′,6,6′-tetrachloro-1,1′-dioctylbenzimidacarbocyanine (C8S3).^22–25 Cryo-TEM reveals a hierarchy of supramolecular architectures, including double-walled nanotubes; under certain conditions, bundles of nanotubes arise.²⁶ Thus, this system allows for the occurrence of electronic excitation energy transport at various levels: within one wall, between walls of one tube, and between different tubes, similar to the situation in natural systems.^27,28To understand how such supramolecular systems work, as well as propose design rules for new materials, it is essential to determine the relationship between molecular structure and optical properties. Current experimental techniques, however, are unable to resolve the structure at the molecular level. This, in combination with the sensitivity of spectral properties to the details of the molecular packing, leads to a crucial role for theoretical modeling.²⁹ For example, molecular dynamics (MD) simulations have been used to predict the molecular packing within a variety of supramolecular assemblies.^30–34 However, synthetic amphiphiles with aromatic groups, such as cyanine dyes—often used to prepare aggregates with optical functionality—tend to fall into kinetic traps during spontaneous self-assembly simulations and the packing of the aromatic chromophores remains highly disordered on the accessible time scale, leading to predicted (optical) spectra that are not consistent with experimental data.³⁵ This problem can be overcome by building assemblies based upon proposed architectures and assessing their stability in relatively short MD simulations.^36–38 The drawback of this approach is the requirement of a thorough understanding of what to use as a starting point and how to validate the structure. In any case, proper validation requires the modeling of the optical spectra of the obtained structure, and finally, comparing it to the experiment. The demanding character of such methods explains why an important role is played by phenomenological modeling, in which a molecular packing is guessed and the optics is obtained from parametrizing an exciton model that describes the collective excited states of the assembly with interactions dictated by the guessed packing. By comparing the calculated spectra to experimental ones, the structure and exciton model may be fine-tuned. While this method has been successful in describing spectra,^23,39 it is limited in its predictive power and also lacks access to essential microscopic parameters, such as tuning of the optical excitation energies imposed by the environment, disorder in these energies and structural heterogeneity.In this work, we use an advanced multiscale approach to determine structure–optical property relationships for the C8S3 double-walled nanotubes, guided by comparison to experiments. The optical spectrum of these aggregates, in which multiple exciton peaks may be discerned, suggests a rather complex underlying molecular packing. This fact, combined with their sheer size going up to many thousands of molecules, makes these systems exceptionally challenging to resolve and leaves important questions concerning structure–function relationships unanswered or under debate, for instance the origin of the splitting between the two lowest-energy spectral bands.^23,38 Here, we answer these questions by iteratively combining MD simulations to capture the details of molecular packing and structural disorder, an exciton Hamiltonian approach to calculate optical signatures, and explicit microelectrostatic calculations to estimate energetic disorder and solvent shifts. Previous attempts to reveal the structure of cyanine-based nanotubes were limited to small-scale system sizes,^37,38 modeling optical features phenomenologically rather than using atomistic information³⁸ or featuring simpler, single-walled systems.³⁷ In addition to answering important questions for the C8S3 double-walled nanotubes, our study opens the way to explain and predict at an unprecedented level of detail the functional properties of other highly complex molecular materials.