Source author record

Peter Boyle

Peter Boyle 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
11topics
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)

preprint2022arXiv

Discovering new physics in rare kaon decays

The decays and mixing of $K$ mesons are remarkably sensitive to the weak interactions of quarks and leptons at high energies. They provide important tests of the standard model at both first and second order in the Fermi constant $G_F$ and offer a window into possible new phenomena at energies as high as 1,000 TeV. These possibilities become even more compelling as the growing capabilities of lattice QCD make high-precision standard model predictions possible. Here we discuss and attempt to forecast some of these capabilities.

preprint2022arXiv

Grid: OneCode and FourAPIs

We discuss a substantial update to the Grid software library for Lattice QCD, enabling it to port to multiple GPU architectures while retaining CPU vectorisation and SIMD execution within OpenMP threads. The GPU environments supported include vendor specific Nvidia CUDA and AMD HIP environments and a (mostly) standards based SYCL implementation. This is performed by an internal abstraction interface giving single source cross-platform performance portability across all number of planned Exascale architectures, and all those planned by the US Department of Energy.

preprint2022arXiv

Lattice QCD and the Computational Frontier

The search for new physics requires a joint experimental and theoretical effort. Lattice QCD is already an essential tool for obtaining precise model-free theoretical predictions of the hadronic processes underlying many key experimental searches, such as those involving heavy flavor physics, the anomalous magnetic moment of the muon, nucleon-neutrino scattering, and rare, second-order electroweak processes. As experimental measurements become more precise over the next decade, lattice QCD will play an increasing role in providing the needed matching theoretical precision. Achieving the needed precision requires simulations with lattices with substantially increased resolution. As we push to finer lattice spacing we encounter an array of new challenges. They include algorithmic and software-engineering challenges, challenges in computer technology and design, and challenges in maintaining the necessary human resources. In this white paper we describe those challenges and discuss ways they are being dealt with. Overcoming them is key to supporting the community effort required to deliver the needed theoretical support for experiments in the coming decade.

preprint2022arXiv

Semileptonic $D \rightarrow π\ell ν$, $D \rightarrow K \ell ν$ and $D_s \rightarrow K \ell ν$ decays with 2+1f domain wall fermions

We present the status of our project to calculate $D \to π\ell ν$, $D \to K \ell ν$ and $D_s \to K \ell ν$ semileptonic form factors using domain wall fermions for both heavy and light quarks. Our computations are performed using RBC/UKQCD's set of 2+1 flavour domain wall fermion and Iwasaki gauge field ensembles. We plan to calculate three-point functions covering the full, physically allowed kinematic range. Given that the signal decays faster than the noise, unambiguously and reliably extracting the ground state is critical for success. We include an analysis of operator diagonalisation within several possible $2 \times 2$ operator bases and find an admixture of gauged fixed wall and $\mathbb{Z} \left( 2 \right)$ wall sources to be acceptable at both zero and non-zero momentum. Initial results for semileptonic form factors are presented for first ensembles.

preprint2021arXiv

Comparison of Domain Wall Fermion Multigrid Methods

We present a detailed comparison of several recent and new approaches to multigrid solver algorithms suitable for the solution of 5d chiral fermion actions such as Domain Wall fermions in the Shamir formulation, and also for the Partial Fraction and Continued Fraction overlap. Our focus is on the acceleration of gauge configuration sampling, and a compact nearest neighbour stencil is required to limit the calculational cost of obtaining a coarse operator. This necessitates the coarsening of a nearest neighbour operator to preserve sparsity in coarsened grids, unlike HDCG. We compare the approaches of HDCR and the Multigrid algorithm and also several new hybrid schemes. In this work we introduce a new recursive Chebyshev polynomial based setup scheme. We find that the HDCR approach, can both setup, and solve standard Shamir Domain Wall Fermions faster than a single solve with red-black preconditioned Conjugate Gradients on large volumes and for modern GPU systems such as the Summit supercomputer. This is promising for the acceleration of HMC, particularly if setup costs are shared across multiple Hasenbusch determinant factors. The setup scheme is likely generally applicable to other Fermion actions.

preprint2021arXiv

Proton decay matrix elements on the lattice at physical pion mass

Proton decay is a major prediction of Grand-Unified Theories (GUT) and its observation would indicate baryon number violation that is required for baryogenesis. Many decades of searching for proton decay have constrained its rate and ruled out some of the simplest GUT models. Apart from the baryon number-violating interactions, this rate also depends on transition amplitudes between the proton and mesons or leptons produced in the decay, which are matrix elements of three-quark operators. We report nonperturbative calculation of these matrix elements for the most studied two-body decay channels into a meson and antilepton done on a lattice with physical light and strange quark masses and lattice spacings $a\approx0.14$ and 0.20 fm. We perform nonperturbative renormalization and excited state analysis to control associated systematic effects. Our results largely agree with previous lattice calculations done with heavier quark masses and thus remove ambiguity in ruling out some simple GUT theories due to quark mass dependence of hadron structure.

preprint2020arXiv

Retinal vessel segmentation by probing adaptive to lighting variations

We introduce a novel method to extract the vessels in eye fun-dus images which is adaptive to lighting variations. In the Logarithmic Image Processing framework, a 3-segment probe detects the vessels by probing the topographic surface of an image from below. A map of contrasts between the probe and the image allows to detect the vessels by a threshold. In a lowly contrasted image, results show that our method better extract the vessels than another state-of the-art method. In a highly contrasted image database (DRIVE) with a reference , ours has an accuracy of 0.9454 which is similar or better than three state-of-the-art methods and below three others. The three best methods have a higher accuracy than a manual segmentation by another expert. Importantly, our method automatically adapts to the lighting conditions of the image acquisition.

preprint2019arXiv

Magritte, a modern software library for 3D radiative transfer: I. Non-LTE atomic and molecular line modelling

Radiative transfer is a key component in almost all astrophysical and cosmological simulations. We present Magritte: a modern open-source software library for 3D radiative transfer. It uses a deterministic ray-tracer and formal solver, i.e. it computes the radiation field by tracing rays through the model and solving the radiative transfer equation in its second-order form along a fixed set of rays originating from each point. Magritte can handle structured and unstructured input meshes, as well as smoothed-particle hydrodynamics (SPH) particle data. In this first paper, we describe the numerical implementation, semi-analytic tests and cross-code benchmarks for the non-LTE line radiative transfer module of Magritte. This module uses the radiative transfer solver to self-consistently determine the populations of the quantised energy levels of atoms and molecules using an accelerated Lambda iteration (ALI) scheme. We compare Magritte with the established radiative transfer solvers Ratran (1D) and Lime (3D) on the van Zadelhoff benchmark and present a first application to a simple Keplerian disc model. Comparing with Lime, we conclude that Magritte produces more accurate and more precise results, especially at high optical depth, and that it is faster.

preprint2016arXiv

A massive momentum-subtraction scheme

We introduce a new massive renormalization scheme, denoted mSMOM, as a modification of the existing RI/SMOM scheme. We use SMOM for defining renormalized fermion bilinears in QCD at non-vanishing fermion mass. This scheme has properties similar to those of the SMOM scheme, such as the use of non-exceptional symmetric momenta, while in contrast to SMOM, it defines the renormalized fields away from the chiral limit. Here we discuss some of the properties of mSMOM, and present non-perturbative arguments for deriving some renormalization constants. The results of a 1-loop calculation in dimensional regularization are briefly summarised to illustrate some properties of the scheme.

preprint2016arXiv

An exploratory study of heavy domain wall fermions on the lattice

We report on an exploratory study of domain wall fermions (DWF) as a lattice regularisation for heavy quarks. Within the framework of quenched QCD with the tree-level improved Symanzik gauge action we identify the DWF parameters which minimise discretisation effects. We find the corresponding effective 4$d$ overlap operator to be exponentially local, independent of the quark mass. We determine a maximum bare heavy quark mass of $am_h\approx 0.4$, below which the approximate chiral symmetry and O(a)-improvement of DWF are sustained. This threshold appears to be largely independent of the lattice spacing. Based on these findings, we carried out a detailed scaling study for the heavy-strange meson dispersion relation and decay constant on four ensembles with lattice spacings in the range $2.0-5.7\,\mathrm{GeV}$. We observe very mild $a^2$ scaling towards the continuum limit. Our findings establish a sound basis for heavy DWF in dynamical simulations of lattice QCD with relevance to Standard Model phenomenology.

preprint2016arXiv

Charm Physics with Domain Wall Fermions and Physical Pion Masses

We present RBC/UKQCD's charm project using $N_f=2+1$ flavour ensembles with inverse lattice spacings in the range $1.73-2.77\,\mathrm{GeV}$ and two physical pion mass ensembles. Domain wall fermions are used for the light as well as the charm quarks. We discuss our strategy for the extraction of the decay constants $f_D$ and $f_{D_s}$ and their extrapolation to the continuum limit, physical pion masses and the physical heavy quark mass. Our preliminary results are $f_D=208.7(2.8)\,\mathrm{MeV}$ and $f_{D_s}=246.4(1.9)\,\mathrm{MeV}$ where the quoted error is statistical only. We outline our current approach to extend the reach in the heavy quark mass and present preliminary results.

preprint2016arXiv

Electromagnetic Corrections to Meson Masses and the HVP

We present an exploratory study of the electromagnetic corrections to meson masses and the hadronic vacuum polarization using $N_f=2+1$ Domain Wall fermions. These corrections are estimated with two different approaches, a stochastic approach using $U(1)$ gauge configurations for the photon fields, and a perturbative approach through a QED perturbative expansion of the QCD+QED path integral. We compare results and statistical errors from both methods.

preprint2016arXiv

Hierarchically deflated conjugate residual

We present a progress report on a new class of multigrid solver algorithm suitable for the solution of 5d chiral fermions such as Domain Wall fermions and the Continued Fraction overlap. Unlike HDCG \cite{Boyle:2014rwa}, the algorithm works directly on a nearest neighbour fine operator. The fine operator used is Hermitian indefinite, for example $Γ_5 D_{dwf}$, and convergence is achieved with an indefinite matrix solver such as outer iteration based on conjugate residual. As a result coarse space representations of the operator remain nearest neighbour, giving an 8 point stencil rather than the 81 point stencil used in HDCG. It is hoped this may make it viable to recalculate the matrix elements of the little Dirac operator in an HMC evolution.

preprint2015arXiv

$N_f=2+1+1$ renormalisation of four-quark operators

When several four-quark operators are allowed to mix through renormalisation, this can considerably amplify the problems coming from perturbative truncation and discretisation effects. In this work we investigate whether our previous $B_K$ strategy can conveniently be generalised to a wider set of operators, corresponding to the so-called "SUSY $B_K$" basis of four-quark operators. We show that the discretisation effects, when plotted as a function of $ap$, are surprisingly independent of the lattice spacing. They appear reasonably under control up to very large energy scales. This allows us to discuss the effect of varying the intermediate scale on which the perturbative matching is done, and therefore the prospects of future high-precision studies with a Rome-Southampton renormalisation.

preprint2015arXiv

Charm physics with Moebius Domain Wall Fermions

We present results showing that Domain Wall fermions are a suitable discretisation for the simulation of heavy quarks. This is done by a continuum scaling study of charm quarks in a Möbius Domain Wall formalism using a quenched set-up. We find that discretisation effects remain well controlled by the choice of Domain Wall parameters preparing the ground work for the ongoing dynamical $2+1f$ charm program of RBC/UKQCD.

preprint2015arXiv

Domain Wall Charm Physics with Physical Pion Masses: Decay constants, bag and $ξ$ parameters

We provide an overview of RBC/UKQCD's charm project on 2+1 flavour physical pion mass ensembles using Möbius Domain Wall Fermions for the light as well as for the charm quark. We discuss the analysis strategy in detail and present results at the different stages of the analysis for $D$ and $D_s$ decay constants as well as the bag and $ξ$ parameters. We also discuss future approaches to extend the reach in the heavy quark mass.

preprint2015arXiv

Grid: A next generation data parallel C++ QCD library

In this proceedings we discuss the motivation, implementation details, and performance of a new physics code base called Grid. It is intended to be more performant, more general, but similar in spirit to QDP++\cite{QDP}. Our approach is to engineer the basic type system to be consistently fast, rather than bolt on a few optimised routines, and we are attempt to write all our optimised routines directly in the Grid framework. It is hoped this will deliver best known practice performance across the next generation of supercomputers, which will provide programming challenges to traditional scalar codes. We illustrate the programming patterns used to implement our goals, and advances in productivity that have been enabled by using new features in C++11.

preprint2014arXiv

Clover Action for Blue Gene-Q and Iterative solvers for DWF

In Lattice QCD, a major challenge in simulating physical quarks is the computational complexity of these simulations. In this proceeding, we describe the optimisation of Clover fermion action for Blue gene-Q architecture and how different iterative solvers behave for Domain Wall Fermion action. We find that the optimised Clover term achieved a maximum efficiency of 29.1% and 20.2% for single and double precision respectively for iterative Conjugate Gradient solver. For Domain Wall Fermion action (DWF) we found that Modified Conjugate Residual(MCR) as the most efficient solver compared to CG and GCR. We have developed a new multi-shift MCR algorithm that is 18.5% faster compared to multi-shift CG for the evaluation of rational functions in RHMC.

preprint2014arXiv

NPR step-scaling across the charm threshold

Matching Non-Perturbative Renormalisation on the lattice and perturbative renormalisation would benefit from higher matching scales, which are needed for observables entering their permile era such as BK. In this work we lay down a strategy, within the Rome-Southampton framework, to push this scale higher across the charm threshold, and apply it to an exploration of the BK running from 3 GeV to 9 GeV. This is done on Nf = 2+1+1 ensembles generated by the RBC-UKQCD collaboration, and features a close study of the discretisation effects.

preprint2011arXiv

Lattice Determination of the Hadronic Contribution to the Muon $g-2$ using Dynamical Domain Wall Fermions

We present a calculation of the leading order hadronic contribution to the anomalous magnetic moment of the muon for a dynamical simulation of 2+1 flavour QCD using domain wall fermions. The electromagnetic 2-point function is evaluated on the RBC-UKQCD lattice gauge configurations and this is fitted to a continuous form motivated by models of vector dominance. We determine a robust and reliable technique for performing this fit, allowing us to extract the most accurate results possible from our ensembles. This combined with data at very light quark masses produces the result a_μ^{(2)had}=641(33)(32) x 10^{-10} at the physical point, where the first uncertainty is statistical, and the second is an estimate of systematics, which is in agreement with previous results. We outline various methods by which this calculation can and will be improved in order to compete with the accuracy of alternative techniques of deducing this quantity from experimental scattering data.

preprint2011arXiv

Non-perturbative renormalization of kaon four-quark operators with nf=2+1 Domain Wall fermions

We present our strategy and some preliminary results for the renormalization of four-quark operators relevant for kaon physics. We follow the non-perturbative Rome-Southampton method, with both exceptional and non-exceptional kinematics. We also implement momentum sources and twisted boundary conditions. We use an (almost) unitary setup: Domain-Wall valence on nf=2+1 Domain-Wall sea and Iwasaki gauge action, at two values of the lattice spacing corresponding to approximately 0.086 fm and 0.114 fm . The chiral properties of these fermions play a crucial role in this compuation and are studied in detail in this work.

preprint2001arXiv

The Static Potential to O(α^2) in Lattice Perturbation Theory

We present a calculation of Wilson loops, and the static inter-quark potential to $O(α^2)$ in lattice perturbation theory. This is carried out with the Wilson, Symanzik-Weisz, and Iwasaki gauge actions and the Wilson, Sheikholeslami-Wohlert, and Kogut-Susskind dynamical fermion action for small Wilson loops, and with the Wilson gauge action and each of the dynamical quark actions in the case of the static potential.

preprint1998arXiv

A lattice potential investigation of quark mass and volume dependence of the $Υ$ spectrum

We investigate bottomonia splittings by solving a Schrodinger-Pauli-type equation with parametrisations of QCD potentials around those that have been determined previously in lattice simulations. This is done both, in the continuum and on finite lattices with resolutions ranging from a=0.2 fm down to a=0.025 fm and extent of up to 12 fm or 144^3 lattice points. We find a strong dependence of some splittings, in particular the 2S-1S and 1P-1S splittings, on both the quark mass and the short range form of the static potential in the neighbourhood of the bottom quark mass, while splittings such as 3S-2S and 2P-2S show reduced dependence on the short distance potential. We conclude that the quenched quarkonium spectrum cannot be matched to experiment with a simple redefinition of the lattice spacing. We investigate the size of relativistic corrections as a function of the quark mass. Finite size effects are shown to die out rather rapidly as the volume is increased, and we demonstrate the restoration of rotational symmetry as the continuum limit is taken.