{"id":"3f768c01-9405-4233-96ee-650b46844b78","arxiv_id":"2508.14270","paper_version":1,"verdict":"CONDITIONAL","confidence":"LOW","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A real-space first-principles method for computing orbital Hall transport and orbital accumulation in metallic multilayers is developed, predicting substantial orbital accumulation even in centrosymmetric systems.","lead":"A new computational method combines real-space density functional theory with a Chebyshev expansion to compute orbital and spin transport in metal multilayers directly in real space. It targets orbitronics, an emerging field that aims to use the electrons' orbital motion for information processing.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quantitative orbital-transport results rest on an undefined orbital current operator and an intra-atomic L approximation; without a bulk benchmark against established OHE calculations, the central numbers are unsecured.","rationale":"The reader identified the intra-atomic definition of the orbital angular momentum operator and the missing explicit orbital current operator as the weakest assumption. I agree: the central claim's quantitative content is entirely determined by this choice, and the supplied manuscript text does not state how the position operator, the velocity operator, and the orbital current operator are constructed in the real-space LMTO-ASA basis. Ref. [53] is cited, showing awareness of the anomalous-position issue, but the visible methodology does not explain whether or how that contribution is included. The paper has real supporting infrastructure—the RS-LMTO-ASA method is established, and the KPM/Kubo-Bastin route is well referenced—so the concern is not about the framework as a whole but about the unvalidated operator representation. A single bulk benchmark would settle this. Since the reader's verdict was already CONDITIONAL with low confidence, my analysis does not move the verdict; it sharpens the condition that must be met before acceptance.","tokens_in":9990,"tokens_out":5945,"duration_ms":81272,"concrete_test":"Compute the orbital Hall conductivity sigma_xy^L for a bulk centrosymmetric transition metal used in the paper (e.g., hcp Ti or fcc Pt/W) with the same RS-LMTO-ASA real-space Kubo-Bastin implementation and the actual orbital-current operator employed in the code. Compare the value at the same chemical potential and broadening with converged Wannier-interpolation results from Salemi & Oppeneer (PRM 6, 095001) and Go et al. (PRB 109, 174435). If sigma_xy^L differs by more than ~25% or changes sign, the intra-atomic/anomalous-position treatment is quantitatively inadequate and the multilayer accumulation predictions must be revised; if it matches, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that the RS-LMTO-ASA real-space method computes orbital Hall transport and orbital accumulation in transition-metal multilayers. In the visible text, the orbital angular momentum and orbital current operators are never explicitly defined. Only the Hamiltonian, Eq. (1), and the Chebyshev LDOS expansion, Eq. (2), are shown. If the implementation uses an intra-atomic, site-diagonal orbital angular momentum operator and neglects the inter-atomic 'anomalous position' contributions highlighted in Ref. [53], then the orbital Hall conductivity and the resulting accumulation values can be quantitatively wrong, potentially even in sign, especially for 5d metals where orbital texture is strongly nonlocal. This is not a proof of error, but it is the weakest load-bearing link: the new physical predictions—substantial orbital/spin accumulation in these heterostructures—depend directly on this operator choice, and no benchmark against established Wannier/Kubo orbital Hall conductivity results is shown in the supplied text. The absence of the explicit current operator and of any validation against known OHE values makes the central numerical claims conditional.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a real-space DFT-based method for orbitronic transport in metallic multilayers. It combines the RS-LMTO-ASA Hamiltonian with a Chebyshev polynomial expansion of the Green's functions and evaluates Kubo-Bastin linear response directly in real space. The method is applied to Co/Ti, Fe/Ti, Ni/Ti, Co/W, Fe/W, and Ni/W heterostructures, and the central claim is that substantial orbital (and spin) accumulation can appear even in centrosymmetric systems, driven by band-structure asymmetries and interfacial scattering. The authors also claim linear scaling with system size and the natural inclusion of disorder, finite-size effects, and interface roughness.","tokens_in":10262,"tokens_out":2508,"duration_ms":31163,"significance":"If fully validated, this approach would provide a valuable real-space alternative to Wannier/Kubo methods for orbital transport in large, inhomogeneous heterostructures, with no fitted parameters entering the transport calculation. The methodological skeleton is credible: Kubo-Bastin plus KPM is a well-established combination, and the RS-LMTO-ASA Hamiltonian is an independent first-principles construction. However, the key physical observables--the orbital current operator, the orbital accumulation definition, and the orbital angular momentum operator in the chosen basis--are never explicitly given, and no benchmark against established bulk orbital Hall conductivity values is shown. Without these, the quantitative results, including the central 'centrosymmetric accumulation' claim, remain unsecured.","major_comments":[{"comment":"The orbital Hall conductivity and orbital accumulation are the central outputs, but the manuscript never defines the orbital angular momentum operator L, the orbital current operator J_orb, or the accumulation operator used in the Kubo-Bastin trace. Eq. (1) gives the Hamiltonian and Eq. (2) gives the LDOS expansion, but no formula such as J_orb^γ = {L^α, v^β}_+/2 or the corresponding Kubo-Bastin expression is displayed. The numerical results in Figs. 4-5 depend directly on this operator choice. Please provide these definitions explicitly.","section":"Sec. II"},{"comment":"No validation against known orbital Hall conductivity values for the same elemental materials is presented. For example, bulk Ti, W, and Pt have published first-principles OHE results from Wannier/Kubo methods; a single comparison for one of these would anchor the method and calibrate the magnitude and sign of the results. As written, the quantitative values of the orbital conductivities and accumulations rest on an unverified implementation of the current operator.","section":"Sec. III / bulk validation"},{"comment":"The manuscript cites Ref. [53] on the role of inter-atomic 'anomalous position' contributions to the orbital Hall effect, but does not state how that issue is handled in the present real-space framework. If the implementation uses an intra-atomic, site-diagonal L and neglects inter-atomic contributions, the resulting OHE values can be quantitatively inaccurate, especially for 5d metals where orbital texture is nonlocal. The authors should either justify this approximation for their specific systems or provide a concrete test of its magnitude.","section":"Sec. II, Ref. [53]"},{"comment":"The orbital accumulation results are reported as layer-resolved in-plane components L_x and L_y, but no definition of 'orbital accumulation' is given. It is not stated whether this is a non-equilibrium expectation value of an intra-atomic L operator, how it is projected onto atomic layers, or what units are used (e.g., ħ per atom). This makes the central quantitative claim about substantial accumulation in centrosymmetric systems impossible to interpret or reproduce.","section":"Sec. III, Fig. 5"}],"minor_comments":[{"comment":"Several mathematical symbols appear as placeholder glyphs in the displayed equations (e.g., the Hamiltonian and LDOS formulas). The typesetting needs to be fixed so that the equations are readable.","section":"Eqs. (1)-(3)"},{"comment":"The abstract states that the method 'naturally incorporates disorder, finite-size effects, and interface roughness,' but the visible results and discussions do not include any disordered system or explicit interface roughness calculation. Either add such a calculation or temper the claim.","section":"Abstract / Sec. I"},{"comment":"The Kubo-Bastin formula is cited as [54-60], but the actual linear-response expression used for the conductivity and accumulation is not written out. Including the explicit Kubo-Bastin formula would also clarify how the Chebyshev moments enter.","section":"Sec. II"},{"comment":"The convergence of the Chebyshev expansion is reported only as a fixed 'N = 500'. A brief convergence test with respect to N would strengthen confidence in the numerical results.","section":"Sec. II"},{"comment":"There is a typographical error: 'layered TM systems an a series of TM-based heterostructures' should read 'and a series.'","section":"Introduction"}],"recommendation":"major_revision","confidential_remarks":"The methodological direction is promising and I do not see circularity: the transport calculation uses the self-consistent DFT Hamiltonian with no fitted parameters. The main risk is that the missing operator definitions and lack of benchmark leave the quantitative claims unverified. The authors should be asked to provide the explicit current/accumulation operators and at least one bulk OHE benchmark. I would also encourage them to consider making the code or representative data available, since the method paper's value depends on reproducibility."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a new and non-routine methodological step—RS-LMTO-ASA real-space DFT with Chebyshev-based Kubo-Bastin for orbital and spin Hall transport—and it deserves a serious referee. But the central quantitative claim, substantial orbital accumulation in centrosymmetric multilayers, is not yet supported by the visible evidence.\n\nWhat's new: the integration itself. The RS-LMTO-ASA Hamiltonian is used directly in a Kubo-Bastin transport calculation, so you get linear scaling, open boundaries, and natural treatment of disorder and interfaces. That's a real addition to the orbitronics toolbox, which currently mostly relies on Wannier or tight-binding fits from momentum-space DFT. The paper is honest about the ingredients being established; the new part is putting them together.\n\nWhat's good: the technical foundation is sound. The KPM recursion for the LDOS is standard, and using it for self-consistency is a sensible extension. There are no fitted parameters in the transport calculation, and the paper cites the relevant literature, including Go et al. on the anomalous position.\n\nSoft spots, in proportion: the visible text never defines the orbital current operator or the orbital angular momentum operator. This is not cosmetic. The orbital Hall conductivity in 5d metals is known to be sensitive to whether you use intra-atomic L or include the inter-atomic anomalous position terms. Since the paper's headline result—accumulation in centrosymmetric systems—depends directly on this operator choice, the absence of a bulk benchmark against established OHE values (e.g., from Wannier calculations) is a real gap. There is also no convergence study for the N=500 Chebyshev truncation or system size. These are all fixable in a revision, but they are not present in this version.\n\nOne thing the reader's report got right: this is not circular. RS-LMTO-ASA is independently validated and the transport calculation is parameter-free.\n\nWho this is for: anyone working on orbitronics or real-space transport methods. The paper would benefit from a referee who knows the orbital-current literature, because the key issue is exactly the operator definition.\n\nRecommendation: send it to peer review. It deserves referee time. Ask the authors to (1) state the orbital current and L operator definitions explicitly, (2) validate against known OHE values for bulk Ti, W, etc., and (3) provide convergence data. With those, it could become a solid methods paper.","headline":"A genuinely new real-space method for orbital/spin Hall transport, but the central numbers are unsecured until the orbital current operator is defined and benchmarked.","tokens_in":10737,"tokens_out":3231,"would_cite":true,"duration_ms":35455,"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 real-space DFT method computes orbital and spin Hall transport and accumulation directly in transition-metal multilayers, showing substantial orbital accumulation even in centrosymmetric systems.","keywords":["orbitronics","orbital Hall effect","orbital accumulation","real-space DFT","Kubo-Bastin formula","Chebyshev polynomial expansion","transition-metal multilayers","spin Hall effect"],"falsifier":"Compare the real-space intra-atomic orbital Hall conductivities for bulk transition metals such as Ti, Cr, or W against converged momentum-space DFT calculations that include the orbital dependence of the anomalous position; a mismatch beyond numerical precision would show the omitted inter-atomic contributions are not negligible. A second test is to measure layer-resolved orbital accumulation in a centrosymmetric Co/Ti or Ni/W stack with magneto-optical detection: the paper's central prediction is a clear in-plane orbital/spin signal, so a null result would refute it.","tokens_in":9929,"feed_emoji":"🧲","tokens_out":7824,"duration_ms":81353,"temperature":0.7,"pith_summary":"This paper develops a first-principles method for computing orbitronic transport in complex metallic heterostructures without ever leaving real space. It combines the RS-LMTO-ASA Hamiltonian with a Chebyshev polynomial expansion of the Green's function, so that orbital Hall conductivities and orbital/spin accumulations are evaluated directly from the real-space DFT Hamiltonian via the Kubo-Bastin formula. The method scales linearly with the number of nonequivalent atoms and can incorporate disorder, interface roughness, and open boundaries without extra approximations. Applied to transition-metal multilayers such as Co/Ti, Fe/Ti, Ni/Ti, Co/W, Fe/W, and Ni/W, it predicts substantial orbital (spin) accumulation even in centrosymmetric systems, driven by band-structure asymmetries and interfacial scattering. If correct, it provides a scalable way to simulate and design orbital-torque devices in realistic, disordered heterostructures.","feed_headline":"Even symmetric metal stacks build up orbital currents","feed_subtitle":"Real-space DFT plus Chebyshev expansions compute orbital Hall transport in realistic, disordered multilayers.","key_machinery":"The load-bearing objects are (i) the RS-LMTO-ASA real-space Hamiltonian in its orthogonal tight-binding representation (Eq. 1), built from potential parameters and structure constants, and (ii) the Chebyshev polynomial expansion of the Green's function used to evaluate the Kubo-Bastin formula. The Chebyshev recursion constructs moments and reconstructs LDOS and response functions with numerical stability and linear scaling; each nonequivalent atom needs one recursion. This replaces momentum-space DFT followed by Wannier or pseudo-atomic projection, keeping disorder, interfaces, and open boundaries as explicit real-space features. The orbital angular momentum current is defined within an intr","core_discovery":"The central claim is that a fully real-space first-principles calculation can capture orbital Hall transport and non-equilibrium orbital/spin accumulation in layered transition-metal systems. The RS-LMTO-ASA orthogonal Hamiltonian is kept in real space, and the Kubo-Bastin formula is evaluated with a Chebyshev expansion of the Green's functions. Layer-resolved results for ferromagnet/normal-metal stacks show substantial in-plane orbital and spin accumulation even in centrosymmetric geometries, driven by band-structure asymmetry and interfacial scattering. The paper further argues this accumulation can be steered by structural confinement to engineer orbital torque in transition-metal bilayer","pith_inferences":["If the centrosymmetric accumulation result survives inclusion of inter-atomic anomalous-position terms, it would open a design route for orbital-memory and torque devices based on symmetric stacks, which are easier to grow; the paper itself only notes the possibility of torque engineering.","The same Chebyshev/Kubo-Bastin real-space pipeline could be applied to other linear-response functions, such as the orbital Edelstein effect, orbital Nernst response, or frequency-dependent conductivities, without changing the Hamiltonian, though the paper does not compute these.","A direct numerical comparison of the intra-atomic orbital-current definition with calculations including the orbital dependence of the anomalous position would map the regime where the real-space method is quantitatively reliable.","Varying interface roughness or alloy disorder systematically in these multilayers and comparing layer-resolved accumulation would help separate band-structure-driven from scattering-driven contributions."],"forward_implications":["Orbital and spin Hall conductivities and accumulations can be computed for large, disordered, finite multilayers without Wannier projection, directly from the DFT Hamiltonian.","Because the cost scales linearly with the number of nonequivalent atoms, full heterostructures with interface roughness or compositional disorder become tractable at first-principles level.","Centrosymmetric metallic stacks are predicted to host substantial orbital and spin accumulation, so inversion symmetry alone does not suppress orbitronic effects.","The interplay between orbital Hall generation and structural confinement in transition-metal bilayers is a handle for engineering orbital torque.","The same Green's-function machinery yields layer-resolved profiles, tying transport quantities to individual atomic layers in the stack."],"supporting_citations":[{"why":"Supplies the real-space DFT electronic-structure method (RS-LMTO-ASA) whose Hamiltonian is the starting point for transport.","marker":"[49–52]"},{"why":"Provides the underlying LMTO-ASA linear band-theory formalism on which the real-space Hamiltonian is built.","marker":"[72]"},{"why":"Establishes the real-space Chebyshev approach to the Kubo conductivity that the paper adapts for orbital and spin Hall currents.","marker":"[54]"},{"why":"Provides the kernel polynomial method used to regularize and truncate the Chebyshev expansion.","marker":"[85]"},{"why":"Defines the Wannier-interpolation route and the orbital-dependent anomalous position, the comparison point for the paper's direct real-space intra-atomic approach.","marker":"[53]"},{"why":"Gives quantitative references for spin and orbital polarization in heavy-metal/3d-metal bilayers for comparison with the accumulation results.","marker":"[86]"}],"fun_headline_variants":["Symmetric metal stacks still accumulate orbital currents","Orbital currents arise even in symmetric metal stacks","Centrosymmetric stacks yield orbital accumulation","Symmetry can't stop orbital currents in metal layers","Real-space DFT finds orbital currents in symmetric stacks"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The calculation assumes each atom's orbital angular momentum lives entirely inside that atom's sphere, ignoring the inter-atomic pieces of the orbital motion that some analyses show contribute to the orbital Hall effect; if those inter-atomic pieces are significant, the predicted numbers change.","fun_headline_variants_meta":{"raw":{"variants":["Symmetric metal stacks still accumulate orbital currents","Orbital currents arise even in symmetric metal stacks","Centrosymmetric stacks yield orbital accumulation","Symmetry can't stop orbital currents in metal layers","Real-space DFT finds orbital currents in symmetric stacks"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001376,"raw_usage":{"total_tokens":5354,"prompt_tokens":631,"completion_tokens":4723,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":375,"completion_tokens_details":{"reasoning_tokens":4667}},"tokens_in":375,"tokens_out":4723,"duration_ms":35552,"temperature":1.0,"reasoning_tokens":4667,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T18:39:54.925780+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compare the real-space intra-atomic orbital Hall conductivities for bulk transition metals such as Ti, Cr, or W against converged momentum-space DFT calculations that include the orbital dependence of the anomalous position; a mismatch beyond numerical precision would show the omitted inter-atomic contributions are not negligible. A second test is to measure layer-resolved orbital accumulation in a centrosymmetric Co/Ti or Ni/W stack with magneto-optical detection: the paper's central prediction is a clear in-plane orbital/spin signal, so a null result would refute it.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the underlying LMTO-ASA linear band-theory formalism on which the real-space Hamiltonian is built."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the real-space Chebyshev approach to the Kubo conductivity that the paper adapts for orbital and spin Hall currents."},{"cited_title":"Weiße, G","cited_arxiv_id":null,"evidence_quote":"Provides the kernel polynomial method used to regularize and truncate the Chebyshev expansion."},{"cited_title":"Go, H.-W","cited_arxiv_id":null,"evidence_quote":"Defines the Wannier-interpolation route and the orbital-dependent anomalous position, the comparison point for the paper's direct real-space intra-atomic approach."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives quantitative references for spin and orbital polarization in heavy-metal/3d-metal bilayers for comparison with the accumulation results."}],"review_version":1}