Ising model has been successful in describing ferromagnetism and its phase transition to paramagnet. At the critical point, the free energy density function and its derivatives diverge. Their behaviour near the critical point are described by power-laws with associated critical exponents. In many cases, the critical exponents can be determined from analytic solutions via conventional renormalization groups methods or from Monte Carlo simulations. However, for quantum many-body systems, very few are tractable to analytical solutions. The quantum many-body wavefunction belongs to large dimensional Hilbert space that increases exponentially with system size. If the Hamiltonian is gapped and only local interaction is considered, then the wavefunction can be efficiently truncated. Tensor Network formalism provides a scheme to truncate the less important degrees of freedom via Singular Value Decomposition (SVD) of the density matrix. In this study, we investigated the thermodynamic properties and phase transition of one-dimensional transverse-field quantum Ising model (1D-tQIM) under the finite-size effect and random coupling strength. Starting with Matrix Product States (MPS) as a wavefunction ansatz, the Density Matrix Renormalization Group algorithm is applied to the MPS. The variational algorithm, which iteratively performs SVDs and truncation at each bond, approximates the ground state MPS wavefunction. All quantum observables are calculated from the contraction of the resultant ground state MPS. Although theoretically, divergence at critical points only happen in an infinite system, one can obtain the critical exponents through simulation of finite-size 1D-tQIM. Using the analytic solution as a benchmark, we compared the finite-size effects of the system using finite-size scaling analysis and MPS methods. The critical exponents of 1D-tQIM are independently calculated and compared with the analytical results. Thermodynamic quantities such as magnetization, susceptibility and correlation function are calculated for system sizes of 20, 40, 60, 80, 100 and 120 spins. We determined the respective critical exponents:
The random transverse-field Ising ferromagnet (RTFIF) is a highly disordered quantum system which contains randomness in the coupling strengths as well as in the transverse-field strengths. In one dimension, the critical properties are governed by an infinite-randomness fixed point (IRFP), and renormalization-group studies argue that the two-dimensional (2D) model is also governed by an IRFP. However, even the location of the critical point remains unsettled among quantum Monte Carlo (QMC) studies. In this work, we perform extensive QMC simulations to locate the quantum critical point and attempt a finite-size scaling analysis to observe the critical behavior. We estimate the critical field strength of the 2D RTFIF as ${\mathrm{\ensuremath{\Gamma}}}_{c}=7.52(2)$, together with critical exponents such as $\ensuremath{\beta}=1.5(3), \ensuremath{\nu}=1.6(3)$, and
z=3.3(3) or $\ensuremath{\psi}=0.50(3)$. We have also considered the McCoy-Wu model, which has randomness in the ferromagnetic coupling strengths but not in the transverse-field strength. Our QMC calculation shows that the critical behavior of the 2D McCoy-Wu model is closer to that of the 2D transverse-field Ising spin glass than to that of the 2D RTFIF. These numerical findings enhance our understanding of disordered 2D quantum systems.
We study the ferromagnetic transverse-field Ising model with quenched disorder at T = 0 T = 0 in one and two dimensions by means of stochastic series expansion quantum Monte Carlo simulations using a rigorous zero-temperature scheme. Using a sample-replication method and averaged Binder ratios, we determine the critical shift and width exponents \nu_\mathrm{s} ν s and \nu_\mathrm{w} ν w as well as unbiased critical points by finite-size scaling. Further, scaling of the disorder-averaged magnetisation at the critical point is used to determine the order-parameter critical exponent \beta β and the critical exponent \nu_{\mathrm{av}} ν a v of the average correlation length. The dynamic scaling in the Griffiths phase is investigated by measuring the local susceptibility in the disordered phase and the dynamic exponent z' z ′ is extracted. By applying various finite-size scaling protocols, we provide an extensive and comprehensive comparison between the different approaches on equal footing. The emphasis on effective zero-temperature simulations resolves several inconsistencies in existing literature.
We propose a scheme to perform tensor network based finite-size scaling analysis for two-dimensional classical models. In the tensor network representation of the partition function, we use higher-order tensor renormalization group (HOTRG) method to coarse grain the weight tensor. The renormalized tensor is then used to construct the approximated transfer matrix of an infinite strip of finite width. By diagonalizing the transfer matrix we obtain the correlation length, the magnetization, and the energy density, which are used in finite-size scaling analysis to determine the critical temperature and the critical exponents. As a benchmark we study the two-dimensional classical Ising model. We show that the critical temperature and the critical exponents can be accurately determined. With HOTRG bond dimension
D=70, the absolute errors of the critical temperature
Tc and the critical exponent $\ensuremath{\nu}$, $\ensuremath{\beta}$ are at the order of ${10}^{\ensuremath{-}7},{10}^{\ensuremath{-}5}$, ${10}^{\ensuremath{-}4}$, respectively. Furthermore, the results can be systematically improved by increasing the bond dimension of the HOTRG method. Finally, we study the length scale induced by the finite cutoff in bond dimension and elucidate its physical meaning in this context.
In this work, we benchmark the well-controlled and numerically accurate exponential thermal tensor renormalization group (XTRG) in the simulation of interacting spin models in two dimensions. Finite temperature introduces a finite thermal correlation length $\ensuremath{\xi}$, such that for system sizes $L\ensuremath{\gg}\ensuremath{\xi}$ finite-size calculations actually simulate the thermodynamic limit. In this paper, we focus on the square lattice Heisenberg antiferromagnet (SLH) and quantum Ising models (QIM) on open and cylindrical geometries up to width
W=10. We explore various one-dimensional mapping paths in the matrix product operator (MPO) representation, whose performance is clearly shown to be geometry dependent. We benchmark against quantum Monte Carlo (QMC) data, yet also the series-expansion thermal tensor network results. Thermal properties including the internal energy, specific heat, and spin structure factors, etc. are computed with high precision, obtaining excellent agreement with QMC results. XTRG also allows us to reach remarkably low temperatures. For SLH, we obtain an energy per site ${u}_{g}^{*}\ensuremath{\simeq}\ensuremath{-}0.6694(4)$ and a spontaneous magnetization ${m}_{S}^{*}\ensuremath{\simeq}0.30(1)$ already consistent with the ground-state properties, which is obtained from extrapolated low-
T thermal data on $W\ensuremath{\le}8$ cylinders and $W\ensuremath{\le}10$ open strips, respectively. We extract an exponential divergence versus
T of the structure factor
S(M), as well as the correlation length $\ensuremath{\xi}$, at the ordering wave vector $M=(\ensuremath{\pi},\ensuremath{\pi})$, which represents the renormalized classical behavior and can be observed over a narrow but appreciable temperature window, by analyzing the finite-size data by XTRG simulations. For the QIM with a finite-temperature phase transition, we employ several thermal quantities, including the specific heat, Binder ratio, as well as the MPO entanglement to determine the critical temperature
Tc.
We investigate the quantum phase transition in the transverse-field Ising model on the Sierpiński gasket using finite-size scaling (FSS) and numerical renormalization group (NRG). Since next generations of the fractal lattice contain exponentially more spins, which in turn increase exponentially the Hilbert space dimension, we challenge and prove usefulness of small systems in FSS. We identified a quantum critical point at
λc≈2.63−2.93, with critical exponents $ ν\approx 0.64 - 0.71, β\approx 0.30, γ\approx 1.67$ and
z≈1.33. The numerical renormalization group method produced results consistent with finite-size scaling approach (
λc=2.766,β=0.316), supporting our findings. Compared to the values reported so far in the literature, critical field is in a strong disagreement, while exponents are generally similar excluding
β and
γ. However, it should be noted, that the lattice investigated in previous works is different from ours, while the latter is in our opinion the standard Sierpiński gasket.
Phase transitions in dissipative quantum systems have been investigated using various analytical approaches, particularly in the mean-field (MF) limit. However, analytical results often depend on specific methodologies. For instance, the Keldysh formalism shows that the dissipative transverse Ising (DTI) model exhibits a discontinuous transition at the upper critical dimension,
dc=3; however, the fluctuationless MF approach shows that a continuous transition occurs in infinite dimensions (
d∞). Thus, analytical results must be confirmed by numerical verification. However, this may not be feasible owing to the exponential increase in computational complexity as
O(22N) with system size
N. However, if we assume that the spins are fully connected in infinite dimensions (
d∞), the spin indices are permutation-symmetric. This symmetry significantly reduces the complexity, to
O(N3). Using this symmetry, we simulate the DTI model on fully connected networks with large system sizes. Using finite-size scaling analysis, we show that the upper critical dimension of the DTI model is not
3 but a non-integer dimension
3.5; thus, the discontinuous transition that occurs at high dissipation strength is not an MF result, but is induced by fluctuation effects below
dc. We conclude that the DTI model with strong dissipation exhibits a continuous transition above
dc, as predicted using the fluctuationless MF approach, and we determine the critical exponents.