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φ→y→j 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φ→y→j-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,N→∞N/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 ϖX≠1 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<ϖX≤2-1≤ζX≤1-∞<C<∞0≤D. | | (3.7) |
[68.1.0.5] If the limit exists, and D≠0, the constants DN must
have the form
where ΛN is a slowly varying function [28], i.e.
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 ϖX≠1 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
ϖX≠1) 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 ζX≠1 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 x≤01-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,N→∞N/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
PXMNx≈Hx-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 ϖX≠1-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]
pXMNx≈hx;ϖ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¯MNx≈L/ξ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.