Source author record

R. R. Horgan

R. R. Horgan 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

24works
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

24 published item(s)

preprint2015arXiv

Bottomonium hyperfine splittings from lattice NRQCD including radiative and relativistic corrections

We present a calculation of the hyperfine splittings in bottomonium using lattice Nonrelativistic QCD. The calculation includes spin-dependent relativistic corrections through O(v^6), radiative corrections to the leading spin-magnetic coupling and, for the first time, non-perturbative 4-quark interactions which enter at alpha_s^2 v^3. We also include the effect of u,d,s and c quark vacuum polarisation. Our result for the 1S hyperfine splitting is M(Upsilon,1S) - M(eta_b,1S)= 60.0(6.4) MeV. We find the ratio of 2S to 1S hyperfine splittings (M(Upsilon,2S) - M(eta_b,2S))/ (M(Upsilon,1S) - M(eta_b,1S)) = 0.445(28).

preprint2015arXiv

Radiative improvement of the lattice NRQCD action using the background field method and application to the hyperfine splitting of quarkonium states

We present the first application of the background field method to Non-Relativistic QCD (NRQCD) on the lattice in order to determine the one-loop radiative corrections to the coefficients of the NRQCD action in a manifestly gauge-covariant manner. The coefficient of the $σ\cdot B$ term in the NRQCD action is computed at the one-loop level; the resulting shift of the hyperfine splitting of bottomonium is found to bring the lattice predictions in line with experiment.

preprint2015arXiv

Radiative improvement of the lattice NRQCD action using the background field method with applications to quarkonium spectroscopy

We apply the background field (BF) method to Non-Relativistic QCD (NRQCD) on the lattice in order to determine the one-loop radiative corrections to the coefficients of the NRQCD action in a manifestly gauge-covariant manner by matching the NRQCD prediction for particular on-shell processes with those of relativistic continuum QCD. We explain how the BF method is implemented in automated perturbation theory and discuss the technique for matching the relativistic and non-relativistic theories. We compute the one-loop radiative corrections to the sigma.B and Darwin terms for the NRQCD action currently used in simulations, as well as the one-loop coefficients of the spin-dependent O(alpha^2) four-fermion contact terms. The effect of the corrections on the hyperfine splitting of bottomonium is estimated using earlier simulation results; the corrected lattice prediction is found to be in agreement with experiment. Agreement of the hyperfine splitting of bottomonium and the B-meson system is confirmed by recent simulation studies (Dowdall et al.) which include our NRQCD radiative corrections for the first time.

preprint2015arXiv

Rare $B$ decays using lattice QCD form factors

In this write-up we review and update our recent lattice QCD calculation of $B \to K^*$, $B_s \to ϕ$, and $B_s \to K^*$ form factors [arXiv:1310.3722]. These unquenched calculations, performed in the low-recoil kinematic regime, provide a significant improvement over the use of extrapolated light cone sum rule results. The fits presented here include further kinematic constraints and estimates of additional correlations between the different form factor shape parameters. We use these form factors along with Standard Model determinations of Wilson coefficients to give Standard Model predictions for several observables [arXiv:1310.3887]. The modest improvements to the form factor fits lead to improved determinations of $F_L$, the fraction of longitudinally polarized vector mesons, but have little effect on most other observables.

preprint2014arXiv

B-meson mixing from full lattice QCD with physical u, d, s and c quarks

We present the first lattice QCD calculation of the $B_s$ and $B_d$ mixing parameters with physical light quark masses. We use MILC gluon field configurations that include $u$, $d$, $s$ and $c$ sea quarks at 3 values of the lattice spacing and with 3 values of the $u/d$ quark mass going down to the physical value. We use improved NRQCD for the valence $b$ quarks. Preliminary results show significant improvements over earlier values.

preprint2013arXiv

B-meson decay constants from improved lattice NRQCD and physical u, d, s and c sea quarks

We present the first lattice QCD calculation of the decay constants f_B and f_{B_s} with physical light quark masses. We use configurations generated by the MILC collaboration including the effect of u, d, s and c HISQ sea quarks at three lattice spacings and with three u/d quark mass values going down to the physical value. We use improved NRQCD for the valence b quarks. Our results are f_B = 0.186(4) GeV, f_{B_s} = 0.224(5) GeV, f_{B_s}/f_B = 1.205(7) and M_{B_s}-M_B=85(2) MeV, superseding earlier results with NRQCD b quarks. We discuss the implications of our results for the Standard Model rates for B_s -> mu+ mu- and B -> tau nu.

preprint2013arXiv

B, Bs, K and pi weak matrix elements with physical light quarks

Calculations of pseudoscalar decay constants of B, Bs, K and pi mesons with physical light quarks are presented. We use HISQ ensembles that include u,d,s and c sea quarks at three lattice spacings. HISQ is used for the valence light quarks and a radiatively improved NRQCD action for the heavy quarks. The key results are f_{B^+}=0.184(4)$ GeV, f_{B_s}=0.224(4) GeV, f_{B_s}/f_{B^+}=1.217(8), f_{K^+}/f_{pi^+}=1.1916(21), f_{K^+}=155.37(34) MeV, giving a significant improvement over previous results that required chiral extrapolation. We also calculate the Wilson flow scale w_0, finding w_0=0.1715(9) fm.

preprint2013arXiv

The mass of the b-quark from lattice NRQCD and lattice perturbation theory

We present a determination of the b-quark mass accurate through O(α_s^2) in perturbation theory and including partial contributions at O(α_s^3). Nonperturbative input comes from the calculation of the Upsilon and B_s energies in lattice QCD including the effect of u, d and s sea quarks. We use an improved NRQCD action for the b-quark. This is combined with the heavy quark energy shift in NRQCD determined using a mixed approach of high-beta simulation and automated lattice perturbation theory. Comparison with experiment enables the quark mass to be extracted: in the MS bar scheme we find m_b(m_b) = 4.166(43) GeV.

preprint2012arXiv

Precise heavy-light meson masses and hyperfine splittings from lattice QCD including charm quarks in the sea

We present improved results for the B and D meson spectrum from lattice QCD including the effect of u/d,s and c quarks in the sea. For the B mesons the Highly Improved Staggered Quark action is used for the sea and light valence quarks and NonRelativistic QCD for the b quark including O(α_s) radiative corrections to many of the Wilson coefficients for the first time. The D mesons use the Highly Improved Staggered Quark action for both valence quarks on the same sea. We find M_{B_s}-M_B=84(2) MeV, M_{B_s}=5.366(8) GeV, M_{B_c}=6.278(9) GeV, M_{D_s}=1.9697(33) GeV, and M_{D_s}-M_{D}=101(3) MeV. Our results for the B meson hyperfine splittings are M_{B^*}-M_{B}=50(3) MeV, M_{B_s^*}-M_{B_s}=52(3) MeV, in good agreement with existing experimental results. This demonstrates that our perturbative improvement of the NRQCD chromo-magnetic coupling works for both heavyonium and heavy-light mesons. We predict M_{B_c^*}-M_{B_c}=54(3) MeV. We also present first results for the radially excited B_c states as well as the orbitally excited scalar B_c0^* and axial vector B_c1 mesons.

preprint2012arXiv

Radiative improvement of spin and Darwin terms in the NRQCD action

We present updated results for the radiative improvement of the σ.B term and the spin-dependent four-fermion terms in the lattice NRQCD action, and first results for the radiative corrections to the NRQCD Darwin term and spin-independent four-fermion terms. The spin-dependent terms have significant impact on getting the correct hyperfine splitting for both bottomonium and heavy-light mesons, while the spin-independent terms suffer from a conspiracy between lattice artifacts and severe IR divergences that complicates their evaluation.

preprint2011arXiv

Renormalization of heavy-light currents in moving NRQCD

Heavy-light decays such as $B \to π\ell ν$, $B \to K^{*} γ$ and $B \to K^{(*)} \ell \ell$ can be used to constrain the parameters of the Standard Model and in indirect searches for new physics. While the precision of experimental results has improved over the last years this has still to be matched by equally precise theoretical predictions. The calculation of heavy-light form factors is currently carried out in lattice QCD. Due to its small Compton wavelength we discretize the heavy quark in an effective non-relativistic theory. By formulating the theory in a moving frame of reference discretization errors in the final state are reduced at large recoil. Over the last years the formalism has been improved and tested extensively. Systematic uncertainties are reduced by renormalizing the m(oving)NRQCD action and heavy-light decay operators. The theory differs from QCD only for large loop momenta at the order of the lattice cutoff and the calculation can be carried out in perturbation theory as an expansion in the strong coupling constant. In this paper we calculate the one loop corrections to the heavy-light vector and tensor operator. Due to the complexity of the action the generation of lattice Feynman rules is automated and loop integrals are solved by the adaptive Monte Carlo integrator VEGAS. We discuss the infrared and ultraviolet divergences in the loop integrals both in the continuum and on the lattice. The light quarks are discretized in the ASQTad and highly improved staggered quark (HISQ) action; the formalism is easily extended to other quark actions.

preprint2011arXiv

Stochastic Processes, Slaves and Supersymmetry

We extend the work of Tanase-Nicola and Kurchan on the structure of diffusion processes and the associated supersymmetry algebra by examining the responses of a simple statistical system to external disturbances of various kinds. We consider both the stochastic differential equations (SDEs) for the process and the associated diffusion equation. The influence of the disturbances can be understood by augmenting the original SDE with an equation for {\it slave variables}. The evolution of the slave variables describes the behaviour of line elements carried along in the stochastic flow. These line elements together with the associated surface and volume elements constructed from them provide the basis of the supersymmetry properties of the theory. For ease of visualisation, and in order to emphasise a helpful electromagnetic analogy, we work in three dimensions. The results are all generalisable to higher dimensions and can be specialised to one and two dimensions. The electromagnetic analogy is a useful starting point for calculating asymptotic results at low temperature that can be compared with direct numerical evaluations. We also examine the problems that arise in a direct numerical simulation of the stochastic equation together with the slave equations. We pay special attention to the dependence of the slave variable statistics on temperature. We identify in specific models the critical temperature below which the slave variable distribution ceases to have a variance and consider the effect on estimates of susceptibilities.

preprint2011arXiv

The Upsilon spectrum and the determination of the lattice spacing from lattice QCD including charm quarks in the sea

We give results for the Upsilon spectrum from lattice QCD using an improved version of the NRQCD action for b quarks which includes radiative corrections to kinetic terms at O(v^4) in the velocity expansion. We also include for the first time the effect of up, down, strange and charm quarks in the sea using 'second generation' gluon field configurations from the MILC collaboration. Using the Υ 2S - 1S splitting to determine the lattice spacing, we are able to obtain the 1P - 1S splitting to 1.4% and the 3S - 1S splitting to 2.4%. Our improved result for M (Upsilon) - M (eta_b) is 70(9) MeV and we predict M(Upsilon) - M(eta_b') = 35(3) MeV. We also calculate pi, K and eta_s correlators using the Highly Improved Staggered Quark action and perform a chiral and continuum extrapolation to give values for M(eta_s) (0.6893(12) GeV) and f_{eta_s} (0.1819(5) GeV) that allow us to tune the strange quark mass as well as providing an independent and consistent determination of the lattice spacing. Combining the NRQCD and HISQ analyses gives mb/ms = 54.7(2.5) and a value for the heavy quark potential parameter of r_1 = 0.3209(26) fm.

preprint2010arXiv

Improved automated lattice perturbation theory in background field gauge

We present an algorithm to automatically derive Feynman rules for lattice perturbation theory in background field gauge. Vertices with an arbitrary number of both background and quantum legs can be derived automatically from both gluonic and fermionic actions. The algorithm is a generalisation of our earlier algorithm based on prior work by Lüscher and Weisz. We also present techniques allowing for the parallelisation of the evaluation of the often rather complex lattice Feynman rules that should allow for efficient implementation on GPUs, but also give a significant speed-up when calculating the derivatives of Feynman diagrams with respect to external momenta.

preprint2010arXiv

The b quark mass from lattice nonrelativistic QCD

We present the first two-loop calculation of the heavy quark energy shift in lattice nonrelativistic QCD (NRQCD). This calculation allow us to extract a preliminary prediction of $m_b(m_b, n_f = 5) = 4.25(12)$ GeV for the mass of the b quark from lattice NRQCD simulations performed with a lattice of spacing $a=0.12$fm. Our result is an improvement on a previous determination of the b quark mass from unquenched lattice NRQCD simulations, which was limited by the use of one-loop expressions for the energy shift. Our value is in good agreement with recent results of $m_b(m_b) = 4.163(16)$ GeV from QCD sum rules and $m_b(m_b, n_f = 5) = 4.170(25)$ GeV from realistic lattice simulations using highly-improved staggered quarks. We employ a mixed strategy to simplify our calculation. Ghost, gluon and counterterm contributions to the energy shift and mass renormalisation are extracted from quenched high-beta simulations whilst fermionic contributions are calculated using automated lattice perturbation theory. Our results demonstrate the effectiveness of such a strategy.

preprint2009arXiv

Automated generation of lattice QCD Feynman rules

The derivation of the Feynman rules for lattice perturbation theory from actions and operators is complicated, especially for highly improved actions such as HISQ. This task is, however, both important and particularly suitable for automation. We describe a suite of software to generate and evaluate Feynman rules for a wide range of lattice field theories with gluons and (relativistic and/or heavy) quarks. Our programs are capable of dealing with actions as complicated as (m)NRQCD and HISQ. Automated differentiation methods are used to calculate also the derivatives of Feynman diagrams.

preprint2009arXiv

Thermal Casimir effect between random layered dielectrics

We study the thermal Casimir effect between two thick slabs composed of plane-parallel layers of random dielectric materials interacting across an intervening homogeneous dielectric. It is found that the effective interaction at long distances is self averaging and is given by a description in terms of effective dielectric functions. The behavior at short distances becomes random (sample dependent) and is dominated by the local values of the dielectric function proximal to each other across the dielectrically homogeneous slab.

preprint2008arXiv

Perturbative calculations for the HISQ action: the gluon action at $O(N_fα_sa^2)$

We present a new (and general) algorithm for deriving lattice Feynman rules which is capable of handling actions as complex as the Highly Improved Staggered Quark (HISQ) action. This enables us to perform a perturbative calculation of the influence of dynamical HISQ fermions on the perturbative improvement of the gluonic action in the same way as we have previously done for asqtad fermions. We find the fermionic contributions to the radiative corrections in the Lüscher-Weisz gauge action to be somewhat larger for HISQ fermions than for asqtad.

preprint2003arXiv

Field theoretic calculation of the surface tension for a model electrolyte system

We carry out the calculation of the surface tension for a model electrolyte to first order in a cumulant expansion about a free field theory equivalent to the Debye-Hückel approximation. In contrast with previous calculations, the surface tension is calculated directly without recourse to integrating thermodynamic relations. The system considered is a monovalent electrolyte with a region at the interface, of width h, from which the ionic species are excluded. In the case where the external dielectric constant epsilon_0 is smaller than the electrolyte solution's dielectric constant epsilon we show that the calculation at this order can be fully regularized. In the case where h is taken to be zero the Onsager-Samaras limiting law for the excess surface tension of dilute electrolyte solutions is recovered, with corrections coming from a non-zero value of epsilon_0/epsilon.

preprint2003arXiv

The contribution of O(alpha) radiative corrections to the renormalised anisotropy and application to general tadpole improvement schemes: addendum to "One loop calculation of the renormalised anisotropy for improved anisotropic gluon actions on a lattice" [hep-lat/0208010]

General O(alpha) radiative corrections to lattice actions may be interpreted as counterterms that give additive contributions to the one-loop renormalisation of the anisotropy. The effect of changing the radiative coefficients is thus easily calculable. In particular, the results obtained in a previous paper for Landau mean link improved actions apply in any tadpole improvement scheme. We explain how this method can be exploited when tuning radiatively improved actions. Efficient methods for self-consistently tuning tadpole improvement factors are also discussed.

preprint2002arXiv

Semileptonic B Decays from an NRQCD/D234 Action

Semileptonic, B --> pi, l nubar, decays are studied on quenched anisotropic lattices using tree-level tadpole improved Symanzik glue, NRQCD heavy quark and D234 light quark actions. Constrained fitting methods are applied to extract groundstate contributions to two-point and three-point correlators. We agree with previous lattice determinations of the form factors. The major source of systematic error here, as in previous work, comes from the chiral extrapolation to the physical pion mass. Future calculations must work at lighter quark masses to resolve this.

preprint1997arXiv

Upsilon Spectrum from NRQCD with Improved Action

We explore the effect of higher order operators in the non-relativistic formulation of QCD (NRQCD). We calculated masses in the b \bar b spectrum using quenched gauge configurations at beta = 6.0 and two different NRQCD actions which have been corrected to order mv^4 and mv^6. The two-point functions are calculated in a gauge invariant fashion. We find the general structure of the spectrum to be the same in the two cases. Using the P - 3S1 splitting we determine the inverse lattice spacings to be 2.44(4) GeV and 2.44(5) GeV for the mv^4 and mv^6 actions, respectively. We do observe shifts in the spin splittings. The hyperfine splitting is reduced by approximately 4 MeV, while the fine splitting is down by up to 10 MeV, albeit with large statistical errors.

preprint1993arXiv

Radial Excited States for Heavy Quark Systems in NRQCD

Following the Non-Relativistic QCD approach we use a gauge invariant smearing method with factorization to measure the excitation energies for a heavy $Q\bar{Q}$ system on a $24^3\times 48$ lattice at $β=6.2$. The results come from averaging over an ensemble of 60 QCD configurations. In order to enhance the signal from each configuration we use wall sources for quark propagators. The quark Hamiltonian contains only the simplest non-relativistic kinetic energy term. The results are listed for a range of bare quark masses. The mass splittings are insensitive to this variable though there are a slight trends with increasing quark mass. For an appropriate choice of UV cut-off ($a^{-1}=3.2$Gev) the mass spectrum compares reasonably well with the experimental values for the spin-averaged energy gaps of the $Υ$ system. We also present results for the $DE$ and $DT$ waves for the lowest bare quark mass. The results are consistent with degeneracy between the two types of $D$ wave. This encourages the idea that even with our simple quark Hamiltonian the departure from rotational invariance is not great.