Work distribution for the driven harmonic oscillator with time-dependent strength: exact solution and slow driving

We study the work distribution of a single particle moving in a harmonic potential with time-dependent strength. This simple system has a non-Gaussian work distribution with exponential tails. The time evolution of the corresponding moment generating function is given by two coupled ordinary differential equations that are solved numerically. Based on this result we study the behavior of the work distribution in the limit of slow but finite driving and show that it approaches a Gaussian distribution arbitrarily well.