BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Emanuel Gull

Emanuel Gull contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.str-elphysics.comp-phcond-mat.mtrl-scicond-mat.supr-conphysics.chem-phcond-mat.dis-nncond-mat.mes-hallcond-mat.stat-mechquant-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
22works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

22 published item(s)

preprint2023arXiv

Stark-Many body localization in interacting infinite dimensional systems

We study bulk particle transport in a Fermi-Hubbard model on an infinite-dimensional Bethe lattice, driven by a constant electric field. Previous numerical studies showed that one dimensional analogs of this system exhibit a breakdown of diffusion due to Stark many-body localization (Stark-MBL) at least up to time which scales exponentially with the system size. Here, we consider systems initially in a spin density wave state using a combination of numerically exact and approximate techniques. We show that for sufficiently weak electric fields, the wave's momentum component decays exponentially with time in a way consistent with normal diffusion. By studying different wavelengths, we extract the dynamical exponent and the generalized diffusion coefficient at each field strength. Interestingly, we find a non-monotonic dependence of the dynamical exponent on the electric field. As the field increases towards a critical value proportional to the Hubbard interaction strength, transport slows down, becoming sub-diffusive. At large interaction strengths, however, transport speeds up again with increasing field, exhibiting super-diffusive characteristics when the electric field is comparable to the interaction strength. Eventually, at the large field limit, localization occurs and the current through the system is suppressed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Interaction expansion inchworm Monte Carlo solver for lattice and impurity models

Multi-orbital quantum impurity models with general interaction and hybridization terms appear in a wide range of applications including embedding, quantum transport, and nanoscience. However, most quantum impurity solvers are restricted to a few impurity orbitals, discretized baths, diagonal hybridizations, or density-density interactions. Here, we generalize the inchworm quantum Monte Carlo method to the interaction expansion and explore its application to typical single- and multi-orbital problems encountered in investigations of impurity and lattice models. Our implementation generically outperforms bare and bold-line quantum Monte Carlo algorithms in the interaction expansion. So far, for the systems studied here, it remains inferior to the more specialized hybridization expansion and auxiliary field algorithms. The problem of convergence to unphysical fixed points, which hampers so-called bold-line methods, is not encountered in inchworm Monte Carlo.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mechanism of Superconductivity in the Hubbard Model at Intermediate Interaction Strength

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 both in 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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantifying the role of antiferromagnetic fluctuations in the superconductivity of the doped Hubbard model

We study the contribution of the electron-spin fluctuation coupling to the superconducting state of the two dimensional Hubbard model within dynamical cluster approximation (DCA) using a numerical exact continuous time Monte Carlo solver. By analyzing the frequency dependence of the self energy, we show that only about half of the superconductivity can be attributed to a "pairing glue" arising from treating spin fluctuations as a pairing boson in the standard one-loop theory.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Reduced Dynamics of Full Counting Statistics

We present a theory of modified reduced dynamics in the presence of counting fields. Reduced dynamics techniques are useful for describing open quantum systems at long emergent timescales when the memory timescales are short. However, they can be difficult to formulate for observables spanning the system and its environment, such as those characterizing transport properties. A large variety of mixed system--environment observables, as well as their statistical properties, can be evaluated by considering counting fields. Given a numerical method able to simulate the field-modified dynamics over the memory timescale, we show that the long-lived full counting statistics can be efficiently obtained from the reduced dynamics. We demonstrate the utility of the technique by computing the long-time current in the nonequilibrium Anderson impurity model from short-time Monte Carlo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Relativistic Self-Consistent $GW$: Exact Two-Component Formalism with One-Electron Approximation for Solids

We present a formulation of relativistic self-consistent $GW$ for solids based on the exact two-component formalism with one-electron approximation (X2C1e) and non-relativistic Coulomb interactions. Our theory allows us to study scalar relativistic effects, spin-orbit coupling, and the interplay of relativistic effects with electron correlation without adjustable parameters. Our all-electron implementation is fully $ab$ $initio$ and does not require a pseudopotential constructed from atomic calculations. We examine the effect of the X2C1e approximation by comparison to the established four-component formalism and reach excellent agreement. The simplicity of X2C1e enables the construction of higher order theories, such as embedding theories, on top of perturbative calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Phase transitions in partial summation methods: Results from the 3D Hubbard model

The accurate determination of magnetic phase transitions in electronic systems is an important task of solid state theory. While numerically exact results are readily available for model systems such as the half-filled 3D Hubbard model, the complexity of real materials requires additional approximations, such as the restriction to certain classes of diagrams in perturbation theory, that reduce the precision with which magnetic properties are described. In this work, we examine the description of magnetic properties in second order perturbation theory, GW, FLEX, and two TMatrix approximations to numerically exact CT-QMC reference data. We assess finite-size effects and compare periodic lattice simulations to cluster embedding. We find that embedding substantially improves finite size convergence. However, by analyzing different partial summation methods we find no systematic improvement in the description of magnetic properties, with most methods considered in this work predicting first-order instead of continuous transitions, leading us to the conclusion that non-perturbative methods are necessary for the accurate determination of magnetic properties and phase transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Hubbard model: A computational perspective

The Hubbard model is the simplest model of interacting fermions on a lattice and is of similar importance to correlated electron physics as the Ising model is to statistical mechanics or the fruit fly to biomedical science. Despite its simplicity, the model exhibits an incredible wealth of phases, phase transitions, and exotic correlation phenomena. While analytical methods have provided a qualitative description of the model in certain limits, numerical tools have shown impressive progress in achieving quantitative accurate results over the last years. This article gives an introduction to the model, motivates common questions, and illustrates the progress that has been achieved over the last years in revealing various aspects of the correlation physics of the model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

$T$-linear resistivity in models with local self-energy

A theoretical understanding of the enigmatic linear-in-temperature ($T$) resistivity, ubiquitous in strongly correlated metallic systems, has been a long sought-after goal. Furthermore, the slope of this robust $T$-linear resistivity is also observed to stay constant through crossovers between different temperature regimes: a phenomenon we dub "slope invariance". Recently, several solvable models with $T$-linear resistivity have been proposed, putting us in an opportune moment to compare their inner workings in various explicit calculations. We consider two strongly correlated models with local self-energies that demonstrate $T$-linearity: a lattice of coupled Sachdev-Ye-Kitaev (SYK) models and the Hubbard model in single-site dynamical mean-field theory (DMFT). We find that the two models achieve $T$-linearity through distinct mechanisms at intermediate temperatures. However, we also find that these mechanisms converge to an identical form at high temperatures. Surprisingly, both models exhibit "slope invariance" across the two temperature regimes. We thus not only reveal some of the diversity in the theoretical inner workings that can lead to $T$-linear resistivity, but we also establish that different mechanisms can result in "slope invarance".

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Ab-Initio self-energy embedding for the photoemission spectra of NiO and MnO

The accurate ab-initio simulation of periodic solids with strong correlations is one of the grand challenges of condensed matter. While mature methods exist for weakly correlated solids, the ab-initio description of strongly correlated systems is an active field of research. In this work, we show results for the single particle spectral function of the two correlated $d$-electron solids NiO and MnO from self-energy embedding theory. Unlike earlier work, the theory does not use any adjustable parameters and is fully ab-initio, while being able to treat both the strong correlation and the non-local screening physics of these materials. We derive the method, discuss aspects of the embedding and choices of physically important orbitals, and compare our results to x-ray and angle-resolved photoemission spectroscopy as well as bremsstrahlung-isochromat spectroscopy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Diagrammatic Monte Carlo Method for Impurity Models with General Interactions and Hybridizations

We present a diagrammatic Monte Carlo method for quantum impurity problems with general interactions and general hybridization functions. Our method uses a recursive determinant scheme to sample diagrams for the scattering amplitude. Unlike in other methods for general impurity problems, an approximation of the continuous hybridization function by a finite number of bath states is not needed, and accessing low temperature does not incur an exponential cost. We test the method for the example of molecular systems, where we systematically vary temperature, interatomic distance, and basis set size. We further apply the method to an impurity problem generated by a self-energy embedding calculation of correlated antiferromagnetic NiO. We find that the method is ideal for quantum impurity problems with a large number of orbitals but only moderate correlations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Error propagation in the fully self-consistent stochastic second-order Green's function method

We present an implementation of a fully self-consistent finite temperature second order Green's function perturbation theory (GF2) within the diagrammatic Monte Carlo framework. In contrast to the previous implementations of stochastic GF2 ({\it J. Chem. Phys.},{\bf 151}, 044144 (2019)), the current self-consistent stochastic GF2 does not introduce a systematic bias of the resulting electronic energies. Instead, the introduced implementation accounts for the stochastic errors appearing during the solution of the Dyson equation. We present an extensive discussion of the error handling necessary in a self-consistent procedure resulting in dressed Green's function lines. We test our method on a series of simple molecular examples.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Legendre-spectral Dyson equation solver with super-exponential convergence

Quantum many-body systems in thermal equilibrium can be described by the imaginary time Green's function formalism. However, the treatment of large molecular or solid ab inito problems with a fully realistic Hamiltonian in large basis sets is hampered by the storage of the Green's function and the precision of the solution of the Dyson equation. We present a Legendre-spectral algorithm for solving the Dyson equation that addresses both of these issues. By formulating the algorithm in Legendre coefficient space, our method inherits the known faster-than-exponential convergence of the Green's function's Legendre series expansion. In this basis, the fast recursive method for Legendre polynomial convolution, enables us to develop a Dyson equation solver with quadratic scaling. We present benchmarks of the algorithm by computing the dissociation energy of the helium dimer He$_2$ within dressed second-order perturbation theory. For this system, the application of the Legendre spectral algorithm allows us to achieve an energy accuracy of $10^{-9} E_h$ with only a few hundred expansion coefficients.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multiorbital Quantum Impurity Solver for General Interactions and Hybridizations

We present a numerically exact Inchworm Monte Carlo method for equilibrium multiorbital quantum impurity problems with general interactions and hybridizations. We show that the method, originally developed to overcome the dynamical sign problem in certain real-time propagation problems, can also overcome the sign problem as a function of temperature for equilibrium quantum impurity models. This is shown in several cases where the current method of choice, the continuous-time hybridization expansion, fails due to the sign problem. Our method therefore enables simulations of impurity problems as they appear in embedding theories without further approximations, such as the truncation of the hybridization or interaction structure or a discretization of the impurity bath with a set of discrete energy levels, and eliminates a crucial bottleneck in the simulation of ab initio embedding problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Nevanlinna Analytical Continuation

Simulations of finite temperature quantum systems provide imaginary frequency Green's functions that correspond one-to-one to experimentally measurable real-frequency spectral functions. However, due to the bad conditioning of the continuation transform from imaginary to real frequencies, established methods tend to either wash out spectral features at high frequencies or produce spectral functions with unphysical negative parts. Here, we show that explicitly respecting the analytic `Nevanlinna' structure of the Green's function leads to intrinsically positive and normalized spectral functions, and we present a continued fraction expansion that yields all possible functions consistent with the analytic structure. Application to synthetic trial data shows that sharp, smooth, and multi-peak data is resolved accurately. Application to the band structure of silicon demonstrates that high energy features are resolved precisely. Continuations in a realistic correlated setup reveal additional features that were previously unresolved. By substantially increasing the resolution of real frequency calculations our work overcomes one of the main limitations of finite-temperature quantum simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Short-range charge fluctuations in the two-dimensional Hubbard model

We investigate charge fluctuations in the two-dimensional Hubbard model as a function of doping, interaction strength, next-nearest-neighbor hopping, and temperature within the eight-site dynamical cluster approximation. In the regime of intermediate interaction strengths, we find that $d$ density wave fluctuations, which had previously been postulated using analytical arguments, are present and strong but cannot be interpreted as the cause of the pseudogap in the model, due to their evolution with doping and interaction strength. For all parameters away from half filling, the charge fluctuations investigated, including $d$ density wave fluctuations, are weaker than $d$-wave superconducting fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Direct comparison of many-body methods for realistic electronic Hamiltonians

A large collaboration carefully benchmarks 20 first principles many-body electronic structure methods on a test set of 7 transition metal atoms, and their ions and monoxides. Good agreement is attained between the 3 systematically converged methods, resulting in experiment-free reference values. These reference values are used to assess the accuracy of modern emerging and scalable approaches to the many-electron problem. The most accurate methods obtain energies indistinguishable from experimental results, with the agreement mainly limited by the experimental uncertainties. Comparison between methods enables a unique perspective on calculations of many-body systems of electrons.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Ferromagnetic spin correlations in the two-dimensional Hubbard model

We analyze the dynamical nearest-neighbor and next-nearest-neighbor spin correlations in the 4-site and 8-site dynamical cluster approximation to the two-dimensional Hubbard model. Focusing on the robustness of these correlations at long imaginary times, we reveal enhanced ferromagnetic correlations on the lattice diagonal, consistent with the emergence of composite spin-1 moments at a temperature scale that essentially coincides with the pseudo-gap temperature $T^*$. We discuss these results in the context of the spin-freezing theory of unconventional superconductivity.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Magnetic and charge susceptibilities in the half-filled triangular lattice Hubbard model

We study magnetic and charge susceptibilities in the half-filled two-dimensional triangular Hubbard model within the dual fermion approximation in the metallic, Mott insulating, and crossover regions of parameter space. In the \textcolor{black}{insulating state}, we find strong spin fluctuations at the K point at low energy corresponding to the \textcolor{black}{120$^{\circ}$} antiferromagnetic order. These spin fluctuations persist into the metallic phase and move to higher energy. We also present data for simulated neutron spectroscopy and \textcolor{black}{spin-lattice} relaxation times, and perform direct comparisons to inelastic neutron spectroscopy experiments on the triangular material Ba$_8$CoNb$_6$O$_{24}$ and to the relaxation times on $κ$-(ET)$_2$Cu$_2$(CN)$_3$. Finally, we present charge susceptibilities in different areas of parameter space, which should correspond to momentum-resolved electron-loss spectroscopy measurements on triangular compounds.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Sparse sampling and tensor network representation of two-particle Green's functions

Many-body calculations at the two-particle level require a compact representation of two-particle Green's functions. In this paper, we introduce a sparse sampling scheme in the Matsubara frequency domain as well as a tensor network representation for two-particle Green's functions. The sparse sampling is based on the intermediate representation basis and allows an accurate extraction of the generalized susceptibility from a reduced set of Matsubara frequencies. The tensor network representation provides a system independent way to compress the information carried by two-particle Green's functions. We demonstrate efficiency of the present scheme for calculations of static and dynamic susceptibilities in single- and two-band Hubbard models in the framework of dynamical mean-field theory.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Sparse sampling approach to efficient ab initio calculations at finite temperature

Efficient ab initio calculations of correlated materials at finite temperature require compact representations of the Green's functions both in imaginary time and Matsubara frequency. In this paper, we introduce a general procedure which generates sparse sampling points in time and frequency from compact orthogonal basis representations, such as Chebyshev polynomials and intermediate representation (IR) basis functions. These sampling points accurately resolve the information contained in the Green's function, and efficient transforms between different representations are formulated with minimal loss of information. As a demonstration, we apply the sparse sampling scheme to diagrammatic $GW$ and GF2 calculations of a hydrogen chain, of noble gas atoms and of a silicon crystal.

0 citations0 reviewsNo code linkedOpen access