BZPEERReading, trust and collaboration for research
Menu

Discover

SearchFind papers, people, topics and opportunities from one place.GraphOpen maps around papers, topics and saved views.

Workspaces

HomeOpen the ranked research feed.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.

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

Tomonori Shirakawa

Tomonori Shirakawa contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.str-elquant-phcond-mat.mtrl-scicond-mat.supr-concond-mat.stat-mechMachine Learningphysics.chem-phphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - Emerging
19works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0
Open graphBrowse works

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 map preview

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

Published work

19 published item(s)

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

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
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
preprint2020arXiv

General correlated geminal ansatz for electronic structure calculations: exploiting Pfaffians in place of determinants

We propose here a single Pfaffian correlated variational ansatz, that dramatically improves the accuracy with respect to the single determinant one, while remaining at a similar computational cost. A much larger correlation energy is indeed determined by the most general two electron pairing function, including both singlet and triplet channels, combined with a many-body Jastrow factor, including all possible spin-spin spin-density and density-density terms. The main technical ingredient to exploit this accuracy is the use of the Pfaffian for antisymmetrizing an highly correlated pairing function, thus recovering the Fermi statistics for electrons with an affordable computational cost. Moreover the application of the Diffusion Monte Carlo, within the fixed node approximation, allows us to obtain very accurate binding energies for the first preliminary calculations reported in this study: C$_2$, N$_2$ and O$_2$ and the benzene molecule. This is promising and remarkable, considering that they represent extremely difficult molecules even for computationally demanding multi-determinant approaches, and opens therefore the way for realistic and accurate electronic simulations with an algorithm scaling at most as the fourth power of the number of electrons.

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

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
preprint2016arXiv

Density matrix renormalization group study in energy space for a single-impurity Anderson model and an impurity quantum phase transition

The density matrix renormalization group method is applied to obtain the ground state phase diagram of the single impurity Anderson model on the honeycomb lattice at half filling. The calculation of local static quantities shows that the phase diagram contains two distinct phases, the local moment (LM) phase and the asymmetric strong coupling (ASC) phase. These results are supported by the local spin and charge excitation spectra, which exhibit qualitatively different behavior in these two phases and also reveal the existence of the valence fluctuating point at the phase boundary. For comparison, we also study the low-energy effective pseudogap Anderson model. Although the high-energy excitations are obviously different, we find that the ground state phase diagram and the asymptotically low-energy excitations are in good quantitative agreement with those for the single impurity Anderson model on the honeycomb lattice, thus providing the first quantitative justification for the previous studies based on low-energy approximate approaches. Furthermore, we find that the lowest entanglement level is doubly degenerate for the LM phase, whereas it is singlet for the ASC phase and is accidentally three fold degenerate at the valence fluctuating point. Our results therefore clearly demonstrate that the low-lying entanglement spectrum can be used to determine with high accuracy the phase boundary of the impurity quantum phase transition.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Emergence of massless Dirac quasiparticles in correlated hydrogenated graphene with broken sublattice symmetry

Using the variational cluster approximation (VCA) and the cluster perturbation theory, we study the finite temperature phase diagram of a half-depleted periodic Anderson model on the honeycomb lattice at half filling for a model of graphone, i.e., single-side hydrogenated graphene. The ground state of this model is found to be ferromagnetic (FM) semi-metal. The spin wave dispersion in the FM state is linear in momentum at zero temperature but becomes quadratic at finite temperatures, implying that the FM state is fragile against thermal fluctuations. Indeed, our VCA calculations find that the paramagnetic (PM) state dominates the finite temperature phase diagram. More surprisingly, we find that massless Dirac quasiparticles with the linear energy dispersion emerge at the Fermi energy upon introducing the electron correlation $U$ at the impurity sites in the PM phase. The Dirac Fermi velocity is found to be highly correlated to the quasiparticle weight of the emergent massless Dirac quasiparticles at the Fermi energy and monotonically increases with $U$. These unexpected massless Dirac quasiparticles are also examined with the Hubbard-I approximation and the origin is discussed in terms of the spectral weight redistribution involving a large energy scale of $U$. Considering an effective quasiparticle Hamiltonian which reproduces the single-particle excitations obtained by the Hubbard-I approximation, we argue that the massless Dirac quasiparticles are protected by the electron correlation. Our finding is therefore the first example of the emergence of massless Dirac quasiparticles due to dynamical electron correlations without breaking any spatial symmetry. The experimental implications of our results for graphone as well as a graphene sheet on transition metal substrates are also briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

From a quasimolecular band insulator to a relativistic Mott insulator in $t_{2g}^5$ systems with a honeycomb lattice structure

The $t_{2g}$ orbitals of an edge-shared transition-metal oxide with a honeycomb lattice structure form dispersionless electronic bands when only hopping mediated by the edge-sharing oxygens is accessible. This is due to the formation of isolated quasimolecular orbitals (QMOs) in each hexagon, introduced recently by Mazin et al. [Phys. Rev. Lett. 109, 197201 (2012)], which stabilizes a band insulating phase for $t_{2g}^5$ systems. However, with help of the exact diagonalization method to treat the electron kinetics and correlations on an equal footing, we find that the QMOs are fragile against not only the spin-orbit coupling (SOC) but also the Coulomb repulsion. We show that the electronic phase of $t_{2g}^5$ systems can vary from a quasimolecular band insulator to a relativistic $J_{\rm eff}=1/2$ Mott insulator with increasing the SOC as well as the Coulomb repulsion. The different electronic phases manifest themselves in electronic excitations observed in optical conductivity and resonant inelastic x-ray scattering. Based on our calculations, we assert that the currently known Ru$^{3+}$- and Ir$^{4+}$-based honeycomb systems are far from the quasimolecular band insulator but rather the relativistic Mott insulator.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Slave rotor theory of the Mott transition in the Hubbard model: a new mean field theory and a new variational wave function

A new mean field theory is proposed for the Mott transition in the Hubbard model based on the slave rotor representation of the electron operator. This theory provides a better description of the role of the long range charge correlation in the Mott insulating state and offers a good estimation of the critical correlation strength for the Mott transition. We have constructed a new variational wave function for the Mott insulating state based on this new slave rotor mean field theory. We find this new variational wave function outperforms the conventional Jastrow type wave function with long range charge correlator in the Mott insulating state. It predicts a continuous Mott transition with non-divergent quasiparticle mass at the transition point. We also show that the commonly used on-site mean field decoupling for the slave rotor corresponds to the Gutzwiller approximation for the Gutzwiller projected wave function with only on-site charge correlator, which can not describe the Mott transition in any finite dimensional system.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spin-orbit-induced exotic insulators in a three-orbital Hubbard model with $(t_{2g})^5$ electrons

On the basis of the multi-orbital dynamical mean field theory, a three-orbital Hubbard model with a relativistic spin-orbit coupling (SOC) is studied at five electrons per site. The numerical calculations are performed by employing the continuous-time quantum Monte Carlo (CTQMC) method based on the strong coupling expansion. We find that appropriately choosing bases, i.e., the maximally spin-orbit-entangled bases, drastically improve the sign problem in the CTQMC calculations, which enables us to treat exactly the full Hund's coupling and pair hopping terms. This improvement is also essential to reach at low temperatures for a large SOC region where the SOC most significantly affects the electronic structure. We show that a metal-insulator transition is induced by the SOC for fixed Coulomb interactions. The insulating state for smaller Coulomb interactions is antiferromagnetically ordered with the local effective total angular momentum $j=1/2$, in which the $j=1/2$ based band is essentially half-filled while the $j=3/2$ based bands are completely occupied. More interestingly, for larger Coulomb interactions, we find that an excitonic insulating state emerges, where the condensation of an electron-hole pair in the $j=1/2$ and $j=3/2$ based bands occurs. The origin of the excitonic insulator as well as the experimental implication is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Block Lanczos density-matrix renormalization group method for general Anderson impurity models: Application to magnetic impurity problems in graphene

We introduce a block Lanczos (BL) recursive technique to construct quasi-one-dimensional models, suitable for density-matrix renormalization group (DMRG) calculations, from single- as well as multiple-impurity Anderson models in any spatial dimensions. This new scheme, named BL-DMRG method, allows us to calculate not only local but also spatially dependent static and dynamical quantities of the ground state for general Anderson impurity models without losing elaborate geometrical information of the lattice. We show that the BL-DMRG method can be easily extended to treat a multi-orbital Anderson impurity model. We also show that the symmetry adapted BL bases can be utilized, when it is appropriate, to reduce the computational cost. As a demonstration, we apply the BL-DMRG method to three different models for graphene: (i) a single adatom on the honeycomb lattice, (ii) a substitutional impurity in the honeycomb lattice, and (iii) an effective model for a single carbon vacancy in graphene. Our analysis reveals that, for the particle-hole symmetric case at half filling of electron density, the ground state of model (i) behaves as an isolated magnetic impurity with no Kondo screening while the ground state of the other two models forms a spin singlet state. We also calculate the real-space dependence of the spin-spin correlation functions between the impurity site and the conduction sites for these three models. Our results clearly show that, reflecting the presence of absence of unscreened magnetic moment at the impurity site, the spin-spin correlation functions decay as $r^{-3}$, differently from the non-interacting limit ($r^{-2}$), for model (i) and as $ r^{-4}$, exactly the same as the non-interacting limit, for models (ii) and (iii) in the asymptotic $r$, where $r$ is the distance between the impurity site and the conduction site.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Theoretical study of insulating mechanism in multi-orbital Hubbard models with a large spin-orbit coupling: Slater versus Mott scenario in Sr$_2$IrO$_4$

To examine the insulating mechanism of 5$d$ transition metal oxide Sr$_2$IrO$_4$, we study the ground state properties of a three-orbital Hubbard model with a large relativistic spin-orbit coupling on a square lattice. Using a variational Monte Carlo method, we find that the insulating state appearing in the ground state phase diagram for one hole per site varies from a weakly correlated to a strongly correlated antiferromagnetic (AF) state with increasing Coulomb interactions. This crossover is characterized by the different energy gain mechanisms of the AF insulating state, i.e., from an interaction-energy driven Slater-type insulator to a band-energy driven Mott-type insulator with increasing Coulomb interactions. Our calculations reveal that Sr$_2$IrO$_4$ is a "moderately correlated" AF insulator located in the intermediate coupling region between a Slater-type and a Mott-type insulators.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Rock-salt SnS and SnSe: Native Topological Crystalline Insulators

Unlike time-reversal topological insulators, surface metallic states with Dirac cone dispersion in the recently discovered topological crystalline insulators (TCIs) are protected by crystal symmetry. To date, TCI behaviors have been observed in SnTe and the related alloys Pb$_{1-x}$Sn$_{x}$Se/Te, which incorporate heavy elements with large spin-orbit coupling (SOC). Here, by combining first-principles and {\it ab initio} tight-binding calculations, we report the formation of a TCI in the relatively lighter rock-salt SnS and SnSe. This TCI is characterized by an even number of Dirac cones at the high-symmetry (001), (110) and (111) surfaces, which are protected by the reflection symmetry with respect to the ($\bar{1}$10) mirror plane. We find that both SnS and SnSe have an intrinsically inverted band structure and the SOC is necessary only to open the bulk band gap. The bulk band gap evolution upon volume expansion reveals a topological transition from an ambient pressure TCI to a topologically trivial insulator. Our results indicate that the SOC alone is not sufficient to drive the topological transition.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Large spectral weight transfer in optical conductivity of SrTiO$_{3}$ induced by intrinsic vacancies

The optical conductivity ($σ_{1}$) of SrTiO$_{3}$ for various vacancies has been systematically studied using a combination of ultraviolet - vacuum ultraviolet (UV-VUV) reflectivity and spectroscopic ellipsometry. For cation (Ti) vacancies, $σ_{1}$ shows large spectral weight transfer over a wide range of energy from as high as 35 eV to as low as 0.5 eV and the presence of mid-gap states, suggesting that strong correlations play an important role. While for anion (O) vacancies $σ_{1}$ shows changes from 7.4 eV up to 35 eV. These unexpected results can be explained in terms of orbital reconstruction.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Monte Carlo study of an unconventional superconducting phase in Ir-oxide J_{eff}=1/2 Mott insulators induced by carrier doping

Based on a microscopic theoretical study, we show that novel superconductivity is induced by carrier doping in layered perovskite Ir oxides where a strong spin-orbit coupling causes an effective total angular momentum J_{eff}=1/2 Mott insulator. Using a variational Monte Carlo method, we find an unconventional superconducting state in the ground state phase diagram of a t_{2g} three-orbital Hubbard model on the square lattice. This superconducting state is characterized by a d_{x^2-y^2}-wave "pseudospin singlet" formed by the J_{eff}=1/2 Kramers doublet, which thus contains inter-orbital as well as both singlet and triplet components of t_{2g} electrons. The superconducting state is found stable only by electron doping, but not by hole doping, for the case of carrier doped Sr_2IrO_4. We also study an effective single-orbital Hubbard model to discuss the similarities to high-T_c cuprate superconductors and the multi-orbital effects.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Microscopic Study of Electronic and Magnetic Properties for Ir Oxide

The exact diagonalization and the variational cluster approximation (VCA) are used to study the nature of a novel Mott insulator induced by a strong spin-orbit coupling for a two-dimensional three-band Hubbard model consisting of the $t_{2g}$ manifold of $5d$ orbitals. To characterize the ground state, we introduce a local Kramer's doublet which can represent a state with effective angular momentum $J_{\rm eff}=|{\bm S}-{\bm L}|=1/2$ as well as spin $S=1/2$. Our systematic study of the pseudo-spin structure factor defined by the Kramer's doublet shows that the $J_{\rm eff}=1/2$ Mott insulator is smoothly connected to the $S=1/2$ Mott insulator. Using the Kramer's doublet as a variational state for the VCA, we examine the one-particle excitations for the Mott insulating phase. These results are compared with recent experiments on Sr$_2$IrO$_4$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Microscopic study of spin-orbit-induced Mott insulator in Ir oxides

Motivated by recent experiments of a novel 5$d$ Mott insulator in Sr$_2$IrO$_4$, we have studied the two-dimensional three-orbital Hubbard model with a spin-orbit coupling $λ$. The variational Monte Carlo method is used to obtain the ground state phase diagram with varying a on-site Coulomb interaction $U$ as well as $λ$. It is found that the transition from a paramagnetic metal to an antiferromagnetic (AF) insulator occurs at a finite $U=U_{\mathrm{MI}}$, which is greatly reduced by a large $λ$, characteristic of 5$d$ electrons, and leads to the "spin-orbit-induced" Mott insulator. It is also found that the Hund's coupling induces the anisotropic spin exchange and stabilizes the in-plane AF order. We have further studied the one-particle excitations using the variational cluster approximation, and revealed the internal electronic structure of this novel Mott insulator. These findings are in agreement with experimental observations on Sr$_2$IrO$_4$, and qualitatively different from those of extensively studied 3$d$ Mott insulators.

0 citations0 reviewsNo code linkedOpen access