Source author record

Benjamin J. Morgan

Benjamin J. Morgan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

6works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2026arXiv

Accurate Estimation of Diffusion Coefficients and their Uncertainties from Computer Simulation

Self-diffusion coefficients, $D^*$, are routinely estimated from molecular dynamics simulations by fitting a linear model to the observed mean-squared displacements (MSDs) of mobile species. MSDs derived from simulation exhibit statistical noise that causes uncertainty in the resulting estimate of $D^*$. An optimal scheme for estimating $D^*$ minimises this uncertainty, i.e., it will have high statistical efficiency, and also gives an accurate estimate of the uncertainty itself. We present a scheme for estimating $\D$ from a single simulation trajectory with high statistical efficiency and accurately estimating the uncertainty in the predicted value. The statistical distribution of MSDs observable from a given simulation is modelled as a multivariate normal distribution using an analytical covariance matrix for an equivalent system of freely diffusing particles, which we parameterise from the available simulation data. We use Bayesian regression to sample the distribution of linear models that are compatible with this multivariate normal distribution, to obtain a statistically efficient estimate of $D^*$ and an accurate estimate of the associated statistical uncertainty.

preprint2020arXiv

An introduction to classical molecular dynamics simulation for experimental scattering users

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.

preprint2015arXiv

Density Functional Theory screening of gas-treatment strategies for stabilization of high energy-density lithium metal anodes

To explore the potential of molecular gas treatment of freshly cut lithium foils in non-electrolyte based passivation of high energy-density Li anodes, density functional theory (DFT) has been used to study the decomposition of molecular gases on metallic lithium surfaces. By combining DFT geometry optimization and Molecular Dynamics, the effects of atmospheric (N2, O2, CO2) and hazardous (F2, SO2) gas decomposition on Li(bcc) (100), (110), and (111) surfaces on relative surface energies, work functions, and emerging electronic and elastic properties are investigated. The simulations suggest that exposure to different molecular gases can be used to induce and control reconstructions of the metal Li surface and substantial changes (up to over 1 eV) in the work function of the passivated system. Contrary to the other considered gases, which form metallic adlayers, SO2 treatment emerges as the most effective in creating an insulating passivation layer for dosages <= 1 mono-layer. The substantial Li->adsorbate charge transfer and adlayer relaxation produce marked elastic stiffening of the interface, with the smallest change shown by nitrogen-treated adlayers.

preprint2015arXiv

Sparse Cyclic Excitations Explain the Low Ionic Conductivity of Stoichiometric Li$_{7}$La$_{3}$Zr$_{2}$O$_{12}$

We have performed long time-scale molecular dynamics simulations of the cubic and tetragonal phases of the solid lithium-ion--electrolyte Li$_{7}$La$_{3}$Zr$_{2}$O$_{12}$ (LLZO), using a first-principles parameterised interatomic potential. Collective lithium transport was analysed by identifying dynamical excitations; persistent ion displacements over distances comparable to the separation between lithium sites, and string-like clusters of ions that undergo cooperative motion. We find that dynamical excitations in c-LLZO are frequent, with participating lithium numbers following an exponential distribution, mirroring the dynamics of fragile glasses. In contrast, excitations in t-LLZO are both temporally and spatially sparse, consisting preferentially of highly concerted lithium motion around closed loops. This qualitative difference is explained as a consequence of lithium ordering in t-LLZO, and provides a mechanistic basis for the much lower ionic conductivity of t-LLZO compared to c-LLZO.

preprint2014arXiv

Molecular Dynamics Simulation of Coherent Interfaces in Fluorite Heterostructures

The standard model of enhanced ionic conductivities in solid electrolyte heterostructures follows from a continuum mean-field description of defect distributions that makes no reference to crystalline structure. To examine ionic transport and defect distributions while explicitly accounting for ion-ion correlations and lattice effects, we have performed molecular dynamics simulations of a model coherent fluorite heterostructure without any extrinsic defects, with a difference in standard chemical potentials of mobile fluoride ions between phases induced by an external potential. Increasing the offset in fluoride ion standard chemical potentials across the internal interfaces decreases the activation energies for ionic conductivity and diffusion and strongly enhances fluoride ion mobilities and defect concentrations near the heterostructure interfaces. Non-charge-neutral "space-charge" regions, however, extend only a few atomic spacings from the interface, suggesting a continuum model may be inappropriate. Defect distributions are qualitatively inconsistent with the predictions of the continuum mean-field model, and indicate strong lattice-mediated defect-defect interactions. We identify an atomic-scale "Frenkel polarisation" mechanism for the interfacial enhancement in ionic mobility, where preferentially oriented associated Frenkel pairs form at the interface and promote local ion mobility via concerted diffusion processes.

preprint2014arXiv

Relationships Between Atomic Diffusion Mechanisms and Ensemble Transport Coefficients in Crystalline Polymorphs

Ionic transport in conventional ionic solids is generally considered to proceed via independent diffusion events or "hops''. This assumption leads to well-known Arrhenius expressions for transport coefficients, and is equivalent to assuming diffusion is a Poisson process. Using molecular dynamics simulations of the low-temperature B1, B3, and B4 AgI polymorphs, we have compared rates of ion-hopping with corresponding Poisson distributions to test the assumption of independent hopping in these common structure-types. In all cases diffusion is a non-Poisson process, and hopping is strongly correlated in time. In B1 the diffusion coefficient can be approximated by an Arrhenius expression, though the physical significance of the parameters differs from that commonly assumed. In low temperature B3 and B4 diffusion is characterised by concerted motion of multiple ions in short closed loops. Diffusion coefficients can not be expressed in a simple Arrhenius form dependent on single-ion free-energies, and intrinsic diffusion must be considered a many-body process.