Pith. sign in

REVIEW 3 major objections 6 minor 36 references

Thermodynamic Properties and Superstatistics of Graphene under a Constant Magnetic Field

T0 review · 3 major / 6 minor · reviewed 2026-08-07 · deepseek-v4-flash

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

Pith's one-line read For graphene in a constant magnetic field, superstatistical temperature fluctuations yield larger entropy and smoother thermodynamic response than the canonical ensemble.

desk verdict The superstatistical partition function in Eq. (52) does not follow from the defining integral; the paper's central 'additional disorder' conclusion is unsupported. read the letter →

arxiv 2506.00709 v1 pith:3JWBWMBP submitted 2025-05-31 quant-ph

Yarou M. Assimiou , Daniel S. Takou , Boukari Amidou , Guingarey Issoufou , Finagnon A. Dossa , Gabriel Y. H. Avossevou This is my paper · ORCID
classification quant-ph
keywords grapheneDirac-WeylequationLandaulevelsconstantmagneticfieldsuperstatisticsdeformedBoltzmannfactorpartitionfunctionthermodynamicproperties
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 solves the Dirac-Weyl equation for graphene in a constant magnetic field, obtaining the relativistic Landau spectrum \(E_n=\pm\hbar v_F\sqrt{2nD}\). It feeds that spectrum into a partition function and derives mean energy, specific heat, entropy, and free energy in two statistical descriptions: the ordinary canonical ensemble and superstatistics with a q-deformed Boltzmann factor. The central claim is that local fluctuations of the inverse temperature add genuine disorder, so the superstatistical entropy lies above the canonical entropy at the same temperature and the free energy is smoothed out. If the claim is right, superstatistics offers a more robust thermodynamic description of graphene under magnetic confinement, which matters for interpreting experiments on driven or non-equilibrium graphene.

What carries the argument

The central objects are the Landau-level spectrum \(E_n=\pm\hbar v_F\sqrt{2nD}\), obtained by reducing the Dirac-Weyl equation to a harmonic-oscillator form through the Nikiforov-Uvarov method, and the q-deformed Boltzmann factor \($B_E^{{(q)}}$(\$\beta$)=$e^{{-\beta E}}$\left(1+\frac{q}{2}\$beta^{2}$ $E^{2}$\right)\), which encodes superstatistical averaging over inverse-temperature fluctuations. The spectrum supplies the energy input to the partition function; the deformed Boltzmann factor, integrated over the level index, produces a closed-form superstatistical partition function \(Z_s=(2+6q)/($a^{2}$\$beta^{2}$)\), from which all thermodynamic quantities follow.

What would settle it

Recompute the superstatistical partition function with the same finite upper limit \(\$\lambda$\) used in the canonical case, \(Z_s=\int_0^\$\lambda$ $e^{{-\beta E}}$(1+\frac{q}{2}\$beta^{2}$ $E^{2}$)\,dn\), and evaluate \(S_s\) at \(q=0.5\). If the superstatistical entropy no longer exceeds the canonical entropy \(S(\$\beta$)\), then the reported additional disorder is an artifact of the integration domain rather than a physical effect.

Watch

Extended reading notes

Core claim

On the paper's own terms, the discovery is that magnetized graphene's thermodynamic response is not fully described by the canonical ensemble once local inverse-temperature fluctuations are present. Starting from the Landau spectrum \(E_n=\pm\hbar v_F\sqrt{2nD}\) and the deformed Boltzmann factor \($B_E^{{(q)}}$(\$\beta$)=$e^{{-\beta E}}$\left(1+\frac{q}{2}\$beta^{2}$ $E^{2}$\right)\), the paper computes closed forms for the partition function, mean energy, specific heat, entropy, and free energy. The key numerical finding is that the superstatistical entropy exceeds the canonical entropy at the same \(\$\beta$\), while the superstatistical mean energy falls as \(2/\$\beta$\) and the specific heat is constant at \(2k\). The paper reads this as evidence that fluctuations introduce additional disorder and that superstatistics smooths out the extreme thermal effects seen in the canonical framework.

Load-bearing premise

The comparison that yields the entropy excess uses a finite cutoff \(\$\lambda$\) for the canonical partition function but an infinite integration domain for the superstatistical partition function; if that asymmetry is removed, the conclusion that fluctuations add disorder may not follow.

Editorial extensions

If this is right

  • The spectrum \(E_n=\pm\hbar v_F\sqrt{2nD}\) gives relativistic Landau quantization in graphene, with the \(n=0\) level shared by electrons and holes and no dependence on the wave number \(k_1\).
  • At a fixed \(\beta\), the superstatistical entropy \(S_s=k\left(2+\ln\frac{2+6q}{a^2\beta^2}\right)\) lies above the canonical entropy, so local temperature fluctuations add measurable disorder.
  • Superstatistical mean energy \(U_s=2/\beta\) and heat capacity \(C_s=2k\) are independent of the cutoff and of \(q\), whereas the canonical counterparts depend on the finite cutoff \(\lambda\), making the superstatistical description structurally simpler.
  • In the limit \(q\to 0\), the deformed Boltzmann factor reduces to the ordinary Boltzmann factor, so the canonical results are recovered as a special case.
  • The smoother superstatistical free energy implies that thermal fluctuations suppress the sharp low-temperature features of the canonical free energy.

Reading between the lines

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

  • A direct check is to compare the two frameworks with identical integration domains, for example both truncated at the same cutoff \(\lambda\); this would isolate whether the entropy excess comes from physics or from the choice of domain.
  • The q-deformed Boltzmann factor is one superstatistical model; other inverse-temperature distributions in Eq. (49) would modify the deformed factor, so the size of the entropy excess is model-dependent even if the qualitative direction of the effect is not.
  • Because the spectrum is the standard relativistic Landau spectrum, the same formalism should apply to other two-dimensional Dirac materials by replacing \(v_F\), an extension the paper does not pursue.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 6 minor

Summary. The manuscript solves the Dirac-Weyl equation for graphene in a constant magnetic field, obtaining the standard Landau-level spectrum En = ±ℏvF √(2nD) with D = eB0/(cℏ). From this spectrum it constructs a canonical partition function with a finite upper cutoff λ, computes the mean energy, heat capacity, entropy, and free energy, and then repeats the calculation in a superstatistical framework based on a q-deformed Boltzmann factor. The central physical claim is that superstatistics introduces additional disorder, evidenced by a larger entropy than in the canonical treatment. The spectrum derivation is standard and the canonical thermodynamic formulas are internally consistent, but the superstatistical partition function in Eq. (52) is neither derived nor correct, and the comparison underlying the entropy-excess claim is inconsistent in its integration domains and parameter choices.

Significance. If the superstatistical analysis were correct, the paper would offer a useful illustration of how local thermal fluctuations modify the thermodynamics of Dirac fermions in graphene. The strengths of the paper are the explicit Nikiforov-Uvarov solution for the Landau levels and the closed-form canonical partition function with a finite cutoff, which are straightforward and reproducible. However, the claimed new result — that superstatistics yields additional disorder — is unsupported: it follows from an incorrect superstatistical partition function and from an apples-to-oranges comparison between a finite-cutoff canonical integral and an infinite-domain superstatistical integral. The paper therefore does not provide a reliable new physical conclusion beyond the standard Landau-level spectrum.

major comments (3)
  1. [Section 4.2, Eqs. (51) and (52)] The superstatistical partition function Z_s in Eq. (52) is asserted rather than derived. Substituting the spectrum (37) and the deformed Boltzmann factor (50) into the defining integral (51) and evaluating directly gives Z_s = (1+3q)/(β² b), where b = D ℏ² v_F² as defined in Eq. (45); this result has no dependence on the cutoff parameter a and does not reduce to the canonical Z(β) of Eq. (45) when q=0. The printed expression (2+6q)/(a²β²) is therefore not obtainable from (51). Since Eqs. (53)–(57) are all built on Eq. (52), the entire superstatistical thermodynamic analysis is unsupported by the paper's own definitions.
  2. [Section 5, Fig. 5] The 'additional disorder' conclusion compares the superstatistical entropy S_s (Eq. (56)) with the canonical entropy S (Eq. (48)) under unequal conditions: the canonical partition function in Eqs. (43)–(45) uses a finite upper cutoff λ (encoded in a), whereas Z_s in (51) integrates n to infinity; furthermore, the plotted curves use different parameter sets (a = 1, b = 2 for the canonical case versus a = 1, q = 0.5 for the superstatistical case, with b absent from Eq. (52)). The entropy excess in Fig. 5 is therefore an artifact of the inconsistent domains and parameter choices, not a physical effect of local temperature fluctuations.
  3. [Section 4.2, Eq. (50)] The q-deformed Boltzmann factor B_E^(q)(β) = e^{-βE}(1 + q β² E²/2) is introduced without specifying the fluctuation distribution f(β',β) in Eq. (49) that would produce it; no derivation is given, and q is treated as an adjustable parameter (taken as 0.5 in Section 5). Because any positive q increases Z_s and hence S_s, the finding that superstatistics 'introduces additional disorder' is built in by the choice of q rather than being a falsifiable prediction of the model.
minor comments (6)
  1. [Figures 5 and 6 captions] The captions of Figures 5 and 6 refer to 'equation (49)' for the entropy and free energy, but the entropy is given in Eq. (48) and the free energy is given immediately after Eq. (48) without a number.
  2. [Section 3, Eq. (21)] Equation (21) writes E_n = ±ℏv_F √(E−); the argument of the square root should be ϵ⁻_n (or a clearly subscripted E⁻), since E− is used as the eigenvalue of the auxiliary Schrödinger problem in Eq. (26).
  3. [Section 4.1, Eq. (47)] Equation (47) contains an unbalanced parenthesis in the printed expression for C(β); the bracket structure should be corrected for readability.
  4. [References] References [29] and [30] both cite the same work by Beck (2004); one of them should be removed or replaced with the intended distinct source.
  5. [Figure 1 caption] The left panel of Figure 1 labels the horizontal axis as 'Hn', which should read 'n'; the caption phrase 'Hn = 0 to 5' is also inconsistent with the text.
  6. [Section 4.2, Eq. (52)] The parameter b is defined only in Eq. (45), but Z_s in Eq. (52) involves only a and β; since direct integration of (51) yields a b-dependent result, the omission of b in (52) obscures the inconsistency.

Circularity Check

2 steps flagged · score 6.0 of 10

Superstatistical 'additional disorder' claim is a restatement of the q-deformed Boltzmann factor ansatz; Eq. (52) is asserted, not derived.

  1. self definitional [Section 4.2, Eq. (50); Section 5, Fig. 5 discussion]
    "When f (β′, β) is modeled by a Dirac delta function δ(β −β′), the integral simplifies, and a deformation parameter q can be introduced to yield the generalized Boltzmann factor: B(q) E (β) = e−βE (1 + q 2 β2E2). Here, q ∈ [0, 1] is a deformation parameter that quantifies the departure from classical Boltzmann-Gibbs statistics. ... The higher value of Ss(β) compared to S(β) indicates that fluctuations introduce additional disorder."

    The entropy increase S_s > S in Fig. 5 is a direct algebraic consequence of the positive qβ²E² term inserted by hand in Eq. (50). q is not derived from graphene physics or from the delta-function average; it is set to q=0.5 in the figures. Thus the paper's central claim that fluctuations 'introduce additional disorder' is not a new result deduced from the spectrum, but a restatement of the q-deformation already assumed. A true delta-function f(β′,β) in Eq. (49) would give q=0, so the 'fluctuations' are an input, not an output.

  2. other [Section 4.2, Eq. (52)]
    "Using equations (50), (51) and Mathematica, the superstatistics partition function equation is given as : Zs(a, b, β, q) = 2 + 6q a2β2 . (52)"

    This step is asserted rather than derived: evaluating the defining integral (51) with E = ℏvF√(2Dn) and the deformed Boltzmann factor (50) yields a b-dependent result, not the printed a-only expression (2+6q)/(a²β²). The displayed Z_s also fails to reduce to the canonical partition function when q→0. Since the superstatistical entropy (56) and the 'additional disorder' comparison are built on Eq. (52), the central conclusion is a property of this un-derived formula, not a consequence of the model's own definition. This is not circularity in the fit sense, but it means the derivation chain breaks exactly where the claimed result enters.

full rationale

The paper independently derives the Dirac-Weyl Landau-level spectrum (Eq. 37) and canonical thermodynamic quantities from that spectrum, so those parts are not circular. The circularity resides in the superstatistical part. The q-deformed Boltzmann factor (Eq. 50) is introduced as an ansatz, with q chosen as 0.5 rather than derived. The abstract and Section 5 then present 'fluctuations introduce additional disorder' as a discovery, but that conclusion is essentially the same positive q term rewritten as an entropy excess; the superstatistical framework already contains the fluctuation by construction. In addition, Eq. (52) is not obtained from the defining integral (51), so the entropy curve that drives the headline claim is built on an asserted closed form rather than on the model's own equations. The inconsistent integration domains (finite λ for the canonical Z, infinite for Z_s) further make the entropy comparison an artifact of the setup. Overall, the central superstatistics claim reduces to the q-deformation ansatz and to an un-derived partition function, while the rest of the paper is a standard textbook calculation.

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

The paper's central thermodynamic results rest on several assumptions: the Dirac-Weyl model for graphene, a specific deformed Boltzmann factor, and different integration ranges for the canonical versus superstatistical partition functions. The q parameter is chosen ad hoc, and the cutoff λ is arbitrary.

free parameters (3)
  • λ (upper cutoff in canonical partition function) = not given; figures use a=1, b=2 which imply λ = a²/(2b) = 0.25
    The canonical partition function truncates the integral over Landau levels at an arbitrary highest quantum number λ (Eq. 43). This cutoff is not derived from physics and affects all canonical thermodynamic quantities.
  • q (superstatistical deformation parameter) = 0.5 in figures
    The q parameter in the deformed Boltzmann factor (Eq. 50) is chosen as 0.5 for the plots without any derivation from a physical temperature-fluctuation distribution.
  • a and b (combinations of physical constants and λ) = a=1, b=2 in figures
    The figures use a=1 and b=2, but these are not connected to actual values of the Fermi velocity, magnetic field, or temperature. They are effectively arbitrary scale choices.
assumptions (4)
  • domain assumption Graphene electrons obey the Dirac-Weyl equation with Fermi velocity vF = c/300 and minimal coupling in Gaussian units (Eqs. 14-15).
    This is the standard model for low-energy electrons in graphene, cited to Refs. [2,3,11,12].
  • ad hoc to paper The deformed Boltzmann factor in Eq. (50), B_E(q)(β) = e^{-βE}(1 + q β²E²/2), is the correct superstatistical generalization.
    Eq. (50) is presented without derivation or citation of the specific temperature-fluctuation distribution that produces it. The q parameter is free and unphysical as used here.
  • domain assumption The sum over Landau levels can be approximated by an integral with an upper cutoff λ in the canonical calculation (Eq. 43).
    This classical-limit approximation is common in statistical mechanics but introduces a free cutoff that is not physically motivated for graphene.
  • ad hoc to paper The superstatistical partition function integrates over n from 0 to ∞ (Eq. 51), unlike the canonical finite-cutoff calculation.
    The different integration domains in the two frameworks are not justified and are the source of the apparent entropy excess.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Thermodynamic Properties and Superstatistics of Graphene under a Constant Magnetic Field." pith.science (2026). https://pith.science/paper/3JWBWMBP

@misc{pith2026250600709,
  author       = {Pith},
  title        = {Pith review of: Thermodynamic Properties and Superstatistics of Graphene under a Constant Magnetic Field},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3JWBWMBP}},
  note         = {Machine review of arXiv:2506.00709}
}
read the original abstract

In this paper, we present the solutions of the Dirac-Weyl equation for graphene under a constant magnetic field. The resulting spectrum is used to determine the partition function, a key quantity in the study of thermodynamic properties. From this function, we analyze the mean energy, specific heat, entropy, and free energy in two different frameworks: the canonical ensemble and the superstatistical approach. The study confirms the relativistic nature of electron transport in graphene under a magnetic field. It also reveals that fluctuations introduce additional disorder in the system. The obtained results are in good agreement with those already reported in the literature.

Figures

Figures reproduced from arXiv: 2506.00709 by the authors.

Figure 1
Figure 1. Variation of the energy spectrum from equation (37) on the left and that of [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. Variation of the partition function in the classical framework on the left [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Variation of the mean energy in the classical framework on the left (equation [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: Variation of the specific heat in the classical regime (equation (48)) on the [PITH_FULL_IMAGE:figures/full_fig_p010_4.png]
Figure 5
Figure 5. Figure 5: Behavior of the entropy in the classical regime (equation (49)) on the left and [PITH_FULL_IMAGE:figures/full_fig_p011_5.png]
Figure 6
Figure 6. Figure 6: Variation of the free energy in the canonical framework (equation (49)) on [PITH_FULL_IMAGE:figures/full_fig_p011_6.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

36 extracted references · 36 canonical work pages

  1. [1]

    Novoselov, A.K

    K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, S.V. Dobonos and I.V. Grigorieva, Science 306 (2004) 666–669

    work page 2004

  2. [2]

    Castro Neto, F

    A.H. Castro Neto, F. Guinea, N.M.R. Peres, K.S. Novoselov, A.K. Giem, Rev. Mod. Phys. 81 (2009) 109–162

    work page 2009

  3. [3]

    Kotov, B

    V.N. Kotov, B. Uchoa, V.M. Pereira, F. Guinea, A.H. Castro Neto, Rev. Mod. Phys. 84 (2012) 1067–1125

    work page 2012

  4. [4]

    Peres, A.H

    N.M.R. Peres, A.H. Castro Neto, F. Guinea, Phys. Rev. B 73 (2006) 241403(R)

    work page 2006

  5. [5]

    De Martino, L.D

    A. De Martino, L.D. Anna, R. Egger, Phys. Rev. Lett. 98 (2007) 066802

    work page 2007

  6. [6]

    Myoung, G

    N. Myoung, G. Ihm, Physica E 42 (2009) 70

    work page 2009

  7. [7]

    Ghosh, J

    T.K. Ghosh, J. Phys. Condens. Matter 21 (2009) 045505

    work page 2009

  8. [8]

    Bliokh, V

    Y.P. Bliokh, V. Freilikher, F. Nori, Phys. Rev. B 81 (2010) 075410

    work page 2010

Show all 36 references
  1. [9]

    Milpas, M

    E. Milpas, M. Torres, G. Murguia, J. Phys. Condens. Matter 23 (2011) 245304

    2011

  2. [10]

    De Silva Leite, C

    L.G. De Silva Leite, C. Filgueiras, D. Cogollo, Edilberto O. Silva, Phys. Lett. A 379 (2015) 907–911

    2015

  3. [11]

    S. Kuru, J. Negro, L. M. Nieto, J. Phys.: Condens. Matter, 21 (2009) 455305

    2009

  4. [12]

    David, C

    J. David, C. Fern´ andez, D. Juan Garc ´ ıa-Mu˜ noz, Eur. Phys. J. Plus 137 (2022) 1013

    2022

  5. [13]

    Eshghi, H

    M. Eshghi, H. Mehraban, C. R., Physique 18 (2017) 47–56

    2017

  6. [14]

    M. R. Setare, O. Hatami, CHIN.PHYS.LETT. 11 (2008) 3848

    2008

  7. [15]

    Hosseini, H

    M. Hosseini, H. Hassanabadi, S. Hassanabadi, Eur. Phys. J. Plus 134 (2019) 6

    2019

  8. [16]

    D. K. Demir, S. Kuru, J. Negro, Physica. E. 118 (2020) 113926

    2020

  9. [17]

    Sun W., Liu Y., Li M., Cheng Q., Zhao L., Energy., 2023, 269, 127001

    2023

  10. [18]

    Bai B., Zhou R., Cai G., Hu, W., Yang G., 2021, 137, 104272

    2021

  11. [19]

    Kuang W., Wang H., Li X., Zhang J., Zhou Q., Zhao Y., Acta Materialia., 2018, 159, 16

    2018

  12. [20]

    F., Uvarov V

    Nikiforov A. F., Uvarov V. B., Special Functions of Mathematical Physics, Birkhauser Basel, 1988, 10.1007/978-1-4757-1595-8

    1988 doi

  13. [21]

    Castro Neto A. H,. Guinea F., Peres N. M., Novoselov K. S. Geim A. K., Rev. Mod. Phys. 81 (2009) 109

    2009

  14. [22]

    E˘ grifes H., Demirhan D., B¨ uy¨ ukkili¸ c F., Physica Scripta., 1999, 60, 195

    1999

  15. [23]

    B¨ uy¨ ukkili¸ c F., E˘ grifes H., Demirhan D., Theor. Chem. Acc., 1997, 98, 192

    1997

  16. [24]

    Infeld L., Hull T. E., Rev. Mod. Phys. 23 (1951) 21

    1951

  17. [25]

    Demircioglu, S

    B. Demircioglu, S. Kuru, ¨Onder M., Ver¸ cin A., J. Math. Phys. 43 (2002) 2133

    2002

  18. [26]

    Kasap, B

    E. Kasap, B. G¨ on¨ ul, M. S ¸im¸ sek, Chem. Phys. Lett., 172 (1990) 499

    1990

  19. [27]

    N., Okon I

    Isonguyo C. N., Okon I. B., Antia A. D., Oyewumi K. J., Omugbe E., Onate C. A., Joshua R. U., Udoh M. E., Ituen Eno. E., Aruajo J. P., Frontiers in Physics., 2022, 10, 962717

    2022

  20. [28]

    M., Daniel S

    Assimiou Y. M., Daniel S. T., Issoufou G., Anselme D. F., and Avossevou G. Y. H., Condensed Matter Physics., 2024, Vol. 27, No. 4, 43301

    2024

  21. [29]

    Beck, Superstatistics: theory and application, Contin

    C. Beck, Superstatistics: theory and application, Contin. Mech. Thermod. 16, 293, (2004)

    2004

  22. [30]

    Superstatistics: theory and applications,

    C. Beck, “Superstatistics: theory and applications,” Continuum Mechanics and Thermodynamics, vol. 16, no. 3, pp. 293–304, 2004

    2004

  23. [31]

    C. O. Edet, P. O. Amadi, U. S. Okorie, A. Tas, A. N. Ikot, and G. Rampho, Revis. Mexic. de fisca. 66, 824, (2020)

    2020

  24. [32]

    Superstatistics,

    C. Beck and E. G. Cohen, “Superstatistics,” Physica A, vol. 322, pp. 267–275, 2003

    2003

  25. [33]

    Ikot ,U.S

    A.N. Ikot ,U.S. Okorie, G. Osobonye,P.O. Amadi ,C.O. Edet, M.J. Sithole, G.J. Rampho, R. Sever, Heliyon 6 (2020) e03738

    2020

  26. [34]

    C. O. Edet, P. O. Amadi, U. S. Okorie, A. Tas, A. N. Ikot, and G. Rampho, Revista Mexicana defisca, vol. 66, pp. 824–839, 2020

    2020

  27. [35]

    Sargolzaeipor, H

    S. Sargolzaeipor, H. Hassanabadi,W.S. Chung, Eur. Phys. J. Plus 133 (2018)5

    2018

  28. [36]

    T., Assimiou Y

    Daniel S. T., Assimiou Y. M., and Avossevou G. Y. H., arXiv:submit/6464474 [quant-ph] 21 May 2025

    2025

Pith tools

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