Self-consistent approach for magnetic ordering and excited site occupation processes in a two-level system

Ferromagnetic ordering in a two-level partially excited system is studied in details. Magnitudes of magnetization (magnetic order parameter) and lattice ordering (excited level occupation number) are calculated self-consistently. The influence of an external magnetic field and excited level gap on t...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Datum:2015
Hauptverfasser: Korynevskii, M.A., Solovyan, V.B.
Format: Artikel
Sprache:English
Veröffentlicht: Інститут фізики конденсованих систем НАН України 2015
Schriftenreihe:Condensed Matter Physics
Online Zugang:http://dspace.nbuv.gov.ua/handle/123456789/155267
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Zitieren:Self-consistent approach for magnetic ordering and excited site occupation processes in a two-level system / M.A. Korynevskii, V.B. Solovyan // Condensed Matter Physics. — 2015. — Т. 18, № 3. — С. 33706: 1–12. — Бібліогр.: 34 назв. — англ.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-155267
record_format dspace
spelling irk-123456789-1552672019-06-17T01:28:31Z Self-consistent approach for magnetic ordering and excited site occupation processes in a two-level system Korynevskii, M.A. Solovyan, V.B. Ferromagnetic ordering in a two-level partially excited system is studied in details. Magnitudes of magnetization (magnetic order parameter) and lattice ordering (excited level occupation number) are calculated self-consistently. The influence of an external magnetic field and excited level gap on the ferromagnetic phase transition is discussed. Детально вивчено процес феромагнiтного впорядкування у дворiвневiй частково збудженiй системi. Величина намагнiченостi (магнiтний параметр впорядкування) та граткове впорядкування (заповнення збуджених рiвнiв) розрахованi самоузгоджено. Обговорено вплив зовнiшнього магнiтного поля i ширини енергетичної щiлини мiж основним i збудженим рiвнями на феромагнiтний фазовий перехiд. 2015 Article Self-consistent approach for magnetic ordering and excited site occupation processes in a two-level system / M.A. Korynevskii, V.B. Solovyan // Condensed Matter Physics. — 2015. — Т. 18, № 3. — С. 33706: 1–12. — Бібліогр.: 34 назв. — англ. 1607-324X DOI:10.5488/CMP.18.33706 arXiv:1510.06907 PACS: 75.10.Hk, 75.30.Cr http://dspace.nbuv.gov.ua/handle/123456789/155267 en Condensed Matter Physics Інститут фізики конденсованих систем НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language English
description Ferromagnetic ordering in a two-level partially excited system is studied in details. Magnitudes of magnetization (magnetic order parameter) and lattice ordering (excited level occupation number) are calculated self-consistently. The influence of an external magnetic field and excited level gap on the ferromagnetic phase transition is discussed.
format Article
author Korynevskii, M.A.
Solovyan, V.B.
spellingShingle Korynevskii, M.A.
Solovyan, V.B.
Self-consistent approach for magnetic ordering and excited site occupation processes in a two-level system
Condensed Matter Physics
author_facet Korynevskii, M.A.
Solovyan, V.B.
author_sort Korynevskii, M.A.
title Self-consistent approach for magnetic ordering and excited site occupation processes in a two-level system
title_short Self-consistent approach for magnetic ordering and excited site occupation processes in a two-level system
title_full Self-consistent approach for magnetic ordering and excited site occupation processes in a two-level system
title_fullStr Self-consistent approach for magnetic ordering and excited site occupation processes in a two-level system
title_full_unstemmed Self-consistent approach for magnetic ordering and excited site occupation processes in a two-level system
title_sort self-consistent approach for magnetic ordering and excited site occupation processes in a two-level system
publisher Інститут фізики конденсованих систем НАН України
publishDate 2015
url http://dspace.nbuv.gov.ua/handle/123456789/155267
citation_txt Self-consistent approach for magnetic ordering and excited site occupation processes in a two-level system / M.A. Korynevskii, V.B. Solovyan // Condensed Matter Physics. — 2015. — Т. 18, № 3. — С. 33706: 1–12. — Бібліогр.: 34 назв. — англ.
series Condensed Matter Physics
work_keys_str_mv AT korynevskiima selfconsistentapproachformagneticorderingandexcitedsiteoccupationprocessesinatwolevelsystem
AT solovyanvb selfconsistentapproachformagneticorderingandexcitedsiteoccupationprocessesinatwolevelsystem
first_indexed 2025-07-14T07:19:33Z
last_indexed 2025-07-14T07:19:33Z
_version_ 1837606645138456576
fulltext Condensed Matter Physics, 2015, Vol. 18, No 3, 33706: 1–12 DOI: 10.5488/CMP.18.33706 http://www.icmp.lviv.ua/journal Self-consistent approach to magnetic ordering and excited site occupation processes in a two-level system M.A. Korynevskii1,2,3, V.B. Solovyan1 1 Institute for Condensed Matter Physics of the National Academy of Sciences of Ukraine, 1 Svientsitskii St., 79011 Lviv, Ukraine 2 Lviv Polytechnic National University, 12 Bandera St., 79013 Lviv, Ukraine 3 Institute of Physics of the University of Szczecin, 15 Wielkopolska St., 70451 Szczecin, Poland Received March 26, 2015, in final form July 1, 2015 Ferromagnetic ordering in a two-level partially excited system is studied in detail. Magnitudes of magnetiza- tion (magnetic order parameter) and lattice ordering (excited level occupation number) are calculated self- consistently. The influence of an external magnetic field and excited level gap on the ferromagnetic phase transition is discussed. Key words:magnetics, two-level system, phase transitions PACS: 75.10.Hk, 75.30.Cr 1. Introduction In numerous instances of real physical systems, the problem of ferromagnetic (or ferroelectric) or- dering strongly depends on the density of magnetic (or dipole) particles in a medium. Such a situation is observed in different types of materials, namely in solids and liquids. Particularly for solids, the irreg- ularity in the position of magnetic particles is determined by a type of the sample preparation (slow or rapid cooling), and requires the use of equilibrium or nonequilibriummethods of statistics to investigate them. On the other hand, the irregularity of a system can be created by the excitation radiation or by heat- ing in such a way that only a part of the total number of particles moves to the overlying energy levels. In case of constant external radiation (temperature), such a system is equilibrium and stable. The degree of occupation of the ground and excited states depends on their mutual energy distance and the intensity of radiation (thermostat temperature). The exchange interaction between the magnetic particles (in the ground or excited states) enables magnetic ordering processes in a system. Their intensity essentially de- pends on the correlations between ground-ground, ground-excited and excited-excited states of particles. Both an external magnetic field and radiation (thermostat temperature) form self-consistent thermo- dynamic states of a system with a certain magnetization and an excited state occupation. For physically- critical applications, the intensity of external radiation (pump) or thermostat temperature may be con- sidered as a constant. A two-level model system can be used to investigate such a system and calculate its characteristics. Two-level systems are relatively simple, they are generally well studied (see [1–4]) and form the basis for the microscopic mechanisms description of different physical phenomena and processes, such as non-linear optics, lasers, Josephson junctions, Kondo scattering, glasses, etc. [5–9]. We shall regard the capability of spin particles of occupying one of the two levels: the ground and excited one with a certain probability dependent on the temperature and interaction between particles. A special attention is paid ©M.A. Korynevskii, V.B. Solovyan, 2015 33706-1 http://dx.doi.org/10.5488/CMP.18.33706 http://www.icmp.lviv.ua/journal M.A. Korynevskii, V.B. Solovyan herein to the analysis of a correlation between spin orientation and occupation processes, as well as to the role of the external magnetic field and interlevel distance in the stable state of a system formation. 2. Hamiltonian and free energy We consider the N -site lattice system of spin-like particles in each site which is capable of occupying one of the two quantum states, i.e., ground state and the excited one. The Hamiltonian of such a system in the external field γ, considering the pair exchange interactions between the particles, can be presented in the form: Ĥ =−γ N∑ i=1 Ŝi − 1 2 N∑ i, j=1 2∑ λ,ν=1 Jλνi j Ĉ (λ) i Ĉ (ν) j Ŝi Ŝ j + N∑ i=1 2∑ λ=1 ελĈ (λ) i . (2.1) Here, Ŝi = Ŝz i is z Pauli matrix; Ĉ (λ) i is the operator of the number of particles in the i -th site in the quantum state λ; Jλνi j is an integral of the exchange interaction between spin particles in the i -th and j -th sites of the crystalline lattice; ελ is a configuration part of the energy per one particle in the λ state. Since the total number of particles coincides with N , N∑ i=1 2∑ λ=1 Ĉ (λ) i = N Î , (2.2) Ĉ (1) i = ( 1 0 0 0 ) , Ĉ (2) i = ( 0 0 0 1 ) , Î = ( 1 0 0 1 ) , Ŝi = ( 1 0 0 −1 ) , we can pass from two operators Ĉ (1) i and Ĉ (2) i to the unique operator Ĉi = Ĉ (2) i , which coincides with the excited particle number operator. The eigenvalues of operator Ĉi are: Ci = { 0, when a particle in the i -th site is in the ground state, 1, when a particle in the i -th site is in the excited state. (2.3) The orthonormal set of the wave functions for each site of the lattice is of the form of the product of configuration and spin components in the following way:{ ψ(λ,s) i } ≡ { ϕ(λ) i } · { χ(s) i } , (2.4) ϕ(1) i = ( 1 0 ) c , ϕ(2) i = ( 0 1 ) c ; χ(1) i = ( 1 0 ) s , χ(−1) i = ( 0 1 ) s . Thus, in reality we have a four-state situation. It is well-known [10–12] that for atoms, thewave functions of excited electron states aremuch broader than the other ones. As a result, the exchange integrals J 11 i j between a pair of the nearest electrons being in the ground state (similarly between a particle in the ground state and the second one in the excited state J 12 i j ) can be regarded as small compared with the exchange integral J 22 i j between both particles in excited states. As a zero approximation for exchange interaction in expression (2.1), we take J 11 i j and J 12 i j negligi- bly small compared with J 22 i j . This is not a principal limitation, but makes it possible to use only one interaction parameter J 22 i j . Such an approach is very close to Vonsovskii s–d exchange model in the the- ory of magnetism [13, 14], according to which s-electrons are responsible for electric conductivity while d -electrons form magnetic properties of solids. The comparison of exchange integrals for ground and excited states for rare-earth metals can be also found in [15]. Following the above mentioned restrictions and taking the configurational part of energy of particles in the ground state equal to zero, the model Hamiltonian (2.1) takes the form: Ĥ =−γ N∑ i=1 Ŝi − 1 2 N∑ i, j=1 J 22 i j Ĉi Ĉ j Ŝi Ŝ j +ε0 N∑ i=1 Ĉi . (2.5) 33706-2 Self-consistent approach to magnetic ordering Here, ε0 ≡ ε2 − ε1 plays the role of energy gap between excited and ground states for a particle in each site, J 22 i j > 0, which corresponds to ferromagnetic exchange coupling. By certain features, the model described by the Hamiltonian (2.5) is close to a well-known Blume- Emery-Griffits (BEG) model [16] introduced to simulate the thermodynamic behaviour of 3 He- 4 He mix- tures and to numerous modifications of BEG model used for interpretation of diluted magnets properties [17–23]. The papers in this field of investigation are generally based on the use of the Ising Hamiltonian with a spin equal to unity as well as on some assumption concerning the character of randomness for magnetic particles in a lattice. Briefly speaking, those models are three-state site models due to the three possible projections of Sz -operator. Quite different physical situation is observed in our investigation, namely a physical interpretation of expression (2.5) comes to two-level, four-state site model (two spin- orientations for the particle in the ground state and two spin-orientations for this particle in the excited state). Temperature helps to increase the excited state occupation but destroys spontaneous magnetiza- tion caused by the external field and the exchange interaction between spin particles in the excited states. Such an approach, in our opinion, is new and quite interesting. Moreover, both quantities c and s being physically different are calculated here self-consistently. Themethod of self-consistent field approximation is used to calculate statistical and thermodynamical properties of the system described by the Hamiltonian (2.5). From the physical point of view, this approx- imation corresponds to the neglect of square fluctuations of the identical physical values for a pair of different sites in a crystalline lattice, the same being true for different physical values in the common site. It is well-known that this approximation is valid in a wide range of temperatures and external fields except for the immediate phase transition point neighbourhood [24, 25]. Thus, we accept the following expressions for the products of the operators: Ĉi Ĉ j = 〈Ĉi 〉Ĉ j + Ĉi 〈Ĉ j 〉−〈Ĉi 〉〈Ĉ j 〉, Ŝi Ŝ j = 〈Ŝi 〉Ŝ j + Ŝi 〈Ŝ j 〉−〈Ŝi 〉〈Ŝ j 〉, (2.6) Ĉi Ŝi = 〈Ĉi 〉Ŝi + Ĉi 〈Ŝi 〉−〈Ĉi 〉〈Ŝi 〉, where 〈Ĉi 〉 ≡ c and 〈Ŝi 〉 ≡ s are mean values of Ĉi and Ŝi operators, correspondingly, according to the Gibbs distribution based on the Hamiltonian (2.5) and expressions (2.6). Being limited in the summation over the sites i , j only to the nearest neighbours ( ∑ j J 22 i j = Xv J , Xv is the number of the nearest neighbours, J is an interaction constant), we obtain for Ĥ the following expression: Ĥ = N∑ i=1 Ĥi = N Ĥi , Ĥi =−γŜi − csXv J ( cŜi + sĈi )+ε0Ĉi + 3 2 c2s2Xv J . (2.7) Thus, the Hamiltonian (2.7) describes the N -particle system with one spin particle in each site. At the same time, each particle staying at one of the two energy levels can take up one of the four states. The total number of particles in a system is constant, so canonical Gibbs ensemble may be used. Contrary to the well-known situations for the instance of independent subsystems: 1) pure or diluted magnetic system (c = const), or 2) simple lattice-gas system (s = const), the Hamiltonian (2.7) describes both binding subsystems, and the values of c and s are in a close connection depending on temperature T , external field γ and energy gap ε0. A partition function of the system described by the Hamiltonian (2.7) is presented as: Z = Tr[exp(−βĤ) ]= { Tr [ exp(−βĤi ) ]}N . (2.8) Here, β= (kBT )−1 , kB is the Boltzmann constant, T is the absolute temperature, and thus the Tr operation should be performed based on the wave functions (2.4). After performing this operation we obtain the following expression for the free energy F per one site: F =− 1 β ln { 2cosh [ β ( c2sXv J +γ)]}− 1 β ln { 1+exp [ β(cs2Xv J −ε0) ]}+ 3 2 c2s2Xv J . (2.9) 33706-3 M.A. Korynevskii, V.B. Solovyan 3. Order parameters c and s The free energy in the form of (2.9) is the function of three independent parameters: temperature T , external field γ and excited energy gap ε0. Using a self-consistent field approximation, we introduced two order parameters c and s, which as the functions of T , γ, ε0 must satisfy the conditions of thermodynamic stability in the form [26, 27]: ∂F ∂ai = 0, i = 1,2,3, . . . , (3.1) where a1 = c and a2 = s in our instance. On the other hand, a system entropy S and magnetization per one site s can be calculated as partial derivatives: S =−∂F ∂T , (3.2) s =−∂F ∂γ . (3.3) Undoubtedly, equation (3.1) for i = 2 and equation (3.3) coincide. Based on the equations (3.1), we get: 3cs −2c tanh [ β ( c2sXv J +γ)]− s { 1+exp [−β(cs2Xv J −ε0) ]}−1 = 0, 3cs − c tanh [ β ( c2sXv J +γ)]−2s { 1+exp [−β(cs2Xv J −ε0) ]}−1 = 0. (3.4) Combining these two equations we obtain for c and s the system of equations in a more convenient form: c = 1 2 { 1+ tanh [ β 2 ( cs2Xv J −ε0 )]} , s = tanh [ β ( c2sXv J +γ)] . (3.5) It can be shown that for arbitrary values of Xv J , the following relationship between c and s takes place: exp [ 2β(cε0 + sγ) ]= ( 1+ s 1− s )s ( 1− c c )2c , (3.6) which at γ= 0 and ε0 = 0 transforms into a simpler form:( 1− s 1+ s )s = ( 1− c c )2c . (3.7) Equations (3.6), (3.7) show an unambiguous relation between c and s for different values of external parameters T , γ, ε0. We start the analysis of equations (3.5) from the situation where the occupation process andmagnetic ordering are completely independent. In other words, the occupation of the excited level c is determined by temperature only (the effect of magnetic ordering is missing, s = 0). In this instance, according to (3.5), we obtain: c = [ 1+exp ( βε0 )]−1 , (3.8) and, correspondingly, s = tanh [ β ( c2 0 sXv J +γ)] , (3.9) where c0 is independent of temperature and has got a fixed (for example by external radiation) con- centration of particles in the excited states. The reduced temperature t = T (Xv J )−1 dependencies of c [figure 1 (a)] and s [figure 1 (b)] for different values of the reduced energy gap ε = ε0(Xv J )−1 and the reduced external field h = γ(Xv J )−1 are presented in figure 1. One can observe a smooth behaviour of both c and s for any h, ε. The curves for magnetization are characteristic of the instance of a second order phase transition. 33706-4 Self-consistent approach to magnetic ordering (a) (b) Figure 1. Temperature dependencies of independent occupation number [figure 1 (a), curves: 1 — ε = −0.5, 2 — ε = −0.25, 3 — ε = −0.1, 4 — ε = 0, 5 — ε = 0.25, 6 — ε = 0.45, 7 — ε = 0.55] and magnetic order parameter [figure 1 (b), curves: 1 — h = 0, 2 — h = 0.07, 3 — h = 0.115. . . , 4 — h = 0.192. . . , 5 — h = 0.318. . . , 6— h = 0.5, (c0 = 1)]. Quite different is the situation when considering bounded system equations (3.5) (see figure 2). Now the correlation between two order parameters (c and s) becomes significant and the role of “external” independent factors (h,ε) is non-trivial. For each value of energy gap ε < 0.5, the curves for c and s demonstrate a smooth temperature be- haviour for h > hc (hc dependent on ε). At h < hc, c and s show jumps at a certain temperature. The values of those jumps decrease with an increase of h. At h = hc, this jump turns into zero [curves 3 in figures 2 (a),(e), curves 4 in figures 2 (b),(f), curves 5 in figures 2 (c),(g)]. Consequently, for h < hc, the first order magnetic phase transition takes place in the investigated system, at h = hc the order of phase transition changes into a second order (tricritical point) and for h > hc no magnetic phase transition is possible. The point hc divides the external fields into “weak” and “strong” ones [27]. Those results are presented in figures 2 (a)–(c) and figures 2 (e)–(g). For ε > 0.5 instance, only smooth curves c and s are observed. The role of the external field h is quantitative, no critical hc was found [figures 2 (d),(h)] and spontaneous magnetization is impossible. A special role of ε parameter is connected with a strong correlation between the processes of the excited level occupation and spontaneous ordering of magnetic particles at this level. Heating facilitates the filling of the excited level while the orderedmagnetic state decreases. The measure of those processes is the value of the system energy, i.e., the ratio between its growth in the excitation and a decrease due to a magnetic ordering. We observe a situation close to the phenomenon of percolation in the systems having a constant concentration of magnetic particles [28]. However, the phenomenon of percolation in the system studied occurs self-consistently where the excited level occupation by magnetic particles depends on the degree of magnetization. That is why for the first order magnetic phase transition, a sharp decrease (dip) of the excited level occupation is observed above the phase transition temperature point Tc [see figures 2 (a)–(d)]. Taking into account that magnetization of the system is represented by the following formula: ~M =~µN s, (3.10) where ~µ is a magnetic moment of an individual particle, based on the equations (3.5), we have obtained the expression for generalized magnetic susceptibility (as a function of temperature and external mag- netic field): χ = |~µ|N Y1 Xv J (Y1Y2 −4c2s2) , (3.11) Y1 = 1 βXv J (1− c)c − s2, (3.12) Y2 = 1 βXv J (1− s2) − c2. 33706-5 M.A. Korynevskii, V.B. Solovyan (a) (e) (b) (f) (c) (g) (d) (h) Figure 2. Temperature dependencies for occupation number [figures 2 (a)–(d)] and magnetic order pa- rameter [figures 2 (e)–(h)] at different energy gaps ε and external fields h (curves: 1— h = 0, 2— h = 0.07, 3— h = 0.115. . . , 4— h = 0.192. . . , 5— h = 0.318. . . , 6— h = 0.50). 33706-6 Self-consistent approach to magnetic ordering (a) (b) (c) (d) Figure 3. Temperature dependencies of generalized magnetic susceptibility at different energy gaps ε and external fields h (numbers of curves correspond to the same values of h as in figure 2). Dependencies of magnetic susceptibility χ (in arbitrary units) on the temperature t for different val- ues of magnetic fields h and energy gap ε are presented in figure 3. Naturally, for h < hc, its behaviour is typical of the systems with the first order phase transition point. At h = hc, the curves transform into the characteristic for a system with the second order phase transition point. For h > hc, the magnetic susceptibility behaviour is characteristic of paramagnetic systems. For ε> 0.5 instance, only one curve (h = 0) demonstrates the behaviour typical of the systems with a possible phase transition in the point T = 0 K [figure 3 (d)]. Spontaneousmagnetization here is impossible at any finite temperature, but the system is sensitive to an infinitesimally small external field at T → 0 K. 4. Thermodynamic functions Using the expression (3.5) we can rewrite the equation (2.9) for the free energy in the following form: F =− 1 β ln 2 (1− c) p 1− s2 + 3 2 c2s2Xv J . (4.1) This representation for F takes into account the dependencies of c and s on t , h, ε and will be useful for the calculation of different thermodynamic functions. Having differentiated free energy (4.1) on temperature, we obtained for entropy the following expres- sion: S = N k [ ln 2 (1− c) p 1− s2 + ( 3cs2βXv J − 1 1− c ) X1 + ( 3c2sβXv J − s 1− s2 ) X2 ] , (4.2) where X1 = 2cs(c2s +h)+ (cs2 −ε0)Y2 Y1Y2 −4c2s2 , X2 = (c2s +h)Y1 +2cs(cs2 −ε0) Y1Y2 −4c2s2 , (4.3) 33706-7 M.A. Korynevskii, V.B. Solovyan (a) (b) (c) (d) Figure 4. Temperature dependencies of the entropy at different energy gaps ε and external fields h (num- bers of curves correspond to the same values of h as in figure 2). Y1,Y2 are determined by expressions (3.12). The dependencies of entropy S (in arbitrary units) on the reduced temperature at different values of external field h and energy gap ε are presented in figure 4. All of them demonstrate significant jumps at some temperature points dependent on the field for all h < hc [see figures 4 (a)–(c)]. This confirms the existence of the first order magnetic phase transition. For h = hc there is a noticeable kink in the entropy behaviour in the corresponding temperature point. Such a behaviour is in accordance with the thermodynamic classification of the phase transition order [29, 30]. A specific role of the zero external field at the ε> 0.5 instance is demonstrated in figure 4 (d). At T → 0 K, the total entropy of the system in zero field (h = 0) does not turn into zero, becausemagnetic ordering is not realized. At the same time, all particles remain at the ground level (c = 0) and the configuration part of entropy is equal to zero. Comparison of numerical values of entropy to the right and to the left of the transition point makes it possible to find the latent heat of the first order magnetic phase transition: q = Tc (Sr −Sl ) , (4.4) where Sr = S+ is a limit value of S when temperature falls to the phase transition point Tc, and Sl = S− is a limit value of S when temperature increases to this point. The latent heat (in arbitrary units) behaviour under the external field h is presented in figure 5. The q is not equal to zero for h < hc and vanishes at h = hc for all ε. Every curve possesses a slight maximum connected, in our opinion, with the different temperature behaviours of c and s parameters. Based on the well-known thermodynamic relation: CV =−T ( ∂2F ∂T 2 ) (4.5) 33706-8 Self-consistent approach to magnetic ordering Figure 5. External field dependencies of latent heat at different energy gaps c (curves: 1 — ε = 0, 2 — ε= 0.25, 3— ε= 0.45, 4— ε= 0.495). the heat capacity (specific heat) has been calculated. Being cumbersome the expression forCV is omitted here. The temperature dependencies of heat capacity (in arbitrary units) at different values of external field h at some ε are presented in figure 6. All of them demonstrate an infinite growth at temperatures of the first and the second order phase transitions and a smooth behaviour for h > hc (ε < 0,5) [see figures 6 (a)–(c)]. A specific place of the ε > 0.5 instance [figure 6 (d)] manifests itself in the appearance of an addi- tional maximum of CV for small (curve 2) and large (curve 6) fields. As can be seen by comparison with (a) (b) (c) (d) Figure 6. Temperature dependencies of the heat capacity at different energy gaps ε and external fields h (numbers of curves correspond to the same values of h as in figure 2). 33706-9 M.A. Korynevskii, V.B. Solovyan figures 2 (d) and 2 (h), the right-hand maximum in curve 2 and the left-hand maximum in curve 6 cor- respond to a rapid growth of occupation of the excited level by particles (c), but both supplementary maxima correspond to a rapid growth of magnetization (s). Thus, for ε> 0.5, the temperature behaviours of c and s are noticeably independent. Such a tendency is also visible for ε = 0.45 [see figures 6 (c)]. It should be underlined that the capacity calculated here corresponds only to the site excited and spin ori- entation degrees of freedom. Other important degrees of freedom, vibrational or rotational, for example, are not taken into account here. The points of a visible change of thermodynamic functions behaviour (jumps for order parameters c and s) are connected with a transition of the free energy (4.1) from a certain branch to another one. The numerical analysis shows the existence of 1–3 different branches of F in the whole temperature range from 0 to∞. Two of them realized a minimum value of F and the third one realized its maximum value. Naturally, the real behaviour of the system investigated is determined by the branches with minimum F . The temperature points of critical behaviour are found from the continuity condition of free energy (4.1) in those points [F (Tc,c1, s1) = F (Tc,c2, s2)]: ln (1− c1) √ 1− s2 1 (1− c2) √ 1− s2 2 = 3 2 βcXv J ( c2 2 s2 2 − c2 1 s2 1 ) , (4.6) where c1, s1 and c2, s2 are roots of equations (3.5) at β→ βc for two different branches of F and βc = (kTc)−1 . The analysis of all physical properties of the system investigated shows their strict connection with the behaviour of c and s parameters. Rapid changes of thermodynamic functions take place when the occupation number c jumps from a value above 0.5 to a value below 0.5. According to equations (3.5), this corresponds to a relation: cs2 −ε= 0. (4.7) Consequently, the excitation energy gap ε= 0.5 determines the upper limit of ε, above which all changes in a system have a continuous character. It is hard to compare the theoretical results obtained here with the other ones due to the lack of an exact coincidence of our model with the one well-known in the literature. However, experimental in- vestigations of magnetic phenomena induced by thermal or photo excited electrons in complex metal or metal-organic species are now intensively carried out [31–34]. Photo-induced magnetic structures based on the organic matrix intercalated by Rb and Cs atoms demonstrate thermodynamic properties above the Curie temperature close to the ones presented in figure 3 and figure 6. Those structures in general are stable at low temperatures and only in V-Cr Prussian Blue Analogy compounds instance [34], a critical temperature for magnetization state is about 350 K. Molecular based magnets of this type are a relatively new class of materials which have good prospects to be used as spin-electronics memory cells. The phys- ical mechanisms of their operation can be quite fully explained within the framework of the two-level four state model discussed here. 5. Conclusions A self-consistent theory for a description of the excited energy levels occupation and spin-orientation ordering in the crystalline lattice is proposed. The Ising-like magnetic interaction between particles in the excited states are taken into account. When this interaction is omitted, the occupation number c and the spin-order parameter s become continuous functions of temperature and an external field, which is characteristic of two independent subsystems: noninteracting two-level site system and diluted magnetic system. The exchange interaction being taken into account significantly changes the situation. Depending on the excited energy gap and on the magnitude of the external magnetic field, the transition to a sponta- neous magnetically ordered state of the system is possible only for the fields less than the critical one. Moreover, at the field less than critical, there holds the first order phase transition, and when this field turns into a critical value, the transition becomes the second order. For the fields bigger than the critical 33706-10 Self-consistent approach to magnetic ordering one, the spontaneous magnetization of a system is suppressed by an external field (an induced magneti- zation is proportional to the field magnitude). Those effects take place for the reduced excited energy gap ε≡ ε0(Xv J )−1 less than 0.5. For εÊ 0.5, no spontaneous magnetic ordering is possible, total magnetization is completely determined by the external magnetic field. The generalized (dependent on the external field) magnetic susceptibility χ behaviour is in accor- dance with the microscopic parameters c and s temperature dependencies. Its infinite peak at the phase transition points, for fields less than and equal to the critical fields, moves to a significant peak at higher fields. The width of χ curve near Tc is much smaller for the first order phase transition instance (h < hc) compared to the second order phase transition instance (h = hc). The latent heat q of the first order phase transition in the system investigated is strongly field depen- dent, vanishing for h Ê hc in accordance with the change of the order of the phase transition. The heat capacity (specific heat) behaves in a usual manner and for h Ê hc remains finite for all temperatures. All those peculiarities for χ, q , CV depend on the ε value in accordance with the above mentioned objections. A specific place of the excited level gap ε = 0.5 is established as a limit value between the jump-like and continuous behaviour of thermodynamic functions of the investigated system. Based on all of the research, we may conclude that in the investigated system, correlation between the excited level occupation and magnetic ordering plays a decisive role. Specific percolation effects con- cerning the magnetic properties of a system take place. They manifest themselves in a self-consistent way. Disappearance of magnetic ordering entails a sharp decrease in the probability of occupation of the excited levels by particles. The additional stimuli (optical radiation for example) being taken into account for particle excitation are quite promising in the study of the mechanisms of photo-induced magnetic phenomena, which are candidates for future applications in spin-electronics. The proposedmodel also possesses other peculiarities based on the complex ratio between the excited level gap (ε) and the magnitude of the external field (h). A detailed study of these peculiarities will be a matter of future investigations. Acknowledgements Authors thank Prof. I.V. Stasyuk for his interest to the present paper and for useful discussions of the subject of this research. References 1. Feynman R.P., Leighton R.B., Sands M.L., The Feynman Lectures on Physics, Addison-Wesley, Reading, 2005. 2. Griffiths D., Introduction to Quantum Mechanics, Pearsen-Prentice Hall, New Jersey, 2004. 3. Capitelli M., Colonna G., D’Angola A., Fundamental Aspects of Plasma Chemical Physics, Springer, New York, 2012. 4. Messina A., Nakazato H., J. Phys. A: Math. Theor., 2014, 47, 445302; doi:10.1088/1751-8113/47/44/445302. 5. Tritschler T., Mücke O.D., Wegener M., Phys. Rev. A, 2003, 68, 033404; doi:10.1103/PhysRevA.68.033404. 6. Bonifacio R., Lugiato L.A., Phys. Rev. A, 1975, 11, 1507; doi:10.1103/PhysRevA.11.1507. 7. Zagoskin A.M., Ashhab S., Johansson J.R., Nori F., Phys. Rev. Lett., 2006, 97, 077001; doi:10.1103/PhysRevLett.97.077001. 8. Ralph D.C., Buhrman R.A., Phys. Rev. Lett., 1992, 69, 2118; doi:10.1103/PhysRevLett.69.2118. 9. Arnold W., Hunklinger S., Solid State Commun., 1975, 17, 883; doi:10.1016/0038-1098(75)90743-7. 10. Foresman J.B., Head-Gordon M.S., Nakano H., J. Phys. Chem., 1992, 96, 135; doi:10.1021/j100180a030. 11. Dreuw A., Head-Gordon M.S., Chem. Rev., 2005, 105, 4009; doi:10.1021/cr0505627. 12. Yukhnovskii I., Principles of Quantum Mechanics, Lybid’, Kyiv, 1995 (in Ukrainian). 13. Schubin S., Wonsowsky S., Proc. R. Soc. Lond. A, 1934, 145, 159; doi:10.1098/rspa.1934.0089. 14. Vonsovskii S.V., Magnetism, Wiley, New York, 1974. 15. Tyablikov S.V., Methods in the Quantum Theory of Magnetism, Springer, New York, 1995. 16. Blume M., Emery V.J., Griffits R.B., Phys. Rev. A, 1971, 4, 1071; doi:10.1103/PhysRevA.4.1071. 17. Capel H.W., Physica, 1966, 32, 966; doi:10.1016/0031-8914(66)90027-9. 33706-11 http://dx.doi.org/10.1088/1751-8113/47/44/445302 http://dx.doi.org/10.1103/PhysRevA.68.033404 http://dx.doi.org/10.1103/PhysRevA.11.1507 http://dx.doi.org/10.1103/PhysRevLett.97.077001 http://dx.doi.org/10.1103/PhysRevLett.69.2118 http://dx.doi.org/10.1016/0038-1098(75)90743-7 http://dx.doi.org/10.1021/j100180a030 http://dx.doi.org/10.1021/cr0505627 http://dx.doi.org/10.1098/rspa.1934.0089 http://dx.doi.org/10.1103/PhysRevA.4.1071 http://dx.doi.org/10.1016/0031-8914(66)90027-9 M.A. Korynevskii, V.B. Solovyan 18. Chen H.H., Levy P.M., Phys. Rev. Lett., 1971, 27, 1383; doi:10.1103/PhysRevLett.27.1383. 19. Lajzerowicz J., Sivardiere J., Phys. Rev. A, 1975, 11, 2079; doi:10.1103/PhysRevA.11.2079. 20. Sivardiere J., J. Phys. (France), 1990, 41, 1081; doi:10.1051/jphys:0198000410100108100. 21. Ausloos M., Clippe P., Kowalski J.M., Pekalska J., Pekalski A., Phys. Rev. A, 1983, 28, 3080; doi:10.1103/PhysRevA.28.3080. 22. Sokolovskii R.O., Phys. Rev. B, 2000, 61, 36; doi:10.1103/PhysRevB.61.36. 23. Stasyuk I.V., Tovstyuk K.D., Gera O.B., Velychko O.V., Preprint of the Institute for Condensed Matter Physics, ICMP-02-09U, Lviv, 2002 (in Ukrainian). 24. Ma S.K., Modern Theory of Critical Phenomena, Benjamin, Reading, 1976. 25. Stanley H.E., Introduction to Phase Transitions and Critical Phenomena, Oxford University Press, Oxford, 1987. 26. Landau L.D., Zh. Eksp. Teor. Fiz., 1937, 7, 19. 27. Landau L.D., Lifshytz E.M., Statistical Physics, Elsevier, Oxford, 2005. 28. Stauffer D., Ahariny A. Introduction to Percolation Theory, Taylor and Francis, London, 1994. 29. Ehrenfest P., Proc. Kon. Ned. Akad. Wetensch. te Amsterdam, 1933, 36, 153. 30. Reichl L.E., A Modern Course in Statistical Physics, 2nd Edn., Wiley, New York, 1998. 31. Pashechko V.A., Borisenko V.G., Khats’ko E.N., Cherny A.S., Paulsen C., Feher A., Orendac M., Orendacova A., Jansen A.G.M., Wyder P., J. Phys.: Condens. Matter, 2002, 14, 9693; doi:10.1088/0953-8984/14/41/324. 32. Varres F., Goujon A., Boukheddaden K., Nogues M., Bleuzen A., Verdaguer M., Mol. Cryst. Liq. Cryst., 2002, 379, 333; doi:10.1080/713738677. 33. Raghunathan R., Ramasesha S., Mathoniere C., Marvaud V., Phys. Rev. B, 2006, 73, 045131; doi:10.1103/PhysRevB.73.045131. 34. Bozdag K.D., Yoo J.-W., Raju N.P., McConnell A.C., Miller J.S., Epstein A.J., Phys. Rev. B, 2010, 82, 094449; doi:10.1103/PhysRevB.82.094449. Самоузгоджений пiдхiд до процесiв магнiтного впорядкування та заповнення збуджених станiв у дворiвневiй системi М.А. Кориневський1,2,3, В.Б. Солов’ян1 1 Iнститут фiзики конденсованих систем НАН України, вул. I. Свєнцiцького, 1, 79011 Львiв, Україна 2 Нацiональний унiверситет “Львiвська полiтехнiка”, вул. С. Бандери, 12, 79013 Львiв, Україна 3 Iнститут фiзики Щецiнського унiверситету, вул. Вєлькопольська, 15, 70451Щецiн, Польща Детально вивчено процес феромагнiтного впорядкування у дворiвневiй частково збудженiй системi. Величина намагнiченостi (магнiтний параметр впорядкування) та граткове впорядкування (заповнення збуджених рiвнiв) розрахованi самоузгоджено. Обговорено вплив зовнiшнього магнiтного поля i шири- ни енергетичної щiлини мiж основним i збудженим рiвнями на феромагнiтний фазовий перехiд. Ключовi слова: магнетики, дворiвнева система, фазовий перехiд 33706-12 http://dx.doi.org/10.1103/PhysRevLett.27.1383 http://dx.doi.org/10.1103/PhysRevA.11.2079 http://dx.doi.org/10.1051/jphys:0198000410100108100 http://dx.doi.org/10.1103/PhysRevA.28.3080 http://dx.doi.org/10.1103/PhysRevB.61.36 http://dx.doi.org/10.1088/0953-8984/14/41/324 http://dx.doi.org/10.1080/713738677 http://dx.doi.org/10.1103/PhysRevB.73.045131 http://dx.doi.org/10.1103/PhysRevB.82.094449 Introduction Hamiltonian and free energy Order parameters c and s Thermodynamic functions Conclusions