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

Seiji Yunoki

Seiji Yunoki appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.str-elcond-mat.mtrl-sciquant-phcond-mat.stat-mechcond-mat.supr-conhep-latphysics.comp-phcond-mat.mes-hallcond-mat.otherMachine LearningPerformancephysics.chem-ph

Catalog footprint

What is connected

64works
12topics
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

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

Brillouin-zone integration scheme for many-body density of states: Tetrahedron method combined with cluster perturbation theory

By combining the tetrahedron method with the cluster perturbation theory (CPT), we present an accurate method to numerically calculate the density of states of interacting fermions without introducing the Lorentzian broadening parameter $η$ or the numerical extrapolation of $η\to 0$. The method is conceptually based on the notion of the effective single-particle Hamiltonian which can be subtracted in the Lehmann representation of the single-particle Green's function within the CPT. Indeed, we show the general correspondence between the self-energy and the effective single-particle Hamiltonian which describes exactly the single-particle excitation energies of interacting fermions. The detailed formalism is provided for two-dimensional multi-orbital systems and a benchmark calculation is performed for the two-dimensional single-band Hubbard model. The method can be adapted straightforwardly to symmetry broken states, three-dimensional systems, and finite-temperature calculations.

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

Electron-hole doping asymmetry of Fermi surface reconstructed in a simple Mott insulator

It is widely recognised that the effect of doping into a Mott insulator is complicated and unpredictable, as can be seen by examining the Hall coefficient in high $T_{\rm c}$ cuprates. The doping effect, including the electron-hole doping asymmetry, may be more straightforward in doped organic Mott insulators owing to their simple electronic structures. Here we investigate the doping asymmetry of an organic Mott insulator by carrying out electric-double-layer transistor measurements and using cluster perturbation theory. The calculations predict that strongly anisotropic suppression of the spectral weight results in the Fermi arc state under hole doping, while a relatively uniform spectral weight results in the emergence of a non-interacting-like Fermi surface in the electron-doped state. In accordance with the calculations, the experimentally observed Hall coefficients and resistivity anisotropy correspond to the pocket formed by the Fermi arcs under hole doping and to the non-interacting Fermi surface under electron doping.

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

Nearly Free Electron States in MXenes

Using a set of first-principles calculations, we studied the electronic structures of two-dimensional transition metal carbides and nitrides, so called MXenes, functionalized with F, O, and OH. Our projected band structures and electron localization function analyses reveal the existence of nearly free electron (NFE) states in variety of MXenes. The NFE states are spatially located just outside the atomic structure of MXenes and are extended parallel to the surfaces. Moreover, we found that the OH-terminated MXenes offer the NFE states energetically close to the Fermi level. In particular, the NFE states in some of the OH-terminated MXenes, such as Ti2C(OH)2, Zr2C(OH)2, Zr2N(OH)2, Hf2C(OH)2, Hf2N(OH)2, Nb2C(OH)2, and Ta2C(OH)2, are partially occupied. This is in remarkable contrast to graphene, graphane, and MoS2, in which their NFE states are located far above the Fermi level and thus they are unoccupied. As a prototype of such systems, we investigated the electron transport properties of Hf2C(OH)2 and found that the NFE states in Hf2C(OH)2 provide almost perfect transmission channels without nuclear scattering for electron transport. Our results indicate that these systems might find applications in nanoelectronic devices. Our findings provide new insights into the unique electronic band structures of MXenes.

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
preprint2016arXiv

Topological insulators in ordered double transition metals M$'_2$M$"$C$_2$ (M$'$= Mo, W; M$"$= Ti, Zr, Hf) MXenes

The family of two-dimensional transition metal carbides, so called MXenes, has recently found new members with ordered double transition metals M$'_2$M$"$C$_2$, where M$'$ and M$"$ stand for transition metals. Here, using a set of first-principles calculations, we demonstrate that some of the newly added members, oxide M$'_2$M$"$C$_2$ (M$'$= Mo, W; M$"$= Ti, Zr, Hf) MXenes, are topological insulators. The nontrivial topological states of the predicted MXenes are revealed by the $Z_2$ index, which is evaluated from the parities of the occupied bands below the Fermi energy at time reversal invariant momenta, and also by the presence of the edge states. The predicted M$'_2$M$"$C$_2$O$_2$ MXenes show nontrivial gaps in the range of 0.041 -- 0.285 eV within the generalized gradient approximation and 0.119 -- 0.409 eV within the hybrid functional. The band gaps are induced by the spin-orbit coupling within the degenerate states with $d_{x^2-y^2}$ and $d_{xy}$ characters of M$'$ and M$"$, while the band inversion occurs at the $Γ$ point among the degenerate $d_{x^2-y^2}$/$d_{xy}$ orbitals and a non-degenerate $d_{3z^2-r^2}$ orbital, which is driven by the hybridization of the neighboring orbitals. The phonon dispersion calculations find that the predicted topological insulators are structurally stable. The predicted W-based MXenes with large band gaps might be suitable candidates for many topological applications at room temperature. In addition, we study the electronic structures of thicker ordered double transition metals M$'_2$M$"_2$C$_3$O$_2$ (M$'$= Mo, W; M$"$= Ti, Zr, Hf) and find that they are nontrivial topological semimetals.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Universal quantum criticality in the metal-insulator transition of two-dimensional interacting Dirac electrons

The metal-insulator transition has been a subject of intense research since Nevil Mott has first proposed that the metallic behavior of interacting electrons could turn to the insulating one as electron correlations increase. Here, we consider electrons with massless Dirac-like dispersion in two spatial dimensions, described by the Hubbard models on two geometrically different lattices, and perform numerically exact calculations on unprecedentedly large systems that, combined with a careful finite size scaling analysis, allow us to explore the quantum critical behavior in the vicinity of the interaction-driven metal-insulator transition. We find thereby that the transition is continuous and determine the quantum criticality for the corresponding universality class, which is described in the continuous limit by the Gross-Neveu model, a model extensively studied in quantum field theory. We furthermore discuss a fluctuation-driven scenario for the metal-insulator transition in the interacting Dirac electrons: the metal-insulator transition is triggered only by the vanishing of the quasiparticle weight but not the Dirac Fermi velocity, which instead remains finite near the transition. This important feature cannot be captured by a simple mean-field or Gutzwiller-type approximate picture, but is rather consistent with the low energy behavior of the Gross-Neveu model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A variational Monte Carlo study of exciton condensation

Exciton condensation in a two-band Hubbard model on a square lattice is studied with variational Monte Carlo method. We show that the phase transition from an excitonic insulator to a band insulator is induced by increasing the interband Coulomb interaction. To examine the character of the exciton condensation, the exciton pair amplitudes both in $\bm{k}$-space and in real space are calculated. Using these quantities, we discuss the BCS-BEC crossover within the excitonic insulator phase.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Charge-density wave induced by combined electron-electron and electron-phonon interactions in 1$T$-TiSe$_2$: A variational Monte Carlo study

To clarify the origin of a charge-density wave (CDW) phase in 1$T$-TiSe$_2$, we study the ground state property of a half-filled two-band Hubbard model in a triangular lattice including electron-phonon interaction. By using the variational Monte Carlo method, the electronic and lattice degrees of freedom are both treated quantum mechanically on an equal footing beyond the mean-field approximation. We find that the cooperation between Coulomb interaction and electron-phonon interaction is essential to induce the CDW phase. We show that the "pure" exciton condensation without lattice distortion is difficult to realize under the poor nesting condition of the underlying Fermi surface. Furthermore, by systematically calculating the momentum resolved hybridization between the two bands, we examine the character of electron-hole pairing from the viewpoint of BCS-BEC crossover within the CDW phase and find that the strong-coupling BEC-like pairing dominates. We therefore propose that the CDW phase observed in 1$T$-TiSe$_2$ originates from a BEC-like electron-hole pairing.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Enhanced charge excitations in electron-doped cuprates by resonant inelastic x-ray scattering

Resonant inelastic x-ray scattering (RIXS) tuned for the Cu L edge is a possible tool to detect charge excitations in cuprate superconductors. We theoretically investigate the possibility for observing a collective charge excitation by the RIXS. The RIXS process via the intermediate state inevitably makes the spectral weight of charge excitation stronger in electron doping than in hole doping. Electron-hole asymmetry also appears in the dynamical charge structure factor, showing a new enhanced small-momentum low-energy mode in electron doping. These facts indicate a possibility of detecting the new charge mode by RIXS in electron-doped systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Giant switchable Rashba effect in oxide heterostructures

One of the most fundamental phenomena and a reminder of the electron's relativistic nature is the Rashba spin splitting for broken inversion symmetry. Usually this splitting is a tiny relativistic correction, hardly discernible in experiment. Interfacing a ferroelectric, BaTiO$_3$, and a heavy 5$d$ metal with a large spin-orbit coupling, Ba(Os,Ir)O$_3$, we show that giant Rashba spin splittings are indeed possible and even fully controllable by an external electric field. Based on density functional theory and a microscopic tight binding understanding, we conclude that the electric field is amplified and stored as a ferroelectric Ti-O distortion which, through the network of oxygen octahedra, also induces a large Os-O distortion. The BaTiO$_3$/BaOsO$_3$ heterostructure is hence the ideal test station for studying the fundamentals of the Rashba effect. It allows intriguing application such as the Datta-Das transistor to operate at room temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Large-Gap Two-Dimensional Topological Insulator in Oxygen Functionalized MXene

Two-dimensional (2D) topological insulator (TI) have been recognized as a new class of quantum state of matter. They are distinguished from normal 2D insulators with their nontrivial band-structure topology identified by the $Z_2$ number as protected by time-reversal symmetry (TRS). 2D TIs have intriguing spin-velocity locked conducting edge states and insulating properties in the bulk. In the edge states, the electrons with opposite spins propagate in opposite directions and the backscattering is fully prohibited when the TRS is conserved. This leads to quantized dissipationless "two-lane highway" for charge and spin transportation and promises potential applications. Up to now, only very few 2D systems have been discovered to possess this property. The lack of suitable material obstructs the further study and application. Here, by using first-principles calculations, we propose that the functionalized MXene with oxygen, M$_2$CO$_2$ (M=W, Mo and Cr), are 2D TIs with the largest gap of 0.194 eV in W case. They are dynamically stable and natively antioxidant. Most importantly, they are very likely to be easily synthesized by recent developed selective chemical etching of transition-metal carbides (MAX phase). This will pave the way to tremendous applications of 2D TIs, such as "ideal" conducting wire, multifunctional spintronic device, and the realization of topological superconductivity and Majorana modes for quantum computing.

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

0-$π$ transition driven by magnetic proximity effect in a Josephson junction

We theoretically study the Josephson effect in a superconductor/normal metal/superconductor ({\it S}/{\it N}/{\it S}) Josephson junction composed of $s$-wave {\it S}s with {\it N} which is sandwiched by two ferromagnetic insulators ({\it F}s), forming a spin valve, in the vertical direction of the junction. We show that the 0-$π$ transition of the Josephson critical current occurs with increasing the thickness of {\it N} along the junction. This transition is due to the magnetic proximity effect (MPE) which induces ferromagnetic magnetization in the {\it N}. Moreover, we find that, even for fixed thickness of {\it N}, the proposed Josephson junction with the spin valve can be switched from $π$ to 0 states and vice versa by varying the magnetization configuration (parallel or antiparallel) of two {\it F}s. We also examine the effect of spin-orbit scattering on the Josephson critical current and argue that the 0-$π$ transition found here can be experimentally observed within the current nanofabrication techniques, thus indicating a promising potential of this junction as a 0-$π$ switching device operated reversibly with varying the magnetic configuration in the spin valve by, e.g., applying an external magnetic field. %with the magnetization configuration in the spin valve. Our results not only provide possible applications in superconducting electronics but also suggest the importance of a fundamental concept of MPE in nanostructures of multilayer {\it N}/{\it F} systems.

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

Density-matrix renormalization group study of third harmonic generation in one-dimensional Mott insulator coupled with phonon

We examine the third-order non-linear optical response of a one-dimensional Mott insulator coupled with phonons. The Mott insulator is described by an extended Hubbard-Holstein model. The third harmonic generation (THG) of the model is calculated by dynamical density-matrix renormalization group. We find that the electron-phonon interaction enhances the intensity of THG if the on-site Coulomb interaction is comparable to the band width, but the enhancement is small for a realistic parameter value of a typical Mott insulator Sr$_2$CuO$_3$. The low-energy spin excitation does not contribute to the optical response if the electron-phonon coupling is neglected. We find that the THG spectrum can detect the spin excitations even for no electron-phonon coupling. The introduction of the electron-phonon coupling leads to a slight increase of low-energy weight in THG without peak structure through phonon-assisted spin excitation. To fully understand the phonon-assisted spin excitation, we additionally calculate the linear susceptibility of a doped Hubbard-Holstein model, showing a spectral distribution different from that at half filling.

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
preprint2014arXiv

Unexpectedly high pressure for molecular dissociation in liquid hydrogen by a reliable electronic simulation

The study of the high pressure phase diagram of hydrogen has continued with renewed effort for about one century as it remains a fundamental challenge for experimental and theoretical techniques. Here we employ an efficient molecular dynamics based on the quantum Monte Carlo method, which can describe accurately the electronic correlation and treat a large number of hydrogen atoms, allowing a realistic and reliable prediction of thermodynamic roperties. We find that the molecular liquid phase is unexpectedly stable and the transition towards a fully atomic liquid phase occurs at much higher pressure than previously believed. The old standing problem of low temperature atomization is, therefore, still far from experimental reach.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Study of Parallelizing O(N) Green-Function-Based Monte Carlo Method for Many Fermions Coupled with Classical Degrees of Freedom

Models of fermions interacting with classical degrees of freedom are applied to a large variety of systems in condensed matter physics. For this class of models, Weiße [Phys. Rev. Lett. {\bf 102}, 150604 (2009)] has recently proposed a very efficient numerical method, called O($N$) Green-Function-Based Monte Carlo (GFMC) method, where a kernel polynomial expansion technique is used to avoid the full numerical diagonalization of the fermion Hamiltonian matrix of size $N$, which usually costs O($N^3$) computational complexity. Motivated by this background, in this paper we apply the GFMC method to the double exchange model in three spatial dimensions. We mainly focus on the implementation of GFMC method using both MPI on a CPU-based cluster and Nvidia's Compute Unified Device Architecture (CUDA) programming techniques on a GPU-based (Graphics Processing Unit based) cluster. The time complexity of the algorithm and the parallel implementation details on the clusters are discussed. We also show the performance scaling for increasing Hamiltonian matrix size and increasing number of nodes, respectively. The performance evaluation indicates that for a $32^3$ Hamiltonian a single GPU shows higher performance equivalent to more than 30 CPU cores parallelized using MPI.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamical spin structure factors of quantum spin nematic states

Dynamical spin structure factors of quantum spin nematic states are calculated in a spin-1/2 square-lattice J1-J2 model with ferromagnetic J1 and competing antiferromagnetic J2 interactions. To this end, we use a fermion representation, generalizing it to N flavors. We begin with a spin-triplet pairing state of fermion fields, called Z2 planar state, which is a stable saddle-point solution in the large-N limit in a finite parameter range where the couplings J1 and J2 compete strongly [R. Shindou and T. Momoi, Phys. Rev. B 80, 064410 (2009)]. Using a large-N expansion, we take into account fluctuations around this saddle point up to corrections of order 1/N. The dynamical spin structure factors thus obtained signify the existence of gapless q-linear director-wave (spin-wave) modes at q=(0,0) and gapped `gauge-field' like collective modes at q=(pi,pi), whose spectral weight vanishes as a linear and quadratic function of the momentum respectively. The low-energy collective modes contain fluctuations of nematic-director, spin, and gauge degrees of freedom. Associated with the gapless q-linear modes, we evaluate the temperature dependence of the nuclear spin relaxation rate 1/T1 in the low-temperature regime as 1/T1 \propto T^{2d-1}, where d is the effective spatial dimension.

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

Absence of a Spin Liquid Phase in the Hubbard Model on the Honeycomb Lattice

A spin liquid is a novel quantum state of matter with no conventional order parameter where a finite charge gap exists even though the band theory would predict metallic behavior. Finding a stable spin liquid in two or higher spatial dimensions is one of the most challenging and debated issues in condensed matter physics. Very recently, it has been reported that a model of graphene, i.e., the Hubbard model on the honeycomb lattice, can show a spin liquid ground state in a wide region of the phase diagram, between a semi-metal (SM) and an antiferromagnetic insulator (AFMI). Here, by performing numerically exact quantum Monte Carlo simulations, we extend the previous study to much larger clusters (containing up to 2592 sites), and find, if any, a very weak evidence of this spin liquid region. Instead, our calculations strongly indicate a direct and continuous quantum phase transition between SM and AFMI.

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

Magnetic and Orbital Order in ($R$MnO$_3$)$_n$/($A$MnO$_3$)$_{2n}$ Superlattices Studied via a Double-Exchange Model with Strain

The two-orbital double-exchange model is employed for the study of the magnetic and orbital orders in ($R$MnO$_3$)$_n$/($A$MnO$_3$)$_{2n}$ ($R$: rare earths; $A$: alkaline earths) superlattices. The A-type antiferromagnetic order is observed in a broad region of parameter space for the case of SrTiO$_3$ as substrate, in agreement with recent experiments and first-principles calculations using these superlattices. In addition, also a C-type antiferromagnetic state is predicted to be stabilized when using substrates like LaAlO$_3$ with smaller lattice constants than SrTiO$_3$, again in agreement with first principles results. The physical mechanism for the stabilization of the A- and C- magnetic transitions is driven by the orbital splitting of the $x^2-y^2$ and $3z^2-r^2$ orbitals. This splitting is induced by the $Q_3$ mode of Jahn-Teller distortions created by the strain induced by the substrates. In addition to the special example of (LaMnO$_3$)$_n$/(SrMnO$_3$)$_{2n}$, our phase diagrams can be valuable for the case where the superlattices are prepared employing narrow bandwidth manganites. In particular, several non-homogenous magnetic profiles are predicted to occur in narrow bandwidth superlattices, highlighting the importance of carrying out investigations in this mostly unexplored area of research.

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
preprint2012arXiv

Strain-engineered magnetic order in (LaMnO$_{3}$)$_n$/(SrMnO$_{3}$)$_{2n}$ superlattices

Using first-principles calculations based on the density functional theory, we show a strong strain dependence of magnetic order in (LaMnO$_{3}$)$_n$/(SrMnO$_{3}$)$_{2n}$ (001) superlattices with $n=1,2$. The epitaxial strain lifts the degeneracy of Mn $e_{g}$ orbitals, thus inducing an inherent orbital order, which in turn strongly affects the ferromagnetic double exchange of itinerant $e_{g}$ electrons, competing with the antiferromagnetic superexchange of localized $t_{2g}$ electrons. For the case of tensile strain induced by SrTiO$_3$ (001) substrate, we find that the ground state is A-type antiferromagnetic and $d_{x^2-y^2}$ orbital ordered, which is in excellent agreement with recent experiments [S. J. May {\it et al.}, Nature Materials {\bf 8}, 892 (2009)]. Instead, for the case of compressive strain induced by LaAlO$_3$ (001) substrate, we predict that the ground state is C-type antiferromagnetic and $d_{3z^2-r^2}$ orbital ordered.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Ab Initio Study of the Intrinsic Exchange Bias at the SrRuO$_3$/SrMnO$_3$ Interface

In a recent publication (S. Dong et al., Phys. Rev. Lett.103, 127201 (2009)), two (related) mechanisms were proposed to understand the intrinsic exchange bias present in oxides heterostructures involving G-type antiferromagnetic perovskites. The first mechanism is driven by the Dzyaloshinskii-Moriya interaction, which is a spin-orbit coupling effect. The second is induced by the ferroelectric polarization, and it is only active in heterostructures involving multiferroics. Using the SrRuO$_3$/SrMnO$_3$ superlattice as a model system, density-functional calculations are here performed to verify the two proposals. This proof-of-principle calculation provides convincing evidence that qualitatively supports both proposals.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

First-principles study of the spin-mixing conductance in Pt/Ni$_{81}$Fe$_{19}$ junctions

Based on the spin-pumping theory and first-principles calculations, the spin-mixing conductance (SMC) is theoretically studied for Pt/Permalloy (Ni$_{81}$Fe$_{19}$, Py) junctions. We evaluate the SMC for ideally clean Pt/Py junctions and examine the effects of interface randomness. We find that the SMC is generally enhanced in the presence of interface roughness as compared to the ideally clean junctions. Our estimated SMC is in good quantitative agreement with the recent experiment for Pt/Py junctions. We propose possible routes to increase the SMC in Pt/Py junctions by depositing a foreign magnetic metal layer in Pt, offering guidelines for designing the future spintronic devices.

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
preprint2011arXiv

Multiferroic Properties of CaMn$_7$O$_{12}$

We report that CaMn$_7$O$_{12}$ is a new magnetic multiferroic material. The appearance of a ferroelectric polarization coinciding with the magnetic phase transition ($\sim90$ K) suggests the presence of ferroelectricity induced by magnetism, further confirmed by its strong magnetoelectric response. With respect to other known magnetic multiferroics, CaMn$_7$O$_{12}$ displays attractive multiferroic properties, such as a high ferroelectric critical temperature and large polarization. More importantly, these results open a new avenue to search for magnetic multiferroics in the catalogue of doped oxides.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Performance Acceleration of Kernel Polynomial Method Applying Graphics Processing Units

The Kernel Polynomial Method (KPM) is one of the fast diagonalization methods used for simulations of quantum systems in research fields of condensed matter physics and chemistry. The algorithm has a difficulty to be parallelized on a cluster computer or a supercomputer due to the fine-gain recursive calculations. This paper proposes an implementation of the KPM on the recent graphics processing units (GPU) where the recursive calculations are able to be parallelized in the massively parallel environment. This paper also illustrates performance evaluations regarding the cases when the actual simulation parameters are applied, the one for increased intensive calculations and the one for increased amount of memory usage. Finally, it concludes that the performance on GPU promises very high performance compared to the one on CPU and reduces the overall simulation time.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Projective studies of spin nematics in a quantum frustrated ferromagnet

We study the ground state properties of the spin-1/2 frustrated ferromagnetic J1-J2 Heisenberg model on the square lattice, employing projected BCS wavefunctions with spin-triplet pairings of the spinon fields as trial wavefunctions. Based on the variational Monte Carlo analysis, we argue that, in the competing coupling regime, a certain type of the projected BCS wavefunction, dubbed the projected Z2 planar state, achieves the best optimal energy among the other competing states such as the ferromagnetic state and collinear antiferromagnetic state. Like in quantum spin liquids, the projected Z2 planar state preserves the translational symmetry of the square lattice. However, it is also accompanied by a d-wave ordering of the quadrupole moments, breaking the spin rotational symmetry. The state thus describes a quantum spin analogue of the nematic liquid crystals. The calculated static correlation functions also reveal that the projected Z2 planar state has a strong collinear antiferromagnetic fluctuation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Strain driven onset of non-trivial topological insulating states in Zintl Sr$_2X$ compounds ($X$=Pb, Sn)

We explore the topological behavior of the binary Zintl phase of the alkaline earth metals based compounds Sr$_2$Pb and Sr$_2$Sn using both standard and hybrid density functional theory. It is found that Sr$_2$Pb lies on the verge of a topological instability which can be suitably tuned through the application of a small uniaxial expansion strain ($>$ 3%). The resulting non-trivial topologically insulating state display well-defined metallic states in the Sr$_{2}$Pb(010) surface, whose evolution is studied as a function of the film thickness.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Highly anisotropic resistivities in the double-exchange model for strained manganites

The highly anisotropic resistivities in strained manganites are theoretically studied using the two-orbital double-exchange model. At the nanoscale, the anisotropic double-exchange and Jahn-Teller distortions are found to be responsible for the robust anisotropic resistivities observed here via Monte Carlo simulations. An unbalanced in the population of orbitals caused by strain is responsible for these effects. In contrast, the anisotropic superexchange is found to be irrelevant to explain our results. Our model study suggests that highly anisotropic resistivities could be present in a wide range of strained manganites, even without (sub)micrometer-scale phase separation. In addition, our calculations also confirm the formation of anisotropic clusters in phase-separated manganites, which magnifies the anisotropic resistivities.

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
preprint2009arXiv

Double-exchange model study of multiferroic $R$MnO$_3$ perovskites

In this proceeding, recent theoretical investigations by the authors on the multiferroic $R$MnO$_3$ perovskites are briefly reviewed at first. Using the double-exchange model, the realistic spiral spin order in undoped manganites such as TbMnO$_3$ and DyMnO$_3$ is well reproduced by incorporating a weak next-nearest neighbor superexchange ($\sim10%$ of nearest neighbor superexchange) and moderate Jahn-Teller distortion. The phase transitions from the A-type antiferromagnet (as in LaMnO$_3$), to the spiral phase (as in TbMnO$_3$), and finally to the E-type antiferromagnet (as in HoMnO$_3$), with decreasing size of the $R$ ions, were also explained. Moreover, new results of phase diagram of the three-dimensional lattice are also included. The ferromagnetic tendency recently discovered in the LaMnO$_3$ and TbMnO$_3$ thin films is explained by considering the substrate stress. Finally, the relationship between our double-exchange model and a previously used $J_1$-$J_2$-$J_3$ model is further discussed from the perspective of spin wave excitations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Electronic and Magnetic Properties of $\rm{RMnO}_3/\rm{AMnO}_3$ Heterostructures

In this paper, we investigate the ground state properties of the $\rm{RMnO}_3/\rm{AMnO}_3$ ($\rm{RMO/AMO}$) heterostructures (R=trivalent cation, and A=divalent cation) by using a two-orbital double-exchange model supplemented by the Poisson's equation. We find that the state stabilized near the interface of the heterostructure is similar to the state of the bulk compound $\rm{(R,A)MO}$ at electronic density close to 0.5. Depending on the bandwidth it will be a charge and orbital ordered CE state or an A-AF state with $x^2-y^2$ orbital order. Our results can explain some properties of long-period superlattices. As another interesting result, we find exotic intermediate states stabilized in between the interface and the bulk-like regimes of the heterostructure. For instance, a spin "canted CE" state and others. They may not have an analog in experimentally known bulk phase diagrams, but provide a natural interpolation between magnetically-ordered states that are stable in the bulk at different electronic densities.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

One-Particle Anomalous Excitations of Gutzwiller Projected BCS Superconductors and Bogoliubov Quasi-Particle Characteristics

Low-lying one-particle anomalous excitations are studied for Gutzwiller projected strongly correlated BCS states. It is found that the one-particle anomalous excitations are highly coherent, and the numerically calculated spectrum can be reproduced quantitatively by a renormalized BCS theory, thus strongly indicating that the nature of low-lying excitations described by the projected BCS states is essentially understood within a renormalized Bogoliubov quasi-particle picture. This finding resembles a well known fact that a Gutzwiller projected Fermi gas is a Fermi liquid. The present results are consistent with numerically exact calculations of the two-dimensional t-J model as well as recent photoemission experiments on high-T_{\rm C} cuprate superconductors.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Two spin liquid phases in the spatially anisotropic triangular Heisenberg model

The quantum spin-1/2 antiferromagnetic Heisenberg model on a two dimensional triangular lattice geometry with spatial anisotropy is relevant to describe materials like ${\rm Cs_2 Cu Cl_4}$ and organic compounds like {$κ$-(ET)$_2$Cu$_2$(CN)$_3$}. The strength of the spatial anisotropy can increase quantum fluctuations and can destabilize the magnetically ordered state leading to non conventional spin liquid phases. In order to understand these intriguing phenomena, quantum Monte Carlo methods are used to study this model system as a function of the anisotropic strength, represented by the ratio $J'/J$ between the intra-chain nearest neighbor coupling $J$ and the inter-chain one $J'$. We have found evidence of two spin liquid regions. The first one is stable for small values of the coupling $J'/J \alt 0.65$, and appears gapless and fractionalized, whereas the second one is a more conventional spin liquid with a small spin gap and is energetically favored in the region $0.65\alt J'/J \alt 0.8$. We have also shown that in both spin liquid phases there is no evidence of broken translation symmetry with dimer or spin-Peirls order or any broken spatial reflection symmetry of the lattice. The various phases are in good agreement with the experimental findings, thus supporting the existence of spin liquid phases in two dimensional quantum spin-1/2 systems.

0 citations0 reviewsNo code linkedOpen access