{"id":"b8505641-153c-43e5-a933-f4497e05689f","arxiv_id":"2501.01286","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"First numerical study of color-restoration phase transitions in a radiative neutrino mass leptoquark model, predicting LISA/BBO/DECIGO-detectable gravitational waves for leptoquark masses near 1.5 TeV.","lead":"This paper studies gravitational waves from early-universe phase transitions in a leptoquark model that also explains neutrino masses. It predicts that a detectable signal in the millihertz to 0.1 hertz band would come with leptoquark masses near 1.5 TeV, testable at future space-based detectors and colliders.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The LISA-reachable region is populated by (h,s,r)->0 transitions with mUS/mu ~ 2, but the only scale-variation test (Fig. 7) is performed on r->h benchmarks, leaving the headline detection claim without a validity check in the relevant vacuum sector.","rationale":"The reader identifies the reliability of the high-temperature EFT at mUS/mu ~ 2 as the weakest assumption, and I agree that this is the most load-bearing point: the central detection claim explicitly relies on relaxing Eq. (4.11), and the paper's own Sec. V states that values around mUS/mu ~ 2 are 'sufficient' while the precise perturbativity boundary is left undetermined. My stress-test adds a concrete specificity: the two benchmark points used for the scale-variation test (Fig. 7) are r -> h transitions, while the LISA-reachable region is dominated by (h,s,r) -> 0 and (h,s,r) -> h transitions, with the strongest (h,s,r) -> 0 points at mUS/mu ~ 2. Therefore the existing validation does not cover the vacuum sector responsible for the headline prediction. This does not prove the prediction is wrong, and the paper is honest about the approximation and its limitations, but it does mean the detectability claim currently rests on an untested extrapolation. Other concerns, such as the unshown neutrino-oscillation compatibility scan and the unreleased Dratopi code, are real but secondary: the former is an existence claim supported by the structure of the model and prior work, and the latter affects reproducibility rather than correctness. The reader's CONDITIONAL verdict is appropriate; my analysis does not change it, but it strengthens the specific condition that must be checked before the central claim can be accepted.","tokens_in":21513,"tokens_out":7900,"duration_ms":83222,"concrete_test":"Take the scan points of type (h,s,r) -> 0 that lie above the LISA PISC or have mUS/mu ~ 2, and repeat the scale-variation test of Fig. 7 for these points, varying C4d = mu/(pi T) over [0.5, 2] and recomputing the peak amplitude h^2 Omega_GW. Also re-run the scan with the strict condition mUS/mu < 1 enforced and count how many points remain above the LISA/BBO/DECIGO PISCs. If the (h,s,r) -> 0 points show large amplitude shifts under the scale variation, or if enforcing Eq. (4.11) removes all LISA-reachable points, then the abstract's 'detectable signal' claim is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract's 'detectable signal' claim depends on points that violate the high-temperature perturbativity bound stated in Eq. (4.11), mUS/mu < 1, with the majority of the strong-signal points clustering at mUS/mu ~ 2 (Sec. IV B, Figs. 5-6). The authors explicitly relax this bound, calling deviations by a factor of a few 'acceptable' and deferring a finer analysis. The only quantitative evidence for this relaxation is the scale-variation test in Fig. 7, which varies mu/(pi T) over [0.5, 2] for two benchmarks of the r -> h transition pattern. However, Sec. IV C and Fig. 11a show that the strongest LISA-reachable signals are dominated by (h,s,r) -> 0 and (h,s,r) -> h transitions, and that the strongest (h,s,r) -> 0 points specifically have mUS/mu ~ 2. Thus the scale-variation test is performed on a vacuum configuration that is not the one carrying the headline sensitivity, and it does not check the convergence of the NLO dimensional-reduction expansion when the ultrasoft mass exceeds the hard matching scale by roughly 2 pi T. Because the effective potential (3.5) and the resulting alpha and beta/H, and hence the predicted h^2 Omega_GW, are computed within this EFT, the detection claim is an extrapolation outside the demonstrated validity domain of the computation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies cosmological first-order phase transitions in a scalar leptoquark extension of the SM (an R2 doublet and an S1 singlet) that radiatively generates Majorana neutrino masses. The authors match the model to the SM, construct the finite-temperature effective potential by NLO dimensional reduction using DRalgo/Dratopi, scan roughly O(10^5) parameter points, and compute the stochastic gravitational wave background from sound waves using the LISA Cosmology Working Group templates. They identify color-breaking 'color-restoration' transitions as the strongest sources and claim that a detectable SGWB in the mHz–0.1 Hz band is correlated with LQ masses near 1.5 TeV. The paper is candid about the main caveats: the high-temperature perturbativity parameter m_US/mu reaches about 2 in the detectable region, the neutrino-data verification is not described, and the bubble wall velocity is fixed at 0.95.","tokens_in":21863,"tokens_out":6941,"duration_ms":73207,"significance":"If the detection claim survives closer scrutiny, the paper is valuable: it gives one of the first quantitative GW analyses of color-restoring FOPTs in a radiative neutrino-mass model and yields concrete, falsifiable correlations between GW observability and LQ masses/couplings. The use of a dimensionally reduced NLO effective potential, a large Monte-Carlo scan, explicit detector PISC curves, and candid acknowledgement of the EFT validity limitation are strengths. The workflow is not circular: the GW parameters are computed from the effective potential after matching to the SM, not fitted to the claimed result. The main risk is not internal inconsistency but extrapolation beyond the demonstrated regime of the EFT, together with two auxiliary assumptions (neutrino-data reproduction and supersonic detonation) that would benefit from dedicated support.","major_comments":[{"comment":"The paper's detectability claim is carried by (h,s,r)->0 and (h,s,r)->h transitions, whose strongest points have m_US/mu approximately 2, yet the scale-variation check in Fig. 7 is performed only for r->h benchmarks. Since Eq. (4.11) defines the EFT validity condition and the effective potential (3.5) is computed inside that EFT, the factor-of-two excursion is an extrapolation for exactly the transitions that populate the LISA-reachable region. Please either (i) repeat the scale-variation test on the dominant vacuum channels, ideally with several benchmarks spanning the LISA-reachable region, (ii) estimate the size of the missing NNLO or higher-order terms for those channels, or (iii) explicitly restrict the headline claim to points satisfying Eq. (4.11).","section":"Sec. IV B, Eq. (4.11); Sec. IV C, Fig. 11a"},{"comment":"The text states 'We have verified that for every viable FOPT... there is at least one solution where the entries of Theta and Omega are small and simultaneously reproduce neutrino oscillation data,' but no procedure, observable definition, or statistical criterion is given, and no reference to a companion calculation is supplied. Because the abstract advertises the model as explaining neutrino oscillation data and the scan accepts or rejects points based on this verification, the claim needs to be reproducible: specify the fit to oscillation parameters, the allowed ranges, and how 'viable FOPT' points were enumerated. Otherwise the model interpretation of the GW predictions is unsupported.","section":"Sec. III A"},{"comment":"All GW amplitudes are computed with xi_w fixed to 0.95, which enters both the mean bubble separation R* and the sound-shell width Delta_w. The justification via large alpha and Ref. [52] is plausible, but since the quantitative claim is that a specific region is observable, the authors should show how the peak amplitude and the LISA/BBO/DECIGO reach change when xi_w is varied over a conservative range (e.g., 0.7-0.99) for a few benchmark points in the detection region. This would separate the EFT-validity uncertainty from the hydrodynamic-model uncertainty.","section":"Sec. IV A and Eq. (3.28)"}],"minor_comments":[{"comment":"The caption of Fig. 10(b) describes the example as an (r,h)->h transition, while Sec. IV C cites (r,h)->0 as a typical color-breaking example; please make the notation consistent.","section":"Fig. 10, Sec. IV C"},{"comment":"The numerical factor appears as '10.^-5' in Eq. (3.28); it should read 10^{-5}.","section":"Eq. (3.28)"},{"comment":"The text states 'beta/H is less than or similar to 105' where the upper limit should evidently read 10^5; please correct the superscript to match the subsequent discussion.","section":"Sec. IV A"},{"comment":"Dratopi is cited as 'To appear' [49], and the numerical results depend on this package; the authors should release it or provide a public artifact (e.g., a repository with the model export and scan scripts) before or at acceptance to make the analysis reproducible.","section":"Sec. IV, Ref. [49]"},{"comment":"The Monte-Carlo improvement weights a and b are set to 10 without a convergence or bias check; since the histograms in Fig. 11 are based on this scan distribution, a brief statement that the scan is exploratory and not intended as a Bayesian posterior would help.","section":"Appendix B, Eq. (B1)"}],"recommendation":"major_revision","confidential_remarks":"I am sympathetic to the authors' candid treatment of the EFT-validity caveat, and I do not regard the m_US/mu approximately 2 excursions as automatic grounds for rejection. However, the central observable claim is only as strong as the validity of the dimensional-reduction expansion in the dominant vacuum channel; without the requested checks, the headline sentence in the abstract outruns the evidence. The neutrino-data verification and the unreleased Dratopi package are secondary but should be addressed before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a competent first pass at GW predictions from color-restoring FOPTs in the minimal radiative-neutrino-mass leptoquark model, and the central find—LISA/DECIGO/BBO-reachable signals with LQ masses around 1.5 TeV—is plausible but sits outside the strict domain of the EFT calculation in exactly the region that matters.\n\nWhat's new: the paper applies dimensional reduction (DRalgo) plus CosmoTransitions to a specific LQ model with three VEVs, classifies eleven transition patterns, and identifies the (h,s,r)->0 and (h,s,r)->h channels as the strongest. It is candid: Sec. IV B and V explicitly flag mUS/mu ~ 2 and say the perturbativity criterion is not a strict threshold. That transparency is real and should be credited. The scan is large (~1e5 points) and the GW pipeline follows the LISA Cosmology Working Group templates, so the machinery is sound.\n\nSoft spots: first, the stress-test point is correct. The scale-variation test in Fig. 7 covers two r->h benchmarks, one at mUS/mu ~ 2, but the headline signals are dominated by (h,s,r)->0 points with mUS/mu ~ 2. That leaves the central detectability claim without a direct validity check in the relevant vacuum sector. The line about deviations by a factor of a few being acceptable is an opinion, not a demonstration. Second, Dratopi is unreleased; the numbers cannot be reproduced, though the public pieces (DRalgo, CosmoTransitions) are standard. Third, the statement that every viable FOPT can reproduce neutrino data is asserted, not shown; that is a minor point, but a plot or a footnote with the scan would help. Fourth, the abstract's 'detectable signal' wording is stronger than the body's caveats, which is common but worth tightening.\n\nNet: no fatal flaw. The central argument—that this model has strong color-restoring FOPTs with potentially observable GWs—holds up as a plausible claim, but the precision on the LISA-reachable boundary is not as solid as the abstract suggests. The paper deserves peer review; a referee should ask for either an extended scale-variation test in the (h,s,r)->0 sector or a revised claim that the mUS/mu ~ 2 region is indicative rather than definitive. I would cite it if I worked on BSM phase transitions, and it is a good reading-group case for the dangers of EFT-validity boundaries.","headline":"Solid, candid first numerical study of color-restoring FOPTs in a radiative-neutrino-mass LQ model with a plausible LISA-reachable region, but the headline claim rests on points at the edge of the EFT's validity and the one robustness check does not cover that sector.","tokens_in":22447,"tokens_out":2952,"would_cite":true,"duration_ms":30504,"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 paper argues that a minimal scalar leptoquark model of radiative neutrino masses predicts a stochastic gravitational-wave background from first-order color-restoring phase transitions in the mHz–0.1 Hz band, with leptoquark masses…","keywords":["gravitational waves","stochastic gravitational wave background","first-order phase transition","color restoration","leptoquark","radiative neutrino mass","dimensional reduction","LISA"],"falsifier":"Measure the stochastic gravitational-wave background in the mHz–0.1 Hz band with LISA, BBO, or DECIGO and simultaneously search for scalar leptoquarks near 1.5 TeV at the high-luminosity LHC; null results at the peak amplitudes shown in Figs. 5 and 6 would rule out the paper's detectable parameter region. A dedicated high-order (NNLO) or non-perturbative calculation of the effective potential at $m_{\\rm US}/\\mu\\approx 2$ could also check whether the strong first-order transition survives in that regime.","tokens_in":21227,"feed_emoji":"🌊","tokens_out":8868,"duration_ms":83136,"temperature":0.7,"pith_summary":"This paper tries to show that a minimal extension of the Standard Model—one pair of scalar leptoquarks, hypothetical particles that carry both quark and lepton quantum numbers—can do three things at once: radiatively generate neutrino masses, pass existing flavor-physics constraints, and drive a first-order phase transition in the early Universe in which color symmetry is briefly broken and then restored. If that is right, the transition would emit a stochastic gravitational-wave background in the mHz–0.1 Hz band that planned observatories LISA, BBO, and DECIGO could detect. The paper's central, testable correlation is that detectable signals require color restoration and leptoquark masses near 1.5 TeV, putting the scenario within reach of both gravitational-wave astronomy and collider searches. The analysis also presents the first application of a new automated dimensional-reduction pipeline, Dratopi, for phase-transition studies.","feed_headline":"Color restoration may ring in mHz gravitational waves","feed_subtitle":"A radiative neutrino mass model predicts LISA-visible signals from first-order color-restoring transitions with 1.5 TeV leptoquarks.","key_machinery":"The load-bearing object is the next-to-leading-order, dimensionally reduced finite-temperature effective potential of the three-scalar (Higgs plus two leptoquarks) model, computed through the paper's Dratopi pipeline. Dimensional reduction integrates out heavy high-temperature modes to build a three-dimensional effective theory at the ultrasoft scale; from this potential the paper computes the three-dimensional bounce action, nucleation and percolation temperatures, transition strength $\\alpha$, inverse duration $\\beta/H$, and the sound-wave gravitational-wave spectrum. A validity check tracks the ratio $m_{\\rm US}/\\mu$ of the largest ultrasoft boson mass to the hard matching scale, and points with $m_{\\rm US}/\\mu \\lesssim 2$ are used for the detectability claims.","core_discovery":"The paper's central claim is that the strongest first-order phase transitions in this leptoquark model are not ordinary electroweak transitions but color-breaking ones: at high temperature one or more of the scalar leptoquark fields develops a vacuum expectation value that breaks $SU(3)_C$, and as the Universe cools the color symmetry is restored. These color-restoring transitions produce gravitational-wave spectra with peak frequencies between roughly $10^{-3}$ and $0.1$ Hz, and a sizable part of the model's parameter space falls within the peak-integrated sensitivity of LISA, BBO, and DECIGO. In the detectable region the three physical leptoquark masses sit near $1.5$ TeV, the trilinear coupling $a_1$ is of order $10^3$ GeV, the mixed quartic couplings are of order $10^{-2}$, and the doublet self-coupling $\\lambda_R$ is of order $1$; the ratio $m_{\\rm US}/\\mu$ characterizing high-temperature perturbativity is close to $2$. The paper also catalogs eleven viable phase-transition patterns, eight of which involve color breaking in the high-temperature phase.","pith_inferences":["If the $m_{\\rm US}/\\mu\\approx 2$ assumption is relaxed, the predicted peak amplitudes in the most sensitive region could shift substantially; a systematic scan that enforces $m_{\\rm US}/\\mu<1$ would show how much of the LISA-accessible region survives.","The color-restoration mechanism is not tied to the specific $S_1+\\tilde R_2$ matter content; other radiative neutrino mass models containing colored scalars could exhibit the same correlation between mHz gravitational waves and TeV-scale scalar masses, making the predicted band a generic target.","An observed stochastic background in this band could act as a coarse 'mass spectrometer' for colored scalars: the peak frequency correlates with the transition temperature, which is set by leptoquark masses, so combined collider exclusions and gravitational-wave observations would tightly constrain the model."],"forward_implications":["If the prediction is correct, the stochastic background searched for by LISA, BBO, and DECIGO in the mHz–0.1 Hz band is a direct probe of color restoration in this model, not just of electroweak-scale physics.","A detection would single out leptoquark masses near 1.5 TeV and a particular corner of the scalar potential, giving LHC Run-3 and HL-LHC searches concrete targets.","The same phase transition is tied to the radiative neutrino mass mechanism through the trilinear coupling $a_1$, so a gravitational-wave signal would indirectly support the loop-level origin of neutrino masses in this construction.","Transitions that preserve color throughout are predicted to be too weak to detect ($h^2\\Omega_{\\rm GW}^{\\rm peak}\\approx 10^{-21}$), so a future detection would favor color-breaking histories and exclude those color-preserving alternatives."],"supporting_citations":[{"why":"Defines the minimal $S_1+\\tilde R_2$ leptoquark model, its radiative neutrino mass mechanism, and flavor-physics constraints that set the scanned parameter ranges.","marker":"[16]"},{"why":"Provides the automated dimensional-reduction computation that builds the NLO 3d effective potential used for all phase-transition results.","marker":"[41]"},{"why":"Supplies the bounce-action solver that the paper modifies to trace phases and compute tunneling actions in three-scalar-field space.","marker":"[50]"},{"why":"Gives the sound-wave and bubble-collision gravitational-wave spectral templates used to compute peak amplitudes and frequencies.","marker":"[6]"},{"why":"Provides the peak-integrated sensitivity curves for LISA, DECIGO, and BBO against which the predicted signals are judged.","marker":"[48]"},{"why":"Motivates the percolation-temperature choice and supplies the false-vacuum fraction, energy budget, and percolation condition formulas.","marker":"[42]"},{"why":"Provides the semi-analytic efficiency-coefficient fits used to estimate the sound-wave energy budget $\\kappa$.","marker":"[47]"},{"why":"Defines the region of the $\\alpha$–$\\beta/H$ plane where the effective-field-theory treatment gives genuine first-order transitions, used to impose the $\\beta/H<10^5$ cut.","marker":"[12]"},{"why":"Earlier study of finite-temperature color breaking that motivates the color-restoration scenario examined here.","marker":"[28]"}],"fun_headline_variants":["Color-restoring phase transitions may ring in mHz GWs","Leptoquark color restoration may produce mHz GWs","LISA-visible GWs from color-restoring leptoquark transitions","Color-restoring first-order transitions emit mHz GWs","GWs ring out when color symmetry is restored in leptoquark model"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The most detectable predictions assume the calculation stays reliable when the relevant particle masses are about twice the scale set by the high-temperature matching procedure, a regime reached only by relaxing the conventional validity condition.","fun_headline_variants_meta":{"raw":{"variants":["Color-restoring phase transitions may ring in mHz GWs","Leptoquark color restoration may produce mHz GWs","LISA-visible GWs from color-restoring leptoquark transitions","Color-restoring first-order transitions emit mHz GWs","GWs ring out when color symmetry is restored in leptoquark model"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00076,"raw_usage":{"total_tokens":3379,"prompt_tokens":956,"completion_tokens":2423,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":572,"completion_tokens_details":{"reasoning_tokens":2334}},"tokens_in":572,"tokens_out":2423,"duration_ms":15972,"temperature":1.0,"reasoning_tokens":2334,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:31:16.230794+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the stochastic gravitational-wave background in the mHz–0.1 Hz band with LISA, BBO, or DECIGO and simultaneously search for scalar leptoquarks near 1.5 TeV at the high-luminosity LHC; null results at the peak amplitudes shown in Figs. 5 and 6 would rule out the paper's detectable parameter region. A dedicated high-order (NNLO) or non-perturbative calculation of the effective potential at $m_{\\rm US}/\\mu\\approx 2$ could also check whether the strong first-order transition survives in that regime.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the sound-wave and bubble-collision gravitational-wave spectral templates used to compute peak amplitudes and frequencies."},{"cited_title":"First Order Color Symmetry Breaking and Restoration Triggered by Electroweak Symmetry Non-restoration","cited_arxiv_id":"2112.13580","evidence_quote":"Motivates the percolation-temperature choice and supplies the false-vacuum fraction, energy budget, and percolation condition formulas."},{"cited_title":"improvement","cited_arxiv_id":null,"evidence_quote":"Defines the region of the $\\alpha$–$\\beta/H$ plane where the effective-field-theory treatment gives genuine first-order transitions, used to impose the $\\beta/H<10^5$ cut."},{"cited_title":"On interplay between flavour anomalies and neutrino properties","cited_arxiv_id":"2206.01674","evidence_quote":"Earlier study of finite-temperature color breaking that motivates the color-restoration scenario examined here."}],"review_version":1}