_{1}

^{*}

The critical limit principle of maximum entropy is put forward, it’s a sufficient condition to obtain accurate critical points, and ensure that the new phase system is still in the maximum entropy state. Two representations for the phase transition of Ising models are found; the universal formula of critical points is explained by thermodynamics. From the point of view of fractal geometry and the correspondence between symmetry and conservation, the scaling laws are reinterpreted. The self consistence equations for the universal class are set up, by which and the scaling laws higher accurate exponents to date are obtained. The temperature where the self similar transformation disappears is calculated.

As a fundamental issue, the critical phenomena theory has pervaded modern physics, the approach to the critical point, the critical fluctuation, and the critical exponents help us understand critical laws deeply. The precise calculation of critical points and critical exponents is the common goal of the researchers. The symmetry analysis that is usually applied in particle physics makes our research to a new state, and has got more quantitative results [

The focus of both views is the lattice correlation length. Whether the correlation length is infinite or finite, we can respond in only one of the two ways. In fact, the concept of infinite correlation length is merely a theoretical assumption. We can’t measure the infinite length; we can only measure the exponent υ that is associated with the length in term of the conventional theory. According to the assumption, the infinite correlation length occurs at the critical temperature T c f , and disappears at the temperature T that is in the vicinity of T c f , and T < T c f . The idea that the ordered state is derived from the lattice correlation is reasonable, but the argument that the correlation length disappears suddenly in the vicinity of T c f is hard to convince us. Let a system have N m lattices in total, a block contain N lattices, the maximum hierarchy number of the SS transformation be r m , we have

r m = ln N m / ln N (1)

where the symbol “ln” represents natural logarithm. If the system is infinite, N m → + ∞ , with the deceasing of temperature the lattice correlation length n increases from n * , which is the block side length at the critical pint, to infinity, so does N at T c s , where the correlation length becomes infinity and the whole system is a block [

r m ∝ | T − T c f | − υ (2)

For the reducible lattice system [

This paper is formed by the above point of view. In Section 2, we put forward the critical limit principle of maximum entropy. We prove again the uniqueness of the critical point. From the dual relation of J and S 2 , we find two representations for the phase transition: spin representation and coupling constant representation. We clarify, in the thermodynamics sense, the meaning of the critical point. From the point of view of the SS transformation, we reinterpret the scaling laws that they are conserved laws corresponding to the symmetry transformation. We set up universal class self consistence constant and the relevant equations for the exponents, and point out that the exponents must satisfy simultaneously the scaling laws and the equations, and the scaling laws are available for the whole renormalization region T c s < T ≤ T c f rather than in the vicinity of T c f .. The structural meaning of the disorder-order phase transition is that the original lattice spin system has become the self-similar block spin system. In Section 3, we calculate T c s and exponents, get a group of exponents, which is the most accurate and self consistent datum to date. Section 4 is conclusion remark. For simplicity, both blocks and sub-blocks are uniformly called blocks in this paper, and the relation between the sub-block’s dimension and the block’s dimension for the reducible lattice system is given by the formula (6.2) of reference [

The conventional theory thinks that since lattice spins are correlated over scales up to the correlation length ζ , it may be plausible to regard the spins with regions up to ζ in size as behaving like a single block spin of side length n, and n ≤ ζ , i.e., the block formation originates from the infinite correlation length. It’s said that one can get accurate critical point only if let n take ζ , if ζ → + ∞ . In fact, this is not true [

We find that there exist simultaneously two sorts of blocks: n + -blocks and n − -blocks, only one of them can become ordered, another disordered [

F = ∑ i = 1 ∞ P i F i (3)

where P i is the occurrence probability of the i-th subsystem. The value of P i can’t change the influence of the singularity of F i on the singularity of F. If P i is determined by Gaussian distribution, F is also the Gaussian distribution free energy, and represented in the form as the Equation (11) of reference [

Critical limit principle of maximum entropy: A thermodynamic equilibrium system keeps itself in the maximum entropy state before disorder-order phase transition, and tends to the critical point in this manner.

According to the fractal theory that a system containing a large plenty of lattice spins becomes ordered as the same as a point spin is a kind of spin contraction mapping, which ensures the uniqueness and existence of self similar transformation [

( 2 − D min ) n 2 + 3 ( 1 − D min ) n − 2 D min = 0 (4)

Equation (4) has only a positive finite real solution n * due to 1 < D min < 2 . For the tetrahedron lattice system, N = ( n 3 + 6 n 2 + 11 n + 6 ) / 6 , see [

( 3 − D min ) n 3 + ( 12 − 6 D min ) n 2 + ( 11 − 11 D min ) n − 6 D min = 0 (5)

This equation has a unique finite real root n * . As the same reason, there is a unique finite side length n * at the critical point for any one system, which accords with the result of numerical calculation [

The energy conservation of transformation produces a dual relation of J and S 2 [

Z J S 2 = 2 D j s 2 (6)

where the lattice spin s 2 = 1 , its coupling constant j, the block fractal dimension is D, its spin square is S 2 , coupling constant J, and coordinate number Z. For a given system Z is constant. Let J = j , we get spin representation:

S 2 = ( 2 D / Z ) s 2 (7)

If let S 2 = s 2 , we have coupling constant representation from Equation (6):

J = ( 2 D / Z ) j (8)

The two representations are equivalent. In the early renormalization group (RG) theory the coupling constant representation is adopted [

The singularity of the free energy gives rise to

Z J S min 2 = k B T c f (9)

This is a universal relation for the critical points of Ising models [

F = 〈 E 〉 − T S e (10)

Only if the entropy S e is maximum can the free energy F be minimized when the T is constant. That the block spin takes the minimum when the system reaches to T c f is the best selection to satisfy the condition. We may view Z J S min 2 as an ordering tendency, k B T c f as a disordering tendency caused by thermal motion. The critical point not only is a balance point for the two tendencies, but also a starting point of the ordering of the system. The fact that the block spin formation is the self-organization behavior and the cooperation phenomena of the lattice spins. Discussing the atomic system, Dirac pointed out that the observed specific heats at ordinary temperatures are given fairly well by a theory that takes into account merely the motion of each atom as a whole and assigns no internal motion to it at all [

Z J S 2 > k B T (11)

It’s the lattice correlation length that determines the size of the block space, i.e., the block side length. The lower the temperature, the longer the correlation length, and the larger the block side length (the block spin) and the higher the ordering degree. By means of boundary condition, Onsager carried the plane square lattice system onto a torus, and obtained a lager critical point than ours [

The T c f is just the critical temperature T c defined by the conventional phase transition theory. We notice that there is no singularity behavior below T c f , i.e., the block structure exists until T c s and the block spins replace the lattice spins. We call the region T c s < T ≤ T c f as renormalization region, rather than merely the criticality one.

The scaling laws are based on the scaling hypothesis, which is purely a conjecture. Taking Ising model as an example, we interpret it by Euclidean geometry and fractal geometry. From the point of view of Euclidean geometry, the reason why the GR theory works is that the lattice correlation length becomes infinity in the vicinity of T c f , such that the influence of all finite microscopic characteristic lengths are wiped out, and the correlation length is a unique characteristic quantity. The singularities of thermodynamic parameters are attributed to the correlation length singularity, and any size transformation can’t change the form of the free energy function, only changes the parameter scaling. On the other hand, the transformation hierarchy is emphasized in terms of fractal geometry: The formation of blocks with finite side length n begins at T c f , they can go through r-hierarchical SS transformations at the same temperature, and the original lattice spin system is replaced by the block spin system. For any finite value of r, one can find the same form of the free energy function represented by the block spins, no matter what value of r will be. Clearly, both interpretations describe the same phenomena. It’s impossible to explain in terms of Euclidean geometry the uniqueness of the block side length n * at the critical point and the coexistence of n + -blocks and n − -blocks. In fact, the renormalization transformation is just the self similar one. In a word, the geometrical structure and the physical parameters function forms on the (r + 1)-th hierarchy maintain the original ones on the r-hierarchy. By the forms, we are unable to recognize which belong to the r-hierarchy and what is the actual value of r. The indistinguishability means symmetry [

For a universal class of space dimension d and order parameter dimension m, there are two independent conserved constants: p and q, they can be considered as the magnitudes of two orthogonal constant vectors p and q, respectively. According to Wilson theory [

p ∝ ln ( 2 D / Z ) / ln N , q ∝ ln ( 2 D / Z ) / ln N (12)

where the elementary excitation of block spins is badly neglected [

p = σ ( d , m ) q (13)

where σ ( d , m ) is a universal class self consistence constant. Three parameters in this equation, only two of them are independent. p and q control six exponents:

α = 2 − 1 / p , β = ( 1 − q ) / p , γ = ( 2 q − 1 ) / p (14)

δ = q / ( 1 − q ) , η = d ( 1 − 2 q ) + 2 , υ = 1 / p d (15)

They determine four scaling laws:

α + 2 β + γ = 2 , γ = β ( δ − 1 ) (16)

γ = υ ( 2 − η ) , υ d = 2 − α (17)

Using Equations (13)-(17), there should theoretically be 15 self consistence equations for σ ( d , m ) by the combinations of any two different exponents. But, the combinations of α and υ , η and δ have no consequence, 13 equations are actually established:

σ 1 = 1 / ( 2 − α − β ) , σ 2 = 1 / ( β + γ ) , σ 3 = 2 d / [ ( 2 − α ) ( 2 + d − η ) ] (18)

σ 4 = ( δ + 1 ) / δ ( 2 − α ) , σ 5 = 2 / ( γ + 2 − α ) , (19)

σ 6 = 1 / ( β δ ) , σ 7 = 1 / ( d υ − β ) , σ 8 = ( δ − 1 ) / ( δ γ ) , σ 9 = 2 / ( γ + d υ ) (20)

σ 10 = ( 4 − 2 η ) / γ ( d + 2 − η ) , σ 11 = 2 d / [ υ d ( d + 2 − η ) ] (21)

σ 12 = ( 1 + δ ) / ( υ d δ ) , σ 13 = ( d − 2 + η ) / β ( d + 2 − η ) (22)

The conservation requires that all exponents should not only obey the scaling laws, but also satisfy the self consistence equations, and σ i = σ ( d , m ) for all i, i = 1 , 2 , ⋯ , 13 . Our experience indicates that the scaling laws are not sensitive to the errors of the exponents: if the deviation between the exponents values and a law is 0.0001, the difference between the high accurate value of the exponent and its crude value is usually greater than 0.0001. As this reason, the exponents must strictly obey the laws. In most cases, the self consistence equations have higher sensitivity to the exponents values: the exponents that can comply with the laws are not necessarily accord with them. Only those exponents that satisfy both the laws and the self consistence equations have higher accuracy.

The experimental datum illustrate that all order parameters curves are smoothly continuous when the temperature is below T c f [

( T y − T ) ± x (23)

where T y is in the normalization region, T y > T , and T in the vicinity of T y , x = α , β , γ , δ , η , υ , their positive and negative signs are the same as the ones’ at T c f − 0 , which is in the vicinity of T c f , and T c f − 0 < T c f . For example, by Equations ((1), (2), and (23)), we get

ln N ∝ ( ln N m ) ( T y − T ) υ (24)

N increases with the decreasing of T. Equation (24) verifies that the power law is a characteristic of fractals [

Whether it is irreducible system or reducible system [

〈 E 〉 = ∫ 0 x D x 2 e x − 1 d x = 0 (25)

where x = ℏ ω / k B T , x D = ℏ ω D / k B T c s , ℏ is Plank constant, ω the magnetic phonon frequency, ω D the Debye frequency. Generally, we may get T c s from the integral equation. If e x ≫ 1 , e x − 1 ≈ e x . Equation (25) simplifies as

∫ 0 x D x 2 e − x d x = 0 (26)

This leads to an algebraic equation

x D 2 e − x D + 2 x D e − x D + 2 e − x D − 2 = 0 (27)

We’ll get T c s from the solution. If e x ≈ 1 , Equation (25) becomes

∫ 0 x D x d x = 0 (28)

x D = 0 , the solution is unreasonable.

The elementary excitation of block spins also contributes to the free energy nonsingularity [

d,m | σ ( d , m ) | q | p | α | β | γ | δ | η | υ |
---|---|---|---|---|---|---|---|---|---|

2.1 | 8/15 | 15/16 | 1/2 | 0 | 1/8 | 7/4 | 15 | 1/4 | 1 |

3.1 | 0.638613 | 0.827204 | 0.528263 | 0.1070 | 0.3271 | 1.2388 | 4.7872 | 0.0368 | 0.6310 |

3.2 | 0.610231 | 0.817227 | 0.498697 | −0.0052 | 0.3665 | 1.2722 | 4.4713 | 0.0966 | 0.6684 |

3.3 | 0.583511 | 0.826243 | 0.482122 | −0.0742 | 0.3604 | 1.3534 | 4.7552 | 0.0425 | 0.6914 |

constant σ i accurate to the fourth decimal place, so as to let the extra errors lie as behind as possibly the fourth decimal point. In principle, if we want the exponents and the constant σ i to be accurate to the x-th decimal place, the constant σ ( d , m ) , p, and q should be to the (x + 2)-th decimal place. For example, the following exponents and σ i for the universal class of d = 3 and m = 1 are accurate to the sixth decimal point, x = 6: α = 0.107063 , β = 0.327093 , γ = 1.238571 , δ = 4.787153 , η = 0.036779 , υ = 0.630979 , σ i = 0.638633 ( i = 1 , 2 , ⋯ , 13 ) ; and the others to the eighth decimal place: σ ( 3 , 1 ) = 0.63863311 , p = 0.52827948 , q = 0.82720347 . This group of exponents is the most accurate and self consistent datum to date. More accurate exponents will be obtained by the same method.

The critical limit principle of maximum entropy introduces a correct way to obtain accurate critical points. The reinterpretation of the scaling laws from the correspondence between symmetry and conservation gives us a new understanding of scaling invariance. Based on the experimental and theoretical datum given by other authors, and applying golden section method and numerical approximation in the scaling laws and the self consistence equations, we can get higher accurate exponents.

The author declares no conflicts of interest regarding the publication of this paper.

Feng, Y.G. (2018) Critical Limit Principle of Maximum Entropy and Reinterpretation of Scaling Laws. Journal of Modern Physics, 9, 2035-2044. https://doi.org/10.4236/jmp.2018.911128