Abstract: Heavy-ion collisions, such as Pb-Pb or p-Pb, produce extreme conditions in temperature and density that make the hadronic matter transition to a new state, called quark-gluon plasma (QGP). Simulations of heavy-ion collisions provide a way to improve our understanding of the QGP's properties. These simulations are composed of a hybrid description that results in final observables in agreement with accelerators like LHC and RHIC. However, recent works pointed out that these hydrodynamic simulations can display acausal behavior during the evolution in certain regions, indicating a deviation from a faithful representation of the underlying QCD dynamics. To pursue a better understanding of this problem and its consequences, this work simulated two different collision systems, Pb-Pb and p-Pb at root sNN = 5.02 TeV. In this context, our results show that causality violation, even though always present, typically occurs on a small part of the system, quantified by the total energy fraction residing in the acausal region. In addition, the acausal behavior can be reduced with changes in the prehydrodynamic factors and the definition of the bulk-viscous relaxation time. Since these aspects are fairly arbitrary in current simulation models, without solid guidance from the underlying theory, it is reasonable to use the disturbing presence of acausal behavior in current simulations to guide improvements towards more realistic modeling. While this work does not solve the acausality problem, it sheds more light on this issue and also proposes a way to solve this problem in simulations of heavy-ion collisions.

Abstract: Background: The mercury isotopes around N = 104 are a well-known example of nuclei exhibiting shape coex-istence. Mixing of configurations can be studied by measuring the monopole strength rho^2(E0), however, currently the experimental information is scarce and lacks precision, especially for the I^pi -> I^pi (I not = 0) transitions. Purpose: The goals of this study were to increase the precision of the known branching ratios and internal conversion coefficients, to increase the amount of available information regarding excited states in 182,184,186Hg, and to interpret the results in the framework of shape coexistence using different models. Method: The low-energy structures in 182,184,186Hg were populated in the & beta; decay of 182,184,186Tl, produced at ISOLDE, CERN and purified by laser ionization and mass separation. The & gamma;-ray and internal conversion electron events were detected by five germanium clover detectors and a segmented silicon detector, respectively, and correlated in time to build decay schemes.Results: In total, 193, 178, and 156 transitions, including 144, 140, and 108 observed for the first time in a & beta;-decay experiment, were assigned to 182,184,186Hg, respectively. Internal conversion coefficients were determined for 23 transitions, out of which 12 had an E0 component. Extracted branching ratios allowed the sign of the interference term in 182Hg as well as & rho;2(E 0; 0+2 & RARR; 0+1 ) and B(E2; 0+2 & RARR; 2+1 ) in 184Hg to be determined. By means of electron-electron coincidences, the 0+3 state was identified in 184Hg. The experimental results were qualitatively reproduced by five theoretical approaches, the interacting boson model with configuration mixing with two different parametrizations, the general Bohr Hamiltonian, the beyond mean-field model, and the symmetry-conserving configuration-mixing model. However, a quantitative description is lacking. Conclusions: The presence of shape coexistence in neutron-deficient mercury isotopes was confirmed and evidence for the phenomenon existing at higher energies was found. The new experimental results provide important spectroscopic input for future Coulomb excitation studies.