{"id":"2752d73d-abdc-45f2-a4dd-7dad2e12e1ea","arxiv_id":"2603.28799","paper_version":3,"verdict":"REJECT","confidence":"HIGH","novelty_score":2.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"PID control of tritium breeding in liquid-lithium fusion systems is shown to be expressible, after linearization, as a localized Bessel-type differential equation with an explicit map between PID gains and Bessel parameters.","lead":"This paper proposes that PID controllers for tritium breeding and liquid-lithium jet cooling can be mathematically recast as localized Bessel-type operators. The claimed benefit is a compact analytical framework for designing feedback control in fusion blanket and jet systems.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The PID-to-Bessel map depends on an asserted error dynamics: Eq. (21) is not derived from the plant, so Eq. (23) does not connect PID tuning to Bessel physics.","rationale":"The reader's weakest_assumption is exactly the load-bearing point: Eq. (21) is not derived from the jet/blanket physics. The central claim is a formal correspondence between PID gains and Bessel parameters. The mapping (23) follows algebraically if Eq. (21) holds, but Eq. (21) itself is an assumption. The paper's only validation (Section VI) fits a second-order error model to simulated closed-loop data and then maps the fitted coefficients; it does not independently derive those coefficients from the plant equations. The paper's limitations (Section VII) concede the models are simplified and the mapping is local, but do not acknowledge that the specific coefficient identification (21) is unjustified. A concrete first-order plant model with the PID law yields closed-loop error coefficients containing the plant time constant, and they match Eq. (21) only under degenerate gain/plant constraints. Therefore the central physical claim is unsupported; the individual algebraic steps after Eq. (21) are fine. I agree with the reader's REJECT verdict, so no verdict change is needed.","tokens_in":11581,"tokens_out":7390,"duration_ms":70083,"concrete_test":"Using the lumped jet model from Section VI (first-order relaxation τ dT/dt + T = q_in + u), close the loop with the PID law u = aE + b dE/dt + c∫E, E = T_target - T. Symbolically derive the exact closed-loop homogeneous error equation and compare its coefficients with Eq. (21). For the canonical first-order plant the coefficients are (1+a)/(τ+b) and c/(τ+b), which coincide with a/b and c/b only for c=0 and τ=b/a; if the published simulation uses a different plant, use that plant and recompute the fitted a1, a0. If they do not equal a/b and c/b, Eq. (21) is imposed rather than derived and Eq. (23) does not establish the PID-Bessel correspondence.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. (21) is the load-bearing step: it converts the PID increment (15) into the second-order error model (20). But (15) is a control signal, not a dynamical law for E(t); no plant model is specified, so a1 = a/b and a0 = c/b are asserted. For the first-order lumped jet model used in Section VI (τ dT/dt + T = q_in + u, with u = aE + bE' + c∫E and E = T_target - T), the closed-loop homogeneous error equation is (τ+b)E'' + (1+a)E' + cE = 0, i.e. E'' + [(1+a)/(τ+b)]E' + [c/(τ+b)]E = 0, not Eq. (21). Agreement with Eq. (21) would require c=0 and τ=0 simultaneously (or τ=b/a when c=0), so the correspondence (23) is not a generic consequence of PID control of the jet. The numerical experiment in Section VI cannot rescue this: it injects Bessel-modulated sources and fits the assumed second-order template, validating the fitting procedure rather than the derivation of (21). Thus the central claim — that PID-tuned error dynamics can be interpreted as localized Bessel modes — rests on an error equation chosen to make the map work, not on the jet/blanket physics. The paper's own limitation statements acknowledge reduced fidelity but do not flag this missing derivation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a unified analytical framework connecting tritium breeding in lithium-based fusion blankets, thermal response of a liquid-lithium jet, and PID feedback control. The authors derive a reduced jet thermal-expansion model from mass and heat transport equations, express PID control as a linear operator, and define a local correspondence between the PID increment and a Bessel-type differential operator acting on the tritium-inventory error. The central claim is that any PID-tuned second-order error dynamics can be interpreted as a localized Bessel mode via the explicit map x0 = b/a, ν² = (1 − c/b) b²/a² (Eq. 23). The paper illustrates this mapping with a numerical jet/blanket–PID simulation and discusses limitations.","tokens_in":12039,"tokens_out":2648,"duration_ms":28132,"significance":"If the central claim were valid, the paper would provide a genuinely compact analytical bridge between neutronics, thermohydraulics, and control in Li-based fusion systems, complementing purely numerical PID-tuning studies. The manuscript is clearly structured, uses explicit operator notation, and the algebraic steps from Eqs. (15)–(23) are internally consistent. The numerical experiments in Section VI are described in sufficient detail that the fits in Figs. 4–5 are plausible and reproducible. However, the central identification is not actually derived from the fusion or jet physics: it rests on an assumed second-order error model (Eq. 21) whose coefficients are set to PID gain ratios, and the numerical validation fits the same second-order template. The paper's own limitation section acknowledges reduced fidelity of the thermohydraulic and TBR models, but it does not flag the missing derivation of the error equation. Because the core contribution depends on this unproven and, in the specific plant model used, contradicted assumption, the claimed unification is not established.","major_comments":[{"comment":"Equation (21) is assumed, not derived. It asserts that the closed-loop error obeys E'' + (a/b)E' + (c/b)E = 0, but the PID increment (15) is a control signal, not a dynamical law for E(t). For the first-order lumped jet model actually used in Section VI, τ dT/dt + T = q_in + u with u = aE + bE' + c∫E and E = T_target − T, the homogeneous closed-loop error equation is (τ+b)E'' + (1+a)E' + cE = 0, i.e., E'' + [(1+a)/(τ+b)]E' + [c/(τ+b)]E = 0, not Eq. (21). Agreement with Eq. (21) would require c = 0 and τ = 0 simultaneously (or τ = b/a when c = 0). Thus Eq. (23) is not a generic consequence of PID control of a jet; it is an imposed ansatz.","section":"§IV.B, Eq. (21)"},{"comment":"The PID-to-Bessel map is circular. Equation (20) postulates a second-order error model with coefficients a0 and a1, and Eq. (21) identifies those coefficients with PID gain ratios. Equation (23) then derives (x0, ν) from that same assumed model. The claimed 'prediction' of Bessel-type behavior is therefore built into the construction. The numerical example in Section VI and Figs. 4–5 does not break the circularity: it generates error trajectories from a Bessel-modulated source and then fits E'' + a1E' + a0E ≈ 0, confirming that the simulated data can be fit by the assumed second-order template. This validates the fitting procedure, not the derivation of Eq. (21) from the jet/blanket physics.","section":"§IV.B, Eqs. (20)–(23)"},{"comment":"The TBR connection depends on the linear surrogate TBR(t) ≈ αR(t) + β, where α and β are free fit parameters. The paper acknowledges that α and β are local effective parameters, but it does not provide any neutronics calculation or data supporting the chosen values (α = 0.08, β = 1.0). Since the central claim involves 'lithium-based tritium breeding,' the absence of a concrete neutronics link weakens the claim that the PID–Bessel dynamics are tied to TBR physics rather than to a generic linear model.","section":"§V, Eq. (28)"},{"comment":"The jet thermal-expansion model is not actually solved. Section III derives Eqs. (5)–(6) and presents Figure 3, but no closed-form or numerical solution of those partial differential equations is given. Section VI substitutes a lumped-parameter jet model with a Bessel-modulated heat source q_in(t), which is not derived from the PDEs and does not test the earlier reduction chain. The numerical experiment therefore does not provide independent support for the claimed coupling between jet thermohydraulics and the PID–Bessel correspondence.","section":"§III and §VI"}],"minor_comments":[{"comment":"The sentence 'Our results indicate that Li-based breeding and heat-removal systems exhibit low-order.' is grammatically incomplete. Also, 'fast' is typo for 'first' and 'Ttarge' for 'T_target' in the caption of Figure 5.","section":"Abstract / Introduction"},{"comment":"The caption 'Different values for the maximum velocity vz...' is vague; the figure lacks axis labels and a description of which parameter is varied, making it difficult to interpret.","section":"Figure 3"},{"comment":"Some references are in nonstandard format (e.g., [2], [3], [14]) and several URLs are incomplete or missing access dates. The manuscript would benefit from a uniform citation style.","section":"References"},{"comment":"The discussion of the analogy between PID and Bessel shift operators is purely structural; it may help the reader to explicitly state that no physical equivalence is being claimed beyond the local second-order identification.","section":"§IV.A"}],"recommendation":"reject","confidential_remarks":"The paper attempts to bridge two interesting areas, but the central identification is an assumed ansatz rather than a derived result, and the numerical validation is circular. The load-bearing error in Eq. (21) and the circular fit in Section VI cannot be fixed by minor revision; they require either a genuine plant-derived derivation of the second-order error dynamics or a rewording of the claim to a much weaker 'one can always choose a second-order error model' statement, which would eliminate the novelty. I therefore recommend rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague —\n\nThis is a short, clearly written paper that works the analogy between PID control and localized Bessel operators in the context of lithium-based fusion systems. The algebraic steps are consistent, and the authors are upfront that the blanket/jet model is a reduced-order conceptual demonstration, not a realistic design. The literature review is serviceable, and the operator parallel (PID as a linear combination of I, identity, and D; Bessel recurrences as linear combinations of d/dx and multiplication by ν/x) is drawn neatly. So there is something here, but it is not a new result.\n\nThe central mapping is Eq. (23): x0 = b/a, ν² = (1 – c/b)b²/a². It comes from equating the second-order error model (20) to the localized Bessel equation (22). The problem is that (20) is assumed, not derived. The paper says, without a plant model, that the physical system responds approximately proportionally to the PID increment, and the integral term is a slow drift. That is load-bearing. For the lumped jet model actually used in Sec. VI, the closed-loop homogeneous error equation is (τ+b)E'' + (1+a)E' + cE = 0, which reduces to (20) only under special conditions (e.g., τ = 0 and c = 0) that are not the generic case. So the correspondence does not follow from PID control of the jet; it is imposed by choosing the error equation that makes the map work.\n\nThe numerical validation in Sec. VI does not fix this. The authors inject a Bessel-modulated heat source, run the lumped model, then fit the same second-order template to the error. The close overlap confirms the fit, not the derivation. The TBR connection adds fitted linear parameters α and β, which are acknowledged to be local surrogates. So the whole construction is a re-parameterization: any constant-coefficient second-order ODE can be written in the localized Bessel form. That is true by coefficient matching, but it carries no new physics and no new control insight.\n\nWhat the paper does well is scope itself honestly. The discussion lists the right limitations: linearization, local validity, scalar variables, and the need for higher-fidelity checks. The authors are not claiming more than a conceptual framework. But the conceptual framework is too thin to support the title's promise. The mapping is an exercise in relabeling unless the error dynamics are derived from the plant.\n\nWho should read this? Someone thinking about abstract operator analogies for control might find it a neat illustration; someone working on fusion blanket control will not take practical guidance from it. It could serve as a cautionary example of circular validation in a reading group, but I would not cite it.\n\nFor peer review, I would desk reject. The core result is not new and the load-bearing assumption is exactly where the physics should be. The authors would need to derive the error dynamics from a specific plant model, or clearly present this as a purely formal analogy without claiming physical content.","headline":"A cleanly written but essentially circular exercise: the PID-to-Bessel map rests on an assumed error equation, so the central claim is a re-parameterization rather than a result.","tokens_in":12447,"tokens_out":3123,"would_cite":false,"duration_ms":31802,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["34H05","33C10"],"pacs":[],"model":"deepseek-v4-flash","headline":"A PID feedback law on tritium-inventory error is shown to be locally equivalent to a Bessel-type differential operator, giving a unified analytic view of lithium-based breeding and heat-removal control.","keywords":["tritium breeding","PID control","Bessel functions","lithium blanket","liquid lithium jet","feedback control","operator theory","fusion heat removal"],"falsifier":"Run a high-fidelity simulation or experiment of a lithium jet with a PID controller on tritium inventory, perturb the setpoint, and record the error. If the error trajectory is not well fit by a second-order linear ODE with coefficients a/b and c/b, or if the fitted Bessel parameters deviate beyond tolerance from (23), the claim is falsified. In particular, if the residual of the second-order fit grows with time or with perturbation amplitude, the local equivalence fails.","tokens_in":11529,"feed_emoji":"🌀","tokens_out":6232,"duration_ms":58771,"temperature":0.7,"pith_summary":"The paper claims that a standard continuous-time PID controller applied to tritium-inventory error in lithium-based fusion systems can be embedded as a special case of a Bessel-type differential operator after localizing around a reference operating point. The explicit parameter map x0=b/a, ν²=(1−c/b)(b/a)² turns PID gains (a,b,c) into effective Bessel parameters (x0,ν). If correct, every PID-tuned second-order error trajectory is a localized Bessel mode, connecting tritium breeding, jet thermal response, and controller design in one modal language. The authors support this with a reduced liquid-lithium jet thermal model and numerical fits showing the error dynamics follow a second-order LTI form. The contribution is an analytical bridge rather than a full reactor model.","feed_headline":"PID tuning becomes Bessel-mode selection in lithium fusion control","feed_subtitle":"A single modal picture links tritium breeding, jet heat removal, and controller design.","key_machinery":"The argument rests on comparing the PID operator C=K(I+(1/Ti)∫+Td d/dt) with the Bessel shift/recurrence operator and the full Bessel differential operator Lν=x²(d²/dx²)+x(d/dx)+(x²−ν²). The load-bearing step is localizing Lν around a reference point x0 (freezing coefficients) and rescaling x=αt, which converts it into a constant-coefficient second-order operator. Combining this with a second-order error model yields the explicit map from PID gains (a,b,c) to Bessel parameters (x0,ν), realizing PID as a localized Bessel mode.","core_discovery":"The central claim is that the PID increment Δr∞(t)=aE(t)+bE'(t)+c∫E(τ)dτ, acting on the error trajectory E(t), can be identified with a localized Bessel-type operator. Around a reference point x0, the Bessel operator x²(d²/dx²)+x(d/dx)+(x²−ν²) approximates to a constant-coefficient second-order operator. Under the time scaling x=αt and matching coefficients with the closed-loop error equation E''+(a/b)E'+(c/b)E=0, the authors derive the explicit correspondence x0=b/a and ν²=(1−c/b)(b²/a²). Thus the PID controller's gains select an effective Bessel order and radius, so the closed-loop error behaves locally like a Bessel mode of order ν near x0. The equivalence is local, not global: the spectr","pith_inferences":["A natural extension is to ask whether a multivariable PID controller on several errors (tritium inventory, temperature, impurity level) maps to a product or tensor of Bessel operators; if so, mode-decoupling may be possible.","The local nature suggests the identification may fail during large transients or actuator saturation; a gain-scheduled PID that preserves the map would be a testable follow-up.","Since the mapping only uses the second-order error model, it likely generalizes to any plant with approximately second-order linearized error dynamics, not just lithium fusion systems; this could be checked on standard benchmark control problems.","The Bessel-order parameter ν acts like a tunable damping/frequency knob; linking it to physical quantities (e.g., 6Li enrichment sensitivity) could yield engineering guidelines for choosing a and b."],"forward_implications":["PID tuning for tritium inventory can be reinterpreted as selecting an effective Bessel order and radius, enabling mode-based insight into error dynamics.","The local equivalence means different PID gain sets correspond to different Bessel modes, so controller design can be guided by Bessel-function properties.","The reduced jet model suggests beam-induced thermal perturbations are captured by a few scalar observables, so low-order controllers suffice near an operating point.","If the equivalence holds, the same analytical language applies to tritium breeding and heat-removal control, possibly simplifying multi-loop design in lithium systems.","The mapping provides a testable prediction: closed-loop error should follow E''+(a/b)E'+(c/b)E=0, with the Bessel parameters fixed by (23)."],"fun_headline_variants":["PID gains map to Bessel order and radius in fusion blanket","Fusion control: PID gains pick Bessel mode parameters","Lithium fusion control: PID tuning equals Bessel selection","Bessel modes from PID gains in lithium fusion systems","Linking PID control to Bessel operators for lithium fusion"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The entire correspondence rests on the unproven assumption that the closed-loop error obeys E'' + (a/b)E' + (c/b)E = 0, justified by asserting the plant responds proportionally to the PID increment on a short timescale and the integral term is a slow drift; this error equation is not derived from the jet or blanket physics.","fun_headline_variants_meta":{"raw":{"variants":["PID gains map to Bessel order and radius in fusion blanket","Fusion control: PID gains pick Bessel mode parameters","Lithium fusion control: PID tuning equals Bessel selection","Bessel modes from PID gains in lithium fusion systems","Linking PID control to Bessel operators for lithium fusion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000471,"raw_usage":{"total_tokens":2234,"prompt_tokens":851,"completion_tokens":1383,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":595,"completion_tokens_details":{"reasoning_tokens":1302}},"tokens_in":595,"tokens_out":1383,"duration_ms":10792,"temperature":1.0,"reasoning_tokens":1302,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T05:38:08.673163+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a high-fidelity simulation or experiment of a lithium jet with a PID controller on tritium inventory, perturb the setpoint, and record the error. If the error trajectory is not well fit by a second-order linear ODE with coefficients a/b and c/b, or if the fitted Bessel parameters deviate beyond tolerance from (23), the claim is falsified. In particular, if the residual of the second-order fit grows with time or with perturbation amplitude, the local equivalence fails.","supporting_citations":[],"review_version":1}