  The holographic phenomena of pole skipping have been studied in the presence
of scalar-Gauss-Bonnet interaction in the four-dimensional Anti-de
Sitter-Schwarzchild black hole background. Pole skipping points are special
points in phase space where the bulk linearised differential equations have
multiple ingoing solutions. Those special points are claimed to be connected to
chaos. In this paper, we initiated a novel study on understanding the response
of those special points under the application of external sources. The source
is identified with the holographic dual operator of the bulk scalar field with
its non-normalizable solutions. We analyze in detail the dynamics of pole
skipping points in both sound and shear channels, considering linear
perturbation in bulk. In the perturbative regime, characteristic parameters for
chaos, namely Lyapunov exponent and butterfly velocity, remain unchanged.
However, the diffusion coefficient has evolved non-trivially under the external
source.
