arxiv:cond-mat/ v1 [cond-mat.stat-mech] 12 Jun 2006

Transcription

1 Continuous phase transitions with a convex dip in the microcanonica entropy Hans Behringer Fakutät für Physik, Universität Bieefed, D Bieefed, Germany arxiv:cond-mat/66283v1 [cond-mat.stat-mech] 12 Jun 26 Miche Peiming Institut für Theoretische Physik I, Universität Erangen-Nürnberg, D 9158 Erangen, Germany The appearance of a convex dip in the microcanonica entropy of finite systems usuay signas a first order transition. However, a convex dip aso shows up in some systems with a continuous transition as for exampe in the Baxter-Wu mode and in the four-state Potts mode in two dimensions. We demonstrate that the appearance of a convex dip in those cases can be traced back to a finitesize effect. The properties of the dip are markedy different from those associated with a first order transition and can be understood within a microcanonica finite-size scaing theory for continuous phase transitions. Resuts obtained from numerica simuations corroborate the predictions of the scaing theory. PACS numbers: 5.5.+q, 5.7.Jk, 75.1.Hk I. INTRODUCTION The microcanonica anaysis of phase transitions has gained growing interest in recent years due to the inequivaence of the canonica and the microcanonica descriptions of finite systems with short-range interactions as we as of systems with ong-range interactions [1, 2, 3, 4, 5, 6, 7, 8, 9, 1, 11, 12, 13, 14, 15, 16, 17]. Whereas the discrimination of the order of a phase transition from finite-size data can be a difficut task in the canonica ensembe, microcanonica quantities usuay yied cear signatures which revea the continuous or the discontinuous character of a phase transition [2]. In sma microcanonica systems which exhibit a continuous transition in the thermodynamic imit the character of the transition is signaed by typica features of symmetry breaking, as e. g. the abrupt onset of a non-zero order parameter when a pseudo-critica point is approached from above. These features are accompanied by singuar physica quantities [18, 19, 2, 21, 22, 23, 24, 25, 26, 27]. The associated singuarities, however, are not the consequence of a non-anaytic microcanonica entropy and can therefore be characterized by cassica critica exponents [19, 22]. Simiary, intriguing features are aso reveaed in the microcanonica entropy of sma systems with a discontinuous transition in the infinite voume imit [5, 28, 29, 3, 31, 32, 33]. A typica back-bending of the microcanonica caoric curve is observed, eading to a negative heat capacity. This property, which has aso been observed in recent experiments [34, 35], has been used as a too for the determination of the first order character of phase transitions in microcanonica systems. In the canonica context the first order character is then signaed by a doube-peak structure in certain histograms (e. g. [36, 37, 38]). However, some systems (as for exampe the Baxter- Wu mode [39, 4]), with a continuous phase transition in the thermodynamic imit, have been reported to show a simiar doube-peak structure in canonica histograms or, equivaenty, a negative microcanonica specific heat [41, 42]. It foows from this observation that the mere existence of a convex dip in the entropy is not sufficient for the identification of a discontinuous transition [43]. In this work we examine more cosey the finite-size behaviour of the microcanonica specific heat and of the convex dip in the entropy function of this kind of systems, thereby taking the Baxter-Wu mode and the four-state Potts mode in two dimensions as exampes. In particuar we demonstrate that the microcanonica finite-size scaing theory presented in Reference [26] for continuous phase transitions reates the appearance of these features to a pecuiar finite-size effect. In addition we deduce properties of the dip in the microcanonica entropy which we test in numerica simuations of both the Baxter-Wu and the four-state Potts mode. The rest of the artice is organized as foows. In Section II a brief introduction to the microcanonica anaysis of physica systems is given. We particuary discuss the properties of the microcanonica specific heat both for systems with continuous and discontinuous phase transitions in the infinite voume imit. The microcanonica finite-size scaing theory of Reference [26] is used in Section III to deduce the scaing properties of the convex intruder in the microcanonica specific heat of systems undergoing a continuous phase transition in the thermodynamic imit. In Section IV we present numerica resuts both for the Baxter-Wu and for the two-dimensiona fourstate Potts mode which beong to the same universaity cass. These numerica resuts corroborate the findings obtained from the finite-size scaing theory as the properties of the convex dip are consistent with the theoretica predictions. Finay, Section V summarizes our resuts. II. SIGNATURES OF PHASE TRANSITIONS IN THE MICROCANONICAL SPECIFIC HEAT In finite systems the thermostatic quantities of a system undergoing a phase transition in the thermodynamic

2 2 imit show characteristic features. The canonica specific heat of finite systems, for exampe, is a reguar function that exhibits a maximum at a certain temperature. The continuous or discontinuous character of the phase transition is reveaed by the size-dependence of the increase of the specific heat maximum which is we understood within canonica finite-size scaing theory [44, 45]. In the case of a continuous phase transition the increase of the specific heat maximum c max with the system size L is given by c max L α/ν where α describes the divergence of the specific heat of the infinite system when approaching the critica temperature (here and in the foowing we ony consider the case α > ), whereas ν is the usua correation ength critica exponent. For a discontinuous phase transition, however, one finds c max L d [46] where d is the dimensionaity of the system. The thermodynamic behaviour can aso be investigated in the microcanonica ensembe which describes systems which are isoated from any environment [1, 2]. The starting point in the microcanonica anaysis is the density of states Ω which, for a discrete spin system, measures the degeneracy of the macrostate (E), i.e. Ω(E) is the number of configurations (or microstates) with energy E. For a d-dimensiona system with inear extension L physica quantities are then cacuated from the microcanonica entropy density s(e, L 1 ) := 1 L d n Ω(Ld e, L 1 ) (1) which depends on the energy density e := E/L d and the system size L. In the foowing the dependence on the system size is indicated by the inverse system size := L 1. Note that the specific entropy aso depends on the boundary conditions which have to be specified for finite systems. This dependence, however, is not indicated in the notation here. The inverse microcanonica temperature β µ shows up as the conjugate variabe to the natura variabe energy e and is defined by β µ (e, ) := d s(e, ). (2) de The corresponding microcanonica specific heat is given by ( ) 1 c(e, ) = (β µ (e, )) 2 dβµ (e, ). (3) de Simiar to the canonica case the finite-size behaviour of the microcanonica specific heat shows aso certain characteristics if the system has a phase transition in the infinite attice. In the case of a continuous transition the entropy density is concave everywhere. The microcanonica specific heat (3) is thus positive and dispays a maximum whose vaue increases with growing system size L [23, 26]. As the derivative of β µ (e, ) is the second order derivative of the entropy, it measures the curvature of the entropy density. This curvature is negative and has a maximum at a certain energy which can be identified as the pseudo-critica energy e pc of the finite system. The vaue of the derivative of β µ at this pseudo-critica energy scaes ike β µ(e pc, ) := d de β µ(e, ) e=epc = d2 de 2 s(e, ) e=epc αε/νε (4) in the regime. The microcanonica critica exponent α ε is thereby reated to the canonica critica exponent α by α ε = α/(1 α) [18, 19]. Generay, a canonica critica exponent κ transates to the microcanonica exponent κ ε = κ/(1 α) if α is positive, whereas for non-positive α the critica exponents of the microcanonica description are identica to the canonica ones [47]. Reation (4) is obtained from a decomposition of the microcanonica entropy into a reguar part and a singuar part which obeys a certain homogeneity reation in the vicinity of the pseudo-critica point. Note that the size-dependence of the increase of the specific heat maximum aows the determination of the corresponding microcanonica critica exponent of the specific heat of the infinite system. For more detais see Section III and in particuar Reference [26]. Note aso that other microcanonica quantities ike the spontaneous magnetization and the susceptibiity aso exhibit typica features of continuous phase transitions in finite systems [19, 2, 21, 22, 25, 27]. These aspects are however not studied in the present work. The behaviour of the specific heat is competey different if the system undergoes a discontinuous phase transition. In this case the microcanonica entropy density of a finite system has a convex dip which eads to a negative microcanonica specific heat for a certain energy interva [5, 28, 29, 3, 31, 32, 33]. This is different in the canonica ensembe where the specific heat, which is reated to the variance of the energy, is aways positive for any finite system. Of course in the thermodynamic imit the convex dip in the microcanonica entropy density disappears in systems with finite-range interactions and the specific heat is positive everywhere. The curvature properties of the entropy density of a finite system aow the determination of an inverse pseudo-transition temperature β t by the doube-tangent (or Maxwe) construction. This is depicted schematicay in Figure 1. The discontinuous character of the transition is reveaed by a certain scaing behaviour of the convex dip. Definig the auxiiary function s(e, ) := s(e, ) β t ()e, a size-dependent width ω() and depth δ() of the dip can be measured, as shown in Figure 1. With increasing system size, i. e. for the regime, the depth scaes ike δ() and the width ike ω() ω [48, 49, 5]. Here ω is the ength of the inear section of the function s of the infinite system and is in fact the atent heat of the transition. The asymptotic behaviour inear in and the appearance of a non-zero atent heat ω signa the discontinuous character of the transition in the infinite system. This can be iustrated by investigating the three-state Potts mode in three dimensions. The genera q-state Potts mode is

3 3 s(e) δ () e..5.1 s(e).25 ω ().2 δ ω FIG. 1: Schematic depiction of the convex dip in the entropy of a finite system whose corresponding infinite system undergoes a discontinuous phase transition. The doube-tangent with sope β t is indicated by a dashed ine (top). For the function s = s β t e the definitions of the width ω and of the depth δ of the dip are shown (bottom). defined by the Hamitonian H(σ) = ij e δ σi,σ j (5) where the sum runs over nearest neighbour sites of the attice and the spins σ i can take on the vaues σ i = 1,...,q [4, 51, 52]. The couping constant is set to unity. The three-state mode in three dimensions undergoes a (reativey weak) first order transition and therefore one expects that the characteristic quantities of the convex dip evove ineary in the inverse system size (compare [53] for a canonica and [29] for a microcanonica investigation of the mode). This is indeed confirmed by numerica data as shown in Figure 2. III. MICROCANONICAL FINITE-SIZE SCALING THEORY In this Section the impications of the microcanonica finite-size scaing theory [26] are investigated for systems with a continuous phase transition where α ε ν ε (canonicay α ν). Note that the condition α ε ν ε corresponds to a rather restricted cass of mode systems. Nevertheess important spin modes ike the Baxter-Wu FIG. 2: The depth δ() and width ω() of the convex dip in the microcanonica entropy of the three-state Potts mode in three dimensions as a function of the inverse system size. The dashed ines are inear fits to the data. Error bars are comparabe to the size of the symbos. and the four-state Potts mode in two dimensions have critica exponents α ε = ν ε and therefore beong to this cass. These modes are coser studied in the next Section. Other modes are reported in the iterature to have exponents satisfying α ε > ν ε [54, 55]. Consider the entropy density s(e, ) of a finite d- dimensiona microcanonica system near the pseudocritica energy e pc. Reca that e pc is defined by the maximum of the second derivative of s. The microcanonica finite-size scaing theory for a system with a continuous transition [26] assumes that the entropy can be spit up into two parts, namey into a so-caed reguar part s r and into a part s s which is caed singuar as it becomes nonanaytic in the infinite voume imit. This is in contrast to the reguar part which is anaytic aso in the infinite system. The singuar part has to satisfy the scaing reation s s (ε, ) = d φ( 1/νε ε) (6) with ε := e e pc and φ being a scaing function with no further expicit size-dependence. This decomposition s(ε, ) = s r (ε, ) + s s (ε, ) is ony vaid for asymptoticay arge system sizes and in the regime 1/νε ε, in particuar for energies cose to the pseudo-critica one. Note that the singuar part and thus the function φ in reation (6) is aso anaytic for finite systems and must fufi the required scaing property which states that it is a homogeneous function. Due to the anayticity of the microcanonica entropy the decomposition together with the required scaing reation for the singuar part states

4 4 that the entropy is of the form s(ε, ) = s ()+β pc ()ε+ 1 2 A 2()ε A 4()ε (7) for sma ε [26]. Here β pc () = β µ (e pc, ). The coefficients of the second and fourth order terms are of the form and A 2 () = B 2 αε νε + C 2 () (8) A 4 () = B 4 αε 2 νε + C 4 () (9) in terms of the critica exponents of the microcanonica system. The coefficient A 4 has to be negative, otherwise higher order terms in the energy deviation ε have to be taken into account. The size-dependence contained in the coefficients C 2 and C 4 stems from the reguar part s r of the entropy density whereas the dependencies invoving the critica exponents α ε and ν ε have their origin in the singuar part s s (or equivaenty in the scaing function φ in reation (6)) of the entropy. In particuar, the coefficient B 2 has to be negative in the case of a continuous transition and, due to the reguarity of s r, the coefficient C 2 is of the form C 2 () = v 1 + v (1) for sma whereas the coefficient C 4 has the expansion C 4 () = z +z in the inverse system size (see reference [26] for further detais). Note that α ε has to satisfy α ε 2 for C 4 to be the subdominant term in the coefficient A 4. This is assumed in the foowing considerations. In certain mode systems the critica exponents satisfy the reation α ε = ν ε. Exampes incude the Baxter-Wu mode and the four-state Potts mode in two dimensions. In this situation the eading behaviour of the coefficient A 2 () in reation (7) is the inear term in the inverse system size, i. e. A 2 () a for where the ampitude a is the sum B 2 +v 1. If the ampitude a is positive, i. e. v 1 is positive and arge enough compared to the moduus of the negative ampitude B 2, the entropy density deveops a convex dip for a finite system sizes. Note that the corresponding condition A 2 () > is a necessary condition for the appearance of a convex dip in the microcanonica entropy. Thus, in the framework of the microcanonica finite-size scaing theory a convex dip can appear even so the transition in the infinite system is continuous. Furthermore, from expression (7) one can then deduce scaing reations for the width and the depth of the convex dip in such a case. Appying the doube-tangent construction to the entropy (7) with α ε = ν ε and a positive A 2 one finds that and ω() = 2 δ() = (A 2()) 2 A 4 () A 2 () 1/αε (11) A 4 () 1+2/αε (12) for sma. This scaing behaviour of the convex dip in the microcanonica entropy is therefore different from the inear behaviour in the inverse system size encountered for discontinuous transitions and discussed in Section II [54, 55]. For a negative ampitude a, on the other hand, the signatures of a continuous transition with the scaing behaviour β µ(e pc, ) A 2 () (= αε/νε ) are observed. Whether the predicted convex dip reay appears or whether one encounters typica features of continuous transitions thus depends on the sign and reative size of the ampitudes entering the scaing form (7) of the entropy near the pseudo-critica energy. Let us end this Section by considering the situation where α ε is stricty arger than ν ε. In this case the eading behaviour of the coefficient A 2 is sti inear in provided the ampitude v 1 is non-zero. For positive v 1 one then has again a convex dip with the scaing reations ω() 1/2 (αε 2)/(2νε) and δ() 2 (αε 2)/νε for its characteristic extensions. For negative v 1 instead one observes features of a continuous transition with the asymptotic behaviour β µ (e pc, ), which does not aow to uncover the actua critica exponents from the eading size-dependence of the evoution of the microcanonica specific heat, in contrast to the corresponding situation for α ε ν ε. In concusion, a convex dip in the microcanonica entropy is in genera possibe for systems with a continuous phase transition if the critica exponents satisfy α ε ν ε. For systems with α ε > ν ε such a dip shows up if the ampitude v 1 of the inear term in the expansion of C 2 in (1) is positive. In a situation with α ε = ν ε the ampitude v 1 has to satisfy the additiona requirement that v 1 > B 2. In both cases we end up with the necessary condition A 2 () >, see Eq. (8). IV. MICROCANONICAL SPECIFIC HEAT OF THE BAXTER-WU AND OF THE FOUR-STATE POTTS MODEL A. The Baxter-Wu mode The Baxter-Wu mode is a cassica spin mode whose free energy can be cacuated exacty in zero magnetic fied [39, 4]. The mode is defined on a two-dimensiona trianguar attice, with the Hamitonian H(σ) = ijk σ i σ j σ k (13) where the strength of the couping constant has been set to unity. The spins σ i are Ising spins with the possibe vaues σ i = ±1. The sum in (13) extends over a trianges of the trianguar attice (here denoted by ijk with i, j and k specifying the vertices of the trianges). The system has a continuous phase transition in the infinite system at the critica temperature T c = 2/ n(1 + 2). The critica exponents which characterize the singuar

5 5 behaviour of physica quantities near the critica point are known exacty. For the specific heat and the correation ength one has for exampe α = ν = 2/3 or α ε = ν ε = 2 microcanonicay. In this Section the microcanonica specific heat of the Baxter-Wu mode is investigated numericay [56]. To this end the microcanonica entropies of systems with inear sizes ranging from 16 to 96 are cacuated using the very efficient transition observabe method [2, 27] for both open and periodic boundary conditions. As the Baxter-Wu mode exhibits a continuous phase transition in the infinite attice, one expects at first sight that the microcanonica specific heat shows a maximum with the associated increase for increasing system sizes. For finite systems with open boundary conditions this is indeed the case. The derivative β µ(e pc, ) of the microcanonica temperature at its maximum vaue, i. e. at the pseudo-critica energy, is negative and its moduus decreases, eading to an increasing specific heat maximum when the system size is increased. For the Baxter-Wu mode the critica exponents α ε and ν ε are equa and therefore one expects from (4) and the discussion in Section III to observe the scaing reation β µ(e pc, ) for sma. This expectation is indeed confirmed by the data shown in Figure 3. dβ µ (e pc, ) /de FIG. 3: The derivative β µ(e pc, ) = dβ µ(e pc, )/de of the Baxter-Wu mode as a function of the inverse system size for open boundary conditions. The errors are approximatey as arge as the symbo size. The dashed ine shows a inear fit. The data extrapoate to zero in the imit, corresponding to a divergent specific heat in the infinite system. The situation is different for finite Baxter-Wu systems with periodic boundary conditions. In Figure 4 the function s is shown for a finite system with inear extension L = 6. A convex dip is ceary visibe. This observation apparenty suggests at first sight that the mode seems to undergo a discontinuous transition in contradiction to the known continuous character of the phase transition. Note that recenty a doube-peak structure in canonica histograms of the energy was reported in the iterature [42]. In the previous Section III we have shown that the microcanonica finite-size scaing theory permits the appearance of such a convex dip for the case α ε = ν ε. With the reations (12) and (11) and the known vaue α ε = 2 s(e) - β t e e FIG. 4: The function s for the finite Baxter-Wu system with L = 6 and periodic boundary conditions (here β t =.43945). Note that s is ony obtained from a numerica simuation up to an additive constant. The data show a doube-peak structure typica of finite systems with a discontinuous transition in the macroscopic imit. for the Baxter-Wu mode we expect to observe the scaing behaviour ω() for the width and δ() 2 for the depth of the convex dip. These predictions are indeed confirmed by our numerica data obtained for the Baxter-Wu mode as shown in Figure 5. The fu ines are fitted functions containing the expected asymptotic behaviour in the inverse system size. For the dashed ines the eading correction term coming from the reguar part s r of the microcanonica entropy (which is proportiona to 3 for the depth and to 3/2 for the width) has been incuded. In agreement with the continuous character of the phase transition the dip scaes away in the macroscopic imit. The size-dependences of both the depth and the width of the convex dip foow the theoretica predictions. This is especiay cear from the inset of Figure 5 where we compare ω() with the expected behaviour in a doube-ogarithmic pot. This agreement gives strong evidence to the genera scaing theory presented in the previous Section III. It shoud be stressed again that the different behaviour of the microcanonica entropy of the Baxter-Wu mode for different boundary conditions is reated to different properties of the expansion coefficients in (7). In particuar the coefficient A 2 seems to be deepy affected by the boundary conditions. A microscopic expanation of this observation is however sti acking. The findings of this section have direct impications for the behaviour of canonica energy histograms, which have been discussed for the Baxter-Wu mode in Reference [42]. The energy histograms are directy reated to the canonica distribution function P β (e) which is the probabiity of finding a microstate with energy density e at the inverse canonica temperature β. In terms of the microcanonica entropy this distribution function is basicay given by P β (e) exp ( βl d e L d s(e, ) ). Thus for the inverse pseudo-transition temperature β t obtained from the Maxwe construction appied to the microcanonica entropy the canonica energy histogram exhibits two peaks at equa height with the peaks appearing at the same

6 6.2 B. The four-state Potts mode δ () ω () n(ω) n() FIG. 5: The depth δ() and width ω() as functions of the inverse system size for Baxter-Wu systems with periodic boundary conditions. The errors are of the order of the symbo size. The fu ines show fits to the theoreticay predicted behaviour (11) and (12) using the six argest system sizes. The dashed ines incude in addition the first correction term coming from the reguar part s r of the microcanonica entropy. The inset shows a doube-ogarithmic pot of the width together with a straight ine of sope 1/2. From the six argest systems one gets a sope of.481(3). energies as the peaks in s = s β t e. It foows that in the canonica histograms the peaks at the energies e 1 and e 2 approach each other according to the scaing aw e 1 e 2. An important quantity which infuences the kinetics of first order phase transitions is the interface tension Σ defined by 2β t Σ = 1 L d 1 n P β t (e max ) P βt (e min ) (14) where the doube-peaked distribution P βt has one of the two peaks of equa height at energy e max whereas the minimum between the two peaks appears at the energy e min [57]. From the above discussion one gets Σ δ()/ and thus the interface tension vanishes in the imit of asymptoticay arge systems. Consequenty a coexistence between an ordered phase and a disordered phase at temperature β t is hardy detectabe in this case athough the canonica distribution function P βt exhibits a doube-peak structure. For moderatey sma systems, however, the dip represents sti a considerabe barrier between the order phase and the disordered one. In some sense one can speak of a phase coexistence for not too arge finite systems. The four-state Potts mode in two dimensions beongs to the same universaity cass as the Baxter-Wu mode. The character of the phase transition of the q-state Potts mode in two dimensions changes exacty at q = 4 from continuous for ow vaues of q to discontinuous for arge q. This makes the four-state Potts mode particuary interesting. Due to the presence of margina scaing fieds for q = 4 the singuar behaviour of the four-state Potts mode is not simpy described by pure power aws, but physica quantities ike the specific heat or the susceptibiity acquire additiona ogarithmic factors [58, 59]. Logarithmic corrections then aso show up in the finite-size scaing form of the specific heat [44]. This is in contrast to the Baxter-Wu mode where these marginaity effects are absent and the singuar behaviour of the specific heat does not exhibit ogarithmic corrections [58, 6]. dβ µ (e pc, ) /de FIG. 6: The derivative β µ(e pc, ) = dβ µ(e pc, )/de of the fourstate Potts mode as a function of the inverse system size for open boundary conditions. The errors are approximatey of the order of the symbo size. The dashed ine shows a fit with an asymptotic inear term and corrections incuding a ogarithmic term. The data extrapoate to zero for the imit, yieding a divergent specific heat in the infinite system. The inset shows a fit without the ogarithmic correction term (see main text for further detais). Using the same numerica technique as for the Baxter- Wu mode, we have computed the microcanonica entropy density of the four-state Potts mode on a square attice with open and periodic boundary conditions. For open boundary conditions system sizes ranging from L = 12 to L = 8 have been considered, whereas for periodic boundary conditions system sizes from L = 14 to L = 5 have been simuated. As shown in Figure 6 the derivative β µ(e pc, ) of the microcanonica temperature evauated at the pseudocritica energy is found to be negative for open boundary conditions. As its moduus decreases as a function of the inverse system size, the resuting specific heat maximum increases when the system size is increased, thus yieding the typica finite-size features of a continuous transition in the thermodynamic imit. The critica exponents α ε and ν ε being equa for the four-state Potts mode, one expects (see Eq. (4)) the scaing reation β µ (e pc, ) for

7 7 sma. However, as Figure 6 demonstrates, this asymptotic regime is not yet fuy reached for the system sizes considered here. Presumaby, this has its origin in the presence of ogarithmic correction terms. Nevertheess, the data are consistent with the expectation of a eading inear dependence in the inverse system size in the asymptotic imit. To demonstrate this a function of the form (β µ (e pc, )) fit = a + b/ n() + c 2 has been fitted to the data (dashed ine in Figure 6). To get a rough impression of the presence of a ogarithmic correction, a fit function with correction terms b c 1 3, i. e. without a ogarithmic term but with the same number of fit parameters, has aso been fitted to the numerica data, yieding a much poorer resut (see inset of Figure 6). For the Potts mode with periodic boundary conditions our findings are in quaitative agreement with the considerations of Section III and with the resuts obtained for the Baxter-Wu mode with periodic boundaries. The microcanonica entropy dispays the same convex dip for finite systems. Note that a doube-peak structure in energy histograms obtained from canonica Monte Caro simuations was aready reported in the iterature for this mode [41]. The scaing behaviours of the depth and the width of this dip are again consistent with the predicted scaing reations (11) and (12) as depicted in Figure 7, even so we were forced to incude a ogarithmic correction term in addition to the eading correction term coming from the reguar part. δ () ω () FIG. 7: The depth δ() and width ω() as functions of the inverse system size for four-state Potts systems with periodic boundary conditions. The errors are approximatey as arge as the symbo size. The dashed ines show fits of the expected scaing reations (11) and (12) to the data with a ogarithmic correction term. V. CONCLUSIONS The determination of the order of a phase transition presents a basic issue in the statistica anaysis of physica systems. As there are ony few mode systems that can be tacked anayticay this question is usuay addressed through numerica simuations. A widey used indicator of a first order transition is the metastabiity between two coexisting phases [36, 37, 41]. This eads to a convex dip in the microcanonica entropy and thus to a negative microcanonica specific heat in the transition region. In the canonica ensembe the convex dip corresponds to a doube-peak structure in histograms of suitabe observabes. However, the appearance of a doube-peak structure in canonica histograms or a convex dip in the microcanonica entropy is not sufficient to identify a first order transition as such features have been reported in the iterature for certain mode systems that are known to undergo a continuous transition in the thermodynamic imit. To identify the first order character additiona scaing properties of the convex dip in the microcanonica entropy (or the doube-peak structure in canonica histograms [36, 37]) have to be satisfied. In this work we anayzed the Baxter-Wu mode and the two-dimensiona four-state Potts mode, which are known to have continuous transitions in the thermodynamic imit, but for which doube-peak structures in canonica histograms are observed [41, 42]. Using a recenty deveoped phenomenoogica finite-size scaing theory for systems with continuous phase transitions [26] the possibe properties of the microcanonica entropy of modes where α ν have been investigated. It turns out that finite-size effects may indeed ead to a convex dip in the entropy of any finite system athough the modes show continuous transitions. The finite-size scaing reations for the microcanonica entropy predict in this case certain scaing properties of the convex dip that are different from those encountered at a first order transition. The numerica studies of the Baxter-Wu system and the four-state Potts mode with periodic boundary conditions indeed revea an entropy with a convex dip. The properties of the convex dip are in agreement with the expected scaing reations deduced from the microcanonica finite-size scaing theory. This gives strong support to the vaidity of the scaing theory proposed in [26] for finite microcanonica systems. The systems with open boundaries on the other hand do not show a convex dip in the corresponding entropy. For this case too, the computed specific heat satisfies the scaing reation of the microcanonica scaing theory. The different behaviour for different boundary conditions coud be traced back to the quaitativey different properties of the expansion coefficient A 2 in (7). Future studies have to carify in more detais the infuence of different boundary conditions on the finite-size behaviour of the microcanonica entropy function. Even so the appearance of a convex dip in the microcanonica entropy usuay signas a first order transition,

8 8 it may aso show up in systems with a continuous transition. We have shown in this work that its appearance in those cases (and therefore aso the reated appearance of a doube-peak structure in canonica histograms) can be traced back to a finite-size effect and that the properties of the dip, which are markedy different from those associated with a first order transition, can be understood within a scaing theory deveoped for the description of continuous phase transitions in the microcanonica ensembe. [1] D. H. E. Gross, Microcanonica Thermodynamics: Phase Transitions in Sma Systems (Lecture Notes in Physic 66) (Word Scientific, Singapore 21). [2] M. Peiming and H. Behringer, Phase Transitions 78, 787 (25). [3] W. Thirring, Z. Phys. 235, 339 (197). [4] M. Bixon and J. Jortner, J. Chem.Phys. 91, 1631 (1989). [5] A. Hüer, Z. Phys. B 93, 41 (1994). [6] R. S. Eis, K. Haven and B. Turkington, J. Stat. Phys. 11, 999 (2). [7] T. Dauxois, P. Hodsworth and S. Ruffo, Eur. Phys. J. B 16, 659 (2). [8] D. H. E. Gross and E. V. Votyakov, Eur. Phys. J. B 15, 115 (2). [9] M. Peiming and A. Hüer, J. Stat. Phys. 14, 971 (21). [1] J. Barré, D. Mukame and S. Ruffo, Phys. Rev. Lett. 87, 361 (21). [11] I. Ispoatov and E. G. D. Cohen, Physica A 295, 475 (21). [12] F. Guminei and Ph. Chomaz, Phys. Rev. E 66, 4618 (22). [13] F. Leyvraz and S. Ruffo, J. Phys. A 35, 285 (22). [14] F. Bouchet and J. Barré, J. Stat. Phys. 118, 173 (25). [15] M. Costeniuc, R. S. Eis, H. Touchette, and B. Turkington, J. Stat. Phys. 119, 1283 (25). [16] M. Costeniuc, R. S. Eis, H. Touchette, and B. Turkington, Phys. Rev. E 67, 2615 (26). [17] M. Kastner, Physica A 359, 447 (26). [18] M. Promberger and A. Hüer, Z. Phys. B 97, 341 (1995). [19] M. Kastner, M. Promberger, and A. Hüer, J. Stat. Phys. 99, 1251 (2). [2] A. Hüer and M. Peiming, Int. J. Mod. Phys. C 13, 947 (22). [21] H. Behringer, J. Phys. A 36, 8739 (23). [22] H. Behringer, J. Phys. A 37, 1443 (24). [23] J. Hove, Phys.Rev.E 7, 5677 (24). [24] J. Naudts, Europhys. Lett. 69, 719 (25). [25] M. Peiming, H. Behringer, and A. Hüer, Phys. Lett. A 328, 432 (24). [26] H. Behringer, M. Peiming, and A. Hüer, J.Phys.A 38, 973 (25). [27] A. Richter, M. Peiming, and A. Hüer, Phys. Rev. E 71, 5675 (25). [28] F. R. Brown and A. Yeguap, Phys.Lett.A 155, 252 (1991). [29] M. Schmidt, Z.Phys.B 95, 327 (1994). [3] D. J. Waes and R. S. Berry, Phys.Rev.Lett. 73, 2875 (1994). [31] D. H. E. Gross, A. Ecker, and X. Z. Zhang, Ann. Phys. Lpz (1996). [32] M. Deserno, Phys. Rev. E 56, 524 (1997). [33] S. B. Ota and S. Ota, J. Phys.: Condens.Matter 12, 2233 (2). [34] M. D Agostino, F. Guminei, Ph. Chomaz, M. Bruno, F. Cannata, R. Bougaut, F. Gramegna, I. Iori, N. Le Neindre, G. V. Margagiotti, A. Moroni and G. Vannini, Phys. Lett. B 473, 219 (2). [35] M. Schmidt, R. Kusche, T. Hipper, J. Donges, W. Kronmüer, B. von Issendorff, and H. Haberand, Phys. Rev. Lett. 86, 1191 (21). [36] J. Lee and J. M. Kosteritz, Phys. Rev. Lett. 65, 137 (199). [37] J. Lee and J. M. Kosteritz, Phys. Rev. B 43, 3265 (1991). [38] W. Janke, Nuc. Phys. B (Proc. Supp) 63A-C, 631 (1998). [39] R. J. Baxter and F. Y. Wu, Phys. Rev. Lett. 31, 1294 (1973). [4] R. J. Baxter, Exacty Soved Modes in Statistica Mechanics (London, Academic Press, 1982). [41] M. Fukugita, H. Mino, M. Okawa, and A. Ukawa, J. Phys. A 23, L561 (199). [42] N. Schreiber and J. Ader, J. Phys. A 38, 7253 (25). [43] A. Stemmer and A. Hüer, Phys. Rev. B 58, 887 (1998). [44] M. N. Barber, in Phase transition and Critica Phenomena Vo. 8, edited by C. Domb and J. L. Lebowitz (London, Academic Press, 1983). [45] V. Privman, Finite Size Scaing and Numerica Simuation of Statistica Systems (Singapore, Word Scientific, 199) [46] M. S. S. Chaa, D. P. Landau, and K. Binder, Phys. Rev. B 34, 1841 (1986). [47] H. Behringer J. Stat. Mech.: Theory Exp. 25, 614 (25). [48] C. Borgs and R. Kotecký, J. Stat. Phys. 61, 79 (199). [49] C. Borgs and R. Kotecký, Phys. Lev. Lett. 68, 1734 (1992). [5] C. Borgs, R. Kotecký, and S. Mirace-Soé, J. Stat. Phys. 62, 529 (1991). [51] R. B. Potts, Proc. Camb. Phi. Soc. 48, 16 (1952). [52] F. Y. Wu, Rev. Mod. Phys. 54, 235 (1982). [53] W. Janke and R. Vianova, Nuc. Phys. B 489, 679 (1997). [54] R. Badke, P. Reinicke, and V. Rittenberg, J. Phys. A 18, 73 (1985). [55] R. Badke, P. Reinicke, and V. Rittenberg, J. Phys. A 18, 653 (1985). [56] S. S. Martinos, A. Maakis, and I. Hadjiagapiou, Physica A 352, 447 (25); in this work the density of states of the Baxter-Wu mode has been determined in an approximate way and has then been used in a standard canonica anaysis. [57] K. Binder, Rep. Prog. Phys. 5, 783 (1987). [58] M. Nauenberg and D. J. Scaapino, Phys. Rev. Lett. 44, 837 (198). [59] J. Saas and A.D. Soka, J. Stat. Phys. 88, 567 (1997). [6] M. A. Novotny, D. P. Landau, and R. H. Swendsen, Phys. Rev. B 26, 33 (1982).