
The emergent dynamical processes associated with the magnetic excitations in
heavy fermion SmB$_6$ are investigated. By imposing multiorbital 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
spinspin 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 slavespin method is introduced to systematically investigate the
solutions of the Hubbard model including the symmetrybroken phases. In this
method, the electron operator is factorized into a fermioninc spinon describing
the physical spin and a slavespin 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 selfconsistent equations to explore various
symmetrybroken density wave states can be constructed in general with a
cluster of multiple slavespin sites. We employ this method to study the
antiferromagnetic (AFM) state in the single band Hubbard model with the two and
foursite clusters of slave spins. While the Hubbard gap, the charge gap due to
the doublyoccupied 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 strongcoupling behaviors as $U$
increases. Our cluster slavespin method reproduces not only the traditional
meanfield behavior of $\Delta \sim U$ in the weakcoupling limit, but also the
behavior of $\Delta \sim t^2/U$ predicted by the superexchange mechanism in the
strong coupling limit. In addition, the holondoublon correlator as functions
of $U$ and doping $x$ is also computed, which exhibits a strong tendency toward
the holondoublon binding in the strong coupling regime. We further show that
the quasiparticle weight obtained by the cluster slavespin 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 slavespin 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 gaplike feature in the
particlehole 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 densitydensity correlation function at q~0 and q~Q_(CDW).
Based on our analysis, we propose that measuring these resonant plasmonaxion
excitations around q~0 and q~Q_(CDW) by the momentumresolved electron energy
loss spectroscopy (MEELS) 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 nonsuperconducting 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 nonsuperconducting 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 tunnelingbased probes. Moreover, we will show how
the SchwingerKadanoffBaymKeldysh (SKBK) formalism together with the
appropriate modeling of the nanoscale 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 $tJU$ model within a slavespin
formalism. We show that the BCS meanfield 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 halffilling, $Z$ develops a strong doping dependence,
leading to a dopingdriven 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 meanfield
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 caxis 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 straininduced modulation of electron correlation
effects and increased orbital anisotropy in the rutile phase of epitaxial
VO$_2$/TiO$_2$ films from hard xray photoelectron spectroscopy and soft V
Ledge xray 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 noninteger, but close to the halffilling. Because the
overlaps of wave functions between $d$ orbitals are modified by the strain,
orbitallydependent 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 noninteger filling number near the halffilling such as for
VO$_2$, certain orbitals could reach an occupation number closer to
halffilling 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 strainengineering of
correlated electron systems.

The quasiparticle interferences (QPIs) of a Mott insulator are investigated
using the $T$matrix formalism implemented with the dynamical meanfield 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 energydependent 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 noninteracting
Fermi surfaces appears at this critical bias voltage $eV_c$ in Mott insulating
state. We propose that this reentry of noninteracting 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 midinfrared scenario for hightemperature superconductivity
proposed by Leggett, the effects of Umklapp processes on the densitydensity
correlation function in the presence of longrange 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 midinfrared
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 pointcontact conductance between a
metalic electrode and a strongly correlated material using the nonequilibrium
SchwingerKadanoffBaymKeldysh 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 inplane momentum. We find that the conductance $dI/dV$
is proportional to the {\it effective density of states}, that is, the
integrated singleparticle spectral function $A(\omega=eV)$ over the whole
Brillouin zone. From this conclusion, we are able to establish the conditions
under which a nonFermi liquid metal exhibits a zerobias peak in the
conductance. This finding is discussed in the context of recent pointcontact
spectroscopy on the iron pnictides and chalcogenides which has exhibited a
zerobias conductance peak.

Using a degenerate doubleexchange 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 randomphase
approximation calculation on both spinons and itinerant electrons. We find that
in the $\left(\pi,0\right)$ magneticallyordered state, the spinwave
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 HolsteinPrimakoff transformation. In the paramagnetic state, the
particlehole 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 particlehole 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 dwave 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$functionlike 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
ironbased 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 spinbased scenarios have been
proposed as the driving force. In this review, we focus on theories from the
orbitalbased 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 longstanding debate between orbital and
spinbased theories.

The critical behavior of a twoorbital 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
twoorbital model and becomes the dominant fluctuation in the vicinity of the
orbitalordering quantum critical point. Since the very existence of this $z=3$
overdamped collective mode induces nonFermi liquid behavior near the quantum
critical point, we conclude that a twoorbital model generally has a sizable
area in the phase diagram showing nonFermi liquid behavior. Furthermore, we
show that the bosonic theory resembles the continuous model near the dwave
Pomeranchuk instability, suggesting that orbital order in a twoorbital 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 ironbased
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 ironpnictide superconductors using a five orbital tight
binding model with generalized Hubbard onsite interactions. Within a oneloop
treatment, we find that an overdamped collective mode develops at low frequency
in channels associated with quasi1D $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 nonFermi 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 nonFermi liquid behavior is responsible for the recently observed
zerobias enhancement in the tunneling signal in quantum point contact
spectroscopy. A key experimental test of this proposal is the absence of the
nonFermi liquid behaviour in the holedoped materials. Our result suggests
that quantum criticality plays an important role in understanding the normal
state properties of ironpnictide superconductors.

We show here that orbital degrees of freedom produce a distinct signature in
the magnetic excitation spectrum of ironbased 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 fiveorbital 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 longrange orbital order is absent,
Gaussian fluctuations beyond the standard randomphase 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 incommensuratetocommensurate 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
twoband superconductor with s+ type of pairing. This mechanism exploits the
proximity effect with a conventional swave superconductor which favors the
alignment of the phases of the twoband superconductor which, in the case of
s+ pairing, are pishifted 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 solitonfree 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 ironbased 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$_{2x}$Se$_2$
using a threeorbital ($t_{2g}$) tightbinding 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 highspin to lowspin states
accompanied by a crossover from orbitallydisordered to orbitallyordered
states will be driven by doping the parent compound with electrons, which can
be verified by neutron scattering and soft Xray measurements.

We develop here a general formalism for multiorbital 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 multiorbital structure. As a consequence
certainly dynamically generated symmetries obtain at lower values of doping
than in the singleband Hubbard model. For example, in the atomic limit, the
particlehole symmetric condition in the lower band shifts from the oneband
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 lowenergy
theories which forbid double occupancy of bare electrons, such as the
multiorbital tJ generalization, should all be be sensitive to this
particlehole symmetric condition. Away from the atomic limit, the dynamical
contributions increase the transferred spectral weight. Consequently, any
phenomena which are sensitive to an emergent particlehole 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 spinpolaron 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 highspin 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 nonbipartite lattices which leads to frustrated supersolid
states with nonuniform distributions of the Cooper pair phase and density.
This exotic pairing state naturally occurs in the $p$orbital band in optical
lattices with ultracold 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 nontrivial 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
sixfold 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 tightbinding model with the
onsite interactions treated at the meanfield level. Due to the strong
spinorbit 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
tightbinding model is constructed accounting for features such as spinorbit
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 angularresolved photoemission spectroscopy experiment.
The calculated quasiparticle interference at zero magnetic field exhibits a
hollow squarelike feature arising from the nesting of the quasi1d $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 quasiparticle
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 Tmatrix formalism. Both the
scalar potential scattering and the spinorbit scattering on the warped
hexagonal isoenergy contour are considered. While backscatterings are forbidden
by timereversal 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.