1 / 29

Topology conserving gauge action and the overlap Dirac operator

Topology conserving gauge action and the overlap Dirac operator. Hidenori Fukaya Yukawa Institute, Kyoto Univ. Collaboration with S.Hashimoto (KEK,Sokendai), T.Hirohashi (Kyoto Univ.), H.Matsufuru(KEK), K.Ogawa(Sokendai), T.Onogi(YITP). Contents. Introduction Lattice simulations

elan
Download Presentation

Topology conserving gauge action and the overlap Dirac operator

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Topology conserving gauge actionand the overlap Dirac operator Hidenori Fukaya Yukawa Institute, Kyoto Univ. Collaboration with S.Hashimoto (KEK,Sokendai), T.Hirohashi (Kyoto Univ.), H.Matsufuru(KEK), K.Ogawa(Sokendai), T.Onogi(YITP)

  2. Contents • Introduction • Lattice simulations • The static quark potential • Stability of the topological charge • The overlap Dirac operator • Conclusion and outlook

  3. 1. Introduction Lattice regularization of the gauge theory is a very powerful tool to analyze strong coupling regime but it spoils a lot of symmetries… • Translational symmetry • Lorentz invariance • Chiral symmetry or topology • Supersymmetry…

  4. Nucl.Phys.B185,20 (‘81),Nucl.Phys.B193,173 (‘81) • Nielsen-Ninomiya theorem Any local Dirac operator satisfying chiral symmetry; has unphysical poles (doublers). Example – free naïve fermion – • Continuum has no doubler. • Lattice has unphysical poles at . • Wilson Dirac operator Doublers are decoupled (heavy) but no chiral symmetry.

  5. Phys.Lett.B417,141(‘98) • Ginsparg-Wilson relation The Neuberger’s overlap operator: satisfying the Ginsparg-Wilson relation: realizes ‘modified’ exact chiral symmetry on the lattice; the action is invariant under NOTE • Expansion in Wilson Dirac operator ⇒No doubler. • Fermion measure is not invariant; ⇒chiral anomaly, index theorem locality? or smoothness? Phys.Rev.D25,2649(‘82)

  6. The locality of the overlap operator If the gauge fields are smooth; with a (small) fixed number ε , it gives a lower bound to ⇒ (exponential) locality and smoothness of the Dirac operator; ⇒ Stability of the index (topological charge); P.Hernandez,K.Jansen,M.Luescher,Nucl.Phys.B552,363(‘99)

  7. M.Luescher,Phys.Lett.B428,342(‘98),Nucl.Phys.B549,295(‘99)M.Luescher,Phys.Lett.B428,342(‘98),Nucl.Phys.B549,295(‘99) • Topology conserving gauge action The bound is automatically satisfied with the action • Hybrid Monte Carlo (HMC) algorithm • Heat bath algorithm = local and large updates • HMC algorithm = global and small updates ⇒Q may be fixed along the simulations

  8. - Motivations - • We studied the Nf=2 massive Schwinger model • with a theta term using the topology conserving action and the domain-wall fermion action; • HF,T.Onogi,Phys.Rev.D68,074503(‘03),Phys.Rev.D70,054508(’04) • Good chiral behavior (consistent with the bosonization approach). • Perfect stability of Q with ε= 1.0. • Efficient sampling of the configurations with fixed Q. ⇒ ε- regime with fixed Q • Improvement of the locality of the overlap operator. • Lower numerical costs of the overlap operator ⇒dynamical overlap fermions • Anyway, it’s interesting…

  9. - Our goals - • Scaling studies (to show validity in 4D quenched QCD.) • determination of the lattice spacing • scaling violations of the quark potential • How stable Q ? • Locality and numerical cost of GW fermion improved ? c.f. S.Shcheredin, W.Bietenholz, K.Jansen, K.I.Nagai, S.Necco and L.Scorzato, hep-lat/0409073 , hep-lat/0412017 .

  10. 2. Lattice simulations • Action: Luescher action (quenched) with 1/ε= 1.0, 2/3, 0.0 (=plaquette action) . • Algorithm: The standard HMC method. • Lattice size : 124,164,204 . • 1 trajectory = 20 - 40 molecular dynamics steps with stepsize Δτ= 0.01 - 0.02. The simulations were done on the Alpha work station at YITP and SX-5 at RCNP.

  11. The overlap Dirac operator • We use the implicit restarted Arnoldi method (ARPACK) to calculate the eigenvalues of . • To compute , we use the Chebyshev polynomial approximation after subtracting 10 lowest eigenmodes exactly. • We set s=0.6 ARPACK, available from http://www.caam.rice.edu/software/

  12. Initial configuration For topologically non-trivial initial configuration, we use a discretized version of instanton solution on 4D torus; which gives constant field strength with arbitrary Q. A.Gonzalez-Arroyo,hep-th/9807108, M.Hamanaka,H.Kajiura,Phys.Lett.B551,360(‘03)

  13. New cooling method to measure Q We “cool” the configuration smoothly by performing HMC steps with exponentially increasing (The bound is always satisfied along the cooling). ⇒ We obtain a “cooled ” configuration close to the classical background at very high β～106, (after 40-50 steps) then gives a number close to the index of the overlap operator. NOTE: 1/εcool= 2/3 is useful for 1/ε= 0.0 .

  14. Eigen function （chirality +） Eigen function （chirality －） Eigen values of γ５D Cooling -Example: 2D QED- - instanton anti-instanton pair on 162 lattice - Let’s cool it !

  15. →1 instanton Eigenvalues Cooling -Example: 2D QED- - 2 instanton and 1 anti-instanton - Let’s cool it !

  16. Cooling -4D QCD results- The agreement of Q with cooling and the index of overlap D is roughly (with only 20-80 samples) • ~ 90-95% for 1/ε= 1.0 and 2/3. • ~ 60-70% for 1/ε=0.0 (plaquette action)

  17. 3. The static quark potential In the following, we assume Q does not affect the Wilson loops. ( initial Q=0 ) • We measure the Wilson loops, in 6 different spatial direction, using smearing. G.S.Bali,K.Schilling,Phys.Rev.D47,661(‘93) • The potential is extracted as . • From results, we calculate the force following S.Necco,R.Sommer,Nucl.Phys.B622,328(‘02) • Sommer scales are determined by

  18. Assuming r0 ~ 0.5 fm, we obtain (preliminary)

  19. The quark potential itself is good.

  20. Perturbative renormalization of the coupling At 1-loop level, renormalized couplings have similar values but the convergence of the perturbative series is much worse for 1/ε= 1, 2/3. Tadpole improvement at 2-loop level may be required or non-perturbative approach should be done. R.K.Ellis,G.Martinelli, Nucl.Phys.B235,93(‘84) Erratum-ibid.B249,750(‘85)

  21. SG ε= 1.0 Q=0 Q=1 If the barrier is high enough, Q may be fixed. SG ε< 1/30 Q=0 ε=∞Q=1 4. Stability of the topological charge The stability of Q for 4D QCD is proved only when ε＜ εmax～1/30 ,which is not practical… Let us search the parameters; β,ε,V which can fix Q.

  22. M.Luescher, hep-lat/0409106 Appendix E. We measure Q using cooling per 20 trajectories : auto correlation for the plaquette : total number of trajectories : (lower bound of ) number of topology changes We define “stability” by the ratio of topology change rate ( ) over the plaquette autocorrelation( ). Note that this gives only the upper bound of the stability, since there may be topology changes during every 20 trajectories.

  23. Stability is better for smaller ε, smaller L , smaller a. • A significantly good stability can be seen • for • ⇒maybe valid in ε- regime.

  24. Q dependence of the quark potential seems week as we expected.

  25. 5. The overlap Dirac operator It is known that the lowest eigen value of is related to • Localityis guaranteed if never touches zero. • The order of Chebyshev polynomial can be reduced; : the condition number : a random vector

  26. The locality For should exponentially decay. 1/a~0.08fm (with 4 samples), no remarkable improvement of locality is seen… ⇒lower beta? ＋ : beta = 1.42, 1/e=1.0 × : beta = 2.7, 1/e=2/3 ＊ : beta = 6.13, 1/e=0.0

  27. and are nor independent. • For 1/ε=1.0, is 1.6~2.5 times larger ( is 1.2~1.4 times larger) than that with the plaquette action. We use 10 different Npol and 4 different confs for each parameter set. • The condition number and Chebyshev approximation

  28. 6. Conclusion and outlook Topology conserving gauge action may be practical. • New cooling method does work. • The lattice spacing can be determined in a conventional manner, ant the quark potential show no large deviation from the continuum limit. • Q can be stable for . • No improvement of the locality (for high beta). • The numerical cost of Chebyshev approximation would be 1.2-2.5 times better than that with plaquette action.

  29. Outlooks • Accept/reject for topology changes to fix Q completely (per every 100trj would be enough.). • Quenched ε-regime with fixed Q. • Nf=2 full QCD in the ε-regime with the overlap fermion.

More Related