CC BY
12
УДК 537.6
Features of the Non-equilibrium Critical Dynamics in 3D Pure and Diluted Ising-like Ferromagnets
Vyacheslav V. Krizhanovskiy* Peter N. Malyarenko^ Pavel V. Prudnikov* Vladimir V. Prudnikov§
Department of Theoretical Physics Dostoevsky Omsk State University Mira, 55a, Omsk, 644077 Russia
Received 20.08.2016, received in revised form 10.10.2016, accepted 12.11.2016 The non-equilibrium critical behaviour of the three-dimensional pure and site-diluted spin systems described by Ising model is studied with different spin concentrations for evolution from various initial magnetic states. The universal scaling functions are determined for the magnetization and autocorrelation function. The aging and superaging effects are revealed in autocorrelation function behaviour for diluted systems with evolution from low-temperature initial states.
Keywords: non-equilibrium critical behaviour, Ising model, initial state, structure defects, correlation function.
DOI: 10.17516/1997-1397-2016-9-4-463-468.
Introduction
The investigation of critical behaviour of disordered systems remains one of the main problems in condensed-matter physics and excites a great interest because all real solids contain structural defects. The structural disorder breaks the translational symmetry of the crystal and thus greatly complicates the theoretical description of the material. Then, to describe the influence of disorder on the critical behaviour the methods of computer calculations are used from which Monte-Carlo methods are main. The development of computer technology and increase of their capacity contributes to this direction of investigations.
Critical dynamics provides a simple instance of slow collective evolution. Statistical systems near the critical temperature are displayed by dramatic slowing down of relaxation processes, hysteresis, memory effects and etc. The influence of disorder is particularly important near critical point Tc where the behaviour of a system is characterized by an anomalous large relaxation time described by the law trel ~ \T— Tc\-zv, where 2 is the dynamic critical exponent and v is the correlation length exponent. At a time t C trel, the effects of aging manifest themselves, namely, at the two-time correlation and response functions on the waiting time tw and observation time t — tw. An effect of initial states influence on the non-equilibrium critical behaviour is displayed within the t,tw C trel time range [1, 2].
* FehtunKrizh@gmail.com
t petr.malyarenko@yandex.ru
^ prudnikp@univer.omsk.su
§ prudnikv@univer.omsk.su © Siberian Federal University. All rights reserved
1. Characteristics of non-equilibrium critical behaviour and results of simulations
In investigations of influence of the initial states on characteristics of non-equilibrium critical behaviour, there are differed the high-temperature initial state, when a system has been prepared at a temperature To ^ Tc with the initial magnetization m0 = 0 before it was quenched at Tc, and the low-temperature magnetized initial state with preparation at T0 < Tc with m0 = 0. The non-equilibrium evolution starts when the system after initial preparation is placed to the thermostat with T = Tc. At the waiting time tw, we have begun to measure two-time quantities such as response and correlation functions. We have considered the ferromagnetic Ising model site-diluted by nonmagnetic impurity atoms. The Hamiltonian of system is given by
H = —J PiPjSiSj, (1)
<i,j>
where J > 0 is short-range exchange interaction integral between spins Si = ±1 fixed at the lattice sites . The random occupation numbers pi take the values 0 or 1, and pi equals 1 if the site of the lattice contains spin and 0 otherwise. We considered the cubic lattice with periodic boundary conditions. Let us denote Ns = pL3 as the number of spins in a lattice with linear size L = 64 and 128, where p is the spin concentration.
To investigate the non-equilibrium evolution, we computed such characteristics as the magnetization M(t) and the two-time dependent autocorrelation function C(t,tw), which are defined by relations
M (t) = ddx [{S(x,t))] =
C (t,tw) = ddx[{S(x,t)S(x,tw)) — (S(x,t)) (S (x,tw )>],
where angle brackets stands for the mean over the different realizations of initial state, the square brackets are for averaging over the different impurity configurations.
At the present time, it becomes established that the time dependence for magnetization and autocorrelation function for systems starting from a low-temperature initial state with m0 = 0 satisfies the following scaling forms [1]:
M (t,tm) = Amt-ii/zvFm(t/tm), (3)
C (t,tw ,tm) = Ac (t — tw )a+1-d/z (t/tw )9-1 Fc (tw/t,t/tm),
where tm = Bmm0 1/k is a new timescale with exponent k = 9' + ¡3/(vz) > 0, a = (2 — n — z)/z, d = 9' — a, 9' is the initial slip exponent. The scaling functions FM(t/tm) and FC(tw/t,t/tm) are finite at limits tw ^ 0 and t/tm ^ 0, AM and AC are non-universal amplitudes.
At the beginning, we realized investigation of critical relaxation of magnetization M(t) for pure (p = 1) and diluted systems (p = 0.8, 0.6, and 0.5) with various initial states m0 (Fig. 1a), which demonstrate that the relaxation time is increased with growth of the impurity concentration cimp = 1 — p. The curves in Fig. 1b and c for pure and diluted systems with p = 0.6 demonstrate an essential differences in relaxation from the high-temperature with m0 ^ 1 and low-temperature with m0 = 1 initial states, and from the intermediate initial states with 0.2 ^ m0 ^ 0.6. So, a remarkable property of the non-equilibrium critical relaxation from the high-temperature initial state is the increase of magnetization M(t) ~ te with 9' > 0. The initial rise of magnetization is changed to the well known decay M(t) ~ t-ii/zv for t > tcr ~ m0 1/k. The critical relaxation from the completely ordered initial state with m0 = 1 is immediately
Ns
—J2 piSi(t)
(2)
characterized by power-law of magnetization decay M(t) ~ t-l3/zu. The intermediate cases with 0.2 < m0 < 0.6 and 0.4 are characterized by a short stage of magnetization rise under the power-law M(t) ~ te which is changed to much long-continued stage of relaxation with M(t) ~ t-l3/zu.
s
o.i
1 10 100 1000 10000 t, MCS/s
Fig. 1. Time dependence of the magnetization M(t) for various values of the initial magnetization m0 and different spin concentrations p (a), for pure model, p =1.0 (b) and diluted system, p = 0.6 (c)
To check the scaling prediction for M(t) as function of the initial magnetization m0 given by relation (3), we plot the dependence of t^/zvM(t) versus x = tm0/k in Fig. 2 with the use of exponent values for the pure 3D Ising model: ¡/v = 0.516(2), 2 = 2.024(6), and diluted systems: 3/v = 0.508(16), 2 = 2.191(21) for p = 0.8 and ¡3/v = 0.462(40), 2 = 2.663(30) for p = 0.6 [3]. We can see in Fig. 2a, b and c the collapse of curves FM (x) for different m0 with fixed spin concentration p into a single universal curve with linear initial curve piece (in log - log scale) with FM(x) ~ xk. We obtained the following values of exponent k: for pure case k = 0.361(3), for systems with p = 0.8 k = 0.353(4) and with p = 0.6 k = 0.325(7).
The dependence of the autocorrelation function on the initial magnetization m0 in concordance with relation (3) can be determined when the waiting time is taken as a quantity proportional to the observation time, as example, tw = t/3. Then, the scaling dependence for the autocorrelation function in (3) takes on form:
C(t,tw = t/3, tm) ~ t-2^/zvGc(tm0/k). (4)
We presented in Fig. 3a and b the time dependence of function C(t,tw = t/3,tm) for pure and diluted with p = 0.8 models for various values of the initial magnetization m0, and in Fig. 3c
Fig. 2. Dependence of the magnetization scaling functions FM(t/tm) = M(t)t^/ziJ on x = tm])/k for different initial magnetizations m0 and spin concentrations p = 1.0 (a), p = 0.8 (b) and p = 0.6 (c)
and d the dependence of the scaling functions Gc(tm0/k) = t2i3/zvC(t,tw = t/3,tm) on variable x = tm])/k with k = 0.361 for pure model and k = 0.353 for system with p = 0.8. The curves in Fig. 3a and b demonstrate the aging effects through deflection of the curves for time dependence of C(t,t/3) from power law, which must look as straight lines at log-log scale, and characterize the slowing down of time correlations with increasing system age tw ~ t. One can see that the increase of m0 leads to the amplification of aging effects.
The graph in Fig. 3c demonstrates the scaling collapse of function Gc (x) for various values of m0 into an universal curve that confirms the scaling forms (3) and (4) for pure Ising model. However, the graph in Fig. 3d shows the violation of scaling form (4) for lack of the data collapse for dilute system. We connect this phenomenon with essential influence of defects on non-equilibrium dynamics of magnetic domain structure and therefore on non-equilibrium correlation properties of diluted system, which lead to superaging effects [2]. When we use the presentation of autocorrelation function in the form t2l3/^zvC(t,t1/^,tm) on x = t^m^/k it gives possibility to get with the use of superaging exponent n = 2.30(6) (p = 0.8) [2] the data collapse into a single curve for values of the initial magnetization m0 > 0.25 (Fig. 3f). Thus, the autocorrelation function is characterized in the aging regime with t > tW by more complicated scaling dependence in the form C(t,tw,tm) ~ t—2l3/zvFc(t/tW,t/tU^) for the site-diluted Ising model.
Investigations were supported by Russian Science Foundation, project 14-12-00562.
ioJ 10J 10"' 10° io' io5 io! 104 10" 10s io' io"
Fig. 3. Dependence of the autocorrelation function C(t,tw = t/3,tm) on time t with different initial magnetizations m0 for systems with p =1.0 (a) and p = 0.8 (b) and its scaling functions Gc(t/tm) = C(t,tw = t/3,tm)t2l3/zv on x = tm0/k for systems with p = 1.0 (c) and p = 0.8 (d). Superaging effects for scaling function t2^/^zvC(t,t1/tm) on x = t^m0/k for system with p = 0.8 (f)
References
[1] P.Calabrese, A.Gambassi, F.Krzakala, J. Stat. Mech., (2006), P06016.
[2] V.V.Prudnikov, P.V.Prudnikov, E.A.Pospelov, J. Stat. Mech., (2016), 043303.
[3] P.V.Prudnikov, V.V.Prudnikov, E.A.Pospelov, P.N. Malyarenko, A.N.Vakilov, Prog. Theor. Exp. Phys., (2015), 053A01.
Особенности неравновесной критической динамики в чистых и разбавленных трехмерных изингоподобных ферромагнетиках
Вячеслав В. Крижановский Петр Н. Маляренко Павел В. Прудников Владимир В. Прудников
Физический факультет Омского государственного университета Мира, ббя, Омск, 644077 Россия
Представлены результаты численного исследования особенностей неравновесного критического поведения в трехмерной чистой и структурно неупорядоченной модели Изинга с широким спектром изменения спиновой концентрации при её эволюции из различных начальных магнитных состояний. Выявлены эффекты старения и сверхстарения в двухвременном поведении автокорреляционной функции для структурно неупорядоченных систем при эволюции из низкотемпературных начальных состояний.
Ключевые слова: неравновесное критическое поведение, модель Изинга, начальное состояние, дефекты структуры, корреляционная функция.
