{"id":"e53a18f1-5b2d-4f5d-875e-557600750f0a","arxiv_id":"2501.16631","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Classical atomic polarizability is recalculated for spheroidally distorted and non-uniform electron clouds, with parameters adjusted to match hydrogen's quantum value.","lead":"This paper derives updated classical formulas for how an atom's electron cloud distorts in an electric field, modeling the cloud as a stretched spheroid or as a mix of dense and thin charge regions. The work is mainly educational, offering intuitive corrections to the textbook atomic polarizability model.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The non-uniform model's claimed (1+βγ^3) enhancement is an electrostatics artifact: the dense core is omitted from the restoring field in Eq. (10) but included in the dipole moment p=qd, so a consistent calculation gives α/(4πε0)=a0^3.","rationale":"Sec. IV's non-uniform model is the source of the numerical match to hydrogen (β=8, γ=0.76), and Eq. (12) multiplies the distortion result by the same factor. The error above removes that factor: a consistent treatment of the core either adds to the restoring field or cancels from the dipole. Either way the non-uniform model does not improve on the simple cloud, so the central claim that non-uniform density yields improved polarizability collapses. The distortion model itself survives: Eq. (7)'s g(ε) matches the standard prolate-spheroid depolarization factor, so that part is internally consistent, though it is a known result. The reader's weakest assumption identified the rigid-core assumption; our concern sharpens it to a quantitative contradiction between Eq. (10) and p=qd. This is not an aesthetic objection: it is an internal inconsistency in the derivation of Eq. (11).","tokens_in":3272,"tokens_out":10691,"duration_ms":115809,"concrete_test":"Recompute p_ind = q r_nucleus + ∫ r ρ_e(r)dV for the Sec. IV configuration under the two possible assumptions for the dense sphere: (i) centered on the displaced nucleus; (ii) fixed at the original center. In both cases derive E_e and p_ind separately; if the calculation yields α/(4πε0)=a0^3 rather than (1+βγ^3)a0^3, Eq. (11) is invalid.","verdict_should_be":"REJECT","load_bearing_attack":"In Sec. IV the field at the displaced nucleus is written as E_e=ρV d/(3ε0) (Eq. (10)). This is valid only if the dense inner sphere (ρN, radius R_N=γa0) is centered on the nucleus and moves with it; if it stayed at the original center, its field would enter Eq. (10). Under that same assumption, the total induced dipole is not qd. Taking the nucleus at x=d and the dense sphere also centered at x=d, the nuclear contribution +qd is partly cancelled by the dense negative sphere, whose contribution is −ρN(4π/3)R_N^3 d. The background ρV sphere centered at x=0 contributes zero. Hence p_ind = [q − ρN(4π/3)R_N^3]d = ρV(4π/3)a0^3 d. Combining this with Eq. (10) gives α/(4πε0)=a0^3, not the (1+βγ^3)a0^3 of Eq. (11). If instead the dense sphere is held fixed at the original center, its field must be added to Eq. (10) and, for d<R_N, the same a0^3 result is recovered. Thus the claimed improvement is an artifact of counting the core charge in p=qd while excluding it from the restoring field, and Eq. (12) inherits the error.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":3627,"tokens_out":25483,"duration_ms":232343,"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":[{"comment":"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.","section":"III, Eqs. (5)-(7)"},{"comment":"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.","section":"IV, Eqs. (9)-(11)"},{"comment":"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.","section":"IV and Conclusion"}],"minor_comments":[{"comment":"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.","section":"IV"},{"comment":"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.","section":"V"},{"comment":"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.","section":"IV"},{"comment":"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.","section":"Abstract and Introduction"}],"recommendation":"reject","confidential_remarks":"The manuscript is clearly written and the author is transparent about the fitted parameters, but the two main derivations contain load-bearing errors. The distortion model's Eq. (5) omits the field-direction factor, and the non-uniform model counts the dense core inconsistently; when treated self-consistently, the non-uniform enhancement vanishes. These are not presentation issues and cannot be repaired by local rewriting without changing the main results. I therefore recommend rejection in the current form. A substantially revised manuscript that redoes the derivations and restates the claims as a corrected classical model could be considered, but the corrected results would no longer support the abstract's central claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","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.\n\nThe 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.\n\nCredit 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.\n\nThe 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.","headline":"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.","tokens_in":4114,"tokens_out":12940,"would_cite":false,"duration_ms":121301,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Classical atomic polarizability improves when the electron cloud is allowed to deform into a spheroid or carry a denser core.","keywords":["atomic polarizability","classical electrodynamics","spheroidal electron cloud","non-uniform charge density","hydrogen atom","Gauss's law","electric dipole","Bohr radius"],"falsifier":"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.","tokens_in":3051,"feed_emoji":"⚛️","tokens_out":11241,"duration_ms":100331,"temperature":0.7,"pith_summary":"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.","feed_headline":"Two classical tweaks hit the hydrogen polarizability target","feed_subtitle":"Uniform-sphere model falls short by a factor 4.5; the paper's formulas close the gap classically.","key_machinery":"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)]$.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the simple uniform-sphere model and its polarizability $a_0^3$, the baseline this paper modifies.","marker":"1"},{"why":"Gives the standard textbook treatment plus a rigid-cloud calculation with quantum density that yields no improvement, motivating the shape and density refinements.","marker":"2"},{"why":"Provides the semi-classical hydrogen value $21a_0^3/4$ that the improved classical models are compared against.","marker":"3"},{"why":"Supplies the quantum mechanical hydrogen polarizability $9a_0^3/2$ that the non-uniform model reproduces with $\\beta=8$, $\\gamma=0.76$.","marker":"4"}],"fun_headline_variants":["Spheroid electron cloud matches quantum hydrogen polarizability","Distortion plus non-uniform charge nails polarizability","Classical model hits quantum target for hydrogen polarizability","Shape and density: classical polarizability refined","Hydrogen polarizability: beyond the rigid sphere"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Spheroid electron cloud matches quantum hydrogen polarizability","Distortion plus non-uniform charge nails polarizability","Classical model hits quantum target for hydrogen polarizability","Shape and density: classical polarizability refined","Hydrogen polarizability: beyond the rigid sphere"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000228,"raw_usage":{"total_tokens":1479,"prompt_tokens":956,"completion_tokens":523,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":572,"completion_tokens_details":{"reasoning_tokens":451}},"tokens_in":572,"tokens_out":523,"duration_ms":15212,"temperature":1.0,"reasoning_tokens":451,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T11:49:30.059822+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the simple uniform-sphere model and its polarizability $a_0^3$, the baseline this paper modifies."},{"cited_title":"The simple model underestimates the quantum result by a factor of 4.5","cited_arxiv_id":null,"evidence_quote":"Gives the standard textbook treatment plus a rigid-cloud calculation with quantum density that yields no improvement, motivating the shape and density refinements."},{"cited_title":"The modified expression for atomic polariza bility with distortion offers an improvement over the simple model","cited_arxiv_id":null,"evidence_quote":"Provides the semi-classical hydrogen value $21a_0^3/4$ that the improved classical models are compared against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the quantum mechanical hydrogen polarizability $9a_0^3/2$ that the non-uniform model reproduces with $\\beta=8$, $\\gamma=0.76$."}],"review_version":1}