Mineral crystallisationMinerals are of great significance to everything from health, to the environment and economy. Despite their ubiquitous importance, our understanding of how minerals form remains incomplete. This in turn limits our ability to control and manipulate their formation. Being able to exert such control would have benefits in reducing scale production in industrial processes through to processes such as long term carbon sequestration. Our research focuses on the family of biominerals that are pivotal to both these aforementioned areas, as well as responsible for the minerals that occur in the human body. The objective of our project is to understand how such minerals assemble at the atomic level – something that is hard to achieve in the laboratory, but is becoming possible using computer models.
Principal investigatorJulian Gale J.Gale@curtin.edu.au
Area of scienceChemistry
Applications usedLAMMPS, OpenMM, PSI4, Plumed
Crystals are the basis of much of the solid matter around us. Understanding how they form from their component molecules and/or ions in solution has long been a challenge. If we were able to control this process then we would be able to limit the formation of scale that blocks pipes and vessels that costs industry in lost productivity, as well as being able to design new materials such as pharmaceuticals. In the context of minerals specifically, the formation of crystals is one of the processes by which systems are purified in industry and in nature this goes to form everything from sand on the beach, to teeth. At the moment we have models of the rate of crystal growth based on fitting macroscopic experimental data. However this requires many assumptions about the component steps that are happening. Even worse, there is even less known regarding the initial “birth” of crystals, referred to as nucleation.
One way that we can determine the individual steps take by ions and molecules to form a crystal is through using computer simulation. Here we can model what happens when atoms are placed in a given arrangement and how this changes with time through the use of a technique known as molecular dynamics. Unfortunately the amount of real time that can observed in this way is limited to typically less than a millionth of a second, which is slower than most experiments, especially geological processes. However, by applying a bias to the simulations that allows them to jump over large barriers we can dramatically accelerate the evolution of the model. Although simple knowledge of time is lost, we gain a wealth of information regarding the thermodynamics, which is the driving force for chemical and physical reactions. Nevertheless, the computational demands remain high as billions of steps are still required – something that is only possible with state of the art supercomputers
Although we know how to overcome the problem of slow processes in our simulation model, the process still has to be solved in discrete time steps that are of the order of 10^-15 of a second due to the fundamental nature of chemical interactions. Therefore each model requires up to a billion iterations to acquire sufficient data to be useful. Each iteration can also be expensive when using the most accurate methods available. Only through the use of many thousands of processors working simultaneously on the same model is this feasible. Access to world-class supercomputers, such as the Magnus system with tens of thousands of cores, is therefore critical to our research. In future, new technologies, such as graphical processing units, may also allow us to speed up the evaluation of each step – something that we have also been exploring using the resources of the Pawsey Centre.
The CrystalGrower software is being developed by the University of Manchester as a spin out company and we are working with this team to advance the method. We have developed an automated protocol that allows the initial parameters required to run a CrystalGrower simulation to be run for molecular crystalline systems and this has been embedded in a version of the CrystalExplorer code. This creates the potential for the commercial use of our method via this spin out company.
We have also developed a periodic implementation of the GFN-FF universal force field. While the parameter generation method will be made available under LGPL due to the dependence on existing code, the use of the parameters will form part of the GULP code which is distributed via Biovia (Dassault Systems) to clients.
NB: Development also linked to pawsey0184, f97
Commercial Advantage of this Project
Most of the developments are not yet at the commercial stage yet, though there is the potential for industry to benefit from these developments. At this stage we have formed a partnership with BASF along with the University of Manchester to explore the development of the CrystalGrower technology for industry application. There is the potential for further commercial use in the pharmaceutical sector in particular. NB: Development also linked to pawsey0184, f97
List of Publications
1) J.D. Gale, L. LeBlanc, P.R. Spackman, A. Silvestri and P. Raiteri, “A universal force field for materials, periodic GFN-FF: Implementation and examination”, J. Chem. Theory Comput., 17, 7827-7849 (2021)
2) A. Schuitemaker, J. Aufort, K.B. Koziara, R. Demichelis, P. Raiteri and J.D. Gale, “Simulating the bind-ing of key organic functional groups to aqueous calcium carbonate species”, PhysChemChemPhys, 23, 27253-27265 (2021)
3) P. Kraus, “Extrapolating DFT toward the complete basis set limit: Lessons from the PBE family of functionals”, J. Chem. Theory Comput., 17, 5651-5660 (2021)
4) J. Aufort, A. Schuitemaker, R. Green, R. Demichelis, P. Raiteri and J.D. Gale, “Determining the adsorp-tion free energies of small organics and intrinsic ions at the terrace and steps of calcite”, Crystal Growth and Design, (accepted December 29th 2021)