Topological Susceptibility under Gradient Flow
Abstract
We study the impact of the Gradient Flow on the topology in various models of lattice field theory. The topological susceptibility is measured directly, and by the slab method, which is based on the topological content of subvolumes (“slabs”) and estimates even when the system remains trapped in a fixed topological sector. The results obtained by both methods are essentially consistent, but the impact of the Gradient Flow on the characteristic quantity of the slab method seems to be different in 2flavour QCD and in the 2d O(3) model. In the latter model, we further address the question whether or not the Gradient Flow leads to a finite continuum limit of the topological susceptibility (rescaled by the correlation length squared, ). This ongoing study is based on direct measurements of in lattices, at .
Topological Susceptibility under Gradient Flow
Héctor MejíaDíaz^{}^{}Speaker, email: he_mejia@yahoo.com.mx, Wolfgang Bietenholz, Krzysztof Cichy, Philippe de Forcrand,
Arthur Dromard, Urs Gerber, and Ilya Orson Sandoval
1Instituto de Ciencias Nucleares,
Universidad Nacional Autónoma de México
A.P. 70543, C.P. 04510 Ciudad de México, Mexico
\@textsuperscript2GoetheUniversität Frankfurt am Main, Institut für Theoretische
Physik
MaxvonLaueStraße 1, 60438 Frankfurt am Main, Germany
\@textsuperscript3Faculty of Physics, Adam Mickiewicz University,
Umultowska 85, 61614 Poznán, Poland
\@textsuperscript4Institut für Theoretische Physik, ETH Zürich,
WolfgangPauliStrasse 27, CH–8093 Zürich, Switzerland
\@textsuperscript5CERN, Theory Division, CH1211 Genève 23, Switzerland
\@textsuperscript6Institut für Theoretische Physik, Universität
Regensburg, D93040 Regensburg, Germany
\@textsuperscript7Instituto de Física y Matemáticas,
Universidad Michoacana de San Nicolás de Hidalgo
Edificio C3, Apdo. Postal 282, C.P. 58040,
Morelia, Michoacán, Mexico

Abstract.
1 Introduction
In some quantum field theories, the set of configurations is divided into topological sectors, labelled by a topological charge . This is the case in QCD, and in dimensional O() models (with periodic boundary conditions for the gluon and spin fields), due to and . Hence this class of models includes 2flavour QCD, as well as the 1d O(2) and the 2d O(3) model, which we are going to deal with.
For usual lattice actions, all configurations can be continuously deformed into one another, at finite action, hence there are no topological sectors in a strict sense. Exceptions are topological lattice actions, with a sharp cutoff for the angles between nearest neighbour spin variables [1], or for each plaquette variable [2] (see also Ref. [3]) in spin models and gauge theories, respectively. However, even for conventional lattice actions there are established ways to divide the configurations into sectors, which turn into topological sectors in the continuum limit.
Here we consider 2flavour QCD with twistedmass quarks [4] (at full twist) and the Wilson gauge action. For the O() models we employ the standard lattice action,
(1) 
where the sum runs over all nearest neighbour lattice sites. For these O() models we apply the geometric definition of the topological charge density on the lattice [5], which leads to integer charges . In QCD we use a clover discretisation of , where is the field strength tensor.
In all cases under consideration, parity symmetry implies , hence the topological susceptibility takes the form
(2) 
2 The slab method to measure the topological susceptibility
Once we have fixed a formulation of the topological charge on the lattice, it is straightforward to measure by means of Monte Carlo simulations, if the Markov chain frequently changes , such that the sectors are sampled correctly. In practice, however, such simulations are often confronted with the severe problem of “topological freezing”: in particular, the algorithms, which proceed in small update steps, tend to get stuck in one topological sector for a huge number of steps, since the topological sectors are effectively separated by high potential barriers. The autocorrelation time with respect to increases with a high power of the inverse lattice spacing as we approach the continuum limit (“topological slowing down”), see e.g. Ref. [6].
A variety of approaches to handle this problem is reviewed in Ref. [7]. One strategy aims at extracting physical observables even from a Markov chain which is entirely trapped in a single topological sector. For general observables such a method was suggested in Ref. [8], and tested and extended in Refs. [9, 10, 11, 12]. More specifically, a procedure to measure within a fixed topological sector was proposed in Ref. [13] and tested in Refs. [14, 9, 15, 12]. Here we consider the slab method as another way to evaluate from data obtained at fixed (actually data from can be combined). The idea was mentioned in Ref. [16], implemented in Ref. [17], and further explored in Refs. [18, 12, 19]. A different variant was applied in Ref. [20], and there are similarities with the approach in Ref. [21].
We briefly review the simplest version of the slab method, which assumes the statistical distribution of the topological charges to be Gaussian [17], . We split the volume into two subvolumes (“slabs”) and (). By summing up the topological charge density in each of them, in a configuration of topological charge , we obtain the slab charges (they do not need to be integer, since the slabs do not have periodic boundaries). At fixed , and , the corresponding slab probability distributions and obey
(3) 
Measuring yields a value for . A sequence of such measurements, at different parameters , enables a fit to the prediction
(4) 
which provides a result for . In practice, the most reliable fitting regime is around the center, , because the size of the slabs should be large compared to that of topological excitations; one should not include or (these regions involve very small slabs, where the Gaussian distribution is not a good approximation).
2.1 Results by the slab method under Gradient Flow (GF)
As a smoothing procedure for lattice field configurations, the GF corresponds to a renormalisation group scheme. When the GF proceeds, the distinction between topological sectors becomes more marked, approaching the continuum feature of a total separation [22, 23].
The slab method has been tested in 2flavour QCD, with twistedmass quarks and the Wilson gauge action, in a volume , at and bare mass , which corresponds to a pion mass of and a lattice spacing [18, 19].
The GF was implemented with the RungeKutta method; the results with time step and agree. The GF flow time unit was fixed to , based on the criterion proposed in Refs. [22, 23]: , where is the mean energy density. The effect of the GF on the curves to be fitted in the slab method is shown in Figure 1 (left). As the GF proceeds, the fit has to be restricted to a narrower interval centered at . Moreover, the fitting function (4) has to be extended to
(5) 
where is a constant (with respect to ), which increases
roughly like [19].
The fits at different instances of the GF time,
, yield very stable results
for , which are compatible both with a direct measurement
(,
and with the method of Ref. [13]
(, see Figure 1
(right). In the interval we consistently
obtain . This is a
success of the slab method, but the rôle of the subtractive
constant is not obvious.
In the continuum, the Gradient Flow in O() models takes the form [24]
(6) 
where is the GF time (of dimension [length]) and is the Laplace operator, which we handle by standard lattice discretisation. In order to proceed in discrete flow time steps, we apply the RungeKutta method. For a given configuration we first compute the gradients of all spin variables (at all times required by the RungeKutta 4point scheme), then all spins are modified simultaneously, with a time step ; afterwards the spins are normalised again.
In this manner, we considered lattices, for instance with , , where the correlation length amounts to [25]. Figure 2 (left) shows the (approximate) slab parabolae obtained for in the sector , in even multiples of the flow time unit (which obeys , cf. Section 3). We see a qualitative difference from the QCD result in Figure 1: here the fits do not require any subtractive constant, i.e. the original formula (4) can be used, and keeps on decreasing as the GF proceeds (the curvature of the parabola is reduced). This feature was also observed in all data sets to be reported in Section 3.
The fitting results for — obtained by the slab
method, separately in the sectors — are
close to the directly measured values, as we see in Figure
2 (right).
In this case, the direct measurement is not problematic, since our
simulations were carried out with the Wolff cluster algorithm,
which proceeds in nonlocal update steps, thus suppressing the
effect of topological freezing [26].
In order to investigate further these qualitatively different behaviours, we tested the 1d O(2) model, or quantum rotor, as a toy model. Here we refer to , as an example. In infinite volume we can compute analytically [27] (still for the standard lattice action, in lattice units)
(7) 
in agreement with our simulation results for size (we used again the cluster algorithm). Regarding the slab method under GF, Figure 3 shows the results for in the sector at various flow times , and as a function of . For the scaling quantity we obtain numerically at , and at ; thus we gradually approach the continuum value of .
Regarding the different features in QCD and in the 2d O(3) model, the quantum rotor at short flow times () seems compatible with the latter, but at a nonnegligible constant has to be subtracted for a successful fit; at it takes values in the sectors .
In the 2d O(3) model our flow times are short so far, cf. Section 3; the question whether the same behaviour sets in after a long GF is under investigation.
3 Topological scaling in the 2d O(3) model
In a dimensional quantum field theory with topological sectors, the dimensionless term is supposed to be a scaling quantity, which converges to a finite value in the continuum limit. For the 1d O(2) model the continuum value is attained without problems [1, 11, 27]. In QCD, and in SU() YangMills theories (), a straight approach based on the expression , where is the lattice topological charge density, faces problems. In conventional formulations one encounters a divergence due to the point , which prevents the continuum scaling. In these cases, there are known solutions to this problem, in particular the application of the GF [22, 23].
In the 2d O(3) model, this question has been controversial in the 1980s and 1990s. The consensus is now that the quantity seems to diverge in the continuum limit, as first observed in [5], and later underpinned by semiclassical studies, e.g. in Refs. [28, 29]. Again the problem can be traced back to the topological density correlation at distance zero (see e.g. Ref. [1]), and the semiclassical picture suggests an abundance of very small topological windings, socalled dislocations, as the continuum limit is approached. Ref. [30] constructed and applied a sophisticated (truncated) classically perfect lattice action, with a host of couplings beyond nearest neighbour lattice sites, which suppress such dislocations. Nevertheless, the simulation results with this action still suggest a logarithmic divergence of the term in the continuum limit ( in lattice units). Hence the continuum limit of this popular model is generally assumed to be illdefined, at least with respect to its topology (the reason is again the contribution , see e.g. Refs. [31, 32, 1]).
However, the question remains whether or not this divergence could be overcome by the GF; we gave preliminary results in Ref. [25]. In the following we summarise the status of this study. So far it involves nine lattices, in the range of , where in each volume has been tuned such that . Therefore, increasing corresponds to a controlled step towards the continuum limit, at a fixed and large physical box size. The statistics in each volume are configurations (generated by the Wolff cluster algorithm, both in the singlecluster and the multicluster version). We repeat that we perform the GF with the RungeKutta 4point method, with a time step of , which is simultaneously applied to all spin variables.
Figure 4 (left) shows the early GF time evolution of the product . At longer flow times it increases up to some maximum, before monotonically decreasing again. For increasing volume, the value of the maximum of decreases, hence the reference value has to be sufficiently small, such that it is attained in all volumes under consideration. We are in the process of extending this study up to and , where for instance is not attained anymore. Hence we chose the reference value , which works up to [33]. So we refer to the definition
(8) 
as we anticipated in Figure 2.
Figure 4 (right) shows an example for the zeromomentum spinspin correlation function at . At a fixed distance between Euclidean time layers, increases as the GF proceeds, but the correlation length remains unchanged within the errors. We measured by a fit to a coshfunction in the range . We see that the GF — up to — hardly affects this longrange property.
Figure 5 (left) illustrates the GF time evolution of the term , which is supposed to be the scaling quantity. We see a rapid decrease at an early stage of the GF flow, in particular at large and . This observation is compatible with an increasing dominance of small dislocations on fine lattices: these topological windings are destroyed even by a short GF flow.
Finally, Figure 5 (right) shows as a function of , at various multiples of the flow time unit (the lines are drawn to guide the eye). At this stage, no trend towards a stabilisation is visible. The observed behaviour at is well compatible with a logarithmically divergent function of the form ; at the quality of this fit is somewhat worse, but it still follows roughly this behaviour (although a powerlaw can be fitted with a similar quality) [25]. In any case, the present data do not provide a basis for revising the standard lore of a topologically illdefined continuum limit in this model.
4 Summary and outlook
In Section 2 we have discussed the slab method, which enables
a reliable measurement of the topological susceptibility within a
fixed topological sector, i.e. from a Markov chain at fixed
. This method does still work quite well when the GF is applied,
but in some cases — in particular in QCD — we observed the necessity
to subtract a constant from the expected fitting function.
In the 2d O(3) model, the data presented in Section 3 do not suggest a continuum convergence of the term after application of the GF. Further details, including a table with numerical results, are given in Ref. [25].
However, the impact range of the GF is short in these examples. It can be estimated based on the heat kernel ; in dimensions we obtain
(9) 
In our 2d model it attains at most at this stage of our study; this refers to our largest volume, , with .
In order to arrive at conclusive results, we are now going to fix
an extended GF time unit by a condition
, and investigate flow times up
to an impact range of .
This study is in progress [33],
and it should finally reveal whether or not the GF leads
to a continuum scaling of the quantity .
Acknowledgements We thank Martin Lüscher for attracting our interest to the subject of Section 3, and for advice regarding the strategy towards conclusive results. We further thank Marc Wagner for helpful discussions about the slab method, which we discussed in Section 2, and the organisers of the 35th International Symposium on Lattice Field Theory. This work was supported by DGAPAUNAM, grant IN107915, by the Consejo Nacional de Ciencia y Tecnología (CONACYT) through project CB2013/222812, and by the Helmholtz International Center for FAIR within the framework of the LOEWE program launched by the State of Hesse. K.C. was supported by the Deutsche Forschungsgemeinschaft (DFG), project nr. CI 236/11, and A.D. by the Emmy Noether Programme of the DFG, grant WA 3000/11.
References
 [1] W. Bietenholz, U. Gerber, M. Pepe, U.J. Wiese, JHEP 12, 020 (2010), 1009.2146
 [2] O. Akerlund, P. de Forcrand, JHEP 06, 183 (2015), 1505.02666
 [3] M. Lüscher, Nucl. Phys. B549, 295 (1999), heplat/9811032
 [4] R. Frezzotti, P.A. Grassi, S. Sint, P. Weisz (Alpha), JHEP 08, 058 (2001), heplat/0101001
 [5] B. Berg, M. Lüscher, Nucl. Phys. B190, 412 (1981)
 [6] S. Schaefer, R. Sommer, F. Virotta (ALPHA), Nucl. Phys. B845, 93 (2011), 1009.5228
 [7] M.G. Endres, PoS LATTICE2016, 014 (2016), 1612.01609
 [8] R. Brower, S. Chandrasekharan, J.W. Negele, U.J. Wiese, Phys. Lett. B560, 64 (2003), heplat/0302005
 [9] W. Bietenholz, I. Hip, S. Shcheredin, J. Volkholz, Eur. Phys. J. C72, 1938 (2012), 1109.2649
 [10] A. Dromard, M. Wagner, Phys. Rev. D90, 074505 (2014), 1404.0247
 [11] W. Bietenholz, C. Czaban, A. Dromard, U. Gerber, C.P. Hofmann, H. MejíaDíaz, M. Wagner, Phys. Rev. D93, 114516 (2016), 1603.05630
 [12] A. Dromard, Ph.D. thesis, GoetheUniversität Frankfurt am Main (2016)
 [13] S. Aoki, H. Fukaya, S. Hashimoto, T. Onogi, Phys. Rev. D76, 054508 (2007), 0707.0396
 [14] S. Aoki et al. (TWQCD, JLQCD), Phys. Lett. B665, 294 (2008), 0710.1130
 [15] I. Bautista, W. Bietenholz, A. Dromard, U. Gerber, L. Gonglach, C.P. Hofmann, H. MejíaDíaz, M. Wagner, Phys. Rev. D92, 114510 (2015), 1503.06853
 [16] P. de Forcrand, M. García Pérez, J.E. Hetrick, E. Laermann, J.F. Lagae, I.O. Stamatescu, Nucl. Phys. Proc. Suppl. 73, 578 (1999), heplat/9810033
 [17] W. Bietenholz, P. de Forcrand, U. Gerber, JHEP 12, 070 (2015), 1509.06433
 [18] A. Dromard, W. Bietenholz, K. Cichy, M. Wagner, Acta Phys. Polon. Supp. 9, 635 (2016), 1605.08637
 [19] W. Bietenholz, K. Cichy, P. de Forcrand, A. Dromard, U. Gerber, PoS LATTICE2016, 321 (2016), 1610.00685
 [20] S. Aoki, G. Cossu, H. Fukaya, S. Hashimoto, T. Kaneko (JLQCD) (2017), 1705.10906
 [21] R.C. Brower et al. (LSD), Phys. Rev. D90, 014503 (2014), 1403.2761
 [22] M. Lüscher, JHEP 08, 071 (2010), [Erratum: JHEP03,092(2014)], 1006.4518
 [23] M. Lüscher, PoS LATTICE2010, 015 (2010), 1009.5877
 [24] H. Makino, H. Suzuki, PTEP 2015, 033B08 (2015), 1410.7538
 [25] I.O. Sandoval, W. Bietenholz, P. de Forcrand, U. Gerber, H. MejíaDíaz, Topology in the 2d Heisenberg Model under Gradient Flow, in 31st Annual Meeting of the Division of Particles and Fields (DPyC) of the Mexican Physical Society (DPyCSMF2017) Mexico City, Mexico, May 2426, 2017 (2017), 1709.06180, http://inspirehep.net/record/1624423/files/arXiv:1709.06180.pdf
 [26] U. Wolff, Phys. Rev. Lett. 62, 361 (1989)
 [27] W. Bietenholz, R. Brower, S. Chandrasekharan, U.J. Wiese, Phys. Lett. B407, 283 (1997), heplat/9704015
 [28] F. Farchioni, A. Papa, Nucl. Phys. B431, 686 (1994), heplat/9407026
 [29] C. Michael, P.S. Spencer, Phys. Rev. D50, 7570 (1994), heplat/9404001
 [30] M. Blatter, R. Burkhalter, P. Hasenfratz, F. Niedermayer, Phys. Rev. D53, 923 (1996), heplat/9508028
 [31] J. Balog, M. Niedermaier, Phys. Rev. Lett. 78, 4151 (1997), hepth/9701156
 [32] J. Balog, M. Niedermaier, Nucl. Phys. B500, 421 (1997), hepth/9612039
 [33] W. Bietenholz, P. de Forcrand, U. Gerber, H. MejíaDíaz, I.O. Sandoval, in preparation (2017)