Sie sind hier: ICP » R. Hilfer » Publikationen

III Finite ensemble scaling

[67.2.1.1] The quantity of main interest for finite ensemble scaling [21, 22, 23] is the macroscopic ensemble sum XMNφ given by (2.19). [67.2.1.2] The idea is to neglect completely its microscopic definition (2.18) in terms of cell variables, and to consider the mesoscopic block variables XMNφy as a starting point. [67.2.1.3] The univariate probability distribution of the ensemble variable is defined as

PXMN(x)=Prob{XMN(φ)x}.(3.1)

[67.2.1.4] Because the ensemble limit automatically generates independent and identically distributed block variables XMNφyj the standard theory of stable laws [27, 28] can be applied. [67.2.1.5] It yields the existence and uniqueness of limiting distributions for the linearly renormalized ensemble sums

ZMNφ=XMNφ-CNDN=j=1NXMNφyj-CNDN(3.2)

where DN>0 and CN are real numbers. [67.2.1.6] Remember that this holds for sums of arbitrary block variables independent of their microscopic definition. [67.2.1.7] The index M serves only as a reminder for the fact that the ensemble limit is used.

[page 68, §0]    [68.1.0.1] The distribution function PZMNx for ZMNφ is given in terms of PXMNx as PXMNDNx+CN and it is thus sufficient to consider PXMNx. [68.1.0.2] The (weak) ensemble limit of these probability distribution functions

limM,NN/M=cPXMNDNx+CN=Hx;ϖXc,ζXc,Cc,Dc(3.3)

exists if and only if Hx;ϖXc,ζXc,Cc,Dc is a stable distribution function whose characteristic function

hk=eikX=-eikxdHx(3.4)

has the form

hk;ϖX,ζX,C,D=expiCk-DkϖXeiπ21-1-ϖXζXsgnk(3.5)

for ϖX1 and

hk;1,ζX,C,D=expiCk-Dk1-iζX2πsgnklogk(3.6)

for ϖX=1. [68.1.0.3] The c-dependence of the parameters ϖXc,ζXc,Cc,Dc has been suppressed to shorten the notation. [68.1.0.4] The parameters ϖX,ζX,C,D obey

    0<ϖX2-1ζX1-<C<0D.(3.7)

[68.1.0.5] If the limit exists, and D0, the constants DN must have the form

DN=NΛN1/ϖX(3.8)

where ΛN is a slowly varying function [28], i.e.

limxΛbxΛx=1(3.9)

for all b>0.

[68.1.0.6] The forms (3.5) and (3.6) of the limiting characteristic functions imply the following scaling relations for the stable probability densities hx;ϖX,ζX,C,D. [68.1.0.7] If ϖX1 then

hx;ϖX,ζX,C,D=D-1/ϖXhx-CD-1/ϖX;ϖX,ζX,0,1(3.10)

holds, while for ϖX=1 one has

hx;ϖX,ζX,C,D=D-1hx-CD-1-2ζXπlogD;ϖX,ζX,0,1.(3.11)

[68.2.0.1] The parameters C and D correspond to the centering and the width of the distribution.

[68.2.1.1] Strictly stable probability densities (i.e. those with ϖX1) are conveniently written in terms of Mellin transforms [29, 30]. [68.2.1.2] This representation is useful for computations and involves the general class of H-functions [31, 32]. [68.2.1.3] For 1<ϖX<2 corresponding to equilibrium phase transitions two cases are distinguished. [68.2.1.4] If ζX1 then [30, 22]

h(x;ϖX,ζX,0,1)=1ϖXH2211(x|1-1/ϖX,1/ϖX1-ϱ,ϱ0,11-ϱ,ϱ)(3.12)

where ϱ=12-ζXϖX+ζX2 and the definition of the general H-function HPQmn is given in the appendix. [68.2.1.5] If ζX=1 then for 1<ϖX<2

h(x;ϖX,±1,0,1)=1ϖXH1110(x|1-1/ϖX,1/ϖX0,1)(3.13)

[68.2.1.6] Similar expressions hold for 0<ϖX<1[30, 22]. [68.2.1.7] The special case ϖX=2 of the general limit theorem (3.3) is the central limit theorem [28] and in this case the stable probability density

hx;2,ζX,C,D=14Dπe-x-C2/4D(3.14)

is the Gaussian distribution with mean C and variance σ2=4D. [68.2.1.8] Note that the right hand side is independent of ζX in this case. [68.2.1.9] Another special case expressible in terms of elementary functions is ϖX=1,ζX=0 where

hx;1,0,C,D=1πDD2D2+x-C2(3.15)

is the Cauchy distribution centered at C and having width D.

[68.2.2.1] For sufficiently large but finite N=L/ξd equation (3.3) implies that the distribution function of ensemble variables may be approximately written as

PXMNx=Rx,M,N,cHx-CNDN;ϖX,ζX,C,D    :    for  x01-Rx,M,N,c1-Hx-CNDN;ϖX,ζX,C,D    :    for  x>0(3.16)

[page 69, §0]    involving a nonuniversal cutoff functionRx,M,N,c such that R0,M,N,c=1 andlimx±Rx,M,N,c=0 for all M,N<. In the ensemble limit the cutoff function must obey

limM,NN/M=cRx;M,N,c=1,(3.17)

for all x and c as a result of equation (3.3). [69.1.0.1] Note that equation (3.17) does not hold for the finite size scaling limit. [69.1.0.2] Instead Table I implies that for the finite size scaling limit

limL,ξL/ξ=cRx;M,N,N/M=Rx;,cd,0(3.18)

if the limit exists, and where now c=L/ξ. [69.1.0.3] The function Rx;,L/ξd,0 may in general differ from unity, and thus the finite size scaling limit may involve a nonuniversal cutoff function which is absent in the finite ensemble limit.

[69.1.1.1] Wherever possible equation (3.16) will be abbreviated as

PXMNxHx-CNDN;ϖX,ζX,C,D.(3.19)

to shorten the equations. [69.1.1.2] If the centering constants are now chosen as

CN=-DNC    :    for  ϖX1-DNC+2πζXDlogD    :    for  ϖX=1(3.20)

then using equations (3.10),(3.11) and (3.8) the basic finite ensemble scaling result [21, 22]

pXMNxhx;ϖX,ζX,0,DNΛN(3.21)

is obtained for the probability density function pXMNx of suitably centered and renormalized ensemble sums. [69.1.1.3] The approximate result (3.21) has formed the basis for the statistical mechanical classification of phase transitions [21, 22].

[69.1.2.1] From the basic result (3.21) the scaling form for the probability density of ensemble averaged block variables X¯MNφ=XMNφ/MN is readily obtained using eq. (3.10) as

p¯X¯MNxL/ξd1-1/ϖXDΛL/ξd1/ϖXhxL/ξd1-1/ϖXDΛL/ξd1/ϖX;ϖX,ζX,0,1.(3.22)

[69.1.2.2] Setting X=Ψ this result is found to be distinctly different from equation (1.3). This shows that finite ensemble scaling (3.22) and finite size scaling (1.3) are not equivalent.