Professors
Readers
Lecturers
Academic Fellows
Research Fellows
Postgraduate Students

         > Title: Universal tight binding model for chemical reactions in solution and at surfaces. III. Stoichiometric and reduced surfaces of titania and the adsorption of water Author(s): Lozovoi1 A.Y., Pashov D.L., Sheppard T.J., Kohanoff J.J., Paxton A.T. Journal of Chemical Physics, 141, pp. 044505 (2014) doi:
We demonstrate a model for stoichiometric and reduced titanium dioxide intended for use in molecular dynamics and other atomistic simulations and based in the polarizable ion tight binding theory. This extends the model introduced in two previous papers from molecular and liquid applications into the solid state, thus completing the task of providing a comprehensive and unified scheme for studying chemical reactions, particularly aimed at problems in catalysis and electrochemistry. As before, experimental results are given priority over theoretical ones in selecting targets for model fitting, for which we used crystal parameters and band gaps of titania bulk polymorphs, rutile and anatase. The model is applied to six low index titania surfaces, with and without oxygen vacancies and adsorbed water molecules, both in dissociated and nondissociated states. Finally, we present the results of molecular dynamics simulation of an anatase cluster with a number of adsorbed water molecules and discuss the role of edge and corner atoms of the cluster. Title: Universal tight binding model for chemical reactions in solution and at surfaces. II. Water Author(s): Lozovoi1 A.Y., Pashov D.L., Sheppard T.J., Kohanoff J.J., Paxton A.T. Journal of Chemical Physics, 141, pp. 044504 (2014) doi:
A revised water model intended for use in condensed phase simulations in the framework of the self consistent polarizable ion tight binding theory is constructed. The model is applied to water monomer, dimer, hexamers, ice, and liquid, where it demonstrates good agreement with theoretical results obtained by more accurate methods, such as DFT and CCSD(T), and with experiment. In particular, the temperature dependence of the self diffusion coefficient in liquid water predicted by the model, closely reproduces experimental curves in the temperature interval between 230 K and 350 K. In addition, and in contrast to standard DFT, the model properly orders the relative densities of liquid water and ice. A notable, but inevitable, shortcoming of the model is underestimation of the static dielectric constant by a factor of two. We demonstrate that the description of inter and intramolecular forces embodied in the tight binding approximation in quantum mechanics leads to a number of valuable insights which can be missing from ab initio quantum chemistry and classical force fields. These include a discussion of the origin of the enhanced molecular electric dipole moment in the condensed phases, and a detailed explanation for the increase of coordination number in liquid water as a function of temperature and compared with ice—leading to insights into the anomalous expansion on freezing. The theory holds out the prospect of an understanding of the currently unexplained density maximum of water near the freezing point. Title: Universal tight binding model for chemical reactions in solution and at surfaces. I. Organic molecules Author(s): Lozovoi1 A.Y., Pashov D.L., Sheppard T.J., Kohanoff J.J., Paxton A.T. Journal of Chemical Physics, 141, pp. 044503 (2014) doi:
As is now well established, a first order expansion of the Hohenberg–Kohn total energy density functional about a trial input density, namely, the Harris–Foulkes functional, can be used to rationalize a non self consistent tight binding model. If the expansion is taken to second order then the energy and electron density matrix need to be calculated self consistently and from this functional one can derive a charge self consistent tight binding theory. In this paper we have used this to describe a polarizable ion tight binding model which has the benefit of treating charge transfer in point multipoles. This admits a ready description of ionic polarizability and crystal field splitting. It is necessary in constructing such a model to find a number of parameters that mimic their more exact counterparts in the density functional theory. We describe in detail how this is done using a combination of intuition, exact analytical fitting, and a genetic optimization algorithm. Having obtained model parameters we show that this constitutes a transferable scheme that can be applied rather universally to small and medium sized organic molecules. We have shown that the model gives a good account of static structural and dynamic vibrational properties of a library of molecules, and finally we demonstrate the model's capability by showing a real time simulation of an enolization reaction in aqueous solution. In two subsequent papers, we show that the model is a great deal more general in that it will describe solvents and solid substrates and that therefore we have created a self consistent quantum mechanical scheme that may be applied to simulations in heterogeneous catalysis. Title: An ignition key for atomicscale engines Author(s): Dundas D., Cunningham B., Buchanan C., Terasawa A., Anthony T Paxton A.T., Todorov T.N. Journal of Physics: Condensed Matter, 24, pp. 40220314022036 (2012) doi:
A currentcarrying resonant nanoscale device, simulated by nonadiabatic
molecular dynamics, exhibits sharp activation of nonconservative
currentinduced forces with bias. The result, above the critical bias,
is generalized rotational atomic motion with a large gain in kinetic energy.
The activation exploits
sharp features in the electronic structure, and constitutes, in effect, an
ignition key for atomicscale motors. A controlling factor for the effect
is the nonequilibrium dynamical response matrix for smallamplitude atomic motion
under current. This matrix can be found from the steadystate electronic structure
by a simpler static calculation, providing a way to detect the likely appearance, or
otherwise, of nonconservative dynamics, in advance of realtime modelling.
Title: Is the pinning of ordinary dislocations in γTiAl intrinsic or extrinsic in nature? A combined atomistic and kinetic Monte Carlo approach. Author(s): Katzarov I.H., Paxton A.T. Acta Materialia, 59, No. 3, pp. 12811290 (2011) doi:
We address the question of the observed pinning of 1/2 <110] ordinary screw dislocations in γTiAl which leads to the characteristic trailing of dipoles in the microstructure. While it has been proposed that these may be variously intrinsic or extrinsic in nature, we are able to rule out the former mechanism. We do this by means of very large scale, three dimensional atomistic simulations using the quantum mechanical bond order potential. We find that the kinkpair formation energy is large: 6eV, while the single kink migration energy is conversely very small: 0.13eV. Using these, and other atomistically derived data, we make kinetic Monte Carlo simulations at realistic time and length scales to simulate dislocation mobility as a function of stress and temperature. In the temperature range of the stress anomaly in γTiAl, we determine whether one or several of the pinning and unzipping processes associated with generation of jogs are observed during our simulations. We conclude that the pinning of ordinary dislocations and anomalous mechanical behaviour in γTiAl must be attributed to a combination of extrinsic obstacles and extensive crossslip in a crystal containing impurities. Title: A tight binding model for water Author(s): Paxton A.T., Kohanoff J.J. Journal of Chemical Physics, 134, No. 4, Art. No. 044130 (2011) doi:
We demonstrate for the first time a tight binding model for water incorporating polarizable oxygen atoms. A novel aspect is that we adopt a "ground up" approach in that properties of the monomer and dimer only are fitted. Subsequently we make predictions of the structure and properties of hexamer clusters, iceXI and liquid water. A particular feature, missing in current tight binding and semiempirical hamiltonians, is that we reproduce the almost twofold increase in molecular dipole moment as clusters are built up towards the limit of bulk liquid. We concentrate on properties of liquid water, particularly dielectric constant and self diffusion coefficient, which are very well rendered in comparison with experiment. Finally we comment on the question of the contrasting densities of water and ice which is central to an understanding of the subtleties of the hydrogen bond. Title: Nonconservative currentinduced forces: A physical interpretation
Author(s): Todorov T.N., Dundas D., Paxton A.T., Horsfield A.P. Beilstein Journal of Nanotechnology, 2, pp. 727733 (2011) doi:
We give a physical interpretation of the recently demonstrated nonconservative nature of interatomic forces in currentcarrying nanostructures. We start from the analytical expression for the curl of these forces, and evaluate it for a point defect in a currentcarrying system. We obtain a general definition of the capacity of electrical current flow to exert a nonconservative force, and thus do net work around closed paths, by a formal noninvasive test procedure. Second, we show that the gain in atomic kinetic energy over time, generated by nonconservative currentinduced forces, is equivalent to the uncompensated stimulated emission of directional phonons. This connection with electronâphonon interactions quantifies explicitly the intuitive notion that nonconservative forces work by angular momentum transfer.
Title: Microscopic Origin of Channeled Flow in Lamellar Titanium Aluminide Author(s): Katzarov I.H., Paxton A.T. Physical Review Letters, 104, No. 22, Art. No. 225502 (2010) doi:
We employ a quantum mechanical bond order potential in an atomistic simulation of channeled flow. We show that the original hypothesis that this is achieved by a cooperative deployment of slip and twinning is correct, first because a twin is able to “protect” a 60° ordinary dislocation from becoming sessile, and second because the two processes are found to be activated by Peierls stresses of similar magnitude. In addition we show an explicit demonstration of the lateral growth of a twin, again at a similar level of stress. Thus these simultaneous processes are shown to be capable of channeling deformation into the observed state of plane strain in socalled “A”oriented mechanical testing of titanium aluminide superalloy. Title: Electronic structure and total energy of interstitial hydrogen in iron: Tight binding models Author(s): Paxton A.T., Elsässer C. Physical Review B, 82, No. 23, Art. No. 23512 (2010) doi:
An application of the tight binding approximation is presented
for the description of electronic structure and interatomic force
in magnetic iron, both pure and containing hydrogen
impurities. We assess the simple canonical dband description in
comparison to a non orthogonal model including s and d bands.
The transferability of our models is tested against
known properties including the segregation energies of hydrogen
to vacancies and to surfaces of iron. In many cases agreement is
remarkably good, opening up the way to quantum mechanical
atomistic simulation of the effects of hydrogen on mechanical
properties.
Title: Ring currents in azulene
Author(s): Paxton A.T., Todorov T.N., Elena A.M. Chemical Physics Letters, 483, No. 13, pp. 154158 (2009) doi:
We propose a self consistent polarisable ion tight binding theory for the study of push–pull processes in
aromatic molecules. We ﬁnd that the method quantitatively reproduces ab initio calculations of dipole
moments and polarisability. We apply the scheme in a simulation which solves the time dependent
Schrödinger equation to follow the relaxation of azulene from the second excited to the ground states.
We observe rather spectacular oscillating ring currents which we explain in terms of interference
between the HOMO and LUMO states. Title: Atomistic studies of <101] screw dislocation core structures and glide in γTiAl Author(s): Katzarov I.H., Paxton A.T. Philosophical Magazine, 89, No. 21, pp. 17311750 (2009) doi:
The core structure of superdislocations in L10 TiAl was investigated with a view to clarifying their dissociation abilities and the mechanisms by which they may become sessile by selflocking. A detailed knowledge of the fine structure of dislocations is essential in analysing the origin of the various deformation features. Atomistic simulation of the core structure and glide of the screw superdislocation was carried out using a bond order potential for γTiAl. The core structure of the screw superdislocation was examined, starting with initial unrelaxed configurations corresponding to various dislocation dissociations discussed in the literature. The superdislocation was found to possess in the screw orientation either planar (glissile) or nonplanar (sessile) core structures. The response of the core configurations to externally applied shear stress was studied. Some implications were considered of the dissociated configurations and their response to externally applied stress on dislocation dynamics, including the issue of dislocation decomposition, the mechanism of locking and the orientation dependence of the dislocation substructure observed in singlephase γTiAl. An unexpectedly rich and complex set of candidate core structures, both planar and nonplanar, was found, the cores of which may transform under applied stress with consequent violation of Schmid's law. Title: Atomistic studies of interactions between the dominant
lattice dislocations and γ/γlamellar boundaries in lamellar γTiAl Author(s): Katzarov I.H., Paxton A.T. Acta Materialia, 57, No. 11, pp. 33493366 (2009) doi:
This paper reports on atomistic simulations of the interactions between the dominant lattice dislocations in γTiAl (ordinary screw 1/2<101]
and <101] superdislocations) with all three kinds of γ/γlamellar boundaries in polysynthetically twinned (PST) TiAl. The purpose of this study is to clarify the early stage of lamellar boundary controlled plastic deformation in PST TiAl. The interatomic inter
actions in our simulations are described by a bond order potential for L10TiAl which provides a proper quantum mechanical description
of the bonding. We are interested in the dislocation core geometries that the lattice produces in proximity to lamellar boundaries and the
way in which these cores are affected by the elastic and atomistic effects of dislocationlamellar boundary interaction. We study the way in
which the interfaces aﬀect the activation of ordinary dislocation and superdislocation slip inside the γlamellae and transfer of plastic
deformation across lamellar boundaries. We ﬁnd three new phenomena in the atomicscale plasticity of PST TiAl, particularly due to
elastic and atomic mismatch associated with the 60° and 120° γ/γinterfaces: (i) two new roles of the γ/γinterfaces, i.e. decomposition
of superdislocations within 120° and 60° interfaces and subsequent detachment of a single ordinary dislocation and (ii) blocking of ordin
ary dislocations by 60° and 120° interfaces resulting in the emission of a twinning dislocation. Title: An introduction to the tight binding approximationimplementation by diagonalisation Author(s): Paxton A.T. Winter School: Multiscale Simulation Methods in Molecular Sciences, Forschungszentrum Jülich (MAR 2009)
Title: Magnetic tight binding and the ironchromium enthalpy anomaly Author(s): Paxton A.T., Finnis M.W. Physical Review B, 77, No. 2, Art. No. 024428 (2008) doi:
We describe a self consistent magnetic tightbinding theory based
in an expansion of the HohenbergKohn 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 a fixed moment model. We employ a simple set of tightbinding
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 fixed moment 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.
Title: Boron in copper: A perfect misfit in the bulk and cohesion enhancer at a grain boundary Author(s): Lozovoi A.Y., Paxton A.T. Physical Review B, 77, No. 16, Art. No. 165413 (2008) doi:
Using first principles electronic structure methods, we calculate the effects of boron impurities in bulk copper and at surfaces and grain boundaries. We find that boron segregation to the Sigma5(310)[001] grain boundary should strengthen the boundary up to 1.5 ML coverage (15.24 at./nm2). The maximal effect is observed at 0.5 ML and corresponds to boron atoms filling exclusively grain boundary interstices. In copper bulk, B causes significant distortion both in interstitial and regular lattice sites, for which boron atoms are either too big or too small. The distortion is compensated to a large extent when the interstitial and substitutional boron combine together to form a strongly bound dumbbell. Our prediction is that bound boron impurities should appear in a sizable proportion if not dominate in most experimental conditions. A large discrepancy between calculated heats of solution and experimental terminal solubility of B in Cu is found, indicating either a significant failure of the density functional approach or, more likely, strongly overestimated solubility limits in the existing BCu phase diagram. Title: Atomistic study of ordinary 1/2 < 110] screw dislocations in singlephase and lamellar gammaTiAl Author(s): Katzarov I.H., Cawkwell M.J., Paxton A.T., Finnis M.W. Philosophical Magazine, 87, No. 12, pp. 17951809 (2007) doi:
Computer simulation of the core structure and glide of ordinary 1/2 [110] screw dislocations in singlephase L1(0) TiAl and in two lamellae forming a twin gamma/gammainterface has been performed using recently constructed BondOrder Potentials (BOPs). BOPs represent a semiempirical, numerically efficient scheme that works within the orthogonal tightbinding approximation and is able to capture the directionality of bonding. We have studied dislocation glide in perfect L10 TiAl and along a twin interface, transmission of an ordinary screw dislocation between lamellae, and the core structure, mobility and detachment of an interfacial 1/2 [110] screw dislocation from a twin boundary under applied shear stresses in directions parallel and perpendicular to a (111) plane. Our results show that the glide of ordinary 1/2 [110] straight screw dislocations under applied stresses in L1(0) TiAl is characterized by zigzag movement on two conjugated {111} planes. The nonplanar core of the 1/2 (110] screw dislocation is distorted asymmetrically when the elastic centre of the dislocation is close to a twin gamma/gammainterface and the dislocation moves on one of the (111) planes, depending on the magnitude of the corresponding Schmid factor. Ordinary dislocations become ordinary interfacial dislocations when they reach the interface. With increasing applied stress they can glide into the adjacent lamella, leaving no remnant interfacial dislocation. Title: Structural and chemical embrittlement of grain boundaries by impurities: A general theory and firstprinciples calculations for copper Author(s): Lozovoi A.Y., Paxton A.T., Finnis M.W. Physical Review B, 74, No. 15, Art. No. 155416 (2006) doi:
Firstprinciples 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 semimetal 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 unraveled. 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 net 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. Title: Stability of Sr adatom model structures for SrTiO3(001) surface reconstructions Author(s): Liborio L.M., Sanchez C.G., Paxton A.T., Finnis M.W. Journal Of PhysicsCondensed Matter, 17, No. 23, pp. L223L230 (2005) doi:
We report results of firstprinciples calculations on the thermodynamic stability of different Sr adatom structures that have been proposed to explain some of the observed reconstructions of the (001) surface of strontium titanate (Kubo and Nozoye 2003 Surf Sci. 542 177). From surface free energy calculations, a phase diagram is constructed indicating the range of conditions over which each structure is most stable. These results are compared with Kubo and Nozoye's experimental observations. It is concluded that low Sr adatom coverage structures can only be explained if the surface is far from equilibrium. Intermediate coverage structures are stable only if the surface is in or very nearly in equilibrium with the strontium oxide. Title: Theory of the near Kedge structure in electron energy loss spectroscopy Author(s): Paxton A.T. Journal Of Electron Spectroscopy And Related Phenomena, 143, No. 23, pp. 5164 (2005) doi:
Arguments are given that lead to a formalism for calculating near Kedge structure in electron energy loss spectroscopy (EELS). This is essentially a one electron picture, while many body effects may be introduced at different levels, such as the local density approximation to density functional theory or the GW approximation to the electron selfenergy. Calculations are made within the all electron LMTO scheme in crystals with complex atomic and electronic structures, and these are compared with experiment. (c) 2004 Elsevier B.V. All rights reserved. Title: Bismuth embrittlement of copper is an atomic size effect Author(s): Schweinfest R., Paxton A.T., Finnis M.W. Nature, 432, No. 7020, pp. 10081011 (DEC 23 2004) doi:
Embrittlement by the segregation of impurity elements to grain boundaries is one of a small number of phenomena that can lead to metallurgical failure by fast fracture(1). Here we settle a question that has been debated for over a hundred years(2): how can minute traces of bismuth in copper cause this ductile metal to fail in a brittle manner? Three hypotheses for Bi embrittlement of Cu exist: two assign an electronic effect to either a strengthening(3) or weakening(4) of bonds, the third postulates a simple atomic size effect(5). Here we report first principles quantum mechanical calculations that allow us to reject the electronic hypotheses, while supporting a size effect. We show that upon segregation to the grain boundary, the large Bi atoms weaken the interatomic bonding by pushing apart the Cu atoms at the interface. The resolution of the mechanism underlying grain boundary weakening should be relevant for all cases of embrittlement by oversize impurities. Title: On the solvation of Laspartic acid Author(s): Paxton A.T., Harper J.B. Molecular Physics, 102, No. 910, pp. 953958 (2004) doi:
We use molecular statics and dynamics to study the stability of Laspartic acid both in vacuo and solvated by polar and nonpolar molecules using density functional theory in the generalized gradient approximation. We find that structures stable in vacuo are unstable in aqueous solution and vice versa. From our simulations we are able to come to some conclusions about the mechanism of stabilisation of zwitterions by polar protic solvents, water and methanol. Title: SrTiO3(001)(2x1) reconstructions: Firstprinciples calculations of surface energy and atomic structure compared with scanning tunneling microscopy images Author(s): Johnston K., Castell M.R., Paxton A.T., Finnis M.W. Physical Review B, 70, No. 8, Art. No. 085415 (2004) doi:
(1x1) and (2x1) reconstructions of the (001) SrTiO3 surface were studied using the firstprinciples fullpotential linear muffintin orbital method. Surface energies were calculated as a function of TiO2 chemical potential, oxygen partial pressure and temperature. The (1x1) unreconstructed surfaces were found to be energetically stable for many of the conditions considered. Under conditions of very low oxygen partial pressure the (2x1) Ti2O3 reconstruction [Martin R. Castell, Surf. Sci. 505, 1 (2002)] is stable. The question as to why STM images of the (1x1) surfaces have not been obtained was addressed by calculating charge densities for each surface. These suggest that the (2x1) reconstructions would be easier to image than the (1x1) surfaces. The possibility that the presence of oxygen vacancies would destabilise the (1x1) surfaces was also investigated. If the (1x1) surfaces are unstable then there exists the further possibility that the (2x1) DLTiO2 reconstruction [Natasha Erdman Nature (London) 419, 55 (2002)] is stable in a TiO2rich environment and for p(O2)>10(18) atm. Title: Bandstructure approach to near edge structure Author(s): Paxton A.T., Craven A.J., Gregg J.M., McComb D.W. Journal Of MicroscopyOxford, 210, No. 1, pp. 3544 (2003) doi:
We review the current state of the art in EELS fingerprinting by computer simulation, focusing on the bandstructure approach to the problem. Currently calculations are made using a one electron theory, but we describe in principle the way to go beyond this to include final state effects. We include these effects within the one electron framework using the Slater transition state formula and assess the errors involved. Two examples are then given which illustrate the use of the one electron approximation within density functional theory. Our approach is to combine predicted atomic structure with predicted electronic structure to assist in fingerprinting of complex crystal structures. Title: A stabilization mechanism of zirconia based on oxygen vacancies only Author(s): Fabris S., Paxton A.T., Finnis M.W. Acta Materialia, 50, No. 20, pp. 51715178 (2002) doi:
The microscopic mechanism leading to stabilization of cubic and tetragonal forms of zirconia (ZrO2) is analyzed by means of a selfconsistent tightbinding 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 (Y2O3) 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. (C) 2002 Acta Materialia Inc. Published by Elsevier Science Ltd. All rights reserved. Title: Electron energyloss nearedge shape as a probe to investigate the stabilization of yttriastabilized zirconia Author(s): Ostanin S., Craven A.J., McComb D.W., Vlachos D., Alavi A., Paxton A.T., Finnis M.W. Physical Review B, 65, No. 22, Art. No. 224109 (2002) doi:
The electron energyloss nearedge structure (ELNES) at the O K edge has been studied in yttriastabilized zirconia (YSZ). The electronic structure of YSZ for compositions between 3 and 15 mol % Y2O3 has been computed using a pseudopotentialbased technique to calculate the local relaxations near the O vacancies. The results showed phase transition from the tetragonal to cubic YSZ at 10 mol % of Y2O3, reproducing experimental observations. Using the relaxed defect geometry, calculation of the ELNES was carried out using the fullpotential linear muffintin orbital method. The results show very good agreement with the experimental O Kedge signal, demonstrating the power of using ELNES to probe the stabilization mechanism in doped metal oxides. Title: Material effects on stressinduced defect generation in trenched silicononinsulator structures Author(s): Nevin W.A., Somasundram K., McCann P., Magee S., Paxton A.T. Journal Of The Electrochemical Society, 148, No. 11, pp. G649G654 (2001) doi:
We have investigated the influence of the material properties of the silicon device layer on the generation of defects, and in particular slip dislocations, in trenched and refilled fusionbonded silicononinsulator structures. A strong dependence of the ease of slip generation on the type of dopant species was observed, with the samples falling into three basic categories; heavily borondoped silicon showed ready slip generation, arsenic and antimonydoped material was fairly resistant to slip, while silicon moderately or lightly doped with phosphorous or boron gave intermediate behavior. The observed behavior appears to be controlled by differences in the dislocation generation mechanism rather than by dislocation mobility. The introduction of an implanted buried layer at the bonding interface was found to result in an increase in slip generation in the silicon, again with a variation according to the dopant species. Here, the greatest slip occurred for both boron and antimonyimplanted samples. The weakening of the implanted material may be related to the presence of a band of precipitates observed in the silicon near the bonding interface. (C) 2001 The Electrochemical Society. Title: Free energy and molecular dynamics calculations for the cubictetragonal phase transition in zirconia Author(s): Fabris S., Paxton A.T., Finnis M.W. Physical Review B, 63, No. 9, Art. No. 094101 (2001) doi:
The hightemperature cubictetragonal 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, selfconsistent, orthogonal tightbinding 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 freeenergy 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 secondorder phase transition but with unusual, very anharmonic behavior above the transition temperature. Title: Effect of relaxation on the oxygen Kedge electron energyloss nearedge structure in yttriastabilized zirconia Author(s): Ostanin S., Craven A.J., McComb D.W., Vlachos D., Alavi A., Finnis M.W., Paxton A.T. Physical Review B, 62, No. 22, pp. 1472814735 (2000) doi:
The electron energyloss nearedge structure (ELNES) at the oxygen Kedge has been investigated in a range of yttriastabilized zirconia (YSZ) materials. The electronic structure of the three polymorphs of pure ZrO2 and of the doped YSZ structure close to the 33 mol %Y2O3 composition have been calculated using a fullpotential linear muffintin orbital method (NFPLMTO) as well as a pseudopotential based technique. Calculations of the ELNES dipole transition matrix elements in the framework of the NFPLMTO scheme and inclusion of core hole screening within Slater's transition state theory enable the ELNES to be computed. Good agreement between the experimental and calculated ELNES is obtained for pure monoclinic ZrO2. The agreement is less good with the ideal tetragonal and cubic structures. This is because the inclusion of defects is essential in the calculation of the YSZ ELNES. If the model used contains ordered defects such as vacancies and metal Y planes, agreement between the calculated and experimental O Kedges is significantly improved. The calculations show how the five different O environments of Zr,Y,O, are connected with the features observed in the experimental spectra and demonstrate clearly the power of using ELNES to probe the stabilization mechanism in doped metal oxides. Title: Relative energetics and structural properties of zirconia using a selfconsistent tightbinding model Author(s): Fabris S., Paxton A.T., Finnis M.W. Physical Review B, 61, No. 10, pp. 66176630 (2000) doi:
We describe an empirical, selfconsistent, orthogonal tightbinding 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 selfconsistency and polarizabilities are exactly obtained by straightforward electrostatics, by analogy with the HellmannFeynman theorem as applied in firstprinciples calculations. The model correctly orders the zero temperature energies of all zirconia polymorphs. The ZrO 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 firstprinciples calculations. We suggest that the model will be useful for studying finite temperature effects by means of molecular dynamics. Title: The nearedge structure in energyloss spectroscopy: manyelectron and magnetic effects in transition metal nitrides and carbides Author(s): Paxton A.T., van Schilfgaarde M., MacKenzie M., Craven A.J. Journal Of PhysicsCondensed Matter, 12, No. 5, pp. 729750 (2000) doi:
We investigate the ability of the local density approximation (LDA) in density functional theory to predict the nearedge structure in electron energyloss spectroscopy in the dipole approximation. We include screening of the core hole within the LDA using Slater's transition state theory. We find that anion Kedge threshold energies are systematically overestimated by 4.22 +/ 0.44 eV in twelve transition metal carbides and nitrides in the rocksalt (B1) structure. When we apply this 'universal' manyelectron correction to energyloss spectra calculated within the transition state approximation to LDA, we find quantitative agreement with experiment to within one or two eV for TiC, TiN and VN. We compare our calculations to a simpler approach using a projected Mulliken density which honours the dipole selection rule, in place of the dipole matrix element itself. We find remarkably close agreement between these two approaches. Finally, we show an anomaly in the nearedge structure in CrN to be due to magnetic structure. In particular, we find that the N K edge in fact probes the magnetic moments and alignments of ther sublattice.
