Exchange bias and coercivity of ferromagnetic/antiferromagnetic multilayers. Klaus D. Usadel

Similar documents
Monte Carlo simulation study on magnetic hysteresis loop of Co nanowires

Exchange bias in core/shell magnetic nanoparticles: experimental results and numerical simulations

Phenomenology and Models of Exchange Bias in Core /Shell Nanoparticles

Conservation of Angular Momentum = "Spin"

Quantum Mechanics I Problem set No.1

PHYSICS - CLUTCH CH 28: INDUCTION AND INDUCTANCE.

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices

SUPPLEMENTARY INFORMATION

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

A quote of the week (or camel of the week): There is no expedience to which a man will not go to avoid the labor of thinking. Thomas A.

Ising-like dynamics and frozen states in systems of ultrafine magnetic particles

Susceptibility and Inverted Hysteresis Loop of Prussian Blue Analogs with Orthorhombic Structure

Spin-rotation coupling of the angularly accelerated rigid body

A Hybrid Variational Iteration Method for Blasius Equation

Monte Carlo study of the exchange bias effects in magnetic nanoparticles with core/shell morphology

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

11. Dynamics in Rotating Frames of Reference

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2

Oguchi Approximation of a mixed spin-2 and spin-5/2 Blume-Capel Ising ferrimagnetic system

PHYSICS - CLUTCH 1E CH 28: INDUCTION AND INDUCTANCE.

SUPPLEMENTARY INFORMATION

The exchange interaction between FM and AFM materials

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.

Dipolar Interaction Effects in the Magnetic and Magnetotransport Properties of Ordered Nanoparticle Arrays

Parametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

PHYS 1101 Practice problem set 12, Chapter 32: 21, 22, 24, 57, 61, 83 Chapter 33: 7, 12, 32, 38, 44, 49, 76

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Lecture 20: Noether s Theorem

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Mass Transfer Processes

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Physics 207 Lecture 6

Relaxation in water /spin ice models

8 Derivation of Network Rate Equations from Single- Cell Conductance Equations

Physics 240: Worksheet 30 Name:

PHYS 215C: Quantum Mechanics (Spring 2017) Problem Set 3 Solutions

Quantum spin system with on-site exchange in a magnetic field

Physics 181. Particle Systems

Thermo-Calc Software. Modelling Multicomponent Precipitation Kinetics with CALPHAD-Based Tools. EUROMAT2013, September 8-13, 2013 Sevilla, Spain

OFF-AXIS MECHANICAL PROPERTIES OF FRP COMPOSITES

CHAPTER II THEORETICAL BACKGROUND

Classical Mechanics ( Particles and Biparticles )

Lab 2e Thermal System Response and Effective Heat Transfer Coefficient

Supplemental Material: Causal Entropic Forces

Tensor Smooth Length for SPH Modelling of High Speed Impact

Chapter Newton s Method

Lecture 3. Interaction of radiation with surfaces. Upcoming classes

Thermodynamics General

...Thermodynamics. If Clausius Clapeyron fails. l T (v 2 v 1 ) = 0/0 Second order phase transition ( S, v = 0)

Mathematical Preparations

and Statistical Mechanics Material Properties

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Uniqueness of Weak Solutions to the 3D Ginzburg- Landau Model for Superconductivity

Non-interacting Spin-1/2 Particles in Non-commuting External Magnetic Fields

Physics 111: Mechanics Lecture 11

8 Derivation of Network Rate Equations from Single- Cell Conductance Equations

CHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)

Why BP Works STAT 232B

Computing π with Bouncing Balls

2 Finite difference basics

Resource Allocation with a Budget Constraint for Computing Independent Tasks in the Cloud

8.1 Arc Length. What is the length of a curve? How can we approximate it? We could do it following the pattern we ve used before

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force.

Physical Chemistry I for Biochemists. Chem340. Lecture 16 (2/18/11)

The Geometry of Logit and Probit

Magnon pairing in quantum spin nematic

Lecture Note 3. Eshelby s Inclusion II

Wilbur and Ague 4 WILBUR AND AGUE; APPENDIX DR1. Two-dimensional chemical maps as well as chemical profiles were done at 15 kv using

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

APPENDIX F A DISPLACEMENT-BASED BEAM ELEMENT WITH SHEAR DEFORMATIONS. Never use a Cubic Function Approximation for a Non-Prismatic Beam

Supplementary Materials for

V.C The Niemeijer van Leeuwen Cumulant Approximation

FTCS Solution to the Heat Equation

Percolation and the Potts Model

Relaxation laws in classical and quantum long-range lattices

Chapter 4 The Wave Equation

A Solution of the Harry-Dym Equation Using Lattice-Boltzmannn and a Solitary Wave Methods

Chapter 3. r r. Position, Velocity, and Acceleration Revisited

Errors for Linear Systems

Week 9 Chapter 10 Section 1-5

Solutions Review Worksheet

ERROR RESEARCH ON A HEPA FILTER MEDIA TESTING SYSTEM OF MPPS(MOST PENETRATION PARTICLE SIZE) EFFICIENCY

DC Circuits. Crossing the emf in this direction +ΔV

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015

International Journal of Mathematics Trends and Technology (IJMTT) Volume 48 Number 2 August 2017

Effect of loading frequency on the settlement of granular layer

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018

modeling of equilibrium and dynamic multi-component adsorption in a two-layered fixed bed for purification of hydrogen from methane reforming products

Investigation of a New Monte Carlo Method for the Transitional Gas Flow

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model

Chapter Eight. Review and Summary. Two methods in solid mechanics ---- vectorial methods and energy methods or variational methods

Appendix B. The Finite Difference Scheme

Brownian-Dynamics Simulation of Colloidal Suspensions with Kob-Andersen Type Lennard-Jones Potentials 1

Perfect Fluid Cosmological Model in the Frame Work Lyra s Manifold

Lecture 21: Numerical methods for pricing American type derivatives

Transcription:

Exchange bas and coercvty of ferromagnetc/antferromagnetc multlayers Klaus D. Usadel Unverstät Dusburg-Essen, Germany Contents: exchange bas: ntroducton doman state model: results from MC smulatons mean feld approach for Isng AFM: coercvty and exchange bas generalsaton to local vector spn models: dynamcal consequences conclusons Coworkers: U. Nowak, A. Msra, B. Beckmann, R. Stamps

Exchange bas exchange bas (EB) = shfted hysteress loop frst observaton n Co/CoO partcles (Meklejohn and Bean, Phys. Rev. 12,1413 (1956)) typcal for FM/AFM compounds lke multlayers, nanopartcles ntal procedure: coolng the system n an external feld from above to below Neel temperature of the AFM loop s shfted upwards and asymmetrc; enhancement of the coercvty B Hyst. FC T N T

Applcatons: sensors 1.5 F(Sx) magnetzaton -.5 H 2 < H 1 H = 1 < -1 1 S x -1-1 -.5.5 1 external feld F(S x ) = DSx 2 BS x The felds B and B + at whch the magnetzaton of the FM swtches are obtaned from F = at S x = 1 and S x = 1, respectvely: B = 2D, B + = 2D.

Understandng EB: problems external feld and exchange feld of the FM nterface layer polarzes the AFM nterface layer why s the polarsaton stable durng a hysteress cycle? why s the symmetry broken? uncompensated compensated perfect nterface rough nterface Doman wall formaton

Soluton: doman walls n the AFM breakng the symmetry necessary for EB by ntroducng frozen domans perpendcular doman walls n the AFM due to nterface roughness (Malozemoff, Phys. Rev. B 35, 3679 (1987)) but: doman wall formaton unlkely for large AFM thcknesses dea: ntroducng defects n the AFM to stablze domans: doman state model (Mlteny, Gerlngs, Keller, Beschoten, Güntherodt, Nowak, and Usadel, PRL 84, 4224 (2))

Experments wth dluted AFM Idea: generate defects n AFM: CoO Co 1 x Mg x O nterface layer wthout defects vary bulk dluton FM: Co AFM: CoO CoMgO bulk dluton enhances exchange bas(eb) assocated wth EB s an enhancement of the coercvty see also: Mewes et al., APL 76, 157 (2) Sh et al., JAP 91, 7763 (22) Mlteny, Gerlngs, Keller, Beschoten, Güntherodt, Nowak, Usadel, PRL 84, 4224 (2)

Local spn model FM M AFM z y EB x B H = H F + H AF + H nt H F = J F S S j,j H AF = J AF ǫ ǫ j σ σ j,j H nt = J nt ( ) DS 2 x + S x B ǫ S x σ, ǫ σ B D > ; J AF = J F /2; J nt = +/ J AF ; Monte Carlo smulaton, system sze up to 128 128 (9 + 1) up to 136 MCS per hysteress, average over up to 1 runs local mean feld theory (nt)

Hysteress of FM and AFM nterface: MC smulatons system cooled n a feld B c =.25J F down to 1.5 B + k B T =.1J F AFM frozen n a doman state wth nterface magnetzaton m (FM) -.5-1 B_ ts rreversble part leads to exchange bas B EB.3J nt D =.1J F m (AFM Interface).1 -.1 m I -.3 -.2 -.1.1.2.3 B

Structure of the AFM domans B FC Hyst. T N T Staggered AFM magnetzaton after feld coolng above: large dluton, p =.5, small domans, below: small dluton, p =.3, larger domans.

Magnetzaton reversal of the FM layer 1 θ = 8 4 Magnetzaton my -1 3 2 1 coherent rotaton of the FM magnetzaton no asymmetry for θ local FM spns S = S replaced by a macro spn -1 1 Magnetzaton m z F(S x ) = NlDSx 2 NlBS x k B T Tr e β(h AF+H nt ) F (S x ) = 2NlDS x NlB J nt σ F (S x ) = 2NlD βj 2 nt nt ( ) 2 (σ σ ) nt m nt = nt σ l number of FM monolayers N spns per FM monolayer

Ñ ¼ ½µ Â ÒØ Ñ ¼ Ñ ¼ ½µ Â ÒØ Swtchng felds ½ F(Sx) À H 2 < H 1 H = 1 < -1 1 S x ½ F (S x ) = 2NlDS x NlB J nt m nt Felds B and B + for magnetzaton swtchng from F = at S x = 1 and S x = 1, respectvely: À B = 2D J nt m nt ( B, S x = 1 ) /l, B + = 2D J nt m nt ( B+, S x = 1 ) /l, m nt = m + m rev m rev = χ (1) AF J nts x + χ (2) AF B

Hysteress of FM and AFM nterface: MC smulatons system cooled n a feld B c =.25J F down to 1.5 B + k B T =.1J F AFM frozen n a doman state wth nterface magnetzaton m (FM) -.5-1 B_ ts rreversble part leads to exchange bas B EB.3J nt D =.1J F m (AFM Interface).1 -.1 m I -.3 -.2 -.1.1.2.3 B

Lnear approxmaton m nt (B, S x = ±1): AFM nterface magnetzaton determnes coercvty and EB B ± = ±2D J nt m nt (B ±, S x = 1)/l m nt = m + m rev m rev = χ (1) AF J nts x + χ (2) AF B B ± = ±2D J ntm /l±j 2 nt χ(1) AF /l 1+J nt χ (2) AF /l B eb = 1 2 (B + + B ) = J ntm /l 1+J nt χ (2) AF /l B c = 1 2 (B + B ) = 2D+J2 nt χ(1) AF /l 1+J nt χ (2) AF /l dependence on the sgn of J INT H c /J AF.14.1.6.2 H c ln H c num.5 1 1.5 2 (1) T/T N mean feld approach

MF approxmaton; no dluton σ = m = tanh (β ( J AF j m j + J nt S x + B )).18 (a).18 (b).14.14 H c /J AF.1 H c /J AF.1.6.2 L=1 L=2 L=3 L=8.6.2 L=1 L=2 L=3 L=8.5 1 1.5 2 (1) T/T N.5 1 1.5 2 (1) T/T N Coercve feld as functon of reduced temperature for dfferent AFM layer thcknesses l. left: J nt = J AF rght: J nt = J AF D/J F =.5 large ncrease of the coercvty n the vcnty of the Neel temperature T N (Scholten, Usadel, Nowak, Phys. Rev. B 71, 64413 (25))

Dluted systems Local mean feld equatons: m = ǫ tanh (β ( J AF j ǫ j m j + J nt S x + B )) Coolng: Iteraton at a fxed temperature untl a (metastable) self-consstent soluton s obtaned, then reducng temperature n small steps.

Coercvty and bas felds: dfferent AFM dluton p.18.18.14.14 H c /J AF.1 H c /J AF.1.6.6.2.2.5.5 H eb /J AF.4.3.2 p=. p=.1 p=.3 p=.6 H eb /J AF.4.3.2 p=. p=.1 p=.3 p=.6.1.1.5 1 1.5 2 (1,) T/T N l = 3 and J nt = J AF.5 1 1.5 2 (1,) T/T N l = 3 and J Int = J AF maxmum of the coercvty ndependent of dluton: χ max = x/t N (x); x concentraton of magnetc stes

Generalzaton to AFM vector spns H o = BS x () D S x () 2 + H ex h α () = S α () H o + J nt σ α () h α () = S α () H o + J nt m α () m α () = χ (1) αβ B β + J nt χ (2) αβ S β() + m,α () Effectve (free) energy of the FM layer: separaton of tme scales and/or slow varaton n space: thermal average restrcted to the AFM AFM eqlbrum susceptbltes χ (1) αβ and χ (2) αβ and exchange feld as response to external feld F = H o J nt S α ()χ (1) αβ B β 1 2 J2 nt S α ()χ (2) αβ S β() S α ()m,α () enhanced moment; enhanced ansotropy

F = H o J nt S α ()χ (1) αβ B β 1 2 J2 nt S α ()χ (2) αβ S β() S α ()m,α () maxmum lowerng of symmetry occurs when only one component of the susceptblty tensor, χ x,x, s nonzero: Isng antferromagnet H o = BS x () DS 2 x() + H ex S α ()m,α () B = B[1 + (J nt /l)χ] D = D + [Jnt/(2l)]χ 2 = B c /2 strong dependence on temperature relatvely weak ansotropy n the ferromagnet: maxmum value of h c /J nt.1 for 2D/J nt =.2. Ths corresponds to ( B/B) max 1.1 and ( D/D) max 5. Dynamcal consequences: doman wall wdth and doman wall energy E DW : ( ) max (ẼDW) max J/ D; J D wall velocty v DW : ṽ DW B. Stamps, Usadel,Europhys. Lett., 74, 512 (26)

Conclusons Doman state model explans exchange bas and many effects assocated wth t wthout explctely assumng some net AFM nterface magnetzaton frozen AFM nterface magnetzaton leads to EB reversble part of the AFM nterface magnetzaton leads to enhanced coercvty provdng the AFM s ansotropc for slow FM dynamcs and slow spatal varaton of ts magnetzaton an effectve FM energy can be obtaned after ntegratng out the AFM degrees of freedom