{"id":"97d0918e-65fa-4035-9b2d-1701685d0a35","arxiv_id":"1908.05231","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Numerical simulations show that two equal-amplitude drives on an optomechanical oscillator lock its stabilized average energy onto discrete, noise-robust levels.","lead":"This paper reports that a classical mechanical oscillator in an optomechanical cavity, driven by two equally strong lasers at matched frequencies, can settle onto one of several discrete amplitude levels, resembling energy quantization. The authors argue these 'classical energy levels' form through a new synchronization mechanism and could enable sensitive switching detectors or multistable mechanical memories.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Numerically observed energy levels are not verified as true attractors: Appendix B admits that for n≥4 a different computation precision leads to a different final level, so the discrete-level claim may reflect solver artifacts rather than invariant dynamics of Eq. (1).","rationale":"The reader's conditional verdict is appropriate. The paper asserts a new synchronization-induced discretization of average mechanical energy, but the evidence is entirely numerical: Fig. 2A-2B, Fig. 4, and Appendices B-D. The strongest independent support is that the lower levels (n≤3) are stated to be precision-independent, and Appendix C checks the single-mode ansatz by comparing the linearized Eq. (C2) with the full Eq. (1). That is a real consistency check, but it is conditional on An already being read from the same numerical contours, and it does not prove existence of the levels as dynamical invariants. The manuscript's own Appendix B is the crux: for n≥4, changing computation precision changes the level reached. This is not disqualifying by itself, because multistable systems can select different basins under small numerical perturbations, but it places the burden on the authors to show that the observed plateaus are true periodic orbits with finite basins. No code, data, or basin/orbit analysis is provided, so an independent reproduction cannot currently be performed. Therefore the central claim is plausible but conditional: if the proposed concrete test shows that the level shifts with solver or precision, or that no periodic orbit with the claimed amplitude exists, the claim should be rejected; if it passes, the phenomenon is established at least for the tested drive range, though the extrapolation to arbitrary drive intensity would still need separate support.","tokens_in":13010,"tokens_out":5088,"duration_ms":55358,"concrete_test":"At the Fig. 4B operating point E=(10^8+6.01×10^-4)κ, integrate Eq. (1) from the stated initial conditions using (i) at least two independent high-order ODE solvers (e.g., an adaptive RK with tight tolerance and a fixed-step variable-precision integrator) and (ii) 32-, 64-, and 128-digit precision. Record the final level ⟨Em⟩ and the stabilized waveform. Then test whether the final waveform is a true periodic orbit by checking self-recurrence Xm(t+T)=Xm(t) and Pm(t+T)=Pm(t) with T=2π/ωm to 1e-10 relative tolerance over many periods, and compute a small basin around it. If the selected level changes across solvers/precisions at fixed physical tolerance, or if no period-2π/ωm orbit with the claimed A_n exists, the discrete-level structure is not established as a property of Eq. (1).","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim requires that the plateaus in ⟨Em⟩ correspond to genuine invariant trajectories of Eq. (1), independent of the numerical method. The paper's own Appendix B states this condition is not met in the regime where the claim is most dramatic: 'in the simulations involving the higher energy levels (n≥4), a different computation precision will lead to a different energy level being reached in the end.' The authors reinterpret this as sensitivity to initial conditions, since different precisions are 'realized by the different iteration step sizes.' That interpretation is plausible only if the set of candidate levels is independently known to be a set of stable limit cycles of Eq. (1). But no such independent verification appears: the table in Appendix C reads An and dn from the same numerical contours it is meant to validate, and the linearized Eq. (C2) is checked after those amplitudes are already inserted, so it cannot establish existence. With no code, data, or basin/orbit analysis, the possibility remains that the n≥4 levels, or the specific level reached, are artifacts of a particular integrator, step size, and integration horizon. Since the abstract extrapolates the claim to 'whatever drive intensity beyond a threshold,' the weakest load-bearing step is exactly the identification of numerical plateaus with dynamical invariants.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports numerical evidence that a classical optomechanical oscillator driven by two equal-amplitude fields, one red-detuned by the mechanical frequency and one resonant, can synchronize onto a discrete set of stabilized oscillatory states whose time-averaged mechanical energy ⟨Em⟩ takes discrete values resembling energy levels. The central claim is that the amplitude and phase of the mechanical motion become locked to one of several discrete levels, that the levels are robust to noise after stabilization, and that the transition between levels is sensitive to initial conditions and drive amplitude. The paper presents simulations of Eq. (1), a single-mode approximation for the locked motion, a comparison with linearized cavity dynamics in Appendix C, and numerical robustness checks under drive and thermal noise in Section VII.","tokens_in":13230,"tokens_out":3123,"duration_ms":33996,"significance":"If the central claim is correct, the paper reports a genuinely striking classical analogue of discrete energy levels, with potential implications for nonlinear dynamics, synchronization, and precision measurement. The manuscript gives credit for a number of careful numerical checks: the single-mode approximation is compared against the full nonlinear dynamics in Fig. C1, the noise robustness of the stabilized levels is demonstrated in Fig. 5, and the dependence on the two drives is explored systematically in Appendices C and D. However, the significance is conditional because the existence of the discrete levels, especially for n≥4, is established only by numerical integration of the ordinary differential equations, with no independent verification that the plateaus are true invariant attractors of Eq. (1). The paper also contains no analytic derivation of the level positions or of the claimed power law, so the result currently rests on numerical observation.","major_comments":[{"comment":"The admitted precision dependence for the higher energy levels is a load-bearing gap. Appendix B states that for n≥4, a different computation precision will lead to a different energy level being reached, and the authors reinterpret this as sensitivity to initial conditions caused by different iteration step sizes. That interpretation is plausible only if the set of candidate levels is independently known to consist of stable limit cycles of Eq. (1). The paper provides no Floquet stability analysis, no basin-of-attraction computation, and no convergence study demonstrating that the level values are invariant under changes of integrator, tolerance, and integration horizon. Since the abstract extrapolates the phenomenon to \"whatever drive intensity beyond a threshold,\" the n≥4 levels are not yet established as properties of the differential equations rather than of a particular numerical solver.","section":"Appendix B"},{"comment":"The linearized validation in Appendix C is partly circular. Table C1 reads the amplitudes An and displacements dn from the contours of the same numerical simulations that the linearized Eq. (C2) is meant to validate, and the agreement in Fig. C1 then confirms the single-mode ansatz but not the existence or discreteness of the levels. To substantiate the central claim, the candidate amplitudes should be obtained independently, for example by solving the algebraic conditions for a periodic orbit of the full system or by continuing the n=1 branch in drive amplitude and checking for bifurcations. Without such an independent determination, the plateaus in ⟨Em⟩ could be an artifact of the way the simulations are initialized and integrated.","section":"Appendix C"},{"comment":"The power law ⟨Em⟩(n)∼n^2.2 is presented without derivation, uncertainty estimate, or a statement of the fitted range. With only five levels, irregular horizontal spacing, and overlapping levels as shown in Figs. 2B1 and 2B2, the exponent is not a quantitatively supported claim. Please provide the fit residuals, the dependence of the exponent on system parameters, or an analytic derivation; otherwise the claim should be explicitly labeled as an empirical observation over the computed range.","section":"Section III"},{"comment":"The claim that the sensitivity to initial conditions \"has nothing to do with chaos\" is supported only by the observation that final differences appear bounded in two example trajectories. A bounded difference over a finite integration time does not exclude transient chaos or a weak positive Lyapunov exponent, especially in a regime where the manuscript itself reports extreme sensitivity to precision. A quantitative Lyapunov-exponent calculation or a return-map analysis of the stabilized motion is needed to make the non-chaos claim persuasive.","section":"Section VI"}],"minor_comments":[{"comment":"The abstract contains typographical errors: \"continues values\" should be \"continuous values,\" and \"against intuition\" would read more naturally as \"contrary to intuition.\"","section":"Abstract"},{"comment":"The caption contains a typo: \"different enenrgy level\" should be \"different energy level.\"","section":"Fig. 4 caption"},{"comment":"In the sentence \"after magnifying the scales on the horizontal axis as in Figs. B1 and B2,\" the references should be to Figs. 2B1 and 2B2 in the main text.","section":"Section III"},{"comment":"The reference to \"Fig. S2\" in the caption of Fig. D2 appears to be a leftover from a supplementary file and should refer to Fig. D2 or another figure in the paper.","section":"Appendix D"},{"comment":"The formal time-ordered exponential expansion in Eq. (C3) is not used in the comparison; the text says the linearized equations are integrated numerically. The expansion could be removed or its purpose clarified to avoid confusion.","section":"Appendix C"},{"comment":"The manuscript would benefit from a description of the numerical integration method, tolerances, and convergence criteria, as well as a statement on code or data availability, to make the numerical results reproducible.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"I see no ethical concerns about the manuscript itself. The main issue is whether the reported discrete levels are genuine invariant sets of Eq. (1) or numerical artifacts; the authors' own Appendix B makes this the central open point. The paper is potentially interesting, but the numerical evidence needs to be complemented by a convergence study, a stability analysis, or an independent determination of the level amplitudes before the claims can be accepted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Alex—quick read of He et al. 1908.05231. The claim: with two equal-strength drives (one red-detuned by ω_m, one resonant), a classical optomechanical oscillator stabilizes onto one of several discrete average-energy plateaus, like energy levels, robust to noise once settled. That is genuinely new as far as I can tell; the cited literature on two-drive optomechanics uses one strong and one weak field, and single-drive multistability doesn't give these sharp plateaus.\n\nWhat's done well: the paper is honest about the numerical basis. They check the single-mode ansatz X_m(t)=A sin(ω_m t)+d against full simulations, compare the linearized cavity equations to the nonlinear ones, and get good agreement for the first few levels. The lower levels (n=1-3) look solid—the plateaus are crisp across a range of drive amplitudes, and the noise robustness is demonstrated with delayed perturbations.\n\nThe soft spots are real but not fatal. First, no analytic derivation of why the levels exist or why the energies scale like n^2.2. The 'predictions' in Appendix C use amplitudes A_n read from the same simulations they are meant to validate—that's circular, not an independent check. Second, Appendix B admits that for n≥4 different computation precision leads to different final levels. The authors interpret this as sensitivity to initial conditions, which is plausible if the set of level positions is actually invariant. They claim that any stabilized result lands on one of the fixed levels regardless of precision—that's the key test, and it partially addresses the concern. But without code, data, or a basin/orbit analysis, I can't fully rule out that the higher levels are integrator artifacts. The extrapolation 'for whatever drive intensity beyond a threshold' is also beyond the simulated range.\n\nWho's this for? Nonlinear dynamics and optomechanics researchers, especially anyone looking for switchable multistable states. It deserves serious peer review, but conditional on the authors providing code/data and ideally an analytic or experimental test of the level structure. I wouldn't cite it yet as established fact, but I'd want it in the discussion.\n\nSend it to a competent referee, but the authors should be asked to share the numerics.","headline":"A novel numerical observation of discrete mechanical energy levels in a doubly driven optomechanical system, but the evidence is almost entirely computational and the higher-level claims need independent reproduction and an analytic anchor.","tokens_in":13757,"tokens_out":2558,"would_cite":false,"duration_ms":26274,"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":"Two equally strong, properly detuned drives can lock a classical mechanical oscillator onto discrete, noise-robust energy levels.","keywords":["synchronization","mode locking","optomechanics","classical energy levels","nonlinear dynamics","drive-induced stabilization","amplitude locking","noise robustness"],"falsifier":"Take the parameters of Fig. 1B ($\\omega_m=50\\kappa$, $g/\\sqrt{2}=10^{-5}\\kappa$, $\\gamma_m=10^{-5}\\kappa$) and sweep $E$ continuously through the claimed threshold $E\\approx 5\\times 10^5\\kappa$ up to $E\\approx 10^8\\kappa$, using two independent high-precision integrators with step sizes differing by, say, a factor of two. If the stabilized time-averaged mechanical energy $\\langle E_m\\rangle$ is not constant over finite intervals but instead varies continuously, or if the set of stabilized values depends on the integrator's step size rather than converging to the tabulated levels $A_n$, then the discrete energy levels are numerical artifacts. The same check can be performed experimentally by measuring $\\langle E_m\\rangle$ versus drive amplitude in a high-finesse cavity with the two matched drives.","tokens_in":12774,"feed_emoji":"🔒","tokens_out":7374,"duration_ms":76080,"temperature":0.7,"pith_summary":"This paper aims to show that a classical macroscopic mechanical oscillator can acquire discrete energy levels through synchronization, contrary to the intuition that classical energies vary continuously. The mechanism is two equal-amplitude drives on an optomechanical cavity: one red-detuned by the mechanical frequency (the cooling field) and one resonant. When both drives exceed a threshold, the oscillator's stabilized oscillation amplitude $A$ freezes at a series of discrete values $A_n$, so the time-averaged mechanical energy jumps between plateaus $\\langle E_m\\rangle^{(n)}\\approx \\tfrac{1}{2}A_n^2$ even as drive amplitude is swept continuously. The selected level is highly sensitive to initial conditions and early fluctuations, but once stabilized the locked oscillation is robust to noise and drive fluctuations. These results matter because they suggest that quantization-like energy levels could be produced purely classically and detected in existing optomechanical setups.","feed_headline":"Two synchronized drives freeze an oscillator on discrete energy levels","feed_subtitle":"Matched drives freeze a macroscopic oscillator's average energy onto fixed plateaus.","key_machinery":"The central object is the pair of cooperating drives on the optomechanical system described by Eq. (1): a cooling field detuned by $\\omega_m$ from the cavity and a resonant field, with equal amplitudes $E_1=E_2$. What carries the argument is the nonlinear coupling $\\tfrac{\\sqrt{2}}{4}gX_m(X_c^2+P_c^2)$ between the cavity field quadratures and the mechanical displacement. Under matched frequencies this nonlinearity turns the resonant drive's intensity into a staircase response of the mechanical amplitude, while the cooling field lowers the resulting levels and locks their phases; the single-mode approximation $X_m(t)=A\\sin(\\omega_m t)+d$ (Eq. C1) then lets the discrete $A_n$ be read off the stabilized trajectories and connects each level to the average energy $\\langle E_m\\rangle^{(n)}\\approx \\tfrac{1}{2}A_n^2$.","core_discovery":"Under the frequency-matching condition $\\Delta_1=\\omega_m$ and $\\Delta_2=0$ with equal drive amplitudes $E_1=E_2=E$, the mechanical displacement approaches a single-frequency locked oscillation $X_m(t)=A\\sin(\\omega_m t+\\varphi)+d$. The amplitude $A$ takes one of a discrete set of values $A_n$ rather than growing continuously with $E$; the time-averaged mechanical energy therefore sits on levels $\\langle E_m\\rangle^{(n)}\\approx \\tfrac{1}{2}A_n^2$, with the empirical power law $\\langle E_m\\rangle^{(n)}\\sim n^{2.2}$. The cavity field develops a characteristic oscillation pattern whose number of peaks in half a period identifies the mechanical level, and all frequency components of the two oscillators are phase-synchronized ($n{:}m$ synchronization). The resonant drive alone creates the discrete amplitude response through nonlinear saturation; adding the equal cooling drive lowers the levels and synchronizes phases across different $E$. The level ultimately reached depends sensitively on initial conditions (a difference of $10^{-18}$ in initial energy can change the final level), but the final locked state is immune to noise applied after stabilization. The paper's claim is that this discrete-level structure is a genuine property of the underlying differential equations, validated by accepting only numerical results that survive refinement of precision.","pith_inferences":["Inference: The irregular, nearly random pattern of which level is selected as the drive amplitude changes suggests that basin boundaries between neighboring levels are complicated, possibly fractal; if so, this two-drive system could serve as a deterministic yet unpredictable level selector, a candidate physical source of random bits.","Inference: The same mechanism may generalize beyond optomechanics to any two nonlinearly coupled oscillators with one resonant and one cooling-like drive; a direct test would be to look for amplitude plateaus in driven micromechanical or optoelectronic oscillator experiments.","Inference: The empirical power law $\\langle E_m\\rangle^{(n)}\\sim n^{2.2}$ is close to quadratic; if higher levels tighten toward $n^2$, the classical level spacings would mimic a harmonic oscillator, a connection the paper does not make.","Inference: The robustness of stabilized levels suggests a possible control protocol: apply a short pulse during the transient to select a desired level, then let the system settle; the paper's square-pulse tests hint at this but do not demonstrate deterministic level selection."],"forward_implications":["If the claim is right, a large classical object can have discrete, noise-robust energy levels, so quantization-like behavior can be mimicked by classical nonlinear dynamics without invoking quantum mechanics.","The mechanical energy level can be read out directly from the cavity field's oscillation pattern (peak count in a half period), giving a practical non-invasive level indicator.","Because the final level is highly sensitive to perturbations during the transient but immune to them after stabilization, the system could act as a sensor for small early disturbances or tiny differences in initial conditions.","The strict condition $\\Delta_1=\\omega_m$ with $\\Delta_2=0$ turns the two-drive setup into a way to measure the mechanical frequency precisely: energy-level locking occurs only at exact matching.","The first energy level is predicted at effective cooling intensity $J\\approx 0.1$, within reach of current optomechanical systems, so the phenomenon should be testable now."],"supporting_citations":[{"why":"Supplies the standard optomechanical model, the damping hierarchy $\\gamma_m\\ll\\kappa$, and the quadratic coupling terms used in Eq. (1).","marker":"[10]"},{"why":"Defines synchronization and mode locking as the phenomena being generalized from frequency and phase locking to amplitude locking.","marker":"[2]"},{"why":"Shapiro steps provide the canonical mode-locking precedent that motivates asking whether oscillation amplitude can also be locked to fixed values.","marker":"[8]"},{"why":"Establishes the previous two-drive optomechanics setting with unequal amplitudes, against which the equal-amplitude cooperating drives are the new situation.","marker":"[16]"},{"why":"Provides the effective cooling intensity $J=gE/\\omega_m$ used to argue that the first energy level is experimentally attainable at $J\\approx 0.1$.","marker":"[26]"},{"why":"Gives prior numerically studied radiation-pressure nonlinearity where the stabilized amplitude varies continuously, highlighting the contrast with the discrete plateaus reported here.","marker":"[28–33]"}],"fun_headline_variants":["Classical oscillator locked to discrete energy levels by twin drives","Two drives trap a classical oscillator on quantum-like energy steps","Matched drives freeze a macroscopic oscillator onto level plateaus","Synchronized equal drives quantize a classical oscillator's motion"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole phenomenon rests on the numerical solutions of Eq. (1) faithfully representing the dynamics; for levels $n\\ge 4$ the paper reports that changing computation precision can switch which level is reached, so if the discrete plateaus are artifacts of the solver rather than true properties of the equations, the central claim fails.","fun_headline_variants_meta":{"raw":{"variants":["Classical oscillator locked to discrete energy levels by twin drives","Two drives trap a classical oscillator on quantum-like energy steps","Matched drives freeze a macroscopic oscillator onto level plateaus","Synchronized equal drives quantize a classical oscillator's motion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000818,"raw_usage":{"total_tokens":3584,"prompt_tokens":950,"completion_tokens":2634,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":566,"completion_tokens_details":{"reasoning_tokens":2566}},"tokens_in":566,"tokens_out":2634,"duration_ms":18039,"temperature":1.0,"reasoning_tokens":2566,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:18:03.475067+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the parameters of Fig. 1B ($\\omega_m=50\\kappa$, $g/\\sqrt{2}=10^{-5}\\kappa$, $\\gamma_m=10^{-5}\\kappa$) and sweep $E$ continuously through the claimed threshold $E\\approx 5\\times 10^5\\kappa$ up to $E\\approx 10^8\\kappa$, using two independent high-precision integrators with step sizes differing by, say, a factor of two. If the stabilized time-averaged mechanical energy $\\langle E_m\\rangle$ is not constant over finite intervals but instead varies continuously, or if the set of stabilized values depends on the integrator's step size rather than converging to the tabulated levels $A_n$, then the discrete energy levels are numerical artifacts. The same check can be performed experimentally by measuring $\\langle E_m\\rangle$ versus drive amplitude in a high-finesse cavity with the two matched drives.","supporting_citations":[{"cited_title":"Aspelmeyer, T","cited_arxiv_id":null,"evidence_quote":"Supplies the standard optomechanical model, the damping hierarchy $\\gamma_m\\ll\\kappa$, and the quadratic coupling terms used in Eq. (1)."},{"cited_title":"Pikovsky, M","cited_arxiv_id":null,"evidence_quote":"Defines synchronization and mode locking as the phenomena being generalized from frequency and phase locking to amplitude locking."},{"cited_title":"Shapiro, Josephson Currents in Superconducting Tunneling: The Effect of Microwaves and Other Observations","cited_arxiv_id":null,"evidence_quote":"Shapiro steps provide the canonical mode-locking precedent that motivates asking whether oscillation amplitude can also be locked to fixed values."},{"cited_title":"Kronwald, F","cited_arxiv_id":null,"evidence_quote":"Establishes the previous two-drive optomechanics setting with unequal amplitudes, against which the equal-amplitude cooperating drives are the new situation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the effective cooling intensity $J=gE/\\omega_m$ used to argue that the first energy level is experimentally attainable at $J\\approx 0.1$."}],"review_version":1}