Numerical and analytical approaches to strongly correlated electron systems

In this thesis we consider three different models for strongly correlated electrons, namely a multi-band Hubbard model as well as the spinless Falicov-Kimball model, both with a semi-elliptical density of states in the limit of infinite dimensions d, and the attractive Hubbard model on a square lattice in d=2. In the first part, we study a two-band Hubbard model with unequal bandwidths and anisotropic Hund's rule coupling (J_z-model) in the limit of infinite dimensions within the dynamical mean-field theory (DMFT). Here, the DMFT impurity problem is solved with the use of quantum Monte Carlo (QMC) simulations. Our main result is that the J_z-model describes the occurrence of an orbital-selective Mott transition (OSMT), in contrast to earlier findings. We investigate the model with a high-precision DMFT algorithm, which was developed as part of this thesis and which supplements QMC with a high-frequency expansion of the self-energy. The main advantage of this scheme is the extraordinary accuracy of the numerical solutions, which can be obtained already with moderate computational effort, so that studies of multi-orbital systems within the DMFT+QMC are strongly improved. We also found that a suitably defined Falicov-Kimball (FK) model exhibits an OSMT, revealing the close connection of the Falicov-Kimball physics to the J_z-model in the OSM phase. In the second part of this thesis we study the attractive Hubbard model in two spatial dimensions within second-order self-consistent perturbation theory. This model is considered on a square lattice at finite doping and at low temperatures. Our main result is that the predictions of first-order perturbation theory (Hartree-Fock approximation) are renormalized by a factor of the order of unity even at arbitrarily weak interaction (U->0). The renormalization factor q can be evaluated as a function of the filling n for 00, the q-factor vanishes, signaling the divergence of self-consistent perturbation theory in this limit. Thus we present the first asymptotically exact results at weak-coupling for the negative-U Hubbard model in d=2 at finite doping.

Theory of strongly correlated electron systems. I. Intersite coulomb interaction and the approximation of renormalized fermions in total energy calculations

The diagrammatic strong-coupling perturbation theory (SCPT) for correlated electron systems is developed for intersite Coulomb interaction and for a nonorthogonal basis set. The construction is based on iterations of exact closed equations for many - electron Green functions (GFs) for Hubbard operators in terms of functional derivatives with respect to external sources. The graphs, which do not contain the contributions from the fluctuations of the local population numbers of the ion states, play a special role: a one-to-one correspondence is found between the subset of such graphs for the many - electron GFs and the complete set of Feynman graphs of weak-coupling perturbation theory (WCPT) for single-electron GFs. This fact is used for formulation of the approximation of renormalized Fermions (ARF) in which the many-electron quasi-particles behave analogously to normal Fermions. Then, by analyzing: (a) Sham's equation, which connects the self-energy and the exchange- correlation potential in density functional theory (DFT); and (b) the Galitskii and Migdal expressions for the total energy, written within WCPT and within ARF SCPT, a way we suggest a method to improve the description of the systems with correlated electrons within the local density approximation (LDA) to DFT. The formulation, in terms of renormalized Fermions LIDA (RF LDA), is obtained by introducing the spectral weights of the many electron GFs into the definitions of the charge density, the overlap matrices, effective mixing and hopping matrix elements, into existing electronic structure codes, whereas the weights themselves have to be found from an additional set of equations. Compared with LDA+U and self-interaction correction (SIC) methods, RF LDA has the advantage of taking into account the transfer of spectral weights, and, when formulated in terms of GFs, also allows for consideration of excitations and nonzero temperature. Going beyond the ARF SCPT, as well as RF LIDA, and taking into account the fluctuations of ion population numbers would require writing completely new codes for ab initio calculations. The application of RF LDA for ab initio band structure calculations for rare earth metals is presented in part 11 of this study (this issue). (c) 2005 Wiley Periodicals, Inc.

Theory of strongly correlated electron systems. II. Including correlation effects into electronic structure calculations

We have previously shown that a division of the f-shell into two subsystems gives a better understanding of the cohesive properties as well the general behavior of lanthanide systems. In this article, we present numerical computations, using the suggested method. We show that the picture is consistent with most experimental data, e.g., the equilibrium volume and electronic structure in general. Compared with standard energy band calculations and calculations based on the self-interaction correction and LIDA + U, the f-(non-f)-mixing interaction is decreased by spectral weights of the many-body states of the f-ion. (c) 2005 Wiley Periodicals, Inc.

Quantum Properties in Transport Through Nanoscopic Rings:Charge-Spin Separation and Interference Effects

Many of the most intriguing quantum effects are observed or could be measured in transport experiments through nanoscopic systems such as quantum dots, wires and rings formed by large molecules or arrays of quantum dots. In particular, the separation of charge and spin degrees of freedom and interference effects have important consequences in the conductivity through these systems. Charge-spin separation was predicted theoretically in one-dimensional strongly inter-acting systems (Luttinger liquids) and, although observed indirectly in several materials formed by chains of correlated electrons, it still lacks direct observation. We present results on transport properties through Aharonov-Bohmrings (pierced by a magnetic flux) with one or more channels represented by paradigmatic strongly-correlated models. For a wide range of parameters we observe characteristic dips in the conductance as a function of magnetic flux which are a signature of spin and charge separation. Interference effects could also be controlled in certain molecules and interesting properties could be observed. We analyze transport properties of conjugated molecules, benzene in particular, and find that the conductance depends on the lead configuration. In molecules with translational symmetry, the conductance can be controlled by breaking or restoring this symmetry, e.g. by the application of a local external potential. These results open the possibility of observing these peculiar physical properties in anisotropic ladder systems and in real nanoscopic and molecular devices.

Double-electron above threshold ionization of helium

We present calculations of intense-field multiphoton ionization processes in helium at XUV wavelengths. The calculations are obtained from a full-dimensional integration of the two-electron time-dependent Schrödinger equation. A momentum-space analysis of the ionizing two-electron wavepacket reveals the existence of double-electron above threshold ionization (DATI). In momentum-space two distinct forms of DATI are resolved, namely non-sequential and sequential. In non-sequential DATI correlated electrons resonantly absorb and share energy in integer units of Ïlaser.

Analyse des propriétés électroniques de supraconducteurs à l’aide de la théorie de la fonctionnelle de la densité

Cette thèse traite de la structure électronique de supraconducteurs telle que déterminée par la théorie de la fonctionnelle de la densité. Une brève explication de cette théorie est faite dans l’introduction. Le modèle de Hubbard est présenté pour pallier à des problèmes de cette théorie face à certains matériaux, dont les cuprates. L’union de deux théories donne la DFT+U, une méthode permettant de bien représenter certains systèmes ayant des électrons fortement corrélés. Par la suite, un article traitant du couplage électron- phonon dans le supraconducteur NbC1−xNx est présenté. Les résultats illustrent bien le rôle de la surface de Fermi dans le mécanisme d’appariement électronique menant à la supraconductivité. Grâce à ces résultats, un modèle est développé qui permet d’expliquer comment la température de transition critique est influencée par le changement des fré- quences de vibration du cristal. Ensuite, des résultats de calcul d’oscillations quantiques obtenus par une analyse approfondie de surfaces de Fermi, permettant une comparaison directe avec des données expérimentales, sont présentés dans deux articles. Le premier traite d’un matériau dans la famille des pnictures de fer, le LaFe2P2. L’absence de su- praconductivité dans ce matériau s’explique par la différence entre sa surface de Fermi obtenue et celle du supraconducteur BaFe2As2. Le second article traite du matériau à fermions lourds, le YbCoIn5. Pour ce faire, une nouvelle méthode efficace de calcul des fréquences de Haas-van Alphen est développée. Finalement, un dernier article traitant du cuprate supraconducteur à haute température critique YBa2Cu3O6.5 est présenté. À l’aide de la DFT+U, le rôle de plusieurs ordres magnétiques sur la surface de Fermi est étudié. Ces résultats permettent de mieux comprendre les mesures d’oscillations quan- tiques mesurées dans ce matériau.

Interacting fermionic atoms in optical lattices - A quantum simulator for condensed matter physics

This thesis reports on the creation and analysis of many-body states of interacting fermionic atoms in optical lattices. The realized system can be described by the Fermi-Hubbard hamiltonian, which is an important model for correlated electrons in modern condensed matter physics. In this way, ultra-cold atoms can be utilized as a quantum simulator to study solid state phenomena. The use of a Feshbach resonance in combination with a blue-detuned optical lattice and a red-detuned dipole trap enables an independent control over all relevant parameters in the many-body hamiltonian. By measuring the in-situ density distribution and doublon fraction it has been possible to identify both metallic and insulating phases in the repulsive Hubbard model, including the experimental observation of the fermionic Mott insulator. In the attractive case, the appearance of strong correlations has been detected via an anomalous expansion of the cloud that is caused by the formation of non-condensed pairs. By monitoring the in-situ density distribution of initially localized atoms during the free expansion in a homogeneous optical lattice, a strong influence of interactions on the out-of-equilibrium dynamics within the Hubbard model has been found. The reported experiments pave the way for future studies on magnetic order and fermionic superfluidity in a clean and well-controlled experimental system.

Interacting bosons and fermions in three-dimensional optical lattice potentials

This thesis reports on the realization, characterization and analysis of ultracold bosonic and fermionic atoms in three-dimensional optical lattice potentials. Ultracold quantum gases in optical lattices can be regarded as ideal model systems to investigate quantum many-body physics. In this work interacting ensembles of bosonic 87Rb and fermionic 40K atoms are employed to study equilibrium phases and nonequilibrium dynamics. The investigations are enabled by a versatile experimental setup, whose core feature is a blue-detuned optical lattice that is combined with Feshbach resonances and a red-detuned dipole trap to allow for independent control of tunneling, interactions and external confinement. The Fermi-Hubbard model, which plays a central role in the theoretical description of strongly correlated electrons, is experimentally realized by loading interacting fermionic spin mixtures into the optical lattice. Using phase-contrast imaging the in-situ size of the atomic density distribution is measured, which allows to extract the global compressibility of the many-body state as a function of interaction and external confinement. Thereby, metallic and insulating phases are clearly identified. At strongly repulsive interaction, a vanishing compressibility and suppression of doubly occupied lattice sites signal the emergence of a fermionic Mott insulator. In a second series of experiments interaction effects in bosonic lattice quantum gases are analyzed. Typically, interactions between microscopic particles are described as two-body interactions. As such they are also contained in the single-band Bose-Hubbard model. However, our measurements demonstrate the presence of multi-body interactions that effectively emerge via virtual transitions of atoms to higher lattice bands. These findings are enabled by the development of a novel atom optical measurement technique: In quantum phase revival spectroscopy periodic collapse and revival dynamics of the bosonic matter wave field are induced. The frequencies of the dynamics are directly related to the on-site interaction energies of atomic Fock states and can be read out with high precision. The third part of this work deals with mixtures of bosons and fermions in optical lattices, in which the interspecies interactions are accurately controlled by means of a Feshbach resonance. Studies of the equilibrium phases show that the bosonic superfluid to Mott insulator transition is shifted towards lower lattice depths when bosons and fermions interact attractively. This observation is further analyzed by applying quantum phase revival spectroscopy to few-body systems consisting of a single fermion and a coherent bosonic field on individual lattice sites. In addition to the direct measurement of Bose-Fermi interaction energies, Bose-Bose interactions are proven to be modified by the presence of a fermion. This renormalization of bosonic interaction energies can explain the shift of the Mott insulator transition. The experiments of this thesis lay important foundations for future studies of quantum magnetism with fermionic spin mixtures as well as for the realization of complex quantum phases with Bose-Fermi mixtures. They furthermore point towards physics that reaches beyond the single-band Hubbard model.

Drinfeld twists and algebraic bethe ansatz of the supersymmetric model associated with U-q(gl(m vertical bar n))

We construct the Drinfeld twists (or factorizing F-matrices) of the supersymmetric model associated with quantum superalgebra U-q(gl(m vertical bar n)), and obtain the completely symmetric representations of the creation operators of the model in the F-basis provided by the F-matrix. As an application of our general results, we present the explicit expressions of the Bethe vectors in the F-basis for the U-q(gl(2 vertical bar 1))-model (the quantum t-J model).

Bosonisation as The Hubbard-Stratonovich transformation

The problem of strongly correlated electrons in one dimension attracted attention of condensed matter physicists since early 50’s. After the seminal paper of Tomonaga [1] who suggested the first soluble model in 1950, there were essential achievements reflected in papers by Luttinger [2] (1963) and Mattis and Lieb [3] (1963). A considerable contribution to the understanding of generic properties of the 1D electron liquid has been made by Dzyaloshinskii and Larkin [4] (1973) and Efetov and Larkin [5] (1976). Despite the fact that the main features of the 1D electron liquid were captured and described by the end of 70’s, the investigators felt dissatisfied with the rigour of the theoretical description. The most famous example is the paper by Haldane [6] (1981) where the author developed the fundamentals of a modern bosonisation technique, known as the operator approach. This paper became famous because the author has rigourously shown how to construct the Fermi creation/anihilation operators out of the Bose ones. The most recent example of such a dissatisfaction is the review by von Delft and Schoeller [7] (1998) who revised the approach to the bosonisation and came up with what they called constructive bosonisation.

Correlated electron-ion dynamics: the excitation of atomic motion by energetic electrons

Correlated electron-ion dynamics (CEID) is an extension of molecular dynamics that allows us to introduce in a correct manner the exchange of energy between electrons and ions. The formalism is based on a systematic approximation: small amplitude moment expansion. This formalism is extended here to include the explicit quantum spread of the ions and a generalization of the Hartree-Fock approximation for incoherent sums of Slater determinants. We demonstrate that the resultant dynamical equations reproduce analytically the selection rules for inelastic electron-phonon scattering from perturbation theory, which control the mutually driven excitations of the two interacting subsystems. We then use CEID to make direct numerical simulations of inelastic current-voltage spectroscopy in atomic wires, and to exhibit the crossover from ionic cooling to heating as a function of the relative degree of excitation of the electronic and ionic subsystems.

Long-range Coulomb interactions and nanoscale electronic inhomogeneities in correlated oxides

Electronic, magnetic, or structural inhomogeneities ranging in size from nanoscopic to mesoscopic scales seem endemic and are possibly generic to colossal magnetoresistance manganites and other transition metal oxides. They are hence of great current interest and understanding them is of fundamental importance. We show here that an extension, to include long-range Coulomb interactions, of a quantum two-fluid l-b model proposed recently for manganites [Phys. Rev. Lett. 92, 157203 (2004)] leads to an excellent description of such inhomogeneities. In the l-b model two very different kinds of electronic states, one localized and polaronic (l) and the other extended or broad band (b) coexist. For model parameters appropriate to manganites and even within a simple dynamical mean-field theory (DMFT) framework, it describes many of the unusual phenomena seen in manganites, including colossal magnetoresistance (CMR), qualitatively and quantitatively. However, in the absence of long-ranged Coulomb interaction, a system described by such a model would actually phase separate, into macroscopic regions of l and b electrons, respectively. As we show in this paper, in the presence of Coulomb interactions, the macroscopic phase separation gets suppressed and instead nanometer scale regions of polarons interspersed with band electron puddles appear, constituting a kind of quantum Coulomb glass. We characterize the size scales and distribution of the inhomogeneity using computer simulations. For realistic values of the long-range Coulomb interaction parameter V-0, our results for the thresholds for occupancy of the b states are in agreement with, and hence support, the earlier approach mentioned above based on a configuration averaged DMFT treatment which neglects V-0; but the present work has features that cannot be addressed in the DMFT framework. Our work points to an interplay of strong correlations, long-range Coulomb interaction, and dopant ion disorder, all inevitably present in transition metal oxides as the origin of nanoscale inhomogeneities rather than disorder frustrated phase competition as is generally believed. As regards manganites, it argues against explanations for CMR based on disorder frustrated phase separation and for an intrinsic origin of CMR. Based on this, we argue that the observed micrometer (meso) scale inhomogeneities owe their existence to extrinsic causes, e.g., strain due to cracks and defects. We suggest possible experiments to validate our speculation.

Electronic Structure Evolution across the Peierls Metal-Insulator Transition in a Correlated Ferromagnet

Transition metal compounds often undergo spin-charge-orbital ordering due to strong electron-electron correlations. In contrast, low-dimensional materials can exhibit a Peierls transition arising from low-energy electron-phonon-coupling-induced structural instabilities. We study the electronic structure of the tunnel framework compound K2Cr8O16, which exhibits a temperature-dependent (T-dependent) paramagnetic-to-ferromagnetic- metal transition at T-C = 180 K and transforms into a ferromagnetic insulator below T-MI = 95 K. We observe clear T-dependent dynamic valence (charge) fluctuations from above T-C to T-MI, which effectively get pinned to an average nominal valence of Cr+3.75 (Cr4+:Cr3+ states in a 3:1 ratio) in the ferromagnetic-insulating phase. High-resolution laser photoemission shows a T-dependent BCS-type energy gap, with 2G(0) similar to 3.5(k(B)T(MI)) similar to 35 meV. First-principles band-structure calculations, using the experimentally estimated on-site Coulomb energy of U similar to 4 eV, establish the necessity of strong correlations and finite structural distortions for driving the metal-insulator transition. In spite of the strong correlations, the nonintegral occupancy (2.25 d-electrons/Cr) and the half-metallic ferromagnetism in the t(2g) up-spin band favor a low-energy Peierls metal-insulator transition.

Density functional theory embedding for correlated wavefunctions

Methods that exploit the intrinsic locality of molecular interactions show significant promise in making tractable the electronic structure calculation of large-scale systems. In particular, embedded density functional theory (e-DFT) offers a formally exact approach to electronic structure calculations in which the interactions between subsystems are evaluated in terms of their electronic density. In the following dissertation, methodological advances of embedded density functional theory are described, numerically tested, and applied to real chemical systems.

First, we describe an e-DFT protocol in which the non-additive kinetic energy component of the embedding potential is treated exactly. Then, we present a general implementation of the exact calculation of the non-additive kinetic potential (NAKP) and apply it to molecular systems. We demonstrate that the implementation using the exact NAKP is in excellent agreement with reference Kohn-Sham calculations, whereas the approximate functionals lead to qualitative failures in the calculated energies and equilibrium structures.

Next, we introduce density-embedding techniques to enable the accurate and stable calculation of correlated wavefunction (CW) in complex environments. Embedding potentials calculated using e-DFT introduce the effect of the environment on a subsystem for CW calculations (WFT-in-DFT). We demonstrate that WFT-in-DFT calculations are in good agreement with CW calculations performed on the full complex.

We significantly improve the numerics of the algorithm by enforcing orthogonality between subsystems by introduction of a projection operator. Utilizing the projection-based embedding scheme, we rigorously analyze the sources of error in quantum embedding calculations in which an active subsystem is treated using CWs, and the remainder using density functional theory. We show that the embedding potential felt by the electrons in the active subsystem makes only a small contribution to the error of the method, whereas the error in the nonadditive exchange-correlation energy dominates. We develop an algorithm which corrects this term and demonstrate the accuracy of this corrected embedding scheme.

Hydrogen contaminant and its correlation with background electrons in GaN

GaN epilayers grown by molecular beam epitaxy using NH3 as the nitrogen source were found to contain hydrogen. We further notice that the background electron concentration in GaN can be correlated with the amount of hydrogen contaminant. X-ray photoelectron spectroscopy (XPS) measurements of the N Is peak reveal that hydrogen is bound to nitrogen. This will make the corresponding Ga atom see insufficient N counterpart, as can be inferred from the XPS Ga 3d spectrum. We then think that nitrogen in the lattice terminated by hydrogen is an effective nitrogen vacancy and hence a donor accounting for the background electrons.