diff --git a/src/MBPT/dynamic_Tmatrix_A.f90 b/src/MBPT/dynamic_Tmatrix_A.f90
index 469ebe8..2120a3a 100644
--- a/src/MBPT/dynamic_Tmatrix_A.f90
+++ b/src/MBPT/dynamic_Tmatrix_A.f90
@@ -68,12 +68,12 @@ subroutine dynamic_Tmatrix_A(eta,nBas,nC,nO,nV,nR,nS,nOO,nVV,lambda,eGT,Omega1,O
           chi = 0d0
 
           do cd=1,nVV
-            eps = + OmBSE - Omega1(cd) - (eGT(a) - eGT(j))
+            eps = + OmBSE - Omega1(cd) + (eGT(i) + eGT(j))
             chi = chi + rho1(i,j,cd)*rho1(a,b,cd)*eps/(eps**2 + eta**2) 
           end do
 
           do kl=1,nOO
-            eps = + OmBSE - Omega2(kl) - (eGT(b) - eGT(i))
+            eps = + OmBSE + Omega2(kl) - (eGT(a) + eGT(b))
             chi = chi + rho2(i,j,kl)*rho2(a,b,kl)*eps/(eps**2 + eta**2)
           end do
 
diff --git a/src/MBPT/dynamic_Tmatrix_ZA.f90 b/src/MBPT/dynamic_Tmatrix_ZA.f90
index fe58db4..ad7cbee 100644
--- a/src/MBPT/dynamic_Tmatrix_ZA.f90
+++ b/src/MBPT/dynamic_Tmatrix_ZA.f90
@@ -54,12 +54,12 @@ subroutine dynamic_Tmatrix_ZA(eta,nBas,nC,nO,nV,nR,nS,nOO,nVV,lambda,eGT,Omega1,
  
           chi = 0d0
           do cd=1,nVV
-            eps = + OmBSE - Omega1(cd) - (eGT(a) - eGT(j))
+            eps = + OmBSE - Omega1(cd) + (eGT(i) + eGT(j))
             chi = chi + rho1(i,j,cd)*rho1(a,b,cd)*(eps**2 - eta**2)/(eps**2 + eta**2)**2
           end do
 
           do kl=1,nOO
-            eps = + OmBSE - Omega2(kl) - (eGT(b) - eGT(i))
+            eps = + OmBSE + Omega2(kl) - (eGT(a) + eGT(b))
             chi = chi + rho2(i,j,kl)*rho2(a,b,kl)*(eps**2 - eta**2)/(eps**2 + eta**2)**2
           end do