{"id":"9d53c169-c198-4892-b59f-621f1903cd9b","arxiv_id":"2505.02776","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A Hermitian parametric pairing model with local damping exhibits non-Hermitian topology and skin modes in the continuum, with a bulk-boundary correspondence established via winding numbers and non-Bloch theory.","lead":"This paper shows that a Hermitian, nonlocal parametric pairing process, combined with uniform local loss, can produce non-Hmitian skin modes in a continuum bosonic system. The result offers a simpler route to non-Hermitian topological phases than engineering complex baths, which may ease experimental realization in photonic, atomic, or acoustic platforms.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The spectral winding number in Eq. (4) is not well defined for the compactified continuum: det(D_k−E) ~ −k^4 as |k|→∞, so the arctan compactification places a pole on the integration contour and the claimed W=±1 regions are not established.","rationale":"The reader's weakest assumption identifies the same load-bearing step: the compactification of continuum momentum used to define the spectral winding number. My analysis sharpens the concern: it is not merely that convergence or boundary terms are unproven; for the given dynamical matrix the determinant grows as −k^4, so the compactified loop passes through a fourth-order pole of det(D_k−E). Consequently the winding integral of Eq. (4) is not a well-defined loop integral in C* without a regularization that the paper does not provide. A direct cutoff evaluation suggests the naive integral would vanish in the continuum limit, making the claimed W=±1 regions dependent on an unspecified prescription. This is significant because the abstract and the bulk-boundary discussion present the spectral winding number as one of the two tools establishing a robust bulk-boundary correspondence. However, the paper contains independent and much more secure evidence: the open-boundary numerical diagonalization shows exponentially localized modes, and the non-Bloch calculation, which does not rely on the compactification, reproduces the open-boundary spectrum and the localization lengths analytically. Those checks are reasonably convincing and support the central physical claim of a non-Hermitian skin effect from parametric pairing with local damping. The compactification gap therefore warrants a conditional verdict rather than rejection: the topological-index interpretation is not established as stated, but the main phenomenon and the non-Bloch correspondence are supported. The verdict should remain CONDITIONAL, as the reader already concluded, so no change is needed.","tokens_in":12238,"tokens_out":14536,"duration_ms":184740,"concrete_test":"Take the parameters of Fig. 3(a) (k0=0.3, g=γ=1) and a complex energy E in the claimed W=+1 region, e.g., E≈0.23−1.83i from Appendix B. Compute W_Λ(E)=(1/2πi)∫_{−Λ}^{Λ} dk d/dk log det(D_k−E) for increasing Λ, with a fixed branch of the logarithm and with k=tan(k̃/2). If W_Λ does not converge to an integer independent of Λ and of the branch choice, or if it converges to 0, then Eq. (4) as written does not produce the stated W=±1. Then repeat with the natural regularization det(D_k−E)/(1+k^2)^2 and state whether the ±1 regions survive; this would identify the missing prescription needed to make the compactified topological index well defined.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The compactification used for Eq. (4) is the load-bearing step for the topological classification. For the dynamical matrix (2) with damping, det(D_k−E) is a quartic polynomial in k with leading behavior −k^4 as |k|→∞. Under k̃=arctan(k), or equivalently k=tan(k̃/2), the matrix function is therefore not a smooth finite periodic function on the Brillouin circle: it has a pole of order four at the identified point k=±∞, and log det diverges there. The integral in Eq. (4) is then not defined as an ordinary winding integral, and the statement that compactification 'guarantees a well quantized topological index' is unsupported. In fact, evaluating the same winding expression over k∈[−Λ,Λ] gives log det at ±Λ, whose leading −Λ^4 terms cancel, so the direct continuum limit yields 0 rather than the claimed ±1; any nonzero result requires an additional regularization or normalization, e.g., dividing det by (1+k^2)^2, that is not specified. The claimed regions W(E)=±1 in Fig. 3(a) are therefore cutoff- and branch-dependent until a prescription is provided. This does not call into question the numerically observed skin modes or the non-Bloch spectrum, which agree independently, but it does remove the bulk spectral-winding argument as evidence for the bulk-boundary correspondence.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a one-dimensional continuum bosonic model with momentum-shifted parametric pairing, Eq. (1), and shows that the Hermitian pairing produces level attraction between shifted particle and hole branches. With uniform local damping added, the bulk spectrum contains a tilted diabolical line, and finite systems display non-Hermitian skin modes. The authors support this with exact numerical diagonalization of a discretized real-space model and with non-Bloch theory, obtaining analytical localization lengths that match numerics. They also attempt to classify the skin effect through a spectral winding number computed after compactifying the continuum momentum via k~=arctan(k), and claim W=±1 regions for energies inside the diabolical-line area.","tokens_in":12547,"tokens_out":16851,"duration_ms":210756,"significance":"The proposed mechanism—Hermitian parametric pairing generating non-Hermitian skin modes without complex bath engineering—is conceptually appealing and potentially impactful for bosonic platforms. The non-Bloch analysis is a genuine strength: no parameter is fitted, the localization lengths (|Im kbar|≈0.345 and 0.339 for g=γ=1, k0=0.3) are analytical predictions that agree with exact diagonalization, and the open-boundary spectrum from Eq. (5) matches numerics. However, the spectral-winding classification, one of the two advertised routes to the bulk-boundary correspondence, is not currently well defined. The significance of the paper depends on whether that index can be given a rigorous regularization; the non-Bloch route alone already establishes the main boundary physics.","major_comments":[{"comment":"The compactified spectral winding in Eq. (4) is not a well-defined integral for this model. With the damped dynamical matrix D~k = σz Hk − iγ1, one has det(D~k − E) = −(k^2−k0^2)^2 + 4(E+iγ)k0 k + (E+iγ)^2 + g^2, which behaves as −k^4 + O(k^2) as |k|→∞. Under the substitution k=tan(k~/2) used in the compactification, the integrand d/dk~ log det(D~k − E) has a fourth-order pole at the identified point k~=π/2. Consequently the integral in Eq. (4) is not defined as an ordinary winding integral, and the statement that this compactification 'guarantees a well quantized topological index' is unsupported. A symmetric cutoff k∈[−Λ,Λ] followed by Λ→∞ gives (1/2πi)[log det(D~Λ−E)−log det(D~−Λ−E)] → 0, i.e., the direct continuum limit is 0 rather than the claimed ±1, while nonsymmetric cutoffs give cutoff-dependent results. The W(E)=±1 regions in Fig. 3(a) therefore require an explicit regularization (for example, dividing det by (1+k^2)^2) that is not specified. Please supply such a prescription or remove the spectral-winding claim and rely on the non-Bloch argument, which is independent and appears sound.","section":"Eq. (4) and the paragraph introducing k~=arctan(k)"},{"comment":"The bulk-anomaly argument is presented as a consequence of Eq. (3), but the current Jk used to color the bands and to define JΣ is never defined, and the statement that diabolical-line modes 'cannot act as counter-propagating modes' is qualitative. In particular, the claim that the unscreened current JΣ≠0 'requires' a boundary effect is an inference, not a derived result. Since the boundary modes are already established by exact numerical diagonalization and non-Bloch theory, this section can be reframed as heuristic motivation; as written, it overstates the logical status of the anomaly. Please define Jk explicitly and either justify the screening statement or label it as a conjecture.","section":"Eq. (3) and the paragraph defining JΣ"}],"minor_comments":[{"comment":"The expression 'det(D~kbar−1 E)=0' appears to contain a typesetting artifact; it should read det(D~kbar − E)=0 or det(D~kbar − 1·E)=0.","section":"Eq. (5)"},{"comment":"The figure references in the discussion of the three regimes (see Fig. 1(a), Fig. 1(c), Fig. 1(d)) should be to Fig. 2(a), Fig. 2(c), and Fig. 2(d), since the relevant spectra and diabolical lines are displayed there.","section":"Paragraph after Eq. (2)"},{"comment":"The discretized real-space Hamiltonian used for the open-boundary numerical diagonalization is not described; please provide the lattice Hamiltonian, the boundary conditions, and the discretization scheme, or state that the code is available, so that the numerical results are reproducible.","section":"Numerical diagonalization, Fig. 3"},{"comment":"The phrase 'We usek0 =0.3' is missing a space; it should read 'We use k0 =0.3'.","section":"Fig. 1 caption"}],"recommendation":"major_revision","confidential_remarks":"The central mechanism and the non-Bloch results appear solid, and the spectral-winding problem is likely fixable in revision; I would not recommend rejection. The authors should be asked to supply a rigorous compactification prescription or to present the spectral winding as an informal indicator, and to make the status of the anomaly argument explicit."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague],\n\nQuick take: the paper demonstrates a genuinely new mechanism for non-Hermitian skin modes in a continuum bosonic system, backed by clean analytics and numerics, but the spectral winding number used to claim topological classification is not well defined as written. The core physics survives.\n\nThe new thing is the model itself: a 1D continuum of bosons with a Hermitian, momentum-shifted pairing term (i g a†_k a†_{-k-2k0} + h.c.) plus uniform local damping. The pairing is Hermitian and nonlocal in momentum, and the damping is just simple loss. The authors work out the spectrum, find a tilted diabolical line, and argue for a bulk anomaly from unbalanced current. On open boundaries they find skin modes numerically, and they reproduce the open-boundary spectrum and specific localization lengths (0.345 and 0.339 for g=γ=1, k0=0.3) using non-Bloch theory. Those analytical localization lengths matching numerics is the strongest evidence in the paper, and it is independent of the topological-index claim.\n\nThe soft spot is the winding number. In Eq. (4) they write the usual spectral winding after compactifying k with k̃ = arctan(k). But det(D_k - E) is a quartic polynomial in k with leading behaviour -k^4, so under compactification it has a fourth-order pole at the identified point k=±∞. The integral is not defined as an ordinary winding integral there; any result depends on a regularization that is not specified. In fact, computing the integral over k ∈ [-Λ,Λ] and taking Λ→∞ gives zero for generic E because the leading contributions from ±Λ cancel. The regions W=±1 in Fig. 3(a) are therefore not established by this calculation. This does not undermine the numerically observed skin modes or the non-Bloch spectrum, but it does mean the claimed bulk-boundary correspondence via spectral winding is not proven as written. A normalization like det under division by (1+k^2)^2 would likely fix it; the authors need to provide that.\n\nMinor things: the figure cross-referencing is sloppy (e.g., calling bulk spectra Fig. 1 when they are in Fig. 2), and the 'anomaly' argument is heuristic—useful as a picture, not a proof.\n\nVerdict: worth refereeing. The central result is likely correct and the winding issue is reparable. People working on non-Hermitian bosonic systems, parametric drives, or continuum topology will get value from this. I would send it to reviewers, expecting a revision that either repairs the winding-number definition or drops the topological-index claim and relies on the non-Bloch correspondence.","headline":"A genuinely new mechanism for non-Hermitian skin modes in a continuum bosonic system, with solid non-Bloch evidence, but the claimed spectral winding number is not well defined as written.","tokens_in":13071,"tokens_out":9707,"would_cite":true,"duration_ms":105454,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that a purely Hermitian, nonlocal parametric pairing process can induce non-Hermitian topology and skin modes in a continuous bosonic medium, with only uniform local loss added for stability.","keywords":["non-Hermitian topology","non-Hermitian skin effect","parametric driving","bosonic systems","continuum model","exceptional points","non-Bloch band theory","spectral winding number"],"falsifier":"Evaluate $W(E)=\\oint d\\tilde{k}\\,\\partial_{\\tilde{k}}\\log\\det(\\tilde{D}_{\\tilde{k}}-E)/(2\\pi i)$ on truncated domains $[-L,L]$ for fixed $E$ in the predicted $\\pm1$ regions; if the winding changes with cutoff or depends on a regularization of the $k\\to\\pm\\infty$ phase, the topological index is not well defined. Alternatively, realize the model in a driven bosonic gas with local loss and measure the steady-state density: absence of exponential accumulation at the predicted edge for $|k_0|<\\sqrt{|g|}$ would falsify the bulk-boundary claim.","tokens_in":12016,"feed_emoji":"🌀","tokens_out":5279,"duration_ms":63213,"temperature":0.7,"pith_summary":"This paper claims that a purely Hermitian, nonlocal parametric pairing term—bosons created and annihilated in pairs with opposite momenta plus a fixed momentum shift—can generate non-Hermitian topology and the non-Hermitian skin effect in a one-dimensional continuum, with only uniform local loss added for stability. The pairing produces level attraction between particle and hole branches, and beyond a threshold it creates a tilted diabolical line in the complex spectrum. The tilt supplies a net current that cannot be screened by the remaining bulk modes, a bulk anomaly that forces boundary modes when the system is opened. The authors compute a spectral winding number after compactifying the unbounded momentum axis and confirm the resulting skin modes and their localization lengths using non-Bloch theory. If the claim holds, parametric drives become a minimal, lattice-free route to non-Hermitian topological behavior in bosonic systems.","feed_headline":"A Hermitian pairing process produces non-Hermitian skin modes","feed_subtitle":"Momentum-shifted pair creation plus local loss localizes modes at edges, with a topological index to back it up.","key_machinery":"The machinery is the momentum-shifted Nambu dynamical matrix $\\tilde{D}_k=\\sigma_z H_k-i\\gamma\\mathbb{1}$, with off-diagonal pairing $ig$ that couples $k-k_0$ particles to $-k-k_0$ holes. The imaginary pairing causes level attraction, and for $k_0\\neq 0$ the resulting diabolical line—the one-dimensional set in the complex spectrum where the two bands touch in real part, bounded by exceptional points—is tilted, creating the unpaired current. The spectral winding number is made well defined by compactifying momentum via $\\tilde{k}=\\arctan(k)$, and non-Bloch theory with complex wave number $\\bar{k}$ provides the generalized Brillouin zone and analytic localization lengths. Together these convert a local continuity anomaly into a quantized bulk index and edge modes.","core_discovery":"The central result is that the non-Hermitian skin effect in this model is not caused by engineered non-Hermitian couplings but by the Hermitian coupling between momentum-shifted particle and hole modes. In the shifted Nambu basis, the spectrum is $E=-2kk_0\\pm\\sqrt{(k^2+k_0^2)^2-g^2}$; for $|k_0|<\\sqrt{|g|}$ the particle and hole bands coalesce into a tilted diabolical line terminated by exceptional points. Uniform loss $\\gamma$ moves all lifetimes to the stable side, $\\gamma\\ge\\sqrt{g^2-k_0^4}$, while the tilt of the diabolical line produces an unscreened bulk current. Opening the boundaries yields exponentially localized edge modes, classified by a spectral winding number $W(E)=\\pm1$ computed after mapping $k\\in(-\\infty,\\infty)$ to a circle with $\\tilde{k}=\\arctan(k)$; non-Bloch theory reproduces the open-boundary spectrum and localization lengths exactly.","pith_inferences":["A natural extension is to read the pair creation as a squeezed-drive term: a Josephson or optical parametric oscillator with a finite-momentum pump should show the same edge accumulation, testable in existing tabletop setups.","Since the skin effect here comes from Hermitian pairing rather than loss engineering, it may appear in effectively Hermitian Bogoliubov descriptions of interacting condensates whenever a small local loss is added; that would move the effect from specially built lattices to ordinary superfluids.","One could test the anomaly directly by measuring the bulk current $J_\\Sigma$ as a function of $k_0$: the predicted non-monotonic signal and its sign change should correlate with the onset of edge modes.","The paper does not address interactions; strong interactions could either stabilize or destroy the tilted diabolical line, so a generalized non-Bloch analysis with a Hartree term is a concrete next test."],"forward_implications":["For $|k_0|<\\sqrt{|g|}$ and damping $\\gamma\\ge\\sqrt{g^2-k_0^4}$, the system is stable and every bulk mode has positive lifetime, so the skin effect can be observed without gain.","The spectral winding number predicts right- or left-edge localization from the sign of $k_0$, giving directional control by reversing the momentum shift.","Non-Bloch theory yields analytic localization lengths, so edge-mode decay rates are quantitatively predictable from bulk parameters.","Because the construction is continuum and lattice-free, parametric drives in ultracold atoms, polaritonic condensates, or acoustic metamaterials are candidate platforms.","The same logic extends the bulk-boundary correspondence for non-Hermitian systems from tight-binding lattices to continuum models."],"supporting_citations":[{"why":"Supplies the finite-momentum pairing concept that shifts particle and hole momenta by $2k_0$.","marker":"[52]"},{"why":"Provides the level-attraction mechanism from which the diabolical line emerges.","marker":"[54]"},{"why":"Justifies compactifying unbounded momentum to apply topological invariants in continuous media.","marker":"[55]"},{"why":"Defines the spectral winding number used to classify the skin effect.","marker":"[26]"},{"why":"Supplies non-Bloch band theory used to restore the bulk-boundary correspondence.","marker":"[29]"},{"why":"Extends non-Bloch theory to continuum systems, grounding its application here.","marker":"[31]"}],"fun_headline_variants":["Non-Hermitian skin modes from Hermitian pairing","Hermitian pairing spawns non-Hermitian skin modes","Skin modes from Hermitian pairing, no bath engineering","Hermitian pairing induces non-Hermitian skin effect","No complex baths: Hermitian pairing yields skin modes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that compactifying the unbounded momentum axis by $\\tilde{k}=\\arctan(k)$ makes the spectral winding integral well defined: the phase of $\\det(\\tilde{D}_k-E)$ must approach a controlled value as $k\\to\\pm\\infty$, with no boundary contributions, even though that convergence is not proven in the paper.","fun_headline_variants_meta":{"raw":{"variants":["Non-Hermitian skin modes from Hermitian pairing","Hermitian pairing spawns non-Hermitian skin modes","Skin modes from Hermitian pairing, no bath engineering","Hermitian pairing induces non-Hermitian skin effect","No complex baths: Hermitian pairing yields skin modes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001427,"raw_usage":{"total_tokens":5716,"prompt_tokens":864,"completion_tokens":4852,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":480,"completion_tokens_details":{"reasoning_tokens":4775}},"tokens_in":480,"tokens_out":4852,"duration_ms":35066,"temperature":1.0,"reasoning_tokens":4775,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:41:04.985089+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate $W(E)=\\oint d\\tilde{k}\\,\\partial_{\\tilde{k}}\\log\\det(\\tilde{D}_{\\tilde{k}}-E)/(2\\pi i)$ on truncated domains $[-L,L]$ for fixed $E$ in the predicted $\\pm1$ regions; if the winding changes with cutoff or depends on a regularization of the $k\\to\\pm\\infty$ phase, the topological index is not well defined. Alternatively, realize the model in a driven bosonic gas with local loss and measure the steady-state density: absence of exponential accumulation at the predicted edge for $|k_0|<\\sqrt{|g|}$ would falsify the bulk-boundary claim.","supporting_citations":[{"cited_title":"Gardin, G","cited_arxiv_id":null,"evidence_quote":"Provides the level-attraction mechanism from which the diabolical line emerges."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Justifies compactifying unbounded momentum to apply topological invariants in continuous media."},{"cited_title":"Hu, Y.-Q","cited_arxiv_id":null,"evidence_quote":"Extends non-Bloch theory to continuum systems, grounding its application here."}],"review_version":1}