## Valence transition in topological Kondo insulator

Jia-Tao Zhuang, 1 Xiao-Jun Zheng, 1 Zhi-Yong Wang, 1 Xing Ming, 1 Huan Li, 1, ∗ Yu Liu, 2,3, † and Hai-Feng Song 2,3

1

College of Science, Guilin University of Technology, Guilin 541004, China 2 Institute of Applied Physics and Computational Mathematics, Beijing 100088, China 3 Software Center for High Performance Numerical Simulation, China Academy of Engineering Physics, Beijing 100088, China

(Dated: August 5, 2019)

We investigate the valence transition in three-dimensional topological Kondo insulator through slave-boson analysis of periodic Anderson model. By including the effect of intra-atomic Coulomb correlation U fc between conduction and local electrons, we find a first-order valence transition from Kondo region to mixed valence upon ascending of local level above a critical U fc , and this valence transition usually occurs very close to or simultaneously with a topological transition. Near the parameter region of zero-temperature valence transition, rise of temperature can generate a thermal valence transition from mixed valence to Kondo region, accompanied by a first-order topological transition. Remarkably, above a critical U fc which is considerable smaller than that generating paramagnetic valence transition, the original continuous antiferromagnetic transition is shifted to first order one, at which a discontinuous valence shift takes place. Upon increased U fc , the paramagnetic valence transition approaches then converges with the first-order antiferromagnetic transition, leaving an significant valence shift on the magnetic boundary. The continuous antiferromagnetic transition, first-order antiferromagnetic transition, paramagnetic valence transition and topological transitions are all summarized in a global phase diagram. Our proposed exotic transition processes can help to understand the thermal valence variation as well as the valence shift around the pressure-induced magnetic transition in topological Kondo insulator candidates and in other heavy-fermion systems.

PACS numbers: 75.30.Mb, 75.70.Tj, 75.30.Kz

## I. INTRODUCTION

Since the proposal of topological intrinsic in some special Kondo insulators (KI) 1,2 , renewed attention has been attracted on these old materials, now are known as 'topological Kondo insulators'(TKI) represented by SmB 6 3-18 . In these TKIs, the strong spin-orbit coupled hybridization between d and f electrons guarantees time-reversal symmetry (TRS) and generates a band inversion between d and f orbits at certain time-reversal-invariant momenta (TRIM), leading to a topologically protected state classified by Z 2 invariants 1,2,4,19,20 . Via the spin- and angle-resolved photoemission spectroscopy (SARPES), the metallic surface states with Dirac points have been observed in SmB , confirming its topological character6 istic 11,21 . A Variety of theoretical works have been carried out for TKI to reveal its rich topological phases, topological transitions, surface states, and magnetic transitions, through first-principle calculations as well as model calculations for periodic Anderson model (PAM) 3,5-9,12,22-24 .

Recently, the high-pressure studies of SmB 6 demonstrate a magnetic transition at 6 ∼ 8 GPa 25-31 , around which the mean valence ν of Sm ion (Sm ν + ) displays a quite rapid variation from 2.5 to saturated value 3 31 , in addition, the mean valence also increases as temperature rises 32,33 , similar to other TKI candidates such as YbB 12 and pressured golden SmS 34-36 . In the hole representation for the f shell (in the filled 4 f 6 base) of Sm atom, the valence ν is related to the mean occupation of f holes n f by ν = 2 + n f 4,26 , so the valence shift of Sm ion indicates a variation of f -occupation with pressure, from mixed valence (MV), to Kondo region (or local moment region) in which the f electrons are nearly localized to create a magnetic order. In this context, the possible valence transi- tion or valence crossover in TKI, as well as the relation to the magnetic transition deserve further theoretical investigation.

For Ce- and Yb- based heavy-fermion compounds, the firstorder valence transition (FOVT) was observed decades ago in pressure studies 37-39 , in which the valences of Ce and Yb ion increase abruptly at a critical temperature, and under enhanced pressure, FOVT can be suppressed and terminates at a critical end point (CEP) to become a valence crossover. The valence of Ce (Yb) is manipulated by the electron (hole) occupation number n f of the 4 f shell 40 , so the discontinuous valence jump indicates a first-order transition from MV (with small n f ) to Kondo region (nearly localized f electrons with n f ∼ 1 ). Pressure can drive magnetic transitions in some heavy-fermion compounds, and it's found that FOVT reaches the magnetic boundary (e.g. in YbInCu 4 38,41 ), implying a strong interplay between FOVT and magnetic transition. The zero-temperature FOVT in heavy-fermion systems can be understood through PAM by considering the on-site Coulomb repulsion U fc between local f and conducting c electrons 42,43 . Upon ascent of the energy level /epsilon1 f of f orbit by strengthened pressure, the influence of U fc causes a much rapider ascent of the renormalized f level (relative to the chemical potential), pours the electrons into conduction band then consequently drives a decrease of f occupation number n f , leading to the crossover behavior of the valence. Above a critical U fc , the valence crossover is strengthened then finally changed into FOVT from Kondo region to MV, showing a abrupt fall of n f at a critical /epsilon1 f 40,44 .

Beside /epsilon1 f , pressure applied on heavy-fermion compounds can lead to variations of other model parameters such as the hybridization strength, which all can affect the FOVT 40 , so the experimentally observed FOVT is a combined outcome

within model description. In YbInCu , the enhanced pres4 sure drives a FOVT from MV to Kondo region 38 , contrary to Ce-based systems, therefore the pressure-induced FOVT cannot be simply attributed to the ascent of /epsilon1 f 42 . By contrast, the FOVT from MV to Kondo region in Ce- and Yb- compounds by rising temperature is a purer effect thus can be interpreted more straightforwardly, but such thermal FOVT is still lacking of theoretical investigation. On the other hand, for SmB 6 and TKI candidate golden SmS, near the pressure-driven magnetic boundary, the valence of Sm ion shows an active increase 30,31,45 , similarly, at low temperatures, the pressure-induced magnetic transition in YbInCu 4 holds simultaneously a FOVT 38 , therefore, the relation of FOVT or valence crossover in TKI (and also in other heavy-fermion systems) to the magnetic transition should be clarified in an unified framework. Particularly, in TKI, variation of model parameters can produce various topological phases and distinct transition processes among them, and can driven magnetic transition as well 46,47 , so interest questions arise, as how the FOVT or valence crossover appears in TKI? what is the relation between FOVT, topological transitions and magnetic transition? can any unrevealed novel transition process takes place in TKI?

This work is devoted to answer above questions. We first study the zero-temperature valence transition in threedimensional TKI, modeling by spin-orbit coupled PAM with U fc interactions. Similar to other typical KIs, we find the rise of f energy level /epsilon1 f can induce a rapid decrease of f electron number n f , causing a valence crossover of f orbit, then above a critical U fc , the valence crossover is shifted to a FOVT from Kondo region to MV. We also find that the U fc interaction has an significant impact on the topological boundaries of TKI. Remarkably, FOVT generally takes place very close to a topological boundary, and with slightly greater U fc than the critical one, FOVT can simultaneously cause a first-order topological transition, which should be clarified more rigorously in future studies. We also propose a thermal FOVT of TKI from MV to Kondo region by rising temperature in a narrow parameter regime, simultaneously with a first-order topological transition. Furthermore, we find the effect of U fc interaction leads to a strong variation of f valence near the continuous antiferromagnetic (AF) transition in TKI, and further enhancement of U fc can drive the AF transition from continuous one to first order one, meanwhile, the continuous valence variation is shifted to a FOVT at the AF boundary. In a narrow U fc window, we find a gradual approaching then convergence between paramagnetic (PM) FOVT and FOVT-associated firstorder AF transition on U fc -/epsilon1 f plane. In addition, such FOVTassociated first-order AF transition can also be generated by increased temperature in some parameter regime. The PM FOVT, two classes of AF transitions and the topological transitions are all summarized in a global phase diagram.

Our work provides new insight into the exotic transition precesses in TKI, and can be used to qualitatively understand the observed thermal valence variation as well as the valence change around the magnetic transition in SmB 6 , other TKIs and heavy-fermion systems.

TABLE I: Two sets of EHA used in this work

|         |   t d |   t ′ d |   t ′′ d |   t f |   t ′ f |   t ′′ f |
|---------|-------|---------|----------|-------|---------|----------|
| EHA(I)  |     1 |   0.15  |    0     |  -0.2 |   -0.02 |     0    |
| EHA(II) |     1 |  -0.375 |   -0.375 |  -0.2 |    0.09 |     0.09 |

## II. VALENCE TRANSITION AND TOPOLOGICAL TRANSITIONS

We use the spin-1/2 half-filled PAM in cubic lattice with a spin-orbit coupled c f -hybridization to describe TKI 2,5 . This model is adopted frequently in the literature and successfully reveals the topological aspects of TKI 6,12,46,48 . Besides, the on-site c f -Coulomb interaction U fc is the crucial driving force of FOVT in heavy-fermion systems thus should be included 40,42,44 . The model Hamiltonian of PAM reads:

$$& \text{lis} \text{-} \quad \text{cluded} ^ { 0, 4 2, 4 4 }. \text{The model Hamiltonian of PAM reads:} \\ & \text{ow} \quad \text{$\mathcal{H} = \sum_{i,j,\sigma}(-t_{ij}^{c}c_{ij}^{\dagger}c_{ij}^{\sigma}-t_{ij}^{f_{ij}^{\dagger}f_{ij}^{\sigma})+\epsilon_{f}\sum_{i,\sigma}f_{ij}^{\dagger}f_{ij}^{\sigma}f_{ij}^{\sigma} \\ & \text{tic} \\ & \text{kes} \quad \text{+} U \sum_{i} n_{i\uparrow} n_{i\downarrow}^{f_{ij}^{n_{ij}^{f_{ij}^{n_{ij}^{f_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{ny_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{N_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{nan_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{[n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}${n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}[[n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}
n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}`
n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}}
n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}\n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}<n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}{\n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}_{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}){n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}\n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij}^{n_{ij$$

in which half-filling of total electrons n t = n c + n f = 2 is fixed by adjusting the chemical potential µ , /vector σ is threedimensional vector formed by three Pauli matrices. We choose the electron hopping amplitude (EHA) up to nextnext-nearest-neighbor, with the two sets of EHAs shown in Tab.I, both retaining an insulating bulk gap 46 . The spin- and oriental- dependent hybridization between c and f electrons on neighboring sites linking by coordination vector /vector l guarantees the TRS of above PAM. The on-site f f -Coulomb repulsion U is general large thus is set to infinite for simplicity 4,12,46,48 , appropriate for applying the standard slave-boson mean-field technique, through which we arrive at the effective Hamiltonian in momentum space as

$$\begin{smallmatrix} p \alpha \, & \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \text{-} \, \end{smallmatrix} \\ \text{global} \, & \, H _ { \text{MF} } = N [ \lambda ( b ^ { 2 } - 1 ) - U _ { f c } n _ { c } n _ { f } ] + \\ \text{san-} \, & \, \text{$\sum_{k,\alpha,\beta}$} \, ( c _ { k \alpha } ^ { \dagger }, f _ { k \alpha } ^ { \dagger } ) \left ( \begin{smallmatrix} \tilde { \epsilon } _ { k } ^ { c } \delta _ { \alpha \beta } & \tilde { \ V S } _ { k } \cdot \vec { \sigma } _ { \alpha \beta } \\ \tilde { \ V S } _ { k } \cdot \vec { \sigma } _ { \alpha \beta } & \tilde { \epsilon } _ { k } ^ { f } \delta _ { \alpha \beta } \end{smallmatrix} \right ) \left ( \begin{smallmatrix} c _ { k \beta } \\ f _ { k \beta } \end{smallmatrix} \right ), \, \text{(2)} \\ \text{one} \, & \, \text{$\dots} \, & \, \text{$\dots} \, & \, \text{$\dots} \, & \, \text{$\dots} \, & \, \text{$\dots} \, & \, \text{$\dots} \, & \, \text{$\dots} \, & \, \text{$\dots} \, & \, \text{$\dots} \, & \, \text{$\dots} \, & \, \text{$\dots} \, & \, \text{$\dots} \, & \, \text{$\dots} \, & \, \text{$\dots} \, & \, \text{$\dots} \, & \, \text{$\dots} \, &                                                                                                                                                                                                         \end{smallmatrix}$$

in which α and β represent spin orientations, λ is the lagrange multiplier, b is the mean expectation value of slave bosons obeying the relation n f = 1 -b 2 . The effective hybridization is renormalized as ˜ V = V b . U fc term in Eq. 1 has been decoupled via Hatree-Fock approximation. The renormalized c and f dispersions are ˜ /epsilon1 c k = /epsilon1 c k + U fc n f -µ and ˜ /epsilon1 f k = b 2 /epsilon1 f k + /epsilon1 f + λ + U fc n c -µ , respectively, in which the tight-binding dispersions /epsilon1 c k and /epsilon1 f k are determined by their corresponding EHAs 46 . S k = (sin k · a 1 , sin k · a 2 , sin k · a 3 ) , 12 where a 1 , a 2 , a 3 are the element vectors of cubic lattice.

The quasi-particle dispersions are derived by diagonalizing the Hamiltonian matrix in Eq.2 (in its modified form) as E ± k = 1 2 [˜ /epsilon1 c k + ˜ /epsilon1 f k ± √ (˜ /epsilon1 c k -˜ ) /epsilon1 f k 2 +4 ˜ V 2 S 2 k ] , both doubly degenerated, then the ground state energy is expressed by

/s32

/s32

FIG. 1: (Color online) Occupation number of f orbit n f and bulk insulating gap ∆ as functions of local level /epsilon1 f . From (a) to (d), U fc =4.1, 4.31, 4.6, 5. U fc =4.31 is the critical value, above which a FOVT occurs at a certain /epsilon1 f , denoted by red dashed vertical lines. When U fc is greater than 4.38, a first-order topological transition takes place simultaneously at FOVT ((c) and (d)). Black dashed vertical lines denote the conventional topological transitions. Parameters: EHA(II) and V = 1 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00913v1-figure-6d3ae407-f220-4f20-85dd-a5ccb4c897b2.jpg)

E g = N λ b [ ( 2 -1) -U fc n n c f ] + 2 ∑ k , ± θ ( -E ± k ) E ± k . The mean-field parameters λ , b , and µ are then solved by saddle point approximation of E g via ∂E /∂µ g = -n t , ∂E /∂λ g = 0 and ∂E /∂b g = 0 , to obtain the following set of equations

$$\text{and} \, \partial \hat { E } _ { g } ^ { \mu } / \partial b = 0, \, \text{to obtain the following set of equations} \\ n _ { t } = \frac { 2 } { N } \sum _ { k, \pm } \theta ( - E _ { k } ^ { \pm } ), & \text{fig. 2} \text{ energy} \\ n _ { f } = \frac { 1 } { N } \sum _ { k, \pm } \theta ( - E _ { k } ^ { \pm } ) [ 1 \mp \frac { \tilde { \epsilon } _ { k } ^ { c } - \tilde { \epsilon } _ { k } ^ { f } } { \sqrt { ( \tilde { \epsilon } _ { k } ^ { c } - \tilde { \epsilon } _ { k } ^ { f } ) ^ { 2 } + 4 V ^ { 2 } b ^ { 2 } S _ { k } ^ { 2 } } } ], & \text{and} \, \text{a zoon} \, \text{eq} _ { f } \, \text{of} \, \text{(a) and} \\ \lambda = - \frac { 1 } { N } \sum _ { k, \pm } \theta ( - E _ { k } ^ { \pm } ) [ \tilde { e } _ { k } ^ { f } \pm \frac { 2 V ^ { 2 } S _ { k } ^ { 2 } - \tilde { e } _ { k } ^ { f } ( \tilde { \epsilon } _ { k } ^ { c } - \tilde { e } _ { k } ^ { f } ) } { \sqrt { ( \tilde { \epsilon } _ { k } ^ { c } - \tilde { e } _ { k } ^ { f } ) ^ { 2 } + 4 V ^ { 2 } b ^ { 2 } S _ { k } ^ { 2 } } } ], & \text{and} \, \\ \text{where} \, \theta ( - E _ { k } ^ { \pm } ) \text{ is an step function. } \text{ the above equations} \, \text{cendir}$$

where θ ( -E ± k ) is an step function. The above equations should be solved by numerical iteration.

In SmB , the itinerant 5d band and local 4f band both lo6 cate near the Fermi level 22 , causing the emergence of intraatomic Coulomb repulsion U fc . Similarly, in Ce and Yb systems, such interactions are also non-ignorable and play a crucial role in their valence transitions, although the magnitude of U fc may differs from each other 40 . In this sense, the pressure- or temperature-induced valence variation of SmB 6 and other TKIs 31-36 should be understood on the basis of U fc interaction. On the other hand, the pressure applied on heavyfermion systems enhances the hybridization strength V , and elevates the local f level /epsilon1 f as well, therefore, to simulate the valence transition in TKI, we first test the valence change (embodied by variation of n f ) as a function of /epsilon1 f under various magnitudes of U fc , then study the evolution of valence transition point with V .

/s32

/s32

/s102

FIG. 2: (Color online) FOVT (red solid lines) driven by shift of local energy level /epsilon1 f at presence of c f -intra-atomic Coulomb interaction U fc . Above a critical U fc =4.31, a FOVT occurs between Kondo state ( n f ∼ 1) and MV state (small n f ). Variation of U fc also remarkably shifts the topological boundaries (blue dashed lines). (b) is a zoom in of the dashed area in (a). (c) illustrates the critical U fc and /epsilon1 f of QCP as functions of hybridization V . Parameters: EHA(II); in (a) and (b), V = 1 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00913v1-figure-54360d38-b377-4f96-a933-a17f706e1f0c.jpg)

The numerical results are illustrated in Fig.1 with EHA(II) and V = 1 . With weak U fc , n f decreases gradually with ascending /epsilon1 f , showing a valence crossover behavior, see Fig.1a. While U fc is enhanced, n f shows a stronger decline. At a critical value U fc = 4 31 . , n f drops sharply at /epsilon1 f = -2 89 . , indicating the emergence of a quantum critical point (QCP), see Fig.1b. While U fc is further increased, n f shows an abrupt jump from n f ∼ 1 to n f /lessmuch 1 at a critical /epsilon1 f , generating a FOVT from Kondo region in which local moments weakly hybridize with c electrons, to MV region in which c and f electrons are strongly coupled, see Fig.1c-d. The FOVT is generated by strong U fc interaction which forces the f electrons to pour into c band when the effective f level is suddenly lifted during the ascent of /epsilon1 f .

Previous works have shown that the variation of /epsilon1 f or V can produce successive topological transitions driven by closing and reopening of bulk insulating gap in TKI 3,5,47 , and through present calculations, we find that these topological boundaries can be shifted sensitively by U fc towards lower /epsilon1 f direction,

/s50/s46/s53

/s32

FIG. 3: (Color online) Phase diagram at non-zero temperatures for (a) U fc = 4 1 . and (b) U fc = 4 4 . . Thick red lines denote FOVT, thin red lines denote continuous decoupling of heavy-fermion state, blue dashed lines denote topological transitions, red dots represent CEP of FOVT. Inset of (b) shows the thermal FOVTs among MV, Kondo and DM states. Parameters: EHA(II) and V = 1 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00913v1-figure-ef2835d5-156e-4665-8242-bd1ab8b32843.jpg)

see the black dashed lines in Fig.1a-c and blue dashed lines in Fig.2. Remarkably, attributing to the rapid variation of bulk gap around the QCP of FOVT, the QCP usually locates very close to a certain topological transition, see Fig.1b. Upon further increased U fc from the QCP, the topological boundary first approaches then finally converges with FOVT, leaving a first-order topological transition characterized by discontinuous change of insulating gap. The unconventional first-order topological transition is a direct outcome of FOVT in TKI, and should be verified more rigorously in future work.

## III. THERMALVALENCE TRANSITIONS

The QCP, FOVT and topological transitions are summarized on U fc -/epsilon1 f plane from U fc = 0 in Fig.2a. In Fig.2b, the phase diagram around the QCP is shown in enlarged pattern, in which the FOVT successively converges with the topological boundaries at two points, one at U fc = 4 38 . and /epsilon1 f = -2 93 . , the other at U fc = 4 82 . and /epsilon1 f = -3 19 . , separating the FOVT into three segments: (1) when 4.31 &lt; U fc &lt; 4.38, FOVT occurs near conventional STI ¯ ¯ Γ X -WTI ¯ X transition ( the subscripts denote the Dirac points on the surface Brillouin zone 47 ); (2) when 4.38 &lt; U fc &lt; 4.82, FOVT occurs simultaneously with first-order STI ¯ ¯ Γ X -WTI ¯ X transition; (3) when U fc &gt; 4.82, a first-order STI ¯ ¯ Γ X -STI ¯ M transition takes place at FOVT. Below the QCP, namely when U fc &lt; 4.31, U fc is insufficient to produce FOVT, but rather a valence crossover between Kondo region and MV. Therefore, whether a heavyfermion system undergoes a FOVT mainly depends on the strength of U fc , e.g., in SmB 6 and SmS, the pressure- and temperature- induced valence shift in PM phase seem like a crossover behavior 30-32 , which may imply a weak U fc . It is also found that the critical U fc of QCP mainly depends on the hybridization strength V . Due to the enhanced coherence between localized and itinerant electrons upon the increase of V , stronger U fc is required to develop a QCP of valence transition, and the critical /epsilon1 f is pushed to deeper direction, see Fig.2c.

The ground-state FOVT discussed in above section may be experimentally realizable through external pressure, e.g. in Ce and Yb systems 37-39 , however, the alteration of external factors will lead to a series of changes for model parameters in PAM which are hard to track, hence it is difficult to compare above theoretical results with experiments. Therefore, in this section, we turn to the possible valence transition caused by temperature variation in TKI, which can be more clearly interpreted theoretically. Such 'thermal valence transition' has been detected in Ce and Yb systems decades ago, in which the increase of temperature causes a FOVT from MV to Kondo state 37,39,40 , but the mechanism and transition process of thermal valence transition are still lacking of theoretical verification. In SmB 6 and TKI candidates such as golden SmS, YbB 6 and YbB 12 , the valences of Sm and Yb both increase with temperature at ambient pressure, and in some systems, the valence even reaches the saturated value 3 30,32-36 , but it seems no discontinuous valence shift was observed in PM phase. We expect that the thermal valence behavior in TKI can be dramatically alternated under different magnitudes of pressure, since the pressure applied can affect the interaction strengths in heavy-fermion systems and lead to significant change of valence in the ground state 30,31,34,36 .

At finite temperatures, the quasi-particle spectrums share the same formulas as those of the ground state, except that the mean-field parameters are now temperature-dependent, and are determined by the saddle-point equations derived from the free energy F = N λ b [ ( 2 -1) -U fc n n c f ] -2 T ∑ k , ± ln(1 + e -E ± k /T ) by ∂F/∂µ = -n t , ∂F/∂λ = 0 and ∂F/∂b = 0 , giving rise to the equation set the same as the ground state one (Eq.3), except that now the step function θ ( -E ± k ) should be replaced by Fermi distribution function f ± k = 1 (1 + / e E ± k /T ) at non-zero temperatures.

/s32

/s32

FIG. 4: (Color online) Thermal MV-Kondo and Kondo-DM FOVTs in the inset of Fig.3(b), with /epsilon1 f = -2 9376 . . (a) Energy comparison of three solutions. (b) f -electron number n f with increasing temperature, in which the dashed vertical lines denote two thermal FOVTs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00913v1-figure-f4daea6c-190b-46ec-b3c0-09e7ad7656aa.jpg)

The T /epsilon1 -f phase diagrams are given in Fig.3a and 3b, with two magnitudes of U fc = 4 1 . and 4 4 . , which are slightly smaller and greater than the QCP value U fc = 4 31 . , respectively. As temperature increases, b is reduced gradually, but only when /epsilon1 f exceeds a certain value (-0.6 and -0.525 for U fc = 4 1 . , 4 4 . , respectively), b (consequently ˜ V ) gradually vanishes at a critical temperature, results in an second-order insulator-metal transition from heavy-fermion MV insulating state to a decoupled metallic state (DM) in which n f = 1 with completely localized f electrons. Surprisingly, when /epsilon1 f is further reduced, the MV-DM transition becomes first-order, indicating a thermal FOVT, at which b and n f undergo an abrupt jump to 0 and 1, respectively, Such discontinuous thermal decoupling of heavy-fermion state and MV-DM thermal FOVT accompanied by an insulator-metal transition deserves further verification.

Remarkably, when U fc is greater than that of the QCP, we find that with /epsilon1 f slightly above that triggers the zerotemperature FOVT, there are two successive FOVTs with increasing temperature, first MV-to-Kondo FOVT, then Kondoto-DM FOVT, see Fig.4, with the transition boundaries on T /epsilon1 -f plane shown in the inset of Fig.3b. Similar to zerotemperature case, the FOVT between MV and Kondo states also holds a first-order WTI ¯ X -STI ¯ ¯ Γ X topological transition. As /epsilon1 f rises, the Kondo region gradually shrinks, then two thermal FOVTs converge, leaving only a MV-DM transition. Moreover, with /epsilon1 f lower than that of zero-temperature FOVT, only Kondo-DM FOVT occurs with increased temperature. For U fc less than the QCP value, no thermal MVKondo FOVT appears, since it is insufficient to support zerotemperature FOVT between MV and Kondo states.

In Ce systems, pressure can suppress the thermal FOVT, forces it to terminate at a CEP, leading to a thermal valence crossover behavior ? , this feature corresponds to a shift of thermal behavior from first-order to continuous one. In this sense, since pressure elevates /epsilon1 f , the /epsilon1 f -driven order change of thermal transitions may explain the appearance of CEP in Ce systems, which is marked by the red dots in Fig.3. In contrast to Ce systems, for SmB 6 and TKI candidates SmS, YbB 6 and YbB 12 , the Sm and Yb valences in PM phases seem to increase continuously with temperature 30,32-36 , showing no clear signal of thermal FOVT. The absence of thermal FOVT in these TKIs may be ascribed to their weak U fc interactions or their high f level /epsilon1 f which exceeds the CEP in Fig.3, leading to a continuous valence increase with temperature. Previous works suggested that magnetic field can shift the CEP of FOVT in Ce and Yb systems towards lower U fc and higher /epsilon1 f 40,40 , similarly, we expect magnetic field can also shift the CEP of TKI to the right side on T /epsilon1 -f plane in Fig.3, and consequently recover a thermal FOVT in TKI under ambient pressure, which requires experimental verification.

## IV. VALENCE VARIATION AROUND THE MAGNETIC TRANSITION

The high-pressure experiments of SmB 6 display a firstorder transition to magnetically ordered state at 6 ∼ 8 Gpa 25,26 , around which the Sm valence shows a rapid increase 31 , afterwards, the Sm valence in the magnetic phase increases smoothly with pressure then approaches pure trivalent above 10 Gpa. For TKI candidate golden SmS 35 , the pressureinduced magnetic transition at about 2 Gpa also exhibits firstorder character 45 , but around this transition, the valence increase of Sm in SmS is much sharper than in SmB 6 30 . In YbInCu , with enhanced pressure, FOVT reaches the mag4 netic boundary at low temperatures 38,41 , implying that the pressure-induced magnetic transition can hold a FOVT simultaneously. The active valence change at and around the magnetic transition indicates an intimate relation between these two transitions in TKI and in other heavy-fermion systems.

In order to study the valence variation around the magnetic transition in TKI, which is mostly likely an antiferromagnetic (AF) transition 48,49 , we adopt Kotliar-Ruckenstein (K-R) slave-boson mean-filed method 50-53 . This technique produces very close results of TKI to those of the conventional slave-boson approach, and has the advantage to include the magnetic order conveniently 46,47 . The main feature of KR representation for AF phase is that each f -creation and annihilation operator in c f -hybridization and f f -hopping is multiplied by a factor Z + or Z -with Z ± = √ 2(1 -n f ) 2 -n f ∓ m f , depending on the spin and sublattice of the f operator, in which m f is the staggered magnetization of f electrons. Secondly, /epsilon1 f is shifted by a parameter η . Z ± and η are analogue to b and λ in conventional slave-boson method, respectively. The K-R mean-field treatment leads to the following effective Hamiltonian

$$\mathcal { H } = N ( h m _ { f } - \eta n _ { f } - U _ { f c } n _ { c } n _ { f } ) + \sum _ { \mathbf k } \Psi _ { \mathbf k } ^ { \dagger } \mathbf H _ { \mathbf k } \Psi _ { \mathbf k }, \quad ( 4 )$$

in which the summation of k is restricted to the magnetic Brillouin zone (MBZ), h is a staggered parameter associated with m f , Ψ † k = ( d † k A ↑ , d † k A ↓ , d † k B ↑ , d † k B ↓ , f † k A ↑ , f † k A ↓ , f † k B ↑ , f † k B ↓ ) is eight-component creation operator for c and f electrons in sublattice and spin spaces, the Hamiltonian matrix is given

$$H _ { k } = \left ( \begin{array} { c c } H _ { k } ^ { d } & V _ { k } \\ V _ { k } ^ { + } & H _ { k } ^ { f } \end{array} \right ), \quad \quad \ ( 5 ) \quad \ \ \ \ 0. 6 \right ]$$

in which

$$\mathbf V _ { \mathbf k } = V \left ( \begin{array} { c c c c } 0 & 0 & Z _ { - } s _ { z } & Z _ { + } \Gamma _ { \mathbf k } \\ 0 & 0 & Z _ { - } \Gamma _ { \mathbf k } ^ { * } & - Z _ { + } s _ { z } \\ Z _ { + } s _ { z } & Z _ { - } \Gamma _ { \mathbf k } & 0 & 0 \\ Z _ { + } \Gamma _ { \mathbf k } ^ { * } & - Z _ { - } s _ { z } & 0 & 0 \end{array} \right ) \quad ( 6 ) & \quad \ \ 0. 8 \right ) \\.$$

hybridizes c and f electrons on neighboring sites.

$$\mathbf H _ { \mathbf k } ^ { d } = \left ( \begin{array} { c c } ( t _ { c } ^ { \prime } \gamma _ { \mathbf k } + U _ { f c } n _ { f } - \mu ) & u _ { \mathbf k } ^ { c } \\ u _ { \mathbf k } ^ { c } & ( t _ { c } ^ { \prime } \gamma _ { \mathbf k } + U _ { f c } n _ { f } - \mu ) \end{array} \right ) \otimes I _ { 2 }, \quad \begin{array} { c c } 1. 0 \\ \infty \\ ( 7 ) & \infty \\ \infty. 4 \end{array} \right )$$

$$\mathbf H _ { k } ^ { f } = \ \left ( \begin{array} { c c c c } e _ { k } ^ { + } & 0 & Z _ { + } Z _ { - } u _ { k } ^ { f } & 0 \\ 0 & e _ { k } ^ { - } & 0 & Z _ { + } Z _ { - } u _ { k } ^ { f } \\ Z _ { + } Z _ { - } u _ { k } ^ { f } & 0 & e _ { k } ^ { - } & 0 \\ 0 & Z _ { + } Z _ { - } u _ { k } ^ { f } & 0 & e _ { k } ^ { + } \end{array} \right ), \quad \mathbf F I G. \, 5 \colon \\ \text{(s)} \quad \text{(s)} \quad \text{from (i} \\ \text{contin} }$$

are the sublattice versions of tight-binding Hamiltonian for c and f , respectively, in which the f hopping amplitudes are renormalized by Z ± . I 2 is a two-order unit matrix, Γ k = s x -i s y , γ k = -4( c x c y + c x c z + c y c z ) , e ± k = /epsilon1 f + η + U fc n c -µ ∓ h + t ′ f Z γ 2 ± k , u c f ( ) k = t c f ( ) λ k + t ′′ c f ( ) g k with λ k = -2( c x + c y + c z ) and g k = -8 c x c y c z , where we denote c i = cos k i and s i = sin k i ( i = x, y, z ) for simplicity.

The free energy can be expressed as the expectation value of the effective Hamiltonian by

$$F = & N ( h m _ { f } - \eta n _ { f } - U _ { f c } n _ { c } n _ { f } + \mu n _ { t } ) \\ & + \sum _ { \mathbf k, n, m } ( \mathbf H _ { \mathbf k } ) _ { n m } \langle ( \Psi _ { \mathbf k } ^ { \dagger } ) _ { n } ( \Psi _ { \mathbf k } ) _ { m } \rangle, \quad \ ( 9 )$$

where subscripts n and m are row or column numbers. Use the unitary transformation matrix U k extracted from numerical diagonalization of H k , we have

$$\langle ( \Psi _ { k } ^ { \dagger } ) _ { n } ( \Psi _ { k } ) _ { m } \rangle = \sum _ { i = 1 } ^ { 8 } ( U _ { \mathbf k } ) _ { n i } ^ { * } ( U _ { \mathbf k } ) _ { m i } f _ { \mathbf k } ^ { ( i ) }, \quad ( 1 0 ) \quad \begin{array} { c } \text{in w} \end{array} \\ \mu, \, f _ { \mathbf l } \, \\ \text{tho c} \end{array}$$

where f ( ) i k = 1 (1 + / e E ( ) i k /T ) is the Fermi distribution of quasi-particles, and equals the step function θ ( -E ( ) i k ) at zero temperature. The set of equations determining the parameters n f , m f , h η , , µ can be obtained through the zero point of the derivation of free energy (Eq.9) with respect to them, by using the formulas of matrix elements ( H k ) nm in Eq.5 46 .

In order to give a global phase diagram including both PM and AF phases through K-R solution, we also solve the K-R saddle-point equations for n f , η and µ in PM phase, which

/s32

/s32

/s32

/s32

/s32

FIG. 5: (Color online) Staggered magnetization m f and f -electron occupation n f vs local level /epsilon1 f . Parameters: EHA(I), V = 1 , and from (a) to (f), U fc =2, 2.92, 3.4, 3.906, 3.94, 4.2. (a) and (b) show continuous AF transition without valence shift. For (c) and (d), the AF transition is shifted to first-order and is accompanied by a valence jump. For (e), the FOVT-associated AF transition and paramagnetic FOVT take place at different /epsilon1 f , while in (f), two transitions converge, leaving a FOVT-associated AF transition with significant valence jump.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00913v1-figure-e788ced1-7d0e-440b-bba8-68e54eca194f.jpg)

are similar to those in section II and III:

$$\text{salt} \quad & \text{are similar to those in section II and III:} \\ \text{lue} \quad & n _ { t } = \frac { 2 } { N } \sum _ { k, \pm } f _ { k } ^ { \pm }, \\ n _ { f } = \frac { 1 } { N } \sum _ { k, \pm } f _ { k } ^ { \pm } [ 1 \mp \frac { \tilde { \epsilon } _ { k } ^ { c } - \tilde { \epsilon } _ { k } ^ { f } } { \sqrt { ( \tilde { \epsilon } _ { k } ^ { c } - \tilde { \epsilon } _ { k } ^ { f } ) ^ { 2 } + 4 V ^ { 2 } Z ^ { 2 } S _ { k } ^ { 2 } } } ], \\ \text{jse} \quad & \eta = \frac { 2 Z } { N } \frac { \partial Z } { \partial n _ { f } } \sum _ { k, \pm } f _ { k } ^ { \pm } [ \epsilon _ { k } ^ { f } \pm \frac { 2 V ^ { 2 } S _ { k } ^ { 2 } - \epsilon _ { k } ^ { f } ( \tilde { \epsilon } _ { k } ^ { c } - \tilde { \epsilon } _ { k } ^ { f } ) } { \sqrt { ( \tilde { \epsilon } _ { k } ^ { c } - \tilde { \epsilon } _ { k } ^ { f } ) ^ { 2 } + 4 V ^ { 2 } Z ^ { 2 } S _ { k } ^ { 2 } } } ], \\ \text{i.e. which } \tau = \sqrt { 2 ( 1 - n _ { f } ) } \ \tilde { \epsilon } _ { f } = \frac { 2, f + \cdot \cdot \cdot \cdot \cdot \cdot \cdot T _ { f } } { \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot }$$

in which Z = √ 2(1 -n f ) 2 -n f , ˜ /epsilon1 f k = Z /epsilon1 2 f k + /epsilon1 f + η + U fc n c -µ , f ± k = 1 (1 + / e E ± k /T ) for finite temperatures and equals the step function θ ( -E ± k ) at zero temperature, where E ± k = 1 2 [˜ /epsilon1 c k +˜ /epsilon1 f k ± √ (˜ /epsilon1 c k -˜ ) /epsilon1 f k 2 +4 V 2 Z S 2 2 k ]

We first study the zero-temperature phase diagram on U fc -/epsilon1 f plane. For convenience, we choose EHA(I) in Tab.I, set various values of U fc from zero, then calculate the corresponding evolution of m f and n f vs /epsilon1 f to locate the AF transitions. Above the AF boundaries, the evolution of PM phases is computed via Eq.11 to determine the FOVT and topological boundaries. The results are shown with increasing U fc from Fig.5a to 5f. For small U fc , m f variation indicates an

/s32

/s32

/s32

/s102

FIG. 6: (Color online) U fc -/epsilon1 f magnetic phase diagram, showing the continuous AF transition (thin red line), first-order AF transition (thick red line) associating with a FOVT, and FOVT in PM phase (thick blue line). Topological transitions are denoted by dashed green lines. (b) and (c) are enlarged view of the dashed boxes in (a) and (b), respectively. Parameters: EHA(I), V =1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00913v1-figure-ee06bbfc-c45e-442e-90c1-6b050477f6bd.jpg)

second-order magnetic transition, at which n f changes continuously with a kink; when U fc exceeds 2.92, we find a firstorder magnetic transition at which n f shows an abrupt jump, indicating a FOVT-associated first-order AF transition; in the PM phase with /epsilon1 f above the first-order AF boundary, n f decreases rapidly with /epsilon1 f , and when U fc &gt; 3.906 gives rise to a FOVT at a critical /epsilon1 f above the first-order AF transition; with further increased U fc , the two transitions approach each other and then converge, leaving a large n f jump at the first-order AF transition.

The resulting ground-state phase diagram on U fc -/epsilon1 f plane is summarized in Fig.6, including the topological transitions in PM phase. In Fig.6a, it can be clearly seen that a critical point C 1 at ( U fc = 2 92 . , /epsilon1 f = -1 92 . ) separates the AF boundaries into continuous part and first-order part, which is below and above C , respectively. The FOVT in PM phase is 1 denoted by blue solid lines in Fig.6a, and the enlarged view near the QCP at ( U fc = 3 906 . , /epsilon1 f = -2 348 . ) is demonstrated in Fig.6c, in which the relative position between FOVT and topological boundaries behaves similarly to Fig.2b. Here we should point out that in K-R solution, the critical U fc and /epsilon1 f for the QCP of FOVT in PM phase is slightly greater and lower than the conventional slave-boson solution, respectively, meanwhile the FOVT and topological transition processes in PM phase remain essential the same. The interplay between first-order AF transition and PM FOVT is clearly illustrated in Fig.6b: with increasing U fc from QCP, the two boundaries approach each other and finally converge at C 2 ( U fc = 4 01 . , /epsilon1 f = -2 45 . ), leading to an single first-order AF boundary accompanied by a FOVT with significant valence jump. In addition, the coexistence of FOVT-associated AF transition and PM FOVT on /epsilon1 f axis only occurs in a narrow U fc regime from 3.906 to 4.01. Therefore, the FOVTassociated first-order AF transitions are also classified into two classes: one with weak valence shift below C 2 , and the other with strong valence jump above C 2 . Consequently, the order of AF transition, as well as the 'strength' of first-order AF transition both depend on the magnitude of U fc , in this sense, the weak first-order magnetic transitions in pressured SmB 6 and SmS may be ascribed to their relatively small d -f intra-atomic Coulomb repulsion as analysed in above sections 25,26,45 , while the remarkable valence change at the magnetic transition in YbInCu 4 may be attributed to its stronger U fc 38,41 . It is worth noting that the critical U fc = 2 92 . of C 1 which triggers a first-order AF transition is much smaller than that drives a FOVT in PM phase ( U fc = 3 906 . of QCP), explaining why Sm valence varies continuously with temperature at ambient pressure in SmB 6 and SmS, meanwhile firstorder AF transitions appear under high pressure.

Then we study the thermal phase transitions on T /epsilon1 -f plane, the results are displayed in Fig.7. For U fc smaller than the critical point C 1 , the zero-temperature AF transition with /epsilon1 f is continuous, while /epsilon1 f descends, the thermal AF transition is shifted to first-order with valence jump (see Fig.7a), implying that the critical point C 1 separating these two types of AF transitions is temperature dependent and is pushed by increasing temperature towards lower /epsilon1 f and weaker U fc . By comparison, for U fc stronger than zero-temperature C 1 , thermal AF transition remains first-order, because U fc keeps greater than C 1 when temperature rises. For temperature above the AF transitions, Z factor reduces with temperature, leading to the increase of n f then finally a decoupling of heavy-fermion state to DM state with n f = 1 at a critical temperature, giving rise to HF-DM transition. By contrast to AF transitions, the order of this HF-DM transition is only weakly temperaturedependent and mainly depends on U fc : when U fc is less or larger than zero-temperature C 1 in Fig.6a, HF-DM transition is continuous or first-order, respectively, and the latter coex-

/s84

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00913v1-figure-708e6aed-001f-441a-9a95-4d173133f179.jpg)

/s32

FIG. 7: (Color online) AF transitions (blue lines) and transitions between heavy-fermion state and DM state (green lines). Thin lines denote continuous transitions, while thick lines denote first-order transitions with valence shift. Parameter: EHA(I), V = 1 . U fc = 2 and 3 4 . in (a) and (b), respectively.

ists a FOVT. As /epsilon1 f descends, the AF boundary and HF-DM boundary gradually converge, leaving a first-order AF-DM transition with abrupt valence shift. The thermal first-order AF transition coexisting a FOVT deserves more rigorous theoretical and experimental examinations.

## V. CONCLUSION AND DISCUSSION

To summarize, we have proposed a series of exotic valence transition processes in TKI, driven by the effect of on-site Coulomb repulsion U fc between conduction and local electrons in the PAM. For the ground state, we located a QCP on the U fc -/epsilon1 f plane, above which a FOVT boundary separates the Kondo state with MV state, and holds an abrupt valence jump upon ascending /epsilon1 f . The QCP is very close to a topological transition, and the FOVT boundary gradually approaches then converges with the topological boundary, leading to a discontinuous topological transition. Near the parameter regime of the zero-temperature FOVT, increase in temperature can lead to a series of thermal FOVTs among the Kondo state, MV state and DM state. In addition, we found a CEP on T /epsilon1 -f plane which terminates the thermal FOVT

∗ Electronic address: lihuan@glut.edu.cn

† Electronic address: liu yu@iapcm.ac.cn

1 Maxim Dzero, Kai Sun, Victor Galitski, and Piers Coleman, Phys. Rev. Lett. 104 , 106408 (2010).

2 Maxim Dzero, Kai Sun, Piers Coleman, and Victor Galitski, Phys. Rev. B 85 , 045130 (2012).

3 Minh-Tien Tran, Tetsuya Takimoto, and Ki-Seok Kim, Phys. Rev. B 85 , 125128 (2012).

4 Victor Alexandrov, Maxim Dzero, and Piers Coleman, Phys. Rev.

from heavy-fermion state to DM state, and shift the thermal decoupling of heavy-fermion state into continuous one under higher /epsilon1 f above the CEP. Further descent of /epsilon1 f can generate an AF order in TKI, and above a critical U fc which is considerably smaller than that of the QCP which triggers a PM FOVT, the /epsilon1 f -driven continuous AF transition is shifted to first-order, accompanied by a FOVT. Further increasing of U fc can push the first-order AF transition to approach then to join the PM FOVT boundary, giving rise to a first-order AF transition with large valence shift. The thermal AF transitions can also be classified into continuous one or first-order one, depending mainly on the magnitude of U fc . The distinct valence-variation processes with temperature and pressure observed in TKI candidates SmB 6 , golden SmS and other heavy-fermion compounds, as well as their valence behaviors near the magnetic transitions can be qualitatively understood by our results in terms of their different magnitudes of U fc . This work has used an simplified PAM with single conduction and f band to describe TKI, actually in real cubic TKI systems such as SmB , dominant hybridization channel in6 volves spin-degenerated 5 d x 2 -y 2 , 5 d 3 z 2 -r 2 orbits, and Γ (1) 8 , Γ (2) 8 4 f quartet, resulting in much complicated tight-binding and hybridization terms 4,7,54 . Although the multiple parameters can all influence the detailed valence-variation process in TKI, we believe that the intra-atomic Coulomb interaction U fc between d and f orbits still provides a crucial factor controlling the valence transition, and the magnitudes of U fc in TKI candidates may be examined through first-principle simulations 22 .

## Acknowledgments

H. Li is supported by NSFC (No. 11764010) and Guangxi Natural Science Foundation (No. 2017GXNSFAA198169). Y. Liu and H. F. Song thank the Science Challenge Project (Grant No. TZ2018002). Y. Liu is also supported by Open Research Fund Program of the State Key Laboratory of LowDimensional Quantum Physics (No. KF201702) and the SPC-Lab Research Fund (No. XKFZ201605). Z. Y. Wang is supported by NSFC (No. 11564008). X. J. Zheng is supported by NSFC (No. 11704084) and Guangxi NSF (No. 2017GXNSFBA198115). X. Ming is supported by NSFC (No. 11864008) and Guangxi NSF (No. 2018GXNSFAA138185).

Lett. 111 , 226403 (2013).

5 Markus Legner, Andreas R¨ uegg, and Manfred Sigrist, Phys. Rev. B 89 , 085110 (2014).

6 Jan Werner and Fakher F. Assaad, Phys. Rev. B 89 , 245119 (2014).

7 Bitan Roy, Jay D. Sau, Maxim Dzero, and Victor Galitski ,Phys. Rev. B 90 , 155314 (2014).

8 Junwon Kim, Kyoo Kim, Chang-Jong Kang, Sooran Kim, Hong Chul Choi, J.-S. Kang, J. D. Denlinger, and B. I. Min, Phys. Rev.

- B 90 , 075131 (2014).
- 9 Pier Paolo Baruselli and Matthias Vojta, Phys. Rev. B 90 , 201106(R) (2014).
- 10 G. Li, Z. Xiang, F. Yu, T. Asaba, B. Lawson, P. Cai, C. Tinsman, A. Berkley, S. Wolgast, Y. S. Eo, Dae-Jeong Kim, C. Kurdak, J. W. Allen, K. Sun, X. H. Chen, Y. Y. Wang, Z. Fisk, Lu Li, Science 346 , 1208 (2014).
- 11 N. Xu, P.K. Biswas, J.H. Dil1, R.S. Dhaka, G. Landolt, S. Muff, C.E. Matt, X. Shi, N.C. Plumb, M. Radovic, E. Pomjakushina, K. Conder, A. Amato, S.V. Borisenko, R. Yu, H.-M. Weng, Z. Fang, X. Dai, J. Mesot1, H. Ding, and M. Shi, Nat. Commun. 5 , 4566 (2014).
- 12 Victor Alexandrov, Piers Coleman, and Onur Erten, Phys. Rev. Lett. 114 , 177202 (2015).
- 13 Zhen Zheng, Xu-Bo Zou and Guang-Can Guo, New J. Phys. 20 , 023039 (2018).
- 14 Harris Pirie, Yu Liu, A. Soumyanarayanan, Pengcheng Chen, Yang He, M.M. Yee, P.F.S. Rosa, J.D. Thompson, Dae-Jeong Kim, Z. Fisk, Xiangfeng Wang, J. Paglione, Dirk K. Morr, M. H. Hamidian, and Jennifer E. Hoffman, arXiv:1810.13419v1.
- 15 Jianjun Ying, Lingyun Tang, Fei Chen, Xianhui Chen, and Viktor V. Struzhkin, Phys. Rev. B 97 , 121101(R) (2018).
- 16 J.J van den Broekey, S.N. Kempkesy, A. Quelley, X.F. Wangx, J. Paglionex, and C. Morais Smith, arXiv: 1803.03553v1.
- 17 Yoshiyuki Ohtsubo, Yuki Yamashita, Kenta Hagiwara, Shinichiro Ideta, Kiyohisa Tanaka, Ryu Yukawa, Koji Horiba, Hiroshi Kumigashira, Koji Miyamoto, Taichi Okuda, Wataru Hirano, Fumitoshi Iga, and Shin-ichi Kimura1, arXiv:1803.09433v3.
- 18 W. T. Fuhrman, J. C. Leiner, J. W. Freeland, M. van Veenendaal, S. M. Koohpayeh, W. Adam Phelan, T. M. McQueen, and C. Broholm, Phys. Rev. B 97 , 121101(R) (2018).
- 19 Maxim Dzero and Victor Galitski, J. Exp. Theor. Phys. 117 , 499 (2013).
- 20 Maxim Dzero, Jing Xia, Victor Galitski, and Piers Coleman, Annu. Rev. Condens. Matter Phys. 7 249 (2016)
- 21 Nan Xu, Hong Ding, and Ming Shi, J. Phys.: Condens. Matter 28 , 363001 (2016).
- 22 Feng Lu, JianZhou Zhao, Hongming Weng, Zhong Fang, and Xi Dai, Phys. Rev. Lett. 110 , 096401 (2013).
- 23 Rui Yu, Hongming Weng, XiaoHu, Zhong Fang, and Xi Dai, New J. Phys. 17 , 023012 (2015).
- 24 Markus Legner, Andreas R¨ uegg, and Manfred Sigrist, Phys. Rev. Lett. 115 , 156405 (2015).
- 25 A. Barla, J. Derr, J. P. Sanchez, B. Salce, G. Lapertot, B. P. Doyle, R. R¨ uffer, R. Lengsdorf, M. M. Abd-Elmeguid, and J. Flouquet, Phys. Rev. Lett. 94 , 166401 (2005).
- 26 J. Derr, G. Knebel, G. Lapertot, B. Salce, M-A M´ easson and J. Flouquet, J. Phys.: Condens. Matter 18 , 2089 (2006).
- 27 J. Derr, G. Knebel, D. Braithwaite, B. Salce, J. Flouquet, K. Flachbart, S. Gab´ ani, and N. Shitsevalova, Phys. Rev. B 77 , 193107 (2008).
- 28 Parisiades Paraskevas, Bremholm Martin, and Mezouar Mohamed, Europhys. Letts. 110 , 66002 (2015).
- 29 N. Emi, N. Kawamura, M. Mizumaki, T. Koyama, N. Ishimatsu, G. Prist´ a˘, s T. Kagayama, K. Shimizu, Y. Osanai, F. Iga, and T. Mito, Phys. Rev. B 97 , 161116(R) (2018).
- 30 Nicholas P. Butch, Johnpierre Paglione, Paul Chow, Yuming Xiao, Chris A. Marianetti, Corwin H. Booth, and Jason R. Jeffries, Phys. Rev. Lett. 116 , 156401 (2016).
- 31 Yazhou Zhou, Qi Wu, Priscila F. S. Rosa, Rong Yu, Jing Guo, Wei Yi, Shan Zhang, Zhe Wang, Honghong Wang, Shu Cai, Ke Yang, Aiguo Li, Zheng Jiang, Suo Zhang, Xiangjun Wei, Yuying Huang, Peijie Sun, Yi-feng Yang, Zachary Fisk, Qimao Si, Zhongxian Zhao, Liling Sun, Science Bulletin 62 , 1439 (2017).
- 32 Masaichiro Mizumaki, Satoshi Tsutsui, and Fumitoshi Iga, J. Phys.: Conf. Ser. 176 , 012034 (2009).
- 33 Qi Wu and Liling Sun, Rep. Prog. Phys. 80 112501 (2017).
- 34 E. Annese, A. Barla, C. Dallera, G. Lapertot, J-P. Sanchez, and G. Vank´ o, Phys. Rev. B 73 , 140409(R) (2006).
- 35 Zhi Li, Jin Li, Peter Blaha, and Nicholas Kioussis, Phys. Rev. B 89 , 121117 (2014).
- 36 Kenta Hagiwara, Yusuke Takeno, Yoshiyuki Ohtsubo, Ryu Yukawa, Masaki Kobayashi, Koji Horiba, Hiroshi Kumigashira, Julien Rault, Patrick Le F` evre, Franc ¸ois Bertran, Amina TalebIbrahimi, Fumitoshi Iga, Shin-ichi Kimura, J. Phys.: Conf. Series 807 , 012003 (2017).
- 37 I. Felner and I. Nowik, Phys. Rev. B 33 , 617 (1985).
- 38 T. Mito, T. Koyama, M. Shimoide, S. Wada, T. Muramatsu, T. C. Kobayashi, and J. L. Sarrao, Phys. Rev. B 67 , 224409 (2003).
- 39 T. Park, V. A. Sidorov, J. L. Sarrao, and J. D. Thompson, Phys. Rev.Lett. 96 ,046405 (2006).
- 40 Shinji Watanabe, Atsushi Tsuruta, Kazumasa Miyake, and Jacques Flouquet, J. Phys. Soc. Jpn 78 , 104706 (2009).
- 41 C. D. Immer, J. L. Sarrao, Z. Fisk, A. Lacerda, C. Mielke, and J. D. Thompson, Phys. Rev. B 56 , 71 (1997).
- 42 A. V. Goltsev and G. Bruls, Phys. Rev. B 63 , 155109 (2001).
- 43 Shinji Watanabe, Masatoshi Imada, and Kazumasa Miyake, J. Phys. Soc. Jpn 75 , 043710 (2006).
- 44 Shinji Watanabe, Atsushi Tsuruta, Kazumasa Miyake, and Jacques Flouquet, Phys. Rev. Lett. 100 , 236401 (2008).
- 45 A. Barla, J. P. Sanchez, Y. Haga, G. Lapertot, B. P. Doyle, O. Leupold, R. R¨ uffer, M.M. Abd-Elmeguid, R. Lengsdorf, and J. Flouquet, Phys. Rev. Lett. 92 , 066401 (2004).
- 46 Huan Li, Yin Zhong, Yu Liu, Hong-Gang Luo and Hai-Feng Song, Journal of Physics: Condensed Matter 30 , 435601 (2018).
- 47 Huan Li,Zhi-Yong Wang, Xiao-Jun Zheng, Yu Liu and Yin Zhong,Chin. Phys. Lett. 35 , 127501 (2018).
- 48 Robert Peters, Tsuneya Yoshida, and Norio Kawakami, Phys. Rev. B 98 , 075104 (2018).
- 49 Kai-Wei Chang and Peng-Jen Chen, Phys. Rev. B 97 , 195145 (2018).
- 50 Gabriel Kotliar and Andrei E. Ruckenstein, Phys. Rev. Lett. 57 , 1362 (1986).
- 51 Min-Fong Yang, Shih-Jye Sun, and Tzay-Ming Hong, Phys. Rev. B 48 , 16123 (1993).
- 52 Shih-Jye Sun, Min-Fong Yang, and Tzay-Ming Hongt, Phys. Rev. B 48 , 16127 (1993).
- 53 Shih-Jye Sun, Tzay-Ming Hong, Min Fong Yang, Physica B 216 , 111 (1995).
- 54 Tetsuya Takimoto, J. Phys. Soc. Jpn 80 , 123710 (2011).
Jia-Tao Zhuang Xiao-Jun Zheng Zhi-Yong Wang Xing Ming Huan Li Yu Liu Hai-Feng Song
Valence transition in topological Kondo insulator
We investigate the valence transition in three-dimensional topological Kondo
insulator through slave-boson analysis of periodic Anderson model. By including
the effect of intra-atomic Coulomb correlation $U_{fc}$ between conduction and
local electrons, we find a first-order valence transition from Kondo region to
mixed valence upon ascending of local level above a critical $U_{fc}$, and this
valence transition usually occurs very close to or simultaneously with a
topological transition. Near the parameter region of zero-temperature valence
transition, rise of temperature can generate a thermal valence transition from
mixed valence to Kondo region, accompanied by a first-order topological
transition. Remarkably, above a critical $U_{fc}$ which is considerable smaller
than that generating paramagnetic valence transition, the original continuous
antiferromagnetic transition is shifted to first order one, at which a
discontinuous valence shift takes place. Upon increased $U_{fc}$, the
paramagnetic valence transition approaches then converges with the first-order
antiferromagnetic transition, leaving an significant valence shift on the
magnetic boundary. The continuous antiferromagnetic transition, first-order
antiferromagnetic transition, paramagnetic valence transition and topological
transitions are all summarized in a global phase diagram. Our proposed exotic
transition processes can help to understand the thermal valence variation as
well as the valence shift around the pressure-induced magnetic transition in
topological Kondo insulator candidates and in other heavy-fermion systems.