Source author record

Tsuyoshi Okubo

Tsuyoshi Okubo 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

20works
7topics
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

20 published item(s)

preprint2026arXiv

Clifford Circuits Augmented Grassmann Matrix Product States

Recent advances in combining Clifford circuits with tensor-network (TN) methods have shown that classically simulable disentanglers can suppress substantial portions of the entanglement structure, effectively alleviating the bond-dimension bottleneck in TN simulations. In this work, we develop a variational TN framework based on Grassmann tensor networks, which natively encode fermionic statistics while preserving locality. By incorporating locally defined Clifford circuits within the fermionic formalism, we simulate benchmark models including the tight-binding and $t$-$V$ models. Our results show that Clifford disentangling removes the classically simulable component of entanglement, leading to a reduced bond dimension and improved accuracy in ground-state energy estimates. Interestingly, once the natural Grassmann-evenness requirement of the fermionic formulation is taken into account and Clifford gates with identical entanglement action are grouped together, the original set of 11520 two-qubit Clifford gates reduces to only 12 distinct gates. This strong reduction leads to a more efficient disentangling scheme within the fermionic framework. These findings highlight the potential of Clifford-augmented Grassmann TNs as a scalable and accurate tool for studying strongly correlated fermionic systems, particularly in higher dimensions.

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.

preprint2022arXiv

TeNeS: Tensor Network Solver for Quantum Lattice Systems

TeNeS (Tensor Network Solver) is a free/libre open-source software program package for calculating two-dimensional many-body quantum states based on the tensor network method and the corner transfer matrix renormalization group (CTMRG) method. This package calculates ground-state wavefunctions for user-defined Hamiltonians and evaluates user-defined physical quantities such as magnetization and correlation functions. For certain predefined models and lattices, there is a tool that makes it easy to generate input files. TeNeS uses an OpenMP/MPI hybrid parallelized tensor operation library and thus can perform large-scale calculations using massively parallel machines.

preprint2020arXiv

Abelian and non-Abelian chiral spin liquids in a compact tensor network representation

We provide new insights into the Abelian and non-Abelian chiral Kitaev spin liquids on the star lattice using the recently proposed loop gas (LG) and string gas (SG) states [H.-Y. Lee, R. Kaneko, T. Okubo, N. Kawashima, Phys. Rev. Lett. 123, 087203 (2019)]. Those are compactly represented in the language of tensor network. By optimizing only one or two variational parameters, accurate ansatze are found in the whole phase diagram of the Kitaev model on the star lattice. In particular, the variational energy of the LG state becomes exact(within machine precision) at two limits in the model, and the criticality at one of those is analytically derived from the LG feature. It reveals that the Abelian CSLs are well demonstrated by the short-ranged LG while the non-Abelian CSLs are adiabatically connected to the critical LG where the macroscopic loops appear. Furthermore, by constructing the minimally entangled states and exploiting their entanglement spectrum and entropy, we identify the nature of anyons and the chiral edge modes in the non-Abelian phase with the Ising conformal field theory.

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.

preprint2020arXiv

Magnetic field induced quantum phases in a tensor network study of Kitaev magnets

Recent discovery of the half quantized thermal Hall conductivity in $α$-RuCl$_3$, a candidate material for the Kitaev spin liquid, suggests the presence of a highly entangled quantum state in external magnetic fields. This field induced phase appears between the low field zig-zag magnetic order and the high field polarized state. Motivated by this experiment, we study possible field induced quantum phases in theoretical models of the Kitaev magnets, using the two dimensional tensor network approach or infinite tensor product states. We find various quantum ground states in addition to the chiral Kitaev spin liquid occupying a small area in the phase diagram. They form a band of emergent quantum phases in an intermediate window of external magnetic fields, somewhat reminiscent of the experiment. We discuss the implications of these results in view of the experiment and previous theoretical studies.

preprint2020arXiv

Zero-energy excitation in the classical kagome antiferromagnet NaBa$_{2}$Mn$_{3}$F$_{11}$

We performed inelastic neutron scattering measurements on a polycrystalline sample of a classical kagome antiferromagnet NaBa$_{2}$Mn$_{3}$F$_{11}$ to investigate the possibility of a dispersionless zero-energy excitation associated with rotation of spins along the chains. The observed spectra indeed exhibit such an excitation with strong intensity at low energy, as well as dispersive excitations with weak intensity at high energy. Combining the measurements with calculations from linear spin-wave theory reveals that NaBa$_{2}$Mn$_{3}$F$_{11}$ is a good realization of the classical kagome antiferromagnet which exhibits a dispersionless mode lifted by the magnetic dipole-dipole interaction.

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.

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.

preprint2014arXiv

A Scaling Relation for Dangerously Irrelevant Symmetry-Breaking Fields

We propose a scaling relation for critical phenomena in which a symmetry-breaking field is dengerously irrelevant. We confirm its validity on the 6-state clock model in three and four dimensions by numerical simulation. In doing so, we point out the problem in the previously-used order parameter, and present an alternative evidence based on the mass-dependent fluctuation.

preprint2014arXiv

Kagome-Triangular Lattice Antiferromagnet NaBa2Mn3F11

A new type of magnetically frustrated lattice is found in the layered fluoride NaBa2Mn3F11. A kagome-type array of regular triangles composed of Mn2+ ions (spin 5/2) deforms so as to generate the next-nearest-neighbor interaction J2 between three out of six spins in the hexagon of a normal kagome lattice, in addition to the nearest-neighbor interaction J1. As a function of |J2/J1|, this lattice can interconnect the kagome (J2 = 0) and the triangular (J2 = J1) lattices and thus is called the kagome-triangular (KT) lattice. Magnetic susceptibility and heat capacity measurements performed on a polycrystalline sample of NaBa2Mn3F11 show an intensive short-range antiferromagnetic correlation below 14 K probably due to the specific magnetic frustration of the KT lattice. In addition, a long-range order at 2.0 K is observed, which is significantly low compared with the antiferromagnetic Weiss temperature of 32.3 K. Theoretical considerations of Heisenberg spins in the KT lattice reveal unique non-coplanar magnetic orders in the case of ferromagnetic J1 and antiferromagnetic J2. NaBa2Mn3F11 may actually exist in this regime according to the results of our analysis based on classical Monte-Carlo simulation.

preprint2014arXiv

Snapshot spectrum and critical phenomenon for two-dimensional classical spin systems

We investigate the eigenvalue distribution of the snapshot density matrix (SDM) generated by Monte Carlo simulation for two-dimensional classical spin systems. We find that the distribution in the high-temperature limit is well explained by the random-matrix theory, while that in the low-temperature limit can be characterized by the zero-eigenvalue condensation. At the critical point, we obtain the power-law distribution with a nontrivial exponent $α\equiv(2-η)/(1-η)$ and the asymptotic form of the snapshot entropy, on the basis of the relationship of the SDM with the correlation function matrix. The aspect-ratio dependence of the SDM spectrum is also mentioned.

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

Finite-temperature transition of the distorted kagome-lattice Heisenberg antiferromagnet

Motivated by the recent experiment on kagome-lattice antiferromagnets, we study the zero-field ordering behavior of the antiferromagnetic classical Heisenberg model on a uniaxially distorted kagome lattice by Monte Carlo simulations. A first-order transition, which has no counterpart in the corresponding undistorted model, takes place at a very low temperature. Origin of the transition is ascribed to a cooperative proliferation of topological excitations inherent to the model.

preprint2011arXiv

Cubic and non-cubic multiple-q states in the Heisenberg antiferromagnet on the pyrochlore lattice

The ordering of the classical Heisenberg model on the pyrochlore lattice with the antiferromagnetic nearest-neighbor interaction J1 and the ferromagnetic next-nearest-neighbour interaction J2 is investigated by means of a mean-field analysis and a Monte Carlo simulation. For a moderate J2/J1-value, the model exhibits a first-order transition into an incommensurate multiple-q ordered state where multiple Bragg peaks coexist in the spin structure factor. We show that there are two types of metastable multiple-q states, a cubic symmetric sextuple-q state and a non-cubic symmetric quadruple-q state. Based on a Monte Carlo simulation, we find that the cubic sextuple-q state appears just below the first-order transition temperature, while another transition from the cubic sextuple-q state to the non-cubic quadruple-q state occurs at a lower temperature.

preprint2011arXiv

Multiple-q states and skyrmion lattice of the triangular-lattice Heisenberg antiferromagnet under magnetic fields

Ordering of the frustrated classical Heisenberg model on the triangular-lattice with an incommensurate spiral spin structure is studied under magnetic fields by means of a mean-field analysis and a Monte Carlo simulation. Several types of multiple-q states including the "skyrmion-lattice" state is observed in addition to the standard single-q state. In contrast to the Dzyaloshinskii-Moriya interaction driven system, the present model allows both skyrmions and anti-skyrmions, together with a new thermodynamic phase where skyrmion and anti-skyrmion lattices form a domain state.

preprint2010arXiv

Novel spin-liquid states in the frustrated Heisenberg antiferromagnet on the honeycomb lattice

Recent experiment on a honeycomb-lattice Heisenberg antiferromagnet (AF) Bi$_3$Mn$_4$O$_{12}$(NO$_3$) revealed a novel spin-liquid-like behavior down to low temperature, which was ascribed to the frustration effect due to the competition between the AF nearest- and next-nearest-neighbor interactions $J_1$ and $J_2$. Motivated by the experiment, we study the ordering of the $J_1$ -$J_2$ frustrated classical Heisenberg AF on a honeycomb lattice both by a low-temperature expansion and a Monte Carlo simulation. The model has been known to possess a massive degeneracy of the ground state, which, however, might be lifted due to thermal fluctuations leading to a unique ordered state, the effect known as 'order-by-disorder'. We find that the model exhibits an intriguing ordering behavior, particularly near the AF phase boundary. The energy scale of the order-by-disorder is suppressed there down to extremely low temperatures, giving rise to exotic spin-liquid states like a "ring-liquid" or a "pancake-liquid" state accompanied by the characteristic spin structure factor and the field-induced antiferromagnetism. We argue that the recent experimental data are explicable if the system is in such exotic spin-liquid states.

preprint2010arXiv

Phase transition of the three-dimensional chiral Ginzburg-Landau model -- search for the chiral phase

Nature of the phase transition of regularly frustrated vector spin systems in three dimensions is investigated based on a Ginzburg-Landau-type effective Hamiltonian. On the basis of the variational analysis of this model, Onoda et al recently suggested the possible occurrence of a chiral phase, where the vector chirality exhibits a long-range order without the long-range order of the spin [Phys. Rev. Lett. 99, 027206 (2007)]. In the present paper, we elaborate their analysis by considering the possibility of a first-order transition which was not taken into account in their analysis. We find that the first-order transition indeed occurs within the variational approximation, which significantly reduces the stability range of the chiral phase, while the chiral phase still persists in a restricted parameter range. Then, we perform an extensive Monte Carlo simulation focusing on such a parameter range. Contrary to the variational result, however, we do not find any evidence of the chiral phase. The range of the chiral phase, if any, is estimated to be less than 0.1% in the temperature width.

preprint2010arXiv

Signature of a Z_2 vortex in the dynamical correlations of the triangular-lattice Heisenberg antiferromagnet

Dynamical properties of the classical Heisenberg antiferromagnet on the triangular lattice are investigated by means of a spin-dynamics simulation and an analytical calculation. While the model was suggested to exhibit a topological transition driven by the Z_2-vortex binding-unbinding, weakness of the associated thermodynamic singularity has made it difficult to observe the evidence of the Z_2 vortex experimentally so far, only some indirect support. Here, we show that the signature of the Z_2-vortex excitation and the vortex-induced topological transition can be captured from the dynamical spin correlations. In particular, the dynamical spin structure factor exhibits a characteristic central peak around the Z_2-vortex transition, and this central peak will be a fingerprint of the Z_2-vortex excitation.

preprint2010arXiv

Z_2-vortex ordering of the triangular-lattice Heisenberg antiferromagnet

Ordering of the classical Heisenberg antiferromagnet on the triangular lattice is studied by means of a mean-field calculation, a scaling argument and a Monte Carlo simulation, with special attention to its vortex degree of freedom. The model exhibits a thermodynamic transition driven by the Z_2-vortex binding-unbinding, at which various thermodynamic quantities exhibit an essential singularity. The low-temperature state is a "spin-gel" state with a long but finite spin correlation length where the ergodicity is broken topologically. Implications to recent experiments on triangular-lattice Heisenberg antiferromagnets are discussed.