[multiorbital] fix move_move trace ratio

c++/triqs_ctseg/moves/move_segment.cpp

@@ -71,6 +71,13 @@ namespace moves {
     double ln_trace_ratio =
        (flipped ? -1 : 1) * (wdata.mu(dest_color) - wdata.mu(origin_color)) * double(origin_segment.length());
 
+    for (auto const &[c, slist] : itertools::enumerate(config.seglists)) {
+      if (c != dest_color && c != origin_color) {
+        ln_trace_ratio += -wdata.U(dest_color, c) * overlap(slist, origin_segment) * (flipped ? -1 : 1);
+        ln_trace_ratio -= -wdata.U(origin_color, c) * overlap(slist, origin_segment) * (flipped ? -1 : 1);
+      }
+    }
+
     if (wdata.has_Dt) {
       auto tau_c    = origin_segment.tau_c;
       auto tau_cdag = origin_segment.tau_cdag;
@@ -80,6 +87,7 @@ namespace moves {
         ln_trace_ratio += K_overlap(slist, tau_c, tau_cdag, wdata.K, dest_color, c);
         ln_trace_ratio -= K_overlap(slist, tau_c, tau_cdag, wdata.K, origin_color, c);
       }
+      // Correct double counting
       auto Kl = wdata.K(double(origin_segment.length())); // matrix
       ln_trace_ratio -= real(Kl(origin_color, origin_color));
       ln_trace_ratio -= real(Kl(dest_color, dest_color));

paper/paper.md

@@ -68,7 +68,7 @@ a metal at a doping-induced quantum critical point (QCP). Dumitrescu et al.
 obtained solutions at inverse temperatures up to $\beta = 65$, limited by the fermionic sign problem 
 of their interaction expansion solver. The hybridization expansion carried out by `CTSEG` is 
 sign-problem-free for the $t-J-U$ model, allowing us to reach $\beta = 300$ and to obtain a 
-more accurate localization of the QCP at doping $p \approx 0.16$ (Fig. 2).
+more accurate localization of the QCP at doping $p \approx 0.16$ (Fig. 1b).
 
 # Acknowledgements