{"id":"263eacc7-bd34-4f0a-ad4b-fe4cb7f02cdc","arxiv_id":"2507.06320","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Adding 2HDM-type two-loop diagrams to the NMSSM electron EDM calculation tightens constraints on the CP-violating phase phi4 and changes the allowed parameter space for phi3, phi4, and phi5.","lead":"This paper finds that previously neglected two-loop Feynman diagrams, known from two-Higgs-doublet models, contribute to the electron electric dipole moment predicted by the general NMSSM, a supersymmetric extension of the Standard Model. If correct, these contributions change which NMSSM parameters are allowed by the current electron EDM bound, with implications for whether this model can explain the matter-antimatter asymmetry of the universe.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The electron EDM bounds in Figs. 4-6 are evaluated at CP-violating parameter points that are never validated for vacuum stability or collider exclusions, so the inferred phase limits may not correspond to viable NMSSM configurations.","rationale":"The paper's central quantitative claims are the new constraints on phi3, phi4, phi5 (Fig. 6) and the per-phase sensitivities (Figs. 4-5). These are computed at CP-violating phase values. The text explicitly limits its NMSSMTools validation to the CP-conserving limit, where all CP-violating phases are 0 or pi. This is not a subtle formal issue: the same input parameters with nonzero sin(phi) can yield a different stationary point, a tachyonic scalar, or a vacuum that is not the global minimum. Since the paper solves tadpole equations but does not check global stability or renormalization-group running at those points, the plotted EDM values and the extracted sin(phi4) bound are not yet connected to viable NMSSM vacua. A concrete re-run of the scanned points with a CP-violating spectrum and vacuum checker would settle this. I considered whether the unprinted analytic expressions are a more fundamental risk; however, the substitution logic in Appendix A is plausible and can in principle be checked against Ref. [46], whereas the validation gap is explicitly admitted in the text and directly affects the numerical conclusions. Thus the reader's weakest_assumption identifies the same load-bearing concern, and the conditional verdict remains appropriate: the diagrammatic observation is credible, but the quantitative phase limits are provisional until the CP-violating parameter space is validated.","tokens_in":16895,"tokens_out":13136,"duration_ms":143754,"concrete_test":"Take the benchmark BP and a random subset of the scanned points used in Fig. 6 with sin(phi3,4,5) = +/-0.01, +/-0.05, +/-0.1. For each, solve the full tree-level tadpole equations, compute the scalar mass-squared eigenvalues, and scan the potential along the h_d, h_u, h_s, a, a_s, and charged directions for deeper minima; then run a public spectrum generator that supports complex NMSSM phases (e.g., CPsuperH or SPheno with NMSSM) and NMSSMTools/HiggsBounds-style checks. If any of those points have a tachyonic Higgs, a deeper minimum, or are excluded by LEP/Tevatron/LHC constraints, exclude them and redraw Figs. 4-6. If the allowed region for phi4 changes by more than a factor of about 2, or if most high-|phi4| points are invalid, the paper's quantitative phase limits do not hold for physical NMSSM vacua.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3 states: 'we confirm the viability of our benchmark points with NMSSMTools in the CP-conserving limit, ensuring the absence of Landau poles, false vacua, and conflicts with LEP, Tevatron, and LHC exclusion data.' But the electron EDM constraints in Figs. 4-6 and the scan over phi3, phi4, phi5 are evaluated at CP-violating values where these checks are not repeated. At such points the tadpole equations are solved for m^2_Hd, m^2_Hu, m^2_S, phi7, phi8, but a stationary solution is not sufficient to guarantee a physical vacuum: the tree-level potential may develop tachyonic scalar masses, a deeper charge-breaking or singlet-direction minimum, or Landau poles under renormalization-group running. If a sizeable fraction of the scanned points with sin(phi4) ~ 0.05 fail such checks, then the claimed |sin(phi4)| <= 0.05 bound, and the qualitative statement that phi4 is an order of magnitude more constrained than phi3 or phi5, would be anchored to unphysical parameter points rather than to valid NMSSM configurations. Because the central conclusion is precisely that the new two-loop diagrams alter the allowed electroweak-baryogenesis parameter space, this validation gap directly undermines the quantitative content of the paper, even if the diagrammatic identification itself is correct.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript revisits the electron EDM in the general (non-Z3) NMSSM with explicit CP violation. It constructs a basis of 14 rephasing-invariant CP-violating phases, identifies three classes of two-loop diagrams (charged-Higgs loop, W-boson loop, and Kite diagrams) that were not included in earlier NMSSM EDM calculations, and adapts the analytic 2HDM results of Ref. [46] to the NMSSM via a set of substitution rules. The numerical analysis, based on a benchmark point and a scan over the phases phi3, phi4, phi5, shows that the new diagrams dominate the electron EDM for these phases, tightening the bound on phi4 by roughly an order of magnitude relative to phi3 and phi5, while for phi'_0 and phi6 they partially cancel the previously known MSSM two-loop contributions. The paper concludes that the allowed NMSSM parameter space for electroweak baryogenesis is materially affected by these new contributions.","tokens_in":17151,"tokens_out":10628,"duration_ms":112163,"significance":"If the numerical results are correct, the paper provides a nontrivial correction to the NMSSM electron EDM literature and has direct implications for electroweak baryogenesis model-building. The identification of the missing diagrams is topologically well-motivated: the NMSSM scalar sector has tree-level CP violation, so the 2HDM-type diagrams must be present, and the adaptation from the type-II 2HDM is plausible. The paper also gives a clear basis of invariant phases. The main limitation is that the quantitative constraints are derived at parameter points whose vacuum stability and collider viability are verified only in the CP-conserving limit, so the specific phase limits in Figs. 5 and 6 should be treated with caution until that gap is closed. The paper is not circular: the EDM bound is an external constraint, and no parameter is tuned to reproduce it.","major_comments":[{"comment":"The paper states that 'we confirm the viability of our benchmark points with NMSSMTools in the CP-conserving limit,' but the electron EDM constraints in Figs. 4–6 are evaluated at CP-violating values of the phases (e.g., sin(phi)=0.01 in Fig. 4 and the scan over sin(phi3,4,5) in Fig. 6). No check of Landau poles, false vacua, or collider exclusions is reported for those CP-violating points. Since the central quantitative conclusion is the revised bound |sin(phi4)| ≲ 0.05 and the claim that phi4 is an order of magnitude more constrained than phi3 or phi5, this gap is load-bearing: if a substantial fraction of the scan points have tachyonic scalars or deeper minima, the allowed regions shown in Fig. 6 would not correspond to physical NMSSM vacua. Please repeat the NMSSMTools checks (or an equivalent vacuum-stability and Landau-pole check) at the CP-violating points, at minimum for the benchmark and for a representative sample of the accepted and rejected scan points, and either restrict the scan to viable points or quantify how many points are affected.","section":"Section 3, first paragraph and Figs. 4–6"},{"comment":"The paper's central novelty claim is that the three two-loop diagram classes (d_E^f)^{H±}, (d_E^f)^{W±}, and (d_E^f)^{Kite} 'have been overlooked in the literature.' However, the paper does not explicitly inventory which of these diagrams are present or absent in the most complete previous NMSSM EDM calculations, in particular Ref. [37], which is cited only in the context of cancellations in the electron EDM. To make the novelty claim verifiable, please add a short comparison (either a table or a few sentences) listing the diagram topologies included in Refs. [30,36,37] and any other relevant prior work, and state explicitly which of the three classes are new relative to each.","section":"Introduction and Section 2, with references to Refs. [30,36,37]"}],"minor_comments":[{"comment":"The counting of the tadpole conditions, especially the statement that ∂V/∂a_d and ∂V/∂a_u are degenerate, is presented in a single sentence; a short explanation of how the three imaginary-component tadpole equations reduce to two independent conditions and how φ7 and φ8 are determined would make the basis of 14 phases more transparent.","section":"Section 2.1, paragraph after Eq. (3)"},{"comment":"The phase convention in the numerical study is ambiguous: Eq. (8) sets φ'_0 = φ3 = φ5 = π, but the text says 'setting each phase to a non-zero value of sin(phi)=0.01.' Please state explicitly whether this means φ_i = arcsin(0.01) with cos(φ_i) > 0 or φ_i = π + δ with sin(δ) = 0.01 (with either sign), because the sign of cos(φ_i) affects the Higgs couplings and therefore the EDM.","section":"Section 3, paragraph on Fig. 5 and Eq. (8)"},{"comment":"The derivation of the adapted couplings is compressed: in particular, the definition of the neutral-scalar rotation matrix O and the treatment of the VEV phases e^{iφ_u}, e^{iφ_d}, e^{iφ_s} in the Yukawa couplings are not spelled out. A brief statement of these conventions would make the substitution rule (A.5) checkable and would help readers apply the results in other bases.","section":"Appendix A, Eqs. (A.3)–(A.5)"},{"comment":"The caption says 'The lengths of bars are shown in unit of dExp_e,' but the horizontal axis has values from -10 to 10; please state explicitly that the bar lengths are in units of the experimental bound dExp_e and clarify whether the 'Total' bar is the algebraic sum of the three displayed contributions.","section":"Fig. 4 caption"},{"comment":"The caption says 'The gray region denotes the current experimental bound,' but the plot shows signed values of d_e; please specify that the bound applies to |d_e| and give the numerical value (4.1 × 10^{-30} e cm).","section":"Fig. 5 caption"},{"comment":"There are several typographical artifacts in the text, such as 'o ffers' in the abstract and 'e ffective' in Section 1; please run a spell-check and correct these and similar LaTeX/OCR artifacts.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the core physics—the existence of 2HDM-type two-loop contributions in the NMSSM—is well motivated. The main technical gap is the missing validation of the CP-violating parameter points used for the quantitative EDM constraints; this is fixable by rerunning the public NMSSMTools checks at those points. The novelty claim should also be sharpened by an explicit comparison with the most complete prior calculations, especially Ref. [37]. With these revisions, the paper could become publishable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my read. This paper does a genuine service: previous NMSSM EDM analyses stopped at the MSSM-type two-loop graphs, but this one correctly notes that the NMSSM's scalar sector has tree-level CP violation, so the charged-Higgs, W, and Kite two-loop diagrams from the complex 2HDM must be present. The substitution rules in Appendix A look sound, and the topological case for the diagrams is airtight. The new 14-phase basis is a useful organizing tool, and the phase classification in Table 1 should help future work. The qualitative finding—phi4 gets much tighter after including these diagrams, while phi0' and phi6 see partial cancellations—is well supported by the structure of the couplings.\n\nThe soft spot flagged in the stress test is real. The authors validate benchmark points with NMSSMTools only in the CP-conserving limit, then evaluate the EDM constraints and the Fig. 6 scan at CP-violating phases without rechecking vacuum stability, Landau poles, or collider bounds. That means the headline number, |sin(phi4)| <~ 0.05, is provisional. I don't think this is fatal—the phases involved are small, so many of those points are probably fine—but it is a gap in evidence, and it should be fixable in revision. The paper also doesn't print the full analytic expressions, which makes independent verification harder, and it ships no code or data. Those are addressable issues.\n\nThe central claim holds up. Nothing here is tuned to reproduce the measured EDM; the formulas come from an independent 2HDM calculation. The citation pattern looks fine. I'd send this to a serious referee. The missing phase-space validation and the unpublished analytic details are exactly what a referee should ask for. This is a paper I'd cite when working on NMSSM EDMs, and I'd bring it to reading group as a useful example of transferring 2HDM results into the NMSSM.","headline":"Correct and important observation—NMSSM needs 2HDM-type two-loop EDMs—but the numerical bounds rest on scan points that were never validated for vacuum stability or collider constraints.","tokens_in":17756,"tokens_out":2906,"would_cite":true,"duration_ms":32019,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["11.30.Er","12.60.Jv","13.40.Em"],"model":"deepseek-v4-flash","headline":"The electron's electric dipole moment in the NMSSM receives three overlooked two-loop contributions — charged-Higgs, W-boson, and Kite diagrams — that tighten the bound on the phase φ4 by about an order of magnitude while partly relaxing…","keywords":["electron electric dipole moment","NMSSM","CP violation","two-loop EDM diagrams","Barr-Zee diagrams","charged Higgs","electroweak baryogenesis","two-Higgs-doublet model"],"falsifier":"Recompute the electron EDM over the $\\phi_4$–$\\phi_5$ scan used for Fig. 6 with an independent implementation of the full one- and two-loop diagram set: the paper predicts that including the three new classes excludes points with $|\\sin\\phi_4|\\gtrsim 0.05$ while $\\phi_3$ and $\\phi_5$ remain allowed across the full $\\pm 0.1$ range, and that the exclusion disappears when the new diagrams are dropped. Separately, run the benchmark and scan points through a spectrum generator that accepts complex NMSSM parameters, checking for Landau poles, tachyonic scalars, and deeper minima: if many of the points used for the limits are not valid vacua once CP phases are switched on, the bounds in Figs. 4–6 would not correspond to physical NMSSM configurations.","tokens_in":16633,"feed_emoji":"🧲","tokens_out":21056,"duration_ms":194611,"temperature":0.7,"pith_summary":"The paper argues that the electron electric dipole moment (EDM) predicted by the Next-to-Minimal Supersymmetric Standard Model (NMSSM) receives three classes of two-loop contributions — charged-Higgs loops, $W$-boson loops, and a 'Kite' diagram, all of the type familiar from the two-Higgs-doublet model (2HDM) — that every earlier NMSSM EDM analysis overlooked because those analyses carried over Minimal Supersymmetric Standard Model (MSSM) results without accounting for tree-level CP violation in the extended Higgs sector. For the NMSSM-specific phases $\\phi_3$, $\\phi_4$, and $\\phi_5$ these new diagrams dominate the electron EDM and reverse its sign relative to the previously known terms, tightening the constraint on $\\phi_4$ by roughly an order of magnitude beyond what earlier studies concluded, while for $\\phi_0'$ and $\\phi_6$ the new terms partially cancel the known two-loop contributions and relax the bounds. The paper constructs a basis of 14 independent CP-violating phases and transfers the gauge-invariant 2HDM formulas to the NMSSM by a short set of substitution rules. A sympathetic reader would care because the revised limits reshape the parameter space in which NMSSM electroweak baryogenesis could explain the observed matter–antimatter asymmetry of the universe.","feed_headline":"Three overlooked diagrams tighten NMSSM electron EDM bounds","feed_subtitle":"Including them tightens the constraint on φ4 by an order of magnitude while loosening two others.","key_machinery":"The load-bearing object is the set of three two-loop diagram classes taken over from the two-Higgs-doublet model — the charged-Higgs loop $(d_E^f)^{H^\\pm}$, the $W$-boson loop $(d_E^f)^{W^\\pm}$, and the 'Kite' diagram $(d_E^f)^{\\mathrm{Kite}}$ — each evaluated in the background-field gauge so each group is separately gauge invariant. They carry the argument because they are the only diagrams sensitive to tree-level CP violation in the NMSSM Higgs sector, and they are ported from the complex 2HDM calculation by the substitution $q_{i1}\\to a_i$, $c_f\\,\\mathrm{Im}(q_{i2})\\to g^P_{H_i\\bar f f}$, and $\\lambda_{iH^+H^-}\\to\\lambda_i$, with all phase dependence encoded in the neutral-scalar rotation matrix $O_{ij}$ and the tree-level Higgs couplings of Appendix A. The paper's organizing device is a basis of 14 independent CP-violating invariants — built by assigning spurion charges under the global $U(1)$ symmetries, including a new $U(1)_S$ — together with a table assigning each phase to the one-loop MSSM, two-loop MSSM, or 2HDM diagram groups; this classification is what lets the authors say which previously unconstrained NMSSM phases are actually probed by the electron EDM.","core_discovery":"On the paper's own terms, previous NMSSM fermion-EDM calculations are incomplete: because the NMSSM, unlike the MSSM, has tree-level CP violation in its scalar Higgs sector, the three 2HDM-type two-loop diagram classes $(d_E^f)^{H^\\pm}$, $(d_E^f)^{W^\\pm}$, and $(d_E^f)^{\\mathrm{Kite}}$ must be added to the familiar one-loop and Barr–Zee two-loop sets. Working in the general NMSSM with 14 independent CP-violating phases, the paper adapts the gauge-invariant complex-2HDM electron-EDM results through the substitutions $q_{i1}\\to a_i$, $c_f\\,\\mathrm{Im}(q_{i2})\\to g^P_{H_i\\bar f f}$, and $\\lambda_{iH^+H^-}\\to\\lambda_i$, and finds that for $\\phi_3,\\phi_4,\\phi_5$ the new diagrams dominate, reversing the predicted EDM sign; in a numerical scan $\\phi_4$ is limited to $|\\sin\\phi_4|\\lesssim 0.05$ while $\\phi_3$ and $\\phi_5$ remain allowed across the full scanned range, so satisfying the bound $d_e\\le 4.1\\times 10^{-30}$ e cm demands a fine-tuned cancellation between $\\phi_4$ and $\\phi_5$. For $\\phi_0'$ and $\\phi_6$, conversely, the new contributions are smaller than and opposite in sign to the MSSM two-loop terms, so those phases stay weakly constrained. The upshot is that the electron EDM excludes or permits a qualitatively different set of NMSSM configurations than the standard formulas indicated, which is the paper's reason for revisiting electroweak-baryogenesis studies that relied on the older expressions.","pith_inferences":["If the sign reversal the paper finds for $\\phi_3,\\phi_4,\\phi_5$ holds, then earlier NMSSM electroweak-baryogenesis studies were tested against an incomplete EDM observable; re-deriving the baryogenesis transport equations with the revised CP-violating Higgs couplings could shift the viable parameter region rather than merely shrink it.","The partial cancellation for $\\phi_0'$ and $\\phi_6$ is a concrete accidental-suppression mechanism: an electron-EDM measurement reaching $\\sim 10^{-31}$ e cm would distinguish a genuinely small underlying coupling from the cancellation the paper identifies.","The same diagrams feed quark chromo-EDMs and the Weinberg operator through the heavy internal loops already in the formulas, so a neutron- and atomic-EDM analysis built on these expressions could constrain the NMSSM phases differently than the electron alone.","Because the ordering of the bounds tracks ratios such as $|\\kappa v_s/A_\\lambda|\\sim 0.1$ and $|\\beta/\\kappa v_s|\\sim O(1)$, the pattern is checkable by design: parameter regions with larger $A_\\lambda$ should show a weaker $\\phi_4$ bound, and repeating the scan while varying that coupling would test the prediction."],"forward_implications":["The current electron-EDM bound $d_e\\le 4.1\\times 10^{-30}$ e cm restricts $\\phi_4$ to $|\\sin\\phi_4|\\lesssim 0.05$ in the scan, roughly an order of magnitude tighter than $\\phi_3$ and $\\phi_5$, which remain viable across nearly the full scanned range and evade the bound at the benchmark with $\\sin\\phi\\lesssim 10^{-2}$.","Satisfying the bound requires fine-tuned interplay between $\\phi_4$ and $\\phi_5$; the allowed points in the $\\phi_4$–$\\phi_5$ plane cluster in the first and third quadrants, where the two phases generate opposite-sign EDM contributions that cancel.","For $\\phi_0'$ and $\\phi_6$ the new diagrams partially cancel the previously known two-loop contributions, so these NMSSM-specific phases remain weakly constrained and can take values of $O(10^{-1})$ without violating the bound.","The same three diagram classes are required in any MSSM extension with an enlarged Higgs sector and tree-level CP violation — the paper names the $\\mu\\nu$SSM — and the revised formulas supersede the standard NMSSM expressions used in earlier EDM and baryogenesis studies."],"supporting_citations":[{"why":"Supplies the gauge-invariant two-loop electron-EDM results for the charged-Higgs, W-boson, and Kite diagrams in the complex 2HDM that the paper adapts to the NMSSM via the substitutions of Appendix A.","marker":"[46]"},{"why":"The previous NMSSM fermion-EDM analysis whose diagram set is shown to be incomplete; its analytic expressions and phase conventions are the baseline the paper revisits.","marker":"[30]"},{"why":"The experimental electron-EDM upper bound $d_e\\le 4.1\\times 10^{-30}$ e cm against which every phase constraint in the paper is set.","marker":"[16]"},{"why":"Provides the two-loop MSSM-type EDM expressions (the $\\gamma H$, $WH$, $WW$, $ZH$ diagrams) that make up the MSSM two-loop group in Table 1.","marker":"[35]"},{"why":"Gauge-invariant Barr–Zee contributions to fermion EDMs in the 2HDM, the origin of the charged-Higgs and W-boson upper-loop diagrams.","marker":"[44]"},{"why":"Defines the general (non-$Z_3$-symmetric) NMSSM superpotential and soft-breaking terms from which the CP-violating phase basis is built.","marker":"[28]"},{"why":"Earlier EDM-based study of the CP-violating NMSSM parameter space whose conclusions are altered once the 2HDM-type diagrams are included.","marker":"[36]"},{"why":"The NMSSM spectrum code used to confirm that benchmark points are free of Landau poles, false vacua, and LEP/Tevatron/LHC exclusions in the CP-conserving limit.","marker":"[49, 50]"}],"fun_headline_variants":["Overlooked two-loop diagrams shift NMSSM electron EDM limits","New two-loop diagrams either relax or tighten NMSSM EDM bounds","Missing diagrams flip NMSSM electron EDM constraints","Revisited electron EDM: new NMSSM diagrams matter","Two-loop diagrams change NMSSM electron EDM predictions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the parameter points at which the EDM limits are computed are genuine, stable NMSSM vacua, but the paper's spectrum check validates them only in the CP-conserving limit and at tree level, so the CP-violating points that actually set the phase bounds are not individually verified against false minima, Landau poles, or collider exclusions.","fun_headline_variants_meta":{"raw":{"variants":["Overlooked two-loop diagrams shift NMSSM electron EDM limits","New two-loop diagrams either relax or tighten NMSSM EDM bounds","Missing diagrams flip NMSSM electron EDM constraints","Revisited electron EDM: new NMSSM diagrams matter","Two-loop diagrams change NMSSM electron EDM predictions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000887,"raw_usage":{"total_tokens":3872,"prompt_tokens":1035,"completion_tokens":2837,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":651,"completion_tokens_details":{"reasoning_tokens":2750}},"tokens_in":651,"tokens_out":2837,"duration_ms":21776,"temperature":1.0,"reasoning_tokens":2750,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T19:07:43.240109+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the electron EDM over the $\\phi_4$–$\\phi_5$ scan used for Fig. 6 with an independent implementation of the full one- and two-loop diagram set: the paper predicts that including the three new classes excludes points with $|\\sin\\phi_4|\\gtrsim 0.05$ while $\\phi_3$ and $\\phi_5$ remain allowed across the full $\\pm 0.1$ range, and that the exclusion disappears when the new diagrams are dropped. Separately, run the benchmark and scan points through a spectrum generator that accepts complex NMSSM parameters, checking for Landau poles, tachyonic scalars, and deeper minima: if many of the points used for the limits are not valid vacua once CP phases are switched on, the bounds in Figs. 4–6 would not correspond to physical NMSSM configurations.","supporting_citations":[{"cited_title":"Higgs Mediated EDMs in the Next-to-MSSM: An Application to Electroweak Baryogenesis","cited_arxiv_id":"1102.5679","evidence_quote":"The previous NMSSM fermion-EDM analysis whose diagram set is shown to be incomplete; its analytic expressions and phase conventions are the baseline the paper revisits."}],"review_version":1}