From 8a3e115d0077db1116ca04943ef747747ff628f9 Mon Sep 17 00:00:00 2001 From: Jaron Kent-Dobias Date: Wed, 20 Oct 2021 21:18:12 +0200 Subject: More writing. --- data/phi_comparison.dat | 18 +++++++------- data/yl_comparison.dat | 8 +++++++ ising_scaling.tex | 62 +++++++++++++++++++++++++++++++++++++------------ 3 files changed, 63 insertions(+), 25 deletions(-) create mode 100644 data/yl_comparison.dat 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} -- cgit v1.2.3-70-g09d2