diff --git a/data/simulations/lif-mu011-D0010-chi2.npz b/data/simulations/lif-mu011-D0010-chi2.npz
index 67b74f2..e7976c1 100644
Binary files a/data/simulations/lif-mu011-D0010-chi2.npz and b/data/simulations/lif-mu011-D0010-chi2.npz differ
diff --git a/lifsuscept.py b/lifsuscept.py
index 3b94aa8..dbf079f 100644
--- a/lifsuscept.py
+++ b/lifsuscept.py
@@ -84,7 +84,8 @@ def load_lifdata(mu, D, vT=1, vR=0, tau_ref=0):
     file_path = sims_path / f'lif-mu{10*mu:03.0f}-D{10000*D:04.0f}-chi2.npz'
     if not file_path.exists():
         freqs1 = np.linspace(0.01, 1.0, 500)
-        freqs2 = np.linspace(0.01, 1.0, 200)
+        #freqs2 = np.linspace(0.01, 1.0, 200)
+        freqs2 = np.linspace(-1.01, 1.0, 400)
         r0, chi1, chi2 = susceptibilities(freqs1, freqs2, mu, D,
                                           tau_ref, vR, vT)
         np.savez(file_path, mu=mu, D=D, vT=vT, vR=vR,
@@ -120,13 +121,18 @@ def plot_chi2(ax, s, r0, freqs, chi2):
     vmax = 300
     pc = ax.pcolormesh(freqs, freqs, chi2, vmin=0, vmax=vmax,
                        rasterized=True)
+    ax.plot([0, 0, 1], [1, 0, 0], **s.lsDiag)
     ax.set_aspect('equal')
     ax.set_xlabel('$f_1$')
     ax.set_ylabel('$f_2$', labelpad=6)
-    ax.set_xlim(0, 1)
-    ax.set_ylim(0, 1)
-    ax.set_xticks_delta(0.2)
-    ax.set_yticks_delta(0.2)
+    #ax.set_xlim(0, 1)
+    #ax.set_ylim(0, 1)
+    #ax.set_xticks_delta(0.2)
+    #ax.set_yticks_delta(0.2)
+    ax.set_xlim(-1, 1)
+    ax.set_ylim(-1, 1)
+    ax.set_xticks_delta(0.5)
+    ax.set_yticks_delta(0.5)
     cax = ax.inset_axes([1.04, 0, 0.05, 1])
     cax.set_spines_outward('lrbt', 0)
     cb = fig.colorbar(pc, cax=cax)