1 / 42

Magnetic and charge order phase transition in YBaFe 2 O 5 (Verwey transition)

Magnetic and charge order phase transition in YBaFe 2 O 5 (Verwey transition). Peter Blaha , Ch. Spiel , K.Schwarz Institute of Materials Chemistry TU Wien Thanks to P.Karen (Univ. Oslo, Norway). E.Verwey, Nature 144, 327 (1939). Fe 3 O 4 , magnetite.

marlie
Download Presentation

Magnetic and charge order phase transition in YBaFe 2 O 5 (Verwey transition)

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. Magnetic and charge order phase transition in YBaFe2O5 (Verwey transition) Peter Blaha, Ch. Spiel, K.Schwarz Institute of Materials Chemistry TU Wien Thanks to P.Karen (Univ. Oslo, Norway)

  2. E.Verwey, Nature 144, 327 (1939) Fe3O4, magnetite phase transition between a mixed-valence and a charge-ordered configuration 2 Fe2.5+ Fe2+ + Fe3+ cubic inverse spinel structure AB2O4 Fe2+A (Fe3+,Fe3+)B O4 Fe3+A (Fe2+,Fe3+)B O4 B A  small, but complicated coupling between lattice and charge order

  3. Double-cell perovskites: RBaFe2O5 ABO3 O-deficient double-perovskite Ba Y (R) square pyramidal coordination Antiferromagnet with a 2 step Verwey transition around 300 K Woodward&Karen, Inorganic Chemistry 42, 1121 (2003)

  4. structural changes in YBaFe2O5 • above TN (~430 K): tetragonal (P4/mmm) • 430K: slight orthorhombic distortion (Pmmm) due to AFM • all Fe in class-III mixed valence state +2.5; • ~334K: dynamic charge order transition into class-II MV state, • visible in calorimetry and • Mössbauer, but not with X-rays • 308K: complete charge order into • class-I MV state (Fe2+ + Fe3+) • large structural changes (Pmma) • due to Jahn-Teller distortion; • change of magnetic ordering: • direct AFM Fe-Fe coupling vs. • FM Fe-Fe exchange above TV

  5. structural changes CO structure: Pmma VM structure: Pmmm a:b:c=2.09:1:1.96 (20K) a:b:c=1.003:1:1.93 (340K) • Fe2+ and Fe3+ chains along b • contradicts Anderson charge-ordering conditions with minimal electrostatic repulsion (checkerboard like pattern) • has to be compensated by orbital ordering and e--lattice coupling c b a

  6. antiferromagnetic structure CO phase: G-type AFM VM phase: • AFM arrangement in all directions, AFM for all Fe-O-Fe superexchange paths also across Y-layer FM across Y-layer (direct Fe-Fe exchange) • Fe moments in b-direction 4 8 independent Fe atoms

  7. Theoretical methods: WIEN2K An Augmented Plane Wave Plus Local Orbital Program for Calculating Crystal Properties Peter Blaha Karlheinz Schwarz Georg Madsen Dieter Kvasnicka Joachim Luitz http://www.wien2k.at • WIEN2k (APW+lo) calculations • Rkmax=7, 100 k-points • spin-polarized, various spin-structures • + spin-orbit coupling • based on density functional theory: • LSDA or GGA (PBE) • Exc≡ Exc(ρ, ∇ρ) • description of “highly correlated electrons” using “non-local” (orbital dep.) functionals • LDA+U, GGA+U • hybrid-DFT (only for correlated electrons) • mixing exact exchange (HF) + GGA 1400 registered groups 2000 mailinglist users

  8. GGA-results: • Metallic behaviour/No bandgap • Fe-dn t2g states not splitted at EF • overestimated covalency between O-p and Fe-eg • Magnetic moments too small • Experiment: • CO: 4.15/3.65 (for Tb), 3.82 (av. for Y) • VM: ~3.90 • Calculation: • CO: 3.37/3.02 • VM: 3.34 • no significant charge order • charges of Fe2+ and Fe3+ sites nearly identical • CO phase less stable than VM • LDA/GGA NOT suited for this compound! Fe-eg t2g eg*t2g eg

  9. “Localized electrons”: GGA+U • Hybrid-DFT • ExcPBE0 [r] = ExcPBE [r] + a (ExHF[Fsel] – ExPBE[rsel]) • LDA+U, GGA+U • ELDA+U(r,n) = ELDA(r) + Eorb(n) – EDCC(r) • separate electrons into “itinerant” (LDA) and localized e-(TM-3d, RE 4f e-) • treat them with “approximate screened Hartree-Fock” • correct for “double counting” • Hubbard-U describes coulomb energy for 2e- at the same site • orbital dependent potential

  10. Determination of U • Take Ueff as “empirical” parameter (fit to experiment) • Estimate Ueff from constraint LDA calculations • constrain the occupation of certain states (add/subtract e-) • switch off any hybridization of these states (“core”-states) • calculate the resulting Etot • we used Ueff=7eV for all calculations

  11. DOS: GGA+U vs. GGA GGA+U GGA single lower Hubbard-band in VM splits in CO with Fe3+ states lower than Fe2+ insulator, t2g band splits metallic

  12. magnetic moments and band gap • magnetic moments in very good agreement with exp. • LDA/GGA: CO: 3.37/3.02 VM: 3.34 mB • orbital moments small (but significant for Fe2+) • band gap: smaller for VM than for CO phase • exp: semiconductor (like Ge); VM phase has increased conductivity • LDA/GGA: metallic

  13. Charge transfer • Charges according to Baders “Atoms in Molecules” theory • Define an “atom” as region within a zero flux surface • Integrate charge inside this region

  14. Structure optimization (GGA+U) O2a • CO phase: • Fe2+:shortest bond in y (O2b) • Fe3+: shortest bond in z (O1) • VM phase: • all Fe-O distances similar • theory deviates along z !! • Fe-Fe interaction • different U ?? • finite temp. ?? O2b O3 O1 O1

  15. Can we understand these changes ? • Fe2+ (3d6) CO Fe3+ (3d5) VM Fe2.5+ (3d5.5) • majority-spin fully occupied • strong covalency effects very localized states at lower energy than Fe2+ in eg and d-xz orbitals • minority-spin states • d-xz fully occupied (localized) empty d-z2 partly occupied • short bond in y short bond in z (one O missing) FM Fe-Fe; distances in z ??

  16. Difference densities Dr=rcryst-ratsup • CO phase VM phase Fe2+: d-xz Fe3+: d-x2 O1 and O3: polarized toward Fe3+ Fe: d-z2 Fe-Fe interaction O: symmetric

  17. dxz spin density (rup-rdn) of CO phase • Fe3+: no contribution • Fe2+: dxz • weak p-bond with O • tilting of O3 p-orbital

  18. Mössbauer spectroscopy: • Isomer shift: d = a (r0Sample – r0Reference); a=-.291 au3mm s-1 • proportional to the electron density r at the nucleus • Magnetic Hyperfine fields: Btot=Bcontact + Borb + Bdip • Bcontact = 8p/3 mB [rup(0) – rdn(0)] … spin-density at the nucleus … orbital-moment … spin-moment S(r) is reciprocal of the relativistic mass enhancement

  19. Electric field gradients (EFG) • Nuclei with a nuclear quantum number I≥1 have an electrical quadrupole moment Q • Nuclear quadrupole interaction (NQI) between “non-spherical” nuclear charge Q times the electric field gradient F • Experiments • NMR • NQR • Mössbauer • PAC EFG traceless tensor with traceless |Vzz| |Vyy| |Vxx| EFG Vzz asymmetry parameter principal component

  20. theoretical EFG calculations: EFG is tensor of second derivatives of VC at the nucleus: Cartesian LM-repr. EFG is proportial to differences of orbital occupations

  21. Mössbauer spectroscopy Hyperfine fields: Fe2+ has large Borb and Bdip Isomer shift: charge transfer too small in LDA/GGA EFG: Fe2+ has too small anisotropy in LDA/GGA CO VM

  22. magnetic interactions • CO phase: • magneto-crystalline anisotropy: moments point into y-direction in agreement with exp. • experimental G-type AFM structure (AFM direct Fe-Fe exchange) is 8.6 meV/f.u. more stable than magnetic order of VM phase (direct FM) • VM phase: • experimental “FM across Y-layer” AFM structure (FM direct Fe-Fe exchange) is 24 meV/f.u. more stable than magnetic order of CO phase (G-type AFM)

  23. Exchange interactions Jij • Heisenberg model: H = Si,j JijSi.Sj • 4 different superexchange interactions (Fe-Fe exchange interaction mediated by an O atom) • J22b : Fe2+-Fe2+ along b • J33b : Fe3+-Fe3+ along b • J23c : Fe2+-Fe3+ along c • J23a : Fe2+-Fe3+ along a • 1 direct Fe-Fe interaction • Jdirect: Fe2+-Fe3+ along c • Jdirect negative (AFM) in CO phase • Jdirect positive (FM) in VM phase

  24. Inelastic neutron scattering • S.Chang etal., PRL 99, 037202 (2007) • J33b = 5.9 meV • J22b = 3.4 meV • J23 = 6.0 meV • J23 = (2J23a + J23c)/3

  25. Theoretical calculations of Jij • Total energy of a certain magnetic configuration given by: ni … number of atoms i zij … number of atoms j which are neighbors of i Si = 5/2 (Fe3+); 2 (Fe2+) si = ±1 • Calculate E-diff when a spin on atom i (Di) or on two atoms i,j (Dij) are flipped • Calculate a series of magnetic configurations and determine Jij by least-squares fit

  26. Investigated magnetic configurations

  27. Calculated exchange parameters

  28. Summary • Standard LDA/GGA methods cannot explain YBaFe2O5 • metallic, no charge order (Fe2+-Fe3+), too small moments • Needs proper description of the Fe 3d electrons (GGA+U, …) • CO-phase: Fe2+: high-spin d6, occupation of a single spin-dn orbital (dxz) • Fe2+/Fe3+ ordered in chains along b, cooperative Jahn-Teller distortion and strong e--lattice coupling which dominates simple Coulomb arguments (checkerboard structure of Fe2+/Fe3+) • VM phase: small orthorhombic distortion (AFM order, moments along b) • Fe d-z2 spin-dn orbital partly occupied (top of the valence bands) leads to direct Fe-Fe exchange across Y-layer and thus to ferromagnetic order (AFM in CO phase). • Quantitative interpretation of the Mössbauer data • Calculated exchange parameters Jij in reasonable agreement with exp.

  29. Thank you for your attention !

  30. tight-binding MO-schemes: too simple? VM phase: Fe2.5+ CO phase: Fe2+ CO phase: Fe3+ Woodward&Karen, Inorganic Chemistry 42, 1121 (2003)

  31. Concepts when solving Schrödingers-equation Treatment of spin Form of potential “Muffin-tin” MT atomic sphere approximation (ASA) pseudopotential (PP) Full potential : FP Non-spinpolarized Spin polarized (with certain magnetic order) Relativistic treatment of the electrons exchange and correlation potential Hartree-Fock (+correlations) Density functional theory (DFT) Local density approximation (LDA) Generalized gradient approximation (GGA) Beyond LDA: e.g. LDA+U non relativistic semi-relativistic fully-relativistic Schrödinger - equation Basis functions Representation of solid plane waves : PW augmented plane waves : APW atomic oribtals. e.g. Slater (STO), Gaussians (GTO), LMTO, numerical basis non periodic (cluster, individual MOs) periodic (unit cell, Blochfunctions, “bandstructure”)

  32. APW Augmented Plane Wave method • The unit cell is partitioned into: • atomic spheres • Interstitial region unit cell Rmt Basisset: PW: ul(r,e) are the numerical solutions of the radial Schrödinger equation in a given spherical potential for a particular energy e AlmKcoefficients for matching the PW join Atomic partial waves

  33. APW based schemes • APW (J.C.Slater 1937) • Non-linear eigenvalue problem • Computationally very demanding • LAPW (O.K.Andersen 1975) • Generalized eigenvalue problem • Full-potential (A. Freeman et al.) • Local orbitals (D.J.Singh 1991) • treatment of semi-core states (avoids ghostbands) • APW+lo (E.Sjöstedt, L.Nordstörm, D.J.Singh 2000) • Efficience of APW + convenience of LAPW • Basis for K.Schwarz, P.Blaha, G.K.H.Madsen, Comp.Phys.Commun.147, 71-76 (2002)

  34. variational methods (L)APW + local orbitals - basis set n…50-100 PWs /atom Trial wave function Variational method: minimum upper bound Generalized eigenvalue problem: H C=E S C Diagonalization of (real symmetric or complex hermitian) matrices ofsize 100 to 50.000 (up to 50 Gb memory)

  35. Quantum mechanics at work

  36. WIEN2k software package An Augmented Plane Wave Plus Local Orbital Program for Calculating Crystal Properties Peter Blaha Karlheinz Schwarz Georg Madsen Dieter Kvasnicka Joachim Luitz November 2001 Vienna, AUSTRIA Vienna University of Technology http://www.wien2k.at WIEN97: ~500 users WIEN2k: ~1150 users mailinglist: 1800 users

  37. Development of WIEN2k • Authors of WIEN2k P. Blaha, K. Schwarz, D. Kvasnicka, G. Madsen and J. Luitz • Other contributions to WIEN2k • C. Ambrosch-Draxl (Univ. Graz, Austria), optics • T. Charpin (Paris), elastic constants • R. Laskowski (Vienna), non-collinear magnetism, parallelization • L. Marks (Northwestern, US) , various optimizations, new mixer • P. Novák and J. Kunes (Prague), LDA+U, SO • B. Olejnik (Vienna), non-linear optics, • C. Persson (Uppsala), irreducible representations • M. Scheffler (Fritz Haber Inst., Berlin), forces • D.J.Singh (NRL, Washington D.C.), local oribtals (LO), APW+lo • E. Sjöstedt and L Nordström (Uppsala, Sweden), APW+lo • J. Sofo and J. Fuhr (Barriloche), Bader analysis • B. Yanchitsky and A. Timoshevskii (Kiev), spacegroup • and many others ….

  38. International co-operations • More than 1000 user groups worldwide • Industries(Canon, Eastman, Exxon, Fuji, A.D.Little, Mitsubishi, Motorola, NEC, Norsk Hydro, Osram, Panasonic, Samsung, Sony). • Europe: (EHT Zürich, MPI Stuttgart, Dresden, FHI Berlin, DESY, TH Aachen, ESRF, Prague, Paris, Chalmers, Cambridge, Oxford) • America:ARG, BZ, CDN, MX, USA (MIT, NIST, Berkeley, Princeton, Harvard, Argonne NL, Los Alamos Nat.Lab., Penn State, Georgia Tech, Lehigh, Chicago, SUNY, UC St.Barbara, Toronto) • far east: AUS, China, India, JPN, Korea, Pakistan, Singapore,Taiwan (Beijing, Tokyo, Osaka, Sendai, Tsukuba, Hong Kong) • Registration at www.wien2k.at • 400/4000 Euro for Universites/Industries • code download via www (with password), updates, bug fixes, news • usersguide, faq-page, mailing-list with help-requests

  39. w2web GUI (graphical user interface) • Structure generator • spacegroup selection • import cif file • step by step initialization • symmetry detection • automatic input generation • SCF calculations • Magnetism (spin-polarization) • Spin-orbit coupling • Forces (automatic geometry optimization) • Guided Tasks • Energy band structure • DOS • Electron density • X-ray spectra • Optics

  40. Program structure of WIEN2k • init_lapw • initialization • symmetry detection (F, I, C-centering, inversion) • input generation with recommended defaults • quality (and computing time) depends on k-mesh and R.Kmax (determines #PW) • run_lapw • scf-cycle • optional with SO and/or LDA+U • different convergence criteria (energy, charge, forces) • save_lapw tic_gga_100k_rk7_vol0 • cp case.struct and clmsum files, • mv case.scf file • rm case.broyd* files

  41. Advantage/disadvantage of WIEN2k + robust all-electron full-potential method (new effective mixer) + unbiased basisset, one convergence parameter (LDA-limit) + all elements of periodic table (equal expensive), metals + LDA, GGA, meta-GGA, LDA+U, spin-orbit + many properties and tools (supercells, symmetry) + w2web (for novice users) ? speed + memory requirements + very efficient basis for large spheres (2 bohr) (Fe: 12Ry, O: 9Ry) - less efficient for small spheres (1 bohr) (O: 25 Ry) - large cells, many atoms (n3, but new iterative diagonalization) - full H, S matrix stored  large memory required + effective dual parallelization (k-points, mpi-fine-grain) + many k-points do not require more memory - no stress tensor - no linear response

  42. Magnetite ferri-magnetic natural mineral, TN=850 K  early “Compass”  proof of earth magnetic field flips Structure below TV ~ 120 K: charge order along (001) planes (1A and 2 B sites) is too simple Mössbauer: 1 A and 4 B Fe sites NMR: 8 A and 16 B sites, Cc symmetry single crystal diffraction, synchrotron diffraction …  small distortions ???  small, but complicated coupling between lattice and charge order

More Related