From 8e67a5a75e4ae916a7649c36144fef7cf685feb7 Mon Sep 17 00:00:00 2001
From: Valentin Bruch <valentin.bruch@rwth-aachen.de>
Date: Sun, 12 Feb 2023 23:06:03 +0100
Subject: [PATCH] improved TikZ plots

---
 final_plots.py                   |  26 +++-
 kondo.py                         |   2 +-
 plot.py                          |  27 ++++
 tikz/asymmetry.tex               |  16 +--
 tikz/bruhat18.tex                | 231 ++++++++++++++-----------------
 tikz/commutator_RG2.tex          |   1 +
 tikz/contour.tex                 |  29 ++--
 tikz/convergence-j.tex           |  11 +-
 tikz/convergence-o2.tex          |  11 +-
 tikz/convergence-o3.tex          |  11 +-
 tikz/convergence-o3a.tex         |  11 +-
 tikz/convergence-o3a_ispline.tex |  11 +-
 tikz/current_convergence.tex     |   1 +
 tikz/current_convergence_bad.tex |   1 +
 tikz/floquet_matrices.tex        |   3 +-
 tikz/g_approximations.tex        |  15 +-
 tikz/ga.tex                      |   5 +-
 tikz/harmonic_current.tex        |   3 +-
 tikz/harmonic_modes.tex          |   3 +-
 tikz/kogan04.tex                 |  98 +++++--------
 tikz/omega5-3d.tex               |   2 +-
 tikz/overview.tex                |   1 +
 tikz/preamble.tex                |  90 ++++++++++++
 tikz/pulse_charge.tex            |  15 +-
 tikz/pulse_current.tex           |  15 +-
 25 files changed, 353 insertions(+), 286 deletions(-)
 create mode 100644 tikz/preamble.tex

diff --git a/final_plots.py b/final_plots.py
index 171e5dc..bfbf17a 100644
--- a/final_plots.py
+++ b/final_plots.py
@@ -886,35 +886,49 @@ def prepare_RGconvergence():
     data_vb7 = np.load("figdata/omega5_interp_vb7.npz")
     data_rel = np.load("figdata/omega5_interp_deviation.npz")
 
+    file = open("tikz/colorscales.tex", "w")
+
     # O2
     omega = 16.5372
     omega_o2 = 48.695054
     vdc = data["vdc"][0] / omega
     vac = data["vac"][:,0] / omega
 
+    print(f"extent: xmin={1.5*vdc[0]-0.5*vdc[1]}, xmax={1.5*vdc[-1]-0.5*vdc[-2]}, ymin={1.5*vac[0]-0.5*vac[1]}, ymax={1.5*vac[-1]-0.5*vac[-2]}")
+    print(f"""\\newcommand\\extentxmin{{{1.5*vdc[0]-0.5*vdc[1]}}}
+\\newcommand\\extentxmax{{{1.5*vdc[-1]-0.5*vdc[-2]}}}
+\\newcommand\\extentymin{{{1.5*vac[0]-0.5*vac[1]}}}
+\\newcommand\\extentymax{{{1.5*vac[-1]-0.5*vac[-2]}}}""", file=file)
+
     for method, u, ref in (("o2","mu","o3p"), ("o3","mu","o3p"), ("o3a","mu","o3p"), ("o3a_ispline","mu","o3a"), ("o3a","J","o3a"), ("o3a_p","J","o3a")):
         gdiff = data[f"gdc_{u}_{method}"] - data[f"gdc_mu_{ref}"]
         if u != "mu":
             method += "_" + u
         normalize = np.nanmax(np.abs(gdiff))
+        tex_name = method.replace('_', 'x').replace('o2','ol').replace('o3','ox').replace('vb7','vb')
         export_img(f"convergence_{method}_diff", gdiff/normalize)
+        print(f"\\newcommand\\scale{tex_name}diff{{{np.pi*normalize}}}% convergence_{method}_diff.png", file=file)
         print(f"Color scale for {method} difference: {np.pi*normalize}")
         grel = gdiff/data[f"gdc_mu_{ref}"]
         normalize = np.nanmax(np.abs(grel))
         export_img(f"convergence_{method}_relative", grel/normalize)
         print(f"Color scale for {method} relative: {normalize}")
+        print(f"\\newcommand\\scale{tex_name}rel{{{normalize}}}% convergence_{method}_relative.png", file=file)
 
     for method, u, sign in (("o3p","mu",-1), ("o3a_vb7","mu",-1), ("o3a","J",1), ("o3a_p","J",1)):
         gdiff = sign*data_rel[f"gdc_{u}_{method}"]
         if u != "mu":
             method += "_" + u
         normalize = np.nanmax(np.abs(gdiff))
+        tex_name = method.replace('_', 'x').replace('o2','ol').replace('o3','ox').replace('vb7','vb')
         export_img(f"convergence_{method}_diff_fromdiff", gdiff/normalize)
         print(f"Color scale for {method} difference (interpolate in the end): {np.pi*normalize}")
+        print(f"\\newcommand\\scale{tex_name}diffX{{{np.pi*normalize}}}% convergence_{method}_diff_fromdiff.png", file=file)
         grel = gdiff/data[f"gdc_mu_o3a"]
         normalize = np.nanmax(np.abs(grel))
         export_img(f"convergence_{method}_relative_fromdiff", grel/normalize)
         print(f"Color scale for {method} relative (interpolate in the end): {normalize}")
+        print(f"\\newcommand\\scale{tex_name}relX{{{normalize}}}% convergence_{method}_relative_fromdiff.png", file=file)
 
     export_img(f"convergence_o3a_J_relative_fromdiff_pnorm", data_rel[f"gdc_J_o3a"]/normalize)
 
@@ -925,23 +939,27 @@ def prepare_RGconvergence():
     g1img = np.array(0xffff*real_cmap(g1rel[::-1]/normalize), dtype=np.uint16)
     g2img = np.array(0xffff*real_cmap(g2rel[::-1]/normalize), dtype=np.uint16)
     g2img[-16:,:8] = 0x8000
-    imwrite(f"figdata/convergence_o3a_J_relative_fromdiff_maskednorm.png", g1img, format="PNG-FI")
-    imwrite(f"figdata/convergence_o3a_J_p_relative_fromdiff_masked.png", g2img, format="PNG-FI")
+    imwrite("figdata/convergence_o3a_J_relative_fromdiff_maskednorm.png", g1img, format="PNG-FI")
+    imwrite("figdata/convergence_o3a_J_p_relative_fromdiff_masked.png", g2img, format="PNG-FI")
     print(f"Color scale for o3a J relative (interpolate in the end, masked): {normalize}")
+    print(f"\\newcommand\\scaleJMaskedRelX{{{normalize}}}% convergence_o3a_J_relative_fromdiff_maskednorm.png, convergence_o3a_J_p_relative_fromdiff_masked.png", file=file)
 
     grel = -data_rel["gdc_mu_o3a_vb7"]/data["gdc_mu_o3a"]
     normalize = max(np.nanmax(np.abs(grel[16:])), np.nanmax(np.abs(grel[:16,8:])))
-    export_img(f"convergence_o3a_vb7_relative_fromdiff_masked", grel/normalize)
+    export_img("convergence_o3a_vb7_relative_fromdiff_masked", grel/normalize)
     print(f"Color scale for vb7 masked relative (interpolate in the end): {normalize}")
+    print(f"\\newcommand\\scalevbMaskedRelX{{{normalize}}}% convergence_o3a_vb7_relative_fromdiff_masked.png", file=file)
 
     gdiff = data["gdc_mu_o3a"] - data_vb7["gdc_mu_o3a"]
     normalize = np.nanmax(np.abs(gdiff))
     export_img("convergence_o3a_vb7_diff", gdiff/normalize)
     print(f"Color scale for o3a_vb7 difference: {np.pi*normalize}")
+    print(f"\\newcommand\\scalevbMaskedDiff{{{np.pi*normalize}}}% convergence_o3a_vb7_diff.png", file=file)
     grel = gdiff/data_vb7["gdc_mu_o3a"]
     normalize = np.nanmax(np.abs(grel))
     export_img("convergence_o3a_vb7_relative", grel/normalize)
     print(f"Color scale for o3a_vb7 relative: {normalize}")
+    print(f"\\newcommand\\scalevbMaskedRel{{{normalize}}}% convergence_o3a_vb7_relative.png", file=file)
 
     for method, u  in (("o2","mu"), ("o3","mu"), ("o3a","mu"), ("o3p","mu"), ("o3a_ispline","mu")):
         if u != "mu":
@@ -964,6 +982,8 @@ def prepare_RGconvergence():
     print("DC voltage: " + " ".join("%.6g"%vdc[i] for i in dc_indices))
     print("AC voltage: " + " ".join("%.6g"%vac[i] for i in ac_indices))
 
+    file.close()
+
 
 def adjust_argument_array(start, end, *funcs, initial_resolution=51, s1=0.01, maxit=5, scale_coef=-0.5, yscales=1):
     x = np.linspace(start, end, initial_resolution)
diff --git a/kondo.py b/kondo.py
index 69e8663..8d785b6 100644
--- a/kondo.py
+++ b/kondo.py
@@ -535,7 +535,7 @@ class Kondo:
             assert not solve_integral_exactly
         if solve_integral_exactly:
             if padding > 0:
-                settings.logger.warn("solve_integral_exactly and padding>0: will lead not handle padding consistenly")
+                settings.logger.warn("solve_integral_exactly and padding>0: will not handle padding consistenly")
             if integral_method == -15:
                 self.integral_method = -1
                 settings.logger.warn("setting integral_method = -1")
diff --git a/plot.py b/plot.py
index a9a3ba3..10a0447 100644
--- a/plot.py
+++ b/plot.py
@@ -2456,5 +2456,32 @@ def kondo_temperature(dm, **parameters):
     print(f"G/G0 = {np.pi*shift:.5g} + {np.pi*prefactor:.5g}/(1 + {cv:.5g} (Vdc/Tkrg)^2)")
 
 
+def compare_resonant(dm, **parameters):
+    if parameters.get("omega", None) is None:
+        parameters["omega"] = 16.5372
+    parameters.pop("voltage_branches", None)
+    parameters.pop("resonant_dc_shift", None)
+    data = dm.list(**parameters)
+    data = data.loc[(((data.vdc + 1e-4)%data.omega) < 2e-4) & (data.vdc > 0)]
+    data.omega = np.round(data.omega, 8)
+    data.vdc = np.round(data.vdc, 8)
+    data.vac = np.round(data.vac, 8)
+    data_resonant = data.loc[(data.method == "J") & (data.padding >= (data.nmax>0)) & (data.resonant_dc_shift > 0) & (data.voltage_branches == 0)]
+    data_ref = data.loc[(data.method == "mu") & (data.padding == 0) & (data.resonant_dc_shift == 0) & (data.voltage_branches == 4)]
+    merged = pd.merge(
+            data_resonant,
+            data_ref,
+            how="inner",
+            on=["vdc","vac","d","omega","energy_im","energy_re","solver_tol_rel","solver_tol_abs","xL","lazy_inverse_factor"],
+            suffixes=("_res", "_ref"))
+    print(data_resonant.shape, data_ref.shape, merged.shape)
+    diff_rel = (merged.dc_conductance_res - merged.dc_conductance_ref)/merged.dc_conductance_res
+    diff_rel_max = np.nanmax(np.abs(diff_rel[merged.vac > 0]))
+    norm = plt.Normalize(-diff_rel_max, diff_rel_max)
+    s = plt.scatter(merged.vdc, merged.vac, c=diff_rel, cmap=plt.cm.seismic, norm=norm)
+    plt.colorbar(s)
+    plt.show()
+
+
 if __name__ == '__main__':
     main()
diff --git a/tikz/asymmetry.tex b/tikz/asymmetry.tex
index 441b27e..9e61920 100644
--- a/tikz/asymmetry.tex
+++ b/tikz/asymmetry.tex
@@ -1,5 +1,4 @@
-%\tikzsetnextfilename{asymmetry}%
-\definecolor{darkgreen}{HTML}{00c000}
+\tikzsetnextfilename{asymmetry}%
 \begin{tikzpicture}
   \begin{axis}[
       name = GvsVdc,
@@ -11,6 +10,7 @@
       xlabel = {$\vdc~(\Omega)$},
       ylabel = {$\gdc/(4x_Lx_R)\quad(2e^2/h)$},
       scale only axis,
+      axis on top,
       every axis plot post/.append style = {
         line width = 0.6pt,
         line cap = round,
@@ -27,12 +27,12 @@
         $x_R=0.5$,
       },
     ]
-    \addplot[red] table [x=vdc, y={xL0.001}] {../figdata/asymmetry.dat};
-    \addplot[violet] table [x=vdc, y={xL0.1}] {../figdata/asymmetry.dat};
-    \addplot[blue] table [x=vdc, y={xL0.2}] {../figdata/asymmetry.dat};
-    \addplot[darkgreen] table [x=vdc, y={xL0.3}] {../figdata/asymmetry.dat};
-    \addplot[orange] table [x=vdc, y={xL0.4}] {../figdata/asymmetry.dat};
-    \addplot[black] table [x=vdc, y={xL0.5}] {../figdata/asymmetry.dat};
+    \addplot[color=c24] table [x=vdc, y={xL0.001}] {../figdata/asymmetry.dat};
+    \addplot[color=c20] table [x=vdc, y={xL0.1}] {../figdata/asymmetry.dat};
+    \addplot[color=c14] table [x=vdc, y={xL0.2}] {../figdata/asymmetry.dat};
+    \addplot[color=c10] table [x=vdc, y={xL0.3}] {../figdata/asymmetry.dat};
+    \addplot[color=c06] table [x=vdc, y={xL0.4}] {../figdata/asymmetry.dat};
+    \addplot[color=black] table [x=vdc, y={xL0.5}] {../figdata/asymmetry.dat};
     \node[black,anchor=south west] at (axis description cs:0.04,0.04) {\parbox{5.9em}{\raggedright\noindent$\vac=1.6\Omega$\\$\hphantom{\vac}\mathllap{\Omega}=4.8\tkv$}};
   \end{axis}
 \end{tikzpicture}
diff --git a/tikz/bruhat18.tex b/tikz/bruhat18.tex
index 8c8ab18..832ced2 100644
--- a/tikz/bruhat18.tex
+++ b/tikz/bruhat18.tex
@@ -1,32 +1,5 @@
-% Configure figures
-\definecolor{frtrg}{rgb}{.95,0,0}%
-\definecolor{adiabatic}{rgb}{0,0,0}%
-\definecolor{phenomenological}{rgb}{0,0,1}%
-\definecolor{analytic}{rgb}{0,.6,0}%
-\definecolor{experiment}{rgb}{0,0,0}%
-\definecolor{ac}{rgb}{.5,.5,.5}%
-\definecolor{green}{rgb}{0,.67,0}%
-\pgfplotsset{adiabatic style/.style={densely dotted, thick, color=adiabatic}}%
-\pgfplotsset{frtrg style/.style={solid, line width=0.75pt, color=frtrg, no marks}}%
-\pgfplotsset{phenomenological style/.style={color=phenomenological}}%
-\pgfplotsset{ac style/.style={thick, dashdotted, color=ac}}%
-\pgfplotsset{analytic style/.style={thick, densely dashed}}%
-\definecolor{c02}{rgb}{.5,.5,.5}%
-\definecolor{c04}{rgb}{.7,.25,.35}%
-\definecolor{c06}{rgb}{.95,.1,0}%
-\definecolor{c08}{rgb}{.9,.5,0}%
-\definecolor{c10}{rgb}{.67,.54,0}%
-\definecolor{c12}{rgb}{.5,.75,0}%
-\definecolor{c14}{rgb}{0,.62,0}%
-\definecolor{c16}{rgb}{0,.8,.5}%
-\definecolor{c18}{rgb}{0,.6,.7}%
-\definecolor{c20}{rgb}{.25,.3,1}%
-\definecolor{c22}{rgb}{.6,.35,.9}%
-\definecolor{c24}{rgb}{.9,0,1}%
-\definecolor{c30}{rgb}{.8,.25,.6}%
-\definecolor{lightgray}{rgb}{.6,.6,.6}%
-\newcommand\plotsep{0.04}%
-%\tikzsetnextfilename{bruhat18}%
+\def\plotsep{0.04}%
+\tikzsetnextfilename{bruhat18}%
 \begin{tikzpicture}
   \begin{axis}[
       name = 19GHz,
@@ -42,6 +15,7 @@
       clip marker paths = true,
       xtick distance = 2,
       minor x tick num = 1,
+      axis on top,
       legend entries = {
         {\clap{$\vac~(\mu V)$}}, {\null},
         20, 40, 60, 80, 100, 120, 140, 160, 180, 200, 220, 240, {\null}, 300
@@ -68,74 +42,74 @@
     \addlegendimage{empty legend};
     \addlegendimage{empty legend};
     %
-    \addplot[frtrg style, on layer=main, color=c02] coordinates {(9,0) (10,0)};
-    \addplot[frtrg style, on layer=main, color=c04] coordinates {(9,0) (10,0)};
-    \addplot[frtrg style, on layer=main, color=c06] coordinates {(9,0) (10,0)};
+    \addplot[exp style, on layer=main, color=c02] coordinates {(9,0) (10,0)};
+    \addplot[exp style, on layer=main, color=c04] coordinates {(9,0) (10,0)};
+    \addplot[exp style, on layer=main, color=c06] coordinates {(9,0) (10,0)};
     %
     \addplot[only marks, on layer=main, color=c08, forget plot] table[x expr={12.7263*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac80}  + 9*\plotsep}] {../exp_data/KondoAC_Freq19GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g0.906} + 9*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c08] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g0.906} + 9*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g0.906} + 9*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, on layer=main, color=c08] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g0.906} + 9*\plotsep}] {../figdata/bruhat18_19GHz.dat};
     \addplot[only marks, on layer=main, color=c10, forget plot] table[x expr={12.7263*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac100} + 8*\plotsep}] {../exp_data/KondoAC_Freq19GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.13}  + 8*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c10] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.13}  + 8*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.13}  + 8*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, on layer=main, color=c10] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.13}  + 8*\plotsep}] {../figdata/bruhat18_19GHz.dat};
     \addplot[only marks, on layer=main, color=c12, forget plot] table[x expr={12.7263*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac120} + 7*\plotsep}] {../exp_data/KondoAC_Freq19GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.36}  + 7*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c12] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.36}  + 7*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.36}  + 7*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, on layer=main, color=c12] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.36}  + 7*\plotsep}] {../figdata/bruhat18_19GHz.dat};
     \addplot[only marks, on layer=main, color=c14, forget plot] table[x expr={12.7263*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac140} + 6*\plotsep}] {../exp_data/KondoAC_Freq19GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.59}  + 6*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c14] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.59}  + 6*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.59}  + 6*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, on layer=main, color=c14] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.59}  + 6*\plotsep}] {../figdata/bruhat18_19GHz.dat};
     \addplot[only marks, on layer=main, color=c16, forget plot] table[x expr={12.7263*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac160} + 5*\plotsep}] {../exp_data/KondoAC_Freq19GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.81}  + 5*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c16] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.81}  + 5*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.81}  + 5*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, on layer=main, color=c16] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.81}  + 5*\plotsep}] {../figdata/bruhat18_19GHz.dat};
     \addplot[only marks, on layer=main, color=c18, forget plot] table[x expr={12.7263*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac180} + 4*\plotsep}] {../exp_data/KondoAC_Freq19GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.04}  + 4*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c18] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.04}  + 4*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.04}  + 4*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, on layer=main, color=c18] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.04}  + 4*\plotsep}] {../figdata/bruhat18_19GHz.dat};
     \addplot[only marks, on layer=main, color=c20, forget plot] table[x expr={12.7263*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac200} + 3*\plotsep}] {../exp_data/KondoAC_Freq19GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.27}  + 3*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c20] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.27}  + 3*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.27}  + 3*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, on layer=main, color=c20] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.27}  + 3*\plotsep}] {../figdata/bruhat18_19GHz.dat};
     \addplot[only marks, on layer=main, color=c22, forget plot] table[x expr={12.7263*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac220} + 2*\plotsep}] {../exp_data/KondoAC_Freq19GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.49}  + 2*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c22] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.49}  + 2*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.49}  + 2*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, on layer=main, color=c22] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.49}  + 2*\plotsep}] {../figdata/bruhat18_19GHz.dat};
     \addplot[only marks, on layer=main, color=c24, forget plot] table[x expr={12.7263*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac240} + 1*\plotsep}] {../exp_data/KondoAC_Freq19GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.72}  + 1*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c24] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.72}  + 1*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.72}  + 1*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, on layer=main, color=c24] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.72}  + 1*\plotsep}] {../figdata/bruhat18_19GHz.dat};
     \addlegendimage{empty legend};
     \addplot[only marks, on layer=main, color=c30, forget plot] table[x expr={12.7263*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac300} + 0*\plotsep}] {../exp_data/KondoAC_Freq19GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g3.4}   + 0*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c30] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g3.4}   + 0*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g3.4}   + 0*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, on layer=main, color=c30] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g3.4}   + 0*\plotsep}] {../figdata/bruhat18_19GHz.dat};
     %
-    \addplot[frtrg style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg0.906} + 9*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg1.13}  + 8*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg1.36}  + 7*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg1.59}  + 6*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg1.81}  + 5*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg2.04}  + 4*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg2.27}  + 3*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg2.49}  + 2*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg2.72}  + 1*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg3.4}   + 0*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg0.906} + 9*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg1.13}  + 8*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg1.36}  + 7*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg1.59}  + 6*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg1.81}  + 5*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg2.04}  + 4*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg2.27}  + 3*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg2.49}  + 2*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg2.72}  + 1*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg3.4}   + 0*\plotsep}] {../figdata/bruhat18_19GHz.dat};
     %
-    \addplot[frtrg style, color=black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg0.906} + 9*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, color=black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg1.13}  + 8*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, color=black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg1.36}  + 7*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, color=black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg1.59}  + 6*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, color=black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg1.81}  + 5*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, color=black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg2.04}  + 4*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, color=black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg2.27}  + 3*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, color=black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg2.49}  + 2*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, color=black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg2.72}  + 1*\plotsep}] {../figdata/bruhat18_19GHz.dat};
-    \addplot[frtrg style, color=black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg3.4}   + 0*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, color=black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg0.906} + 9*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, color=black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg1.13}  + 8*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, color=black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg1.36}  + 7*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, color=black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg1.59}  + 6*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, color=black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg1.81}  + 5*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, color=black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg2.04}  + 4*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, color=black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg2.27}  + 3*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, color=black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg2.49}  + 2*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, color=black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg2.72}  + 1*\plotsep}] {../figdata/bruhat18_19GHz.dat};
+    \addplot[exp style, color=black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg3.4}   + 0*\plotsep}] {../figdata/bruhat18_19GHz.dat};
     %
-    %\addplot[frtrg style, cyan, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g0.906} + 9*\plotsep}] {../figdata/bruhat18_19GHz_asym.dat};
-    %\addplot[frtrg style, cyan, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.13}  + 8*\plotsep}] {../figdata/bruhat18_19GHz_asym.dat};
-    %\addplot[frtrg style, cyan, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.36}  + 7*\plotsep}] {../figdata/bruhat18_19GHz_asym.dat};
-    %\addplot[frtrg style, cyan, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.59}  + 6*\plotsep}] {../figdata/bruhat18_19GHz_asym.dat};
-    %\addplot[frtrg style, cyan, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.81}  + 5*\plotsep}] {../figdata/bruhat18_19GHz_asym.dat};
-    %\addplot[frtrg style, cyan, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.04}  + 4*\plotsep}] {../figdata/bruhat18_19GHz_asym.dat};
-    %\addplot[frtrg style, cyan, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.27}  + 3*\plotsep}] {../figdata/bruhat18_19GHz_asym.dat};
-    %\addplot[frtrg style, cyan, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.49}  + 2*\plotsep}] {../figdata/bruhat18_19GHz_asym.dat};
-    %\addplot[frtrg style, cyan, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.72}  + 1*\plotsep}] {../figdata/bruhat18_19GHz_asym.dat};
-    %\addplot[frtrg style, cyan, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g3.4}   + 0*\plotsep}] {../figdata/bruhat18_19GHz_asym.dat};
+    %\addplot[exp style, cyan, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g0.906} + 9*\plotsep}] {../figdata/bruhat18_19GHz_asym.dat};
+    %\addplot[exp style, cyan, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.13}  + 8*\plotsep}] {../figdata/bruhat18_19GHz_asym.dat};
+    %\addplot[exp style, cyan, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.36}  + 7*\plotsep}] {../figdata/bruhat18_19GHz_asym.dat};
+    %\addplot[exp style, cyan, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.59}  + 6*\plotsep}] {../figdata/bruhat18_19GHz_asym.dat};
+    %\addplot[exp style, cyan, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.81}  + 5*\plotsep}] {../figdata/bruhat18_19GHz_asym.dat};
+    %\addplot[exp style, cyan, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.04}  + 4*\plotsep}] {../figdata/bruhat18_19GHz_asym.dat};
+    %\addplot[exp style, cyan, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.27}  + 3*\plotsep}] {../figdata/bruhat18_19GHz_asym.dat};
+    %\addplot[exp style, cyan, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.49}  + 2*\plotsep}] {../figdata/bruhat18_19GHz_asym.dat};
+    %\addplot[exp style, cyan, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.72}  + 1*\plotsep}] {../figdata/bruhat18_19GHz_asym.dat};
+    %\addplot[exp style, cyan, on layer=main, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g3.4}   + 0*\plotsep}] {../figdata/bruhat18_19GHz_asym.dat};
   \end{axis}
   \begin{axis}[
       name = 12GHz,
@@ -154,67 +128,68 @@
       yticklabels = {},
       mark size = 0.4pt,
       clip marker paths = true,
+      axis on top,
     ]
     \node[black] at (axis description cs:0.25,0.9) {$f=12\,\mathrm{GHz}$};
     %
     \addplot[only marks, on layer=main, color=c02, forget plot] table[x expr={20.1499*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac20} + 10*\plotsep}] {../exp_data/KondoAC_Freq12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g0.431} +10*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c02] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g0.431} +10*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g0.431} +10*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=c02] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g0.431} +10*\plotsep}] {../figdata/bruhat18_12GHz.dat};
     \addplot[only marks, on layer=main, color=c04, forget plot] table[x expr={20.1499*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac40}  + 9*\plotsep}] {../exp_data/KondoAC_Freq12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g0.862} + 9*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c04] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g0.862} + 9*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g0.862} + 9*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=c04] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g0.862} + 9*\plotsep}] {../figdata/bruhat18_12GHz.dat};
     \addplot[only marks, on layer=main, color=c06, forget plot] table[x expr={20.1499*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac60}  + 8*\plotsep}] {../exp_data/KondoAC_Freq12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.29}  + 8*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c06] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.29}  + 8*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.29}  + 8*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=c06] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.29}  + 8*\plotsep}] {../figdata/bruhat18_12GHz.dat};
     \addplot[only marks, on layer=main, color=c08, forget plot] table[x expr={20.1499*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac80}  + 7*\plotsep}] {../exp_data/KondoAC_Freq12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.72}  + 7*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c08] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.72}  + 7*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.72}  + 7*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=c08] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g1.72}  + 7*\plotsep}] {../figdata/bruhat18_12GHz.dat};
     \addplot[only marks, on layer=main, color=c10, forget plot] table[x expr={20.1499*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac100} + 6*\plotsep}] {../exp_data/KondoAC_Freq12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.16}  + 6*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c10] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.16}  + 6*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.16}  + 6*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=c10] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.16}  + 6*\plotsep}] {../figdata/bruhat18_12GHz.dat};
     \addplot[only marks, on layer=main, color=c12, forget plot] table[x expr={20.1499*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac120} + 5*\plotsep}] {../exp_data/KondoAC_Freq12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.59}  + 5*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c12] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.59}  + 5*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.59}  + 5*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=c12] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g2.59}  + 5*\plotsep}] {../figdata/bruhat18_12GHz.dat};
     \addplot[only marks, on layer=main, color=c14, forget plot] table[x expr={20.1499*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac140} + 4*\plotsep}] {../exp_data/KondoAC_Freq12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g3.02}  + 4*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c14] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g3.02}  + 4*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g3.02}  + 4*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=c14] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g3.02}  + 4*\plotsep}] {../figdata/bruhat18_12GHz.dat};
     \addplot[only marks, on layer=main, color=c16, forget plot] table[x expr={20.1499*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac160} + 3*\plotsep}] {../exp_data/KondoAC_Freq12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g3.45}  + 3*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c16] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g3.45}  + 3*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g3.45}  + 3*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=c16] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g3.45}  + 3*\plotsep}] {../figdata/bruhat18_12GHz.dat};
     \addplot[only marks, on layer=main, color=c18, forget plot] table[x expr={20.1499*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac180} + 2*\plotsep}] {../exp_data/KondoAC_Freq12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g3.88}  + 2*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c18] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g3.88}  + 2*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    %\addplot[frtrg style, on layer=main, color=c20] coordinates {(9,0) (10,0)};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g3.88}  + 2*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=c18] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g3.88}  + 2*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    %\addplot[exp style, on layer=main, color=c20] coordinates {(9,0) (10,0)};
     \addplot[only marks, on layer=main, color=c22, forget plot] table[x expr={20.1499*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac220} + 1*\plotsep}] {../exp_data/KondoAC_Freq12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g4.74}  + 1*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c22] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g4.74}  + 1*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    %\addplot[frtrg style, on layer=main, color=c24] coordinates {(9,0) (10,0)};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g4.74}  + 1*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=c22] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g4.74}  + 1*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    %\addplot[exp style, on layer=main, color=c24] coordinates {(9,0) (10,0)};
     \addplot[only marks, on layer=main, color=c30, forget plot] table[x expr={20.1499*\thisrow{Vsd_mV}}, y expr={\thisrow{Vac300} + 0*\plotsep}] {../exp_data/KondoAC_Freq12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g6.47}  + 0*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=c30] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g6.47}  + 0*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g6.47}  + 0*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=c30] table[x expr={\thisrow{vdc}}, y expr={\thisrow{g6.47}  + 0*\plotsep}] {../figdata/bruhat18_12GHz.dat};
     %
-    \addplot[frtrg style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg0.431} +10*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg0.862} + 9*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg1.29}  + 8*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg1.72}  + 7*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg2.16}  + 6*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg2.59}  + 5*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg3.02}  + 4*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg3.45}  + 3*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg3.88}  + 2*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg4.74}  + 1*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg6.47}  + 0*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg0.431} +10*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg0.862} + 9*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg1.29}  + 8*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg1.72}  + 7*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg2.16}  + 6*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg2.59}  + 5*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg3.02}  + 4*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg3.45}  + 3*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg3.88}  + 2*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg4.74}  + 1*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=lightgray, line width=0.4pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gpatbg6.47}  + 0*\plotsep}] {../figdata/bruhat18_12GHz.dat};
     %
-    \addplot[frtrg style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg0.431} +10*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg0.862} + 9*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg1.29}  + 8*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg1.72}  + 7*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg2.16}  + 6*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg2.59}  + 5*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg3.02}  + 4*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg3.45}  + 3*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg3.88}  + 2*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg4.74}  + 1*\plotsep}] {../figdata/bruhat18_12GHz.dat};
-    \addplot[frtrg style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg6.47}  + 0*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg0.431} +10*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg0.862} + 9*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg1.29}  + 8*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg1.72}  + 7*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg2.16}  + 6*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg2.59}  + 5*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg3.02}  + 4*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg3.45}  + 3*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg3.88}  + 2*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg4.74}  + 1*\plotsep}] {../figdata/bruhat18_12GHz.dat};
+    \addplot[exp style, on layer=main, color=black, line width=0.6pt, forget plot] table[x expr={\thisrow{vdc}}, y expr={\thisrow{gbg6.47}  + 0*\plotsep}] {../figdata/bruhat18_12GHz.dat};
   \end{axis}
 \end{tikzpicture}
diff --git a/tikz/commutator_RG2.tex b/tikz/commutator_RG2.tex
index d53708a..1d7cbe9 100644
--- a/tikz/commutator_RG2.tex
+++ b/tikz/commutator_RG2.tex
@@ -1,3 +1,4 @@
+\tikzsetnextfilename{commutator_RG2}%
 \begin{tikzpicture}
   \pgfplotsset{set layers = axis on top}
   \begin{axis}[
diff --git a/tikz/contour.tex b/tikz/contour.tex
index a8ecaeb..7cdc29e 100644
--- a/tikz/contour.tex
+++ b/tikz/contour.tex
@@ -1,32 +1,33 @@
-\definecolor{darkgreen}{HTML}{00c000}
+\tikzsetnextfilename{contour}%
 \begin{tikzpicture}
   \begin{axis}[
-      width = 12cm,
-      height = 10cm,
+      width = 13cm,
+      height = 11cm,
       xmin = 0, xmax = 10,
       ymin = 0, ymax = 10,
       xlabel = {$\vdc~(\Omega)$},
       ylabel = {$\vac~(\Omega)$},
+      axis on top,
       every axis plot post/.append style = {
-        line width=0.2pt,
+        line width=0.4pt,
         line cap=round,
         line join=round
       },
       legend cell align = left,
       legend style = {
         anchor = south east,
-        at = {(0.99,0.01)}
+        at = {(0.975,0.025)}
       }
     ]
-    \addplot[darkgreen] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_mu_o2.dat};
-    \addplot[darkgreen, forget plot] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_J_o2.dat};
-    \addplot[darkgreen, 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[c20] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_mu_o2.dat};
+    \addplot[c20, forget plot] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_J_o2.dat};
+    \addplot[c20, forget plot] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_J_o2_p.dat};
+    \addplot[c14] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_mu_o3.dat};
+    \addplot[c14, forget plot] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_J_o3.dat};
+    \addplot[c14, forget plot] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_J_o3_p.dat};
+    \addplot[c06] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_mu_o3a.dat};
+    \addplot[c06, forget plot] table[x=vdc, y=vac, z=gdc]{../figdata/contour_gdc_J_o3a.dat};
+    \addplot[c06, 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};
diff --git a/tikz/convergence-j.tex b/tikz/convergence-j.tex
index 8456c84..770bed7 100644
--- a/tikz/convergence-j.tex
+++ b/tikz/convergence-j.tex
@@ -1,3 +1,4 @@
+\tikzsetnextfilename{convergence_J}%
 \begin{tikzpicture}
   \pgfplotsset{set layers = axis on top}
   \begin{axis}[
@@ -31,15 +32,17 @@
       yticklabel = \empty,
       clip = true,
     ]
-    \addplot graphics[xmin=-.03333333, xmax=10.03333333, ymin=-.05, ymax=10.05] {../figdata/convergence_o3a_J_p_relative_fromdiff_masked.png};
+    \addplot graphics[xmin=\extentxmin, xmax=\extentxmax, ymin=\extentymin, ymax=\extentymax] {../figdata/convergence_o3a_J_p_relative_fromdiff_masked.png};
     \node[anchor=north west] at (axis description cs:0.05,0.95) {(b)};
   \end{axis}
   \begin{axis}[
       name = colorbar,
       at = (Jp.east),
       anchor = west,
-      ymin = -0.0004348761377387133,
-      ymax = 0.0004348761377387133,
+      %ymin = -0.0004348761377387133,
+      %ymax = 0.0004348761377387133,
+      ymin = -\scaleoxaxJrelX,
+      ymax = \scaleoxaxJrelX,
       xshift = 8mm,
       enlargelimits = false,
       axis on top,
@@ -57,6 +60,6 @@
       ylabel = {relative deviation},
       %ylabel shift = -5mm,
     ]
-    \addplot graphics[xmin=0, xmax=1, ymin=-0.0004348761377387133, ymax=0.0004348761377387133] {../figdata/colorbar_seismic_vertical.png};
+    \addplot graphics[xmin=0, xmax=1, ymin=-\scaleoxaxJrelX, ymax=\scaleoxaxJrelX] {../figdata/colorbar_seismic_vertical.png};
   \end{axis}
 \end{tikzpicture}
diff --git a/tikz/convergence-o2.tex b/tikz/convergence-o2.tex
index 0b96d0a..a2d1699 100644
--- a/tikz/convergence-o2.tex
+++ b/tikz/convergence-o2.tex
@@ -1,3 +1,4 @@
+\tikzsetnextfilename{convergence_o2}%
 \begin{tikzpicture}
   \pgfplotsset{set layers = axis on top}
   \begin{axis}[
@@ -17,7 +18,7 @@
       },
       clip = false,
     ]
-    \addplot graphics[xmin=-.03333333, xmax=10.03333333, ymin=-.05, ymax=10.05] {../figdata/convergence_o2_diff.png};
+    \addplot graphics[xmin=\extentxmin, xmax=\extentxmax, ymin=\extentymin, ymax=\extentymax] {../figdata/convergence_o2_diff.png};
     \node at (2,9.2) {(b)};
     % vertical lines
     \addplot[no marks, black] coordinates {(0,1.6) (10,1.6)};
@@ -133,8 +134,10 @@
   \begin{axis}[
       name = colorbar,
       at = (dc.north west),
-      xmin = -0.0314149933,
-      xmax = 0.0314149933,
+      %xmin = -0.0314149933,
+      %xmax = 0.0314149933,
+      xmin = -\scaleoldiff,
+      xmax = \scaleoldiff,
       xshift = 5mm,
       yshift = 11mm,
       anchor = south west,
@@ -152,6 +155,6 @@
       xlabel = {deviation $(2e^2/h)$},
       xlabel shift = -15mm,
     ]
-    \addplot graphics[ymin=0, ymax=1, xmin=-0.0314149933, xmax=0.0314149933] {../figdata/colorbar_seismic.png};
+    \addplot graphics[ymin=0, ymax=1, xmin=-\scaleoldiff, xmax=\scaleoldiff] {../figdata/colorbar_seismic.png};
   \end{axis}
 \end{tikzpicture}
diff --git a/tikz/convergence-o3.tex b/tikz/convergence-o3.tex
index 280378b..bf7d672 100644
--- a/tikz/convergence-o3.tex
+++ b/tikz/convergence-o3.tex
@@ -1,3 +1,4 @@
+\tikzsetnextfilename{convergence_o3}%
 \begin{tikzpicture}
   \pgfplotsset{set layers = axis on top}
   \begin{axis}[
@@ -17,7 +18,7 @@
       },
       clip = false,
     ]
-    \addplot graphics[xmin=-.03333333, xmax=10.03333333, ymin=-.05, ymax=10.05] {../figdata/convergence_o3_diff.png};
+    \addplot graphics[xmin=\extentxmin, xmax=\extentxmax, ymin=\extentymin, ymax=\extentymax] {../figdata/convergence_o3_diff.png};
     \node at (2,9.2) {(b)};
     % vertical lines
     \addplot[no marks, black] coordinates {(0,1.6) (10,1.6)};
@@ -133,8 +134,10 @@
   \begin{axis}[
       name = colorbar,
       at = (dc.north west),
-      xmin = -0.02863684157,
-      xmax = 0.02863684157,
+      %xmin = -0.02863684157,
+      %xmax = 0.02863684157,
+      xmin = -\scaleoxdiff,
+      xmax = \scaleoxdiff,
       xshift = 5mm,
       yshift = 11mm,
       anchor = south west,
@@ -152,6 +155,6 @@
       xlabel = {deviation $(2e^2/h)$},
       xlabel shift = -15mm,
     ]
-    \addplot graphics[ymin=0, ymax=1, xmin=-0.02863684157, xmax=0.02863684157] {../figdata/colorbar_seismic.png};
+    \addplot graphics[ymin=0, ymax=1, xmin=-\scaleoxdiff, xmax=\scaleoxdiff] {../figdata/colorbar_seismic.png};
   \end{axis}
 \end{tikzpicture}
diff --git a/tikz/convergence-o3a.tex b/tikz/convergence-o3a.tex
index 789f577..90ce622 100644
--- a/tikz/convergence-o3a.tex
+++ b/tikz/convergence-o3a.tex
@@ -1,3 +1,4 @@
+\tikzsetnextfilename{convergence_o3a}%
 \begin{tikzpicture}
   \pgfplotsset{set layers = axis on top}
   \begin{axis}[
@@ -17,7 +18,7 @@
       },
       clip = false,
     ]
-    \addplot graphics[xmin=-.03333333, xmax=10.03333333, ymin=-.05, ymax=10.05] {../figdata/convergence_o3p_diff_fromdiff.png};
+    \addplot graphics[xmin=\extentxmin, xmax=\extentxmax, ymin=\extentymin, ymax=\extentymax] {../figdata/convergence_o3p_diff_fromdiff.png};
     \node at (2,9.2) {(b)};
     % vertical lines
     \addplot[no marks, black] coordinates {(0,1.6) (10,1.6)};
@@ -133,8 +134,10 @@
   \begin{axis}[
       name = colorbar,
       at = (dc.north west),
-      xmin = -0.002963806506594897,
-      xmax = 0.002963806506594897,
+      %xmin = -0.002963806506594897,
+      %xmax = 0.002963806506594897,
+      xmin = -\scaleoxpdiffX,
+      xmax = \scaleoxpdiffX,
       xshift = 5mm,
       yshift = 11mm,
       anchor = south west,
@@ -151,6 +154,6 @@
       xlabel = {deviation $(2e^2/h)$},
       xlabel shift = -15.85mm,
     ]
-    \addplot graphics[ymin=0, ymax=1, xmin=-0.002963806506594897, xmax=0.002963806506594897] {../figdata/colorbar_seismic.png};
+    \addplot graphics[ymin=0, ymax=1, xmin=-\scaleoxpdiffX, xmax=\scaleoxpdiffX] {../figdata/colorbar_seismic.png};
   \end{axis}
 \end{tikzpicture}
diff --git a/tikz/convergence-o3a_ispline.tex b/tikz/convergence-o3a_ispline.tex
index fdf693a..53542c1 100644
--- a/tikz/convergence-o3a_ispline.tex
+++ b/tikz/convergence-o3a_ispline.tex
@@ -1,3 +1,4 @@
+\tikzsetnextfilename{convergence_o3a_ispline}%
 \begin{tikzpicture}
   \pgfplotsset{set layers = axis on top}
   \begin{axis}[
@@ -17,7 +18,7 @@
       },
       clip = false,
     ]
-    \addplot graphics[xmin=-.03333333, xmax=10.03333333, ymin=-.05, ymax=10.05] {../figdata/convergence_o3a_ispline_relative.png};
+    \addplot graphics[xmin=\extentxmin, xmax=\extentxmax, ymin=\extentymin, ymax=\extentymax] {../figdata/convergence_o3a_ispline_relative.png};
     \node at (2,9.2) {(b)};
     % vertical lines
     \addplot[no marks, black] coordinates {(0,1.6) (10,1.6)};
@@ -133,8 +134,10 @@
   \begin{axis}[
       name = colorbar,
       at = (dc.north west),
-      xmin = -0.06062738883,
-      xmax = 0.06062738883,
+      %xmin = -0.06062738883,
+      %xmax = 0.06062738883,
+      xmin = -\scaleoxaxisplinerel,
+      xmax = \scaleoxaxisplinerel,
       xshift = 5mm,
       yshift = 11mm,
       anchor = south west,
@@ -152,6 +155,6 @@
       xlabel = {relative deviation},
       xlabel shift = -14mm,
     ]
-    \addplot graphics[ymin=0, ymax=1, xmin=-0.06062738883, xmax=0.06062738883] {../figdata/colorbar_seismic.png};
+    \addplot graphics[ymin=0, ymax=1, xmin=-\scaleoxaxisplinerel, xmax=\scaleoxaxisplinerel] {../figdata/colorbar_seismic.png};
   \end{axis}
 \end{tikzpicture}
diff --git a/tikz/current_convergence.tex b/tikz/current_convergence.tex
index cbbaa54..b4acc9e 100644
--- a/tikz/current_convergence.tex
+++ b/tikz/current_convergence.tex
@@ -1,3 +1,4 @@
+\tikzsetnextfilename{current_convergence}%
 \begin{tikzpicture}
   \begin{axis}[
       name = idc,
diff --git a/tikz/current_convergence_bad.tex b/tikz/current_convergence_bad.tex
index 26f1a60..7c554a5 100644
--- a/tikz/current_convergence_bad.tex
+++ b/tikz/current_convergence_bad.tex
@@ -1,3 +1,4 @@
+\tikzsetnextfilename{current_convergence_bad}%
 \begin{tikzpicture}
   \begin{axis}[
       name = idc,
diff --git a/tikz/floquet_matrices.tex b/tikz/floquet_matrices.tex
index f05d5a1..b6b1c3a 100644
--- a/tikz/floquet_matrices.tex
+++ b/tikz/floquet_matrices.tex
@@ -1,4 +1,5 @@
-\newcommand\matrixsize{24mm}
+\tikzsetnextfilename{floquet_matrices}%
+\newcommand\matrixsize{24mm}%
 \begin{tikzpicture}
   \begin{axis}[
       name = z_mu_small,
diff --git a/tikz/g_approximations.tex b/tikz/g_approximations.tex
index 0472313..83d256a 100644
--- a/tikz/g_approximations.tex
+++ b/tikz/g_approximations.tex
@@ -1,17 +1,4 @@
-% Configure figures
-\definecolor{frtrg}{rgb}{.95,0,0}%
-\definecolor{adiabatic}{rgb}{0,0,0}%
-\definecolor{phenomenological}{rgb}{0,0,1}%
-\definecolor{analytic}{rgb}{0,.6,0}%
-\definecolor{experiment}{rgb}{0,0,0}%
-\definecolor{ac}{rgb}{.5,.5,.5}%
-\definecolor{green}{rgb}{0,.67,0}%
-\pgfplotsset{adiabatic style/.style={densely dotted, thick, color=adiabatic}}%
-\pgfplotsset{frtrg style/.style={solid, thick, color=frtrg}}%
-\pgfplotsset{phenomenological style/.style={color=phenomenological}}%
-\pgfplotsset{ac style/.style={thick, dashdotted, color=ac}}%
-\pgfplotsset{analytic style/.style={thick, densely dashed}}%
-%\tikzsetnextfilename{g_approximations}%
+\tikzsetnextfilename{g_approximations}%
 \begin{tikzpicture}
   \begin{axis}[
       name = GvsVdc,
diff --git a/tikz/ga.tex b/tikz/ga.tex
index ee7ed0e..99d19af 100644
--- a/tikz/ga.tex
+++ b/tikz/ga.tex
@@ -1,3 +1,4 @@
+\tikzsetnextfilename{ga}%
 \begin{tikzpicture}
   \pgfplotsset{set layers = axis on top}
   \begin{axis}[
@@ -13,7 +14,7 @@
       ylabel = {$\vac~(\Omega)$},
       clip = true,
     ]
-    \addplot graphics[xmin=-.03333333, xmax=10.03333333, ymin=-.05, ymax=10.05] {../figdata/Ga_max_abs.png};
+    \addplot graphics[xmin=-0.0166667, xmax=10.0167, ymin=-0.025, ymax=10.025] {../figdata/Ga_max_abs.png};
     \node[white,anchor=north west] at (axis description cs:0.05,0.95) {(a)};
   \end{axis}
   \begin{axis}[
@@ -55,7 +56,7 @@
       ylabel = {$\vac~(\Omega)$},
       clip = true,
     ]
-    \addplot graphics[xmin=-.03333333, xmax=10.03333333, ymin=-.05, ymax=10.05] {../figdata/Ga_max_rel.png};
+    \addplot graphics[xmin=-0.0166667, xmax=10.0167, ymin=-0.025, ymax=10.025] {../figdata/Ga_max_rel.png};
     \node[white,anchor=north west] at (axis description cs:0.05,0.95) {(b)};
   \end{axis}
   \begin{axis}[
diff --git a/tikz/harmonic_current.tex b/tikz/harmonic_current.tex
index d11b0a4..98a554a 100644
--- a/tikz/harmonic_current.tex
+++ b/tikz/harmonic_current.tex
@@ -1,5 +1,4 @@
-\definecolor{green}{rgb}{0,.67,0}%
-%\tikzsetnextfilename{harmonic_current}%
+\tikzsetnextfilename{harmonic_current}%
 \begin{tikzpicture}
   \begin{axis}[
       name = current,
diff --git a/tikz/harmonic_modes.tex b/tikz/harmonic_modes.tex
index 1be66a6..2ecb39d 100644
--- a/tikz/harmonic_modes.tex
+++ b/tikz/harmonic_modes.tex
@@ -1,5 +1,4 @@
-\definecolor{green}{rgb}{0,.67,0}%
-%\tikzsetnextfilename{harmonic_modes}%
+\tikzsetnextfilename{harmonic_modes}%
 \begin{tikzpicture}
   \begin{axis}[
       name = omega1,
diff --git a/tikz/kogan04.tex b/tikz/kogan04.tex
index 56310b3..432511b 100644
--- a/tikz/kogan04.tex
+++ b/tikz/kogan04.tex
@@ -1,37 +1,10 @@
-% Configure figures
-\definecolor{frtrg}{rgb}{.95,0,0}%
-\definecolor{adiabatic}{rgb}{0,0,0}%
-\definecolor{phenomenological}{rgb}{0,0,1}%
-\definecolor{analytic}{rgb}{0,.6,0}%
-\definecolor{experiment}{rgb}{0,0,0}%
-\definecolor{ac}{rgb}{.5,.5,.5}%
-\definecolor{green}{rgb}{0,.67,0}%
-\pgfplotsset{adiabatic style/.style={densely dotted, thick, color=adiabatic}}%
-\pgfplotsset{frtrg style/.style={solid, line width=0.75pt, color=frtrg, no marks}}%
-\pgfplotsset{phenomenological style/.style={color=phenomenological}}%
-\pgfplotsset{ac style/.style={thick, dashdotted, color=ac}}%
-\pgfplotsset{analytic style/.style={thick, densely dashed}}%
-\definecolor{c02}{rgb}{.5,.5,.5}%
-\definecolor{c04}{rgb}{.7,.25,.35}%
-\definecolor{c06}{rgb}{.95,.1,0}%
-\definecolor{c08}{rgb}{.9,.5,0}%
-\definecolor{c10}{rgb}{.67,.54,0}%
-\definecolor{c12}{rgb}{.5,.75,0}%
-\definecolor{c14}{rgb}{0,.62,0}%
-\definecolor{c16}{rgb}{0,.8,.5}%
-\definecolor{c18}{rgb}{0,.6,.7}%
-\definecolor{c20}{rgb}{.25,.3,1}%
-\definecolor{c22}{rgb}{.6,.35,.9}%
-\definecolor{c24}{rgb}{.9,0,1}%
-\definecolor{c30}{rgb}{.8,.25,.6}%
-\definecolor{lightgray}{rgb}{.6,.6,.6}%
-\newcommand\plotsep{0.02}%
-%\tikzsetnextfilename{kogan04}%
+\def\plotsep{0.02}%
+\tikzsetnextfilename{kogan04}%
 \begin{tikzpicture}
   \begin{axis}[
       %ymin = 0.037, ymax = 0.205,
       ymin = 0.037, ymax = 0.188,
-      xmin = -7.3, xmax = 6.7,
+      xmin = -7, xmax = 7,
       mark = *,
       width = 9cm,
       height = 8cm,
@@ -45,6 +18,7 @@
       clip marker paths = true,
       xtick distance = 2,
       minor x tick num = 1,
+      axis on top,
       legend entries = {
         {\hspace*{-1.8em}{$\vac~(\mu V)$}},
         \hphantom{1}29, \hphantom{1}45, \hphantom{1}60, \hphantom{1}67, 144,
@@ -70,46 +44,46 @@
     ]
     \addlegendimage{empty legend};
     \addplot[only marks, on layer=main, color=c06, forget plot] table[x expr={17.95*\thisrowno{1}}, y expr={\thisrowno{0}/2 + 4*\plotsep}] {../exp_data/d764n766_didv_trace0++.txt};
-    \addplot[only marks, on layer=main, color=c08, forget plot] table[x expr={17.95*\thisrowno{1}}, y expr={\thisrowno{0}/2 + 3*\plotsep}] {../exp_data/d764n766_didv_trace10++.txt};
+    \addplot[only marks, on layer=main, color=c10, forget plot] table[x expr={17.95*\thisrowno{1}}, y expr={\thisrowno{0}/2 + 3*\plotsep}] {../exp_data/d764n766_didv_trace10++.txt};
     \addplot[only marks, on layer=main, color=c14, forget plot] table[x expr={17.95*\thisrowno{1}}, y expr={\thisrowno{0}/2 + 2*\plotsep}] {../exp_data/d764n766_didv_trace21++.txt};
     \addplot[only marks, on layer=main, color=c20, forget plot] table[x expr={17.95*\thisrowno{1}}, y expr={\thisrowno{0}/2 + 1*\plotsep}] {../exp_data/d764n766_didv_trace26++.txt};
     \addplot[only marks, on layer=main, color=c24, forget plot] table[x expr={17.95*\thisrowno{1}}, y expr={\thisrowno{0}/2 + 0*\plotsep}] {../exp_data/d764n766_didv_trace40++.txt};
     %
-    %\addplot[frtrg style, cyan, on layer=main] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g0.754asym}+ 4*\plotsep}] {../figdata/kogan04.dat};
-    %\addplot[frtrg style, cyan, on layer=main] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.17asym} + 3*\plotsep}] {../figdata/kogan04.dat};
-    %\addplot[frtrg style, cyan, on layer=main] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.56asym} + 2*\plotsep}] {../figdata/kogan04.dat};
-    %\addplot[frtrg style, cyan, on layer=main] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.74asym} + 1*\plotsep}] {../figdata/kogan04.dat};
-    %\addplot[frtrg style, cyan, on layer=main] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g3.74asym} + 0*\plotsep}] {../figdata/kogan04.dat};
+    %\addplot[exp style, cyan, on layer=main] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g0.754asym}+ 4*\plotsep}] {../figdata/kogan04.dat};
+    %\addplot[exp style, cyan, on layer=main] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.17asym} + 3*\plotsep}] {../figdata/kogan04.dat};
+    %\addplot[exp style, cyan, on layer=main] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.56asym} + 2*\plotsep}] {../figdata/kogan04.dat};
+    %\addplot[exp style, cyan, on layer=main] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.74asym} + 1*\plotsep}] {../figdata/kogan04.dat};
+    %\addplot[exp style, cyan, on layer=main] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g3.74asym} + 0*\plotsep}] {../figdata/kogan04.dat};
     % asymmetry:
-    \addplot[frtrg style, black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g0.729}+ 4*\plotsep}] {../figdata/kogan04.dat};
-    \addplot[frtrg style, black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.13} + 3*\plotsep}] {../figdata/kogan04.dat};
-    \addplot[frtrg style, black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.51} + 2*\plotsep}] {../figdata/kogan04.dat};
-    \addplot[frtrg style, black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.68} + 1*\plotsep}] {../figdata/kogan04.dat};
-    \addplot[frtrg style, black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g3.62} + 0*\plotsep}] {../figdata/kogan04.dat};
+    \addplot[exp style, black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g0.729}+ 4*\plotsep}] {../figdata/kogan04.dat};
+    \addplot[exp style, black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.13} + 3*\plotsep}] {../figdata/kogan04.dat};
+    \addplot[exp style, black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.51} + 2*\plotsep}] {../figdata/kogan04.dat};
+    \addplot[exp style, black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.68} + 1*\plotsep}] {../figdata/kogan04.dat};
+    \addplot[exp style, black, line width=0.6pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g3.62} + 0*\plotsep}] {../figdata/kogan04.dat};
     %
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g0.729asym}+ 4*\plotsep}] {../figdata/kogan04.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.13asym} + 3*\plotsep}] {../figdata/kogan04.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.51asym} + 2*\plotsep}] {../figdata/kogan04.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.68asym} + 1*\plotsep}] {../figdata/kogan04.dat};
-    \addplot[frtrg style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g3.62asym} + 0*\plotsep}] {../figdata/kogan04.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g0.729asym}+ 4*\plotsep}] {../figdata/kogan04.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.13asym} + 3*\plotsep}] {../figdata/kogan04.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.51asym} + 2*\plotsep}] {../figdata/kogan04.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.68asym} + 1*\plotsep}] {../figdata/kogan04.dat};
+    \addplot[exp style, on layer=main, color=white, opacity=0.8, line width=1.65pt, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g3.62asym} + 0*\plotsep}] {../figdata/kogan04.dat};
     %
-    %\addplot[frtrg style, on layer=main] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g0.754}+ 4*\plotsep}] {../figdata/kogan04.dat};
-    %\addplot[frtrg style, on layer=main] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.17} + 3*\plotsep}] {../figdata/kogan04.dat};
-    %\addplot[frtrg style, on layer=main] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.56} + 2*\plotsep}] {../figdata/kogan04.dat};
-    %\addplot[frtrg style, on layer=main] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.74} + 1*\plotsep}] {../figdata/kogan04.dat};
-    %\addplot[frtrg style, on layer=main] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g3.74} + 0*\plotsep}] {../figdata/kogan04.dat};
-    \addplot[frtrg style, on layer=main, color=c06] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g0.729asym}+ 4*\plotsep}] {../figdata/kogan04.dat};
-    \addplot[frtrg style, on layer=main, color=c08] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.13asym} + 3*\plotsep}] {../figdata/kogan04.dat};
-    \addplot[frtrg style, on layer=main, color=c14] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.51asym} + 2*\plotsep}] {../figdata/kogan04.dat};
-    \addplot[frtrg style, on layer=main, color=c20] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.68asym} + 1*\plotsep}] {../figdata/kogan04.dat};
-    \addplot[frtrg style, on layer=main, color=c24] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g3.62asym} + 0*\plotsep}] {../figdata/kogan04.dat};
+    %\addplot[exp style, on layer=main] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g0.754}+ 4*\plotsep}] {../figdata/kogan04.dat};
+    %\addplot[exp style, on layer=main] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.17} + 3*\plotsep}] {../figdata/kogan04.dat};
+    %\addplot[exp style, on layer=main] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.56} + 2*\plotsep}] {../figdata/kogan04.dat};
+    %\addplot[exp style, on layer=main] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.74} + 1*\plotsep}] {../figdata/kogan04.dat};
+    %\addplot[exp style, on layer=main] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g3.74} + 0*\plotsep}] {../figdata/kogan04.dat};
+    \addplot[exp style, on layer=main, color=c06] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g0.729asym}+ 4*\plotsep}] {../figdata/kogan04.dat};
+    \addplot[exp style, on layer=main, color=c10] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.13asym} + 3*\plotsep}] {../figdata/kogan04.dat};
+    \addplot[exp style, on layer=main, color=c14] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.51asym} + 2*\plotsep}] {../figdata/kogan04.dat};
+    \addplot[exp style, on layer=main, color=c20] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.68asym} + 1*\plotsep}] {../figdata/kogan04.dat};
+    \addplot[exp style, on layer=main, color=c24] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g3.62asym} + 0*\plotsep}] {../figdata/kogan04.dat};
     %
-    \addplot[frtrg style, black, line width=0.6pt] coordinates {(8,0) (9,0)};
+    \addplot[exp style, black, line width=0.6pt] coordinates {(8,0) (9,0)};
     %
-    %\addplot[frtrg style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g0.729pat}+ 4*\plotsep}] {../figdata/kogan04.dat};
-    %\addplot[frtrg style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.13pat} + 3*\plotsep}] {../figdata/kogan04.dat};
-    %\addplot[frtrg style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.51pat} + 2*\plotsep}] {../figdata/kogan04.dat};
-    %\addplot[frtrg style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.68pat} + 1*\plotsep}] {../figdata/kogan04.dat};
-    %\addplot[frtrg style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g3.62pat} + 0*\plotsep}] {../figdata/kogan04.dat};
+    %\addplot[exp style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g0.729pat}+ 4*\plotsep}] {../figdata/kogan04.dat};
+    %\addplot[exp style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.13pat} + 3*\plotsep}] {../figdata/kogan04.dat};
+    %\addplot[exp style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.51pat} + 2*\plotsep}] {../figdata/kogan04.dat};
+    %\addplot[exp style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g1.68pat} + 1*\plotsep}] {../figdata/kogan04.dat};
+    %\addplot[exp style, color=lightgray, line width=0.4pt, on layer=main, forget plot] table[x expr={\thisrow{vdc}+0.3}, y expr={\thisrow{g3.62pat} + 0*\plotsep}] {../figdata/kogan04.dat};
   \end{axis}
 \end{tikzpicture}
diff --git a/tikz/omega5-3d.tex b/tikz/omega5-3d.tex
index 099c435..32233c3 100644
--- a/tikz/omega5-3d.tex
+++ b/tikz/omega5-3d.tex
@@ -1,4 +1,4 @@
-%\tikzsetnextfilename{gacdc_vdc_vac}%
+\tikzsetnextfilename{omega5_3d}%
 \begin{tikzpicture}
   \begin{axis}[
       name = G3d,
diff --git a/tikz/overview.tex b/tikz/overview.tex
index d526711..5d2bf2e 100644
--- a/tikz/overview.tex
+++ b/tikz/overview.tex
@@ -1,3 +1,4 @@
+\tikzsetnextfilename{overview}%
 \begin{tikzpicture}
   \begin{axis}[
       width=9.5cm,
diff --git a/tikz/preamble.tex b/tikz/preamble.tex
new file mode 100644
index 0000000..e3508d4
--- /dev/null
+++ b/tikz/preamble.tex
@@ -0,0 +1,90 @@
+\newcommand\tikzsetnextfilename[1]{}
+\newcommand\tkv{T_\mathrm{K}}
+\newcommand\tkrg{T_\mathrm{K}'}
+\newcommand\vdc{V_\mathrm{avg}}
+\newcommand\vac{V_\mathrm{osc}}
+\newcommand\gdc{G_\mathrm{avg}}
+\newcommand\gac{G_\mathrm{osc}}
+\newcommand\idc{I_\mathrm{avg}}
+\newcommand\iac{I_\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
+
+
+
+\definecolor{frtrg}{rgb}{.95,0,0}%
+\definecolor{adiabatic}{rgb}{0,0,0}%
+\definecolor{phenomenological}{rgb}{0,0,1}%
+\definecolor{analytic}{rgb}{0,.6,0}%
+\definecolor{experiment}{rgb}{0,0,0}%
+\definecolor{ac}{rgb}{.5,.5,.5}%
+\definecolor{green}{rgb}{0,.67,0}%
+\pgfplotsset{adiabatic style/.style={densely dotted, thick, color=adiabatic}}%
+\pgfplotsset{frtrg style/.style={solid, thick, color=frtrg, no marks}}%
+\pgfplotsset{phenomenological style/.style={color=phenomenological}}%
+\pgfplotsset{ac style/.style={thick, dashdotted, color=ac}}%
+\pgfplotsset{analytic style/.style={thick, densely dashed}}%
+\pgfplotsset{exp style/.style={solid, line width=0.75pt, no marks}}%
+
+\definecolor{c02}{rgb}{.5,.5,.5}%
+\definecolor{c04}{rgb}{.7,.25,.35}%
+\definecolor{c06}{rgb}{.95,.1,0}%
+\definecolor{c08}{rgb}{.9,.5,0}%
+\definecolor{c10}{rgb}{.67,.54,0}%
+\definecolor{c12}{rgb}{.5,.75,0}%
+\definecolor{c14}{rgb}{0,.62,0}%
+\definecolor{c16}{rgb}{0,.8,.5}%
+\definecolor{c18}{rgb}{0,.6,.7}%
+\definecolor{c20}{rgb}{.25,.3,1}%
+\definecolor{c22}{rgb}{.6,.35,.9}%
+\definecolor{c24}{rgb}{.9,0,1}%
+\definecolor{c30}{rgb}{.8,.25,.6}%
+\definecolor{lightgray}{rgb}{.6,.6,.6}%
+\definecolor{darkgreen}{HTML}{00c000}%
diff --git a/tikz/pulse_charge.tex b/tikz/pulse_charge.tex
index 4585eef..2147a79 100644
--- a/tikz/pulse_charge.tex
+++ b/tikz/pulse_charge.tex
@@ -1,17 +1,4 @@
-% Configure figures
-\definecolor{frtrg}{rgb}{.95,0,0}%
-\definecolor{adiabatic}{rgb}{0,0,0}%
-\definecolor{phenomenological}{rgb}{0,0,1}%
-\definecolor{analytic}{rgb}{0,.6,0}%
-\definecolor{experiment}{rgb}{0,0,0}%
-\definecolor{ac}{rgb}{.5,.5,.5}%
-\definecolor{green}{rgb}{0,.67,0}%
-\pgfplotsset{adiabatic style/.style={densely dotted, thick, color=adiabatic}}%
-\pgfplotsset{frtrg style/.style={solid, thick, color=frtrg}}%
-\pgfplotsset{phenomenological style/.style={color=phenomenological}}%
-\pgfplotsset{ac style/.style={thick, dashdotted, color=ac}}%
-\pgfplotsset{analytic style/.style={thick, densely dashed}}%
-%\tikzsetnextfilename{pulse_current}%
+\tikzsetnextfilename{pulse_charge}%
 \begin{tikzpicture}
   \begin{axis}[
       name = qp,
diff --git a/tikz/pulse_current.tex b/tikz/pulse_current.tex
index 2728f6c..536eefc 100644
--- a/tikz/pulse_current.tex
+++ b/tikz/pulse_current.tex
@@ -1,17 +1,4 @@
-% Configure figures
-\definecolor{frtrg}{rgb}{.95,0,0}%
-\definecolor{adiabatic}{rgb}{0,0,0}%
-\definecolor{phenomenological}{rgb}{0,0,1}%
-\definecolor{analytic}{rgb}{0,.6,0}%
-\definecolor{experiment}{rgb}{0,0,0}%
-\definecolor{ac}{rgb}{.5,.5,.5}%
-\definecolor{green}{rgb}{0,.67,0}%
-\pgfplotsset{adiabatic style/.style={densely dotted, thick, color=adiabatic}}%
-\pgfplotsset{frtrg style/.style={solid, thick, color=frtrg}}%
-\pgfplotsset{phenomenological style/.style={color=phenomenological}}%
-\pgfplotsset{ac style/.style={thick, dashdotted, color=ac}}%
-\pgfplotsset{analytic style/.style={thick, densely dashed}}%
-%\tikzsetnextfilename{pulse_current}%
+\tikzsetnextfilename{pulse_current}%
 \begin{tikzpicture}
   \begin{axis}[
       name = izoom,
-- 
GitLab