1 / 19

Electronic Structure of PuO2: DFT+U Calculations

This study presents DFT+U calculations of the electronic structure of perfect and defective PuO2, a material widely used in commercial nuclear fuel. The results show the accurate description of radiation-induced defects and the reproduction of UPS spectra. Defect formation energies are also analyzed using GGA+U and SIC methods. The findings provide valuable insights for understanding and improving nuclear fuel behavior.

jmclain
Download Presentation

Electronic Structure of PuO2: DFT+U Calculations

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. DFT+U calculations of the electronic structure of perfect and defective PuO2Eugene Kotomin and Denis GryaznovLaboratory of Theoretical Physics and Computer Modelling at the Institute of Solid State Physics, University of Latvia Riga, Latvia CECAM workshop on Actinides, Manchester, June 2010

  2. Motivation • Solid solution (U,Pu)O2 (known as MOX) is nowadays widely used as a commercial nuclear fuel. This is why understanding of the basic physical and chemical properties of MOX and parent materials (UO2 and PuO2) including radiation-induced defects is of great importance. Accurate theoretical calculations of these 5f-electron materials is a challenge due to strong electron correlation effects. CECAM workshop on Actinides, Manchester, June 2010

  3. PuO2 first principles calculations • L.Petit et al, Science 301, 498 (2003) SIC LSD method, including defects • R.Martin et al, Phys Rev B 76, 033101 (2007): Gaussian code + Periodic Boundary Conditions, HSE (Heyd,Scuseria, Enzerhof) hybrid functionals • Jomard et al, Phys Rev B 78, 075125 (2008) GGA+U (PBE functional), ABINIT code CECAM workshop on Actinides, Manchester, June 2010

  4. Basic experimental data • Cubic fluorite structure, a=5.396 A • semiconductor, with a gap of 1.8 eV • Pu shows no magnetic moment (?) • 5f states dominate at the Fermi level – not reproduced in previous studies CECAM workshop on Actinides, Manchester, June 2010

  5. UPS spectraT. Gouder et al, Surf Sci Lett. 601, L 77 (2007) CECAM workshop on Actinides, Manchester, June 2010

  6. Present study • Spin-polarized GGA+U VASP calculations • PAW pseudopotentials, 78 e (Pu), 2e (O) in a core • a comparison of Dudarev and Liechtenstein exchange-correlation functionals: rotationally invariant Ueff= U – J vs independent U, J parameters • Supercells 3, 6, 12 (conventional cell) and 24 atoms, FM and AFM magnetic solutions CECAM workshop on Actinides, Manchester, June 2010

  7. 24 atom supercell O Pu CECAM workshop on Actinides, Manchester, June 2010

  8. DOS for FM PuO2 Dudarev functional for a) the cubic structure and b) the tetragonal structure. Spin-up electrons are only shown and the J-exchange parameter fixed at 0.5 eV. The Fermi level is taken as zero CECAM workshop on Actinides, Manchester, June 2010

  9. Dudarev Exc vs Liechtenstein:AFM structure AFM PuO2 a) the Dudarev functional at U = 6.0 eV, J=0.5 eV, b) Liechtenstein functional for different J with fixed U = 3.0 eV. Spin-up electrons are only shown. The Fermi level is taken as zero. CECAM workshop on Actinides, Manchester, June 2010

  10. Optimal parameters • Liechtenstein functional, U=3 eV, J=1.5 eV • a=5.51 A, c=5.41 A (AFM, tetragonal str.) • DOS corresponds to UPS spectra • The effective (Bader) charges: 2.48 e (Pu) -1.24 e (O)  considerable covalency! CECAM workshop on Actinides, Manchester, June 2010

  11. Hybrid (PBE0) calculations CECAM workshop on Actinides, Manchester, June 2010

  12. Defects in PuO2 • Change of supercell volume due to Pu and O vacancies (Å**3). Positive sign means increase of volume, negative its decrease. The values correspond to the complete local structure and lattice parameter optimisation. Numbers in brackets show the defect concentration. • Defect 12 atom cell 24 atom • Pu vacancy -12.60 (25%) -14.43 (12.5%) • O vacancy 4.58 (12.5%) Lattice constant variation • 0.8 % (O) vs -2.6% (Pu) • 0.9 % (SIC)-3.4% (SIC) (12% of defects) CECAM workshop on Actinides, Manchester, June 2010

  13. Local lattice, charge perturbations • Pu vacancy NN O ions displaced 1.1% charge -0.95e NNN Pu 0.0 2.56 e • O vacancy NN Pu ions displaced 1.2 % charge 2.20 e NNN O 0.0 -1.25e Charge redistribution is very local, only NN CECAM workshop on Actinides, Manchester, June 2010

  14. Defect formation energies, eV O vacancy GGA+U SIC O-rich 3.97 (3.60-24 at.) Pu-rich -1.08 -- O interstitial O-rich 1.78 1.90 Pu-rich 6.83 6.70 Pu vacancy O-rich 3.07 9.0 Pu-rich 13.18 18.5 CECAM workshop on Actinides, Manchester, June 2010

  15. Thermodynamic reference • Pu-rich conditions correspond to the chemical potential for excess Pu • μ(Pu) = - Ecoh(δ-Pu), • μ(O)=½( - Ecoh(PuO2)- μ(Pu)); • O-rich condition implies • μ(O) = - ½Ecoh(O2) • μ(Pu)= - Ecoh(PuO2)- 2μ(O). • Defect formation energy definition: • Ef = - Ecoh(PunOm) - n μ(Pu) – m μ(O). CECAM workshop on Actinides, Manchester, June 2010

  16. Analysis • Defect formation energy deacreases 10% as concentration increases from 12 to 25% • Reasonable agreement of GGA+U and SIC, discrepancy is due to difference in the reference states-cohesive energies of pure Pu, PuO2 and O2 molecule. CECAM workshop on Actinides, Manchester, June 2010

  17. CONCLUSIONS • Only use of the Liechtenstein functional within GGA+U permits to reproduce the UPS spectra • Defect calculations using GGA+U and SIC agree quite well • Hybrid functionals in VASP and GAUSSIAN-PBC give similar results, but this does not help! • Defect migration calculations is the next step in multi-scale study of nuclear fuel long-time performance (TRANSURANUS code) CECAM workshop on Actinides, Manchester, June 2010

  18. Many thanks to L.Petit, A.Svane, M.Freyss, R.A.Evarestov for many stumulating discussions, as well as EC FP7 F-Bridge project for a financial support CECAM workshop on Actinides, Manchester, June 2010

  19. Thank you for your attention! CECAM workshop on Actinides, Manchester, June 2010

More Related