BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Tanmoy Bhattacharya

Tanmoy Bhattacharya appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
hep-lathep-phnucl-thquant-phhep-thcond-mat.str-elMachine Learningnucl-exComputation and LanguageCryptography and Securityhep-exphysics.soc-phPopulations and EvolutionQuantitative Methods

Catalog footprint

What is connected

59works
14topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

59 published item(s)

preprint2026arXiv

The Spectrum and Scale Setting on 2+1-flavor NME Lattices

This paper describes the thirteen ensembles, named NME, generated with 2+1-flavor Wilson-clover fermions by the JLab/W\&M/LANL/MIT/Marseille collaborations, and presents an analysis of the meson and baryon spectrum, decay constants $f_π$ and $f_K$, flow scales $t_0$ and $w_0$, and time histories of the $Θ$ and Weinberg operators under gradient flow. Using these quantities, the physical point values of the two flow scales, ${t_0^{\rm Phy}}$ and ${w_0^{\rm Phy}}$, and the ratio $\mathop{f_K / f_π}^{\rm Phy}$ are determined. The masses of the octet and decuplet baryons are analyzed using both the next-to-leading order (NLO) and the next-next-to-leading order (NNLO) ansatz from heavy baryon chiral perturbation theory (HB$χ$PT). The NNLO fit to the octet baryons, $M_N$, $M_Σ$, $M_Λ$ and $M_Ξ$, is preferred while the corresponding fits to the decuplet Omega mass, $M_Ω$, are not distinguished. We also present a study of the autocorrelations in the data and show that there is no evidence, even at large flow time, of the freezing of the topological charge or the Weinberg three-gluon operator.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Update on flavor diagonal nucleon charges from clover fermions

We present a summary of the full calculation of the axial, scalar and tensor flavor diagonal charges of the nucleon carried out using Wilson-clover fermions on eight ensembles generated using 2+1+1-flavors of highly improved staggered quarks (HISQ) by the MILC collaboration. We also give results for the $3\times 3$ matrix of renormalization factors between the RI-sMOM and $\overline{\rm MS}$ scheme for the 2+1 flavor theory that include flavor mixing. Preliminary results for $g_{A,S,T}^{u,d,s}$ are presented in the $\overline{\rm MS}$ scheme at scale 2 GeV.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

nEDM from the theta-term and chromoEDM operators

In a previous work, we showed that unresolved excited state contaminations provide a major source of systematic uncertainty in the calculation of the nucleon electric dipole moment due to the QCD topological term theta. Here we extend the calculation to the quark chromo-electric dipole moment operator (qcEDM). We also show quantitatively the impact of the mixing of the qcEDM with lower-dimensional operators on the lattice. Finally, we present preliminary results from a unitary clover-on-clover calculation for the QCD topological term.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Nucleon form factors and the pion-nucleon sigma term

This talk summarizes the progress made since Lattice 2021 in understanding and controlling the contributions of towers of multihadron excited states with mass gaps starting lower than of radial excitations, and in increasing our confidence in the extraction of ground state nucleon matrix elements. The most clear evidence for multihadron excited state contributions (ESC) is in axial/pseudoscalar form factors that are required to satisfy the PCAC relation between them. The talk examines the broader question--which and how many of the theoretically allowed positive parity states $N(\textbf p)π(-\textbf p)$, $N(\textbf 0)π(\textbf 0)π(\textbf 0)$, $N(\textbf p)π(\textbf 0)$, $N(\textbf 0)π(\textbf p),\ \ldots$ make significant contributions to a given nucleon matrix element? New data for the axial, electric and magnetic form factors are presented. They continue to show trends observed in Ref[1]. The N${}^2$LO $χ$PT analysis of the ESC to the pion-nucleon sigma term, $σ_{πN}$, has been extended to include the $Δ$ as an explicit degree of freedom [2]. The conclusion reached in Ref [3] that $N π$ and $N ππ$ states each contribute about 10 MeV to $σ_{πN}$, and the consistency between the lattice result with $N π$ state included and the phenomenological estimate is not changed by this improvement.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Phases of 2d massless QCD with qubit regularization

We investigate the possibility of reproducing the continuum physics of 2d SU(N) gauge theory coupled to a single flavor of massless Dirac fermions using qubit regularization. The continuum theory is described by N free fermions in the ultraviolet (UV) and a coset Wess-Zumino-Witten (WZW) model in the infrared (IR). In this work, we explore how well these features can be reproduced using the Kogut-Susskind Hamiltonian with a finite-dimensional link Hilbert space and a generalized Hubbard coupling. Using strong coupling expansions, we show that our model exhibits a gapped dimer phase and another phase described by a spin-chain. Furthermore, for N=2, using tensor network methods, we show that there is a second-order phase transition between these two phases. The critical theory at the transition can be understood as an SU(2)_1 WZW model, using which we determine the phase diagram of our model quantitatively. Using the confinement properties of the model we argue how the UV physics of free fermions could also emerge, but may require further modifications to our model.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Synergy between NP and HEP research goals and efforts in fundamental symmetries and interactions

The aim of this white paper is to highlight several areas for which the Department of Energy's Office of Nuclear Physics has primary stewardship or significant investment and expertise, and for which there is also significant interest and expertise within the HEP community. These areas of overlap offer exciting opportunities for collaboration.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Update on Flavor Diagonal Nucleon Charges

This talk provides an update on the calculation of matrix elements of flavor diagonal axial, scalar and tensor quark bilinear operators between the nucleon ground state. The simulations are done using Wilson-clover fermions on a sea of eight 2+1+1-flavor HISQ ensembles generated by the MILC collaboration. We discuss the signal in the sum of the connected and disconnected contributions for the up, down and strange quarks, control over fits to remove excited state contamination, and the simultaneous chiral-continuum fit used to extract the charges.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

AI and Theoretical Particle Physics

Theoretical particle physicists continue to push the envelope in both high performance computing and in managing and analyzing large data sets. For example, the goals of sub-percent accuracy in predictions of quantum chromodynamics (QCD) using large scale simulations of lattice QCD and in finding signals of rare events and new physics in exabytes of data produced by experiments at the high luminosity large hadron collider (LHC) require new tools beyond just developments in hardware. Machine learning and artificial intelligence offer the promise of dramatically reducing the computational cost and time. This chapter reviews selected areas where AI/ML tools could have a major impact, provides an overview of the challenges, and discusses how new ideas such as normalizing flows can speed up the generation of gauge configurations needed in lattice QCD calculations; the growth of ML in surrogate models and pattern matching to reduce the cost of event generators and in the analysis of experimental data; and in the search for viable vacua in the landscape of string theories. While such approaches transform aspects of particle theory into computational problems, and thus black boxes, we argue that physics-aware development of these tools combined with algorithms that ensure that the results are bias free will continue to require a deep understanding of the physics. We see this broader transformation as akin to formulating and extracting observables from simulations of lattice QCD, a numerical integration of the path integral formulation of QCD that nevertheless requires a deep understanding of the underlying quantum field theory, the standard model of particle physics and effective field theory methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Calculation of neutron electric dipole moment due to the QCD topological term, Weinberg three-gluon operator and the quark chromoelectric moment

We present results for the neutron electric dipole moment due to the dimension 4 and dimension 6 gluonic CP violation, and the isovector quark chromoelectric dipole moment using clover valence quarks on HISQ dynamical ensembles generated by the MILC Collaboration. For the gluonic operators, we use the gradient flow scheme to obtain divergence-free continuum results. For the chromoelectric dipole moment operator, we use the unflowed local operator but discuss how the quadratically divergent mixing with the pseudoscalar operator can be controlled nonperturbatively.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Excited states and precision results for nucleon charges and form factors

The exponentially falling signal-to-noise ratio in all nucleon correlation functions, and the presence of towers of multihadron excited states with relatively small mass gaps makes extraction of matrix elements of various operators within the ground state nucleon challenging. Theoretically, the allowed positive parity states with the smallest mass gaps are the $N(\bm p)π(-\bm p)$, $N(\bm 0)π(\bm 0)π(\bm 0)$, $N(\bm p)π(\bm 0)$, $N(\bm 0)π(\bm p),\ \ldots$, states. A priori, the contribution of these states arises at one loop in chiral perturbation theory ($χ$PT), however, in many cases the contributions are enhanced. In this talk, I will review four such cases: the correlation functions from which the axial form factors, electric and magnetic form factors, the $Θ$-term contribution to neutron electric dipole moment (nEDM), and the pion-nucleon sigma term are extracted. Including appropriate multihadron states in the analysis can lead to significantly different results compared to standard analyses with the mass gaps taken from fits to 2-point functions. The $χ$PT case for $N π$ states is the most clear in the axial/pseudoscalar form factors which need to satisfy the PCAC relation between them. Our analyses, supported by $χ$PT, suggests similarly large effects in the calculations of the $Θ$-term and the pion-nucleon sigma term that have significant phenomenological implications.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Flavor Diagonal Nucleon Charges

This talk provides an update on the calculation of matrix elements of flavor diagonal axial, scalar and tensor quark bilinear operators between the nucleon ground state. The simulations are done using Wilson-clover fermions on a sea of eight 2+1+1-flavor HISQ ensembles generated by the MILC collaboration. We discuss the signal in the connected and disconnected contributions, calculation of the renormalization constants and mixing in the RI-sMOM scheme, and control over the simultaneous chiral-continuum-finite-volume fit used to extract the final charges.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Precision Nucleon Charges and Form Factors Using 2+1-flavor Lattice QCD

We present high statistics results for the isovector nucleon charges and form factors using seven ensembles of 2+1-flavor Wilson-clover fermions. The axial and pseudoscalar form factors obtained on each ensemble satisfy the PCAC relation once the lowest energy $Nπ$ excited state is included in the spectral decomposition of the correlation functions used for extracting the ground state matrix elements. Similarly, we find evidence that the $Nππ$ excited state contributes to the correlation functions with the vector current, consistent with the vector meson dominance model. The resulting form factors are consistent with the Kelly parameterization of the experimental electric and magnetic data. Our final estimates for the isovector charges are $g_{A}^{u-d} = 1.31(06)(05)_{sys}$, $g_{S}^{u-d} = 1.06(10)(06)_{sys}$, and $g_{T}^{u-d} = 0.95(05)(02)_{sys}$, where the first error is the overall analysis uncertainty and the second is an additional combined systematic uncertainty. The form factors yield: (i) the axial charge radius squared, ${\langle r_A^2 \rangle}^{u-d}=0.428(53)(30)_{sys}\ {\rm fm}^2$, (ii) the induced pseudoscalar charge, $g_P^\ast=7.9(7)(9)_{sys}$, (iii) the pion-nucleon coupling $g_{π{\rm NN}} = 12.4(1.2)$, (iv) the electric charge radius squared, ${\langle r_E^2 \rangle}^{u-d} = 0.85(12)(19)_{sys} \ {\rm fm}^2$, (v) the magnetic charge radius squared, ${\langle r_M^2 \rangle}^{u-d} = 0.71(19)(23)_{\rm sys} \ {\rm fm}^2$, and (vi) the magnetic moment $μ^{u-d} = 4.15(22)(10)_{\rm sys}$. All our results are consistent with phenomenological/experimental values but with larger errors. Lastly, we present a Padé parameterization of the axial, electric and magnetic form factors over the range $0.04< Q^2 <1$ GeV${}^2$ for phenomenological studies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The pion-nucleon sigma term from Lattice QCD

We summarize recent evidence, both from lattice QCD and chiral perturbation theory, that suggests that larger-than-expected excited-state contamination could be the reason for the tension between phenomenological determinations and previous direct lattice-QCD calculations of the pion--nucleon sigma term $σ_{πN}$. In addition, we extend the $χ$PT analysis by calculating the corrections due to including the $Δ(1232)$ resonance as an explicit degree of freedom. This correction is found to be small, thereby corroborating the excited-state effects found in the $Δ$-less calculation and the result for $σ_{πN}$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Why I'm not Answering: Understanding Determinants of Classification of an Abstaining Classifier for Cancer Pathology Reports

Safe deployment of deep learning systems in critical real world applications requires models to make very few mistakes, and only under predictable circumstances. In this work, we address this problem using an abstaining classifier that is tuned to have $>$95% accuracy, and then identify the determinants of abstention using LIME. Essentially, we are training our model to learn the attributes of pathology reports that are likely to lead to incorrect classifications, albeit at the cost of reduced sensitivity. We demonstrate an abstaining classifier in a multitask setting for classifying cancer pathology reports from the NCI SEER cancer registries on six tasks of interest. For these tasks, we reduce the classification error rate by factors of 2--5 by abstaining on 25--45% of the reports. For the specific task of classifying cancer site, we are able to identify metastasis, reports involving lymph nodes, and discussion of multiple cancer sites as responsible for many of the classification mistakes, and observe that the extent and types of mistakes vary systematically with cancer site (e.g., breast, lung, and prostate). When combining across three of the tasks, our model classifies 50% of the reports with an accuracy greater than 95% for three of the six tasks\edit, and greater than 85% for all six tasks on the retained samples. Furthermore, we show that LIME provides a better determinant of classification than measures of word occurrence alone. By combining a deep abstaining classifier with feature identification using LIME, we are able to identify concepts responsible for both correctness and abstention when classifying cancer sites from pathology reports. The improvement of LIME over keyword searches is statistically significant, presumably because words are assessed in context and have been identified as a local determinant of classification.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Moments of nucleon isovector structure functions in $2+1+1$-flavor QCD

We present results on the isovector momentum fraction, $\langle x \rangle_{u-d}$, helicity moment, $\langle x \rangle_{Δu-Δd}$, and the transversity moment, $\langle x \rangle_{δu-δd}$, of the nucleon obtained using nine ensembles of gauge configurations generated by the MILC collaboration using $2+1+1$-flavors of dynamical highly improved staggered quarks (HISQ). The correlation functions are calculated using the Wilson-Clover action and the renormalization of the three operators is carried out nonperturbatively on the lattice in the RI${}^\prime$-MOM scheme. The data have been collected at lattice spacings $a \approx 0.15,\ 0.12,\ 0.09,$ and 0.06 fm and $M_π\approx 310,\ 220$ and 135 MeV, which are used to obtain the physical values using a simultaneous chiral-continuum-finite-volume fit. The final results, in the $\overline{MS}$ scheme at 2 GeV, are $\langle x \rangle_{u-d} = 0.173(14)(07)$, $\langle x \rangle_{Δu-Δd} = 0.213(15)(22)$ and $\langle x \rangle_{δu-δd} = 0.208(19)(24)$, where the first error is the overall analysis uncertainty and the second is an additional systematic uncertainty due to possible residual excited-state contributions. These results are consistent with other recent lattice calculations and phenomenological global fit values.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Nucleon isovector momentum fraction, helicity and transversity moment using Lattice QCD

We present our recent high precision calculations (Phys. Rev. D102 (2020) no.5, 054512 and JHEP 04 (2021) 044, JHEP 21 (2020) 004) of the first moment of nucleon isovector polarized, unpolarized and transversity distributions, i.e., momentum fraction, helicity and transversity moment, respectively. We use the standard method for the calculation of these moments (via matrix elements of twist two operators), and carry out a detailed analysis of the sources of systematic uncertainty, in particular of excited state contributions. Our calculations have been performed using two different lattice setups (Clover-on-HISQ and Clover-on-Clover), each with several ensembles. They give consistent results that are in agreement with global fit analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Space-time symmetric qubit regularization of the asymptotically free two-dimensional O(4) model

We explore if space-time symmetric lattice field theory models with a finite Hilbert space per lattice site can reproduce asymptotic freedom in the two-dimensional $O(4)$ model. We focus on a simple class of such models with a five dimensional local Hilbert space. We demonstrate how even the simplest model reproduces asymptotic freedom within the D-theory formalism but at the cost of increasing the size of the Hilbert space through coupling several layers of a two-dimensional lattice. We then argue that qubit regularization can be viewed as an effective field theory (EFT) even if the continuum limit cannot be reached, as long as we can tune the model close enough to the continuum limit where perturbation theory, or other analytical techniques, become viable. We construct a simple lattice model on a single layer with a four dimensional local Hilbert space that acts like an excellent EFT of the original theory.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Neutron Electric Dipole Moments with Clover Fermions

We present preliminary results for the contributions to the neutron EDM arising from the QCD $θ$-term, the Weinberg three-gluon and the quark chromo-EDM operators from our ongoing lattice calculations using clover valence quarks on the MILC HISQ lattices. We use the gradient-flow technique to smooth the lattices and renormalize the gluonic operators, and use the Schwinger source method to incorporate the quark chromo-EDM interactions in the quark propagator. For the QCD $θ$-term and the Weinberg three-gluon operator, we report results in the gradient-flow scheme from 8 ensembles at four lattice spacings and three pion masses, including 2 physical pion mass ensembles. For the quark chromo-EDM, unrenormalized results are presented at two lattice spacings, $a=0.12$ and $0.09$ fm, and two pion masses, $M_π= 310$ MeV and $220$ MeV.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Nucleon Axial Form Factors from Clover Fermion on 2+1+1-flavor HISQ Lattice

The nucleon axial form factors -- axial $G_A$, induced pseudoscalar $\widetilde{G}_P$ and pseudoscalar $G_P$ -- have displayed large systematics in lattice QCD calculations. The major symptoms were the violation of the partially conserved axial current (PCAC) relation between the three form factors, and the underestimation of the induced pseudoscalar coupling $g_P^\ast$ and the axial charge radius $r_A$ compared to phenomenological estimates. The small $g_P^\ast$ was a consequence of the failure of the pion-pole dominance (PPD) hypothesis, especially at low $M_π^2$. The small charge radius $r_A$ and the underestimate of $g_A$ were related. The dominant systematic responsible is the lack of inclusion of low-energy ($N π$) states that are not manifest in the multiexponential fit to the nucleon two-point correlator. We show that this low-energy state can be determined from the three-point correlator $\langle N A_4 N \rangle $ with the insertion of the temporal component of the axial current $A_4$ within the nucleon state, ie, the strategy labeled $S_{A4}$ [1]. Including this low-energy state in fits to control excited-state contamination (ESC) gives results for $g_A$, $r_A$, and $g_P^\ast$ that are consistent with experimental/phenomenological values. However, the systematic uncertainties, especially in data at small $Q^2$, are now much larger.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Nucleon charges and form factors using clover and HISQ ensembles

We present high statistics ($\mathcal{O}(2\times 10^5)$ measurements) preliminary results on (i) the isovector charges, $g^{u-d}_{A,S,T}$, and form factors, $G^{u-d}_E(Q^2)$, $G^{u-d}_M(Q^2)$, $G^{u-d}_A(Q^2)$, $\widetilde G^{u-d}_P(Q^2)$, $G^{u-d}_P(Q^2)$, on six 2+1-flavor Wilson-clover ensembles generated by the JLab/W&M/LANL/MIT collaboration with lattice parameters given in Table 1. Examples of the impact of using different estimates of the excited state spectra are given for the clover-on-clover data, and as discussed in [1], the biggest difference on including the lower energy (close to $Nπ$ and $Nππ$) states is in the axial channel. (ii) Flavor diagonal axial, tensor and scalar charges, $g^{u,d,s}_{A,S,T}$, are calculated with the clover-on-HISQ formulation using nine 2+1+1-flavor HISQ ensembles generated by the MILC collaboration [2] with lattice parameters given in Table 2. Once finished, the calculations of $g^{u,d,s}_{A,T}$ will update the results given in Refs.[3,4]. The estimates for $g^{u,d,s}_{S}$ and $σ_{Nπ}$ are new. Overall, a large part of the focus is on understanding the excited state contamination (ESC), and the results discussed provide a partial status report on developing defensible analyses strategies that include contributions of possible low-lying excited states to individual nucleon matrix elements.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On Mixup Training: Improved Calibration and Predictive Uncertainty for Deep Neural Networks

Mixup~\cite{zhang2017mixup} is a recently proposed method for training deep neural networks where additional samples are generated during training by convexly combining random pairs of images and their associated labels. While simple to implement, it has been shown to be a surprisingly effective method of data augmentation for image classification: DNNs trained with mixup show noticeable gains in classification performance on a number of image classification benchmarks. In this work, we discuss a hitherto untouched aspect of mixup training -- the calibration and predictive uncertainty of models trained with mixup. We find that DNNs trained with mixup are significantly better calibrated -- i.e., the predicted softmax scores are much better indicators of the actual likelihood of a correct prediction -- than DNNs trained in the regular fashion. We conduct experiments on a number of image classification architectures and datasets -- including large-scale datasets like ImageNet -- and find this to be the case. Additionally, we find that merely mixing features does not result in the same calibration benefit and that the label smoothing in mixup training plays a significant role in improving calibration. Finally, we also observe that mixup-trained DNNs are less prone to over-confident predictions on out-of-distribution and random-noise data. We conclude that the typical overconfidence seen in neural networks, even on in-distribution data is likely a consequence of training with hard labels, suggesting that mixup be employed for classification tasks where predictive uncertainty is a significant concern.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Recent results of nucleon structure & matrix element calculations

A review of recent lattice calculations of nucleon structure and matrix elements of operators in nucleons is presented. It primarily covers developments in the calculation of the matrix elements of the scalar, tensor, pseudo-scalar, axial-vector and vector operators in the ground state of neutrons and protons in the isospin symmetric limit. Some preliminary calculations of the electric dipole moment, the gravitational moments and stress-energy distribution, and the magnetic polarizability are briefly described.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

State preparation and measurement in a quantum simulation of the O(3) sigma model

Recently, Singh and Chandrasekharan showed that fixed points of the non-linear O(3) sigma model can be reproduced near a quantum phase transition of a spin model with just two qubits per lattice site. In a paper by the NuQS collaboration, the proposal is made to simulate such field theories on a quantum computer using the universal properties of a similar model. In this paper, following that direction, we demonstrate how to prepare the ground state of the model from and measure a dynamical quantity of interest, the O(3) Noether charge, on a quantum computer. In particular, we apply Trotter methods to obtain results for the complexity of adiabatic ground state preparation in both the weak-coupling and quantum-critical regimes and use shadow tomography to measure the dynamics of local observables. We then present and analyze a quantum algorithm based on non-unitary randomized simulation methods that may yield an approach suitable for intermediate-term noisy quantum devices.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Axial Vector Form Factors from Lattice QCD that Satisfy the PCAC Relation

Previous lattice QCD calculations of axial vector and pseudoscalar form factors show significant deviation from the partially conserved axial current (PCAC) relation between them. Since the original correlation functions satisfy PCAC, the observed deviations from the operator identity cast doubt on whether all the systematics in the extraction of form factors from the correlation functions are under control. We identify the problematic systematic as a missed excited state, whose energy as a function of the momentum transfer squared, $Q^2$, is determined from the analysis of the 3-point functions themselves. Its mass is much smaller than those of the excited states previously considered and including it impacts the extraction of all the ground state matrix elements. The form factors extracted using these mass/energy gaps satisfy PCAC and other consistency conditions, and validate the pion-pole dominance hypothesis. We also show that the extraction of the axial charge $g_A$ is very sensitive to the value of the mass gaps of the excited states used and current lattice data do not provide an unambiguous determination of these, unlike the $Q^2 \neq 0$ case. To highlight the differences and improvement between the conventional versus the new analysis strategy, we present a comparison of results obtained on a physical pion mass ensemble at $a\approx 0.0871\,\mathrm{fm}$. With the new strategy, we find $g_A = 1.30(6)$. A very significant improvement over previous lattice results is found for the axial charge radius $r_A = 0.74(6)\,\mathrm{fm}$, extracted using the $z$-expansion to parameterize the $Q^2$ behavior of $G_A(Q^2)$, and $g_P^\ast = 8.06(44)$ obtained using the pion pole-dominance ansatz to fit the $Q^2$ behavior of the induced pseudoscalar form factor $\widetilde{G}_P(Q^2)$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Nucleon Electromagnetic Form Factors in the Continuum Limit from 2+1+1-flavor Lattice QCD

Results are presented for the nucleon isovector electromagnetic form factors using 11 ensembles generated by the MILC collaboration using the 2+1+1-flavors HISQ action. They span 4 lattice spacings $a \sim$ 0.06, 0.09, 0.12 and 0.15~fm and 3 values of $M_π\sim 135, 225$ and 315 MeV. High-statistics estimates are used to perform a simultaneous extrapolation in the lattice spacing, lattice volume and light-quark masses. The $Q^2$ dependence over the range 0.05-1.4 ${\rm GeV}^2$ is investigated using both the $z$-expansion and the dipole form. Final $z$-expansion estimates for the isovector r.m.s. radius are $r_E = 0.769(27)(30)$ fm $r_M = 0.671(48)(76)$ fm and $μ^{p-n} = 3.939(86)(138)$ Bohr magneton. The first error is the combined uncertainty from the leading-order analysis, and the second is an estimate of the additional uncertainty due to using the leading order chiral-continuum-finite-volume fits. The dipole estimates, $r_E = 0.765(11)(8)$ fm, $r_M = 0.704(21)(29)$ fm and $μ^{p-n} = 3.975(84)(125)$, are consistent with those from the $z$-expansion but with smaller errors. Our analysis highlights three points. First, all data from the eleven ensembles and existing lattice data on, or close to, physical mass ensembles from other collaborations collapses more clearly onto a single curve when plotted versus $Q^2/M_N^2$ as compared to $Q^2$ with the scale set by quantities other than $M_N$. The difference between these two analyses is indicative of discretization errors, some of which presumably cancel when the data are plotted versus $Q^2/M_N^2$. Second, the size of the remaining deviation of this common curve from the Kelly curve is small and can be accounted for by statistical and possible systematic uncertainties. Third, to improve lattice estimates, high statistics data for $Q^2 < 0.1$ ${\rm GeV}^2$ are needed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Likelihood Ratio Detector for Identifying Within-Perimeter Computer Network Attacks

The rapid detection of attackers within firewalls of enterprise computer net- works is of paramount importance. Anomaly detectors address this problem by quantifying deviations from baseline statistical models of normal network behav- ior and signaling an intrusion when the observed data deviates significantly from the baseline model. However, many anomaly detectors do not take into account plausible attacker behavior. As a result, anomaly detectors are prone to a large number of false positives due to unusual but benign activity. This paper first in- troduces a stochastic model of attacker behavior which is motivated by real world attacker traversal. Then, we develop a likelihood ratio detector that compares the probability of observed network behavior under normal conditions against the case when an attacker has possibly compromised a subset of hosts within the network. Since the likelihood ratio detector requires integrating over the time each host be- comes compromised, we illustrate how to use Monte Carlo methods to compute the requisite integral. We then present Receiver Operating Characteristic (ROC) curves for various network parameterizations that show for any rate of true posi- tives, the rate of false positives for the likelihood ratio detector is no higher than that of a simple anomaly detector and is often lower. We conclude by demon- strating the superiority of the proposed likelihood ratio detector when the network topologies and parameterizations are extracted from real-world networks.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Iso-vector and Iso-scalar Tensor Charges of the Nucleon from Lattice QCD

We present results for the iso-vector and flavor diagonal tensor charges $g^{u-d}_T$, $g^{u}_T$, $g^{d}_T$, and $g^{s}_T$ needed to probe novel tensor interactions at the TeV scale in neutron and nuclear $β$-decays and the contribution of the quark electric dipole moment (EDM) to the neutron EDM. The lattice QCD calculations were done using nine ensembles of gauge configurations generated by the MILC collaboration using the HISQ action with 2+1+1 dynamical flavors. These ensembles span three lattice spacings $a \approx 0.06, 0.09$ and $0.12 $ fm and three quark masses corresponding to the pion masses $M_π\approx 130, 220$ and $310 $ MeV. Using estimates from these ensembles, we quantify all systematic uncertainties and perform a simultaneous extrapolation in the lattice spacing, volume and light quark masses for the connected contributions. The final estimates of the connected nucleon (proton) tensor charge for the iso-vector combination is $g_T^{u-d} = 1.020(76) $ in the $\bar{\text{MS}}$ scheme at $2$ GeV. The additional disconnected quark loop contributions needed for the flavor-diagonal matrix elements are calculated using a stochastic estimator employing the truncated solver method with the all-mode-averaging technique. We find that the size of the disconnected contribution is smaller than the statistical error in the connected contribution. This allows us to bound the disconnected contribution and include it as an additional uncertainty in the flavor-diagonal charges. After a continuum extrapolation, we find $g_T^{u} = 0.774(66) $, $g_T^{d} = -0.233(28) $ and $g_T^{u+d} = 0.541(67) $. The strangeness tensor charge, that can make a significant contribution to the neutron EDM due to the large ratio $m_s/m_{u,d}$, is $g_T^{s}=0.008(9)$ in the continuum limit.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Neutron Electric Dipole Moment and Tensor Charges from Lattice QCD

We present Lattice QCD results on the neutron tensor charges including, for the first time, a simultaneous extrapolation in the lattice spacing, volume, and light quark masses to the physical point in the continuum limit. We find that the "disconnected" contribution is smaller than the statistical error in the "connected" contribution. Our estimates in the $\bar{\text{MS}}$ scheme at $2$ GeV, including all systematics, are $g_T^{d-u}=1.020(76)$, $g_T^d = 0.774(66)$, $g_T^u = - 0.233(28)$, and $g_T^s = 0.008(9)$. The flavor diagonal charges determine the size of the neutron electric dipole moment (EDM) induced by quark EDMs that are generated in many new scenarios of CP-violation beyond the Standard Model (BSM). We use our results to derive model-independent bounds on the EDMs of light quarks and update the EDM phenomenology in split Supersymmetry with gaugino mass unification, finding a stringent upper bound of $d_n < 4 \times 10^{-28} \, e$ cm for the neutron EDM in this scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Neutron Electric Dipole Moment from quark Chromoelectric Dipole Moment

The connection between a regularization-independent symmetric momentum substraction (RI-$\tilde{\rm S}$MOM) and the $\overline{\rm MS}$ scheme for the quark chromo EDM operators is discussed. A method for evaluating the neutron EDM from quark chromoEDM is described. A preliminary study of the signal in the matrix element using clover quarks on a highly improved staggered quark (HISQ) ensemble is shown.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nucleon Charges, Form-factors and Neutron EDM

We present an update of our analysis of statistical and systematic errors in the calculation of iso-vector scalar, axial and tensor charges of the nucleon. The calculations are done using $N_f=2+1+1$ flavor HISQ ensembles generated by the MILC Collaboration at three values of the lattice spacing ($a=0.12,\ 0.09,$ and $0.06$ fm) and three values of the quark mass ($M_π\approx 310,\ 220$ and $130$ MeV); and clover fermions for calculating the correlation functions, i.e., we use a clover-on-HISQ lattice formulation. The all-mode-averaging method allows us to increase the statistics by a factor of eight for the same computational cost leading to a better understanding of and control over excited state contamination. Our current results, after extrapolation to the continuum limit and physical pion mass are $g_A^{u-d} = 1.21(3)$, $g_T^{u-d} = 1.005(59)$ and $g_S^{u-d} = 0.95(12) $. Further checks of control over all systematic errors, especially in $g_A^{u-d}$, are still being performed. Using results for the flavor-diagonal charges, $g_T^{u} $, $g_T^{d} $ and $g_T^{s}$, we analyze contributions of the quark electric dipole moment to the neutron EDM and the consequences for split SUSY model.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Quark Chromoelectric Dipole Moment Contribution to the Neutron Electric Dipole Moment

The quark chromo-electric dipole moment operator and the pseudo-scalar fermion bilinear with which it mixes under renormalization can both be included in a calculation of the electromagnetic form factor of the nucleon using the Schwinger source method. A preliminary calculation of these operators using clover quarks on HISQ lattices generated by MILC collaboration will be presented showing the quality of the signal in the correlators necessary for calculating the neutron electric dipole moment.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Tuning the hopping parameter in the Oktay-Kronfeld action for charm and bottom quarks on a MILC HISQ ensemble

The first step in the calculation of semi-leptonic form factors in the decay of heavy mesons is the tuning of the hopping parameter $κ$ for the charm and bottom quark masses. Results for the Oktay-Kronfeld (OK) action are presented for one $N_f=2+1+1$ HISQ ensemble generated by the MILC collaboration at $a\approx 0.12\,\mathrm{fm}$ and $M_π\approx 310$ MeV. Estimates of hyperfine splitting of heavy-light and heavy-heavy mesons are presented and the inconsistency parameter is evaluated.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Worldline approach for numerical computation of electromagnetic Casimir energies. I. Scalar field coupled to magnetodielectric media

We present a worldline method for the calculation of Casimir energies for scalar fields coupled to magnetodielectric media. The scalar model we consider may be applied in arbitrary geometries, and it corresponds exactly to one polarization of the electromagnetic field in planar layered media. Starting from the field theory for electromagnetism, we work with the two decoupled polarizations in planar media and develop worldline path integrals, which represent the two polarizations separately, for computing both Casimir and Casimir-Polder potentials. We then show analytically that the path integrals for the transverse-electric (TE) polarization coupled to a dielectric medium converge to the proper solutions in certain special cases, including the Casimir-Polder potential of an atom near a planar interface, and the Casimir energy due to two planar interfaces. We also evaluate the path integrals numerically via Monte-Carlo path-averaging for these cases, studying the convergence and performance of the resulting computational techniques. While these scalar methods are only exact in particular geometries, they may serve as an approximation for Casimir energies for the vector electromagnetic field in other geometries.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Disconnected Quark Loop Contributions to Nucleon Structure

We calculate the disconnected contribution to isoscalar nucleon charges for scalar, axial and tensor channels of light and strange quarks. The calculation has been done with the Clover valence quarks on the MILC $N_f=2+1+1$ HISQ lattices whose light quark masses corresponding to the pion masses of 305 MeV and 217 MeV at $a \approx 0.12$ fm and 312 MeV at $a \approx 0.09$ fm. All-mode-averaging technique is used for the evaluation two-point functions. Disconnected quark loops are estimated by using the truncated solver method with Gaussian random noise sources. Contamination from the excited states is removed by fitting the results of various source-sink separations and operator insertions to the formula including up to the first excited state, simultaneously.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Longitudinal Antigenic Sequences and Sites from Intra-Host Evolution (LASSIE) Identifies Immune-Selected HIV Variants

Within-host genetic sequencing from samples collected over time provides a dynamic view of how viruses evade host immunity. Immune-driven mutations might stimulate neutralization breadth by selecting antibodies adapted to cycles of immune escape that generate within-subject epitope diversity. Comprehensive identification of immune-escape mutations is experimentally and computationally challenging. With current technology, many more viral sequences can readily be obtained than can be tested for binding and neutralization, making down-selection necessary. Typically, this is done manually, by picking variants that represent different time-points and branches on a phylogenetic tree. Such strategies are likely to miss many relevant mutations and combinations of mutations, and to be redundant for other mutations. Longitudinal Antigenic Sequences and Sites from Intrahost Evolution (LASSIE) uses transmitted-founder loss to identify virus "hot-spots" under putative immune selection and chooses sequences that represent recurrent mutations in selected sites. LASSIE favors earliest sequences in which mutations arise. With well-characterized longitudinal Env sequences, we confirmed selected sites were concentrated in antibody contacts and selected sequences represented diverse antigenic phenotypes. Practical applications include rapidly identifying immune targets under selective pressure within a subject, selecting minimal sets of reagents for immunological assays that characterize evolving antibody responses, and for immunogens in polyvalent "cocktail" vaccines.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the universal structure of human lexical semantics

How universal is human conceptual structure? The way concepts are organized in the human brain may reflect distinct features of cultural, historical, and environmental background in addition to properties universal to human cognition. Semantics, or meaning expressed through language, provides direct access to the underlying conceptual structure, but meaning is notoriously difficult to measure, let alone parameterize. Here we provide an empirical measure of semantic proximity between concepts using cross-linguistic dictionaries. Across languages carefully selected from a phylogenetically and geographically stratified sample of genera, translations of words reveal cases where a particular language uses a single polysemous word to express concepts represented by distinct words in another. We use the frequency of polysemies linking two concepts as a measure of their semantic proximity, and represent the pattern of such linkages by a weighted network. This network is highly uneven and fragmented: certain concepts are far more prone to polysemy than others, and there emerge naturally interpretable clusters loosely connected to each other. Statistical analysis shows such structural properties are consistent across different language groups, largely independent of geography, environment, and literacy. It is therefore possible to conclude the conceptual structure connecting basic vocabulary studied is primarily due to universal features of human cognition and language use.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Precision calculations of nucleon charges $g_A$, $g_S$, $g_T$

We present a detailed analysis of statistical and systematic errors in the calculation of matrix elements of iso-vector scalar, axial and tensor charges between a neutron and a proton state. These analyses are being done on dynamical $N_f=2+1+1$ HISQ configurations generated by the MILC Collaboration using valence clover fermions. Using ensembles at three values of the lattice spacing ($a=0.12,\ 0.09,$ and $0.06$ fm) and three values of the quark mass ($M_π\approx 310,\ 220$ and $130$ MeV) we find that the estimates of the tensor charge are stable and it can be extracted with $5\%$ precision with O(10,000) measurements. We also find that higher statistics are needed to resolve the various uncertainties in the calculation of $g_A$ and improve the signal in $g_S$, which with present data has large errors. A brief status report on the mixing and renormalization of novel operators contributing to nEDM is also given.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Neutron Electric Dipole Moments from Beyond the Standard Model Physics

Neutron Electric Dipole Moment (nEDM), a generic feature of CP-violation, is predicted to be very small in the Standard Model, but can be much larger in most extensions of the model. In this talk, I will discuss the classification of the CP violating operators up to dimension 6 that can give rise to nEDM, and then describe the mixing and renormalization structure of the operators of dimension 5 and lower in both dimensional and cutoff regularizations in general terms. Finally I will describe how to connect the dimension 5 operators, in particular, the Chromoelectric Dipole Moment of the quarks, between MSbar scheme and a Regularization Independent prescription in the chiral limit.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nucleon Charges and Electromagnetic Form Factors from 2+1+1-Flavor Lattice QCD

We present lattice-QCD results on the nucleon isovector axial, scalar and tensor charges, the isovector electromagnetic Dirac and Pauli form factors, and the connected parts of the isoscalar charges. The calculations have been done using two ensembles of HISQ lattices generated by the MILC Collaboration with 2+1+1 dynamical flavors at a lattice spacing of 0.12 fm and with light-quark masses corresponding to pions with masses 310 and 220 MeV. We perform a systematic study including excited-state degrees of freedom and examine the dependence of the extracted nucleon matrix elements on source-sink separation. This study demonstrates with high-statistics data that including excited-state contributions and generating data at multiple separations is necessary to remove contamination that would otherwise lead to systematic error. We also determine the renormalization constants of the associated quark bilinear operators in the RI-sMOM scheme and make comparisons of our renormalized results with previous dynamical-lattice calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Probing TeV scale physics in precision UCN decays

We present the calculation of matrix elements of iso-vector scalar, axial and tensor charges between a neutron and a proton state on dynamical $N_f=2+1+1$ HISQ configurations generated by the MILC Collaboration using valence clover fermions. These matrix elements are needed to probe novel scalar and tensor interactions in neutron beta-decay that can arise in extensions to the Standard Model at the TeV scale. Results are presented at one value of the lattice spacing, $a=0.12$ fm, and two values of light quarks corresponding to $M_π=310$ and $220$ MeV. We discuss two sources of systematic errors, contribution of excited states to these matrix elements and the renormalization constants, and the efficacy of methods used to control them.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The equation of state in (2+1)-flavor QCD

We present results for the equation of state in (2+1)-flavor QCD using the highly improved staggered quark action and lattices with temporal extent $N_τ=6,~8,~10$, and $12$. We show that these data can be reliably extrapolated to the continuum limit and obtain a number of thermodynamic quantities and the speed of sound in the temperature range $(130-400)$ MeV. We compare our results with previous calculations, and provide an analytic parameterization of the pressure, from which other thermodynamic quantities can be calculated, for use in phenomenology. We show that the energy density in the crossover region, $145~ {\rm MeV} \leq T \leq 163$ MeV, defined by the chiral transition, is $ε_c=(0.18-0.5)~{\rm GeV}/{\rm fm}^3$, $i.e.$, $(1.2-3.1)\ ε_{\rm nuclear}$. At high temperatures, we compare our results with resummed and dimensionally reduced perturbation theory calculations. As a byproduct of our analyses, we obtain the values of the scale parameters $r_0$ from the static quark potential and $w_0$ from the gradient flow.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The QCD phase transition with physical-mass, chiral quarks

We report on the first lattice calculation of the QCD phase transition using chiral fermions at physical values of the quark masses. This calculation uses 2+1 quark flavors, spatial volumes between (4 fm$)^3$ and (11 fm$)^3$ and temperatures between 139 and 196 MeV . Each temperature was calculated using a single lattice spacing corresponding to a temporal Euclidean extent of $N_t=8$. The disconnected chiral susceptibility, $χ_{\rm disc}$ shows a pronounced peak whose position and height depend sensitively on the quark mass. We find no metastability in the region of the peak and a peak height which does not change when a 5 fm spatial extent is increased to 10 fm. Each result is strong evidence that the QCD ``phase transition'' is not first order but a continuous cross-over for $m_π=135$ MeV. The peak location determines a pseudo-critical temperature $T_c = 155(1)(8)$ MeV. Chiral $SU(2)_L\times SU(2)_R$ symmetry is fully restored above 164 MeV, but anomalous $U(1)_A$ symmetry breaking is non-zero above $T_c$ and vanishes as $T$ is increased to 196 MeV.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Fluctuations and Correlations of net baryon number, electric charge, and strangeness: A comparison of lattice QCD results with the hadron resonance gas model

We calculate the quadratic fluctuations of net baryon number, electric charge and strangeness as well as correlations among these conserved charges in (2+1)-flavor lattice QCD at zero chemical potential. Results are obtained using calculations with tree level improved gauge and the highly improved staggered quark (HISQ) actions with almost physical light and strange quark masses at three different values of the lattice cut-off. Our choice of parameters corresponds to a value of 160 MeV for the lightest pseudo scalar Goldstone mass and a physical value of the kaon mass. The three diagonal charge susceptibilities and the correlations among conserved charges have been extrapolated to the continuum limit in the temperature interval 150 MeV <T < 250 MeV. We compare our results with the hadron resonance gas (HRG) model calculations and find agreement with HRG model results only for temperatures T<= 150 MeV. We observe significant deviations in the temperature range 160 MeV < T < 170 MeV and qualitative differences in the behavior of the three conserved charge sectors. At T < 160 MeV quadratic net baryon number fluctuations in QCD agree with HRG model calculations while, the net electric charge fluctuations in QCD are about 10% smaller and net strangeness fluctuations are about 20% larger. These findings are relevant to the discussion of freeze-out conditions in relativistic heavy ion collisions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Probing novel TeV physics through precision calculations of scalar and tensor charges of the nucleon

We present an update on the calculation of matrix elements of iso-vector scalar, axial and tensor charges between a neutron and a proton state. These matrix elements are needed to probe novel scalar and tensor interactions in neutron beta-decay that can arise in extensions of the Standard Model at the TeV scale. Our calculations are being done using valence clover fermions on dynamical N_f=2+1+1 HISQ configurations generated by the MILC Collaboration. We provide preliminary estimates of the dependence of these matrix elements on the light quark masses, lattice spacing, and the time separation between the source and sink of the nucleons. We also find that the renormalization constants calculated using the RI-sMOM scheme are close to unity for the HYP smeared HISQ lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Probing TeV scale physics via ultra cold neutron decays and calculating non-standard baryon matrix elements

We motivate undertaking precision analyses of neutron decays to look for signatures of new scalar and tensor interactions that can arise in extensions of the Standard Model at the TeV scale. The key ingrediant needed to connect experimental data with theoretical analysis are high-precision calculations of matrix elements of isovector bilinear operators between the decaying neutron and final state proton. We describe the status of our Lattice QCD program of using valence clover fermions on dynamical N_f=2+1+1 HISQ configurations generated by the MILC Collaboration. On the theoretical side we use the effective field theory method and provide both model independent and dependent analyses to obtain bounds on possible scalar and tensor interactions, both from low energy experiments and LHC data.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The chiral transition and U(1)_A symmetry restoration from lattice QCD using Domain Wall Fermions

We present results on both the restoration of the spontaneously broken chiral symmetry and the effective restoration of the anomalously broken U(1)_A symmetry in finite temperature QCD at zero chemical potential using lattice QCD. We employ domain wall fermions on lattices with fixed temporal extent N_τ= 8 and spatial extent N_σ= 16 in a temperature range of T = 139 - 195 MeV, corresponding to lattice spacings of a \approx 0.12 - 0.18 fm. In these calculations, we include two degenerate light quarks and a strange quark at fixed pion mass m_π= 200 MeV. The strange quark mass is set near its physical value. We also present results from a second set of finite temperature gauge configurations at the same volume and temporal extent with slightly heavier pion mass. To study chiral symmetry restoration, we calculate the chiral condensate, the disconnected chiral susceptibility, and susceptibilities in several meson channels of different quantum numbers. To study U(1)_A restoration, we calculate spatial correlators in the scalar and pseudo-scalar channels, as well as the corresponding susceptibilities. Furthermore, we also show results for the eigenvalue spectrum of the Dirac operator as a function of temperature, which can be connected to both U(1)_A and chiral symmetry restoration via Banks-Casher relations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Theoretical Bounds on New Four-Fermion Interactions and TeV Scale Physics

The standard model weak interactions can be described by four-fermion V-A operators at low energies. New physics at the TeV scale can, however, generate the other Lorentz structures. In this talk, we review the constraints on such interactions from nuclear and hadronic decays, as well as from collider searches. Currently the most stringent bounds come from the analysis of the 0+ to 0+ nuclear and the pi to e nu gamma radiative pion decays. In the near future, the ultracold neutron beta decay experiments and the direct LHC measurements will compete in setting the most stringent bounds, provided, however, that the neutron-to-proton non-perturbative transition matrix elements can be calculated to a level of 10-20% accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Probing Novel Scalar and Tensor Interactions from (Ultra)Cold Neutrons to the LHC

Scalar and tensor interactions were once competitors to the now well-established V-A structure of the Standard Model weak interactions. We revisit these interactions and survey constraints from low-energy probes (neutron, nuclear, and pion decays) as well as collider searches. Currently, the most stringent limit on scalar and tensor interactions arise from 0+ -> 0+ nuclear decays and the radiative pion decay pi -> e nu gamma, respectively. For the future, we find that upcoming neutron beta decay and LHC measurements will compete in setting the most stringent bounds. For neutron beta decay, we demonstrate the importance of lattice computations of the neutron-to-proton matrix elements to setting limits on these interactions, and provide the first lattice estimate of the scalar charge and a new average of existing results for the tensor charge. Data taken at the LHC is currently probing these interactions at the 10^-2 level (relative to the standard weak interactions), with the potential to reach the 10^-3 level. We show that, with some theoretical assumptions, the discovery of a charged spin-0 resonance decaying to an electron and missing energy implies a lower limit on the strength of scalar interactions probed at low energy.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Nonlinear Quantum Dynamics

The vast majority of the literature dealing with quantum dynamics is concerned with linear evolution of the wave function or the density matrix. A complete dynamical description requires a full understanding of the evolution of measured quantum systems, necessary to explain actual experimental results. The dynamics of such systems is intrinsically nonlinear even at the level of distribution functions, both classically as well as quantum mechanically. Aside from being physically more complete, this treatment reveals the existence of dynamical regimes, such as chaos, that have no counterpart in the linear case. Here, we present a short introductory review of some of these aspects, with a few illustrative results and examples.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Testing improved staggered fermions with $m_s$ and $B_K$

We study the improvement of staggered fermions using hypercubically smeared (HYP) links. We calculate the strange quark mass and the kaon B-parameter, $B_K$, in quenched QCD on a $16^3 \times 64$ lattice at $β=6.0$. We find $m_s(\bar{\rm MS},2 {\rm GeV})=101.2\pm1.3\pm4 $MeV and $B_K(\bar{\rm MS},2 {\rm GeV}) = 0.578 \pm 0.018\pm 0.042$, where the first error is from statistics and fitting, and the second from using one-loop matching factors. The scale ($1/a=1.95$GeV) is set by $M_ρ$, and $m_s$ is determined using the kaon mass. Comparing to quenched results obtained using unimproved staggered fermions and other discretizations, we argue that the size of discretization errors in $B_K$ is substantially reduced by improvement.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Scaling behavior of improvement and renormalization constants

This talk summarizes results for all the scale independent renormalization constants for bilinear currents ($Z_A$, $Z_V$, and $Z_S/Z_P$), the improvement constants ($c_A$, $c_V$, and $c_T$), the quark mass dependence of $Z_O$, and the coefficients of the equation of motion operators for $O(a)$ improved lattice QCD. Using data at $β=6.0$, 6.2 and 6.4 we study the scaling behavior of these quantities and quantify residual discretization errors.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Non-perturbative Renormalization Constants using Ward Identities

We extend the application of axial Ward identities to calculate $b_A, b_P$ and $b_T$, coefficients that give the mass dependence of the renormalization constants of the corresponding bilinear operators in the quenched theory. The extension relies on using operators with non-degenerate quark masses. It allows a complete determination of the O(a) improvement coefficients for bilinears in the quenched approximation using Ward Identities alone. Only the scale dependent normalization constants $Z_P^0$ (or $Z_S^0$) and $Z_T$ are undetermined. We present results of a pilot numerical study using hadronic correlators.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

A pot-pourri of results in QCD from large lattice simulations on the CM5

We present a status report on simulations being done on $32^3 \times 64$ lattices at $β= 6.0$ using quenched Wilson fermions. Results for the spectrum, decay constants, the kaon B-parameter $B_K$, and semi-leptonic form factors are given based on the current statistical sample of 28 configurations. These ``Grand Challenge'' calculations are being done on the CM5 at Los Alamos in collaboration with G. Kilcup, S. Sharpe and P. Tamayo. We end with a brief statement of code performance.

0 citations0 reviewsNo code linkedOpen access