Relative energetics and structural properties of zirconia using a self-consistent tight-binding model

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.

A tight binding model for water

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, ice-XI and liquid water. A particular feature, missing in current tight binding and semiempirical hamiltonians, is that we reproduce the almost two-fold 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.

Universal tight binding model for chemical reactions in solution and at surfaces. II. Water

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. 

Universal tight binding model for chemical reactions in solution and at surfaces. III. Stoichiometric and reduced surfaces of titania and the adsorption of water

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 non-dissociated 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. (C) 2014 AIP Publishing LLC.

Universal tight binding model for chemical reactions in solution and at surfaces. I. Organic molecules

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.

Tight-binding simulation of current-carrying nanostructures

The tight-binding (TB) approach to the modelling of electrical conduction in small structures is introduced. Different equivalent forms of the TB expression for the electrical current in a nanoscale junction are derived. The use of the formalism to calculate the current density and local potential is illustrated by model examples. A first-principles time-dependent TB formalism for calculating current-induced forces and the dynamical response of atoms is presented. An earlier expression for current-induced forces under steady-state conditions is generalized beyond local charge neutrality and beyond orthogonal TB. Future directions in the modelling of power dissipation and local heating in nanoscale conductors are discussed.

Magnetic tight binding and the iron-chromium enthalpy anomaly

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-polarized 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 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 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.

Electronic structure and total energy of interstitial hydrogen in iron: Tight-binding models

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 d-band 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.

Time-dependent tight binding

Starting from a Lagrangian mean-field theory, a set of time-dependent tight-binding equations is derived to describe dynamically and self-consistently an interacting system of quantum electrons and classical nuclei. These equations conserve norm, total energy and total momentum. A comparison with other tight-binding models is made. A previous tight-binding result for forces on atoms in the presence of electrical current flow is generalized to the time-dependent domain and is taken beyond the limit of local charge neutrality.

A simple model of atomic interactions in noble metals based explicitly on electronic structure

A total energy tight-binding model with a basis of just one s state per atom is introduced. It is argued that this simplest of all tight-binding models provides a surprisingly good description of the structural stability and elastic constants of noble metals. By assuming inverse power scaling laws for the hopping integrals and the repulsive pair potential, it is shown that the density matrix in a perfect primitive crystal is independent of volume, and structural energy differences and equations of state are then derived analytically. The model is most likely to be of use when one wishes to consider explicitly and self-consistently the electronic and atomic structures of a generic metallic system, with the minium of computation expense. The relationship to the free-electron jellium model is described. The applicability of the model to other metals is also considered briefly.

A stabilization mechanism of zirconia based on oxygen vacancies only

The microscopic mechanism leading to stabilization of cubic and tetragonal forms of zirconia (ZrO2) 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 (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.

Free energy and molecular dynamics calculations for the cubic-tetragonal phase transition in zirconia

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 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 behavior above the transition temperature.

Inelastic current-voltage spectroscopy of atomic wires

A tight-binding model is developed to describe the electron-phonon coupling in atomic wires under an applied voltage and to model, their inelastic current-voltage spectroscopy. Particular longitudinal phonons are found to have greatly enhanced coupling to the electronic states of the system. This leads to a large drop in differential conductance at threshold energies associated with these phonons. It is found that with increasing tension these energies decrease, while the size of the conductance drops increases, in agreement with experiment.

Power dissipation in nanoscale conductors

A previous tight-binding model of power dissipation in a nanoscale conductor under an applied bias is extended to take account of the local atomic topology and the local electronic structure. The method is used to calculate the power dissipated at every atom in model nanoconductor geometries: a nanoscale constriction, a one-dimensional atomic chain between two electrodes with a resonant double barrier, and an irregular nanowire with sharp corners. The local power is compared with the local current density and the local density of states. A simple relation is found between the local power and the current density in quasiballistic geometries. A large enhancement in the power at special atoms is found in cases of resonant and anti-resonant transmission. Such systems may be expected to be particularly unstable against current-induced modifications.