Loading [MathJax]/jax/output/HTML-CSS/jax.js

Feb 22, 2013

Gauge Parallel Transport

the Covariant Derivative in Differential Geometry

Let V(x) be a field of tangent vector (or what mathematicians called section) in some smooth manifold M. In mathematical jargon, it is a section of the tangent vector bundle of M. It may be expressed by local basis V(x)=VμμVμ(x)eμ. In general, the basis vectors eμ depend on the coordinates x (See the 2-sphere S2 for example).

Fig. 1 A manifold S2

So when we try to compute the derivative of V, we have to differentiate the basis as well: μV=μVνeν+Vνμeν. But the problem is, the differentiation of the basis vectors may not be tangent vectors on the manifold (See Fig. 1). Seen from the embedded space (R3 in the case of S2), the differentiation of the basis vectors has components perpendicular to the manifold. So in general μeν=Γαμνeα+baνna, where naMμV=(μVν+ΓνμαVα)eν+Vνbaνna. Thus, μV is not a tangent vector on M.

In order to introduce a proper tangent vector derivative, we can take the tangential part of  μV and call it the covariant derivative (meaning it is a tangent vector) DμV(μVν+ΓνμαVα)eν. It can be shown that the covariant derivative satisfies the rules for a derivative:

  • (linearity): D(V+W)=DV+DW;
  • (Leibniz rule): D(fV)=dfV+fDV,
where f(x) is a scalar function. In fact it is more elegant to define the covariant derivative from these properties. Nevertheless, our description from the embedded space provided an intuitive picture.

the Parallel Transport in Differential Geometry

As we have seen in last section, the partial derivative of a vector may not be a vector. Differentiation is essentially vector addition (subtraction). The closeness of the vector addition, for vectors at different coordinates, is violated because on the manifold, the vector space defined at two different points are different. This idea leads to another solution of the problem. If we can move (transport) the second vector "parallelly" to the position of the first vector and then do differentiation, the result will be a tangent vector.

However, on the manifold, we can only guarantee local parallelism. As such, a finite parallel transport may depend on the path. We denote a parallel transport along curve γ:IM from γ(s) to γ(t) as T[γts]. After the transportation, vector at γ(s), V(γ(s))T[γts]V(γ(s)) becomes a vector at γ(t). Or, in terms of the components, Vμ(γ(s))Tμν[γts]Vν(γ(s)).

Fig. 2 shows a vector (yellow) parallelly transported along a curve (black) and resulted a vector (blue) differently from the original vector.

Fig. 2 parallel transport on S2
Given such a parallel transport, we can define a covariant derivative along a curve γ as DsV(γ(s))|sslimϵ0V(γ(s+ϵ))T[γs+ϵs]V(γ(s))ϵ. For an infinitesimal interval, the covariant derivative only depends on the tangent vector of the curve at s, X˙γ(s). This gives us the directional covariant derivative: XDV(γ(s))DXV(γ(s))DsV[γ(s)]|ss=limϵ0V(γ(s)+ϵX)T[γs+ϵs]V(γ(s))ϵ It's components are denoted as (DμV)ν(DμV)νVν;μ. Compare to our finding in last section, Vν;μ=μVν+ΓνμαVα. Expand the parallel transport around γ(s), T[γs+ϵs]=1+ϵddtT[γs+ts]. Therefore, XμΓνμαVα=ddt(T[γs+ts])ναVα. To put it in a nicer form, let's define vμ(x)=Tμα[γts]Vα(x0), x=γ(t),x0=γ(s), then ddtvμ(x)=˙xρΓμρν(x)vα(x).
Similarly, the parallel transport along curve γ satisfies ddtTμν[γts]+˙γρΓμρσTσν[γts]=0.
The operator equations have Dyson-series solution: vμ(x)=Pexp{γdxρΓμρν}vν(x0)Tμν[γts]=Pexp{γtsdt˙γρΓμρν}. Because v(x) and v(x0) are two vectors defined at different point x and x, the parallel transport under the coordinate transformation should transforms as Tμν[γts]Λμα(γ(t))Λβν(γ(s))Tαβ[γ] If γ is a closed path, Tμν transforms like a tensor.

We can construct an invariant, the holonomy (physicists call it Wilson loop after Kenneth G. Wilson, who first studied the similar object in gauge theory): WCTμμ[C]=trPexp{CdxρΓρ} where C is a closed curve. In the flat space, Tμν=δμν, and Wd.

the SU(N) Gauge Symmetry

We have a similar problem in the gauge theories, where the tangent-vector fields are replaced by  quantized color-vector field. We don't have the geometric visualization, yet under a gauge transformation, a color-vector field still transforms covariantly just like the tangent vector. Let ϕi(x),i=1,2,N be an SU(N) color-vector field. Again, the partial derivative of ϕi(x) ξααϕi(x)=limϵ0ϕi(x+ϵξ)ϕi(x)ϵ, is not covariant. Indeed, let V(x) be a gauged SU(N) transformation, i.e. ϕi(x)Vij(x)ϕj(x). Then the derivative becomes ξααϕi(x)ξαα(Vij(x)ϕj(x))=Vij(x)ξααϕj(x)+ξααVij(x)ϕj(x). Similarly, the covariant derivative can be defined as, ξαDαϕi(x)=limϵ0ϕi(x+ϵξ)Uij(x+ϵξ,x)ϕj(x)ϵ,  where Uij(y,x) is a parallel transport (also named comparator, gauge link, Wilson line) in the "color" space. It satisfies the following properties,
  • It's path dependent: U(y,x)=Uγ(y,x)U[γ], where γ:[0,1]M is a curve with γ(0)=x, γ(1)=y;
  • Uγ(x,x)=1;
  • Uγ(z,y)Uσ(y,x)=Uγσ(z,x), if γ(0)=σ(1);
  • Uγ[γ1]=U1[γ];
  • U(y,x) transforms as U(y,x)V(y)U(y,x)V(x); In this way, the covariant derivative transforms covariantly: ξμDμϕ(x)limϵ0V(x+ϵξ)ϕ(x+ϵξ)V(x+ϵξ)U(x+ϵξ,x)V(x)V(x)ϕ(x)ϵ=limϵ0V(x+ϵξ)ϕ(x+ϵξ)U(x+ϵξ,x)ϕ(x)ϵ=V(x)ξμDμϕ(x)
Let γ be some curve in M. Define a curve γst along γ,
γst(u)=γ(t+u(st)),t,s[0,1].
Then U[γst]=Uγst(γ(s),γ(t)), U[γst]U[γtp]=U[γsp]. Suppose the covariant derivative is, Dμ=μigAμU(x+ϵξ,x)=1+igϵξμAμ(x)+O(ϵ2). If an operator ϕ is parallelly transported (start from x) along γ, then ϕ(γ(s))=U[γst]ϕ(γ(t)). This transport induces a field along γ. Its covariant derivative vanishes, ˙γμ(s)Dμϕ(γ(s))=ϕ(γ(s+ds))U[γs+dss]ϕ(γ(s))ds=0 Here ˙γ(s) is the tangent vector at γ(s). So ˙γμμf(x)=ddsf(γ(s)). Therefore, we have defined an initial value problem, ddsϕ(γ(s))=igdγμdsAμ(γ(s))ϕ(γ(s));(1) where Aμ(x) is an operator. s,t[0,1],st,[A(γ(s)),A(γ(t))]0. Recall initial value problem of Schroedinger's equation, ddt|ψ(t)=iH|ψ(t); The solution of this problem is |ψ(t)=Texp{it0dτH(τ)}|ψ(0),where T is the time-ordering operator.
Similarly, The solution of Eq. (1) is ϕ(γ(1))=Pexp{igγds˙γμ(s)Aμ(γ(s))}ϕ(γ(0)) where P is the path-ordering operator. P{A(γ(s1))A(γ(s2))}=θ(s1s2)A(γ(s1))A(γ(s2))+(1)Aθ(s2s1)A(γ(s2))A(γ(s1)).
Apparently, U[γ]=Pexp{igγds˙γμ(s)Aμ(γ(s))}Pexp{igγdxμAμ(x)} If γ is a closed loop, i.e. γ(0)=γ(1), trU[γ])tr{V(γ(0))U[γ]V(γ(0))}=trU[γ], is a gauge invariant. This observable is called a Wilson loop. W[γ]trPexp{igγdxμAμ(x)} In abelian gauge theories, Stokes theorem implies, W[S]=exp{igSdxμdxν(μAννAμ)}. Apparently, μAννAμ=Fμν is the field tensor. This result can be generalized to non-abelian case and again yields field tensor, but the corresponding Stokes theorem is more complicated (See [1]).


Update:
Wilson loops are fundamental gauge invariants. In quantum field theory, the vacuum expectation value (VEV) of a Wilson loop is essential to study the property of gauge fields. W[γ1]W[γ2]W[γn]=DAW[γ1]W[γ2]W[γn]exp{iS[A]} One particular intriguing theory is the Chern-Simons theory (more broadly topological quantum field theory). Edward Witten (Witten The Magnificent) proved that in Chern-Simons theory, only the topology of the Wilson loops is important. Such a observable is called a knot invariant. The study of Wilson loops becomes the study of knot invariant. In fact, the EVE of Wilson loops satisfies the skein relation. This makes topological QFT fun.

Update: Aharonov-Bohm effect
Consider a quantum mechanical particle with minimal coupling to electromagnetism,S[γ]=τfτidτ(p2+m2+epμAμ), where p=dγ/dτ˙γ is the particle momentum.
The Feynman propagator is defined, K(xf;xi)=Dγexp{iS[γ]} with γ(τi)=xi,γ(τf)=xf.
Apparently, it can be written as, K(xf;xi)=Dγexp{ieτ1τ0dxμAμ}exp{iS0[γ]}
In a double-slit experiment, the (stationary) paths reduce to two, each acquiring a phase factor due to the non-vanishing magnetic vector potential outside the solenoid. exp(ieφ1+ieφ2)=exp(ieΔφ)=exp(ieBS) equals a (the) Wilson loop enclosing the solenoid. Therefore, the amplitude (propagator) is shifted by a phase factor leading to a shift in the interference pattern, even though there is no magnetic field outside the solenoid. This is the famous Aharonov-Bohm effect.
Aharonov-Bohm effect is a result of presence of gauge link instead of tensor field in the action.

Update: Gauge Invariants
Our introduction of Wilson lines (Wilson loops as well) does not depend on specific dynamics (given by an action). In reality, a dynamical theory is important for discussion of quantization, hence for quantum Wilson lines. The action of the theory has to be gauge invariant. We have already had gauge covariant quantities, the "colored" field ϕ(x), and the covariant derivative Dμ. But we want gauge invariants/covariants involving only the gauge field to close the theory. [Dμ,Dν]=igFμν is such a covariant quantity. To get gauge invariants, one simply takes 12trFμνFμν. This term is called Yang-Mills.

Wilson loop is another gauge invariants involving pure gauge field. It can also be used to construct an action. But, just as we have shown above with Stokes theorem, Wilson loops is equivalent to Yang-Mills plus other non-linear terms of it. Non-linear terms (Yang-Mills is non-linear in non-abelian gauge symmetries) are plausible in non-linear optics for example. In lattice gauge theory, Wilson loop is actually used as the action.

˜Fρλ12ϵμνρλFμν is also gauge covariant. But tr˜Fμν˜Fμν=trFμνFμν. tr˜FμνFμν=12ϵμνρλFμνFρλ is a gauge invariant. Such a term is called θ-term in QCD.

In (2+1)-dimension (not limited to Minkowski space), we also have, k8πMd3xϵαβγtr{Aα(βAγγAβ)+23Aα[Aβ,Aγ]} This term is called Chern-Simons. k is an integer for quantized theory. Chern-Simons is an important theory, because it is metric free. Such a theory is called a topological quantum field theory. Wilson loops, also metric free, are the primary observables in Chern-Simons theory.





[1]: N. E. Bralic, Phys. Rev. D 22 (1980) 3090

Some notes on non-abelian gauge theory

Non-Abelian Gauge Theory Lecture Notes, Based on Srednicki §24 & §69 

This is the lecture notes I used for this week's lecture of PHYS 682. We ended up to Yang-Mills Lagrangian, LYM=14FμνaFaμν+ˉψi(iγμDμijmδij)ψj where Dμij=μδijigAμaτaij is the gauge covariant derivative; Fμνa=μAνaνAμa+gfabcAμbAνc is the field tensor.

Feb 14, 2013

The Secret Formula for Love


Mathematica code:

(* a xkcd-style graph converter from http://mathematica.stackexchange.com/questions/11350/xkcd-style-graphs 
by Simon Woods *)

xkcdStyle = {FontFamily -> "Comic Sans MS", 16};

xkcdLabel[{str_, {x1_, y1_}, {xo_, yo_}}] := 
  Module[{x2, y2}, x2 = x1 + xo; y2 = y1 + yo;
   {Inset[
     Style[str, xkcdStyle], {x2, y2}, {1.2 Sign[x1 - x2], 
      Sign[y1 - y2] Boole[x1 == x2]}], Thick, 
    BezierCurve[{{0.9 x1 + 0.1 x2, 0.9 y1 + 0.1 y2}, {x1, y2}, {x2, 
       y2}}]}];

xkcdRules = {EdgeForm[ef : Except[None]] :> 
    EdgeForm[Flatten@{ef, Thick, Black}], 
   Style[x_, st_] :> Style[x, xkcdStyle], 
   Pane[s_String] :> 
    Pane[Style[s, xkcdStyle]], {h_Hue, l_Line} :> {Thickness[0.02], 
     White, l, Thick, h, l}, 
   Grid[{{g_Graphics, s_String}}] :> Grid[{{g, Style[s, xkcdStyle]}}],
    Rule[PlotLabel, lab_] :> Rule[PlotLabel, Style[lab, xkcdStyle]]};

xkcdShow[p_] := 
 Show[p, AxesStyle -> Thick, LabelStyle -> xkcdStyle] /. xkcdRules

xkcdShow[Labeled[p_, rest__]] := 
 Labeled[Show[p, AxesStyle -> Thick, LabelStyle -> xkcdStyle], 
   rest] /. xkcdRules

(* Mathematica Version 8.0+ *)
xkcdDistort[p_] := 
  Module[{r, ix, iy, rand}, 
   r = ImagePad[Rasterize@p, 10, Padding -> White];
   {ix, iy} = 
    Table[RandomImage[{-1, 1}, ImageDimensions@r]~ImageConvolve~
      GaussianMatrix[10], {2}];
   ImagePad[
    ImageTransformation[
     r, # + 15 {ImageValue[ix, #], ImageValue[iy, #]} &, 
     DataRange -> Full], -5]];

(* for Mathematica version < 7 *)
xkcdDistort[p_] := 
 Module[{r, id, ix, iy, samplepoints, funcs, channels},
  r = ImagePad[Rasterize@p, 10, Padding -> White]; 
  id = Reverse@ImageDimensions[r];
  {ix, iy} = Table[ListInterpolation[ImageData[
      Image@RandomReal[{-1, 1}, id]~ImageConvolve~GaussianMatrix[10]]], {2}]; 
  samplepoints = Table[{x + 15 ix[x, y], y + 15 iy[x, y]}, {x, id[[1]]}, {y, id[[2]]}]; 
  funcs = ListInterpolation[ImageData@#] & /@ ColorSeparate[r]; 
  channels = Apply[#, samplepoints, {2}] & /@ funcs; 
  ImagePad[ColorCombine[Image /@ channels], -10]]


xkcdConvert[x_] := xkcdDistort[xkcdShow[x]]


(* the heart plot *)

xkcdConvert[
 ContourPlot[x^2 + (y - (x^2)^(1/3))^2 == 1, {x, -1, 1}, {y, -1, 2}, 
  ContourStyle -> {Red, Thick}, ContourShading -> Red, 
  AspectRatio -> 1, ImageSize -> 400, Axes -> True, Frame -> False, 
  AxesLabel -> Map[Text@Style[#, 20, Italic] &, {"x", "y"}], 
  AxesStyle -> FontSize -> 15, 
  Epilog -> {Inset[
     Style[x^2 + (y - (x^2)^(1/3))^2 == 1, Medium, Italic], 
     Offset[{0, 0}], Scaled[{-0.1, -5.5}], 1], 
    Inset[Style["xkcd-style graph", xkcdStyle, FontSize -> 12 ], 
     Scaled[{0.15, 0.05}], Automatic]}, 
  PlotLabel -> 
   Style["The Secret Formula for Love", xkcdStyle, 
    FontSize -> 24]]]

Feb 4, 2013

Cisco AnyConnect Client was not able to establish connection ...


I installed and used Cisco AnyConnect for a while and then suddenly one day, I cannot make a connection. The prompted error information is as following,


The Cisco troubleshooting website just told me to retry, not helpful at all. After some googling, it seems the trouble is caused by incompatibility of AnyConnect client and Windows network and file sharing service. But I have my network and file sharing option disabled already. This post suggests, not only the network and file sharing option in network and sharing center has to be disabled, but also the Internet Connection Sharing (ICS) service has to be shut down. It can be done via the following steps: find
start -> services.msc -> Internet Connection Sharing 
and disable it. (The following prompt windows is in Chinese, but similar to ones of other languages)