
Physical properties of multiorbital materials depend not only on the
strength of the effective interactions among the valence electrons but also on
their type. Strong correlations are caused by either Mott physics that captures
the Coulomb repulsion among charges, or Hund physics that aligns the spins in
different orbitals. We identify four energy scales marking the onset and the
completion of screening in orbital and spin channels. The differences in these
scales, which are manifest in the temperature dependence of the local spectrum
and of the charge, spin and orbital susceptibilities, provide clear signatures
distinguishing Mott and Hund physics. We illustrate these concepts with
realistic studies of two archetypal strongly correlated materials, and
corroborate the generality of our conclusions with a model Hamiltonian study.

We study the fundamental question of the lattice dynamics of a metallic
ferromagnet in the regime where the static long range magnetic order is
replaced by the fluctuating local moments embedded in a metallic host. We use
the \textit{ab initio} Density Functional Theory(DFT)+embedded Dynamical
MeanField Theory(eDMFT) functional approach to address the dynamic stability
of iron polymorphs and the phonon softening with increased temperature. We show
that the nonharmonic and inhomogeneous phonon softening measured in iron is a
result of the melting of the long range ferromagnetic order, and is unrelated
to the first order structural transition from the BCC to the FCC phase, as is
usually assumed. We predict that the BCC structure is dynamically stable at all
temperatures at normal pressure, and is only thermodynamically unstable between
the BCC$\alpha$ and the BCC$\delta$ phase of iron.

Twoparticle Green's functions and the vertex functions play a critical role
in theoretical frameworks for describing strongly correlated electron systems.
However, numerical calculations at twoparticle level often suffer from large
computation time and massive memory consumption. We derive a general expansion
formula for the twoparticle Green's functions in terms of an overcomplete
representation based on the recently proposed "intermediate representation"
basis. The expansion formula is obtained by decomposing the spectral
representation of the twoparticle Green's function. We demonstrate that the
expansion coefficients decay exponentially, while all highfrequency and
longtail structures in the Matsubarafrequency domain are retained. This
representation therefore enables efficient treatment of twoparticle quantities
and opens a route to the application of modern manybody theories to realistic
strongly correlated electron systems.

Iridates abound with interesting magnetic behaviours because of their strong
spinorbit coupling. Sr$_3$NiIrO$_6$ brings together the spinorbital
entanglement of the Ir$^{4+}$ ion with a 3$d$ Ni cation and a onedimensional
crystal structure. It has a ferrimagnetic ground state with a 55 T coercive
field. We perform a theoretical study of the magnetic interactions in this
compound, and elucidate the role of anisotropic symmetric exchange as the
source of its strong magnetic anisotropy. Our firstprinciples calculations
reproduce the magnon spectra of this compound and predict a signature in the
cross sections that can differentiate the anisotropic exchange from singleion
anisotopy.

We report infrared studies of AFe$_{2}$As$_{2}$ (A = Ba, Sr), two
representative parent compounds of ironarsenide superconductors, at magnetic
fields (B) up to 17.5 T. Optical transitions between Landau levels (LLs) were
observed in the antiferromagnetic states of these two parent compounds. Our
observation of a $\sqrt{B}$ dependence of the LL transition energies, the
zeroenergy intercepts at B = 0 T under the linear extrapolations of the
transition energies and the energy ratio ($\sim$ 2.4) between the observed LL
transitions, combined with the linear band dispersions in twodimensional (2D)
momentum space obtained by theoretical calculations, demonstrates the existence
of massless Dirac fermions in antiferromagnetic BaFe$_{2}$As$_{2}$. More
importantly, the observed dominance of the zerothLLrelated absorption
features and the calculated bands with extremely weak dispersions along the
momentum direction $k_{z}$ indicate that massless Dirac fermions in
BaFe$_{2}$As$_{2}$ are 2D. Furthermore, we find that the total substitution of
the barium atoms in BaFe$_{2}$As$_{2}$ by strontium atoms not only maintains 2D
massless Dirac fermions in this system, but also enhances their Fermi velocity,
which supports that the Dirac points in ironarsenide parent compounds are
topologically protected.

The metalinsulator transition (MIT) remains among the most thoroughly
studied phenomena in solid state physics, but the complexity of the phenomena,
which usually involves cooperation of many degrees of freedom including
orbitals, fluctuating local moments, magnetism, and the crystal structure, have
resisted predictive abinitio treatment. Here we develop abinitio theoretical
method for correlated electron materials, based on Dynamical Mean Field Theory,
which can predict the change of the crystal structure across the MIT at finite
temperature. This allows us to study the coupling between electronic, magnetic
and orbital degrees of freedom with the crystal structure across the MIT in
rareearth nickelates. We predict the free energy profile of the competing
states, and the theoretical magnetic ground state configuration, which is in
agreement with neutron scattering data, but is different from the magnetic
models proposed before. The resonant elastic Xray response at the Kedge,
which was argued to be a direct probe of the charge order, is theoretically
modelled within the Dynamical Mean Field Theory, including the corehole
interaction. We show that the lineshape of the measured resonant elastic Xray
response can be explained with the "siteselective" Mott scenario without real
charge order on Ni sites.

We investigate the evolution of magnetic properties as a function of hydrogen
doping in iron based superconductor LaFeAsO$_{1x}$H$_x$ using the dynamical
meanfield theory combined with the densityfunctional theory. We find that two
independent consequences of the doping, the increase of the electron occupation
and the structural modification, have the opposite effects on the strength of
electron correlation and magnetism, resulting in the minimum of the calculated
magnetic moment around the intermediate doping level as a function of $x$. Our
result provides a natural explanation for the puzzling recent experimental
discovery of the two separated antiferromagnetic phases at low and high doping
limits. Furthermore, the increase of orbital occupation and correlation
strength with the doping results in reduced orbital polarization of $d_{xz/yz}$
orbitals and the enhanced role of $d_{xy}$ orbital in the magnetism at high
doping levels, and their possible implications to the superconductivity are
discussed in line with the essential role of the magnetism.

We implemented combination of DMFT and $GW$ in its fully selfconsistent way,
one shot $GW$ approximation, and quasiparticle selfconsistent scheme, and
studied how well these combined methods perform on H$_2$ molecule as compared
to more established methods such as LDA+DMFT. We found that most flavors of
$GW$+DMFT break down in strongly correlated regime due to causality violation.
Among $GW$+DMFT methods, only the selfconsistent quasiparticle $GW$+DMFT with
static doublecounting, and a new method with causal doublecounting, correctly
recover the atomic limit at large Hatom separation. While some flavors of
$GW$+DMFT improve the singleelectron spectra as compared to LDA+DMFT, the
total energy is best predicted by LDA+DMFT, for which the exact doublecounting
is known, and is static.

We use inelastic neutron scattering to study temperature and energy
dependence of spin excitations in optimally Pdoped BaFe2(As0.7P0.3)2
superconductor (Tc = 30 K) throughout the Brillouin zone. In the undoped state,
spin waves and paramagnetic spin excitations of BaFe2As2 stem from
antiferromagnetic (AF) ordering wave vector QAF= (1/1,0) and peaks near zone
boundary at (1/1,1/1) around 180 meV. Replacing 30% As by smaller P to induce
superconductivity, lowenergy spin excitations of BaFe2(As0.7P0.3)2form a
resonance in the superconducting state and highenergy spin excitations now
peaks around 220 meV near (1/1,1/1). These results are consistent with
calculations from a combined density functional theory and dynamical mean field
theory, and suggest that the decreased average pnictogen height in
BaFe2(As0.7P0.3)2 reduces the strength of electron correlations and increases
the effective bandwidth of magnetic excitations.

The theoretical understanding of emergent phenomena in quantum materials is
one of the greatest challenges in condensed matter physics. In contrast to
simple materials such as noble metals and semiconductors, macroscopic
properties of quantum materials cannot be predicted by the properties of
individual electrons. One of the examples of scientific importance is strongly
correlated electron system. Neither localized nor itinerant behaviors of
electrons in partially filled 3d, 4f, and 5f orbitals give rise to rich physics
such as Mott insulators, hightemperature superconductors, and superior
thermoelectricity, but hinder quantitative understanding of lowlying
excitation spectrum. Here we present a new firstprinciples approach to
strongly correlated solids. It is Q4 based on a combination of the
quasiparticle selfconsistent GW approximation and the dynamical meanfield
theory. The sole input in this method is the projector to the set of correlated
orbitals for which all local Feynman graphs are being evaluated. For that
purpose, we choose very localized quasiatomic orbitals spanning large energy
window, which contains most strongly hybridized bands, as well as upper and
lower Hubbard bands. The selfconsistency is carried out on the Matsubara axis.
This method enables the firstprinciples study of Mott insulators in both their
paramagnetic and antiferromagnetic phases. We illustrate the method on the
archetypical charge transfer correlated insulators La2CuO4 and NiO, and obtain
spectral properties and magnetic moments in good agreement with experiments.

The spin state transition in LaCoO$_3$ has eluded description for decades
despite concerted theoretical and experimental effort. In this study, we
approach this problem using charge consistent Density Functional Theory +
Embedded Dynamical Mean Field Theory (DFT+DMFT). We show, from first
principles, that LaCoO$_3$ cannot be described by a single, pure spin state at
any temperature, but instead shows a gradual change in the population of higher
spin multiples as temperature is increased. We explicitly elucidate the
critical role of the lattice expansion and oxygen octahedral rotations in the
spin state transition. We also show that the spin state transition and the
metalinsulator transition in LaCoO$_3$ occur at different temperature scales.
In addition, our results shed light on the importance of electronic entropy,
which has so far been ignored in all first principles studies of this material.

We introduce a cluster DMFT (Dynamical Mean Field Theory) approach to study
the normal state of the iron pnictides and chalcogenides. In the regime of
moderate mass renormalizations, the selfenergy is very local, justifying the
success of single site DMFT for these materials and for other Hunds metals. We
solve the corresponding impurity model with CTQMC (Continuous Time Quantum
MonteCarlo) and find that the minus sign problem is not severe in regimes of
moderate mass renormalization.

We use neutron scattering to study spin excitations in single crystals of
LiFe$_{0.88}$Co$_{0.12}$As, which is located near the boundary of the
superconducting phase of LiFe$_{1x}$Co$_{x}$As and exhibits nonFermiliquid
behavior indicative of a quantum critical point. By comparing spin excitations
of LiFe$_{0.88}$Co$_{0.12}$As with a combined density functional theory (DFT)
and dynamical mean field theory (DMFT) calculation, we conclude that
wavevector correlated low energy spin excitations are mostly from the $d_{xy}$
orbitals, while highenergy spin excitations arise from the $d_{yz}$ and
$d_{xz}$ orbitals. Unlike most iron pnictides, the strong orbital selective
spin excitations in LiFeAs family cannot be described by anisotropic Heisenberg
Hamiltonian. While the evolution of lowenergy spin excitations of
LiFe$_{1x}$Co$_x$As are consistent with electronhole Fermi surface nesting
condition for the $d_{xy}$ orbital, the reduced superconductivity in
LiFe$_{0.88}$Co$_{0.12}$As suggests that Fermi surface nesting conditions for
the $d_{yz}$ and $d_{xz}$ orbitals are also important for superconductivity in
iron pnictides.

We implemented the derivative of the free energy functional with respect to
the atom displacements, so called force, within the combination of Density
Functional Theory and the Embedded Dynamical Mean Field Theory. We show that in
combination with the numerically exact quantum Monte Carlo (MC) impurity
solver, the MC noise cancels to a great extend, so that the method can be used
very efficiently for structural optimization of correlated electron materials.
As an application of the method, we show how strengthening of the fluctuating
moment in FeSe superconductor leads to a substantial increase of the anion
height, and consequently to a very large effective mass, and also strong
orbital differentiation.

Magnetoelectric (ME) effect is recognized for its utility for lowpower
electronic devices.Largest ME coefficients are often associated with phase
transitions in which ferroelectricity is induced by magnetic order.
Unfortunately, in these systems, large ME response is revealed only upon
elaborate poling procedures. These procedures may become unnecessary in
singlepolardomain crystals of polar magnets. Here we report giant ME effects
in a polar magnet Fe2Mo3O8 at temperatures as high as 60 K.Polarization jumps
of 0.3 {\mu}C/cm2, and repeated mutual control of ferroelectric and magnetic
moments with differential ME coefficients on the order of 10$^4$ ps/m are
achieved.Importantly, no electric or magnetic poling is needed, as necessary
for applications.The sign of the ME coefficients can be switched by changing
the applied "bias" magnetic field. The observed effects are associated with a
hidden ferrimagnetic order unveiled by application of a magnetic field.

Dynamical mean field theory (DMFT) combined with the local density
approximation (LDA) is widely used in solids to predict properties of
correlated systems. In this paper, its application to one of the simplest
strongly correlated systems, the hydrogen molecule H$_2$, is demonstrated to
develop a parameterfree LDA+DMFT framework. We propose a method to calculate
the exact intersection of LDA and DMFT that leads to highly accurate
subtraction of the doubly counted correlation in both methods. When the exact
doublecounting treatment and a good projector to the correlated subspace are
used, LDA+DMFT yields very accurate total energy and excitation spectrum of the
H$_2$ molecule. We also discuss how this doublecounting scheme can be extended
to solid state calculations.

Combining density functional theory (DFT) and embedded dynamical meanfield
theory (DMFT) methods, we study the metalinsulator transition in
$R_2$Ir$_2$O$_7$ ($R$=Y, Eu, Sm, Nd, Pr, and Bi) and the topological nature of
the insulating compounds. Accurate free energies evaluated using the charge
selfconsistent DFT+DMFT method reveal that the metalinsulator transition
occurs for an Acation radius between that of Nd and Pr, in agreement with
experiments. The allinallout magnetic phase, which is stable in the Nd
compound but not the Pr one, gives rise to a small Ir$^{4+}$ magnetic moment of
$\approx 0.5\mu_B$ and opens a sizable correlated gap. We demonstrate that
within this stateoftheart theoretical method, the insulating bulk pyrochlore
iridates are topologically trivial.

We present a systematical theoretical study on the transport properties of an
archetypal family of Hund's metals, Sr$_2$RuO$_4$,Sr$_3$Ru$_2$O$_7$, SrRuO$_3$
and CaRuO$_3$, within the combination of first principles density functional
theory and dynamical mean field theory. The agreement between theory and
experiments for optical conductivity and resistivity is good, which indicates
that electronelectron scattering dominates the transport of ruthenates. We
demonstrate that in the singlesite dynamical mean field approach the transport
properties of Hund's metals fall into the scenario of "resilient
quasiparticles". We explains why the single layered compound Sr$_2$RuO$_4$ has
a relative weak correlation with respect to its siblings, which corroborates
its good metallicity.

Discovery of new transition metal compounds with large spin orbit coupling
(SOC) coexisting with strong electronelectron correlation among the $d$
electrons is essential for understanding the physics that emerges from the
interplay of these two effects. In this study, we predict a novel class of
$J_{eff}=1/2$ Mott insulators in a family of fluoride compounds that are
previously synthesized, but not characterized extensively. First principles
calculations in the level of all electron Density Functional Theory + Dynamical
Mean Field Theory (DFT+DMFT) indicate that these compounds have large Mott gaps
and some of them exhibit unprecedented proximity to the ideal, $SU(2)$
symmetric $J_{eff}=1/2$ limit.

The stationary functional of the allelectron density functional plus
dynamical mean field theory (DFT+DMFT) formalism to perform free energy
calculations and structural relaxations is implemented for the first time.
Here, the first order error in the density leads to a much smaller, second
order error in the free energy. The method is applied to several well known
correlated materials; metallic SrVO$_3$, Mott insulating FeO, and elemental
Cerium, to show that it predicts the lattice constants with very high accuracy.
In Cerium, we show that our method predicts the isostructural transition
between the $\alpha$ and $\gamma$ phases, and resolve the long standing
controversy in the driving mechanism of this transition.

We propose a continuum representation of the Dynamical Mean Field Theory, in
which we were able to derive an exact overlap between the Dynamical Mean Field
Theory and band structure methods, such as the Density Functional Theory. The
implementation of this exact doublecounting shows improved agreement between
theory and experiment in several correlated solids, such as the transition
metal oxides and lanthanides. Previously introduced nominal doublecounting is
in much better agreement with the exact doublecounting than most widely used
fully localized limit formula.

Layered 5d transition metal dichalcogenide (TMD) IrTe2 is distinguished from
the traditional TMDs (such as NbSe2) by the existence of multiple CDWlike
stripe phases and superconductivity at low temperatures. Despite of intensive
studies, there is still no consensus on the physical origin of the stripe
phases or even the ground state modulation for this 5d material. Here, we
present atomicscale evidence from scanning tunneling microscopy and
spectroscopy (STM/STS), that the ground state of IrTe2 is a q=1/6 stripe phase,
identical to that of the Sedoped compound. Furthermore, our data suggest that
the multiple transitions and stripe phases are driven by the intralayer IrIr
dimerization that competes against the interlayer TeTe bonding. The
competition results in a unified phase diagram with a series of hierarchical
modulated stripe phases, strikingly similar to the renowned "devil's staircase"
phenomena.

The solution of a generalized impurity model lies at the heart of electronic
structure calculations with dynamical meanfield theory (DMFT). In the
stronglycorrelated regime, the method of choice for solving the impurity model
is the hybridization expansion continuous time quantum Monte Carlo(CTHYB).
Enhancements to the CTHYB algorithm are critical for bringing new physical
regimes within reach of current computational power. Taking advantage of the
fact that the bottleneck in the algorithm is a product of hundreds of matrices,
we present optimizations based on the introduction and combination of two
concepts of more general applicability: a) skip lists and b) fast rejection of
proposed configurations based on matrix bounds. Considering two very different
test cases with $d$ electrons, we find speedups of $\sim 25$ up to $\sim 500$
compared to the direct evaluation of the matrix product. Even larger speedups
are likely with $f$ electron systems and with clusters of correlated atoms.

A combination of dynamical mean field theory and density functional theory,
as implemented in Phys. Rev. B 81, 195107 (2010), is applied to both the early
and late transition metal oxides. For fixed value of the local Coulomb
repulsion, without fine tuning, we obtain the main features of these series,
such as the metallic character of SrVO$_3$ and the the insulating gaps of
LaVO$_3$, LaTiO$_3$ and La$_2$CO$_4$ which are in good agreement with
experiment. The study highlights the importance of local physics and high
energy hybridization in the screening of the Hubbard interaction and how
different low energy behaviors can emerge from the a unified treatment of the
transition metal series.

Using tightbinding models and firstprinciples calculations, we demonstrate
the possibility to achieve a quantum anomalous Hall (QAH) phase on a
twodimensional square lattice, which can be realized in monolayers of double
perovskites. We show that effective intersite spinorbit coupling between eg
orbitals can be induced perturbatively, giving rise to a QAH state. Moreover,
the effective spinorbit coupling can be enhanced by octahedral rotations.
Based on firstprinciples calculations, we propose that this type of QAH state
could be realized in La2MnIrO6 monolayers, with the size of the gap as large as
26 meV in the ideal case. We observe that the electronic structure is sensitive
to structural distortions, and that an enhanced Hubbard U tends to stabilize
the nontrivial gap.