%  Ising model abrupt transition.
%
%  Created by Jaron Kent-Dobias on Thu Apr 20 12:50:56 EDT 2017.
%  Copyright (c) 2017 Jaron Kent-Dobias. All rights reserved.
%
\documentclass[aps,prl,reprint,fleqn]{revtex4-1}

\usepackage[utf8]{inputenc}
\usepackage[T1]{fontenc}
\usepackage{amsmath,amssymb,latexsym,concmath,mathtools,xifthen,mfpic}

\mathtoolsset{showonlyrefs=true}

\def\[{\begin{equation}}
\def\]{\end{equation}}

\def\re{\mathop{\mathrm{Re}}\nolimits}
\def\im{\mathop{\mathrm{Im}}\nolimits}
\def\dd{\mathrm d}
\def\O{\mathcal O}
\def\o{\mathcal o}
\def\ei{\mathop{\mathrm{Ei}}\nolimits}
\def\b{\mathrm b}
\def\c{\mathrm c}

\newcommand\pd[3][]{
  \ifthenelse{\isempty{#1}}
    {\def\tmp{}}
    {\def\tmp{^#1}}
  \frac{\partial\tmp#2}{\partial#3\tmp}
}

\begin{document}

\title{Essential Singularity in the Ising Abrupt Transition}
\author{Jaron Kent-Dobias}
\author{James P.~Sethna}
\affiliation{Cornell University}

\date{April 20, 2017}

\begin{abstract}
  Test
\end{abstract}

\maketitle

The Ising model is the canonical example of a system with a continuous phase
transition, and the study of its singular properties marked the first success
of the renormalization group ({\sc rg}) method in statistical physics
\cite{wilson.1971.renormalization}. This status makes sense: it's a simple
model whose phase transition admits {\sc rg} methods in a straightforward way,
and has exact solutions in certain dimensions and for certain parameter
restrictions. However, in one respect the Ising critical point is not simply a
continuous transition: it ends the line of abrupt phase transitions at zero
field below the critical temperature. Though typically neglected in {\sc rg}
scaling analysis of the critical point, we demonstrate that there are
numerically measurable contributions to scaling due to the abrupt transition
line that cannot be accounted for by analytic changes of control or
thermodynamic variables.

{\sc Rg} analysis predicts that the singular part of the free energy per site
$F$ as a function of reduced temperature $t=1-\frac{T_c}T$ and field $h=H/T$ in the vicinity of the critical point takes the scaling form
$F(t,h)=g_t^{2-\alpha}\mathcal F(g_hg_t^{-\Delta})$, where
$\Delta=\beta\delta$ and $g_t$, $g_h$ are analytic functions of $t$, $h$ that
transform exactly linearly under {\sc rg}. When studying the properties of the
Ising critical point, it is nearly always assumed that $\mathcal F(X)$, the
universal scaling function, is an analytic function of $X$. However, it has
long been known that there exists an essential singularity in $\mathcal F$ at
$X=0$, though its effects have long been believed to be unobservable
\cite{fisher.1967.condensation}. With careful analysis, we have found that
assuming the presence of the essential singularity is predictive of the
scaling form of e.g. the susceptibility.

The providence of the essential singularity can be understood using the
methods of critical droplet theory for the decay of an Ising system in a
metastable state, i.e., an equilibrium Ising state for $T<T_c$, $H>0$
subjected to a small negative external field $H<0$.
It's long been known that the decay rate $\Gamma$ of metastable states in
statistical mechanics is often related to the metastable free energy $F$ by
$\Gamma\propto\im F$
\cite{langer.1967.condensation,langer.1969.metastable,gaveau.1989.analytic}.
`Metastable free energy' can be thought of as either an analytic continuation of the free energy
through the abrupt phase transition, or restriction of the partition function
trace to states in the vicinity of the local free energy minimum that
characterizes the metastable state. In any case, the free energy develops a
nonzero imaginary part in the metastable region. Heuristically, this can be
thought of as similar to what happens in quantum mechanics with a non-unitary
Hamiltonian: the imaginary part describes loss of probability in the system
that corresponds to decay. 

In critical droplet theory, the metastable state decays when a domain of the
equilibrium state forms that is sufficiently large to grow and envelop the
whole system. Assuming the free energy cost of the surface of the droplet
scales with the number of spins $N$ like $\Sigma N^\sigma$ and that of its
bulk scales like $-MHN$, the critical droplet size scales like
$N_\c\sim(MH/\Sigma)^{-1/(1-\sigma)}$ and the free energy of the critical
droplet scales like $\Delta F_\c\sim\Sigma^{1/(1-\sigma)}(HM)^{-\sigma/(1-\sigma)}$.
Assuming domains have minimal surfaces, $\sigma=1-\frac1d$ and
$\Delta F_\c\sim\Sigma^d(HM)^{-(d-1)}$. Assuming the scaling forms
$\Sigma=g_t^\mu\mathcal S(g_hg_t^{-\Delta})$ and $M=g_t^\beta\mathcal
M(g_hg_t^{-\Delta})$ and using known hyperscaling relations
\cite{widom.1981.interface}, this implies a scaling form
\begin{align}
  \Delta F_c&
  \sim\mathcal S^d(g_hg_t^{-\Delta})(g_hg_t^{-\Delta}\mathcal
    M(g_hg_t^{-\Delta}))^{-(d-1)}\\
  &\sim\mathcal G^{-(d-1)}(g_h g_t^{-\Delta})
\end{align}
Since both surface tension and magnetization are finite and nonzero for $H=0$
at $T<T_c$, $\mathcal G(X)=\O(X)$ for small $X$.
The decay rate of the metastable state will be roughly given by the Boltzmann
factor for the creation of a critical droplet, or $\Gamma\sim e^{-\beta\Delta
F_c}$, so that
\[
  \im F\sim e^{-\mathcal G(g_hg_t^{-\Delta})^{-(d-1)}}
\]
For $d>1$ this function has an essential singularity in the invariant
combination $g_hg_t^{-\Delta}$.

This form of $\im F$ for small $h$ is known. Henceforth we will assume $h$ and
$t$ are sufficiently small that $g_t\simeq t$, $g_h\simeq h$, and all
functions of both variables can be truncated at lowest order. We make the scaling ansatz that
the imaginary part of the metastable free energy has the same singular
behavior as the real part of the equilibrium free energy, and that for small
$t$, $h$, $\im F(t,h)=t^{2-\alpha}\mathcal H(ht^{-\Delta})$ for
\[
  \mathcal H(X)=A\Theta(-X)X^\zeta e^{-1/(BX)^{d-1}}
\]
where $\Theta$ is the Heaviside function and with $\zeta=-(d-3)d/2$ for $d=2,4$ and $\zeta=-7/3$ for $d=3$
\cite{houghton.1980.metastable,gunther.1980.goldstone}. Assuming that $F$ is
analytic in the upper complex-$h$ plane, the real part of $F$ in the
equilibrium state can be extracted from this imaginary metastable free energy
using the Kramers--Kronig relation
\[
  \re F(t,h)=\frac1\pi\int_{-\infty}^\infty\frac{\im F(t,h')}{h'-h}\,\dd h'
\]
In {\sc 3d} and {\sc 4d} this can be computed explicitly given our scaling
ansatz, yielding
\[
  \begin{aligned}
    \mathcal F(X)&=
    \frac{AB^{1/3}}{12\pi X^2}e^{-1/(BX)^2}
    \bigg[\Gamma(\tfrac16)E_{7/6}((BX)^{-2})\\
    &\hspace{10em}-4BX\Gamma(\tfrac23)E_{5/3}((BX)^{-2})\bigg]
  \end{aligned}
  \notag
\]
\[
  \begin{aligned}
    \mathcal F(X)&=
    \frac{A}{9\pi X^2}e^{1/(BX)^3}
    \Big[3\Gamma(0,(BX)^{-3})\\
    &\hspace{2em}-3\Gamma(\tfrac23)\Gamma(\tfrac13,(BX)^{-3})
    -\Gamma(\tfrac13)\Gamma(-\tfrac13,(BX)^{-3})\Big]
  \end{aligned}
  \notag
\]
At the level of truncation we are working at, the Kramers--Kronig relation
does not converge in {\sc 2d}. However, the higher moments can still be
extract, e.g., the susceptibility, by taking
\[
  \chi\propto\pd[2]Fh
  =\frac2\pi\int_{-\infty}^\infty\frac{\im F(t,h')}{(h'-h)^3}\,\dd h'
\]
This yields
\[
  \chi=|t|^{-\gamma}\frac{C}{2(BX)^3}\big[BX(BX-1)-e^{1/BX}\ei(-1/BX)\big]
\]
Scaling forms for the free energy can then be extracted by integration and
comparison with known exact results at zero field, yielding
\[
  \mathcal M(X)=\frac{D}{BX}(BX-1)e^{1/BX}\ei(-1/BX)-D+\mathcal M(0)
\]
with $\mathcal
M(0)=\big(2(\sqrt2-1)\big)^{1/4}\big((4+3\sqrt2)\sinh^{-1}1\big)^{1/8}$, and
\[
  \mathcal F(X)=\mathcal F(0)+EX(\mathcal M(0)+De^{1/BX}\ei(-1/BX))
\]

How predictive are these scaling forms in the proximity of the critical point
and the abrupt transition line? We used a form of the Wolff algorithm modified
to remain efficient in the presence of an external field by incorporating the
field as another spin with coupling $|H|$ to all others
\cite{dimitrovic.1991.finite}. Data was then taken for susceptibility and
magnetization for $|t|,h\leq0.1$


\bibliographystyle{apsrev4-1}
\bibliography{essential-ising}

\end{document}