Researcher profile

Christopher J. Monahan

Christopher J. Monahan contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2022arXiv

A lattice QCD perspective on weak decays of b and c quarks Snowmass 2022 White Paper

Lattice quantum chromodynamics has proven to be an indispensable method to determine nonperturbative strong contributions to weak decay processes. In this white paper for the Snowmass community planning process we highlight achievements and future avenues of research for lattice calculations of weak $b$ and $c$ quark decays, and point out how these calculations will help to address the anomalies currently in the spotlight of the particle physics community. With future increases in computational resources and algorithmic improvements, percent level (and below) lattice determinations will play a central role in constraining the standard model or identifying new physics.

preprint2022arXiv

A novel nonperturbative renormalization scheme for local operators

The gradient flow exponentially suppresses ultraviolet field fluctuations and removes ultraviolet divergences (up to a multiplicative fermionic wavefunction renormalization). It can be used to describe real-space Wilsonian renormalization group transformations and determine the corresponding beta function. We propose a new nonperturbative renormalization scheme for local composite fermionic operators that uses the gradient flow and is amenable to lattice QCD calculations. We present preliminary nonperturbative results for the running of quark bilinear operators in this scheme and outline the calculation of perturbative matching to the MS-bar scheme.

preprint2022arXiv

Detailed analysis of excited state systematics in a lattice QCD calculation of $g_A$

Excited state contamination remains one of the most challenging sources of systematic uncertainty to control in lattice QCD calculations of nucleon matrix elements and form factors: early time separations are contaminated by excited states and late times suffer from an exponentially bad signal-to-noise problem. High-statistics calculations at large time separations $\gtrsim1$ fm are commonly used to combat these issues. In this work, focusing on $g_A$, we explore the alternative strategy of utilizing a large number of relatively low-statistics calculations at short to medium time separations (0.2--1 fm), combined with a multi-state analysis. On an ensemble with a pion mass of approximately 310 MeV and a lattice spacing of approximately 0.09 fm, we find this provides a more robust and economical method of quantifying and controlling the excited state systematic uncertainty. A quantitative separation of various types of excited states enables the identification of the transition matrix elements as the dominant contamination. The excited state contamination of the Feynman-Hellmann correlation function is found to reduce to the 1% level at approximately 1 fm while for the more standard three-point functions, this does not occur until after 2 fm. Critical to our findings is the use of a global minimization, rather than fixing the spectrum from the two-point functions and using them as input to the three-point analysis. We find that the ground state parameters determined in such a global analysis are stable against variations in the excited state model, the number of excited states, and the truncation of early-time or late-time numerical data.

preprint2022arXiv

One-loop matching for quark dipole operators in a gradient-flow scheme

The quark chromoelectric dipole (qCEDM) operator is a CP-violating operator describing, at hadronic energies, beyond-the-standard-model contributions to the electric dipole moment of particles with nonzero spin. In this paper we define renormalized dipole operators in a regularization-independent scheme using the gradient flow, and we perform the matching at one loop in perturbation theory to renormalized operators of the same and lower dimension in the more familiar MS scheme. We also determine the matching coefficients for the quark chromomagnetic dipole operator (qCMDM), which contributes, for example, to matrix elements relevant to CP-violating and CP-conserving kaon decays. The calculation provides a basis for future lattice QCD computations of hadronic matrix elements of the qCEDM and qCMDM operators.

preprint2021arXiv

The transversity parton distribution function of the nucleon using the pseudo-distribution approach

We present a determination of the non-singlet transversity parton distribution function (PDF) of the nucleon, normalized with respect to the tensor charge at $μ^2=2$ GeV$^2$ from lattice quantum chromodynamics. We apply the pseudo-distribution approach, using a gauge ensemble with a lattice spacing of 0.094 fm and the light quark mass tuned to a pion mass of 358 MeV. We extract the transversity PDF from the analysis of the short-distance behavior of the Ioffe-time pseudo-distribution using the leading-twist next-to-leading order (NLO) matching coefficients calculated for transversity. We reconstruct the $x$-dependence of the transversity PDF through an expansion in a basis of Jacobi polynomials in order to reduce the PDF ansatz dependence. Within the limitations imposed by a heavier-than-physical pion mass and a fixed lattice spacing, we present a comparison of our estimate for the valence transversity PDF with the recent global fit results based on single transverse spin asymmetry. We find the intrinsic nucleon sea to be isospin symmetric with respect to transversity.

preprint2020arXiv

$F_K / F_π$ from Möbius domain-wall fermions solved on gradient-flowed HISQ ensembles

We report the results of a lattice quantum chromodynamics calculation of $F_K/F_π$ using Möbius domain-wall fermions computed on gradient-flowed $N_f=2+1+1$ highly-improved staggered quark (HISQ) ensembles. The calculation is performed with five values of the pion mass ranging from $130 \lesssim m_π\lesssim 400$ MeV, four lattice spacings of $a\sim 0.15, 0.12, 0.09$ and $0.06$ fm and multiple values of the lattice volume. The interpolation/extrapolation to the physical pion and kaon mass point, the continuum, and infinite volume limits are performed with a variety of different extrapolation functions utilizing both the relevant mixed-action effective field theory expressions as well as discretization-enhanced continuum chiral perturbation theory formulas. We find that the $a\sim0.06$ fm ensemble is helpful, but not necessary to achieve a subpercent determination of $F_K/F_π$. We also include an estimate of the strong isospin breaking corrections and arrive at a final result of $F_{K^\pm}/F_{π^\pm} = 1.1942(45)$ with all sources of statistical and systematic uncertainty included. This is consistent with the Flavour Lattice Averaging Group average value, providing an important benchmark for our lattice action. Combining our result with experimental measurements of the pion and kaon leptonic decays leads to a determination of $|V_{us}|/|V_{ud}| = 0.2311(10)$.

preprint2020arXiv

Notes on lattice observables for parton distributions: nongauge theories

We review recent theoretical developments concerning the definition and the renormalization of equal-time correlators that can be computed on the lattice and related to Parton Distribution Functions (PDFs) through a factorization formula. We show how these objects can be studied and analyzed within the framework of a nongauge theory, gaining insight through a one-loop computation. We use scalar field theory as a playground to revise, analyze and present the main features of these ideas, to explore their potential, and to understand their limitations for extracting PDFs. We then propose a framework that would allow to include the available lattice QCD data in a global analysis to extract PDFs.

preprint2020arXiv

Short flow-time coefficients of CP-violating operators

Measurements of a permanent neutron electric dipole moment (EDM) potentially probe Beyond-the-Standard Model (BSM) sources of CP-violation. At low energy the CP-violating BSM interactions are parametrized by flavor-conserving CP-violating operators of dimension higher than four. QCD calculations of the nucleon matrix elements of these operators are required to fully reconstruct the sources and magnitudes of the different CP-violating contributions to the nucleon EDM. Herein we study the quark-chromo electric dipole moment (qCEDM) operator and the three-gluon Weinberg operator. The non-perturbative determination, using lattice QCD, of the nucleon matrix elements of these CP-violating operators is hampered by their short-distance behavior. Under renormalization these operators mix with lower dimensional operators, which induces power divergences in the lattice spacing, as the continuum limit is approached. We study the short-distance behavior of the qCEDM and the Weinberg operators using the gradient flow. We perform a short flow time expansion and determine, in perturbation theory, the expansion coefficients of the linearly-divergent terms stemming from the mixing with the pseudoscalar density and the topological charge, confirming the expectations of the operator product expansion. We introduce a new method to perform calculations at non-zero flow-time for arbitrary values of the external momenta. This method allows us to work in four dimensions for most of the calculations described in this paper, avoiding the complications associated with defining $γ_5$ in generic d dimensions. We show that leading contributions in the external momenta can be reproduced by defining $γ_5$ using the 't Hooft-Veltman-Breitenlohner-Maison scheme.

preprint2019arXiv

Lattice QCD Determination of $g_A$

The nucleon axial coupling, $g_A$, is a fundamental property of protons and neutrons, dictating the strength with which the weak axial current of the Standard Model couples to nucleons, and hence, the lifetime of a free neutron. The prominence of $g_A$ in nuclear physics has made it a benchmark quantity with which to calibrate lattice QCD calculations of nucleon structure and more complex calculations of electroweak matrix elements in one and few nucleon systems. There were a number of significant challenges in determining $g_A$, notably the notorious exponentially-bad signal-to-noise problem and the requirement for hundreds of thousands of stochastic samples, that rendered this goal more difficult to obtain than originally thought. I will describe the use of an unconventional computation method, coupled with "ludicrously'" fast GPU code, access to publicly available lattice QCD configurations from MILC and access to leadership computing that have allowed these challenges to be overcome resulting in a determination of $g_A$ with 1% precision and all sources of systematic uncertainty controlled. I will discuss the implications of these results for the convergence of $SU(2)$ Chiral Perturbation theory for nucleons, as well as prospects for further improvements to $g_A$ (sub-percent precision, for which we have preliminary results) which is part of a more comprehensive application of lattice QCD to nuclear physics. This is particularly exciting in light of the new CORAL supercomputers coming online, Sierra and Summit, for which our lattice QCD codes achieve a machine-to-machine speed up over Titan of an order of magnitude.

preprint2019arXiv

Lattice QCD matrix elements for the ${B_s^0-\bar{B}_s^0}$ width difference beyond leading order

Predicting the $B_s^0-\bar{B}_s^0$ width difference $ΔΓ_s$ relies on the heavy quark expansion and on hadronic matrix elements of $ΔB=2$ operators. We present the first lattice QCD results for matrix elements of the dimension-7 operators $R_{2,3}$ and linear combinations $\tilde{R}_{2,3}$ using nonrelativistic QCD for the bottom quark and a highly improved staggered quark (HISQ) action for the strange quark. Computations use MILC ensembles of gauge field configuations with $2+1+1$ flavors of sea quarks with the HISQ discretization, including lattices with physically light up/down quark masses. We discuss features unique to calculating matrix elements of these operators and analyze uncertainties from series truncation, discretization, and quark mass dependence. Finally we report the first Standard Model determination of $ΔΓ_s$ using lattice QCD results for all hadronic matrix elements through $\mathcal{O}(1/m_b)$. The main result of our calculations yields the $1/m_b$ contribution $ΔΓ_{1/m_b} = -0.022(10)~\mathrm{ps}^{-1}$. Adding this to the leading order contribution, the Standard Model prediction is $ΔΓ_s = 0.092(14)~\mathrm{ps}^{-1}$.