{"id":"e329aa0c-7cca-4f1e-8320-180428d8e37f","arxiv_id":"1909.01997","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"An exact ground-state purity formula for three coupled harmonic oscillators is derived via SU(3) rotation matrices, recovering known two-oscillator results in decoupling limits.","lead":"This paper diagonalizes a Hamiltonian of three coupled harmonic oscillators and derives the ground-state entanglement, expressed as a purity function. The result generalizes a previously published two-oscillator formula to the three-mode case using SU(3) rotations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The exact solution and purity formula are valid only when the potential matrix R in Eq (6) is positive definite; the paper neither states this nor keeps the strong-coupling limit inside the stable domain.","rationale":"After rederiving the diagonalization, I find the SU(3)-based rotation is a valid SO(3) parameterization and the expressions (12)-(17) are consistent; the purity formula (40) also reduces correctly in the three decoupling limits, including the known two-oscillator result (60). The one load-bearing gap is the unstated stability condition: the ground state and all subsequent formulas exist only for R>0. The paper's strong-coupling limit is the concrete casualty because it lets the auxiliary differences (43) go to +-infinity without specifying a physical path that keeps R positive definite. This is a domain restriction rather than an algebraic error, so the reader's CONDITIONAL verdict is appropriate. I do not see a separate flaw that would justify rejection.","tokens_in":10919,"tokens_out":19681,"duration_ms":184963,"concrete_test":"Choose the equal-mass, equal-frequency case with m_i=1, omega_i=1, and J_12=J_13=J_23=2, so that Eq (6) is R = I + 2(1-delta_ij). Show that the eigenvalues of R are 5, -1, -1 and attempt to insert the resulting normal-mode data into Eq (19) and Eq (40): the quantity e^{rho-kappa} is not a positive real number, so the purity formula cannot be evaluated as written. If the authors instead restrict to R>0, verify Eq (40) numerically for a positive-definite example, e.g. m_i=1, omega=(1,2,3), J_ij chosen to satisfy |J_ij|<omega_i omega_j, by comparing with the purity obtained from direct diagonalization of the covariance matrix.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is an exact ground-state solution and a purity function for 'arbitrary masses and couplings.' For that claim to hold, the potential matrix R in Eq (6) must be positive definite: only then do the normal-mode frequencies Sigma_i exist as real numbers, the parameters in Eq (19) make sense, and the ground state in Eq (23) be normalizable. The manuscript never states or verifies this condition, and the strong-coupling limit in Section 4 is where the issue matters. Taking the couplings J_ij to infinity while holding the diagonal frequencies fixed drives R indefinite; for example, with m_i=1, omega_i=1, and J_12=J_13=J_23=2, R has eigenvalues 5, -1, -1. In that regime one of the Sigma_i^2 in Eq (11) is negative, the logarithms in Eq (19) are not real, and Eq (40) is not a physical purity. The statement in Section 6 that the results are 'general and derived without making use of any assumption or approximation' is therefore too strong: the derivation is conditional on R>0, and the strong-coupling limit must be taken along a path that preserves positive definiteness.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents an exact treatment of three coupled harmonic oscillators with arbitrary masses, frequencies, and quadratic position couplings. After a mass rescaling, the potential matrix R is diagonalized by an orthogonal transformation M parameterized by three angles and expressed through Gell-Mann matrices; the resulting normal modes define three decoupled oscillators. The ground-state wavefunction is written in the original coordinates, and after tracing out two oscillators the reduced state is a single-mode Gaussian. The paper's central result is the closed-form purity Eq. (40) as a function of the three normal-mode ratios and three angles. The authors discuss weak- and strong-coupling limits and show, in three decoupling limits, that the result reduces to the two-oscillator formula of Ref. [22].","tokens_in":11170,"tokens_out":29542,"duration_ms":244931,"significance":"If correct, Eq. (40) is a useful closed-form expression for bipartite entanglement of a three-mode harmonic system with arbitrary quadratic interactions, expressed through normal-mode data. The diagonalization is self-contained, the Gaussian reduction is standard, and the recovery of the known two-oscillator result in Section 5 provides an independent check. The paper includes no fitting parameters, and the formula is directly falsifiable by independent Gaussian-state calculations. The main weaknesses are the missing stability condition on R and the incorrect limiting formulas, which currently prevent the consistency check from being conclusive.","major_comments":[{"comment":"The purity formula and the strong-coupling limit require the potential matrix R in Eq. (6) to be positive definite, since otherwise the normal-mode frequencies Sigma_i in Eq. (11) are not real, the parameters in Eq. (19) are not real, and the ground state in Eq. (23) is not normalizable. The manuscript never states or verifies this condition. For example, with m_i=1, omega_i=1, and J_12=J_13=J_23=2, R has eigenvalues 5, -1, -1, so the Hamiltonian has no ground state and Eq. (40) is not a physical purity. The limit in Eq. (43) must be taken along a path that keeps R positive definite; the sentence in Section 6 claiming the results are derived without any assumption or approximation is therefore too strong.","section":"Section 4, Eqs. (40)-(44)"},{"comment":"The definitions in Eq. (55) are inconsistent with the normalization of J_12 given in Eq. (4). Since J_12 is the off-diagonal element of R, the determinant of the two-oscillator block is omega_1^2 omega_2^2 - J_12^2, not omega_1^2 omega_2^2 - J_12^2/4, and the eigenvalues are [omega_1^2+omega_2^2 +/- sqrt((omega_1^2-omega_2^2)^2 + 4J_12^2)]/2. The correct relations are k_12 = sqrt(omega_1^2 omega_2^2 - J_12^2) and e^{+-2eta_12} = [omega_1^2+omega_2^2 +/- sqrt((omega_1^2-omega_2^2)^2 + 4J_12^2)]/(2k_12). As printed, substituting Eq. (55) into Eqs. (54) and (60) does not reproduce the known two-oscillator purity, so the consistency check with Ref. [22] fails unless these formulas are corrected.","section":"Section 5, Eq. (55)"},{"comment":"The central purity formula Eq. (40) is obtained by a 'straightforward calculation' that is not shown. Because Eq. (40) is the main quantitative result of the paper, the authors should include the intermediate algebra, or an appendix, demonstrating how the Gaussian reduction from Eqs. (37)-(38) leads to the factorized form (40). The limiting checks in Section 5 are helpful but do not replace the derivation.","section":"Section 4, Eqs. (37)-(40)"},{"comment":"The paper claims to solve the problem for arbitrary masses and couplings, but Eqs. (12)-(17) give the physical parameters (omega_i, J_ij) in terms of the normal-mode data (Sigma_i, angles), not the inverse. To use Eq. (40) for a specified Hamiltonian, one still needs to diagonalize R. The authors should state this explicitly or provide the forward map from the original parameters to the angles and normal-mode frequencies.","section":"Sections 2-3, Eqs. (12)-(17)"}],"minor_comments":[{"comment":"The potential term in Eq. (5) should be (m/2) sum X_i R_ij X_j, not (1/(2m)) sum X_i R_ij X_j; the factor as printed makes the equation dimensionally inconsistent with Eq. (2) and with the later expression in Eq. (18).","section":"Eq. (5)"},{"comment":"In Eq. (37), the factor printed as B2 should read B, so that the denominator is B - Gamma_23^2/C; otherwise the Gaussian integration leading to Eq. (38) is not reproduced.","section":"Eq. (37)"},{"comment":"The phrase 'tow rather than three free parameters' contains a typo and should read 'two rather than three free parameters'.","section":"Section 3"},{"comment":"Ref. [24] cites the Wikipedia page for Gell-Mann matrices; this should be replaced by a standard textbook reference.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper's use of SU(3) is essentially a parametrization of SO(3) by Euler angles; the representation-theory content is light. The main value is the explicit closed formula, which is useful as a reference, but the manuscript would benefit from a more careful derivation of Eq. (40) and from correcting the limiting formulas. Novelty is moderate; I would not reject on that basis if the technical issues are fixed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this paper delivers what it says in the title—an exact diagonalization of three coupled harmonic oscillators and an explicit ground-state purity function. The purity formula (40) is new as far as I know and is a genuine extension of the authors' two-oscillator result; it is parameter-free, derived from the Hamiltonian data, and it reduces correctly to [22] in three decoupling limits. The SU(3) rotation is a rephrasing of a standard orthogonal diagonalization, so the novelty is not the method; the value is the closed-form expression for the reduced single-mode purity in terms of the six parameters. That is useful for three-mode Gaussian states, and the check against earlier work gives me real confidence the algebra is right.\n\nThe soft spots are real but not fatal. The derivation assumes the potential matrix R in Eq (6) is positive definite: only then are the normal-mode frequencies real and the ground state normalizable. The paper never states this. The stress-test example is concrete: with unit masses, unit frequencies and all three couplings equal to 2, R has a negative eigenvalue, so the purity formula is not a physical purity. The strong-coupling limit in Section 4 takes (ς−ρ, κ−ς, ρ−κ) to infinity in a way that, if interpreted as sending the bare couplings to infinity with fixed diagonal frequencies, leaves the stable domain. That limit needs to be reformulated along a path that keeps R>0, or at least with the condition stated.\n\nMinor issues: the step from Eqs (37)–(38) to Eq (40) is not shown, and it is not a one-liner; there is a typo in Eq (37), where the denominator should be B, not B^2; and the conclusion's claim that everything was derived 'without making use of any assumption or approximation' is too strong.\n\nBottom line: the central result is honest progress on a standard model. It should be sent to peer review; a good referee will ask for the missing derivation and a clear statement of the stability condition. I would not desk-reject it.","headline":"Exact three-oscillator purity formula is real progress but rests on an unstated stability condition; the strong-coupling limit as written leaves the valid domain.","tokens_in":11673,"tokens_out":2510,"would_cite":true,"duration_ms":25087,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["03.67.Bg","03.65.-w","02.20.Sv"],"model":"deepseek-v4-flash","headline":"This paper derives the exact ground-state purity function for three coupled harmonic oscillators with arbitrary quadratic couplings, showing that the entanglement is fully determined by three normal-mode ratios and three rotation angles.","keywords":["three coupled harmonic oscillators","SU(3) representation theory","reduced density matrix","purity function","Gaussian states","entanglement","normal-mode diagonalization","quadratic couplings"],"falsifier":"Take a concrete positive-definite potential matrix, for instance unit bare frequencies with couplings $J_{12}=0.4$, $J_{13}=0.2$, $J_{23}=0.3$, and numerically diagonalize $R$ to obtain the ground-state covariance matrix; trace out oscillators 2 and 3 and compute the purity, comparing it with Eq. (40) using the angles and ratios read from that same diagonalization. A boundary check exposing the hidden assumption is the symmetric case of unit frequencies with all couplings equal to $J>1$, where $R$ has a negative eigenvalue and Eq. (40) is formally evaluated outside its domain of validity.","tokens_in":10739,"feed_emoji":"🔗","tokens_out":20398,"duration_ms":186113,"temperature":0.7,"pith_summary":"The paper aims to solve exactly the quantum mechanics of three harmonic oscillators of arbitrary masses and frequencies coupled by arbitrary quadratic coordinate couplings. By diagonalizing the three-by-three potential matrix with a three-angle rotation built from SU(3) generators, the Hamiltonian becomes three independent oscillators, and the spectrum and eigenfunctions follow in closed form. From the ground state the authors compute the reduced density matrix after tracing out two oscillators and obtain an explicit formula for the purity, a number that measures how entangled the remaining oscillator is with the other two. The formula gives purity 1 for zero couplings (a product state) and drives the purity to 0 for strong couplings (maximal entanglement), and it reproduces the known two-oscillator purity in three decoupling limits.","feed_headline":"Three coupled oscillators: exact ground-state entanglement formula","feed_subtitle":"The purity depends on six parameters and recovers the two-oscillator result in three decoupling limits.","key_machinery":"The load-bearing object is the real orthogonal three-angle rotation $M$ of Eq. (9), built from the imaginary Gell-Mann matrices, which diagonalizes the real symmetric potential matrix $R$ of Eq. (6) into $\\mathrm{diag}(\\Sigma_1^2,\\Sigma_2^2,\\Sigma_3^2)$ and turns the interacting Hamiltonian into three decoupled oscillators. The normal-mode frequencies are reparametrized by their geometric mean and three exponential ratios, and the purity formula is the result of a Gaussian integration over the two traced oscillators expressed entirely in those ratios and the three angles. The same machinery supplies the full energy spectrum and the exact eigenfunctions, so the entanglement calculation rests on the same diagonalization that solves the spectrum.","core_discovery":"The central claim is that the ground-state purity of the three-oscillator system is exactly $$P = \\prod_{i=1}^3 \\left( M_{i1}^2 $e^{{\\rho-\\varsigma}}$ + M_{i2}^2 $e^{{\\varsigma-\\kappa}}$ + M_{i3}^2 $e^{{\\kappa-\\rho}}$ \\right)^{-1/2},$$ where $M$ is the three-angle rotation of Eq. (9) and $(\\rho,\\varsigma,\\kappa)$ are the exponential normal-mode parameters of Eq. (19); expanding the rows of $M$ gives the explicit six-parameter expression in Eq. (40). In the weak-coupling limit the rotation reduces to the identity and the exponentials become the bare frequency ratios, so the formula returns $P=1$; in the strong-coupling limit the exponentials become extreme and the formula gives $P\\to 0$. The same expression reduces to the earlier two-oscillator purity $P_{0,0}(\\eta_{12},\\varphi)$ when the third oscillator is decoupled, and to its two partner expressions in the other decoupling limits.","pith_inferences":["The product-of-three-square-roots structure of Eq. (40) suggests that the total purity factorizes into three directional factors associated with the rows of the rotation matrix; the paper does not state this, but it implies the tripartite entanglement can be tuned by controlling the parameters of a single spatial direction in the diagonalizing frame.","The same SU(3) diagonalization should extend to the time-dependent Hamiltonian by allowing the three angles and the exponential ratios to depend on time, giving an exact time-dependent purity; the authors list this as future work without carrying it out.","Because the formula is only valid where the potential matrix is positive definite, plotting the purity over coupling space would reveal a stability boundary where the ground state ceases to exist; the paper does not map this boundary."],"forward_implications":["For zero couplings the purity returns 1, so the ground state is fully separable; the formula therefore quantifies how much the quadratic couplings create entanglement.","In the strong-coupling limit the purity goes to 0, which the paper interprets as maximal entanglement of the ground state.","The exact eigenfunctions expressed back in the original coordinates give a complete basis for computing any observable of the three-body system, not only the ground-state purity.","When one oscillator is decoupled, the purity formula reduces to the two-oscillator purity in each of the three possible decoupling limits, so the three-body result contains the two-body result.","Taking the special case where one normal-mode frequency vanishes recovers the known propagator result for a three-dimensional three-body system, indicating the diagonalization also covers a free-particle limit."],"supporting_citations":[{"why":"Provides the two-oscillator purity result that the three-body formula must reproduce in three decoupling limits, serving as the consistency benchmark.","marker":"[22]"},{"why":"Demonstrates the use of Lie-algebra representation theory for three-mode coupled systems, the method the paper generalizes via SU(3).","marker":"[13]"},{"why":"The quantum-propagator solution for three-dimensional three-body systems that the paper recovers in the limit of one zero normal-mode frequency, validating the diagonalization's reach.","marker":"[23]"}],"fun_headline_variants":["Exact entanglement formula for three coupled oscillators","Ground-state purity: exact for three oscillators","Three coupled modes: exact purity solved","Six-parameter purity exact for triple oscillator","Three oscillator entanglement: exact ground-state"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation assumes the potential matrix $R$ in Eq. (6) is positive definite so that the normal-mode frequencies are real; the paper neither states nor verifies this condition, and its strong-coupling limit approaches the boundary of that domain.","fun_headline_variants_meta":{"raw":{"variants":["Exact entanglement formula for three coupled oscillators","Ground-state purity: exact for three oscillators","Three coupled modes: exact purity solved","Six-parameter purity exact for triple oscillator","Three oscillator entanglement: exact ground-state"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000261,"raw_usage":{"total_tokens":1547,"prompt_tokens":855,"completion_tokens":692,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":471,"completion_tokens_details":{"reasoning_tokens":626}},"tokens_in":471,"tokens_out":692,"duration_ms":6146,"temperature":1.0,"reasoning_tokens":626,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:05:18.456067+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a concrete positive-definite potential matrix, for instance unit bare frequencies with couplings $J_{12}=0.4$, $J_{13}=0.2$, $J_{23}=0.3$, and numerically diagonalize $R$ to obtain the ground-state covariance matrix; trace out oscillators 2 and 3 and compute the purity, comparing it with Eq. (40) using the angles and ratios read from that same diagonalization. A boundary check exposing the hidden assumption is the symmetric case of unit frequencies with all couplings equal to $J>1$, where $R$ has a negative eigenvalue and Eq. (40) is formally evaluated outside its domain of validity.","supporting_citations":[{"cited_title":"Jellal, F","cited_arxiv_id":null,"evidence_quote":"Provides the two-oscillator purity result that the three-body formula must reproduce in three decoupling limits, serving as the consistency benchmark."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Demonstrates the use of Lie-algebra representation theory for three-mode coupled systems, the method the paper generalizes via SU(3)."},{"cited_title":"de Souza Dutra, Ann","cited_arxiv_id":null,"evidence_quote":"The quantum-propagator solution for three-dimensional three-body systems that the paper recovers in the limit of one zero normal-mode frequency, validating the diagonalization's reach."}],"review_version":1}