Source author record

Andrew J. Millis

Andrew J. Millis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.str-el cond-mat.mtrl-sci cond-mat.mes-hall cond-mat.supr-con cond-mat.stat-mech quant-ph cond-mat.other Machine Learning cond-mat.dis-nn physics.optics physics.chem-ph physics.data-an

Catalog footprint

What is connected

107works

12topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

Dynamical instability in a Floquet-Driven Dissipative System

We analyse the magnon spectrum and distribution function of the antiferromagnetic phase of the Floquet-driven Hubbard model. Above a critical drive strength, we find a dynamical instability, resulting from a change in sign of the magnon damping at a non-zero wavevector. The change in sign means that infinitesimal fluctuations grow with time, corresponding to an instability of the driven state. Implications for the nonequilibrium distribution function and the strong drive nonlinear dynamics are discussed.

preprint2025arXiv

Superfluid stiffness bounds in time-reversal symmetric superconductors

Quantum geometry has been shown to make an important contribution to the superfluid stiffness of superconductors, especially for flat-band systems such as moiré materials. In this work we use mean-field theory to derive an expression for the superfluid stiffness of time-reversal symmetric superconductors at zero temperature by computing the energy of the mean-field ground state as a function of pairing momentum. We show that the quantum geometric contribution to superfluid stiffness is a consequence of broken Galilean invariance in the interaction Hamiltonian, arising from momentum-dependent form factors related to the momentum dependence of Bloch states. The effects of broken Galilean invariance are not fully parametrized by the quantum metric considered in previous work. We obtain general lower and upper bounds that apply to both continuum and lattice models and present numerical calculations of the precise value in several important cases. The superfluid stiffness of superconductivity in a Landau level saturates the lower bound and the superfluid stiffness of the other cases we consider is close to the general lower bound we derive. In multilayer rhombohedral graphene the geometric contribution is shown not to be the dominant contribution to the superfluid stiffness, despite the flat band behavior in the vicinity of the Fermi level. Finally, assuming contact interaction and uniform pairing, we show that the superfluid stiffness is proportional to the ``minimal quantum metric" introduced in previous work. We provide a continuum version of the minimal quantum metric and explain its physical origin.

preprint2022arXiv

Atomically imprinted graphene plasmonic cavities

Plasmon polaritons in van der Waals (vdW) materials hold promise for next-generation photonics. The ability to deterministically imprint spatial patterns of high carrier density in cavities and circuitry with nanoscale features underlies future progress in nonlinear nanophotonics and strong light-matter interactions. Here, we demonstrate a general strategy to atomically imprint low-loss graphene plasmonic structures using oxidation-activated charge transfer (OCT). We cover graphene with a monolayer of WSe$_2$, which is subsequently oxidized into high work-function WOx to activate charge transfer. Nano-infrared imaging reveals low-loss plasmon polaritons at the WOx/graphene interface. We insert WSe$_2$ spacers to precisely control the OCT-induced carrier density and achieve a near-intrinsic quality factor of plasmons. Finally, we imprint canonical plasmonic cavities exhibiting laterally abrupt doping profiles with single-digit nanoscale precision via programmable OCT. Specifically, we demonstrate technologically appealing but elusive plasmonic whispering-gallery resonators based on free-standing graphene encapsulated in WOx. Our results open avenues for novel quantum photonic architectures incorporating two-dimensional materials.

preprint2022arXiv

Common microscopic origin of the phase transitions in Ta$_2$NiS$_5$ and the excitonic insulator candidate Ta$_2$NiSe$_5$

The structural phase transition in Ta$_2$NiSe$_5$ has been envisioned as driven by the formation of an excitonic insulating phase. However, the role of structural and electronic instabilities on crystal symmetry breaking has yet to be disentangled. Meanwhile, the phase transition in its complementary material Ta$_2$NiS$_5$ does not show any experimental hints of an excitonic insulating phase. We present a microscopic investigation of the electronic and phononic effects involved in the structural phase transition in Ta$_2$NiSe$_5$ and Ta$_2$NiS$_5$ using extensive first-principles calculations. In both materials the crystal symmetries are broken by phonon instabilities, which in turn lead to changes in the electronic bandstructure also observed in experiment. A total energy landscape analysis shows no tendency towards a purely electronic instability and we find that a sizeable lattice distortion is needed to open a bandgap. We conclude that an excitonic instability is not needed to explain the phase transition in both Ta$_2$NiSe$_5$ and Ta$_2$NiS$_5$.

preprint2022arXiv

Dynamical Mean Field Theory of Moiré Bilayer Transition Metal Dichalcogenides: Phase Diagram, Resistivity, and Quantum Criticality

We present a comprehensive dynamical mean field study of the triangular lattice moiré Hubbard model, which is believed to represent the physics of moiré bilayer transition metal dichalcogenides. In these materials, important aspects of the band structure including the bandwidth and the order and location of van Hove singularities can be tuned by varying the interlayer potential. We present a magnetic and metal-insulator phase diagram and a detailed study of the dependence of the resistivity on temperature, band filling and interlayer potential. We find that transport displays Fermi liquid, strange metal and quantum critical behaviors in distinct regions of the phase diagram. Specifically, we find that the cube-root van Hove singularity ($ρ(ε) \sim|ε|^{-1 / 3}$) gives a strange metal behavior with a $T$-linear scattering rate and $ω/T$ scaling. We show how magnetic order affects the resistivity. Our results elucidate the physics of the correlated states and the metal-insulator continuous transition recently observed in twisted homobilayer WSe$_2$ and heterobilayer MoTe$_2$/WSe$_2$ experiments.

preprint2022arXiv

Machine Learning for Optical Scanning Probe Nanoscopy

The ability to perform nanometer-scale optical imaging and spectroscopy is key to deciphering the low-energy effects in quantum materials, as well as vibrational fingerprints in planetary and extraterrestrial particles, catalytic substances, and aqueous biological samples. The scattering-type scanning near-field optical microscopy (s-SNOM) technique has recently spread to many research fields and enabled notable discoveries. In this brief perspective, we show that the s-SNOM, together with scanning probe research in general, can benefit in many ways from artificial intelligence (AI) and machine learning (ML) algorithms. We show that, with the help of AI- and ML-enhanced data acquisition and analysis, scanning probe optical nanoscopy is poised to become more efficient, accurate, and intelligent.

preprint2022arXiv

Nonequilibrium phase transition in a driven-dissipative quantum antiferromagnet

A deeper theoretical understanding of driven-dissipative interacting systems and their nonequilibrium phase transitions is essential both to advance our fundamental physics understanding and to harness technological opportunities arising from optically controlled quantum many-body states. This paper provides a numerical study of dynamical phases and the transitions between them in the nonequilibrium steady state of the prototypical two-dimensional Heisenberg antiferromagnet with drive and dissipation. We demonstrate a nonthermal transition that is characterized by a qualitative change in the magnon distribution, from subthermal at low drive to a generalized Bose-Einstein form including a nonvanishing condensate fraction at high drive. A finite-size analysis reveals static and dynamical critical scaling at the transition, with a discontinuous slope of the magnon number versus driving field strength and critical slowing down at the transition point. Implications for experiments on quantum materials and polariton condensates are discussed.

preprint2022arXiv

Optical Conductivity of the two dimensional Hubbard model: vertex corrections, emergent Galilean invariance and the accuracy of the single-site dynamical mean field approximation

We compute the frequency dependent conductivity of the two dimensional square lattice Hubbard model at zero temperature as a function of density to second order in the interaction strength, and compare the results to the predictions of single-site dynamical mean field theory computed at the same order. We find that despite the neglect of vertex corrections, the single site dynamical mean field approximation produces semiquantitatively accurate results for most carrier concentrations, but fails qualitatively for the nearly empty or nearly filled band cases where the theory exhibits an emergent Galilean invariance. The theory also becomes qualitatively inaccurate very near half filling if nesting is important.

preprint2022arXiv

Second-order magnetic responses in quantum magnets: Magnetization under ac magnetic fields

We investigate second-order magnetic responses of quantum magnets against ac magnetic fields. We focus on the case where the $z$ component of the spin is conserved in the unperturbed Hamiltonian and the driving field is applied in the $xy$ plane. We find that linearly polarized driving fields induce a second-harmonic response, while circularly polarized fields generate only a zero-frequency response, leading to a magnetization with a direction determined by the helicity. Employing an unbiased numerical method, we demonstrate the nonlinear magnetic effect driven by the circularly polarized field in the XXZ model and show that the magnitude of the magnetization can be predicted by the dynamical spin structure factor in the linear response regime.

preprint2022arXiv

Superconductivity and Antiferromagnetism in NdNiO$_2$ and CaCuO$_2$: A Cluster DMFT Study

We perform a comparative $2 \times 2$ real space cluster DMFT study on minimal models for NdNiO$_2$ and CaCuO$_2$ obtained from downfolding DFT states, using a Nambu formalism that allows for both superconducting and antiferromagnetic order. We produce a phase diagram in temperature and doping. We find that for the nickelate, like the cuprate, the stoichiometric compound is antiferromagnetic. We find superconductivity in a doping range bounded, with a small coexistence region, by the onset of antiferromagnetism at low doping and with transition temperature becoming immeasurably small at high doping. Superconductivity emerges at around the same hole doping for both compounds, but requires a larger deviation from half filling for the nickelate. Both antiferromagnetic and superconducting order lead to a partial gapping of the $d_{x^2-y^2}$ Fermi surface sheet. Our similar results for the cuprate and nickelate suggest that nickelate superconductivity is cupratelike. We compare our results to the experimental phase diagram.

preprint2022arXiv

Transient trapping into metastable states in systems with competing orders

The quench dynamics of a system involving two competing orders is investigated using a Ginzburg-Landau theory with relaxational dynamics. We consider the scenario where a pump rapidly heats the system to a high temperature, after which the system cools down to its equilibrium temperature. We study the evolution of the order parameter amplitude and fluctuations in the resulting time dependent free energy landscape. Exponentially growing thermal fluctuations dominate the dynamics. The system typically evolves into the phase associated with the faster-relaxing order parameter, even if it is not the global free energy minimum. This theory offers a natural explanation for the widespread experimental observation that metastable states may be induced by laser induced collapse of a dominant equilibrium order parameter.

preprint2021arXiv

Exploring nonequilibrium phases of photo-doped Mott insulators with Generalized Gibbs ensembles

Photo-excited strongly correlated systems can exhibit intriguing non-thermal phases, but the theoretical investigation of them poses significant challenges. In this work, we introduce a generalized Gibbs ensemble type description for long-lived photo-doped states in Mott insulators. This framework enables systematic studies of photo-induced phases based on equilibrium methods, as demonstrated here for the one-dimensional extended Hubbard model. We determine the nonequilibrium phase diagram, which features $η$-pairing and charge density wave phases in a wide doping range, and reveal physical properties of these phases. We show that the peculiar kinematics of photo-doped carriers, and the interaction between them, play an essential role in the formation of the non-thermal phases, and we clarify the differences between photo-doped Mott insulators, chemically-doped Mott insulators and photo-doped semiconductors. Our results demonstrate a new path for the systematic exploration of nonequilibrium strongly correlated systems and show that photo-doped Mott insulators host different phases than conventional semiconductors.

preprint2021arXiv

Moiré heterostructures as a condensed matter quantum simulator

Twisted van der Waals heterostructures have latterly received prominent attention for their many remarkable experimental properties, and the promise that they hold for realising elusive states of matter in the laboratory. We propose that these systems can, in fact, be used as a robust quantum simulation platform that enables the study of strongly correlated physics and topology in quantum materials. Among the features that make these materials a versatile toolbox are the tunability of their properties through readily accessible external parameters such as gating, straining, packing and twist angle; the feasibility to realize and control a large number of fundamental many-body quantum models relevant in the field of condensed-matter physics; and finally, the availability of experimental readout protocols that directly map their rich phase diagrams in and out of equilibrium. This general framework makes it possible to robustly realize and functionalize new phases of matter in a modular fashion, thus broadening the landscape of accessible physics and holding promise for future technological applications.

preprint2021arXiv

Strongly correlated ladders in K-doped $p$-terphenyl crystals

Potassium-doped terphenyl has recently attracted attention as a potential host for high-transition-temperature superconductivity. Here, we elucidate the many-body electronic structure of recently synthesized potassium-doped terphenyl crystals. We show that this system may be understood as a set of weakly coupled one-dimensional ladders. Depending on the strength of the inter-ladder coupling the system may exhibit spin-gapped valence-bond solid or antiferromagnetic phases, both of which upon hole doping may give rise to superconductivity. This terphenyl-based ladder material serves as a new platform for investigating the fate of ladder phases in presence of three-dimensional coupling as well as for novel superconductivity.

preprint2021arXiv

Topological charge pumping in excitonic insulators

We show that in excitonic insulators with $s$-wave electron-hole pairing, an applied electric field (either pulsed or static) can induce a $p$-wave component to the order parameter, and further drive it to rotate in the $s+ip$ plane, realizing a Thouless charge pump. In one dimension, each cycle of rotation pumps exactly two electrons across the sample. Higher dimensional systems can be viewed as a stack of one dimensional chains in momentum space in which each chain crossing the fermi surface contributes a channel of charge pumping. Physics beyond the adiabatic limit, including in particular dissipative effects is discussed.

preprint2021arXiv

Trigonal Symmetry Breaking and its Electronic Effects in Two-Dimensional Dihalides and Trihalides

We study the consequences of the approximately trigonal ($D_{3d}$) point symmetry of the transition metal (M) site in two-dimensional van der Waals MX$_2$ dihalides and MX$_3$ trihalides. The trigonal symmetry leads to a 2-2-1 orbital splitting of the transition metal $d$ shell, which may be tuned by the interlayer distance, and changes in the ligand-ligand bond lengths. Orbital order coupled to various lower symmetry lattice modes may lift the remaining orbital degeneracies, and we explain how these may support unique electronic states using ZrI$_2$ and CuCl$_2$ as examples, and offer a brief overview of possible electronic configurations in this class of materials. By building and analysing Wannier models adapted to the appropriate symmetry we examine how the interplay among trigonal symmetry, electronic correlation effects, and $p$-$d$ orbital charge transfer leads to insulating, orbitally polarized magnetic and/or orbital-selective Mott states. Our work establishes a rigorous framework to understand, control, and tune the electronic states in low-dimensional correlated halides. Our analysis shows that trigonal symmetry and its breaking is a key feature of the 2D halides that needs to be accounted for in search of novel electronic states in materials ranging from CrI$_3$ to $α$-RuCl$_3$.

preprint2020arXiv

Bardasis-Schrieffer polaritons in excitonic insulators

Bardasis-Schrieffer modes in superconductors are fluctuations in subdominant pairing channels, e.g., d-wave fluctuations in an s-wave superconductor. This Rapid Communication shows that these modes also generically occur in excitonic insulators. In s-wave excitonic insulators, a p-wave Bardasis-Schrieffer mode exists below the gap energy, is optically active and hybridizes strongly with photons to form Bardasis-Schrieffer polaritons, which are observable in both far-field and near-field optical experiments.

preprint2020arXiv

Charge stiffness and long-range correlation in the optically induced $η$-pairing state of the one-dimensional Hubbard model

We show that optical excitation of the Mott insulating phase of the one-dimensional Hubbard model can create a state possessing two of the hallmarks of superconductivity: a nonvanishing charge stiffness and long-ranged pairing correlation. By employing the exact diagonalization method, we find that the superposition of the $η$-pairing eigenstates induced by the optical pump exhibits a nonvanishing charge stiffness and a pairing correlation that decays very slowly with system size in sharp contrast to the behavior of an ensemble of thermally excited eigenstates, which has a vanishing charge stiffness and no long-ranged pairing correlations. We show that the charge stiffness is indeed directly associated with the $η$-pairing correlation in the Hubbard model. Our finding demonstrates that optical pumping can actually lead to superconducting-like properties on the basis of the $η$-pairing states.

preprint2020arXiv

Collective modes and terahertz near-field response of superconductors

We theoretically study the low energy electromagnetic response of BCS type superconductors focusing on propagating collective modes that are observable with THz near field optics. The interesting frequency and momentum range is $ω< 2Δ$ and $q < 1/ξ$ where $Δ$ is the gap and $ξ$ is the coherence length. We show that it is possible to observe the superfluid plasmons, amplitude (Higgs) modes, Bardasis-Schrieffer modes and Carlson-Goldman modes using THz near field technique, although none of these modes couple linearly to far field radiation. Coupling of THz near field radiation to the amplitude mode requires particle-hole symmetry breaking while coupling to the Bardasis-Schrieffer mode does not and is typically stronger. For parameters appropriate to layered superconductors of current interest, the Carlson-Goldman mode appears in the near field reflection coefficient as a weak feature in the sub-THz frequency range. In a system of two superconducting layers with nanometer scale separation, an acoustic phase mode appears as the antisymmetric density fluctuation mode of the system. This mode produces well defined resonance peaks in the near-field THz response and has strong anticrossings with the Bardasis-Schrieffer and amplitude modes, enhancing their response. In a slab consisting of many layers of quasi-two dimensional superconductors, realized for example in samples of high T$_c$ cuprate compounds, many branches of propagating Josephson plasmon modes are found to couple to the THz near field radiation.

preprint2020arXiv

Collective Modes in Excitonic Insulators: Effects of Electron-Phonon Coupling and Signatures in Optical Response

We consider a two-band spinless model describing an excitonic insulator (EI) on the two-dimensional square lattice with anisotropic hopping parameters and electron-phonon (el-ph) coupling, inspired by the EI candidate Ta$_2$NiSe$_5$. We systematically study the nature of the collective excitations in the ordered phase which originates from the interband Coulomb interaction and the el-ph coupling. When the ordered phase is stabilized only by the Coulomb interaction (pure EI phase), its collective response exhibits a massless phase mode in addition to the amplitude mode. We show that in the BEC regime, the signal of the amplitude mode becomes less prominent and that the anisotropy in the phase mode velocities is relaxed compared to the model bandstructure. Through coupling to the lattice, the phase mode acquires a mass and the signal of the amplitude mode becomes less prominent. Importantly, character of the softening mode at the boundary between the normal semiconductor phase and the ordered phase depends on the parameter condition. In particular, we point out that even for el-ph coupling smaller than the Coulomb interaction the mode that softens to zero at the boundary can have a phonon character. We also discuss how the collective modes can be observed in the optical conductivity. Furthermore, we study the effects of nonlocal interactions on the collective modes and show the possibility of realizing a coexistence of an in-gap mode and an above-gap mode split off from the single amplitude mode in the system with the local interaction only.

preprint2020arXiv

Ground-state properties of the hydrogen chain: insulator-to-metal transition, dimerization, and magnetic phases

Accurate and predictive computations of the quantum-mechanical behavior of many interacting electrons in realistic atomic environments are critical for the theoretical design of materials with desired properties, and require solving the grand-challenge problem of the many-electron Schrodinger equation. An infinite chain of equispaced hydrogen atoms is perhaps the simplest realistic model for a bulk material, embodying several central themes of modern condensed matter physics and chemistry, while retaining a connection to the paradigmatic Hubbard model. Here we report a combined application of cutting-edge computational methods to determine the properties of the hydrogen chain in its quantum-mechanical ground state. Varying the separation between the nuclei leads to a rich phase diagram, including a Mott phase with quasi long-range antiferromagnetic order, electron density dimerization with power-law correlations, an insulator-to-metal transition and an intricate set of intertwined magnetic orders.

preprint2020arXiv

Nature of symmetry breaking at the excitonic insulator transition: Ta$_2$NiSe$_5$

Ta$_2$NiSe$_5$ is one of the most promising materials for hosting an excitonic insulator ground state. While a number of experimental observations have been interpreted in this way, the precise nature of the symmetry breaking occurring in Ta$_2$NiSe$_5$, the electronic order parameter, and a realistic microscopic description of the transition mechanism are, however, missing. By a symmetry analysis based on first-principles calculations, we uncover the \emph{discrete} lattice symmetries which are broken at the transition. We identify a purely electronic order parameter of excitonic nature that breaks these discrete crystal symmetries and contributes to the experimentally observed lattice distortion from an orthorombic to a monoclinic phase. Our results provide a theoretical framework to understand and analyze the excitonic transition in Ta$_2$NiSe$_5$ and settle the fundamental questions about symmetry breaking governing the spontaneous formation of excitonic insulating phases in solid-state materials.

preprint2020arXiv

Negative absolute conductivity in photoexcited metals

We show that in a model of a metal photoexcited by a transient pump pulse resonant with a phonon mode the absolute dc conductivity may become negative, depending on the interplay between the electronic structure, the phonon frequency and the pump intensity. The analysis includes the effects of inelastic scattering and thermal relaxation. Results for the time evolution of the negative conductivity state are presented; the associated non-equilibrium physics may persist for long times after the pulse. Our findings provide a theoretical justification for previously proposed phenomenology and indicate new routes to the generation and exploration of intrinsically non-equilibrium states.

preprint2020arXiv

Nonlinear spectroscopy of collective modes in excitonic insulator

The nonlinear optical response of an excitonic insulator coupled to lattice degrees of freedom is shown to depend in strong and characteristic ways on whether the insulating behavior originates primarily from electron-electron or electron-lattice interactions. Linear response optical signatures of the massive phase mode and the amplitude (Higgs) mode are identified. Upon nonlinear excitation resonant to the phase mode, a new in-gap mode at twice the phase mode frequency is induced, leading to a huge second harmonic response. Excitation of in-gap phonon modes leads to different and much smaller effects. A Landau-Ginzburg theory analysis explain these different behavior and reveals that a parametric resonance of the strongly excited phase mode is the origin of the photo-induced mode in the electron-dominant case. The difference in the nonlinear optical response serve as a measure of the dominant mechanism of the ordered phase.

preprint2020arXiv

Pump induced motion of an interface between competing orders

We study the motion of an interface separating two regions with different electronic orders following a short duration pump that drives the system out of equilibrium. Using a generalized Ginzburg-Landau approach and assuming that the main effect of the nonequilibrium drive is to transiently heat the system we address the question of the direction of interface motion; in other words, which ordered region expands and which contracts after the pump. Our analysis includes the effects of differences in free energy landscape and in order parameter dynamics and identifies circumstances in which the drive may act to increase the volume associated with the subdominant order, for example when the subdominant order has a second order free energy landscape while the dominant order has a first order one.

preprint2020arXiv

The spontaneous symmetry breaking in Ta$_2$NiSe$_5$ is structural in nature

The excitonic insulator is an electronically-driven phase of matter that emerges upon the spontaneous formation and Bose condensation of excitons. Detecting this exotic order in candidate materials is a subject of paramount importance, as the size of the excitonic gap in the band structure establishes the potential of this collective state for superfluid energy transport. However, the identification of this phase in real solids is hindered by the coexistence of a structural order parameter with the same symmetry as the excitonic order. Only a few materials are currently believed to host a dominant excitonic phase, Ta$_2$NiSe$_5$ being the most promising. Here, we test this scenario by using an ultrashort laser pulse to quench the broken-symmetry phase of this transition metal chalcogenide. Tracking the dynamics of the material's electronic and crystal structure after light excitation reveals surprising spectroscopic fingerprints that are only compatible with a primary order parameter of phononic nature. We rationalize our findings through state-of-the-art calculations, confirming that the structural order accounts for most of the electronic gap opening. Not only do our results uncover the long-sought mechanism driving the phase transition of Ta$_2$NiSe$_5$, but they also conclusively rule out any substantial excitonic character in this instability.

preprint2019arXiv

Metal-Insulator and Magnetic Phase Diagram of Ca$_2$RuO$_4$ from Auxiliary Field Quantum Monte Carlo and Dynamical Mean Field Theory

Layered perovskite ruthenium oxides exhibit a striking series of metal-insulator and magnetic-nonmagnetic phase transitions easily tuned by temperature, pressure, epitaxy, and nonlinear drive. In this work, we combine results from two complementary state of the art many-body methods, Auxiliary Field Quantum Monte Carlo and Dynamical Mean Field Theory, to determine the low-temperature phase diagram of Ca$_2$RuO$_4$. Both methods predict a low temperature, pressure-driven metal-insulator transition accompanied by a ferromagnetic-antiferromagnetic transition. The properties of the ferromagnetic state vary non-monotonically with pressure and are dominated by the ruthenium $d_{xy}$ orbital, while the properties of the antiferromagnetic state are dominated by the $d_{xz}$ and $d_{yz}$ orbitals. Differences of detail in the predictions of the two methods are analyzed. This work is theoretically important as it presents the first application of the Auxiliary Field Quantum Monte Carlo method to an orbitally-degenerate system with both Mott and Hunds physics, and provides an important comparison of the Dynamical Mean Field and Auxiliary Field Quantum Monte Carlo methods.

preprint2016arXiv

Antisite defects at oxide interfaces

We use \textit{ab initio} calculations to estimate formation energies of cation (transition metal) antisite defects at oxide interfaces and to understand the basic physical effects that drive or suppress the formation of these defects. Antisite defects are found to be favored in systems with substantial charge transfer across the interface, while Jahn-Teller distortions and itinerant ferromagnetism can prevent antisite defects and help stabilize atomically sharp interfaces. Our results enable identification of classes of systems that may be more and less susceptible to the formation of antisite defects and motivate experimental studies and further theoretical calculations to elucidate the local structure and stability of oxide interface systems.

preprint2016arXiv

Charge transfer and electron-phonon coupling in monolayer FeSe on Nb doped SrTiO$_3$

Monolayer films of FeSe grown on SrTiO$_3$ substrates are electron doped relative to bulk and exhibit a significantly higher superconducting transition temperatures. We present density functional calculations and a modified Schottky model incorporating the strong paraelectricity of SrTiO$_3$ which demonstrate that the doping may be due to charge transfer from SrTiO$_3$ impurity bands driven by work function mismatch. Physically relevant levels of Nb doping are shown to lead to doping of the FeSe compatible with observation. The coupling of electrons in FeSe to polar phonons in the depletion regime of the SrTiO$_3$ is calculated. A $λ\sim 0.4$ is found; the coupling to long-wavelength phonons is found to be dominant.

preprint2016arXiv

Hybrid quantum-classical approach to correlated materials

Recent improvements in control of quantum systems make it seem feasible to finally build a quantum computer within a decade. While it has been shown that such a quantum computer can in principle solve certain small electronic structure problems and idealized model Hamiltonians, the highly relevant problem of directly solving a complex correlated material appears to require a prohibitive amount of resources. Here, we show that by using a hybrid quantum-classical algorithm that incorporates the power of a small quantum computer into a framework of classical embedding algorithms, the electronic structure of complex correlated materials can be efficiently tackled using a quantum computer. In our approach, the quantum computer solves a small effective quantum impurity problem that is self-consistently determined via a feedback loop between the quantum and classical computation. Use of a quantum computer enables much larger and more accurate simulations than with any known classical algorithm, and will allow many open questions in quantum materials to be resolved once a small quantum computer with around one hundred logical qubits becomes available.

preprint2016arXiv

Influence of quantum confinement and strain on orbital polarization of four-layer LaNiO$_3$ superlattices: a DFT+DMFT study

Atomically precise superlattices involving transition metal oxides provide a unique opportunity to engineer correlated electron physics using strain (modulated by choice of substate) and quantum confinement (controlled by layer thickness). Here we use the combination of density functional theory and dynamical mean field theory (DFT+DMFT) to study Ni E$_g$ $d$-orbital polarization in strained LaNiO$_3$/LaAlO$_3$ superlattices consisting of four layers of nominally metallic NiO$_2$ and four layers of insulating AlO$_2$ separated by LaO layers. The layer-resolved orbital polarization is calculated as a function of strain and analyzed in terms of structural, quantum confinement, and correlation effects. The effect of strain is determined from the dependence of the results on the Ni-O bond-length ratio and the octahedral rotation angles; quantum confinement is studied by comparison to bulk calculations with similar degrees of strain; correlation effects are inferred by varying interaction parameters within our DFT+DMFT calculations. The calculated dependence of orbital polarization on strain in superlattices is qualitatively consistent with recent X-ray absorption spectroscopy and resonant reflectometry data. However, interesting differences of detail are found between theory and experiment. Under tensile strain, the two inequivalent Ni ions display orbital polarization similar to that calculated for strained bulk LaNiO$_3$ and observed in experiment. Compressive strain produces a larger dependence of orbital polarization on Ni position and even the inner Ni layer exhibits orbital polarization different from that calculated for strained bulk LaNiO$_3$.

preprint2016arXiv

Phase diagram of Sr$_{1-x}$Ba$_x$MnO$_3$ as a function of chemical doping, epitaxial strain and external pressure

We use \textit{ab initio} calculations to systematically study the phase diagram of multiferroic Sr$_{1-x}$Ba$_x$MnO$_3$ ($0 \leq x \leq 1$) as a function of chemical doping, epitaxial strain and external pressure. We find that by replacing Sr with Ba in cubic SrMnO$_3$ and imposing epitaxial strain, the material can be tuned to the vicinity of a first order transition between two multiferroic phases, one antiferromagnetic with a smaller polarization and one ferromagnetic with a larger polarization. A giant effective magneto-electric coupling and cross-field control (electric field control of magnetism or magnetic field control of polarization) can be achieved in the vicinity of the transition. The dependence of the theoretically computed transition point on the choice of exchange correlation functionals is determined and is found to be non-negligible. We also show that the perovskite structure of BaMnO$_3$ can be stabilized relative to its hexagonal polymorphs at pressures larger than 20 GPa.

preprint2016arXiv

Projected Regression Methods for Inverting Fredholm Integrals: Formalism and Application to Analytical Continuation

We present a machine learning approach to the inversion of Fredholm integrals of the first kind. The approach provides a natural regularization in cases where the inverse of the Fredholm kernel is ill-conditioned. It also provides an efficient and stable treatment of constraints. The key observation is that the stability of the forward problem permits the construction of a large database of outputs for physically meaningful inputs. We apply machine learning to this database to generate a regression function of controlled complexity, which returns approximate solutions for previously unseen inputs; the approximate solutions are then projected onto the subspace of functions satisfying relevant constraints. We also derive and present uncertainty estimates. We illustrate the approach by applying it to the analytical continuation problem of quantum many-body physics, which involves reconstructing the frequency dependence of physical excitation spectra from data obtained at specific points in the complex frequency plane. Under standard error metrics the method performs as well or better than the Maximum Entropy method for low input noise and is substantially more robust to increased input noise. We expect the methodology to be similarly effective for any problem involving a formally ill-conditioned inversion, provided that the forward problem can be efficiently solved.

preprint2016arXiv

The Anderson--Holstein Model in Two Flavors of the Non--Crossing Approximation

The dynamical interplay between electron-electron interactions and electron-phonon coupling is investigated within the Anderson-Holstein model, a minimal model for open quantum systems that embody these effects. The influence of phonons on spectral and transport properties is explored in equilibrium, for non-equilibrium steady state and for transient dynamics after a quench. Both the particle-hole symmetric and the more generic particle-hole asymmetric cases are studied. The treatment is based on two complementary non-crossing approximations, the first of which is constructed around the weak-coupling limit and the second around the polaron limit. In general, the two methods disagree in nontrivial ways, indicating that more reliable approaches to the problem are needed. The frameworks used here can form the starting point for numerically exact methods based on bold-line continuous-time quantum Monte Carlo algorithms capable of treating open systems simultaneously coupled to multiple fermionic and bosonic baths.

preprint2015arXiv

Band Structure and Terahertz Optical Conductivity of Transition Metal Oxides: Theory and Application to CaRuO$_3$

Density functional plus dynamical mean field calculations are used to show that in transition metal oxides, rotational and tilting (GdFeO$_3$-type) distortions of the ideal cubic perovskite structure produce a multiplicity of low-energy optical transitions which affect the conductivity down to frequencies of the order of $1$ or $2$~mV (terahertz regime), mimicking non-Fermi-liquid effects even in systems with a strictly Fermi-liquid self-energy. For CaRuO$_3$, a material whose measured electromagnetic response in the terahertz frequency regime has been interpreted as evidence for non-Fermi-liquid physics, the combination of these band structure effects and a renormalized Fermi-liquid self-energy accounts for the low frequency optical response which had previously been regarded as a signature of exotic physics. Signatures of deviations from Fermi-liquid behavior at higher frequencies ($\sim 100$~meV) are discussed.

preprint2015arXiv

Charge transfer across transition metal oxide interfaces: emergent conductance and new electronic structure

We perform density functional theory plus dynamical mean field theory calculations to investi- gate internal charge transfer in an artificial superlattice composed of alternating layers of vanadate and manganite perovskite and Ruddlesden-Popper structure materials. We show that the elec- tronegativity difference between vanadium and manganese causes moderate charge transfer from VO2 to MnO2 layers in both perovskite and Ruddlesden-Popper based superlattices, leading to hole doping of the VO2 layer and electron doping of the MnO2 layer. Comparison of the perovskite and Ruddlesden-Popper based heterostructures provides insights into the role of the apical oxy- gen. Our first principles simulations demonstrate that the combination of internal charge transfer and quantum confinement provided by heterostructuring is a powerful approach to engineering electronic structure and tailoring correlation effects in transition metal oxides.

preprint2015arXiv

Density functional versus spin-density functional and the choice of correlated subspace in multi-variable effective action theories of electronic structure

Modern extensions of density functional theory such as the density functional theory plus U and the density functional theory plus dynamical mean-field theory require choices, including selection of variable (charge vs spin density) for the density functional and specification of the correlated subspace. This paper examines these issues in the context of the "plus U" extensions of density functional theory, in which additional correlations on specified correlated orbitals are treated using a Hartree-Fock approximation. Differences between using charge-only or spin-density-dependent exchange-correlation functionals and between Wannier and projector-based definitions of the correlated orbitals are considered on the formal level and in the context of the structural energetics of the rare earth nickelates. It is demonstrated that theories based on spin-dependent exchange-correlation functionals can lead to large and in some cases unphysical effective on-site exchange couplings. Wannier and projector-based definitions of the correlated orbitals lead to similar behavior near ambient pressure, but substantial differences are observed at large pressures. Implications for other beyond density functional methods such as the combination of density functional and dynamical mean field theory are discussed.

preprint2015arXiv

Electronic correlations, magnetism and Hund's rule coupling in the ruthenium perovskites SrRuO$_3$ and CaRuO$_3$

A comparative density functional plus dynamical mean field theory study of the pseudocubic ruthenate materials CaRuO$_3$ and SrRuO$_3$ is presented. Phase diagrams are determined for both materials as a function of Hubbard repulsion $U$ and Hund's rule coupling $J$. Metallic and insulating phases are found, as are ferromagnetic and paramagnetic states. The locations of the relevant phase boundaries are determined. Based on the computed phase diagrams, Mott-dominated and Hund's dominated regimes of strong correlation are distinguished. Comparison of calculated properties to experiments indicates that the actual materials are in the Hund's coupling dominated region of the phase diagram so can be characterized as Hund's metals, in common with other members of the ruthenate family. Comparison of the phase diagrams for the two materials reveals the role played by rotational and tilt (GdFeO$_3$-type) distortions of the ideal perovskite structure. The presence of magnetism in SrRuO$_3$ and its absence in CaRuO$_3$ despite the larger mass and larger tilt/rotational distortion amplitude of CaRuO$_3$ can be understood in terms of density of states effects in the presence of strong Hund's coupling. Comparison of the calculated low-$T$ properties of CaRuO$_3$ to those of SrRuO$_3$ provides insight into the effects of magnetic order on the properties of a Hund's metal. The study provides a simultaneous description of magnetism and correlations and explicates the roles played by band theory and Hubbard and Hund's interactions.

preprint2015arXiv

Imaginary-time matrix product state impurity solver for dynamical mean-field theory

We present a new impurity solver for dynamical mean-field theory based on imaginary-time evolution of matrix product states. This converges the self-consistency loop on the imaginary-frequency axis and obtains real-frequency information in a final real-time evolution. Relative to computations on the real-frequency axis, required bath sizes are much smaller and less entanglement is generated, so much larger systems can be studied. The power of the method is demonstrated by solutions of a three band model in the single and two-site dynamical mean-field approximation. Technical issues are discussed, including details of the method, efficiency as compared to other matrix product state based impurity solvers, bath construction and its relation to real-frequency computations and the analytic continuation problem of quantum Monte Carlo, the choice of basis in dynamical cluster approximation, and perspectives for off-diagonal hybridization functions.

preprint2015arXiv

Machine learning for many-body physics: efficient solution of dynamical mean-field theory

Machine learning methods for solving the equations of dynamical mean-field theory are developed. The method is demonstrated on the three dimensional Hubbard model. The key technical issues are defining a mapping of an input function to an output function, and distinguishing metallic from insulating solutions. Both metallic and Mott insulator solutions can be predicted. The validity of the machine learning scheme is assessed by comparing predictions of full correlation functions, of quasi-particle weight and particle density to values directly computed. The results indicate that with modest further development, machine learning approach may be an attractive computational efficient option for real materials predictions for strongly correlated systems.

preprint2015arXiv

Metal-insulator transitions in GdTiO3/SrTiO3 superlattices

The density functional plus U method is used to obtain the electronic structure, lattice relaxation and metal-insulator phase diagram of superlattices consisting of $m$ layers of Gadolinium Titanate (GdTiO$_{3}$) alternating with $n$ layers of Strontium Titanate (SrTiO$_{3}$). Metallic phases occur when the number of SrTiO$_3$ layers is large or the interaction $U$ is small. In metallic phases, the mobile electrons are found in the SrTiO$_3$ layers, with near-interface electrons occupying $xy$-derived bands, while away from the interface the majority of electrons reside in $xz/yz$ bands. As the thickness of the SrTiO$_3$ layers decreases or the on-site interaction U increases a metal-insulator transition occurs. Two different insulating states are found. When the number of SrTiO$_{3}$ layers is larger than one, we find an insulating state with two sublattice charge and orbital ordering and associated Ti-O bond length disproportionations. When the number of SrTiO$_{3}$ units per layer is one, a different insulating phase occurs driven by orbital ordering within the quasi one-dimensional $xz/yz$ bonding bands connecting Ti atoms across the SrO layer. In this phase there is no sublattice charge ordering or bond disproportionation. The critical U for the single-layer insulator is $\sim$ 2.5 eV, much less than critical U $\sim$ 3.5 eV required to drive the metal-insulator transition when the number of SrTiO$_3$ is larger than one. Inconsistencies between the calculation and the experiment suggest that many-body correlations may be important. A local inversion symmetry breaking around Ti atoms suggests the possibility of in-plane ferroelectric polarization in the insulating phase.

preprint2015arXiv

Photoinduced phase transitions in narrow-gap Mott insulators: the case of VO$_2$

We study the nonequilibrium dynamics of photoexcited electrons in the narrow-gap Mott insulator VO$_2$. The initial stages of relaxation are treated using a quantum Boltzmann equation methodology, which reveals a rapid ($\sim$ femtosecond time scale) relaxation to a pseudothermal state characterized by a few parameters that vary slowly in time. The long-time limit is then studied by a Hartree-Fock methodology, which reveals the possibility of nonequilibrium excitation to a new metastable $M_1$ metal phase that is qualitatively consistent with a recent experiment. The general physical picture of photoexcitation driving a correlated electron system to a new state that is not accessible in equilibrium may be applicable in similar materials.

preprint2015arXiv

Quasiparticle Properties of the Superconducting State of the Two Dimensional Hubbard Model

Cluster dynamical mean field methods are used to calculate the normal and anomalous components of the electron self energy of the two dimensional Hubbard model. From these the evolution of the superconducting gap and the momentum dependent photoemission and inverse photoemission spectra across the phase diagram are determined. In the pseudogap regime, decreasing the temperature into the superconducting state leads to a decrease in the energy gap and the formation of a `peak-dip-hump' structure in the electronic density of states. The peak feature disperses very weakly. The calculated spectral functions are in good qualitative agreement with published data. The mathematical origin of the behavior is found to be the effect of the superconductivity on the pole structure giving rise to the normal state pseudogap. In particular the "hump" feature is found to arise from a zero crossing of the real part of the electron self energy rather than from an onset of scattering. The effect of superconductivity on the zone diagonal spectra is presented.

preprint2015arXiv

Solutions of the Two Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms

Numerical results for ground state and excited state properties (energies, double occupancies, and Matsubara-axis self energies) of the single-orbital Hubbard model on a two-dimensional square lattice are presented, in order to provide an assessment of our ability to compute accurate results in the thermodynamic limit. Many methods are employed, including auxiliary field quantum Monte Carlo, bare and bold-line diagrammatic Monte Carlo, method of dual fermions, density matrix embedding theory, density matrix renormalization group, dynamical cluster approximation, diffusion Monte Carlo within a fixed node approximation, unrestricted coupled cluster theory, and multi-reference projected Hartree-Fock. Comparison of results obtained by different methods allows for the identification of uncertainties and systematic errors. The importance of extrapolation to converged thermodynamic limit values is emphasized. Cases where agreement between different methods is obtained establish benchmark results that may be useful in the validation of new approaches and the improvement of existing methods.

preprint2015arXiv

Spatial Correlations and the Insulating Phase of the High-$T_c$ Cuprates: Insights from a Configuration-Interaction-Based Solver for Dynamical Mean Field Theory

A recently proposed configuration-interaction based impurity solver is used in combination with the single-site and four-site cluster dynamical mean field approximations to investigate the three-band copper oxide model believed to describe the electronic structure of high transition temperature copper-oxide superconductors. Use of the configuration interaction solver enables verification of the convergence of results with respect to the number of bath orbitals. The spatial correlations included in the cluster approximation substantially shift the metal-insulator phase boundary relative to the prediction of the single-site approximation and increase the predicted energy gap of the insulating phase by about 1 eV above the single-site result. Vertex corrections occurring in the four-site approximation act to dramatically increase the value of the optical conductivity near the gap edge, resulting in a better agreement with the data. The calculations reveal two distinct correlated insulating states: the `magnetically correlated insulator', in which nontrivial intersite correlations play an essential role in stabilizing the insulating state, and the strongly correlated insulator, in which local physics suffices. Comparison of the calculations to the data place the cuprates in the magnetically correlated Mott insulator regime.

preprint2015arXiv

Spin-density functional theories and their $+U$ and $+J$ extensions: a comparative study of transition metals and transition metal oxides

Previous work on the physical content of exchange correlation functionals that depend on both charge and spin densities is extended to elemental transition metals and a wider range of perovskite transition metal oxides. A comparison of spectra and magnetic moments calculated using exchange correlation functionals depending on charge density only or on both charge and spin densities, as well as the $+U$ and $+J$ extensions of these methods confirms previous conclusions that the spin-dependent part of the exchange correlation functional provides an effective Hund's interaction acting on the transition metal $d$ orbitals. For the local spin density approximation and spin-dependent Perdew-Burke-Ernzerhof generalized gradient approximation, the effective Hund's exchange is found to be larger than 1 eV. The results indicate that at least as far as applications to transition metals and their oxides are concerned, $+U$, $+J$ and +dynamical mean field theory extensions of density functional theory should be based on exchange-correlation functionals of charge density only.

preprint2015arXiv

Strain Control of Electronic Phase in Rare Earth Nickelates

We use density functional plus $U$ methods to study the effects of a tensile or compressive substrate strain on the charge-ordered insulating phase of LuNiO$_3$. The numerical results are analyzed in terms of a Landau energy function, with octahedral rotational distortions of the perovskite structure included as a perturbation. Approximately 4% tensile or compressive strain leads to a first-order transition from an insulating structure with large amplitude breathing mode distortions of the NiO$_6$ octahedra to a metallic state in which breathing mode distortions are absent but Jahn-Teller distortions in which two Ni-O bonds become long and the other four become short are present. Compressive strain produces uniform Jahn-Teller order with the long axis aligned perpendicular to the substrate plane while tensile strain produces a staggered Jahn-Teller order in which the long bond lies in the plane and alternates between two nearly orthogonal in-plane directions forming a checkerboard pattern. In the absence of the breathing mode distortions and octahedral rotations, the tensile strain-induced transition to the staggered Jahn-Teller state would be of second order.

preprint2015arXiv

Taming the dynamical sign problem in real-time evolution of quantum many-body problems

Current nonequilibrium Monte Carlo methods suffer from a dynamical sign problem that makes simulating real-time dynamics for long times exponentially hard. We propose a new `Inchworm Algorithm', based on iteratively reusing information obtained in previous steps to extend the propagation to longer times. The algorithm largely overcomes the dynamical sign problem, changing the scaling from exponential to quadratic. We use the method to solve the Anderson impurity model in the Kondo and mixed valence regimes, obtaining results both for quenches and for spin dynamics in the presence of an oscillatory magnetic field.

preprint2014arXiv

Computing total energies in complex materials using charge self-consistent DFT+DMFT

We have formulated and implemented a fully charge-self-consistent density functional theory plus dynamical mean field theory methodology which enables an efficient calculation of the total energy of realistic correlated electron systems. The density functional portion of the calculation uses a plane wave basis set within the projector augmented wave method enabling study of systems with large, complex unit cells. The dynamical mean field portion of the calculation is formulated using maximally localized Wannier functions, enabling a convenient implementation which is independent of the basis set used in the density functional portion of the calculation. The importance of using a correct double counting term is demonstrated. A generalized form of the standard double counting correction, which we refer to as the $U^\prime$ form, is described in detail and used. For comparison the density functional plus U method is implemented within the same framework including the generalized double counting. The formalism is validated via a calculation of the metal-insulator and structural phase diagrams of the rare-earth nickelate perovskites as functions of applied pressure and A-site rare-earth ions. The calculated density functional plus dynamical mean field results are found to be consistent with experiment. The density functional plus U method is shown to grossly overestimate the tendency for bond-disproportionation and insulating behavior.

preprint2014arXiv

Covalency and the metal-insulator transition in titanate and vanadate perovskites

A combination of density functional and dynamical mean-field theory is applied to the perovskites SrVO$_3$, LaTiO$_3$ and LaVO$_3$. We show that DFT+DMFT in conjunction with the standard fully localized-limit (FLL) double-counting predicts that LaTiO$_3$ and LaVO$_3$ are metals even though experimentally they are correlation-driven ("Mott") insulators. In addition, the FLL double counting implies a splitting between oxygen $p$ and transition metal $d$ levels which differs from experiment. Introducing into the theory an \textit{ad hoc} double counting correction which reproduces the experimentally measured insulating gap leads also to a $p$-$d$ splitting consistent with experiment if the on-site interaction $U$ is chosen in a relatively narrow range ($\sim 6\pm 1$ eV). The results indicate that these early transition metal oxides will serve as critical test for the formulation of a general \textit{ab initio} theory of correlated electron metals.

preprint2014arXiv

Density functional plus dynamical mean field theory of the metal-insulator transition in early transition metal oxides

The combination of density functional theory and single-site dynamical mean-field theory, using both Hartree and full continuous-time quantum Monte Carlo impurity solvers, is used to study the metal-insulator phase diagram of perovskite transition-metal oxides of the form $AB$O$_3$ with a rare-earth ion $A$=Sr, La, Y and transition metal $B$=Ti, V, Cr. The correlated subspace is constructed from atomiclike $d$ orbitals defined using maximally localized Wannier functions derived from the full $p$-$d$ manifold; for comparison, results obtained using a projector method are also given. Paramagnetic DFT+DMFT computations using full charge self-consistency along with the standard "fully localized limit" (FLL) double counting are shown to incorrectly predict that LaTiO$_3$, YTiO$_3$, LaVO$_3$ and SrMnO$_3$ are metals. A more general examination of the dependence of physical properties on the mean $p$-$d$ energy splitting, the occupancy of the correlated $d$ states, the double-counting correction, and the lattice structure demonstrates the importance of charge-transfer physics even in the early transition-metal oxides and elucidates the factors underlying the failure of the standard approximations. If the double counting is chosen to produce a $p$-$d$ splitting consistent with experimental spectra, single-site dynamical mean-field theory provides a reasonable account of the materials properties. The relation of the results to those obtained from "$d$-only" models in which the correlation problem is based on the frontier orbital $p$-$d$ antibonding bands is determined. It is found that if an effective interaction $U$ is properly chosen the $d$-only model provides a good account of the physics of the $d^1$ and $d^2$ materials.

preprint2014arXiv

Determination of Intrinsic Magnetic Response from Local Measurements of Fringing Fields

Micron-sized Hall bars and micro-SQUIDs are now used routinely to measure the local static and dynamic magnetic response with micron-scale spatial resolution. While this provides a powerful new tool, determining the intrinsic magnetization presents new challenges, as it requires correcting for demagnetization fields that vary widely with position on a sample. In this paper we develop a method to correct for the demagnetization effect at local points of a rectangular prism shaped sample using a finite element analysis of Maxwell's equation applied to local Hall sensor measurements calibrated by bulk measurements of the magnetization. This method can be generalized to other geometric shapes to analyze data obtained with local magnetic probes.

preprint2014arXiv

Effect of dilute, strongly pinning impurities on charge density waves

We study theoretically the effects of strong pinning centers on a charge density wave in the limit that the charge density wave coherence length is shorter than the average inter-impurity distance. An analysis based on a Ginzburg-Landau model shows that long range forces arising from the elastic response of the charge density wave induce a kind of collective pinning which suppresses impurity-induced phase fluctuations leading to a long ranged ordered ground state. The effective correlations among impurities are characterized by a length scale parametrically longer than the average inter-impurity distance. The thermal excitations are found to be gapped implying the stability of the ground state. We also present Monte Carlo simulations that confirm the basic features of the analytical results.

preprint2014arXiv

Emergent properties hidden in plane view: Strong electronic correlations at oxide interfaces

Finding new collective electronic states in materials is one of the fundamental goals of condensed matter physics. Atomic-scale superlattices formed from transition metal oxides are a particularly appealing hunting ground for new physics. In bulk form, transition metal oxides exhibit a remarkable range of magnetic, superconducting, and multiferroic phases that are of great scientific interest and are potentially capable of providing innovative energy, security, electronics and medical technology platforms. In superlattices new states may emerge at the interfaces where dissimilar materials meet. Here we illustrate the essential features that make transition metal oxide-based heterostructures an appealing discovery platform for emergent properties with a few selected examples, showing how charge redistributes, magnetism and orbital polarization arises and ferroelectric order emerges from heterostructures comprised of oxide components with nominally contradictory behavior with the aim providing insight into the creation and control of novel behavior at oxide interfaces by suitable mechanical, electrical or optical boundary conditions and excitations.

preprint2014arXiv

Experimental evidence for a Bragg glass density wave phase in a transition-metal dichalcogenide

Analysis of the spatial dependence of current-voltage characteristics obtained from scanning tunneling microscopy experiments indicates that the charge density wave (CDW) occurring in NbSe$_2$ is subject to locally strong pinning by a non-negligible density of defects, but that on the length scales accessible in this experiment the material is in a "Bragg glass" phase where dislocations and anti-dislocations occur in bound pairs and free dislocations are not observed. A Landau theory-based analysis is presented showing how a strong local modulation may produce only a weak long range effect on the CDW phase.

preprint2014arXiv

Green's functions from real-time bold-line Monte Carlo

We present two methods for computing two-time correlation functions or Green's functions from real time bold-line continuous time quantum Monte Carlo. One method is a formally exact generalized auxiliary lead formalism by which spectral properties may be obtained from single-time observables. The other involves the evaluation of diagrams contributing to two-time observables directly on the Keldysh contour. Additionally, we provide a detailed description of the bold-line Monte Carlo method. Our methods are general and numerically exact, and able to reliably resolve high-energy features such as band edges. We compare the spectral functions obtained from real time methods to analytically continued spectral functions obtained from imaginary time Monte Carlo, thus probing the limits of analytic continuation.

preprint2014arXiv

Green's functions from real-time bold-line Monte Carlo: spectral properties of the nonequilibrium Anderson impurity model

The nonequilibrium spectral properties of the Anderson impurity model with a chemical potential bias are investigated within a numerically exact real time quantum Monte Carlo formalism. The two-time correlation function is computed in a form suitable for nonequilibrium dynamical mean field calculations. Additionally, the evolution of the model's spectral properties are simulated in an alternative representation, defined by a hypothetical but experimentally realizable weakly coupled auxiliary lead. The voltage splitting of the Kondo peak is confirmed and the dynamics of its formation after a coupling or gate quench are studied. This representation is shown to contain additional information about the dot's population dynamics. Further, we show that the voltage-dependent differential conductance gives a reasonable qualitative estimate of the equilibrium spectral function, but significant qualitative differences are found including incorrect trends and spurious temperature dependent effects.

preprint2014arXiv

Machine learning for many-body physics: The case of the Anderson impurity model

Machine learning methods are applied to finding the Green's function of the Anderson impurity model, a basic model system of quantum many-body condensed-matter physics. Different methods of parametrizing the Green's function are investigated; a representation in terms of Legendre polynomials is found to be superior due to its limited number of coefficients and its applicability to state of the art methods of solution. The dependence of the errors on the size of the training set is determined. The results indicate that a machine learning approach to dynamical mean-field theory may be feasible.

preprint2014arXiv

Pairing Glue in the Two Dimensional Hubbard Model

Cluster dynamical mean field calculations are used to construct the superconducting gap function of the two dimensional Hubbard model. The frequency dependence of the imaginary part of the gap function indicates that the pairing is dominated by fluctuations at two characteristic frequencies: one at the scale of the hopping matrix element $t$ and one at a much lower scale. The lower frequency component becomes more important as the doping is reduced into the pseudogap regime. Comparison to available information on the spin fluctuation spectrum of the model suggests that the superconductivity arises from exchange of spin fluctuations. The inferred pairing glue function is in remarkable qualitative consistency with the pairing function inferred from time-resolved optical conductivity data.

preprint2014arXiv

Selectively Localized Wannier Functions

Since the seminal work of Marzari and Vanderbilt, maximally localized Wannier functions have become widely used as a real-space representation of the electronic structure of periodic materials. In this paper we introduce selectively localized Wannier functions which extend the method of Marzari and Vanderbilt in two important ways. First, our method allows us to focus on localizing a subset of orbitals of interest. Second, our method allows us to fix centers of these orbitals, and ensure the preservation of the point-group symmetry. These characteristics are important when Wannier functions are used in methodologies that go beyond density functional theory by treating a local subspace of the Hamiltonian more effectively. Application of our method to GaAs, SrMnO$_3$, and Co demonstrates that selectively localized Wannier functions can offer improvements over the maximally localized Wannier function technique.

preprint2014arXiv

Three Dimensional Metallic and Two Dimensional Insulating Behavior in Tantalum Dichalcogenides

Using density functional theory with added on-site interactions (DFT+U), we study the electronic structure of bulk, monolayer, and bilayer of the layered transition-metal dichalcogenide $1T-TaS_2$. We show that a two-dimensional spin--$\frac{1}{2}$ Mott-phase exists for the monolayer in the charge-density wave state (CDW) and that such a phase is systematically destroyed by packing of the distorted layers leading to a one dimensional metal for bulk, CDW-distorted TaS$_2$. The latter finding is in contrast with previous DMFT predictions --disagreement we explain by the weak effective interaction felt by the electrons in the CDW state. Experimental observations of insulating behavior may arise from disorder due to stacking faults.

preprint2014arXiv

Total energy calculations using DFT+DMFT: computing the pressure phase diagram of the rare earth nickelates

A full implementation of the $ab$ $initio$ density functional plus dynamical mean field theory (DFT+DMFT) formalism to perform total energy calculations and structural relaxations is proposed and implemented. The method is applied to the structural and metal-insulator transitions of the rare earth nickelate perovskites as a function of rare earth ion, pressure, and temperature. In contrast to previous DFT and DFT+$U$ theories, the present method accounts for the experimentally observed structure of $La$NiO$_3$ and the insulating nature of the other perovskites, and quantitatively reproduces the metal-insulator and structural phase diagram in the plane of pressure and rare earth element. The temperature dependence of the energetics of the phase transformation indicates that the thermal transition is driven by phonon entropy effects.

preprint2013arXiv

Charge density distribution and optical response of the LaAlO3/SrTiO3 interface

We present calculations of the charge density profile, subband occupancy and ellipsometry spectra of the electron gas at the LaAlO3/SrTiO3 interface. The calculations employ self-consistent Hartree and random phase approximations, a tight binding parametrization of the band structure and a model for the optical phonon of SrTiO3. The dependence of the spatial structure and occupancy of subbands on the magnitude of the polarization charge at the interface and the dielectric function is determined. The interface-confined subbands may be labelled by the symmetry (xy, xz, or yz) of the Ti d-orbitals from which they mainly derive. The xy-derived band nearest the interface contains the major proportion of the electronic charge, but a large number of more distant, slightly occupied xy-derived bands are also found. Depending on the magnitude of the polarization charge, zero, one, or two xz/yz derived subbands are found. When present, these xz/yz bands give the dominant contribution to the long-distance tail of the interface charge. The response to applied ac electric fields polarized parallel and perpendicular to the interface is calculated and the results are presented in terms of ellipsometry angles. Two features are found: a dip in the spectrum near the LO feature of the STO phonon and a peak at the higher energy. We show that the form and magnitude of the dip is related to the Drude response of carriers moving in the plane of the interface while the peak arises from the plasmon excitation of the xz and yz electrons. The relation of the features of the subband occupancies and the in-plane conductivities is given.

preprint2013arXiv

Designing ferromagnetism in vanadium-oxide based superlattices

Motivated by recent reports (Phys. Rev. B\textbf{80}, 241102) of room-temperature ferromagnetism in vanadium-oxide based superlattices, a single-site dynamical mean field study of the dependence of the paramagnetic-ferromagnetic phase boundary on superlattice geometry was performed. An examination of variants of the experimentally determined crystal structure indicate that ferromagnetism is found only in a small and probably inaccessible region of the phase diagram. Design criteria for increasing the range over which ferromagnetism might exist are proposed.

preprint2013arXiv

Engineering Correlation Effects via Artificially Designed Oxide Superlattices

Ab initio calculations are used to predict that a superlattice composed of layers of LaTiO3 and LaNiO3 alternating along the [001] direction is a S=1 Mott insulator with large magnetic moments on the Ni sites, negligible moments on the Ti sites and a charge transfer gap set by the energy difference between Ni d and Ti d states, distinct from conventional Mott insulators. Correlation effects are enhanced on the Ni sites via filling the oxygen p states and reducing the Ni-O-Ni bond angle. Small hole (electron) doping of the superlattice leads to a two-dimensional single-band situation with holes (electrons) residing on Ni dx2-y2 (Ti dxy) orbital and coupled to antiferromagnetically correlated spins in the NiO2 layer.

preprint2013arXiv

Nematicity as a probe of superconducting pairing in iron-based superconductors

In several families of iron-based superconducting materials, a d-wave pairing instability may compete with the leading s-wave instability. Here we show that when both states have comparable free energies, superconducting and nematic degrees of freedom are strongly coupled. While nematic order causes a sharp non-analytic increase in $T_{c}$, nematic fluctuations can change the character of the s-wave to d-wave transition, favoring an intermediate state that does not break time-reversal symmetry but does break tetragonal symmetry. The coupling between superconductivity and nematicity is also manifested in the strong softening of the shear modulus across the superconducting transition. Our results show that nematicity can be used as a diagnostic tool to search for unconventional pairing states in iron pnictides and chalcogenides.

preprint2013arXiv

Sum rule violation in self-consistent hybridization expansions

We show that for multi-orbital quantum impurity models the non-crossing approximation and one-crossing approximation versions of the self-consistent hybridization expansions violate the sum rules relating the coefficients of the high-frequency expansion of the self energy and the product of the self energy and Green function to thermodynamic expectation values. Comparison of non-crossing/one-crossing results to numerically exact quantum Monte-Carlo calculations shows that the consistency with sum rules provides a useful estimate of the reliability of the approximations. The sum rule violations are more pronounced, and therefore the quality of the non-crossing/one-crossing approximation is poorer, in situations with multiple orbitals and away from particle-hole symmetry but becomes less severe as the correlation strength increases. The one crossing approximation is markedly superior to the non-crossing approximation.

preprint2013arXiv

Superconducting and Pseudogap effects on the interplane conductivity and Raman scattering cross section in the two dimensional Hubbard Model

Cluster dynamical mean field methods are used to calculate the superconductivity-induced changes in the interplane conductivity and Raman scattering cross section of the two dimensional Hubbard model. When superconductivity emerges from the pseudogap, the superconducting response is found to be diminished in amplitude, broadened and, in the case of the interplane conductivity, shifted to higher frequency. The results are in agreement with data on high temperature copper-oxide superconductors indicating that the Hubbard model contains the essential low energy physics of the pseudogap and its interplay with superconductivity in the cuprates.

preprint2013arXiv

Superconductivity and the Pseudogap in the two-dimensional Hubbard model

Recently developed numerical methods have enabled the explicit construction of the superconducting state of the Hubbard model of strongly correlated electrons in parameter regimes where the model also exhibits a pseudogap and a Mott insulating phase. $d_{x^2-y^2}$ symmetry superconductivity is found to occur in proximity to the Mott insulator, but separated from it by a pseudogapped nonsuperconducting phase. The superconducting transition temperature and order parameter amplitude are found to be maximal at the onset of the normal-state pseudogap. The emergence of superconductivity from the normal state pseudogap leads to a decrease in the excitation gap. All of these features are consistent with the observed behavior of the copper-oxide superconductors.

preprint2013arXiv

Suppression of superconductivity by Neel-type magnetic fluctuations in the iron pnictides

Motivated by recent experimental detection of Neel-type ($(π,π)$) magnetic fluctuations in some iron pnictides, we study the impact of competing $(π,π)$ and $(π,0)$ spin fluctuations on the superconductivity of these materials. We show that, counter-intuitively, even short-range, weak Neel fluctuations strongly suppress the $s^{+-}$ state, with the main effect arising from a repulsive contribution to the $s^{+-}$ pairing interaction, complemented by low frequency inelastic scattering. Further increasing the strength of the Neel fluctuations leads to a low-$T_{c}$ d-wave state, with a possible intermediate $s+id$ phase. The results suggest that the absence of superconductivity in a series of hole-doped pnictides is due to the combination of short-range Neel fluctuations and pair-breaking impurity scattering, and also that $T_{c}$ of optimally doped pnictides could be further increased if residual $(π,π)$ fluctuations were reduced.

preprint2013arXiv

The failure of DFT-based computations for a stepped-substrate-supported correlated Co wire

Density functional theory (DFT) has been immensely successful in its ability to predict physical properties, and, in particular, structures of condensed matter systems. Here, however, we show that DFT qualitatively fails to predict the dimerized structural phase for a monatomic Co wire that is self-assembled on a vicinal, i.e. stepped, Cu(111) substrate. To elucidate the nature of this failure, we compute the energetics of a Co chain on a Cu surface, step, notch, and embedded in bulk. The results demonstrate that increasing Co coordination extinguishes the dimerization, indicating that the failure of DFT for Co on the Cu step arises from excessive hybridization, which both weakens the ferromagnetic correlations that drive the dimerization and increases the bonding that opposes dimerization. Additionally, we show that including local interactions via DFT+U or DFT+DMFT does not restore the dimerization for the step-substrate supported wire, though the Co wire does dimerize in DFT+DMFT for the isolated vacuum case. This system can serve as a benchmark for future electronic structure methods.

preprint2013arXiv

Theory of ferromagnetism in vanadium-oxide based perovskites

The conditions under which ferromagnetism may occur in transition metal oxides with partially filled $t_{2g}$ shells such as vanadium-based perovskites are studied using a combination of density functional and single-site dynamical mean field methods. For reasonable values of the correlation strength, rotations of the VO$_6$ octahedra play an important role in enabling ferromagnetism, with ferromagnetism typically occurring for rotations larger than a nonzero critical value. Ferromagnetism is suppressed near the Mott insulating phase but the phase boundary is otherwise only weakly dependent on carrier concentration. Design rules are suggested for new oxide systems exhibiting ferromagnetism.

preprint2012arXiv

Covalency, double-counting and the metal-insulator phase diagram in transition metal oxides

Dynamical mean field theory calculations are used to show that for late transition-metal-oxides a critical variable for the Mott/charge-transfer transition is the number of d-electrons, which is determined by charge transfer from oxygen ions. Insulating behavior is found only for a narrow range of d-occupancy, irrespective of the size of the intra-d Coulomb repulsion. The result is useful in interpreting 'density functional +U' and 'density functional plus dynamical mean field' methods in which additional correlations are applied to a specific set of orbitals and an important role is played by the 'double counting correction' which dictates the occupancy of these correlated orbitals. General considerations are presented and are illustrated by calculations for two representative transition metal oxide systems: layered perovskite Cu-based "high-Tc" materials, an orbitally non-degenerate electronically quasi-two dimensional systems, and pseudocubic rare earch nickelates, an orbitally degenerate electronically three dimensional system. Density functional calculations yield d-occupancies very far from the Mott metal-insulator phase boundary in the nickelate materials, but closer to it in the cuprates, indicating the sensitivity of theoretical models of the cuprates to the choice of double counting correction and corroborating the critical role of lattice distortions in attaining the experimentally observed insulating phase in the nickelates.

preprint2012arXiv

Energetics of superconductivity in the two dimensional Hubbard model

The energetics of the interplay between superconductivity and the pseudogap in high temperature superconductivity is examined using the eight-site dynamical cluster approximation to the two dimensional Hubbard model. Two regimes of superconductivity are found: a weak coupling/large doping regime in which the onset of superconductivity causes a reduction in potential energy and an increase in kinetic energy, and a strong coupling regime in which superconductivity is associated with an increase in potential energy and decrease in kinetic energy. The crossover between the two regimes is found to coincide with the boundary of the normal state pseudogap, providing further evidence of the unconventional nature of superconductivity in the pseudogap regime. However the absence, in the strongly correlated but non-superconducting state, of discernibly nonlinear response to an applied pairing field, suggests that resonating valence bond physics is not the origin of the kinetic-energy driven superconductivity.

preprint2012arXiv

Numerically Exact Long Time Magnetization Dynamics at the Nonequilibrium Kondo Crossover of the Anderson Impurity Model

We investigate the dynamical and steady-state spin response of the nonequilibrium Anderson model to magnetic fields, bias voltage, and temperature using a numerically exact method combining a bold-line quantum Monte Carlo technique with the memory function formalism. We obtain converged results in a range of previously inaccessible regimes, in particular the crossover to the Kondo domain. We provide detailed predictions for novel nonequilibrium phenomena, including non-monotonic temperature dependence of observables at high bias voltage and oscillatory quench dynamics at high magnetic fields.

preprint2012arXiv

Site-selective Mott transition in rare earth nickelates

A combination of density functional and dynamical mean field theory calculations are used to show that the remarkable metal-insulator transition in the rare earth nickelate perovskites arise from a site-selective Mott phase, in which the $d$-electrons on a half of the Ni ions are localized to form a fluctuating moment while the $d$-electrons on other Ni ions form a singlet with holes on the surrounding oxygen ions. The calculation reproduces key features observed in the nickelate materials, including an insulating gap in the paramagnetic state, a strong variation of static magnetic moments among Ni sites and an absence of "charge order". A connection between structure and insulating behavior is documented. The site-selective Mott transition may be a more broadly applicable concept in the description of correlated materials.

preprint2012arXiv

Spin-exchange-induced dimerization of an atomic 1-D system

Using low-temperature scanning tunneling microscopy, we demonstrate an unambiguous 1-D system that surprisingly undergoes a CDW instability on a metallic substrate. Our ability to directly and quantitatively measure the structural distortion of this system provides an accurate reference for comparison with first principles theory. In comparison to previously proposed physical mechanisms, we attribute this particular 1-D CDW instability to a ferromagnetic state. We show that though the linear arrayed dimers are not electronically isolated, they are magnetically independent, and hence can potentially serve as a binary spin-memory system.

preprint2012arXiv

Theory of the magnetic and metal-insulator transitions in RNiO3 bulk and layered structures

A slave rotor-Hartree Fock formalism is presented for studying the properties of the p-d model describing perovskite transition metal oxides, and a flexible and efficient numerical formalism is developed for its solution. The methodology is shown to yield, within an unified formulation, the significant aspects of the rare earth nickelate phase diagram, including the paramagnetic metal state observed for the LaNiO3 and the correct ground-state magnetic order of insulating compounds. It is then used to elucidate ground state changes occurring as morphology is varied from bulk to strained and un-strained thin-film form. For ultrathin films, epitaxial strain and charge-transfer to the apical out-of-plane oxygen sites are shown to have significant impact on the phase diagram.

preprint2012arXiv

Two-particle response in Cluster Dynamical Mean-Field Theory: Formalism and application to the Raman Response of High-temperature Superconductors

A method is presented for the unbiased numerical computation of two-particle response functions of correlated electron materials via a solution of the dynamical mean-field equations in the presence of a perturbing field. The power of the method is demonstrated via a computation of the Raman $B_{1g}$ and $B_{2g}$ scattering intensities of the two dimensional Hubbard model, in parameter regimes believed to be relevant to high-temperature superconductivity. The theory reproduces the `two-magnon' peak characteristic of the Raman intensity of the insulating parent compounds of high-$T_c$ copper oxide superconductors and shows how it evolves to a quasiparticle response as carriers are added. The method can be applied in any situation where a solution of the equilibrium dynamical mean-field equations is feasible.

preprint2011arXiv

Continuous-time Monte Carlo methods for quantum impurity models

Quantum impurity models describe an atom or molecule embedded in a host material with which it can exchange electrons. They are basic to nanoscience as representations of quantum dots and molecular conductors and play an increasingly important role in the theory of "correlated electron" materials as auxiliary problems whose solution gives the "dynamical mean field" approximation to the self energy and local correlation functions. These applications require a method of solution which provides access to both high and low energy scales and is effective for wide classes of physically realistic models. The continuous-time quantum Monte Carlo algorithms reviewed in this article meet this challenge. We present derivations and descriptions of the algorithms in enough detail to allow other workers to write their own implementations, discuss the strengths and weaknesses of the methods, summarize the problems to which the new methods have been successfully applied and outline prospects for future applications.

preprint2011arXiv

Current driven defect unbinding transition in an XY ferromagnet

A Keldysh-contour effective field theory is derived for magnetic vortices in the presence of current flow. The effect of adiabatic and non-adiabatic spin transfer torques on vortex motion is highlighted. Similarities to and differences from the superconducting case are presented and explained. Current flow across a magnetically ordered state is shown to lead to a defect-unbinding phase transition which is intrinsically nonequilibrium in the sense of not being driven by a variation in effective temperature. The dependence of the density of vortices on the current density is determined.

preprint2011arXiv

d_{3z^2-r^2} Orbital in high-Tc cuprates: Excitonic spectrum, metal-insulator phase diagram, optical conductivity and orbital character of doped holes

The single-site dynamical mean-field approximation is used to solve a model of high-Tc cuprate superconductors which includes both d_{x^2-y^2} and d_{3z^2-r^2} orbitals on the Cu as well as the relevant oxygen states. Both T (with apical oxygen) and T' (without apical oxygen) crystal structures are considered. In both phases, inclusion of the d_{3z^2-r^2} orbital is found to broaden the range of stability of the charge transfer insulating phase. For equal charge transfer energies and interaction strengths, the T' phase is found to be less strongly correlated than the T phase. For both structures, d-d excitons are found within the charge-transfer gap. However, for all physically relevant dopings the Fermi surface is found to have only one sheet and the admixture of d_{3z^2-r^2} into ground state wave function remains negligible (<5%). Inclusion of the extra orbitals is found not to resolve the discrepancy between computed and observed conductivity in the insulating state.

preprint2011arXiv

Dynamical Mean-Field Theory for Quantum Chemistry

The dynamical mean-field concept of approximating an unsolvable many-body problem in terms of the solution of an auxiliary quantum impurity problem, introduced to study bulk materials with a continuous energy spectrum, is here extended to molecules, i.e., finite systems with a discrete energy spectrum. The application to small clusters of hydrogen atoms yields ground state energies which are competitive with leading quantum chemical approaches at intermediate and large interatomic distances as well as good approximations to the excitation spectrum.

preprint2011arXiv

High-frequency asymptotic behavior of self-energies in quantum impurity models

We present explicit expressions for the high-frequency asymptotic behavior of electron self-energy of general quantum impurity models, which may be useful for improving the convergence of dynamical mean-field calculations and for the analytic continuation of the electron self-energy. We also give results, expressed in more physical terms, for the two-orbital and three-orbital rotationally invariant Slater-Kanamori interactions, in order to facilitate calculations of transition metal oxides.

preprint2011arXiv

Numerically Exact Long Time Behavior of Nonequilibrium Quantum Impurity Models

A Monte Carlo sampling of diagrammatic corrections to the non-crossing approximation is shown to provide numerically exact estimates of the long-time dynamics and steady state properties of nonequilibrium quantum impurity models. This `bold' expansion converges uniformly in time and significantly ameliorates the sign problem that has heretofore limited the power of real-time Monte Carlo approaches to strongly interacting real-time quantum problems. The new approach enables the study of previously intractable problems ranging from generic long time nonequilibrium transport characteristics in systems with large onsite repulsion to the direct description of spectral functions on the real frequency axis in Dynamical Mean Field Theory.

preprint2011arXiv

Role of oxygen-oxygen hopping in the three-band copper-oxide model: quasiparticle weight, metal insulator and magnetic phase boundaries, gap values and optical conductivity

We investigate the effect of oxygen-oxygen hopping on the three-band copper-oxide model relevant to high-$T_c$ cuprates, finding that the physics is changed only slightly as the oxygen-oxygen hopping is varied. The location of the metal-insulator phase boundary in the plane of interaction strength and charge transfer energy shifts by $\sim 0.5$eV or less along the charge transfer axis, the quasiparticle weight has approximately the same magnitude and doping dependence and the qualitative characteristics of the electron-doped and hole-doped sides of the phase diagram do not change. The results confirm the identification of La$_2$CuO$_4$ as a material with intermediate correlation strength. However, the magnetic phase boundary as well as higher-energy features of the optical spectrum are found to depend on the magnitude of the oxygen-oxygen hopping. We compare our results to previously published one-band and three-band model calculations.

preprint2011arXiv

Satellites and large doping- and temperature-dependence of electronic properties in hole-doped BaFe2As2

Over the last years, superconductivity has been discovered in several families of iron-based compounds. Despite intense research, even basic electronic properties of these materials, such as Fermi surfaces, effective electron masses, or orbital characters are still subject to debate. Here, we address an issue that has not been considered before, namely the consequences of dynamical screening of the Coulomb interactions among Fe-d electrons. We demonstrate its importance not only for correlation satellites seen in photoemission spectroscopy, but also for the low-energy electronic structure. From our analysis of the normal phase of BaFe2As2 emerges the picture of a strongly correlated compound with strongly doping- and temperature-dependent properties. In the hole overdoped regime, an incoherent metal is found, while Fermi-liquid behavior is recovered in the undoped compound. At optimal doping, the self-energy exhibits an unusual square-root energy dependence which leads to strong band renormalizations near the Fermi level.

preprint2010arXiv

Bold Line Diagrammatic Monte Carlo Method: General formulation and application to expansion around the Non-Crossing Approximation

We present a general framework for performing "bold-line" diagrammatic Monte Carlo calculations based on using an analytical partial resummation as a starting point for stochastically summing all diagrams. As an example we solve the one orbital Anderson model using the non-crossing approximation as a starting point. We assess the accuracy of the method by using it to solve the equations of single-site dynamical mean field theory. We establish the validity of the starting approximations and show that the bold method provides a very accurate treatment of the Mott-insulating phase.

preprint2010arXiv

Chemical control of orbital polarization in artificially structured transition-metal oxides: La2NiXO6 (X=B, Al, Ga, In) from first principles

The application of modern layer-by-layer growth techniques to transition-metal oxide materials raises the possibility of creating new classes of materials with rationally designed correlated electron properties. An important step toward this goal is the demonstration that electronic structure can be controlled by atomic composition. In compounds with partially occupied transition-metal d shells, one important aspect of the electronic structure is the relative occupancy of different d orbitals. Previous work has established that strain and quantum confinement can be used to influence orbital occupancy. In this paper we demonstrate a different modality for orbital control in transition-metal oxide heterostructures, using density-functional band calculations supplemented by a tight-binding analysis to show that the choice of nontransition-metal counterion X in transition-metal oxide heterostructures composed of alternating LaNiO3 and LaXO3 units strongly affects orbital occupancy, changing the magnitude and in some cases the sign of the orbital polarization.

preprint2010arXiv

Dynamical Screening in Correlated Electron Materials

We present an efficient method for incorporating the dynamical effects of the screening of the Hubbard U by electronic degrees of freedom in the solid into the single site dynamical mean field approximation. The formalism is illustrated by model system calculations which capture the essential features of the frequency dependent interactions proposed for Gd, Ni, SrVO_3 and other compounds. Screening leads to shifts in the metal-insulator phase boundary, changes in the spectral function near the Mott-Hubbard gap edge and to a renormalization of the quasiparticle weight. Hubbard bands are generically neither separated by the screened nor the unscreened interaction energy, implying that the common practice of extracting the Hubbard U from the energies of features in photoemission and inverse photoemission spectra requires reexamination.

preprint2010arXiv

Effect of a Coulombic dot-lead coupling on the dynamics of a quantum dot

The effect of a Coulombic coupling on the dynamics of a quantum dot hybridized to leads is determined. The calculation treats the interaction between charge fluctuations on the dot and the dynamically generated image charge in the leads. A formally exact solution is presented for a dot coupled to a Luttinger liquid and an approximate solution, equivalent to treating the lead dynamics within a random phase approximation, is given for a dot coupled to a two- or three-dimensional metallic lead. The leading divergences arising from the long-ranged Coulomb interaction are found to cancel, so that in the two- and three-dimensional cases the quantum-dot dynamics is equivalent to that obtained by neglecting both the dot-lead Coulomb coupling and the Coulomb renormalization of the lead electrons, while in the one-dimensional case the dot-lead mixing is enhanced relative to the non-interacting case. Explicit results are given for the short-time dynamics.

preprint2010arXiv

Electron-Hole Liquids in Transition Metal Oxide Heterostructures

Appropriately designed transition metal oxide heterostructures involving small band gap Mott insulators are argued to support spatially separated electron and hole gasses at equilibrium. Spatial separations and carrier densities favoring the formation of excitonic states are achievable. The excitonic states may exhibit potentially novel properties. Energetic estimates are given, candidate material systems are discussed, and the possibility of large photvoltaic effects is mentioned

preprint2010arXiv

Momentum space anisotropy and pseudogaps: a comparative cluster dynamical mean field analysis of the doping-driven metal-insulator transition in the two dimensional Hubbard model

Cluster dynamical mean field calculations based on 2, 4, 8 and 16 site clusters are used to analyze the doping-driven metal-insulator transition in the two dimensional Hubbard model. Comparison of results obtained on different clusters enables a determination of those aspects of the physics that are common to all clusters and permits identification of artifacts associated with particular cluster geometries. A modest particle-hole asymmetry in the underlying band structure is shown to lead to qualitatively different behavior on the hole doped side than on the electron doped side. For particle-hole asymmetry of the sign and magnitude appropriate to high-$T_c$ cuprates, the approach to the insulator from the hole-doping side is found to proceed in two stages from a high-doping region where the properties are those of a Fermi liquid with moderately renormalized parameters and very weak momentum dependence. As doping is reduced the system first enters an intermediate doping regime where the Fermi liquid renormalizations are larger and the electron self energy varies significantly around the Fermi surface and then passes to a small doping regime characterized by a gap in some regions of the Fermi surface but gapless behavior in other regions. On the electron doped side the partially gapped regime does not occur, and the momentum dependence of the electron self energy is less pronounced. Implications for the high-$T_c$ cuprates and for the use of cluster dynamical mean field methods in wider classes of problems are discussed.

preprint2010arXiv

Numerically exact path integral simulation of nonequilibrium quantum transport and dissipation

We develop an iterative, numerically exact approach for the treatment of nonequilibrium quantum transport and dissipation problems that avoids the real-time sign problem associated with standard Monte Carlo techniques. The method requires a well-defined decorrelation time of the non-local influence functional for proper convergence to the exact limit. Since finite decorrelation times may arise either from temperature or from a voltage drop at zero temperature, the approach is well suited for the description of the real-time dynamics of single-molecule devices and quantum dots driven to a steady-state via interaction with two or more electron leads. We numerically investigate two non-trivial models: the evolution of the nonequilibrium population of a two-level system coupled to two electronic reservoirs, and quantum transport in the nonequilibrium Anderson model. For the latter case, two distinct formulations are described. Results are compared to those obtained by other techniques.

preprint2010arXiv

Physics of the Pseudogap in 8-site Cluster Dynamical Mean Field Theory: photoemission, Raman scattering, in-plane and c-axis conductivity

Cluster dynamical mean field and maximum entropy analytical continuation methods are used to obtain theoretical estimates for the many-body density of states, electron self-energy, in-plane and c-axis optical conductivity and the $B_{1g}$ and $B_{2g}$ Raman scattering spectra of the two dimensional square lattice Hubbard model at intermediate interaction strengths and carrier concentrations near half filling. The calculations are based on an 8-site cluster approximation which gives access to both the zone-diagonal and zone-face portions of the Fermi surface. At low dopings the zone-face regions exhibit a 'pseudogap', a suppression of the many-body density of states for energies near the Fermi surface. The pseudogap magnitude is largest near half filling and decreases smoothly with doping, but as temperature is increased the gap fills in rather than closes. The calculated response functions bear a strong qualitative resemblance to data taken in the pseudogap regime of high-$T_c$ cuprates, strongly suggesting that the intermediate coupling Hubbard model accounts for much of the exotic behavior observed in high-$T_c$ materials.

preprint2010arXiv

Theory of charged impurities in correlated electron materials: Application to muon spectroscopy of high-Tc superconductors

Single-site and cluster dynamical mean-field methods are used to estimate the response of a doped Mott insulator to a charged impurity. The effect of correlations on the Thomas-Fermi screening properties is determined. The charge density, the on-site and near-neighbor spin-spin correlations in the vicinity of the impurity are compared to those far from it. The theory is used to address the question of the effect of the density perturbation induced by the muon charge on the local response functions of a high-temperature superconductor. For reasonable values of the background dielectric constant and basic correlation strength, a muon is shown to lead to an observable perturbation of the local spin dynamics, raising questions about the interpretation of muon-spin-rotation experiments in metallic high-temperature superconductors.

preprint2010arXiv

Transport through a quantum dot with excitonic dot-lead coupling

We study the effect of a dot-lead interaction on transport through a quantum dot hybridized to two semi-infinite Luttinger-liquid leads. A bosonization approach is applied to treat the interaction between charge fluctuations on the dot and the dynamically generated image charge in the leads. The nonequilibrium distribution function of the dot and the tunneling current are computed within a master-equation approach. The presence of the excitonic dot-lead coupling is found to enhance transport in the vicinity of the Coulomb-blockade threshold. This behavior is in contrast to the usual power-law suppression of electronic tunneling which is found if this interaction is ignored.

preprint2009arXiv

GW approach to Anderson model out of equilibrium: Coulomb blockade and false hysteresis in the I-V characteristics

The Anderson model for a single impurity coupled to two leads is studied using the $GW$ approximation in the strong electron-electron interaction regime as a function of the alignment of the impurity level relative to the chemical potentials in the leads. We employ a non-equilibrium Green's function technique to calculate the electron self-energy, the spin density and the current as a function of bias across the junction. In addition we develop an expression for the change in the expectation value of the energy of the system that results when the impurity is coupled to the leads, including the role of Coulomb interactions through the electron self energy in the region of the junction. The current-voltage characteristics calculated within the GW approximation exhibit Coulomb blockade. Depending on the gate voltage and applied bias, we find that there can be more than one steady-state solution for the system, which may give rise to a hysteresis in the I-V characteristics. We show that the hysteresis is an artifact of the $GW$ approximation and would not survive if quantum fluctuations beyond the $GW$ approximation are included.

preprint2009arXiv

Magnetism and orbital-ordering in an interacting three band model: a dynamical mean field study

Single-site dynamical mean field theory is used to determine the magnetic and orbital-ordering phase diagram for a model of electrons moving on a lattice with three orbital states per site and with the fully rotationally invariant Slater-Kanamori on-site interactions. The model captures important aspects of the physics of transition metal oxides with partially filled $t_{2g}$ shells, and of electron-doped C$_{60}$. We introduce an unbiased, computationally simple and inexpensive method for estimating the presence of two sublattice order, determine the regimes in which spatially uniform and two-sublattice spin and orbital orderings are present and give physical arguments for the origins of the different phases. Guidelines are determined for optimizing the presence of ferromagnetism, which may be desirable in applications

preprint2009arXiv

Quantum criticality and non-Fermi-liquid behavior in a two-level two-lead quantum dot

Analytical and continuous-time quantum Monte Carlo methods are used to investigate the possibility of occupation switching and quantum criticality in a model of two quantum impurities coupled to two leads. A general discussion of potential occupancy-switching related quantum critical points is given, and a detailed analysis is made of a specific model which has been recently discussed. For spinless electrons, no phase transition is found. For electrons with spin, a critical value of the interaction strength separates a weak coupling regime in which all properties vary smoothly with parameters from a strong coupling phase in which occupation numbers vary discontinuously as level energies are changed. The discontinuity point is characterized by non-Fermi-liquid behavior. Results for self-energies and correlation functions are given. Phase diagrams are presented.

preprint2009arXiv

Weak-coupling quantum Monte Carlo calculations on the Keldysh contour: theory and application to the current-voltage characteristics of the Anderson model

We present optimized implementations of the weak-coupling continuous-time Monte Carlo method defined for nonequilibrium problems on the Keldysh contour. We describe and compare two methods of preparing the system before beginning the real-time calculation: the "interaction quench" and the "voltage quench", which are found to be suitable for large and small voltage biasses, respectively. We also discuss technical optimizations which increase the efficiency of the real-time measurements. The methods allow the accurate simulation of transport through quantum dots over wider interaction ranges and longer times than have heretofore been possible. The current-voltage characteristics of the particle-hole symmetric Anderson impurity model is presented for interactions U up to 10 times the intrinsic level width Gamma. We compare the Monte Carlo results to fourth order perturbation theory, finding that perturbation theory begins to fail at U/Gamma>4. Within the parameter range studied we find no evidence for a splitting of the Kondo resonance due to the applied voltage. The interplay of voltage and temperature and the Coulomb blockade conductance regime are studied.

preprint2008arXiv

The metal insulator transition in cluster dynamical mean field theory: intersite correlation, cluster size, interaction strength, and the location of the transition line

To gain insight into the physics of the metal insulator transition and the effectiveness of cluster dynamical mean field theory (DMFT) we have used one, two and four site dynamical mean field theory to solve a polaron model of electrons coupled to a classical phonon field. The cluster size dependence of the metal to polaronic insulator phase boundary is determined along with electron spectral functions and cluster correlation functions. Pronounced cluster size effects start to occur in the intermediate coupling region in which the cluster calculation leads to a gap and the single-site approximation does not. Differences (in particular a sharper band edge) persist in the strong coupling regime. A partial density of states is defined encoding a generalized nesting property of the band structure; variations in this density of states account for differences between the dynamical cluster approximation and the cellular-DMFT implementations of cluster DMFT, and for differences in behavior between the single band models appropriate for cuprates and the multiband models appropriate for manganites. A pole or strong resonance in the self energy is associated with insulating states; the momentum dependence of the pole is found to distinguish between Slater-like and Mott-like mechanisms for metal insulator transition. Implications for the theoretical treatment of doped manganites are discussed.

preprint2007arXiv

Doping-driven Mott transition in the one-band Hubbard model

A powerful new impurity solver is shown to permit a systematic study of the doping driven Mott transition in a one-band Hubbard model within the framework of single-site dynamical mean field theory. At small dopings and large interaction strengths we are able access low enough temperatures that a reliable extrapolation to temperature T=0 may be performed, and ground state energies of insulating and metallic states may be compared. We find that the T=0 doping-driven transition is of second order and is characterized by an interaction-strength dependent electronic compressibility, which vanishes at the critical interaction strength of the half filled model. Over wide parameter ranges the compressibility is substantially reduced relative to the non-interacting system. The metal insulator transition is characterized by the appearance of in-gap states, but these are relevant only for very low dopings of less than 3%.

preprint2007arXiv

Non-Equilibrium Quantum Dissipation

Dissipative processes in non-equilibrium many-body systems are fundamentally different than their equilibrium counterparts. Such processes are of great importance for the understanding of relaxation in single molecule devices. As a detailed case study, we investigate here a generic spin-fermion model, where a two-level system couples to two metallic leads with different chemical potentials. We present results for the spin relaxation rate in the nonadiabatic limit for an arbitrary coupling to the leads, using both analytical and exact numerical methods. The non-equilibrium dynamics is reflected by an exponential relaxation at long times and via complex phase shifts, leading in some cases to an "anti-orthogonality" effect. In the limit of strong system-lead coupling at zero temperature we demonstrate the onset of a Marcus-like Gaussian decay with {\it voltage difference} activation. This is analogous to the equilibrium spin-boson model, where at strong coupling and high temperatures the spin excitation rate manifests temperature activated Gaussian behavior. We find that there is no simple linear relationship between the role of the temperature in the bosonic system and a voltage drop in a non-equilibrium electronic case. The two models also differ by the orthogonality-catastrophe factor existing in a fermionic system, which modifies the resulting lineshapes. Implications for current characteristics are discussed. We demonstrate the violation of pair-wise Coulomb gas behavior for strong coupling to the leads. The results presented in this paper form the basis of an exact, non-perturbative description of steady-state quantum dissipative systems.

preprint2007arXiv

Performance analysis of continuous-time solvers for quantum impurity models

Impurity solvers play an essential role in the numerical investigation of strongly correlated electrons systems within the "dynamical mean field" approximation. Recently, a new class of continuous-time solvers has been developed, based on a diagrammatic expansion of the partition function in either the interactions or the impurity-bath hybridization. We investigate the performance of these two complementary approaches and compare them to the well-established Hirsch-Fye method. The results show that the continuous-time methods, and in particular the version which expands in the hybridization, provide substantial gains in computational efficiency.

preprint2005arXiv

Electronic reconstruction in correlated electron heterostructures

Electronic phase behavior in correlated-electron systems is a fundamental problem of condensed matter physics. We argue here that the change in the phase behavior near the surface and interface, i.e., {\em electronic reconstruction}, is the fundamental issue of the correlated-electron surface or interface science. Beyond its importance to basic science, understanding of this behavior is crucial for potential devices exploiting the novel properties of the correlated systems. % We present a general overview of the field, and then illustrate the general concepts by theoretical studies of the model heterostructures comprised of a Mott-insulator and a band-insulator, which show that spin (and orbital) orderings in thin heterostructures are generically different from the bulk and that the interface region, about three-unit-cell wide, is always metallic, demonstrating that {\em electronic reconstruction} generally occurs. % Predictions for photoemission experiments are made to show how the electronic properties change as a function of position, and the magnetic phase diagram is determined as a function of temperature, number of layers, and interaction strength. Future directions for research are also discussed.

preprint2004arXiv

Electron-lattice coupling, orbital stability and the phase diagram of Ca$_{2-x}$Sr$_x$RuO$_4$

Hartree-Fock calculations are presented of a theoretical model describing the Sr/CaRuO$_4$ family of compounds. Both commensurate and incommensurate magnetic states are considered, along with orbital ordering and the effect of lattice distortions. For reasonable parameter values, interactions disfavor orbital disproportionation. A coherent description of the observed phase diagram is obtained.

Andrew J. Millis

What is connected

Connect this record

See the researcher in context

Building this map preview

107 published item(s)

Dynamical instability in a Floquet-Driven Dissipative System

Superfluid stiffness bounds in time-reversal symmetric superconductors

Atomically imprinted graphene plasmonic cavities

Common microscopic origin of the phase transitions in Ta$_2$NiS$_5$ and the excitonic insulator candidate Ta$_2$NiSe$_5$

Dynamical Mean Field Theory of Moiré Bilayer Transition Metal Dichalcogenides: Phase Diagram, Resistivity, and Quantum Criticality

Machine Learning for Optical Scanning Probe Nanoscopy

Nonequilibrium phase transition in a driven-dissipative quantum antiferromagnet

Optical Conductivity of the two dimensional Hubbard model: vertex corrections, emergent Galilean invariance and the accuracy of the single-site dynamical mean field approximation

Second-order magnetic responses in quantum magnets: Magnetization under ac magnetic fields

Superconductivity and Antiferromagnetism in NdNiO$_2$ and CaCuO$_2$: A Cluster DMFT Study

Transient trapping into metastable states in systems with competing orders

Exploring nonequilibrium phases of photo-doped Mott insulators with Generalized Gibbs ensembles

Moiré heterostructures as a condensed matter quantum simulator

Strongly correlated ladders in K-doped $p$-terphenyl crystals

Topological charge pumping in excitonic insulators

Trigonal Symmetry Breaking and its Electronic Effects in Two-Dimensional Dihalides and Trihalides

Bardasis-Schrieffer polaritons in excitonic insulators

Charge stiffness and long-range correlation in the optically induced $η$-pairing state of the one-dimensional Hubbard model

Collective modes and terahertz near-field response of superconductors

Collective Modes in Excitonic Insulators: Effects of Electron-Phonon Coupling and Signatures in Optical Response

Ground-state properties of the hydrogen chain: insulator-to-metal transition, dimerization, and magnetic phases

Nature of symmetry breaking at the excitonic insulator transition: Ta$_2$NiSe$_5$

Negative absolute conductivity in photoexcited metals

Nonlinear spectroscopy of collective modes in excitonic insulator

Pump induced motion of an interface between competing orders

The spontaneous symmetry breaking in Ta$_2$NiSe$_5$ is structural in nature

Metal-Insulator and Magnetic Phase Diagram of Ca$_2$RuO$_4$ from Auxiliary Field Quantum Monte Carlo and Dynamical Mean Field Theory

Antisite defects at oxide interfaces

Charge transfer and electron-phonon coupling in monolayer FeSe on Nb doped SrTiO$_3$

Hybrid quantum-classical approach to correlated materials

Influence of quantum confinement and strain on orbital polarization of four-layer LaNiO$_3$ superlattices: a DFT+DMFT study

Phase diagram of Sr$_{1-x}$Ba$_x$MnO$_3$ as a function of chemical doping, epitaxial strain and external pressure

Projected Regression Methods for Inverting Fredholm Integrals: Formalism and Application to Analytical Continuation

The Anderson--Holstein Model in Two Flavors of the Non--Crossing Approximation

Band Structure and Terahertz Optical Conductivity of Transition Metal Oxides: Theory and Application to CaRuO$_3$

Charge transfer across transition metal oxide interfaces: emergent conductance and new electronic structure

Density functional versus spin-density functional and the choice of correlated subspace in multi-variable effective action theories of electronic structure

Electronic correlations, magnetism and Hund's rule coupling in the ruthenium perovskites SrRuO$_3$ and CaRuO$_3$

Imaginary-time matrix product state impurity solver for dynamical mean-field theory

Machine learning for many-body physics: efficient solution of dynamical mean-field theory

Metal-insulator transitions in GdTiO3/SrTiO3 superlattices

Photoinduced phase transitions in narrow-gap Mott insulators: the case of VO$_2$

Quasiparticle Properties of the Superconducting State of the Two Dimensional Hubbard Model

Solutions of the Two Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms

Spatial Correlations and the Insulating Phase of the High-$T_c$ Cuprates: Insights from a Configuration-Interaction-Based Solver for Dynamical Mean Field Theory

Spin-density functional theories and their $+U$ and $+J$ extensions: a comparative study of transition metals and transition metal oxides

Strain Control of Electronic Phase in Rare Earth Nickelates

Taming the dynamical sign problem in real-time evolution of quantum many-body problems

Computing total energies in complex materials using charge self-consistent DFT+DMFT

Covalency and the metal-insulator transition in titanate and vanadate perovskites

Density functional plus dynamical mean field theory of the metal-insulator transition in early transition metal oxides

Determination of Intrinsic Magnetic Response from Local Measurements of Fringing Fields

Effect of dilute, strongly pinning impurities on charge density waves

Emergent properties hidden in plane view: Strong electronic correlations at oxide interfaces

Experimental evidence for a Bragg glass density wave phase in a transition-metal dichalcogenide

Green's functions from real-time bold-line Monte Carlo

Green's functions from real-time bold-line Monte Carlo: spectral properties of the nonequilibrium Anderson impurity model

Machine learning for many-body physics: The case of the Anderson impurity model

Pairing Glue in the Two Dimensional Hubbard Model

Selectively Localized Wannier Functions

Three Dimensional Metallic and Two Dimensional Insulating Behavior in Tantalum Dichalcogenides

Total energy calculations using DFT+DMFT: computing the pressure phase diagram of the rare earth nickelates

Charge density distribution and optical response of the LaAlO3/SrTiO3 interface

Designing ferromagnetism in vanadium-oxide based superlattices

Engineering Correlation Effects via Artificially Designed Oxide Superlattices

Nematicity as a probe of superconducting pairing in iron-based superconductors

Sum rule violation in self-consistent hybridization expansions

Superconducting and Pseudogap effects on the interplane conductivity and Raman scattering cross section in the two dimensional Hubbard Model

Superconductivity and the Pseudogap in the two-dimensional Hubbard model

Suppression of superconductivity by Neel-type magnetic fluctuations in the iron pnictides

The failure of DFT-based computations for a stepped-substrate-supported correlated Co wire

Theory of ferromagnetism in vanadium-oxide based perovskites

Covalency, double-counting and the metal-insulator phase diagram in transition metal oxides

Energetics of superconductivity in the two dimensional Hubbard model