{"id":"423894df-5c2d-4b74-a4e8-18fc6b4aec3f","arxiv_id":"2507.21365","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A new COSY-Infinity code, GES, for electrostatic toroidal bends with fringe fields agrees with GIOS and with the author's hard-edge theory, and shows that COSY's ESP routine misses a curvature-derivative fringe-field term.","lead":"A TRIUMF physicist reports that a new COSY-Infinity routine, GES, computes higher-order beam maps for electrostatic toroidal bends with fringe fields, and that this routine matches the independent GIOS code but not COSY's built-in ESP routine. The likely reason is that ESP omits a term involving the derivative of the bending curvature in fringe regions, a second-order effect that matters for beam optics design.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The four 'mysteriously' disagreeing second-order dispersive terms between GES and GIOS at finite aperture leave the joint validation of the hard-edge map incomplete.","rationale":"The h''-zeroing experiment in Section 4.4 is strong: when h'' is removed from GES, its second-order map becomes virtually identical to ESP, so the specific claim that ESP misses the curvature-derivative term is well supported even without inspecting ESP internals. The theory behind Eqs. (8)-(9) is also consistent with the fringe-field-alone test in Section 4.2. The soft spot is the finite-aperture comparison with GIOS: the four unexplained dispersive terms mean the two 'correct' codes do not fully agree where the finite fringe shape matters. This is exactly the regime where the paper's universality claim for the hard-edge map is needed and least tested. I also note a secondary algebraic issue in Section 2.2: Eq. (7) for a12 does not follow from the recursion relation as written and appears inconsistent with the neighboring third-order coefficients; this does not affect the code or the hard-edge map, but it lowers confidence in the presentation of the theory. On balance, the central diagnosis is likely correct, so the reader's CONDITIONAL verdict remains appropriate.","tokens_in":11080,"tokens_out":34545,"duration_ms":407338,"concrete_test":"Re-run the Section 4.1 comparison with GES and GIOS using the same finite aperture, the same Enge falloff coefficients, and the same effective fringe length, and inspect all second-order coefficients, especially (X,AD), (A,XD), (Y,BD), (B,YD). Independently, build a 3D Poisson-solver model of the same toroidal bender ends and extract entry/exit fringe maps. If the four terms converge to the Eq. (8)-(9) hard-edge values at the few-percent level or better, the concern is resolved; if they persist or the 3D maps contain extra y-dependent second-order terms, the hard-edge map is not the complete second-order fringe-field map and the GES/GIOS/ESP diagnosis needs qualification.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that GES and GIOS agree with the hard-edge map while ESP lacks the h'' term. Section 4.1's finite-aperture comparison is the main joint validation, but it contains four second-order dispersive terms, (X,AD), (A,XD), (Y,BD), (B,YD), that the author calls 'mysteriously in disagreement.' The few-percent differences in the other terms are attributed to GIOS's default fringe integrals, but these four are not explained, and no matched-parameter rerun is provided. The later zero-aperture agreement only shows that the hard-edge limit is common; it removes exactly the finite-fringe effects that could generate shape-dependent second-order terms. Thus the paper's strongest external confirmation that the hard-edge map is the complete irreducible second-order effect is incomplete at realistic apertures, and the claim that ESP disagrees with 'GES and GIOS' is fully established only in the zero-aperture limit.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents GES, a new COSY-Infinity element for electrostatic toroidal benders, and uses it to argue that the hard-edge fringe-field map for such a bender has second-order terms Δx = x²/(2ρ) and ΔPx = -xPx/ρ at entry, with signs reversed at exit, arising from the h'' curvature-derivative term in the potential. The paper compares GES with the built-in ESP/ECL elements and with the external code GIOS at finite and zero aperture, tests the map in the fringe-field-only and thin-lens limits, and performs a controlled experiment in which the h'' term is removed from GES, showing that the resulting map matches ESP. It concludes that ESP omits the curvature-derivative term, which explains the long-standing discrepancy.","tokens_in":11255,"tokens_out":9882,"duration_ms":100318,"significance":"If the claims hold, the paper resolves a practical puzzle in charged-particle optics and supplies a simple, parameter-free hard-edge map that should be included by default in transport codes. Strengths include the reproducible COSY code in the appendix, the comparison against the independent GIOS code, the clean thin-lens-limit checks, and the controlled h''-removal experiment. The paper also makes a falsifiable prediction: the hard-edge map is independent of fringe-field falloff shape. However, because four dispersive second-order terms in the finite-aperture GIOS comparison remain unexplained, the central claim is not yet fully established at realistic apertures.","major_comments":[{"comment":"The finite-aperture comparison of GES and GIOS lists four second-order dispersive terms, (X,AD), (A,XD), (Y,BD), and (B,YD), as 'mysteriously in disagreement.' These are not small residuals: for example, GES (X,AD) = -0.1466 while GIOS (X,AD) = 0.2071. Since this comparison is the main independent check at a realistic aperture, leaving these discrepancies unexplained means the claim that GES and GIOS agree substantially is unsupported for exactly the terms that test dispersive fringe-field coupling. Please either match the GIOS fringe-field integrals to GES's Enge parameters and rerun, or explain the source of these terms and state what the agreement claim does and does not cover.","section":"Section 4.1"},{"comment":"The conclusion that ESP lacks the h'' curvature-derivative term is inferred indirectly: GES with hpp set to zero reproduces ESP to 10^-4. This is a strong controlled experiment, but it would be more conclusive to identify the omission directly in the POTXZ/ESP potential or Hamiltonian. Please either inspect the ESP source and point to the absent h''-dependent term, or explicitly label this as an inference from a numerical null experiment.","section":"Section 4.4"},{"comment":"The derivation of the x-shift in Eq. (9), Δx = x²/(2ρ), is referenced to Eq. (20) of the author's earlier note [7] but not shown. Since this hard-edge map is the central theoretical result and is used to interpret all later comparisons, please include the canonical transformation or a short derivation so the paper is self-contained.","section":"Section 2.3"},{"comment":"The paper's premise that the lowest-order irreducible fringe-field effect is independent of the falloff shape is taken from the author's earlier Snowmass talk [3] and is not tested here. The zero-aperture test only isolates the hard-edge limit; it does not probe shape dependence at finite aperture. Given that the finite-aperture GIOS comparison has unresolved terms (Major comment 1), the shape-independence assumption should be tested, for example by varying the Enge coefficients in GES and checking that the second-order map changes only by terms beyond the hard-edge map.","section":"Section 1 / Ref. [3]"}],"minor_comments":[{"comment":"The phrase 'theCOSY code' should read 'the COSY code'.","section":"Abstract / page 2"},{"comment":"The asterisks in the fringe-field-alone output are not explained in the text; state explicitly that the two starred entries correspond to the theoretical values of Δx/x² and ΔPx/(xPx).","section":"Section 4.2"},{"comment":"The reference to Valetov et al. [9] would be more useful if it named the specific claim or section that the present results contradict.","section":"Section 4.1"},{"comment":"The last expression, Px/(1 ± x/ρ), appears to have a sign opposite to the series expansion just above it; please check the sign convention.","section":"Eq. (10)"},{"comment":"The use of the Enge function in HVBEND, form(s) = enge(1,1,2,-10*s+5), is not documented; adding the meaning of the four arguments and the role of lenfr would aid reproducibility.","section":"Appendix A"}],"recommendation":"major_revision","confidential_remarks":"The paper is largely self-referential in that the theory, the code, and several key references are all by the same author, but the comparison with the external code GIOS and the controlled h''-removal experiment provide independent grounding. The main risk is the four unexplained finite-aperture discrepancies; please ask the author to address them before acceptance. The paper may also need a statement about the scope of the agreement with GIOS."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The punchline: this note gives an explicit second-order hard-edge fringe-field map for electrostatic toroidal benders, Eqs. (8)-(9), and traces a long-standing COSY/ESP vs GIOS discrepancy to ESP's omission of curvature-derivative terms. That diagnosis is well supported by the numerical tests, especially the h''-removal test in Sec. 4.4, where GES without h'' matches ESP almost exactly. The thin-lens limits also agree with the theory. The code is shipped in the appendix, and the symplectic error checks are reported, so the evidence is reproducible.\n\nThe main soft spot is the four second-order dispersive terms that disagree between GES and GIOS at finite aperture, which the author himself calls mysterious. The zero-aperture and thin-lens tests cleanly validate the hard-edge limit, but they remove exactly the finite-fringe shape effects that could hide a shape-dependent term. So the claim that the hard-edge map is the complete irreducible second-order effect for realistic apertures is not fully nailed down. That said, the few-percent agreement with GIOS on all other terms and the h''-removal test mean this is a minor-to-moderate caveat, not a load-bearing flaw.\n\nOne more thing: the paper is self-referential, but the core comparison against GIOS and the internal h''-removal test give independent grounding. No fitted parameters are in play. The citation pattern is appropriate.\n\nWho is this for? Practitioners using COSY or GIOS for electrostatic bend design, and anybody who has hit the ESP/GIOS mismatch. It is a focused technical note, not a broad review article, but it is worth a serious referee because it resolves a real, long-standing code inconsistency and provides a simple map that could be added to transport codes. My recommendation: send it to peer review, ask the author to either explain or resolve the four finite-aperture dispersive discrepancies, and make the zero-aperture comparison more prominent as the definitive validation.","headline":"A short, well-grounded technical note that likely resolves why COSY's ESP disagrees with GIOS on electrostatic bender fringe fields; the hard-edge map is new and useful, but a few unexplained second-order terms remain.","tokens_in":11763,"tokens_out":835,"would_cite":true,"duration_ms":12002,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Electrostatic toroidal benders carry an irreducible hard-edge fringe-field map, $\\Delta x = x^2/(2\\rho)$ and $\\Delta P_x = -xP_x/\\rho$, that the standard built-in element omits.","keywords":["electrostatic toroidal bender","fringe-field effects","hard-edge transfer map","second-order aberrations","curvature derivative terms","symplectic transfer maps","beam optics","GES element"],"falsifier":"Track a charged particle through a numerically computed or measured electrostatic field of a toroidal bender with a smooth, finite-length fringe region, and extract the entry-map coefficients $(x|xx)$ and $(x|xp)$; if they deviate from $1/(2\\rho)$ and $-1/\\rho$ in a shape-dependent way beyond numerical error, the hard-edge universality claim fails.","tokens_in":10868,"feed_emoji":"⚡","tokens_out":11412,"duration_ms":118254,"temperature":0.7,"pith_summary":"An electrostatic toroidal bender has an intrinsic second-order fringe-field effect that no amount of end-face shaping can remove. In the hard-edge limit, the entry fringe field shifts the transverse coordinate and momentum by $\\Delta x = x^2/(2\\rho)$ and $\\Delta P_x = -xP_x/\\rho$, with both signs reversed at the exit. The paper traces this to the curvature-derivative term $-h'' x^3/6$ in the potential expansion, and shows that the built-in ESP element omits it while the new GES element reproduces it. If the claim is right, beam-optics codes that ignore these maps misplace second-order bend aberrations even when they idealize the fringe field.","feed_headline":"Bender fringe fields kick x by x²/2ρ and Px by −xPx/ρ","feed_subtitle":"A hard-edge map shows the second-order bend aberration is intrinsic—and the standard built-in element misses it.","key_machinery":"The load-bearing object is the curvature-derivative term $-h'' x^3/6$ in the third-order toroidal potential, $V_T = h x - h(h+k) x^2/2 + [2h(h^2+hk+k^2) - h''] x^3/6$. In the hard-edge limit this term becomes singular in the Hamiltonian, and removing it by a canonical transformation produces exactly the entry/exit jumps above. In the code, GES (a general electrostatic bender element written for the same map-code environment) feeds $h(s)$ and this potential directly into Runge-Kutta integration, and in zero-aperture mode applies precomputed entry and exit matrices so the second-order fringe-field map comes from two matrix applications rather than an integration.","core_discovery":"The central discovery is a two-line hard-edge fringe-field map, not a detailed field integration. At entry, the second-order transfer effect is $\\Delta x = x^2/(2\\rho)$ and $\\Delta P_x = -xP_x/\\rho$; at exit the signs are reversed. The singular term $h'' x^3/6$ in the Frenet-Serret potential expansion is responsible: it must be eliminated by a canonical transformation, and two integrations by parts turn it into exactly these jumps in $x$ and $P_x$. Because the singular coefficient depends only on the bend-plane curvature $h = 1/\\rho$ and not on the vertical curvature $k$, the result holds for spherical, cylindrical, and intermediate electrode geometries. The paper's new GES element, which supplies $h(s)$ and the derived potential directly to the equations of motion, agrees with the theoretical map and with an independent ion-optics program, while the built-in element agrees only when fringe fields are turned off and departs at second order when they are on, because it drops curvature derivatives in the fringe regions.","pith_inferences":["The same canonical-transformation treatment should produce hard-edge fringe maps for magnetic bends or for any element whose curvature falls off longitudinally; deriving the analogous entry/exit matrices is a direct extension.","If the built-in element's omission is systematic, previously published optics designs and emittance-growth estimates that used it for electrostatic bends may need rechecking wherever second-order bend aberrations are important.","The four second-order dispersive coefficients that still disagree between GES and the independent benchmark provide a clean test: run both codes with identical smooth fringe-field shapes and identical aperture to isolate whether the remaining difference is physical or a modeling convention."],"forward_implications":["Any beamline design that omits fringe-field maps carries an irreducible second-order bend aberration that cannot be designed away by shaping element ends.","The zero-aperture GES mode gives the exact second-order fringe-field effect by applying the entry and exit matrices directly, making the correction cheap enough to include in matching and optimization.","The built-in element's second-order maps are missing the curvature-derivative contribution and therefore disagree with GES and the independent calculation whenever fringe fields are enabled.","The hard-edge correction applies to every toroidal electrode geometry, spherical, cylindrical, and intermediate, because only the bend-plane curvature enters.","Adding such fringe-field maps by default would keep transport calculations symplectic and conservation-law-respecting even when the detailed falloff is idealized."],"supporting_citations":[{"why":"Derives the toroidal-coordinate potential expansion and the canonical transformation whose hard-edge limit gives the entry/exit shifts.","marker":"[7]"},{"why":"Establishes the general theorem that the lowest-order irreducible fringe-field effect is independent of falloff shape, justifying the hard-edge map for real electrodes.","marker":"[3]"},{"why":"Presents the earlier cross-validation claiming the built-in element is correct; this paper repeats the case and reaches the opposite conclusion.","marker":"[9]"},{"why":"Supplies an independent ion-optics transfer-map calculation used as the benchmark; GES agrees with it to about one percent.","marker":"[1]"},{"why":"Provides the map-code environment, the coordinate expansion and recursion relation used to build the potential, and the built-in element whose omission is found.","marker":"[2]"}],"fun_headline_variants":["Hard-edge map pins down fringe kick x²/2ρ","Built-in bend element missing curvature derivative kick","Singular term in potential becomes bend fringe kick","Fringe kick holds for spherical, cylindrical electrodes","30-year COSY discrepancy resolved by hard-edge map"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the lowest-order fringe-field effect is independent of the shape of the falloff, so the hard-edge map captures the full irreducible second-order effect for any real electrode geometry.","fun_headline_variants_meta":{"raw":{"variants":["Hard-edge map pins down fringe kick x²/2ρ","Built-in bend element missing curvature derivative kick","Singular term in potential becomes bend fringe kick","Fringe kick holds for spherical, cylindrical electrodes","30-year COSY discrepancy resolved by hard-edge map"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000716,"raw_usage":{"total_tokens":3154,"prompt_tokens":817,"completion_tokens":2337,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":433,"completion_tokens_details":{"reasoning_tokens":2263}},"tokens_in":433,"tokens_out":2337,"duration_ms":24900,"temperature":1.0,"reasoning_tokens":2263,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T12:50:40.321054+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Track a charged particle through a numerically computed or measured electrostatic field of a toroidal bender with a smooth, finite-length fringe region, and extract the entry-map coefficients $(x|xx)$ and $(x|xp)$; if they deviate from $1/(2\\rho)$ and $-1/\\rho$ in a shape-dependent way beyond numerical error, the hard-edge universality claim fails.","supporting_citations":[{"cited_title":"Electrostatic Bender Fields, Optics, Aberrations, with Application to the Proton EDM Ring","cited_arxiv_id":"1508.00157","evidence_quote":"Derives the toroidal-coordinate potential expansion and the canonical transformation whose hard-edge limit gives the entry/exit shifts."},{"cited_title":"Baartman, End Effects of Beam Transport Elements, talk at Snowmass (July 2001)","cited_arxiv_id":null,"evidence_quote":"Establishes the general theorem that the lowest-order irreducible fringe-field effect is independent of falloff shape, justifying the hard-edge map for real electrodes."},{"cited_title":"Valetov, M","cited_arxiv_id":null,"evidence_quote":"Presents the earlier cross-validation claiming the built-in element is correct; this paper repeats the case and reaches the opposite conclusion."},{"cited_title":"Wollnik, J","cited_arxiv_id":null,"evidence_quote":"Supplies an independent ion-optics transfer-map calculation used as the benchmark; GES agrees with it to about one percent."},{"cited_title":"Berz, Computational aspects of optics design and simulation: COSY INFINITY, Nuclear Instruments and Methods in Physics Research 298 (1990) 473–479","cited_arxiv_id":null,"evidence_quote":"Provides the map-code environment, the coordinate expansion and recursion relation used to build the potential, and the built-in element whose omission is found."}],"review_version":1}