{"id":"0452c312-48f4-4876-ba85-500d137feb6e","arxiv_id":"2501.02217","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"In the TNMSSM, the new CPV couplings χ_d and χ_t give significant contributions to EDMs, with neutron EDM data constraining their phases and b, c quark EDMs reaching 10^-22 and 10^-23 e·cm.","lead":"This paper calculates electric dipole moments for the neutron, electron, mercury, and b and c quarks in a supersymmetric model with extra triplet and singlet fields, focusing on two new CP-violating couplings. It finds that the neutron EDM bound strongly restricts those couplings, while the predicted b and c quark EDMs could be large enough to interest future experiments.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Neutron EDM hadronic relation (Eq. 2) is the bottleneck: coefficient fixed points are used to derive quark EDM bounds, and the paper never propagates the stated ±1 and ±10 MeV uncertainties, so the claimed constraints on χ_d, χ_t and the d_b, d_c ceiling are not robustly established.","rationale":"The reader identified essentially the same weakness: the paper fixes the hadronic relation's coefficients at 0.5 and 12 MeV without propagating the stated ±1 and ±10 uncertainties. My independent read of the text confirms this is the most load-bearing step, because the entire quantitative story (large χ_d/χ_t contributions, strict constraints, and observable d_b/d_c) is calibrated through Eq. (2). The loop amplitudes and model setup are standard and internally coherent; no formal verification or reproducible code is provided, but that is a reproducibility issue rather than an internal contradiction. The reader's additional concern about an internal four-order-of-magnitude inconsistency in the d_b/d_c predictions is real: the abstract and summary quote 10^-22 and 10^-23 e·cm, while Section IV states 10^-26 and 10^-27 e·cm for the χ_d/χ_t-only contributions. That is an internal inconsistency, but it concerns presentation of which parameter region produces which prediction; the hadronic uncertainty concern is more fundamental because it undermines the quantitative calibration of the constraints themselves. I therefore agree with the reader's identified weakest assumption, and retain the CONDITIONAL verdict: the claim is plausible and the framework is reasonable, but the specific numbers are not yet trustworthy without propagating the hadronic uncertainties and resolving the internal prediction discrepancy.","tokens_in":20503,"tokens_out":1953,"duration_ms":16323,"concrete_test":"Recompute the χ_d and χ_t exclusion contours and the maximum d_b/d_c values using the full uncertainty ranges in Eq. (2), i.e. take the coefficient of the first bracket as 1±0.5 at both endpoints, vary the C5 prefactor from 12 to 32 MeV, and propagate these through the RG running to the inferred |dg_c| and |dg_b| limits. If the allowed χ_d/χ_t regions and the reported d_b/d_c maxima shift by more than a factor of two, the central quantitative claims need revision.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central quantitative claim is that the new TNMSSM CPV sources χ_d and χ_t are tightly constrained by the neutron EDM upper bound and that d_b and d_c can still reach 10^-22 e·cm and 10^-23 e·cm. The chain of reasoning passes through Eq. (2), d_n = (1±0.5)[1.4(dγ_d − 0.25 dγ_u) + 1.1 e(dg_d + 0.5 dg_u)] ± (22 ± 10) MeV C5, which the authors explicitly evaluate at the fixed coefficient values 0.5 and 12 MeV (Section III.B). They then use the resulting bound to translate the d_n limit into limits on d_b and d_c, citing Ref. [22]. Two problems make this the most load-bearing step. First, the hadronic uncertainties on the two terms and on C5 are order-100% and are dropped without discussion; a factor-of-√[(0.5)²+(0.5)²] or larger shift in the d_n prediction directly changes the allowed regions of χ_d and χ_t and the inferred upper bounds on d_b and d_c by a comparable factor. Second, the paper does not show the intermediate mapping from d_n to |dg_c| and |dg_b|, so the reader cannot check whether the claimed bounds and the quoted figures actually follow. The loop calculation itself is conventional and internally consistent in structure, so the weakest point is not the model amplitude but the hadronic bridge from a quark-level dipole to the neutron observable.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper computes electric dipole moments of the neutron, electron, mercury, b quark, and c quark in the TNMSSM, including one-loop diagrams, two-loop gluino corrections, Barr-Zee-type diagrams, and the Weinberg operator. It studies the dependence of these EDMs on the gluino phase θ3, the μ phase θμ, and the TNMSSM-specific CPV sources χd and χt, and it compares the results with MSSM predictions. The central claims are that the neutron EDM upper bound imposes strict constraints on χd and χt and that the b and c quark EDMs can reach about 10^-22 e·cm and 10^-23 e·cm respectively, making them potentially observable in future experiments.","tokens_in":20899,"tokens_out":6692,"duration_ms":63128,"significance":"If the numerical results are correct, the paper provides a useful phenomenological study of CPV in a well-motivated SUSY extension and gives concrete, falsifiable predictions for b and c quark EDMs. The analytic expressions and the explicit TNMSSM-versus-MSSM comparison are strengths, and I see no circularity: the EDM predictions are computed from the model Lagrangian and then compared with experimental bounds. However, the headline numbers contain a four-order-of-magnitude internal inconsistency, and the hadronic uncertainty in the neutron-EDM relation is not propagated; both issues directly affect the claimed constraints and the quoted reach of d_b and d_c.","major_comments":[{"comment":"The text after Fig. 7 states that d_b and d_c can reach 10^-26 e·cm and 10^-27 e·cm respectively from the new CPV sources alone, but Fig. 8(c)-(f) has axes labeled db×10^22 and dc×10^23 with order-unity values, and both the Abstract and Sec. V quote 10^-22 e·cm and 10^-23 e·cm. This is a four-order-of-magnitude discrepancy in the central quantitative claim of the paper. Please correct the statement or clearly distinguish the scenarios that produce the two sets of numbers.","section":"Sec. IV (after Fig. 7); Abstract; Sec. V"},{"comment":"Equation (2) contains the hadronic uncertainties (1±0.5) and (22±10) MeV, but the text states that the coefficients are fixed to 0.5 and 12 MeV, and the numerical analysis never propagates these uncertainties. Because the paper's headline result—that |d_n| imposes strict constraints on χ_d and χ_t, and hence sets the allowed d_b and d_c ceiling—passes through this relation, the order-100% uncertainty on the hadronic coefficients should be reflected in the constraints. Please provide the resulting uncertainty band or justify why fixing the central values is sufficient.","section":"Sec. III.B, Eq. (2)"},{"comment":"The paper derives or adopts bounds on d_b and d_c from a mapping through |dg_c| and dg_b, citing Ref. [22], but it does not show the intermediate steps connecting the neutron bound to these quark-level bounds. Since the claimed reach of d_b ~ 10^-22 e·cm and d_c ~ 10^-23 e·cm depends on this chain, please display the explicit mapping or provide the numerical inputs so the reader can verify the quoted limits.","section":"Sec. I, Eqs. (1)-(3); Sec. IV"}],"minor_comments":[{"comment":"The text says 'the fact can be seed explicitly'; this should read 'seen'. Similar typos appear elsewhere, such as 'in in Fig. 6'.","section":"Sec. IV"},{"comment":"The caption of Fig. 7 does not specify the units of the plotted quantities; please add the axis labels and units so that the numbers quoted in the text can be checked against the figures.","section":"Fig. 7 caption"},{"comment":"Ref. [64] is incomplete ('B. Yan, S. M. Zhao, and T. F. Feng et al. [Authors]'), and Ref. [76] contains a stray bibliographic entry ('975 (2022) 115671, doi: ...').","section":"References"},{"comment":"The mercury EDM formula involves dg_s, but no strange-quark CEDM expression is given in Sec. III; please state how dg_s is obtained.","section":"Eq. (27)"},{"comment":"The argument of ℑ and the H function in Eq. (10) are not cleanly delimited by parentheses; please fix the notation for readability.","section":"Eq. (10)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a conventional SUSY EDM calculation, and the main obstacle is the internal inconsistency in the quoted b and c quark EDM values. Once that is corrected and the hadronic uncertainty issue is addressed, the paper could be suitable for publication. I do not see evidence of circularity or of deliberate overclaiming, but the headline numbers need to be reconciled."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a workmanlike TNMSSM EDM calculation. The genuinely new piece is the focus on the TNMSSM-specific CPV phases χ_d and χ_t and their effect on d_b and d_c, plus the claim that the neutron EDM bound constrains those phases tightly. That part is plausible and worth a careful look.\n\nWhat the paper does well: the effective-operator formalism is standard but competently applied, with one-loop, two-loop, and Weinberg-operator contributions all included. The observation that χ_t dominates over χ_d because of tan β enhancement is sensible, and the comparison to the MSSM in the same parameter space is a nice way to isolate the new effects. The prior literature, including the authors' own TNMSSM EDM papers, is cited appropriately.\n\nNow the soft spots. First and most serious: the paper states two different predictions for d_b and d_c. In Section IV, after the Fig. 7 scan, it says d_b and d_c can reach 10^-26 and 10^-27 e·cm from the new CPV sources alone. But the Fig. 6 axes, the Fig. 8 axes, and the Summary all quote 10^-22 and 10^-23 e·cm. That is four orders of magnitude. It cannot be a minor typo in one place because the figures are labeled with the larger values and the Summary repeats them; if the 10^-26 statement is the real one, then the \"observable in future\" claim collapses. Resolving this is essential.\n\nSecond, the hadronic relation in Eq. (2) is the load-bearing bridge from the neutron EDM to the quark-level constraints, and the paper fixes the coefficients at 0.5 and 12 MeV without propagating the ±1 and ±10 MeV uncertainties. The stress-test note is right: the allowed regions for χ_d and χ_t, and therefore the inferred ceilings on d_b and d_c, shift by order-one factors under those hadronic uncertainties. Even a short sensitivity check would help.\n\nThird, there is no code or numerical table, and the mapping from d_n to |dg_c| and |dg_b| is asserted via Ref. [22] rather than shown. A referee can't verify the chain without that intermediate step.\n\nThe loop calculation itself looks structurally sound; I don't see a fatal error in the amplitudes. But as submitted, the quantitative claims are not internally consistent, and the hadronic uncertainty issue makes the bounds less robust than the text implies. This paper deserves peer review — a referee can sort out the inconsistency — but it should be a conditional accept at best. I wouldn't cite it until the numbers are settled.","headline":"A workmanlike TNMSSM EDM calculation with a four-order-of-magnitude internal inconsistency between the quoted d_b and d_c predictions that needs to be resolved before the headline numbers can be trusted.","tokens_in":21470,"tokens_out":2598,"would_cite":false,"duration_ms":24244,"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":"The TNMSSM's new CP-violating couplings drive observable EDMs, with the neutron bound the tightest constraint.","keywords":["electric dipole moment","CP violation","TNMSSM","supersymmetry","neutron EDM","b quark EDM","c quark EDM","Weinberg operator"],"falsifier":"Recompute the bounds on $\\chi_d$ and $\\chi_t$ and the maxima of $d_b$ and $d_c$ using the full $\\pm0.5$ and $\\pm10\\,\\text{MeV}$ uncertainty band of Eq. (2) instead of the fixed central values the paper adopts; the claimed $10^{-22}\\,\\text{e}\\!\\cdot\\!\\text{cm}$ reach for $d_b$ survives only if it changes by less than an order of magnitude under that propagation.","tokens_in":20286,"feed_emoji":"⚛️","tokens_out":9217,"duration_ms":86659,"temperature":0.7,"pith_summary":"This paper tries to establish that the TNMSSM—a supersymmetric model extended by two scalar triplets and a singlet—contains new CP-violating couplings, $\\chi_d$ and $\\chi_t$, that generate sizable electric dipole moments for the neutron, electron, mercury, and the $b$ and $c$ quarks. The authors compute these EDMs from one-loop diagrams, two-loop gluino and Barr-Zee diagrams, and the Weinberg operator, then confront them with current experimental upper bounds. They find that the neutron EDM bound is the sharpest constraint on $\\chi_d$ and $\\chi_t$, while the $b$- and $c$-quark EDMs can still reach about $10^{-22}\\,\\text{e}\\!\\cdot\\!\\text{cm}$ and $10^{-23}\\,\\text{e}\\!\\cdot\\!\\text{cm}$, respectively, which future experiments might observe. If correct, EDM experiments can directly probe the CP-violating sector of this model and help decide whether it can explain the matter-antimatter asymmetry.","feed_headline":"b-quark EDM could reach 10^-22 e·cm under neutron-EDM bounds","feed_subtitle":"A triplet-singlet SUSY extension keeps heavy-quark EDMs observable while the neutron bound squeezes the new CP phases.","key_machinery":"The analysis is organized by a five-operator effective Lagrangian: the quark EDM operator, the quark CEDM operator, and the purely gluonic Weinberg operator, with renormalization-group running from the matching scale down to the quark mass scale. The new couplings $\\chi_d$ and $\\chi_t$ enter the superpotential through $\\chi_d H_d \\cdot T H_d$ and $\\chi_u H_u \\cdot \\bar{T} H_u$, shift the squark and neutralino/chargino mass matrices, and supply the CP-violating phases that feed the loop amplitudes. Two-loop Barr-Zee diagrams, two-loop gluino self-energy diagrams, and the Weinberg operator are all needed to match the measured observables, and the neutron and mercury EDMs are then assembled from quark EDMs and CEDMs through hadronic relations.","core_discovery":"The central claim is that the TNMSSM's own CP-violating couplings, $\\chi_d$ and $\\chi_t$, produce electric and chromoelectric dipole moments that are not merely present but experimentally relevant. Through one-loop diagrams, two-loop gluino and Barr-Zee diagrams, and the Weinberg operator, the authors compute $d_n$, $d_e$, $d_{\\text{Hg}}$, $d_b$, and $d_c$ and compare them with current bounds. They find that $\\chi_t$, enhanced by $\\tan\\beta$ and acting through up-type squarks, dominates the EDMs; the neutron upper bound then sharply restricts $\\chi_d$ and $\\chi_t$, while $d_b$ and $d_c$ can still be as large as about $10^{-22}\\,\\text{e}\\!\\cdot\\!\\text{cm}$ and $10^{-23}\\,\\text{e}\\!\\cdot\\!\\text{cm}$, respectively. In the same parameter region the TNMSSM predictions are larger than the MSSM ones, making the model distinguishable through EDM measurements.","pith_inferences":["The paper stops at quoting reachable values; a natural next step is to map the same parameter space onto CP asymmetries in $b \\to s\\gamma$ or $B$-meson decays, which would give an independent handle on $d_b$ without a dedicated EDM measurement.","If a future EDM experiment sees a signal, the correlation structure among $d_n$, $d_e$, $d_{\\text{Hg}}$, $d_b$, and $d_c$—all driven by $\\chi_t$ and the $\\mu$ phase—would help disentangle the TNMSSM phases from MSSM phases, something the paper does not explicitly work out.","The large $\\mu$ phase preferred for electroweak baryogenesis is the same phase most tightly bounded by $d_e$ and $d_{\\text{Hg}}$; the paper notes cancellations can hide it, but does not quantify how much fine-tuning that requires, leaving the model's baryogenesis viability as an open question."],"forward_implications":["The neutron EDM bound excludes most of the $\\chi_t$ parameter space, leaving only $|\\chi_t| \\lesssim 0.2$ when other phases are fixed, so neutron EDM experiments act as a direct filter on the TNMSSM's new CPV sources.","Because $\\chi_t$ is enhanced by $\\tan\\beta$ and enters through up-type squark and chargino/neutralino loops, the TNMSSM predictions for all five EDMs exceed the MSSM predictions in the same parameter region.","Within current constraints, $d_b$ reaches about $10^{-22}\\,\\text{e}\\!\\cdot\\!\\text{cm}$ and $d_c$ reaches about $10^{-23}\\,\\text{e}\\!\\cdot\\!\\text{cm}$, values that future heavy-quark EDM searches could probe.","The electron and mercury EDMs constrain the phase of $\\mu$ more strongly than the neutron EDM does, so combining all four observables gives a tighter test of the model than any single bound."],"supporting_citations":[{"why":"compiles the experimental upper bounds on $d_n$, $d_e$, $d_{\\text{Hg}}$, $d_b$, and $d_c$ quoted in Eq. (1).","marker":"[16–22]"},{"why":"provides the improved upper bounds on $b$- and $c$-quark EDMs and the hadronic running used to convert neutron limits into heavy-quark constraints.","marker":"[22]"},{"why":"supplies the neutron-EDM relation in Eq. (2) connecting quark EDMs, CEDMs, and the Weinberg operator.","marker":"[23]"},{"why":"provides the mercury-EDM formula $d_{\\text{Hg}}$ from quark CEDMs used in the numerical analysis.","marker":"[65]"},{"why":"gives the two-loop Barr-Zee type contributions to the quark EDM.","marker":"[63]"},{"why":"gives the Weinberg operator Wilson coefficient $C_5$ from gluino-squark loops.","marker":"[50]"},{"why":"defines the loop functions $F_3$, $F_4$, $F_5$ used in the two-loop gluino self-energy contributions.","marker":"[62]"},{"why":"defines the TNMSSM field content, superpotential, and soft-breaking terms, including the triplet and singlet extensions.","marker":"[13–15]"},{"why":"generates the interaction vertices $C^{L,R}$ used in the loop amplitudes from the model's mass matrices.","marker":"[57–61]"}],"fun_headline_variants":["Heavy-quark EDMs up to 10^-22 e·cm in SUSY","Neutron EDM squeezes new SUSY CP phases","TNMSSM: b-quark EDM could hit 10^-22 e·cm","SUSY triplet-singlet boosts quark EDMs","New CP phases in TNMSSM face neutron EDM limits"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire constraint analysis rests on one hadronic formula, Eq. (2), that turns quark electric and chromoelectric dipole moments into the neutron's EDM; the paper fixes that formula's coefficients at their central values (0.5 and 12 MeV) and does not propagate the $\\pm0.5$ and $\\pm10\\,\\text{MeV}$ uncertainties, so a change in this hadronic input would shift the bounds on the new CP phases and the predicted heavy-quark EDMs.","fun_headline_variants_meta":{"raw":{"variants":["Heavy-quark EDMs up to 10^-22 e·cm in SUSY","Neutron EDM squeezes new SUSY CP phases","TNMSSM: b-quark EDM could hit 10^-22 e·cm","SUSY triplet-singlet boosts quark EDMs","New CP phases in TNMSSM face neutron EDM limits"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000374,"raw_usage":{"total_tokens":2026,"prompt_tokens":1006,"completion_tokens":1020,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":622,"completion_tokens_details":{"reasoning_tokens":921}},"tokens_in":622,"tokens_out":1020,"duration_ms":8738,"temperature":1.0,"reasoning_tokens":921,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:14:14.971195+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the bounds on $\\chi_d$ and $\\chi_t$ and the maxima of $d_b$ and $d_c$ using the full $\\pm0.5$ and $\\pm10\\,\\text{MeV}$ uncertainty band of Eq. (2) instead of the fixed central values the paper adopts; the claimed $10^{-22}\\,\\text{e}\\!\\cdot\\!\\text{cm}$ reach for $d_b$ survives only if it changes by less than an order of magnitude under that propagation.","supporting_citations":[{"cited_title":"Sala, JHEP 1403, 061 (2014)","cited_arxiv_id":null,"evidence_quote":"provides the improved upper bounds on $b$- and $c$-quark EDMs and the hadronic running used to convert neutron limits into heavy-quark constraints."},{"cited_title":"Pospelov and A","cited_arxiv_id":null,"evidence_quote":"supplies the neutron-EDM relation in Eq. (2) connecting quark EDMs, CEDMs, and the Weinberg operator."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the mercury-EDM formula $d_{\\text{Hg}}$ from quark CEDMs used in the numerical analysis."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"gives the two-loop Barr-Zee type contributions to the quark EDM."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"gives the Weinberg operator Wilson coefficient $C_5$ from gluino-squark loops."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"defines the loop functions $F_3$, $F_4$, $F_5$ used in the two-loop gluino self-energy contributions."}],"review_version":1}