Dispersion-Driven Cooperativity in Alkyl Perylene Diimide Oligomers: Insights from Density Functional Theory

The cooperative mechanism is of paramount importance in the synthesis of supramolecular polymers with desired characteristics, including molecular mass, polydispersity, and morphology. It is primarily driven by the presence of intermolecular interactions, which encompass strong hydrogen bonding, metal-ligand interactions, and dipole-dipole interactions. In this study, we utilize density functional theory and energy decomposition analysis to investigate the cooperative behavior of perylene diimide (PDI) oligomers with alkyl chains at their imide positions, which lack the previously mentioned interactions. Our systematic examination reveals that dispersion interactions originating from the alkyl side-chain substituents play an important role in promoting cooperativity within these PDIs. This influence becomes even more pronounced for alkyl chain lengths beyond hexyl groups. The energy decomposition analysis reveals that the delicate balance between dispersion energy and Pauli repulsion energy is the key driver of cooperative behavior in PDIs. Additionally, we have developed a mathematical model capable of predicting the saturated binding energies for PDI oligomers of varying sizes and alkyl chain lengths. Overall, our findings emphasize the previously undervalued significance of dispersion forces in cooperative supramolecular polymerization, enhancing our overall understanding of the cooperative mechanism.

Breaking Radial Dipole Symmetry in Planar Macrocycles Modulates Edge-to-Edge Packing and Disrupts Cofacial Stacking

Dipolar interactions are ever-present in supramolecular architectures, though their impact is typically revealed by making dipoles stronger. While it is also possible to assess the role of dipoles by altering their orientations by using synthetic design, doing so without altering the molecular shape is not straightforward. We have now done this by flipping one triazole unit in a rigid macrocycle, tricarb. The macrocycle is composed of three carbazoles (2 Debye) and three triazoles (5 Debye) defining an array of dipoles aligned radially but organized alternately in and out. These dipoles are believed to dictate edge-to-edge tiling and face-to-face stacking. We modified our synthesis to prepare isosteric macrocycles with the orientation of one triazole dipole rotated 40°. The new dipole orientation guides edge-to-edge contacts to reorder the stability of two surface-bound 2D polymorphs. The impact on dipole-enhanced π stacking, however, was unexpected. Our stacking model identified an unchanged set of short-range (3.4 Å) anti-parallel dipole contacts. Despite this situation, the reduction in self-association was attributed to long-range (~6.4 Å) dipolar repulsions between π-stacked macrocycles. This work highlights our ability to control the build-up and symmetry of macrocyclic skeletons by synthetic design, and the work needed to further our understanding of how dipoles control self-assembly.

Molecular mechanisms underlying nanowire formation in pristine phthalocyanine

Understanding the molecular processes of nanowire self-assembly is crucial for designing and controlling nanoscale structures that could lead to breakthroughs in functional materials. In this work, we focus on pristine phthalocyanines as a representative example of mesogenic supramolecular assemblies and have analyzed the formation of nanowires using classical molecular dynamics simulations. In the simulations, the molecules spontaneously form multi-columnar structures resembling supramolecular polymers that subsequently grow into more ordered aggregates. These self-assemblies are concentration dependent, leading to the formation of multi-columnar, dynamic aggregates at higher concentrations and nanowires at lower concentrations. The multi-columnar assemblies on a whole are more disordered than the nanowires, but have locally ordered domains of parallel facing molecules that can fluctuate while maintaining their overall shape. The nanowire formation at lower concentrations involves the initial interaction and clustering of randomly oriented phthalocyanine molecules, followed by the merging of small clusters into elongated segments and the eventual formation of a stable nanowire. We observe three main conformers in these self-assemblies, the parallel, T-shaped and edge-to-edge stacking of the phthalocyanine dimers. We calculate the underlying free energy landscape and show that the parallel conformers form the most stable configuration which is followed by the T-shaped and edge-to-edge dimer configurations. The findings provide insights into the mechanisms and pathways of nanowire formation and a step towards the understanding of self-assembly processes in supramolecular mesogens.

Columnar Aggregates of Azobenzene Stars: Exploring Intermolecular Interactions, Structure, and Stability in Atomistic Simulations

We present a simulation study of supramolecular aggregates formed by three-arm azobenzene (Azo) stars with a benzene-1,3,5-tricarboxamide (BTA) core in water. Previous experimental works by other research groups demonstrate that such Azo stars assemble into needle-like structures with light-responsive properties. Disregarding the response to light, we intend to characterize the equilibrium state of this system on the molecular scale. In particular, we aim to develop a thorough understanding of the binding mechanism between the molecules and analyze the structural properties of columnar stacks of Azo stars. Our study employs fully atomistic molecular dynamics (MD) simulations to model pre-assembled aggregates with various sizes and arrangements in water. In our detailed approach, we decompose the binding energies of the aggregates into the contributions due to the different types of non-covalent interactions and the contributions of the functional groups in the Azo stars. Initially, we investigate the origin and strength of the non-covalent interactions within a stacked dimer. Based on these findings, three arrangements of longer columnar stacks are prepared and equilibrated. We confirm that the binding energies of the stacks are mainly composed of π-π interactions between the conjugated parts of the molecules and hydrogen bonds formed between the stacked BTA cores. Our study quantifies the strength of these interactions and shows that the π-π interactions, especially between the Azo moieties, dominate the binding energies. We clarify that hydrogen bonds, which are predominant in BTA stacks, have only secondary energetic contributions in stacks of Azo stars but remain necessary stabilizers. Both types of interactions, π-π stacking and H-bonds, are required to maintain the columnar arrangement of the aggregates.

Tuning Supramolecular Polymer Assembly through Stereoelectronic Interactions

The supramolecular polymerization of 2,11-dithia[3.3]paracyclophanes through self-complementary intermolecular and transannular amide hydrogen bonding is presented. An n → π* interaction between the amide hydrogen bonding units and the central bridging atom results from the single-point exchange of a carbon atom for a sulfur atom. This orbital donor-acceptor interaction can be strengthened by oxidizing the sulfide to a sulfone which acts to shorten the donor···acceptor distance and increase orbital overlap. Experimental signatures of the increased n → π* interaction include larger isodesmic polymerization elongation constants in solution, changes in characteristic bond stretching frequencies, and geometric/structural changes evaluated by X-ray crystallography. The experimental data are supported by extensive computational investigations of both assembling and nonassembling 2,11-dithia[3.3]paracyclophanes as well as a rationally designed model system to confirm the role of stereoelectronic effects on supramolecular polymer assembly.

Dipolar relaxation in thin films of supramolecular stacks of benzenecarboxamides and insights to enhance their ferroelectric characteristics

The relationship between molecular structure and ferroelectric behaviour of thin films is explored in an all-organic supramolecular polymer material based on benzenecarboxamides, using atomistic molecular dynamics simulations. While increasing the number of amide groups around the phenyl core increases the dipole density of a molecule, increasing the length of the corresponding alkyl groups decreases the same. The interplay between these two contributions displays a rich behaviour on key material characteristics, in particular, the polarisation retention time. The latter is shown to be inversely proportional to the alkyl chain length, a consequence of weaker interactions between macrodipoles of stacks. Polarisation retention time was observed to be the highest in a molecule with five amide groups around the aromatic phenyl core which is explained as due to the large barrier for amide group rotation, which is one of the crucial channels for dipolar relaxation. Simulations also demonstrate that the barrier, however, does not affect the switchability of polarization, upon field reversal.

Stacks of Azobenzene Stars: Self-Assembly Scenario and Stabilising Forces Quantified in Computer Modelling

In this paper, the columnar supramolecular aggregates of photosensitive star-shaped azobenzenes with benzene-1,3,5-tricarboxamide core and azobenzene arms are analyzed theoretically by applying a combination of computer simulation techniques. Without a light stimulus, the azobenzene arms adopt the trans-state and build one-dimensional columns of stacked molecules during the first stage of the noncovalent association. These columnar aggregates represent the structural elements of more complex experimentally observed morphologies-fibers, spheres, gels, and others. Here, we determine the most favorable mutual orientations of the trans-stars in the stack in terms of (i) the π - π distance between the cores lengthwise the aggregate, (ii) the lateral displacements due to slippage and (iii) the rotation promoting the helical twist and chirality of the aggregate. To this end, we calculate the binding energy diagrams using density functional theory. The model predictions are further compared with available experimental data. The intermolecular forces responsible for the stability of the stacks in crystals are quantified using Hirshfeld surface analysis. Finally, to characterize the self-assembly mechanism of the stars in solution, we calculate the hydrogen bond lengths, the normalized dipole moments and the binding energies as functions of the columnar length. For this, molecular dynamics trajectories are analyzed. Finally, we conclude about the cooperative nature of the self-assembly of star-shaped azobenzenes with benzene-1,3,5-tricarboxamide core in aqueous solution.

Do Columns of Azobenzene Stars Disassemble under Light Illumination?

The clustering properties of star-shaped molecules comprising three photochromic azobenzene-containing arms are investigated with specific focus on the influence of light on these structures. Previous experimental works report self-assembly of azobenzene stars in aqueous solution into long columnar clusters that are detectable using optical microscopy. These clusters appear to vanish under UV irradiation, which is known to induce trans-to-cis photoisomerization of the azobenzene groups. We have performed MD simulations, density functional theory, and density functional tight binding calculations to determine conformational properties and binding energies of these clusters. Our simulation data suggest that the binding strength of the clusters is large enough to prevent a breaking along their main axis. We conclude that very likely other mechanisms lead to the apparent disappearance of the clusters.