More writing.

3 files changed, 63 insertions, 25 deletions

diff --git a/data/phi_comparison.dat b/data/phi_comparison.dat
index 9193688..d58e96c 100644
--- a/data/phi_comparison.dat
+++ b/data/phi_comparison.dat
@@ -1,10 +1,8 @@
-2	0.002860955318525926	0.004496459219585747	0.0025781014469987568	0.0004361990091461404
-3	0.0005720429508622171	0.0010847239134089692	0.000805365486839224	0.00006427186448818359
-4	0.00003961608489011503	0.0001278039316774393	0.00018174532718064074	0.00013408467605927413
-5	0.0000622987443403833	0.00016919055775577174	0.0002085264051783775	0.0001300161350704411
-6	0.00005016392362722222	0.00014150435296594877	0.00016732830408854038	0.00007562595035311148
-7	0.000015703311988302104	0.00005762691693961264	0.00009304388663239349	0.00007579125219901034
-8	8.173890766238756e-6	0.000021607891761421527	0.000022883111337773654	6.9772437447900015e-6
-9	2.4873158455118727e-6	7.768088409076945e-6	9.816343265717231e-6	3.88602466879287e-6
-10	3.7198724605058686e-6	0.000011866786064407275	0.00001594714070687897	8.677026311975505e-6
-11	0.00006257848162638524	0.00006629563696586294	0.0000134012475765527	0.0000386060655095978
\ No newline at end of file
+2	0.002860955318525926	0.004496459219585747	0.0025781014469987568	0.0004361990091461404	0.0015990766131468608
+3	0.0005720429508622171	0.0010847239134089692	0.000805365486839224	0.00006427186448818359	0.0003917599113948043
+4	0.00003961608489011503	0.0001278039316774393	0.00018174532718064074	0.00013408467605927413	0.00002358497641938859
+5	0.0000622987443403833	0.00016919055775577174	0.0002085264051783775	0.0001300161350704411	3.732946060521912e-6
+6	0.00005016392362722222	0.00014150435296594877	0.00016732830408854038	0.00007562595035311148	0.00005202003334785774
+7	0.000015703311988302104	0.00005762691693961264	0.00009304388663239349	0.00007579125219901034	0.000010561158237231718
+8	8.173890766238756e-6	0.000021607891761421527	0.000022883111337773654	6.9772437447900015e-6	0.000010861362342928695
+9	2.4873158455118727e-6	7.768088409076945e-6	9.816343265717231e-6	3.88602466879287e-6	5.648983756909563e-6
\ No newline at end of file
diff --git a/data/yl_comparison.dat b/data/yl_comparison.dat
new file mode 100644
index 0000000..f1ece63
--- /dev/null
+++ b/data/yl_comparison.dat
@@ -0,0 +1,8 @@
+2	0.0156377454168716	0.00005
+3	0.0026092903452304694	0.00005
+4	0.00287653656882067	0.00005
+5	0.0030257297041298703	0.00005
+6	0.0006769621913844392	0.00005
+7	0.0005281357188132441	0.00005
+8	0.0000685310321167365	0.00005
+9	0.000035091061848013805	0.00005
\ No newline at end of file
diff --git a/ising_scaling.tex b/ising_scaling.tex
index b148eec..e3461fa 100644
--- a/ising_scaling.tex
+++ b/ising_scaling.tex
@@ -381,7 +381,7 @@ simplest form of the imaginary part to be fixed later by the real part.
 \end{tikzcd}
 \]
 We require that, for $\theta\in\mathbb R$
-\begin{equation}
+\begin{equation} \label{eq:imaginary.abrupt}
   \operatorname{Im}\mathcal F(\theta+0i)=\operatorname{Im}\mathcal F_0(\theta+0i)=C_0[\Theta(\theta-\theta_0)\mathcal I(\theta)-\Theta(-\theta-\theta_0)\mathcal I(-\theta)]
 \end{equation}
 where
@@ -493,37 +493,40 @@ abrupt transition and the Yang--Lee point, the coefficients in the analytic
 part $G$ of $\mathcal F$, and the coefficients in the undetermined function
 $g$.  Other parameters are determined by known properties.
 
-For $\theta>\theta_0$,
+For $\theta>\theta_0$, the form \eqref{eq:essential.singularity} can be
+expanded around $\theta=\theta_0$ to yield
 \begin{equation}
   \begin{aligned}
     \operatorname{Im}u_f
     &\simeq A_0 u_t(\theta)^{D\nu}\xi(\theta)\exp\left\{\frac1{b\xi(\theta)}\right\} \\
-    &=A_0R^{D\nu}(\theta_0-1)^{D\nu}\xi'(\theta_0)(\theta-\theta_0)
+    &=A_0R^{D\nu}(\theta_0^2-1)^{D\nu}\xi'(\theta_0)(\theta-\theta_0)
     \exp\left\{\frac1{b\xi'(\theta_0)}\left(\frac1{\theta-\theta_0}
       -\frac{\xi''(\theta_0)}{2\xi'(\theta_0)}\right)
       \right\}\left(1+O[(\theta-\theta_0)^2]\right)
   \end{aligned}
 \end{equation}
+Comparing this with the requirement \eqref{eq:imaginary.abrupt}, we find that
 \begin{equation}
-  B=-b\xi'(\theta_0)=-b\frac{h'(\theta_0)}{(\theta_0^2-1)^{1/\beta\delta}}
+  B=-b\xi'(\theta_0)=-b\frac{g'(\theta_0)}{(\theta_0^2-1)^{1/\beta\delta}}
 \end{equation}
+and
 \begin{equation}
   \begin{aligned}
-    C_0&=At(\theta_0)^{D\nu}\xi'(\theta_0)\exp\left\{
-    -\frac{\xi''(\theta_0)}{2\tilde B\xi'(\theta_0)^2}
+    C_0&=A_0t(\theta_0^2-1)^{D\nu}\xi'(\theta_0)\exp\left\{
+    -\frac{\xi''(\theta_0)}{2b\xi'(\theta_0)^2}
   \right\} \\
        &=
-       A|t(\theta_0)|^{D\nu-\Delta}h'(\theta_0)
-       \exp\left\{-\frac1{\tilde B}\left(\frac{|t(\theta_0)|^\Delta h''(\theta_0)}{2h'(\theta_0)^2}+\frac{\Delta|t(\theta_0)|^{\Delta - 1} t'(\theta_0)}{h'(\theta_0)}
-       \right)\right\} 
+       A_0(\theta_0^2-1)^{D\nu-\Delta}g'(\theta_0)
+       \exp\left\{-\frac1b\left(\frac{(\theta_0^2-1)^\Delta g''(\theta_0)}{2g'(\theta_0)^2}-\frac{2\Delta(\theta_0^2-1)^{\Delta - 1}\theta_0}{g'(\theta_0)}
+       \right)\right\}
   \end{aligned}
 \end{equation}
-fixing $B$ and $F_c$. Since $A$ and $\tilde B$ are known exactly, these forms can be substituted.
+fixing $B$ and $C_0$.
 
 This leaves as unknown variables the positions $\theta_0$ and
 $\theta_{\mathrm{YL}}$ of the abrupt transition and Yang--Lee edge singularity,
-the amplitude $A_\mathrm{YL}$ of the latter, and the unknown functions $F$ and
-$h$. We determine these approximately by iteration in the polynomial order at
+the amplitude $C_\mathrm{YL}$ of the latter, and the unknown functions $G$ and
+$g$. We determine these approximately by iteration in the polynomial order at
 which the free energy and its derivative matches known results. We write as a
 cost function the difference between the known series coefficients of the
 scaling functions $\mathcal F_\pm$ and the series coefficients of our
@@ -671,22 +674,51 @@ accuracy of the fit results can be checked against the known values here.
 \end{table}
 
 \begin{figure}
-  \begin{gnuplot}[terminal=epslatex, terminaloptions={size 8.65cm,5.35cm}]
+  \begin{gnuplot}[terminal=epslatex]
     dat = 'data/phi_comparison.dat'
 
     set xlabel '$n$'
-    set ylabel '$|\mathcal F_0^{(n)}-|$'
+    set ylabel '$|\Delta\mathcal F_0^{(m)}(0)|$'
 
+    set format y '$10^{%T}$'
     set style data linespoints
     set logscale y
+    set key title '\raisebox{0.5em}{$m$}' bottom left
 
     plot \
       dat using 1:2 title '0', \
       dat using 1:3 title '1', \
       dat using 1:4 title '2', \
-      dat using 1:5 title '3'
+      dat using 1:5 title '3', \
+      dat using 1:6 title '4'
+  \end{gnuplot}
+  \caption{
+    The error in the $m$th derivative of the scaling function $\mathcal F_0$
+    with respect to $\eta$ evaluated at $\eta=0$, as a function of the
+    polynomial order $n$ at which the scaling function was fit.
+  }
+\end{figure}
+
+\begin{figure}
+  \begin{gnuplot}[terminal=epslatex]
+    dat = 'data/yl_comparison.dat'
+
+    set xlabel '$n$'
+    set ylabel '$|\Delta\xi_\mathrm{YL}|$'
+    set xrange [1.5:9.5]
+    set yrange [0.000005:0.05]
+
+    set format y '$10^{%T}$'
+    set style data yerrorlines
+    set logscale y
+    unset key
+
+    plot dat using 1:2:3
   \end{gnuplot}
   \caption{
+    The error in the location of the Yang--Lee edge singularity as a function
+    of the polynomial order $n$ at which the scaling function was fit. Error
+    bars denote the uncertainty in the known location of the singularity.
   }
 \end{figure}