{"id":"d11ce23b-c30f-40cb-b096-1b7fa6a6746f","arxiv_id":"2411.19469","paper_version":4,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"With optimized selection and an optimal observable, STCF could reach |d_tau| below 3.9e-18 e cm at 68% confidence, improving the Belle limit by about a factor of three.","lead":"This simulation study projects that the planned Super Tau-Charm Facility in China could measure the tau lepton's electric dipole moment with a sensitivity of about 4 x 10^-18 e cm after ten years, roughly three times better than the current best limit from Belle. It shows how machine-learning event selection and an optimized observable could extract this small signal from background.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quoted 3.89e-18 sensitivity rests on Sec. IV.C MadGraph samples preserving full tau spin correlations in rho-rho decays, but the paper never states how decays and polarimeter vectors are generated; this is the load-bearing unverified link.","rationale":"The paper is a careful feasibility study with internally consistent selection numbers and a standard optimal-observable method. The stress-test pass identifies one unverified procedural premise in the calibration chain: the MadGraph samples of Sec. IV.C must encode the same spin correlations as the analytic squared matrix element used to define O_Re. The text does not describe the tau decay generation or spin-density-matrix handling, so the numerical value of the sensitivity cannot be independently reconstructed from the paper. This is exactly the reader's weakest assumption, and it does not require changing the CONDITIONAL verdict; it does require the requested closure test or code release. I do not see an internal inconsistency in the selection or kinematic reconstruction; the issue is the missing link between the generator-level EDM amplitude and the reconstructed observable.","tokens_in":13538,"tokens_out":4824,"duration_ms":44710,"concrete_test":"Obtain the LHE event record or run cards for the samples behind Fig. 8 and Table III and check the tau decay provenance. Then generate two MadGraph/UFO samples with identical setup but (a) decays handled by a spin-correlation-preserving generator such as Tauola or MadSpin and (b) decays with tau spins averaged or decayed isotropically; recompute <O_Re> for d_tau = 0 and d_tau = 2e-16 e cm. If the fitted a_Re in Table III changes by more than about 10% between cases, or if the d_tau dependence in case (b) is not linear with the quoted slope, the 3.89e-18 projection is not supported. A minimal closure check is to compare truth-level MC values of <O_Re> against Eq. (7) using the analytic a_Re and b_Re.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The linear calibration in Table III and the resulting |d_tau| < 3.89e-18 e cm require the optimal observable O_Re = |M_inter|^2/|M_SM|^2 to be computed from reconstructed tau-decay products in events that carry the full spin-correlated production amplitude. The analytic framework in Eqs. (24)-(31) is explicit that the polarimeter vectors h_+ and h_- and the c_i^{mu nu} spin-correlation term enter the observable. Section IV.C describes only MadGraph with custom UFO model files to simulate tau-pair production, followed by Delphes detector simulation; it does not state how tau decays to pi+pi0 nu and pi-pi0 nu are generated, whether Tauola or MadSpin preserves spin correlations, or how the two-fold tau momentum ambiguity and the average over solutions are implemented in the calibration. If the MadGraph samples were generated with tau spins averaged at production, or with isotropic or uncorrelated decays, the d_tau dependence of <O_Re> would be altered or vanish entirely, so the fitted slope a_Re would not calibrate the EDM. The central numerical claim is therefore conditional on an unstated simulation procedure, not on disagreement with consensus.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a Monte Carlo sensitivity study of the tau lepton electric dipole moment (EDM) at the proposed Super Tau-Charm Facility (STCF). The authors simulate e+e- -> tau+tau- events with tau decays to rho nu (tau+ -> pi+ pi0 anti-nu_tau, tau- -> pi- pi0 nu_tau), develop a BDTG-based event selection achieving 80.0% signal purity with 6.3% efficiency, use a joint kinematic fit for particle pairing, reconstruct tau momenta analytically, and construct the optimal observable O_Re = |M_inter|^2 / |M_SM|^2. The central claim, stated in the abstract and Section V, is that after ten years of STCF operation the tau EDM can be constrained to |d_tau| < 3.89e-18 e cm at 68% confidence, based on a linear calibration of the mean optimal observable against d_tau using MadGraph samples with custom UFO model files and Delphes fast simulation.","tokens_in":13819,"tokens_out":4897,"duration_ms":45354,"significance":"If the quoted sensitivity is correct, this would be the most sensitive tau EDM projection to date, about a factor of three stronger than the current Belle limit near 1e-17 e cm, and it would be a useful input for the STCF physics program. The paper has genuine strengths: the optimal-observable framework is analytically derived and internally consistent, the selection chain is described in detail, the comparison with Belle yields a plausible statistical improvement, and the use of a linear calibration from MC samples is a standard and testable procedure. The main weakness is that the central numerical result depends on an incompletely documented simulation chain, and the quoted number is not directly derivable from the tabulated fit results as presented. The study is also statistical-only, with no systematic uncertainty assessment, which is acceptable for a first projection but should be stated explicitly.","major_comments":[{"comment":"The calibration of the optimal observable response rests on the MadGraph samples with custom UFO files reproducing the full spin correlations between tau production and tau decay. The analytic expressions in Eqs. (24)-(31) show that the observable depends on the polarimeter vectors h_+ and h_- and on the spin-correlation terms c_i^{mu nu}, which are constructed from the tau decay matrix elements. However, Section IV.C states only that MadGraph and custom UFO files are used to simulate tau-pair production, followed by Delphes detector simulation; it does not state how the tau decays to pi+ pi0 nu and pi- pi0 nu are generated, whether Tauola or MadSpin or another package preserves the spin correlations, or how the two-fold tau momentum ambiguity and the averaging over the two solutions are implemented in the calibration. If the MadGraph samples were generated with tau spins averaged at production or with uncorrelated decays, the dependence of <O_Re> on d_tau would be altered or could vanish, so the fitted slope a_Re in Table III would not calibrate the EDM. This unstated link is load-bearing for the central claim and must be documented and validated.","section":null},{"comment":"The quoted sensitivity |d_tau| < 3.89e-18 e cm is not derived in the text and is not directly reproducible from Table III. Taking the statistical uncertainty on the mean observable as delta<O_Re> = 1.088e-4 and the fitted slope a_Re = 3.06e-5 (10^-18 e cm)^-1, the simple ratio delta<O_Re>/a_Re gives 3.55e-18 e cm, not 3.89e-18 e cm. The additional factor of about 1.09 presumably comes from an unstated propagation of the fit parameter uncertainties or from another statistical formula, but this is not shown. The paper should present the exact formula used for the 68% confidence bound, including how the uncertainties on a_Re, b_Re, and the number of selected events are propagated.","section":null},{"comment":"The study provides no systematic uncertainty assessment. The selection has non-negligible background contamination (about 20% background, dominated by tau->pi pi0 pi0 nu), the tau momentum reconstruction has a two-fold ambiguity that is averaged over, and the photon and track resolutions are folded in through fast simulation. Any of these could bias the mean of the optimal observable or dilute the slope a_Re. Since the paper claims an 'estimated sensitivity' of 3.89e-18 e cm, it should at minimum state explicitly that this is a statistical-only projection and list the dominant expected sources of systematic uncertainty, with rough estimates where possible. Without this, the numerical claim is presented with unwarranted precision.","section":null}],"minor_comments":[{"comment":"The quantity v in the amplitude Mprod = M_SM + (Re(d_tau)/v) M_Re is not defined; it should be stated explicitly (presumably the Higgs vacuum expectation value) since it enters the definition of the slope in Eq. (7).","section":null},{"comment":"The phrase 'we only keep the real part of Re(d_tau)' is confusing because Re(d_tau) is already the real part; the authors should say they assume d_tau is real and retain only the linear interference term.","section":null},{"comment":"The sign conventions in Eq. (31) for the a_i, b_i, and c_i coefficients are not defined; they should be related to the quantities introduced in Eq. (4) or defined in the text.","section":null},{"comment":"The legend in the left panel of Fig. 8 is described as 'Red: d_tau<0; Green: d_tau>0' in the caption, but the colored curves and the definition of the ratio r are not clearly identified in the figure itself; please clarify the plot and state where r is used.","section":null},{"comment":"The paper does not state how many reconstructed signal events are used in the calibration fits for each injected d_tau value, nor the corresponding statistical uncertainty on each mean <O_Re>; providing this would improve reproducibility.","section":null},{"comment":"There is a typo in the text before Eq. (31): 'produciton' should be 'production'.","section":null}],"recommendation":"major_revision","confidential_remarks":"The paper fits the scope of a high-energy physics experiment journal and is a reasonable first sensitivity study. The main reason for major revision is not disagreement with the physics but the incomplete documentation of the MC calibration chain and the unreproducible numerical result. The self-citation in Refs. [14] and [15], which share a coauthor with the present paper, is mild and not problematic, but the authors may wish to add a brief note on the relation to those works. No data or code are provided, which is not unusual for this type of study, but the analysis would be stronger if the custom UFO model files and the decay-generation configuration were made available."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Sun, Wu, and Zhou have produced a serious feasibility study for a tau EDM measurement at STCF. The headline number, |d_tau| < 3.89e-18 e cm at 68% CL after ten years, is about three times better than Belle's existing limit. The paper deserves a referee's time.\n\nWhat is actually new: this is the first STCF-specific projection using a modern selection chain (BDTG photon and event selection), a joint kinematic fit for particle pairing, and an analytic tau momentum reconstruction. The event selection tables are internally consistent, and the kinematic fitting raises the correct-pairing rate to 82.5%. The spin density matrix and optimal observable treatment follow the established Belle framework but are re-derived cleanly. Credit is due for the careful step-by-step selection accounting.\n\nThe soft spots are real but mostly about documentation, not about the physics. The largest one: Section IV.C says only that MadGraph with a custom UFO file 'simulates the production' of tau pairs, and that Delphes performs fast detector simulation. It never states how the taus are decayed, whether Tauola or MadSpin preserves spin correlations, or how the two-fold momentum ambiguity is averaged in the calibration. The observed linear shift of <ORe> with d_tau (Fig. 8, right) shows the samples must carry some spin information, so the stress-test concern does not land as a fatal flaw. But the slope's magnitude, and hence the 3.89e-18 projection, cannot be independently checked without that information. This is a reproducibility gap, not a contradiction.\n\nA smaller but annoying issue: the quoted limit does not match the table. With a_Re = 3.06e-5 and sigma_<ORe> = 1.088e-4, the sensitivity is 3.55e-18, not 3.89e-18. The paper never shows the extra propagation (uncertainty on aRe? a scale factor?). That needs a sentence or two.\n\nThere are no systematics, which is forgivable in a projection but should be listed. And no code release, which would settle the spin-correlation question.\n\nWho should read this: anyone planning tau physics at STCF, and theorists who want a concrete target for tau EDM predictions. It is not a measurement paper, and it does not settle any BSM question. It is a carefully done projection with one undocumented step.\n\nRecommendation: send it to peer review. The referee should ask for the simulation chain details and the error calculation, but the core work is solid and the result is useful for the STCF community.","headline":"A careful STCF tau EDM projection that deserves refereeing, but the key calibration step is under-documented and the quoted limit does not trace from the tables.","tokens_in":14354,"tokens_out":5079,"would_cite":true,"duration_ms":43578,"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":"A planned tau-charm collider could measure the tau lepton EDM to 3.9e-18 e cm, about three times tighter than today's best limit.","keywords":["tau lepton","electric dipole moment","CP violation","e+e- collisions","optimal observable","spin correlations","Monte Carlo sensitivity study","tau-charm collider"],"falsifier":"Generate $e^+e^-\\to\\tau^+\\tau^-$ events with the EDM operator included in a custom model, and compute the optimal observable distribution twice: once with tau decays carrying the full spin information from the production amplitude, and once with decays treated without spin correlations. If the mean shift of $\\mathcal{O}_{\\mathrm{Re}}$ as a function of $d_\\tau$ differs between the two treatments, or deviates from the linear relation with the slope quoted in the paper's Table III, then the projected sensitivity is not reproducible.","tokens_in":13327,"feed_emoji":"⚛️","tokens_out":11289,"duration_ms":90760,"temperature":0.7,"pith_summary":"The paper argues that a future high-luminosity electron-positron collider in the tau-charm energy region can measure the tau lepton's electric dipole moment to $|d_\\tau| < 3.89\\times 10^{-18}\\,e\\cdot\\mathrm{cm}$ at 68% confidence after ten years of data taking. This would be about three times tighter than the current best experimental limit. To get there, the paper constructs a complete analysis chain: machine-learning photon and event selection, joint kinematic fitting to pair decay products, analytic tau momentum reconstruction, and an optimal observable derived from the squared spin density matrix. The observable's mean value is calibrated against Monte Carlo samples with different EDM values, giving a linear relation that converts a measured shift in the observable into an EDM value. The projected sensitivity would probe CP-violating new physics in the tau sector beyond the Standard Model.","feed_headline":"Tau EDM sensitivity projected to 3.9e-18 e cm","feed_subtitle":"Ten years of data at a planned tau-charm collider could beat today's best tau EDM limit by three times.","key_machinery":"The load-bearing object is the optimal observable $\\mathcal{O}_{\\mathrm{Re}} = |M_{\\mathrm{inter}}|^2/|M_{\\mathrm{SM}}|^2$, where $|M_{\\mathrm{inter}}|^2$ is the SM-EDM interference term in the squared production amplitude and $|M_{\\mathrm{SM}}|^2$ is the Standard Model term. Constructed from the reconstructed tau momenta and the polarimeter vectors of the hadronic decays, its mean over phase space is $\\langle\\mathcal{O}_{\\mathrm{Re}}\\rangle = a_{\\mathrm{Re}}\\,\\mathrm{Re}(d_\\tau)/v + b_{\\mathrm{Re}}$, so the EDM is read off from a measured shift in the mean. The squared spin density matrix decomposition $|M_i|^2 = |M_i|^2(1+h_i^+\\cdot s_+ + h_i^-\\cdot s_- + c_i\\cdot s_+ s_-)$ connects the production spin correlations to the decay polarimeters; the two-fold ambiguity in the reconstructed tau momentum is handled by averaging the observable over the two solutions, following the earlier experimental approach.","core_discovery":"The central claim is that using the $\\tau^\\pm\\to\\pi^\\pm\\pi^0\\nu_\\tau$ (rho-rho) decay channel, the planned collider can constrain the tau EDM to $|d_\\tau|<3.89\\times10^{-18}\\,e\\cdot\\mathrm{cm}$ at 68% confidence after ten years, roughly three times stronger than the limit from the previous B-factory experiment. The analysis chain achieves 80.0% signal purity with 6.3% efficiency, and the joint kinematic fitting pairs the final-state particles correctly in 82.5% of signal events. The EDM sensitivity comes from the optimal observable $\\mathcal{O}_{\\mathrm{Re}} = |M_{\\mathrm{inter}}|^2/|M_{\\mathrm{SM}}|^2$, whose phase-space mean is a linear function of $\\mathrm{Re}(d_\\tau)/v$; Monte Carlo samples with several input EDM values fix the slope of that line, and the statistical uncertainty on the mean observable translates into the quoted limit.","pith_inferences":["The paper leaves open how tau decays and spin correlations are realized in the calibration samples; a truth-level comparison between samples with spin correlations included and excluded would settle whether the fitted slope $a_{\\mathrm{Re}}$ is biased.","Resolving the two-fold tau-momentum ambiguity with vertex-detector impact parameters, which the paper identifies as a possible future improvement, could push the limit below the quoted value by removing the averaging that dilutes the observable.","The same optimal-observable pipeline could, in principle, be applied to data already collected at existing tau-charm experiments, giving an intermediate EDM measurement before the new collider turns on.","Because the dominant remaining background is $\\tau^\\pm\\to\\pi^\\pm\\pi^0\\pi^0\\nu_\\tau$ at about 14%, a further selection optimized against extra neutral pions would raise purity and may improve sensitivity beyond the quoted projection."],"forward_implications":["A null result after ten years would set $|d_\\tau| < 3.89\\times10^{-18}\\,e\\cdot\\mathrm{cm}$ at 68% CL, roughly three times stronger than the current best limit.","The selection chain yields 80.0% signal purity at 6.3% efficiency, corresponding to about $1.4\\times10^7$ rho-rho signal events per year, two orders of magnitude more than the earlier B-factory sample.","Because the mean optimal observable is linear in $\\mathrm{Re}(d_\\tau)$, any future measurement of $\\langle \\mathcal{O}_{\\mathrm{Re}}\\rangle$ with the quoted uncertainty can be converted directly into an EDM value or limit.","The reconstructed tau momentum, after averaging the two analytical solutions, has a transverse-momentum relative deviation with FWHM 0.10 and a direction RMS of about 10 degrees, enough to keep the optimal-observable calibration stable."],"supporting_citations":[{"why":"Provides the current experimental EDM limit and the rho-rho optimal-observable analysis whose method this study adapts and extends.","marker":"[11]"},{"why":"Supplies the EDM Lagrangian and the earlier experimental search formalism used to define the squared matrix element.","marker":"[14]"},{"why":"Establishes the tau-factory EDM search framework, including the optimal observable approach for $e^+e^-\\to\\tau^+\\tau^-$, that this study applies.","marker":"[15]"},{"why":"Introduces the optimal observable construction that maximizes sensitivity to a form-factor-like parameter.","marker":"[16]"},{"why":"Generates the $e^+e^-\\to\\tau^+\\tau^-$ production sample with radiative corrections used as the study's Monte Carlo dataset.","marker":"[17]"},{"why":"Provides the tau decay model and polarimeter-vector formalism used to connect tau spin to final-state momenta.","marker":"[18]"},{"why":"Simulates EDM-modified production with a custom model file, producing the samples whose optimal-observable means calibrate the linear fit.","marker":"[24]"},{"why":"Performs the fast detector simulation that turns generated events into reconstructed quantities used in the selection.","marker":"[25]"},{"why":"Defines the planned collider's luminosity, energy range, and detector performance assumptions from which the ten-year event yield is derived.","marker":"[13]"}],"fun_headline_variants":["STCF projected to probe tau EDM to 3.9e-18","Tau EDM limit cut threefold by STCF study","ML boosts tau EDM sensitivity to 3.9e-18","STCF study sets tau EDM target of 3.9e-18"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The projection assumes the Monte Carlo samples used to calibrate the optimal observable reproduce the full spin correlation between tau production and tau decay exactly as in the analytic matrix element; the paper does not describe how tau decay and spin are handled in those samples, so if the simulation misses those correlations the fitted slope and the quoted limit would not hold.","fun_headline_variants_meta":{"raw":{"variants":["STCF projected to probe tau EDM to 3.9e-18","Tau EDM limit cut threefold by STCF study","ML boosts tau EDM sensitivity to 3.9e-18","STCF study sets tau EDM target of 3.9e-18"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000677,"raw_usage":{"total_tokens":3105,"prompt_tokens":997,"completion_tokens":2108,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":613,"completion_tokens_details":{"reasoning_tokens":2027}},"tokens_in":613,"tokens_out":2108,"duration_ms":15621,"temperature":1.0,"reasoning_tokens":2027,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T10:11:46.135774+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Generate $e^+e^-\\to\\tau^+\\tau^-$ events with the EDM operator included in a custom model, and compute the optimal observable distribution twice: once with tau decays carrying the full spin information from the production amplitude, and once with decays treated without spin correlations. If the mean shift of $\\mathcal{O}_{\\mathrm{Re}}$ as a function of $d_\\tau$ differs between the two treatments, or deviates from the linear relation with the slope quoted in the paper's Table III, then the projected sensitivity is not reproducible.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the current experimental EDM limit and the rho-rho optimal-observable analysis whose method this study adapts and extends."},{"cited_title":"04 (2022) 110","cited_arxiv_id":null,"evidence_quote":"Supplies the EDM Lagrangian and the earlier experimental search formalism used to define the squared matrix element."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the tau-factory EDM search framework, including the optimal observable approach for $e^+e^-\\to\\tau^+\\tau^-$, that this study applies."},{"cited_title":"Achasov, X.C","cited_arxiv_id":null,"evidence_quote":"Introduces the optimal observable construction that maximizes sensitivity to a form-factor-like parameter."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Generates the $e^+e^-\\to\\tau^+\\tau^-$ production sample with radiative corrections used as the study's Monte Carlo dataset."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the tau decay model and polarimeter-vector formalism used to connect tau spin to final-state momenta."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Performs the fast detector simulation that turns generated events into reconstructed quantities used in the selection."},{"cited_title":"Pospelov, A","cited_arxiv_id":null,"evidence_quote":"Defines the planned collider's luminosity, energy range, and detector performance assumptions from which the ten-year event yield is derived."}],"review_version":1}