From 38743e76d8cf16a7f45e3f65b14b26eb51c480a8 Mon Sep 17 00:00:00 2001
From: Jarvist Moore Frost <jarvist@gmail.com>
Date: Fri, 20 Oct 2017 14:26:17 +0100
Subject: [PATCH] Additional figure for SI of https://arxiv.org/abs/1708.04158

V and W coplotted with Schultz's Polaron radius r_f
---
 HalidePerovskites/HalidePerovskites.gpt | 20 ++++++++++++--------
 HalidePerovskites/ThermalPathways.gpt   | 19 +++++++++++++++++++
 2 files changed, 31 insertions(+), 8 deletions(-)
 create mode 100644 HalidePerovskites/ThermalPathways.gpt

diff --git a/HalidePerovskites/HalidePerovskites.gpt b/HalidePerovskites/HalidePerovskites.gpt
index 72ec4e39..9d323b5c 100644
--- a/HalidePerovskites/HalidePerovskites.gpt
+++ b/HalidePerovskites/HalidePerovskites.gpt
@@ -106,16 +106,20 @@ p 'MAPI-electron.dat' u 1:4 w lp title "MAPI-electron",\
 load "gnuplot-render.gpt"
 
 # This one is now a bit of a mess, but leaving it here to mine the axes setup
-#prefix=FILE."-variational"
+FILE="MAPI-electron"
+prefix=FILE."-VW-vs-rf"
 
-#set y2range [0:1]
+set y2range [0:50]
+set yrange [0:25]
+set xrange [0:400] # beyond here just boring linear. Including 300K point.
+set ylabel "Variational parameters"
+set y2label "Schultz r_f (Å)"
 #set y2label "Mass / Tau (reduced units)"
-#set y2tics
-#p FILE.'.dat' u 1:15 w lp title "Polaron radius", \
+set y2tics
+p FILE.'.dat' u 1:($15*10**10) w lp axes x1y2 title "Polaron radius", \
+ '' u 1:6 w lp  title "v", '' u 1:7 w lp  title "w"
+load "gnuplot-render.gpt"
+
 # '' u 1:6 w lp axes x1y2 title "v (variational)", '' u 1:7 w lp axes x1y2 title "w (variational)", \
 # '' u 1:9 w lp axes x1y2 title "M (Effective mass)", '' u 1:14 w lp axes x1y2 title "Scattering Tau (ps)"
 
-#load "gnuplot-render.gpt"
-#unset y2tics
-#unset y2label
-
diff --git a/HalidePerovskites/ThermalPathways.gpt b/HalidePerovskites/ThermalPathways.gpt
new file mode 100644
index 00000000..57f7e1e0
--- /dev/null
+++ b/HalidePerovskites/ThermalPathways.gpt
@@ -0,0 +1,19 @@
+# Thermal Pathways - vw vs rf, for MAPI-electron 
+set terminal dumb
+
+FILE="MAPI-electron"
+prefix=FILE."-VW-vs-rf"
+
+set y2range [0:50]
+set yrange [0:25]
+set xrange [0:400] # beyond here just boring linear. Including 300K point.
+
+set xlabel "Temperature (K)"
+set ylabel "Variational parameters"
+set y2label "Schultz r_f (Å)"
+set y2tics
+
+p FILE.'.dat' u 1:($15*10**10) w lp axes x1y2 title "Polaron radius", \
+ '' u 1:6 w lp  title "v", '' u 1:7 w lp  title "w"
+load "gnuplot-render.gpt"
+