This is supporting data for "Mass-Assisted Local Deconfinement in a Confined Z2 Lattice Gauge Theory". DOI: https://doi.org/10.48550/arXiv.2404.11645 The data is provided under CC BY-NC 4.0 license. This contains all the data plotted in the main text and Appendices as well as code (in the form of Python Jupyther notebooks) to reproduce these plots. The data and code are found in the file zip file "Data+Code.zip", which contains two folders names "Data" and "Code". The data is found in the folder "Data" in numpy compressed format (npz) files. Python code to load the data and plot all figures is found in the folder "Code", with one Jupyter notebook per figure. There is also a Python file called custom_funcs.py which contains a few functions shared between different notebooks. Crucially, it contains the function "get_dirp" which returns the relative path from the Code fodler to the Data folder. If you want to modify where the Data and Code folders are relative to each other you will need to modify this function accordingly. The code used to produce the data is available upon reasonable request to the first author. Content of files: Z2_grid_res_rf_N{0}.npz : data for level statistics for {0} sites and for various values of mu and h res_r: 2D array, full-spectrum value res_r_half: 2D array, in the lower half of the spectrum mul: 1D array, values of mu (along first dimension of res_r and res_r_half) hl: 1D array, values of h (along second dimension of res_r and res_r_half) N: number of gauge sites Z2_stag_hole_dyn_N60.npz : data for quenches in a system of 60 fermionic sites with two holes in the middle t_vec: 1D array, time points mul: 1D array, list of values of mu used hl: 1D array, list of values of h used J: value of J used N: number of sites res_nh: 3D array, expectation value of the fermionic number operator n for each set of parameter values (as given in mul and hl, along axis 0) for each site (along axis 1) and at all times (along axis 2) Hkh: 2D array, eigenvalues of the system (along axis 1) for each set of parameter values (as given in mul and hl, along axis 0) olaph: 2D array, overlap of eigenstates of the system (along axis 1) with the initial state for each set of parameter values (as given in mul and hl, along axis 0) Z2_stag_part_dyn_N60.npz : data for quenches in a system of 60 fermionic sites with two particles in the middle t_vec: 1D array, time points mul: 1D array, list of values of mu used hl: 1D array, list of values of h used J: value of J used N: number of sites res_np: 3D array, expectation value of the fermionic number operator n for each set of parameter values (as given in mul and hl, along axis 0) for each site (along axis 1) and at all times (along axis 2) Hkp: 2D array, eigenvalues of the system (along axis 1) for each set of parameter values (as given in mul and hl, along axis 0) olapp: 2D array, overlap of eigenstates of the system (along axis 1) with the initial state for each set of parameter values (as given in mul and hl, along axis 0) Z2_all_stats_h{0}_mu{1}_N{2}.npz : ED data in the relevant symmetry sectors for the states Psi_2, Psi_4, Phi and Psi'_2 for a system with {2} sites (at half-filling of domain walls) with J=1, h={0} and mu={1} with PBC N: number of gauge sites mu: value of mu h: value of h w: value of alpha used for the Psi'_2 state kp_lim: 1D array, indices marking the limits between the different symmetry sectors in the eigenvalue and eigenstate data Hk: 1D array, eigenvalues olap: 2D array, overlap of eigenstates of the system (along axis 0) with the various initial states (Psi_2, Psi_4, Phi and Psi'_2 along axis 1) res_H0: 1D array, expectation value of the zeroth-order Hamiltonian H0 (corresponding to the classical energy) for all eigenstates res_S: 1D array, entanglement entropy of eigenstates in the fully symmetric sector t_vec: 1D array, time points res_fid: 2D array, fidelity after a quench from various initial states (Psi_2, Psi_4, Phi and Psi'_2 along axis 0) at all times (along axis 1) Z2_all_stats_Z3_h{0}_mu{1}_N{2}.npz : ED data in the relevant symmetry sectors for the state Psi_3 for a system with {2} sites (at half-filling of domain walls) with J=1, h={0} and mu={1} with PBC N: number of gauge sites mu: value of mu h: value of h w: value of alpha used for the Psi'_2 state kp_lim: 1D array, indices marking the limits between the different symmetry sectors in the eigenvalue and eigenstate data Hk: 1D array, eigenvalues olap: 1D array, overlap of eigenstates of the system with the Psi_3 state res_H0: 1D array, expectation value of the zeroth-order Hamiltonian H0 (corresponding to the classical energy) for all eigenstates res_S: 1D array, entanglement entropy of eigenstates in the fully symmetric sector t_vec: 1D array, time points res_fid: 1D array, fidelity after a quench from Psi_3 Z2_full_stag_proj_mu{0}_N{1}.npz : Full dynamics data from Phi and Psi_4 in a system with {1} sites at the resonance mu=h={0} N: number of gauge sites mu: value of mu and h t_vec: 1D array, time points res_proj: 2D array, expectation value of the projector on the H_1 sector of the initial states over time (along axis 1) and for both initial states (Phi and Psi_4, along axis 0) res_projH0: 2D array, expectation value of the projector on the H_0 sector of the initial states over time (along axis 1) and for both initial states (Phi and Psi_4, along axis 0) Z2_stag_Heff2_spectrum_mu{0}_N{1}.npz : ED data in the relevant symmetry sectors for the states Psi_2, Psi_4 and Phi for the effective Hamiltonian up to 2nd order in a system with {1} sites (at half-filling of domain walls) with J=1, h=mu={0} with PBC N: number of gauge sites r_val: gap ratio in fully symmetric sector kp_lim: 1D array, indices marking the limits between the different symmetry sectors in the eigenvalue and eigenstate data Hk: 1D array, eigenvalues olap: 1D array, overlap of eigenstates of the system with the Psi state eig_S: 1D array, entanglement entropy of eigenstates in the fully symmetric sector t_vec: 1D array, time points res_ff: 3D array, overlap of the time-evolved state after a quench from various initial states (Psi_2, Psi_4 and Phi, along axis 0) on the states Psi_2, Psi_4 and Phi (along axis 1) at all times (along axis 2) Z2_spec_SW2_N{0}.npz: data for gap-ratio in the effective model up to second order in a system with {0} sites and various values of mu N: number of gauge sites mul: 1D array, values of mu=h considered res_r: 1D array, gap ratio in fully symmetric sector for all values of mu Z2_Heff2_proj_dyn_N{0}.npz : Dynamics data from Psi_4 and Phi for the effective Hamiltonian up to 2nd order in a system with {0} sites at the resonance mu=h=1 N: number of gauge sites mu: value of mu h: value of h t_vec: 1D array, time points res_proj: 2D array, expectation value of the projector on the H_1 sector of the initial states over time (along axis 1) and for both initial states (Psi_4 and Phi, along axis 0) res_proj2: 2D array, expectation value of the projector on H_2 neighboring states of the H_1 sector of the initial states over time (along axis 1) and for both initial states (Psi_4 and Phi, along axis 0) res_ff: 3D array, overlap of the time-evolved state after a quench from various initial states (Psi_4 and Phi, along axis 0) on the states Psi_4, Phi, Psi_2, beta and Psi'_2 (along axis 1) at all times (along axis 2) Z2_stag_dyn_neignbor_states_N{0}.npz: Dynamics data from Phi for the effective Hamiltonian up to 2nd order in a system with {0} sites at the resonance mu=h=1 N: number of gauge sites mu: value of mu h: value of h t_vec: 1D array, time points res_Krylov: 2D array, overlap of the wavefunction over time (along axis 0) with all symmetry-resolved neighbor states (along axis 1) Z2_revs_alpha_N{0}.npz: Dynamics data for quenches from the Psi'_2 state in the full Hamiltonian a system with {0} sites at the resonance for various values of mu and alpha N: number of gauge sites mul: 1D array, values of mu=h used wl: 1D array, values of alpha used t_vec: 1D array, time points res_fid: 3D array, fidelity over time (along axis 0) for quenches from the Psi'_2 with various values of mu (as given in mul, along aixs 1) and various values of alpha (given in wl, along axis 2) max_fid: 1D array, maximum revival fidelity overall all alpha for each value of mu max_w: 1D array, for each value of mu, value of alpha for which the maximum revival fidelity is obtained Z2_stag_matter_dyn_mu{0}_N{1}.npz: Data for the evolution of fermionic occupation from the Psi_2 and Psi_4 states in the full system with {1} sites (at half-filling of domain walls) with J=1, h=mu={0}, with PBC N: number of gauge sites mu: value of mu h: value of h J: value of J t_vec: 1D array, time points res_n: 3D array, expectation value of the occupation of each fermionic site (along axis 1) over time (along axis 2) for both initial states (Psi_2 and Psi_4, along axis 0) Z2_stag_matter_dyn_Psi3_mu{0}_N{1}.npz: Data for the evolution of fermionic occupation from the Psi_3 state in the full system with {1} sites (at half-filling of domain walls) with J=1, h=mu={0}, with PBC N: number of gauge sites mu: value of mu h: value of h J: value of J t_vec: 1D array, time points res_n: 2D array, expectation value of the occupation of each fermionic site (along axis 0) over time (along axis 1) after a quench from the Psi_3 state Z2_stag_matter_dyn_imb_mu{0}_N{1}.npz: Data for the evolution of fermionic occupation imbalance from the Psi_2 and Psi_4 states in the full system with {1} sites (at half-filling of domain walls) with J=1, h=mu={0}, with PBC N: number of gauge sites mu: value of mu h: value of h J: value of J t_vec: 1D array, time points res_imb: 2D array, expectation value of the fermionic occupaion imbalance over time (along axis 1) after a quench from both initial states (Psi_2 and Psi_4, along axis 0) CE: prediction of imbalance from the canonical ensemble beta: effective temperature of the initial states CEr: prediction of imbalance from the canonical ensemble restricted to the relevant symmetry sectors betar: effective temperature of the initial states among states restricted to the relevant symmetry sectors Z2_grid_res_spec_N{0}.npz: Data for the full spectrum in the most symmetric sector of the full model with {0} sites with PCB, h=-0.25712973861329 and mu=0.34641016151377546. We note that the spectrum is identical for a positive value of h. mu: value of mu h value of h N: number of gauge sites Hk: 1D array, eigenvalues HK_unf: 1D array, eigenvalues after unfolding Level_stats.npz: Non-model-specific data showing the expected distribution of unfolded level-spacing according to the Poisson distribution and to the Wigner-Dyson ensemble for the GOE. We emphasize that here we use a numerically computed Wigner-Dyson distribution for an *infinite* GOE matrix instead of the simple Wigner surmise for a 2 by 2 matrix. s: 1D array, level spacing values sampled Poisson: 1D array, P(s) according to a Poisson distribution WD: 1D array, P(s) according to the Wigner-Dyson distribution for an infinite GOE matrix Z2_stag_exact_QMBSN{0}.npz: Data for exact scars in the full model with {0} sites, J=1, mu=0.74 and h=0.3, with PBC N: number of gauge sites mu: value of mu h: value of h J: value of J Hk: 1D array, energy eigenvalues at half-filling in the most symmetric sector res_S: 1D array, entanglement entropy of eigenstates at half-filling in the most symmetric sector idx_dw: index pointing to the scarred states of each kind which live in the same domain-wall sector as the rest of the data EEs: energy of first kind of scarred states SSs: entanglement entropy of first kind of scarred states EEsb: energy of second kind of scarred states SSsb: entanglement entropy of second kind of scarred states Z2_n_detuning_mu{0}_N{1}.npz: data for dynamics in the full model with {1} sites with detuning, for a value mu={0}, after quenches from Psi_2 and Psi_4 N: number of gauge sites mu: value of mu dl: 1D array, list of values of detuning tvl: 2D array, time steps (along axis 1) for each value of detuning (given in dl, along axis 0) res_fid: 3D array, fidelity over time (along axis 2) after a quench from both initial states (Psi_2 and Psi_4, along axis 1), for various values of detuning (given in dl, along axis 0) nl: 1D array, values of n to which the detuning in dl corresponds to according to the resonance condition dl=(nl-1)/sqrt(2n) Z2_stag_detuning_N{0}.npz: Data for quenches from the Psi_4 state in the full model with N={0} sites with detuning N: number of sites mul: 1D array, values of mu used hl: 1D array, values of h used t_vec: time date res_fid: 3D array, for each mu (given in mul, along axis 0), fidelity over time (along axis 1) with either no detuning (i.e h=mu, slice 0 of axis 2) or with detuning (i.e h given in hl, slice 1 of axis 2)