pith. machine review for the scientific record. sign in

arxiv: 2605.08107 · v1 · submitted 2026-04-27 · ❄️ cond-mat.str-el · cond-mat.stat-mech· physics.optics

Criticality in optical properties of the Drude and Drude-Sommerfeld metals around the plasma frequencies for high carrier concentrations

Pith reviewed 2026-05-12 02:04 UTC · model grok-4.3

classification ❄️ cond-mat.str-el cond-mat.stat-mechphysics.optics
keywords Drude modelplasma frequencyattenuation constantcritical exponentsoptical propertiesDrude-Sommerfeld modelThomas-Fermi screening
0
0 comments X

The pith

For Drude metals with high carrier concentrations, the attenuation constant simplifies to an absolute-value expression around the plasma frequency, inducing criticality in optical properties.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper derives an analytical expression for the attenuation constant of electromagnetic waves propagating through a Drude metal that also possesses background linear permittivity and permeability. In the limit of high carrier density where the product of plasma frequency and relaxation time greatly exceeds unity, this expression reduces to a compact form involving the absolute value of the difference between plasma frequency squared and wave frequency squared. This non-analytic form produces critical behavior in the attenuation constant, group velocity, and dielectric function exactly at the plasma frequency, complete with specific critical exponents. The analysis is extended to include a quantum correction using the Drude-Sommerfeld model and Thomas-Fermi screening.

Core claim

In the Drude model for a conductor with linear dielectric and magnetic responses from bound charges and currents, and in the high carrier concentration regime ω_p τ ≫ 1, the attenuation constant takes the approximate form k_- ≃ +√(μϵ/2) √(ω_p² - ω² + |ω_p² - ω²|) for frequencies both below and above the plasma frequency ω_p. This leads to criticality near ω = ω_p, with critical exponents obtained for the attenuation constant, group velocity, and complex dielectric constant. A quantum correction to these optical properties is derived within the Drude-Sommerfeld model incorporating Thomas-Fermi screening.

What carries the argument

the absolute-value expression inside the square root for the attenuation constant k_- that creates the non-analyticity at the plasma frequency

Load-bearing premise

The conductor must have linear dielectric and magnetic properties due to bound charges and currents, and must satisfy the high-carrier-density condition ω_p τ ≫ 1 so that damping can be neglected over a wide frequency range.

What would settle it

Measure the frequency dependence of the attenuation constant or the real and imaginary parts of the dielectric function in a high-density conductor such as a dense electron gas or doped semiconductor around its plasma frequency and check whether the data follow the predicted square-root absolute-value form and the associated critical exponents.

Figures

Figures reproduced from arXiv: 2605.08107 by Bikram Keshari Behera, Rhitabrata Bhattacharyya, Shyamal Biswas.

Figure 1
Figure 1. Figure 1: FIG. 1: The solid line follows Eqn. ( [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
read the original abstract

We have analytically determined the attenuation constant of the Drude metal for the entire range of frequency ($0<\omega<\infty$) of an electromagnetic (plane) wave incident on it within a single framework of classical electrodynamics. Here, by the Drude metal, we mean an electrical conductor that obeys the Drude model for the conduction electrons. We further consider the conductor to have linear dielectric and magnetic properties (i.e. permittivity $\epsilon>\epsilon_0$ and permeability $\mu>\mu_0$) due to the bound charges and bound currents in the background. Interestingly, for such a conductor with a high carrier concentration ($\omega_p\tau\gg1$), we have obtained a simple form of the attenuation constant $k_-\simeq+\sqrt{\frac{\mu\epsilon}{2}}\sqrt{\omega_p^2-\omega^2+|\omega_p^2-\omega^2|}$ for a wide range of high frequencies below and above plasma frequency $\omega_p$. Such a result gives rise to criticality in the conductor's optical properties, such as -- the attenuation constant, group velocity, and complex dielectric constant near around $\omega=\omega_p$. We have obtained the critical exponents for these quantities. We also have obtained a quantum correction to the optical properties within the Drude-Sommerfeld model with the Thomas-Fermi screening.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript analytically derives the attenuation constant k_- for plane electromagnetic waves in a Drude metal possessing background linear permittivity ε > ε0 and permeability μ > μ0. Under the high-carrier-density limit ω_p τ ≫ 1, it obtains the simplified form k_- ≃ √(μϵ/2) √(ω_p² - ω² + |ω_p² - ω²|) valid for a wide range of frequencies below and above the plasma frequency ω_p. This absolute-value expression is used to identify non-analytic critical behavior in the attenuation constant, group velocity, and complex dielectric function near ω = ω_p, with explicit critical exponents extracted. The work also presents a quantum correction within the Drude-Sommerfeld model incorporating Thomas-Fermi screening.

Significance. If the central approximation remains valid in the claimed frequency window, the paper supplies a compact analytical expression that renders the critical non-analyticity at ω_p explicit and derives concrete exponents, offering a useful classical-electrodynamics perspective on optical response in high-density conductors. The single-framework derivation across all frequencies and the inclusion of a quantum correction are constructive elements. The overall significance is reduced by the approximation's breakdown precisely where criticality is asserted.

major comments (2)
  1. [Abstract and derivation of attenuation constant] Abstract and the derivation leading to Eq. for k_-: the simplified absolute-value form is obtained by neglecting the 1/τ damping term in the Drude dielectric function on the basis of ω_p τ ≫ 1. However, in the critical region |ω - ω_p| ≲ 1/τ the difference |ω_p² - ω²| becomes comparable to or smaller than the damping scale, so the imaginary part of ε(ω) cannot be dropped; the absolute-value expression and the associated √|ω_p - ω| critical behavior therefore do not hold inside the very window defining the criticality. The manuscript must either retain the full damped expression near ω_p or demonstrate that the claimed exponents survive when damping is restored.
  2. [Abstract] Abstract claim of validity 'for a wide range of high frequencies below and above plasma frequency ω_p': this statement is inconsistent with the fact that the critical region itself lies outside the regime where the damping-free approximation is justified. The frequency window of applicability of the simple form versus the window used for the criticality analysis requires explicit delineation, including a quantitative estimate of the width 1/τ relative to the claimed 'wide range'.
minor comments (1)
  1. The definitions of the real and imaginary parts of the wave vector (k_+ and k_-) and their relation to the complex dielectric function should be stated explicitly from Maxwell's equations at the outset, before the high-frequency approximation is introduced.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and valuable comments on the regimes of validity of our approximations. We have revised the manuscript to explicitly delineate the frequency windows, include a discussion of the full damped Drude expression near ω_p, and demonstrate that the critical exponents are recovered in the intermediate regime 1/τ ≪ |ω − ω_p| ≪ ω_p. Point-by-point responses to the major comments follow.

read point-by-point responses
  1. Referee: [Abstract and derivation of attenuation constant] Abstract and the derivation leading to Eq. for k_-: the simplified absolute-value form is obtained by neglecting the 1/τ damping term in the Drude dielectric function on the basis of ω_p τ ≫ 1. However, in the critical region |ω - ω_p| ≲ 1/τ the difference |ω_p² - ω²| becomes comparable to or smaller than the damping scale, so the imaginary part of ε(ω) cannot be dropped; the absolute-value expression and the associated √|ω_p - ω| critical behavior therefore do not hold inside the very window defining the criticality. The manuscript must either retain the full damped expression near ω_p or demonstrate that the claimed exponents survive when damping is restored.

    Authors: We agree that the damping term cannot be neglected when |ω_p² − ω²| ≲ ω/τ (i.e., |ω − ω_p| ≲ 1/τ), as the imaginary part of the Drude dielectric function then becomes comparable to the real part. However, because ω_p τ ≫ 1 by assumption, this interval is narrow. In the intermediate window 1/τ ≪ |ω − ω_p| ≪ ω_p the damping-free approximation remains valid and the √|δω| form is recovered. We have added a new subsection deriving the full damped expression for k_- near ω_p and showing analytically that the critical exponents for k_-, group velocity, and dielectric response are asymptotically restored in the stated window as τ → ∞. The abstract has been updated to state the precise applicability condition. revision: partial

  2. Referee: [Abstract] Abstract claim of validity 'for a wide range of high frequencies below and above plasma frequency ω_p': this statement is inconsistent with the fact that the critical region itself lies outside the regime where the damping-free approximation is justified. The frequency window of applicability of the simple form versus the window used for the criticality analysis requires explicit delineation, including a quantitative estimate of the width 1/τ relative to the claimed 'wide range'.

    Authors: We have revised the abstract to read: 'for a wide range of high frequencies below and above the plasma frequency ω_p, provided |ω_p² − ω²| ≫ ω/τ'. This excludes only a narrow interval of width ∼1/τ around ω_p. In the revised text we provide the quantitative estimate that the relative width of the excluded region is ∼1/(ω_p τ) ≪ 1, so that the simple form still covers a wide range around ω_p. The criticality analysis is performed in the broad intermediate window 1/τ ≪ |ω − ω_p| ≪ ω_p, which is fully consistent with the high-carrier-density limit. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation follows directly from standard Drude model under explicit approximation

full rationale

The central expression for the attenuation constant is obtained by substituting the Drude dielectric function (with background ε, μ) into the plane-wave dispersion relation and taking the high-carrier-density limit ω_p τ ≫ 1 to drop the damping term over a broad range. The absolute-value form and subsequent critical-exponent analysis are algebraic consequences of that limit applied to the real part of the wave vector; ω_p and τ are independently defined material parameters, not fitted to the optical quantities being derived. No load-bearing self-citations, no parameter fitting renamed as prediction, and no self-definitional loops appear in the chain. The approximation's validity near ω = ω_p is a separate question of correctness, not circularity.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on the classical Drude conductivity model, linear electromagnetic response of the background medium, and the high-carrier-density limit that suppresses damping; these are standard domain assumptions rather than new postulates.

free parameters (2)
  • plasma frequency ω_p
    Material parameter set by carrier density; appears explicitly in the simplified attenuation expression
  • relaxation time τ
    Material parameter; used only to define the high-density regime ω_p τ ≫ 1
axioms (2)
  • domain assumption The metal obeys the Drude model for conduction electrons
    Explicitly stated as the definition of the Drude metal under study
  • domain assumption The conductor has linear dielectric and magnetic properties (ε > ε0, μ > μ0) due to bound charges and currents
    Invoked to justify the background permittivity and permeability in the wave propagation analysis

pith-pipeline@v0.9.0 · 5569 in / 1813 out tokens · 88601 ms · 2026-05-12T02:04:15.216085+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

47 extracted references · 47 canonical work pages

  1. [1]

    and (2) inside the conductor, under the normal incidence as mentioned above, take the form ⃗E(⃗ r, t) = ⃗E0ei[⃗k·⃗ r− ωt ] (3) for the electric field and ⃗B(⃗ r, t) = ⃗B0ei[⃗k·⃗ r− ωt ] (4) for the magnetic field with the same frequency ω but a different (complex) wave vector ⃗k = (k+ + ik− )ˆk (where k+ > 0 and k− > 0) due to the Drude model AC electrical c...

  2. [2]

    I n that case, m∗ e should be replaced with the effective mass of a hole m∗ h, which can be negative according to the band theory of solids

    and ( 4)) is known as 1 For a p-type semiconductor, holes are the charge carriers. I n that case, m∗ e should be replaced with the effective mass of a hole m∗ h, which can be negative according to the band theory of solids. For this reason, in the formula of the DC conductivit y, we have written |m∗ e | rather than only m∗ e so that the same formula can al...

  3. [3]

    (or ( 4)) to differential equation ( 1) (or ( 2)) with the complex wave-vector ⃗k = ( k+ + ik− )ˆk and the complex conduc- tivity σ (as mentioned in Eqn. ( 5)) results in k2 = iµσω + µǫω 2 which yields k2 + − k2 − = ǫµω 4τ2 + (µǫ − µσ 0τ)ω 2 1 + ω 2τ2 (6) and 2k+k− = µσ 0ω 1 + ω 2τ2 (7) once the real and imaginary parts of the complex wave- number ( k = k+...

  4. [4]

    and ( 7), we get the phase constant, as k+ = √ µǫω 2 2(1 + ω 2τ2) [ √ (1 + ω 2τ2) [ 1 + (ω 2 − ω 2p)2τ2/ω 2] + [ 1 + (ω 2 − ω 2 p)τ2] ] 1/ 2 (8) and the attenuation constant, as k− = √ µǫω 2 2(1 + ω 2τ2) [ √ (1 + ω 2τ2) [ 1 + (ω 2 − ω 2p)2τ2/ω 2] − [ 1 + (ω 2 − ω 2 p)τ2] ] 1/ 2 , (9) where ω p = √ nee2 |m∗ e |ǫ [ 6, 20] is the plasma frequency of the cond...

  5. [5]

    and (9) are the unified results for the phase con- stant and the attenuation constant, respectively, as they are true for the entire range of frequencies (0 < ω < ∞ ) of the incident waves and the bound charges and currents (with ǫ > ǫ0 & µ > µ 0) in the background of the conduc- tion electrons in the Drude metal. Such a unification not only allows us to st...

  6. [6]

    A case of the high carrier concentration in the quasi-static regime The quality of an electrical conductor is decided by the value of its skin-depth ( δ) or attenuation constant (k− = 1 /δ ). The Drude metal would be good with a large skin-depth (or poor with a small skin depth) if the very low-frequency condition ωτ ≪ ω 2 pτ2 (or a low carrier concentrat...

  7. [7]

    In such a case, the attenuation constant takes the form Eqn

    A case of the high carrier concentration in the non-quasi-static regime: a criticality at ω = ω p The quality of the Drude metal with the bound charges and currents in the background has not been analysed for a high carrier concentration ( ω pτ ≫ 1) in the non-quasi- static regime ( ωτ ≳ 1), in particular, around the plasma frequency ( ω ∼ ω p). In such a...

  8. [8]

    for ω < ω p and ω pτ → ∞ , however, takes the very simple form k− = √ µǫ √ ω 2p − ω 2 which, of course, is not available in the literature. The two completely different behaviours result in a criticality in the attenuation constant k− ≃ √ 2µǫω p|1 − ω/ω p|ν (13) near around ω = ω p with the critical exponent ν = 1 / 2 for ω < ω p and ν → ∞ for ω > ω p. We ...

  9. [9]

    and ( 9)) for the phase consonant ( k+) and the attenuation constant ( k− ) also unify other optical properties of the Drude metal, in particular, skin-depth ( δ = 1 /k − ), absorption coeffi- cient ( α = 2 k− ), complex refractive index (˜ n = c ω k = c ω [k+ + ik− ]2), complex dielectric constant (˜ ǫ/ǫ 0 = ˜n2), phase speed ( vp = ω/k +), group velocity ...

  10. [10]

    is the chemical po- tential of the conduction electrons at a low temperature (kBT ≪ ǫF ). Eqn. (

  11. [11]

    ( 19), of course, is an out- come of the Sommerfeld asymptotic expansion 7 of the right-hand side of Eqn

    along with Eqn. ( 19), of course, is an out- come of the Sommerfeld asymptotic expansion 7 of the right-hand side of Eqn. ( 17). This makes the model for the optical properties as the Drude-Sommerfeld model, and the conductor, as the Drude-Sommerfeld metal. It is interesting to note that, while the 2nd term of Eqn. ( 18) is of purely quantum many-body scr...

  12. [12]

    We have determined the critical ex- ponents for these quantities

    to (16). We have determined the critical ex- ponents for these quantities. The Drude model, however, is a classical model which is applicable at a high temper- ature ( T ≫ TF ). For a low temperature ( T ≪ TF ), we have employed the Drude-Sommerfeld model to get quan- tum corrections on the optical properties, especially with the plasma frequency squared,...

  13. [13]

    The criticality we have explored is not the criticality often discussed in connection with the continuous phase transitions in statistical mechanics or thermodynamics [ 32]

    and ( 19). The criticality we have explored is not the criticality often discussed in connection with the continuous phase transitions in statistical mechanics or thermodynamics [ 32]. It is rather the criticality of the optical properties in the frequency domain around the plasma frequency. Here, we haven’t considered the permittivity and per- meability ...

  14. [14]

    This is an approximation we have made by further considering the conductor to be a linear medium

    and (2). This is an approximation we have made by further considering the conductor to be a linear medium. Such an approximation is applicable for a long time ( t ≫ ǫ/σ 0) scale and for all frequencies ω > 0. The static equilib- rium, of course, is not achieved for such an oscillatory case of ω > 0. In case of the static equilibrium ( ω = 0, and t ≫ ǫ/σ 0...

  15. [15]

    D. J. Griffiths, Introduction to Electrodynamics, 4th ed., Sec. 9.4.1, pp. 412-415 (Pearson, Boston, 2013)

  16. [16]

    Born and E

    M. Born and E. Wolf, Principles of Optics , 7th ed., Ch. 14, pp. 735-789 (Cambridge University Press, Cam- bridge, 1999)

  17. [17]

    Roberts, Optical properties of copper , Phys

    S. Roberts, Optical properties of copper , Phys. Rev. 118, 1509 (1960)

  18. [18]

    Linz, Jr

    A. Linz, Jr. and R. E. Newnham, Ultraviolet absorption spectra in ruby , Phys. Rev. 123, 500 (1961)

  19. [19]

    Wooten, Optical Properties of Solids , Ch

    F. Wooten, Optical Properties of Solids , Ch. 2-4, pp. 15- 84 (Academic Press, New York, 1972)

  20. [20]

    Fox, Optical Properties of Solids , Ch

    M. Fox, Optical Properties of Solids , Ch. 2, pp. 25-84 & Ch. 7. pp. 143-164 (Oxford University Press, Oxford, 2001)

  21. [21]

    Murakami, D

    Y. Murakami, D. Goleˇ z, M. Eckstein, and P. Werner, Photoinduced nonequilibrium states in Mott insulators , Rev. Mod. Phys. 97, 035001 (2025)

  22. [22]

    Cardona and M

    M. Cardona and M. L. W. Thewalt, Isotope effects on the optical spectra of semiconductors , Rev. Mod. Phys. 77, 1173 (2005)

  23. [23]

    T. P. Devereaux and R. Hackl, Inelastic light scattering from correlated electrons, Rev. Mod. Phys. 79, 175 (2007)

  24. [24]

    Giustino, Electron-phonon interactions from first prin- ciples, Rev

    F. Giustino, Electron-phonon interactions from first prin- ciples, Rev. Mod. Phys. 89, 015003 (2017)

  25. [25]

    E. D. Russo, T. Verstijnen, P. Koenraad, K. Pantzas, G. Patriarche, and L. Rigutti, Order and disorder at the atomic scale: Microscopy applied to semiconductors , Rev. Mod. Phys. 97, 025006 (2025)

  26. [26]

    de Vries, J

    N. de Vries, J. Chen, E. Hoglund, X. Guo, D. Chaudhuri, J. Hachtel, and P. Abbamonte, Comparative analysis of plasmon modes in layered Lindhard metals and strange metals, Phys. Rev. B 112, 165145 (2025) 7

  27. [27]

    Arafune, H

    R. Arafune, H. Ishida, C.-L. Lin, and N. Takagi, Probing moir´ e Bloch bands of photoexcited electrons on graphene/Ir(111), Phys. Rev. B 112, L161408 (2025)

  28. [28]

    C. M. Tonauer, E.-M. K¨ ock, R. Henn, C. Kappacher, C. W. Huck, and T. Loerting, Near-Infrared spectroscopic sensing of hydrogen order in ice XIII , Phys. Rev. Lett. 135, 018002 (2025)

  29. [29]

    Tiwari, R

    V. Tiwari, R. Korol, and I. Franco, Robust purely opti- cal signatures of Floquet states in laser-dressed crystals , Phys. Rev. Lett. 135, 186901 (2025)

  30. [30]

    Huang and Z

    T. Huang and Z. Sun, Universal phase transitions of mat- ter in optically driven cavities , Phys. Rev. Lett. 136, 036901 (2026)

  31. [31]

    E. J. Rothwell and M. J. Cloud, Electromagnetics, edn. 3, Sec. 3.6.2.3, pp. 253-255, CRC Press (Taylor and Francis Group), Boca Raton (2018)

  32. [32]

    Drude, Zur elektronentheorie der metalle , Ann

    P. Drude, Zur elektronentheorie der metalle , Ann. der Physik 306, 566 (1900) & Zur elektronentheorie der met- alle; II. Teil. Galvanomagnetische und thermomagnetis- che effecte , Ann. der Physik 308, 369 (1900)

  33. [33]

    J. D. Jackson, Classical Electrodynamics, 3rd ed., Sec. 7.5, pp. 309-316, (Wiley, New York, 1999)

  34. [34]

    Zangwill, Modern Electrodynamics, Sec

    A. Zangwill, Modern Electrodynamics, Sec. 14.7 pp. 472- 481 & Sec. 18.5-6, pp. 630-648, Cambridge University Press, Cambridge (2013)

  35. [35]

    S. Biswas, Revisiting the integral form of Gauss’ law for a generic case of electrodynamics with arbitrarily mov- ing Gaussian surface , arXiv:2412.13221v3 [Accepted in Resonance (2025).]

  36. [36]

    K. E. Oughstun, Electromagnetic and Optical Pulse Propagation, Vol. 1: Spectral Representations in Tempo- rally Dispersive Media Second Edition , edn. 2, Sec. 5.2.6, pp. 300-306 (Springer Nature, Cham, 2019)

  37. [37]

    Boschini, M

    F. Boschini, M. Zonno and A. Damascelli, Time-resolved ARPES studies of quantum materials , Rev. Mod. Phys. 96, 015003 (2024)

  38. [38]

    Sommerfeld, Zur Elektronentheorie der Metalle auf Grund der Fermischen Statistik , Z

    A. Sommerfeld, Zur Elektronentheorie der Metalle auf Grund der Fermischen Statistik , Z. f¨ ur Physik 47, 1 (1928)

  39. [39]

    L. H. Thomas, The calculation of atomic fields , Math. Proc. Cam. Philos. Soc. 23, 542 (1927) ; E. Fermi, Eine statistische methode zur bestimmung einiger eigen- schaften des atoms und ihre anwendung auf die theorie des periodischen systems der elemente , Z. f¨ ur Physik48, 73 (1928)

  40. [40]

    H. W. Verleur, Determination of optical constants from reflectance or transmittance measurements on bulk crys- tals or thin films , J. Opt. Soc. Am. 58, 1356 (1968)

  41. [41]

    B. K. Behera, S. K. Gour, and S. Biswas, Refractive in- dex for the mechanical refraction of a relativistic particl e, Eur. Phys. J. D 78, 60 (2024)

  42. [42]

    Chaturvedi and S

    S. Chaturvedi and S. Biswas, Fermi-Dirac statistics, Re- son. 19, 45 (2014)

  43. [43]

    Biswas, D

    S. Biswas, D. Jana, and R. K. Manna, Excess energy of an ultracold Fermi gas in a trapped geometry , Eur. Phys. J. D 66, 217 (2012)

  44. [44]

    N. W. Ashcroft and N. D. Mermin, Solid State Physics , Ch. 17, pp. 340-344 (Harcourt Asia, Singapore, 1976)

  45. [45]

    Biswas and D

    S. Biswas and D. Jana, Thermodynamics of quantum gases for the entire range of temperature , Eur. J. Phys. 33, 1527 (2012)

  46. [46]

    M. E. Fisher, The renormalization group in the theory of critical behavior, Rev. Mod. Phys. 46, 597 (1974) ; Erra- tum: Rev. Mod. Phys. 47, 543 (1975)

  47. [47]

    J. F. Annett, Superconductivity, Superfluids and Con- densates, Sec. 3.2, pp. 47-49 (Oxford University Press, Oxford, 2004)