arxiv:cond-mat/ v1 [cond-mat.stat-mech] 17 Sep 1999

Transcription

1 Shape Effects of Finite-Size Scalin Functions for Anisotropic Three-Dimensional Isin Models Kazuhisa Kaneda and Yutaka Okabe Department of Physics, Tokyo Metropolitan University, Hachioji, Tokyo , Japan arxiv:cond-mat/999262v [cond-mat.stat-mech] 7 Sep 999 Macoto Kikuchi Department of Physics, Osaka University, Osaka 56-43, Japan (Received 22 July 999) The finite-size scalin functions for anisotropic threedimensional Isin models of size L L al (a: anisotropy parameter) are studied by Monte Carlo simulations. We study the a dependence of finite-size scalin functions of the Binder parameter and the manetization distribution function p(m). We have shown that the finite-size scalin functions for p(m) at the critical temperature chane from a twopeak structure to a sinle-peak one by increasin or decreasin a from. We also study the finite-size scalin near the critical temperature of the layered square-lattice Isin model, when the systems have a lare two-dimensional anisotropy. We have found the three-dimensional and two-dimensional finite-size scalin behavior dependin on the parameter which is fixed; a unified view of 3D and 2D finite-size scalin behavior has been obtained for the anisotropic 3D Isin models. I. INTRODUCTION Finite-size scalin (FSS) has been increasinly important in the study of critical phenomena [,2]. This is partly due to the proress in the theoretical understandin of finite-size effects, and partly due to the application of FSS in the analysis of simulational results. Recently, more attention has been paid to the universality of FSS functions [3] for both percolation models [4 7]andIsinmodels[8 ]. Ithasbeenshownthatseveral quantities includin distribution functions on various lattices have universal FSS functions by choosin appropriate metric factors. It should be noted that universal FSS functions depend on boundary conditions and the shape of finite systems. In the percolation problem, it was considered until a quite recent time that there exists only one percolatin cluster on two-dimensional (2D) lattices at the percolatin threshold. But the importance of the number of percolatin clusters for anisotropic systems, which was pointed out by Hu and Lin [5], has captured current interest [5, 3]. Aizenman [] derived the upper and lower bounds of the probability for the appearance of n percolatin clusters, W n, for 2D percolation at criticality. His result was later confirmed by Monte Carlo simulations [2]. By use of conformal field theory, Cardy [3] proposed an exact formula for W n at the percolation threshold for the systems with lare aspect ratios [3]. The aspect ratio L /L 2 is an important quantity in the FSS functions for anisotropic 2D systems of size L L 2. The nonuniversal scalin of the low-temperature conductance peak heihts for Corbino disks in the quantum Hall effect has been discussed in terms of the number of the percolatin clusters [4]. For the Isin model, the aspect ratio dependence of Binder parameter [5] at the critical temperature was studied by Kamieniarz and Blöte [6]. Binder and Wan discussed anisotropic FSS [7]. The aspect ratio dependence of the universal FSS functions for Binder parameter and manetization distribution function was discussed by Okabe et al. [] for the 2D Isin model with tilted boundary conditions. For fixed set of the aspect ratio and the tilt parameter, the FSS functions were shown to be universal []. Quite recently, based on the connection between the Isin model and a correlated percolation model, Tomita et al. [8] studied the FSS properties of distribution functions for the fraction of lattice sites in percolatin clusters in subraphs with n percolatin clusters, f n (c), and the distribution function for manetization in subraphs with n percolatin clusters, p n (m). They studied the chane of the structure of the manetization distribution function for the 2D system with lare aspect ratio in terms of percolatin clusters. Almost all results above for the shape effects on FSS functions are for 2D systems except for some works on percolation problem [6,7]. It is quite important to extend these aruments for hiher-dimensional systems. Conformal invariance plays a role in 2D systems [9], but is not so powerful for three-dimensional (3D) systems as for 2D ones. Moreover, there are two limitin cases for anisotropic 3D systems, that is, the one-dimensional(d) limit and the 2D limit. It is interestin to study various types of scalin for anisotropic 3D systems. In this paper, we study the FSS functions for anisotropic 3D systems by Monte Carlo simulations. We are concerned with the ferromanetic Isin model on the L L al simple cubic lattices with the periodic boundary conditions, where a is rearded as the anisotropy parameter. Attention is mainly paid to the effect of anisotropy, and we study the a dependence of the FSS functions for quantities such as the Binder parameter and the manetization distribution function.

2 At the critical temperature for the 3D Isin model, we find that FSS functions for the manetization distribution function chane from a two-peak structure to a sinle-peak one, when the anisotropy parameter a is varied from. This behavior is observed both for the systems with 2D anisotropy and for those with D anisotropy. When a system has a lare 2D anisotropy, we find that finite-size systems with fixed anisotropy parameter show ood FSS behavior as 3D systems near the critical temperatures for the layered square-lattice Isin models. In contrast, when we fix a number of layers and apply the FSS analysis for layered systems, these systems are scaled as 2D ones. We oranize the rest of the paper as follows. In Sec. 2, we define the FSS functions and describe the quantities we treat in this paper. In Sec. 3, we present our simulational results. In Sec. 4, our results are summarized and discussed. II. FINITE-SIZE SCALING If a quality Q has a sinularity of the form Q(t) t ω (t = T T c ) near the criticality t =, then the correspondin quantity Q(L, t) for finite systems with the linear size L has a followin scalin form, Q(L,t) L ω/ν f(tl /ν ), () where ν is the correlation-lenth critical exponent and f(x) is the scalin function. Of course, the corrections to FSS are not neliible for smaller L. The FSS is also applicable to the distribution function of Q. At the criticality t =, the FSS function has a followin form, p(q;l,t = ) L ω/ν F(QL ω/ν ). (2) In this paper we focus on the Binder parameter and the manetization distribution function. The Binder parameter [5] is iven by = 2 (3 m4 m 2 2), (3) and serves as a measure of the non-gaussian nature of the distribution function. The FSS functions for these quantities are iven by (L,t) (tl /ν ), (4) p(m;l,t = ) L β/ν p(ml β/ν ), (5) where β is the manetization exponent. simple cubic lattices with different values of L and a. We use the periodic boundary conditions x6x8 24x24x2 32x32x6 6x6x4 24x24x6 32x32x T a=/4 6x6x8 24x24x2 32x32x6 6x6x4 24x24x6 32x32x8 a=/2 a= 5 5 (T T c ) L.6 (a) 6x6x6 (b) 6x6x6 FIG.. (a) Temperature dependence of for several lattices with a=, /2, and /4. The temperature is represented in unit of J. (b) Plot of as a function of (T T c)l /ν. First, we show the results for a, that is, a 2D anisotropic case. We plot the temperature dependence of the Binder parameter for a=, /2, and /4 of various sizes in Fi. (a). Error-barsare within the size of marks unless specified. We ive the FSS plots in Fi. (b); is plotted as a function of (T T c )L /ν, where L is iven by (L L al ) /3. For T c, β and ν, we use the numericallyestimatedvalues forthe 3DIsinmodel [2], that is, T c = 4.54, β =.32 and ν =.625. From now on, we represent the temperature in unit of J. From Fi. (b) we see that the FSS functions have a dependence in 3D systems, which is similar to the case of 2D systems. III. RESULTS We use the Metropolis Monte Carlo method to simulate the ferromanetic Isin model on L L al 2

3 p(m) L x6x6 a= 6x6x8 24x24x2 32x32x6 a=/2 a=/4 6x6x4 24x24x6 32x32x m L.52 FIG. 2. p(m)l β/ν at T = T c as a function of ml β/ν for several lattices with a=, /2, and /4. Wecanetmoreinformationontheadependencefrom the manetization distribution function p(m). In Fi. 2, we show the scalin plot of the manetization distribution function p(m) at T = T c for various sizes with different anisotropy parameters (a ). From the fiure we find ood FSS behavior and also a lare a dependence. By decreasin the anisotropy parameter a from, the FSS functions for the manetization distribution function chane from a two-peak structure to a sinle-peak one. Next turnto the caseofa, that is, a Danisotropic case. We plot the temperature dependence of for a=, 2 and 4 of various sizes in Fi. 3(a). The FSS plots are iven in Fi. 3(b). We see similar behavior of the FSS functions for a, that is, take smaller values if we chane a from. We show the a dependence of the FSS functions of the manetization distribution functions at T = T c for a in Fi. 4. By increasin a from, FSS functions for p(m) chane from a two-peak structure to a sinle-peak one aain..5 a=4 2x2x24 6x6x32 24x24x48 2x2x48 6x6x64 24x24x96 a=2 a= 5 5 (T T c ) L.6 (b) 6x6x6 FIG. 3. (a) Temperature dependence of for several lattices with a=, 2, and 4. (b) Plot of as a function of (T T c)l /ν. If we compare the results of Fis. (b) and 3(b), and also those of Fis. 2 and 4, it seems that there could be a set of the anisotropy parameters a < and a > which ive the same FSS functions. But it needs more careful calculations to clarify this universality. p(m) L x2x2 6x6x6 a= 2x2x24 6x6x32 24x24x48 a=2 a=4 2x2x48 6x6x64 24x24x96 (a) 6x6x m L.52 FIG. 4. p(m)l β/ν at T = T c as a function of ml β/ν for several lattices with a=, 2, and x2x24 6x6x32 24x24x48 2x2x48 6x6x64 24x24x T Let us consider the oriin of the a dependence of FSS functions. This problem is related to the multiple percolatin clusters. The connection between critical phenomena of spin models and percolation problems has been studied for a lon time [2 23]. Quite recently, Tomita et al. [8] have used the cluster formalism to investiate the percolatin properties of the 2D Isin model. They elucidated that the existence of several percolatin clusters for anisotropic finite systems leads to the chane of the structure of p(m). That is, the combination of up-spin clusters and down-spin clusters ives the contribution of m in p(m) for anisotropic case. The a dependence of 3

4 p(m) for the 3D systems can be understood in the same way as the 2D systems. We can apply this arument both cases for 2D anisotropy (a < ) and D anisotropy (a > ). The value of Binder parameter at T = T c, of course, reflects upon the structure of the distribution function p(m)..5 s=4 s=6 s=8 c.5 64x64x4 96x96x4 28x28x4 64x64x6 28x28x T 64x64x8 96x96x8 28x28x8 FIG. 6. Temperature dependence of for several systems. The numbers of layers are 4, 6, and a FIG. 5. at T = T c, c, as a function of a. We use a semi-loarithmic scale. a=/6 In order to see the anisotropy parameter dependence more clearly, we plot the a dependence of the Binder parameter at criticality, c, in Fi. 5. Here, we use the semi-loarithmic scale for plot. We add several data other than those iven in Fis. and 3. For lare enouh a or small enouh a, the corrections to FSS become larer. From Fi. 5, we find that c takes the maximum value at a = and decreases radually in both directions, a < and a >. For very lare a or very small a, c tends to vanish, which indicates that the distribution function p(m) approaches the Gaussian distribution. This behavior is consistent with above-mentioned multiple-percolatin-cluster arument. Such an a dependence of at criticality for D anisotropy of 2D Isin model was already discussed for 2D Isin model [6]. It is also interestin to consider the relation to the layered square-lattice Isin model for the case of lare 2D anisotropy. The critical temperatures and the shift exponent for the layered square-lattice Isin model were studied by Kitatani et al. [24]. In Fi. 6, we ive the temperature dependence of for several lattices, where the number of layers, s, is 4, 6, and 8. If we make the system size lare with fixin s, the system becomes the layered square-lattice Isin model. We can estimate the s-dependent critical temperature T cs from the crossin point of the curves for fixed s in Fi. 6. Our estimates, T c4 = (4), T c6 = 4.37(4), and T c8 = (4), are consistent with those of Ref. [24]..5 64x64x4 28x28x8 5 5 (T T cs ) L.6 FIG. 7. Scalin plot of for fixed a(= /6). The scaled variable is (T T cs)l /ν for each s. The values of T cs are T c4 = , T c6 = 4.37, and T c8 = We use the 3D exponent, ν =.625. Pickin up the data for a = /6 from Fi. 6, we plot as a function of (T T cs )L /ν in Fi. 7. Here we treat T T cs for each s, and ν is chosen as.625, that is, the 3D value. We have ood FSS behavior in Fi. 7, but this FSS plot is different from that iven in Fi. (b), where T T c has been treated. We should note that T T cs can be rewritten as T T cs = (T T c ) + (T c T cs ), and T c T cs s λ. Here λ is the shift exponent and is shown to be /ν 3D [24]. Thus, for fixed a, we have 3D FSS for T T cs. In this case, the value of c at T = T cs is.92, which is nothin but the c value for the isotropic 2D Isin model. It is different from.4 for c value of anisotropic 3D Isin model for a = /6, which is iven in Fi. 5. In Fi. 8, we also show the scalin plot of the manetization distribution function p(m) at T = T cs for each s, fixin a as /6. We have ood FSS behavior 4

5 usin 3D exponents. There are two distinct peaks in p(m), which is a typical character of p(m) at T c for the isotropic 2D Isin model. In contrast, if we pick up the data for fixed s, for example, s = 6, these systems are expected to be scaled by usin 2D critical exponents. Since we fix s(= al ), the scaled variable is iven by (T T cs )L /ν. Usin the 2D exponent, ν 2D =, we et very ood scalin plot, shown in Fi. 9. Of course, the critical value c at T = T c6 is.92, that is, the c value for the isotropic 2D system. p(m) L x64x4 28x28x8 a=/ m L.52 FIG. 8. Scalin plotofp(m)at T = T cs for fixeda(= /6); p(m)l β/ν as a function of ml β/ν. The values of T cs are T c4 = , T c6 = 4.37, and T c8 = We use the 3D exponents, β =.32 and ν =.625. s=6 IV. SUMMARY AND DISCUSSIONS To summarize, we have studied the FSS functions for anisotropic 3D Isin model by Monte Carlo simulations. The anisotropy parameter (a) dependence of the FSS functions for the Binder parameter and the manetization distribution function has been investiated. We have shown that the manetization distribution functions chane from two-peak structures to sinle-peak ones as a increases or decreases from. This chane of distribution functions may be attributed to the combination of up-spin percolatin clusters and down-spin percolatin clusters for anisotropic systems, which is the same as the case of 2D systems [8]. We should also note that the distribution function p(m) approaches the Gaussian distribution for lare anisotropy, which is another indication of the multiple-percolation-cluster arument. It is interestin to calculate the probability for the appearance of n percolatin clusters, W n, for the anisotropic 3D Isin model and to make cluster analysis, which will be left to a future study. Moreover, we have considered the relation to the layered square-lattice Isin models for lare 2D anisotropy. We have obtained 3D FSS behavior near the critical temperature of the layered square-lattice Isin models, T cs, for fixed a. In contrast, when we fix the number oflayers, we have shown 2D FSS behavior near T cs. Thus, we have obtained a unified view of 3D and 2D FSS behavior for the anisotropic 3D Isin models. We have studied the shape effects of FSS functions for 3D Isin models. FSS functions also depend on boundary conditions. Extension of the 2D results for various boundary conditions, for example, tilted boundary conditions [], to 3D systems will be interestin. This study is now in proress. ACKNOWLEDGMENTS.5 64x64x6 28x28x6 2 2 (T T c6 ) L We would like to thank C.-K. Hu and Y. Tomita for valuable discussions. The computation in this work has been done usin the facilities of the Supercomputer Center, Institute for Solid State Physics, University of Tokyo. This work was supported by a Grant-in-Aid for Scientific Research from the Ministry of Education, Science, Sports and Culture, Japan. FIG. 9. Scalin plot of for fixed s(= 6). The scaled variable is (T T c6)l /ν, and T c6 = We use the 2D exponent, ν =. Electronic address: kkaneda@phys.metro-u.ac.jp Electronic address: okabe@phys.metro-u.ac.jp [] Fisher M E 97 Proc. 97 E. Fermi Int. School of Physics vol 5, ed M S Green (New York: Academic) p 5

6 [2] Cardy J L (ed) 988 Finite-size Scalin (New York: North-Holland) [3] Privman V and Fisher M E 984 Phys. Rev. B [4] Hu C-K, Lin C-Y and Chen J-A 995 Phys. Rev. Lett and 2786E [5] Hu C-K and Lin C-Y 996 Phys. Rev. Lett [6] Lin C-Y, Hu C-K and Chen J-A 998 J. Phys. A 3 L [7] Lin C-Y and Hu C-K 998 Phys. Rev. E [8] Okabe Y and Kikuchi M 996 Int. J. Mod. Phys. C [9] Wan F-G and Hu C-K 997 Phys. Rev. E [] Okabe Y, Kaneda K, Kikuchi M and Hu C-K 999 Phys. Rev. E [] Aizenman M 997 Nucl. Phys. B [FS] [2] Shchur L N and Kosyakov S S 997 Int. J. Mod. Phys. C [3] Cardy J L 998 J. Phys. A 3 L5 [4] Cooper N R, Halperin B I, Hu C-K and Ruzin I M 997 Phys. Rev. B [5] Binder K 98 Z. Phys. B 43 9 [6] Kamieniarz G and Blöte H W J 993 J. Phys. A 26 2 [7] Binder K and Wan J-S 989 J. Stat. Phys [8] Tomita Y, Okabe Y and Hu C-K 999 Phys. Rev. E [9] Cardy J L 987 Phase Transitions and Critical Phenomena vol, ed C Domb and J L Lebowitz (New York: Academic) p 55 [2] Ferrenber A M and Landau D P 99 Phys. Rev. B [2] Kasteleyn P W and Fortuin C M 969 J. Phys. Soc. Jpn. Suppl. 26 Fortuin C M and Kasteleyn P W 972 Physica [22] Conilio A and Klein W 98 J. Phys. A [23] Hu C-K 984 Phys. Rev. B and 59 [24] Kitatani H, Ohta M and Ito N 996 J. Phys. Soc. Jpn