Source author record

Martin J. Savage

Martin J. Savage 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

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

65 published item(s)

preprint2023arXiv

Multi-Body Entanglement and Information Rearrangement in Nuclear Many-Body Systems

We examine how effective-model-space (EMS) calculations of nuclear many-body systems rearrange and converge multi-particle entanglement. The generalized Lipkin-Meshkov-Glick (LMG) model is used to motivate and provide insight for future developments of entanglement-driven descriptions of nuclei. The effective approach is based on a truncation of the Hilbert space together with a variational rotation of the qubits (spins), which constitute the relevant elementary degrees of freedom. The non-commutivity of the rotation and truncation allows for an exponential improvement of the energy convergence throughout much of the model space. Our analysis examines measures of correlations and entanglement, and quantifies their convergence with increasing cut-off. We focus on one- and two-spin entanglement entropies, mutual information, and $n$-tangles for $n=2,4$ to estimate multi-body entanglement. The effective description strongly suppresses entropies and mutual information of the rotated spins, while being able to recover the exact results to a large extent with low cut-offs. Naive truncations of the bare Hamiltonian, on the other hand, artificially underestimate these measures. The $n$-tangles in the present model provide a basis-independent measures of $n$-particle entanglement. While these are more difficult to capture with the EMS description, the improvement in convergence, compared to truncations of the bare Hamiltonian, is significantly more dramatic. We conclude that the low-energy EMS techniques, that successfully provide predictive capabilities for low-lying observables in many-body systems, exhibit analogous efficacy for quantum correlations and multi-body entanglement in the LMG model, motivating future studies in nuclear many-body systems and effective field theories relevant to high-energy physics and nuclear physics.

preprint2023arXiv

Quantum Simulations of SO(5) Many-Fermion Systems using Qudits

The structure and dynamics of quantum many-body systems are the result of a delicate interplay between underlying interactions, which leads to intricate entanglement structures. Despite this apparent complexity, symmetries emerge and have long been used to determine the relevant degrees of freedom and simplify classical descriptions of these systems. In this work, we explore the potential utility of quantum computers with arrays of qudits in simulating interacting fermionic systems, when the qudits can naturally map these relevant degrees of freedom. The Agassi model of fermions is based on an underlying $so(5)$ algebra, and the systems it describes can be partitioned into pairs of modes with five basis states, which naturally embed in arrays of $d=5$ qudits (qu5its). Classical noiseless simulations of the time evolution of systems of fermions embedded in up to twelve qu5its are performed using Google's cirq software. The resource requirements of the qu5it circuits are analyzed and compared with two different mappings to qubit systems, a physics-aware Jordan-Wigner mapping and a state-to-state mapping. We find advantages in using qudits, specifically in lowering the required quantum resources and reducing anticipated errors that take the simulation out of the physical space. A previously unrecognized sign problem has been identified from Trotterization errors in time evolving high-energy excitations. This has implications for quantum simulations in high-energy and nuclear physics, specifically of fragmentation and highly inelastic, multi-channel processes.

preprint2022arXiv

Snowmass Computational Frontier: Topical Group Report on Quantum Computing

Quantum computing will play a pivotal role in the High Energy Physics (HEP) science program over the early parts of the 21$^{st}$ Century, both as a major expansion of our capabilities across the Computational Frontier, and in synthesis with quantum sensing and quantum networks. This report outlines how Quantum Information Science (QIS) and HEP are deeply intertwined endeavors that benefit enormously from a strong engagement together. Quantum computers do not represent a detour for HEP, rather they are set to become an integral part of our discovery toolkit. Problems ranging from simulating quantum field theories, to fully leveraging the most sensitive sensor suites for new particle searches, and even data analysis will run into limiting bottlenecks if constrained to our current computing paradigms. Easy access to quantum computers is needed to build a deeper understanding of these opportunities. In turn, HEP brings crucial expertise to the national quantum ecosystem in quantum domain knowledge, superconducting technology, cryogenic and fast microelectronics, and massive-scale project management. The role of quantum technologies across the entire economy is expected to grow rapidly over the next decade, so it is important to establish the role of HEP in the efforts surrounding QIS. Fully delivering on the promise of quantum technologies in the HEP science program requires robust support. It is important to both invest in the co-design opportunities afforded by the broader quantum computing ecosystem and leverage HEP strengths with the goal of designing quantum computers tailored to HEP science.

preprint2021arXiv

Hierarchical Qubit Maps and Hierarchical Quantum Error Correction

We consider hierarchically implemented quantum error correction (HI-QEC), in which the fidelities of logical qubits are differentially optimized to enhance the capabilities of quantum devices in scientific applications. By employing qubit representations that propagate hierarchies in simulated systems to those in logical qubit noise sensitivities, heterogeneity in the distribution of physical-to-logical qubits can be systematically structured. For concreteness, we estimate HI-QEC's impact on surface code resources in computing low-energy observables to fixed precision, finding up to $\sim 60\%$ reductions in qubit requirements plausible in early error corrected simulations. Hierarchical qubit maps are also possible without error correction in qubit and qudit systems where fidelities are non-uniform, either unintentionally or by design. Hierarchical optimizations are another element in the co-design process of quantum simulations for nuclear and particle physics.

preprint2021arXiv

Standard Model Physics and the Digital Quantum Revolution: Thoughts about the Interface

Advances in isolating, controlling and entangling quantum systems are transforming what was once a curious feature of quantum mechanics into a vehicle for disruptive scientific and technological progress. Pursuing the vision articulated by Feynman, a concerted effort across many areas of research and development is introducing prototypical digital quantum devices into the computing ecosystem available to domain scientists. Through interactions with these early quantum devices, the abstract vision of exploring classically-intractable quantum systems is evolving toward becoming a tangible reality. Beyond catalyzing these technological advances, entanglement is enabling parallel progress as a diagnostic for quantum correlations and as an organizational tool, both guiding improved understanding of quantum many-body systems and quantum field theories defining and emerging from the Standard Model. From the perspective of three domain science theorists, this article compiles thoughts about the interface on entanglement, complexity, and quantum simulation in an effort to contextualize recent NISQ-era progress with the scientific objectives of nuclear and high-energy physics.

preprint2020arXiv

Systematically Localizable Operators for Quantum Simulations of Quantum Field Theories

Correlations and measures of entanglement in ground state wavefunctions of relativistic quantum field theories are spatially localized over length scales set by the mass of the lightest particle. We utilize this localization to design digital quantum circuits for preparing the ground states of lattice scalar quantum field theories. Controlled rotations that are exponentially localized in their position-space extent are found to provide exponentially convergent wavefunction fidelity. These angles scale with the correlation between sites and the classical two-point correlation function, as opposed to the more localized mutual information or the hyper-localized negativity. We anticipate that further investigations will uncover quantum circuit designs with controlled rotations dictated by the measures of entanglement. This work is expected to impact quantum simulations of systems of importance to nuclear physics, high-energy physics, and basic energy sciences research.

preprint2020arXiv

Two Nucleon Systems at $m_π\sim 450~{\rm MeV}$ from Lattice QCD

Nucleon-nucleon systems are studied with lattice quantum chromodynamics at a pion mass of $m_π\sim 450~{\rm MeV}$ in three spatial volumes using $n_f=2+1$ flavors of light quarks. At the quark masses employed in this work, the deuteron binding energy is calculated to be $B_d = 14.4^{+3.2}_{-2.6} ~{\rm MeV}$, while the dineutron is bound by $B_{nn} = 12.5^{+3.0}_{-5.0}~{\rm MeV}$. Over the range of energies that are studied, the S-wave scattering phase shifts calculated in the 1S0 and 3S1-3D1 channels are found to be similar to those in nature, and indicate repulsive short-range components of the interactions, consistent with phenomenological nucleon-nucleon interactions. In both channels, the phase shifts are determined at three energies that lie within the radius of convergence of the effective range expansion, allowing for constraints to be placed on the inverse scattering lengths and effective ranges. The extracted phase shifts allow for matching to nuclear effective field theories, from which low energy counterterms are extracted and issues of convergence are investigated. As part of the analysis, a detailed investigation of the single hadron sector is performed, enabling a precise determination of the violation of the Gell-Mann--Okubo mass relation. [An Erratum to the published version is included as an appendix. It details the impact of an error discovered in 2020 and corrects typographical errors.]

preprint2019arXiv

Minimally-Entangled State Preparation of Localized Wavefunctions on Quantum Computers

Initializing a single site of a lattice scalar field theory into an arbitrary state with support throughout the quantum register requires ${\cal O}(2^{n_Q})$ entangling gates on a quantum computer with $n_Q$ qubits per site. It is conceivable that, instead, initializing to functions that are good approximations to states may have utility in reducing the number of required entangling gates. In the case of a single site of a non-interacting scalar field theory, initializing to a symmetric exponential wavefunction requires $n_Q-1$ entangling gates, compared with the $2^{n_Q-1} + n_Q-3 + δ_{n_Q,1}$ required for a symmetric Gaussian wavefunction. In this work, we explore the initialization of 1-site ($n_Q=4$), 2-site ($n_Q=3$) and 3-site ($n_Q=3$) non-interacting scalar field theories with symmetric exponential wavefunctions using IBM's quantum simulators and quantum devices (Poughkeepsie and Tokyo). With the digitizations attainable with $n_Q = 3,4$, these tensor-product wavefunctions are found to have large overlap with a Gaussian wavefunction, and provide a suitable low-noise initialization for improvement and \emph{Somma Inflation}. In performing these simulations, we have employed a workflow that interleaves calibrations to mitigate systematic errors in production. The calibrations allow tolerance cuts on gate performance including the fidelity of the symmetrizing Hadamard gate, both in vacuum ($|{\bf 0}\rangle^{\otimes n_Q}$) and in medium ($n_Q-1$ qubits initialized to an exponential function). The results obtained in this work are relevant to systems beyond scalar field theories, such as the deuteron radial wavefunction, 2- and 3-dimensional cartesian-space wavefunctions, and non-relativistic multi-nucleon systems built on a localized eigenbasis.

preprint2019arXiv

SU(2) non-Abelian gauge field theory in one dimension on digital quantum computers

An improved mapping of one-dimensional SU(2) non-Abelian gauge theory onto qubit degrees of freedom is presented. This new mapping allows for a reduced unphysical Hilbert space. Insensitivity to interactions within this unphysical space is exploited to design more efficient quantum circuits. Local gauge symmetry is used to analytically incorporate the angular momentum alignment, leading to qubit registers encoding the total angular momentum on each link. The results of a multi-plaquette calculation on IBM's quantum hardware are presented.

preprint2016arXiv

Axial-Current Matrix Elements in Light Nuclei from Lattice QCD

I present results from the first lattice QCD calculations of axial-current matrix elements in light nuclei, performed by the NPLQCD collaboration. Precision calculations of these matrix elements, and the subsequent extraction of multi-nucleon axial-current operators, are essential in refining theoretical predictions of the proton-proton fusion cross section, neutrino-nucleus cross sections and $ββ$-decay rates of nuclei. In addition, they are expected to shed light on the phenomenological quenching of $g_A$ that is required in nuclear many-body calculations.

preprint2016arXiv

Nuclear Physics

Lattice QCD is making good progress toward calculating the structure and properties of light nuclei and the forces between nucleons. These calculations will ultimately refine the nuclear forces, particularly in the three- and four-nucleon sector and the short-distance interactions of nucleons with electroweak currents, and allow for a reduction of uncertainties in nuclear many-body calculations of nuclei and their reactions. After highlighting their importance, particularly to the Nuclear Physics and High-Energy Physics experimental programs, I discuss the progress that has been made toward achieving these goals and the challenges that remain.

preprint2016arXiv

Unitary Limit of Two-Nucleon Interactions in Strong Magnetic Fields

Two-nucleon systems are shown to exhibit large scattering lengths in strong magnetic fields at unphysical quark masses, and the trends toward the physical values indicate that such features may exist in nature. Lattice QCD calculations of the energies of one and two nucleons systems are performed at pion masses of $m_π\sim 450$ and 806 MeV in uniform, time-independent magnetic fields of strength {\bf B}| \sim 10^{19}$-$10^{20}$ Gauss to determine the response of these hadronic systems to large magnetic fields. Fields of this strength may exist inside magnetars and in peripheral relativistic heavy ion collisions, and the unitary behavior at large scattering lengths may have important consequences for these systems.

preprint2015arXiv

Ab initio calculation of the $np \to d γ$ radiative capture process

Lattice QCD calculations of two-nucleon systems are used to isolate the short-distance two-body electromagnetic contributions to the radiative capture process $np \to dγ$, and the photo-disintegration processes $γ^{(\ast)} d \to np$. In nuclear potential models, such contributions are described by phenomenological meson-exchange currents, while in the present work, they are determined directly from the quark and gluon interactions of QCD. Calculations of neutron-proton energy levels in multiple background magnetic fields are performed at two values of the quark masses, corresponding to pion masses of $m_π\sim 450$ and 806 MeV, and are combined with pionless nuclear effective field theory to determine these low-energy inelastic processes. Extrapolating to the physical pion mass, a cross section of $σ^{lqcd}(np\to dγ)=332.4({\tiny \begin{array}{l}+5.4 \\ - 4.7\end{array}})\ mb$ is obtained at an incident neutron speed of $v=2,200\ m/s$, consistent with the experimental value of $σ^{expt}(np \to dγ) = 334.2(0.5)\ mb$.

preprint2015arXiv

Magnetic properties of light nuclei from lattice QCD

After a short review of Lattice QCD methodology and techniques, I summarize recent results of Lattice QCD calculations of the interactions of nucleons and light nuclei with magnetic fields at pion masses of 805 MeV and 450 MeV. Interestingly, the magnetic moments are found to be consistent with the experimental values when given in terms of natural nuclear magnetons. The very low-energy cross section for $np\rightarrow dγ$ is calculated and found to agree with the experimental measurement. First calculations of the magnetic polarizabilities of light nuclei are presented, with a large isovector polarizability observed for the nucleon at these heavier pion masses.

preprint2015arXiv

Nuclear Physics from Lattice Quantum Chromodynamics

Quantum Chromodynamics and Quantum Electrodynamics, both renormalizable quantum field theories with a small number of precisely constrained input parameters, dominate the dynamics of the quarks and gluons - the underlying building blocks of protons, neutrons, and nuclei. While the analytic techniques of quantum field theory have played a key role in understanding the dynamics of matter in high energy processes, they encounter difficulties when applied to low-energy nuclear structure and reactions, and dense systems. Expected increases in computational resources into the exascale during the next decade will provide the ability to determine a range of important strong interaction processes directly from QCD using the numerical technique of Lattice QCD. This will complement the nuclear physics experimental program, and in partnership with new thrusts in nuclear many-body theory, will enable unprecedented understanding and refinement of nuclear forces and, more generally, the visible matter in our universe. In this presentation, I will discuss the state-of-the-art Lattice QCD calculations of quantities of interest in nuclear physics, progress that is expected in the near future, and the anticipated impact.

preprint2015arXiv

The Magnetic Structure of Light Nuclei from Lattice QCD

Lattice QCD with background magnetic fields is used to calculate the magnetic moments and magnetic polarizabilities of the nucleons and of light nuclei with $A\le4$, along with the cross-section for the $M1$ transition $np\rightarrow dγ$, at the flavor SU(3)-symmetric point where the pion mass is $m_π\sim 806$ MeV. These magnetic properties are extracted from nucleon and nuclear energies in six uniform magnetic fields of varying strengths. The magnetic moments are presented in a recent Letter. For the charged states, the extraction of the polarizability requires careful treatment of Landau levels, which enter non-trivially in the method that is employed. The nucleon polarizabilities are found to be of similar magnitude to their physical values, with $β_p=5.22(+0.66/-0.45)(0.23) \times 10^{-4}$ fm$^3$ and $β_n=1.253(+0.056/-0.067)(0.055) \times 10^{-4}$ fm$^3$, exhibiting a significant isovector component. The dineutron is bound at these heavy quark masses and its magnetic polarizability, $β_{nn}=1.872(+0.121/-0.113)(0.082) \times 10^{-4}$ fm$^3$ differs significantly from twice that of the neutron. A linear combination of deuteron scalar and tensor polarizabilities is determined by the energies of the $j_z=\pm 1$ deuteron states, and is found to be $β_{d,\pm 1}=4.4(+1.6/-1.5)(0.2) \times 10^{-4}$ fm$^3$. The magnetic polarizabilities of the three-nucleon and four-nucleon systems are found to be positive and similar in size to those of the proton, $β_{^{3}\rm He}=5.4(+2.2/-2.1)(0.2) \times 10^{-4}$ fm$^3$, $β_{^{3}\rm H}=2.6(1.7)(0.1) \times 10^{-4}$ fm$^3$, $β_{^{4}\rm He}=3.4(+2.0/-1.9)(0.2) \times 10^{-4}$ fm$^3$. Mixing between the $j_z=0$ deuteron state and the spin-singlet $np$ state induced by the background magnetic field is used to extract the short-distance two-nucleon counterterm, ${\bar L}_1$, of the pionless effective theory for $NN$ systems (equivalent to the meson-exchange current contribution in nuclear potential models), that dictates the cross-section for the $np\to dγ$ process near threshold. Combined with previous determinations of NN scattering parameters, this enables an ab initio determination of the threshold cross-section at these unphysical masses.

preprint2014arXiv

Finite-Volume Electromagnetic Corrections to the Masses of Mesons, Baryons and Nuclei

Now that Lattice QCD calculations are beginning to include QED, it is important to better understand how hadronic properties are modified by finite-volume QED effects. They are known to exhibit power-law scaling with volume, in contrast to the exponential behavior of finite-volume strong interaction effects. We use non-relativistic effective field theories describing the low-momentum behavior of hadrons to determine the finite-volume QED corrections to the masses of mesons, baryons and nuclei out to O(1/L^4) in a volume expansion, where L is the spatial extent of the cubic volume. This generalizes the previously determined expansion for mesons, and extends it by two orders in 1/L to include contributions from the charge radius, magnetic moment and polarizabilities of the hadron. We make an observation about direct calculations of the muon g-2 in a finite volume.

preprint2014arXiv

Two-Particle Elastic Scattering in a Finite Volume Including QED

The presence of long-range interactions violates a condition necessary to relate the energy of two particles in a finite volume to their S-matrix elements in the manner of Luscher. While in infinite volume, QED contributions to low-energy charged particle scattering must be resummed to all orders in perturbation theory (the Coulomb ladder diagrams), in a finite volume the momentum operator is gapped, allowing for a perturbative treatment. The leading QED corrections to the two-particle finite-volume energy quantization condition below the inelastic threshold, as well as approximate formulas for energy eigenvalues, are obtained. In particular, we focus on two spinless hadrons in the A1+ irreducible representation of the cubic group, and truncate the strong interactions to the s-wave. These results are necessary for the analysis of Lattice QCD+QED calculations of charged-hadron interactions, and can be straightforwardly generalized to other representations of the cubic group, to hadrons with spin, and to include higher partial waves.

preprint2014arXiv

Uncertainty Quantification in Lattice QCD Calculations for Nuclear Physics

The numerical technique of Lattice QCD holds the promise of connecting the nuclear forces, nuclei, the spectrum and structure of hadrons, and the properties of matter under extreme conditions with the underlying theory of the strong interactions, quantum chromodynamics. A distinguishing, and thus far unique, feature of this formulation is that all of the associated uncertainties, both statistical and systematic can, in principle, be systematically reduced to any desired precision with sufficient computational and human resources. We review the sources of uncertainty inherent in Lattice QCD calculations for nuclear physics, and discuss how each is quantified in current efforts.

preprint2013arXiv

Nuclear Forces from Lattice Quantum Chromodynamics

A century of coherent experimental and theoretical investigations have uncovered the laws of nature that underly nuclear physics. The standard model of strong and electroweak interactions, with its modest number of input parameters, dictates the dynamics of the quarks and gluons - the underlying building blocks of protons, neutrons, and nuclei. While the analytic techniques of quantum field theory have played a key role in understanding the dynamics of matter in high energy processes, they encounter difficulties when applied to low-energy nuclear structure and reactions, and dense systems. Expected increases in computational resources into the exa-scale during the next decade will provide the ability to numerically compute a range of important strong interaction processes directly from QCD with quantifiable uncertainties using the technique of Lattice QCD. These calculations will refine the chiral nuclear forces that are used as input into nuclear many-body calculations, including the three- and four-nucleon interactions. I discuss the state-of-the-art Lattice QCD calculations of quantities of interest in nuclear physics, progress that is expected in the near future, and the impact upon nuclear physics.

preprint2013arXiv

Two-Baryon Systems with Twisted Boundary Conditions

We explore the use of twisted boundary conditions in extracting the nucleon mass and the binding energy of two-baryon systems, such as the deuteron, from Lattice QCD calculations. Averaging the results of calculations performed with periodic and anti-periodic boundary conditions imposed upon the light-quark fields, or other pair-wise averages, improves the volume dependence of the deuteron binding energy from ~exp(-kappa*L)/L to ~exp(-sqrt(2)kappa*L)/L. However, a twist angle of pi/2 in each of the spatial directions improves the volume dependence from ~exp(-kappa*L)/L to ~exp(-2kappa*L)/L. Twist averaging the binding energy with a random sampling of twist angles improves the volume dependence from ~exp^(-kappa*L)/L to ~exp(-2kappa*L)/L, but with a standard deviation of ~exp(-kappa*L)/L, introducing a signal-to-noise issue in modest lattice volumes. Using the experimentally determined phase shifts and mixing angles, we determine the expected energies of the deuteron states over a range of cubic lattice volumes for a selection of twisted boundary conditions.

preprint2013arXiv

Two-Nucleon Systems in a Finite Volume: (II) 3S1-3D1 Coupled Channels and the Deuteron

The energy spectra of two nucleons in a cubic volume provide access to the two phase shifts and one mixing angle that define the S-matrix in the 3S1-3D1 coupled channels containing the deuteron. With the aid of recently derived energy quantization conditions for such systems, and the known scattering parameters, these spectra are predicted for a range of volumes. It is found that extractions of the infinite-volume deuteron binding energy and leading scattering parameters, including the S-D mixing angle at the deuteron pole, are possible from Lattice QCD calculations of two-nucleon systems with boosts of |P| <= 2pi sqrt{3}/L in volumes with 10 fm <~ L <~ 14 fm. The viability of extracting the asymptotic D/S ratio of the deuteron wavefunction from Lattice QCD calculations is discussed.

preprint2012arXiv

Constraints on the Universe as a Numerical Simulation

Observable consequences of the hypothesis that the observed universe is a numerical simulation performed on a cubic space-time lattice or grid are explored. The simulation scenario is first motivated by extrapolating current trends in computational resource requirements for lattice QCD into the future. Using the historical development of lattice gauge theory technology as a guide, we assume that our universe is an early numerical simulation with unimproved Wilson fermion discretization and investigate potentially-observable consequences. Among the observables that are considered are the muon g-2 and the current differences between determinations of alpha, but the most stringent bound on the inverse lattice spacing of the universe, b^(-1) >~ 10^(11) GeV, is derived from the high-energy cut off of the cosmic ray spectrum. The numerical simulation scenario could reveal itself in the distributions of the highest energy cosmic rays exhibiting a degree of rotational symmetry breaking that reflects the structure of the underlying lattice.

preprint2012arXiv

Restoration of Rotational Symmetry in the Continuum Limit of Lattice Field Theories

We explore how rotational invariance is systematically recovered from calculations on hyper-cubic lattices through the use of smeared lattice operators that smoothly evolve into continuum operators with definite angular momentum as the lattice-spacing is reduced. Perturbative calculations of the angular momentum violation associated with such operators at tree-level and at one-loop are presented in phi^4 theory and QCD. Contributions from these operators that violate rotational invariance occur at tree-level, with coefficients that are suppressed by O(a^2) in the continuum limit. Quantum loops do not modify this behavior in phi^4, nor in QCD if the gauge-fields are smeared over a comparable spatial region. Consequently, the use of this type of operator should, in principle, allow for Lattice QCD calculations of the higher moments of the hadron structure functions.

preprint2011arXiv

Extracting Scattering Phase-Shifts in Higher Partial-Waves from Lattice QCD Calculations

Lüscher's method is routinely used to determine meson-meson, meson-baryon and baryon-baryon s-wave scattering amplitudes below inelastic thresholds from Lattice QCD calculations - presently at unphysical light-quark masses. In this work we review the formalism and develop the requisite expressions to extract phase-shifts describing meson-meson scattering in partial-waves with angular-momentum l<=6 and l=9. The implications of the underlying cubic symmetry, and strategies for extracting the phase-shifts from Lattice QCD calculations, are presented, along with a discussion of the signal-to-noise problem that afflicts the higher partial-waves.

preprint2011arXiv

Improving the Volume Dependence of Two-Body Binding Energies Calculated with Lattice QCD

Volume modifications to the binding of two-body systems in large cubic volumes of extent L depend upon the total momentum and exponentially upon the ratio of L to the size of the boosted system. Recent work by Bour et al determined the momentum dependence of the leading volume modifications to nonrelativistic systems with periodic boundary conditions imposed on the single-particle wavefunctions, enabling them to numerically determine the scattering of such bound states using a low-energy effective field theory and Luscher's finite-volume method. The calculation of bound nuclear systems directly from QCD using Lattice QCD has begun, and it is important to reduce the systematic uncertainty introduced into such calculations by the finite spatial extent of the gauge-field configurations. We extend the work of Bour et al from nonrelativistic quantum mechanics to quantum field theory by generalizing the work of Luscher and of Gottlieb and Rummukainen to boosted two-body bound states. The volume modifications to binding energies can be exponentially reduced from ~ e^{-kappa L}/L to ~ e^{-2 kappa L}/L in nonrelativistic systems (where kappa is the binding momentum of the state) by forming particular combinations of the binding energies determined in the four lowest-lying boosted systems. Relativistic corrections to this combination, and others, that violate the exponential reduction are determined. An analysis of what can be expected from Lattice QCD calculations of the deuteron is performed, the results of which are representative of a generic loosely bound system.

preprint2010arXiv

A method to study complex systems of mesons in Lattice QCD

Finite density systems can be explored with Lattice QCD through the calculation of multi-hadron correlation functions. Recently, systems with up to 12 $π^+$'s or $K^+$'s have been studied to determine the 3-$π^+$ and 3-$K^+$ interactions, and the corresponding chemical potentials have been determined as a function of density. We derive recursion relations between correlation functions that allow this work to be extended to systems of arbitrary numbers of mesons and to systems containing many different types of mesons, such as $π^+$'s, $K^+$'s, $\bar{D}^0$'s and $B^+$'s. These relations allow for the study of finite-density systems in arbitrary volumes, and for the study of high-density systems.

preprint2009arXiv

High Statistics Analysis using Anisotropic Clover Lattices: (I) Single Hadron Correlation Functions

We present the results of high-statistics calculations of correlation functions generated with single-baryon interpolating operators on an ensemble of dynamical anisotropic gauge-field configurations generated by the Hadron Spectrum Collaboration using a tadpole-improved clover fermion action and Symanzik-improved gauge action. A total of 292,500 sets of measurements are made using 1194 gauge configurations of size 20^3 x 128 with an anisotropy parameter ξ= b_s/b_t = 3.5, a spatial lattice spacing of b_s=0.1227\pm 0.0008 fm, and pion mass of m_π~ 390 MeV. Ground state baryon masses are extracted with fully quantified uncertainties that are at or below the ~0.2%-level in lattice units. The lowest-lying negative-parity states are also extracted albeit with a somewhat lower level of precision. In the case of the nucleon, this negative-parity state is above the Nπthreshold and, therefore, the isospin-1/2 πN s-wave scattering phase-shift can be extracted using Luescher's method. The disconnected contributions to this process are included indirectly in the gauge-field configurations and do not require additional calculations. The signal-to-noise ratio in the various correlation functions is explored and is found to degrade exponentially faster than naive expectations on many time-slices. This is due to backward propagating states arising from the anti-periodic boundary conditions imposed on the quark-propagators in the time-direction. We explore how best to distribute computational resources between configuration generation and propagator measurements in order to optimize the extraction of single baryon observables.

preprint2009arXiv

High Statistics Analysis using Anisotropic Clover Lattices: (III) Baryon-Baryon Interactions

Low-energy baryon-baryon interactions are calculated in a high-statistics lattice QCD study on a single ensemble of anisotropic clover gauge-field configurations at a pion mass of m_π~ 390 MeV, a spatial volume of L^3 ~ (2.5 fm)^3, and a spatial lattice spacing of b~0.123 fm. Lüscher's method is used to extract nucleon-nucleon, hyperon-nucleon and hyperon-hyperon scattering phase shifts at one momentum from the one- and two-baryon ground-state energies in the lattice volume. The isospin-3/2 NΣinteractions are found to be highly spin-dependent, and the interaction in the ^3S_1 channel is found to be strong. In contrast, the NΛinteractions are found to be spin-independent, within the uncertainties of the calculation, consistent with the absence of one-pion-exchange. The only channel for which a negative energy-shift is found is ΛΛ, indicating that the ΛΛinteraction is attractive, as anticipated from model-dependent discussions regarding the H-dibaryon. The NN scattering lengths are found to be small, clearly indicating the absence of any fine-tuning in the NN-sector at this pion mass. This is consistent with our previous Lattice QCD calculation of NN interactions. The behavior of the signal-to-noise ratio in the baryon-baryon correlation functions, and in the ratio of correlation functions that yields the ground-state energy splitting is explored. In particular, focus is placed on the window of time slices for which the signal-to-noise ratio does not degrade exponentially, as this provides the opportunity to extract quantitative information about multi-baryon systems.

preprint2009arXiv

Meson-Baryon Scattering Lengths from Mixed-Action Lattice QCD

The $π^+Σ^+$, $π^+Ξ^0$, $K^+p$, $K^+n$, and $\overline{K}{}^0 Ξ^0$ scattering lengths are calculated in mixed-action Lattice QCD with domain-wall valence quarks on the asqtad-improved coarse MILC configurations at four light-quark masses, and at two light-quark masses on the fine MILC configurations. Heavy Baryon Chiral Perturbation Theory with two and three flavors of light quarks is used to perform the chiral extrapolations. To the order we work in the three-flavor chiral expansion, the kaon-baryon processes that we investigate show no signs of convergence. Using the two-flavor chiral expansion for extrapolation, the pion-hyperon scattering lengths are found to be $a_{π^+Σ^+}=-0.197\pm0.017$ fm, and $a_{π^+Ξ^0}=-0.098\pm0.017$ fm, where the comprehensive error includes statistical and systematic uncertainties.

preprint2007arXiv

Multi-Pion Systems in Lattice QCD and the Three-Pion Interaction

The ground-state energies of 2, 3, 4 and 5 π^+'s in a spatial volume V (2.5 fm)^3 are computed with lattice QCD. By eliminating the leading contribution from three-π^+ interactions, particular combinations of these n-π^+ ground-state energies provide precise extractions of the π^+π^+ scattering length in agreement with that obtained from calculations involving only two π^+'s. The three-π^+ interaction can be isolated by forming other combinations of the n-π^+ ground-state energies. We find a result that is consistent with a repulsive three-π^+ interaction for m_π< 352 MeV.

preprint2005arXiv

Extrapolation Formulas for Neutron EDM Calculations in Lattice QCD

Lattice QCD is rapidly progressing toward being able to reliably compute the electric dipole moment of the neutron as a function of the strong CP-violating parameter theta. Present day calculations are performed at unphysical values of the light quark masses, in volumes that are not exceptionally large and at lattice spacings that are not exceptionally small. We use chiral perturbation theory to determine the leading contributions to the neutron electric dipole moment at finite volume, and in partially-quenched calculations.

preprint2004arXiv

DVCS-Dissociation of the Deuteron and the EMC Effect

The break-up of the deuteron during deeply-virtual Compton scattering, gamma* d --> gamma(*) n p, is explored. In the effective field theory describing nucleon dynamics at momenta below the pion mass, the EMC effect results from four-nucleon interactions with the twist-2 operators, appropriate for describing forward, and near-forward, matrix elements in the two-nucleon system. We point out that the break-up of the deuteron to low-energy final states during deeply-virtual Compton scattering is a process with which to explore strong-interaction physics closely related to that responsible for the EMC effect. The single-nucleon contribution to the break-up depends on the moments of the spin-dependent structure functions and contributions from local four-nucleon operators. Experimental deviations from the single-nucleon prediction would provide a probe of strong interactions complimentary to the EMC effect.

preprint2003arXiv

Goldstone Bosons in the 3P2 Superfluid Phase of Neutron Matter and Neutrino Emission

At the high densities present in the interior of neutron stars, the neutrons are condensed into the 3P2 superfluid phase. While this condensation has little impact on the equation of state, it can have an important role in determining the low-temperature energy-momentum transport properties. The spontaneous breaking of baryon number by the condensate gives rise to the familiar Goldstone boson, but in addition, the spontaneous breaking of rotational invariance by the condensate gives rise to three Goldstone bosons, in general, one for each broken generator of rotations. These Goldstone bosons, which couple to the Z0, provide a new mechanism for neutrino emission. Using a low-energy effective field theory to describe the dynamics of these Goldstone bosons we estimate the neutrino emissivity of dense neutron matter and show that their annihilation is the dominant energy-loss mechanism over a range of temperatures.

preprint2002arXiv

Hadronic Parity-Violation on the Lattice

We motivate lattice QCD studies of the parity-violating pion-nucleon coupling constant and extend flavor-conserving hadronic parity-violation from QCD to partially-quenched QCD. The parity-violating pion-nucleon coupling and the anapole form factor (and moment) of the proton are computed to one-loop order in the partially-quenched chiral expansion. For the parity-violating pion-nucleon interaction, we include the contributions from total derivative operators necessary to match the kinematics that will be used in lattice simulations.

preprint2002arXiv

Variation of Fundamental Couplings and Nuclear Forces

The dependence of the nuclear force on standard model parameters plays an important role in bounding time and space variations of fundamental couplings over cosmological time scales. We discuss the quark-mass dependence of deuteron and di-neutron binding in a systematic chiral expansion. The leading quark-mass dependence of the nuclear force arises from one-pion exchange and from local quark-mass dependent four-nucleon operators with coefficients that are presently unknown. By varying these coefficients while leaving nuclear observables at the physical values of the quark masses invariant, we find scenarios where two-nucleon physics depends both weakly and strongly on the quark masses. While the determination of these coefficients is an exciting future opportunity for lattice QCD, we conclude that, at present, bounds on time and space variations of fundamental parameters from the two-nucleon sector are much weaker than previously claimed. This brings into question the reliability of coupling-constant bounds derived from more complex nuclei and nuclear processes.

preprint2001arXiv

Parity Violation in Low-Energy np->d gamma and the Deuteron Anapole Moment

Parity violation in low-energy nuclear observables is included in the pionless effective field theory. The model-independent relation between the parity-violating asymmetry in polarized np -> d gamma and the non-nucleon part of the deuteron anapole moment is discussed. The asymmetry in np -> d gamma computed with KSW power-counting, and recently criticized by Desplanques, is discussed.

preprint2000arXiv

Extra dimensions, SN1987a, and nucleon-nucleon scattering data

One of the strongest constraints on the existence of large, compact, "gravity-only" dimensions comes from SN1987a. If the rate of energy loss into these putative extra dimensions is too high, then the neutrino pulse from the supernova will differ from that actually seen. The dominant mechanism for the production of Kaluza-Klein gravitons and dilatons in the supernova is via gravistrahlung and dilastrahlung from the nucleon-nucleon system. In this paper we compute the rates for these processes in a model-independent way using low-energy theorems which relate the emissivities to the measured nucleon-nucleon cross section. This is possible because for soft gravitons and dilatons the leading contribution to the energy-loss rate is from graphs in which the gravitational radiation is produced from external nucleon legs. Previous calculations neglected these mechanisms. We re-evaluate the bounds on toroidally-compactified "gravity-only" dimensions (GODs), and find that consistency with the observed SN1987a neutrino signal requires that if there are two such dimensions then their radius must be less than 1 micron.

preprint2000arXiv

From Hadrons to Nuclei: Crossing the Border

The study of nuclei predates by many years the theory of quantum chromodynamics. More recently, effective field theories have been used in nuclear physics to ``cross the border'' from QCD to a nuclear theory. We are now entering the second decade of efforts to develop a perturbative theory of nuclear interactions using effective field theory. This work describes the current status of these efforts.

preprint2000arXiv

Rearranging Pionless Effective Field Theory

We point out a redundancy in the operator structure of the pionless effective field theory which dramatically simplifies computations. This redundancy is best exploited by using dibaryon fields as fundamental degrees of freedom. In turn, this suggests a new power counting scheme which sums range corrections to all orders. We explore this method with a few simple observables: the deuteron charge form factor, n p -> d gamma, and Compton scattering from the deuteron. Higher dimension operators involving electroweak gauge fields are not renormalized by the s-wave strong interactions, and therefore do not scale with inverse powers of the renormalization scale. Thus, naive dimensional analysis of these operators is sufficient to estimate their contribution to a given process.

preprint1999arXiv

Improving the Convergence of NN Effective Field Theory

We study a low-energy effective field theory (EFT) describing the NN system in which all exchanged particles are integrated out. We show that fitting the residue of the 3S1 amplitude at the deuteron pole, rather than the 3S1 effective range, dramatically improves the convergence of deuteron observables in this theory. Reproducing the residue ensures that the tail of the deuteron wave function, which is directly related to NN scattering data via analytic continuation, is correctly reproduced in the EFT at next-to-leading order. The role of multi-nucleon-electroweak operators which produce deviations from effective-range theory can then be explicitly separated from the physics of the wave function tail. Such an operator contributes to the deuteron quadrupole moment, mu_Q, at low order, indicating a sensitivity to short-distance physics. This is consistent with the failure of impulse approximation calculations in NN potential models to reproduce mu_Q. The convergence of NN phase shifts in the EFT is unimpaired by the use of this new expansion.

preprint1997arXiv

Non-1/m_b^n Power Suppressed Contributions to Inclusive b->s l+ l- Decays

We compute non-perturbative contributions to b->s l+ l- that are not explicitly suppressed by powers of the b-quark mass. They are proportional to lambda_2 and arise from an interference between the free-quark amplitude and higher order terms in the matrix element of a four-quark operator. This correction is found to be small over most of the dalitz plot except near the charm threshold. Unfortunately, the perturbative computation we have performed is invalid near charm threshold and we do not except to see the structure found at lowest order reproduced in the data. We conclude that these non-perturbative contributions do not significantly modify the previous analysis of b->s l+ l-.

preprint1996arXiv

Heavy Baryon Mixing in Chiral Perturbation Theory

We discuss the SU(3) and heavy quark spin-symmetry breaking mixing between the Xi_c and Xi'_c charmed baryons. Chromomagnetic hyperfine interactions are the leading source of spin-symmetry breaking and together with the SU(3) breaking mass differences between the lightest pseudo-Goldstone bosons gives the leading contribution to the mixing. Such contributions are computed in chiral perturbation theory and compared to quark model expectations. We also compute the leading contribution to the semileptonic decay Xi_b -> Xi'_c l nu at zero recoil, and find that it is an order of magnitude smaller than naive power counting would suggest. It appears that Xi_b -> Xi'_c l nu is dominated by incalculable counterterms, and we discuss the implications for quark models based on the essential role of hyperfine interactions.

preprint1996arXiv

Higher Order QCD Corrections to $b\to c\overline{c}s$

We calculate the ${\cal O}(α_s^2 β_0 )$ corrections to the decay rate $b\rightarrow c\bar c s$. For reasonable values of $m_c/m_b$ this term is of the same order as both the one-loop and $\CO(α_s^2 \log^2 m_W/m_b)$ corrections to the decay rate. For $m_c/m_b=0.3$ the $\CO(α_s^2β_0)$ corrections enhance the rate by $\sim 18\%$ . We also discuss the $\CO(α_s^2β_0)$ corrections to $R_τ$, the $B$ semileptonic branching fraction and the charm multiplicity.

preprint1996arXiv

Phenomenology of the $1/m_Q$ Expansion in Inclusive $B$ and $D$ Meson Decays

We apply a recent theoretical analysis of hadronic observables in inclusive semileptonic heavy hadron decays to the phenomenology of $B$ and $D$ mesons. Correlated bounds on the nonperturbative parameters $\barΛ$ and $λ_1$ are derived by considering data from $B$ decays and, independently, data from $D$ decays. The two sets of bounds are found to be consistent with each other. The data from $B$ decays are then are used to extract a lower limit on the CKM matrix element $|V_{cb}|$. We address the issue of the convergence of the perturbative expansions used in the analysis, and compare our bounds on $\barΛ$ and $λ_1$ to lattice and QCD sum rule results. Finally, we argue that a comparison of the analyses of $D$ and $D_s$ decays provides evidence for the applicability of parton-hadron duality in the semileptonic decay of charmed hadrons.

preprint1995arXiv

Hadron Spectra for Semileptonic Heavy Quark Decay

We calculate the leading perturbative and power corrections to the hadronic invariant mass and energy spectra in semileptonic heavy hadron decays. We apply our results to the $B$ system. Moments of the invariant mass spectrum, which vanish in the parton model, probe gluon bremsstrahlung and nonperturbative effects. Combining our results with recent data on $B$ meson branching ratios, we obtain a lower bound $\barΛ>410\,{\rm MeV}$ and an upper bound $m_b^{\rm pole}<4.89\,$GeV. The Brodsky-Lepage-Mackenzie scale setting procedure suggests that higher order perturbative corrections are small for bottom decay, and even tractable for charm decay.

preprint1995arXiv

K_L to π^+ π^- e^+ e^-

We calculate all of the form factors for the one-photon, $K_L \rightarrow π^+ π^- γ^* \rightarrow π^+ π^- e^+e^-$ contribution to the $K_L \rightarrow π^+ π^- e^+e^-$ decay amplitude at leading order in chiral perturbation theory. These form factors depend on one unknown constant that is a linear combination of coefficients of local ${\cal O}(p^4)$ operators in the chiral lagrangian for weak radiative kaon decay. We determine the differential rate for $\K$ and also the magnitude of two CP violating observables.

preprint1995arXiv

The Radiative Charmed Baryon Decay $Ξ_{c2}^{0*} \to Ξ_{c1}^0 γ$

V-spin symmetry ($s \leftrightarrow d$ symmetry) forbids the radiative decay $Ξ_{c2}^{0*} \rightarrow Ξ_{c1}^0 γ$ in the SU(3) limit. The quark mass term breaks V-spin symmetry and the leading nonanalytic contribution to the radiative decay amplitude is computable in heavy baryon chiral perturbation theory. The radiative decay branching ratio is determined by the coupling constant $g_2$ and at leading order in chiral perturbation theory is given by $Br(Ξ_{c2}^{0*} \rightarrow Ξ_{c1}^0 γ) = 1.0\times 10^{-3} g_2^2$. Measurement of this branching fraction will determine $|g_2|$.

preprint1994arXiv

Charm Mass Dependence of the $\CO(α_s^2 n_f)$ Correction to Inclusive $B\rightarrow X_c e\barν_e$ Decay

We compute the $α_s^2 n_f$ perturbative QCD contribution to semileptonic B decay, including the finite mass of the charm quark. This result provides an estimate of the size of the two-loop correction, which is found to be about 50\% of the one-loop correction. We use these results to set the scale for the one-loop correction using the scheme of Brodsky, Lepage and Mackenzie and find a BLM scale of $μ_{\rm BLM} = 0.13\, m_b$, when the inclusive semileptonic rate is expressed in terms of the $b$ and $c$ quark pole masses and the \msbar strong coupling. The two loop correction lies roughly midway between that obtained at $m_c=0$ and that obtained in the Shifman-Voloshin limit $m_b, m_c\gg m_b-m_c\ggΛ_{\rm QCD}$ while the corresponding BLM scale is somewhat closer to that obtained in the former case.

preprint1994arXiv

E2 Strength in the Radiative Charmed Baryon Decay

The radiative decay $Σ_Q^* \rightarrow Λ_Qγ$ can have both magnetic dipole (M1) and electric quadrupole (E2) components. In the heavy quark limit $M_Q\rightarrow\infty$ the transition arises from the spin of the light degrees of freedom changing from $s_l=1$ to $s_l=0$ and hence the E2 contribution vanishes. We compute the leading contribution to the E2 strength in chiral perturbation theory and find that the amplitude is enhanced by a small energy denominator in the chiral limit. This enhancement essentially compensates for the $1/M_c$ suppression that is present in the charm system. We find a mixing ratio of order a few percent dependent upon the $Σ_c^*--Σ_c$ spin symmetry breaking mass difference. The analogous quantity in the b-baryon sector is smaller by a factor of $\sim M_c/M_b$.

preprint1994arXiv

Perturbative Strong Interaction Corrections to the Heavy Quark Semileptonic Decay Rate

We calculate the part of the order $α_s^2$ correction to the semileptonic heavy quark decay rate proportional to the number of light quark flavors, and use our result to set the scale for evaluating the strong coupling in the order $α_s$ term according to the scheme of Brodsky, Lepage and Mackenzie. Expressing the decay rate in terms of the heavy quark pole mass $m_Q$, we find the scale for the $\overline{MS}$ strong coupling to be $0.07\, m_Q$. If the decay rate is expressed in terms of the $\overline{MS}$ heavy quark mass $\overline m_Q(m_Q)$ then the scale is $0.12\, m_Q$. We use these results along with the existing calculations for hadronic $τ$ decay to calculate the BLM scale for the nonleptonic decay width and the semileptonic branching ratio. The implications for the value of $|V_{bc}|$ extracted from the inclusive semileptonic $B$ meson decay rate are discussed.

preprint1994arXiv

Strong $Λπ$ Phase Shifts for CP Violation in Weak $Ξ\rightarrow Λπ$ Decay

Strong interaction $Λπ$ phase shifts relevant for the weak nonleptonic decay $Ξ\rightarrow Λπ$ are calculated using baryon chiral perturbation theory. We find in leading order that the S-wave phase shift vanishes and the $J={1 \over 2}$ P-wave phase shift is $-1.7 ^{\rm o} $. The small phase shifts imply that CP violation in this decay will be difficult to observe. Our results follow from chiral $SU(2)_L\times SU(2)_R$ symmetry.

preprint1993arXiv

Extracting $|V_{bc}|$, $m_c$ and $m_b$ from Inclusive $D$ and $B$ Decays

Using recent results for nonperturbative contributions to the $B$ and $D$ meson inclusive semileptonic widths, a model independent extraction of $\vbc$, $m_c$ and $m_b$ is made from the experimentally measured $B$ and $D$ lifetimes and semileptonic branching ratios. Constraining the parameters of the HQET at $\CO(1/m_Q^2)$ by the $D$ semileptonic width, $\vbc$ is found to lie in the range $.040<\vbc< 0.057$. The $c$ and $b$ quark masses are not well constrained due to uncertainty in the relevant scale of $α_s$. These results assume the validity of perturbative QCD at the low scales relevant to semileptonic charm decay. Without making this assumption, somewhat less stringent bounds on $V_{bc}$ from $B$ decay alone may be obtained.

preprint1993arXiv

Heavy Quark Fragmentation to Baryons Containing Two Heavy Quarks

We discuss the fragmentation of a heavy quark to a baryon containing two heavy quarks of mass $m_Q\ggΛ_{\rm QCD}$. In this limit the heavy quarks first combine perturbatively into a compact diquark with a radius small compared to $1/Λ_{\rm QCD}$, which interacts with the light hadronic degrees of freedom exactly as does a heavy antiquark. The subsequent evolution of this $QQ$ diquark to a $QQq$ baryon is identical to the fragmentation of a heavy antiquark to a meson. We apply this analysis to the production of baryons of the form $ccq$, $bbq$, and $bcq$.

preprint1992arXiv

A QCD Calculation of the Interaction of Quarkonium with Nuclei

The interaction of quarkonium with nuclei is studied in the $m_Q\rightarrow \infty$ limit of QCD, where the binding energy is found to be exactly computable. The dominant contribution to the interaction is from two-gluon operators. The forward matrix elements of these two-gluon operators can be determined from the QCD scale anomaly, and from deep inelastic scattering. We apply our results to the $Υ$ and $J/ψ$, treating the $\qqbar$ interaction as purely Coulombic. We find the $Υ$ binds in nuclear matter with a binding energy of a few $\mev$, while for the $J/ψ$ binding is of order 10 $\mev$. For the $J/ψ$ in particular we expect confinement effects to produce large corrections to this result.

preprint1992arXiv

Semileptonic $B_c$ Decay and Heavy Quark Spin Symmetry

Semileptonic decay of the $B_c$ meson is studied in the heavy quark limit. The six possible form factors for $B_c \rightarrow B_s (B^0),B_s^* (B^{*0})$ semileptonic decay are determined by two invariant functions. Only one of these functions contributes at zero recoil, where it is calculable to lowest order in an operator product expansion in terms of the meson decay constant $f_B$ and the $B_c$ wavefunction. A similar result is found for $B_c \rightarrow D^0,D^{*0}$ and for $B_c\rightarrowη_c,J/ψ$ semileptonic decay for a restricted kinematic region. Semileptonic $B_c$ decay provides a means for determining the KM mixing angle $|V_{ub}|$.

preprint1992arXiv

Two Photon Contribution to Polarization in $K^+ \rightarrow π^+ μ^+ μ^-$

Short distance physics involving virtual top and charm quarks contributes to $μ^+$ (and $μ^-$) polarization in the decay $K^+ \rightarrow π^+ μ^+ μ^-$. Measurement of the parity violating asymmetry $(Γ_R - Γ_L)/(Γ_R + Γ_L)$, where $Γ_R$ and $Γ_L$ are the rates to produce right and left-handed $μ^+$, may provide valuable information on the unitarity triangle. The parity violating asymmetry also gets a contribution from Feynman diagrams with two photon intermediate states. We estimate this two photon contribution to the asymmetry and discuss briefly the two photon contribution to time reversal odd asymmetries that involve both the $μ^+$ and $μ^-$ polarizations.