{"id":"0aeb9661-8820-442a-bd9b-6f2bf2121ae3","arxiv_id":"2511.11272","paper_version":4,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"In a p-orbital square-lattice model, chiral phonons induce a static orbital moment at second order in the lattice displacement via coupling to orbital quadrupoles.","lead":"Chiral phonons—atoms moving in circles—are predicted to create a static electronic orbital angular momentum without spin-orbit coupling or magnetism. The work offers a new mechanism for controlling orbital degrees of freedom in simple metals, relevant to orbitronics.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (20) is derived from a long-wavelength Hamiltonian (|k|,|q|<<1/d), then evaluated over the full BZ; for acoustic phonons the q integral is UV-dominated, so the quantitative predictions are uncontrolled.","rationale":"The reader's conditional verdict is appropriate. The strongest claim is Eq. (20), and the weakest point is not the diagonal coupling, which appears benign in the q_y=0 setup because kappa^y=0 and no y-z off-diagonal channel survives at second order; rather, the mismatch between the small-momentum Hamiltonian used to obtain Eq. (20) and the full-BZ sums used to evaluate it is the load-bearing issue. The authors honestly cross-checked the algebra with three methods, and the existence of a second-order chiral-phonon orbital response is plausible (analogous to inverse Faraday and phonon spin magnetization in Ref. [31]). However, the numerical magnitudes and the claimed comparability to photon-driven orbital moments depend on a regime in which the simplified Hamiltonian (10) is not valid. The proposed recomputation with the full matrix elements would settle whether Eq. (20) overestimates, underestimates, or even changes sign relative to the actual tight-binding model. Since the qualitative effect is not in doubt but the quantitative claims are, the verdict should remain CONDITIONAL.","tokens_in":26566,"tokens_out":18641,"duration_ms":176065,"concrete_test":"Re-run the numerical evaluation of Sec. IV replacing the approximate coupling in Eq. (10) with the full k- and q-dependent matrix elements of Eq. (7)/SM Eq. (37), i.e. C_q -> 2i(t_sigma-t_pi)/(dN) kappa^x_{k,q} (and the analogous kappa^y terms), keeping the same tight-binding dispersions and phonon distributions. Repeat for both the coherent phonon case (Eq. 8) and the thermal phonon case (Eq. 24). If the resulting <Lx_0> differs from Eq. (20) by more than ~50% or changes sign, the quoted efficiencies and the '0.4 micro-eV seems detectable' claim are not supported by the paper's derivation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central formula (20) is not a controlled evaluation of the full model. The simplified H_ep in Eq. (10) is obtained in the limit |k|,|q|<<1/d, q_y=q_z<<q_x (main text before Eq. (10); SM Eq. (42)), where kappa^x_{k,q}=sin(k_x d)-sin(k_x+q_x)d is approximated by -q_x d, independent of k. Eq. (20) then sums A_k over the entire Brillouin zone using the full tight-binding dispersions (3)-(5), and the thermal calculation (24) integrates q to the zone boundary. The omitted k-dependence of kappa^x is not small away from Gamma, and the acoustic-phonon integrand |C'_q|^2(n+_q-n-_q) hbar omega_q grows like |q| for small q, so the q sum is dominated by the boundary of the BZ, where the approximation used to derive Eq. (10) is invalid. Consequently, the numerical estimates in Figs. 3-4, the 'efficiency comparable to photons' statement, and the mu_l ~ 0.4 micro-eV estimate are not consequences of the derived formula; they rely on an uncontrolled extrapolation. The effect may survive, but Eq. (20) as applied is a long-wavelength asymptotic expression, not a quantitative prediction.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a tight-binding model of p-orbitals on a square lattice with orbital-dependent electron-phonon coupling. For a chiral phonon mode propagating along x, the authors derive a static orbital accumulation ⟨L_x⟩ at second order in the lattice displacement. The central result is Eq. (20) with A_k defined in Eq. (21), obtained from a Berry-curvature argument and cross-checked by linear response, Keldysh, and equations-of-motion methods. The authors extend the result to nonequilibrium phonons in Eq. (24) and provide numerical estimates, including a claim that the phonon-to-orbital-moment conversion efficiency is comparable to that of circularly polarized light and that the effect is enhanced near orbital degeneracies.","tokens_in":26948,"tokens_out":13379,"duration_ms":131762,"significance":"The qualitative claim — that chiral phonons can produce a static orbital accumulation without spin-orbit coupling, with the sign set by phonon chirality — is conceptually interesting and, if correct, relevant to orbitronics. The derivation is self-contained: it starts from a microscopic tight-binding model and is cross-checked by four independent analytic methods, with no parameter fitted to the target effect. The prediction of enhancement at orbital-degeneracy hot spots is falsifiable. However, the quantitative estimates in Figs. 3-4 and the efficiency comparison to light rest on an uncontrolled extrapolation from a long-wavelength Hamiltonian to the full Brillouin zone. The qualitative mechanism may survive, but the numerical claims are not yet established.","major_comments":[{"comment":"The simplified H_ep in Eq. (10) is derived only in the limit |k|,|q|≪1/d and q_y,q_z≪q_x. In particular, κ^x_{k,q}=sin(k_x d)-sin((k_x+q_x)d) is replaced by -q_x d, dropping a cos(k_x d) factor that is not small over most of the Brillouin zone. Equation (20) then sums A_k over the full BZ, and Eq. (24) integrates q up to the zone boundary. For the acoustic-phonon injection model, |C'_q|^2 (n_+ - n_-) ℏω_q grows with q (linearly in the high-temperature limit), so the q-sum is ultraviolet-dominated and the largest contributions come from the region where the long-wavelength approximation is invalid. The numerical estimates in Figs. 3-4, the μ_l ~ 0.4 μeV value, and the efficiency comparison to photons are therefore not controlled consequences of the derived formula. Please redo the numerical evaluation with the full κ(k,q), or alternatively restrict the calculation to q≪1/d with an explici","section":"Eqs. (10), (20)-(21), (24); SM Eq. (42)"},{"comment":"The statement that the orbital accumulation originates from the entire Fermi sea is inconsistent with Eq. (21). The quantity A_k vanishes identically when f_x=f_y=f_z at the same k, so the nonzero contribution comes only from k points where there is an occupation imbalance between the px, py, and pz orbitals. These are regions where the Fermi level lies between the orbital-resolved bands, not 'all occupied energy eigenstates.' This also affects the comparison in footnote [44] between the present effect and Fermi-surface effects. Please revise the interpretation and make the statement consistent with Eq. (21).","section":"Main text after Eq. (21); Fig. 3 caption"},{"comment":"The diagonal electron-phonon coupling s_αα is dropped with the assertion that it does not contribute to the orbital accumulation. The statement that it 'does not change the orbital degree of freedom' is not by itself a proof, because at second order in the displacement a diagram with one diagonal vertex and one off-diagonal vertex could in principle contribute. Please provide an explicit calculation, or at least a clear symmetry selection-rule argument, showing that such mixed second-order terms vanish. This is load-bearing because if they do not vanish, Eq. (20) is incomplete.","section":"Eq. (7); SM after Eq. (22)"}],"minor_comments":[{"comment":"The relative sign of the u^-_q Q^+_{-q} term differs between Eq. (10) in the main text and the corresponding expression in SM Eq. (42). Since the sign of the chirality response is a central prediction, please reconcile this discrepancy.","section":"Eq. (10) vs. SM Eq. (42)"},{"comment":"The same symbol T is used for both temperature and transmission coefficient in the nonequilibrium phonon distribution. This is confusing and should be changed.","section":"Numerical Estimation"},{"comment":"The vertical axis is labeled in Bohr magnetons per site, but ⟨L_x⟩ as defined is an angular momentum. Please state explicitly how the conversion to μ_B is made (e.g., whether ℏ is set to unity).","section":"Fig. 3 caption"},{"comment":"There are minor typographical issues: 'Kas hiwa' in the affiliation, 'Bohr magnetron' should be 'Bohr magneton,' and 'SA W' in Fig. 2. Please proofread.","section":"Affiliations and general"},{"comment":"The sentence noting that 'the zeroth order terms about ℏω_q are finite, but it vanishes when we assume |q|≪|k|' is cryptic. Please clarify which terms are being referred to and why the q≪k assumption removes them.","section":"SM Keldysh formalism after Eq. (130)"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern about the long-wavelength-to-full-BZ extrapolation is legitimate and lands on the quantitative claims. The qualitative mechanism, however, is plausible and the analytic cross-checks are a real strength. I would not reject the paper: the uncontrolled numerical integration can in principle be fixed by using the full κ(k,q) or by cutting off q and testing convergence. If the numerical estimates cannot be salvaged, the authors should withdraw the quantitative efficiency comparison and reframe the paper around the symmetry-based mechanism. The concurrent work in Ref. [50] is acknowledged; the overlap should be monitored but does not by itself undermine this manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Brief for you: the paper derives a second-order static orbital moment from chiral phonons in a p-orbital square lattice without spin-orbit coupling, and it does so with an unusually thorough toolset — Berry curvature, linear response, Keldysh, and equations of motion all give the same answer. That cross-check is real, and the authors honestly flag that their formula is the orbital analogue of Yao–Murakami's spin result (Ref. 31), plus a concurrent honeycomb paper. For anyone working in orbitronics, the conceptual point — chiral phonons couple to orbital quadrupoles and can rectify into a static orbital dipole — is worth knowing.\n\nWhere it gets soft: the quantitative claims are not on the same footing as the algebra. Eq. (10) is derived assuming |k|,|q| << 1/d and q_y,q_z << q_x, but Eq. (20) is then summed over the whole Brillouin zone, and the thermal phonon sum integrates q to the zone boundary. For acoustic phonons the integrand grows ~q at small q, so the q sum is UV-dominated; the numerical estimates in Figs. 3-4, the \"efficiency comparable to photons\" line, and the 0.4 μeV orbital chemical potential all depend on this uncontrolled extrapolation. The effect may still exist, but those numbers are not consequences of the derived formula. A controlled calculation — keeping the full k-dependence in κ or introducing a proper ultraviolet cutoff or phonon lifetime — is needed before quoting magnitudes.\n\nTwo smaller points. First, the \"entire Fermi sea\" phrasing is wrong in a strict sense: Eq. (21) vanishes when f_x=f_y=f_z, so the effect is not a property of the occupied sea as a whole; it rests on inter-band asymmetries. That is a wording problem, not a fatal one. Second, dropping the diagonal electron-phonon coupling with a one-line assertion is a real gap; if that channel affects the occupations or the energy denominators in second order, Eq. (20) could be incomplete. The authors should prove or at least properly discuss it. And \"0.4 μeV seems detectable\" overstates things.\n\nWorth a serious referee. The mechanism is novel within the subfield, the derivation is internally consistent and multi-method, and the flaws are quantitative overreach and a few unproved shortcuts — all fixable with revision. It is not a desk reject. My recommendation: send to peer review, ask for a controlled evaluation of the q and k sums and a more careful statement of what the formula actually proves.","headline":"A clean multi-method derivation of a plausible but poorly-quantified orbital accumulation from chiral phonons; the numbers should not be trusted, but the mechanism may well survive.","tokens_in":27394,"tokens_out":4750,"would_cite":true,"duration_ms":44463,"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":"Chiral phonons can create a static orbital moment in a nonmagnetic metal without spin-orbit coupling.","keywords":["chiral phonons","orbital accumulation","orbital angular momentum","electron-phonon coupling","orbital quadrupole moment","second-order response","orbitronics","Berry curvature"],"falsifier":"Compute ⟨Lx_0⟩^(2) with the diagonal coupling s_αα retained alongside the off-diagonal channel: any nonzero contribution of order |u|^2 would invalidate Eq. (20); experimentally, reversing the chirality of a surface-acoustic-wave drive in a nonmagnetic, spin-orbit-free film must exactly reverse the sign of the static orbital accumulation.","tokens_in":26439,"feed_emoji":"🌀","tokens_out":7345,"duration_ms":59557,"temperature":0.7,"pith_summary":"The paper claims that chiral phonons—lattice vibrations whose atoms move in circular or helical paths—can leave behind a static, chirality-dependent orbital angular momentum in a conducting electron system, with no magnetism and no spin-orbit coupling required. Working on a square-lattice tight-binding model of p orbitals, the authors show that the orbital-dependent electron-phonon coupling acts on the electrons' orbital quadrupole moments; the linear response oscillates at the phonon frequency and time-averages to zero, but the second-order response rectifies into a uniform orbital accumulation. The central formula, ⟨Lx_0⟩^(2) ≈ Σ_k 4λℏω_q |C_q|^2 U_q^2 A_k, sets the sign by the phonon chirality λ and the magnitude by the occupied band structure. The authors estimate the conversion efficiency per energy flux to be comparable to circularly polarized light—about 10^-17 vs 10^-16 m²/W—and argue the effect grows sharply near orbital degeneracies in the band structure. If correct, this gives a phonon-based route to orbitronic control that does not rely on heavy elements.","feed_headline":"Chiral phonons produce static orbital moments with no spin-orbit coupling","feed_subtitle":"A second-order rectification turns circular lattice motion into a lasting orbital polarization; chirality sets the sign.","key_machinery":"The load-bearing object is the orbital-dependent electron-phonon coupling, which the two-center approximation reduces to H_ep = Σ_q C_q [u+_q Q-_-q + u-_q Q+_-q] in the long-wavelength limit. Here u±_q = uy_q ± i uz_q encode the chiral circular motion of the lattice, Q±_q = Q_xy_q ± i Q_zx_q are orbital quadrupole operators built from the p-orbital angular momentum, and C_q ∝ -(tσ - tπ) q_x/N measures the difference between σ and π hopping. This coupling is what allows chirality to be transferred from the lattice to the orbital sector; the quadrupole operators act as the intermediate that converts the oscillating linear response into a static dipole moment at second order. The same result is","core_discovery":"The central discovery is that a static orbital dipole moment can be induced purely by chiral lattice dynamics, at second order in the lattice displacement, in a system with neither magnetism nor spin-orbit coupling. In their two-center tight-binding model, the electron-phonon coupling between different p orbitals (px, py, pz) is fixed by the difference tσ - tπ of hopping integrals; after a long-wavelength expansion it takes the form H_ep = Σ_q C_q (u+_q Q-_-q + u-_q Q+_-q), where u± are chiral phonon coordinates and Q± = Q_xy ± i Q_zx are orbital quadrupole operators. The linear response of the orbital moment ⟨L_q⟩ oscillates with the phonon phase and vanishes on time average; the second-ord","pith_inferences":["The same rectification should have a reciprocal counterpart: a static orbital polarization in the electrons should couple back to the chiral phonon modes, allowing electrical or optical detection of phonon chirality.","Since the effect scales with ℏω_q and with the Fermi-sea occupation, applying the formula to optical phonons or to band crossings such as Dirac/Weyl points could raise the accumulation well beyond the ~10^-7 Bohr magneton per site estimated here.","A direct test can be made in a nonmagnetic metal driven by a surface acoustic wave: time-resolved X-ray circular dichroism or orbital-sensitive probes should see a static orbital moment that reverses sign when the wave's handedness is reversed, with no accompanying spin signal."],"forward_implications":["A static orbital accumulation can be produced in a nonmagnetic, spin-orbit-free electron system, with its sign controlled simply by the handedness of the chiral phonons.","The energy-flux-normalized conversion efficiency is comparable to that of circularly polarized light (~10^-17 vs ~10^-16 m²/W), so phonon-driven orbitronics can work with ordinary materials rather than heavy elements.","The effect is enhanced near orbital degeneracies—for example, where the pz band touches the px and py bands—and along high-symmetry directions in the Brillouin zone, offering a band-structure design principle.","For thermal chiral phonons, the accumulation is proportional to the imbalance n+_q - n-_q and grows linearly with temperature above the Debye temperature, making a temperature gradient across a chiral material a practical drive."],"fun_headline_variants":["Chiral phonons rectify into static orbital dipoles","Orbital order from chiral phonons without spin-orbit","Static orbital moments via chiral phonon rectification","Chiral lattice motion induces lasting orbital polarization","No spin-orbit needed: chiral phonons make orbital dipoles"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The result leans on the stated-but-unproved assertion that the diagonal electron-phonon coupling s_αα makes no second-order contribution to orbital accumulation, together with the long-wavelength limits |k|,|q| ≪ 1/d and q_y,q_z ≪ q_x.","fun_headline_variants_meta":{"raw":{"variants":["Chiral phonons rectify into static orbital dipoles","Orbital order from chiral phonons without spin-orbit","Static orbital moments via chiral phonon rectification","Chiral lattice motion induces lasting orbital polarization","No spin-orbit needed: chiral phonons make orbital dipoles"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000188,"raw_usage":{"total_tokens":1119,"prompt_tokens":647,"completion_tokens":472,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":391,"completion_tokens_details":{"reasoning_tokens":404}},"tokens_in":391,"tokens_out":472,"duration_ms":4447,"temperature":1.0,"reasoning_tokens":404,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T22:13:31.072826+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute ⟨Lx_0⟩^(2) with the diagonal coupling s_αα retained alongside the off-diagonal channel: any nonzero contribution of order |u|^2 would invalidate Eq. (20); experimentally, reversing the chirality of a surface-acoustic-wave drive in a nonmagnetic, spin-orbit-free film must exactly reverse the sign of the static orbital accumulation.","supporting_citations":[],"review_version":1}