diff --git a/src/GW/GW_ppBSE_dynamic_kernel_D.f90 b/src/GW/GW_ppBSE_dynamic_kernel_D.f90
index 1497061..0070269 100644
--- a/src/GW/GW_ppBSE_dynamic_kernel_D.f90
+++ b/src/GW/GW_ppBSE_dynamic_kernel_D.f90
@@ -106,7 +106,6 @@ subroutine GW_ppBSE_dynamic_kernel_D(ispin,eta,nBas,nC,nO,nV,nR,nS,nOO,lambda,eG
             kl = kl + 1
   
             do m=1,nS
-
               dem = - OmBSE + eGW(k) - Om(m) + eGW(j)
               num = rho(i,k,m)*rho(j,l,m)
 
diff --git a/src/GW/GW_ppBSE_dynamic_perturbation.f90 b/src/GW/GW_ppBSE_dynamic_perturbation.f90
index 011b91f..4f5b96e 100644
--- a/src/GW/GW_ppBSE_dynamic_perturbation.f90
+++ b/src/GW/GW_ppBSE_dynamic_perturbation.f90
@@ -121,7 +121,7 @@ subroutine GW_ppBSE_dynamic_perturbation(ispin,dTDA,eta,nBas,nC,nO,nV,nR,nS,nOO,
 
       call GW_ppBSE_dynamic_kernel_D(ispin,eta,nBas,nC,nO,nV,nR,nS,nOO,1d0,eGW,OmRPA,rho_RPA,Om2(ij),KD_dyn,ZD_dyn)
  
-      Z2_dyn(kl)  = + dot_product(Y2(:,ij),matmul(ZD_dyn,Y2(:,ij)))
+      Z2_dyn(kl)  = dot_product(Y2(:,ij),matmul(ZD_dyn,Y2(:,ij)))
       Om2_dyn(kl) = - dot_product(Y2(:,ij),matmul(KD_dyn - KD_sta,Y2(:,ij)))
 
     else