Featured Researches

High Energy Physics Lattice

A lattice investigation of exotic tetraquark channels

We perform annf=2+1lattice study of a number of channels where past claims exist in the literature for the existence of strong-interaction-stable light-heavy tetraquarks. We find no evidence for any such deeply-bound states, beyond theJP=1+,I=0udb¯b¯andI=1/2lsb¯b¯states already identified in earlier lattice studies. We also describe a number of systematic improvements to our previous lattice studies, including working with largermπLto better suppress possible finite volume effects, employing extended sinks to better control excited-state contamination, and expanding the number of operators used in the GEVP analyses. Our results also allow us to rule out several phenomenological models which predict significant tetraquark binding in channels where no such binding is found.

Read more
High Energy Physics Lattice

A model-independent framework for determining finite-volume effects of spatially nonlocal operators

We present a model-independent framework to determine finite-volume corrections of matrix elements of spatially-separated current-current operators. We define these matrix elements in terms of Compton-like amplitudes, i.e. amplitudes coupling single-particle states via two current insertions. We show that the infrared behavior of these matrix elements is dominated by the single-particle pole, which is approximated by the elastic form factors of the lowest-lying hadron. Therefore, given lattice data on the relevant elastic form factors, the finite-volume effects can be estimated non-perturbatively and without recourse to effective field theories. For illustration purposes, we investigate the implications of the proposed formalism for a class of scalar theories in two and four dimensions.

Read more
High Energy Physics Lattice

A multigrid accelerated eigensolver for the Hermitian Wilson-Dirac operator in lattice QCD

Eigenvalues of the Hermitian Wilson-Dirac operator are of special interest in several lattice QCD simulations, e.g., for noise reduction when evaluating all-to-all propagators. In this paper we present a Davidson-type eigensolver that utilizes the structural properties of the Hermitian Wilson-Dirac operatorQto compute eigenpairs of this operator corresponding to small eigenvalues. The main idea is to exploit a synergy between the (outer) eigensolver and its (inner) iterative scheme which solves shifted linear systems. This is achieved by adapting the multigrid DD-αAMG algorithm to a solver for shifted systems involving the Hermitian Wilson-Dirac operator. We demonstrate that updating the coarse grid operator using eigenvector information obtained in the course of the generalized Davidson method is crucial to achieve good performance when calculating many eigenpairs, as our study of the local coherence shows. We compare our method with the commonly used software-packages PARPACK and PRIMME in numerical tests, where we are able to achieve significant improvements, with speed-ups of up to one order of magnitude and a near-linear scaling with respect to the number of eigenvalues. For illustration we compare the distribution of the small eigenvalues ofQon a64×323lattice with what is predicted by the Banks-Casher relation in the infinite volume limit.

Read more
High Energy Physics Lattice

A novel computational paradigm for a precise determination of the hadronic contribution to(gμ??)from lattice QCD

The hadronic contribution to the muon anomalous magnetic momentaμ=(gμ??)/2has to be determined at the per-mille level for the Standard Model prediction to match the expected final uncertainty of the ongoing E989 experiment. That is 3 times better than the current precision from the dispersive approach, and 5-15 times smaller than the uncertainty based on the purely theoretical determinations from lattice QCD. So far the stumbling-block is the large statistical error in the Monte Carlo evaluation of the required correlation functions which can hardly be tamed by brute force. In this talk we present our proposal to solve this problem by multi-level Monte Carlo integration, a technique which reduces the variance of correlators exponentially in the distance of the fields. We report the results of our feasibility tests for the computation of the Hadronic Vacuum Polarization on a lattice with a linear extension of 3~fm, a spacing of 0.065 fm, and a pion mass of 270 MeV. Indeed the two-level integration makes the contribution to the statistical error from long-distances de-facto negligible by accelerating its inverse scaling with the cost of the simulation. These findings establish multi-level Monte Carlo as a solid and efficient method for a precise lattice determination of the hadronic contribution toaμ.

Read more
High Energy Physics Lattice

A physicist-friendly reformulation of the Atiyah-Patodi-Singer index and its mathematical justification

The Atiyah-Patodi-Singer index theorem describes the bulk-edge correspondence of symmetry protected topological insulators. The mathematical setup for this theorem is, however, not directly related to the physical fermion system, as it imposes on the fermion fields a non-local and unnatural boundary condition known as the "APS boundary condition" by hand. In 2017, we showed that the same integer as the APS index can be obtained from theηinvariant of the domain-wall Dirac operator. Recently we gave a mathematical proof that the equivalence is not a coincidence but generally true. In this contribution to the proceedings of LATTICE 2019, we try to explain the whole story in a physicist-friendly way.

Read more
High Energy Physics Lattice

Ab-initio calculation of the proton and the neutron's scalar couplings for new physics searches

Many low-energy, particle-physics experiments seek to reveal new fundamental physics by searching for very rare scattering events on atomic nuclei. The interpretation of their results requires quantifying the non-linear effects of the strong interaction on the spin-independent couplings of this new physics to protons and neutrons. Here we present a fully-controlled, ab-initio calculation of these couplings to the quarks within those constituents of nuclei. We use lattice quantum chromodynamics computations for the four lightest species of quarks and heavy-quark expansions for the remaining two. We determine each of the six quark contributions with an accuracy better than 15%. Our results are especially important for guiding and interpreting experimental searches for our universe's dark matter.

Read more
High Energy Physics Lattice

Adding machine learning within Hamiltonians: Renormalization group transformations, symmetry breaking and restoration

We present a physical interpretation of machine learning functions, opening up the possibility to control properties of statistical systems via the inclusion of these functions in Hamiltonians. In particular, we include the predictive function of a neural network, designed for phase classification, as a conjugate variable coupled to an external field within the Hamiltonian of a system. Results in the two-dimensional Ising model evidence that the field can induce an order-disorder phase transition by breaking or restoring the symmetry, in contrast with the field of the conventional order parameter which causes explicit symmetry breaking. The critical behavior is then studied by proposing a Hamiltonian-agnostic reweighting approach and forming a renormalization group mapping on quantities derived from the neural network. Accurate estimates of the critical point and of the critical exponents related to the operators that govern the divergence of the correlation length are provided. We conclude by discussing how the method provides an essential step toward bridging machine learning and physics.

Read more
High Energy Physics Lattice

Alternating Descent Method for Gauge Cooling of Complex Langevin Simulations

We study the gauge cooling technique for the complex Langevin method applied to the computation in lattice quantum chromodynamics. We propose a new solver of the minimization problem that optimizes the gauge, which does not include any parameter in each iteration, and shows better performance than the classical gradient descent method especially when the lattice size is large. Two numerical tests are carried out to show the effectiveness of the new algorithm.

Read more
High Energy Physics Lattice

Alternative derivation of the relativistic three-particle quantization condition

We present a simplified derivation of the relativistic three-particle quantization condition for identical, spinless particles described by a generic relativistic field theory satisfying aZ2symmetry. The simplification is afforded by using a three-particle quasilocal K matrix that is not fully symmetrized,K˜(u,u)df,3, and makes extensive use of time-ordered perturbation theory (TOPT). We obtain a new form of the quantization condition. This new form can then be related algebraically to the standard quantization condition, which depends on a fully symmetric three-particle K matrix,Kdf,3. The new derivation is fully explicit, allowing, for example, a closed-form expression forKdf,3to be given in terms of TOPT amplitudes. The new form of the quantization condition is similar in structure to that obtained in the "finite-volume unitarity" approach, and in a companion paper we make this connection concrete. Our simplified approach should also allow a more straightforward generalization of the quantization condition to nondegenerate particles, and perhaps also to more than three particles.

Read more
High Energy Physics Lattice

An analysis of systematic effects in finite size scaling studies using the gradient flow

We propose a new strategy for the determination of the step scaling functionσ(u)in finite size scaling studies using the Gradient Flow. In this approach the determination ofσ(u)is broken in two pieces: a change of the flow time at fixed physical size, and a change of the size of the system at fixed flow time. Using both perturbative arguments and a set of simulations in the pure gauge theory we show that this approach leads to a better control over the continuum extrapolations. Following this new proposal we determine the running coupling at high energies in the pure gauge theory and re-examine the determination of theΛ-parameter, with special care on the perturbative truncation uncertainties.

Read more

Ready to get started?

Join us today