Molecular dynamics calculations of oxide/water interfaces


Computer simulations hold tremendous promise to accelerate the molecular engineering of green energy technologies, such as new systems for electrical energy storage and solar energy usage, as well as carbon dioxide capture from the environment. However, the predictive power of these simulations depends on having a means to confirm that they do indeed describe the real world.

Such confirmation is no simple task. Many assumptions enter the setup of these simulations. As a result, the simulations must be carefully checked by using an appropriate validation protocol involving experimental measurements.

To address this challenge, a team of scientists at the U.S. Department of Energy’s (DOE) Argonne National Laboratory, the University of Chicago and the University of California, Davis, developed a groundbreaking validation protocol for simulations of the atomic structure of the interface between a solid (a metal oxide) and liquid water. The team was led by Giulia Galli, a theorist with a joint appointment at Argonne and the University of Chicago, and Paul Fenter, an Argonne experimentalist.

The authors focused on a solid/liquid interface because interfaces are ubiquitous in materials, and those between oxides and water are key in many energy applications. To date, most validation protocols have been designed for bulk materials, ignoring interfaces. The authors believed that the atomic-scale structure of surfaces and interfaces in realistic environments would present a particularly sensitive, and therefore challenging, validation approach.

The validation procedure they designed uses high-resolution X-ray reflectivity (XR) measurements as the experimental pillar of the protocol. The team compared XR measurements for an aluminum oxide/water interface, conducted at beamline 33-ID-D at Argonne’s Advanced Photon Source (APS), with results obtained by running high-performance computer simulations at the Argonne Leadership Computing Facility (ALCF).  

Their measurements detected the reflection of very high energy X-ray beams from an oxide/water interface. At the beam energies generated at the APS, the X-ray wavelengths are similar to interatomic distances. This allows the researchers to directly probe the molecular-scale structure of the interface.

XR is an ideal probe to obtain experimental results directly comparable to simulations. The team ran the simulations at the ALCF using the Qbox code, which is designed to study finite temperature properties of materials and molecules using simulations based on quantum mechanics. They were able to test several approximations of the theory. The team compared measured XR intensities with those calculated from several simulated structures. They also investigated how X-rays scattered from the electrons in different parts of the sample would interfere to produce the experimentally observed signal.

The endeavor of the team turned out to be more challenging than anticipated. There was a bit of a trial and error at the beginning when we were trying to understand the right geometry to adopt and the right theory that would give us accurate results. However, their back and forth between theory and experiment paid off, and they were able to set up a robust validation protocol that can now be deployed for other interfaces as well. Indeed, the validation protocol helped quantify the strengths and weaknesses of the simulations, providing a pathway toward building more accurate models of solid/liquid interfaces in the future. The simulations also shed new insight on the XR measurements themselves. In particular, they showed that the data are sensitive not only to the atomic positions, but also to the electron distribution surrounding each atom in subtle and complex ways. These insights will prove beneficial to future experiments on oxide/liquid interfaces.

molecular dynamics calculations of oxide/water interfaces - Advances in Engineering
FIGURE: Pictorial representation of joint experimental and computational study of materials. The study utilized the Advanced Photon Source (upper panel) and Argonne Leadership Computing Facility (lower panel). The team addressed the atomistic structure of interfaces, which are ubiquitous in materials. Credit: Emmanuel Gygi, University of California, San Diego

About the author

Giulia Galli

Liew Family Professor of Molecular Engineering
Professor of Chemistry

Giulia Galli is the Liew Family professor of Electronic Structure and Simulations in the Pritzker School of Molecular Engineering and Professor of Chemistry at the University of Chicago. She also holds a Senior Scientist position at Argonne National Laboratory (ANL) and she is the director of the Midwest Integrated Center for Computational Materials (MICCoM). Prior to joining UChicago and ANL, she was Professor of Chemistry and Physics at UC Davis (2005-2013) and the head of the Quantum Simulations group at the Lawrence Livermore National Laboratory (1998-2005). She holds a Ph.D. in Physics from the International School of Advanced Studies (SISSA) in Trieste, Italy.

She is a Fellow of the American Physical Society, the American Association for the Advancement of Science and the American Academy of Arts of Sciences, as well as a member of the National Academy of Sciences. She is a recipient of numerous awards, including the Lawrence Livermore National Laboratory Science and Technology Award, the US Department of Energy Award of Excellence, the Materials Research Society Theory Award, the APS David Adler Lectureship in Materials Physics, the Feynman Nanotechnology Prize in Theory, the medal of the Schola Physica Romana and the Tomassoni-Chisesi award by the Sapienza University of Rome.

Her research activity is focused on the development and use of theoretical and computational methods to understand and predict the properties and behavior of materials (solids, liquids, and nanostructures) from first principles.


Katherine J. Harmon, Kendra Letchworth-Weaver, Alex P. Gaiduk, Federico Giberti, Francois Gygi, Maria K. Y. Chan, Paul Fenter, and Giulia Galli. Validating first-principles molecular dynamics calculations of oxide/water interfaces with x-ray reflectivity data Phys. Rev. Materials 4, 113805 – 2020

Go To Phys. Rev. Materials

Check Also

A facile mechanochemical technique to visualize polymer chain breakage in hydrogel - Advances in Engineering

A facile mechanochemical technique to visualize polymer chain breakage in hydrogel