  The QCD equation of state at zero baryon chemical potential is the only
element of the standard dynamical framework to describe heavy ion collisions
that can be directly determined from first principles. Continuum extrapolated
lattice QCD equations of state have been computed using 2+1 quark flavors
(up/down and strange) as well as 2+1+1 flavors to investigate the effect of
thermalized charm quarks on QCD thermodynamics. Lattice results have also
indicated the presence of new strange resonances that not only contribute to
the equation of state of QCD matter but also affect hadronic afterburners used
to model the later stages of heavy ion collisions. We investigate how these new
developments obtained from first principles calculations affect multiparticle
correlations in heavy ion collisions. We compare the commonly used equation of
state S95n-v1, which was constructed using what are now considered outdated
lattice results and hadron states, to the current state-of-the-art lattice QCD
equations of state with 2+1 and 2+1+1 flavors coupled to the most up-to-date
hadronic resonances and their decays. New hadronic resonances lead to an
enhancement in the hadronic spectra at intermediate $p_T$. Using an outdated
equation of state can directly affect the extraction of the shear viscosity to
entropy density ratio, $\eta/s$, of the quark-gluon plasma and results for
different flow observables. The effects of the QCD equation of state on
multiparticle correlations of identified particles are determined for both AuAu
$\sqrt{s_{NN}}=200$ GeV and PbPb $\sqrt{s_{NN}}=5.02$ TeV collisions. New
insights into the $v_2\{2\}$ to $v_3\{2\}$ puzzle in ultracentral collisions
are found. Flow observables of heavier particles exhibit more non-linear
behavior regardless of the assumptions about the equation of state, which may
provide a new way to constrain the temperature dependence of $\eta/s$.
