Dimitar Pashov

Fully quantum mechanical calculation of the diffusivity of hydrogen in iron using the tight-binding approximation and path integral theory

We present calculations of free energy barriers and diffusivities as functions of temperature for the diffusion of hydrogen in bcc-Fe. This is a fully quantum mechanical approach since the total energy landscape is computed using a new self consistent, transferable tight binding model for interstitial impurities in magnetic iron. Also the hydrogen nucleus is treated quantum mechanically and we compare here two approaches in the literature both based in the Feynman path integral formulation of statistical mechanics. We find that the quantum transition state theory which admits greater freedom for the proton to explore phase space gives result in better agreement with experiment than the alternative which is based on fixed centroid calculations of the free energy barrier. We also find results in better agreement compared to recent centroid molecular dynamics (CMD) calculations of the diffusivity which employed a classical interatomic potential rather than our quantum mechanical tight binding theory. In particular we find first that quantum effects persist to higher temperatures than previously thought, and conversely that the low temperature diffusivity is smaller than predicted in CMD calculations and larger than predicted by classical transition state theory. This will have impact on future modeling and simulation of hydrogen trapping and diffusion.

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

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.

Accurate optical properties from first principles: a Quasiparticle Self consistent GW plus Bethe-Salpeter Equation approach

We present an approach to calculate the optical absorption spectra that combines the quasiparticle selfconsistent GW method [Phys. Rev. B, 76 165106 (2007)] for the electronic structure with the solution of the ladder approximation to the Bethe-Salpeter equation for the macroscopic dielectric function. The solution of the Bethe-Salpeter equation has been implemented within an all-electron framework, using a linear muffin-tin orbital basis set, with the contribution from the non-local self-energy to the transition dipole moments (in the optical limit) evaluated explicitly. This approach addresses those systems whose electronic structure is poorly described within the standard perturbative GW approaches with as a starting point density-functional theory calculations. The merits of this approach have been exemplified by calculating optical absorption spectra of a strongly correlated transition metal oxide, NiO, and a narrow gap semiconductor, Ge. In both cases, the calculated spectrum is in good agreement with the experiment. It is also shown that for systems whose electronic structure is well-described within the standard perturbative GW , such as Si, LiF and h-BN, the performance of the present approach is in general comparable to the standard GW plus Bethe-Salpeter equation. It is argued that both vertex corrections to the electronic screening and the electron-phonon interaction are responsible for the observed systematic overestimation of the fundamental bandgap and spectrum onset.

Theoretical evaluation of the role of crystal defects on local equilibrium and effective diffusivity of hydrogen in iron

ABSTRACT Hydrogen diffusion and trapping in ferrite is evaluated by quantum mechanically informed kinetic Monte Carlo simulations in defective microstructures. We find that the lattice diffusivity is attenuated by two to four orders of magnitude due to the presence of dislocations. We also find that pipe diffusivity is vanishingly small along screw dislocations and demonstrate that dislocations do not provide fast diffusion pathways for hydrogen as is sometimes supposed. We make contact between our simulations and the predictions of Orianis theory of ‘effective diffusivity’, and find that local equilibrium is maintained between lattice and trap sites. We also find that the predicted effective diffusivity is in agreement with our simulated results in cases where the distribution of traps is spatially homogeneous; in the trapping of hydrogen by dislocations where this condition is not met, the Oriani effective diffusivity is in agreement with the simulations to within a factor of two. This paper is part of a thematic issue on Hydrogen in Metallic Alloys