  Heavy-ion collisions at varying collision energies provide access to
different regions of the QCD phase diagram. In particular collisions at
intermediate energies are promising candidates to experimentally identify the
postulated first order phase transition and critical end point.
  While heavy-ion collisions at low and high collision energies are
theoretically well described by transport approaches and
hydrodynamics+transport hybrid approaches, respectively, intermediate energy
collisions remain a challenge.
  In this work, a modular hybrid approach, the SMASH-vHLLE-hybrid coupling 3+1D
viscous hydrodynamics (vHLLE) to hadronic transport (SMASH), is introduced. It
is validated and subsequently applied in Au+Au/Pb+Pb collisions between
$\sqrt{s_\mathrm{NN}}$ = 4.3 GeV and $\sqrt{s_\mathrm{NN}}$ = 200.0 GeV to
study the rapidity and transverse mass distributions of identified particles as
well as excitation functions for $\mathrm{dN}/\mathrm{d}y|_{y = 0}$ and
$\langle p_\mathrm{T} \rangle$. A good agreement with experimental measurements
is obtained, including the baryon stopping dynamics. The transition from a
Gaussian rapidity spectrum of protons at lower energies to the double-hump
structure at high energies is reproduced. The centrality and energy dependence
of charged particle $v_2$ is also described reasonably well. This work serves
as a basis for further studies, e.g. systematic investigations of different
equations of state or transport coefficients.
