Researcher profile

Tsuyoshi Okubo

Tsuyoshi Okubo contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 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.

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.