2011
DOI: 10.1103/physrevb.83.075122
|View full text |Cite
|
Sign up to set email alerts
|

Submatrix updates for the continuous-time auxiliary-field algorithm

Abstract: We present a submatrix 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 ≈ 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 i… Show more

Help me understand this report
View preprint versions

Search citation statements

Order By: Relevance

Paper Sections

Select...
1
1
1

Citation Types

2
89
1

Year Published

2013
2013
2022
2022

Publication Types

Select...
6
2

Relationship

5
3

Authors

Journals

citations
Cited by 91 publications
(92 citation statements)
references
References 44 publications
2
89
1
Order By: Relevance
“…We describe only a brief outline of CT-INT in this section. For efficient computation, we also employ the submatrix update algorithm 34,35) extended to CT-INT. 33) We divide the impurity Hamiltonian into two parts, H = H 0 + H i .…”
Section: Continuous Time Quantum Monte Carlo Methodsmentioning
confidence: 99%
See 1 more Smart Citation
“…We describe only a brief outline of CT-INT in this section. For efficient computation, we also employ the submatrix update algorithm 34,35) extended to CT-INT. 33) We divide the impurity Hamiltonian into two parts, H = H 0 + H i .…”
Section: Continuous Time Quantum Monte Carlo Methodsmentioning
confidence: 99%
“…Given this background, we analyze the one-particle spectrum of a single band model of a cuprate superconductor near the Fermi level using dynamical mean field theory (DMFT) 21) with two kinds of impurity solvers : iterated perturbation theory (IPT) [22][23][24][25][26][27][28][29] and continuous time quantum Monte Carlo (CT-QMC) [30][31][32][33][34][35] We find that the electron-hole asymmetry can exist even under common interaction strengths between the hole-and the electron-doped systems. The asymmetry of the spectrum already exists in the non-interacting case, but it is drastically enhanced, especially when the interaction is strong enough to make the non-doped case a Mott insulator.…”
Section: Introductionmentioning
confidence: 99%
“…More recent developments [18] have enabled researchers to access clusters large enough to obtain a reasonable picture of the N → ∞ limit [14,16,[47][48][49][50][51]. It has been found [16] that in DCA clusters of size N > 4 the Mott transition is multistaged, with the fully gapped Mott insulating state being separated from the Fermi liquid state by an intermediate phase, in which regions of momentum space near the (0, π)/(π, 0) point are gapped and regions of momentum space near (±π/2, ±π/2) are not.…”
mentioning
confidence: 99%
“…Specifics of our methods are given in the Supplemental Material; here we briefly note that a key aspect of our study is the use of recently developed "submatrix update" numerical techniques [17][18][19] which enable access to couplings strong enough to produce a pseudogap at temperatures low enough to construct the superconducting state for cluster size N large enough to reasonably represent the N → ∞ limit. Our key results are that the pseudogap and superconductivity are competing phases and that, remarkably, the onset of superconductivity within the pseudogap phase leads to a decrease in the excitation gap, in sharp contrast to conventional situations where the onset of superconductivity increases the gap.…”
mentioning
confidence: 99%
“…With the introduction of continuous-time Monte Carlo solvers [10][11][12][13][14], it has become possible to directly measure the self-energy on the Matsubara axis with unprecedented speed [15,16] and accuracy [17]. This motivated us to investigate the possibility of an analytic continuation of the self-energy directly in frequency space with the relationship (z) = 0 + 1 2π…”
Section: S = − Dωa(ω) − M(ω) − A(ω) Log[a(ω)/m(ω)] (3)mentioning
confidence: 99%