  The majority of Monte-Carlo (MC) simulation campaigns for future $e^+e^-$
colliders has so far been based on the leading-order (LO) matrix elements
provided by Whizard 1.95, followed by parton shower and hadronization in
Pythia6, using the tune of the OPAL experiment at LEP. In this contribution, we
test and develop the interface between Whizard3 and Pythia8. As a first step,
we simulate the $e^+e^-\to q\bar{q}$ process with LO matrix elements, and
compare three tunes in Pythia8: the standard Pythia8 tune, the OPAL tune and
the ALEPH tune. At stable-hadron level, predictions of charged and neutral
hadron multiplicities of these tunes are compared to LEP data, since they are
strongly relevant to the performance of particle flow algorithms.
  The events are used to perform a full detector simulation and reconstruction
of the International Large Detector concept (ILD) as an example for a
particle-flow-optimised detector. At reconstruction level, a comparison of the
jet energy resolution in these tunes is presented. We found good agreement with
previous results that were simulated by Whizard1+Pythia6. In addition, the
preliminary next-to-leading order (NLO) results are also presented. This modern
MC simulation chain, with matched NLO matrix elements in the future, should be
introduced to ILC or other future $e^+e^-$ colliders.
