Sie sind hier: ICP » R. Hilfer » Publikationen

IV Finite size scaling

[69.1.3.1] This section discusses the implications of finite ensemble scaling for finite size scaling at a critical point. [69.1.3.2] Contrary to finite ensemble scaling the theory of finite size scaling includes the strongly correlated microscopic cell variables into the theoretical consideration. [69.2.0.1] This can be done in two ways. [69.2.0.2] Thermodynamic finite size scaling concentrates on the thermodynamic fluctuations within the ensemble, while statistical mechanical (or fieldtheoretical) finite size scaling focusses on the correlation functions on the block level. [69.2.0.3] The distinction appears already in equations (1.1) and (1.2). [69.2.0.4] The general identification of thermodynamics as the infinite volume limit of statistical mechanics implies a relation between the two parts which is at the origin of hyperscaling relations.

IV.A Thermodynamic finite size scaling

[69.2.1.1] The thermodynamic method of reintroducing the strongly correlated cell variables is to to use the definition of block variables (2.18) and to define

ZMNφ=XMNφ-CMNDMN=j=1Ni=1MXMNφyj+xi-CMNDMN(4.1)

as a double sum over correlated microscopic cell variables. [69.2.1.2] Although the microscopic variables are strongly correlated inside the blocks they remain uncorrelated at separations larger than ξ. [69.2.1.3] Therefore the property of strong mixing [33, 34] continues to hold in the ensemble limit. [69.2.1.4] Therefore the same considerations as in the previous section can also be applied to the double sums (4.1) to give the finite size scaling result

PXMNxHx;ϖX,ζX,0,DDMNϖX(4.2)

where now

DNM=MNΛMN1/ϖX(4.3)

similar to equation (3.8).

[69.2.1.5] To exhibit the relation of the result (4.2) with the usual thermodynamic finite size scaling Ansatz (1.3) [8] for the order parameter distribution it is first necessary to rewrite the results in terms of the probability density for the ensemble averages X¯MNφ=XMNφ/MN. [69.2.1.6] This gives the thermodynamic finite size scaling result

p¯X¯MNxL/ad-d/ϖXΛL/ad1/ϖXhxL/ad-d/ϖXΛL/ad1/ϖX;ϖX,ζX,0,D.(4.4)

[69.2.1.7] Setting X=Ψ and comparing with [5] yields the identification [21, 22]

ϖΨ=min2,γΨΨ+2βΨγΨΨ+βΨ=min2,λΨ(4.5)

where γΨΨ is the order parameter susceptibility exponent, βΨ is the order parameter exponent, and λΨ is the [page 70, §0]    generalized Ehrenfest order [19] in the conjugate field direction. [70.1.0.1] The appearance of the min-function results from the general inequality (3.7). [70.1.0.2] Similarly for the energy density X=E the result

ϖE=min2,2-αE=min2,λE(4.6)

is obtained with αE=α the specific heat exponent. [70.1.0.3] In general the identification is given as ϖX=min2,2-αX=min2,λX where αX is the thermodynamic fluctuation exponent [9] defined in terms of derivatives of the free energy. [70.1.0.4] Equation (4.4) in combination with equations (4.5) and (3.16) determines the thermodynamic finite size scaling function for the order parameter distribution in (1.3) explicitly as

p~Ψx,y=Rx,yhΨx+HΨxRx,yx    :    for  x0Rx,yhΨx-1-HΨxRx,yx    :    for  x>0(4.7)

where

hΨx,y=hx;γΨΨ+2βΨγΨΨ+βΨ,ζΨ,0,D=dHΨxdx(4.8)

and h is defined through the H-functions in equations (3.12), (3.13) and the appendix. [70.1.0.5] Note that the thermodynamic finite size scaling function depends on y only through the nonuniversal cutoff function Rx;M,y,c. [70.1.0.6] It will be seen below that the dependence on y in the universal function h reappears in fieldtheoretical finite size scaling. [70.1.0.7] Note also that the general inequalities βΨ>0 and γΨΨ>0 imply ϖΨ<2.

IV.B Fieldtheoretical finite size scaling

[70.1.1.1] The fieldtheoretical or statistical mechanical method of reintroducing the microscopic cell variables uses the same uncorrelated block sums as in finite ensemble scaling (3.2), but multiplies them with the M-dependent field theoretic renormalization factor for block sums DM/M from (2.17) which has to be calculated independently. [70.1.1.2] In this case the renormalized ensemble sums are defined as

ZMNφ=DM/MYMNφ-CNDN=DM/Mj=1NYMNφyj-CNDN(4.9)

where CN and DN are constants as in (3.2). [70.1.1.3] The composite operators YMNφyj have been denoted differently from the the thermodynamic case to indicate that the variables of interest in mesoscopic fieldtheoretic or statistical mechanical calculations (block level) may in general differ from those accessible to macroscopic thermodynamic experiments (ensemble level). [70.2.0.1] Particular examples are the staggered magnetization for antiferromagnets or the quantum mechanical wave function. [70.2.0.2] Applying the same limit theorem as in the previous section now gives the fieldtheoretic finite size scaling result

PYMNxHx;ϖY,ζY,0,DMDNDMϖY(4.10)

for the limiting probability distribution function of ensemble sums in the ensemble limit. [70.2.0.3] Using (2.13) and going over to averages the finite size scaling form for the probability density of ensemble averages is found as

p¯Y¯MNxLadYhxLadY;ϖY,ζY,0,DLξd-ϖYd-dYΛLξ(4.11)

which is exactly of the form (1.3) with d*=0. [70.2.0.4] Thus the validity of (2.13), which has to be established by independent calculation, implies the validity of hyperscaling. [70.2.0.5] Note that the fieldtheoretic finite size scaling result (4.11) appears to be different from the thermodynamic one (4.4) in that it depends on L/ξ. [70.2.0.6] It will be seen below however that the two forms are generally identical except for ϖX=2.

IV.C Hyperscaling and the Structure of the Gaussian Fixed Point

[70.2.1.1] To establish the connection between thermodynamic fluctuation exponents ϖX and fieldtheoretic correlation exponents ϖY it is necessary to compare the scaling results (4.4) and (4.11). [70.2.1.2] Note that (4.4) holds generally by virtue of the ensemble limit while the validity of (4.11) depends upon the validity of (2.13). [70.2.1.3] The connection between thermodynamics and statistical mechanics is generally given by identifying-logZ with the free energy or, in the microcanonical ensemble, by inverting the logarithm of the density of states to give the internal energy as function of entropy. [70.2.1.4] Thus the identification rests upon the identification of microscopic and macroscopic energies. [70.2.1.5] In fact the energy is the only observable which will always exist microscopically and macroscopically for thermal systems because it is a defining property of the system, and generates the thermal fluctuations of interest. [70.2.1.6] Thus the connection between thermodynamics and statistical mechanics in the present probabilistic approach is provided by identifying equation (4.4) for X=H with equation (4.11) for Y=H. [70.2.1.7] This yields the algebraic form

DM=M1-1/ϖEDΛNDΛMN1/ϖE(4.12)

for the energy renormalization. [70.2.1.8] Comparison with (2.13) and (4.6) gives the identification (first obtained in [21, 22]) [page 71, §0]

ϖE=min2,dd-dE=min2,dν=min2,2-α(4.13)

where ν=νE is the correlation length exponent, and α=αE is the specific heat exponent. [71.1.0.1] Thus equation (4.12) and (4.13) combined with the general relation [9]

2-αXνX=2-αYνY(4.14)

establish the general validity of hyperscaling for all microscopically and macroscopically accessible observables whenever the specific heat exponent is positive. [71.1.0.2] Therefore the hyperscaling relation

ϖX=min2,dd-dX=min2,dνX=min2,2-αX(4.15)

holds for all phase transitions with α>0. [71.1.0.3] This result is a direct consequence of identifying thermal fluctuations in thermodynamics with those in statistical mechanics or field theory.

[71.1.1.1] The violation of hyperscaling above four dimensions in field theory is now a simple consequence of the renormalization group eigenvalues yE=1/νE=2 and yΨ=1/νΨ=d+2/2 for the Gaussian fixed point. [71.1.1.2] Equation (4.13) implies ϖE=2 at d=4.

[71.1.2.1] Of course the present theory does not allow to conclude that hyperscaling is generally violated for α0. [71.1.2.2] In fact very often hyperscaling continues to be valid in such cases. [71.1.2.3] To see how this is possible it is instructive to consider the domains of attraction for the stable laws appearing in the finite size and finite ensemble scaling formulas. [71.1.2.4] Within the present approach the fact that only stable distributions have nonempty domains of attraction [28] is the reason for the existence of fixed points in the renormalization group picture and for universality of critical behaviour [35]. [71.1.2.5] It is well known [27, 28] that the domain of attraction is very different for gaussian and nongaussian fixed points.

[71.1.3.1] The existence of the limit distribution in (4.2) for the correlated ensemble sums implies by virtue of (2.18) and (2.19) that the limiting distribution of the correlated block sums

PXMNj(x)=Prob{XMN(φ(yj))x}(4.16)

must approach a distribution within the domain of attraction of the stable distribution (4.2) for all blocks j=1,,N. [71.1.3.2] In order that a distribution PXMNjx belongs to the domain of attraction of the stable law with index 0<ϖX<2 and parameters ζX,D it is necessary and sufficient [28] that, as x,

PXMNjx=c-Λ-x-x-ϖX    :    for  x<01-c+Λxx-ϖX    :    for  x>0(4.17)

where Λx is slowly varying and the constants c-,c+0,c-+c+>0 are related to the parameters ϖX,ζX,D by

c±=Dcosω1ζX2Γ1-ϖXcosω11cotω1tanω1ζX    :    for  0<ϖx<1DπcosπζX/21cotω1tanπζX/2    :    for  ϖX=1D1-ϖXcosω2ζX2Γ2-ϖXcosω11cotω1tanω2ζX    :    for  1<ϖx<2(4.18)

with ω1=πϖX/2 and ω2=π2-ϖX/2. [71.2.0.1] For ϖX=2 on the other hand the domain of attraction is much larger. A distribution PXMNjx belongs to the domain of attraction of the Gaussian if it has a finite variance or if, for x>0,

1-PXMNjx+PXMNj-x=x-2Λx(4.19)

where Λx is slowly varying.

[71.2.1.1] Equation (4.17) implies that for ϖX<2 the generalized susceptibility which is proportional to the second moment of the renormalized block variables

=limM,N-x2dPXMNjxχXX(4.20)

diverges in each block j=1,,N. [71.2.1.2] For ϖX=2 on the other hand the second moment may either diverge or else it is finite and nonzero. [71.2.1.3] (A zero value occurs only away from the critical point). [71.2.1.4] This result underlines the general validity of the algebraic form (2.13) derived in (4.12) for nongaussian fixed points, i.e. ϖE<2, which then implies the validity of hyperscaling. [71.2.1.5] The Gaussian fixed point ϖE=2 on the other hand has a much larger domain of attraction. [71.2.1.6] In particular it contains both distribution functions with algebraic tails and distributions without algebraic tails. [71.2.1.7] No general conclusion about the validity or violation of hyperscaling can be drawn in the present approach for the Gaussian fixed point.