Source author record

Hidemaro Suwa

Hidemaro Suwa 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

9works
10topics
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

9 published item(s)

preprint2022arXiv

Nematicity and fractional magnetization plateaus induced by spin-lattice coupling in the classical kagome-lattice Heisenberg antiferromagnet

We investigate the effect of spin-lattice coupling (SLC) on the magnetic properties of the classical kagome-lattice Heisenberg antiferromagnet (KHAF) using improved Monte Carlo updates. The lattice modes are represented by Einstein site phonons, which introduce effective further-neighbor spin interactions in addition to the nearest-neighbor biquadratic interactions. In the weak SLC, the macroscopically degenerate coplanar ground state remains at zero field, while a $\sqrt{3} \times \sqrt{3}$ ordered phase accompanied by a 1/3-magnetization plateau appears in external magnetic fields. In the strong SLC, we find a nematic order at zero field and a 1/9-magnetization plateau associated with a $3 \times 3$ collinear order. Near the phase transition between the 1/9- and 1/3-plateau states, the ergodicity in the single spin flip is practically broken, and slow dynamics appear. We propose that relevant KHAFs with strong SLC would be realized in spinel-based materials.

preprint2021arXiv

Geometric allocation approach to accelerating directed worm algorithm

The worm algorithm is a versatile technique in the Markov chain Monte Carlo method for both classical and quantum systems. The algorithm substantially alleviates critical slowing down and reduces the dynamic critical exponents of various classical systems. It is crucial to improve the algorithm and push the boundary of the Monte Carlo method for physical systems. We here propose a directed worm algorithm that significantly improves computational efficiency. We use the geometric allocation approach to optimize the worm scattering process: worm backscattering is averted, and forward scattering is favored. Our approach successfully enhances the diffusivity of the worm head (kink), which is evident in the probability distribution of the relative position of the two kinks. Performance improvement is demonstrated for the Ising model at the critical temperature by measurement of exponential autocorrelation times and asymptotic variances. The present worm update is approximately 25 times as efficient as the conventional worm update for the simple cubic lattice model. Surprisingly, our algorithm is even more efficient than the Wolff cluster algorithm, which is one of the best update algorithms. We estimate the dynamic critical exponent of the simple cubic lattice Ising model to be $z \approx 0.27$ in the worm update. The worm and the Wolff algorithms produce different exponents of the integrated autocorrelation time of the magnetic susceptibility estimator but the same exponent of the asymptotic variance. We also discuss how to quantify the computational efficiency of the Markov chain Monte Carlo method. Our approach can be applied to a wide range of physical systems, such as the $| ϕ|^4$ model, the Potts model, the O($n$) loop model, and lattice QCD.

preprint2016arXiv

Level spectroscopy in a two-dimensional quantum magnet: Linearly dispersing spinons at the deconfined quantum critical point

We study the level structure of excitations at the "deconfined" critical point separating antiferromagnetic and valence-bond-solid phases in two-dimensional quantum spin systems using the $J$-$Q$ model as an example. Energy gaps in different spin ($S$) and momentum (${\bf k}$) sectors are extracted from imaginary-time correlation functions obtained in quantum Monte Carlo simulations. We find strong quantitative evidence for deconfined linearly dispersing spinons with gapless points at ${\bf k}=(0,0)$, $(π,0)$, $(0,π)$, and $(π,π)$, as inferred from two-spinon excitations ($S=0$ and $S=1$ states) around these points. We also observe a duality between singlet and triplet excitations at the critical point and inside the ordered phases, in support of an enhanced symmetry, possibly SO(5).

preprint2015arXiv

Generalized Moment Method for Gap Estimation and Quantum Monte Carlo Level Spectroscopy

We formulate a convergent sequence for the energy gap estimation in the worldline quantum Monte Carlo method. The ambiguity left in the conventional gap calculation for quantum systems is eliminated. Our estimation will be unbiased in the low-temperature limit and also the error bar is reliably estimated. The level spectroscopy from quantum Monte Carlo data is developed as an application of the unbiased gap estimation. From the spectral analysis, we precisely determine the Kosterlitz-Thouless quantum phase-transition point of the spin-Peierls model. It is stablished that the quantum phonon with a finite frequency is essential to the critical theory governed by the antiadiabatic limit, i.e., the $k=1$ SU(2) Wess-Zumino-Witten model.

preprint2015arXiv

Stochastic approximation of dynamical exponent at quantum critical point

We have developed a unified finite-size scaling method for quantum phase transitions that requires no prior knowledge of the dynamical exponent $z$. During a quantum Monte Carlo simulation, the temperature is automatically tuned by the Robbins-Monro stochastic approximation method, being proportional to the lowest gap of the finite-size system. The dynamical exponent is estimated in a straightforward way from the system-size dependence of the temperature. As a demonstration of our novel method, the two-dimensional $S=1/2$ quantum $XY$ model in uniform and staggered magnetic fields is investigated in the combination of the world-line quantum Monte Carlo worm algorithm. In the absence of the uniform magnetic field, we obtain the fully consistent result with the Lorentz invariance at the quantum critical point, $z=1$, i.e., the three-dimensional classical $XY$ universality class. Under a finite uniform magnetic field, on the other hand, the dynamical exponent becomes two, and the mean-field universality with effective dimension $(2+2)$ governs the quantum phase transition.

preprint2015arXiv

Velocity of excitations in ordered, disordered and critical antiferromagnets

We test three different approaches, based on quantum Monte Carlo simulations, for computing the velocity $c$ of triplet excitations in antiferromagnets. We consider the standard $S=1/2$ one- and two-dimensional Heisenberg models, as well as a bilayer Heisenberg model at its critical point. Computing correlation functions in imaginary time and using their long-time behavior, we extract the lowest excitation energy versus momentum using improved fitting procedures and a generalized moment method. The velocity is then obtained from the dispersion relation. We also exploit winding numbers to define a cubic space-time geometry, where the velocity is obtained as the ratio of the spatial and temporal lengths of the system when all winding number fluctuations are equal. The two methods give consistent results for both ordered and critical systems, but the winding-number estimator is more precise. For the Heisenberg chain, we accurately reproduce the exactly known velocity. For the two-dimensional Heisenberg model, our results are consistent with other recent calculations, but with an improved statistical precision; $c=1.65847(4)$. We also use the hydrodynamic relation $c^2=ρ_s/χ_\perp(q\to 0)$ between $c$, the spin stiffness $ρ_s$, and the transversal susceptibility $χ_\perp$, using the smallest non-zero momentum $q=2π/L$. This method also is well controlled in two dimensions, but the cubic criterion for winding numbers delivers better numerical precision. In one dimension the hydrodynamic relation is affected by logarithmic corrections which make accurate extrapolations difficult. As an application of the winding-number method, for the quantum-critical bilayer model our high-precision determination of the velocity enables us to quantitatively test, at an unprecedented level, field-theoretic predictions for low-temperature scaling forms where $c$ enters.

preprint2013arXiv

Geometric Allocation Approaches in Markov Chain Monte Carlo

The Markov chain Monte Carlo method is a versatile tool in statistical physics to evaluate multi-dimensional integrals numerically. For the method to work effectively, we must consider the following key issues: the choice of ensemble, the selection of candidate states, the optimization of transition kernel, algorithm for choosing a configuration according to the transition probabilities. We show that the unconventional approaches based on the geometric allocation of probabilities or weights can improve the dynamics and scaling of the Monte Carlo simulation in several aspects. Particularly, the approach using the irreversible kernel can reduce or sometimes completely eliminate the rejection of trial move in the Markov chain. We also discuss how the space-time interchange technique together with Walker's method of aliases can reduce the computational time especially for the case where the number of candidates is large, such as models with long-range interactions.

preprint2012arXiv

Geometric Allocation Approach for Transition Kernel of Markov Chain

We introduce a new geometric approach that constructs a transition kernel of Markov chain. Our method always minimizes the average rejection rate and even reduce it to zero in many relevant cases, which cannot be achieved by conventional methods, such as the Metropolis-Hastings algorithm or the heat bath algorithm (Gibbs sampler). Moreover, the geometric approach makes it possible to find not only a reversible but also an irreversible solution of rejection-free transition probabilities. This is the first versatile method that can construct an irreversible transition kernel in general cases. We demonstrate that the autocorrelation time (asymptotic variance) of the Potts model becomes more than 6 times as short as that by the conventional Metropolis-Hastings algorithm. Our algorithms are applicable to almost all kinds of Markov chain Monte Carlo methods and will improve the efficiency.

preprint2010arXiv

Markov Chain Monte Carlo Method without Detailed Balance

We present a specific algorithm that generally satisfies the balance condition without imposing the detailed balance in the Markov chain Monte Carlo. In our algorithm, the average rejection rate is minimized, and even reduced to zero in many relevant cases. The absence of the detailed balance also introduces a net stochastic flow in a configuration space, which further boosts up the convergence. We demonstrate that the autocorrelation time of the Potts model becomes more than 6 times shorter than that by the conventional Metropolis algorithm. Based on the same concept, a bounce-free worm algorithm for generic quantum spin models is formulated as well.