How to pick a crystal thickness in 50 simulations ✛ ✎

Original: "Scalable Bayesian Optimization of Composite Functions for Image-Based Inverse Problems in Materials Characterization"
arXiv:2609.02126v1 · 2026-09-02 · arXiv non-exclusive · independent analysis · 2 min
The new SBOCF method reduces 24,649 diffraction pattern pixels to 11 numbers and finds sample thickness and tilt with errors up to 290 times smaller.

A translator at the negotiating table doesn't listen to every word; they grasp the meaning of whole phrases. Electron microscopy poses a similar challenge: you need to figure out the crystal's thickness and its tilt in the beam, and all you have is a diffraction pattern, a mosaic of thousands of bright and dark spots. Comparing every pixel of that pattern with a simulation is too costly, so researchers proposed listening not to individual words but to the meaning of fragments.

SBOCF cuts the image into nine squares, like a page into paragraphs, and for each it sums the brightness. Instead of 24,649 pixels, only 11 numbers remain: nine sums plus two correction coefficients that capture residual discrepancy. These numbers are modeled by independent Gaussian processes, cheap stand-ins for expensive simulations. Bayesian optimization, like a seasoned negotiator, picks the next thickness and tilt pair that promises the greatest improvement.

On a synthetic SrTiO3 sample with a thickness of 380 Å and a tilt of 1.5 mrad, the method reduced the median final error by a factor of 290 compared with standard Bayesian optimization. For a thin sample of 100 Å, the gain was 47 times; for 200 Å, 37 times.

Previously, for each new material you had to either run thousands of simulations or train a neural network on a large library. SBOCF requires no training set: it uses only the structure of the objective function and a budget of 50 simulations. On experimental data, the estimated thickness and tilt agree with previously published values, and the time to choose the next point is about 0.2 minutes per iteration.

A translator who can hear the meaning of whole phrases works faster and more accurately. Likewise, SBOCF reduces the problem to 11 numbers, finds the crystal parameters in 50 simulations, and improves the subsequent reconstruction of the atomic structure. The method can already be integrated into electron microscopy workflows, and in the future it could be applied to X-ray diffraction or neutron scattering.

The number of modeled outputs dropped from 24,649 to 11, a reduction of more than 2000 times, while the final error decreased by up to 290 times. If each simulation took an hour, a full grid search over a million points would take 114 years, but SBOCF does it in 50 hours.

f(x) = \sum_{k=1}^{M}\sum_{\ell=1}^{M} (y_{k,\ell}(x) - y^{true}_{k,\ell})^2
Objective function to be minimized for parameter estimation
f(x) = \epsilon(x) f_P(x) + \delta(x)
Relation between pixel-level and patch-level errors
Concepts
single pixel imaging image reconstruction x ray imaging optimization synthetic images Stochastic modeling information criterion image classification optimization problem 3d reconstruction signal decomposition bayesian hierarchical model least squares fitting photographic body analysis