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

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

Insight into nucleon structure from generalized parton distributions

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.Comment: 3 pages, Lattice2003(Theoretical developments

Moments of nucleon spin-dependent generalized parton distributions

We present a lattice measurement of the first two moments of the spin-dependent GPD H-tilde(x,xi,t). From these we obtain the axial coupling constant and the second moment of the spin-dependent forward parton distribution. The measurements are done in full QCD using Wilson fermions. In addition, we also present results from a first exploratory study of full QCD using Asqtad sea and domain-wall valence fermions.Comment: Lattice2003(Theory), 3 pages, 3 figures, to appear in the Proceedings of Lattice 200

Higher Order Processes in Electromagnetic Production of Electron Positron Pairs in Relativistic Heavy Ion Collisions

We study higher-order effects in the electromagnetic production of electron-positron pairs in relativistic heavy ion collisions. Treating the field of the heavy ions as an external field and neglecting the interaction among electrons and positrons, we show that the $N$-pair creation amplitude is the antisymmetrised product of $N$ one-pair creation amplitudes and the vacuum amplitude. Neglecting contributions coming from exchange terms, we show that the total probability for $N$ pairs is approximately a Poisson distribution. We investigate further the structure of the reduced one-pair amplitude, concentrating especially on multiple-particle corrections. We calculate the first of these corrections in second order Magnus theory based on our previous result in second-order Born approximation for impact parameter $b$ zero. Explicit calculations show that the total probability is increased up to 10 \% by this correction for realistic collider parameters. The calculations can also be used to confirm the use of the Poisson distribution for the total probability.Comment: 29 pages RevTeX and 12 uuencoded figures (compressed postscript

Numerical atomic orbitals for linear scaling

The performance of basis sets made of numerical atomic orbitals is explored in density-functional calculations of solids and molecules. With the aim of optimizing basis quality while maintaining strict localization of the orbitals, as needed for linear-scaling calculations, several schemes have been tried. The best performance is obtained for the basis sets generated according to a new scheme presented here, a flexibilization of previous proposals. The basis sets are tested versus converged plane-wave calculations on a significant variety of systems, including covalent, ionic and metallic. Satisfactory convergence (deviations significantly smaller than the accuracy of the underlying theory) is obtained for reasonably small basis sizes, with a clear improvement over previous schemes. The transferability of the obtained basis sets is tested in several cases and it is found to be satisfactory as well.Comment: 9 pages with 2 encapsulated postscript figures, submitted to Phys. Rev.

Systematic generation of finite-range atomic basis sets for linear-scaling calculations

Basis sets of atomic orbitals are very efficient for density functional calculations but lack a systematic variational convergence. We present a variational method to optimize numerical atomic orbitals using a single parameter to control their range. The efficiency of the basis generation scheme is tested and compared with other schemes for multiple zeta basis sets. The scheme shows to be comparable in quality to other widely used schemes albeit offering better performance for linear-scaling computations

Multicanonical Hybrid Monte Carlo: Boosting Simulations of Compact QED

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.Comment: 15 pages, 8 figure

Chiral extrapolation of lattice moments of proton quark distributions

We present the resolution of a long-standing discrepancy between the moments of parton distributions calculated from lattice QCD and their experimental values. We propose a simple extrapolation formula for the moments of the nonsinglet quark distribution u-d, as a function of quark mass, which embodies the general constraints imposed by the chiral symmetry of QCD. The inclusion of the leading nonanalytic behavior leads to an excellent description of both the lattice data and the experimental values of the moments.Comment: 9 pages, 1 figure, to appear in Physical Review Letter

Unquenched QCD with Light Quarks

We present recent results in unquenched lattice QCD with two degenerate light sea quarks using the truncated determinant approximation (TDA). In the TDA the infrared modes contributing to the quark determinant are computed exactly up to some cutoff in quark off-shellness (typically 2$\Lambda_{QCD}$). This approach allows simulations to be performed at much lighter quark masses than possible with conventional hybrid MonteCarlo techniques. Results for the static energy and topological charge distributions are presented using a large ensemble generated on very coarse (6$^4$) but physically large lattices. Preliminary results are also reported for the static energy and meson spectrum on 10$^3$x20 lattices (lattice scale $a^{-1}$=1.15 GeV) at quark masses corresponding to pions of mass $\leq$ 200 MeV. Using multiboson simulation to compute the ultraviolet part of the quark determinant the TDA approach becomes an exact with essentially no increase in computational effort. Some preliminary results using this fully unquenched algorithm are presented.Comment: LateX, 39 pages, 16 eps figures, 1 ps figur