Source author record

Masatoshi Imada

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

63works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

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

63 published item(s)

preprint2022arXiv

$Ab$ $initio$ material design of Ag-based oxides for high-$T_c$ superconductor

We propose silver-based oxides with layered perovskite structure as candidates of exhibiting intriguing feature of strongly correlated electrons. The compounds show unique covalence between Ag d and O p orbitals with the strong electron correlation of their antibonding orbital similar to the copper oxide high-temperature superconductors, but stronger covalency and slightly smaller correlation strength. We examine $A_2$AgO$_2X_2$ with $A$ =Sr and Ba and $X$ =F, I and Cl in detail. Among them, Sr$_2$AgO$_2$F$_2$ has the largest effective onsite Coulomb repulsion and shows an antiferromagnetic insulating ground state, which competes with correlated metals. It offers features both similar and distinct from the copper oxides, and paves a new route. The possibility of superconductivity in doped systems is discussed.

preprint2022arXiv

Optimized Implementation for Calculation and Fast-Update of Pfaffians Installed to the Open-Source Fermionic Variational Solver mVMC

In this article, we present a high performance, portable and well templated implementation for computing and fast-updating Pfaffian and inverse of an even-ranked skew-symmetric (antisymmetric) matrix. It is achieved with a skew-symmetric, blocked variant of the Parlett-Reid algorithm and a blocked update scheme based on the Woodbury matrix identity. Installation of this framework into the geminal-wavefunction-based many-variable Variational Monte Carlo (mVMC) code boosts sampling performance to up to more than $6$ times without changing Markov chain's behavior. The implementation is based on an extension of the BLAS-like instantiation software (BLIS) framework which has optimized kernel for many state-of-the-art processors including Intel Skylake-X, AMD EPYC Rome and Fujitsu A64FX.

preprint2022arXiv

Quantum mechanically driven structural-spin glass in two dimensions at finite temperature

In magnetic materials, spins sometimes freeze into spatially disordered glassy states. Glass forming liquids or structural glasses are found very often in three dimensions. However, in two dimensions(2D) it is believed that both spin glass and structural glass can never exist at a finite temperature because they are destroyed by thermal fluctuations. Using a large-scale quantum Monte Carlo simulation, we discover a quantum-mechanically driven 2D glass phase at finite temperatures. Our platform is an Ising spin model with a quantum transverse field on a frustrated triangular lattice. How the present glass phase is formed is understood by the following three steps. First, by the interplay of geometrical frustration and quantum fluctuation, part of the spins spontaneously form an antiferromagnetic honeycomb spin-superstructure. Then, small randomness in the bond interaction works as a relevant perturbation to this superstructure and breaks it up into {\it domains}, making it a structural glass. The glassiness of the superstructure, in turn, generates an emergent random magnetic field acting on the remaining fluctuating spins and freezes them. The shape of domains thus formed depends sensitively on the quenching process, which is one of the characteristic features of glass, originating from a multi-valley free-energy landscape. The present system consists only of {\it a single} bistable Ising degree of freedom, which naturally does not become a structural glass alone nor a spin glass alone. Nevertheless, a glass having both types of nature emerges in the form of coexisting two-component glasses, algebraic structural-glass and long-range ordered spin-glass. This new concept of glass-forming mechanism opens a way to realize functional glasses even in low dimensional systems.

preprint2022arXiv

Unconventional dual 1D-2D quantum spin liquid revealed by $ab$ $initio$ studies on organic solids family

Organic solids host various electronic phases. Especially, a milestone compound of organic solid, $β'$-$X$[Pd(dmit)$_2$]$_2$ with $X$=EtMe$_3$Sb shows quantum spin-liquid (QSL) properties suggesting a novel state of matter. However, nature of the QSL has been largely unknown. Here, we computationally study five compounds comprehensively with different $X$ using 2D $ab$ $initio$ Hamiltonians and correctly reproduce experimental phase diagram with antiferromagnetic order for $X$=Me$_4$P, Me$_4$As, Me$_4$Sb, Et$_2$Me$_2$As and a QSL for $X$=EtMe$_3$Sb without adjustable parameters. We find that the QSL for $X$=EtMe$_3$Sb exhibits 1D nature characterized by algebraic decay of spin correlation along one direction, while exponential decay in the other direction, indicating dimensional reduction from 2D to 1D. The 1D nature indeed accounts for the experimental specific heat, thermal conductivity and magnetic susceptibility. The identified QSL, however, preserves 2D nature as well consistently with spin fractionalization into spinon with Dirac-like gapless excitations and reveals duality bridging the 1D and 2D QSLs.

preprint2020arXiv

Charge dynamics of correlated electrons: Variational description with inclusion of composite fermions

We propose a method to calculate the charge dynamical structure factors for the ground states of correlated electron systems based on the variational Monte Carlo method. Our benchmarks for the one- and two-dimensional Hubbard models show that inclusion of composite-fermion excitations in the basis set greatly improves the accuracy, in reference to the exact charge dynamical structure factors for clusters. Together with examination for larger systems beyond tractable sizes by the exact diagonalization, our results indicate that the variational Monte Carlo method is a promising way for studies on the nature of charge dynamics in correlated materials such as the copper oxide superconductors, if the composite-fermion excitations are properly included in the restricted Hilbert space of intermediate states in the linear response theory. Our results are consistent with the particle-hole excitations inferred from the single-particle spectral function $A(\mathbf{k}, ω)$ in the literature. We also discuss the importance of incorporating nonlocal composite fermion for more accurate description. Future issues for further improvements are also discussed.

preprint2019arXiv

Ab Initio Study on Superconductivity and Inhomogeneity in Hg-based Cuprate Superconductor

Understanding physics of high-$T_c$ cuprate superconductors remains one of the important problems in materials science. Though a number of diverse theories argue about the superconductivity and competing orders, ab initio and quantitative understanding is lacking. Here, we reproduce the experimental phase diagram of HgBa$_2$CuO$_{4+y}$ by solving its ab initio low-energy effective Hamiltonian without adjustable parameters. It shows a superconducting phase in a wide range of hole density $δ$, and its competition with charge period-4 plus spin period-8 stripe order near $δ\sim 0.1$, in agreement with experimental results including recent X-ray scattering. Then a crucial role of off-site interactions in stabilizing the superconductivity is elucidated with emphasis on charge fluctuations. It also clarifies the condensation energy mainly contributed from the onsite Coulomb interaction. The present achievement will enable deeper, predictable understanding on open issues of the high-$T_c$ superconducting mechanism and promote ab initio studies on strongly correlated electrons beyond parametrized model studies.

preprint2016arXiv

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

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

preprint2016arXiv

Finite-Temperature Variational Monte Carlo Method for Strongly Correlated Electron Systems

A new computational method for finite-temperature properties of strongly correlated electrons is proposed by extending the variational Monte Carlo method originally developed for the ground state. The method is based on the path integral in the imaginary-time formulation, starting from the infinite-temperature state that is well approximated by a small number of certain random initial states. Lower temperatures are progressively reached by the imaginary-time evolution. The algorithm follows the framework of the quantum transfer matrix and finite-temperature Lanczos methods, but we extends them to treat much larger system sizes without the negative sign problem by optimizing the truncated Hilbert space on the basis of the time-dependent variational principle (TDVP). This optimization algorithm is equivalent to the stochastic reconfiguration (SR) method that has been frequently used for the ground state to optimally truncate the Hilbert space. The obtained finite-temperature states allow an interpretation based on the thermal pure quantum (TPQ) state instead of the conventional canonical-ensemble average. Our method is tested for the one- and two-dimensional Hubbard models and its accuracy and efficiency are demonstrated.

preprint2016arXiv

Hidden-Fermion Representation of Self-energy in Pseudogap and Superconducting States of Two-Dimensional Hubbard Model

We study the frequency-dependent structure of electronic self-energy in the pseudogap and superconducting states of the two-dimensional Hubbard model. We present the self-energy calculated with the cellular dynamical mean-field theory systematically in the space of temperature, electron density, and interaction strength. We show that the low-frequency part of the self-energy is well represented by a simple equation, which describes the transitions of an electron to and from a hidden fermionic state. By fitting the numerical data with this simple equation, we determine the parameters characterizing the hidden fermion and discuss its identity. The simple expression of the self-energy offers a way to organize numerical data of this uncomprehended superconducting and pseudogap states, as well as a useful tool to analyze spectroscopic experimental results. The successful description by the simple two-component fermion model supports the idea of "dark" and "bright" fermions emerging from a bare electron as bistable excitations in doped Mott insulators.

preprint2016arXiv

Real-space renormalized dynamical mean field theory

We propose real-space renormalized dynamical mean field theory (rr-DMFT) to deal with large clusters in the framework of a cluster extension of the DMFT. In the rr-DMFT, large clusters are decomposed into multiple smaller clusters through a real-space renormalization. In this work, the renormalization effect is taken into account only at the lowest order with respect to the intercluster coupling, which nonetheless reproduces exactly both the noninteracting and atomic limits. Our method allows us large cluster-size calculations which are intractable with the conventional cluster extensions of the DMFT with impurity solvers, such as the continuous-time quantum Monte Carlo and exact diagonalization methods. We benchmark the rr-DMFT for the two-dimensional Hubbard model on a square lattice at and away from half filling, where the spatial correlations play important roles. Our results on the spin structure factor indicate that the growth of the antiferromagnetic spin correlation is taken into account beyond the decomposed cluster size. We also show that the self-energy obtained from the large-cluster solver is reproduced by our method better than the solution obtained directly for the smaller cluster. When applied to the Mott metal-insulator transition, the rr-DMFT is able to reproduce the reduced critical value for the Coulomb interaction comparable to the large cluster result.

preprint2016arXiv

Self-Optimized Superconductivity Attainable by Interlayer Phase Separation at Cuprate Interfaces

Stabilizing superconductivity at high temperatures and elucidating its mechanism have long been major challenges of materials research in condensed matter physics. Meanwhile, recent progress in nanostructuring offers unprecedented possibilities for designing novel functionalities. Above all, thin films of cuprate and iron-based high-temperature superconductors exhibit remarkably better superconducting characteristics (e.g. higher critical temperatures) than in the bulk, but the underlying mechanism is still not understood. Solving microscopic models suitable for cuprates, here, we demonstrate that at an interface between a Mott insulator and an overdoped non-superconducting metal, the superconducting amplitude stays always pinned at the optimum achieved in the bulk, independently of the carrier concentration in the metal. This is in contrast to the dome-like dependence in bulk superconductors, but consistent with the astonishing independence of the critical temperature from the carrier density x observed at interfaces of La2CuO4 and La2-xSrxCuO4. We furthermore identify a self-organization mechanism as responsible for the pinning at the optimum amplitude: An emergent electronic structure induced by interlayer phase separation eludes bulk phase separation and inhomogeneities that would kill superconductivity in the bulk. Thus, interfaces provide an ideal tool to enhance and stabilize superconductivity. This interfacial example opens up further ways of shaping superconductivity by suppressing competing instabilities, with direct perspectives for designing devices.

preprint2016arXiv

Volume-wise destruction of the antiferromagnetic Mott insulating state through quantum tuning

Metal-to-insulator transitions (MITs) are a dramatic manifestation of strong electron correlations in solids1. The insulating phase can often be suppressed by quantum tuning, i.e. varying a nonthermal parameter such as chemical composi- tion or pressure, resulting in a zero-temperature quantum phase transition (QPT) to a metallic state driven by quantum fluctuations, in contrast to conventional phase transitions driven by thermal fluctuations. Theories of exotic phenomena known to occur near the Mott QPT such as quantum criticality and high-temperature superconductivity often assume a second-order QPT, but direct experimental evidence for either first- or second-order behavior at the magnetic QPT associated with the Mott transition has been scarce and further masked by the superconducting phase in unconventional superconductors. Most measurements of QPTs have been performed by volume-integrated probes, such as neutron scattering, magnetization, and transport, in which discontinuous behavior, phase separation, and spatially inhomogeneous responses are averaged and smeared out, leading at times to misidentification as continuous second-order transitions. Here, we demonstrate through muon spin relaxation/rotation (MuSR) experiments on two archetypal Mott insulating systems, composition-tuned RENiO3 (RE=rare earth element) and pressured-tuned V2O3, that the QPT from antiferromagnetic insulator to paramagnetic metal is first-order: the magnetically ordered volume fraction decreases to zero at the QPT, resulting in a broad region of intrinsic phase separation, while the ordered magnetic moment retains its full value across the phase diagram until it is suddenly destroyed at the QPT. These findings call for further investigation into the role of inelastic soft modes and the nature of dynamic spin and charge fluctuations underlying the transition.

preprint2015arXiv

Ab initio Studies of Magnetism in the Iron Chalcogenides FeTe and FeSe

The iron chalcogenides FeTe and FeSe belong to the family of iron-based superconductors. We study the magnetism in these compounds in the normal state using the ab initio downfolding scheme developed for strongly correlated electron systems. In deriving ab initio low-energy effective models, we employ the constrained GW method to eliminate the double counting of electron correlations originating from the exchange correlations already taken into account in the density functional theory. By solving the derived ab initio effective models, we reveal that the elimination of the double counting is important in reproducing the bicollinear antiferromagnetic order in FeTe, as is observed in experiments. We also show that the elimination of the double counting induces a unique degeneracy of several magnetic orders in FeSe, which may explain the absence of the magnetic ordering. We discuss the relationship between the degeneracy and the recently found puzzling phenomena in FeSe as well as the magnetic ordering found under pressure.

preprint2015arXiv

Exciton Lifetime Paradoxically Enhanced by Dissipation and Decoherence - Toward Efficient Energy Conversion of Solar Cell

Energy dissipation and decoherence are at first glance harmful to acquiring long exciton lifetime desired for efficient photovoltaics. In the presence of both optically forbidden (namely, dark) and allowed (bright) excitons, however, they can be instrumental as suggested in photosynthesis. By simulating quantum dynamics of exciton relaxations, we show that the optimized decoherence that imposes a quantum-to-classical crossover with the dissipation realizes a dramatically longer lifetime. In an example of carbon nanotube, the exciton lifetime increases by nearly two orders of magnitude when the crossover triggers stable high population in the dark exciton.

preprint2015arXiv

Hidden fermionic excitation in the superconductivity of the strongly attractive Hubbard model

We scrutinize the real-frequency structure of the self-energy in the superconducting state of the attractive Hubbard model within the dynamical mean-field theory. Within the strong-coupling superconducting phase which has been understood in terms of the Bose-Einstein condensation in the literature, we find two qualitatively different regions crossing over each other. In one region close to zero temperature, the self-energy depends on the frequency only weakly at low energy. On the other hand, in the region close to the critical temperature, the self-energy shows a pole structure. The latter region becomes more dominant as the interaction becomes stronger. We reveal that the self-energy pole in the latter region is generated by a coupling to a hidden fermionic excitation. The hidden fermion persists in the normal state, where it yields a pseudogap. We compare these properties with those of the repulsive Hubbard model relevant for high-temperature cuprate superconductors, showing that hidden fermions are a key common ingredient in strongly correlated superconductivity.

preprint2015arXiv

Modulated Helical Metals at Magnetic Domain Walls of Pyrochlore Iridium Oxides

Spontaneous symmetry breakings, metal-insulator transitions, and transport properties of magnetic-domain-wall states in pyrochlore iridium oxides are studied by employing a symmetry adapted effective hamiltonian with a slab perpendicular to the (111) direction of the pyrochlore lattice. Emergent metallic domain wall, which has unconventional topological nature with a controllable and mobile metallic layer, is shown to host Fermi surfaces with modulated helical spin textures resembling Rashba metals. The helical nature of the domain-wall Fermi surfaces is experimentally detectable by anomalous Hall conductivity, circular dichroism, and optical Hall conductivity under external magnetic fields. Possible applications of the domain-wall metals to spin-current generation and "half-metallic" conduction are also discussed.

preprint2015arXiv

Nonequilibrium Pump-Probe Photoexcitation as a Tool for Analyzing Unoccupied Equilibrium States of Correlated Electrons

Relaxation of electrons in a Hubbard model coupled to a dissipative bosonic bath is studied to simulate the pump-probe photoemission measurement. From this insight, we propose an experimental method of eliciting unoccupied part of the single-particle spectra at the equilibrium of doped-Mott insulators. We reveal first that effective temperatures of distribution functions and electronic spectra are different during the relaxation, which makes the frequently employed thermalization picture inappropriate. Contrary to the conventional analysis, we show that the unoccupied spectra at equilibrium can be detected as the states that relax faster.

preprint2015arXiv

Quantum Spin Liquid in Spin 1/2 J1-J2 Heisenberg Model on Square Lattice: Many-Variable Variational Monte Carlo Study Combined with Quantum-Number Projections

The nature of quantum spin liquids is studied for the spin-$1/2$ antiferromagnetic Heisenberg model on a square lattice containing exchange interactions between nearest-neighbor sites, $J_1$, and those between next-nearest-neighbor sites, $J_2$. We perform variational Monte Carlo simulations together with the quantum-number-projection technique and clarify the phase diagram in the ground state together with its excitation spectra. We obtain the nonmagnetic phase in the region $0.4< J_2/J_1\le 0.6$ sandwiched by the staggered and stripe antiferromagnetic (AF) phases. Our direct calculations of the spin gap support the notion that the triplet excitation from the singlet ground state is gapless in the range of $0.4 < J_2/J_1 \le 0.5$, while the gapped valence-bond-crystal (VBC) phase is stabilized for $0.5 < J_2/J_1 \le 0.6$. The VBC order is likely to have the columnar symmetry with a spontaneous symmetry breaking of the $C_{4v}$ symmetry. The power-law behaviors of the spin-spin and dimer-dimer correlation functions in the gapless region are consistent with the emergence of the algebraic quantum-spin-liquid phase (critical phase). The exponent of the spin correlation $\langle S(0)S(r)\rangle \propto 1/r^{z+η}$ at a long distance $r$ appears to increase from $z+η\sim 1$ at $J_2/J_1\sim0.4$ toward the continuous transition to the VBC phase at $J_1/J_1\sim0.5$. Our results, however, do not fully exclude the possibility of a direct quantum transition between the staggered AF and VBC phases with a wide critical region and deconfined criticality.

preprint2015arXiv

Stabilization of Topological Insulator Emerging from Electron Correlations on Honeycomb Lattice and Its Possible Relevance in Twisted Bilayer Graphene

Realization and design of topological insulators emerging from electron correlations, called topological Mott insulators (TMIs), is pursued by using mean-field approximations as well as multi-variable variational Monte Carlo (MVMC) methods for Dirac electrons on honeycomb lattices. The topological insulator phases predicted in the previous studies by the mean-field approximation for an extended Hubbard model on the honeycomb lattice turn out to disappear, when we consider the possibility of a long-period charge-density-wave (CDW) order taking over the TMI phase. Nevertheless, we further show that the TMI phase is still stabilized when we are able to tune the Fermi velocity of the Dirac point of the electron band. Beyond the limitation of the mean-field calculation, we apply the newly developed MVMC to make accurate predictions after including the many-body and quantum fluctuations. By taking the extrapolation to the thermodynamic and weak external field limit, we present realistic criteria for the emergence of the topological insulator caused by the electron correlations. By suppressing the Fermi velocity to a tenth of that of the original honeycomb lattice, the topological insulator emerges in an extended region as a spontaneous symmetry breaking surviving competitions with other orders. We discuss experimental ways to realize it in a bilayer graphenesystem.

preprint2015arXiv

Tensor network algorithm by coarse-graining tensor renormalization on finite periodic lattices

We develop coarse-graining tensor renormalization group algorithms to compute physical properties of two-dimensional lattice models on finite periodic lattices. Two different coarse-graining strategies, one based on the tensor renormalization group and the other based on the higher-order tensor renormalization group, are introduced. In order to optimize the tensor-network model globally, a sweeping scheme is proposed to account for the renormalization effect from the environment tensors under the framework of second renormalization group. We demonstrate the algorithms by the classical Ising model on the square lattice and the Kitaev model on the honeycomb lattice, and show that the finite-size algorithms achieve substantially more accurate results than the corresponding infinite-size ones.

preprint2015arXiv

Time-dependent many-variable variational Monte Carlo method for nonequilibrium strongly correlated electron systems

We develop a time-dependent variational Monte Carlo (t-VMC) method for quantum dynamics of strongly correlated electrons. The t-VMC method has been recently applied to bosonic systems and quantum spin systems. Here, we propose a time-dependent trial wave function with many variational parameters, which is suitable for nonequilibrium strongly correlated electron systems. As the trial state, we adopt the generalized pair-product wave function with correlation factors and quantum-number projections. This trial wave function has been proven to accurately describe ground states of strongly correlated electron systems. To show the accuracy and efficiency of our trial wave function in nonequilibrium states as well, we present our benchmark results for relaxation dynamics during and after interaction quench protocols of fermionic Hubbard models. We find that our trial wave function well reproduces the exact results for the time evolution of physical quantities such as energy, momentum distribution, spin structure factor, and superconducting correlations. These results show that the t-VMC with our trial wave function offers an efficient and accurate way to study challenging problems of nonequilibrium dynamics in strongly correlated electron systems.

preprint2015arXiv

Variational Monte Carlo Method in the Presence of Spin-Orbit Interaction and Its Application to the Kitaev and Kitaev-Heisenberg Model

We propose an accurate variational Monte Carlo method applicable in the presence of the strong spin-orbit interaction. Our variational wave functions consist of generalized Pfaffian-Slater wave functions that involve mixtures of singlet and triplet Cooper pairs, Jastrow-Gutzwiller-type pro- jections, and quantum number projections. The generalized wave functions allow any symmetry breaking states, ranging from magnetic and/or charge ordered states to superconducting states and their fluctuations, on equal footing without any ad hoc ansatz for variational wave functions. We de- tail our optimization scheme for the generalized Pfaffian-Slater wave functions with complex-number variational parameters. Generalized quantum number projections are also introduced, which im- poses the conservation of not only spin and momentum quantum numbers but also Wilson loops. As a demonstration of the capability of the present variational Monte Carlo method, the accuracy and efficiency of the present method is tested for the Kitaev and Kitaev-Heisenberg models. The Kitaev model serves as a critical benchmark of the present method: The exact ground state of the model is a typical gapless quantum spin liquid far beyond the applicability range of simple mean-field wave functions. The newly introduced quantum number projections precisely reproduce the ground state degeneracy of the Kitaev spin liquids, in addition to their ground state energy. Our framework offers accurate solutions for the systems where strong electron correlation and spin-orbit interaction coexist.

preprint2014arXiv

Electron Correlation Effects on Topological Phases

Topological insulators are found in materials that have elements with strong spin orbit interaction. However, electron Coulomb repulsion also potentially generates the topological insulators as well as Chern insulators by the mechanism of spontaneous symmetry breaking, which is called topological Mott insulators. The quantum criticality of the transition to the topological Mott insulators from zero-gap semiconductors follows unconventional universality distinct from the Landau-Ginzburg-Wilson scenario. On the pyrochlore lattice, the interplay of the electron correlation and the spin orbit interaction provides us in a rich phase diagram not only with simple topological insulators but also with Weyl semimetal and topologically distinct antiferromagnetic phases. Magnetic domain wall of the all-in-all-out type antiferromagnetic order offers a promising candidate of magnetically controlled transport, because, even when the Weyl points disappears, the domain wall maintains robust gapless excitations with Fermi surfaces around it embedded in the bulk insulator and bears uniform magnetization simultaneously. The ingap state is protected by a mechanism similar to the solitons in polyacetylene. Puzzling experimental results of pyrochlore iridates are favorably compared with the prediction of the domain wall theory.

preprint2014arXiv

Evidence for hidden fermion that triggers high-temperature superconductivity in cuprates

In superconductors, electrons bound into Cooper pairs conduct a dissipationless current. The strength of the Cooper pairs scales with the value of the critical transition temperature (Tc). In cuprate high-Tc superconductors, however, the pairing mechanism is still unexplained. Here we unveil why in the cuprates the Cooper pairs are so strongly bound to work out the extraordinary high Tc. From one-to-one correspondence between numerical simulation on a microscopic cuprate model and a simple two-component fermion model, we show that hidden fermions emerge from the strong electron correlation and give birth to the strongly bound Cooper pairs. This mechanism is distinct from a conventional pairing mediated by some bosonic glue, such as phonons in conventional superconductors. The hidden fermions survive even above Tc and generate the strange-metal pseudogap phase. This reveals an unprecedented direct relationship between the pseudogap phase and superconductivity in the cuprates.

preprint2014arXiv

First-Principles Study of the Honeycomb-Lattice Iridates Na2IrO3 in the Presence of Strong Spin-Orbit Interaction and Electron Correlations

An effective low-energy Hamiltonian of itinerant electrons for iridium oxide Na2IrO3 is derived by an ab initio downfolding scheme. The model is then reduced to an effective spin model on a honeycomb lattice by the strong coupling expansion. Here we show that the ab initio model contains spin-spin anisotropic exchange terms in addition to the extensively studied Kitaev and Heisenberg exchange interactions, and allows to describe the experimentally observed zigzag magnetic order, interpreted as the state stabilized by the antiferromagnetic coupling of the ferromagnetic chains. We clarify possible routes to realize quantum spin liquids from existing Na2IrO3.

preprint2014arXiv

Metallic Interface Emerging at Magnetic Domain Wall of Antiferromagnetic Insulator---Fate of Extinct Weyl Electrons

Topological insulators, in contrast to ordinary semiconductors, accompany protected metallic surfaces described by Dirac-type fermions. Here, we theoretically show another emergent two-dimensional metal embedded in the bulk insulator is realized at a magnetic domain wall. The domain wall has long been studied as ingredients of both old-fashioned and leading-edge spintronics. The domain wall here, as an interface of seemingly trivial antiferromagnetic insulators, emergently realizes a functional interface preserved by zero modes with robust two-dimensional Fermi surfaces, where pyrochlore iridium oxides proposed to host condensed-matter realization of Weyl fermions offer such examples at low temperatures. The existence of ingap states pinned at domain walls, theoretically resembling spin/charge solitons in polyacetylene, and protected as the edge of hidden one-dimensional weak Chern insulators characterized by a zero-dimensional class A topological invariant, solves experimental puzzles observed in R2Ir2O7 with rare earth elements R. The domain wall realizes a novel quantum confinement of electrons and embosses a net uniform magnetization, which enables magnetic control of electronic interface transports beyond semiconductor paradigm.

preprint2014arXiv

Origin of high-Tc superconductivity in doped Hubbard models and their extensions: Roles of uniform charge fluctuations

Doped Hubbard model is a simple model for the high-Tc cuprate superconductors, while its ground state remains a challenge. Here, by performing state-of-the-art variational Monte Carlo calculations for the strong-coupling Hubbard model, we find evidences that the d-wave superconducting phase emerges always near the phase separation region and the superconducting order has one-to-one correspondence with the enhancement of charge compressibility. The order as well as the phase separation are vulnerable to realistic intersite Coulomb interaction while the superexchange interaction enhances both. An appropriate combination of these two widens the stable superconducting phase.

preprint2014arXiv

Superconductivity and its mechanism in an ab initio model for electron-doped LaFeAsO

Two families of high temperature superconductors whose critical temperatures are higher than 50K are known. One is the copper oxides and the other is the iron-based superconductors. Comparisons of mechanisms between these two in terms of common ground as well as distinctions will greatly help in searching for higher Tc superconductors. However, studies on mechanisms for the iron family based on first principles calculations are few. Here we first show that the superconductivity emerges in the state-of-the-art numerical calculations for an ab initio multi-orbital model of an electron-doped iron-based superconductor LaFeAsO, in accordance with experimental observations. Then the mechanism of the superconductivity is identified as enhanced uniform density fluctuations by one-to-one correspondence with the instability toward inhomogeneity driven by first-order antiferromagnetic and nematic transitions. Despite many differences, certain common features with the copper oxides are also figured out in terms of the underlying orbital selective Mottness found in the iron family.

preprint2014arXiv

Universal departure from Johnson-Nyquist relation caused by limited resolution

Exploiting the two-point measurement statistics, we propose a quantum measurement scheme of current with limited resolution of electron counting. Our scheme is equivalent to the full counting statistics in the long-time measurement with the ideal resolution, but is theoretically extended to take into account the resolution limit of actual measurement devices. Applying our scheme to a resonant level model, we show that the limited resolution of current measurement gives rise to a positive excess noise, which leads to a deviation from the Johnson-Nyquist relation. The deviation exhibits universal single-parameter scaling with the scaling variable $Q\equiv S_{\rm{}M}/S_0$, which represents the degree of the insufficiency of the resolution. Here, $S_0$ is the intrinsic noise, and $S_{\rm{}M}$ is the positive quantity that has the same dimension as $S_0$ and is defined solely by the measurement scheme. For the lack of the ideal resolution, the deviation emerges for $Q<1$ as $2\exp[-(2π)^2/Q]$ having an essential singularity at $Q=0$, which followed by the square root dependence $\sqrt{Q/4π}$ for $Q\gg1$. Our findings offer an explanation for the anomalous enhancement of noise temperature observed in Johnson noise thermometry.

preprint2014arXiv

Variational Monte Carlo Method for Electron-Phonon Coupled Systems

We develop a variational Monte Carlo (VMC) method for electron-phonon coupled systems. The VMC method has been extensively used for investigating strongly correlated electrons over the last decades. However, its applications to electron-phonon coupled systems have been severely restricted because of its large Hilbert space. Here, we propose a variational wave function with a large number of variational parameters which is suitable and tractable for systems with electron-phonon coupling. In the proposed wave function, we implement an unexplored electron-phonon correlation factor which takes into account the effect of the entanglement between electrons and phonons. The method is applied to systems with diagonal electron-phonon interactions, i.e. interactions between charge densities and lattice displacements (phonons). As benchmarks, we compare VMC results with previous results obtained by the exact diagonalization, the Green function Monte Carlo and the density matrix renormalization group for the Holstein and Holstein-Hubbard model. From these benchmarks, we show that the present method offers an efficient way to treat strongly coupled electron-phonon systems.

preprint2013arXiv

Improved Multi-Variable Variational Monte Carlo Method Examined by High-Precision Calculations of One-Dimensional Hubbard Model

We revisit the accuracy of the variational Monte Carlo (VMC) method by taking an example of ground state properties for the one-dimensional Hubbard model. We start from the variational wave functions with the Gutzwiller and long-range Jastrow factor introduced by Capello et al. [Phys. Rev. B 72, 085121 (2005)] and further improve it by considering several quantum-number projections and a generalized one-body wave function. We find that the quantum spin projection and total momentum projection greatly improve the accuracy of the ground state energy within 0.5% error, for both small and large systems at half filling. Besides, the momentum distribution function n(k) at quarter filling calculated up to 196 sites allows us direct estimate of the critical exponents of the charge correlations from the power-law behavior of n(k) near the Fermi wave vector. Estimated critical exponents well reproduce those predicted by the Tomonaga-Luttinger theory.

preprint2013arXiv

Phase Diagram Structure of Topological Mott Transition for Zero-gap Semiconductors beyond Conventional Landau-Ginzburg-Wilson Scenario

We show that a wide class of unconventional quantum criticality emerges when orbital currents cause quantum phase transitions from zero-gap semiconductors such as Dirac fermions to topological insulator (TI) or Chern insulator (CI). Changes in Fermi surface topology concomitant with (SU(2) or time reversal) symmetry breakings generate quantum critical lines (QCL) even beyond the quantum critical point. This QCL running at temperature $T=0$ separates two distinct topological phases. This is in contrast to the simple termination of the finite temperature critical line at the quantum critical point without any extension of it at $T=0$. Topology change causes the unconventionality beyond the concept of simple spontaneous symmetry breaking assumed in the conventional Landau-Ginzburg-Wilson (LGW) scenario. The unconventional universality implied by mean-field critical exponents beta>1/2 and delta<3 is protected by the existence of the quantum critical line. It emerges for several specific lattice models including the honeycomb, kagome, diamond and pyrochlore lattices. We also clarify phase diagrams of the topological phases in these lattices at finite temperatures.

preprint2013arXiv

Quantum Monte Carlo Simulations for Stacked Spin Ladder Systems Containing Low Concentrations of Non-Magnetic Impurities; Application to the Low Temperature Broadening of NMR-Spectra in SrCu2O3

We present a Quantum Monte Carlo study for Heisenberg spin-1/2 two leg ladder systems doped with non magnetic impurities. The simulations are applied to the doped spin ladder compound SrCu2O3 doped with Zn, where a large broadening of the 65Cu NMR lines has been observed in experiment at low temperatures but far above the Neel temperature. We find that interladder-couplings with a sizeable coupling in stacking direction are required to describe the line broadening, which can not be explained by considering a single ladder only. Around a single impurity, spin correlations cause an exponentially decaying antiferromagnetic local magnetization in a magnetic field. We develop an effective model for the local magnetization of systems with many randomly distributed impurities, with few parameters which can be extracted out of QMC calculations with a single impurity. The broadening arises from a drag effect, where the magnetization around an impurity works as an effective field for spins on the neighboring ladders, causing a non-exponentially decaying magnetization cloud around the impurity. Our results show that even for impurity concentrations as small as x = 0.001 and x = 0.0025, the broadening effect is large, in good quantitative agreement with experiment. We also develop a simple model for the effective interaction of two impurity spins.

preprint2012arXiv

Ab initio two-dimensional multiband low-energy models of EtMe_3Sb[Pd(dmit)_2]_2 and κ-(BEDT-TTF)_2Cu(NCS)_2 with comparisons to single-band models

We present ab initio two-dimensional extended Hubbard-type multiband models for EtMe_3Sb[Pd(dmit)_2]_2 and κ-(BEDT-TTF)_2Cu(NCS)_2, after a downfolding scheme based on the constrained random phase approximation (cRPA) and maximally-localized Wannier orbitals, together with the dimensional downfolding. In the Pd(dmit)_2 salt, the antibonding state of the highest occupied molecular orbital (HOMO) and the bonding/antibonding states of the lowest unoccupied molecular orbital (LUMO) are considered as the orbital degrees of freedom, while, in the κ-BEDT-TTF salt, the HOMO-antibonding/bonding states are considered. Accordingly, a three-band model for the Pd(dmit)_2 salt and a two-band model for the κ-(BEDT-TTF) salt are derived. We derive single band models for the HOMO-antibonding state for both of the compounds as well.

preprint2012arXiv

Ab-initio Low-Energy Model of Transition-Metal-Oxide Heterostructure LaAlO3/SrTiO3

We develop the multi-scale ab-initio scheme for correlated electrons (MACE) for transitionmetal-oxide heterostructures, and determine the parameters of the low-energy effective model. By separating Ti t2g bands near the Fermi level from the global Kohn-Sham (KS) bands of LaAlO3/SrTiO3 which are highly entangled with each other, we are able to calculate the parameters of the low-energy effective model of the interface with the help of constrained random phase approximation (cRPA). The on-site energies of the Ti t2g orbitals in the 1stlayer is about 650 meV lower than those in the 2nd-layer. In the 1st-layer, the transfer integral of the Ti t2g orbital is nearly the same as that of the bulk SrTiO3, while the effective screened Coulomb interaction becomes about 10 percent larger than that of the bulk SrTiO3. The differences of the parameters from the bulk SrTiO3 reduce rapidly with increasing distance from the interface. Our present versatile method makes it possible to derive effective ab-initio low-energy models and allows studying interfaces of strongly correlated electron systems from first principles.

preprint2012arXiv

Derivation of Static Low-Energy Effective Models by ab initio Downfolding Method without Double Counting of Coulomb Correlations: Application to SrVO3, FeSe and FeTe

Derivation of low-energy effective models by a partial trace summation of the electronic degrees of freedom far away from the Fermi level, called downfolding, is reexamined. We propose an improved formalism free from the double-counting of electron correlation in the low-energy degrees of freedom. In this approach, the exchange-correlation energy in the local density approximation (LDA) is replaced with the constrained self-energy corrections defined by the sum of the contribution from eliminated high-energy degrees of freedom, Sigma H and that from the frequency-dependent part of the partially screened interaction, Delta Sigma L. We apply the formalism to SrVO3 as well as to two iron-based superconductors, FeSe and FeTe. The resultant bandwidths of the effective models are nearly the same as those of the previous downfolding formalism because of striking cancellations of Sigma H and Delta Sigma L. In SrVO3, the resultant bandwidth of the effective low-energy model is 2.56 eV (in comparison to LDA: 2.58 eV, full GW approximation: 2.19 eV). However, in the non-degenerate multi-band materials such as FeSe and FeTe, the momentum dependent self-energy effects yield substantial modifications of the band structures and relative shifts of orbital-energy levels of the effective models. The FeSe model indicates a substantial downward shift of the x2-y2 orbital level, which may lead to an increase in the filling of the x2-y2 orbital above half filling. It suggests a destruction of the orbital selective Mott phase of the x2-y2 orbital in agreement with the experimental absence of the antiferromagnetic phase.

preprint2012arXiv

Effective Onsite Interaction for Dynamical Mean-Field Theory

A scheme to incorporate non-local polarizations into the dynamical mean-field theory (DMFT) and a tailor-made way to determine the effective interaction for the DMFT are systematically investigated. Applying it to the two-dimensional Hubbard model, we find that non-local polarizations induce a non-trivial filling-dependent anti-screening effect for the effective interaction. The present scheme combined with density functional theory offers an ab initio way to derive effective onsite interactions for the impurity problem in DMFT. We apply it to SrVO3 and find that the anti-screening competes with the screening caused by the off-site interaction.

preprint2012arXiv

High-temperature Superconductivity in Layered Nitrides β-Li$_x$MNCl (M = Ti, Zr, Hf): Insights from Density-functional Theory for Superconductors

We present an ab initio analysis with density functional theory for superconductors (SCDFT) to understand the superconducting mechanism of doped layered nitrides β-Li$_x$MNCl (M=Ti, Zr, and Hf). The current version of SCDFT is based on the Migdal-Eliashberg theory and has been shown to reproduce accurately experimental superconducting-transition temperatures Tc of a wide range of phonon-mediated superconductors. In the present case, however, our calculated Tc$\leq$4.3 K (M=Zr) and $\leq$10.5 K (M=Hf) are found to be less than a half of the experimental Tc. In addition, Tc obtained in the present calculation increases with the doping concentration x, opposite to that observed in the experiment. Our results indicate that we need to consider some elements missing in the present SCDFT based on the Migdal-Eliashberg theory.

preprint2012arXiv

Impurity Effects on Superconductivity on Surfaces of Topological Insulators

A two-dimensional superconductor (SC) on surfaces of topological insulators (TIs) is a mixture of s-wave and helical p-wave components when induced by s-wave interactions, since spin and momentum are correlated. On the basis of the Abrikosov-Gor'kov theory, we reveal that unconventional SCs on the surfaces of TIs are stable against time-reversal symmetric (TRS) impurities within a region of small impurity concentration. Moreover, we analyze the stability of the SC on the surfaces of TIs against impurities beyond the perturbation theory by solving the real-space Bogoliubov-de Gennes equation for an effective tight-binding model of a TI. We find that the SC is stable against strong TRS impurities. The behaviors of bound states around an impurity suggest that the SC on the surfaces of TIs is not a topological SC.

preprint2012arXiv

Mott Transition and Phase Diagram of $κ$-(BEDT-TTF)2Cu(NCS)2 Studied by Two-Dimensional Model Derived from Ab initio Method

We present an ab initio analysis for the ground-state properties of a correlated organic compound $κ$-(BEDT-TTF)2Cu(NCS)2. First, we derive an effective two-dimensional low-energy model from first principles, having short-ranged transfers and short-ranged Coulomb and exchange interactions. Then, we perform many-variable variational Monte Carlo calculations for this model and draw a ground-state phase diagram as functions of scaling parameters for the onsite and off-site interactions. The phase diagram consists of three phases; a paramagnetic metallic phase, an antiferromagnetic (Mott) insulating phase, and a charge-ordered insulating phase. In the phase diagram, the parameters for the real compound are close to the first-order Mott transition, being consistent with experiments. We show that the off-site Coulomb and exchange interactions affect the phase boundary; (i) they appreciably stabilize the metallic state against the Mott insulating phase and (ii) enhance charge fluctuations in a wide parameter region in the metallic phase. We observe arc-like structure in Fermi surface around the region where the charge fluctuations are enhanced. Possible relevance of the charge fluctuations to the experimentally observed dielectric anomaly in the $κ$-BEDT-TTF family compounds is also pointed out.

preprint2011arXiv

Ab initio Evidence for Strong Correlation Associated with Mott Proximity in Iron-based Superconductors

We predict that iron-based superconductors discovered near d6 configuration (5 Fe 3d orbitals filled by 6 electrons) is located on the foot of an unexpectedly large dome of correlated electron matter centered at the Mott insulator at d5 (namely, half filling). This is based on the many-variable variational Monte-Carlo results for ab initio low-energy models derived by the downfolding. The d5 Mott proximity extends to subsequent emergence of incoherent metals, orbital differentiations due to the Mott physics and Hund's-rule coupling, followed by antiferromagnetic quantum criticality, in quantitative accordance with available experiments.

preprint2011arXiv

Cluster-size dependence in cellular dynamical mean-field theory

We examine the cluster-size dependence of the cellular dynamical mean-field theory (CDMFT) applied to the two-dimensional Hubbard model. Employing the continuous-time quantum Monte Carlo method as the solver for the effective cluster model, we obtain CDMFT solutions for 4-, 8-, 12-, and 16-site clusters at a low temperature. Comparing various periodization schemes, which are used to construct the infinite-lattice quantities from the cluster results, we find that the cumulant periodization yields the fastest convergence for the hole-doped Mott insulator where the most severe size dependence is expected. We also find that the convergence is much faster around (0,0) and (pi/2,pi/2) than around (pi,0) and (pi,pi). The cumulant-periodized self-energy seems to be close to its thermodynamic limit already for a 16-site cluster in the range of parameters studied. The 4-site results remarkably agree well with the 16-site results, indicating that the previous studies based on the 4-site cluster capture the essence of the physics of doped Mott insulators.

preprint2011arXiv

Cofermion Theory for Pseudogap Phenomena and Superconducting Mechanism of Underdoped Cuprate Superconductors

We study pseudogap phenomena and Fermi-arc formation experimentally observed in typical two dimensional doped Mott insulators, namely, underdoped cuprate superconductors. To develop a physically unequivocal theory, we start from the slave-boson mean-field theory for the Hubbard model on a square lattice. Our crucial step is to further take into account the charge dynamics and fluctuations. The extra charge fluctuations seriously modify low-energy single-particle spectra of doped Mott insulators near the Fermi level: An electron added around an empty site (or a hole added around a doubly occupied site) constitutes composite fermion (cofermion), called holo-electron (or doublo-hole) at low energy in distinction from the normal quasiparticles. These unexplored composite fermions substantiate the extra charge fluctuation. We show that the quasiparticles hybridize with the holo-electrons and doublo-holes. The resultant hybridization gap is identified as the pseudogap observed in the underdoped region of the high-T_c cuprates. Because the Fermi level crosses the top (bottom) of the low-energy band formed just below (above) the hybridization gap in the hole-doped (electron-doped) case, it causes a Fermi-surface reconstruction, namely, a topological change in the Fermi surface forced by the penetration of zeros of the quasiparticle Green function. This reconstruction signals the emergence of a non-Fermi-liquid phase. The pseudogap, and the resultant formation of pocket or arc of the Fermi surface reproduce the experimental results for the cuprate superconductors in the underdoped region. The pairing channel opens not only between two quasiparticles but also between a quasiparticle and a cofermion. This pairing solves the puzzle of the dichotomy between the d-wave superconductivity and the precursors of the the insulating gap in the antinodal region.

preprint2011arXiv

Stability of Unconventional Superconductivity on Surfaces of Topological Insulators

Superconductivity on the surface of topological insulators is known to be anisotropic and unconventional in that the symmetry is the mixture of s-wave and nodeless p-wave component. In contrast to Anderson's theorem for the insensitivity of the s-wave superconducting critical temperature to the nonmagnetic (time-reversal symmetric (TRS)) impurities, anisotropic superconductors including nodeless p-wave one are in general fragile even with small concentration of the TRS impurities. By employing the Abrikosov-Gor'kov theory, we clarify that this type of unconventional superconductivity emergent on the surface state of the strong topological insulators robustly survive against TRS impurities.

preprint2011arXiv

Theory of pseudogap and superconductivity in doped Mott insulators

Underdoped Mott insulators provide us with a challenge of many-body physics. Recent renewed understanding is discussed in terms of the evolution of pole and zero structure of the single-particle Green's function. Pseudogap as well as Fermi arc/pocket structure in the underdoped cuprates is well reproduced from the recent cluster extension of the dynamical mean-field theory. Emergent coexisting zeros and poles set the underdoped Mott insulator apart from the Fermi liquid, separated by topological transitions. Cofermion proposed as a generalization of exciton in the slave-boson framework accounts for the origin of the zero surface formation. The cofermion-quasiparticle hybridization gap offers a natural understanding of the pseudogap and various unusual Mottness. Furthermore the cofermion offers a novel pairing mechanism, where the cofermion has two roles: It reinforces the Cooper pair as a pair partner of the quasiparticle and acts as a glue as well. It provides a strong insight for solving the puzzle found in the dichotomy of the gap structure.

preprint2011arXiv

Topological Insulators from Spontaneous Symmetry Breaking Induced by Electron Correlation on Pyrochlore Lattices

We study an extended Hubbard model with the nearest-neighbor Coulomb interaction on the pyrochlore lattice at half filling. An interaction-driven insulating phase with nontrivial Z_2 invariants emerges at the Hartree-Fock mean-field level in the phase diagram. This topological insulator phase competes with other ordered states and survives in a parameter region surrounded by a semimetal, antiferromagnetic and charge ordered insulators. The symmetries of these phases are group-theoretically analyzed. We also show that the ferromagnetic interaction enhances the stability of the topological phase.

preprint2010arXiv

Ab initio Low-Dimensional Physics Opened Up by Dimensional Downfolding: Application to LaFeAsO

An ab initio downfolding method is formulated to construct low-dimensional models for correlated electrons. In addition to the band downfolding by constrained random phase approximation formulated for 3D models, screening away from the target layer (chain) is further involved. Eliminating the off-target degrees of freedom, namely, dimensional downfolding yields ab initio low-dimensional models. The method is applied to derive a 2D model for a layered superconductor LaFeAsO, where the interlayer screening crucially makes the effective interaction short ranged and reduces the onsite Coulomb interactions by 10-20 \% from the 3D model for the 5 iron-3d orbitals.

preprint2010arXiv

Comparison of Ab initio Low-Energy Models for LaFePO, LaFeAsO, BaFe2As2, LiFeAs, FeSe and FeTe: Electron Correlation and Covalency

Effective low-energy Hamiltonians for several different families of iron-based superconductors are compared after deriving them from the downfolding scheme based on first-principles calculations. Systematic dependences of the derived model parameters on the families are elucidated, many of which are understood from the systematic variation of the covalency between Fe-3d and pnictogen-/chalcogen-p orbitals.

preprint2010arXiv

Composite-Fermion Theory for Pseudogap, Fermi Arc, Hole Pocket, and Non-Fermi-Liquid of Underdoped Cuprate Superconductors

We propose that an extension of the exciton concept to doped Mott insulators offers a fruitful insight into challenging issues of the copper oxide superconductors. In our extension, new fermionic excitations called cofermions emerge in conjunction to generalized excitons. The cofermions hybridize with conventional quasiparticles. Then a hybridization gap opens, and is identified as the pseudogap observed in the underdoped cuprates. The resultant Fermi-surface reconstruction naturally explains a number of unusual properties of the underdoped cuprates, such as the Fermi arc and/or pocket formation.

preprint2010arXiv

Doped high-Tc cuprate superconductors elucidated in the light of zeros and poles of electronic Green's function

We study electronic structure of hole- and electron-doped Mott insulators in the two-dimensional Hubbard model to reach a unified picture for the normal state of cuprate high-Tc superconductors. By using a cluster extension of the dynamical mean-field theory, we demonstrate that structure of coexisting zeros and poles of the single-particle Green's function holds the key to understand Mott physics in the underdoped region. We show evidence for the emergence of non-Fermi-liquid phase caused by the topological quantum phase transition of Fermi surface by analyzing low-energy charge dynamics. The spectra calculated in a wide range of energy and momentum reproduce various anomalous properties observed in experiments for the high-Tc cuprates. Our results reveal that the pseudogap in hole-doped cuprates has a d-wave-like structure only below the Fermi level, while it retains non-d-wave structure with a fully opened gap above the Fermi energy even in the nodal direction due to a zero surface extending over the entire Brillouin zone. In addition to the non-d-wave pseudogap, the present comprehensive identifications of the spectral asymmetry as to the Fermi energy, the Fermi arc, and the back-bending behavior of the dispersion, waterfall, and low-energy kink, in agreement with the experimental anomalies of the cuprates, do not support that these originate from (the precursors of) symmetry breakings such as the preformed pairing and the d-density wave fluctuations, but support that they are direct consequences of the proximity to the Mott insulator. Several possible experiments are further proposed to prove or disprove our zero mechanism.

preprint2010arXiv

Electronic and magnetic properties of metallic phases under coexisting short-range interaction and diagonal disorder

We study a three-dimensional Anderson-Hubbard model under the coexistence of short-range interaction and diagonal disorder within the Hartree-Fock approximation. We show that the density of states at the Fermi energy is suppressed in the metallic phases near the metal-insulator transition as a proximity effect of the soft Hubbard gap in the insulating phases. The transition to the insulator is characterized by a vanishing DOS in contrast to formation of a quasiparticle peak at the Fermi energy obtained by the dynamical mean field theory in pure systems. Furthermore, we show that there exist frozen spin moments in the paramagnetic metal.

preprint2010arXiv

Electronic Structure Calculation by First Principles for Strongly Correlated Electron Systems

Recent trends of ab initio studies and progress in methodologies for electronic structure calculations of strongly correlated electron systems are discussed. The interest for developing efficient methods is motivated by recent discoveries and characterizations of strongly correlated electron materials and by requirements for understanding mechanisms of intriguing phenomena beyond a single-particle picture. A three-stage scheme is developed as renormalized multi-scale solvers (RMS) utilizing the hierarchical electronic structure in the energy space. It provides us with an ab initio downfolding of the global band structure into low-energy effective models followed by low-energy solvers for the models. The RMS method is illustrated with examples of several materials. In particular, we overview cases such as dynamics of semiconductors, transition metals and its compounds including iron-based superconductors and perovskite oxides, as well as organic conductors of kappa-ET type.

preprint2010arXiv

Magnetic Properties of Ab initio Model for Iron-Based Superconductors LaFeAsO

By using variational Monte Carlo method, we examine an effective low-energy model for LaFeAsO derived from an ab initio downfolding scheme. We show that quantum and many-body fluctuations near a quantum critical point largely reduce the antiferromagnetic (AF) ordered moment and the model not only quantitatively reproduces the small ordered moment in LaFeAsO, but also explains the diverse dependence on LaFePO, BaFe2As2 and FeTe. We also find that LaFeAsO is under large orbital fluctuations, sandwiched by the AF Mott insulator and weakly correlated metals. The orbital fluctuations and Dirac-cone dispersion hold keys for the diverse magnetic properties.

preprint2010arXiv

Mott Physics on Helical Edges of 2D Topological Insulators

We study roles of electron correlations on topological insulators on the honeycomb lattice with the spin-orbit interaction. Accurate variational Monte Carlo calculations show that the increasing on-site Coulomb interactions cause a strong suppression of the charge Drude weight in the helical-edge metallic states leading to a presumable Mott transition from a conventional topological insulator to an edge Mott insulator before a transition to a bulk antiferromagnetic insulator. The intermediate bulk-topological and edge-Mott-insulator phase has a helical spin-liquid character with the protected time-reversal symmetry.

preprint2010arXiv

Theoretical evidence for strong correlations and incoherent metallic state in FeSe

The role of electronic Coulomb correlations in iron-based superconductors is an important open question. We provide theoretical evidence for strong correlation effects in the FeSe compound, based on dynamical mean field calculations. A lower Hubbard band is found in the spectral properties. Moreover, together with significant orbital-dependent mass enhancements, we find that the normal state is a bad metal over an extended temperature range, implying a non-Fermi liquid. Predictions for angle-resolved photoemission spectroscopy are made.

preprint2010arXiv

Theory of Electron Transport near Anderson-Mott Transitions

We present a theory of the DC electron transport in insulators near Anderson-Mott transitions under the influence of coexisting electron correlation and randomness. At sufficiently low temperatures, the DC electron transport in Anderson-Mott insulators is determined by the single-particle density of states (DOS) near the Fermi energy. Anderson insulators, caused by randomness, are characterized by a nonzero DOS at the Fermi energy. However, recently, the authors proposed that coexisting randomness and short-ranged interaction in insulators open a soft Hubbard gap in the DOS, and the DOS vanishes only at the Fermi energy. Based on the picture of the soft Hubbard gap, we derive a formula for the critical behavior for the temperature dependence of the DC resistivity. Comparisons of the present theory with experimental results of electrostatic carrier doping into an organic conductor kappa-(BEDT-TTF)_2Cu[N(CN)_2]Br demonstrate the evidence for the present soft-Hubbard scaling.

preprint2010arXiv

Unconventional quantum criticality emerging as a new common language of transition-metal compounds, heavy-fermion systems, and organic conductors

We analyze and overview several different unconventional quantum criticalities. One origin of the unconventionality is the proximity to first-order transitions. The border between the first-order and continuous transitions is described by a quantum tricritical point (QTCP) for symmetry-breaking transitions. One of the characteristic features is the concomitant divergence of order-parameter and uniform fluctuations in contrast to the conventional quantum critical point (QCP). Several puzzling non-Fermi-liquid properties are referred to be accounted for as in the cases of YbRh2Si2, CeRu2Si2 and beta-YbAlB4. Another more dramatic unconventionality appears in this case for topological transitions such as metal-insulator and Lifshitz transitions. This border, the marginal quantum critical point (MQCP), belongs to an unprecedented universality class with diverging uniform fluctuations at zero temperature. The MQCP has a unique feature by a combined character of symmetry-breaking and topological transitions. The theoretical results are supported by experimental indications for V2-xCrxO3 and an organic conductor kappa-(ET)2Cu[N(CN)2]Cl. Identifying topological transitions also reveals how non-Fermi liquid appears as a phase in metals. The theory also accounts for the criticality of a metamagnetic transition in ZrZn2, by interpreting it as an interplay of Lifshitz transition and correlation effects. We discuss common underlying physics in these examples.

preprint2009arXiv

Ab initio Derivation of Low-Energy Model for $κ$-ET Type Organic Conductors

We derive effective Hubbard-type Hamiltonians of $κ$-(ET)$_2X$, using an {\em ab initio} downfolding technique, for the first time for organic conductors. They contain dispersions of the highest occupied Wannier-type molecular orbitals with the nearest neighbor transfer $t$$\sim$0.067 eV for a metal $X$=Cu(NCS)$_2$ and 0.055 eV for a Mott insulator $X$=Cu$_2$(CN)$_3$, as well as screened Coulomb interactions. It shows unexpected differences from the conventional extended Hückel results, especially much stronger onsite interaction $U$$\sim$0.8 eV ($U/t$$\sim$12-15) than the Hückel estimates ($U/t$$\sim$7-8) as well as an appreciable longer-ranged interaction. Reexamination on physics of this family of materials is required from this realistic basis.

preprint2009arXiv

Ab-initio procedure for effective models of correlated materials with entangled band structure

We present a first-principles method for deriving effective low-energy models of electrons in solids having entangled band structure. The procedure starts with dividing the Hilbert space into two subspaces, the low-energy part ("$d$ space'') and the rest of the space ("$r$ space''). The low-energy model is constructed for the $d$ space by eliminating the degrees of freedom of the $r$ space. The thus derived model contains the strength of electron correlation expressed by a partially screened Coulomb interaction, calculated in the constrained random-phase-approximation (cRPA) where screening channels within the $d$ space, $P_d$, are subtracted. One conceptual problem of this established downfolding method is that for entangled bands it is not clear how to cut out the $d$ space and how to distinguish $P_d$ from the total polarization. Here, we propose a simple procedure to overcome this difficulty. In our scheme, the $d$ subspace is cut out from the Hilbert space of the Kohn Sham eigenfunctions with the help of a procedure to construct a localized Wannier basis. The $r$ subspace is constructed as the complementary space orthogonal to the $d$ subspace. After this disentanglement, $P_d$ becomes well defined. Using the disentangled bands, the effective parameters are uniquely determined in the cRPA. The method is successfully applied to 3$d$ transition metals.

preprint2009arXiv

Spin Fluctuation Theory for Quantum Tricritical Point Arising in Proximity to First-Order Phase Transitions: Applications to Heavy-Fermion Systems, YbRh2Si2, CeRu2Si2, and beta-YbAlB4

We propose a phenomenological spin fluctuation theory for antiferromagnetic quantum tricritical point (QTCP), where the first-order phase transition changes into the continuous one at zero temperature. Under magnetic fields, ferromagnetic quantum critical fluctuations develop around the antiferromagnetic QTCP in addition to antiferromagnetic ones, which is in sharp contrast with the conventional antiferromagnetic quantum critical point. For itinerant electron systems,} we show that the temperature dependence of critical magnetic fluctuations around the QTCP are given as chiQ \propto T^{-3/2} (chi0\propto T^{-3/4}) at the antiferromagnetic ordering (ferromagnetic) wave number q=Q (q=0). The convex temperature dependence of chi0^{-1} is the characteristic feature of the QTCP, which is never seen in the conventional spin fluctuation theory. We propose that the general theory of quantum tricriticality that has nothing to do with the specific Kondo physics itself, solves puzzles of quantum criticalities widely observed in heavy-fermion systems such as YbRh2Si2, CeRu2Si2, and beta-YbAlB4. For YbRh2Si2, our theory successfully reproduces quantitative behaviors of the experimental ferromagnetic susceptibility and the magnetization curve by choosing the phenomenological parameters properly. The quantum tricriticality is also consistent with singularities of other physical properties such as specific heat, nuclear magnetic relaxation time 1/T_1T, and Hall coefficient. For CeRu2Si2 and beta-YbAlB4, we point out that the quantum tricriticality is a possible origin of the anomalous diverging enhancement of the uniform susceptibility observed in these materials.

preprint2009arXiv

Systematic Control of Carrier Doping without Disorder at Interface of Oxide Heterostructures

We propose a method to systematically control carrier densities at the interface of transition-metal oxide heterostructures without introducing disorders. By inserting non-polar layers sandwiched by polar layers, continuous carrier doping into the interface can be realized. This method enables us to control the total carrier densities per unit cell systematically up to high values of the order unity.

preprint2009arXiv

Transition from Band insulator to Bose-Einstein Condensate superfluid and Mott State of Cold Fermi Gases with Multiband Effects in Optical Lattices

We study two models realized by two-component Fermi gases loaded in optical lattices. We clarify that multi-band effects inevitably caused by the optical lattices generate a rich structure, when the systems crossover from the region of weakly bound molecular bosons to the region of strongly bound atomic bosons. Here the crossover can be controlled by attractive fermion interaction. One of the present models is a case with attractive fermion interaction, where an insulator-superfluid transition takes place. The transition is characterized as the transition between a band insulator and a Bose-Einstein condensate (BEC) superfluid state. Differing from the conventional BCS superfluid transition, this transition shows unconventional properties. In contrast to the one particle excitation gap scaled by the superfluid order parameter in the conventional BCS transition, because of the multi-band effects, a large gap of one-particle density of states is retained all through the transition although the superfluid order grows continuously from zero. A reentrant transition with lowering temperature is another unconventionality. The other model is the case with coexisting attractive and repulsive interactions. Within a mean field treatment, we find a new insulating state, an orbital ordered insulator. This insulator is one candidate for the Mott insulator of molecular bosons and is the first example that the orbital internal degrees of freedom of molecular bosons appears explicitly. Besides the emergence of a new phase, a coexisting phase also appears where superfluidity and an orbital order coexist just by doping holes or particles. The insulating and superfluid particles show differentiation in momentum space as in the high-Tc cuprate superconductors.

preprint2006arXiv

Quantum criticality around metal-insulator transitions of strongly correlated electrons

Quantum criticality of metal-insulator transitions in correlated electron systems is shownto belong to an unconventional universality class with violation of Ginzburg-Landau-Wilson(GLW) scheme formulated for symmetry breaking transitions. This unconventionality arises from an emergent character of the quantum critical point, which appears at the marginal point between the Ising-type symmetry breaking at nonzero temperatures and the topological transition of the Fermi surface at zero temperature. We show that Hartree-Fock approximations of an extended Hubbard model on square latticesare capable of such metal-insulator transitions with unusual criticality under a preexisting symmetry breaking. The obtained universality is consistent with the scaling theory formulated for Mott transition and with a number of numerical results beyond the mean-field level, implying that the preexisting symmetry breaking is not necessarily required for the emergence of this unconventional universality. Examinations of fluctuation effects indicate that the obtained critical exponents remain essentially exact beyond the mean-field level. Detailed analyses on the criticality, containing diverging carrier density fluctuations around the marginal quantum critical point, are presented from microscopic calculations and reveal the nature as quantum critical "opalescence". Analyses on crossovers between GLW type at nonzero temperature and topological type at zero temperature show that the critical exponents observed in (V,Cr)2O3 and kappa-ET-type organic conductor provide us with evidences for the existence of the present marginal quantum criticality.