materials-chemistry/ai produced the result/arXiv 2026 · v2
Robot laser-heats thin films while software picks each next heating condition
Researchers built a loop that heats patches of thin film with a laser, reads the resulting crystal structure by X-ray, and lets software decide the next condition to try. Machine learning identified the phases and chose the experiments.
spectrum · one line per step, placed by what the step does · bright lines used AI
Autonomous Materials Exploration by Integrating Automated Phase Identification and AI-Assisted Human Reasoning
arXiv, 2026
doi:10.48550/arxiv.2601.08185 · record aix-00139 v2 · checked 2026-10-09
- AI was for
- Experimental design, Classification, Property prediction
- Model family
- Gaussian process, Clustering, Probabilistic graphical model
- Checked by
- Replication
- Code
- not reported
The finding the paper is about came from the AI.
What this research was about
Many useful materials can exist in more than one crystalline arrangement, or phase. The same atoms, packed differently, can behave quite differently. Which phase you get depends not only on the recipe but on how the material is processed: how hot it gets and for how long. Charting that is slow work. Each combination of composition, temperature and heating time is a separate experiment, and the number of combinations grows quickly. Reading the result is slow too. The standard tool is X-ray diffraction, which fires X-rays at a sample and records the pattern they scatter into. Turning that pattern back into a list of phases normally takes a trained eye.
The researchers set out to run this exploration as a closed loop. Thin films were laid down in a non-crystalline, glassy state, then heated in stripes by a scanning laser, each stripe carrying a spread of temperatures along its length. Synchrotron X-rays mapped each stripe. The loop then chose where to point the laser next. They ran it on bismuth oxide, tin oxide and a bismuth-titanium-oxide film whose composition varied across the wafer, over 120 heating steps for the last of these, and compared the resulting maps with earlier work that had sampled the space exhaustively.
Where AI came in
Software did two jobs. First, reading the X-ray data: a factorisation method broke each stripe's many patterns down into a handful of representative ones, a background model let the glassy, non-crystalline ones be set aside, and an algorithm called CrystalShift matched the rest against candidate crystal structures from a database, returning probabilities rather than a single answer. That matching stood in for the expert judgement usually needed to label diffraction patterns. Second, choosing experiments: a Gaussian process, a statistical model that predicts a quantity and says how unsure it is, was fitted to the phase amounts seen so far across composition, temperature and heating time. An acquisition function then picked the condition expected to be most informative, in place of sampling the space exhaustively.
A human stayed in the loop. Experts looked at the emerging maps mid-run and added or removed phases from the list the software was aiming for, which changed what the selection step was optimising. The reported processing maps are products of this machinery: the phase labels come from the diffraction analysis, and the filled-in regions between measurements from the Gaussian process that also chose the measurements.
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
SARA-H couples robotic lateral-gradient laser spike annealing of thin films with on-the-fly X-ray diffraction analysis: non-negative matrix factorization and the CrystalShift probabilistic labeling algorithm turn each annealed stripe into expected phase fractions, which update a Gaussian process model of the processing space, and a logEI acquisition function selects the next anneal condition. A human expert can add or remove candidate phases during a campaign, which changes the objective the active learning cycle optimises. The framework was run on Bi2O3, SnOx and a Bi-Ti-O composition spread; the Bi-Ti-O campaign is reported over 120 iterations divided into 45 exploration, 25 refined exploration and 50 targeted exploitation anneals. The resulting processing phase diagrams map extensive domains of delta-Bi2O3 and Bi2Ti2O7 and regions of Bi-substituted anatase, which the authors read as evidence that Bi doping inhibits the anatase-to-rutile transformation, and no sillenite Bi12TiO20 was observed.
How AI was used
Each iteration of the closed loop begins with an lg-LSA anneal of an amorphous thin film at a condition expressed as laser power and scan velocity, converted to peak temperature and dwell time through a thermoreflectance calibration. Spatially resolved synchrotron diffraction patterns from across the stripe are integrated with PyFAI, then reduced by a non-negative matrix factorization based on an extreme-ray finding algorithm (rank k=4), whose bases are fitted with a kernel ridge regression background model so that amorphous bases can be discarded. The retained bases are passed to CrystalShift, which performs a best-first tree search with lattice refinement and Laplace-approximated Bayesian model comparison over a user-supplied candidate phase set, with a temperature-scaling parameter calibrated on 10,000 synthetic XRD patterns; the top candidates per basis are combined with the NMF activation coefficients into expected phase activations as a function of temperature. These activations, with measurement and temperature uncertainties propagated, train a Gaussian process regression model with a Matérn-5/2 kernel over composition, log dwell time and peak temperature, which serves as the surrogate for Bayesian optimisation. A logEI acquisition function, optimised with BFGS and optionally summed over the temperature range of a stripe, selects the next anneal; the targeted phase is either held fixed (target sampling) or cycled round-robin over the candidate pool (cycle sampling). Human experts inspect the evolving phase fields mid-campaign and add or remove candidate phases, shifting the objective between exploration and exploitation. Sampling strategies were also compared on a synthetic model built from GPs fitted to hand-labeled Bi-Ti-O data, with random and uncertainty sampling as baselines.
The shape of the work
Structural · the record, drawn
no AI
Deposit amorphous thin-film libraries
Physical execution, by hand or by robot.
Thin films were then deposited on the wafers using radio frequency (RF) magnetron sputteringwhere the paper describes this · verbatim
no AI
Laser spike anneal a stripe at the chosen condition
Physical execution, by hand or by robot.
For each anneal, the stage was first accelerated to the desired velocity corresponding to the desired dwell timewhere the paper describes this · verbatim
no AI
Map the stripe with synchrotron X-ray diffraction
Obtaining raw data, whether by measurement, download or retrieval.
A total of 151 two-dimensional X-ray diffraction images were collected for each anneal stripewhere the paper describes this · verbatim
AI
Factorize the XRD map into representative bases
Encoding data into features, descriptors, embeddings or graphs.
nonnegative matrix factorization (NMF) was performed to reduce the dimensionality of the datawhere the paper describes this · verbatim
AI
Label phases probabilistically with CrystalShift
Running a trained model over new data to predict, classify or score. The AI stood in for expert judgement.
CrystalShift, a rapid probabilistic multiphase labeling algorithm was employed to label phases and estimate the phase fractionswhere the paper describes this · verbatim
AI
Fit the Gaussian process phase map
Fitting model parameters, including fine-tuning an existing model. The AI stood in for physical experiment.
The expected phase activation of the annealed stripes was used to create and continuously update a Gaussian process (GP) modelwhere the paper describes this · verbatim
AI
Choose the next anneal by maximizing logEI
Iterative search over a space. The AI stood in for exhaustive search. Its result feeds back into an earlier step.
We optimized LogEI with the BFGS algorithm to find the optimal 𝐱.where the paper describes this · verbatim
no AI
Expert reviews phase fields and edits the target pool
Extracting understanding from model behaviour. Its result feeds back into an earlier step.
The human can examine the data during the run and decide to add or remove candidate phase(s)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.
The reported processing phase diagrams are produced by the autonomous loop itself: phase labels come from the NMF/CrystalShift analysis and the phase fields from the Gaussian process model that also chose which anneals were performed.
We utilized an AL scheme based on Bayesian optimization, which includes a surrogate model and an acquisition function for experimental decision making.where the paper describes this · verbatim
A comparison with the corresponding figures from the exhaustive sampling in Ref. demonstrates excellent agreement.where the paper describes this · verbatim
The calculations are performed on an HPC cluster to parallelize the workload.where 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.
- DataWhether the data are available is not stated.
- How many were testedThe paper gives no count of what was tested.
- Version of Gaussian process surrogate of expected phase activation (Matérn-5/2 kernel)Which version of the model was used is not stated.
- Version of CrystalShift probabilistic multiphase labeling algorithmWhich version of the model was used is not stated.
- Version of Nonnegative matrix factorization with extreme-ray basis selectionWhich version of the model was used is not stated.
- Version of Kernel ridge regression background model (RBF kernel) for amorphous-basis screeningWhich version of the model was used is not stated.
- Version of Gaussian process interpolation of XRF composition map (RBF kernel)Which version of the model was used is not stated.
- What step 4 replacedThe paper gives no basis for what the AI stood in for.
About this article
Record aix-00139, 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