Pith. sign in

REVIEW 3 major objections 4 minor 21 references

Analytic expressions for Debye functions and the heat capacity of a solid

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

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

Pith's one-line read The paper claims an integral-free polylogarithm formula for Debye functions that makes both low- and high-temperature limits immediate.

desk verdict The new Debye-function identity (3.22) is correct and useful, but the method-of-brackets derivation is asserted rather than derived—worth refereeing, not desk-rejecting. read the letter →

arxiv 1908.08667 v1 pith:PQJEXPI5 submitted 2019-08-23 math-ph math.MP

Ivan Gonzalez , Igor Kondrashuk , Victor H. Moll , Alfredo Vega This is my paper · ORCID
classification math-phmath.MP MSC 33E2033F10
keywords methodofbracketsDebyefunctionspolylogarithmheatcapacityasymptoticexpansionssolidstatephysicsincompletegammafunctionbracketseries
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

Debye functions $D_N(X)$ describe the heat capacity of a crystalline solid, but their standard form is an integral that must be evaluated numerically or expanded separately in the low- and high-temperature regimes. This paper derives a new, integral-free closed form for the extended Debye function $D_N(\alpha,X)$ using the method of brackets. The central identity writes $D_N(\alpha,X)$ as $N\Gamma(N+1)/(X^N\alpha)$ times the difference of two polylogarithm terms, one evaluated at $\alpha$ and one at $\alpha e^{-X}$ summed with finite $X$ powers. Because the second term is exponentially small in $X$, the low-temperature limit becomes immediate, and the paper shows the same formula yields the standard high-temperature expansion as well. If the identity is right, Debye-function evaluation and Debye-solid thermodynamics reduce to evaluating polylogarithms.

What carries the argument

The load-bearing object is the bracket series produced by expanding $1/(e^t-\alpha)$ and integrating termwise with the bracket rule $\langle a\rangle = \int_0^\infty x^{a-1}\,dx$. The paper evaluates this five-index series by solving the linear system obtained from setting brackets to zero, which yields four candidate series $S_1$ through $S_4$; $S_4$ is discarded as divergent and $S_3$ is treated as a large-$X$ asymptotic term. The central identity comes from $S_2$, after the inner hypergeometric sum is rewritten as an incomplete gamma function and then as a finite combination of polylogarithms using the gamma-function identity for integer $N$.

What would settle it

Evaluate the right-hand side of the claimed identity numerically for $N=3$, $\alpha=1$, and several temperatures, say $X=0.01$, $X=1$, and $X=10$, and compare with high-precision quadrature of $D_3(X) = \frac{3}{X^3}\int_0^X \frac{t^3}{e^t-1}\,dt$; a mismatch beyond machine precision at any of these points would falsify the claimed identity.

Watch

Extended reading notes

Core claim

The paper claims that the method of brackets, applied to the integral $D_N(\alpha,X) = \frac{N}{X^N}\int_0^X \frac{t^N}{e^t-\alpha}\,dt$, produces a new analytic representation of the Debye function: $D_N(\alpha,X) = \frac{N\Gamma(N+1)}{X^N\alpha}\left[\mathrm{Li}_{N+1}(\alpha) - \sum_{k=0}^N \mathrm{Li}_{N+1-k}(\alpha e^{-X})\frac{X^k}{k!}\right]$. This is presented as a new result for general nonnegative integer $N$ and positive $\alpha$, including the physical case $\alpha=1$ after taking a limit. The paper states that this formula is equivalent to an earlier expression obtained from a companion series, which recovers the known integral-free formula of the literature, but that the new form is better suited to asymptotic analysis, making both $T\to 0$ and $T\to\infty$ limits of the Debye function, internal energy, and heat capacity follow directly from polylogarithm asymptotics.

Load-bearing premise

The derivation rests on the extension of the method of brackets that assigns values to bracket series when the number of summation indices exceeds the number of brackets; the paper cites that extension rather than proving it, and if that extension is invalid the central identity is unsupported.

Editorial extensions

If this is right

  • If the identity is correct, every Debye function for integer $N$ and positive $\alpha$ is a finite combination of polylogarithms, so no numerical quadrature is needed for tabulation.
  • The low-temperature limit $D_3(u)\sim 18\zeta(4)/u^3$ follows because polylogarithms $\mathrm{Li}_n(e^{-u})$ are exponentially small, making the Debye $T^3$ law a direct consequence.
  • The high-temperature expansion is recovered by expanding the polylogarithms near $\alpha=1$, reproducing the familiar $1-\frac38 u+\frac1{20}u^2-\cdots$ series for the Debye function.
  • The representation yields closed-form expressions for internal energy and heat capacity at arbitrary temperature, not only in the two limiting regimes.
  • Because the formula is manifestly smooth in $X$ for $X>0$, it provides an analytic continuation of the Debye function away from the integral's original domain.

Reading between the lines

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

  • The identity suggests a natural recurrence in $N$: differentiating with respect to $X$ expresses $dD_N/dX$ again in polylogarithms, potentially yielding new identities among Debye functions of consecutive dimensions.
  • The same bracket-series route could be applied to Fermi-Dirac analogues ($\alpha=-1$) or to higher-order quantum-statistical integrals, provided the extension of the method used here remains valid.
  • One testable extension is to check whether the formula holds for non-integer $N$ by analytic continuation of the gamma and polylogarithm functions, which would go beyond the paper's integer-$N$ statement.
  • The finite sum over $k$ is numerically stable for large $X$, but for very small $X$ the prefactor $X^{-N}$ may amplify rounding errors; a compensated evaluation scheme would be a practical follow-up.
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 / 4 minor

Summary. The paper uses the method of brackets (MoB) to obtain analytic, integral-free expressions for the N-dimensional Debye function, working with the extension DN(α,X) = N/X^N ∫_0^X t^N/(e^t - α) dt. It reproduces a known polylogarithmic representation (3.19) and claims a new representation (3.22), and then applies these expressions to recover the standard high- and low-temperature limits of the internal energy and heat capacity in the Debye model.

Significance. The central identity (3.22) is correct: expanding 1/(e^t - α) = Σ_{m≥1} α^{m-1} e^{-mt} and using the finite incomplete-gamma formula (3.16) reproduces it directly. This gives a compact, integral-free representation for Debye functions that is useful for asymptotic analysis, and the high- and low-temperature limits in §4.1 are correctly stated. However, the derivation of (3.22) via the method of brackets is not fully presented, and two algebraic slips appear in working formulas. The paper has merit but needs a substantial revision before it is publishable.

major comments (3)
  1. [§3.1, Eqs. (3.6)–(3.10)] The transition from the five-index bracket series (3.6) to the four series S1–S4 is the load-bearing step of the paper, but it is not derived. Rule 4 in Eq. (2.5) covers only the case where the number of summation indices equals the number of brackets; here there are five indices and three brackets, and the text merely cites the extension in [11,13]. The reader is not shown the solution of the underdetermined linear system, the choice of free indices, or the regularization of the negative-integer Gamma values that occur. Because the derivation of (3.22) rests on this step, the manuscript should supply the complete MoB reduction or, more economically, verify (3.22) directly by expanding 1/(e^t - α) = Σ_{m≥1} α^{m-1} e^{-mt} and using Eq. (3.16). I have checked that this direct verification works, so the final identity is correct; the gap is in the presented derivation, not in the result.
  2. [§3.2.2, Eq. (3.21)] Equation (3.21) contains an algebraic slip that breaks the displayed transition to (3.22). The second sum in brackets is written with [α e^{-X}]^{n1}; multiplied by the prefactor N Γ(N+1)/(X^N α), this produces α^{n1-1} e^{-n1 X}, whereas the preceding line (3.20), combined with (3.16), requires α^{n1} e^{-(n1+1)X}. Replacing [α e^{-X}]^{n1} by [α e^{-X}]^{n1+1} in that sum makes the step consistent.
  3. [§4, Eq. (4.2)] Equation (4.2) is inconsistent with the cited representation (3.19). For N=3, (3.19) gives the coefficient of ζ(4)/u^3 as 18 and the coefficients of Li_4(e^u)/u^3, Li_3(e^u)/u^2, Li_2(e^u)/u, and Li_1(e^u) as 18, 18, 9, and 3 (with signs), but (4.2) uses 24, 24, 24, 12, and 4. The new expression (4.3) has the correct coefficients, so the asymptotic limits in §4.1, which are quoted correctly, should be re-derived from the corrected (4.2) or from (4.3). This is a local error, but as written (4.2) is not an analytic expression for D_3(u).
minor comments (4)
  1. [§3.1, Eq. (3.11)] The label S4 in Eq. (3.11) should be S3; S4 has already been discarded, and the truncated series described in the preceding bullet is S3.
  2. [Throughout] There are numerous typos and formatting issues, including 'ANAL YTIC' in the title, 'his corresponding power series' in Rule 1, and inconsistent capitalization of 'the method of brackets'.
  3. [§3.1, Eq. (3.20)] The convergence of the intermediate double series (3.20) at α=1 is not discussed; the n2=0 layer is individually divergent, and cancellations are essential. A brief convergence statement would clarify the status of the intermediate series.
  4. [§3, Eq. (3.2)] The domain of α in (3.2) is not specified. For α>1 and X>ln α the defining integral diverges, while the polylogarithmic expressions are analytic continuations; a sentence clarifying the validity region of (3.19) and (3.22) would be helpful.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the polylog Debye formula (3.22) follows from the defining integral by independent expansions; the cited method of brackets is a heuristic, not a hidden input.

full rationale

The derivation chain starts from the integral definition (3.2) of the extended Debye function. The method-of-brackets series (3.6) is obtained by power-series expansions of the integrand. From (3.6), the paper states that the method yields the four series S1–S4 (Eqs. 3.7–3.10) without displaying the underdetermined-system reduction; this is a rigor gap, but it is not a circular reduction. S1 is rearranged to the known polylog representation (3.19), and S2 is rearranged via the incomplete-gamma identity (3.14)/(3.16) to the new polylog formula (3.22). Neither (3.19) nor (3.22) is assumed as an input. In fact, (3.22) can be derived independently from the same integral definition by expanding 1/(e^t−α)=Σ_{m≥1} α^{m−1}e^{−mt} and using ∫_0^X t^N e^{−mt}dt = m^{−(N+1)}γ(N+1,mX) together with (3.16). No fitted parameters appear, and no quantity is predicted from a fit to itself. The citation of the method-of-brackets extension [11,13] is a methodological self-citation, but it is not load-bearing in the circularity sense: the extension is a general procedure for underdetermined bracket systems and does not encode the Debye result. The asymptotic limits in Section 4 are obtained from the polylog formulas and reproduce standard T→0 and T→∞ results. Overall, the central identity has independent mathematical content and is not equivalent to its inputs by construction; the unsupported transition (3.6)→(3.7)–(3.10) is a correctness/rigor concern, not circularity.

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

The paper introduces no free parameters or new physical entities. The central identity depends on the formal validity of the method of brackets and on analytic continuation of polylogarithms, both methodological or standard mathematical assumptions.

assumptions (3)
  • domain assumption The method of brackets rules, including the extension for non-invertible A (more sums than brackets), correctly evaluate the integral.
    Invoked in Section 3.1 to produce S1-S4 from Eq (3.6); the extension is cited to [11,13] but not proven.
  • domain assumption The expansion 1/(e^t - α) = Σ_{m≥1} α^{m-1} e^{-mt} with interchange of sum and integral is valid for all α > 0 through analytic continuation.
    Used implicitly in deriving the polylog identities; for α > 1 the series does not converge at t=0, so the identity relies on analytic continuation.
  • standard math Polylogarithm analytic continuation for arguments > 1 gives real results after cancellation of imaginary parts in Eq (3.19) and the heat capacity formulas.
    Used when evaluating Li_s(e^u) for u>0 in Section 4.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Analytic expressions for Debye functions and the heat capacity of a solid." pith.science (2026). https://pith.science/paper/PQJEXPI5

@misc{pith2026190808667,
  author       = {Pith},
  title        = {Pith review of: Analytic expressions for Debye functions and the heat capacity of a solid},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/PQJEXPI5}},
  note         = {Machine review of arXiv:1908.08667}
}
abstract

Analytic expressions for the $N$-dimensional Debye function are obtained by the method of brackets. The new expressions are suitable for the analysis of the asymptotic behavior of this function, both in the high and low temperature limits.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

21 extracted references · 21 canonical work pages

  1. [6]

    A. E. Dubinov and A. A. Dubinova. Exact integral-free exp ressions for the integral Debye functions. Tech. Phys. Let. , 34(12):999–1001, 2008

    work page 2008

  2. [1]

    Abramowitz and I

    M. Abramowitz and I. Stegun. Handbook of Mathematical Functions with Formulas, Graphs and Mathematical Tables . Dover, New York, 1972

    work page 1972

  3. [2]

    Amdeberhan, O

    T. Amdeberhan, O. Espinosa, I. Gonzalez, M. Harrison, V. Moll, and A. Straub. Ramanujan Master Theorem. The Ramanujan Journal , 29:103–120, 2012

    work page 2012

  4. [3]

    Anastasiou, E

    C. Anastasiou, E. W. N. Glover, and C. Oleari. Applicatio n of the negative-dimension ap- proach to massless scalar box integrals. Nucl. Phys. B , 565:445–467, 2000. ANALYTIC EXPRESSIONS FOR DEBYE FUNCTIONS 9

    work page 2000

  5. [4]

    Anastasiou, E

    C. Anastasiou, E. W. N. Glover, and C. Oleari. Scalar one- loop integrals using the negative- dimension approach. Nucl. Phys. B , 572:307–360, 2000

    work page 2000

  6. [5]

    P. Debye. Zur Theorie der spezifischen W¨ armen. Ann. Phys. , 39:789, 1912

    work page 1912

  7. [7]

    G. V. Dunne and I. G. Halliday. Negative dimensional inte gration. 2. Path integrals and fermionic equivalence. Phys. Lett. B , 193:247, 1987

    work page 1987

  8. [8]

    Einstein

    A. Einstein. Die Plancksche Theorie der Strahlung und di e Theorie der spezifischen w¨ arme. Ann. Phys. (Leipzig) , 22:180–190, 1907

    work page 1907

Show all 21 references
  1. [9]

    Gonzalez, K

    I. Gonzalez, K. Kohl, and V. Moll. Evaluation of entries i n Gradshteyn and Ryzhik employing the method of brackets. Scientia, 25:65–84, 2014

    2014

  2. [10]

    Gonzalez and M

    I. Gonzalez and M. Loewe. Feynman diagrams and a combina tion of the Integration by Parts (IBP) and the Integration by Fractional Expansion (IBFE) Te chniques. Physical Review D , 81:026003, 2010

    2010

  3. [11]

    Gonzalez and V

    I. Gonzalez and V. Moll. Definite integrals by the method of brackets. Part 1. Adv. Appl. Math., 45:50–73, 2010

    2010

  4. [12]

    Gonzalez, V

    I. Gonzalez, V. Moll, and I. Schmidt. Ramanujan’s Maste r Theorem applied to the evaluation of Feynman diagrams. Adv. Applied Math. , 63:214–230, 2015

    2015

  5. [13]

    Gonzalez, V

    I. Gonzalez, V. Moll, and A. Straub. The method of bracke ts. Part 2: Examples and appli- cations. In T. Amdeberhan, L. Medina, and Victor H. Moll, edi tors, Gems in Experimental Mathematics, volume 517 of Contemporary Mathematics , pages 157–172. American Mathe- matical Society, 2010

    2010

  6. [14]

    Gonzalez and I

    I. Gonzalez and I. Schmidt. Optimized negative dimensi onal integration method (NDIM) and multiloop Feynman diagram calculation. Nuclear Physics B , 769:124–173, 2007

    2007

  7. [15]

    Gonzalez and I

    I. Gonzalez and I. Schmidt. Modular application of an in tegration by fractional expansion (IBFE) method to multiloop Feynman diagrams. Phys. Rev. D , 78:086003, 2008

    2008

  8. [16]

    Gonzalez and I

    I. Gonzalez and I. Schmidt. Modular application of an in tegration by fractional expansion (IBFE) method to multiloop Feynman diagrams II. Phys. Rev. D , 79:126014, 2009

    2009

  9. [17]

    I. S. Gradshteyn and I. M. Ryzhik. Table of Integrals, Series, and Products . Edited by D. Zwillinger and V. Moll. Academic Press, New York, 8th editio n, 2015

    2015

  10. [18]

    I. G. Halliday and R. M. Ricotta. Negative dimensional i ntegrals. I. Feynman graphs. Phys. Lett. B , 193:241, 1987

    1987

  11. [19]

    K. Huang. Statistical Mechanics. Harper & Row, New York, 1975

    1975

  12. [20]

    D. A. McQuarrie. Statistical Mechanics. Harper & Row, New York, 1975

    1975

  13. [21]

    A. T. Suzuki and A. G. M. Schmidt. Massless and massive on e-loop three-point functions in negative dimensional approach. Eur. Phys. J. , C-26:125–137, 2002. Departmento de F ´ısica y Astronomia, Universidad de V alparaiso, V alparaiso, C hile E-mail address : ivan.gonzalez@uv...

    2002

Pith tools

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