{"id":"adad4cac-e65e-4821-8ffc-a624898114e3","arxiv_id":"2608.12483","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"New singlet and doublet fermions, a scalar singlet, and right-handed neutrinos can simultaneously explain Dirac neutrino masses, dark matter, and a two-step inverse electroweak phase transition with observable gravitational waves.","lead":"This paper proposes a Standard Model extension with new fermions and a scalar that can explain dark matter, tiny neutrino masses, and a two-step early-universe phase transition that emits gravitational waves. It maps the parameter space that survives current constraints and identifies which future experiments could detect the predicted signals.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The daisy-resummed effective potential omits the large singlet-doublet fermion contributions to the Higgs and Goldstone thermal masses, which for y2 ~ 3 dominate the SM terms and may eliminate the two-step inverse FOPT.","rationale":"The reader's verdict is CONDITIONAL, focused on the large Yukawa couplings and the Landau pole. My pass identifies a more specific, internal issue: the daisy-resummed potential omits the singlet-doublet fermion contributions to the Higgs and Goldstone Debye masses. For y2 ~ 2.8–3.1, these contributions are several times the SM thermal masses, so omitting them is not a negligible approximation. The Arnold–Espinosa procedure requires including all species that couple to the bosons; the authors themselves include the SD fermions in the CW and thermal potentials (Eqs. 4.6–4.10), so their absence in Eqs. (4.12)–(4.15) is inconsistent. The concrete test would recompute the benchmark transitions with the full Debye masses; if the barrier disappears, the central GW prediction fails. I do not claim the model is impossible, but the presented evidence for the two-step FOPT is not yet reliable. The paper has other strengths: the radiative seesaw and DM relic calculation are standard, and the authors are transparent about the Landau pole. Still, the phase-transition result needs this correction. The verdict should remain conditional, with the added condition that the daisy resummation be corrected.","tokens_in":23675,"tokens_out":15378,"duration_ms":136088,"concrete_test":"Recompute Veff(h,T) for BP1–BP3 with the daisy-resummed thermal masses augmented by the singlet-doublet fermion self-energy: add to Πh,g(T) in Eq. (4.12) the term (T^2/24)(d^2(m_X^2+m_Y^2)/dh^2|_0) ≈ (T^2/12) y_2^2 per generation (or the exact finite-T expression). If, with this term, the benchmark points no longer exhibit two first-order transitions at Tc1,Tc2 and the GW spectra of Fig. 12 are lost, the central claim is invalidated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result—the two inverse first-order transitions and the GW spectra (Figs. 11–12, Tables 2–3)—rests on the finite-temperature effective potential Veff in Eq. (4.16). In the daisy resummation, Eqs. (4.11)–(4.15), the thermal masses Πh,g, ΠW,Z,γ are taken from the SM alone. The singlet-doublet fermions with large Yukawa couplings y2 ≈ 2.8–3.1 also couple to the Higgs doublet through y_i Ψ̄_i H̃ χ_i. Their contribution to the Higgs/Goldstone thermal self-energy is of order (Σ_i y_i^2/3) T^2 in the high-T limit; for the benchmark points this is ≳ 2.5 T^2, several times larger than the SM value ≈ 0.4 T^2. Eq. (4.12) contains no such term. Because the daisy correction enters as (m_i^2(h)+Π_i)^{3/2}, omitting a dominant thermal-mass contribution can change the shape and minima of Veff and, plausibly, remove the barrier needed for the first-order transitions. The same omission affects the Goldstones (with identical Yukawa couplings) and therefore the ring-resummed potential. This is an internal inconsistency in the very calculation that produces the two-step FOPT, independent of the separate concern about the Landau pole at ~26 TeV noted in Appendix A. If the full thermal masses restore the usual high-T symmetric phase without a second minimum, the benchmark GW signals do not arise.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends the singlet-doublet fermion dark matter framework by adding two generations of vector-like singlet-doublet fermions, one complex scalar singlet, and three right-handed neutrinos, with the aim of simultaneously explaining radiative Dirac neutrino masses, the observed dark matter relic density, and an inverse two-step first-order electroweak phase transition. The authors reconstruct the neutrino-portal Yukawa couplings from oscillation data via a Casas-Ibarra style parameterization, impose dark matter relic and direct-detection constraints, and then compute the finite-temperature effective potential with singlet-doublet fermions to obtain three benchmark points with two first-order phase transitions and gravitational-wave spectra that they project to be observable by future experiments. The paper also discusses \\Delta N_eff and collider signatures, and it acknowledges in Appendix A that the large fermion-Higgs couplings hit a Landau pole near 26 TeV.","tokens_in":1620,"tokens_out":1674,"duration_ms":230006,"significance":"If the phase-transition calculation survives scrutiny, the model would be a compact and phenomenologically rich framework connecting dark matter, radiative neutrino mass, and gravitational-wave astronomy, with falsifiable predictions for cLFV, \\Delta N_eff, displaced-vertex searches, and direct detection. The paper uses standard tools (micrOMEGAs, FindBounce, SARAH) and the neutrino Yukawa reconstruction makes the parameter scan predictive rather than ad hoc. However, the central FOPT and GW claims rest on a one-loop effective potential with a load-bearing technical gap: the daisy-resummed Debye masses omit the singlet-doublet fermion contributions, and the mass-eigenvalue formulas contain an internal factor inconsistency. These issues must be fixed before the benchmark results can be considered reliable.","major_comments":[{"comment":"The daisy-resummed thermal masses for the Higgs and Goldstone bosons contain only SM contributions, yet the singlet-doublet fermions with y2 \\approx 2.8-3.1 couple to the Higgs with O(1) Yukawa couplings and contribute a leading high-T term to Pi_{h,g} of order y2^2 T^2 times an O(1) group-theory coefficient. For the benchmark values this contribution is at least comparable to, and plausibly larger than, the entire SM contribution in Eq. (4.12), and it is simply absent from the calculation. Because the daisy correction enters as (m_i^2 + Pi_i)^{3/2}, omitting this term can change the shape of V_eff, remove or create the barrier, and therefore alter the two-step inverse-FOPT structure in Figs. 11-12 and the GW spectra in Fig. 12. The authors should include the singlet-doublet fermion contributions to the Debye masses, or demonstrate numerically that the benchmark transitions survive their inclusion.","section":"Sec. 4, Eqs. (4.11)-(4.15)"},{"comment":"The mass-matrix formulas are mutually inconsistent. From the interaction -y_i \\bar Psi_i \\tilde H chi_i, the off-diagonal entry after EWSB is y_i v / sqrt(2), which gives eigenvalues containing sqrt((M_{Psi_i}-M_{chi_i})^2 + 2 y_i^2 v^2) and tan(2 theta_i) = sqrt(2) y_i v / (M_{Psi_i}-M_{chi_i}). This matches Eq. (2.7), but Eqs. (2.8) and (4.6) instead contain 4 y_i^2 h^2 inside the square root, a factor-of-two error. Since the field-dependent masses m_X(h) and m_Y(h) enter V_CW and V_T and hence V_eff in Eq. (4.16), the FOPT results and the benchmark points in Tables 2-3 can shift. The authors must correct the factor or explicitly clarify the definition of y_i used in the numerical calculation.","section":"Sec. 2 and Sec. 4, Eqs. (2.7), (2.8), (4.6)"},{"comment":"The paper acknowledges that the singlet-doublet Yukawa coupling exceeds the perturbative limit at about 4.85 TeV and develops a Landau pole near 26 TeV. Although the phase transition occurs at temperatures below 300 GeV, the effective potential is computed at one-loop order with y2 \\approx 3, so higher-order fermionic corrections to V_eff and to the thermal self-energies are not parametrically small. The central claim of an inverse FOPT with detectable GW signals should be accompanied by a quantitative estimate of the sensitivity of the barrier and of the transition parameters to the truncation of the perturbative expansion, for example by varying the renormalization scale or by estimating the size of two-loop contributions.","section":"Appendix A and Sec. 7"}],"minor_comments":[{"comment":"The sums are written over i=1 to 3, but there are only two generations of singlet-doublet fermions; the sums should run over i=1,2 unless a third generation is intended.","section":"Eqs. (2.18) and (2.19)"},{"comment":"The quantities M_chi, Delta M_1, M_chi', and M_psi' should be explicitly labeled as physical mass-eigenstate masses or as Lagrangian masses; the text appears to use both notions, and the distinction matters for Eqs. (2.7)-(2.8).","section":"Table 2 and related text"},{"comment":"The DarkSide reference is rendered as 'Darkside binom.' and the DARWIN reference as 'DAR WIN'; these should be corrected to 'DarkSide-50' and 'DARWIN'.","section":"References [84] and [85]"},{"comment":"The legend entries for the benchmark points are garbled; they should read BP1, BP2, BP3.","section":"Fig. 12"},{"comment":"There are several typographical errors, including 'natural occurance' in the introduction and 'equilibrium comiving densities' in Section 3; these should be corrected in a final proofreading.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper is well within the scope of the journal and the novelty is not in question. The main risk is technical: the FOPT benchmark results depend on a one-loop effective potential whose daisy resummation omits the dominant new-fermion thermal-mass contributions, and the mass formulas contain an internal factor inconsistency. I would ask the authors to rerun the phase-transition and GW analysis with the corrected thermal masses and formulas before publication. The absence of any shipped code makes independent verification more difficult, though this is not by itself disqualifying."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper's core new claim—two-step inverse first-order EWPT from singlet-doublet fermions with y2~3, with double-peaked GW and correlated Delta Neff—rests on an effective potential that omits the new fermion contributions to the Higgs/Goldstone thermal masses. That looks like a load-bearing error. If I'm right, the benchmarks in Tables 2–3 and Figs. 11–12 are not reliable as they stand.\n\nWhat's genuinely new: earlier work (refs [35,56,58]) did the same dark sector and Dirac seesaw, and refs [66,68] did fermion-driven inverse FOPT, but the joint analysis with two FOPTs, GW spectra, Delta Neff, and collider constraints is new. The model is clearly laid out, the neutrino mass parametrization via Casas-Ibarra is standard, relic density and direct detection are handled with micrOMEGAs, and the authors are upfront about the Landau pole at ~26 TeV in Appendix A. That honesty is real.\n\nThe problem: Eq. (4.12) for Pi_{h,g} contains only SM contributions—Higgs, Z, W, top—even though the heavy singlet-doublet fermions carry Yukawa couplings y2~2.8–3.1 to the same Higgs doublet. Their thermal self-energy contribution scales like y^2 T^2, which for these benchmark points is several times larger than the SM terms. In the Arnold-Espinosa daisy resummation, leaving that out can shift the minima and plausibly remove the barrier that the two FOPTs rely on. You cannot compensate with V_T and V_CW; the ring-resummed zero-mode part is missing a term of the same parametric size. This is independent of the Landau pole worry. If the omitted terms restore the usual high-T symmetric phase, the GW signals disappear.\n\nOther soft spots are secondary: no code or data shipped, so the micrOMEGAs and FindBounce numbers are not independently reproducible; and the O(1) Yukawas with a Landau pole at ~26 TeV mean the effective potential itself is on shaky ground without a UV completion. But the missing daisy term is the one that decides the verdict.\n\nWho it's for: model-builders working on radiative seesaw plus GW phenomenology will want to know about this combination, but they should treat the phase-transition predictions as needing a corrected calculation. It deserves a serious referee—the error is a straightforward fix and the question is worth asking—but I'd only accept after the authors recompute with the full thermal masses. My own verdict is conditional, leaning no on the central claim until it's fixed.","headline":"An interesting new combination of single-doublet dark matter, radiative Dirac neutrino mass, and a fermion-driven inverse first-order EWPT, but the two-step FOPT and GW predictions look unreliable because the daisy-resummed potential omits the new fermions' large contributions to the Higgs thermal mass.","tokens_in":24718,"tokens_out":3614,"would_cite":false,"duration_ms":32874,"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":"One extension with singlet-doublet fermions gives dark matter, neutrino mass, and a two-step electroweak phase transition.","keywords":["singlet-doublet fermion dark matter","radiative Dirac seesaw","inverse first-order electroweak phase transition","gravitational waves","neutrino mass","dark matter relic density","dark radiation","collider signatures"],"falsifier":"Compute the finite-temperature effective potential beyond one-loop, or on the lattice, for the benchmark points of Table 2: if the barrier between the symmetric and broken minima vanishes at $y_2 \\simeq 2.8$–$3.1$, the inverse two-step transition is not real. Observationally, a null search for the predicted double-peaked stochastic gravitational-wave background—with the second transition strengths and durations of Table 3 (i.e. $\\alpha_2$ between 0.03 and 0.25 and $\\beta/H$ between 180 and 1358)—would exclude the two-transition claim.","tokens_in":23441,"feed_emoji":"🌌","tokens_out":16472,"duration_ms":128797,"temperature":0.7,"pith_summary":"This paper argues that a single compact extension of the standard model—two generations of singlet-doublet fermions, one complex scalar singlet, and three right-handed neutrinos—can account for the two measured neutrino mass splittings, the thermal dark-matter relic density, and a two-step inverse first-order electroweak phase transition. In that phase history the Universe first tunnels from a broken Higgs phase to a symmetric phase and then, at a lower temperature, tunnels back to the broken phase, giving two first-order transitions instead of the standard crossover. The lightest singlet-doublet admixture, stabilized by the residual $Z_2$ symmetry, is dark matter; the heavier generation couples to the Higgs with order-one Yukawa couplings and drives the transition. The same particle content keeps all existing constraints satisfied and predicts gravitational-wave signals in the band of planned observatories, a dark-radiation shift, and distinct collider signatures.","feed_headline":"One model gives dark matter, neutrino mass, and a two-stage phase jump","feed_subtitle":"The same fermions that set the dark-matter relic also drive two gravitational-wave-producing phase transitions.","key_machinery":"The carrying mechanism is the singlet-doublet fermion pair $(\\chi_i,\\Psi_i)$ and its Yukawa coupling $y_i\\,\\overline{\\Psi_i}\\tilde{H}\\chi_i$ to the Higgs. After electroweak symmetry breaking this pair mixes by an angle $\\theta_i$, producing the light dark-matter state $\\chi$ and a heavier state $\\psi$; the large coupling $y_2$ feeds the field-dependent fermion masses in the one-loop finite-temperature effective potential, and the competition with the standard-model thermal potential creates the barrier that makes the two-step inverse transition possible. The same two generations run in the one-loop radiative Dirac seesaw, where the soft $Z_4$-breaking term $\\mu_\\phi^2$ splits the singlet scalar into $\\phi_1,\\phi_2$ and generates the neutrino mass: each generation contributes the factor $(M_Y-M_X)\\sin 2\\theta_i$ times a loop function. Thus one portal, the SD fermion-Higgs coupling, drives the phase transition, while the residual $Z_2$ supplies dark-matter stability and the loop seesaw supplies neutrino mass.","core_discovery":"The central claim is that two vector-like singlet-doublet fermion generations, a complex scalar singlet with a soft $Z_4$-breaking mass term, and three Dirac right-handed neutrinos solve three problems at once. The neutral state $\\chi$, a mix of singlet and doublet after electroweak symmetry breaking, provides the observed dark matter relic density; the one-loop diagram with the split scalar singlet and the two fermion generations produces two non-zero Dirac neutrino masses consistent with oscillation data; and the heavier singlet-doublet generation, with Yukawa couplings $y_2 \\simeq 2.8$–$3.1$, turns the electroweak transition into an inverse first-order transition with two first-order phase transitions (at $T_{c1} \\simeq 215$–$283$ GeV and $T_{c2} \\simeq 80$–$89$ GeV for the benchmark points). The resulting gravitational-wave spectrum has two peaks with distinct frequencies and strengths, within reach of future detectors. The paper thereby demonstrates that dark matter, neutrino mass, and a nonstandard electroweak thermal history can share one minimal origin.","pith_inferences":["A consequence left implicit is that during the intermediate symmetric phase the electroweak gauge symmetry is restored, so any primordial magnetic fields or baryon-number-violating processes active before the second transition would be reprocessed by the thermal plasma; this is a testable cosmological stamp of the scenario.","The perturbative-control issue noted in the paper (the singlet-doublet Yukawa coupling exceeds the perturbative limit around 4.85 TeV and hits a Landau pole near 26 TeV) suggests the benchmark couplings sit near a strong-dynamics threshold; a UV completion at that scale could alter the free-energy barrier and shift the gravitational-wave frequencies, so the GW prediction doubles as a probe of the ","The same two-generation structure could be converted to a Majorana neutrino seesaw by adjusting the symmetry, and the Dirac choice (motivated by dark radiation and the absence of neutrinoless double beta decay) would lose those observational channels while keeping the phase-transition mechanism; comparing the two variants would show how generic the inverse-transition prediction is."],"forward_implications":["If the central claim holds, the same benchmark points that satisfy dark matter and neutrino constraints predict a double-peaked stochastic gravitational-wave background, measurable with planned space- and ground-based detectors.","The two-transition thermal history replaces the single electroweak crossover, so any computation of electroweak-scale out-of-equilibrium processes in the early Universe must be redone in a universe that passes through an intermediate symmetric phase.","Because the neutrinos are Dirac, neutrinoless double beta decay is absent and the right-handed neutrinos contribute an extra radiation component $\\Delta N_{\\rm eff}$; part of the otherwise allowed dark-matter parameter space is already excluded by current CMB bounds and more will be probed by future CMB experiments.","The heavier charged doublet fermion either decays promptly to multilepton final states (parts now excluded by prompt-decay searches) or is long-lived, producing displaced vertices testable at colliders and at future long-lived-particle detectors; if no such signatures appear, the preferred mixing region shrinks.","When the heavier generation mass exceeds about 2215 GeV with perturbative couplings, the inverse transition no longer occurs, so the phase-transition requirement directly bounds the mass spectrum of the second generation."],"supporting_citations":[{"why":"Establishes the radiative seesaw template in which the loop particles are $Z_2$-odd and the lightest one is dark matter, the structure this model builds on.","marker":"[24]"},{"why":"Provides the singlet-doublet Dirac radiative seesaw with $(g-2)_\\mu$ and $\\Delta N_{\\rm eff}$ that this paper extends to include the inverse phase transition.","marker":"[35]"},{"why":"Earlier gravitational-wave study of singlet-doublet dark matter with radiative neutrino mass; the comparison case in which singlet scalars, not fermions, drive the transition.","marker":"[58]"},{"why":"Shows that new TeV-scale fermions can generate multistep strongly first-order phase transitions, the fermion-driven mechanism adopted here.","marker":"[66]"},{"why":"Connects the inverse seesaw with the inverse electroweak phase transition, the two-step transition path used in this model.","marker":"[68]"},{"why":"Supplies the two-sector Boltzmann equations and numerical solver used to compute the dark-matter relic density including coannihilation and conversion processes.","marker":"[74]"},{"why":"Provides the zero-temperature one-loop radiative correction to the scalar potential used in the effective potential.","marker":"[87]"},{"why":"Supplies the high-temperature resummation of bosonic zero modes used for the infrared-corrected effective potential.","marker":"[92]"},{"why":"Computes the Euclidean bounce action $S_3$ used to set nucleation temperatures and gravitational-wave parameters.","marker":"[109]"}],"fun_headline_variants":["One fermion sector: dark matter, neutrino mass, two phase jumps","Singlet-doublet fermions: DM, Dirac neutrinos, double EWPT","One model solves DM, neutrino mass, and a two-step phase change","Triple fix: dark matter, neutrino mass, twin phase transitions","Two fermion generations: DM, neutrino mass, and a double phase jump"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the one-loop finite-temperature effective potential with daisy resummation, evaluated at the benchmark singlet-doublet Yukawa couplings $y_2 \\simeq 2.8$–$3.1$, correctly gives the free-energy barrier and minimum ordering; if higher-order or non-perturbative effects destroy that barrier, the inverse first-order transition and its gravitational-wave spectrum disappear.","fun_headline_variants_meta":{"raw":{"variants":["One fermion sector: dark matter, neutrino mass, two phase jumps","Singlet-doublet fermions: DM, Dirac neutrinos, double EWPT","One model solves DM, neutrino mass, and a two-step phase change","Triple fix: dark matter, neutrino mass, twin phase transitions","Two fermion generations: DM, neutrino mass, and a double phase jump"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000219,"raw_usage":{"total_tokens":1495,"prompt_tokens":1052,"completion_tokens":443,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":668,"completion_tokens_details":{"reasoning_tokens":344}},"tokens_in":668,"tokens_out":443,"duration_ms":4147,"temperature":1.0,"reasoning_tokens":344,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:08:06.618473+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the finite-temperature effective potential beyond one-loop, or on the lattice, for the benchmark points of Table 2: if the barrier between the symmetric and broken minima vanishes at $y_2 \\simeq 2.8$–$3.1$, the inverse two-step transition is not real. Observationally, a null search for the predicted double-peaked stochastic gravitational-wave background—with the second transition strengths and durations of Table 3 (i.e. $\\alpha_2$ between 0.03 and 0.25 and $\\beta/H$ between 180 and 1358)—would exclude the two-transition claim.","supporting_citations":[],"review_version":1}