Source author record

Waseem Kamleh

Waseem Kamleh appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

63works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

63 published item(s)

preprint2022arXiv

Centre Vortex Structure of QCD-Vacuum Fields and Confinement

The non-trivial ground-state vacuum fields of QCD form the foundation of matter. Using modern visualisation techniques, this presentation examines the microscopic structure of these fields. Of particular interest are the centre vortices identified within the ground-state fields of lattice QCD. Our current focus is on understanding the manner in which dynamical fermions in the QCD vacuum alter the centre-vortex structure. The impact of dynamical fermions is significant and provides new insights into the role of centre vortices in underpinning both confinement and dynamical chiral symmetry breaking in QCD.

preprint2022arXiv

Centre Vortices in the Presence of Dynamical Fermions

The behaviour of centre vortices in the presence of dynamical fermions is studied for the first time in the context of the static quark potential and the Landau-gauge gluon propagator. These results indicate that in the presence of dynamical fermions, centre vortices are able to better encapsulate the long-range behaviour of QCD, compared to previous studies on pure Yang-Mills lattices. This work provides strong evidence that centre vortices are responsible for the long-range linear potential and the enhanced infrared gluon propagator, which are indicative of confinement in QCD.

preprint2022arXiv

Disconnected contributions to the magnetic polarisability of the neutral pion

The magnetic polarisability of the neutral pion has been calculated using the background field method in lattice QCD. These early results do not consider the effect of the disconnected loop contractions arising from the breaking of charge symmetry in a background magnetic field. Recent work in chiral perturbation theory has shown that these quark-self-annihilation contractions provide the leading loop-order contributions to the magnetic polarisability. A first investigation of these contractions in a background magnetic field is presented.

preprint2022arXiv

Dynamical fermions, centre vortices, and emergent phenomena

The non-trivial ground-state vacuum fields of QCD form the foundation of matter. Here we examine the centre vortices identified within the ground-state fields of lattice QCD. We aim to understand the manner in which dynamical fermions in the QCD vacuum alter the centre-vortex structure. Using modern visualisation techniques, the centre-vortex structure of pure-gauge and dynamical-fermion fields is quantified and compared. We then explore the impact this modified structure has on measures of confinement and dynamical mass generation. The string tension of the static quark potential, positivity-violation in the gluon propagator, and dynamical mass generation in the overlap quark propagator are of particular interest. The impact of dynamical fermions is significant and provides new insights into the role of centre vortices in underpinning both confinement and dynamical chiral symmetry breaking in QCD.

preprint2022arXiv

Emergent phenomena from centre vortices

Quark confinement is perhaps the most important emergent property of the theory of quantum chromodynamics. Herein we review some key aspects of centre vortices in SU(3) lattice gauge theory. Starting from the original Monte Carlo gauge fields, a vortex identification procedure yields vortex-removed and vortex-only backgrounds. The comparison between the original `untouched' Monte Carlo gauge fields and these so called vortex-modified ensembles has provided a variety of results that support the notion that centre vortices are fundamental to confinement in pure gauge theory. For the first time we perform direct numerical tests of the response of centre vortices to the presence of dynamical quarks in SU(3).

preprint2022arXiv

Impact of centre vortex removal on the Landau-gauge quark propagator in dynamical QCD

The important role of centre vortices in dynamical chiral symmetry breaking and corresponding dynamical mass generation has been demonstrated in pure gauge studies of the Landau-gauge quark propagator. We present initial results of our investigation into the impact of centre vortex removal on the Landau-gauge propagator on dynamical gauge fields. To ensure sensitivity to the topology of the gauge fields, the propagator is computed with overlap fermions. Upon removal of centre vortices we find that dynamical mass generation is suppressed. The quark renormalisation function remains remarkably flat except in the deep infrared where it exhibits significant suppression. These effects become more prominent at lighter quark masses.

preprint2022arXiv

Impact of Dynamical Fermions on Centre Vortex Structure

This presentation examines the centre vortex structure of Monte-Carlo generated gauge-field configurations using modern visualisation techniques. This time, the manner in which light dynamical fermion degrees of freedom impact the centre-vortex structure is explored. Focusing on the thin vortices identified by plaquettes having a non-trivial centre phase, the vortex structure is illustrated through 3D renderings of oriented plaquettes. The impact of light dynamical fermions is not subtle, changing both the density of vortices and the complexity of the vortex structures observed. The role of vortex branching points in full QCD is highlighted in the survey of results presented.

preprint2022arXiv

Impact of Dynamical Fermions on the Centre Vortex Gluon Propagator

The impact of $SU(3)$ centre vortices on the Landau-gauge gluon propagator is calculated in the presence of dynamical fermions and compared to the pure Yang-Mills case. The presence of dynamical fermions is found to alter the behaviour of the centre vortex propagator when compared to the established pure-gauge result. The gluon spectral representation is also explored from the centre vortex perspective, where centre vortices are shown to exhibit clear signs of positivity violation, which is an indicator of confinement. Vortex removal subsequently restores positivity, demonstrating the crucial role centre vortices play in the confinement of gluons.

preprint2022arXiv

Smoothing algorithms for projected center-vortex gauge fields

We study the application of $\mathrm{SU}(3)$ gauge field smoothing methods to $\mathrm{Z}(3)$-projected center-vortex gauge fields. Due to the proportionality of the vortex links to the identity, naive applications of these methods are either ineffectual or limited in scope, containing subtle issues which are not obviously manifest. To overcome these issues we introduce centrifuge preconditioning, a novel method applied prior to smoothing that rotates the links away from the center while preserving the fundamental structure of the vortex field. Additionally, the concept of vortex-preserved annealed smoothing is formulated to ensure the smoothing procedure maintains the underlying vortex structure. The application of these new methods in the context of annealed smoothing applied to vortex fields is shown to successfully achieve the desired smoothness condition required for the study of more advanced operators.

preprint2022arXiv

Static quark potential from centre vortices in the presence of dynamical fermions

For the first time, centre vortices are identified on SU(3) lattice ensembles that include dynamical fermions. Using a variational method, the static quark potential is calculated on untouched, vortex-removed, and vortex-only fields. Two dynamical ensembles and one pure gauge ensemble are studied, allowing for an exploration of the impact of dynamical fermions on the centre-vortex vacuum. Novel modifications to the standard Coulomb term are introduced to describe the long range behaviour of the vortex-removed potential. These modifications remove a source of systematic error in the fitted string tension on the original ensembles. Our pure Yang-Mills result is consistent with previous studies, where projected centre-vortex fields only reproduce approximately two thirds of the string tension. Remarkably, we find that the vortex-only fields on both dynamical lattices are able to fully reproduce the respective untouched string tensions.

preprint2021arXiv

Towards high partial waves in lattice QCD with a dumbbell-like operator

An extended two-hadron operator is developed to extract the spectra of irreducible representations (irreps) in the finite volume. The irreps of the group for the finite volume system are projected using a coordinate-space operator. The correlation function of this operator is computationally efficient to extract lattice spectra of the specific irrep. In particular, this new formulation only requires propagators to be computed from two distinct source locations, at fixed spatial separation. We perform a proof-of-principle study on a $24^3 \times 48$ lattice volume with $m_π\approx 900$ MeV by isolating various spectra of the $ππ$ system with isospin-2 including a range of total momenta and irreps. By applying the Lüscher formalism, the phase shifts of $S$-, $D$- and $G$-wave $ππ$ scattering with isospin-2 are extracted from the spectra.

preprint2020arXiv

Elastic Form Factors of Nucleon Excitations in Lattice QCD

First principles calculations of the form factors of baryon excitations are now becoming accessible through advances in Lattice QCD techniques. In this paper, we explore the utility of the parity-expanded variational analysis (PEVA) technique in calculating the Sachs electromagnetic form factors for excitations of the proton and neutron. We study the two lowest-lying odd-parity excitations and demonstrate that at heavier quark masses, these states are dominated by behaviour consistent with constituent quark models for the $N^*(1535)$ and $N^*(1650)$, respectively. We also study the lowest-lying localised even-parity excitation, and find that its form factors are consistent with a radial excitation of the ground state nucleon. A comparison of the results from the PEVA technique with those from a conventional variational analysis exposes the necessity of the PEVA approach in baryon excited-state studies.

preprint2020arXiv

Emergent Structure in QCD

The structure of the SU(3) gauge-field vacuum is explored through visualisations of centre vortices and topological charge density. Stereoscopic visualisations highlight interesting features of the vortex vacuum, especially the frequency with which singular points appear and the important connection between branching points and topological charge. This work demonstrates how visualisations of the QCD ground-state fields can reveal new perspectives of centre-vortex structure.

preprint2020arXiv

Magnetic polarisability of the nucleon using a Laplacian mode projection

Conventional hadron interpolating fields, which utilise gauge-covariant Gaussian smearing, are ineffective in isolating ground state nucleons in a uniform background magnetic field. There is evidence that residual Landau mode physics remains at the quark level, even when QCD interactions are present. In this work, quark-level projection operators are constructed from the $SU(3) \times U(1)$ eigenmodes of the two-dimensional lattice Laplacian operator associated with Landau modes. These quark-level modes are formed from a periodic finite lattice where both the background field and strong interactions are present. Using these eigenmodes, quark-propagator projection operators provides the enhanced hadronic energy-eigenstate isolation necessary for calculation of nucleon energy shifts in a magnetic field. The magnetic polarisability of both the proton and neutron is calculated using this method on the $32^3 \times 64$ dynamical QCD lattices provided by the PACS-CS Collaboration. A chiral effective-field theory analysis is used to connect the lattice QCD results to the physical regime, obtaining magnetic polarisabilities of $β^p = 2.79(22)({}^{+13}_{-18}) \times 10^{-4}$ fm$^3$ and $β^n = 2.06(26)({}^{+15}_{-20}) \times 10^{-4}$ fm$^3$, where the numbers in parantheses describe statistical and systematic uncertainties.

preprint2020arXiv

Pion in a uniform background magnetic field with clover fermions

Background field methods provide an important nonperturbative formalism for the determination of hadronic properties which are complementary to matrix-element calculations. However, new challenges are encountered when utilising a fermion action exposed to additive mass renormalisations. In this case, the background field can induce an undesired field-dependent additive mass renormalisation that acts to change the quark mass as the background field is changed. For example, in a calculation utilising Wilson fermions in a uniform background magnetic field, the Wilson term introduced a field-dependent renormalisation to the quark mass which manifests itself in an unphysical increase of the neutral-pion mass for large magnetic fields. Herein, the clover fermion action is studied to determine the extent to which the removal of $\mathcal{O}(a)$ discretisation errors suppresses the field-dependent changes to the quark mass. We illustrate how a careful treatment of nonperturbative improvement is necessary to resolve this artefact of the Wilson term. Using the $32^3 \times 64$ dynamical-fermion lattices provided by the PACS-CS Collaboration we demonstrate how our technique suppresses the unphysical mass renormalisation over a broad range of magnetic field strengths.

preprint2020arXiv

Role of chiral symmetry in the nucleon excitation spectrum

The origin of the low-lying nature of the $N$*(1440), or Roper resonance, has been the subject of significant interest for many years, including several investigations using lattice QCD. The majority of lattice studies have not observed a low-lying excited state energy level in the region of the Roper resonance. However, it has been claimed that chiral symmetry could play an important role in our understanding of this resonance. The purpose of this study is to systematically examine the role of chiral symmetry in the low-lying nucleon spectrum by directly comparing the clover and overlap fermion actions. To ensure any differences in results are attributable to the choice of fermion action, simulations are performed on the same set of gauge field configurations at matched pion masses. Correlation matrix techniques are employed to determine the excitation energy of the first positive-parity excited state for each action. The clover and overlap actions show a remarkable level of agreement. We do not find any evidence that fermion action chiral symmetry plays a significant role in understanding the Roper resonance on the lattice.

preprint2020arXiv

Structure and transitions of nucleon excitations via parity-expanded variational analysis

The recently-introduced Parity Expanded Variational Analysis (PEVA) technique allows for the isolation of baryon eigenstates on the lattice at finite momentum free from opposite-parity contamination. We find that this technique introduces a statistically significant correction in extractions of the electromagnetic form factors of the ground state nucleon. It also allows first extractions of the elastic and transition form factors of nucleon excitations on the lattice. We present the electromagnetic elastic form factors and helicity amplitudes of two odd-parity excitations of the nucleon. These results provide valuable insight into the structure of these states, and allow for a connection to be made to quark-model states in this energy region.

preprint2020arXiv

Visualisation of Centre Vortex Structure

The centre vortex structure of the $SU(3)$ gauge field vacuum is explored through the use of novel visualisation techniques. The lattice is partitioned into 3D time slices, and vortices are identified by locating plaquettes with nontrivial centre phases. Vortices are illustrated by rendering vortex lines that pierce these nontrivial plaquettes. Nontrivial plaquettes with one dimension in the suppressed time direction are rendered by identifying the visible spatial link. These visualisations highlight the frequent presence of singular points and reveal an important role for branching points in $SU(3)$ gauge theory in creating high topological charge density regimes. Visualisations of the topological charge density are presented, and an investigation into the correlation between vortex structures and topological charge density is conducted. The results provide new insight into the mechanisms by which centre vortices generate nontrivial gauge field topology. This work demonstrates the utility of visualisations in conducting centre vortex studies, presenting new avenues with which to investigate this perspective of the QCD vacuum.

preprint2020arXiv

Visualisations of Centre Vortices

The centre vortex structure of the vacuum is visualised through the use of novel 3D visualisation techniques. These visualisations allow for a hands-on examination of the centre-vortex matter present in the QCD vacuum, and highlights some of the key features of the centre-vortex model. The connection between topological charge and singular points is also explored. This work highlights the useful role visualisations play in the exploration of the QCD vacuum.

preprint2016arXiv

Hamiltonian effective field theory study of the $\mathbf{N^*(1535)}$ resonance in lattice QCD

Drawing on experimental data for baryon resonances, Hamiltonian effective field theory (HEFT) is used to predict the positions of the finite-volume energy levels to be observed in lattice QCD simulations of the lowest-lying $J^P=1/2^-$ nucleon excitation. In the initial analysis, the phenomenological parameters of the Hamiltonian model are constrained by experiment and the finite-volume eigenstate energies are a prediction of the model. The agreement between HEFT predictions and lattice QCD results obtained on volumes with spatial lengths of 2 and 3 fm is excellent. These lattice results also admit a more conventional analysis where the low-energy coefficients are constrained by lattice QCD results, enabling a determination of resonance properties from lattice QCD itself. Finally, the role and importance of various components of the Hamiltonian model are examined.

preprint2016arXiv

Nucleon matrix elements using the variational method in lattice QCD

The extraction of hadron matrix elements in lattice QCD using the standard two- and three-point correlator functions demands careful attention to systematic uncertainties. One of the most commonly studied sources of systematic error is contamination from excited states. We apply the variational method to calculate the axial vector current $g_{A}$, the scalar current $g_{S}$ and the quark momentum fraction $\left<x\right>$ of the nucleon and we compare the results to the more commonly used summation and two-exponential fit methods. The results demonstrate that the variational approach offers a more efficient and robust method for the determination of nucleon matrix elements.

preprint2015arXiv

Improved determination of hadron matrix elements using the variational method

The extraction of hadron form factors in lattice QCD using the standard two- and three-point correlator functions has its limitations. One of the most commonly studied sources of systematic error is excited state contamination, which occurs when correlators are contaminated with results from higher energy excitations. We apply the variational method to calculate the axial vector current gA and compare the results to the more commonly used summation and two-exponential fit methods. The results demonstrate that the variational approach offers a more efficient and robust method for the determination of nucleon matrix elements.

preprint2015arXiv

Instanton contributions to the low-lying hadron mass spectrum

The role of instanton-like objects in the QCD vacuum on the mass spectrum of low-lying light hadrons is explored in lattice QCD. Using over-improved stout-link smearing, tuned to preserve instanton-like objects in the QCD vacuum, the evolution of the mass spectrum under smearing is examined. The calculation is performed using a $20^3\times40$ dynamical fat-link-irrelevant-clover (FLIC) fermion action ensemble with lattice spacing 0.126 fm. Through the consideration of a range of pion masses, the effect of the vacuum instanton content is compared at a common pion mass. While the qualitative features of ground-state hadrons are preserved on instanton-dominated configurations, the excitation spectrum experiences significant changes. The underlying physics revealed shows little similarity to the direct-instanton interaction predictions of the instanton liquid model.

preprint2015arXiv

Lattice baryon spectroscopy with multi-particle interpolators

In 2 + 1 flavour lattice QCD the spectrum of the nucleon is presented for both parities using local meson-baryon type interpolating fields in addition to the standard three-quark nucleon interpolators. The role of local five-quark operators in extracting the nucleon excited state spectrum via correlation matrix techniques is explored on dynamical gauge fields with $m_π$ = 293 MeV, leading to the observation of a state consistent with the S-wave $Nπ$ scattering threshold in the negative-parity sector. Furthermore, the robustness of the variational technique is examined by studying the spectrum on a variety of operator bases. Fitting a single-state ansatz to the eigenstate-projected correlators provides robust energies for the low-lying spectrum that are essentially invariant despite being extracted from qualitatively different bases.

preprint2015arXiv

Lattice QCD Evidence that the Lambda(1405) Resonance is an Antikaon-Nucleon Molecule

For almost 50 years the structure of the Lambda(1405) resonance has been a mystery. Even though it contains a heavy strange quark and has odd parity, its mass is lower than any other excited spin-1/2 baryon. Dalitz and co-workers speculated that it might be a molecular state of an antikaon bound to a nucleon. However, a standard quark-model structure is also admissible. Although the intervening years have seen considerable effort, there has been no convincing resolution. Here we present a new lattice QCD simulation showing that the strange magnetic form factor of the Lambda(1405) vanishes, signaling the formation of an antikaon-nucleon molecule. Together with a Hamiltonian effective-field-theory model analysis of the lattice QCD energy levels, this strongly suggests that the structure is dominated by a bound antikaon-nucleon component. This result clarifies that not all states occurring in nature can be described within a simple quark model framework and points to the existence of exotic molecular meson-nucleon bound states.

preprint2015arXiv

Light Meson Form Factors at near Physical Masses

The ability for most hadrons to decay via strong interactions prevents the direct measurement of their electromagnetic properties. However, a detailed understanding of how these resonant states feature in scattering processes can allow one to disentangle such information from photo production processes. In particular, there has been increasing interest in the determination of magnetic dipole moments using such methods. In a recent study, Gudino et al. provide the first experimental determination of the magnetic dipole moment of the rho meson. To facilitate a comparison with this experimental determination, we present a calculation of the rho meson and pion electromagnetic form factors calculated in the framework of Lattice QCD. Using the PACS-CS 2+1 flavour full QCD gauge field configurations, we are able to access low $Q^2$ values at near-physical quark masses. Through the use of variational techniques, we control excited state systematics in the matrix elements of the lowest-lying states and gain access to the matrix elements of the first excited state. Our determination of the rho meson $g$-factor $g_ρ = 2.21(8)$ is in excellent agreement with this experimental determination, but with a significantly smaller uncertainty.

preprint2015arXiv

N* Spectroscopy from Lattice QCD: The Roper Explained

This brief review focuses on the low-lying even- and odd-parity excitations of the nucleon obtained in recent lattice QCD calculations. Commencing with a survey of the 2014-15 literature we'll see that results for the first even-parity excitation energy can differ by as much as 1 GeV, a rather unsatisfactory situation. Following a brief review of the methods used to isolate excitations of the nucleon in lattice QCD, and drawing on recent advances, we'll see how a consensus on the low-lying spectrum has emerged among many different lattice groups. To provide insight into the nature of these states we'll review the wave functions and electromagnetic form factors that are available for a few of these states. Consistent with the Luscher formalism for extracting phase shifts from finite volume spectra, the Hamiltonian approach to effective field theory in finite volume can provide guidance on the manner in which physical quantities manifest themselves in the finite volume of the lattice. With this insight, we will address the question; Have we seen the Roper in lattice QCD?

preprint2015arXiv

Transition of $ρ\rightarrow πγ$ in Lattice QCD

With the ongoing experimental interest in exploring the excited hadron spectrum, evaluations of the matrix elements describing the formation and decay of such states via radiative processes provide us with an important connection between theory and experiment. In particular, determinations obtained via the lattice allow for a direct comparison of QCD-expectation with experimental observation. Here we present the first light quark determination of the $ρ\rightarrow πγ$ transition form factor from lattice QCD using dynamical quarks. Using the PACS-CS 2+1 flavour QCD ensembles we are able to obtain results across a range of masses, to the near physical value of $m_π= 157$ MeV. An important aspect of our approach is the use of variational methods to isolate the desired QCD eigenstate. For low-lying states, such techniques facilitate the removal of excited state contributions. In principle the method enables one to consider arbitrary eigenstates. We find our results are in accord with the non-relativistic quark model for heavy masses. In moving towards the light-quark regime we observe an interesting quark mass dependence, contrary to the quark model expectation. Comparison of our light-quark result with experimental determinations highlights a significant discrepancy suggesting that disconnected sea-quark loop contributions may play a significant role in fully describing this process.

preprint2014arXiv

Electromagnetic matrix elements for negative parity nucleons

Here we present preliminary results for the evaluation of the electromagnetic form factors for the lowest-lying negative-parity, spin-$\frac{1}{2}$ nucleons, namely the $S_{11}(1535)$ and $S_{11}(1650)$, through the use of the variational method. We find that the characteristics of the electric form factor, $G_{E}$, are similar between these states, however significant differences are observed between the quark-sector contributions to the magnetic form factor, $G_{M}$. Within simple constituent quark models, these states are understood to be admixtures of $s=\frac{1}{2}$ and $s=\frac{3}{2}$ states coupled to orbital angular momentum $\ell = 1$. Our results reveal a qualitative difference in the manner in which the singly-represented quark sector contributes to these baryon magnetic form factors.

preprint2014arXiv

Nucleon spectroscopy using multi-particle operators

The role of 5-quark operators in extracting the nucleon excited state spectrum via correlation matrix techniques is explored. In particular, the coupling of meson-baryon operators to nucleon resonance states and scattering states is studied. Results are presented for 2+1 flavour dynamical ensembles in both the positive and negative parity channels. Fitting a single-state ansatz to the eigenstate-projected correlators provides robust energies for the low-lying spectrum that are essentially invariant across a variety of different operator bases. In particular, the resonant state energies obtained both with and without the use of meson-baryon operators agree, demonstrating that resonance energies can be reliably extracted solely using 3-quark operators.

preprint2014arXiv

On the Structure of the Lambda 1405

For almost 50 years the structure of the Lambda 1405 resonance has been a mystery. Recently, a new lattice QCD simulation showing that its strange magnetic form factor vanishes, together with a comprehensive Hamiltonian analysis of the lattice QCD energy levels, has unambiguously established that the structure is dominated by a bound anti-kaon--nucleon component [1]. Here we present supplementary information for Ref. [1] including a presentation of the relevant Hamiltonian effective field theory and an illustration of the volume dependence of the results and their connection to the infinite volume limit of Nature.

preprint2014arXiv

Searching for low-lying multi-particle thresholds in lattice spectroscopy

We explore the Euclidean-time tails of odd-parity nucleon correlation functions in a search for the S-wave pion-nucleon scattering-state threshold contribution. The analysis is performed using 2+1 flavor 32^3 x 64 PACS-CS gauge configurations available via the ILDG. Correlation matrices composed with various levels of fermion source/sink smearing are used to project low-lying states. The consideration of 25,600 fermion propagators reveals the presence of more than one state in what would normally be regarded as an eigenstate-projected correlation function. This observation is in accord with the scenario where the eigenstates contain a strong mixing of single and multi-particle states but only the single particle component has a strong coupling to the interpolating field. Employing a two-exponential fit to the eigenvector-projected correlation function, we are able to confirm the presence of two eigenstates. The lower-lying eigenstate is consistent with a N-pi scattering threshold and has a relatively small coupling to the three-quark interpolating field. We discuss the impact of this small scattering-state contamination in the eigenvector projected correlation function on previous results presented in the literature.

preprint2014arXiv

Visualisations of coherent centre domains in local Polyakov loops

Quantum Chromodynamics exhibits a hadronic confined phase at low to moderate temperatures and, at a critical temperature $T_C$, undergoes a transition to a deconfined phase known as the quark-gluon plasma. The nature of this deconfinement phase transition is probed through visualizations of the Polyakov loop, a gauge independent order parameter. We produce visualizations that provide novel insights into the structure and evolution of center clusters. Using the HMC algorithm the percolation during the deconfinement transition is observed. Using 3D rendering of the phase and magnitude of the Polyakov loop, the fractal structure and correlations are examined. The evolution of the center clusters as the gauge fields thermalize from below the critical temperature to above it are also exposed. We observe deconfinement proceeding through a competition for the dominance of a particular center phase. We use stout-link smearing to remove small-scale noise in order to observe the large-scale evolution of the center clusters. A correlation between the magnitude of the Polyakov loop and the proximity of its phase to one of the center phases of SU(3) is evident in the visualizations.

preprint2013arXiv

Electromagnetic Form Factors for the $Λ$(1405)

Building on our successful technique to isolate the otherwise-elusive $Λ$(1405) using correlation matrix techniques and multiple source and sink smearings, we present calculations of the quark sector contributions to the electric form factors of the $Λ$(1405). Using the PACS-CS $(2+1)$-flavour full-QCD ensembles available through the ILDG, our calculations reveal behaviour consistent with the development of a non-trivial molecular $\overline{K}N$ bound-state component as one approaches the physical values of the $u$ and $d$ quark masses.

preprint2013arXiv

Magnetic properties of the nucleon in a uniform background field

We present results for the magnetic moment and magnetic polarisability of the neutron and the magnetic moment of the proton. These results are calculated using the uniform background field method on 32^3 x 64 dynamical QCD lattices provided by the PACS-CS collaboration as part of the ILDG. We use a uniform background magnetic field quantised by the periodic spatial volume. We investigate ways to improve the effective energy plots used to calculate magnetic polarisabilities, including the use of correlation matrix techniques with various source smearings.

preprint2013arXiv

Multi-Particle Baryon Spectroscopy

In Nature the excited states of the hadron spectrum appear as resonances. Consequently, there has been significant interest in studying the excited baryon spectrum using lattice QCD. With this in mind we perform spectroscopic calculations with five-quark interpolating fields. Stochastic estimation techniques are used in order to calculate the loop propagators, with dilution in spin, colour and time implemented as a means of variance reduction. We present effective mass plots extracted from these five-quark interpolators, and examine the contributions from fully-connected and loop-containing pieces of the correlation function, keeping in mind their use in future corre- lation matrix studies.

preprint2013arXiv

Nucleon Excited State Wave Functions from Lattice QCD

We apply the eigenvectors from a variational analysis to successfully extract the wave functions of even-parity excited states of the nucleon, including the Roper. We explore the first four states in the spectrum excited by the standard nucleon interpolating field. We find that the states exhibit a structure qualitatively consistent with a constituent quark model, where the ground, first-, second- and third-excited states have 0, 1, 2, and 3 nodes in the radial wave function of the d-quark about two $u$ quarks at the origin. Moreover the radial amplitude of the probability distribution is similar to that predicted by constituent quark models. We present a detailed examination of the quark-mass dependence of the probability distributions for these states, searching for a nontrivial role for the multi-particle components mixed in the finite-volume QCD eigenstates. Finally we examine the dependence of the d-quark probability distribution on the positions of the two u quarks. The results are fascinating, with the underlying S-wave orbitals governing the distributions even at rather large u-quark separations.

preprint2013arXiv

Structure and Flow of the Nucleon Eigenstates in Lattice QCD

A determination of the excited energy eigenstates of the nucleon, $s=c{1}{2}$, $I={1}{2}$, $N^{\pm}$, is presented in full QCD using 2+1 flavor PACS-CS gauge configurations. The correlation-matrix method is used and is built using standard nucleon interpolators employing smearings at the fermion sources and sinks. We develop and demonstrate a new technique that allows the eigenvectors obtained to be utilized to track the propagation of the intrinsic nature of energy-states from one quark mass to the next. This approach is particularly useful for larger dimension correlation matrices where more near-degenerate energy-states can appear in the spectrum.

preprint2013arXiv

SU(3) centre vortices underpin confinement and dynamical chiral symmetry breaking

The mass function of the nonperturbative quark propagator in SU(3) gauge theory shows only a weak dependence on the vortex content of the gauge configurations. Of particular note is the survival of dynamical mass generation on vortex-free configurations having a vanishing string tension. This admits the possibility that mass generation associated with dynamical chiral symmetry breaking persists without confinement. In this presentation, we examine the low-lying ground-state hadron spectrum of the pi, rho, N and Delta and discover that while dynamical mass generation persists in the vortex-free theory, it is not connected to dynamical chiral symmetry breaking. In this way, centre vortices in SU(3) gauge theory are intimately linked to both confinement and dynamical chiral symmetry breaking. We conclude that centre vortices are the essential underlying feature of the QCD vacuum.

preprint2013arXiv

Variational Approach to the Calculation of gA

We present a lattice QCD calculation of the nucleon axial charge, gA, using a variational approach. After a brief outline of how the variational method is applied to the calculation of form factors, we present results for gA using this method. We find that ground state dominance is rapid, evident in the early onset of a clear plateau in the correlation function ratio proportional to gA. Through a comparison with results obtained via traditional methods, we show how excited state effects can suppress gA by as much as 8% if sources are not properly tuned.

preprint2013arXiv

Wave Function of the Roper from Lattice QCD

We apply the eigenvectors from a variational analysis in lattice QCD to successfully extract the wave function of the Roper state, and a higher mass P_11 state of the nucleon. We use the 2+1 flavour 32^3x64 PACS-CS configurations at a near physical pion mass of 156 MeV. We find that both states exhibit a structure consistent with a constituent quark model. The Roper d-quark wave function contains a single node consistent with a 2S state, and the third state wave function contains two, consistent with a 3S state. A detailed comparison with constituent quark model wave functions is carried out, obtained from a Coulomb plus ramp potential. These results validate the approach of accessing these states by constructing a variational basis composed of different levels of fermion source and sink smearing. Furthermore, significant finite volume effects are apparent for these excited states which mix with multi-particle states, driving their masses away from physical values and enabling the extraction of resonance parameters from lattice QCD simulations.

preprint2012arXiv

Accessing High Momentum States In Lattice QCD

Two measures are defined to evaluate the coupling strength of smeared interpolating operators to hadronic states at a variety of momenta. Of particular interest is the extent to which strong overlap can be obtained with individual high-momentum states. This is vital to exploring hadronic structure at high momentum transfers on the lattice and addressing interesting phenomena observed experimentally. We consider a novel idea of altering the shape of the smeared operator to match the Lorentz contraction of the probability distribution of the high-momentum state, and show a reduction in the relative error of the two-point function by employing this technique. Our most important finding is that the overlap of the states becomes very sharp in the smearing parameters at high momenta and fine tuning is required to ensure strong overlap with these states.

preprint2012arXiv

GPUs: An Oasis in the Supercomputing Desert

A novel metric is introduced to compare the supercomputing resources available to academic researchers on a national basis. Data from the supercomputing Top 500 and the top 500 universities in the Academic Ranking of World Universities (ARWU) are combined to form the proposed "500/500" score for a given country. Australia scores poorly in the 500/500 metric when compared with other countries with a similar ARWU ranking, an indication that HPC-based researchers in Australia are at a relative disadvantage with respect to their overseas competitors. For HPC problems where single precision is sufficient, commodity GPUs provide a cost-effective means of quenching the computational thirst of otherwise parched Lattice practitioners traversing the Australian supercomputing desert. We explore some of the more difficult terrain in single precision territory, finding that BiCGStab is unreliable in single precision at large lattice sizes. We test the CGNE and CGNR forms of the conjugate gradient method on the normal equations. Both CGNE and a modified form of CGNR (with restarts) provide reliable convergence for quark propagator calculations in single precision.

preprint2012arXiv

Low-lying Odd-parity States of the Nucleon in Lattice QCD

The world's first examination of the odd-parity nucleon spectrum at light quark masses in 2+1 flavor lattice QCD is presented. Configurations generated by the PACS-CS collaboration and made available through the ILDG are used, with the lightest pion mass at 156 MeV. A novel method for tracking the individual energy eigenstates as the quark mass changes is introduced. The success of this approach reveals the flow of the states towards the physical masses. Using the correlation matrix method, the two lowest-energy states revealed are found to be in accord with the physical spectrum of Nature.

preprint2012arXiv

Magnetic properties of the neutron in a uniform background field

We present calculations of the magnetic moment and magnetic polarisability of the neutron from the background field method. The calculations are performed on $32^3\times64$ dynamical lattices generated by the PACS-CS collaboration and made available via the ILDG. We consider uniform fields quantised by the periodic spatial volume. We explore different approaches for improving the quality of the fits used in the results. Also included are initial results for the magnetic moment of the lowest lying negative parity nucleon states.

preprint2012arXiv

Roper Resonance in 2+1 Flavor QCD

The low-lying even-parity states of the nucleon are explored in lattice QCD using the PACS-CS collaboration 2+1-flavor dynamical-QCD gauge-field configurations made available through the International Lattice Datagrid (ILDG). The established correlation-matrix approach is used, in which various fermion source and sink smearings are utilized to provide an effective basis of interpolating fields to span the space of low-lying energy eigenstates. Of particular interest is the nature of the first excited state of the nucleon, the $N{1/2}^{+}$ Roper resonance of $P_{11}$ pion-nucleon scattering. The Roper state of the present analysis approaches the physical mass, displaying significant chiral curvature at the lightest quark mass. These full QCD results, providing the world's first insight into the nucleon mass spectrum in the light-quark regime, are significantly different from those of quenched QCD and provide interesting insights into the dynamics of QCD.

preprint2011arXiv

A Novel Multiple-Time Scale Integrator for the Hybrid Monte Carlo Algorithm

Hybrid Monte Carlo simulations that implement the fermion action using multiple terms are commonly used. By the nature of their formulation they involve multiple integration time scales in the evolution of the system through simulation time. These different scales are usually dealt with by the Sexton-Weingarten nested leapfrog integrator. In this scheme the choice of time scales is somewhat restricted as each time step must be an exact multiple of the next smallest scale in the sequence. A novel generalisation of the nested leapfrog integrator is introduced which allows for far greater flexibility in the choice of time scales, as each scale now must only be an exact multiple of the smallest step size.

preprint2011arXiv

Isolating the Λ(1405) in Lattice QCD

The negative-parity ground state of the Λ baryon lies surprisingly low in mass. At 1405.1 MeV, it lies lower than the negative-parity ground state nucleon, even though it has a valence strange quark. Using the PACS-CS (2+1)-flavour full-QCD ensembles available through the ILDG, we employ a variational analysis using source and sink smearing to isolate this elusive state. We find three low-lying odd-parity states, and for the first time reproduce the correct level ordering with respect to the nearby scattering thresholds.

preprint2011arXiv

Polynomial Filtered HMC -- an algorithm for lattice QCD with dynamical quarks

Polynomial approximations to the inverse of the fermion matrix are used to filter the dynamics of the upper energy scales in HMC simulations. The use of a multiple time-scale integration scheme allows the filtered pseudofermions to be evolved using a coarse step size. We introduce a novel generalisation of the nested leapfrog which allows for far greater flexibility in the choice of time scales. We observe a reduction in the computational expense of the molecular dynamics integration of between 3--5 which improves as the quark mass decreases.

preprint2011arXiv

Wave Functions of the Proton Ground State in the Presence of a Uniform Background Magnetic Field in Lattice QCD

We calculate the probability distributions of quarks in the ground state of the proton, and how they are affected in the presence of a constant background magnetic field. We focus on wave functions in the Landau and Coulomb gauges. We observe the formation of a scalar u-d diquark clustering. The overall distortion of the quark probability distribution under a very large magnetic field, as demanded by the quantisation conditions on the field, is quite small. The effect is to elongate the distributions along the external field axis while localizing the remainder of the distribution.

preprint2010arXiv

Ordering of Spin-$\frac{1}{2}$ Excitations of the Nucleon in Lattice QCD

We present results for the negative parity low-lying state of the nucleon, $N{\frac{1}{2}}^{-}$ (1535 MeV) ${\rm S}_{11}$, from a variational analysis method. The analysis is performed in quenched QCD with the FLIC fermion action. The principal focus of this paper is to explore the level ordering between the Roper (${\rm{P}}_{11}$) and the negative parity ground (${\rm{S}}_{11}$) states of the nucleon. Evidence of the physical level ordering is observed at light quark masses. A wide variety of smeared-smeared correlation functions are used to construct correlation matrices. A comprehensive correlation matrix analysis is performed to ensure an accurate isolation of the $N{\frac{1}{2}}^{-}$ state.

preprint2010arXiv

Positive-parity Excited-states of the Nucleon in Quenched Lattice QCD

Positive-parity spin-1/2 excitations of the nucleon are explored in lattice QCD. The variational method is used in this investigation and several correlation matrices are employed. As our focus is on the utility and methodology of the variational approach, we work in the quenched approximation to QCD. Various sweeps of Gaussian fermion-field smearing is applied at the source and at the sink of $χ_{1}\barχ_{1}$ and $χ_{1}χ_{2}$ correlation functions to obtain a large basis of operators. Using several different approaches for constructing basis interpolators, we demonstrate how improving the basis can split what otherwise might be interpreted as a single state into multiple eigenstates. Consistency of the extracted excited energy states are explored over various dimensions of the correlation matrices. The use of large correlation matrices is emphasized for the reliable extraction of the excited eigenstates of QCD.

preprint2009arXiv

Isolating Excited States of the Nucleon in Lattice QCD

We discuss a robust projection method for the extraction of excited-state masses of the nucleon from a matrix of correlation functions. To illustrate the algorithm in practice, we present results for the positive parity excited states of the nucleon in quenched QCD. Using eigenvectors obtained via the variational method, we construct an eigenstate-projected correlation function amenable to standard analysis techniques. The method displays its utility when comparing results from the fit of the projected correlation function with those obtained from the eigenvalues of the variational method. Standard nucleon interpolators are considered, with $2\times 2$ and $3\times 3$ correlation matrix analyses presented using various combinations of source-smeared, sink-smeared and smeared-smeared correlation functions. Using these new robust methods, we observe a systematic dependency of the nucleon excited-state masses on source- and sink-smearing levels. To the best of our knowledge, this is the first clear indication that a correlation matrix of standard nucleon interpolators is insufficient to isolate the eigenstates of QCD.

preprint2009arXiv

Low-lying positive-parity excited states of the nucleon

We present an overview of the correlation-matrix methods developed recently by the CSSM Lattice Collaboration for the isolation of excited states of the nucleon. Of particular interest is the first positive-parity excited-state of the nucleon known as the Roper resonance. Using eigenvectors of the correlation matrix we construct parity and eigenstate projected correlation functions which are analysed using standardized methods. The robust nature of this approach for extracting the eigenstate energies is presented. We report the importance of using a variety of source and sink smearings in achieving this. Ultimately the independence of the eigenstate energies from the interpolator basis is demonstrated. In particular we consider $4\times 4$ correlation matrices built from a variety of interpolators and smearing levels. Using FLIC fermions to access the light quark mass regime, we explore the curvature encountered in the energy of the states as the chiral limit is approached. We report a low-lying Roper state contrasting earlier results using correlation matrices. To the best of our knowledge, this is the first time a low-lying Roper resonance has been found using correlation matrix methods. Finally, we present our results in the context of the Roper results reported by other groups.

preprint2008arXiv

Preconditioning Maximal Center Gauge with Stout Link Smearing in SU(3)

Center vortices are studied in SU(3) gauge theory using Maximal Center Gauge (MCG) fixing. Stout link smearing and over-improved stout link smearing are used to construct a preconditioning gauge field transformation, applied to the original gauge field before fixing to MCG. We find that preconditioning successfully achieves higher gauge fixing maxima. We observe a reduction in the number of identified vortices when preconditioning is used, and also a reduction in the vortex-only string tension.

preprint2004arXiv

Dynamical fat link fermions

The use of APE smearing or other blocking techniques in fermion actions can provide many advantages. There are many variants of these fat link actions in lattice QCD currently, such as FLIC fermions. Frequently, fat link actions make use of the APE blocking technique in combination with a projection of the blocked links back into the special unitary group. This reunitarisation is often performed using an iterative maximisation of a gauge invariant measure. This technique is not differentiable with respect to the gauge field and thus prevents the use of standard Hybrid Monte Carlo simulation algorithms. The use of an alternative projection technique circumvents this difficulty and allows the simulation of dynamical fat link fermions with standard HMC and its variants.

preprint2004arXiv

Hybrid Monte Carlo with Fat Link Fermion Actions

The use of APE smearing or other blocking techniques in lattice fermion actions can provide many advantages. There are many variants of these fat link actions in lattice QCD currently, such as FLIC fermions. The FLIC fermion formalism makes use of the APE blocking technique in combination with a projection of the blocked links back into the special unitary group. This reunitarisation is often performed using an iterative maximisation of a gauge invariant measure. This technique is not differentiable with respect to the gauge field and thus prevents the use of standard Hybrid Monte Carlo simulation algorithms. The use of an alternative projection technique circumvents this difficulty and allows the simulation of dynamical fat link fermions with standard HMC and its variants. The necessary equations of motion for FLIC fermions are derived, and some initial simulation results are presented. The technique is more general however, and is straightforwardly applicable to other smearing techniques or fat link actions.

preprint2004arXiv

Improved chiral properties of FLIC fermions

The chiral properties of the fat-link irrelevant clover (FLIC) fermion action are examined. The improved chiral properties of fermion actions incorporating smoothed links are realized in the FLIC action where only the irrelevant operators of the fermion action are constructed with smoothed links. In particular, the histogram of the additive mass renormalization encountered in chiral-symmetry breaking Wilson-type fermion actions is seen to narrow upon introducing fat-links in the irrelevant operators. The exceptional configuration problem of quenched QCD is reduced, enabling access to the light quark mass regime of $m_π / m_ρ \sim 0.35$. In particular, quenched chiral non-analytic behavior is revealed in the light quark mass dependence of the $Δ$-baryon mass. FLIC fermions offer a promising approach to revealing the properties of full QCD at light quark masses.

preprint2004arXiv

The FLIC Overlap Quark Propagator

FLIC overlap fermions are a variant of the standard (Wilson) overlap action, with the FLIC (Fat Link Irrelevant Clover) action as the overlap kernel rather than the Wilson action. The structure of the FLIC overlap fermion propagator in momentum space is studied, and a comparison against previous studies of the Wilson overlap propagator in quenched QCD is performed. To explore the scaling properties of the propagator for the two actions, numerical calculations are performed in Landau Gauge across three lattices with different lattice spacing $a$ and similar physical volumes. We find that at light quark masses the acti ons agree in both the infrared and the ultraviolet, but at heavier masses some disagreement in the ultraviolet appears. This is attributed to the two action s having different discretisation errors with the FLIC overlap providing superior performance in this regime. Both actions scale reasonably, but some scaling violations are observed.

preprint2003arXiv

Low-lying eigenmodes of the Wilson-Dirac operator and correlations with topological objects

The probability density of low-lying eigenvectors of the hermitian Wilson-Dirac operator is examined. Comparisons in position and size between eigenvectors, topological charge and action density are made. We do this for standard Monte-Carlo generated SU(3) background fields and for single instanton background fields. Both hot and cooled SU(3) background fields are considered. An instanton model is fitted to eigenmodes and topological charge density and the sizes and positions of these are compared.

preprint2002arXiv

Generalised Spin Projection for Fermion Actions

The majority of compute time doing lattice QCD is spent inverting the fermion matrix. The time that this takes increases with the condition number of the matrix. The FLIC(Fat Link Irrelevant Clover) action displays, among other properties, an improved condition number compared to standard actions and hence is of interest due to potential compute time savings. However, due to its two different link sets there is a factor of two cost in floating point multiplications compared to the Wilson action. An additional factor of two has been attributed due to the loss of the so-called spin projection trick. We show that any split-link action may be written in terms of spin projectors, reducing the additional cost to at most a factor of two. Also, we review an efficient means of evaluating the clover term, which is additional expense not present in the Wilson action.