arxiv: v3 [cond-mat.mtrl-sci] 15 Jul 2016

Similar documents
On the magnetization process of ferromagnetic materials

MAC-solutions of the nonexistent solutions of mathematical physics

Physics 202, Lecture 14

Classical Mechanics. From Molecular to Con/nuum Physics I WS 11/12 Emiliano Ippoli/ October, 2011

Consequently, the temperature must be the same at each point in the cross section at x. Let:

SUMMER KNOWHOW STUDY AND LEARNING CENTRE

DETERMINATION OF MECHANICAL PROPERTIES OF NANOSTRUCTURES WITH COMPLEX CRYSTAL LATTICE USING MOMENT INTERACTION AT MICROSCALE

KINEMATICS OF RIGID BODIES

New Expansion and Infinite Series

4 7x =250; 5 3x =500; Read section 3.3, 3.4 Announcements: Bell Ringer: Use your calculator to solve

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Measuring Electron Work Function in Metal

1 Bending of a beam with a rectangular section

7.2 The Definite Integral

Read section 3.3, 3.4 Announcements:

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

Factors affecting the phonation threshold pressure and frequency

Lesson 1: Quadratic Equations

Review of Calculus, cont d


Math 8 Winter 2015 Applications of Integration

CHM Physical Chemistry I Chapter 1 - Supplementary Material

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

Analytically, vectors will be represented by lowercase bold-face Latin letters, e.g. a, r, q.

13.4 Work done by Constant Forces

Casimir-Polder interaction in the presence of parallel walls

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

5.7 Improper Integrals

Week 10: Line Integrals

Reversible magnetization processes in scalar Preisachtype models of hysteresis

1.1. Linear Constant Coefficient Equations. Remark: A differential equation is an equation

8 Laplace s Method and Local Limit Theorems

Physics 3323, Fall 2016 Problem Set 7 due Oct 14, 2016

Green s functions. f(t) =

1.9 C 2 inner variations

Higher Checklist (Unit 3) Higher Checklist (Unit 3) Vectors

Lecture 1. Functional series. Pointwise and uniform convergence.

Jackson 2.26 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

MATHS NOTES. SUBJECT: Maths LEVEL: Higher TEACHER: Aidan Roantree. The Institute of Education Topics Covered: Powers and Logs

Part I: Basic Concepts of Thermodynamics

Predict Global Earth Temperature using Linier Regression

MASKING OF FERROMAGNETIC ELLIPTICAL SHELL IN TRANSVERSE MAGNETIC FIELD

New data structures to reduce data size and search time

- 5 - TEST 2. This test is on the final sections of this session's syllabus and. should be attempted by all students.

On the Uncertainty of Sensors Based on Magnetic Effects. E. Hristoforou, E. Kayafas, A. Ktena, DM Kepaptsoglou

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams

Intro to Nuclear and Particle Physics (5110)

ARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac

Physics 201 Lab 3: Measurement of Earth s local gravitational field I Data Acquisition and Preliminary Analysis Dr. Timothy C. Black Summer I, 2018

Rates of chemical reactions

Conducting Ellipsoid and Circular Disk

The Wave Equation I. MA 436 Kurt Bryan

p-adic Egyptian Fractions

THE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM

The Properties of Stars

Chapter 1. Chapter 1 1

Thomas Whitham Sixth Form

Physics 9 Fall 2011 Homework 2 - Solutions Friday September 2, 2011

Energy Bands Energy Bands and Band Gap. Phys463.nb Phenomenon

potentials A z, F z TE z Modes We use the e j z z =0 we can simply say that the x dependence of E y (1)

The Regulated and Riemann Integrals

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations

Alg. Sheet (1) Department : Math Form : 3 rd prep. Sheet

This chapter will show you. What you should already know. 1 Write down the value of each of the following. a 5 2

Duality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.

A5682: Introduction to Cosmology Course Notes. 4. Cosmic Dynamics: The Friedmann Equation. = GM s

Math 1B, lecture 4: Error bounds for numerical methods

Vadose Zone Hydrology

Special Relativity solved examples using an Electrical Analog Circuit

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes

6. Intermediate state

Lesson 8. Thermomechanical Measurements for Energy Systems (MENR) Measurements for Mechanical Systems and Production (MMER)

Summary of equations chapters 7. To make current flow you have to push on the charges. For most materials:

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah

CBE 291b - Computation And Optimization For Engineers

Review of basic calculus

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION

The practical version

Recitation 3: More Applications of the Derivative

CHEMICAL KINETICS

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS

than 1. It means in particular that the function is decreasing and approaching the x-

INTERACTION BETWEEN THE NUCLEONS IN THE ATOMIC NUCLEUS. By Nesho Kolev Neshev

DIRECT CURRENT CIRCUITS

Math 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008

CAPACITORS AND DIELECTRICS

Module 6: LINEAR TRANSFORMATIONS

UNIT 1 FUNCTIONS AND THEIR INVERSES Lesson 1.4: Logarithmic Functions as Inverses Instruction

Energy creation in a moving solenoid? Abstract

13: Diffusion in 2 Energy Groups

The International Association for the Properties of Water and Steam. Release on the Ionization Constant of H 2 O

BME 207 Introduction to Biomechanics Spring 2018

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).

Partial Derivatives. Limits. For a single variable function f (x), the limit lim

Quantum Physics III (8.06) Spring 2005 Solution Set 5

R. I. Badran Solid State Physics

MIXED MODELS (Sections ) I) In the unrestricted model, interactions are treated as in the random effects model:

200 points 5 Problems on 4 Pages and 20 Multiple Choice/Short Answer Questions on 5 pages 1 hour, 48 minutes

Transcription:

On the mgnetiztion process in ferromgnetic mterils Ruben Khchturyn nd Vhrm ekhitrin Institute for Physicl Reserch, NA of Armeni, Ashtrk, Armeni (Dted: July 18, 2016) rxiv:1506.01805v [cond-mt.mtrl-sci] 15 Jul 2016 The present rticle concludes tht ferromgnetic smple could be considered s prmgnetic system where roles of mgnetic moments ply mgnetic domins. sed on this conclusion nd tking into ccount presence of n nisotropy field the formul which describes mgnetiztion dependence on the externl mgnetic field is derived. xpressions for remnent mgnetiztion nd coercive force re presented. The new prmeter to chrcterize mgnetic stiffness of mteril is introduced. A physicl expression for dynmic mgnetic susceptibility s function of mterils chrcteristics, externl mgnetic field, nd temperture is given. I. INTRODUCTION A physicl theory permits correctly involve ll interctions in mgnetiztion process nd to revel reltionship between structure nd physicl properties of mgnetic mteril. Applicble mthemticl model could be derived from such theory. This model will give possibility to investigte rel physicl nd structurl properties of mgnetic mterils from experimentl dt. It is essentil for synthesis of new mterils with desired properties. Nowdys there re severl models for describing mgnetiztion of ferromgnetic mterils. ore detiled description nd nlyzes of dvntges nd disdvntges of these models one cn find in works [1]-[4]. In the present pper, there suggests new theory of mgnetiztion nd n ttempt to derive n pplicble generl mthemticl model to describe mgnetiztion curve for soft nd stiff mgnetic mterils. uch generic specil points in mgnetiztion curve re lso elicited from the model. A formul for dynmic mgnetic susceptibility is derived from mgnetiztion field dependence. A new prmeter which cn numericlly chrcterizes the mgnetic stiffness of mteril is introduced nd its physicl interprettion is given. II. AIN IDA Two competing interctions could be distinguished in ferromgnetic mterils: n exchnge interction (exch) which tends to orient mgnetic moments in the sme direction nd by this mgnetizes the system nd dipoledipole interction (dip-dip) which tends to orient mgnetic moments ntiprllel to ech other nd by this demgnetize the system. A relevnt difference between these interctions is tht exch cts between nerest toms nd its energy is independent of mgnetic moment of rubenftf@gmil.com vhrm.mekhitrin@gmil.com system. On the contrry, dip-dip energy rises s mgnetic moment increses. Incresing mteril size dipdip energy cn overcome exch energy. In this cse, two nd more domins structures become more fvorble. It is schemticlly shown in figure 1. As exch energy is much bigger thn dip-dip energy between nerest toms mgnetic moments inside domin re firmly connected in the sme direction. The mgnetic shell where dip-dip energy becomes equl to exch energy could be ccepted s border of the domin nd by this defining size of domin [5]. o, fter domin ws completed the next shell of mgnetic moments will recline in the opposite direction, figure 1. It should be mentioned tht mgnetic moments inside domin re not strictly directed in the sme direction. They recline under n ngle to ech other from shell to shell until the domin would not be finished. The trnsition lyers behviour is detiled described in the work of Lndu-Lifshitz [6]. ecuse of it, mgnetic moment of domin is less thn sum of mgnetic moments of toms inside it. size exch const dipdip ( V ) FIG. 1. Digrm representtion how exch energy nd dip-dip energy behve with incresing of size. On the ssumption of foregoing we could consider domins like solitry mgnetic prticles. The seprted qusiprticles we would cll supermgneton (sm) nlogiclly to R.rrison [7]-[9] where domins re lso supposed to be unit mgnetic moments in the mgnetiztion processes. We know tht exch is compensted by dip-dip between sms. It mens tht sms mgnetic moments re exempt from exch. o, the problem of ferromgnetic mterils is

2 brought to problem of prmgnetic mterils where sms ply the role of mgnetic moments. eeing tht there re xes of esiest mgnetiztion in ferromgnetic mterils [10] sms re distributed in field of nisotropy ccording to the oltzmnn distribution. For simplicity, we would consider cse of unixil nisotropy. ms with positive projection on ny selected direction long the nisotropy xis seprted from sms with negtive projection by nisotropy brrier, figure 2. θ As result of the energetic shift mgnetiztion of the smple is ppered. To estimte the mgnetiztion one needs to clculte difference between mgnetic moments with positive nd negtive projections. ere nd in further positive nd negtive directions would be considered reltive to the direction in which externl mgnetic field is pplied. Tking into ccount tht t ny field vlue distribution of sms is obey to the oltzmnn sttistics we cn clculte difference between mount of smswith negtie nd positive projections: N1 N2 Z g ( + )e d θ +π 0 π π 2 (1) Z θ FIG. 2. ms seprted by nisotropy brrier. g ( 2m)e kt e e N1 + N 2 e d 2m kt 2m2 (2) 2m2 () y the sme wy we get: III. III. AGNTIZATION PROC y pplying mgnetic field to ferromgnetic mteril mgnetic field is induced inside. ms with positive projection on the mgnetic field direction obtin energy -m+, nd sms with negtive projection obtin energy m, where m+ nd m re mgnetic moments of sms with positive nd negtive projection respectively. Tking into ccount tht m+ + m 2m, m is mgnetic moment of sms in zero field, we cn conclude tht m + m+ 2m. It mens tht domins could be effectively replced by sms mgnetic moment of which remins unchnged nd equl to mgnetic moment of domin in nonmgnetized stte. o, potentil wells shift on vlue 2m, s shown on figure. N1 + N 2 e kt Dividing 2 on we get: N1 -N2 1e 2(m+ ) 2(m+ ) (4) 2 (m ) (N1 +N2 ) tnh ultiplying 4 on m nd tking into ccount tht by deffinition: m (N1 + N2 ) (5) m (N1 - N2 ) finlly we obtin: tnh 2 (m ) (6) where is sturtion mgnetiztion nd is mgnetiztion in field. To explore mgnetiztion process more thoroughly it is necessry to understnd field distribution inside mteril during mgnetiztion. As direction of externl field is given we will concentrte on distribution of dipolr field. The mgnetic field which is induced by mgnetic moment in ny point could be clculted in first pproximtion s [12]: FIG.. hift of potentil wells in the induced field. ~ n ) m ~ ~ n (m r (7)

where η is dipolr term, is distnce between nerest domins domins (liner size of sm), η is coefficient which represents difference between positive nd negtive dipolr influences. As n exmple, we will consider two-domin rod s shown on figure 5. Lets compre processes of mgnetiztion of such rod when the externl mgnetic field is directed long longitude (x xis) nd when the externl mgnetic field is directed long width (y or z xes). z FIG. 4. Dipolr field distribution round the mgnetic moment. nd b re lines where mgnetic field chnge their projection sing on mgnetiztion direction. ψ is n ngle between mgnetiztion direction nd lines nd b. 1, re spce regions where dipolr field hs positive projection on mgnetic moment direction nd 2, 4 re spce regions where dipolr field hs negtive projection on mgnetic moment direction. where ~n is unit vector in the direction to the point, nd r is distnce between mgnetic moment nd the point where the field is clculted. It is seen from 7 tht hs negtive or positive projection on the direction of mgnetic moment in different points. It is not difficult to found points where chnges its projection sign on direction of the mgnetiztion from 7. _ ~ _ ~ 0 n n (8) x y FIG. 5. Two domin rod. When the rod is mgnetized long x xis the dipolr field of one domin will mgnetize the second domin becuse it belongs to the region of the spce where dipolr field hs positive projection on mgnetiztion direction. The digrm of the process is shown in figure 6. y this dipolr field will help to mgnetize the system. In this cse s-shped hysteresis loop would be observed [7]-[8]. When the rod is mgnetized long z xis one domin will x 2 ~_ ~_ n0 n ~_ n cos (ψ) (9) cos (ψ) ψ 55 (10) (11) These points belong to the lines nd b which decline under ngles ψ to mgnetiztion direction s shown on figure 4. Lines nd b demrcte spce round mgnetic moment on four regions: 1, 2,, 4. o, hs positive projection in ny point which belongs to regions 1 nd with biggest vlue when or nd hs negtive projection in ny pint of region 2 nd 4 with ~ ~ 2m biggest vlue or 180. r when ψ 0 y this, dipolr field plys both mgnetize (positive) nd demgnetize (negtive) roles. Due to positive influence of depolr field it is possible to mgnetize bulk smples. Thus mgnetic field inside smple could be represented s: η (12) FIG. 6. The rod is mgnetized long longitude. direct the second domin in the opposite direction becuse it belongs to the region of the spce where dipolr field hs negtive projection on mgnetiztion direction s shown in figure 8. In this cse η 1. ecuse of ll written bove, the mgnetic energy term in 6 should be replced by ηm, where η is prmeter which depend on difference between demgnetizing nd mgnetizing prts of dip-dip influences, by this η depends on should depend on surfce nd surfce volume rtio. At ech certin vlues of externl mgnetic field the certin distribution of field inside mteril exists. This field distribution chnges with externl field, nd by this η chnge with externl field s well. m η m + (γ) tnh (1)

4 z R FIG. 7. The rod is mgnetized long width. As eqution 1 is trncendent it could be rewritten in more convenient form s done in [7]-[9]: η + k ( ) T 2m ln + (γ) (14) + m It should be noticed tht the lst term in 14 is independent on the vlue of mgnetic moment s both mgnetic moment of the domins nd nisotropy energy term re both volume dependnt. IV. RANNT AGNTIZATION AND CORCIV FORC After the externl mgnetic field ws bolished dipdip tends to demgnetize the smple. ecuse of this sms would turn from positive projection to negtive through nisotropy brrier. ut not ll sms cn overcome nisotropy brrier nd the prt of them would remin with positive projection. o, fter the externl mgnetic field ws bolished ferromgnetic smple would steel remin in mgnetize condition. This mgnetiztion is clled remnent mgnetiztion. It is represented in figure 9. In order to get expression for remnent mgnetiztion it is necessry to put 0 into 1: [ η m ] R tnh R + (γ) (15) k T where R is remnent mgnetiztion. As is known coercive force is mgnetic field which should be pplied to the smple to demgnetize it. o to get n expression of coercive force one need to put 0 into 14: C (T ) (γ) m (16) It is seen tht coercive force depends on the field direction (ngle between field nd nisotropy xis). V. AGNTIC TIFFN The next prmeter cn chrcterize mgnetic stiffness of ferromgnetic mterils: k e 2 ηm R (γ) k T (17) FIG. 8. sms overcome nisotropy brrier when smple chnge its stte from sturted mgnetiztion to remnent mgnetiztion. Dividing 2 on it is possible to show tht k N2 N 1 nd represents the rtio of the mount of sms which overcme nisotropy brrier to the mount of sms which remin with positive projection fter the externl mgnetic field ws bolished, figure??. Vlues of k could chnge in the rnge from 0 to 1. The bigger the vlue of k the softer mgnetic mteril is nd vice verse. For exmple, in cse of strong nisotropy no of sms re ble to overcome nisotropy brrier, it mens tht there re no sms with negtive projection (k 0) nd R rectngulr-like hysteresis loop, figure 12 ). In the cse when ll sms were ble to overcome nisotropy brrier, the mount of sms with negtive projection is equl to the mount of positive projection (k 1 ) nd consequently R 0, figure 12 c). ) b) FIG. 9. Three possible cse of sms distribution in stte of remnent mgnetiztion: )k 0; b)k 1;; Cses ) nd b) corresponds to mgneto stiff mteril nd mgneto soft mterils correspondingly. VI. DYNAIC AGNTIC UCPTIILITY A low which describes mgnetiztion dependence on externl mgnetic field gives possibility to find out lw for mgnetic susceptibility: χ d [ ] d d m ( d tnh η ) + (γ) k T ( ) m 1 dη d 4π η χ + d (γ) /d [ ] k T ch 2 m(η )+ (γ) k T (18)

χ ( ) [m 1 dη d + d ] (γ) /d k T ch 2 [ m(η )+ (γ) k T ] + η m (19) It should be noted tht η decreses with incresing of thus dη / d < 0. This term is one of fctors which re responsible for high vlue of mgnetic susceptibility of ferromgnetic mterils. It is seen tht mgnetic susceptibility depends on the externl mgnetic field, the dipolr interction, the temperture of the smple nd the nisotropy energy, s it ws expected. In cse of high temperture or smll externl mgnetic field (tnh (x) x; η 1) from 19 one cn get: ( m + d ) (γ) /d χ k T + m / (20) VII. CONCLUION There ws shown tht mgnetic smple could be considered s set of mgnetic prticles, clled superprmgnetons (sms). gnetic moment of sm is equl to mgnetic moment of domin plus mgnetic moment of the domin wll, before externl mgnetic field ws pplied. The strk difference of the sm is tht it s mgnetic moment doesn t chnge during mgnetistion precess nd sms re free from exchnge interction wht gives possibility to pply the oltzmnn sttistics for them s it done for prmgnetic smples. sed on this ssumption there ws derived the nlyticl excretion to describe dependence of mgnetiztion of ferromgnetic mteril on n internl 6. Considering dipolr field distribution inside mteril n nlyticl expression of mgnetiztion dependence on n externl mgnetic field 1 is lso deduced. It is importnt tht ll energies tht tke plce in the process of mgnetiztion re included dditively. It mens tht dditionl energies like energies on domin wlls pinning could be esily dded in the formul in 1. There were derived expressions for remnent mgnetiztion 15 nd coercive force 16 s specil points on ( ) curve. It is seen tht temperture dependence of remnent mgnetiztion bers exponentil chrcter nd depend on dip-dip in smple nd on n nisotropy brrier. It is lso seen how coercive force depend on nisotropy brrier nd mgnetic moment of domin. A new prmeter which chrcterizes mgnetic stiffness of mteril nd its temperture dependence is introduced 17. From ( ) function there ws derived n expression for mgnetic susceptibility 19. It is shown tht in extreml cses, like high temperture or low field, mgnetiztion depends on field linerly, nd mgnetic susceptibility is field independent like in the prmgnetic cse. It is essentil to note tht expressions 19 20 re pplicble t ll temperture regions. 5 [1] D.C Jiles, X. Fng, W. Zhng, ndbook of Advnced gnetic terils, (2006). [2] F.Liorzou,. Phelps, nd D.L. Atherton, croscopic models of mgnetiztion,i trnsctions on mgnetics, vol.6, No.2., (2000). [] ergey. Zirk, Yuriy L. oroz, Robert G. rrison, nd Krzysztof Chwstek, On the physicl spects of the Jiles- Atherton hysteresis models, J. Appl. Phys. 112, (2012). [4] D.C. Jiles nd D.L. Atherton, Theory of ferromgnetic hysteresis, J. gn. gn. ter. 61, 48, (1986). [5] J. Frnkel nd J. Dorfmn; pontneous nd Induced gnetistion in Ferromgnetic odies, (190). [6] L. Lndu,. Lifshits, On the theory of the dispersion of mgnetic permebillity in ferromgnetic bodies, Phys. Zeitsch. der ow. 8, pp. 15-169, (195). [7] R.G. rrison, Physicl model of spin ferromgnetism, (200). [8] R.G. rrison, Vrible-Domin-ize theory of spin ferromgnetism, (2004). [9] R.G. rrison, Physicl theory of ferromgnetic firstorder return curves, I TRANITION ON AG- NTIC, Vol. 45 NO, 4, (2009). [10] L. D. Lndu nd.. Lifshitz, Course of Theoreticl Physics, Vol. 8, (2005). [11] C. Kittel, Theory of the structure of ferromgnetic domins in films nd smll prticles, Phys.Rev. vol.70, NO. 11 nd 12, (1946). [12] L. D. Lndu nd.. Lifshitz, Course of Theoreticl Physics, Vol. 2: The Clssicl Theory of Fields (Nuk, oscow, 1988; Pergmon, Oxford, 1975). [1] Alberto P. Gimres, Principles of Nnomgnetism, pringer, 2009.