  We study the impact of pions in simulations of neutron star mergers and
explore the impact on gravitational-wave observables. We model charged and
neutral pions as a non-interacting Boson gas with a chosen, constant effective
mass. We add the contributions of pions, which can occur as a condensate or as
a thermal population, to existing temperature and composition dependent
equations of state. Compared to the models without pions, the presence of a
pion condensate decreases the characteristic properties of cold, non-rotating
neutron stars such as the maximum mass, the radius and the tidal deformability.
We conduct relativistic hydrodynamical simulations of neutron star mergers for
these modified equations of state models and compare to the original models,
which ignore pions. Generally, the inclusion of pions leads to a softening of
the equation of state, which is more pronounced for smaller effective pion
masses. We find a shift of the dominant postmerger gravitational-wave frequency
by up to 150~Hz to higher frequencies and a reduction of the threshold binary
mass for prompt black-hole formation by up to 0.07~$M_\odot$. We evaluate
empirical relations between the threshold mass or the dominant postmerger
gravitational-wave frequency and stellar parameters of nonrotating neutron
stars. These relations are constructed to extract these stellar properties from
merger observations and are built based on large sets of equation of state
models which do not include pions. Comparing to our calculations with pions, we
find that these empirical relations remain valid to good accuracy, which
justifies their use although they neglect a possible impact of pions. We also
address the mass ejection by neutron star mergers and observe a moderate
enhancement of the ejecta mass by a few ten per cent. (abridged)
