BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
K. Schilling appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
The lowest three moments of generalized parton distributions are calculated in full QCD and provide new insight into the behavior of nucleon electromagnetic form factors, the origin of the nucleon spin, and the transverse structure of the nucleon.
This work presents the first calculation in lattice QCD of three moments of spin-averaged and spin-polarized generalized parton distributions in the proton. It is shown that the slope of the associated generalized form factors decreases significantly as the moment increases, indicating that the transverse size of the light-cone quark distribution decreases as the momentum fraction of the struck parton increases.
The fomalism is developed to express nucleon matrix elements of the electromagnetic current in terms of form factors consistent with the translational, rotational, and parity symmetries of a cubic lattice. We calculate the number of these form factors and show how appropriate linear combinations approach the continuum limit.
Moments of light cone quark density, helicity, and transversity distributions are calculated in unquenched lattice QCD at $β= 5.5$ and $β= 5.3$ using Wilson fermions on $ 16^3 \times 32 $ lattices. These results are combined with earlier calculations at $β= 5.6$ using SESAM configurations to study the continuum limit.
Moments of the quark density, helicity, and transversity distributions are calculated in unquenched lattice QCD. Calculations of proton matrix elements of operators corresponding to these moments through the operator product expansion have been performed on $16^3 \times 32$ lattices for Wilson fermions at $β= 5.6$ using configurations from the SESAM collaboration and at $β= 5.5$ using configurations from SCRI. One-loop perturbative renormalization corrections are included. At quark masses accessible in present calculations, there is no statistically significant difference between quenched and full QCD results, indicating that the contributions of quark-antiquark excitations from the Dirac Sea are small. Close agreement between calculations with cooled configurations containing essentially only instantons and the full gluon configurations indicates that quark zero modes associated with instantons play a dominant role. Naive linear extrapolation of the full QCD calculation to the physical pion mass yields results inconsistent with experiment. Extrapolation to the chiral limit including the physics of the pion cloud can resolve this discrepancy and the requirements for a definitive chiral extrapolation are described.
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.
We present a report on the status of the GRAL project (Going Realistic And Light), which aims at simulating full QCD with two dynamical Wilson quarks below the vector meson decay threshold, m_ps/m_v < 0.5, making use of finite-size-scaling techniques.
In order to construct improved overlap fermions, we start from a short ranged approximate Ginsparg-Wilson fermion and insert it into the overlap formula. We show that its polynomial evaluation is accelerated considerably compared to the standard Neuberger fermion. In addition the degree of locality is strongly improved.
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.
We demonstrate that substantial progress can be achieved in the study of the phase structure of 4-dimensional compact QED by a joint use of hybrid Monte Carlo and multicanonical algorithms, through an efficient parallel implementation. This is borne out by the observation of considerable speedup of tunnelling between the metastable states, close to the phase transition, on the Wilson line. Our approach leads to a general parallelization scheme for the efficient stochastic sampling of systems where (a part of) the Hamiltonian involves the total action or energy in each update step.
It is shown that the 21264 Alpha processor can reach about 20% sustained efficiency for the inversion of the Wilson-Dirac operator. Since fast ethernet is not sufficient to get balancing between computation and communication on reasonable lattice- and system-sizes, an interconnection using Myrinet is discussed. We find a price/performance ratio comparable with state-of-the-art SIMD-systems for lattice QCD.
Considerable progress has been achieved recently in the observation of string breaking within non-Abelian Higgs models, by use of multi-channel methods allowing for broken string states. Similarly, in pure gauge theory this approach has been shown to reveal string breaking for color charges in the adjoint represaentation. For QCD with dynmical fermions, one needs substantial progress in noise reduction, however, to render such techniques viable.
We present the final analysis of the light and strange hadron spectra from a full QCD lattice simulation with two degenerate dynamical sea quark flavours at $β= 5.6$ on a $16^3 \times 32$ lattice. Four sets of sea quark masses corresponding to the range $.69 \leq m_π/m_ρ\leq .83$ are investigated. For reference we also ran a quenched simulation at $β_{\sf eff} = 6.0$, which is the point of equal lattice spacing, $a_ρ^{-1}$. In the light sector, we find the chiral extrapolation to physical u- and d- masses to present a major source of uncertainty, comparable to the expected size of unquenching effects. From linear and quadratic fits we can estimate the errors on the hadron masses made from light quarks to be on a 15 % level prior to the continuum extrapolation. For the hadrons with strange valence quark content, the $N_F = 2$ approximation to QCD appears not to cure the well-known failure of quenched QCD to reproduce the physical $K-K^*$ splitting.
We demonstrate that substantial progress can be achieved in the study of the phase structure of 4-dimensional compact QED by a joint use of hybrid Monte Carlo and multicanonical algorithms, through an efficient parallel implementation. This is borne out by the observation of considerable speedup of tunnelling between the metastable states, close to the phase transition, on the Wilson line. We estimate that the creation of adequate samples (with order 100 flip-flops) becomes a matter of half a year's runtime at 2 Gflops sustained performance for lattices of size up to 24^4.
We calculate connected and disconnected contributions to the flavour singlet scalar density amplitude of the nucleon in a full QCD lattice simulation with $n_f=2$ dynamical Wilson fermions at $β=5.6$ on a $16^3 \times 32$ lattice. We find that both contributions are of similar size at the light quark mass. We arrive at the estimate $σ_{πN} = 18(5)$MeV. Its smallness is directly related to the apparent decrease of $u$, $d$ quark masses when unquenching QCD lattice simulations. The $y$ parameter can be estimated from a semi-quenched analysis, in which there are no strange quarks in the sea, the result being $y=0.59(13)$.}
We investigate the critical dynamics of the Hybrid Monte Carlo algorithm approaching the chiral limit of standard Wilson fermions. Our observations are based on time series of lengths O(5000) for a variety of observables. The lattice sizes are 16^3 x 32 and 24^3 x 40. We work at beta=5.6, and kappa=0.156, 0.157, 0.1575, 0.158, with 0.83 > m_pi/m_rho > 0.55. We find surprisingly small integrated autocorrelation times for local and extended observables. The dynamical critical exponent $z$ of the exponential autocorrelation time is compatible with 2. We estimate the total computational effort to scale between V^2 and V^2.25 towards the chiral limit.
We present results for spectroscopy, quark masses and decay constants obtained from SESAM's and TkL's large statistics simulations of QCD with two dynamical Wilson fermions.
In a full QCD lattice study with $N_f = 2$ Wilson fermions, we seek to optimize the signals for the disconnected contributions to the matrix element of flavour-singlet operators between nucleon states, which are indicative for sea quark effects. We demonstrate, in form of a fluctuation analysis to the noisy estimator technique, that -- in order to achieve a tolerable signal to noise-ratio in full QCD -- it is advantageous to work with a $Z_2$-noise source rather than to rely only on gauge invariance to cancel non-gauge-invariant background. In the case of the $π$N $σ$-term, we find that 10 $Z_2$-noise sources suffice on our sample ( about 150 independent QCD configurations at $β= 5.6$ on $16^3\times32$ with $κ_{sea} = 0.157$, equivalent to $M_π/M_ρ = 0.76(1)$), to achieve decent signals and adequate fluctuations, rather than 300 such sources as recently used in quenched simulations.
The investigation of light sea-quark effects in lattice QCD with dynamical Wilson fermions requires both larger physical volumes and finer lattice resolutions than achieved previously. As high-end supercomputers like the 512-node APE Tower provide the compute power to perform a major step towards the chiral limit (T-chi-L), we have launched a feasibility study on a 24^3 x 40 lattice. We approach the chiral limit--while refining the resolution--, using the standard Wilson fermion action. Following previous work, our Hybrid Monte Carlo simulation runs at beta=5.6 and two kappa-values, 0.1575 and 0.158. From our study, we are confident that, for the APE Tower, a realistic working point has been found corresponding to a volume of 2 fm^3, with chirality characterized by 1/(a m_pi) = 5.6.
We present recent results of SESAM's large scale lattice simulation of QCD with two dynamical flavours of Wilson fermions. The emphasis is on future prospects in the extraction of flavour-singlet matrix elements, i.e the pi-nucleon sigma term, with dynamical configurations.
Pure gauge compact QED on hypercubic lattices is considered with periodically closed monopole currents suppressed. We compute observables on sublattices which are nested around the centre of the lattice in order to locate regions where translation symmetry is approximately recovered. Our Monte Carlo simulations on $24^4$-lattices give indications for a first-order nature of the U(1) phase transition.
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 feasable to investigate the thermodynamics of a full gravitating n-body problem with O(10000) particles using the new method on a QH4 system.
We present evidence that the existence of a first order phase transition in compact U(1) with Wilson action is not related to monopole loops wrapping around the toroidal lattice, as has been previously suggested. Our analysis is based on the suppression of such loops by `soft boundary conditions' that correspond to an infinitely large chemical potential for the monopoles on the boundary, during the updating process. It is observed that the double peak structure characteristic for the first order phase transition reappears at sufficiently large lattice sizes and separations from the lattice boundary.
We present intermediate results on our ongoing investigation concerning semileptonic decays of heavy pseudoscalar mesons into pseudoscalar and vector mesons. The corresponding formfactors are evaluated at several momenta and appropriate combinations of four light and four heavy quarks, which are chosen to allow for an extrapolation into the B Meson region. In order to obtain clear groundstate signals we apply gauge invariant ``Wuppertal'' smearing to the quarks. The analysis is based on 32 quenched gauge configurations of size $24^3 \times 64$ at $β=6.3$, with Wilson fermions.
The stabilized biconjugate gradient algorithm BiCGStab recently presented by van der Vorst is applied to the inversion of the lattice fermion operator in the Wilson formulation of lattice Quantum Chromodynamics. Its computational efficiency is tested in a comparative study against the conjugate gradient and minimal residual methods. Both for quenched gauge configurations at beta= 6.0 and gauge configurations with dynamical fermions at beta=5.4, we find BiCGStab to be superior to the other methods. BiCGStab turns out to be particularly useful in the chiral regime of small quark masses.
We present results on spin dependent potentials from lattice simulations of SU(2) gauge theory. The Coulomb like short range part of the central potential is identified as a mixed vector-scalar exchange while the linear long range part is pure scalar.
We present an extensive study on the direct determination of the running coupling alpha_s from the static quark antiquark force at short distances, in quenched QCD. We find from our high statistics potential analysis that alpha_qq exhibits two-loop asymptotic behaviour for momenta as low as .5 GeV. As a result, we determine the zero flavour Lambda-parameter to be Lambda^0_MSBAR = 0.630(38)\sqrtσ = 293(18)^{+25}_{-63} MeV. A rough estimate of full QCD effects leads to the five flavour value alpha_MSBAR(m_Z) = .102^{+6}_{-11}. A comparison with other lattice results is made.
A detailed investigation of the temperature dependence of the spatial string tension $σ_s$ in $SU(2)$ gauge theory is presented. A sustained performance of 3~GFLOPS on a 64K Connection Machine CM-2 equivalent has been achieved. Scaling of $σ_s$ between $β=2.5115$ and $β=2.74$, on large lattices, is demonstrated. Below the critical temperature, $T_c$, $σ_s$ remains constant. For temperatures larger than $2T_c$ the temperature dependence can be parametrized by $σ_s(T) = (0.369\pm 0.014)^2 g^4(T)T^2$, where $g(T)$ is a 2-loop running coupling constant with the scale parameter determined as $Λ_T = (0.076\pm 0.013)T_c$.
We describe the Wuppertal university pilot project in applied parallel computing. We report on a comprehensive high statistics determination of the static quark-antiquark potential and related quantities from quenched quantum chromodynamics. New data for the string tension and the plaquette action for the region 5.5 < beta < 6.8 is presented.