-
We investigate the normal state of the superconducting compound PuCoGa$_5$
using the combination of density functional theory (DFT) and dynamical mean
field theory (DMFT), with the continuous time quantum Monte Carlo (CTQMC) and
the vertex-corrected one-crossing approximation (OCA) as the impurity solvers.
Our DFT+DMFT(CTQMC) calculations suggest a strong tendency of Pu-5$f$ orbitals
to differentiate at low temperatures. The renormalized 5$f_{5/2}$ states
exhibit a Fermi-liquid behavior whereas one electron in the 5$f_{7/2}$ states
is at the edge of a Mott localization. We find that the orbital differentiation
is manifested as the removing of 5$f_{7/2}$ spectral weight from the Fermi
level relative to DFT. We corroborate these conclusions with DFT+DMFT(OCA)
calculations which demonstrate that 5$f_{5/2}$ electrons have a much larger
Kondo scale than the 5$f_{7/2}$.
-
Here we use angle-resolved photoemission spectroscopy to study
superconductivity that emerges in two extreme cases, from a Fermi liquid phase
(LiFeAs) and an incoherent bad-metal phase (FeTe0.55Se0.45). We find that
although the electronic coherence can strongly reshape the single-particle
spectral function in the superconducting state, it is decoupled from the
maximum superconducting pairing amplitude, which shows a universal scaling that
is valid for all FeSCs. Our observation excludes pairing scenarios in the BCS
and the BEC limit for FeSCs and calls for a universal strong coupling pairing
mechanism for the FeSCs.
-
We report a study of Co-doped berthierite Fe$_{1-x}$Co$_x$Sb$_2$S$_4$ ($x$ =
0, 0.1, and 0.2). The alloy series of Fe$_{1-x}$Co$_x$Sb$_2$S$_4$ crystallize
in an orthorhombic structure with the Pnma space group, similar to FeSb$_2$,
and show semiconducting behavior. The large discrepancy between activation
energy for conductivity, $E_\rho$ (146 $\sim$ 270 meV), and thermopower, $E_S$
(47 $\sim$ 108 meV), indicates the polaronic transport mechanism. Bulk
magnetization and heat capacity measurements of pure FeSb$_2$S$_4$ ($x$ = 0)
exhibit a broad antiferromagnetic (AFM) transition ($T_N$ = 46 K) followed by
an additional weak transition ($T^*$ = 50 K). Transition temperatures ($T_N$
and $T^*$) slightly decrease with increasing Co content $x$. This is also
reflected in the thermal conductivity measurement, indicating strong
spin-lattice coupling. Fe$_{1-x}$Co$_x$Sb$_2$S$_4$ shows relatively high value
of thermopower (up to $\sim$ 624 $\mu$V K$^{-1}$ at 300 K) and thermal
conductivity much lower when compared to FeSb$_{2}$, a feature desired for
potential applications based on FeSb$_{2}$ materials.
-
Bismuthates were the first family of oxide high-temperature superconductors,
exhibiting superconducting transition temperatures (Tc) up to 32K, but the
superconducting mechanism remains under debate despite more than 30 years of
extensive research. Our angle-resolved photoemission spectroscopy studies on
Ba$_{0.51}$K$_{0.49}$BiO$_3$ reveal an unexpectedly 34% larger bandwidth than
in conventional density functional theory calculations. This can be reproduced
by calculations that fully account for long-range Coulomb interactions --- the
first direct demonstration of bandwidth expansion due to the Fock exchange
term, a long-accepted and yet uncorroborated fundamental effect in many body
physics. Furthermore, we observe an isotropic superconducting gap with
2\Delta$_0$/k$_B$ T$_c$ = 3.51 $\pm$ 0.05, and strong electron-phonon
interactions with a coupling constant \lambda$\sim$ 1.3 $\pm$ 0.2. These
findings solve a long-standing mystery --- Ba$_{0.51}$K$_{0.49}$BiO$_3$ is an
extraordinary Bardeen-Cooper-Schrieffer (BCS) superconductor, where long-range
Coulomb interactions expand the bandwidth, enhance electron-phonon coupling,
and generate the high Tc. Such effects will also be critical for finding new
superconductors.
-
A repulsive Coulomb interaction between electrons in different orbitals in
correlated materials can give rise to bound quasiparticle states. We study the
non-hybridized two-orbital Hubbard model with intra (inter)-orbital interaction
$U$ ($U_{12}$) and different band widths using an improved dynamical mean field
theory numerical technique which leads to reliable spectra on the real energy
axis directly at zero temperature. We find that a finite density of states at
the Fermi energy in one band is correlated with the emergence of well defined
quasiparticle states at excited energies $\Delta=U-U_{12}$ in the other band.
These excitations are inter-band holon-doublon bound states. At the symmetric
point $U=U_{12}$, the quasiparticle peaks are located at the Fermi energy,
leading to a simultaneous and continuous Mott transition settling a
long-standing controversy.
-
Resonant elastic X-ray scattering is a powerful technique for measuring
multipolar order parameters. In this paper, we theoretically and experimentally
study the possibility of using this technique to detect the proposed multipolar
order parameters in URu$_2$Si$_2$ at the U-$L_{3}$ edge with the electric
quadrupolar transition. Based on an atomic model, we calculate the azimuthal
dependence of the quadrupolar transition at the U-$L_{3}$ edge. The results
illustrate the potential of this technique for distinguishing different
multipolar order parameters. We then perform experiments on ultra-clean single
crystals of URu$_2$Si$_2$ at the U-$L_{3}$ edge to search for the predicted
signal, but do not detect any indications of multipolar moments within the
experimental uncertainty. We theoretically estimate the orders of magnitude of
the cross-section and the expected count rate of the quadrupolar transition and
compare them to the dipolar transitions at the U-$M_4$ and U-$L_3$ edges,
clarifying the difficulty in detecting higher order multipolar order parameters
in URu$_2$Si$_2$ in the current experimental setup.
-
In this work we present a comparative investigation of the electronic
structures of NbO$_2$ and VO$_2$ obtained within the combination of density
functional theory and cluster-dynamical mean field theory calculations. We
investigate the role of dynamic electronic correlations on the electronic
structure of the metallic and insulating phases of NbO$_2$ and VO$_2$, with
focus on the mechanism responsible for the gap opening in the insulating
phases. For the rutile metallic phases of both oxides, we obtain that
electronic correlations lead to strong renormalization of the $t_{2g}$
subbands, as well as the emergence of incoherent Hubbard subbands, signaling
that electronic correlations are also important in the metallic phase of
NbO$_2$. Interestingly, we find that nonlocal dynamic correlations do play a
role in the gap formation of the (bct) insulating phase of NbO$_2$, by a
similar physical mechanism as that recently proposed by us in the case of the
(M$_1$) dimerized phase of VO$_2$ (\textit{Phys. Rev. Lett. 117, 056402
(2016)}). Although the effect of nonlocal dynamic correlations in the gap
opening of bct phase is less important than in the (M$_1$ and M$_2$) monoclinic
phases of VO$_2$, their presence indicates that the former is not a purely
Peierls-type insulator, as it was recently proposed.
-
A few approximate schemes to solve the Hedin equations self-consistently
introduced in (Phys. Rev. B 94, 155101 (2016)) are explored and tested for the
3D electron gas at metallic densities. We calculate one electron spectra,
dielectric properties, compressibility, and correlation energy. Considerable
reduction in the calculated band width (as compared to the self consistent GW
result) has been found when vertex correction was used for both polarizability
and self energy. Generally, it is advantageous to obtain the diagrammatic
representation of polarizability from the definition of this quantity as a
functional derivative of the electronic density with respect to the total field
(external plus induced). For self energy, the first order vertex correction
seems to be sufficient for the range of densities considered. Whenever it is
possible, we compare the accuracy of our vertex-corrected schemes with the
accuracy of the self-consistent quasi-particle GW approximation (QSGW), which
is less expensive computationally. We show that QSGW approach performs poorly
and we relate this poor performance with an inaccurate description of the
screening in the QSGW method (with an error comprising a factor 2-3 in the
physically important range of momenta).
-
We report on the emergence of robust superconducting order in single crystal
alloys of 2H-TaSe$_{2-x}$S$_{x}$ (0$\leq$x$\leq$2) . The critical temperature
of the alloy is surprisingly higher than that of the two end compounds
TaSe$_{2}$ and TaS$_{2}$. The evolution of superconducting critical temperature
T$_{c} (x)$ correlates with the full width at half maximum of the Bragg peaks
and with the linear term of the high temperature resistivity. The conductivity
of the crystals near the middle of the alloy series is higher or similar than
that of either one of the end members 2H-TaSe$_{2}$ and/or 2H-TaS$_{2}$. It is
known that in these materials superconductivity (SC) is in close competition
with charge density wave (CDW) order. We interpret our experimental findings in
a picture where disorder tilts this balance in favor of superconductivity by
destroying the CDW order.
-
The rich variety of iron-based superconductors and their complex electronic
structure lead to a wide range of possibilities for gap symmetry and pairing
components. Here we solve in the 2-Fe Brillouin zone the full
frequency-dependent linearized Eliashberg equations for LiFeAs with
spin-fluctuation mediated pairing interactions. The magnetic excitations are
calculated with the random phase approximation on a correlated electronic
structure obtained with density functional theory and dynamical mean field
theory. Correlations induce long-lived local moments with orbital-dependent
dynamics. The interaction between electrons through Hund's coupling promotes
inter-orbital magnetic susceptibility. As a consequence, the leading pairing
channel, conventional $s^{+-}$, acquires sizeable inter-orbital
$d_{xy}-d_{xz(yz)}$ singlet pairing with odd parity under glide-plane symmetry.
These components reduce the superconducting gap magnitude induced by the
intra-orbital components of the gap function at the electron pockets
intersection where the Fe-d $t_{2g}$ orbitals strongly mix. This in turn makes
the results consistent with available experiments on the angular dependence of
the gaps observed on the different Fermi surfaces.
-
In the iron-based superconductors (FeSCs), orbital differentiation is an
important phenomenon, whereby correlations stronger on the dxy orbital than on
the dxz/yz orbital yield quasi-particles with dxy} orbital character having
larger mass renormalization and abnormal temperature evolution. However, the
physical origin of this orbital differentiation is debated between the Hund's
coupling induced unbinding of spin and orbital degrees of freedom and the
Hubbard interaction instigated orbital selective Mott transition. Here we use
angle-resolved photoemission spectroscopy to identify an orbital-dependent
correlation-induced quasi-particle (QP) anomaly in LiFeAs. The excellent
agreement between our photoemission measurements and first-principles many-body
theory calculations shows that the orbital-differentiated QP lifetime anomalies
in LiFeAs are controlled by the Hund's coupling.
-
We present a theoretical investigation of the electronic structure of rutile
(metallic) and M$_1$ and M$_2$ monoclinic (insulating) phases of VO$_2$
employing a fully self-consistent combination of density functional theory and
embedded dynamical mean field theory calculations. We describe the electronic
structure of the metallic and both insulating phases of VO$_2$, and propose a
distinct mechanism for the gap opening. We show that Mott physics plays an
essential role in all phases of VO$_2$: undimerized vanadium atoms undergo
classical Mott transition through local moment formation (in the M$_2$ phase),
while strong superexchange within V-dimers adds significant dynamic intersite
correlations, which remove the singularity of self-energy for dimerized
V-atoms. The resulting transition from rutile to dimerized M$_1$ phase is
adiabatically connected to Peierls-like transition, but is better characterized
as the Mott transition in the presence of strong intersite exchange. As a
consequence of Mott physics, the gap in the dimerized M$_1$ phase is
temperature dependent. The sole increase of electronic temperature collapses
the gap, reminiscent of recent experiments.
-
We present optical transmission measurements that reveal a charge gap of 0.86
eV in the local moment antiferromagnetic insulator BaMn$_2$As$_2$, an order of
magnitude larger than previously reported. Density functional theory plus
dynamical mean field theory (DFT+DMFT) calculations correctly reproduce this
charge gap only when a strong Hund's coupling is considered. Thus,
BaMn$_2$As$_2$ is a member of a wider class of Mn pnictide compounds that are
Mott-Hund's insulators. We also present optical reflectance for metallic 2% K
doped BaMn$_2$As$_2$ that we use to extract the optical conductivity at
different temperatures. The optical conductivity $\sigma_1$($\omega$) exhibits
a metallic response that is well described by a simple Drude term. Both
$\sigma$($\omega$$\rightarrow$0, T) and $\rho$(T) exhibit Fermi liquid
temperature dependencies. From these measurements, we argue that a more
strongly correlated Hund's metal version of the parent compounds of the iron
pnictide superconductors has not yet been realized by doping this class of
Hund's insulators.
-
We analyze enhancements in the magnitude of the self-energy for electrons far
away from the Fermi surface in doped Mott insulators using the dynamical
cluster approximation to the Hubbard model. For large onsite repulsion, U, and
hole doping, the magnitude of the self-energy for imaginary frequencies at the
top of the band (k=(pi,pi) is enhanced with respect to the self-energy
magnitude at the bottom of the band (k=(0,0)). The self-energy behavior at
these two k-points is switched for electron doping. Although the hybridization
is much larger for (0,0) than for (pi,pi), we demonstrate that this is not the
origin of this difference. Isolated clusters under a downward shift of the
chemical potential, mu<U/2, at half-filling reproduce the overall self-energy
behavior at (0,0) and (pi,pi) found in low hole doped embedded clusters. This
happens although there is no change in the electronic structure of the isolated
clusters. Our analysis shows that a downward shift of the chemical potential
which weakly hole dopes the Mott insulator can lead to a large enhancement of
the (pi,pi) self-energy which is not necessarily associated with electronic
correlation effects, even in embedded clusters.
-
We show that the numerical renormalization group (NRG) is a viable multi-band
impurity solver for Dynamical Mean Field Theory (DMFT), offering unprecedent
real-frequency spectral resolution at arbitrarily low energies and
temperatures. We use it to obtain a numerically exact DMFT solution to the
Hund's metal problem for a three-orbital model with filling factor $n_d=2$. The
ground state is a Fermi liquid. The one-particle spectral function has a strong
particle-hole asymmetry, with a clear apparent power law for positive
frequencies only. With increasing temperature it shows a coherence-incoherence
crossover with spectral weight transfered from low to high energies and evolves
qualitatively differently from a doped Mott insulator. The spin and orbital
spectral functions show "spin-orbital separation": spin screening occurs at
much lower energies than orbital screening. The renormalization group flows
clearly reveal the relevant physics at all energy scales.
-
Ab initio calculations of the electric polarization of correlation-driven
insulating materials, namely Mott insulators, have not been possible so far.
Using a combination of density functional theory and dynamical-mean-field
theory we study the electric polarization of the Mott insulator
Sr$_{0.5}$Ba$_{0.5}$MnO$_{3}$. We predict a ferroelectric polarization of
$\simeq 16.5 \mu C/cm^2$ in the high temperature paramagnetic phase and recover
the measured value of $\simeq 13.3 \mu C/cm^2$ in the low temperature
antiferromagnetic phase. Our calculations reveal that the driving force for the
ferroelectricity, the hybridization between Mn e$_g$ and O p orbitals, is
suppressed by correlations, in particular by the Hund coupling and by the onset
of magnetic order. They also confirm that the half-filled Mn $t_{2g}$ orbitals
give rise to the antiferromagnetic Mott phase. This magnetic ordering leads to
changes in the ionic polar displacement and in turn to the electronic
polarization. In addition, for a fixed ionic displacement, we find that there
is a reduction in the electronic contribution due to partial magnetic
polarization of the Mn e$_g$ orbitals. The reduction of the polarization due to
ionic displacement dominates over the additional electronic part, hence the net
magneto-electric coupling is negative.
-
The electronic structure and numerous magnetic properties of MnBi magnetic
systems are investigated using local spin density approximation (LSDA) with
on-cite Coulomb correlations (LSDA+$U$) included. We show that the inclusion of
Coulomb correlations provides a much better description of equilibrium magnetic
moments on Mn atom as well as the magnetic anisotropy energy (MAE) behavior
with temperature and magneto-optical effects. We found that the inversion of
the anisotropic pairwise exchange interaction between Bi atoms is responsible
for the observed spin reorientation transition at 90 K. This interaction
appears as a result of strong spin orbit coupling on Bi atoms, large magnetic
moments on Mn atoms, significant $p-d$ hybridization between Mn and Bi atoms,
and it depends strongly on lattice constants. A better agreement with the
magneto-optical Kerr measurements at higher energies is obtained. We also
present the detailed investigation of the Fermi surface, the de Haas-van Alphen
(dHvA) effect and the X-ray magnetic circular dichroism in MnBi.
-
Spin-orbit coupling introduces chirality into electronic structure. This can
have profound effects on the magnetization induced by orbital motion of
electrons. Here we derive a formula for the orbital magnetization of
interacting electrons in terms of the full Green's function and vertex
functions. The formula is applied within dynamical mean-field theory to the
Kane-Mele-Hubbard model that allows both topological and trivial insulating
phases. We study the insulating and metallic phases in the presence of an
exchange magnetic field. In the presence of interactions, the orbital
magnetization of the quantum spin Hall insulating phase with inversion symmetry
is renormalized by the bulk quasi-particle weight. The orbital magnetization
vanishes for the in-plane antiferromagnetic phase with trivial topology. In the
metallic phase, the enhanced effective spin-orbit coupling due to the
interaction sometimes leads to an enhancement of the orbital magnetization.
However, at low doping, magnetization is suppressed at large interaction
strengths.
-
We report the first comprehensive study of the high temperature form
($\alpha$-phase) of iron disilicide. Measurements of the magnetic
susceptibility, magnetization, heat capacity and resistivity were performed on
well characterized single crystals. With a nominal iron $d^6$ configuration,
and a quasi-two dimensional crystal structure that strongly resembles that of
LiFeAs, $\alpha$-FeSi$_2$ is a potential candidate for unconventional
superconductivity. Akin to LiFeAs, $\alpha$-FeSi$_2$ does not develop any
magnetic order, and we confirm its metallic state down to the lowest
temperatures ($T$=1.8 K). However, our experiments reveal that paramagnetism
and electronic correlation effects in $\alpha$-FeSi$_2$ are considerably weaker
than in the pnictides. Band theory calculations yield small Sommerfeld
coefficients of the electronic specific heat $\gamma=C_e/T$ that are in
excellent agreement with experiment. Additionally, realistic many-body
calculations further corroborate that quasi-particle mass enhancements are only
modest in $\alpha$-FeSi$_{2}$ .
Remarkably, we find that the natural tendency to vacancy formation in the
iron sublattice has little influence on the iron valence and the density of
states at the Fermi level. Moreover, Mn doping does not significantly change
the electronic state of the Fe ion. This suggests that the iron valence is
protected against hole doping, and indeed the substitution of Co for Fe causes
a rigid-band like response of the electronic properties. As a key difference
from the pnictides, we identify the smaller inter-iron layer spacing, which
causes the active orbitals near the Fermi level to be of a different symmetry
in $\alpha$-FeSi$_2$. This change in orbital character might be responsible for
the lack of superconductivity in this system, providing constraints on pairing
theories in the iron based pnictides and chalcogenides.
-
We investigate signatures of electronic correlations in the narrow-gap
semiconductor FeGa$_3$ by means of electrical resistivity and thermodynamic
measurements performed on single crystals of FeGa$_3$, Fe$_{1-x}$Mn$_x$Ga$_3$
and FeGa$_{3-y}$Zn$_y$, complemented by a study of the 4$d$ analog material
RuGa$_3$. We find that the inclusion of sizable amounts of Mn and Zn dopants
into FeGa$_3$ does not induce an insulator-to-metal transition. Our study
indicates that both substitution of Zn onto the Ga site and replacement of Fe
by Mn introduces states into the semiconducting gap that remain localized even
at highest doping levels. Most importantly, using neutron powder diffraction
measurements, we establish that FeGa$_3$ orders magnetically above room
temperature in a complex structure, which is almost unaffected by the doping
with Mn and Zn. Using realistic many-body calculations within the framework of
dynamical mean field theory (DMFT), we argue that while the iron atoms in
FeGa$_3$ are dominantly in an $S=1$ state, there are strong charge and spin
fluctuations on short time scales, which are independent of temperature.
Further, the low magnitude of local contributions to the spin susceptibility
advocates an itinerant mechanism for the spin response in FeGa$_3$. Our joint
experimental and theoretical investigations classify FeGa$_3$ as a correlated
band insulator with only small dynamical correlation effects, in which
non--local exchange interactions are responsible for the spin gap of 0.4 eV and
the antiferromagnetic order. We show that hole doping of FeGa$_3$ leads, within
DMFT, to a notable strengthening of many--body renormalizations.
-
Using a combination of density functional theory and dynamical mean field
theory we show that electric polarization and magnetism are strongly
intertwined in (TMTTF)$_2$-$X$ (X$=$PF$_6$, As$F_6$, and SbF$_6$) organic
crystals and they originate from short-range Coulomb interactions. Electronic
correlations induce a charge-ordered state which, combined with the molecular
dimerization, gives rise to a finite electronic polarization and to a
ferroelectric state. We predict that the value of the electronic polarization
is enhanced by the onset of antiferromagnetism showing a sizable
magnetoelectric leading to a multiferroic behavior of (TMTTF)$_2$-$X$
compounds.
-
We present inelastic neutron scattering and magnetization measurements of the
antiferromagnetic insulator LaMnPO that are well described by a Heisenberg spin
model. These measurements are consistent with the presence of two-dimensional
magnetic correlations up to a temperature T$_{max}$ $\approx$ 700 K >> T$_{N}$
= 375 K, the N\'eel temperature. Optical transmission measurements show the T =
300 K direct charge gap $\Delta$ = 1 eV has decreased only marginally by 500 K
and suggest it decreases by only 10% at T$_{max}$. Density functional theory
and dynamical mean field theory calculations reproduce a direct charge gap in
paramagnetic LaMnPO only when a strong Hund's coupling J$_{H}$ = 0.9 eV is
included, as well as onsite Hubbard U = 8 eV. These results show the direct
charge gap in LaMnPO is rather insensitive to antiferromagnetic exchange
coupling and instead is a result of the local physics governed by U and
J$_{H}$.
-
The nature of the elemental cerium phases, undergoing an isostructural volume
collapse transition, cannot be understood using conventional solid state
concepts. Using the combination of density functional theory and dynamical
mean-field theory, we study the magnetic properties of both the $\alpha$ and
the $\gamma$ phases. We compute the magnetic form factor, and show that it is
very close to free ion behavior in both the local moment $\gamma$ phase as well
as the more itinerant $\alpha$ phase, in agreement with neutron scattering
experiments. In sharp contrast, the dynamic local magnetic susceptibility of
the two phases is strikingly different. In the $\gamma$ phase, the sharp low
energy peak due to local moment formation and consequently low Kondo
temperature dominates the spectra. In the $\alpha$ phase two broad peaks can be
identified, the first is due to Kondo screening, and the second is due to
Hund's coupling. This shows that hybridization plays a central role in the
$\alpha-\gamma$ transition in cerium, and that from the point of view of
magnetic properties, the $4f$ electrons are strongly correlated in both phases.
-
An understanding of the phase diagram of elemental plutonium (Pu) must
include both the effects of the strong directional bonding and the high density
of states of the Pu 5f electrons, as well as how that bonding weakens under the
influence of strong electronic correlations. We present for the first time
electronic-structure calculations of the full 16-atom per unit cell alpha-phase
structure within the framework of density functional theory (DFT) together with
dynamical mean-field theory (DMFT). Our calculations demonstrate that Pu atoms
sitting on different sites within the alpha-Pu crystal structure have a
strongly varying site dependence of the localization-delocalization correlation
effects of their 5f electrons and a corresponding effect on the bonding and
electronic properties of this complicated metal. In short, alpha-Pu has the
capacity to simultaneously have multiple degrees of electron
localization/delocalization of Pu 5f electrons within a pure single-element
material.
-
The symmetry of the wave function describing the Cooper pairs is one of the
most fundamental quantities in a superconductor but its measurement in the
iron-based superconductors has proved to be very difficult. The complex
multi-band nature of these materials makes the interplay of superconductivity
with spin and orbital dynamics very intriguing, leading to very material
dependent magnetic excitations, and pairing symmetries. Here we use
first-principles many-body method, including ab initio determined two-particle
vertex function, to study the spin dynamics and superconducting pairing
symmetry in a large number of iron-based superconductors. In iron compounds
with high transition temperature, we find both the dispersive high-energy spin
excitations, and very strong low energy commensurate or nearly commensurate
spin response, suggesting that these low energy spin excitations play the
dominate role in cooper pairing. We find three closely competing types of
pairing symmetries, which take a very simple form in the space of active iron
$3d$ orbitals, and differ only in the relative quantum mechanical phase of the
$xz$, $yz$ and $xy$ orbital contributions. The extensively discussed s$^{+-}$
symmetry appears when contributions from all orbitals have equal sign, while
the opposite sign in $xz$ and $yz$ orbitals leads to the $d$ wave symmetry. A
novel orbital antiphase $s^{+-}$ symmetry emerges when $xy$ orbital has
opposite sign to $xz$ and $yz$ orbitals. We propose that this orbital-antiphase
pairing symmetry explains the puzzling variation of the experimentally observed
superconducting gaps on all the Fermi surfaces of LiFeAs. This novel symmetry
of the order parameter may be realized in other iron superconductors.