{"id":"52563b6c-84a7-4a33-9121-a0ac2809d37f","arxiv_id":"2412.08546","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A repulsively coupled impurity in a finite-temperature Bose gas develops an attractive spectral line when it binds to vortices and density holes, providing a potential local probe of vortex proliferation across the BKT transition in 2D and below the BEC transition in 3D.","lead":"This paper simulates an impurity particle in a warm, weakly interacting gas of bosons and finds that the impurity's spectrum develops extra low-energy lines when vortices appear in the gas. The result suggests a new way to spot vortex formation in flat superfluids, such as excitonic insulators, where vortices are hard to detect directly.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Impurity couples only to density, not phase: the attractive branch is a generic density-hole/Lifshitz-tail effect, so the vortex-specific interpretation and the T_BKT minimum are not yet established.","rationale":"I read the conditional verdict as appropriate. The numerical implementation is careful, finite-size checks are provided, and the comparison with static vortex bound-state energies gives independent support for the binding mechanism. The central claim, however, is about detecting vortices specifically, and that is the least secure part of the argument. The impurity Hamiltonian is a density probe; without a control that separates vortex cores from non-topological density holes, the attractive branch is equally explained by impurity localization in the deepest density fluctuations. The cutoff dependence of n_inc compounds this by making the apparent BKT minimum sensitive to a numerical parameter. Both issues are testable and would not necessarily overturn the core numerical finding, so I would keep the verdict CONDITIONAL (unchanged).","tokens_in":16977,"tokens_out":7620,"duration_ms":85138,"concrete_test":"Compute the 2D impurity spectrum at T ~ 1.1 T_inf_BKT for two modified baths: (i) the same equilibrated |psi|^2 but with vortex cores smoothly filled in (an approximately vortex-free density field), and (ii) a truly vortex-free positive field with matched density power spectrum. If the attractive branch persists with comparable weight and position, the vortex-specific interpretation is unsupported. Separately, rerun with epsilon_cut=mu+2T log2 (and corresponding n_inc) and check whether the branch minimum shifts by more than 5% of g_BI n_0; a shift would confirm the cutoff artifact.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. (4) makes the impurity respond exclusively to n_B(x,t)=|psi|^2+n_inc; the phase of psi never enters. Vortices therefore affect the impurity only through the density depletion at their cores, and the paper's own snapshots (Fig. S4/S5) show that many density holes are not vortices. The low-energy branch is thus a bound state in generic low-density regions; no calculation isolates the vortex contribution or shows that the branch weight/position tracks vortex density rather than, say, the depth distribution of thermal density holes. The minimum near T_inf_BKT is also not a clean observable: below T_BKT the branch redshifts as deep holes appear, while above T_BKT it blueshifts because n_inc=(mT/2pi)log2 grows linearly with T. The cutoff epsilon_cut=mu+T log2 is conventional but arbitrary; a different cutoff changes n_inc and therefore the location/strength of the minimum. So the quantitative claim that spectroscopy can estimate T_BKT is cutoff-dependent until shown otherwise.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses the stochastic projected Gross-Pitaevskii equation (SPGPE) for the weakly interacting Bose bath and a full Schrödinger evolution for the impurity to compute the zero-momentum injection spectrum of a repulsively coupled impurity at finite temperature. In 2D, the authors find a low-energy attractive branch whose weight grows near the BKT transition and whose minimum lies close to T_BKT; they attribute this branch to impurity binding to vortices and density holes. In 3D, they report an attractive branch associated with vortex rings below the BEC transition and a single Lorentzian above it. The numerical setup includes density fixing via the chemical potential, finite-size checks, and a comparison of the attractive branch with bound-state energies in an isolated static vortex.","tokens_in":17103,"tokens_out":3892,"duration_ms":46782,"significance":"If the vortex-specific interpretation is established, the paper would provide a concrete local spectroscopic route to detecting vortex proliferation in neutral superfluids, with potential relevance to excitonic insulators and to finite-temperature Bose polaron physics. The numerical work is carefully designed: the bath is equilibrated within SPGPE, the density is fixed, finite-size convergence of local observables and of the spectra is reported, and the bound-state comparison uses an independent static vortex profile (Pitaevskii-Berloff). No parameter is fitted to the target spectral data. The main weakness is that the central attribution of the attractive branch to vortices, as opposed to generic density holes, is not directly demonstrated.","major_comments":[{"comment":"The impurity Schrödinger equation (4) depends on the bosonic density n_B(x,t)=|psi|^2+n_inc and not on the phase of the classical field psi. Consequently, a vortex affects the impurity only through the density depletion at its core, and any low-density region produces the same kind of attractive potential. The paper itself states that 'not all density holes correspond to vortices', and Figs. S4-S5 show density deeps without any phase winding. The attractive branch is therefore equally consistent with a generic density-hole or Lifshitz-tail mechanism. To support the central claim that spectroscopy detects vortex proliferation, the authors should provide a diagnostic that isolates the vortex contribution, for example by conditioning the spectrum on the impurity being initialized at detected vortex cores, or by recomputing the spectrum after removing phase windings while preserving the density field. Without such a test, the vortex-specific interpretation is underdetermined.","section":"Model and methods, Eq. (4); Numerical results in 2D; Figs. S4-S5"},{"comment":"The presence and location of the minimum of the attractive branch near T_inf_BKT depend in part on the incoherent density n_inc=(mT/2pi)log2, whose numerical value is set by the cutoff epsilon_cut=mu+T log2. The authors note that the blueshift above T_BKT is 'in part due to the increase of the incoherent fraction', which is a uniform contribution. A different, equally conventional choice of cutoff would change n_inc and therefore shift the position and strength of the apparent minimum. The quantitative claim that the spectra can be used to estimate T_BKT is thus cutoff-dependent unless the authors demonstrate robustness by varying epsilon_cut or by subtracting the uniform n_inc contribution before extracting the minimum.","section":"Numerical results in 2D, Fig. 1(b); Eq. (5)"},{"comment":"For the L=20 micron box used in the main text, the algebraic correlation exponent alpha reaches the BKT critical value 0.25 at T approx 1.25 T_inf_BKT, which the authors attribute to the known finite-size shift of the BKT transition. However, the attractive-branch minimum is reported near T_inf_BKT, not near the finite-size crossover visible in the same simulations. This discrepancy needs an explicit explanation. If the minimum tracks the thermodynamic T_BKT while the actual finite-size vortex proliferation occurs at higher temperature, the vortex-binding mechanism becomes harder to reconcile with the observed minimum; if instead the minimum is controlled by generic density fluctuations, that also weakens the central claim.","section":"Numerical results in 2D, Fig. 1(b); Supplemental Fig. S6(b)"},{"comment":"The same density-only coupling issue applies to the 3D vortex-ring interpretation. The vortex-ring density in Fig. 3(a) grows below TBEC, but thermal density fluctuations also grow in the same temperature range, and the impurity potential (4) does not distinguish a ring core from any other low-density tubular region. The observation that the spectrum becomes a single Lorentzian above TBEC is expected for any nearly homogeneous thermal gas and does not by itself prove that the sub-TBEC attractive branch is caused by vortex rings. A control calculation that removes or masks the phase information, or a correlation analysis between the spectrum and ring-core locations, is needed before the 3D claim can be considered established.","section":"Numerical results in 3D, Fig. 3(b)"}],"minor_comments":[{"comment":"There is a typo: 'irrelavant' should be 'irrelevant'.","section":"Model and methods, Eq. (1)"},{"comment":"The phrase 'above the cutoff blueepsilon_cut' contains an apparent editing artifact ('blue'); it should simply read 'above the cutoff epsilon_cut'.","section":"Model and methods, after Eq. (5)"},{"comment":"The number of stochastic realizations Mstat is not specified in the main text, and the spectral plots do not show statistical error estimates; this information should be provided or referenced to the Supplemental Material.","section":"Model and methods"},{"comment":"Reference [68] is given as a placeholder string rather than a complete citation; the actual Supplemental Material URL should be supplied.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The numerical study is competently executed and the raw spectral phenomenology is interesting. However, the paper's headline claim that the attractive branch provides a direct spectroscopic signature of vortex proliferation is not yet isolated from a generic density-hole effect, because the impurity couples only to density. The finite-size shift of the BKT exponent versus the location of the spectral minimum is an additional inconsistency that the authors should address explicitly. I would encourage the editor to send the manuscript back for a major revision rather than reject it, since the missing diagnostic (conditioning on vortices, removing phase information, or varying the cutoff) appears feasible within the manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the goods: this is a careful numerical study. The authors use SPGPE for the bath and the full Schrödinger equation for the impurity, equilibrate properly, fix density, and check finite-size convergence. The comparison of the attractive branch position to a T=0 bound state in an isolated vortex (via Pitaevskii-Berloff profile) is a nice independent check. The 3D result — vortex rings well below T_BEC, visible in projected vortex detection — is new and interesting. The work is worth serious attention.\n\nThe soft spot is the interpretation. The impurity equation (Eq. 4) couples only to the total density n_B = |psi|^2 + n_inc; the phase never appears. Vortices affect the impurity only through the density dip at their cores, and the paper's own snapshots show many density holes with no vortex. So the low-energy branch is a bound state in generic low-density regions. Nothing in the calculation isolates the vortex contribution or shows that the branch weight tracks vortex density rather than the distribution of thermal hole depths. The abstract says 'direct spectroscopic signatures of vortex proliferation' — that overstates what is shown.\n\nThe T_BKT minimum is also less clean than it looks. Below T_BKT the branch redshifts as deep holes appear; above T_BKT it blueshifts partly because n_inc = (mT/2pi) log 2 grows linearly with T. The energy cutoff epsilon_cut = mu + T log 2 is conventional but not unique; changing it changes n_inc and therefore the position and strength of the minimum. So the quantitative claim that this spectroscopy can estimate T_BKT is cutoff-dependent until the authors show otherwise. The lack of error bars on the spectra doesn't help, though the finite-size checks soften that.\n\nNone of this sinks the central numerical finding: a repulsively coupled impurity develops a low-energy branch as thermal density holes appear, and the branch energy is consistent with binding to vortices. It is a good paper, but the interpretation needs sharpening. A referee should ask for a calculation that separates vortices from non-vortex holes — e.g., masking out vortices or correlating the bound-state weight with vortex density versus hole depth statistics — and for a scan over the cutoff to quantify the robustness of the T_BKT minimum.\n\nI'd send it to peer review; it deserves referee time, but it needs revision before publication.","headline":"Careful SPGPE-plus-Schrödinger numerics give a real temperature-dependent attractive branch in impurity spectra, but the vortex-specific reading and the BKT-estimation claim outrun what the calculation actually shows.","tokens_in":17701,"tokens_out":2771,"would_cite":true,"duration_ms":28995,"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":"A repulsive impurity's injection spectrum develops a redshifted branch that tracks vortex proliferation across the BKT and BEC transitions, offering a local spectroscopic probe of these phase transitions.","keywords":["Bose polaron","injection spectroscopy","Berezinskii-Kosterlitz-Thouless transition","vortex proliferation","stochastic Gross-Pitaevskii equation","classical-field methods","excitonic insulator","Bose-Einstein condensation"],"falsifier":"Measure the injection spectrum of a repulsive impurity in a uniform 2D Bose gas while independently imaging the vortex distribution (e.g., by matter-wave interference or stirring); if the low-energy sideband does not appear at the temperature where vortices proliferate, or if its peak position disagrees with the bound-state energy computed from the measured vortex density profile by more than the linewidth, the vortex-binding mechanism is falsified.","tokens_in":16706,"feed_emoji":"🌀","tokens_out":11091,"duration_ms":99408,"temperature":0.7,"pith_summary":"This paper claims that a single repulsive impurity injected into a weakly interacting Bose gas at finite temperature can act as a local spectroscopic detector of vortex proliferation. In two dimensions, the injection spectrum develops a low-energy attractive branch whose weight grows as the temperature crosses the Berezinskii-Kosterlitz-Thouless transition, with a minimum around $T_{\\rm BKT}$; the branch is attributed to the impurity binding to vortices and density holes, where the reduced boson density lowers the repulsive Hartree energy. In three dimensions, an attractive branch appears once vortex rings nucleate well below the BEC temperature and disappears above it, leaving a single Lorentzian from the uniform thermal gas. The authors propose this mechanism as a spectroscopic estimator of the BKT transition in excitonic insulators and as a route to understanding finite-temperature Bose polaron spectra.","feed_headline":"An impurity's red-shifted line detects vortices near BKT transition","feed_subtitle":"A repulsive impurity binds to density holes, giving a spectroscopic marker of vortex proliferation in 2D and 3D.","key_machinery":"The central machinery is the stochastic projected Gross-Pitaevskii equation (SPGPE) for the Bose bath, coupled to a full Schrödinger equation for the impurity moving in the stochastic potential $g_{BI} n_B(\\mathbf{x},t)$. The SPGPE describes the highly occupied modes of the gas up to an energy cutoff $\\epsilon_{\\rm cut} = \\mu + T \\log 2$, with a projected white noise enforcing thermal fluctuations. For each stochastic realization, the impurity wavefunction is evolved from a uniform initial state, and the injection spectrum is obtained from the Fourier transform of the averaged overlap function. The physical mechanism is that the impurity is repelled from high-density regions and can become quasi-bound in the density-depleted core of a vortex; the binding energies are computed by exact diagonalization in the Padé-approximated vortex density profile, and they match the observed spectral sidebands.","core_discovery":"The authors show that the injection spectrum of a zero-momentum impurity repulsively coupled to a weakly interacting Bose gas carries direct, local signatures of vortex proliferation. At low temperature the spectrum is a delta-like peak at the Hartree shift $g_{BI} n_0$. As temperature rises and vortices (or vortex rings in 3D) appear, the impurity can bind in the density-depleted core of a vortex, producing a red-shifted spectral branch whose minimum lies at the BKT transition in 2D and at the threshold of vortex-ring nucleation in 3D. Above the BEC transition in 3D the condensate is gone, vortices vanish, and the spectrum collapses to a single Lorentzian. The energies of the observed attractive lines match the s-wave bound states of the impurity in an isolated vortex computed from the Padé-approximated Pitaevskii profile.","pith_inferences":["If this mechanism is confirmed, injection spectroscopy could be used as a non-destructive, local sensor of topological order, complementing global probes such as the superfluid density jump.","A natural testable extension is to measure the spatial map of the attractive-branch weight in a trapped gas and check that it tracks the local vortex density.","One could apply the same stochastic approach to systems where vortices proliferate out of equilibrium (driven-dissipative polariton condensates or quenched gases), extending the detection scheme beyond equilibrium transitions.","The neglect of impurity back-action sets a practical limit: at higher impurity densities or stronger repulsion, self-localization or bath deformation would renormalize the binding energy, and the simple single-impurity picture would need correction."],"forward_implications":["A repulsive impurity can serve as a local spectroscopic probe of vortex proliferation, providing an estimate of the BKT transition temperature in 2D without imaging phase winding or measuring superfluid density.","The same mechanism could be used to detect superfluidity and vortices in excitonic insulators in transition-metal dichalcogenide heterostructures.","In 3D, the impurity spectrum distinguishes the vortex-ring regime below $T_{\\rm BEC}$ from the normal-gas regime above it, potentially allowing a local thermometer for the condensation transition.","The mass dependence of the attractive sidebands (one vs two lines) reflects the number of vortex-bound states, making the spectroscopy sensitive to the vortex density profile.","Finite-temperature Bose polaron injection spectra generically feature a split repulsive branch and an attractive branch whenever vortices or density holes are present."],"supporting_citations":[{"why":"Defines the SPGPE framework and the energy cutoff $\\epsilon_{\\rm cut} = \\mu + T \\log 2$ used to equilibrate the finite-temperature bath.","marker":"[47]"},{"why":"Introduces the stochastic Gross-Pitaevskii equation that underlies the bath dynamics.","marker":"[44]"},{"why":"Provides the weakly-interacting BKT temperature formula used to place the critical point in 2D.","marker":"[66]"},{"why":"Supplies the ideal vortex density profile whose bound states explain the attractive spectral branch.","marker":"[81]"},{"why":"Gives the Padé approximation of the vortex profile used to compute impurity-vortex binding energies.","marker":"[82]"},{"why":"Provides the vortex detection algorithm used to count vortices and correlate them with spectra.","marker":"[69]"},{"why":"Extends the vortex detection method used in the 3D vortex-ring analysis.","marker":"[70]"},{"why":"Establishes the BKT vortex-unbinding picture that the 2D probe is designed to detect.","marker":"[43]"},{"why":"Prior SPGPE study of the 2D BKT transition, validating the classical-field approach for this regime.","marker":"[58]"}],"fun_headline_variants":["Quantum impurity spies vortices via red-shifted spectral line","Impurity binds to vortices, marking BKT transition in 2D","Vortex proliferation read out by impurity spectra in Bose gases","Red-shift in impurity spectrum signals vortex nucleation","Probing BKT transition with a single impurity's spectrum"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the stochastic classical-field description of the Bose gas, together with the neglect of the impurity's back-action on the bath, faithfully reproduces the density holes and vortex configurations that the impurity binds to; if those fluctuations are captured incorrectly, or if back-action becomes significant at the interaction strengths used, the predicted attractive sideband would shift or disappear.","fun_headline_variants_meta":{"raw":{"variants":["Quantum impurity spies vortices via red-shifted spectral line","Impurity binds to vortices, marking BKT transition in 2D","Vortex proliferation read out by impurity spectra in Bose gases","Red-shift in impurity spectrum signals vortex nucleation","Probing BKT transition with a single impurity's spectrum"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000444,"raw_usage":{"total_tokens":2233,"prompt_tokens":920,"completion_tokens":1313,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":536,"completion_tokens_details":{"reasoning_tokens":1230}},"tokens_in":536,"tokens_out":1313,"duration_ms":10666,"temperature":1.0,"reasoning_tokens":1230,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T17:44:11.304753+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the injection spectrum of a repulsive impurity in a uniform 2D Bose gas while independently imaging the vortex distribution (e.g., by matter-wave interference or stirring); if the low-energy sideband does not appear at the temperature where vortices proliferate, or if its peak position disagrees with the bound-state energy computed from the measured vortex density profile by more than the linewidth, the vortex-binding mechanism is falsified.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the stochastic Gross-Pitaevskii equation that underlies the bath dynamics."},{"cited_title":"Prokof’ev and B","cited_arxiv_id":null,"evidence_quote":"Provides the weakly-interacting BKT temperature formula used to place the critical point in 2D."},{"cited_title":"Pitaevskii, Vortex lines in an imperfect Bose gas, Sov","cited_arxiv_id":null,"evidence_quote":"Supplies the ideal vortex density profile whose bound states explain the attractive spectral branch."},{"cited_title":"Quantum impurities in finite-temperature Bose gases: Detecting vortex proliferation across the BKT and BEC transitions","cited_arxiv_id":null,"evidence_quote":"Gives the Padé approximation of the vortex profile used to compute impurity-vortex binding energies."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Extends the vortex detection method used in the 3D vortex-ring analysis."}],"review_version":1}