{"id":"bc4cefb0-df2a-4852-b683-ffc7b572bc0e","arxiv_id":"2504.15840","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A linear molecule in a superfluid can carve a density dimple that supports a bound attractive angulon state, and rotational spectroscopy should show a repulsive-to-attractive crossover as bath density rises.","lead":"A theory paper predicts that a molecule spinning inside superfluid helium can carve a small dimple in the surrounding fluid, and that this dimple can trap a new kind of rotational quasiparticle called the attractive angulon. The same calculation predicts a density-driven crossover between two distinct angulon branches that should be visible in rotational spectroscopy.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed B*<B resolution rests on E_att from the single-excitation ansatz, but the paper's own spectroscopy puts the AA peak above E_att; no calculation shows the observable peak or a correlation-corrected E1-E0 remains below 2B.","rationale":"The reader's conditional verdict is appropriate. The most load-bearing weakness is not merely the absence of a convergence check, but the internal inconsistency between the quantity used for the headline B*<B statement (E_att from the single-excitation ansatz) and the quantity the paper itself identifies as the spectroscopically observed attractive-angulon feature (a higher-frequency AA peak). Since the conclusion explicitly claims to resolve the anomalous rotational constant problem, this gap matters. A second, related issue is that B* is a difference of two variational energies from different approximation levels; the variational upper-bound property of each energy individually does not protect the difference, so quantum depletion could in principle push E1-E0 above 2B. These are addressable with a more complete variational or numerically exact calculation, and the paper contains no benchmark against Monte Carlo or experiment. The existence of an attractive angulon state below the single-phonon continuum is supported by the variational calculation, so the concern does not warrant rejection; it warrants conditional acceptance pending a concrete convergence or benchmark test. The reader's rationale already noted the E_att/spectrum discrepancy and the truncation issue, so my assessment agrees partially rather than fully with the stated weakest assumption.","tokens_in":19135,"tokens_out":15624,"duration_ms":165298,"concrete_test":"Perform a systematic improvement of the variational space in the same model: include two-phonon states (|1,m> b†b†|Φ0>) and allow the condensate wavefunction f_M to be optimized independently in the J=0 and J=1 sectors, for the same parameters and box size. Compute ΔE = E1 - E0 and locate the spectral AA peak. If ΔE >= 2B or the AA peak lies at or above 2B at intermediate densities (e.g., ρ = 3.5e-2), the central B*<B claim fails. A complementary check is exact diagonalization for small boson number N with the same pseudopotential to benchmark the energy difference.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section VI defines B*_1 = (E_att - E_GS)/2 and concludes B*_1 < B because E_att - E_GS < 2B for all densities. However, E_att is the lowest eigenvalue of the single-Bogoliubov-excitation effective Hamiltonian of Eq. (12), while Section V explicitly states that the AA peak in the multireference rotational spectrum appears at a frequency higher than this E_att, attributing the discrepancy to interactions among multiple Bogoliubov excitations. The paper never demonstrates that this observable AA peak, or the energy of the true J=1 ground state, remains below 2B. Moreover, B* is a difference of two variational energies: E_GS is obtained from a single coherent state, and E_att from a one-phonon ansatz built on that coherent state. Even though E_att is an upper bound to the exact J=1 ground-state energy and E_GS is an upper bound to the exact J=0 energy, these individual bounds do not bound their difference. If quantum depletion lowers the J=0 energy more than correlations lower the J=1 energy, the true E1-E0 could exceed 2B and the anomalous rotational constant would not be resolved. No convergence study, Monte Carlo benchmark, or independent spectrum is supplied to rule out this possibility.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a linear rotor immersed in a bosonic superfluid, using a Lee-Low-Pines-like transformation followed by a multireference variational ansatz in which each rotor projection state is entangled with a coherent bath state; excited states are obtained by diagonalizing an effective Hamiltonian in a single-Bogoliubov-excitation subspace built on the self-consistently deformed condensate. The central claims are: (i) the rotor's backaction creates a superfluid density defect; (ii) at intermediate and high densities this defect supports a bound 'attractive angulon' below the Bogoliubov continuum, while dilute systems exhibit a long-lived 'repulsive angulon'; (iii) rotational spectral functions show a density-driven crossover between these states; and (iv) the effective rotational constant B*_1, defined through the attractive angulon energy, is always below the bare rotor value, thereby 'resolving' the anomalous moment-of-inertia problem in earlier Chevy-ansatz treatments. Appendices A and B provide the variational equations of motion and the Bogoliubov effective Hamiltonian used in the numerics.","tokens_in":19349,"tokens_out":6251,"duration_ms":62509,"significance":"If established, the attractive angulon would be a qualitatively new quasiparticle: a rotor bound to a localized superfluid density wave, with useful analogies to photonic bound states and YSR states. The combination of the LLP transformation with a deformed, nonuniform condensate background goes beyond the uniform-background Chevy ansatz of Ref. [50] and yields explicit dynamical predictions (spectral peak positions, linewidths, and density evolution) that are in principle falsifiable in ultracold-gas experiments. A clear strength is the level of detail in Appendices A and B, which makes the variational computation reproducible. The main caveats are that the conclusions rest on a restricted variational manifold with no reported convergence checks, and the B*_1 < B resolution uses an energy from a smaller Hilbert space than the spectroscopic peak to which it is compared. These gaps are load-bearing for the paper's central claims, but they are addressable in a revision.","major_comments":[{"comment":"The proof that B*_1 < B is not established. Equation (16) defines B*_J through the lowest energy in the J sector, but the argument in Sec. VI uses E_att from the single-Bogoliubov-excitation effective Hamiltonian of Eq. (12), while Sec. V explicitly states that the AA peak in the multireference rotational spectrum lies at a higher frequency than that E_att, attributing the difference to interactions among multiple Bogoliubov excitations. Because E_GS and E_att are separately variational upper bounds, their difference is not an upper bound on the exact E_1 - E_0, and the paper never shows that the observable AA peak or a correlation-corrected J=1 ground-state energy remains below 2B. Additionally, in the dilute regime the manuscript evaluates B*_1 using the metastable repulsive angulon energy E_rep, which is not the lowest J=1 eigenstate; using a resonance energy in a ground-state definition of the rotational constant requires justification. The authors should compute E_1 - E_0 within the multireference ansatz (or an extended manifold) and show that it remains below 2B, or redefine B* through the observed spectroscopic peak and discuss the finite-lifetime caveat.","section":"Sec. VI, Eq. (16); Sec. V"},{"comment":"The admitted boundary artifact undermines the dilute-regime part of the claims. The text states that in the dilute regime the density distribution rho_1m(r) is localized at the hard-wall edge and 'is merely a boundary effect.' Nevertheless, Fig. 2(c) reports E_att and Z_att in this regime as physical, and Sec. VI uses E_att - E_GS < 2B 'for all densities' to conclude B*_1 < B. Thus the low-density attractive angulon is not established as a state localized within the density defect, and its use in the crossover scenario or in the all-density B* statement is unsupported. A finite-size study (varying R and the boundary condition) and a clear separation of physical bound states from edge-localized states are needed before the dilute-regime spectroscopy can be interpreted as repulsive-angulon physics with an attractive angulon present but dark.","section":"Sec. IV, Fig. 2(d); Sec. VI"},{"comment":"No convergence checks are provided for the variational truncations, and the central existence of the attractive angulon depends on this truncation. The ground state is represented by one coherent state per rotor projection M (Eq. (7)), and the J>0 states are obtained in a subspace with a single Bogoliubov excitation on top of that coherent state (Eq. (12)); the explicit Bogoliubov matrices in Appendix B are written for angular-momentum cutoff l_c = 1. The manuscript does not report tests against larger l_c, higher momentum cutoffs k_c, larger system radii R, or multi-phonon sectors, nor any benchmark against Monte Carlo data. Since the density defect and the claimed shape-resonance bound state are consequences of this truncated manifold, the authors should show that E_att, Z_att, and the spectral features are stable when the variational space is enlarged; otherwise the attractive angulon may be an artifact of the ansatz.","section":"Sec. III; Sec. IV; Appendix B, Eqs. (7), (12)"}],"minor_comments":[{"comment":"Equation (15) contains an unbalanced parenthesis: rho_lm(r,t) is defined with a trailing ')' after |psi_J(t)> that has no matching opening parenthesis.","section":"Eq. (15)"},{"comment":"The phrase 'the healing timethealing~1/(rho abb)' is missing a space and the quantity is not defined elsewhere; please define the healing time and use consistent notation.","section":"Sec. V, Fig. 4 discussion"},{"comment":"The notation in Eq. (13) is opaque: k0 = alpha_1^(1)/R uses the spherical-Bessel zero introduced later in Appendix B, but the equation does not define the normalization of the state or the integral measure over the rotor Euler angles; please clarify.","section":"Eq. (13) and Appendix B"},{"comment":"The sentence 'It is important to note that ˆJ 2 = ˆJ2' is typographically confusing; the intended identity should be typeset clearly (e.g., hat J^2 = hat J^2 or a comment on the body-fixed vs laboratory frame) to avoid ambiguity.","section":"Sec. II, after Eq. (6)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a good fit for a quantum-gas/many-body journal and the variational framework is competently presented. My main concern is that the headline conceptual claims (attractive angulon, B*_1 < B resolution) are supported only within a narrow variational manifold and a specific definition of the energy that is not the one seen in the reported spectroscopy. The authors should be pushed to provide convergence checks and to reconcile Eq. (16) with the multireference spectral peaks. The low-density boundary artifact needs to be quantified, not just acknowledged. These are fixable within the manuscript's scope, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nHere's my take on 2504.15840. The genuinely new thing is the self-consistent density defect: a linear rotor deforming its own condensate, and a bound \"attractive angulon\" state living in that defect. Earlier angulon work relied on a uniform condensate with a Chevy-like ansatz; this paper lets the bath adjust and finds a bound state below the Bogoliubov continuum, with a crossover in the rotational spectrum as density increases. The time-dependent calculation and the spectroscopy are also more than static ground-state stuff—the multiphonon emission and relaxation into the defect-bound state is a substantial piece of work.\n\nThe soft spots are real, though not fatal. The biggest one is the claimed resolution of the anomalous rotational constant. The paper defines B*_1 from E_att, the lowest eigenvalue of a single-phonon ansatz, and notes that E_att - E_GS < 2B. But the spectroscopy they compute with the multireference ansatz puts the AA peak at a higher frequency than E_att, and they never show that the observable peak—or the true J=1 ground-state energy in the larger variational space—remains below 2B. Since E_GS and E_att are variational upper bounds on the J=0 and J=1 ground states, their difference is not a bounded estimate of the exact gap. If quantum depletion lowers the J=0 energy substantially, the true gap could exceed 2B. This is a load-bearing gap in the anomaly-resolution claim, not a cosmetic issue.\n\nSecond, there are no convergence checks for the angular-momentum cutoff, momentum cutoff, or box size; no code or data; and no benchmark against Monte Carlo or experimental rotational spectra. The dilute-regime boundary artifact is admitted, which is honest, but it means that part of the phase diagram is not yet under control.\n\nThese are all addressable. The central mechanism—the density defect hosting a bound state—is physically plausible and new, and the calculation is coherent. But the paper oversells the B* resolution relative to what its own spectroscopy shows.\n\nWho is this for? People working on angulons, helium-droplet spectroscopy, and Bose polarons. It deserves a serious referee; I'd send it to review but with a clear request to fix the B* argument and add convergence checks.","headline":"The attractive angulon and the density-defect mechanism are a real step beyond earlier angulon work, but the claimed B*<B resolution is not backed by the paper's own spectroscopy.","tokens_in":19955,"tokens_out":5644,"would_cite":true,"duration_ms":54704,"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":"A rotating molecule in a superfluid creates a density defect that binds a new class of rotor-bath states, the attractive angulon, while dilute baths host long-lived repulsive angulons.","keywords":["angulon","superfluid density defect","rotational spectroscopy","quantum impurity","Lee-Low-Pines transformation","multireference variational ansatz","Bogoliubov excitations","bound states"],"falsifier":"Compute the $J=1$ spectral function with a larger variational manifold (two or more Bogoliubov excitations, or a quantum Monte Carlo treatment): the claim stands if a bound state below the Bogoliubov continuum survives at intermediate and high densities and the spectrum shows the predicted sequence — one narrow line at low density, a broad line plus a sharp below-continuum line with satellites at intermediate density, and only the red-shifted line at high density. If the bound state disappears or the sharp attractive line is absent, the central claim is wrong.","tokens_in":18813,"feed_emoji":"🌀","tokens_out":10920,"duration_ms":93455,"temperature":0.7,"pith_summary":"The paper argues that the backaction of a rotating molecule on a superfluid bath is not a perturbative detail: in the ground state the molecule depletes the condensate and creates a density defect, and this defect is the trap that holds a new class of bound rotor-bath states. Two such states are identified. At intermediate and high densities, the ground state of a fixed-angular-momentum sector is an attractive angulon, a rotor entangled with a Bogoliubov excitation localized inside the defect and lying below the excitation continuum. In dilute environments, the long-lived excited state is instead a repulsive angulon, a quasibound state immersed in the continuum. Rotational spectroscopy of molecules in superfluids should then show a density-driven crossover between the two states, and the effective rotational constant satisfies $B^*_J < B$ at all densities, resolving the earlier anomalous renormalization.","feed_headline":"A molecule binds to the density hole it digs in a superfluid","feed_subtitle":"Rotational spectroscopy maps a density-driven crossover between repulsive and attractive angulon states.","key_machinery":"The load-bearing machinery is a Lee-Low-Pines-like body-fixed frame transformation that decouples the rotor's total angular momentum, combined with a multireference variational ansatz for each sector: $|\\psi_J\\rangle = \\sum_M c_M |J,M\\rangle|f_M\\rangle$, one coherent bath state $|f_M\\rangle$ per rotor projection $M$. Imaginary-time evolution of these parameters self-consistently determines the deformed condensate $|\\Phi_0\\rangle$ with its density defect. On top of that deformed background, the paper builds a Bogoliubov-excitation subspace and diagonalizes an effective Hamiltonian, so the density defect acts as a trapping potential for finite-angular-momentum excitations; when the defect is deep enough, it binds one such excitation around the rotor, producing the attractive angulon below the continuum.","core_discovery":"The central discovery is that the ground state of a sector with nonzero total angular momentum $J$ is not a bare rotor dressed by delocalized bath excitations, but a rotor bound to a localized density wave. The molecule's repulsive interaction with the bosons depletes the condensate over a region of the order of the healing length; once the bath density is high enough that this defect deepens past a scattering resonance, a Bogoliubov excitation with angular momentum $l=1$ localizes inside the defect. The resulting attractive angulon has energy $E_{att} < 2B$ (equivalently $B^*_1 < B$) and a quasiparticle residue $Z_{att}$ that grows with density. In the dilute regime the rotor instead sits inside the Bogoliubov continuum, hybridizing weakly and leaving a sharp, long-lived repulsive angulon close to the bare rotor energy. Real-time rotational spectroscopy computed from the variational equations shows the crossover: a single narrow line at low density, a broadened repulsive line plus a sharp attractive-angulon peak with emission satellites at intermediate density, and only the red-shifted attractive line at high density.","pith_inferences":["If the defect really binds attractive angulons, tuning the bath density or the molecule-boson scattering length should continuously tune the bound-state depth and localization length, making the crossover a dialable spectroscopic feature in ultracold molecular gases.","The same defect-binding mechanism should apply to other angular-momentum-carrying impurities and to anisotropic or microwave-dressed molecular gases, where the rotor-bath interaction can be engineered.","The satellite peaks between the attractive-angulon line and $2B$ encode interactions among multiple Bogoliubov excitations; fitting their spacing could give a direct measure of angulon-cloud interactions, something the paper does not extract.","A natural stress test is whether the bound state survives in a two-excitation or Monte Carlo treatment; if depletion fills the defect, the predicted crossover would shift or disappear."],"forward_implications":["In dilute superfluids, the narrow line seen in rotational spectroscopy is a long-lived repulsive angulon, not the true ground state of the $J=1$ sector.","At intermediate densities the spectrum should show both a broadened repulsive peak and a sharp attractive-angulon line below $2B$, with satellite peaks from sequential emission of Bogoliubov excitations.","At high densities the repulsive angulon disappears and only the red-shifted attractive-angulon line remains.","The effective rotational constant obeys $B^*_J < B$ in all regimes, removing the nonphysical $B^*_1 > B$ found by the single-excitation ansatz.","The attractive angulon is the $J>0$ analogue of an attractive polaron: an impurity bound to a localized bath excitation below the Bogoliubov continuum."],"supporting_citations":[{"why":"Supplies the rotor-bath model and the angulon framework that this paper extends with backaction.","marker":"[48]"},{"why":"Prior LLP treatment of a rotating impurity deforming the bath; the deformed background here becomes the density defect.","marker":"[49]"},{"why":"The prior variational angulon theory whose anomalous $B^*_1 > B$ result this paper identifies as the artifact to resolve.","marker":"[50]"},{"why":"Provides the attractive-polaron concept and dynamics that the attractive angulon is modeled on.","marker":"[42]"},{"why":"Gives the general theory of bound states in boson impurity models used as the analogy for defect-localized angulons.","marker":"[56]"},{"why":"Supplies the variational imaginary- and real-time projection equations used to evolve the multireference state.","marker":"[62]"},{"why":"Supplies the multireference configuration-approach idea behind the variational ansatz.","marker":"[59]"}],"fun_headline_variants":["Molecule binds to its own density hole in superfluid","Angulons switch from repulsive to attractive as density rises","Superfluid rotor binds to its self-dug density crater","Rotational impurity latches onto density defect in superfluid"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The variational space is truncated to one coherent state per rotor projection together with at most one Bogoliubov excitation on the deformed condensate; if multi-excitation correlations or quantum depletion soften or fill the density defect, the claimed attractive angulon bound state may not survive.","fun_headline_variants_meta":{"raw":{"variants":["Molecule binds to its own density hole in superfluid","Angulons switch from repulsive to attractive as density rises","Superfluid rotor binds to its self-dug density crater","Rotational impurity latches onto density defect in superfluid"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000169,"raw_usage":{"total_tokens":1275,"prompt_tokens":966,"completion_tokens":309,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":582,"completion_tokens_details":{"reasoning_tokens":241}},"tokens_in":582,"tokens_out":309,"duration_ms":3347,"temperature":1.0,"reasoning_tokens":241,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:17:18.314684+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the $J=1$ spectral function with a larger variational manifold (two or more Bogoliubov excitations, or a quantum Monte Carlo treatment): the claim stands if a bound state below the Bogoliubov continuum survives at intermediate and high densities and the spectrum shows the predicted sequence — one narrow line at low density, a broad line plus a sharp below-continuum line with satellites at intermediate density, and only the red-shifted line at high density. If the bound state disappears or the sharp attractive line is absent, the central claim is wrong.","supporting_citations":[{"cited_title":"Rotation of quantum impurities in the presence of a many-body environment,","cited_arxiv_id":null,"evidence_quote":"Supplies the rotor-bath model and the angulon framework that this paper extends with backaction."},{"cited_title":"Deformation of a quantum many-particle system by a rotating impurity,","cited_arxiv_id":null,"evidence_quote":"Prior LLP treatment of a rotating impurity deforming the bath; the deformed background here becomes the density defect."},{"cited_title":"Variational theory of angulons and their rotational spectroscopy,","cited_arxiv_id":null,"evidence_quote":"The prior variational angulon theory whose anomalous $B^*_1 > B$ result this paper identifies as the artifact to resolve."},{"cited_title":"Bound states in boson impurity models,","cited_arxiv_id":null,"evidence_quote":"Gives the general theory of bound states in boson impurity models used as the analogy for defect-localized angulons."},{"cited_title":"Variational approach for many-body systems at finite temperature,","cited_arxiv_id":null,"evidence_quote":"Supplies the variational imaginary- and real-time projection equations used to evolve the multireference state."},{"cited_title":"A complete active space scf method (casscf) using a density matrix formulated super-ci approach,","cited_arxiv_id":null,"evidence_quote":"Supplies the multireference configuration-approach idea behind the variational ansatz."}],"review_version":1}