“…For large g, which we define for all protocols as g ≡ τ a /τ I , the excitation number decays as a power law. For any g and N , our analytically solvable model (Protocol 1) always outperforms the other protocols, although all of them show scaling similar to 1/g for large g. In numerous other tests (not shown), we found that all non-powerlaw schedules, e.g., with g(t) decaying exponentially, had much worse performance for the same values of τ I , τ a , and T , in agreement with [22]. Figure 3 also shows the data that we used to extrapolate the results to larger N .…”