Bayesian analysis of a (3+1)D hybrid approach with initial conditions from hadronic transport

This study aims to apply statistical learning, specifically Bayesian inference, to the (3+1)D SMASH-vHLLE-hybrid model using initial conditions generated by the SMASH transport code itself, with the objective of constraining model parameters and gaining deeper insight on the temperature and baryochemical potential dependence of both the shear and the bulk viscosity. This study is performed in the hybrid approach SMASH-vHLLE, composed of the hadronic transport approach SMASH and the (3+1)D viscous hydrodynamic code vHLLE. A Bayesian framework is employed, utilizing Markov Chain Monte Carlo (MCMC) sampling to explore the parameter space. The analysis compares model predictions against experimental observables, including particle yields, momentum and flow coefficients both at midrapidity as well as in forward and backward direction. We find that the SMASH-vHLLE-hybrid framework, using hadronic initial conditions for Au+Au collisions at different beam energies, can reproduce a variety of experimental observables at midrapidity and forward/backward rapidities. Notably, the preferred posterior distribution suggests a near-vanishing specific shear viscosity in the high-temperature QGP phase, combined with moderate-to-large bulk viscosity around the phase transition region, although the constraints on baryochemical potential dependence are weak. Our findings reveal that a hadronic initial condition constrains the evolution more strictly at intermediate energies, making parameters such as the hydrodynamic onset time highly sensitive. Intriguingly, the extracted shear viscosity differs substantially from previous Bayesian analyses, motivating further systematic studies with higher-statistics data sets and refined modeling assumptions.