diff --git a/plot.py b/plot.py
index 49df0cd72f9dc6c722b8f9204cdcdc5dd45b923f..17deef0492efd785d13a35674843488a6f325462 100644
--- a/plot.py
+++ b/plot.py
@@ -17,6 +17,7 @@ from numbers import Number
 from final_plots import filter_grid_data, photon_assisted_tunneling, gen_photon_assisted_tunneling, interp_vac0
 from scipy.interpolate import bisplrep, bisplev, splrep, BSpline, griddata
 from scipy.optimize import curve_fit, leastsq
+from PIL import Image
 import settings
 from logging import _levelToName
 from data_management import DataManager, KondoImport
@@ -1725,9 +1726,12 @@ def RGeq_comparison_contours(
 
 
 def contours_to_coordinates(contour):
-    for level, collection in zip(contour.levels, contour.collections):
-        for path in collection.get_paths():
-            yield level, path.vertices
+    for level, segments in zip(contour.levels, contour.allsegs):
+        for seg in segments:
+            yield level, seg
+    #for level, collection in zip(contour.levels, contour.collections):
+    #    for path in collection.get_paths():
+    #        yield level, path.vertices
 
 def contours_to_coordinates_unified(contour):
     for level, vertices in contours_to_coordinates(contour):
@@ -2175,7 +2179,7 @@ def compare_limits(
     plt.show()
 
 
-def overview_3d_contours(dm, show_contours=False, show_surfaces=True, offset=0, **trash):
+def overview_3d_contours(dm, show_contours=False, show_surfaces=True, save_images=True, offset=0, **trash):
     import mpl_toolkits.mplot3d
     data_vdc0 = np.load("figdata/vdc0_interp.npz")
     data_omega5 = np.load("figdata/omega5_interp.npz")
@@ -2184,6 +2188,7 @@ def overview_3d_contours(dm, show_contours=False, show_surfaces=True, offset=0,
     fill_levels = 1.1/np.linspace(1, 6.12, 256)[::-1] - 0.1
     colors = plt.cm.viridis(np.linspace(0, 1, 256))
     levels = 1/np.arange(1, 50, 0.5)[::-1]
+    real_cmap = lambda x: plt.cm.viridis(1 - (1.1/(x+0.1)-1)/5.12)[...,:3]
     vdc_resolution = 400
     omega_resolution = 400
     omega_max = 16.5372
@@ -2200,7 +2205,7 @@ def overview_3d_contours(dm, show_contours=False, show_surfaces=True, offset=0,
     cmap_ax.set_xlim(0, 1)
     cmap_ax.set_ylim(min_m, 1)
 
-    vdc0_selection = data_vdc0["omega"][0] <= omega_max
+    vdc0_selection = data_vdc0["omega"][0] <= omega_max + 2e-3
     vdc0_y = data_vdc0["vac_omega"][:,vdc0_selection]
     vdc0_x = np.zeros_like(vdc0_y)
     vdc0_z = data_vdc0["omega"][:,vdc0_selection]
@@ -2209,7 +2214,12 @@ def overview_3d_contours(dm, show_contours=False, show_surfaces=True, offset=0,
     if show_surfaces:
         ax.contourf(vdc0_m, vdc0_y, vdc0_z/TK_VOLTAGE, zdir="x", levels=fill_levels, offset=0, zorder=-10, colors=colors)
     if show_contours:
-        ax.contour(vdc0_m, vdc0_y, vdc0_z/TK_VOLTAGE, zdir="x", levels=levels, offset=-offset, zorder=10, colors="black")
+        contour = ax.contour(vdc0_m, vdc0_y, vdc0_z/TK_VOLTAGE, zdir="x", levels=levels, offset=-offset, zorder=10, colors="black")
+        save_contour_coordinates(contour, "figdata/overview_contour_x.dat", "vdc", "omega", "gdc")
+    if save_images:
+        array = np.array(256*real_cmap(vdc0_m.T[::-1,::-1]), dtype=np.ubyte)
+        img = Image.fromarray(array)
+        img.save("figdata/overview_x.png")
 
     omega5 = data_omega5["omega"]
     omega5_x = data_omega5["vdc"]/omega5
@@ -2218,7 +2228,12 @@ def overview_3d_contours(dm, show_contours=False, show_surfaces=True, offset=0,
     if show_surfaces:
         ax.contourf(omega5_x, omega5_y, omega5_m, zdir="z", levels=fill_levels, offset=omega5/TK_VOLTAGE, colors=colors)
     if show_contours:
-        ax.contour(omega5_x, omega5_y, omega5_m, zdir="z", levels=levels, offset=omega5/TK_VOLTAGE+offset, colors="black")
+        contour = ax.contour(omega5_x, omega5_y, omega5_m, zdir="z", levels=levels, offset=omega5/TK_VOLTAGE+offset, colors="black")
+        save_contour_coordinates(contour, "figdata/overview_contour_z.dat", "vac", "vdc", "gdc")
+    if save_images:
+        array = np.array(256*real_cmap(omega5_m[::-1]), dtype=np.ubyte)
+        img = Image.fromarray(array)
+        img.save("figdata/overview_z.png")
 
     vac0_x = np.linspace(0, 10, vdc_resolution)
     vac0_z = np.linspace(1e-3, omega_max, omega_resolution)
@@ -2227,9 +2242,18 @@ def overview_3d_contours(dm, show_contours=False, show_surfaces=True, offset=0,
     if show_surfaces:
         ax.contourf(vac0_x, vac0_m, vac0_z/TK_VOLTAGE, zdir="y", levels=fill_levels, offset=0, zorder=-10, colors=colors)
     if show_contours:
-        ax.contour(vac0_x, vac0_m, vac0_z/TK_VOLTAGE, zdir="y", levels=levels, offset=-offset, zorder=10, colors="black")
+        contour = ax.contour(vac0_x, vac0_m, vac0_z/TK_VOLTAGE, zdir="y", levels=levels, offset=-offset, zorder=10, colors="black")
+        save_contour_coordinates(contour, "figdata/overview_contour_y.dat", "omega", "vac", "gdc")
+    if save_images:
+        array = np.array(256*real_cmap(vac0_m.T[::-1]), dtype=np.ubyte)
+        img = Image.fromarray(array)
+        img.save("figdata/overview_y.png")
 
     #ax.plot_surface(vdc0_x, vdc0_y, vdc0_z, rstride=1, cstride=1, facecolors=plt.cm.viridis(vdc0_m), shade=False, zorder=-10)
+    #ax.set_axis_off()
+    #fig.set_size_inches(14/2.56, 12/2.56)
+    #fig.savefig("/tmp/figure.png", dpi=512, transparent=True)
+    #ax.set_axis_on()
     plt.show()
 
 
diff --git a/tikz/contour.tex b/tikz/contour.tex
new file mode 100644
index 0000000000000000000000000000000000000000..0a6ca5e4449bb1820e79499607b12aab0fb17638
--- /dev/null
+++ b/tikz/contour.tex
@@ -0,0 +1,39 @@
+\begin{tikzpicture}
+  \begin{axis}[
+      width = 12cm,
+      height = 10cm,
+      xmin = 0, xmax = 10,
+      ymin = 0, ymax = 10,
+      xlabel = {$\vac~(\Omega)$},
+      ylabel = {$\vdc~(\Omega)$},
+      every axis plot post/.append style = {
+        line width=0.2pt,
+        line cap=round,
+        line join=round
+      },
+      legend style = {
+        anchor = south east,
+        at = {(0.99,0.01)}
+      }
+    ]
+    \addplot[green] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_mu_o2.dat};
+    \addplot[green, forget plot] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_J_o2.dat};
+    \addplot[green, forget plot] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_J_o2_p.dat};
+    \addplot[blue] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_mu_o3.dat};
+    \addplot[blue, forget plot] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_J_o3.dat};
+    \addplot[blue, forget plot] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_J_o3_p.dat};
+    \addplot[red] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_mu_o3a.dat};
+    \addplot[red, forget plot] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_J_o3a.dat};
+    \addplot[red, forget plot] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_J_o3a_p.dat};
+    \addplot[black] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_mu_o3p.dat};
+    \addplot[black, forget plot] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_J_o3p.dat};
+    \addplot[black, forget plot] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_J_o3p_p.dat};
+
+    \legend{
+      $O(J^2)$,
+      {$O(J^3)$, $G^a=0$},
+      {$O(J^3)$, approx.~$I[X]$},
+      $O(J^3)$,
+    }
+  \end{axis}
+\end{tikzpicture}
diff --git a/tikz/contour_wrapper.tex b/tikz/contour_wrapper.tex
new file mode 100644
index 0000000000000000000000000000000000000000..fd51447292cb401e7b74d4b6c8425c94a76b268f
--- /dev/null
+++ b/tikz/contour_wrapper.tex
@@ -0,0 +1,59 @@
+\documentclass{standalone}
+\newcommand\tikzsetnextfilename[1]{}
+\newcommand\tkv{T_\mathrm{K}}
+\newcommand\vdc{V_\mathrm{avg}}
+\newcommand\vac{V_\mathrm{osc}}
+\newcommand\gdc{G_\mathrm{avg}}
+\newcommand\gac{G_\mathrm{osc}}
+\usepackage{pgfplots}
+\pgfplotsset{compat=1.18}
+\pgfplotsset{surf shading/precision=pdf}
+\usepackage{fontspec}
+\defaultfontfeatures{Ligatures=TeX}
+\setsansfont[
+  UprightFont = *-Regular,
+  BoldFont = *-Bold,
+  ItalicFont = *-Italic,
+  %BoldItalicFont = LinBiolinum_K, % TODO: Check this font, use a better one
+  SmallCapsFont = *-Regular,
+  SmallCapsFeatures = {Letters=SmallCaps},
+  SlantedFont = *-Regular,
+  SlantedFeatures = {FakeSlant=0.15},
+]{LibertinusSans}
+\setmainfont[
+  %UprightFont = *Display-Regular,
+  UprightFont = *-Regular,
+  BoldFont = *-Semibold,
+  ItalicFont = *-Italic,
+  BoldItalicFont = *-SemiboldItalic,
+  SmallCapsFont = *-Regular,
+  SmallCapsFeatures = {Letters=SmallCaps},
+  %SlantedFont = *-Italic,
+  SlantedFont = *-Regular,
+  SlantedFeatures = {FakeSlant=0.15},
+]{LibertinusSerif}
+%\setmonofont[
+%  UprightFont = *-Regular
+%]{LibertinusMono}
+\setmonofont[
+  UprightFont = *-Regular,
+  BoldFont = *-Medium,
+  ItalicFont = *-It,
+  BoldItalicFont = *-MediumIt,
+]{SourceCodePro}
+
+\usepackage{mathtools}
+%\usepackage{amsmath}
+\usepackage{amssymb}
+%\usepackage{amsfonts}
+\usepackage[math-style=ISO, bold-style=TeX, partial=upright, nabla=upright]{unicode-math}
+\setmathfont{Libertinus Math}
+\setmathfont{STIXTwoMath-Regular.otf}[range=cal]
+\setmathfont{STIXTwoMath-Regular.otf}[range={\mathscr},StylisticSet=1]
+\setmathfont{STIXTwoMath-Regular.otf}[range={"02022, "02218, "02200, "0223C, "02248}] % Slightly larger \bullet and \circ; larger \forall; different \sim, \approx
+%\setmathfont{Latin Modern Math}[range={"027F8-"027FA}] % \iff, \implies, \impliedby
+%\setmathfont{STIX Two Math}[range={"027E8, "027E9}] % \langle, \rangle
+\setmathfont{STIXTwoMath-Regular.otf}[range={"0221A-"0221C, "023B7}] % root (sqrt, cbrt, \dots) symbols
+\begin{document}
+\input{contour}
+\end{document}
diff --git a/tikz/g_overview_vdc_vac.tex b/tikz/g_overview_vdc_vac.tex
new file mode 100644
index 0000000000000000000000000000000000000000..f419e8065bf63e999790bc49d10b3e09904227c7
--- /dev/null
+++ b/tikz/g_overview_vdc_vac.tex
@@ -0,0 +1,62 @@
+\tikzsetnextfilename{gacdc_vdc_vac}%
+\begin{tikzpicture}
+  \begin{axis}[
+      name = G3d,
+      view={110}{22},
+      width = 0.84\linewidth,
+      height = 0.72\linewidth,
+      zmin = 0, zmax = 1,
+      ymin = 0, ymax = 10,
+      xmin = 0, xmax = 10,
+      xlabel = {$\vac~(\Omega)$},
+      ylabel = {$\vdc~(\Omega)$},
+      zlabel = {$\gdc~(2e^2/h)$},
+      x label style = {sloped},
+      y label style = {sloped},
+      label style = {font=\small},
+      every tick label/.append style = {font=\small},
+      colormap/viridis,
+      clip = false,
+      scale only axis,
+    ]
+    \addplot3[
+      surf,
+      shader = faceted interp,
+      point meta min = 0,
+      point meta max = 1,
+      mesh/rows = 101,
+      point meta = \thisrow{gdc},
+      line width = 0.1pt,
+    ] table [x=vac, y=vdc, z=gdc]{../figdata/vdc_vac_omega16.5372.dat};
+    \addplot3[color=black         ] table [x expr=10, y=vdc, z=gdc]{../figdata/vdc_vac0_omega16.5372.dat};
+    \addplot3[color=black!40!red  ] table [x expr=10, y=vdc, z=gdc]{../figdata/vdc_vac1_omega16.5372.dat};
+    \addplot3[color=red           ] table [x expr=10, y=vdc, z=gdc]{../figdata/vdc_vac2_omega16.5372.dat};
+    \addplot3[color=red!50!blue   ] table [x expr=10, y=vdc, z=gdc]{../figdata/vdc_vac3_omega16.5372.dat};
+    \addplot3[color=blue          ] table [x expr=10, y=vdc, z=gdc]{../figdata/vdc_vac4_omega16.5372.dat};
+    \addplot3[color=blue!50!green ] table [x expr=10, y=vdc, z=gdc]{../figdata/vdc_vac5_omega16.5372.dat};
+    \addplot3[color=green!80!black] table [x expr=10, y=vdc, z=gdc]{../figdata/vdc_vac6_omega16.5372.dat};
+    \addplot3[color=green!40!red  ] table [x expr=10, y=vdc, z=gdc]{../figdata/vdc_vac7_omega16.5372.dat};
+    \addplot3[color=green!40!black] table [x expr=10, y=vdc, z=gdc]{../figdata/vdc_vac8_omega16.5372.dat};
+    \addplot3[color=blue!50!black ] table [x expr=10, y=vdc, z=gdc]{../figdata/vdc_vac9_omega16.5372.dat};
+    %
+    \addplot3[color=black         ] table [x=vac, y expr=10, z=gdc]{../figdata/vac_vdc0_omega16.5372.dat};
+    \addplot3[color=black!40!red  ] table [x=vac, y expr=10, z=gdc]{../figdata/vac_vdc1_omega16.5372.dat};
+    \addplot3[color=red           ] table [x=vac, y expr=10, z=gdc]{../figdata/vac_vdc2_omega16.5372.dat};
+    \addplot3[color=red!50!blue   ] table [x=vac, y expr=10, z=gdc]{../figdata/vac_vdc3_omega16.5372.dat};
+    \addplot3[color=blue          ] table [x=vac, y expr=10, z=gdc]{../figdata/vac_vdc4_omega16.5372.dat};
+    \addplot3[color=blue!50!green ] table [x=vac, y expr=10, z=gdc]{../figdata/vac_vdc5_omega16.5372.dat};
+    \addplot3[color=green!80!black] table [x=vac, y expr=10, z=gdc]{../figdata/vac_vdc6_omega16.5372.dat};
+    \addplot3[color=green!40!red  ] table [x=vac, y expr=10, z=gdc]{../figdata/vac_vdc7_omega16.5372.dat};
+    \addplot3[color=green!40!black] table [x=vac, y expr=10, z=gdc]{../figdata/vac_vdc8_omega16.5372.dat};
+    \addplot3[color=blue!50!black ] table [x=vac, y expr=10, z=gdc]{../figdata/vac_vdc9_omega16.5372.dat};
+    %\addplot3[
+    %  surf,
+    %  shader = faceted interp,
+    %  point meta min = 0,
+    %  point meta max = 1,
+    %  point meta = \thisrow{gdc},
+    %  line width = 0.1pt,
+    %] table [x=vac, y=vdc, z=gdc]{../figdata/vdc_vac_omega16.5372_zoom.dat};
+    \node[black] at (axis description cs:0.75,0.82) {$\Omega=5\tkv$};
+  \end{axis}
+\end{tikzpicture}
diff --git a/tikz/g_overview_vdc_vac_old.tex b/tikz/g_overview_vdc_vac_old.tex
new file mode 100644
index 0000000000000000000000000000000000000000..ad72a2c930c20000a176558a00d4dd79f5214af7
--- /dev/null
+++ b/tikz/g_overview_vdc_vac_old.tex
@@ -0,0 +1,60 @@
+\tikzsetnextfilename{gacdc_vdc_vac}%
+\begin{tikzpicture}
+  \begin{axis}[
+      name = G3d,
+      view={110}{22},
+      width = 0.84\linewidth,
+      height = 0.72\linewidth,
+      zmin = 0.075, zmax = 80,
+      ymin = 0, ymax = 10,
+      xmin = 0, xmax = 10,
+      zmode = log,
+      xlabel = {$\vac~(\Omega)$},
+      ylabel = {$\vdc~(\Omega)$},
+      zlabel = {$\quad\gac~(2e^2/h)\qquad\quad\gdc~(2e^2/h)$},
+      x label style = {sloped},
+      y label style = {sloped},
+      label style = {font=\small},
+      every tick label/.append style = {font=\small},
+      label shift = {-.8ex},
+      zlabel shift = {-1.5ex},
+      ztick = {0.1, 1, 8, 80},
+      zticklabels = {$0.1$, $1$, $0.1$, $1$},
+      minor ztick = {0.08,0.09,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9, 5.6,6.4,7.2,16,32,40,48,56,64,72},
+      %zlabel style = {at = (ticklabel cs:0.65)},
+      xtick = {0,4,8},
+      xticklabels = {$\hspace*{-.5em}0$, $4$, $8$},
+      ytick = {0,2,4,6,8},
+      yticklabels = {$\hspace*{.5em}0$, $2$, $4$, $6$, $8\hspace*{.5em}$},
+      colormap/viridis,
+      clip = false,
+      scale only axis,
+    ]
+    \addplot3[gray, no marks, line width=0.2pt] coordinates {(10,0,0.3) (0,0,0.3) (0,10,0.3)};
+    \addplot3[
+      line width = 0.18pt,
+      surf,
+      point meta min = -1.15,
+      point meta max = -0.4,
+      point meta = {log10(pi*\thisrow{gac})},
+    ] table [x=vac, y=vdc, z expr=gac]{../figdata/vdc_vac_omega10.dat};
+    \addplot3[gray, no marks, line width=0.2pt] coordinates {(10,0,10) (0,0,10) (0,10,10)};
+    \addplot3[
+      line width = 0.18pt,
+      surf,
+      point meta min = -1.15,
+      point meta max = -0.4,
+      point meta = {log10(pi*\thisrow{gdc_g})},
+    ] table [x=vac, y=vdc, z expr={80*gdc_g}]{../figdata/vdc_vac_omega10.dat};
+    %\addplot3[
+    %  red,
+    %  line width = 0.5pt,
+    %] table [x expr=4.8, y expr={\thisrow{vdc}/10}, z expr={80*\thisrow{gdc}}]{../figdata/vdc_vac120_omega10.dat};
+    %\addplot3[
+    %  red,
+    %  line width = 0.5pt,
+    %] table [y expr=4, x expr={\thisrow{vac}/10}, z expr={80*\thisrow{gdc}}]{../figdata/vac_vdc100_omega10_cutoff.dat};
+    \addplot3[gray, no marks, line width=0.2pt] coordinates {(8,0,8) (8,8,8) (0,8,8)};
+    \node[black] at (axis description cs:0.75,0.82) {$\Omega=7.55\tkv$};
+  \end{axis}
+\end{tikzpicture}
diff --git a/tikz/g_overview_vdc_vac_wrapper.tex b/tikz/g_overview_vdc_vac_wrapper.tex
new file mode 100644
index 0000000000000000000000000000000000000000..c8580d48db6f2ced7f4f5300955956d5e7ec7f7a
--- /dev/null
+++ b/tikz/g_overview_vdc_vac_wrapper.tex
@@ -0,0 +1,13 @@
+\documentclass{standalone}
+\newcommand\tikzsetnextfilename[1]{}
+\newcommand\tkv{T_\mathrm{K}}
+\newcommand\vdc{V_\mathrm{dc}}
+\newcommand\vac{V_\mathrm{ac}}
+\newcommand\gdc{G_\mathrm{dc}}
+\newcommand\gac{G_\mathrm{ac}}
+\usepackage{pgfplots}
+\pgfplotsset{compat=1.18}
+\pgfplotsset{surf shading/precision=pdf}
+\begin{document}
+\input{g_overview_vdc_vac}
+\end{document}
diff --git a/tikz/omega5-3d.tex b/tikz/omega5-3d.tex
new file mode 100644
index 0000000000000000000000000000000000000000..7492d2ce1e578dabbd020ac81b9de1a194fb5a8e
--- /dev/null
+++ b/tikz/omega5-3d.tex
@@ -0,0 +1,52 @@
+%\tikzsetnextfilename{gacdc_vdc_vac}%
+\begin{tikzpicture}
+  \begin{axis}[
+      name = G3d,
+      view={110}{22},
+      width = 7cm,
+      height = 7cm,
+      zmin = 0.075, zmax = 80,
+      ymin = 0, ymax = 10,
+      xmin = 0, xmax = 10,
+      zmode = log,
+      xlabel = {$\vac~(\Omega)$},
+      ylabel = {$\vdc~(\Omega)$},
+      zlabel = {$\quad\gac~(2e^2/h)\qquad\quad\gdc~(2e^2/h)$},
+      x label style = {sloped},
+      y label style = {sloped},
+      label style = {font=\small},
+      every tick label/.append style = {font=\small},
+      label shift = {-.8ex},
+      zlabel shift = {-1.5ex},
+      ztick = {0.1, 1, 8, 80},
+      zticklabels = {$0.1$, $1$, $0.1$, $1$},
+      minor ztick = {0.08,0.09,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9, 5.6,6.4,7.2,16,32,40,48,56,64,72},
+      %zlabel style = {at = (ticklabel cs:0.65)},
+      xtick = {0,5,10},
+      xticklabels = {$\hspace*{-.5em}0$, $5$, $10$},
+      ytick = {0,2,4,6,8,10},
+      yticklabels = {$\hspace*{.5em}0$, $2$, $4$, $6$, $8$, $10\hspace*{.5em}$},
+      colormap/viridis,
+      clip = false,
+      scale only axis,
+    ]
+    \addplot3[gray, no marks, line width=0.2pt] coordinates {(10,0,0.3) (0,0,0.3) (0,10,0.3)};
+    \addplot3[
+      line width = 0.18pt,
+      surf,
+      point meta min = -1.15,
+      point meta max = -0.4,
+      point meta = {log10(\thisrow{gac})},
+    ] table [x expr={\thisrow{vac}/5.000007155}, y expr={\thisrow{vdc}/5.000007155}, z expr={\thisrow{gac}}]{../figdata/vdc_vac_omega16.5372_interp.dat};
+    \addplot3[gray, no marks, line width=0.2pt] coordinates {(10,0,8) (0,0,8) (0,10,8)};
+    \addplot3[
+      line width = 0.18pt,
+      surf,
+      point meta min = -1.15,
+      point meta max = -0.4,
+      point meta = {log10(\thisrow{gdc_g})},
+    ] table [x expr={\thisrow{vac}/5.000007155}, y expr={\thisrow{vdc}/5.000007155}, z expr={80*\thisrow{gdc_g}}]{../figdata/vdc_vac_omega16.5372_interp.dat};
+    \addplot3[gray, no marks, line width=0.2pt] coordinates {(10,0,8) (10,10,8) (0,10,8)};
+    \node[black] at (axis description cs:0.75,0.82) {$\Omega=4.8\tkv$};
+  \end{axis}
+\end{tikzpicture}
diff --git a/tikz/omega5-3d_wrapper.tex b/tikz/omega5-3d_wrapper.tex
new file mode 100644
index 0000000000000000000000000000000000000000..6eeb7c4a107c93d22a63e19c32a1144be7076f34
--- /dev/null
+++ b/tikz/omega5-3d_wrapper.tex
@@ -0,0 +1,59 @@
+\documentclass{standalone}
+\newcommand\tikzsetnextfilename[1]{}
+\newcommand\tkv{T_\mathrm{K}}
+\newcommand\vdc{V_\mathrm{avg}}
+\newcommand\vac{V_\mathrm{osc}}
+\newcommand\gdc{G_\mathrm{avg}}
+\newcommand\gac{G_\mathrm{osc}}
+\usepackage{pgfplots}
+\pgfplotsset{compat=1.18}
+\pgfplotsset{surf shading/precision=pdf}
+\usepackage{fontspec}
+\defaultfontfeatures{Ligatures=TeX}
+\setsansfont[
+  UprightFont = *-Regular,
+  BoldFont = *-Bold,
+  ItalicFont = *-Italic,
+  %BoldItalicFont = LinBiolinum_K, % TODO: Check this font, use a better one
+  SmallCapsFont = *-Regular,
+  SmallCapsFeatures = {Letters=SmallCaps},
+  SlantedFont = *-Regular,
+  SlantedFeatures = {FakeSlant=0.15},
+]{LibertinusSans}
+\setmainfont[
+  %UprightFont = *Display-Regular,
+  UprightFont = *-Regular,
+  BoldFont = *-Semibold,
+  ItalicFont = *-Italic,
+  BoldItalicFont = *-SemiboldItalic,
+  SmallCapsFont = *-Regular,
+  SmallCapsFeatures = {Letters=SmallCaps},
+  %SlantedFont = *-Italic,
+  SlantedFont = *-Regular,
+  SlantedFeatures = {FakeSlant=0.15},
+]{LibertinusSerif}
+%\setmonofont[
+%  UprightFont = *-Regular
+%]{LibertinusMono}
+\setmonofont[
+  UprightFont = *-Regular,
+  BoldFont = *-Medium,
+  ItalicFont = *-It,
+  BoldItalicFont = *-MediumIt,
+]{SourceCodePro}
+
+\usepackage{mathtools}
+%\usepackage{amsmath}
+\usepackage{amssymb}
+%\usepackage{amsfonts}
+\usepackage[math-style=ISO, bold-style=TeX, partial=upright, nabla=upright]{unicode-math}
+\setmathfont{Libertinus Math}
+\setmathfont{STIXTwoMath-Regular.otf}[range=cal]
+\setmathfont{STIXTwoMath-Regular.otf}[range={\mathscr},StylisticSet=1]
+\setmathfont{STIXTwoMath-Regular.otf}[range={"02022, "02218, "02200, "0223C, "02248}] % Slightly larger \bullet and \circ; larger \forall; different \sim, \approx
+%\setmathfont{Latin Modern Math}[range={"027F8-"027FA}] % \iff, \implies, \impliedby
+%\setmathfont{STIX Two Math}[range={"027E8, "027E9}] % \langle, \rangle
+\setmathfont{STIXTwoMath-Regular.otf}[range={"0221A-"0221C, "023B7}] % root (sqrt, cbrt, \dots) symbols
+\begin{document}
+\input{omega5-3d}
+\end{document}
diff --git a/tikz/overview.tex b/tikz/overview.tex
new file mode 100644
index 0000000000000000000000000000000000000000..cf5a4481af99e6e22a94b1837db7537e4fb7d67d
--- /dev/null
+++ b/tikz/overview.tex
@@ -0,0 +1,36 @@
+\begin{tikzpicture}
+  \begin{axis}[
+      width=11cm,
+      height=11cm,
+      scale only axis,
+      xmin=0, xmax=10,
+      ymin=0, ymax=10,
+      zmin=0, zmax=4.804,
+      x label style = {sloped},
+      y label style = {sloped},
+      xlabel = {$\vdc/\Omega$},
+      ylabel = {$\vac/\Omega$},
+      zlabel = {$\Omega~(\tkv)$},
+    ]
+    \addplot3[surf] graphics[
+      points={
+        %(10,10,4.804) => (757.5,1644)
+        %%(0,10,4.804) => (0,1151)
+        %(10,0,4.804) => (1515,1151)
+        %(10,0,0) => (1471,549)
+        %(0,0,0) => (757.5,0)
+        (10,10,4.804) => (400,1000)
+        %(10,0,4.804) => (800,700)
+        (10,0,0) => (800,300)
+        (0,0,0) => (400,0)
+        (0,0,4.804) => (400,400)
+      }
+    ] {overview_pixmap.png};
+    \addplot3[black] table[x expr={0}, y=vdc, z=omega]{../figdata/overview_contour_x.dat};
+    \addplot3[black] table[x=vac, y expr={0}, z=omega]{../figdata/overview_contour_y.dat};
+    \addplot3[black] table[x=vac, y=vdc, z expr={4.804}]{../figdata/overview_contour_z.dat};
+    \addplot3[gray,thin] coordinates {(0,0,0) (0,0,4.804)};
+    \addplot3[gray,thin] coordinates {(10,0,4.804) (0,0,4.804) (0,10,4.804)};
+    %\addplot3[surf] graphics[points={(0,0,0)=>(0,0)}] {../figdata/overview_x.png};
+  \end{axis}
+\end{tikzpicture}
diff --git a/tikz/overview_wrapper.tex b/tikz/overview_wrapper.tex
new file mode 100644
index 0000000000000000000000000000000000000000..4efc8fc54b9621be2cd63bf7375fbd85fccb1fc2
--- /dev/null
+++ b/tikz/overview_wrapper.tex
@@ -0,0 +1,59 @@
+\documentclass{standalone}
+\newcommand\tikzsetnextfilename[1]{}
+\newcommand\tkv{T_\mathrm{K}}
+\newcommand\vdc{V_\mathrm{avg}}
+\newcommand\vac{V_\mathrm{osc}}
+\newcommand\gdc{G_\mathrm{avg}}
+\newcommand\gac{G_\mathrm{osc}}
+\usepackage{pgfplots}
+\pgfplotsset{compat=1.18}
+\pgfplotsset{surf shading/precision=pdf}
+\usepackage{fontspec}
+\defaultfontfeatures{Ligatures=TeX}
+\setsansfont[
+  UprightFont = *-Regular,
+  BoldFont = *-Bold,
+  ItalicFont = *-Italic,
+  %BoldItalicFont = LinBiolinum_K, % TODO: Check this font, use a better one
+  SmallCapsFont = *-Regular,
+  SmallCapsFeatures = {Letters=SmallCaps},
+  SlantedFont = *-Regular,
+  SlantedFeatures = {FakeSlant=0.15},
+]{LibertinusSans}
+\setmainfont[
+  %UprightFont = *Display-Regular,
+  UprightFont = *-Regular,
+  BoldFont = *-Semibold,
+  ItalicFont = *-Italic,
+  BoldItalicFont = *-SemiboldItalic,
+  SmallCapsFont = *-Regular,
+  SmallCapsFeatures = {Letters=SmallCaps},
+  %SlantedFont = *-Italic,
+  SlantedFont = *-Regular,
+  SlantedFeatures = {FakeSlant=0.15},
+]{LibertinusSerif}
+%\setmonofont[
+%  UprightFont = *-Regular
+%]{LibertinusMono}
+\setmonofont[
+  UprightFont = *-Regular,
+  BoldFont = *-Medium,
+  ItalicFont = *-It,
+  BoldItalicFont = *-MediumIt,
+]{SourceCodePro}
+
+\usepackage{mathtools}
+%\usepackage{amsmath}
+\usepackage{amssymb}
+%\usepackage{amsfonts}
+\usepackage[math-style=ISO, bold-style=TeX, partial=upright, nabla=upright]{unicode-math}
+\setmathfont{Libertinus Math}
+\setmathfont{STIXTwoMath-Regular.otf}[range=cal]
+\setmathfont{STIXTwoMath-Regular.otf}[range={\mathscr},StylisticSet=1]
+\setmathfont{STIXTwoMath-Regular.otf}[range={"02022, "02218, "02200, "0223C, "02248}] % Slightly larger \bullet and \circ; larger \forall; different \sim, \approx
+%\setmathfont{Latin Modern Math}[range={"027F8-"027FA}] % \iff, \implies, \impliedby
+%\setmathfont{STIX Two Math}[range={"027E8, "027E9}] % \langle, \rangle
+\setmathfont{STIXTwoMath-Regular.otf}[range={"0221A-"0221C, "023B7}] % root (sqrt, cbrt, \dots) symbols
+\begin{document}
+\input{overview}
+\end{document}