{"id":"cac6f5e0-5088-4912-beea-970ce3e01f9e","arxiv_id":"2602.21122","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"SM+S has four and SM+2S has eleven inequivalent maximal scalar Lie point symmetry algebras of the field equations, with parameter-based algorithms to identify them.","lead":"This paper classifies all scalar Lie point symmetries of the Standard Model extended by one or two real gauge-singlet scalars, listing the inequivalent symmetry algebras for each model. It also gives parameter-based algorithms for reading the symmetry algebra off any specific potential without solving the determining equations.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The SVS/DS/NVS classification rests on the unverified 'no polynomial consequences' assumption in Theorem 1/Corollary 1; if it fails, the g_var and g_svar lists could change.","rationale":"The paper's central mathematical framework for distinguishing symmetry types is Corollary 1, whose proof depends on the no-polynomial-consequences condition. The reader identified exactly this as the weakest assumption, and my reading agrees: it is load-bearing for the g_var and g_svar classifications, it is explicitly assumed rather than proved, and it is not self-evident from the text because the field equations are never written out and analyzed for algebraic constraints. At the same time, the concern is not that the assumption is false—for the SM+S/SM+2S Lagrangians, with canonical kinetic terms for all scalars, the Euler–Lagrange equations are nondegenerate second-order PDEs, so no polynomial consequence is to be expected. The issue is a gap in rigor: the paper states the hypothesis and then applies it without verifying it. A short argument (or a reference to a standard Cauchy–Kovalevskaya result) would close the gap. Since this is a verifiable technical condition rather than a demonstrated error, the appropriate verdict remains CONDITIONAL: accept after the assumption is proved or explicitly checked for each branch. The g_EL completeness is also computationally dependent on unprovided SYM/Reduce outputs, but that is a reproducibility issue secondary to the logical gap in the theorem's application. I therefore agree with the reader's assessment and recommend no change to the verdict.","tokens_in":34963,"tokens_out":10041,"duration_ms":102749,"concrete_test":"For each reduced branch (SM+S cases and SM+2S Leaves 1–31), verify that the Euler–Lagrange system E(L)=0 has no polynomial consequences. Concretely: (1) Compute the principal symbol of the system; for the scalar equations it is the identity matrix (for SM+S a 5×5 or for SM+2S a 6×6 block with □ for each scalar), and the gauge-field equations are also second-order, so the system is nondegenerate. (2) Invoke the Cauchy–Kovalevskaya theorem for analytic initial data: for any prescribed values (φ(x0), s(x0)) = c of the scalar fields at a point, some local solution of the full system exists (choose the gauge-field initial data to satisfy the Gauss-law constraints, which is always possible locally). Then any polynomial p(φ,s) vanishing on all solutions would vanish at x0 for every c, forcing p≡0. If this check succeeds, the 'no polynomial consequences' assumption is satisfied; a failure would","verdict_should_be":"UNCHANGED","load_bearing_attack":"Theorem 1 and Corollary 1 characterize strict variational, divergence, and non-variational symmetries under the explicit hypothesis that the system E(L)=0 'has no polynomial consequences' (Section 2.4, before Theorem 1). This hypothesis is used to conclude that any polynomial in the scalar fields that vanishes on-shell must be identically zero, yielding prX(V)=a_i α_i as a constant. Sections 3.3 and 4.5 apply this theorem to SM+S and SM+2S without ever verifying the hypothesis. If some branch's field equations implied a nonzero polynomial relation among s1,s2 (or s), then a generator that is actually a divergence symmetry could be misclassified as non-variational, or vice versa, because the key step 'on-shell vanishing implies polynomial identity' would fail. This would affect the central classifications (3.28), (5.3)–(5.5), at least for the variational and strict-variational lists; the g_EL list itself is derived directly from determining equations and may survive. For these physical models the assumption is very plausible—every scalar field has a nondegenerate kinetic term, so the Euler–Lagrange equations are second-order PDEs rather than algebraic constraints—but the paper provides no proof or explicit verification. The reader's conditional verdict correctly targets this gap.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper classifies scalar Lie point symmetries for the Standard Model extended by one or two real gauge-singlet scalars (SM+S and SM+2S), distinguishing strict variational, divergence, and non-variational (Euler–Lagrange) symmetries. For SM+S it finds four inequivalent realizable algebras; for SM+2S eleven. It provides parameter-based algorithms to determine the symmetry algebra of any numerical instance without solving determining equations, and it proves general results (Theorem 1, Proposition 1, Corollary 1, Proposition 2) characterizing symmetry types for Lagrangians with potentials. The classification is obtained by solving determining equations with the Mathematica packages SYM/Reduce over a reduction tree of 31 parameter branches.","tokens_in":35275,"tokens_out":15201,"duration_ms":138209,"significance":"If the results hold, the paper delivers a useful catalogue of scalar Lie point symmetries for two popular BSM models and an efficient algorithm for symmetry identification in parameter scans. The general characterization in Corollary 1 and the explicit proof of Proposition 3 are valuable extensions of earlier work. The g_EL classification is derived directly from the determining equations and is supported by explicit examples, e.g., the Leaf 6 analysis and the calculations around Eqs. (4.113)–(4.115). However, the SVS/DS/NVS classification rests on a structural hypothesis that is stated but never verified for the physical models, which limits confidence in the variational and strict-variational lists.","major_comments":[{"comment":"The theorem and corollary assume that the system E(L)=0 'has no polynomial consequences'. This hypothesis is never verified for SM+S or SM+2S. It is essential for the step concluding that prX(V) is constant from on-shell vanishing of its partial derivatives, and thereby underpins the SVS/DS/NVS classification of every generator. If some branch of the reduction tree admitted a nonzero polynomial relation among the scalar fields, a generator classified as non-variational could actually be a divergence symmetry (or vice versa), changing the g_var and g_svar lists in Eqs. (5.4)–(5.5). The authors should either prove the property for these field equations (e.g., by a local-solvability/Cauchy–Kovalevskaya argument using the nondegenerate kinetic terms) or provide a direct per-leaf computation of prX(L) and E(prX(L)) for each generator in the set ℵ, thereby avoiding the hypothesis altogether.","section":"§2.4 (Theorem 1, Corollary 1) and §§3.3, 4.5"},{"comment":"The exhaustiveness of the computational classification is not fully transparent. The text states that SYM/Reduce solves the determining equations for each of the 31 leaves, but the solving steps and intermediate differential-elimination results are not shown. To support the claimed completeness of the reduction tree and the resulting classification, the authors should make the code available or include a detailed appendix (or supplementary material) with, for each leaf, the reduced determining equations and the final solutions, together with a justification that the branching in Figures 1–3 covers all possible parameter values.","section":"§4.4 and Appendix A"},{"comment":"The summary list contains 'sc_2⊕u(1)_Y' as one of the eleven inequivalent algebras, but Section 4.5.1 shows that sc_1≂O sc_2≂O sc, so sc_2⊕u(1)_Y is not inequivalent to sc⊕u(1)_Y, which is already listed. The intended entry is presumably 'sc_1⊕sc_2⊕u(1)_Y', which appears as one of the three (2+1)d algebras in (4.157). As written, the list contradicts the 'inequivalent' claim and the count of eleven; please correct this.","section":"Eq. (5.3)"}],"minor_comments":[{"comment":"The notation 'sc_2' is ambiguous: it is used for a one-dimensional algebra acting on the second singlet, but in (5.3) it appears to represent a two-dimensional direct sum. Recommend writing 'sc_1⊕sc_2' explicitly when that is intended.","section":"Notation in §5 and §4.5"},{"comment":"The reduction tree is very dense. A legend distinguishing red leaves, intermediate nodes, and '∅' nodes, and perhaps a higher-resolution or simplified version, would improve readability.","section":"Figure 2"},{"comment":"The 'no cancellations' argument between Eqs. (2.55) and (2.56) is somewhat terse. A sentence explaining why, under the assumptions on T, E_i(prX(T)) cannot contain a pure polynomial term in the φ variables would help the reader follow the proof.","section":"§2.4, proof of Proposition 1"},{"comment":"Several standard results are cited to [19], the author's previous work. Since the present paper aims to be self-contained, consider restating the most important definitions (e.g., the equivalence notation ≂O and the affine-vector-field notation) with equation numbers, rather than relying on the earlier reference.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The main barrier is the unverified 'no polynomial consequences' assumption. This is likely fixable—either by a short proof using the nondegenerate kinetic structure or by direct per-leaf checks—but it is load-bearing for the variational and strict-variational classifications. If the authors supply that, the paper could become acceptable. The computational completeness would also be strengthened by sharing code or detailed solving logs. Minor: correct the (5.3) list."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"One thing to know up front: this paper is a solid, useful classification of scalar Lie point symmetries for the Standard Model plus one or two real singlets. It is not a breakthrough, but it is a careful piece of work that should be taken seriously by model builders and by people doing symmetry analysis of scalar field theories.\n\nWhat's new: the SM+2S classification, the parameter-decision algorithm of Section 4.6, the equivalence analysis under affine reparametrizations, and Corollary 1—a general characterization of strict variational, divergence, and non-variational generators for Lagrangians with polynomial potentials. The paper is honest about what it is doing: it separates the Euler-Lagrange algebra from the variational subalgebras, gives explicit parameter conditions for every branch of the reduction tree, and does not fit anything to data. The internal consistency checks, like the explicit computation for the kinetic terms, are well done.\n\nThe soft spots are real. The biggest is the 'no polynomial consequences' assumption in Theorem 1 and Corollary 1. The paper states it, uses it to conclude that any polynomial in the fields that vanishes on-shell must vanish identically, and then applies the theorem to SM+S and SM+2S without verifying the hypothesis. For these models, the assumption is almost certainly true—the scalar fields have nondegenerate kinetic terms, so the Euler-Lagrange equations are second-order PDEs, not algebraic constraints—but the paper does not actually show it. A referee should insist on a proof or an explicit verification. If the assumption failed, the g_var and g_svar lists could change, though the g_EL classification, which comes directly from the determining equations, would likely survive.\n\nThe second issue is reproducibility. The SM+2S classification rests on solving large systems of determining equations with SYM/Mathematica. The parameter conditions are listed, but the actual solving steps are not. This makes the completeness claim hard to audit. It is fixable—supply the notebooks or an appendix with the linear systems—but it is a real gap for a paper whose central claim is completeness.\n\nOn the citation pattern and the literature engagement: everything is in order. The paper builds on the author's earlier 2HDM work, cites the standard references on Lie symmetry analysis and singlet extensions, and is careful to distinguish what is proved here from what is assumed.\n\nWho is this for? People who work on the physics of Higgs-portal models, especially singlet extensions, and anyone doing symmetry classification of scalar field theories. I would not put it in every reading group, but it would be a good topic for a group with a taste for computational Lie symmetry analysis.\n\nRecommendation: send to serious peer review. It should not be desk rejected. The right outcome is 'major revision'—the author needs to close the no-polynomial-consequences gap and make the computer algebra auditable. With that, this will be a reference worth citing.","headline":"Careful and useful classification, held back by an unproved technical step and hidden computer algebra; worth reviewing, not rejecting.","tokens_in":35728,"tokens_out":3355,"would_cite":true,"duration_ms":34515,"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":"For the Standard Model extended by one or two real gauge-singlet scalars, all scalar Lie point symmetries — strict variational, divergence, and non-variational — are classified into a short list of inequivalent algebras, with algorithms tha","keywords":["scalar Lie point symmetries","Standard Model singlet extensions","variational symmetries","divergence symmetries","non-variational symmetries","symmetry classification","Higgs portal"],"falsifier":"Pick a concrete parameter point, e.g. the SM+S free massive singlet (µ_s² ≠ 0, all other singlet couplings zero), write out the five scalar Euler–Lagrange equations, and test whether the system implies any nonzero polynomial identity p(ϕ_1,...,ϕ_5)=0. If such an identity exists, Theorem 1's conclusion that pr X(V) is constant fails, and the s_c⊕u(1)_Y classification of symmetry types would need revision.","tokens_in":34852,"feed_emoji":"⚛️","tokens_out":6724,"duration_ms":59253,"temperature":0.7,"pith_summary":"This paper aims to fully classify the scalar Lie point symmetries — continuous symmetries of the scalar field equations — of the Standard Model extended by one or two real gauge-singlet scalars (SM+S and SM+2S). It claims that, up to affine reparametrization, only four inequivalent symmetry algebras occur for SM+S and eleven for SM+2S, and it separates them into strict variational, divergence, and non-variational types. It also provides a direct parameter-based algorithm that, for any numerical potential, identifies the maximal symmetry algebra of each type without solving the determining equations. A complete classification matters because variational symmetries give conserved currents and constrain model parameters, while non-variational symmetries (such as the free-singlet scaling) are symmetries of the equations of motion only. If correct, the paper gives model-builders a complete catalogue of the continuous scalar symmetries available to singlet-extension models.","feed_headline":"All scalar symmetries of singlet-extended SM classified","feed_subtitle":"Parameter recipes identify which symmetry algebra any SM+S or SM+2S potential realizes, without solving PDEs.","key_machinery":"The load-bearing machinery is Corollary 1, which characterizes the three symmetry types for Lagrangians of the form L = T − V. For a scalar point symmetry X = η^i(φ) ∂_{φ^i} with constant part a^i and linear-potential coefficients α_i, X is a strict variational symmetry if and only if pr X(T)=0 and a^i α_i = 0; a divergence symmetry if pr X(T) is a total divergence and (a^i α_i ≠ 0 or the divergence is nonzero); and non-variational if E(pr X(T)) ≠ 0. The proof rests on the identity pr X(E(L)) = E(pr X(L)) − D*_Q E(L) and on the hypothesis that E(L)=0 has no polynomial consequences, which forces pr X(V) to be a constant. Affine reparametrizations (orthogonal changes of singlet basis and const","core_discovery":"The paper establishes that the scalar Lie point symmetries of the SM+S and SM+2S are completely classified. For the SM+S, the realizable, inequivalent Euler–Lagrange symmetry algebras are exactly a(1)⊕u(1)_Y (a free, massless singlet), s_h⊕u(1)_Y (massless singlet with a linear term), s_c⊕u(1)_Y (free massive singlet), and u(1)_Y (the hypercharge rephasing, always present). For the SM+2S, the list grows to 11 inequivalent Euler–Lagrange algebras, including so(2) rotations corresponding to a complex singlet with a global U(1), the affine algebra a(2) of the free two-singlet kinetic theory, and mixed shift/scaling algebras, with the variational and strict variational subalgebras separately enu","pith_inferences":["The 'no polynomial consequences' hypothesis is the main unverified link; checking it explicitly for representative parameter points via elimination of derivatives (e.g., by computing a Gröbner basis of the differential ideal of the field equations) would confirm the classification's scope.","The parameter-based algorithms could be turned into automated scans that flag which symmetry algebra a given benchmark point in the singlet-extension parameter space realizes, useful for dark matter and baryogenesis studies.","The classification's approach — start from the symmetry algebra of the kinetic terms and reduce by reparametrizations — suggests a route to SM+KS with K>2, though the reduction tree complexity grows rapidly; the paper explicitly notes the boost algebra k⊕u(1)_Y with eigenvalues ±1 is not realizable in SM+2S, a caution for naive symmetry searches."],"forward_implications":["For any numerical SM+S or SM+2S potential, the symmetry algebras g_EL, g_var, and g_svar can be read off from parameter conditions alone, making systematic parameter scans feasible without solving determining equations.","The theorems give a general criterion for when a scalar field symmetry is variational (and hence yields a Noether current) versus merely a symmetry of the field equations, applicable beyond these two models to any Lagrangian with a potential of the assumed form.","The variational symmetries of SM+S reduce to the singlet shift (strict when α=0, divergence otherwise) plus hypercharge; in SM+2S the only variational generators are the two singlet shifts, the rotation s_2∂_{s1} − s_1∂_{s2}, and hypercharge, which tells model-builders which continuous symmetries can be imposed consistently.","Models realizing the a(2)⊕u(1)_Y algebra correspond to a free, massless two-singlet sector; models realizing so(2)⊕u(1)_Y are equivalent to a complex singlet with a global U(1), matching known pseudo-Goldstone dark matter setups."],"fun_headline_variants":["All scalar symmetries of SM+S and SM+2S classified","Parameter recipes classify all scalar symmetries of SM+S and SM+2S","No PDEs needed: new algorithms classify SM+S and SM+2S symmetries","Exact symmetry algebras for singlet-extended SM now known","SM+S and SM+2S: complete classification of scalar symmetries"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The classification of strict variational, divergence, and non-variational symmetries assumes that the SM+S and SM+2S field equations imply no nonzero polynomial relation among the scalar fields, a hypothesis the paper invokes but never verifies.","fun_headline_variants_meta":{"raw":{"variants":["All scalar symmetries of SM+S and SM+2S classified","Parameter recipes classify all scalar symmetries of SM+S and SM+2S","No PDEs needed: new algorithms classify SM+S and SM+2S symmetries","Exact symmetry algebras for singlet-extended SM now known","SM+S and SM+2S: complete classification of scalar symmetries"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000711,"raw_usage":{"total_tokens":3036,"prompt_tokens":743,"completion_tokens":2293,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":487,"completion_tokens_details":{"reasoning_tokens":2196}},"tokens_in":487,"tokens_out":2293,"duration_ms":13817,"temperature":1.0,"reasoning_tokens":2196,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T21:07:12.961021+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Pick a concrete parameter point, e.g. the SM+S free massive singlet (µ_s² ≠ 0, all other singlet couplings zero), write out the five scalar Euler–Lagrange equations, and test whether the system implies any nonzero polynomial identity p(ϕ_1,...,ϕ_5)=0. If such an identity exists, Theorem 1's conclusion that pr X(V) is constant fails, and the s_c⊕u(1)_Y classification of symmetry types would need revision.","supporting_citations":[],"review_version":1}