corrected some typos. put localup and wormup in a subsection.

[phys/proceedings_lattice13.git] / proceed.tex

diff --git a/proceed.tex b/proceed.tex
index 5d1181d33d4e0be1efbe4cbf284c133bad48620f..12e77681f304f33968b88b05394b66881a4a7fb8 100644 (file)
--- a/proceed.tex
+++ b/proceed.tex
@@ -6,7 +6,7 @@
 \usepackage{dsfont}
 \usepackage{subfigure}
 
 \usepackage{dsfont}
 \usepackage{subfigure}
 

-\title{Solving the sign problem of two-flavored scalar electrodynamics at finite chemical potential}
+\title{Solving the sign problem of two flavored scalar electrodynamics at finite chemical potential}

 
 \ShortTitle{Solving the sign problem of two-falvored scalar electrodynamics at finite chemical potential}
 
 
 \ShortTitle{Solving the sign problem of two-falvored scalar electrodynamics at finite chemical potential}
 


@@ -27,11 +27,10 @@
 We explore two-flavored scalar electrodynamics on the lattice, which has a complex phase problem 
 at finite chemical potential. By rewriting the action in terms of dual variables 
 this complex phase problem can be solved exactly. The dual variables are links and plaquettes, subject to non-trivial 
 We explore two-flavored scalar electrodynamics on the lattice, which has a complex phase problem 
 at finite chemical potential. By rewriting the action in terms of dual variables 
 this complex phase problem can be solved exactly. The dual variables are links and plaquettes, subject to non-trivial 

-constraints, which have to be respected by the Monte Carlo algorithm. 
-bvFor the simulation we use a local update that always obeys the constraints and the surface worm algorithm (SWA). 
+constraints, which have to be respected by the Monte Carlo algorithm. For the simulation we use a local update that always obeys the constraints and the surface worm algorithm (SWA). 

 The SWA is a generalization of the Prokof'ev Svistunov 
 worm algorithm concept to simulate the dual representation of abelian Gauge-Higgs models on a lattice. 
 The SWA is a generalization of the Prokof'ev Svistunov 
 worm algorithm concept to simulate the dual representation of abelian Gauge-Higgs models on a lattice. 

-We also assess the performance of the SWA and compare it with a local update in the dual representation. 
+We also assess the performance of the SWA and compare it with the local update algorithm in the dual representation. 

 Finally, we determine the full phase diagram of the model.
 }
 
 Finally, we determine the full phase diagram of the model.
 }
 


@@ -64,7 +63,7 @@ After mapping the degrees of freedom of the system to its dual variables, the we
 partition sum is positive and real and usual Monte Carlo techniques can be applied.  However, 
 the dual variables, links and plaquettes for this model, are subject to non-trivial constraints.
 Therefore one has to choose a proper algorithm in order to sample the system efficiently.  In our case, we have
 partition sum is positive and real and usual Monte Carlo techniques can be applied.  However, 
 the dual variables, links and plaquettes for this model, are subject to non-trivial constraints.
 Therefore one has to choose a proper algorithm in order to sample the system efficiently.  In our case, we have

-used two different Monte Carlo algorithms:  A local update (LMA) \cite{z3} and an extension \cite{swa} of the
+used two different Monte Carlo algorithms:  A local update algorithm (LMA) \cite{z3} and an extension \cite{swa} of the

 Prokof'ev Svistunov worm algorithm \cite{worm}.   Here we present
 some technical comparison of both algorithms in addition to the physics of the model.
  
 Prokof'ev Svistunov worm algorithm \cite{worm}.   Here we present
 some technical comparison of both algorithms in addition to the physics of the model.
  


@@ -186,13 +185,14 @@ flux variables as the new degrees of freedom.
 \vspace*{-1mm}
 \noindent
 Because the dual variables are subject to non-trivial constraints, they cannot be modified randomly during the update.
 \vspace*{-1mm}
 \noindent
 Because the dual variables are subject to non-trivial constraints, they cannot be modified randomly during the update.

-An straight forward way to sample the system is to change allowed surfaces.  
+A straight forward way to sample the system is to change allowed surfaces.  

 Thus we choose the smallest possible structures in order to
 increase the acceptance rate.  This algorithm is called local update
 Thus we choose the smallest possible structures in order to
 increase the acceptance rate.  This algorithm is called local update

-(LMA) and was used in \cite{z3,swa,prl}.  Other possibility is to use an extension of the worm
+(LMA) and was used in \cite{z3,swa,prl}.  Another possibility is to use an extension of the worm

 algorithm \cite{worm}, the so called surface worm algorithm \cite{swa}.  For this model we use both algorithms and
 algorithm \cite{worm}, the so called surface worm algorithm \cite{swa}.  For this model we use both algorithms and

-asses their performance.
+assess their performance.

 
 

+\subsection{Local update algorithm}

 Let us begin by describing the LMA. It consists of the following updates:
 \begin{itemize}
 \vspace*{-1mm}
 Let us begin by describing the LMA. It consists of the following updates:
 \begin{itemize}
 \vspace*{-1mm}


@@ -252,6 +252,7 @@ the directions $1 \leq \nu_1 < \nu_2 < \nu_3 \leq 4$.} \label{cube}
 \vspace*{-2mm}
 \end{figure}
 
 \vspace*{-2mm}
 \end{figure}
 

+\subsection{Worm algorithm}

 \noindent
 Instead of the plaquette and cube updates we can use the worm algorithm.
 Here we will shortly describe the SWA (see \cite{swa} for a detailed description).
 \noindent
 Instead of the plaquette and cube updates we can use the worm algorithm.
 Here we will shortly describe the SWA (see \cite{swa} for a detailed description).


@@ -309,7 +310,7 @@ $L_V$ where the constraints are violated.} \label{segments}
 \noindent
 In this section we describe the numerical analysis.  We first show the assessment of both algorithms
 and then the physics of the model.  In both cases we use thermodynamical observables and their fluctuations.
 \noindent
 In this section we describe the numerical analysis.  We first show the assessment of both algorithms
 and then the physics of the model.  In both cases we use thermodynamical observables and their fluctuations.

-We study in particular three observables:  the first and second derivatives with respect to the inverse
+We study in particular three observables:  The first and second derivatives with respect to the inverse

 gauge coupling $\beta$, i.e., the plaquette expectation value and its susceptibility,
 
 \begin{equation}
 gauge coupling $\beta$, i.e., the plaquette expectation value and its susceptibility,
 
 \begin{equation}


@@ -336,7 +337,7 @@ n  = \frac{1}{N_s^3 N_t}\frac{\partial}{\partial \mu} \ln\ Z\quad , \quad
 \subsection{Algorithm assessment}
 \noindent
 For the comparison of both algorithms we considered the U(1) gauge-Higgs model coupled
 \subsection{Algorithm assessment}
 \noindent
 For the comparison of both algorithms we considered the U(1) gauge-Higgs model coupled

-with two (as described previously) and with only one scalar field \cite{swa}.  
+with two (as described above) and with only one scalar field \cite{swa}.  

 First we checked the correctness of the SWA comparing the results for different 
 lattices sizes and parameters.  Examples for the one flavor model are shown in \cite{swa}.
 Fig.~\ref{obs} shows two observables for the two flavor case.
 First we checked the correctness of the SWA comparing the results for different 
 lattices sizes and parameters.  Examples for the one flavor model are shown in \cite{swa}.
 Fig.~\ref{obs} shows two observables for the two flavor case.


@@ -427,7 +428,7 @@ The dual reformulation of a problem makes it possible to look at the same physic
 
 \noindent
 In Fig.~\ref{occutrans_plaq} we plot the plaquette expectation value $\langle U \rangle$ and the corresponding susceptibility $\chi_U$ as function of the chemical potential, for two different volumes $12^3\times60$ and $16^3\times60$. We see that for the larger volume the transition is shifted slightly towards lower chemical potential, but the volume dependence seems to be reasonably small. The parameters $\beta$ and $M^2$ are fixed to $\beta=0.75$ and $M^2=5.73$. Increasing the chemical potential takes us from the confining- to the Higgs-phase where we cross the phase boundary at some critical value of $\mu$, which is $\mu\simeq2.65$ for the larger and $\mu\simeq2.7$ for the smaller lattice, telling us that the Higgs phase is tilted towards the confining phase in $\mu$-direction. Below the critical value of the chemical potential both 
 
 \noindent
 In Fig.~\ref{occutrans_plaq} we plot the plaquette expectation value $\langle U \rangle$ and the corresponding susceptibility $\chi_U$ as function of the chemical potential, for two different volumes $12^3\times60$ and $16^3\times60$. We see that for the larger volume the transition is shifted slightly towards lower chemical potential, but the volume dependence seems to be reasonably small. The parameters $\beta$ and $M^2$ are fixed to $\beta=0.75$ and $M^2=5.73$. Increasing the chemical potential takes us from the confining- to the Higgs-phase where we cross the phase boundary at some critical value of $\mu$, which is $\mu\simeq2.65$ for the larger and $\mu\simeq2.7$ for the smaller lattice, telling us that the Higgs phase is tilted towards the confining phase in $\mu$-direction. Below the critical value of the chemical potential both 

-$\langle U \rangle$ and $\chi_U$ are independent of $\mu$, which is typical for a Silverblaze type transition \cite{cohen}.
+$\langle U \rangle$ and $\chi_U$ are independent of $\mu$, which is characteristic for a Silverblaze type transition \cite{cohen}.

 
 \noindent
 Then in Fig.~\ref{occutrans} we show the occupation numbers of all dual link variables $\bar{j}$, $\bar{l}$, $j$, $l$ and dual plaquette variables $p$ just below (top) and above (bottom) the critical chemical potential $\mu_c$. Here blue links/plaquettes depict positive occupation numbers, green links/plaquettes depict negative occupation numbers and links/plaquettes with $0$-occupation are spared out. It can be seen that below $\mu_c$ links and plaquettes are hardly occupied, while above $\mu_c$ they are highly occupied. In that sense the Silverblaze transition shown in Fig.~\ref{occutrans_plaq} can be understood as condensation phenomenon, which is a new perspective on the underlying physics we gained from the dual reformulation of the problem.
 
 \noindent
 Then in Fig.~\ref{occutrans} we show the occupation numbers of all dual link variables $\bar{j}$, $\bar{l}$, $j$, $l$ and dual plaquette variables $p$ just below (top) and above (bottom) the critical chemical potential $\mu_c$. Here blue links/plaquettes depict positive occupation numbers, green links/plaquettes depict negative occupation numbers and links/plaquettes with $0$-occupation are spared out. It can be seen that below $\mu_c$ links and plaquettes are hardly occupied, while above $\mu_c$ they are highly occupied. In that sense the Silverblaze transition shown in Fig.~\ref{occutrans_plaq} can be understood as condensation phenomenon, which is a new perspective on the underlying physics we gained from the dual reformulation of the problem.