{"id":"7cbc80c3-3319-4be3-b92a-64be7f2f3fc2","arxiv_id":"2412.16733","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A general closed-form expression is presented for the expansion of real and complex solid harmonics into Cartesian monomials, with coefficients listed up to ℓ=10.","lead":"The authors derive a closed-form formula for the coefficients that convert spherical harmonic Gaussian orbitals into Cartesian Gaussians, and tabulate values up to ℓ=10. They demonstrate the formula in a multipole expansion of the Hartree potential for the La2 molecule.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No internal flaw found in Eq. (29); the remaining load-bearing caveat is the single-point validation of the multipole Hartree-potential application.","rationale":"After re-deriving Eqs. (12)–(29), I found no algebraic error: the round-down upper limit σ is forced by the vanishing of the derivative when 2k > ℓ−|m|, and the parity factor M and Heaviside constraints correctly encode the delta-function conditions. Spot checks of Tables I and II against direct expansion of the defining solid harmonics (e.g., Y_5^5, Y_4^2, Y_2^1) reproduce the listed D coefficients exactly. Thus the reader's strongest_claim about Eq. (29) is sound. The CONDITIONAL verdict is justified by the application, not by the transformation formula: the multipole expansion is demonstrated at only one point, and the non-monotone errors at ℓ=18–20 cast doubt on the claim of an 'arbitrary-order convergent' expansion for Gaussian densities. The proposed multi-point convergence test would directly settle this. I therefore leave the verdict at CONDITIONAL.","tokens_in":11783,"tokens_out":22099,"duration_ms":186747,"concrete_test":"Rerun the La2 Hartree-potential comparison at a grid of evaluation points (e.g., 3, 4, 6, and 10 Å from the nearest La nucleus, both along and perpendicular to the bond), using the same density, multipole expansion, and Boys-function target; if the ℓ-truncated error does not decrease monotonically to the target precision and then grow at any of these points, the 'arbitrary-order multipole expansion' claim should be softened to 'truncated asymptotic multipole expansion.'","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (29) is derived consistently from Eq. (14) through a derivative-at-origin extraction; spot checks against known solid harmonics for ℓ=1,2,4,5 and the values in Tables I and II reproduce the exact polynomial expansions, including the imaginary-unit factors arising from (x+iy)^|m|. I therefore find no mathematical defect in the central transformation formula. The load-bearing weakness lies in the application claim: Table III and Fig. 1 test the multipole expansion of Eq. (32) at a single point, 4 Å from the nearest La nucleus, with no specification of the expansion center R or the direction of r. Moreover, the tabulated values are not monotone in ℓ (ℓ=16 error ≈5.4e-8, ℓ=18 ≈8.4e-8), which is consistent with an asymptotic expansion for a non-compact Gaussian density rather than a convergent arbitrary-order expansion. The abstract's 'arbitrary-order multipole expansion' is therefore stronger than what the single-point demonstration supports.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives a closed-form expression, Eq. (29), for the coefficients C^{ℓ|m|}_{tuv} in the Cartesian monomial expansion of unnormalized complex solid harmonics Y^{|m|}_ℓ(r), starting from Rodrigues' formula for associated Legendre functions and extracting coefficients by differentiation at the origin. Real solid-harmonic coefficients D^{ℓm}_{tuv} are then obtained by taking real or imaginary parts, and numerical values are tabulated up to ℓ=6 in the main text and ℓ=10 in the Supplementary Material, with a Fortran source code on GitHub. As an application, the coefficients are used in a multipole expansion of the Hartree potential, which is tested on the La2 molecule at a single evaluation point 4 Å from the nearest La nucleus.","tokens_in":11922,"tokens_out":11209,"duration_ms":95882,"significance":"If Eq. (29) is correct, it is a useful and general closed-form alternative to the Schlegel–Frisch transformation, with immediate applicability to electronic-structure codes that use real solid harmonic Gaussians together with Cartesian integral algorithms. The derivation is transparent, the formula is parameter-free, and the paper ships machine-checkable artifacts: Fortran source, GitHub repository, and tabulated coefficients up to ℓ=10. The Hartree-potential demonstration is illustrative rather than systematic; its main role is to show that the coefficients can be embedded in a standard multipole workflow, but the convergence claims attached to it are stronger than the single test supports.","major_comments":[{"comment":"The multipole Hartree-potential validation is limited to a single evaluation point 4 Å from the nearest La nucleus, and the manuscript does not specify the expansion center R or the direction of r relative to the La2 axis. Because a Gaussian density has infinite support, no finite sphere centered at R contains the whole charge distribution, so the multipole series in Eq. (32) is asymptotic rather than convergent; the non-monotone errors in Table III (ℓ=16 and ℓ=18 overshoot the target) are consistent with that behavior. The abstract's phrase 'arbitrary-order multipole expansion' is therefore stronger than what is demonstrated. The authors should either specify R and r, add convergence tests at several geometries and evaluation points, and discuss the asymptotic character, or temper the application claim to a formal capability to generate multipole moments up to arbitrary order.","section":"Section III, Eq. (32), Table III, Fig. 1"},{"comment":"As printed, Eq. (29) contains the factor (-1)^{(|m|-t+2q)/2} for values of |m|-t+2q that are odd, where the exponent is a half-integer and the expression is not well-defined. This is not an empty case: for example, the entry (ℓ,m,t,u,v)=(1,-1,0,1,0) from Table I corresponds to |m|-t+2q=1. The underlying complex coefficient is genuinely imaginary in such cases, so the formula should use i^{|m|-t+2q} explicitly, or the real and imaginary parts should be written out separately without fractional powers of -1. This is a formal defect in the central equation, even though the tabulated real coefficients are recovered by taking the appropriate branch.","section":"Eq. (29)"},{"comment":"The statement that Eq. (32) applies when r is external to a sphere centered at R 'which contains the whole charge distribution' is not strictly satisfiable for a Gaussian-basis density, which has non-vanishing support everywhere. This assumption is load-bearing for the claimed convergence guarantee of the multipole expansion, and the manuscript should state that the expansion is used in an asymptotic sense for exponentially decaying (but not compact) densities.","section":"Section III, text before Eq. (32)"}],"minor_comments":[{"comment":"The derivative rule uses m!/(m-n)! together with H[m-n]; for n>m the factorial of a negative integer is undefined. This is harmless only if the Heaviside factor is interpreted as forcing the term to zero before evaluating the factorial, but the convention should be stated explicitly.","section":"Eq. (11)"},{"comment":"The transition from Eq. (22) to Eq. (25) is terse: the elimination of the summation indices s and p via the Kronecker deltas is not shown, and the reader must reconstruct the intermediate algebra. A brief explanatory sentence or an intermediate display would improve accessibility.","section":"Eqs. (16)–(29)"},{"comment":"The abstract contains a grammatical slip: 'a linear combinations' should be 'linear combinations'.","section":"Abstract"},{"comment":"The illustrative sentence 'For instance: X3 3 =15x3 −45xy2' is incomplete; it should be completed with the full monomial expression or removed.","section":"Table I, caption"},{"comment":"The formula '(ℓ+2 choose 2)' for the number of CGTOs is written without the binomial notation being typeset, and the sentence ends with a stray comma and period ('2ℓ + 1,.'). This should be cleaned up.","section":"Ref. 28 footnote"}],"recommendation":"major_revision","confidential_remarks":"The central transformation formula appears to be correct, and the reproducible code and tables are strong assets. The main obstacle to acceptance is the mismatch between the strong application claim in the abstract and the single-point, unspecified-center Hartree-potential test. I recommend major revision focused on the Hartree-potential claims and on fixing the phase-factor notation in Eq. (29)."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper's Eq. (29) is the real deal. The closed-form expression for the spherical-to-Cartesian coefficients with the floor upper limit is new as far as I know, and the derivation from Rodrigues plus derivative-at-origin extraction is transparent. I spot-checked Tables I and II against known solid harmonics for ℓ=1..4 and they pass, including the imaginary-unit factors. The Fortran source and tables up to ℓ=10 are provided, which makes this reproducible without much effort. The critique of Schlegel-Frisch lists concrete problems (z-derivative applied to z-independent function, index outside sum, noninteger upper bound), and that part reads plausible.\n\nThe weak part is Section III. The Hartree potential application is more of a sketch than a benchmark. Only one evaluation point, 4 Å from the nearest La, is tested. The multipole expansion center R and direction of r are not given. The La2 density matrix and basis set details are not shipped, so the numbers in Table III cannot be reproduced independently. The convergence is real but slow (ℓ=14 for high accuracy), and the errors are not monotone in ℓ. That is consistent with truncating an asymptotic expansion of a non-compact Gaussian density, not evidence against the formula. The abstract's \"arbitrary-order multipole expansion\" is defensible as a statement about the exact expansion, but the demonstration only shows one favorable geometry.\n\nThe central math is solid; the application section is underdocumented. This belongs in the literature as a methods reference. A referee should ask for more details on the Hartree test, or soften the claims, but nothing here is load-bearing flawed.\n\nMy recommendation: send to peer review. I would cite it if I worked on GTO integral transformations.","headline":"Solid closed-form formula for solid-harmonic coefficients; the Hartree application section is a sketch, not a full benchmark.","tokens_in":12485,"tokens_out":2093,"would_cite":true,"duration_ms":17618,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["33C55","33C45"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper derives a closed-form formula for the coefficients that convert solid-harmonic Gaussian and Slater orbitals into Cartesian monomials, and demonstrates it by building an arbitrary-order multipole Hartree potential for the La$_2$…","keywords":["solid harmonics","spherical-to-Cartesian transformation","Gaussian-type orbitals","transformation coefficients","multipole expansion","Hartree potential","real solid harmonics","Slater functions"],"falsifier":"Evaluate Eq. (29) for a high angular momentum, say $\\ell=9$, $m=4$, and a nontrivial exponent triple such as $(t,u,v)=(3,2,4)$, and compare the polynomial with a direct numerical evaluation of $Y^4_9(r)$ at a point with all three coordinates nonzero; any mismatch would falsify the coefficient formula. For the Hartree application, compute the multipole potential at several off-axis distances and for a second molecule and check whether the $\\ell=14$ partial sum still matches the reference value from the direct integral.","tokens_in":11559,"feed_emoji":"🧮","tokens_out":9528,"duration_ms":76853,"temperature":0.7,"pith_summary":"Solid harmonics are the angular building blocks of atomic orbitals; quantum-chemistry programs often prefer them because fewer functions are needed, but integral routines work with Cartesian monomials. This paper derives a single closed-form expression, Eq. (29), for the coefficients $C^{\\ell|m|}_{tuv}$ in the expansion of a solid harmonic into monomials $x^t y^u z^v$ with $t+u+v=\\ell$, including the round-down upper limit $\\sigma_{\\ell|m|}=\\lfloor(\\ell-|m|)/2\\rfloor$ that earlier treatments mishandled. The authors tabulate real coefficients up to $\\ell=10$ and supply code for arbitrary $\\ell$. They then use the coefficients to evaluate the Hartree potential of the La$_2$ molecule by multipole expansion, finding that moments up to $\\ell=14$ reproduce the reference value. If the formula is correct, any spherical-harmonic basis function can be converted exactly to Cartesian form at any angular momentum, removing a long-standing bottleneck in high-$\\ell$ calculations.","feed_headline":"Exact formula converts any spherical harmonic to Cartesian monomials","feed_subtitle":"Closed coefficients for all angular momenta, tested in a multipole Hartree potential that converges at ℓ=14 for La2.","key_machinery":"The load-bearing mechanism is the derivative identity $(\\partial/\\partial x)^n (x+a)^m = \\frac{m!}{(m-n)!}(x+a)^{m-n}H[m-n]$ with $H[0]=1$, applied twice: once to move from the derivative in $z/r$ to powers of $r$ and $z$, and once to extract monomial coefficients by differentiating at the origin. Around this sit the standard derivative formula for the associated Legendre function $P^{|m|}_\\ell$, the binomial expansion of $(x+iy)^{|m|}$ and of $(x^2+y^2+z^2)^k$, and the Kronecker-delta conditions that pin each coefficient to a specific triple $(t,u,v)$. The parity factor $M_{tu|m|}$ selects monomials whose total degree in $x$ and $y$ has the same parity as $|m|$, and the Heaviside product $\\Theta_{tu|m|kq}=H[k-\\eta_{tu|m|}]H[t-2q]H[|m|-t+2q]$ with $\\eta_{tu|m|}=(t+u-|m|)/2$ keeps all summation indices inside their valid ranges.","core_discovery":"Starting from the standard derivative formula for the associated Legendre function, the paper rewrites $Y^{|m|}_\\ell(r)$ as a derivative of a polynomial in $z/r$, trades $e^{i|m|\\varphi}$ for $(x+iy)^{|m|}$, and expands $r^{2k}=(x^2+y^2+z^2)^k$ by the binomial theorem. Differentiating both sides of the monomial expansion and evaluating at the origin isolates each coefficient through Kronecker deltas. The resulting Eq. (29) expresses $C^{\\ell|m|}_{tuv}$ as a double sum over $k$ and $q$ of binomial coefficients, factorials, and Heaviside steps, with a parity factor $M_{tu|m|}=[1-\\mathrm{mod}(t+u-|m|,2)]\\,H[t+u-|m|]$ that forces the correct symmetry and a step function $\\Theta_{tu|m|kq}$ that enforces the summation bounds. Real solid-harmonic coefficients are obtained by taking the real part for $m\\ge0$ and the imaginary part for $m<0$. The paper's central claim is that this formula is exact for all $\\ell$ and $m$, with the floor operation on $\\sigma_{\\ell|m|}$ resolving the noninteger upper bound that plagued an earlier transformation formula.","pith_inferences":["The coefficient formula is purely algebraic and exact, so it should extend to analytic derivatives of the basis functions; forces and stress tensors for spherical-basis calculations could be generated by differentiating the same sums.","The multipole validation is limited to one molecule at one evaluation point 4 Å from the nearest nucleus; testing additional distances, bond lengths, and molecular shapes would show whether $\\ell=14$ is a general rule or a property of this particular density.","Because the coefficients are rational for integer exponents, they could be precomputed once and stored, making the spherical-to-Cartesian conversion essentially free inside an integral code.","Fast multipole methods in periodic electronic-structure codes need solid-harmonic moments of Gaussian densities at high $\\ell$; this closed form supplies those moments directly and may simplify their implementation."],"forward_implications":["Any spherical-harmonic Gaussian or Slater function can be rewritten exactly as a finite sum of Cartesian Gaussians, so Cartesian integral engines can process spherical bases without approximation at any $\\ell$.","Real solid-harmonic coefficients follow directly from real and imaginary parts, enabling real-algebra implementations that avoid complex arithmetic.","The Hartree potential can be expanded to arbitrarily high multipole order with these coefficients; in the tested La$_2$ case the series reaches the reference value with moments through $\\ell=14$.","The same transformation coefficients apply unchanged to Slater functions, since only the radial factor is replaced.","The floor operation in Eq. (29) removes the noninteger summation bound of an earlier transformation formula, giving a well-defined prescription for every $\\ell$ and $m$."],"supporting_citations":[{"why":"Supplies the standard derivative formula for the associated Legendre function, the starting identity for the whole derivation.","marker":"Ref. 12"},{"why":"Gives the definition and properties of solid harmonics that set up the spherical-to-Cartesian expansion.","marker":"Ref. 11"},{"why":"An earlier transformation formula that this paper revisits and whose noninteger summation bound it fixes.","marker":"Ref. 32"},{"why":"The semi-analytical reference method used as the target for validating the multipole Hartree potential.","marker":"Ref. 19"},{"why":"Establishes the multipole expansion of the Coulomb potential for Gaussian charge densities, which the Hartree-potential application relies on.","marker":"Ref. 5"}],"fun_headline_variants":["Exact closed formula for spherical-to-Cartesian harmonic transformation","All-ℓ exact formula maps spherical harmonics to Cartesian monomials","Exact solid-harmonic to Cartesian conversion, any ℓ","One exact formula for all spherical harmonic orders to Cartesian","Closed coefficients: spherical harmonics to Cartesian monomials exactly"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The broad practical claim that the method produces accurate arbitrary-order multipole Hartree potentials for general molecular densities rests on a single validation point, 4 Å from the nearest La nucleus in La$_2$; if the multipole series converges differently for other geometries, the application claim would exceed what is shown.","fun_headline_variants_meta":{"raw":{"variants":["Exact closed formula for spherical-to-Cartesian harmonic transformation","All-ℓ exact formula maps spherical harmonics to Cartesian monomials","Exact solid-harmonic to Cartesian conversion, any ℓ","One exact formula for all spherical harmonic orders to Cartesian","Closed coefficients: spherical harmonics to Cartesian monomials exactly"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000847,"raw_usage":{"total_tokens":3654,"prompt_tokens":883,"completion_tokens":2771,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":499,"completion_tokens_details":{"reasoning_tokens":2687}},"tokens_in":499,"tokens_out":2771,"duration_ms":17440,"temperature":1.0,"reasoning_tokens":2687,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T10:15:35.476685+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate Eq. (29) for a high angular momentum, say $\\ell=9$, $m=4$, and a nontrivial exponent triple such as $(t,u,v)=(3,2,4)$, and compare the polynomial with a direct numerical evaluation of $Y^4_9(r)$ at a point with all three coordinates nonzero; any mismatch would falsify the coefficient formula. For the Hartree application, compute the multipole potential at several off-axis distances and for a second molecule and check whether the $\\ell=14$ partial sum still matches the reference value from the direct integral.","supporting_citations":[],"review_version":1}