{"id":"72806fe6-4acf-436e-914b-2d08380dcb7b","arxiv_id":"2608.03578","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"In a Coulomb-coupled electro-optomechanical ring cavity, the Coulomb interaction is required to produce entanglement, and an intracavity optical parametric amplifier can significantly enhance both bipartite and tripartite entanglement.","lead":"This paper studies a triangular optical cavity in which light pushes on mechanical mirrors, two of which are electrically coupled, and an optical parametric amplifier boosts the coupling. It reports numerical evidence that the Coulomb coupling creates both two- and three-way entanglement, and that the amplifier can dramatically increase it, at the cost of a smaller stable operating range.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The numerical results are built on an invalid steady-state amplitude: Eq. (19) treats α* as α, so the α used in the linearized couplings does not solve Eq. (16); all Section 5 entanglement curves are suspect.","rationale":"The reader's weakest assumption identified the α-real simplification as the main risk; I agree that this is the load-bearing point, but the problem is sharper than a phase-convention caveat. Eq. (19) is not merely a phase-restricted version of the correct steady state; it is algebraically inconsistent with Eq. (16) whenever G≠0 and Δ≠2G sinθ. The exact complex solution has a different magnitude and phase, and the linearized drift matrix (22) omits the imaginary-part couplings that a complex α would generate. Since every figure in Section 5 is produced from this drift matrix and this α, the specific quantitative maxima quoted in the abstract and results (EM1M2≈0.52, RM1M2C≈0.04) are unsupported. The qualitative claim that Coulomb coupling is necessary and that OPA can enhance correlations may survive a corrected calculation, but the current manuscript's numbers cannot be trusted. The paper does have positive features: a standard Gaussian-covariance pipeline, stability checks at every parameter point, and an explicit statement that SymPy was used for the commutation algebra. However, symbolic verification of the algebra does not catch the error because the incorrect steady-state equation was coded. The missing geometric angle cos²φ and the absence of code/data are secondary but should also be fixed. Therefore I recommend REJECT for the current version, with the expectation that a revised manuscript solving the exact complex steady state and rerunning the full linearization could be reconsidered.","tokens_in":14490,"tokens_out":15425,"duration_ms":136780,"concrete_test":"Recompute Fig. 2(a) from the exact fixed point: solve (κ+iΔ)α − 2G e^{iθ}α* = EL for complex α, set G1,3 = √2 g1,3 cos²φ α (complex), build the 8×8 drift matrix including both Re(α)δx and Im(α)δy couplings, and solve AV+VA^T=−D for the same parameters (θ=π, λ=0.9ωm, PL=60 mW, T=1 mK). If the peak EM1M2 at G/κ=0.4 differs from ≈0.471 by more than 10%, or if the stability boundary shifts, the quantitative claims in the abstract and Section 5 are unsupported.","verdict_should_be":"REJECT","load_bearing_attack":"Section 2.3.2 solves the steady-state of Eq. (16) by writing α = EL/[κ−2G cosθ + i(Δ−2G sinθ)] (Eq. 19). This is the solution of ((κ+iΔ)−2G e^{iθ})α = EL, i.e. it replaces α* by α. The actual steady-state condition is (κ+iΔ)α − 2G e^{iθ} α* = EL. Solving it with α = u+iv gives u = EL(κ+2G cosθ)/(κ²+Δ²−4G²) and v = −EL(Δ−2G sinθ)/(κ²+Δ²−4G²), which is not Eq. (19). For θ=π, for example, Eq. (19) gives α = EL[(κ+2G)−iΔ]/[(κ+2G)²+Δ²], whereas the exact solution is α = EL[(κ−2G)−iΔ]/[κ²+Δ²−4G²]; the two agree only for G=0 or Δ=2G sinθ. The paper then forms G1,3 = √2 g1,3 cos²φ α and, in Section 2.3.3, takes α real, omitting the δy-quadrature couplings that a complex α produces in the linearized optomechanical interaction. The drift matrix (22) therefore is not the linearization of Eq. (16) at its fixed point. Because the Section 5 parameter scans vary Δ and θ independently without enforcing Δ=2G sinθ, the covariance matrices and quantitative entanglement maxima (EM1M2≈0.52, RM1M2C≈0.04) are computed from a state that is not a steady state. This is an internal algebraic inconsistency, independent of any external consensus.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript proposes a hybrid electro-optomechanical ring-cavity system containing two Coulomb-coupled charged mechanical resonators, a third uncharged mechanical resonator, and an intracavity optical parametric amplifier. The authors derive the rotating-frame Hamiltonian, write nonlinear quantum Langevin equations, linearize around steady-state solutions, and analyze Gaussian fluctuations through the covariance matrix and the Lyapunov equation. Bipartite entanglement is quantified by logarithmic negativity and tripartite entanglement by the minimum residual contangle, computed as functions of detuning, OPA gain and phase, Coulomb coupling strength, laser power, and temperature. The central claims are that the Coulomb interaction is indispensable for entanglement generation and that the OPA dramatically enhances all bipartite and tripartite correlations, subject to a stability trade-off.","tokens_in":14897,"tokens_out":11873,"duration_ms":99377,"significance":"If the quantitative results were correct, the paper would provide a useful parameter map for entanglement engineering in a hybrid optomechanical platform and a concrete example of OPA-based control of Gaussian multipartite entanglement. The work has genuine strengths: the optomechanical coupling is derived with oblique incidence, the covariance-matrix formalism is applied in a standard way, stability is checked eigenvalue-by-eigenvalue, and the parameter values are physically motivated. However, the numerical foundation of the paper is compromised by an algebraic error in the steady-state optical amplitude, and because every numerical curve in Section 5 is generated from that amplitude, the claimed enhancement factors and quantitative maxima are not currently supported.","major_comments":[{"comment":"The steady-state optical amplitude in Eq. (19) is not the solution of the steady-state condition obtained from Eq. (16). Setting ⟨ȧ⟩ = 0 in Eq. (16) gives (κ + iΔ)α − 2G e^{iθ} α* = E_L, whereas Eq. (19) solves (κ + iΔ − 2G e^{iθ})α = E_L, i.e., it replaces α* by α. The correct solution is α = E_L[(κ + 2G cosθ) − i(Δ − 2G sinθ)]/(κ² + Δ² − 4G²), which differs from Eq. (19) except in the special cases G = 0 or Δ = 2G sinθ. Because Section 4 constructs G1,3 from this α and every subsequent stability check, Lyapunov solution, and entanglement curve depends on it, all numerical results in Section 5 are computed for a point that is not a steady state of the model. The authors must correct Eq. (19), rerun the numerical pipeline, and verify whether the claimed maxima (e.g., EM1M2 ≈ 0.52 and RM1M2C ≈ 0.04) survive.","section":"§2.3.2, Eq. (19)"},{"comment":"The linearization in Section 2.3.3 assumes α can be taken real 'without loss of generality' and defines G1,3 = √2 g1,3 cos²φ α, while Section 4 states G1,3 = √2 g1,3 |α|. For a genuinely complex α, the optomechanical interaction linearizes into both δx and δy quadrature couplings, and the drift matrix entries involving G1,3 change; the δy couplings are absent from Eq. (20) and Eq. (22). A phase rotation of the cavity mode can make α real, but only at the cost of also rotating the effective detuning and the OPA phase. Since Section 5 scans Δ and θ independently as if no such rotation had been performed, the drift matrix used in the numerics is not the linearization of Eq. (16) at the actual steady state. The authors should either carry the complex α through the linearization or define and apply the full phase rotation consistently before scanning parameters.","section":"§2.3.3 and §4"},{"comment":"Because of the steady-state error, the quantitative claims in Section 5 are not reliable in their present form. This includes the reported maxima EM1M2 ≈ 0.52 at λ = 0.95ωm and Δ = 0.65ωm, RM1M2C ≈ 0.04 at λ = 0.94ωm and Δ = 0.31ωm, and the OPA-enhancement statements for EM3C. In addition, the text in Fig. 5(a) reports a peak logarithmic negativity of 0.81 for G/κ = 0.3, which is inconsistent with the nearby values in Figs. 5 and 6 for the same quantity and is likely a typographical error. After correcting the steady state, all figures, peak values, and the conclusions drawn from them must be regenerated and re-evaluated.","section":"§5"}],"minor_comments":[{"comment":"The formula for ϱ_{r|st} is written as a minimum eigenvalue, but the residual contangle requires the minimum symplectic eigenvalue of the partially transposed covariance matrix; please correct Eq. (30) and specify the quadrature ordering for the partial-transposition matrices P_{r|st}, P_{s|rt}, and P_{t|sr}.","section":"§3.2, Eq. (30)"},{"comment":"The reported peak EM3C ≈ 0.81 at G/κ = 0.3 should be checked against Fig. 6, where the same quantity with G/κ = 0.3 is approximately 0.08; either the value, the axis label, or the figure caption contains an error.","section":"§5, Fig. 5(a)"},{"comment":"The geometric angle φ is never given a numerical value, although every optomechanical coupling strength G1,3 contains cos²φ; please state the value used or explain why the results are insensitive to it.","section":"§5"},{"comment":"The symbol Δ denotes the effective detuning in Eq. (19) but is used as the scan variable in Section 5; please clarify whether the plotted Δ/ωm is the effective detuning including the mechanical steady-state shifts or the bare detuning, and if necessary redefine the variable to avoid ambiguity.","section":"§2.3.2–§5"},{"comment":"Some cross-references are imprecise, e.g., the text refers to 'Fig. 6(a)' although the caption of Fig. 6 does not label panels; please correct these references.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The decisive issue is the incorrect steady-state amplitude in Eq. (19), which invalidates all quantitative results as currently presented. This is fixable within the scope of a revision: the authors must correct the steady-state solution, rerun the stability and Lyapunov calculations, and regenerate every figure and quantitative claim in Section 5. I also note that the introduction cites a long series of nonlinear-wave papers (Refs. 37–45) that appear tangential to the optomechanical content; the authors should either connect these references explicitly to their model or remove them."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hi —\n\nBottom line: the stress-test note is right, and the issue is load-bearing. Eq. (19) is not the solution of Eq. (16) for a real drive. Solving (κ+iΔ)α − 2G e^{iθ} α* = E_L gives α = E_L(κ+2G cosθ − i(Δ−2G sinθ))/(κ²+Δ²−4G²), which differs from Eq. (19) at second order in G/κ. The paper never imposes the condition that would make the two coincide, and the scans in Section 5 vary Δ and θ independently. The linearization also uses G1,3 = sqrt(2) g1,3 cos²φ |α| with only δx couplings; that is a valid simplification only after a phase rotation that changes the OPA phase and the drive phase. As written, the drift matrix is not the linearization at the fixed point, so the covariance matrices and all Section 5 curves are not steady-state results.\n\nWhat the paper does well: the architecture is new — two Coulomb-coupled mechanical oscillators plus a third oscillator in a triangular ring with an intracavity OPA — and the qualitative logic is plausible: Coulomb coupling creates entanglement and the OPA amplifies it. The Hamiltonian derivation, the covariance-matrix formalism, and the stability analysis are standard and clearly described. The symbolic-algebra checks are a good practice.\n\nThe other problems are minor in comparison but worth listing. The geometric angle φ is never given numerically, so cos²φ is an unspecified scale factor in every plot. There is a likely typo in the Fig. 5(a) discussion (a peak of 0.81 that should be around 0.081) and a caption/label mismatch in Fig. 9(b). No code or data are provided. The citation list includes several nonlinear-wave papers that are not used in the derivation; they read as padding.\n\nMy recommendation: do not publish this version. The authors need to solve the complex fixed point correctly or state the phase convention explicitly, then re-run the parameter scans and report the actual numbers. If they do, the paper could be a reasonable contribution for a specialized journal. A referee familiar with OPA optomechanics should be asked to check the fixed-point algebra before anything else. But in its current state the central quantitative claims do not survive.","headline":"Interesting new optomechanical architecture, but the numerical results use a steady-state amplitude that does not solve the stated equations; the quantitative entanglement claims are unreliable as written.","tokens_in":15446,"tokens_out":9930,"would_cite":false,"duration_ms":88303,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["03.67.-a","03.67.Bg","42.50.Lc","42.50.Wk","42.65.Yj"],"model":"deepseek-v4-flash","headline":"This paper claims that Coulomb coupling generates entanglement and an optical parametric amplifier amplifies it in an electro-optomechanical ring cavity.","keywords":["optomechanics","ring cavity","bipartite entanglement","tripartite entanglement","parametric nonlinearity","nonlinear optics","Coulomb coupling","optical parametric amplifier"],"falsifier":"Set the Coulomb coupling to zero ($\\lambda=0$) and run the same covariance-matrix calculation: the paper predicts every entanglement measure is exactly zero for all detunings and gains; observing a nonzero logarithmic negativity in that case would refute the claim that the Coulomb interaction is indispensable.","tokens_in":14309,"feed_emoji":"🔗","tokens_out":8308,"duration_ms":66346,"temperature":0.7,"pith_summary":"This paper claims that in a triangular electro-optomechanical ring cavity—two charged mechanical resonators coupled by the Coulomb force, a third uncharged mechanical oscillator, and an optical parametric amplifier (OPA)—entanglement is created only by the Coulomb interaction, and the OPA then acts as a tunable amplifier for that entanglement. Working with the linearized quantum Langevin equations and Gaussian states, the authors compute steady-state covariance matrices and quantify bipartite entanglement by logarithmic negativity $E_N$ and tripartite entanglement by the minimum residual contangle $R_{\\min}^{\\tau}$. They report that increasing the OPA gain $G$ and choosing its phase $\\theta$ can push the mechanical–mechanical entanglement from roughly $0.29$ to about $0.52$, while tripartite measures grow from about $0.018$ to $0.04$, with the caveat that high gain or strong Coulomb coupling shrinks the stable operating window. If correct, the result identifies two practical control knobs—charge voltage and parametric drive—for distributing stronger quantum correlations across mechanical and optical modes in a single hybrid device.","feed_headline":"OPA gain lifts mechanical entanglement to 0.52","feed_subtitle":"Two Coulomb-coupled resonators supply the correlations; a parametric amplifier turns them into a tunable quantum link.","key_machinery":"The load-bearing object is the $8\\times 8$ drift matrix $A$ of the linearized quantum Langevin equations, written in the fluctuation basis $(\\delta q_1,\\delta p_1,\\delta q_2,\\delta p_2,\\delta q_3,\\delta p_3,\\delta x,\\delta y)$. Its entries contain the Coulomb coupling $\\lambda$, the linearized optomechanical couplings $G_1=\\sqrt{2}g_1\\cos^2\\phi\\,\\alpha$ and $G_3$, and the OPA terms $\\pm 2G\\cos\\theta$, $\\pm 2G\\sin\\theta$; stability requires all its eigenvalues to have negative real parts. For each stable parameter set, the steady-state covariance matrix $V$ solves the Lyapunov equation $AV+VA^T=-D$, and bipartite entanglement is read from the smallest symplectic eigenvalue of the partially transposed reduced state, while tripartite entanglement comes from the minimum residual contangle obeying the Coffman–Kundu–Wootters monogamy inequality.","core_discovery":"On the paper's own terms, the central discovery is a division of labour between two nonlinear couplings. The Coulomb interaction $\\lambda q_1 q_2$ between the two charged membranes is the indispensable generator: at $\\lambda=0$ every computed entanglement measure, including entanglement between the optical mode and the third mechanical oscillator, vanishes. The intracavity OPA, with gain $G$ and phase $\\theta$, does not create correlations from nothing but dramatically amplifies the ones the Coulomb coupling seeds; by scanning detuning $\\Delta$, gain, phase, Coulomb strength, laser power, and temperature, the paper maps out stable parameter regions and reports maxima such as $E_{M_1M_2}\\approx 0.52$ at $\\lambda=0.95\\,\\omega_m$ and $\\Delta=0.65\\,\\omega_m$, and $R_{M_1M_2C}^{\\min}\\approx 0.04$ at $\\lambda=0.94\\,\\omega_m$, $\\Delta=0.31\\,\\omega_m$. The dependence on $\\lambda$ is non-monotonic for the optical–mechanical bipartition, and all correlations degrade with temperature, so the paper's conclusion is that low-temperature operation and careful parameter selection are required to exploit the OPA enhancement.","pith_inferences":["A natural extension the paper does not pursue is to compute Gaussian quantum steering or Bell-nonlocality measures from the same covariance matrix; if the OPA amplifies logarithmic negativity, it likely also amplifies steering, which would make the device more directly useful for one-sided quantum communication.","Because the third oscillator M3 is not directly charged, the nonzero M3–cavity entanglement at intermediate $\\lambda$ suggests the Coulomb pair acts as a nonlinear mediator; this could be tested by replacing the Coulomb coupling with a different quadratic interaction and checking whether the same 'seeding plus OPA amplification' pattern survives.","The reported stability threshold at $G/\\kappa\\simeq 0.5$ is a sharp, testable prediction: an experiment should see the output field transition from a stationary state to an unstable or oscillatory regime as the parametric pump crosses that gain, at fixed $\\theta=\\pi$.","The paper fixes the geometric angle $\\phi$ in the equations but never specifies its value in the numerics; since the effective couplings scale as $\\cos^2\\phi$, a quantitative comparison with experiment requires $\\phi$ to be reported or measured. This is an editorial caution, not a claim in the paper."],"forward_implications":["At $\\lambda=0$ none of the predicted bipartite or tripartite entanglement survives, so any experiment aiming to realize this scheme must first establish a nonzero charge-mediated coupling between the two mechanical resonators.","The OPA converts a modest baseline ($E_{M_1M_2}\\simeq 0.29$ at $G=0$) into a substantially stronger correlation ($\\simeq 0.52$ near $G=0.5\\kappa$, $\\theta=\\pi$, $\\lambda=0.95\\omega_m$), giving experimentalists a single optical knob to tune entanglement strength without changing the mechanical setup.","Detuning $\\Delta$ controls where the maximum sits: increasing OPA gain shifts the optimal detuning downward, so the operating point must be re-optimized whenever $G$ or $\\theta$ is changed.","Thermal noise degrades every correlation measure, and the OPA only partially compensates; the paper's numbers therefore presuppose cryogenic operation (order 1 mK) for the reported magnitudes.","Strong coupling and high gain trade against stability: near $\\lambda\\simeq\\omega_m$ or $G/\\kappa>0.5$ the stable detuning window narrows, so maximizing entanglement requires operating close to an instability boundary."],"supporting_citations":[{"why":"Supplies the realistic Coulomb coupling value $\\lambda\\approx 0.32\\omega_m$ and the experimental parameter set used to argue the coupling is tunable and attainable.","marker":"[48]"},{"why":"Provides the Lyapunov-equation method linking the drift matrix to the steady-state covariance matrix, the backbone of all entanglement numbers.","marker":"[52]"},{"why":"Defines logarithmic negativity, the bipartite entanglement measure used for every two-mode subsystem.","marker":"[53]"},{"why":"Gives the PPT criterion that makes $\\varrho<1/2$ necessary and sufficient for bipartite entanglement in Gaussian states.","marker":"[55]"},{"why":"Defines the residual contangle and the monogamy framework used to quantify genuine tripartite entanglement.","marker":"[56]"},{"why":"Supplies the stability condition (all drift-matrix eigenvalues with negative real parts) that defines the allowed parameter regions.","marker":"[51]"},{"why":"Shows that an OPA can improve entanglement in cavity optomechanical systems, the precedent the paper extends to the Coulomb-coupled ring cavity.","marker":"[34]"},{"why":"Another prior demonstration of OPA-enhanced optomechanical entanglement, used as a baseline for the enhancement claim.","marker":"[35]"}],"fun_headline_variants":["Coulomb seeds, OPA amplifies: tunable ring-cavity entanglement","Without Coulomb coupling, no quantum link; OPA tunes it up","OPA boost to 0.52 for Coulomb-paired mechanical modes","Trade-off: max entanglement vs stable operation in ring cavity","Parametric drive amplifies Coulomb-generated multipartite entanglement"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's numbers depend on assuming the cavity field settles into a stable steady state whose phase can be rotated away; if that rotation changes the meaning of the detuning and amplifier phase, the reported maxima may not be physically reachable.","fun_headline_variants_meta":{"raw":{"variants":["Coulomb seeds, OPA amplifies: tunable ring-cavity entanglement","Without Coulomb coupling, no quantum link; OPA tunes it up","OPA boost to 0.52 for Coulomb-paired mechanical modes","Trade-off: max entanglement vs stable operation in ring cavity","Parametric drive amplifies Coulomb-generated multipartite entanglement"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000307,"raw_usage":{"total_tokens":1827,"prompt_tokens":1088,"completion_tokens":739,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":704,"completion_tokens_details":{"reasoning_tokens":648}},"tokens_in":704,"tokens_out":739,"duration_ms":6925,"temperature":1.0,"reasoning_tokens":648,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T14:50:12.175877+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Set the Coulomb coupling to zero ($\\lambda=0$) and run the same covariance-matrix calculation: the paper predicts every entanglement measure is exactly zero for all detunings and gains; observing a nonzero logarithmic negativity in that case would refute the claim that the Coulomb interaction is indispensable.","supporting_citations":[{"cited_title":"et al,Sci","cited_arxiv_id":null,"evidence_quote":"Supplies the realistic Coulomb coupling value $\\lambda\\approx 0.32\\omega_m$ and the experimental parameter set used to argue the coupling is tunable and attainable."},{"cited_title":"Vitali, S","cited_arxiv_id":null,"evidence_quote":"Provides the Lyapunov-equation method linking the drift matrix to the steady-state covariance matrix, the backbone of all entanglement numbers."},{"cited_title":"Vidal and R","cited_arxiv_id":null,"evidence_quote":"Defines logarithmic negativity, the bipartite entanglement measure used for every two-mode subsystem."},{"cited_title":"Simon,Phys","cited_arxiv_id":null,"evidence_quote":"Gives the PPT criterion that makes $\\varrho<1/2$ necessary and sufficient for bipartite entanglement in Gaussian states."},{"cited_title":"Adesso and F","cited_arxiv_id":null,"evidence_quote":"Defines the residual contangle and the monogamy framework used to quantify genuine tripartite entanglement."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the stability condition (all drift-matrix eigenvalues with negative real parts) that defines the allowed parameter regions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows that an OPA can improve entanglement in cavity optomechanical systems, the precedent the paper extends to the Coulomb-coupled ring cavity."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Another prior demonstration of OPA-enhanced optomechanical entanglement, used as a baseline for the enhancement claim."}],"review_version":2}