{"id":"ef4920b6-0971-4619-9d0d-2f06184ab0b0","arxiv_id":"hep-ph/0404220","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The lightest CP-odd Higgs of the NMSSM can be as light as about 10 MeV without conflicting with rare K, B, and Upsilon decay data, and its b-physics signatures differ qualitatively from the MSSM.","lead":"A one-loop calculation of how the lightest pseudoscalar Higgs of the next-to-minimal supersymmetric Standard Model would appear in rare K, B, and Upsilon decays, finding that masses as low as about 10 MeV are not excluded by current data. It matters because it identifies concrete search channels and two qualitative ways the NMSSM differs from the MSSM in flavor physics.","discovery_kind":"unclear","skeptic_critique":{"model":"selfhosted/deepseek-v4.1-flash","headline":"The light-A_1^0 contribution to B-Bbar mixing (Eq. 26) grows as 1/m_{A_1^0}^2 and is most restrictive in the O(10) MeV window the paper claims is allowed; Figure 2 shows the Delta m_d constraint only for m_A >= 500 GeV, leaving the small-m_{A_1^0} region untested.","rationale":"The paper's qualitative MSSM-vs-NMSSM differences (RG mixing to dipoles, breakdown of the Bs -> mu+mu- / Delta m_s correlation) are durable and argued cleanly. The O(10) MeV claim is a scan result, and the reader correctly notes the scan priors and soft exclusions. My concern adds a specific quantitative gap: the 1/m_{A_1^0}^2 growth of the light-A_1^0 Delta m_d contribution is not displayed for small m_{A_1^0}. This is checkable and would settle whether the O(10) MeV window survives. The reader's lifetime concern is related but distinct; both point to the same conclusion that the O(10) MeV number should be read as 'not excluded by this scan' rather than established. The verdict remains CONDITIONAL.","tokens_in":17298,"tokens_out":20182,"duration_ms":342677,"concrete_test":"Recompute f_d from Eqs. (26)-(28) at m_{A_1^0} = 10 MeV for the parameter points in Figure 2 with |delta_- v/x| ~ 0.1-0.5 and m_A = 500 GeV, using the same C_A from Eq. (9). If f_d < -0.6 for all such points, the O(10) MeV window is excluded by Delta m_d and the headline claim weakens.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's Eq. (26) gives C_SRR1 = -(1/4pi^2) C_A^2 m_W^2/m_{A_1^0}^2. Since C_A (Eq. 9) is independent of m_{A_1^0}, this contribution grows without bound as m_{A_1^0} -> 0. The constraint f_d > -0.6 (Eq. 28) therefore becomes most restrictive precisely at small m_{A_1^0}. Yet Figure 2 displays the Delta m_d constraint only for m_A >= 500 GeV and m_A >= 1000 GeV, and Section 4.2 says 'We find constraints for larger values of |delta_- v/x| and m_A >= 500 GeV.' The O(10) MeV region sits at small |delta_- v/x| and small m_{A_1^0}; the paper does not show whether the 1/m_{A_1^0}^2 growth of Eq. (26) is compensated by a sufficiently small C_A there. If it is not, the O(10) MeV window is excluded by Delta m_d, not allowed. This is a specific, checkable gap in the scan, distinct from the reader's lifetime concern.","agreement_with_reader":"partial"},"referee_report":{"model":"selfhosted/deepseek-v4.1-flash","summary":"The manuscript computes the b → s A_1^0 amplitude at leading order in tan β in the NMSSM, where A_1^0 is the lightest CP-odd Higgs, and uses it to derive constraints from rare K, B, and radiative Υ decays, B − B̄ mixing, and B_s → μ+μ−. The central claim (abstract:s4) is that A_1^0 masses down to O(10) MeV survive all these constraints, opening a low-mass window and identifying b → s τ+τ− as a search channel with branching ratios around 10^-3. The paper also argues that the NMSSM breaks the MSSM correlation between B_s → μ+μ− and B_s − B̄_s mixing, and that neutral-Higgs renormalization of the photon and gluon dipole operators can be sizeable in the NMSSM but not in the MSSM.","tokens_in":35050,"tokens_out":44266,"duration_ms":790243,"significance":"The paper is a careful, comprehensive phenomenology study. The one-loop matching calculation is standard and the loop functions are given explicitly in Appendix B, so the amplitude is checkable. The identification of b → s τ+τ− as a discovery channel and the demonstration that the MSSM B_s → μ+μ− / Δm_s correlation is not a generic SUSY prediction are genuine, useful results. The paper is a solid contribution to the NMSSM phenomenology literature.","major_comments":[{"comment":"The A_1^0 contribution to B − B̄ mixing, C_SRR1 = −(1/4π²) C_A² m_W²/m_A1², scales as 1/m_A1² and is therefore most restrictive at small m_A1, precisely the region where the paper claims O(10) MeV masses are allowed. Figure 2 displays the Δm_d constraint only for m_A ≥ 500 GeV and m_A ≥ 1000 GeV, and Section 4.2 states only that constraints are found 'for larger values of |δ− v/x| and m_A ≥ 500 GeV.' The manuscript does not show whether the 1/m_A1² growth is compensated by a sufficiently small C_A in the small-m_A1 region, nor does it explain why the constraint is truncated at m_A ≥ 500 GeV. Since the headline claim (abstract:s4) depends on the O(10) MeV window surviving, the paper should either display the Δm_d constraint across the full m_A1 range or give an explicit argument (e.g., that at small m_A1 the constraint lies inside the rare-K excluded region of Eq. (21), or that C_A is GIM-suppressed there) showing that it does not close the window.","section":"Section 4.2, Eq. (26), Eq. (28), Figure 2","authors_acknowledge":null,"authors_acknowledge_p":null}],"minor_comments":[{"comment":"The caption and text should state explicitly what the Δm_d curves represent (e.g., the f_d = −0.6 contour) and why they are truncated at m_A ≥ 500 GeV, so that the reader can judge whether the small-m_A1 region is constrained or merely not plotted.","section":"Figure 2","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The statement that 'all allowed points are at the GIM boundary θ_t ≃ π/2, which is set by our value of the cut-off ε' suggests that the allowed region in this mass range is an artifact of the scan's GIM cut. The paper should quantify the sensitivity of the allowed region to ε.","section":"Section 4.1.2","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The bound |δ− v/x| < 3.7 for m_A1 < 23 MeV/|δ− v/x| < 2m_μ is self-referential, with |δ− v/x| appearing on both sides. The paper should state the self-consistent solution explicitly.","section":"Section 4.4, Eq. (35)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The MSSM comparison uses FeynHiggs v2.02, but the scan range and number of points are not given. Stating these would let the reader judge the robustness of the quoted bounds |C_S + C_P|/|C_S − C_P| < 0.2 and |C^b_R| < 1.3 × 10^-3.","section":"Section 5.1","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The statement that masses in the range of Eq. (20) are 'disfavored since the B → K A_1^0 decay would happen too rapidly for most of the parameter space, although cannot rigorously be excluded' is vague. The paper should state what fraction of the scanned parameter space is excluded and under what assumptions.","section":"Section 4.1.1","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"major_revision"},"author_rebuttal":null,"desk_editor":{"model":"selfhosted/deepseek-v4.1-flash","letter":"Hiller's paper is a solid, useful contribution to NMSSM flavor phenomenology, and the durable parts are not the headline number. What's actually new: the b → s A_1^0 amplitude at leading order in tan β in the NMSSM (Eqs. 8-11). The MSSM analogue was known; the NMSSM singlet admixture changes the coupling structure, and Hiller works it out cleanly with explicit loop functions in Appendix B. The second durable point is qualitative: the MSSM relation Eq. (1) fails in the NMSSM, so neutral Higgs contributions to the photon and gluon dipole operators don't cancel, and the Bs → μμ / Δm_s correlation breaks down. Both observations are argued well and are worth citing. The headline — A_1^0 masses down to O(10) MeV are allowed — is a scan result, not a theorem. It depends on priors for stop and scharm masses (~1 TeV), a GIM cutoff ε = 0.05, and assumed MSSM ranges for Δm_d. The paper is honest that all FCNC bounds can be evaded by tuning the super-GIM, but the flip side is that the allowed region is partly an artifact of where you put the cutoff. Read the O(10) MeV number as 'not excluded by this scan' rather than 'established.' The stress-test note flags Eq. (26): the B-Bbar mixing contribution grows as 1/m_{A_1^0}², so Δm_d should be most restrictive at small m_{A_1^0}. The scaling is real, but the specific claim that Figure 2 leaves the small-mass region untested doesn't hold up — the Δm_d curves span the full m_{A_1^0} range; the 'm_A ≥ 500 GeV' label refers to the heavy Higgs mass, not the light one. And at the δ_- values required by the K → π and Υ bounds (δ_- v/x ~ 0.5-2 at m_{A_1^0} = 10 MeV), C_A is small enough that Δm_d doesn't bite. Still, Section 4.2's statement that constraints appear 'for larger values of |δ_- v/x|' is loose: C_A has a δ_- independent piece (the λ term in Eq. 9), so the constraint doesn't vanish as δ_- → 0. The text should say this more clearly. Smaller things: the B → K missing-energy bound is shaded as excluded in Figure 2 but called 'disfavored' in the text — minor inconsistency. The allowed points in the 2m_μ to 2m_τ window sit at the GIM boundary set by ε, which is a scan artifact. And the matching contributions to C_7,8 from neutral Higgs loops are left for future work, so the RG mixing analysis is an estimate, not a full calculation. Who should read it: anyone working on NMSSM flavor or light pseudoscalar searches. The amplitude and the qualitative MSSM-vs-NMSSM differences are worth citing. The exclusion plot is useful but should be read with the priors in mind. I'd send it to referees — the calculation is checkable and the qualitative points are durable.","headline":"A clean NMSSM calculation with durable qualitative points, but the O(10) MeV headline is a scan result that should be read as 'not excluded' rather than 'established.'","tokens_in":18268,"tokens_out":10946,"duration_ms":192200},"pith_extraction":{"msc":[],"pacs":[],"model":"selfhosted/deepseek-v4.1-flash","headline":"In the next-to-minimal supersymmetric Standard Model, the lighter CP-odd Higgs can be as light as about 10 MeV without conflicting with rare kaon, B-meson, or Upsilon data, because a softly broken global symmetry keeps its mass small.","keywords":[],"falsifier":"Placeholder.","tokens_in":18547,"feed_emoji":"","tokens_out":2021,"duration_ms":36948,"temperature":0.7,"pith_summary":"Supersymmetry's minimal version ties the CP-odd Higgs mass to the electroweak scale, but the next-to-minimal version adds a gauge-singlet superfield whose Higgs potential carries an approximate global $U(1)_R$ symmetry, broken only by small trilinear soft terms, so one pseudoscalar $A_1^0$ can naturally sit far below the weak scale. Hiller asks whether existing $b$-physics data already exclude such a state. She computes the one-loop $b$ to $s A_1^0$ amplitude at leading order in $1/|1 - 2kx/A_1^0|$... more precisely at leading order in $1/|1 - 2kx/A_1^0|$","feed_headline":"","feed_subtitle":"","key_machinery":"Placeholder.","core_discovery":"Placeholder.","pith_inferences":[],"forward_implications":[],"supporting_citations":[],"fun_headline_variants":["The lightest NMSSM pseudoscalar may weigh only about 10 MeV","b-physics allows an NMSSM pseudoscalar down to 10 MeV","Rare K and B decays constrain the NMSSM light pseudoscalar","b → s τ⁺τ⁻ searches can probe the NMSSM pseudoscalar","NMSSM and MSSM differ in B_s mixing and B_s → μ⁺μ⁻"],"cache_read_input_tokens":0,"weakest_assumption_plain":"Placeholder.","fun_headline_variants_meta":{"raw":{"variants":["The lightest NMSSM pseudoscalar may weigh only about 10 MeV","b-physics allows an NMSSM pseudoscalar down to 10 MeV","Rare K and B decays constrain the NMSSM light pseudoscalar","b → s τ⁺τ⁻ searches can probe the NMSSM pseudoscalar","NMSSM and MSSM differ in B_s mixing and B_s → μ⁺μ⁻"]},"model":"selfhosted/deepseek-v4.1-flash","effort":"low","attempts":1,"cost_usd":0.0,"raw_usage":{"total_tokens":1590,"prompt_tokens":789,"completion_tokens":801,"prompt_tokens_details":null,"completion_tokens_details":{"reasoning_tokens":699}},"tokens_in":789,"tokens_out":801,"duration_ms":12844,"image_pages":[],"temperature":1.0,"reasoning_budget":50,"reasoning_tokens":699,"transport_provider":"selfhosted","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-11-13T08:46:05.547542+00:00","model_set":{"reader":"selfhosted/deepseek-v4.1-flash"},"pipeline_version":"pith-pipeline@v1.5.0-archive-selfhost","input_kind":null,"input_chars":null,"input_total_chars":null,"passes_missing":null,"published_standing":null,"published_confidence":null,"objection_audit":null,"run_record":null,"basis":{"kind":"not_recorded","source":null,"chars_given":null,"chars_total":null},"receipt":null,"falsifier":"Placeholder.","supporting_citations":[],"lane_accuracy":{"status":"not yet measured","lane":{"model":"selfhosted/deepseek-v4.1-flash","input_kind":null},"url":"https://pith.science/accuracy"},"review_version":1}