Charge order and broken rotational symmetry in magic-angle twisted bilayer graphene.

Bilayer graphene can be modified by rotating (twisting) one layer with respect to the other. The interlayer twist gives rise to a moiré superlattice that affects the electronic motion and alters the band structure1-4. Near a 'magic angle' of twist2,4, where the emergence of a flat band causes the charge carriers to slow down3, correlated electronic phases including Mott-like insulators and superconductors were recently discovered5-8 by using electronic transport. These measurements revealed an intriguing similarity between magic-angle twisted bilayer graphene and high-temperature superconductors, which spurred intensive research into the underlying physical mechanism9-14. Essential clues to this puzzle, such as the symmetry and spatial distribution of the spectral function, can be accessed through scanning tunnelling spectroscopy. Here we use scanning tunnelling microscopy and spectroscopy to visualize the local density of states and charge distribution in magic-angle twisted bilayer graphene. Doping the sample to partially fill the flat band, we observe a pseudogap phase accompanied by a global stripe charge order that breaks the rotational symmetry of the moiré superlattice. Both the pseudogap and the stripe charge order disappear when the band is either empty or full. The close resemblance to similar observations in high-temperature superconductors15-21 provides new evidence of a deeper link underlying the phenomenology of these systems.

Dynamic Response of an Electron Gas: Towards the Exact Exchange-Correlation Kernel.

Precise calculations of dynamics in the homogeneous electron gas (jellium model) are of fundamental importance for design and characterization of new materials. We introduce a diagrammatic Monte Carlo technique based on algorithmic Matsubara integration that allows us to compute frequency and momentum resolved finite temperature response directly in the real frequency domain using a series of connected Feynman diagrams. The data for charge response at moderate electron density are used to extract the frequency dependence of the exchange-correlation kernel at finite momenta and temperature. These results are as important for development of the time-dependent density functional theory for materials dynamics as ground state energies are for the density functional theory.

Interfacial charge-transfer Mott state in iridate-nickelate superlattices.

We investigate [Formula: see text]/[Formula: see text] superlattices in which we observe a full electron transfer at the interface from Ir to Ni, triggering a massive structural and electronic reconstruction. Through experimental characterization and first-principles calculations, we determine that a large crystal field splitting from the distorted interfacial [Formula: see text] octahedra surprisingly dominates over the spin-orbit coupling and together with the Hund's coupling results in the high-spin (S = 1) configurations on both the Ir and Ni sites. This demonstrates the power of interfacial charge transfer in coupling lattice, charge, orbital, and spin degrees of freedom, opening fresh avenues of investigation of quantum states in oxide superlattices.

Mott Metal-Insulator Transitions in Pressurized Layered Trichalcogenides.

Transition metal phosphorous trichalcogenides, MPX_{3} (M and X being transition metal and chalcogen elements, respectively), have been the focus of substantial interest recently because they are unusual candidates undergoing Mott transition in the two-dimensional limit. Here we investigate material properties of the compounds with M=Mn and Ni employing ab initio density functional and dynamical mean-field calculations, especially their electronic behavior under external pressure in the paramagnetic phase. Mott metal-insulator transitions (MIT) are found to be a common feature for both compounds, but their lattice structures show drastically different behaviors depending on the relevant orbital degrees of freedom, i.e., t_{2g} or e_{g}. Under pressure, MnPS_{3} can undergo an isosymmetric structural transition within monoclinic space group by forming Mn-Mn dimers due to the strong direct overlap between the neighboring t_{2g} orbitals, accompanied by a significant volume collapse and a spin-state transition. In contrast, NiPS_{3} and NiPSe_{3}, with their active e_{g} orbital degrees of freedom, do not show a structural change at the MIT pressure or deep in the metallic phase within the monoclinic symmetry. Hence NiPS_{3} and NiPSe_{3} become rare examples of materials hosting electronic bandwidth-controlled Mott MITs, thus showing promise for ultrafast resistivity switching behavior.

Phonon Softening due to Melting of the Ferromagnetic Order in Elemental Iron.

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 ab initio density functional theory + embedded dynamical mean-field theory functional approach to address the dynamic stability of iron polymorphs and the phonon softening with an 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 thermodynamically unstable only between the bcc-α and the bcc-δ phases of iron.

Metal-Insulator Transition and Topological Properties of Pyrochlore Iridates.

Combining density functional theory (DFT) and embedded dynamical mean-field theory (DMFT) methods, we study the metal-insulator 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 self-consistent DFT+DMFT method reveal that the metal-insulator transition occurs for an A-cation radius between that of Nd and Pr, in agreement with experiments. The all-in-all-out magnetic phase, which is stable in the Nd compound but not the Pr one, gives rise to a small Ir^{4+} magnetic moment of ≈0.4 µ_{B} and opens a sizable correlated gap. We demonstrate that within this state-of-the-art theoretical method, the insulating bulk pyrochlore iridates are topologically trivial.

Two-Dimensional Massless Dirac Fermions in Antiferromagnetic AFe_{2}As_{2} (A=Ba,Sr).

We report infrared studies of AFe_{2}As_{2} (A=Ba, Sr), two representative parent compounds of iron-arsenide 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 zero-energy intercepts at B=0 T under the linear extrapolations of the transition energies and the energy ratio (∼2.4) between the observed LL transitions, combined with the linear band dispersions in two-dimensional (2D) momentum space obtained by theoretical calculations, demonstrates the existence of massless Dirac fermions in the antiferromagnet BaFe_{2}As_{2}. More importantly, the observed dominance of the zeroth-LL-related 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 iron-arsenide parent compounds are topologically protected.

Transport Properties of Metallic Ruthenates: A DFT+DMFT Investigation.

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 electron-electron scattering dominates the transport of ruthenates. We demonstrate that in the single-site dynamical mean field approach the transport properties of Hund's metals fall into the scenario of "resilient quasiparticles." We explain why the single layered compound Sr_{2}RuO_{4} has a relative weak correlation with respect to its siblings, which corroborates its good metallicity.

Orbital Selective Spin Excitations and their Impact on Superconductivity of LiFe_{1-x}Co_{x}As.

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_{1-x}Co_{x}As and exhibits non-Fermi-liquid 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 and dynamical mean field theory calculation, we conclude that wave-vector correlated low energy spin excitations are mostly from the d_{xy} orbitals, while high-energy spin excitations arise from the d_{yz} and d_{xz} orbitals. Unlike most iron pnictides, the strong orbital selective spin excitations in the LiFeAs family cannot be described by an anisotropic Heisenberg Hamiltonian. While the evolution of low-energy spin excitations of LiFe_{1-x}Co_{x}As is consistent with the electron-hole Fermi surface nesting conditions 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.

Exact Double Counting in Combining the Dynamical Mean Field Theory and the Density Functional Theory.

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; double counting. The implementation of this exact double counting shows improved agreement between the theory and experiment in several correlated solids, such as the transition metal oxides and lanthanides. Previously introduced nominal double counting is in much better agreement with the exact double counting than the most widely used fully localized limit formula.

Free Energy from Stationary Implementation of the DFT+DMFT Functional.

The stationary functional of the density functional plus embedded dynamical mean field theory 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 good accuracy. In cerium, we show that our method predicts the isostructural transition between the α and γ phases, and resolve the long-standing controversy in the driving mechanism of this transition.

J(eff)=1/2 Mott-insulating state in Rh and Ir fluorides.

Discovery of new transition metal compounds with large spin orbit coupling coexisting with strong electron-electron 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 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.

Signatures of electronic correlations in iron silicide.

The intermetallic FeSi exhibits an unusual temperature dependence in its electronic and magnetic degrees of freedom, epitomized by the cross-over from a low-temperature nonmagnetic semiconductor to a high-temperature paramagnetic metal with a Curie-Weiss-like susceptibility. Many proposals for this unconventional behavior have been advanced, yet a consensus remains elusive. Using realistic many-body calculations, we here reproduce the signatures of the metal-insulator cross-over in various observables: the spectral function, the optical conductivity, the spin susceptibility, and the Seebeck coefficient. Validated by quantitative agreement with experiment, we then address the underlying microscopic picture. We propose a new scenario in which FeSi is a band insulator at low temperatures and is metalized with increasing temperature through correlation induced incoherence. We explain that the emergent incoherence is linked to the unlocking of iron fluctuating moments, which are almost temperature independent at short timescales. Finally, we make explicit suggestions for improving the thermoelectric performance of FeSi based systems.

Shining light on transition-metal oxides: unveiling the hidden Fermi liquid.

We use low energy optical spectroscopy and first principles local density approximation plus dynamical mean field theory calculations to test the hypothesis that the anomalous transport properties of strongly correlated metals originate in the strong temperature dependence of their underlying resilient quasiparticles. We express the resistivity in terms of an effective plasma frequency ω(p)* and an effective scattering rate 1/τ(tr)*. We show that in the archetypal correlated material V2O3, ω(p)* increases with increasing temperature, while the plasma frequency from the partial sum rule exhibits the opposite trend. 1/τ(tr)* has a more pronounced temperature dependence than the scattering rate obtained from the extended Drude analysis. The theoretical calculations of these quantities are in quantitative agreement with experiment. We conjecture that these are robust properties of all strongly correlated metals, and test the conjecture by carrying out a similar analysis on thin film NdNiO3 on a LaAlO3 substrate.

Quantum critical phase of FeO spans conditions of Earth's lower mantle.

Seismic and mineralogical studies have suggested regions at Earth's core-mantle boundary may be highly enriched in FeO, reported to exhibit metallic behavior at extreme pressure-temperature (P-T) conditions. However, underlying electronic processes in FeO remain poorly understood. Here we explore the electronic structure of B1-FeO at extreme conditions with large-scale theoretical modeling using state-of-the-art embedded dynamical mean field theory (eDMFT). Fine sampling of the phase diagram reveals that, instead of sharp metallization, compression of FeO at high temperatures induces a gradual orbitally selective insulator-metal transition. Specifically, at P-T conditions of the lower mantle, FeO exists in an intermediate quantum critical state, characteristic of strongly correlated electronic matter. Transport in this regime, distinct from insulating or metallic behavior, is marked by incoherent diffusion of electrons in the conducting t2g orbital and a band gap in the eg orbital, resulting in moderate electrical conductivity (~105 S/m) with modest P-T dependence as observed in experiments. Enrichment of solid FeO can thus provide a unifying explanation for independent observations of low seismic velocities and elevated electrical conductivities in heterogeneities at Earth's mantle base.

Hidden Fermi liquid, scattering rate saturation, and Nernst effect: a dynamical mean-field theory perspective.

We investigate the transport properties of a correlated metal within dynamical mean-field theory. Canonical Fermi liquid behavior emerges only below a very low temperature scale T(FL). Surprisingly the quasiparticle scattering rate follows a quadratic temperature dependence up to much higher temperatures and crosses over to saturated behavior around a temperature scale T(sat). We identify these quasiparticles as constituents of the hidden Fermi liquid. The non-Fermi-liquid transport above T(FL), in particular the linear-in-T resistivity, is shown to be a result of a strongly temperature dependent band dispersion. We derive simple expressions for the resistivity, Hall angle, thermoelectric power and Nernst coefficient in terms of a temperature dependent renormalized band structure and the quasiparticle scattering rate. We discuss possible tests of the dynamical mean-field theory picture of transport using ac measurements.

Plutonium hexaboride is a correlated topological insulator.

We predict that plutonium hexaboride (PuB(6)) is a strongly correlated topological insulator, with Pu in an intermediate valence state of Pu(2.7+). Within the combination of dynamical mean field theory and density functional theory, we show that PuB(6) is an insulator in the bulk, with nontrivial Z(2) topological invariants. Its metallic surface states have a large Fermi pocket at the X[over ¯] point and the Dirac cones inside the bulk derived electronic states, causing a large surface thermal conductivity. PuB(6) has also a very high melting temperature; therefore, it has ideal solid state properties for a nuclear fuel material.

Effective J=1/2 insulating state in Ruddlesden-Popper iridates: an LDA+DMFT study.

Using ab initio methods for correlated electrons in solids, we investigate the metal-insulator transition across the Ruddlesden-Popper (RP) series of iridates and explore the robustness of the Jeff=1/2 state against band effects due to itineracy, tetragonal distortion, octahedral rotation, and Coulomb interaction. We predict the effects of epitaxial strain on the optical conductivity, magnetic moments, and Jeff=1/2 ground-state wave functions in the RP series. To describe the solution of the many-body problem in an intuitive picture, we introduce a concept of energy-dependent atomic states, which strongly resemble the atomic Jeff=1/2 states but with coefficients that are energy or time dependent. We demonstrate that the deviation from the ideal Jeff=1/2 state is negligible at short time scales for both single- and double-layer iridates, while it becomes quite significant for Sr3Ir2O7 at long times and low energy. Interestingly, Sr2IrO4 is positioned very close to the SU(2) limit, with only ∼3% deviation from the ideal Jeff=1/2 situation.

γ-α isostructural transition in cerium.

We present zero-temperature first-principles calculations of elemental cerium and we compute its pressure-volume phase diagram within a theoretical framework able to describe simultaneously both the α and the γ phases. A surprising result revealed by our study is the presence of a clear signature of the transition at zero temperature and that this signature can be observed if and only if the spin-orbit coupling is taken into account. Our calculations indicate that the transition line in the pressure-temperature phase diagram of this material has a low-T critical point at negative pressures, placed very close to zero temperature. This suggests that cerium is very close to being "quantum critical," in agreement with recent experiments.

Ground state of the parallel double quantum dot system.

We resolve the controversy regarding the ground state of the parallel double quantum dot system near half filling. The numerical renormalization group predicts an underscreened Kondo state with residual spin-1/2 magnetic moment, ln2 residual impurity entropy, and unitary conductance, while the Bethe ansatz solution predicts a fully screened impurity, regular Fermi-liquid ground state, and zero conductance. We calculate the impurity entropy of the system as a function of the temperature using the hybridization-expansion continuous-time quantum Monte Carlo technique, which is a numerically exact stochastic method, and find excellent agreement with the numerical renormalization group results. We show that the origin of the unconventional behavior in this model is the odd-symmetry "dark state" on the dots.