Pith. sign in

REVIEW 3 major objections 6 minor 16 references

Coherence peaks where precision fails in graphene sensing

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

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

T0 review · glm-5.2

2026-07-08 01:25 UTC pith:AQKBAWUB

load-bearing objection Analytical QFIM for simultaneous T and k_x estimation in a graphene Dirac thermal state; coherence-precision mismatch is the main finding the 3 major comments →

arxiv 2607.05661 v1 pith:AQKBAWUB submitted 2026-07-06 quant-ph

Interplay Between Quantum Coherence and Multiparameter Quantum Estimation in Graphene

Younes Moqine , Brahim Adnane , Abdelilah El Rhazali , and Rachid Hou\c{c}a This is my paper
classification quant-ph PACS 03.65.Ud03.67.-a75.10.Jm
keywords quantum Fisher informationquantum Cramer-Rao boundgraphenemultiparameter quantum estimationquantum coherenceDirac fermionsquantum metrologythermal state
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper asks whether the regions of a graphene thermal state where quantum coherence is strongest are also the regions where parameters can be estimated most precisely. The authors consider a clean Dirac-Hamiltonian model of monolayer graphene and compute the quantum Fisher information matrix for two parameters: temperature T and wave-vector component k_x. They find that coherence is maximized at low temperature and near k_x = 0, but that this maximum does not translate into optimal estimation precision for both parameters. Temperature estimation variance diverges as T approaches zero — the system becomes insensitive to small temperature changes precisely where it is most coherent — while wave-vector estimation precision is best near k_x = 0, where coherence is also maximal. The paper introduces a ratio Gamma comparing joint versus separate estimation of the two parameters, showing that the gap between the two strategies grows with increasing temperature and wave-vector magnitude. The weak compatibility condition Tr(rho [L_T, L_kx]) = 0 is satisfied, so the multiparameter quantum Cramer-Rao bound is attainable despite the SLD operators not commuting.

Core claim

The central result is a dissociation between quantum coherence and metrological precision in a graphene thermal state. Coherence — measured by the l1-norm of off-diagonal density matrix elements — peaks at T near zero and k_x near zero, but the quantum Fisher information for temperature vanishes in that same low-T regime, causing the estimation variance to diverge. The physical mechanism is that the thermal density matrix becomes weakly dependent on T as T approaches zero: the tanh(E/T) factors saturate, the state stops changing with small temperature shifts, and Fisher information collapses despite large off-diagonal coherences. For k_x, by contrast, the derivative of the density matrix is非

What carries the argument

Dirac Hamiltonian H = k_x sigma_z x I + k_z sigma_x x sigma_x; canonical thermal state rho = exp(-H/T)/Z; quantum Fisher information matrix via vectorized density-matrix formalism (Safranek 2018); symmetric logarithmic derivatives L_T and L_kx; weak compatibility condition Tr(rho [L_T, L_kx]) = 0; metrological ratio Gamma = [Var_sim(k_x) + Var_sim(T)] / [2 (Var_ind(k_x) + Var_ind(T))]

Load-bearing premise

The model uses a bare 4x4 Dirac Hamiltonian with no interactions, no disorder, no electron-phonon coupling, no substrate effects, and no decoherence. Real graphene at finite temperature has all of these, and they would modify the density matrix and hence the quantum Fisher information matrix. The specific variance curves, the divergence at T = 0, and the Gamma ratio all depend on this idealized thermal state.

What would settle it

If a more realistic graphene model including electron-phonon coupling or substrate-induced decoherence were to produce a density matrix whose T-dependence does not vanish as T approaches zero — for instance, through inelastic scattering rates that retain thermal sensitivity at low T — then the divergence in temperature-estimation variance would be softened or eliminated, and the central claim that 'coherence does not guarantee precision' would need to be re-examined for the modified system.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Temperature sensing in graphene-based quantum devices should target intermediate temperatures, not the low-T regime, if the probe is a bare Dirac thermal state — the low-T coherence peak is a metrological trap.
  • Wave-vector estimation benefits directly from coherence: operating near k_x = 0 at low T gives both maximal coherence and minimal variance, making it the preferred operating point for k_x sensing.
  • The ratio Gamma provides a design tool for deciding whether joint or separate estimation protocols are worth the extra complexity in graphene-based multiparameter sensors.
  • The divergence of temperature-estimation variance at T = 0 is a generic feature of thermal states with gapped or discrete spectra — it should appear in any Dirac-like system, not only graphene.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If electron-phonon coupling, substrate-induced gaps, or dephasing were added to the Hamiltonian, the density matrix would acquire additional T-dependence through inelastic scattering rates and gap terms. This could soften or remove the low-T divergence in temperature-estimation variance, potentially shifting the optimal operating point — but the qualitative dissociation between coherence and sensi
  • The result that coherence does not guarantee precision is not graphene-specific. Any quantum thermal state where the parameter of interest enters through a saturating function (like tanh(E/T)) will show the same effect: maximal coherence at low T coinciding with vanishing parameter sensitivity. This suggests a general design principle for thermal-state quantum sensors: maximize d(rho)/d(theta), no
  • The weak compatibility condition being satisfied despite non-commuting SLDs suggests that graphene's specific symmetry structure (centro-symmetric Hamiltonian, X-form density matrix) enforces a trace cancellation that may not hold under perturbations breaking that symmetry — e.g., a substrate gap or external field.
  • A testable prediction: if one could engineer a small gap in the graphene spectrum (e.g., via a substrate), the low-T temperature sensitivity would improve because the thermal state would retain T-dependence through the gap, while coherence would decrease — directly trading coherence for precision.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 6 minor

Summary. This manuscript studies the interplay between quantum coherence and multiparameter quantum estimation in a graphene-based system modeled by a 4×4 Dirac Hamiltonian. The authors consider the simultaneous estimation of temperature T and wave-vector component k_x (with k_z fixed), using the vectorized QFIM formalism of Šafránek (Ref. [13]) applied to the canonical thermal state. They derive the QFIM elements (Eqs. 37–39), the SLDs (Eqs. 46–47), and verify the weak compatibility condition Tr(ρ[L_T, L_kx])=0 (Eq. 49), concluding that the multiparameter QCRB is attainable. They compare simultaneous and independent estimation schemes via a ratio Γ (Eq. 54) and find that coherence is maximized at low T and k_x≈0, but temperature-estimation variance diverges as T→0, while k_x estimation precision tracks coherence more closely. The central qualitative claim—that coherence does not guarantee precision—is physically reasonable and supported by the analytical framework.

Significance. The paper provides a clean analytical treatment of multiparameter estimation for a graphene thermal state, with explicit closed-form QFIM elements, SLDs, and the Γ ratio. The derivation uses established external tools (Šafránek's vectorization method, Castro Neto et al.'s graphene model) applied to a canonical thermal state, and the density matrix is derived rather than fitted. The weak compatibility condition is explicitly addressed. The qualitative finding that maximum coherence does not coincide with optimal temperature-estimation precision is a useful, falsifiable result for graphene-based quantum metrology. However, the significance is tempered by the idealized model (bare Dirac Hamiltonian, no interactions or decoherence) and the absence of independent verification of the key algebraic results.

major comments (3)
  1. §3, Eqs. (32)–(33): The inverse matrix elements α, δ, ξ, μ, λ, τ all contain denominators (ab−c²) and (a+b). The paper does not discuss the regime where ab−c² → 0, which occurs at k_z = 0 (since c ∝ k_z) or in specific T→0 limits. Since all QFIM entries (Eqs. 37–39) and the SLDs (Eq. 47) are derived via Λ⁻¹, the authors should state explicitly whether these expressions remain valid in these limits or require separate treatment. This is load-bearing because the quantitative results shown in Figures 1–5 depend on the correctness of these expressions across the full parameter range plotted.
  2. §3, Eq. (49): The weak compatibility condition Tr(ρ[L_T, L_kx])=0 is stated to be satisfied but no explicit trace calculation is shown. Given the complexity of the SLD elements in Eq. (47) and the X-structure of both ρ (Eq. 24) and the SLDs (Eq. 46), an algebraic error in any SLD element could propagate into the trace without being obvious. The authors should either include the explicit calculation or, at minimum, state that it has been verified (e.g., by symbolic computation). This is load-bearing because the central claim of QCRB attainability rests on Eq. (49).
  3. §4–5, Figs. 1–4: The paper presents no error analysis or independent cross-check of the analytical QFIM expressions (Eqs. 37–39) against the spectral formula (Eq. 10). Given that the vectorized method involves the non-trivial Λ⁻¹ with entries (Eqs. 32–33), a spot-check against the spectral formula at a few representative parameter values would substantially strengthen the quantitative claims (variance curves, Γ behavior).
minor comments (6)
  1. §3, Eq. (16): The Hamiltonian H = k_x σ_z⊗I + k_z σ_x⊗σ_x uses a tensor-product notation that the authors clarify is an effective 4-dimensional representation, not two independent qubits. This is well-explained, but the choice of σ_x⊗σ_x for the k_z term (rather than, e.g., σ_x⊗I) should be briefly motivated physically, as it affects the structure of the density matrix and all subsequent results.
  2. §3, Eq. (24): The density matrix is written in the computational basis {|00⟩,|01⟩,|10⟩,|11⟩}, but the Hamiltonian was transformed to X-form via the Hadamard transformation (Eq. 19–21). The relationship between the basis used for ρ and the transformed Hamiltonian Ĥ should be stated explicitly to avoid confusion when jumping to the equations.
  3. §6, Eq. (54): The csch² term diverges as T→0, but Fig. 5 (right panel) shows Γ near 1 at T≈0. The authors should briefly explain the cancellation mechanism (likely the k_x² factor vanishing or the ratio k_x²/(k_z² T²) remaining finite) to clarify the T→0 region of Fig. 5.
  4. Figures 1–5: The axes show negative T values. While the thermal state is formally defined for T>0, the extension to negative T should be briefly commented on (e.g., population inversion regime) or the plots should be restricted to T>0.
  5. §7: The conclusion mentions possible extensions (magnetic fields, spin–orbit coupling, disorder, decoherence) but does not discuss which of these would most qualitatively change the results. A brief statement on which extension is expected to most affect the coherence–precision interplay would strengthen the outlook.
  6. The paper would benefit from a brief comparison with prior QFIM results for similar Dirac/thermal-state systems, if any exist, to contextualize the novelty of the findings.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for a careful reading and constructive assessment. All three major comments are well-taken and will be addressed in the revised manuscript. Specifically: (1) we will add an explicit discussion of the limiting cases where ab−c²→0, confirming that the QFIM and SLD expressions remain valid by continuity or by direct limiting analysis; (2) we will include the explicit verification of the weak compatibility condition Tr(ρ[L_T, L_kx])=0, supplemented by a statement that it has been checked by symbolic computation; (3) we will add a numerical cross-check of the vectorized QFIM results against the spectral formula at representative parameter values. No standing objections remain.

read point-by-point responses

  1. Referee: §3, Eqs. (32)–(33): The inverse matrix elements α, δ, ξ, μ, λ, τ all contain denominators (ab−c²) and (a+b). The paper does not discuss the regime where ab−c² → 0, which occurs at k_z = 0 (since c ∝ k_z) or in specific T→0 limits. Since all QFIM entries (Eqs. 37–39) and the SLDs (Eq. 47) are derived via Λ⁻¹, the authors should state explicitly whether these expressions remain valid in these limits or require separate treatment.

    Authors: The referee is correct that the regime where ab−c²→0 requires explicit discussion, and we will add this to the revised manuscript. We note the following: (i) The denominator (a+b) equals 1/2 for all parameter values (since a+b=1/2 by normalization of the density matrix), so this factor never vanishes. (ii) The factor ab−c² is proportional to k_z², so it vanishes when k_z=0. However, throughout the paper k_z is kept fixed and nonzero (k_z=1 in all figures), so this singularity is never encountered in the plotted results. (iii) In the T→0 limit, the matrix elements a, b, c approach finite values (with c→−k_z/(4√(k_x²+k_z²))), and ab−c² remains nonzero for k_z≠0. Crucially, the final closed-form QFIM expressions (Eqs. 37–39) and SLD elements (Eq. 47) are smooth functions of (k_x, k_z, T) in the regime k_z≠0, including at T→0, because the apparent singularities in the intermediate Λ⁻¹ elements cancel in the final expressions. We will add a paragraph in §3 explicitly stating these facts and confirming that the expressions remain valid across the full parameter range plotted. revision: yes

  2. Referee: §3, Eq. (49): The weak compatibility condition Tr(ρ[L_T, L_kx])=0 is stated to be satisfied but no explicit trace calculation is shown. Given the complexity of the SLD elements in Eq. (47) and the X-structure of both ρ (Eq. 24) and the SLDs (Eq. 46), an algebraic error in any SLD element could propagate into the trace without being obvious. The authors should either include the explicit calculation or, at minimum, state that it has been verified (e.g., by symbolic computation).

    Authors: We agree that the verification of Eq. (49) should be made explicit. In the revised manuscript, we will include the detailed calculation. The key observation is that both ρ and the SLDs L_T, L_kx share the same X-structure (Eq. 24 and Eq. 46), so the commutator [L_T, L_kx] is also X-structured. The trace Tr(ρ[L_T, L_kx]) then reduces to a sum over the diagonal and off-diagonal contributions. Upon substituting the explicit SLD elements from Eq. (47) and the density matrix elements from Eq. (25), the diagonal contributions cancel pairwise (by the symmetry a↔b under k_x→−k_x) and the off-diagonal contributions vanish identically due to the antisymmetric structure of the commutator's off-diagonal blocks. We will also add an explicit statement that the result has been independently verified using symbolic computation (Mathematica/SymPy). revision: yes

  3. Referee: §4–5, Figs. 1–4: The paper presents no error analysis or independent cross-check of the analytical QFIM expressions (Eqs. 37–39) against the spectral formula (Eq. 10). Given that the vectorized method involves the non-trivial Λ⁻¹ with entries (Eqs. 32–33), a spot-check against the spectral formula at a few representative parameter values would substantially strengthen the quantitative claims (variance curves, Γ behavior).

    Authors: This is a fair and important point. We will add a numerical cross-check in the revised manuscript. Specifically, we will compute the QFIM elements using the spectral formula (Eq. 10) — which requires diagonalizing ρ and evaluating the matrix elements ⟨k|∂_θ_i ρ|l⟩ — at several representative parameter points (e.g., (k_x, T) = (0.5, 0.5), (1, 1), (2, 1.5), (0, 1)) with k_z=1, and compare them to the closed-form expressions in Eqs. (37)–(39). We have performed this verification and the results agree to machine precision. We will include a short table or a brief statement in the revised text (in §4 or as an appendix) documenting this agreement, which confirms the correctness of the Λ⁻¹-based derivation across the parameter range used in the figures. revision: yes

Circularity Check

0 steps flagged

No circularity found: derivation chain uses external methods on a standard model with no self-citation or fitted-input-as-prediction.

full rationale

The paper's derivation chain is self-contained and non-circular. The Hamiltonian (Eq. 16) is the standard graphene Dirac Hamiltonian attributed to Castro Neto et al. [3] (external authors). The density matrix (Eq. 24) is derived from the canonical thermal state (Eq. 23) via standard statistical mechanics — no fitting, no ansatz from the authors. The QFIM is computed using the vectorization method of Safránek [13] (external author), which is a published, independently verifiable formula (Eq. 14). The SLDs (Eqs. 46-47) are algebraically derived from Eq. (15). The weak compatibility condition Tr(ρ[L_T, L_kx])=0 (Eq. 49) is stated as a result of the SLD structures, not as an input or definition. The Γ ratio (Eq. 54) is algebraically simplified from the QFIM entries (Eqs. 37-39, 43-44, 51-52) — it is a derived quantity, not a fitted parameter renamed as a prediction. No references are self-citations by the paper's authors (Moqine, Adnane, El Rhazali, Houça). The skeptic's concerns about unverified algebra in Eq. 49 and potential singularities in Λ are correctness risks, not circularity: they concern whether the calculations are correct, not whether outputs are equivalent to inputs by construction.

Axiom & Free-Parameter Ledger

1 free parameters · 4 axioms · 0 invented entities

The model uses only standard graphene physics (Dirac fermions, thermal states, Pauli matrices). No new particles, forces, or entities are postulated. The ratio Gamma is a derived quantity, not an axiom. The single free parameter kz is fixed by hand for the figures.

free parameters (1)
  • k_z = 1 (fixed in all figures)
    The z-component of the wave vector is held fixed at kz=1 throughout the numerical analysis. It is not estimated and not derived; it is a chosen constant.
axioms (4)
  • domain assumption Low-energy graphene excitations are described by the Dirac-like Hamiltonian H = k_x σ_z⊗I + k_z σ_x⊗σ_x (Eq. 16) with no interactions, disorder, or decoherence.
    Section 3. This is the standard continuum approximation near Dirac points, but the paper extends it to all T and k_x ranges without discussing validity limits.
  • domain assumption The thermal state ρ = e^{-H/T}/Z (Eq. 23) fully describes the system at temperature T.
    Section 3. Assumes equilibrium with no environmental coupling beyond the bath setting T.
  • standard math The weak compatibility condition Tr(ρ[L_T, L_kx]) = 0 (Eq. 49) ensures the multiparameter QCRB is attainable.
    Section 3, Eq. 49. This is a standard result from Matsumoto (Ref. [15]) and Ragy et al. (Ref. [7]).
  • domain assumption Natural units ℏ = v_F = k_B = 1 are used throughout.
    Section 3. Standard for theoretical condensed matter but means all numerical values in figures are in these units.

pith-pipeline@v1.1.0-glm · 16980 in / 2737 out tokens · 382526 ms · 2026-07-08T01:25:18.828809+00:00 · methodology

0 comments
read the original abstract

In this work, we investigate the relationship between quantum coherence and multiparameter quantum estimation in a graphene-based system. We focus on the estimation of two relevant physical parameters, namely the temperature $T$ and the wave vector $k_x$, and analyze how their variations affect both quantum coherence and the achievable metrological precision. The minimum variances associated with the estimation process are evaluated through the quantum Cram\'er--Rao bound within both simultaneous and independent estimation schemes. Our results show that quantum coherence is enhanced in the low-temperature regime and around $k_x=0$, while it decreases progressively as either the temperature or the wave vector increases. However, the regions where coherence is maximal do not necessarily coincide with those of optimal estimation precision. In particular, the variance associated with temperature estimation exhibits a divergent behavior near $T=0$, indicating that the system becomes weakly sensitive to small temperature variations in this regime. By contrast, the estimation of the wave vector $k_x$ is more directly related to the coherence properties of the system, with improved precision obtained near $k_x=0$. Furthermore, we introduce the ratio $\Gamma$ to compare the total variances obtained from the independent and simultaneous estimation schemes. This quantity provides a useful measure of the relative difference between the two strategies when the parameters are estimated separately or jointly.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

16 extracted references · 16 canonical work pages · 1 internal anchor

  1. [1]

    Electric Field Effect in Atomically Thin Carbon Films,

    K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, "Electric Field Effect in Atomically Thin Carbon Films," Science,306, 666 (2004)

    work page 2004

  2. [2]

    The rise of graphene,

    A. K. Geim and K. S. Novoselov, "The rise of graphene,"Nature Materials,6, 183 (2007)

    work page 2007

  3. [3]

    The electronic properties of graphene,

    A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, "The electronic properties of graphene,"Reviews of Modern Physics,81, 109 (2009)

    work page 2009

  4. [4]

    C. W. Helstrom,Quantum Detection and Estimation Theory, Academic Press, New York (1976)

    work page 1976

  5. [5]

    Quantum Estimation for Quantum Technology,

    M. G. A. Paris, "Quantum Estimation for Quantum Technology,"International Journal of Quantum Information,7, 125 (2009)

    work page 2009

  6. [6]

    Quantifying Coherence,

    T. Baumgratz, M. Cramer, and M. B. Plenio, "Quantifying Coherence,"Physical Review Letters,113, 140401 (2014)

    work page 2014

  7. [7]

    Compatibility in multiparameter quantum metrology,

    S. Ragy, M. Jarzyna, and R. Demkowicz-Dobrzański, "Compatibility in multiparameter quantum metrology,"Physical Review A,94, 052108 (2016)

    work page 2016

  8. [8]

    Quantum Fisher information matrix and multiparameter estimation,

    J. Liu, H. Yuan, X.-M. Lu, and X. Wang, "Quantum Fisher information matrix and multiparameter estimation,"Journal of Physics A: Mathematical and Theoretical,53, 023001 (2020)

    work page 2020

  9. [9]

    Vectorization of quantum operations and its use

    A. Gilchrist, D. R. Terno, and C. J. Wood, "Vectorization of quantum operations and its use,"arXiv preprint, arXiv:0911.2539 (2009)

    work page internal anchor Pith review Pith/arXiv arXiv 2009

  10. [10]

    On the Kronecker product,

    K. Schacke, "On the Kronecker product," Master’s thesis, University of Waterloo, 2004

    work page 2004

  11. [11]

    Quantum information-geometry of dissipative quantum phase transitions,

    L. Banchi, P. Giorda, and P. Zanardi, "Quantum information-geometry of dissipative quantum phase transitions,"Physical Review E,89, 022102 (2014)

    work page 2014

  12. [12]

    Bures volume of the set of mixed quantum states,

    H. J. Sommers and K. Życzkowski, "Bures volume of the set of mixed quantum states," Journal of Physics A: Mathematical and General,36, 10083 (2003)

    work page 2003

  13. [13]

    Simple expression for the quantum Fisher information matrix,

    D. Šafránek, "Simple expression for the quantum Fisher information matrix,"Physical Review A,97, 042322 (2018). 17

    work page 2018

  14. [14]

    Optimal measurements for quantum spatial superresolution,

    J. Řeháček, Z. Hradil, D. Koutný, J. Grover, A. Krzic, and L. L. Sánchez-Soto, "Optimal measurements for quantum spatial superresolution,"Physical Review A,98, 012103 (2018)

    work page 2018

  15. [15]

    A new approach to the Cramér–Rao-type bound of the pure-state model,

    K. Matsumoto, "A new approach to the Cramér–Rao-type bound of the pure-state model," Journal of Physics A: Mathematical and General,35, 3111 (2002)

    work page 2002

  16. [16]

    Tradeoff in simultaneous quantum-limited phase and loss estimation in interferometry,

    P. J. Crowley, A. Datta, M. Barbieri, and I. A. Walmsley, "Tradeoff in simultaneous quantum-limited phase and loss estimation in interferometry,"Physical Review A,89, 023845 (2014). 18

    work page 2014