Source author record

Kenji Harada

Kenji Harada appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2020arXiv

Finite-$m$ scaling analysis of Berezinskii-Kosterlitz-Thouless phase transitions and entanglement spectrum for the six-state clock model

We investigate the Berezinskii-Kosterlitz-Thouless transitions for the square-lattice six-state clock model with the corner-transfer matrix renormalization group (CTMRG). Scaling analyses for effective correlation length, magnetization, and entanglement entropy with respect to the cutoff dimension $m$ at the fixed point of CTMRG provide transition temperatures consistent with a variety of recent numerical studies. We also reveal that the fixed point spectrum of the corner transfer matrix in the critical intermediate phase of the six-state clock model is characterized by the scaling dimension consistent with the $c=1$ boundary conformal field theory associated with the effective $Z_6$ dual sine-Gordon model.

preprint2020arXiv

Universal spectrum structure at nonequilibrium critical points in the (1+1)-dimensional directed percolation

Using a tensor renormalization group method with oblique projectors for an anisotropic tensor network, we confirm that the rescaled spectrum of transfer matrices at nonequilibrium critical points in the (1+1)-dimensional directed percolation, a canonical model of nonequilibrium critical phenomena, is scale-invariant and its structure is universal.

preprint2015arXiv

Kernel method for corrections to scaling

Scaling analysis, in which one infers scaling exponents and a scaling function in a scaling law from given data, is a powerful tool for determining universal properties of critical phenomena in many fields of science. However, there are corrections to scaling in many cases, and then the inference problem becomes ill-posed by an uncontrollable irrelevant scaling variable. We propose a new kernel method based on Gaussian process regression to fix this problem generally. We test the performance of the new kernel method for some example cases. In all cases, when the precision of the example data increases, inference results of the new kernel method correctly converge. Because there is no limitation in the new kernel method for the scaling function even with corrections to scaling, unlike in the conventional method, the new kernel method can be widely applied to real data in critical phenomena.

preprint2015arXiv

SU(N) Heisenberg model with multi-column representations

The $\mathrm{SU}(N)$ symmetric antiferromagnetic Heisenberg model with multi-column representations on the two-dimensional square lattice is investigated by quantum Monte Carlo simulations. For the representation of Young diagram with two columns, we confirm that a valence-bond solid order appears as soon as the Néel order disappears at $N = 10$ indicating no intermediate phase. In the case of the representation with three columns, there is no evidence for both of the Néel and the valence-bond solid ordering for $N\ge 15$. This is actually consistent with the large-$N$ theory, which predicts that the VBS state immediately follows the Néel state, because the expected spontaneous order is too weak to be detected.

preprint2015arXiv

Thermal Phase Transition of Generalized Heisenberg Models for SU(N) Spins on Square and Honeycomb Lattices

We investigate thermal phase transitions to a valence-bond solid phase in SU(N) Heisenberg models with four- or six-body interactions on a square or honeycomb lattice, respectively. In both cases, a thermal phase transition occurs that is accompanied by rotational symmetry breaking of the lattice. We perform quantum Monte Carlo calculations in order to clarify the critical properties of the models. The estimated critical exponents indicate that the universality classes of the square- and honeycomb-lattice cases are identical to those of the classical $XY$ model with a $Z_4$ symmetry-breaking field and the 3-state Potts model, respectively. In the square-lattice case, the thermal exponent, $ν$, monotonically increases as the system approaches the quantum critical point, while the values of the critical exponents, $η$ and $γ/ν$, remain constant. From a finite-size scaling analysis, we find that the system exhibits weak universality, because the $Z_4$ symmetry-breaking field is always marginal. In contrast, $ν$ in the honeycomb-lattice case exhibits a constant value, even in the vicinity of the quantum critical point, because the $Z_3$ field remains relevant in the SU(3) and SU(4) cases.

preprint2014arXiv

Parallelized Quantum Monte Carlo Algorithm with Nonlocal Worm Updates

Based on the worm algorithm in the path-integral representation, we propose a general quantum Monte Carlo algorithm suitable for parallelizing on a distributed-memory computer by domain decomposition. Of particular importance is its application to large lattice systems of bosons and spins. A large number of worms are introduced and its population is controlled by a fictitious transverse field. For a benchmark, we study the size-dependence of the Bose-condensation order parameter of the hardcore Bose-Hubbard model with $L\times L\times βt = 10240\times 10240\times 16$, using 3200 computing cores, which shows good parallelization efficiency.

preprint2014arXiv

Symmetry-protected topological order and negative-sign problem for SO(N) bilinear-biquadratic chains

Using a generalized Jordan-Wigner transformation combined with the defining representation of the SO(N) spin, we map the SO(N) bilinear-biquadratic(BLBQ) spin chain into the N-color bosonic particle model. We find that, when the Jordan-Wigner transformation disentangles the symmetry-protected topological entanglement, this bosonic model becomes negative-sign free in the context of quantum Monte-Carlo simulation. For the SO(3) case, moreover, the Kennedy-Tasaki transformation for the S=1 BLBQ chain, which is also a topological disentangler, derives the same bosonic model through the dimer-R bases. We present the temperature dependence of the energy, entropy and string order parameter for the SO(N=3, 4, 5) BLBQ chains by a world-line Monte-Carlo simulation for the N-color bosonic particle model.

preprint2013arXiv

Possibility of Deconfined Criticality in SU(N) Heisenberg Models at Small N

To examine the validity of the scenario of the deconfined critical phenomena, we carry out a quantum Monte Carlo simulation for the SU($N$) generalization of the Heisenberg model with four-body and six-body interactions. The quantum phase transition between the SU($N$) Néel and valence-bond solid phases is characterized for $N=2,3,$ and $4$ on the square and honeycomb lattices. While finite-size scaling analysis works well up to the maximum lattice size ($L=256$) and indicates the continuous nature of the phase transition, a clear systematic change towards the first-order transition is observed in the estimates of the critical exponent $y \equiv 1/ν$ as the system size increases. We also confirm the relevance of a squared valence-bond solid field $Ψ^2$ for the SU(3) model.

preprint2012arXiv

Numerical study of incommensurability of the spiral state on spin-1/2 spatially anisotropic triangular antiferromagnets using entanglement renormalization

The ground state of an S=1/2 antiferromagnetic Heisenberg model on a spatially anisotropic triangular lattice, which is an effective model of Mott insulators on a triangular layer of organic charge transfer salts or Cs2CuCl4, is numerically studied. We apply a numerical variational method by using a tensor network with entanglement renormalization, which improves the capability of describing a quantum state. Magnetic ground states are identified for 0.7 <= J2/J1 <= 1 in the thermodynamic limit, where J1 and J2 denote the inner-chain and inter-chain coupling constants, respectively. Except for the isotropic case (J1=J2), the magnetic structure is spiral with an incommensurate wave vector that is different from the classical one. The quantum fluctuation weakens the effective coupling between chains, but the magnetic order remains in the thermodynamic limit. In addition, the incommensurate wave number is in good agreement with that of the series expansion method.

preprint2012arXiv

Study of the Shastry Sutherland Model Using Multi-scale Entanglement Renormalization Ansatz

We performed variational calculation based on the multi-scale entanglemnt renormalization ansatz, for the antiferromagnetic Heisenberg model on a Shastry Sutherland lattice (SSL). Our results show that at coupling ratio J'/J= 0.687(3), the system undergoes a quantum phase transition from the orthogonal dimer order to the plaquette valence bond solid phase, which then transits into the antiferromagnetic order above J'/J=0.75. In the presence of an external magnetic field, our calculations show clear evidences of various magnetic plateaux in systems with different coupling ratios range from 0.5 to 0.69. Our calculations are not limited to the small coupling ratio region, and we are able to show strong evidence of the presence of several supersolid phases, including ones above 1/2 and 1/3 plateaux. Such supersolid phases, which feature the coexistence of compressible superfluidity and crystalline long range order in triplet excitations, emerge at relatively large coupling ratio (J'/J>0.5). A schematic phase diagram of the SSL model in the presence of magnetic field is provided.

preprint2011arXiv

Bayesian Inference in the Scaling Analysis of Critical Phenomena

To determine the universality class of critical phenomena, we propose a method of statistical inference in the scaling analysis of critical phenomena. The method is based on Bayesian statistics, most specifically, the Gaussian process regression. It assumes only the smoothness of a scaling function, and it does not need a form. We demonstrate this method for the finite-size scaling analysis of the Ising models on square and triangular lattices. Near the critical point, the method is comparable in accuracy to the least-square method. In addition, it works well for data to which we cannot apply the least-square method with a polynomial of low degree. By comparing the data on triangular lattices with the scaling function inferred from the data on square lattices, we confirm the universality of the finite-size scaling function of the two-dimensional Ising model.

preprint2008arXiv

Diffusion in the Continuous-Imaginary-Time Quantum World-Line Monte Carlo Simulations with Extended Ensembles

The dynamics of samples in the continuous-imaginary-time quantum world-line Monte Carlo simulations with extended ensembles are investigated. In the case of a conventional flat ensemble on the one-dimensional quantum S=1 bi-quadratic model, the asymmetric behavior of Monte Carlo samples appears in the diffusion process in the space of the number of vertices. We prove that a local diffusivity is asymptotically proportional to the number of vertices, and we demonstrate the asymmetric behavior in the flat ensemble case. On the basis of the asymptotic form, we propose the weight of an optimal ensemble as $1/\sqrt{n}$, where $n$ denotes the number of vertices in a sample. It is shown that the asymmetric behavior completely vanishes in the case of the proposed ensemble on the one-dimensional quantum S=1 bi-quadratic model.

preprint2006arXiv

Dimer-Quadrupolar Quantum Phase Transition in the Quasi-One-Dimensional Heisenberg Model with Biquadratic Interaction

The quasi-one-dimensional S=1 Heisenberg antiferromagnet with a biquadratic term is investigated at zero temperature by quantum Monte Carlo simulation. As the magnitude of the inter-chain coupling is increased, the system undergoes a phase transition from a spontaneously dimerized phase to a Néel ordered or spin nematic phase. The numerical results suggest the possibility of an unconventional second-order transition in which the symmetry group characterizing one phase is not a subgroup of the other.

preprint2006arXiv

Quantum World-line Monte Carlo Method with Non-binary Loops and Its Application

A quantum world-line Monte Carlo method for high-symmetrical quantum models is proposed. Firstly, based on a representation of a partition function using the Matsubara formula, the principle of quantum world-line Monte Carlo methods is briefly outlined and a new algorithm using non-binary loops is given for quantum models with high symmetry as SU(N). The algorithm is called non-binary loop algorithm because of non-binary loop updatings. Secondary, one example of our numerical studies using the non-binary loop updating is shown. It is the problem of the ground state of two-dimensional SU(N) anti-ferromagnets. Our numerical study confirms that the ground state in the small N <= 4 case is a magnetic ordered Neel state, but the one in the large N >= 5 case has no magnetic order, and it becomes a dimer state.