From 786ca95479f4f03756e453959aae065f761cbd5b Mon Sep 17 00:00:00 2001 From: Jaron Kent-Dobias Date: Mon, 22 Oct 2018 12:25:41 -0400 Subject: renamed the cluster figure and added hyperlink line breaks --- .gitignore | 1 + figs/fig_cluster-scaling.gplot | 159 ++++++++++++++++++++++++++++++++++++++++ figs/fig_clusters_ising2d.gplot | 159 ---------------------------------------- makefile | 2 +- monte-carlo.tex | 5 +- 5 files changed, 164 insertions(+), 162 deletions(-) create mode 100644 figs/fig_cluster-scaling.gplot delete mode 100644 figs/fig_clusters_ising2d.gplot diff --git a/.gitignore b/.gitignore index 5e54a7f..aa5e894 100644 --- a/.gitignore +++ b/.gitignore @@ -1,4 +1,5 @@ *.pdf +*.out *.aux *.log *.bbl diff --git a/figs/fig_cluster-scaling.gplot b/figs/fig_cluster-scaling.gplot new file mode 100644 index 0000000..d22be31 --- /dev/null +++ b/figs/fig_cluster-scaling.gplot @@ -0,0 +1,159 @@ + +# data and constants, sorted by column + +isingRange = "[0.0000005:5000000]" + +dataIsing2D = "data/correlation-times/ising-2d.dat" +betaIsing2D = 1.0 / 8.0 +nuIsing2D = 1.0 +deltaIsing2D = 15.0 +gammaIsing2D = 7.0 / 4.0 + +dataIsing3D = "data/correlation-times/ising-3d.dat" +betaIsing3D = 0.326419 +nuIsing3D = 0.629971 +deltaIsing3D = 4.78984 +gammaIsing3D = 1.237075 + +pottsRange="[0.000001:5000000]" + +data3Potts2D = "data/correlation-times/potts-3.dat" +beta3Potts2D = 1.0 / 9.0 +nu3Potts2D = 5.0 / 6.0 +delta3Potts2D = 14.0 +gamma3Potts2D = 13.0 / 9.0 + +data4Potts2D = "data/correlation-times/potts-4.dat" +beta4Potts2D = 1.0 / 12.0 +nu4Potts2D = 2.0 / 3.0 +delta4Potts2D = 15.0 +gamma4Potts2D = 7.0 / 6.0 + +vectorRange="[0.000001:10000000]" + +data2Vector3D = "data/correlation-times/planar-pert.dat" +beta2Vector3D = 0.3470 +nu2Vector3D = 0.6703 +delta2Vector3D = 4.79539 +gamma2Vector3D = 1.3169 + +data3Vector3D = "data/correlation-times/heisenberg-pert.dat" +beta3Vector3D = 0.3662 +nu3Vector3D = 0.7073 +delta3Vector3D = 4.79465 +gamma3Vector3D = 1.3895 + +# define colors + +cc1 = "#5e81b5" +cc2 = "#e19c24" +cc3 = "#8fb032" +cc4 = "#eb6235" +cc5 = "#8778b3" +cc6 = "#c56e1a" +cc7 = "#5d9ec7" +cc8 = "#ffbf00" + +set palette model RGB defined ( 4 cc1, 8 cc2, 16 cc3, 32 cc4, 64 cc5, 128 cc6, 256 cc7 , 512 cc8 ) +set cbrange [4:512] + +# global plotting settings + +set terminal epslatex size 6.834, 3.3 +set output "fig_cluster-scaling.tex" + +set multiplot layout 2,3 margins 0.05, 0.995, 0.15, 0.9 spacing 0.05, 0.015 +set size ratio 1 / 1.61803398875 +set key autotitle columnhead +unset colorbox +set nokey +set logscale xy + +scaledylabel = '"$\\avg{s_{\\text{\\sc 1c}}}L^{-D\\gamma/\\nu}$" offset 2,0' + +# plotting the first row +set format x "" +unset xlabel +set ylabel @scaledylabel + +# ploting Ising 2D +set xrange @isingRange +set yrange [0.8:6] + +set title '\twodee Ising' offset -6.5,-2.5 + +plot dataIsing2D using (column("H")/column("T") * exp(log(column("L")) * betaIsing2D * deltaIsing2D / nuIsing2D)):(column("S") * exp(-log(column("L")) * (gammaIsing2D / nuIsing2D))):(column("dS") * exp(-log(column("L")) * (gammaIsing2D / nuIsing2D))):1 with yerrorbars pt 0 lw 2 palette,\ + dataIsing2D using (column("H")/column("T") * exp(log(column("L")) * betaIsing2D * deltaIsing2D / nuIsing2D)):(column("S") * exp(-log(column("L")) * (gammaIsing2D / nuIsing2D))):1 with lines palette,\ + 1.4 * exp(log(x) * 2 / deltaIsing2D) with lines lc rgb "black" + +unset ylabel + +# plotting 3 Potts 2D +set xrange @pottsRange +set yrange [0.9:6] + +set title '\twodee 3-State Potts' offset -4,-2.5 + +set key above center vertical maxrows 1 + +plot data3Potts2D using (column("H")/column("T") * exp(log(column("L")) * beta3Potts2D * delta3Potts2D / nu3Potts2D)):(column("S") * exp(-log(column("L")) * (gamma3Potts2D / nu3Potts2D))):(column("dS") * exp(-log(column("L")) * (gamma3Potts2D / nu3Potts2D))):1 with yerrorbars pt 0 lw 2 palette notitle,\ + data3Potts2D using (column("H")/column("T") * exp(log(column("L")) * beta3Potts2D * delta3Potts2D / nu3Potts2D)):(column("S") * exp(-log(column("L")) * (gamma3Potts2D / nu3Potts2D))):1 with lines palette notitle,\ + 1.4 * exp(log(x) * 2 / delta3Potts2D) with lines lc rgb "black",\ + NaN with lines lc rgb cc1 title "4",\ + NaN with lines lc rgb cc2 title "8",\ + NaN with lines lc rgb cc3 title "16",\ + NaN with lines lc rgb cc4 title "32",\ + NaN with lines lc rgb cc5 title "64",\ + NaN with lines lc rgb cc6 title "128",\ + NaN with lines lc rgb cc7 title "256",\ + NaN with lines lc rgb cc8 title "512" + +unset key + +# plotting 2 Vector 3D +set xrange @vectorRange +set yrange [0.2:80] + +set title '\threedee O(2)' offset -6.5,-2.5 + +plot data2Vector3D using (column("H")/column("T") * exp(log(column("L")) * beta2Vector3D * delta2Vector3D / nu2Vector3D)):(column("S") * exp(-log(column("L")) * (gamma2Vector3D / nu2Vector3D))):(column("dS") * exp(-log(column("L")) * (gamma2Vector3D / nu2Vector3D))):1 with yerrorbars pt 0 lw 2 palette,\ + data2Vector3D using (column("H")/column("T") * exp(log(column("L")) * beta2Vector3D * delta2Vector3D / nu2Vector3D)):(column("S") * exp(-log(column("L")) * (gamma2Vector3D / nu2Vector3D))):1 with lines palette,\ + 1 * exp(log(x) * 2 / delta2Vector3D) with lines lc rgb "black" + +# plotting the second row +set format x '$10^{%T}$' +set xlabel '$hL^{\beta\delta/\nu}$' +set ylabel @scaledylabel + +# plotting Ising 3D +set xrange @isingRange +set yrange [0.5:200] + +set title '\threedee Ising' offset -6.5,-2.5 + +plot dataIsing3D using (column("H")/column("T") * exp(log(column("L")) * betaIsing3D * deltaIsing3D / nuIsing3D)):(column("S") * exp(-log(column("L")) * (gammaIsing3D / nuIsing3D))):(column("dS") * exp(-log(column("L")) * (gammaIsing3D / nuIsing3D))):1 with yerrorbars pt 0 lw 2 palette,\ + dataIsing3D using (column("H")/column("T") * exp(log(column("L")) * betaIsing3D * deltaIsing3D / nuIsing3D)):(column("S") * exp(-log(column("L")) * (gammaIsing3D / nuIsing3D))):1 with lines palette,\ + 4 * exp(log(x) * 2 / deltaIsing3D) with lines lc rgb "black" + +unset ylabel + +# plotting 4 Potts 2D +set xrange @pottsRange +set yrange [0.8:6] + +set title '\twodee 4-State Potts' offset -4,-2.5 + +plot data4Potts2D using (column("H")/column("T") * exp(log(column("L")) * beta4Potts2D * delta4Potts2D / nu4Potts2D)):(column("S") * exp(-log(column("L")) * (gamma4Potts2D / nu4Potts2D))):(column("dS") * exp(-log(column("L")) * (gamma4Potts2D / nu4Potts2D))):1 with yerrorbars pt 0 lw 2 palette,\ + data4Potts2D using (column("H")/column("T") * exp(log(column("L")) * beta4Potts2D * delta4Potts2D / nu4Potts2D)):(column("S") * exp(-log(column("L")) * (gamma4Potts2D / nu4Potts2D))):1 with lines palette,\ + 1.4 * exp(log(x) * 2 / delta4Potts2D) with lines lc rgb "black" + +# plotting 3 Vector 3D +set xrange @vectorRange +set yrange [0.3:100] + +set title '\threedee O(3)' offset -6.5,-2.5 + +plot data3Vector3D using (column("H")/column("T") * exp(log(column("L")) * beta3Vector3D * delta3Vector3D / nu3Vector3D)):(column("S") * exp(-log(column("L")) * (gamma3Vector3D / nu3Vector3D))):(column("dS") * exp(-log(column("L")) * (gamma3Vector3D / nu3Vector3D))):1 with yerrorbars pt 0 lw 2 palette,\ + data3Vector3D using (column("H")/column("T") * exp(log(column("L")) * beta3Vector3D * delta3Vector3D / nu3Vector3D)):(column("S") * exp(-log(column("L")) * (gamma3Vector3D / nu3Vector3D))):1 with lines palette,\ + 1 * exp(log(x) * 2 / delta3Vector3D) with lines lc rgb "black" + diff --git a/figs/fig_clusters_ising2d.gplot b/figs/fig_clusters_ising2d.gplot deleted file mode 100644 index 2a4c3fb..0000000 --- a/figs/fig_clusters_ising2d.gplot +++ /dev/null @@ -1,159 +0,0 @@ - -# data and constants, sorted by column - -isingRange = "[0.0000005:5000000]" - -dataIsing2D = "data/correlation-times/ising-2d.dat" -betaIsing2D = 1.0 / 8.0 -nuIsing2D = 1.0 -deltaIsing2D = 15.0 -gammaIsing2D = 7.0 / 4.0 - -dataIsing3D = "data/correlation-times/ising-3d.dat" -betaIsing3D = 0.326419 -nuIsing3D = 0.629971 -deltaIsing3D = 4.78984 -gammaIsing3D = 1.237075 - -pottsRange="[0.000001:5000000]" - -data3Potts2D = "data/correlation-times/potts-3.dat" -beta3Potts2D = 1.0 / 9.0 -nu3Potts2D = 5.0 / 6.0 -delta3Potts2D = 14.0 -gamma3Potts2D = 13.0 / 9.0 - -data4Potts2D = "data/correlation-times/potts-4.dat" -beta4Potts2D = 1.0 / 12.0 -nu4Potts2D = 2.0 / 3.0 -delta4Potts2D = 15.0 -gamma4Potts2D = 7.0 / 6.0 - -vectorRange="[0.000001:10000000]" - -data2Vector3D = "data/correlation-times/planar-pert.dat" -beta2Vector3D = 0.3470 -nu2Vector3D = 0.6703 -delta2Vector3D = 4.79539 -gamma2Vector3D = 1.3169 - -data3Vector3D = "data/correlation-times/heisenberg-pert.dat" -beta3Vector3D = 0.3662 -nu3Vector3D = 0.7073 -delta3Vector3D = 4.79465 -gamma3Vector3D = 1.3895 - -# define colors - -cc1 = "#5e81b5" -cc2 = "#e19c24" -cc3 = "#8fb032" -cc4 = "#eb6235" -cc5 = "#8778b3" -cc6 = "#c56e1a" -cc7 = "#5d9ec7" -cc8 = "#ffbf00" - -set palette model RGB defined ( 4 cc1, 8 cc2, 16 cc3, 32 cc4, 64 cc5, 128 cc6, 256 cc7 , 512 cc8 ) -set cbrange [4:512] - -# global plotting settings - -set terminal epslatex size 6.834, 3.3 -set output "fig_clusters_ising2d.tex" - -set multiplot layout 2,3 margins 0.05, 0.995, 0.15, 0.9 spacing 0.05, 0.015 -set size ratio 1 / 1.61803398875 -set key autotitle columnhead -unset colorbox -set nokey -set logscale xy - -scaledylabel = '"$\\avg{s_{\\text{\\sc 1c}}}L^{-D\\gamma/\\nu}$" offset 2,0' - -# plotting the first row -set format x "" -unset xlabel -set ylabel @scaledylabel - -# ploting Ising 2D -set xrange @isingRange -set yrange [0.8:6] - -set title '\twodee Ising' offset -6.5,-2.5 - -plot dataIsing2D using (column("H")/column("T") * exp(log(column("L")) * betaIsing2D * deltaIsing2D / nuIsing2D)):(column("S") * exp(-log(column("L")) * (gammaIsing2D / nuIsing2D))):(column("dS") * exp(-log(column("L")) * (gammaIsing2D / nuIsing2D))):1 with yerrorbars pt 0 lw 2 palette,\ - dataIsing2D using (column("H")/column("T") * exp(log(column("L")) * betaIsing2D * deltaIsing2D / nuIsing2D)):(column("S") * exp(-log(column("L")) * (gammaIsing2D / nuIsing2D))):1 with lines palette,\ - 1.4 * exp(log(x) * 2 / deltaIsing2D) with lines lc rgb "black" - -unset ylabel - -# plotting 3 Potts 2D -set xrange @pottsRange -set yrange [0.9:6] - -set title '\twodee 3-State Potts' offset -4,-2.5 - -set key above center vertical maxrows 1 - -plot data3Potts2D using (column("H")/column("T") * exp(log(column("L")) * beta3Potts2D * delta3Potts2D / nu3Potts2D)):(column("S") * exp(-log(column("L")) * (gamma3Potts2D / nu3Potts2D))):(column("dS") * exp(-log(column("L")) * (gamma3Potts2D / nu3Potts2D))):1 with yerrorbars pt 0 lw 2 palette notitle,\ - data3Potts2D using (column("H")/column("T") * exp(log(column("L")) * beta3Potts2D * delta3Potts2D / nu3Potts2D)):(column("S") * exp(-log(column("L")) * (gamma3Potts2D / nu3Potts2D))):1 with lines palette notitle,\ - 1.4 * exp(log(x) * 2 / delta3Potts2D) with lines lc rgb "black",\ - NaN with lines lc rgb cc1 title "4",\ - NaN with lines lc rgb cc2 title "8",\ - NaN with lines lc rgb cc3 title "16",\ - NaN with lines lc rgb cc4 title "32",\ - NaN with lines lc rgb cc5 title "64",\ - NaN with lines lc rgb cc6 title "128",\ - NaN with lines lc rgb cc7 title "256",\ - NaN with lines lc rgb cc8 title "512" - -unset key - -# plotting 2 Vector 3D -set xrange @vectorRange -set yrange [0.2:80] - -set title '\threedee O(2)' offset -6.5,-2.5 - -plot data2Vector3D using (column("H")/column("T") * exp(log(column("L")) * beta2Vector3D * delta2Vector3D / nu2Vector3D)):(column("S") * exp(-log(column("L")) * (gamma2Vector3D / nu2Vector3D))):(column("dS") * exp(-log(column("L")) * (gamma2Vector3D / nu2Vector3D))):1 with yerrorbars pt 0 lw 2 palette,\ - data2Vector3D using (column("H")/column("T") * exp(log(column("L")) * beta2Vector3D * delta2Vector3D / nu2Vector3D)):(column("S") * exp(-log(column("L")) * (gamma2Vector3D / nu2Vector3D))):1 with lines palette,\ - 1 * exp(log(x) * 2 / delta2Vector3D) with lines lc rgb "black" - -# plotting the second row -set format x '$10^{%T}$' -set xlabel '$hL^{\beta\delta/\nu}$' -set ylabel @scaledylabel - -# plotting Ising 3D -set xrange @isingRange -set yrange [0.5:200] - -set title '\threedee Ising' offset -6.5,-2.5 - -plot dataIsing3D using (column("H")/column("T") * exp(log(column("L")) * betaIsing3D * deltaIsing3D / nuIsing3D)):(column("S") * exp(-log(column("L")) * (gammaIsing3D / nuIsing3D))):(column("dS") * exp(-log(column("L")) * (gammaIsing3D / nuIsing3D))):1 with yerrorbars pt 0 lw 2 palette,\ - dataIsing3D using (column("H")/column("T") * exp(log(column("L")) * betaIsing3D * deltaIsing3D / nuIsing3D)):(column("S") * exp(-log(column("L")) * (gammaIsing3D / nuIsing3D))):1 with lines palette,\ - 4 * exp(log(x) * 2 / deltaIsing3D) with lines lc rgb "black" - -unset ylabel - -# plotting 4 Potts 2D -set xrange @pottsRange -set yrange [0.8:6] - -set title '\twodee 4-State Potts' offset -4,-2.5 - -plot data4Potts2D using (column("H")/column("T") * exp(log(column("L")) * beta4Potts2D * delta4Potts2D / nu4Potts2D)):(column("S") * exp(-log(column("L")) * (gamma4Potts2D / nu4Potts2D))):(column("dS") * exp(-log(column("L")) * (gamma4Potts2D / nu4Potts2D))):1 with yerrorbars pt 0 lw 2 palette,\ - data4Potts2D using (column("H")/column("T") * exp(log(column("L")) * beta4Potts2D * delta4Potts2D / nu4Potts2D)):(column("S") * exp(-log(column("L")) * (gamma4Potts2D / nu4Potts2D))):1 with lines palette,\ - 1.4 * exp(log(x) * 2 / delta4Potts2D) with lines lc rgb "black" - -# plotting 3 Vector 3D -set xrange @vectorRange -set yrange [0.3:100] - -set title '\threedee O(3)' offset -6.5,-2.5 - -plot data3Vector3D using (column("H")/column("T") * exp(log(column("L")) * beta3Vector3D * delta3Vector3D / nu3Vector3D)):(column("S") * exp(-log(column("L")) * (gamma3Vector3D / nu3Vector3D))):(column("dS") * exp(-log(column("L")) * (gamma3Vector3D / nu3Vector3D))):1 with yerrorbars pt 0 lw 2 palette,\ - data3Vector3D using (column("H")/column("T") * exp(log(column("L")) * beta3Vector3D * delta3Vector3D / nu3Vector3D)):(column("S") * exp(-log(column("L")) * (gamma3Vector3D / nu3Vector3D))):1 with lines palette,\ - 1 * exp(log(x) * 2 / delta3Vector3D) with lines lc rgb "black" - diff --git a/makefile b/makefile index 8a25388..c3d9764 100644 --- a/makefile +++ b/makefile @@ -1,6 +1,6 @@ DOC=monte-carlo -FIGS=fig_correlation_collapse-hL fig_correlation fig_metastable fig_correlation-temp fig_clusters_ising2d fig_correlation-times fig_harmonic-susceptibilities fig_generator-times +FIGS=fig_cluster-scaling fig_correlation-times fig_harmonic-susceptibilities fig_generator-times DATA=correlation correlation_temp metastable_raw metastable_func all: ${DOC}.pdf diff --git a/monte-carlo.tex b/monte-carlo.tex index a0a26fc..52a9e4a 100644 --- a/monte-carlo.tex +++ b/monte-carlo.tex @@ -4,7 +4,8 @@ \documentclass[aps,pre,reprint]{revtex4-1} -\usepackage{amsmath,amssymb,latexsym,mathtools,algorithm,listings} +\usepackage{amsmath,amssymb,latexsym,mathtools,algorithm,listings,hyperref} +\hypersetup{breaklinks=true} % uncomment to label only equations that are referenced in the text @@ -680,7 +681,7 @@ proportional to their size, or \] \begin{figure*} - \input{fig_clusters_ising2d} + \input{fig_cluster-scaling} \caption{Collapses of rescaled average Wolff cluster size $\avg s_{\text{\sc 1c}}L^{-\gamma/\nu}$ as a function of field scaling variable $hL^{\beta\delta/\nu}$ for a variety of models. Critical exponents -- cgit v1.2.3-54-g00ecf