Pith. Sign in

REVIEW

The one-dimensional Hubbard model with open ends: Universal divergent contributions to the magnetic susceptibility

Strongly Correlated Electrons · Submitted 22 Feb 2006 · 41 references

Pith's verdict: The T=0 boundary divergence in h is solid; the T>0 boundary divergence is asserted from the Heisenberg analogy and never checked, so the paper's universality claim needs to be read with that split.

T0 · AI review

How this review was made: An AI referee (deepseek-v4-flash) read the paper on August 28, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 4 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDcond-mat/0602511 (version 2)
Pith IDLETUPD55, Pith's permanent record of this paper and its review
CitationJ. Sirker, M. Bortz. "The one-dimensional Hubbard model with open ends: Universal divergent contributions to the magnetic susceptibility." arXiv:cond-mat/0602511 (2006). https://pith.science/paper/LETUPD55
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

The boundary contribution to the magnetic susceptibility of the one-dimensional Hubbard model with open ends diverges as χB(h)=1/(4h ln²(h0/h)) at small field and χB(T)=−1/(12T ln T) at low temperature, with prefactors independent of…

The reading

The paper establishes that in the one-dimensional Hubbard model with open ends, the boundary contribution χB to the magnetic susceptibility diverges at small magnetic field h and low temperature T as χB(h)=1/(4h ln²(h0/h)) and χB(T)=−1/(12T ln T), respectively, up to next-to-leading logarithms. These leading prefactors are independent of the electron filling and of the interaction strength, and coincide with the spin-1/2 Heisenberg chain result. For the bulk susceptibility, the paper derives the logarithmic correction with the scale h0 determined exactly from the Bethe ansatz at arbitrary filling. It also predicts Friedel oscillations in the local magnetization and susceptibility near the boundary, confirmed numerically with the density-matrix renormalization group applied to transfer matrices. The result matters because a Curie-like boundary term can appear in clean one-dimensional conductors with non-magnetic defects and could be misread as magnetic-impurity physics.

What carries the argument

The load-bearing object is the function P2(v) (and its companion P1(v)) defined by a Wiener-Hopf integral equation with kernel κ(1)(v), the same auxiliary function that controls the boundary susceptibility of the spin-1/2 Heisenberg chain. In the Hubbard model with open boundaries, the Bethe-ansatz integral equations for the root densities reduce, in the low-field limit, to the same equation; the paper takes P2 from the Heisenberg-chain solution, so the entire boundary divergence is carried by this function. The charge sector enters only through renormalized velocities and scales, and the SU(2) symmetry pins the spin Luttinger parameter Ks=1, which is why the prefactor 1/4 survives.

Extended reading notes

Core claim

At low energies the Hubbard Hamiltonian separates into commuting charge and spin sectors, and the SU(2) symmetry fixes the Luttinger parameter in the spin sector to Ks=1. Because the boundary contribution to the susceptibility in the Heisenberg chain depends only on the solution P2(v) of a Wiener-Hopf integral equation and not on the spin velocity, the same P2(v) carries over to the Hubbard model; charge degrees of freedom only renormalize the spin velocity and the logarithmic scales h0 and T0. Consequently, χB(h)=1/(4h ln²(h0/h))(1−ln ln(h0/h)/ln(h0/h)+...) at T=0 and χB(T)=−1/(12T ln(T0/T))(1−ln ln(T0/T)/(2 ln(T0/T))+...) at h=0, with the leading small-field and small-temperature forms fully universal and identical to the Heisenberg result. For the bulk, χbulk(h)=χ0(1−1/(2 ln(h/h0))−ln|ln(h/h0)|/(4 ln²(h/h0))+5/(16 ln²(h/h0))) with χ0=1/(2πvs) and h0 determined exactly from the Bethe ansatz, Eq. (76).

Load-bearing premise

The calculation assumes that cutting the chain at its ends does not introduce any new local physics coupling spin to charge, so the same boundary response function that solves the Heisenberg chain applies to the Hubbard model; if cutting creates such a coupling, the universal prefactor 1/4 in χB(h) would change.

Editorial extensions

If this is right

  • A measurement of the low-temperature or low-field magnetic response of a clean one-dimensional conductor with open ends should show a Curie-like divergent boundary term whose logarithmic prefactor is set solely by the field or temperature.
  • Non-magnetic impurities or structural kinks that effectively cut a one-dimensional chain will produce the same universal boundary susceptibility as magnetic impurities would, so Curie tails cannot be used by themselves as evidence of magnetic scatterers.
  • The bulk susceptibility at zero temperature is now pinned down to next-to-leading logarithms with an exactly determined scale h0 for arbitrary filling, so quantitative fits of χ(h) data have no free scale parameter.
  • The conformal-field-theory formulas for Friedel oscillations in the local magnetization and susceptibility can be used to extract velocities and exponents from local probes such as NMR Knight shift measurements.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If this universality is robust, it should hold for any one-dimensional itinerant system whose low-energy theory is an SU(2)-symmetric spin-1/2 Heisenberg sector plus a gapless charge sector; a direct numerical test in a non-integrable chain would separate the role of integrability from that of spin-charge separation.
  • The logarithmic temperature dependence of the effective Curie constant implies that standard Curie-tail analysis should be modified: plotting χB(T)·T versus ln T should reveal a straight line with slope −1/12 rather than a constant.
  • The paper leaves the next-to-leading constant γ_B^(h) undetermined; until it is fixed by a higher-order Wiener-Hopf calculation, fits of χB(h) at intermediate fields contain one free scale parameter h0, which is independent of the exact bulk scale h0.
Original abstract

The magnetic susceptibility of the one-dimensional Hubbard model with open boundary conditions at arbitrary filling is obtained from field theory at low temperatures and small magnetic fields, including leading and next-leading orders. Logarithmic contributions to the bulk part are identified as well as algebraic-logarithmic divergences in the boundary contribution. As a manifestation of spin-charge separation, the result for the boundary part at low energies turns out to be independent of filling and interaction strength and identical to the result for the Heisenberg model. For the bulk part at zero temperature, the scale in the logarithms is determined exactly from the Bethe ansatz. At finite temperature, the susceptibility profile as well as the Friedel oscillations in the magnetisation are obtained numerically from the density-matrix renormalisation group applied to transfer matrices. Agreement is found with an exact asymptotic expansion of the relevant correlation function.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 4,962 input tokens, $0.0002 reported.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.16 · human participants 0.01 · controlled comparison 0.07 · randomized assignment 0.01 · states quantitative result 0.14 · compares with prior results 0.89 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v2. Results and the main claim are underlined; .

The magnetic susceptibility of the one-dimensional Hubbard model with open boundary conditions at arbitrary filling is obtained from field theory at low temperatures and small magnetic fields, including leading and next-leading orders. Logarithmic contributions to the bulk part are identified as well as algebraic-logarithmic divergences in the boundary contribution. As a manifestation of spin-charge separation, the result for the boundary part at low energies turns out to be independent of filling and interaction strength and identical to the result for the Heisenberg model. For the bulk part at zero temperature, the scale in the logarithms is determined exactly from the Bethe ansatz. At finite temperature, the susceptibility profile as well as the Friedel oscillations in the magnetisation are obtained numerically from the density-matrix renormalisation group applied to transfer matrices. Agreement is found with an exact asymptotic expansion of the relevant correlation function.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 The magnetic susceptibility of the one-dimensional Hubbard model with open boundary conditions at arbitrary filling is obtained from field theory at low temperatures and small magnetic fields, including leading and next-leading orders. contribution 0.94 0.93 0.22 asserted 1.00
s1 Logarithmic contributions to the bulk part are identified as well as algebraic-logarithmic divergences in the boundary contribution. result 0.99 0.90 0.16 asserted 1.00
s2 As a manifestation of spin-charge separation, the result for the boundary part at low energies turns out to be independent of filling and interaction strength and identical to the result for the Heisenberg model. interpretation 0.67 0.89 0.42 asserted 1.00
s3 For the bulk part at zero temperature, the scale in the logarithms is determined exactly from the Bethe ansatz. result 0.71 0.79 0.07 asserted 1.00
s4 At finite temperature, the susceptibility profile as well as the Friedel oscillations in the magnetisation are obtained numerically from the density-matrix renormalisation group applied to transfer matrices. method 0.81 0.94 0.05 asserted 1.00
s5 Agreement is found with an exact asymptotic expansion of the relevant correlation function. result 0.99 0.89 0.08 asserted 1.00

What kind of paper the classifier reports this is: theoretical result 1.00 · evidence from proof or derivation 0.90.

What was examined, and what was not

sources read
  • arXiv listing abstract v2
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. Spin- (black point-dashed lines) and charge- (black full lines) velocities at fillings n = 1/4, 1/2, 3/4 (pairs from bottom to top). The dashed lines on the left are the low-u results Eqs. (10), (13). The horizontal bars on the right indicate the limiting value vc|u→∞ = 2 sin(πn), [33]. Note that this asymptotic value is the same for both n = 1/4, 3/4. By comparing vc|u=0 = vF = 2 sin(πn/2) with vc|u→∞, one concludes that vc(u) is maximal at a finite u for 2/3 < n < 1. The susceptibility χc = ∂µn in turn can be expressed by ǫ ′ (k)|k=B , ρ(B) and a related function [32]. Although the
Figure 2. Luttinger parameter Kc at fillings n = 1/4, 1/2, 3/4. The blue dashed lines are the low-u results Eq. (12). Substitute Eq. (29) into Eq. (30) to obtain σ(v) = 1 L (g (0) Q (v) + SQ(v, 0)) + g (0) Q (v) − Z B −B SQ(v, v′ )σ(v ′ )dv ′ (39) SQ(v, v′ ) := a2(v − v ′ ) − Z Q −Q a1(v − sin k)a1(sin k − v ′ ) cos k dk. Here g (0) Q is the ν = 0-case of g (ν) Q , defined by g (ν) Q (v) := Z Q −Q 1 π a1(v − sin k) cosν k dk. (40) Furthermore, σ (ν) (v) := 1 L (g (ν) Q (v) + SQ(v, 0)) + g (ν) Q (v) − Z ∞ −∞ SQ(v, v′ )σ (ν) (v ′ )dv ′ . (41) Eq
Figure 3. The susceptibility χbulk over u at densities n = 0.2, 0.3, 0.4, 0.5, 0.7, 1 (from top to bottom at u = 0). A similar figure has been shown by Shiba [9]. Here, we additionally compare with the field-theory result at small u (A.15) (blue dashed lines) and with the XXX-limit (for better comparison the inset shows the same figure on a larger scale; the n = 1-line is printed bold here). that this way of obtaining χ0 is equivalent to calculating vs from Eq. (34) and then using Eq. (15). The scale h0 as a function of u at different fillings is depicted in
Figure 4. The scale h0 over u for n = 0.2, 0.3, 0.4, 0.5, 0.7, 1.0 (from bottom to top). In a), the large-coupling result Eq. (82) is verified; in b), the small-coupling formula Eq. (87) is visualised (dashed lines). 10-6 10-5 10-4 10-3 10-2 10-1 100 h 0 0.4 0.8 1.2 1.6 χbulk(h) 10-7 10-6 10-5 10-4 10-3 10-2 h 0 2 4 6 8 10 χbulk(h) a) b)
Figure 5. The bulk susceptibility for n = 0.2 in a) and for n = 1 in b) at u = 1, 2, 10 (from bottom to top). The dashed blue curves are the analytical result (69). The diamonds indicate χ0, showing the decrease of the scale h0 with increasing u
Figure 6. The boundary susceptibility χB(h) for u = 1 in a) and for u = 10 in b) with n = 0.2, 1 (from top to bottom). The dots are the numerical data, the lines the asymptotic form (17) where h (B) 0 has been determined by a fit. The insets show χB · h · ln2 h, the red horizontal line denotes the asymptotic value 1/4
Figure 7. TMRG data (black circles) for the local magnetisation hs z (x)i where u = 1.0, T = 0.131, s = 0.037 and n = 0.886. In comparison the field theory result (93) is shown (blue squares) where the amplitudes and phases have been determined by a fit. parameters. The agreement is very good. In particular, the exponential decay of the correlation function is correctly described by the exponents and velocities obtained by Bethe ansatz. The surprisingly rich structure of hs z (x)i is caused by a competition between the first three terms in (93) which oscillate with different wave vectors but…
Figure 8. TMRG data (black circles) for the local susceptibility χ(x) where u = 1.0, T = 0.131, s = 0.0 and n = 0.886. In comparison the field theory result (95) is shown (blue squares) where the amplitudes and phases have been determined by a fit. important and the bulk susceptibility becomes χbulk = 1 2πvs  1 − 1 2 ln E − ln | ln E| ln2 E  . For the boundary part we even find that the result for T = 0, h ≪ h0 χB = 1 4h ln2 h  1 + ln | ln h| ln h  as well as the result for h = 0, T ≪ T0 χB = − 1 12T ln T  1 + ln | ln T | 2 ln T  are completely universal. In pa

Figures reproduced from arXiv:cond-mat/0602511 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.