Modern Physics Letters B c

- © World Scientific Publishing Company

## LATTICE GAUGE THEORY FOR CONDENSED MATTER PHYSICS: FERROMAGNETIC SUPERCONDUCTIVITY AS ITS EXAMPLE

## IKUO ICHINOSE

Department of Applied Physics, Nagoya Institute of Technology, Nagoya, 466-8555, Japan

## TETSUO MATSUI

Department of Physics, Kinki University, Higashi-Osaka, 577-8502, Japan

Received (Day Month Year) Revised (Day Month Year)

Recent theoretical studies of various strongly-correlated systems in condensed matter physics reveal that the lattice gauge theory(LGT) developed in high-energy physics is quite a useful tool to understand physics of these systems. Knowledges of LGT are to become a necessary item even for condensed matter physicists. In the first part of this paper, we present a concise review of LGT for the reader who wants to understand its basics for the first time. For illustration, we choose the abelian Higgs model, a typical and quite useful LGT, which is the lattice verison of the Ginzburg-Landau model interacting with a U(1) gauge field (vector potential). In the second part, we present an account of the recent progress in the study of ferromagnetic superconductivity (SC) as an example of application of LGT to topics in condensed matter physics, . As the ferromagnetism (FM) and SC are competing orders with each other, large fluctuations are expected to take place and therefore nonperturbative methods are required for theoretical investigation. After we introduce a LGT describing the FMSC, we study its phase diagram and topological excitations (vortices of Cooper pairs) by Monte-Carlo simulations.

Keywords : Lattice gauge theory; Strongly-correlated systems; Ginzburg-Landau theory; Ferromagnetic superconductivity; Monte-Carlo simulation; Path-integral

## 1. Introduction

A strongly-correlated system is a set of either fermions or bosons which interact strongly each other. Such a system may exhibit collective behaviors that cannot be expected from the behaviors of each single particle, and has been the subject of great interest in various fields of condensed matter physics. For example, the t -J model, which is a typical model of electrons interacting via strong Coulomb repulsion, is taken as a standard model of high-temperature superconductors 1 . Another example is the Bose Hubbard model 3 and the related bosonic t -J model 4 , which are to be used to study cold bosonic atoms put on an optical lattice 5 and tuned to have

2

large interaction parameters. In the last two decades, theoretical approach using gauge theory and the associated concepts has proved itself a convincing way to understand essentials of their physics 6 7 . ,

The gauge theory put on a spatial lattice (and even on the imaginary discrete time in the path-integral formalism) is known as lattice gauge theory (LGT), which is introduced by Wilson in 1974 as a model of quark confinement 8 . Since then, it has contributed to increase our knowledges not only of quantum chromodynamics (QCD) itself, but also of nonperturbative aspects of general quantum field theory 9 . The reason to consider a lattice instead of the usual continuum space is to reduce the degrees of freedom from uncountable infinity to countable infinity, which allows us a concrete definition of the field theory in the nonperturbative region and also ways to examine it such as the strong-coupling (high-temperature) expansion. Concerning to the global phase structure, LGT shows, in addition to the ordinary Coulomb phase and Higgs phase of gauge dynamics (they are sometimes called deconfinement phases), the existence of so called confinement phase. This phase was crucial to explain confinement of quarks 8 , one of the longstanding problem in high-energy physics, because, in the confinement phase, the potential energy V ( r ) between a pair of static quark and anti-quark separated by the distance r increases linearly in r and costs infinite energy to separate them ( r →∞ ). These three phases are classified by the magnitude of fluctuations of gauge field (See Table I below). The gauge field put on the lattice can be treated as periodic variables (one calls this case a compact gauge field; see Sect.2.1), and the gauge-field configurations reflecting this periodicity may involve large spatial variations beyond the conventional perturbative ones. Monopoles are a typical and important example of such configurations (monopoles in the continuum gauge theory cost infinite energy and so not possible). In the confinement phase, such monopoles condense 10 , and the resulting strong fluctuations of gauge field are argued to generate squeezed one-dimensional line of electric flux (variables conjugate the the gauge field itself). Then, in the confinement phase, such an electric flux gives rise to the linearly-rising confining potential. This phenomenon of formation of electric fluxes is sometimes called the 'dual' Meissner effect in contrast to the Meissner effect in ordinary superconductors where magnetic field is squeezed into magnetic fluxes due to Cooper-pair condensation. One way to recognize the relevance of monopoles among other possible (topological) excitations is to perform the duality transformation 11 , which rewrites the system as an ensemble of monopoles or monopole loops.

Table 1. Three phases of gauge dynamics: V ( r ) is the potential energy stored between a pair of point charges with the opposite signs and separated by the distance r .

| hase        | fluctuation of gauge field   | potential energy V ( r )   |
|-------------|------------------------------|----------------------------|
| Higgs       | small                        | ∝ exp( - mr ) r            |
| Coulomb     | medium                       | ∝ 1 r                      |
| Confinement | large                        | ∝ r                        |

Also, in some restricted fields of condensed matter physics, LGT has been applied successfully. One example is the t -J model of high-T c superconductivity (SC) 6 . Here the strong correlations are implemented by excluding the double occupancy state of electrons, the state of two electrons with the opposite spin directions occupying the same site, from the physical states. This may be achieved by viewing an electron as a composite of fictitious particles, a so-called holon carrying the charge degree of freedom of the electron and a spinon carrying the spin degree of freedom 12 . By the method of auxiliary field in path-integral, one may introduce U(1) gauge field which binds holon and spinon. Then the possible phases may be classified according to the dynamics of this U(1) gauge field. In the confinement phase, the relevant quasiparticles are the original electrons in which holons and spinons are confined. This corresponds to the state at the overdoped high temperature region. In the deconfinement Coulomb phase, quasiparticles are holons and spinons in which charge and spin degrees of freedom behave independently. This phenomenon is called the charge-spin separation 13 and are to describe the anomalous behavior of the normal state such as the liner-T resistivity. Further correspondence between the gauge dynamics and the observed phases are also satisfactory as discussed in Refs. 1 6 13 . In particular, possibility of experimental observation of , , spinons is discussed in Ref. 2 .

Another application of LGT to relate the dynamics of auxiliary gauge field and the observed phases is the (fractional) Hall effect 14 . Among a couple of parallel formulations 15 , one may view an electron as a composite of a so called bosonic fluxon carrying fictitious (Chern-Simons) gauge flux and a so called fermionic chargon carrying the electric charge 7 . In the deconfinement phase, these fluxons and chargons are separated, i.e., particle-flux separation takes place. If fluxons Bose condense at sufficiently low temperatures, the resulting system is fermionic chargons in a reduced uniform magnetic field. These chargons are nothing but the Jain's composite fermions, and the integer Hall effect of them explain the fractional Hall effect of the original electron system.

In this paper, we shall see yet another example of LGT approach to condensed matter physics. Namely, we consider a lattice model of ferromagnetic SC, which exhibits the ferromagnetism (FM) and/or SC, two typical and important phenomena in condensed matter physics. We note here that the introduction of a spatial lattice in the original LGT in high-energy theory is traced back, as mentioned, to the necessity to define a nonperturbative field theory. The lattice spacing there should be taken as a running cutoff(scale) parameter in the renormalization-group theory 18 and the continuum limit should be taken carefully. However, in condensed matter physics, introduction of a lattice has often a practical support independent of the above meaning, i.e., the system itself has a lattice structure. One may naturally regard the lattice of the model as the real lattice of material in question. In this paper, we shall not touch the continuum limit in the sense of renormalization group.

The layout of the paper is as follows: In Sect. 2, we review LGT. The reader who is not familiar with this subject may acquire necessary knowledge to understand

## 4 I.Ichinose and T.Matsui

Fig. 1. i -j plane of the 3d lattice and ˆ φ r defined on the site r and ˆ θ ri and ˆ E ri defined on the link ( ri ) = ( r, r + ). i

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0089v2-figure-c7b4f648-ede0-42f9-a104-49c219d6ec45.jpg)

the successive sections. In Sect. 3, we explain the lattice model of ferromagnetic superconductivity (FMSC). It is based on the known Ginzburg-Landau theory defined in the continuum space. In Sect. 4, we present numerical results of Monte Carlo simulations for the lattice model in Sect. 3. In Sect. 5, we give conclusion.

## 2. Introduction to lattice gauge thoery

Formulating field-theoretical models by using space(-time) lattice is quite useful for studying their properties nonperturbatively. In particular, because of the lattice formulation itself, the high-temperature expansion (strong-coupling expansion) can be performed analytically even for the case in which fluctuations of the fields are very strong. Furthermore, numerical studies like the Monte-Carlo simulation can be applicable for wide range of the lattice field theories. In this section, we review LGT consisting of a U(1) gauge field and a charged bosonic matter field.

## 2.1. Lattice gauge theory in Hamiltonian formalism

In this subsection, we explain the Hamiltonian formalism of LGT 16 . The reasons are two fold: (i) it gives rise to an intuitive picture of the relevant states in terms of electric field and magnetic field, and (ii) it plays an important role in recent development in quantum simulations of the dynamical gauge field by using ultracold atomic systems put on an optical lattice.

Let us start with an example of Hamiltonian system defined on a threedimensional ( d = 3) spatial cubic lattice with lattice spacing a . Its sites are labeled by r , and its links (bonds) are labeled as ( r, i ) where i is the direction index i = 1 2 3. We sometimes use , , i also for the unit vector in the i th-direction, such as ( r, i ) = ( r, r + ) (See Fig.1). i

One may put matter fields on each site r , which may be bosons/fermions composing materials in question, e.g., electrons, atoms in an optical lattice, Coope-pair field of SC, spins, or even more complicated ones. For definiteness, we put below a one-component canonical nonrelativistic boson field 17 , whose annihilation and

creation operators ˆ φ , φ r ˆ † r satisfy

$$[ \hat { \phi } _ { r }, \hat { \phi } _ { r ^ { \prime } } ^ { \dagger } ] = \delta _ { r r ^ { \prime } }.$$

Gauge fields are put on each link ( r, i ) in LGT (See Fig.1). They may be associated with the gauge group U(1), SU(N), etc, and meditate interactions between matters. They are viewed as connections in differential geometry that relate two local frames of the internal space of matter field at neighboring lattice sites r and r + . i For definiteness, we consider below a U(1) gauge field operator ˆ θ ri (= ˆ θ † ri ). This is viewed as the lattice version of well-known vector potential field ˆ ( A r i ) (we use the same letter r for the coordinate of the continuum space) describing the electromagnetic(EM) interaction, but may have different origins, e.g., mediating interactions between holons and spinons, as explained in Section 1. The momentum operator ˆ E ri conjugate to ˆ θ ri describes the electric field and satisfy

$$[ \hat { \theta } _ { r i }, \hat { E } _ { r ^ { \prime } j } ] = i \delta _ { r r ^ { \prime } } \delta _ { i j }, \ [ \hat { \theta } _ { r i }, \hat { \theta } _ { r ^ { \prime } j } ] = [ \hat { E } _ { r i }, \hat { E } _ { r ^ { \prime } j } ] = 0.$$

To understand the LGT Hamiltonian, it is helpful to start with a model Hamiltonian ˆ H c defined in the continuum space. Let us consider the following GinzburgLandau(GL)-type model:

$$H _ { \text{c} } & = \int \stackrel { \ell } { d ^ { 3 } r } \Big [ \sum _ { i } | D _ { i } \hat { \phi } ( r ) | ^ { 2 } + m ^ { 2 } \hat { \rho } ( r ) + \lambda \hat { \rho } ( r ) ^ { 2 } + \frac { 1 } { 2 } \sum _ { i } \hat { E } _ { i } ( r ) ^ { 2 } + \frac { 1 } { 4 } \sum _ { i, j } \hat { F } _ { i j } ( r ) ^ { 2 } \Big ], \\ D _ { i } & = \partial _ { i } + i g \hat { A } _ { i } ( r ), \ \hat { \rho } ( r ) = \hat { \phi } ^ { \dagger } ( r ) \hat { \phi } ( r ), \ \hat { F } _ { i j } ( r ) = \partial _ { i } \hat { A } _ { j } ( r ) - \partial _ { j } \hat { A } _ { i } ( r ),$$

where D i is the covariant derivative expressing the minimal coupling of φ -particle with the EM field, g is the coupling constant (charge of φ -particle 20 ), m, λ are GL parameters whose meaning is explained later, and ˆ F ij ( r ) is the magnetic field B r i ( ) = /epsilon1 ijk F jk ( r ). One may check that the each term of H c is invariant under the local gauge transformation,

$$\hat { \phi } ( r ) \to e ^ { i \Lambda ( r ) } \hat { \phi } ( r ), \ \hat { A } _ { i } ( r ) \to \hat { A } _ { i } ( r ) - \frac { 1 } { g } \partial _ { i } \Lambda ( r ), \ \hat { E } _ { i } ( r ) \to \hat { E } _ { i } ( r ), \quad ( 2. 4 )$$

for an arbitrary function Λ( ). r By using the canonical commutation relations [ ˆ ( φ r , φ ) ˆ ( † r ′ )] = δ (3) ( r -r ′ ) and [ ˆ ( A r ,E i ) ˆ ( j r ′ )] = iδ ij δ (3) ( r -r ′ ) , one may check that its generator ˆ ( Q r ) is given by

$$\hat { Q } ( r ) = \sum _ { i } \partial _ { i } \hat { E } _ { i } ( r ) - g \hat { \rho } ( r ),$$

and ˆ ( Q r ) commutes with ˆ H c , [ ˆ ( Q r , H ) ˆ c ] = 0. So one may define the physical states | phys 〉 so that they are eigenstates of ˆ ( Q r ) (superselection rule). In the usual case, i.e., when there are no external electric field, one imposes the local constraint,

$$\hat { Q } ( r ) | \text{phys} \rangle = 0,$$

which leads to the lattice version of the Gauss' law ˆ Q r = 0 [See Eq.(2.17].

Now we consider the Hamiltonian on the cubic lattice. In LGT, the operator ˆ θ ri defined on link ( r, i ) is assumed to approach as ˆ θ ri → agA ˆ ( i r ) as the lattice

## 6 I.Ichinose and T.Matsui

spacing a becomes small. Furthermore, as we shall see, the natural U(1) gauge field operator in LGT is not ˆ θ ri but its exponential,

$$\hat { U } _ { r i } \equiv \exp ( i \hat { \theta } _ { r i } ), \, \hat { U } _ { r i } ^ { \dagger } \hat { U } _ { r i } = \hat { 1 },$$

and so ˆ θ ri is viewed as an angle operator (we discuss this point in detail later). Then, by using Eq.(2.2) we have ˆ E ri → a E 2 i ( r ) as a → 0 (note ˆ E ri is dimensionless and δ (3) ( r ) /similarequal δ r 0 a -3 ) and

$$[ \hat { E } _ { r i }, \hat { U } _ { r ^ { \prime } j } ] = \delta _ { i j } \delta _ { r r ^ { \prime } } \hat { U } _ { r i }.$$

(Note that Eq.(2.2) implies the replacement ˆ E ri → -i∂/∂θ ri .) To construct a set of eigenstates satisfying the completeness conditions, let us start with a pair of operators ˆ U and ˆ sitting on certain link (we suppress their link index). Then we E have

$$\overset { \sim } { \hat { U } } | \theta \rangle = e ^ { i \theta } | \theta \rangle, \ \theta \in [ 0, 2 \pi \}, \ \hat { E } | n \rangle = n | n \rangle, \ n \in \mathbf Z, \\ \langle \theta | n \rangle = \frac { e ^ { i n \theta } } { \sqrt { 2 \pi } }, \ \hat { 1 } _ { U } = \int _ { 0 } ^ { 2 \pi } d \theta | \theta \rangle \langle \theta |, \ \hat { 1 } _ { E } = \sum _ { n = - \infty } ^ { \infty } | n \rangle \langle n |. \quad \ ( 2. 9 ) \\ \hat { \theta } \text{ and } \hat { E } \text{ correspond to the position and momentum operators in the ordinary}$$

So ˆ and θ ˆ E correspond to the position and momentum operators in the ordinary quantum mechanics respectively, but the significant difference is that the theory is taken to be compact , i.e., the physical state | phys 〉 is imposed to be periodic in θ ,

$$\langle \theta | \text{phys} \rangle = \langle \theta + 2 \pi | \text{phys} \rangle,$$

which implies that the eigenvalues are an angle θ ∈ [0 , 2 π ) defined by mod 2 π , and integers n , instead of two real numbers. Because one may write

$$& | n \rangle = \hat { 1 } _ { U } | n \rangle = \int \frac { d \theta } { \sqrt { 2 \pi } } e ^ { i n \theta } | \theta \rangle = e ^ { i n \hat { \theta } } \cdot \int \frac { d \theta } { \sqrt { 2 \pi } } | \theta \rangle = e ^ { i n \hat { \theta } } | n = 0 \rangle = \hat { U } ^ { n } | n = 0 \rangle, ( 2. 1 1 ) \\ & \hat { U } ^ { ( \dagger ) } \text{ is just the creation (annihilation)} \text{ operator of the electric field of unit electric}$$

ˆ U ( † ) is just the creation (annihilation) operator of the electric field of unit electric flux (the signature of n distinguishes the direction of the flux), which is called a 'string bit'. These relations remind us the analogy with the angular momentum operators /vector L as ˆ E ↔ ˆ L z ( n ↔ L / z /planckover2pi1 ) , U ˆ ↔ ˆ L x + iL , θ ˆ y ↔ ϕ (azimuthal angle) with the limit of large angular momentum /lscript → ∞ . We note that there holds the uncertainty principle,

$$\Delta \theta \cdot \Delta n \geq \frac { 1 } { 2 }.$$

The Hamiltonian ˆ H of LGT should (i) reduce to ˆ H c of Eq.(2.3) in the naive continuum limit a → 0, and (ii) respect the U(1) local gauge symmetry as in the continuum. A simple example of ˆ H satisfying these two conditions is given by

$$\text{continuation}. A \, \text{on} p r e \, \text{campaign} \, \text{or} \, \text{in} \, \text{saving} \, \text{use} \, \text{when} \, \text{in} \, \text{given} \, \text{by} \\ \hat { H } & = t \sum _ { r, i } ( \hat { \phi } _ { r + i } ^ { \dagger } - \hat { U } _ { r i } \hat { \phi } _ { r } ^ { \dagger } ) ( \hat { \phi } _ { r + i } - \hat { U } _ { r i } ^ { \dagger } \hat { \phi } _ { r } ) + \sum _ { r } V ( \hat { \phi } _ { r } ) + \hat { H } _ { g }, \\ \hat { H } _ { g } & = \frac { g ^ { 2 } } { 2 a } \sum _ { r, i } ( \hat { E } _ { r i } ) ^ { 2 } + \frac { 1 } { 2 g ^ { 2 } a } \sum _ { r } \sum _ { i < j } \left ( \hat { U } _ { r j } \hat { U } _ { r + j, i } \hat { U } _ { r + i, j } ^ { \dagger } \hat { U } _ { r i } ^ { \dagger } + \text{H.c.} - 2 \right ).$$

Here V ( ˆ φ r ) represents the self-interaction of ˆ . The naive continuum limit of φ r ˆ H is checked as ˆ H → ˆ H c + O a ( 2 ) with a suitable choice of V ( ˆ φ r ) by using the following relations as a → 0:

$$\iota \iota \iota \iota \iota \sigma \, \alpha \to \, \iota. \quad & \hat { \phi } _ { r } \to a ^ { \frac { 3 } { 2 } } \hat { \phi } ( r ), \ \hat { U } _ { r i } \to e ^ { i g a \hat { A } _ { i } ( r ) } \simeq 1 + i g a \hat { A } _ { i } ( r ) - \frac { a ^ { 2 } g ^ { 2 } } { 2 } ( \hat { A } _ { i } ( r ) ) ^ { 2 }, \\ \hat { E } _ { r i } \to \frac { a ^ { 2 } } { g } \hat { E } _ { i } ( r ), \ \sum _ { r } \to \frac { 1 } { a ^ { 3 } } \int d ^ { 3 } r. \\ \iota \iota \iota \iota \iota \iota \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \, \alpha \to \}$$

In particular, the last UUUU -term reduces to the magnetic term FF as

$$& \ U U U + c. c. \to \exp [ i g a ( A + A - A - A ) ] + c. c. \\ & = \exp [ i g a ^ { 2 } ( \partial _ { i } A _ { j } ( r ) - \partial _ { j } A _ { i } ( r ) ) ] + c. c. \to 2 \cos ( g a ^ { 2 } F _ { i j } ( r ) ) \to 2 - g ^ { 2 } a ^ { 4 } F _ { i j } ^ { 2 } ( r ). \ ( 2. 1 5 )$$

For gauge invariance, one may define the gauge transformation on the lattice as

$$\hat { U } _ { r i } \rightarrow V _ { r + i } ^ { * } \hat { U } _ { r i } V _ { r }, \ V _ { r } = e ^ { i \Lambda _ { r } } \in \mathbf U ( 1 ), \ \hat { \phi } _ { r } \rightarrow V _ { r } \hat { \phi } _ { r },$$

which reduces to Eq.(2.4) by scaling Λ r → Λ( ) as r a → 0, and check that each term of ˆ H is invariant under Eq.(2.16). We define the physical states in the same manner as in the continuum case,

$$Q _ { r } \equiv \sum _ { i } \nabla _ { i } \hat { E } _ { r i } - \hat { \phi } _ { r } ^ { \dagger } \hat { \phi } _ { r }, \ [ \hat { H }, \hat { Q } _ { r } ] = 0, \ \hat { Q } _ { r } | \text{phys} \rangle = 0,$$

where the forward difference operator ∇ i is defined as ∇ i f r ≡ f r + i -f r . Both ∇ i ˆ E ri and ˆ φ φ † r ˆ r have integer eigenvalues. We note that ˆ H of Eq.(2.13) has the periodicity ˆ θ ri → ˆ θ ri +2 , which may be regarded as a gauge symmetry under a special gauge π transformation Λ r = 2 πr i .

## 2.2. Physical properties of LGT

Let us discuss the expected properties of the model (2.13). First, we focus on the case of pure gauge theory ˆ H g without matter field. In te continuum theory, the corresponding system,

$$\hat { H } _ { g } \rightarrow \frac { 1 } { 2 } \int d ^ { 3 } r \sum _ { i } \left ( \hat { E } _ { i } ( r ) \hat { E } _ { i } ( r ) + \hat { B } _ { i } ( r ) \hat { B } _ { i } ( r ) \right ),$$

is well known to describe the ensemble of free photons. In contrast, ˆ H g on the lattice is an interacting theory. One can confirm it from cos( ga F 2 ij ( r )) in Eq.(2.15). It contains, in addition to the leading F ij ( r ) 2 term, F ij ( r ) 4 and higher interaction terms describing scattering of gauge bosons. These terms are traced back to the compactness of the model (the periodicity (2.10)).

For large g /greatermuch 1, which is called the strong coupling region, one may take the electric term as the unperturbed Hamiltonian and the magnetic term as a perturbation. By recalling Eq.(2.17), the unperturbed eigenstate | eigen 〉 is the eigenstate of ˆ E ri with their eigenvalues satisfying the divergenceless condition. A simple example

## 8 I.Ichinose and T.Matsui

is an electric flux of strength n along a closed loop C on the 3D lattice. The general eigenstate may be composed of such loops. So one may write

$$| \text{eigen} \rangle = \prod _ { k } \prod _ { \ell _ { k } \in C _ { k } } ( \hat { U } _ { \ell _ { k } } ) ^ { n _ { k } } | 0 \rangle,$$

where C k denotes the k -th flux loop, /lscript k denotes links composing C k , and n k is the strength of the k -th flux. If there are two external sources of charge g at r 1 and -g at r 2 , the lowest-energy state is an electrix flux state forming a straight line ˜ C 12 (on the lattice) connecting these two charges and written as

$$\prod _ { \ell \in \tilde { C } _ { 1 2 } } \hat { U } _ { \ell } | 0 \rangle,$$

(See Fig.2a). This state has an energy g 2 / (2 a ) × number of links in C 12 , thus proportional to the distance r . If these two charges are quark and antiquark, isolation of each quark from the other implies r → ∞ and costs infinite amount of energy, being impossible to be realized. This is Wilson's explanation of quark confinement 8 .

On the other hand, for the weak-coupling region, g /lessmuch 1, the magnetic term is leading. Up to the gauge transformation, the energy eigenstates in the limit g → 0 should be the eigenstate of ˆ θ ri in contrast to the strong-coupling limit. By the uncertainty principle (2.12), the eigenstate there should be a superposition of electric fluxes of various strength and locations as in the ordinary classical Coulombic electromagnetism shown in Table.1. So the physical state corresponding to the state considered in Eq.(2.20) is replaced by

$$\sum _ { C _ { 1 2 } } \Gamma ( C _ { 1 2 } ) \prod _ { \ell \in C _ { 1 2 } } \hat { U } _ { \ell } | 0 \rangle,$$

where the coefficient Γ( C 12 ) is the weight for the path C 12 connecting the sites 1 and 2. For the lowest-energy state, Γ( C 12 ) is to be determined by the minimum-energy condition, and the energy stored in the two charges here should be proportional to 1 /r as in the Coulomb potential energy (See Fig.2b).

The above qualitative consideration may indicate that there exists a phase boundary separating the strong and weak-coupling regions. Dynamics of the pure

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0089v2-figure-d7bcb552-e008-4a96-a9e9-2e4456be83ea.jpg)

a

b

c

/negationslash

Fig. 2. Illustration of electric flux lines between two oppositely charged particles in the three phases. (a) confinement phase; (b) Coulomb phase; (c) Higgs phase. Dashed lines in (c) imply that the strength of fluxes are not conserved as Eq.(2.17) with 〈 ˆ φ φ † r ˆ r 〉 = const = 0 shows.

gauge system H g in (2.13) is quite nontrivial and careful investigation such as numerical study is required to obtain the phase diagram. At present, it is known that the 3D pure U(1) gauge system H g has two phases; The phase at strong coupling region is called the confinement phase, and the phase at weak coupling region is called Coulomb phase. The phase transition is said to be weak first order 26 . These two phases are characterized by the behavior of the potential energy V ( r ) for a pair of static charges of opposite signs and separated by a distance r ; ∝ r or ∝ 1 /r . Also they are distinguished by the fluctuation of gauge field ˆ θ ri , ∆ , as ∆ θ θ /greatermuch 1 in the confinement phase and ∆ θ /lessmuch 1 in Coulomb phase. We note that, for the 2D spatial lattice in contrast, there are no phase transitions, and only the confinement phase exists 21 .

Now let us include the matter part and consider the full Hamiltonian ˆ H of Eq.(2.13). By referring to Eq.(2.3), one may write V ( ˆ φ r ) as

$$V ( \hat { \phi } _ { r } ) = \lambda ( \hat { \phi } _ { r } ^ { \dagger } \hat { \phi } _ { r } - v ) ^ { 2 }.$$

/negationslash

Let us recall that the GL theory with the gauge interaction being turned off by setting g = 0 and for sufficiently large and positive λ , the system exhibits a secondorder phase transition by varying v . It is signaled by the order parameter φ = 〈 ˆ φ r 〉 , (or more strictly defined as | φ | 2 ≡ lim | r -r ′ |→∞ 〈 ˆ φ φ † r ˆ r ′ 〉 ), and separates the ordered (condensed) phase 〈 φ 〉 = 0 and the disordered phase 〈 φ 〉 = 0. In the mean-field theory, the transition temperature T c is determined as the vanishing point of v as v ∝ ( T c -T ). However, the fluctuations of the phase degree of freedom ϕ r of ˆ φ r = | ˆ φ r | exp( iϕ r ), which are controlled by the hopping t -term of Eq.(2.13), play an essentially important role for realization of a condensation. In fact, these two phases may be characterized by the magnitude of fluctuation ∆ ϕ of ϕ r as ∆ ϕ /lessmuch 1 (ordered phase) and ∆ ϕ /greatermuch 1(disordered phase). To obtain the correct transition temperature T c , more detailed study such as MC simulation is required as we discuss later.

For the full fledged gauge-matter system ˆ H of Eq.(2.13), one may naively consider four combinations such as (∆ θ /lessmuch 1 or /greatermuch 1) and (∆ ϕ /lessmuch 1 or /greatermuch 1) for a possible phase. However, one may easily recognize that the combination ∆ ϕ /lessmuch 1 and ∆ θ /greatermuch 1 is impossible. This is because ˆ ϕ r appears in the hopping term with the combination ∇ i ˆ ϕ r -ˆ θ ri and the condition ∆ θ /greatermuch 1 destroys the phase coherence ∇ i ˆ ϕ r /similarequal 0. In fact, the numerical study and the mean field theory predict the three phases listed in Table 1; (1) confinement phase ∆ θ /greatermuch 1 ∆ , ϕ /greatermuch 1, (2) Coulomb phase ∆ θ /lessmuch 1 ∆ , ϕ /greatermuch 1, (3) Higgs phase ∆ θ /lessmuch 1 ∆ , ϕ /lessmuch 1.

Let us see this point in some details because it helps us to understand the property of each phase appeared in an effective LGT of various condensed-matter systems.

- (1) The confinement phase appears for the parameter region of small t [ t is the hopping amplitude in Eq.(2.13)] and large g regardless of the values of λ and v . In a sense, this phase corresponds to the high-temperature ( T ) region and

## 10 I.Ichinose and T.Matsui

there exist no long-range orders. The parameter region of the confinement phase is enlarged by decreasing a value of v . For pure U(1) gauge theory without matter field ˆ , a potential φ r V ( r ) for a pair of oppositely charged static sources separated with distance r behaves as V ( r ) ∝ r as explained.

- (2) The Coulomb phase appears for small g and small t . Fluctuations of the gauge field is suppressed by the plaquette terms in Eq.(2.13) and the compactness of ˆ θ ri plays no role. The potential of a pair of two static charges is given by the Coulomb law, which is V ( r ) ∝ 1 /r for 3D space.
- (3) In the Higgs phase, a coherent condensation of the boson field ˆ φ r takes place. Then the hopping term in Eq.(2.13) gives the following term of U xµ ,

$$t \hat { \phi } _ { r + i } ^ { \dagger } \hat { U } _ { r i } ^ { \dagger } \hat { \phi } _ { r } + \text{H.c.} \to t v \ ( \hat { U } _ { r i } ^ { \dagger } + U _ { r i } ) + \dots \simeq - t v a ^ { 2 } g ^ { 2 } \hat { \theta } _ { r i } ^ { 2 } + \dots \,. \quad ( 2. 2 3 )$$

As a result, the gauge field acquires a mass m θ = √ 2 tv and a force mediated by the gauge field ˆ θ ri becomes short range like 〈 ˆ θ r + Ri,j ˆ θ tj 〉 = e -m R θ (See Fig.2c). The Higgs phase is realized for large value of t and a finite v (in addition to small g ). Superconducting (SC) phase is a kind of the Higgs phase in which the boson field ˆ φ r describes Cooper pairs, and the gauge field becomes short-ranged by the Meissner effect.

So far, we studied the compact LGT. In some literatures, so called noncompact LGT is considered 22 . The noncompact U(1) LGT is given by the following Hamiltonian for the gauge part,

$$\hat { H } _ { g } ^ { \prime } = \frac { g ^ { 2 } } { 2 a } \sum _ { r, i } ( \hat { E } _ { r i } ) ^ { 2 } + \frac { 1 } { 2 g ^ { 2 } a } \sum _ { r } \sum _ { i < j } ( \nabla _ { i } \hat { \theta } _ { r j } - \nabla _ { j } \hat { \theta } _ { r i } ) ^ { 2 }.$$

The eigenvalue of ˆ θ ri now runs in ( -∞ ∞ , ) and ˆ H ′ g loses the periodicity and describes a set of free photons. For small fluctuations ∆ θ /lessmuch 1, the compact pure gauge theory (2.13) reduces to above ˆ H ′ g . However, for large fluctuations, these two models behave quite differently and they generally have different phase diagrams. In the noncompact gauge theories with the gauge part (2.24), the Coulomb and Higgs phases are realized but the confinement phase is impossible because of the suppression of large fluctuations. In other words, the large fluctuations of ˆ θ ri in the confinement phase are achieved by large-field configurations called topological excitations, the typical examples of which are instantons and monopoles. They are allowed only for a system with periodicity such as compact LGT. To describe a SC phase transition with the ordinary EM interactions, one should use the noncompact U(1) gauge theory. The gauge field ˆ θ ri there is nothing but the vector potential for the electro-magnetism, and the single-component boson field ˆ φ r corresponds to the s -wave Cooper-pair field as mentioned. Therefore, as we have explained above, the Higgs phase is nothing but the SC phase. In order to describe a multi-component SC state such as the p -wave SC, introduction of a multi-component boson field is necessary. This system will be discussed rather in detail in the subsequent sections.

## 2.3. Partition function in the path-integral representation on the Euclidean lattice

Let us consider the partition function Z for ˆ H of Eq.(2.13) at the temperature T ,

$$Z = \text{Tr} \hat { P } \exp ( - \beta \hat { H } ), \ \beta \equiv \frac { 1 } { k _ { B } T }, \ \hat { P } \equiv \prod _ { r } \delta _ { \hat { Q } _ { r }, 0 },$$

where Tr is over the space of all the values of ˆ Q r , and ˆ is the projection operator P to the physical states. As usual, we start by factorizing the Boltzmann factor into N factors but with care of ˆ as P

$$\hat { P } \exp ( - \beta \hat { H } ) = [ \hat { P } \exp ( - \Delta \beta \hat { H } ) ] ^ { N }, \ \Delta \beta \equiv \frac { \beta } { N },$$

where we used ˆ P 2 = ˆ [ ˆ P, P,H ˆ ] = 0, and then insert the complete set between every successive factors.

For ˆ , we use the following coherent states φ r | φ 〉 :

$$\hat { \phi } | \phi \rangle = \phi | \phi \rangle, | \phi \rangle = \exp ( - \bar { \phi } \hat { \phi } + \hat { \phi } ^ { \dagger } \phi ) | 0 \rangle, \ \hat { \phi } | 0 \rangle = 0, \ \langle \phi ^ { \prime } | \phi \rangle = \exp ( - \frac { \bar { \phi } ^ { \prime } \phi ^ { \prime } } { 2 } - \frac { \bar { \phi } \phi } { 2 } + \bar { \phi } ^ { \prime } \phi ^ { \prime } ), \\ | \{ \phi \} \rangle = \prod _ { r } | \phi _ { r } \rangle, \ \int d ^ { 2 } \phi = \prod _ { r } \int \frac { d ^ { 2 } \phi _ { r } } { \pi }, \ \hat { 1 } _ { \phi } = \int d ^ { 2 } \phi | \{ \phi \} \rangle \langle \{ \phi \} |,$$

where φ is a complex number and ¯ is its complex conjugate (The bar denotes the φ complex-conjugate quantity).

The eigenstates and the completeness of gauge field for the entire lattice are written by using Eqs.(2.9) as

$$\hat { 1 } _ { \theta } = \int d \theta | \{ \theta \} \rangle \langle \{ \theta \} |, \ \hat { 1 } _ { n } = \sum _ { n } | \{ n \} \rangle \langle \{ n \} |, \ | \{ \theta \} \rangle = \prod _ { r, i } | \theta _ { r i } \rangle, \ | \{ n \} \rangle = \prod _ { r, i } | n _ { r i } \rangle, \\ \int d \theta = \prod _ { r, i } \int d \theta _ { r i }, \ \sum _ { n } = \prod _ { r, i } \sum _ { n _ { r i } \in \mathbf z }.$$

(We abbreviate the symbol ⊗ of tensor product in |{ θ }〉 and |{ n }〉 .) For example, the matrix element of ˆ is calculated as P

$$& \delta _ { \hat { Q } _ { r }, 0 } = \int _ { 0 } ^ { 2 \pi } \frac { d \theta _ { r 0 } } { 2 \pi } \exp ( i \theta _ { r 0 } \hat { Q } _ { r } ), \ \langle \{ \phi ^ { \prime }, n ^ { \prime } \} | \hat { P } | \{ \phi, n \} \rangle = \prod _ { r } \int \frac { d \theta _ { r 0 } } { 2 \pi } \prod _ { i } \delta _ { n _ { r i } ^ { \prime } n _ { r i } } \,. \\ & \prod _ { r } \exp ( - \frac { 1 } { 2 } \bar { \phi } _ { r } ^ { \prime } \phi _ { r } ^ { \prime } - \frac { 1 } { 2 } \bar { \phi } _ { r } \phi _ { r } + \bar { \phi } _ { r } ^ { \prime } e ^ { - i \theta _ { r 0 } } \phi _ { r } ) \exp ( i \theta _ { r 0 } \sum _ { i } \nabla _ { i } E _ { r i } ).$$

The Lagrange multiplier θ r 0 will be interpreted as the imaginary-time component of gauge field. Then we obtain the following expression of Z for sufficiently large

N :

$$& N \colon \\ & Z = \int [ d U ] \int [ d \phi ] \sum _ { \{ n \} } \exp ( A _ { H } ), \ [ d U ] = \prod _ { x, \mu } d U _ { x \mu }, \ [ d \phi ] = \prod _ { x } d \phi _ { x }, \sum _ { \{ n \} } = \prod _ { x, i } \sum _ { n _ { x i } }, \\ & A _ { H } = \Delta \beta \sum _ { x } \left [ \frac { i } { \Delta \beta } \sum _ { i } ( E _ { x i } \nabla _ { 0 } \theta _ { x i } + \theta _ { x 0 } \nabla _ { i } E _ { x i } ) - \frac { 1 } { \Delta \beta } \bar { \phi } _ { x + 0 } ( \phi _ { x + 0 } - e ^ { - i \theta _ { x 0 } } \phi _ { x } ) - V ( \phi _ { x } ) \\ & - t \sum _ { i } | \phi _ { x + i } - \bar { U } _ { x i } \phi _ { x } | ^ { 2 } - \frac { g ^ { 2 } } { 2 a } \sum _ { i } E _ { x i } ^ { 2 } + \frac { 1 } { 2 g ^ { 2 } a } \sum _ { i < j } ( U _ { x j } U _ { x + j, i } \bar { U } _ { x + i, j } \bar { U } _ { x i } + \text{c.c.} ) \right ].$$

$$- t \sum _ { i } | \phi _ { x + i } - \bar { U } _ { x i } \phi _ { x } | ^ { 2 } - \frac { g ^ { 2 } } { 2 a } \sum _ { i } E _ { x i } ^ { 2 } + \frac { 1 } { 2 g ^ { 2 } a } \sum _ { i < j } ( U _ { x j } U _ { x + j, i } \bar { U } _ { x + i, j } \bar { U } _ { x i } + c. c. ) \Big ].$$

where we assume that V ( ˆ φ x ) is normal ordered. In Eq.(2.30) we introduced the imaginary-time coordinate x 0 (= 0 , · · · , N -1), the site index of 4D Euclidean lattice x ≡ ( x , r 0 ) with r = ( x , x 1 2 , x 3 ), and its direction index/unit vector µ = 0 1 2 3. , , , In path integral, every field φ , θ x xµ , n xi is defined on each time slice x 0 together with the spatial site r . We note that φ , θ x xµ are periodic under x 0 → x 0 + N . The 4-component U(1) gauge field U xµ = exp( iθ xµ ) and its complex conjugate ¯ U xµ are sitting on the link ( x, x + µ ).

Next, we carry out the summation over n xi to obtain 23

$$c$$

$$\text{Next, we carry out the summation over $n_{xi}$ to obtain $\vdots$} \\ Z & = \int [ d U ] [ d \phi ] \ \exp ( A ), \ A = \sum _ { x } \left [ - \bar { \phi } _ { x + 0 } ( \phi _ { x + 0 } - \bar { U } _ { x 0 } \phi _ { x } ) - \Delta \beta V ( \phi _ { x } ) \\ & \quad - t \Delta \beta \sum _ { i } | \phi _ { x + i } - \bar { U } _ { x i } \phi _ { x } | ^ { 2 } + \frac { 1 } { 2 } \sum _ { \mu < \nu } c _ { 2 \mu \nu } ( U _ { x \nu } U _ { x + \nu, \mu } \bar { U } _ { x + \mu, \nu } \bar { U } _ { x \mu } + c. c. ) \right ], \\ c _ { 2 i j } & = \frac { \Delta \beta } { g ^ { 2 } a }, \ c _ { 2 0 i } = \frac { a } { g ^ { 2 } \Delta \beta }, \ U _ { x 0 } \equiv \exp ( i \theta _ { x 0 } ). \\ \pi \dots \dots \dots \lambda _ { x } - e. \pi \dots \dots \dots \dots \lambda _ { x } - \dots \dots \dots \dots \lambda _ { x } - \dots \dots \dots \dots \lambda _ { x }.$$

$$- \, \overline { g ^ { 2 } a }, \, \lambda _ { 2 0 n } - \overline { g ^ { 2 } \Delta \beta }, \, \mathfrak { c }$$

Here we made the following replacement in ∑ n xi respecting the periodicity in θ :

$$\sum _ { n } e ^ { - a n ^ { 2 } + i n \theta } \simeq \int _ { - \infty } ^ { \infty } d n \ e ^ { - a n ^ { 2 } + i n \theta } = C \exp ( - \frac { \theta ^ { 2 } } { 4 a } ) \to C ^ { \prime } \exp ( \frac { 1 } { 2 a } \cos \theta ).$$

The action A and the measure [ dU ][ dφ ] in Eq.(2.31 are invaraiant under the U(1) gauge transformation rotating the local internal coordinate of complex field at each site x on the 4-dimensional lattice,

$$U _ { x \mu } \to \bar { V } _ { x + \mu } U _ { x \mu } V _ { x }, \ V _ { x } = e ^ { i \Lambda _ { x } } \in \mathbf U ( 1 ), \ \phi _ { x } \to V _ { x } \phi _ { x },.$$

## 2.4. Phase structure and Correlation functions

Phase structure of various gauge-field models has been investigated by both analytic and numerical methods. In particular, for the gauge-Higgs model, which describes the SC phase transitions, interesting phase diagrams have been obtained. To this end, the numerical study by the Monte-Carlo (MC) simulations plays a very important role as they provide us reliable results including all nonperturbative effects. Here we present some known schematic phase diagrams restricting the system to the relativistic Higgs coupling and the frozen radial degrees of freedom of φ x field,

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0089v2-figure-8a99af17-ded6-46c9-8eae-f951d33c694d.jpg)

$$c _ { 1 } ( \underset { x } { \overset { \phi } { x } } \bar { \underline { U } } _ { x \mu } \bar { \phi } _ { x + \mu } + \underset { x + \mu } { \overset { x + \nu } { x + \mu } } ) + c _ { 2 } \left ( \underset { x } { \overset { x + \nu } { \underline { U } } _ { x \mu } } + \underset { x } { \overset { x + \nu } { \underline { U } } } \right ) \\ \text{ illustration of the action A of the compact LGT (2.34). Open (filled) circle denot}$$

Fig. 3. Illustration of the action A of the compact LGT (2.34). Open (filled) circle denotes φ x = exp( iϕ x )( ¯ φ x ) and straight line with an arrow in the negative (positive) µ -direction denotes U xµ = exp( iθ xµ )( ¯ U xµ )

i.e., φ x = exp( iϕ x ), but with spatial dimensions d = 2 and 3 (dimension of the corresponding Euclidean lattice D = d +1 is 3 and 4). Explicitly, we consider the following partition function;

$$& \text{sources} \text{,} \\ & Z = \int [ d \varphi ] [ d \theta ] _ { G } \exp ( A ), \ A = c _ { 1 } \sum _ { x, \mu } \cos ( \varphi _ { x + \mu } + \theta _ { x \mu } - \varphi _ { x } ) + A _ { g } ^ { ( ^ { \prime } ) }, \\ & A _ { g } = c _ { 2 } \sum _ { x, \mu < \nu } \cos \theta _ { x \mu \nu }, \ A _ { g } ^ { \prime } = - \frac { c _ { 2 } } { 2 } \sum _ { x, \mu < \nu } \theta _ { x \mu \nu } ^ { 2 }, \ \theta _ { x \mu \nu } \equiv \theta _ { x \nu } + \theta _ { x + \nu, \mu } - \theta _ { x + \mu, \nu } - \theta _ { x \mu }, ( 2. 3 4 ) \\ & \text{s.s.s. \dots - \infty } \quad \ A \ \text{s.s.s. \dots - \infty } \quad \text{s.s.s. \dots - \infty } \quad \text{s.s.s. \dots - \infty } \quad \text{s.s.s. \dots - \infty } \quad \text{s.s.s. \dots - \infty } \quad \text{s.s.s. \dots - \infty } \quad \text{s.s.s. \dots - \infty } \quad \text{s.s.s. \dots - \infty } \quad \text{s.s.s. \dots - \infty } \quad \text{s.s.s. \dots - \infty } \quad \text{s.s.s. \dots - \infty } \quad \text{s.s.s. C }$$

where µ = 0 , · · · , d, and c 1 and c 2 are the parameters. A g is for the compact LGT and A ′ g is for the noncompact LGT. Each term of the action A is depicted in Fig.3. This model is sometimes called the Abelian Higgs model 11 . Note that the Higgs coupling (the c 1 term) along the imaginary time µ = 0 direction has the same form as the spatial directions, i.e., every direction has couplings both in the positive and negative directions. This is in contrast to the nonrelativistic model (2.31) in which the coupling for µ = 0 is only in the positive direction. In short, the coupling in the negative µ = 0 direction describes the existence of antiparticles which are allowed in the relativistic theory but not in nonrelativistic theory 27 .

Here we comment on the integration measure [ dθ ] G and gauge fixing. The path integral (2.34) involves the same contributions from a set of configurations of { ϕ, θ } (so called gauge copies) that are connected with each other by gauge transformations. So one may reduce this redundancy by choosing only one representative from each set of gauge copies without changing physical contents. This procedure is called gauge fixing and achieved by inserting the gauge fixing function G ( { ϕ, θ } ), which is not gauge invariant under Eq.(2.33), to the measure, [ dϕ dθ ][ ] → [ dϕ dθ ][ ] G ≡ G ( { ϕ, θ } )[ dϕ dθ ][ ]. The difference between the gauge fixed case [ dθ ] G and the nonfixed case [ dθ ] appears as a multiplicative factor in Z ,

$$Z _ { \text{nonfixed} } = C _ { \text{gv} } Z _ { \text{fixed} }, \ C _ { \text{gv} } = \prod _ { x, \mu } \int d \Lambda _ { x \mu },$$

where the constant C gv is called gauge volume. In the compact case, the region of θ xµ and Λ xµ are both compact [0 2 , π ). Therefore C gv is finite and the gauge fixing is irrelevant (either choice of fixing gauge or not will do). In the noncompact case, the region of θ xµ and Λ xµ is ( -∞ ∞ , ), so C gv diverges. Thus, one must fix the

gauge in formal argument in the noncompact case, [ dϕ dθ ][ ] G . We said 'in formal argument' here because the choice G ( { ϕ, θ } ) = 1 works even in the noncompact case as long as one calculate the average of gauge-invariant quantities, e.g., in the MC simulations. This is because the average does not suffer from the overall factor in Z . We shall present in Sect.4.1. another reason supporting this point, which is special in MC simulations.

Fig. 4. Schematic phase diagrams in the c 2 -c 1 plane of the 3D and 4D LGT with compact and noncompact action. All the Higgs-Coulomb transitions are of second order. The Higgs-confinement one in (d) is first order and the confinemnet-Coulomb one in (d) is weak first order. The partition function along the line c 2 = 0 is evaluated exactly by single link integral. It is an analytic function of c 1 and so no transitions exist along this line.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0089v2-figure-f1aa2d11-7fa2-4d6c-8481-b2a56d375c3e.jpg)

In Fig.4 we present the schematic phase diagrams in the c 2 -c 1 plane. The noncompact 3D and 4D models have the Coulomb phase and the Higgs phase. The compact 3D model has only the confinement phase. The compact 4D model has all the three phases. These points are understood as follows; (i) the Higgs phase may exist at large c 1 above the Coulomb phase in 3D and 4D, and above the confinement phase in 4D, (ii) the Coulomb phase may exist at sufficiently large c 2 in 4D case, (iii) the confinement phase exists in the compact case, but not in the noncompact phase.

We note that, as c 2 decreases, all the existing confinement-Higgs transition curves terminate at points with c 2 &gt; 0. This can be understood by calculating Z for c 2 = 0 exactly. Because the integration over θ xµ there is factorized as

$$Z ( c _ { 1 }, c _ { 2 } = 0 ) = \prod _ { x, \mu } \int d \theta _ { x \mu } ^ { \prime } \ e ^ { c _ { 1 } \cos \theta _ { x \mu } ^ { \prime } } = I _ { 0 } ( c _ { 1 } ) ^ { N _ { \ell } }, \ \theta _ { x \mu } ^ { \prime } \equiv \varphi _ { x + \mu } + \theta _ { x \mu } - \varphi _ { x }, \ ( 2. 3 6 )$$

where I 0 ( c 1 ) is the modified Bessel function and N /lscript ( ≡ ∑ x,µ 1) is the total number of links. Z c , c ( 1 2 = 0) is then an analytic function of c 1 and no phase transition can exists along the c 1 axis. This fact is known as complementarity and implies that the confinement phase and the Higgs phase are connected with each other without crossing a phase boundary.

Another limiting case c 2 → ∞ may imply ∆ θ → 0. Therefore the resulting system becomes

$$Z ( c _ { 1 }, c _ { 2 } \rightarrow \infty ) = \prod _ { x, \mu } \int d \varphi _ { x } ^ { \prime } \exp [ c _ { 1 } \sum _ { x, \mu } \cos ( \varphi _ { x + \mu } ^ { \prime } - \varphi _ { x } ^ { \prime } ) ],$$

by setting as θ xµ = ∇ µ Λ x (pure gauge configuration) with ϕ ′ x ≡ ϕ x + Λ . This x system is just the 3D and 4D XY spin models with the XY spin (cos ϕ , x sin ϕ x ), which are known to exhibit a second-order transition. This corresponds to the Higgs-Coulomb transition in Fig.4.

Let us discuss a general correlation function 〈 O ( { } { θ , ϕ } 〉 ) (the expectation value of the function O ( { } { θ , ϕ } )) of the model (2.34) defined as

$$\langle O ( \{ \theta \}, \{ \varphi \} ) \rangle = \frac { 1 } { Z } \int [ d \varphi ] [ d \theta ] O ( \{ \theta \}, \{ \varphi \} ) \exp ( A ).$$

First, 〈 O ( { } { θ , ϕ } 〉 ) is known to satisfy the Elitzur theorem 24 , which says

$$\langle O ( \{ \theta \}, \{ \varphi \} ) \rangle = 0 \text{ if } O ( \{ \theta \}, \{ \varphi \} ) \text{ is not gauge invariant.}$$

Its proof is rather simple and given by making a change of variables in a form of gauge transformation and noting the measure [ dθ ][ dϕ ] is invariant under it. A famous example of gauge-invariant correlation function is the Wilson loop defined on a closed loop C along successive links as

$$W ( C ) = \langle \prod _ { \ell \in C } U _ { \ell } \rangle.$$

It is used by Wilson as a nonlocal 'order parameter' to distinguish a confinement phase and a deconfinement phase for pure gauge system (i.e., involving no dynamical charged particles). If W C ( ) satisfies the 'area law' as C becomes large, W C ( ) ∼ exp( -aS C ( )) where a &gt; 0 and S C ( ) is the minimum area bound by C , the system is in the confinement phase. On the other hand, if it satisfies the perimeter law, W C ( ) ∼ exp( -a L C ′ ( )) where L C ( ) is the length of C , it is in the deconfinement phase. This interpretation comes from the fact that for the rectangular C of the horizontal size R and the vertical size T [The coordinate ( x , x 0 1 ) of four corners in the 0-1 plane, e.g., are (0,0), (0,R), (T,0), (T,R)], V ( R ) defined through W C ( ) ∼ exp( -V ( R T ) ) at sufficiently large T expresses the lowest potential energy stored by a pair of external charges separated by the distance R . The area law implies the confining potential V ( R ) = aR .

Second, one may estimate Z and the correlation functions of gauge-invariant objects by fixing the gauge because they are gauge invariant. Fixing implies to fix one degree of freedom per site x corresponding to Λ x whatever one wants.

For example, the choice ϕ x = 0 is called the unitary gauge, and another choice θ x 0 = 0 is called the temporal gauge. In the mean field theory for the compact LGT, one choose gauge variant quantities 〈 exp( iθ xµ ) 〉 and 〈 exp( iϕ x ) 〉 as a set of order parameters and predicts the phase structure involving the above mentioned three phases. They should be considered as a result after gauge fixing. If one wants to make it compatible with the Elitzur theorem, one simply need to superpose over the gauge copies of these order parameters as Drouffe suggested 25 . For 〈 exp( iθ xµ ) , 〉 it implies

$$\prod _ { x } \int _ { 0 } ^ { 2 \pi } \frac { d \Lambda _ { x } } { 2 \pi } \cdot \langle \exp ( i \theta _ { x \mu } ^ { \prime } ) \rangle = 0, \ \theta _ { x \mu } ^ { \prime } = \theta _ { x \mu } - \nabla _ { \mu } \Lambda _ { x }, \quad \quad ( 2. 4 1 )$$

without changing the gauge-invariant results such as phase diagram, etc..

## 3. Ferromagnetic Superconductivity and its Lattice Gauge Model

In this section and Sect.4, we shall consider a theoretical model for the FMSC, which can be regarded as a kind of LGT. In some materials such as UGe , it has 2 been observed that SC can coexist with a FM long-range order 29 . In UGe 2 and URhGe, a SC appears only within the FM state in the pressure-temperature ( P T -) phase diagram, whereas, in UCoGe, the SC exists both in the FM and paramagnetic states. Soon after their discovery, phenomenological models of FMSC materials were proposed 30 . In those studies, the FMSC state is characterized as a spin-triplet p -wave state of electron pairs as suggested experimentally 31 . In this section, we introduce the GL theory of FMSC materials defined on the 3d spatial lattice (4D Euclidean lattice). As we see, the FM order naturally induces nontrivial vector potential (gauge field) and it competes with another order, the SC order, therefore a nonperturbative study is required in order to clarify the phase diagram etc.

## 3.1. Ginzburg Landau theory of FMSC in the continuum space

Before considering a lattice model of FMSC itself, let us start with the GL theory of FMSC in the 3d continuum space 30 . It was introduced phenomenologically for the case with strong spin-orbit coupling, although the origin of SC, etc. had not been clarified. Its free energy density f GL contains a pair of basic 3d-vector fields, /vector ψ r ( ) and /vector m r ( ), which are to be regarded as c-numbers appearing in the pathintegral expression for the partition function of the underlying quantum theory as in Eq.(2.31) [We discuss the partition function in detail later; see Eq.(3.16)].

The three-component complex field /vector ψ = ( ψ , ψ 1 2 , ψ 3 ) t is the SC order parameter, describing the degrees of freedom of electrons participating in SC. More explicitly, /vector ψ r ( ) describes Cooper-pairs of electrons with total spin s = 1 (triplet) and the relative angular momentum /lscript = 1 ( p -wave) (Note that only odd /lscript 's are allowe for s = 1 by Pauli principle). The complex Cooper-pair field Φ ( q r ) for a s = 1 and /lscript = 1 state generally has 3 × 3 = 9 components ( q = 1 , · · · , 9). The strong spinorbit coupling forces the spin /vector S r ( ) and the angular momentum /vector L r ( ) made of Φ ( q r )

parallel each other, and reduces their degrees of freedom down to 3. Such a Cooper pair is directly expressed in terms of the amplitudes defined as

$$\Delta _ { \sigma \sigma ^ { \prime } } ( r ) \equiv \int d r _ { 0 } \Gamma ( r _ { 0 } ) \langle \hat { C } _ { \sigma } ( r + \frac { r _ { 0 } } { 2 } ) \hat { C } _ { \sigma ^ { \prime } } ( r - \frac { r _ { 0 } } { 2 } ) \rangle.$$

ˆ C σ ( r ) is the annihilation operator of electron at r and spin σ = ↑ ↓ , . Γ( r 0 ) with r 0 being the relative coordinate of electrons is the weight reflecting the attractive interaction and the p -wave nature. Because we assign /vector ψ r ( ) as a 3d vector, it is convenient to use the so called /vector d -vector, which transforms also as a vector in the 3d space;

$$\cdot \quad \vec { \psi } ( r ) \in \vec { d } ( r ) = \begin{pmatrix} d _ { x } ( r ) \\ d _ { y } ( r ) \\ d _ { z } ( r ) \end{pmatrix} \equiv \begin{pmatrix} - \frac { 1 } { 2 } ( \Delta _ { \uparrow \uparrow } ( r ) - \Delta _ { \downarrow \downarrow } ( r ) ) \\ - \frac { 1 } { 2 } ( \Delta _ { \uparrow \uparrow } ( r ) + \Delta _ { \downarrow \downarrow } ( r ) ) \\ \Delta _ { \uparrow \downarrow } ( r ) ( = \Delta _ { \downarrow \uparrow } ( r ) ) \end{pmatrix}. \quad \quad \ ( 3. 2 ) \\ \text{The road vector field $\vec{m}/r\}$ is the FM order parameter, discovery of}$$

The real vector field /vector m r ( ) is the FM order parameter, describing the degrees of freedom of electrons participating in the normal state. More explicitly, it describes their magnetization,

$$\vec { m } ( r ) = \langle \hat { C } ^ { \prime } _ { \sigma } ( r ) \Big ( - \, i \delta _ { \sigma \sigma ^ { \prime } } \vec { r } \times \vec { \nabla } + \frac { 1 } { 2 } \vec { \sigma } _ { \sigma \sigma ^ { \prime } } \Big ) \hat { C } ^ { \prime } _ { \sigma ^ { \prime } } ( r ) \rangle, \quad \quad ( 3. 3 ) \\ \hat { \sim } \, \omega \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdots$$

where ˆ C ′ σ ( r ) denotes the annihilation operator of electrons not participating in the SC state. Because there holds div /vector m r ( ) = 0 as in the usual magnetic field, /vector m r ( ) is expressed by using the vector potential /vector A r ( ) (gauge field) as

$$\vec { m } ( r ) = \text{rot} \vec { A } ( r ).$$

Therefore, the fundamental fields may be /vector ψ r ( ) and /vector A r ( ), The GL free energy density f GL is then given by ?

$$\text{$i$ me $ or $i$ree energy energy $JGL$ is $i$} \text{$given $y$} \\ f _ { \text{GL} } & = f _ { \psi } + f _ { m } + f _ { \text{Z} }, \\ f _ { \psi } & = K \sum _ { i } ( D _ { i } \vec { \psi } ) ^ { * } \cdot ( D _ { i } \vec { \psi } ) + \alpha ( T - T _ { \text{SC} } ^ { 0 } ) \vec { \psi } ^ { * } \cdot \vec { \psi } + \lambda ( \vec { \psi } ^ { * } \cdot \vec { \psi } ) ^ { 2 }, \\ f _ { m } & = K ^ { \prime } \sum _ { i } ( \partial _ { i } \vec { m } ) ^ { 2 } + \alpha ^ { \prime } ( T - T _ { \text{FM} } ^ { 0 } ) \vec { m } ^ { 2 } + \lambda ^ { \prime } ( \vec { m } ^ { 2 } ) ^ { 2 }, \\ f _ { \text{Z} } & = - J \vec { m } \cdot \vec { S }, \ D _ { i } = \partial _ { i } - 2 i e A _ { i }, \ \vec { S } = - i \vec { \psi } ^ { * } \times \vec { \psi }. \\ K ^ { ( \prime ) } \cdot \alpha ^ { ( \prime ) } \cdot \lambda ^ { ( \prime ) } \text{ are $GL$ numerates: real positive numerates characterizing each $ma$}.$$

K ( ) ′ , α ( ) ′ , λ ( ) ′ are GL parameters; real positive parameters characterizing each material. /vector S r ( ) is a real vector field describing the spin of Cooper pairs. For example, S 3 = 1 2 ( ∆ | ↑↑ | 2 -| ∆ ↓↓ | 2 ). f Z with J &gt; 0 is nothing but the Zeeman coupling of /vector S and /vector m , which enhances the coexistence of the FM and SC orders as one easily expects. Although the existence of f Z term is supported phenomenologically, its microscopic origin should be clarified by detailed study of the interactions of electrons (and nuclei) in each material.

f GL is gauge invariant under the following U(1) gauge transformation,

$$\vec { \psi } ( r ) \rightarrow \vec { \psi } ^ { \prime } ( r ) = e ^ { i \lambda ( r ) } \vec { \psi } ( r ), \ \vec { A } ( r ) \rightarrow \vec { A } ^ { \prime } ( r ) = \vec { A } ( r ) + \frac { 1 } { 2 e } \vec { \nabla } \lambda ( r ). \quad ( 3. 6 )$$

So /vector m r ( ) and /vector S r ( ) are gauge invariant. This gauge invariance reflects the original U(1) gauge invariance of the EM interaction of electrons. For example, the K term with the covariant derivative D i describes the interaction of SC Cooper pairs of charge -2 e with the vector potential /vector A made of the normal electrons. Its microscopic origin should be traced back to the repulsive Coulombic interactions between electrons.

Before introducing the lattice GL theory, let us list up some main characteristics suggested by the continuum GL theory (3.5).

- · In the mean-field approximation with ignoring the Zeeman coupling f Z , T 0 FM and T 0 SC are critical temperatures of the FM and SC phase transitions, respectively.
- · Existence of a finite magnetization, 〈 /vector m 〉 /negationslash = 0, means a nontrivial configuration of A r i ( ). This induces nontrivial spatial dependence of the Cooper pair /vector ψ through the kinetic term K D /vector ( i ψ ) ∗ · ( D /vector i ψ ) in f GL , because this term favors D /vector i ψ /similarequal 0. A typical example of such configurations is a vortex configuration characterized by a nonvanishing vorticity /vector v r ( ) of /vector ψ . The 3d vector /vector v r ( ) of a complex field φ r ( ) = | φ r ( ) | exp( iϕ r ( )) is defined generally as the winding number of its phase ϕ r ( ). Explicitly, the component of /vector v r ( ) along the direction of a normal vector /vector n is defined by the circle integral over ϕ r ( ),

$$\vec { n } \cdot \vec { v } ( r ^ { \prime } ) = \frac { 1 } { 2 \pi } \oint _ { C _ { \vec { n } } ( r ^ { \prime } ) } d \varphi ( r ) \ ( \in \mathbf Z ),$$

where C /vector n ( r ′ ) is a circle lying in the plane parpendicular to /vector n with its center at r ′ . /vector n · /vector v r ( ) takes integers due to the single-valuedness of φ r ( ). In the later sections, we shall consider two vorticities /vector v ± ( r ) corresponding to ψ 1 ( r ) ± iψ 2 ( r ) respectively. Then the simple mean-field like estimation of the SC critical temperature has to be reexamined by more reliable analysis.

- · The Zeeman coupling f Z obviously enhances the coexisting phase of the FM and SC orders, since it favors a set of antiparallel /vector S and /vector m with large | /vector ψ | and | /vector m | .

/negationslash

- · In the FM phase, /vector A produces a nonvanishing magnetic field /vector m r ( ) inside the materials. Because the K term in f ψ of Eq.(3.5) disfavors the coherent (spatially uniform) condensation of /vector ψ ∂ /vector ( i ψ = 0), due to rot A = 0 (favoring nontrivial ones such as vortices as explained above), the SC in a magnetic field is disfavored. In this sense, the system (3.5) is a kind of a frustrated system (the K term disfavors coexistence although the Zeeman term does).

Then a detailed study by using numerical methods is needed to obtain the correct phase diagram. To this end, we shall introduce a lattice version of the GL theory (3.5) that is suitable for the investigation by using MC simulations.

## 3.2. GL theory on the lattice

In this subsection we explain the lattice GL theory introduced in Ref. 32 . In constructing the theory, we start with the following two simplifications:

- (1) We consider the 'London' limit of the SC, i.e., we assume /vector ψ ∗ ( r ) · /vector ψ r ( ) = const . ignoring the fluctuations of the radial degrees of freedom | /vector ψ r ( ) . | This assumption is legitimate as the phase degrees of freedom /vector ψ r ( ) themselves play an essentially important role in the SC transition (Recall the account of phase coherence for a Bose-Einstein condensate) 28 .
- (2) We assume that the third component ψ 3 ( r ) ∝ ∆ ↑↓ ( r ) is negligibly small compared to the remaining ones ψ 1 2 , ( r ),

$$\psi _ { \pm } \equiv \frac { 1 } { \sqrt { 2 } } ( \psi _ { 1 } \pm i \psi _ { 2 } ), \, \psi _ { \uparrow \uparrow } \equiv \psi _ { - } \, \infty \, \Delta _ { \uparrow \uparrow }, \, \psi _ { \downarrow \downarrow } \equiv - \psi _ { + } \, \infty \, \Delta _ { \downarrow \downarrow }. \quad \ \ ( 3. 8 ) \\ \text{This simplification is consistent with the fact that the real materials exhibit}$$

This simplification is consistent with the fact that the real materials exhibit FM orders of Ising-type with the i = 3-direction as the easy axis. In fact, /vector S = -i /vector ψ ∗ × /vector ψ is calculated with ψ ↑↓ ≡ ψ 3 as

$$\vec { S } = ( - \sqrt { 2 } \, \text{Re} ( \psi _ { \uparrow \uparrow } + \psi _ { \downarrow } ) \psi _ { \uparrow \downarrow } ^ { * }, \sqrt { 2 } \, \text{Im} ( \psi _ { \uparrow \uparrow } - \psi _ { \downarrow } ) \psi _ { \uparrow \downarrow } ^ { * }, | \psi _ { \uparrow \uparrow } | ^ { 2 } - | \psi _ { \downarrow \downarrow } | ^ { 2 } ) ^ { t }. \quad ( 3. 9 ) \\ \intertext { \Sigma \, \text{the 7\,\text{common column $f$} \, \text{for\,\text{our\,as} \,\text{more\,as} } \, \text{more\,and} \, \text{to\,as} \, \text{for\,\text{our} } \, \vec { m } \, \sim }$$

So the Zeeman coupling f Z requires large ψ 1 2 , compared to ψ 3 to favor /vector m ∝ (0 , 0 , m ).

In summary, we parametrize the Cooper-pair field in the London limit, | ψ r ( ) | 2 = [ α/ (2 λ )]( T 0 SC -T ), at low T &lt; T ( 0 SC ) with the third easy axis as

$$\binom { \psi _ { 1 } ( r ) } { \psi _ { 2 } ( r ) } = \sqrt { \frac { \alpha } { 2 \lambda } ( T _ { \text{SC} } ^ { 0 } - T ) } \times \binom { z _ { 1 } ( r ) } { z _ { 2 } ( r ) } \,, \ | z _ { 1 } ( r ) | ^ { 2 } + | z _ { 2 } ( r ) | ^ { 2 } = 1. \quad ( 3. 1 0 ) \\ \tau _ { 1 } \, - \, \alpha _ { 1 } \, \dots \, \tau _ { 2 } \, - \, \dots \, \tau _ { 3 } \, - \, \tau _ { 4 } \, - \, \tau _ { 5 } \, - \, \tau _ { 6 } \, - \, \tau _ { 7 } \, - \, \tau _ { 8 } \, - \, \tau _ { 9 } \, - \, \tau _ { 9 } \, - \, \tau _ { 9 } \, - \, \tau _ { 9 } \, - \, \tau _ { 9 } \, - \, \tau _ { 9 } \, - \, \tau _ { 9 } \, - \, \tau _ { 9 } \, - \, \tau _ { 9 } \, - \, \tau _ { 9 } \, - \, \tau _ { 9 } \, - \, \tau _ { 9 } \, - \, \tau _ { 9 } \, - \, \tau _ { 9 } \, - \, \tau _ { 9 } \, - \, \tau _ { 9 } \, - \, \tau _ { 9 } \, - \, \ t }$$

The 'normalized' two-component complex field z a ( r ) ( a = 1 2) satisfying the con-, straint (3.10) is called CP 1 variable (CP stands for complex projective group). In terms of the SC order-parameter field z a ( r ), the first K term of (3.5) is rewritten as Kα (2 λ ) -1 ( T 0 SC -T ) i a D z i a · D z i a .

∑ ∑ Now let us introduce the lattice GL theory of FMSC defined on the 3d cubic lattice. Its free energy density per site f r is given by 32

$$\text{www.w3.org/wiki/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/h3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3//w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w1/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/g3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3\w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w4/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/W3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3/w3:/
f_r = - \frac { c _ { 1 } } { 2 } \sum _ { i = 1 } ^ { 3 } \sum _ { a = 1 } ^ { 2 } \frac { c _ { 3 } } { c _ { 4 } } \sum _ { r = 1 } ^ { 2 } \frac { c _ { 5 } } { c _ { 6 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 7 } } { c _ { 8 } } \sum _ { a = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac < c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { [ } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _{ i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 ] } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _{{ i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1. } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 3 } \frac { c _ { 9 } } { c _ { 9 } } \sum _ { i = 1 } ^ { 2 } \frac { c _ { 9$$

The five coefficients c i ( i = 1 ∼ 5) in (3.11) are real nonnegative parameters that are to distinguish various materials in various environments. z ra is the CP 1 variable put on the site r and plays the role of SC order-parameter field. U ri is the exponentiated

vector potential, θ ri , put on the link ( r, r + ). i /vector m r = ( m ,m r 1 r 2 , m r 3 ) t is the magnetic field made out of θ ri as

$$m _ { r i } \equiv \sum _ { j, k = 1 } ^ { 3 } \epsilon _ { i j k } \nabla _ { j } \theta _ { r k }, \ \nabla _ { j } \theta _ { r k } \equiv \theta _ { r + j, k } - \theta _ { r k }.$$

/vector m r serves as the FM order-parameter field. Following Eq.(3.5), the spin /vector S r of the Cooper pair is defined as

$$\vec { S } _ { r } = ( 0, 0, S _ { r 3 } ), \ S _ { r 3 } \equiv - i ( \bar { z } _ { r 1 } z _ { r 2 } - \bar { z } _ { r 2 } z _ { r 1 } ) \ ( \in | \psi _ { \uparrow \uparrow } | ^ { 2 } - | \psi _ { \downarrow \downarrow } | ^ { 2 } ), \quad ( 3. 1 3 )$$

where we absorbed the normalization of /vector S r ( ) into the coefficient c 3 . /vector m r , /vector S r and f r are invariant under the following gauge transformation,

$$z _ { r a } \to z _ { r a } ^ { \prime } = e ^ { i \lambda _ { r } } z _ { r a }, \ U _ { r i } \to U _ { r i } ^ { \prime } = e ^ { - i \lambda _ { r + i } } \ U _ { r i } \ e ^ { i \lambda _ { r } }. \quad \quad ( 3. 1 4 )$$

Here we comment on the way of putting a 3d continuum vector field on the 3d lattice. It sounds natural that a vector field /vector B r ( ) should be put on the link ( r, r + ) as i B ri . However, it is too naive. In fact, the gauge field /vector A r ( ) has the nature of connection as explained and the connection between two points separated by finite distance such as nearest-neighbor pair of sites cannot be implemented by θ ri itself but by its exponentiated form such as U ri . This is the reason why one has ¯ z r + i,a ¯ U ri z ra term in Eq.(3.11). On the other hands, concerning to the SC order field /vector ψ r ( ), its suitable lattice version depends on the coherence length ξ . If ξ is of the same order of lattice spacing a , its vector nature should be respected and the link field ψ ri is adequate. If ξ /greatermuch a , the detailed lattice structure is irrelevant and a simpler version ψ ra ( a = 1 2 3) on the site (or its Ising counterpart , , z ra ( a = 1 2)) , will do. Eq.(3.11) corresponds to the latter case ξ /greatermuch a as some materials show.

Let us see the meaning of each term in f r . The c 1 -term describes a hopping of Cooper pairs from site r to r + i (and from r + i to r ). As the Cooper pair has the electric charge -2 , it minimally couples with the vector potential e θ ri via U ri as explained above. It is important to observe the relation between the hopping parameter c 1 and T . From the expression given in the paragraph below Eq. (3.5), we have

$$c _ { 1 } \sim K \alpha \lambda ^ { - 1 } ( T _ { S C } ^ { 0 } - T ) a.$$

At sufficiently large βc 1 = c 1 /T (we set k B = 1) that corresponds to low T 's, the c 1 -term stabilizes the combination of phases of ¯ z r + i,a ¯ U ri z ra , and then, if U ri is stabilized already by other mechanism, a coherent condensation of the phase degrees of freedom of z r is realized inducing the superconductivity . The c 2 and c 4 -terms are the quartic GL potential of /vector m r , favoring a finite amount of local magnetization 〈 /vector m r 〉 /negationslash = 0 (note that we take c 2 &gt; 0). We should remark that these terms controls intrinsic magnetization, which is different from the fluctuating but external magnetic field. The c 5 -term represents the NN coupling of the magnetization and corresponds to the K ′ term of Eq.(3.5). It enhances uniform configurations of /vector m r , i.e., a FM longrange order signaled by a finite magnetization m m , 2 ≡ lim | r -r ′ |→∞ 〈 /vector m r · /vector m r ′ 〉 /negationslash = 0.

The c 3 -term is the Zeeman coupling, which favors collinear configurations of /vector m r and /vector S r , i.e., enhances the coexistence of FM and SC orders.

The partition function Z at T for the energy density f r of Eq.(3.11) is given by the integral over a set of two fundamental fields z ra and θ ri as

$$\text{integral over a set or two runamental neas } z _ { r a } \text{ and } \nu _ { r i } \text{ as} \\ Z & = \int [ d z ] [ d \theta ] \exp ( - \beta F ), \ \beta = T ^ { - 1 }, \ F = \sum _ { r } f _ { r }, \\ [ d z ] & = \prod _ { r } d ^ { 2 } z _ { r 1 } d ^ { 2 } z _ { r 2 } \ \delta ( | z _ { r 1 } | ^ { 2 } + | z _ { r 2 } | ^ { 2 } - 1 ), \\ [ d \theta ] & = G ( \{ \theta \} ) \prod _ { r, i } d \theta _ { r i }, \ \theta _ { r i } \in ( - \infty, \infty ). \\ \text{re} C ( I \theta ) \, \text{ is a square fusion function} [ \succeq the naraarh contain} \, F ( 9 3 ) ]$$

where G ( { } θ ) is a gauge fixing function[see the paragraph containing Eq.(2.35)].

We stress here that the free energy F and the integration variables in Eq.(3.16) have no dependence on the imaginary time x 0 in contrast with the expression (2.34). The ordinary GL theory explained in the literature also shares this properties. This implies that we ignore the x 0 -dependent modes of would-be ϕ x and θ xµ as A ( { ϕ x } { , θ xµ } ) → βF ( { ϕ r } { , θ ri } ). This is an approximation applicable for small β 's, i.e., for high T 's. This procedure allows us to make use of MC method that is based on the probabilistic process. In fact, the form (3.16) has a real function F and the Boltzmann factor exp( -βF /Z ) is able to be interpreted as a probability. This is in contrast with the original action A corresponding to F ; A is certainly a complex function as Eq.(2.34) giving rise to a complex probability. We come back to this point later.

The coefficients c i ( i = 1 , · · · , 5) in f r may have nontrivial T -dependence as Eqs.(3.5) and (3.15) suggest. However, in the present study we consider the response of the system by varying the 'temperature' T ≡ 1 /β defined by β , an overall prefactor in Eq.(3.16), while keeping c i fixed. This method corresponds to well-known studies such as the FM transition by means of the O (3) nonlinearσ model 33 and the lattice gauge-Higgs models discussed in Sec.2.1, and is sufficient to determine the critical temperature [See, e.g., Eq.(4.7)].

In the following section, we shall show the results obtained by MC numerical evaluation of Eq.(3.16). Physical quantities O ( z, θ ) like the internal energy, specific heat, correlation functions, etc. are also calculated by the MC simulations as

$$\langle \mathcal { O } ( z, \theta ) \rangle = \frac { 1 } { Z } \int [ d z ] [ d \theta ] \mathcal { O } ( z, \theta ) \exp ( - \beta F ).$$

From the obtained results, we shall clarify the phase diagram of the system and physical properties of each phase.

## 4. Results of Monte Carlo simulation

In the present section, we introduce and discuss the results of numerical calculations 32 34 , for the lattice GL model (3.11), which include the phase diagram

and Meissner effect in the SC phase. Figs.5 ∼ 9 presented below are from Ref. 32 and Fig.10 ∼ 13 are from Ref. 34 .

## 4.1. Phase diagram and Meissner effect

We first seek for a suitable boundary condition (BC) of θ ri as the magnetization is expressed by θ ri as Eq.(3.12). A familiar condition is the periodic boundary condition such as θ r + Li,j = θ rj for the linear system size L . However, this condition necessarily gives rise to a vanishing net mean magnetization due to the lattice Storkes theorem. So we consider the 3d lattice of the size (2 + L +2) 2 × L and take a 'free BC' on z r in the i = 1 2 directions defined by ,

$$z _ { r + i, a } - \bar { U } _ { r i } z _ { r a } = 0 \text{ for } r = \begin{cases} ( 0, r _ { 2 }, r _ { 3 } ), & i = 1, \\ ( L + 2, r _ { 2 }, r _ { 3 } ), i = 1, \\ ( r _ { 1 }, 0, r _ { 3 } ), & i = 2, \\ ( r _ { 1 }, L + 2, r _ { 3 } ), i = 2, \end{cases} & ( 4. 1 ) \\ \text{hereas we impose the free boundary condition on } \theta _ { r i }. \text{ It is easily shown that the amount } i \text{$\mathbb{C}$}$$

whereas we impose the free boundary condition on θ ri . It is easily shown that the BC (4.1) implies that the suppercurrent j SC ri ,

$$j _ { r i } ^ { \text{SC} } \in \text{Im} \ ( \sum _ { a } \bar { z } _ { r + i, a } \bar { U } _ { r i } z _ { r a } ),$$

satisfies j SC r 1(2) = 0 on the boundary surfaces in the 2(1) -3 planes. That is, the supercurrent does not leak out of the SC material. Another possible BC is the one that imposes vanishing Cooper-pair amplitude z ra , but we expect, as usual, that the qualitative bulk properties of the results are not affected by the BC seriously. [The above BC (4.1) means that the magnetization is just like the Ising type and 〈 /vector m x 〉 = (0 0 , , m 3 ), which describes the real materials properly.]

In the simulation made in Refs. 32 34 , we use the standard Metropolis algorithm 35 for the lattice size up to L = 30. The typical number of sweeps for a measurement is (30000 ∼ 50000) × 10 and the acceptance ratio is 40% ∼ 50%. The error is calculated as the standard deviation of 10 samples obtained by dividing one measurement run into 10 pieces. Concerning to the gauge-fixing term of Eq.(2.35), even for a noncompact gauge theory such as the present one, gauge fixing is not necessary in practical MC simulations. This is because, in addition to the reasons given at the end of paragraph of Eq.(2.35), the MC updates using random numbers almost never generate exactly gauge-equivalent configurations.

We calculate the internal energy U , the specific heat C of the central region R of the size L 3 , the magnetization m i , which are defined as follows,

$$U = \frac { 1 } { L ^ { 3 } } \langle F _ { L } \rangle, \ C = \frac { 1 } { L ^ { 3 } } \langle ( F _ { L } - \langle F _ { L } \rangle ) ^ { 2 } \rangle, \ F _ { L } \equiv \sum _ { r \in R } f _ { r }, \ m _ { i } \equiv \frac { 1 } { L ^ { 3 } } \langle | \sum _ { r \in R } m _ { r i } | \rangle, \ \ ( 4. 3 )$$

and the normalized correlation functions,

$$G _ { m } ( r - r _ { 0 } ) = \frac { \langle \vec { m } _ { r } \cdot \vec { m } _ { r _ { 0 } } \rangle } { \langle \vec { m } _ { r _ { 0 } } \cdot \vec { m } _ { r _ { 0 } } \rangle }, \, G _ { S } ( r - r _ { 0 } ) = \frac { \langle S _ { r 3 } S _ { r _ { 0 }, 3 } \rangle } { \langle S _ { r _ { 0 }, 3 } S _ { r _ { 0 }, 3 } \rangle }, \quad \quad ( 4. 4 )$$

where r 0 is chosen on the boundary of R such as (3 , 2+ L/ , z 2 ). Singular behaviors of U and/or C indicate the existence of phase transitions, and the magnetization m i and the correlation functions in Eqs.(4.4) clarify the physical meaning of the observed phase transitions. As the present system is a frustrated system as we explained above, G m and G S may exhibit some peculiar behavior in the FMSC state (the state in which FM and SC orders coexist).

To verify that the model (3.16) actually exhibits a FM phase transition as T is lowered, we put c 1 = c 3 = 0 and ( c 2 , c 4 , c 5 ) = (0 5 4 0 1 0), . , . , . and measured C increasing β in the Boltzmann factor of Eq.(3.16). In Fig.5 we show C and m i , which obviously indicate that a second-order phase transition to the FM state takes place at β FM = 1 /T FM /similarequal 2 0. We observed that other cases with various values of . c 2 4 5 , , exhibit similar FM phase transitions.

Let us see the SC phase transition and how it coexists with the FM. We recall that the case of all c i = 0 except for c 1 was studied in the previous paper Ref. 19 . There it was found that the phase transition from the confinement phase to the Higgs phase takes place at c 1 /similarequal 2 85. This result suggests that the SC state exists . at sufficiently large c 1 also in the present system with c i (=1) = 0.

/negationslash

/negationslash

To study this possibility, we use ( c 2 , c 4 , c 5 ) = (0 5 . , 4 0 . , 1 0) as in the FM sector . above and ( c 1 , c 3 ) = (0 2 . , 0 2). In Fig.6a, we show . C vs β , which exhibits a large and sharp peak at β /similarequal 2 1 and a small and broad one at . β /similarequal 4 5. To understand the . physical meaning of the second broad peak, it is useful to measure 'partial specific heat' C i for each term F i in the free energy (3.11) defined by

$$C _ { i } = \frac { 1 } { L ^ { 3 } } \langle ( F _ { i } - \langle F _ { i } \rangle ) ^ { 2 } \rangle, \ F _ { i } = \sum _ { r \in R } f _ { i r },$$

where f ir is the c i -term in f r of Eq.(3.11).

Figs.6b,c show that the partial specific heat C 2 4 5 , , have a sharp peak at β /similarequal 2 1. . Thus the peak of total C there should indicate the FM phase transition. The analysis of Meissner effect (see the discussion below using Fig.8) supports this point. On the other hand, C 1 of the c 1 -term has a relatively large and broad peak at β /similarequal 4 5. . C 3 also shows a broader peak there. Then we judge that the SC phase transition

Fig. 5. (a) Specific heat for ( c 2 , c 4 , c 5 ) = (0 5 . , 4 0 . , 1 0) and . c 1 = c 3 = 0. At β /similarequal 2 0, . C exhibits a sharp peak indicating a second-order FM phase transition. (b) Each component of magnetization m µ vs β . For T &lt; T FM , m 3 develops, whereas m 1 and m 2 are zero within the errors as expected.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0089v2-figure-d06e8cda-3a84-4c88-b501-c1d6e941e020.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0089v2-figure-e9b7d581-31c2-4a95-be11-cb3f52eae1b1.jpg)

## 24 I.Ichinose and T.Matsui

Fig. 6. (a) Specific heat vs β for ( c 1 , c 3 ) = (0 2 0 2) and ( . , . c 2 , c 4 , c 5 ) = (0 5 4 0 1 0) ( . , . , . L = 20). There are a large peak at β /similarequal 2 1 and a small one at . β /similarequal 4 5. . (b,c) Partial specific heat C i of Eq.(4.5) vs β . The small and broad peak at β /similarequal 4 5 in . C is related to fluctuations of c 1 -term.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0089v2-figure-d0481072-c341-468c-be6b-34b3147c842d.jpg)

/negationslash takes place at β SC /similarequal 4 5. To support these conclusions, we show . G m ( r ) and G S ( r ) in Fig.7. At β = 2 5, . G m ( r ) exhibits a finite amount of the FM order, whereas G S ( r ) decreases very rapidly to vanish. This means that, as T is decreased, the FM transition takes place first and then the SC transition does. Therefore, for β ≥ β SC /similarequal 4 5, the FM and SC orders coexist. In this way, the partial specific heat . may be used to judge the nature of each transition found by the peak of full specific heat C (We note C = i C i in general due to interference).

∑ As we explained in the previous section, the bare transition temperature T 0 SC in Eq.(3.5) and the genuine transition temperature T SC are different. Then it is interesting to clarify the relation between them. From Eq.(3.16), any physical quantity is a function of βc i . In the numerical simulations, we fix the values of c i and vary β as explained. Then the result β SC /similarequal 4 5 means .

Fig. 7. Correlation functions G m ( r ) and G S ( r ) at various T 's for L = 20. c i 's are the same as in Fig.6.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0089v2-figure-aac26530-69c5-430a-a8a7-b5cbe8f0a747.jpg)

Fig. 8. Gauge-boson mass M G of the external magnetic field propagating in x y -plane vs β for the same c i as in Fig.6 and c 2 ′ = 3 0 ( . L = 16). At the SC phase transition point β SC /similarequal 4 5 (indicated . by an arrow) determined by the peak of C M , G starts to increase from small values.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0089v2-figure-92fe004d-ebb6-465d-b53f-bd57cd29e661.jpg)

By using Eq.(3.15), this gives the following relation;

$$\frac { 1 } { T _ { \text{SC} } } \frac { K \alpha ( T _ { \text{SC} } ^ { 0 } - T _ { \text{SC} } ) a } { \lambda } = 0. 9 0 \ \to T _ { \text{SC} } = \left ( 1 + \frac { 0. 9 0 \ \lambda } { K \alpha \ a } \right ) ^ { - 1 } T _ { \text{SC} } ^ { 0 }. \quad ( 4. 7 ) \\ \quad \ \times \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \dots \ \delta \ \ \end{array}$$

Eq.(4.7) shows that the transition temperature is lowered by the fluctuations of the phase degrees of freedom of Cooper pairs. We expect that a relevant contribution to lowering the SC transition temperature comes from vortices that are generated spontaneously in the FMSC as we shall show in Sec.4.2.

After having confirmed that the genuine critical temperature can be calculated by the critical value of β with fixed c i , we use the word temperature in the rest of the paper just as the one defined by T ≡ 1 /β while c i are T -independent parameters.

One of the most important phenomena of the SC is the appearance of a finite mass of the electromagnetic field, i.e., the Meissner effect. In theoretical study of the SC that includes the effect of fluctuations of the Cooper-pair wave function as in the present one, a finite mass of the photon is the genuine order parameter of the SC. To study it, we follow the following steps 19 ; (i) introduce a vector potential θ ex ri for an external magnetic field, (ii) couple it to Cooper pairs by replacing U ri → U ri exp( iθ ex ri ) in the c 1 term of f x and add its magnetic term f ex r = + c 2 ′ ( /vector m ex r ) 2 ( c ′ 2 &gt; 0) to f r with /vector m ex r defined in the same way as (3.12) by using θ ex ri , (iii) let θ ex ri fluctuate together with z xa and θ ri and measure an effective mass M G of θ ex ri via the decay of correlation functions of /vector m ex r . The result of /vector m ex r propagating in the 1-2 plane is shown in Fig.8. It is obvious that the mass M G starts to develop at the SC phase transition point, and we conclude that Meissner effect takes place in the SC state.

## 4.2. Vortices in FMSC state

In this section, we study the FMSC state observed in the previous section in detail. In particular, we are interested in whether there exist vortices and their density if any. There are two kinds of vortices as the present SC state contains two gaps described by ψ + r ∝ z + r and ψ -r ∝ z -r , where z ± r are defined as

$$\cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad\cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \ \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cd. \quad \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cd. \quad \quad \quad$$

Corresponding to the above Cooper pair fields, one may define the following two kinds of gauge-invariant vortex densities V + r and V -r in the 1-2 plane;

$$& \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \alpha \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xright ], \quad ( 4. 9 ) \\ & \quad - \mod \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad$$

where mod( x ) ≡ mod( x, 2 π ). In short, V ± r describes vortices of electron pairs with the amplitude z ± r .

To study the behaviors of vortices step by step, we first consider the case of constant magnetization (magnetic field), i.e., we set θ ri by hand so that

$$m _ { r 1 } = m _ { r 2 } = 0, \ m _ { r 3 } = f = \text{constant},$$

freezing their fluctuations. As the above numerical studies show that a typical magnitude of the magnetization 〈 m r 3 〉 = 0 8 . · · · , we put f = π 4 . In this case, the SC phase transition takes place at β = 4 8. Snapshots of vortices are shown in Fig.9. . It is obvious that at lowerT , i.e., β = 7 0, densities of vortices are low compared to . the higherT ( β = 3 0) case but they are still nonvanishing. This result obviously . comes from the existence of the finite magnetization, i.e., the internal magnetic field.

Fig. 9. Snapshots of vortex densities V ± r at ( c 1 , c 2 , c 3 , c 4 , c 5 ) = (0 2 0 5 0 2 . , . , . , 4 0 . , 1 0) for a . fixed magnetization /vector m r = (0 0 , , π/ 4) ( L = 16) from Ref. 32 . Black dots; V ± r = 0.875, Dark gray dots; -1.125, Light gray dots; -0.125. (A) V + r at β = 3 0, (B) . V -r at β = 3 0, (C) . V + r at β = 7 0, (D) . V -r at β = 7 0. The average magnitude . 〈| V r ± |〉 is (A) 0.387, (B) 0.380, (C) 0.331 and (D) 0.335. The points V ± r = -0 125 = . -m / 3 (2 π ) reflect /vector m itself, and corresponds to the state without genuine vortices.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0089v2-figure-e09f0687-c034-475c-83c7-07143fb7a70a.jpg)

Fig. 10. Snapshot at β = 2 0 . with dynamically fluctuating /vector m r . The values of c i are ( c 1 , c 2 , c 3 , c 4 , c 5 ) = (0 2 . , 0 5 . , 0 2 . , 4 0 . , 1 0). .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0089v2-figure-781e57fa-e0c5-45ab-8980-c73484188333.jpg)

It is interesting to study behavior of vortices and the magnetization in the FMSC state rather in detail. We now allow for fluctuations of θ ri . In this state, the magnetization m r 3 fluctuates around its mean value. As we showed in the previous section, the FM and SC phase transitions take place at β FM /similarequal 2 1 and . β SC /similarequal 4 5, . respectively. We are interested in the β -dependence of various quantities such as the magnetic field 〈 m r 3 〉 , the vortex densities V ± r , and the Cooper-pair densities 〈 ρ ± r 〉 .

In Figs. 10 ∼ 13 we show these quantities as snapshots at β = 2 0 3 0 4 0 . , . , . , 5 2 for . the set of parameters ( c 1 , c 2 , c 3 , c 4 , c 5 ) = (0 2 . , 0 5 . , 0 2 . , 4 0 . , 1 0). At . β = 2 0 (Fig.10), . the local magnetization m r 3 fluctuates rather strongly, and also the density of Cooper pairs ρ + r and ρ -r are almost equal. Densities of vortices V ± r are also fairly large. At β = 3 0(Fig.11), . m r 3 tends to have a positive value. This corresponds to the fact that the system has the FM long-range order. At β = 4 0(Fig.12), . ρ + r

Fig. 11. Snapshot at β = 3 0 with . dynamically fluctuating /vector m r . The values of c i are the same as in Fig.10.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0089v2-figure-f2d19732-99c7-4169-961c-574ac76ec4cf.jpg)

## 28 I.Ichinose and T.Matsui

- density

Fig. 12. Snapshot at β = 4 0 with . dynamically fluctuating /vector m r . The values of c i are the same as in Fig.10.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0089v2-figure-224da6a3-b298-4409-998a-4ed9c6148763.jpg)

is getting larger compared to ρ -r . Similarly, V -r is sightly larger than V + r . These results reflect the Zeeman coupling. At β = 5 2(Fig.13), vortex-antivortex pairs are . formed, and V -r is larger than V + . As in the case of the constant magnetic field, r the densities of vortices are nonvanishing even in the SC phase.

## 5. Conclusion

In this paper, we present a brief review of LGT with Higgs matter field, and then apply it to the GL lattice model of FMSC. The detailed MC simulations of the model reveal the global structure such as the phase diagram and the physical characteristics of each phase in an explicit manner. The concepts and knowledges of LGT considerably help us to build an appropriate lattice model as well as to interpret the results of MC simulations. We want to add this example to the rather

Fig. 13. Snapshot at β = 5 2 with . dynamically fluctuating /vector m r . The values of c i are the same as in Fig.10.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0089v2-figure-a64bfce2-4453-4dec-a881-c9cf10a8630f.jpg)

long list of studies that have utilized the following common three-step method: (i) Establishing a lattice gauge model of the phenomenon of one's interest, (ii) Performing MC simulations of the model, (iii) Intepreting the MC results according to LGT. Here we list up a couple of tips to remove obstacles in carrying out this program.

- · MC simulations become more effective with less variables. To avoid unnecessary inflation of degrees of freedom, one should identify what are the relevant variables. Neglecting the radial degrees of freedom of the Higgs field φ r is an example 28 .
- · MC simulations require positive measure to make use of a probabilistic (Markoff) process. Fermionic systems in their original forms suffer from the well known negative-sign problem, which generates negative probabilities. By restricting the region of parameters of the model, it could be possible to avoid this problem. For example, the critical region around a continuous phase transition, the order-parameters are small and the GL expansion may assure us a positive provability 6 .
- · Another example to avoid nonpositive provability is the case of bosonic system at finite temperatures. The path integral of bosonic variables usually contain the imaginary kinetic term like i ∫ dτφ dφ/dτ ¯ ( ) in the action, where τ ∈ (0 , β = 1 / k ( B T )) is the imaginary time. By focusing in the high-temperature region, the τ -dependent modes of φ τ ( ) may be ignored and so is this kinetic term. This is the case of FMSC studied in Sect. 3 and 4, and some of the cases studied in Ref. 6 7 . , This procedure has an extra merit that the dimension of the effective lattice system is reduce to d -dimensions instead of D = d +1 dimensions due to the absence of the x 0 -direction.

Nowadays, it is well known that a quite wide variety of systems of cold atoms put on an optical lattice may serve as a quantum simulator of certain definite quantum model with almost arbitrary values of interaction parameters. Such a simulator is to clarify the dynamical (time-dependent) properties of the quantum model. Therefore, the importance of theoretical analysis of such models become increased as a guide for experiments. For example, information on the global phase structure of a quantum model is quite useful to perform experiments in an effective manner, although the former is confined to the static (time-independent) properties in thermal equilibrium. We expect that the above path (i-iii) with the successive tips may be useful for this purpose, in particular to simulate the model including compact or noncompact gauge fields. We hope as many readers as possible make use of this practical but powerful tool to strengthen and expand their researches.

## Acknowledgments

We thank Mr. T. Noguchi, our coauthor of Ref. 34 , for his help and discussions. We also thank Dr. K. Kasamatsu for reading the manuscript and presenting comments and suggestions. This work was supported by JSPS KAKENHI Grantnumbers 23540301, 26400246, 26400412.

## References

- 1. See, e.g. a review; P. A. Lee, N. Nagaosa, and X-G Wen, Rev. Mod. Phys. 78 , 17 (2006).
- 2. P. A. Lee, Science 321, 1306 (2008).
- 3. M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. S. Fisher, Phys. Rev. B 40 , 546 (1989); D. Jaksch, C. Bruder, J. I. Cirac, C. W. Cardiner, and P. Zoller, Phys. Rev. Lett. 81 , 3108 (1998).
- 4. M. W. Long and X. Zotos, Phys. Rev. B 45 , 9932 (1992); M. Boninsegni, Phys. Rev. Lett. 87 , 087201 (2001); Phys. Rev. B 65 , 134403 (2002); L.-M. Duan, E. Demler, and M. D. Lukin, Phys. Rev. Lett. 91 , 090402 (2003); E. Altman,W. Hofstetter, E. Demler, andM. D. Lukin, New J. Phys. 5 , 113 (2003). For the relation between the Bose Hubbard model and the bosonic t -J model, see Appendix A of the last reference of Ref. 7 .
- 5. See, e.g., M. Lewenstein, A. Sanpera, and V. Ahufinger, ' Ultracold Atoms in Optical Lattices: Simulating Quantum Many-body Systems ', (Oxford University Press, 2012).
- 6. For the t -J model of electrons, see G. Baskaran and P. W. Anderson, Phys. Rev. B 37 , 580 (1988); A. Nakamura and T. Matsui, Phys. Rev. B 37 , 7940 (1988); I. Ichinose and T. Matsui, Phys. Rev. B 45 , 9776 (1992); Nucl. Phys. B 394 , 281 (1993); Phys. Rev. B 51 , 11860 (1995); P. A. Lee and N. Nagaosa, Phys. Rev. B 46 , 5621 (1992); N. Nagaosa and P. A. Lee, Phys. Rev. Lett. 64 , 2450 (1990); Phys. Rev. B 43 , 1233 (1991). I. Ichinose, T. Matsui and M. Onoda, Phys. Rev. B 64 , 104516 (2001).
- 7. K. Aoki, K. Sakakibara, I. Ichinose, and T. Matsui, Phys. Rev. B 80 , 144510 (2009). Y. Nakano, T. Ishima, N. Kobayashi, K. Sakakibara, I. Ichinose, and T. Matsui, Phys. Rev. B 83 , 235116 (2011).
- 8. K. G. Wilson, Phys. Rev. D 10 , 2445 (1974).
- 9. For earlier reviews and books of LGT, see, e.g., J. B. Kogut, Rev. Mod. Phys. 51 , 659 (1979); M. Creutz, ' Quarks, gluons and lattices ', (Cambridge University Press 1985); I. Montvay and G. M¨nster, ' u Quantum Fields on a Lattice ', (Cambridge University Press 1997).
- 10. G. 't Hooft, Nucl. Phys. B 138 ,1 (1978); M. E. Peskin, Ann. Phys. 113 , 122 (1978).
- 11. M. B. Einhorn and R. Savit, Phys. Rev. D 17 , 2583 (1978).
- 12. G. Kotliar and A . Ruckenstein, Phys. Rev. Lett. 57 , 2790 (1987); Z. Zou and P. W. Anderson, Phys. Rev. B 37 , 627 (1988).
- 13. P. W. Anderson, Phys. Rev. Lett. 64 , 1839 (1990).
- 14. See, e.g., J. K. Jain, ' Composite fermions ', (Cambridge University Press, 2007) and references cited therein.
- 15. For composite boson formalism, see M. P. A. Fisher and D. H. Lee, Phys. Rev. B 39 , 2756 (1989); S. C. Zhang, Int. J. Mod. Phys. B 6 , 25 (1992). For composite fermion, formalism, see J. Jain, Phys. Rev. Lett. 63 , 199 (1989), Phys. Rev. B 40 , 8079 (1989); 41 , 7653 (1990).
- 16. J. Kogut and L. Susskind, Phys. Rev. D 11 , 3376 (1974).
- 17. For inclusion of relativistic fermions, see Refs. 8 9 16 , , For relativistic bosons, see Sect.2.4 and Ref. 11 . For nonrelativistic fermions, see I. Ichinose and T. Matsui, Phys.

## Rev. B 45 , 9776 (1992).

- 18. K. G. Wilson and J. Kogut, Phys. Rep. 12 , 75 (1974).
- 19. S. Takashima, I. Ichinose and T. Matsui, Phys. Rev. B 72 , 075112 (2005).
- 20. In the relativistic theory, the charge-conjugate pair of φ -particle (antiparticle) having the charge -g should appear on an equal footing.
- 21. A. M. Polyakov, Phys. Lett. B 59 , 82 (1975).
- 22. Sometimes the noncompact model is studied just for the purpose of comparison with the compact model. In general, one must ask for every system and situation under consideration which model (compact vs. noncompact) is more suitable.
- 23. In some literatures, the term | φ x + µ -U xµ φ x | 2 appears instead of | φ x + µ -¯ U xµ φ x | 2 . They are related by φ x ↔ ¯ φ x or U xµ ↔ ¯ U xµ and equivalent each other.
- 24. S. Elitzur, Phys. Rev. D 12 , 3978 (1975).
- 25. J. M. Drouffe, Nucl. Phys. B 170 , 211 (1980).
- 26. See, e.g., E. Sa'nchez-Velasco, Phys. Rev. E 54 , 5819 (1996), and references cited therein.
- 27. In the relativistic theory, Eq.(2.34) is obtained by starting with a set of annihilation operator, ˆ a r of particle and ˆ b r of antiparticle. Then, the charge operator is written as q a (ˆ † r ˆ a r -ˆ ˆ b † r b r ). One writes the path-integral expression of Z by using the canonical pair ( ˆ φ ± r , π ˆ ± r ) with ˆ a r = ˆ φ + r + ˆ iπ + r , ˆ b r = ˆ φ -r + ˆ iπ -r . Then by integrating over π ± r ∈ C one arrives at Eq.(2.34) with φ r = φ + r + iφ -r ∈ C . This φ r is a canonical field in contrast with the coherent-state field φ r in Eq.(2.31) although the same symbol φ r is used.
- 28. Of course, the fluctuations of radial component | φ x | may give rise to certain effects. For example, in the abelian Higgs model, fluctuations of | φ x | may change the secondorder SC transition to a first-order one if the coefficient λ of the | φ x | 4 coupling [See Eq.(2.22)] becomes sufficiently small. See, e.g., K. Jansen, J. Jers´k, C. B. Lang, T. a Neuhaus, G. Vones, Nucl. Phys. B 265 , 129 (1986).
- 29. For UGe , see S. S. Saxena 2 et al. , Nature (London) 406 , 587 (2000); A. Huxley et al. , Phys. Rev. B 63 , 144519 (2001); N. Tateiwa et al. , J. Phys. Condens. Matter 13 , L17 (2001). For URhGe, see D. Aoki et al. , Nature (London) 413 , 613 (2001). For UCoGe, see N. T. Huy et al. , Phys. Rev. Lett. 99 , 067006 (2007); E. Slooten et al. Phys. Rev. Lett. 103 , 097003 (2009).
- 30. K. Machida and T. Ohmi, Phys. Rev. Lett. 86 , 850 (2001); M. B. Walker and K. V. Samokhin, Phys. Rev. Lett. 88 , 207001 (2002).
- 31. F. Hardy and A. D. Huxley, Phys. Rev. Lett. 94 , 247006 (2005); Harada et al. , Phys. Rev. B 75 , 140502 (2007).
- 32. A. Shimizu, H. Ozawa, I. Ichinose, and T. Matsui, Phys. Rev. B 85 , 144524 (2012). Here, we adopted the alternative notation explained in Ref. 23 .
- 33. See, e.g., J. Cardy, ' Scaling and Renormalization in Statistical Physics ', (Cambridge University Press, 1996).
- 34. T. Noguchi, I. Ichinose and T. Matsui in preparation.
- 35. N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. M. Teller, and E. Teller, J. Chem. Phys. 21 , 1087 (1953); J. M. Thijssen, ' Computational Physics ', (Cambridge University Press, 1999).
Ikuo Ichinose Tetsuo Matsui
Lattice Gauge Theory for Condensed Matter Physics: Ferromagnetic Superconductivity as its Example
Recent theoretical studies of various strongly-correlated systems in
condensed matter physics reveal that the lattice gauge theory(LGT) developed in
high-energy physics is quite a useful tool to understand physics of these
systems. Knowledges of LGT are to become a necessary item even for condensed
matter physicists. In the first part of this paper, we present a concise review
of LGT for the reader who wants to understand its basics for the first time.
For illustration, we choose the abelian Higgs model, a typical and quite useful
LGT, which is the lattice verison of the Ginzburg-Landau model interacting with
a U(1) gauge field (vector potential). In the second part, we present an
account of the recent progress in the study of ferromagnetic superconductivity
(SC) as an example of application of LGT to topics in condensed matter physics,
. As the ferromagnetism (FM) and SC are competing orders with each other, large
fluctuations are expected to take place and therefore nonperturbative methods
are required for theoretical investigation. After we introduce a LGT describing
the FMSC, we study its phase diagram and topological excitations (vortices of
Cooper pairs) by Monte-Carlo simulations.