More is different.
God is in the details.

NEWS

Title
Neural network approach to scaling analysis of critical phenomena
Reference
Physical Review E 107, 044128 (2023)
DOI
10.1103/PhysRevE.107.044128
Author
Ryosuke Yoneda and Kenji Harada
Abstract
Determining the universality class of a system exhibiting critical phenomena is one of the central problems in physics. There are several methods to determine this universality class from data. As methods to collapse plots onto scaling functions, polynomial regression, which is less accurate, and Gaussian process regression, which provides high accuracy and flexibility but is computationally expensive, have been proposed. In this paper, we propose a regression method using a neural network. The computational complexity is linear only in the number of data points. We demonstrate the proposed method for the finite-size scaling analysis of critical phenomena in the two-dimensional Ising model and bond percolation problem to confirm the performance. This method efficiently obtains the critical values with accuracy in both cases.
Comments
10 pages, 10 figures
Preprint
arXiv.2209.01777

Title
Quantum critical dynamics in the two-dimensional transverse Ising model
Reference
Physical Review Research 5, 013186 (2023)
DOI
10.1103/PhysRevResearch.5.013186
Author
Chisa Hotta, Tempei Yoshida, and Kenji Harada
Abstract
In the vicinity of the quantum critical point (QCP), thermodynamic properties diverge toward zero temperature governed by universal exponents. Although this fact is well known, how it is reflected in quantum dynamics has not been addressed. The QCP of the transverse Ising model on a triangular lattice is an ideal platform to test the issue, since it has an experimental realization, the dielectrics being realized in an organic dimer Mott insulator, κ−ET2X, where a quantum electric dipole represents the Ising degrees of freedom. We track the Glauber-type dynamics of the model by constructing a kinetic protocol based on the quantum Monte Carlo method. The dynamical susceptibility takes the form of the Debye function and shows a significant peak narrowing in approaching a QCP due to the divergence of the relaxation timescale. It explains the anomaly of dielectric constants observed in the organic materials, indicating that the material is very near the ferroelectric QCP. We disclose how the dynamical and other critical exponents develop near QCP beyond the simple field theory.
Comments
12 pages, 8 figures
Preprint
arXiv:2209.11599

Title
Quantum Circuit Simulation by SGEMM Emulation on Tensor Cores and Automatic Precision Selection
Reference
ISC 2023
Author
Hiryuki Ootomo, Hidetaka Manabe, Kenji Harada, and Rio Yokota
Abstract
Quantum circuit simulation provides the foundation for the development of quantum algorithms and the verification of quantum supremacy. Among the various methods for quantum circuit simulation, tensor network contraction has been increasing in popularity due to its ability to simulate a larger number of qubits. During tensor contraction, the input tensors are reshaped to matrices and computed by a GEMM operation, where these GEMM operations could reach up to 90\% of the total calculation time. GEMM throughput can be improved by utilizing mixed-precision hardware such as Tensor Cores, but straightforward implementation results in insufficient fidelity for deep and large quantum circuits. Prior work has demonstrated that compensated summation with special care of the rounding mode can fully recover the FP32 precision of SGEMM even when using TF32 or FP16 Tensor Cores. The exponent range is a critical issue when applying such techniques to quantum circuit simulation. While TF32 supports almost the same exponent range as FP32, FP16 supports a much smaller exponent range. In this work, we use the exponent range statistics of input tensor elements to select which Tensor Cores we use for the GEMM. We evaluate our method on Random Circuit Sampling (RCS), including Sycamore's quantum circuit, and show that the throughput is 1.86 times higher at maximum while maintaining accuracy.
Preprint
arXiv:2303.08989 [quant-ph]

Title
Automatic structural optimization of tree tensor networks
Reference
Physical Review Research 5, 013031 (2023)
DOI
10.1103/PhysRevResearch.5.013031
Author
Toshiya Hikihara, Hiroshi Ueda, Kouichi Okunishi, Kenji Harada, and Tomotoshi Nishino
Abstract
Tree tensor network (TTN) provides an essential theoretical framework for the practical simulation of quantum many-body systems, where the network structure defined by the connectivity of the isometry tensors plays a crucial role in improving its approximation accuracy. In this paper, we propose a TTN algorithm that enables us to automatically optimize the network structure by local reconnections of isometries to suppress the bipartite entanglement entropy on their legs. The algorithm can be seamlessly implemented to such a conventional TTN approach as density-matrix renormalization group. We apply the algorithm to the inhomogeneous antiferromagnetic Heisenberg spin chain having a hierarchical spatial distribution of the interactions. We then demonstrate that the entanglement structure embedded in the ground-state of the system can be efficiently visualized as a perfect binary tree in the optimized TTN. Possible improvements and applications of the algorithm are also discussed.
Comments
11 pages, 10 figures, 2 tables
Preprint
arXiv:2209.03196

Date
Aug 26, 2022
Conference (invited talk)
The 15th Asia Pacific Physics Conference (APPC15), Korea (online)
Title
Tensor renormalization group study of the non-equilibrium critical fixed point of the one-dimensional contact process
Abstract
The steady-state of many stochastic systems is non-equilibrium. We studied the phase of non-equilibrium systems and the transition similar to equilibrium systems. In particular, the critical phase transition is interesting because we can define the non-equilibrium universality class. To confirm the existence of a non-equilibrium critical fixed point, we study the time evolution operator of one-dimensional contact processes by using a tensor renormalization group technique. The time evolution operators converge to universal critical tensors in the tensor renormalization group flow. The spectrums of critical tensors are strongly anisotropic but share the intrinsic structure each for the universality class. The integer structure for the universality class of compact-directed percolation in the time direction is consistent with the exact spectrum structure of the diffusion-annihilation process.

We will hold the oneline workshop, Tensor Network States: Algorithms and Applications (TNSAA) 2021-2022, from Jan. 17 to Jan. 21, 2022. The series of workshop, TNSAA, has been organized for the purposes of exchanging new developments, having discussions toward future studies, and providing intruductory talks for new generation of researchers about tensor networks.

TOPICS

Toolkit of Bayesian Scaling Analysis

Reference application software of a new scaling analysis method of critical phenomena based on Bayesian inference.

To demo To details
Monte Carlo simulations

This demonstration shows a Monte Carlo simulation of the two-dimensional Ising model by three algorithms: Metropolis, Swendsen-Wang, and Wolff algorithms.

To demo

ABOUT

Kenji Harada

Kenji Harada ( 原田健自 )
Assistant Professor, Graduate School of Informatics, Kyoto University, Japan.
harada.kenji.8e@kyoto-u.ac.jp
Room 203, Research Bldg. No.8, Yoshida Campus, Kyoto Univ., Kyoto, 606-8501, Japan. Map (No.59)

LINKS