  Gravitational wave (GW) astronomy has consolidated its role as a new
observational window to reveal the properties of compact binaries in the
Universe. In particular, the discovery of the first binary neutron star
coalescence, GW170817, led to a number of scientific breakthroughs as the
possibility to place constraints on the equation of state of cold matter at
supranuclear densities. These constraints and all scientific results based on
them require accurate models describing the GW signal to extract the source
properties from the measured signal.
  In this article, we study potential systematic biases during the extraction
of source parameters using different descriptions for both, the point-particle
dynamics and tidal effects. We find that for the considered cases the mass and
spin recovery show almost no systematic bias with respect to the chosen
waveform model. However, the extracted tidal effects can be strongly biased,
where we find generally that Post-Newtonian approximants predict neutron stars
with larger deformability and radii than numerical relativity tuned models.
Noteworthy, an increase in the Post-Newtonian order in the tidal phasing does
not lead to a monotonic change in the estimated properties.
  We find that for a signal with strength similar to GW170817, but observed
with design sensitivity, the estimated tidal parameters can differ by more than
a factor of two depending on the employed tidal description of the waveform
approximant. This shows the current need for the development of better waveform
models to extract reliably the source properties from upcoming GW detections.
