{"id":"60ed326b-a348-4647-91b2-ba59d9c9fc07","arxiv_id":"1908.08134","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Numerical evidence shows a quantum analogue of the Neimark-Sacker bifurcation in a driven open bosonic dimer: a change from a single-peaked to bagel-shaped asymptotic distribution together with Floquet eigenvalues approaching the unit circle.","lead":"A periodically driven open quantum dimer develops a doughnut-shaped state distribution, the quantum counterpart of a classical Neimark-Sacker bifurcation. The result adds a new type of transition to the small family of bifurcations already observed in open quantum systems.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No scaling analysis links the finite-N bagel onset to the mean-field Neimark-Sacker point, so the claimed quantum-classical correspondence is not established.","rationale":"The reader identified the semiclassical mean-field truncation as the weakest assumption. My concern sharpens this: the paper's own results show that the quantum system does not simply mirror the mean-field bifurcation at the largest simulated N. At N=500, the bagel appears at U=0.1 while the mean-field fixed point is still stable, and an extra quantum bifurcation occurs at U∈[0.6,0.7] without a classical counterpart. These discrepancies suggest that the mean-field reduction may not capture the relevant physics, or that the quantum signatures are finite-size effects that do not converge to the classical bifurcation. The paper's statement that the spectral gap decreases with N is not backed by a scaling analysis, and the D(U) curves in Fig. 6 are not used to extract a size-dependent bifurcation point. A concrete test is to simulate larger N (1000, 2000) and check whether the onset U*(N) approaches the mean-field value; if not, the central claim of a quantum Neimark-Sacker bifurcation analogous to the classical one lacks quantitative support. This does not invalidate the numerical observations but makes the interpretation conditional on additional validation, matching the reader's CONDITIONAL verdict.","tokens_in":9293,"tokens_out":9088,"duration_ms":83897,"concrete_test":"Recompute the stroboscopic Husimi distributions for N = 1000 and N = 2000 (same parameters, J=1, γ=0.1, A=3.4, T=2π) and extract the bagel onset U*(N) from the D(U) curves defined in Section III. If U*(N) does not approach the mean-field Neimark-Sacker value U_c≈0.11 as N increases (within numerical error), the quantum transition cannot be identified as the counterpart of the classical bifurcation. As a complementary check, fit the spectral gap 1-|μ_{2,3}| at U=0.1 versus N to a power law; a nonvanishing gap as N→∞ would indicate no sharp bifurcation in the classical limit.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the quantum system's transition at U≈0.1 is the quantum counterpart of the classical Neimark-Sacker bifurcation in the mean-field model (Eq. 7). This requires that the quantum signatures (bagel-shaped Husimi distribution, conjugate Floquet eigenvalues approaching the unit circle, rotation number ω≈0.58) are the finite-N precursors of the classical bifurcation. The paper, however, reports quantitative discrepancies that directly challenge this: at N=500 the bagel is already present at U=0.1 (Fig. 2a) while the mean-field still has a fixed point, and an additional quantum bifurcation appears for U∈[0.6,0.7] with no mean-field counterpart (Section III, Fig. 1). Furthermore, the spectral gap 1-|μ_{2,3}| is stated to decrease with N but no scaling law is given, and the D(U) curves in Fig. 6 are said to depend on N without a finite-size scaling analysis. Since the mean-field equations themselves are derived by dropping higher-order correlators and subleading-in-N dissipative terms (Section II), the validity of the classical reference for N=50..500 is not quantified. Without evidence that the quantum onset U*(N) converges to the mean-field bifurcation value U≈0.11 as N→∞, the observed bagel could be a finite-size crossover unrelated to the classical Neimark-Sacker mechanism.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a periodically modulated open bosonic dimer with N bosons subject to Lindblad dissipation, and reports a qualitative transition in the stroboscopic asymptotic state as the interaction strength U increases through U about 0.1: the Husimi distribution and the histograms of observables change from unimodal to bagel-shaped, a conjugate pair of Floquet eigenvalues approaches the unit circle, and trajectory-based rotation numbers become localized near omega about 0.58. The authors identify this transition as the quantum counterpart of the classical Neimark-Sacker bifurcation in the mean-field equations (Eq. 7), and they further show that the transition depends on the particle number N, which can itself serve as a bifurcation parameter.","tokens_in":9610,"tokens_out":4660,"duration_ms":46087,"significance":"If established, this result would extend the catalog of dissipative quantum bifurcations to the Neimark-Sacker (torus-birth) case, in a model that is experimentally relevant for cavity and circuit QED systems. The paper's strengths are its multi-probe numerical evidence: Husimi distributions, observable histograms, Floquet spectra, and rotation-number statistics all point to the same transition, and the quantum results are compared with the mean-field model without parameter fitting. The quantum-trajectory unraveling additionally resolves dynamics on the quantum attractor, which goes beyond earlier static Husimi pictures. The main weakness is that the quantitative link between the finite-N quantum transition and the mean-field bifurcation point is not established by a finite-size scaling analysis.","major_comments":[{"comment":"The central claim of a quantum Neimark-Sacker bifurcation requires the finite-N onset to converge to the mean-field bifurcation at U approximately 0.11 as N goes to infinity. The manuscript explicitly reports that at N=500 the bagel is already present at U=0.1 while the mean-field map still has a fixed point, and it states only that the spectral gap 1-|mu_{2,3}| decreases with N, without giving a quantitative scaling law. No finite-size scaling of U*(N), of the bagel diameter D(U), or of the spectral gap is provided, so the observed bagel could be a finite-size crossover rather than a precursor of the classical Neimark-Sacker bifurcation. Please provide estimates of U*(N) for the reported N values (for example, from a threshold in D(U) or in 1-|mu_{2,3}|), an extrapolation to N to infinity, and a comparison with U approximately 0.11, together with the scaling exponent of the spectral gap with N.","section":"Section III, Figs. 2 and 6 and the Floquet spectral paragraph"},{"comment":"The quantum bifurcation diagram in Fig. 1(b) shows an additional quantum transition within U in [0.6,0.7] that has no counterpart in the mean-field model. Because the mean-field model is the classical reference used to identify the bifurcation, this unexplained feature leaves the claimed correspondence incomplete. The authors should characterize this transition, for instance by testing whether its location and N-dependence indicate a separate finite-size instability or a second quantum Neimark-Sacker bifurcation, and by stating explicitly whether the feature persists as N increases.","section":"Section III, Fig. 1(b)"},{"comment":"The mean-field equations are obtained by truncating the cumulant hierarchy at the level of first-order expectation values and neglecting subleading-in-N dissipative terms, but the accuracy of this truncation for N=50..500 is not quantified. Since the mean-field bifurcation point U approximately 0.11 is used as the classical reference to which the quantum data are compared, the validity of this reference is load-bearing. A consistency check should be reported, for example the magnitude of the leading neglected second-order cumulants in the Heisenberg equations or a comparison with a next-order truncation, so that the reader can assess how much the mean-field bifurcation value may shift at finite N.","section":"Section II, Eq. (6)"}],"minor_comments":[{"comment":"The word 'qunntum torus' appears in the Introduction; it should be 'quantum torus'.","section":"Introduction"},{"comment":"The phrase 'The maximal element for each value of U is normalized to 1' is ambiguous; please specify the plotted quantity (for example, a density-matrix diagonal element, a histogram count, or a probability) and the binning or normalization procedure used.","section":"Fig. 1 caption"},{"comment":"Please specify how the average rotation number (solid line) is computed, including whether it is a time average or an ensemble average over trajectories, and provide error bars or a standard deviation to quantify the localization of the omega distribution.","section":"Fig. 4"},{"comment":"The statement that the largest eigenvalue 'is always unity' should be qualified: for a trace-preserving completely positive stroboscopic map with a unique asymptotic state there is an eigenvalue equal to 1, and numerical diagonalization of the finite-dimensional Floquet map will yield values close to but not exactly equal to 1.","section":"Section III, Floquet paragraph"},{"comment":"Reference [17] gives the publisher location as 'Brlin'; this should be 'Berlin'.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript presents a plausible and well-illustrated numerical observation, but the central quantum-classical correspondence is not quantitatively supported without a finite-size scaling analysis. The missing analysis is achievable with additional numerical work, so I recommend major revision rather than rejection. There are no concerns about citation practices or novelty disclosure in my reading."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is the first paper I know of that reports a quantum analogue of the Neimark-Sacker bifurcation, the torus-birth transition, in a driven open bosonic dimer. The previous work by this community covered pitchfork, saddle-node, and period-doubling; this adds the fourth type. The numerical evidence is coherent: as U crosses about 0.1 for N=500, the stroboscopic Husimi distribution changes from unimodal to bagel-shaped, the histograms of n and e do the same, a conjugate pair of Floquet eigenvalues approaches the unit circle with phase consistent with rotation number ω≈0.58, and quantum trajectories on the 'torus' show rational locking at 3/5 later. That is a solid package.\n\nThe paper is also honest about the discrepancies with the mean-field model: at U=0.1 the quantum bagel is already present while the classical fixed point is still stable, and an extra quantum bifurcation appears for U in [0.6,0.7] with no classical counterpart. Those are stated plainly, not hidden.\n\nWhere the paper is weaker: the connection to the classical Neimark-Sacker bifurcation is qualitative. The mean-field equations are derived by neglecting higher-order correlators and subleading-N dissipative terms, and the error of that truncation for N=50..500 is never quantified. The stress-test note says the finite-N onset U*(N) is not shown to converge to the classical U≈0.11; that is a fair criticism, and the paper does not provide the scaling analysis. The D(U) curves in Fig. 6 are said to depend on N without a finite-size collapse or even error bars. The spectral gap 1-|μ_{2,3}| is stated to decrease with N but no scaling law is given. These are all missing quantitative support, not contradictions of the main observation.\n\nThe central claim, though, survives the missing analysis. The paper does not assert that the quantum transition is exactly the classical bifurcation; it claims a quantum counterpart, and the multiple independent diagnostics (Husimi, observables, Floquet, rotation number) all point to the same qualitative change. The unexplained extra bifurcation is interesting and deserves follow-up, but it does not invalidate the main result.\n\nWho this is for: anyone working on driven open quantum systems, dissipative phase transitions, or quantum chaos analogues. It deserves a serious referee. I would send it to review and ask for the finite-size scaling and error quantification; without those the quantitative correspondence remains open, but the phenomenon is new and worth reporting.\n\nYours,","headline":"First numerical evidence for a quantum Neimark-Sacker bifurcation in a driven open dimer; the case is coherent but the link to the classical mean-field limit is asserted rather than demonstrated.","tokens_in":10096,"tokens_out":2448,"would_cite":true,"duration_ms":22425,"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":"The paper establishes a quantum analog of the Neimark-Sacker bifurcation in a periodically driven open bosonic dimer.","keywords":["quantum Neimark-Sacker bifurcation","open quantum dimer","Lindblad master equation","Floquet map","Husimi distribution","quantum trajectories","rotation number","driven dissipative bosons"],"falsifier":"Exact diagonalization of the Floquet map at $N=500$, $U=0.1125$, should show a bagel-shaped Husimi distribution and a conjugate eigenvalue pair with $|\\mu_{2,3}|$ close to 1 and phase $\\theta_0\\approx 2\\pi\\times 0.58$; if the distribution remains unimodal or no such pair exists at that parameter set, the claimed quantum Neimark-Sacker bifurcation is not present there.","tokens_in":9124,"feed_emoji":"🥯","tokens_out":8278,"duration_ms":67603,"temperature":0.7,"pith_summary":"The paper reports that a periodically modulated open quantum dimer undergoes a qualitative transition that mirrors the classical Neimark-Sacker bifurcation: as the boson interaction strength $U$ crosses about $0.1$ (for $N\\ge 25$ particles), the stroboscopic asymptotic state changes from a point-like unimodal distribution to a bagel-shaped one, the signature of a torus in the classical Poincaré section. The same bagel appears in histograms of directly measurable observables sampled by quantum trajectories, not only in the Husimi quasiprobability. A conjugate pair of Floquet eigenvalues approaches the unit circle at a phase consistent with a rotation number $\\omega\\approx 0.58$, and at larger $U$ the rotation number locks to $3/5$, giving a five-periodic structure. If correct, this extends the catalogue of quantum bifurcations from pitchfork, saddle-node, and period doubling to the birth of a torus, and it makes the particle number itself a usable bifurcation parameter.","feed_headline":"Bagel-shaped states signal a quantum Neimark-Sacker bifurcation","feed_subtitle":"In a driven open dimer, a stable point blossoms into a torus-like bagel at interaction strength U≈0.1.","key_machinery":"The central objects are the stroboscopic Floquet map $P_F=\\mathcal{T}\\exp[\\int_0^T \\mathcal{L}\\,dt]$, the one-period evolution operator of the Lindblad master equation, and the mean-field Bloch-sphere equations obtained by replacing bosonic operators with expectation values. The Floquet eigenvalues play the role of classical multipliers, with a conjugate pair approaching the unit circle at the bifurcation; the Husimi distribution built from SU(2) coherent states visualizes the attractor shape; and the Monte-Carlo wave-function unraveling supplies individual trajectories whose polar angle in the $(n,e)$ plane defines the rotation number $\\omega_m$.","core_discovery":"For the open dimer with hopping $J=1$, dissipation $\\gamma=0.1$, modulation amplitude $A=3.4$, and period $T=2\\pi$, the paper shows that the quantum stroboscopic state reproduces the mean-field Neimark-Sacker scenario: the fixed point of the Poincaré map loses stability and an invariant curve is born, while the quantum Husimi distribution and the stroboscopic distributions of particle number and energy become bagel-shaped. The Floquet map's subleading eigenvalues $\\mu_{2,3}$ approach the unit circle with phase $\\theta_0\\approx 2\\pi\\omega$, where the rotation number $\\omega\\approx 0.58$ is measured from the winding of individual quantum trajectories on the attractor. Increasing $U$ leads to frequency locking at $\\omega=3/5$ and a period-5 structure on the torus, and increasing $N$ at fixed $U$ converts a unimodal distribution into a bagel, so the bifurcation is controlled both by interaction strength and by system size.","pith_inferences":["If the finite-$N$ shift of the apparent bifurcation point scales as a power of $1/N$, measuring that shift in the dimer would quantify how much the mean-field reference overestimates or underestimates the true quantum threshold; the paper reports the shift but does not extract this scaling.","The same stroboscopic rotation-number diagnostic could be applied to other driven-dissipative bosonic systems, such as an open Dicke model, to search for torus bifurcations from time-series data alone.","The observed relation between the Floquet eigenvalue phase and the trajectory rotation number suggests that rational plateaus of the rotation number could be predicted directly from spectral properties without simulating trajectories.","Comparing the spread of quantum trajectories on the bagel with the width predicted from the spectral gap would allow a test of whether quantum fluctuations, rather than finite sampling, set the thickness of the quantum torus."],"forward_implications":["At $N=500$ the quantum bifurcation diagram reproduces the classical sequence of torus birth, period-6 cycle, chaos, and crisis, with the quantum transition occurring slightly earlier ($U\\approx 0.1$).","The rotation-number distribution is well localized near $\\omega\\approx 0.58$ just after the bifurcation, clearly distinct from the period-doubling value $1/2$, and locks to the rational $3/5$ at $U\\approx 0.15$.","The spectral gap $1-|\\mu_{2,3}|$ decreases as $N$ grows, so the relaxation time to the asymptotic state, estimated as $t\\sim(1-|\\mu_{2,3}|)^{-1}$, can become very large near the quantum bifurcation.","The number of bosons acts as a bifurcation parameter: at $U=0.1125$ the Husimi distribution is unimodal for $N=50$ but develops a bagel by $N=500$.","The diameter of the quantum bagel grows with $U$ in a way that depends on $N$, so the finite-size scaling near the bifurcation is not universal."],"supporting_citations":[{"why":"It defines the classical Neimark-Sacker bifurcation as the birth of an invariant curve from a fixed point and the corresponding crossing of the unit circle by a conjugate pair of multipliers.","marker":"[2]"},{"why":"It established the quantum period-doubling bifurcation in the same driven dimer and supplied the Floquet-map spectral method used here.","marker":"[24]"},{"why":"It showed that quantum bifurcations can be detected directly in the structure of the asymptotic density matrix, the approach the paper extends to observables.","marker":"[25]"},{"why":"It provided the treatment of two-time correlations in the asymptotic regime that the paper adapts to relate the Floquet phase to the rotation number.","marker":"[30]"},{"why":"It introduced the Monte-Carlo wave-function method that yields individual quantum trajectories on the quantum torus.","marker":"[35]"},{"why":"It developed the stochastic unraveling of the Lindblad equation used to sample the stroboscopic observables.","marker":"[36]"},{"why":"It supplies the high-performance realization of quantum trajectories that makes large-$N$ simulations feasible.","marker":"[50]"}],"fun_headline_variants":["Quantum torus emerges from dimer instability","Bagel-shape reveals quantum Neimark-Sacker","Quantum dimers birth a torus at critical U","From point to bagel: quantum bifurcation","Quantum Neimark-Sacker seen in driven dimer"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole analogy rests on the assumption that the mean-field equations obtained by replacing quantum operators with expectation values are accurate for $N$ between 50 and 500, and the paper does not quantify the error of that semiclassical truncation.","fun_headline_variants_meta":{"raw":{"variants":["Quantum torus emerges from dimer instability","Bagel-shape reveals quantum Neimark-Sacker","Quantum dimers birth a torus at critical U","From point to bagel: quantum bifurcation","Quantum Neimark-Sacker seen in driven dimer"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000173,"raw_usage":{"total_tokens":1282,"prompt_tokens":949,"completion_tokens":333,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":565,"completion_tokens_details":{"reasoning_tokens":274}},"tokens_in":565,"tokens_out":333,"duration_ms":3516,"temperature":1.0,"reasoning_tokens":274,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:48:02.998274+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Exact diagonalization of the Floquet map at $N=500$, $U=0.1125$, should show a bagel-shaped Husimi distribution and a conjugate eigenvalue pair with $|\\mu_{2,3}|$ close to 1 and phase $\\theta_0\\approx 2\\pi\\times 0.58$; if the distribution remains unimodal or no such pair exists at that parameter set, the claimed quantum Neimark-Sacker bifurcation is not present there.","supporting_citations":[{"cited_title":"Kuznetsov, Elements of Applied Bifurcation The- ory, Springer, 3rd edition (2004)","cited_arxiv_id":null,"evidence_quote":"It defines the classical Neimark-Sacker bifurcation as the birth of an invariant curve from a fixed point and the corresponding crossing of the unit circle by a conjugate pair of multipliers."},{"cited_title":"Hartmann, D","cited_arxiv_id":null,"evidence_quote":"It established the quantum period-doubling bifurcation in the same driven dimer and supplied the Floquet-map spectral method used here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It showed that quantum bifurcations can be detected directly in the structure of the asymptotic density matrix, the approach the paper extends to observables."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It provided the treatment of two-time correlations in the asymptotic regime that the paper adapts to relate the Floquet phase to the rotation number."},{"cited_title":"Dum, A.S","cited_arxiv_id":null,"evidence_quote":"It introduced the Monte-Carlo wave-function method that yields individual quantum trajectories on the quantum torus."},{"cited_title":"Volokitin, A","cited_arxiv_id":null,"evidence_quote":"It supplies the high-performance realization of quantum trajectories that makes large-$N$ simulations feasible."}],"review_version":1}