-
Grain boundary segregation and phase separation in ceria-zirconia from atomistic simulation
Authors:
Tom L. Underwood,
Susanna Vigorito,
Marco Molinari,
John Purton,
Nigel B. Wilding,
John T. S. Irvine,
Stephen C. Parker
Abstract:
Doping is the most common strategy employed in the development of new and improved materials. However, predicting the effects of doping on the atomic-scale structure of a material is often difficult or limited to high-end experimental techniques. Doping can induce phase separation in a material, undermining the material's stability. A further complication is that dopant atoms can segregate to inte…
▽ More
Doping is the most common strategy employed in the development of new and improved materials. However, predicting the effects of doping on the atomic-scale structure of a material is often difficult or limited to high-end experimental techniques. Doping can induce phase separation in a material, undermining the material's stability. A further complication is that dopant atoms can segregate to interfaces in a material such as grain boundaries (GBs), with consequences for key macroscopic properties of the material such as its conductivity. Here, we describe a computational methodology based on semi-grand canonical Monte Carlo which can be used to probe these phenomena at the atomic scale for metal oxide solid solutions. The methodology can provide precise predictions of the thermodynamic conditions at which phase separation occurs. It can also provide the segregation patterns exhibited by GBs at given conditions. We apply the methodology to one of the most important catalytic materials, ceria-zirconia. Our calculations reveal an interesting richness in the GB segregation in this system. Most GBs we examined exhibited continuous increases in Zr segregation upon Zr doping, with a concomitant reduction in the formation enthalpies of the GBs. However, a few GBs exhibited no segregation at low temperatures. We also observed evidence of first-order complexion transitions in some GBs.
△ Less
Submitted 7 October, 2023; v1 submitted 29 September, 2023;
originally announced September 2023.
-
dlmontepython: A Python library for automation and analysis of Monte Carlo molecular simulations
Authors:
T. L. Underwood,
J. A. Purton,
J. R. H. Manning,
A. V. Brukhno,
K. Stratford,
T. Düren,
N. B. Wilding,
S. C. Parker
Abstract:
We present an open source Python 3 library aimed at practitioners of molecular simulation, especially Monte Carlo simulation. The aims of the library are to facilitate the generation of simulation data for a wide range of problems; and to support data analysis methods which enable one to make the most of previously generated data. The library contains a framework for automating the task of measuri…
▽ More
We present an open source Python 3 library aimed at practitioners of molecular simulation, especially Monte Carlo simulation. The aims of the library are to facilitate the generation of simulation data for a wide range of problems; and to support data analysis methods which enable one to make the most of previously generated data. The library contains a framework for automating the task of measuring target physical properties (e.g. density) over a range of thermodynamic parameters (e.g. temperature) calculated using a molecular simulation program, in particular the Monte Carlo program DL_MONTE. The library also supports analysis methods including block averaging, equilibration detection and histogram reweighting. Here we describe the library and provide examples to demonstrate its key functionality: we use the library to automatically calculate isotherms to a specified precision; and to calculate the surface tension and liquid-vapour coexistence properties of methane.
△ Less
Submitted 8 April, 2021;
originally announced April 2021.
-
Accuracy of hybrid functionals with non-self-consistent Kohn-Sham orbitals for predicting the properties of semiconductors
Authors:
Jonathan M. Skelton,
David S. D. Gunn,
Sebastian Metz,
Stephen C. Parker
Abstract:
Accurately modeling the electronic structure of materials is a persistent challenge to high-throughput screening. A promising means of balancing accuracy against computational cost are non-self-consistent calculations with hybrid density-functional theory, where the electronic band energies are evaluated using a hybrid functional from orbitals obtained with a less demanding (semi-)local functional…
▽ More
Accurately modeling the electronic structure of materials is a persistent challenge to high-throughput screening. A promising means of balancing accuracy against computational cost are non-self-consistent calculations with hybrid density-functional theory, where the electronic band energies are evaluated using a hybrid functional from orbitals obtained with a less demanding (semi-)local functional. We have quantified the performance of this technique for predicting the physical properties of sixteen tetrahedral semiconductors with bandgaps from 0.2-5.5 eV. Provided the base functional predicts a non-metallic electronic structure, bandgaps within 5 % of the PBE0 and HSE06 gaps can be obtained with an order of magnitude reduction in computing time. The positions of the valence and conduction band extrema and the Fermi level are well reproduced, further enabling calculation of the band dispersion, density of states, and dielectric properties using Fermi's Golden Rule. While the error in the non-self-consistent total energies is ~50 meV atom-1, the energy-volume curves are reproduced accurately enough to obtain the equilibrium volume and bulk modulus with minimal error. We also test the dielectric-dependent scPBE0 functional and obtain bandgaps and dielectric constants to within 2.5 % of the self-consistent results, which amount to a significant improvement over self-consistent PBE0 for a similar computational cost. We identify cases where the non-self-consistent approach is expected to perform poorly, and demonstrate that partial self-consistency provides a practical and efficient workaround. Finally, we perform proof-of-concept calculations on CoO and NiO to demonstrate the applicability of the technique to strongly-correlated open-shell transition-metal oxides.
△ Less
Submitted 1 May, 2020;
originally announced May 2020.
-
An introduction to classical molecular dynamics simulation for experimental scattering users
Authors:
Andrew R. McCluskey,
James Grant,
Adam R. Symington,
Tim Snow,
James Doutch,
Benjamin J. Morgan,
Stephen C. Parker,
Karen J. Edler
Abstract:
Classical molecular dynamics simulations are a common component of multi-modal analyses from scattering measurements, such as small-angle scattering and diffraction. Users of these experimental techniques often have no formal training in the theory and practice of molecular dynamics simulation, leading to the possibility of these simulations being treated as a "black box" analysis technique. In th…
▽ More
Classical molecular dynamics simulations are a common component of multi-modal analyses from scattering measurements, such as small-angle scattering and diffraction. Users of these experimental techniques often have no formal training in the theory and practice of molecular dynamics simulation, leading to the possibility of these simulations being treated as a "black box" analysis technique. In this article, we describe an open educational resource (OER) designed to introduce classical molecular dynamics to users of scattering methods. This resource is available as a series of interactive web pages, which can be easily accessed by students, and as an open source software repository, which can be freely copied, modified, and redistributed by educators. The topic covered in this OER includes classical atomistic modelling, parameterising interatomic potentials, molecular dynamics simulations, typical sources of error, and some of the approaches to using simulations in the analysis of scattering data.
△ Less
Submitted 11 May, 2020; v1 submitted 4 February, 2019;
originally announced February 2019.
-
Assessing molecular simulation for the analysis of lipid monolayer reflectometry
Authors:
Andrew R. McCluskey,
James Grant,
Andrew J. Smith,
Jonathan L. Rawle,
David J. Barlow,
M. Jayne Lawrence,
Stephen C. Parker,
Karen J. Edler
Abstract:
Using molecular simulation to aid in the analysis of neutron reflectometry measurements is commonplace. However, reflectometry is a tool to probe large-scale structures, and therefore the use of all-atom simulation may be irrelevant. This work presents the first direct comparison between the reflectometry profiles obtained from different all-atom and coarse-grained molecular dynamics simulations.…
▽ More
Using molecular simulation to aid in the analysis of neutron reflectometry measurements is commonplace. However, reflectometry is a tool to probe large-scale structures, and therefore the use of all-atom simulation may be irrelevant. This work presents the first direct comparison between the reflectometry profiles obtained from different all-atom and coarse-grained molecular dynamics simulations. These are compared with a traditional model layer structure analysis method to determine the minimum simulation resolution required to accurately reproduce experimental data. We find that systematic limits reduce the efficacy of the MARTINI potential model, while the Berger united-atom and Slipids all-atom potential models agree similarly well with the experimental data. The model layer structure gives the best agreement, however, the higher resolution simulation-dependent methods produce an agreement that is comparable. Finally, we use the atomistic simulation to advise on possible improvements that may be offered to the model layer structures, creating a more realistic monolayer model.
△ Less
Submitted 11 May, 2020; v1 submitted 16 January, 2019;
originally announced January 2019.
-
Bayesian determination of the effect of a deep eutectic solvent on the structure of lipid monolayers
Authors:
Andrew R. McCluskey,
Adrian Sanchez-Fernandez,
Karen J. Edler,
Stephen C. Parker,
Andrew J. Jackson,
Richard A. Campbell,
Tom Arnold
Abstract:
In this work, we present the first example of the self-assembly of phospholipid monolayers at the interface between air and an ionic solvent. Deep eutectic solvents are a novel class of environmentally friendly non-aqueous room temperature liquids with tunable properties, that have wide ranging potential applications and are capable of promoting the self-assembly of surfactant molecules. We use a…
▽ More
In this work, we present the first example of the self-assembly of phospholipid monolayers at the interface between air and an ionic solvent. Deep eutectic solvents are a novel class of environmentally friendly non-aqueous room temperature liquids with tunable properties, that have wide ranging potential applications and are capable of promoting the self-assembly of surfactant molecules. We use a chemically-consistent Bayesian modelling of X-ray and neutron reflectometry measurements to show that these monolayers broadly behave as they do on water. This method allows for the monolayer structure to be determined, alongside the molecular volumes of the individual monolayer components without the need for water-specific constraints to be introduced. Furthermore, using this method we are able to better understand the correlations present between parameters in the analytical model. This example of a non-aqueous phospholipid monolayer has important implications for the potential uses of these solvents and for our understanding of how biomolecules behave in the absence of water.
△ Less
Submitted 11 May, 2020; v1 submitted 17 October, 2018;
originally announced October 2018.
-
Electronic and Phonon Instabilities in Bilayer Graphene under Applied External Bias
Authors:
E. Lora da Silva,
M. C. Santos,
J. M. Skelton,
Tao Yang,
T. Santos,
S. C. Parker,
A. Walsh
Abstract:
We have performed electronic-structure and lattice-dynamics calculations on the AB and AA structures of bilayer graphene. We study the effect of external electric fields and compare results obtained with different levels of theory to existing theoretical and experimental results. Application of an external field to the AB bilayer alters the electronic spectrum, with the bands changing under bias f…
▽ More
We have performed electronic-structure and lattice-dynamics calculations on the AB and AA structures of bilayer graphene. We study the effect of external electric fields and compare results obtained with different levels of theory to existing theoretical and experimental results. Application of an external field to the AB bilayer alters the electronic spectrum, with the bands changing under bias from a parabolic to a "Mexican hat" double-well structure. This results in a semi-metal-to-semiconductor phase transition, with the size of the induced electronic band-gap being tuneable through the field strength. A reduction of continuous symmetry from a hexagonal to a triangular lattice is also evidenced through in-plane electronic charge inhomogeneities between the sublattices. When spin-orbit coupling is turned on for the AB system, we find that the bulk gap decreases, gradually increasing for larger intensities of the bias. Under large bias the energy dispersion recovers the Mexican hat structure, since the energy interaction between the layers balances the coupling interaction. We find that external bias perturbs the harmonic phonon spectra and leads to anomalous behaviour of the out-of-plane flexural ZA and layer-breathing ZO' modes. For the AA system, the electronic and phonon dispersions both remain stable under bias, but the phonon spectrum exhibits zone-center imaginary modes due to layer-sliding dynamical instabilities.
△ Less
Submitted 18 September, 2018;
originally announced September 2018.
-
Particle Morphology and Lithium Segregation to Surfaces of the Li$_7$La$_3$Zr$_2$O$_{12}$ Solid Electrolyte
Authors:
Pieremanuele Canepa,
James A. Dawson,
Gopalakrishnan Sai Gautam,
Joel M. Statham,
Stephen C. Parker,
M. Saiful Islam
Abstract:
Solid electrolytes for solid-state Li-ion batteries are stimulating considerable interest for next-generation energy storage applications. The Li$_7$La$_3$Zr$_2$O$_{12}$ garnet-type solid electrolyte has received appreciable attention as a result of its high ionic conductivity. However, several challenges for the successful application of solid-state devices based on Li$_7$La$_3$Zr$_2$O$_{12}$ rem…
▽ More
Solid electrolytes for solid-state Li-ion batteries are stimulating considerable interest for next-generation energy storage applications. The Li$_7$La$_3$Zr$_2$O$_{12}$ garnet-type solid electrolyte has received appreciable attention as a result of its high ionic conductivity. However, several challenges for the successful application of solid-state devices based on Li$_7$La$_3$Zr$_2$O$_{12}$ remain, such as dendrite formation and maintaining physical contact at interfaces over many Li intercalation/extraction cycles. Here, we apply first-principles density functional theory to provide insights into the Li$_7$La$_3$Zr$_2$O$_{12}$ particle morphology under various physical and chemical conditions. Our findings indicate Li segregation at the surfaces, suggesting Li-rich grain boundaries at typical synthesis and sintering conditions. On the basis of our results, we propose practical strategies to curb Li segregation at the Li$_7$La$_3$Zr$_2$O$_{12}$ interfaces. This approach can be extended to other Li-ion conductors for the design of practical energy storage devices.
△ Less
Submitted 14 April, 2018;
originally announced April 2018.
-
Anharmonicity in the high-temperature Cmcm phase of SnSe: soft modes and three-phonon interactions
Authors:
Jonathan M. Skelton,
Lee A. Burton,
Stephen C. Parker,
Aron Walsh,
Chang-Eun Kim,
Aloysius Soon,
John Buckeridge,
Alexey A. Sokol,
C. Richard A. Catlow,
Atsushi Togo,
Isao Tanaka
Abstract:
The layered semiconductor SnSe is one of the highest-performing thermoelectric materials known. We demonstrate, through a first-principles lattice-dynamics study, that the high-temperature Cmcm phase is a dynamic average over lower-symmetry minima separated by very small energetic barriers. Compared to the low-temperature Pnma phase, the Cmcm phase displays a phonon softening and enhanced three-ph…
▽ More
The layered semiconductor SnSe is one of the highest-performing thermoelectric materials known. We demonstrate, through a first-principles lattice-dynamics study, that the high-temperature Cmcm phase is a dynamic average over lower-symmetry minima separated by very small energetic barriers. Compared to the low-temperature Pnma phase, the Cmcm phase displays a phonon softening and enhanced three-phonon scattering, leading to an anharmonic damping of the low-frequency modes and hence the thermal transport. We develop a renormalisation scheme to quantify the effect of the soft modes on the calculated properties, and confirm that the anharmonicity is an inherent feature of the Cmcm phase. These results suggest a design concept for thermal insulators and thermoelectric materials, based on displacive instabilities, and highlight the power of lattice-dynamics calculations for materials characterization.
△ Less
Submitted 30 July, 2016; v1 submitted 11 February, 2016;
originally announced February 2016.
-
Thermal physics of the lead chalcogenides PbS, PbSe, and PbTe from first principles
Authors:
Jonathan Michael Skelton,
Stephen C. Parker,
Atsushi Togo,
Isao Tanaka,
Aron Walsh
Abstract:
The lead chalcogenides represent an important family of functional materials, in particular due to the benchmark high-temperature thermoelectric performance of PbTe. A number of recent investigations, experimental and theoretical, have aimed to gather insight into their unique lattice dynamics and electronic structure. However, the majority of first-principles modelling has been performed at fixed…
▽ More
The lead chalcogenides represent an important family of functional materials, in particular due to the benchmark high-temperature thermoelectric performance of PbTe. A number of recent investigations, experimental and theoretical, have aimed to gather insight into their unique lattice dynamics and electronic structure. However, the majority of first-principles modelling has been performed at fixed temperatures, and there has been no comprehensive and systematic computational study of the effect of temperature on the material properties. We report a comparative lattice-dynamics study of the temperature dependence of the properties of PbS, PbSe and PbTe, focussing particularly on those relevant to thermoelectric performance, viz. phonon frequencies, lattice thermal conductivity, and electronic band structure. Calculations are performed within the quasi-harmonic approximation, with the inclusion of phonon-phonon interactions from many-body perturbation theory, which are used to compute phonon lifetimes and predict the lattice thermal conductivity. The results are critically compared against experimental data and other calculations, and add new insight to on-going research on the PbX compounds in relation to the off-centring of Pb at high temperatures, which is shown to be related to phonon softening. The agreement with experiment suggests that this method could serve as a straightforward, powerful and generally-applicable means of investigating the temperature dependence of material properties from first principles.
△ Less
Submitted 24 May, 2014;
originally announced May 2014.