{"id":"afb16892-4186-42c1-94a1-061dd432bb7b","arxiv_id":"2411.17368","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A finite-length R-CPE transmission line is solved exactly, yielding the anomalous finite-length Warburg impedance s^{-alpha/2} coth(s^{alpha/2}) plus its step response and relaxation-time distribution.","lead":"This paper solves the time-fractional diffusion equation for a finite-length transmission line made of resistors and constant phase elements, and derives the impedance, step response, and relaxation-time distribution for porous electrodes under galvanostatic charging. The work gives a mathematical basis to an empirical impedance formula widely used to fit supercapacitor and battery data.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The DRT derivation in Eq. 52 contains an algebraic error: the substitution u=b^{-1/alpha}v gives a denominator b^2+2b v^alpha cos(alpha pi)+v^{2alpha}, not the b^{1/alpha+1}, b^{1/alpha}, b^{1/alpha-1} form used in Eqs. 52-53, so the claimed DRT and its Fig. 4(c) agreement are suspect.","rationale":"The paper's central claim as emphasized by the reader is Eq. 44, and my independent re-derivation confirms that Eq. 44 follows from Eq. 8 under the uniform-pore r=r0 assumption, so I do not object to the impedance formula itself. The reader's flagged weakest assumption (uniformity, constant CPE parameters, blocking end) is a legitimate modeling limitation, but it is not an internal error and does not threaten the mathematical derivation within the stated model. The more load-bearing concern I found is the DRT derivation: Eq. 52 is algebraically incorrect by a factor b^{1/alpha-1}, which propagates into Eq. 53 and invalidates the plotted g(tau) as the DRT of the derived impedance. Since the paper explicitly claims to have derived the DRT and uses the agreement between Eq. 53 and Eq. 57 as supporting evidence, this is a correctness risk in a claimed contribution, not merely an unclear normalization. My recommended verdict remains CONDITIONAL rather than REJECT because the impedance formula is correct and the DRT section is local and fixable; this does not change the reader's overall verdict, so I mark the verdict unchanged. I partially agree with the reader because the reader noted the DRT normalization was unclear, but our specific algebraic diagnosis is different from the reader's stated weakest assumption.","tokens_in":9,"tokens_out":27350,"duration_ms":304693,"concrete_test":"Recompute Eq. 52 step by step: substitute u=b^{-1/alpha} v into Eq. 51, multiply numerator and denominator by b^2, and cancel the common b^{1-1/alpha} factor; if the denominator is b^2+2b v^alpha cos(alpha pi)+v^{2alpha}, then Eq. 52 is wrong. Then rederive Eq. 53 from the corrected expression, recompute g(tau) for alpha=0.5, 0.75, 0.90, and compare with Eq. 57 and with the impedance integral Eq. 54 evaluated numerically. If the corrected g(tau) differs visibly from Fig. 4(c), the DRT section must be revised before the paper can be accepted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central impedance result, Eq. 44, follows correctly from Eq. 8 under the stated r=r0 uniformity assumption, so the main impedance claim is not the place I would attack. The load-bearing flaw is in the distribution-of-relaxation-times derivation. Starting from Eq. 51 and substituting x=b s^alpha, u=b^{-1/alpha} v, the integrand becomes v^alpha e^{-s v} divided by 1+2b^{-1}v^alpha cos(alpha pi)+b^{-2}v^{2alpha}. Multiplying numerator and denominator by b^2 and cancelling the common factor b^{1-1/alpha} on both sides gives exactly s^{alpha-1} E_{alpha,alpha}(-b s^alpha) = (sin(alpha pi)/pi) integral v^alpha e^{-s v}/(b^2+2b v^alpha cos(alpha pi)+v^{2alpha}) dv. Eq. 52 instead has denominator b^{1/alpha+1}+2b^{1/alpha}v^alpha cos(alpha pi)+b^{1/alpha-1}v^{2alpha}, which is b^{1/alpha-1} times the correct denominator; the extra b^{1/alpha-1} factor is spurious. Eq. 53 therefore does not give the DRT of Eq. 44/46, and the 'excellent agreement' with Eq. 57 in Fig. 4(c) cannot be a valid check as written. This is an internal algebraic inconsistency, not a disagreement with prior literature. The reader's uniformity concern is a modeling limitation, but the concrete DRT error is the more urgent correctness risk.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives the analytical solution of the time-fractional diffusion equation governing a finite-length resistor-CPE transmission line under galvanostatic charging. The main claimed result is the reduced impedance Z_TL(s_n)/(r_0 L)=s_n^{-alpha/2} coth(s_n^{alpha/2}) (Eq. 44), which the authors present as a physical derivation of Cabanel et al.'s empirical fractional finite-length Warburg impedance. The paper also gives an infinite-series expression for the voltage, a step-response formula, and a distribution of relaxation times g(tau) claimed to correspond to this impedance.","tokens_in":13234,"tokens_out":12072,"duration_ms":102345,"significance":"If the central derivation is correct, the paper provides a useful physical justification for a widely used empirical impedance model and supplies explicit closed-form series for the time-domain response of porous electrodes. The impedance part is based on standard Laplace and finite Fourier cosine transform steps and is likely sound once the r0=r issue is resolved. The distribution-of-relaxation-times part, however, contains an algebraic error that invalidates Eq. (53) as written. The paper's contribution is therefore significant but needs a substantive correction before the DRT claims can be accepted.","major_comments":[{"comment":"The substitution stated in the text does not produce Eq. (52). Starting from Eq. (51) with x=b s^alpha and u=b^{-1/alpha}v gives b^{1-1/alpha}s^{alpha-1}E_{alpha,alpha}(-b s^alpha) = (sin(alpha pi)/pi) integral_0^infty [b v^alpha e^{-s v}]/[b^2+2b v^alpha cos(alpha pi)+v^{2alpha}] dv, so that s^{alpha-1}E_{alpha,alpha}(-b s^alpha) = (sin(alpha pi)/pi) integral_0^infty [b^{1/alpha} v^alpha e^{-s v}]/[b^2+2b v^alpha cos(alpha pi)+v^{2alpha}] dv. The denominator in Eq. (52) is b^{1/alpha-1} times this correct denominator, and the numerator lacks the b^{1/alpha} factor; the two expressions agree only when b=1. Consequently Eq. (53) is not the distribution of relaxation times corresponding to Eq. (44)/(46), and the agreement with Eq. (57) shown in Fig. 4(c) cannot validate Eq. (53). Please correct the substitution and Eq. (53), or remove the DRT claim.","section":"Section II.C, Eq. (52)"},{"comment":"The model introduces r0 in the boundary condition Eq. (9) while the current-voltage relation Eq. (39) uses r, and the two are never related. Combining Eq. (40) at x=0 with Eq. (9) gives i(0,t)=(r0/r)i0, so the physical current at the pore mouth equals the applied current i0 only when r0=r. Since the impedance in Eq. (41) is defined through i0/s, the derivation of Eqs. (43)-(44) implicitly requires this equality. The paper should state r0=r explicitly or use a single resistance-per-unit-length parameter throughout.","section":"Section II.A, Eqs. (9), (39), (40)"}],"minor_comments":[{"comment":"For integer k, sin(k pi)=0, so the displayed expression can be simplified to the single term -r0 i0 L^2/(k^2 pi^2); presenting the simplification would make the k=0 case less opaque.","section":"Eq. (26)"},{"comment":"The caption groups magnitude and phase in panel (a) and labels the Nyquist plot as (b) and the DRT plot as (c), while the text refers to the phase plot as Fig. 4(b); these panel labels should be aligned.","section":"Fig. 4 caption and text"},{"comment":"The notation H_{p,q}^{m,n} is used without specifying p, q, m, and n; either define these parameters or explicitly state that the notation follows Ref. [39].","section":"Eq. (57)"},{"comment":"Reference [38] is a bare URL without a full bibliographic entry; please supply author, title, and access date.","section":"Reference [38]"},{"comment":"The title as posted contains an OCR artifact ('transm ission') that should be corrected to 'transmission'.","section":"Title"}],"recommendation":"major_revision","confidential_remarks":"The algebraic error in the DRT derivation is internal and fixable, and the impedance result appears sound after the r0=r assumption is made explicit, so the manuscript should be revisable within its scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe useful part of this paper is the time-domain solution. The authors solve the Caputo time-fractional diffusion equation for a finite-length R-CPE transmission line under galvanostatic charging and obtain the voltage series (Eq. 34), the current profile (Eq. 40), and the step response (Eq. 46). These are cleanly derived, the alpha=1 limit checks out, and the plots are sensible. If you work on impedance spectroscopy of porous electrodes, Eq. 46 is a practical closed form you might use.\n\nThe impedance result itself, Z(s_n)=s_n^{-alpha/2} coth(s_n^{alpha/2}), is not new. It is the anomalous finite-length Warburg impedance already derived by Bisquert and Compte (ref. 1), and the voltage series was given by Luchko (refs. 34-35). The paper cites these works, so it is honest, but the claim that this provides a 'physical derivation' of Cabanel et al.'s empirical form is overstated; that derivation already exists in the anomalous diffusion literature. What is new is framing it for the R-CPE transmission line and adding the time-domain formulas.\n\nTwo soft spots. First, the paper never states that r=r0. The boundary condition (Eq. 9) uses r0, but the current at x=0 from Eq. 40 is r0 i0/r. It equals i0 only if r=r0. The verification sentence mentions this in passing, but it should be a stated assumption in the model.\n\nSecond, and more serious: the DRT derivation in Section II-C contains an algebraic error. Substituting u=b^{-1/alpha}v into Eq. 51 and canceling gives the denominator b^2+2b v^alpha cos(alpha pi)+v^{2alpha}. Equation 52 instead has the denominator multiplied by b^{1/alpha-1}. So Eq. 53 does not follow from Eq. 51, and the 'excellent agreement' with Eq. 57 in Fig. 4(c) cannot hold as written. This is an internal inconsistency, not a disagreement with prior literature. The H-function expression (Eq. 57) may well be the correct DRT, but the Mittag-Leffler form (Eq. 53) is wrong.\n\nBottom line: the paper is worth sending to a referee because the time-domain solution is a legitimate contribution and the impedance derivation is correct. But the DRT section needs to be redone, and the r=r0 assumption needs to be explicit. As written, I would not cite Eq. 53, and I would treat Fig. 4(c)'s agreement claim with caution. With the algebra fixed, the paper would be a reasonable contribution to impedance spectroscopy.\n\nRecommendation: send to peer review, major revision required.","headline":"The time-domain solution is solid, but the impedance is already known and the DRT derivation has a correctable algebraic error.","tokens_in":13736,"tokens_out":4347,"would_cite":true,"duration_ms":36323,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35R11","26A33","44A10"],"pacs":[],"model":"deepseek-v4-flash","headline":"The bounded resistor-CPE transmission line has exact impedance $s_n^{-\\alpha/2}\\coth(s_n^{\\alpha/2})$, derived from the time-fractional diffusion equation.","keywords":["transmission line model","constant phase element","time-fractional diffusion","finite-length Warburg impedance","galvanostatic charging","Mittag-Leffler function","distribution of relaxation times","porous electrodes"],"falsifier":"Take a well-characterized uniform single pore with a blocked far end, measure its impedance over a wide frequency range, and independently fix $\\alpha$ from the low-frequency slope; if the full spectrum deviates from $s_n^{-\\alpha/2}\\coth(s_n^{\\alpha/2})$, or if the low-frequency real-axis intercept of the normalized impedance is not $1/3$, the central derivation is falsified. A galvanostatic step measured on the same pore should also match the Mittag-Leffler series solution for the voltage at $x=0$.","tokens_in":12646,"feed_emoji":"⚡","tokens_out":13031,"duration_ms":109425,"temperature":0.7,"pith_summary":"This paper derives the exact impedance of a finite-length transmission line in which each unit length of an electrode pore carries a constant resistance in series with a constant-phase element (CPE) rather than an ideal capacitor. Solving the resulting time-fractional diffusion equation under galvanostatic charging, with a perfectly blocking end at $x=L$, gives the reduced impedance $Z_{\\mathrm{TL}}(s_n)/(r_0 L)=s_n^{-\\alpha/2}\\coth(s_n^{\\alpha/2})$. This is exactly the fractional modification of the reflective finite-length Warburg impedance that has long been used, on empirical grounds, to fit porous electrodes, supercapacitors, and insertion electrodes. The derivation supplies that empirical formula with a physical starting point, and it also yields the voltage and current step responses and an analytical distribution of relaxation times. If the derivation is right, one dispersion exponent $\\alpha$ and the pore's resistance and length set the entire frequency- and time-domain response of a uniform blocking pore.","feed_headline":"Fractional diffusion yields Warburg's empirical impedance","feed_subtitle":"Derived from the time-fractional diffusion equation, it justifies an empirical Warburg model used in porous electrodes.","key_machinery":"The load-bearing object is the time-fractional diffusion equation for the voltage $v(x,t)$, with the Caputo fractional derivative of order $\\alpha$ replacing the ordinary time derivative; the constant-phase element (CPE, impedance $1/(c_\\alpha s^\\alpha)$) is the energy-storage element per unit length. The equation is solved by applying a finite Fourier cosine transform in space and a Laplace transform in time, which turns the problem into algebraic equations whose inversion uses Mittag-Leffler inverse Laplace formulas. The impedance is then the ratio of the Laplace-transformed surface voltage to the Laplace-transformed input current, and the reflection at the blocked end $x=L$ produces the hyperbolic cotangent, while the fractional CPE law produces the exponent $\\alpha/2$.","core_discovery":"The central claim is that the impedance of a bounded, one-dimensional resistor-CPE transmission line is exactly $$\\frac{Z_{\\mathrm{TL}}(s_n)}{r_0 L}=$s_n^{{-\\alpha/2}}$\\coth\\!\\left($s_n^{{\\alpha/2}}$\\right),$$ with dimensionless frequency $s_n=s\\left(r c_\\alpha L^2\\right)^{1/\\alpha}$, obtained by solving the Caputo time-fractional diffusion equation $${}_0D_t^\\$\\alpha$ v=\\frac{1}{r c_\\$\\alpha$}\\frac{\\$partial^{2}$ v}{\\partial $x^{2}$},\\quad 0<x<L,$$ with the galvanostatic boundary condition $-r_0^{-1}\\partial_x v|_{x=0}=i_0$ and the blocking condition $\\partial_x v|_{x=L}=0$. The same solution gives the voltage as an infinite Mittag-Leffler series, the step response through one further time derivative, and an RC relaxation-time distribution $g(\\tau)$ such that the impedance is the integral of $g(\\tau)/(1+s\\tau)$. For $\\alpha=1$ the formulas reduce to the classical reflective finite-length Warburg impedance of an ideal RC line, and the paper notes that this limit corrects missing terms in earlier galvanostatic solutions.","pith_inferences":["A direct cross-check the paper does not perform: fix $\\alpha$ from the low-frequency CPE slope of a uniform pore, then predict the whole impedance and step response with no additional free parameters and compare; agreement would be a strong test of the derivation.","The same fractional diffusion equation with potentiostatic or permeable-end boundary conditions should produce analogous $\\tanh$- or $\\sinh$-type impedance forms, so the derivation likely extends beyond the perfectly blocking galvanostatic case.","The natural time scale $t/(r c_\\alpha L^2)^{1/\\alpha}$ predicts that galvanostatic voltage curves for pores of different lengths or electrolyte conductivities collapse onto one master curve when rescaled, a testable scaling law.","The model suggests interpreting $\\alpha$ as a property of the distributed electrode/electrolyte interface rather than a pure fit parameter; one could test this by comparing $\\alpha$ from impedance fits with $\\alpha$ obtained from independent time-domain fractional-capacitance measurements on the same material."],"forward_implications":["The empirical fractional finite-length Warburg impedance now follows from a first-principles fractional diffusion equation, so porous-electrode fits using $s_n^{-\\alpha/2}\\coth(s_n^{\\alpha/2})$ have a stated physical model behind them.","The same dispersion exponent $\\alpha$ governs both the high-frequency half-order CPE branch and the low-frequency $s_n^{-\\alpha}$ branch, linking the two regions of a Nyquist plot to one parameter.","The voltage solution gives the full galvanostatic charging curve, so chronopotentiometry and impedance spectroscopy can be analyzed with the same set of parameters $r$, $c_\\alpha$, $L$, and $\\alpha$.","The analytical distribution of relaxation times makes it possible to compute the RC time-constant spectrum of a bounded CPE pore directly, without numerical deconvolution schemes.","For $\\alpha=1$, all derived formulas reduce to the classical reflective finite-length Warburg/RC transmission-line results, with the earlier ideal-capacitor galvanostatic solution corrected."],"supporting_citations":[{"why":"Supplies the standard reflective finite-length Warburg impedance form that the fractional model generalizes.","marker":"[1]"},{"why":"Proposes the empirical $s_n^{-\\alpha/2}\\coth(s_n^{\\alpha/2})$ modification that this paper derives from a physical model.","marker":"[11]"},{"why":"Provides the transmission-line diffusion equation and the impedance definition via Laplace-transformed voltage and current at the pore mouth.","marker":"[14]"},{"why":"Gives the semi-infinite R-CPE transmission-line solution, including the half-order CPE impedance, that the bounded-domain analysis extends.","marker":"[15]"},{"why":"Supplies the combined finite Fourier cosine transform and Laplace transform method used to solve the fractional diffusion equation.","marker":"[26]"},{"why":"Gives the ideal RC bounded transmission-line galvanostatic solution whose $\\alpha=1$ limit the paper generalizes and corrects.","marker":"[27]"},{"why":"Provides the inverse Laplace transform formula used to obtain Mittag-Leffler series solutions from the transformed equations.","marker":"[32]"},{"why":"Supplies the procedure for deriving an analytical distribution of relaxation times from impedance functions via Fox H-functions.","marker":"[39]"}],"fun_headline_variants":["Fractional diffusion yields exact Warburg impedance","Time-fractional diffusion solves reflective Warburg model","Exact impedance for bounded resistor-CPE line","Resistor-CPE line solved exactly via fractional diffusion","New analytical proof for fractal Warburg impedance"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation assumes the electrode pore is a single uniform tube whose electrical resistance per unit length and whose capacitive response, described by a constant-phase element with a fixed exponent $\\alpha$, do not vary along its length, and whose far end blocks all current; if any of these vary with position or state of charge, the exact impedance formula and time-domain solution no longer apply.","fun_headline_variants_meta":{"raw":{"variants":["Fractional diffusion yields exact Warburg impedance","Time-fractional diffusion solves reflective Warburg model","Exact impedance for bounded resistor-CPE line","Resistor-CPE line solved exactly via fractional diffusion","New analytical proof for fractal Warburg impedance"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001016,"raw_usage":{"total_tokens":4324,"prompt_tokens":1011,"completion_tokens":3313,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":627,"completion_tokens_details":{"reasoning_tokens":3242}},"tokens_in":627,"tokens_out":3313,"duration_ms":22061,"temperature":1.0,"reasoning_tokens":3242,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:17:01.444297+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a well-characterized uniform single pore with a blocked far end, measure its impedance over a wide frequency range, and independently fix $\\alpha$ from the low-frequency slope; if the full spectrum deviates from $s_n^{-\\alpha/2}\\coth(s_n^{\\alpha/2})$, or if the low-frequency real-axis intercept of the normalized impedance is not $1/3$, the central derivation is falsified. A galvanostatic step measured on the same pore should also match the Mittag-Leffler series solution for the voltage at $x=0$.","supporting_citations":[{"cited_title":"Bisquert and A","cited_arxiv_id":null,"evidence_quote":"Supplies the standard reflective finite-length Warburg impedance form that the fractional model generalizes."},{"cited_title":"Cabanel, G","cited_arxiv_id":null,"evidence_quote":"Proposes the empirical $s_n^{-\\alpha/2}\\coth(s_n^{\\alpha/2})$ modification that this paper derives from a physical model."},{"cited_title":"Pedersen, T","cited_arxiv_id":null,"evidence_quote":"Provides the transmission-line diffusion equation and the impedance definition via Laplace-transformed voltage and current at the pore mouth."},{"cited_title":"Allagui and E","cited_arxiv_id":null,"evidence_quote":"Gives the semi-infinite R-CPE transmission-line solution, including the half-order CPE impedance, that the bounded-domain analysis extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the combined finite Fourier cosine transform and Laplace transform method used to solve the fractional diffusion equation."},{"cited_title":"Posey and T","cited_arxiv_id":null,"evidence_quote":"Gives the ideal RC bounded transmission-line galvanostatic solution whose $\\alpha=1$ limit the paper generalizes and corrects."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the inverse Laplace transform formula used to obtain Mittag-Leffler series solutions from the transformed equations."},{"cited_title":"Allagui and A","cited_arxiv_id":null,"evidence_quote":"Supplies the procedure for deriving an analytical distribution of relaxation times from impedance functions via Fox H-functions."}],"review_version":1}