Advanced Monte-Carlo simulation techniques for classical crystalline solids, including confined crystals

Lead Research Organisation: University of Strathclyde
Department Name: Chemical and Process Engineering

Abstract

All chemicals form crystalline solids, so naturally the variety of crystals and crystal properties is immense. Much of the technology of our modern society is based on crystalline solids. For example, inorganic semi-conductor crystals have given us computers, the internet and just about every electronic device there is. Molecular crystal engineering will very likely be just as important in the future. There is currently great interest in polymer crystals for optical and electronic applications (LEDs, lasers, and solar cells for example), in molecular organic crystals for high-technology devices such as photocopiers and computer displays, in gas-hydrates (which are a form of ice with encapsulated gas) which could well be important for providing our future energy needs (it is thought that there is more methane in the worlds natural gas-hydrate fields than in its oil and gas fields) and for locking-away the resulting greenhouse gases that are thought to be responsible for global warming, in advanced porous materials (or special crystal 'frameworks') which could be important for future gas separation processes (allowing carbon dioxide to be separated from exhausts, for example), in biologically related crystals (pharmaceuticals and proteins, for example) that are used in medicines and to understand how the human body works, and so on.All these developments will depend on being able to understand why a particular crystalline solid exists in a particular form, and how its properties depend on factors such as its dimension, or geometry, the shape of its molecules, their interactions with other molecules, and how molecules assemble into larger building blocks.Modern computers and methods now allow scientists to model, or simulate, many important molecular systems quickly, efficiently and at low cost. This is particularly true when modelling extremely small, or nanoscale molecular systems (where the system's dimensions are not that much larger than its molecules), because experiments on this scale are very difficult, and therefore expensive, to setup up and analyse. For example, molecular simulation is used routinely to model fluids, especially those confined in nanoscale pores. However, current methods for simulating crystalline solids, particularly those within these pores, are not entirely satisfactory. At best they are unwieldy and in the case of nanoscale crystals they are often flawed. The aim of the work proposed here is to develop new and better methods for simulating molecular crystals, particularly nanoscale crystals. I expect the molecular simulation methods developed during this work will become the standard method for simulating crystals because they will be the most straightforward, robust and efficient methods available. So these methods could have a direct impact on our ability to simulate molecular crystals correctly and efficiently, and therefore provide a 'step-change' in the prospects for engineering these crystals in the future.

Publications

10 25 50
publication icon
Sweatman M (2009) New techniques for simulating crystals in Molecular Simulation

publication icon
Sweatman MB (2008) The self-referential method combined with thermodynamic integration. in The Journal of chemical physics

 
Description Free energy calculations are important for determining the relative stability of different competing phases, and hence map phase diagrams of model systems of interest. For example, they can be used to understand the phase behaviour materials at high pressure or low temperature. This is useful for a wide range of technologies, e.g. pharmaceuticals. This work concerns a new method for calculating the free energy of crystalline materials. The aim was to further develop a novel 'self-referential' technique so that it can be used quite efficiently and generally for crystals. The first, and key, outcome of this work[1] demonstrated how the self-referential technique could be made much more efficient by using a mathematical procedure called 'thermodynamic integration'. The second and third outcomes demonstrates that the self-referential technique can be applied quite generally to molecular crystals with internal degrees of freedom[2,3], i.e. essentially arbitrary bulk molecular crystals. For example, the method has recently been applied to calculate the free energy of gas hydrates used to store hydrogen[4].

However, it is now expected that the self-referential method is not restricted to molecular crystals, but can be used to calculate the free energy of any equilibrium system, including confined cyrystals, fluids, liquid crystals, etc, and work is ongoing to demonstrate this. If true, this would yield a very valuable tool for the whole molecular simulation community.

[1] The self-referential method combined with thermodynamic integration, Sweatman, M. B. , Atamas, A. A. & Leyssale, J-M. Feb-2008 In : Journal of Chemical Physics. 128, 6, 10.
[2] The self-referential method for linear rigid bodies: Application to hard and Lennard-Jones dumbbells, Sweatman, M. B. , Atamas, A. & Leyssale, J-M. Jan-2009 In : Journal of Chemical Physics. 130, 2, p. 024101.
[3] New techniques for simulating crystals, Sweatman, M. B. 2009 In : Molecular Simulation. 35, 10-11, p. 897.
[4] Monte Carlo calculations of the free energy of ice-like structures using the self-referential method, Atamas, A. , Koudriachova, M. V. , de Leeuw, S. W. & Sweatman, M. B. 2011 In : Molecular Simulation. 37, 4, p. 284.