diff --git a/src/GW/GW_plot_self_energy.f90 b/src/GW/GW_plot_self_energy.f90
index 85b36da..f8df205 100644
--- a/src/GW/GW_plot_self_energy.f90
+++ b/src/GW/GW_plot_self_energy.f90
@@ -43,8 +43,8 @@ subroutine GW_plot_self_energy(nBas,eta,nC,nO,nV,nR,nS,eHF,eGW,Om,rho)
 
 ! Minimum and maximum frequency values
 
-  wmin = -5d0
-  wmax = +5d0
+  wmin = -10d0
+  wmax = +0d0
   dw = (wmax - wmin)/dble(ngrid)
 
   do g=1,nGrid
diff --git a/src/GW/RGWC.f90 b/src/GW/RGWC.f90
index 18d9154..167d6d4 100644
--- a/src/GW/RGWC.f90
+++ b/src/GW/RGWC.f90
@@ -29,7 +29,7 @@ subroutine RGWC(dotest,eta,nBas,nC,nO,nV,nR,nS,Om,rho,eHF,e,eGW,Z)
   integer                       :: p,q,i,a,m
   integer                       :: iSat
   double precision              :: num,eps
-  double precision,parameter    :: cutoff = 1d-2
+  double precision,parameter    :: cutoff = 0d-3
 
   logical,parameter             :: do_hole_branch = .true.
   logical,parameter             :: do_electron_branch = .false.
@@ -190,7 +190,7 @@ subroutine RGWC(dotest,eta,nBas,nC,nO,nV,nR,nS,Om,rho,eHF,e,eGW,Z)
 
 ! Construct grid
 
-  nGrid = 1000
+  nGrid = 5000
   allocate(w(nGrid),AGWC(nBas,nGrid))
 
 ! Minimum and maximum frequency values