Cookies

We use cookies to ensure that we give you the best experience on our website. You can change your cookie settings at any time. Otherwise, we'll assume you're OK to continue.

Department of Chemistry

Prof. Mark R. Wilson


Researcher ID profile

Personal web page

Professor in the Department of Chemistry
Telephone: +44 (0) 191 33 42144

(email at mark.wilson@durham.ac.uk)

Research Interests

We use classical molecular dynamics and Monte Carlo simulation methods to study problems in soft matter chemistry. We concentrate on materials that are ordered on the nanoscale, including liquid crystals, colloids, polymers and amphiphilic molecules.
Read about our latest simulations at my research group web pages.

Simulation of chromonic liquid crystals

Sunset yellow structure, chromonic column of sunset yellow molecules and chromonic liquid crystal phase

We have recently carried out the first detailed study of the structure and dynamics of a chromonic liquid crystal (edicol - the food dye sunset yellow).1 The simulations show spontaneous self-assembly in water at atomistic detail, and can be used to determine the free energy of association in solution. Sunset yellow is seen to form single molecule stacks with a preference for head-to-tail ordering, which at higher concentrations form liquid crystal phases in water (as shown in the picture above). Our work shows that simulation can be used to understand, predict and control molecular self-assembly with a view to the future use of chromonic materials for applications in self-assembled molecular electronics.

Simulation of liquid crystal molecules and calculation of material properties

Molecular simulations provide a way of understanding the structure and dynamics3 of complex materials such as liquid crystals. We have recently used4 atomistic simulation to study the structure of the first low molecular weight biaxial nematic liquid crystal, ODBP-Ph-C7.

We have also been developing techniques to simulate liquid crystal molecules at the atomistic level and predict material properties in liquid crystal phases from simulation.5 For example, the rotational viscosity determines how quickly a nematic liquid crystal can reorientate in a liquid crystal display; and can be determined from our simulations by monitoring the liquid crystal director in the course of a molecular dynamics simulation.5 To carry out this work, we need accurate representations of intramolecular interactions within a molecule and intermolecular interactions between molecules. This is made possible by a new force field for use in materials chemistry, which we have been developing.6

Growth of a biaxial phase from a uniaxial starting point showing also the formation of ferroelectric domains (blue and red colours)

The picture shows a sequence of time frames showing the slow growth of a biaxial phase from a uniaxial starting point. As the phase grows it develops ferroelectric domains, which are shown colour-coded. Molecules which are blue and red have transverse dipoles pointing in opposite directions.

A snapshot from a molecular dynamics simulation showing the liquid crystal molecule PCH5 in a nematic phase.

Dynamics in soft solids

motion in an octafluronaphthalene crystal

We have been using atomistic simulation to study dynamics in molecular solids.2 Our recent work with octafluoronaphthalene and urea inclusion compounds has indicated that its is now possible to look at dynamical processes over tens of ns, and study correlated motion in solids. This provides the possiblity of using simulation to complement and help interpret experimental solid state NMR data, and raises intriguing questions as to whether we can use simulation to  help engineer tractable molecular machines for future applications.

Simulation studies of amphiphilic molecules in water

We are interested in the properties of amphiphilic molecules and have recently used simulation to understand the behaviour of amphiphilic polymers at a water/air interface. The simulations can be used to interpret the results of complementary neutron diffraction studies by calculating neutron reflectivity curves to match with experiment.7

Similar modelling techniques are being used to study lipid - amino acid interactions. Our recent atomistic work in this area seeks to understand how these interactions (specifically the balance of hydrogen bonding and cation-Pi interactions) can be perturbed by small chemical modifications of the amino acid structure. Such studies seek to aid our understanding of binding between proteins and cell membranes.

Snapshots from a simulation study of an amphiphilic polymer at a water-air interface at low surface concentration. The hydrophobic polynorbornene backbone is shown in blue and the hydrophilic poly(ethylene oxide) (PEO) grafts are shown in red. The dots represent water molecules. At higher surface concentrations the PEO chains are forced down into the aqueous subphase to form a polymer brush.

Coarse-grained models for the simulation of polymers, dendrimers, liquid crystals and lipids

We have been developing coarse-grained simulation methods that allow us understand structure and organisation in a range of molecular materials. We have studied liquid crystalline polymers and dendrimers8 and are working on “molecular engineering” tools to design molecules to have the desired organisation at the nanoscale. This research is being extended to lipid bilayers and vesicles. An important part of our work has involved developing efficient simulation methods that make use of parallel computers, so we are able to simulate extremely large systems composed of many molecules and take advantage of today’s supercomputers.9
The figure shows microphase separation in a liquid crystalline side chain polymer to give polymer-rich and liquid crystal-rich regions. Blue: polymer backbone (methylsiloxane), white: liquid crystal group, red: flexible spacer groups (alkyl chain).

More on Soft Matter Simulations

Read about our latest simulations at my research group web pages

Picture of a polymer nanowire

References

  1. Chami F. and Wilson M. R. Molecular Order in a Chromonic Liquid Crystal: A Molecular Simulation Study of the Anionic Azo Dye Sunset Yellow. J. Am. Chem. Soc., 2010, 132, 7794. 
  2. Ilott A. J., Palucha S., Batsanov A. S., Wilson M. R. and Hodgkinson P.  Elucidation of Structure and Dynamics in Solid Octafluoronaphthalene from Combined NMR, Diffraction, And Molecular Dynamics Studies.   J. Am. Chem. Soc., 2010, 132, 5179.
  3. Oganesyan, V. S., Kuprusevicius, E., Gopee, H., Cammidge, A. N., Wilson, M. R. Simulation of EPR spectra directly from molecular dynamics trajectories of a liquid crystal with doped paramagnetic spin probe. Phys. Rev. Lett., 2009, 102, 013005(1)-013005(4).
  4. Pelaez J. and Wilson M. R. Atomistic Simulations of a Thermotropic Biaxial Liquid Crystal. Phys. Rev. Lett., 2006, 97., 267801.
  5. Cheung D. L., Clark S. J, Wilson M. R. Calculation of the rotational viscosity of a nematic liquid crystal. Chem. Phys. Lett., 2002, 356, 140.
  6. Cheung D. L., Clark S. J., Wilson M. R. Parametrization and validation of a force field for liquid-crystal forming molecules. Phys. Rev. E, 2002, 65,051709.
  7. Anderson P. M. and Wilson M. R. Molecular Dynamics Simulations of an Amphiphilic Graft Copolymer at a Water/Air Interface., 2004, J. Chem. Phys, 121, 8503.
  8. Wilson M. R., Ilnytskyi J. M., Stimson L. M., Computer simulations of a liquid crystalline dendrimer in liquid crystalline solvents. J. Chem. Phys., 2003, 119, 3509.
  9. Ilnytskyi J, Wilson M. R.. A domain decomposition molecular dynamics program for the simulation of flexible molecules with an arbitrary topology of Lennard-Jones and/or Gay-Berne sites. Comput. Phys. Comm., 2001, 134, 23.
  10. Earl D. J., M., Wilson M. R. Predictions of molecular chirality and helical twisting powers: A theoretical study. J. Chem. Phys., 2003, 119, 10280.


Read our latest work at our publications page.