1 files changed, 59 insertions, 34 deletions
diff --git a/fig-stiffnesses.gplot b/fig-stiffnesses.gplot
index e132869..9fa686e 100644
--- a/fig-stiffnesses.gplot
+++ b/fig-stiffnesses.gplot
@@ -1,60 +1,85 @@
 
+# colors
 cc1 = "#5e81b5"
 cc2 = "#e19c24"
-cc3 = "#8fb032"
+cc3 = "#6fa012"
 cc4 = "#eb6235"
 
-set terminal epslatex size 8.68cm, 8.2cm standalone
+Tc = 17.26
+
+a1 = 71.1212
+b1 = 0.0104105
+c1 = 0.00378087
+d1 = 6.2662
+
+C10(T) = a1 - b1 * T
+C1(T) = C10(T) / (1 + d1 / (1 + c1 * abs(T - Tc)) / C10(T))
+
+a2 = 144.345
+b2 = 0.019492
+c2 = 120.462 
+
+C20(T) = a2 - b2 * T**2 / (c2 + T)
+
+set terminal epslatex size 8.68cm, 6.5cm standalone header \
+  '\usepackage{xcolor}\definecolor{mathc3}{HTML}{6fa012}\definecolor{mathc4}{HTML}{eb6235}'
 set output "fig-stiffnesses.tex"
 
-set multiplot  layout 3, 2  margins 0.13, 0.9, 0.1, 0.998  spacing 0.01, 0.01
+set multiplot  layout 2, 2  margins 0.1, 0.88, 0.1, 0.99  spacing 0.01, 0.01
 
 set nokey
-set xrange [0:315]
-set arrow 1 from 17.26,graph 0 to 17.26,graph 1 nohead lw 4 lc rgb cc2
 
-lam(T) = 71.1212 - 0.0104105 * T
 
-set format x ""
-set xtics 0, 100, 300 offset 0,0.5
 set mxtics 2
-set ylabel '\tiny $C / \mathrm{GPa}$' offset 4.5
-set format y '\tiny $%.1f$'
+set mx2tics 2
+set mytics 2
+set my2tics 2
+set format y '\tiny $%1.f$'
 set format y2 '\tiny $%.1f$'
+set format x '\tiny $%.0f$'
+#set format x2 '\tiny $%.0f$'
+set format x2 ''
 
+set arrow 1 from Tc,graph 0 to Tc,graph 1 nohead lw 5 lc rgb cc2
+unset xtics
+set x2tics 50, 50, 250 offset 0,-0.5 mirror
+#set x2label '\tiny $T / \mathrm K$' offset 0,-1.0
+set xrange [0:300]
+
+set ylabel '\scriptsize $C_{\mathrm{B_{2\mathrm g}}} / \mathrm{GPa}$' offset 3.5
 set yrange [140:145]
 set ytics 141,1,144 offset 0.5
-set title '\tiny (e) $C_{\mathrm{B_{2\mathrm g}}}$' offset -1,-6
 plot "data/c66.dat" using 1:(100 * $2) with lines lw 3 lc rgb cc3, \
-     144.345 - 0.019492 * x**2 / (120.462 + x) dt 3 lw 4 lc rgb cc4
-
-set ylabel ''
+     C20(x) dt 3 lw 4 lc rgb cc4
 
+unset ylabel
+set y2label '\scriptsize $C_{\mathrm{B_{1\mathrm g}}} / \mathrm{GPa}$' offset -5 rotate by -90
 set yrange [65.05:65.7]
 set y2tics 62.1,0.1,65.6 offset -0.5 mirror
-set title '\tiny (c) $C_{\mathrm{B_{1\mathrm g}}}$' offset -1,-6
-plot "data/c11mc12.dat" using 1:(100 * $2) with lines lw 3 lc rgb cc3
+plot "data/c11mc12.dat" using 1:(100 * $2) with lines lw 3 lc rgb cc3, \
+     C1(x) dt 3 lw 4 lc rgb cc4
 
-set format x '\tiny $%.0f$'
-set ylabel '\tiny $C / \mathrm{GPa}$'
-set xlabel '\tiny $T / \mathrm K$' offset 0,1
-unset y2tics
 
-set yrange [65:71.5]
+unset x2label
+unset x2tics
+
+set xlabel '\scriptsize $T / \mathrm K$' offset 0,1.25
+set xtics 50, 50, 250 offset 0,0.5 mirror
+unset y2tics
 set ytics 63,1,72
-set title '\tiny (c) $C_{\mathrm{B_{1\mathrm g}}}$' offset -1,-6
+unset y2label
+set ylabel '\scriptsize $\textcolor{mathc3}{C_{\mathrm{B_{1\mathrm g}}}},\textcolor{mathc4}{C_{\mathrm{B_{1\mathrm g}}}^0} / \mathrm{GPa}$' offset 3
+
+set yrange [64.5:71.5]
 plot "data/c11mc12.dat" using 1:(100 * $2) with lines lw 3 lc rgb cc3, \
-     lam(x) dt 3 lw 4 lc rgb cc4
+     C10(x) dt 3 lw 4 lc rgb cc4
 
 set ylabel ''
-set y2label '\tiny $\log[(C^0 - C) / \mathrm{GPa}]$' offset -5 rotate by -90
-set xlabel '\tiny $\log(T / \mathrm K+D_\perp q_*^4/a)$' offset 0,1
-
-set yrange [1:1.85]
-set format x '\tiny $%0.1f$'
-set xrange [5.5:6.35]
-set y2tics 1,0.1,1.9 offset -0.7 mirror
-set xtics 5.6,0.2,6.4
-set title '\tiny (c) $C_{\mathrm{B_{1\mathrm g}}}$' offset -1,-6
-plot "data/c11mc12.dat" using (log($1 + 1 / 0.00378087)):(log(lam($1) - 100 * $2)) with lines lw 3 lc rgb cc3, \
-     7.3 - x
+set y2label '\scriptsize $[C^0(C^0/C - 1)]^{-1} \cdot \mathrm{GPa}$' offset -5.5 rotate by -90
+
+set format y2 '\tiny $%0.2f$'
+set format y '\tiny $%0.2f$'
+set yrange [0.12:0.345]
+set y2tics 0.15,0.05,0.39 offset -0.7 mirror
+plot "data/c11mc12.dat" using 1:(1 / (C10($1)*(C10($1) / (100 * $2) - 1))) with lines lw 3 lc rgb cc3, \
+     1/(C10(x) * (C10(x) / C1(x) - 1)) dt 3 lw 4 lc rgb cc4