materials-chemistry/ai produced the result/Nature Communications 2026 · v2
Neural network designs five hard metallic glasses, all confirmed by experiment
Researchers trained a neural network on 673 published hardness measurements of metallic glasses, then used it in reverse to invent new alloy recipes. Five were cast and tested, and their measured hardness matched the model's predictions.
spectrum · one line per step, placed by what the step does · bright lines used AI
Attention-enhanced variational learning for physically informed discovery of exceptionally hard multicomponent bulk metallic glasses
Nature Communications, 2026
doi:10.1038/s41467-026-73008-0 · record aix-00206 v2 · checked 2026-10-09
- AI was for
- Property prediction, Candidate generation, Simulation surrogate
- Model family
- Autoencoder, Transformer, Multilayer perceptron, Graph neural network, Gaussian process, Random forest, Gradient-boosted trees, Linear model, Clustering
- Checked by
- Experimental5 tested, 5 worked
- Code
- not reported
The finding the paper is about came from the AI.
What this research was about
Most metals are crystalline: their atoms sit in a repeating grid. Cool certain mixtures of metals fast enough and the atoms freeze in a jumbled arrangement instead, more like glass than like ordinary steel. These bulk metallic glasses can be very hard, meaning they resist being dented. Hardness is usually measured by pressing a diamond tip into the surface under a known load and measuring the dimple left behind. The trouble is choosing what to mix. An alloy can draw on dozens of elements in any proportion, so the number of possible recipes is effectively endless, and each one has to be melted and cast before anyone knows whether it worked.
The researchers wanted a way to go backwards: instead of guessing a recipe and measuring its hardness, start from the hardness wanted and get a recipe out. They gathered 673 published records of metallic glass compositions, the load used in the test and the resulting hardness, and built a model called VIBANN to learn the link between the two.
Where AI came in
The network learned to squeeze each alloy recipe into a compact numerical summary, a sixteen-number description, and from that summary to predict hardness and to rebuild the recipe. Because the compression works in both directions, the researchers could search in the space of summaries rather than the space of recipes. They sampled promising points in that space, nudged them further using the model's own sense of which direction raises hardness, then decoded them back into real compositions. A set of filters discarded candidates the model was unsure about or that sat too far from anything it had seen. This stood in for trial-and-error searching and expert intuition about which mixtures are worth trying.
Five surviving compositions, based on boron, niobium and iron, were melted and cast into thin rods. X-ray and electron measurements showed they were glassy rather than crystalline, and the measured hardness fell within the ranges the model had predicted, with one alloy reaching 2447 hardness units under a 0.5 newton load. A separate, ready-made machine-learning model of how atoms push on one another was then used to simulate the five glasses atom by atom, standing in for far costlier physics calculations, and to connect their hardness ordering to how tightly the atoms pack.
Written by AIxSci from the checked record below, to give context for readers outside the field. It is not part of the record.
The work
Technical · from the record
A neural network called VIBANN, combining a self-attention encoder with a variational information bottleneck, was trained on 673 literature records of bulk metallic glass composition, indentation load and Vickers hardness (538 training, 135 held-out test; test R about 0.943 and mean absolute error about 55.6 HV). Its latent space was sampled with Gaussian-mixture-guided MCMC and refined by gradient-based optimisation of an uncertainty-penalised hardness objective, and decoded compositions were filtered on latent support, novelty and uncertainty. Five B-Nb-Fe-based compositions were arc melted and suction cast into 2 mm rods that X-ray diffraction and TEM showed to be amorphous, with B68Nb24Fe4W4 measuring 2447 ± 44 HV at 0.5 N; measured hardness across alloys and loads lay within the model's prediction intervals. Machine-learning-potential molecular dynamics was then used to relate the hardness ordering across the five alloys to local packing, coordination and bond-angle statistics.
How AI was used
A curated literature dataset of bulk metallic glass compositions, indentation loads and Vickers hardness values was encoded as 56-dimensional simplex-normalised composition vectors with a standardised scalar load and split by composition-based clustering. An attention network with a variational information bottleneck was trained from scratch to predict standardised hardness while a decoder head reconstructed the composition from the 16-dimensional latent vector, with the bottleneck weight set by a feedback controller targeting a fixed KL per dimension, hyperparameters chosen by Optuna, and epistemic uncertainty obtained from Monte Carlo dropout passes with batch normalisation frozen. For inverse design at a fixed 0.5 N load, a three-component Gaussian mixture model was fitted to the training latents to define distributional support; mixture-aware multi-chain MCMC proposed latent points under a risk-aware utility combining predictive mean and uncertainty, selected seeds were refined by backpropagating an uncertainty-aware objective into the latent coordinate and re-decoding, and surviving candidates had to pass latent-likelihood, novelty-distance, uncertainty-cap and lower-confidence-bound gates before synthesis. The trained model was also run to predict hardness for the synthesised alloys across loads, to produce attention scores, integrated-gradient attributions and latent traversals, and separately a pre-trained SevenNet machine learning interatomic potential drove LAMMPS molecular dynamics quenches of the five compositions for structural analysis.
The shape of the work
Structural · the record, drawn
no AI
Curate and encode hardness dataset
Cleaning, filtering, normalising or labelling data already obtained.
A curated dataset of 673 bulk MG compositions with reported Vickers hardness (HV) values and indentation loads was compiled from peer-reviewed literature.where the paper describes this · verbatim
AI
Train VIBANN surrogate
Fitting model parameters, including fine-tuning an existing model. The AI stood in for statistical model.
Models were trained using the Adam optimizer for up to 1000 epochs with early stopping (patience 100)where the paper describes this · verbatim
AI
Sample candidate latents and decode compositions
Producing candidate objects that did not previously exist. The AI stood in for exhaustive search.
Sampling of the new points was carried out using Markov-chain Monte Carlo (MCMC), with the GMM density as the priorwhere the paper describes this · verbatim
AI
Refine latent seeds by gradient optimisation
Iterative search over a space. The AI stood in for expert judgement. Its result feeds back into an earlier step.
We then refined a subset of ten high-quality latent seeds using gradient-based optimization.where the paper describes this · verbatim
no AI
Gate candidates on support, novelty and uncertainty
Reducing a candidate set by filtering or ranking, in a single pass.
Candidates that passed this acceptance mask were ranked and reported for experimental validation.where the paper describes this · verbatim
no AI
Synthesise and characterise selected alloys
Physical execution, by hand or by robot.
All five compositions were fabricated by arc melting followed by suction casting, yielding 2 mm diameter rodswhere the paper describes this · verbatim
AI
Atomistic simulation of designed glasses
Numerical or physics simulation, including where a learned surrogate replaces it. The AI stood in for simulation.
molecular dynamics (MD) simulations were performed using the SevenNet machine learning interatomic potentialwhere the paper describes this · verbatim
AI
Attribute hardness to elements and latent directions
Extracting understanding from model behaviour.
Integrated Gradients (IG) was used to compute feature attributions for composition and load with respect to the hardness output.where the paper describes this · verbatim
What the record says
Technical · every part carries its own basis
+ in the paper~ our reading− not reported
How to read the quotations. A quotation shows where the paper describes something. It does not quote every value beside it: one passage locates a part of the work, and values without their own quotation are our reading of that passage.
All five synthesized compositions were produced by the model's latent-space generation and refinement; the reported alloys would not exist without it
We develop VIBANN, a variational information bottleneck-augmented attention-based neural network frameworkwhere the paper describes this · verbatim
Across all alloys and loads, the measured values follow the predicted trends closely and remain within the prediction intervals.where the paper describes this · verbatim
representative LAMMPS input files and supporting metadata, have been deposited in Zenodowhere the paper describes this · verbatim
What this paper did not report
Technical · absence is published deliberately
- CodeWhether the code is available is not stated.
- Trained model weightsWhether the trained model is available is not stated.
- ComputeThe hardware or time used is not stated.
- Version of VIBANN (variational information bottleneck-augmented attention neural network)Which version of the model was used is not stated.
- Version of VIBANN ablation with deterministic bottleneck replacing the VIBWhich version of the model was used is not stated.
- Version of Gaussian process surrogate for Bayesian optimisation baselineWhich version of the model was used is not stated.
- Version of Random forest regression baselineWhich version of the model was used is not stated.
- Version of Gradient boosting regression baselineWhich version of the model was used is not stated.
- Version of Multilayer perceptron regression baselineWhich version of the model was used is not stated.
- Version of Linear regression baselineWhich version of the model was used is not stated.
- Version of Nearest-neighbours regression baselineWhich version of the model was used is not stated.
- Version of Gaussian mixture model of the latent distributionWhich version of the model was used is not stated.
- What step 8 replacedThe paper gives no basis for what the AI stood in for.
About this article
Record aix-00206, version 2, checked by a person on 2026-10-09. The record describes the paper; it does not assess whether the paper's findings are right. The paper is published under CC-BY; quotations are at most 25 words. How we work · Report an error