Search

We introduce sparse-ir, a collection of libraries to efficiently handle imaginary-time propagators, a central object in finite-temperature quantum many-body calculations. We leverage two concepts: firstly, the intermediate representation (IR), an optimal compression of the propagator with robust a-priori error estimates, and secondly, sparse sampling, near-optimal grids in imaginary time and imaginary frequency from which the propagator can be reconstructed and on which diagrammatic equations can be solved. IR and sparse sampling are packaged into stand-alone, easy-to-use Python, Julia and Fortran libraries, which can readily be included into existing software. We also include an extensive set of sample codes showcasing the library for typical many-body and ab initio methods., 8 pages, 4 figures

Analytic continuation (AC) from the imaginary-time Green's function to the spectral function is a crucial process for numerical studies of the dynamical properties of quantum many-body systems. This process, however, is an ill-posed problem; that is, the obtained spectrum is unstable against the noise of the Green's function. Though several numerical methods have been developed, each of them has its own advantages and disadvantages. The sparse modeling (SpM) AC method, for example, is robust against the noise of the Green's function but suffers from unphysical oscillations in the low-energy region. We propose a new method that combines the SpM AC with the Pad\'{e} approximation. This combination, called SpM-Pad\'{e}, inherits robustness against noise from SpM and low-energy accuracy from Pad\'{e}, compensating for the disadvantages of each. We demonstrate that the SpM- Pad\'{e} method yields low-variance and low-biased results with almost the same computational cost as that of the SpM method., Comment: 8 pages, 6 figures (proofread; changed text color, fixed some words)

We present a new open-source program, DCore, that implements dynamical mean-field theory (DMFT). DCore features a user-friendly interface based on text and HDF5 files. It allows DMFT calculations of tight-binding models to be performed on predefined lattices as well as \textit{ab initio} models constructed by external density functional theory codes through the Wannier90 package. Furthermore, DCore provides interfaces to many advanced quantum impurity solvers such as quantum Monte Carlo and exact diagonalization solvers. This paper details the structure and usage of DCore and shows some applications., Comment: 33 pages in single column format. Based on DCore v3.0.0

This lecture note reviews recently proposed sparse-modeling approaches for efficient ab initio many-body calculations based on the data compression of Green's functions. The sparse-modeling techniques are based on a compact orthogonal basis, an intermediate representation (IR) basis, for imaginary-time and Matsubara Green's functions. A sparse sampling method based on the IR basis enables solving diagrammatic equations efficiently. We describe the basic properties of the IR basis, the sparse sampling method and its applications to ab initio calculations based on the GW approximation and the Migdal-Eliashberg theory. We also describe a numerical library for the IR basis and the sparse sampling method, sparse-ir, and provide its sample codes. This lecture note follows the Japanese review article with major revisions [H. Shinaoka et al., Solid State Physics 56(6), 301 (2021)]., Comment: 29 pages, 11 figures

In this work, we investigate models for bulk, bi- and multilayers containing half-metallic ferromagnets (HMFs), at zero and at finite temperature, in order to elucidate the effects of strong electronic correlations on the spectral properties (density of states). Our focus is on the evolution of the finite-temperature many-body induced tails in the half-metallic gap. To this end, the dynamical mean-field theory (DMFT) is employed. For the bulk, a Bethe lattice model is solved using a matrix product states based impurity solver at zero temperature and a continuous-time quantum Monte Carlo (CT-QMC) solver at finite temperature. We demonstrate numerically, in agreement with the analytical result, that the tails vanish at the Fermi level at zero temperature. In order to study multilayers, taken to be square lattices within the layers, we use the real-space DMFT extension with the CT-QMC impurity solver. For bilayers formed by the HMF with a band or correlated insulator, we find that charge fluctuations between the layers enhance the finite temperature tails. In addition, in the presence of inter-layer hopping, a coherent quasiparticle peak forms in the otherwise correlated insulator. In the multilayer heterostructure setup, we find that by suitably choosing the model parameters, the tails at the HMF/Mott insulator interface can be reduced significantly, and that a high spin polarization is conceivable, even in the presence of long-ranged electrostatic interactions., Comment: to be published in Physical Review Research

We formulate an efficient scheme to perform Migdal-Eliashberg calculation considering the retardation effect from first principles. While the conventional approach requires a huge number of Matsubara frequencies, we show that the intermediate representation of the Green's function [H. Shinaoka et al., Phys. Rev. B 96, 035147 (2017)] dramatically reduces the numerical cost to solve the linearized gap equation. Without introducing any empirical parameter, we demonstrate that we can successfully reproduce the experimental superconducting transition temperature of elemental Nb ($\sim 10$ K) very accurately. The present result indicates that our approach has a superior performance for many superconductors for which $T_{\rm c}$ is lower than ${\mathcal O}(10)$ K, Comment: 8 pages, 4 figures

Many-body calculations at the two-particle level require a compact representation of two-particle Green's functions. In this paper, we introduce a sparse sampling scheme in the Matsubara frequency domain as well as a tensor network representation for two-particle Green's functions. The sparse sampling is based on the intermediate representation basis and allows an accurate extraction of the generalized susceptibility from a reduced set of Matsubara frequencies. The tensor network representation provides a system independent way to compress the information carried by two-particle Green's functions. We demonstrate efficiency of the present scheme for calculations of static and dynamic susceptibilities in single- and two-band Hubbard models in the framework of dynamical mean-field theory., 27 pages in single column format, 12 pages (added missing references)

Computing momentum-dependent susceptibilities in the dynamical mean-field theory (DMFT) requires solving the Bethe-Salpeter equation, which demands large computational cost. Exploiting the strong-coupling feature of local fluctuations, we derive a simplified formula that can be solved at a considerably lower cost. The validity and the physical meaning of the formula are confirmed by deriving the effective intersite interactions in the strong-coupling limit, such as the kinetic exchange and RKKY interactions. Furthermore, numerical calculations for single-orbital and multiorbital models demonstrate surprisingly wider applicability including weak-coupling region. Based on this formula, we propose three levels of practical approximations that can be chosen depending on complexity of problems. Simpler evaluations of spin and orbital susceptibilities in multiorbital systems thus become possible within DMFT., Comment: 18 pages, 20 figures

Optical responses in an excitonic insulating (EI) system with strong electron correlation are studied. We adopt the two-orbital Hubbard model with a finite energy difference between the two orbitals where the spin state degree of freedom exists. This model is analyzed by the variational cluster approach. In order to include the local electron correlation effect, the vertex correction is taken into account in the formulation of the optical conductivity spectra. We calculate a finite-temperature phase diagram, in which an EI phase appears between a low-spin band insulating state and a high-spin Mott insulating state. Characteristic components of the optical conductivity spectra consisting of a sharp peak and continuum appear in the EI phase. Integrated intensity almost follows the order parameter of the EI state, suggesting that this component is available to identify the EI phases and transitions., Comment: 9 pages, 11 figures

The imaginary-time Green's function is a building block of various numerical methods for correlated electron systems. Recently, it was shown that a model-independent compact orthogonal representation of the Green's function can be constructed by decomposing its spectral representation. We investigate the performance of this so-called \textit{intermedaite representation} (IR) from several points of view. First, for two simple models, we study the number of coefficients necessary to achieve a given tolerance in expanding the Green's function. We show that the number of coefficients grows only as $O(\log \beta)$ for fermions, and converges to a constant for bosons as temperature $T=1/\beta$ decreases. Second, we show that this remarkable feature is ascribed to the properties of the physically constructed basis functions. The fermionic basis functions have features in the spectrum whose width is scaled as $O(T)$, which are consistent with the low-$T$ properties of quasiparticles in a Fermi liquid state. On the other hand, the properties of the bosonic basis functions are consistent with those of spin/orbital susceptibilities at low $T$. These results demonstrate the potential wide application of the IR to calculations of correlated systems., Comment: 12 pages, 12 figures

In the present paper, we present an efficient continuous-time quantum Monte Carlo impurity solver with high acceptance rate at low temperature for multi-orbital quantum impurity models with general interaction. In this hybridization expansion impurity solver, the imaginary time evolution operator for the high energy multiplets, which decays very rapidly with the imaginary time, is approximated by a probability normalized $\delta$-function. As the result, the virtual charge fluctuations of $f^{n}\rightarrow f^{n\pm1}$ are well included on the same footing without applying Schrieffer-Wolff transformation explicitly. As benchmarks, our algorithm perfectly reproduces the results for both Coqblin-Schriffeer and Kondo lattice models obtained by CT-J method developed by Otsuki {\it et al}. Furthermore, it allows capturing low energy physics of heavy-fermion materials directly without fitting the exchange coupling $J$ in the Kondo model., Comment: 14 pages, 9 figures

The open-source library, irbasis, provides easy-to-use tools for two sets of orthogonal functions named intermediate representation (IR). The IR basis enables a compact representation of the Matsubara Green's function and efficient calculations of quantum models. The IR basis functions are defined as the solution of an integral equation whose analytical solution is not available for this moment. The library consists of a database of pre-computed high-precision numerical solutions and computational code for evaluating the functions from the database. This paper describes technical details and demonstrates how to use the library., Comment: 12 pages, 5 figures

Two-particle Green's functions and the vertex functions play a critical role in theoretical frameworks for describing strongly correlated electron systems. However, numerical calculations at two-particle level often suffer from large computation time and massive memory consumption. We derive a general expansion formula for the two-particle Green's functions in terms of an overcomplete representation based on the recently proposed "intermediate representation" basis. The expansion formula is obtained by decomposing the spectral representation of the two-particle Green's function. We demonstrate that the expansion coefficients decay exponentially, while all high-frequency and long-tail structures in the Matsubara-frequency domain are retained. This representation therefore enables efficient treatment of two-particle quantities and opens a route to the application of modern many-body theories to realistic strongly correlated electron systems., Comment: 17 pages, 12 figures (minor modifications)

Superconductivity in the two-dimensional Kondo lattice model is derived by means of the dual-fermion approach. This framework takes full account of local correlations based on the dynamical mean-field theory and further includes influence of critical antiferromagnetic fluctuations. It is found that pairing with dx2 −y2 symmetry similar to the Hubbard model appears only in the weak-coupling regime. As the coupling is increased, we observe transition to p-wave spin-singlet superconductivity. The emergence of the unconventional p-wave pairing is ascribed to formation of large Fermi surface, which embodies itinerancy of f electrons.

This review paper describes the basic concept and technical details of sparse modeling and its applications to quantum many-body problems. Sparse modeling refers to methodologies for finding a small number of relevant parameters that well explain a given dataset. This concept reminds us physics, where the goal is to find a small number of physical laws that are hidden behind complicated phenomena. Sparse modeling extends the target of physics from natural phenomena to data, and may be interpreted as "physics for data". The first half of this review introduces sparse modeling for physicists. It is assumed that readers have physics background but no expertise in data science. The second half reviews applications. Matsubara Green's function, which plays a central role in descriptions of correlated systems, has been found to be sparse, meaning that it contains little information. This leads to (i) a new method for solving the ill-conditioned inverse problem for analytical continuation, and (ii) a highly compact representation of Matsubara Green's function, which enables efficient calculations for quantum many-body systems., Comment: 24 pages, 32 figures, to be published in J. Phys. Soc. Jpn. as an invited review paper

A new approach of solving the ill-conditioned inverse problem for analytical continuation is proposed. The root of the problem lies in the fact that even tiny noise of imaginary-time input data has a serious impact on the inferred real-frequency spectra. By means of a modern regularization technique, we eliminate redundant degrees of freedom that essentially carry the noise, leaving only relevant information unaffected by the noise. The resultant spectrum is represented with minimal bases and thus a stable analytical continuation is achieved. This framework further provides a tool for analyzing to what extent the Monte Carlo data need to be accurate to resolve details of an expected spectral function., 7 pages, 5 figures

Extensions of dynamical-mean-field-theory (DMFT) make use of quantum impurity models as non-perturbative and exactly solvable reference systems which are essential to treat the strong electronic correlations. Through the introduction of retarded interactions on the impurity, these approximations can be made two-particle self-consistent. This is of interest for the Hubbard model, because it allows to suppress the antiferromagnetic phase transition in two-dimensions in accordance with the Mermin-Wagner theorem, and to include the effects of bosonic fluctuations. For a physically sound description of the latter, the approximation should be conserving. In this paper we show that the mutual requirements of two-particle self-consistency and conservation lead to fundamental problems. For an approximation that is two-particle self-consistent in the charge- and longitudinal spin channel, the double occupancy of the lattice and the impurity are no longer consistent when computed from single-particle properties. For the case of self-consistency in the charge- and longitudinal as well as transversal spin channels, these requirements are even mutually exclusive so that no conserving approximation can exist. We illustrate these findings for a two-particle self-consistent and conserving DMFT approximation., 17 pages, 9 figures

New model-independent compact representations of imaginary-time data are presented in terms of the intermediate representation (IR) of analytical continuation. This is motivated by a recent numerical finding by the authors [J. Otsuki et al., arXiv:1702.03056]. We demonstrate the efficiency of the IR through continuous-time quantum Monte Carlo calculations of an Anderson impurity model. We find that the IR yields a significantly compact form of various types of correlation functions. The present framework will provide general ways to boost the power of cutting-edge diagrammatic/quantum Monte Carlo treatments of many-body systems., Comment: 9 pages, 7 figures

A quantum critical behavior of the magnetic susceptibility was observed in a quasicrystal containing ytterbium. At the same time, a mixed-valence feature of Yb ions was reported, which appears to be incompatible with the magnetic instability. We derive the magnetic susceptibility by expressing the quasiperiodicity as the distributed hybridization strength between Yb 4f and conduction electrons. Assuming a wide distribution of the hybridization strength, the most f electrons behave as renormalized paramagnetic states in the Kondo or mixed-valence regime, but a small number of f moments remain unscreened. As a result, the bulk magnetic susceptibility exhibits a nontrivial power-law-like behavior, while the average f-electron occupation is that of mixed-valence systems. This model thus resolves two contradictory properties of Yb quasicrystals., 4 pages, 5 figures

A modified version of the spinless Anderson model is studied by means of the continuous-time quantum Monte Carlo method. This study is motivated by the peculiar heavy-fermion behavior observed in certain Samarium compounds, which is insensitive to magnetic field. The model involves M channels for conduction electrons, all of which interact with local f electron via the Coulomb repulsion Ufc, while only one channel has hybridization with the local state. The effective hybridization is reduced by the Anderson orthogonality effect, and a quantum critical point occurs with increasing M and/or increasing Ufc. The numerical results at finite temperature of the local charge susceptibility are well fitted by a simple scaling theory for all M. However, the single-particle spectrum is described by a double Lorentzian for M>1, in contrast with the single Lorentzian with M=1. A quasi-particle perturbation theory is presented that reproduces the quantum critical point for large M. The quasi-particle theory gives not only the renormalized energy scale, but its extrapolation toward higher energies being consistent with the double Lorentzian spectrum., 10 pages, 10 figures

The Kondo lattice model describes a quantum phase transition between the antiferromagnetic state and heavy-fermion states. Applying the dual-fermion approach, we explore possible superconductivities emerging due to the critical antiferromagnetic fluctuations. The d-wave pairing is found to be the leading instability only in the weak-coupling regime. As the coupling is increased, we observe a change of the pairing symmetry into a p-wave spin-singlet pairing. The competing superconductivities are ascribed to crossover between small and large Fermi surfaces, which occurs with the formation of heavy quasiparticles.

The dual-fermion approach offers a way to perform diagrammatic expansion around the dynamical mean-field theory. Using this formalism, the influence of antiferromagnetic fluctuations on the self-energy is taken into account through ladder-type diagrams in the particle-hole channel. The resulting phase diagram for the (quasi-)two-dimensional Hubbard model exhibits antiferromagnetism and d-wave superconductivity. Furthermore, a uniform charge instability, i.e., phase separation, is obtained in the low doping regime around the Mott insulator. We also examine spin/charge density wave fluctuations including d-wave symmetry. The model exhibits a tendency towards an unconventional charge density-wave, but no divergence of the susceptibility is found., Comment: 14 pages, 13 figures, 2 tables

To explore the nature of the metallic state near the transition to a Mott insulator, we investigate the $t\mathrm{\text{\ensuremath{-}}}J$ model with random exchange interaction in $d=\ensuremath{\infty}$ dimensions. A numerically exact solution is obtained by an extension of the continuous-time quantum Monte Carlo method to the case of a vector bosonic field coupled to a local spin. We show that the paramagnetic solution near the Mott insulator describes an incoherent metal with a residual moment, and that single-particle excitations produce an additional band, which is separated from the Mott-Hubbard band.

A vector bosonic field coupled to the electronic spin is treated by means of the continuous-time quantum Monte Carlo method. In the Bose Kondo model with a sub-Ohmic density of states $\rho_{B}(\omega) \propto \omega^{s}$ with s=0.2, two contributions to the spin susceptibility, the Curie term T^{-1} and the term T^{-s} due to bosonic fluctuations, are observed separately. This result indicates the existence of a residual moment and a hidden critical behavior. By including hybridization with itinerant electrons, a quantum critical point is identified between this local-moment state and the Kondo singlet state. It is demonstrated that the energy scale of the bosonic fluctuations is not affected by the quantum phase transition., Comment: 8 pages, 7 figures

The continuous-time quantum Monte Carlo method is applied to the interacting resonant level model (IRLM) using double expansion with respect to Coulomb interaction Ufc and hybridization V. Thermodynamics of the IRLM without spin is equivalent to the anisotropic Kondo model in the low-energy limit. Exact dynamics and thermodynamics of the IRLM are derived numerically for a wide range of Ufc with a given value of V. For negative Ufc, excellent agreement including a quantum critical point is found with a simple scaling formula that deals with V in the lowest-order, and Ufc up to infinite order. As Ufc becomes positive and large, lower order scaling results deviate from exact numerical results. Possible relevance of the results is discussed to certain Samarium compounds with unusual heavy-fermion behavior., Comment: 9 pages, 8 figures

Peculiar property of electronic order is clarified for the two-channel Kondo lattice. With two conduction electrons per site, the order parameter is a composite quantity involving both local and itinerant degrees of freedom. In contrast to the ordinary Kondo lattice, a heavy electron band is absent above the transition temperature, but is rapidly formed below it. The change of entropy associated with the ordering is found to be close to $\ln \sqrt 2$ per site. This entropy corresponds to the residual entropy in a two-channel Kondo impurity, which has been regarded as due to localized free Majorana particles. The present composite order is interpreted as instability of Majorana particles toward non-Kramers conduction electrons plus heavy fermions that involve localized electrons., Comment: 9 pages, 11 figures

To explore the nature of the metallic state near the transition to a Mott insulator, we investigate the t-J model with random exchange interaction in d=∞ dimensions. A numerically exact solution is obtained by an extension of the continuous-time quantum Monte Carlo method to the case of a vector bosonic field coupled to a local spin. We show that the paramagnetic solution near the Mott insulator describes an incoherent metal with a residual moment, and that single-particle excitations produce an additional band, which is separated from the Mott-Hubbard band.

A nonperturbative approach to unconventional superconductivity is developed based on the idea of two-particle self-consistent (TPSC) theory by Vilk and Tremblay. A sum rule, which the momentum-dependent pairing susceptibility satisfies, is derived. An effective pairing interaction between quasiparticles is determined so that the susceptibility should fulfill this exact sum rule, in which fluctuations belonging to different symmetries couple at finite momentum. It is demonstrated that the mode coupling between $d$-wave and $s$-wave pairing fluctuations leads to suppression of the $d$-wave fluctuation near the Mott insulator.

We show by quantum Monte Carlo simulations of realistic Kondo lattice models derived from electronic--structure calculations that multiple quantum critical points can be realized in Plutonium--based materials. We place representative systems including PuCoGa5 on a realistic Doniach phase diagram and identify the regions where the magnetically mediated superconductivity could occur. Solution of an inverse problem to restore the quasiparticle renormalization factor for f-electrons is shown to be sufficiently good to predict the trends among Sommerfeld coefficients and magnetism. Suggestion on the possible experimental verification for this scenario is given for PuAs., 4 pages, 3 figures, 2 tables, references updated

A novel type of symmetry breaking is reported for the two-channel Kondo lattice where conduction electrons have spin and orbital (channel) degrees of freedom. Using the continuous-time quantum Monte Carlo and the dynamical mean-field theory, a spontaneous breaking of the orbital symmetry is observed. The tiny breakdown of orbital occupation number, however, vanishes if the conduction electrons have the particle-hole symmetry. The proper order parameter instead is identified as a composite quantity representing the orbital-selective Kondo effect. The single-particle spectrum of the selected orbital shows insulating property, while the other orbital behaves as Fermi liquid. This composite order is the first example of odd-frequency order other than off-diagonal order (superconductivity), and is a candidate of hidden order in $f$-electron systems., Comment: 5pages, 4figures

A microscopic model is proposed for the scalar order in PrFe4P12 where f2 crystalline electric field (CEF) singlet and triplet states interact with two conduction bands. By combining the dynamical mean-field theory and the continuous-time quantum Monte Carlo, we obtain an electronic order with staggered Kondo and CEF singlets with the total conduction number being unity per site. The ground state becomes semimetallic provided that the two conduction bands have different occupation numbers. This model naturally explains experimentally observed properties in the ordered phase of PrFe4P12 such as the scalar order parameter, temperature dependence of the resistivity, field-induced staggered moment, and inelastic features in neutron scattering. The Kondo effect plays an essential role for ordering, in strong contrast with ordinary magnetic orders by the RKKY interaction., 4 pages, 4figures

We present realistic Kondo-lattice simulation results for the recently-discovered heavy-fermion antiferromagnet CePt2In7 comparing with its three-dimensional counterpart CeIn3 and the less two-dimensional ones, Ce-115’s. We find that the distance to the magnetic quantum critical point is the largest for CeIn3 and the smallest for Ce-115’s, and CePt2In7 falls in between. We argue that thetrend in quasi-two-dimensional materials stems from the frequency dependence of the hybridization between Cerium 4f-electrons and the conduction bands.

The impurity Green's function Gf in the local non-Fermi liquid state is evaluated by means of the continuous-time quantum Monte Carlo method extended to the multichannel Anderson model. For N=M (where N and M are numbers of spin components and channels, respectively), Gf is expressed as -Im Gf(w+i0) = c - b |w|^{1/2}, and the zero-frequency value c depends only on N (=M). A corresponding impurity self-energy at low frequencies is composed of two parts: a resonance term related to c, and a non-Fermi liquid term proportional to |w|^{1/2}. The characteristic energy scale is discussed in terms of the non-Fermi liquid term in the self-energy., 4 pages, 6 figures, proceedings of ICHE2010

Novel electronic order is found theoretically for a system where even number of localized electrons per site are coupled with conduction electrons. For precise quantitative study, a variant of the Kondo lattice model is taken with crystalline electric field (CEF) singlet and triplet states for each site. Using the dynamical mean-field theory combined with the continuous-time quantum Monte Carlo method, a staggered order with alternating Kondo and CEF singlets is identified for a case with one conduction electron per site being distributed in two conduction bands each of which is quarter-filled. This electronic order accompanies a charge density wave (CDW) of conduction electrons that accumulate more on Kondo-singlet sites than on CEF-singlet sites. Possible relevance of the present order to the scalar order in PrFe$_4$P$_{12}$ is discussed., Comment: 11 pages, 17 figures

We propose a new framework for first-principle calculations of heavy-fermion materials. These are described in terms of the Kondo lattice Hamiltonian with the parameters extracted from a realistic density functional based calculation which is then solved using continuous-time quantum Monte Carlo method and dynamical mean field theory. As an example, we show our results for the Neel temperatures of Cerium-122 compounds (CeX2Si2 with X=Ru, Rh, Pd, Cu, Ag, and Au) where the general trend around the magnetic quantum critical point is successfully reproduced. Our results are organized on a universal Doniach phase diagram in a semi-quantitative way., Comment: revised to fit in 4 pages, to appear in Phys. Rev. Lett

Single-particle spectrum of the Kondo lattice model is derived with use of the continuous-time quantum Monte Carlo method, combined with the dynamical mean-field theory. Crossover behavior is traced quantitatively either to a heavy Fermi-liquid state or to a magnetically ordered state from the local-moment state at high temperatures. The momentum distribution in the low-temperature limit acquires a discontinuity at the location that involves the local-spin degrees of freedom. Even without the charge degrees of freedom for local electrons, the excitation spectra exhibit hybridized bands similar to those in the Anderson lattice. Temperature dependence in the zero-energy component of the self-energy is crucial in forming the Fermi-liquid state with the large Fermi surface., 4 pages, 4 figures, to appear in Phys. Rev. Lett

Highly accurate numerical results for single-particle spectrum and order parameter are obtained for the magnetically ordered Kondo lattice by means of the dynamical mean-field theory combined with the continuous-time quantum Monte Carlo method. Hybridized energy bands involving local spins are identified in the N\'{e}el state as a hallmark of itinerant antiferromagnetism. At the boundary of the reduced Brillouin zone, the two-fold degeneracy remains in spite of the doubled unit cell. This degeneracy results if the molecular field felt by localized spins has identical magnitude and reversed direction with that of conduction electrons. The persistent Kondo effect is responsible for the behavior. The antiferromagnetic quantum transition occurs inside the itinerant regime, and does not accompany the itinerant-localized transition., Comment: 4 pages, 5 figures

Dynamical properties are studied numerically for a variant of the Kondo model with singlet and triplet crystalline electric field (CEF) levels where Kondo and CEF singlets compete for the ground state. Using the continuous-time quantum Monte Carlo method, we derive the $t$-matrix of conduction electrons and dynamical susceptibilities of local electrons without encountering the negative sign problem. When the CEF splitting is comparable to the Kondo temperature, the dynamical response has only a quasi-elastic peak. Nevertheless, the local single-particle spectrum shows an energy gap in strong contrast with the ordinary Kondo model., Comment: 9 pages, 10 figures

Magnetic and charge susceptibilities in the Kondo lattice are derived by the continuous-time quantum Monte Carlo (CT-QMC) method combined with the dynamical mean-field theory. For a weak exchange coupling J and near half filling of the conduction band, antiferromagnetic transition occurs as signalled by divergence of the staggered magnetic susceptibility with lowering temperature. With increasing J, the Kondo effect suppresses the divergence, and the critical value of J agrees well with Doniach's estimate which considers the RKKY interaction as competing with the Kondo effect. For low density of conduction electrons, a ferromagnetic ordering is observed where Doniach's estimate does not work. Around quarter filling, a charge-density-wave (CDW) transition is found. The CDW is interpreted from the strong-coupling limit in terms of effective repulsion between Kondo singlets., 10 pages, 16 figures

A Green-function formalism for the Kondo lattice model is presented, which is designed to be combined with the dynamical mean-field theory. With use of Wick's theorem only for conduction electrons, dynamical quantities are represented in terms of the t-matrix and its two-particle generalizations. By taking the high-frequency limit of the t-matrix with respect to a part of the fermion frequencies, one obtains dynamical correlation functions of localized moments. Explicit examples of calculational steps are provided with use of the continuous-time quantum Monte Carlo method., 13 pages, 10 figures

An impurity solver based on a continuous-time quantum Monte Carlo method is developed for the Coqblin-Schrieffer model. The Monte Carlo simulation does not encounter a sign problem for antiferromagnetic interactions, and accurately reproduces the Kondo effect. Our algorithm can deal with an arbitrary number N of local degrees of freedom, becomes more efficient for larger values of N, and is hence suitable for models with orbital degeneracy. The dynamical susceptibility and the impurity t-matrix are derived with the aid of the Pad\'e approximation for various values of N, and good agreement is found with other methods and available exact results. We point out that the Korringa-Shiba relation needs correction for a finite value of the exchange interaction., Comment: 12 pages, 15 figures

Characteristics of hybridization and multipole orders of $4f$-electrons in Pr skutterudites are explained in terms of a pseudo-quartet composed of crystalline electric field (CEF) singlet and the triplet. It is shown that the contrasting behaviors observed in PrFe$_4$P$_{12}$ and PrOs$_4$Sb$_{12}$ are ascribed to the difference in triplet wave functions and the CEF splittings.Since a macroscopic degeneracy remains even in the ordered phase with the $\Gamma_3$-type antiferro-quadrupole (AFQ) order,a model with strong AFQ fluctuation and static monopole/hexadecapole order is proposed for PrFe$_4$P$_{12}$.The remaining cubic symmetry explains qualitatively the behavior of staggered magnetization observed by NMR and neutron scattering. For identification of possible hexadecapole orders in PrFe$_4$P$_{12}$ andPrRu$_4$P$_{12}$, eightfold intensity pattern is predicted in the azimuthal angle scan of resonant X-ray scattering. In PrOs$_4$Sb$_{12}$, the ferromagnetic exchange coupling with the conduction band does not lead to magnetic Kondo effect, but a momentum-dependent quadrupole coupling can give rise to enhanced effective mass, which should be sensitive to disorder., Comment: Invited paper to be published in Proc. Int. Symp. on Advances in the Physics and Chemistry of Actinide Compounds (Suppl. J. Phys. Soc. Jpn)