{"id":"62dcbf1d-9400-446e-8070-f7b0b3b51a9d","arxiv_id":"2412.12267","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"A Z4-symmetric Dirac seesaw makes the imaginary part of a seesaw scalar a stable dark matter candidate and links its phenomenology to gravitational wave and CMB observables.","lead":"The paper constructs a neutrino-mass and dark matter model from the same small set of new particles, and predicts gravitational wave signals from the decay of domain walls created when a discrete symmetry breaks. A generalist might read it because the model also makes testable predictions for dark matter searches and cosmic microwave background measurements of extra radiation.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Zero-bias v1–v3 domain walls remain exactly degenerate under the CP-even bias term, so the paper's explicit exclusion of them in Sec. 3 leaves the central DW-annihilation claim unproven; a surviving network would be cosmologically problematic.","rationale":"The reader's weakest_assumption identifies the same gap, and our stress test confirms it is the most load-bearing issue. All other components (seesaw formula, DM relic, Dirac leptogenesis, Neff) are internally consistent parameter scans; the GW prediction is the one place where a qualitative cosmological statement rests on an unproven assumption. The paper itself flags the assumption in the quoted sentence of Sec. 3, so this is not a manufactured objection. A false-vacuum decay calculation or lattice simulation could resolve it; hence CONDITIONAL is the right verdict, with a request for either a CP-violating bias (which would require a DM-stability reanalysis) or a dynamical demonstration that v1–v3 walls are eliminated before they dominate.","tokens_in":36422,"tokens_out":11344,"duration_ms":108074,"concrete_test":"Run a 3D lattice simulation (e.g., with CosmoLattice) of the scalar potential Eq. (2.1) plus the bias Eq. (3.6), starting from random initial fields after the Z4-breaking transition, with benchmark values λρ=0.1, λ1=λρ/72, µ1=0.01 GeV, vρ=10^5 GeV, and ε ∈ [10^-26, 10^-21], and track the energy density stored in walls separating the v1 and v3 minima. If this component scales as a^-1 and survives past the biased-wall annihilation time or dominates before BBN, the central claim fails; if it decays via v0 interposition, the exclusion in Sec. 3 is justified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central cosmological claim is that the Z4-breaking bias V_bias = -√2 ε vρ^4 cosθ (Eq. 3.6) makes all catastrophic domain walls annihilate. This fails for the v1–v3 walls: Eq. (3.7) lists δV13 = 0, and the section explicitly states: 'The domain wall between v1 and v3 is excluded from our analysis since the associated bias, δV13 = 0.' Because V_bias is CP-even and the model is taken in the CP-conserving limit, the minima v1 and v3 remain exactly degenerate by the same C_dark/CP symmetry that stabilizes the DM candidate χ. A percolating network of such walls is the standard Z2 scaling solution and has no pressure difference to drive annihilation; if the v1/v3 metastable regions decay to the global minimum v0 only via bubble nucleation, the timescale can exceed the age of the universe for the tiny ε values used (e.g., 10^-26–10^-21). The QG origin of the bias discussed in Appendix D is also CP-even (e.g., c1(H†H)^2ρ/Λ), so it does not lift the degeneracy. The paper provides no dynamical argument or simulation for this wall type; simply excluding it leaves open a surviving domain-wall network that would spoil the advertised GW scenario, neutrino/DM cosmology, and BBN.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a minimal Type-I Dirac seesaw extension of the SM containing heavy Dirac fermions N_L,N_R, right-handed neutrinos nu_R, a real scalar eta, and a complex scalar rho, with a global Z4 symmetry and an unbroken U(1)_L. The imaginary component of rho, chi, is a stable scalar dark matter candidate protected by a C_dark/CP symmetry. Spontaneous Z4 breaking generates a Dirac neutrino mass through the induced VEV of eta, and also produces domain walls. A small bias term V_bias = -sqrt(2) epsilon v_rho^4 cos theta (Eq. 3.6) is introduced to make the walls annihilate and emit gravitational waves within reach of LISA, BBO, mu-ARES, and other experiments. The paper studies DM relic density, direct and indirect detection, Dirac leptogenesis, Delta N_eff, and presents numerical scans showing parameter regions consistent with neutrino oscillation data, baryon asymmetry, and DM constraints, together with projected GW sensitivities.","tokens_in":36779,"tokens_out":8787,"duration_ms":79900,"significance":"If valid, the model would be a minimal combined solution to neutrino mass, dark matter, and baryogenesis, with correlated observational signals in gravitational waves, Delta N_eff, neutrinoless double-beta decay, and DM searches. The economy of using the seesaw field content also for DM, stabilized by C_dark, is a genuine virtue. The numerical work uses standard public tools (SARAH, SPheno, micrOMEGAs) with clearly stated scan ranges and benchmark points, making the analysis reproducible in principle. The GW spectra are compared with a broad set of current and future experiments, and the Delta N_eff discussion is quantitative. However, the central GW claim depends on the fate of the zero-bias v1-v3 domain walls, which is not established, and the predictive content of the GW reach is weakened because the bias epsilon is a free parameter scanned over many orders of magnitude.","major_comments":[{"comment":"The paper excludes the domain wall between minima v1 and v3 on the grounds that delta V13 = 0 and that its collapse dynamics is more complex, but no dynamical proof or simulation is provided. Since the bias term in Eq. (3.6) is even under rho -> rho*, the minima v1 and v3 remain exactly degenerate, so a wall separating them has no volume pressure to drive annihilation. If such walls form a percolating network, they would survive and could dominate the universe, invalidating the BBN and Delta N_eff constraints and the GW signal calculation. The authors must either prove by a lattice simulation or a controlled analytic argument that these walls are destabilized by the surrounding network, introduce a small C_dark-breaking bias that lifts delta V13 and show that the DM lifetime constraint (Appendix D) is still satisfied, or demonstrate that v1/v3 regions do not percolate for the relevant initial conditions.","section":"Section 3, Eq. (3.7) and the paragraph following Eq. (3.12)"},{"comment":"The GW spectrum is computed using the instantaneous-annihilation approximation and formulas calibrated for Z2 domain walls from Ref. [165], while the Z4 network has two distinct wall types with different tensions and biases. The paper acknowledges that delayed annihilation shifts the peak frequency and modifies the spectral shape (Refs. [187-189]) but does not incorporate these effects. Since the peak amplitude scales as sigma^4/delta V^2 and the detectability in Fig. 13 is quantified by SNR > 10, the projected reach could shift substantially; a quantitative estimate of this systematic uncertainty is needed before the GW claims can be considered robust.","section":"Section 3, Eqs. (3.10)-(3.12) and Fig. 8"},{"comment":"The bias parameter epsilon is treated as a free parameter scanned in the range 10^-26 to 10^-21, and the colored GW-sensitive region is the locus of SNR > 10 for that scan. The QG discussion in Section 3 and Appendix D only establishes that epsilon values of this order can be compatible with DM lifetime; it does not predict epsilon. The abstract's statement that the model generates GWs 'within reach' of future experiments is therefore a conditional projection, not a parameter-free prediction. This should be stated explicitly in the conclusions, and the claim of verifiability should be tempered accordingly.","section":"Section 5, Fig. 13"},{"comment":"The baryon asymmetry is computed with the simplified efficiency kappa_f = 1 for K <= 1 and kappa_f = 0.12/K^1.1 for K > 1. For weak washout with thermal initial abundance, kappa_f is not generally unity but depends on the initial conditions and on washout processes; the adopted approximation can overestimate eta_B. Since the lower bound M1 >~ 1.6 x 10^9 (v_eta/100 GeV) GeV and the colored regions in Fig. 13 depend on this efficiency, a more complete solution of the Boltzmann equations is needed to confirm the claimed parameter space for Dirac leptogenesis.","section":"Section 2.2, Eqs. (2.23)-(2.24)"}],"minor_comments":[{"comment":"The manuscript contains numerous typos and formatting errors, including 'Additionnaly', 'magnituded', 'mestasble', 'T able 1', 'a-as-well-as', and inconsistent table/figure captions; a careful proofread is needed.","section":"Throughout"},{"comment":"The statement that 'approximately 4/5th of the domain walls evolve under the adjacent bias and 1/5th under the non-adjacent bias' is not derived from the actual population of wall types in the Z4 network; weighting by the number of bias types rather than by the produced wall fractions is an oversimplification that affects the total GW spectrum in Fig. 8.","section":"Section 3, text after Eq. (3.12)"},{"comment":"The scan imposes the resonance condition m_chi = m_h3/2, so the plotted DM mass ranges are projections onto the resonance, not independent predictions of the model; this should be clarified when summarizing the 'predicted' m_chi ranges.","section":"Section 5"},{"comment":"The caption does not specify whether the y-axis is epsilon itself or epsilon times v_rho^4; the axis label should be defined explicitly.","section":"Fig. 13 caption"}],"recommendation":"major_revision","confidential_remarks":"The zero-bias v1-v3 domain wall issue is likely to be the main point of contention from other referees and should be made a mandatory condition for acceptance. The paper's presentation of the GW signal as a prediction of the model, rather than as a consequence of choosing epsilon, should also be revised. The leptogenesis efficiency approximation is standard in the literature but should be checked or at least its limitations acknowledged."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper adds a Z4-symmetric Type-I Dirac seesaw where the same complex scalar ρ that breaks Z4 gives a stable scalar DM through the unbroken dark-CP symmetry. That combination—and the two-wall-type Z4 domain-wall spectra in Figs. 8-9—is a genuine extension of the earlier Z2 treatments. The neutrino-mass, relic-density, and direct-detection scans are done with standard packages and are internally consistent. The authors also flag correctly that some \"predictions\" are projections onto chosen parameters rather than independent derivations.\n\nThe problem is in the domain-wall section. The bias term V_bias = -√2 ε vρ^4 cosθ is CP-even, so the minima v1 and v3 remain exactly degenerate. The paper simply says they are excluded \"since the associated bias δV13 = 0,\" without a dynamical argument. That's the load-bearing step: if those walls form in the network, the bias gives them no pressure to collapse, and a percolating v1-v3 network would survive and dominate before BBN. The QG operators in Appendix D are also CP-even, so they don't break this degeneracy. The stress-test concern is correct, and it's not a minor gap.\n\nThe instantaneous-annihilation formula and the tuned ε values are additional soft spots, but the paper acknowledges those approximations. The ΔNeff statement is conditional on a thermalization band, which the scan supports.\n\nThe paper is worth a serious referee: a concrete model with a clean, identifiable flaw is exactly what referees are for. But the advertised result—that the walls all annihilate and produce LISA/BBO signals—is not yet established. A revision should add an imaginary bias component to lift the v1-v3 degeneracy or provide a lattice or analytic proof that those walls do not survive, and should publish the scan inputs. I'd bring it to a reading group to discuss the wall question, but I wouldn't cite the GW claim in my own work as it stands.","headline":"Viable Dirac seesaw with DM, but the v1-v3 domain-wall exclusion is a load-bearing gap that leaves the advertised GW signal unsupported.","tokens_in":37371,"tokens_out":4785,"would_cite":false,"duration_ms":42695,"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 minimal Z4-symmetric Type-I Dirac seesaw can simultaneously supply a stable thermal dark-matter candidate, generate light Dirac neutrino masses, produce the baryon asymmetry through Dirac leptogenesis, and leave gravitational-wave and…","keywords":["Dirac seesaw","scalar dark matter","Z4 symmetry","domain walls","gravitational waves","Dirac leptogenesis","effective neutrino number Neff","dark charge conjugation"],"falsifier":"A lattice simulation of the Z4-breaking phase transition with the bias term of Eq. (3.6) would settle the claim: if the degenerate v1–v3 wall channel percolates and fails to annihilate, the cosmological history is broken. A null result from planned gravitational-wave observatories in the predicted peak band for the scanned benchmarks vρ ∈ [$10^{5}$, 2×$10^{8}$] GeV and ε ∈ [$10^{-26}$, $10^{-21}$] would exclude the central parameter region.","tokens_in":36183,"feed_emoji":"🌊","tokens_out":8661,"duration_ms":71228,"temperature":0.7,"pith_summary":"This paper tries to establish that one minimal field content—the fields already required for the Type-I Dirac seesaw, arranged under a Z4 symmetry and with lepton number conserved by a global U(1)_L—can explain neutrino mass, dark matter, and the baryon asymmetry together. The imaginary part of the complex scalar ρ is a stable thermal scalar dark-matter candidate protected by an effective dark charge-conjugation symmetry, so no extra dark-sector fields are needed. The same spontaneous Z4 breaking that generates tiny Dirac neutrino masses also produces domain walls, and a small explicit bias term V_bias = -√2 ε $v_ρ^{4}$ cos θ makes those walls annihilate, producing a stochastic gravitational-wave background within reach of planned detectors. The parameter space that fits neutrino oscillations and the observed baryon asymmetry via Dirac leptogenesis also predicts enhanced ΔNeff and resonant dark-matter masses that future dark-matter and CMB experiments can probe.","feed_headline":"One seesaw yields dark matter, neutrino mass, and gravitational waves","feed_subtitle":"Stable scalar dark matter plus gravitational waves and ΔNeff in reach of planned detectors.","key_machinery":"The central mechanism is the Z4-symmetric scalar potential with an induced VEV hierarchy, where the coupling λ1 both fixes the dark-matter mass through the resonance condition mχ = mh3/2 and controls the domain-wall profile. The imaginary component χ of ρ is the stable scalar dark matter, and the four degenerate minima of the potential give adjacent and non-adjacent walls with different tensions, σ_adj = mθ $v_ρ^{2}$/2 and σ_non-adj = (2√2/3)√λρ $v_ρ^{3}$. The load-bearing bias term V_bias = -√2 ε $v_ρ^{4}$ cos θ lifts the degeneracy so walls annihilate when vacuum pressure equals tension; its magnitude, together with the wall tensions, sets the gravitational-wave spectrum via the peak formulas of Eqs. (3.10)-(3.11). The same parameter set fixes v_η, the seesaw scale M_1, and the Yukawa couplings Y_L and Y_R that determine the neutrino mass, Dirac leptogenesis, and the right-handed-neutrino decoupling temperature controlling ΔNeff.","core_discovery":"In this model, the neutrino mass formula is mν = Y_L $M_N^{{-1}}$ Y_R v v_η / 2, with the η VEV vη induced by the ρ VEV vρ, while the dark-matter candidate χ is the imaginary component of ρ with mass $mχ^{2}$ ≃ 8λ1 $v_ρ^{2}$ + 2√2 μ1 v_η. The authors show that after spontaneous Z4-breaking there are four degenerate minima and two classes of domain walls; the bias in Eq. (3.6) gives adjacent walls a potential difference √2 ε $v_ρ^{4}$ and non-adjacent walls 2√2 ε $v_ρ^{4}$, making the walls annihilate before BBN and emit gravitational waves whose peak amplitude and frequency are set by wall tensions and δV. A numerical scan restricted to the resonant regime mχ = mh3/2 finds points that simultaneously satisfy relic density, direct and indirect detection bounds, electroweak precision data, and vacuum stability up to the Planck scale, while reproducing neutrino oscillation data and a baryon asymmetry ηB in the observed range. The same scan places the additional relativistic degrees of freedom from right-handed neutrinos at ΔNeff = 0.14 in the thermalized region, and predicts GW signals with SNR > 10 for planned observatories across a range of dark-matter masses and seesaw scales.","pith_inferences":["An extension the paper leaves implicit is deriving the bias parameter ε from explicit Planck-scale-suppressed operators and checking that the same operators do not make the dark matter decay faster than the age of the Universe.","The v1–v3 wall channel is excluded rather than simulated; a dedicated lattice study would either strengthen the model by showing that channel annihilates through neighbouring-wall pressure or overturn the gravitational-wave prediction, so it is the clearest next step.","The reported ΔNeff value is 0.14 when right-handed neutrinos thermalize and decouple; coupling the decoupling temperature to leptogenesis efficiency could produce a continuum of smaller values, giving a testable correlation that the paper does not explore.","Because the gravitational-wave peak frequency and the dark-matter mass are linked by the resonance condition, a future GW measurement plus a direct scalar mass measurement would overdetermine vρ and λ1, offering a multi-messenger discriminator among models of this type."],"forward_implications":["A detection of the predicted gravitational-wave peak would directly measure the Z4-breaking scale vρ and, through the resonance condition, the dark-matter mass mχ.","The same parameter space ties the seesaw and leptogenesis scale to vη via M1 ≳ 10^9 GeV (vη/100 GeV), so the mechanism is not pushed to arbitrarily high scales.","Regions with thermalized right-handed neutrinos predict ΔNeff = 0.14, which future CMB surveys will be able to confirm or exclude.","In the resonant regime, dark-matter masses range from about 10^2 to 10^5 GeV, with direct-detection lower bounds between 161 and 593 GeV for the sampled portal couplings; these are testable by next-generation dark-matter experiments.","A future positive signal in neutrinoless double beta decay would falsify the Dirac nature of neutrinos assumed here and break the model's connection to the observed baryon asymmetry."],"supporting_citations":[{"why":"Establishes the earlier Z2-symmetric Dirac seesaw with domain-wall gravitational waves that this paper extends to the Z4 case and to the scalar dark-matter candidate.","marker":"[80]"},{"why":"Supplies the complex-scalar dark-matter mechanism, where the imaginary component is stabilized by CP or dark charge conjugation, and the associated domain-wall collapse analysis.","marker":"[103]"},{"why":"Provides the adjacent-wall profile and tension formula σ_adj = mθ vρ²/2 used for the Z4 walls.","marker":"[146]"},{"why":"Supplies the area parameter A ≈ 1.5 and simulation-based estimates for gravitational-wave spectra from long-lived domain walls.","marker":"[148]"},{"why":"Gives the 4/5 adjacent and 1/5 non-adjacent wall population fractions and the pV = pT annihilation condition used for the weighted gravitational-wave spectrum.","marker":"[153]"},{"why":"Provides the peak amplitude and peak frequency formulas for gravitational waves from collapsing domain walls used in Eqs. (3.10)-(3.11).","marker":"[165]"},{"why":"Supplies the Dirac leptogenesis CP asymmetry and Boltzmann-equation formalism used to compute the baryon asymmetry.","marker":"[111]"},{"why":"Supplies the ΔNeff = 0.14 contribution for three thermalized right-handed neutrinos decoupling before the electroweak scale.","marker":"[59]"}],"fun_headline_variants":["One seesaw, three cosmic gifts: DM, neutrino mass, and GWs","Minimal seesaw solves DM, neutrino mass, and GWs","Seesaw triple play: dark matter, neutrinos, gravitational waves","Seesaw yields dark matter, neutrinos, and spacetime ripples","Single seesaw seeds DM, neutrino mass, and GW signals"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the one domain-wall channel that the bias does not split—the wall between the two opposite minima that remain exactly degenerate—can be safely ignored and will not survive as a percolating network that dominates the universe.","fun_headline_variants_meta":{"raw":{"variants":["One seesaw, three cosmic gifts: DM, neutrino mass, and GWs","Minimal seesaw solves DM, neutrino mass, and GWs","Seesaw triple play: dark matter, neutrinos, gravitational waves","Seesaw yields dark matter, neutrinos, and spacetime ripples","Single seesaw seeds DM, neutrino mass, and GW signals"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001709,"raw_usage":{"total_tokens":6807,"prompt_tokens":1031,"completion_tokens":5776,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":647,"completion_tokens_details":{"reasoning_tokens":5681}},"tokens_in":647,"tokens_out":5776,"duration_ms":35209,"temperature":1.0,"reasoning_tokens":5681,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T14:16:05.822120+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A lattice simulation of the Z4-breaking phase transition with the bias term of Eq. (3.6) would settle the claim: if the degenerate v1–v3 wall channel percolates and fails to annihilate, the cosmological history is broken. A null result from planned gravitational-wave observatories in the predicted peak band for the scanned benchmarks vρ ∈ [$10^{5}$, 2×$10^{8}$] GeV and ε ∈ [$10^{-26}$, $10^{-21}$] would exclude the central parameter region.","supporting_citations":[],"review_version":1}