Meson and Baryon dispersion relations with Brillouin fermions

We study the dispersion relations of mesons and baryons built from Brillouin quarks on one N_f=2 gauge ensemble provided by QCDSF. For quark masses up to the physical strange quark mass, there is hardly any improvement over the Wilson discretization, if either action is link-smeared and tree-level clover improved. For quark masses in the range of the physical charm quark mass, the Brillouin action still shows a perfect relativistic behavior, while the Wilson action induces severe cut-off effects. As an application we determine the masses of the \Omega_c^0, \Omega_{cc}^+ and \Omega_{ccc}^{++} baryons on that ensemble.Comment: 16 pages, 9 figures, 4 tables; v2: one Reference added, matches published versio

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

Cost of dynamical quark simulations: O(a) improved Wilson fermions

I report on cost estimates and algorithmic performance in simulations using 2 flavours of non-perturbatively O(a) improved Wilson quarks together with the Wilson plaquette action.Comment: Contribution to Lattice2001 (panel discussion), 2 pages, 2 figure

Electron transport through an interacting region: The case of a nonorthogonal basis set

The formula derived by Meir and Wingreen [Phys. Rev. Lett. {\bf 68}, 2512 (1992)] for the electron current through a confined, central region containing interactions is generalized to the case of a nonorthogonal basis set. As in the original work, the present derivation is based on the nonequilibrium Keldysh formalism. By replacing the basis functions of the central region by the corresponding elements of the dual basis, the lead- and central region-subspaces become mutually orthogonal. The current formula is then derived in the new basis, using a generalized version of second quantization and Green's function theory to handle the nonorthogonality within each of the regions. Finally, the appropriate nonorthogonal form of the perturbation series for the Green's function is established for the case of electron-electron and electron-phonon interactions in the central region.Comment: Added references. 8 pages, 1 figur

Proposal for an interference experiment to test the applicability of quantum theory to event-based processes

We analyze a single-particle Mach-Zehnder interferometer experiment in which the path length of one arm may change (randomly or systematically) according to the value of an external two-valued variable $x$, for each passage of a particle through the interferometer. Quantum theory predicts an interference pattern that is independent of the sequence of the values of $x$. On the other hand, corpuscular models that reproduce the results of quantum optics experiments carried out up to this date show a reduced visibility and a shift of the interference pattern depending on the details of the sequence of the values of $x$. The proposed experiment will show that: (1) it can be described by quantum theory, and thus not by the current corpuscular models, or (2) it cannot be described by quantum theory but can be described by the corpuscular models or variations thereof, or (3) it can neither be described by quantum theory nor by corpuscular models. Therefore, the proposed experiment can be used to determine to what extent quantum theory provides a description of observed events beyond the usual statistical level.Comment: Accepted for publication in J. Phys. Soc. Jp

A comparative study of numerical methods for the overlap Dirac operator--a status report

Improvements of various methods to compute the sign function of the hermitian Wilson-Dirac matrix within the overlap operator are presented. An optimal partial fraction expansion (PFE) based on a theorem of Zolotarev is given. Benchmarks show that this PFE together with removal of converged systems within a multi-shift CG appears to approximate the sign function times a vector most efficiently. A posteriori error bounds are given.Comment: 3 pages, poster contribution to Lattice2001(algorithms

Improved Upsilon Spectrum with Dynamical Wilson Fermions

We present results for the b \bar b spectrum obtained using an O(M_bv^6)-correct non-relativistic lattice QCD action, where M_b denotes the bare b-quark mass and v^2 is the mean squared quark velocity. Propagators are evaluated on SESAM's three sets of dynamical gauge configurations generated with two flavours of Wilson fermions at beta = 5.6. These results, the first of their kind obtained with dynamical Wilson fermions, are compared to a quenched analysis at equivalent lattice spacing, beta = 6.0. Using our three sea-quark values we perform the ``chiral'' extrapolation to m_eff = m_s/3, where m_s denotes the strange quark mass. The light quark mass dependence is found to be small in relation to the statistical errors. Comparing the full QCD result to our quenched simulation we find better agreement of our dynamical data with experimental results in the spin-independent sector but observe no unquenching effects in hyperfine-splittings. To pin down the systematic errors we have also compared quenched results in different ``tadpole'' schemes as well as using a lower order action. We find that spin-splittings with an O(M_bv^4) action are O(10%) higher compared to O(M_bv^6) results. Relative to the results obtained with the plaquette method the Landau gauge mean link tadpole scheme raises the spin splittings by about the same margin so that our two improvements are opposite in effect.Comment: 24 pages (latex file, Phys Rev D style file, uses epsf-style

String breaking with dynamical Wilson fermions

We present results of our ongoing determination of string breaking in full QCD with N_f=2 Wilson fermions. Our investigation of the fission of the static quark-antiquark string into a static-light meson-antimeson system is based on dynamical configurations of size 24^3 x 40 produced by the TxL collaboration. Combining various optimization methods we determine the matrix elements of the two-by-two system with so far unprecedented accuracy. The all-to-all light quark propagators occurring in the transition element are computed from eigenmodes of the Hermitian Wilson-Dirac matrix complemented by stochastic estimates in the orthogonal subspace. We observe a clear signature for level-splitting between ground state and excited potential. Thus, for the first time, string breaking induced by sea quarks is observed in a simulation of 4-dimensional lattice-QCD.Comment: 3 pages, 3 figures, contribution to Lattice 200

Compression domain volume rendering for distributed environments

This paper describes a method for volume data compression and rendering which bases on wavelet splats. The underlying concept is especially designed for distributed and networked applications, where we assume a remote server to maintain large scale volume data sets, being inspected, browsed through and rendered interactively by a local client. Therefore, we encode the server‘s volume data using a newly designed wavelet based volume compression method. A local client can render the volumes immediately from the compression domain by using wavelet footprints, a method proposed earlier. In addition, our setup features full progression, where the rendered image is refined progressively as data comes in. Furthermore, framerate constraints are considered by controlling the quality of the image both locally and globally depending on the current network bandwidth or computational capabilities of the client. As a very important aspect of our setup, the client does not need to provide storage for the volume data and can be implemented in terms of a network application. The underlying framework enables to exploit all advantageous properties of the wavelet transform and forms a basis for both sophisticated lossy compression and rendering. Although coming along with simple illumination and constant exponential decay, the rendering method is especially suited for fast interactive inspection of large data sets and can be supported easily by graphics hardware