Source author record

Thomas Blum

Thomas Blum 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

30works
8topics
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

30 published item(s)

preprint2022arXiv

$K\toππ$ decay matrix elements at the physical point with periodic boundary conditions

We calculate $K\toππ$ matrix elements using periodic boundary conditions as an independent calculation from our previous study with G-parity boundary conditions. We present our preliminary results for $K\toππ$ three-point functions and matrix elements on a $24^3, a^{-1} = 1$~GeV, $2+1$-flavor Möbius DWF ensemble at physical pion and kaon masses generated by the RBC and UKQCD collaborations and discuss the prospect for high-precision computation of $\varepsilon'$ with periodic boundary conditions.

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

Electric dipole moments and the search for new physics

Static electric dipole moments of nondegenerate systems probe mass scales for physics beyond the Standard Model well beyond those reached directly at high energy colliders. Discrimination between different physics models, however, requires complementary searches in atomic-molecular-and-optical, nuclear and particle physics. In this report, we discuss the current status and prospects in the near future for a compelling suite of such experiments, along with developments needed in the encompassing theoretical framework.

preprint2020arXiv

Light quark vacuum polarization at the physical point and contribution to the muon $g-2$

We report on the computation of the connected light quark vacuum polarization with 2+1+1 flavors of HISQ fermions at the physical point and its contribution to the muon anomalous magnetic moment. Three ensembles, generated by the MILC collaboration, are used to take the continuum limit. The finite volume correction to this result is computed in the (Euclidean) time-momentum representation to NNLO in chiral perturbation theory. We find $a_μ^{ll}({\rm HVP})=(659\pm 20\pm 5\pm 5\pm 4)\times 10^{-10}$, where the errors are statistical and estimates of residual uncertainties from taking the continuum limit, scale setting, and truncation of chiral perturbation theory at NNLO. We compare our results with recent ones in the literature.

preprint2020arXiv

Nucleon mass and isovector couplings in 2+1-flavor dynamical domain-wall lattice QCD near physical mass

We report nucleon mass, isovector vector and axial-vector charges, and tensor and scalar couplings, calculated using two recent 2+1-flavor dynamical domain-wall fermions lattice-QCD ensembles generated jointly by the RIKEN-BNL-Columbia and UKQCD collaborations. These ensembles were generated with Iwasaki $\times$ dislocation-suppressing-determinant-ratio gauge action at inverse lattice spacing of 1.378(7) GeV and pion mass values of 249.4(3) and 172.3(3) MeV. The nucleon mass extrapolates to a value $m_N = 0.950(5)$ GeV at physical point. The isovector vector charge renormalizes to unity in the chiral limit, narrowly constraining excited-state contamination in the calculation. The ratio of the isovector axial-vector to vector charges shows a deficit of about ten percent. The tensor coupling no longer depends on mass and extrapolates to 1.04(5) in $\overline {\rm MS}$ 2-GeV renormalization at physical point, in a good agreement with the value obtained at the lightest mass in our previous calculations and other calculations that followed. The scalar charge, though noisier, does not show mass dependence and is in agreement with other calculations.

preprint2020arXiv

The hadronic light-by-light scattering contribution to the muon anomalous magnetic moment from lattice QCD

We report the first result for the hadronic light-by-light scattering contribution to the muon anomalous magnetic moment with all errors systematically controlled. Several ensembles using 2+1 flavors of physical mass Möbius domain-wall fermions, generated by the RBC/UKQCD collaborations, are employed to take the continuum and infinite volume limits of finite volume lattice QED+QCD. We find $a_μ^{\rm HLbL} = 7.87(3.06)_\text{stat}(1.77)_\text{sys}\times 10^{-10}$. Our value is consistent with previous model results and leaves little room for this notoriously difficult hadronic contribution to explain the difference between the Standard Model and the BNL experiment.

preprint2020arXiv

Time-lapse reconstruction of the fracture front from diffracted waves arrivals in laboratory hydraulic fracture experiments

4D acoustic imaging via an array of 32 sources / 32 receivers is used to monitor hydraulic fracture propagating in a 250~mm cubic specimen under a true-triaxial state of stress. We present a method based on the arrivals of diffracted waves to reconstruct the fracture geometry (and fluid front when distinct from the fracture front). Using Bayesian model selection, we rank different possible fracture geometries (radial, elliptical, tilted or not) and estimate model error. The imaging is repeated every 4 seconds and provide a quantitative measurement of the growth of these low velocity fractures. We test the proposed method on two experiments performed in two different rocks (marble and gabbro) under experimental conditions characteristic respectively of the fluid lag-viscosity (marble) and toughness (gabbro) dominated hydraulic fracture propagation regimes. In both experiments, about 150 to 200 source-receiver combinations exhibit clear diffracted wave arrivals. The results of the inversion indicate a radial geometry evolving slightly into an ellipse towards the end of the experiment when the fractures feel the specimen boundaries. The estimated modelling error with all models is of the order of the wave arrival picking error. Posterior estimates indicate an uncertainty of the order of a millimeter on the fracture front location for a given acquisition sequence. The reconstructed fracture evolution from diffracted waves is shown to be consistent with the analysis of $90^{\circ}$ incidence transmitted waves across the growing fracture.

preprint2016arXiv

Finite-volume effects in the muon anomalous magnetic moment on the lattice

We investigate finite-volume effects in the hadronic vacuum polarization, with an eye toward the corresponding systematic error in the muon anomalous magnetic moment. We consider both recent lattice data as well as lowest-order, finite-volume chiral perturbation theory, in order to get a quantitative understanding. Even though leading-order chiral perturbation theory does not provide a good description of the hadronic vacuum polarization, it turns out that it gives a good representation of finite-volume effects. We find that finite-volume effects cannot be ignored when the aim is a few percent level accuracy for the leading-order hadronic contribution to the muon anomalous magnetic moment, even when using ensembles with $m_πL> 4$ and $m_π\sim 200$ MeV.

preprint2016arXiv

Project X: Physics Opportunities

Part 2 of "Project X: Accelerator Reference Design, Physics Opportunities, Broader Impacts". In this Part, we outline the particle-physics program that can be achieved with Project X, a staged superconducting linac for intensity-frontier particle physics. Topics include neutrino physics, kaon physics, muon physics, electric dipole moments, neutron-antineutron oscillations, new light particles, hadron structure, hadron spectroscopy, and lattice-QCD calculations. Part 1 is available as arXiv:1306.5022 [physics.acc-ph] and Part 3 is available as arXiv:1306.5024 [physics.acc-ph].

preprint2016arXiv

The connected and leading disconnected diagrams of the hadronic light-by-light contribution to muon $g - 2$

We report our recent lattice calculation of hadronic light-by-light contribution to muon $g-2$ using our recently developed moment method. The connected diagrams and the leading disconnected diagrams are included. The calculation is performed on a $48^3 \times 96$ lattice with physical pion mass and 5.5 fm box size. We expect sizable finite volume and finite lattice spacing corrections to the results of these calculations which will be estimated in calculations to be carried out over the next 1-2 years.

preprint2015arXiv

Covariant approximation averaging

We present a new class of statistical error reduction techniques for Monte-Carlo simulations. Using covariant symmetries, we show that correlation functions can be constructed from inexpensive approximations without introducing any systematic bias in the final result. We introduce a new class of covariant approximation averaging techniques, known as all-mode averaging (AMA), in which the approximation takes account of contributions of all eigenmodes through the inverse of the Dirac operator computed from the conjugate gradient method with a relaxed stopping condition. In this paper we compare the performance and computational cost of our new method with traditional methods using correlation functions and masses of the pion, nucleon, and vector meson in $N_f=2+1$ lattice QCD using domain-wall fermions. This comparison indicates that AMA significantly reduces statistical errors in Monte-Carlo calculations over conventional methods for the same cost.

preprint2015arXiv

Finite-volume effects in the hadronic vacuum polarization

We investigate finite-volume effects in the hadronic vacuum polarization, with an eye toward the corresponding systematic error in the muon anomalous magnetic moment. While it is well known that leading-order chiral perturbation theory does not provide a good description of the hadronic vacuum polarization, it turns out that it gives a much better representation of finite-volume effects. Indications are that finite-volume effects cannot be ignored when the aim is a few percent level accuracy for the hadronic contribution to the muon anomalous magnetic moment, even when $m_πL \sim 4$ and $m_π\sim 200$ MeV.

preprint2015arXiv

Hadronic Light by Light Contributions to the Muon Anomalous Magnetic Moment With Physical Pions

The current measurement of muonic $g - 2$ disagrees with the theoretical calculation by about 3 standard deviations. Hadronic vacuum polarization (HVP) and hadronic light by light (HLbL) are the two types of processes that contribute most to the theoretical uncertainty. The current value for HLbL is still given by models. I will describe results from a first-principles lattice calculation with a 139 MeV pion in a box of 5.5 fm extent. Our current numerical strategies, including noise reduction techniques, evaluating the HLbL amplitude at zero external momentum transfer, and important remaining challenges, in particular those associated with finite volume effects, will be discussed.

preprint2015arXiv

Lattice Calculation of Hadronic Light-by-Light Contribution to the Muon Anomalous Magnetic Moment

The quark-connected part of the hadronic light-by-light scattering contribution to the muon's anomalous magnetic moment is computed using lattice QCD with chiral fermions. We report several significant algorithmic improvements and demonstrate their effectiveness through specific calculations which show a reduction in statistical errors by more than an order of magnitude. The most realistic of these calculations is performed with a near-physical, $171$ MeV pion mass on a $(4.6\;\mathrm{fm})^3$ spatial volume using the $32^3\times 64$ Iwasaki+DSDR gauge ensemble of the RBC/UKQCD Collaboration.

preprint2015arXiv

Lattice Calculation of the Connected Hadronic Light-by-Light Contribution to the Muon Anomalous Magnetic Moment

The anomalous magnetic moment of muon, $g-2$, is a very precisely measured quantity. However, the current measurement disagrees with standard model by about 3 standard deviations. Hadronic vacuum polarization and hadronic light by light are the two types of processes that contribute most to the theoretical uncertainty. I will describe how lattice methods are well-suited to provide a first-principle's result for the hadronic light by light contribution, the various numerical strategies that are presently being used to evaluate it, our current results and the important remaining challenges which must be overcome.

preprint2015arXiv

Neutron and proton electric dipole moments from $N_f=2+1$ domain-wall fermion lattice QCD

We present a lattice calculation of the neutron and proton electric dipole moments (EDM's) with $N_f=2+1$ flavors of domain-wall fermions. The neutron and proton EDM form factors are extracted from three-point functions at the next-to-leading order in the $θ$ vacuum of QCD. In this computation, we use pion masses 0.33 and 0.42 GeV and 2.7 fm$^3$ lattices with Iwasaki gauge action and a 0.17 GeV pion and 4.6 fm$^3$ lattice with I-DSDR gauge action, all generated by the RBC and UKQCD collaborations. The all-mode-averaging technique enables an efficient and high statistics calculation. Chiral behavior of lattice EDM's is discussed in the context of baryon chiral perturbation theory. In addition, we also show numerical evidence on relationship of three- and two-point correlation function with local topological distribution.

preprint2015arXiv

On calculating disconnected-type hadronic light-by-light scattering diagrams from lattice QCD

For reliable comparison of the standard model prediction to the muon g-2 with its experimental value, the hadronic light-by-light scattering (HLbL) contribution must be calculated by lattice QCD simulation. HLbL contribution has many types of disconnected-type diagrams. Here, we start with recalling the point that must be taken care of in every method to calculate them by lattice QCD, and present one concrete method called nonperturbative QED method.

preprint2014arXiv

Hadronic light-by-light scattering contribution to the muon anomalous magnetic moment from lattice QCD

The form factor that yields the light-by-light scattering contribution to the muon anomalous magnetic moment is computed in lattice QCD+QED and QED. A non-perturbative treatment of QED is used and is checked against perturbation theory. The hadronic contribution is calculated for unphysical quark and muon masses, and only the diagram with a single quark loop is computed. Statistically significant signals are obtained. Initial results appear promising, and the prospect for a complete calculation with physical masses and controlled errors is discussed.

preprint2013arXiv

Non-degenerate light quark masses from 2+1f lattice QCD+QED

We report on a calculation of the effects of isospin breaking in Lattice QCD+QED. This involves using Chiral Perturbation Theory with Electromagnetic corrections to find the renormalized, non-degenerate, light quark masses. The calculations are carried out on QCD ensembles generated by the RBC and UKQCD collaborations using Domain Wall Fermions and the Iwasaki and Iwasaki+DSDR Gauge Actions with unitary pion masses down to 170 MeV. Non-compact QED is treated in the quenched approximation. The simulations use a $32^3$ lattice size with $a^{-1}=2.28(3)$ GeV (Iwasaki) and 1.37(1) (Iwasaki+DSDR). This builds on previous work from the RBC/UKQCD collaboration with lattice spacing $a^{-1}=1.78(4)$ GeV.

preprint2013arXiv

The hadronic vacuum polarization with twisted boundary conditions

The leading-order hadronic contribution to the anomalous magnetic moment of the muon is given by a weighted integral over the subtracted hadronic vacuum polarization. This integral is dominated by euclidean momenta of order the muon mass which are not available on current lattice volumes with periodic boundary conditions. Twisted boundary conditions can in principle help access momenta of any size even in a finite volume. We investigate the implementation of twisted boundary conditions both numerically (using all-mode averaging for improved statistics) and analytically, and present our initial results.

preprint2013arXiv

The Muon (g-2) Theory Value: Present and Future

This White Paper briefly reviews the present status of the muon (g-2) Standard-Model prediction. This value results in a 3 - 4 standard-deviation difference with the experimental result from Brookhaven E821. The present experimental uncertainty is $\pm 63 \times 10^{-11}$ (0.54~ppm), and the Standard-Model uncertainty is $\simeq \pm 49 \times 10^{-11}$. Fermilab experiment E989 has the goal to reduce the experimental error to $\pm 16 \times 10^{-11}$. Improvements in the Standard-Model value, which should be achieved between now and when the first results from Fermilab E989 could be available, should lead to a Standard-Model uncertainty of $\sim \,\pm 35 \times 10^{-11}$. These improvements would halve the uncertainty on the difference between experiment and theory, and should clarify whether the current difference points toward New Physics, or to a statistical fluctuation. At present, the (g-2) result is arguably the most compelling indicator of physics beyond the Standard Model and, at the very least, it represents a major constraint for speculative new theories such as supersymmetry, dark gauge bosons or extra dimensions.

preprint2012arXiv

A new class of variance reduction techniques using lattice symmetries

We present a general class of unbiased improved estimators for physical observables in lattice gauge theory computations which significantly reduces statistical errors at modest computational cost. The error reduction techniques, referred to as covariant approximation averaging, utilize approximations which are covariant under lattice symmetry transformations. We observed cost reductions from the new method compared to the traditional one, for fixed statistical error, of 16 times for the nucleon mass at $M_π\sim 330$ MeV (Domain-Wall quark) and 2.6-20 times for the hadronic vacuum polarization at $M_π\sim 480$ MeV (Asqtad quark). These cost reductions should improve with decreasing quark mass and increasing lattice sizes.

preprint2012arXiv

Error reduction technique using covariant approximation and application to nucleon form factor

We demonstrate the new class of variance reduction techniques for hadron propagator and nucleon isovector form factor in the realistic lattice of $N_f=2+1$ domain-wall fermion. All-mode averaging (AMA) is one of the powerful tools to reduce the statistical noise effectively for wider varieties of observables compared to existing techniques such as low-mode averaging (LMA). We adopt this technique to hadron two-point functions and three-point functions, and compare with LMA and traditional source-shift method in the same ensembles. We observe AMA is much more cost effective in reducing statistical error for these observables.

preprint2012arXiv

Full QED+QCD Low-Energy Constants through Reweighting

The effect of sea quark electromagnetic charge on meson masses is investigated, and first results for full QED+QCD low-energy constants are presented. The electromagnetic charge for sea quarks is incorporated in quenched QED+full QCD lattice simulations by a reweighting method. The reweighting factor, which connects quenched and unquenched QED, is estimated using a stochastic method on 2+1 flavor dynamical domain-wall quark ensembles.

preprint2012arXiv

Model-independent parametrization of the hadronic vacuum polarization and g-2 for the muon on the lattice

The leading hadronic contribution to the muon anomalous magnetic moment is given by a weighted integral over euclidean momentum of the hadronic vacuum polarization. This integral is dominated by momenta of order the muon mass. Since the finite volume in lattice QCD makes it difficult to compute the vacuum polarization at a large number of low momenta with high statistics (combined with the fact that one cannot compute it at zero momentum), a parametrization of the vacuum polarization is required to extrapolate the data. A much used functional form is based on vector meson dominance, which introduces model dependence into the lattice computation of the magnetic moment. Here we introduce a model-independent extrapolation method, and present a few first tests of this new method.

preprint2012arXiv

Pade approximants and g-2 for the muon

The leading hadronic contribution to the muon anomalous magnetic moment is given by a weighted euclidean momentum integral of the hadronic vacuum polarization. This integral is dominated by momenta of order the muon mass. Since in lattice QCD it is difficult to compute the vacuum polarization at a large number of low momenta, a parametrization of the vacuum polarization is required to extrapolate the data. Most fits to date are based on vector meson dominance, which introduces model dependence into the lattice computation of the magnetic moment. Here we introduce a model-independent extrapolation method, and present a few first tests of this new method.

preprint2009arXiv

Controlling Residual Chiral Symmetry Breaking in Domain Wall Fermion Simulations

At stronger gauge-field couplings, the domain wall fermion (DWF) residual mass, a measure of chiral symmetry breaking, grows rapidly. This measure is largely due to near zero fermion eigenmodes of logarithm of the 4D transfer matrix along the fifth dimension, and these eigenmodes increase rapidly at strong coupling. To suppress these eigenmodes, we have added to the DWF path integral a multiplicative weighting factor consisting of a ratio of determinants of Wilson-Dirac fermions having a chirally twisted mass with a large negative real component and a small imaginary chiral component. Numerical results show that this weighting factor with an appropriate choice of twisted masses significantly suppresses the residual mass while allowing adequate topological tunneling.

preprint2005arXiv

Nucleon structure functions with domain wall fermions

We present a quenched lattice QCD calculation of the first few moments of the polarized and un-polarized structure functions of the nucleon. Our calculations are done using domain wall fermions and the DBW2 gauge action with inverse lattice spacing ~1.3GeV, physical volume approximatelly (2.4 fm)^3, and light quark masses down to about 1/4 the strange quark mass. Values of the individual moments are found to be significantly larger than experiment, as in past lattice calculations, but interestingly the chiral symmetry of domain wall fermions allows for a precise determination of the ratio of the flavor non-singlet momentum fraction to the helicity distribution, which is in very good agreement with experiment. We discuss the implications of this result. Next, we show that the chiral symmetry of domain wall fermions is useful in eliminating mixing of power divergent lower dimensional operators with twist-3 operators. Finally, we find the isovector tensor charge at renormalization scale 2 GeV in the MS bar scheme to be 1.192(30), where the error is the statistical error only.