REVIEW 3 major objections 4 minor 52 references
Strange diffusivity of incoherent metal in half-filled two-dimensional Hubbard model
T0 review · 3 major / 4 minor · reviewed 2026-08-05 · deepseek-v4-flash
Pith's one-line read The half-filled 2D Hubbard model's incoherent metal has a charge diffusion constant D ~ 1/√T across all anomalous-resistivity regimes down to the metal-insulator crossover, plus a Pseudogap Metal state where insulating compressibility coexi
desk verdict D~T^-1/2 diffusion in the half-filled Hubbard model is the new result; the analytic continuation is the main caveat, but the paper deserves serious refereeing. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central computational object is the imaginary-time current-current correlator Λ(iωn), continued via its spectral representation to the optical conductivity σ(ω) using the SOCC/MCC method with a multiscale high-frequency tail from Bold4. The central physical identity is the Nernst-Einstein relation D = σDC/κ, which particle-hole symmetry at half-filling makes exact because thermoelectric response vanishes. Its diagnostic power lies in decomposing anomalous resistivity into a universal diffusion law (1/√T) and a non-universal compressibility; the diagrammatic spin-resolved decomposition separates bubble, same-spin, and opposite-spin vertex contributions that explain the Drude-to-continuum
What would settle it
A measurement of the optical conductivity in the half-filled 2D Hubbard model that does not rest on maximum-smoothness analytic continuation—for example a real-time dynamical quantum Monte Carlo calculation or a cold-atom density-response measurement at half-filling—should find σDC at U=4,T=0.4-4 within the stretch-test window; if σDC deviates significantly, then D=σ/κ will not follow the reported 1/√T law.
Extended reading notes
Core claim
For the half-filled 2D Hubbard model, the authors compute the current-current correlation function directly in the thermodynamic limit with diagrammatic Monte Carlo, splice the high-frequency tail with the self-consistent Bold4 diagrammatic theory, and continue to real frequencies with a stretch-test-controlled stochastic method. They find that the DC resistivity's anomalous scaling between high and low temperatures is the product of a near-universal diffusivity D∼1/√T multiplied by a compressibility that is non-universal; at lower temperatures the compressibility turns over and the system enters a Pseudogap Metal state with ∂κ/∂T>0 but metallic ρDC(T), until the metal-insulator crossover. D
Load-bearing premise
The numerical analytic continuation must recover the true optical conductivity near zero frequency; the smoothness-constrained spectrum is validated by a stretch test but is not uniquely determined by the Matsubara data.
Editorial extensions
If this is right
- The 1/√T diffusivity extends through Diffusive Metal I, Diffusive Metal II, and the Pseudogap Metal, so in that window the anomalous resistivity exponent α is controlled mostly by the compressibility, not by the diffusion constant.
- Because particle-hole symmetry makes the Nernst-Einstein extraction exact at half-filling, the reported D is a genuine charge diffusion constant rather than a thermoelectric mixture.
- The Pseudogap Metal state is defined by dκ/dT>0 together with metallic transport, showing that insulating charge response and metallic conduction can coexist without an actual gap in the conductivity.
- The opposite-spin vertex corrections deplete the low-frequency Drude peak and feed a high-frequency continuum, explaining why the central optical peak is suppressed before the system becomes insulating.
Reading between the lines
- One implication beyond the paper's claims is that D∼1/√T may be universal across doping, with all material-specific and model-specific variation residing in the compressibility; the cold-atom data at substantial doping already point in this direction.
- A testable extension would be a real-time measurement of charge diffusion in an ultracold-atom realization at half-filling: if the Pseudogap Metal picture is right, density-response (compressibility) and cloud-expansion (transport) probes should show opposite temperature trends.
- The spin-resolved vertex mechanism suggests an optical sum-rule-style diagnostic: the same-spin vertex adds spectral weight in a window set by temperature, while the opposite-spin vertex removes it from the Drude peak; this separation could be probed by spin-resolved or polarized light experiments in analog systems.
- If the 1/√T law is the universal incoherent-metal signature, then future analytic-continuation studies should report D rather than only ρDC, since D is the quantity that exposes the regularity behind the non-universal resistivity exponents.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies charge transport across the metal-insulator crossover in the half-filled two-dimensional Hubbard model. The authors compute the current-current correlation function directly in the thermodynamic limit with connected determinant diagrammatic Monte Carlo (CDet) up to order 8-10, combine it with a self-consistent Bold4 high-frequency tail, and use numerical analytic continuation (SOCC/MCC with a 'stretch test') to obtain the optical conductivity. Combining the DC conductivity with the compressibility computed earlier by the same group, they extract the charge diffusion constant D = sigma_DC / kappa. The central claim is that, across a broad temperature range where the DC resistivity has anomalous scaling rho_DC ~ T^alpha with 0<alpha<1, the diffusion constant displays a robust ~1/sqrt(T) 'strange metal' behavior. They also identify a 'Pseudogap Metal' regime characterized by insulating charge compressibility coexisting with metallic transport, and analyze the diagrammatic origin of this behavior, finding that opposite-spin vertex corrections transfer Drude weight to a high-frequency continuum.
Significance. If the reported 1/sqrt(T) law holds, it would provide a striking universal diffusive signature of incoherent transport in the half-filled 2D Hubbard model, with a direct connection to cold-atom measurements in the doped system. The numerical work is careful: the CDet series are pushed to high order and extrapolated with Pade approximants, error bars are propagated into observables, and the analytic continuation is guarded by a stretch test and cross-checked against SOM at one representative point. The diagrammatic decomposition also benefits from an exact order-by-order selection rule for spin-resolved vertex contributions. These strengths make the paper a serious candidate for publication, provided the principal scaling claim is robust against the known ill-posedness of the analytic continuation.
major comments (3)
- [SM Sec. III, Eq. (4)] The stretch test defines the sigma(0) error bar as the range over which the MCC spectrum does not develop visible wiggles; this is a smoothness-prior sensitivity test, not a uniqueness bound. Because the kernel K(i omega_n, omega) has small eigenvalues, non-smooth spectra with different sigma(0) can fit Lambda(i omega_n) within the DiagMC error bars. The main claim D(T) ~ T^{-1/2} is obtained from sigma_DC(T)/kappa(T) across a range of T, so a T-dependent continuation bias could fake or distort the power law. Please cross-validate with an independent continuation method (e.g., MaxEnt, SOM, sparse modeling, or Nevanlinna continuation) at enough (U,T) points to map D(T), and/or demonstrate with synthetic spectra that non-smooth solutions consistent with Lambda(i omega_n) change D(T) by less than the quoted error.
- [SM Sec. III, Fig. S2(a)] The high-frequency tail for omega_n > 15 is replaced by Bold4 values with an ad hoc O(1/omega_n^3) uncertainty. This tail enters as data in Eq. (4) and can influence the low-frequency inversion through overall normalization and constraints. The systematic difference between Bold4 and the exact CCF is acknowledged but is not propagated into the final D(T) error bars in a well-controlled way. Please test the sensitivity of sigma(0) to the substitution threshold and to the assumed tail-error model, and include the resulting uncertainty in the reported diffusion constant.
- [Fig. 2(c), Table I] The central claim of a robust ~1/sqrt(T) scaling is supported visually by reference lines, but no fitting procedure, exponent uncertainty, or temperature-window definition is reported. Provide a power-law fit D(T) = A T^{-gamma} for each U with statistical errors, state the criterion for the fitting window (e.g., DM1 through PGM), and report gamma and its uncertainty. This is needed to substantiate the universality statement and to make the claimed T^{-1/2} exponent falsifiable.
minor comments (4)
- [Fig. 3] The panel labels and caption are confusing: both panels appear to be labeled 'a' and 'b', and the inset showing sigma_DC * Gamma_tr and D * Gamma_tr is not fully described. Please clarify the figure layout and the definition of Gamma_tr in the caption.
- [SM Sec. IV] The temperature T_II is defined through the inflection point of a 5th-order polynomial fit to d kappa / d log T. Please comment on the stability of T_II with respect to the fit order and the fitting range.
- [Table I] The quoted ranges of the resistivity exponent alpha (0 <~ alpha <~ 0.5 and 0.5 <~ alpha <~ 1) are not tied to explicit fits. These ranges should be obtained from and reported with the same fitting procedure used for the diffusion constant.
- [References] Reference [48] is an arXiv preprint; if a published version exists, please cite it.
Circularity Check
No significant circularity: the D(T) scaling is read off independently computed σDC and κ; analytic-continuation systematics are a robustness concern, not a circular reduction.
full rationale
The derivation chain is: Λ(iωn) is computed by DiagMC directly in the thermodynamic limit; the high-frequency tail is taken from Bold4 with an explicitly assigned O(1/ω_n^3) uncertainty; σ(ω) is obtained by numerical analytic continuation; σDC=σ(0); and D=σDC/κ with κ computed by DiagMC (with first-order analytic complements). Each step is an independent numerical evaluation. No quantity along this chain is defined in terms of the claimed D(T)∼1/√T scaling, and no parameter is fitted to the ρDC(T) or D(T) behavior that is then presented as a prediction. The paper explicitly frames the NAC error using a stretch test around the "as-smooth-as-possible" solution, which is a sensitivity statement about smoothness-constrained spectra, not a uniqueness certificate; the possibility that a less smooth true spectrum could bias σ(0) is a correctness/robustness risk, not a circularity of the paper's own equations. The uses of the authors' prior results (Ref. [13] for compressibility and Ref. [15] for self-energy anisotropy) are citations to separate computations of different observables, not inputs that already contain the transport conclusion; they provide independent, externally falsifiable information. Thus no load-bearing step reduces to its own input by construction.
Assumptions & free parameters
assumptions (5)
- domain assumption Particle-hole symmetry at half-filling gives zero thermoelectric response, so the Nernst-Einstein relation D=σ/κ holds in its simple form.
- domain assumption The diagrammatic expansion of the current-current correlator converges and the Padé extrapolation to infinite order is reliable within the quoted error bars.
- ad hoc to paper The numerical analytic continuation (SOCC/MCC) with the stretch test yields the true optical conductivity σ(ω), particularly σ(0), within the stated error bars.
- ad hoc to paper The Bold4 self-consistent diagrammatic results accurately describe the high-frequency tail of the current correlator, with the introduced O(1/ω_n^3) error bounding the uncertainty.
- domain assumption Perfect nesting and Mermin-Wagner theorem ensure the finite-temperature transport is incoherent and the phase diagram is as described.
Cite this review
Pith. "Pith review of Strange diffusivity of incoherent metal in half-filled two-dimensional Hubbard model." pith.science (2026). https://pith.science/paper/NQAANOZB
@misc{pith2026250900281,
author = {Pith},
title = {Pith review of: Strange diffusivity of incoherent metal in half-filled two-dimensional Hubbard model},
year = {2026},
howpublished = {\url{https://pith.science/paper/NQAANOZB}},
note = {Machine review of arXiv:2509.00281}
}
abstract
We study charge transport across the metal-insulator crossover in the half-filled two-dimensional Hubbard model, with particular emphasis on precision control. The dynamic current-current correlation function is obtained directly in the thermodynamic limit, and the optical conductivity is extracted using numerical analytic continuation. To achieve this, we develop a multiscale approach: the non-perturbative low-frequency behavior is computed using the unbiased diagrammatic Monte Carlo technique, while the high-frequency physics is captured via a self-consistent (semi-)analytic diagrammatic theory. We found that across a broad temperature range where the DC resistivity displays anomalous scaling, $\sim T^\alpha$ with $0<\alpha\lesssim 1$, the Nernst-Einstein relation implies the diffusion constant with the characteristic $\sim 1/\sqrt{T}$ "strange metal" behavior. It was also revealed that the insulating regime is entered through a peculiar non-Fermi liquid state-which we call a Pseudogap Metal-characterized by insulating charge compressibility coexisting with metallic transport. Diagrammatically, the high-temperature incoherent transport is captured by the dressed polarization bubble, whereas near the metal-insulator crossover, the effective interaction vertex between opposite-spin particles is responsible for transferring the Drude weight to a high-frequency continuum.
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[1]
J. A. N. Bruin, H. Sakai, R. S. Perry, and A. P. Macken- zie, Similarity of scattering rates in metals showing t- linear resistivity, Science 339, 804 (2013)
work page 2013
-
[2]
A. F. Ioffe and A. R. Regel, Non-crystalline semiconduc- tors, Prog. Semicond. 4, 237 (1960)
work page 1960
-
[3]
N. F. Mott, Conduction in non-crystalline materials, Phi- los. Mag. 26, 1015 (1972)
work page 1972
-
[4]
O. Gunnarsson, M. Calandra, and J. E. Han, Colloquium: Saturation of electrical resistivity, Rev. Mod. Phys. 75, 1085 (2003). 6
work page 2003
-
[5]
N. E. Hussey, K. Takenaka, and H. Takagi, Universal- ity of the mott–ioffe–regel limit in metals, Philosophical Magazine 84, 2847 (2004)
work page 2004
-
[6]
S. A. Hartnoll, Theory of universal incoherent metallic transport, Nat. Phys. 11, 54 (2015)
2015
-
[7]
P. T. Brown, D. Mitra, E. Guardado-Sanchez, R. Nourafkan, A. Reymbaut, C.-D. H´ ebert, S. Bergeron, A. M. S. Tremblay, J. Kokalj, D. A. Huse, P. Schauß, and W. S. Bakr, Bad metallic transport in a cold atom fermi-hubbard system, Science 363, 379–382 (2019)
work page 2019
-
[8]
P. A. Lee and N. Read, Why is t c of the oxide supercon- ductors so low?, Phys. Rev. Lett. 58, 2691 (1987)
work page 1987
Show all 52 references
-
[9]
Virosztek and J
A. Virosztek and J. Ruvalds, Nested-fermi-liquid theory, Phys. Rev. B 42, 4064 (1990)
1990
-
[10]
T. G. Kiely and E. J. Mueller, Transport in the two- dimensional fermi-hubbard model: Lessons from weak coupling, Phys. Rev. B 104, 165143 (2021)
2021
-
[11]
Vuˇ ciˇ cevi´ c, S
J. Vuˇ ciˇ cevi´ c, S. Predin, and M. Ferrero, Charge fluctua- tions, hydrodynamics, and transport in the square-lattice hubbard model, Phys. Rev. B 107, 155140 (2023)
2023
-
[12]
Beni and C
G. Beni and C. F. Coll, Thermoelectric power in half- filled bands, Phys. Rev. B 11, 573 (1975)
1975
-
[13]
A. J. Kim, F. Simkovic, and E. Kozik, Spin and charge correlations across the metal-to-insulator crossover in the half-filled 2d hubbard model, Phys. Rev. Lett. 124, 117602 (2020)
2020
-
[14]
Kokalj, Bad-metallic behavior of doped mott insula- tors, Phys
J. Kokalj, Bad-metallic behavior of doped mott insula- tors, Phys. Rev. B 95, 041110 (2017)
2017
-
[15]
ˇSimkovic, J
F. ˇSimkovic, J. P. F. LeBlanc, A. J. Kim, Y. Deng, N. V. Prokof’ev, B. V. Svistunov, and E. Kozik, Ex- tended crossover from a fermi liquid to a quasiantiferro- magnet in the half-filled 2d hubbard model, Phys. Rev. Lett. 124, 017003 (2020)
2020
-
[16]
Lenihan, A
C. Lenihan, A. J. Kim, F. ˇSimkovic Iv, and E. Kozik, Entropy in the non-fermi-liquid regime of the doped 2D hubbard model, Phys. Rev. Lett. 126, 105701 (2021)
2021
-
[17]
Schaefer, N
T. Schaefer, N. Wentzell, F. Simkovic, Y. Y. He, C. Hille, M. Klett, C. J. Eckhardt, B. Arzhang, V. Harkov, F. M. L. Regent, A. Kirsch, Y. Wang, A. J. Kim, E. Kozik, E. A. Stepanov, A. Kauch, S. Andergassen, P. Hansmann, D. Rohe, Y. M. Vilk, J. P. F. LeBlanc, S. Zhang, A.-M. S...
2021
-
[18]
N. V. Prokof’ev and B. V. Svistunov, Polaron problem by diagrammatic quantum monte carlo, Phys. Rev. Lett. 81, 2514 (1998)
1998
-
[19]
Kozik, K
E. Kozik, K. V. Houcke, E. Gull, L. Pollet, N. Prokof’ev, B. Svistunov, and M. Troyer, Diagrammatic monte carlo for correlated fermions, EPL 90, 10004 (2010)
2010
-
[20]
K. V. Houcke, E. Kozik, N. Prokof’ev, and B. Svistunov, Diagrammatic monte carlo, Phys. Procedia 6, 95 (2010)
2010
-
[21]
Rossi, Determinant diagrammatic monte carlo algo- rithm in the thermodynamic limit, Phys
R. Rossi, Determinant diagrammatic monte carlo algo- rithm in the thermodynamic limit, Phys. Rev. Lett. 119, 045701 (2017)
2017
-
[22]
See Supplemental Material at [URL will be inserted by publisher] for more details on the current–current cor- relation function, diagrammatic series convergence, the stretch test in numerical analytic continuation, the ef- fects of the van Hove singularity, as well as vertex c...
-
[23]
A. S. Mishchenko, N. V. Prokof’ev, A. Sakamoto, and B. V. Svistunov, Diagrammatic quantum monte carlo study of the fr¨ ohlich polaron, Phys. Rev. B 62, 6317 (2000)
2000
-
[24]
A. S. Mishchenko, in Correlated Electrons: From Mod- els to Materials Modeling and Simulation, edited by E. Pavarini, E. Koch, F. Anders, and M. Jarrell (Verlag des Forschungszentrum, Julich, 2012), Vol. 2
2012
-
[25]
Goulko, A
O. Goulko, A. S. Mishchenko, L. Pollet, N. Prokof’ev, and B. Svistunov, Numerical analytic continuation: An- swers to well-posed questions, Phys. Rev. B 95, 014102 (2017)
2017
-
[26]
ˇSimkovic, Y
F. ˇSimkovic, Y. Deng, N. V. Prokof’ev, B. V. Svistunov, I. S. Tupitsyn, and E. Kozik, Magnetic correlations in the two-dimensional repulsive fermi-hubbard model, Phys. Rev. B 96, 081117 (2017)
2017
-
[27]
K. S. D. Beach, Identifying the maximum entropy method as a special limit of stochastic analytic continu- ation, arXiv:cond-mat/0403055
-
[28]
A. W. Sandvik, Stochastic method for analytic contin- uation of quantum monte carlo data, Phys. Rev. B 57, 10287 (1998)
1998
-
[29]
Fuchs, M
S. Fuchs, M. Jarrell, and T. Pruschke, Application of bayesian inference to stochastic analytic continua- tion, Journal of Physics: Conference Series 200, 012041 (2010)
2010
-
[30]
Yoshimi, J
K. Yoshimi, J. Otsuki, Y. Motoyama, M. Ohzeki, and H. Shinaoka, Spm: Sparse modeling tool for analytic con- tinuation of imaginary-time green’s function, Computer Physics Communications 244, 319 (2019)
2019
-
[31]
Motoyama, K
Y. Motoyama, K. Yoshimi, and J. Otsuki, Robust ana- lytic continuation combining the advantages of the sparse modeling approach and the pad´ e approximation, Phys. Rev. B 105, 035139 (2022)
2022
-
[32]
Yoon, J.-H
H. Yoon, J.-H. Sim, and M. J. Han, Analytic continuation via domain knowledge free machine learning, Phys. Rev. B 98, 245101 (2018)
2018
-
[33]
Fournier, L
R. Fournier, L. Wang, O. V. Yazyev, and Q. Wu, Artifi- cial neural network approach to the analytic continuation problem, Phys. Rev. Lett. 124, 056401 (2020)
2020
-
[34]
Fei, C.-N
J. Fei, C.-N. Yeh, and E. Gull, Nevanlinna analytical con- tinuation, Phys. Rev. Lett. 126, 056402 (2021)
2021
-
[35]
Fei, C.-N
J. Fei, C.-N. Yeh, D. Zgid, and E. Gull, Analytical con- tinuation of matrix-valued functions: Carath´ eodory for- malism, Phys. Rev. B 104, 165111 (2021)
2021
-
[36]
Huang, E
Z. Huang, E. Gull, and L. Lin, Robust analytic continua- tion of green’s functions via projection, pole estimation, and semidefinite relaxation, Phys. Rev. B 107, 075151 (2023)
2023
-
[37]
Zhang and E
L. Zhang and E. Gull, Minimal pole representation and controlled analytic continuation of matsubara response functions, Phys. Rev. B 110, 035154 (2024)
2024
-
[38]
Zhang, Y
L. Zhang, Y. Yu, and E. Gull, Minimal pole representa- tion and analytic continuation of matrix-valued correla- tion functions, Phys. Rev. B 110, 235131 (2024)
2024
-
[39]
G. A. Baker and P. Graves-Morris, Pade Approxi- mants: Encyclopedia of Mathematics and It’s Applica- tions (Cambridge University Press, Cambridge, 1996)
1996
-
[40]
Sch¨ ott, I
J. Sch¨ ott, I. L. M. Locht, E. Lundin, O. Gr ˚ an¨ as, O. Eriks- son, and I. Di Marco, Analytic continuation by averaging pad´ e approximants, Phys. Rev. B93, 075104 (2016)
2016
-
[41]
R. K. Bryan, Maximum entropy analysis of oversampled data problems, Eur. Biophys. J. 18, 165 (1990)
1990
-
[42]
Jarrell and J
M. Jarrell and J. E. Gubernatis, Bayesian inference and the analytic continuation of imaginary-time quantum 1 monte carlo data, Phys. Rep. 269, 133 (1996)
1996
-
[43]
Perepelitsky, A
E. Perepelitsky, A. Galatas, J. Mravlje, R. ˇZitko, E. Khatami, B. S. Shastry, and A. Georges, Transport and optical conductivity in the hubbard model: A high- temperature expansion perspective, Phys. Rev. B 94, 235115 (2016)
2016
-
[44]
Technically, TI can also be defined by the condition ∂2κ/∂ log T 2 = 0, in agreement with the definition in the main text within the error bars [22]
-
[45]
E. W. Huang, R. Sheppard, B. Moritz, and T. P. Dev- ereaux, Strange metallicity in the doped hubbard model, Science 366, 987 (2019)
2019
-
[46]
Since no time scale is supposed to be larger than the temperature, in this case, the HWHM should be inter- preted as a signal of the strong renormalization of other quantities
-
[47]
Kauch, P
A. Kauch, P. Pudleiner, K. Astleithner, P. Thunstr¨ om, T. Ribic, and K. Held, Generic optical excitations of cor- related systems: π-tons, Phys. Rev. Lett. 124, 047401 (2020)
2020
-
[48]
Dasari, H
N. Dasari, H. U. R. Strand, M. Eckstein, A. I. Lichten- stein, and E. A. Stepanov, Non-Local Correlation Effects in DC and Optical Conductivity of the Hubbard Model, arXiv:2507.16673
-
[49]
ˇSimkovic and M
F. ˇSimkovic and M. Ferrero, Fast principal minor algo- rithms for diagrammatic monte carlo, Phys. Rev. B 105, 125104 (2022)
2022
-
[50]
N. V. Prokof’ev and B. V. Svistunov, Spectral analysis by the method of consistent constraints, JETP Letters 97, 649 (2013)
2013
-
[51]
W. H. Furry, A symmetry theorem in the positron theory, Phys. Rev. 51, 125 (1937)
1937
-
[52]
W. Wu, M. Ferrero, A. Georges, and E. Kozik, Control- ling feynman diagrammatic expansions: Physical nature of the pseudogap in the two-dimensional hubbard model, Physical Review B 96, 041105 (2017). Supplemental Materials: Strange diffusivity of incoherent metal in half-fille...
2017
Reviewed August 5, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.