arxiv: v1 [cond-mat.stat-mech] 9 Jan 2012

Transcription

1 arxiv: v1 [con-mat.stat-mech] 9 Jan 2012 Externally riven one-imensional Ising moel Amir Aghamohammai a 1, Cina Aghamohammai b 2, & Mohamma Khorrami a 3 a Department of Physics, Alzahra University, Tehran , Iran b Department of Electrical Engineering, Sharif University of Technology, , Tehran, Iran Abstract A one imensional kinetic Ising moel at a finite temperature on a semiinfinite lattice with time varying bounary spins is consiere. Exact expressions for the expectation values of the spin at each site are obtaine, in terms of the time epenent bounary conition an the initial conitions. The solution consists of a transient part which is ue to the initial conition, an a part riven by the bounary. The latter is an evanescent wave when the bounary spin is oscillating harmonically. Low- an high-frequency limits are investigate with greater etail. The total magnetization of the lattice is also obtaine. It is seen that for any arbitrary rapily varying bounary conitions, this total magnetization is equal to the the bounary spin itself, plus essentially the time integral of the bounary spin. A nonuniform moel is also investigate. 1 mohamai@alzahra.ac.ir 2 c aghamohammai@yahoo.com 3 mamwa@mailaps.org 1

2 1 Introuction Dynamical spin systems have playe a central role in non-equilibrium statistical moels. The Ising moel is wiely stuie in statistical mechanics, as it is simple an allows one to unerstan many features of phase transitions. The non-equilibrium properties of the Ising moel follow from the spin ynamics. In his article [1], Glauber introuce a ynamical moel formulating the ynamics of spins, base on the rates coming from a etaile balance analysis. It is a simple non-equilibrium moel of interacting spins with spin flip ynamics. An extension of the kinetic Ising moel with nonuniform coupling constants on a one-imensional lattice was introuce in[2]. In[3], a amage spreaing metho was use to stuy the sensitivity of the time evolution of a kinetic Ising moel with Glauber ynamics against the initial conitions. The full time epenence of the space-epenent magnetization an of the equal time spin-spin correlation functions were stuie in[4]. Non-equilibrium two-time correlation an response functions for the ferromagnetic Ising chain with Glauber ynamics have been stuie in [5, 6]. The ynamics of a left-right asymmetric Ising chain has been stuie in [7]. The response function to an infinitesimal magnetic fiel for the Ising-Glauber moel with arbitrary exchange couplings was aresse in [8]. In [9], a Glauber moel on a one-imensional lattice with bounaries was stuie, for both ferromagnetic an anti-ferromagnetic couplings. The largetime behavior of the one-point function was stuie. It was shown that the system exhibits a ynamical phase transition, which is controlle by the rate of spin flip at the bounaries. It was shown in [10, 11] that for a nonuniform extension of the kinetic Ising moel, there are cases where the system exhibits static an ynamical phase transitions. Using a transfer matrix metho, it was shown that there are cases where the system exhibits a static phase transition, which is a change of behavior of the static profile of the expectation values of the spins near en points [10]. Using the same metho, it was shown in [11] that a ynamic phase transition coul occur as well: there is a fast phase where the relaxation time is inepenent of the reaction rates at the bounaries, an a slow phase where the relaxation time oes epen on the reaction rates at the bounaries. Most of the stuies on reaction iffusion moels have been on cases where the bounary conitions are constant in time. Among the few moels with time epenent bounary conitions, is the asymmetric simple exclusion process on a semi-infinite chain couple at the en to a reservoir with a particle ensity that changes perioically in time [12]. The situation is similar regaring the case of the kinetic Ising moel as well. Among the exceptions are the stuy of the ynamical response of a two-imensional Ising moel subject to a square-wave external fiel [13], an the stuy of a harmonic oscillator linearly couple with a linear chain of Ising spins [14, 15]. In this article a one imensional Ising moel at temperature T on a semiinfinite lattice with time varying bounary spin is investigate. The paper is organize as follows. In section 2 a brief review of the formalism is presente, mainly to introuce the notation. In section 3, a semi-infinite lattice with 2

3 oscillating bounary spin is stuie. The exact solution for the expectation values of the spin at any site is obtaine. It is shown that there the bounary prouces an evanescent wave in the lattice. The low an high frequency limits are stuie in greater etail. The total magnetization of the lattice, M(t), is also obtaine. It is shown that for rapily changing bounary conitions, the total magnetization is equal to the the bounary spin itself, plus something proportional to the time integral of the bounary spin. A nonuniform moel in also investigate. It is shown that its evolution operator eigenvalues are real. For the specific case of a two-part lattice with each part being homogeneous, the reflection an transmission coefficients corresponing to a harmonic source at the en of the lattice are calculate. Finally, section 4 is evote to the concluing remarks. 2 One-imensional Ising moel with nonuniform coupling constants Consier an Ising moel on a one-imensional lattice with L sites, labele from 1 to L. At each site of the lattice there is a spin interacting with its nearest neighboring sites accoring to the Ising Hamiltonian. At the bounaries there arefixe magneticfiels. Denoting the spin at the site j by s j, an the magnetic fiel at the sites 1 an L by B 1 an B L, one has for the Ising Hamiltonian H = L µ α=1+µ where J α is the coupling constant in the link α, an J α s α µ s α+µ B 1 s 1 B L s L. (1) µ = 1 2. (2) The link α links the sites α µ an α + µ, so that α ± µ are integers, an α runs from µ up to (L µ). Throughout this paper, sites are enote by Latin letters which represent integers, while links are enote by Greek letters which represent integers plus one half (µ). The spin variable s j takes the values +1 for spin up ( ), or 1 for spin own ( ). Define βj α, 1 < α < L K α := βb 1, α = µ (3) βb L, α = L+µ where β := 1 k B T, (4) an k B is the Boltzmann s constant, an T is the temperature. Denoting the reaction rate from the configuration A to the configuration B by ω(a B), 3

4 an assuming that in each step only one spin flips, etaile balance emans the following for the reaction rates. ω[(s,s j ) (S, s j )] = Γ j [1 s j tanh(k j µ s j 1 +K j+µ s j+1 )], 1 < j < L, (5) ω[(s,s 1 ) (S, s 1 )] = Γ 1 [1 s 1 tanh(k µ +K 1+µ s 2 )], (6) ω[(s,s L ) (S, s L )] = Γ L [1 s L tanh(k L µ s L 1 +K L+µ )]. (7) Γ j s are inepenent of the configurations. For simplicity, we take them to be inepenent of the site. Then, rescaling the time they are set equal to one. So the evolution equation for the expectation value of the spin in the site j is t s j = 2 s j +[tanh(k j µ +K j+µ )+tanh(k j µ K j+µ )] s j 1 +[tanh(k j µ +K j+µ ) tanh(k j µ K j+µ )] s j+1, 1 < j < L, t s 1 = 2 s 1 +[tanh(k µ +K 1+µ )+tanh(k µ K 1+µ )] +[tanh(k µ +K 1+µ ) tanh(k µ K 1+µ )] s 2, t s L = 2 s L +[tanh(k L µ +K L+µ )+tanh(k L µ K L+µ )] s L 1 +[tanh(k L µ +K L+µ ) tanh(k L µ K L+µ )]. (8) These can be written in the form t s j = 2 s j +[tanh(k j µ +K j+µ )+tanh(k j µ K j+µ )] s j 1 +[tanh(k j µ +K j+µ ) tanh(k j µ K j+µ )] s j+1, 1 j L, (9) s 0 = 1, (10) s L+1 = 1. (11) 3 Time varying bounary conitions on a semiinfinite lattice Consier a lattice for which the bounary spins (s 0 an s L+1 ) are externally controlle, but the reactions at the internal sites satisfy etaile balance. The evolution equation is then the same as (9), but combine with bounary conitions ifferent from (10) an (11). A semi-infinite lattice the bounary of which is externally controlle, is obtaine by letting L ten to infinity, an using the following bounary conitions s 0 = f(t), (12) s j oes not blow up as j tens to infinity, (13) 4

5 instea of (10) an (11). A general solution of (9), combine with (12) an (13), can be written as the sum of a particular solution plus a general solution of (9), combine with the homogeneous bounary conitions. 3.1 Semi-infinite lattice with uniform couplings: the homogeneous solution For a lattice with uniform couplings, K α s are enote by K. The solution to the homogenous equation (vanishing f) is enote by s j h. One arrives at Defining one arrives at t s j h = 2 s j h +[tanh(2k)]( s j 1 h + s j+1 h ), 0 < j. (14) s j h := s j h, j < 0, (15) t s j h = 2 s j h +[tanh(2k)]( s j 1 h + s j+1 h ), (16) which hols for all integers j. Denoting the linear operator acting on s l s in the right-han sie of (16) by h, the above equation is of the form t s j h = h l j s l h, (17) where h l j s are the matrix elements of h. Defining the generating function G through G(z,t) := z j s j h (t), (18) one arrives at resulting in j= G t = [ 2+(z +z 1 ) tanh(2k)]g, (19) G(z,t) = exp{[ 2+(z +z 1 ) tanh(2k)]t}g(z,0), = exp( 2t) z k I k [2t tanh(2k)]g(z,0), = exp( 2t) k= j= z j l= I j l [2t tanh(2k)] s l h (0), (20) 5

6 where I k is the moifie Bessel function of first kin of orer k. (20) results in s j h (t) = exp( 2t) = exp( 2t) l= I j l [2t tanh(2k)] s l h (0), {I j l [2t tanh(2k)] I j+l [2t tanh(2k)]} s l h (0). l=1 (21) Using the large argument behavior of the moifie Bessel functions, it is seen that s j h (t) exp{ 2[1 tanh(2k)]t}, j > 0, (22) showing that the homogeneous solution tens to zero at large times. 3.2 Semi-infinite lattice with uniform couplings: the particular solution corresponing to harmonic bounary conitions The harmonic bounary conition is s 0 = Re[σ 0 exp( iωt)] (23) The following ansatz for a particular solution s j p to equations (9) an (12) results in This has a solution of the form where z satisfies s j p = Re[σ j exp( iωt)] (24) (iω 2)σ j +[tanh(2k)](σ j+1 +σ j 1 ) = 0. (25) σ j = cz j, (26) z +z 1 = iω+2 tanh(2k). (27) It is obvious that changing the sign of K results in changing the sign of z, while changing the sign of ω results in changing z to its complex conjugate. So it is sufficient to consier only nonegative values of K an ω. From now on, it is assume that K an ω are nonnegative. (27) has two solution for z, which are inverse of each other, an none are unimoular. The bounary conition at infinity imposes that of the two solutions of type (26), only that solution is acceptable which correspons to the root of (27) with moulus less than one. From now, only this root is enote by z: z := r exp(iθ) (28) 6

7 where r an θ are real an r is positive an less that one. The solution to (25) is then σ j = σ 0 z j. (29) As z is less than one, the particular solution (24) escribes an evanescent wave. Obviously, the rate of ecay length an the phase spee, l an v respectively, satisfy l = 1 lnr, v = θ ω. (30) As the homogenous solution (21) tens to zero for large times, the particular solution (24) is in fact the large times solution to the problem of harmonic bounary conition. Defining a := ω 2, b := tanh2k, the real an imaginary parts of (27) rea u := r +r 1, (31) 2 u cosθ = 1 b, u2 1 sinθ = a b. (32) So u satisfies b 2 u 4 (a 2 +b 2 +1)u 2 +1 = 0, (33) from which one arrives, for the solution which is larger than one, at ( 1+a 2 +b 2 + ) 1/2 (1+a u = 2 +b 2 ) 2 4b 2 2b 2. (34) This is increasing with respect to a, an ecreasing with respect to b. Noting that u r = 1 (1 1r ) 2 2, (35) which shows that u is ecreasing with respect to r, it is seen that r is ecreasing with respect to ω, an increasing with respect to K. One also has r = u u 2 1. (36) 7

8 Regaring θ, ifferentiating the first equation in (32)with respect to a, one has cosθ u θ u sinθ = 0, (37) a a resulting in θ a = cosθ u u sinθ a, b 2 u 2 b 2 = bu sinθ (1+a 2 +b 2 ) 2 4b 2 a, = Equation (34) shows that b 2 u 2 b 2 sinθ bu(b 2 u 2 u 2 ) a. (38) bu 1, (39) from which it is seen that 0 < θ a sinθ a. (40) The first inequality shows θ is an increasing function of a, so it is an increasing function of ω. The secon inequality results in θ a θ a, (41) which shows that (θ/a) is a ecreasing function of a. So (a/θ) is an increasing function of a, or (ω/θ) is an increasing function of ω. One also has (2b 2 u 2 ) a 2 +b 2 1 b 2 = 1+ (42) (1+a2 +b 2 ) 2 4b2, an as (1+a2 +b 2 ) 2 4b 2 1 b 2, (43) it turns out that (bu) is increasing with b, so that θ is increasing with b. Hence (ω/θ) is ecreasing with K. The asymptotic behavior of r an θ is summarize as tanhk, ω 1, r = tanh(2 K), 1 ω, ω tanh(2 K) 4+ω 2, K 1,, (44) 8+ω2 +ω 16+ω 2 ω 2 +ω 16+ω 2, 1 K 8 8

9 r tanh(2 K) ω Figure 1: The plot of r versus ω for ifferent values of tanh(2k) an ω cosh(2k), ω 1, 2 π 2, 1 ω, θ = tan 1 ω, (45) 2, K 1, cos 1 8+ω 2 ω 16+ω 2, 1 K 8 Among other things, it is seen that the phase spee, at low frequencies approaches the constant value 2/[cosh(2 K)], while at high frequencies varies like (2ω/π). Figure 1 is a plot of r versus ω for ifferent values of tanh(2k) from 0.1 to 0.9. Figure 2 is a plot of the phase spee (ω/θ) versus ω for ifferent values of tanh(2k) from 0.1 to 0.9. The total magnetization, efine as the sum of the expectation values of the spins, is enote by M. At large times only the particular solution contributes to the magnetization. So, [ ] σ0 M = Re 1 z exp( iωt), t. (46) For the time-inepenent bounary conition, this leas to σ 0 M =, (t, ω = 0). (47) 1 tanhk 9

10 (ω/θ) 0.1 tanh(2 K) ω Figure 2: The plot of the phase spee (ω/θ) versus ω for ifferent values of tanh(2k) For high frequencies, [ M = Re {σ 0 1 tanh(2k) ] 1 exp( iωt)}, (t, ω ). (48) iω This can be simplifie to [ M = Re {σ 0 1+ tanh(2k) ] iω where } exp( iωt), = s 0 (t)+[tanh(2k)][s 0 (t) S 0 ], (t, ω ), (49) t S 0 (t) := t s 0 (t ), 0 [ ] 1 S 0 := lim t S 0 (t). (50) T T One then arrives at a similar result for the magnetization when the bounary conition is any arbitrary rapily varying function of time (so that its low frequency components are negligible): M = s 0 (t)+[tanh(2k)][s 0 (t) S 0 ], T (t, rapily varying bounary conitions). (51) 10

11 3.3 Semi-infinite lattice with two parts of uniform couplings: the particular solution corresponing to harmonic bounary conitions Consier a semi-infinite lattice consisting of two parts, so that { K 1, α < N K α = K 2, α > N, (52) The time evolution equations for the expectation values of the spins are where ṡ j = 2 s j +[tanh(2k 1 )]( s j 1 + s j+1 ), 0 < j < N, (53) ṡ N = 2 s N +κ s N 1 +κ + s N+1, (54) ṡ j = 2 s j +[tanh(2k 2 )]( s j 1 + s j+1 ), N < j. (55) κ := tanh(k 1 +K 2 )+tanh(k 1 K 2 ), κ + := tanh(k 1 +K 2 ) tanh(k 1 K 2 ). (56) Applying a harmonic bounary conition (23), one has for the particular solution of the kin (24), { (A 1 z j 1 σ j = +B 1z j 1 ), 0 j N A 2 z j 2, N j, (57) where z l +z 1 l = iω+2, l = 1,2 (58) tanh(2k l ) an z l is smaller than one. The bounary conition results in From (57) for j = N, one arrives at Finally, (54) results in A 1 +B 1 = σ 0. (59) A 1 z N 1 +B 1z N 1 = A 2 z N 2. (60) κ (A 1 z N 1 1 +B 1 z N+1 1 )+κ + A 2 z N+1 2 = ( iω +2)A 2 z N 2. (61) Equations (59) through (61) give A 1 = B 1 = A 2 = (κ z 1 +κ + z 2 +iω 2)z1 N σ 0 κ (z1 N+1 z1 N 1 )+(κ + z 2 +iω 2)(z1 N z1 N), (κ z1 1 +κ + z 2 +iω 2)z1 N σ 0 κ (z1 N+1 z1 N 1 )+(κ + z 2 +iω 2)(z1 N z1 N), κ (z 1 z1 1 )z N 2 σ 0 κ (z1 N+1 z1 N 1 )+(κ + z 2 +iω 2)(z1 N (62) z1 N ). 11

12 For large N, these simplify to A 1 = σ 0, 1 σ 0 B 1 = (κ z1 1 +κ + z 2 +iω 2)z 2N κ z 1 +κ + z 2 +iω 2 A 2 = κ (z 1 z 1 1 )z N 2 z N 1 σ 0 κ z 1 +κ + z 2 +iω 2. (63) It can be easily shown that for the nonuniform lattice an at high frequencies, up to first term in ω 1 the magnetization is similar to the case of the uniform lattice. 3.4 Semi-infinite lattice with nonuniform couplings: the relaxation times The general solution of (9) with (12) an (13) is the sum of a particular solution an the general solution to (9) an (12) an (13) with vanishing f. The latter (the homogeneous solution) satisfies t s j h = 2 s j h +[tanh(k j µ +K j+µ )+tanh(k j µ K j+µ )] s j 1 h +[tanh(k j µ +K j+µ ) tanh(k j µ K j+µ )] s j+1 h, 1 j L, s 0 h = 0, s L+1 h = 0, (64) which can be written as (17). Denoting an eigenvalue of h by E, an the corresponing eigenvector by ψ E, it is seen that there are solutions to (64) of the form s j h (t) = ψ Ej exp(et). (65) These solutions ecay with a relaxation time τ satisfying τ = 1 Re(E). (66) One can see that the eigenvalues of the operator h are real. To see this, one notices that equations (64) are the same as the equations corresponing to the homogeneous solution of (8). So the homogeneous solution to the Ising chain externally riven at ens, is the same as the homogenous solution to the Ising chain with magnetic fiels at bounaries. The evolution equation for the latter satisfies the etaile balance. For any evolution satisfying etaile balance, the eigenvalues of the evolution operator are real. To see this, one notices that the criterion of the etaile balance is ω(b A) = Y A B exp[β(e B E A )], (67) where A an B are two ifferent state, YB A s are real nonnegative numbers (for B A) satisfying YA B = Y B A, (68), 12

13 an E is the energy of the system in the state A. So the matrix Y is Hermitian. Equation (67) means that the evolution matrix H, the components of which are ω(b A) s, isasimilarity-transformeofy. AsY ishermitian,theeigenvalues of Y are real. As H is a similarity-transforme of Y, the eigenvalues of H are the same as the eigenvalues of Y. So the eigenvalues of H are real ([16] for example). The eigenvalues of h are eigenvalues of H as well. So the eigenvalues of h are real. 4 Concluing remarks A one imensional kinetic Ising moel at temperature T, with time varying bounary conitions was stuie, for the case the lattice is semi-infinite. The evolution equation for the expectation values of the spins was investigate. For the case of harmonic bounary conitions, with uniform couplings, exact particular solutions were obtaine for the expectation values of the spins, as well as the total magnetization. The low- an high-frequency behaviors were stuie in more etail. Moels for which the coupling constant is nonuniform were also stuie. Physically, such nonuniform couplings coul arise when either the interaction between spins or the temperature epens on the position. As a specific example, the harmonic solution on a semi-infinite lattice consisting of two homogeneous parts stuie. Finally, it was shown that for a general (nonuniform) lattice, the eigenvalues corresponing to the evolution operator are real. Acknowlegement: This work was supporte by the research council of the Alzahra University. 13

14 References [1] R. J. Glauber; J. Math. Phys. 4 (1963) 294. [2] M. Droz, J. Kamphorst Leal Da Silva & A. Malaspinas; Phys. Lett. A115 (1986) 448 [3] T. Vojta; Phys. Rev E55 (1997) [4] R. B. Stinchcombe, J. E. Santos, & M. D. Grynberg; J. Phys. A31 (1998) 541. [5] C. Gorèche & J. M. Luck; J. Phys. A33 (2000) [6] M. Droz, Z. Rácz, & J. Schmit; Phys. Rev. A39 (1989) [7] C. Gorèche; J. Stat. Mech. (2011) P [8] C. Chatelain; J. Phys. A A36 (2003) [9] M. Khorrami & A. Aghamohammai, 2002 Phys. Rev. E [10] M. Khorrami & A. Aghamohammai, 2008 arxiv: [11] A. Aghamohammai & M. Khorrami; J. Stat. Mech. (2010) P1001. [12] V. Popkov, M. Salerno & G. M. Schutz; Phys. Rev E78 (2008) [13] G. M. Bueni & P. A. Rikvol; Phys. Rev. E78 (2008) [14] A. Praos, L. L. Bonilla & A. Carpio; J. Stat. Mech. (2010) P [15] L. L. Bonilla, A. Praos, & A. Carpio; J. Stat. Mech. (2010) P [16] Geoffrey R. Grimmett & Davi R. Stirzaker; Probability an ranom processes (Oxfor University Press, 2001) thir eition, section