{"id":"56505b48-9e2d-4cad-a006-9aaf2068aa0c","arxiv_id":"2607.00756","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":6.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":0,"one_line_summary":"BiP-PRISM extends the PRISM method with beam partitioning and on-atom interpolation to cut memory use by 5x while controlling error via a locality result for core-loss EELS simulations.","lead":"The BiP-PRISM algorithm speeds up core-loss STEM-EELS simulations by partitioning beams in the probe-forming and detector matrices and using local natural-neighbor interpolation at ionized atoms. If it works, this would let researchers run full-resolution elemental maps and 4D data on ordinary GPUs instead of large clusters.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"Central claim depends on locality result that total error is governed entirely by on-atom reconstruction error","rationale":"The reader's weakest_assumption matches the load-bearing premise exactly. The abstract-only basis leaves the locality result unexamined, so the verdict moves from UNVERDICTED to CONDITIONAL pending verification of that result; no other internal inconsistency is visible from the given material.","tokens_in":1744,"tokens_out":299,"duration_ms":27257,"concrete_test":"From the paper's error-analysis or methods section, extract the locality result and re-derive it without assuming the on-atom dominance; then run a control simulation on the FePt nanoparticle case with doubled thickness and compare whether total error remains equal to the on-atom reconstruction error within 5%.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The abstract states that a locality result shows total error is governed entirely by on-atom reconstruction error from natural-neighbor interpolation on sparse parent beams. This premise is what justifies removing per-scan exit-wave propagation for S1 and S2 matrices and the claimed 5x memory reduction while retaining accuracy for elemental maps and qEELS. No derivation, assumptions (e.g., probe delocalization bounds, neglect of long-range dynamical scattering), or quantitative error bound appear in the supplied abstract; if the locality fails outside the characterized regime (thick samples, high-q transfer), the error is not controlled solely by the local interpolation step.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The paper introduces the BiP-PRISM algorithm as an extension of PRISM for core-loss STEM-EELS simulations. It partitions the probe-forming (S1) and detector-propagating (S2) matrices onto sparse parent beams, reconstructs them at ionized atoms via natural-neighbor interpolation, and invokes a locality result asserting that total simulation error is controlled solely by this local reconstruction error. This permits removal of per-scan exit-wave propagation, yielding ~5x memory reduction while enabling full-resolution elemental maps, 4D cubes, and qEELS on consumer GPUs. The validity regime is characterized and the method is demonstrated on a multimodal five-edge oxide interface and an FePt nanoparticle Fe-L edge map, reporting high accuracy relative to reference calculations.","tokens_in":1870,"tokens_out":557,"duration_ms":27516,"significance":"If the locality result is rigorously established and the error remains localized under the stated conditions, the work would meaningfully advance quantitative EELS modeling by making previously intractable large-scale simulations practical. The reported 5x memory savings and GPU demonstrations directly address a recognized computational bottleneck in the field, with potential impact on routine analysis of interfaces and nanoparticles. The explicit partitioning and interpolation steps constitute a clear algorithmic contribution.","major_comments":[{"comment":"The locality result (abstract and method section) is load-bearing for the central claim that per-scan S1/S2 propagation can be eliminated while retaining controlled error. The manuscript must supply the full derivation, including all assumptions (probe delocalization bounds, neglect of long-range dynamical scattering, thickness limits) and a quantitative error bound showing that off-atom contributions are negligible; without this, the justification for skipping full propagation and the 5x memory claim cannot be verified.","section":"Abstract and locality-result section"},{"comment":"Demonstration cases (oxide-interface and FePt maps): the reported accuracy must be accompanied by explicit error metrics (e.g., pixel-wise L2 or edge-integrated intensity differences) versus a full PRISM reference, together with the precise beam sparsity and interpolation parameters used, so that the claimed “high accuracy” can be assessed against the locality premise.","section":"Results section"}],"minor_comments":[{"comment":"Notation for the partitioned matrices S1 and S2 should be defined once with explicit dimensions and beam indices before the interpolation step is introduced.","section":null},{"comment":"The validity-regime characterization would benefit from a single summary table or plot showing error versus sample thickness, probe convergence angle, and q-transfer range.","section":null}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for their careful reading and constructive comments on our manuscript. We address each major comment point by point below.","responses":[{"response":"We agree that the locality result is central and that its full derivation, assumptions, and quantitative bound must be provided for the claims to be verifiable. The current manuscript sketches the result and states the governing principle but does not contain the complete derivation. In the revised version we will expand the methods section to include the full mathematical derivation together with explicit statements of all assumptions (probe delocalization bounds, neglect of long-range dynamical scattering, and thickness limits) and a quantitative error bound demonstrating negligibility of off-atom contributions.","revision_made":"yes","referee_comment":"[Abstract and locality-result section] The locality result (abstract and method section) is load-bearing for the central claim that per-scan S1/S2 propagation can be eliminated while retaining controlled error. The manuscript must supply the full derivation, including all assumptions (probe delocalization bounds, neglect of long-range dynamical scattering, thickness limits) and a quantitative error bound showing that off-atom contributions are negligible; without this, the justification for skipping full propagation and the 5x memory claim cannot be verified."},{"response":"We agree that explicit quantitative error metrics and parameter values are necessary to allow readers to assess the reported accuracy against the locality premise. The revised manuscript will add pixel-wise L2 differences and edge-integrated intensity differences relative to full PRISM references for both demonstration cases, together with the precise beam-sparsity factors and natural-neighbor interpolation parameters used in each simulation.","revision_made":"yes","referee_comment":"[Results section] Demonstration cases (oxide-interface and FePt maps): the reported accuracy must be accompanied by explicit error metrics (e.g., pixel-wise L2 or edge-integrated intensity differences) versus a full PRISM reference, together with the precise beam sparsity and interpolation parameters used, so that the claimed “high accuracy” can be assessed against the locality premise."}],"tokens_in":1430,"tokens_out":437,"duration_ms":28401,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The new piece here is applying beam partitioning to both the S1 probe-forming and S2 detector-propagating PRISM matrices, then reconstructing only at the ionized atom with natural-neighbor interpolation. That combination is not in the earlier PRISM papers referenced, and it lets the method drop per-scan exit-wave propagation while claiming a 5x memory cut.\n\nThe demonstrations are the strongest part. The oxide-interface five-edge map and the FePt nanoparticle Fe-L map both run at full resolution on consumer GPUs with reported high accuracy. Those are concrete, relevant test cases for people who actually need elemental or momentum-resolved EELS data.\n\nThe soft spot is the locality result. The abstract states that total error is governed entirely by the on-atom reconstruction error, which is what justifies skipping full propagation. No derivation, assumptions about probe delocalization, or bounds on long-range scattering appear in the supplied text, so it is hard to judge how far the validity regime extends beyond the two examples. If that step does not hold for thicker specimens or high-q transfer, the error control claim weakens.\n\nThis paper is for groups already running or wanting to run quantitative dynamical EELS simulations at scale. The algorithmic steps are explicit enough that a reader could implement and test them. It deserves a serious referee because the practical payoff is clear and the demos are on real materials problems, even though the central error argument will need close checking in review.","headline":"BiP-PRISM adds dual-matrix beam partitioning plus local natural-neighbor interpolation to PRISM, which could cut memory and propagation costs for core-loss EELS, but the locality claim that controls all error needs a clear derivation.","tokens_in":2342,"tokens_out":378,"would_cite":false,"duration_ms":23274,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"Beam partitioning in PRISM matrices lets core-loss STEM-EELS simulations skip per-scan propagation while keeping error local to each atom.","keywords":["STEM-EELS","PRISM algorithm","core-loss simulation","beam partitioning","dynamical diffraction","GPU acceleration","atomic-resolution mapping"],"falsifier":"Run both the full per-scan PRISM calculation and the BiP-PRISM version on the same atomic model for one probe position and check whether the difference in the core-loss spectrum exceeds the bound given by the on-atom interpolation error alone.","tokens_in":2630,"feed_emoji":"🔬","tokens_out":685,"duration_ms":23441,"temperature":0.7,"pith_summary":"The paper develops the BiP-PRISM algorithm to make quantitative dynamical simulations of atomic-resolution core-loss electron energy loss spectroscopy practical on ordinary hardware. It splits the probe-forming and post-loss scattering matrices into sparse parent beams that are computed once and then interpolated locally around each ionized atom. A locality result shows that the overall error is controlled only by how well those local interpolations match the full matrices at the atom sites. This change removes the need to propagate the exit wave separately for every probe position and cuts memory use by a factor of five. The authors demonstrate the approach on a five-edge oxide interface and an FePt nanoparticle map while staying within the validity regime they map out.","feed_headline":"Beam partitioning cuts memory use 5x in core-loss EELS simulations","feed_subtitle":"BiP-PRISM interpolates sparse matrices locally at atoms and removes per-scan propagation for full 4D maps on consumer GPUs.","key_machinery":"Beam partitioning of the probe-forming (S1) and detector-propagating (S2) PRISM matrices with natural-neighbor interpolation on sparse parent beams at each ionized atom.","core_discovery":"By calculating the S1 and S2 PRISM matrices only on a sparse set of parent beams and reconstructing the values at each ionized atom via natural-neighbor interpolation, the BiP-PRISM method removes per-scan exit-wave propagation; the total error is governed entirely by the on-atom reconstruction error, which the authors bound and validate on multimodal five-edge and Fe-L maps at 5x lower memory cost.","pith_inferences":["The same locality principle could be tested on thicker specimens or different acceleration voltages to map the practical range further.","Integration with experimental workflows would let users decide scan density on the fly based on the reported error bound.","The approach might reduce the barrier to simulating momentum-resolved data for materials where only average spectra were previously affordable."],"forward_implications":["Full-resolution elemental mapping, 4D data cubes, and momentum-resolved qEELS become feasible on consumer-grade GPUs.","Memory footprint drops by a factor of five while accuracy remains high enough for the demonstrated oxide-interface and nanoparticle cases.","The validity regime of the approximation is characterized so users can decide when the method applies.","Multimodal simulations involving multiple core-loss edges can be performed without prohibitive cost."],"fun_headline_variants":["BiP-PRISM cuts core-loss EELS memory 5x with beam partitioning","Sparse S1 S2 interpolation enables 4D EELS maps on GPUs","BiP-PRISM skips per-scan propagation in scalable EELS","Error bound validates BiP-PRISM for multimodal EELS mapping"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The total error in the simulated EELS signal is governed entirely by the reconstruction error at the ionized atom sites.","fun_headline_variants_meta":{"raw":{"variants":["BiP-PRISM cuts core-loss EELS memory 5x with beam partitioning","Sparse S1 S2 interpolation enables 4D EELS maps on GPUs","BiP-PRISM skips per-scan propagation in scalable EELS","Error bound validates BiP-PRISM for multimodal EELS mapping"]},"model":"grok-4.3","cost_usd":0.006441,"raw_usage":{"total_tokens":3006,"prompt_tokens":644,"num_sources_used":0,"completion_tokens":76,"cost_in_usd_ticks":64412000,"prompt_tokens_details":{"text_tokens":644,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":2286,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":644,"tokens_out":76,"duration_ms":20155,"temperature":1.0,"reasoning_tokens":2286,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-02T10:14:00.308733+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"Run both the full per-scan PRISM calculation and the BiP-PRISM version on the same atomic model for one probe position and check whether the difference in the core-loss spectrum exceeds the bound given by the on-atom interpolation error alone.","supporting_citations":[],"review_version":1}