Table of Contents
Fetching ...

Bayesian Optimization for Function-Valued Responses under Min-Max Criteria

Pouya Ahadi, Reza Marzban, Ali Adibi, Kamran Paynabar

TL;DR

This paper tackles optimizing function-valued responses under a worst-case (min–max) criterion by introducing MM-FBO, which combines FPCA-based dimensionality reduction with Gaussian process surrogates for the principal component scores. The authors derive a tractable acquisition function that blends worst-case exploitation with integrated uncertainty and prove theoretical guarantees: a discretization bound connecting the continuous and discrete objectives and a consistency result ensuring convergence of the acquisition to the true min–max objective. Empirical results on synthetic benchmarks and physics-inspired case studies in metasurface design and vapor-phase infiltration show MM-FBO consistently achieves faster, more stable convergence and significantly stronger worst-case performance than standard BO baselines. The work highlights the value of explicitly modeling functional uncertainty in design optimization and provides a practical framework for robust function-valued optimization in engineering and physics-oriented applications.

Abstract

Bayesian optimization is widely used for optimizing expensive black box functions, but most existing approaches focus on scalar responses. In many scientific and engineering settings the response is functional, varying smoothly over an index such as time or wavelength, which makes classical formulations inadequate. Existing methods often minimize integrated error, which captures average performance but neglects worst case deviations. To address this limitation we propose min-max Functional Bayesian Optimization (MM-FBO), a framework that directly minimizes the maximum error across the functional domain. Functional responses are represented using functional principal component analysis, and Gaussian process surrogates are constructed for the principal component scores. Building on this representation, MM-FBO introduces an integrated uncertainty acquisition function that balances exploitation of worst case expected error with exploration across the functional domain. We provide two theoretical guarantees: a discretization bound for the worst case objective, and a consistency result showing that as the surrogate becomes accurate and uncertainty vanishes, the acquisition converges to the true min-max objective. We validate the method through experiments on synthetic benchmarks and physics inspired case studies involving electromagnetic scattering by metaphotonic devices and vapor phase infiltration. Results show that MM-FBO consistently outperforms existing baselines and highlights the importance of explicitly modeling functional uncertainty in Bayesian optimization.

Bayesian Optimization for Function-Valued Responses under Min-Max Criteria

TL;DR

This paper tackles optimizing function-valued responses under a worst-case (min–max) criterion by introducing MM-FBO, which combines FPCA-based dimensionality reduction with Gaussian process surrogates for the principal component scores. The authors derive a tractable acquisition function that blends worst-case exploitation with integrated uncertainty and prove theoretical guarantees: a discretization bound connecting the continuous and discrete objectives and a consistency result ensuring convergence of the acquisition to the true min–max objective. Empirical results on synthetic benchmarks and physics-inspired case studies in metasurface design and vapor-phase infiltration show MM-FBO consistently achieves faster, more stable convergence and significantly stronger worst-case performance than standard BO baselines. The work highlights the value of explicitly modeling functional uncertainty in design optimization and provides a practical framework for robust function-valued optimization in engineering and physics-oriented applications.

Abstract

Bayesian optimization is widely used for optimizing expensive black box functions, but most existing approaches focus on scalar responses. In many scientific and engineering settings the response is functional, varying smoothly over an index such as time or wavelength, which makes classical formulations inadequate. Existing methods often minimize integrated error, which captures average performance but neglects worst case deviations. To address this limitation we propose min-max Functional Bayesian Optimization (MM-FBO), a framework that directly minimizes the maximum error across the functional domain. Functional responses are represented using functional principal component analysis, and Gaussian process surrogates are constructed for the principal component scores. Building on this representation, MM-FBO introduces an integrated uncertainty acquisition function that balances exploitation of worst case expected error with exploration across the functional domain. We provide two theoretical guarantees: a discretization bound for the worst case objective, and a consistency result showing that as the surrogate becomes accurate and uncertainty vanishes, the acquisition converges to the true min-max objective. We validate the method through experiments on synthetic benchmarks and physics inspired case studies involving electromagnetic scattering by metaphotonic devices and vapor phase infiltration. Results show that MM-FBO consistently outperforms existing baselines and highlights the importance of explicitly modeling functional uncertainty in Bayesian optimization.

Paper Structure

This paper contains 13 sections, 2 theorems, 56 equations, 9 figures, 5 tables, 1 algorithm.

Key Result

Proposition 1

Assume $e(\theta,\lambda)$ admits a modulus of continuity $\omega$ uniformly in $\theta$ as above. Then for every $\theta\in\Theta$, If, in addition, $e(\theta,\lambda)$ is $L$ Lipschitz in $\lambda$ uniformly in $\theta$, then In the one-dimensional case with $\Lambda=[a,b]\subset\mathbb R$ and an equally spaced grid of $T$ points with spacing $\Delta=(b-a)/(T-1)$, the fill distance is $h_T=\De

Figures (9)

  • Figure 1: Representative trajectories drawn from each simulator. Each panel shows multiple random draws of the response over time to illustrate functional variability across the design space.
  • Figure 2: Regret versus iteration with median and interquartile bands.
  • Figure 3: Normalized regret trajectories showing convergence speed and stability.
  • Figure 4: Box plots of final regret distributions across replications.
  • Figure 5: AUOC distributions summarizing full regret trajectories.
  • ...and 4 more figures

Theorems & Definitions (2)

  • Proposition 1
  • Theorem 2: Consistency of the acquisition under vanishing uncertainty