1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
|
cc1 = "#5e81b5"
cc2 = "#e19c24"
cc3 = "#8fb032"
cc4 = "#eb6235"
set terminal epslatex size 8.68cm, 8.2cm standalone
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 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 format y2 '\tiny $%.1f$'
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 ''
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
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]
set ytics 63,1,72
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, \
lam(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
|
