Supramolecular Polymerization of N,N′,N″,N‴-tetra-(Tetradecyl)-1,3,6,8-pyrenetetracarboxamide: A Computational Study

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.

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.

Differentiating the mechanism of self-assembly in supramolecular polymers through computation

The mechanism by which monomers in solution, beyond a certain concentration or below a certain temperature, self-assemble to form one dimensional supramolecular polymers determines much of the bulk properties of the polymer.

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.

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.

Hierarchical computational modelling approaches for the study of supramolecular polymers is reviewed.