No Arabic abstract
We have made substantial advances in elucidating the properties of the susceptibility of the square lattice Ising model. We discuss its analyticity properties, certain closed form expressions for subsets of the coefficients, and give an algorithm of complexity O(N^6) to determine its first N coefficients. As a result, we have generated and analyzed series with more than 300 terms in both the high- and low-temperature regime. We quantify the effect of irrelevant variables to the scaling-amplitude functions. In particular, we find and quantify the breakdown of simple scaling, in the absence of irrelevant scaling fields, arising first at order |T-T_c|^{9/4}, though high-low temperature symmetry is still preserved. At terms of order |T-T_c|^{17/4} and beyond, this symmetry is no longer present. The short-distance terms are shown to have the form (T-T_c)^p(log|T-T_c|)^q with p ge q^2. Conjectured exact expressions for some correlation functions and series coefficients in terms of elliptic theta functions also foreshadow future developments.
We investigate the role of a transverse field on the Ising square antiferromagnet with first-($J_1$) and second-($J_2$) neighbor interactions. Using a cluster mean-field approach, we provide a telltale characterization of the frustration effects on the phase boundaries and entropy accumulation process emerging from the interplay between quantum and thermal fluctuations. We found that the paramagnetic (PM) and antiferromagnetic phases are separated by continuous phase transitions. On the other hand, continuous and discontinuous phase transitions, as well as tricriticality, are observed in the phase boundaries between PM and superantiferromagnetic phases. A rich scenario arises when a discontinuous phase transition occurs in the classical limit while quantum fluctuations recover criticality. We also find that the entropy accumulation process predicted to occur at temperatures close to the quantum critical point can be enhanced by frustration. Our results provide a description for the phase boundaries and entropy behavior that can help to identify the ratio $J_2/J_1$ in possible experimental realizations of the quantum $J_1$-$J_2$ Ising antiferromagnet.
The competition between interactions and dissipative processes in a quantum many-body system can drive phase transitions of different order. Exploiting a combination of cluster methods and quantum trajectories, we show how the systematic inclusion of (classical and quantum) nonlocal correlations at increasing distances is crucial to determine the structure of the phase diagram, as well as the nature of the transitions in strongly interacting spin systems. In practice, we focus on the paradigmatic dissipative quantum Ising model: in contrast to the non-dissipative case, its phase diagram is still a matter of debate in the literature. When dissipation acts along the interaction direction, we predict important quantitative modifications of the position of the first-order transition boundary. In the case of incoherent relaxation in the field direction, our approach confirms the presence of a second-order transition, while does not support the possible existence of multicritical points. Potentially, these results can be tested in up-to date quantum simulators of Rydberg atoms.
We calculate very long low- and high-temperature series for the susceptibility $chi$ of the square lattice Ising model as well as very long series for the five-particle contribution $chi^{(5)}$ and six-particle contribution $chi^{(6)}$. These calculations have been made possible by the use of highly optimized polynomial time modular algorithms and a total of more than 150000 CPU hours on computer clusters. For $chi^{(5)}$ 10000 terms of the series are calculated {it modulo} a single prime, and have been used to find the linear ODE satisfied by $chi^{(5)}$ {it modulo} a prime. A diff-Pade analysis of 2000 terms series for $chi^{(5)}$ and $chi^{(6)}$ confirms to a very high degree of confidence previous conjectures about the location and strength of the singularities of the $n$-particle components of the susceptibility, up to a small set of ``additional singularities. We find the presence of singularities at $w=1/2$ for the linear ODE of $chi^{(5)}$, and $w^2= 1/8$ for the ODE of $chi^{(6)}$, which are {it not} singularities of the ``physical $chi^{(5)}$ and $chi^{(6)},$ that is to say the series-solutions of the ODEs which are analytic at $w =0$. Furthermore, analysis of the long series for $chi^{(5)}$ (and $chi^{(6)}$) combined with the corresponding long series for the full susceptibility $chi$ yields previously conjectured singularities in some $chi^{(n)}$, $n ge 7$. We also present a mechanism of resummation of the logarithmic singularities of the $chi^{(n)}$ leading to the known power-law critical behaviour occurring in the full $chi$, and perform a power spectrum analysis giving strong arguments in favor of the existence of a natural boundary for the full susceptibility $chi$.
An important problem in statistical physics concerns the fascinating connections between partition functions of lattice models studied in equilibrium statistical mechanics on the one hand and graph theoretical enumeration problems on the other hand. We investigate the nature of the relationship between the number of spanning trees and the partition function of the Ising model on the square lattice. The spanning tree generating function $T(z)$ gives the spanning tree constant when evaluated at $z=1$, while giving he lattice green function when differentiated. It is known that for the infinite square lattice the partition function $Z(K)$ of the Ising model evaluated at the critical temperature $K=K_c$ is related to $T(1)$. Here we show that this idea in fact generalizes to all real temperatures. We prove that $ ( Z(K) {rm sech~} 2K ~!)^2 = k expbig[ T(k) big] $, where $k= 2 tanh(2K) {rm sech}(2K)$. The identical Mahler measure connects the two seemingly disparate quantities $T(z)$ and $Z(K)$. In turn, the Mahler measure is determined by the random walk structure function. Finally, we show that the the above correspondence does not generalize in a straightforward manner to non-planar lattices.
Energy eigenvalues and order parameters are calculated by exact diagonalization for the transverse Ising model on square lattices of up to 6x6 sites. Finite-size scaling is used to estimate the critical parameters of the model, confirming universality with the three-dimensional classical Ising model. Critical amplitudes are also estimated for both the energy gap and the ground-state energy.