{"id":"54f6768b-1280-4fba-9f29-9d74d7239c6e","arxiv_id":"2507.22571","paper_version":2,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":2,"one_line_summary":"Phonons are predicted to enhance spin polarization and charge transport in β12-borophene nanoribbons, but the supporting self-energy derivation is inconsistent.","lead":"This paper models electron-phonon coupling in borophene nanoribbons and claims that lattice vibrations can strengthen spin-polarized currents in zigzag-edged ribbons. The claim is built on a tight-binding model with a Holstein phonon self-energy, but the inelastic-current derivation in the appendix contains an inconsistent step.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Appendix D's phonon lesser self-energy omits the ±ω0 energy shifts; the claimed spin-polarization enhancement is not supported by the equations as written.","rationale":"The reader's final verdict, REJECT, is well supported, and the decisive issue is the same one identified in the reader's rationale: the phonon self-energy derivation in Appendix D is internally inconsistent. However, the reader's explicitly labeled \"weakest assumption\" was the rigid exchange splitting Δ = 0.20 eV, not the Appendix D error. I regard the Appendix D error as more load-bearing because it invalidates the computed phonon-assisted current, which is the mechanism underlying the paper's central claim. The exact convolution in Eq. D1 produces energy-shifted spectral terms; Eq. D5 collapses them to a single inelastic term at the same energy. Given the sharp DOS features and the sizable optical-phonon coupling reported in the paper, this is not a benign approximation. A corrected derivation might still show phonon-enhanced spin polarization, so the paper's physical idea is not ruled out, but the present manuscript does not establish it. I therefore recommend no change to the reader's verdict; the paper should be revised with a corrected self-energy evaluation before the central claim can be assessed.","tokens_in":27608,"tokens_out":3576,"duration_ms":50582,"concrete_test":"Recompute the phonon-assisted current for the 4-ZZ device at 300 K by evaluating Eq. D1 numerically with the exact delta-function phonon propagator of Eq. D2, instead of using Eq. D5, and then feeding the corrected Σph into Eqs. 11 and 14. If the resulting Iph and spin polarization differ by more than about 20% from Figures 4(b)–4(d), or if the bias onset of spin polarization shifts, the central claim is not supported as written.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim depends on the phonon-assisted current Iph (Eq. 14), which is controlled by the lesser/greater phonon self-energies. Appendix D begins with a convolutional expression, Eq. D1, in which the phonon propagator D<(ω−E) in Eq. D2 enforces energy conservation at E = ω ∓ ω0. However, Eq. D5 replaces that convolution by a single term proportional to fC(ω)AC(ω)(2nB(ω0)+1). The exact evaluation of Eq. D1 gives [nB(ω0)fC(ω−ω0)AC(ω−ω0) + (1+nB(ω0))fC(ω+ω0)AC(ω+ω0)] times the coupling sum; the two terms coincide with Eq. D5 only if fC and AC are constant on the scale of the phonon energy. The paper's own DOS figures show features with widths of order 50–70 meV, and the dominant optical phonon mode contributes a large coupling, so the energy shift is not negligible. Because the claimed spin-polarization enhancement arises precisely from phonon-assisted access to sharp, spin-split conduction states, dropping the ±ω0 shifts can create or eliminate the effect. The central claim is therefore unsupported by the derivation as written. A second concern, the rigid exchange splitting Δ = 0.20 eV, also deserves scrutiny, but the Appendix D error alone blocks evaluation of the paper's main conclusion.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a multiscale study of electron-phonon coupling (EPC) in β12-borophene nanoribbons (BNRs). The authors combine DFT-derived tight-binding parameters and EPC matrix elements with a Holstein-model self-energy and a nonequilibrium Green's function/Landauer-Büttiker transport calculation. Their central claim is that phonon-assisted inelastic processes substantially modify both charge and spin transport: in zigzag-edged BNRs, EPC enhances the weak intrinsic spin polarization at lower bias than in the ballistic case, while in nonmagnetic edge configurations, EPC increases the total charge current. The conclusions are based on spin-resolved current-voltage characteristics, differential conductance, and density-of-states comparisons for several ribbon widths and lengths.","tokens_in":27796,"tokens_out":7441,"duration_ms":98450,"significance":"If the central claim were established, the paper would identify lattice vibrations as a tunable knob for spin filtering in borophene-based nanoribbons, which is a potentially interesting result for 2D spintronics. The work also has genuine strengths: the EPC constants are obtained from DFT rather than fitted to the transport target, and the paper attempts a transparent multiscale workflow. However, the main quantitative conclusion rests on a phonon self-energy derivation in Appendix D that is not correct as written, and the spin-polarization mechanism is built on a rigid exchange-splitting assumption that is not validated at the device level. The significance of the paper therefore cannot be assessed until these load-bearing points are resolved.","major_comments":[{"comment":"Evaluating the convolution in Eq. (D1) with D< from Eq. (D2) and g< from Eq. (D3) gives, up to the common prefactor i/(2N)∑g², the expression n_B(ω0) f_C(ω−ω0) A_C(ω−ω0) + [1+n_B(ω0)] f_C(ω+ω0) A_C(ω+ω0). Equation (D5) instead states f_C(ω) A_C(ω)[2 n_B(ω0)+1], which is valid only if the spectral function and Fermi function are approximately constant on the scale of the phonon energy. The paper's own DOS features (Sec. III.A, Fig. 3) have widths of order 50–70 meV, comparable to the phonon energies relevant here, and the claimed mechanism relies on phonon-assisted access to sharp spin-split states. Since Eq. (14) builds Iph directly from Σ<,>ph, the central quantitative claim that EPC enhances spin polarization is not supported by the equations as written.","section":"Appendix D, Eqs. (D1)–(D6)"},{"comment":"The spin polarization in the transport model is introduced as a fixed, rigid exchange splitting Δ = 0.20 eV taken from prior DFT work and applied uniformly to the nonmagnetic Hamiltonian for all ribbon widths, lengths, biases, and temperatures. The claimed phonon enhancement of spin polarization is therefore conditional on this rigid-band approximation being quantitatively valid for the finite devices studied. This is an assumption that can be tested: a comparison with spin-resolved DFT density of states for the relevant ribbon widths, or a self-consistent treatment of the exchange splitting under bias, would show whether the splitting remains rigid and of the assumed magnitude. I do not regard the use of Δ as circular, since the EPC parameters are computed independently, but the paper should provide such a test before claiming that EPC enhances intrinsic polarization.","section":"Sec. III.B and Supplemental Eq. (S1)"},{"comment":"It is not clear which phonon self-energy is actually used in the numerical calculations. The main text presents the retarded self-energy in Eq. (10) with the expected ±ωqλ energy shifts, while Appendix D derives lesser and greater self-energies in Eqs. (D5)–(D6) that omit these shifts. If the latter are analytically continued and used in Eq. (11), the peak-shift error propagates into the renormalized DOS and the coherent current I0 as well as into Iph. If, instead, the exact Eq. (10) is used for the retarded self-energy, then the manuscript uses two inconsistent phonon self-energies, and the origin of the numerical results is not reproducible. The authors should specify the exact self-energies used and provide a corrected derivation.","section":"Sec. II.C, Sec. II.D, and Appendix D"}],"minor_comments":[{"comment":"The tight-binding block notation is inconsistent: the main text and the first version of Appendix B use α, β, βw, βl, while the duplicated version of the appendix uses ε, ω, ωw, ωl for the same objects. This makes it unnecessarily difficult to reconstruct the device Hamiltonian.","section":"Appendix B and main text"},{"comment":"Several paragraphs and figure captions appear twice with different wording and corrupted symbols (e.g., the description of Fig. 3 near the start of Sec. III.A and the repeated Fig. 3/Fig. 4 captions). These artifacts must be removed and replaced by a single consistent text.","section":"General presentation"},{"comment":"The text refers to 'Fig. 3(e)' when discussing the decomposition of the differential conductance into coherent and phonon-assisted components, but the corresponding panel appears to be Fig. 4(e). Please correct cross-references between figures and text.","section":"Sec. III.B, Fig. 4"},{"comment":"The transfer-matrix expansion in Eq. (C2) uses quantities t0, ¯t0, t1, etc., but the recursion relations in Eq. (C3) are written for tn and ¯tn without explicitly defining the relationship between these sequences; adding a short definition would make the algorithm reproducible.","section":"Appendix C, Eqs. (C2)–(C3)"}],"recommendation":"reject","confidential_remarks":"The Appendix D error is central: the phonon-assisted current Iph, which carries the claimed spin-polarization enhancement, is computed from lesser/greater self-energies that omit the ±ω0 energy conservation shifts. Correcting this requires recomputing essentially all transport results, so the present manuscript cannot be accepted as is; a substantially revised version with corrected derivation and new numerics would be needed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the mode-resolved electron-phonon coupling for β12-borophene nanoribbons computed from DFT, and the spin-resolved transport with and without Holstein phonons for different edge terminations. Those numbers are new, and the qualitative trends — phonons broaden the DOS, smooth the conductance steps, and increase the current — are plausible and consistent with the machinery. There is no sign of curve-fitting; the EPC inputs come from independent DFT calculations, which is to the authors' credit.\n\nThe trouble is that the central quantitative claim rests on a self-energy that is not actually derived correctly. In Appendix D, Eq. (D1) defines a convolution of the electronic Green's function with the phonon propagator. Eq. (D5) then evaluates that convolution by pulling f_C(ω) and A_C(ω) out of the integral. The exact evaluation of (D1) with Eq. (D2) gives the Bose-weighted combination [n_B(ω₀) f_C(ω−ω₀) A_C(ω−ω₀) + (1+n_B(ω₀)) f_C(ω+ω₀) A_C(ω+ω₀)]. Eq. (D5) replaces this with f_C(ω) A_C(ω) (2n_B+1), which is only valid if the spectral function is flat on the scale of the phonon energy. The paper's own Figure 3 shows peaks with widths of tens of meV, and the dominant phonon modes are optical modes with substantial ω₀. This is not a cosmetic slip: the spin-polarization enhancement is claimed to arise precisely because phonon-assisted processes access sharp spin-split states, and shifting those states by ±ω₀ can create or destroy the effect. As the equations stand, the main claim is not supported.\n\nBeyond that, the rigid exchange splitting Δ = 0.20 eV taken from bulk DFT and applied uniformly to every ribbon width, length, and bias is a justified simplification only up to a point; in nanoribbons with edge magnetism the splitting can be width- and position-dependent. The manuscript also has presentational problems: duplicated passages, inconsistent notation (ω12 vs. β12), a dangling \"Figure??\", and a placeholder SM URL. None of this is fatal by itself, but it compounds the review burden.\n\nWho should read this? People working on borophene spintronics and phonon-assisted transport will want to know the mode-resolved EPC numbers and the general idea that lattice vibrations can tune spin transport in these ribbons. As submitted, though, the paper does not establish that idea quantitatively.\n\nI would send this to peer review rather than desk reject, because the error is specific, checkable, and fixable, and the qualitative physics is not obviously wrong. But a serious referee should expect major revision: the self-energy must be corrected and the transport recalculated. If the effect survives, it becomes a solid paper; if not, the phonon-tuning claim goes away. Either way, the current version is not publishable.","headline":"The phonon-assisted transport derivation in Appendix D drops the ±ω₀ energy shifts, and the paper's own DOS makes that error load-bearing; the manuscript deserves a demanding referee, not a desk reject.","tokens_in":28426,"tokens_out":3363,"would_cite":false,"duration_ms":38501,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["73.63.-b","72.10.Di","71.38.-k"],"model":"deepseek-v4-flash","headline":"This paper argues that lattice vibrations can tune both charge and spin transport in β12-borophene nanoribbons: phonon-assisted processes strengthen the weak spin polarization of zigzag edges and enhance current in nonmagnetic edge…","keywords":["β12-borophene nanoribbons","electron-phonon coupling","Holstein model","spin polarization","phonon-assisted transport","Landauer-Büttiker formalism","zigzag edge states","nonmagnetic edge configurations"],"falsifier":"Measure the current–voltage curve of a zigzag β12-borophene nanoribbon at low temperature and again at 300 K: if spin polarization appears at the same bias and the conductance peaks near 1.1 V and 1.7 V do not grow with temperature, the phonon-assisted mechanism is not what controls the transport. A complementary calculation would recompute the spin splitting self-consistently for each ribbon width and show whether it stays at 0.20 eV.","tokens_in":27311,"feed_emoji":"🧲","tokens_out":8113,"duration_ms":87403,"temperature":0.7,"pith_summary":"The paper claims that electron-phonon coupling (EPC) can be used as a tuning knob for transport in β12-borophene nanoribbons. In zigzag-edged ribbons, where intrinsic spin polarization is weak, lattice vibrations strengthen the spin polarization and make it appear at lower bias, while in nonmagnetic edge configurations they boost the charge current. The effect shows up as marked changes in current–voltage characteristics, including conductance peaks at biases near 1.1 V and 1.7 V, so it is open to direct experimental probing. The authors build a tight-binding model of the boron pz orbitals, add local Holstein phonons with coupling constants computed from first principles, and solve the resulting device currents with the Landauer–Büttiker Green's function formalism.","feed_headline":"Phonons amplify spin polarization in borophene nanoribbons","feed_subtitle":"Thermal lattice motion opens new conduction channels, so zigzag ribbons polarize spin at lower voltage.","key_machinery":"The central object is the Holstein-model phonon self-energy $\\Sigma_{\\mathrm{ph}}$ added to the retarded Green's function of the central scattering region. It produces a total current $I = I_0 + I_{\\mathrm{ph}}$, where $I_0$ is the coherent Landauer–Büttiker current and $I_{\\mathrm{ph}}$ is the phonon-assisted inelastic current carried by electrons that emit or absorb a phonon. The coupling constants $g_{q\\lambda}$ that feed this self-energy are Brillouin-zone-averaged first-principles electron-phonon matrix elements, and spin enters through a rigid exchange splitting applied to the tight-binding Hamiltonian. This machinery converts lattice motion into bias-dependent extra conduction channels and predicts the voltages at which those channels switch on.","core_discovery":"The central discovery is that phonon-assisted inelastic processes open additional conduction channels in β12-borophene nanoribbons, and these channels can favor one spin direction. In zigzag-edged ribbons, where edge magnetism is weak and modeled as a rigid splitting Δ = 0.20 eV, adding EPC raises the current in both spin channels but selectively strengthens the spin-up channel in wider ribbons, enough to reverse the spin-conductance order found without phonons and to make spin polarization appear at lower bias. The extra current is dominated by the inelastic component $I_{\\mathrm{ph}}$, while the coherent part barely changes. In nonmagnetic line-line, AC-AC, and AC-line edges, EPC smooths quantized conductance steps into a more ohmic-like response and raises the current; for line-line ribbons it reverses the length trend, so the longest ribbon carries the most current. The microscopic source is traced to two phonon modes—a long-wavelength acoustic mode and a bond-stretching optical mode—that together contribute more than 90 percent of the total coupling strength of about 0.55 eV.","pith_inferences":["If the rigid-exchange-splitting model is replaced by a self-consistent magnetic treatment, the qualitative phonon-induced enhancement is likely to survive, but the bias at which spin polarization appears and the size of the I–V features could shift; this is a natural next calculation.","A direct experimental test would be to measure the temperature dependence of the differential conductance: the predicted resonances near 1.1 V and 1.7 V should grow as the phonon population increases, and should weaken toward the ballistic low-temperature limit.","The same Holstein-plus-NEGF machinery could be applied to other zigzag-edged 2D ribbons with weak edge magnetism, where phonon-assisted channels might similarly alter the dominant spin channel.","Because the two dominant phonon modes are identified, isotope substitution or selective strain could shift their frequencies and thereby move the conductance resonances, effectively phonon-engineering the spin filter."],"forward_implications":["A zigzag β12-borophene nanoribbon device operated at finite temperature should show spin polarization at lower bias than its ballistic prediction, with the spin-down channel carrying more current at low bias and the spin-up channel catching up in wider ribbons.","Phonon-assisted transport lifts the exponential length suppression of coherent tunneling, so in nonmagnetic line-line ribbons longer channels can carry more current than shorter ones.","Differential-conductance measurements should reveal phonon-specific resonances near 1.1 V and 1.7 V in zigzag ribbons, effectively turning transport measurements into phonon spectroscopy.","Edge configuration becomes a design variable: line-line edges maximize phonon-enhanced charge current, while zigzag edges give phonon-enhanced spin polarization, so a single material can be patterned for either function.","Because the inelastic component dominates at 300 K, moderate heating acts as a controllable conductance boost rather than merely a scattering loss."],"supporting_citations":[{"why":"Supplies the prior prediction that β12-borophene has medium-strength EPC (dimensionless coupling near 0.89), motivating the inclusion of phonons.","marker":"[33]"},{"why":"Provides DFT-derived tight-binding parameters and the rigid exchange splitting Δ = 0.20 eV used for the magnetic zigzag ribbon.","marker":"[22]"},{"why":"Provides the Dirac-fermion characterization of β12-borophene and the exchange-splitting value the spin model is built on.","marker":"[24]"},{"why":"Defines the Holstein molecular-crystal model used for local electron-phonon coupling.","marker":"[37]"},{"why":"Defines the small-polaron Holstein model underlying the phonon self-energy treatment.","marker":"[38]"},{"why":"Supplies the first-principles method for computing electron-phonon coupling matrix elements used to obtain $g_{q\\lambda}$.","marker":"[44]"},{"why":"Supplies the complementary first-principles formalism for phonon-limited transport used for the coupling matrix elements.","marker":"[45]"},{"why":"Gives the Matsubara-summed Holstein self-energy expression used to renormalize the electronic Green's function.","marker":"[47]"},{"why":"Supplies the Landauer–Büttiker formula and NEGF transport framework for the coherent current $I_0$.","marker":"[35]"},{"why":"Supplies the iterative surface Green's function scheme used to build the semi-infinite lead self-energies.","marker":"[58]"}],"fun_headline_variants":["Phonons turn up spin in borophene nanoribbons","Lattice vibrations sharpen spin polarization in borophene","Thermal motion boosts spin current in borophene ribbons","Phonon-assisted spin filtering in borophene nanoribbons","Vibrations enhance spin transport in borophene edges"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculations assume the magnetism of every zigzag ribbon can be captured by one fixed rigid energy splitting between the two spin directions, 0.20 eV, no matter how wide, long, biased, or hot the ribbon is.","fun_headline_variants_meta":{"raw":{"variants":["Phonons turn up spin in borophene nanoribbons","Lattice vibrations sharpen spin polarization in borophene","Thermal motion boosts spin current in borophene ribbons","Phonon-assisted spin filtering in borophene nanoribbons","Vibrations enhance spin transport in borophene edges"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000931,"raw_usage":{"total_tokens":4020,"prompt_tokens":1014,"completion_tokens":3006,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":630,"completion_tokens_details":{"reasoning_tokens":2920}},"tokens_in":630,"tokens_out":3006,"duration_ms":23983,"temperature":1.0,"reasoning_tokens":2920,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T11:32:05.656296+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the current–voltage curve of a zigzag β12-borophene nanoribbon at low temperature and again at 300 K: if spin polarization appears at the same bias and the conductance peaks near 1.1 V and 1.7 V do not grow with temperature, the phonon-assisted mechanism is not what controls the transport. A complementary calculation would recompute the spin splitting self-consistently for each ribbon width and show whether it stays at 0.20 eV.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the prior prediction that β12-borophene has medium-strength EPC (dimensionless coupling near 0.89), motivating the inclusion of phonons."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides DFT-derived tight-binding parameters and the rigid exchange splitting Δ = 0.20 eV used for the magnetic zigzag ribbon."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Dirac-fermion characterization of β12-borophene and the exchange-splitting value the spin model is built on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Holstein molecular-crystal model used for local electron-phonon coupling."},{"cited_title":"Gao, Q.-Z","cited_arxiv_id":null,"evidence_quote":"Defines the small-polaron Holstein model underlying the phonon self-energy treatment."},{"cited_title":"Gunst, T","cited_arxiv_id":null,"evidence_quote":"Supplies the first-principles method for computing electron-phonon coupling matrix elements used to obtain $g_{q\\lambda}$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the complementary first-principles formalism for phonon-limited transport used for the coupling matrix elements."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the Matsubara-summed Holstein self-energy expression used to renormalize the electronic Green's function."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Landauer–Büttiker formula and NEGF transport framework for the coherent current $I_0$."},{"cited_title":"Chakrabarty, A","cited_arxiv_id":null,"evidence_quote":"Supplies the iterative surface Green's function scheme used to build the semi-infinite lead self-energies."}],"review_version":1}