Effect of impurities on supersolid condensate: a Ginzburg-Landau approach

Alexander V. Balatsky; Elihu Abrahams

REVIEW

Effect of impurities on supersolid condensate: a Ginzburg-Landau approach

Alexander V. Balatsky, Elihu Abrahams

Superconductivity · Submitted 22 Feb 2006 · 17 references

Pith's verdict: A speculative but honest GL note with a clean wave-function core; the Tc enhancement mechanism is plausible, but the stiffness suppression in Eq. (15) is internally inconsistent with the model's own coupling once self-consistency is imposed.

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.

T0 AI review An AI read the paper and answered a fixed list of questions.
T1 Audit Automatic checks that the references and files the paper relies on exist and match.
T2 Compute The paper's code or calculations re-run with the same results.
T3 Formal statement The claims written out in a formal math language a computer can read.
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/0602530 (version 1)

Pith IDLMS7M56O, Pith's permanent record of this paper and its review

CitationAlexander V. Balatsky, Elihu Abrahams. "Effect of impurities on supersolid condensate: a Ginzburg-Landau approach." arXiv:cond-mat/0602530 (2006). https://pith.science/paper/LMS7M56O

Cite this reviewBibTeX for Pith's review of this paper

Claim paper

Why Pith says this

Highlight any text to comment on it.

Editor's note

read the letter

Dear colleague,

The one thing you should know: this is a short, openly speculative Ginzburg-Landau note with a clean formal core—the tilted wave function gives a two-order-parameter functional with a legitimate stiffness term—and a genuinely new physical idea, that 3He acts as an attractive potential for vacancies. But the headline suppression of superfluid stiffness is not consistent with the model's own coupling once you enforce self-consistency.

What is actually new: the specific disorder coupling F3 = -λ n3 |ψ_V|^2 and the resulting linear Tc enhancement. That part is well motivated: a repulsive potential for 4He atoms is attractive for vacancies, and shifting the quadratic GL coefficient is a natural way to raise Tc. The authors are also transparent about the fits: they use the unpublished Chan data to estimate λ and g, and they flag the assumption that nV is constant. The stiffness derivation from the gradient term is the most solid piece.

The soft spots, in order. First, the two linear effects are not predictions in any strong sense: λ and g are fitted to the same data the model claims to explain, so Eqs. (14) and (15) partly restate the fit. The linear functional form is a prediction, but the coefficients are not. Second, and more serious, the stress-test concern is right. The same F3 term that shifts Tc also shifts the equilibrium amplitude: minimizing the GL functional at T=0 gives δ|ψ_V|^2/|ψ_V|^2 = λ n3/a, exactly equal to δTc/Tc0. With the authors' numbers (δTc/Tc0 ≈ 0.75, g n3/nV ≈ 0.25), the product (1 - 0.25)(1 + 0.75) ≈ 1.31, so the stiffness would increase, not decrease. Eq. (15) is therefore not just an unproven assumption; it is in tension with the model's own potential term. The conclusion's caveat—"we have no evidence" that nV is constant—undersells this: the model predicts nV changes unless an additional pinning mechanism is added.

Citation pattern is fine; the relevant vacancy literature is cited and the discussion of Anderson's objection is honest. This is a paper for supersolid specialists. It deserves a serious referee, but the referee should push the authors to either derive g, add a mechanism that keeps nV fixed, or reframe the stiffness drop as an independent phenomenological assumption rather than a consequence of the GL functional.

Referee report

Significance

The wave-function construction in Eqs. (3)-(10) is clean and gives a simple, physically transparent expression for the superfluid stiffness in terms of the vacancy condensate amplitude; this part is a useful contribution. If the disorder effects were derived rather than fitted, the linear dependence on 3He concentration would be an attractive explanation of the strong impurity sensitivity seen in torsional-oscillator experiments. However, the quantitative claims are not predictions: λ is fixed by the Tc shift and g by the stiffness drop in the same dataset, and the stiffness suppression is posited rather than obtained from the model. The manuscript is therefore best read as a speculative framework whose central numerical statements still need a self-consistent derivation.

Major objections (3)

[EFFECT OF DISORDER, Eq. (14)] Equation (14) is not an independent prediction. The parameter λ/a is fixed by the same Chan et al. data it is used to explain: the authors insert δTc/Tc0 ≈ 0.75 at δn3 ≈ 10 ppm and conclude λδn3/a ≈ 0.75. Consequently the linear Tc-dependence is a restatement of the fitted input, not a falsifiable consequence of the model. A testable statement would require λ to be computed from a microscopic model or determined from a different observable.
[EFFECT OF DISORDER, Eq. (15); Conclusion] The stiffness suppression in Eq. (15) is internally inconsistent with the model's own potential term. The same coupling F3 that shifts Tc also changes the equilibrium vacancy amplitude: minimizing Eq. (12) at low T gives |ψ_V|^2 = (a + λn3)/(2b), so δnV/nV = λn3/a = δTc/Tc0. With the authors' values, this is about 0.75 at 10 ppm, not negligible. Since ρs ∝ t|ψ_V|^2, combining t = t0(1 − g n3/nV) with the self-consistent amplitude yields ρs/ρs0 ≈ (1 − g n3/nV)(1 + λn3/a) ≈ 0.75 × 1.75 ≈ 1.31, an increase that contradicts Eq. (15). The conclusion's caveat that 'we have no evidence' for constant nV is therefore not merely an empirical uncertainty; the GL functional used to derive Eq. (14) predicts a large nV increase unless an additional mechanism is introduced to pin nV.
[EFFECT OF DISORDER, t = t0(1 − g n3/nV)] The stiffness softening is asserted, not derived. The text says 'we expect' that substitutional impurities inhibit 4He exchanges and reduce t, but no calculation, estimate, or model is given for g. Because the stiffness decrease is one of the two central claims, this needs at least a microscopic argument or a self-consistent treatment of the order-parameter amplitude.

Minor objections (4)

[Introduction and later sections] There are several typos: 'The main purpose of purpose of this note', 'discusssion', 'degreees', and 'particpate' should be corrected.
[Paragraph after Eq. (13)] The text says 'the effective hopping t in the GL free energy density of Eq. (13)' but Eq. (13) is the disorder coupling; the reference should be to Eq. (12).
[Effect of disorder, Fig. 1] The data from Chan et al. [2] are described as reproduced in Fig. 1, but the numerical values δTc/Tc0 ≈ 0.75 and δρs/ρs0 ≈ 1/4 should be listed explicitly in the text so the estimates of λ and g can be checked without access to the unpublished figure.
[Footnote [15]] The authors note in the footnote that if the Anderson-Brinkman-Huse objection is correct, the existing data cannot be used to estimate λ; this undercuts the numerical agreement advertised in the main text and should be stated more prominently.

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

3He impurities increase supersolid Tc and decrease superfluid stiffness linearly in concentration.

The reading

This paper proposes a Ginzburg-Landau description of a 4He supersolid in which solid order and superfluidity are carried by two order parameters: the density modulation n_B(r) and a vacancy Bose condensate psi_V. It argues that each 3He impurity acts as an attractive potential for vacancies, creating local patches of higher transition temperature, while the same defects impede vacancy hopping and therefore reduce the global superfluid stiffness. The result is a linear rise in Tc and a linear drop in stiffness for small 3He concentration. If correct, the model accounts for the striking experimental observation that about 10 ppm of 3He nearly doubles the supersolid onset temperature while cutting the stiffness by about a quarter.

What carries the argument

The machinery is a two-order-parameter Ginzburg-Landau free energy density, F(r) = [a(T/Tc0 - 1) - u(r)/$\Omega$] |psi_V|^2 + (t/a0)|psi_V|^2 |grad $\alpha$|^2 + b |psi_V|^4, where psi_V is the vacancy condensate field and n_B is the solid density. The vacancy field is constructed from a hard-core boson (spin-1/2) wave function rotated at every site, which produces psi_V and shows that a potential repulsive to 4He atoms is attractive to vacancies. Disorder from 3He enters as an annealed random potential -$\lambda$ n_3 |psi_V|^2 in the quadratic term (the source of the Tc shift) and as a suppression of the hopping t (the source of the stiffness drop).

Extended reading notes

Core claim

The central claim is that adding 3He impurities to a vacancy-based 4He supersolid changes the bulk superfluid properties in a simple, linear way because a 3He atom is simultaneously an attractor for vacancies and a barrier to their motion. The paper encodes both effects in a single Ginzburg-Landau functional with a random attractive potential -$\lambda$ n_3 |psi_V|^2 in the quadratic coefficient, which raises Tc according to Tc = (1 + $\lambda$ n_3/a)Tc0, and with a renormalized hopping t = t0(1 - g n_3/n_V), which lowers the stiffness according to rho_s = rho_s0(1 - g n_3/n_V). The functional is motivated by a microscopically derived tilted-spin wave function that yields the vacancy order parameter psi_V = (1/2) sin phi $e^{{-i alpha}}$.

Load-bearing premise

The load-bearing premise is that adding 3He does not change the average vacancy density nV, so the stiffness decrease is purely kinetic; the paper states it has no evidence for this.

Editorial extensions

If this is right

At low 3He concentration, the transition temperature should rise linearly with slope lambda/a per 3He atom, independent of the vacancy density.
The superfluid stiffness should fall linearly and would vanish at n3 = nV/g, indicating a disorder-driven superfluid-to-normal transition at sufficient 3He doping.
The coexistence of enhanced local condensate around each 3He impurity with globally suppressed phase rigidity implies a granular, puddle-like supersolid state.
Using the reported data, the model's parameters are lambda/a about 0.75 per 10 ppm and g about 25 for nV about 0.1%, implying an impurity sensitivity stronger than any known superfluid.
Because the two linear responses have opposite signs, measuring both Tc and rho_s as functions of n3 provides an internal consistency check of the model without additional fitting.
If the measured onset is not a true thermodynamic transition, only the relative Tc increase from the attractive-potential term survives; the stiffness suppression remains a separate prediction.

Reading between the lines

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

If the model is correct, the sharpest test is to measure the vacancy concentration independently: a rise in nV with added 3He would transfer the Tc increase from the attractive-potential mechanism to a vacancy-density effect and invalidate Eq. (15).
The same Ginzburg-Landau mechanism should apply to other light impurities in supersolid candidates, predicting a general linear enhancement of ordering by any substitutional impurity that repels host atoms.
The paper's proximity-effect remark suggests a concrete extension: treating 3He-doped regions as islands of higher Tc in a lower-Tc matrix would produce percolation-like onset behavior as a function of n3, not merely a linear shift.
Path-integral simulations of solid 4He with a single 3He impurity could test the microscopic picture directly by computing the local vacancy density and the superfluid density around the impurity.

Original abstract

Basing our arguments on a wave function that contains both positional and superfluid order, we propose a Ginzburg-Landau functional for a supersolid with the two order parameters necessary to describe such a phase: density n_B(r) and supersolid order parameter \psi_V. We argue that adding lighter $^3$He atoms to a 4He supersolid produces attractive regions for vacancies, leading to patches of higher T_c. On the other hand, the supersolid stiffness decreases in this granular state with increased 3He disorder. Both effects are linear in 3He concentration.

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, 3,618 input tokens, $0.0002 reported.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.30 · human participants 0.01 · controlled comparison 0.11 · randomized assignment 0.01 · states quantitative result 0.08 · compares with prior results 0.06 · mentions code or data availability 0.01.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; show every role.

Basing our arguments on a wave function that contains both positional and superfluid order, we propose a Ginzburg-Landau functional for a supersolid with the two order parameters necessary to describe such a phase: density n_B(r) and supersolid order parameter \psi_V. We argue that adding lighter $^3$He atoms to a 4He supersolid produces attractive regions for vacancies, leading to patches of higher T_c. On the other hand, the supersolid stiffness decreases in this granular state with increased 3He disorder. Both effects are linear in 3He concentration.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#	sentence	role	the authors' own work	states the main conclusion	strength
s0	Basing our arguments on a wave function that contains both positional and superfluid order, we propose a Ginzburg-Landau functional for a supersolid with the two order parameters necessary to describe such a phase: density n_B(r) and supersolid order parameter \psi_V.	contribution 0.60	0.86	0.08	asserted 0.74
s1	We argue that adding lighter $^3$He atoms to a 4He supersolid produces attractive regions for vacancies, leading to patches of higher T_c.	result 0.88	0.76	0.34	asserted 0.65
s2	On the other hand, the supersolid stiffness decreases in this granular state with increased 3He disorder.	result 1.00	0.80	0.20	asserted 1.00
s3	Both effects are linear in 3He concentration.	result 1.00	0.81	0.11	asserted 1.00

What kind of paper the classifier reports this is: theoretical result 0.99 · evidence from argument 0.84.

What was examined, and what was not

sources read

arXiv listing abstract v1

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.** FIG. 1: Dependence of transition temperature and superfluid

**Figure 2.** FIG. 2: When included in the

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

References

Works cited (17)

[1]

The presence of lighter 3He atom defects increases the zero point ﬂuctuations locally; this repels nearby host atoms. As we have seen, a repulsion for 4He atoms is an attraction for vacancies; this increases the local vacancy density as seen in Fig. 2. This leads to an increase in Tc. To ﬁrst approximation, the eﬀect would be linear in the 3He concentrati...

View on Pith
[2]

M.H.W. Chan et al , unpublished

View on Pith
[3]

E. Kim and M.H.W. Chan, Nature 427, 225 (2004); Sci- ence 305, 1941 (2005)

View on Pith 2004
[4]

L. Reatto, Phys. Rev. 183, 334 (1969); G.V. Chester and L. Reatto, Phys. Rev. 155, 88 (1967)

View on Pith 1969
[5]

A.F. Andreev and I.M. Lifshitz, Sov. Phys. JETP 29, 1107 (1969)

View on Pith 1969
[6]

A.J. Leggett, Phys. Rev. Lett. 25, 2543 (1970)

View on Pith 1970
[7]

G.V. Chester, Phys. Rev. A 2, 256 (1970)

View on Pith 1970
[8]

P.W. Anderson, cond-mat/0504731

View on Pith Pith review arXiv
[9]

have discussed the problem. Recent developments have been reported at a KITP workshop and are available online [10] The main purpose of purpose of this note is to de- velop an approach that allows a discusssion of the eﬀect of 3He impurities on a 4He supersolid. In what follows, we show how the approach leads to a Ginzburg-Landau functional for the supers...

View on Pith
[10]

P.W. Anderson, Basic Notions of Condensed Matter Physics, Ch. 4 pp. 143 et seq., Benjamin, Menlo Park CA (1984)

View on Pith 1984
[11]

P.W. Anderson, W.F. Brinkman and D.A. Huse, Science 310, 1164 (2005)

View on Pith 2005
[12]

Online at http://online.itp.ucsb.edu/online/smatter_m06/

View on Pith
[13]

Furthermore, we believe that the activation energy for in- terstitials is larger than that for vacancies and the band- width of the former is likely to be smaller than that of the latter

View on Pith
[14]

This representation has been used for discussion of liq - uid 4He in a number of works. See, for example, R.T. Whitlock and P.R. Zilsel, Phys. Rev. 131, 2409 (1963); P.R. Zilsel, Phys. Rev. Lett. 15, 476 (1965). A quantum lattice gas model also has a history in the literature. In particular, Liu and Fisher, J. Low Temp. Phys 10, 655 (1973) discussed super...

View on Pith 1963
[15]

We are aware that simple nearest neighbor hopping may not be adequate to describe the multiparticle exchanges responsible for vacancy motion

View on Pith
[16]

The kinetic energy operator in the spin language is (−t/ 2) ∑ [σ i ·σ j − σ z iσ z j ]

View on Pith
[17]

According to Anderson, Brinkman and Huse [9], the Tc determined in the experiments is actually not a true tran- sition temperature but some vortex dynamics tempera- ture and that the true Tc cannot be determined directly from the torsional oscillator experiments. If this objecti on is correct then we cannot use existing data to estimate the λ coeﬃcient. N...

View on Pith

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.

[1] [1]

The presence of lighter 3He atom defects increases the zero point ﬂuctuations locally; this repels nearby host atoms. As we have seen, a repulsion for 4He atoms is an attraction for vacancies; this increases the local vacancy density as seen in Fig. 2. This leads to an increase in Tc. To ﬁrst approximation, the eﬀect would be linear in the 3He concentrati...

View on Pith

[2] [2]

M.H.W. Chan et al , unpublished

View on Pith

[3] [3]

E. Kim and M.H.W. Chan, Nature 427, 225 (2004); Sci- ence 305, 1941 (2005)

View on Pith 2004

[4] [4]

L. Reatto, Phys. Rev. 183, 334 (1969); G.V. Chester and L. Reatto, Phys. Rev. 155, 88 (1967)

View on Pith 1969

[5] [5]

A.F. Andreev and I.M. Lifshitz, Sov. Phys. JETP 29, 1107 (1969)

View on Pith 1969

[6] [6]

A.J. Leggett, Phys. Rev. Lett. 25, 2543 (1970)

View on Pith 1970

[7] [7]

G.V. Chester, Phys. Rev. A 2, 256 (1970)

View on Pith 1970

[8] [8]

P.W. Anderson, cond-mat/0504731

View on Pith Pith review arXiv

[9] [9]

have discussed the problem. Recent developments have been reported at a KITP workshop and are available online [10] The main purpose of purpose of this note is to de- velop an approach that allows a discusssion of the eﬀect of 3He impurities on a 4He supersolid. In what follows, we show how the approach leads to a Ginzburg-Landau functional for the supers...

View on Pith

[10] [10]

P.W. Anderson, Basic Notions of Condensed Matter Physics, Ch. 4 pp. 143 et seq., Benjamin, Menlo Park CA (1984)

View on Pith 1984

[11] [11]

P.W. Anderson, W.F. Brinkman and D.A. Huse, Science 310, 1164 (2005)

View on Pith 2005

[12] [12]

Online at http://online.itp.ucsb.edu/online/smatter_m06/

View on Pith

[13] [13]

Furthermore, we believe that the activation energy for in- terstitials is larger than that for vacancies and the band- width of the former is likely to be smaller than that of the latter

View on Pith

[14] [14]

This representation has been used for discussion of liq - uid 4He in a number of works. See, for example, R.T. Whitlock and P.R. Zilsel, Phys. Rev. 131, 2409 (1963); P.R. Zilsel, Phys. Rev. Lett. 15, 476 (1965). A quantum lattice gas model also has a history in the literature. In particular, Liu and Fisher, J. Low Temp. Phys 10, 655 (1973) discussed super...

View on Pith 1963

[15] [15]

We are aware that simple nearest neighbor hopping may not be adequate to describe the multiparticle exchanges responsible for vacancy motion

View on Pith

[16] [16]

The kinetic energy operator in the spin language is (−t/ 2) ∑ [σ i ·σ j − σ z iσ z j ]

View on Pith

[17] [17]

According to Anderson, Brinkman and Huse [9], the Tc determined in the experiments is actually not a true tran- sition temperature but some vortex dynamics tempera- ture and that the true Tc cannot be determined directly from the torsional oscillator experiments. If this objecti on is correct then we cannot use existing data to estimate the λ coeﬃcient. N...

View on Pith