Dataset: Probabilistic work extraction on a classical oscillator beyond the second law

Data set for the article "Probabilistic work extraction on a classical oscillator beyond the second law" , currently published as a preprint on arXiv, and accepted in Physical Review Letters.

The preprint and SM have been uploaded for convenience.

Inside the compressed file, you will find several folders

It contains:

/Exemple trajectories/

 10 folders named Protocol 0 , Protocol 1 .... each containing the raw data of the 100 first experiments for each protocol. Each protocol was repeated 3000 time to increase the statistics for the article. We provide a representative subset.

The Python analysis code "readexperimentalfile.py" . It reads the data in the aforementioned files, and computes the initial potential, the final potential, and the work associated to each protocols. Optionnaly, it saves trajectories used create Fig3 and Fig 4.

/Figures/

Fig 1 : The pdf file for Fig1. (stochastic2.pdf) as presented in the article. Other formats can be provided on demand.

Fig 2: The pdf file for Fig2, Setup.pdf

Fig 3:  The python code used to create Fig3 , named "plotpotentiel.py". To use it, one need the parameters tuning the initial potential used in the fifth protocol 'fitdepart5.npy', the final potential 'fitfin5.npy' as well as the set of trajectories "zfin5.npy" and "zinit5.npy". All those files were created using the provided python code "readexperimentalfile.py" , and are already provided in the same folder.

Fig4: The python code used to create Fig4, named "plottrajex.py". "ztotsubset5.npy" and "lambdasubset5.npy" are provided to directly plot z(t) and lambda(t) during the protocol. "fitdepart5.npy" and 'fitfin5.npy", also created using "readexperimentalfile.py", are needed. 

Fig5: a) Everything needed to plot the first part of Fig 5. using the python file "workdpdf.py".  All the datas of the other experiments are provided, and one can also choose to plot the work probability density function for other protocols (by changing the 5 t o another number. b) "mainfig.py" creates the bottom part of Fig5.