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

Naoki Kawashima

Naoki Kawashima 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-mechcond-mat.str-elcond-mat.quant-gascond-mat.mtrl-sciphysics.comp-phcond-mat.dis-nncond-mat.supr-conmath.NAquant-phhep-lathep-thmath-phmath.MPNumerical Analysis

Catalog footprint

What is connected

46works
14topics
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

46 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
preprint2016arXiv

Clues and criteria for designing Kitaev spin liquid revealed by thermal and spin excitations of honeycomb iridates Na$_2$IrO$_3$

Contrary to the original expectation, Na$_2$IrO$_3$ is not a Kitaev's quantum spin liquid (QSL) but shows a zig-zag-type antiferromagnetic order in experiments. Here we propose experimental clues and criteria to measure how a material in hand is close to the Kitaev's QSL state. For this purpose, we systematically study thermal and spin excitations of a generalized Kitaev-Heisenberg model studied by Chaloupka $et$ $al$. in Phys. Rev. Lett. 110, 097204 (2013) and an effective ab initio Hamiltonian for Na$_2$IrO$_3$ proposed by Yamaji $et$ $al$. in Phys. Rev. Lett. 113, 107201 (2014), by employing a numerical diagonalization method. We reveal that closeness to the Kitaev's QSL is characterized by the following properties, besides trivial criteria such as reduction of magnetic ordered moments and Neel temperatures: (1) Two peaks in the temperature dependence of specific heat at $T_{\ell}$ and $T_h$ caused by the fractionalization of spin to two types of Majorana fermions. (2) In between the double peak, prominent plateau or shoulder pinned at $(R/2)\ln 2$ in the temperature dependence of entropy, where $R$ is the gas constant. (3) Failure of the linear spin wave approximation at the low-lying excitations of dynamical structure factors. (4) Small ratio $T_{\ell}/T_h$ close to or less than 0.03. According to the proposed criteria, Na$_2$IrO$_3$ is categorized to a compound close to the Kitaev's QSL, and is proven to be a promising candidate for the realization of the QSL if the relevant material parameters can further be tuned by making thin film of Na$_2$IrO$_3$ on various substrates or applying axial pressure perpendicular to the honeycomb networks of iridium ions. Applications of these characterization to (Na$_{1-x}$Li$_x$)$_2$IrO$_3$ and other related materials are also discussed.

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

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

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.

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

Phase Transitions with Discrete Symmetry Breaking in Antiferromagnetic Heisenberg Models on a Triangular Lattice

We study phase transition behavior of the Heisenberg model on a distorted triangular lattice with competing interactions. The ground-state phase diagram indicates that underlying symmetry can be changed by tuning parameters. We focus on two cases in which a phase transition with discrete symmetry breaking occurs. The first is that the order parameter space is SO(3)$\times C_3$. In this case, a first-order phase transition, with threefold symmetry breaking, occurs. The second has the order parameter space SO(3)$\times Z_2$. In this case, a second-order phase transition occurs with twofold symmetry breaking. To investigate finite-temperature properties of these phase transitions from a microscopic viewpoint, we introduce a method to make the connection between continuous frustrated spin systems and the Potts model with invisible states.

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
preprint2013arXiv

Second-order phase transition in the Heisenberg model on a triangular lattice with competing interactions

We discover an example where the dissociation of the Z2 vortices occurs at the second-order phase transition point. We investigate the nature of phase transition in a classical Heisenberg model on a distorted triangular lattice with competing interactions. The order parameter space of the model is SO(3)xZ2. The dissociation of the Z2 vortices which comes from SO(3) and a second-order phase transition with Z2 symmetry breaking occur at the same temperature. We also find that the second-order phase transition belongs to the universality class of the two-dimensional Ising model.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Method to Control Order of Phase Transition: Invisible States in Discrete Spin Models

It is an important topic to investigate nature of the phase transition in wide area of science such as statistical physics, materials science, and computational science. Recently it has been reported the efficiency of quantum adiabatic evolution/quantum annealing for systems which exhibit a phase transition, and we cannot obtain a good solution in such systems. Thus, to control the nature of the phase transition has been also attracted attention in quantum information science. In this paper we review nature of the phase transition and how to control the order of the phase transition. We take the Ising model, the standard Potts model, the Blume-Capel model, the Wajnflasz-Pick model, and the Potts model with invisible states for instance. Until now there is no general method to avoid the difficulty of annealing method in systems which exhibit a phase transition. It is a challenging problem to propose a method how to erase or how to control the nature of the phase transition in the target system.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Commensurate Supersolid of Three-Dimensional Lattice Bosons

Using an unbiased quantum Monte Carlo method, we obtain convincing evidence of the existence of a checkerboard supersolid at a {\it commensurate} filling factor 1/2 (commensurate supersolid) in the soft-core Bose-Hubbard model with nearest-neighbor repulsions on a cubic lattice. In conventional cases, supersolids are realized at incommensurate filling factors by a doped-defect-condensation mechanism, where particles (holes) doped into a perfect crystal act as interstitials (vacancies) and delocalize in the crystal order. However, in the above model, a supersolid state is stabilized even at the commensurate filling factor 1/2 {\it without doping}. By performing grand canonical simulations, we obtain a ground-state phase diagram that suggests the existence of a supersolid at a commensurate filling. To obtain direct evidence of the commensurate supersolid, we next perform simulations in canonical ensembles at a particle density $ρ=1/2$ and exclude the possibility of phase separation. From the obtained snapshots, we discuss the microscopic structure and observe that interstitial-vacancy pairs are unbound in the crystal order.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Entanglement Spectra of the 2D AKLT Model: VBS/CFT Correspondence

We investigate the entanglement properties of the valence-bond-solid (VBS) state defined on two-dimensional lattices, which is the exact ground state of the Affleck-Kennedy-Lieb-Tasaki model. It is shown that the entanglement entropy obeys an area law and the non-universal prefactor of the leading term is strictly less than $\ln 2$. The analysis of entanglement spectra for various lattices reveals that the reduced density matrix associated with the VBS state is closely related to a thermal density matrix of a {\it holographic} spin chain, whose spectrum is reminiscent of that of the spin-1/2 Heisenberg chain. This correspondence is further supported by comparing the entanglement entropy in the holographic spin chain with conformal field theory predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ground-State Phase Diagram of the Two-Dimensional Extended Bose-Hubbard Model

We investigate the ground-state phase diagram of the soft-core Bose-Hubbard model with the nearest-neighbor repulsion on a square lattice by using an unbiased quantum Monte Carlo method. In contrast to the previous study[P. Sengupta {\it et. al.}, Phys. Rev. Lett. {\bf 94}, 207202 (2005)], we present the ground-state phase diagrams up to large hopping parameters. As a result, in addition to the known supersolid above half-filling, we find supersolid even below and at half-filling for large hopping parameters. Furthermore, for the strong nearest-neighbor repulsion, we show that the supersolid phase occupies a remarkably broad region in the phase diagram. The results are in qualitative agreement with that obtained by the Gutzwiller mean-field approximation[M. Iskin, Phys. Rev. A {\bf 83}, 051606(R) (2011) and T. Kimura, Phys. Rev. A {\bf 84}, 063630 (2011)].

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quantum phases of hard-core bosons on 2D lattices with anisotropic dipole-dipole interaction

By using an unbiased quantum Monte Carlo method, we investigate the hard-core Bose-Hubbard model on a square lattice with anisotropic dipole-dipole interaction. To study the effect of the anisotropy, dipole moments are assumed to be polarized in the $y$ direction on the two-dimensional (2D) $xy$ plane. To perform efficient simulations of long-range interacting systems, we use the worm algorithm with an O(N) Monte Carlo method. We obtain the ground-state phase diagram that includes a superfluid phase and a striped solid phase at half-filling as two main phases. In addition to these two main phases, we find a small region where there are multi plateaus in the particle density for small hopping amplitudes. In this region, the number of plateaus increases as the system size increases. This indicates the appearance of numerous competing ground states due to frustrated interactions. In our simulations, we find no evidence of a supersolid phase.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Study of the Shastry Sutherland Model Using Multi-scale Entanglement Renormalization Ansatz

We performed variational calculation based on the multi-scale entanglemnt renormalization ansatz, for the antiferromagnetic Heisenberg model on a Shastry Sutherland lattice (SSL). Our results show that at coupling ratio J'/J= 0.687(3), the system undergoes a quantum phase transition from the orthogonal dimer order to the plaquette valence bond solid phase, which then transits into the antiferromagnetic order above J'/J=0.75. In the presence of an external magnetic field, our calculations show clear evidences of various magnetic plateaux in systems with different coupling ratios range from 0.5 to 0.69. Our calculations are not limited to the small coupling ratio region, and we are able to show strong evidence of the presence of several supersolid phases, including ones above 1/2 and 1/3 plateaux. Such supersolid phases, which feature the coexistence of compressible superfluidity and crystalline long range order in triplet excitations, emerge at relatively large coupling ratio (J'/J>0.5). A schematic phase diagram of the SSL model in the presence of magnetic field is provided.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Double-$\it q$ Order in a Frustrated Random Spin System

We use the three-dimensional Heisenberg model with site randomness as an effective model of the compound Sr(Fe$_{1-x}$Mn$_x$)O$_2$. The model consists of two types of ions that correspond to Fe and Mn ions. The nearest-neighbor interactions in the ab-plane are antiferromagnetic. The nearest-neighbor interactions along the c-axis between Fe ions are assumed to be antiferromagnetic, whereas other interactions are assumed to be ferromagnetic. From Monte Carlo simulations, we confirm the existence of the double-$\boldsymbol{q}$ ordered phase characterized by two wave numbers, $(πππ)$ and $(π\pi0)$. We also identify the spin ordering pattern in the double-$\boldsymbol{q}$ ordered phase.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

First-Order Phase Transition with Breaking of Lattice Rotation Symmetry in Continuous-Spin Model on Triangular Lattice

Using a Monte Carlo method, we study the finite-temperature phase transition in the two-dimensional classical Heisenberg model on a triangular lattice with or without easy-plane anisotropy. The model takes account of competing interactions: a ferromagnetic nearest-neighbor interaction $J_1$ and an antiferromagnetic third nearest-neighbor interaction $J_3$. As a result, the ground state is a spiral spin configuration for $-4 < J_1/J_3 < 0$. In this structure, global spin rotation cannot compensate for the effect of 120-degree lattice rotation, in contrast to the conventional 120-degree structure of the nearest-neighbor interaction model. We find that this model exhibits a first-order phase transition with breaking of the lattice rotation symmetry at a finite temperature. The transition is characterized as a $Z_2$ vortex dissociation in the isotropic case, whereas it can be viewed as a $Z$ vortex dissociation in the anisotropic case. Remarkably, the latter is continuously connected to the former as the magnitude of anisotropy decreases, in contrast to the recent work by Misawa and Motome [J. Phys. Soc. Jpn. \textbf{79} (2010) 073001.] in which both the transitions were found to be continuous.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Novel Mechanism of Supersolid of Ultracold Polar Molecules in Optical Lattices

We study the checkerboard supersolid of the hard-core Bose-Hubbard model with the dipole-dipole interaction. This supersolid is different from all other supersolids found in lattice models in the sense that superflow paths through which interstitials or vacancies can hop freely are absent in the crystal. By focusing on repulsive interactions between interstitials, we reveal that the long-range tail of the dipole-dipole interaction have the role of increasing the energy cost of domain wall formations. This effect produces the supersolid by the second-order hopping process of defects. We also perform exact quantum Monte Carlo simulations and observe a novel double peak structure in the momentum distribution of bosons, which is a clear evidence for supersolid. This can be measured by the time-of-flight experiment in optical lattice systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Phase Transition of Generalized Ferromagnetic Potts Model - Effect of Invisible States -

We investigate the nature of the phase transition of the ferromagnetic Potts model with invisible states. The ferromagnetic Potts model with invisible states can be regarded as straightforward extension of the standard ferromagnetic Potts model. The invisible states contribute the entropy, however they do not affect the internal energy. They also do not change the symmetry which breaks at the transition temperature. The invisible states stimulate a first-order phase transition. We confirm that the first-order phase transition with spontaneous $q$-fold symmetry breaking for $q=2,3$, and 4 takes place even on two-dimensional lattice by Monte Carlo simulation. We also find that the transition temperature decreases and the latent heat increases as the number of invisible states increases.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Random Fan-Out State Induced by Site-Random Interlayer Couplings

We study the low-temperature properties of a classical Heisenberg model with site-random interlayer couplings on the cubic lattice. This model is introduced as a simplified effective model of Sr(Fe$_{1-x}$Mn$_{x}$)O$_2$, which was recently synthesized. In this material, when $x=0.3$, $(πππ)$ and $(π\pi0)$ mixed ordering is observed by neutron diffraction measurements. By Monte Carlo simulations, we find an exotic bulk spin structure that explains the experimentally obtained results. We name this spin structure the "random fan-out state". The mean-field calculations provide an intuitive understanding of this phase being induced by the site-random interlayer couplings. Since Rietveld analysis assuming the random fan-out state agrees well with the neutron diffraction pattern of Sr(Fe$_{0.7}$Mn$_{0.3}$)O$_2$, we conclude that the random fan-out state is reasonable for the spin-ordering pattern of Sr(Fe$_{0.7}$Mn$_{0.3}$)O$_2$ at the low-temperature phase.

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

Finite-Temperature Transition in the Spin-Dimer Antiferromagnet BaCuSi2O6

We consider a classical XY-like Hamiltonian on a body-centered tetragonal lattice, focusing on the role of interlayer frustration. A three-dimensional (3D) ordered phase is realized via thermal fluctuations, breaking the mirror-image reflection symmetry in addition to the XY symmetry. A heuristic field-theoretical model of the transition has a decoupled fixed point in the 3D XY universality, and our Monte Carlo simulation suggests that there is such a temperature region where long-wavelength fluctuations can be described by this fixed point. However, it is shown using scaling arguments that the decoupled fixed point is unstable against a fluctuation-induced biquadratic interaction, indicating that a crossover to nontrivial critical phenomena with different exponents appears as one approaches the critical point beyond the transient temperature region. This new scenario clearly contradicts the previous notion of the 3D XY universality.

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

Quantum Monte Carlo method for pairing phenomena: Super-counter-fluid of two-species Bose gases in optical lattices

We study the super-counter-fluid(SCF) states in the two-component hardcore Bose-Hubbard model on the square lattice, using the quantum Monte Carlo method based on the worm(directed loop) algorithm. Since the SCF state is a state of a pair-condensation characterized by $< a^{\dagger} b > \neq 0, < a > = 0$, and $< b > = 0$, where $a$ and $b$ are the order parameters of the two components, it is important to study behaviors of the pair-correlation function $< a_{i} b_{i}^{\dagger} a_{j}^{\dagger} b_{j} >$. For this purpose, we propose a choice of the worm head for calculating the pair-correlation function. From this pair-correlation, we confirm the Kosterlitz-Thouless(KT) charactor of the SCF phase. The simulation efficiency is also improved in the SCF phase.

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

Direct mapping of the finite temperature phase diagram of strongly correlated quantum models

Optical lattice experiments, with the unique potential of tuning interactions and density, have emerged as emulators of nontrivial theoretical models that are directly relevant for strongly correlated materials. However, so far the finite temperature phase diagram has not been mapped out for any strongly correlated quantum model. We propose a remarkable method for obtaining such a phase diagram for the first time directly from experiments using only the density profile in the trap as the input. We illustrate the procedure explicitly for the Bose Hubbard model, a textbook example of a quantum phase transition from a superfluid to a Mott insulator. Using "exact" quantum Monte Carlo simulations in a trap with up to $10^6$ bosons, we show that kinks in the local compressibility, arising from critical fluctuations, demarcate the boundaries between superfluid and normal phases in the trap. The temperature of the bosons in the optical lattice is determined from the density profile at the edge. Our method can be applied to other phase transitions even when reliable numerical results are not available.

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
preprint2009arXiv

Magnetic properties of the spin-1/2 XXZ model on the Shastry-Sutherland lattice: Effect of long-range interactions

We study magnetic properties of the $S=1/2$ Ising-like XXZ model on the Shastry-Sutherland lattices with long-range interactions, using the quantum Monte Carlo method. This model shows magnetization plateau phases at one-half and one-third of the saturation magnetization when additional couplings are considered. We investigate the finite temperature transition to one-half and one-third plateau phases. The obtained results suggest that the former case is of the first order and the latter case is of the second order. We also find that the system undergoes two successive transitions with the 2D Ising model universality, although there is a single phase transition in the Ising limit case. Finally, we estimate the coupling ratio to explain the magnetization process observed in ${\rm TmB_4}$

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Successive phase transitions and phase diagrams of the quasi-two-dimensional triangular antiferromagnet Rb4Mn(MoO4)3

Comprehensive experimental studies by magnetic, thermal and neutron measurements have clarified that Rb4Mn(MoO4)3 is a model system of a quasi-2D triangular Heisenberg antiferromagnet with an easy-axis anisotropy, exhibiting successive transitions across an intermediate collinear phase. As a rare case for geometrically frustrated magnetism, quantitative agreement between experiment and theory is found for complete, anisotropic phase diagrams as well as magnetic properties.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Dimer-Quadrupolar Quantum Phase Transition in the Quasi-One-Dimensional Heisenberg Model with Biquadratic Interaction

The quasi-one-dimensional S=1 Heisenberg antiferromagnet with a biquadratic term is investigated at zero temperature by quantum Monte Carlo simulation. As the magnitude of the inter-chain coupling is increased, the system undergoes a phase transition from a spontaneously dimerized phase to a Néel ordered or spin nematic phase. The numerical results suggest the possibility of an unconventional second-order transition in which the symmetry group characterizing one phase is not a subgroup of the other.

0 citations0 reviewsNo code linkedOpen access