无需代理模型,实现中子反射率实时贝叶斯反演
Towards real-time surrogate-free Bayesian inversion for neutron reflectometry
- 直接计算反射率梯度,跳过传统代理模型
- 哈密顿蒙特卡洛使采样效率大幅提升,速度比MCMC快数倍
- 支持秒级不确定性量化,适合复杂多层器件分析
中子反射率(NR)是诸多科学领域的重要技术。尽管正向反射模型已明确,但从数据反推样品物理性质需解逆问题。随着束线科学家在快速动力学和复杂结构探测中广泛应用NR,不确定性日益显著。现有不确定性量化方法如马尔可夫链蒙特卡洛(MCMC)样本效率低、收敛慢。虽有机器学习代理模型被提出,但会丢失物理规律。本文提出一种快速、无代理的贝叶斯反演方法:首次实现反射率精确梯度计算,支持高性能梯度推理,哈密顿蒙特卡洛相比MCMC显著提升采样效率;变分推断可在秒级完成近似不确定性量化。在厚氧化硅石英膜上达到当前最优性能,在有机LED多层器件高复杂度场景中表现稳健。此外,提供开源的Python反射率核函数库。
原文摘要 · Abstract (English)
Neutron reflectometry (NR) is a key enabling technology for many areas of scientific development. Although the forward reflectivity model is well-known, inferring the physical properties of a sample from NR data requires the solution of an inverse problem. Increasingly, beamline scientists are using NR in fast kinetic configurations and probing highly-complex structures and interfaces, introducing significant uncertainty. Existing uncertainty quantification (UQ) approaches in NR, such as Markov-Chain Monte-Carlo (MCMC), suffer from poor sample efficiency and slow convergence times. Recently, surrogate machine learning models have been proposed as an alternative. However, physical intuition is lost when replacing governing equations with fast surrogates. Instead, we propose a rapid, surrogate-free Bayesian inversion for NR. Our approach offers a step-change in inference speed and efficiency. For the first time in NR, exact gradients through the reflectivity are computed, enabling highly performant gradient-based inference schemes: Hamiltonian Monte-Carlo offers significant advances in sample efficiency compared to MCMC. Variational inference enables approximate UQ on the order of seconds rather than hours. We demonstrate state-of-the-art performance on a thick oxide quartz film, and robust co-fitting performance in the high complexity regime of organic LED multilayer devices. Additionally, we provide an open-source library of reflectometry kernels in the python language.
Thank you to arXiv for use of its open access interoperability. PaperDance 不是 arXiv 官方产品;中文卡片由大模型生成,请以原文为准。