  In a viscoelastic environment, the diffusion of a particle becomes
non-Markovian due to the memory effect. An open question is to quantitatively
explain how self-propulsion particles with directional memory diffuse in such a
medium. Based on simulations and analytic theory, we address this issue with
active viscoelastic systems where an active particle is connected with multiple
semi-flexible filaments. Our Langevin dynamics simulations show that the active
cross-linker displays super- and sub-diffusive athermal motion with a
time-dependent anomalous exponent $\alpha$. In such viscoelastic feedback, the
active particle always has superdiffusion with $\alpha=3/2$ at times shorter
than the self-propulsion time ($\tau_A$). At times greater than $\tau_A$, the
subdiffusion emerges with $\alpha$ bounded between $1/2$ and $3/4$. Remarkably,
the active subdiffusion is reinforced as the active propulsion (Pe) is more
vigorous. In the high-Pe limit, the athermal fluctuation in the stiff filament
eventually leads to $\alpha=1/2$, which can be misinterpreted with the thermal
Rouse motion in a flexible chain. We demonstrate that the motion of active
particles cross-linking a network of semi-flexible filaments can be governed by
a fractional Langevin equation combined with fractional Gaussian noise and an
Ornstein-Uhlenbeck noise. We analytically derive the velocity autocorrelation
function and mean-squared displacement of the model, explaining their scaling
relations as well as the prefactors. We find that there exist the threshold Pe
($\mathrm{Pe}^*$) and cross-over times ($\tau^*$ and $\tau^\dagger$) above
which the active viscoelastic dynamics emerge on the timescales of $\tau^*
\lesssim t \lesssim \tau^\dagger$. Our study may provide a theoretical insight
into various nonequilibrium active dynamics in intracellular viscoelastic
environments.
