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

Seiji Yunoki

Seiji Yunoki contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.str-elquant-phcond-mat.stat-mechcond-mat.supr-conhep-latcond-mat.mtrl-sciMachine Learningphysics.comp-ph

Trust snapshot

Quick read

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

22 published item(s)

preprint2026arXiv

Dissipative ground-state preparation of a quantum spin chain on a trapped-ion quantum computer

We demonstrate a dissipative protocol for ground-state preparation of a quantum spin chain on a trapped-ion quantum computer. As a first step, we derive a Kraus representation of a dissipation channel for the protocol recently proposed by Ding et al. [Phys. Rev. Res. 6, 033147 (2024)] that still holds for arbitrary temporal discretization steps, extending the analysis beyond the Lindblad dynamics regime. The protocol guarantees that the fidelity with the ground state monotonically increases (or remains unchanged) under repeated applications of the channel to an arbitrary initial state, provided that the ground state is the unique steady state of the dissipation channel. Using this framework, we implement dissipative ground-state preparation of a transverse-field Ising chain for up to 19 spins on the trapped-ion quantum computer Reimei provided by Quantinuum. Despite the presence of hardware noise, the dynamics consistently converges to a low-energy state far away from the maximally mixed state even when the corresponding quantum circuits contain as many as 4110 entangling gates, demonstrating the intrinsic robustness of the protocol. By applying zero-noise extrapolation, the resulting energy expectation values are systematically improved to agree with noiseless simulations within statistical uncertainties.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Extending the Handover-Iterative VQE to Challenging Strongly Correlated Systems: $N_2$ and Fe-S Cluster

Accurately describing strongly correlated electronic systems remains a central challenge in quantum chemistry, as electron-electron interactions give rise to complex many-body wavefunctions that are difficult to capture with conventional approximations. Classical wavefunction-based approaches, such as the Semistochastic Heat-bath Configuration Interaction (SHCI) and the Density Matrix Renormalization Group (DMRG), currently define the state of the art, systematically converging toward the Full Configuration Interaction (FCI) limit, but at a rapidly increasing computational cost. Quantum computing algorithms promise to alleviate this scaling bottleneck by leveraging entanglement and superposition to represent correlated states more compactly. We introduced the Handover-Iterative Variational Quantum Eigensolver (HI-VQE) as a practical quantum computing algorithm with an iterative "handover" mechanism that dynamically exchanges information between quantum and classical computers, even using Noisy Intermediate-Scale Quantum (NISQ) computers. In this work, we extend the HI-VQE to benchmark two prototypical strongly correlated systems, the nitrogen molecule $N_2$ and iron-sulfur (Fe-S) cluster, which serve as stringent tests for both classical and quantum electronic-structure methods. By comparing HI-VQE results against Heat-bath Configuration Interaction (HCI) benchmarks, we assess its accuracy, scalability, and ability to capture multireference correlation effects. Achieving quantitative agreement on these canonical systems demonstrates a viable pathway toward quantum-enhanced simulations of complex bioinorganic molecules, catalytic mechanisms, and correlated materials.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Low-energy photoexcitations inside the Mott gap in doped Hubbard and t-J ladders

We investigate changes in the optical conductivity of doped Mott insulators by tuning ultrashort pump pulses to target either the Drude or low-energy absorption regions. Using a hole-doped two-leg Hubbard ladder and a four-leg t-J ladders, we calculate the optical conductivity after pump by employing the time-dependent density matrix renormalization group. We find that a monocycle electric field pulse tuned to the Drude absorption reduces the Drude weight, accompanied by a slight enhancement in the mid-infrared (mid-IR) spectral weight. However, this enhancement diminishes as the pulse intensity increases. In contrast, a pump pulse tuned to the mid-IR absorption only affects the Drude weight. This behavior arises because the mid-IR absorption originates from magnetic excitations that do not couple directly to photons. These predictions can be tested experimentally by applying ultrashort low-energy pump pulses to cuprate materials.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Tensor Computing Interface: An Application-Oriented, Lightweight Interface for Portable High-Performance Tensor Network Applications

Tensor networks (TNs) are a central computational tool in quantum science and artificial intelligence. However, the lack of unified software interface across tensor-computing frameworks severely limits the portability of TN applications, coupling algorithmic development to specific hardware and software back ends. To address this challenge, we introduce the Tensor Computing Interface (TCI) -- an application-oriented, lightweight application programming interface designed to enable framework-independent, high-performance TN applications. TCI provides a well-defined type system that abstracts tensor objects together with a minimal yet expressive set of core functions covering essential tensor manipulations and tensor linear-algebra operations. Through numerical demonstrations on representative tensor-network applications, we show that codes written against TCI can be migrated seamlessly across heterogeneous hardware and software platforms while achieving performance comparable to native framework implementations. We further release an open-source implementation of TCI based on \textit{Cytnx}, demonstrating its practicality and ease of integration with existing tensor-computing frameworks.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gutzwiller wave function on a quantum computer using a discrete Hubbard-Stratonovich transformation

We propose a quantum-classical hybrid scheme for implementing the nonunitary Gutzwiller factor using a discrete Hubbard-Stratonovich transformation, which allows us to express the Gutzwiller factor as a linear combination of unitary operators involving only single-qubit rotations, at the cost of the sum over the auxiliary fields. To perform the sum over the auxiliary fields, we introduce two approaches that have complementary features. The first approach employs a linear-combination-of-unitaries circuit, which enables one to probabilistically prepare the Gutzwiller wave function on a quantum computer, while the second approach uses importance sampling to estimate observables stochastically, similar to a quantum Monte Carlo method in classical computation. The proposed scheme is demonstrated with numerical simulations for the half-filled Fermi-Hubbard model. Furthermore, we perform quantum simulations using a real quantum device, demonstrating that the proposed scheme can reproduce the exact ground-state energy of the two-site Fermi-Hubbard model within error bars.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Local multiplet formation around a single vacancy in graphene: an effective Anderson model analysis based on the block-Lanczos DMRG method

To better understand the electronic structure of a single vacancy in graphene, we study the ground state property of an effective Anderson model, consisting of three dangling $sp^2$ orbitals of the surrounding carbon atoms around the vacancy and the $π$ orbitals of carbon atoms that form the honeycomb lattice with a single vacancy. Employing the block-Lanczos density-matrix renormalization group method, we show that there are two phases in the relevant parameter space, i.e., a nonmagnetic phase in the weak coupling region and a free magnetic moment phase in the realistic parameter region. The systematic analysis finds that, in the free magnetic moment phase, local multiplets of the doubly degenerate irreducible representation of the $\mathcal{C}_3$ point group with spin 1 become dominant in the ground state, and approximately half of this local spin 1 is screened by electrons in the surrounding $π$ orbitals, indicating the emergence of the residual spin-1/2 free magnetic moment. Furthermore, we find that the emergence of the free magnetic moment is robust against carrier doping, which is in sharp contrast to the case of graphene with an adatom, thus explaining the qualitative difference observed experimentally in these two classes of systems. We also find the enhancement of the spin correlation function between $π$ electrons around the vacancy and those in the conduction band away from the vacancy in the undoped case, as compared to that in the doped case, while the spin correlation function between the $σ$ electrons in the $sp^2$ dangling orbitals around the vacancy and the $π$ electrons in the conduction band remains large in both undoped and doped cases. Our calculations thus support qualitatively the previous experiment that suggests the emergence of free magnetic moment with two distinct origins.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Spatial, spin, and charge symmetry projections for a Fermi-Hubbard model on a quantum computer

We propose an extended version of the symmetry-adapted variational-quantum-eigensolver (VQE) and apply it to a two-component Fermi-Hubbard model on a bipartite lattice. In the extended symmetry-adapted VQE method, the Rayleigh quotient for the Hamiltonian and a parametrized quantum state in a properly chosen subspace is minimized within the subspace and is optimized among the variational parameters implemented on a quantum circuit to obtain variationally the ground state and the ground-state energy. The corresponding energy derivative with respect to a variational parameter is expressed as a Hellmann-Feynman-type formula of a generalized eigenvalue problem in the subspace, which thus allows us to use the parameter-shift rules for its evaluation. The natural-gradient-descent method is also generalized to optimize variational parameters in a quantum-subspace-expansion approach. As a subspace for approximating the ground state of the Hamiltonian, we consider a Krylov subspace generated by the Hamiltonian and a symmetry-projected variational state, and therefore the approximated ground state can restore the Hamiltonian symmetry that is broken in the parametrized variational state prepared on a quantum circuit. We show that spatial symmetry operations for fermions in an occupation basis can be expressed as a product of the nearest-neighbor fermionic swap operations on a quantum circuit. We also describe how the spin and charge symmetry operations, i.e., rotations, can be implemented on a quantum circuit. By numerical simulations, we demonstrate that the spatial, spin, and charge symmetry projections can improve the accuracy of the parametrized variational state, which can be further improved systematically by expanding the Krylov subspace without increasing the number of variational parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Discretized quantum adiabatic process for free fermions and comparison with the imaginary-time evolution

Motivated by recent progress of quantum technologies, we study a discretized quantum adiabatic process for a one-dimensional free fermion system described by a variational wave function, i.e., a parametrized quantum circuit. The wave function is composed of $M$ layers of two elementary sets of time-evolution operators, each set being decomposed into commutable local operators. The evolution time of each time-evolution operator is treated as a variational parameter so as to minimize the expectation value of the energy. We show that the exact ground state is reached by applying the layers of time-evolution operators as many as a quarter of the system size. This is the minimum number $M_B$ of layers set by the limit of speed, i.e., the Lieb-Robinson bound, for propagating quantum entanglement via the local time-evolution operators. Quantities such as the energy $E$ and the entanglement entropy $S$ of the optimized variational wave function with $M < M_B$ are independent of the system size $L$ but fall into some universal functions of $M$. The development of the entanglement in these ansatz is further manifested in the progressive propagation of single-particle orbitals in the variational wave function. We also find that the optimized variational parameters show a systematic structure that provides the optimum scheduling function in the quantum adiabatic process. We also investigate the imaginary-time evolution of this variational wave function, where the causality relation is absent due to the non-unitarity of the imaginary-time evolution operators, thus the norm of the wave function being no longer conserved. We find that the convergence to the exact ground state is exponentially fast, despite that the system is at the critical point, suggesting that implementation of the non-unitary imaginary-time evolution in a quantum circuit is highly promising to further shallow the circuit depth.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

QMC study of the chiral Heisenberg Gross-Neveu universality class

We investigate a quantum criticality of an antiferromagnetic phase transition in the Hubbard model on a square lattice with a $d$-wave pairing field by large-scale auxiliary-field quantum Monte Carlo simulations. Since the $d$-wave pairing filed induces Dirac cones in the non-interacting single-particle spectrum, the quantum criticality should correspond to the chiral Heisenberg universality class in terms of the Gross-Neveu theory, which is the same as those expected in the Hubbard model on the honeycomb lattice, despite the unit cells being different (e.g., they contain one and two sites, respectively). We show that both the two phase transitions, expected to occur on the square and on the honeycomb lattices, indeed have the same quantum criticality. We also argue that details of the models, i.e., the way of counting the total number $N$ of fermion components and the anisotropy of the Dirac cones, do not change the critical exponents.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Quantum Power Method by a Superposition of Time-Evolved States

We propose a quantum-classical hybrid algorithm of the power method, here dubbed as quantum power method, to evaluate $\hat{\cal H}^n |ψ\rangle$ with quantum computers, where $n$ is a nonnegative integer, $\hat{\cal H}$ is a time-independent Hamiltonian of interest, and $|ψ\rangle$ is a quantum state. We show that the number of gates required for approximating $\hat{\cal H}^n$ scales linearly in the power and the number of qubits, making it a promising application for near term quantum computers. Using numerical simulation, we show that the power method can control systematic errors in approximating the Hamiltonian power ${\hat{\cal H}^n}$ for $n$ as large as 100. As an application, we combine our method with a multireference Krylov-subspace-diagonalization scheme to show how one can improve the estimation of ground-state energies and the ground-state fidelities found using a variational-quantum-eigensolver scheme. Finally, we outline other applications of the quantum power method, including several moment-based methods. We numerically demonstrate the connected-moment expansion for the imaginary-time evolution and compare the results with the multireference Krylov-subspace diagonalization.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Quantum spin solver near saturation: QS$^3_{~}$

We develop a program package named QS$^{3}$ [\textipa{kjú:-és-kjú:b}] based on the (thick-restart) Lanczos method for analyzing spin-1/2 XXZ-type quantum spin models on spatially uniform/non-uniform lattices near fully polarized states, which can be mapped to dilute hardcore Bose systems. All calculations in QS$^{3}$, including eigenvalue problems, expectation values for one/two-point spin operators, and static/dynamical spin structure factors, are performed in the symmetry-adapted bases specified by the number $N_{\downarrow}$ of down spins and the wave number $\boldsymbol{k}$ associated with the translational symmetry without using the bit representation for specifying spin configurations. Because of these treatments, QS$^{3}$ can support large-scale quantum systems containing more than 1000 sites with dilute $N_{\downarrow}$. We show the benchmark results of QS$^{3}$ for the low-energy excitation dispersion of the isotropic Heisenberg model on the $10\times10\times10$ cubic lattice, the static and dynamical spin structure factors of the isotropic Heisenberg model on the $10\times10$ square lattice, and the open-MP parallelization efficiency on the supercomputer (Ohtaka) based on AMD Epyc 7702 installed at the Institute for the Solid State Physics (ISSP). Theoretical backgrounds and the user interface of QS$^{3}$ are also described.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Quantum-inspired search method for low-energy states of classical Ising Hamiltonians

We develop a quantum-inspired numerical procedure for searching low-energy states of a classical Hamiltonian composed of two-body fully-connected random Ising interactions and a random local longitudinal magnetic field. In this method, we introduce infinitesimal quantum interactions that do not commute with the original Ising Hamiltonian, and repeatedly generate and truncate direct product states, inspired by the Krylov subspace method, to obtain the low-energy states of the original classical Ising Hamiltonian. The computational cost is controlled by the form of infinitesimal quantum interactions (e.g., one-body or two-body interactions) and the numbers of infinitesimal interaction terms introduced, different initial states considered, and low-energy states kept during the iteration. For a demonstrate of the method, here we introduce as the infinitesimal quantum interactions pair products of Pauli $X$ operators acting on different sites and on-site Pauli $X$ operators into the random Ising Hamiltonian, in which the numerical cost is $O(N^3)$ per iteration with the system size $N$. We consider 120 instances of the random coupling realizations for the random Ising Hamiltonian with $N$ up to 600 and search the 120 lowest-energy states for each instance. We find that the time-to-solution by the quantum-inspired method proposed here, with parallelization in terms of the different initial states, for searching the ground state of the random Ising Hamiltonian scales approximately as $N^5$ for $N$ up to 600. We also examine the basic physical properties such as the ensemble-averaged ground-state and first-excited energies and the ensemble-averaged number of states in the low-energy region of the random Ising Hamiltonian.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Characterization of photoexcited states in the half-filled one-dimensional extended Hubbard model assisted by machine learning

Photoinduced nonequilibrium states can provide new insight into dynamical properties of strongly correlated electron systems. One of the typical and extensively studied systems is the half-filled one-dimensional extended Hubbard model (1DEHM). Here, we propose that the supervised machine learning (ML) can provide useful information for characterizing photoexcited states in 1DEHM. Using entanglement spectra as a training dataset, we construct neural network. Judging from the trained network, we find that bond-spin-density wave (BSDW) order can be enhanced in photoexcited states if the frequency of a driving pulse nearly resonates with gap. We separately calculate the time evolution of local and non-local order parameters and confirm that the correlation functions of BSDW are actually enhanced by photoexcitation as predicted by ML. The successful prediction of BSDW demonstrates the advantage of ML to assist characterizing photoexcited quantum states.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Charge stiffness and long-range correlation in the optically induced $η$-pairing state of the one-dimensional Hubbard model

We show that optical excitation of the Mott insulating phase of the one-dimensional Hubbard model can create a state possessing two of the hallmarks of superconductivity: a nonvanishing charge stiffness and long-ranged pairing correlation. By employing the exact diagonalization method, we find that the superposition of the $η$-pairing eigenstates induced by the optical pump exhibits a nonvanishing charge stiffness and a pairing correlation that decays very slowly with system size in sharp contrast to the behavior of an ensemble of thermally excited eigenstates, which has a vanishing charge stiffness and no long-ranged pairing correlations. We show that the charge stiffness is indeed directly associated with the $η$-pairing correlation in the Hubbard model. Our finding demonstrates that optical pumping can actually lead to superconducting-like properties on the basis of the $η$-pairing states.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Corner transfer matrix renormalization group analysis of the two-dimensional dodecahedron model

We investigate the phase transition of the dodecahedron model on the square lattice. The model is a discrete analogue of the classical Heisenberg model, which has continuous $O(3)$ symmetry. In order to treat the large on-site degree of freedom $q = 20$, we develop a massively parallelized numerical algorithm for the corner transfer matrix renormalization group method, incorporating EigenExa, the high-performance parallelized eigensolver. The scaling analysis with respect to the cutoff dimension reveals that there is a second-order phase transition at $T^{~}_{\rm c}=0.4398(8)$ with the critical exponents $ν=2.88(8)$ and $β=0.21(1)$. The central charge of the system is estimated as $c=1.99(6)$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Finite-$m$ scaling analysis of Berezinskii-Kosterlitz-Thouless phase transitions and entanglement spectrum for the six-state clock model

We investigate the Berezinskii-Kosterlitz-Thouless transitions for the square-lattice six-state clock model with the corner-transfer matrix renormalization group (CTMRG). Scaling analyses for effective correlation length, magnetization, and entanglement entropy with respect to the cutoff dimension $m$ at the fixed point of CTMRG provide transition temperatures consistent with a variety of recent numerical studies. We also reveal that the fixed point spectrum of the corner transfer matrix in the critical intermediate phase of the six-state clock model is characterized by the scaling dimension consistent with the $c=1$ boundary conformal field theory associated with the effective $Z_6$ dual sine-Gordon model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Photoinduced $η$ pairing in the Kondo lattice model

The previous theoretical study has shown that pulse irradiation to the Mott insulating state in the Hubbard model can induce the enhancement of superconducting correlation due to the generation of $η$ pairs. Here, we show that the same mechanism can be applied to the Kondo lattice model, an effective model for heavy electron systems, by demonstrating that the pulse irradiation indeed enhances the $η$-pairing correlation. As in the case of the Hubbard model, the non-linear optical process is essential to increase the number of photoinduced $η$ pairs and thus the enhancement of the superconducting correlation. We also find the diffusive behavior of the spin dynamics after the pulse irradiation, suggesting that the increase of the number of $η$ pairs leads to the decoupling between the conduction band and the localized spins in the Kondo lattice model, which is inseparably related to the photodoping effect.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Photoinduced $η$-pairing at finite temperatures

We numerically prove photoinduced $η$-pairing in a half-filled fermionic Hubbard chain at both zero and finite temperature. The result, obtained by combining the matrix-product-state based infinite time-evolving block decimation technique and the purification method, applies to the thermodynamic limit. Exciting the Mott insulator by a laser electric field docked on via the Peierls phase, we track the time-evolution of the correlated many-body system and determine the optimal parameter set for which the nonlocal part of the $η$-pair correlation function becomes dominant during the laser pump at zero and low temperatures. These correlations vanish at higher temperatures and long times after pulse irradiation. In the high laser frequency strong Coulomb coupling regime we observe a remnant enhancement of the Brillouin-zone boundary pair-correlation function also at high temperatures, if the Hubbard interaction is about a multiple of the laser frequency, which can be attributed to an enhanced double occupancy in the virtual Floquet state.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Photoinduced $η$-pairing in One-dimensional Mott Insulators

Employing the density-matrix renormalization group technique in the matrix-product-state representation, we investigate the photoexcited superconducting correlations induced by the $η$-pairing mechanism in the half-filled Hubbard chain. We estimate the characteristic pair correlation function and verify the accuracy of our numerical results by comparison with exact-diagonalization data for small systems. The optimal parameter set of the pump that most enhances the $η$-pair correlations, is calculated in the strong-coupling regime. For such a pump, we explore the possibility of quasi-long-range order.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Superconductivity and charge density wave under a time-dependent periodic field in the one-dimensional attractive Hubbard model

We investigate the competition between superconductivity (SC) and charge density wave (CDW) under a time-dependent periodic field in the attractive Hubbard model. By employing the time-dependent exact diagonalization method, we show that the driving frequency and amplitude of the external field can control the enhancement of either the superconducting pair or the CDW correlations, which are degenerate in the ground state of the half-filled attractive Hubbard model in the absence of the field. In the strong-coupling limit of the attractive Hubbard interaction, the controllability is characterized by the anisotropic interaction of the effective model. The anisotropy is induced by the external field and lifts the degeneracy of SC and CDW. We find that the enhancement or suppression of the superconducting pair and CDW correlations in the periodically-driven attractive Hubbard model can be well interpreted by the quench dynamics of the effective model derived in the strong-coupling limit.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Symmetry-adapted variational quantum eigensolver

We propose a scheme to restore spatial symmetry of Hamiltonian in the variational-quantum-eigensolver (VQE) algorithm for which the quantum circuit structures used usually break the Hamiltonian symmetry. The symmetry-adapted VQE scheme introduced here simply applies the projection operator, which is Hermitian but not unitary, to restore the spatial symmetry in a desired irreducible representation of the spatial group. The entanglement of a quantum state is still represented in a quantum circuit but the nonunitarity of the projection operator is treated classically as postprocessing in the VQE framework. By numerical simulations for a spin-$1/2$ Heisenberg model on a one-dimensional ring, we demonstrate that the symmetry-adapted VQE scheme with a shallower quantum circuit can achieve significant improvement in terms of the fidelity of the ground state and has a great advantage in terms of the ground-state energy with decent accuracy, as compared to the non-symmetry-adapted VQE scheme. We also demonstrate that the present scheme can approximate low-lying excited states that can be specified by symmetry sectors, using the same circuit structure for the ground-state calculation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Thermodynamic properties of an $S=1/2$ ring-exchange model on the triangular lattice

By using a numerically exact diagonalization technique and a block-extended version of the finite-temperature Lanczos method, we study thermodynamic properties of an $S=1/2$ Heisenberg model on the triangular lattice with an antiferromagnetic nearest-neighbor interaction $J$ and a four-spin ring-exchange interaction $J_{\rm c}$. Calculations are performed on small clusters under the periodic-boundary conditions. In contrast to the purely triangular case with $J_{\rm c}=0$, the specific heat exhibits a characteristic double-peak structure for $J_{\rm c}/J \gtrsim 0.04$. From the calculation of the entropy and the uniform magnetic susceptibility, it is shown that non-magnetic excitations exist below the magnetic excitation for $J_{\rm c}/J \gtrsim 0.04$.

0 citations0 reviewsNo code linkedOpen access