Hybrid density functional theory applied to magnetite: Crystal structure, charge order, and phonons

Andrew D. Rowan, Charles H. Patterson, and L. V. Gasparov
Phys. Rev. B 79, 205103 – Published 5 May 2009
PDFHTMLExport Citation

Abstract

The electronic structure and equilibrium structure of magnetite (Fe3O4) in the high temperature cubic Fd3¯m and low temperature monoclinic P2/c unit cells have been computed using the Perdew-Wang generalized gradient approximation (GGA) to density functional theory (DFT) and the B3LYP hybrid density functional. The ground state for the GGA-DFT is an itinerant electron metallic state in the cubic unit cell and the ground state for the B3LYP functional is a charge ordered semiconducting state in the monoclinic unit cell. The equilibrium structure predicted by the B3LYP functional for Fe3O4 in the P2/c unit cell has been calculated with lattice parameters fixed at values obtained in recent x-ray diffraction work and with the lattice fully relaxed. Bond lengths obtained with lattice parameters fixed at experimental values are in excellent agreement with x-ray measurements [J. P. Wright et al., Phys. Rev. B 66, 214422 (2002)]. The degree of charge order, measured as disproportionation of charge on octahedral B sites, is considerably less than unity and in reasonable agreement with values from resonant x-ray diffraction measurements. However, conduction electrons are found to be fully localized on B1 and B4 sites in orbitally ordered t2g states. This shows that they are formally Fe2+ ions while Fe B2 and B3 sites are formally Fe3+ sites. Therefore Verwey’s original conjecture regarding charge localization in Fe3O4 applies, even though the specific pattern of charge order is different. GGA-DFT and B3LYP density functionals were used to calculate phonons at the Γ point of the Brillouin zone. Phonon frequencies predicted for these crystal structures are compared to frequencies from infrared conductivity and Raman scattering experiments. Charge ordering causes symmetry breaking of force constants on symmetry lowering from the cubic Fd3¯m unit cell to the P2/c unit cell. This produces frequency splitting of modes which are degenerate in the cubic unit cell and concentration of ion displacements in phonon eigenvectors on particular Fe octahedral B site chains, especially in the highest frequency bands.

  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
9 More
  • Received 18 November 2008

DOI:https://doi.org/10.1103/PhysRevB.79.205103

©2009 American Physical Society

Authors & Affiliations

Andrew D. Rowan and Charles H. Patterson

  • School of Physics, Trinity College Dublin, Dublin 2, Ireland

L. V. Gasparov

  • Department of Chemistry and Physics, University of North Florida, St. John’s Bluff Road South, Jacksonville, Florida 32224, USA

Article Text (Subscription Required)

Click to Expand

Supplemental Material (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand
Issue

Vol. 79, Iss. 20 — 15 May 2009

Reuse & Permissions
Access Options
Author publication services for translation and copyediting assistance advertisement

Authorization Required


×
×

Images

×

Sign up to receive regular email alerts from Physical Review B

Log In

Cancel
×

Search


Article Lookup

Paste a citation or DOI

Enter a citation
×