We caution that for comparison with dynamo simulations in a collapsing system,
given the limited numerical resolution,
it is very challenging to numerically resolve the entire turbulent cascade down to the smallest turbulent eddies.
As the fastest dynamo rate is determined by the turnover rate of the smallest eddies,
the dynamo growth is affected and
the kinematic dynamo in simulations is prolonged
[1]}
(also see Fig. REF ).
In addition, here we adopt a Kolmogorov turbulent cascade by assuming the turbulent cascade is faster than the
gravitational contraction.
In the central collapsing region, when the local contraction becomes faster than the turbulent cascade,
additional turbulence can be driven with the gravitational energy converted to turbulent energy,
which can have a different scaling from the Kolmogorov one
[2]}, [3]}, [4]}.
