Next Article in Journal
Design, Synthesis, and Structural Characterization of Novel Diazaphenothiazines with 1,2,3-Triazole Substituents as Promising Antiproliferative Agents
Next Article in Special Issue
The Self-Aggregation of Porphyrins with Multiple Chiral Centers in Organic/Aqueous Media: The Case of Sugar- and Steroid-Porphyrin Conjugates
Previous Article in Journal
Development of a Microfluidic Array to Study Drug Response in Breast Cancer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Vladyslav Savchenko
1,2,†,
Markus Koch
2,†,
Aleksander S. Pavlov
3,
Marina Saphiannikova
1,2 and
Olga Guskova
1,2,*
1
Dresden Center for Computational Materials Science (DCMS), Technische Universität Dresden, 01062 Dresden, Germany
2
Institute Theory of Polymers, Leibniz Institute of Polymer Research Dresden, Hohe Str. 6, 01069 Dresden, Germany
3
Department of Physical Chemistry, Faculty of Chemistry and Technology, Tver State University, Sadovyj per. 35, Tver 170002, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Molecules 2019, 24(23), 4387; https://doi.org/10.3390/molecules24234387
Submission received: 12 November 2019 / Revised: 25 November 2019 / Accepted: 28 November 2019 / Published: 30 November 2019
(This article belongs to the Special Issue Self-Aggregation in Supramolecular Systems)

Abstract

:
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.

1. Introduction

In recent years, azobenzene-containing molecules (azos), which are reversibly responding to the light stimulus, have been widely applied in science and engineering to construct materials with various photosensitive characteristics [1,2,3,4,5]. Among others, star-shaped azos, i.e., multi-photochromic compounds with azobenzene moieties attached to the central unit as arms, have received considerable attention in the research of the past decade. On one hand, such molecules combine mesomorphism and photoswitching [6,7,8], they can more effectively reply to the light stimulus (for example, during inscribing a surface relief grating) as compared to monomeric azos. On the other hand, they can be readily synthesized and purified, be either evaporated or cast from solution, build stable amorphous phases and exhibit more uniform physical properties. There is no dilution of the chromophores or any effects coming from entanglements unlike azo–photochromic polymers and composite polymer systems [9,10]. Star-shaped azos are also different from linear or cyclic azo n-mers with the same number of azo-fragments (trimers, tetramers, etc.). First of all, this is manifested in different supramolecular structures the stars and the linear or cyclic [11] n-mers can build. Moreover, the branching results in lower melting temperatures, and lower viscosity [6,12,13], meaning that these properties are topology-dependent.
There are several classifications of star-shaped azos: they may differ in (i) the number of arms (three [9,10,14,15,16,17,18,19,20], four [7,21,22], six [23,24,25], nine [26]), (ii) in their centres (nitrogen [9,19], phosphorus [27], silicon [28] or carbon atoms [21,22], benzene [14,16,20,29,30,31,32], polycyclic aromatic hydrocarbons [25,33], heterocyclic hydrocarbons [34,35,36,37], some chiral groups [6] or bio-active residues [38]), (iii) in the way azobenzenes are attached to the centre (covalently or noncovalently [39]), (iv) in their conformational rigidity (flexible or rigid), and (v) in their overall geometry (quasi-planar or 3D). Both the conformational rigidity and the shape of the star are closely related to the chemical nature of the star core. For instance, the planar conjugated fragment in the absence of any fatty linkers between the centre and the azobenzene arms supports the planarity of the star and its rigidity, and vice versa, long hydrocarbon linkers provoke its conformational flexibility [8,30,31].
In this work, we investigate aqueous solutions of azobenzene star molecules with a planar benzene-1,3,5-tricarboxamide (BTA) core and three azo arms (azo–BTA star, see Figure 1a). This star is rigid and quasi planar in its all-trans state, as recently confirmed by us [14,29]. Its BTA core is a well-known supramolecular entity in colloid science, used for the fabrication of dynamic and responsive materials [40,41,42,43,44,45,46]. The interactions between BTA molecules are inherently dynamic and involve hydrogen bonding between three amide groups, π π stacking interactions between benzene cores, and hydrophobic forces originating from the arms attached to the centre. Tuning these interactions and applying a light stimulus, one can control the colloidal assembly of such molecules in solution.
As literature survey shows, azo–BTA stars that differ from the one shown in Figure 1a only in the para-substituents at the star periphery, demonstrate light-controlled self-assembly and actuation [17,26,47]. For example, Choi et al. observed the construction of nanocolumns and further lateral self-assembly into the low-ordered hexagonal columnar liquid crystal phase or lamellar-columnar crystals [26]. Not only the shape of supramolecular assemblies but also their color are shown to be remote-controlled by light [26]. Lee et al. have studied the solvent dependent generation of microfibers, gels, and hollow spheres in aqueous organic solvent systems and their photoinduced and reversible phase transitions upon UV-Vis light illumination [17]. One new study demonstrated the perpetual columnar self-assembly of azo–BTA discotic mesogenic blocks at room temperature which hold together despite light stimulation [48]. Another interesting example of self-assembly of nonchiral azo–BTA stars into mesophases with photo-switchable chirality was demonstrated at a recent time by Choi et al. [49].
Light-sensitivity of azo–BTA stars has been used for the fabrication of lyotropic liquid crystalline alignment media [50], for the light-controlled synthesis of gold nanoparticles [51], for the inscription of surface relief [10,52,53] and holographic [32] gratings, as well as for patterning and molecular alignment [30,31,54]. In a recently published review [55], it has been highlighted that crystalline phases consisting of azo–BTA stars show remarkable promise to be applied as photomechanical crystalline nanowire actuator utilized for the gripper of a micropipette [18]. Moreover, azo–BTA stars, owing to their chemical structure and depending on the experimental conditions, show the ability to form homogeneous molecular glasses [10,32,52,53] (an amorphous phase), to build a liquid-crystalline phase [26,30,31,54], and to self-assemble into ordered crystalline arrays [17,18].
As pointed out in numerous publications devoted to the self-assembly of azo–BTA stars, the hydrogen bonding and π π stacking between the neighboring star centers along the columns play an essential role in this system. These two interactions are responsible for assembling a single structure from very first aggregates, that appear at the initial stage of the ordering process. The question of quantifying these stabilizing forces is open and is, therefore, one of the goals of this study.
Our attention will also be given to the self-assembly mechanism of azo–BTA stars. According to the latest scientific ideas in this field, three scenarios of the supramolecular assembly are distinguished: cooperative, isodesmic and ring-chain mechanisms [56]. Two of them, cooperative [43,44] and isodesmic [40,41,42], were identified for BTA-centered molecules. In a recent review, the “structure–mechanism” correlation has been established and the ordering behaviors of BTA-containing derivatives were placed on the so-called cooperativity scale [45]. Thus, the analysis of the signatures of self-assembly is the second aim of this study.
The remainder of this article is divided into three sections. The next chapter contains the discussion of the results and consists of three sub-chapters: The binding energy diagrams calculated for the azo–BTA dimers are discussed first, followed by the Hirshfeld surface analysis of the crystalline azo–BTA star phase found in the literature [17], and then concluded by the quantification of the signatures of cooperative growth of the star-shape azobenzenes stacks. This section is followed by conclusions and how the results can be interpreted in the perspective of previous studies and of the working hypotheses. Lastly, we give an overview of the models and computer simulation techniques applied in the study.

2. Results and Discussion

2.1. Binding Energy Diagrams

We begin the description of the results with a discussion of the dimers made of azo–BTA stars (Figure 1a). The orientation of the molecules in dimers is parallel, and hence their structure can be quantified by their segregation meaning that they are placed one above another. The binding energy (EB in kJ mol−1) for the molecular dimers and the assemblies of larger N (see below) is estimated as
E B = E N N E 1 N 1 ,
where N is the assembly size, EN is the energy of the column of size N, and E1 is the energy of the monomer, i.e., the azo–BTA star. More details can be found in Section 4.
The varied parameter for calculation of the binding strength is the distance d between the centers of the BTA cores, as schematically shown in Figure 2. The selected d values are ranging from 3 to 6 A ˚ with a step of 0.1 A ˚ which is near experimentally observed π π distances. The task here is to determine the most energetically favorable distance between the star cores in the column. This task sounds rather trivial, but on one side, the knowledge of this distance is necessary for subsequent simulations, and on another side, the experimental benchmarking using data from the literature allows us to learn about the accuracy of the computer-based research methods.
The analysis of the experimental data [17,48,57,58,59,60] on the ordered phases of various BTA–cored three-armed stars reveals that the intermolecular separations in the stacks range from 3.51 to 3.89 A ˚ for the molecules with different periphery around the star center. Apparently, there is no clear relationship between the distance d and the bulkiness of the star arms. For example, for almost the same azobenzene [17] and alkoxy azobenzene [48] BTA–stars, the values differ by 0.47 A ˚ , whereas the π π stacking separations are nearly equal to 3.6 A ˚ for shorter (2-methoxyethyl [57] and carboxymethyl [60]) and longer/bulkier arms (linear and branched C8–C11 alkyls [58], Glu-based oligopeptides [59]). The reason may lie in the many possibilities the molecules can be assembled in the ordered phases: The stars can build segregated columns, be rotated around the core in the star plane or shifted laterally along the columns.
We would like to stress here that the values found in experimental investigations of similar compounds belong to the range of negative binding strengths. The strongest coupling between the azo–BTA stars in computer simulations corresponds to d = 3.6 A ˚ , which is in the range of the experimental values discussed above [17,48]. The energetics of the binding is ca. −200 kJ mol−1. For example, the binding of methyl-substituted BTA proceeds with the energy of ca. −150 kJ mol−1, as evaluated in DFT calculations with counterpoise correction by Filot et al. [44]. Taking into account additional dispersion interactions coming from azobenzene arms, the estimated energy for the azo–BTA star seems reasonable. Increasing d up to 3.9 A ˚ has only little effect on the binding of the stars in the column, indicating that the separations 3.5 ≤ d ≤ 3.9 A ˚ could be realized without significant energy losses (the variations in EB are 1–20 kJ mol−1).
As mentioned above, not only the segregated stacks were found in experiments or predicted theoretically, but also molecules rotated around the star center [59,60] leading to helically twisted arrangements [42,44,49,58,61] in the columns or shifted laterally along the stack [17]. Therefore, as the next step of the study, diagrams are constructed showing the binding energy EB depending on the lateral shifts l, k and the twist angle Ω in star dimer at d = 3.6 A ˚ (Figure 3). The l-shift assumes the displacement of one star in dimer along one of the arms, whereas the k-shift occurs in between two arms perpendicular to the l-shift. The angle Ω , defined as the rotation angle of the upper star, is shown schematically in Figure 3. Such a systematic rotation of the consecutive stars along the column leads to the formation of the helical structures found in the literature for BTA-centered molecules [44,58,62,63,64]. It was discussed that this twisting allows the binding of the stars through hydrogen bonding, in which all three amide groups of the BTA are involved. Moreover, it reduces possible steric clashes of the arms.
The diagrams provide a graphical visualization of the effects of the variables l, k and Ω on the intermolecular coupling, i.e. the binding strength. Both panels of Figure 3 illustrate that the most stable stacking is achieved if the stars are perfectly segregated, or slightly shifted/rotated. Concerning the rotation of the stars in the stack (Figure 3a), which is favorable for the triple helical hydrogen bonding, we can conclude that rotation is strongly suppressed in the perfectly stacked columns at the Ω values of 5–10°, where the binding proceeds with positive energy. The coupling becomes more developed for the stacks with shifted molecules along one of the arms and for rotation angles Ω > 10°, meaning that the energy is always negative. For example, for small l-shifts and Ω = 0° there is a relatively large area with energies ranging from −50 to −30 kJ mol−1 (marked with red color), and another big area for angles 10 < Ω 60° with energies from −120 to −70 kJ mol−1 for l-shifts from four to eight A ˚ (marked with green color).
Let us now compare the literature data with our predictions. The most common rotational angle between consecutive unmoved (l = k = 0) BTA-centred monomers is 60–65° [44,58,62,63,64]. The calculations show that the binding energy corresponding to these coordinates is about −73 kJ mol−1 for the azo–BTA star (marked with yellow color). Interestingly, around Ω = 60° a local minimum on this symmetrical diagram is observed. Smaller rotational angles are also found in literature, for instance, 13–16° for the consecutive stars [40].
The second panel (Figure 3b) shows a wide range of binding energies as well. A combination of small k and l up to 3–4 A ˚ is characterized by the strongest binding. The experimentally found coordinates [17] belong to this global minimum, which coincides with our predictions. The most disadvantageous mutual orientation of stars in the stack is at k = 8–10 A ˚ and l = 0–2 A ˚ , which is caused by the steric repulsion of the molecular centre and dimethylamino groups at the arm ends. In both panels, the dramatic shifts weaken the coupling considerably.
Note here, that the binding diagrams are calculated with the assumption that the molecular conformation (in our case, almost planar geometry of the star after its optimization [14]) does not alter upon dimerization, translation or rotation. On the contrary, the formation of hydrogen bonds requires not only the rotation of the stars but also the alteration of the planarity of the BTA cores. More realistic binding models will be considered further when we will discuss the mechanisms of self-assembly of azo–BTA stars in aqueous solution.

2.2. Analysis of the Trans-Azo–BTA Star Crystal

In this section, we move from the simplified dimer representation to a more complex organization of the azo–stars, as experimentally investigated by Lee et al. [17]. The crystalline phase found in literature [17] is shown in Figure 1d.
We are most interested in studying the intermolecular contacts, holding together or destabilizing the stacks of stars in the crystal. To probe them, we use the Hirshfeld surface analysis [65] (Figure 4, left column) and plot the de-di diagrams (Figure 4, right column). The color scheme of the surface codes the sign of the normalized contact distance dnorm defined in terms of de, di and the van der Waals (vdW) radii of the atoms (see Section 4 for more detail).
Here, d e and d i are the distances from the Hirshfeld surface to the nearest atom outside (external) and inside (internal) the surface, correspondingly. The color of the Hirshfeld surface changes from red (showing the distances shorter than the sum of vdW radii) through white (distances are equal to the sum of vdW radii) to blue (distances are longer than the sum of vdW radii). The fraction of the corresponding atomic contacts is also given in percent.
For example, the contacts between oxygen and hydrogen atoms (Figure 4a) are shorter than the sum of vdW radii, indicating the hydrogen-bonded arms along the molecular stack. The distance of the separation between two participating atoms is 2.24 A ˚ . Two sharp lateral spikes with a minimum de + di ≈ 2.24 A ˚ correspond to the O⋯H contacts (Figure 4a). The O⋯H contacts comprise 5.0% of the surface.
The Hirshfeld surface provides evidence for C–H π interactions (between hydrogen atoms of the phenyl ring of an azobenzene arm and the carbon atoms of the ring from the neighboring star) beyond dispersion, which further links the stars into columns. These contacts have the largest fraction of 13.3% from the intermolecular interactions in the crystal. The separation of these contacts is 2.57 A ˚ (Figure 4c).
A very small fraction of close contacts making up only <1% of the surface occurs between two hydrogen atoms of the phenyl ring at the star periphery with distances of 2.04 or 2.30 A ˚ . This type of close contact is on the spikes of the fingerprint plots as well (Figure 4d). The N⋯H, C⋯C, N⋯N contacts are marked as blue, meaning that these atoms are relatively far from each other in the crystals.
There are two important conclusions from the Hirshfeld surface analysis. First of all, the visual inspection of these charts shows the molecular asymmetry in terms of interactions with the surrounding molecules in the crystal. Secondly, the strongest interactions are realized along the stack due to the hydrogen bonds and C–H π contacts. The last stabilizing interaction contributes to the lateral coupling between the columns.
The considered intermolecular interactions only qualitatively describe the forces holding together the ordering structures. Our next task is the quantification of the interactions. Partially, this has been solved recently by us [29] for the same object of the study. In that article, we have shown that the binding along the column is much stronger (up to a factor of six) than the intercolumnar coupling. The next section will be focused on the columnar stacks (Figure 1b,e) in order to establish which self-assembly mechanism the azo–BTA stars obey.

2.3. On the Cooperative Nature of Self-Assembly

The noncovalent association into one-dimensional columns is the first step in the appearance of complex supramolecular structures. Discovering the factors that affect the growth from a single azo–BTA star to the stacks and columns is a step forward to understanding the multi-structured morphologies found experimentally. One of these factors is the assembly mechanism, which is determined by the dependence of the association constant for monomer addition (in our case, azo–BTA star) on the assembly size N [56]. Here, the growth is isodesmic if the binding of each new monomer to the growing chain is independent of N, and cooperative if the binding shows enhanced propensity for monomer addition as N increases [44].
Our working hypothesis is that azo–BTA stars gather together cooperatively and that this growth in binding strength is due to electronic effects. According to literature data, these effects are manifested primarily based on pairwise electrostatic interactions, including long-range dipole-dipole interaction, non-pairwise short-range polarization, and resonance-assisted hydrogen bonding, and caused by the redistribution of the electron density along the column. As descriptors of cooperativity, the following features are analyzed as functions of N or depending on the position of the star in the stack: binding energies, the length of hydrogen bonds, and macrodipoles. These descriptors have been used by many researchers for the classification of the self-assembly mechanisms in BTA-containing systems [44,61,64,67]. The stacks of azo–BTA stars are analyzed using data sets from molecular dynamics (MD) simulations [29] and from experiments [17] applying density functional theory and classical modelling techniques.
In Figure 5, we compare the binding energies per pair for azo–star columns of various length from experiments, i.e., for the stack shown in Figure 1e [17] (panel a), and modelled in MD simulations (panel b). Both data sets show an obvious N-dependence. This behaviour is analogous to results reported earlier [44,61,64,67]. The interaction strength increases with the number of molecules, proving that the formation of longer columns is energetically increasingly favorable, or in other words cooperative, as compared to shorter associates. On both panels, the binding energy decreases toward an asymptotic value for the infinitely long supramolecule. The values for the stacks with N = 20 are comparable to each other with −213 and −296 kJ mol−1 for Lee’s column (Figure 1e) and stacks from our recent publication [29], respectively. Here, it is important to underline that for such long columns of molecules with an all-atom representation of the star, the self-assembly rules are analyzed for the first time. In previous works, the objects of the simulations were either very short columns, or coarse-grainedd modelled stacks.
To rationalize the origin of this cooperative effect inherent to the longer azo–BTA stacks, the cumulative parameter Σ is calculated (Figure 6). First, the data from this plot indicate that the local interactions become stronger as the stack grows: The values for longer stacks are systematically shifted to more negative energies. This is a result of the binding to a larger number of molecules that are located beyond the nearest neighbors. Correspondingly, the circumstance that the outer stars only have binding partners to one side, leads to a weakening of their attachment. Hence, we observe approximately halved energy values when comparing the outer stars to centrally located molecules within the same cluster. Secondly, the coupling closer to the middle of the column is stronger than the interactions near the edges of the stack. Consider, however, that the irregular shape of curves in Figure 6 could hint at weak points and irregularities in the cluster themselves caused by effects coming from solvent molecules and thermal motion simulated explicitly in MD runs. The mentioned irregularities are locations where a cluster may be destabilized in ensuing photoisomerization events when the clusters are exposed to light.
The next feature pointing out the cooperativity of the self-assembly is the shortening of the hydrogen bonds upon stack growth (Figure 7a). Such a behavior has been observed in computational studies of smaller BTA-cored molecules [44,61]. In contrast to this, the close contacts do not show any clear dependence on the assembly size, and they cannot be considered as one of the descriptors of cooperativity (Figure 7b). As follows from Table 1, the strength of hydrogen bonding also changes depending on the position inside the stack. Although the values do not give a monotonic dependence, the trend is clearly seen: The shortest hydrogen bonds are found close to the center of the stack. Similar conclusions have been drawn in a computational study by Albuquerque et al. [61].
To examine the long-range dipole–dipole interactions, the average dipole moment μ N and normalized macrodipole along the stacking direction μ N / μ max are calculated (Figure 8). The latter one is defined as the ratio of the electric dipole moment of an N-mer to the largest value of a dipole moment observed in the calculations, as described by Korlepara et al. [64]. In our calculations, μ max corresponds to N = 20. The dipole of the stack progressively grows (Figure 8a) due to the formation of hydrogen bonds, which requires the distortion of the planarity in the BTA-core: Weak intermolecular contacts force the rotation of the amide groups so that they do not remain in the plane of the benzene ring anymore. Most vividly it is seen for the stack from the crystalline structure [17]: The molecules have identical geometries, the dihedral angle Θ between the benzene core and amide group is around 30° for HB-bonded arms and ca. 8° otherwise, and the polar C = O groups point in the same direction increasing the long-axis component of the dipole.
A growing dipole of the stacks simulated in MD is predicted using different methods. However, the scattering of the data is substantial, which is related to a large number of molecular geometries that are realized at room temperature in computer simulations. For example, the mean value of the dihedral angle is Θ = 32.1 ± 17.8° for dimers. The carbonyl groups point in the same direction most of the time but not always, since thermal motion affects the molecular conformations. The values of the macrodipoles (Figure 8b) grow as well, which is in line with findings of other theoretical studies [61,67]. Thus, the amplification of electrostatic interactions along the stacks is the main factor leading to stronger hydrogen bonding between azo–BTA stars governing the cooperative growth of ordered columns.

3. Conclusions

Employing a number of computer simulation techniques, we have analyzed the intermolecular interactions between the azo–BTA stars in columns in experimentally resolved crystals [17] and as obtained in our recent molecular dynamics study [29]. Our results suggest that during the formation of ordered phases, these molecules prefer to build 1D columns with the equilibrium distance between the star cores of 3.6 A ˚ , in which stars are held together via weak interactions. These noncovalent interactions include among others the hydrogen bonding, whose length shortens with increasing the column length.
The binding energy, which accounts for all the noncovalent through-space interactions, grows with increasing size of the aggregate, as our calculations demonstrate. This phenomenon called cooperative self-assembly is proven for azo–BTA stars in aqueous phases. The simulations clearly attribute the cooperativity of the self-assembly to the development of a macrodipole, in other words, to the redistribution of the electronic density. This, in turn, enhances the strength of the hydrogen bonds, and subsequently the intermolecular coupling. We find that the predictions of binding energies and the molecular coordinates in the stacks coincide with data for the photosensitive azobenzene stars investigated in seminal works by Lee et al. [17,18,47].
The next steps of this work could include, for instance, the calculation of the cooperativity factor σ [45] in order to place the modelled azobenzene star on the cooperativity scale for BTA-cored derivatives. Another interesting topic is related to the recently published paper by Bochicchio et al. [68] on the light-induced evolution of an azobenzene-containing supramolecular tubule in view of our data showed in Figure 6, where the defects are characterized quantitatively for the azo–BTA column. Finally, here we considered the stacks of azobenzenes in the dark, i.e., the stars in all-trans state. It would be interesting to elucidate the behavior of all-cis azo–BTA stars and the energetics of their self-assembly, since, as we showed recently, the stacks made by all-cis azobenzene stars maintain their clusters upon light stimulus, even though their shape is more disordered, as compared to the stacks in the dark.

4. Materials and Methods

4.1. Objects of the Study and Their Models

In this paper, we consider three objects. The first one is a dimer constructed by doubling the all-trans azo–BTA star shown in Figure 1a in its optimized geometry, as obtained in our previous study [14]. The second one is the column of N = 2–20 stars formed in MD simulations in a recent paper [29], where N is the assembly size. One example with N = 20 is depicted in Figure 1b. Finally, the third object is the crystal of azo–BTA star experimentally characterized by Lee et al. [17] (Figure 1c–e). Here, we study the crystal itself, as well as the molecular stacks obtained by multiplying the original unit cell in b-direction 2–20 times. In this case, we get the columnar assemblies of the same size N as in MD simulations, but having identical molecular geometries, as revealed in experiments in ref. [17]. In all the cases, we use a full-atomistic representation of the stars, whereas the aqueous environment is modelled applying two schemes: In density functional theory (DFT) and density functional based tight-binding (DFTB+) calculations, the effects coming from the solvent are taken into account via the polarizable continuum model [69], whereas for the simulation of assemblies in molecular dynamics, water is simulated explicitly.
The columns, which vary in their size between two and 20 molecules, are manually assembled for the subsequent MD simulations. The columns are built by repeatedly placing centrally aligned stars at a distance of 5.0 A ˚ from each other. Later, during the minimization process and MD simulations, the molecules find their preferred distances. Using the Forcite module of Materials Studio [70], the preformed clusters are placed into a rectangular box with side lengths of Lx = LB = 81 A ˚ and Lz = 140 A ˚ for the initial setup. Subsequently, to achieve an overall density of 1.0 g cm−3, a total of 28978 water molecules are randomly and homogeneously distributed to solvate the columnar cluster. Finally, the configurations are transferred to LAMMPS to simulate each system for 10 ns. All simulations are performed in the NPT ensemble (Nosé–Hoover barostat and thermostat) at the temperature T = 298 K and pressure P = 1 atm while using a time step of 1 fs. More details on the MD protocols, polymer consistent force field (PCFF) parameters, etc. can be found in ref. [29].

4.2. Computation of the Binding Energy Diagrams and the Binding Energy as a Function of Assembly Size N

The definition of the binding energy is given by Equation (1). The diagrams for the dimers (Figure 2 and Figure 3) are obtained using the M06-L/aug-cc-pVTZ functional [71] and basis set with a semi-empirical counterpoise-type correction for basis set superposition error, as implemented in Gaussian 09 [72]. The binding of larger stacks from either MD simulations or experimental data (Figure 5) is quantified applying density functional based tight-binding with parameter set mio-1-1 [73].
The cumulative parameter Σ (Figure 6) is calculated as a function of the molecule index i within the stack by using the group/group command of LAMMPS [74]. For each molecule in the column-shaped cluster the pairwise interactions with all the other azobenzene stars are calculated, excluding any self-interactions within molecule i. More precisely, considering a molecule with index i, all pairwise Coulombic and Lennard–Jones energy contributions coming from the N-1 surrounding molecules (with indices ≠ i) are summed up to Σ (i). The parameter Σ (i) therefore describes the strength of binding of the i-th azo star with the entire rest of the cluster. Please note, that the sum of Σ (i) over all molecules i in the stack is clearly not equal to the binding energy E B of a cluster as defined above. This occurs, since for any molecules i and j in the stack there are reappearing energy contributions to Σ (i) or Σ (j), respectively, for example from another nearby molecule k. The results for the cumulative parameter are obtained by averaging over an ensemble of three independent simulations of 10 ns for each cluster size.

4.3. Hirshfeld Surface Analysis of the azo–BTA Star Crystal

To more closely examine the intermolecular interactions in the crystals of azo-stars experimentally characterized by Lee et al. [17], the Hirshfeld surfaces [65,66] are generated using CrystalExplorer (Version 3.1) [75]. This method has been successfully utilized to characterize aminoazobenzene polymorphs [76], to investigate the intermolecular interactions in azobenzenes revealing photoinduced crystal–melt transition [77], and to explain the molecular origin of the color changes of guest signaling azobenzene [78,79], to name a few.
The Hirshfeld surface is a three-dimensional isosurface effectively determining the volume of the molecule in the crystal. The surface itself is defined so that for each point on this surface, one-half of the electron density is due to the spherically averaged noninteracting atoms comprising the molecule, and the other half is due to those comprising the rest of the crystal [66]. The quantity which is mapped onto the surface is d norm (normalized contact distance), which is the negative or positive sum of d i and d e when the contacts are shorter or longer than the sum of the van der Waals radii of the interacting atoms, respectively:
d norm = ( d i r i vdW ) r i vdW + ( d e r e vdW ) r e vdW .
The fingerprint plots (de-di diagrams) are constructed using d e and d i distances from the Hirshfeld surface to the nearest atom outside (external) and inside (internal) the surface, correspondingly. These plots provide a summary of the intermolecular contacts in the crystal.

4.4. Quantification of Hydrogen Bonds, Close Contacts and Macrodipoles

The length of the hydrogen bond lHB is considered as a separation between the hydrogen and oxygen atoms of the amide groups from the neighboring stars in the stack applying the structural definitions of the bond suggested by Luzar and Chandler [80]. The close contacts lCC are defined in terms of the sum of the specified atomic radii. Finally, the dipole moments μ N of the stars and stacks with N≤ 8 are calculated using density functional theory with M06-L/aug-cc-pVTZ hybrid functional, as justified by Zawada et al. for hydrogen-bonded molecular chains [81]. For comparison and for larger aggregates, the point-like charges on atoms from PCFF parametrization are used to calculate the normalized dipole moments as μ N / μ max , where μ max is the maximal value of the dipole molecule of the stack with N = 20. All the aforementioned properties are analyzed using the Forcite module of BIOVIA Materials Studio [70].

Author Contributions

Conceptualization, O.G. and A.S.P.; methodology, M.K., A.S.P. and O.G.; formal analysis, M.K. and A.S.P.; investigation, V.S., M.K. and O.G.; data curation, V.S. and M.K.; writing—original draft preparation, O.G.; writing—review and editing, O.G. and M.S.; visualization, V.S. and O.G.; supervision, M.S. and O.G.; project administration, O.G.; funding acquisition, O.G.

Funding

This research was funded by German Research Foundation (DFG), grant numbers GU1510/3-1, SA1657/13-1 and GU1510/5-1. V.S. and O.G. thank the European Regional Development Fund and the Free State of Saxony for the service contract within the European Social Fund project 100231947 (Young Investigators Group “Computer Simulations for Materials Design—CoSiMa”, Dresden Center for Computational Materials Science, Technische Universität Dresden).

Acknowledgments

We would like to acknowledge high-performance computing support from the Center for Information Services and High Performance Computing (ZIH) provided by Technische Universität Dresden.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BTABenzene-1,3,5-tricarboxamide
azo–BTA starThree-arm star-shaped molecule with BTA core and azobenzene arms
DFTDensity functional theory
DFTB+Density functional based tight-binding
MDMolecular dynamics
PCFFPolymer consistent force field
LAMMPSLarge-scale atomic/molecular massively parallel simulator

References

  1. Ikeda, T.; Tsutsumi, O. Optical switching and image storage by means of azobenzene liquid-crystal films. Science 1995, 268, 1873–1875. [Google Scholar] [CrossRef] [PubMed]
  2. Beharry, A.A.; Wooley, G.A. Azobenzene photoswitches for biomolecules. Chem. Soc. Rev. 2011, 40, 4422–4437. [Google Scholar] [CrossRef]
  3. Gillono, M.; Roppolo, I.; Frascella, F.; Scaltrito, L.; Pirri, C.F.; Chiappone, A. CO2 permeability control in 3D printed light responsive structures. Appl. Mater. Today 2019, 100470. [Google Scholar] [CrossRef]
  4. Serak, S.V.; Tabiryan, N.V.; White, T.J.; Bunning, T.J. Azobenzene liquid crystal polymer-based membrane and cantilever optical systems. Opt. Express 2009, 17, 15736–15746. [Google Scholar] [CrossRef] [PubMed]
  5. Jelken, J.; Pandiyarajan, C.K.; Genzer, J.; Lomadze, N.; Santer, S. Fabrication of flexible hydrogel sheets featuring periodically spaced holes with continuously adjustable size in real time. ACS Appl. Mater. Interfaces 2018, 10, 30844–30851. [Google Scholar] [CrossRef] [PubMed]
  6. Wang, D.; Huang, Y.; Lv, J.-M.; Min, Y.; Quan, Y.-Y.; Ge, L.-N.; Tian, M.; Yao, D.-S. Multi-arm azobenzene liquid crystal based on cholic acid: Synthesis and mesophase properties. Liq. Cryst. 2018, 45, 1813–1824. [Google Scholar] [CrossRef]
  7. Vialetto, J.; Groppi, J.; La Rosa, M.; Silvi, S.; Credi, A.; Baroncini, M. Solution and solid state photochromism in a family of shape persistent azobenzene tetramers functionalized with alkyloxy substituents. Photochem. Photobiol. Sci. 2019, 18, 2281–2286. [Google Scholar] [CrossRef]
  8. Ooi, Y.-H.; Yeap, G.-Y.; Han, C.-C.; Lin, H.-C.; Shinomiya, T.; Ito, M.M. Non-conventional three-armed star-shaped mesogens based on 1, 3, 5-trisubstituted benzene with azobenzene moieties at the periphery: Synthesis, and mesomorphic behaviour. Liq. Cryst. 2014, 41, 1017–1033. [Google Scholar] [CrossRef]
  9. Takahashi, T.; Tanino, T.; Ando, H.; Nakano, H.; Shirota, Y. Surface relief grating formation using a novel azobenzene-based photochromic amorphous molecular material, tris [4-(phenylazo) phenyl] amine. Mol. Cryst. Liq. Cryst. 2005, 430, 9–14. [Google Scholar] [CrossRef]
  10. Yadavalli, N.S.; Saphiannikova, M.; Santer, S. Photosensitive response of azobenzene containing films towards pure intensity or polarization interference patterns. Appl. Phys. Lett. 2014, 105, 051601. [Google Scholar] [CrossRef]
  11. Jiang, X.; Lu, J.; Zhou, F.; Zhang, Z.; Pan, X.; Zhang, W.; Wang, Y.; Zhou, N.; Zhu, X. Molecularly-defined macrocycles containing azobenzene main-chain oligomers: Modular stepwise synthesis, chain-length and topology-dependent properties. Polym. Chem. 2016, 7, 2645–2651. [Google Scholar] [CrossRef]
  12. Imrie, C.T.; Henderson, P.A. Liquid crystal dimers and higher oligomers: Between monomers and polymers. Chem. Soc. Rev. 2007, 36, 2096–2124. [Google Scholar] [CrossRef] [PubMed]
  13. Imrie, C.T.; Henderson, P.; Yeap, G.-Y. Liquid crystal oligomers: Going beyond dimers. Liq. Cryst. 2009, 36, 755–777. [Google Scholar] [CrossRef]
  14. Koch, M.; Saphiannikova, M.; Santer, S.; Guskova, O. Photoisomers of azobenzene star with a flat core: Theoretical insights into multiple states from DFT and MD perspective. J. Phys. Chem. B 2017, 121, 8854–8867. [Google Scholar] [CrossRef]
  15. Galanti, A.; Santoro, J.; Mannancherry, R.; Duez, Q.; Diez-Cabanes, V.; Valášek, M.; De Winter, J.; Cornil, J.; Gerbaux, P.; Mayor, M.; et al. A new class of rigid multi (azobenzene) switches featuring electronic decoupling: Unravelling the isomerization in individual photochromes. J. Am. Chem. Soc. 2019, 141, 9273–9283. [Google Scholar] [CrossRef]
  16. Han, M.; Nakanishi, H. Light-sensitive microspheres based on spherical assembly of star-shaped chromophores. J. Photopolym. Sci. Technol. 2018, 31, 527–531. [Google Scholar] [CrossRef]
  17. Lee, S.; Oh, S.; Lee, J.; Malpani, Y.; Jung, Y.-S.; Kang, B.; Lee, J.Y.; Ozasa, K.; Isoshima, T.; Lee, S.Y.; et al. Stimulus-responsive azobenzene supramolecules: Fibers, gels, and hollow spheres. Langmuir 2013, 29, 5869–5877. [Google Scholar] [CrossRef]
  18. Lee, J.; Oh, S.; Pyo, J.; Kim, J.-M.; Je, J.H. A light-driven supramolecular nanowire actuator. Nanoscale 2015, 7, 6457–6461. [Google Scholar] [CrossRef]
  19. Bahrenburg, J.; Sievers, C.M.; Schönborn, J.B.; Hartke, B.; Renth, F.; Temps, F.; Näther, C.; Sönnichsen, F.D. Photochemical properties of multi-azobenzene compounds. Photochem. Photobiol. Sci. 2013, 12, 511–518. [Google Scholar] [CrossRef]
  20. Adachi, H.; Hirai, Y.; Ikeda, T.; Maeda, M.; Hori, R.; Kutsumizu, S.; Haino, T. Photoresponsive toroidal nanostructure formed by self-assembly of azobenzene-functionalized tris(phenylisoxazolyl) benzene. Org. Lett. 2016, 18, 924–927. [Google Scholar] [CrossRef]
  21. Baroncini, M.; d’Agostino, S.; Bergamini, G.; Ceroni, P.; Comotti, A.; Sozzani, P.; Bassanetti, I.; Grepioni, F.; Hernande, T.M.; Silvi, S.; et al. Photoinduced reversible switching of porosity in molecular crystals based on star-shaped azobenzene tetramers. Nat. Chem. 2015, 7, 634. [Google Scholar] [CrossRef] [PubMed]
  22. Nacci, C.; Baroncini, M.; Credi, A.; Grill, L. Reversible photoswitching and isomer-dependent diffusion of single azobenzene tetramers on a metal surface. Angew. Chem. Int. Ed. 2018, 130, 15254–15259. [Google Scholar] [CrossRef]
  23. Akiyama, H. Photochemically reversible liquefaction/solidification of sugar–alcohol derivatives. In Functional Organic Liquids; Wiley–VCH Verlag GmbH Co.: Weinheim, Germany, 2019; pp. 75–86. [Google Scholar]
  24. Liao, L.-X.; Junge, D.M.; McGrath, D.V. Photochromic dendrimers containing six azobenzenes. Macromolecules 2002, 35, 319–322. [Google Scholar] [CrossRef]
  25. Tanaka, D.; Ishiguro, H.; Shimizu, Y.; Uchida, K. Thermal and photoinduced liquid crystalline phase transitions with a rod–disc alternative change in the molecular shape. J. Mater. Chem. 2012, 22, 25065–25071. [Google Scholar] [CrossRef]
  26. Choi, Y.-J.; Kim, D.-Y.; Park, M.; Yoon, W.-J.; Lee, Y.; Hwang, J.-K.; Chiang, Y.-W.; Kuo, S.-W.; Hsu, C.-H.; Jeong, K.-U. Self-assembled hierarchical superstructures from the benzene-1, 3, 5-tricarboxamide supramolecules for the fabrication of remote-controllable actuating and rewritable films. ACS Appl. Mater. Interfaces 2016, 8, 9490–9498. [Google Scholar] [CrossRef] [PubMed]
  27. Segarra-Maset, M.D.; van Leeuwen, P.W.N.M.; Freixa, Z. Light switches the ligand! Photochromic azobenzene–phosphanes. Eur. J. Inorg. Chem. 2010, 2010, 2075–2078. [Google Scholar] [CrossRef]
  28. Liu, J.; Zhang, Q.; Zhang, J.; Hou, W. Photochemical behaviours of star-like liquid crystal with azobenzene terminal groups. J. Mater. Sci. 2005, 40, 4517–4521. [Google Scholar] [CrossRef]
  29. Koch, M.; Saphiannikova, M.; Guskova, O. Do columns of azobenzene stars disassemble under light illumination? Langmuir 2019, 35, 14659–14669. [Google Scholar] [CrossRef] [PubMed]
  30. Wolfer, P.; Audorff, H.; Kreger, K.; Kador, L.; Schmidt, H.-W.; Stingelin, N.; Smith, P. Photo-induced molecular alignment of trisazobenzene derivatives. J. Mater. Chem. 2011, 21, 4339–4345. [Google Scholar] [CrossRef]
  31. Wolfer, P.; Kreger, K.; Schmidt, H.-W.; Stingelin, N.; Smith, P. Photo-oriented trisazobenzene layers for patterned liquid-crystal alignment. Mol. Cryst. Liq. Cryst. 2012, 562, 133–140. [Google Scholar] [CrossRef]
  32. Kreger, K.; Wolfer, P.; Audorff, H.; Kador, L.; Stingelin-Stutzmann, N.; Smith, P.; Schmidt, H.-W. Stable holographic gratings with small-molecular trisazobenzene derivatives. J. Am. Chem. Soc. 2010, 132, 509–516. [Google Scholar] [CrossRef] [PubMed]
  33. Wang, J.; Wang, S.; Zhou, Y.; Wang, X.; He, Y. Fast photoinduced large deformation of colloidal spheres from a novel 4-arm azobenzene compound. ACS Appl. Mater. Interfaces 2015, 7, 16889–16895. [Google Scholar] [CrossRef] [PubMed]
  34. Sperner, M.; Tober, N.; Detert, H. Tristriazolotriazines with azobenzene arms-acidochromic dyes and discotic liquid crystals. Eur. J. Org. Chem. 2019, 2019, 4688–4693. [Google Scholar] [CrossRef] [Green Version]
  35. Goldmann, D.; Janietz, D.; Schmidt, C.; Wendorff, J.H. Liquid crystalline 1, 3, 5-triazines incorporating rod-like azobenzene sub-units. Liq. Cryst. 1998, 25, 711–719. [Google Scholar] [CrossRef]
  36. Pokladek, Z.; Dudek, M.; Mongin, O.; Metivier, R.; Mlynarz, P.; Samoc, M.; Matczyszyn, K.; Paul, F. Linear and third-order nonlinear optical properties of triazobenzene-1, 3, 5-triazinane-2, 4, 6-trione (isocyanurate) derivatives. ChemPlusChem 2017, 82, 1372–1383. [Google Scholar] [CrossRef]
  37. Wu, S.L.; Zhao, F.Y.; Luo, S.W.; Zhang, W.J.; Xu, S.D.; Zeng, Q.D. Synthesis and scanning tunneling microscopy observation of 1, 3, 5-triazine bearing azo-carboxylate. J. Nanosci. Nanotechnol. 2017, 17, 725–728. [Google Scholar] [CrossRef]
  38. Jin, H.; Jian, T.; Ding, Y.-H.; Chen, Y.; Mu, P.; Wang, L.; Chen, C.-L. Solid-phase synthesis of three-armed star-shaped peptoids and their hierarchical self-assembly. Biopolymers 2019, 110, e23258. [Google Scholar] [CrossRef]
  39. Vera, F.; Tejedor, R.M.; Romero, P.; Barberá, J.; Ros, M.B.; Serrano, J.L.; Sierra, T. Light-driven supramolecular chirality in propeller-like hydrogen-bonded complexes that show columnar mesomorphism. Angew. Chem. Int. Ed. 2007, 46, 1873–1877. [Google Scholar] [CrossRef]
  40. Metzroth, T.; Hoffmann, A.; Martín-Rapún, R.; Smulders, M.M.J.; Pieterse, K.; Palmans, A.R.A.; Vekemans, J.A.J.M.; Meijer, E.W.; Spiess, H.W.; Gauss, J. Unravelling the fine structure of stacked bipyridine diamine-derived C3-discotics as determined by X–ray diffraction, quantum-chemical calculations, fast-MAS NMR and CD spectroscopy. Chem. Sci. 2011, 2, 69–76. [Google Scholar] [CrossRef]
  41. Van Herrikhuyzen, J.; Jonkheijm, P.; Schenning, A.P.H.J.; Meijer, E.W. The influence of hydrogen bonding and ππ stacking interactions on the self-assembly properties of C3–symmetrical oligo (p-phenylenevinylene) discs. Org. Biomol. Chem. 2006, 4, 1539–1545. [Google Scholar] [CrossRef]
  42. Narayan, B.; Kulkarni, C.; George, S.J. Synthesis and self-assembly of a C3–symmetric benzene-1, 3, 5-tricarboxamide (BTA) anchored naphthalene diimide disc. J. Mater. Chem. C 2013, 1, 626–629. [Google Scholar] [CrossRef]
  43. Smulders, M.M.J.; Schenning, A.P.H.J.; Meijer, E.W. Insight into the mechanisms of cooperative self-assembly: The “sergeants-and-soldiers” principle of chiral and achiral C3–symmetrical discotic triamides. J. Am. Chem. Soc. 2008, 130, 606–611. [Google Scholar] [CrossRef] [PubMed]
  44. Filot, I.A.W.; Palmans, A.R.A.; Hilbers, P.A.J.; van Santen, R.A.; Pidko, E.A.; de Greef, T.F.A. Understanding cooperativity in hydrogen-bond-induced supramolecular polymerization: A density functional theory study. J. Phys. Chem. B 2010, 114, 13667–13674. [Google Scholar] [CrossRef] [PubMed]
  45. Kulkarni, C.; Meijer, E.W.; Palmans, A.R.A. Cooperativity scale: A structure–mechanism correlation in the self-assembly of benzene-1, 3, 5-tricarboxamides. Acc. Chem. Res. 2017, 50, 1928–1936. [Google Scholar] [CrossRef] [PubMed]
  46. Bochicchio, D.; Pavan, G.M. Molecular modelling of supramolecular polymers. Adv. Phys. X 2018, 3, 316–338. [Google Scholar] [CrossRef] [Green Version]
  47. Malpani, Y.R.; Oh, S.; Lee, S.; Jung, Y.-S.; Kim, J.-M. Photoinduced phase transition of azobenzene-coupled benzenetricarboxamide. Bull. Korean Chem. Soc 2014, 35, 2563–2566. [Google Scholar] [CrossRef] [Green Version]
  48. Devi, S.; Bala, I.; Gupta, S.P.; Kumar, P.; Pal, S.K.; Venkataramani, S. Reversibly photoswitchable alkoxy azobenzenes connected benzenetricarboxamide discotic liquid crystals with perpetual long range columnar assembly. Org. Biomol. Chem. 2019, 17, 1947–1954. [Google Scholar] [CrossRef]
  49. Choi, Y.-J.; Yoon, W.-J.; Park, M.; Kang, D.-G.; Bang, G.; Koo, J.; Lim, S.-I.; Park, S.; Jeong, K.-U. Construction of light-responsive phase chirality from an achiral macrogelator. J. Mater. Chem. C 2019, 7, 3231–3237. [Google Scholar] [CrossRef]
  50. Leyendecker, M. Supramolecular Lyotropic Liquid Crystalline Alignment Media Based on Benzene-1,3,5-tricarboxamides; Technische Universität Darmstadt: Darmstadt, Germany, 2019. [Google Scholar]
  51. Huang, Y.; Kim, D.-H. Light-controlled synthesis of gold nanoparticles using a rigid, photoresponsive surfactant. Nanoscale 2012, 4, 6312–6317. [Google Scholar] [CrossRef]
  52. Goldenberg, L.M.; Kulikovsky, L.; Kulikovska, O.; Tomczyk, J.; Stumpe, J. Thin layers of low molecular azobenzene materials with effective light-induced mass transport. Langmuir 2010, 26, 2214–2217. [Google Scholar] [CrossRef]
  53. Yadavalli, N.S. Advances in Experimental Methods to Probe Surface Relief Grating Formation Mechanism in Photosensitive Materials; Universität Potsdam: Potsdam, Germany, 2014. [Google Scholar]
  54. Wolfer, P. Ordering Molecular Organics into Functionality; ETH Zürich: Zürich, Germany, 2011. [Google Scholar]
  55. Taniguchi, T.; Asahi, T.; Koshima, H. Photomechanical azobenzene crystals. Crystals 2019, 9, 437. [Google Scholar] [CrossRef] [Green Version]
  56. De Greef, T.F.A.; Smulders, M.M.J.; Wolffs, M.; Schenning, A.P.H.J.; Sijbesma, R.P.; Meijer, E.W. Supramolecular polymerization. Chem. Rev. 2009, 109, 5687–5754. [Google Scholar] [CrossRef] [PubMed]
  57. Lightfoot, M.P.; Mair, F.S.; Pritchard, R.G.; Warren, J.E. New supramolecular packing motifs: π-stacked rods encased in triply-helical hydrogen bonded amide strands. Chem. Commun. 1999, 1945–1946. [Google Scholar] [CrossRef]
  58. Nakano, Y.; Hirose, T.; Stals, P.J.M.; Meijer, E.W.; Palmans, A.R.A. Conformational analysis of supramolecular polymerization processes of disc-like molecules. Chem. Sci. 2012, 3, 148–155. [Google Scholar] [CrossRef]
  59. Ranganathan, D.; Kurur, S.; Gilardi, R.; Karle, I.L. Design and synthesis of AB3-type (A = 1,3,5-benzenetricarbonyl unit; B = Glu diOMe or Glu7 octa OMe) peptide dendrimers: Crystal structure of the first generation. Biopolym. Orig. Res. Biom. 2000, 54, 289–295. [Google Scholar] [CrossRef]
  60. Gong, B.; Zheng, C.; Yan, Y. Structure of N,N′,N″-tris(carboxymethyl)-1,3,5-benzenetricarboxamide trihydrate. J. Chem. Crystallogr. 1999, 29, 649–652. [Google Scholar] [CrossRef]
  61. Albuquerque, R.Q.; Timme, A.; Kress, R.; Senker, J.; Schmidt, H.-W. Theoretical investigation of macrodipoles in supramolecular columnar stackings. Chem. Eur. J. 2013, 19, 1647–1657. [Google Scholar] [CrossRef]
  62. Wegner, M.; Dudenko, D.; Sebastiani, D.; Palmans, A.R.A.; de Greef, T.F.A.; Graf, R.; Spiess, H.W. The impact of the amide connectivity on the assembly and dynamics of benzene-1,3,5-tricarboxamides in the solid state. Chem. Sci. 2011, 2, 2040–2049. [Google Scholar] [CrossRef]
  63. Greciano, E.E.; Calbo, J.; Buendía, J.; Cerdá, J.; Aragó, J.; Ortí, E.; Sánchez, L. Decoding the consequences of increasing the size of self-assembling tricarboxamides on chiral amplification. J. Am. Chem. Soc. 2019, 141, 7463–7472. [Google Scholar] [CrossRef]
  64. Korlepara, D.B.; Henderson, W.R.; Castellano, R.K.; Balasubramanian, S. Differentiating the mechanism of self-assembly in supramolecular polymers through computation. Chem. Commun. 2019, 55, 3773–3776. [Google Scholar] [CrossRef]
  65. Spackman, M.A.; Jayatilaka, D. Hirshfeld surface analysis. CrystEngComm 2009, 11, 19–32. [Google Scholar] [CrossRef]
  66. Spackman, M.; McKinnon, J.J. Fingerprinting intermolecular interactions in molecular crystals. CrystEngComm 2002, 4, 378–392. [Google Scholar] [CrossRef]
  67. Garzoni, M.; Baker, M.B.; Leenders, C.M.A.; Voets, I.K.; Albertazzi, L.; Palmans, A.R.A.; Meijer, E.W.; Pavan, G.M. Effect of H–bonding on order amplification in the growth of a supramolecular polymer in water. J. Am. Chem. Soc. 2016, 138, 13985–13995. [Google Scholar] [CrossRef] [PubMed]
  68. Bochicchio, D.; Kwangmettatam, S.; Kudernac, T.; Pavan, G.M. How defects control the out-of-equilibrium dissipative evolution of a supramolecular tubule. ACS Nano 2019, 13, 4322–4334. [Google Scholar] [CrossRef] [PubMed]
  69. Tomasi, J.; Mennicci, B.; Cammi, R. Quantum mechanical continuum solvation models. Chem. Rev. 2005, 105, 2999–3093. [Google Scholar] [CrossRef] [PubMed]
  70. BIOVIA Materials Studio. Dassault Systéms; Version 8; BIOVIA Materials Studio: San Diego, CA, USA, 2014. [Google Scholar]
  71. Zhao, Y.; Truhlar, D.G. A new local density functional for main group thermochemistry, transition metal bonding, thermochemical kinetics, and noncovalent interactions. J. Chem. Phys. 2006, 125, 194101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09, Revision D.01; Gaussian, Inc.: Wallingford, UK, 2009. [Google Scholar]
  73. Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, T.; Suhai, S.; Seifert, G. Self-consistent-charge density-functional tight-binding method for simulations of complex materials properties. Phys. Rev. B Condens. Matter Phys. 1998, 58, 7260–7268. [Google Scholar] [CrossRef]
  74. Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 1995, 117, 1–19. [Google Scholar] [CrossRef] [Green Version]
  75. Wolff, S.K.; Grimwood, D.J.; McKinnon, J.J.; Turner, M.J.; Jayatilaka, D.; Spackman, M.A. CrystalExplorer (Version 3.1); University of Western Australia: Perth, WA, Australia, 2012. [Google Scholar]
  76. Cruickshank, D.L.; Hendon, C.H.; Verbeek, M.J.R.; Walsh, A.; Wilson, C.C. Polymorphism of the azobenzene dye compound methyl yellow. CrystEngComm 2016, 18, 3456–3461. [Google Scholar] [CrossRef] [Green Version]
  77. Hoshino, M.; Uchida, E.; Norikane, Y.; Azumi, R.; Nozawa, S.; Tomita, A.; Sato, T.; Adachi, S.-I.; Koshihara, S.-Y. Crystal melting by light: X-ray crystal structure analysis of an azo crystal showing photoinduced crystal–melt transition. J. Am. Chem. Soc. 2014, 136, 9158–9164. [Google Scholar] [CrossRef]
  78. Scott, J.L.; Almesaker, A.; Sumi, Y.; Tanaka, K. Guest signaling compound: Trans-3,3′-bis(diphenylhydroxymethyl)azobenzene. Cryst. Growth Des. 2007, 7, 1049–1054. [Google Scholar] [CrossRef]
  79. Scott, J.L.; Tanaka, K. Signalling by modulation of intermolecular interactions. In Engineering of Crystalline Materials Properties. NATO Science for Peace and Security Series B: Physics and Biophysics; Springer: Dordrecht, The Netherlands, 2008; pp. 429–447. [Google Scholar]
  80. Luzar, A.; Chandler, D. Hydrogen-bond kinetics in liquid water. Nature 1996, 379, 55–57. [Google Scholar] [CrossRef]
  81. Zawada, A.; Kaczmarek-Kędziera, A.; Bartkowiak, W. On the potential application of DFT methods in predicting the interaction-induced electric properties of molecular complexes. Molecular H-bonded chains as a case of study. J. Mol. Model. 2012, 18, 3073–3086. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. The object of the study: (a) the azo–benzene-1,3,5-tricarboxamide (BTA) star with benzene-1,3,5-tricarboxamide core and three azobenzene arms. Each arm carries the N, N-dimethylamino group. Here, the optimized geometry is shown [14]. A snapshot of the columnar aggregate with N = 20 from MD trajectory [29] is illustrated in panel (b). The azobenzene star studied experimentally by Lee et al. [17] is depicted in panel (c). The conformation of this star is the same as in the crystalline structure, shown in panel (d). Here, a = 59.37 A ˚ , b = 9.42 A ˚ and c = 39.95 A ˚ ; α = β = 90° and γ = 103.16°. The original unit cell [17] is multiplied twice along the b-direction to illustrate the stacking of the stars. The panel (e) shows the stack made of 20 such molecules. The carbon, hydrogen, oxygen and nitrogen atoms are colored dark grey, light grey, red and blue spheres, respectively.
Figure 1. The object of the study: (a) the azo–benzene-1,3,5-tricarboxamide (BTA) star with benzene-1,3,5-tricarboxamide core and three azobenzene arms. Each arm carries the N, N-dimethylamino group. Here, the optimized geometry is shown [14]. A snapshot of the columnar aggregate with N = 20 from MD trajectory [29] is illustrated in panel (b). The azobenzene star studied experimentally by Lee et al. [17] is depicted in panel (c). The conformation of this star is the same as in the crystalline structure, shown in panel (d). Here, a = 59.37 A ˚ , b = 9.42 A ˚ and c = 39.95 A ˚ ; α = β = 90° and γ = 103.16°. The original unit cell [17] is multiplied twice along the b-direction to illustrate the stacking of the stars. The panel (e) shows the stack made of 20 such molecules. The carbon, hydrogen, oxygen and nitrogen atoms are colored dark grey, light grey, red and blue spheres, respectively.
Molecules 24 04387 g001
Figure 2. Binding energy EB depending on the π π distance d between the star cores in the dimer. The literature values [17,44,48,57,58,59,60] of the π π distances of BTA-cored stars (in some of these publications called dendrimers) with various arms are shown as symbols with ordinate EB = −100 kJ mol−1 for convenience. In case of the slipped stacks from the experimental data [17], the stars were shifted or rotated so that they segregated in order to evaluate a true π π distance between their centres. The definition of d is given by the sketch.
Figure 2. Binding energy EB depending on the π π distance d between the star cores in the dimer. The literature values [17,44,48,57,58,59,60] of the π π distances of BTA-cored stars (in some of these publications called dendrimers) with various arms are shown as symbols with ordinate EB = −100 kJ mol−1 for convenience. In case of the slipped stacks from the experimental data [17], the stars were shifted or rotated so that they segregated in order to evaluate a true π π distance between their centres. The definition of d is given by the sketch.
Molecules 24 04387 g002
Figure 3. Binding energy EB depending on the l shift and the rotation angle Ω in star dimer (d = 3.6 A ˚ , k = 0 A ˚ ) (a) and as a function of l, k shifts (d = 3.6 A ˚ , Ω = 0°) (b). The experimental coordinates [17] are shown as a black bullet on the top projection (b). The definitions of d, k, l and Ω are given on schemes at the bottom.
Figure 3. Binding energy EB depending on the l shift and the rotation angle Ω in star dimer (d = 3.6 A ˚ , k = 0 A ˚ ) (a) and as a function of l, k shifts (d = 3.6 A ˚ , Ω = 0°) (b). The experimental coordinates [17] are shown as a black bullet on the top projection (b). The definitions of d, k, l and Ω are given on schemes at the bottom.
Molecules 24 04387 g003
Figure 4. Hirshfeld surface analysis [65] of the azo–BTA star crystal studied by Lee et al. [17] and the fingerprint plots [66] (de-di diagrams) of intermolecular interactions between oxygen (a), nitrogen (b), carbon (c), hydrogen (d) and other atoms in the molecular crystal. The percentage of the corresponding interaction is given in parentheses. The Hirshfeld surface is shown as a semitransparent surface with a reference molecule inside. The contacting atoms from the neighboring stars in the crystal are shown as beads. For the color scheme see Figure 1.
Figure 4. Hirshfeld surface analysis [65] of the azo–BTA star crystal studied by Lee et al. [17] and the fingerprint plots [66] (de-di diagrams) of intermolecular interactions between oxygen (a), nitrogen (b), carbon (c), hydrogen (d) and other atoms in the molecular crystal. The percentage of the corresponding interaction is given in parentheses. The Hirshfeld surface is shown as a semitransparent surface with a reference molecule inside. The contacting atoms from the neighboring stars in the crystal are shown as beads. For the color scheme see Figure 1.
Molecules 24 04387 g004
Figure 5. Binding energy EB per pair as a function of assembly size N calculated for the star stack from the experimental crystal cell [17] (a) and for the columns modelled in MD simulations [29] (b). The dashed red line on panel (b) guides for the eye.
Figure 5. Binding energy EB per pair as a function of assembly size N calculated for the star stack from the experimental crystal cell [17] (a) and for the columns modelled in MD simulations [29] (b). The dashed red line on panel (b) guides for the eye.
Molecules 24 04387 g005
Figure 6. Cumulative parameter Σ calculated for each molecule in a stack depending on its position along the column of size N.
Figure 6. Cumulative parameter Σ calculated for each molecule in a stack depending on its position along the column of size N.
Molecules 24 04387 g006
Figure 7. Average length of the intermolecular hydrogen bonds lHB (a) and average length of the close contacts lCC (b) as functions of assembly size N. The MD data from ref. [29] are taken for the analysis. The blue line in panel (a) shows the length of the hydrogen bonds as measured from the crystal cell studied by Lee et al. [17]. The dashed black line on panel (a) guides for the eye. The red line in panel (b) shows an average value for the close contacts as measured from the crystal cell studied by Lee et al. [17] and the red hatched field on the graph depicts the error bar above and below the average value.
Figure 7. Average length of the intermolecular hydrogen bonds lHB (a) and average length of the close contacts lCC (b) as functions of assembly size N. The MD data from ref. [29] are taken for the analysis. The blue line in panel (a) shows the length of the hydrogen bonds as measured from the crystal cell studied by Lee et al. [17]. The dashed black line on panel (a) guides for the eye. The red line in panel (b) shows an average value for the close contacts as measured from the crystal cell studied by Lee et al. [17] and the red hatched field on the graph depicts the error bar above and below the average value.
Molecules 24 04387 g007
Figure 8. Average dipole moment μ N (a) and normalized dipole moment μ N / μ max (b) as functions of assembly size N. The results shown on both panels as triangles refer to the structural data from ref. [17]. The other two data sets (red rhombi and open blue squares) are taken from snapshots of the MD trajectories [29]. The dashed blue line on panel (b) guides for the eye. In panel (a), error bars are depicted for the two MD data sets.
Figure 8. Average dipole moment μ N (a) and normalized dipole moment μ N / μ max (b) as functions of assembly size N. The results shown on both panels as triangles refer to the structural data from ref. [17]. The other two data sets (red rhombi and open blue squares) are taken from snapshots of the MD trajectories [29]. The dashed blue line on panel (b) guides for the eye. In panel (a), error bars are depicted for the two MD data sets.
Molecules 24 04387 g008
Table 1. Mean intermolecular hydrogen bond length lHB in A ˚ as a function of the position inside the stack. The indices i and j of Nij (1–2, 2–3, etc.) denote the location of the hydrogen bonds, i.e., the index numbers of hydrogen-bonded molecules in the stack of size N = 20.
Table 1. Mean intermolecular hydrogen bond length lHB in A ˚ as a function of the position inside the stack. The indices i and j of Nij (1–2, 2–3, etc.) denote the location of the hydrogen bonds, i.e., the index numbers of hydrogen-bonded molecules in the stack of size N = 20.
NijN1–2N2–3N3–4N4–5N5–6N6–7N7–8N8–9N9–10
lHB2.342.142.162.222.151.932.152.002.02
NijN10–11N11–12N12–13N13–14N14–15N15–16N16–17N17–18N18–19N19–20
lHB2.182.082.172.051.992.042.112.082.212.25

Share and Cite

MDPI and ACS Style

Savchenko, V.; Koch, M.; Pavlov, A.S.; Saphiannikova, M.; Guskova, O. Stacks of Azobenzene Stars: Self-Assembly Scenario and Stabilising Forces Quantified in Computer Modelling. Molecules 2019, 24, 4387. https://doi.org/10.3390/molecules24234387

AMA Style

Savchenko V, Koch M, Pavlov AS, Saphiannikova M, Guskova O. Stacks of Azobenzene Stars: Self-Assembly Scenario and Stabilising Forces Quantified in Computer Modelling. Molecules. 2019; 24(23):4387. https://doi.org/10.3390/molecules24234387

Chicago/Turabian Style

Savchenko, Vladyslav, Markus Koch, Aleksander S. Pavlov, Marina Saphiannikova, and Olga Guskova. 2019. "Stacks of Azobenzene Stars: Self-Assembly Scenario and Stabilising Forces Quantified in Computer Modelling" Molecules 24, no. 23: 4387. https://doi.org/10.3390/molecules24234387

APA Style

Savchenko, V., Koch, M., Pavlov, A. S., Saphiannikova, M., & Guskova, O. (2019). Stacks of Azobenzene Stars: Self-Assembly Scenario and Stabilising Forces Quantified in Computer Modelling. Molecules, 24(23), 4387. https://doi.org/10.3390/molecules24234387

Article Metrics

Back to TopTop