We use cookies to ensure that we give you the best experience on our website. By continuing to browse this repository, you give consent for essential cookies to be used. You can read more about our Privacy and Cookie Policy.

Durham Research Online
You are in:

Thermodynamics of the self-assembly of non-ionic chromonic molecules using atomistic simulations. The case of TP6EO2M in aqueous solution.

Akinshina, A. and Walker, M. and Wilson, M. R. and Tiddy, G. J. T. and Masters, A. J. and Carbone, P. (2015) 'Thermodynamics of the self-assembly of non-ionic chromonic molecules using atomistic simulations. The case of TP6EO2M in aqueous solution.', Soft matter., 11 (4). pp. 680-691.


Atomistic molecular dynamic simulations have been performed for the non-ionic chromonic liquid crystal 2,3,6,7,10,11-hexa-(1,4,7-trioxa-octyl)-triphenylene (TP6EO2M) in aqueous solution. TP6EO2M molecules consist of a central poly-aromatic core (a triphenylene ring) functionalized by six hydrophilic ethyleneoxy (EO) chains, and have a strong tendency to aggregate face-to-face into stacks even in very dilute solution. We have studied self-assembly of the molecules in the low concentration range corresponding to an isotropic solution of aggregates, using two force fields GAFF and OPLS. Our results reveal that the GAFF force field, even though it was successfully used previously for modelling of ionic chromonics, overestimates the attraction of TP6EO2M molecules in water. This results in an aggregation free energy which is too high, a reduced hydration of EO chains and, therefore, molecular self-assembly into compact disordered clusters instead of stacks. In contrast, use of the OPLS force field, leads to self-assembly into ordered stacks in agreement with earlier experimental studies of triphenylene-based chromonics. The free energy of association follows a “quasi-isodesmic” pattern, where the binding free energy of two molecules to form a dimer is of the order of 2.5 RT larger than the corresponding energy of addition of a molecule into a stack. The obtained value for the binding free energy, ΔG0agg = −12 RT, is found to be in line with the published values for typical ionic chromonics (−7 to −12 RT), and agrees reasonably well with the experimental results for this system. The calculated interlayer distance between the molecules in a stack is 0.37 nm, which is at the top of the range found for typical chromonics (0.33–0.37 nm). We suggest that the relatively large layer spacing can be attributed to the repulsion between EO side chains.

Item Type:Article
Full text:(AM) Accepted Manuscript
Download PDF
Publisher Web site:
Date accepted:24 November 2014
Date deposited:24 February 2015
Date of first online publication:25 November 2014
Date first made open access:No date available

Save or Share this output

Look up in GoogleScholar