
where \(C\)  denotes a deviation from the analytical –exactly adiabatic and harmonic– estimations above, to be computed numerically. For our setup, we find numerically \(C\simeq 2\) , which confirms the result first obtained in Ref. [1]}; this correction will be included in all figures below for the kinetic+trapped mechanism.
