{"id":"52c2455b-7584-4ac1-a428-a59d29eb5ad7","arxiv_id":"2607.28843","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A differentiable forward model for far-, near-, and point-focused high-energy diffraction microscopy matches the MIDAS C simulator to pixel precision and enables joint gradient-based refinement of grain and detector parameters.","lead":"A new PyTorch simulator for high-energy diffraction microscopy computes X-ray spot positions and their derivatives, matching the established MIDAS C code to pixel precision. It enables one-shot gradient-based fitting of grain orientation, strain, and detector geometry across far-field, near-field, and point-focused setups.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"NF pixel-exact agreement rests on unspecified Gaussian-splatting parameters; if tuned to the C reference output, the match is circular and the headline evidence is compromised.","rationale":"The reader's weakest assumption focuses precisely on the unspecified NF splatting parameters, and I agree that this is the most load-bearing concern. The paper's strongest evidence for the differentiable physics is the pixel-exact NF agreement, but if the splatting kernel is not specified, the agreement cannot be assessed or reproduced. The concern is not that the physics is wrong—the FF and pf agreements are on spot coordinates and are well-specified—but that the NF pixel-exact claim, as reported, is not independently verifiable. This does not invalidate the whole paper; it is an addressable documentation/reproducibility gap, consistent with a CONDITIONAL verdict. I would not escalate to REJECT because the underlying forward model is likely correct based on the coordinate-level agreements and real-data round-trip, but the NF claim needs clarification. The proposed concrete test directly determines whether the match is a tuning artifact or robust.","tokens_in":31595,"tokens_out":2996,"duration_ms":37452,"concrete_test":"Inspect the repository script paper/scripts/run_side_by_side_nf.py and the fwd_sim/hedm_forward.py model defaults to recover the exact sigma and radius used for the NF comparison, and check whether they are explicitly listed in the paper. Then re-run the NF cross-code comparison with sigma set to a near-zero (delta-like) splat and with a range of sigma values (e.g., 0.1, 0.5, 1.0 px). If the 2304/2304 exact match persists only for the originally used sigma, the agreement is a fitted artifact; if it holds for all reasonable splatting kernels (because the C reference also renders delta-like spots), the concern is resolved. Additionally, instrument simulateNF on a synthetic single-spot input to characterize its actual pixel point-spread function and compare against the Gaussian splatting kernel.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of the paper is that the differentiable forward model achieves pixel-exact agreement with the established MIDAS reference simulators. The NF-HEDM agreement (2304/2304 pixels identical) is the strongest quantitative evidence for this claim, yet it depends on the Gaussian-splatting output head (Section 4.2.5) whose sigma and radius are described only as 'configurable', with no default or chosen value reported. The C reference rendering in simulateNF (Section S1.2) is described only as 'ray-traces the diffracted beam through the detector plane', with no pixel-level point-spread function or discretization rule. Without knowing the exact splatting parameters used in the comparison, the claim that the PyTorch model reproduces the C image 'pixel-exact' could be an artifact of having tuned those parameters against the reference output. This would make the agreement circular: the rendering head is fitted to match the reference, so the match is not an independent test of the physics. The concern is concrete: if the splatting sigma/radius is not specified and not justified, the reproducibility and evidentiary value of the NF agreement collapse. The reader's weakest-assumption flag on this point is well placed; the paper does not document the parameter values anywhere in the main text or supplement, and the 'configurable' wording leaves open the possibility of post-hoc tuning.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents HEDMForwardModel, a PyTorch implementation of an end-to-end differentiable forward model for far-field (FF), near-field (NF), and point-focused (pf) HEDM. The model is validated against the MIDAS C simulators (ForwardSimulationCompressed and simulateNF), reporting pixel-exact agreement for FF and NF and 99.27% for pf, with analytic gradients verified against finite differences. Three optimization demonstrations are given: joint NF orientation–strain–position recovery, a real-data round-trip on a 214-grain α-Ti dataset, and joint multi-panel FF geometry–grain refinement. The software is released as the open-source midas-diffract package.","tokens_in":31988,"tokens_out":6763,"duration_ms":73452,"significance":"If the validation holds, this is a significant methodological contribution: it opens gradient-based joint refinement, regularized reconstruction, and Bayesian inference for HEDM, with a unified framework across all three HEDM geometries. The paper includes extensive cross-code comparisons, gradient-correctness checks, basin-of-convergence sweeps, and a reproducible-script pipeline. However, the NF pixel-exact agreement rests on undocumented Gaussian-splatting parameters, and the validation references are the authors' own C simulators; these limit the independent evidentiary value of the headline claims.","major_comments":[{"comment":"The NF 'predict_images' output head uses a Gaussian splatting kernel with 'configurable sigma and radius' (§4.2.5), but the values used for the 2304/2304 pixel-exact comparison (§5.3) are never reported. Supplementary §S1.2 describes the C simulateNF reference only as 'ray-traces the diffracted beam through the detector plane', without specifying the pixel-level point-spread function or discretization rule. If the sigma/radius were chosen to match the C output, the agreement is circular and not an independent validation. Please specify the exact splatting parameters, how they were set, and document the C rendering kernel. Without this, the NF agreement is not reproducible and the strongest quantitative evidence for the pixel-exact claim collapses.","section":"§4.2.5, §5.3, §S1.2"},{"comment":"The abstract's claim that residuals match the production fit 'to 0.3%' is not supported by the body. Section 6.4 reports median ω residuals 'essentially tied' (0.0901° vs 0.0904°, a 0.33% difference), while η residuals are 14.7% tighter and position residuals 39% tighter. The 0.3% figure appears to refer only to the ω residual, not to an overall residual match, and it conflates the perturb-and-recover experiment (100% recovery) with the forward-model agreement test. Similarly, the abstract's '~6 nm precision' for NF-HEDM (Section 6.3, Table 2) is the noise-free result; the noisy case gives 66 nm. Please revise the abstract and contribution 3 to state precisely what is measured and under what noise conditions.","section":"Abstract; §6.4, Test 1"},{"comment":"The validation references are the MIDAS C simulators (ForwardSimulationCompressed and simulateNF), which are authored by the same group and, as stated in §9, are documented in the peer-reviewed literature for the first time in this paper. The PyTorch model is a direct port of those codes. The pixel-exact agreement therefore primarily demonstrates that the port is faithful, not that the physics is independently correct. The manuscript should qualify the phrase 'established reference simulators' and clarify that this is a self-consistency check. An independent cross-check against another HEDM code (e.g., ImageD11 or HEXRD spot positions) would materially strengthen the external validity of the central claim.","section":"§3, §5.1, §9"}],"minor_comments":[{"comment":"'~6 nm precision' should be explicitly labeled as the noise-free result; the noisy result is 66 nm. As written, the abstract implies robustness that the body does not support.","section":"Abstract; §6.3, Table 2"},{"comment":"The straight-through estimator for pixel rounding is mentioned but not described. Please provide the exact gradient rule or a reference, since it is essential for the NF image-space gradients.","section":"§4.2.5"},{"comment":"The tilt sweep is reported only as 'up to 1° per axis'. Please state the number of tilt configurations, the step size, and the exact comparison criterion (e.g., max pixel difference or exact index match) so the 2322/2322 result can be reproduced.","section":"§5.3"},{"comment":"The description of simulateNF should include the pixel intensity distribution rule (e.g., bilinear interpolation, as in ForwardSimulationCompressed) and how the detector plane is discretized, to make the NF comparison fully reproducible.","section":"§S1.2"},{"comment":"The text says 'below the diagonal on 179/214 grains for η and position' and '120/214 for ω'. Please ensure the counts and percentages are internally consistent and clearly tied to the specific residual definitions.","section":"§6.4, Test 1"}],"recommendation":"major_revision","confidential_remarks":"The claim to be 'the first end-to-end differentiable HEDM forward model' should be checked carefully against the literature and concurrent arXiv postings. The validation is entirely against the authors' own MIDAS C simulators; an external independent benchmark would substantially increase confidence in the pixel-exact agreement claim. The NF splatting-parameter omission is the key point to resolve before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is the first end-to-end differentiable HEDM forward model covering FF, NF, and pf geometries, and it ships with code, validation scripts, and a real-data demonstration. The core is sound; I'd send it to peer review. But two things need fixing before I'd trust the headline claims: the NF pixel-exact match depends on Gaussian-splatting parameters that aren't reported, and the abstract's \"0.3% residual\" line doesn't match the body's numbers.\n\nWhat's new: the PyTorch implementation makes every input (orientation, position, strain, detector geometry) a leaf in the autograd graph, which unlocks gradient-based joint refinement, Bayesian UQ, and physics-informed regularisation. The validation is extensive: FF 162/162 spots with sub-10^-3 px residuals, NF 2304/2304 pixels matched exactly under zero and non-zero tilts, pf 1088/1096 spots, finite-difference gradient checks at 10^-8, a 250-grain noisy cross-code test, and a perturb-and-recover on a real 214-grain alpha-Ti dataset. The multi-panel joint geometry-and-grain refinement is a nice concrete demonstration. They also release the code and reproducibility scripts, which is real.\n\nSoft spots, in proportion. First, the NF pixel-exact claim: section 4.2.5 says Gaussian splatting uses \"configurable sigma and radius\" but never reports the values, and the C reference's rendering is described only as \"ray-traces\" in S1.2. If those parameters were chosen to force the match, the agreement is circular. This doesn't sink the paper—the FF and pf agreements don't involve splatting—but it needs a clear statement of the chosen parameters and a justification that they're not tuned. Second, the abstract says residuals match the production fit to 0.3%; the body reports median omega tied and eta/position tighter by 14.7% and 39%. That's an overstatement. Third, the validation reference is the group's own MIDAS C code; the agreement is partly self-consistency. Given the C code's production history that's reasonable, but the paper should say so more explicitly. Minor: the pf match is 99.27%, so \"pixel-exact\" should be qualified.\n\nBottom line: the central claim—that a differentiable HEDM forward model can be built without sacrificing agreement with the reference physics—holds up. This is a paper for the synchrotron materials-science audience and for method developers; it deserves a serious referee, with the splatting documentation and abstract fix as required revisions.","headline":"Genuinely useful first differentiable HEDM forward model covering FF/NF/pf, with solid core validation, but the NF pixel-exact claim needs the splatting parameters reported and the abstract's 0.3% residual line corrected.","tokens_in":32441,"tokens_out":2717,"would_cite":true,"duration_ms":29206,"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":"This paper claims that high-energy diffraction microscopy can be made fully differentiable end to end—across far-field, near-field, and point-focused geometries—without sacrificing pixel-level agreement with the established C reference simu","keywords":["high-energy diffraction microscopy","differentiable forward model","automatic differentiation","grain orientation refinement","detector geometry calibration","near-field diffraction","far-field diffraction","point-focused HEDM"],"falsifier":"Inspect the test configuration for the near-field comparison and confirm that the splatting sigma and radius are set from an independent physical or numerical criterion; then re-run the 2304-pixel comparison after varying sigma across a small range. If the exact match holds only at the chosen sigma, the pixel-exact claim is parameter-dependent rather than a robust property of the physics port.","tokens_in":31479,"feed_emoji":"🔬","tokens_out":4619,"duration_ms":50467,"temperature":0.7,"pith_summary":"The paper sets out to close a structural gap in high-energy diffraction microscopy (HEDM): the forward models that map grain state and instrument geometry to detector images are fast but not differentiable, so gradient-based joint refinement, physics-informed regularization, and Bayesian uncertainty quantification are unavailable. It presents a differentiable forward model built from standard tensor operations, covering all three HEDM geometries, and validates it against the production C simulators to pixel-level agreement. With exact gradients available, it demonstrates joint recovery of orientation, strain, and position in near-field geometry to about six nanometers, and joint refinement of all per-detector geometry parameters plus per-grain state on a four-panel far-field setup. If the claims hold, HEDM reconstruction moves from sequential, coordinate-descent-style fitting to single gradient flows over the entire experimental state.","feed_headline":"X-ray grain mapping now differentiable end to end","feed_subtitle":"Far-, near-, and point-focused geometries match C references to pixel level, unlocking joint gradient refinement.","key_machinery":"The load-bearing object is the omega quadratic: the rotating-crystal diffraction condition reduced to a quadratic equation in cos(omega), solved analytically, with branch selection handled by masked tensor operations rather than branching. Around it, four substitutions preserve differentiability: masked selection between two fully evaluated paths in place of conditional branching, clamped inverse trigonometric functions, Newton–Schulz projection onto SO(3) instead of SVD, and binary validity masks multiplied into losses. Near-field images are produced by Gaussian splatting with a straight-through pixel estimator so gradients flow through sub-pixel spot positions.","core_discovery":"The central discovery is that every step of the HEDM forward simulation—rotation-matrix construction, reciprocal-lattice strain, the quadratic rotating-crystal diffraction condition, azimuthal angle, and detector projection—can be re-expressed as differentiable tensor primitives without changing the physics. The resulting model reproduces the validated C simulators to within floating-point precision on the vast majority of reflections (162/162 far-field, 2304/2304 near-field pixels, 1088/1096 point-focused), with the only mismatches at a known branch-selection boundary near the azimuthal axis. Because every input is a leaf of the autograd graph, scalar losses on predicted spot positions prop","pith_inferences":["The pixel-exact near-field claim may be sensitive to the Gaussian-splatting bandwidth, which the paper does not specify; a fair test would measure whether the exact match survives when the splatting sigma is fixed by an independent point-spread calibration.","The eight unmatched point-focused spots sit at a genuine physical degeneracy; this suggests that a continuous branch mixture rather than a hard choice would make gradients stable across that boundary.","The demonstrated 10–15 degree basin could potentially be widened by soft spot-matching losses, a direction the paper mentions but does not explore.","If these results replicate, HEDM becomes a testbed for differentiable experimental design: because detector geometry is differentiable, one could optimise detector placement or the omega sweep to maximise information gain before the experiment runs."],"forward_implications":["HEDM inverse problems can be solved by joint gradient descent over orientation, strain, position, and detector geometry in one flow, replacing sequential fits.","The same gradients enable physics-informed regularisation and Bayesian uncertainty quantification via Hamiltonian Monte Carlo or variational inference.","Detector auto-calibration becomes possible in near-field geometry, where no powder ring exists, recovering tilt angles to sub-microradian precision.","The basin of convergence for orientation refinement is about 10–15 degrees, matching typical production indexer accuracy, so existing indexers can feed this refiner.","The model scales to thousands of grains and can be coupled to differentiable finite-element solvers, closing the experiment–simulation loop for crystal plasticity."],"fun_headline_variants":["End-to-end differentiable HEDM forward model released","HEDM simulation becomes differentiable across all geometries","Differentiable X-ray grain mapping matches MIDAS to pixel level","Gradient-based HEDM refinement now possible with new model","First differentiable forward model for HEDM in PyTorch"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The near-field 'pixel-exact' agreement depends on the Gaussian-splatting output head reproducing exactly the same pixel values as the C reference, but the paper does not state the splatting sigma and radius used; if those were tuned to force the match, that agreement is not independent evidence.","fun_headline_variants_meta":{"raw":{"variants":["End-to-end differentiable HEDM forward model released","HEDM simulation becomes differentiable across all geometries","Differentiable X-ray grain mapping matches MIDAS to pixel level","Gradient-based HEDM refinement now possible with new model","First differentiable forward model for HEDM in PyTorch"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000187,"raw_usage":{"total_tokens":1188,"prompt_tokens":790,"completion_tokens":398,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":534,"completion_tokens_details":{"reasoning_tokens":317}},"tokens_in":534,"tokens_out":398,"duration_ms":5113,"temperature":1.0,"reasoning_tokens":317,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T01:30:59.618641+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Inspect the test configuration for the near-field comparison and confirm that the splatting sigma and radius are set from an independent physical or numerical criterion; then re-run the 2304-pixel comparison after varying sigma across a small range. If the exact match holds only at the chosen sigma, the pixel-exact claim is parameter-dependent rather than a robust property of the physics port.","supporting_citations":[],"review_version":1}