diff --git a/latex/images/burn_in_time_energy_1_0.pdf b/latex/images/burn_in_time_energy_1_0.pdf new file mode 100644 index 0000000..c58334f Binary files /dev/null and b/latex/images/burn_in_time_energy_1_0.pdf differ diff --git a/latex/images/burn_in_time_energy_2_4.pdf b/latex/images/burn_in_time_energy_2_4.pdf new file mode 100644 index 0000000..e23072e Binary files /dev/null and b/latex/images/burn_in_time_energy_2_4.pdf differ diff --git a/latex/images/burn_in_time_magnetization_1_0.pdf b/latex/images/burn_in_time_magnetization_1_0.pdf new file mode 100644 index 0000000..eac658c Binary files /dev/null and b/latex/images/burn_in_time_magnetization_1_0.pdf differ diff --git a/latex/images/burn_in_time_magnetization_2_4.pdf b/latex/images/burn_in_time_magnetization_2_4.pdf new file mode 100644 index 0000000..0f90aef Binary files /dev/null and b/latex/images/burn_in_time_magnetization_2_4.pdf differ diff --git a/latex/images/phase_transition/fox/narrow/10M/energy.pdf b/latex/images/phase_transition/fox/narrow/10M/energy.pdf new file mode 100644 index 0000000..54df904 Binary files /dev/null and b/latex/images/phase_transition/fox/narrow/10M/energy.pdf differ diff --git a/latex/images/phase_transition/fox/narrow/10M/heat_capacity.pdf b/latex/images/phase_transition/fox/narrow/10M/heat_capacity.pdf new file mode 100644 index 0000000..25bc859 Binary files /dev/null and b/latex/images/phase_transition/fox/narrow/10M/heat_capacity.pdf differ diff --git a/latex/images/phase_transition/fox/narrow/10M/linreg.pdf b/latex/images/phase_transition/fox/narrow/10M/linreg.pdf new file mode 100644 index 0000000..422b6af Binary files /dev/null and b/latex/images/phase_transition/fox/narrow/10M/linreg.pdf differ diff --git a/latex/images/phase_transition/fox/narrow/10M/magnetization.pdf b/latex/images/phase_transition/fox/narrow/10M/magnetization.pdf new file mode 100644 index 0000000..2d2e1e3 Binary files /dev/null and b/latex/images/phase_transition/fox/narrow/10M/magnetization.pdf differ diff --git a/latex/images/phase_transition/fox/narrow/10M/susceptibility.pdf b/latex/images/phase_transition/fox/narrow/10M/susceptibility.pdf new file mode 100644 index 0000000..1d64a4a Binary files /dev/null and b/latex/images/phase_transition/fox/narrow/10M/susceptibility.pdf differ diff --git a/latex/images/phase_transition/fox/wide/10M/energy.pdf b/latex/images/phase_transition/fox/wide/10M/energy.pdf new file mode 100644 index 0000000..b546179 Binary files /dev/null and b/latex/images/phase_transition/fox/wide/10M/energy.pdf differ diff --git a/latex/images/phase_transition/fox/wide/10M/heat_capacity.pdf b/latex/images/phase_transition/fox/wide/10M/heat_capacity.pdf new file mode 100644 index 0000000..3c0242c Binary files /dev/null and b/latex/images/phase_transition/fox/wide/10M/heat_capacity.pdf differ diff --git a/latex/images/phase_transition/fox/wide/10M/linreg.pdf b/latex/images/phase_transition/fox/wide/10M/linreg.pdf new file mode 100644 index 0000000..5b85e49 Binary files /dev/null and b/latex/images/phase_transition/fox/wide/10M/linreg.pdf differ diff --git a/latex/images/phase_transition/fox/wide/10M/magnetization.pdf b/latex/images/phase_transition/fox/wide/10M/magnetization.pdf new file mode 100644 index 0000000..5e43e16 Binary files /dev/null and b/latex/images/phase_transition/fox/wide/10M/magnetization.pdf differ diff --git a/latex/images/phase_transition/fox/wide/10M/susceptibility.pdf b/latex/images/phase_transition/fox/wide/10M/susceptibility.pdf new file mode 100644 index 0000000..9bd2a32 Binary files /dev/null and b/latex/images/phase_transition/fox/wide/10M/susceptibility.pdf differ diff --git a/latex/images/phase_transition/fox/wide/1M/energy.pdf b/latex/images/phase_transition/fox/wide/1M/energy.pdf new file mode 100644 index 0000000..42bb99c Binary files /dev/null and b/latex/images/phase_transition/fox/wide/1M/energy.pdf differ diff --git a/latex/images/phase_transition/fox/wide/1M/heat_capacity.pdf b/latex/images/phase_transition/fox/wide/1M/heat_capacity.pdf new file mode 100644 index 0000000..b86622d Binary files /dev/null and b/latex/images/phase_transition/fox/wide/1M/heat_capacity.pdf differ diff --git a/latex/images/phase_transition/fox/wide/1M/linreg.pdf b/latex/images/phase_transition/fox/wide/1M/linreg.pdf new file mode 100644 index 0000000..3876087 Binary files /dev/null and b/latex/images/phase_transition/fox/wide/1M/linreg.pdf differ diff --git a/latex/images/phase_transition/fox/wide/1M/magnetization.pdf b/latex/images/phase_transition/fox/wide/1M/magnetization.pdf new file mode 100644 index 0000000..ba4c903 Binary files /dev/null and b/latex/images/phase_transition/fox/wide/1M/magnetization.pdf differ diff --git a/latex/images/phase_transition/fox/wide/1M/susceptibility.pdf b/latex/images/phase_transition/fox/wide/1M/susceptibility.pdf new file mode 100644 index 0000000..af2b3c5 Binary files /dev/null and b/latex/images/phase_transition/fox/wide/1M/susceptibility.pdf differ diff --git a/latex/images/phase_transition/hp/energy.pdf b/latex/images/phase_transition/hp/energy.pdf new file mode 100644 index 0000000..bf64478 Binary files /dev/null and b/latex/images/phase_transition/hp/energy.pdf differ diff --git a/latex/images/phase_transition/hp/heat_capacity.pdf b/latex/images/phase_transition/hp/heat_capacity.pdf new file mode 100644 index 0000000..ad90c26 Binary files /dev/null and b/latex/images/phase_transition/hp/heat_capacity.pdf differ diff --git a/latex/images/phase_transition/hp/linreg.pdf b/latex/images/phase_transition/hp/linreg.pdf new file mode 100644 index 0000000..e8ca142 Binary files /dev/null and b/latex/images/phase_transition/hp/linreg.pdf differ diff --git a/latex/images/phase_transition/hp/magnetization.pdf b/latex/images/phase_transition/hp/magnetization.pdf new file mode 100644 index 0000000..7d1cfb4 Binary files /dev/null and b/latex/images/phase_transition/hp/magnetization.pdf differ diff --git a/latex/images/phase_transition/hp/susceptibility.pdf b/latex/images/phase_transition/hp/susceptibility.pdf new file mode 100644 index 0000000..a70bf10 Binary files /dev/null and b/latex/images/phase_transition/hp/susceptibility.pdf differ diff --git a/latex/images/test_burn_in/energy.pdf b/latex/images/test_burn_in/energy.pdf new file mode 100644 index 0000000..d3d431c Binary files /dev/null and b/latex/images/test_burn_in/energy.pdf differ diff --git a/latex/images/test_burn_in/heat_capacity.pdf b/latex/images/test_burn_in/heat_capacity.pdf new file mode 100644 index 0000000..eb8c051 Binary files /dev/null and b/latex/images/test_burn_in/heat_capacity.pdf differ diff --git a/latex/images/test_burn_in/magnetization.pdf b/latex/images/test_burn_in/magnetization.pdf new file mode 100644 index 0000000..9d5de5d Binary files /dev/null and b/latex/images/test_burn_in/magnetization.pdf differ diff --git a/latex/images/test_burn_in/susceptibility.pdf b/latex/images/test_burn_in/susceptibility.pdf new file mode 100644 index 0000000..7bd2d37 Binary files /dev/null and b/latex/images/test_burn_in/susceptibility.pdf differ diff --git a/python_scripts/burn_in_time.py b/python_scripts/burn_in_time.py index ef49d3f..429f26d 100644 --- a/python_scripts/burn_in_time.py +++ b/python_scripts/burn_in_time.py @@ -48,13 +48,13 @@ def plot_1_0(): ax1.plot(t, energy, label=f"{label}", color=color) ax2.plot(t, magnetization, label=f"{label}", color=color) - ax1.set_xlabel("t") - ax1.set_ylabel(r"$\langle \epsilon \rangle$") + ax1.set_xlabel("t (MC cycles)") + ax1.set_ylabel(r"$\langle \epsilon \rangle$ $(J)$") ax1.legend(title="Initial state", loc="upper right") - fig1.savefig("./latex/images/burn_in_time_energy_1_0.pdf", bbox_inches="tight") - ax2.set_ylabel(r"$\langle |m| \rangle$") - ax2.set_xlabel("t") + + ax2.set_ylabel(r"$\langle |m| \rangle$ (unitless)") + ax2.set_xlabel("t (MC cycles)") ax2.legend(title="Initial state", loc="upper right") fig2.savefig( "./latex/images/burn_in_time_magnetization_1_0.pdf", bbox_inches="tight" @@ -94,13 +94,13 @@ def plot_2_4(): ax1.plot(t, energy, label=f"{label}", color=color) ax2.plot(t, magnetization, label=f"{label}", color=color) - ax1.set_xlabel("t") - ax1.set_ylabel(r"$\langle \epsilon \rangle$") + ax1.set_xlabel("t (MC cycles)") + ax1.set_ylabel(r"$\langle \epsilon \rangle$ $(J)$") ax1.legend(title="Initial state", loc="upper right") - fig1.savefig("./latex/images/burn_in_time_energy_2_4.pdf", bbox_inches="tight") - ax2.set_ylabel(r"$\langle |m| \rangle$") - ax2.set_xlabel("t") + + ax2.set_ylabel(r"$\langle |m| \rangle$ (unitless)") + ax2.set_xlabel("t (MC cycles)") ax2.legend(title="Initial state", loc="upper right") fig2.savefig( "./latex/images/burn_in_time_magnetization_2_4.pdf", bbox_inches="tight" diff --git a/python_scripts/pd_estimate.py b/python_scripts/pd_estimate.py index b1508c6..400e92c 100644 --- a/python_scripts/pd_estimate.py +++ b/python_scripts/pd_estimate.py @@ -40,7 +40,7 @@ def plot(infile, outfile): ec="white", lw=0.2, ) - ax1.set_xlabel(r"$\epsilon$") + ax1.set_xlabel(r"$\epsilon$ $(J)$") ax1.set_ylabel("Density") mu, sigma = np.mean(eps), np.std(eps) diff --git a/python_scripts/phase_transition.py b/python_scripts/phase_transition.py index c8186ca..d6f5b28 100644 --- a/python_scripts/phase_transition.py +++ b/python_scripts/phase_transition.py @@ -4,6 +4,20 @@ from pathlib import Path import matplotlib.pyplot as plt import numpy as np from scipy.stats import linregress +import seaborn as sns + +sns.set_theme() +params = { + "font.family": "Serif", + "font.serif": "Roman", + "text.usetex": True, + "axes.titlesize": "large", + "axes.labelsize": "large", + "xtick.labelsize": "large", + "ytick.labelsize": "large", + "legend.fontsize": "medium", +} +plt.rcParams.update(params) def plot_phase_transition_alt(indir, outdir): @@ -65,20 +79,40 @@ def plot_phase_transition_alt(indir, outdir): x = np.linspace(0, 1 / 20, 1001) regression = linregress(inv_L, Tc) f = lambda x: regression[0] * x + regression[1] + stats = ( + f"$\\beta_{0}$ = {regression[1]:.4f}\n" + f"$\\beta_{1}$ = {regression[0]:.4f}" + ) + bbox = dict(boxstyle="round", pad=0.3, fc="white", ec="white", alpha=0.5) + ax5.text( + 0.6, 0.6, stats, bbox=bbox, transform=ax1.transAxes, ha="right", va="center" + ) ax5.scatter(inv_L, Tc) - ax5.plot(x, f(x), label=f"m = {regression[0]}, i = {regression[1]}") + ax5.plot(x, f(x)) - figure1.legend() - figure2.legend() - figure3.legend() - figure4.legend() - figure5.legend() + ax1.set_xlabel(r"T $(J/k_{B})$") + ax1.set_ylabel(r"$\langle \epsilon \rangle$ $(J)$") + ax1.legend(title="Lattice size", loc="upper right") - figure1.savefig(Path(outdir, "energy.pdf")) - figure2.savefig(Path(outdir, "magnetization.pdf")) - figure3.savefig(Path(outdir, "heat_capacity.pdf")) - figure4.savefig(Path(outdir, "susceptibility.pdf")) - figure5.savefig(Path(outdir, "linreg.pdf")) + ax2.set_xlabel(r"T $(J/k_{B})$") + ax2.set_ylabel(r"$\langle |m| \rangle$ $(unitless)$") + ax2.legend(title="Lattice size", loc="upper right") + + ax3.set_xlabel(r"T $(J/k_{B})$") + ax3.set_ylabel(r"$C_{V}$ $(k_{B})$") + ax3.legend(title="Lattice size", loc="upper right") + + ax4.set_xlabel(r"T $(J/k_{B})$") + ax4.set_ylabel(r"$\chi$ $(1/J)$") + ax4.legend(title="Lattice size", loc="upper right") + + # ax5.legend() + + figure1.savefig(Path(outdir, "energy.pdf"), bbox_inches="tight") + figure2.savefig(Path(outdir, "magnetization.pdf"), bbox_inches="tight") + figure3.savefig(Path(outdir, "heat_capacity.pdf"), bbox_inches="tight") + figure4.savefig(Path(outdir, "susceptibility.pdf"), bbox_inches="tight") + figure5.savefig(Path(outdir, "linreg.pdf"), bbox_inches="tight") plt.close(figure1) plt.close(figure2) @@ -99,11 +133,11 @@ def plot_phase_transition(indir, outdir): "size_100.txt", ] labels = [ - "L = 20", - "L = 40", - "L = 60", - "L = 80", - "L = 100", + "20", + "40", + "60", + "80", + "100", ] figure1, ax1 = plt.subplots() @@ -151,20 +185,40 @@ def plot_phase_transition(indir, outdir): x = np.linspace(0, 1 / 20, 1001) regression = linregress(inv_L, Tc) f = lambda x: regression[0] * x + regression[1] + stats = ( + f"$\\beta_{0}$ = {regression[1]:.4f}\n" + f"$\\beta_{1}$ = {regression[0]:.4f}" + ) + bbox = dict(boxstyle="round", pad=0.3, fc="white", ec="white", alpha=0.5) + ax5.text( + 0.6, 0.6, stats, bbox=bbox, transform=ax1.transAxes, ha="right", va="center" + ) ax5.scatter(inv_L, Tc) - ax5.plot(x, f(x), label=f"m = {regression[0]}, i = {regression[1]}") + ax5.plot(x, f(x)) - figure1.legend() - figure2.legend() - figure3.legend() - figure4.legend() - figure5.legend() + ax1.set_xlabel(r"T $(J/k_{B})$") + ax1.set_ylabel(r"$\langle \epsilon \rangle$ $(J)$") + ax1.legend(title="Lattice size", loc="upper right") - figure1.savefig(Path(outdir, "energy.pdf")) - figure2.savefig(Path(outdir, "magnetization.pdf")) - figure3.savefig(Path(outdir, "heat_capacity.pdf")) - figure4.savefig(Path(outdir, "susceptibility.pdf")) - figure5.savefig(Path(outdir, "linreg.pdf")) + ax2.set_xlabel(r"T $(J/k_{B})$") + ax2.set_ylabel(r"$\langle |m| \rangle$ $(unitless)$") + ax2.legend(title="Lattice size", loc="upper right") + + ax3.set_xlabel(r"T $(J/k_{B})$") + ax3.set_ylabel(r"$C_{V}$ $(k_{B})$") + ax3.legend(title="Lattice size", loc="upper right") + + ax4.set_xlabel(r"T $(J/k_{B})$") + ax4.set_ylabel(r"$\chi$ $(1/J)$") + ax4.legend(title="Lattice size", loc="upper right") + + # ax5.legend() + + figure1.savefig(Path(outdir, "energy.pdf"), bbox_inches="tight") + figure2.savefig(Path(outdir, "magnetization.pdf"), bbox_inches="tight") + figure3.savefig(Path(outdir, "heat_capacity.pdf"), bbox_inches="tight") + figure4.savefig(Path(outdir, "susceptibility.pdf"), bbox_inches="tight") + figure5.savefig(Path(outdir, "linreg.pdf"), bbox_inches="tight") plt.close(figure1) plt.close(figure2) diff --git a/python_scripts/test_burn_in.py b/python_scripts/test_burn_in.py index 6360cde..25708f5 100644 --- a/python_scripts/test_burn_in.py +++ b/python_scripts/test_burn_in.py @@ -2,6 +2,20 @@ from os import makedirs from pathlib import Path import matplotlib.pyplot as plt +import seaborn as sns + +sns.set_theme() +params = { + "font.family": "Serif", + "font.serif": "Roman", + "text.usetex": True, + "axes.titlesize": "large", + "axes.labelsize": "large", + "xtick.labelsize": "large", + "ytick.labelsize": "large", + "legend.fontsize": "medium", +} +plt.rcParams.update(params) def plot(indir, outdir): @@ -13,8 +27,12 @@ def plot(indir, outdir): "burn_in.txt", ] labels = [ - "Without burn-in time", - "With burn-in time", + "Without", + "With", + ] + colors = [ + "darkred", + "seagreen", ] figure1, ax1 = plt.subplots() @@ -22,7 +40,7 @@ def plot(indir, outdir): figure3, ax3 = plt.subplots() figure4, ax4 = plt.subplots() - for file, label in zip(files, labels): + for file, label, color in zip(files, labels, colors): t = [] e = [] m = [] @@ -39,20 +57,31 @@ def plot(indir, outdir): CV.append(float(l[3])) X.append(float(l[4])) - ax1.plot(t, e, label=label) - ax2.plot(t, m, label=label) - ax3.plot(t, CV, label=label) - ax4.plot(t, X, label=label) + ax1.plot(t, e, label=label, color=color) + ax2.plot(t, m, label=label, color=color) + ax3.plot(t, CV, label=label, color=color) + ax4.plot(t, X, label=label, color=color) - figure1.legend() - figure2.legend() - figure3.legend() - figure4.legend() + ax1.set_xlabel(r"T $(J/k_{B})$") + ax1.set_ylabel(r"$\langle \epsilon \rangle$ $(J)$") + ax1.legend(title="Burn-in time", loc="upper right") - figure1.savefig(Path(outdir, "energy.pdf")) - figure2.savefig(Path(outdir, "magnetization.pdf")) - figure3.savefig(Path(outdir, "heat_capacity.pdf")) - figure4.savefig(Path(outdir, "susceptibility.pdf")) + ax2.set_xlabel(r"T $(J/k_{B})$") + ax2.set_ylabel(r"$\langle |m| \rangle$ $(unitless)$") + ax2.legend(title="Burn-in time", loc="upper right") + + ax3.set_xlabel(r"T $(J/k_{B})$") + ax3.set_ylabel(r"$C_{V}$ $(k_{B})$") + ax3.legend(title="Burn-in time", loc="upper right") + + ax4.set_xlabel(r"T $(J/k_{B})$") + ax4.set_ylabel(r"$\chi$ $(1/J)$") + ax4.legend(title="Burn-in time", loc="upper right") + + figure1.savefig(Path(outdir, "energy.pdf"), bbox_inches="tight") + figure2.savefig(Path(outdir, "magnetization.pdf"), bbox_inches="tight") + figure3.savefig(Path(outdir, "heat_capacity.pdf"), bbox_inches="tight") + figure4.savefig(Path(outdir, "susceptibility.pdf"), bbox_inches="tight") plt.close(figure1) plt.close(figure2)