Source author record

Sara Collins

Sara Collins 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

33works
4topics
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

33 published item(s)

preprint2023arXiv

Doubly charmed tetraquark: isospin channels and diquark-antidiquark interpolators

We perform a lattice simulation to investigate the doubly charmed tetraquark $T^+_{cc}$ observed by the LHCb collaboration, slightly below the $D^{*+}D^0$ threshold, with flavor content $cc\bar{u}\bar{d}$ and isospin-$0$. Two-meson interpolators are implemented to explore the isospin quantum numbers $I=0$ and $I=1$. We observe attraction near the $DD^*$ threshold for $I=0$ and repulsion for $I=1$. Moreover, we also include diquark-antidiquark interpolators to study their effect on the energy spectrum. There is no significant shift in the ground state energy when adding diquark-antidiquark interpolators to the interpolator basis when the heavy quark mass is close to the physical charm quark mass. However, we observe a non-negligible shift in the second energy level. This effect has to be taken into account to extract the scattering amplitude of the $T^+_{cc}$. Finally, with a higher mass (close to the bottom quark), the ground state is shifted down significantly. The simulation is performed on $N_f=2+1$ CLS ensembles with $m_π\simeq 280$ MeV.

preprint2023arXiv

Sigma terms of the baryon octet in $N_\mathrm{f} = 2+1$ QCD with Wilson quarks

A lot of progress has been made in the direct determination of nucleon sigma terms. Using similar methods, we consider the sigma terms of the other octet baryons as well. These are determined on CLS gauge field ensembles employing the Lüscher-Weisz gluon action and the Sheikholeslami-Wohlert fermion action with $N_\mathrm{f} = 2 + 1$. The ensembles have pion masses ranging from ${410}\,\mathrm{MeV}$ down to the physical value and lattice spacings covering a range between ${0.098}\,\mathrm{fm}$ and ${0.039}\,\mathrm{fm}$. We present some preliminary results for the pion and strange sigma terms and compare to indirect determinations. To do so, we discuss multi-state fits to tackle the well-known problem of excited state contamination comparing the ratio and summation methods also including priors.

preprint2022arXiv

Determination of $m_c$ from $N_f = 2+1$ QCD with Wilson fermions

We present preliminary results for the charm quark mass in the $N_f=4$ RGI scheme. These were obtained using $N_f=2+1$ CLS ensembles with $\mathcal{O}(a)$ non-perturbatively improved Wilson fermions. We employed five different lattice spacings, ranging down to $a\lesssim 0.04$ fm and realized approximately physical pion and kaon masses, with ensembles spread out along three different trajectories in the quark mass plane, enabling a thorough study of the dependence on the lattice spacing and the light and strange sea quark masses. We sketch our analysis strategy and find that the dominant errors at present are due to the renormalization and scale setting uncertainties.

preprint2022arXiv

Leading order mesonic and baryonic SU(3) low energy constants from $N_f = 3$ lattice QCD

We determine the leading order mesonic~($B_0$ and $F_0$) and baryonic~($m_0$, $D$ and $F$) SU(3) chiral perturbation theory low energy constants from lattice QCD. We employ gauge ensembles with $N_f=3$ (i.e., $m_u=m_d=m_s$) non-perturbatively improved Wilson fermions at six distinct values of the lattice spacing in the range $a\approx (0.039 - 0.098)$ fm, which constitute a subset of the Coordinated Lattice Simulations (CLS) gauge ensembles. The pseudoscalar meson mass $M_π$ ranges from around $430$ MeV down to $240$ MeV and the linear spatial lattice extent $L$ from $6.4\,M_π^{-1}$ to $3.3\,M_π^{-1}$, where $ L M_π\geq 4$ for the majority of the ensembles. This allows us to perform a controlled extrapolation of all the low energy constants to the chiral, infinite volume and continuum limits. We find the SU(3) chiral condensate and $F_0$ to be smaller than their SU(2) counterparts while the Gell-Mann--Oakes--Renner parameters $B_0\approx B$ are similar. Regarding baryonic LECs, we obtain $F/D = 0.612^{(14)}_{(12)}$.

preprint2020arXiv

Nucleon axial structure from lattice QCD

We present a new analysis method that allows one to understand and model excited state contributions in observables that are dominated by a pion pole. We apply this method to extract axial and (induced) pseudoscalar nucleon isovector form factors, which satisfy the constraints due to the partial conservation of the axial current up to expected discretization effects. Effective field theory predicts that the leading contribution to the (induced) pseudoscalar form factor originates from an exchange of a virtual pion, and thus exhibits pion pole dominance. Using our new method, we can recover this behavior directly from lattice data. The numerical analysis is based on a large set of ensembles generated by the CLS effort, including physical pion masses, large volumes (with up to $96^3 \times 192$ sites and $L m_π= 6.4$), and lattice spacings down to $0.039 \, \text{fm}$, which allows us to take all the relevant limits. We find that some observables are much more sensitive to the choice of parametrization of the form factors than others. On the one hand, the $z$-expansion leads to significantly smaller values for the axial dipole mass than the dipole ansatz ($M_A^{\text{$z$-exp}}=1.02(10) \, \text{GeV}$ versus $M_A^{\text{dipole}} = 1.31(8) \, \text{GeV}$). On the other hand, we find that the result for the induced pseudoscalar coupling at the muon capture point is almost independent of the choice of parametrization ($g_P^{\star \ \text{$z$-exp}} = 8.68(45)$ and $g_P^{\star \ \text{dipole}} = 8.30(24)$), and is in good agreement with both, chiral perturbation theory predictions and experimental measurement via ordinary muon capture. We also determine the axial coupling constant $g_A$.

preprint2016arXiv

$ρ$ and $K^*$ resonances on the lattice at nearly physical quark masses and $N_f=2$

Working with a pion mass $m_π\approx 150$ MeV, we study $ππ$ and $Kπ$ scattering using two flavours of non-perturbatively improved Wilson fermions at a lattice spacing $a\approx 0.071$ fm. Employing two lattice volumes with linear spatial extents of $N_s=48$ and $N_s=64$ points and moving frames, we extract the phase shifts for p-wave $ππ$ and $Kπ$ scattering near the $ρ$ and $K^*$ resonances.Comparing our results to those of previous lattice studies, that used pion masses ranging from about 200 MeV up to 470 MeV, we find that the coupling $g_{ρππ}$ appears to be remarkably constant as a function of $m_π$.

preprint2016arXiv

Direct determinations of the nucleon and pion $σ$ terms at nearly physical quark masses

We present a high statistics study of the pion and nucleon light and strange quark sigma terms using $N_f=2$ dynamical non-perturbatively improved clover fermions with a range of pion masses down to $m_π\sim 150$ MeV and several volumes, $Lm_π=3.4$ up to $6.7$, and lattice spacings, $a=0.06-0.08$ fm, enabling a study of finite volume and discretisation effects for $m_π\gtrsim 260$ MeV. Systematics are found to be reasonably under control. For the nucleon we obtain $σ_{πN}=35(6)$ MeV and $σ_s=35(12)$ MeV, or equivalently in terms of the quark fractions, $f_{T_u}=0.021(4)$, $f_{T_d}=0.016(4)$ and $f_{T_s}=0.037(13)$, where the errors include estimates of both the systematic and statistical uncertainties. These values, together with perturbative matching in the heavy quark limit, lead to $f_{T_c}=0.075(4)$, $f_{T_b}=0.072(2)$ and $f_{T_t}=0.070(1)$. In addition, through the use of the (inverse) Feynman-Hellmann theorem our results for $σ_{πN}$ are shown to be consistent with the nucleon masses determined in the analysis. For the pion we implement a method which greatly reduces excited state contamination to the scalar matrix elements from states travelling across the temporal boundary. This enables us to demonstrate the Gell-Mann-Oakes-Renner expectation $σ_π=m_π/2$ over our range of pion masses.

preprint2016arXiv

Nucleon generalized form factors from lattice QCD with nearly physical quark masses

We determine generalized form factors of the nucleon from lattice simulations with $N_f = 2$ mass-degenerate non-perturbatively improved Wilson-Sheikholeslami-Wohlert fermions down to a pion mass of 150 MeV. We also present the resulting isovector quark angular momentum. Possible excited-state contaminations are investigated with correlated simultaneous fits.

preprint2016arXiv

Testing the hadro-quarkonium model on the lattice

Recently the LHCb experiment found evidence for the existence of two exotic resonances consisting of $c\bar{c}uud$ quarks. Among the possible interpretations is the hadro-charmonium model, in which charmonium is bound "within" a light hadron. We test this idea on CLS $N_f$=2+1 lattices using the static formulation for the heavy quarks. We find that the static potential is modified by the presence of a hadron such that it becomes more attractive. The effect is of the order of a few MeV.

preprint2015arXiv

(Approximate) Low-Mode Averaging with a new Multigrid Eigensolver

We present a multigrid based eigensolver for computing low-modes of the Hermitian Wilson Dirac operator. For the non-Hermitian case multigrid methods have already replaced conventional Krylov subspace solvers in many lattice QCD computations. Since the $γ_5$-preserving aggregation based interpolation used in our multigrid method is valid for both, the Hermitian and the non-Hermitian case, inversions of very ill-conditioned shifted systems with the Hermitian operator become feasible. This enables the use of multigrid within shift-and-invert type eigensolvers. We show numerical results from our MPI-C implementation of a Rayleigh quotient iteration with multigrid. For state-of-the-art lattice sizes and moderate numbers of desired low-modes we achieve speed-ups of an order of magnitude and more over PARPACK. We show results and develop strategies how to make use of our eigensolver for calculating disconnected contributions to hadronic quantities that are noisy and still computationally challenging. Here, we explore the possible benefits, using our eigensolver for low-mode averaging and related methods with high and low accuracy eigenvectors. We develop a low-mode averaging type method using only a few of the smallest eigenvectors with low accuracy. This allows us to avoid expensive exact eigensolves, still benefitting from reduced statistical errors.

preprint2015arXiv

$D_s \rightarrow η, η'$ semileptonic decay form factors with disconnected quark loop contributions

We calculate for the first time the form factors of the semi-leptonic decays of the $D_s$ meson to $η$ and $η^\prime$ using lattice techniques. As a by-product of the calculation we obtain the masses and leading distribution amplitudes of the $η$ and $η^\prime$ mesons. We use $N_f=2+1$ non-perturbatively improved clover fermions on configurations with a lattice spacing $a\sim 0.075$ fm. We are able to obtain clear signals for relevant matrix elements, using several noise reduction techniques, both for the connected and disconnected contributions. This includes a new method for reducing the variance of pseudoscalar disconnected two-point functions. At zero momentum transfer, we obtain for the scalar form factors, $|f_0^{D_s\to η}|=0.564(11)$ and $|f_0^{D_s\to η'}|=0.437(18)$ at $M_π\approx 470\, {\rm MeV}$, as well as $|f_0^{D_s\to η}|=0.542(13)$ and $|f_0^{D_s\to η'}|=0.404(25)$ at $M_π\approx 370\, {\rm MeV}$, where the errors are statistical only.

preprint2015arXiv

Charmed baryon spectroscopy and light flavour symmetry from lattice QCD

We determine the ground state and first excited state masses of singly and doubly charmed spin 1/2 and 3/2 baryons with positive and negative parity. Configurations with $N_f=2+1$ non-perturbatively improved Wilson-clover fermions were employed, with the same quark action also being used for the valence quarks, including the charm. The spectrum is calculated for pion masses in the range $M_π\sim 259-460$ MeV at a lattice spacing $a\sim 0.075$ fm. Finite volume effects are studied comparing lattices with two different linear spatial extents ($1.8\,{\rm fm}$ and $2.4\,{\rm fm}$). The physical point is approached from the SU(3) limit keeping the flavour averaged light quark mass fixed. The baryon masses are extrapolated using expansions in the strange-light quark mass difference. Most particles fall into the expected SU(3) multiplets with well constrained extrapolations, the exceptions having a possibly more complex internal structure. Overall agreement is found with experiment for the masses and splittings of the singly charmed baryons. As part of the calculation an analysis of the lower lying charmonium, $D$ and $D_s$ spectra was performed in order to assess discretisation errors. The gross spectra are reproduced, including the $D^*_{s0}$, $D_{s1}$ and $D_1$ mesons, while at this single lattice spacing hyperfine splittings come out $10-20$ MeV too low.

preprint2015arXiv

Nucleon isovector couplings from $N_f=2$ lattice QCD

We compute the axial, scalar, tensor and pseudoscalar isovector couplings of the nucleon as well as the induced tensor and pseudoscalar charges in lattice simulations with $N_f=2$ mass-degenerate non-perturbatively improved Wilson-Sheikholeslami-Wohlert fermions. The simulations are carried out down to a pion mass of 150 MeV and linear spatial lattice extents of up to 4.6 fm at three different lattice spacings ranging from approximately 0.08 fm to 0.06 fm. Possible excited state contamination is carefully investigated and finite volume effects are studied. The couplings, determined at these lattice spacings, are extrapolated to the physical pion mass. In this limit we find agreement with experimental results, where these exist, with the exception of the magnetic moment. A proper continuum limit could not be performed, due to our limited range of lattice constants, but no significant lattice spacing dependence is detected. Upper limits on discretization effects are estimated and these dominate the error budget.

preprint2014arXiv

The moment $\langle x\rangle_{u-d}$ of the nucleon from $N_f=2$ lattice QCD down to nearly physical quark masses

We present an update of our analysis [1] which includes additional ensembles at different quark masses, lattice spacings and volumes, all with high statistics. We use $N_f=2$ mass-degenerate quark flavours, employing the non-perturbatively improved clover action. The lattice matrix elements are converted to the $\overline{\rm MS}$ scheme via renormalization factors determined non-perturbatively in the RI$^\prime$-MOM scheme. We have systematically investigated excited state contributions, in particular, at the smallest, near physical, pion mass. While our results~(with much increased precision) are consistent with Ref.~[1], comparing with previous determinations we find that excited state contributions can be significant if the quark smearing is not suitably optimized, in agreement with other recent studies. The difference with respect to the value for $\langle x\rangle_{u-d}$ extracted from experimental data is reduced but not resolved. Using lattice sizes in the range $L m_π\sim 3.4-6.7$, no significant finite volume effects were observed. Performing a controlled continuum limit that may remove the discrepancy will require simulations at lattice spacings $a< 0.06$ fm.

preprint2013arXiv

Mesons in large-N QCD

We present the results of a systematic, first-principles study of the spectrum and decay constants of mesons for different numbers of color charges N, via lattice computations. We restrict our attention to states in the non-zero isospin sector, evaluating the masses associated with the ground-state and first excitation in the pseudoscalar, vector, scalar, and axial vector channels. Our results are based on a new set of simulations of four dimensional SU(N) Yang-Mills theories with the number of colors ranging from N=2 to N=17; the spectra and the decay constants are computed in the quenched approximation (which becomes exact in the 't Hooft limit) using Wilson fermions. After discussing the extrapolations to the chiral and large-N limits, we present a comparison of our results to some of the numerical computations and analytical predictions available in the literature - including, in particular, those from holographic computations.

preprint2013arXiv

Nucleon structure from stochastic estimators

Using stochastic estimators for connected meson and baryon three-point functions has successfully been tried in the past years. Compared to the standard sequential source method we trade the freedom to compute the current-to-sink propagator independently of the hadron sink for additional stochastic noise in our observables. In the case of the nucleon we can use this freedom to compute many different sink-momentum/polarization combinations, which grants access to more virtualities. We will present preliminary results on the scalar, electro-magnetic and axial form factors of the nucleon in $N_f=2+1$ lattice QCD and contrast the performance of the stochastic method to the sequential source method. We find the stochastic method to be competitive in terms of errors at fixed cost.

preprint2012arXiv

<x>_{u-d} from lattice QCD at nearly physical quark masses

We determine the second Mellin moment of the isovector quark parton distribution function <x>_{u-d} from lattice QCD with N_f=2 sea quark flavours, employing the non-perturbatively improved Wilson-Sheikholeslami-Wohlert action at a pseudoscalar mass of 157(6) MeV. The result is converted non-perturbatively to the RI'-MOM scheme and then perturbatively to the MSbar scheme at a scale mu = 2 GeV. As the quark mass is reduced we find the lattice prediction to approach the value extracted from experiments.

preprint2012arXiv

Strangeness Contribution to the Proton Spin from Lattice QCD

We compute the strangeness and light-quark contributions Delta s, Delta u and Delta d to the proton spin in nf=2 lattice QCD at a pion mass of about 285 MeV and at a lattice spacing a approx 0.073 fm, using the non-perturbatively improved Sheikholeslami-Wohlert Wilson action. We carry out the renormalization of these matrix elements which involves mixing between contributions from different quark flavours. Our main result is the small negative value MSbar scheme value Delta s = -0.020(10)(4), at a scale mu = 7.4 GeV, of the strangeness contribution to the nucleon spin. The second error is an estimate of the uncertainty, due to the missing extrapolation to the physical point.

preprint2012arXiv

The strange and light quark contributions to the nucleon mass from Lattice QCD

We determine the strangeness and light quark fractions of the nucleon mass by computing the quark line connected and disconnected contributions to the matrix elements m_q <N|qbar q|N> in lattice QCD, using the non-perturbatively improved Sheikholeslami-Wohlert Wilson Fermionic action. We simulate n_F=2 mass degenerate sea quarks with a pion mass of about 285 MeV and a lattice spacing a approx 0.073 fm. The renormalization of the matrix elements involves mixing between contributions from different quark flavours. The pion-nucleon sigma-term is extrapolated to physical quark masses exploiting the sea quark mass dependence of the nucleon mass. We obtain the renormalized values σ_{πN} = 38(12) MeV at the physical point and f_{T_s}=σ_s/m_N= 0.012(14)^{+10}_{-3} for the strangeness contribution at our larger than physical sea quark mass.

preprint2011arXiv

Charmonium spectroscopy and mixing with light quark and open charm states from nF=2 lattice QCD

We study the charmonium spectrum including higher spin and gluonic excitations. We determine an upper limit on the mixing of the eta_c ground state with light pseudoscalar flavour-singlet mesons and investigate the mixing of charmonia near open charm thresholds with pairs of (excited) D and anti-D mesons. For charm and light valence quarks and nF=2 sea quarks, we employ the non-perturbatively improved Sheikholeslami-Wohlert (clover) action. Excited states are accessed using the variational technique, starting from a basis of suitably optimised operators. For some aspects of this study, the use of improved stochastic all-to-all propagators was essential.

preprint2010arXiv

Effective noise reduction techniques for disconnected loops in Lattice QCD

Many Lattice QCD observables of phenomenological interest include so-called all-to-all propagators. The computation of these requires prohibitively large computational resources, unless they are estimated stochastically. This is usually done. However, the computational demand can often be further reduced by one order of magnitude by implementing sophisticated unbiased noise reduction techniques. We combine both well known and novel methods that can be applied to a wide range of problems. We concentrate on calculating disconnected contributions to nucleon structure functions, as one realistic benchmark example. In particular we determine the strangeness contributions to the nucleon, <N|ss|N>, and to the spin of the nucleon, Delta s.

preprint2010arXiv

Improved Semileptonic Form Factor Calculations in Lattice QCD

We investigate the computational efficiency of two stochastic based alternatives to the Sequential Propagator Method used in Lattice QCD calculations of heavy-light semileptonic form factors. In the first method, we replace the sequential propagator, which couples the calculation of two of the three propagators required for the calculation, with a stochastic propagator so that the calculations of all three propagators are independent. This method is more flexible than the Sequential Propagator Method but introduces stochastic noise. We study the noise to determine when this method becomes competitive with the Sequential Propagator Method, and find that for any practical calculation it is competitive with or superior to the Sequential Propagator Method. We also examine a second stochastic method, the so-called ``one-end trick", concluding it is relatively inefficient in this context. The investigation is carried out on two gauge field ensembles, using the non-perturbatively improved Wilson-Sheikholeslami-Wohlert action with N_f=2 mass-degenerate sea quarks. The two ensembles have similar lattice spacings but different sea quark masses. We use the first stochastic method to extract ${\mathcal O}(a)$-improved, matched lattice results for the semileptonic form factors on the ensemble with lighter sea quarks, extracting f_+(0).

preprint2010arXiv

Meson and baryon masses with low mode averaging

We describe and test a method known in the literature as low mode averaging to improve Euclidean two-point functions in lattice QCD using the low-lying eigenmodes of the Wilson-Dirac operator D. The contribution from the low modes is averaged over all positions of the quark sources while the contribution from high modes is calculated in the traditional way using one source point per lattice. We apply this method to different baryon and meson two-point functions and we compare the improvements using the eigenmodes of the non-hermitian operator D and the eigenmodes of the hermitian operator Q=gamma5 D. The convergence strongly depends on the parity of the states.

preprint2009arXiv

A Stochastic Method for Semileptonic Form Factor Calculations on the Lattice

We investigate an alternative to the Sequential Propagator Method used in Lattice QCD calculations of semileptonic form factors. We replace the sequential propagator with a stochastic propagator so that, in principle, all momentum and sink smearing combinations are available with only a single spin-color inversion. Practically, the stochastic noise is significant and must be reduced at the cost of more inversions. We study the behavior of the stochastic noise and compare the computational costs of this stochastic technique and the Sequential Propagator Method. We also present preliminary semileptonic form factor results using the stochastic technique on N_f=2 configurations with a non-perturbatively improved Sheikoleslami-Wohlert action generated by the QCDSF collaboration. At a fixed cost, measured in terms of the number of heavy-quark inversions, the method provides more correlators for the extraction of the form factors at various q^2's than the Sequential Propagator Method. These additional correlators reduce the total statistical errors of certain kinematic points, although the stochastic error is still comparable to the gauge error at other points.

preprint2009arXiv

Strangeness and charm content of the nucleon

We present results on the scalar strangeness and charm contents of the nucleon and of the disconnected contributions to the nucleon spin. These are obtained on two flavour non-perturbatively improved Sheikholeslami-Wilson configurations at a pseudoscalar mass of 290 MeV. We quote f_{T_s} = m_s<N|ss|N>/m_N< = 0.070(22) as our preliminary value for the strange quark fraction of the nucleon mass and a modified minimal subtraction scheme value Delta s = -0.015(10), with as yet unknown systematics, for the strangeness contribution to the spin.

preprint1997arXiv

SCRI Results With the Tadpole-Improved Clover Action

We compare light hadron spectroscopy using the Wilson and Clover fermionic actions. We show that a Clover coefficient chosen using tadpole-improved tree-level perturbation theory effectively eliminates the O(a) discretization errors present in the Wilson action. We find that discretization errors in light spectroscopy for both the Wilson and Clover actions are characterized by an energy scale mu of about 200-300 MeV, indicating that these errors can be reduced to the 5% level by using the Clover action at an inverse lattice spacing of about 1.3 GeV.

preprint1996arXiv

Comparing Wilson and Clover Quenched $SU(3)$ Spectroscopy with an Improved Gauge Action

We present results of quenched $SU(3)$ hadron spectroscopy comparing $\order(a)$ improved Wilson (Clover) fermions with conventional Wilson fermions. The configurations were generated using an $\order(a^2)$ improved 6-link $SU(3)$ pure gauge action at $β$'s corresponding to lattice spacings of $0.15$, $0.18$, $0.20$, $0.33$, and $0.43$ fm. We find evidence that fermionic scaling violations are consistent with $\order(a^2)$ for Clover and $\order(a)$ with a nonnegligible $\order(a^2)$ term for standard Wilson fermions. This latter mixed ansatz makes a reliable continuum extrapolation problematic for Wilson fermions. We also find that the slope of the scaling violations is roughly $250 MeV$ for both Wilson and Clover fermions.

preprint1995arXiv

A Comparison of Clover and Wilson Spectroscopy in the Presence of Dynamical Quarks

We present preliminary results of light hadron spectroscopy using valence, tadpole-improved, Clover fermions on an ensemble of gauge configurations generated with 2 flavors of staggered fermions at a beta of 5.6. We compare the slope and intercept of the curve M_V vs. M_PS^2 for Clover and Wilson fermions. We show that a higher order chiral perturbation theory ansatz works very well for chiral extrapolations.

preprint1995arXiv

Quenched $SU(3)$ hadron spectroscopy using improved fermionic and gauge actions

We present results of quenched $SU(3)$ hadron spectroscopy using $\order(a)$ improved Wilson fermions. The configurations were generated using an $\order(a^2)$ improved 6-link $SU(3)$ pure gauge action at $β$'s corresponding to lattice spacings of $0.43$, $0.25$, $0.20$, $0.18$, and $0.15$ fm. We find evidence that fermionic scaling violations are consistent with $\order(a^2)$ errors.