arxiv:cond-mat/ v1 [cond-mat.stat-mech] 10 Jul 2001

Transcription

1 Random-luster mult-hstogram samplng for the q-state Potts model arxv:ond-mat/010701v1 [ond-mat.stat-meh] 10 Jul 001 Martn Wegel and Wolfhard Janke Insttut für Theoretshe Physk, Unverstät Lepzg, Augustusplatz 10/11, Lepzg, Germany Chn-Kun Hu Insttute of Physs, Aadema Sna, ankang, Tape, Tawan 1159, Republ of Chna (Dated: July 06, 001) Usng the random-luster representaton of the q-state Potts models we onsder the poolng of data from luster-update Monte Carlo smulatons for dfferent thermal ouplngs K and number of states per spn q. Proper ombnaton of hstograms allows for the evaluaton of thermal averages n a broad range of K and q values, nludng non-nteger values of q. Due to restrtons n the samplng proess proper normalzaton of the ombned hstogram data s nontrval. We dsuss the dfferent possbltes and analyze ther respetve ranges of applablty. PACS numbers: Ln, q, Hk I. ITRODUCTIO Durng the last deade the queston of how to make most effent use of the data sampled durng a Monte Carlo (MC) smulaton has reeved an nreasng amount of attenton. The dea of reweghtng [1] of tme-seres data from a sngle anonal smulaton at a gven fxed value of a ouplng parameter (.e., most ommonly temperature or magnet feld) to nearby regons of the ouplng-parameter spae, allows for the analyss of thermal averages as ontnuous funtons of external parameters and thus a muh more prese determnaton of extremal, pseudo-rtal ponts. As an extenson of ths, the ombnaton of data from smulatons at dfferent ponts n the ouplng-parameter spae, ommonly known as mult-hstogram tehnque [], n prnple allows to get aurate estmates for thermal averages over a marosopal regon of ouplngs from a relatvely small number of smulaton (that, however, generally has to be nreased wth the sze of the system). The bas problem wth ollapsng data from dfferent smulatons s that of fndng the orret relatve normalzaton of the sngle hstograms. Consder the sampled energy hstogram ĤK (E) of, e.g., an Isng model smulaton at the ouplng K = Jβ, onsstng of energymeasurements. The thermal average of an observable A(E) at K s just gven by the tme-seres average n the mportane-samplng sheme and thus nsenstve to the value of the partton funton at that pont. Combnng two hstograms, however, amounts to the ombnaton of the temperature-ndependent expressons Z K Ĥ K (E)/ e KE (1) Eletron address: wegel@tp.un-lepzg.de Eletron address: janke@tp.un-lepzg.de Eletron address: huk@phys.sna.edu.tw for dfferent smulatons, where the partton funton Z K appears as a normalzaton onstant. Thus, for orret relatve normalzaton of the hstograms to be ombned, one has to know the rato of partton funtons Z K1 /Z K or, equvalently, dfferenes n the free energy denstes f K1 f K at the smulated ouplngs K. In Ref. [] ths problem has been solved by an teratve soluton of self-onssteny equatons for the free energy dfferenes at adjaent smulaton ouplngs K. Sne the ombnaton gven n Eq. (1) s nothng but an estmator for the densty of (energy) states Ω(E), mult-hstogramng data analyss amounts to estmatng the densty of states of the varable that s thermodynamally onjugate to the onsdered ouplng parameter. Gong to the random-luster representaton of the Potts model,.e., ts nterpretaton as orrelated perolaton model [3, 4, 5], the relevant densty of states s gven by the number g(b, n) of bond onfguratons wth b bonds and n lusters on the latte. Apart from ganng ontrol over two parameters, the thermal ouplng K and the number of states q, ths language suggests the use of luster estmators for thermal averages lke orrelaton funtons, whh are known to yeld a varane-reduton n ertan stuatons [6]. One of us [7, 8] has proposed a mult-hstogram tehnque for the q-state Potts model and smulatons at dfferent temperatures makng use of the samplng of luster deompostons of the latte as they our n the Swendsen-Wang luster-update algorthm [9]. There, the relatve normalzaton of the ndvdual hstograms at ouplngs K s aomplshed by makng use of the known absolute number of onfguratons wth b atve bonds on the latte, whh s just gven by the bnomal ( E b), E beng the total number of bonds of the latte. Whle ths method appears advantageous at frst sght and gves ne results from the ases of perolaton (q 1) [10] and the Isng model (q = ) [8], we fnd that ths proedure s not the best hoe of normalzaton for smulatons of Potts models wth q

2 larger than or 3 and propose a dfferent approah for normalzaton to rumvent ths problem. The outlne of the paper s as follows. In Seton II we restate the mult-hstogram approah of Ref. [8] n the random-luster representaton ( RC hstogramng ), whh was orgnally formulated for smulatons at fxed q only, and generalze t to smulatons of multple q values. Applyng t to the q = 10 Potts model n two dmensons we fnd large devatons from the expeted results. As an alternatve, n Seton III we propose an adaptvely normalzed RC mult-hstogramng ansatz. We dsuss detals of ts mplementaton and present a omparatve reweghtng analyss for the q = 10 ase. For energyrelated observables we also ompare hstogramng n the random-luster language to the samplng of the energy densty of states n the well-known framework of hstogramng n the energy/magnetzaton language ( EM hstogramng ) []. Comparng both methods, n Seton IV we trak down the observed devatons wth the frst ansatz to be a result of the applaton of the above mentoned normalzaton ondton. Ths problem an thus be resolved by the seond ansatz. Fnally, Seton V ontans our onlusons. II. RC HISTOGRAMS AD ABSOLUTE ORMALIZATIO Consder the Hamltonan of the q-state Potts model n zero magnet feld, H = J,j δ(σ, σ j ), σ = 1,..., q, () on a general graph G wth stes and E bonds. Transformng to the random-luster representaton [3], the partton funton beomes Z e K,j δ(σ,σj) = (e K 1) b(g ) q n(g ), (3) {σ } G G where the sum runs over all bond onfguratons G on the graph (subgraphs), and K = denotes the thermal ouplng. ote that the formulaton (3) n ontrast to that of Eq. () allows for a natural ontnuaton of the model to non-nteger values of the parameter q. Usng the subgraph expanson of the q-state Potts n external feld, one of us [5] has shown that the q-state Potts model an be onsdered as a bond-orrelated perolaton model (BCPM) wth bond oupaton probablty p = 1 e K. Eq. (3) an be rewrtten as: Z p,q (G) = e KE p b(g ) (1 p) E b(g ) q n(g ) G G E = e KE g(b, n)p b (1 p) E b q n, (4) b=0 n=1 where g(b, n) denotes the number of subgraphs of G wth b atvated bonds and n lusters resultng therefrom. Ths purely ombnatoral quantty orresponds to the densty of states of the BCPM. The Swendsen-Wang luster-update algorthm generates bond onfguratons drawn from the equlbrum anonal dstrbuton of ths model. Thus, the probablty for the ourrene of a subgraph wth b bonds and n lusters s gven by P p,q (b, n) = W 1 p,q (G)g(b, n)pb (1 p) E b q n, (5) whh n turn s the expetaton value of the normalzed sampled hstogram of bond onfguratons,.e., P p,q (b, n) = Ĥp,q(b, n)/, where denotes the length of the tme seres of measurements. Here, we separated the ommon fator exp(ke) from the partton funton: Z p,q (G) = e KE W p,q (G). (6) An estmator for the densty of states g(b, n) s therefore gven by Ĥ p,q (b, n) ĝ(b, n) = W p,q (G) p b (1 p) E b q n. (7) Sne the redued partton funton W p,q (G) s a pror unknown, the orret normalzaton of ths estmator s not known at the begnnng. Probably the most obvous way of fxng the normalzaton would be to estmate the redued partton funton W p,q (G) dretly from Eq. (7). One an do better than that, however, by onsderng the aumulated densty g(b), whh s obvously just a bnomal [7], g(b) = ( ) E g(b, n) =. (8) b n Imposng ths restrton on the estmate ĝ(b, n) also, one arrves at ( E Ŵ p,q (G) Ĉ p,q (b) b) p b (1 p) E b = n Ĥp,q(b,, (9) n)q n so that the absolute values of ĝ(b, n) are now fxed by E ndependent normalzaton ondtons, one for eah number of atve bonds b. Thus we have the followng estmate for the densty of states [8]: ĝ(b, n) = Ĉp,q(b)Ĥp,q(b, n)q n. (10) ow, we want to ombne the estmates ĝ () (b, n) from several smulatons at dfferent parameters (p, q ),.e., we want to do mult-hstogramng n both parameters, p and q. Then, we have ĝ(b, n) = α (b, n)ĝ () (b, n), α (b, n) = 1. (11) Sne we want to mnmze the varane ˆσ [ĝ] of the fnal estmate and the dfferent smulatons are statstally ndependent, the orret hoe of the weghts α obvously s gven by α (b, n) = 1/ˆσ [ĝ () (b, n)] 1/ˆσ [ĝ () (b, n)]. (1)

3 3 From Eq. (10) we have ĝ(b, n) = ˆσ [ĝ () (b, n)] = Ĉ p,q (b)q n ˆσ [Ĥp,q (b, n)] Ĉ p,q (b)q n Ĥ p,q (b, n), (13) suh that the varane-optmzed estmate for g(b, n) beomes ( E b) µ Ĥp,q (b, µ)q µ q n j [ ν Ĥp j,q j (b, ν)qj ν ]. qj n [Ĥp j,q j (b, n)] 1 (14) In wrtng ths expresson we allow for several approxmatons: frst, we treat Ĉp,q (b) as a parameter n Eq. (13) nstead of takng ts own varane nto aount; ths s justfed by the lear suppresson of varane of ths quantty as a sum as ompared to the the varane of ts summands Ĥ p,q (b, n). Seondly, we take ˆσ [Ĥ(b, n)] = Ĥ(b, n),.e., we treat the ndvdual bns (b, n) as ndependently dstrbuted aordng to an unorrelated 1/ statsts, whh wll n general not be exatly fulflled. Sne those assumptons only affet the varane of the fnal estmate, however, and do not ntrodue a bas, we onsder them justfed. Fnally, we do not take autoorrelatons between suessve measurements (b, n) nto aount,.e., we assume here and n the followng that measurements n the samplng proess are taken wth a frequeny around 1/τ nt, where τ nt denotes the ntegrated autoorrelaton tme, resultng n an effetvely unorrelated tme seres. From the partton funton Eq. (4) we nfer the followng luster-language estmators for the free energy densty f = F/, the per-ste nternal energy u = U/ and spef heat v = C v /, ˆf = 1 K lnẑp,q(g), û = 1 p b ĝ, ĉ v = K p [ ] (b b ĝ) ĝ (1 p) b ĝ, (15) f. Appendx A. Here, the estmated expetaton value of an observable O(b, n) s defned as O(b, n) ĝ Ẑ 1 p,q(g)e KE E b=0 n=1 ĝ(b, n) p b (1 p) E b q n O(b, n). (16) For the evaluaton of magnet observables one has to dstngush perolatng lusters, denoted by ndes π, from non-perolatng, fnte lusters, denoted by ndes φ. Let {(G )} be the set of lusters of a subgraph G of the latte and n (G ) the number of stes n luster of G. Then, onsder the followng mroanonal averages: m π 1 (b, n) 1 g(b, n) G G,b(G )=b, n(g )=n { π (G )} n π (G ), (17) - F/ (C-C ex )/C ex 10 3 (F-F ex )/F ex RC hstograms EM hstograms exat RC hstograms EM hstograms exat RC hstograms EM hstograms exat FIG. 1: Results from the EM and RC mult-hstogram analyses of tme seres from 9 luster update smulatons of the two-dmensonal q = Potts model on a = 16 latte. (a) Free energy densty (left sale) and spef heat (rght sale) as a funton of the ouplng K = as ompared to the exat soluton of Ref. [11]. (b) Relatve devaton of the results for the spef heat from the exat soluton for both methods. () Relatve devaton for the free energy. All results shown are re-saled from the q = Potts model to the Isng model formulaton to ft the results from Ref. [11]..e., the average number of stes n perolatng lusters (a) (b) () C V / hallo hallo

4 4 6 5 smulaton exat, =16 exat, = (a) EM hstograms RC hstograms σ RC (U)/σ EM (U) 4 3 -U/ σ RC (M)/σ EM (M) FIG. : Rato of standard devatons for estmates of (a) the nternal energy U and (b) the spontaneous magnetzaton M of the q = Potts model on a = 16 latte and RC and EM hstogram analyses as a funton of the ouplng K =. The varanes are estmated by a jakknfe tme seres analyss [1]. The sold lne of (a) shows the exat result of Eq. (4) and Ref. [11], the dashed lne that for from Ref. [13]. In (b) we use the defnton (A9) for K < 0.8 and (A8) for K 0.8. for subgraphs wth b atve bonds and n lusters, m π (b, n) 1 g(b, n) G G,b(G )=b, n(g )=n { π (G )} (b) n π (G ), (18).e., the mean square number of stes n non-perolatng lusters for those subgraphs, and m φ 3 (b, n) 1 g(b, n) G G,b(G )=b, n(g )=n { φ (G )} [n φ (G )], (19).e., the mean squared sum of the sze of nonperolatng lusters. These mroanonal averages obvously an be estmated by addng { π (G )} nπ (G ) FIG. 3: Internal energy of the two-dmensonal q = 10 Potts model on a = 16 square latte wth perod boundary ondtons as gven by the random-luster (RC) and energy representaton (EM) mult-hstogram analyses from smulatons for dfferent thermal ouplngs K. The transton pont of the nfnte system s gven by K t = ln(1 + 10) 1.46 [14]. to ˆMπ 1,p,q (b, n) for eah bond onfguraton G wth {b(g ) = b, n(g ) = n}, where n π (G ) should be taken 0 for non-perolatng bond onfguratons, addng [ { π (G )} nπ (G )] to ˆM,p,q(b, π n), and addng { φ (G )} [nφ (G )] to ˆM φ 3,p,q (b, n) for eah suh observed onfguraton. Then, f we defne Ĥ(b, n) = p,q Ĥ p,q (b, n), ˆM 1 π(b, n) = p ˆMπ,q 1,p,q (b, n), ˆM π(b, n) = p ˆMπ,q,p,q (b, n), ˆM φ 3 (b, n) = p ˆMφ,q 3,p,q (b, n), (0) we have the followng estmates for m π 1 (b, n), mπ (b, n), and m φ 3 (b, n), ˆm π/φ 1//3 ˆMπ/φ 1//3 (b, n) (b, n) =, (1) Ĥ(b, n) whh, fnally, result n the followng expressons for the (zero-feld) magnetzaton m and the magnet suseptblty χ, ˆ m = q 1 ˆm π 1 q (b, n) ĝ + 1 q, ( ) q 1 [ ] ˆ χ = m π ĝ m π 1 ĝ q + q 1 q m φ 3, () ĝ f. Appendx A. ote, that we smply add up hstograms from dfferent smulatons n Eqs. (0) and (1) wthout usng any reweghtng fators n p and q. Ths s orret sne the ondtonal probablty of the ourrene

5 5 of, say, a gven number of stes n perolatng lusters n a subgraph wth b atve bonds and n lusters does no longer depend on p and q. The order parameter m of the Potts model s gven by [14] m = q m 1 q 1, (3) and the orrespondng, resaled suseptblty s χ = [(q 1)/q] χ. As a frst omparatve test for the method we performed a Swendsen-Wang luster MC smulaton for the Isng model ase (q = ) on a small, smple-ub = 16 latte wth perod boundary ondtons. We gathered hstograms from 9 dfferent smulatons at the ouplngs K = 0.1, 0.,..., 0.8 and K = K = 1 ln(1 + ), where the ouplngs are gven n the language of the Isng model n ths ase,.e., are half of the ouplngs of the orrespondng q = Potts model. Eah run sampled 17 = bond onfguratons resultng n orrespondng tme seres of (b, n) and of the energy/magnetzaton pars (E, M) for omparson wth the EM hstogramng method. Thus, any dfferenes n the results must be solely due to the method of data analyss, the underlyng smulaton data beng exatly dental. For the EM hstograms throughout ths paper we use a mult-hstogram analyss aordng to Ref. [] very smlar to that presented for the RC hstograms n Se. III. Fgure 1 shows the results n omparson to the exat expresson for F and C v on square lattes as gven by Kaufman [15] and analysed by Ferdnand and Fsher [11]. Statstal errors for both analyss shemes were evaluated usng the jakknfe error estmaton tehnque [1]. The relatve devatons (Ĉv C v )/C v from the exat result are noteably larger for the RC hstogram analyss, however n agreement wth statstal errors n both ases, p. Fg. 1(b). The same holds true for the nternal energy U gven by the dfferent estmates, whh s not shown n Fg. 1. ote that for the energy related observables from Eq. (15) only the lmtng dstrbuton P p,q (b) s needed, whh, n general, has a dfferent wdth than the dstrbuton of energes P p,q (E), thus leadng to dfferent varanes. In fat, by nspeton of the random-luster expresson for the spef heat Eq. (15) and omparson wth ts defnton n the energy language as C v = K ( E E ) we an nfer the followng relaton between the varanes of energy estmates n the RC and EM shemes: σrc (U) 1 p = 1 K (U) p σ EM U C v 1. (4) Thus, energy estmates from RC hstograms are always less prese than those from EM hstograms, regardless of the temperature. Fgure (a) shows the rato of jakknfeestmated varanes of the two dfferent estmates of nternal energy, ompared to the result from Eq. (4) wth the exat expressons for U and C v for the q = ase nserted [11]. ote, that from Fg. (a) ths quantty seems to have extremely small fnte-sze orretons. As a remnder, ths shows learly that luster estmators are not always mproved estmators [16], but sometmes deterorated estmators. ote, however, that ths effet wll derease wth nreasng number of states q, at least n the transton regon, sne the sngularty n C v sharpens n ths lmt, whereas the energes U always stay n the range 0 U/. For the q = 10 model s has been observed that at the transton pont P p,q (b) s almost ndstngushable from P p,q (E), when sutably resaled [17]. The mnmum of the exat urve of Fg. at the rtal pont s somewhat n ontrast to the usual noton that luster estmators work best off rtalty [6]; ths result, however, apples to the spn-spn orrelaton funton at medum and long dstanes and to magnet observables lke the suseptblty, whh s the ntegral of the orrelaton funton, whereas the nternal energy U onsttutes the extreme short dstane lmt of ths quantty. For the magnet observables m and χ the stuaton s reversed, the varane of the RC estmators beng strongly redued as ompared to the EM estmators, p. Fg. (b). ote that n ontrast to the EM ase, the RC estmators provde a sngle onsstent defnton of m and χ for both, the broken and unbroken phases, f. Appendx A. For the free energy, on the other hand, devatons for the RC method are by far smaller than those of the EM method, p. Fg. 1(). Moreover, devatons are not overed by statstal errors n the latter ase, a fat we wll omment on later n the next Seton. In the EM ase F s beng fxed by makng ontat wth the non-nteratng lmt K = 0 resp. p = 0, where Z p=0,q (G) = E Ω(E) = q, (5) so that KF(p = 0)/ = lnq. Ths equaton orresponds to the normalzaton ondton Eq. (8). It s obvous that havng a normalzaton ondton for eah number b of atve bonds and therefore, mpltly, for eah (mroanonal) temperature n the RC ase allows for aurate estmaton of the free energy even far away from K = 0, whereas for EM hstograms the results deterorate wth the dstane from the only normalzaton pont K = 0. Thus, for samplng free energes the RC multhstogram tehnque normalzed by Eq. (8) seems to be a good hoe. As a slghtly less trval example we performed smulatons for the q = 10 Potts model on the same latte, whh exhbts a strongly frst-order phase transton. It s well known that luster algorthms are not effent to redue the super-rtal (exponentally strong) slowng down of the loal MC dynams at frst-order transtons. For the small latte under onsderaton, however, autoorrelaton tmes are stll qute moderate, so that one gets relable results wthout havng to resort to more sophstated methods lke mult-anonal smulatons [17, 18]. We gathered data from 11 sngle-hstogram smulatons at ouplngs K = 0.8, 0.9,..., 1.8 wth 0 = measurements eah. Fgure 3 shows the qute astonsh-

6 EM, smulatons at q=10 RC, smulatons at q= U/ ln Ω(E) =0.4 =0.8 =1. exat FIG. 4: Internal energy of the q = 10 Potts model on a = 16 square latte, reweghted from the q = model smulatons shown n Fg. 1 usng multple RC hstograms aordng to Eq. (14). The q = 10 results from the EM multhstogram analyss are shown for omparson E FIG. 5: Densty of states for the q = Potts model n two dmensons on a = 16 latte from sngle-hstogram lusterupdate smulatons at ouplng K = 0.4, 0.8 and 1. usng the estmator Eq. (6). The sold lne shows the exat result of Ref. [0]. ng results for the nternal energy from ths smulaton data usng the analyses n the RC and EM languages, respetvely. aturally, we do not have exat results to ompare wth n ths ase; nevertheless, the results from the EM analyss are ompletely n agreement wth our expetatons and also well ompatble wth results from prevous smulatons [19]. So, obvously, the results from the RC hstogram analyss are strkngly wrong and n a way that s learly not overed by the present statstal errors. Obvously, the results for the spef heat, whh are not shown, look even worse, wth a pronouned, unphysal double-peak resultng from the devatons n nternal energy shown n Fg. 3. Sne as one of ts major strengths n the RC approah we have the possblty of reweghtng n the parameter q also, as a frst lue to the reason for ths onspuous falure we show the outome of usng the q = smulaton data from above for determnng the nternal energy of the q = 10 ase, p. Fg. 4. The agreement wth the dret EM analyss of the q = 10 smulatons s remarkably good onsderng the large dstane n q between the smulaton and analyss ponts. Comparng Fgs. 3 and 4 t s qute natural to suspet that the applaton of the normalzaton ondton (8) s not a proper hoe for smulaton data from larger q models. III. RC HISTOGRAMS AD ADAPTIVE ORMALIZATIO To understand ths normalzaton problem let us shortly go bak to the samplng of the energy densty of states Ω(E) for the ase of the two-dmensonal q = (Isng) model. Here, exat results are not only avalable for thermal averages, but for Ω(E) tself [0]. Usng the K = 0 normalzaton ondton (5), a sngle-hstogram estmator for the densty of states n the energy language would be gven by ˆΩ(E) = E Ĥ K (E)e KE ĤK(E)e KE. (6) Ths works qute well n the hgh-temperature phase and for small lattes. For lower temperatures, however the hstogram loses ontat wth the normalzaton pont K = 0, resultng n large devatons from the orret normalzaton, p. Fg. 5. Clearly, eah smulaton samples only a rather small wndow of energy spae; from the exponental n the denomnator of Eq. (6), however, onfguratons near the maxmal energy E = reeve the largest weght n the sum, so that mssng those onfguratons, whh s the ase for large K, results n an exponentally wrong normalzaton fator (lnear n ln Ω(E)). In other words, the absolute normalzaton ondton (5) reweghts the hstogram data to the pont K = 0, whh wll have no relable outome f the overlap between the hstograms at the smulaton ouplng and at K = 0 s too small or even vanshes. ote also, that the statstal error bars gven n Fg. 5 do not reflet ths fundamental falure, although t s statstal n nature. Ths s due to the fat that the usual mplementaton of error estmaton shemes for hstograms take the error of hstogram bns wthout entres to be zero, whereas aordng to 1/ statsts t should n some sense be onsdered nfntely large. Expetng a smlar samplng-related normalzaton falure for the RC hstograms normalzed by the ondton (8) let us relax ths absolute normalzaton and apply an adaptve normalzaton sheme as n the orgnal EM mult-hstogramng formulaton of Ref. []. Consder sngle-hstogram estmates of the partton funton from

7 7 - F/ RC hstograms EM hstograms U/ s obvous that for the teraton to onverge, one needs some overlap between the Ĥ(b, n) hstograms between adjaent smulatons,.e., at least parwse overlap. Apart from ths restrton, however, we fnd ths teratve sheme to be very well-behaved, onvergng rapdly n every ase that fulfls the overlap-ondton. To get started, we use frst-guess values of the F from thermodynam ntegraton. Assume that the smulaton ponts () = (P, q ) are ordered suh that the hstograms of () and ( + 1) have reasonable overlap; then FIG. 6: Internal energy and free energy of the q = 10 Potts model on a = 16 square latte as omputed by the adaptvely normalzed RC mult-hstogramng sheme aordng to Eqs. (9) and (30). The results from the EM mult-hstogram analyss of the same data are shown for omparson. smulatons at (p, q ) and defne redued free energes F as F = 1 K lnz p,q (G) = E F K. (7) From Eq. (7) estmates of the densty g(b, n) from the sngle hstograms Ĥp,q (b, n) Ĥ() (b, n) are gven by: ĝ () Ĥ () (b, n) (b, n) = e F p b (1 p ) E b q n. (8) One agan ombnng these estmates n a varaneoptmzed way as above n Seton II, treatng the F as parameters wth zero varane, we arrve at the followng expresson ĝ(b, n) = e F p b (1 p ) E b q n j j e Fj p b j (1 p j) (E b) q n j [Ĥ() (b, n)] 1. (9) ow, from ths estmate one has the followng a posteror relaton for omputaton of the parameters F : e F = E b=0 n=1 0.5 ĝ(b, n)p b (1 p ) E b q n. (30) Eqs. (9) and (30) form a par of equatons to be solved self-onsstently for the determnaton of the parameters F, whh an be straghtforwardly terated by pluggng n the results for F from Eq. (9) nto Eq. (30) and ve versa. One an mprove on that by applyng more sophstated teraton shemes lke, e.g., the ewton- Raphson teraton [1]. We fnd, however, that the radus of onvergene of ths method s qute small; therefore, we adaptvely revert to the smple teraton f the proedure leaves the ewton-raphson onvergene regon. It E Ĥ ( 1) (b, n) p b e F = (1 p ) E b q n b=0 n=1 1 p b 1 (1 p 1) E b q 1 n e F 1 (31) s a good startng pont for the desrbed teraton sheme. F 1 an be hosen arbtrarly, sne the gven par of equatons s obvously nvarant under a global shft F F F 1. Thus, we have determned the fnal estmate ĝ(b, n) only up to a global fator. To fx ths last normalzaton we propose two dfferent possbltes; on the one hand, we an use the free model lmt,.e., evaluate F p=0,q from Eq. (30) and use Eq. (5) for any q: F p=0,q = lnz p=0,q (G) KE = lnq. (3) On the other hand, also the q = 1 partton funton s trval, [ E ( ] E F p,q=1 = ln )p b (1 p) E b = 0, (33) b b=0 and an serve as normalzaton pont for arbtrary p. In prate the best hoe depends on the set of smulated ouplngs (p, q ): for large-q smulatons one mght want to resort to Eq. (3), whle otherwse Eq. (33) should be the better hoe. ow, we an re-onsder the nternal energy of the q = 10 ase from above wth the new, adaptvely normalzed RC mult-hstogramng sheme. Fgure 6 shows nternal energy and free energy from ths analyss as ompared to the EM mult-hstogram approah. As far as the error estmates are onerned, we apply the jakknfe proess to the whole teraton run,.e., the teraton sheme for fxng the weghts F s done for eah jakknfe blok of data separately, takng full aount of statstal errors. Clearly, now the results from both approahes perfetly agree, the devatons of Fg. 4 have vanshed. As antpated n Se. II, also the luster estmator for the nternal energy performs noteably better than n the q = ase, suh that at least n the rtal regon t s qute omparable n preson to the EM estmator, p. Fg. 7. For the free energy t s obvous that wth the adaptve normalzaton sheme of RC hstograms we lose the espeally hgh preson throughout the whole K regon obtaned by the applaton of the sum rule (8) n Fg. 1 for the Isng model. To amend ths, havng fxed the relatve normalzaton of the sngle hstograms adaptvely,

8 I II III 3 00 σ RC (U)/σ EM (U) n FIG. 7: Rato of standard devatons for estmates of the nternal energy U of the two-dmensonal q = 10 Potts model on a = 16 latte and RC and EM mult-hstogram analyses as a funton of the ouplng K =. The varanes are estmated by a jakknfe tme-seres analyss [1]. The RC hstograms are normalzed aordng to Eqs. (9) and (30). one mght onsder applyng the sum rule (8) to the fnal result ĝ(b, n) nstead of usng the normalzatons Eq. (33) or Eq. (3). Ths, however, gves results lookng almost dental to those shown above n Fg. 3,.e., the large devatons reappear, whh learly reveals the soure they are resultng from b FIG. 8: Densty plot of the relatve dfferenes [ĝ (abs) ĝ (rel) ]/ĝ (rel) of the densty of states as sampled from the q = 10 Potts model on a = 16 latte by the absolutely normalzed RC hstogramng sheme of Eq. (14) (ĝ (abs) ) and the adaptvely normalzed sheme of Eq. (9) (ĝ (rel) ), respetvely. Dark shadng ndates that ĝ (abs) (b, n) > ĝ (rel) (b, n) and ve versa. bonds, or, for the square latte, n 1 (E b), (36) E IV. COMPARISO OF THE METHODS n b + 1. (37) The effet of ths normalzaton problem should also be learly seen n the fnal estmates for the densty of states g(b, n) from the two RC hstogramng methods. In Fg. 8 we show a densty plot of the relatve dfferenes of the estmated densty of states ĝ(b, n) from the absolutely normalzed hstogramng sheme of Eq. (14), ĝ (abs) (b, n), and of the adaptvely normalzed sheme of Eqs. (9) and (30), ĝ (rel) (b, n),.e., the quantty ˆ (b, n) ĝ(abs) (b, n) ĝ (rel) (b, n). (34) ĝ (rel) (b, n) ote that the range of possble value pars (b, n) s restrted by two smple bounds n the (b, n) plane. Frst, startng from the pont (b = 0, n = ) eah added bond an at most redue the number of lusters by one, namely by onnetng two prevously unonneted lusters,.e., one has n b. (35) On the other hand, startng from the opposte pont (b = E, n = 1) one has b/e bonds per ste, so that for produng a new luster one must at least remove /E Apart from sngle ponts near those bounds, all onfguratons wthn ths trangle an atually appear n a Potts model smulaton wth non-vanshng probablty. ow, from Fg. 8 t s obvous, gven that the estmate ĝ (rel) (b, n) s orret up to an overall fator, that the absolutely normalzed hstogramng estmate ĝ (abs) (b, n) gves too large estmates for b values near the entre b = as ompared to the other regons of b (dark shadng n Fg. 8). Then, onsderng agan the devaton n nternal energy shown n Fg. 3, ts orgn beomes lear: the hstogram Ĥp,q(b, n) for a smulaton somewhat below the transton pont wll be entred around the lne b = b I n Fg. 8; then, usng the densty of states estmate ĝ (abs) (b, n) for evaluatng U, the parts of the hstogram lyng to the rght of b = b I wll have too large weght as ompared to the values b < b I, thus by Eq. (15) resultng n a too large estmate for the nternal energy U. On the other hand, for ouplngs above the transton pont the hstogram wll be entred around b = b III wth too large weghts for b < b III, leadng to estmates for U that are too low. Dretly n the vnty of the transton pont, devatons n normalzaton are symmetr wth respet to the hstogram, whh wll be entred around b = b II, thus leadng to an unbased estmate for U. Ths s exatly the behavor found n Fg. 3. Fnally, ontemplatng on the reason for the devatons n normalzaton

9 9 σ RC (F)/σ RC,norm (F) FIG. 9: Rato of standard devatons of estmates of the free energy from the adaptvely normalzed RC hstogramng sheme of Eq. (9) wth (σ RC,norm) and wthout (σ RC) a fnal applaton of the sum rule Eq. (8) to the densty of states after determnng the weghts. The used tme seres data nludes both, the q = and the q = 10 smulatons reported above. shown n Fg. 8 n the frst plae, t beomes obvous that they have the same orgn as those shown n Fg. 5. The exponental fator q n from the sum rule Eq. (9) attahes large weght to the onfguratons wth small number of lusters n; f, however, hstograms mss entres for small n, as s the ase for hstograms n the transton regon b of Fg. 8, the sum n Ĥp,q(b, n)q n wll beome too small, resultng n too large normalzaton fators Ĉp,q(b). Thus, for the applaton of the sum rule (8) to the fnal result from the adaptvely normalzed RC hstogramng sheme to work relably, one always has to nlude hstograms from small-q smulatons, lke perolaton (q 1) or the Isng model (q = ), that produe onfguratons wth relatvely small numbers of lusters n. To llustrate ths, we ombned the data from the q = and q = 10 smulatons reported above, used the hstogramng sheme Eqs. (9) and (30) to get results for q = 10 and appled the sum rule Eq. (8) afterwards,.e., the normalzaton of the sngle hstograms was found adaptvely, whereas the total hstogram was normalzed by the sum rule (8); ths yelds results for the nternal and free energes ndstngushable from those of the pure q = 10 results of Fg. 6. For the free energy, however, the sze of statstal errors s largely affeted by the fnal normalzaton, f. Fg. 9. For most of the ouplngs shown, the estmate from the fnally sum-rule-normalzed densty of states s up to about 10 tmes more aurate n terms of the statstal errors. The presene and sze of suh a gan for a gven ouplng s not manly physally motvated, but rather depends on the relaton of the smulaton ponts (q, p ) to the ponts of data analyss. V. COCLUSIOS We have onsdered mult-hstogram data analyses of tme seres from luster-update Monte Carlo smulatons of the q-state Potts model n the random-luster language. Generalzng the orgnal formulaton of Ref. [8] to the ase of smulatons of dfferent number of states q, we found the orgnal ansatz of absolutely normalzng the ndvdual hstograms wth a geometral sum rule for fnte-length tme seres to produe large devatons from the expeted behavor when appled to ases q larger then about 3 or 4 n two dmensons. We trak ths error down to a msmath between exponental suppresson of a part of the state spae (b, n) and a smultaneous exponental enhanement of ths regon n the sum rule Eq. (9). To rumvent ths problem, we propose a dfferent ansatz, normalzng the hstograms adaptvely va a set of self-onssteny equatons amng at the mnmzaton of the varane of the fnal estmate of the densty of states g(b, n). Absolute normalzaton over the whole temperature regon an stll be mantaned by makng ontat wth the trval partton funton of the perolaton lmt q 1 or by ombnng large- and small-q data and applyng the sum rule (8) after the adaptve normalzaton. Ths new approah does not exhbt the lmtatons of the absolutely normalzed ansatz to small-q smulatons. Comparng the newly ntrodued, adaptvely normalzed random-luster RC mult-hstogram tehnque wth mult-hstogramng n the energy/magnetzaton EM language, we an make the followng statements: (a) The luster varables (b, n) form the natural state spae for the analyss of the Potts model. Usng the Swendsen-Wang luster-update algorthm, these numbers are automatally known as a by-produt of the update-steps; no addtonal measurement steps are needed. (b) The RC representaton allows for reweghtng n both parameters, the thermal ouplng p resp. K, and the number of states q, wthout systematal errors as n the partal transformaton of Ref. []. Espeally, the model an be onsdered for the ase of non-nteger q. It s easy to ombne data from dfferent-q smulatons to enhane the auray for large q. () Cluster estmators our naturally n the RC language. Although we found that short-dstane observables lke the nternal energy and spef heat are sampled systematally less aurate by luster estmators, ths stuaton s reversed for observables senstve to long-range order lke the magnetzaton, suseptblty and orrelaton funtons. Also, even short-range luster estmators perform omparable to EM language estmators for larger q values, at least n the transton regon. In the RC language, the magnet observables an be defned onsstently throughout the broken and unbroken phases, f. Appendx A. Apart from that, the ombnaton of data from small and large q models an serve as a new method to ope wth the super-rtal slowng down at the frst-order transtons for large q: for suffently large lattes sm-

10 10 ulaton runs wll entrely stay n one of the pure phases dependng on the ntal onfguraton and boundary ondtons due to ergodty breakng at the transton pont. However, ombnng suh data wth smaller-q smulatons from the seond-order or weak frst-order regme allows the adaptve normalzaton sheme to stll fnd the orret normalzaton of the pure-phase hstograms wthout real tunnelng events. Ths approah s smlar n sprt to the smulated temperng tehnque [3, 4]. One mght also thnk of applyng the mult-anonal [18] resp. mult-bond [17] smulaton approah or one of the related tehnques to the samplng of the densty of states g(b, n). Espeally, applaton of the absolute normalzaton Eq. (8) to ths ase mght be of nterest. Ths approah s urrently under nvestgaton. However, samplng the omplete range of possble values n the (b, n) plane wth suffent auray s found to be a omputatonally very demandng problem. In ontrast, n the urrent approah, we stll stk to the physally sensble approah of mportane-samplng,.e., samplng the phase-spae aordng to the loal anonal weghts. Furthermore, we are able to take full advantage of the omputatonal gan of luster algorthms, whereas the generalzed-ensemble algorthms put forward so far employ loal updates (apart from the mult-bond algorthm of Ref. [17]). As an nterestng applaton of our ansatz we suggest the analyss of the tr-rtal pont q, where the order of the thermal transtons hanges from seond to frst order, n three dmensons. There has been qute some debate about the loaton of ths pont, estmates rangng from q =.15 [5] to q =.6 [6]. Furthermore, the unversalty lass, rtal exponents et. of ths transton have not yet been properly analyzed. A test for the q = 4 ase n two dmensons shows that our method s well suted for suh an analyss. Ths problem wll be onsdered n a forthomng publaton. ACKOWLEDGMETS The ollaboraton of the authors has been supported by DAAD-SC grant o. D/ MW aknowledges support by the DFG through the Graduertenkolleg QFT. APPEDIX A: CLUSTER ESTIMATORS Consder the Potts model oupled to an external magnet feld H wth Hamltonan [7] H = J,j δ(σ, σ j ) H δ(σ, 1). (A1) Then, the random-luster representaton of the partton funton on a graph G onsstng of stes and E edges (bonds) s gven by Z p,q (G, B) = G G = e KE G G (e K 1) b(g ) p b(g ) (1 p) E b(g ) [(q 1) + e Bn ], [(q 1) + e Bn ] (A) where the produt runs over the set of lusters {} of the subgraph G, n s the number of stes n luster, K = s the thermal ouplng parameter and B = βh denotes the redued magnet feld. p = 1 e K s the probablty for the atvaton of bonds. The zero-feld, per-ste nternal energy u s then gven by u = [ ] lnzp,q (G, B = 0) K = 1 K [(ek 1) b ] (e K 1) b = 1 p b, (A3) whh shows the lose onneton between the b and E dstrbutons. The zero-feld spef heat v follows from v = K K [ lnzp,q (G, B = 0) ] = K [ 1 p b + K [(e K 1) b ] (e K 1) b = K [ b p b (1 p) b ]. (A4) In the thermodynam lmt, the per-ste, zero-feld magnetzaton m = δ(σ, 1)/ s gven by [ ] m = lm lm lnzp,q (G, B) B 0+ B ]

11 = lm B 0+ = lm B 0+ lm lm e KE p b(g ) (1 p) E b(g ) [(q 1) + e Bn ] n e Bn Z p,q (G, B) (q 1) + e Bn G G [ ] nπ 1 + nφ 1. (A5) (q 1)e Bnπ + 1 (q 1)e Bnφ + 1 π φ 11 Here, we splt the luster ontrbutons of the subgraph nto perolatng lusters π and non-perolatng, fnte lusters φ. In the ndated order of takng the lmts, frst and then B 0, the fators exp( Bn ) take the values 0 and 1 for perolatng and non-perolatng lusters, respetvely. Thus, we arrve at m = nπ + 1 q π φ = q 1 q π nφ nπ + 1 q, (A6) whh expltly reflets the symmetry-breakng nature of the perolatng onfguratons. For the order parameter m, whh vares between 0 for the ompletely dsordered state and 1 for the ground states, we fnd m q m 1 q 1 = π nπ. (A7) For fnte lattes one an retan ths defnton sne the noton of perolatng and non-perolatng lusters s stll well-defned. ote, that ths gves a onsstent defnton of the order parameter throughout the dsordered and broken phases. In ontrast, n the EM language one has to expltly break symmetry n the low-temperature phase, whh s usually done by defnng m K>Kt = max 1 j q δ(σ, j), (A8) whereas for the unbroken phase one uses m K Kt = δ(σ, 1). (A9) Obvously, for fnte lattes, the expetaton values of the RC and EM defntons wll not onde exatly; rtal exponents, however, wll of ourse agree. The zero-feld suseptblty χ s gven by χ = lm B 0+ lm B = m + lm B 0+ ( ) lnz(b) lm ( ) ( q 1 = q π = lm B 0+ n ) nπ π lm 1 (q 1)e Bn + 1 [ ( ) 1 Z(B) Z(B) + B ( n + + nπ + q 1 q 1 Z(, B) 1 (q 1)e Bn + 1 ] Z(B) B ) ] n 1 [(q 1)e Bn + 1] nφ. (A10) φ From Eq. (A10) one reognzes the wdely used mproved luster estmator for the hgh-temperature phase, namely the last term. ote, however, that the orgnal mproved estmator nludes all lusters here nstead of only the non-perolatng ones, whh makes a dfferene for fnte lattes. For fnte lattes, one agan, from Eq. (A10) we have a sngle defnton for both, the unbroken and broken phases. Alternatvely defnng the suseptblty orrespondng to the order parameter m we get ( ) q χ = χ = 1 nφ q 1 q 1 φ ( ) + nπ nπ.(a11) π π

12 1 REFERECES [1] A. M. Ferrenberg and R. H. Swendsen, Phys. Rev. Lett. 61, 635 (1988); Phys. Rev. Lett. 63, 1658(E) (1989). [] A. M. Ferrenberg and R. H. Swendsen, Phys. Rev. Lett. 63, 1195 (1989). [3] C. M. Fortun and P. W. Kasteleyn, Physa 57, 536 (197). [4] A. Conglo and W. Klen, J. Phys. A 13, 775 (1980). [5] C.-K. Hu, Phys. Rev. B 9, 5103 and 5109 (1984); J. Phys. A 19, 3067 (1986). [6] U. Wolff, ul. Phys. B 300, 501 (1988). [7] C.-K. Hu, Phys. Rev. Lett. 69, 739 (199). [8] J.-A. Chen and C.-K. Hu, Phys. Rev. B 50, 660 (1994). [9] R. H. Swendsen and J.-S. Wang, Phys. Rev. Lett. 58, 86 (1987). [10] C.-K. Hu, Phys. Rev. B 46, 659 (199). [11] A. E. Ferdnand and M. E. Fsher, Phys. Rev. 185, 83 (1969). [1] B. Efron, The Jakknfe, the Bootstrap and Other Resamplng Plans (Soety for Industral and Appled Mathemats [SIAM], Phladelpha, 198). [13] L. Onsager, Phys. Rev. 65, 117 (1944). [14] F. Y. Wu, Rev. Mod. Phys. 54, 35 (198). [15] B. Kaufman, Phys. Rev. 76, 13 (1949). [16] U. Wolff, ul. Phys. B 334, 581 (1990). [17] W. Janke and S. Kappler, Phys. Rev. Lett. 74, 1 (1995). [18] B. A. Berg and T. euhaus, Phys. Rev. Lett. 68, 9 (199). [19] W. Janke and S. Kappler, J. Phys. I (Frane) 7, 663 (1997). [0] P. D. Beale, Phys. Rev. Lett. 76, 78 (1996). [1] W. H. Press, S. A. Teukolsky, W. T. Vetterlng, and B. P. Flannery, umeral Repes n C The Art of Sentf Computng (Cambrdge Unversty Press, Cambrdge, 199), nd ed. [] J. Lee and J. M. Kosterltz, Phys. Rev. B 43, 168 (1991). [3] E. Marnar and G. Pars, Europhys. Lett. 19, 451 (199). [4] W. Kerler and A. Weber, Phys. Rev. B 47, (1993). [5] S. Grollau, M. L. Rosnberg, and G. Tarjus, ondmat/ preprnt (000). [6] B. enhus, E. K. Redel, and M. Shk, Phys. Rev. B 3, 6055 (1981). [7] ote that n Ref. [5] a dfferent ouplng between magnet feld and spns s employed.