Dataset for "Measurement-Induced Crossover of Quantum Jump Statistics in Postselection-Free Many-Body Dynamics"

Dataset for the Python language.

The "np" stands for the abbreviation for "numpy".

The folder "Variance" includes the two-dimensional array for the set [SFQJ, variance in the total system] for the measurement strengths \gamma=0.05, 0.1, ..., 1.

If you load one of them like

$ variance_array_8=np.load("variance_8.npy")

and print it as

$ variance_array_8[:, 0]

you can obtain SFQJ for \gamma=0.05, 0.1,..., 1 for L=8, and the same for the variance in the total system as

$ variance_array_8[:, 1]

In the folder "singlesite" inside "Variance", you can obtain the variance for a single site (e.g., for L=8) as

$ variance_array_8=np.load("variance_8.npy")[:, 0]

The folder "LiouvillianGap" includes the two-dimensional array for the Liouvillian gap.

You can access them by loading the files as

#### gamma = [0.05, 0.1, ..., 1] ####

physical_quantities_2D_8=np.load("physical_quantities_2D_8.npy")[:, 0]

physical_quantities_2D_10=np.load("physical_quantities_2D_10.npy")[:, 0]

physical_quantities_2D_12=np.load("physical_quantities_2D_12.npy")[:, 0]

#### gamma = [0.5, 0.55, ..., 1] (γ is different from the other data due to the convergence problem) #### 

physical_quantities_2D_14=np.load("physical_quantities_2D_14.npy")[:, 0]

for the Heisenberg model.

For the XX model, you should load them as

#### gamma = [0.05, 0.1, ..., 1] ####

physical_quantities_2D_10=np.load("physical_quantities_2D_10.npy")[:, 0]
#### gamma = [0.5, 0.55, ..., 1] (γ is different from the other data due to the convergence problem) ####

physical_quantities_2D_8=np.load("physical_quantities_2D_8.npy")[:, 0]

physical_quantities_2D_12=np.load("physical_quantities_2D_12.npy")[:, 0]

physical_quantities_2D_14=np.load("physical_quantities_2D_14.npy")[:, 0]

instead.

The folder "Autocorrelation" includes the two-dimensional array for the integrated autocorrelation function.

If you load one of them like

$ physical_quantities_2D_4=np.load("physical_quantities_2D_4.npy")

and print it as

$ physical_quantities_2D_4[:, 0]

you can obtain the half-chain integrated autocorrelation functions for \gamma=0.05, 0.1,..., 1 for L=4, and the same for that in the single site as

$ physical_quantities_2D_4[:, 1]