• The emergent dynamical processes associated with the magnetic excitations in heavy fermion SmB$_6$ are investigated. By imposing multi-orbital interactions on a first principles tight binding model, we find the interplay between magnetic and orbital fluctuations in the localized $f$-levels is highly sensitive to the Hubbard interaction $U$ and the Hund's coupling $J$. The magnetic phase diagram constructed at zero temperature reveals novel quantum critical features with the existence of several competing magnetic phases. Within the random phase approximation, we perform a comprehensive study of the spin-spin correlation function, and our results agree with the neutron scattering experiments. Furthermore, analysis of the spectral weight of the magnetic excitations shows that some collective spin excitations are accompanied by orbital fluctuations, indicating an interesting entanglement between the spin and the orbital degrees of freedom in SmB$_6$. Our work underscores the importance of the orbital degree of freedom in the physics of heavy fermion systems.
  • The cluster slave-spin method is introduced to systematically investigate the solutions of the Hubbard model including the symmetry-broken phases. In this method, the electron operator is factorized into a fermioninc spinon describing the physical spin and a slave-spin describing the charge fluctuations. Following the $U(1)$ formalism derived by Yu and Si [Phys. Rev. B 86, 085104 (2012)], it is shown that the self-consistent equations to explore various symmetry-broken density wave states can be constructed in general with a cluster of multiple slave-spin sites. We employ this method to study the antiferromagnetic (AFM) state in the single band Hubbard model with the two and four-site clusters of slave spins. While the Hubbard gap, the charge gap due to the doubly-occupied states, scales with the Hubbard interaction $U$ as expected, the AFM gap $\Delta$, the gap in the spinon dispersion in the AFM state, exhibits a crossover from the weak to strong-coupling behaviors as $U$ increases. Our cluster slave-spin method reproduces not only the traditional mean-field behavior of $\Delta \sim U$ in the weak-coupling limit, but also the behavior of $\Delta \sim t^2/U$ predicted by the superexchange mechanism in the strong coupling limit. In addition, the holon-doublon correlator as functions of $U$ and doping $x$ is also computed, which exhibits a strong tendency toward the holon-doublon binding in the strong coupling regime. We further show that the quasiparticle weight obtained by the cluster slave-spin method is in a good agreement with the generalized Gutzwiller approximation in both AFM and paramagnetic states, and the results can be improved beyond the generalized Gutzwiller approximation as the cluster is enlarged from a single site to 4 sites. Our results demonstrate that the cluster slave-spin method can be a powerful tool to systematically investigate the strongly correlated system.
  • We investigate the charge excitations of a Weyl semimetal in the axionic charge density wave (axionic CDW) state. While it has been shown that the topological response (anomalous Hall conductivity) is protected against the CDW state, we find that the long wavelength plasmon excitation is radically influenced by the dynamics of the CDW order parameter. In the normal state, we show that an undamped collective mode should exist at q~Q_(CDW) if there is an attractive interaction favoring the formation of the CDW state. The undamped nature of this collective mode is attributed to a gap-like feature in the particle-hole continuum at q~Q_(CDW) due to the chirality of the Weyl nodes, which is not seen in other materials with CDW instability. In the CDW state, the long wavelength plasmon excitations become more dispersive due to the additional interband scattering not allowed in the normal state. Moreover, because the translational symmetry is spontaneously broken, Umklapp scattering, the process conserving the total momentum only up to nQ_(CDW) with n an integer and Q_(CDW) the ordering wave vector, emerges in the CDW state. We find that the plasmon excitation couples to the phonon mode of the CDW order via the Umklapp scattering, leading to two branches of resonant collective modes observable in the density-density correlation function at q~0 and q~Q_(CDW). Based on our analysis, we propose that measuring these resonant plasmon-axion excitations around q~0 and q~Q_(CDW) by the momentum-resolved electron energy loss spectroscopy (M-EELS) could serve as a reliable way to detect the axionic CDW state in Weyl semimetals.
  • We review recent progress in point contact spectroscopy (PCS) to extract spectroscopic information out of correlated electron materials, with the emphasis on non-superconducting states. PCS has been used to detect bosonic excitations in normal metals, where signatures (e.g. phonons) are usually less than 1$\%$ of the measured conductance. In the superconducting state, point contact Andreev reflection (PCAR) has been widely used to study properties of the superconducting gap in various superconductors. In the last decade, there have been more and more experimental results suggesting that the point contact conductance could reveal new features associated with the unusual single electron dynamics in non-superconducting states, shedding a new light on exploring the nature of the competing phases in correlated materials. We will summarize the theories for point contact spectroscopy developed from different approaches and highlight these conceptual differences distinguishing point contact spectroscopy from tunneling-based probes. Moreover, we will show how the Schwinger-Kadanoff-Baym-Keldysh (SKBK) formalism together with the appropriate modeling of the nano-scale point contacts randomly distributed across the junction leads to the conclusion that the point contact conductance is proportional to the {\it effective density of states}, a physical quantity that can be computed if the electron self energy is known. The experimental data on iron based superconductors and heavy fermion compounds will be analyzed in this framework. These recent developments have extended the applicability of point contact spectroscopy to correlated materials, which will help us achieve a deeper understanding of the single electron dynamics in strongly correlated systems.
  • We investigate the superconductivity in the $t-J-U$ model within a slave-spin formalism. We show that the BCS mean-field theory implemented with the slave spin formalism naturally predicts two distinct gaps which are the pairing gap of the spinons $\Delta_f$ and the Cooper pairing gap of the electrons $\Delta_{SC} = Z\Delta_f$, where $Z$ is the quasiparticle weight. Furthermore, we find that the nature of the superconducting state depends crucially on the interaction providing the pairing mechanism. For spin interactions, the bandwidth of the spinon hopping term is renormalized by $Z$ but the its pairing term is not. As a result, if $U$ exceeds the critical value for the Mott insulating state at half-filling, $Z$ develops a strong doping dependence, leading to a doping-driven crossover from strong to weak pairing states. In the strong pairing state, while $\Delta_f$ is enhanced as $x\to 0$, $\Delta_{SC}\sim x$ due to the renormalization of $Z$. In the weak pairing state, $Z$ does not change with $x$ significantly. Therefore, $\Delta_{SC}$ is mainly controlled by $\Delta_f$, and both of them go to zero at larger doping. The crossover from strong to weak pairing states is well captured by the slave spin formalism within reasonable range of parameters just at the mean-field level, indicating the slave spin formalism is a powerful tool to study correlated materials. For charge interactions, we find that the bandwidth of the spinon hopping term and its pairing term are renormalized by $Z$ and $Z^2$ respectively. Consequently, both $\Delta_f$ and $\Delta_{SC}$ are suppressed at small $x$ in large U, and no crossover will occur. The implication of our results for the superconducting states in correlated materials will be discussed.
  • We investigate the evolution of magnetic excitations in Sr$_3$Ru$_2$O$_7$ using a three band tight binding model that takes into account the influence of Mn and Ti dopant ions. The effect of dopant ions on the Sr$_3$Ru$_2$O$_7$ band structure has been included by taking into account the dopant induced suppression of the oxygen octahedral rotation in the tight binding band structure and changes in electron occupation. We find that the low energy spin fluctuations are dominated by three wave vectors around Q=$((0,0),(\pi/2,\pi/2))$, and $(\pi,0)$ which compete with each other. As the octahedral rotation is suppressed with increasing doping, the three wave vectors evolve differently. In particular, the undoped compound has dominant wavevectors at Q=$((0,0),(\pi/2,\pi/2))$, but doping Sr$_3$Ru$_2$O$_7$ leads to a significant enhancement in the spin susceptibility at the Q=$(\pi,0)$ wavevector bringing the system closer to a magnetic instability. All the features calculated from our model are in agreement with neutron scattering experiments. We have also studied the effect of a c-axis Zeeman field on the low energy spin fluctuations. We find that an increasing magnetic field suppresses the AFM fluctuations and leads to stronger competition between the AFM and FM spin fluctuations. The magnetic field dependence observed in our calculations therefore supports the scenario that the observed nematic phase in the metamagnetic region in Sr$_3$Ru$_2$O$_7$ is intimately related to the presence of a competing ferromagnetic instability.
  • We present evidence of strain-induced modulation of electron correlation effects and increased orbital anisotropy in the rutile phase of epitaxial VO$_2$/TiO$_2$ films from hard x-ray photoelectron spectroscopy and soft V L-edge x-ray absorption spectroscopy, respectively. By using the U(1) slave spin formalism, we further argue that the observed anisotropic correlation effects can be understood by a model of orbital selective Mott transition at a filling that is non-integer, but close to the half-filling. Because the overlaps of wave functions between $d$ orbitals are modified by the strain, orbitally-dependent renormalizations of the bandwidths and the crystal fields occur with the application of strain. These renormalizations generally result in different occupation numbers in different orbitals. We find that if the system has a non-integer filling number near the half-filling such as for VO$_2$, certain orbitals could reach an occupation number closer to half-filling under the strain, resulting in a strong reduction in the quasiparticle weight $Z_{\alpha}$ of that orbital. Moreover, an orbital selective Mott transition, defined as the case with $Z_{\alpha} = 0$ in some, but not all orbitals, could be accessed by epitaxial strain-engineering of correlated electron systems.
  • The quasiparticle interferences (QPIs) of a Mott insulator are investigated using the $T$-matrix formalism implemented with the dynamical mean-field theory ($T$-DMFT). In Mott insulating state, because DMFT predicts a singularity in the real part of the electron self energy $\Re \Sigma(\omega) \sim \eta/\omega$ at low frequency, where $\eta$ can be considered as the 'order parameter' for Mott insulating state, QPIs are completely washed out at the small bias voltage. However, a further analysis shows that $\Re \Sigma(\omega)$ in fact serves as an energy-dependent chemical potential shift. As a result, the effective bias voltage seen by the system is $eV' = eV - \Re \Sigma(eV)$. Due to the singular behavior of $\Re \Sigma(\omega)$, a critical bias voltage $eV_c = \sqrt{\eta}$ satisfying $eV'=0$ exists if and only if in Mott insulating state. Consequently, the same QPI patterns produced by the non-interacting Fermi surfaces appears at this critical bias voltage $eV_c$ in Mott insulating state. We propose that this reentry of non-interacting QPI patterns at $eV_c$ could serve as an experimental signature of Mott insulating state, and the 'order parameter' can be experimentally measured as $\eta \sim (eV_c)^2$.
  • Motivated by the mid-infrared scenario for high-temperature superconductivity proposed by Leggett, the effects of Umklapp processes on the density-density correlation function in the presence of long-range Coulomb interaction have been investigated on a microscopic model. We show that because Umklapp processes enable scatterings that conserve total momentum only up to $n\vec{K}$, where $n$ is an integer and $\vec{K}$ is the reciprocal wave vector, significant amounts of spectral weight in the plasmonic excitations at long wavelength are transferred into lower frequency around mid-infrared regime. We further find that regardless of the gap symmetry, superconductivity generally suppresses the Umklapp scatterings due to the formation of the electron pairing. This suppression is unique for the superconductivity due to the interplay between electron pairing and the odd parity of the matrix elements associated with Umklapp channels, which usually does not occur in other known competing orders. Specific predictions for the experimental signatures in optical conductivity and electron energy loss spectroscopy (EELS) will be discussed.
  • We develop a microscopic theory for the point-contact conductance between a metalic electrode and a strongly correlated material using the non-equilibrium Schwinger-Kadanoff-Baym-Keldysh formalism. We explicitly show that in the classical limit, contact size shorter than the scattering length of the system, the microscopic model can be reduced to an effective model with transfer matrix elements that conserve in-plane momentum. We find that the conductance $dI/dV$ is proportional to the {\it effective density of states}, that is, the integrated single-particle spectral function $A(\omega=eV)$ over the whole Brillouin zone. From this conclusion, we are able to establish the conditions under which a non-Fermi liquid metal exhibits a zero-bias peak in the conductance. This finding is discussed in the context of recent point-contact spectroscopy on the iron pnictides and chalcogenides which has exhibited a zero-bias conductance peak.
  • Using a degenerate double-exchange model, we investigate the spin excitation spectra of iron pnictides. The model consists of local spin moments on each Fe site, as well as itinerant electrons from the degenerate $d_{xz}$ and $d_{yz}$ orbitals. The local moments interact with each other through antiferromagnetic $J_{1}$-$J_{2}$ Heisenberg interactions, and they couple to the itinerant electrons through a ferromagnetic Hund coupling. We employ the fermionic spinon representation for the local moments and perform a generalized random-phase approximation calculation on both spinons and itinerant electrons. We find that in the $\left(\pi,0\right)$ magnetically-ordered state, the spin-wave excitation at $\left(\pi,\pi\right)$ is pushed to a higher energy due to the presence of itinerant electrons, which is consistent with a previous study using the Holstein-Primakoff transformation. In the paramagnetic state, the particle-hole continuum keeps the collective spin excitation near $\left(\pi,\pi\right)$ at a higher energy even without any $C_{4}$ symmetry breaking. The implications for recent high temperature neutron scattering measurements will be discussed.
  • We show that the fluctuations associated with ferro orbital order in the $d_{xz}$ and $d_{yz}$ orbitals can develop a sharp resonance mode in the superconducting state with a nodeless gap on the Fermi surface. This orbital resonance mode appears below the particle-hole continuum and is analogous to the magnetic resonance mode found in various unconventional superconductors. If the pairing symmetry is $s_{\pm}$, a dynamical coupling between the orbital ordering and the d-wave subdominant pairing channels is present by symmetry. Therefore the nature of the resonance mode depends on the relative strengths of the fluctuations in these two channels, which could vary significantly for different families of the iron based superconductors. The application of our theory to a recent observation of a new $\delta$-function-like peak in the B$_{1g}$ Raman spectrum of Ba$_{0.6}$K$_{0.4}$Fe$_2$As$_2$ is discussed.
  • One central issue under intense debate in the study of the iron based superconductors is the origin of the structural phase transition that changes the crystal lattice symmetry from tetragonal to orthorhombic. This structural phase transition, occurring universally in almost every family of the iron-based superconductors, breaks the lattice $C_4$ rotational symmetry and results in an anisotropy in a number of physical properties. Due to the unique topology of the Fermi surface, both orbital- and spin-based scenarios have been proposed as the driving force. In this review, we focus on theories from the orbital-based scenario and discuss several related experiments. It is pointed out that although both scenarios lead to the same macroscopic phases and are not distinguishable in bulk measurements of the thermodynamic properties, the elementary excitations could be fundamentally different, and provide us with the possibility to resolve this long-standing debate between orbital- and spin-based theories.
  • The critical behavior of a two-orbital model with degenerate $d_{xz}$ and $d_{yz}$ orbitals is investigated by multidimensional bosonization. We find that the corresponding bosonic theory has an overdamped collective mode with dynamical exponent $z=3$, which appears to be a general feature of a two-orbital model and becomes the dominant fluctuation in the vicinity of the orbital-ordering quantum critical point. Since the very existence of this $z=3$ overdamped collective mode induces non-Fermi liquid behavior near the quantum critical point, we conclude that a two-orbital model generally has a sizable area in the phase diagram showing non-Fermi liquid behavior. Furthermore, we show that the bosonic theory resembles the continuous model near the d-wave Pomeranchuk instability, suggesting that orbital order in a two-orbital model is identical to nematic order in a continuous model. Our results can be applied to systems with degenerate $d_{xz}$ and $d_{yz}$ orbitals such as iron-based superconductors and bilayer strontium ruthenates Sr$_3$Ru$_2$O$_7$.
  • We study the influence of quantum fluctuations on the electron self energy in the normal state of iron-pnictide superconductors using a five orbital tight binding model with generalized Hubbard on-site interactions. Within a one-loop treatment, we find that an overdamped collective mode develops at low frequency in channels associated with quasi-1D $d_{xz}$ and $d_{yz}$ bands. When the critical point for the $C_4$ symmetry broken phase (structural phase transition) is approached, the overdamped collective modes soften, and acquire increased spectral weight, resulting in non-Fermi liquid behavior at the Fermi surface characterized by the frequency dependence of the imaginary part of electron self energy of the form $\omega^\lambda$, $0<\lambda<1$. We argue that this non-Fermi liquid behavior is responsible for the recently observed zero-bias enhancement in the tunneling signal in quantum point contact spectroscopy. A key experimental test of this proposal is the absence of the non-Fermi liquid behaviour in the hole-doped materials. Our result suggests that quantum criticality plays an important role in understanding the normal state properties of iron-pnictide superconductors.
  • We show here that orbital degrees of freedom produce a distinct signature in the magnetic excitation spectrum of iron-based superconductors above the magnetic ordering temperature. Because $d_{xz}$ and $d_{yz}$ orbitals are strongly connected with the Fermi surface topology, the nature of magnetic excitations can be modified significantly due to the presence of either static or fluctuating orbital correlations. Within a five-orbital itinerant model, we show that static orbital order generally leads to an enhancement of commensurate magnetic excitations even when the original Fermi surface lacks nesting at $(\pi,0)$ or $(0,\pi)$. When long-range orbital order is absent, Gaussian fluctuations beyond the standard random-phase approximation (RPA) capture the effects of fluctuating orbital correlations on the magnetic excitations. We find that commensurate magnetic excitations can also be enhanced if the orbital correlations are strong. Our results offer a natural explanation for the incommensurate-to-commensurate transformation observed in a recent neutron scattering measurement (Z. Xu, et. al., arXiv:1201.4404), and we propose that this unusual transformation is an important signature to distinguish orbital from spin physics in the normal state of the pnictides. Implications for the magnetic and superconducting states are discussed.
  • We suggest a mechanism which promotes the existence of a phase soliton -- topological defect formed in the relative phase of superconducting gaps of a two-band superconductor with s+- type of pairing. This mechanism exploits the proximity effect with a conventional s-wave superconductor which favors the alignment of the phases of the two-band superconductor which, in the case of s+- pairing, are pi-shifted in the absence of proximity. In the case of a strong proximity such effect can be used to reduce soliton's energy below the energy of a soliton-free state thus making the soliton thermodynamically stable. Based on this observation we consider an experimental setup, applicable both for stable and metastable solitons, which can be used to distinguish between s+- and s++ types of pairing in the iron-based multiband superconductors.
  • We investigate the effects of the $\sqrt{5}\times\sqrt{5}$ Fe vacancy ordering on the orbital and magnetic order in (K,Tl,Cs)$_y$Fe$_{2-x}$Se$_2$ using a three-orbital ($t_{2g}$) tight-binding Hamiltonian with generalized Hubbard interactions. We find that vacancy order enhances electron correlations, resulting in the onset of a block antiferromagnetic phase with large moments at smaller interaction strengths. In addition, vacancy ordering modulates the kinetic energy differently for the three $t_{2g}$ orbitals. This results in a breaking of the degeneracy between the $d_{xz}$ and $d_{yz}$ orbitals on each Fe site, and the onset of orbital order. Consequently, we obtain a novel inverse relation between orbital polarization and the magnetic moment. We predict that a transition from high-spin to low-spin states accompanied by a crossover from orbitally-disordered to orbitally-ordered states will be driven by doping the parent compound with electrons, which can be verified by neutron scattering and soft X-ray measurements.
  • We develop here a general formalism for multi-orbital Mott systems which can be used to understand dynamical and static spectral weight transfer. We find that the spectral weight transferred from the high energy scales is greatly increased as a result of the multi-orbital structure. As a consequence certainly dynamically generated symmetries obtain at lower values of doping than in the single-band Hubbard model. For example, in the atomic limit, the particle-hole symmetric condition in the lower band shifts from the one-band result of $x=1/3$ to $x=1/(2n_o+1)$, where $n_o$ is the number of orbitals with an unpaired spin. Transport properties computed from effective low-energy theories which forbid double occupancy of bare electrons, such as the multi-orbital t-J generalization, should all be be sensitive to this particle-hole symmetric condition. Away from the atomic limit, the dynamical contributions increase the transferred spectral weight. Consequently, any phenomena which are sensitive to an emergent particle-hole symmetry should obtain at $x<(1/(2n_o+1)$.
  • In a recent Letter (Phys. Rev. Lett. 106, 036401 (2011)), Lau, et al. investigated the single hole problem in an effective model containing Cu $d$ and O $2p$ orbitals by exact diagonalization on clusters composed of 32 CuO$_2$ unit cells. With full quantum fluctuations due to the antiferromagnetic (AFM) backgroud, they found that spin-polaron solutions with spin 3/2 become lowest energy solutions in certain region of Brillouin zone (BZ), which they claimed no similiar solutions were obtained in other models or approximations before. In this Comment, we would like to point out that such high-spin polaron solutions have been found in our previous Letter (Phys. Rev. Lett. 91, 057001, (2003)) and many of features described by Lau et al. have in fact been seen in our Letter.
  • Geometric frustration in quantum magnetism refers to that magnetic interactions on different bonds cannot be simultaneously minimized. The usual Cooper pairing systems favor the uniform distribution of the pairing phase among lattice sites without frustration. In contrast, we propose "frustrated Cooper pairing" in non-bipartite lattices which leads to frustrated supersolid states with non-uniform distributions of the Cooper pair phase and density. This exotic pairing state naturally occurs in the $p$-orbital band in optical lattices with ultra-cold spinless fermions. In the triangular lattice, it exhibits an unconventional supersolid state with the $f$-wave symmetry.
  • The tremendous development of cold atom physics has opened up a whole new opportunity to study novel states of matter which are not easily accessible in solid state systems. Here we propose to realize the $f$-wave pairing superfluidity of spinless fermions in the $p_{x,y}$-orbital bands of the two dimensional honeycomb optical lattices. The non-trivial orbital band structure rather than strong correlation effects gives rise to the unconventional pairing with the nodal lines of the $f$-wave symmetry. With a confining harmonic trap, zero energy Andreev bound states appear around the circular boundary with a six-fold symmetry. The experimental realization and detection of this novel pairing state are feasible.
  • The thermodynamic properties of the bilayer Sr$_3$Ru$_2$O$_7$ at very low temperatures are investigated by a realistic tight-binding model with the on-site interactions treated at the mean-field level. Due to the strong spin-orbit coupling, the band structure undergoes a significant change in Fermi surface topology as the external magnetic field is applied, invalidating the rigid band picture in which the Zeeman energy only causes chemical potential shifts. In addition, since Sr$_3$Ru$_2$O$_7$ is a $t_{2g}$ active system with unquenched orbital moments, the orbital Zeeman energy is not negligible and plays an important role in the phase diagram on the magnetic field orientation. We find that both the total density of states at the Fermi energy and the entropy exhibit a sudden increase near the critical magnetic field for the nematic phase, echoing the experimental findings. Our results suggest that extra cares are necessary to isolate the contributions due to the quantum criticality from the band structure singularity in this particular material. The effects of quantum critical fluctuations are briefly discussed.
  • Quasiparticle interference (QPI) in spectroscopic imaging scanning tunneling microscopy provides a powerful method to detect orbital band structures and orbital ordering patterns in transition metal oxides. We use the $T$-matrix formalism to calculate the QPI spectra for the unconventional metamagnetic system of Sr$_3$Ru$_2$O$_7$ with a $t_{2g}$-orbital band structure. A detailed tight-binding model is constructed accounting for features such as spin-orbit coupling, bilayer splitting, and the staggered rotation of the RuO octahedra. The band parameters are chosen by fitting the calculated Fermi surfaces with those measured in the angular-resolved photo-emission spectroscopy experiment. The calculated quasiparticle interference at zero magnetic field exhibits a hollow square-like feature arising from the nesting of the quasi-1d $d_{xz}$ and $d_{yz}$ orbital bands, in agreement with recent measurements by J. Lee {\it et al.} (Nature Physics {\bf 5}, 800 (2009)). Rotational symmetry breaking in the nematic metamagnetic state also manifests in the quasi-particle interference spectra.
  • The quasiparticle interference of the spectroscopic imaging scanning tunneling microscopy has been investigated for the surface states of the large gap topological insulator Bi$_2$Te$_3$ through the T-matrix formalism. Both the scalar potential scattering and the spin-orbit scattering on the warped hexagonal isoenergy contour are considered. While backscatterings are forbidden by time-reversal symmetry, other scatterings are allowed and exhibit strong dependence on the spin configurations of the eigenfunctions at k points over the isoenergy contour. The characteristic scattering wavevectors found in our analysis agree well with recent experiment results.