Chemical short-range order (CSRO) in Co-Ni-V alloys will modulate the adsorption energy distributions (AEDs) of CO2 and H2 on alloy nanocatalyst surfaces, directly altering catalytic selectivity in CO2 hydrogenation by >=20% as predicted by machine-learned force fields and validated via high-throughput computational campaigns.
Adversarial Debate Score
66% survival rate under critique
Expert panel critique
Independent views, each critiquing the hypothesis on its own — the score rewards genuine disagreement and discounts consensus.
Supporting Research Papers
- Chemical Short-Range Order Regulates Hydrogen Energetics and Hydrogen-Dislocation Interactions in CoNiV
Chemical short-range order (CSRO) has emerged as a critical structural feature in concentrated alloys, yet its coupling with hydrogen remains an active discussion. Here, we develop a machine-learning ...
- Chemical short-range order controls deformation pathways in a complex concentrated alloy
Chemical short-range order (CSRO) is an intrinsic feature of complex concentrated alloys (CCAs), yet its influence on deformation mechanisms is controversial because of the inconclusive state of concu...
- Selectivity- and Activity-Aware Catalyst Descriptors for CO₂ Hydrogenation on Alloy Nanocatalysts using Machine-Learned Force Fields
Adsorption energy distributions (AEDs) have emerged as a powerful and increasingly adopted descriptor for catalytic performance in high-entropy alloys and, more recently, in conventional metallic allo...
- Interaction energies of H₂ and CO on transition-metal surfaces computed by a range-separated hybrid van der Waals density functional
Dissociative chemisorption (DC) of H₂ on the Cu(111) surface is a prototypical problem for understanding elements of heterogeneous catalysis [Science 326, 832 (2009)]. The challenge lies in modeling t...
Computational Result
An LLM's reading of the literature — not computational verification.
CSRO may influence catalytic selectivity, but evidence is mixed.
Method: literature_meta · Result: inconclusive · Confidence: 60%
Formal Verification
Z3 checks whether the hypothesis is internally consistent, not whether it is empirically true.
This discovery has a Claude-generated validation package with a full experimental design.
Precise Hypothesis
For Co-Ni-V ternary alloy nanocatalysts (compositions within Co(30-50)Ni(30-50)V(10-30) at.%), varying the degree of chemical short-range order (CSRO, quantified via Warren-Cowley parameters α_ij computed over 1st-3rd coordination shells) from a quenched/random state (α≈0) to an annealed/ordered state (|α|>0.15 for at least one pairwise correlation) will shift the adsorption energy distribution (AED) means of CO2 and H2 on the dominant exposed facets (111, 100, 211) by ≥0.15 eV and/or broaden AED standard deviation by ≥30%, as computed by a machine-learned interatomic potential (MLIP, e.g., MACE/NequIP-class model trained on DFT data) combined with a hybrid Monte Carlo/molecular dynamics (MC/MD) sampling campaign over ≥500 distinct local surface motifs per CSRO state. This AED shift will produce a computed change in CO2 hydrogenation selectivity (CO vs. CH4 vs. CH3OH branching ratio, estimated via microkinetic modeling using the MLIP-derived energetics) of ≥20% relative selectivity change between the disordered and ordered CSRO states, at fixed T=473K, P=20 bar, H2:CO2=3:1.
- If AED mean shift between ordered and disordered CSRO states is <0.05 eV and AED std-dev change is <10% across all three facets, for both CO2 and H2, the CSRO→AED link is disproven at the stated composition/size window.
- If microkinetic-model-predicted selectivity change is <10% despite a confirmed AED shift ≥0.15 eV (i.e., AED changes do not propagate to selectivity), the "directly alters selectivity" claim is disproven even if the AED-modulation sub-claim holds.
- If MLIP validation against held-out DFT fails accuracy thresholds (MAE >50 meV/atom), results are inconclusive rather than disproving — must be flagged as a methodology failure, not a hypothesis failure.
- If Warren-Cowley parameters cannot be pushed beyond |α|<0.10 via any accessible annealing/quench protocol in the simulated composition space (i.e., no experimentally realizable CSRO contrast exists), the hypothesis is untestable as stated and must be narrowed.
Spine & Adversarial Read
- highAdsorption energy shifts computed from a mean-field microkinetic model built on scaling relations are notoriously insensitive to input energy perturbations — the 20% selectivity threshold may never be reachable regardless of true AED shift, making the hypothesis's causal chain (AED→selectivity) methodologically fragile.Protocol includes an explicit disproof criterion for AED-selectivity decoupling and an abort checkpoint at Day 90; however, the EVP does not yet justify why a mean-field microkinetic model (vs. kinetic Monte Carlo, which captures site-heterogeneity/AED effects more faithfully) was chosen — this is a real methodology justification gap that should be addressed by running a KMC cross-check on the MVP composition before full-scale commitment.
- highCSRO in real nanoparticles synthesized under practical annealing/quenching conditions may never achieve the |Δα|≥0.15 contrast assumed achievable in simulation, especially at 2-10nm particle sizes where surface energy minimization can override bulk-like SRO ordering tendencies.Partially addressed via disproof criterion (Checkpoint 2) and the boundary condition on metastability, but no experimental literature citation is provided (search unavailable) to confirm CSRO magnitudes are experimentally realizable in this specific ternary at this size regime — this remains an open validation gap pending literature/EXAFS confirmation.
- mediumWhy MACE/NequIP-class MLIPs specifically, and why this particular DFT functional (PBE+D3/RPBE) rather than a hybrid functional or higher-level method, given that adsorption energetics on transition-metal alloys are known to be functional-sensitive by 0.1-0.3 eV — comparable to the claimed AED shift itself?Not resolved in this EVP: functional sensitivity is of the same order as the effect being measured, meaning results could be an artifact of DFT functional choice rather than a genuine CSRO effect. A required addition is a functional-sensitivity control (re-running a subset with a second functional, e.g., RPBE-D3 vs PBE-D3, or BEEF-vdW ensemble) to bound this error source before claiming a true positive.
Experimental Protocol
Minimum viable test (single composition, e.g., Co40Ni40V20):
- Generate 2 bulk/slab configurations: CSRO-random (Monte Carlo quench, high-T anneal + rapid quench in lattice MC) and CSRO-ordered (simulated anneal to low-T equilibrium with pairwise cluster expansion or MLIP-driven MC).
- Confirm CSRO contrast via Warren-Cowley parameters (target Δα ≥0.15 between states).
- Build (111), (100), (211) slabs (4-6 layers, 3x3-4x4 supercells) from each bulk state; relax with MLIP.
- Sample ≥500 adsorption sites per facet per CSRO state for CO2 and H2 (top/bridge/hollow, multiple local compositional environments) using MLIP single-point + local relaxation.
- Build AEDs (histogram + KDE) per adsorbate/facet/CSRO state.
- Feed AED-derived site-averaged energetics into a mean-field microkinetic model (CO2 hydrogenation network: RWGS, methanation, methanol synthesis pathways) to predict selectivity at 473K/20bar/H2:CO2=3:1.
- Compare selectivity distributions (bootstrap over AED sampling, N=1000 resamples) between ordered vs. disordered states; report % relative change with 95% CI.
- DFT training set for MLIP: ≥5,000-10,000 single-point energies/forces on Co-Ni-V bulk, surface slabs, and adsorbate configurations (VASP/Quantum ESPRESSO, PBE+D3 or RPBE functional), spanning composition sweep and CSRO states.
- Existing open MLIP architectures: MACE, NequIP, or Allegro (pretrained or fine-tuned).
- Reference binary alloy CSRO datasets (e.g., published Co-Ni, Ni-V short-range order DFT/MC studies) for cross-validation of Warren-Cowley computation pipeline.
- Microkinetic modeling framework: CatMAP, MKMCXX, or custom mean-field solver with CO2 hydrogenation reaction network (≥15 elementary steps).
- Adsorbate reference database: CO2, H2, H, CO, HCOO*, CH3O*, OH*, H2O binding energies on Co/Ni/V surfaces (from Catalysis-Hub / OC20-style datasets if available, for baseline sanity-checking).
- MLIP validated: energy MAE <30 meV/atom, force MAE <100 meV/Å on held-out test set.
- CSRO contrast achieved: |Δα| ≥0.15 for at least one first-shell pair correlation between ordered/disordered states.
- AED shift: mean adsorption energy shift ≥0.15 eV and/or std-dev change ≥30% for at least one adsorbate-facet combination, statistically significant (KS test p<0.01).
- Selectivity shift: ≥20% relative change in at least one product branching ratio (CO vs CH4 vs CH3OH) between CSRO states, with bootstrap 95% CI excluding zero.
- (Stretch) Experimental cross-validation confirms directionality (not necessarily magnitude) of the selectivity shift.
- AED shift <0.05 eV and std-dev change <10% across all facets/adsorbates despite confirmed CSRO contrast ≥0.15 → CSRO does not measurably affect AEDs at this composition/size.
- Selectivity shift <10% despite adequate AED shift → AED changes don't propagate through microkinetics (scaling-relation buffering).
- MLIP fails accuracy validation and cannot be improved within 2x additional training data/compute → methodology inconclusive, must halt and reassess architecture.
- CSRO states cannot be experimentally realized/quenched-in for the studied compositions (literature/thermodynamic check) → hypothesis lacks physical relevance even if computationally "true."
100
GPU hours
30d
Time to result
$1,000
Min cost
$10,000
Full cost
ROI Projection
Implementation Sketch
# Stage 1: MLIP training dft_data = generate_or_curate_dft(Co-Ni-V, bulk+slab+adsorbates, n=8000) mlip = train_MACE(dft_data, val_split=0.15) assert mlip.energy_MAE < 0.03 # eV/atom assert mlip.force_MAE < 0.10 # eV/A # Stage 2: CSRO structure generation for state in ["random", "ordered"]: bulk[state] = hybrid_MC_anneal(mlip, composition="Co40Ni40V20", target_alpha=CSRO_TARGETS[state]) verify_warren_cowley(bulk[state], shells=[1,2,3]) # Stage 3: Slab + adsorption sampling for state in ["random","ordered"]: for facet in ["111","100","211"]: slab = cut_and_relax(bulk[state], facet, mlip) sites = enumerate_sites(slab, n_min=500) for adsorbate in ["CO2","H2"]: AED[state][facet][adsorbate] = [ mlip.adsorption_energy(slab, site, adsorbate) for site in sites ] # Stage 4: Statistics + microkinetics for state in ["random","ordered"]: stats = summarize(AED[state]) # mean, std, KS-test vs other state rates = scaling_relations_to_rates(AED[state]) selectivity[state] = microkinetic_solve(rates, T=473, P=20, H2_CO2=3) delta_selectivity = relative_change(selectivity["ordered"], selectivity["random"]) report(delta_selectivity, bootstrap_CI=1000)
- Checkpoint 1 (Day 20): MLIP validation fails accuracy threshold after 2 retraining iterations → abort/redesign.
- Checkpoint 2 (Day 35): Cannot generate CSRO contrast |Δα|≥0.15 within accessible MC anneal protocols → abort or narrow hypothesis to smaller Δα claim.
- Checkpoint 3 (Day 60): AED shift <0.05 eV across all facets/adsorbates with confirmed CSRO contrast → stop, report negative result (disproof), do not proceed to microkinetics.
- Checkpoint 4 (Day 90): Microkinetic selectivity shift <10% despite adequate AED shift → stop, report partial disproof (AED-selectivity decoupling), do not proceed to experimental validation stage.
NAMED_EXPERTS: []
CLOSEST_EXISTING_WORK: []
NOVELTY_NARROWING_REQUIRED: false
SPINE_STATEMENT: This hypothesis tests whether tuning chemical short-range order in Co-Ni-V alloys, independent of nominal composition, causes a machine-learned-force-field-computed shift in CO2/H2 adsorption energy distributions large enough to change predicted CO2 hydrogenation selectivity by at least 20%.