Source author record

Tim Harris

Tim Harris 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

14works
7topics
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

14 published item(s)

preprint2022arXiv

Non-perturbative thermal QCD at all temperatures: the case of mesonic screening masses

We present a strategy based on the step-scaling technique to study non-perturbatively thermal QCD up to very high temperatures. As a first concrete application, we compute the flavour non-singlet meson screening masses at 12 temperatures covering the range from $T \sim 1$ GeV up to $\sim 160$ GeV in the theory with three massless quarks. The calculation is carried out by Monte Carlo simulations on the lattice by considering large spatial extensions in order to have negligible finite volume effects. For each temperature we have simulated 3 or 4 values of the lattice spacing, so as to perform the continuum limit extrapolation with confidence at a few permille accuracy. Chiral symmetry restoration manifests itself in our results through the degeneracy of the vector and the axial vector channels and of the scalar and the pseudoscalar ones. In the entire range of temperatures explored, the meson screening masses deviate from the free theory result, $2 πT$, by at most a few percent. These deviations, however, cannot be explained by the known leading term in the QCD coupling constant $g$ up to the highest temperature, where other contributions are still very relevant. In particular the vector-pseudoscalar mass splitting turns out to be of $O(g^4)$ in the entire range explored, and it remains clearly visible up to the highest temperature, where the two screening masses are still significantly different within our numerical precision. The pattern of different contributions that we have found explains why it has been difficult in the past to match non-perturbative lattice results at $T \sim 1$ GeV with the analytic behaviour at asymptotically high temperatures.

preprint2021arXiv

A novel computational paradigm for a precise determination of the hadronic contribution to $(g_μ-2)$ from lattice QCD

The hadronic contribution to the muon anomalous magnetic moment $a_μ=(g_μ-2)/2$ has 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 to $a_μ$.

preprint2021arXiv

Vacuum correlators at short distances from lattice QCD

Non-perturbatively computing the hadronic vacuum polarization at large photon virtualities and making contact with perturbation theory enables a precision determination of the electromagnetic coupling at the $Z$ pole, which enters global electroweak fits. In order to achieve this goal ab initio using lattice QCD, one faces the challenge that, at the short distances which dominate the observable, discretization errors are hard to control. Here we address challenges of this type with the help of static screening correlators in the high-temperature phase of QCD, yet without incurring any bias. The idea is motivated by the observations that (a) the cost of high-temperature simulations is typically much lower than their vacuum counterpart, and (b) at distances $x_3$ far below the inverse temperature $1/T$, the operator-product expansion guarantees the thermal correlator of two local currents to deviate from the vacuum correlator by a relative amount that is power-suppressed in $(x_3\:T)$. The method is first investigated in lattice perturbation theory, where we point out the appearance of an O$(a^2 \log(1/a))$ lattice artifact in the vacuum polarization with a prefactor that we calculate. It is then applied to non-perturbative lattice QCD data with two dynamical flavors of quarks. Our lattice spacings range down to 0.049 fm for the vacuum simulations and down to 0.033 fm for the simulations performed at a temperature of 250 MeV.

preprint2020arXiv

Frequency-splitting estimators for single-propagator traces

In these proceedings we address the computation of quark-line disconnected diagrams in lattice QCD. The evaluation of these diagrams is required for many phenomenologically interesting observables, but suffers from large statistical errors due to the vacuum and random-noise contributions to their variances. Motivated by a theoretical analysis of the variances, we introduce a new family of stochastic estimators of single-propagator traces built upon a frequency splitting combined with a hopping expansion of the quark propagator, and test their efficiency in two-flavour QCD with pions as light as 190 MeV. The use of these estimators reduces the cost of the computation by one to two orders of magnitude over standard estimators depending on the fermion bilinear. As a concrete application, we show the impact of these findings on the computation of the hadronic vacuum polarization contribution to the muon anomalous magnetic moment.

preprint2016arXiv

A variational method for spectral functions

The Generalized Eigenvalue Problem (GEVP) has been used extensively in the past in order to reliably extract energy levels from time-dependent Euclidean correlators calculated in Lattice QCD. We propose a formulation of the GEVP in frequency space. Our approach consists of applying the model-independent Backus-Gilbert method to a set of Euclidean two-point functions with common quantum numbers. A GEVP analysis in frequency space is then applied to a matrix of estimators that allows us, among other things, to obtain particular linear combinations of the initial set of operators that optimally overlap to different local regions in frequency. We apply this method to lattice data from NRQCD. This approach can be interesting both for vacuum physics as well as for finite-temperature problems.

preprint2016arXiv

The End of a Myth: Distributed Transactions Can Scale

The common wisdom is that distributed transactions do not scale. But what if distributed transactions could be made scalable using the next generation of networks and a redesign of distributed databases? There would be no need for developers anymore to worry about co-partitioning schemes to achieve decent performance. Application development would become easier as data placement would no longer determine how scalable an application is. Hardware provisioning would be simplified as the system administrator can expect a linear scale-out when adding more machines rather than some complex sub-linear function, which is highly application specific. In this paper, we present the design of our novel scalable database system NAM-DB and show that distributed transactions with the very common Snapshot Isolation guarantee can indeed scale using the next generation of RDMA-enabled network technology without any inherent bottlenecks. Our experiments with the TPC-C benchmark show that our system scales linearly to over 6.5 million new-order (14.5 million total) distributed transactions per second on 56 machines.

preprint2015arXiv

Non-perturbative improvement of the vector current in Wilson lattice QCD

Many observables of interest in lattice QCD are extracted from correlation functions involving the vector current. If Wilson fermions are used, it is therefore of practical importance that, besides the action, the current be O($a$) improved in order to remove the leading discretization errors from the observables. Here we introduce and apply a new method to determine the improvement coefficient for the two most widely used discretizations of the current.

preprint2015arXiv

Nucleon electromagnetic form factors and axial charge from CLS $N_\mathrm{f}=2+1$ ensembles

We present preliminary results on the electromagnetic form factors and axial charge of the nucleon from ensembles generated by the CLS effort with $N_\mathrm{f}=2+1$ flavours of non-perturbatively $\mathrm{O}(a)$-improved Wilson fermions and open temporal boundary conditions. Systematic effects due to excited-state contamination are accounted for using both two-state fits and the method of summed operator insertions. This exploratory analysis demonstrates the viability of obtaining precision baryon observables with $N_\mathrm{f}=2+1$ flavours of Wilson fermions on fine lattices, aiming towards controlled chiral and continuum limits in the future.

preprint2015arXiv

The Influence of Malloc Placement on TSX Hardware Transactional Memory

The hardware transactional memory (HTM) implementation in Intel's i7-4770 "Haswell" processor tracks the transactional read-set in the L1 (level-1), L2 (level-2) and L3 (level-3) caches and the write-set in the L1 cache. Displacement or eviction of read-set entries from the cache hierarchy or write-set entries from the L1 results in abort. We show that the placement policies of dynamic storage allocators -- such as those found in common "malloc" implementations -- can influence the L1 conflict miss rate in the L1. Conflict misses -- sometimes called mapping misses -- arise because of less than ideal associativity and represent imbalanced distribution of active memory blocks over the set of available L1 indices. Under transactional execution conflict misses may manifest as aborts, representing wasted or futile effort instead of a simple stall as would occur in normal execution mode. Furthermore, when HTM is used for transactional lock elision (TLE), persistent aborts arising from conflict misses can force the offending thread through the so-called "slow path". The slow path is undesirable as the thread must acquire the lock and run the critical section in normal execution mode, precluding the concurrent execution of threads in the "fast path" that monitor that same lock and run their critical sections in transactional mode. For a given lock, multiple threads can concurrently use the transactional fast path, but at most one thread can use the non-transactional slow path at any given time. Threads in the slow path preclude safe concurrent fast path execution. Aborts rising from placement policies and L1 index imbalance can thus result in loss of concurrency and reduced aggregate throughput.

preprint2014arXiv

The bottomonium spectrum at finite temperature from $N_f=2+1$ lattice QCD

We present results on the bottomonium spectrum at temperatures above and below the deconfinement crossover temperature, $T_c$, from dynamical lattice QCD simulations. The heavy quark is treated with a non-relativistic effective field theory on the lattice and serves as a probe of the hot medium. Ensembles with a finer spatial lattice spacing and a greater range of temperatures below $T_c$ than those previously employed by this collaboration are used. In addition, there are $N_f=2+1$ flavours of Wilson clover quark in the sea with $M_π\approx400$ MeV and we perform a more careful tuning of the bottom quark mass in this work. We calculate the spectral functions of S and P wave bottomonium states using the maximum entropy method and confirm earlier findings on the survival of the ground state S wave states up to at least $2T_c$ and the immediate dissociation of the P wave states above $T_c$.

preprint2013arXiv

Bottomonium spectrum at finite temperature

We investigate the modification of S and P wave states in the bottomonium spectrum above and below the deconfinement crossover temperature through their spectral functions obtained from the maximum entropy method. Anisotropic ensembles with $N_f=2+1$ Wilson clover fermions with tadpole-improvement are used while the bottom quark is treated with an improved non-relativistic action. We observe significant modifications of the P wave $χ_{b1}$ ground state directly above the crossover temperature, $T_c$, while the S wave $Υ$ ground state remains relatively unchanged up to temperatures of almost $2T_c$. This evidence supports earlier conclusions from our FASTSUM collaboration of the immediate dissociation of the P wave states above $T_c$ and the survival of the S wave ground states up to $2T_c$.

preprint2013arXiv

Quark-gluon plasma phenomenology from the lattice

The FASTSUM Collaboration has calculated several quantities relevant for QCD studies at non-zero temperature using the lattice technique. We report here our results for the (i) interquark potential in charmonium; (ii) bottomonium spectral functions; and (iii) electrical conductivity. All results were obtained with 2+1 flavours of dynamical fermions on an anisotropic lattice which allows greater resolution in the temporal direction.

preprint2010arXiv

Mechanism of ferroelectric instabilities in non d^0 perovskites: LaCrO_3 versus CaMnO_3

The incompatibility of partial d occupation on the perovskite B-site with the standard charge transfer mechanism for ferroelectricity has been a central paradigm in multiferroics research. Nevertheless, it was recently shown by density functional theory calculations that CaMnO_3 exhibits a polar instability that even dominates over the octahedral tilting for slightly enlarged unit cell volume. Here, we present similar calculations for LaCrO_3, which has the same d^3 B-site electron configuration as CaMnO_3. We find that LaCrO_3 exhibits a very similar, albeit much weaker, polar instability as CaMnO_3. In addition, while the Born effective charge (BEC) of the Mn^{4+} cation in CaMnO_3 is highly anomalous, the BEC of Cr^{3+} in LaCrO_3 is only slightly enhanced. By decomposing the BECs into contributions of individual Wannier functions we show that the ferroelectric instabilities in both systems can be understood in terms of charge transfer between TM d and O p states, analogously to the standard d^0 perovskite ferroelectrics.