Source author record

Anthony G. Williams

Anthony G. Williams 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

61works
10topics
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

61 published item(s)

preprint2022arXiv

Global fit of 2HDM with future collider results

In this work, we summarize a global fit study of Type-II two Higgs doublet models (2HDM), and explore the impact of future SM-like Higgs and Z-pole precision measurements on the allowed parameter space. The work is based on the study results of a global fit of 2HDMs with the tool GAMBIT, utilising various current constraints including theoretical constraints (unitarity, perturbativity and vacuum stability), Higgs searches at colliders, electroweak physics and flavour constraints. We further investigate the ability of future facilities, such as the HL-LHC, CEPC, ILC and FCC-ee to explore the 2HDM parameter space.

preprint2021arXiv

Convergent Bayesian Global Fits of 4D Composite Higgs Models

Models in which the Higgs boson is a composite pseudo-Nambu-Goldstone boson offer attractive solutions to the Higgs mass naturalness problem. We consider three such models based on the minimal $SO(5) \rightarrow SO(4)$ symmetry breaking pattern, and perform convergent global fits on the models under a Bayesian framework in order to find the regions of their parameter spaces that best fit a wide range of constraints, including recent Higgs measurements. We use a novel technique to analyse the fine-tuning of the models, quantifying the tuning as the Kullback-Leibler divergence from the prior to the posterior probability on the parameter space. Each model is found to be able to satisfy all constraints at the $3σ$ level simultaneously. As a by-product of the fits, we analyse the collider phenomenology of our models in these viable regions. In two of the three models, we find that the $g g \rightarrow H \rightarrow γγ$ cross section is less than ${\sim}90$% that predicted by the SM, which is already in slight tension with experiment and could potentially be ruled out in the future high-luminosity run of the LHC. In addition, the lightest fermions $F$ arising from the new strong dynamics in these models are seen in general to lie above ${\sim}1.1$ TeV, with the $F \rightarrow tW^{+}$ and $F \rightarrow \bar{b}W^{+}$ decays offering particularly promising channels for probing these models in future collider searches.

preprint2020arXiv

Exploring the 2HDM with Global Fits in $\texttt{GAMBIT}$

In this work, we present preliminary results of a global fit of the type-II two-Higgs-doublet model (2HDM) with the tool $\texttt{GAMBIT}$. Our study includes various constraints, including the theoretical constraints (unitarity, perturbativity and vacuum stability), Higgs searches at colliders, electroweak physics and flavour constraints. With the latest experimental results, our results not only confirm past studies but also go further in probing the model. We find, for example, that the measurements of $B\rightarrow K^*μ^+μ^-$ angular observables cannot be explained in the type-II 2HDM.

preprint2020arXiv

Global analyses of Higgs portal singlet dark matter models using GAMBIT

We present global analyses of effective Higgs portal dark matter models in the frequentist and Bayesian statistical frameworks. Complementing earlier studies of the scalar Higgs portal, we use GAMBIT to determine the preferred mass and coupling ranges for models with vector, Majorana and Dirac fermion dark matter. We also assess the relative plausibility of all four models using Bayesian model comparison. Our analysis includes up-to-date likelihood functions for the dark matter relic density, invisible Higgs decays, and direct and indirect searches for weakly-interacting dark matter including the latest XENON1T data. We also account for important uncertainties arising from the local density and velocity distribution of dark matter, nuclear matrix elements relevant to direct detection, and Standard Model masses and couplings. In all Higgs portal models, we find parameter regions that can explain all of dark matter and give a good fit to all data. The case of vector dark matter requires the most tuning and is therefore slightly disfavoured from a Bayesian point of view. In the case of fermionic dark matter, we find a strong preference for including a CP-violating phase that allows suppression of constraints from direct detection experiments, with odds in favour of CP violation of the order of 100:1. Finally, we present DDCalc 2.0.0, a tool for calculating direct detection observables and likelihoods for arbitrary non-relativistic effective operators.

preprint2020arXiv

Minimal 4D Composite Higgs Models Under Current LHC Constraints

We present preliminary results of the first convergent global fits of several minimal composite Higgs models. Our fits are performed using the differential evolution optimisation package $\texttt{Diver}$. A variety of physical constraints are taken into account, including a wide range of exclusion bounds on heavy resonance production from Run 2 of the LHC. As a by-product of the fits, we analyse the collider phenomenology of the lightest new up-type and down-type resonances in the viable regions of our models, finding some low-mass resonances that can be probed in future collider searches.

preprint2020arXiv

Model-independent approach for incorporating interference effects in collider searches for new resonances

The presence of large-mass resonances in the data collected at the Large Hadron Collider would provide direct evidence of physics beyond the Standard Model. A key challenge in current resonance searches at the LHC is the modelling of signal--background interference effects, which can severely distort the shape of the reconstructed invariant mass distribution relative to the case where there is no interference. Such effects are strongly dependent on the beyond the Standard Model theory that must be considered as unknown if one aims to minimise any theoretical bias on the search results. In this paper, we describe a procedure which employs a physically-motivated, model-independent template functional form that can be used to model interference effects, both for the characterisation of positive discoveries, and in the presentation of null results. We illustrate the approach with the example of a scalar resonance decaying into a pair of photons.

preprint2019arXiv

Higgs Quark Flavor Violation: Simplified Models and Status of General Two-Higgs-Doublet Model

We study quark flavor violating interactions mediated by the Higgs boson $h$. We consider observables involving a third generation quark, of both the up and the down quark sectors, like $h\rightarrow bs$ and $t\rightarrow ch$. Using an effective field theory approach we systematically list all the possible tree-level ultraviolet completions, which comprise models with vector-like quarks and/or extra scalars. We provide upper bounds on the flavor violating transitions allowed by current limits stemming from low energy processes, such as meson mixing and $b \rightarrow s γ$. We find that scenarios with vector-like quarks always have very suppressed flavor-violating transitions, while a general two Higgs doublet model may have a sizeable rate. To study the latter case in detail, we perform a full numerical simulation taking into account all relevant theoretical and phenomenological constraints. Our results show that ${\rm BR}(t\rightarrow ch)$ [${\rm BR}(h\rightarrow bs)$] are still allowed at the subpercent [percent] level, which are being [may be] explored at the LHC [future colliders]. Finally, we have found that the mild mass-splitting discrepancy with respect to the SM in the $B_s$ meson system can be accommodated in the Two-Higgs-Doublet Model. If confirmed, it yields the prediction ${\rm BR}(h\rightarrow bs)\simeq 10^{-4}$, if the new contribution to the mass-splitting is dominated by tree-level Higgs boson exchange.

preprint2016arXiv

Boosting the charged Higgs search prospects using jet substructure at the LHC

Charged Higgs bosons are predicted in variety of theoretically well-motivated new physics models with extended Higgs sectors. In this study, we focus on a type-II two Higgs doublet model (2HDM-II) and consider a heavy charged Higgs with its mass ranging from 500 GeV to 1 TeV as dictated by the $b\to sγ$ constraints which render $M_{H^\pm}>480$ GeV. We study the dominant production mode $H^\pm t$ associated production with $H^\pm \to W^\pm A$ being the dominant decay channel when the pseudoscalar $A$ is considerably lighter. For such a heavy charged Higgs, both the decay products $W^\pm$ and $A$ are relatively boosted. In such a scenario, we apply the jet substructure analysis of tagging the fat pseudoscalar and $W$ jets in order to eliminate the standard model background efficiently. We perform a detailed detector simulation for the signal and background processes at the 14 TeV LHC. We introduce various kinematical cuts to determine the signal significance for a number of benchmark points with charged Higgs boson mass from 500 GeV to 1 TeV in the $W^\pm A$ decay channel. Finally we perform a multivariate analysis utilizing a boosted decision tree algorithm to optimize these significances.

preprint2016arXiv

Combined analysis of effective Higgs portal dark matter models

We combine and extend the analyses of effective scalar, vector, Majorana and Dirac fermion Higgs portal models of dark matter (DM), in which DM couples to the Standard Model (SM) Higgs boson via an operator of the form $\mathcal{O}_{\textrm{DM}}\, H^\dagger H$. For the fermion models, we take an admixture of scalar $\overlineψ ψ$ and pseudoscalar $\overlineψ iγ_5 ψ$ interaction terms. For each model, we apply constraints on the parameter space based on the Planck measured DM relic density and the LHC limits on the Higgs invisible branching ratio. For the first time, we perform a consistent study of the indirect detection prospects for these models based on the WMAP7/Planck observations of the cosmic microwave background, a combined analysis of 15 dwarf spheroidal galaxies by Fermi-LAT and the upcoming Cherenkov Telescope Array (CTA). We also perform a correct treatment of the momentum-dependent direct search cross section that arises from the pseudoscalar interaction term in the fermionic DM theories. We find, in line with previous studies, that current and future direct search experiments such as LUX and XENON1T can exclude much of the parameter space, and we demonstrate that a joint observation in both indirect and direct searches is possible for high mass weakly interacting massive particles. In the case of a pure pseudoscalar interaction of a fermionic DM candidate, future gamma-ray searches are the only class of experiment capable of probing the high mass range of the theory.

preprint2016arXiv

Gamma rays from muons from WIMPs: Implementation of radiative muon decays for dark matter analyses

Dark matter searches in gamma ray final states often make use of the fact that photons can be produced from final state muons. Modern Monte Carlo generators and DM codes include the effects of final state radiation from muons produced in the dark matter annihilation process itself, but neglect the O(1%) radiative correction that arises from the subsequent muon decay. After implementing this correction we demonstrate the effect that it can have on dark matter phenomenology by considering the case of dark matter annihilation to four muons via scalar mediator production. We first show that the AMS-02 positron excess can no longer easily be made consistent with this final state once the Fermi-LAT dwarf limits are calculated with the inclusion of radiative muon decays, and we next show that the Fermi-LAT galactic centre gamma excess can be improved with this final state after inclusion of the same effect. We provide code and tables for the implementation of this effect in the popular dark matter code micrOMEGAs, providing a solution for any model producing final state muons.

preprint2016arXiv

Investigating light NMSSM pseudoscalar states with boosted ditau tagging

We study a class of realizations of the Next-to-Minimal Supersymmetric Standard Model that is motivated by dark matter and Higgs data, and in which the lightest pseudoscalar Higgs boson mass is smaller than twice the bottom quark mass and greater than twice the tau lepton mass. In such scenarios, the lightest pseudoscalar Higgs boson can be copiously produced at the LHC from the decay of heavier superpartners and will dominantly further decay into a pair of tau leptons that is generally boosted. We make use of a boosted object tagging technique designed to tag such a ditau jet, and estimate the sensitivity of the LHC to the considered supersymmetric scenarios with 20 to 50~fb$^{-1}$ of proton-proton collisions at a center-of-mass energy of 13~TeV.

preprint2016arXiv

Tau reconstruction methods at an electron-positron collider in the search for new physics

By exploiting the relatively long lifetime of the tau lepton, we propose several novel methods for searching for new physics at an electron-positron collider. We consider processes that involve final states consisting of a tau lepton pair plus two missing particles. The mass and spin of the new physics particles can be measured in 3-prong tau decays. The tau polarization, which reflects the coupling to new physics, can be measured from the $τ\to πν$ decay channel using the impact parameter distribution of the charged pion. We also discuss the corresponding backgrounds for these measurements, the next-to-leading order (NLO) effects, and the implications of finite detector resolution.

preprint2015arXiv

LUX likelihood and limits on spin-independent and spin-dependent WIMP couplings with LUXCalc

We present LUXCalc, a new utility for calculating likelihoods and deriving WIMP-nucleon coupling limits from the recent results of the LUX direct search dark matter experiment. After a brief review of WIMP-nucleon scattering, we derive LUX limits on the spin-dependent WIMP-nucleon couplings over a broad range of WIMP masses, under standard assumptions on the relevant astrophysical parameters. We find that, under these and other common assumptions, LUX excludes the entire spin-dependent parameter space consistent with a dark matter interpretation of DAMA's anomalous signal, the first time a single experiment has been able to do so. We also revisit the case of spin-independent couplings, and demonstrate good agreement between our results and the published LUX results. Finally, we derive constraints on the parameters of an effective dark matter theory in which a spin-1 mediator interacts with a fermionic WIMP and Standard Model fermions via axial-vector couplings. A detailed appendix describes the use of LUXCalc with standard codes to place constraints on generic dark matter theories.

preprint2015arXiv

NMSSM explanations of the Galactic center gamma ray excess and promising LHC searches

The Galactic Center Excess (GCE) is explained in the framework of the Next-to-Minimal Supersymmetric Standard Model (NMSSM) with a $Z_3$ discrete symmetry. We show that a resonant CP-odd Higgs boson with mass twice that of the Dark Matter (DM) candidate is favoured. Meanwhile, the DM candidate is required to have relatively large coupling with the $Z$ boson through its Higgsino component in order to obtain correct DM relic density. Its LHC discovery potential via four signatures is discussed in detail. We find that the most sensitive signals are provided by the Higgsino-like chargino and neutralino pair production with their subsequent decays into $W$ bosons, $Z$ bosons, and DM. The majority of the relevant parameter space can be probed at the Large Hadron Collider with a centre-of-mass energy of 14 TeV and an integrated luminosity 1000 fb$^{-1}$.

preprint2015arXiv

Spin glass behavior of the antiferromagnetic Heisenberg model on scale free network

Randomness and frustration are considered to be the key ingredients for the existence of spin glass (SG) phase. In a canonical system, these ingredients are realized by the random mixture of ferromagnetic (FM) and antiferromagnetic (AF) couplings. The study by Bartolozzi {\it et al.} [Phys. Rev. B{\bf 73}, 224419 (2006)] who observed the presence of SG phase on the AF Ising model on scale free network (SFN) is stimulating. It is a new type of SG system where randomness and frustration are not caused by the presence of FM and AF couplings. To further elaborate this type of system, here we study Heisenberg model on AF SFN and search for the SG phase. The canonical SG Heisenberg model is not observed in $d$-dimensional regular lattices for ($d \leq 3$). We can make an analogy for the connectivity density ($m$) of SFN with the dimensionality of the regular lattice. It should be plausible to find the critical value of $m$ for the existence of SG behaviour, analogous to the lower critical dimension ($d_l$) for the canonical SG systems. Here we study system with $m=2,3,4$ and $5$. We used Replica Exchange algorithm of Monte Carlo Method and calculated the SG order parameter. We observed SG phase for each value of $m$ and estimated its corersponding critical temperature.

preprint2015arXiv

Two-color QCD at high density

QCD at high chemical potential has interesting properties such as deconfinement of quarks. Two-color QCD, which enables numerical simulations on the lattice, constitutes a laboratory to study QCD at high chemical potential. Among the interesting properties of two-color QCD at high density is the diquark condensation, for which we present recent results obtained on a finer lattice compared to previous studies. The quark propagator in two-color QCD at non-zero chemical potential is referred to as the Gor'kov propagator. We express the Gor'kov propagator in terms of form factors and present recent lattice simulation results.

preprint2014arXiv

Dynamical mass generation in unquenched QED using the Dyson--Schwinger equations

We present a comprehensive numerical study of dynamical mass generation for unquenched QED in four dimensions, in the absence of four-fermion interactions, using the Dyson--Schwinger approach. We begin with an overview of previous investigations of criticality in the quenched approximation. To this we add an analysis using a new fermion-antifermion-boson interaction ansatz, the Kizilersu-Pennington (KP) vertex, developed for an unquenched treatment. After surveying criticality in previous unquenched studies, we investigate the performance of the KP vertex in dynamical mass generation using a renormalized fully unquenched system of equations. This we compare with the results for two hybrid vertices incorporating the Curtis--Pennington vertex in the fermion equation. We conclude that the KP vertex is as yet incomplete, and its relative gauge-variance is due to its lack of massive transverse components in its design.

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.

preprint2013arXiv

Strongly-Coupled Unquenched QED4 Propagators Using Schwinger-Dyson Equations

We study unquenched QED in four dimensions using renormalised Schwinger-Dyson equations and focus on the behaviour of the fermion and photon propagators. For this purpose we use an improved Kizilersu-Pennington (KP) vertex which respects gauge invariance, multiplicative renormalizability for the massless case, agrees with perturbation theory in the weak coupling regime and is free of kinematic singularities. We find that the KP vertex performs very well as expected specially in comparison with other vertex choices. We find that the Landau pole problem familiar from perturbative QED persists in the nonperturbative case with the renormalised inverse photon propagator having zero crossing.

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.

preprint2012arXiv

Bosonic stringlike behavior and the Ultraviolet filtering of QCD

The gluonic action density is calculated in static mesons at finite temperature just below the deconfinement point. Our focus is to elucidate the role of vacuum ultraviolet fluctuations which are filtered using an improved smearing algorithm. In the intermediate source separation distance, where the free string picture poorly describes the flux tube width profile, we find upon reducing the vacuum action towards the classical instanton vacuum, the characteristics of the flux tube converge and compare favorably with the predictions of the free bosonic string. This result establishes a connection between the free string action and vacuum gauge fields and reveals the important role of ultraviolet physics in understanding the lattice data at this temperature scale. As a by-product of these calculations, we find the broadening of the QCD flux tube to be independent of the ultraviolet filtering at large distances. Our results exhibit a linearly divergent pattern in agreement with the string picture predictions.

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

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

Gluonic profile of the static baryon at finite temperature

The gluon flux distribution of a static three quark system has been revealed at finite temperature in the pure SU(3) Yang-Mills theory. An action density operator is correlated with three Polyakov loops representing the baryonic state at a temperatures near the end of the QCD plateau, T/T_c = 0.8, and another just before the deconfinement point, T/T_c = 0.9. The flux distributions at short distance separations between the quarks display an action-density profile consistent with a filled Delta shape iso-surface. However the Delta-shaped action iso-surface distributions are found to persist even at large inter-quark separations for both temperatures. The action density distribution in the quark plane exhibits a nonuniform pattern for all quark separations considered. This result contrasts the well-known Y-shaped uniform action density gluonic-flux profile obtained using the Wilson-loop as a quark source operator at zero temperature. We systematically measure and compare the main aspects of the profile of the flux distribution at the two considered temperature scales for three sets of isosceles triangle quark configurations. The radii, amplitudes and rate of change of the width of the flux distribution are found to reverse their behavior as the temperature increase from the end of the QCD plateau region towards the deconfinement point. Remarkably, we find the mean square width of the flux distribution shrinks and localizes for quark separations larger than 1.0 fm at T/T_c = 0.8 which results in an identifiable Y-shaped radius profile. Near the deconfinement point, the action-density delocalizes and the width broadens linearly at large quark separations.

preprint2011arXiv

On the ground state of Yang-Mills theory

We investigate the overlap of the ground state meson potential with sets of mesonic-trial wave functions corresponding to different gluonic distributions. We probe the transverse structure of the flux tube through the creation of non-uniform smearing profiles for the string of glue connecting two color sources in Wilson loop operator. The non-uniformly UV-regulated flux-tube operators are found to optimize the overlap with the ground state and display interesting features in the ground state overlap.

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.

preprint2010arXiv

Role of center vortices in chiral symmetry breaking in SU(3) gauge theory

We study the behavior of the AsqTad quark propagator in Landau gauge on SU(3) Yang-Mills gauge configurations under the removal of center vortices. In SU(2) gauge theory, center vortices have been observed to generate chiral symmetry breaking and dominate the infrared behavior of the quark propagator. In contrast, we report a weak dependence on the vortex content of the gauge configurations, including the survival of dynamical mass generation on configurations with vanishing string tension.

preprint2009arXiv

Impact of stout-link smearing in lattice fermion actions

The impact of stout-link smearing in lattice fermion actions is examined through the consideration of the mass and renormalization functions of the overlap quark propagator over a variety of smeared configurations. Up to six sweeps of stout-link smearing are investigated. For heavy quark masses, the quark propagator is strongly affected by the smearing procedure. For moderate masses, the effect appears to be negligible. A small effect is seen for light quark masses, where dynamical mass generation is suppressed through the smearing procedure.

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.

preprint2009arXiv

Stout-link smearing in lattice fermion actions

The properties of the momentum space quark propagator in Landau gauge are studied for the overlap quark action in quenched lattice QCD. Numerical calculations are performed over four ensembles of gauge configurations, where three are smeared using either 1, 3, or 6 sweeps of stout-link smearing. We calculate the non-perturbative wave function renormalization function $Z(p)$ and the non-perturbative mass function $M(p)$ for a variety of bare quark masses. We find that the wave-function renormalization function is slightly sensitive to the number of stout-link smearing sweeps. For the mass function we find the effect of the stout-link smearing algorithm to be small for moderate to light bare quark masses. For a heavy bare quark mass we find a strong dependence on the number of smearing sweeps.

preprint2007arXiv

Modified Lattice Landau Gauge

We propose a modified lattice Landau gauge based on stereographically projecting the link variables on the circle S^1 -> R for compact U(1) or the 3-sphere S^3 -> R^3 for SU(2) before imposing the Landau gauge condition. This can reduce the number of Gribov copies exponentially and solves the Gribov problem in compact U(1) where it is a lattice artifact. Applied to the maximal Abelian subgroup this might be just enough to avoid the perfect cancellation amongst the Gribov copies in a lattice BRST formulation for SU(N), and thus to avoid the Neuberger 0/0 problem. The continuum limit of the Landau gauge remains unchanged.

preprint2006arXiv

Scaling behavior of quark propagator in full QCD

We study the scaling behavior of the quark propagator on two lattices with similar physical volume in Landau gauge with 2+1 flavors of dynamical quarks in order to test whether we are close to the continuum limit for these lattices. We use configurations generated with an improved staggered (``Asqtad'') action by the MILC collaboration. The calculations are performed on $28^3\times 96$ lattices with lattice spacing $a = 0.09$ fm and on $20^3\times 64$ lattices with lattice spacing $a = 0.12$ fm. We calculate the quark mass function, $M(q^2)$, and the wave-function renormalization function, $Z(q^2)$, for a variety of bare quark masses. Comparing the behavior of these functions on the two sets of lattices we find that both $Z(q^2)$ and $M(q^2)$ show little sensitivity to the ultraviolet cutoff.

preprint2005arXiv

Quark propagator in Landau and Laplacian gauges with overlap fermions

The properties of the momentum space quark propagator in Landau gauge and Gribov copy free Laplacian gauge are studied for the overlap quark action in quenched lattice QCD. Numerical calculations are done on two lattices with different lattice spacing $a$ and the same physical volume. We have calculated the nonperturbative wave function renormalization function $Z(q)$ and the nonperturbative mass function $M(p)$ for a variety of bare quark masses and perform a simple linear extrapolation to the chiral limit. We focus on the comparison of the behavior of $Z(q)$ and $M(p)$ in the chiral limit in the two gauge fixing schemes as well as the behavior on two lattices with different lattice spacing $a$. We find that the mass functions $M(p)$ are very similar for the two gauges while the wave-function renormalization function $Z(q)$ is more strongly infrared suppressed in the Laplacian gauge than in the Landau gauge on the finer lattice. For Laplacian gauge, it seems that the finite $a$ error is large on the coarse lattice which has a lattice spacing $a$ of about 0.124 fm.

preprint2005arXiv

Unquenched quark propagator in Landau gauge

We present an unquenched calculation of the quark propagator in Landau gauge with 2+1 flavors of dynamical quarks. We use configurations generated with an improved staggered (``Asqtad'') action by the MILC collaboration. This quark action has been seen to have excellent rotational symmetry and scaling properties in the quenched quark propagator. Quenched and dynamical calculations are performed on a $20^3\times 64$ lattice with a nominal lattice spacing of $a = 0.125$ fm. The matched quenched and dynamical lattices allow us to investigate the relatively subtle sea quark effects, and even in the quenched case the physical volume of these lattices gives access to lower momenta than our previous study. We calculate the quark mass function and renormalization function for a variety of valence and sea quark masses.

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

Infrared and ultraviolet properties of the Landau gauge quark propagator

We present a current summary of a program to study the quark propagator using lattice QCD. We use the Overlap and ``Asqtad'' quark actions on a number of lattice ensembles to assess systematic errors. We comment on the place of this work amongst studies of QCD Green's functions in other formulations. A preliminary calculation of the running quark mass is presented.

preprint2004arXiv

Quark-gluon vertex in arbitrary kinematics

We compute the quark-gluon vertex in quenched lattice QCD, in the Landau gauge using an off-shell mean-field O(a)-improved fermion action. The complete vertex is computed in two specific kinematical limits, while the Dirac-vector part is computed for arbitrary kinematics. We find a nontrivial and rich tensor structure, including a substantial infrared enhancement of the interaction strength regardless of kinematics.

preprint2004arXiv

Scaling behavior of the overlap quark propagator in Landau gauge

The properties of the momentum space quark propagator in Landau gauge are examined for the overlap quark action in quenched lattice QCD. Numerical calculations are done on three lattices with different lattice spacings and similar physical volumes to explore the approach of the quark propagator toward the continuum limit. We have calculated the nonperturbative momentum-dependent wave function renormalization function Z(p) and the nonperturbative mass function M(p) for a variety of bare quark masses and perform an extrapolation to the chiral limit. We find the behavior of Z(p) and M(p) are in reasonable agreement between the two finer lattices in the chiral limit, however the data suggest that an even finer lattice is desirable. The large momentum behavior is examined to determine the quark condensate.

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.

preprint2004arXiv

Topological Charge Evolution in the Markov-Chain of QCD

The topological charge is studied on lattices of large physical volume and fine lattice spacing. We illustrate how a parity transformation on the SU(3) link-variables of lattice gauge configurations reverses the sign of the topological charge and leaves the action invariant. Random applications of the parity transformation are proposed to traverse from one topological charge sign to the other. The transformation provides an improved unbiased estimator of the ensemble average and is essential in improving the ergodicity of the Markov chain process.

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.

preprint2003arXiv

Nonperturbative QCD, gauge-fixing, Gribov copies, and the lattice

Perturbative QCD uses the Faddeev-Popov gauge-fixing procedure, which leads to ghosts and the local BRST invariance of the gauge-fixed perturbative QCD action. In the asymptotic regime, where perturbative QCD is relevant, Gribov copies can be neglected. In the nonperturbative regime, one must adopt either a nonlocal Gribov-copy free gauge (e.g., Laplacian gauge) or attempt to maintain local BRST invariance at the expense of admitting Gribov copies. These issues are explored. In addition, we discuss the relationship between recent Dyson-Schwinger based model calculations of the infrared behavior of QCD Green's functions and the lattice calculation of these quantities.

preprint2003arXiv

The Hamiltonian limit of (3+1)D SU(3) lattice gauge theory on anisotropic lattices

The extreme anisotropic limit of Euclidean SU(3) lattice gauge theory is examined to extract the Hamiltonian limit, using standard path integral Monte Carlo (PIMC) methods. We examine the mean plaquette and string tension and compare them to results obtained within the Hamiltonian framework of Kogut and Susskind. The results are a significant improvement upon previous Hamiltonian estimates, despite the extrapolation procedure necessary to extract observables. We conclude that the PIMC method is a reliable method of obtaining results for the Hamiltonian version of the theory. Our results also clearly demonstrate the universality between the Hamiltonian and Euclidean formulations of lattice gauge theory. It is particularly important to take into account the renormalization of both the anisotropy, and the Euclidean coupling $ β_E $, in obtaining these results.

preprint2002arXiv

Improved Smoothing Algorithms for Lattice Gauge Theory

The relative smoothing rates of various gauge field smoothing algorithms are investigated on ${\cal O}(a^2)$-improved $\suthree$ Yang--Mills gauge field configurations. In particular, an ${\cal O}(a^2)$-improved version of APE smearing is motivated by considerations of smeared link projection and cooling. The extent to which the established benefits of improved cooling carry over to improved smearing is critically examined. We consider representative gauge field configurations generated with an ${\cal O}(a^2)$-improved gauge field action on $\1$ lattices at $β=4.38$ and $\2$ lattices at $β=5.00$ having lattice spacings of 0.165(2) fm and 0.077(1) fm respectively. While the merits of improved algorithms are clearly displayed for the coarse lattice spacing, the fine lattice results put the various algorithms on a more equal footing and allow a quantitative calibration of the smoothing rates for the various algorithms. We find the relative rate of variation in the action may be succinctly described in terms of simple calibration formulae which accurately describe the relative smoothness of the gauge field configurations at a microscopic level.

preprint2002arXiv

Overlap Quark Propagator in Landau Gauge

The properties of the quark propagator in Landau gauge in quenched QCD are examined for the overlap quark action. The overlap quark action satisfies the Ginsparg-Wilson relation and as such provides an exact lattice realization of chiral symmetry. This in turn implies that the quark action is free of ${\cal O}(a)$ errors. We present results using the standard Wilson fermion kernel in the overlap formalism on a $12^3\times 24$ lattice at a spacing of 0.125 fm. We obtain the nonperturbative momentum-dependent wavefunction renormalization function $Z(p)$ and the nonperturbative mass function $M(p)$ for a variety of bare masses. We perform a simple extrapolation to the chiral limit for these functions. We clearly observe the dynamically generated infrared mass and confirm the qualitative behavior found for the Landau gauge quark propagator in earlier studies. We attempt to extract the quark condensate from the asymptotic behavior of the mass function in the chiral limit.

preprint2001arXiv

General Algorithm For Improved Lattice Actions on Parallel Computing Architectures

Quantum field theories underlie all of our understanding of the fundamental forces of nature. The are relatively few first principles approaches to the study of quantum field theories [such as quantum chromodynamics (QCD) relevant to the strong interaction] away from the perturbative (i.e., weak-coupling) regime. Currently the most common method is the use of Monte Carlo methods on a hypercubic space-time lattice. These methods consume enormous computing power for large lattices and it is essential that increasingly efficient algorithms be developed to perform standard tasks in these lattice calculations. Here we present a general algorithm for QCD that allows one to put any planar improved gluonic lattice action onto a parallel computing architecture. High performance masks for specific actions (including non-planar actions) are also presented. These algorithms have been successfully employed by us in a variety of lattice QCD calculations using improved lattice actions on a 128 node Thinking Machines CM-5. {\underline{Keywords}}: quantum field theory; quantum chromodynamics; improved actions; parallel computing algorithms.

preprint2001arXiv

Nonperturbative improvement and tree-level correction of the quark propagator

We extend an earlier study of the Landau gauge quark propagator in quenched QCD where we used two forms of the O(a)-improved propagator with the Sheikholeslami-Wohlert quark action. In the present study we use the nonperturbative value for the clover coefficient c_sw and mean-field improvement coefficients in our improved quark propagators. We compare this to our earlier results which used the mean-field c_sw and tree-level improvement coefficients for the propagator. We also compare three different implementations of tree-level correction: additive, multiplicative, and hybrid. We show that the hybrid approach is the most robust and reliable and can successfully deal even with strong ultraviolet behavior and zero-crossing of the lattice tree-level expression. We find good agreement between our improved quark propagators when using the appropriate nonperturbative improvement coefficients and hybrid tree-level correction. We also present a simple extrapolation of the quark mass function to the chiral limit.

preprint2001arXiv

Quark propagator in the Landau gauge

The Landau gauge quark propagator in momentum space is investigated using the O(a)-improved Sheikholeslami--Wohlert (SW) quark action with a tree-level mean-field improved coefficient c_sw. We have studied the unimproved definition of the quark propagator, as well as two different tree-level O(a)-improved propagators. The ultraviolet behavior of the free lattice propagator is studied for each of these in order to establish which of them provides the most reliable description of the quark propagator up to the medium momentum regime. A general method of tree-level correction is introduced. This exploits asymptotic freedom and removes much of the trivial lattice artifacts at medium to high momenta. We obtain results for the quark propagator which are qualitatively similar to those typically used in quark models. A simple extrapolation of the infrared quark mass M(p^2=0) to the chiral limit gives (298+/-8+/-30) MeV, which is consistent with phenomenological expectations.

preprint2000arXiv

Calibration of Smearing and Cooling Algorithms in SU(3)-Color Gauge Theory

The action and topological charge are used to determine the relative rates of standard cooling and smearing algorithms in pure SU(3)-color gauge theory. We consider representative gauge field configurations on $16^3\times 32$ lattices at $β=5.70$ and $24^3\times 36$ lattices at $β=6.00$. We find the relative rate of variation in the action and topological charge under various algorithms may be succinctly described in terms of simple formulae. The results are in accord with recent suggestions from fat-link perturbation theory.

preprint1997arXiv

Solving the Bethe-Salpeter equation for bound states of scalar theories in Minkowski space

We apply the perturbation theory integral representation (PTIR) to solve for the bound state Bethe-Salpeter (BS) vertex for an arbitrary scattering kernel, without the need for any Wick rotation. The results derived are applicable to any scalar field theory (without derivative coupling). It is shown that solving directly for the BS vertex, rather than the BS amplitude, has several major advantages, notably its relative simplicity and superior numerical accuracy. In order to illustrate the generality of the approach we obtain numerical solutions using this formalism for a number of scattering kernels, including cases where the Wick rotation is not possible.

preprint1996arXiv

On Renormalized Strong-Coupling Quenched QED in Four Dimensions

We study renormalized quenched strong-coupling QED in four dimensions in arbitrary covariant gauge. Above the critical coupling leading to dynamical chiral symmetry breaking, we show that there is no finite chiral limit. This behaviour is found to be independent of the detailed choice of photon-fermion proper vertex in the Dyson-Schwinger equation formalism, provided that the vertex is consistent with the Ward-Takahashi identity and multiplicative renormalizability. We show that the finite solutions previously reported lie in an unphysical regime of the theory with multiple solutions and ultraviolet oscillations in the mass functions. This study supports the assertion that in four dimensions strong coupling QED does not have a continuum limit in the conventional sense.

preprint1996arXiv

Renormalization and Chiral Symmetry Breaking in Quenched QED in Arbitrary Covariant Gauge

We extend a previous Landau-gauge study of subtractive renormalization of the fermion propagator Dyson-Schwinger equation (DSE) in strong-coupling, quenched QED_4 to arbitrary covariant gauges. We use the fermion-photon proper vertex proposed by Curtis and Pennington with an additional correction term included to compensate for the small gauge-dependence induced by the ultraviolet regulator. We discuss the chiral limit and the onset of dynamical chiral symmetry breaking in the presence of nonperturbative renormalization. We extract the critical coupling in several different gauges and find evidence of a small residual gauge-dependence in this quantity.

preprint1995arXiv

SOLVING THE BETHE--SALPETER EQUATION IN MINKOWSKI SPACE: SCALAR THEORIES

The Bethe-Salpeter (BS) equation for scalar-scalar bound states in scalar theories without derivative coupling is formulated and solved in Minkowski space. This is achieved using the perturbation theory integral representation (PTIR), which allows these amplitudes to be expressed as integrals over weight functions and known singularity structures and hence allows us to convert the BS equation into an integral equation involving weight functions. We obtain numerical solutions using this formalism for a number of scattering kernels to illustrate the generality of the approach. It applies even when the na\"ıve Wick rotation is invalid. As a check we verify, for example, that this method applied to the special case of the massive ladder exchange kernel reproduces the same results as are obtained by Wick rotation.

preprint1995arXiv

Solving the Bethe-Salpeter Equation for Scalar Theories in Minkowski Space

The Bethe-Salpeter (BS) equation for scalar-scalar bound states in scalar theories without derivative coupling is formulated and solved in Minkowski space. This is achieved using the perturbation theory integral representation (PTIR), which allows these amplitudes to be expressed as integrals over weight functions and known singularity structures and hence allows us to convert the BS equation into an integral equation involving weight functions. We obtain numerical solutions using this formalism for a number of scattering kernels to illustrate the generality of the approach. It applies even when the na\"ıve Wick rotation is invalid. As a check we verify, for example, that this method applied to the special case of the massive ladder exchange kernel reproduces the same results as are obtained by Wick rotation.

preprint1994arXiv

Chiral Symmetry Breaking in Quenched Massive Strong-Coupling QED$_4$

We present results from a study of subtractive renormalization of the fermion propagator Dyson-Schwinger equation (DSE) in massive strong-coupling quenched QED$_4$. Results are compared for three different fermion-photon proper vertex {\it Ansätze\/}: bare $γ^μ$, minimal Ball-Chiu, and Curtis-Pennington. The procedure is straightforward to implement and numerically stable. This is the first study in which this technique is used and it should prove useful in future DSE studies, whenever renormalization is required in numerical work.

preprint1994arXiv

Dynamical chiral symmetry breaking and confinement with an infrared-vanishing gluon propagator?

We study a model Dyson-Schwinger equation for the quark propagator closed using an {\it Ansatz} for the gluon propagator of the form \mbox{$D(q) \sim q^2/[(q^2)^2 + b^4]$} and two {\it Ansätze} for the quark-gluon vertex: the minimal Ball-Chiu and the modified form suggested by Curtis and Pennington. Using the quark condensate as an order parameter, we find that there is a critical value of $b=b_c$ such that the model does not support dynamical chiral symmetry breaking for $b>b_c$. We discuss and apply a confinement test which suggests that, for all values of $b$, the quark propagator in the model {\bf is not} confining. Together these results suggest that this Ansatz for the gluon propagator is inadequate as a model since it does not yield the expected behaviour of QCD.