Modeling Atomic Polarizability: From Charge Distortion to Non-Uniform Distributions

O. F. de Alcantara Bonfim

REVIEW 3 major objections 4 minor 5 references

Modeling Atomic Polarizability: From Charge Distortion to Non-Uniform Distributions

T0 review · 3 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash

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

Pith's one-line read Classical atomic polarizability improves when the electron cloud is allowed to deform into a spheroid or carry a denser core.

desk verdict A clean teaching note whose non-uniform model's key result is an artifact of mixing two inconsistent assumptions; the distortion model is standard electrostatics with a sign error in the derivation. read the letter →

arxiv 2501.16631 v2 pith:FAB4TXWO submitted 2025-01-28 physics.class-ph physics.atom-ph

O. F. de Alcantara Bonfim This is my paper · ORCID

classification physics.class-phphysics.atom-ph

keywords atomicpolarizabilityclassicalelectrodynamicsspheroidalelectroncloudnon-uniformchargedensityhydrogenatomGauss'slawelectricdipoleBohrradius

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 tries to improve the textbook classical model of atomic polarizability, in which an applied electric field shifts a rigid spherical electron cloud relative to the nucleus. The author replaces the rigid sphere with a uniformly charged spheroid that has the nucleus at one focus, and derives $\alpha/(4\pi\epsilon_0)=2a_0^3/[3g(\epsilon)]$, where the eccentricity factor $g(\epsilon)$ encodes how the distortion changes the restoring field at the nucleus. A second model adds a denser inner sphere to the uniform cloud and gives $\alpha/(4\pi\epsilon_0)=(1+\beta\gamma^3)a_0^3$, and the two effects combine into Eq. (12). The point of these refinements is to close part of the gap between the simple model's $a_0^3$ and the quantum hydrogen value $9a_0^3/2$ while keeping the picture purely classical and calculable.

What carries the argument

The central object is the spheroidal electron cloud, an ellipsoid of revolution, with the nucleus at one focus; its boundary is $x^2/a^2+r^2/b^2=1$ with $b=a\sqrt{1-\epsilon^2}$. The calculation slices this spheroid into thin disks perpendicular to the field and integrates the axial field of each disk at the nucleus, packaging the result into the geometric factor $g(\epsilon)=2-2/\epsilon^2+(1/\epsilon)(1-1/\epsilon^2)\ln[(1-\epsilon)/(1+\epsilon)]$. This factor carries all distortion dependence: $g(0)=2/3$ gives back the uniform-sphere field, and the polarizability formula $\alpha/(4\pi\epsilon_0)=2a_0^3/[3g(\epsilon)]$ follows once the density is normalized by $ab^2=a_0^3$. In the non-uniform model, the machinery is charge conservation together with the assumption that only the uniform background density $\rho_V$ creates a restoring field; the dense inner core rides rigidly with the nucleus. Combining the two mechanisms yields Eq. (12), $\alpha/(4\pi\epsilon_0)=(1+\beta\gamma^3)\,2a_0^3/[3g(\epsilon)]$.

What would settle it

A direct numerical integration of the field at the focus of a uniformly charged spheroid with density $\rho=3q/(4\pi ab^2)$ should reproduce $E_e=\rho d\,g(\epsilon)/(2\epsilon_0)$; if the computed field deviates from this expression, the central geometric factor is wrong and the polarizability formula falls with it.

Watch

Extended reading notes

Core claim

The central claim is that the restoring field at the nucleus depends on the shape and internal charge layout of the electron cloud, not just its radius. For a uniformly charged spheroid with the nucleus at a focus and fixed density, the field at the focus is $E_e=\rho d\,g(\epsilon)/(2\epsilon_0)$ with $g(\epsilon)$ given by Eq. (7); balancing this against the applied field yields $\alpha/(4\pi\epsilon_0)=2a_0^3/[3g(\epsilon)]$, which reduces to the rigid-sphere result at $\epsilon=0$ and grows with distortion. In the non-uniform model, a dense inner sphere of radius $\gamma a_0$ and density $\beta\rho_V$ adds charge but, being rigidly attached to the nucleus, contributes no restoring field, so $\alpha/(4\pi\epsilon_0)=(1+\beta\gamma^3)a_0^3$; choosing $\beta=8$ and $\gamma=0.76$ reproduces the quantum hydrogen value $9a_0^3/2$. The combined expression, Eq. (12), multiplies the two correction factors. The paper presents these formulas as classical refinements that identify the separate roles of geometry and charge distribution in polarization.

Load-bearing premise

The load-bearing premise is that the applied field changes the electron cloud's shape but not its internal charge density; if the field actually rearranges the charge, the restoring field computed from the fixed density is no longer valid and the polarizability formulas do not follow.

Editorial extensions

If this is right

Any spheroidal distortion with $\epsilon>0$ raises the predicted polarizability above the rigid-sphere value $a_0^3$, so geometric distortion alone can move classical estimates toward the larger quantum values.
For small distortion, $\alpha/(4\pi\epsilon_0)\approx a_0^3(1+2\epsilon^2/5)$, giving a simple one-line upgrade to the textbook formula.
In the non-uniform model, polarizability grows linearly with the fractional core charge $\beta\gamma^3$, so a modest dense inner region is an efficient way to increase $\alpha$.
Because the distortion and non-uniformity factors multiply in Eq. (12), moderate amounts of each can together reach the quantum hydrogen value without requiring the extreme $\epsilon\approx0.97$ distortion that the distortion-only model needs.

Reading between the lines

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

The paper does not explore how sensitive the factor $g(\epsilon)$ is to the assumed placement of the nucleus; applying the same disk integration to an ellipsoid with the nucleus at the center or at other interior points would show how much of the enhancement is due to the focus geometry itself.
Relaxing the rigid-core assumption in the non-uniform model, so that the dense inner sphere can deform or lag behind the nucleus, would add a restoring-field contribution from $\rho_N$ and break the clean multiplicative form $(1+\beta\gamma^3)$; this is a natural next test of the model's assumptions.
Fitting the shape parameters $\epsilon$ or $(\beta,\gamma)$ to measured or computed polarizabilities of real atoms would turn the formulas into a diagnostic for how much of an atom's response is geometric versus charge-redistributional, a use the paper leaves implicit.

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.

Desk Editor's Note

read the letter

Here's my read. The paper is a short teaching note with two toy models for hydrogen polarizability. The non-uniform model, which gives the main quantitative improvement, is internally inconsistent. Eq. (10) for the field at the nucleus includes only the background density ρV. That's valid only if the dense core moves with the nucleus, but in that case the induced dipole is qV d, not qd, and the polarizability goes back to a0^3. If the core stays fixed, then its field has to be in Eq. (10), and you get a different expression. No consistent reading yields (1+βγ^3)a0^3; the enhancement is an artifact of counting the core charge in p=qd but not in Ee. The stress-test's conclusion holds, though its second alternative value isn't right; the core point stands.

The distortion model is cleaner in spirit but has a sign problem. Eq. (5) as written uses |d-x| without a direction factor, so the integrand isn't antisymmetric around d=0 and the integral won't give the linear-in-d field of Eq. (6). The final formula matches the standard ellipsoid depolarization factor, so the result is likely correct, but the derivation needs fixing.

Credit where due: the paper is honest that its parameters are fitted, the writing is clear, and the idea of connecting distortion to depolarization is a useful classroom illustration. The ε≈0.97 required to match hydrogen is frankly acknowledged as unrealistic.

The load-bearing flaw is in Sec. IV. That's not a minor gap. For a teaching journal like Am. J. Phys. I'd still send it to a referee, because the distortion part is salvageable and the topic is well suited to that venue, but the referee should require a consistent derivation of the non-uniform model or a removal of it. I wouldn't cite this in my own work as it stands. If the author fixes the sign and clarifies the physical picture, it could become a decent pedagogical note.

Referee Report

3 major / 4 minor

Summary. The manuscript proposes two classical refinements to the textbook uniformly-charged-sphere model of atomic polarizability. In the distortion model (Section III), the electron cloud is taken to be a uniformly charged prolate spheroid with the nucleus at one focus, and the paper derives alpha/(4 pi epsilon_0) = [2/(3 g(epsilon))] a0^3, Eq. (8). In the non-uniform model (Section IV), the cloud is a uniform background plus a denser central sphere, giving alpha/(4 pi epsilon_0) = (1 + beta gamma^3) a0^3, Eq. (11), and the parameters beta = 8, gamma = 0.76 are chosen to reproduce the quantum hydrogen value 4.5 a0^3. Section V combines the two mechanisms into Eq. (12). The central claim is that these models improve the classical estimate and provide a better pedagogical account of atomic polarizability.

Significance. If the derivations were correct, the paper would offer a simple classroom extension of the classical polarizability model. The spheroid-distortion idea is physically motivated, and the manuscript is transparent about the parameter choices and about the large distortion required to match quantum results. These are real strengths in presentation. However, the quantitative claims are not supported: the distortion derivation omits the field-direction factor, and the non-uniform result is an artifact of inconsistent charge accounting. When the non-uniform model is treated self-consistently, the enhancement (1 + beta gamma^3) disappears and the model returns the unimproved result alpha/(4 pi epsilon_0) = a0^3. The agreement with the quantum value is therefore fitted rather than predictive. The corrected derivations would not support the abstract's central claim.

major comments (3)

[III, Eqs. (5)-(7)] Equation (5) is not the x-component of the field of a disk. For a disk located at x, the field at x = d points in the +x direction when d > x and in the -x direction when d < x; this requires an extra factor sign(d - x) (equivalently (d - x)/|d - x|). Without that factor, the integrand in Eq. (5) is an even function of d, so its integral over x in (-a, a) has no term linear in d and cannot yield the linear restoring field in Eq. (6). When the sign factor is restored, the uniform-sphere limit gives the standard field rho d/(3 epsilon_0), and for the spheroid the leading correction is different from the expansion of g(epsilon) in Eq. (7): the coefficient of epsilon^2 in the restoring field is -12/5, not -2/5, relative to the rho d/(3 epsilon_0) limit. Equation (8) is therefore not established by the derivation as written.
[IV, Eqs. (9)-(11)] The derivation of Eq. (11) is internally inconsistent. It uses p = q d, with q given by Eq. (9), while the restoring field in Eq. (10) is computed only from the background density rho_V, with the dense sphere omitted because it is assumed centered on the nucleus. If the dense sphere of charge -rho_N (4 pi/3) R_N^3 is centered on the nucleus at x = d, its dipole moment relative to the origin is -rho_N (4 pi/3) R_N^3 d. The total induced dipole is then p_ind = [q - rho_N (4 pi/3) R_N^3] d = rho_V (4 pi/3) a0^3 d. Combining this p_ind with Eq. (10) gives alpha/(4 pi epsilon_0) = a0^3, not Eq. (11). If instead the dense sphere is held fixed at the origin, its field must be added to Eq. (10), and again Eq. (11) is not obtained. The factor (1 + beta gamma^3) is an artifact of counting the core charge in p but excluding it from the restoring field. Equation (12) inherits this error.
[IV and Conclusion] The quantitative agreement with the quantum hydrogen value is fitted, not predicted. In Section IV the values beta = 8 and gamma = 0.76 are chosen precisely so that Eq. (11) equals 4.5 a0^3; with two free parameters and one target number, the match is guaranteed by construction. The same is true for the distortion model when epsilon is solved from Eq. (8) to reproduce the same target. The manuscript acknowledges this in the text, but the Conclusion's statement that the non-uniform model 'achieves closer agreement with quantum mechanical results' overstates the evidential value of a fitted match.

minor comments (4)

[IV] The manuscript does not specify whether the background sphere and the dense sphere are fixed in space or move with the nucleus; this ambiguity is central to the inconsistency described in the major comment and should be clarified in any revision.
[V] Equation (12) multiplies the distortion factor and the non-uniform factor without justifying that the two mechanisms act independently; a brief argument for independence, or a statement that the combination is a heuristic, would be needed.
[IV] The notation q is used for both the total electron charge and the charge appearing in p = q d without a clear sign convention; since the electron charge is negative, the text should state whether q denotes a magnitude and should use consistent signed densities.
[Abstract and Introduction] The abstract says the models 'aim to refine the classical approximation,' but the fitted nature of the parameters should be acknowledged in the abstract as well, so that readers are not left with the impression that the quantum value is obtained from first principles.

Circularity Check

2 steps flagged · score 6.0 of 10

The quantitative match to the quantum hydrogen polarizability is obtained by choosing model parameters (β, γ, and effectively ε) to reproduce that target; the second model's coefficient is a free knob and the distortion model's large eccentricity is selected by the desired answer.

fitted input called prediction [Sec. IV, Eq. (11) and the following sentence]
"imposing E = Ee, and using Eq. (9) along with the definition of polarizability, we obtain: α/4πε0 = (1 + βγ^3)a0^3. (11) ... By choosing β = 8 and γ = 0.76, we find α/4πε0 = 4.5a0^3, which reproduces the quantum mechanical result for the hydrogen atom."

Equation (11) makes α/(4πε0) equal to (1+βγ^3)a0^3. With β and γ free parameters, 1+βγ^3 = 1+8(0.76)^3 ≈ 4.5, so the advertised 'reproduction' of the quantum hydrogen value is inserted by the parameter choice. The model provides no independent determination of β or γ, so the agreement with α/(4πε0)=4.5a0^3 is an algebraic consequence of the chosen inputs, not a prediction emerging from the model.
fitted input called prediction [Sec. III, after Eq. (8)]
"As expected, for ǫ → 0, g(0) = 2/3, and the classical result for the polarizability with no distortion is recovered. ... However, significant distortion (ǫ ≈ 0.97) is required to match quantum mechanical results, which may not represent a physically realistic distortion."

Equation (8) gives α/(4πε0) = [2/(3g(ǫ))] a0^3. The quantum value 4.5a0^3 is not obtained by an independent calculation of ǫ; instead ǫ ≈ 0.97 is the value that makes 2/(3g(ǫ)) = 4.5. In the weak-field limit ǫ→0 the formula returns the original a0^3, so the claimed improvement over the simple model is conditional on imposing an eccentricity chosen to force the known target.

full rationale

The paper contains no load-bearing self-citation and no imported uniqueness theorem; the derivation is algebraic and self-contained. The central circularity is parametric: the non-uniform model's coefficient (1+βγ^3) is a free knob, and the paper turns it to β=8, γ=0.76 to reproduce the known quantum hydrogen polarizability. The distortion model similarly reaches the quantum value only by requiring ǫ≈0.97, a value selected by the target rather than derived from an independent mechanical balance; its small-ǫ limit is the old a0^3 result. Thus the headline quantitative agreements are fitted rather than predicted, which warrants a score of 6 rather than 0. A separate electrostatics inconsistency exists — depending on whether the dense core moves with the nucleus or stays fixed, a consistent treatment of Eqs. (9)-(11) gives α/(4πε0)=a0^3 rather than (1+βγ^3)a0^3 — but that is a correctness problem, not an additional circularity, and I do not count it in the score.

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

The central formulas depend on several ad hoc modeling assumptions: a uniform density that stays uniform under distortion, a fixed spheroidal shape with the nucleus at a focus, and an inner dense sphere that moves rigidly with the nucleus. The parameters beta, gamma, and epsilon are fitted to the quantum hydrogen value, so they are free parameters rather than derived constants.

free parameters (3)

beta = 8
Chosen so that (1+beta*gamma^3)=4.5 matches the quantum hydrogen polarizability in Sec. IV.
gamma = 0.76
Chosen together with beta to reproduce alpha/(4 pi epsilon0)=4.5 a0^3 for hydrogen.
eccentricity epsilon = about 0.97
Selected so the distortion model alone matches the quantum value; the paper admits this is physically unrealistic.

assumptions (4)

domain assumption The electron cloud's charge density remains uniform and is not altered by the applied field during distortion.
Stated in Sec. III and used to compute the field from the spheroidal cloud.
ad hoc to paper The distorted electron cloud takes a prolate spheroidal shape with the nucleus at one focal point.
Postulated in Sec. III; no physical derivation for this shape.
ad hoc to paper The inner dense charge region in the non-uniform model moves with the nucleus and contributes no electric field at the nucleus.
Assumed in Sec. IV; this simplifies the restoring field to come only from the outer uniform sphere.
domain assumption Classical electrostatics (Gauss's law, disk field formula) applies to the atomic electron cloud.
The entire paper relies on continuous charge distributions and classical field equations.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Modeling Atomic Polarizability: From Charge Distortion to Non-Uniform Distributions." pith.science (2026). https://pith.science/paper/FAB4TXWO

@misc{pith2026250116631,
  author       = {Pith},
  title        = {Pith review of: Modeling Atomic Polarizability: From Charge Distortion to Non-Uniform Distributions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/FAB4TXWO}},
  note         = {Machine review of arXiv:2501.16631}
}

read the original abstract

In the classical model of atomic polarizability, atomic charges are displaced by an applied electric field, assuming the electron cloud remains spherically symmetric but with its center shifted from the nucleus, thereby inducing an electric dipole. In this work, we propose that the applied electric field distorts the initially spherical electron cloud into a spheroidal shape, with the nucleus positioned at one of its focal points. We assume that the electron cloud's charge density remains uniform and is not altered by the applied field. We derive a modified expression for the polarizability that accounts for the distortion factor, represented by the eccentricity of the spheroid, and analyze an additional model incorporating non-uniform charge distribution. Furthermore, we present an expression for polarizability that combines the effects of both distortion and non-uniform charge density. These models aim to refine the classical approximation and offer deeper insights into the mechanisms underlying atomic polarizability.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

5 extracted references · 5 canonical work pages

[1]

(3) 2 For the hydrogen atom, semi-classical calculations 3 give α/4πε0 = 21 4 a3 0, while quantum mechanical calculations 4 yield α/4πε0 = 9 2a3

work page
[2]

The simple model underestimates the quantum result by a factor of 4.5

Here a0 is the Bohr radius. The simple model underestimates the quantum result by a factor of 4.5. Attempts to improve the model include using the quantum mechanica l charge density ρ(r) = q πa 3 0 e− 2r/a 0 and a rigid spherical cloud 2. However, calculations assuming d ≪ a0 yield α/4πε0 = 3 4a3 0, which is no improvement. Here, we propose a modiﬁed clas...

work page
[3]

The modiﬁed expression for atomic polariza bility with distortion oﬀers an improvement over the simple model

(8) 3 As expected, for ǫ → 0, g(0) = 2 /3, and the classical result for the polarizability with no distortion is recovered. The modiﬁed expression for atomic polariza bility with distortion oﬀers an improvement over the simple model. However, signiﬁcant dis tortion ( ǫ ≈ 0.97) is required to match quantum mechanical results, which may not rep resent a phy...

work page
[4]

(11) This expression improves upon the simple model. By choosing β = 8 and γ = 0.76, we ﬁnd α/4πε0 = 4 .5a3 0, which reproduces the quantum mechanical result for the hydrog en atom, though the parameters used may not be entirely realistic. Comparin g the two models intro- duced here, we ﬁnd that the model with a non-uniform charge distr ibution provides a...

work page
[5]

The Classical Polarizability of the Hydro gen Atom,

(12) This combined model captures the fundamental mechanisms behind the induction of an electric dipole in an atom. Any further reﬁnements would be quantita tive in nature, involving adjustments to the shape and charge distribution of the electron cloud in response to the external electric ﬁeld. VI. CONCLUSION We proposed a classical model for atomic pola...

work page 2004

Pith tools

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

[1] [1]

(3) 2 For the hydrogen atom, semi-classical calculations 3 give α/4πε0 = 21 4 a3 0, while quantum mechanical calculations 4 yield α/4πε0 = 9 2a3

work page

[2] [2]

The simple model underestimates the quantum result by a factor of 4.5

Here a0 is the Bohr radius. The simple model underestimates the quantum result by a factor of 4.5. Attempts to improve the model include using the quantum mechanica l charge density ρ(r) = q πa 3 0 e− 2r/a 0 and a rigid spherical cloud 2. However, calculations assuming d ≪ a0 yield α/4πε0 = 3 4a3 0, which is no improvement. Here, we propose a modiﬁed clas...

work page

[3] [3]

The modiﬁed expression for atomic polariza bility with distortion oﬀers an improvement over the simple model

(8) 3 As expected, for ǫ → 0, g(0) = 2 /3, and the classical result for the polarizability with no distortion is recovered. The modiﬁed expression for atomic polariza bility with distortion oﬀers an improvement over the simple model. However, signiﬁcant dis tortion ( ǫ ≈ 0.97) is required to match quantum mechanical results, which may not rep resent a phy...

work page

[4] [4]

(11) This expression improves upon the simple model. By choosing β = 8 and γ = 0.76, we ﬁnd α/4πε0 = 4 .5a3 0, which reproduces the quantum mechanical result for the hydrog en atom, though the parameters used may not be entirely realistic. Comparin g the two models intro- duced here, we ﬁnd that the model with a non-uniform charge distr ibution provides a...

work page

[5] [5]

The Classical Polarizability of the Hydro gen Atom,

(12) This combined model captures the fundamental mechanisms behind the induction of an electric dipole in an atom. Any further reﬁnements would be quantita tive in nature, involving adjustments to the shape and charge distribution of the electron cloud in response to the external electric ﬁeld. VI. CONCLUSION We proposed a classical model for atomic pola...

work page 2004