BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Emanuel Gull

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

ResearcherUnclaimed source record
cond-mat.str-elphysics.comp-phcond-mat.supr-concond-mat.mtrl-scicond-mat.quant-gascond-mat.mes-hallphysics.chem-phcond-mat.stat-mechcond-mat.dis-nnquant-ph

Catalog footprint

What is connected

65works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

65 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
preprint2016arXiv

Diagrammatic Monte Carlo for Dual Fermions

We introduce a numerical algorithm to stochastically sample the dual fermion perturbation series around the dynamical mean field theory, generating all topologies of two-particle interaction vertices. We show results in the weak and strong coupling regime of the half-filled Hubbard model in two dimensions, illustrating that the method converges quickly where dynamical mean field theory is a good approximation, and show that corrections are large in the strong correlation regime at intermediate interaction. The fast convergence of dual corrections to dynamical mean field results illustrates the power of the approach and opens a practical avenue towards the systematic inclusion of non-local correlations in correlated materials simulations. An analysis of the frequency scale shows that only low-frequency propagators contribute substantially to the diagrams, putting the inclusion of higher order vertices within reach.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Efficient temperature-dependent Green's function methods for realistic systems: using cubic spline interpolation to approximate Matsubara Green's functions

The popular, stable, robust and computationally inexpensive cubic spline interpolation algorithm is adopted and used for finite temperature Green's function calculations of realistic systems. We demonstrate that with appropriate modifications the temperature dependence can be preserved while the Green's function grid size can be reduced by about two orders of magnitude by replacing the standard Matsubara frequency grid with a sparser grid and a set of interpolation coefficients. We benchmarked the accuracy of our algorithm as a function of a single parameter sensitive to the shape of the Green's function. Through numerous examples, we confirmed that our algorithm can be utilized in a systematically improvable, controlled, and black-box manner and highly accurate one- and two-body energies and one-particle density matrices can be obtained using only around 5% of the original grid points. Additionally, we established that to improve accuracy by an order of magnitude, the number of grid points needs to be doubled, whereas for the Matsubara frequency grid an order of magnitude more grid points must be used. This suggests that realistic calculations with large basis sets that were previously out of reach because they required enormous grid sizes may now become feasible.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Continuous-time Quantum Monte Carlo using Worm Sampling

We present a worm sampling method for calculating one- and two-particle Green's functions using continuous-time quantum Monte Carlo simulations in the hybridization expansion (CT-HYB). Instead of measuring Green's functions by removing hybridization lines from partition function configurations, as in conventional CT-HYB, the worm algorithm directly samples the Green's function. We show that worm sampling is necessary to obtain general two-particle Green's functions which are not of density-density type and that it improves the sampling efficiency when approaching the atomic limit. Such two-particle Green's functions are needed to compute off-diagonal elements of susceptibilities and occur in diagrammatic extensions of the dynamical mean field theory and efficient estimators for the single-particle self-energy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mechanisms of finite temperature magnetism in the three-dimensional Hubbard model

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$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

opendf - an implementation of the dual fermion method for strongly correlated systems

The dual fermion method is a multiscale approach for solving lattice problems of interacting strongly correlated systems. In this paper, we present the \texttt{opendf} code, an open-source implementation of the dual fermion method applicable to fermionic single-orbital lattice models in dimensions $D=1,2,3$ and $4$. The method is built on a dynamical mean field starting point, which neglects all local correlations, and perturbatively adds spatial correlations. Our code is distributed as an open-source package under the GNU public license version 2.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Superconducting Fluctuations in the Normal State of the Two-Dimensional Hubbard Model

We compute the two-particle quantities relevant for superconducting correlations in the two-dimensional Hubbard model within the dynamical cluster approximation. In the normal state we identify the parameter regime in density, interaction, and second-nearest-neighbor hopping strength that maximizes the $d_{x^2-y^2}$ superconducting transition temperature. We find in all cases that the optimal transition temperature occurs at intermediate coupling strength, and is suppressed at strong and weak interaction strengths. Similarly, superconducting fluctuations are strongest at intermediate doping and suppressed towards large doping and half-filling. We find a change in sign of the vertex contributions to $d_{xy}$ superconductivity from repulsive near half filling to attractive at large doping. $p$-wave superconductivity is not found at the parameters we study, and $s$-wave contributions are always repulsive. For negative second-nearest-neighbor hopping the optimal transition temperature shifts towards the electron-doped side in opposition to the van Hove singularity which moves towards hole doping. We surmise that an increase of the local interaction of the electron-doped compounds would increase $T_c$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Systematically improvable multi-scale solver for correlated electron systems

The development of numerical methods capable of simulating realistic materials with strongly correlated electrons, with controllable errors, is a central challenge in quantum many-body physics. Here we describe how a hybrid between self-consistent second order perturbation theory and exact diagonalization can be used as a multi-scale solver for such systems. Using a quantum impurity model, generated from a cluster dynamical mean field approximation to the 2D Hubbard model, as a benchmark, we show that our method allows us to obtain accurate results at a fraction of the cost of typical Monte Carlo calculations. We test the behavior of our method in multiple regimes of interaction strengths and doping of the model. The algorithm avoids difficulties such as double counting corrections, frequency dependent interactions, or vertex functions. As it is solely formulated at the level of the single-particle Green's function, it provides a promising route for the simulation of realistic materials that are currently difficult to study with other methods.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Thermodynamics of the Hubbard model on stacked honeycomb and square lattices

We present a numerical study of the Hubbard model on simply stacked honeycomb and square lattices, motivated by a recent experimental realization of such models with ultracold atoms in optical lattices. We perform simulations with different interlayer coupling and interaction strengths and obtain Néel transition temperatures and entropies. We provide data for the equation of state to enable comparisons of experiments and theory. We find an enhancement of the short-range correlations in the anisotropic lattices compared to the isotropic cubic lattice, in parameter regimes suitable for the interaction driven adiabatic cooling.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Voltage quench dynamics of a Kondo system

We examine the dynamics of a correlated quantum dot in the mixed valence regime. We perform numerically exact calculations of the current after a quantum quench from equilibrium by rapidly applying a bias voltage in a wide range of initial temperatures. The current exhibits short equilibration times and saturates upon the decrease of temperature at all times, indicating Kondo behavior both in the transient regime and in steady state. The time-dependent current saturation temperature matches the Kondo temperature at small times or small voltages; a substantially increased value is observed outside of linear response. These signatures are directly observable by experiments in the time-domain.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A comparative DMFT study of the eg-orbital Hubbard model in thin films

Heterostructures of transition-metal oxides emerged as a new route to engineer electronic systems with desired functionalities. Motivated by these developments, we study a two-orbital Hubbard model in a thin-film geometry confined along the cubic [001] direction using the dynamical mean-field theory. We contrast the results of two approximate impurity solvers (exact diagonalization and one-crossing approximation) to the results of the numerically exact continuous-time quantum Monte Carlo solver. Consistent with earlier studies, we find that the one-crossing approximation performs well in the insulating regime, while the advantage of the exact-diagonalization based solver is more pronounced in the metallic regime. We then investigate various aspects of strongly correlated eg-orbital systems in thin film geometries. In particular, we show how the interfacial orbital polarization dies off quickly a few layers from the interface and how the film thickness affects the location of the interaction-driven Mott transition. In addition, we explore the changes in the electronic structure with varying carrier concentration and identify large variations of the orbital polarization in the strongly correlated regime.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Critical Exponents of Strongly Correlated Fermion Systems from Diagrammatic Multi-Scale Methods

Self-consistent dynamical approximations for strongly correlated fermion systems are particularly successful in capturing the dynamical competition of local correlations. In these, the effect of spatially extended degrees of freedom is usually only taken into account in a mean field fashion or as a secondary effect. As a result, critical exponents associated with phase transitions have mean field character. Here, we demonstrate that diagrammatic multi-scale methods anchored around local approximations are indeed capable of capturing the non mean-field nature of the critical point of the lattice model encoded in a non-vanishing anomalous dimension, and to correctly describe the transition to mean field like behavior as the number of spatial dimensions increases.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Pairing Glue in the Two Dimensional Hubbard Model

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Efficient implementation of the continuous-time hybridization expansion quantum impurity solver

Strongly correlated quantum impurity problems appear in a wide variety of contexts ranging from nanoscience and surface physics to material science and the theory of strongly correlated lattice models, where they appear as auxiliary systems within dynamical mean-field theory. Accurate and unbiased solutions must usually be obtained numerically, and continuous-time quantum Monte Carlo algorithms, a family of algorithms based on the stochastic sampling of partition function expansions, perform well for such systems. With the present paper we provide an efficient and generic implementation of the hybridization expansion quantum impurity solver, based on the segment representation. We provide a complete implementation featuring most of the recently developed extensions and optimizations. Our implementation allows one to treat retarded interactions and provides generalized measurement routines based on improved estimators for the self-energy and for vertex functions. The solver is embedded in the ALPS-DMFT application package.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Equation of State of the Fermionic 2D Hubbard Model

We present results for the equation of state of the two-dimensional Hubbard model on an isotropic square lattice as obtained from a controlled and numerically exact large-cluster dynamical mean field simulation. Our results are obtained for large but finite systems and are extrapolated to infinite system size using a known finite size scaling relation. We present the energy, entropy, double occupancy and nearest-neighbour spin correlations extrapolated to the thermodynamic limit and discuss the implications of these calculations on pseudogap physics of the 2D-Hubbard model away from half filling. We find a strong behavioural shift in energy below a temperature $T^*$ which becomes more pronounced for larger clusters. Finally, we provide reference calculations and tables for the equation of state for values of doping away from half filling which are of interest to cold atom experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Sum rule violation in self-consistent hybridization expansions

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Superconductivity and the Pseudogap in the two-dimensional Hubbard model

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Thermodynamics and magnetic properties of the anisotropic 3D Hubbard model

We study the 3D Hubbard model with anisotropic nearest neighbor tunneling amplitudes using the dynamical cluster approximation and compare the results with a quantum simulation experiment using ultracold fermions in an optical lattice, focussing on magnetic correlations. We find that the short-range spin correlations are significantly enhanced in the direction with stronger tunneling amplitudes. Our results agree with the experimental observations and show that the experimental temperature is lower than the strong tunneling amplitude. We characterize the system by examining the spin correlations beyond neighboring sites and determine the distribution of density, entropy and spin correlation in the trapped system. We furthermore investigate the dependence of the critical entropy at the Néel transition on anisotropy.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Energetics of superconductivity in the two dimensional Hubbard model

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Truncated Configuration Interaction expansions as solvers for correlated quantum impurity models and dynamical mean field theory

The development of polynomial cost solvers for correlated quantum impurity models, with controllable errors, is a central challenge in quantum many-body physics, where these models find applications ranging from nano-science to the dynamical mean-field theory (DMFT). Here we describe how configuration interaction (CI) approximations to exact diagonalization (ED) may be used as solvers in DMFT. CI approximations retain the main advantages of ED, such as the ability to treat general interactions and off-diagonal hybridizations and to obtain real spectral information, but are of polynomial cost. Furthermore, their errors can be controlled by monitoring the convergence of physical quantities as a function of the CI hierarchy. Using benchmark DMFT applications, such as single-site DMFT of the 1D Hubbard model and $2\times 2$ cluster DMFT of the 2D Hubbard model, we show that CI approximations allow us to obtain near-exact ED results for a tiny fraction of the cost. This is true over the entire range of interaction strengths including "difficult" regimes, such as in the pseudogap phase of the 2D Hubbard model. We use the ability of CI approximations to treat large numbers of orbitals to demonstrate convergence of the bath representation in the $2\times 2$ cluster DMFT using a 24 bath orbital representation. CI approximations thus form a promising route to extend ED to problems that are currently difficult to study using other solvers such as continuous-time quantum Monte Carlo, including impurity models with large numbers of orbitals and general interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Charge and spin criticality for the continuous Mott transition in a two-dimensional organic conductor

We study the continuous bandwidth-controlled Mott transition in the two-dimensional single-band Hubbard model with a focus on the critical scaling behavior of charge and spin degrees of freedom. Using plaquette cluster dynamical mean-field theory, we find charge and spin criticality consistent with experimental results for organic conductors. In particular, the charge degree of freedom measured via the local density of states at the Fermi level shows a smoother transition than expected for the Ising universality class and in single-site dynamical mean-field theory, revealing the importance of short-ranged nonlocal correlations in two spatial dimensions. The spin criticality measured via the local spin susceptibility agrees quantitatively with nuclear magnetic resonance measurements of the spin-lattice relaxation rate.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Continuous-time Monte Carlo methods for quantum impurity models

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamical mean-field theory for bosons

We discuss the recently developed bosonic dynamical mean-field (B-DMFT) framework, which maps a bosonic lattice model onto the selfconsistent solution of a bosonic impurity model with coupling to a reservoir of normal and condensed bosons. The effective impurity action is derived in several ways: (i) as an approximation to the kinetic energy functional of the lattice problem, (ii) using a cavity approach, and (iii) by using an effective medium approach based on adding a one-loop correction to the selfconsistently defined condensate. To solve the impurity problem, we use a continuous-time Monte Carlo algorithm based on a sampling of a perturbation expansion in the hybridization functions and the condensate wave function. As applications of the formalism we present finite temperature B-DMFT phase diagrams for the bosonic Hubbard model on a 3d cubic and 2d square lattice, the condensate order parameter as a function of chemical potential, critical exponents for the condensate, the approach to the weakly interacting Bose gas regime for weak repulsions, and the kinetic energy as a function of temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Impurity Effects in Highly Frustrated Diamond Lattice Antiferromagnets

We consider the effects of local impurities in highly frustrated diamond lattice antiferromagnets, which exhibit large but non-extensive ground state degeneracies. Such models are appropriate to many A-site magnetic spinels. We argue very generally that sufficiently dilute impurities induce an ordered magnetic ground state, and provide a mechanism of degeneracy breaking. The states which are selected can be determined by a "swiss cheese model" analysis, which we demonstrate numerically for a particular impurity model in this case. Moreover, we present criteria for estimating the stability of the resulting ordered phase to a competing frozen (spin glass) one. The results may explain the contrasting finding of frozen and ordered ground states in CoAl2O4 and MnSc2S4, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerically Exact Long Time Behavior of Nonequilibrium Quantum Impurity Models

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Spectral properties of the three-dimensional Hubbard model

We present momentum resolved single-particle spectra for the three-dimensional Hubbard model for the paramagnetic and antiferromagnetically ordered phase obtained within the dynamical cluster approximation. The effective cluster problem is solved by continuous-time Quantum Monte Carlo simulations. The absence of a time discretization error and the ability to perform Monte Carlo measurements directly in Matsubara frequencies enable us to analytically continue the self-energies by maximum entropy, which is essential to obtain momentum resolved spectral functions for the N'eel state. We investigate the dependence on temperature and interaction strength and the effect of magnetic frustration introduced by a next-nearest neighbor hopping. One particular question we address here is the influence of the frustrating interaction on the metal insulator transition of the three-dimensional Hubbard model.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Sub-matrix updates for the Continuous-Time Auxiliary Field algorithm

We present a sub-matrix update algorithm for the continuous-time auxiliary field method that allows the simulation of large lattice and impurity problems. The algorithm takes optimal advantage of modern CPU architectures by consistently using matrix instead of vector operations, resulting in a speedup of a factor of $\approx 8$ and thereby allowing access to larger systems and lower temperature. We illustrate the power of our algorithm at the example of a cluster dynamical mean field simulation of the Néel transition in the three-dimensional Hubbard model, where we show momentum dependent self-energies for clusters with up to 100 sites.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Superconductivity and Pairing Fluctuations in the Half-Filled Two-Dimensional Hubbard Model

The two-dimensional Hubbard model exhibits superconductivity with d-wave symmetry even at half-filling in the presence of next-nearest neighbor hopping. Using plaquette cluster dynamical mean-field theory with a continuous-time quantum Monte Carlo impurity solver, we reveal the non-Fermi liquid character of the metallic phase in proximity to the superconducting state. Specifically, the low-frequency scattering rate for momenta near (π, 0) varies non-monotonously at low temperatures, and the dc conductivity is T-linear at elevated temperatures with an upturn upon cooling. Evidence is provided that pairing fluctuations dominate the normal-conducting state even considerably above the superconducting transition temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dynamical mean field solution of the Bose-Hubbard model

We present the effective action and self-consistency equations for the bosonic dynamical mean field (B-DMFT) approximation to the bosonic Hubbard model and show that it provides remarkably accurate phase diagrams and correlation functions. To solve the bosonic dynamical mean field equations we use a continuous-time Monte Carlo method for bosonic impurity models based on a diagrammatic expansion in the hybridization and condensate coupling. This method is readily generalized to bosonic mixtures, spinful bosons, and Bose-Fermi mixtures.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Optimized broad-histogram simulations for strong first-order phase transitions: Droplet transitions in the large-Q Potts model

The numerical simulation of strongly first-order phase transitions has remained a notoriously difficult problem even for classical systems due to the exponentially suppressed (thermal) equilibration in the vicinity of such a transition. In the absence of efficient update techniques, a common approach to improve equilibration in Monte Carlo simulations is to broaden the sampled statistical ensemble beyond the bimodal distribution of the canonical ensemble. Here we show how a recently developed feedback algorithm can systematically optimize such broad-histogram ensembles and significantly speed up equilibration in comparison with other extended ensemble techniques such as flat-histogram, multicanonical or Wang-Landau sampling. As a prototypical example of a strong first-order transition we simulate the two-dimensional Potts model with up to Q=250 different states on large systems. The optimized histogram develops a distinct multipeak structure, thereby resolving entropic barriers and their associated phase transitions in the phase coexistence region such as droplet nucleation and annihilation or droplet-strip transitions for systems with periodic boundary conditions. We characterize the efficiency of the optimized histogram sampling by measuring round-trip times tau(N,Q) across the phase transition for samples of size N spins. While we find power-law scaling of tau vs. N for small Q \lesssim 50 and N \lesssim 40^2, we observe a crossover to exponential scaling for larger Q. These results demonstrate that despite the ensemble optimization broad-histogram simulations cannot fully eliminate the supercritical slowing down at strongly first-order transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quenched dynamics in interacting one-dimensional systems: Appearance of current carrying steady states from initial domain wall density profiles

We investigate dynamics arising after an interaction quench in the quantum sine-Gordon model for a one-dimensional system initially prepared in a spatially inhomogeneous domain wall state. We study the time-evolution of the density, current and equal time correlation functions using the truncated Wigner approximation (TWA) to which quantum corrections are added in order to set the limits on its validity. For weak to moderate strengths of the back-scattering interaction, the domain wall spreads out ballistically with the system within the light cone reaching a nonequilibrium steady-state characterized by a net current flow. A steady state current exists for a quench at the exactly solvable Luther-Emery point. The magnitude of the current decreases with increasing strength of the back-scattering interaction. The two-point correlation function of the variable canonically conjugate to the density reaches a spatially oscillating steady state at a wavelength inversely related to the current.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Performance analysis of continuous-time solvers for quantum impurity models

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

0 citations0 reviewsNo code linkedOpen access