BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Synge Todo

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

ResearcherUnclaimed source record
cond-mat.stat-mechphysics.comp-phcond-mat.str-elcond-mat.dis-nnquant-phmath-phmath.MPmath.NAComputationcond-mat.softNumerical Analysisq-fin.CP

Catalog footprint

What is connected

24works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

24 published item(s)

preprint2022arXiv

Nested Iterative Shift-invert Diagonalization for Many-body Localization in the Random-field Heisenberg Chain

We study the many-body localization of the random-field Heisenberg chain using the nested shift-invert Lanczos method with an iterative linear solver. We use the minimum residual method (MINRES) inside each Lanczos iteration. The memory consumption of the proposed method is only proportional to the dimension of the Hilbert space. We also introduce a preconditioner that takes into account the effects of disorder and interaction in the random-field Heisenberg chain. As a probe of many-body localization transition, we propose a unitary operator called the twist operator, which has a clear interpretation in the real space. We discuss its behavior for thermal and localized eigenstates. We demonstrate the efficiency of the nested iterative shift-invert diagonalization method with the proposed preconditioner for the many-body localization problem and estimate the transition point of the random-field Heisenberg chain more precisely based on the finite-size analysis of the expectation value of the twist operator.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Optimized Implementation for Calculation and Fast-Update of Pfaffians Installed to the Open-Source Fermionic Variational Solver mVMC

In this article, we present a high performance, portable and well templated implementation for computing and fast-updating Pfaffian and inverse of an even-ranked skew-symmetric (antisymmetric) matrix. It is achieved with a skew-symmetric, blocked variant of the Parlett-Reid algorithm and a blocked update scheme based on the Woodbury matrix identity. Installation of this framework into the geminal-wavefunction-based many-variable Variational Monte Carlo (mVMC) code boosts sampling performance to up to more than $6$ times without changing Markov chain's behavior. The implementation is based on an extension of the BLAS-like instantiation software (BLIS) framework which has optimized kernel for many state-of-the-art processors including Intel Skylake-X, AMD EPYC Rome and Fujitsu A64FX.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Randomized-gauge test for machine learning of Ising model order parameter

Recently, machine learning has been applied successfully for identifying phases and phase transitions of the Ising models. The continuous phase transition is characterized by spontaneous symmetry breaking, which can not be detected in general from a single spin configuration. To investigate if neural networks can extract correlations among spin snapshots, we propose a new test using the random-gauge Ising model. We show that neural networks can extract the order parameter or the energy of the random-gauge model as in the ferromagnetic case. We also discuss how and where the information of random gauge is encoded in neural networks and attempt to reconstruct the gauge from the neural network parameters. We find that the fully connected network encodes the effect of random gauge to its weights naturally. In contrast, the convolutional network copes with the randomness by assigning different network parts to local gauge patterns. This observation indicates that although the latter demonstrates higher performance than the former for the present randomized-gauge test, the former is more effective and suitable for dealing with models with spatial randomness.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Bond-weighted Tensor Renormalization Group

We propose an improved tensor renormalization group (TRG) algorithm, the bond-weighted TRG (BTRG). In BTRG, we generalize the conventional TRG by introducing bond weights on the edges of the tensor network. We show that BTRG outperforms the conventional TRG and the higher-order tensor renormalization group with the same bond dimension, while its computation time is almost the same as that of TRG. Furthermore, BTRG can have non-trivial fixed-point tensors at an optimal hyperparameter. We demonstrate that the singular value spectrum obtained by BTRG is invariant under the renormalization procedure in the case of the two-dimensional Ising model at the critical point. This property indicates that BTRG performs the tensor contraction with high accuracy while keeping the scale-invariant structure of tensors.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

K$ω$ -- Open-source library for the shifted Krylov subspace method of the form $(zI-H)x=b$

We develop K$ω$, an open-source linear algebra library for the shifted Krylov subspace methods. The methods solve a set of shifted linear equations $(z_k I-H)x^{(k)}=b\, (k=0,1,2,...)$ for a given matrix $H$ and a vector $b$, simultaneously. The leading order of the operational cost is the same as that for a single equation. The shift invariance of the Krylov subspace is the mathematical foundation of the shifted Krylov subspace methods. Applications in materials science are presented to demonstrate the advantages of the algorithm over the standard Krylov subspace methods such as the Lanczos method. We introduce benchmark calculations of (i) an excited (optical) spectrum and (ii) intermediate eigenvalues by the contour integral on the complex plane. In combination with the quantum lattice solver $\mathcal{H} Φ$, K$ω$ can realize parallel computation of excitation spectra and intermediate eigenvalues for various quantum lattice models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multithreaded event-chain Monte Carlo with local times

We present a multithreaded event-chain Monte Carlo algorithm (ECMC) for hard spheres. Threads synchronize at infrequent breakpoints and otherwise scan for local horizon violations. Using a mapping onto absorbing Markov chains, we rigorously prove the correctness of a sequential-consistency implementation for small test suites. On x86 and ARM processors, a C++ (OpenMP) implementation that uses compare-and-swap primitives for data access achieves considerable speed-up with respect to single-threaded code. The generalized birthday problem suggests that for the number of threads scaling as the square root of the number of spheres, the horizon-violation probability remains small for a fixed simulation time. We provide C++ and Python open-source code that reproduces all our results.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Anisotropic Tensor Renormalization Group

We propose a new tensor renormalization group algorithm, Anisotropic Tensor Renormalization Group (ATRG), for lattice models in arbitrary dimensions. The proposed method shares the same versatility with the Higher-Order Tensor Renormalization Group (HOTRG) algorithm, i.e., it preserves the lattice topology after the renormalization. In comparison with HOTRG, both of the computation cost and the memory footprint of our method are drastically reduced, especially in higher dimensions, by renormalizing tensors in an anisotropic way after the singular value decomposition. We demonstrate the ability of ATRG for the square lattice and the simple cubic lattice Ising models. Although the accuracy of the present method degrades when compared with HOTRG of the same bond dimension, the accuracy with fixed computation time is improved greatly due to the drastic reduction of the computation cost.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Dynamical phase transition in Floquet optical bistable systems: An approach from finite-size quantum systems

We study a dynamical phase transition in optical bistable systems subject to a time-periodic driving field. The phase transition occurs in the structure of limit cycle as a function of the frequency of the driving field. In the thermodynamic limit, a single limit cycle is divided into two separated limit cycles at the transition point. In finite-size systems, however, there is always a single limit cycle due to the quantum tunneling effect. We use a Floquet dissipative map, which is a time-evolution operator over one period in a dynamics given by a quantum master equation, and discuss the decay rate of relaxation dynamics into the limit cycle based on the dominant eigenvalue of the map. We found that the decay rate exhibits qualitatively different system-size dependence before and after the phase transition, and it shows a finite-size scaling of spinodal phenomena around the transition point. The present work provides a systematic way of studying dynamical phase transition observed in time-periodically driven open systems in terms of the Floquet dissipative map.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Monte Carlo simulation with aspect ratio optimization: Anomalous anisotropic scaling in dimerized antiferromagnet

We present a method that optimizes the aspect ratio of a spatially anisotropic quantum lattice model during the quantum Monte Carlo simulation, and realizes the virtually isotropic lattice automatically. The anisotropy is removed by using the Robbins-Monro algorithm based on the correlation length in each direction. The method allows for comparing directly the value of critical amplitude among different anisotropic models, and identifying the universality more precisely. We apply our method to the staggered dimer antiferromagnetic Heisenberg model and demonstrate that the apparent non-universal behavior is attributed mainly to the strong size correction of the effective aspect ratio due to the existence of the cubic interaction.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Path-integral Monte Carlo for the Gauge-fixed Berry Connection and the Local Z_2 Berry Phase

We present a general formula of the gauge-fixed Berry connection which can be evaluated by path- integral Monte Carlo method. We also propose that the gauge-fixed local Berry connection can be an- other effective tool to estimate precisely the quantum critical point. For a demonstration, we calculate the gauge-fixed Berry connection and the local Z2 Berry phase of the antiferromagnetic Heisenberg model on a staggered bond-alternating ladder, and estimated quantum critical point is consistent with other methods.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Path-integral Monte Carlo method for the local Z_2 Berry phase

We present a loop cluster algorithm Monte Carlo method for calculating the local Z_2 Berry phase of the quantum spin models. The Berry connection, which is given as the inner product of two ground states with different local twist angles, is expressed as a Monte Carlo average on the worldlines with fixed spin configurations at the imaginary-time boundaries. The "complex weight problem" caused by the local twist is solved by adopting the meron cluster algorithm. We present the results of simulation on the antiferromagnetic Heisenberg model on an out-phase bond-alternating ladder to demonstrate that our method successfully detects the change in the valence bond pattern at the quantum phase transition point. We also propose that the gauge-fixed local Berry connection can be an effective tool to estimate pricisely the quantum critical point.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Critical temperature and correlation length of an elastic interaction model for spin-crossover materials

It has previously been pointed out that the coexistence of infinite-range and short-range interactions causes a system to have a phase transition of the mean-field universality class, in which the cluster size is finite even at the critical point. In the present paper, we study this property in a model of bistable molecules, whose size changes depending on the bistable states. The molecules can move in space, interacting via an elastic interaction. It is known that due to the different sizes, an effective long-range interaction between the spins appears, and thus this model has a mean-field type of phase transition. It is found that the scaling properties of the shift of the critical temperature from the pure short-range limit in the model with infinite-range and short-range interactions hold also in the present model, regarding the ratio of the size of the two states as a control parameter for the strength of the long-range interaction. By studying the structure factor, it is shown that the dependence of the cluster size at the critical temperature also shows the same scaling properties as a previously studied model with both infinite-range and short-range interactions. We therefore conclude that these scaling relations hold universally in hybrid models with both short-range and weak long-range interactions.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ordering and Excitation in Orbital Compass Model on a Checkerboard Lattice

We study an orbital compass model on a checkerboard lattice where orbital degree of freedom is represented by the pseudo-spin operator. Competition arises from an Ising interaction for the $z$ component of pseudo-spins along the vertical/horizontal bonds and an Ising interaction for the $x$ component along diagonal bonds. Classical and quantum compass models are analyzed by utilizing several analytical methods and numerical simulations. At a fully frustrated point where the two Ising interactions compete with each other, a macroscopic number of orbital configurations are degenerate in a classical ground state. This degeneracy is lifted by thermal and quantum fluctuations, and a staggered long-range order of the $z$ component of the pseudo-spin is realized. A tricritical point for this order appears due to competition between the bond dependent Ising interactions. Roles of geometrical frustration on excitation dynamics are also examined.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Orbital ordering in e_g orbital systems: Ground states and thermodynamics of the 120 degree model

Orbital degrees of freedom shape many of the properties of a wide class of Mott insulating, transition metal oxides with partially filled 3d-shells. Here we study orbital ordering transitions in systems where a single electron occupies the e_g orbital doublet and the spatially highly anisotropic orbital interactions can be captured by an orbital-only model, often called the 120 degree model. Our analysis of both the classical and quantum limits of this model in an extended parameter space shows that the 120 degree model is in close proximity to several T=0 phase transitions and various competing ordered phases. We characterize the orbital order of these nearby phases and their associated thermal phase transitions by extensive numerical simulations and perturbative arguments.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Bond-Dilution-Induced Quantum Phase Transitions in Heisenberg Antiferromagnets

Bond-dilution effects on the ground state of the square-lattice antiferromagnetic Heisenberg model, consisting of coupled bond-alternating chains, are investigated by means of the quantum Monte Carlo simulation. It is found that, when the ground state of the non-diluted system is a non-magnetic state with a finite spin gap, a sufficiently weak bond dilution induces a disordered state with a mid gap in the original spin gap, and under a further stronger bond dilution an antiferromagnetic long-range order emerges. While the site-dilution-induced long-range order is induced by an infinitesimal concentration of dilution, there exists a finite critical concentration in the case of bond dilution. We argue that this essential difference is due to the occurrence of two types of effective interactions between induced magnetic moments in the case of bond dilution, and that the antiferromagnetic long-range-ordered phase does not appear until the magnitudes of the two interactions become comparable.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Transfer-Matrix Study of Hard-Core Singularity for Hard-Square Lattice Gas

A singularity on the negative fugacity axis of the hard-square lattice gas is investigated in terms of numerical diagonalization of transfer matrices. The location of the singular point $z_c^-$ and the critical exponent $ν$ are accurately determined by the phenomenological renormalization technique as -0.11933888188(1) and 0.416667(1), respectively. It is also found that the central charge $c$ and the dominant scaling dimension $x_σ$ are -4.399996(8) and -0.3999996(7), respectively. These results strongly support that this singularity belongs to the same universality class as the Yang-Lee edge singularity ($c=-22/5, x_σ=-2/5 and ν=5/12$).

0 citations0 reviewsNo code linkedOpen access