Durham University Logo           

Soft Matter Simulations in the Wilson Group


Contact Details
Prof. Mark Wilson
Department of Chemistry
University Science Laboratories
South Road
Durham
DH1 3LE
United Kingdom
Tel: +44 (0) 191 334 2144
Fax: +44 (0) 191 384 4737
Email: mark.wilson@durham.ac.uk
                      Jump to Liquid Crystal Simulations
Jump to Amphiphilic Molecule Simulations
Jump to Polymer and Dendrimer Simulations
Jump to Coarse-grained Simulations
Jump to Chiral System Simulations
Jump to Parallel Simulations
      Link to Wilson Group Home Page
Link to Mark Wilson Home Page
Link to Mark Wilson Publications Page
Link to Mark Wilson Publications Page
Link to Chemistry Department Home Page
Link to Durham University Home Page

Research Summary

We use classical molecular dynamics and Monte Carlo simulation methods to study problems in soft matter chemistry, concentrating on materials that are ordered on the nanoscale. Our simulations look at: This page contains a brief summary of some of our recent research projects, with links to more details. See also our publications page to find out more about the research work we do.

Simulation of liquid crystal molecules

ODBP-Ph-C7
Chemical structure of the biaxial nematic liquid crystal molecule ODBP=Ph-C7
The first "true" low molecular weight biaxial nematic.






We have recently1 used atomistic simulation to study the structure of the first low molecular weight biaxial liquid crystal, ODBP-Ph-C7. The sequence 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. Read More.

snapshots showing a sequence with time of the growth of ferroelectric biaxial domains



chemical structure of the liquid crystal PCH5

We have also (with our experimental colleagues) been developing techniques to understand dynamics in liquid crystals and to use simulation as a tool to help interpret experimental measurements of dynamical processes.2 Read More.

We have developed new methods for calculating material properties in liquid crystal phases from atomistic simulations. 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.3 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.4

Read More.
chemical structure of the liquid crystal PCH5

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





Simulation studies of amphiphilic molecules in water

side view of amphiphlic polymer at water surface

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.5 The simulations can be used to interpret the results of complementary neutron diffraction studies by calculating neutron reflectivity curves to match with experiment.

The snapshot (left )is from a simulation study of an amphiphilic polymer at a water-air interface (see cartoon below). At higher surface concentrations the PEO chains are forced down into the aqueous subphase to form a polymer brush.

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.

diagram of an amphiphilic polymer at water surface


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

picture of a side chain liquid crystal polymer

We have been developing coarse-grained simulation methods that allow us to understand structure and organisation in a range of molecular materials. We have recently studied liquid crystalline polymers and dendrimers 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.

The figure (left) 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).



Pictures of Mesophases formed by multiblock copolymers composed of rod and coil segments

The pictures (left) show the structure of a series of mesophases formed by a multiblock copolymer as studied by coarse-grained modelling.7 The phases are nematic (top left), lamellar (top right and middle right for different lengths of rods), gyroid (middle left), micellar (bottom). The bottom left structure show elongated micelles formed from polymer rods.

Picture of a polymer nanowire

The second picture shows the formation of a polymer nanowire within a simulation.

Such wires may well have future applications in the area of nanoelectronics. Read More.








Calculations of helical twisting power - chirality mediated by intermolecular forces

When a nematic liquid crystal phase is doped with a small amount of a chiral molecule, the whole phase becomes chiral, and the liquid crystal director (which describes the preferred direction of  alignment for the molecules) is twisted through space to form a helix. Some molecules are better than others at inducing this twist, and the effectiveness of an individual molecule is measured by its Helical Twisting Power (HTP). There is considerable industrial interest in making new molecules with very large HTPs, to use as chiral dopants for chiral polymer films. However, it is difficult to relate chemical structure to HTP.

We developed new methods for predicting HTP prior to synthesis.

Our work suggests that each separate molecular conformation has a different HTP, and we can show how the overall HTP value changes with increases in temperature (as higher energy conformations with different HTPs start to become populated). This mechanism explains the temperature-induced reversal in HTP that occurs with some materials. It also explains the fact that some molecules can have different HTPs in different solvents. The latter is caused by some conformations being preferentially selected in certain solvents.

Recently, we have applied out methods to "banana molecules" and the results help explain the remarkable observation that  bananas can behave as if they are chiral, even though the molecules themselves are achiral.

Read More.


Parallel Simulation Techniques

A key part of our work involves the design and implementation of parallel simulation techniques. These techniques allow computational problems to be tackled by an array of powerful processors working together. We have recently developed new computer codes which allow molecular dynamics and Monte Carlo simulations to be carried out in parallel. These techniques are readily translated to other parallel systems, paving the way for a new generation of large-scale simulations of chemical systems. Our current work in this area involves in the design of new methods for the large-scale simulation of polymer systems and the development of parallel tempering methods. More information about our programs is available from our software page.
Read More.


References

  1. Pelaez J. and Wilson M. R. Atomistic Simulations of a Thermotropic Biaxial Liquid Crystal. Phys. Rev. Lett., 2006, 97., 267801.

  2. 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).

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

  4. 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, art. no. 051709, 1.

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

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

  7. Lintuvuori, J. S., Wilson M. R. A coarse-grained simulation study of mesophase formation in a series of rod-coil multiblock copolymers. , Phys. Chem. Chem. Phys., 2009,
    DOI:10.1039/b818616b.

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

  9. Earl D. J., M., Wilson M. R. Predictions of molecular chirality and helical twisting powers: A theoretical study.  J. Chem. Phys., 2003, 119, 10280.


Link to Full Publication List
Page maintained by Mark Wilson
Document updated: