{"id":"5129c0e9-b3e5-4559-b287-b34a6c541c16","arxiv_id":"2412.04240","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"ESQPTs in constrained bosonic systems can be found from stationary points of a Lagrange function, and the boundary singularities of Holstein-Primakoff mappings disappear when a complete atlas of mappings is used.","lead":"This paper develops a Lagrange-multiplier method to locate and classify excited-state quantum phase transitions in quantum systems with conserved quantities, and shows that standard Holstein-Primakoff mappings hide some of these transitions unless a full atlas of mappings is used. A reader interested in algebraic many-body models gets a more direct route to ESQPTs and a resolution of the long-standing boundary ESQPT puzzle.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Generality claim overreaches: Appendix A's proof of the Lagrange-multiplier ESQPT theorem assumes each constraint is paired with a cyclic action-angle coordinate, i.e., first-class/Abelian constraints; non-commuting conserved quantities are not covered.","rationale":"The reader's weakest_assumption identifies exactly the same gap: the proof of Eqs. (14) and (16) silently assumes the constraints are in involution and admit an action-angle separation. This is the most load-bearing concern because the abstract and Section 2.3 explicitly promise applicability to an arbitrary number of integrals of motion, and the Lagrange method is presented as a general replacement for explicit canonical elimination. If the involution condition is not stated, a reader could apply the method to non-commuting conserved charges (e.g., in superintegrable systems with non-Abelian symmetries) and obtain stationary-point indices that do not correspond to ESQPT singularities. The numerical and semiclassical demonstrations in Sections 3.4–3.5 are convincing for the u(3) model and for the two commuting constraints used there, and the HP-atlas discussion is a useful resolution of the boundary problem. These results are independent support, but they do not test the non-Abelian case. The issue is addressable by adding the involution/first-class condition to the theorem statement or by providing a separate derivation using Dirac brackets. Thus the reader's CONDITIONAL verdict is appropriate; our stress-test does not change it.","tokens_in":19420,"tokens_out":14403,"duration_ms":145958,"concrete_test":"A direct test: take the two-dimensional isotropic harmonic oscillator with two independent conserved charges that do not commute, e.g., Φ1 = H_x − a (the first oscillator's energy) and Φ2 = J − b (the angular momentum), for which {Φ1, Φ2} = 2H_x ≠ 0. Attempt to find a canonical transformation of the form (43) with both φ1 and φ2 cyclic; such a transformation fails for non-commuting charges. Then apply the Lagrange-multiplier equations (11) and the restricted-Hessian index formula (14) to the common level set, and compare the resulting stationary-point indices with the exact quantum level-density singularities in the simultaneous eigenspace of H_x and J (which is generically empty because the two observables do not commute). If the correspondence fails, the theorem's involution assumption is essential.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim (Sections 2.2–2.3, Eqs. (10)–(16)) is that the Lagrange-multiplier method classifies ESQPTs in systems with an arbitrary number of integrals of motion. The proof in Appendix A, however, rests on the construction of a local canonical transformation (43) to coordinates (x, φ, Φ) in which the Hamiltonian is independent of c cyclic coordinates φα. Such a simultaneous action-angle separation exists only when the constraints are in involution (Poisson-commuting) and generate a first-class constraint surface; this is precisely the Liouville–Arnold condition. The paper nowhere states this condition: Section 2.2 only requires the gradients ∇Φα to be nonzero and linearly independent, and calls the constraints arbitrary integrals of motion. For conserved quantities that do not commute, no common action-angle pair exists, the decomposition of the phase-space measure in Eq. (16) is unjustified, and the index identity (14) is not established. The u(3) demonstrations use commuting constraints (Φ_N and Φ_l with {Φ_N, Φ_l}=0), so they provide no evidence for the non-Abelian case. The generality claim is therefore broader than the proof supports.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends the semiclassical theory of excited-state quantum phase transitions (ESQPTs) to systems with constraints induced by conserved quantities. It proposes to find and classify the stationary points of the constrained classical Hamiltonian by solving the Lagrange equations of L = H + Σ_α λ_α Φ_α, claiming that each stationary point's energy and Hessian index coincide with those obtained after explicit reduction to the unconstrained phase space. The method is demonstrated on a u(3) boson model with one constraint (fixed total boson number N) and with two constraints (additional conserved O(2) Casimir). The paper also analyzes Holstein-Primakoff (HP) mappings for fixed-N bosonic systems, showing that a single HP mapping is singular at the boundary of the reduced phase space and that a complete atlas of HP mappings is needed to reveal all ESQPTs. The Lagrange-multiplier results are verified against numerical diagonalization and against the HP atlas.","tokens_in":19616,"tokens_out":16141,"duration_ms":164709,"significance":"If the claimed correspondence is valid, the paper provides a practical method for identifying ESQPTs in constrained systems without constructing explicit canonical transformations, and it resolves a long-standing issue with boundary stationary points in Holstein-Primakoff mappings by showing they are coordinate singularities rather than a distinct type of ESQPT. The numerical verification in the u(3) model with one and two constraints is convincing, and the method contains no fitted parameters: all semiclassical predictions are derived from the Hamiltonian and constraints. The demonstration that an additional constraint can move ESQPT singularities to lower derivatives, change their indices, or remove them entirely is a useful contribution to the classification of ESQPTs. The construction of a complete atlas of HP mappings is a new and clearly explained technical development.","major_comments":[{"comment":"The proof of the central correspondence (13)-(14) relies on the existence of a local canonical transformation to coordinates (x, φ_1,...,φ_c, Φ_1,...,Φ_c) in which the Hamiltonian is independent of the cyclic coordinates φ_α (Eq. (43)). Such a simultaneous action-angle separation exists only when the constraint functions are in involution (Poisson-commuting) and define a first-class constraint surface in the sense of Dirac. The manuscript does not state this condition: Section 2.2 requires only nonzero, linearly independent gradients, and the abstract claims 'an arbitrary number of integrals of motion'. For conserved quantities that do not Poisson-commute, no common action-angle representation is guaranteed, the measure factorization in Eq. (16) is not justified, and the index identity (14) is not established by the given proof. The demonstrations in Sections 3.4 and 3.5 use Poisson-commuting constraints (Φ_N and Φ_l), so they provide no evidence for the non-commuting case. Please either restrict the claim to commuting integrals of motion or supply a proof that avoids the action-angle separation, for example by analyzing the Hessian of H restricted to the constraint surface directly.","section":"Section 2.2 and Appendix A"},{"comment":"The quantum formulation of multiple constraints implicitly requires a common zero eigenspace of the operators ʈΦ_α. For constraints of the form (8), ʈΦ_α = ʈI_α − I_α, this means the conserved operators ʈI_α must possess a joint eigenspace, which generally requires the ʈI_α to commute (or at least to act as scalars on H_c). The paper does not state this commutativity requirement, despite claiming in the abstract that the method applies to an arbitrary number of integrals of motion. This is not merely a technicality: without it, the physical meaning of imposing several non-commuting constraints is unclear. The authors should state explicitly that the method is intended for commuting (Abelian) integrals of motion, or justify that the non-commuting case is physically meaningful and covered by the proof.","section":"Section 2.2, Eq. (7)"}],"minor_comments":[{"comment":"The symbol l^2 is used both for the operator ʈl^2 and for its eigenvalue; the classical counterpart in Eq. (41) uses ℓ^2, which is clearer. Please unify the notation for the eigenvalue throughout Section 3.5.","section":"Eq. (40) and Sec. 3.5"},{"comment":"The block decomposition of the Hessian would be easier to follow if the rows and columns were labeled (x, φ, Φ) explicitly; the current display has three rows of blocks but the second row contains only zero blocks, which is initially confusing.","section":"Eq. (47)"},{"comment":"The expression 'j ≠ k ≠ j′' is ambiguous; it should read 'k ≠ j and k ≠ j′'.","section":"Appendix B, Eq. (57)"},{"comment":"The statement that the ESQPT at E = −0.2 'can be considered as a result of an interplay of all the subspaces' is vague; a brief explanation of how the finite-dimensional subspaces H_l conspire to produce the semiclassical singularity would strengthen the argument.","section":"Section 3.5, final paragraph"}],"recommendation":"major_revision","confidential_remarks":"The paper is well within the scope of J. Phys. A and the numerical verification is convincing. The main issue is the unstated involutivity/commutativity assumption; if the authors are willing to either restrict the generality claim or prove the correspondence under weaker assumptions, the paper should be publishable. The heavy citation of the authors' own prior work is appropriate given the topic."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this paper is worth engaging with. The Lagrange-multiplier method for finding and classifying ESQPT stationary points is a genuine, practical extension of the standard semiclassical machinery, and the Holstein-Primakoff atlas argument is the strongest part: it convincingly shows that boundary ESQPTs in the u(3) vibron model are artifacts of a single singular mapping, not a separate phenomenon. The numerical checks against Gaussian-smoothed quantum spectra are convincing, with stationary-point energies and indices matching the observed nonanalytic structure. No parameters are fitted to the target spectra. The citation pattern is normal; the heavy use of the authors' own prior classification is background, not circularity.\n\nWhat is genuinely new: the constrained stationary-point and index theorem, the demonstration that a complete HP atlas resolves the boundary singularity, and the multiply-constrained u(3) example showing that additional integrals of motion can change the ESQPT count and character. The paper also gives a useful practical argument that Lagrange multipliers avoid the square-root clutter and boundary issues of HP mappings.\n\nThe main soft spot is exactly what the stress-test identifies. Appendix A proves the theorem only under the assumption that the constraints can be paired with cyclic coordinates via a local action-angle canonical transformation, which requires the constraints to Poisson-commute and be independent. The paper never states this involution condition; Section 2.2 only requires nonzero, linearly independent gradients and calls the constraints arbitrary integrals of motion. That overstates the generality. The proof of the level-density equality (16) also relies on the same measure factorization, so the unqualified claim of \"arbitrary number of constraints\" is not established. The demonstrated examples all use commuting constraints, so the gap does not infect the numerics, but it is a real gap in the general theorem. A second, smaller gap: the claim that a complete HP atlas always suffices is argued and demonstrated for f=2, but not proven for general f or higher-rank algebras. That is minor relative to the main result.\n\nI think the central argument holds for the cases it actually treats. The paper is clearly written and the authors know the ESQPT literature well. The right fix is to make the involution/separability assumption explicit and qualify the generality claim, not to redo the analysis.\n\nWho is this for? Anyone working on ESQPTs in bosonic algebraic models—vibron, IBM, spinor condensates—will get direct value. It deserves a serious referee. I would send it to review with the request that the authors state and prove their theorem under the correct hypotheses, or explicitly restrict the claims to commuting integrals of motion. I would cite it in my own work on ESQPTs.","headline":"A useful, mostly correct extension of ESQPT semiclassics whose advertised generality outruns its proof.","tokens_in":20162,"tokens_out":2276,"would_cite":true,"duration_ms":30064,"reading_group":"yes","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 claims that all excited-state quantum phase transitions in constrained Hamiltonian systems can be found and classified directly from the Hamiltonian plus Lagrange multipliers, without constructing the canonical transformation…","keywords":["excited-state quantum phase transitions","Lagrange multipliers","constrained systems","Holstein-Primakoff mapping","level density","stationary points","algebraic boson models","u(3) model"],"falsifier":"Choose a Hamiltonian with two independent, regular constraints whose Poisson brackets do not vanish, compute the quantum level density at large size, and compare the locations and derivative-singularity types of its nonanalytic features with predictions from $\\nabla L=0$ and $D^2L|_{\\Sigma}$; a mismatch would falsify the claimed generality. Within the u(3) model, a numerical check near $\\xi\\approx0.42$ should show the derivative of the smoothed level density developing singular kinks exactly at the new stationary energies for $\\xi>\\xi_e$.","tokens_in":19211,"feed_emoji":"⚛️","tokens_out":7376,"duration_ms":71481,"temperature":0.7,"pith_summary":"The paper extends the standard semiclassical theory of excited-state quantum phase transitions (ESQPTs) to systems whose dynamics is restricted by conserved quantities. Its central claim is that all stationary points of the reduced classical Hamiltonian—the objects that produce ESQPT singularities in the level density—can be found and classified directly from the original Hamiltonian together with Lagrange multipliers for the constraints. This removes the need to construct the canonical transformation that eliminates the constrained degrees of freedom. The paper proves the correspondence on general grounds and demonstrates it on a u(3) boson model with one and two constraints. It also shows that boundary stationary points that appear in a single Holstein-Primakoff mapping are an artifact of that singular coordinate choice: a complete atlas of mappings reveals them as ordinary ESQPTs.","feed_headline":"One formula exposes every hidden phase transition in constrained systems","feed_subtitle":"Excited-state singularities in the level density are classified straight from the Hamiltonian, with no coordinate surgery.","key_machinery":"The Lagrange function $L(X,\\lambda)=H(X)+\\sum_\\alpha\\lambda_\\alpha\\Phi_\\alpha(X)$ is the primary object: solving its full gradient gives stationary points $\\{X_{\\rm st},\\lambda_{\\rm st}\\}$. The second object is the restricted Hessian $D^2L|_{\\Sigma}$, obtained by evaluating the Hessian of $L$ only in directions tangent to the constraint surface $\\Sigma$; its negative-eigenvalue count supplies the index $r_\\Sigma$. Appendix A shows that a local canonical transformation $X\\mapsto(x,\\varphi,\\Phi)$ separates the constraints into conjugate pairs $(\\Phi_\\alpha,\\varphi_\\alpha)$ with cyclic angles, which turns the restricted Hessian into the Hessian of the reduced Hamiltonian plus harmless zero eigenvalues. The Holstein-Primakoff map $M^{(j)}$ is the third object: a singular projection from the constraint sphere $\\Sigma$ onto the compact ball $\\sigma^{(j)}$, and the atlas $\\{M^{(0)},\\dots,M^{(f)}\\}$ is what guarantees that every stationary point appears in the interior of at least one chart.","core_discovery":"For a constrained system with $f+c$ degrees of freedom and $c$ independent constraints $\\Phi_\\alpha$, define $L(X,\\lambda)=H(X)+\\sum_\\alpha\\lambda_\\alpha\\Phi_\\alpha(X)$. The paper's claim is that the solutions of $\\nabla_{(X,\\lambda)}L=0$ stand in one-to-one correspondence with the stationary points $x_{\\rm st}$ of the reduced Hamiltonian $H^{(\\sigma)}$ on the physical phase space, with equal energies $E_{\\rm st}=H(X_{\\rm st})=L(X_{\\rm st},\\lambda_{\\rm st})$ and equal singularity indices $r(x_{\\rm st})=r_\\Sigma(X_{\\rm st})$, the latter read off from the Hessian of $L$ restricted to directions tangent to the constraint surface. This is proven in Appendix A using local action-angle coordinates in which each constraint is a momentum conjugate to a cyclic angle. Applied to the u(3) boson model, the method reproduces the full set of ESQPT energies and indices, including the transition at $E=1-\\xi$ that a single Holstein-Primakoff chart misses. The companion claim about the Holstein-Primakoff mapping is that its classical singular boundary is a coordinate artifact: boundary stationary points move to the interior of some chart of a complete atlas of HP mappings, so they are not a separate class of ESQPT.","pith_inferences":["The same Lagrange construction should apply to systems whose constraints are not integrals of motion but are imposed by the physical setup, such as gauge constraints, provided the involutive and action-angle conditions hold.","The appearance of complex stationary points near the real axis suggests a general precursor effect: before a pair of real stationary points emerges, the smoothed level density should show a smooth but non-monotonic wiggle, which could serve as an early spectroscopic warning of an approaching ESQPT.","For boson systems with more than two degrees of freedom, the atlas of HP maps grows with the number of boson types, so the Lagrange method likely becomes the only practical route; this could be checked by extending the u(3) analysis to u(4) vibron models.","The u(3) model with two constraints is integrable; applying the Lagrange method to other integrable limits should predict which invariant subspaces carry ESQPTs and which do not."],"forward_implications":["ESQPT energies and singularity types in any constrained system follow from solving $\\nabla L=0$ and diagonalizing $D^2L|_{\\Sigma}$; the explicit canonical transformation that removes the constraints is never needed.","Boundary stationary points seen in a single Holstein-Primakoff map are ordinary ESQPTs, not a separate class; a complete atlas of HP maps finds every one in the interior of some chart.","Adding extra conserved constraints lowers the effective number of degrees of freedom, moving ESQPT singularities to lower derivatives of the level density and changing the indices.","Additional constraints can also delete ESQPTs: some stationary points live only in subspaces with particular values of the conserved quantities, and at least one singularity in the doubly constrained u(3) model arises from the interplay of all invariant subspaces.","For polynomial algebraic Hamiltonians the Lagrange equations are polynomial, so stationary points can be found analytically and their count bounds the number of ESQPTs."],"supporting_citations":[{"why":"Supplies the classification of ESQPTs by stationary-point index for arbitrary numbers of degrees of freedom, the framework Eq. (6) extends to constrained systems.","marker":"[27]"},{"why":"Provides the Holstein-Primakoff treatment of interacting boson systems whose boundary singularities the paper reinterprets via a complete atlas.","marker":"[39]"},{"why":"Review of ESQPT semiclassical theory that frames the level-density singularity analysis used throughout.","marker":"[22]"},{"why":"Gives the Weyl formula and phase-space integration used to equate the constrained and reduced level densities in Eq. (16).","marker":"[46]"},{"why":"Supplies the Lagrange multiplier method for constrained optimization that the paper adopts to find stationary points.","marker":"[50]"},{"why":"Gives the coherent-state classical limit of the u(3) vibron model and the QPT at $\\xi=1/5$ used as a reference case.","marker":"[43]"},{"why":"Reports stationary points hidden at the phase-space boundary in atom-field systems and conjectures a different nature, which the paper addresses.","marker":"[14]"},{"why":"Identifies the ESQPT in the $\\ell=0$ subspace of the u(3) model, the doubly constrained result reproduced by the Lagrange method.","marker":"[62]"}],"fun_headline_variants":["One formula reveals all ESQPTs in constrained systems","Lagrange multipliers expose every ESQPT in constrained systems","Complete atlas of HP mappings uncovers all ESQPTs","New theory classifies constrained ESQPTs without coordinate surgery","All ESQPTs in constrained systems from a single formula"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument relies on every constraint being pairable with a cyclic coordinate through a local action-angle canonical transformation; if the conserved quantities fail to commute or no such separation exists, the equality of energies and indices is not proven.","fun_headline_variants_meta":{"raw":{"variants":["One formula reveals all ESQPTs in constrained systems","Lagrange multipliers expose every ESQPT in constrained systems","Complete atlas of HP mappings uncovers all ESQPTs","New theory classifies constrained ESQPTs without coordinate surgery","All ESQPTs in constrained systems from a single formula"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001136,"raw_usage":{"total_tokens":4733,"prompt_tokens":974,"completion_tokens":3759,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":590,"completion_tokens_details":{"reasoning_tokens":3675}},"tokens_in":590,"tokens_out":3759,"duration_ms":25582,"temperature":1.0,"reasoning_tokens":3675,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T21:36:45.252375+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Choose a Hamiltonian with two independent, regular constraints whose Poisson brackets do not vanish, compute the quantum level density at large size, and compare the locations and derivative-singularity types of its nonanalytic features with predictions from $\\nabla L=0$ and $D^2L|_{\\Sigma}$; a mismatch would falsify the claimed generality. Within the u(3) model, a numerical check near $\\xi\\approx0.42$ should show the derivative of the smoothed level density developing singular kinks exactly at the new stationary energies for $\\xi>\\xi_e$.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the classification of ESQPTs by stationary-point index for arbitrary numbers of degrees of freedom, the framework Eq. (6) extends to constrained systems."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Holstein-Primakoff treatment of interacting boson systems whose boundary singularities the paper reinterprets via a complete atlas."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Review of ESQPT semiclassical theory that frames the level-density singularity analysis used throughout."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the Weyl formula and phase-space integration used to equate the constrained and reduced level densities in Eq. (16)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Lagrange multiplier method for constrained optimization that the paper adopts to find stationary points."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the coherent-state classical limit of the u(3) vibron model and the QPT at $\\xi=1/5$ used as a reference case."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports stationary points hidden at the phase-space boundary in atom-field systems and conjectures a different nature, which the paper addresses."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Identifies the ESQPT in the $\\ell=0$ subspace of the u(3) model, the doubly constrained result reproduced by the Lagrange method."}],"review_version":1}