{"id":"11e239d9-094d-47c8-b19a-88096a6f8265","arxiv_id":"2412.19094","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A KSVZ-like model with a fermion dark matter candidate produced by UV freeze-in through the axion portal is shown to have viable parameter space at TeV mass and axion decay constant near 10^10 to 10^11 GeV.","lead":"Dark matter in this model is a new fermion produced out of equilibrium through the axion portal, with the axion scale set by Peccei-Quinn symmetry breaking. The same symmetry gives tiny masses to neutrinos and solves the strong CP problem, and the paper maps out where the fermion alone can be all of the dark matter.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Isocurvature bound is satisfied only by an unjustified fI = 10^17 GeV; with the model's physical fI = fa ~ 10^10 GeV, the claimed parameter space violates Eq. (3.6).","rationale":"The reader's weakest assumption was the exclusion of the post-inflationary PQ-breaking scenario. My stress-test identifies a different, more direct flaw inside the pre-inflationary scenario they do consider: the isocurvature constraint is apparently satisfied only by decoupling fI from fa, which is not possible in the single-PQ-scalar model described in Section 2. This is more load-bearing than the post-inflationary concern because it affects the claimed 'not excluded' statement even after accepting the paper's own cosmological assumption. It is also concrete and checkable with a back-of-the-envelope calculation. I do not think this changes the overall verdict: the paper remains conditionally plausible if the fI issue is resolved (e.g., by lowering H_I or adding a justified mechanism for fI > fa), but the current presentation is not yet self-consistent. The reader's CONDITIONAL verdict therefore stands, and I set verdict_should_be to UNCHANGED. I only partially agree with the reader's weakest_assumption because they located the soft spot in the pre/post-inflation dichotomy, whereas the more serious issue is the unjustified fI value within the chosen pre-inflation trajectory.","tokens_in":15765,"tokens_out":9167,"duration_ms":84617,"concrete_test":"Take the benchmark point mψ = 1 TeV, fa = 10^10 GeV, TRH = 10^8 GeV, θ_i = 1, H_I = 10^14 GeV. Evaluate Eq. (3.6) with f_I = fa rather than 10^17 GeV. If P_a/P_r > 0.04, the benchmark point is excluded by Planck isocurvature data. Then repeat for all points in Fig. 3; if none survive, the paper must either (i) introduce additional PQ-breaking fields or a non-canonical axion kinetic term to justify f_I >> fa, or (ii) adopt a lower H_I (e.g., H_I ≲ 10^9 GeV) and redo the isocurvature check. This algebraic test settles whether the claimed parameter space is actually allowed under the paper's own Eq. (3.6).","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.3 claims the FIMP parameter space is not excluded by isocurvature, using Eq. (3.6) with H_I = 10^14 GeV, θ_i = {0.1, 1}, and f_I = 10^17 GeV. The problem is that this model contains exactly one PQ-breaking complex scalar Φ whose VEV determines the axion decay constant fa ≈ vΦ ≈ 10^10–10^11 GeV. There is no second PQ field or non-minimal kinetic term to make the inflation-era decay constant f_I three orders of magnitude larger than fa. If the physically consistent identification f_I = fa is imposed, the isocurvature ratio becomes (Ωa/0.12)^2 × 4(H_I/2π)^2/(fa θ_i)^2. For fa = 10^10 GeV, H_I = 10^14 GeV, θ_i = 1, this is about 4 × (1.6×10^13/10^10)^2 ≈ 10^7, exceeding the Planck bound P_a/P_r ≤ 0.04 by many orders of magnitude. Even for θ_i = 0.1 it remains huge. Thus the claimed allowed region in Fig. 3 is, as presented, excluded by isocurvature constraints. The paper says only that 'a larger fI (> fa) is important to suppress these fluctuations' and then sets fI by hand, but offers no mechanism consistent with the field content in Table 1. This is an internal consistency problem, not a disagreement with external consensus: the same scalar VEV sets both fa and fI, so they cannot be chosen independently.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a KSVZ-like extension of the Standard Model with Peccei-Quinn (PQ) symmetry, adding a complex scalar Φ, a vector-like heavy quark Q, a second Higgs doublet H2, a Dirac fermion ψ, and three right-handed neutrinos. PQ symmetry generates the QCD axion, gives Dirac masses to neutrinos via the small H2 VEV, and, together with a Z2 symmetry, stabilizes ψ. The dark matter candidate ψ is produced by UV freeze-in through axion-portal interactions. The author solves the coupled Boltzmann equations for ψ and the axion, identifies parameter space in (mψ, fa) around (1–10 TeV, 10^10–10^11 GeV) for TRH = 10^8 GeV that reproduces Ωψh2 ≈ 0.12, and compares with axion and direct-detection bounds. The paper also estimates the non-thermal (misalignment) axion relic and claims the combined two-DM scenario respects isocurvature constraints.","tokens_in":16013,"tokens_out":20157,"duration_ms":187132,"significance":"If the results hold, the model offers an economical simultaneous solution to strong CP, neutrino mass, and dark matter, and it provides a concrete UV-freeze-in setup with analytic cross sections. The paper correctly uses standard Boltzmann equations and lists explicit cross-section expressions in the appendices; it also transparently states its assumption that PQ breaking occurs before/during inflation. However, the cosmological viability of the claimed parameter space is not established because the isocurvature analysis relies on an unjustified inflation-era decay constant f_I, and the quoted isocurvature inequality in Eq. (3.6) is not the standard expression for the axion isocurvature relative to curvature fluctuations. These issues are load-bearing for the central claim that the FIMP parameter space is not excluded.","major_comments":[{"comment":"The isocurvature constraint is not satisfied within the model as written. The axion decay constant is fixed by the single PQ-charged scalar Φ in Eq. (2.12), fa ≈ x_Φ v_Φ, and the FIMP relic-density region of Fig. 3 requires v_Φ ≈ 10^10–10^11 GeV. The paper nevertheless sets f_I = 10^17 GeV in Eq. (3.6) with the comment that 'a larger fI (> fa) is important to suppress these fluctuations,' but no field content or coupling in Table 1 can produce an inflation-era decay constant three orders of magnitude larger than fa. If the physical identification f_I = fa is imposed, the standard isocurvature ratio P_a/P_r = (Ω_a/0.12)^2 (H_I/(π f_a θ_i))^2 / P_r, with Ω_a from Eq. (3.5), exceeds the Planck bound P_a/P_r ≤ 0.04 by many orders of magnitude (for fa = 10^10 GeV, H_I = 10^14 GeV, θ_i = 1, and Ω_a ≈ 5×10^-4, the ratio is O(10^10)). Moreover, Eq. (3.6) as written is not the standard definition of P_a/P_r; the usual expression divides by the curvature power amplitude P_r ≈ 2×10^-9, which makes the discrepancy even larger. The statement that the FIMP parameter space is 'not excluded' by isocurvature is therefore unsupported unless a concrete mechanism producing f_I ≫ fa is supplied.","section":"Section 3.3, Eqs. (3.5)-(3.6); Table 1"}],"minor_comments":[{"comment":"The paper explicitly excludes the post-inflationary PQ-breaking scenario (fa < TRH) because axions could thermalize and affect the FIMP analysis. This is a genuine limitation of the claimed parameter space, and it should be stated more prominently as an assumption limiting the validity of all conclusions to pre-inflationary PQ breaking.","section":"Section 3.3"},{"comment":"The caption lists five curves ('black, green, blue, pink, and red') but only three process categories are named ('DM - axion, DM - gluon, and axion - gluon'). Each curve (Hubble, aa→ψψ, ψψ→gg, aa→gg, GG→Ga) should be identified explicitly by color.","section":"Figure 1"},{"comment":"The numerical solution of the coupled Boltzmann equations is not described beyond stating that they are solved numerically; providing the integration method, precision criteria, and initial conditions would improve reproducibility.","section":"Sections 3.1–3.2"},{"comment":"The FIMP contour in the (mψ, |gaγ|) plane is overlaid on the axion-mass plane in Fig. 4 via Eq. (2.13); this mapping should be stated in the caption so that the reader can follow the conversion.","section":"Figures 3 and 4"},{"comment":"There are several typographical errors: 'isocuravture' should be 'isocurvature', 'Similalry' should be 'Similarly', 'boltzmann' should be 'Boltzmann', and 'T able 1' in the Table 1 caption should be 'Table 1'.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope and the freeze-in calculation is competently presented, but the isocurvature section is a serious gap. The f_I = 10^17 GeV input is an ad hoc assumption with no support in the model, and Eq. (3.6) appears non-standard. I recommend major revision: the author should either introduce a concrete mechanism for f_I > fa (e.g., a second PQ-charged field or an inflationary coupling) or reanalyze the parameter space under f_I = fa, in which case the claimed region is likely excluded. The paper would also benefit from a clearer statement of the scope limitation to pre-inflationary PQ breaking."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version: this is a standard KSVZ-plus-freeze-in model paper, competently done, but the claim that the parameter space is not excluded by isocurvature rests on an input the model does not actually have. The paper sets fI = 10^17 GeV to suppress the axion fluctuations, but there is only one PQ-breaking scalar in Table 1, and its VEV sets both fa and fI. With fI = fa ≈ 10^10–10^11 GeV, Eq. (3.6) gives a ratio roughly seven orders of magnitude above the Planck bound for the stated HI and θi. The text only says 'a larger fI (> fa) is important' and then treats fI as free. That is the load-bearing soft spot.\n\nWhat the paper does well: it assembles the known ingredients (KSVZ axion, Dirac neutrino mass from a PQ-charged Higgs doublet, axion-portal freeze-in of a Z2-stable fermion) into a concrete model, and the freeze-in machinery is standard. The yield curves in Fig. 2 behave as expected, and the allowed band in the (mψ, fa) plane is a reasonable outcome of solving the coupled Boltzmann equations. The direct-detection suppression (q^4 and 1/fa^2) is correctly noted.\n\nThe other soft spots are less severe. The post-inflationary PQ-breaking scenario is discarded in a sentence, although a proper treatment of strings and domain walls could change the axion contribution. There is no code or detailed scan recipe, so the numerics are not independently reproducible. Several parameters (TRH, θi, HI, vH2, xΦ) are fixed by hand without a sensitivity study.\n\nAll that said, the paper is not incoherent. The standard freeze-in calculation is legitimate, and the combination of three problems in one PQ framework is a sensible exercise. If the fI issue is fixed by either presenting a mechanism for fI >> fa (e.g., a second PQ field or a non-minimal kinetic term) or restricting the axion relic so that isocurvature is safely small, the FIMP-only scenario may survive. As written, that central claim needs revision.\n\nThe right reader is someone working on axion-portal freeze-in or combined models of DM, neutrino mass, and the strong CP problem. They will find the parameter scan useful, but they should be warned about the fI assumption. I would send it to a serious referee, but the report should emphasize the isocurvature handling and ask for a response on the physical status of fI. A quick revision could make this paper much more solid.","headline":"A workmanlike KSVZ-plus-freeze-in model paper, but the isocurvature bound is handled by setting fI = 10^17 GeV when the single PQ field in the model gives fI = fa ~ 10^10 GeV.","tokens_in":16746,"tokens_out":5367,"would_cite":false,"duration_ms":153772,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["95.35.+d","14.80.Va","14.60.Pq"],"model":"deepseek-v4-flash","headline":"A Dirac fermion produced by UV freeze-in through the axion portal can supply the full dark matter relic density in an unexcluded window of mass and axion scale.","keywords":["dark matter","FIMP","UV freeze-in","QCD axion","Peccei-Quinn symmetry","Dirac neutrino mass","KSVZ model","axion portal"],"falsifier":"A decisive test would be to establish the timing of PQ breaking: if axion dark matter is found with $f_a$ below about $10^{10}$ GeV, or if axion minicluster or string-wall signals show that PQ broke after inflation, the assumed $f_a > T_{\\rm RH}$ regime is wrong and the FIMP-only relic calculation collapses. Alternatively, a measurement of the axion-photon coupling $|g_{a\\gamma}|$ in the model's predicted range around $10^{-14}$ to $10^{-12}$ GeV$^{-1}$ at the corresponding axion masses would test the plotted parameter space directly.","tokens_in":15374,"feed_emoji":"🌌","tokens_out":6937,"duration_ms":69297,"temperature":0.7,"pith_summary":"This paper builds a single extension of the Standard Model that tackles three open problems at once: the strong CP problem, the smallness of neutrino masses, and the particle identity of dark matter. The proposal is a KSVZ-like model in which a Peccei-Quinn symmetry generates the QCD axion, forbids Majorana neutrino masses so neutrinos are Dirac, and opens an axion portal to a $\\mathbb{Z}_2$-stable Dirac fermion. The central result is that this fermion, produced by UV freeze-in, can supply the full observed relic abundance $\\Omega_\\psi h^2 \\approx 0.12$ for masses around 1 to 10 TeV with axion decay constant around $10^{10}$ to $10^{11}$ GeV and reheating temperature $10^8$ GeV, while staying clear of current axion, direct-detection, and isocurvature limits. The paper also shows that the axion itself can contribute through misalignment, making a two-component dark matter scenario possible. If correct, the model would connect three otherwise separate mysteries to one symmetry-breaking scale.","feed_headline":"Axion-portal fermion can be all of dark matter","feed_subtitle":"KSVZ-like model with PQ-protected Dirac neutrinos leaves a viable 1–10 TeV freeze-in window.","key_machinery":"The engine of the model is the axion portal: after PQ breaking, the pseudoscalar axion $a$ emerges from the phase of $\\Phi$ and couples derivatively to the dark fermion $\\psi$, the heavy quark $Q$, and neutrinos with strength set by $1/f_a$. These derivative couplings, combined with $f_a \\gtrsim 10^{10}$ GeV, keep $\\psi$ out of thermal equilibrium with the quark-gluon plasma, while axion production channels such as $G G \\to G a$ can thermalize axions at temperatures above about $10^9$ GeV; the coupled Boltzmann equations then track both yields. The same PQ symmetry forbids Majorana neutrino masses while allowing the Dirac Yukawa coupling $\\bar{\\ell}_L \\tilde{H}_2 \\nu_R$, giving $m_\\nu = y_\\nu v_{H_2}/\\sqrt{2}$ with $v_{H_2} = 10^{-9}$ GeV. Direct detection is suppressed by the $\\gamma_5$ derivative structure of the portal, which produces a $q^4$ momentum suppression on top of the $f_a^{-2}$ factor.","core_discovery":"The central claim is that a KSVZ-like extension containing a Dirac fermion $\\psi$, three right-handed neutrinos, a PQ-charged second Higgs doublet, a complex scalar $\\Phi$, and a vector-like heavy quark can solve the strong CP problem, give small Dirac neutrino masses, and supply dark matter all at once. The paper argues that $\\psi$, stabilized by a $\\mathbb{Z}_2$ symmetry, is produced out of equilibrium through the axion portal via derivative couplings suppressed by $1/f_a$, and that solving the coupled Boltzmann equations for $\\psi$ and the axion yields $\\Omega_\\psi h^2 \\approx 0.12$ for $m_\\psi$ around 1 to 10 TeV and $f_a$ around $10^{10}$ to $10^{11}$ GeV with reheating temperature $T_{\\rm RH} = 10^8$ GeV. It further argues that this FIMP-only solution satisfies current axion-photon, direct-detection, and isocurvature constraints, and that adding misalignment-produced axions with $\\theta_i = 0.1$ or $1$, $H_I = 10^{14}$ GeV, and $f_I = 10^{17}$ GeV admits a two-component dark matter scenario whose total relic abundance also fits the observation.","pith_inferences":["One consequence the paper leaves implicit is that because the ratio $E/N$ vanishes in this model, the axion-photon coupling is suppressed relative to standard KSVZ, making near-future haloscopes and helioscopes less likely to see this axion; a positive axion-photon detection would point away from this specific construction.","If PQ symmetry broke after inflation, axion production from strings and domain walls would likely exceed the FIMP abundance in the 1 to 10 TeV window, so observing axion minicluster or string-wall signatures would disfavor the single-FIMP interpretation.","The paper's neutrino masses are purely Dirac; an observation of neutrinoless double beta decay would require adding a Majorana source, breaking the link the paper draws between PQ symmetry and the absence of Majorana masses."],"forward_implications":["If the claim holds, a single PQ-breaking sector accounts for the strong CP problem, small Dirac neutrino masses, and the observed dark matter abundance without WIMP-scale couplings.","The FIMP mass is predicted to lie near 1 to 10 TeV when $f_a$ is around $10^{10}$ to $10^{11}$ GeV and $T_{\\rm RH} = 10^8$ GeV; lower $f_a$ requires a lighter $\\psi$ and is more tightly constrained by axion searches.","The axion can also be dark matter via misalignment, so the model admits a two-component dark matter sector; for $\\theta_i = 1$ and the chosen inflationary parameters, the axion relic can rival the FIMP relic.","Direct detection experiments such as LUX and XENON1T cannot see this dark matter because the scattering cross section is suppressed by $q^4/f_a^2$, leaving indirect or axion-mediated signatures as the main observational channels."],"supporting_citations":[{"why":"Supplies the Planck relic-density constraint $\\Omega_{\\rm DM} h^2 = 0.12 \\pm 0.001$ that defines the allowed FIMP parameter space.","marker":"[7]"},{"why":"Introduces the FIMP and freeze-in production mechanism that the paper applies to the Dirac fermion.","marker":"[23]"},{"why":"Defines the KSVZ invisible axion model whose field content and PQ charge assignments the paper extends.","marker":"[37, 38]"},{"why":"Establishes the Peccei-Quinn solution to the strong CP problem and the axion as the resulting Nambu-Goldstone boson.","marker":"[34–36]"},{"why":"Provides the thermal axion production rates from the quark-gluon plasma used in the coupled Boltzmann equations.","marker":"[64–66]"},{"why":"Gives the misalignment axion relic-density formula for PQ breaking before or during inflation.","marker":"[84, 85]"},{"why":"Yields the isocurvature bound $P_a/P_r \\le 0.04$ used to choose $f_I = 10^{17}$ GeV and $H_I = 10^{14}$ GeV.","marker":"[86]"},{"why":"Sets the LUX and XENON1T direct-detection limits that the paper's $q^4/f_a^2$-suppressed cross section must evade.","marker":"[15, 17]"}],"fun_headline_variants":["Axion portal freeze-in fermion alone explains dark matter","KSVZ model with Dirac neutrino: freeze-in DM and axion","UV-freeze-in dark fermion via axion portal survives all bounds","Single axion-portal fermion is viable dark matter candidate","1-10 TeV freeze-in fermion dark matter from axion portal"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The parameter space presupposes that PQ symmetry broke before or during inflation, with $f_a > T_{\\rm RH} = 10^8$ GeV, $H_I = 10^{14}$ GeV, and $f_I = 10^{17}$ GeV; the paper explicitly discards the post-inflationary case where strings, domain walls, and thermalized axions could shift the dark matter budget and invalidate the FIMP-only claim.","fun_headline_variants_meta":{"raw":{"variants":["Axion portal freeze-in fermion alone explains dark matter","KSVZ model with Dirac neutrino: freeze-in DM and axion","UV-freeze-in dark fermion via axion portal survives all bounds","Single axion-portal fermion is viable dark matter candidate","1-10 TeV freeze-in fermion dark matter from axion portal"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000384,"raw_usage":{"total_tokens":2054,"prompt_tokens":989,"completion_tokens":1065,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":605,"completion_tokens_details":{"reasoning_tokens":973}},"tokens_in":605,"tokens_out":1065,"duration_ms":9188,"temperature":1.0,"reasoning_tokens":973,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T00:59:31.122336+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test would be to establish the timing of PQ breaking: if axion dark matter is found with $f_a$ below about $10^{10}$ GeV, or if axion minicluster or string-wall signals show that PQ broke after inflation, the assumed $f_a > T_{\\rm RH}$ regime is wrong and the FIMP-only relic calculation collapses. Alternatively, a measurement of the axion-photon coupling $|g_{a\\gamma}|$ in the model's predicted range around $10^{-14}$ to $10^{-12}$ GeV$^{-1}$ at the corresponding axion masses would test the plotted parameter space directly.","supporting_citations":[],"review_version":1}