Ferro- and ferrielectricity and negative piezoelectricity in thioamide-based supramolecular organic discotics

Amide-based discotic supramolecular organic materials are of interest for fundamental understanding of cooperative self-assembly and collective dipole switching mechanisms as well as for practically relevant ferroelectric and piezoelectric properties. Here, we show how replacing amides (dipole moment of ∼3.5 D) with thioamides (∼5.1 D) as dipolar moieties in the archetypal C₃-symmetric discotic molecule BTA leads to ferroelectric materials with a higher remnant polarization and lower coercive field. The thioamide-based materials also demonstrate a rare negative piezoelectricity and a previously predicted, yet never experimentally observed, polarization reversal via asymmetric intermediate states, that is, ferrielectric switching.

Cyclic Photoisomerization of Azobenzene in Atomistic Simulations: Modeling the Effect of Light on Columnar Aggregates of Azo Stars

This computational study investigates the influence of light on supramolecular aggregates of three-arm azobenzene stars. Every star contains three azobenzene (azo) moieties, each able to undergo reversible photoisomerization. In solution, the azo stars build column-shaped supramolecular aggregates. Previous experimental works report severe morphological changes of these aggregates under UV-Vis light. However, the underlying molecular mechanisms are still debated. Here we aim to elucidate how light affects the structure and stability of the columnar stacks on the molecular scale. The system is investigated using fully atomistic molecular dynamics (MD) simulations. To implement the effects of light, we first developed a stochastic model of the cyclic photoisomerization of azobenzene. This model reproduces the collective photoisomerization kinetics of the azo stars in good agreement with theory and previous experiments. We then apply light of various intensities and wavelengths on an equilibrated columnar stack of azo stars in water. The simulations indicate that the aggregate does not break into separate fragments upon light irradiation. Instead, the stack develops defects in the form of molecular shifts and reorientations and, as a result, it eventually loses its columnar shape. The mechanism and driving forces behind this order-disorder structural transition are clarified based on the simulations. In the end, we provide a new interpretation of the experimentally observed morphological changes.

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.

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.

Kinetic Monte Carlo simulations of organic ferroelectrics

Ferroelectrics find broad applications, e.g. in non-volatile memories, but the switching kinetics in real, disordered, materials is still incompletely understood. Here, we develop an electrostatic model to study ferroelectric switching using 3D Monte Carlo simulations. We apply this model to the prototypical small molecular ferroelectric trialkylbenzene-1,3,5-tricarboxamide (BTA) and find good agreement between the Monte Carlo simulations, experiments, and molecular dynamics studies. Since the model lacks any explicit steric effects, we conclude that these are of minor importance. While the material is shown to have a frustrated antiferroelectric ground state, it behaves as a normal ferroelectric under practical conditions due to the large energy barrier for switching that prevents the material from reaching its ground state after poling. We find that field-driven polarization reversal and spontaneous depolarization have orders of magnitude different switching kinetics. For the former, which determines the coercive field and is relevant for data writing, nucleation occurs at the electrodes, whereas for the latter, which governs data retention, nucleation occurs at disorder-induced defects. As a result, by reducing the disorder in the system, the polarization retention time can be increased dramatically while the coercive field remains unchanged.

Resistive switching in an organic supramolecular semiconducting ferroelectric

The combination of switchable dipolar side groups and the semiconducting core of the newly synthetized C3-symmetric benzotrithiophene molecule (BTTTA) leads to an ordered columnar material showing continuous tunability from injection- to bulk-limited conductivity modulation.

Molecular modelling of supramolecular one dimensional polymers

Supramolecular polymers exemplify the need to employ several computational techniques to study processes and phenomena occuring at varied length and time scales. Electronic processes, conformational and configurational excitations of small aggregates of chromophoric molecules, solvent effects under realistic thermodynamic conditions and mesoscale morphologies are some of the challenges which demand hierarchical modelling approaches. This review focusses on one-dimensional supramolecular polymers, the mechanism of self-assembly of monomers in polar and non-polar solvents and properties they exhibit. Directions for future work are as well outlined.

How the Dynamics of a Supramolecular Polymer Determines Its Dynamic Adaptivity and Stimuli-Responsiveness: Structure-Dynamics-Property Relationships From Coarse-Grained Simulations

The rational design of supramolecular polymers that can adapt or respond in time to specific stimuli in a controlled way is interesting for many applications, but this requires understanding the molecular factors that make the material faster or slower in responding to the stimulus. To this end, it is necessary to study the dynamic adaptive properties at submolecular resolution, which is difficult at an experimental level. Here we show coarse-grained molecular dynamics simulations (<5 Å resolution) demonstrating how the dynamic adaptivity and stimuli responsiveness of a supramolecular polymer is controlled by the intrinsic dynamics of the assembly, which is in turn determined by the structure of the monomers. As a representative case, we focus on a water-soluble 1,3,5-benzenetricarboxamide (BTA) supramolecular polymer incorporating (charged) receptor monomers, experimentally seen to undergo dynamic clustering following the superselective binding to a multivalent recruiter. Our simulations show that the dynamic reorganization of the supramolecular structure proceeds via monomer diffusion on the dynamic fiber surface (exchange within the fiber). Rationally changing the structure of the monomers to make the fiber surface more or less dynamic allows tuning the rate of response to the stimulus and of supramolecular reconfiguration. Simple in silico experiments draw a structure-dynamics-property relationship revealing the key factors underpinning the dynamic adaptivity and stimuli-responsiveness of these supramolecular polymers. We come out with clear evidence that to master the bioinspired properties of these fibers, it is necessary to control their intrinsic dynamics, while the high-resolution of our molecular models permits us to show how.

Communication: Self-assembly of a model supramolecular polymer studied by replica exchange with solute tempering

Conventional molecular-dynamics (cMD) simulation has a well-known limitation in accessible time and length scales, and thus various enhanced sampling techniques have been proposed to alleviate the problem. In this paper, we explore the utility of replica exchange with solute tempering (REST) (i.e., a variant of Hamiltonian replica exchange methods) to simulate the self-assembly of a supramolecular polymer in explicit solvent and compare the performance with temperature-based replica exchange MD (T-REMD) as well as cMD. As a test system, we consider a relatively simple all-atom model of supramolecular polymerization (namely, benzene-1,3,5-tricarboxamides in methylcyclohexane solvent). Our results show that both REST and T-REMD are able to predict highly ordered polymer structures with helical H-bonding patterns, in contrast to cMD which completely fails to obtain such a structure for the present model. At the same time, we have also experienced some technical challenge (i.e., aggregation-dispersion transition and the resulting bottleneck for replica traversal), which is illustrated numerically. Since the computational cost of REST scales more moderately than T-REMD, we expect that REST will be useful for studying the self-assembly of larger systems in solution with enhanced rearrangement of monomers.

Unveiling the Role of Macrodipolar Interactions in the Properties of Self-Assembled Supramolecular Materials

Self-assembling of supramolecules composed of benzene and cyclohexane tricarboxamide derivatives can form highly organized 1 D fibers exhibiting macrodipoles. The way fibers pack in the condensed phase governs the final properties of the supramolecular material, in which macrodipoles can be oriented parallel or antiparallel to each other, and their magnitude can be tuned by additional intra-columnar dipole stabilization. X-ray structural elucidation of these materials remains a real challenge due to the difficulty in growing single crystals. This problem can be tackled by using atomistic molecular dynamics to simulate supramolecular materials composed of cyclohexanetricarboxamide derivatives assuming different magnitudes and orientations of macrodipoles in the condensed phase, as we show here. The results provide insight on the isotropization mechanism of the supramolecules and also reveal that the relative orientation between macrodipoles can indeed influence their stability. This work nicely complements X-ray structural characterizations of supramolecular materials, and helps understand structure-property relationships of a range of similar noncovalent materials.

Photoisomers of Azobenzene Star with a Flat Core: Theoretical Insights into Multiple States from DFT and MD Perspective

This study focuses on comparing physical properties of photoisomers of an azobenzene star with benzene-1,3,5-tricarboxamide core. Three azobenzene arms of the molecule undergo a reversible trans-cis isomerization upon UV-vis light illumination giving rise to multiple states from the planar all-trans one, via two mixed states to the kinked all-cis isomer. Employing density functional theory, we characterize the structural and photophysical properties of each state indicating a role the planar core plays in the coupling between azobenzene chromophores. To characterize the light-triggered switching of solvophilicity/solvophobicity of the star, the difference in solvation free energy is calculated for the transfer of an azobenzene star from its gas phase to implicit or explicit solvents. For the latter case, classical all-atom molecular dynamics simulations of aqueous solutions of azobenzene star are performed employing the polymer consistent force field to shed light on the thermodynamics of explicit hydration as a function of the isomerization state and on the structuring of water around the star. From the analysis of two contributions to the free energy of hydration, the nonpolar van der Waals and the electrostatic terms, it is concluded that isomerization specificity largely determines the polarity of the molecule and the solute-solvent electrostatic interactions. This convertible hydrophilicity/hydrophobicity together with readjustable occupied volume and the surface area accessible to water, affects the self-assembly/disassembly of the azobenzene star with a flat core triggered by light.

Unravelling the Pathway Complexity in Conformationally Flexible N-Centered Triarylamine Trisamides

Two families of C3 -symmetrical triarylamine-trisamides comprising a triphenylamine- or a tri(pyrid-2-yl)amine core are presented. Both families self-assemble in apolar solvents via cooperative hydrogen-bonding interactions into helical supramolecular polymers as evidenced by a combination of spectroscopic measurements, and corroborated by DFT calculations. The introduction of a stereocenter in the side chains biases the helical sense of the supramolecular polymers formed. Compared to other C3 -symmetrical compounds, a much richer self-assembly landscape is observed. Temperature-dependent spectroscopy measurements highlight the presence of two self-assembled states of opposite handedness. One state is formed at high temperature from a molecularly dissolved solution via a nucleation-elongation mechanism. The second state is formed below room temperature through a sharp transition from the first assembled state. The change in helicity is proposed to be related to a conformational switch of the triarylamine core due to an equilibrium between a 3:0 and a 2:1 conformation. Thus, within a limited temperature window, a small conformational twist results in an assembled state of opposite helicity.

From Cooperative Self-Assembly to Water-Soluble Supramolecular Polymers Using Coarse-Grained Simulations

Supramolecular polymers, formed via noncovalent self-assembly of elementary monomers, are extremely interesting for their dynamic bioinspired properties. In order to understand their behavior, it is necessary to access their dynamics while maintaining high resolution in the treatment of the monomer structure and monomer-monomer interactions, which is typically a difficult task, especially in aqueous solution. Focusing on 1,3,5-benzenetricarboxamide (BTA) water-soluble supramolecular polymers, we have developed a transferable coarse-grained model that allows studying BTA supramolecular polymerization in water, while preserving remarkable consistency with the atomistic models in the description of the key interactions between the monomers (hydrophobic, H-bonding, etc.), self-assembly cooperativity, and amplification of order into the growing fibers. This permitted us to monitor the amplification of the key interactions between the monomers (including H-bonding) in the BTA fibers during the dynamic polymerization process. Our molecular dynamics simulations provide a picture of a stepwise cooperative polymerization mechanism, where initial fast hydrophobic aggregation of the BTA monomers in water is followed by the slower reorganization of these disordered aggregates into ordered directional oligomers. Supramolecular polymer growth then proceeds on a slower time scale. We challenged our models via comparison with the experimental evidence, capturing the effect of temperature variations and subtle changes in the monomer structure on the polymerization and on the properties of the fibers seen in the real systems. This work provides a multiscale spatiotemporal characterization of BTA self-assembly in water and a useful platform to study a variety of BTA-based supramolecular polymers toward structure-property relationships.

Dynamic self-assembled polymer: HCl responsive inversion of supramolecular polymer handedness

Discotic trisamide formed a self-assembled polymer and exhibits inversion of supramolecular polymer handedness in the presence of HCl.

Influence of Metal, Ligand and Solvent on Supramolecular Polymerizations with Transition-Metal Compounds: A Theoretical Study

The nature of intermolecular interactions governing supramolecular polymerizations is very important for controlling their cooperativity. In order to address this problem, supramolecular columns made of Pt^II and Pd^II complexes of oligo(phenylene ethynylene)-based pyridine (OPE) and tetrazolylpyridine ligands (TEP) were investigated through the dispersion-corrected PM6 method. Aromatic, CH-π, M-Cl and metallophilic interactions helped stabilize the supramolecules studied, and their geometries and associated cooperativities were in excellent agreement with experimental data. The OPE ligand and/or the presence of Pt^II led to stronger metallophilic interactions and also to cooperative supramolecular polymerizations, which clearly suggests that metallophilic interactions are a key factor for controlling cooperativity. The results indicate that sequential monomer addition is in general less spontaneous than the combination of two larger preformed stacks. The present theoretical investigations contribute to the further understanding of the relation between the thermodynamics of supramolecular polymerizations and the nature of different synthons.