{"id":"60b95403-563b-4dc1-99ff-885e15553832","arxiv_id":"2412.04388","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A generalized, unitarity-bounded parametrization for b-hadron decay form factors that incorporates subthreshold and anomalous branch cuts is constructed.","lead":"This paper extends the standard BGL parametrization of hadronic form factors to handle extra branch cuts, including cuts that run into the complex plane. The new method restores rigorous unitarity bounds, which could tighten theoretical predictions for B meson decays and help interpret anomalies in rare decays.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The anomalous-cut unitarity bound is not derived: the outer function and the Δχ bound are explicitly left out, so the abstract's claim overreaches.","rationale":"The reader's verdict of CONDITIONAL is appropriate. The subthreshold construction for f_+^{BK} and f_0^{BK} is a genuine, largely sound advance: the outer functions are explicit, the Blaschke factors are handled, and the coefficient bound follows if χ~ is treated as an upper bound. A secondary concern is that the subthreshold Δχ is estimated, not rigorously bounded; the paper argues it is <1% of χ_OPE even for K ~ O(100), but a strictly rigorous unitarity bound would require inserting a conservative upper bound (e.g., inflating χ~ by a safety factor) rather than setting χ~ = χ_OPE. This is fixable and not the dominant issue. The dominant issue is the anomalous-cut section, which fails to deliver the central advertised result. The reader's weakest_assumption focused on the unvalidated Schwarz–Christoffel mapping; I agree that the numerical mapping needs validation, but the deeper problem is that the derivation itself is incomplete: the outer function and the Δχ bound are explicitly deferred. Even a mathematically certified conformal mapping would not supply the unitarity bound without these ingredients. The abstract overclaims, and the conditional verdict stands. My proposed concrete test—attempting to complete the missing construction and bound—would determine whether the anomalous-cut claim can be salvaged or must be removed from the headline claim.","tokens_in":10575,"tokens_out":13139,"duration_ms":135169,"concrete_test":"Attempt to complete the anomalous-cut derivation: using the provided Python code for the conformal mapping, compute the image of the anomalous-cut arc on the unit circle and construct a candidate outer function by the standard BGL boundary-measure matching procedure. Then compute a rigorous, conservative upper bound for Δχ along the anomalous arc—not merely an order-of-magnitude estimate. If the outer function has zeros inside the disk, or if Δχ cannot be bounded above by a known quantity with a safety factor, then the claimed unitarity-bounded parametrization for anomalous cuts is not established and the abstract should be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of the paper is that a unitarity-bounded parametrization is derived for form factors with subthreshold and/or anomalous cuts. For subthreshold cuts, the construction is explicit: Eq. (15) defines the parametrization and Eq. (18) states the bound ∑|c_n|^2 < 1. For anomalous cuts, however, the paper stops short of the promised derivation. In the section 'ANOMALOUS BRANCH CUTS', after introducing the Schwarz–Christoffel mapping, the text states: 'The derivation of the outer function is standard and is not shown here. The only point requiring further clarification is the calculation of ΔχJ , which is particularly challenging in this case and beyond the scope of this work.' These are exactly the two elements needed to convert the dispersion inequality into a coefficient bound. Without the outer function, one cannot write the analogue of Eq. (15); without a bound on Δχ (the integral over the anomalous-cut arc), one cannot extend the real-axis spectral integral to the full unit circle. The subthreshold trick of adding a positive term to both sides of the inequality does not automatically transfer: the anomalous cut is not a physical cut of the two-point function on the real axis, so there is no known positive spectral integral to add. The numerical Schwarz–Christoffel mapping, although insufficiently validated, is not the limiting ingredient; even a perfect mapping would not produce the bound without the missing outer function and the missing Δχ estimate. The paper's own text therefore concedes that the abstract's assertion to 'derive unitarity bounds in the presence of ... anomalous branch cuts' is not yet substantiated.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a generalization of the BGL parametrization for b-hadron decay form factors that have subthreshold branch cuts or anomalous (complex) branch cuts. For subthreshold cuts, the authors add the positive subthreshold integral to both sides of the standard dispersion inequality, choose a conformal variable z(q^2, s_Gamma) that removes the cut from the unit disk, and construct an outer function so that the expanded coefficients satisfy sum |c_n|^2 < 1. This is worked out explicitly for f_+^{BK} and f_0^{BK}, including a pole-removal trick for f_0. For anomalous cuts, the paper computes a Schwarz-Christoffel mapping of the cut domain to the unit disk and sketches how the same strategy would apply, but it does not provide the outer function or the required Delta-chi bound. The conclusion claims a systematic procedure for bounded parametrizations with both types of cuts.","tokens_in":10825,"tokens_out":4566,"duration_ms":45324,"significance":"If fully realized, the subthreshold construction would be a useful, rigorous alternative to the standard BGL treatment for form factors such as those in B -> D and B -> K decays, and for non-local form factors where subthreshold effects are not negligible. The f_+^{BK} and f_0^{BK} examples are explicit and the pole-removal for f_0 is a clean, potentially transferable trick. The paper also makes a concrete computational contribution by supplying Python code for the conformal mapping with a single branch cut. However, the anomalous-cut part is presented only as a road map: the essential elements that would convert the dispersion inequality into a coefficient bound are explicitly deferred. As a result, the abstract's claim to derive unitarity bounds for anomalous cuts is not yet substantiated, and the subthreshold bound's rigor depends on an estimate that is plausible but not proven.","major_comments":[{"comment":"The derivation in this section stops short of the advertised result. The text states that 'The derivation of the outer function is standard and is not shown here' and that the calculation of Delta chi_J is 'particularly challenging in this case and beyond the scope of this work.' These are precisely the two ingredients needed to write the analogue of Eq. (15) and to extend the spectral integral to the full unit circle as in Eq. (17). Without the outer function, no explicit expansion of f_nl^{BK} with a unitarity bound can be written; without a bound on Delta chi_J, the inequality for the non-local form factor does not close. The abstract and introduction claim that the paper 'derive[s] unitarity bounds in the presence of ... anomalous branch cuts', which is not supported by the present content.","section":"ANOMALOUS BRANCH CUTS (after Eq. (27))"},{"comment":"The identification tilde-chi^1,bs = chi^1,bs_OPE used in the outer function (16) is not a rigorous upper bound but an estimate. Since Delta chi^1,bs as defined in Eq. (12) is positive, replacing tilde-chi by chi_OPE in phi_+ makes the right-hand side of Eq. (17) equal to tilde-chi/chi_OPE > 1, so the claimed bound sum |c_+,n|^2 < 1 in Eq. (18) only holds if Delta chi is truly negligible. The argument based on Eq. (14) and K < 1 (or even K ~ 100) is a reasonable order-of-magnitude estimate, but it is not a proof; a rigorous first-principles bound requires a rigorous upper bound on |f_+^{BK}| on the interval [s_Gamma, s_+]. The same issue applies to f_0^{BK}, where Eq. (21) leaves Delta chi_0,bs unevaluated and no numerical value for tilde-chi_0 is provided.","section":"SUBTHRESHOLD BRANCH CUTS, Eqs. (12)-(18)"},{"comment":"The numerical Schwarz-Christoffel mapping is not validated. The text describes a quasi-Newton determination of the prevertex z1, but provides no numerical error estimate, no convergence test, and no independent verification that g maps the unit disk to the claimed domain Omega in Eq. (26) (e.g., by checking boundary correspondence or by integrating g' along test arcs). Because the entire anomalous-cut parametrization rests on this mapping and its inverse, the mapping's accuracy must be demonstrated before the construction can be considered reliable.","section":"ANOMALOUS BRANCH CUTS, Eq. (27)"}],"minor_comments":[{"comment":"There is a spacing typo in the title: 'forb-hadron decays' should read 'for b-hadron decays'.","section":"Title and abstract"},{"comment":"The notation z(q^2, s_+) is used in Eq. (1) and in the text, but the definition in Eq. (2) also depends on the free parameter s_0; the s_0 dependence should be made explicit throughout, or stated as implicit.","section":"Eq. (2)"},{"comment":"The polynomials p_n are described as orthonormal on an arc, but only p_0, p_1, p_2 are given in Ref. [9]; the general construction of these polynomials should be stated or referenced clearly so the reader can reproduce the expansion.","section":"Eq. (5)"},{"comment":"The captions refer to 'magenta line' and 'magenta arc' inconsistently; the visual representation of the subthreshold cut and the integration arc should be clarified in the captions themselves.","section":"Figs. 1 and 2"}],"recommendation":"major_revision","confidential_remarks":"The paper has a sound and interesting subthreshold construction, but the anomalous-cut section does not deliver the promised unitarity bound; the abstract and introduction overstate the result. The authors should either complete the anomalous-cut derivation (outer function and Delta chi bound) or retitle/reframe the paper to accurately reflect that only subthreshold cuts are fully treated. The numerical validation of the Schwarz-Christoffel mapping should also be strengthened before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe thing to know: the subthreshold-cut half of this paper is a genuine, largely working advance; the anomalous-cut half is not yet a derivation, and the abstract sells it as one.\n\nWhat's actually new: they construct an outer function ϕ_+(q²) using the mapping z(q²,s_Γ) that restores the clean diagonal unitarity bound ∑|c_n|²<1 for f_+^{BK}. That fixes the known problem with the Gubernari–van Dyk–Virto arc-polynomial parametrization, whose orthonormal polynomials blow up inside the disk. The pole-subtraction trick for f_0, where they remove the B_s0 pole by weighting the dispersion integral with (q²−m_{B_s0}²)², is clean and genuinely useful. The estimate that Δχ is negligible for the physical B→K case (even for K~O(100), sub-percent of χ_OPE) is sensible and kills a potential objection. The explicit f_+ and f_0 outer functions are ready to use.\n\nThe soft spots: the anomalous-cut section stops short. The Schwarz–Christoffel mapping is computed numerically with no accuracy or convergence checks, and the promised code wouldn't be enough anyway. The two ingredients needed to turn the dispersion inequality into a coefficient bound — the outer function and an estimate of Δχ over the anomalous arc — are both explicitly left out. The text says the outer function derivation is 'standard' and the Δχ calculation is 'beyond the scope of this work.' Those are the two load-bearing parts of the method. So the abstract's phrase 'derive unitarity bounds in the presence of ... anomalous branch cuts' overreaches. What they have for anomalous cuts is a plausible strategy plus a numerical mapping. That's a meaningful step toward the result, not the result itself. The reader's stress-test note is right on this.\n\nThe subthreshold claims, by contrast, hold up. The bound (18) follows from (13) the same way BGL's does, and the positivity argument is intact. The paper also correctly flags the limitation of the older Boyd–Grinstein–Lebed subthreshold treatment, which relies on numerically cancelling a branch cut in the disk.\n\nBottom line: cite it for the subthreshold parametrization; treat the anomalous-cut claim as a research program. A serious referee should see it, but the authors should be pushed to either move the anomalous material to an explicit outlook section or actually produce the outer function and the Δχ bound. As it stands the published abstract would mislead.\n\nRecommendation: send to peer review, but expect revision on the anomalous-cut claims.","headline":"Solid new subthreshold parametrization; the anomalous-cut bound is a program, not a derivation, despite what the abstract says.","tokens_in":11406,"tokens_out":2102,"would_cite":true,"duration_ms":20803,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper generalizes the Boyd-Grinstein-Lebed (BGL) parametrization so that unitarity bounds hold for hadronic form factors with subthreshold and anomalous branch cuts.","keywords":["unitarity bounds","form factors","BGL parametrization","subthreshold cuts","anomalous cuts","conformal mapping","B meson decays","dispersion relations"],"falsifier":"Compute the conformal map for the anomalous-cut domain using an independent numerical method and compare its boundary with the claimed slit from $4m_D^2$ to $24.1-3.5i$; any mismatch, or a failure of the induced inequality $\\sum_n |c_n|^2<1$ to hold for a test expansion, would show that the anomalous-cut construction is not valid.","tokens_in":1052,"feed_emoji":"📐","tokens_out":8051,"duration_ms":177212,"temperature":0.7,"pith_summary":"This paper tries to establish that the BGL parametrization, which enforces a unitarity bound on hadronic form-factor expansions, can be extended to form factors whose analytic structure contains subthreshold and anomalous branch cuts. The central claim is that adding the same positive subthreshold integral to both sides of the dispersive inequality restores the standard bound, so the coefficients in the $z(q^2,s_\\Gamma)$ expansion satisfy $\\sum_n |c_n|^2 < 1$ and the truncation error is under control. The paper demonstrates this explicitly for $f_+^{BK}$ and $f_0^{BK}$ and presents a conformal-mapping procedure intended to do the same for non-local form factors with anomalous cuts, such as $f_{nl}$ in $B\\to K\\ell^+\\ell^-$. If the construction is right, it removes a known source of model dependence and allows first-principles unitarity bounds to constrain fits to $B\\to D^{(*)}$, $B\\to K^{(*)}$, and $\\Lambda_b\\to\\Lambda$ decays.","feed_headline":"Unitarity bounds survive subthreshold and complex cuts","feed_subtitle":"Adding the subthreshold cut integral restores rigorous truncation error for B→K, B→D and Λ_b→Λ analyses.","key_machinery":"The load-bearing identity is the modified dispersive inequality obtained by adding the positive subthreshold integral $\\Delta\\chi_J$ to both sides of the OPE bound. This converts an integral over $[s_\\Gamma,\\infty)$ into a contour integral on the unit circle in the variable $z(q^2,s_\\Gamma)$, and with the outer function chosen to absorb both thresholds the expansion coefficients satisfy $\\sum_n |c_n|^2 < 1$. For $f_+^{BK}$ the outer function is given in closed form; for $f_0^{BK}$ the pole between the thresholds is removed by a pole-subtracted combination of subtracted dispersion relations. For the anomalous-cut case, the analogous object is the conformal map $\\hat z$ built from a Schwarz-Christoffel map $g$ and a Möbius transformation, whose construction is the computational core that must deliver the same unitarity-bound inequality.","core_discovery":"The paper claims the first model-independent and unitarity-bounded parametrization for form factors with subthreshold cuts, demonstrated explicitly for $f_+^{BK}$. The key step is to add the subthreshold integral $\\Delta\\chi_{1,bs}(2)$ to both sides of the dispersive inequality; after this addition the full integral from $s_\\Gamma$ to $\\infty$ can be written as an integral around the unit circle in the variable $z(q^2,s_\\Gamma)$, and with a suitably chosen outer function $\\phi_+(q^2)$ that knows about both thresholds, the expansion coefficients satisfy $\\sum_n |c_{+,n}|^2 < 1$. For $f_0^{BK}$, whose pole at $m_{B_{s0}}^2$ lies between $s_\\Gamma$ and $s_+$, the same trick is applied with the pole removed by an additional subtraction, again yielding $\\sum_n |c_{0,n}|^2 < 1$. For the non-local form factor $f_{nl}$ with an anomalous cut between $4m_D^2$ and $s_A=24.1-3.5i$, the paper constructs a Schwarz-Christoffel map of the slit domain to the unit disk and states that the same bounded-parametrization procedure applies, while leaving the explicit numerical bound for $f_{nl}$ and the estimation of the corresponding $\\Delta\\chi_J$ as future work.","pith_inferences":["If the anomalous-cut conformal mapping proves numerically stable, the same strategy could be applied to form factors with several intersecting anomalous cuts, although the paper only treats one cut explicitly.","The pole-removal trick used for $f_0^{BK}$ could be iterated to suppress any subthreshold contribution by additional subtractions, which the paper mentions but does not develop into a general algorithm.","A direct test would be to compute the coefficients $c_n$ for a known toy form factor with a subthreshold cut and verify that the tail respects $\\sum |c_n|^2<1$; the paper does not provide such a numerical demonstration."],"forward_implications":["Analyses of $B\\to D^{(*)}$, $B\\to K^{(*)}$, and $\\Lambda_b\\to\\Lambda$ form factors can now use fits whose truncation error is controlled by $\\sum_n |c_n|^2<1$ while accounting for the real subthreshold cut.","The same procedure gives a template for unitarity-bounded parametrizations of non-local form factors in rare $B$ decays, which previously had no bounded parametrization at all.","The method makes explicit that the subthreshold contribution $\\Delta\\chi_J$ must be estimated or bounded; when the cut is short, as for $B\\to K$, the paper argues this contribution is negligible compared with the OPE uncertainty.","Once an upper bound on $\\Delta\\chi_J$ is supplied, the paper's construction would reduce the model dependence in predictions of $B\\to K\\mu^+\\mu^-$ observables."],"supporting_citations":[{"why":"Supplies the original BGL unitarity-bound framework, including the dispersion inequality and the bound $\\sum_n |a_n|^2<1$ that the paper generalizes.","marker":"[7]"},{"why":"Defines the non-local form factors and proposes the arc-polynomial parametrization that the paper argues is not unitarity-bounded.","marker":"[9]"},{"why":"Locates the anomalous branch cut endpoints $s_\\Gamma=4m_D^2$ and $s_A=24.1-3.5i$ for the $B\\to K\\ell^+\\ell^-$ non-local form factor.","marker":"[12]"},{"why":"Gives the subthreshold threshold $s_\\Gamma=(m_{B_s}+m_\\pi)^2$ and the $\\bar{B}K$ contribution used in the dispersive inequality for $f_+^{BK}$.","marker":"[17]"},{"why":"Describes the Schwarz-Christoffel parameter determination method used to construct the anomalous-cut conformal map.","marker":"[39]"}],"fun_headline_variants":["Subthreshold cuts tamed for b-hadron form factors","New unitarity bounds with extra branch cuts","BGL parametrization extended to subthreshold cuts","First-principles form factors with anomalous cuts"],"cache_read_input_tokens":13440,"weakest_assumption_plain":"For the anomalous-cut part, the entire bound rests on the assumption that the numerically computed Schwarz-Christoffel mapping really does map the unit disk onto the claimed slit domain, and that the associated outer function makes the unitarity inequality hold; the paper gives no numerical error estimate or independent check.","fun_headline_variants_meta":{"raw":{"variants":["Subthreshold cuts tamed for b-hadron form factors","New unitarity bounds with extra branch cuts","BGL parametrization extended to subthreshold cuts","First-principles form factors with anomalous cuts"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000647,"raw_usage":{"total_tokens":2995,"prompt_tokens":992,"completion_tokens":2003,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":608,"completion_tokens_details":{"reasoning_tokens":1951}},"tokens_in":608,"tokens_out":2003,"duration_ms":15466,"temperature":1.0,"reasoning_tokens":1951,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T21:24:09.269574+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the conformal map for the anomalous-cut domain using an independent numerical method and compare its boundary with the claimed slit from $4m_D^2$ to $24.1-3.5i$; any mismatch, or a failure of the induced inequality $\\sum_n |c_n|^2<1$ to hold for a test expansion, would show that the anomalous-cut construction is not valid.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Describes the Schwarz-Christoffel parameter determination method used to construct the anomalous-cut conformal map."}],"review_version":1}