X-Git-Url: http://git.treefish.org/~alex/phys/proceedings_lattice13.git/blobdiff_plain/fcf0e96d40c7d11e72e34263dc6b8ebf10257a8b..81a169949997087e2e55fd6d80e5edce62725844:/proceed.tex?ds=inline diff --git a/proceed.tex b/proceed.tex index 5d1181d..12e7768 100644 --- a/proceed.tex +++ b/proceed.tex @@ -6,7 +6,7 @@ \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} @@ -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 -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. -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. } @@ -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 -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. @@ -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. -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 -(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 -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} @@ -252,6 +252,7 @@ the directions $1 \leq \nu_1 < \nu_2 < \nu_3 \leq 4$.} \label{cube} \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). @@ -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. -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} @@ -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 -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. @@ -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 -$\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.