********************************************************************************
Title: Data for "Fractionalized fermionic multicriticality in anisotropic Kitaev spin-orbital liquids"
Authors: Max Fornoville and Lukas Janssen
Reference: 
********************************************************************************


********************************************************************************
Description:
********************************************************************************

The archive contains the data used to construct Figs. 1, 2 and 3 of the paper "Fractionalized fermionic multicriticality in anisotropic Kitaev spin-orbital liquids" by M. Fornoville and L. Janssen. 

All data files result from mean-field iteration and contain sublattice magnetizations (phi_Ax, phi_Ay, phi_Az) and (phi_Bx, phi_By, phi_Bz) and bond parameters (chi_x, chi_y, chi_z) for initial conditions (phi_Ax_0, phi_Ay_0, phi_Az_0), (phi_Bx_0, phi_By_0, phi_Bz_0), (chi_x_0, chi_y_0, chi_z_0), obtained for linear system size 2N at temperature T and couplings K, J_xy and J_z. The amount of iterations is given in the header.

The data is structured in 23 columns: N, T, K, J_xy, J_z, phi_Ax, phi_Ay, phi_Az, phi_Bx, phi_By, phi_Bz, chi_x, chi_y, chi_z, phi_Ax_0, phi_Ay_0, phi_Az_0, phi_Bx_0, phi_By_0, phi_Bz_0, chi_x_0, chi_y_0, chi_z_0.

The data is organized in three folders, corresponding to the three figures.


********************************************************************************
Fig_1:
********************************************************************************

The folder contains three subfolders, corresponding to the different transitions. Each of these folders contain data files that were used to extrapolate the mean-field transition points:

"z2" contains the data for the Z2 symmetry-breaking transition, with fixed couplings K, J_xy and varying J_z, obtained at T=0 and fixed linear system size 2N=48 at 10,000 iterations.

"u1" contains the data for the U(1) symmetry-breaking transition, with fixed couplings K, J_z and varying J_xy, obtained at T=0 and fixed linear system size 2N=48 at 10,000 iterations.

"su2" contains the data for the SU(2) symmetry-breaking transition, with fixed coupling K and varying J_xy = J_z, obtained at T=0 and fixed linear system size 2N=48 at 10,000 iterations.


********************************************************************************
Fig_2:
********************************************************************************

The folder contains four subfolders, corresponding to the Figures (a)-(d):

"z2" contains data for Fig. 2(a), with fixed coupling K=1, J_xy=0.2 and varying J_z between 0 and 2, obtained at T=0 and fixed linear system size 2N=48 at 2,500 iterations.

"u1" contains data for Fig. 2(b), with fixed coupling K=1bi" J_z=0.2 and varying J_xy between 0 and 2, obtained at T=0 and fixed linear system size 2N=48 at 2,500 iterations.

"fo" contains data for Fig. 2(c), with fixed coupling K=1 and varying J_xy between 0 and 2 and J_z=1.4-J_xy, obtained at T=0 and fixed linear system size 2N=48. The data for J_xy in [0, 0.6978] and [0.7022, 1.4] was obtained at 2,500 iterations, for J_xy in [0.698, 0.702] 10,000 were used to ensure convergence of the iteration near the first-order transition.

"bi" contains data for Fig. 2(d), with fixed coupling K=1 and varying J_xy between 0 and 2 and J_z=1.4-J_xy, obtained at T=0 and fixed linear system size 2N=48. The data for J_xy in [0, 0.6181] and [0.6225, 1.4] was obtained at 2,500 iterations, for J_xy in [0.6183, 0.6223] 10,000 were used to ensure convergence of the iteration near the bicritical transition.


********************************************************************************
Fig_3:
********************************************************************************

The folder contains mean-field data at (J_xy,J_z)=(0.5,0.2), (J_xy,J_z)=(0.6,0.2), (J_xy,J_z)=(0.2,0.6), obtained at coupling K=1, T=0 and linear system size 2N=48 at 2,500 iterations. From this data, the single-particle spectra in Fig. 3 can be calculated.
