BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Naoki Kawashima

Naoki Kawashima contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.str-elcond-mat.stat-mechphysics.comp-phcond-mat.quant-gasmath.NAcond-mat.dis-nncond-mat.mtrl-scicond-mat.supr-conhep-lathep-thmath-phmath.MPNumerical Analysis

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
13topics
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

20 published item(s)

preprint2025arXiv

Superfluid dome in the spatially modulated two-dimensional XY model

In strongly correlated electron systems, superconductivity and charge density waves often coexist in close proximity, suggesting a deeper relationship between these competing phases. Recent research indicates that these orders can intertwine, with the superconducting order parameter coupling to modulations in the electronic density. To elucidate this interplay, we study a two-dimensional XY model with a periodic modulation of the coupling strength in one spatial direction. Using a combination of tensor network methods and Monte Carlo simulations, we reveal a non-monotonic, dome-like dependence of $T_c$ on the modulation wavelength, with the peak $T_c$ shifting to longer wavelengths as the modulation strength grows. The origin of this phenomenon is traced back to an effective pinning of vortices in the valleys of the modulation, confirmed by a comparison to modulated $q$-state clock models. These findings shed new light on the phase behavior of intertwined superconducting and charge-ordered states, offering a deeper understanding of their complex interactions.

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

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Multiple magnetization plateaus induced by farther neighbor interaction in an $S = 1$ two-leg Heisenberg spin ladder

We study the magnetization process of the $S=1$ Heisenberg model on a two-leg ladder with farther neighbor spin-exchange interaction. We consider the interaction that couples up to the next-nearest neighbor rungs and find an exactly solvable regime where the ground states become product states. The next-nearest neighbor interaction tends to stabilize magnetization plateaus at multiples of 1/6. In most of the exactly solvable regime, a single magnetization curve shows two series of plateaus with different periodicities.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Universal and Non-Universal Correction Terms of Bose Gases in Dilute Region: a Quantum Monte Carlo Study

We study dilute gases of interacting bosons at zero-temperature in the region where the system is characterized only by the s-wave scattering length. We carry out quantum Monte Carlo simulation of the Bose Hubbard model and a continuous-space hard-core model. Fitting the extended Lee-Huang-Yang formula to the Monte Carlo results establishes the detailed mapping from the lattice model to the continuous field theory characterized only by the s-wave scattering length. Our estimate of the intrinsic s-wave scattering length $a_s$ of the Bose-Hubbard model is $a_s/a_l = 0.316(2)$ where $a_l$ is the lattice constant. It turned out that inclusion of the universal second correction term of $O(n \log n)$ makes the fitting worse while it is restored by inclusion of the non-universal third correction term, of which the existence was predicted analytically.

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Boundary conformal spectrum and surface critical behaviors of the classical spin systems: a tensor network renormalization study

We numerically obtain the conformal spectrum of several classical spin models on a two-dimensional lattice with open boundaries, for every boundary fixed point obtained by the Cardy's derivation [J. L. Cardy, Nucl. Phys. B 324, 581 (1989)]. In order to extract accurate conformal data, we implement the tensor network renormalization algorithm [G. Evenbly and G. Vidal, Phys. Rev. Lett. 115, 180405 (2015)] extended so as to be applicable to a square lattice with open boundaries. We successfully compute the boundary conformal spectrum consistent with the underlying boundary conformal field theories (BCFTs) for the Ising, tri-critical Ising, and 3-state Potts models on the lattice, which allows us to confirm the validity of the BCFT analyses for the surface critical behaviors of those lattice models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Construction of Variational Matrix Product States for the Heisenberg Spin-1 Chain

We propose a simple variational wave function that captures the correct ground state energy of the spin-1 Heisenberg chain model to within 0.04\%. The wave function is written in the matrix product state (MPS) form with the bond dimension $D=8$, and characterized by three fugacity parameters. The proposed MPS generalizes the Affleck-Kennedy-Lieb-Tasaki (AKLT) state by dressing it with dimers, trimers, and general $q$-dimers. The fugacity parameters control the number and the average size of the $q$-mers. Furthermore, the $D=8$ variational MPS state captures the ground states of the entire family of bilinear-biquadratic Hamiltonian belonging to the Haldane phase to high accuracy. The 2-4-2 degeneracy structure in the entanglement spectrum of our MPS state is found to match well with the results of density matrix renormalization group (DMRG) calculation, which is computationally much heavier. Spin-spin correlation functions also find excellent fit with those obtained by DMRG.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Continuous phase transition between Neel and valence bond solid phases in a J-Q-like spin ladder system

We investigate a quantum phase transition between a Neel phase and a valence bond solid (VBS) phase, in each of which a different Z2 symmetry is broken, in a spin-1/2 two-leg XXZ ladder with a four-spin interaction. The model can be viewed as a one-dimensional variant of the celebrated J-Q model on a square lattice. By means of variational uniform matrix product state calculations and an effective field theory, we determine the phase diagram of the model, and present evidences that the Neel-VBS transition is continuous and belongs to the Gaussian universality class with the central charge c=1. In particular, the critical exponents $β, η,$ and, $ν$ are found to satisfy the constraints expected for a Gaussian transition within numerical accuracy. These exponents do not detectably change along the phase boundary while they are in general allowed to do so for the Gaussian class.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Data Assimilation Method for Experimental and First-Principles Data: Finite-Temperature Magnetization of (Nd,Pr,La,Ce)$_{2}$(Fe,Co,Ni)$_{14}$B

We propose a data-assimilation method for evaluating the finite-temperature magnetization of a permanent magnet over a high-dimensional composition space. Based on a general framework for constructing a predictor from two data sets including missing values, a practical scheme for magnetic materials is formulated in which a small number of experimental data in limited composition space are integrated with a larger number of first-principles calculation data. We apply the scheme to (Nd$_{1-α-β-γ}$Pr$_α$La$_β$Ce$_γ$)$_{2}$(Fe$_{1-δ-ζ}$Co$_δ$Ni$_ζ$)$_{14}$B. The magnetization in the whole $(α, β, γ, δ, ζ)$ space at arbitrary temperature is obtained. It is shown that the Co doping does not enhance the magnetization at low temperatures, whereas the magnetization increases with increasing $δ$ above 320 K.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Global optimization of tensor renormalization group using the corner transfer matrix

A tensor network renormalization algorithm with global optimization based on the corner transfer matrix is proposed. Since the environment is updated by the corner transfer matrix renormalization group method, the forward-backward iteration is unnecessary, which is a time-consuming part of other methods with global optimization. In addition, a further approximation reducing the order of the computational cost of contraction for the calculation of the coarse-grained tensor is proposed. The computational time of our algorithm in two dimensions scales as the sixth power of the bond dimension while the higher-order tensor renormalization group and the higher-order second renormalization group methods have the seventh power. We perform benchmark calculations in the Ising model on the square lattice and show that the time-to-solution of the proposed algorithm is faster than that of other methods.

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Tensor network wavefunction of $S=1$ Kitaev spin liquids

Recently there has been a great interest in understanding quantum spin liquid phases with varying spin magnitude, partly due to possible material realizations. A number of recent numerical computations suggest that the ground state of the S=1 Kitaev model may be a quantum spin liquid in analogy to the renowned $S$=$1/2$ model. On the other hand, the nature of the ground state remains elusive since the $S$=$1$ model is not exactly solvable in contrast to the $S$=$1/2$ model. In this work, we construct a tensor network ground state wavefunction for the S=1 Kitaev model, which is explicitly written in terms of physical spin operators. We explain how this class of wavefunctions can be successfully used for variational computations and compare the outcomes to known results on finite size systems. We establish the existence of distinct topological sectors on torus by constructing the minimally entangled states in the degenerate ground state manifold and evaluating topological entanglement entropy. Our results suggest that the ground state of the S=1 Kitaev model is a gapped quantum spin liquid with Z2 gauge structure and Abelian quasiparticles.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Calculation of higher-order moments by higher-order tensor renormalization group

A calculation method for higher-order moments of physical quantities, including magnetization and energy, based on the higher-order tensor renormalization group is proposed. The physical observables are represented by impurity tensors. A systematic summation scheme provides coarse-grained tensors including multiple impurities. Our method is compared with the Monte Carlo method on the two-dimensional Potts model. While the nature of the transition of the $q$-state Potts model has been known for a long time owing to the analytical arguments, a clear numerical confirmation has been difficult due to extremely long correlation length in the weakly first-order transitions, e.g., for $q=5$. A jump of the Binder ratio precisely determines the transition temperature. The finite-size scaling analysis provides critical exponents and distinguishes the weakly first-order and the continuous transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Tensor-Ring Decomposition with Index-Splitting

Tensor-ring decomposition of tensors plays a key role in various applications of tensor network representation in physics as well as in other fields. In most heuristic algorithms for the tensor-ring decomposition, one encounters the problem of local-minima trapping. Particularly, the minima related to the topological structure in the correlation are hard to escape. Therefore, identification of the correlation structure, somewhat analogous to finding matching ends of entangled strings, is the task of central importance. We show how this problem naturally arises in physical applications, and present a strategy for winning this string-pull game.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Crossover Behavior from Decoupled Criticality

We study the thermodynamic phase transition of a spin Hamiltonian comprising two 3D magnetic sublattices. Each sublattice contains XY spins coupled by the usual bilinear exchange, while spins in different sublattices only interact via biquadratic exchange. This Hamiltonian is an effective model for XY magnets on certain frustrated lattices such as body centered tetragonal. By performing a cluster Monte Carlo simulation, we investigate the crossover from the 3D-XY fixed point (decoupled sublattices) and find a systematic flow toward a first-order transition without a separatrix or a new fixed point. This strongly suggests that the correct asymptotic behavior is a first-order transition.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Phase Transition in Potts Model with Invisible States

We study phase transition in the ferromagnetic Potts model with invisible states that are added as redundant states by mean-field calculation and Monte Carlo simulation. Invisible states affect the entropy and the free energy, although they do not contribute to the internal energy. The internal energy and the number of degenerated ground states do not change, if invisible states are introduced into the standard Potts model. A second-order phase transition takes place at finite temperature in the standard $q$-state ferromagnetic Potts model on two-dimensional lattice for $q=2,3$, and 4. However, our present model on two-dimensional lattice undergoes a first-order phase transition with spontaneous $q$-fold symmetry breaking ($q=2,3$, and 4) due to entropy effect of invisible states. We believe that our present model is a fundamental model for analysis of a first-order phase transition with spontaneous discrete symmetry breaking.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Reponse to Comment on "Direct mapping of the finite temperature phase diagram of strongly correlated quantum models"

Pollet et al have questioned our claim staked in Phys. Rev. Lett. 103, 085701 (2009) that from the location of sharp features in the compressibility in the trap, it is possible to map out the phase boundary between two phases within reasonable accuracy. We stand by it and use this reply to sharpen the conditions under which "sharp" features in compressibility reflect the effects of critical fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Finite-Size Scaling for Quantum Criticality above the Upper Critical Dimension: Superfluid-Mott-Insulator Transition in Three Dimensions

Validity of modified finite-size scaling above the upper critical dimension is demonstrated for the quantum phase transition whose dynamical critical exponent is $z=2$. We consider the $N$-component Bose-Hubbard model, which is exactly solvable and exhibits mean-field type critical phenomena in the large-$N$ limit. The modified finite-size scaling holds exactly in that limit. However, the usual procedure, taking the large system-size limit with fixed temperature, does not lead to the expected (and correct) mean-field critical behavior due to the limited range of applicability of the finite-size scaling form. By quantum Monte Carlo simulation, it is shown that the same holds in the case of N=1.

0 citations0 reviewsNo code linkedOpen access