Quantum reservoir computing with classical and nonclassical states in an integrated optical circuit

The files contain data essential to reproduce figures 2-6. A detailed list of file explanations is given below: PP - positive-P Figure 2 cat0PP - photons statistics of cat of alpha=2 and phi = 0 - PP calcualtion (panel a) cat0_expected - photons statistics of cat of alpha=2 and phi = 0 - analytic (panel a) cat1PP - photons statistics of cat of alpha=2 and phi = 1 - PP calcualtion (panel b) cat1_expected - photons statistics of cat of alpha=2 and phi = 1 - analytic (panel b) cat0.5PP - photons statistics of cat of alpha=2 and phi = 0.5 - PP calcualtion (panel c) cat0.5_expected - photons statistics of cat of alpha=2 and phi = 0.5 - analytic (panel c) cohPP - photon statistics of cat of alpha=2 - PP calculation (panel d) coh_expected - photons statistics of cat of alpha=2 - analytic (panel d) Q, W - meshgrids of X1 and X2 quadrature to plot Wigner function (panels e-f) W_cat - cat state Wigner function (panel e) WI, WII, WIII, WIV - partial Wigner functions of cat state components (panel e) W_cat4 - 4-th order cat state Wigner function (panel f) Figure 3 couplings - super-gaussian couplings (panel a) W1_before, W1_after - Wigner functions of mode 1 before and after passing through the interferometer (panel a) W2_before, W2_after - Wigner functions of mode 2 before and after passing through the interferometer (panel a) G_av - average intermode cross-correlation function dynamics (panel b) N_av[0] - n1 , N_av[1] - n2 - average occupation dynamics for the first and second mode Q,P - X1 and X2 meshgrids to calculate Wigner functions (panel b) Figure 4 N_av - PP calculated average occupation <n> dynamics (panel a) N_master_4 - ME calculated average occupation <n> dynamics with 4 cutoff (panel a) N_master_8 - ME calculated average occupation <n> dynamics with 8 cutoff (panel a) G_av - PP calculated g2 dynamics (panel b) G_master_4 - ME calculated g2 dynamics with 4 cutoff (panel b) G_master_8 - ME calculated g2 dynamics with 8 cutoff (panel b) cutofs - array of Fock space cutoffs (panel c) Gerr - error corresponding to those cutoffs (panel c) Figure 5 couplings - super-gaussian couplings (panel a,c) N_av - PP calculated average occupation <n> dynamics (panel b) G_av - PP calculated g2 dynamics (panel d) Q,P - X1 and X2 meshgrids to calculate Wigner functions (panel e) W1_after, W2_after, W3_after,W4_after - Wigner functions of modes 1,2,3,4 after passing through the interferometer (panel e) Figure 6 training - training dataset (x and y features) testing - testing dataset (x and y features) Ntr_coh - average occupations for training set examples for coherent state input Nte_coh - average occupations for testing set examples for coherent state input Ndec_coh - average occupations for a uniform grid of (x,y) for coherent state input - decision boundary calculation n_coh - decision boundary for coherent state input using average occupations Ntr - average occupations for training set examples for cat state input Ntr - average occupations for testing set examples for cat state input Ndec- average occupations for a uniform grid of (x,y) for cat state input - decision boundary calculation n - decision boundary for cat state input using cross-correlation functions Gtr - average cross-correlation functions for training set examples for cat state input Gte - average cross-correlation functions for testing set examples for cat state input Gdec- average cross-correlation functions for a uniform grid of (x,y) for cat state input - decision boundary calculation g - decision boundary for cat state input using average cross-correlation functions