Pith. sign in

REVIEW 3 major objections 5 minor 17 references

A refinement of the Lorentz local field expression with impact on the Clausius-Mossotti and Lorentz-Lorenz models

T0 review · 3 major / 5 minor · reviewed 2026-08-15 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The exact local field in a dielectric is $E\,3\varepsilon_r/(2\varepsilon_r+1)$, not $E(\varepsilon_r+2)/3$, and replacing the old factor changes how dense materials' permittivity and refractive index are modeled.

desk verdict The paper's claimed 'exact' local field is a modeling choice, not a correction, but the densified-glass data comparison is worth a serious look. read the letter →

arxiv 2506.01993 v5 pith:OJCUAWP6 submitted 2025-05-18 physics.optics physics.hist-ph

Jeroen van Duivenbode , Anne-Jans Faber , Reinoud Lavrijsen This is my paper · ORCID
classification physics.opticsphysics.hist-ph
keywords localfieldsphericalcavityClausius-MossottiequationLorentz-LorenzpermittivityrefractiveindexmolecularpolarizabilityMossotticatastrophe
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper argues that the standard Lorentz local field, $E_{\rm local} = E + P/(3\varepsilon_0)$, is an approximation that drops the field a molecule's own induced dipole creates in the surrounding dielectric. When that induced surface charge is included, the local field inside the spherical cavity is $E_{\rm local} = E + P/((2\varepsilon_r+1)\varepsilon_0) = 3\varepsilon_r E/(2\varepsilon_r+1)$, the classical spherical-cavity result. Feeding this field into the standard derivation replaces the Clausius-Mossotti and Lorentz-Lorenz equations with Eqs. (10) and (12), removes the Mossotti catastrophe, and gives fits to densified-glass refractive-index data that keep molecular polarizability constant up to high density. The payoff is practical: for dense materials, recovered polarizabilities are 7–36% higher than the old formulas suggest and remain stable as density changes.

What carries the argument

The working object is the boundary-value problem of a void spherical cavity of radius $R$ in a dielectric of relative permittivity $\varepsilon_r$, with the molecule's dipole represented as a surface charge $\sigma_f = P\cos\theta$ on the cavity wall. The solution has two pieces: inside the cavity a uniform field $E_1 = -\sigma_{f0}/(\varepsilon_0(2\varepsilon_r+1))\,\hat{z}$, outside a pure dipole field $E_2 = (\sigma_{f0}/(\varepsilon_0(2\varepsilon_r+1)))\,d$, with $d$ the dipole field factor. The key step is that the dielectric outside the cavity carries an induced bound charge that weakens the dipole's internal field from its vacuum value by the factor $3/(2\varepsilon_r+1)$; this converts the familiar $1/3$ in the Lorentz field into $1/(2\varepsilon_r+1)$ and makes the local field equal to the classical field of a spherical cavity. The same machinery then converts the old Clausius-Mossotti and Lorentz-Lorenz equations into their refined forms.

What would settle it

Take a single dielectric, such as a compressed noble gas or a silica glass, and measure either $\varepsilon_r$ or $n$ over the widest accessible density range; convert each point to $N\alpha$ using both Eq. (11) and Eq. (12). If the refined model is right, the $\alpha$ recovered from Eq. (12) is flat while the old model's $\alpha$ tilts with density, and the tilt grows as $N\alpha$ approaches 3. A decisive regime is $N\alpha>3$: Clausius-Mossotti predicts $\varepsilon_r$ turning negative through an asymptote, whereas Eq. (10) predicts a finite, monotonically increasing $\varepsilon_r$; measuring a dense dielectric with $N\alpha>3$ settles the issue.

Watch

Extended reading notes

Core claim

The central claim is that the exact local field acting on a spherically modeled molecule in a linear, isotropic dielectric is not $(\varepsilon_r+2)E/3$ but $3\varepsilon_r E/(2\varepsilon_r+1)$. The mechanism is the bound charge that the molecule's induced dipole pulls onto the dielectric surface around the excluded spherical volume. Inside the sphere a $P\cos\theta$ surface-charge distribution creates a uniform field; the dielectric outside responds with bound charge that attenuates this field by the factor $3/(2\varepsilon_r+1)$, entering Eq. (7). The authors therefore identify $E_{\rm local}$ with the field in an empty spherical cavity, and from it derive $N\alpha = (\varepsilon_r-1)(2\varepsilon_r+1)/(3\varepsilon_r)$ and $N\alpha = (n^2-1)(2n^2+1)/(3n^2)$. Refractive-index data on densified silicate glasses then give a better fit and, unlike the old Lorentz-Lorenz model, recover molecular polarizabilities that do not drift with density.

Load-bearing premise

The entire derivation rests on the choice that the field a molecule's own induced dipole creates in the surrounding dielectric is subtracted out and not allowed to act back on that molecule; if that self-induced reaction field participates in polarizing the molecule, Eq. (2) and the revised polarizabilities change.

Editorial extensions

If this is right

  • For dense dielectrics, Eq. (10) replaces Clausius-Mossotti: permittivity continues to rise smoothly with $N\alpha$, and the asymptote at $N\alpha=3$ — the Mossotti catastrophe — disappears.
  • For optics, Eq. (12) replaces Lorentz-Lorenz: refractive index as a function of density matches highly densified silica with sub-percent errors where the old model diverges by tens of percent.
  • Molecular polarizabilities recovered from refractive-index or permittivity data become larger by 7–36% depending on material, and stay constant as density changes, instead of tilting against the constant-$\alpha$ hypothesis.
  • When inverting densified-glass data, the base-material polarizabilities for SiO2, TiO2, and Na2O change, so any model that uses these polarizabilities as inputs would need updated values.

Reading between the lines

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

  • A direct extension the paper does not make: applying the same corrected inversion to polar liquids would put the Debye equation on the same footing, since it inherits the original Clausius-Mossotti local field; deviations from a constant recovered dipole moment would then isolate where the spherical-molecule idealization fails.
  • Because the correction factor depends on $\varepsilon_r$, literature polarizability tables derived from dense-phase refractive indices are likely systematically low by a material-dependent amount; re-deriving them could shift derived quantities such as dispersion forces.
  • The same boundary-value mechanism should reappear in local-field estimates in nonlinear optics and in field-enhancement calculations, where the $3\varepsilon_r/(2\varepsilon_r+1)$ factor could replace the familiar $(\varepsilon_r+2)/3$ enhancement in dense media.
  • An independent check would be to compress a noble gas across a wide density range, where chemical complications are absent, and test whether Eq. (12) recovers a flat $\alpha$ as density changes.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The manuscript argues that the standard Lorentz local field E + P/(3ε0) is a rough approximation and that the exact local field is the empty-cavity field E·3εr/(2εr+1), obtained by including the attenuation of the molecule's self-field by induced charges in the surrounding dielectric. It then derives revised Clausius-Mossotti and Lorentz-Lorenz relations, states that the Mossotti catastrophe disappears, and reports that the revised relations give more constant extracted molecular polarizabilities and better fits to refractive-index data for densified glasses.

Significance. If the central claim were correct, the paper would provide a simple closed-form refinement of two classical relations and would matter for extracting polarizabilities from dense media. The algebraic derivations in Section 2 and Appendix C are clear, and the finite-element verification of the boundary-value solutions is a useful check. However, the physical identification of the local field with the empty-cavity field is a modeling choice rather than an exact statement, and the paper's own mirror-charge argument suggests the opposite choice. The empirical comparison is suggestive but does not independently justify the choice. The contribution may be salvageable as a model comparison, but the title claim of an exact refinement is not supported.

major comments (3)
  1. [Section 2, Eq. (7) and following paragraph] The paper's central move is to subtract E_induced,inside, the reaction field of the molecule's own induced dipole, when forming E_local. The mirror-charge analogy in the same paragraph states that induced fields should act on the object. If the reaction field is included in E_local, then E_local = E + P/(3ε0), and combining p = αε0E_local with P = Np and P = ε0(εr−1)E gives Nα = 3(εr−1)/(εr+2), i.e. the original Clausius-Mossotti relation (9), not Eq. (10). Therefore Eq. (2) is a bookkeeping choice, not an exact solution, and the claims that Eq. (2) is 'the exact expression' (Section 2) and that Eqs. (10) and (12) are refinements are not supported.
  2. [Section 3, Eq. (10) and Section 4, Eq. (12)] The disappearance of the Mossotti catastrophe and the reported 10–36% changes in polarizability are direct consequences of excluding the reaction field. The manuscript should either provide a physical argument for why a molecule's polarizability should not respond to its own reaction field, or reframe Eqs. (10) and (12) as an alternative local-field model. As written, the internal inconsistency between the mirror-charge principle and the subtraction in Eq. (7) is a load-bearing flaw.
  3. [Section 5, Fig. 5] The empirical argument for the model is not independent of the local-field choice. Each model is calibrated by fitting α at the pristine reference point, so the subsequent agreement of the ELF curve with densified data shows that a constant α is consistent with the ELF model; it does not by itself establish that the reaction field should be excluded. Please state this limitation and discuss whether Onsager's reaction-field theory would alter the comparison.
minor comments (5)
  1. [Section 2, paragraph beginning 'Should the field...'] The phrasing is ambiguous: the answer 'it should' appears to contradict the subtraction of that field in Eq. (7). Clarify whether the reaction field is being excluded as part of the molecule's own field.
  2. [Appendix C, first line] The note to the editor and reviewers should be removed before publication.
  3. [Throughout] Please fix typographical errors: 'apporach' in Section 4, 'polarizabity' in the Conclusion, 'MCS2' in the footnote, and the legend string 'unveri-edasymptote' in Fig. 3.
  4. [Footnote 4] The conversion statement '1 Å^3 (Gaussian units) = 4π Å^3 (SI units)' is dimensionally confusing; SI polarizability volume is simply 1e-30 m^3, and the factor 4π belongs to the relation between α in SI and in Gaussian units.
  5. [Eq. (15)] The phrase 'linear polarization' should be 'linear combination' or 'linear mixing'; please also define the vector K and αc more explicitly.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity found: the refined local-field formulas are derived from self-contained boundary-value electrostatics, and the fitted polarizabilities are outputs of model comparison, not inputs.

full rationale

The paper derives Eq. (2) from the boundary-value solution for a void spherical cavity in a linear dielectric, worked out in Appendix C from Gauss's law and Legendre potentials rather than assumed from a fitted parameter. Eq. (7) is a consequence of that boundary-value solution, and combining it with p = αε0Elocal, P = Nαε0Elocal, and P = (εr−1)ε0E yields Eq. (10) by explicit algebra; no step in that chain uses Eq. (10) as an input. The polarizability values in Table 1 and the constancy tests in Figs. 4 and 5 are obtained by applying the derived formulas to measured refractive-index data, so the fitted α values are outputs of the model, not quantities fitted to force the model's prediction. There are no self-citations carrying a load-bearing uniqueness claim. The central physical choice, whether the reaction field of the induced dipole should be added to Elocal rather than subtracted as part of the self-field, is a modeling assumption that may be contested on physical grounds and affects correctness, but it does not make the derivation circular because the paper explicitly identifies the assumption and the conclusion does not require that assumption to be re-imported as an equation. Accordingly, no circular step can be exhibited by reduction of one equation to another by construction.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

The refined equations rest on standard electrostatics plus four domain assumptions. The only fitted quantities are the base-material polarizabilities, which are not used to derive the local field formula but are needed for the empirical comparisons. No new entities are introduced. The most consequential axiom is the identification of the polarizing field with the empty-cavity field, which is a model choice rather than a theorem.

free parameters (3)
  • alpha_SiO2 = 37.0 A^3 (LL), 41.9 A^3 (ELF)
    Fitted to refractive index data from Masuno and Arndt-Hummel; central to the LL versus ELF comparison in Sections 4 and 5.
  • alpha_TiO2 = 75.1 A^3 (LL), 102.5 A^3 (ELF)
    Recovered by least-squares fitting (Eq. 16) from Arndt-Hummel refractive index data for TiO2-containing glasses.
  • alpha_Na2O = 32.2 A^3 (LL), 34.4 A^3 (ELF)
    Recovered by least-squares fitting from Na2O-containing glass data.
assumptions (5)
  • domain assumption Dielectric is linear, isotropic, and homogeneous with P = (eps_r - 1) eps_0 E
    Used to connect polarization to permittivity in Eqs. (4)-(5) and throughout Appendix C.
  • domain assumption The molecule occupies a sphere of volume 4 pi R^3 / 3 = 1/N and is modeled by a P cos(theta) surface charge
    Section 2; poor for non-spherical molecules such as CS2 and CCl4.
  • ad hoc to paper The polarizing local field equals the empty-cavity field, so the molecule's induced-charge reaction field is subtracted rather than included in E_local
    Section 2, Eq. (7); this is the load-bearing model choice behind the claim that Eq. (2) is the exact local field.
  • domain assumption Molecular polarizability alpha is independent of density and molar fraction
    The hypothesis tested in Sections 4 and 5.
  • standard math Standard boundary-value electrostatics: Laplace equation, Legendre polynomials, Gauss's law
    Appendix C, Eqs. (22)-(36), used to derive field components and the corrective factor in Eq. (7).

how reviews work

0 comments
Cite this review

Pith. "Pith review of A refinement of the Lorentz local field expression with impact on the Clausius-Mossotti and Lorentz-Lorenz models." pith.science (2026). https://pith.science/paper/OJCUAWP6

@misc{pith2026250601993,
  author       = {Pith},
  title        = {Pith review of: A refinement of the Lorentz local field expression with impact on the Clausius-Mossotti and Lorentz-Lorenz models},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/OJCUAWP6}},
  note         = {Machine review of arXiv:2506.01993}
}
read the original abstract

In the 19th century Mossotti and Clausius developed an expression linking the electrical permittivity of a dielectric to the product of molecular polarizability and number density. Lorenz and Lorentz later extended this framework to encompass the refractive index of the dielectric. These classical expressions have proven remarkably successful in describing how permittivity and refractive index vary with number density, under the assumption that molecular polarizability remains relatively constant. While these models have stood the test of time and continue to offer valuable insights, their derivation relies on an approximation of the local electric field within a spherical cavity that simulates the molecular environment, excluding the field generated by the molecule or molecules themselves. For regimes of higher number densities, such as those encountered in densified dielectrics, employing an exact solution for the local field becomes increasingly important. This refinement extends the applicability of the Clausius-Mossotti and Lorentz-Lorenz equations and leads to more accurate estimates of molecular polarizability in general.

Figures

Figures reproduced from arXiv: 2506.01993 by the authors.

Figure 1
Figure 1. compares the field factors introduced above and shows the difference is significant. As an example, for a regular value of εr = 4, Eq. (1) overestimates the local field by 50%. Section 2 will explain the origin of this differ￾ence and establish Eq. (2) as the exact expres￾sion for the local field. Sections 3 and 4 recall the original Clausius-Mossotti and Lorentz￾Lorenz models and relates them to the local field exp… view at source ↗
Figure 2
Figure 2. Field algebra diagram with the feedback loop that links all field components inside [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Comparing experimental data with three models for relative permittivity versus the product of number density and molecular polarizability. A second case study is to see how a dielec￾tric’s refractive index relates to the molecu￾lar polarizability of the constituent molecules, based on the local field models of either Eq. (1) or (2). 4 The Lorentz-Lorenz model With Maxwell’s identity for the refractive in￾dex n = √εr… view at source ↗
Figures from the paper (10 more)
Figure 4
Figure 4. Figure 4: validates, with three different molar fractions of TiO2, the assumption of linearity between the compound molecular polarizabil￾ity and that of its constituents. Though the data is affected by some measurement noise, it is clear that the αc numbers recovered using the …
Figure 5
Figure 5. Figure 5: Refractive index versus density with highly densified SiO2 measured by Lobanov[9] and others[16]. The ELF model supports the constant molecular polarizability hypothesis up to significant degrees of densification. quickly diverge from the LL1 model (red line, Eq. (11))…
Figure 6
Figure 6. Figure 6: shows how polarization P can be linked to an alternative version of uniform lo￾cal field Elocal, where the latter is then ex￾perienced by a spherical set of molecules in￾stead of the single molecule modeled in [PITH_FULL_IMAGE:figures/full_fig_p008_6.png]
Figure 7
Figure 7. Figure 7: Void sphere with σf0 cos θ surface charge imposed. See Appendix C, paragraph “Graphs” for interpreting [PITH_FULL_IMAGE:figures/full_fig_p013_7.png]
Figure 8
Figure 8. Figure 8: Dielectric sphere with σf0 cos θ surface charge imposed. [3] Rudolf Clausius. Die mechanische Be￾handlung der Electricit¨at. Springer, 2nd edition, 1879. p.94. [4] I.E. Eremin et al. System modification of the equation Lorenz-Lorentz-Clausius￾Mossotti. International Jo…
Figure 9
Figure 9. Figure 9: Void sphere with surrounding average field [PITH_FULL_IMAGE:figures/full_fig_p014_9.png]
Figure 10
Figure 10. Figure 10: Dielectric sphere with surrounding average field [PITH_FULL_IMAGE:figures/full_fig_p014_10.png]
Figure 11
Figure 11. Figure 11: Conductive sphere with surrounding average field [PITH_FULL_IMAGE:figures/full_fig_p015_11.png]
Figure 12
Figure 12. Figure 12: Generic equations for any dielectric material inside/outside. [PITH_FULL_IMAGE:figures/full_fig_p015_12.png]
Figure 13
Figure 13. Figure 13: Verification with finite element model, same situation as Fig. 9b. The color scale is also used for figures 7 - 12. [8] Charles Kittel. Introduction to Solid State Physics. John Wiley & Sons, 8th edition, 2005. p.460. [9] Sergey Lobanov, Sergio Speziale, et al. Electr…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

17 extracted references · 17 canonical work pages

  1. [4]

    0E z 3 r R ^z ^r ^n ^3 <b;1 =0 <b;2 =0 <=0 1 2

    I.E. Eremin et al. System modification of the equation Lorenz-Lorentz-Clausius- Mossotti. International Journal for Light 13 <f =D <="0E z 3 r R ^z ^r ^n ^3 <b;1 =0 <b;2 =0 <=0 1 2 "1 =1 "2 =1 source (at 'z >R): E =E" ^z V1 =V2 =0;V =!E"z E1 =E2 =0 P1 =P2 =0 D1 =D2 ="0 "E (a) in vacuum <f =D <b =!P <="0E z 3 r R ^z ^r ^n ^3 <b;1 =0 <b;2 =Pcos3" !32"r+1 <=...

    work page 2021

  2. [1]

    Arndt and W

    J. Arndt and W. Hummel. The general re- fractivity formula applied to densified sil- icate glasses. Physics and Chemistry of Minerals 15, pages 363–369, 1988

    work page 1988

  3. [2]

    1 =1 "2 =1 source (on sphere): <f =<f0 cos3 V1 = <f0

    Neil Ashcroft and David Mermin. Solid State Physics. Saunders College Publish- ing, 1st edition, 1976. p.541. 12 z 3 r R ^z ^r ^n ^3 <f =<f0 cos3 <=<f 1 2 "1 =1 "2 =1 source (on sphere): <f =<f0 cos3 V1 = <f0 "0 " 1 3 "z E1 = <f0 "0 " !1 3 " ^z P1 =0 D1 ="0 "E1 V2 = <f0 "0 " 1 3 " R3 r3 "z E2 = <f0 "0 " 1 3 "d P2 =0 D2 ="0 "E2 d = 2cos3" ^r+sin3" ^3 (r=R)...

    work page 1976

  4. [3]

    Die mechanische Be- handlung der Electricit¨ at

    Rudolf Clausius. Die mechanische Be- handlung der Electricit¨ at. Springer, 2nd edition, 1879. p.94

    work page

  5. [5]

    0E z 3 r R ^z ^r ^n ^3 <f;1 =

    Richard Feynman. The Feynman lectures on physics, Chapter 11. Inside Dielectrics, volume II. Addison-Wesley, 13 edition, 1964. 14 <f =D <="0E z 3 r R ^z ^r ^n ^3 <f;1 ="0Ecos3"3 <=<f;1 1 2 "1 =N/A "2 =1 source (at 'z >R): E =E" ^z V +V1 =!E"[1!1]"z =0 E+E1 =E"(1!1) =0 P1 =0 D1 =0 V +V2 =!E" 5 1!R3 r3 6 "z E+E2 =E"[^z+d] P2 =0 D2 ="0 "(E+E2) d = 2cos3" ^r+...

    work page 1964

  6. [6]

    Griffiths

    David J. Griffiths. Introduction to electro- dynamics. Pearson, 4th edition, 2013. Ex. 4.2, problem 4.41

    work page 2013

  7. [7]

    Classical Electrody- namics

    John David Jackson. Classical Electrody- namics. Pearson, 1st edition, 1962. p.115- 116. 15 Figure 13: Verification with finite element model, same situation as Fig. 9b. The color scale is also used for figures 7 - 12

    work page 1962

  8. [8]

    Introduction to Solid State Physics

    Charles Kittel. Introduction to Solid State Physics. John Wiley & Sons, 8th edition,

    work page

Show all 17 references
  1. [9]

    Electronic, structural, and mechanical properties of SiO 2 glass at high pressure inferred from its refractive index

    Sergey Lobanov, Sergio Speziale, et al. Electronic, structural, and mechanical properties of SiO 2 glass at high pressure inferred from its refractive index. Physical Review Letters, 128(077403), 2022

    2022

  2. [10]

    Ueber die Beziehung zwischen der Fortpflanzungs- geschwindigkeit des Lichtes und der K¨ orperdichte.Annalen der Physik , pages 641–665, 1880

    Hendrik Lorentz. Ueber die Beziehung zwischen der Fortpflanzungs- geschwindigkeit des Lichtes und der K¨ orperdichte.Annalen der Physik , pages 641–665, 1880

  3. [11]

    Experimentale og theo- retiske undersøgelser over legemers bryd- ningsforhold

    Ludvig Lorenz. Experimentale og theo- retiske undersøgelser over legemers bryd- ningsforhold. Det kongelige danske Vi- denskabernes Selskabs Skrifter 5(8) , page 203–248, 1869

  4. [12]

    On Faraday’s lines of force

    James Clerk Maxwell. On Faraday’s lines of force. Transactions of the Cambridge Philosophical Society , x(Part 1):70–71, 1858

  5. [13]

    Discussione analytica sull’influenza che l’azione di un mezzo dielettrico ha sulla distribuzione dell’electricit` a alla superficie di pi` u corpi elettrici disseminati in esso

    Ottaviano-Fabrizio Mossotti. Discussione analytica sull’influenza che l’azione di un mezzo dielettrico ha sulla distribuzione dell’electricit` a alla superficie di pi` u corpi elettrici disseminati in esso. Memorie di matematica e di fisica delle Societ` a ital- iana delle sci...

  6. [14]

    Electric Moments of Molecules in Liquids

    Lars Onsager. Electric Moments of Molecules in Liquids. Journal of the American Chemical Society , pages 1486– 1493, 1936

    1936

  7. [15]

    Peter Schwerdtfeger and Jeffey K. Nagle. 2018 Table of static dipole polarizabilities of the neutral elements in the periodic ta- ble. Molecular Physics, 2018

    2018

  8. [16]

    C.Z. Tan, J. Arndt, and H.S. Xie. Op- tical properties of densified silica glasses. Physica B, (252):28–33, 1998

    1998

  9. [17]

    Modern Electrodynam- ics

    Andrew Zangwill. Modern Electrodynam- ics. Cambridge University Press, 1st edi- tion, 2012. p.174-175. 16

    2012

Pith tools

Reviewed August 15, 2026 · model on record in the stance chip above.