try C=2 models

niusen · Feb 1, 2023 · 2965779 · 2965779
1 parent c51fbf3
commit 2965779
Show file tree

Hide file tree

Showing 15 changed files with 78 additions and 488 deletions.
diff --git a/.gitignore b/.gitignore
@@ -2,3 2,4 @@
 *.jld2
 *.mat
 *.jld
 *.jld
diff --git a/ES_CTMRG/spinful_M1/ES_Gutzwiller_y6_CTM_iterative.jl b/ES_CTMRG/spinful_M1/ES_Gutzwiller_y6_CTM_iterative.jl
@@ -15,7 15,7 @@ include("D:\\My Documents\\Code\\Julia_codes\\Tensor network\\GfPEPS_parton\\ES_
 include("D:\\My Documents\\Code\\Julia_codes\\Tensor network\\GfPEPS_parton\\ES_CTMRG\\ES_algorithms.jl")
 
 
-chi=70
 chi=20
 tol=1e-6
 CTM_ite_nums=500;
 CTM_trun_tol=1e-10;
@@ -100,7 100,7 @@ A_origin=deepcopy(A);
 
 
 y_anti_pbc=true;
-boundary_phase_y=0;
 boundary_phase_y=0.5;
 
 if y_anti_pbc
  gauge_gate1=gauge_gate(A,2,2*pi/6*boundary_phase_y);

diff --git a/ES_exact_contraction/spinful_M1/ES_Gutzwiller_y4_krylov.jl b/ES_exact_contraction/spinful_M1/ES_Gutzwiller_y4_krylov.jl
@@ -94,7 94,7 @@ A_origin=deepcopy(A);
 
 
 y_anti_pbc=true;
-boundary_phase_y=0.0;
 boundary_phase_y=0.5;
 
 if y_anti_pbc
  gauge_gate1=gauge_gate(A,2,2*pi/4*boundary_phase_y);
@@ -508,12 508,21 @@ k_phase=vcat(k_phase_set1,k_phase_set2);
 Qn=vcat(Qn_set1,Qn_set2);
 Spin=vcat(Spin_set1,Spin_set2);
 
-matwrite("ES_Gutzwiller_M1_Nv4"*".mat", Dict(
- "k_phase" => k_phase,
- "eu" => eu,
- "Qn"=>Qn,
- "Spin"=>Spin
-); compress = false)
-
 
 
 if y_anti_pbc
  matwrite("ES_Gutzwiller_M1_Nv4_APBC"*".mat", Dict(
  "k_phase" => k_phase,
  "eu" => eu,
  "Qn"=>Qn,
  "Spin"=>Spin
  ); compress = false)
 else
  matwrite("ES_Gutzwiller_M1_Nv4"*".mat", Dict(
  "k_phase" => k_phase,
  "eu" => eu,
  "Qn"=>Qn,
  "Spin"=>Spin
  ); compress = false)
 end
 
diff --git a/ES_exact_contraction/spinful_M1/ES_y4_krylov.jl b/ES_exact_contraction/spinful_M1/ES_y4_krylov.jl
@@ -93,7 93,7 @@ A_origin=deepcopy(A);
 
 
 
-y_anti_pbc=true;
 y_anti_pbc=false;
 boundary_phase_y=0.0;
 
 if y_anti_pbc