We investigate the use of normalizing flow (NF) models as flexible priors in Bayesian inference via Markov chain Monte Carlo (MCMC) sampling for iterative Bayesian calibration. Trained on posteriors from previous analyses, these models can be used as informative priors that capture nontrivial distributions and correlations in subsequent inference tasks. We compare different training strategies and loss functions, finding that training based on Kullback-Leibler divergence and unsupervised learning consistently yield the most accurate reproductions of reference distributions. We apply such a sequential Bayesian workflow to a high-energy nuclear physics problem; MCMC with NF-based priors reproduces the results of one-shot joint inference well, provided the target distributions are unimodal. In cases with pronounced multi-modality or dataset tension, distortions may arise, underscoring the need for caution in multi-stage Bayesian inference. A comparison between the pocoMC MCMC sampler and the standard emcee sampler further demonstrates the importance of advanced and robust algorithms for exploring the posterior space. Overall, our results establish NF-based priors as a practical and efficient tool for sequential Bayesian inference in high-dimensional parameter spaces.

Close

PrevNext

View figure in article

• QCD phenomenology
• Relativistic heavy-ion collisions
• Bayesian methods
• Machine learning