Cumulants of the chiral order parameter in a nonequilibrium Bjorken expansion with a QCD critical point

To understand experimentally obtained net-proton number cumulants in the search for the QCD critical point, we study a dynamical model based on an effective quark-meson Lagrangian with chiral symmetry. We investigate the evolution of the expanding medium created in a heavy-ion collision using a spatially homogeneous fluid and a time-dependent order parameter, the sigma field evolved by a Langevin equation. We extract cumulants of the sigma field along a parametrized freeze-out curve and match the obtained freeze-out points to corresponding beam energies. These cumulants can be related to cumulants of the net-proton number through the sigma-proton coupling to provide a qualitative comparison to experimental data from STAR’s beam energy scan program. We demonstrate that the presence of the spinodal or mixed phase region around the first-order chiral phase transition allows for a wide interval of cumulants at the lowest beam energies.