Search

In this work, we investigate collective electronic fluctuations and, in particular, the possibility of the charge density wave ordering in an infinite-layer NdNiO$_2$. We perform advanced many-body calculations for the ab-initio three-orbital model by taking into account local correlation effects, non-local charge and magnetic fluctuations, and the electron-phonon coupling. We find that in the considered material, electronic correlations are strongly orbital- and momentum-dependent. Notably, the charge density wave and magnetic instabilities originate from distinct orbitals. In particular, we show that the correlation effects lead to the momentum-dependent hybridization between different orbitals, resulting in the splitting and shifting of the flat part of the Ni-$d_{z^2}$ band. This strong renormalization of the electronic spectral function drives the charge density wave instability that is related to the intraband Ni-$d_{z^2}$ correlations. Instead, the magnetic instability stems from the Ni-$d_{x^2-y^2}$ orbital, which remains half-filled through the redistribution of the electronic density between different bands even upon hole doping. Consequently, the strength of the magnetic fluctuations remains nearly unchanged for the considered doping levels. We argue that this renormalization is not inherent to the stoichiometric case but can be induced by hole doping.

Superconductivity emerges from the spatial coherence of a macroscopic condensate of Cooper pairs. Increasingly strong binding and localization of electrons into these pairs compromises the condensate's phase stiffness, thereby limiting critical temperatures -- a phenomenon known as the BCS-BEC crossover in lattice systems. In this study, we demonstrate enhanced superconductivity in a multiorbital model of alkali-doped fullerides (A$_3$C$_{60}$) that goes beyond the limits of the lattice BCS-BEC crossover. We identify that the interplay of strong correlations and multiorbital effects results in a localized superconducting state characterized by a short coherence length but robust stiffness and a domeless rise in critical temperature with increasing pairing interaction. To derive these insights, we introduce a new theoretical framework allowing us to calculate the fundamental length scales of superconductors, namely the coherence length ($\xi_0$) and the London penetration depth ($\lambda_{\mathrm{L}}$), even in presence of strong electron correlations., Comment: main: 12 pages, 5 figures, 1 table | SM: 25(+1), 7 figures, 1 table

The electronic structure and magnetism of individual Dy atom adsorbed on the MgO(001) substrate is investigated using the combination of the density functional theory with the Hubbard-I approximation to the Anderson impurity model (DFT+U(HIA)). The divalent Dy$^{2+}$ adatom in $f^{10}$ configuration is found. The calculated x-ray absorption (XAS) and magnetic circular dichroism (XMCD) spectra are compared to the experimental data. Quantum tunneling between degenerate $|{J=8.0, J_z= \pm 4.0}>$ states leads to formation of $|{J=8.0, J_z= 0.0}>$ ground state with an in-plane orientation of the magnetic moment. It explains absence of remanent magnetization in MgO adatom on the top of Mg(001) substrate. Our studies can provide a viable route for further investigation and prediction of the rare-earth single atom magnets., Comment: 6 pages, 5 figures

The zero-bias anomaly in conductance spectra of single Co atoms on Cu(111) observed at $\approx$ 4 K, which has been interpreted as being due to a Kondo resonance, is strongly modified when the Co atoms are attached to monatomic Cu chains. Scanning tunneling spectra measured at 340 mK in magnetic fields exhibit all characteristics of spin-flip excitations. Their dependence on the magnetic field reveals a magnetic anisotropy and suggests a non-collinear spin state. This indicates that spin-orbit coupling (SOC), which has so far been neglected in theoretical studies of Co/Cu(111), has to be taken into account. According to our density functional theory and multi-orbital quantum Monte Carlo calculations SOC suppresses the Kondo effect for all studied geometries.

We develop a strong-coupling perturbation scheme for a generic Hubbard model around a half-filled particle-hole-symmetric reference system, which is free from the fermionic sign problem. The approach is based on the lattice determinantal Quantum Monte Carlo (QMC) method in continuous and discrete time versions for large periodic clusters in a fermionic bath. Considering the first-order perturbation in the shift of the chemical potential and of the second-neighbour hopping gives an accurate electronic spectral function for a parameter range corresponding to the optimally doped cuprate system for temperature of the order of $T=0.1t$, the region hardly accessible for the straightforward lattice QMC calculations. We discuss the formation of the pseudogap and the nodal-antinodal dichotomy for a doped Hubbard system in a strong-coupling regime with the interaction parameter $U$ equal to the bandwidth and the optimal value of the next-nearest-neighbor hopping parameter $t'$ for high-temperature superconducting cuprates., Comment: 14 pages, 17 figures

In this work we computed the phase diagram as a function of temperature and doping for a system of lead adatoms allocated periodically on a silicon (111) surface. This Si(111):Pb material is characterized by a strong and long-ranged Coulomb interaction, a relatively large value of the spin-orbit coupling, and a structural phase transition that occurs at low temperature. In order to describe the collective electronic behavior in the system, we perform many-body calculations consistently taking all these important features into account. We find that charge- and spin-density wave orderings coexist with each other in several regions of the phase diagram. This result is in agreement with the recent experimental observation of a chiral spin texture in the charge density wave phase in this material. We also find that geometries of the charge and spin textures strongly depend on the doping level. The formation of such a rich phase diagram in the Si(111):Pb material can be explained by a combined effect of the lattice distortion and electronic correlations.

In this report we review the method of explicit calculations of interatomic exchange interactions of magnetic materials. This involves exchange mechanisms normally referred to as Heisenberg exchange, Dzyaloshinskii-Moriya interaction and anisotropic symmetric exchange. The connection between microscopic theories of the electronic structure, such as density functional theory or dynamical mean field theory, and interatomic exchange, is given in detail. The different aspects of extracting information for an effective spin Hamiltonian that involves thousands of atoms, from electronic structure calculations considering significantly fewer atoms (1-50) is highlighted. Examples of exchange interactions of a large group of materials is presented, which involves heavy elements of the 3d period, alloys between transition metals, Heusler compounds, multilayer systems as well as overlayers and adatoms on a substrate, transition metal oxides, 4f elements, magnetic materials in two dimensions and molecular magnets. Where possible, a comparison to experimental data is made, that naturally becomes focused on the magnon dispersion. The influence of relativity is reviewed for a few cases, as is the importance of dynamical correlations. Development to theories that handle out of equilibrium conditions is also described here. The review ends with a short description of extensions of the theories behind explicit calculations of interatomic exchange, to non-magnetic situations, e.g. that describe chemical (charge) order and superconductivity.

Many strongly correlated transition metal insulators are colored, even though they have large fundamental band gaps and no quasi-particle excitations in the visible range. Why such insulators possess the colors they do poses a serious challenge for any many-body theory to reliably pick up the interactions responsible for the color. We pick two archetypal cases as examples: NiO with green color and MnF\textsubscript{2} with pink color. The body of literature around the collective charge transitions (excitons) that are responsible for the color in these and other strongly correlated systems, often fail to disentangle two important factors: what makes them form and what makes them optically bright. An adequate answer requires a theoretical approach able to compute such excitations in periodic crystals, reliably and without free parameters -- a formidable challenge. We employ two kinds of advanced \emph{ab initio} many body Green's function theories to investigate both optical and spin susceptibilities. The first, a perturbative theory based on low-order extensions of the $GW$ approximation, is able to explain the color in NiO, and indeed well describe the dielectric response over the entire frequency spectrum, while the same theory is unable to explain why MnF\textsubscript{2} is pink. We show its color originates from higher order spin-flip transitions that modify the optical response. This phenomenon is not captured by low-order perturbation theory, but it is contained in dynamical mean-field theory (DMFT), which has a dynamical spin-flip vertex that contributes to the charge susceptibility. We show that symmetry lowering mechanisms, such as spin-orbit coupling, odd-parity phonons and Jan-Teller distortions, determine how `bright' these excitons are, but are not fundamental to their existence., Comment: 19 pages and 8 figures

We present the multi-band dual triply irreducible local expansion (D-TRILEX) approach to interacting electronic systems and discuss its numerical implementation. This method is designed for a self-consistent description of multi-orbital systems that can also have several atoms in the unit cell. The current implementation of the D-TRILEX approach is able to account for the frequency- and channel-dependent long-ranged electronic interactions. We show that our method is accurate when applied to small multi-band systems such as the Hubbard-Kanamori dimer. Calculations for the extended Hubbard, the two-orbital Hubbard-Kanamori, and the bilayer Hubbard models are also discussed.

A clear separation of the time scales governing the dynamics of "slow" and "fast" degrees of freedom often serves as a prerequisite for the emergence of an independent low-energy theory. Here, we consider (slow) classical spins exchange coupled to a tight-binding system of (fast) conduction electrons. The effective equations of motion are derived under the constraint that the quantum state of the electron system at any instant of time $t$ lies in the $n$-dimensional low-energy subspace for the corresponding spin configuration at $t$. The effective low-energy theory unfolds itself straightforwardly and takes the form of a non-abelian gauge theory with the gauge freedom given by the arbitrariness of the basis spanning the instantaneous low-energy sector. The holonomic constraint generates a gauge covariant spin-Berry curvature tensor in the equations of motion for the classical spins. In the non-abelian theory for $n>1$, opposed to the $n=1$ adiabatic spin dynamics theory, the spin-Berry curvature is generically nonzero, even for time-reversal symmetric systems. Its expectation value with the representation of the electron state is gauge invariant and gives rise to an additional {\em geometrical} spin torque. Besides anomalous precession, the $n\ge 2$ theory also captures the spin nutational motion, which is usually considered as a retardation effect. This is demonstrated by proof-of-principle numerical calculations for a minimal model with a single classical spin. Already for $n=2$ and in parameter regimes where the $n=1$ adiabatic theory breaks down, we find good agreement with results obtained from the full (unconstrained) theory., Comment: 19 pages, 12 figures, v2 with minor changes

In the study of correlated systems, approximations based on the dynamical mean-field theory (DMFT) provide a practical way to take local vertex corrections into account, which capture, respectively, particle-particle screening at weak coupling and the formation of the local moment at strong coupling. We show that in both limits the local vertex corrections can be efficiently parametrized in terms of single-boson exchange, such that the two-particle physics described by DMFT and its diagrammatic extensions is recovered to good approximation and at a reduced computational cost. Our investigation highlights the importance of the frequency-dependent fermion-boson coupling (Hedin vertex) for local vertex corrections. Namely, at weak coupling the fermion-spin-boson coupling suppresses the N\'eel temperature of the DMFT approximation compared to the static mean-field, whereas for large interaction it facilitates a huge enhancement of local spin-fluctuation exchange, giving rise to the effective-exchange energy scale $4t^2/U$. We find that parametrizations of the vertex which neglect the nontrivial part of the fermion-boson coupling fail qualitatively at strong coupling.

We apply the boson exchange parquet solver for dual fermions to the half-filled Hubbard model on a square lattice at small interaction. Our results establish that, in this regime, nonlocal vertex corrections play an important role in the formation of the pseudogap. Namely, in comparison to the simpler ladder approximation, these additional vertex corrections included in the parquet equations enhance the coupling of spin fluctuations with the quasiparticles. The pseudogap thus opens already at a higher temperature, in quantitative agreement with the numerically exact diagrammatic Monte Carlo. The representation of the parquet diagrams in terms of boson exchange facilitates large lattice sizes and gives rise to an unbiased fluctuation diagnostic of the self-energy, which does not rely on the Fierz ambiguity. The fluctuation diagnostic implies that nodal and antinodal fermions are affected equally by spin fluctuations with the exact commensurate nesting vector $(\pi,\pi)$. However, the antinode couples more efficiently to incommensurate fluctuations than the node, leading to the nodal/antinodal dichotomy. We corroborate this finding in terms of a spin-fermion-like calculation., Comment: 13 pages, 12 figures

Characterizing non-local magnetic fluctuations in materials with strong electronic Coulomb interactions remains one of the major outstanding challenges of modern condensed matter theory. In this work we address the spatial symmetry and orbital structure of magnetic fluctuations in perovskite materials. To this aim, we develop a consistent multi-orbital diagrammatic extension of dynamical mean field theory, which we apply to an anisotropic three-orbital model of cubic $t_{2g}$ symmetry. We find that the form of spatial spin fluctuations is governed by the local Hund's coupling. For small values of the coupling, magnetic fluctuations are anisotropic in orbital space, which reflects the symmetry of the considered $t_{2g}$ model. Large Hund's coupling enhances collective spin excitations, which mixes orbital and spatial degrees of freedom, and magnetic fluctuations become orbitally isotropic. Remarkably, this effect can be seen only in two-particle quantities; single-particle observables remain anisotropic for any value of the Hund's coupling. Importantly, we find that the orbital isotropy can be induced both, at half-filling and for the case of $4$ electrons per lattice site, where the magnetic instability is associated with different, antiferromagnetic and ferromagnetic modes, respectively., Comment: 6 pages + 3 pages of Supplemental Material

We present and implement a parquet approximation within the dual-fermion formalism based on a partial bosonization of the dual vertex function which substantially reduces the computational cost of the calculation. The method relies on splitting the vertex exactly into single-boson exchange contributions and a residual four-fermion vertex, which physically embody respectively long-range and short-range spatial correlations. After recasting the parquet equations in terms of the residual vertex, these are solved using the truncated unity method of Eckhardt et al. [Phys. Rev. B 101, 155104 (2020)], which allows for a rapid convergence with the number of form factors in different regimes. While our numerical treatment of the parquet equations can be restricted to only a few Matsubara frequencies, reminiscent of Astretsov et al. [Phys. Rev. B 101, 075109 (2020)], the one- and two-particle spectral information is fully retained. In applications to the two-dimensional Hubbard model the method agrees quantitatively with a stochastic summation of diagrams over a wide range of parameters., Comment: 13 pages + appendix, 15 figures

We present a purely diagrammatic derivation of the dual fermion scheme [Phys. Rev. B 77 (2008) 033101]. The derivation makes particularly clear that a similar scheme can be developed for an arbitrary reference system provided it has the same interaction term as the original system. Thereby no restrictions are imposed by the locality of the reference problem or by the nature of the original problem as a lattice one. We present new arguments in favour of keeping the dual denominator in the expression for the lattice self-energy independently of the truncation of the dual interaction. As an example we present the computational results for the half-filled 2D Hubbard model with the choice of a $2\times2$ plaquette with periodic boundary conditions as a reference system. We observe that obtained results are in a good agreement with numerically exact lattice quantum Monte Carlo data.

We establish a way to handle main collective fluctuations in correlated quantum systems based on a Fluctuation Local Field concept. This technique goes beyond standard mean-field approaches, such as Hartree-Fock and dynamical mean-field theories (DMFT), as it includes a fluctuating classical field that acts on the leading order parameter of the system. Effective model parameters of this new theory are determined from the variational principle, which allows to resolve the Fierz ambiguity in decoupling of the local interaction term. In the saddle-point approximation for the fluctuating field our method reproduces the mean-field result. The exact numerical integration over this field allows to consider nonlinear fluctuations of the global order parameter of the system while local correlations can be accounted by solving the DMFT impurity problem. We apply our method to the magnetic susceptibility of finite Hubbard systems at half-filling and demonstrate that the introduced technique leads to a superior improvement of results with respect to parental mean-field approaches without significant numerical complications. We show that the Fluctuation Local Field method can be used in a very broad range of temperatures substantially below the N\'eel temperature of DMFT, which remains a major challenge for all existing theoretical approaches.

Floquet engineering of closed quantum systems can lead to the formation of long-lived prethermal states that, in general, eventually thermalize to infinite temperature. Coupling these driven systems to dissipative baths can stabilize such states, establishing a true nonequilibrium steady state. We demonstrate that in a certain strongly interacting lattice model coupled to a bath and driven by an electric field, such steady states can have the remarkable property that the cycle-averaged rate of energy transfer between the lattice and the baths vanishes. Despite this, we show that these states retain a clear nonequilibrium nature.

We present an infinite-dimensional lattice of two-by-two plaquettes, the quadruple Bethe lattice, with Hubbard interaction and solve it exactly by means of the cluster dynamical mean-field theory. It exhibits a $d$-wave superconducting phase that is related to a highly degenerate point in the phase diagram of the isolated plaquette at that the groundstates of the particle number sectors $N=2,3,4$ cross. The superconducting gap is formed by the renormalized lower Slater peak of the correlated, hole-doped Mott insulator. We engineer parts of the interaction and find that pair hoppings between $X/Y$-momenta are the main two-particle correlations of the superconducting phase. The suppression of the superconductivity in the overdoped regime is caused by the diminishing of pair hopping correlations and in the underdoped regime by charge blocking. The optimal doping is $\sim 0.15$ at which the underlying normal state shows a Lifshitz transition. The model allows for different intra- and inter-plaquette hoppings that we use to disentangle superconductivity from antiferromagnetism as the latter requires larger inter-plaquette hoppings., Comment: 16 pages, 19 figures

In this work we introduce a modified real-time continuous-time hybridization-expansion quantum Monte Carlo solver for a time-dependent single-orbital Anderson impurity model: CT-1/2-HYB-QMC. In the proposed method the diagrammatic expansion is performed only for one out of the two spin channels, while the resulting effective single-particle problem for the other spin is solved semi-analytically for each expansion diagram. CT-1/2-HYB-QMC alleviates the dynamical sign problem by reducing the order of sampled diagrams and makes it possible to reach twice as long time scales in comparison to the standard CT-HYB method. We illustrate the new solver by calculating an electric current through impurity in paramagnetic and spin-polarized cases., Comment: 18 pages, 8 figures, Submission to SciPost

A considerable success in phenomenological description of high-T$_{\rm c}$ superconductors has been achieved within the paradigm of Quantum Critical Point (QCP) - a parental state of a variety of exotic phases that is characterized by dense entanglement and absence of well-defined quasiparticles. However, the microscopic origin of the critical regime in real materials remains an open question. On the other hand, there is a popular view that a single-band $t-t'$ Hubbard model is the minimal model to catch the main relevant physics of superconducting compounds. Here, we suggest that emergence of the QCP is tightly connected with entanglement in real space and identify its location on the phase diagram of the hole-doped $t-t'$ Hubbard model. To detect the QCP we study a weighted graph of inter-site quantum mutual information within a four-by-four plaquette that is solved by exact diagonalization. We demonstrate that some quantitative characteristics of such a graph, viewed as a complex network, exhibit peculiar behavior around a certain submanifold in the parametric space of the model. This method allows us to overcome difficulties caused by finite size effects and to identify the transition point even on a small lattice, where long-range asymptotics of correlation functions cannot be accessed., Comment: 17+4 pages, 4+6 figures; v2: title changed, results on different sectors added, figures refined, discussion improved

We consider the interaction-driven Mott transition at zero temperature from the viewpoint of microscopic Fermi liquid theory. To this end, we derive an exact expression for the Landau parameters within the dynamical mean-field theory (DMFT) approximation to the single-band Hubbard model. At the Mott transition the symmetric and the anti-symmetric Landau parameter diverge. The vanishing compressibility at the Mott transition directly implies the divergence of the forward scattering amplitude in the charge sector, which connects the proximity of the Mott phase to a tendency towards phase separation. We verify the expected behavior of the Landau parameters in a DMFT application to the Hubbard model on the triangular lattice at finite temperature. Exact conservation laws and the Ward identity are crucial to capture vertex divergences related to the Mott transition. We furthermore generalize Leggett's formula for the static susceptibility of the Fermi liquid to the static fermion-boson response function. In the charge sector the limits of small transferred momentum and frequency of this response function commute at the Mott transition., Comment: Only 17 pages (+ 6 pages Appendix)

We derive an expression for the effective Josephson coupling from the microscopic Hubbard model. It serves as a starting point for the description of phase fluctuations of local Cooper pairs in $d_{x^2-y^2}$-wave superconductors in the framework of an effective $XY$ model of plaquettes, the Josephson lattice. The expression for the effective interaction is derived by means of the local-force theorem, and it depends on local symmetry-broken correlation functions that we obtain using the cluster dynamical mean-field theory. Moreover, we apply the continuum limit to the Josephson lattice to obtain an expression for the gradient term in the Ginzburg-Landau theory and compare predicted London penetration depths and Kosterlitz-Thouless transition temperatures with experimental data for YBa$_2$Cu$_3$O$_{7-x}$., Comment: 13 pages, 13 figures

We investigate the violation of the first Hund's rule in 4$d$ and 5$d$ transition metal oxides that form solids of dimers. Bonding states within these dimers reduce the magnetization of such materials. We parametrize the dimer formation with realistic hopping parameters and find not only regimes, where the system behaves as a Fermi liquid or as a Peierls insulator, but also strongly correlated regions due to Hund's coupling and its competition with the dimer formation. The electronic structure is investigated using the cluster dynamical mean-field theory for a dimer in the two-plane Bethe lattice with two orbitals per site and $3/8$-filling, that is three electrons per dimer. It reveals dimer-antiferromagnetic order of a high-spin (double exchange) state and a low-spin (molecular orbital) state. At the crossover region we observe the suppression of long-range magnetic order, fluctuation enhancement and renormalization of electron masses. At certain interaction strengths the system becomes an incoherent antiferromagnetic metal with well defined local moments., Comment: 11 pages, 10 figures

In the study of strongly interacting systems, correlations on the two-particle level are receiving more and more attention. In this work, we study a particular two-particle correlation function: the fermion-boson vertex. It describes the response of the Green's function when an external field is applied and is an important ingredient of diagrammatic extensions of Dynamical Mean-Field Theory. We provide several perspectives on this object, using Ward identities, sum rules, perturbative analysis and asymptotic relations. We then use these tools to study the vertex across the doping-driven metal-insulator transition and find a divergence in the imaginary part., Comment: 21 pages

The growing family of two-dimensional (2D) materials that are now available can be used to assemble van der Waals heterostructures with a wide range of properties. Of particular interest are tunnelling heterostructures, which have been used to study the electronic states both in the tunnelling barrier and in the emitter and collector contacts. Recently, 2D ferromagnets have been studied theoretically and experimentally. Here we investigate electron tunnelling through a thin (2-6 layers) ferromagnetic CrBr3 barrier. For devices with non-magnetic barriers, conservation of momentum can be relaxed by phonon-assisted tunnelling or by tunnelling through localised states. In the case of our ferromagnetic barrier the dominant tunnelling mechanisms are the emission of magnons at low temperatures or scattering of electrons on localised magnetic excitations above the Curie temperature. Tunnelling with magnon emission offers the possibility of injecting spin into the collector electrode.

We study the Hubbard model on the honeycomb lattice in the vicinity of the quantum critical point by means of a multiband formulation of the Dual Fermion approach. Beyond the strong local correlations of the dynamical mean field, critical fluctuations on all length scales are included by means of a ladder diagram summation. Analysis of the susceptibility yields an estimate of the critical interaction strength of the quantum phase transition from a paramagnetic semimetal to an antiferromagnetic insulator, in good agreement to other numerical methods. We further estimate the crossover temperature to the renormalized classical regime. Our data imply that, at large interaction strengths, the Hubbard model on the honeycomb lattice behaves like a quantum nonlinear sigma model, while being in a non-Fermi liquid state., Comment: - changed the last sentence in the abstract - rewrote Sec. II.A - added references 11,18 and 29 - changed styles in Fig.4, Fig.5 and Fig.6 to improve coherence - revised Fig.7 promoted insets to individual panels and added polynomial extrapolation of the self-energy data - Rewrote Sec. III.B. - Rewrote parts of the appendix - added Refs. 19-21

We present an efficient exact diagonalization scheme for the extended dynamical mean-field theory and apply it to the extended Hubbard model on the square lattice with nonlocal charge-charge interactions. Our solver reproduces the phase diagram of this approximation with good accuracy. Details on the numerical treatment of the large Hilbert space of the auxiliary Holstein-Anderson impurity problem are provided. Benchmarks with a numerically exact strong-coupling continuous-time quantum-Monte Carlo solver show better convergence behavior of the exact diagonalization in the deep insulator. Special attention is given to possible effects due to the discretization of the bosonic bath. We discuss the quality of real axis spectra and address the question of screening in the Mott insulator within extended dynamical mean-field theory., Comment: 12 pages, 8 figures

The quantum nature of a microscopic system can only be revealed when it is sufficiently decoupled from surroundings. Interactions with the environment induce relaxation and decoherence that turn the quantum state into a classical mixture. Here, we study the timescales of these processes for a qubit encoded in the collective state of a set of magnetic atoms deposited on a metallic surface. For that, we provide a generalization of the commonly used definitions of $T_{1}$ and $T_{2}$ characterizing relaxation and decoherence rates. We calculate these quantities for several atomic structures, including a collective spin, a setup implementing a decoherence-free subspace, and two examples of spin chains. Our work contributes to the comprehensive understanding of the relaxation and decoherence processes and shows the advantages of the implementation of a decoherence free subspace in these setups., Comment: 13 pages, 8 figures

Various intermetallic compounds harbor subtle electronic correlation effects. To elucidate this fact for the Fe-Al system, we perform a realistic many-body investigation based on the combination of density functional theory with dynamical mean-field theory in a charge self-consistent manner. A better characterization and understanding of the phase stability of bcc-based D0$_3$-Fe$_3$Al through an improved description of the correlated charge density and the magnetic-energy is achevied. Upon replacement of one Fe sublattice by V, the Heusler compound Fe$_2$VAl is realized, known to display bad-metal behavior and increased specific heat. We here document a charge-gap opening at low temperatures in line with previous experimental work. The gap structure does not match conventional band theory and is reminiscent of (pseudo)gap charateristics in correlated oxides., Comment: 8 pages, 7 figures

The electronic structure of selected rare-earth atoms adsorbed on a free-standing graphene was investigated using methods beyond the conventional density functional theory (DFT+U, DFT+HIA and DFT+ED). The influence of the electron correlations and the spin-orbit coupling on the magnetic properties has been examined. The DFT+U method predicts both atoms to carry local magnetic moments (spin and orbital) contrary to a nonmagnetic $f^6$ ($J=0$) ground-state configuration of Sm in the gas phase. Application of DFT${}+{}$Hubbard-I (HIA) and DFT${}+{}$exact diagonalization (ED) methods cures this problem, and yields a nonmagnetic ground state with six $f$ electrons and $J=0$ for the Sm adatom. Our calculations show that Nd adatom remains magnetic, with four localized $f$ electrons and $J=4.0$. These conclusions could be verified by STM and XAS experiments., Comment: REVTeX 4-1, 8 pages, 4 figures

We discuss the calculation of the double occupancy using Dynamical Mean-Field Theory (DMFT) in finite dimensions. The double occupancy can be determined from the susceptibility of the auxiliary impurity model or from the lattice susceptibility. The former method underestimates, whereas the latter overestimates the double occupancy. We illustrate this for the square-lattice Hubbard model. We propose an approach for which both methods lead to identical results by construction and which resolves this ambiguity. This self-consistent dual boson scheme results in a double occupancy that is numerically close to benchmarks available in the literature., Comment: 5 pages, 5 figures

We study the dynamical magnetic susceptibility of a strongly correlated electronic system in the presence of a time-dependent hopping field, deriving a generalized Bethe-Salpeter equation which is valid also out of equilibrium. Focusing on the single-orbital Hubbard model within the time-dependent Hartree-Fock approximation, we solve the equation in the non-equilibrium adiabatic regime, obtaining a closed expression for the transverse magnetic susceptibility. From this, we provide a rigorous definition of non-equilibrium (time-dependent) magnon frequencies and exchange parameters, expressed in terms of non-equilibrium single-electron Green functions and self-energies. In the particular case of equilibrium, we recover previously known results., Comment: Revised version

Time-dependent spin phenomena in condensed matter are most often either described in the weakly correlated limit of metallic Stoner/Slater-like magnetism via band theory or in the strongly correlated limit of Heisenberg-like interacting spins in an insulator. However many experimental studies, e.g. of (de)magnetization processes, focus on itinerant local-moment materials such as transition metals and various of their compounds. We here present a general theoretical framework that is capable of addressing correlated spin dynamics, also in the presence of a vanishing charge gap. A real-space implementation of the time-dependent rotational-invariant slave boson methodology allows to treat non-equilibrium spins numerically fast and efficiently beyond linear response as well as beyond the band-theoretical or Heisenberg limit., Comment: 9 pages, 7 figures, revised version

We report on the formation of the high-energy empty-state resonance in the electronic spectrum of the iron adatom on the Pt(111) surface. By using the combination of the first-principles methods and the finite-temperature exact diagonalization approach, we show that the resonance is the result of the valence fluctuations between atomic configurations of the impurity. Our theoretical finding is fully confirmed by the results of the scanning tunneling microscopy measurements [M.F. Crommie et al., Phys. Rev. B 48, 2851 (1993)]. In contrast to the previous theoretical results obtained by using local spin density approximation, the paramagnetic state of the impurity in the experiment is naturally reproduced within our approach. This opens a new way for interpretation of STM data collected earlier for metallic surface nanosystems with iron impurities.

We examine the nature of the transition to the antiferromagnetically ordered state in the half-filled three-dimensional Hubbard model using the dual-fermion multiscale approach. Consistent with analytics, in the weak-coupling regime we find that spin-flip excitations across the Fermi surface are important, and that the strong coupling regime is described by Heisenberg physics. In the intermediate interaction, strong correlation regime we find aspects of both local and non-local correlations. We analyze the critical exponents of the transition in the strong coupling regime and find them to be consistent with Heisenberg physics down to an interaction of $U/t=10$., Comment: 9 pages, 9 figures

We derive a set of equations expressing the parameters of the magnetic interactions characterizing a strongly correlated electronic system in terms of single-electron Green's functions and self-energies. This allows to establish a mapping between the initial electronic system and a spin model including up to quadratic interactions between the effective spins, with a general interaction (exchange) tensor that accounts for anisotropic exchange, Dzyaloshinskii-Moriya interaction and other symmetric terms such as dipole-dipole interaction. We present the formulas in a format that can be used for computations via Dynamical Mean Field Theory algorithms., Comment: Submitted to the Proceedings of the 20th International Conference on Magnetism, Barcelona (ICM2015)

A Mn-porphyrin was contacted on Au(111) in a low-temperature scanning tunneling microscope (STM). Differential conductance spectra show a zero-bias resonance that is due to an underscreened Kondo effect according to many-body calculations. When the Mn center is contacted by the STM tip, the spectrum appears to invert along the voltage axis. A drastic change in the electrostatic potential of the molecule involving a small geometric relaxation is found to cause this observation., Comment: 8 pages, 8 figures (including supplemental material); accepted in Phys. Rev. Lett

We consider the thermodynamic consistency of the charge response function in the (extended) Hubbard model. In DMFT, thermodynamic consistency is preserved. We prove that the static, homogeneous DMFT susceptibility is consistent as long as vertex corrections obtained from the two-particle impurity correlation function are included. In presence of a nonlocal interaction, the problem may be treated within extended DMFT (EDMFT), or its diagrammatic extension, the dual boson approach. We show that here, maintaining thermodynamic consistency requires knowledge of three- and four-particle impurity correlation functions, which are typically neglected. Nevertheless, the dual boson approximation to the response is remarkably close to consistency. This holds even when two-particle vertex corrections are neglected. EDMFT is consistent only in the strongly correlated regime and near half-filling, where the physics is predominantly local., Comment: 11 pages (incl. appendix), 4 figures

The local-density approximation combined with the dynamical mean-field theory (LDA+DMFT) is applied to the paramagnetic phase of light actinide dioxides: UO2, NpO2, and PuO2. The calculated band gaps and the valence-band electronic structure are in a very good agreement with the optical absorption experiments as well as with the photoemission spectra. The hybridization of the actinide 5f shell with the 2p states of oxygen is found to be relatively large, it increases the filling of the 5f orbitals from the nominal ionic configurations with two, three, and four electrons to nearly half-integer values 2.5, 3.4 and 4.4. The large hybridization leaves an imprint also on the core-level photoemission spectra in the form of satellite peaks. It is demonstrated that these satellites are accurately reproduced by the LDA+DMFT calculations., Comment: REVTeX 4-1, 10 pages, 6 figures

We present a technique to map an electronic model with local interactions (a generalized multi-orbital Hubbard model) onto an effective model of interacting classical spins, by requiring that the thermodynamic potentials associated to spin rotations in the two systems are equivalent up to second order in the rotation angles. This allows to determine the parameters of relativistic and non-relativistic magnetic interactions in the effective spin model in terms of equilibrium Green's functions of the electronic model. The Hamiltonian of the electronic system includes, in addition to the non-relativistic part, relativistic single-particle terms such as the Zeeman coupling to an external magnetic fields, spin-orbit coupling, and arbitrary magnetic anisotropies; the orbital degrees of freedom of the electrons are explicitly taken into account. We determine the complete relativistic exchange tensors, accounting for anisotropic exchange, Dzyaloshinskii-Moriya interactions, as well as additional non-diagonal symmetric terms (which may include dipole-dipole interaction). Our procedure provides the complete exchange tensors in a unified framework, including previously disregarded features such as the vertices of two-particle Green's functions and non-local self-energies. We do not assume any smallness in spin-orbit coupling, so our treatment is in this sense exact. Finally, we show how to distinguish and address separately the spin, orbital and spin-orbital contributions to magnetism., Comment: Revised version, including a discussion of the role of symmetry breaking and complete explicit formulas for spin, orbital and spin-orbital parts of the magnetic interactions, both including and excluding vertices of the two-particle Green's functions

The dual-fermion approach offers a way to perform diagrammatic expansion around the dynamical mean-field theory. Using this formalism, the influence of antiferromagnetic fluctuations on the self-energy is taken into account through ladder-type diagrams in the particle-hole channel. The resulting phase diagram for the (quasi-)two-dimensional Hubbard model exhibits antiferromagnetism and d-wave superconductivity. Furthermore, a uniform charge instability, i.e., phase separation, is obtained in the low doping regime around the Mott insulator. We also examine spin/charge density wave fluctuations including d-wave symmetry. The model exhibits a tendency towards an unconventional charge density-wave, but no divergence of the susceptibility is found., Comment: 14 pages, 13 figures, 2 tables

We investigate the excitation spectrum of strongly correlated sodium cobaltate within a realistic many-body description beyond dynamical mean-field theory (DMFT). At lower doping around $x$=0.3, rather close to Mott-critical half-filling, the single-particle spectral function of Na$_x$CoO$_2$ displays an upper Hubbard band which is captured within DMFT. Momentum-dependent self-energy effects beyond DMFT become dominant at higher doping. Around a doping level of $x\sim 0.67$, the incoherent excitations give way to finite-energy spin-polaron excitations in close agreement with optics experiments. These excitations are a direct consequence of the formation of bound states between quasiparticles and paramagnons in the proximity to in-plane ferromagnetic ordering., Comment: 5 pages, 3 figures; supplementary material

The dual boson approach [Ann. Phys. 327, 1320 (2012)] provides a means to construct a diagrammatic expansion around the extended dynamical mean-field theory (EDMFT). In this paper, we present the numerical implementation of the approach and apply it to the extended Hubbard model with nearest-neighbor interaction $V$. We calculate the EDMFT phase diagram and study the effect of diagrams beyond EDMFT on the transition to the charge-ordered phase. Including diagrammatic corrections to the EDMFT polarization shifts the EDMFT phase boundary to lower values of $V$. The approach interpolates between the random phase approximation in the weak coupling limit and EDMFT for strong coupling. Neglecting vertex corrections leads to results reminiscent of the EDMFT+$GW$ approximation. We however find significant deviations from the dual boson results already for small values of the interaction, emphasizing the crucial importance of fermion-boson vertex corrections., Comment: Published version; 24 pages, 32 figures

We present a first-principles many-body analysis of multi-orbital lattice susceptibilities in the metallic phase of the quasi-twodimensional compound LiVS$_2$. We base this on advanced correlated electronic structure methods for the $t_{2g}$-shell to reveal a highly entangled spin-orbital hexagonal ordering (SOHO) bringing about an inherently intersite order parameter for the trimerization transition, which eventually leads to an intriguing insulating phase at low temperature.

We consider the collective, long-wavelength charge excitations in correlated media in presence of short- and long-range forces. As an example for the case of a short-range interaction, we examine the two-dimensional Hubbard model within dynamical mean-field theory (DMFT). It is shown explicitly that the DMFT susceptibility including vertex corrections respects the Ward identity and yields a manifestly gauge-invariant response in finite dimensions. For computing the susceptibility, we use a different expression and establish its formal equivalence to the standard DMFT formula. It allows for a more stable analytical continuation. We find a zero-sound mode expected for short-range forces. The relation between the vertex corrections, gauge invariance, and the appearance of the collective modes is discussed. Long-range forces are treated within extended dynamical mean-field theory. In order to obtain a gauge-invariant response, it is necessary to additionally incorporate some non local vertex corrections into the polarization. In doing so, we obtain plasmons in the three-dimensional Hubbard model. The plasma frequency is determined by the (single-particle) density distribution as a consequence of gauge invariance. We compare this result with the plasma frequency extracted from the analytical continuation of the susceptibility. It is in good agreement with the prediction from the gauge-invariance condition., Comment: Published version; 21 pages, 22 figures