Heteroepitaxial Control of Fermi Liquid, Hund Metal, and Mott Insulator Phases in Single‐Atomic‐Layer Ruthenates
Not Available
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Not Available
The nonequilibrium dynamics of correlated many-particle systems is of interest in connection with pump–probe experiments on molecular systems and solids, as well as theoretical investigations of transport properties and relaxation processes. Nonequilibrium Green’s functions are a powerful tool to study interaction effects in quantum many-particle systems out of equilibrium, and to extract physically relevant information for the interpretation of experiments. Here, we present the open-source software package NESSi (The Non-Equilibrium Systems Simulation package) which allows to perform many-body dynamics simulations based on Green’s functions on the L-shaped Kadanoff–Baym contour. NESSi contains the library libcntr which implements tools for basic operations on these nonequilibrium Green’s functions, for constructing Feynman diagrams, and for the solution of integral and integro-differential equations involving contour Green’s functions. The library employs a discretization of the Kadanoff–Baym contour into time points and a high-order implementation of integration routines. The total integrated error scales up to $\mathcal{O}(N^{-7})$, which is important since the numerical effort increases at least cubically with the simulation time. A distributed-memory parallelization over reciprocal space allows large-scale simulations of lattice systems. We provide a collection of example programs ranging from dynamics in simple two-level systems to problems relevant in contemporary condensed matter physics, including Hubbard clusters and Hubbard or Holstein lattice models. The libcntr library is the basis of a follow-up software package for nonequilibrium dynamical mean-field theory calculations based on strong-coupling perturbative impurity solvers.
Not Available
Abstract Multiorbital systems away from global half-filling host intriguing physical properties promoted by Hund’s coupling. Despite increasing awareness of this regime dubbed Hund’s metal, effect of nonlocal interaction is still elusive. Here we study a three-orbital model with 1/3 filling (two electrons per site) including the intersite Coulomb interaction ( V ). Using the G W plus extended dynamical mean-field theory, the valence-skipping charge order transition is shown to be driven by V . Most interestingly, the instability to this transition is significantly enhanced in the spin-freezing crossover regime, thereby lowering the critical V to the formation of charge order. This behavior is found to be closely related to the population profile of the atomic multiplet states in the spin-freezing regime. In this regime, maximum spin states are dominant in each total charge subspace with substantial amount of one- and three-electron occupations, which leads to almost equal population of one- and the maximum spin three-electron state. Our finding unveils another feature of the Hund’s metal and has potential implications for the broad range of multiorbital systems as well as the recently discovered charge order in iron pnictides.
Abstract The electronic property and magnetic susceptibility of Ce 3 Pd 3 Bi 4 were systemically investigated from 18 to 290 K for varying values of cell volume using dynamic mean-field theory coupled with density functional theory. By extrapolating to zero temperature, the ground state of Ce 3 Pd 3 Bi 4 at ambient pressure is found to be a correlated semimetal due to insufficient hybridization. Upon applying pressure, the hybridization strength increases and a crossover to the Kondo insulator is observed at finite temperatures. The characteristic temperature signaling the formation of Kondo singlet, as well as the characteristic temperature associated with f -electron delocalization–localization change, simultaneously vanishes around a critical volume of 0.992 V 0 , suggesting that such metal–insulator transition is possibly associated with a quantum critical point. Finally, Wilson’s loop calculations indicate that the Kondo insulating side is topologically trivial, thus a topological transition also occurs across the quantum critical point.
Abstract The discovery of two-dimensional (2D) van der Waals (vdW) materials often provides interesting playgrounds to explore novel phenomena. One of the missing components in 2D vdW materials is the intrinsic heavy-fermion systems, which can provide an additional degree of freedom to study quantum critical point (QCP), unconventional superconductivity, and emergent phenomena in vdW heterostructures. Here, we investigate 2D vdW heavy-fermion candidates through the database of experimentally known compounds based on dynamical mean-field theory calculation combined with density functional theory (DFT+DMFT). We have found that the Kondo resonance state of CeSiI does not change upon exfoliation and can be easily controlled by strain and surface doping. Our result indicates that CeSiI is an ideal 2D vdW heavy-fermion material and the quantum critical point can be identified by external perturbations.
The discovery of superconductivity in hole-doped infinite layer nickelates, RNiO 2 (R = Nd, Pr, La) has resulted in sustained interest in the field. A definitive picture of low-energy many-body states has not yet emerged. Here we provide insights into the low-energy physics, based on our embedded dynamical mean-field theory calculations, and propose a correlation (U)-temperature (T) phase diagram. The key features are a low-T Fermi liquid (FL) phase, a high-T Curie-Weiss regime, and an antiferromagnetic phase in a narrow U-T region. We associate the onset of the FL phase with partial screening of Ni-d moments; however, full screening occurs at lower temperatures. This may be related to insufficiency of conduction electrons to effectively screen the Ni-d moments, suggestive of Nozieres Exhaustion Principle. Our results suggest that RNiO 2 are in the paramagnetic state, close to an antiferromagnetic dome, making magnetic fluctuations feasible. This may be consequential for superconductivity.
The intricate nature of magnetism in uranium-based Kondo lattices is a consequence of correlations between U-5 f and conduction electrons. Previously, the source of magnetism has been ascribed to either Mott physics or Ruderman-Kittel-Kasuya-Yosida interaction, both of which are not fully applicable to uranium-based Kondo lattices. Using linearized quasiparticle self-consistent GW plus dynamical mean-field theory, we demonstrate a crossover from incoherent to coherent f - d Kondo cloud in the paramagnetic phase of UTe 2 , USbTe and USbSe. As the transition occurs, we observe an augmented f - d coherence and Pauli-like magnetic susceptibility, with a substantial frozen magnetic moment of U-5 f persisting. We show that momentum dependent f - d hybridization is responsible for the magnetic moments arising from the renormalized f electrons’ van Hove singularity. Our findings provide a perspective to explain the dual nature of magnetism and the long-range magnetic ordering induced by pressure in UTe 2 .
We study the fluctuations responsible for pairing in the d -wave superconducting state of the two-dimensional Hubbard model at intermediate coupling within a cluster dynamical mean-field theory with a numerically exact quantum impurity solver. By analyzing how momentum- and frequency-dependent fluctuations generate the d -wave superconducting state in different representations, we identify antiferromagnetic fluctuations as the pairing glue of superconductivity in both the underdoped and the overdoped regime. Nevertheless, in the intermediate coupling regime, the predominant magnetic fluctuations may differ significantly from those described by conventional spin fluctuation theory.
Electronic flat bands associated with quenched kinetic energy and heavy electron mass have attracted great interest for promoting strong electronic correlations and emergent phenomena such as high-temperature charge fractionalization and superconductivity. Intense experimental and theoretical research has been devoted to establishing the rich nontrivial metallic and heavy fermion phases intertwined with such localized electronic states. Here, in this study, we investigate the transition metal oxide spinel LiV 2 O 4 , an enigmatic heavy fermion compound lacking localized f orbital states. We use angle-resolved photoemission spectroscopy and dynamical mean-field theory to reveal a kind of correlation-induced flat band with suppressed interatomic electron hopping arising from intra-atomic Hund’s coupling. The appearance of heavy quasiparticles is ascribed to a proximate orbital-selective Mott state characterized by fluctuating local moments as evidenced by complementary magnetotransport measurements. The spectroscopic fingerprints of long-lived quasiparticles and their disappearance with increasing temperature further support the emergence of a high-temperature “bad” metal state observed in transport data. This work resolves a long-standing puzzle on the origin of heavy fermion behavior and unconventional transport in LiV 2 O 4 . Simultaneously, it opens a path to achieving flat bands through electronic interactions in d-orbital systems with geometrical frustration, potentially enabling the realization of exotic phases of matter such as the fractionalized Fermi liquids.
Quantum embedding methods, such as dynamical mean-field theory (DMFT), provide a powerful framework for investigating strongly correlated materials. A central computational bottleneck in DMFT is in solving the Anderson impurity model (AIM), whose exact solution is classically intractable for large bath sizes. In this work, we benchmark a quantum-classical hybrid solver tailored for particle-hole symmetric AIMs, using the variational quantum eigensolver to prepare the ground state of the model with shallow quantum circuits. The solver uses shallow quantum ansätze and one set of variational parameters to prepare the ground state and its particle and hole excitations, enabling the construction of the impurity Green’s function through a continued-fraction expansion. We evaluate the performance of this approach across a few bath sizes and interaction strengths under noisy, shot-limited conditions. We compare three optimization routines (COBYLA, Adam, and L-BFGS-B) in terms of convergence and fidelity, assess the benefits of estimating a quantum-computed moment correction to the variational energies, and benchmark the approach by comparing the density of states computed from the impurity Green’s function against that obtained using a classical pipeline. Our results demonstrate the feasibility of Green’s function construction on near-term devices and establish practical benchmarks for quantum impurity solvers embedded within self-consistent DMFT loops.
Rutile ( R ) phase VO 2 is a quintessential example of a strongly correlated bad metal, which undergoes a metal-insulator transition (MIT) concomitant with a structural transition to a V-V dimerized monoclinic ( M 1 ) phase below T MIT ~ 340 K . It has been experimentally illustrated that one can control this transition by doping VO 2 . In particular, doping with oxygen vacancies ( V O ) has been shown to completely suppress this MIT without any structural transition. We explain this suppression by elucidating the influence of oxygen vacancies on the electronic structure of the metallic R phase VO 2 , explicitly treating strong electron-electron correlations using dynamical mean-field theory (DMFT) as well as diffusion Monte Carlo (DMC) flavor of quantum Monte Carlo (QMC) techniques. DMC calculations show a gap closure in the M 1 phase when vacancies are present, suggesting that when vacancies are introduced in the high-temperature rutile phase, the dimerized insulating phase cannot be reached when temperature is lowered. Both DMFT and DMC calculations of nonstoichiometric metallic rutile phase shows that this tendency not to dimerize in the presence of vacancies is because V O 's tend to change the V - 3 d filling away from its nominal half-filled value, with the e g π orbitals competing with the otherwise dominant a 1 g orbital. Loss of this near orbital polarization of the a 1 g orbital is associated with a weakening of electron correlations, especially along the V-V dimerization direction. This removes a charge-density wave (CDW) instability along this direction above a critical doping concentration, which further suppresses the metal-insulator transition. Our research further suggests that the MIT is predominantly driven by a correlation-induced CDW instability along the V-V dimerization direction.
Despite decades of efforts, the origin of metal-insulator transitions (MITs) in strongly correlated materials remains one of the main long-standing problems in condensed-matter physics. An archetypal example is V 2 O 3 , which undergoes simultaneous electronic, structural, and magnetic phase transitions. This remarkable feature highlights the many degrees of freedom at play in this material. In this work, acting solely on the magnetic degree of freedom, we reveal an anomalous feature in the electronic transport of V 2 O 3 : On cooling, the magnetoresistance changes from positive to negative values well above the MIT temperature, and shows divergent behavior at the transition. In this work, the effects are attributed to the magnetic field quenching antiferromagnetic fluctuations above the Néel temperature $T_N$, and preventing long-range antiferromagnetic ordering below $T_N$. In both cases, suppressing the antiferromagnetic ordering prevents the opening of the incipient electronic gap. This interpretation is supported by Hubbard model calculations which fully reproduce the experimental behavior. Our study sheds light on this classic problem providing a clear and physical interpretation of the nature of the metal-insulator transition.
Plutonium is a critically important material as the behavior of its 5f electrons stands midway between the metalliclike itinerant character of the light actinides and localized atomic-core-like character of the heavy actinides. The δ phase of plutonium (δ-Pu), whereas still itinerant, has a large coherent Kondo peak and strong electronic correlations coming from its near-localized character. Here, using sophisticated Gutzwiller wave function and dynamical mean-field theory correlated theories, we study the Fermi surface and associated mass renormalizations of δ-Pu together with calculations of the de Haas–van Alphen frequencies. We find a large (~200%) correlation induced volume expansion in both the hole and the electron pockets of the Fermi surface in addition to an intermediate mass enhancement. All of the correlated electron theories predict, approximately, the same hole pocket placement in the Brillouin zone, which is different from that obtained in conventional density-functional band-structure theory, whereas the electron pockets from all theories are in, roughly, the same place.
We combine transport, angle-resolved photoemission spectroscopy, and scanning tunneling spectroscopy to investigate several low-energy manifestations of the Hund's coupling in a canonical FeSC family Li(Fe,Co)As. We determine the doping dependence of the coherent-incoherent crossover temperature and the quasiparticle effective mass enhancement in the normal state. Furthermore, our tunneling spectroscopy result in the superconducting state supports the idea that superconductivity emerging from Hund's metal state displays a universal maximal superconducting gap vs transition temperature (2Δ max /k B T c ) value, which is independent of doping level and T c .
The dimer Mott insulator kappa - (BEDT-TTF) 2 - Cu 2 (CN) 3 can be tuned into metallic and superconducting states on applying pressure of 1.5 kbar and more. We have performed dielectric measurements (7.5 kHz to 5 MHz) on kappa - (BEDT-TTF) 2 - Cu 2 (CN) 3 single crystals as a function of temperature (down to T = 8 K) and pressure (up to p = 4.3 kbar). In addition to the relaxor-like dielectric behavior seen below 50 K at p = 0, that moves toward lower temperatures with pressure, a second peak emerges in ε 1 (T) around T = 15 K. When approaching the insulator-metal boundary, this peak diverges rapidly reaching ε 1 ≈ 10 5 . Our dynamical mean-field theory calculations substantiate that the dielectric catastrophe at the Mott transition is not caused by closing the energy gap, but due to the spatial coexistence of correlated metallic and insulating regions. We discuss the percolative nature of the first-order Mott insulator-to-metal transition in all details.
The Anderson Impurity Model (AIM) is a canonical model of quantum many-body physics. Here we investigate whether machine learning models, both neural networks (NN) and kernel ridge regression (KRR), can accurately predict the AIM spectral function in all of its regimes, from empty orbital, to mixed valence, to Kondo. To tackle this question, we construct two large spectral databases containing approximately 410 000 and 600 000 spectral functions of the single-channel impurity problem. We show that the NN models can accurately predict the AIM spectral function in all of its regimes, with pointwise mean absolute errors down to 0.003 in normalized units. We find that the trained NN models outperform models based on KRR and enjoy a speedup on the order of 10 5 over traditional AIM solvers. Finally, the required size of the training set of our model can be significantly reduced using farthest point sampling in the AIM parameter space, which is important for generalizing our method to more complicated multichannel impurity problems of relevance to predicting the properties of real materials.
We generalize the formalism of the dynamical vertex approximation (DΓA)—a diagrammatic extension of the dynamical mean-field theory (DMFT)—to treat magnetically ordered phases. To this aim, we start by concisely illustrating the many-electron formalism for performing ladder resummations of Feynman diagrams in systems with broken SU(2) symmetry associated to ferromagnetic (FM) or antiferromagnetic (AF) order. We then analyze the algorithmic simplifications introduced by taking the local approximation of the two-particle irreducible vertex functions in the Bethe-Salpeter equations, which defines the ladder implementation of DΓA for magnetic systems. The relation of this assumption with the DMFT limit of large coordination-number/high dimensions is explicitly discussed. As a last step, we derive the expression for the ladder DΓA self-energy in the FM- and AF-ordered phases of the Hubbard model. The physics emerging in the AF-ordered case is explicitly illustrated by means of approximated calculations based on a static mean-field input for DΓA equations. The results obtained capture fundamental aspects of both metallic and insulating ground states of two-dimensional antiferromagnets, providing a reliable compass for future, more extensive applications of our approach. Furthermore, possible routes to further develop diagrammatic-based treatments of magnetic phases in correlated electron systems are briefly outlined in the Conclusions.