Source author record

A. Schiller

A. Schiller 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

76works
9topics
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

76 published item(s)

preprint2019arXiv

Anomalous magnetic moment of the muon with dynamical QCD+QED

The current $3.5σ$ discrepancy between experimental and Standard Model determinations of the anomalous magnetic moment of the muon $a_μ=(g-2)/2$ can only be extended to the discovery $5σ$ regime through a reduction of both experimental and theoretical uncertainties. On the theory side, this means a determination of the hadronic vacuum polarisation (HVP) contribution to better than 0.5%, a level of precision that demands the inclusion of QCD + QED effects to properly understand how the behaviour of quarks are modified when their electric charges are turned on. The QCDSF collaboration has generated an ensemble of configurations with dynamical QCD and QED fields with the specific aim of studying flavour breaking effects arising from differences in the quark masses and charges in physical quantities. Here we study these effects in a calculation of HVP around the SU(3) symmetric point. Furthermore, by performing partially-quenched simulations we are able to cover a larger range of quark masses and charges on these configurations and then fit the results to an SU(3) flavour breaking expansion. Subsequently, this allows for an extrapolation to the physical point.

preprint2016arXiv

Isospin splittings of meson and baryon masses from three-flavor lattice QCD + QED

Lattice QCD simulations are now reaching a precision where isospin breaking effects become important. Previously, we have developed a program to systematically investigate the pattern of flavor symmetry beaking within QCD and successfully applied it to meson and baryon masses involving up, down and strange quarks. In this Letter we extend the calculations to QCD + QED and present our first results on isospin splittings in the pseudoscalar meson and baryon octets. In particular, we obtain the nucleon mass difference of $M_n-M_p=1.35(18)(8)\,\mbox{MeV}$ and the electromagnetic contribution to the pion splitting $M_{π^+}-M_{π^0}=4.60(20)\,\mbox{MeV}$. Further we report first determination of the separation between strong and electromagnetic contributions in the $\bar{MS}$ scheme.

preprint2016arXiv

Partially conserved axial vector current and applications

We investigate implications of the use of the point-split axial vector current derived from a Wilson like fermionic action. We compute the corresponding renormalization factor nonperturbatively for one beta value. The axial charge gA calculated from this nonlocal current is found to be nearer to the physical value than computed with the local axial vector current -- computed both on the same lattice with the same action.

preprint2016arXiv

QED effects in the pseudoscalar meson sector

We present results on the pseudoscalar meson masses from a fully dynamical simulation of QCD+QED. We concentrate particularly on violations of isospin symmetry. We calculate the $π^+$-$π^0$ splitting and also look at other isospin violating mass differences. We have presented results for these isospin splittings in arXiv:1508.06401 [hep-lat]. In this paper we give more details of the techniques employed, discussing in particular the question of how much of the symmetry violation is due to QCD, arising from the different masses of the $u$ and $d$ quarks, and how much is due to QED, arising from the different charges of the quarks. This decomposition is not unique, it depends on the renormalisation scheme and scale. We suggest a renormalisation scheme in which Dashen's theorem for neutral mesons holds, so that the electromagnetic self-energies of the neutral mesons are zero, and discuss how the self-energies change when we transform to a scheme such as $\bar{MS}$, in which Dashen's theorem for neutral mesons is violated.

preprint2015arXiv

A lattice determination of Sigma - Lambda mixing

Isospin breaking effects in baryon octet (and decuplet) masses are due to a combination of up and down quark mass differences and electromagnetic effects and lead to small mass splittings. Between the Sigma and Lambda this mass splitting is much larger, this being mostly due to their different wavefunctions. However when isospin is broken, there is a mixing between between these states. We describe the formalism necessary to determine the QCD mixing matrix and hence find the mixing angle and mass splitting between the Sigma and Lambda particles due to QCD effects.

preprint2015arXiv

Applications of the Feynman-Hellmann theorem in hadron structure

The Feynman-Hellmann (FH) relation offers an alternative way of accessing hadronic matrix elements through artificial modifications to the QCD Lagrangian. In particular, a FH-motivated method provides a new approach to calculations of disconnected contributions to matrix elements and high-momentum nucleon and pion form factors. Here we present results for the total nucleon axial charge, including a statistically significant non-negative total disconnected quark contribution of around $-5\%$ at an unphysically heavy pion mass. Extending the FH relation to finite-momentum transfers, we also present calculations of the pion and nucleon electromagnetic form factors up to momentum transfers of around 7-8 GeV$^2$. Results for the nucleon are not able to confirm the existence of a sign change for the ratio $\frac{G_E}{G_M}$, but suggest that future calculations at lighter pion masses will provide fascinating insight into this behaviour at large momentum transfers.

preprint2015arXiv

Disconnected contributions to the spin of the nucleon

The spin decomposition of the proton is a long-standing topic of much interest in hadronic physics. Lattice QCD has had much success in calculating the connected contributions to the quark spin. However, complete calculations, which necessarily involve gluonic and strange-quark contributions, still present some challenges. These "disconnected" contributions typically involve small signals hidden against large statistical backgrounds and rely on computationally intensive stochastic techniques. In this work we demonstrate how a Feynman-Hellmann approach may be used to calculate such quantities, by measuring shifts in the proton energy arising from artificial modifications to the QCD action. We find a statistically significant non-zero result for the disconnected quark spin contribution to the proton of about -5% at a pion mass of 470 MeV.

preprint2015arXiv

Improving the lattice axial vector current

For Wilson and clover fermions traditional formulations of the axial vector current do not respect the continuum Ward identity which relates the divergence of that current to the pseudoscalar density. Here we propose to use a point-split or one-link axial vector current whose divergence exactly satisfies a lattice Ward identity, involving the pseudoscalar density and a number of irrelevant operators. We check in one-loop lattice perturbation theory with SLiNC fermion and gauge plaquette action that this is indeed the case including order $O(a)$ effects. Including these operators the axial Ward identity remains renormalisation invariant. First preliminary results of a nonperturbative check of the Ward identity are also presented.

preprint2015arXiv

Reply to "Comment on `Lattice determination of Sigma - Lambda mixing' "

In this Reply, we respond to the above Comment. Our computation [Phys. Rev. D 91 (2015) 074512] only took into account pure QCD effects, arising from quark mass differences, so it is not surprising that there are discrepancies in isospin splittings and in the Sigma - Lambda mixing angle. We expect that these discrepancies will be smaller in a full calculation incorporating QED effects.

preprint2015arXiv

The electric dipole moment of the neutron from 2+1 flavor lattice QCD

We compute the electric dipole moment d_n of the neutron from a fully dynamical simulation of lattice QCD with 2+1 flavors of clover fermions and nonvanishing theta term. The latter is rotated into the pseudoscalar density in the fermionic action using the axial anomaly. To make the action real, the vacuum angle theta is taken to be purely imaginary. The physical value of d_n is obtained by analytic continuation. We find d_n = -3.8(2)(9) x 10^{-16} [theta e cm], which, when combined with the experimental limit on d_n, leads to the upper bound theta < 7.6 x 10^{-11}.

preprint2015arXiv

Wilson flow and scale setting from lattice QCD

We give a determination of the phenomenological value of the Wilson (or gradient) flow scales t0 and w0 for 2+1 flavours of dynamical quarks. The simulations are performed keeping the average quark mass constant, which allows the approach to the physical point to be made in a controlled manner. O(a) improved clover fermions are used and together with four lattice spacings this allows the continuum extrapolation to be taken.

preprint2014arXiv

A Feynman-Hellmann approach to the spin structure of hadrons

We perform a Nf = 2 + 1 lattice QCD simulation to determine the quark spin fractions of hadrons using the Feynman-Hellmann theorem. By introducing an external spin operator to the fermion action, the matrix elements relevant for quark spin fractions are extracted from the linear response of the hadron energies. Simulations indicate that the Feynman-Hellmann method offers statistical precision that is comparable to the standard three-point function approach, with the added benefit that it is less susceptible to excited state contamination. This suggests that the Feynman-Hellmann technique offers a promising alternative for calculations of quark line disconnected contributions to hadronic matrix elements. At the SU(3)-flavour symmetry point, we find that the connected quark spin fractions are universally in the range 55-70% for vector mesons and octet and decuplet baryons. There is an indication that the amount of spin suppression is quite sensitive to the strength of SU(3) breaking.

preprint2014arXiv

A novel approach to nonperturbative renormalization of singlet and nonsinglet lattice operators

A novel method for nonperturbative renormalization of lattice operators is introduced, which lends itself to the calculation of renormalization factors for nonsinglet as well as singlet operators. The method is based on the Feynman-Hellmann relation, and involves computing two-point correlators in the presence of generalized background fields arising from introducing additional operators into the action. As a first application, and test of the method, we compute the renormalization factors of the axial vector current $A_μ$ and the scalar density $S$ for both nonsinglet and singlet operators for $N_f=3$ flavors of SLiNC fermions. For nonsinglet operators, where a meaningful comparison is possible, perfect agreement with recent calculations using standard three-point function techniques is found.

preprint2014arXiv

Connected and disconnected quark contributions to hadron spin

By introducing an external spin operator to the fermion action, the quark spin fractions of hadrons are determined from the linear response of the hadron energies using the Feynman-Hellmann (FH) theorem. At our SU(3)-flavour symmetric point, we find that the connected quark spin fractions are universally in the range 55-70\% for vector mesons and octet and decuplet baryons. There is an indication that the amount of spin suppression is quite sensitive to the strength of SU(3) breaking. We also present first preliminary results applying the FH technique to calculations of quark-line disconnected contributions to hadronic matrix elements of axial and tensor operators. At the SU(3)-flavour symmetric point we find a small negative contribution to the nucleon spin from disconnected quark diagrams, while the corresponding tensor matrix elements are consistent with zero.

preprint2014arXiv

Determining Sigma - Lambda mixing

SU2 isospin breaking effects in baryon octet (and decuplet) masses are due to a combination of up and down quark mass differences and electromagnetic effects. These mass differences are small. Between the Sigma and Lambda the mass splitting is much larger, but this is mostly due to their different wavefunctions. However there is now also mixing between these states. We determine the QCD mixing matrix and hence find the mixing angle and mass splitting.

preprint2014arXiv

Renormalization of local quark-bilinear operators for Nf=3 flavors of SLiNC fermions

The renormalization factors of local quark-bilinear operators are computed non-perturbatively for $N_f=3$ flavors of SLiNC fermions, with emphasis on the various procedures for the chiral and continuum extrapolations. The simulations are performed at a lattice spacing $a=0.074$ fm, and for five values of the pion mass in the range of 290-465 MeV, allowing a safe and stable chiral extrapolation. Emphasis is given in the subtraction of the well-known pion pole which affects the renormalization factor of the pseudoscalar current. We also compute the inverse propagator and the Green's functions of the local bilinears to one loop in perturbation theory. We investigate lattice artifacts by computing them perturbatively to second order as well as to all orders in the lattice spacing. The renormalization conditions are defined in the RI$'$-MOM scheme, for both the perturbative and non-perturbative results. The renormalization factors, obtained at different values of the renormalization scale, are translated to the ${\bar{\rm MS}}$ scheme and are evolved perturbatively to 2 GeV. Any residual dependence on the initial renormalization scale is eliminated by an extrapolation to the continuum limit. We also study the various sources of systematic errors. Particular care is taken in correcting the non-perturbative estimates by subtracting lattice artifacts computed to one loop perturbation theory using the same action. We test two different methods, by subtracting either the ${\cal O}(g^2\,a^2)$ contributions, or the complete (all orders in $a$) one-loop lattice artifacts.

preprint2013arXiv

Perturbatively improving renormalization constants

Renormalization factors relate the observables obtained on the lattice to their measured counterparts in the continuum in a suitable renormalization scheme. They have to be computed very precisely which requires a careful treatment of lattice artifacts. In this work we present a method to suppress these artifacts by subtracting one-loop contributions proportional to the square of the lattice spacing calculated in lattice perturbation theory.

preprint2013arXiv

Perturbatively improving RI-MOM renormalization constants

The determination of renormalization factors is of crucial importance in lattice QCD. They relate the observables obtained on the lattice to their measured counterparts in the continuum in a suitable renormalization scheme. Therefore, they have to be computed as precisely as possible. A widely used approach is the nonperturbative Rome-Southampton method. It requires, however, a careful treatment of lattice artifacts. In this paper we investigate a method to suppress these artifacts by subtracting one-loop contributions to renormalization factors calculated in lattice perturbation theory. We compare results obtained from a complete one-loop subtraction with those calculated for a subtraction of contributions proportional to the square of the lattice spacing.

preprint2013arXiv

SU(3) flavour symmetry breaking and charmed states

By extending the SU(3) flavour symmetry breaking expansion from up, down and strange sea quark masses to partially quenched valence quark masses we propose a method to determine charmed quark hadron masses including possible QCD isospin breaking effects. Initial results for some open charmed pseudoscalar meson states and singly and doubly charmed baryon states are encouraging and demonstrate the potential of the procedure. Essential for the method is the determination of the scale using singlet quantities, and to this end we also give here a preliminary estimation of the recently introduced Wilson flow scales.

preprint2013arXiv

Transitional $γ$ strength in Cd isotopes

The level densities and $γ$-ray strength functions of $^{105,106,111,112}$Cd have been extracted from particle-$γ$ coincidence data using the Oslo method. The level densities are in very good agreement with known levels at low excitation energy. The $γ$-ray strength functions display no strong enhancement for low $γ$ energies. However, more low-energy strength is apparent for $^{105,106}$Cd than for $^{111,112}$Cd. For $γ$ energies above $\approx$ 4 MeV, there is evidence for some extra strength, similar to what has been previously observed for the Sn isotopes. The origin of this extra strength is unclear; it might be due to $E1$ and $M1$ transitions originating from neutron skin oscillations or the spin-flip resonance, respectively.

preprint2012arXiv

Analysis of possible systematic errors in the Oslo method

In this work, we have reviewed the Oslo method, which enables the simultaneous extraction of level density and gamma-ray transmission coefficient from a set of particle-gamma coincidence data. Possible errors and uncertainties have been investigated. Typical data sets from various mass regions as well as simulated data have been tested against the assumptions behind the data analysis.

preprint2012arXiv

Evidence for the ground-state resonance of 26O

Evidence for the ground state of the neutron-unbound nucleus 26O was observed for the first time in the single proton-knockout reaction from a 82 MeV/u 27F beam. Neutrons were measured in coincidence with 24O fragments. 26O was determined to be unbound by 150+50-150 keV from the observation of low-energy neutrons. This result agrees with recent shell model calculations based on microscopic two- and three-nucleon forces.

preprint2012arXiv

Exploring the neutron dripline two neutrons at a time: The first observations of the 26O and 16Be ground state resonances

The two-neutron unbound ground state resonances of $^{26}$O and $^{16}$Be were populated using one-proton knockout reactions from $^{27}$F and $^{17}$B beams. A coincidence measurement of 3-body system (fragment + n + n) allowed for the decay energy of the unbound nuclei to be reconstructed. A low energy resonance, $<$ 200 keV, was observed for the first time in the $^{24}$O + n + n system and assigned to the ground state of $^{26}$O. The $^{16}$Be ground state resonance was observed at 1.35 MeV. The 3-body correlations of the $^{14}$Be + n + n system were compared to simulations of a phase-space, sequential, and dineutron decay. The strong correlations in the n-n system from the experimental data could only be reproduced by the dineutron decay simulation providing the first evidence for a dineutron-like decay.

preprint2012arXiv

Hyperon sigma terms for 2+1 quark flavours

QCD lattice simulations determine hadron masses as functions of the quark masses. From the gradients of these masses and using the Feynman-Hellmann theorem the hadron sigma terms can then be determined. We use here a novel approach of keeping the singlet quark mass constant in our simulations which upon using an SU(3) flavour symmetry breaking expansion gives highly constrained (i.e. few parameter) fits for hadron masses in a multiplet. This is a highly advantageous procedure for determining the hadron mass gradient as it avoids the use of delicate chiral perturbation theory. We illustrate the procedure here by estimating the light and strange sigma terms for the baryon octet.

preprint2012arXiv

Observation of Ground-State Two-Neutron Decay

Neutron decay spectroscopy has become a successful tool to explore nuclear properties of nuclei with the largest neutron-to-proton ratios. Resonances in nuclei located beyond the neutron dripline are accessible by kinematic reconstruction of the decay products. The development of two-neutron detection capabilities of the Modular Neutron Array (MoNA) at NSCL has opened up the possibility to search for unbound nuclei which decay by the emission of two neutrons. Specifically this exotic decay mode was observed in 16Be and 26O.

preprint2012arXiv

Perturbative and Nonperturbative Renormalization in Lattice QCD

We investigate the perturbative and nonperturbative renormalization of composite operators in lattice QCD restricting ourselves to operators that are bilinear in the quark fields (quark-antiquark operators). These include operators which are relevant to the calculation of moments of hadronic structure functions. The nonperturbative computations are based on Monte Carlo simulations with two flavors of clover fermions and utilize the Rome-Southampton method also known as the RI-MOM scheme. We compare the results of this approach with various estimates from lattice perturbation theory, in particular with recent two-loop calculations.

preprint2012arXiv

Perturbative subtraction of lattice artifacts in the computation of renormalization constants

The determination of renormalization factors is of crucial importance. They relate the observables obtained on finite, discrete lattices to their measured counterparts in the continuum in a suitable renormalization scheme. Therefore, they have to be computed as precisely as possible. A widely used approach is the nonperturbative Rome-Southampton method. It requires, however, a careful treatment of lattice artifacts. They are always present because simulations are done at lattice spacings $a$ and momenta $p$ with $ap$ not necessarily small. In this paper we try to suppress these artifacts by subtraction of one-loop contributions in lattice perturbation theory. We compare results obtained from a complete one-loop subtraction with those calculated for a subtraction of $O(a^2)$.

preprint2012arXiv

Primary gamma-ray spectra in 44Ti of astrophysical interest

Primary gamma-ray spectra for a wide excitation-energy range have been extracted for 44Ti from particle-gamma coincidence data of the 46Ti(p,t gamma)44Ti reaction. These spectra reveal information on the gamma-decay pattern of the nucleus, and may be used to extract the level density and radiative strength function applying the Oslo method. Models of the level density and radiative strength function are used as input for cross-section calculations of the 40Ca(alpha,gamma)44Ti reaction. Acceptable models should reproduce data on the 40Ca(alpha,gamma)44Ti reaction cross section as well as the measured primary gamma-ray spectra. This is only achieved when a coherent normalization of the slope of the level density and radiative strength function is performed. Thus, the overall shape of the experimental primary gamma-ray spectra puts a constraint on the input models for the rate calculations.

preprint2012arXiv

Wilson loops to 20th order numerical stochastic perturbation theory

We calculate Wilson loops of various sizes up to 20 loops in SU(3) pure lattice gauge theory at different lattice sizes for Wilson gauge action using the technique of numerical stochastic perturbation theory. This allows us to investigate the perturbative series for various Wilson loops at high loop orders. We observe differences in the behavior of those series as function of the loop order. Up to $n=20$ we do not find evidence for the factorial growth of the expansion coefficients often assumed to characterize an asymptotic series. Based on the actually observed behavior we sum the series in a model parametrized by hypergeometric functions. Alternatively we estimate the total series in boosted perturbation theory using information from the first 14 loops. We introduce generalized ratios of Wilson loops of different sizes. Together with the corresponding Wilson loops from standard Monte Carlo measurements they enable us to assess their non-perturbative parts.

preprint2011arXiv

Flavour blindness and patterns of flavour symmetry breaking in lattice simulations of up, down and strange quarks

QCD lattice simulations with 2+1 flavours (when two quark flavours are mass degenerate) typically start at rather large up-down and strange quark masses and extrapolate first the strange quark mass and then the up-down quark mass to its respective physical value. Here we discuss an alternative method of tuning the quark masses, in which the singlet quark mass is kept fixed. Using group theory the possible quark mass polynomials for a Taylor expansion about the flavour symmetric line are found, first for the general 1+1+1 flavour case and then for the 2+1 flavour case. This ensures that the kaon always has mass less than the physical kaon mass. This method of tuning quark masses then enables highly constrained polynomial fits to be used in the extrapolation of hadron masses to their physical values. Numerical results for the 2+1 flavour case confirm the usefulness of this expansion and an extrapolation to the physical pion mass gives hadron mass values to within a few percent of their experimental values. Singlet quantities remain constant which allows the lattice spacing to be determined from hadron masses (without necessarily being at the physical point). Furthermore an extension of this programme to include partially quenched results is given.

preprint2011arXiv

Neutron knockout of 12Be populating neutron-unbound states in 11Be

Neutron-unbound resonant states of 11Be were populated in neutron knock-out reactions from 12Be and identified by 10Be-n coincidence measurements. A resonance in the decay-energy spectrum at 80(2) keV was attributed to a highly excited unbound state in 11Be at 3.949(2) MeV decaying to the 2+ excited state in 10Be. A knockout cross section of 15(3) mb was inferred for this 3.949(2) MeV state suggesting a spectroscopic factor near unity for this 0p3/2- level, consistent with the detailed shell model calculations.

preprint2011arXiv

Nucleon sigma terms for 2+1 quark flavours

QCD lattice simulations yield hadron masses as functions of the quark masses. From the gradients of the hadron masses the sigma terms can then be determined. We consider here dynamical 2+1 flavour simulations, in which we start from a point of the flavour symmetric line and then keep the singlet or average quark mass fixed as we approach the physical point. This leads to highly constrained fits for hadron masses in a multiplet. The gradient of this path for a hadron mass then gives a relation between the light and strange sigma terms. A further relation can be found from the change in the singlet quark mass along the flavour symmetric line. This enables light and strange sigma terms to be estimated for the baryon octet.

preprint2010arXiv

Exploring the Nucleon Structure from First Principles of QCD

Quantum Chromodynamics (QCD) is generally assumed to be the fundamental theory underlying nuclear physics. In recent years there is progress towards investigating the nucleon structure from first principles of QCD. Although this structure is best revealed in Deep Inelastic Scattering, a consistent analysis has to be performed in a fully non-perturbative scheme. The only known method for this purpose are lattice simulations. We first sketch the ideas of Monte Carlo simulations in lattice gauge theory. Then we comment in particular on the issues of chiral symmetry and operator mixing. Finally we present our results for the Bjorken variable of a single quark, and for the second Nachtmann moment of the nucleon structure functions.

preprint2010arXiv

Flavour symmetry breaking and tuning the strange quark mass for 2+1 quark flavours

QCD lattice simulations with 2+1 flavours typically start at rather large up-down and strange quark masses and extrapolate first the strange quark mass to its physical value and then the up-down quark mass. An alternative method of tuning the quark masses is discussed here in which the singlet quark mass is kept fixed, which ensures that the kaon always has mass less than the physical kaon mass. Using group theory the possible quark mass polynomials for a Taylor expansion about the flavour symmetric line are found, which enables highly constrained fits to be used in the extrapolation of hadrons to the physical pion mass. Numerical results confirm the usefulness of this expansion and an extrapolation to the physical pion mass gives hadron mass values to within a few percent of their experimental values.

preprint2010arXiv

Level densities and $γ$-ray strength functions in Sn isotopes

The nuclear level densities of $^{118,119}$Sn and the $γ$-ray strength functions of $^{116,118,119}$Sn below the neutron separation energy are extracted with the Oslo method using the ($^3$He, \,$αγ$) and ($^3$He,$^3$He$^\primeγ$) reactions. The level density function of $^{119}$Sn displays step-like structures. The microcanonical entropies are deduced from the level densities, and the single neutron entropy of $^{119}$Sn is determined to be $(1.7 \pm 0.2)\,k_B$. Results from a combinatorial model support the interpretation that some of the low-energy steps in the level density function are caused by neutron pair-breaking. An enhancement in all the $γ$-ray strength functions of $^{116-119}$Sn, compared to standard models for radiative strength, is observed for the $γ$-ray energy region of $\simeq (4 -11)$ MeV. These small resonances all have a centroid energy of 8.0(1) MeV and an integrated strength corresponding to $1.7(9)\%$ of the classical Thomas-Reiche-Kuhn sum rule. The Sn resonances may be due to electric dipole neutron skin oscillations or to an enhancement of the giant magnetic dipole resonance.

preprint2010arXiv

NSPT study of the three-loop lattice gluon propagator in Landau gauge

By means of Numerical Stochastic Perturbation Theory (NSPT), we calculate the lattice gluon propagator up to three loops of perturbation theory in the limits of infinite volume and vanishing lattice spacing. Based on known anomalous dimensions and a parametrization of both the hypercubic symmetry group H(4) and finite-size effects, we calculate the non-leading-log and non-logarithmic contributions iteratively, starting with the first-loop expression.

preprint2010arXiv

Renormalisation of composite operators in lattice QCD: perturbative versus nonperturbative

The perturbative and nonperturbative renormalisation of quark-antiquark operators in lattice QCD with two flavours of clover fermions is investigated within the research programme of the QCDSF collaboration. Operators with up to three derivatives are considered. The nonperturbative results based on the RI-MOM scheme are compared with estimates from one- and two-loop lattice perturbation theory.

preprint2010arXiv

Tuning the strange quark mass in lattice simulations

QCD lattice simulations with 2+1 flavours typically start at rather large up-down and strange quark masses and extrapolate first the strange quark mass to its physical value and then the up-down quark mass. An alternative method of tuning the quark masses is discussed here in which the singlet quark mass is kept fixed, which ensures that the kaon always has mass less than the physical kaon mass. It can also take into account the different renormalisations (for singlet and non-singlet quark masses) occurring for non-chirally invariant lattice fermions and so allows a smooth extrapolation to the physical quark masses. This procedure enables a wide range of quark masses to be probed, including the case with a heavy up-down quark mass and light strange quark mass. Results show the correct order for the baryon octet and decuplet spectrum and an extrapolation to the physical pion mass gives mass values to within a few percent of their experimental values.

preprint2010arXiv

Two-point functions of quenched lattice QCD in Numerical Stochastic Perturbation Theory

We summarize the higher-loop perturbative computation of the ghost and gluon propagators in SU(3) Lattice Gauge Theory. Our final aim is to compare with results from lattice simulations in order to expose the genuinely non-perturbative content of the latter. By means of Numerical Stochastic Perturbation Theory we compute the ghost and gluon propagators in Landau gauge up to three and four loops. We present results in the infinite volume and $a \to 0$ limits, based on a general fitting strategy.

preprint2010arXiv

Two-point functions of quenched lattice QCD in Numerical Stochastic Perturbation Theory. (II) The gluon propagator in Landau gauge

This is the second of two papers devoted to the perturbative computation of the ghost and gluon propagators in SU(3) Lattice Gauge Theory. Such a computation should enable a comparison with results from lattice simulations in order to reveal the genuinely non-perturbative content of the latter. The gluon propagator is computed by means of Numerical Stochastic Perturbation Theory: results range from two up to four loops, depending on the different lattice sizes. The non-logarithmic constants for one, two and three loops are extrapolated to the lattice spacing $a \to 0$ continuum and infinite volume $V \to \infty$ limits.

preprint2010arXiv

Very high order lattice perturbation theory for Wilson loops

We calculate perturbative Wilson loops of various sizes up to loop order $n=20$ at different lattice sizes for pure plaquette and tree-level improved Symanzik gauge theories using the technique of Numerical Stochastic Perturbation Theory. This allows us to investigate the behavior of the perturbative series at high orders. We observe differences in the behavior of perturbative coefficients as a function of the loop order. Up to $n=20$ we do not see evidence for the often assumed factorial growth of the coefficients. Based on the observed behavior we sum this series in a model with hypergeometric functions. Alternatively we estimate the series in boosted perturbation theory. Subtracting the estimated perturbative series for the average plaquette from the non-perturbative Monte Carlo result we estimate the gluon condensate.

preprint2009arXiv

A Non-Perturbative Operator Product Expansion

Nucleon structure functions can be observed in Deep Inelastic Scattering experiments, but it is an outstanding challenge to confront them with fully non-perturbative QCD results. For this purpose we investigate the product of electromagnetic currents (with large photon momenta) between quark states (of low momenta). By means of an Operator Product Expansion the structure function can be decomposed into matrix elements of local operators, and Wilson coefficients. For consistency both have to be computed non-perturbatively. Here we present precision results for a set of Wilson coefficients. They are evaluated from propagators for numerous quark momenta on the lattice, where the use of chiral fermions suppresses undesired operator mixing. This over-determines the Wilson coefficients, but reliable results can be extracted by means of a Singular Value Decomposition.

preprint2009arXiv

Gamma-strength functions in 60Ni from two-step cascades following proton capture

The two-step cascade method previously used in neutron capture experiments is now applied to a proton capture reaction. The spectrum of two-step cascades populating the first 2+ level of 60Ni has been measured with 59Co(p,2gamma)60Ni reaction. The simulation technique used for the spectrum analysis allows one to reveal the range of possible shapes of both E1 and M1 gamma-strength functions. The low-energy enhancement previously observed in 3He induced reactions is seen to appear in M1 strength functions of 60Ni.

preprint2009arXiv

Results from 2+1 flavours of SLiNC fermions

QCD results are presented for a 2+1 flavour fermion clover action (which we call the SLiNC action). A method of tuning the quark masses to their physical values is discussed. In this method the singlet quark mass is kept fixed, which solves the problem of different renormalisations (for singlet and non-singlet quark masses) occuring for non-chirally invariant lattice fermions. This procedure enables a wide range of quark masses to be probed, including the case with a heavy up-down quark mass and light strange quark mass. Preliminary results show the correct splittings for the baryon (octet and) decuplet spectrum.

preprint2009arXiv

The lattice ghost propagator in Landau gauge up to three loops using Numerical Stochastic Perturbation Theory

We complete our high-accuracy studies of the lattice ghost propagator in Landau gauge in Numerical Stochastic Perturbation Theory up to three loops. We present a systematic strategy which allows to extract with sufficient precision the non-logarithmic parts of logarithmically divergent quantities as a function of the propagator momentum squared in the infinite-volume and $a\to 0$ limits. We find accurate coincidence with the one-loop result for the ghost self-energy known from standard Lattice Perturbation Theory and improve our previous estimate for the two-loop constant contribution to the ghost self-energy in Landau gauge. Our results for the perturbative ghost propagator are compared with Monte Carlo measurements of the ghost propagator performed by the Berlin Humboldt university group which has used the exponential relation between potentials and gauge links.

preprint2009arXiv

Two-point functions of quenched lattice QCD in Numerical Stochastic Perturbation Theory. (I) The ghost propagator in Landau gauge

This is the first of a series of two papers on the perturbative computation of the ghost and gluon propagators in SU(3) Lattice Gauge Theory. Our final aim is to eventually compare with results from lattice simulations in order to enlight the genuinely non-perturbative content of the latter. By means of Numerical Stochastic Perturbation Theory we compute the ghost propagator in Landau gauge up to three loops. We present results in the infinite volume and $a \to 0$ limits, based on a general strategy that we discuss in detail.

preprint2009arXiv

Wilson loops in very high order lattice perturbation theory

We calculate Wilson loops of various sizes up to loop order $n=20$ for lattice sizes of $L^4 (L=4, 6, 8, 12)$ using the technique of Numerical Stochastic Perturbation Theory in quenched QCD. This allows to investigate the behaviour of the perturbative series at high orders. We discuss three models to estimate the perturbative series: a renormalon inspired fit, a heuristic fit based on an assumed power-law singularity and boosted perturbation theory. We have found differences in the behavior of the perturbative series for smaller and larger Wilson loops at moderate $n$. A factorial growth of the coefficients could not be confirmed up to $n=20$. From Monte Carlo measured plaquette data and our perturbative result we estimate a value of the gluon condensate $<\fracαπGG>$.

preprint2008arXiv

Clover improvement for stout-smeared 2+1 flavour SLiNC fermions: non-perturbative results

We discuss an action in which the fermion matrix has single level stout smearing for the hopping terms together with unsmeared links for the clover term. With the (tree level) Symanzik improved gluon action this constitutes the Stout Link Non-perturbative Clover or SLiNC action. To cancel O(a) terms the clover coefficient, csw, has to be tuned. We present here preliminary results of a non-perturbative determination of csw using the Schrodinger functional and as a by-product also a determination of the critical hopping parameter. A determination of the renormalisation constant for the local vector current is also given. Comparisons of the results are made with lowest order perturbation theory results.

preprint2008arXiv

Clover improvement for stout-smeared 2+1 flavour SLiNC fermions: perturbative results

For the Stout Link Non-perturbative Clover (SLiNC) action we determine in one-loop lattice perturbation theory the critical hopping parameter $κ_c$ and the clover parameter $c_{SW}$ which is needed for $\mathcal{O}(a)$ improvement. Performing this calculation off-shell we are also able to compute the non gauge invariant quark field improvement coefficient $c_{NGI}$. Additionally, we present first results for the renormalization factors of the scalar, pseudoscalar, vector and axial vector currents. We discuss mean field improvement for the SLiNC action.

preprint2008arXiv

Higher-loop gluon and ghost propagators in Landau gauge from numerical stochastic perturbation theory

We present higher loop results for the gluon and ghost propagator in Landau gauge on the lattice calculated in numerical stochastic perturbation theory. We make predictions for the perturbative content of those propagators as function of the lattice momenta for finite lattices. To find out their nonperturbative contributions, the logarithmic definition of the gauge fields and the corresponding Faddeev-Popov operator have to be implemented in the Monte Carlo simulations.

preprint2008arXiv

Renormalisation of Composite Operators in Lattice Perturbation Theory with Clover Fermions: Non-forward Matrix Elements

We consider the renormalisation of lattice QCD operators with one and two covariant derivatives related to the first and second moments of generalised parton distributions and meson distribution amplitudes. Employing the clover fermion action we calculate their non-forward quark matrix elements in one-loop lattice perturbation theory. For some representations of the hypercubic group commonly used in simulations we determine the sets of all possible mixing operators and compute the matrices of renormalisation factors in one-loop approximation. We describe how tadpole improvement is applied to the results.

preprint2006arXiv

Breaking of nucleon Cooper pairs at finite temperature in 93-98Mo

The S shape of the canonical heat capacity is known as a signature of the pairing transition and along an isotopic chain it is significantly more pronounced for nuclei with an even number of neutrons than with an odd number. Although the heat capacities extracted from experimental level densities in 93-98Mo exhibit a clear S shape, they do not show such an odd-even staggering. To understand the underlying physics, we analyze thermal quantities evaluated from the partition function calculated using the static-path plus random-phase approximation (SPA+RPA) in a monopole pairing model with number-parity projection. The calculated level densities reproduce very well the experimental data and they also agree with estimates using the back-shifted Fermi-gas model. We clarify the reason why the heat capacities for Mo isotopes do not show odd-even staggering of the S shape. We also discuss thermal odd-even mass differences in 94-97Mo using the three-, four-, and five-point formula. These thermal mass differences are regarded as indicators of pairing correlations at finite temperature.

preprint2006arXiv

Landau Gauge Gluon and Ghost Propagators from Lattice QCD

We report on recent numerical computations of the Landau gauge gluon and ghost propagators as well as of the ghost-gluon vertex function in pure SU(3) Yang-Mills theory and in full QCD on the lattice. Special emphasis is paid to the low momentum region. In particular, we present new data for the gluon propagator at momenta below 300 MeV. We also discuss different systematic effects as there are finite-size, lattice discretization and Gribov copy but also unquenching effects. A MOM-scheme running coupling α_s(q^2) based on the ghost-gluon vertex is calculated and found to decrease for momenta below 550 MeV, even though the renormalization constant of the vertex deviates only weakly from being constant.

preprint2006arXiv

Level density of $^{56}$Fe and low-energy enhancement of $γ$-strength function

The $^{55}$Mn$(d,n)^{56}$Fe differential cross section is measured at $E_d=7$ MeV\@. The $^{56}$Fe level density obtained from neutron evaporation spectra is compared to the level density extracted from the $^{57}$Fe$(^3$He,$αγ)^{56}$Fe reaction by the Oslo-type technique. Good agreement is found between the level densities determined by the two methods. With the level density function obtained from the neutron evaporation spectra, the $^{56}$Fe $γ$-strength function is also determined from the first-generation $γ$ matrix of the Oslo experiment. The good agreement between the past and present results for the $γ$-strength function supports the validity of both methods and is consistent with the low-energy enhancement of the $γ$ strength below $\sim 4$ MeV first discovered by the Oslo method in iron and molybdenum isotopes.

preprint2006arXiv

Microcanonical entropies and radiative strength functions of $^{50,51}$V

The level densities and radiative strength functions (RSFs) of $^{50,51}$V have been extracted using the ($^3$He,$αγ$) and ($^3$He,$^3$He$^{\prime} γ$) reactions, respectively. From the level densities, microcanonical entropies are deduced. The high $γ$-energy part of the RSF is described by the giant electric dipole resonance. A significant enhancement over the predicted strength in the region of $E_γ \lesssim 3$ MeV is seen, which at present has no theoretical explanation.

preprint2005arXiv

Perturbative renormalisation of quark bilinear operators for overlap fermions with and without stout links and improved gauge action

We calculate lattice renormalisation constants of local and one-link quark operators for overlap fermions and improved gauge actions in one-loop perturbation theory. For the local operators we stout smear the SU(3) links in the fermionic action. Using the popular tadpole improved Lüscher-Weisz actions at $β=8.45$ and $β=8.0$ we present numerical values for the Z factors in the $\bar{MS}$ scheme (partly as function of the stout smearing strength). We compare various levels of mean field (tadpole) improvement which have been applied to our results.

preprint2005arXiv

Studying the infrared region in Landau gauge QCD

We report on the progress we made in studying the infrared behavior of the ghost and gluon dressing functions in Landau gauge. Related to this we also investigate a running coupling given in terms of those functions and compare our results to those coming from the Dyson-Schwinger approach. We present first numerical results for the SU(3) ghost-ghost-gluon vertex renormalization constant. In addition the spectrum of low-lying eigenvalues and eigenfunctions of the Faddeev-Popov operator is determined. The saturation of the ghost propagator in terms of those eigenvalues and eigenmodes is discussed at lower momenta.

preprint2005arXiv

The compact Q=2 Abelian Higgs model in the London limit: vortex-monopole chains and the photon propagator

The confining and topological properties of the compact Abelian Higgs model with doubly-charged Higgs field in three space-time dimensions are studied. We consider the London limit of the model. We show that the monopoles are forming chain-like structures (kept together by ANO vortices) the presence of which is essential for getting simultaneously permanent confinement of singly-charged particles and breaking of the string spanned between doubly-charged particles. In the confinement phase the chains are forming percolating clusters while in the deconfinement (Higgs) phase the chains are of finite size. The described picture is in close analogy with the synthesis of the Abelian monopole and the center vortex pictures in confining non--Abelian gauge models. The screening properties of the vacuum are studied by means of the photon propagator in the Landau gauge.

preprint2004arXiv

Monopole Chains in the Compact Abelian Higgs Model with doubly-charged Matter Field

We study the properties of topological defects in the lattice compact Abelian Higgs Model with charge Q=2 matter field. We find that monopoles and antimonopoles form chain-like structures which are dense in the confinement/symmetric phase. In this phase the mentioned structures explain both the confinement of single-charged and the breaking of strings spanned between doubly-charged test particles. This observation helps to understand how the non-diagonal gluons, once taken into consideration in the Abelian projection of gluodynamics, could reproduce in this framework the string breaking for adjoint charges.

preprint2003arXiv

Decoherence in Disordered Conductors at Low Temperatures, the effect of Soft Local Excitations

The conduction electrons' dephasing rate, $τ_ϕ^{-1}$, is expected to vanish with the temperature. A very intriguing apparent saturation of this dephasing rate in several systems was recently reported at very low temperatures. The suggestion that this represents dephasing by zero-point fluctuations has generated both theoretical and experimental controversies. We start by proving that the dephasing rate must vanish at the $T\to 0$ limit, unless a large ground state degeneracy exists. This thermodynamic proof includes most systems of relevance and it is valid for any determination of $τ_ϕ$ from {\em linear} transport measurements. In fact, our experiments demonstrate unequivocally that indeed when strictly linear transport is used, the apparent low-temperature saturation of $τ_ϕ$ is eliminated. However, the conditions to be in the linear transport regime are more strict than hitherto expected. Another novel result of the experiments is that introducing heavy nonmagnetic impurities (gold) in our samples produces, even in linear transport, a shoulder in the dephasing rate at very low temperatures. We then show theoretically that low-lying local defects may produce a relatively large dephasing rate at low temperatures. However, as expected, this rate in fact vanishes when $T \to 0$, in agreement with our experimental observations.

preprint2003arXiv

New method for calculating helicity amplitudes of jet-like QED processes for high-energy colliders II. Processes with lepton pair production

As continuation of our previous paper we further develop our new method for calculating helicity amplitudes of jet-like QED processes described by tree diagrams, applying it to lepton pair production. This method consists in replacing spinor structures for real and weakly virtual intermediate leptons by simple transition vertices. New vertices are introduced for the pair production case, and previous bremsstrahlung vertices are generalized to include virtual photons inside the considered jet. We present a diagrammatic approach that allows to write down in an efficient way the leading helicity amplitudes, at tree level. The obtained compact amplitudes are particularly suitable for numerical calculations in jet-like kinematics. Several examples with up to four particles in a jet are discussed in detail.

preprint2002arXiv

String Breaking and Monopoles: a Case Study in the 3D Abelian Higgs Model

We study the breaking of the string spanned between test charges in the three dimensional Abelian Higgs model with compact gauge field and fundamentally charged Higgs field at zero temperature. In agreement with current expectations we demonstrate that string breaking is associated with pairing of monopoles. However, the string breaking is not accompanied by an ordinary phase transition.

preprint2001arXiv

Monopoles, confinement and deconfinement in lattice compact QED in (2+1)D with external fields

Finite temperature compact electrodynamics in (2+1) dimensions is studied in the presence of external electromagnetic fields. The deconfinement temperature is found to be insensitive to the external fields. This result corroborates our observation that external fields create additional small--size magnetic dipoles from the vacuum which do not spoil the confining properties of the model at low temperature. However, the Polyakov loop is not an order parameter of confinement. It can vanish in deconfinement in the presence of external field. This does not mean the restoration of confinement for certain external field fluxes. As a next step in the study of (2+1)D QED, the influence of monopoles on the photon propagator is studied. First results are presented showing this connection in the confining phase (without external field).

preprint2001arXiv

New method for calculating helicity amplitudes of jet--like QED processes for high--energy colliders I. Bremsstrahlung processes

Inelastic QED processes, the cross sections of which do not drop with increasing energy, play an important role at high-energy colliders. Such reactions have the form of two-jet processes with the exchange of a virtual photon in the t-channel. We consider them in the region of small scattering angles $m/E \lesssim θ\ll 1$, which yields the dominant contribution to their total cross sections. A new effective method is presented and applied to QED processes with emission of real photons to calculate the helicity amplitudes of these processes. Its basic idea is similar to the well-known equivalent-lepton method. Compact analytical expressions for those amplitudes up to $e^8$ are derived omitting only terms of the order of $m^2/E^2, θ^2$, $θm/E$ and higher order. The helicity amplitudes are presented in a compact form in which large compensating terms are already cancelled. Some common properties for all jet-like processes are found and we discuss their origin.

preprint1999arXiv

Search for fourth generation quarks and leptons at the Fermilab Tevatron and CERN Large Hadron Collider

If next generations of heavy quarks and leptons exist within the standard model (SM), they can manifest themselves in Higgs boson production at the Tevatron and the LHC, before being actually observed. This generation leads to an increase of the Higgs boson production cross section via gluon fusion at hadron colliders by a factor 6-9. So, the study of this process at the Tevatron and LHC can finally fix the number of generations in the SM. Using the $WW^*$ Higgs boson decay channel, the studies at the upgraded Tevatron will answer the question about the next generation for mass values 135 GeV $\lsim M_H\lsim 190$ GeV. Studying the $τ\barτ$ channel we show its large potential for the study of the Higgs boson at the LHC even in the standard case of three generations. At the Tevatron, studies in this channel could explore the mass range 110-140 GeV.

preprint1999arXiv

The visible effect of a very heavy magnetic monopole at colliders

If a heavy Dirac monopole exists, the light-to-light scattering below the monopole production threshold is enhanced due to strong coupling of monopoles to photons. At the next Linear Collider with electron beam energy 250 GeV this photon pair production could be observable at monopole masses less than 2.5-6.4 TeV in the $e^+e^-$ mode or 3.7-10 TeV in the $γγ$ mode, depending on the monopole spin. At the upgraded Tevatron such an effect is expected to be visible at monopole masses below 1-2.5 TeV. The strong dependence on the initial photon polarizations allows to find the monopole spin in experiments at $e^+e^-$ and $γγ$ colliders. We consider the $Zγ$ production and the $3γ$ production at $e^+e^-$ and $pp$ or $p\bar{p}$ colliders via the same monopole loop. The possibility to discover these processes is significantly lower than that of the $γγ$ case.

preprint1999arXiv

Vortex profiles and vortex interactions at the electroweak crossover

Local correlations of Z-vortex operators with gauge and Higgs fields (lattice quantum vortex profiles) as well as vortex two-point functions are studied in the crossover region near a Higgs mass of 100 GeV within the 3D SU(2) Higgs model. The vortex profiles resemble certain features of the classical vortex solutions in the continuum. The vortex-vortex interactions are analogous to the interactions of Abrikosov vortices in a type-I superconductor.

preprint1998arXiv

Embedded Topological Defects in Electroweak Theory: From Percolating Networks to Sphalerons

New defects (Z-vortices and Nambu monopoles) are found to become thermodynamically relevant for the broken phase near to the (weakly first order) electroweak phase transition, and below the crossover for higher Higgs mass. The symmetric phase is characterized by vortex condensation (percolation). The percolation transition persists in the crossover region. The quasiclassical nature of the vortices is supported by correlations with Higgs field and gauge field energy. Sphalerons are shown to be related to monopolium bound states.

preprint1998arXiv

Embedded Topological Defects in Hot Electroweak Theory: a Lattice Study

We study the properties of Nambu monopoles and Z-vortices in the 3D lattice SU(2) Higgs theory which represents the Standard Model at high temperature. We show that the densities of the Nambu monopoles and the Z-vortices are O(1) in the symmetric phase and generically small in the Higgs phase. Near to the critical Higgs mass and in the vicinity of the phase transition the densities are no more negligible in the broken phase. The percolation probability of the Z-vortex lines is found as a new disorder parameter for this phase transition. We conclude that the transition to the symmetric phase is accompanied by Z-vortex condensation. Simulations comparing elementary and extended vortices and monopoles at different β_G values, aiming to show that the density of vortices and monopoles of fixed physical size might have a well-defined continuum limit, gives encouraging but so far inconclusive results.

preprint1998arXiv

Wave Functions and Spectrum in Hot Electroweak Matter for Large Higgs Masses

We present results for the wave functions and the screening mass spectrum for quantum numbers $0^{++}$, $1^{--}$ and $2^{++}$ in the three-dimensional SU(2)-Higgs model near to the phase transition line below the endpoint and in the crossover region. Varying the 3D gauge couplings we study the behaviour along a line of constant physics towards the continuum limit in both phases. In the crossover region the changing spectrum of screening states versus temperature is examined showing the aftermath of the phase transition at lower Higgs mass. Different to smearing concepts we used large sets of operators with various extensions allowing to identify wave functions in position space.

preprint1998arXiv

Z-Vortex Percolation in the Electroweak Crossover Region

We study the statistical properties of Z-vortices and Nambu monopoles in the 3D SU(2) Higgs model for a Higgs mass M_H \approx 100 GeV near and above the crossover temperature, where these defects are thermally excited. Although there is no phase transition at that strong selfcoupling, we observe that the Z-vortices exhibit the percolation transition that has been found recently to accompany the first order thermal transition that exists at smaller Higgs mass. Above the crossover temperature percolating networks of Z-vortex lines are ubiquitous, whereas vortices form a dilute gas of closed vortex loops and (Nambu) monopolium states on the low-temperature side of the crossover. The percolation temperature turns out to be roughly independent of the lattice spacing. We find that the Higgs modulus is smaller (the gauge action is larger) inside the vortices, compared to the bulk average. This correlation becomes very strong on the low-temperature side. The percolation transition is a prerequisite of some string mediated baryon number generation scenarios.

preprint1997arXiv

Where the electroweak phase transition ends

We give a more precise characterisation of the end of the electroweak phase transition in the framework of the effective 3d SU(2)--Higgs lattice model than has been given before. The model has now been simulated at gauge couplings beta_G=12 and 16 for Higgs masses M_H^*=70, 74, 76 and 80 GeV up to lattices 96^3 and the data have been used for reweighting. The breakdown of finite volume scaling of the Lee-Yang zeroes indicates the change from a first order transition to a crossover at lambda_3/g_3^2=0.102(2) in rough agreement with results of Karsch et al (hep-lat/9608087) at β_G=9 and smaller lattices. The infinite volume extrapolation of the discontinuity Delta < phi^+ phi > /g_3^2 turns out to be zero at lambda_3/g_3^2=0.107(2) being an upper limit. We comment on the limitations of the second method.