-
A fast anharmonic free energy method with an application to vacancies in ZrC
Authors:
Thomas A. Mellan,
Andrew I. Duff,
Blazej Grabowski,
Michael W. Finnis
Abstract:
We propose an approach to calculate the anharmonic part of the volumetric-strain and temperature dependent free energy of a crystal. The method strikes an effective balance between accuracy and computational efficiency, showing a $\times10$ speed-up on comparable free energy approaches at the level of density functional theory, with average errors less than 1 meV/atom. As a demonstration we make n…
▽ More
We propose an approach to calculate the anharmonic part of the volumetric-strain and temperature dependent free energy of a crystal. The method strikes an effective balance between accuracy and computational efficiency, showing a $\times10$ speed-up on comparable free energy approaches at the level of density functional theory, with average errors less than 1 meV/atom. As a demonstration we make new predictions on the thermodynamics of substoichiometric ZrC$_x$, including vacancy concentration and heat capacity.
△ Less
Submitted 4 April, 2019;
originally announced April 2019.
-
Spontaneous Frenkel pair formation in Zirconium Carbide
Authors:
Thomas A. Mellan,
Andrew I. Duff,
Michael W. Finnis
Abstract:
With density functional theory we have performed molecular dynamics simulations of ZrC which displayed spontaneous Frenkel pair formation at a temperature of 3200 K, some 500 K below the melting point. To understand this behaviour, rarely seen in equilibrium simulations, we quenched and examined a set of lattices containing a Frenkel pair. Five metastable structures were found, and their formation…
▽ More
With density functional theory we have performed molecular dynamics simulations of ZrC which displayed spontaneous Frenkel pair formation at a temperature of 3200 K, some 500 K below the melting point. To understand this behaviour, rarely seen in equilibrium simulations, we quenched and examined a set of lattices containing a Frenkel pair. Five metastable structures were found, and their formation energies and electronic properties were studied. Their thermal generation was found to be facilitated by a reduction of between 0.7 and 1.5 eV in formation energy due to thermal expansion of the lattice. With input from a quasi-harmonic description of the defect free energy of formation, an ideal solution model was used to estimate lower bounds on their concentration as a function of temperature and stoichiometry. At 3000 K (0.81 of the melting temperature) their concentration was estimated to be 1.2% per mole in a stoichiometric crystal, and 0.3% per mole in a crystal with 10% per mole of constitutional vacancies. Their contribution to heat capacity, thermal expansion and bulk modulus was estimated.
△ Less
Submitted 10 September, 2018;
originally announced September 2018.
-
Anomalous Diffusion along Metal/Ceramic Interfaces
Authors:
Aakash Kumar,
Hagit Barda,
Leonid Klinger,
Michael W. Finnis,
Vincenzo Lordi,
Eugen Rabkin,
David J. Srolovitz
Abstract:
Hole formation in a polycrystalline Ni film on an $α$-Al$_2$O$_3$ substrate coupled with a continuum diffusion analysis demonstrates that Ni diffusion along the Ni/$α$-Al$_2$O$_3$ interface is surprisingly fast. Ab initio calculations demonstrate that both Ni vacancy formation and migration energies at the coherent Ni/$α$-Al$_2$O$_3$ interface are much smaller than in bulk Ni, suggesting that the…
▽ More
Hole formation in a polycrystalline Ni film on an $α$-Al$_2$O$_3$ substrate coupled with a continuum diffusion analysis demonstrates that Ni diffusion along the Ni/$α$-Al$_2$O$_3$ interface is surprisingly fast. Ab initio calculations demonstrate that both Ni vacancy formation and migration energies at the coherent Ni/$α$-Al$_2$O$_3$ interface are much smaller than in bulk Ni, suggesting that the activation energy for diffusion along coherent Ni/$α$-Al$_2$O$_3$ interfaces is comparable to that along (incoherent/high angle) grain boundaries. Based on these results, we develop a simple model for diffusion along metal/ceramic interfaces, apply to a wide range of metal/ceramic systems and validate it with several ab initio calculations. These results suggest that fast metal diffusion along metal/ceramic interfaces should be common, but is not universal.
△ Less
Submitted 12 July, 2018;
originally announced July 2018.
-
A model for time-dependent grain boundary diffusion of ions and electrons through a film or scale, with an application to alumina
Authors:
M. P. Tautschnig,
N. M. Harrison,
M. W. Finnis
Abstract:
A model for ionic and electronic grain boundary transport through thin films, scales or membranes with columnar grain structure is introduced. The grain structure is idealized as a lattice of identical hexagonal cells - a honeycomb pattern. Reactions with the environment constitute the boundary conditions and drive the transport between the surfaces. Time-dependent simulations solving the Poisson…
▽ More
A model for ionic and electronic grain boundary transport through thin films, scales or membranes with columnar grain structure is introduced. The grain structure is idealized as a lattice of identical hexagonal cells - a honeycomb pattern. Reactions with the environment constitute the boundary conditions and drive the transport between the surfaces. Time-dependent simulations solving the Poisson equation self-consistently with the Nernst-Planck flux equations for the mobile species are performed. In the resulting Poisson-Nernst-Planck system of equations, the electrostatic potential is obtained from the Poisson equation in its integral form by summation. The model is used to interpret alumina membrane oxygen permeation experiments, in which different oxygen gas pressures are applied at opposite membrane surfaces and the resulting flux of oxygen molecules through the membrane is measured. Simulation results involving four mobile species, charged aluminum and oxygen vacancies, electrons, and holes, provide a complete description of the measurements and insight into the microscopic processes underpinning the oxygen permeation of the membrane. Most notably, the hypothesized transition between p-type and n-type ionic conductivity of the alumina grain boundaries as a function of the applied oxygen gas pressure is observed in the simulations. The range of validity of a simple analytic model for the oxygen permeation rate, similar to the Wagner theory of metal oxidation, is quantified by comparison to the numeric simulations. The three-dimensional model we develop here is readily adaptable to problems such as transport in a solid state electrode, or corrosion scale growth.
△ Less
Submitted 3 May, 2017; v1 submitted 4 February, 2017;
originally announced February 2017.
-
New Methods for Calculating the Free Energy of Charged Defects in Solid Electrolytes
Authors:
Robert M. Horton,
Andrew J. Haslam,
Amparo Galindo,
George Jackson,
Michael W. Finnis
Abstract:
A methodology for calculating the contribution of charged defects to the configurational free energy of an ionic crystal is introduced. The temperature-independent Wang-Landau Monte Carlo technique is applied to a simple model of a solid electrolyte, consisting of charged positive and negative defects on a lattice. The electrostatic energy is computed on lattices with periodic boundary conditions,…
▽ More
A methodology for calculating the contribution of charged defects to the configurational free energy of an ionic crystal is introduced. The temperature-independent Wang-Landau Monte Carlo technique is applied to a simple model of a solid electrolyte, consisting of charged positive and negative defects on a lattice. The electrostatic energy is computed on lattices with periodic boundary conditions, and used to calculate the density of states and statistical-thermodynamic potentials of this system. The free energy as a function of defect concentration and temperature is accurately described by a regular solution model up to concentrations of 10% of defects, well beyond the range described by the ideal solution theory.
The approach, supplemented by short-ranged terms in the energy, is proposed as an alternative to free-energy methods that require a number of simulations to be carried out over a range of temperatures.
△ Less
Submitted 23 May, 2013;
originally announced May 2013.
-
Oxygen vacancy formation energies in Sr-doped complex perovskites: ab initio thermodynamic study
Authors:
Denis Gryaznov,
Mike W. Finnis,
Robert Evarestov,
Joachim Maier
Abstract:
La1-xSrxCo0.25Fe0.75O3-delta is known as one of the best cathode materials for permeation membranes and solid oxide fuel cells. Optimization of its chemical composition is a challenging problem. One of the key properties is concentration of oxygen vacancies, controlled by their formation energies. Ab initio calculations were employed in order to study the formation of oxygen vacancies in La1-xSrxC…
▽ More
La1-xSrxCo0.25Fe0.75O3-delta is known as one of the best cathode materials for permeation membranes and solid oxide fuel cells. Optimization of its chemical composition is a challenging problem. One of the key properties is concentration of oxygen vacancies, controlled by their formation energies. Ab initio calculations were employed in order to study the formation of oxygen vacancies in La1-xSrxCo0.25Fe0.75O3-delta perovskites by varying the Sr content from x = 12.5% to 50%. The formation energies were obtained for different stoichiometries as functions of temperature and oxygen partial pressure. We have shown that the phonon contribution to the free formation energy becomes increasingly important in La1-xSrxCo0.25Fe0.75O3-delta not only with rising temperature but also with rising Sr content. We have shown that the formation energies are decreased significantly with increasing Sr content due to two effects: charge compensation of Sr2+ ions and phonon contribution. We have suggested a simple explanation to increasing role of phonons for the oxygen vacancies formation energies on the basis of phonon mode changes in comparison to defect-free materials. A careful analysis of the experimental results from the literature is also presented.
△ Less
Submitted 27 June, 2012;
originally announced June 2012.
-
Solid-liquid interface free energy through metadynamics simulations
Authors:
Stefano Angioletti-Uberti,
Michele Ceriotti,
Mike W. Finnis,
Peter D. Lee
Abstract:
The solid-liquid interface free energy γsl is a key parameter controlling nucleation and growth during solidification and other phenomena. There are intrinsic difficulties in obtaining accurate experimental values, and the previous approaches to compute γsl with atomistic simulations are computationally demanding. We propose a new approach, which is to obtain γsl from a free energy map of the phas…
▽ More
The solid-liquid interface free energy γsl is a key parameter controlling nucleation and growth during solidification and other phenomena. There are intrinsic difficulties in obtaining accurate experimental values, and the previous approaches to compute γsl with atomistic simulations are computationally demanding. We propose a new approach, which is to obtain γsl from a free energy map of the phase transition reconstructed by metadynamics. We apply this to the benchmark case of a Lennard-Jones potential and the results confirm the most reliable data obtained previously. We demonstrate several advantages of our new approach: it is simple to implement, robust and free of hysteresis problems, it allows a rigorous and unbiased estimate of the statistical uncertainty and it returns a good estimate of of the thermodynamic limit with system sizes of a just a few hundred atoms. It is therefore attractive for using with more realistic and specific models of interatomic forces.
△ Less
Submitted 19 May, 2010; v1 submitted 13 November, 2009;
originally announced November 2009.
-
Magnetic tight-binding and the iron-chromium enthalpy anomaly
Authors:
A. T. Paxton,
M. W. Finnis
Abstract:
We describe a self consistent magnetic tight-binding theory based in an expansion of the Hohenberg-Kohn density functional to second order, about a non spin polarised reference density. We show how a first order expansion about a density having a trial input magnetic moment leads to the Stoner--Slater rigid band model. We employ a simple set of tight-binding parameters that accurately describes…
▽ More
We describe a self consistent magnetic tight-binding theory based in an expansion of the Hohenberg-Kohn density functional to second order, about a non spin polarised reference density. We show how a first order expansion about a density having a trial input magnetic moment leads to the Stoner--Slater rigid band model. We employ a simple set of tight-binding parameters that accurately describes electronic structure and energetics, and show these to be transferable between first row transition metals and their alloys. We make a number of calculations of the electronic structure of dilute Cr impurities in Fe which we compare with results using the local spin density approximation. The rigid band model provides a powerful means for interpreting complex magnetic configurations in alloys; using this approach we are able to advance a simple and readily understood explanation for the observed anomaly in the enthalpy of mixing.
△ Less
Submitted 24 October, 2007;
originally announced October 2007.
-
Structural and chemical embrittlement of grain boundaries by impurities: a general theory and first principles calculations for copper
Authors:
A. Y. Lozovoi,
A. T. Paxton,
M. W. Finnis
Abstract:
First principles calculations of the Sigma 5 (310)[001] symmetric tilt grain boundary in Cu with Bi, Na, and Ag substitutional impurities provide evidence that in the phenomenon of Bi embrittlement of Cu grain boundaries electronic effects do not play a major role; on the contrary, the embrittlement is mostly a structural or "size" effect. Na is predicted to be nearly as good an embrittler as Bi…
▽ More
First principles calculations of the Sigma 5 (310)[001] symmetric tilt grain boundary in Cu with Bi, Na, and Ag substitutional impurities provide evidence that in the phenomenon of Bi embrittlement of Cu grain boundaries electronic effects do not play a major role; on the contrary, the embrittlement is mostly a structural or "size" effect. Na is predicted to be nearly as good an embrittler as Bi, whereas Ag does not embrittle the boundary in agreement with experiment. While we reject the prevailing view that "electronic" effects (i.e., charge transfer) are responsible for embrittlement, we do not exclude the role of chemistry. However numerical results show a striking equivalence between the alkali metal Na and the semi metal Bi, small differences being accounted for by their contrasting "size" and "softness" (defined here). In order to separate structural and chemical effects unambiguously if not uniquely, we model the embrittlement process by taking the system of grain boundary and free surfaces through a sequence of precisely defined gedanken processes; each of these representing a putative mechanism. We thereby identify three mechanisms of embrittlement by substitutional impurities, two of which survive in the case of embrittlement or cohesion enhancement by interstitials. Two of the three are purely structural and the third contains both structural and chemical elements that by their very nature cannot be further unravelled. We are able to take the systems we study through each of these stages by explicit computer simulations and assess the contribution of each to the nett reduction in intergranular cohesion. The conclusion we reach is that embrittlement by both Bi and Na is almost exclusively structural in origin; that is, the embrittlement is a size effect.
△ Less
Submitted 10 September, 2006; v1 submitted 23 August, 2006;
originally announced August 2006.
-
A Stabilization Mechanism of Zirconia Based on Oxygen Vacancies Only
Authors:
Stefano Fabris,
Anthony T. Paxton,
Michael W. Finnis
Abstract:
The microscopic mechanism leading to stabilization of cubic and tetragonal forms of zirconia (ZrO$_2$) is analyzed by means of a self-consistent tight-binding model. Using this model, energies and structures of zirconia containing different vacancy concentrations are calculated, equivalent in concentration to the charge compensating vacancies associated with dissolved yttria (Y$_2$O$_3$) in the…
▽ More
The microscopic mechanism leading to stabilization of cubic and tetragonal forms of zirconia (ZrO$_2$) is analyzed by means of a self-consistent tight-binding model. Using this model, energies and structures of zirconia containing different vacancy concentrations are calculated, equivalent in concentration to the charge compensating vacancies associated with dissolved yttria (Y$_2$O$_3$) in the tetragonal and cubic phase fields (3.2 and 14.4% mol respectively). The model is shown to predict the large relaxations around an oxygen vacancy, and the clustering of vacancies along the $<111 >$ directions, in good agreement with experiments and first principles calculations. The vacancies alone are shown to explain the stabilization of cubic zirconia, and the mechanism is analyzed.
△ Less
Submitted 6 June, 2002;
originally announced June 2002.
-
Ordering at Solid-Liquid Interfaces Between Dissimilar Materials
Authors:
Adham Hashibon,
Joan Adler,
Michael W. Finnis,
Wayne D. Kaplan
Abstract:
In an earlier report we explored structural correlations at a liquid-solid interface with molecular dynamics simulations of a model aluminium system using the Ercolessi-Adams potential and up to 4320 atoms. Substrate atoms were pinned to their equilibrium fcc crystalline positions while liquid atoms were free to move. A direct correlation between the amount of ordering in the liquid phase and th…
▽ More
In an earlier report we explored structural correlations at a liquid-solid interface with molecular dynamics simulations of a model aluminium system using the Ercolessi-Adams potential and up to 4320 atoms. Substrate atoms were pinned to their equilibrium fcc crystalline positions while liquid atoms were free to move. A direct correlation between the amount of ordering in the liquid phase and the underlying substrate orientation was found. In the present paper we extend this study to the case of a fixed bcc substrate in contact with liquid aluminium. We find surprisingly similar results for the density profiles of both (100) and (110) substrates. However, there is a far greater in-plane ordering in the (100) than for the (110) system. For the (100) substrates we observe adsorption of liquid atoms into the terminating plane of the bcc (100) substrate, effectivel
△ Less
Submitted 30 December, 2001;
originally announced December 2001.
-
Atomistic Study of Structural Correlations at a Liquid-Solid Interface
Authors:
Adham Hashibon,
Joan Adler,
Michael W. Finnis,
Wayne D. Kaplan
Abstract:
Structural correlations at a liquid-solid interface were explored with molecular dynamics simulations of a model aluminium system using the Ercolessi-Adams potential and up to 4320 atoms. Substrate atoms were pinned to their equilibrium crystalline positions while liquid atoms were free to move. The density profile at the interface was investigated for different substrate crystallographic orient…
▽ More
Structural correlations at a liquid-solid interface were explored with molecular dynamics simulations of a model aluminium system using the Ercolessi-Adams potential and up to 4320 atoms. Substrate atoms were pinned to their equilibrium crystalline positions while liquid atoms were free to move. The density profile at the interface was investigated for different substrate crystallographic orientations and temperatures. An exponential decay of the density profile was observed, $ρ(z) \sim \rm{e}^{- κz}$, leading to the definition of $κ$ as a quantitative measure of the ordering at the liquid solid interface. A direct correlation between the amount of ordering in the liquid phase and the underlying substrate orientation was found.
△ Less
Submitted 30 December, 2001;
originally announced December 2001.
-
Free energy and molecular dynamics calculations for the cubic-tetragonal phase transition in zirconia
Authors:
Stefano Fabris,
Anthony T. Paxton,
Michael W. Finnis
Abstract:
The high-temperature cubic-tetragonal phase transition of pure stoichiometric zirconia is studied by molecular dynamics (MD) simulations and within the framework of the Landau theory of phase transformations. The interatomic forces are calculated using an empirical, self-consistent, orthogonal tight-binding (SC-TB) model, which includes atomic polarizabilities up to the quadrupolar level. A firs…
▽ More
The high-temperature cubic-tetragonal phase transition of pure stoichiometric zirconia is studied by molecular dynamics (MD) simulations and within the framework of the Landau theory of phase transformations. The interatomic forces are calculated using an empirical, self-consistent, orthogonal tight-binding (SC-TB) model, which includes atomic polarizabilities up to the quadrupolar level. A first set of standard MD calculations shows that, on increasing temperature, one particular vibrational frequency softens. The temperature evolution of the free energy surfaces around the phase transition is then studied with a second set of calculations. These combine the thermodynamic integration technique with constrained MD simulations. The results seem to support the thesis of a second-order phase transition but with unusual, very anharmonic behaviour above the transition temperature.
△ Less
Submitted 28 February, 2001;
originally announced February 2001.
-
Relative energetics and structural properties of zirconia using a self-consistent tight-binding model
Authors:
Stefano Fabris,
Anthony T. Paxton,
Michael W. Finnis
Abstract:
We describe an empirical, self-consistent, orthogonal tight-binding model for zirconia, which allows for the polarizability of the anions at dipole and quadrupole levels and for crystal field splitting of the cation d orbitals. This is achieved by mixing the orbitals of different symmetry on a site with coupling coefficients driven by the Coulomb potentials up to octapole level. The additional f…
▽ More
We describe an empirical, self-consistent, orthogonal tight-binding model for zirconia, which allows for the polarizability of the anions at dipole and quadrupole levels and for crystal field splitting of the cation d orbitals. This is achieved by mixing the orbitals of different symmetry on a site with coupling coefficients driven by the Coulomb potentials up to octapole level. The additional forces on atoms due to the self-consistency and polarizabilities are exactly obtained by straightforward electrostatics, by analogy with the Hellmann-Feynman theorem as applied in first-principles calculations. The model correctly orders the zero temperature energies of all zirconia polymorphs. The Zr-O matrix elements of the Hamiltonian, which measure covalency, make a greater contribution than the polarizability to the energy differences between phases. Results for elastic constants of the cubic and tetragonal phases and phonon frequencies of the cubic phase are also presented and compared with some experimental data and first-principles calculations. We suggest that the model will be useful for studying finite temperature effects by means of molecular dynamics.
△ Less
Submitted 13 February, 2000;
originally announced February 2000.
-
Equilibrium and adhesion of Nb/sapphire: the effect of oxygen partial pressure
Authors:
Iskander G. Batirev,
Ali Alavi,
Michael W. Finnis
Abstract:
We derive a formula, useful for first-principles calculations, which relates the free energy of an oxide/metal interface to the free energies of surfaces and the work of separation of the interface. We distinguish the latter {\it mechanical} quantity from the {\it thermodynamic} work of adhesion, and we describe explicitly how both may be calculated. Our formulae for interfacial and surface ener…
▽ More
We derive a formula, useful for first-principles calculations, which relates the free energy of an oxide/metal interface to the free energies of surfaces and the work of separation of the interface. We distinguish the latter {\it mechanical} quantity from the {\it thermodynamic} work of adhesion, and we describe explicitly how both may be calculated. Our formulae for interfacial and surface energies are cast in terms of quantities which can be calculated or looked up in tables, and include as additional parameters the ambient temperature and partial pressure of oxygen $P_{O_2}$. From total energy calculations for the Nb(111)/alpha-alumina(0001) interface, free Nb and alumina surfaces, we obtain firstly numerical estimates of the works of separation, which are independent of $P_{O_2}$. We then obtain surface energies, interfacial energies and the equilibrium work of adhesion as a function of $P_{O_2}$.
△ Less
Submitted 12 January, 2000;
originally announced January 2000.