Substituent interference on supramolecular assembly in urea gelators : synthesis , structure prediction and NMR †

Eighteen N-aryl-N0-alkyl urea gelators were synthesised in order to understand the effect of head substituents on gelation performance. Minimum gelation concentration values obtained from gel formation studies were used to rank the compounds and revealed the remarkable performance of 4-methoxyphenyl urea gelator 15 in comparison to 4-nitrophenyl analogue 14, which could not be simply ascribed to substituent effects on the hydrogen bonding capabilities of the urea protons. Crystal structure prediction calculations indicated alternative low energy hydrogen bonding arrangements between the nitro group and urea protons in gelator 14, which were supported experimentally by NMR spectroscopy. As a consequence, it was possible to relate the observed differences to interference of the head substituents with the urea tape motif, disrupting the order of supramolecular packing. The combination of unbiased structure prediction calculations with NMR is proposed as a powerful approach to investigate the supramolecular arrangement in gel fibres and help understand the relationships between molecular structure and gel formation.


Introduction
Supramolecular gels, formed by the self-assembly of molecularscale building blocks, have a wide range of applications. 13][4] The late development of this area can be explained by the serendipitous nature of the discovery of new gelators. 5,6espite much effort, the gelation process of low molecular weight organogelators is still not fully understood. 7In particular, the correlation between precursor composition, assembly conditions and self-assembled structure remains unclear.As a consequence, many researchers have identified new gelators simply by modifying known scaffolds. 8,9However, the approach of taking known gelators and tuning them for a specific purpose is limited by the systems available and solvents required for effective gelation, which may further limit the scope of application. 5 increased understanding of the relationship between molecular and self-assembled structure would not only represent a significant scientific advance, but also enable the design of new gel-based materials. 10t around the same time, Hanabusa, 11 van Esch and Kellogg 7,[12][13][14] have independently exploited urea derivatives to obtain new gelators, where the urea motif represents the molecular equivalent of a Velcro s type hook-and-loop. 15This enabled molecular unidirectional self-association and, for this reason, the urea motif has been extensively used to form supramolecular assemblies.7][18] With this knowledge, functional groups began to be used as a means to manipulate the solubility parameters of gelators (e.g.long alkyl chain in polar solvents), allowing access to this metastable state and taking advantage of the modified solute-solvent interactions. 19These and many other studies during the past twenty years contributed to the formulation of three basic rules 1 for designing low molecular weight gelators: (1) the existence of unidirectional intermolecular interactions to favour one-dimensional self-assembly, (2) intertwining of the aggregates for network formation, and (3) the control of fibre/solvent interfacial energy to manipulate solubility and prevent the crystallisation of the metastable gel.
The original aim of this work was to extend the understanding of molecular structure-gelation relationships by systematically varying some molecular features of a common scaffold, as exemplified in Fig. 1, in order to isolate the contribution of -NO 2 and -OCH 3 head substituents on a phenyl ring.Gelation test outcomes have driven the series of structural modifications that led to the synthesis of a total of eighteen analogues.As a consequence of the unusual observations made and in an effort to understand them, two representative gelators (14 and 15) underwent more comprehensive thermo-mechanical characterisation, molecular modelling and NMR studies.

Synthesis of the gelators
Compounds 1-18 were synthesised and all products precipitated as solids and were isolated in moderate to good yields, ESI.†

Gelation studies
Gel formation studies were performed in fifteen solvents, chosen in order to cover a wide range of polarities from hexane to water.Gelators were tested at nine different concentrations (5, 10, 15, 20, 30, 40, 50, 100, 150 mg mL À1 ).The threshold concentration of 20 mg mL À1 was that from which all gel formation studies commenced.If a gel was formed at this concentration, further tests at lower concentration (15, 10, 5 mg mL À1 ) were performed.If the result was the formation of a partial gel, more tests were conducted at higher concentrations (30, 40, 50, 100, 150 mg mL À1 ) in an attempt to form a gel.No further tests were carried out when no gel or insolubility was observed (ESI †).Although a wide range of solvents were tested, solvent polarity was found not to directly affect (or at least straightforwardly) gelation performance, despite reported evidence on other systems. 20For this reason, it was decided to focus on three solvents: toluene, tetralin and dimethyl sulfoxide (DMSO).Gelation performances of 1-18 were then ranked with respect to the minimum gelation concentration (MGC) values in these three solvents (Table 1).
Gelators 1-3 bearing short alkyl chains did not form any gels due to their high solubility, indicating that solute-solvent interactions prevailed over the solute-solute interactions. 21,22vidently, the relatively short alkyl chain (C 6 ) was insufficient to provide enough solvophobic drive, with a weak contribution to chain packing through van der Waals interactions.
In an effort to optimise gelation performance, the number of carbon atoms in the chain was increased (C 11 -C 14 ) in 4-15.Notably, it was more difficult to obtain gels in toluene in the presence of an odd numbered alkyl chain, while the presence of the even numbered C 14 chain yielded the best performance in both toluene and tetralin.Additionally, no gels were formed in dimethyl sulfoxide until the alkyl chain was sufficiently long.Consequently, C 14 alkyl chain analogues were chosen to study the substituent interference on the self-assembly of the class of molecules investigated in this work.
Since the inclusion of the C 14 apolar moiety in 13-15 was shown to be sufficient to promote either gelation or lower MGC values and literature reports show that long n-alkanes can provide gels of remarkable stability in organic solvents, 23 gelator 18 was designed and found to exhibit low MGC values (15-20 mg mL À1 ).
Finally, to confirm the importance of the urea hydrogen bonding motif for this scaffold, one or both urea hydrogen bonding donor sites were removed in 16 (amide) and 17 (ester).Poor gelation observed was attributed to the either extremely low or high solubility of 16 and 17, respectively, incompatible with the required balance between solute-solute and solutesolvent interactions. 21nspection of the MCG values for gelators 13-15 indicated that the presence of the head substituent -OCH 3 led to better performance (5 mg mL À1 ) than -NO 2 (20 mg mL À1 ).Gelation driven by cooperative hydrogen bonding between urea functional groups in gelators 13-15 was expected to be reinforced by electron-withdrawing substitution of the aromatic head group leading to more acidic NH groups and consequently stronger hydrogen bonding interactions, with improved gelation.However, the reverse was found and gel formation studies on gelators 14 and 15 showed that the presence of electron donating -OCH 3 could improve the gelation performances in two of the three selected solvents.We note that the same effect is seen in comparing 11 to 12, while with shorter chain lengths (5 vs. 6) there is some evidence that the -NO 2 substituent improves gelation (Table 1).These observations indicate that there must be a balance of effects that influence the gelation behaviour, related to alkyl chain length, substituent electron donating or withdrawing ability, and perhaps other, less intuitive, substituent effects.

Thermo-mechanical and morphological characterisation of the gels
Studies were conducted to explore whether a simple correlation existed between the gelation performance of a compound and the thermo-mechanical characteristics of the resulting gel.The gels derived from 14 (-NO 2 ) and 15 (-OCH 3 ) in toluene were selected for more detailed studies.
The gel-sol transition temperatures (T gel ) were obtained by differential scanning calorimetry (DSC) measurements (ESI †).Each gel sample (E20 mg of gels at 20 mg mL À1 ) underwent a heating-cooling-heating cycle from 30 to 180 1C at a rate of 10 1C min À1 in sealed aluminium pans.In both samples, characteristic thermal events were only present in the first heating ramp.Their T gel were found to be 75 1C for gel 14 and 86 1C for gel 15.
Rheological tests were performed (ESI †) with both gels showing a solid-like behaviour, where the storage modulus (G 0 , E10 3 Pa) was systematically larger than the loss modulus (G 00 , E10 2 Pa) by roughly one order of magnitude and parallel to each other with respect to the angular frequency.Their rigid behaviour was further proved in oscillation sweep experiments, where the phase angle d, that measures the delay between stress and strain, was found always around 101 (an ideal elastic solid has d = 01).Additionally, the phase angle showed no significant variation across the frequency range investigated, indicating the ability of all gels to withstand stresses below their yield stress value without showing signs of mechanical fatigue. 24,25The oscillation stress sweep experiments demonstrated that the strength of the material was increasing from gel 14 to gel 15, as confirmed by both relative order and position of the curves.It was not possible to extrapolate the yield stress values at d = 451 since G 0 and G 00 data sets never crossed each other, as the structure broke down suddenly before the final data point could be collected.
Rheology measurements, in accordance with DSC outcomes, showed the existence of a correlation between the gelation performance and the thermo-mechanical characteristics of the gel materials.4-Nitrophenyl urea 14, a gelator with higher MGC values and therefore poorer tendency to self-assemble, also resulted in a less stable network and consequently in a weaker material when compared to 4-methoxy phenyl urea derived gel 15.
Furthermore, environmental scanning electron microscopy (ESEM) was used to investigate morphological differences in dried gel samples, defined as xerogels, obtained from gels 14 and 15 in toluene (ESI †).Two kinds of morphologies were observed in Fig. 2, possibly reflecting the head groups' influence in the supramolecular packing of the corresponding gelators during the formation of the network.Interestingly, xerogel 15 presented a more fibrillar network, when compared to xerogel 14.

Molecular modelling
To help interpret the observed differences in gelation between 14 and 15, electronic-structure calculations were performed on both isolated molecules, as well as crystal structure prediction (CSP) calculations.We applied CSP methods here to take advantage of the developments that have been made in this area over the past few years [26][27][28][29] in predicting the preferred solid-state assembly of organic molecules using a global search of the lattice energy surface.This global optimisation approach provides an unbiased assessment of molecular packing possibilities in a solid state environment, unlike other recent molecular modelling studies of gel structure based on assumed supramolecular arrangements 30,31 conformational sampling of small molecular clusters, 32 or dynamical simulations of pre-arranged fibre models. 33,34hile CSP calculations have been applied to fairly challenging molecules, 35,36 due to the size and possible conformational flexibility of these molecules, the extent to which their crystal packing possibilities could be sampled is limited.Moreover, CSP generates perfect crystalline arrangements of the molecules, which might not correspond to the packing that is present in gel fibres.Therefore, the realistic goal of this study was not to determine the packing of molecules in the gel fibres, but to understand the molecules' preferred intermolecular arrangements, and the influence on solid-state packing of changing the substitution of the aromatic ring.
The initial motivation for investigating different aromatic substitutions stemmed from the potential ability of electron withdrawing or donating substituents to tune the acidity of the urea group.To investigate this inductive effect, single molecule density functional theory (DFT) calculations were performed using DMol3 37 at the B3LYP/DNP level of theory to investigate the electrostatic potential around 14 and 15.These calculations revealed a slightly more positive electrostatic potential at the polar hydrogen atoms of the urea in 14 compared to 15, as shown in Fig. 3, and a more negative electrostatic potential in the p-electron region of the aromatic ring for 15.These differences reflect the expected influence of the electron withdrawing and donating substituents.However, the more conspicuous difference between the molecules was the introduction of a region of negative charge at the -NO 2 substituent in 14.The oxygen atoms in the -NO 2 group of 14 are particularly accessible, so are positioned to easily compete with the urea oxygen as hydrogen bond acceptors for the protons of the urea group.Similarly, the oxygen of the -OCH 3 substituent in gelator 15 also introduced a region of negative charge, but one that was much less strong and was shielded to some extent by the steric bulk of the methyl.These details of the electrostatic potential are significant, as they could affect the disruption of the urea ''zipping'', which was associated with good gelation properties. 15n a study of the effects of anions on urea gelators, Steed et al. has reported that only anions of sufficient hardness can disrupt this pattern. 38In comparing 14 to 15, the increased acidity of the urea hydrogen atoms could enhance the formation of urea tapes in 14, while the presence of competing hydrogen bond acceptors might have disrupted their formation.The purpose of performing CSP calculations was to explore the balance of these competing effects and the overall influence of the aromatic substituents on the preferred arrangement of these molecules.
Putative crystal structures were generated using the Global Lattice Energy Explorer software, 39,40 performing quasirandom crystal structure generation, followed by lattice energy minimisation.An initial study was performed with rigid, gas phase molecular geometries, generating 5000 crystal structures in each of 16 common space groups for organic molecules, for both 14 and 15 (rigid searches were performed with two orientations of the -OCH 3 group on 15, which differed by only 1.2 kJ mol À1 ).All crystal structures were lattice energy minimised using a force field model combining an empirically parameterised model of exchange-repulsion and dispersion interactions with an accurate atomic multipole model for intermolecular electrostatics. 41This rigid molecule approach allowed a fairly rapid exploration of crystal packing phase space, due to the reduced dimensionality of the search space and requirement only to consider intermolecular interactions in assessing the relative stabilities of the predicted structures.
These initial calculations displayed a clear difference in the predicted crystal packing behaviours of gelators 14 and 15.The low energy predicted structures of gelator 15 displayed hydrogen bonding between urea groups, predominantly forming the one-dimensional urea chains or tapes (Fig. 4d and f) that were associated with gel fibre formation.The -OCH 3 did not participate in hydrogen bonding.In contrast, the low energy structures of 14 lacked the N-HÁ Á ÁO urea hydrogen bonding linking urea into one-dimensional chains or tapes; the CSP results suggested that the urea tape formation was disrupted by the presence of the -NO 2 group, which preferentially acted as a hydrogen bond acceptor amongst the low energy predicted crystal structures (Fig. 4b).
To verify that these findings were not biased by the simplification of using rigid molecular geometries in the crystal structure search, further CSP calculations were performed with flexibility of the molecular geometry allowed during structure generation and lattice energy minimisation.These more demanding calculations provide the opportunity for molecules to explore packing arrangements that are inaccessible to the molecules in their gas phase geometries, as well as to optimise their hydrogen bonding networks, which can be very sensitive to small changes in molecular geometry.A further 5000 crystal structures were generated and minimised in each of the same 16 space groups for each molecule, with flexibility in the two dihedral angles around N-C bonds connecting the urea to the aromatic ring and alkyl chain.Flexibility of the dihedral connecting the -NO 2 or -OCH 3 group to the aromatic ring was also included.Molecular distortion corresponding to up 22 kJ mol À1 was allowed during the search, as this is the largest intramolecular strain seen in crystalline molecular geometries of small organic molecules. 42nlike intermolecular interactions, atom-atom force fields to describe intramolecular energies lack the accuracy required for CSP.Modern methods combine atom-atom intermolecular models with a DFT description of intramolecular geometry and energy. 27Here, we applied a newly developed method in which the intramolecular energies and atomic partial charges were modelled with a Gaussian process regression model fitted to DFT data. 43Details of the model, structure generation and energy minimisation are provided as ESI.† All structures generated with the flexible search had their final energies calculated with the same anisotropic atom-atom force field intermolecular energy model used in the rigid structure search.
These flexible-molecule calculations did not change the overall conclusions from the initial calculations: hydrogen bonded urea tapes dominate the low energy structures of 15, while the low energy structures of 14 are dominated by urea hydrogen bonded to -NO 2 .However, the details of these interactions were modified by giving the molecule freedom to distort from its ideal geometry.For 14, a slight reorientation of the alkyl chain allowed the alignment of both -NO 2 oxygen atoms with urea hydrogen atoms, forming a double hydrogen bond, as shown in Fig. 4a.This interaction was found in the overall lowest energy predicted structures of 14 and provided a significant stabilisation.The total energy of the best structure involving the double hydrogen bond (Fig. 4a) is 9 kJ mol À1 more stable than the lowest energy structure involving only one -NO 2 oxygen in hydrogen bonding (Fig. 4b).The lowest energy predicted crystal structure of 14 that displayed the urea hydrogen bond tape motif was a further 13 kJ mol À1 higher in energy (Fig. 4c), outside the energetic range that would be observable in small molecule crystal structures. 44This gives an indication of the strength of the ability of the oxygen atoms on the -NO 2 group to interfere with the urea-tape packing motif.
For gelator 15, the influence of molecular flexibility on the predictions was subtler; here, freedom of the molecule to distort from its gas phase geometry led to some low energy structures with planar urea tapes (Fig. 4e), whereas these interactions were twisted or buckled in the best structures generated with the rigid gas phase geometry, as per Fig. 4d and f.The difference in lattice energies of the most stable predicted structures of each (À198.8kJ mol À1 for 14 vs.À176.7 kJ mol À1 for 15) further highlighted the strength of the ureaÁ Á ÁO 2 N interaction.The calculated lattice energies of the lowest energy urea-tape structures of each molecule, as shown in Fig. 4c vs. Fig.4d, were remarkably similar.
There are various types of interactions in these systems, including hydrogen bonding and van der Waals interactions, which can be maximised by stacking of the aromatic rings and aligning molecular tails.All of the low energy structures in the CSP were densely packed.It was seen, though, that in the systems of gelator 15 all of the types of interactions were associated within a pair of neighbouring molecules, i.e. they had a tendency to align tails, urea groups, and aromatic groups in series, a behaviour that is thought to be associated with the formation of a self-assembled fibrillar network (SAFIN) in gels. 4,45For gelator 14, there was instead a strong interaction between the urea and the -NO 2 group of molecules, but these two molecules could not then align their tails and arrange aromatic groups into p-stack (Fig. 4a and b).These CSP results suggested that, rather than the acidity of the urea group determining the difference in behaviour of 14 and 15, it was such compromises between sources of strong interaction that may cause the differing observations in the gel formation experiments.

Powder X-ray diffraction
Semicrystalline structures attributed to the gel fibres were verified by the presence of diffraction peaks PXRD patterns (ESI †) acquired on gel samples without the need to remove any amount of solvent.
Comparison between the simulated and observed diffraction patterns must take into account the expected discrepancies between the lattice parameters of modelled structures and packing in the gel fibres, which are usually a few percent in lattice parameters, as well as differences in reflection intensities caused by orientation of the gel fibres.Nevertheless, several of the predicted crystal structures for each molecule were identified as giving similar simulated diffraction patterns to those observed from the gels (Fig. 5).These comparisons support the conclusions from the structure prediction: for 15, those structures that give similar PXRD patterns all include the hydrogen bond urea tapes, while urea-nitro hydrogen bonding is present in all structures of 14 that give good agreement with the observed diffraction pattern.The structures that give good agreement with the experimental PXRD patterns include some of the lowest energy predicted structures of each molecule (Fig. 4a and e, À198.8 kJ mol À1 for 14 and À176.8 kJ mol À1 for 15).
NMR chemical shielding calculations were performed on all of the predicted structures that gave good agreement with gel PXRD patterns, vide infra.

NMR spectroscopy
NMR spectroscopy enabled us to get further insight into the molecular level organisation of the gels.This technique has already proven its potential in characterising structure and dynamics of soft materials. 46,47These multiphasic materials required both solution and solid-state NMR spectroscopy to probe mobile entities and rigid components of the gel network, respectively.NMR studies were conducted with gels prepared from dimethyl sulfoxide to avoid overlapping of aromatic peaks with solvent peaks from toluene and tetralin.NOESY solution-state NMR experiments provided information on spatial connectivity between molecular regions in aggregated states. 7Cross-peaks in 2D 1 H-1 H NOESY (ESI †) spectra of gels 14 and 15 displayed the same phase as the diagonal peaks, as shown in Fig. 6, indicative of negative NOE enhancements, which are characteristic of large molecules.Since these systems are exclusively formed of small gelator molecules (o600 Da), these negative enhancements can be attributed to medium-large supramolecular aggregates, As shown in Fig. 6, cross-peaks between the -CH 3 protons of the alkyl chain and the aromatic and urea ones were identified in gel 14, whereas they were not detected for gel 15.The presence of such cross-peaks in gel 14 revealed close proximity between the end of the aliphatic chain, the benzene ring and the urea moiety (Table 2), belonging either to the same or to surrounding molecules.Such an arrangement was suggested by the CSP calculations, whose results showed that the urea HÁ Á ÁO 2 N interaction leads to interdigitation of the aliphatic chains, as shown in Fig. 4a and b, which places the end of the chain near to the aromatic rings of neighbouring molecules.Furthermore, intermolecular cross-peaks between distant NH a and H c or H d (Fig. 6) suggested proximity of the aromatic ring to another molecule with a different orientation.Conversely, the absence of cross-peaks between the alkylic -CH 3 protons, the aromatic, the urea and the -OCH 3 sites in gel 15 suggested a well-aligned network in which the long aliphatic chains are positioned separately from aromatic and urea groups.In addition, the absence of cross-peaks between NH a and H c or H d (Fig. 6) was also in agreement with the predicted well defined stacking of the aromatic rings in the 3D network of gel 15.
1 H- 13 C CP/MAS NMR experiments 48 can probe rigid components of the gel samples and were therefore used to study the molecular packing of the gel fibres.Despite the presence of high contents of solvent which was expected to lead to an increased mobility and a reduced efficiency of polarisation transfer, we were able to acquire 1 H- 13 C CP/MAS NMR spectra of gels 14 and 15.The detection of peaks in 1 H- 13 C CP/MAS NMR spectra of these gels was indicative of a rigid 3D fibrous network, with densely packed molecules as predicted by molecular modelling studies.To the best of our knowledge, this is the first example of supramolecular gels in which 1 H- 13 C CP/MAS NMR peaks were detectable without any physical   modification of the gel samples, i.e. removing any amount of the solvent or lowering the temperature. 49or gels 14 and 15, aliphatic carbons appeared as sharp resonances, while low intensity peaks were observed for the aromatic and carbonyl carbons (and -OCH 3 carbon in gel 15) in spectra acquired at 25 1C and 1 kHz MAS rate, as shown in Fig. 7a (the assignment is provided separately, ESI †).The presence of sharp peaks is in agreement with the semicrystalline gel fibres identified in PXRD experiments.Lowering the temperature to 15, 10, 5 and 0 1C enabled us to record 1 H- 13 C CP/MAS NMR spectra at higher MAS rates, as shown in Fig. 7b-g.This led to significant narrowing of the aromatic and carbonyl carbons (and -OCH 3 carbon in gel 15), with spinning sidebands present in accordance with the specific MAS rate.No significant differences were observed for the chemical shift values in the spectra acquired at lower temperatures.
The 1 H- 13 C CP/MAS NMR spectrum of frozen gel 14 acquired at an MAS rate of 8 kHz (Fig. 7g) presented a single peak per carbon, indicative of one magnetic environment for each carbon site.Peaks at 125.7 and 117.2 ppm resulted from overlapped resonances of carbons C 2 and C 6 , and C 3 and C 5 , respectively (Fig. 8), as predicted by solid-state DFT calculations.The variation of chemical shift values between the reference solid powder 14 (Fig. 7h) and the resulting gel suggested the molecular organisation of the gel fibres was different from the reference solid powder.
In contrast to gel 14, frozen gel 15 acquired at an MAS rate of 8 kHz (Fig. 7g) showed several peaks per carbon in the aromatic region, indicating the same carbon can experience different magnetic environments in different molecules.These results suggested the presence of multiple symmetrically independent gel molecules in the structure.Structures with multiple independent molecules (Z 0 4 1) are not uncommon in molecular crystals and it is not surprising that this phenomenon carries over to gel fibre packing. 50This observation reflects the symmetry of molecular packing and does not imply a less organised structure.
However, this finding does mean that the packing in the gel 15 fibres cannot be perfectly described by any one of the predicted crystal structures, since the molecular modelling studies were conducted considering only one molecule in the asymmetric unit (Z 0 = 1).Furthermore, the similarity in chemical shift values of the reference solid powder 15 (Fig. 7h) and the resulting gel pointed towards a self-assembled network with a 3D organisation similar to that of the crystalline powder.
For comparison, NMR chemical shielding calculations were performed for the predicted crystal structures of gelators 14 and 15 and gave reasonable agreement with the gel X-ray diffraction pattern, ESI † (Fig. 5), using DFT calculations (with the PBE functional and a plane wave basis set, ESI †) and the GIPAW 51 methodology.Despite the similarities in intermolecular interactions, we found a reasonable variation in predicted chemical shifts from the selected structures (see ESI †).Amongst these structures, the simulated 13 C chemical shifts for the lowest energy predicted structures (À198.8kJ mol À1 for 14, reported in Fig. 4a, and À173.8 kJ mol À1 for 15, reported in Fig. 4e) were found to be in best agreement with the 13 C chemical shift values observed for gel fibres (Fig. 9), with RMSD (predicted -observed) of 1.8 ppm for 14 and 1.4 ppm for 15.These deviations are as low as those seen in comparing observed 13 C chemical shifts with those predicted from known single crystal structures of small organic molecules. 52These results strongly support the relevance of the predicted structures in interpreting the properties of the two gelators.The increased resolution obtained for the aromatic and carbonyl carbons (and -OCH 3 carbon in gel 15) upon spinning at higher MAS rates, with no significant changes for the aliphatic carbons, strongly indicated that chemical shift anisotropy is considerable for the aromatic and carbonyl carbons (and -OCH 3 carbon in gel 15), even in the gel state.The solid-state DFT calculations of NMR properties (ESI †) performed for the low energy predicted crystal structures confirmed that the anisotropy of the shielding tensors for the aromatic and carbonyl carbon atoms was much larger (in the range from 108 to 192 ppm) than for the aliphatic carbons (in the range from 19 to 43 ppm) for both gelators, supporting the previous conclusions.The -OCH 3 carbon in gel 15 was narrower than the aromatic carbons, but the peak sharpened as a result of spinning faster, in agreement with the DFT calculated intermediate anisotropy of the shielding tensor (ca.70 ppm).Hence, the significantly broadened peaks detected for aromatic and carbonyl carbons at low MAS rates in 1 H- 13 C CP/ MAS NMR spectra of gels 14 and 15 were most likely due to the high values of chemical shift anisotropy for these carbon sites.
Overall, the combination of calculation, diffraction and NMR studies provides valuable insight into the molecular level organisation, offering likely motifs of intermolecular packing and local structure in the gels.

Conclusions
The aim of this work was to study the influence of the electronwithdrawing (-NO 2 ) and electron-donating (-OCH 3 ) head substituents on the gelation performances of urea gelators.The thermo-mechanical characterisation of gels 14 (20 mg mL À1 ) and 15 (5 mg mL À1 ) confirmed the latter to be a stronger material, expressing higher T gel and more solid-like rheological behaviour.The results of the crystal structure prediction calculations suggested that, rather than changes in the acidity of the urea group, the differences observed in their gelation performance is due to the presence of strong competing interactions between the -NO 2 head substituent and the urea.The computational studies of ideal crystals quantified these interactions, and we therefore suggest that the availability of substituent groups to form competing interactions must be considered in any programme of rationally designed self-assembling systems.The crystalline arrangements produced by the modelling studies and their molecular contacts detected by NOESY solution-state NMR experiments showed that gelator 15 has a tendency to align tails, urea groups, and aromatic rings in series.Conversely, such well stacked supramolecular assemblies were not observed in gelator 14, due to the perturbing influence of the -NO 2 group on the formation of the urea tape motif.These materials were proven to be rigid by 1 H- 13 C CP/MAS solid-state NMR studies and therefore some insight into preferred molecular arrangements was gained, particularly through comparison with simulated chemical shifts of the predicted structures.Despite the limitations of molecular modelling studies of gels, when validated by solid-state NMR experiments, these have a great potential in unravelling structure of supramolecular gels.calculated values for the predicted structures that give best agreement with both NMR and PXRD from the gel.Calculated isotropic chemical shieldings were converted to chemical shifts by matching the calculated and observed chemical shift of the CH 3 carbon in gel.The structures for gelators 14 and 15 correspond to those reported in Fig. 4a and e, respectively.

Fig. 1
Fig. 1 Schematic representation of the urea scaffold investigated.

Fig. 3
Fig. 3 Calculated molecular electrostatic potentials of 14 (a-c) and 15 (d-f).The potentials are plotted on the 0.02 a.u.electronic isodensity surface.The red patches of positive electrostatic potential in (b) and (e) correspond to the urea hydrogen atom positions.(c and f) Show the view down the urea oxygen atom, highlighting its negative electrostatic potential.

Fig. 4
Fig. 4 Molecular packing in representative low energy predicted crystal structures of 14 (a-c) and 15 (d-f).(a and d) Are the global minimum energy crystal structures that have been located for each molecule.Calculated lattice energies of each structure are shown.Hydrogen bonds are shown as dashed blue lines.Hydrogen atoms on the alkyl chains are hidden in (f) for clarity.

Fig. 5
Fig. 5 Powder X-ray diffraction patterns of the gels (a) 14 and (b) 15 in DMSO, compared to simulated diffraction patterns from low energy predicted structures that give similar diffraction patterns to the gels.The simulated patterns from the predicted structures are labelled by their calculated lattice energies.The simulated diffraction patterns of 14 and 15 with lattice energies of À198.8 and À173.8 kJ mol À1 correspond to the predicted crystal structures reported in Fig. 4a and e, respectively.A larger set of comparisons is provided in the ESI.†

Fig. 6
Fig.62D 1 H-1 H NOESY spectra of gels 14 and 15 (60 mg mL À1 in DMSO-d 6 ) acquired at 25 1C with a mixing time of 0.5 s using a 800.23 MHz spectrometer.Relevant cross-peaks have the corresponding spatial connectivity assigned (coloured lines).The cross-peaks observed between the -CH 2 and the aromatic protons in gel 15 do not contain structural information, as they result from indirect NOE enhancements, only observable at long mixing times (ESI †).

aa
Pairs of protons for which cross-peaks were observed only at mixing times longer than 0.25 s, thus preventing accurate calculations of distances.

Fig. 9
Fig.9Experimental13 C chemical shift values for frozen gels 14 and 15 derived from the spectra acquired at 0 1C and an MAS rate of 8 kHz vs. calculated values for the predicted structures that give best agreement with both NMR and PXRD from the gel.Calculated isotropic chemical shieldings were converted to chemical shifts by matching the calculated and observed chemical shift of the CH 3 carbon in gel.The structures for gelators 14 and 15 correspond to those reported in Fig.4a and e, respectively.

Table 2
Intermolecular distances calculated from initial rates of NOE build-up curves with mixing time, using the H c -H d distance as reference