Topological structures and phases. in U(1) gauge theory. Abstract. We show that topological properties of minimal Dirac sheets as well as of

Similar documents
1 Introduction Duality transformations have provided a useful tool for investigating many theories both in the continuum and on the lattice. The term

and B. Taglienti (b) (a): Dipartimento di Fisica and Infn, Universita di Cagliari (c): Dipartimento di Fisica and Infn, Universita di Roma La Sapienza

arxiv:hep-lat/ v1 30 May 1995

of the gauge eld U which one obtains after integrating through K = (8 + 2m) 1. It is easy to show that the lattice Dirac operator

SCATTERING CONFIGURATION SPACES

Estimation of Hottest Spot Temperature in Power Transformer Windings with NDOF and DOF Cooling

SVETLANA KATOK AND ILIE UGARCOVICI (Communicated by Jens Marklof)

arxiv: v1 [hep-lat] 18 Nov 2013

arxiv: v1 [cs.fl] 24 Nov 2017

The exact region of stability for MacCormack scheme. Hoon Hong. Research Institute for Symbolic Computation. Johannes Kepler University

Notes to accompany Continuatio argumenti de mensura sortis ad fortuitam successionem rerum naturaliter contingentium applicata

Technical Report: SRI-CSL July Computer Science Laboratory, SRI International, Menlo Park. also. Technical Report: TR-98-09

hep-lat/ Dec 94

Dual Monte Carlo and Cluster Algorithms. N. Kawashima and J.E. Gubernatis. Abstract. theory and not from the usual viewpoint of a particular model.

arxiv:hep-th/ v2 8 Jun 2000

[15] For details regarding renormalization on the lattice and lattice perturbation theory, see

2 discretized variales approach those of the original continuous variales. Such an assumption is valid when continuous variales are represented as oat

Robot Position from Wheel Odometry

3.14. The model of Haldane on a honeycomb lattice

hep-lat/ May 1996

arxiv: v1 [hep-ph] 30 Dec 2015

A new explanation to the cold nuclear matter effects in heavy ion. collisions

Parallel Recombinative Simulated Annealing: A Genetic Algorithm. IlliGAL Report No July 1993

arxiv:hep-ph/ v1 19 Mar 1998

On the large time behavior of solutions of fourth order parabolic equations and ε-entropy of their attractors

Recall that the action for the GGCS model is given by: S = S YM + S CS + S H (Euclidean) (.) where S YM = ; g Z S CS = ; i g Z S H = g Z d 3 x tr F d

At first numbers were used only for counting, and 1, 2, 3,... were all that was needed. These are called positive integers.

A stochastic method for solving Smoluchowski's. Institute of Computational Mathematics and Mathematical Geophysics

Superluminal Hidden Communication as the Underlying Mechanism for Quantum Correlations: Constraining Models

Simple Examples. Let s look at a few simple examples of OI analysis.

At large distances D-branes [1, 2] interact through the long-range elds of supergravity. At substringy scales, on the other hand, their interaction is

ground state degeneracy ground state energy

Monopole Condensation and Confinement in SU(2) QCD (1) Abstract

Critical value of the total debt in view of the debts. durations

Testing Predictive Ability and Power Robusti cation

Homework 7 Solutions to Selected Problems

Instantaneous Measurement of Nonlocal Variables

arxiv: v1 [hep-ph] 20 Dec 2012

A Continuous Vibration Theory for Beams with a Vertical Edge Crack

EQUILIBRIA AND STABILITY ANALYSIS OF A BRANCHED METABOLIC NETWORK WITH FEEDBACK INHIBITION. F. Grognard Y. Chitour G. Bastin

A comparison between compact and noncompact formulation of the three dimensional lattice QED

138. Last Latexed: April 25, 2017 at 9:45 Joel A. Shapiro. so 2. iψ α j x j

Large N reduction in deformed Yang-Mills theories

condensates and topology fixing action

Exact solutions to plaquette Ising models with free and periodic boundaries

Contents. Learning Outcomes 2012/2/26. Lecture 6: Area Pattern and Spatial Autocorrelation. Dr. Bo Wu

IMSc/98/10/51 hep-th/ The Hagedorn Transition, Deconnement and the AdS/CFT Correspondence S. Kalyana Rama and B. Sathiapalan Institute of Mathe

A general theory of discrete ltering. for LES in complex geometry. By Oleg V. Vasilyev AND Thomas S. Lund

Modifying Shor s algorithm to compute short discrete logarithms

Hot Electron Production in Plasmas Illuminated by Intense Lasers

Explicit constructions of quasi-monte Carlo rules for the numerical integration of high dimensional periodic functions

Technical Report No Pursuit and Evasion on a Ring: An Innite Hierarchy for Parallel Real-Time Systems Stefan D. Bruda and Selim G. Akl Depar

N f = 1. crossover. 2nd order Z(2) m, m

Effect of Uniform Horizontal Magnetic Field on Thermal Instability in A Rotating Micropolar Fluid Saturating A Porous Medium

Konrad-Zuse-Zentrum für Informationstechnik Berlin Heilbronner Str. 10, D-1000 Berlin 31

Exploring the relationship between a fluid container s geometry and when it will balance on edge

CHAPTER 5. Linear Operators, Span, Linear Independence, Basis Sets, and Dimension

Kinetic boundary conditions in the lattice Boltzmann method

arxiv: v1 [hep-lat] 1 Oct 2007

A converse Gaussian Poincare-type inequality for convex functions

arxiv:hep-lat/ v1 10 Oct 1997

THE BALANCED DECOMPOSITION NUMBER AND VERTEX CONNECTIVITY

arxiv:hep-lat/ v6 11 Dec 2003

Calculus and linear algebra for biomedical engineering Week 3: Matrices, linear systems of equations, and the Gauss algorithm

Π( ) Π( ) Π( ) A ε. b b. (d) (c) (a) (b)

ADAM BERGER, JOHN MULVEY, ED ROTHBERG, AND ROBERT VANDERBEI. Abstract. One component of every multi-stage stochastic program is a

An n log n Algorithm for Hyper-Minimizing a (Minimized) Deterministic Automaton

Testing a Fourier Accelerated Hybrid Monte Carlo Algorithm

Determinants of generalized binary band matrices

hep-lat/ Apr 1994

ON AN INFINITE-DIMENSIONAL GROUP OVER A FINITE FIELD. A. M. Vershik, S. V. Kerov

with angular brackets denoting averages primes the corresponding residuals, then eq. (2) can be separated into two coupled equations for the time evol

hep-th/ Mar 94

hep-lat/ Dec 93

hep-ph/ Sep 1998

Structuring Unreliable Radio Networks

Exact Free Vibration of Webs Moving Axially at High Speed

Minimizing a convex separable exponential function subject to linear equality constraint and bounded variables

arxiv:hep-lat/ v1 29 Sep 1997

The High Density Region of QCD from an Effective Model

Constraining the mass distribution of cluster galaxies by weak lensing

ON SYMMETRY NON-RESTORATION AT HIGH TEMPERATURE. G. Bimonte. and. G. Lozano 1. International Centre for Theoretical Physics, P.O.

Demographic noise slows down cycles of dominance in ecological models

Laplacian modes for calorons and as a filter

Luis Manuel Santana Gallego 100 Investigation and simulation of the clock skew in modern integrated circuits. Clock Skew Model

Shape stability of pasta phases: Lasagna case

1 Introduction This work follows a paper by P. Shields [1] concerned with a problem of a relation between the entropy rate of a nite-valued stationary

Travel Grouping of Evaporating Polydisperse Droplets in Oscillating Flow- Theoretical Analysis

cond-mat/ Mar 1995

only nite eigenvalues. This is an extension of earlier results from [2]. Then we concentrate on the Riccati equation appearing in H 2 and linear quadr

Constrained Leja points and the numerical solution of the constrained energy problem

arxiv:hep-ph/ v1 21 May 2001

Quantized 2d Dilaton Gravity and. abstract. We have examined a modied dilaton gravity whose action is separable into the

University of California. Berkeley, CA fzhangjun johans lygeros Abstract

Continuous-Time Gauss-Markov Processes. Alessandro Beghi y. Abstract. Continuing the work started in [11], in this paper we examine

arxiv: v1 [cond-mat.str-el] 3 Apr 2012

arxiv:hep-th/ v1 7 Aug 2001

Depth versus Breadth in Convolutional Polar Codes

arxiv: v3 [math.qa] 11 Jul 2010

Notes on Iterated Expectations Stephen Morris February 2002

Transcription:

BUHEP-94-35 Decemer 1994 Topological structures and phases in U(1) gauge theory Werner Kerler a, Claudio Rei and Andreas Weer a a Fachereich Physik, Universitat Marurg, D-35032 Marurg, Germany Department of Physics, Boston University, Boston, MA 02215, USA hep-lat/9412066 13 Dec 1994 Astract We show that topological properties of minimal Dirac sheets as well as of currents lines characterize the phases unamiguously. We otain the minimal sheets relialy y a suitale simulated-annealing procedure. 1. Introduction We investigate compact U(1) lattice gauge theory in 4 dimensions with the action [1] X X S = cos ;x ) + jm ;x j >;x(1 ;x where M ;x = ( ;x+ ;x )=4. The physical ux ;x 2 [ ; ) is given y [2] ;x = ;x + 2n ;x. We consider periodic oundary conditions (apart from a nal remark on our results for open oundary conditions). The strength of the rst order transition decreases with, the transition ultimately getting of second order [3]. We use this to set up a very ecient algorithm [3, 4] for the Monte Carlo simulations in which ecomes a dynamical variale. Contriution to LATTICE 94, International Symposium on Lattice Field Theory, Bielefeld, Germany, 1994. Research supported in part under DFG grants Ke 250/7-2 and 250/11-1 and under DOE grant DE-FG02-91ER40676

We argue that analyzing topological structures of the congurations (rather than only measuring gloal oservales) is a promising tool. We apply it to monopole currents and Dirac plaquettes. For loops the topological characterization is straightforward. For the actually occurring networks of monopole currents we have recently een ale to present a mathematically sound characterization [3]. This is not only necessary for the unamiguous identication of physical features ut also for the computer analysis of huge networks. Here the analysis of congurations is extended to Dirac sheets and the fact that appropriately specied topological structures signal the phases is conrmed in more detail. In conclusion the general principle, not relying on the particular oundary conditions, is pointed out. 2. Dual-lattice structures On P the dual lattice for the current J ;x = M ;x+ related to links the conservation law (J ;x J ;x ) = 0 holds. The denition of current lines is such that for J ;x = 0 there is no line on the link, for J ;x = 1 there is one line, and for J ;x = 2 there are two lines, in positive or negative direction, respectively. Networks of currents are connected sets of current lines. For a network N P disconnected from the rest one gets the net current ow f ~ with components f 3 = x 0 x 1 x 2 J 3 ;x for J ;x 2 N. The Dirac string content of the plaquettes on the dual lattice is descried y p ;x = 1 2 n ;x++ for which the eld equation P (p ;x p ;x ) = J ;x holds. Dirac plaquettes then are dened such that for p ;x = 0 there is none, for p ;x = 1 there is one, and for p ;x = 2 there are two of them. Dirac sheets are formed y connecting Dirac plaquettes with common edge and appropriate orientation so that J ;x = 0 except at the oundaries of the sheets. Dirac sheet structures are not gauge invariant. However, they elong to equivalence classes which cannot e deformed into each other y gauge transformations. For given oundaries there are topologically distinct possiilities for the sheets. Thus the equivalence classes of sheets carry more information than the current networks. To represent the equivalence classes of sheets y their memers with minimal area we minimize the numer of Dirac plaquettes y gauge transformations. We nd that careful annealing is necessary to get the minimal structures relialy. For the annealing we use Metropolis sweeps for P () exp( ), where = P x;> jn ;xj, with 2

random local gauge transformations keeping the angles in [ ; ). In susequent simulations we increase the parameter in appropriate steps. 3. Topological analysis The elements of the fundamental homotopy group 1 (X; ) of a space X with ase point are equivalence classes of paths starting and ending at which can e deformed continuously into each other. Its generators may e otained emedding a suciently dense network N into X and performing suitale transformations which preserve homotopy. We use the oservation that y such a method only the generators of a sugroup are otained if a given network N does not wrap around in all directions. This provides an unamiguous characterization of networks. Choosing one vertex point of N to e the ase point and considering all paths which start and end at we note that a mapping which shrinks one edge to zero length preserves the homotopy of all of these paths. Therefore, y a sequence of such mappings we can shift all other vertices to to get a ouquet of paths starting and ending at without changing the group content. Descriing a path y a vector which is the sum of oriented steps along the path, for N with K 0 vertices and K 1 edges the ouquet has K = K 1 K 0 + 1 loops represented y vectors ~s i with i = 1; : : : ; K and j-th component s ij = w ij L j where L j is the lattice size. The ouquet matrix w ij then is to e analyzed with respect to the content of generators of 1 (T 4 ; ) = Z 4. Current networks have the additional properties of path orientation and current conservation. The maps reducing the ouquet matrix w ij, therefore, in addition to homotopy have to respect current conservation. This leads to a modied Gauss elimination procedure within which adding of a row to another one requires to sutract it simultaneously from a further row. In this way one arrives at the minimal form with rows ~a 1, : : :, ~a r, ~t, ~0,: : :, ~0 where r 4. The signicance of this gets ovious switching to the pair form with rows ~a 1, ~a 1, : : :, ~a r, ~a r, f, ~ ~0, : : :, ~0 which exhiits the relation to the net current ow f ~ explicitely. It is to e noted that f ~ 6= ~0 occurs only in very rare cases [3]. The numer of independent pairs determines the numer of nontrivial directions. For the topological analysis of minimal Dirac sheets we use the network of plaquettes as a suciently dense auxiliary network. The ouquet matrix then is otained as descried efore. The reduction of the ouquet matrix is simpler here ecause usual Gauss elimination applies, which gives the minimal form with rows ~a 1, : : :, 3

1.0 0.8 P net 0.6 0.4 0.2 0.0-0.3 0.0 0.3 Figure 1: Proaility P net for hot and cold phase as function of for lattice sizes 8 4 (circles) and16 4 (crosses). ~a r, ~0,: : :, ~0. Now the numer of independent vectors corresponds to the numer of nontrivial directions. 4. Results Figure 1 shows our results for the proaility P net to nd a current network which is nontrivial in four directions as function of in the transition region. For larger L or negative, where the peaks of the energy distriution related to the phases get well separated [3, 4], P net is seen to e 1 for the hot (conning) phase as compared to 0 for the cold (Coulom) phase. We thus have the unamiguous topological characterization y the existence of a nontrivial network in the hot phase and its asence in the cold phase. Our analysis of minimal Dirac sheets leads to analogous results as that of current networks. The hot phase is characterized y the existence of a topologically nontrivial sheet and the cold phase y its asence. Because, given a nontrivial current network, the sheets could e trivial or nontrivial, this answers the question aout the location of the sheets. We construct minimal Dirac sheets y making rst the connections where only two Dirac plaquettes meet at an edge with appropriate orientation. Then we consider 4

links 3600 hot(sheet) plaq. 600 2400 400 1200 0 cold hot(network) -0.3 0.0 0.3 200 0 Figure 2: Size of largest current network (circles, links) and largest minimal Dirac sheet (crosses, plaquettes) in hot and cold phase as function of on 8 4 lattice. the places where more than two plaquettes meet at an edge in such a way that locally there is more than one possiility to make a connection. We nd that connecting all and connecting none of the plaquettes at those places makes very little dierence in the results. Thus we get the further remarkale oservation that the minimal sheets do not intersect much. In Figure 2 we present the sizes of the largest current network and of the largest minimal Dirac sheet. Oviously the sheets give a more sensitive signal for the phases. In Figure 3 we demonstrate that our proaility P net is an order parameter superior to n max =n tot, the relative size of the largest network, advocated in [5]. In order to consider our results from a more general point of view we note that a topologically nontrivial structure on a nite lattice with periodic oundary conditions corresponds to an innite structure on an innite lattice. This suggests to characterize the hot phase more generally y the existence of an innite network (or sheet) and the cold phase y its asence, where on nite lattices \innite" is to e interpreted according to the particular oundary conditions. To test this characterization we have also performed simulations with open oundary conditions. They show that outside the transition region (which now is extended due to nite size eects) there is still a clear topological signal for the phases provided that in the analysis the prescription \nontrivial in all directions" is replaced y \touching the oundaries in all directions". This supports the indicated percolation-type picture. 5

Figure 3: Order parameters P net and n max =n tot for = 0 as function of on 8 4 lattice. References [1] J.S. Barer, R.E. Shrock, and R. Schrader, Phys. Lett. B152 (1985) 221. [2] T. DeGrand and D. Toussaint, Phys. Rev. D 22 (1980) 2478. [3] W. Kerler, C. Rei, and A. Weer, Phys. Rev. D 50 (1994) No.11. [4] W. Kerler, C. Rei, and A. Weer, \Dynamical parameter algorithm in U(1) gauge theory", contriution to these proceedings. [5] M. Baig, H. Fort, and J.B. Kogut, Preprint ILL-(TH)-94-10. 6