FFT for the APE Parallel Computer

We present a parallel FFT algorithm for SIMD systems following the `Transpose Algorithm' approach. The method is based on the assignment of the data field onto a 1-dimensional ring of systolic cells. The systolic array can be universally mapped onto any parallel system. In particular for systems with next-neighbour connectivity our method has the potential to improve the efficiency of matrix transposition by use of hyper-systolic communication. We have realized a scalable parallel FFT on the APE100/Quadrics massively parallel computer, where our implementation is part of a 2-dimensional hydrodynamics code for turbulence studies. A possible generalization to 4-dimensional FFT is presented, having in mind QCD applications.Comment: 17 pages, 13 figures, figures include

An Estimate of alpha_S from Bottomonium in Unquenched QCD

We estimate the strong coupling constant from the perturbative expansion of the plaquette. The scale is set by the 2S-1S and 1P-1S splittings in bottomonium which are computed in NRQCD on dynamical gauge configurations with nf=2 degenerate Wilson quarks at intermediate masses. We have increased the statistics of our spectrum calculation in order to reliably extrapolate in the sea-quark mass. We find a value of alpha_MS(m_Z) = 0.1118(26) which is somewhat lower than previous estimates within NRQCD.Comment: LATTICE98(heavyqk

String breaking

We numerically investigate the transition of the static quark-antiquark string into a static-light meson-antimeson system. Improving noise reduction techniques, we are able to resolve the signature of string breaking dynamics for Nf=2 lattice QCD at zero temperature. We discuss the lattice techniques used and present results on energy levels and mixing angle of the static two-state system. We visualize the action density distribution in the region of string breaking as a function of the static colour source-antisource separation. The results can be related to properties of quarkonium systems.Comment: 8 pages, Talk given at the Workshop on Computational Hadron Physics, Nicosia, Cyprus, 14--17 September 200

The Consequences of Non-Normality

The non-normality of Wilson-type lattice Dirac operators has important consequences - the application of the usual concepts from the textbook (hermitian) quantum mechanics should be reconsidered. This includes an appropriate definition of observables and the refinement of computational tools. We show that the truncated singular value expansion is the optimal approximation to the inverse operator D^{-1} and we prove that due to the gamma_5-hermiticity it is equivalent to gamma_5 times the truncated eigenmode expansion of the hermitian Wilson-Dirac operator

Unquenching the Quark Model and Screened Potentials

The low-lying spectrum of the quark model is shown to be robust under the effects of `unquenching'. In contrast, the use of screened potentials is shown to be of limited use in models of hadrons. Applications to unquenching the lattice Wilson loop potential and to glueball mixing in the adiabatic hybrid spectrum are also presented.Comment: 6 pages, 3 ps figures, revtex. Version to appear in J. Phys.

Flavour singlet pseudoscalar masses in N_f = 2 QCD

We perform a lattice mass analysis in the flavour singlet pseudoscalar channel on the SESAM and TXL full QCD vacuum configurations, with 2 active flavours of dynamical Wilson fermions at beta = 5.6. At our inverse lattice spacing, a^-1 = 2.3 GeV, we retrieve by a chiral extrapolation to the physical light quark masses the value m_eta' = 3.7(+8)(-4) m_pi. A crude extrapolation from (N_f = 3) phenomenology would suggest m_eta' \approx 5.1 m_pi for N_f = 2 QCD. we verify that the mass gap between the singlet state eta' and the pi flavour triplt state is due to gauge configurations with non-trivial topology.Comment: 8 pages, 10 figure

Hyper-systolic processing on APE100\/Quadrics: n2^{2}-loop computations

We investigate the performance gains from hyper-systolic implementations of n^2-loop problems on the massively parallel computer Quadrics, exploiting its 3-dimensional interprocessor connectivity. For illustration we study the communication aspects of an exact molecular dynamics simulation of n particles with Coulomb (or gravitational) interactions. We compare the interprocessor communication costs of the standard-systolic and the hyper-systolic approaches for various granularities. We predict gain factors as large as 3 on the Q4 and 8 on the QH4 and measure actual performances on these machine configurations. We conclude that it appears feasile to investigate the thermodynamics of a full gravitating n-body problem with O(10000) particles using the new method on a QH4 system

Instanton dominance of topological charge fluctuations in QCD?

We consider the local chirality of near-zero eigenvectors from Wilson-Dirac and clover improved Wilson-Dirac lattice operators as proposed recently by Horv\'ath et al. We studied finer lattices and repaired for the loss of orthogonality due to the non-normality of the Wilson-Dirac matrix. As a result we do see a clear double peak structure on lattices with resolutions higher than 0.1 fm. We found that the lattice artifacts can be considerably reduced by exploiting the biorthogonal system of left and right eigenvectors. We conclude that the dominance of instantons on topological charge fluctuations is not ruled out by local chirality measurements.Comment: 10 pages, 6 figure

Scanning the Topological Sectors of the QCD Vacuum with Hybrid Monte Carlo

We address a long standing issue and determine the decorrelation efficiency of the Hybrid Monte Carlo algorithm (HMC), for full QCD with Wilson fermions, with respect to vacuum topology. On the basis of five state-of-the art QCD vacuum field ensembles (with 3000 to 5000 trajectories each and m_pi/m_rho-ratios in the regime >0.56, for two sea quark flavours) we are able to establish, for the first time, that HMC provides sufficient tunneling between the different topological sectors of QCD. This will have an important bearing on the prospect to determine, by lattice techniques, the topological susceptibility of the vacuum, and topology sensitive quantities like the spin content of the proton, or the eta' mass.Comment: 5 pages, 4 eps-figure

Multipole Expansion of Bremsstrahlung in Intermediate Energy Heavy Ion Collisions

Using a multipole expansion of the radiated field generated by a classical electric current, we present a way to interprete the bremsstrahlung spectra of low energy heavy ion collisions. We perform the calculation explicitely for the system ^{12}C+ ^{12}C at 84AMeV and compare the result with the experimental data of E. Grosse et al. Using simple model assumptions for the electromagnetic source current we are able to describe the measured data in terms of coherent photon emission. In this context, the information contained in the measured data is discussed.Comment: LaTex, 4 Figure