diff --git a/python/sumk_dft.py b/python/sumk_dft.py
index a2cf6f2e..29c44c0f 100644
--- a/python/sumk_dft.py
+++ b/python/sumk_dft.py
@@ -259,8 +259,9 @@ class SumkDFT:
 
         # Are we including Sigma?
         if with_Sigma:
-            if with_dc: sigma_minus_dc = self.add_dc(iw_or_w)
             Sigma_imp = getattr(self,"Sigma_imp_"+iw_or_w)
+            sigma_minus_dc = Sigma_imp.copy()
+            if with_dc: sigma_minus_dc = self.add_dc(iw_or_w)
             if iw_or_w == "iw":
                 beta = Sigma_imp[0].mesh.beta   # override beta if Sigma_iw is present
                 n_iw = len(Sigma_imp[0].mesh)