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:
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.