Pith. sign in

REVIEW 4 major objections 4 minor 24 references

Entanglement in Three Coupled Harmonic Oscillators

T0 review · 4 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read 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.

desk verdict 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. read the letter →

arxiv 1909.01997 v2 pith:77HEWIUW submitted 2019-09-04 quant-ph hep-th

Abdeldjalil Merdaci , Ahmed Jellal This is my paper · ORCID
classification quant-phhep-th PACS 03.67.Bg03.65.-w02.20.Sv
keywords threecoupledharmonicoscillatorsSU(3)representationtheoryreduceddensitymatrixpurityfunctionGaussianstatesentanglementnormal-modediagonalizationquadraticcouplings
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

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].

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 (4)
  1. [Section 4, Eqs. (40)-(44)] 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.
  2. [Section 5, Eq. (55)] 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.
  3. [Section 4, Eqs. (37)-(40)] 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.
  4. [Sections 2-3, Eqs. (12)-(17)] 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.
minor comments (4)
  1. [Eq. (5)] 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).
  2. [Eq. (37)] 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.
  3. [Section 3] The phrase 'tow rather than three free parameters' contains a typo and should read 'two rather than three free parameters'.
  4. [References] Ref. [24] cites the Wikipedia page for Gell-Mann matrices; this should be replaced by a standard textbook reference.

Circularity Check

0 steps flagged · score 2.0 of 10

No circular derivation: normal modes and purity are computed from the physical parameters; the prior-work comparison is a consistency check, not an input.

full rationale

The derivation chain is self-contained. Eq. (5)-(18) diagonalize the Hamiltonian (1) by an explicit SU(3) rotation M (Eq. 9-10), so the normal-mode frequencies Sigma_i and angles are functions of the masses, frequencies and couplings, not fitted parameters. Eq. (40) is obtained by Gaussian integration over the traced-out variables from the normalized ground state (23), and the parameters entering it are the same diagonalization variables; none of the target quantities (purity limits P=1 and P->0) are inserted as an ansatz. The references to the authors' earlier two-oscillator paper [22] appear only in Section 5 and Appendix B as a limiting consistency check: after taking J13,J23 -> 0, the general formulas reduce to Eqs. (51)-(60), and the paper states the result 'coincides exactly with that obtained in our previous work [22]'. That comparison is not used to establish the general spectrum or purity, so it is a minor self-citation, not load-bearing. The only substantive caveat is a domain condition: the matrix R in Eq. (6) must be positive definite for the Sigma_i^2 to be positive and the ground state to be normalizable, and the statement in Section 6 that the result holds 'without making use of any assumption or approximation' should be read as conditional on that stability condition. This is a mathematical-validity concern, not circularity.

Assumptions & free parameters 0 free parameters · 3 assumptions · 0 invented entities

No parameters are fitted to data. The angles (θ,φ,ϕ) and the normal-mode frequencies (Σ1,Σ2,Σ3) are determined by the Hamiltonian's masses, frequencies, and couplings. No new entities are introduced.

assumptions (3)
  • domain assumption The potential matrix R is positive definite, so all normal-mode frequencies Σ_i are real.
    Eq (11) writes R = M diag(Σ^2) M^{-1} and Eq (19) takes cube roots of Σ_i; if R has a negative eigenvalue, the ground state is not normalizable. The paper never states this stability condition.
  • standard math The rotation matrix M in Eq (10) is orthogonal and provides a valid canonical transformation for all physical parameters.
    This is the Euler-angle decomposition of SO(3); the paper uses it without proof, which is acceptable.
  • standard math The Gaussian integration over x2 and x3 that yields the reduced density matrix (36) is elementary and correct.
    The result is a standard marginal of a multivariate Gaussian; the paper states it without derivation, and we spot-checked the algebra.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Entanglement in Three Coupled Harmonic Oscillators." pith.science (2026). https://pith.science/paper/77HEWIUW

@misc{pith2026190901997,
  author       = {Pith},
  title        = {Pith review of: Entanglement in Three Coupled Harmonic Oscillators},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/77HEWIUW}},
  note         = {Machine review of arXiv:1909.01997}
}
abstract

We develop an approach in solving exactly the problem of three-body oscillators including general quadratic interactions in the coordinates for arbitrary masses and couplings. We introduce a unitary transformation of three independent angles to end up with a diagonalized Hamiltonian. Using the representation theory of the group $SU(3)$, we explicitly determine the solutions of the energy spectrum. Considering the ground state together with reduced density matrix, we derive the corresponding purity function that is giving rise to minimal and maximal entanglement under suitable conditions. The cases of realizing one variable among three is discussed and know results in literature are recovered.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

24 extracted references · 24 canonical work pages

  1. [22]

    Jellal, F

    A. Jellal, F. Madouri and A. Merdaci, J. Stat. Mech. (201 1) P09015

    work page

  2. [1]

    Schr¨ odinger, Naturwissenschaften 23 (1935) 807

    E. Schr¨ odinger, Naturwissenschaften 23 (1935) 807

    work page 1935

  3. [2]

    Einstein, B

    A. Einstein, B. Podolsky and N. Rosen, Phys. Rev. 47 (1935) 777

    work page 1935

  4. [3]

    Bell, ”Speakable and Unspeakable in Quantum Mechanics” (Cambridge University Press, 1987)

    J.S. Bell, ”Speakable and Unspeakable in Quantum Mechanics” (Cambridge University Press, 1987)

    work page 1987

  5. [4]

    Benett and Peter W

    Charles H. Benett and Peter W. Shor, IEEE Transactions on Information Theory 44 (1998) 2724

    work page 1998

  6. [5]

    C. H. Bennett, G. Brassard, C. Cr´ epeau, R. Jozsa, A. Pere s and W. K. Wootters, Phys. Rev. Lett. 70 (1993) 1895

    work page 1993

  7. [6]

    C. H. Bennett and S. J. Wiesner, Phys. Rev. Lett. 69 (1992) 2881

    work page 1992

  8. [7]

    A. K. Ekert, Phys. Rev. Lett. 67 (1991) 661

    work page 1991

Show all 24 references
  1. [8]

    Murao, D

    M. Murao, D. Jonathan, M. B. Plenio and V. Vedral, Phys. Re v. A 59 (1999) 156

    1999

  2. [9]

    C. A. Fuchs, Phys. Rev. Lett. 79 (1997) 1162

    1997

  3. [10]

    Gottesman and I

    D. Gottesman and I. Chuang, Nature 402 (1999) 390. 10

    1999

  4. [11]

    Horne, Harald Weinfurter a nd Marek Åżukowski, Phys

    Anton Zeilinger, Michael A. Horne, Harald Weinfurter a nd Marek Åżukowski, Phys. Rev. Lett. 78 (997) 3031

  5. [12]

    Park, Quantum Inf Process 18 (2019) 282

    D. Park, Quantum Inf Process 18 (2019) 282

    2019

  6. [13]

    M. M. Sebawe Abdalla and M. A. Bashir, Quantum Semiclass . Opt. 10 (1998) 415

    1998

  7. [14]

    Giedke, B

    G. Giedke, B. Kraus, M. Lewenstein and J. I. Cirac, Phys. Rev. A 64 (2001) 052303

    2001

  8. [15]

    Bouwmeester, J.-W

    D. Bouwmeester, J.-W. Pan, M. Daniell, H. Weinfurter an d A. Zeilinger, Phys. Rev. Lett. 82 (1999) 1345

    1999

  9. [16]

    Rauschenbeutel, G

    A. Rauschenbeutel, G. Nogues, S. Osnaghi, P. Bertet, M. Brune, J. M. Raimond and S. Haroche, Science 288 (2000) 2024

    2000

  10. [17]

    Alessandro Ferraro, Matteo G. A. Paris, Maria Bondani, Alessia Allevi, Emiliano Puddu and Alessandra Andreoni, J. Opt. Soc. Am. B 21 (2004) 1241

    2004

  11. [18]

    C. A. Sackett, D. Kielpinski, B. E. King, C. Langer, V. Me yer, C. J. Myatt, M. Rowe, Q. A. Turchette, W. M. Itano, D. J. Wineland and C. Monroe, Nature 404 (2000) 256

    2000

  12. [19]

    J.-W. Pan, M. Daniell, S. Gasparoni, G. Weihs and A. Zeil inger, Phys. Rev. Lett. 86 (2001) 4435

    2001

  13. [20]

    Maria Bondani, Alessia Allevi, Emiliano Puddu, Alessa ndra Andreoni, Alessandro Ferraro and Matteo G. A. Paris, Opt. Lett. 29 (2004) 180

    2004

  14. [21]

    Br iegel and Jian-Wei Pan, Nature 430 (2004) 54

    Zhi Zhao, Yu-Ao Chen, An-Ning Zhang, Tao Yang, Hans J. Br iegel and Jian-Wei Pan, Nature 430 (2004) 54

    2004

  15. [23]

    de Souza Dutra, Ann

    A. de Souza Dutra, Ann. Phys 321 (2006) 1092

    2006

  16. [24]

    Gell-Mann and Y

    M. Gell-Mann and Y. Ne’eman, The Eightfold Way, W. A. Ben jamin (1964). Also see https://en.wikipedia.org/wiki/Gell-Mann matrices. 11

    1964

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.