{"id":"d48bcbc9-5d1b-49c2-96ec-681c61d211d0","arxiv_id":"2607.11209","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.5,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Heteroscedastic GPR recovers morphology-induced, input-dependent stiffness scatter in spinodoids from one-realization data, and only heteroscedastic RBDO meets reliability targets that deterministic and homoscedastic designs miss.","lead":"Identical cone-angle settings for spinodoid metamaterials produce different stiffnesses because Gaussian random fields generate different morphologies. Heteroscedastic Gaussian process models recover that input-dependent scatter from sparse data and show only reliability-based design that uses it meets failure targets.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.5","headline":"Recovery of input-dependent aleatoric variance from one-realization residuals under the parametric noise model is the softest link for the RBDO claim.","rationale":"The reader correctly isolates the residual-to-aleatoric recovery step (distributional continuity + parametric form) as the weakest assumption underlying the framework that produces the designs in Table 3. The post-hoc 1000-realization validation itself is solid and independently confirms the three reported pf intervals, so the specific numerical claim holds. No stronger internal inconsistency or calculation error was found; the physical interpretability of the fitted coefficients and the 1D–3D noise-trend consistency supply modest corroboration. Because the method’s ability to locate reliable designs in general still rests on that recovery step (and on the unmodeled component correlations), the CONDITIONAL verdict and the call for noise-form sensitivity checks remain appropriate. No change to the reader’s assessment is warranted.","tokens_in":20128,"tokens_out":625,"duration_ms":40387,"concrete_test":"Using the same 1000 GRF+homogenization realizations already generated for the Clopper–Pearson intervals in Table 3, compute the empirical standard deviations of E1, E2, E3 (and of the two ratios) at each of the three optimized designs; compare them component-wise to the HOMO and HETERO surrogate-predicted σ̂i(Θ*). If heteroscedastic σ̂ matches empirical std within ~20 % (while constant noise does not) and the independent-ξ MCS ranking of the three designs remains consistent with the true joint pf, the recovery assumption is supported; larger mismatch would show the noise model is unreliable for guiding RBDO.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The claim that only heteroscedastic RBDO meets the true pf target (Table 3) depends on the surrogate correctly locating low-scatter designs during optimization. That location is driven by the inferred noise function σ^{2}(Θ) (Eqs. 6–10 and S1.2), recovered solely from residual patterns of sparse one-realization-per-point data via the distributional-continuity argument (§4.2) and the d=1 polynomial-exponential form chosen by ELPD/parsimony. If the recovered σ systematically mis-ranks high- versus low-scatter regions, the optimizer can return a design whose true morphology-induced failure probability exceeds 0.01 while the surrogate reports compliance; the three-point 1000-realization check then becomes an insufficient guarantee for the method. Separate single-output GPs further ignore correlations among stiffness components that govern the ratio constraints G (Eq. 17), an approximation the authors themselves flag in §5.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript reinterprets cone-angle descriptors of GRF-generated spinodoid metamaterials as stochastic descriptors associated with input-dependent effective-stiffness distributions, rather than deterministic structure–property maps. It trains heteroscedastic Gaussian process surrogates (polynomial-exponential noise, Eqs. 6–10) on sparse one-realization-per-point homogenization data, without empirical variance labels, and shows that inferred aleatoric scatter is component-specific and concentrated at small cone angles. Calibration and uncertainty metrics improve over a homoscedastic baseline in 1D and 3D (Tables 1–2, Figs. 3–4). The calibrated surrogates are then used in DDO and RBDO for a target axial modulus with stiffness-ratio constraints; independent 1000-realization Clopper–Pearson validation (Table 3, §4.4) shows that only the heteroscedastic RBDO design meets pf_target = 0.01, while the deterministic and homoscedastic optima violate it.","tokens_in":20486,"tokens_out":1398,"duration_ms":33431,"significance":"If the results hold, the work makes a clear case that morphology-induced aleatoric uncertainty in spinodoid design is structured, input-dependent, and consequential for reliability-based inverse design—not a constant residual to be averaged out. Strengths include: (i) physically interpretable noise coefficients aligned with mechanically active directions; (ii) consistent d=1 noise structure across 1D and 3D with ELPD-based degree selection (S1.2, Fig. S1); (iii) calibration diagnostics (Fig. 4) and 1D–3D noise-trend consistency (Fig. 6); and (iv) independent multi-realization validation of the three optimized designs with Clopper–Pearson intervals (Table 3), which is stronger evidence than surrogate-only reliability claims. The demonstration that a homoscedastic RBDO can report compliance while failing true reliability is practically important for the spinodoid and broader stochastic-architecture communities.","major_comments":[{"comment":"§3, §4.2 and S1.2: Recovery of input-dependent aleatoric variance σ²(Θ) from residual patterns of sparse one-realization-per-point data (distributional-continuity argument; polynomial-exponential form with d=1 forced by parsimony even when ELPD rises slightly for C1133 and C2323) is load-bearing for the RBDO claim. Fig. 3 reports Emp (n=3) only as a local diagnostic, and Table 3 validates reliability only at three optimized points. Because the optimizer’s preference for low-scatter designs is driven by the inferred σ(Θ), the manuscript should add a denser held-out comparison of multi-realization empirical standard deviations versus predicted σ across the design space (or a clear quantitative caveat that ranking of uncertainty regions is only sparsely validated). Without that, the general claim that heteroscedastic RBDO is essential rests on a soft link between residual-inferred noise and","section":"§3, §4.2, S1.2, Table 3"},{"comment":"§4.4, Eqs. (11)–(17): The limit-state G is defined on stiffness ratios E2/E1 and E3/E1, so the failure probability depends on the joint distribution of (E1, E2, E3). The paper trains independent single-output GPs and samples ξi independently, which neglects physical correlations among tensor components. Section 5 correctly flags multi-output GPR as future work, but this approximation is load-bearing for the present RBDO numbers (own-model pf and, to a lesser extent, interpretation of why HETERO succeeds). A sensitivity check (e.g., copula or residual-correlation sampling) or an explicit statement in §4.4 that reported pf values are under an independence assumption—and may be biased relative to the joint law—should be added before the reliability conclusions are treated as fully quantitative.","section":"§4.4, Eqs. (11)–(17)"}],"minor_comments":[{"comment":"Abstract and §4.4 wording: the abstract states that only the heteroscedastic formulation satisfies the target “under the heteroscedastic uncertainty evaluation,” whereas Table 3’s decisive evidence is independent computational-homogenization validation (pf_CP). Align the abstract with the stronger, model-independent result.","section":"Abstract, §4.4"},{"comment":"Table 2 reports relative improvements (%) without the absolute homoscedastic baselines that Table 1 provides; adding absolute Homo/Hetero values (or a supplementary table) would make the 3D gains easier to interpret, especially where MAE slightly worsens.","section":"Table 2"},{"comment":"S1.3: The power-law warping exponent 1.6 for training LHS is a free design choice that concentrates samples at small cone angles. A brief sensitivity note (or one alternative exponent) would reassure readers that noise-model conclusions are not an artifact of that warping.","section":"S1.3"},{"comment":"Notation: several symbols render inconsistently in the manuscript text (e.g., script C for stiffness, Greek letters appearing as doubled characters). Clean typesetting of ℂ_eff, σ²(Θ), and the cone-angle set S_Θ would improve readability.","section":"Throughout"},{"comment":"Acknowledgements are left as a placeholder (“This work supported by ~”); complete before production.","section":"Acknowledgements"},{"comment":"Fig. 3 caption: “Emp. (n=3)” is useful but the local smoothing method is not specified; one sentence in the caption or S1 would help reproducibility.","section":"Fig. 3"}],"recommendation":"major_revision","confidential_remarks":"Solid fit for a computational mechanics / architected-materials venue. The application of heteroscedastic GPR (building on Ozbayram et al.) to spinodoids plus the RBDO demonstration with independent validation is a genuine contribution, not merely incremental method transfer. The two major points are fixable with additional multi-realization checks and clearer joint-distribution caveats; I would not reject on novelty or scope grounds."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The punchline is practical and well supported: for GRF spinodoids, identical cone-angle sets produce input-dependent stiffness scatter, and only a heteroscedastic RBDO formulation meets a pf=0.01 target under fresh 1000-realization Clopper–Pearson checks. Deterministic and homoscedastic optima both fail that check (Table 3).\n\nWhat is new is not heteroscedastic GPR itself (they cite Ozbayram et al. 2024 correctly) but the application: reinterpreting cone angles as stochastic descriptors, recovering noise coefficients that line up with mechanically active directions, the 1D–3D consistency slice, and the three-way optimization comparison with independent homogenization validation. The diagnostics are careful—calibration curves, noise-vs-residual plots, ELPD degree selection, and the physical reading of the a_rj coefficients against load paths. Mean trends match the expected U-shaped / inverse-U behavior from the wave-vector construction. Circularity is low: variance is inferred from residuals, not labeled by the same quantity later claimed as prediction, and the failure probabilities use new GRF draws.\n\nSoft spots are real but secondary. The noise model is parametric (d=1 forced by parsimony even when ELPD edges up for two components), and the distributional-continuity argument that lets one-realization residuals recover true aleatoric variance is an assumption, not a proof. Separate single-output GPs ignore correlations among stiffness components that enter the ratio constraints; the authors flag this themselves in §5. No code or data release. None of these overturn the central empirical claim on the three optimized designs, but they do mean the method’s ranking of low-scatter regions is only as good as the residual-inferred σ(Θ).\n\nThis is for people who design or optimize GRF-based metamaterials and care about reliability, not for a general ML audience. Math and citations look standard and honest. I would send it to peer review; a serious referee can push on the noise form and multi-output extension without the paper collapsing. Worth engaging if you work in this subfield.","headline":"Solid applied paper: heteroscedastic GPR on spinodoid cone-angle space plus a clean RBDO head-to-head that only the input-dependent model survives independent validation.","tokens_in":21056,"tokens_out":523,"would_cite":true,"duration_ms":5106,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Cone-angle designs for spinodoid metamaterials carry input-dependent stiffness scatter; only heteroscedastic uncertainty models meet reliability targets in inverse design.","keywords":["Spinodoid metamaterials","Gaussian process regression","Heteroscedasticity","Uncertainty quantification","Reliability-based design optimization","Structure-property mapping","Aleatoric uncertainty"],"falsifier":"At many held-out cone-angle points, generate dozens of independent GRF realizations each, compute empirical stiffness standard deviations, and check whether they systematically match the heteroscedastic σ(Θ) inferred from one-realization training; large mismatches in magnitude or in which cone-angle regions are high-scatter would falsify the recovery claim.","tokens_in":21043,"feed_emoji":"⚙️","tokens_out":932,"duration_ms":19309,"temperature":0.7,"pith_summary":"Spinodoid metamaterials are built from Gaussian random fields, so the same cone-angle parameters can produce many different morphologies and different effective stiffnesses. The paper shows that this scatter is not uniform: its width depends on location in design space and on which stiffness component you measure, tracking each component's load-bearing directions. Heteroscedastic Gaussian process regression can recover that input-dependent uncertainty from sparse data with only one realization per design point, without needing empirical variance labels everywhere. When that uncertainty is used in reliability-based design optimization, a deterministic optimum and a constant-noise (homoscedastic) formulation both fail the reliability target under direct multi-realization checks; only the heteroscedastic formulation satisfies it. Readers who design or certify these materials should care because mean-only structure-property maps can look successful while still producing designs that routinely violate constraints once real morphology variability is present.","feed_headline":"Only input-dependent uncertainty yields reliable spinodoid designs","feed_subtitle":"Deterministic and constant-noise optima both violate reliability targets once morphology scatter is measured.","key_machinery":"Heteroscedastic Gaussian process regression with a polynomial-exponential noise model that jointly learns the mean structure-property map and an input-dependent residual variance σ²(Θ) from sparse one-realization-per-point observations, separating aleatoric (realization-induced) from epistemic uncertainty.","core_discovery":"Cone-angle descriptors for GRF-generated spinodoids are stochastic descriptors of input-dependent property distributions, not deterministic maps to single stiffness values. Heteroscedastic GPR infers the morphology-induced aleatoric uncertainty from one-realization-per-point data, and only when that heteroscedastic uncertainty is used in RBDO does the design meet the prescribed reliability target under independent multi-realization validation; deterministic and homoscedastic alternatives do not.","pith_inferences":["The same one-realization residual-pattern approach could transfer to other stochastic generative metamaterial families where full variance labels are expensive.","Joint multi-output models that capture correlations among stiffness components would likely improve failure-probability estimates for ratio constraints beyond independent surrogates.","Active sampling concentrated in high-scatter cone-angle regions could cut the data needed to calibrate the noise map.","Fabrication tolerances on realized cone angles would add a second aleatoric source the same framework could absorb without changing its core structure."],"forward_implications":["Deterministic inverse design that ignores morphology scatter is highly susceptible to constraint violation once variability is accounted for.","Design points with identical mean stiffness can carry very different aleatoric uncertainty, so mean-only selection leaves fabrication risk uncontrolled.","Homoscedastic RBDO can meet its own-model reliability target yet fail when evaluated against true input-dependent scatter.","Uncertainty-aware surrogates are required for reliability-aware inverse design of spinodoid metamaterials.","Scatter magnitude and its decay with cone angle follow the mechanically active directions of each stiffness-tensor component."],"fun_headline_variants":["Heteroscedastic GPR maps spinodoid stiffness scatter for reliable RBDO","Only input-dependent uncertainty meets spinodoid reliability targets","Cone-angle descriptors carry stochastic morphology noise not fixed stiffness","Deterministic and homoscedastic optima violate spinodoid reliability under scatter","Sparse one-realization data suffice for heteroscedastic spinodoid surrogates"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"That residual patterns from single realizations at neighboring design points, under a simple linear-in-exponent noise model, recover the true input-dependent morphology scatter without needing repeated samples as variance labels.","fun_headline_variants_meta":{"raw":{"variants":["Heteroscedastic GPR maps spinodoid stiffness scatter for reliable RBDO","Only input-dependent uncertainty meets spinodoid reliability targets","Cone-angle descriptors carry stochastic morphology noise not fixed stiffness","Deterministic and homoscedastic optima violate spinodoid reliability under scatter","Sparse one-realization data suffice for heteroscedastic spinodoid surrogates"]},"model":"grok-4.5","effort":"low","cost_usd":0.006738,"raw_usage":{"total_tokens":1739,"prompt_tokens":828,"num_sources_used":0,"completion_tokens":96,"cost_in_usd_ticks":67380000,"prompt_tokens_details":{"text_tokens":828,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":815,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":828,"tokens_out":96,"duration_ms":7468,"temperature":1.0,"reasoning_tokens":815,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-14T06:15:41.633668+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"At many held-out cone-angle points, generate dozens of independent GRF realizations each, compute empirical stiffness standard deviations, and check whether they systematically match the heteroscedastic σ(Θ) inferred from one-realization training; large mismatches in magnitude or in which cone-angle regions are high-scatter would falsify the recovery claim.","supporting_citations":[],"review_version":1}