Propagators and topology

Two popular perspectives on the non-perturbative domain of Yang–Mills theories are either in terms of the gluons themselves or in terms of collective gluonic excitations, i.e. topological excitations. If both views are correct, then they are only two different representations of the same underlying physics. One possibility to investigate this connection is by the determination of gluon correlation functions in topological background fields, as created by the smearing of lattice configurations. This is performed here for the minimal Landau gauge gluon propagator, ghost propagator, and running coupling, both in momentum and position space for SU(2) Yang–Mills theory. The results show that the salient low-momentum features of the propagators are qualitatively retained under smearing at sufficiently small momenta, in agreement with an equivalence of both perspectives. However, the mid-momentum behavior is significantly affected. These results are also relevant for the construction of truncations in functional methods, as they provide hints on necessary properties to be retained in truncations.


Introduction
The non-perturbative, low-energy domain of Yang-Mills theory, and in extension of QCD, remains an interesting conceptual challenge, even if approaches like, e.g., lattice gauge theory, functional methods, or chiral perturbation theory already permit to determine many quantities of practical interest. One of the central questions has been for a very long time what the effective degrees of freedom at low energies are. Two perspectives are in terms of the gauge bosons themselves [1][2][3][4][5][6] and in terms of collective, i.e. topological, excitations [7][8][9][10][11][12], like (center) vortices [7], monopoles [7,8,10], instantons [13,14], calorons [15,16], merons [17], and dyons [18,19]. These topological configurations are likely not all independent but intricately related [7,[20][21][22][23]. a e-mail: axelmaas@web.de Provided both views are correct, they are necessarily just two different representations of the same physics. Since for both views a plethora of evidence exists, this appears likely to be the case. In fact, for simpler models such relations are explicitly known [2,24]. Thus, it should be possible to establish this relation explicitly, a challenge which remains so far unsolved [25]. However, since both, gluons and gluonic excitations, are inherently gauge dependent, any such relation could be itself gauge dependent. Still, it would be significant progress to establish the details of this connection at least for one gauge. There have been a number of investigations contributing to this endeavor, in both Coulomb and Landau gauge [1,[26][27][28][29][30][31][32][33][34][35], 1 with respect to different types of topological excitations. These investigations provided evidence that such a link indeed exists, and that the most characteristic low-momentum features of gluonic correlation functions are likely reflecting (or formed by, depending on perspective) topological excitations.
Lattice calculations, using smearing [39][40][41], provide a tool to isolate from field configurations the (self-dual) topological part. This permits another way to determine correlation functions of gluons in a topological field configuration, i.e. how the gluons inside a collective gluonic excitation behave. This possibility has so far only been explored in preliminary investigations [1,32]. Here, this will be extended to a full systematic investigation of this possibility in (minimal) Landau gauge for the gluon, the corresponding ghost, and the running coupling, both in momentum and position space. Because these investigations are computationally much more expensive than calculating just the propagators, this will be only possible over a limited range of lattice settings, and therefore will not yet answer questions about the deep infrared. However, they are complementary to investigations using center projections [27,30,33,34], and therefore they offer a novel perspective on the interplay of collective and single gluon excitations at low energies.
The technical details of these investigations are briefly described in Sect. 2. A discussion of the selection criteria for the background configurations on which to measure the propagators is given in Sect. 3. Results for the gluon are then presented separately for momentum space and position space in Sects. 4.1 and 4.2, respectively, in the results Sect. 4. The ghost propagator is investigated in Sect. 4.3. Section 4 also contains results on further derived quantities. The results are summarized in Sect. 6, after a speculative interpretation in Sect. 5. There are also a few remarks that will be made on how these results can be exploited to ensure that truncations of functional equations [1][2][3][4][5][6] adequately capture the contributions from collective excitations.

Technical details
The lattice configurations have been created using the standard SU(2) Yang-Mills Wilson action [39], using the hybridoverrelaxation algorithm described in [44]. The list of configurations and lattice parameters are given in Table 1.
There are a number of possibilities to isolate the topological content of the generated lattice configurations. All of these algorithms appear to provide qualitatively similar results, but differ quantitatively by 10-30 % [40,45]. Since for this first investigation the qualitative effects are most interesting, these deviations are not too important, and hence only a single method will be used. This will be APE-smearing [41] with a lexicographical update. This procedure is equivalent to the mathematically better defined Wilson flow, provided the smearing levels are suitably selected [45], a subject to be discussed in more detail in Sect. 3.
In the APE-smearing process, a single smearing sweep replaces all link variables by the prescription where "projected to the group" implies that the non-group element U μ found after addition is replaced by the group element U μ closest to the result, where the distance is given by trU μ U μ , with no summation implied. For the present SU(2) case this can be achieved by a multiplicative factor. The parameter α can be used to tune the smearing. It will be set throughout to 0.55, which is a value that is very convenient in many investigations [40]. Typically, only a few APE-smearing sweeps are used in spectroscopy applications. But to smear to an essentially topological content, i.e. one satisfying approximately the (anti-) self-duality equations F μ = ± μνρσ F ρσ , typical for topological (instanton-like) Table 1 The configurations employed. N is the (symmetric) extent of the lattice of total volume N 4 . The lattice spacing at given gauge coupling β has been determined using the data of [42], setting the string tension to (440 MeV) 2 . The number of configurations before and after the slash are the number used for the short and long cooling, respectively; see text. Therm. and sweeps are the number of configurations dropped for thermalization and between two measurements. To reduce correlations, prevalent for topological quantities [43], typical O(100) independent runs have been performed for each lattice setting. Long gives the number of total APE sweeps performed for the self-dual configurations, while Int. gives the measurement interval between APE sweeps in this case configurations [46], several more APE-smearing sweeps are necessary. Depending on the lattice parameters, this can be a few tens, or several hundred. It should be noted that the smearing is done here in all directions, not just in spatial directions, owing to the aim of creating self-dual configurations with full O(4) symmetry. There will be two possibilities considered throughout. The first will be short smearing, with 25 or less smearing sweeps. This will eliminate the ultraviolet fluctuations, but should leave most of the infrared structure intact. The second will be a long smearing, where the configurations are moved far into the self-dual regime. Since this depends on the volume, the number of smearing steps is chosen volume dependent. For the short smearing a measurement will be performed after each smearing sweep, while for the long smearing only in some interval, again depending on the volume, measurements will be performed. These intervals are given in Table 1.
The measurements of the propagators requires to fix the gauge, which is chosen to be the minimal Landau gauge [1], using the algorithm described in [44]. Since the gauge is not preserved under smearing, it was necessary to fix the gauge for each measurement anew. This made this investigation comparatively expensive, and therefore only the limited set of lattice setups, and especially the rather small volumes, listed in Table 1 could be investigated. This also implies that the lowest reachable momenta are still large, even on the largest and coarsest lattice they are O(250) MeV. Thus, the far infrared domain, an actively researched area [1], could not be investigated. However, the domain most relevant to (hadronic) bound states [2,3] is accessible.
The quantities measured are the color-diagonal part 2 of the gluon propagator D and the ghost propagator D G in momentum space and the gluon propagator also in position space, the latter also denoted as the Schwinger function . The measurement is performed with the methods described in [1,44]. Renormalization will be performed for the unsmeared propagators at μ = 2 GeV with μ 2 D(μ 2 ) = μ 2 D G (μ 2 ) = 1, though this is essentially only needed in Fig. 9. Hence, results shown are unrenormalized, except when stated otherwise.
The propagators suffice to calculate [47] the running coupling in the miniMOM scheme [48] as Since this is a renormalization-group-invariant quantity, it can be calculated using the unrenormalized propagators, after fixing α(μ 2 ) [47].
It is finally interesting to see whether the propagators depend on the topological charge. To measure the topological charge, the simplest lattice realization of the continuum 2 The color-off-diagonal part is zero in Landau gauge [1]. topological charge density operator will be used. This will be performed by calculating first the field strength tensor F μν at site x from the link variables, and then calculating the product (1). The full topological charge Q is then obtained by summation Since Q is on a finite lattice usually not exactly an integer, the result of the measurement will be projected to the nearest integer. It should be noted that Q does not count the number of, e.g., instanton-like objects, but it is the net number. Thus, in a configuration with n + objects of topological charge 1 and n − objects of topological charge −1, Q = n + − n − . It therefore cannot give insight on the local structure of a configuration, but merely characterizes the vacuum sector.
The topological charge is also used to monitor the smearing process: In the self-dual regime, the quantity (2) forms (almost) integer plateaus as a function of smearing sweeps, and changes only by (almost) integer jumps [41]. The appearance of these plateaus were used to assure that almost all measurements in the long smearing runs were in the selfdual domain.
The distribution in Q is Gaussian centered around Q = 0. Thus to sample even modestly large Q requires an exponentially large amount of statistics, which restricts this type of investigation only to a very limited set of lattice parameters and/or Q values.

Selection of smearing levels
The primary goal of this work is to understand the behavior of propagators in a topological, i.e. self-dual, background. That such a background is reached is exemplified in Fig. 1. It is visible that at about 300 smearing sweeps the topological charge, even for the very high value of this configuration, has equilibrated and become almost integer. 3 For smaller charger, this state is usually reached earlier.
However, in general for different lattice spacings smearing has quite different (quantitative) effects. The same is also true for different smearing techniques. Hence, a meaningful discussion of the effects requires some comparable quantity [45]. Since the main aim of the present work are the properties of the elementary degrees of freedom, the running coupling  appears to be a suitable measure to compare the impact of smearing. Its definition and other properties will be discussed below in Sect. 4.4. Here, it suffices that it is proportional to some function A( p). Thus, the measure will be given by A smeared ( p)/A( p), at some fixed p. The momentum should be such that A( p) becomes a monotonous function of the smearing sweeps. As will be seen in Sect. 4.4, this is the case at sufficiently large momenta. Since the different β yield a quite different set of momenta, there is no momentum common to all lattices. It is therefore possible to either interpolate the momenta or choose a momentum closest to a reference momentum. Since the results turn out to be only very weakly affected by the choice, here the momentum closest to the reference momentum of 1.65 GeV is chosen. The resulting smearing dependence is shown in Fig. 2.
As is visible, the decrease is for all values of β monotonous, and slower the finer the lattice. This is actually not an effect of the different volumes. The curves are only slightly quantitatively altered by changing the volumes among the available ones, a few-percent-level effect at most. Thus, in the following only the largest volumes will be considered.
It is then visible that the suppression reached for the smearing sweeps required to reach the full self-dual regime according to Fig. 1 is very large. Such a suppression is not reached on the finer lattices. However, as will be seen later the results show a rather smooth development with the number of smearing sweeps, and the behavior deep in the self-dual regime differs very little from the one obtained much earlier in the smearing progress.
Hence, the following smearing levels will be used for comparison of the different lattice spacings employed: • The smallest amount of suppression possible for all lattice spacings is 0.74, corresponding to the set of 1, 2, and 6 The levels of smearing will be used now to compare the results on the different lattices. As noted, however, the development is smooth, and therefore conclusions can be drawn irrespective of where on the smearing trajectories a result isthe development along the trajectory is the truly relevant result. Gluon dressing function, slightly smeared

Fig. 3
The gluon propagator (left panels) and dressing function (right panels) at β = 2.2 for different volumes for different numbers of APE sweeps, being slightly smeared (top panels), moderately smeared (middle panels), or strongly smeared (bottom panels). If no errors bars are visible here and hereafter, then they are smaller than the symbol size. The momenta are always along the x-axis, best suited to reach low momenta [49]

Results
The main aim of this work is to understand how the propagators behave for topological configurations. As noted before, this situation is only reached for the coarsest lattices with the available resources. However, as the following will show, the changes of the propagators are both monotonous and smooth under smearing. Especially, no qualitative change is observed when entering the regime of pure topological field configurations. It appears thus likely that this is also true for finer lattices.
As a consequence, the main observation of this work is the development under smearing. Hence, in the following the development of the correlation functions under smearing The discretization is a = 0.14 fm. All results from 24 4 lattices will be presented, and no distinction will be made between the situation at few levels of smearing and the self-dual regime, except where noted, or indicated in the figure. It will be observed that there a trends under smearing, and these trends are stable. In fact, the results in the topological regime differ from those after only few levels of smearing only quantitatively, but not qualitatively. Also this is an important result in itself.

Lattice artifacts
The simplest object to investigate is the gluon propagator. To assess the results requires to estimate the two potential types of lattice artifacts, which can affect it, the (physical) volume, and the discretization [1]. Without smearing, the dominant qualitative artifact is the volume [1], while discretization artifacts yield only a (sizable) quantitative correction; see e.g. [49][50][51]. To assess the effects, the impact of the volume for a slightly (suppression factor 1.35), moderately (suppression factor 2), and strongly (suppression factor 20) smeared gluon propagator is shown in Fig. 3. In all cases the extent of smearing is not altering the strength of the finite-volume effects significantly. The situation is somewhat different when the discretization is varied, as is visible in Fig. 4. While without smearing the ultraviolet part agrees within a few percent for these discretizations [49], already slightly smearing changes this. Then the ultraviolet tail is the stronger suppressed the coarser the discretization. This effect also increases with increasing number of APE sweeps. This is most visible at the renormalization point μ = 2 GeV, where the dressing functions coincide without smearing, but differ after strong smearing. The effect is much less pronounced at small momenta. Thus the low-momentum regime is not overmuch affected by discretization effects, but the high-momentum tail is. This is not too surprising, since ultraviolet fluctuations are most affected by the smearing operation, and thus discretization effects should become more pronounced at large momenta.

Results for the gluon propagator
The gluon propagator and dressing functions as a function of the number of smearing sweeps is shown in Figs. 5, 6, and 7. While a suppression by a small factor is not a large effect, even for a rather fine discretization, already moderate suppression factors alter the ultraviolet behavior substantially. This is not surprising, as the ultraviolet tail of the propagators are dominated by the hard, perturbative fluctuations, which are suppressed by the smearing process. It is also visible that with increased smearing the suppression starts at smaller and smaller momenta.
The situation is rather different at small momenta, and shows strong dependence on the physical volume. Concentrating on the case of the largest physical volumes, it is found that the peak of the dressing function becomes enhanced and moves to lower momenta. At the same time, the slope towards the peak becomes increased, appearing to aspire to a 1/ p 4 behavior, as shown explicitly in Fig. 8. Indeed, after a sufficient number of sweeps, but already before entering the Smearing progress Fig. 9 The gluon propagator at zero momentum as a function of volume and smearing steps. If for a volume two discretizations are available, both have been included, what especially at small volumes leads to fluctuations. The ridge structure appears due to the substantial impact of various systematic effects at such small volumes and coarse lattices; the relevant information here is merely the common trend. All results are renormalized, as discussed in Sect. 2 self-dual regime, the maximum in the dressing function has moved to so small momenta that it can almost no longer be observed on the present volumes. The finer lattices show the same behavior, but, due to the smaller volumes, the peaks disappear much earlier from sight.
Of course, the lattice volumes are small, and therefore this may change once more when moving towards very large volumes. Especially, it cannot be concluded that the peak will survive forever at smaller and smaller momenta. Nonetheless, it appears that the salient feature of the gluon propagator survives under smearing. The appearance of a 1/ p 4 behavior at typical hadronic scales is, however, interesting. Especially, as such a strong rise coincides with the naive expectation for a linear rising potential [2]. This will be discussed further in Sect. 5.
To trace out this behavior further, a helpful quantity is D(0) as a function of physical volume, which is shown in Fig. 9. Since D(0) more or less continuously increases with the number of smearing steps for all investigated volumes, except perhaps for the largest one, there is not yet any conclusive hint of the fate of the maximum in the dressing function in a topological background field. However, even the largest volumes included have to be considered small with respect to the asymptotic behavior of the gluon propagator [1], and therefore caution is mandatory.

Dependency on the topological charge
Given the observation [43] that lattice simulations algorithms tend to become stuck in a sector of fixed (net) topological charge when moving further and further towards the continuum limit, it is a relevant question whether this affects the correlation functions substantially, so that special care would be required. In general, several APE sweeps are necessary before the topological charge stabilizes itself, so to answer this question for the unsmeared case requires to extrapolate any dependency backwards. One substantial problem in doing so is that the topological charge is Gaussian distributed for configurations, and therefore very large statistics is needed to access large topological charge, especially for highly smeared configurations.
Keeping this statistical limitation in mind, results for several discretizations and smearing sweeps are shown in Figs. 10, 11, and 12. It is visible that there is essentially no effect, no matter the suppression factor nor whether in the self-dual domain or not. Thus, the gluon propagator does not appear to depend too strongly on the topological charge sector.
It should again be stressed that Q is a net-charge. Therefore, this result cannot characterize how the propagator would look, e.g., in configurations with two lumps of topological charge 1 and one with −1 in contrast to a configuration with only one of charge 1. The dependence on such a characterization may be quite different from the one found here.

Gluon propagator: position space
The gluon propagator in position space has been a valuable quantity to indirectly determine generic properties of the analytic structure of the propagator [1,3]. Though by now first direct continuum calculations using functional methods are available for it [52], the systematic uncertainties make indirect information still valuable. Since especially the long-range structure is interesting, the smearing could be expected to improve the quality of the results, similar to what is obtained for bound states [39].
The result for different smearing levels are shown for the three discretizations in Fig. 13. Some qualitative difference is observed for the different discretizations. This is mainly that the finer the discretization, the earlier a behavior is seen which appears on coarser lattice only for a stronger suppression. This effect is of the same size as the difference due to the different volumes of the not smeared case, and hence is probably rather a finite-volume artifact. It is therefore possible to concentrate on the results for the largest physical volumes, and therefore longest accessible times. It is found that the decay becomes slower the more smearing has been performed. This moves the characteristic zero crossing [1] to larger times, but it is even after substantial smearing still observable. Even in the cases where the zero crossing is not visible, and even after the longest amount of smearing, the correlator still curves incorrectly for a physical correlator. Thus, the generic properties of the  Of course, if one is willing to interpret the severely smeared configuration as the relevant structure of the con-figurations, this could also be interpreted as that positivity violations in the gluon propagator, and thus its absence from the physical state spectrum, is caused by topological effects. This interpretation would also imply that the presence of pos-  =2.5, mildly smeared β Gluon dressing function at Fig. 12 The gluon dressing function for mildly (suppression factor 2, top panels), moderately (suppression factor 10, middle panels), and strongly smeared (suppression factor 20, bottom panels) configura-tions, in different fixed topological charge sectors. The lattice spacing is a = 0.087 fm, on a 24 4 lattices QCD [1,3]. Especially, it will be important for the running coupling in the next section [2]. Before analyzing it, the first step is once more to assess the importance of lattice artifacts. The influence of both the lattice volume and the discretization for a smeared ghost dressing func-tion is shown in Fig. 14  Changing the discretization has more effect. One is on the renormalization, which is performed in Fig. 14 with the same renormalization constants as for the unsmeared case. These factors yield even for the smallest suppression no longer coinciding values. At the strongest suppression, the effect is more pronounced as it suppressed the ghost dressing function for the coarsest lattice stronger than it is on the finest lattice. Hence, the results discussed below may be bigger on a finer lattice. Also comparing the values of β = 2.35 and β = 2.5, they move closer together Smearing progress in the far infrared with increasing smearing. Thus, the impact of smearing seems to increase with increasing lattice spacing.
Keeping this in mind, the next step is to proceed to the dependence on the smearing progress. The equivalent to Figs. 5, 6, and 7 for the gluon is shown for the ghost in Figs. 15, 16, and 17. Shown is the unrenormalized ghost dressing function, and that for a particular reason: when increasing the number of APE sweeps, the high-momentum behavior is not removed, like for the gluon propagator. Rather, the behavior of the propagator becomes more and more the one of a free particle, with trivial renormalization. At the same time the qualitative infrared behavior is unaltered, and the ghost dressing function remains enhanced. However, the amount of enhancement is significantly reduced at fixed momenta. The same effect would occur, if the enhancement would be shifted to smaller momenta. Hence, this behavior is in agreement with what is observed for the gluon propagator, where its characteristic infrared behavior sets in for smaller and smaller momenta the more the con-figurations are smeared, and the main impact occurs at midmomentum. There is no distinguishing effect when reaching the self-dual domain: the reduction of infrared strength is a continuous process.
For the same reason as for the gluon propagator, it is interesting to identify the topological-(net-)charge dependence of the ghost propagator, which is shown in Figs. 18, 19, and 20. Similarly to the gluon case, no statistically significant dependency on the topological charge is observed. Therefore, the same conclusion holds as for the gluon propagator, i.e. there is no significant dependency on the topological-(net-)charge sector.
Of course, the same caveat that this is a dependence on a net-charge rather than the number of topological lumps applies here as well.

Running coupling
The running coupling is a useful quantity to assess how the structure of the fields affects the interaction strength. In Lan- where α(μ) normalizes the coupling, making α( p) independent of μ. Since this is an irrelevant (scheme-dependent) overall factor, here only the dimensionless product A = p 6 D 2 G D will be investigated, as the main interest is the change of momentum-dependence under smearing. This dependence on smearing is shown in Fig. 21. The results agree with what could have been inferred from Figs. 5, 6, 7, 15, 16, and 17: At small momenta, the decrease in the ghost propagator cannot compensate for the enhancement of the gluon propagator. Thus, the running coupling becomes a steeply increasing function at mid-momentum, before eventually bending down again, if the volume is sufficiently large to reach small enough momenta. At the same time, the ultraviolet suppression of the gluon propagator dominates the effect at large momenta, since the ghost propagator is almost not altered. As a consequence, the ultraviolet interactions are strongly suppressed, which is reflected in an ultraviolet suppressed running coupling.

Residual configuration results
In principle, it is possible to view the smeared links U s μ just as part of the original links U μ , with the remainder part defined as This implies that for an unsmeared configuration U r μ is the unit matrix, and the corresponding gluon propagator vanishes therefore identical and the ghost propagator is the one of a free particle.
Since matrix multiplication is not commutative, this is not a unique decomposition. Nonetheless, similar decompositions have been used in the past [26,27,30,31,33,34] to characterize the contents of the part of the configuration removed under smearing and/or other operations to isolate the topological content. Based on the assumption that smear- Smearing progress ing removes only the ultraviolet fluctuations, the corresponding propagators have usually been found to be approximately the perturbative or tree-level ones. In Fig. 22 the corresponding propagators and dressing functions for different numbers of APE sweeps are shown. Surprisingly, these function agree within a few percent, after appropriate normalization and after a certain amount of smearing. The overall normalization is a decreasing function of the number of APE sweeps, which appears to tend to a finite value for larger and larger numbers of APE sweeps, though this has not been studied in detail. Also, there is a significant impact, especially for the ghost dressing function, of the discretization. The finer and smaller the lattice, the closer the ghost dressing function at the same amount of suppression moves toward the tree-level behavior.
The results are as expected. The gluon propagator is essentially momentum independent, and thus the one of a random variable, A a μ (x)A a ν (y) = δ ab δ μν δ(x − y). Consequently, the dressing function increases quickly with momentum. Such a behavior is also seen in the Schwinger function: It is zero, within error, except for x = y. The ghost dressing function shows only a mild deviation from the one of a free particle at large momentum for an increasing level of smearing. Otherwise, it is essentially that of a free, massless particle.

A speculative interpretation
The results shown here are for a very limited amount of different volumes and discretizations, due to the computational costs involved. Given the sensitivity of the investigated correlation functions to both types of lattice artifacts [1], any interpretation should only be made with the requirement that it has to be confirmed for larger volumes, and thus remains currently speculative. It should also be kept in mind that any other way of selecting Gribov copies than the minimal Landau gauge used here may yield different results.
Assuming for a moment that there will be no qualitative change for larger volumes leads then to the following interpretation: also for propagators smearing leads to a suppres-  18 The ghost dressing function for mildly (suppression factor 2, top panels), moderately (suppression factor 85, middle panels), and strongly smeared (self-dual regime, bottom panels) configura-tions, in different fixed topological charge sectors. The discretization is a = 0.21 fm. All results from 24 4 lattices sion of ultraviolet physics. Consequently, the gluon propagator is ultraviolet suppressed. The corresponding residual propagator is the one of a free random variable. The ghost propagator shows at large momenta the one of a free field, and it thus decouples from the gluonic fluctuations. At the same time, the infrared and mid-momentum behavior is nontrivial. Thus, the long-distance behavior of the propagators appears to be dominated by self-dual (topological) configurations, in line with the arguments in [28,29,31]. This also confirms the idea that topological scenarios and the physics  19 The ghost dressing function for mildly (suppression factor 2, top panels), moderately (suppression factor 20, middle panels), and strongly smeared (suppression factor 85, bottom panels) configura-tions, in different fixed topological charge sectors. The discretization is a = 0.14 fm. All results from 24 4 lattices captured by the infrared behavior of correlation functions are just two facets of the same underlying physics [25].
A somewhat surprising result is the appearance of a range of momenta at the hadronic scales where the gluon propagator appears to behave like 1/ p 4 . This would be a welcome effect, as such a contribution makes it much easier to understand where the gluon propagator encodes information like the Wilson string. Taking this for granted, this would imply that by smearing the physical irrelevant, gauge-fixingdominated infrared part is shifted to small momenta, and the relevant mid-momentum physics, which is thus generated by the topological excitations, is made evident. This would  The ghost dressing function for mildly (suppression factor 2, top panels), moderately (suppression factor 10, middle panels), and strongly smeared (suppression factor 20, bottom panels) configura-tions, in different fixed topological charge sectors. The discretization is a = 0.087 fm. All results from 24 4 lattices imply, in a very simplified manner, that the unsmeared gluon propagator consists out of three parts: A ultraviolet part, a strongly interacting 1/ p 4 part, and an infrared finite part due to the gauge-fixing process. These are schematically multiplicatively superimposed like where γ is the Gribov parameter [6], and D p is just a perturbative part. Smearing removes the ultraviolet part, and   The top panels are for a discretization of a = 0.21 fm, the middle panels for a = 0.14 fm, and the bottom panels for a = 0.087 fm. All results from 24 4 lattices. Note the scales Besides this conceptual insight, this has also practical implications for non-lattice calculations of correlations functions, especially for functional methods. Any approximation and/or truncation made in such calculations can in principle remove part or all of the physics due to topological configurations. The investigations here now provide a necessary condition for the correct inclusion of this infrared physics. Only if the infrared and the mid-momentum behavior of the correlation functions is correct, the contributions from topology is included. This is not sufficient, at least if only a finite number of correlation functions is determined, since it appears possible that the infrared behavior could also be reproduced without the correct physics.
On the other hand, there is a remarkable implication in return. It has been explicitly tested that for many physical observables the far infrared behavior of correlation function appears to be irrelevant [55][56][57][58][59]. Since this part appears to be retained under smearing, this implies that this unphysical information is still encoded in the topological fluctuations as well, besides any physical information. Given the limited amount of results here, also this statement has to be taken as being speculative.

Conclusions
Summarizing, it is found that for the limited set of lattice settings studied here the low-momentum behavior of correlation functions is qualitatively conserved under APE-smearing. That implies that likely the infrared part of correlation functions carry indeed the imprint of topological configurations with them, emphasizing that both aspects are only two facets of the same underlying physics. Furthermore, no pronounced dependence of the correlation functions on (net-) topological charge is found. This implies that the correlation functions are not strongly affected if lattice simulation algorithms get stuck in a sector of fixed topological charge.
Of course, a precise quantitative investigation will require much larger lattices, especially in the asymptotic domain [1]. However, such an investigation will require substantially more computational resources, and therefore will only be possible when these become available in the future. Still, the present results provide besides the physical insights the necessary requirement that analytical determinations of correlation functions need to be correct at small and intermediate momenta to include topological effects. On the other hand, since the configurations are purely self-dual/topological after a sufficient number of APE sweeps, this implies that correlation functions can capture the physics of such configurations, and no artificial additional introduction of them into, e.g., functional methods is required.