Hard-core Boson Model Research Articles

Supervised learning of an interacting two-dimensional hardcore boson model of a weak topological insulator using correlation functions

Supervised learning of an interacting two-dimensional hardcore boson model of a weak topological insulator using correlation functions

Strong Pairing Originated from an Emergent Z_{2} Berry Phase in La_{3}Ni_{2}O_{7}.

The recent discovery of high-temperature superconductivity in La_{3}Ni_{2}O_{7} offers a fresh platform for exploring unconventional pairing mechanisms. Starting with the basic argument that the electrons in d_{z^{2}} orbitals nearly form local moments, we examine the effect of the Hubbard interaction U on the binding strength of Cooper pairs based on a single-orbital bilayer model with intralayer hopping t_{∥} and interlayer superexchange J_{⊥}. By extensive density matrix renormalization group calculations, we observe a remarkable enhancement in binding energy as much as 10-20 times larger with U/t_{∥} increasing from 0 to 12 at J_{⊥}/t_{∥}∼1. We demonstrate that such a substantial enhancement stems from a kinetic-energy-driven mechanism. Specifically, a Z_{2} Berry phase will emerge at large U due to the Hilbert space restriction (Mottness), which strongly suppresses the mobility of single particle propagation as compared to U=0. However, the kinetic energy of the electrons (holes) can be greatly restored by forming an interlayer spin-singlet pairing, which naturally results in a superconducting state even for relatively small J_{⊥}. An effective hard-core bosonic model is further proposed to estimate the superconducting transition temperature at the mean-field level.

Many-body localization on finite generation fractal lattices

We study many-body localization in a hardcore boson model in the presence of random disorder on finite generation fractal lattices with different Hausdorff dimensions and different local lattice structures. In particular, we consider the Vicsek, T-shaped, Sierpinski gasket, and modified Koch-curve fractal lattices. In the single-particle case, these systems display Anderson localization for arbitrary disorder strength if they are large enough. In the many-body case, the systems available to exact diagonalization exhibit a transition between a delocalized and localized regime, visible in the spectral and entanglement properties of these systems. The position of this transition depends on the Hausdorff dimension of the given fractal, as well as on its local structure.

Many-body entanglement and spectral clusters in the extended hard-core bosonic Hatano-Nelson model

Many-body entanglement and spectral clusters in the extended hard-core bosonic Hatano-Nelson model

Attaining the maximum Bose–Einstein condensation in a finite-size hard-core boson model

We propose a model for hard-core bosons in a lattice which allows to achieve the optimal occupation number predicted by Tennie et al (2017 Phys. Rev. B 96 064502) for a finite number of sites. The model is based on an extension of the Hamiltonian of the so-called Hubbard star, whose quantum properties are studied by means of quantum information descriptors such as the von Neumann entropy and the mutual information. These metrics are analyzed as a function of the one- and two-particle reduced density matrices, allowing to explore the relationship between condensation and entanglement by means of a control parameter that, under a given limit, connects our findings with previous results. All developments comprised in this article have been derived by analytical methods.

Tensor network investigation of the hard-square model

Using the corner-transfer matrix renormalization group to contract the tensor network that describes its partition function, we investigate the nature of the phase transitions of the hard-square model, one of the exactly solved models of statistical physics for which Baxter has found an integrable manifold. The motivation is twofold: assess the power of tensor networks for such models, and probe the 2D classical analog of a 1D quantum model of hard-core bosons that has recently attracted significant attention in the context of experiments on chains of Rydberg atoms. Accordingly, we concentrate on two planes in the 3D parameter space spanned by the activity and the coupling constants in the two diagonal directions. We first investigate the only case studied so far with Monte Carlo simulations, the case of opposite coupling constants. We confirm that, away and not too far from the integrable 3-state Potts point, the transition out of the period-3 phase appears to be unique in the Huse-Fisher chiral universality class, albeit with significantly different exponents as compared to Monte Carlo. We also identify two additional phase transitions not reported so far for that model, a Lifshitz disorder line, and an Ising transition for large enough activity. To make contact with 1D quantum models of Rydberg atoms, we then turn to a plane where the ferromagnetic coupling is kept fixed, and we show that the resulting phase diagram is very similar, the only difference being that the Ising transition becomes first-order through a tricritical Ising point, in agreement with Baxter's prediction that this plane should contain a tricritical Ising point, and in remarkable, almost quantitative agreement with the phase diagram of the 1D quantum version of the model.

Beyond Hard-Core Bosons in Transmon Arrays

Arrays of transmons have proven to be a viable medium for quantum information science and quantum simulations. Despite their widespread popularity as qubit arrays, there remains yet untapped potential beyond the two-level approximation or, equivalently, the hard-core boson model. With the higher excited levels included, coupled transmons naturally realize the attractive Bose-Hubbard model. The dynamics of the model has been difficult to study due to unfavorable scaling of the dimensionality of the Hilbert space with the system size. In this work, we present a framework for describing the effective unitary dynamics of highly-excited states of coupled transmons based on high-order degenerate perturbation theory. This allows us to describe various collective phenomena -- such as bosons stacked onto a single site behaving as a single particle, edge-localization, and effective longer-range interactions -- in a unified, compact and accurate manner. While our examples deal with one-dimensional chains of transmons for the sake of clarity, the theory can be readily applied to more general geometries.

Critical Temperatures of Hard-Core Boson Model on Square Lattice within Bethe Approximation

The short-range correlations are considered for a two-dimensional hard-core boson model on square lattice within Bethe approximation for the clusters consisting of two and four sites. Explicit equations are derived for the critical temperatures of charge and superfluid ordering and their solutions are considered for various ratios of the charge-charge correlation parameter to the transfer integral. It is shown that taking into account short-range correlations for the temperatures of charge ordering results in the appearance of the critical concentration of bosons, which restricts the existence domain of the solutions of charge ordering type. In the case of superfluid ordering with the assumption of short-range correlations, the critical temperature is reduced up to zero values at half-filling. A phase diagram of the hard-core boson model is constructed with the assumption of phase separation within Maxwell's construction and it is shown that taking into account short-range correlations within Bethe approximation quantitatively approximates the form of phase diagram to the results of the quantum Monte Carlo method.

Stable iPEPO Tensor-Network Algorithm for Dynamics of Two-Dimensional Open Quantum Lattice Models

Being able to describe accurately the dynamics and steady-states of driven and/or dissipative but quantum correlated lattice models is of fundamental importance in many areas of science: from quantum information to biology. An efficient numerical simulation of large open systems in two spatial dimensions is a challenge. In this work, we develop a tensor network method, based on an infinite Projected Entangled Pair Operator (iPEPO) ansatz, applicable directly in the thermodynamic limit. We incorporate techniques of finding optimal truncations of enlarged network bonds by optimising an objective function appropriate for open systems. Comparisons with numerically exact calculations, both for the dynamics and the steady-state, demonstrate the power of the method. In particular, we consider dissipative transverse quantum Ising and driven-dissipative hard core boson models in non-mean field limits, proving able to capture substantial entanglement in the presence of dissipation. Our method enables to study regimes which are accessible to current experiments but lie well beyond the applicability of existing techniques.

Extended Coulomb liquid of paired hardcore boson model on a pyrochlore lattice

There is a growing interest in the $U(1)$ Coulomb liquid in both quantum materials in pyrochlore ice and cluster Mott insulators and cold atom systems. We explore a paired hardcore boson model on a pyrochlore lattice. This model is equivalent to the XYZ spin model that was proposed for rare-earth pyrochlores with "dipole-octupole" doublets. Since this model has no sign problem for quantum Monte Carlo (QMC) simulations in a large parameter regime, we carry out both analytical and QMC calculations. We find that the $U(1)$ Coulomb liquid is quite stable and spans a rather large portion of the phase diagram with boson pairing. Moreover, we numerically find thermodynamic evidence that the boson pairing could induce a possible $\mathbb{Z}_2$ liquid in the vicinity of the phase boundary between Coulomb liquid and $\mathbb{Z}_2$ symmetry-broken phase. Besides the materials' relevance with quantum spin ice, we point to quantum simulation with cold atoms on optical lattices.

Worm-algorithm-type simulation of the quantum transverse-field Ising model

We apply a worm algorithm to simulate the quantum transverse-field Ising model in a path-integral representation of which the expansion basis is taken as the spin component along the external-field direction. In such a representation, a configuration can be regarded as a set of non-intersecting loops constructed by "kinks" for pairwise interactions and spin-down (or -up) imaginary-time segments. The wrapping probability for spin-down loops, a dimensionless quantity characterizing the loop topology on a torus, is observed to exhibit small finite-size corrections and yields a high-precision critical point in two dimensions (2D) as $h_c \! =\! 3.044\, 330(6)$, significantly improving over the existing results and nearly excluding the best one $h_c \! =\! 3.044\, 38 (2)$. At criticality, the fractal dimensions of the loops are estimated as $d_{\ell \downarrow} (1{\rm D}) \! = \! 1.37(1) \! \approx \! 11/8 $ and $d_{\ell \downarrow} (2{\rm D}) \! = \! 1.75 (3)$, consistent with those for the classical 2D and 3D O(1) loop model, respectively. An interesting feature is that in one dimension (1D), both the spin-down and -up loops display the critical behavior in the whole disordered phase ($ 0 \! \leq \! h \! < \! h_c$), having a fractal dimension $d_{\ell} \! = \! 1.750 (7)$ that is consistent with the hull dimension $d_{\rm H} \! = \! 7/4$ for critical 2D percolation clusters. The current worm algorithm can be applied to simulate other quantum systems like hard-core boson models with pairing interactions.

Guide to Exact Diagonalization Study of Quantum Thermalization

Exact diagonalization is a powerful numerical method to study isolated quantum many-body systems. This paper provides a review of numerical algorithms to diagonalize the Hamiltonian matrix. Symmetry and the conservation law help us perform the numerical study efficiently. We explain the method to block-diagonalize the Hamiltonian matrix by using particle number conservation, translational symmetry, particle-hole symmetry, and spatial reflection symmetry in the context of the spin-1/2 XXZ model or the hard-core boson model in a one-dimensional lattice. We also explain the method to study the unitary time evolution governed by the Schr\"odinger equation and to calculate thermodynamic quantities such as the entanglement entropy. As an application, we demonstrate numerical results that support that the eigenstate thermalization hypothesis holds in the XXZ model.

Persistent oscillations versus thermalization in the quench dynamics of quantum gases with long-range interactions

Searching for nonthermalized dynamics in interacting quantum systems is not only of fundamental theoretical interest in nonequilibrium quantum physics, but also of immense practical significance in quantum information processing. In this paper, we study quantum quench dynamics in an hard-core bosonic model with infinite-range interactions, which have been realized in recent high-finesse cavity experiments. We show the long-time dynamics of this model can exhibit either undamped oscillations or thermalization depending on the choice of initial states. The long-range nature of the interactions rather than conserved quantities are responsible for such nonergodic dynamical behaviors.

Biexciton Condensation in Electron-Hole-Doped Hubbard Bilayers: A Sign-Problem-Free Quantum MonteCarlo Study.

The bilayer Hubbard model with electron-hole doping is an ideal platform to study excitonic orders due to suppressed recombination via spatial separation of electrons and holes. However, suffering from the sign problem, previous quantum MonteCarlo studies could not arrive at an unequivocal conclusion regarding the presence of phases with clear signatures of excitonic condensation in bilayer Hubbard models. Here, we develop a determinant quantum MonteCarlo algorithm for the bilayer Hubbard model that is sign-problem-free for equal and opposite doping in the two layers and study excitonic order and charge and spin density modulations as a function of chemical potential difference between the two layers, on-site Coulomb repulsion, and interlayer interaction. In the intermediate coupling regime and in proximity to the SU(4)-symmetric point, we find a biexcitonic condensate phase at finite electron-hole doping, as well as a competing (π,π) charge density wave state. We extract the Berezinskii-Kosterlitz-Thouless transition temperature from superfluid density and a finite-size scaling analysis of the correlation functions and explain our results in terms of an effective biexcitonic hard-core boson model.

Pair hopping in systems of strongly interacting hard-core bosons

We have used the Stochastic Series Expansion quantum Monte Carlo method to study interacting hard-core bosons on the square lattice, with pair-hopping processes supplementing the standard single-particle hopping. Such pair hopping arises in effective models for frustrated quantum magnets. Our goal is to investigate the effects of the pair hopping process on the commonly observed superfluid, insulating (Mott), and super-solid ground-state phases in the standard hard-core boson model with various interaction terms. The model is specifically motivated by the observation of finite dispersion of 2-magnon bound states in neutron diffraction experiments SrCu$_2$(BO$_3$)$_2$. Our results show that the pair hopping has different effects on Mott phases at different filling fractions, "melting" them at different critical pair-hopping amplitudes. Thus, it appears that pair hopping may have an important role in determining which out of a potentially large number of Mott phases (stabilized by details of the charge-diagonal interaction terms) actually survive the totality of quantum fluctuations present.

Diagnosing Potts criticality and two-stage melting in one-dimensional hard-core boson models

We investigate a model of hard-core bosons with infinitely repulsive nearest- and next-nearest-neighbor interactions in one dimension, introduced by Fendley, Sengupta and Sachdev in Phys. Rev. B 69, 075106 (2004). Using a combination of exact diagonalization, tensor network, and quantum Monte Carlo simulations, we show how an intermediate incommensurate phase separates a crystalline and a disordered phase. We base our analysis on a variety of diagnostics, including entanglement measures, fidelity susceptibility, correlation functions, and spectral properties. According to theoretical expectations, the disordered-to-incommensurate-phase transition point is compatible with Berezinskii-Kosterlitz-Thouless universal behaviour. The second transition is instead non-relativistic, with dynamical critical exponent $z > 1$. For the sake of comparison, we illustrate how some of the techniques applied here work at the Potts critical point present in the phase diagram of the model for finite next-nearest-neighbor repulsion. This latter application also allows to quantitatively estimate which system sizes are needed to match the conformal field theory spectra with experiments performing level spectroscopy.

Out-of-time-ordered correlators in short-range and long-range hard-core boson models and in the Luttinger-liquid model

We study out-of-time-ordered correlators (OTOC) in hard-core boson models with short-range and long-range hopping and compare the results to the OTOC in the Luttinger liquid model. For the density operator, a related `commutator function' starts at zero and decays back to zero after the passage of the wavefront in all three models, while the wavefront broadens as $t^{1/3}$ in the short-range model and shows no broadening in the long-range model and the Luttinger liquid model. For the boson creation operator, the corresponding commutator function shows saturation inside the light cone in all three models, with similar wavefront behavior as in the density-density commutator function, despite the presence of a nonlocal string in terms of Jordan-Wigner fermions. For the long-range model and the Luttinger liquid model, the commutator function decays as power law outside the light cone in the long time regime when following different fixed-velocity rays. In all cases, the OTOCs approach their long-time values in a power-law fashion, with different exponents for different observables and short-range vs long-range cases. Our long-range model appears to capture exponents in the Luttinger liquid model (which are found to be independent of the Luttinger parameter in the model). This conclusion also bears on the OTOC calculations in conformal field theories, which we propose correspond to long-ranged models.

Numerical study of the chiral Z3 quantum phase transition in one spatial dimension

Recent experiments on a one-dimensional chain of trapped alkali atoms [arXiv:1707.04344] have observed a quantum transition associated with the onset of period-3 ordering of pumped Rydberg states. This spontaneous $\mathbb{Z}_3$ symmetry breaking is described by a constrained model of hard-core bosons proposed by Fendley $et\, \,al.$ [arXiv:cond-mat/0309438]. By symmetry arguments, the transition is expected to be in the universality class of the $\mathbb{Z}_3$ chiral clock model with parameters preserving both time-reversal and spatial-inversion symmetries. We study the nature of the order-disorder transition in these models, and numerically calculate its critical exponents with exact diagonalization and density-matrix renormalization group techniques. We use finite-size scaling to determine the dynamical critical exponent $z$ and the correlation length exponent $\nu$. Our analysis presents the only known instance of a strongly-coupled transition between gapped states with $z \ne 1$, implying an underlying nonconformal critical field theory.

Ground-state phase diagram of the one-dimensional t−J model with pair hopping terms

The $t$-$J$ model is a standard model of strongly correlated electrons, often studied in the context of high-$T_c$ superconductivity. However, most studies of this model neglect three-site terms, which appear at the same order as the superexchange $J$. As these terms correspond to pair-hopping, they are expected to play an important role in the physics of superconductivity when doped sufficiently far from half-filling. In this paper we present a density matrix renormalisation group study of the one-dimensional $t$-$J$ model with the pair hopping terms included. We demonstrate that that these additional terms radically change the one-dimensional ground state phase diagram, extending the superconducting region at low fillings, while at larger fillings, superconductivity is completely suppressed. We explain this effect by introducing a simplified effective model of repulsive hardcore bosons.

Continuous Easy-Plane Deconfined Phase Transition on the Kagome Lattice.

We use large scale quantum MonteCarlo simulations to study an extended Hubbard model of hard core bosons on the kagome lattice. In the limit of strong nearest-neighbor interactions at 1/3 filling, the interplay between frustration and quantum fluctuations leads to a valence bond solid ground state. The system undergoes a quantum phase transition to a superfluid phase as the interaction strength is decreased. It is still under debate whether the transition is weakly first order or represents an unconventional continuous phase transition. We present a theory in terms of an easy plane noncompact CP^{1} gauge theory describing the phase transition at 1/3 filling. Utilizing large scale quantum MonteCarlo simulations with parallel tempering in the canonical ensemble up to 15552 spins, we provide evidence that the phase transition is continuous at exactly 1/3 filling. A careful finite size scaling analysis reveals an unconventional scaling behavior hinting at deconfined quantum criticality.