{"id":"526f3882-cb63-4810-816f-7151eb0cb3c5","arxiv_id":"2411.13936","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Nonlinear modes in a quasiperiodic lattice form via symmetry-breaking pitchfork or saddle-node bifurcations below the mobility edge and via gap-soliton-like localization above it, while asymmetric potentials produce cascades of saddle-node bifurcations.","lead":"Researchers mapped how nonlinear localized modes are born in a one-dimensional Bose-Einstein condensate trapped in a quasiperiodic two-frequency potential. The work classifies bifurcation patterns below and above the mobility edge and may help guide experiments in ultracold atoms and photonic lattices.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Approximant-path robustness is the weakest point: Table I shows band-edge classifications change across rational approximants, and the above-mobility-edge evidence for a pitchfork bifurcation and sqrt-law N(mu) relies on near-degenerate pairs at p/q=89/55 with spacings at numerical precision.","rationale":"The reader's weakest assumption is the approximant-path faithfulness, and the paper's own Table I and Sec. III gap-E discussion expose exactly that weakness. The concern is load-bearing because the central claims are universal statements about quasiperiodic systems, yet every nonlinear family is computed for one approximant. The dichotomy in N(mu) may survive since localized modes below the mobility edge have O(1) IPR while extended modes above it have O(L) IPR, so the slope difference is structural; however, the specific pitchfork/saddle-node taxonomy near a given gap is tied to the band-edge type, which Table I shows is not stable across approximants. Because the manuscript already conditions on this issue and the remedy is a well-defined convergence study, the existing CONDITIONAL verdict is appropriate and no adjustment is needed.","tokens_in":17210,"tokens_out":7423,"duration_ms":77983,"concrete_test":"Compute the same families for p/q=55/34, 144/89, and 233/144 using the paper's stated 4500 Fourier harmonics and v1=v2=0.8. Check (i) whether the gap-E family that bifurcates from the right band edge still has a symmetry-breaking pitchfork when the edge is a single extended state at 144/89 rather than a near-degenerate pair; (ii) whether the N(mu) curves above the mobility edge converge to a common mu^{1/2} asymptote as q grows, while the below-edge slopes dN/dmu remain O(1) and independent of q; (iii) whether the gap-D pitchfork reappears at 233/144 or is replaced by boundary-mode behaviour. If the bifurcation type or the scaling changes with q, the claimed patterns are rational-approximant artifacts rather than properties of the quasiperiodic system.","verdict_should_be":"UNCHANGED","load_bearing_attack":"All nonlinear bifurcation diagrams are computed for a single rational approximant, p/q=89/55, although the paper's claims are about the infinite quasiperiodic potential. The approximant path is not shown to be convergent for the specific objects studied. In Sec. III, Table I lists the right band-edge structures for gaps A-D and shows that the same gap changes type between 's-a/s pair', 'boundary', and 'new s-a/s pair' as q increases; in particular gap D, the flagship below-mobility-edge pitchfork example in Fig. 4, becomes a 'boundary' edge at p/q=144/89. Above the mobility edge the instability of the approximant picture is acknowledged in the text: for gap E the right edge is a pair of virtually coinciding extended eigenfunctions for every approximant except 144/89, with spacings ~10^-8 and ~10^-11, comparable to numerical error. The entire gap-E analysis (Fig. 6), including the symmetry-breaking pitchfork and the square-root-like N(mu) dependence, is computed only at 89/55. Therefore the central claim that the mobility edge dichotomizes the bifurcation structure, and the specific taxonomy of pitchfork vs saddle-node, are not yet established in the infinite-system limit; they may depend on which finite approximant is selected.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies nonlinear steady states of the one-dimensional Gross-Pitaevskii equation with a bichromatic quasiperiodic potential, using rational approximants of the golden-ratio frequency ratio with periodic boundary conditions. The authors classify the right band edges of selected spectral gaps (A–E) and compute bifurcation diagrams of nonlinear modes near those gaps for the approximant p/q = 89/55 with lattice amplitudes v1 = v2 = 0.8. They report symmetry-breaking pitchfork bifurcations for near-degenerate symmetric/antisymmetric pairs below the mobility edge, saddle-node bifurcations due to hybridization of modes with different spatial profiles, a different (square-root-like) N(mu) behavior above the mobility edge, and a cascade of saddle-node bifurcations for a nonzero phase shift. The central claim is that the mobility edge dichotomizes the formation of nonlinear modes: linear N(mu) dependence below, square-root-like above, with the bifurcation taxonomy mimicking symmetric and asymmetric double-well potentials.","tokens_in":17464,"tokens_out":2721,"duration_ms":28761,"significance":"If established, the claimed connection between the fractal linear spectrum, the mobility edge, and the bifurcation structure of nonlinear localized modes would be a useful organizing principle for quasiperiodic BECs and photonic lattices, with clear experimental relevance. The numerical methodology is standard and adequately described: Fourier collocation with 4500 harmonics, Newton's method, and linear-stability analysis via Eq. (8). The perturbation formulas (10)–(11) are standard and Eq. (12) is explicitly checked against numerics. The numerical results are obtained by direct solution of the nonlinear PDE, so there is no circularity in the bifurcation detection. The main weakness is the reliance on a single rational approximant for all nonlinear bifurcation diagrams, while the band-edge classification itself visibly changes across approximants in Table I. This makes the central claims about universal patterns and the role of the mobility edge not yet fully established in the infinite-quasiperiodic limit.","major_comments":[{"comment":"The manuscript's central claim of universal patterns is undermined by the approximant dependence of the band-edge classification. Table I shows that the same gap changes its right-edge structure as p/q increases: gap A switches between 's-a/s pair', 'boundary', and 'new s-a/s pair'; gap B switches between 'boundary' and 's-a/s pair'; and gap D, the flagship below-mobility-edge pitchfork example in Fig. 4, becomes a 'boundary' edge at p/q = 144/89. Since all nonlinear bifurcation diagrams in Sec. IV are computed only for p/q = 89/55, the paper does not demonstrate that the pitchfork-vs-saddle-node taxonomy is the infinite-system behavior rather than an artifact of that particular approximant. The authors should either compute the relevant bifurcation diagrams for at least one more approximant (e.g., 55/34, 144/89, or 233/144) for gaps C, D, and E, or provide a convergence argument showing that the bifurcation type stabilizes for sufficiently large q.","section":"Sec. III, Table I"},{"comment":"The above-mobility-edge analysis for gap E rests on numerically fragile data. The band edge at p/q = 89/55 is formed by a pair of eigenvalues with splitting ~10^-11, which the authors themselves note is comparable to numerical error, and for p/q = 144/89 the edge is instead a single extended eigenfunction well separated from the coincident pair. The symmetry-breaking bifurcation and the square-root-like N(mu) law shown in Fig. 6 are computed only at 89/55. This does not establish the claimed dichotomy between below- and above-mobility-edge bifurcation scenarios in the quasiperiodic limit. A direct convergence test (e.g., increasing the number of Fourier harmonics for the near-degenerate pair, or repeating the gap-E bifurcation analysis at another approximant) is needed before the central claim can be accepted.","section":"Sec. IV D, Fig. 6"},{"comment":"The asserted square-root law for N(mu) above the mobility edge is not actually derived or quantitatively verified. Equation (11) gives a linear dependence N ≈ (mu - mu_n)/chi_n for a bifurcation from an isolated linear mode with fixed IPR. For the extended modes above the mobility edge, chi is not constant along the family, so a different law is expected, but the paper only states that the dependence 'resembles the root-law behavior' and attributes it to similarity with gap solitons in periodic media (Sec. IV D). To make the 'linear below, square-root above' dichotomy load-bearing, the authors should provide a quantitative check (e.g., a log-log fit of N(mu) over the gap, or a reduced-model derivation) and show that it holds for at least one additional approximant.","section":"Sec. IV A and Sec. IV D"}],"minor_comments":[{"comment":"Typographical errors: 'eigenvalue' appears as 'and eigenvalue' and 'the the remaining' appears in the enumeration of spectral properties. These should be corrected.","section":"Sec. II, Eq. (8)"},{"comment":"The phrase 'nonlinear modes theta bifurcate near those bands' in the opening sentence of Sec. IV A should read 'that bifurcate'.","section":"Sec. IV A"},{"comment":"The caption of Fig. 5 is dense and the labels SB1, SB2, SN1, SN2 are not all explicitly explained in the caption; a short list of which line corresponds to which family would improve readability.","section":"Fig. 5"},{"comment":"The caption of Table I states 'the structure of the right band edges adjacent to gaps A–D from Fig. 1' but does not mention that gap E is discussed in the text; adding this would help the reader locate the full classification.","section":"Sec. III, Table I"},{"comment":"The asymmetric case is illustrated for a single arbitrarily chosen phase shift theta ≈ 1.7π in gap B. While the cascade of saddle-node bifurcations is plausible, a brief statement about whether θ = 1.7π is generic (e.g., not accidentally restoring symmetry) would strengthen the presentation.","section":"Sec. IV E"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the scope of the journal and presents a worthwhile numerical study, but the central claims about universal, approximant-independent bifurcation patterns are not yet supported because all nonlinear computations use p/q = 89/55 while Table I shows that the band-edge classification itself varies with q. The authors are clearly aware of the approximant path and even label some structures as 'new s-a/s pair', so the missing piece is a systematic convergence check for the bifurcation diagrams. This is fixable within the manuscript's scope, hence major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: a clear, competent numerical bifurcation study that will be useful, but the claims about the above-mobility-edge regime rest on one rational approximant and need a convergence check before I treat them as established.\n\nWhat's actually new: the paper extends the authors' prior program (Ref. [50], symmetry breaking at the left band edge) to right band edges and to modes above the mobility edge, in a symmetric bichromatic potential with v1=v2=0.8. It classifies right-edge structures in terms of symmetric-antisymmetric pairs, boundary-localized modes, and a single robust symmetric mode, and shows numerically that below the mobility edge N(mu) is linear while above it looks square-root-like. It also shows, for one nonzero phase shift, a cascade of saddle-node bifurcations. The methods are standard (Fourier collocation with 4500 harmonics, Newton, linear stability), the perturbation formulas (10)-(11) and the two-mode estimate (12) are checked against the numerics, and the paper is honest about the gap-E caveats. That is genuinely useful, especially for experimental groups planning quasiperiodic BEC or photonic experiments.\n\nSoft spots. The approximant path is the load-bearing weakness. Table I shows the right-edge structure of gaps A, B, D changing type between 's-a/s pair' and 'boundary' as q increases; gap D, the flagship below-edge pitchfork example, becomes a boundary edge at p/q=144/89. Above the mobility edge, the gap E 'pair of virtually coinciding extended eigenfunctions' has spacings ~1e-8 and ~1e-11 at 55/34 and 89/55, comparable to numerical precision, and at 144/89 the right edge is a single extended eigenfunction. All nonlinear diagrams in gap E are computed only at 89/55. So the most novel claim - that the mobility edge dichotomizes the bifurcation structure and the scaling of N(mu) - is not yet robust in the infinite-system limit; it may be sensitive to which approximant is selected. The generic asymmetric case is one example at one theta. No code, data, or explicit tolerances are provided, which makes independent verification harder.\n\nDoes the central argument hold? The below-edge picture - linear N(mu) and pitchfork/saddle-node from localized double-well pairs - is well supported. The above-edge picture is plausible but not yet established. The paper overclaims slightly with 'general patterns' given the limited approximant data.\n\nWho for: researchers in quasiperiodic lattices, nonlinear waves, BECs. A serious referee could add value by asking for a convergence study and a second approximant for the nonlinear runs. I'd accept it for peer review, and I'd cite it for the classification and below-edge results with a caveat.","headline":"Useful numerical bifurcation study that needs a convergence check before the above-mobility-edge claims (pitchfork, square-root N(mu)) are trusted.","tokens_in":18027,"tokens_out":3873,"would_cite":true,"duration_ms":34856,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q55","37G10","37G15"],"pacs":[],"model":"deepseek-v4-flash","headline":"In a one-dimensional Bose-Einstein condensate with a bichromatic quasiperiodic potential, nonlinear localized modes are born through pitchfork or saddle-node bifurcations dictated by the band-edge structure, with the mobility edge…","keywords":["Bose-Einstein condensate","quasiperiodic lattice","mobility edge","nonlinear modes","pitchfork bifurcation","saddle-node bifurcation","Gross-Pitaevskii equation","approximant path"],"falsifier":"Repeat the bifurcation analysis at the next Fibonacci approximants ($p/q=144/89$, $233/144$, $377/233$) or with a different irrational frequency ratio and check whether the predicted patterns persist: a symmetric-antisymmetric edge that changes to a boundary edge and turns a pitchfork into a different bifurcation, or an above-edge $N(\\mu)$ branch that deviates from the square-root law, would refute the claimed universality.","tokens_in":16939,"feed_emoji":"⚛️","tokens_out":10136,"duration_ms":89855,"temperature":0.7,"pith_summary":"The paper studies how steady nonlinear modes are born from the linear spectrum in a one-dimensional Bose-Einstein condensate with repulsive interactions trapped in a bichromatic quasiperiodic potential. Working with rational approximants of the golden-ratio potential, it finds that the mobility edge splits the formation mechanism in two: below it, the linear mode is already localized and the atom number grows linearly with chemical potential, $N(\\mu)\\approx (\\mu-\\tilde{\\mu}_n)/\\tilde{\\chi}_n$; above it, an extended linear mode gradually localizes and $N(\\mu)$ follows a square-root-like law. In a symmetric potential, nonlinear modes are born through both pitchfork and saddle-node bifurcations that mimic symmetric and asymmetric double wells, while a nonzero phase shift between the lattices produces a cascade of saddle-node bifurcations. If these patterns are robust, they connect the fractal linear spectrum and the mobility edge directly to the multiplicity of nonlinear localized states and give experimental predictions for atom-number scaling.","feed_headline":"Mobility edge separates two ways nonlinear modes are born","feed_subtitle":"Below mobility edge: linear growth; above: square-root. Phase shifts trigger saddle-node cascades.","key_machinery":"The argument is carried by the approximant path: replace the irrational frequency ratio $\\phi$ by a Fibonacci fraction $p/q$, restrict the problem to the periodic cell $I_q=[-\\pi q/2,\\pi q/2)$ with periodic boundary conditions, and solve the Gross-Pitaevskii equation $\\mu\\psi=H\\psi+g\\psi^3$ by Fourier collocation and Newton's method. The organizing device is a classification of right band edges as symmetric-antisymmetric pairs, boundary-localized modes, or isolated symmetric modes. The bifurcation analysis rests on the perturbation formula $\\psi(x)\\approx((\\mu-\\tilde{\\mu}_n)/\\tilde{\\chi}_n)^{1/2}\\tilde{\\psi}_n(x)$, which yields $N(\\mu)\\approx(\\mu-\\tilde{\\mu}_n)/\\tilde{\\chi}_n$, where $\\tilde{\\chi}_n$ is the inverse participation ratio; this formula is what makes below-the-edge branches linear and above-the-edge branches steep and root-like.","core_discovery":"For the symmetric bichromatic potential $V(x)=v_1\\cos(2x)+v_2\\cos(2\\phi x)$ with $\\phi$ the golden ratio, the right edge of a linear band can be a symmetric-antisymmetric pair, a single strongly localized symmetric mode, or a mode pinned to the domain boundary. The paper shows that each edge structure dictates a distinct nonlinear scenario: symmetric-antisymmetric pairs trigger pitchfork symmetry-breaking bifurcations whose particle number is set by the ratio of the eigenvalue splitting to the inverse participation ratio; isolated modes hybridize with nearby states to form in-phase and out-of-phase bound states through saddle-node bifurcations; and above the mobility edge, extended linear states develop into gap-soliton-like families with a root-like $N(\\mu)$ curve. With a nonzero inter-lattice phase shift, the potential is generically asymmetric, the special eigenvalue pairs disappear, and every linear mode produces a nonlinear family that meets the next eigenvalue, yielding a chain of saddle-node bifurcations before localized modes reach the gap. The central claim is that the mobility edge and the Fibonacci-labelled gap structure of the linear spectrum organize which bifurcation pattern applies.","pith_inferences":["Beyond the paper's finite-approximant numerics, the setup suggests that in the true quasiperiodic limit the asymmetric case may exhibit an infinite accumulation of saddle-node bifurcations near a band edge, a property the paper does not directly establish.","The Fibonacci-labelled gaps are special to the golden ratio; for other irrational ratios the universal dichotomy (pitchfork for symmetric, saddle-node cascade for asymmetric) may persist even though the ordering of modes at band edges changes, so testing another ratio would separate the general mechanism from golden-ratio-specific structure.","In the optical realization, the square-root $N(\\mu)$ law above the mobility edge translates into a measurable power-versus-propagation-constant curve for soliton formation in quasiperiodic photonic lattices, and the predicted symmetry-breaking particle number $N_{\\mathrm{SB}}\\approx(\\tilde{\\mu}_{55}-\\tilde{\\mu}_{54})/\\tilde{\\chi}_{55}$ could be checked directly in a cold-atom experiment."],"forward_implications":["Below the mobility edge, atom number grows linearly with chemical potential, $N(\\mu)\\approx (\\mu-\\tilde{\\mu}_n)/\\tilde{\\chi}_n$, so measuring $N(\\mu)$ near a band edge gives the inverse participation ratio of the underlying linear mode.","Above the mobility edge, an initially extended linear mode gradually localizes as $\\mu$ moves into the gap, producing a square-root-like $N(\\mu)$ branch characteristic of gap solitons in periodic media.","In a symmetric bichromatic potential, band edges formed by symmetric-antisymmetric pairs produce pitchfork symmetry breaking at a particle number set by the eigenvalue splitting divided by the inverse participation ratio; isolated or hybridizing edge modes produce saddle-node bifurcations and in-phase/out-of-phase bound states.","A nonzero phase shift between the two lattices removes the symmetry, converting pitchforks into a cascade of saddle-node bifurcations in which each linear mode contributes a localized branch that reaches the gap.","Because the Gross-Pitaevskii equation is mathematically equivalent to the nonlinear Schrödinger equation for paraxial light, the same bifurcation patterns should appear as optical mode formation in quasiperiodic photonic lattices."],"supporting_citations":[{"why":"supplies the approximant path that replaces the irrational frequency ratio by rational Fibonacci fractions, the method's foundation.","marker":"[46]"},{"why":"establishes mobility edges in one-dimensional bichromatic incommensurate potentials, providing the linear-spectrum dichotomy the paper builds on.","marker":"[49]"},{"why":"gives the two-mode approximation and symmetry-breaking bifurcation estimate used to locate pitchfork bifurcations and to compare with prior quasiperiodic-lattice results.","marker":"[50]"},{"why":"supplies the double-well bifurcation paradigm (pitchfork in symmetric wells, saddle-node in asymmetric wells) that the paper maps onto quasiperiodic band edges.","marker":"[67]"},{"why":"recent few-mode analysis predicting saddle-node bifurcations in quasiperiodic potentials, which the paper extends to full cascades.","marker":"[62]"},{"why":"provides the numerical methods (Newton's method, linear-stability eigenvalue problem) and the periodic gap-soliton behavior used to interpret modes above the mobility edge.","marker":"[51]"},{"why":"supplies the perturbation theory for bifurcations from an isolated linear eigenfunction, giving the $N(\\mu)$ approximation used throughout.","marker":"[63]"}],"fun_headline_variants":["Mobility edge dictates nonlinear mode birth patterns","Phase shift triggers cascade of saddle-node bifurcations","Symmetry breaking vs saddle-node: edge matters","Mobility edge splits nonlinear mode formation scenarios","Two paths to nonlinear modes, separated by mobility edge"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusions are drawn from finite rational approximants of the golden-ratio potential with periodic boundary conditions, and the paper assumes this approximant path faithfully represents the true infinite quasiperiodic system even though Table I shows the band-edge structure switching between approximants.","fun_headline_variants_meta":{"raw":{"variants":["Mobility edge dictates nonlinear mode birth patterns","Phase shift triggers cascade of saddle-node bifurcations","Symmetry breaking vs saddle-node: edge matters","Mobility edge splits nonlinear mode formation scenarios","Two paths to nonlinear modes, separated by mobility edge"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000166,"raw_usage":{"total_tokens":1269,"prompt_tokens":973,"completion_tokens":296,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":589,"completion_tokens_details":{"reasoning_tokens":225}},"tokens_in":589,"tokens_out":296,"duration_ms":3555,"temperature":1.0,"reasoning_tokens":225,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:42:40.473972+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the bifurcation analysis at the next Fibonacci approximants ($p/q=144/89$, $233/144$, $377/233$) or with a different irrational frequency ratio and check whether the predicted patterns persist: a symmetric-antisymmetric edge that changes to a boundary edge and turns a pitchfork into a different bifurcation, or an above-edge $N(\\mu)$ branch that deviates from the square-root law, would refute the claimed universality.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the approximant path that replaces the irrational frequency ratio by rational Fibonacci fractions, the method's foundation."},{"cited_title":"Huang, Z","cited_arxiv_id":null,"evidence_quote":"establishes mobility edges in one-dimensional bichromatic incommensurate potentials, providing the linear-spectrum dichotomy the paper builds on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"gives the two-mode approximation and symmetry-breaking bifurcation estimate used to locate pitchfork bifurcations and to compare with prior quasiperiodic-lattice results."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the double-well bifurcation paradigm (pitchfork in symmetric wells, saddle-node in asymmetric wells) that the paper maps onto quasiperiodic band edges."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"recent few-mode analysis predicting saddle-node bifurcations in quasiperiodic potentials, which the paper extends to full cascades."},{"cited_title":"Takahashi, H","cited_arxiv_id":null,"evidence_quote":"provides the numerical methods (Newton's method, linear-stability eigenvalue problem) and the periodic gap-soliton behavior used to interpret modes above the mobility edge."},{"cited_title":"Larcher, F","cited_arxiv_id":null,"evidence_quote":"supplies the perturbation theory for bifurcations from an isolated linear eigenfunction, giving the $N(\\mu)$ approximation used throughout."}],"review_version":1}