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
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.
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
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.
Referee Report
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)
- [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.
- [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)
- 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
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
-
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
-
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
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
free parameters (2)
- plasma frequency ω_p
- relaxation time τ
axioms (2)
- domain assumption The metal obeys the Drude model for conduction electrons
- domain assumption The conductor has linear dielectric and magnetic properties (ε > ε0, μ > μ0) due to bound charges and currents
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
for such a conductor with a high carrier concentration (ω_p τ ≫1), we have obtained a simple form of the attenuation constant k_- ≃ +√(μϵ/2) √(ω_p² - ω² + |ω_p² - ω²|)
-
IndisputableMonolith/Foundation/DimensionForcing.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We have obtained the critical exponents for these quantities.
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
-
[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]
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]
(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]
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]
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]
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]
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]
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]
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]
is the chemical po- tential of the conduction electrons at a low temperature (kBT ≪ ǫF ). Eqn. (
-
[11]
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]
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]
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]
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...
work page 2023
-
[15]
D. J. Griffiths, Introduction to Electrodynamics, 4th ed., Sec. 9.4.1, pp. 412-415 (Pearson, Boston, 2013)
work page 2013
-
[16]
M. Born and E. Wolf, Principles of Optics , 7th ed., Ch. 14, pp. 735-789 (Cambridge University Press, Cam- bridge, 1999)
work page 1999
-
[17]
Roberts, Optical properties of copper , Phys
S. Roberts, Optical properties of copper , Phys. Rev. 118, 1509 (1960)
work page 1960
- [18]
-
[19]
Wooten, Optical Properties of Solids , Ch
F. Wooten, Optical Properties of Solids , Ch. 2-4, pp. 15- 84 (Academic Press, New York, 1972)
work page 1972
-
[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)
work page 2001
-
[21]
Y. Murakami, D. Goleˇ z, M. Eckstein, and P. Werner, Photoinduced nonequilibrium states in Mott insulators , Rev. Mod. Phys. 97, 035001 (2025)
work page 2025
-
[22]
M. Cardona and M. L. W. Thewalt, Isotope effects on the optical spectra of semiconductors , Rev. Mod. Phys. 77, 1173 (2005)
work page 2005
-
[23]
T. P. Devereaux and R. Hackl, Inelastic light scattering from correlated electrons, Rev. Mod. Phys. 79, 175 (2007)
work page 2007
-
[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)
work page 2017
-
[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)
work page 2025
-
[26]
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
work page 2025
-
[27]
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)
work page 2025
-
[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)
work page 2025
- [29]
-
[30]
T. Huang and Z. Sun, Universal phase transitions of mat- ter in optically driven cavities , Phys. Rev. Lett. 136, 036901 (2026)
work page 2026
-
[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)
work page 2018
-
[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)
work page 1900
-
[33]
J. D. Jackson, Classical Electrodynamics, 3rd ed., Sec. 7.5, pp. 309-316, (Wiley, New York, 1999)
work page 1999
-
[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)
work page 2013
- [35]
-
[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)
work page 2019
-
[37]
F. Boschini, M. Zonno and A. Damascelli, Time-resolved ARPES studies of quantum materials , Rev. Mod. Phys. 96, 015003 (2024)
work page 2024
-
[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)
work page 1928
-
[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)
work page 1927
-
[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)
work page 1968
-
[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)
work page 2024
-
[42]
S. Chaturvedi and S. Biswas, Fermi-Dirac statistics, Re- son. 19, 45 (2014)
work page 2014
- [43]
-
[44]
N. W. Ashcroft and N. D. Mermin, Solid State Physics , Ch. 17, pp. 340-344 (Harcourt Asia, Singapore, 1976)
work page 1976
-
[45]
S. Biswas and D. Jana, Thermodynamics of quantum gases for the entire range of temperature , Eur. J. Phys. 33, 1527 (2012)
work page 2012
-
[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)
work page 1974
-
[47]
J. F. Annett, Superconductivity, Superfluids and Con- densates, Sec. 3.2, pp. 47-49 (Oxford University Press, Oxford, 2004)
work page 2004
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.