A new reactive force field to describe proton diffusion within the solid oxide fuel cell material BaZrO3 has been derived. Using a quantum mechanical potential energy surface, the parameters of an interatomic potential model to describe hydroxyl groups within both pure and yttrium-doped BaZrO3 have been determined. Reactivity is then incorporated through the use of the empirical valence bond model. Molecular dynamics simulations (EVB-MD) have been performed to explore the diffusion of hydrogen using a stochastic thermostat and barostat whose equations are extended to the isostress–isothermal ensemble. In the low concentration limit, the presence of yttrium is found not to significantly influence the diffusivity of hydrogen, despite the proton having a longer residence time at oxygen adjacent to the dopant. This lack of influence is due to the fact that trapping occurs infrequently, even when the proton diffuses through octahedra adjacent to the dopant. The activation energy for diffusion is found to be 0.42 eV, in good agreement with experimental values, though the prefactor is slightly underestimated.

Reactive force field simulation of proton diffusion in BaZrO3 using an empirical valence bond approach / Raiteri, P.; Gale, J. D.; Bussi, G.. - In: JOURNAL OF PHYSICS. CONDENSED MATTER. - ISSN 0953-8984. - 23:33(2011), pp. 334213.1-334213.13. [10.1088/0953-8984/23/33/334213]

Reactive force field simulation of proton diffusion in BaZrO3 using an empirical valence bond approach

Bussi, G.
2011-01-01

Abstract

A new reactive force field to describe proton diffusion within the solid oxide fuel cell material BaZrO3 has been derived. Using a quantum mechanical potential energy surface, the parameters of an interatomic potential model to describe hydroxyl groups within both pure and yttrium-doped BaZrO3 have been determined. Reactivity is then incorporated through the use of the empirical valence bond model. Molecular dynamics simulations (EVB-MD) have been performed to explore the diffusion of hydrogen using a stochastic thermostat and barostat whose equations are extended to the isostress–isothermal ensemble. In the low concentration limit, the presence of yttrium is found not to significantly influence the diffusivity of hydrogen, despite the proton having a longer residence time at oxygen adjacent to the dopant. This lack of influence is due to the fact that trapping occurs infrequently, even when the proton diffuses through octahedra adjacent to the dopant. The activation energy for diffusion is found to be 0.42 eV, in good agreement with experimental values, though the prefactor is slightly underestimated.
2011
23
33
1
13
334213
https://arxiv.org/pdf/1104.0773.pdf
Raiteri, P.; Gale, J. D.; Bussi, G.
File in questo prodotto:
File Dimensione Formato  
0953-8984_23_33_334213.pdf

non disponibili

Tipologia: Versione Editoriale (PDF)
Licenza: Non specificato
Dimensione 1.4 MB
Formato Adobe PDF
1.4 MB Adobe PDF   Visualizza/Apri   Richiedi una copia

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/20.500.11767/11900
Citazioni
  • ???jsp.display-item.citation.pmc??? 3
  • Scopus 45
  • ???jsp.display-item.citation.isi??? 42
social impact