Next Article in Journal
Valorization of Anaerobic Digestate from Biowaste to High-Value Bioproducts: A Review
Previous Article in Journal
ALICE Highlights
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

Grassmannization of the 3D Ising Model †

1
Dipartimento di Fisica e Astronomia “Ettore Majorana”, Università di Catania, Via S. Sofia, 64, I-95123 Catania, Italy
2
Scuola Superiore di Catania, Università di Catania, Via Valdisavoia, 9, I-95123 Catania, Italy
3
Istituto Nazionale di Fisica Nucleare (INFN), Sezione di Catania, Via S. Sofia, 64, I-95123 Catania, Italy
4
Consorzio Nazionale Interuniversitario per la Struttura della Materia (CNISM), Unità di Catania, Via S. Sofia, 64, I-95123 Catania, Italy
5
Istituto per la Microelettronica e Microsistemi (IMM-CNR), Zona Industriale, VIII strada, 5, I-95121 Catania, Italy
6
Department of Physics, Arnold Sommerfeld Center for Theoretical Physics, Theresienstrasse 37, D-80333 Munich, Germany
*
Author to whom correspondence should be addressed.
Presented at the 11th Italian Quantum Information Science conference (IQIS2018), Catania, Italy, 17–20 September 2018.
Current address: School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham B152TT, UK.
Proceedings 2019, 12(1), 20; https://doi.org/10.3390/proceedings2019012020
Published: 6 June 2019
(This article belongs to the Proceedings of 11th Italian Quantum Information Science conference (IQIS2018))

Abstract

:
The application of Feynman’s diagrammatic technique to classical link models with local constraints seems impossible due to (i) the absence of a free Gaussian theory on top of which the perturbative expansion can be constructed, and (ii) Dyson’s collapse argument, rendering the perturbative expansion divergent. However, we show for the classical 3D Ising model how both problems can be circumvented using a Grassmann representation. This makes it possible to obtain an expansion of the spin correlation function and the magnetic susceptibility in terms of the inverse temperature in the thermodynamic limit, through which the values for the critical temperature and critical index γ are evaluated within 1.6% and 5.4% of their accepted values, respectively. Our work is a straightforward adaptation of the theory previously developed in an earlier paper.

1. Introduction

The Ising model was introduced in 1920 by Lenz [1] and solved for the first time by Ising in the 1D case [2]. Nearly 20 years later it was solved exactly in the 2D case by Onsager [3] in the thermodynamic limit. For the 3D Ising model no exact solution is known. One typically makes use of Monte Carlo (MC) simulations in order to find its critical temperature and critical indices (although the recently developed conformal bootstrap method [4] can yield results that are two orders of magnitude more precise). MC methods such as the single spin-flip Metropolis algorithm [5] or the superior worm algorithm [6] are applied to finite-size systems. By means of finite size scaling theory one determines the value of the critical indices by demanding that all the rescaled data lines collapse onto each other in the proximity of the critical point (see Ref. [7] for a pedagogical introduction).
Diagrammatic Monte Carlo simulations by contrast simulate the method directly in the thermodynamic limit. The application of diagrammatic Monte Carlo simulations to the 3D Ising model seems at first questionable because (i) there is no free gaussian theory readily available necessary for the application of Wick’s theorem, and (ii) the nature of classical link variables and local constraints is inherently bosonic such that Dyson’s collapse argument, originally formulated for quantum electrodynamics [8], implies a divergent series. One of us (LP) proposed in Ref. [9], together with M. Kiselev, N. V. Prokof’ev and B. V. Svistunov, a theory in which (pairs of) Grassmann variables [10] represent all degrees of freedom of the high-temperature series expansion of the classical partition function (or, more generally, any discrete representation of the system can be used). The resulting expansion is a Feynman series for the spin correlation function, from which the susceptibility can be computed. The susceptibility χ is known to diverge at the critical point as ( T T c ) γ . Since the series is found to be convergent down to the critical temperature T c , the convergence radius of the power series gives therefore access to the critical temperature, T c , and the critical index, γ . The results presented here have been obtained following the Grassmannization recipe developed by Pollet et al. [9] and by generalizing the code written by Pollet from the 2D case to the 3D case.

2. Results

2.1. Grassmannization of the Ising Model

The approach proposed by Pollet et al. [9] consists of representing the partition function of the 3D Ising model in terms of pairs of Grassmann variables. We do not repeat the development of the formalism here but only discuss the differences of the 3D case compared to the 2D case [9], and introduce notation to make the discussion self-contained. Given a simple cubic lattice, one assigns a set of discrete variables to the links (or, bonds, b), { α b } , with a link factor f ( { α b } ) depending on the occupancy number of the link, and one defines a site factor g j = g ( { α b } j ) that depends only on the links attached to the jth site. The partition function is then written as Z = T r b A b j B j , where:
A b = α 0 exp [ ξ ¯ α , b ξ α , b f ( α ) + ξ ¯ α , b ξ α , b f ( α ) ] and B j = { α b } j exp [ λ ( { α b } j ) b { b } j ξ ^ α b , b ξ ^ α b , b ] .
Here, the variables ξ , ξ ¯ , ξ , and ξ ¯ are Grassmann variables living on the bonds of the lattice. The factors λ ( { α b } j ) are related to the site factors in a non-trivial way and have to be computed. The easiest expansion we can use for such a computation is the high temperature series expansion of the Ising model partition function, first obtained by Kramer and Wannier [11], for a lattice with N sites, in an external field h:
Z = Z 0 { s i } i j ( 1 + s i s j v ) i ( 1 + s i w ) ,
where v = tanh β , w = tanh h and Z 0 = ( cosh β ) 2 N ( cosh h ) N . One finds for the link factors that f ( 0 ) = 1 and f ( 1 ) = v and for the site factors that
g ( 0 ) = g ( 2 ) = g ( 4 ) = g ( 6 ) = 1 ; g ( 1 ) = g ( 3 ) = g ( 5 ) = w .
These are related to the λ factors by the following relations,
λ 1 = w ; λ 2 = 1 w 2 ; λ 3 = 2 w + 2 w 3 ; λ 4 = 2 + 8 w 2 6 w 4 ; λ 5 = 16 w 40 w 3 + 24 w 5 ; λ 6 = 16 136 w 2 + 210 w 4 40 w 6 .
We can now take the thermodynamic limit N and set h = 0 . The partition function can then be written as follows:
Z = Z 0 D [ ξ ¯ ξ ξ ¯ ξ ] { α b } { α b } exp ( ξ ¯ α b ξ α b v + ξ ¯ α b ξ α b v ) j exp ( λ α { b } j { b } j ξ ^ α b ξ ^ α b ) .
Equation (2) is the starting point for the formulation of the Feynman rules. By Taylor-expanding the non-Gaussian term, the partition function takes the form
Z = Z ˜ 0 n = 0 x 1 x n ( + 1 ) n n ! V ( x 1 ) V ( x n ) ,
and can be evaluated by means of Wick’s theorem. The Feynman rules for the spin correlator G ( 0 , r ) at the order n + 1 are:
  • Draw all the n topologically different connected diagrams connecting the origin with site r . Any occupied bond contributes 1 to the expansion order n.
  • The n bare propagators live on the links (they are local), and are represented by a pair of Grassmann variables. A pair contributes a factor v in magnitude to the weight of the diagram. Propagation lines have no arrows.
  • The ( n + 1 ) interaction vertices live on the sites of the lattice, and can be of different type: the origin and end vertex belong to the V 1 , V 3 or V 5 class, the ( n 1 ) others belong to V 2 , V 4 or V 6 class, whose weights are in accordance with Equation (1). All j legs of the vertices V j must be connected by propagator lines.
  • If a link is multiply occupied, a minus sign occurs when swapping 2 Grassmann variables. The minus signs can equivalently be inferred by identifying all fermionic loops.
  • The total weight will be ( 1 ) P ( 1 ) q 1 + q 2 ( 2 ) q 3 + q 4 ( 16 ) q 5 + q 6 v n , being P the signature of the exchange permutation, and q j is the sum of al vertices V j that are of type j, j = 1 , , 6 . It follows that the weight of a diagram is an integer number times v n .

2.2. Results of the simulations

The details of the algorithm are discussed in [9]. The output of the code, which generates and evaluates all possible Feynman diagrams for the spin correlator up to order 10, is a table shown in Table 1.
The sum of the number of diagrams weighted by the symmetry factor of the cubic lattice contributes to the series coefficients for the expansion of the susceptibility. The series for the susceptibility, obtained after a calculation lasting about one hour on a laptop equipped with an Intel Core i5 2.5 GHz processor with a maximum resident set size of about 9 Gbytes, is
β 1 χ = n χ n v n = 1 + 6 v + 30 v 2 + 150 v 3 + 726 v 4 + 3510 v 5 + 16710 v 6 + 79494 v 7 + 375174 v 8 + 1769686 v 9 + 8307678 v 10 + O ( v 11 ) ,
which coincides up to the order under consideration with the result of Campostrini [12].
Since the susceptibility diverges at the critical temperature as ( 1 v / v c ) γ , the convergence radius of the power series for χ can provide information about the critical temperature and the critical index γ for the model. In order to obtain these, it is sufficient to compare the power series coefficients of ( 1 v / v c ) γ with χ n , and to estimate the asymptotic behaviour of the ratio χ n / χ n + 1 ,
χ n χ n + 1 = 1 v c + γ 1 v c 1 n .
A linear fit of the obtained coefficient ratios as a function of 1 / n yields the value of the critical temperature and the value of the critical index, as shown in Figure 1.
The results of the fit are such that ( tanh β c ) 1 = 4.53 ( 4 ) and ( γ 1 ) ( tanh β c ) 1 = 1.4 ( 1 ) , and thus T c / J = 4.45 and γ = 1.30 . These values differ from their known values [12] ( T c / J = 4.51149 ( 4 ) and γ = 1.2371( 4 ) ) by 1.6% and 5.4%, respectively, but agree within error bars.
In summary, we have successfully benchmarked our Grassmannization approach for the 3D Ising model and thereby obtained the high-temperature series expansion for the spin correlator.
Acknowledgements. We wish to thank M. Kiselev, N. V. Prokof’ev, and B. V. Svistunov for stimulating discussions. LP is supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – EXC-2111 – 390814868, and FP7/ERC Consolidator grant QSIMCORR No. 771891.

References

  1. Lenz, W. Beitrag zum Verstandnis der magnetischen Erscheinunge in festen Korpern. Phys. Zs. 1920, 21, 613–615. [Google Scholar]
  2. Ising, E. Beitrag zur theorie des ferromagnetismus. Zeitschrift für Physik 1925, 31, 253–258. [Google Scholar] [CrossRef]
  3. Onsager, L. Crystal Statistics. I. A Two-Dimensional Model with an Order-Disorder Transition. Phys. Rev. (Ser. I) 1944, 65. [Google Scholar] [CrossRef]
  4. El-Showk, S.; Paulos, M.F.; Poland, D.; Rychkov, S.; Simmons-Duffin, D.; Vichi, A. Solving the 3D Ising model with the conformal bootstrap. Phys. Rev. D 2012, 86, 025022. [Google Scholar] [CrossRef]
  5. Metropolis, N.; Rosenbluth, A.W.; Rosenbluth, M.N.; Teller, A.H.; Teller, E. Equation of State Calculations by Fast Computing Machines. J. Chem. Phys. 1953, 21. [Google Scholar] [CrossRef]
  6. Prokof’ev, N.V.; Svistunov, B.V. Worm algorithms for classical statistical models. Phys. Rev. Lett. 2001, 87, 160601. [Google Scholar] [CrossRef] [PubMed]
  7. Sandvik, A.W. Computational studies of quantum spin systems. In AIP Conference Proceedings; AIP: Stockholm, Sweden, 2010; Volume 1297, pp. 135–338. [Google Scholar]
  8. Dyson, F.J. Divergence of Perturbation Theory in Quantum Electrodynamics. Phys. Rev. 1952, 85, 631–632. [Google Scholar] [CrossRef]
  9. Pollet, L.; Kiselev, M.N.; Prokof’ev, N.V.; Svistunov, B.V. Grassmannization of classical models. New J. Phys. 2016, 18, 113025. [Google Scholar] [CrossRef]
  10. Berezin, F.A. The Method of Second Quantization, 1st ed.; Pure and Applied Physics 24; Academic Press: Cambridge, MA, USA, 1966. [Google Scholar]
  11. Kramers, H.A.; Wannier, G.H. Statistics of the Two-Dimensional Ferromagnet. Part I. Phys. Rev. (Ser. I) 1941, 60. [Google Scholar] [CrossRef]
  12. Campostrini, M. Linked-cluster expansion of the Ising model. J. Stat. Phys. 2001, 103, 369–394. [Google Scholar] [CrossRef]
Figure 1. Plot of Equation (3). Black dots are the ratios of the expansion coefficients, while the red solid line is a linear fit. The intercept provides information about T c , while the slope is related to γ . Results are in good accordance with the results found in the literature. The red dashed line is a linear fit through the points omitting the one at n = 1 but it leads to results for T c and γ that deviate from the correct answer.
Figure 1. Plot of Equation (3). Black dots are the ratios of the expansion coefficients, while the red solid line is a linear fit. The intercept provides information about T c , while the slope is related to γ . Results are in good accordance with the results found in the literature. The red dashed line is a linear fit through the points omitting the one at n = 1 but it leads to results for T c and γ that deviate from the correct answer.
Proceedings 12 00020 g001
Table 1. The table shows the number of diagrams connecting the origin with the corresponding site. The second column lists the symmetry factor associated with each site on a cubic lattice.
Table 1. The table shows the number of diagrams connecting the origin with the corresponding site. The second column lists the symmetry factor associated with each site on a cubic lattice.
Site S F ν 0 ν 1 ν 2 ν 3 ν 4 ν 5 ν 6 ν 7 ν 8 ν 9 ν 10
(1,0,0)601040400456061000
(1,1,0)12002016017002144030334
(1,1,1)800060540648088400
(4,0,0)60000104001156024136
(4,1,0)240000050202050060
(4,1,1)24000000300936021474
(4,2,0)24000000150748018647
(4,2,1)480000000105035070
(4,2,2)2400000000420013440
(4,3,0)24000000035022190
(4,3,1)4800000000280011060
(4,3,2)4800000000012600
(4,3,3)2400000000004200
(7,0,0)60000000101120
(7,1,0)240000000080802
(7,1,1)24000000000720
(7,2,0)24000000000360
(7,2,1)480000000000360
(7,3,0)240000000000120
(9,0,0)600000000010
(9,1,0)24000000000010
(10,0,0)600000000001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Martello, E.; Angilella, G.G.N.; Pollet, L. Grassmannization of the 3D Ising Model. Proceedings 2019, 12, 20. https://doi.org/10.3390/proceedings2019012020

AMA Style

Martello E, Angilella GGN, Pollet L. Grassmannization of the 3D Ising Model. Proceedings. 2019; 12(1):20. https://doi.org/10.3390/proceedings2019012020

Chicago/Turabian Style

Martello, E., G. G. N. Angilella, and L. Pollet. 2019. "Grassmannization of the 3D Ising Model" Proceedings 12, no. 1: 20. https://doi.org/10.3390/proceedings2019012020

Article Metrics

Back to TopTop