The 33rd International Symposium on Lattice Field Theory (Lattice 2015)

Asia/Tokyo
Kobe International Conference Center

Kobe International Conference Center

6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
Yasumichi Aoki (Nagoya University), Yoshinobu Kuramashi (University of Tsukuba)
Description

The 33rd International Symposium on Lattice Field Theory was held at the
Kobe International Conference Center, Kobe, Japan from Tuesday, July 14
to Saturday, July 18, 2015.

The program includes plenary talks, parallel talks and a poster session on
the following topics:

  • Algorithms and Machines
  • Applications Beyond QCD
  • Chiral Symmetry
  • Hadron Spectroscopy and Interactions
  • Hadron Structure
  • Nonzero Temperature and Density
  • Physics Beyond the Standard Model
  • Standard Model Parameters and Renormalization
  • Theoretical Developments
  • Vacuum Structure and Confinement
  • Weak Decays and Matrix Elements

For further details, please refer to http://www.r-ccs.riken.jp/sympo/lattice2015/.

    • Plenary Session: Conference Opening BF & 1F Main Hall

      BF & 1F Main Hall

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Philippe de Forcrand
      • 09:00
        Conference Opening
      • 1
        Recent progress on chiral symmetry breaking in QCD
        I review recent progress achieved on the lattice in the understanding of chiral symmetry breaking in QCD. Emphasis is given to recent computations of the spectral density of the Dirac operator and of the topological susceptibility.
        Speaker: Prof. Leonardo Giusti (University of Milano Bicocca)
        Slides
      • 2
        Determining the order of chiral phase transition in QCD from conformal bootstrap
        There has been a longstanding debate if the chiral phase transition in two-flavor massless QCD is the first order or the second order. The previous arguments based on epsilon expansions, large N expansions, functional renormalization group, and Monte-Carlo simulations had been all inconclusive with shortcomings. If it were the second order phase transition, there should exist a corresponding three-dimensional conformal field theory which describes the critical phenomenon. The recent development in conformal bootstrap enables us to directly study the (non-)existence of conformal fixed points in a non-perturbative manner. In this talk, I review the conformal bootstrap method and its application to this problem. Our conclusion is that the corresponding conformal fixed point should exist and the phase transition will be the second order if the U(1) chiral anomaly is effectively restored. This means that the original 1-loop prediction by Pisarski and Wilczek would be incorrect. We further provide the most precise prediction of the critical exponent there. We believe future numerical simulations will confirm our prediction.
        Speaker: Dr YU NAKAYAMA (Caltech)
        Slides
    • 10:30
      Coffee Break
    • Plenary Session BF & 1F Main Hall

      BF & 1F Main Hall

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Maarten Golterman
      • 3
        Numerical evaluation of QED contribution to lepton g-2
        The anomalous magnetic moment (g-2) of the electron has played the central role in testing the validity of quantum electrodynamics (QED) as well as the standard model of particle physics. The test has been further improved, which was made possible by the complete evaluation of the tenth-order term in the perturbation theory of QED, together with the latest measurement of the electron g-2 by the Harvard group that has reached the precision of 0.24 ppb. In this talk we will present the numerical approach to the evaluation of QED contribution to lepton g-2 up to the recent developments.
        Speaker: Dr Tatsumi Aoyama (KMI, Nagoya University)
        Slides
      • 4
        Lattice QCD moments - $g-2$ and NEDM -
        There has been much progress in studies of the physics of low energy, precision observables in the flavor neutral sector using Lattice QCD, especially on the lepton anomalous magnetic moments ($g-2$) and nucleon electric dipole moments. The current status of these studies is reviewed, and remaining challenges and new ideas are discussed.
        Speaker: Taku Izubuchi (RIKEN BNL Research Center)
        Slides
    • 12:30
      Lunch Break
    • Algorithms and Machines 406

      406

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Nigel Cundy
      • 5
        Algorithmic improvements for weak coupling simulations of domain wall fermions
        We discuss algorithmic improvements being used in the evolution of new 2+1+1 flavor lattices by the RBC/UKQCD collaboration. A technique similar to Hasenbusch mass splitting and suggested previously by Brower, Neff, and Orginos allows the light quark action to be split into multiple parts with different values of L_s, the size of the fifth dimension. This allows us to reduce chiral symmetry breaking by making L_s large while most inversions are actually done at a cheaper, smaller value of L_s.
        Speaker: Mr Greg McGlynn (Columbia University)
        Slides
      • 6
        zMobius and other recent developments on Domain Wall Fermions
        Recent advances in hardware, combined with various deflation techniques has made it crucial to optimize not only for solver performances, but also for memory and disk footprint to keep the overall efficiency high. Here we report on advances in Domain Wall Fermion formalism(DWF), such as zMobius which achieves a good approximation to the scaled Shamir DWF with smaller extent in 5th dimension, and delayed deflation technique which eliminates the need for the source projection by using predetermined approximation for the inverse.
        Speaker: Chulwoo Jung (Brookhaven National Laboratory)
        Slides
      • 7
        Fermionic twisted boundary conditions with reweighting method
        Imposing twisted boundary conditions on the fermion fields is a procedure extensively used when evaluating, for example, form factors on the lattice. Twisting is usually performed for one flavour and only in the valence, and this causes a breaking of unitarity. In this work we explore the possibility to restore the unitarity through the reweighting method. We evaluate the change in the fermion determinant with different boundary conditions and include that in the expectation values, avoiding in this way the cost of generating new configurations for each choice of the twisting angle, theta. As expected the effect of the reweighting is negligible in the case of large volumes but it is important when the volumes are small and the twisting angles are large. In particular we find a measurable effect for the plaquette and the pion correlation function in the case of theta=pi/2 in a volume 16x8^3, and we found a systematic upward shift in the pion dispersion relation.
        Speaker: Mr Andrea Bussone (CP3-Origins, University of Southern Denmark)
        Slides
      • 8
        Optimizing the domain wall fermion Dirac operator using the R-Stream source-to-source compiler
        The application of the Dirac operator on a spinor field, the Dslash operation, is the most computation-intensive part of the lattice QCD simulations. It is often the key kernel to optimize to achieve maximum performance on various platforms. Here we report on a project to optimize the domain wall fermion Dirac operator in Columbia Physics System (CPS) using the R-Stream source-to-source compiler. Our initial target platform is the PC clusters. We discuss the optimization strategies involved before and after the automatic code generation with R-Stream and present some preliminary benchmark results.
        Speaker: Dr Meifeng Lin (Brookhaven National Laboratory)
      • 9
        A next generation C++ library for data parallel QCD
        We discuss progress and performance from the development of a data parallel interface for cartesian fields of tensor data types. The library is suitable for the simulation of QCD including multiple grids. Wherever appropriate the interface is quite similar to that of QDP++, but the library makes use of C++11 features to reduce the volume of code, compared to QDP++ while bringing greater generality and greater performance. While the library simultaneously targets MPI, OpenMP and SIMD parallelism, the SIMD optimisation is notably flexible. Performance is substantially improved through the use of a data layout transformation and benchmark results are presented for simple Lattice QCD operations and the Wilson operator under AVX, AVX2 and AVX512 instruction set targets.
        Speaker: Prof. Peter Boyle (University of Edinburgh)
        Slides
    • Hadron Spectroscopy and Interactions 402

      402

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Naruhito Ishizuka
      • 10
        Pion-pion interaction from N_f = 2+1 simulations using the stochastic LapH method
        We present results on pion-pion scattering in the isovector channel from N_f = 2+1 CLS simulations with open temporal boundary conditions and m_\pi = 280 MeV. The required correlation functions are computed using the stochastic LapH method that facilitates large-volume calculations in order to maintain m_\pi L > 4. Mapping out the \rho resonance structure is a necessary step towards the extraction of the timelike pion form factor.
        Speaker: Ben Hörz (School of Mathematics, Trinity College Dublin)
        Slides
      • 11
        Large volume calculation of pion-pion scattering phase shifts with the stochastic LapH method on an $N_f=2+1$ anisotropic clover lattice
        Results for pion-pion scattering phase shifts using the stochastic LapH method on a single $N_f = 2+1$ anisotropic clover lattice with $m_{\pi}\equiv 240\mathrm{MeV}$ and $L\equiv 3.8\mathrm{fm}$ will be presented. In particular, for $I=1$ this large volume enables good resolution of the $\rho$ resonance shape. Prospects for the other isospin channels will also be discussed.
        Speaker: Dr John Bulava (Trinity College Dublin)
        Slides
      • 12
        K-pi scattering lengths from physical point ensembles
        I will present preliminary results on K-pi scattering lengths in both I=3/2 and I=1/2 channels obtained by the RBC-UKQCD collaboration. The continuum results are obtained using two of our Moebius domain wall fermion ensembles - 48^3 x 96 and 64^3 x 128 - with physical quark masses. I will demonstrate how the large finite time extent effects which are present with such light quark masses can be controlled.
        Speaker: Mr Tadeusz Janowski (University of Southampton)
        Slides
      • 13
        K-pi and pi-pi scattering close to the physical pion mass
        Results will be presented for pi-pi and K-pi scattering using lattice configurations on which the pion mass is 150 MeV. In this study the focus is on the kinematic region around the rho (for pi-pi scattering) and K* (for K-pi) resonances. Different lattice volumes and moving frames are used to obtain the spectrum for a range of centre of mass energies, which allows the energy dependence of the scattering phase shift to be extracted.
        Speaker: Dr Gordon Donald (University of Regensburg)
        Slides
      • 14
        The mass and leptonic decay constant of rho meson at the physical point
        We use overlap valence quarks on the $L^3\times T=48^3\times 96$ $N_f=2+1$ domain-wall fermion configurations generated by the RBC\&UKQCD Collaboration. The light sea quark mass is at the physical point and the spatial extension is 5.5 fm. The $\rho$ mass is precisely determined to be $m_\rho=778(9)$ MeV at the physical pion mass. We also propose a strategy to reduce the noise of the wall-wall correlation functions of hadrons, through which the leptonic decay constant of $\rho$, $f_\rho$, is determined to be 209(4) MeV at the physical $m_\pi$. This value is in excellent agreement with the experimental value $f_\rho=208$ MeV derived from the process $\tau\rightarrow \rho \nu_\tau$.
        Speaker: Prof. ying chen (institute of high energy physics, chinese academy of sciences, china)
        Slides
      • 15
        Resonance Parameters for the rho-meson from Lattice QCD
        We present a lattice calculation of the mass and width for the rho meson by determining the energy dependent phaseshift for pi-pi scattering in isospin-1, spin-1 channel in the elastic region. We compute the energy spectrum from spatially asymmetric boxes. We use Nf=2 nHYP fermions and generate two sets of ensembles with pion masses of 300MeV and 220MeV. To determine the energy spectrum we perform a variational analysis with interpolating fields including several q-\bar{q} and pi-pi interpolating fields.
        Speaker: Mr Dehua Guo (The George Washington University)
        Slides
    • Hadron Structure 403

      403

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Georg von Hippel
      • 16
        Hadronic Light by Light Contributions to the Muon Anomalous Magnetic Moment With Near Physical Pions
        The current measurement of muonic g-2 disagrees with the theoretical calculation by about 3 standard deviations. Hadronic vacuum polarization (HVP) and hadronic light by light (HLbL) are the two types of processes that contribute most to the theoretical uncertainty. The current value for HLbL is still given by a model . In this talk, I will describe results from a first-principles lattice calculation with a 171MeV pion in a box of 4.6 fm extent. Our current numerical strategies, including noise reduction techniques, evaluating the HLbL amplitude at zero external momentum transfer, and important remaining challenges, in particular those associated with finite volume effects, will be discussed.
        Speaker: Mr Luchang Jin (Columbia University)
        Slides
      • 17
        On calculating disconnected-type hadronic light-by-light scattering diagrams from lattice QCD
        The first-principles calculation of the hadronic light-by-light scattering (HLbL) contribution is the most critical issue in the theory of the muon $g-2$ to confront the standard model prediction with the more accurate experimental result to be provided in the next five years. The feasibility of the lattice QCD computation of the HLbL contribution has been demonstrated for the “connected diagram” where four electromagnetic (EM) vertices appear on a single quark loop. In this talk, we focus on “disconnected diagrams” where the four EM vertices are distributed over more than one quark loop. We note that any method to calculate the disconnected diagrams requires subtraction of unwanted contributions which may affect its feasibility. We present an example of the method with a concrete realization for such a subtraction.
        Speaker: Mr Masashi Hayakawa (Nagoya University)
        Slides
      • 18
        Improving the volume-dependence of lattice QCD+QED simulations
        We explore methods to improve the volume-dependence of lattice QCD+QED simulations. Applications to QED mass-splittings, QED corrections to decay constants, and the computation of the hadronic light-by-light contribution to the muon anomalous magnetic moment are discussed.
        Speaker: Christoph Lehner (BNL)
        Slides
      • 19
        The strange contribution to the anomalous magnetic moment of the muon with physical quark masses using Möbius domain wall fermions
        We present preliminary results for the strange leading-order hadronic contribution to the anomalous magnetic moment of the muon using RBC/UKQCD physical point domain wall fermions ensembles. We discuss various analysis strategies in order to constrain the systematic uncertainty in the final result.
        Speaker: Mr Matt Spraggs (University of Southampton)
        Slides
      • 20
        Finite volume hadronic vacuum polarisation at arbitrary momenta
        The muon anomalous magnetic moment constitute one of the most important tension between experimental measurements and the Standard Model predictions. The dominant theoretical uncertainties is coming from hadronic contributions, therefore it is crucial to aim at determining precisely and reliably this quantity from lattice QCD. Most of the contributions to the leading hadronic contribution comes from the low q^2 sector of the hadronic vacuum polarization (HVP). Because of momentum quantization in finite volume, it is hard to compute the HVP at low momentum and generally one has to use models to describe this region which can introduce significant systematic errors. Here we propose a model-independent method to reconstruct the HVP at continuous momenta. Using sampling theory, we show that this reconstruction is exact up to corrections that decay exponentially fast with the pion mass times the lattice extent. We conclude by presenting preliminary numerical results and we compare them to model-dependent approaches.
        Speaker: Dr Antonin Portelli (University of Southampton)
        Slides
      • 21
        Lattice calculation of the HVP contribution to the anomalous magnetic moment of muon
        We report our (HPQCD) progress with the lattice QCD calculations of both the connected and disconnected hadronic Vacuum Polarisation (HVP) contributions to the anomalous magnetic moment of the muon (g-2). Our method uses Padé approximants to reconstruct the Adler function from its derivatives at q^2=0. These are obtained simply and accurately from time-moments of the vector current-current correlator at zero spatial momentum. We give a full flavor estimate of the total HVP contribution calculated on MILC Collaboration's n_f = 2+1+1 HISQ ensembles at multiple values of the lattice spacing, multiple volumes and multiple light sea quark masses (including physical pion mass configurations) and corrected for finite volume effects.
        Speaker: Ms bipasha Chakraborty (University of Glasgow)
        Slides
    • Nonzero Temperature and Density 401

      401

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Owe Philipsen
      • 22
        The curvature of the chiral phase transition line at small values of the quark chemical potentials
        The QCD chiral phase transition for small values of baryon chemical potentials is expected to be second-order and the dependence of its transition temperature on the chemical potentials is thus controlled by the dependence of $O(N)$ scaling functions on the reduced temperature variables. Though the transition is merely a crossover at physical quark mass, the dependence of the chiral phase transition temperature is nevertheless of interest as it could be related to the freeze-out curve that is obtained from an analysis of hadron multiplicities in heavy-ion collisions. In this talk we will report on our ongoing work on the leading-order dependence of the chiral phase transition temperature as a function of the light and strange quark chemical potentials i.e. the curvature of the chiral phase transition temperature with respect to the light and the strange chemical potentials. We have performed simulations using the Highly Improved Staggered Quark (HISQ) action on lattices with temporal extent $N_\tau=6$ with two degenerate light quarks and a strange quark. The strange quark mass was chosen to be fixed to its physical value and the light quark mass was adjusted in order to have five values of the pseudo-Goldstone pion mass ranging from 160 MeV to about 80 MeV in the continuum limit. The comparison of the resulting curvatures from the current study with those from previous investigations will also be discussed. This study is on behalf of the BNL-Bielefeld-CCNU collaboration.
        Speaker: Dr Prasad Hegde (Central China Normal University)
        Slides
      • 23
        The curvature of the crossover line in the (T, mu)-phase diagram of QCD
        An efficient way to study the QCD phase diagram at small finite density is to extrapolate thermodynamical observables from imaginary chemical potential. The phase diagram features a crossover line starting from the transition temperature already determined at zero chemical potential. In this talk we focus on the curvature of this line at mu = 0. We present the extrapolation of the crossover temperature based on three observables at several lattice spacings. The simulations were performed at zero and at moderate values of the imaginary chemical potential, always in the strangeness neutral point. We used the Symanzik-improved gauge action with four times stout smeared staggered fermions.
        Speaker: Ms Jana Guenther (Wuppertal University)
        Slides
      • 24
        Curvature of the pseudocritical line in (2+1)-flavor QCD with HISQ fermions
        We study QCD with (2 + 1) HISQ fermions at nonzero temperature and nonzero imaginary baryon chemical potential. Monte Carlo simulations are performed using the MILC code along the line of constant physics with a light to strange mass ratio of ml/ms=1/20 on lattices up to 48^3x12 to check for finite cutoff effects. We determine the curvature of the pseudocritical line extrapolated to the continuum limit.
        Speaker: Dr Leonardo Cosmai (INFN)
        Slides
      • 25
        Curvature of the QCD chiral peusdocritical line from analytic continuation
        We present our latest results on the determination of the curvature of the pseudo-critical line of the QCD phase diagram at the physical point, using the method of analytic continuation from an imaginary chemical potential. We also assess the impact of including a non-zero strange quark chemical potential. Our results are obtained with stout improved staggered fermions and the tree level Symanzik improvement for the gauge action.
        Speaker: Michele Mesiti (Universita' di Pisa)
        Slides
      • 26
        Phase structure of Nf=3 QCD at finite temperature and density by Wilson-Clover fermions
        We investigate the phase structure of finite temperature and density QCD with three-flavor by using Wilson-Clover fermions. We focus on locating the critical end point that characterizes the phase structure and extracting the curvature of critical line on the plane of quark chemical potential and pseudo-scalar meson mass.
        Speaker: Dr Shinji TAKEDA (Kanazawa university)
        Slides
      • 27
        Many flavor approach to study the critical point in finite density QCD
        We discuss the QCD critical point at finite density through the study of many-flavor QCD, in which two light flavors and Nf massive flavors exist. Performing simulations of QCD with 2 flavors of improved Wilson fermions, we calculate probability distribution functions in many-flavor QCD at finite temperature and density, where the reweighting technique is used to add the dynamical effect of massive flavors and the chemical potential. From the shape of the distribution functions, we determine the critical surface separating the first order transition and crossover regions in the space spanned by the light and massive quark masses and the chemical potentials. It is found that the critical massive quark mass becomes larger as the chemical potential increases in (2+Nf)-flavor QCD. The indication to the (2+1)-flavor QCD is then discussed.
        Speaker: Mr Ryo IWAMI (Niigata University)
        Slides
    • Physics Beyond the Standard Model 405

      405

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: C.-J. David Lin
      • 28
        Walking and conformal dynamics in many flavor QCD
        In the search for a realistic walking technicolor model, QCD with many flavors is an attractive candidate. From the series of studies by the LatKMI collaboration, we present updated results of the scaling properties of various hadron spectra, including the (pseudo)scalar, vector, and baryon channels, for Nf=8 QCD with the HISQ action. By comparing these with Nf=12 QCD, which has properties consistent with conformality, possible signals of walking dynamics will be discussed.
        Speaker: Dr Hiroshi Ohki (RIKEN)
      • 29
        Topological observables in many-flavour QCD
        SU(3) gauge theory with eight massless flavours is believed to be walking, while the corresponding twelve- and four-flavour appear IR-conformal and confining respectively. Looking at the simulations performed by the LatKMI collaboration of these theories, we use the topological susceptibility as an additional probe of the IR dynamics. By drawing a comparison with SU(3) pure gauge theory, we see a dynamical quenching effect emerge at larger number of flavours, which is suggestive of emerging near-conformal and conformal behaviour.
        Speaker: Dr Ed Bennett (Swansea University)
        Slides
      • 30
        SU(3) gauge theory with four degenerate fundamental fermions on the lattice
        LatKMI Collaboration has been studying SU(3) gauge theories with 4, 8, 12 and 16 fundamental fermions by using HISQ action on the lattice. This talk is dedicated to show results of Nf=4 simulations. We first show the m_rho/ f_pi ratio, from which we see Nf=4 SU(3) gauge theory is quite similar to the actual QCD. Then, we’ll show preliminary results of flavor-singlet scalar bound state and give some discussions.
        Speaker: Dr Masafumi Kurachi (KEK)
      • 31
        Studying near conformal behavior with four light flavors and eight flavors of variable mass.
        After a brief review of the motivation for studying an SU(3) theory with four light flavors and eight flavors of variable mass I will discuss the methodologies we have used to investigate the properties of the system. I will then illustrate some recent results. More detailed results will be given in another talk. We will also illustrate some results of a study of the SU(3) theory with eight staggered flavors carried out by the LSD collaboration.
        Speaker: Prof. Claudio Rebbi (Boston University)
        Slides
      • 32
        Probing near conformal dynamics with 4+8 and 8 flavors: running coupling and the spectrum
        We investigate an SU(3) gauge theory with four light flavors and eight heavy flavors of variable mass as a model of near conformal dynamics. By continuously varying the mass of the eight heavy flavors we can tune between chirally broken and conformal behavior. We report on recent progress with probing this model deeper in the light flavor chiral regime, including updated measurements of a light 0++ isosinglet meson. This will be presented alongside a study of an SU(3) theory with eight light flavors by the LSD collaboration.
        Speaker: Evan Weinberg (Boston University Department of Physics)
        Slides
      • 33
        IR fixed points and conformal window in $SU(3)$ gauge Theories
        We propose a novel RG method to specify the location of the IR fixed point in lattice gauge theories and apply it to the $SU(3)$ gauge theories with $N_f$ fundamental fermions. It is based on the scaling behavior of the propagator through the RG analysis with a finite IR cut-off, which we cannot remove in the conformal field theories in sharp contrast with the confining theories. The method also enables us to estimate the anomalous mass dimension in the continuum limit at the IR fixed point. We perform the program for $N_f=16, 12, 8 $ and $N_f=7$ and indeed identify the location of the IR fixed points in all cases. Our results are consistent with that the conformal window is 7 \le N_f \le 16.
        Speaker: Dr YU NAKAYAMA (Caltech)
        Slides
    • Standard Model Parameters and Renormalization 404

      404

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Francesco Knechtli
      • 34
        A status update on the ALPHA collaboration's project to determine the Lambda-parameter in 3-flavour QCD
        The ALPHA collaboration has set itself the goal to determine alpha_s(m_Z) with a total error substantially below the percent level. A further step towards this goal can be taken by combining results from the recent simulations of 2+1 flavour QCD by the CLS initiative with a number of tools developed over the years: renormalized couplings in finite volume schemes, recursive finite size techniques, two-loop renormalized perturbation theory and the (improved) gradient flow on the lattice. I will sketch the strategy, which involves both the standard SF coupling in the high energy regime and a gradient flow coupling at low energies. This implies the need for matching both schemes at an intermediate switching scale. Our results for the scale evolution at high energies allow for a well-controlled continuum limit and thus yield the Nf=3 Lambda-parameter in units of this switching scale. Simulations at lower energy scales are underway. I will present first results for the scale evolution of a gradient flow coupling with SF boundary conditions.
        Speaker: Stefan Sint (Trinity College Dublin)
        Slides
      • 35
        Prospects and status of quark mass renormalization in three-flavour QCD
        We present the current status of a revised strategy to compute the running of renormalized quark masses in QCD with three flavours of massless O(a) improved Wilson quarks. The strategy employed uses the standard finite-size scaling method in the Schroedinger functional and accommodates for the non-perturbative scheme-switch which becomes necessary at intermediate renormalized couplings as discussed in [1411.7648].
        Speaker: Dr Patrick Fritzsch (Instituto de Física Teoríca UAM/CSIC)
        Slides
      • 36
        Non-perturbative renormalization of tensor bilinears in Schrodinger Functional schemes
        We present preliminary results for the non-perturbative renormalization of the tensor current in QCD for $N_f=0,2$, as well as an ongoing extension to $N_f=2+1$. The renormalization group running is computed in the continuum limit for a wide range of scales, using various Schrodinger Functional (SF) schemes and finite volume recursive techniques. In all these schemes, we have also computed the matching factor to MSbar and RI-MOM, which allows for a NLO matching to the RGI operator at high energies. An example application of the results is the renormalization and matching of the tensor currents entering the effective Hamiltonian for $b \to s$ transitions.
        Speaker: Mr David Preti (Instituto de Física Teoríca UAM/CSIC)
        Slides
      • 37
        Scale determination for the CLS 2+1 ensembles
        During the last two years, the CLS effort has generated ensembles with 2+1 dynamical flavors of non-perturbatively improved Wilson fermions at lattice spacings between 0.05fm and 0.09fm. Most of these ensembles are along lines of constant sum of bare sea quark masses tr(M). We give an overview of these simulations and describe our scale setting procedure using the pseudoscalar decay constants. Since this action is relatively new, particular attention will be given to the size of the discretization effects and the impact of the mistuning of tr(M).
        Speaker: Stefan Schaefer (NIC, DESY)
        Slides
      • 38
        Perturbative renormalization of ∆S = 2 four-fermion operators with the chirally rotated Schrödinger functional
        The chirally rotated Schrödinger functional (χSF) is a renormalization scheme which renders the mechanism of automatic O(a) improvement compatible with the Schrödinger functional (SF) formulation. Here we define a family of renormalization schemes based on the χSF for a complete basis of ∆S = 2 parity-odd four-fermion operators. We compute the scale-dependant renormalization constants of such operators to one-loop in perturbation theory and obtain their NLO anomalous dimensions. After this is done, we compute the cutoff effects in the corresponding step-scaling functions at one-loop. Due to automatic O(a) improvement, once the χSF action is renormalized and O(a) improved, renormalization constants are affected directly by O(a^{2}) effects without the need of operator improvement.
        Speaker: Dr Pol Vilaseca Mainar (Instituto Nazionale de Fisica Nucleare)
        Slides
      • 39
        The chirally rotated Schrödinger functional at work
        The chirally rotated Schrödinger functional has proven to be a powerful tool in addressing non-perturbative renormalization problems in lattice QCD with Wilson-type fermions. In this contribution we consider two novel applications of the method. Firstly, we investigate the renormalization of a complete basis of ΔS=2 four-quark operators relevant for searches of Beyond the Standard Model physics. Preliminary results are presented for the theory with Nf=2 dynamical flavours. Secondly, we discuss the renormalization of several quark-bilinears in the Nf=2+1 theory.
        Speaker: Dr Mattia Dalla Brida (DESY - Zeuthen)
        Slides
    • 16:00
      Coffee Break
    • Hadron Spectroscopy and Interactions 402

      402

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Tomomi Ishikawa
      • 40
        Decay constants and spectroscopy of mesons in lattice QCD using domain-wall fermions
        We report results of masses and decay constants of light and charmed psuedo-scalar mesons using lattice QCD with Moebius domain-wall fermions. Using this formulation we are able to compute psuedoscalar decay constants through the pseudo-scalar density operator as well as with the axial-vector current. Results are shown from several lattice spacings and pion masses between 240 MeV and 500 MeV. We present an analysis of these results at different quark masses to show the chiral properties of the light mesons masses and decay constants.
        Speaker: Dr Brendan Fahy (KEK)
        Slides
      • 41
        2+1 flavor QCD simulation near the physical point on a $96^4$ lattice
        We present our preliminary results on 2+1 flavor QCD with the stout smeared O(a)-improved Wilson quark action and the Iwasaki gauge action. Simulations are carried out at a lattice spacing 0.085fm on a $96^4$ lattice with the quark masses near the physical point. The reweighing technique and ChPT are employed for the extrapolation to the physical point. We will show the preliminary results of the light quark masses, the decay constants and the light hadron spectrum.
        Speaker: Naoya Ukita (University of Tsukuba)
        Slides
      • 42
        Latest Results from RQCD using 2+1f CLS Simulations with Open Boundaries
        We will report on latest results of the RQCD group obtained from 2+1f gauge field configurations with non-perturbatively improved Wilson action and open boundaries generated within the CLS effort. We have extended our simulations to include two chiral trajectories, one keeping the average quark mass fixed and an additional one where the strange quark mass is kept at its physical value. For the latter, we present details of our tuning strategy which enables us to fix the strange quark mass at the percent level (or possibly below). We will focus on spectroscopy, scale setting, and improvement coefficients.
        Speaker: Dr Wolfgang Soeldner (Regensburg University)
        Slides
      • 43
        Extracting the eta-prime meson mass from gluonic correlators in lattice QCD
        Calculation of the eta-prime mass is a notoriously difficult problem, as it requires evaluation of the disconnected diagram which is costly and noisy. In this work, we use a gluonic operator to extract the eta-prime state after smearing the link variables through the Wilson flow. With this choice, one can avoid a large cancellation of pion contribution between the connected and disconnected diagrams. We obtain the eta-prime meson mass on lattices with three different lattice spacings and two physical volumes, which allow us to estimate its continuum and large volume limits.
        Speaker: Dr Hidenori Fukaya (Osaka University)
        Slides
      • 44
        Photon mass term as an IR regularization for QCD+QED on the lattice
        The commonly adopted approach for including QED in lattice QCD simulations introduces power-law finite volume corrections to physical quantities. These effects, which are due to the long-range nature of the electromagnetic interaction, must be removed by performing simulations at multiple lattice volumes, followed by an extrapolation to the infinite volume limit. In this work, we explore the advantages and disadvantages of introducing a photon mass term as an alternative means for gaining control over the finite volume effects associated with the inclusion of QED. We present exploratory findings for hadron mass shifts due to electromagnetic interactions (i.e., for the proton, neutron, charged and neutral kaon) and corresponding mass splittings, and compare them with standard QCD+QED calculations. Preliminary results are reported for numerical studies of three flavor electroquenched QCD using ensembles corresponding to 800 MeV pions and three lattice volumes ranging from 3.4 fm to 6.7 fm.
        Speaker: Dr Michael Endres (Massachusetts Institute of Technology)
        Slides
      • 45
        Dashen's theorem and electromagnetic contributions to pseudoscalar meson masses
        We present results on the pseudoscalar meson masses from a fully dynamical simulation of QCD plus QED. We concentrate particularly on the question of how to divide the measured masses into a QED contribution and a QCD contribution. This decomposition is not unique, it depends on the renormalisation scheme and scale. We suggest a renormalisation scheme in which Dashen's theorem holds, so that the electromagnetic self-energies of the neutral mesons are zero, and discuss how the self-energies change when we transform to a scheme such as MS-bar, in which Dashen's theorem is violated.
        Speaker: Dr Paul Rakow (University of Liverpool)
        Slides
    • Hadron Structure 403

      403

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Masashi Hayakawa
      • 46
        Direct calculation of hadronic light-by-light scattering
        Calculations of hadronic light-by-light scattering amplitudes via lattice QCD evaluation of Euclidean four-point functions of vector currents will be reported. These initial results include only the fully quark-connected contribution. Particular attention will be given to the case of forward scattering, which can be related via dispersion relations to the gamma* gamma* -> hadrons cross section, and thus allows lattice data to be compared with phenomenology. A strategy for computing the hadronic light-by-light contribution to the muon anomalous magnetic moment will also be briefly outlined.
        Speaker: Dr Jeremy Green (Institut für Kernphysik, Johannes Gutenberg-Universität Mainz)
        Slides
      • 47
        Study of the hadronic contributions to the running of the QED coupling
        The electromagnetic coupling receives significant contributions to its running from non-perturbative QCD effects. We present an update of our study of the Adler function and of its application to the determination of leading order hadronic contribution to the running of the QED coupling. We use a high-statistics lattice QCD computation with two flavours of O(a) improved Wilson fermions in a broad range of the momentum transfer Q^2. The running of the QED coupling, including valence contributions from u, d, s and c quarks, is compared to phenomenological results at intermediate Q^2 values.
        Speaker: Gregorio Herdoiza (IFT UAM-CSIC)
        Slides
      • 48
        The hadronic vacuum polarization function with O(a)-improved Wilson fermions - an update
        We present an update of the analysis of the leading order hadronic contribution to the muon g-2 using a high statistics evaluation of the vacuum polarization function (VPF) with O(a)-improved $N_f=2$ Wilson fermions. Partially twisted boundary conditions are used to increase the number of available $Q^2$ points. We employ an extended frequentist's method to study the systematic errors arising from the $Q^2$ dependence of the VPF as well as various ans\"atze for the continuum and chiral extrapolation. We present preliminary results for $a_\mu^\mathrm{HLO}$ including the valence contributions of u,d,s and c quarks.
        Speaker: Hanno Horch (Uni Mainz)
        Slides
      • 49
        Finite volume effects in hadronic vacuum polarization
        We investigate finite volume effects in the hadronic vacuum polarization, with an eye toward the corresponding systematic error in the muon anomalous magnetic moment. While it is well known that leading-order chiral perturbation theory does not provide a good description of the hadronic vacuum polarization, it turns out that it gives a much better representation of finite volume effects. Indications are that finite volume effects cannot be ignored when the aim is a few percent level accuracy for the hadronic contribution to the muon anomalous magnetic moment, even when $m_\pi L \sim 4$ and $m_\pi \sim 200$ MeV.
        Speaker: Prof. Maarten Golterman (San Francisco State University)
        Slides
      • 50
        A derivative based approach for the leading order hadronic contribution to g-2 of the muon
        We describe a lattice approach to calculate the leading-order anomalous magnetic moment of the muon. We employ lattice momentum derivatives to determine the hadronic vacuum polarization scalar at low momenta and construct a smooth, intregrable function in this momentum region. We present preliminary results for hex-smeared Wilson-quark lattice ensembles.
        Speaker: Dr Eric Gregory (Bergische Universitaet Wuppertal & Juelich Supercomputing Centre)
        Slides
      • 51
        Leading isospin breaking correction to the hadronic vacuum polarisation
        Several leading lattice collaborations are investing significant effort to reduce the uncertainty in the lattice computation of the HVP to a sub-percent level. In order to achieve this goal, giving an estimate of the size of the isospin breaking effects becomes relevant. Recently, the RM123 developed a method for separating QCD from QED isospin breaking effects based on expansion of the path–integral in powers of the up and down quark mass difference and the fine structure constant. We discuss how to apply this method to the connected part of the hadronic vacuum polarisation from the lattice.
        Speaker: Dr Marina Marinkovic (CERN)
        Slides
    • Nonzero Temperature and Density 401

      401

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Sourendu Gupta
      • 52
        Many flavor approach to study the nature of chiral phase transition of two-flavor QCD
        We perform lattice numerical simulations to study the phase transition of QCD at finite temperature and the vanishing chemical potential with intent to clarify the nature of the transition of massless two flavor QCD. We investigate QCD with two light and N_f heavy quarks instead of two-flavor QCD, and focus on the light quark mass dependence of the critical heavy mass below which the transition is of first order. The nature of the transition is identified by the shape of the constraint effective potential, constructed from the histogram of the generalized plaquette, at the critical temperature. The heavy quark effects are incorporated in the form of the hopping parameter expansion through the re-weighting technique. Our result indicates that the critical heavy mass remains finite in the chiral limit of the two flavors of light quarks, suggesting the phase transition of massless two-flavor QCD is of second order.
        Speaker: Dr Norikazu Yamada (KEK)
        Slides
      • 53
        The nature of the Roberge-Weiss Transition in N_f=2 QCD with Wilson Fermions on N_t=6 lattices
        The chiral and deconfinement phase transitions at zero density for light and heavy quarks, respectively, have analytic continuations to imaginary chemical potential. At some critical imaginary chemical potential, they meet the high temperature Roberge-Weiss transition between adjacent Z3 sectors. For light and heavy quarks, where the chiral and deconfinement transitions are first order, the transition lines meet in a triple point. For intermediate masses chiral or deconfinement transitions are crossover and the Roberge-Weiss transition ends in a second order point. At the boundary between these regimes the junction is a tricritical point, as shown in studies with N_f=2,3 flavors of staggered and Wilson quarks on Nt=4 lattices. Employing finite size scaling we investigate the nature of this point as a function of quark mass for N_f=2 flavors of Wilson fermions with a temporal lattice extent of Nt=6. In particular we are interested in the change of the location of tricritical points compared to our ealier study on Nt=4.
        Speaker: Mr Christopher Czaban (Goethe University Frankfurt am Main)
        Slides
      • 54
        The $N_f=2$ chiral phase transition from imaginary chemical potential with Wilson Fermions
        The clarification of the order of the thermal transition in the chiral limit of QCD with two dynamical flavours of quarks is a long-standing issue. Still, it is not definitely known if the transition is of first or second order. Which of the two scenarios is realized has important implications for the physical QCD phase diagram, and in particular it is important regarding the existence of a critical endpoint at finite densities. Settling this issue by simulating at successively decreased pion mass was not conclusive yet, primarily because of the increasing demands of the simulations as the pion mass is lowered. Recently, an alternative approach was proposed, which relies on the nontrivial phase structure of QCD at purely imaginary chemical potential $\mu_I$ induced by the Roberge-Weiss symmetry. At critical values of $\mu_I$, the endpoint of the Roberge-Weiss transition is mass-dependent and changes from a triple point at low and high masses to a second order endpoint for intermediate masses. These regions are separated by tricritical points. Leaving the critical $\mu_I$-values, it is known that lines of second order transitions depart from the tricritical points, separating regions of first order transitions from crossover regions. Furthermore, at $\mu_I$ the sign problem is absent and standard simulation algorithms can be applied. Hence, determining the second order line for lower quark masses is possible. The line is governed by tricritical scaling laws, which then allow an extrapolation to the chiral limit. In this way one can clarify the order of the chiral limit at zero chemical potential. Using staggered fermions on $N_t=4$ lattices, it was indeed found that the transition is of first order in the chiral limit. These findings have to be contrasted with other fermion discretizations. In this talk, we report on the status of our simulations with Wilson fermions following the same approach.
        Speaker: Dr Christopher Pinke (Goethe University Frankfurt)
        Slides
      • 55
        Analytic continuation of finite density QCD with heavy quarks in the strong coupling region
        Finite density and temperature QCD with heavy quarks is studied in the complex chemical potential ($¥mu$) plane by use of a mean field method and Monte Carlo simulations, where the former applies to the strong coupling region. We calculate the effective potential as a function of Polyakov line, and study thermodynamic singularities and their associated Stokes lines in the complex $¥mu$ plane. We also perform an explicit analytic continuation of the first order transition and crossover lines appearing on the real $¥mu$ axis in the strong coupling region.
        Speaker: Prof. Hiroshi Yoneyama (Saga University)
        Slides
      • 56
        Topological feature and phase diagram of QCD at complex chemical potential
        One of the important subjects in QCD is understanding properties of deconfinement transition at finite temperature and density. We investigate the pseudo-critical temperature of the deconfinement transition by considering the complex chemical potential. The imaginary chemical potential can be interpreted as the Aharonov-Bohm phase induced by U(1) flux insertions to the fictitious hole of the imaginary-time direction, then the analogy of topological order suggests that we can determine the deconfinement transition temperature by the Roberge-Weiss endpoint. We also present the expected QCD phase diagram based on the perturbative calculation at finite complex chemical potential.
        Speaker: Dr Kouji Kashiwa (Yukawa Institute for Theoretical Physics, Kyoto University)
        Slides
      • 57
        Finite-temperature phase transition of Nf=3 QCD with exact center symmetry
        By imposing flavor-dependent boundary conditions on quarks, we construct three-flavour SU(3) gauge theory reconciling fundamental fermion representation and exact Z3 center symmetry. In this talk, we show the first result of lattice simulation on this QCD-like theory, which we call Z3-QCD, with emphasis on the finite-temperature phase transition with respect to center and chiral symmetries. On the lattice, we formulate Z3-symmetric SU(3) gauge theory with three fundamental Wilson quarks by twisting quark boundary conditions in a compact dimension. We calculate the finite-temperature vacuum expectation value of Polyakov loop and the chiral condensate. We find out the first-order center phase transition where the hysteresis of temperature dependence exists depending on cold and hot starts. We also discuss the chiral crossover transition, and its relation to the center phase transition.
        Speaker: Prof. Tatsuhiro Misumi (Akita University / Keio University)
        Slides
    • Physics Beyond the Standard Model 405

      405

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Claudio Rebbi
      • 58
        A minimal model of the composite Higgs and its Goldstone dynamics
        A low mass scalar state has been observed by the Lattice Higgs Collaboration as a viable candidate for the minimal realization of the composite Higgs from a fermion doublet in the two-index symmetric representation of the SU(3) color gauge group. The difficulty of decoupling Goldstone dynamics from the low mass scalar state in realistic simulations requires the extension of the low energy effective theory for chiral symmetry breaking.
        Speaker: Prof. Julius Kuti (U.C. San Diego)
        Slides
      • 59
        Hadron Spectroscopy with a low-mass composite scalar in the sextet model
        The SU(3) gauge theory with two flavors of two-index symmetric (sextet) fermions has been of great interest recently since a low-mass scalar meson that may serve as a composite Higgs candidate was observed in our previous studies. The hadron spectrum analysis on a more extended dataset will be presented as a followup investigation.
        Speaker: Mr wong chik him (university of wuppertal)
        Slides
      • 60
        Taste symmetry restoration in the sextet model with staggered fermions
        We show the cutoff dependence of taste symmetry breaking on the full staggered Goldstone spectrum in the sextet model. Taste symmetry in the Goldstone valence spectrum is restored using a mixed action setup on the gauge field gradient flow.
        Speaker: Mr santanu Mondal (Eötvös Loránd University)
        Slides
      • 61
        Running coupling of the sextet composite Higgs model
        The renormalized running coupling of SU(3) gauge theory coupled to Nf = 2 flavors of massless Dirac fermions in the 2-index-symmetric (sextet) representation is calculated. This model is of particular interest as a minimal realization of the strongly interacting composite Higgs scenario. A recently proposed finite volume gradient flow scheme is used. The calculations are performed at several lattice spacings and two discretizations allowing for a controlled continuum extrapolation and particular attention is paid to estimating the systematic uncertainties.
        Speaker: Daniel Nogradi (Eotvos University)
        Slides
      • 62
        Conformal symmetry vs. chiral symmetry breaking in the SU(3) sextet model
        We present new results from our simulations of the SU(3) "sextet model" with two flavors in the 2-index symmetric representation of the gauge group. The simulations use unimproved Wilson fermions and we measure the meson and baryon spectrum of the theory for multiple bare quark masses at two different lattice spacings. To address the ongoing issue of whether the model is inside or outside the conformal window, we compare the spectrum to the expectations for a theory with spontaneous chiral symmetry breaking and to those of an IR conformal theory.
        Speaker: Mr Martin Hansen (CP3-Origins, University of Southern Denmark)
        Slides
    • Standard Model Parameters and Renormalization 404

      404

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Stefan Schaefer
      • 63
        Non-perturbative renormalization of the static quark theory in a large volume
        We report on progress to renormalize non-perturbatively the static heavy quark theory on the lattice. In particular, we present first results for position-space renormalization scheme for heavy-light bilinears. We test our approach on RBC's 16^3 x 32 lattice ensemble with pion mass of 420 MeV, Iwasaki gauge action and domain wall light fermions.
        Speaker: Dr Piotr Korcyl (University of Regensburg)
        Slides
      • 64
        Implementation of a non-perturbative matching strategy between heavy-light currents in HQET and QCD
        We outline our strategy to non-perturbatively match all components of the heavy-light axial and vector currents in Heavy Quark Effective Theory (HQET) at O(1/m) to finite-volume lattice QCD. Based on a tree-level study, we propose a set of matching conditions between suitable observables defined in QCD and HQET to fix the parameters of the effective theory, which are required to absorb the power divergences of lattice HQET. These conditions can be evaluated through numerical simulations, and we report on the status of our implementation in two-flavour QCD. The results of this finite-volume matching strategy will enter a HQET computation of the form factors of semi-leptonic B→π decays as a first application.
        Speaker: Mr Christian Wittemeier (University of Muenster, Germany)
        Slides
      • 65
        Perturbative versus non-perturbative decoupling of heavy quarks
        We simulate a theory with N_f=2 heavy quarks of mass M. At energies much smaller than M the heavy quarks decouple and the theory can be described by an effective theory which is a pure gauge theory to leading order in 1/M. We present results for the mass dependence of ratios such as t_0(M)/t_0(0). We compute these ratios from simulations and compare them to the perturbative prediction. The latter relies on a factorization formula for the ratios which is valid to leading order in 1/M in the low energy effective theory.
        Speaker: Prof. Francesco Knechtli (University of Wuppertal)
        Slides
      • 66
        Nonperturbative renormalization in the RI-SMOM scheme and Gribov copies for staggered bilinears
        We present renormalization factors for the bilinear operators obtained using the nonperturbative renormalization method (NPR) in RI-SMOM scheme for improved staggered fermions on the MILC asqtad lattice (Nf=2+1). We compare the RI-SMOM result to the RI-MOM scheme result and the one-loop perturbative result. Since the NPR requires Landau gauge fixing, we study related Gribov copy problem in staggered NPR analysis.
        Speaker: Mr Sungwoo Park (Seoul National University)
        Slides
      • 67
        NPR determination of quark masses from the HISQ action
        I will report on a calculation of quark bilinear Z-factors for the HISQ action using non-perturbative renormalization techniques. One of the main outputs is the mass renormalization factor Z_m relating bare quark masses to the MSbar scheme. This will provide an independent determination of quark masses in addition to other methods being used by the HPQCD collaboration.
        Speaker: Dr Andrew Lytle (University of Glasgow)
        Slides
      • 68
        Electromagnetic effects on the light pseudoscalar mesons and determination of $m_u/m_d$
        The MILC Collaboration has completed production running of electromagnetic effects on light mesons using asqtad improved staggered quarks. In these calculations, the photons are quenched using the non-compact formalism. Four lattice spacings from 0.12 fm to 0.045 fm have been used. Finite volume effects with $a=0.12$ fm have been studied with spacial sizes $L_s=12$, 16, 20, 28, 40, and 48. The chiral-continuum fit of the meson masses allows calculation of corrections to Dashen's theorem and determination of the ratio of up to down quark masses.
        Speaker: Prof. Steven Gottlieb (Indiana University)
        Slides
    • Weak Decays and Matrix Elements 406

      406

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Amarjit Soni
      • 69
        Chiral behavior of light meson form factors in 2+1 flavor QCD with exact chiral symmetry
        We present a study of the chiral behavior of light meson form factors in QCD with three flavors of overlap quarks. Gauge ensembles are generated at single lattice spacing 0.12fm with pion masses down to 300 MeV. The pion and kaon electro-magnetic form factors and the kaon semileptonic form factors are precisely calculated using the all-to-all quark propagator. We discuss their chiral behavior using the NNLO chiral perturbation theory.
        Speaker: Dr Takashi Kaneko (KEK)
        Slides
      • 70
        QED Corrections to Hadronic Processes in Lattice QCD
        We review the method proposed in [1] to compute electromagnetic effects in hadronic processes, such as decays, using lattice simulations. The method can be applied, for example, to the leptonic and semileptonic decays of light or heavy pseudoscalar mesons. For these quantities the presence of infrared divergences in intermediate stages of the calculation makes the procedure much more complicated than is the case for the hadronic spectrum, for which calculations already exist. In order to compute the physical widths, diagrams with virtual photons must be combined with those corresponding to the emission of real photons. Only in this way do the infrared divergences cancel as rst understood by Bloch and Nordsieck in 1937. We present a detailed analysis of the method for the leptonic decays of a pseudoscalar meson and also review the status of exploratory numerical studies. [1] N. Carrasco, V. Lubicz, G. Martinelli, C. T. Sachrajda, N. Tantalo, C. Tarantino and M. Testa, "QED Corrections to Hadronic Processes in Lattice QCD," Phys. Rev. D 91 (2015) 7, 074506 [arXiv:1502.00257 [hep-lat]].
        Speaker: Prof. Christopher Sachrajda (University of Southampton)
        Slides
      • 71
        Neutral Kaon mixing beyond the Standard Model
        We compute the hadronic matrix elements of the four-quark operators needed for the study of $K^0-{\bar K^0}$ mixing beyond the Standard Model (SM). We have used $n_f=2+1$ flavours of domain wall fermion (DWF) at two values of the lattice spacing ($a\approx0.08$ and $a\approx0.11\text{fm}$) and with lightest unitary Pion mass of $\approx300\,\text{MeV}$. Renormalisation is performed non-perturbatively and the impact of several different intermediate momentum schemes is investigated.
        Speaker: Dr Renwick Hudspith (York University, Toronto)
        Slides
      • 72
        Improved Hadronic Matrix Element Determination Using the Variational Method
        Calculating hadronic matrix elements in Lattice QCD that are relevant for e.g. Form Factors and Parton Distribution Functions, provides important insights into the structure of hadrons. However standard 3-point function methods have their limitations. One of the most commonly studied sources of systematic error is excited state contamination, which occurs when correlators are contaminated with results from higher energy excitations. This investigation uses configurations generated by the QCDSF/UKQCD/CSSM collaborations at the SU(3)-symmetric point (pion mass of 460 MeV) a lattice spacing of 0.074 fm with a volume of 32^3 x 64. We apply the variational method to calculate a range of quantities and compare the results to the more commonly used summation and two-state fit methods. The results of this analysis demonstrate that the variational approach offers a more efficient and robust method for the determination of nucleon matrix elements.
        Speaker: Mr Jack Dragos (University of Adelaide)
        Slides
      • 73
        Nuclear Parity Violation from Lattice QCD
        The electroweak interaction on the level of quarks and gluons are well understood from precision measurements in high energy collider experiments. Relating these fundamental parameters to Hadronic Parity Violation in nuclei however remains an outstanding theoretical challenge. One of the most interesting observables in this respect is the parity violating hadronic neutral current: it is hard to measure in collider experiments and is thus the least constrained observable of the Standard Model. Precision measurements of parity violating transitions in nuclei can help to improve these constraints. In these systems however, the weak interaction is masked by effects of the ~7 orders of magnitude stronger non-perturbative strong interaction. Therefore, in order to relate experimental measurements of the parity violating pion-nucleon couplings to the fundamental Lagrangian of the SM, these non-perturbative effects have to be well understood. In this talk, I am going to present a Lattice QCD approach for computing the $\Delta I=2$ parity violating matrix element in proton proton scattering. This process does not involve disconnected diagrams in the isospin symmetric limit and is thus a perfect testbed for studying the feasibility of the more involved calculation of the parity violating pion-nucleon coupling. In this talk, I am going to present our technology for computing this quantity and some preliminary results.
        Speaker: Dr Thorsten Kurth (Lawrence Berkeley National Laboratory)
        Slides
    • Plenary Session BF & 1F Main Hall

      BF & 1F Main Hall

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Steven Gottlieb
      • 74
        Hadron Structure from Lattice QCD
        Here we review the progress made in understanding the internal structure of hadrons in terms of the gluonic and quark constituents. Recent results for standard observables such as the nucleon axial charge, electromagnetic form factors and quark momentum fraction will be summarised, before turning our attention to more challenging quantities, including quark disconnected contributions and so-called parton quasi-distribution functions. With many new experiments whose primary goal is to understand the internal dynamics of hadrons now underway, particular emphasis will be given to lattice calculations that promise to assist and guide these experimental efforts.
        Speaker: Dr James Zanotti (University of Adelaide)
        Slides
      • 75
        Light flavours
        For quite some time now simulations of lattice QCD have allowed for predicting a basic set of light flavour quantities reliably and with increasingly high precision. The field has started to move on: Advances in field theory, algorithms and computing for the first time allow to address more complicated problems like for example hadronic and rare kaon decays, the kaon mass-difference or the conceptually clean inclusion of electromagnetic and isospin effects. This talk aims at providing an overview over the state-of-the-art.
        Speaker: Dr Andreas Juettner (University of Southampton)
        Slides
    • 10:30
      Coffee Break
    • Plenary Session BF & 1F Main Hall

      BF & 1F Main Hall

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Sinya Aoki
      • 76
        Kenneth G. Wilson Award for Excellence in Lattice Field Theory
        Slides
      • 77
        Three-body observables from the lattice
        Scattering and transition amplitudes with three-hadron final states play an important role in nuclear and particle physics. In order to predict such quantities using Lattice QCD, formalism is required to overcome the limitations of Euclidean time and finite volume. In this talk I will focus on extensions of Luescher's work relating the finite-volume energy spectrum to physical scattering amplitudes. I will highlight the challenges that arise in extending the formalism from two- to three-particle states, and will describe how these have been addressed. Finally, I will outline outstanding problems and discuss the prospects of applying the formalism in a numerical Lattice QCD calculation.
        Speaker: Dr Maxwell Hansen (University of Mainz)
        Slides
      • 78
        Standard-model prediction for direct CP violation in K->pi pi decays
        We discuss our recent publication [arXiv:1505.07863] of the first lattice QCD calculation of the complex kaon decay amplitude A_0 with physical kinematics, using a single 32^3 x 64 domain wall ensemble with G-parity spatial boundary conditions. We obtain approximate agreement with the experimental value for Re(A_0), which serves as a test of our method. Our prediction of Im(A_0) can be used to compute the direct CP violating ratio Re(ε′/ε), which we find to be ~2 sigma lower than the experimental value. This result provides a new test of the Standard Model theory of CP violation, one which can be made more accurate with increasing computer capability.
        Speaker: Dr Christopher Kelly (RBRC, Brookhaven National Laboratory)
        Slides
    • 12:30
      Lunch Break
    • Hadron Spectroscopy and Interactions 402

      402

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Michael Endres
      • 79
        Magnetic properties of light nuclei and the $n p \to d \gamma$ transition
        I will discuss the NPLQCD collaboration's recent calculations of the magnetic moments and polarisabilities of light nuclei up to A=4 using background field methods. I will also describe an a calculation of the magnetic transition amplitude that determines the $n p \to d \gamma$ cross section that is of central importance in Big Bang Nucleosynthesis.
        Speaker: Dr William Detmold (MIT)
        Slides
      • 80
        Light nuclei and nucleon form factors in N_f=2+1 lattice QCD
        We present our result of binding energy of light nuclei with the nuclear mass number less than or equal to four at the pion mass m_pi = 0.3 GeV. The simulations are performed in N_f=2+1 QCD with Iwasaki gauge and non-perturbative improved Wilson fermion actions at the lattice spacing of a = 0.09 fm. We discuss the quark mass dependence of the binding energies using our previous results and also a preliminary result at almost physical pion mass m_pi ~ 0.145 GeV with a ~ 0.085 fm. Furthermore, we show a preliminary result of the axial charge and the radii obtained from isovector nucleon form factors at the almost physical pion mass.
        Speaker: Takeshi Yamazaki (University of Tsukuba)
        Slides
      • 81
        The H-dibaryon in two flavor lattice QCD
        We present preliminary results from a lattice QCD calculation of the H-dibaryon using two flavours of O(a) improved Wilson fermions. We employ six-quark interpolating operators with the appropriate quantum numbers of the H-dibaryon and also explore its couplings to two-baryon channels. To improve the overlap to the ground state two smearings are employed and a generalised eigenvalue problem is solved in the aforementioned operator basis. With the application of L\"uscher's finite volume formalism, we explore the nature of the infinite volume interaction of the two baryons. The relevant correlators are projected to three moving frames further enabling the isolation of the infinite volume bound/scattering state. Preliminary results on pion mass of 1 $\mathrm{GeV}$ indicate the H-dibaryon is bound in the infinite volume. Results at a lower pion mass of 451 MeV will also be presented.
        Speaker: Dr Parikshit Junnarkar (Helmholtz Institut Mainz)
        Slides
      • 82
        Two-nucleon scattering in multiple partial waves
        We determine scattering phase shifts for s,p,d, and f partial wave channels in two-nucleon systems using lattice QCD methods. We use a generalization of Luscher's finite volume method to determine infinite volume phase shifts from a set of finite volume ground- and excited-state energy levels on two volumes, V=(3.4 fm)^3 and V=(4.5 fm)^3. The calculations are performed in the SU(3)-flavor limit, corresponding to a pion mass of approximately 800 MeV. From the energy dependence of the phase shifts we are able to extract scattering parameters corresponding to an effective range expansion.
        Speaker: Dr Amy Nicholson (UC Berkeley)
        Slides
      • 83
        Calculation of the decay width of decuplet baryons
        We present an update of our investigation of hadronic decays of decuplet baryons using the transfer matrix method. In addition to the benchmark decay Delta -> pi N we look at other channels, like Sigma^{*} -> pi Sigma and Sigma^{*} -> pi Lambda to further test the reach of the method. Lattice calculations are done using a hybrid setup with domain wall fermions on a staggered sea at pion mass 350 MeV and using a unitary setup with domain wall fermions at pion mass 180 MeV.
        Speaker: Dr Marcus Petschlies (HISKP, University of Bonn)
        Slides
      • 84
        The Roper resonance from spatially large interpolation fields
        We present our calculations of the Roper state obtained from a combination of spatially large interpolation fields. The calculation is carried out with overlap valence fermion on $2 + 1$-flavor domain-wall fermion configurations on the $24^3\times 64$ lattice with $a^{-1} = 1.73\,$GeV. Our result is consistent with that from the Sequential Bayesian method (SBM) on the same lattice, and with the experimental value at $1440$ MeV. We utilize the same method on the anisotropic Clover lattice, and compare with the overlap results. In the end we give an explanation on the nature of the difference.
        Speaker: Mingyang Sun (University of Kentucky)
        Slides
    • Hadron Structure 403

      403

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Paul Rakow
      • 85
        Parton Distribution Function from Hadronic Tensor
        Parton Distribution function can be obtained from calculating the hadronic tensor on the lattice. It involves a 4-point function evaluation with two vector or axial-vector currents. We propose to use the source method for one of the currents and effectively reduce the numerical complexity to that of a 3-point function. The method can be generalized to evaluate other 4-point functions.
        Speaker: Prof. Keh-Fei Liu (University of Kentucky)
        Slides
      • 86
        Nucleon transverse momentum-dependent parton distributions: Comparing Clover and Domain wall fermion results at ~300 MeV pion mass
        We present a lattice QCD calculation of transverse momentum dependent parton distributions (TMDs) in the proton using staple-shaped gauge links. We calculate both naively time-reversal odd (T-odd) observables, namely, the generalized Sivers and Boer-Mulders transverse momentum shifts in SIDIS and DY processes, as well as T-even observables, namely, the transversity related to the tensor charge and the generalized worm-gear shift. The calculation is done on a n_f=2+1 clover ensemble with lattice spacing 0.114 fm and pion mass 317 MeV. The results are compared with a previous calculation on a domain-wall ensemble at 0.084 fm lattice spacing with 297 MeV pion mass.
        Speaker: Dr Boram Yoon (Los Alamos National Laboratory)
        Slides
      • 87
        Nucleon generalized form factors from lattice QCD with nearly physical quark masses
        We show results of generalized form factors of the nucleon from lattice simulations with Nf = 2 mass-degenerate non-perturbatively improved Wilson-Sheikholeslami-Wohlert fermions down to a pion mass of 150 MeV. We also present the resulting isovector quark angular momentum. Possible excited state contaminations are treated with correlated simultaneous fits to all 3pt functions of a given ensemble with fixed momentum transfer squared.
        Speaker: Mr Rudolf Roedl (Universität Regensburg)
        Slides
    • Nonzero Temperature and Density 401

      401

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Atsushi Nakamura
      • 88
        Exploring Complex-Langevin Methods for Finite-Density QCD
        Recent advances in applying complex-Langevin simulations to gauge theories with complex actions, in particular gauge-cooling, have revived interest in using these methods for QCD at finite baryon-number density. QCD at a finite chemical potential $\mu$ for quark number has a complex fermion determinant which prevents use of standard simulation techniques. Complex-Langevin simulations show promise for simulating this theory. We describe our preliminary investigations in applying complex-Langevin simulations to lattice QCD at finite $\mu$. We use a complex extension of the partial second-order approach to Langevin simulations of Fukugita, Ukawa and others, using staggered fermions. Our current studies are of 2-flavour lattice QCD at finite $\mu$ at zero temperature on a $12^4$ lattice. Of special interest is if we can observe the expected transition at $\mu \approx m_N/3$, and whether there is a spurious transition at $\mu \approx m_\pi/2$. Finite-temperature simulations at finite $\mu$ are also planned.
        Speaker: d. k. sinclair (argonne national laboratory)
        Slides
      • 89
        Towards the heavy dense QCD phase diagram using Complex Langevin simulations
        Monte Carlo methods cannot probe far into the QCD phase diagram with a real chemical potential, due to the famous sign problem. Complex Langevin simulations, using adaptive step-size scaling and gauge cooling, are suited for sampling path integrals with complex weights. We report here tests on the deconfinement transition in pure Yang-Mills SU(3) simulations and present an update on the QCD phase diagram in the limit of heavy and dense quarks.
        Speaker: Mr Felipe Attanasio (Swansea University)
        Slides
      • 90
        Insights into the heavy dense QCD phase diagram using Complex Langevin simulations
        Complex Langevin simulations provide an alternative to sample path integrals with complex weights and therefore are suited to determine the phase diagram of QCD from first principles. Adaptive step-size scaling and gauge cooling are used to make correct convergence possible. We present results for the phase diagram of QCD in the limit of heavy quarks and discuss the order of the phase transitions, which are studied by varying the spatial simulation volume.
        Speaker: Dr Benjamin Jaeger (Swansea University)
        Slides
      • 91
        Testing a generalized cooling procedure in the complex Langevin simulation of chiral Random Matrix Theory
        The complex Langevin (CL) method has been attracting much attention as a solution to the sign problem since the method was shown to work in finite density QCD in the deconfined phase by using the so-called gauge cooling procedure. Whether it works also in the confined phase with light quarks is still an open question, though. In order to shed light on this question, we apply the method to the chiral Random Matrix Theory (RMT), which describes the epsilon regime of finite density QCD. Earlier works reported that a naive implementation of the method fails to reproduce the known exact results and that the problem can be solved by choosing a suitable coordinate. In this work we stick to the naive implementation, and show that a generalized gauge cooling procedure can be used to avoid the problem.
        Speaker: Dr Keitaro Nagata (KEK)
        Slides
      • 92
        Understanding the problem with logarithmic singularities in the complex Langevin method
        In recent years, there has been remarkable progress in theoretical justification of the complex Langevin method, which is a promising method for evading the sign problem in the path integral with a complex weight. There still remains, however, an issue concerning occasional failure of this method in the case where the action involves logarithmic singularities such as the one appearing from the fermion determinant in finite density QCD. In this talk, we point out that this failure is due to the breakdown of the relation between the complex weight which satisfies the Fokker-Planck equation and the probability distribution generated by the stochastic process. In fact, this kind of failure can occur in general when the stochastic process involves a singular drift term. We show, however, in simple examples, that there exists a parameter region in which the method works although the standard reweighting method is hardly applicable.
        Speaker: Dr Shinji Shimasaki (KEK)
        Slides
      • 93
        Complex Langevin in low-dimensional QCD: the good and the not-so-good
        We present our latest results on the application of the Complex Langevin method in one- and two-dimensional QCD. Although the method is stable, it unfortunately converges to an incorrect result when applied as such. After applying additional gauge cooling steps, the results agree with the known analytical results in the one-dimensional case. However, in the two-dimensional case the disagreement subsists, even with gauge cooling, when the sign problem is sufficiently large.
        Speaker: Dr Jacques Bloch (University of Regensburg)
        Slides
    • Physics Beyond the Standard Model 405

      405

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Julius Kuti
      • 94
        Lattice study for conformal windows of SU(2) and SU(3) gauge theories with fundamental fermions
        We present results from our study of SU(2) gauge theory with eight flavours, and SU(3) gauge theory with twelve flavours. These two theories may be very close to the lower ends of the conformal windows for the SU(2) and SU(3) gauge groups, respectively, when only the fundamental-representation fermions are present. For the SU(2) theory with eight flavours, we report our investigation for the distribution of the lowest-lying eigenvalues of the Dirac operator. In particular, we compare our numerical results with predictions from Random Matrix Theory and extract the chiral condensate. As for the SU(3) theory with twelve flavours, we show our final analysis for the Gradient-Flow running coupling constant using the step-scaling method. In this presentation, we demonstrate our detailed study of the continuum extrapolation in the step-scaling approach, and comment on the challenges in such computations.
        Speaker: Prof. C.-J. David Lin (National Chiao-Tung University)
        Slides
      • 95
        Approaching the conformal window: systematic study of the particle spectrum in SU(2) field theory with N_f =2,4 and 6.
        It is expected that SU(2) gauge theory with N_f fundamental fermions has an infrared fixed point when N_f is between ~6 and 11. We study the hadron spectrum and scale setting in SU(2) gauge field theory with N_f=2,4,6 using hypercubic stout smeared Wilson-clover (HEX) action. The case N_f=2 is QCD-like, whereas N_f=6 is close to the lower edge of the conformal window. In our study length scales are determined by using gradient flow approach.
        Speaker: Sara Tähtinen (University of Helsinki, Helsinki Institute of Physics)
        Slides
      • 96
        Gradient flow and IR fixed point in SU(2) with Nf=8 flavors
        SU(2) with Nf=8 flavors of fundamental fermions is expected to feature an infrared fixed point (IRFP). We measure the evolution of the coupling constant with Schrödinger functional boundary conditions and gradient flow, using HEX-smeared Wilson-clover action. We observe clear evidence for a fixed point, qualitatively compatible with perturbative results.
        Speaker: Mr Viljami Leino (University of Helsinki)
        Slides
      • 97
        Lattice simulations of technicolour theories with adjoint fermions and supersymmetric Yang-Mills theory
        In this talk I will present our newest results from the lattice simulation of SU(2) adjoint QCD with two Dirac flavours, also known as minimal walking technicolour. These results include the mass spectrum of the theory, in particular mesonic states with disconnected contributions and a specific spin half fermion-gluon state of the theory. I will discuss the results in comparison with SU(2) adjoint QCD with one Dirac flavour and supersymmetric Yang-Mills theory.
        Speaker: Dr Georg Bergner (ITP AEC University of Bern)
        Slides
      • 98
        Wilson Fermions with Four Fermion Interactions
        We present a study of a chirally symmetric model of four fermion interactions, the Nambu Jona-Lasinio model, on the lattice with Wilson fermions. Four fermion interactions are a necessary part of many models of beyond Standard Model physics. In particular, they can couple the Standard Model fermions to a strongly interacting Higgs sector. In technicolor models they may cause spontaneous chiral symmetry breaking and modify the running coupling. As a first step, we study the restoration and spontaneous breaking of chiral symmetry in the lattice NJL model before adding a gauge interactions.
        Speaker: Dr Jarno Rantaharju (CP3-Origins)
        Slides
    • Standard Model Parameters and Renormalization: with Hadron Spectroscopy and Interactions 404

      404

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Patrick Fritzsch
      • 99
        Lattice input on the tau V_us puzzle
        Existing versions of the standard determinations of V_us from flavor-breaking finite energy sum rules with hadronic tau decay data as input yield V_us values 3 sigma or more below the expectations of 3-family unitarity. Slow convergence of the D=2 OPE series employed in these sum rules, however, makes it hard to assess the reliability of the use of the OPE representation. In this paper we use Euclidean Q^2 lattice data for the relevant flavor-breaking polarization function difference to assess the use and reliability of the OPE for this quantity. We then revisit the sum rule determination of V_us with the lessons learned from this study in mind. We show that previously encountered self-consistency problems are solved by the new analysis, and note that, with the the strange spectral distribution modified to account for the somewhat larger preliminary BaBar result for the K^- pi^0 branching fraction, the resulting output V_us is in good agreement with the result obtained from K_{ell 3} using lattice input for f_+(0).
        Speaker: Kim Maltman (York University)
        Slides
      • 100
        Scalar and vector form factors of $D \to \pi \ell \nu$ decays with $N_f=2+1+1$ Twisted fermions
        We present lattice results for the form factors relevant for $D \to \pi \ell \nu$ decays, obtained from simulations performed by the European Twisted Mass Collaboration with $N_f=2+1+1$ flavors of dynamical quarks, at three values of the lattice spacing and pion masses as low as $210$ MeV. We computed both the vector and the scalar form factors, studied their dependance on the momentum transfer and compared our results with the experimental ones. Specifically, by combining our determination of $f_+(0)$ with the experimental result for $ |f_+(0)V_{cd}|$ we are able to determine the CKM matrix element.
        Speaker: Mr Paolo Lami (Università degli Studi Roma Tre & INFN sezione Roma Tre)
        Slides
      • 101
        Nf=2+1+1 renormalisation of four-quark operators
        The phenomenology of several problems of flavour physics require the renormalisation of four-quark effective operators at the electroweak scale. While some part of the running can be done perturbatively, it is important to perform the non-perturbative renormalisation at a scale as high as possible. We will present results at an energy range where four flavours of quarks are present, and show how it improves our control of the perturbative matching.
        Speaker: Dr Julien Frison (University of Edinburgh)
        Slides
      • 102
        The leading hadronic contribution to $\gamma$-$Z$ mixing
        We study the hadronic contribution to the $\gamma$-$Z$ mixing, which determines the leading order hadronic contribution to the running of the electroweak mixing angle $\theta_W$. The required vacuum polarization function $\Pi^{\gamma Z}$ is calculated from the appropriate vector correlation functions in a mixed time-momentum representation. We explicitly calculate the connected and the disconnected contributions to such vector correlators using $N_f=2$ dynamical flavors of non-perturbatively $O(a)$-improved Wilson fermions.
        Speaker: Ms Vera Guelpers (University of Mainz)
        Slides
      • 103
        Determining the scale in Lattice QCD
        Numerical lattice QCD simulations determine mass (or other) ratios but not the scale itself, which has to be determined from experiment. A hadron mass such as the proton mass or decay constant such as the pion decay constant are often used for this purpose. We discuss here the advantages of setting the scale using a flavour-singlet quantity, which in conjunction with simulations keeping the average quark mass constant allow SU(3) flavour breaking expansions to be used. This is illustrated using 2+1 clover fermions, and a determination of the Wilson flow scales t0 and w0 is given.
        Speaker: Roger Horsley (University of Edinburgh)
        Slides
      • 104
        QCD spectroscopy and quark mass renormalisation in external magnetic fields with Wilson fermions
        We study the change of the QCD spectrum of low-lying mesons in the presence of an external magnetic field using Wilson fermions in the quenched approximation. Motivated by qualitative differences observed in the spectra of overlap and Wilson fermions for large magnetic fields, we investigate the dependence of the additive quark mass renormalisation on the magnetic field. To this purpose we derive Ward identities for lattice and continuum QCD+QED from which we can extract the current quark masses. We then compare different strategies of tuning for the quark masses.
        Speaker: Dr Bastian Brandt (University of Regensburg)
        Slides
    • Theoretical Developments 406

      406

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Yannick Meurice
      • 105
        Renormalizability of the Schrödinger Functional.
        Symanzik showed that quantum field theory can be formulated on a space with boundaries by including suitable surface interactions in the action to implement boundary conditions. We show that to all orders in perturbation theory all the divergences induced by these surface interactions can be absorbed by a renormalization of their coefficients.
        Speaker: Prof. Anthony Kennedy (University of Edinburgh)
        Slides
      • 106
        Convergent Perturbation Theory for $\phi^4$ model on lattice
        The small coupling expansions in the lattice field theory are asymptotic and do not converge. It was shown in previous studies, that an appropriate regularization of the integrals or shifting of the initial approximation allows one to construct convergent series. In this work we study the convergent series in application to the lattice $\phi^4$-model and compare the observables calculated using the series with the results of the Monte Carlo simulations. In addition, we consider the Borel resummation of the weak coupling lattice $\phi^4$-perturbation theory.
        Speaker: Mr Aleksandr Ivanov (Faculty of Physics M.V.Lomonosov Moscow State University)
        Slides
      • 107
        Induced YM theory with auxiliary bosons
        We study pure SU(N) lattice gauge theory with a plaquette weight factor given by an inverse determinant which can be written as an integral over auxiliary bosonic fields (modifying a proposal of Budczies and Zirnbauer). We derive conditions for the existence of a continuum limit and its equivalence to Yang-Mills theory. Furthermore, we perturbatively compute the relation between the coupling constants of the 'induced' gauge action and the familiar Wilson gauge action using the background-field technique. The perturbative relation agrees well with numerical results for N=2 in three dimensions.
        Speaker: Robert Lohmayer (University of Regensburg)
        Slides
      • 108
        Renormalization of two-dimensional XQCD
        Recently, Kaplan proposed an interesting extension of QCD : XQCD (=extended QCD) with bosonic auxiliary fields. While its partition function is kept exactly the same as that of QCD, XQCD naturally contains properties of low-energy hadrons. We apply this extension to the two-dimensional QCD in the large N_c limit ('t Hooft model). In this solvable model, it is possible to directly examine the hadronic picture of the 2d XQCD and analyze its renormalization group flow to understand how the auxiliary degrees of freedom behave in the low energy region. We confirm that the additional scalar fields can become dynamical acquiring the kinetic term, and its parity-odd part becomes dominant in the low energy region. This renormalization of XQCD provides an "extension" of the renormalization scheme of QCD, inserting different field variables from those in the original theory, without any changes in physical observables.
        Speaker: Mr Ryo Yamamura (Osaka U.)
        Slides
      • 109
        Combining ordinary and topological finite volume effects for fixed topology simulations
        In lattice simulations at fine lattice spacings typical algorithms tend to freeze topologically. In such cases specific topological finite size effects have to be taken into account, to obtain physical results corresponding to infinite volume or unfixed topology. Moreover, in QCD simulations, where the volume is not that large, it is often also necessary to get rid of ordinary finite volume effects not related to topology freezing. To extract physical results from simulations affected by both types of finite volume effects, we extend a known relation from the literature between hadron masses at fixed and at unfixed topology by incorporating additionally also ordinary finite volume effects. Numerical results for SU(2) Yang-Mills theory are presented.
        Speaker: Mr Arthur Dromard (Goethe University Frankfurt am Main)
        Slides
      • 110
        Three particles in a finite volume
        The spectrum of a bound state of three identical particles with a mass $m$ in a finite cubic box is studied within the effective field theory approach. It is shown that in the limit of a large two-body scattering length, the energy shift of a shallow bound state is given by $\Delta E=c (\kappa^2/m) (\kappa L)^{-3/2}|A|^2 \exp(−2\kappa L/\sqrt{3})$, where $\kappa$ is the bound-state momentum, $L$ is the box size, $|A|^2$ denotes the three-body analog of the asymptotic normalization coefficient of the bound state wave function and $c$ is a numerical constant. The formula is valid for $\kappa L\gg 1$. We further compare these predictions to the results of numerical calculations of the three-body spectrum in a finite volume. Using this approach to study the nature of the three-body bound states on the lattice is discussed.
        Speaker: Dr Akaki Rusetsky (HISKP, University of Bonn)
        Slides
    • 16:00
      Coffee Break
    • Hadron Spectroscopy and Interactions 402

      402

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: William Detmold
      • 111
        First results of baryon interactions from lattice QCD with physical masses (1) -- General overview and two-nucleon forces
        One of the most important frontier in lattice QCD is the first-principles determination of baryon-baryon interactions. Our novel method, the HAL QCD method, has been shown to be effective for this objective, while the calculations so far employed unphysically heavy quark masses. Under these circumstances, we have launched the new project to calculate baryon interactions in lattice QCD, employing quark masses around the physical point on a huge lattice volume of (8fm)^4. Resources such as K computer at AICS, Kobe, Japan are used for this computations. In this talk, we first give the overview of this on-going project, such as the theoretical formalism in HAL QCD method and numerical aspects of the baryon-force computations. We then turn to the numerical results, where first preliminary results for two-nucleon forces will be shown
        Speaker: Dr Takumi Doi (RIKEN)
        Slides
      • 112
        First results of baryon interaction from lattice QCD with physical masses (2) -- S=-3 and S=-4 sectors (XiXi, XiSigma, XiLambda-XiSigma channels) --
        Nucleaon-Nucleon interaction plays an essential role in nuclear physics. In the same way, hyperon-hyperon interactions shold play an important role in hyper nuclear physics. However, unlike the nucleons who are quite stable, hyperons decay quickly so that the direct scattering experiments are difficult. As a result, phenomenological determination of hyperon potentials involves large uncertainly. In this talk, we present our preliminary lattice QCD results on the determination of hyperon potentials in XiXi, XiSigma, and XiLambda-XiSigma channels based on the HAL QCD method by using a 2+1 flavor gauge configuration near the physical point generated by K computer.
        Speaker: Prof. Noriyoshi for HAL QCD Collaboration Ishii (Research Center for Nuclear Physics, Osaka University)
        Slides
      • 113
        First results of baryon interactions from lattice QCD with physical masses (3) -- Strangeness S=-2 two-baryon system
        The determination of baryon-baryon interactions directly from QCD is highly awaited. Especially for the strangeness S=-2 two-baryon system, we can investigate the SU(3) structure and its breaking effect of baryon-baryon interactions because the flavor singlet combination is allowed only in this system. Our approach is deriving a potential from coupled channel Schroedinger equation using Nambu-Bethe-Salpeter wave function measured on the lattice. We will report our latest results of the S=-2 baryon-baryon interactions and the fate of H-dibaryon by lattice QCD simulation employing quark masses around the physical point on a huge lattice volume of L=8fm generated by K computer at AICS.
        Speaker: Dr Kenji Sasaki (CCS, University of Tsukuba)
        Slides
      • 114
        Lattice QCD studies of baryon interactions from HAL QCD method and Luscher's finite volume method
        Both HAL QCD and Luscher's finite volume methods are used to calculate the hadron interactions from lattice QCD. However, there are some systematic deviations between these methods. In this talk, we investigate baryon interactions from both HAL QCD and Luscher's finite volume methods using the same lattice setups. From systematic comparison of these two techniques, we discuss the origin of the deviations between them and possibilities of the improvements of both methods.
        Speaker: Dr Takumi Iritani (Yukawa Institute for Theoretical Physics)
        Slides
      • 115
        Lambda_c-N interaction from lattice QCD
        Lambda_c (2286), the lightest baryon that contains one charm quark, has been found at some experiments, but its detailed properties such as the strength of the interaction with the nucleon, including a possibility of forming a deuteron-like bound state, are not known. Due to the heavy quark symmetry, heavy baryons such as Lambda_c (2286) sometimes show quite different behavior from light baryons, so the direct study of heavy baryons in QCD is important to reveal their properties. Recently, a new approach to investigate hadron interactions in lattice QCD has been proposed and developed extensively by the HAL QCD Collaboration. Employing this method, we have investigated the interaction between Lambda_c and nucleon (N). We first calculate the (single channel) potential between Lambda_c-N, with which we investigate an existence of the bound state in this system. We then consider the coupled channel potentials between Lambda_c-N and Sigma_c-N, in order to see effects from inelastic states to Lambda_c-N interactions. In this talk, we present the current status of our research project on Lambda_c-N interactions as well as future prospects.
        Speaker: Mr Takaya for HAL QCD Collaboration Miyamoto (Yukawa Institute for Theoretical Physics, Kyoto University)
        Slides
      • 116
        Zc(3900) from coupled-channel HAL QCD approach on the lattice
        We study the candidate of charmed tetraquark Zc(3900) from full QCD simulation. The Zc(3900) was first reported by both BESIII and Belle Collaborations in pi J/psi invariant mass spectrum. After that, structure of the Zc(3900) is actively discussed using many phenomenological models. However, due to lack of information of the interactions of the pi J/psi and its coupled systems, the predicted structures of the Zc(3900) are highly dependent on model parameters so that the information on the coupled-channel interactions are necessary to conclude the structure. We perform the coupled-channel analysis for the Zc(3900) through the HAL QCD method. We measure wave functions (NBS wave functions) on the lattice, and extract the potential matrix. Using the interactions obtained from LQCD, we investigate the pi J/psi and DbarD* invariant mass spectra of 2-body scatterings, the pole position of the scattering amplitudes on the complex energy plane, and production reaction of the Zc(3900) from Y(4260) decay. I will report those results for the Zc(3900).
        Speaker: Dr Yoichi Ikeda (RIKEN)
        Slides
    • Hadron Structure 403

      403

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Keh-Fei Liu
      • 117
        Pion electromagnetic form factor from full lattice QCD
        We present the first calculation of the pion electromagnetic form factor at physical light quark masses. This form factor parameterises the deviations from the behaviour of a point-like particle when a photon hits the pion. These deviations result from the internal structure of the pion and can thus be calculated in QCD. We use three sets (different lattice spacings) of n_f = 2+1+1 lattice configurations generated by the MILC collaboration. The Highly Improved Staggered Quark formalism (HISQ) is used for all of the sea and valence quarks. Using lattice configurations with u/d quark masses very close to the physical value is a big advantage, as we avoid the chiral extrapolation. We study the shape of the vector (f_+) form factor in the q^2 range from 0 to -0.15 GeV^2 and extract the mean square radius, <r^2_v>. The shape of the vector form factor and the resulting radius is compared with experiment. We also discuss the scalar form factor and radius extracted from that, which is not directly accessible to experiment. We have also calculated the contributions from the disconnected diagrams to the scalar form factor at small q^2 and discuss their impact on the scalar radius <r^2_s>.
        Speaker: Dr Jonna Koponen (University of Glasgow)
        Slides
      • 118
        SU(3)-breaking effects and induced second-class form factors in hyperon beta decays from 2+1 flavor lattice QCD
        We discuss the effects of SU(3) symmetry breaking measured in hyperon semileptonic decays from fully dynamical lattice QCD. Our calculations are carried out with gauge configurations generated by the RBC and UKQCD collaborations with (2+1)-flavors of dynamical domain-wall fermions and the Iwasaki gauge action at two couplings, ¥beta$ = 2.13 and 2.25. We have estimated the value of the hyperon vector couplings $f_1(0)$ with an accuracy of less than one percent. We then find that lattice results of $f_1(0)$ combined with the best estimate of $|V_{us}|$ with imposing CKM unitarity are slightly deviated from the experimental result of $|V_{us}f_1(0)|$ for $¥Sigma ¥rightarrow N$ decay. This discrepancy can be attributed to an assumption made in the experimental analysis on $|V_{us}f_1(0)|$, where the induced second-class form factor $g_2$ is set to be zero. We will therefore report on this matter and show the preliminary results of $g_2(0)$ evaluated in both indirect and direct ways.
        Speaker: Prof. Shoichi Sasaki (Tohoku University)
        Slides
      • 119
        Nucleon axial and tensor charges with dynamical overlap quarks
        We report on our calculation of the nucleon axial and tensor charges in 2+1-flavor QCD with dynamical overlap quarks. Numerical simulations are carried out at single lattice spacing 0.12 fm with pion masses below 500 MeV. We calculate both of the connected and disconnected contributions by using the all-to-all quark propagator. We present our preliminary results for the isoscalar and isovector charges.
        Speaker: Dr Nodoka Yamanaka (iTHES Group, Riken)
        Slides
      • 120
        High statistic analysis of nucleon form factors and charges in lattice QCD
        I present our recent analysis of isovector nucleon electromagnetic form factors, as well as the axial, scalar and tensor charges with high statistics in lattice QCD. By applying the all-mode-averaging technique to two flavors of Wilson-Clover fermions, we obtain precise lattice results of form factors at various lattice spacings and pion masses in the range 200--500 MeV for \m_{\pi}L>4. An extensive numerical study of excited state contamination suggests that large effects persist for source-sink separations of less than 1.3 fm when the standard plateau method is used. In order to reduce it to less than a few percent uncertainty, separations of more than 1.5 fm are required. I present the results from several analyses including the first excited state and perform a consistency check with the region where the ground state dominates. I also discuss the chiral behaviour and systematic uncertainties arising from finite size effects and lattice artifacts.
        Speaker: Dr Eigo Shintani (Mainz University/AICS)
        Slides
      • 121
        Nucleon-pion-state contributions in the determination of the nucleon axial charge
        The nucleon-pion-state contributions to QCD two- and three-point functions used in the calculation of the nucleon axial charge are studied in chiral perturbation theory. For sufficiently small quark masses and large volumes the nucleon-pion states are expected to have smaller total energy than the single-particle excited states. To leading order in chiral perturbation theory the results do not depend on low-energy constants associated with the interpolating nucleon fields and apply to local as well as smeared interpolators. The nucleon-pion-state contribution is found to be at the few percent level and non-negligible for precise determinations of the axial charge.
        Speaker: Dr Oliver Baer (Humboldt University Berlin)
        Slides
      • 122
        Some nucleon isovector obsesrvables from 2+1-flavor domain-wall QCD at physical mass
        Current status will be reported of the nucleon isovector form factors and low moments of structure functions calculated by RBC and LHP Collaborations using a physical mass 2+1-flavor dynamical domain-wall fermion (DWF) lattice-QCD ensemble generated by UKQCD and RBC Collaborations at momentum cutoff of 1.730(4) GeV and lattice spatial extent of 5.476(12) fm.
        Speaker: Prof. Shigemi Ohta (KEK)
        Slides
    • Nonzero Temperature and Density 401

      401

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Carleton DeTar
      • 123
        Magnetic properties and deconfinement
        We report on our recent results regarding the properties of strongly interacting matter in the presence of strong magnetic background fields below and above the deconfinement transition.
        Speaker: Mr Marco Mariti (University of Pisa)
        Slides
      • 124
        Towards the continuum limit of the critical endline of finite temperature QCD
        We report on preliminary results for the critical endpoint of finite temperature phase transition of $N_f=3$ QCD at $N_t=10, 12$ and the critical endline of $N_f=2+1$ QCD at $N_t=6$ around the SU(3)-flavour symmetric point. We employ the renormalization-group improved Iwasaki gauge action and non-perturbatively $O(a)$-improved Wilson-clover fermion action.
        Speaker: Dr Yoshifumi Nakamura (RIKEN AICS)
        Slides
      • 125
        Chiral phase transition of $N_f$=3 and 2+1 QCD at vanishing baryon chemical potential
        We present updated studies of chiral phase transition in Nf=2+1 and 3 QCD based on the simulations using Highly Improved Staggered Quarks on lattices with temporal extent $N_\tau$ =6 at vanishing baryon chemical potential. In Nf=2+1 QCD we have performed simulations with a strange quark fixed to its physical value and two degenerate light quarks whose values are adjusted to have 5 values of Goldstone pion masses in the region of 160 - 80 MeV in the continuum limit. Datasets with more values of temperature near chiral phase transition temperature are produced compared to our previous study (1312.0119). The universal scaling behavior of chiral condensates as well as chiral susceptibilities and uncertainties of the determination of non-universal parameters, e.g. chiral phase transition temperature Tc, are discussed. Simulations with 3 degenerate quarks have also been performed with more temperature values and with almost doubled statistics compared to our previous investigation (1302.5740). Signatures of first order phase transition and the estimate of the value of critical quark mass where the second order phase transition ends are discussed.
        Speaker: Prof. Heng-Tong DING (Central China Normal University)
        Slides
      • 126
        Thermodynamics and reference scale of SU(3) gauge theory from gradient flow on fine lattices
        We study thermodynamics of SU(3) gauge theory on the basis of the Yang-Mills gradient flow on fine lattices. For this purpose, the lattice spacing of the Wilson gauge action is determined over a wide range of $beta$ with high accuracy. We then measure the flow time and lattice spacing dependences of the expectation values of the energy-momentum tensor. The extrapolation to the continuum limit of these results is performed.
        Speaker: Prof. Masakiyo Kitazawa (Osaka University)
        Slides
      • 127
        Polyakov loop renormalization with gradient flow
        We propose to use gradient flow for the renormalization of Polyakov loops in various representations. We study Polyakov loops in 2+1 flavor QCD using HISQ action and lattices with temporal extents Nt=6,8,10 and 12 in various representations, including fundamental, sextet, adjoint, decuplet, 15-plet and 27-plet. This alternative renormalization procedure allows the renormalization over a large temperature range from T=100 up to T=1000 MeV, with small errors not only for the fundamental but also for higher representations of the Polyakov loop. We discuss the Casimir scaling of the Polyakov loops and also compare with weak coupling results.
        Speaker: Mr Hans-Peter Schadler (University of Graz)
        Slides
      • 128
        Pure SU(3) Topological Susceptibility at Finite Temperature with the Wilson Flow
        We give a summary of our preliminary results on the finite temperature topological susceptibility $\chi$ from pure SU(3) theory. The simulations use a Symanzik improved action and a gluonic definition of the topological charge with cutoff effects at the $a^2$ level. We use the Wilson flow to calculate a properly renormalized topological charge and its susceptibility. Our results suggest a very strong decay of the topological susceptibility above $T_c$ in line with the results in literature.
        Speaker: Dr Simon Mages (Forschungszentrum Jülich)
        Slides
    • Physics Beyond the Standard Model 405

      405

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Daniel Nogradi
      • 129
        Effects of higher order operators on the Standard Model Higgs sector
        We study generic effects of Beyond the Standard Model physics on the Standard Model Higgs sector. Specifically we investigate the effects of higher order operators in the Higgs potential in the Higgs-Yukawa and gauge-Higgs sectors separately. Our main result is the dependence of the critical endpoint of the Electroweak finite temperature phase transition on the scale of new physics.
        Speaker: Mr Oscar Akerlund (ETH Zurich)
        Slides
      • 130
        Lattice study of the Higgs-Yukawa model in and beyond the Standard Model
        We consider a chirally-invariant Higgs-Yukawa model as a limiting case of the Standard Model (SM) and investigate the model using lattice simulation. We present result of the study of the model near the mean-field using Renormalisation Group (RG) improved finite-size scalling. We derived formulae for such scaling behaviour by including the contribuations of the leading logarithms. Results for the magnetisation, susceptibility and Binder's cumulant are tested against the lattice simulation data performed with small bare scalar-quartic and Yukawa couplings. In these tests we found the predicted logarithmic volume dependence of those observables. The technique we establish in this work can be applied to an extensive investigation of the phase structure of the Higgs-Yukawa models. In particular, it can be used to study the possibility of having non-trivial fixed points in such models. In a different aspect of our work, we have also included a $(\phi^{\dagger}\phi)^3$ term to investigate effects of a higher dimensional operator in this model. Such a higher-dimensional operator can originate from physics beyond the SM. With the zero temperature phase structure established in our previous work, it is worthwhile to study the finite temperature properties of such a model.
        Speaker: David Y. -J. Chu (Department of Electrophysics, National Chaio Tung University)
        Slides
      • 131
        Non-Perturbative Gauge-Higgs Unification in Five Dimensions
        We study the phase diagram and mass spectrum of an $SU(2)$ Gauge-Higgs Unification scenario on a five-dimensional orbifold. This theory exhibits spontaneous symmetry breaking, and we observe that a newly discovered phase transition plays an important role in the ability of the theory to produce a standard model-like spectrum. We comment on dimensional reduction and take first steps towards constructing renormalised trajectories along the phase diagram such that physical quantities remain constant.
        Speaker: Dr Graham Moir (Bergische Universität Wuppertal)
        Slides
      • 132
        Five-dimensional Gauge Theories in a Warped Background
        The anisotropic five-dimensional SU(2) lattice gauge theory in a flat space-time background has a rich phase diagram that has been investigated extensively using various techniques. The objective is to understand whether a layered phase exists that could support a four-dimensional brane in the continuum limit, but the results to date provide no evidence for this. We present new results, obtained in mean-field theory, for the phase diagram of this theory when the extra dimension is warped.
        Speaker: Ms Eliana Lambrou (The University of Edinburgh)
        Slides
      • 133
        Gravitational waves from cosmological first order phase transitions
        Future space-based gravitational wave detectors, e.g. the eLISA satellite constellation scheduled for launch in 2034, can potentially observe primordial gravitational radiation, opening a new window to the early Universe. Gravitational radiation can be generated in first order phase transitions, for example the electroweak transition in some extensions of the Standard Model. In order to obtain detailed understanding of the generation of gravitational radiation, we have initiated a research program using large-scale simulations of an effective order parameter field + fluid systems. We find that the dominant source of radiation are the acoustic waves of the fluid, the sound of the transition. This hitherto unappreciated mechanism gives significantly stronger gravitational radiation background than the earlier estimates have indicated.
        Speaker: Prof. Kari Rummukainen (University of Helsinki)
        Slides
    • Standard Model Parameters and Renormalization 404

      404

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Stefan Sint
      • 134
        Analysis of short distance current correlators using OPE
        We investigate correlation functions of light quark flavor non-singlet bilinear operators at short distances, where non-perturbative effects are encoded in the form of Operator Product Expansion (OPE). Comparing lattice results and continuum perturbation theory, we determine renormalization constants and vacuum expectation values appearing in OPE. We use the lattice data obtained with the Mobius domain-wall fermion at three lattice spacings.
        Speaker: Mr Masaaki Tomii (SOKENDAI)
        Slides
      • 135
        Charmonium current-current correlators with Mobius domain-wall fermion
        We calculate the moments of charmonium current-current correlators with the Mobius domain-wall fermion. Following the method adopted by the HPQCD collaboration, we extract the charm quark mass by matching the results with the corresponding perturbative QCD calculations. We use the recently generated domain-wall fermion ensembles by the JLQCD collaboration at lattice spacings a = 0.083 fm, 0.055 fm and 0.044 fm.
        Speaker: katsumasa Nakayama (Nagoya University)
        Slides
      • 136
        Determining the QCD coupling from lattice vacuum polarization
        The QCD coupling appears in the perturbative expansion of the current-current two-point (vacuum polarization) function. Any lattice calculation of vacuum polarization is plagued by several competing non-perturbative effects at small momenta and by discretization errors at large momenta. We work in an intermediate region and compute vacuum polarization for many off-axis momentum directions on the lattice. Having many momentum directions provides a way to monitor and account for lattice artifacts. Our results are competitive with, and have certain systematic advantages over, the alternate phenomenological determination of the strong coupling from the same light quark vacuum polarization produced by sum rule analyses of hadronic tau decay data.
        Speaker: Randy Lewis (York University)
        Slides
      • 137
        Light quark masses from infrared fixed point
        It is argued that low energy parameters, such as the light quark masses, are related by infrared fixed points of the underlying theory. Implications for the u and d quark masses are discussed.
        Speaker: Prof. Gerrit Schierholz (DESY)
        Slides
      • 138
        Clover fermions in Numerical Stochastic Perturbation Theory
        We present an implementation of clover fermions in Numerical Stochastic Perturbation Theory. Two loop computations are already possible, e.g. we already computed critical mass. We stress that any two loop computation only relies on the known value of c_SW to one loop. A two loop determination of c_SW would be of great interest: we report on the status of our computation of the latter.
        Speaker: Dr Michele Brambilla (INFN Milano Bicocca and Università di Parma)
        Slides
    • Theoretical Developments 406

      406

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Anthony Kennedy
      • 139
        Hamiltonian simulation of lattice gauge theories
        We consider the matrix product state formalism for the simulation of Hamiltonian lattice gauge theories. To this end, we define matrix product states which are manifestly gauge invariant. As an application, we study 1+1 dimensional one flavor quantum electrodynamics, also known as the massive Schwinger model. First, we discuss our results on full quantum non-equilibrium dynamics induced by a quench in the form of a uniform background electric field (i.e. the Schwinger pair creation mechanism) . Furthermore we study the effects of charge screening and confinement in the vacuum for integer and non-integer external charges. Finally, we present some results on finite temperature simulations. References: [1] arXiv:1312.6654 [hep-lat] [2] arXiv:1411.0020 [hep-lat]
        Speaker: Mr Boye Buyens (Ghent University)
        Slides
      • 140
        Real-time simulation of dissipation-driven quantum Systems
        We set up a real-time path integral to study the evolution of quantum systems driven in real-time completely through the coupling of the system with the environment. This can also be interpreted as measurements being performed on the system. For a spin-1/2 system, in particular, when the measurement results are averaged over, the resulting sign problem completely disappears, and the system can be simulated with an efficient cluster algorithm.
        Speaker: Dr Debasish Banerjee (DESY Zeuthen)
        Slides
      • 141
        Lefschetz-thimble path integral for solving the mean-field sign problem
        Lefschetz-thimble method is a recently developing tool for solving the sign problem. We showed that the sign problem appearing in the mean-field approximation can be completely solved by applying this technique. The manifest reality of the physical observables is shown to be ensured in spite of the complexification of the field variables. The result is demonstrated through the heavy-quark model, and we also discuss some implications to the lattice QCD.
        Speaker: Mr Yuya Tanizaki (University of Tokyo)
        Slides
      • 142
        Thermal evolution of the 1-flavour Schwinger model with using Matrix Product States
        The Schwinger model, or 1+1 dimensional QED, offers an interesting object of study, both at zero and non-zero temperature, because of its similarities to QCD. In this talk we present the first full calculation of the temperature dependent chiral condensate of this model in the continuum limit using Matrix Product States (MPS). MPS methods, in general tensor networks, constitute a very promising technique for the non-perturbative study of Hamiltonian quantum systems. In the last few years, they have shown their suitability as ansatzes for ground states and low-lying excitations of lattice gauge theories. We show the feasibility of the approach also for finite temperature.
        Speaker: Dr Hana Saito (DESY)
        Slides
      • 143
        Tensor renormalization group analysis of CP(N-1) model in two dimensions
        We apply the tensor renormalization group (TRG) method to the lattice CP(N-1) model in two dimensions. A tensor network representation of this model is derived for arbitrary N. For N=2, we compute the average energy by using the higher-order TRG. We compare it with the result of the O(3) nonlinear sigma model in two dimensions which is analyzed by the same method. Finally, we discuss the tensor network representation in the presence of the theta term.
        Speaker: Mr Hikaru Kawauchi (Kanazawa University)
        Slides
      • 144
        Approaching conformality with the Tensor Renormalization Group method
        We discuss the Tensor Renormalization Group (TRG) method for various spin and gauge lattice models in 1+1 dimensions with emphasis on near gapless/conformal situations. We present recent progress on optimized truncation methods. We describe the entanglement entropy in the superfluid phase of the O(2) model with a chemical potential.
        Speaker: Prof. Yannick Meurice (University of Iowa)
        Slides
    • Poster Session 4th floor

      4th floor

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      • 145
        A new method to calculate the Dirac operator spectral density
        We describe a new method to determine in a stochastic fashion the Dirac operator spectral density over the entire range of eigenvalues. This allows for example studies of the chiral limit of the fermion condensate, the mode number and the anomalous mass dimension. We apply this technique to the 2-flavor sextet SU(3) gauge theory, to investigate this theory's viability as a minimal realization of the composite Higgs scenario.
        Speaker: Prof. Kieran Holland (University of the Pacific)
        Poster
      • 146
        A novel computation of the thermodynamics of SU(3) Yang-Mills theory
        We present an accurate computation of the Equation of State of SU(3) Yang-Mills theory using shifted boundary conditions in the temporal direction. In this framework, the entropy density s(T) can be obtained in a simple way from the v.e.v. of the space-time components T_0k of the energy-momentum tensor. Furthermore, contrary to the standard approach, the Equation of State s(T)/T^3 can be measured in an independent way at any value of the temperature. The extrapolation to the continuum limit shows almost no dependence on the lattice artifacts.
        Speaker: Prof. Leonardo Giusti (University of Milano Bicocca)
      • 147
        Adaptive algebraic multigrid on SIMD architectures
        We present details of our implementation of the Wuppertal adaptive algebraic multigrid code on SIMD architectures, with particular emphasis on the Intel Xeon Phi processor (KNC) used in QPACE 2. As a smoother, the algorithm uses a domain-decomposition-based solver code previously developed for the KNC in Regensburg. We optimized the remaining parts of the multigrid code and conclude that it is a very good target for SIMD architectures. Some of the remaining bottlenecks can be eliminated by a multiple-right-hand-sides setup, which is discussed in the contribution of Daniel Richtmann.
        Speaker: Prof. Tilo Wettig (University of Regensburg)
      • 148
        An application of the hybrid Monte Carlo algorithm for realized stochastic volatility model
        The hybrid Monte Carlo (HMC) algorithm has been widely used for dynamical lattice QCD simulations. One of the advantages of using the HMC algorithm is that it is a global algorithm that can update all link variables simultaneously. In this way we can greatly reduce computational cost concerning the fermionic part. We utilize this advantage for parameter estimations of the realized stochastic volatility model which has been used for modelling time series data. The realized stochastic volatility model includes a number of volatility variables to be updated. We update those variables by the HMC algorithm. It is found that the HMC algorithm de-correlates effectively Monte Carlo samples of volatility variables. We also show that the algorithm can be accelerated by the GPU computing.
        Speakers: Prof. Tetsuya Takaishi (Hiroshima University of Economics), Prof. Yubin Liu (Nankai University)
        Poster
      • 149
        Beating the sign problem in finite density lattice QCD
        At finite density lattice QCD, the famous sign problem has prevented us from studying the QCD phase diagram based on the first-principle calculation. In order to circumvent the problem, we propose a method where we put new wine into old wineskins: i.e., we use an old idea, "canonical approach", but put into it a new ingredient, multi-precision calculations. This is free from the usual limitation due to Taylor expansion. We briefly present the history of the canonical partition function method, describe our formulation, and show outcomes of the first test, which are very promising and encourages us to go to the final stage, i.e., the study of the QCD under extreme conditions by realistic simulations at small quark masses and on large lattices.
        Speaker: Prof. Atsushi Nakamura (RCNP, Osaka Univ.)
        Poster
      • 150
        Bosonization analysis for artificial "atomic collapse" in graphene
        Around a large charge with atomic number Z > 137, the QED vacuum is expected to collapse due to the strong Coulombic force. While the relativistic quantum mechanics fails to make reliable predictions for the fate of the vacuum, the heavy ion collision experiment also does not give clear understanding of this system. Recently,the "atomic collapse" resonances were observed on graphene where an artificial nuclei can be made. In this poster, I will present our non-perturbative study of the vacuum strucutre which contains multi-body effect using bosonization method.
        Speaker: Ms Aya Kagimura (Particle Physics Theory Group, Department of Physics, Osaka university)
        Poster
      • 151
        calculation of strange and light quark condensate using improved staggered fermions and overlap fermions
        We calculate the strange quark condensate and the light quark condensate by using improved staggered fermions and overlap fermions, respectively. Then we report the mass dependence of the quark condensate. We use N_f = 2+1+1 MILC HISQ(highly improved staggered quark) gauge ensembles. We use a normal CG inverter to calculate the quark condensates for higher quark masses, and the Lanczos and the eigCG algorithms for lower quark masses. We also study the multi-mass algorithm.
        Speaker: Mr Hwancheol Jeong (Seoul National University)
        Poster
      • 152
        Determination of $\varepsilon_K$ using lattice QCD inputs
        We present the Standard Model evaluation of the indirect CP violation parameter $\varepislon_K$ determined using lattice QCD inputs: $\hat{B}_K$, $\xi_0$, $V_{us}$, and $V_{cb}$. We find that the Standard Model prediction of $\varepsilon_K$ with exclusive $V_{cb}$ (lattice QCD results) is lower than the experimental value by $3.6\sigma$. However, this tension disappears with inclusive $V_{cb}$ (results of heavy quark expansion).
        Speaker: Prof. Weonjong Lee (Seoul National University)
        Slides
      • 153
        Determination of f_K/f_pi from staggered Nf=2+1+1 ensembles
        Preliminary results are presented from a new determination of f_K/f_pi which is based on a series of staggered Nf=2+1+1 simulations in the vicinity of the physical point. Special emphasis is put on the interpolation to the physical point, the finite-volume corrections, the continuum extrapolation and the breaking of the isospin symmetry.
        Speaker: Dr Stephan Durr (Wuppertal University / JSC Julich)
        Poster
      • 154
        Exploring free-form smearing for bottomonium and B meson spectroscopy
        Free-form smearing was designed as a way to implement source operators of any desired shape. A variation of the method is introduced that reduces the computational cost by reducing the number of link multiplications to its absolute minimum. Practical utility is demonstrated through calculations of bottomonium and B meson masses.
        Speaker: Richard Woloshyn (TRIUMF)
        Poster
      • 155
        Grassmann tensor renormalization group for the lattice Gross-Neveu model with finite chemical potential
        In numerical calculation of the lattice QCD, the sign problem of Monte Carlo methods has been a major issue. A numerical renormalization group called tensor renormalization group (TRG) is one of the candidates for the solution of this problem. Grassmann tensor renormalization group (GTRG) is a generalization of the TRG to fermion systems. In this study, We apply the GTRG to the lattice Gross-Neveu model in the presence of a chemical potential as a benchmark for future works in finite density systems.
        Speaker: Yusuke Yoshimura (RIKEN Advanced Institute for Computational Science)
      • 156
        Grid: A next generation C++ library for data parallel QCD
        We discuss progress and performance from the development of the Grid data parallel interface for cartesian fields of tensor data types. The library is suitable for the simulation of QCD including multiple grids. Wherever appropriate the interface is quite similar to that of QDP++, but the library makes use of C++11 features to reduce the volume of code, compared to QDP++ while bringing greater generality and greater performance. While the library simultaneously targets MPI, OpenMP and SIMD parallelism, the SIMD optimisation is notably flexible. Performance is substantially improved through the use of a data layout transformation and benchmark results are presented for simple Lattice QCD operations and the Wilson operator under AVX, AVX2 and AVX512 instruction set targets. This work is performed as part of the Intel Parallel Computing Centre in the Higgs Centre for Theoretical Physics at Edinburgh. Submitted by Peter Boyle on 14 May 2015 at 23:48 Id: 5 Last modification: 14 May 2015 23:48 Contribution type: Talk Content We discuss progress and performance from the development of a data parallel interface for cartesian fields of tensor data types. The library is suitable for the simulation of QCD including multiple grids. Wherever appropriate the interface is quite similar to that of QDP++, but the library makes use of C++11 features to reduce the volume of code, compared to QDP++ while bringing greater generality and greater performance. While the library simultaneously targets MPI, OpenMP and SIMD parallelism, the SIMD optimisation is notably flexible. Performance is substantially improved through the use of a data layout transformation and benchmark results are presented for simple Lattice QCD operations and the Wilson operator under AVX, AVX2 and AVX512 instruction set targets.
        Speaker: Dr Azusa Yamaguchi (University of Edinburgh)
        Poster
      • 157
        Heavy-heavy current improvement for calculating B ̅ →D^((*) ) lν ̅ semi-leptonic form factors with Oktay-Kronfeld quarks.
        We improve heavy-heavy currents for calculating B ̅ →D^((*) ) lν ̅ semi-leptonic form factors with Oktay-Kronfeld (OK) heavy quarks. The OK action, which has dimension 6 and 7 interaction terms, can control the discretization errors of heavy quarks (b and c quarks). The OK action is improved through third order in HQET power counting. We report work on heavy-heavy currents to get the systematic improvement of the hadronic matrix elements for B ̅ →D^((*) ) lν ̅ processes with the OK action.
        Speaker: Mr JAEHOON LEEM (Seoul National University)
        Slides
      • 158
        Instantaneous Stochastic Perturbation Theory and Gradient flow in $\phi^4$ theory
        L\"uscher recently introduced Instantaneous Stochastic Perturbation Theory (ISPT) as a method for the stochastic evaluation of lattice perturbation theory. We present an exploratory study of its properties in $\phi^4$ theory, and we compare it with both analytical computations and Numerical Stochastic Perturbation Theory. In addition, we propose a automated method based on ISPT for solution of the Gradient flow equations.
        Speaker: Mr Marco garofalo (The University of Edinburgh)
        Poster
      • 159
        Investigating some technical improvements to glueball calculations.
        Preliminary results for the effect of open boundary conditions on the masses of glueballs and the probability distributions of glueball correlators, in quenched QCD, are presented. Initial results from glueball correlators computed on unquenched ensembles downloaded from the ILDG data grid are reported.
        Speaker: Dr Craig McNeile (Plymouth University)
        Slides
      • 160
        Isospin splitting in Wilson chiral perturbation theory for twisted-mass lattice-QCD with three non-degenerate quark flavours
        We present the expansions of the masses of the pseudoscalar mesons in chiral perturbation theory at next-to-leading order for twisted mass lattice QCD with three light quark flavours, taking the mass difference between the up and down quarks into account.
        Speaker: Prof. Gernot Münster (University of Münster)
        Poster
      • 161
        Lattice QCD code set Bridge++ on arithmetic accelerators
        We are developing a code set "Bridge++" for simulations of lattice gauge theories that aims at an extensible, readable, and portable workbench, while achieving high performance. This work concerns the design of Bridge++ to incorporate the accelerator devices, such as GPUs. As generic frameworks for heterogeneous programming using CPU and accelerator devices, we apply OpenCL and OpenACC to Bridge++. OpenCL and OpenACC are based on different policies. The former explicitly controls the devices through API functions, while the latter is a directive-based extension of a programming language. The inversion of a fermion matrix is offloaded to the accelerator devices. From a viewpoint of constructing reusable components based on the object-oriented programming, as well as tuning the code to achieve high performance, we evaluate feasibility of these frameworks.
        Speaker: Dr Shinji Motoki (High Energy Accelerator Research Organization(KEK))
      • 162
        Lattice QCD study of the I=0 scalar channel using four-quark operators
        We study the possible significance of four-quark states in the isosinglet scalar mesons ($J^{PC}=0^{++}$, I=0) by performing two-flavor full lattice QCD simulations on an $8^3 ¥times 16$ lattice using the improved gauge action and the clover-improved Wilson quark action. In particular, we evaluate the propagators of molecular and tetraquark operators together with singly disconnected diagrams. In the computation of the singly disconnected diagrams we employ the $Z_2$-noise method with the truncated eigenmode approach. We show that the quark loops given by the disconnected diagrams play an essential role in propagators of tetraquark and molecular operators.
        Speaker: Mr Masayuki Wakayama (Nagoya University)
      • 163
        Mass and Axial current renormalization in the Schr\"{o}dinger functional scheme for the RG-improved gauge and the stout smeared $O(a)$-improved Wilson quark actions.
        We present the quark mass and axial current renormalization factors for the RG-improved gauge and the stout smeared $O(a)$-improved Wilson quark actions. The $O(a)$ improvement coefficient for the three-flavors of dynamical quarks has been determined previously with the stout-link smearing parameter $\alpha=0.1, n=6$. We employ the Schr\"{o}dinger functional scheme and obtain the renormalization factors at $\beta=1.82$ with three-flavors of quarks where a large scale simulation is being carried out.
        Speaker: Mr Ken-Ichi Ishikawa (Hiroshima University, Department of Physical Science)
        Poster
      • 164
        Mass anomalous dimension of SU2 with Nf=8 using the spectral density method
        SU2 with Nf=8 is believed to have an infrared conformal fixed point. We use the spectral density method to evaluate the coupling constant dependence of the mass anomalous dimension. We use massless theory with Schrödinger functional boundary conditions and HEX smeared Wilson clover fermions.
        Speaker: Mr Joni Suorsa (University of Helsinki, Helsinki Institute of Physics)
        Poster
      • 165
        Monte Carlo studies of dynamical compactification of extra dimensions in a model of nonperturbative string theory
        The IIB matrix model has been proposed as a nonperturbative definition of superstring theory. In this work, we study the Euclidean version of this model in which extra dimensions can be dynamically compactified if a scenario of the spontaneous breakdown of the SO(10) rotational symmetry is realized. Monte Carlo calculations of the Euclidean IIB matrix model suffer from a very strong complex action problem due to the large fluctuations of the complex phase of the Pfaffian after integrating out the fermion. We employ the factorization method to achieve the effective sampling. We report on preliminary results that can be compared with previous studies of the rotational symmetry breakdown using the Gaussian expansion method.
        Speaker: Prof. Takehiro Azuma (Setsunan University)
        Slides
      • 166
        Multiple right-hand side setup for the DD-$\alpha$AMG
        The setup cost of a modern solver such as the DD-$\alpha$AMG (Wuppertal Multigrid) is a significant contribution to the total time spent on solving the Dirac equation, and in HMC it can even be dominant. We present an improved implementation of this algorithm with modified computation order in the setup procedure. By processing multiple right-hand sides simultaneously we can alleviate many of the performance issues of the default single right-hand side setup. The main improvements are as follows: Many matrix-vector products are replaced by matrix-matrix products, leading to better cache reuse. The synchronization overhead inflicted by on-chip parallelization (threading), which is becoming crucial on many-core architectures such as the Intel Xeon Phi, is effectively reduced. By combining multiple right-hand sides the message size for off-chip communication is larger, which leads to better utilization of the network bandwidth. In the parts implemented so far, we observe a speedup of roughly 2x compared to the optimized version of the single right-hand side setup on realistic lattices.
        Speaker: Mr Daniel Richtmann (University of Regensburg)
      • 167
        Non-perturbative Renormalization with RI-MOM scheme for Bilinear Operators on the Fine Lattice
        We present the preliminary results of the wave function renormalization factor Z_q and mass renormalization factor Z_m from the bilinear operators obtained using non-perturbative renormalization method(NPR) in the RI-MOM scheme with improved staggered fermions. We use fine ensembles of MILC asqtad lattices (N_f=2+1) with 28^3 \times 96 geometry and am_\ell / am_s = 0.0062 / 0.031. We also present the dependence of lattice spacing for Z_q and Z_m by comparing the results of coarse and fine lattices.
        Speaker: Dr jangho Kim (Korea Institute of Science and Technology Information)
        Slides
      • 168
        Overlap fermions on GPUs
        We report on our efforts to impliment overlap fermions on NVIDIA GPUs using CUDA, commenting on the algorithms used, implemetation details, and the performance of our code.
        Speaker: Dr Nigel Cundy (Seoul National University)
        Poster
      • 169
        Performance of Maxwell GPUs and Optimization of Non-Perturbative Renormalization codes.
        Matching factors for the four-fermion operators relevant to B_K are obtained using the Non-Perturbative Renormalization(NPR) method in the RI-MOM scheme. Our NPR codes are calculated in Double Precision(DP) and optimized on Fermi GPU. NVIDIA released The GTX Titan X based on the Maxwell architecture. The Single Precision(SP) calculation performance of the GTX Titan X is improved, but the ratio of DP performance to SP performance is low. Hence, We optimize NPR codes on GTX Titan X and measure its scalability on Multi GPUs. And, we measure the performance of GTX Titan X using our conjugate gradient(CG) solver.
        Speaker: Mr Jeonghwan Pak (Seoul National University)
        Slides
      • 170
        Precision study of critical slowing down in lattice simulations of the CP^{N-1} model
        With the aim of studying the relevance and properties of critical slowing down in Monte Carlo simulations of lattice quantum field theories we carried out a high precision numerical study of the discretised two-dimensional CP^{N-1} model at N=10 using an over-heat bath algorithm. We identify critical slowing down in terms of slowly-evolving topological modes and present evidence that other observables couple to these slow modes. This coupling is found to reduce however as we increase the physical volume in which we simulate.
        Speaker: Dr Andreas Juettner (University of Southampton)
      • 171
        Proposal for the Quantum Simulation of the CP(2) Model on Optical Lattices
        The 2d CP(N-1) models share a number of features with QCD, like asymptotic freedom, a dynamically generated mass gap, and topological sectors. They have been formulated and analysed successfully in the formalism of the so-called D-theory. In that framework, we propose an experimental set-up for the quantum simulation of the CP(2) model. It is based on ultracold alkaline-earth-atoms located on the sites of an optical lattice, where the nuclear spins represent the relevant fields. We present results for the correlation length and for tunnelling transitions, to be compared with such a future experiment. The latter would also enable the exploration of theta vacua and the phase diagram at finite chemical potential, since it does not suffer from the sign problem.
        Speaker: Prof. Wolfgang Bietenholz (Universidad Nacional Autonoma de Mexico)
        Poster
      • 172
        Quark Spin in Proton from Anomalous Ward Indentity
        We report a quark spin calculation from the anomalous Ward identity with overlap fermions on 2+1 flavor dynamical fermion configurations. Such a formulation decomposes the divergence of the flavor-singlet axial-vector current into a quark pseudoscalar term and a triangle anomaly term, flavor by flavor. We use the overlap fermion for the valence and the quark loop so that the renormalization constants Zm and ZP cancel in the pseudoscalar operator 2mP. In addition, the overlap Dirac operator is used to calculate the local topological charge in the anomaly so that there is no renormalization for the anomaly term either.
        Speaker: Prof. Ming Gong (Institute of High Energy Physics, CAS)
        Poster
      • 173
        S-parameter and vector decay constant in QCD with eight fundamental fermions
        SU(3) gauge theory with eight massless fundamental fermions seems to be near the conformal boundary, and is a candidate of walking technicolor. Along the series of study by LatKMI collaboration using HISQ fermions, S-parameter and vector decay constant, which provide important constraints in the model building, are calculated for this theory. Use of various volumes allows a systematic investigation of finite volume effect, where a strong sensitivity on the S-parameter is found.
        Speaker: Yasumichi Aoki (KMI, Nagoya University)
      • 174
        Stochastic calculation of the QCD Dirac operator spectrum with Mobius domain-wall fermion
        We calculate the spectral function of the QCD Dirac operator using the four-dimensional effective operator constructed from the Mobius domain-wall implementation. We utilize the eigenvalue filtering technique combined with the stochastic estimate of the mode number. The spectrum in the entire eigenvalue range is obtained with a single set of measurements. Results on 2+1-flavor ensembles with Mobius domain-wall sea quarks at lattice spacing ~0.08 fm are shown.
        Speaker: Prof. Shoji Hashimoto (KEK)
      • 175
        Strong coupling expansion of the generalized t-V model in one dimension
        We employ a strong coupling expansion - similar to the one used in the lattice field theory studies [1] - to solve the one-dimensional extended t-V model of fermions on a lattice [2]. This model is solved for a range of filling factors, including both commensurate - where a charge density wave is present - and incommensurate densities. The first set consists not only of a trivial case of half filling. The method allows us to trace the transition from a Luttinger liquid phase to a Mott insulating phase and calculate the critical parameter K. This simple yet powerful method is not based on Bethe ansatz and it works for both integrable and non-integrable systems. Furthermore, we investigate how tailoring the interaction can introduce other ordered phases of the system [3]. [1] C.J. Hamer, Phys. Lett. B 82, 75-78 (1979); D.P Crewther, C.J. Hamer, Nucl. Phys. B 170, 353-368 (1980). [2] G. Gómez-Santos, Phys. Rev. Lett. 70, 3780 (1993); R.G. Dias, Phys. Rev. B 62, 7791 (2000). [3] P. Schmitteckert, R. Werner, Phys. Rev. B 69, 195115 (2004); T. Mishra et al., Phys. Rev. B 84, 115135 (2011).
        Speaker: Mr Marcin Szyniszewski (Lancaster University)
        Poster
      • 176
        Study of the conformal phase of the SU(3) gauge theory with domain-wall fermions
        We study the vacuum structure of the SU(3) gauge theory with the flavor number within the conformal region. We perform numerical simulation employing the massless Domain-Wall fermions with the one-step stout smearing. We investigate a conjecture based on the idea of "the conformal field theories with an IR cutoff", which was recently studied using the Wilson fermions. Dealing with the Polyakov loops, the pseudo-scalar correlators and the anomalous mass-dimension, we collect further supporting evidence of this conjecture.
        Speaker: Dr Jun Noaki (KEK)
        Poster
      • 177
        SU(2) gauge theory with domain-wall fermions in fundamental and adjoint representations
        We study the SU(2) gauge theory with dynamical domain-wall fermions in fundamental and adjoint representations. Dynamical simulations with fundamental fermions are performed with numbers of flavors Nf=2,4,6, and 8. Nf dependence of the meson masses and decay constants are investigated. The same analyses are applied also to the adjoint fermions, based on the investigation of the Aoki phase structure of the Wilson fermion operator.
        Speaker: Dr Hideo Matsufuru (High Energy Accelerator Research Organization (KEK))
        Poster
      • 178
        The Nonlinear O(3) Model with Chemical Potential in a Dual Representation
        We study the two dimensional nonlinear O(3) model with a chemical potential coupled to the conserved charge. The complex action problem occurring at nonzero chemical potential is solved by means of a duality transform, which yields real and positive weights and makes Monte Carlo simulations feasible. We find a phase transition at the end of a Silver Blaze region, where a net number of particles appears in the system. We show various results for bulk observables and correlators at finite chemical potential and discuss their finite size behavior. Furthermore, an approximate method to obtain the two particle potential is discussed.
        Speaker: Mr Thomas Kloiber (University of Graz)
      • 179
        The one-loop analysis of the beta-function in the Schroedinger Functional for Moebius Domain Wall Fermions
        Last year we reported the property of the the Moebius Domain Wall Fermions (MDWF) in the Schroedinger Functional (SF) scheme with the palindromic fifth dimensional parameters and the appropriate boundary operator at the tree-level, and the non-universal behavior for a small fifth direction extent $N_5$. We find that this behavior comes from the residual mass at the tree-level and the universality of the MDWF in the SF setup is confirmed after removal of the residual mass even if the value of $N_5$ is small. Moreover we discuss the scaling property of the one-loop contribution to the beta function based on this set up.
        Speaker: Ms Yuko Murakami (Hiroshima University)
        Poster
      • 180
        The Rho Resonance from N_f=2+1+1 Twisted Mass Lattice QCD
        We present first results for I=1 pipi scattering with N_f=2+1+1 twisted mass fermions utilizing the sLapH method. The rho mass and decay width are computed for various pion masses and lattice spacings. Special care was taken in the analysis of systematic uncertainties.
        Speaker: Mr Markus Werner (University of Bonn)
        Poster
      • 181
        The static three-quark potential of various quark configurations
        We study the static three-quark potential in SU(3) lattice gauge theory at zero temperature with the Polyakov loop correlation function consisting of three Polyakov loops. By employing the multi-level algorithm we overcome the smallness of the signal to noise ratio and obtain remarkably clean signals. We extend our previous investigation of the three-quark potential and compute the potential for some eccentric cases such that the three quarks are located at the vertices of obtuse triangles, and are located in line. We compare these results with the cases that the three quarks are located at the vertices of the equilateral, the isosceles, and the right triangles.
        Speaker: Dr Miho Koma (Nihon University, College of International Relations)
        Poster
      • 182
        The step scaling function of the SU(3) 2 flavor sextet model with Wilson fermions
        We investigate the discrete $\beta$ function of the 2-flavor SU(3) sextet model using the gradient flow scheme. Staggered fermion investigations suggest that the system is chirally broken, contradicting earlier Schroedinger functional Wilson fermion studies that were consistent with conformality. Our results, using improved Wilson fermions and the gradient flow RG scheme, suggests a step scaling function that differs significantly from the perturbative 2-loop prediction and in tension with the staggered results. Considering the potential phenomenological impact of this model, it is important to resolve this disagreement.
        Speaker: Dr Anna Hasenfratz (University of Colorado)
        Poster
      • 183
        Towards string breaking with 2+1 dynamical fermions using the stochastic LapH-method
        We investigate the use of stochastically estimated light quark propagators in correlation functions involving a static color source. To this end we compute the static-light meson correlation function in the stochastic LapH framework, using N_f= 2+1 gauge configurations generated through the CLS effort. Together with the static potential, we obtain an estimate for the expected string breaking distance.
        Speaker: Vanessa Koch (Trinity College Dublin)
      • 184
        Towards the QCD equation of state at the physical point using Wilson fermion
        We study the $N_f=2+1$ QCD at nonzero temperatures using nonperbatively improved Wilson quarks of the physical masses by the fixed scale approach. We perform physical point simulations at finite temperatures with the coupling parameters which were adopted by the PACS-CS collaboration in their studies using the reweighting technique. Zero temperature values are obtained on the PACS-CS configurations which are open to the public on the ILDG. Finite temperature configurations are generated with the RHMC algorithm. The lattice sizes are $32^3 x N_t$, that $N_t=14,13,\cdots,4$ correspond to $T \sim 140-500$ MeV. We present some basic observables at these temperatures and the status of our calculation of the equation of state.
        Speaker: Dr Takashi Umeda (Hiroshima University)
        Poster
      • 185
        Update on the Heavy-Meson Spectrum Tests of the Oktay--Kronfeld Action
        We present updated results of the heavy-meson spectrum obtained using the Oktay--Kronfeld (OK) action on MILC asqtad lattices. The OK action was designed to improve the heavy-quark action of the Fermilab formulation by including complete dimension six and part of dimension seven bilinear improvement terms. Improvement terms are truncated by power counting of order $\Lambda^3/m^3$ for heavy-light system and $v^6$ for quarkonium. They suffice for tree-level matching to QCD of the given order from the power counting. To assess the improvement, we update previous numerical tests with heavy meson masses by generating new data that covers both charm and bottom quark mass regions on an ensemble of lattice spacing a = 0.12 fm. We update the analyses of the inconsistency parameter and the hyperfine splittings for the rest and kinetic masses.
        Speaker: Mr Yong-Chull Jang (Seoul National University)
        Slides
    • Chiral Symmetry 404

      404

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Robert Mawhinney
      • 186
        Exotic Quantum Critical Points with Staggered Fermions
        We study two flavors of massless staggered fermions interacting via an on-site four- fermion interaction and argue that the model contains an exotic quantum critical point separating the perturbative massless phase from a massive fermion phase at strong cou- plings where the fermion bilinear condensate remains zero. We believe that no spontaneous symmetry breaking occurs at the transition. We have extensive calculations in three Eu- clidian dimensions that are consistent with the existence of a single second order phase transition separating the two phases. Although mean field theory suggests that this tran- sition will turn first order at sufficiently large number of dimensions, preliminary results suggest that the transition remains second order in four-dimensions.
        Speaker: Mr Venkit Ayyar (Duke University)
        Slides
      • 187
        Distribution of the k-th smallest Dirac operator eigenvalue : an update
        Based on the exact relationship to random matrix theory, we present alternative methods of evaluating the probability distribution of the k-th smallest Dirac eigenvalue in the epsilon regime of QCD and related gauge theories. By employing (1) the Nystrom-type evaluation of Fredholm determinants and Pfaffians and/or (2) the interrelationship between tau functions for random matrix ensembles at beta=2,1,4 sharing the weight function, practical troubles and a technical restriction in our previous work [PRD63, 045012 (2001)] are resolved. Especially, this update enables the computation of individual Dirac eigenvalue distributions for a lattice gauge theory with 4n staggered flavors in the pseudo-real representation.
        Speaker: Prof. Shinsuke Nishigaki (Shimane University)
        Slides
      • 188
        The microscopic Twisted Mass Dirac spectrum and the spectral determination of the LECs of Wilson $\chi$-PT
        We present the comparison of the analytical microscopic spectral density for lattice QCD with $N_f=2$ twisted mass fermions with the one obtained on the lattice utilizing configurations produced by the ETM collaboration. We extract estimates for the chiral condensate as well as the low-energy constant $W_8$ of Wilson $\chi$-PT by employing spectral information of the Wilson Dirac operator with fixed index at finite volume.
        Speaker: Dr Savvas Zafeiropoulos (Goethe University Frankfurt)
        Slides
      • 189
        Chiral Symmetry Breaking for Bosonic Partition Functions.
        The behavior of quenched Dirac spectra of two-dimensional lattice QCD is consistent with spontaneous chiral symmetry breaking which is forbidden according to the Coleman-Mermin-Wagner theorem. One possible resolution of this paradox is that because of the bosonic determinant that enters in the partially quenched partition function the conditions of this theorem are violated allowing for spontaneous symmetry breaking in two dimensions or less. This goes back to work by Niedermaier and Seiler on nonamenable symmetries of the hyperbolic spin chain and work with Splittorf on bosonic partition functions at nonzero chemical potential. In this talk we discuss chiral symmetry breaking for a bosonic random matrix theory at imaginary chemical potential for which the chiral symmetry group is noncompact and compare with results from the corresponding fermionic theory with the usual compact symmetries.
        Speaker: Prof. Jacobus Verbaarschot (Stony Brook University)
        Slides
      • 190
        Boundary effects on the chiral condensate from Lattice QCD
        In the presence of a constant electric field the vacuum becomes unstable due to Schwinger pair creation. One way to avoid this instability is to use Dirichlet boundary conditions. In this case the chiral condensate vanishes on the walls and it is important to know how quickly it gets restored to its bulk value. A sigma-model calculation predicts that the region where the condensate differs significantly from its bulk value is rather thick. We show by direct calculation that this region is much thinner than estimated.
        Speaker: Dr Andrei Alexandru (The George Washington University)
        Slides
    • Hadron Spectroscopy and Interactions 402

      402

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Chris Allton
      • 191
        X(3872) and Y(4140) using diquark-antidiquark operators with lattice QCD
        I will discuss recent lattice study of charmonium-like mesons with $J^{PC}=1^{++}$ and three quark contents $\bar cc\bar du$, $\bar cc(\bar uu + \bar dd)$ and $\bar cc\bar ss$, where the latter two can mix with $\bar cc$. In this quantum channel, the long known exotic candidate, X(3872), resides. This simulation employs $N_f=2$, $m_\pi=266~$MeV and a large basis of $\bar cc$, two-meson and diquark-antidiquark interpolating fields, with diquarks in both anti-triplet and sextet color representations. It aims at the possible signatures of four-quark exotic states. Along the way, I will discuss the relations between the diquark-antidiquark operators and the two-meson operators via the Fierz transformations.
        Speaker: Dr Padmanath Madanagopalan (University of Graz)
        Slides
      • 192
        Mass spectrum of mesons containing charm quarks - continuum limit results from twisted mass fermions
        We present results from an on-going computation of masses of D mesons, D_s mesons and charmonium, including both ground states and several parity and angular momentum excitations. We employ 2+1+1 flavours of dynamical maximally twisted mass fermions at three lattice spacings and three u/d quark masses at each lattice spacing. We consider different combinations of valence quarks, with either identical or opposite signs in front of the twisted mass terms. In the end, our setup allows for a good control of different kinds of systematic effects, in particular the quark mass dependence of the resulting meson masses and cut-off effects. We obtain very good agreement with experiment for the well-established states and some of our results are predictions.
        Speaker: Dr Krzysztof Cichy (Goethe University Frankfurt)
        Slides
      • 193
        Ground state charmed meson spectra for $N_f=2+1+1$
        We present a preliminary study of the charmed meson spectra using the electrically neutral subset of the new Budapest-Marseille-Wuppertal $N_f=2+1+1$ gauge configurations that utilise the 3-HEX smeared clover action. The analysis is performed with a focus on the hyperfine splitting.
        Speaker: Dr Thomas Rae (University of Wuppertal)
        Slides
      • 194
        Exploring possibly existing qq-anti-b-anti-b tetraquark states with qq = ud, ss, cc
        We explore the possible existence of qq-anti-b-anti-b tetraquark states, where qq denotes two up/down, strange or charm quarks. To this end we compute potentials of two static antiquarks in the presence of the two quarks qq of finite mass using lattice QCD. Then we check, whether these potentials are sufficiently attractive to host bound states. We find strong indication for a bound four-quark state for qq=(ud-du)/\sqrt{2}, i.e. isospin I=0. At the same time there is clear evidence against the existence of qq-anti-b-anti-b tetraquarks with strange or charm flavor or isospin I=1.
        Speaker: Ms Antje Peters (Goethe-Universität Frankfurt am Main)
        Slides
      • 195
        Computation of correlation matrices for tetraquark candidates with $J^P = 0^+$ and flavor structure $q_1 \bar{q_2} q_3 \bar{q}_3$
        The conjecture that several recently observed mesons have a structure, which is not dominated by an ordinary quark antiquark pair, but by a four-quark structure, is being actively investigated both theoretical and experimentally. Such a state maybe characterised as molecular state or as a diquark-antidiquark pair. Lattice QCD provides a theoretically sound framework to study such states. In order to quantitatively investigate the internal structure of such mesons one needs the precise computation of correlation matrices with several different interpolating fields that include two- and four-quarks. In this talk we present our study of such correlation matrices paying particular attention to the technical aspects involved in the computation of diagrams that include closed quark loops and disconnected parts. We study tetraquark candidates with $J^P = 0^+$ and flavor structure $q_1 \bar{q_2} q_3 \bar{q}_3$ including e.g.\ the $a_0(980)$ meson, the $D_{s0}^\ast$ meson and some of the the charged $c \bar{c}$ $X$ states. Numerical results on the $a_0(980)$ meson are discussed in detail.
        Speaker: Mr Joshua Berlin (University of Frankfurt)
        Slides
    • Hadron Structure 403

      403

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Shigemi Ohta
      • 196
        Applications of the Feynman–Hellmann theorem in hadron structure
        By considering an appropriate modification of the action, the Feynman–Hellman (FH) theorem provides an alternative technique for studying hadron matrix elements. We report on new developments in the extension of the FH technique to the study of non-forward matrix elements, where preliminary results demonstrate the feasibility of accessing relatively large momentum transfers. In the near future, these methods offer the potential to address physics topics such as the transition to the perturbative regime in the pion form factor and the possible zero crossing in the proton’s electric form factor. We also present updated results on the extraction of the quark spin fractions in a range of hadrons, including contributions from disconnected operator insertions. At the quark masses studied, the results suggest a negative contribution to the total quark spin of the nucleon from disconnected insertions.
        Speaker: Dr Ross Young (University of Adelaide)
        Slides
      • 197
        Gluon saturation and gluon densities
        In the talk I will overview status of QCD at region of phase space when the partons form a dens system leading eventually to parton saturation i.e. situation when the partons are not growing any more due to nonlinear dynamics. I will overview recent developments leading to construction of gluon densities which take into account saturation effects. Finally I will present predictions for observables which eventually will test saturation phenomenon in p-Pb collisions at LHC
        Speaker: Dr Krzysztof Kutak (Instytut Fizyki Jadrowej Polskiej Akademii Nauk)
        Slides
      • 198
        Nucleon-Sigma-Terms from Lattice QCD
        The nucleon-sigma-terms are measures of the light-quark- and strange-quark-content of the nucleon. Especially the stangeness-content is of significant interest for dark-matter searches, as it determines the coupling of several dark matter candidates to hadronic matter. While the sigma-terms can not be measured directly they can be determined via lattice QCD from first principles. The sigma-terms are related to the light- and strange-quark mass dependence of the nucleon mass by the Feynman-Hellmann-theorem. To measure this dependence we used $N_f=1+1+1+1$ ensembles generated with tree-level improved Symmanzik gauge action and tree-level improved clover Wilson fermions with three levels of HEX smearing at four values of the lattice spacing.
        Speaker: Mr Lukas Varnhorst (Wuppertal University)
        Slides
      • 199
        A lattice study of the nucleon quark content at the physical point
        After a detailed analysis of possible sources of systematic uncertainty, ab-initio $N_f=2+1$ results for the up-down and strange quark content - with pion masses all the way down to the physical point - will be presented and discussed.
        Speaker: Dr Christian Torrero (on behalf of the Budapest-Marseille-Wuppertal Collaboration (Centre de Physique Théorique and CNRS, Aix-Marseille Université, France))
        Slides
      • 200
        Proton spin decomposition with overlap fermion
        We calculate glue spin with the overlap fermion on 2+1 flavor domain wall fermion configurations as a function of the nucleon momentum. The calculation is carried out on 24^3 x 64 and 32^3 x 64 lattice at two lattice spacings and the light sea quarks correspond to pion mass at 330 MeV and 290 MeV on these two lattices. Other issues related to the proton spin will be discussed.
        Speaker: Dr yibo Yang (University of kentucky)
        Slides
    • Nonzero Temperature and Density 401

      401

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Akira Ohnishi
      • 201
        Study of high density phase transition in lattice QCD with canonical approach
        The canonical partition function is related to the grand canonical one through the fugacity expansion. In this talk we perform the fugacity expansion by a method of the hopping parameter expansion in temporal direction. The canonical partition function is evaluated for Nf=2 QCD upto baryon numbers of nB=30 in a wide range of temperature. After derivation of the canonical partition function we study the chemical potential dependences of hadronic observables like chiral condensate, quark number density and its susceptibility. In this talk we report a phase transition found at real chemical potential and its dependence on the quark mass and the volume.
        Speaker: Prof. Yusuke Taniguchi (University of Tsukuba)
        Slides
      • 202
        Exploring finite density QCD phase transition with canonical approach - Power of multiple precision computation -
        The canonical approach is one of the methods that have been pursued to solve the sign problem. However, this approach has the numerical instability, and most people have thought that this approach does not work for a long time. We found that this instability can be solved by using the multiple precision computation. In this talk, we will study effects of the multiple precision computation and report some physical quantities in the finite density.
        Speaker: Mr Shotaro Oka (Rikkyo Univ.)
        Slides
      • 203
        Validity range of canonical approach to finite density QCD
        We calculate the baryon chemical potential (μ) dependence of thermodynamic observables, i.e., pressure, baryon number density and susceptibility by lattice QCD using the canonical approach with winding number expansion.We compare the results with those obtained by the multi parameter reweighting (MPR) method; Both methods give very consistent values in the regions where errors of the MPR are under control.The canonical method gives reliable results up to μ/T=3.
        Speaker: Mr Ryutaro Fukuda (The University of Tokyo / ETH Zurich)
        Slides
      • 204
        Calculation of high-order cumulants with canonical ensemble method in lattice QCD
        High-order cumulants are numbers characterizing the probability distribution and have a lot of physical information. However, sign problem makes it difficult to numerical calculation of high-order cumulants in finite density lattice QCD. In this study we realize the calculation of high-order cumulants with "canonical ensemble method" in heavy quark region. Also, we study a finite density phase transition from the specific behavior of high-order cumulants.
        Speaker: Mr Asobu Suzuki (University of Tsukuba)
        Slides
      • 205
        Aspects of topological actions on the lattice
        We discuss properties of actions which are invariant under smooth deformations of the lattice fields, which are restricted to a subspace of phase space. In particular we comment on the continuum limit of such actions and show how to directly measure the free energy. We also investigate possible applications to finite density systems.
        Speaker: Philippe de Forcrand (ETH Zurich &amp; CERN)
        Slides
    • Physics Beyond the Standard Model 405

      405

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Hiroshi Ohki
      • 206
        Template Composite Dark Matter : SU(2) gauge theory with Nf=2
        We present a non perturbative study of SU(2) gauge theory with two fundamental Dirac flavours. We discuss how the model can be used as a template for composite Dark Matter (DM). We estimate two types of interactions of the DM candidate with the Standard Model : the interaction through Higgs exchange and the interaction through photon exchange computing the electric polarizability of the DM candidate. Finally, we briefly discuss the viability of the model given the present experimental constraints.
        Speaker: Dr Vincent Drach (CP3 Origins - University of Southern Denmark)
        Slides
      • 207
        Stealth Dark Matter on the lattice
        The Stealth Dark Matter model contains a candidate composite dark matter particle which appears as a bosonic neutral baryon of a new strongly coupled sector. The elementary constituents of this composite state carry electroweak charges. This construction provides a mechanism to naturally reduce the strength of dark matter interactions with standard model particles. However such interactions exist and can allow direct detection, and collider experiments to put constraints on the model. In order to get predictions from this strongly-coupled model, lattice simulations are employed and give definite predictions for the cross-section of the dark matter candidate with standard nuclei in detectors for example. Lattice simulations are shown to be of utmost relevance for some dark matter models and they are increasingly important to obtain phenomenological predictions of experimental interest.
        Speaker: Dr Enrico Rinaldi (Lawrence Livermore National Laboratory)
        Slides
      • 208
        Lattice QCD and Axion Cosmology
        The Strong CP Problem can be resolved by introducing an additional global symmetry known as Peccei-Quinn symmetry. Once PQ symmetry is broken the associated particle, the QCD axion, is a plausible dark matter candidate. Calculating the cosmological energy density of the axion requires nonperturbative QCD input---the high-temperature topological susceptibility. I will show results from a pure-glue calculation and examine the implications for the axion mass and coupling.
        Speaker: Dr Evan Berkowitz (Lawrence Livermore National Laboratory)
        Slides
      • 209
        Beyond the Standard model matrix elements with the gradient flow
        The nucleon EDM receives potential contributions from fermionic and gluonic local operators representing the effective interactions from BSM physics. We have recently proposed to use the gradient flow to calculate the QCD component of these key BSM matrix elements. I present an update of the results we have obtained.
        Speaker: Andrea Shindler (Forschungszentrum Jülich)
        Slides
      • 210
        Neutron EDM from quark chromoEDM
        The presentation will discuss the connection between RI-SMOM and MSbar schemes for the quark chromo EDM operators. A preliminary study of the signal in the matrix element of neutron EDM using clover quarks on an HISQ ensemble will also be described.
        Speaker: Dr Tanmoy Bhattacharya (Los Alamos National Laboratory)
        Slides
    • Weak Decays and Matrix Elements 406

      406

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Weonjong Lee
      • 211
        Neutral B meson mixings and B meson decay constants in the infinite b quark mass limit with domain-wall light quarks
        Neutral B meson mixing matrix elements and B meson decay constants are calculated in the infinite b quark mass limit. For the light quarks, domain-wall fermion (DWF) formalism is employed. The calculations are carried out using all-mode-averaging (AMA) technique on Nf=2+1 dynamical DWF ensembles with lattice spacings of nearly 0.08fm and 0.11fm, where pion masses are relatively heavy (>280MeV). We show statistical improvements from our previous results. We also present preliminary results at physical pion mass using 48^3 dynamical Mobius DWF ensemble mainly for B_s quantities.
        Speaker: Tomomi Ishikawa (RIKEN BNL Research Center)
        Slides
      • 212
        Decay constants $f_B$ and $f_{B_s}$ from HISQ simulations
        We present preliminary, high-precision results for the decay constants $f_B$ and $f_{B_s}$ from simulations with HISQ heavy and light valence and sea quarks. Calculations are carried out with several heavy valence-quark masses on ensembles with 2+1+1 flavors of HISQ sea quarks at five lattice spacings and three light sea quark mass ratios $m_{ud}/m_s$, including approximately physical sea quark masses. This range of parameters provides excellent control of the continuum limit and of heavy-quark discretization errors.
        Speaker: Prof. Carleton DeTar (University of Utah)
        Slides
      • 213
        Neutral $B$-meson and $D$-meson mixing bag parameters from $2+1$ flavor lattice QCD
        We compute hadronic matrix elements for neutral $B$-meson mixing in the Standard Model and beyond as well as the corresponding operators for the the short-distance contribution to $D$-meson mixing. On the same sets of gauge ensembles, we calculate the leptonic decay constants for the $B$- and $D$-meson systems. The matrix elements are computed with asqtad light quarks and clover heavy quarks in the Fermilab interpretation. We use the MILC asqtad gauge ensembles at four lattice spacings between $0.11$ fm and $0.043$ fm. We review the the matrix element calculations and report on the determination of the bag parameters. The bag parameters are computed by combining the mixing matrix elements and the decay constant results in a bootstrap process which accounts for correlations from both statistics and the dominant systematic errors. Some systematic effects, such as discretization errors, and quark mass tunings, tend to cancel in ratios.
        Speaker: Dr James Simone (Fermilab)
        Slides
      • 214
        Semileptonic $B$-meson decay phenomenology with lattice QCD
        We present Standard Model predictions for phenomenologically interesting observables for the rare decays $B \to K \ell^+\ell^-$, $B \to \pi \ell^+\ell^-$, $B \to K \nu \bar{\nu}$, and $B \to \pi \nu \bar{\nu}$, as well as for $B \to\pi\tau\nu$. All of these processes are sensitive to new physics effects, and there are a number of tensions between experimental measurements and Standard Model expectations of these and similar processes. We recently completed lattice calculations of the form factors for the semileptonic $B\to\pi$ and $B\to K$ transitions. Here we use these form factors to explore the phenomenology of these decays with quantitative control over the theoretical uncertainties.
        Speaker: ran zhou (Fermilab)
        Slides
      • 215
        Phenomenology with Lattice NRQCD b Quarks.
        Lattice NRQCD is a relatively inexpensive formalism that can be used to simulate b quarks on the lattice at their physical mass. The HPQCD collaboration has used radiatively-improved NRQCD for b quarks in bottomonium to determine the decay rate to leptons of the ground-sate Upsilon and its first radial excitation in lattice QCD. Using time moments of vector bottomonium current-current correlators, we present our determination of the b quark mass in the MS-bar scheme. We use NRQCD b quarks and HISQ light quarks -- with masses down to their physical values -- in the semileptonic B to pi decay at zero recoil to calculate f^0(q^2_max). Our results are consistent with the soft pion theorem for this decay, which relates the meson decay constants to the form factor by f^0(q^2_max)=f_B/f_pi in the massless pion limit.
        Speaker: Dr Brian Colquhoun (University of Glasgow)
        Slides
    • 10:10
      Coffee Break
    • Algorithms and Machines 405

      405

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Georg Bergner
      • 216
        A Multigrid Based Eigensolver for the Hermitian Wilson Dirac Operator
        In this talk, we present a multigrid based eigensolver for computing low modes of the hermitian Wilson Dirac operator. For the non-hermitian case, multigrid methods have already replaced conventional Krylov subspace solvers in many lattice QCD computations. Since the gamma5 preserving aggregation based interpolation is valid for both, the hermitian and the non-hermitian case, inversions of very ill-conditioned shifted systems with the non-hermitian operator become feasible. This enables the use of multigrid within shift-and-invert type eigensolvers. We show numerical results from our MPI-C implementation of a Rayleigh quotient iteration with multigrid. For state-of-the-art lattice sizes and moderate numbers of desired low modes we achieve speed ups of an order of magnitude and more over PARPACK.
        Speaker: Mr Matthias Rottmann (University of Wuppertal)
        Slides
      • 217
        Multigrid-accelerated Low-Mode Averaging
        Disconnected contributions to hadronic quantities are noisy and still computationally challenging. Here, we explore the possibilities of Multigrid-based eigensolvers for Low-mode averaging and related methods. Using only the few lowest approximate eigenmodes of the Hermitian Dirac operator $\gamma_5 M$ allows us to avoid expensive exact solves and improve the signal efficiently. In this first test we apply the method to pion and eta correlators for $N_f=2$ Wilson-Clover fermions at several quark masses, down to $m_\pi = 150 MeV$ and volumes of up to $64^4$ sites.
        Speaker: Mr Jakob Simeth (University of Regensburg)
        Slides
      • 218
        Eigenspectrum calculation of the non-Hermitian O(a)-improved Wilson-Dirac operator using the Sakurai-Sugiura method
        We are developing a computer code for calculating eigenvalues of the non-Hermitian O(a)-improved Wilson-Dirac operator. We introduce here the Sakurai-Sugiura method, which is an eigensolver algorithm based on a contour integral, allowing us to calculate desired eigenvalues located inside a given contour. We report the test results for low-lying eigenvalues obtained with free-case, quenched and full QCD configurations up to a lattice size of 96^4.
        Speaker: Dr Hiroya Suno (RIKEN Advanced Institute for Computational Science; RIKEN Nishina Center)
        Slides
      • 219
        Accelerating deflation of eigenvalues for fermion matrix inversions on GPUs
        We present a brief overview of deflation algorithms recently implemented in the SciDAC QUDA library for solving the lattice-Dirac equation with multiple right-hand sides on NVIDIA accelerators. In particular, we discuss implementation aspects of the mixed precision technique for the eigenvalue deflation which helps to considerably relax GPU memory requirements while allowing for an appropriate deflation quality in fermion matrix inverters. We analyze the efficiency of mixed precision deflation on the examples of HISQ and Wilson twisted mass ensembles.
        Speaker: Dr Alexei Strelchenko (FNAL)
        Slides
    • Chiral Symmetry 404

      404

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Andrei Alexandru
      • 220
        NLO and NNLO Low Energy Constants for SU(2) Chiral Perturbation Theory
        We have performed global fits of $f_\pi$ and $m_\pi$, from a variety of RBC-UKQCD domain wall fermion ensembles, to $SU(2)$ partially quenched chiral perturbation theory at NNLO. We report values for 9 NLO and 8 linearly independent combinations of NNLO partially quenched low energy constants, which we compare to other lattice and phenomenological determinations. We discuss the convergence of the expansion and use our large set of low energy constants to make predictions for the pion mass splitting due to QCD isospin breaking effects and the S-wave $\pi \pi$ scattering lengths. Finally, we discuss the effects of including lattice results for the $I=2$ $\pi \pi$ scattering length in our chiral fits.
        Speaker: Mr David Murphy (Columbia University)
        Slides
      • 221
        NLO and NNLO Low Energy Constants for SU(3) Chiral Perturbation Theory
        We have performed global fits of pseudoscalar masses and decay constants, from a variety of RBC-UKQCD domain wall fermion ensembles, to $SU(3)$ partially quenched chiral perturbation theory at NNLO. We report values for 9 NLO and 10 linearly independent combinations of NNLO partially quenched low energy constants, which we compare to other lattice and phenomenological determinations. We discuss the convergence of the expansion and use our large set of low energy constants to make predictions for mass and decay constant splittings due to QCD isospin breaking effects. We also discuss S-wave $K \pi$ scattering lengths.
        Speaker: Robert Mawhinney (Columbia University)
        Slides
      • 222
        Evidence for a new SU(4) symmetry with J=2 mesons
        Recently, a new $SU(4)$ symmetry has been established, which appears after removing the quasi-zero modes from the Overlap Dirac operator. Namely, the $\rho, \rho', \omega, \omega', a_1, b_1$ and possibly $f_1$ $J=1$ mesons become mass degenerate after the quasi-zero mode removal. The $SU(4)$ symmetry contains $SU(2)_L \times SU(2)_R \times U(1)_A$ as a subgroup and predicts, among other things, that all isovector states of given spin $J$ become mass degenerate. Here we study isovector $J=2$ mesons and observe the degeneracy of all states after removing the quasi-zero modes. This result gives evidence that the $SU(4)$ symmetry is present for mesons of total spin $J \geq 1$.
        Speaker: Dr Markus Pak (University of Graz)
        Slides
      • 223
        Testing the Witten-Veneziano Formula on the Lattice
        We employ lattice techniques to compute both sides of the Witten-Veneziano formula. For the one side we perform dedicated quenched simulations and use the spectral projector method to determine the topological susceptibility in pure Yang-Mills theory. The other side we determine in lattice QCD with $N_f=2+1+1$ dynamical Wilson twisted mass fermions including for the first time also the flavour singlet decay constant. Taking the continuum and the SU$(2)$ chiral limits we compare both sides and find good agreement within uncertainties.
        Speaker: Dr Konstantin Ottnad (University of Bonn)
        Slides
    • Hadron Structure 403

      403

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Ross Young
      • 224
        Nucleon charges, form-factors and neutron EDM
        Results for the matrix elements of quark bilinear operators within a nucleon state will be presented using nine ensembles of gauge configurations generated by the MILC collaboration using the HISQ action with 2+1+1 dynamical flavors. Excited state contamination and extrapolation to the physical point will be discussed for the charges g_A, g_S and g_T, the vector and axial form factors, and the quark EDM operator. The impact of these results in constraining scalar and tensor interactions at the TeV scale and on new sources of CP violation will be discussed.
        Speaker: Dr Rajan Gupta (Los Alamos National Lab)
        Slides
      • 225
        The electric dipole moment of the neutron from N_f=2+1+1 twisted mass fermions
        We extract the neutron electric dipole moment (nEDM) |d_n| on configurations produced with N_f=2+1+1 twisted mass fermions with lattice spacing of a=0.082fm and a light quark mass that corresponds to M_{\pi}=370 MeV. We do so by evaluating the CP-odd form factor F_3 for small values of the vacuum angle \theta at the limit of zero momentum transfer. This limit is taken using a parametrization of the momentum dependence and performing a fit as well as using the position space methods we refer to as "continuum derivative" and "direct computation". The extraction of the CP-odd form factor F_3 requires the evaluation of the field theoretical topological charge. We measure the topological charge via cooling and gradient flow using the ordinary Wilson, Symanzik tree-level improved and Iwasaki actions for smoothing. We obtain consistent results for all choices of smoothing actions, smoothing procedures and momentum dependence treating techniques. We report an nEDM of |d_n|/ \theta = -0.036(11)(7) e fm.
        Speaker: Dr Andreas Athenodorou (University of Cyprus)
        Slides
      • 226
        Neutron-antineutron oscillation matrix elements with domain wall fermions at the physical point
        Neutron-antineutron oscillations are hypothetical processes changing the baryon number by 2 units. Compared to proton decay, they present another scenario of baryon number violation and, if found, they would have different implications for phenomenology of baryogenesis. If such baryon number violation exists at higher scales beyond the Standard model, it will manifest itself at the hadron scale as effective six-(anti)quark operators turning neutrons into antineutrons and vice versa. Nucleon matrix elements of these operators are important for experiments looking for evidence of such processes. I will present preliminary results for these matrix elements computed with physical Nf=2+1 domain wall quarks. Results are non-perturbatively renormalized and converted to MSbar.
        Speaker: Dr Sergey Syritsyn (RIKEN BNL Research Center)
        Slides
      • 227
        Electromagnetic structure of charmed baryons - extended to spin-3/2
        The electromagnetic form factors of the spin-3/2 $\Omega$ baryons, namely the $\Omega$, $\Omega_c^\ast$, $\Omega_{cc}^\ast$ and $\Omega_{ccc}$, are calculated in full QCD on $32^3$ x 64 PACS-CS lattices with a pion mass of 156(9) MeV. The electric charge radii and magnetic moments from the E0 and M1 multipole form factors are extracted. Results for the electric quadrupole form factor, E2, are also provided. Quark sector contributions are computed individually for each observable and then combined to obtain the baryon properties. Charm quark contributions are observed to be systematically smaller compared to the strange quark contributions in case of the charge radii and magnetic moments. E2 moments of the $\Omega_{cc}^\ast$ and $\Omega_{ccc}$ are estimated significantly enough to show that their electric charge distributions are deformed. Properties of the spin-1/2 counterparts of the $\Omega_c^\ast$ and $\Omega_{cc}^\ast$ baryons are also computed and a comparison is presented.
        Speaker: Mr Utku Can (Tokyo Institute of Technology)
        Slides
    • Nonzero Temperature and Density 401

      401

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Heng-Tong DING
      • 228
        Diagrammatic Monte Carlo simulations of staggered fermions at finite coupling
        The diagrammatic Monte Carlo method has been a very fruitful tool for taming, and in some cases even solving, the sign problem in several lattice models. We have recently proposed a diagrammatic model for simulating lattice gauge theories with staggered fermions at arbitrary coupling, which extends earlier successful efforts to simulate lattice QCD at finite baryon density in the strong-coupling regime. Here we present the first numerical simulations of our model, using a worm algorithm.
        Speaker: Helvio Vairinhos (ETH Zurich)
        Slides
      • 229
        Lattice simulation of the SU(2)-chiral model at zero and non-zero pion density
        We present a flux representation based lattice formulation of the partition function corresponding to the SU(2)-chiral Lagrangian, including a chemical potential and scalar/pseudo-scalar source terms. Lattice simulations are then used to obtain non-perturbative properties of the theory, in particular its mass spectrum at zero and non-zero pion density.
        Speaker: Mr Tobias Rindlisbacher (ETH Zurich)
        Slides
      • 230
        Heavy and dense QCD from an effective lattice theory
        A three-dimensional effective lattice theory of Polyakov loops can be construted by integrating over the spatial links in combined strong coupling and hopping parameter expansions. The effective theory is valid for sufficiently heavy quarks and coarse to intermediate lattice spacings. In its domain of validity, the sign problem is weak enough to enable simulations at finite baryon chemical potential and allows for a description of the onset transition as well as cold nuclear matter directly from QCD. We summarize recent results and test the analytic effective couplings against improved numerically determined versions. We also compare complex Langevin and Monte Carlo simulations of the effective theory.
        Speaker: Owe Philipsen (Goethe-University Frankfurt)
        Slides
      • 231
        Analytic computations of an effective lattice theory for heavy quarks
        We discuss the fermionic contribution to a three dimensional effective lattice theory for heavy quarks. We first study the mathematical structure of higher order contributions of the hopping parameter expansion to the effective action. We then examine how to apply a graphical linked cluster expansion to retrieve analytical results for various lattice quantities such as nucleon number and nucleon binging energy. The free parameters of the analytic results are temperature, chemical potential and the hopping parameter. The results are compared to Monte Carlo simulations of the dimensionally reduced theory.
        Speaker: Mr Jonas Rylund Glesaaen (Goethe University Frankfurt)
        Slides
    • Theoretical Developments 402

      402

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Tatsuhiro Misumi
      • 232
        Pion spectrum for the 2-flavor staggered Wilson fermion
        The 4 flavors of the usual staggered fermion give rise to 16 pions (pseudoscalar mesons). We study the effect on the pion spectrum of turning on the 2-flavor staggered Wilson term in the staggered fermion action. 8 of the pions become heavy while 8 remain light. 6 of the light pions are identified as two copies of the 2-flavor pion triplet. They have the same light quark content but different heavy (doubler) quark content. The pions within each triplet are degenerate. The remaining two light mesons are identified as two copies of the flavor-singlet eta meson, which remains light in our study since the disconnected piece of its propagator is omitted. This verifies that the 2-flavor staggered Wilson formulation works as intended as far as the pion spectrum is concerned.
        Speaker: Mr Christian Zielinski (Nanyang Technological University)
        Slides
      • 233
        Flavor Filtered Fermions
        We present a new discretization of the lattice Dirac operator that preserves chiral symmetry while suppressing the effects of doubler states. The new operator is constructed from the naive Dirac operator by filtering out doubler states while preserving one physical quark state. In the simplest version the filtering is only approximate and the filter would require fine tuning to produce a single flavor in the continuum limit. We will present results from quenched simulations for the doubler suppression, locality and topology of the operator. We also discuss prospects for improving the filter to reduce and possibly eliminate the need for fine tuning.
        Speaker: James Osborn (ANL)
        Slides
      • 234
        Charged particles in QED with C* boundary conditions I
        In order to calculate QED corrections to hadronic quantities by means of lattice simulations, a coherent description of electrically-charged states in finite volume is needed. In the usual periodic setup, Gauss' law and large gauge transformations forbid the propagation of electrically-charged states. A possible solution to this problem, which does not violate the axioms of local quantum field theory, has been proposed by Wiese and Polley, and is based on the use of C* boundary conditions. We discuss the properties and symmetries of QED in isolation and QED coupled to QCD, with C* boundary conditions. We show that a certain class of electrically-charged states can be constructed in this setup in a fully consistent fashion, without relying on gauge fixing. This class of states covers most of the interesting phenomenological applications. We also calculate finite-volume corrections to the mass of stable charged particles and show that these are much smaller than in non-local formulations of QED. This is the first of two consecutive talks on the subject.
        Speaker: Dr Agostino Patella (CERN and Plymouth University)
        Slides
      • 235
        Charged particles in QED with C* boundary conditions II
        In order to calculate QED corrections to hadronic quantities by means of lattice simulations, a coherent description of electrically-charged states in finite volume is needed. In the usual periodic setup, Gauss' law and large gauge transformations forbid the propagation of electrically-charged states. A possible solution to this problem, which does not violate the axioms of local quantum field theory, has been proposed by Wiese and Polley, and is based on the use of C* boundary conditions. We discuss the properties and symmetries of QED in isolation and QED coupled to QCD, with C* boundary conditions. We show that a certain class of electrically-charged states can be constructed in this setup in a fully consistent fashion, without relying on gauge fixing. This class of states covers most of the interesting phenomenological applications. We also calculate finite-volume corrections to the mass of stable charged particles and show that these are much smaller than in non-local formulations of QED. This is the second of two consecutive talks on the subject.
        Speaker: Dr Nazario Tantalo (Rome University and INFN "Tor Vergata")
        Slides
    • Weak Decays and Matrix Elements 406

      406

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Hiroshi Ohno
      • 236
        Charm Physics at the physical point
        This talk will provide an overview of RBC/UKQCD's charm project on their 2+1 flavour physical point ensembles using Moebius Domain Wall Fermions for the light as well as for the charm quark. I will discuss the analysis strategy and present first preliminary results for D and D_s decay constants.
        Speaker: Mr Justus Tobias Tsang (University of Southampton, UK)
        Slides
      • 237
        Neutral D-Meson Mixing near the Charm Mass
        We will report on the RBC/UKQCD charm project with (2+1)-flavor Domain Wall fermions at the physical point, and present progress on the determination of the matrix element of the VV+AA operator, and the Bag and Xi-parameters in the region of the physical charm quark mass. We will address the region of practical use of the domain wall action for heavy quark simulation and outline a strategy for the renormalization, continuum extrapolation and possibly extrapolation in heavy quark mass.
        Speaker: Ms Ava Khamseh (Edinburgh University)
        Slides
      • 238
        D meson semileptonic decays from lattice QCD with chiral fermions
        We report on our study of D meson semileptonic decays in QCD with 2+1 flavors of the Moebius domain-wall fermions. Gauge ensembles are generated at lattice cut-off around 2.4 GeV and pion masses as low as 300 MeV. Relevant form factors are calculated by using the Moebius action also for charm quarks.
        Speaker: Mr Takashi Suzuki (Osaka university)
        Slides
      • 239
        D meson semileptonic form factors at zero momentum transfer in 2+1+1 flavor lattice QCD
        We present a calculation of the $D\to\K,l,\nu$ and $D\to\pi,l,\nu$ semileptonic form factors at $q^2=0$. These form factors are needed for the determination of the CKM matrix elements $\lvert{V_{cs}}\rvert$ and $\lvert{V_{cd}}\rvert$ respectively. This work uses the HISQ action for both valence quarks and sea quarks on MILC $N_f=2+1+1$ configurations, with ensembles including several at the physical pion mass and lattice spacings ranging from $0.12 fm$ to $0.042 fm$. The calculation employs twisted boundary conditions to calculate the form factors at zero momentum transfer directly. We use HM$\chi$PT to obtain preliminary results at the physical point and in the continuum limit.
        Speaker: Dr Thomas Primer (University of Arizona)
        Slides
    • 12:00
      Lunch Break
    • 13:30
      Excursion
    • Plenary Session BF & 1F Main Hall

      BF & 1F Main Hall

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Hiroshi Suzuki
      • 240
        What is QFT? Resurgent trans-series, Lefschetz thimbles, and new exact saddles
        Resurgent trans-series provide a novel mathematical formalism to unify perturbative and non-perturbative physics, leading to new insights into the general structure of quantum theories. I will review the main ideas with some illustrative examples. In many quantum mechanical systems, all orders non-perturbative data is encoded into perturbation theory, and it can be decoded. In QFTs, there are cases in which resurgence provides a new interpretation of IR-renormalon puzzle, reveals the existence of many new saddles (such as magnetic and neutral bions), and potentially provide a non-perturbative continuum definition of QFT in a semi-classical domain. I will also describe a new perspective on path integration, which is intimately tied with resurgence theory and employs some tools of Picard-Lefschetz theory. This perspective leads to many dramatic and surprising results, and implies that the proper framework to study semi-classics in path integral formulation is yet to be developed.
        Speaker: Prof. Mithat Unsal (North Carolina State University)
        Slides
      • 241
        Lattice's Bright, Dark and Safe side
        Circa 95% of the universe is made of unknown forms of matter and energy, while to describe the remaining 5% one needs at least three fundamental forces, i.e. Quantum Electrodynamics, Weak Interactions and Quantum Chromo Dynamics (QCD). Furthermore strong interactions are responsible for creating the bulk of the bright mass, i.e. the 5%. It is therefore natural to expect that to correctly describe the rest of our universe while providing a sensible link to the visible component new forces will soon emerge. For example the standard model sector responsible for breaking the electroweak symmetry spontaneously could be replaced by a new and more natural strongly coupled sector. Another interesting avenue is the possibility that dark matter itself, as the ordinary proton and neutron, is not an elementary particle but rather a state composed by new strongly coupled matter. I will argue that models of composite dynamics are indeed primary candidates for a better understanding of the origin of bright and dark matter in Nature. I will also show by reviewing recent results that Lattice field methods can provide clear and falsifiable predictions for the Large Hadron Collider experiment as well as dark matter direct experiments searches. I will then switch gear and exhibit the first proof of existence of nonsupersymmetric and non-asymptotically free 4D Gauge-Yukawa theories (structurally similar to the standard model) that are UV finite thanks to the existence of an exact interacting quantum UV fixed point in the gauge, Yukawa and scalar self-couplings. These theories are therefore asymptotically safe. Theories with this behaviour have been searched for on the lattice for the past several decades. The results show the critical ingredients needed to construct these new classes of theories and offer a strategy for future lattice studies. I will then comment on the wide phenomenological impact of this discovery. Understanding strong dynamics is therefore crucial to construct phenomenologically relevant extensions of the standard model of particle interactions, and ultimately being able to describe our Universe.
        Speaker: Prof. Francesco Sannino (CP3-Origins & DIAS, University of Southern Denmark)
        Slides
    • 10:30
      Coffee Break
    • Plenary Session BF & 1F Main Hall

      BF & 1F Main Hall

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Rainer Sommer
      • 242
        Physics at the Belle II experiment and Lattice QCD
        In this talk, we discuss the status and future prospect in B physics, focusing on the SuperKEKB/Belle II experiment at KEK, with the target peak luminosity of 8 \times 10^35 cm^-2s^-1. This enables us to study decays of heavy flavor particles, B and D mesons as well as \tau leptons, at the order of O(10^10) per year, and to search for New Physics through processes sensitive to presence of virtual heavy particles. We present physics reach, status and plan of the experiment. We also discuss the inputs necessary from lattice QCD calculations, in order to identify effects of New Physics.
        Speaker: Prof. Toru Iijima (KMI, Nagoya University)
        Slides
      • 243
        Progress and prospects for heavy flavour physics on the lattice
        I will review recent progress in lattice computations relevant for B- and charm physics. Emphasis will be put on the interplay with the upcoming new generation of experimental results.
        Speaker: Prof. Carlos Pena (Universidad Autónoma de Madrid and IFT-UAM/CSIC)
        Slides
    • 12:30
      Lunch Break
    • Applications Beyond QCD 405

      405

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Eigo Shintani
      • 244
        Complex Langevin simulation in condensed matter physics
        The complex Langevin method is one hopeful candidate to tackle the sign problem. This method is applicable not only to QCD but also to nonrelativistic field theory. We present the simulation results of nonrelativistic condensed matter systems.
        Speaker: Arata Yamamoto (The University of Tokyo)
        Slides
      • 245
        Domain-wall/overlap fermion and topological insulators
        Topological insulators are a new class of materials which have gapped spectra in the bulk, but are accompanied by topologically protected gapless excitations at the surface (edge) of the system. These phenomena have a close relationship with symmetry and dimensionality of the system through quantum anomalies. We point out that such a surface state is a physical realization of the domain-wall/overlap fermion. From this point of view, we discuss its implications for experiments of topological insulators. We also discuss an unconventional overlap fermion, which is suggested by the “periodic table" of topological insulators and superconductors.
        Speaker: Dr Taro Kimura (RIKEN)
        Slides
      • 246
        Chiral Magnetic Conductivity in an interacting lattice model of a parity-breaking Weyl semimetal
        With the recent experimental realization of Dirac and Weyl semimetals it has become possible to study the Chiral Magnetic Effect (CME) in clean table-top experiments. We present results of a mean-field study of the static Chiral Magnetic Conductivity (sCMC) in a simple lattice model of a parity-breaking Weyl semimetal. Our model is given by the lattice Wilson-Dirac Hamiltonian with on-site repulsive interaction and a constant chiral chemical potential term. The value of the sCMC for free Dirac fermions is well known and possible corrections due to fermion-interactions are often neglected. We study the model in the parameter space of bare Wilson-Dirac mass, inter-fermion interaction strength and bare chiral chemical potential and find that on-site repulsive interactions affect the sCMC almost exclusively through the enhancement of the renormalized chiral chemical potential. Non-trivial corrections to the sCMC due to inter-fermion interactions seem to be irrelevant in practice, since they become important only in a phase where the CME response is strongly suppressed by a large gap in the energy spectrum.
        Speaker: Mr Matthias Puhr (Regensburg University)
        Slides
      • 247
        Hybrid Monte Carlo simulations of graphene in external magnetic field
        Recent experimental results [Nature Physics, 8 (2012) 550 and Nature 505 (2014) 528] indicate that graphene turns into insulator in sufficiently strong magnetic field. However, the exact nature of this state is still elusive and there are large discrepancies between theoretical predictions and experimental results. To resolve this discrepancies extensive simulations of graphene in external magnetic field were performed using Hybrid Monte Carlo algorithm. Insulating state was observed in agreement with experiment. Mass gap and various order parameters were measured.
        Speaker: Dr Maksim Ulybyshev (Regensburg University)
        Slides
      • 248
        Numerical simulations of graphene conductivity with realistic inter-electron potential
        There was a long time disagreement in literature between theoretical calculations and experimental data about conductor-insulator phase transition in graphene. This contradiction was resolved in the papers [Phys. Rev. Lett. 111, 56801 (2013)] and [Phys. Rev. B 89, 195429 (2014)] where chiral condensate was studied taking into account screening of Coulomb potential by $\sigma$-electrons of carbon atoms. Unfortunately there was no explicit evidence that conductivity drops down simultaneously with appearance of chiral condensate. Therefore it is required to perform direct calculations of graphene conductivity in tight-binding model. We performed Monte-Carlo simulations of graphene tight-biding model with realistic electron-electron interaction and calculated the dependence of conductivity on dielectric permittivity of substrate. According to our calculations suspended graphene has conductivity $\sigma = 0.325\pm0.003 e2/\hbar$ which is very close to experimental value $\sigma = (1.01 \pm 0.04)e^2/4\hbar$. The phase transition in conductivity is also shifted to unphysical region epsilon<1 due to influence of $\sigma$-electrons.
        Speaker: Mr Denis Boyda (Far Eastern Federal University)
        Slides
      • 249
        Lattice gauge theory treatment of strongly correlated Dirac semimetals
        We study the strong correlation effect from Coulomb interaction in two- and three-dimensional Dirac semimetals. Effective field theory for the interacting Dirac semimetals can be constructed in terms of quantum electrodynamics (QED). We propose a U(1) lattice gauge theory formulation with the internal degrees of freedom, such as spin and orbitals (pseudospin), explicitly. In the strong coupling limit, the coupling to the U(1) gauge field leads to local charge neutrality, which can be expressed as a sufficiently strong Hubbard (on-site) repulsion as well, rendering the system into a Mott insulator. We discuss the evolution of charge and spin orders by strong coupling expansion and mean-field analysis, and examine the effect of dimensionality and number of flavors.
        Speaker: Dr Yasufumi Araki (Institute for Materials Research, Tohoku University)
        Slides
    • Nonzero Temperature and Density 403

      403

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Tilo Wettig
      • 250
        Two-Color Lattice QCD with Non-zero Chiral Chemical Potential
        The phase diagram of two-color QCD with non-zero chiral chemical potential is studied by means of lattice simulation. We focus on the influence of a chiral chemical potential on the confinement/deconfinement phase transition and the breaking/restoration of chiral symmetry. The simulation is carried out with dynamical staggered fermions without rooting. The dependences of the Polyakov loop, the chiral condensate and the corresponding susceptibilities on the chiral chemical potential and the temperature are presented. The critical temperature is observed to increase with increasing chiral chemical potential.
        Speaker: Mr Andrey Kotov (ITEP)
        Slides
      • 251
        Lattice simulation of $QC_2D$ with $N_f=2$ at non-zero baryon density
        We present the results of lattice simulation of $QC_2D$ with two flavors of staggered fermions and non-zero chemical potential (\mu_q). Dependencies of the Polyakov loop, chiral condensate and baryon number density on \mu_q were studied. We found, that raising of the baryon chemical potential leads to the chiral symmetry restoration. At small \mu_q our results for the baryon density agree with ChPT predictions.
        Speaker: Mr Alexander Nikolaev (Far Eastern Federal University)
        Slides
      • 252
        Cluster expansions and chiral symmetry at large density in 2-color QCD
        We study lattice $SU(N_c)$ gauge theories with massless staggered fermions at finite quark chemical potential $\mu$. At strong coupling such theories have been studied with a variety of techniques such as mean field, monomer-dimer representations and MC simulations. Here we employ a new cluster expansion which has recently been shown to converge for large $\mu$ and small $\beta$. Extension to weaker coupling, however, does not appear feasible in the presence of complex fermion determinant. For theories with real determinant, however, such as 2-color QCD, we show how such large $\mu$ cluster expansions can be used to obtain information on the behavior of lattice observables in the weak coupling regime. We discuss the intertwined questions of presence of chiral symmetry and quarkyonic vs. superfluidity/superconductivity phases at large $\mu$, and the problem of saturation, a lattice artifact, and means to circumvent it.
        Speaker: Terry E. Tomboulis (UCLA)
        Slides
      • 253
        Thimble regularization at work besides toy models: from Random Matrix Theory to Gauge Theories.
        Thimble regularization as a solution of the sign problem has been successfully put at work for a few toy models. Given the non trivial nature of the method (also from the algorithmic point of view) it is compelling to provide evidences that it works for realistic models. A chiral random matrix theory has been studied in detail. The known analytical solution shows that the model is non-trivial as for the sign problem (in particular, phase quenched results can be very far away from the exact solution). This study gave us the chance to address a couple of key issues: how many thimbles contribute to the solution of a realistic problem? can one devise algorithms which are robust as for staying on the correct manifold? The obvious step forward we are interested in are applications to gauge theories.
        Speaker: Mr Giovanni Eruzzi (University of Parma and INFN)
        Slides
      • 254
        Thimble regularization at work for Gauge Theories: from toy models onwards.
        A final goal for thimble regularization of lattice field theories is the application to lattice QCD and the study of the phase diagram. Gauge theories pose a number of conceptual and algorithmic problems, some of which can be addressed even in the framework of toy models. We report on our progresses in this field, starting in particular from first successes in the study of one link models.
        Speaker: Dr Francesco Di Renzo (University of Parma and INFN)
        Slides
      • 255
        Application of the Lefschetz thimble formulation to the (0+1) dimensional Thirring model at finite density
        We investigate the (0+1) dimensional Thirring model at finite chemical potential by HMC simulations based on the Lefschetz thimble formulation of path-integration. We adopt the lattice model defined with staggered fermion and a compact link field for the auxiliary vector field. We firstly find the critical points (saddle points) under the assumption that the complexified link field is space-time independent, and investigate gradient flows and the Stokes phenomenon to identify the thimbles which contribute to the path-integral. Then, we perform lattice simulations using HMC algorithm with taking one particular thimble and the results are compared to the exact solution. The numerical results are in agreement with the exact ones in small and large chemical potential regions, while we observe a discrepancy between the numerical and exact results around the intermediate transition region in the chemical potential. We also discuss the result in relation to the contributions of the other thimbles.
        Speaker: Dr Syo Kamata (Rikkyo Univ.)
        Slides
    • Nonzero Temperature and Density 401

      401

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Takashi Umeda
      • 256
        Lattice NRQCD study of quarkonium at non-zero temperature
        To study the in-medium modification of quarkonium properties, charmonium correlators at $140.4 (\beta =6.664) \le T \le 221 (\beta = 7.280)$ (MeV) are calculated using the NRQCD formalism on $48^3 \times 12$ gauge configurations with dynamical $N_f = 2 + 1$ flavors of Highly Improved Staggered Quarks (HISQ). To determine the "zero energy shift" for these lattices, we perform a fine zero temperature scan ($\beta =6.664, 6.740, 6.800, 6.880, 6.950, 7.030, 7.150$ and $7.280$). In addition, we investigate the influence of statistical errors on the reconstruction of spectral functions by varying the number of correlators in the study. We find that the temperature dependence of charmonium correlators is stronger than the temperature dependence of bottomonium correlators in a given channel. This fits into the expected pattern of sequential quarkonium melting.
        Speaker: Prof. Seyong Kim (Sejong University)
        Slides
      • 257
        Charmonia and bottomonia at finite temperature on large quenched lattice
        We present our updated study on charmonia and bottomonia at finite temperature with quenched ensembles on large and fine isotropic lattices. Simulations have been performed by using the standard plaquette gauge and the $O(a)$-improved Wilson fermion actions in a region of quark mass for charmonia to bottomonia in order to investigate difference of in-medium behavior between them. We show spectral functions of quarkonia at temperatures in a range between about 0.73$T_c$ and 2.2$T_c$ at both vanishing and finite momenta and discuss about their temperature, quark mass and momentum dependence.
        Speaker: Dr Hiroshi Ohno (Center for Computational Sciences, University of Tsukuba)
        Slides
      • 258
        Transverse and longitudinal spectral functions of charmonia at finite temperature with maximum entropy method
        We study the charmonium spectral functions with finite momentum from lattice Euclidean correlators using the maximum entropy method. In medium, the spectral function of vector channel with finite momentum is decomposed into transverse and longitudinal components because of the lack of Lorentz invariance. We investigate these spectral functions, their residues and the dispersion relations on the quenched lattice below and above the critical temperature.
        Speaker: Mr Atsuro Ikeda (Osaka Universitiy)
        Slides
      • 259
        Polyakov loop correlators and cyclic Wilson loop from lattice QCD
        We calculate the correlation function of Polyakov loops and the cyclic Wilson loops in 2+1 flavor QCD at non-zero temperature. We also study the correlation function of Wilson lines in Coulomb gauge. In our investigations we use the highly improved staggered quark (HISQ) action and lattices with temporal extent Nt=4,6,8,10 and 12. At high temperatures we compare our numerical results with perturbation theory.
        Speaker: Mr Johannes Weber (Technical University Munich)
        Slides
      • 260
        Static quark-antiquark pair free energy and screening masses: continuum results at the QCD physical point
        We study the correlators of Polyakov loops, and the corresponding gauge invariant free energy of a static quark-antiquark pair in 2+1 flavor QCD at finite temperature. Our simulations were carried out on Nt = 6, 8, 10, 12, 16 lattices using a Symanzik improved gauge action and a stout improved staggered action with physical quark masses. The free energies calculated from the Polyakov loop correlators are extrapolated to the continuum limit. For the free energies we use a two step renormalization procedure that only uses data at finite temperature. We also measure correlators with definite Euclidean time reversal and charge conjugation symmetry to extract two different screening masses, one in the magnetic, and one in the electric sector, to distinguish two different correlation lengths in the full Polyakov loop correlator. This conference contribution is based on the paper: JHEP 1504 (2015) 138
        Speaker: Attila Pasztor (University of Wuppertal)
        Slides
      • 261
        Confinement/deconfinement transition temperature from the Polyakov loop potential and gauge-invariant gluon mass
        In this talk we give an analytical derivation for the existence of the confinement/deconfinement phase transition at finite temperature in Yang-Mills theory. Moreover, we give a numerical estimate on the transition temperature $T_d$ in the form of the ratio to a gauge-invariant gluon mass $M(T)$ measured on the lattice. For this purpose, we use the functional renormalization-group equation of the Wetterich type to obtain the effective potential of the Polyakov loop average in the Yang-Mills theory. This result enables us to understand the reason why the phase transition from deconfinement to confinement occurs at a certain temperature and what is the mechanism for confinement at finite temperature.
        Speaker: Prof. Kei-Ichi Kondo (Chiba University)
        Slides
    • Theoretical Developments 402

      402

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Richard Brower
      • 262
        (Dimensional) twisted reduction in large N gauge theories
        We will use the idea of twisted reduction and volume independence to numerically investigate the large N limit of Yang-Mills theories by reducing two of the directions to a single point. We will show that this framework might be convenient to analyze some properties of these theories.
        Speaker: Dr Alberto Ramos (CERN)
        Slides
      • 263
        Large N meson propagators from twisted space-time reduced model
        We propose a new method to calculate meson propagators in the large N limit from twisted space-time reduced model. Preliminary results of meson spectra and decay constants are presented.
        Speaker: Prof. Masanori Okawa (Hiroshima University)
        Slides
      • 264
        Topology and glueballs in $SU(7)$ Yang-Mills with open boundary conditions
        It is well known that the topology of gauge configurations generated in a Markov Monte-Carlo chain freezes as the continuum limit is approached. The corresponding autocorrelation time increases exponentially with the inverse lattice spacing, affecting the ergodicity of the simulation. In SU(N) gauge theories for large N this problem sets in at much coarser lattice spacings than for N=3. This means that its systematics can be studied on lattices that are smaller in terms of the number of lattice sites. It has been shown that using open boundary conditions in time allows instantons to be created and destroyed, restoring topological mobility and ergodicity. However, with open boundary conditions translational invariance is lost and the influence of spurious states propagating from the boundary into the bulk on physical correlators needs to be carefully evaluated. Moreover, while the total topological charge can be changed, the mobility of instantons across the lattice is still reduced. We consider SU(7) Yang-Mills theory and analyse its topological content in the periodic and open boundary condition cases. We also investigate scalar and pseudoscalar glueball correlation functions.
        Speaker: Dr Alessandro Amato (University of Helsinki)
        Slides
      • 265
        Diagrammatic Monte-Carlo algorithms for large-N quantum field theories from Schwinger-Dyson equations
        I present a general framework for constructing Diagrammatic Monte-Carlo algorithms for large-N quantum field theories, which is based on the stochastic solution of the full untruncated hierarchy of Schwinger-Dyson equations and is an extension of the approach proposed in arXiv:1104.3459, arXiv:1009.4033. The algorithms are capable of constructing both weak- or strong-coupling expansions which are convergent at large N and asymptotic 1/N expansions. We illustrate the application of the algorithm on the examples of the planar phi^4 theory, U(N) sigma model and strong-coupling QCD in the Veneziano limit at finite density. In the latter cases of lattice field theories with U(N) degrees of freedom, simulations in the weak-coupling regime are hindered by a sign problem. We discuss possible ways to overcome these problems, e.g. using the "bold" diagrammatic Monte-Carlo.
        Speaker: Dr Pavel Buividovich (Regensburg University)
        Slides
      • 266
        Numerical study of complex instantons in the Gross-Witten U(N) matrix model
        It is known that the 1/N expansion in the Gross-Witten U(N) matrix model is divergent and leading divergence in the weak coupling regime is governed by instanton which is associated with eigenvalue tunneling. However, it was not known what triggers divergences in the strong coupling phase where instantons disappear. To address this problem, we numerically study structure of saddle points in the Gross-Witten model with complexified degrees of freedom. These complex saddle points can be considered as complex instantons. In the strong coupling phase we identify complex instanton which governs leading divergence in the 1/N expansion. Also we find, that in the limit of large N in the strong coupling regime action of the perturbative vacuum has one flat direction. We present an evidence that in the vicinity of the transition point all complex instantons might give a significant contribution to the partition function, which might be important observation for Hybrid Monte-Carlo simulations based on integration on a Lefschetz thimbles.
        Speaker: Mr Semen Valgushev (University of Regensburg)
        Slides
      • 267
        Matrix Geometry and Coherent States
        Some matrix models have been proposed as a non-perturbaive formulation of string theories. In the matrix model formulation of string theories, configurations of strings or D-branes are described by a set of Hermitian matrices. Here, we propose a new class of observables in matrix models, which are made of the Hermitian matrices and encode geometric information of the strings or D-branes. By performing a Monte Carlo simulation and computing those observables for a simple toy model of a bosonic matrix model, we demonstrate how we can see the geometric properties of the strings from the matrix configurations.
        Speaker: Prof. Goro Ishiki (University of Tsukuba)
        Slides
    • Vacuum Structure and Confinement 404

      404

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Akihiro Shibata
      • 268
        Chiral symmetry breaking, instantons, and monopoles
        The purpose of this study is to show that monopoles are related to chiral symmetry breaking and to instantons, using the Overlap fermions as an analytical tool. In order to study these relations, we add one pair of monopoles with different magnetic charges to the quenched SU(3) configurations by a monopole creation operator. We then count the number of zero modes, and compute the average square of the topological charges as the number of instantons. We have shown that one pair of monopoles makes one instanton, and that the chiral condensate which is computed from the eigenvalues and eigenvectors of the Overlap Dirac operator decreases by increasing the charges of the monopoles. We have also found that the fermion spectral density of the eigenvalues increases by increasing the charges of the monopoles. Moreover, we compare the low-lying eigenvalues with the prediction of the random matrix theory. By now these studies have been done using one small lattice ($V = 14^{4}, \beta = 6.00$). Our program is to investigate this subject more precisely using larger lattices and different lattice spacings. In this study, we present results on the chiral condensate, the pseudo scalar mass, and pion decay constants on larger lattice.
        Speaker: Dr Masayasu Hasegawa (Joint Institute for Nuclear Research)
        Slides
      • 269
        Classifying the phases of gauge theories by spectral density of probing chiral quarks
        We describe our recent proposal that distinct phases of vector-like gauge theories with fundamental quarks translate into specific types of low-energy behavior in Dirac spectral density. The resulting scenario is built around new evidence substantiating the existence of a phase characterized by bimodal (anomalous) density, and corresponding to deconfined dynamics with broken valence chiral symmetry. We argue that such anomalous phase occurs quite generically in these theories, including in "real world" QCD above the crossover temperature, and in zero-temperature systems with many light flavors.
        Speaker: Ivan Horvath (University of Kentucky)
        Slides
      • 270
        Dirac spectrum representation of Polyakov loop fluctuations in lattice QCD
        We investigate contribution to the Polyakov loop fluctuations from each Dirac mode. The Polyakov loop fluctuations are sensitive probes for the quark deconfinement even if considering dynamical quarks. We derive analytical relations between the Polyakov loop fluctuations and Dirac modes on the temporally odd-number lattice, where the temporal lattice size is odd, with the normal non-twisted periodic boundary condition for link-variables. These relations indicate that low-lying Dirac modes have little contribution to the Polyakov loop fluctuations. We numerically confirmed that at the quenched level in both confinement and deconfinement phases. Our results suggest that there is no direct one-to-one correspondence between confinement and chiral symmetry breaking in QCD.
        Speaker: Mr Takahiro Doi (Kyoto University)
        Slides
      • 271
        SU(4) symmetry of hadrons upon quasi-zero Dirac mode elimination.
        A large degeneracy of mesons of given spin has been discovered upon reduction of the quasi-zero modes of the Dirac operator in a dynamical simulation. One expects a priori, that upon elimination of the quasi-zero modes of the Dirac operator, the chiral symmetry should be restored, since the quark condensate of the vacuum is connected with the density of the quasi-zero modes. However, it has turned out that not only degeneracy patterns from the SU(2)_R * SU(2)_L and U(1)_A symmetries are observed, but a larger degeneracy that includes all possible chiral multiplets for given spin. This symmetry has been established to be SU(4) that includes both the isospin rotations of quarks of given chirality as well as the rotations chirality itself. This symmetry is higher than the symmetry of the QCD Lagrangian and should be consequently considered as an emergent symmetry that reflects the QCD dynamics without the quasi-zero modes of the Dirac operator. This symmetry implies the absence of the color-magnetic field in the system and might be interpreted as a manifestation of the dynamical QCD string.
        Speaker: Prof. Leonid Glozman (University of Graz)
        Slides
      • 272
        Gauge fixing and the gluon propagator in renormalizable xi gauges
        Covariant $R_\xi$ gauge fixing is notoriously difficult for large lattice volumes, large $\xi$ and small $N_c$. We thoroughly test different convergence techniques, which allows the gauge fixing of lattice configurations with a total volume of (3.25 fm)$^4$, up to $\xi=0.5$. We are able to study the gluon propagator in the infrared region and its dependence on the gauge fixing parameter $\xi$. As expected, the longitudinal gluon dressing functions stay constant at their tree-level value $\xi$. Similar to the Landau gauge, the transverse $R_\xi$ gauge gluon propagators saturate at a non-vanishing value in the deep infrared for all values of $\xi$ studied. We compare with very recent continuum studies and perform a simple analysis of the found saturation with a dynamically generated effective gluon mass.
        Speaker: Pedro Bicudo (IST, Univ. Lisboa)
        Slides
      • 273
        Non-Gaussianity of the topological charge distribution in SU(3) Yang–Mills theory
        We report the results of a study of the distribution of the topological charge of the SU(3) Yang–Mills theory: we estimate the second and fourth cumulant with high precision in order to measure the deviation from the normal distribution. The computation is done on the lattice implementing a naïve discretization of the topological charge evolved with the Yang–Mills gradient flow. A range of high statistics Monte Carlo simulations with different lattice volumes and spacings is used to extrapolate the results to the continuum limit with confidence by keeping finite volume effects negligible with respect to the statistical errors.
        Speaker: Marco Cè (Scuola Normale Superiore, Pisa & INFN, Sezione di Pisa)
        Slides
    • Weak Decays and Matrix Elements 406

      406

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Takashi Kaneko
      • 274
        Kaon semileptonic form factors as functions of the momentum transfer with Nf=2+1+1 Twisted Mass fermions
        We present lattice results for the form factors relevant for K -> pi l nu decays, obtained from simulations performed by the European Twisted Mass Collaboration with Nf=2+1+1 flavors of dynamical quarks, at three values of the lattice spacing and pion masses as low as 250 MeV. Our determination of f_+(0), which in combination with the experimental result for (f_+(0)*V_us) can be used to determine the CKM matrix element, will be presented along with a study of both the vector and the scalar form factors as functions of the momentum transfer comparing our results with the experimental ones.
        Speaker: Dr Lorenzo Riggio (INFN Roma Tre)
        Slides
      • 275
        Long-distance contributions to the rare kaon decay K+ -> pi+ l+ l-
        The rare decays of a kaon into a pion and a charged lepton/antilepton pair proceed via a flavour changing neutral current and therefore may only be induced beyond tree level in the Standard Model. This natural suppression makes these decays sensitive to the effects of potential New Physics. To discern such New Physics one must be able to control the errors on the Standard Model prediction of the decay amplitude. These particular decay channels however are dominated by a single photon exchange; this involves a sizeable long-distance hadronic contribution which represents the current major source of theoretical uncertainty. In this talk I will outline our methodology for the computation of the long distance contributions to these rare decay amplitudes using lattice QCD, and present the numerical results of some exploratory studies using the Domain Wall Fermion ensembles of the RBC and UKQCD collaborations.
        Speaker: Mr Andrew Lawson (University of Southampton)
        Slides
      • 276
        Long-distance contributions to the rare kaon decay K+ -> pi+ nu nu-bar
        As flavor changing neutral current processes, rare kaon decays K -> pi nu nu-bar are highly suppressed in the standard model and thus provide ideal probes for the observation of new physics effects. Given the importance of rare kaon decays, the CERN NA62 experiment aims at an observation of O(100) events of K+ -> pi+ nu nu-bar decay and a 10%-precision measurement of the branching ratio in two years. Another experiment, J-PARC KOTO, is dedicated to search for the CP-violating KL -> pi0 nu nu-bar decay. Recognizing that the standard model predictions will soon be confronted with new experimental results, it is important to determine the long-distance contributions to K+ -> pi+ nu nu-bar with a controlled error. In this talk I will present an exploratory calculation of the long-distance contributions to the K+ -> pi+ nu nu-bar decay amplitude from first principles using lattice QCD.
        Speaker: Dr Xu Feng (Columbia University)
        Slides
      • 277
        Computing the long-distance contributions to $\epsilon_K$
        We will discuss the theoretical framework and initial results for the lattice QCD calculation of the long-distance contributions to $\epsilon_K$, the measure of indirect CP violation in $K_L$ decay. A proof-of-principle calculation has been carried out on a $24^3\times 64$ lattice volume with an inverse lattice spacing of 1.73 GeV and pion mass of 329 MeV. A complete calculation will be described which shows how the logarithmically divergent lattice contribution can be combined with continuum perturbation theory to obtain a physical result.
        Speaker: Prof. Norman Christ (Columbia University)
        Slides
      • 278
        Emerging lattice approach to K-Unitarity Triangle
        It has been been clear for past many years that in low energy observables new physics can only appear as a perturbation. Therefore precise theoretical predictions and precise experimental measurements have become rather mandatory. Here we draw attention to the significant advances that have been made on the lattice in recent years in K=> pi pi, Delta_M_K, long- distance part of epsilon_K and rare K-decays. Thus, in conjunction with experiments, the construction of a unitarity triangle purely from Kaon physics should soon become feasible. In comparison with the B-unitarity triangle, this should allow for more stringent tests of the Standard Model and tighter constraints on new physics.
        Speaker: AMARJIT SONI (Brookhaven National Lab)
        Slides
      • 279
        Lattice QCD calculations of nucleon transverse momentum-dependent parton distributions (TMDs) at 170 MeV pion mass
        An exploration of nucleon TMD observables at a substantially lower pion mass, 170 MeV, than used in previous lattice TMD calculations is presented. On a corresponding RBC/UKQCD DWF ensemble, TMDs are extracted from nucleon matrix elements of a bilocal quark operator containing a staple-shaped gauge link. Appropriate TMD ratios are constructed to cancel divergences associated with the gauge link. In particular, time-reversal odd observables associated with the Sivers and Boer-Mulders effects are considered. The results are compared with previous DWF calculations at 297 MeV pion mass with a view to exploring whether these observables vary strongly as a function of pion mass in the chiral regime.
        Speaker: Michael Engelhardt (NMSU)
        Slides
    • 16:00
      Coffee Break
    • Algorithms and Machines 405

      405

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Hideo Matsufuru
      • 280
        QUDA features, scaling and solvers
        We describe recent developments to the QUDA software framework, a library aimed at deploying lattice QCD computations on GPUs. The library has ever broadening support for various LQCD actions and algorithms, as well as being integrated into many LQCD applications. Recent focus has been on improving strong scaling for multi-GPU calculations and developing an adaptive multigrid solver. We give updates on both of these efforts, and compare performance against other platforms. Lastly, we look to the future and discuss how upcoming technologies such as stacked memory and nvlink, a fast GPU interconnect, will bring disruptive changes to lattice QCD calculations.
        Speaker: Dr M Clark (NVIDIA)
        Slides
      • 281
        Optimization of Lattice QCD with CG and multi-shift CG on Intel Xeon Phi Coprocessor
        We implemented lattice QCD on Xeon Phi coprocessor using intrinsics as vectorization method, and OpenMP and MPI as parallelization method. Our implementation uses double precision conjugate gradient (CG) solver which also supports multi-shift CG. We present our optimization methodology and performance for key steps in CG algorithms.
        Speaker: Dr hirokazu Kobayashi (Intel K. K.)
        Slides
      • 282
        Accelerating twisted mass LQCD with QPhiX
        We present the implementation of twisted mass fermions within the QPhiX code package along with an interface to tmLQCD. We analyze the performance on Intel Xeon Phis and Intel Xeon Haswell CPUs. On the latter we find a speedup of larger than 4x for the mixed precision conjugate gradient inverter.
        Speaker: Dr Mario Schroeck (INFN - Roma Tre)
        Slides
      • 283
        An implementation of hybrid parallel CUDA code for the hyperonic nuclear forces
        We present our recent effort to develop the CUDA code to measure the Nambu-Bethe-Salpeter (NBS) wave function of various baryon-baryon systems on the lattice, which is a primary quantity to make a first-principle calculation of nuclear force and hyperonic-nuclear forces from lattice QCD. Not only the nucleon-nucleon (NN) but also the hyperon-nucleon (YN) and hyperon-hyperon (YY) interactions are fundamental inputs to study the properties of (hyper-)nuclei and the (hyperonic-)matter inside the neutron stars while the present phenomenological YN and YY potentials are not well constrained from experimental data. In the HAL QCD method, which is recently developed to study various interhadron interactions from lattice QCD, the NN, YN and YY potentials can be obtained at the same time by defining the set of interpolating field of baryons. The purpose of this contribution is to present a fast algorithm to calculate a large set of NBS wave functions and its implementation to the heterogeneous system involving GPUs. This is a branched work from the hybrid parallel programming reported in LATTICE 2013 [PoS LATTICE2013 (2014) 426]. A hybrid parallel CUDA code is implemented with utilizing the MPI and OpenMP. The present code works on HA-PACS at University of Tsukuba which comprises 16 CPUs (Intel E5-2670) and 4 GPUs (NVIDIA M2090) per node. We would also like to present preliminary results of YN potentials which is obtained from the 2+1 flavor lattice QCD at almost physical point.
        Speaker: Dr Hidekatsu Nemura (Center for Computational Sciences, University of Tsukuba)
        Slides
      • 284
        Early Performance Evaluation of Lattice QCD on POWER+GPU Cluster
        As supercomputers are shifting from peta-scale to exa-scale, computers with accelerators such as GPUs, MICs and FPGAs have become one of the big trends of supercomputer because of their low energy consumption and high density. Now IBM’s POWER processor has quite new power, Nvidia’s Tesla GPU brings huge computational capability. It is important for us to understand how this new POWER+GPU environment brings power to the actual applications in the early stage. We implemented Wilson-Dirac kernel and BiCGStab solver using CUDA7.0 on the POWER+GPU cluster and evaluated the performance. We will show the performance results of serial and parallel executions.
        Speaker: Mr Jun Doi (IBM Research - Tokyo)
        Slides
    • Nonzero Temperature and Density 403

      403

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Keitaro Nagata
      • 285
        Preweighting method in Monte-Carlo sampling with complex action
        Understanding the QCD phase diagram from first principles is one of the ultimate goals in nuclear and particle physics. The largest obstacle towards this goal is the sign problem in lattice QCD at finite chemical potential. The statistical weight is reduced by the average phase factor and it becomes difficult to obtain precise results in the phase reweighting method. If we know the amount of cancellation beforehand, it would be possible to take account of the reduction during the sampling of configurations. In the case where the phase distribution is Gaussian, the average phase factor is exp(-D^2/2) [1], where D^2 is the variance of the phase. We can approximately include this reduction in the Monte-Carlo sampling by adding "penalty" terms to the action. We examine the usefulness of this "preweighting" method in the auxiliary field Monte-Carlo configurations in the strong coupling limit of lattice QCD [2]. We will also apply the preweighting to the strong coupling lattice QCD with 1/g^2 corrections [3], where the sign problem is more serious. [1] S. Ejiri, Eur. Phys. J. A 49 (2013) 86. [2] T. Ichihara, A. Ohnishi and T. Z. Nakano, Prog. Theor. Exp. Phys. 2014 (2014), 123D02. [3] T. Ichihara, A. Ohnishi, PoS LATTICE2014 (2014), 188.
        Speaker: Prof. Akira Ohnishi (Yukawa Institute for Theoretical Physics, Kyoto University)
        Slides
      • 286
        The density of states approach at finite chemical potential: a numerical study of the Bose gas.
        Recently, a novel algorithm for computing the density of states in statistical systems and quantum field theories has been proposed. The same method can be applied to theories at finite density affected by the notorious sign problem, reducing a high-dimensional oscillating integral to a more tractable one-dimensional one. As an example we applied the method to the relativistic Bose gas.
        Speaker: Dr Roberto Pellegrini (The University of Edinburgh)
        Slides
      • 287
        Solving the complex action problem of the finite density Z3 spin model with the density of states approach using FFA
        In this contribution we apply a variant of the density of states method to the Z3 spin model with a chemical potential. We use a restricted Monte Carlo strategy to compute restricted expectation values and study their dependence on a free parameter \lambda which enters in the Boltzmann weight. When expressed in terms of the density, the expectation values are well known functions of \lambda which we fit to the Monte Carlo data and in this way determine the density of states (functional fit approach (FFA)). We calculate observables related to the particle number and the particle number susceptibility and compare the results to a reference simulation in the dual formulation of the Z3 spin model. We find a good agreement of the methods for a wide range of parameters.
        Speaker: Mr Pascal Toerek (University of Graz)
        Slides
      • 288
        Density of states approach with FFA for an effective Polyakov loop model at finite density
        We study an effective theory for Polyakov loops at finite density (SU(3) spin model) using the density of states (DoS) method. We generalize a recently developed variant of DoS, the so-called functional fit approach (FFA) for systems with continuous degrees of freedom (compare also the related talk by P. Törek). We show that the density of states can be computed with sufficiently high accuracy and we compare our physical observables to the results of a reference simulation of the model in a dual simulation.
        Speaker: Mr Mario Giuliani (University of Graz)
        Slides
      • 289
        Dual representation for massless fermions with chemical potential and U(1) gauge fields
        We show that the massless Schwinger model with staggered fermions can be represented with dual variables which are loops for the fermions and surfaces for the gauge fields. In the dual form the complex action problem of the conventional representation at finite chemical potential or non-vanishing theta-angle is absent, and dual Monte Carlo simulations are possible at arbitrary values of these parameters. We show that the dualization can be generalized to 1+1 dimensional fermions interacting with 4 dimensional U(1) gauge fields, i.e., to a system of wires with relativistic charged fermions.
        Speaker: Prof. Christof Gattringer (University of Graz)
        Slides
    • Nonzero Temperature and Density 401

      401

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Masakiyo Kitazawa
      • 290
        A stochastic approach to the reconstruction of spectral functions in lattice QCD
        All the information of quantum states are encoded in their spectral functions. The reconstruction of the spectral function from the corresponding Euclidean correlation function is, however, a typical ill-posed problem as an analytic continuation from imaginary time to real time is required. The commonly used Maximum Entropy Method (MEM) is based on the Bayesian interference and its output thus may have bias on the prior information provided in the default models. We present results using a stochastic approach for the reconstruction of the spectral functions from the Euclidean correlation functions. In this approach the spectral function is parameterized as a sum of randomly distributed boxes. By varying the width, location and the height of the box stochastically an optimal spectral function can be obtained. Using this approach we can reproduce fairly well the mock spectral functions which contain very sharp resonance peaks, transport peaks and continuum spectra. Application to real lattice QCD charmonium correlator data with N_tau=96, 48, 32 and 24 has also been performed and the same conclusion on the dissociation temperature of ground charmonium states as obtained from MEM is obtained.
        Speaker: Mr Hai-Tao Shu (Key Laboratory of Quark & Lepton Physics (MOE) and Institute of Particle Physics,Central China Normal University)
        Slides
      • 291
        Thermal dilepton rates and electrical conductivity of the QGP
        We present new results on the reconstruction of mesonic continuum spectral functions for three temperatures between 1.1Tc and 1.5Tc in quenched QCD. Making use of non-perturbatively improved clover Wilson valence quarks allows for a clean extrapolation of measured correlator data to the continuum. We found that the vector correlation function is almost temperature independent in the current temperature window. For the case of vanishing momentum the vector spectral function is obtained by fitting the data to a well motivated ansatz. The electrical conductivity of the hot medium, related to the slope of the vector spectral function, is computed from the resulting parameters leading to a temperature independent estimate for all three temperatures. The dilepton rates resulting from the obtained spectral functions are compared to a perturbative expression, obtained by combining the results of different computations in different frequency regimes. Our results qualitatively yield a larger rate in the intermediate frequency regime, and also show no significant temperature dependence. The continuum extrapolated transversally and longitudinally polarized correlation functions at different momenta are compared to the corresponding perturbative expressions.
        Speaker: Mr Florian Meyer (Bielefeld University, Germany)
        Slides
      • 292
        Calculation of free baryon spectral densities at finite temperature
        Following a recent study of nucleon parity doubling at finite temperature from the computation of the two-point nucleon correlators, we study the spectral functions of free baryons at finite temperature. A full analytic review of the spectral densities in the continuum is presented along with a comparison with numerical results. Particular attention will be given to lattice artifacts at higher energies.
        Speaker: Ms Chrisanthi Praki (Swansea University)
        Slides
      • 293
        Nucleons and parity doubling across the deconfinement transition
        The spectrum of nucleons and their parity partners is studied as a function of temperature across the deconfinement transition. We analyse our results using the correlation functions directly, and the Maximum Entropy Method. Both techniques show that there is degeneracy (i.e. parity restoration) in the parity partners' ground states above T_c. This in accordance with the expectation that parity and chiral symmetry are restored in the deconfined phase. We also find that the nucleon ground state is largely independent of the temperature, whereas there are substantial temperature effects in the negative parity (N*) channel, already in the confined phase. All results are obtained using our FASTSUM N_f=2+1 ensembles.
        Speaker: Chris Allton (Swansea University)
        Slides
      • 294
        Study of entropy production in Yang-Mills theory with use of Husimi function
        Understanding the themalization process in a pure quantum system is a challenge in theoretical physics. In this work, we explore possible thermalization mechanism in Yang-Mills(Y-M) theory by using a positive semi-definite quantum distribution function called Husimi function which is given by a coarse graining of Wigner function within the minimal uncertainty. Then entropy is defined in terms of the Husimi function, which is called the Husimi-Wehrl(H-W) entropy. We propose two numerical methods to calculate the H-W entropy. We find that it is feasible to apply the semi-classical approximation with the use of classical Y-M equation. It should be noted that the semi-classical approximation is valid in the systems of physical interest including the early stage of heavy-ion collisions. Using a product ansatz for the Husimi function, which is checked to give only some 10 % over estimate, we succeed in a numerical evaluation of H-W entropy of Y-M theory and show that it surely has a finite value and increases in time.
        Speaker: Hidekazu Tsukiji (YITP, Kyoto univ.)
        Slides
    • Theoretical Developments 402

      402

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Masanori Okawa
      • 295
        Lattice Conformal Field Theory on the Reimann Sphere
        Viable non-perturbative methods for quantum field theories on curved manifold are difficult. By adapting features from both the traditional finite element methods (FEM) and simplicial Regge calculus we are developing a Quantum Finite Element Method(QFEM). To test the QFEM approach, we study the $\lambda \phi^4$ on the simplicial lattice for the Riemann Sphere. To reach the Wilson-Fisher fixed point in the continuum requires modifying the lattice Lagrangian by a counter terms which cancels the ultraviolet distortions of classical FEIM simplicial lattice. In addition Fermions are formulated on the Riemann sphere. Both are compared with the exact solutions to Ising $c = 1/2$ 2D CFT field theory. Future directions and application are entertained.
        Speaker: Prof. Richard Brower (Boston University)
      • 296
        Scalar field theory on a 2-sphere using quantum finite element method (QFEM)
        We demonstrate the quantum finite element method (QFEM) by computing known properties of the 2-D Wilson-Fisher fixed point using scalar field theory on 2-sphere.
        Speaker: Dr George Fleming (Yale University)
        Slides
      • 297
        Spherical Finite Elements for Lattice Radial Quantization
        Recent efforts have been made to study conformal field theories in radial quantization on the lattice. The lattice necessarily breaks spherical symmetry down to icosahedral symmetry, so efforts must be made to restore spherical symmetry — even at the classical level — by clever construction of the action. Linear finite elements have been shown to improve the spectrum of the Laplacian. In this work, we introduce spherical finite elements, which are bounded not by chords but by great circle arcs. We compare the spectra of the linear finite element Laplacian to the spherical finite element Laplacian.
        Speaker: Mr Andrew Gasbarro (Yale University)
        Slides
      • 298
        Encoding field theories into gravities
        We propose a method, which encodes the information of a $d$ dimensional quantum field theory into a $d+1$ dimensional gravity in the $1/N$ expansion. We first construct a $d+1$ dimensional field theory from the $d$ dimensional one via the gradient flow equation, whose flow time $t$ represents the energy scale of the system such that $t¥rightarrow 0$ corresponds to the ultra-violet (UV) while $t¥rightarrow¥infty$ to the infra-red (IR). We then define the induced metric from $d+1$ dimensional field operators. We show that the metric defined in this way becomes classical in the large $N$ limit, in a sense that quantum fluctuations of the metric are suppressed as $1/N$ due to the large $N$ factorization property. As a concrete example, we apply our method to the O(N) non-linear $¥sigma$ model in two dimensions. We calculate the induced metric in three dimensions, which is shown to describe asymptotically De Sitter (dS) or Anti De Sitter (AdS) space in the UV limit. We finally discuss several open issues in future studies.
        Speaker: Prof. Sinya Aoki (Kyoto University)
        Slides
      • 299
        Summary of Super Doubler Approach on Exact Lattice Supersymmetry
        Recently we have proposed a lattice SUSY formulation which we may call super doubler approach, where chiral fermion species doublers are identified as super partners of extended supersymmetry. We claim that the super symmetry is exactly kept on the lattice. However the formulation is nonlocal and breaks lattice translational invariance. We argue that these features cause no fundamental difficulties in the continuum limit. Furthermore a naive formulation breaks associativity of the product of fields and thus cannot be applicable to super Yang-Mills theory. We propose a modified super doubler approach to recover the associativity, which is essentially equivalent to the continuum formulation even for a finite lattice constant. We discuss fundamental meanings of the formulation.
        Speaker: Prof. Noboru KAWAMOTO (Hokkaido University)
    • Vacuum Structure and Confinement 404

      404

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Miho Koma
      • 300
        Study of non-perturbative contributions to surface operator within lattice gauge theory
        The most important probes for the phase states of a four-dimensional gauge field theory are the Wilson and t’Hooft line operators that are defined on one-dimensional curves in the space-time. However, for more detail understanding of four-dimensional gauge field theory dynamics and vacuum topology we need additional probes expressed by operators defined on the subspaces with higher dimensions. Possible candidates are operators that are defined on the two-dimensional surface in the four-dimensional space-time. The corresponding spatial surface operators are sensitive to the magnetic flux through a closed surface, what provides a tool for monopole condensate study. The area and volume dependence of the surface operator can provide a probe for the correlated and uncorrelated monopoles and antimonopoles condensation. In the present work the surface operator dependence on the surface area and volume are studied in the confinement and deconfinement phases within the SU(2) gauge field theory on the lattice. It is shown that spatial and temporal surface operators cannot be defined by the perturbation theory expansion only. The leading non-perturbative correction is obtained. Its dependence on inverse lattice coupling constants and temperature is studied. It is shown that the non-perturbative correction has dimension two and does not depend on temperature.
        Speaker: Prof. Alexander Molochkov (School of Biomedicine, Far Eastern Federal University)
        Slides
      • 301
        Abelian monopole or non-Abelian monopole responsible for quark confinement
        We have pointed out that the SU(3) Yang-Mills theory has a new way of reformulation using new field variables (minimal option), in addition to the conventional option adopted by Cho, Faddeev and Niemi (maximal option). The reformulation enables us to change the original non-Abelian gauge field into the new field variables such that one of them called the restricted field gives the dominant contribution to quark confinement in the gauge-independent way. In the minimal option, especially, the restricted field is non-Abelian U(2) and involves the non-Abelian magnetic monopole. In the preceding lattice conferences, we have accumulated the numerical evidences for the non-Abelian magnetic-monopole dominance in addition to the restricted non-Abelian field dominance for quark confinement supporting the non-Abelian dual superconductivity using the minimal option for the SU(3) Yang-Mills theory. This should be compared with the maximal option which is a gauge invarient version of the Abelian projection in the maximal Abelian gauge: the restricted field is Abelian U(1) x U(1) and involves only the Abelian magnetic monopole, just like the Abelian projection. In this talk, we focus on discriminating between two reformulations, i.e., maximal and minimal options of SU(3) Yang-Mills theory for quark confinement from the viewpoint of dual superconductivity. For this purpose, we measure the distribution of the chromoelectric flux connecting a quark and an antiquark and the induced magnetic-monopole current around the flux tube.
        Speaker: Akihiro Shibata (Computing Research Center, KEK)
      • 302
        Gluon propagators near the phase transition in SU(2) gluodynamics
        We study Landau gauge transverse and longitudinal gluon propagators at non-zero temperature in close vicinity of the phase transition in SU(2) gluodynamics. We compute the screening mass and the electric-magnetic asymmetry of the A**2 condensate. Our goal is to provide results free of finite volume and Gribov copies systematic effects and with small scaling violations.
        Speaker: Dr Vitaly Bornyakov (FEFU, Vladivostok)
        Slides
      • 303
        Anatomy of SU(3) flux tubes at finite temperature
        An attempt to adapt the study of color flux tubes to the case of finite temperature has been made. The field is measured through a connected correlator of Wilson loop and plaquette in the spatial sub-lattice. Still the profile of the flux tube resembles the transverse field distribution around an isolated vortex in an ordinary superconductor. The temperature dependence of all the parameters characterizing the flux tube is investigated.
        Speaker: Francesca Cuteri (University of Calabria)
        Slides
      • 304
        The three-quark potential and perfect Abelian dominance in SU(3) lattice QCD
        We study the static three-quark (3Q) potential with high statistics in SU(3) lattice QCD at the quenched level. For all the distances, the 3Q potential is found to be well described by the Y-Ansatz, i.e., one-gluon-exchange Coulomb plus Y-type linear potential. As a remarkable fact, we find that the quark confinement force in the 3Q system can be perfectly described only with Abelian variables in the maximally Abelian gauge, which we call “perfect Abelian dominance” of the quark confinement. Reference: “Perfect Abelian Dominance of Quark Confinement in SU(3) QCD”, N. Sakumichi and H. Suganuma, Phys. Rev. D90 (Rapid Communication) 111501(R) (2014).
        Speaker: Prof. Hideo Suganuma (Kyoto University)
        Slides
    • Weak Decays and Matrix Elements 406

      406

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: James Simone
      • 305
        Improving our determinations of the decay constant $f_B$ and the $B\to\pi \ell \nu$ semi-leptonic form factors using physical light quarks
        We report on updates of our $B$-physics program using domain-wall light quarks and nonperturbatively tuned, relativistic $b$-quarks by adding measurements obtained with dynamical and physical light quarks. We present progress towards improved determinations of the $B$-meson decay constant $f_B$, the ratio $f_{B_s}/f_B$, and the $B\to \pi \ell\nu$ semi-leptonic form factor. Our results are based on the RBC/UKQCD 2+1 flavor gauge field configurations with (M\"obius) domain-wall fermions and the Iwasaki gauge action at two lattice spacing of 0.086 fm and 0.11 fm.
        Speaker: Dr Taichi Kawanai (Forschungszentrum Jülich GmbH)
        Slides
      • 306
        Hadronic form factors for rare semi-leptonic $B$ decays
        We show first results for the computation of semi-leptonic form factors for rare $B$ decays focusing in particular on $B \to K^{(*)} \ell^+ \ell^-$ and $B_s \to \phi \ell^+ \ell^-$. Our setup uses domain-wall light quarks, relativistic $b$ quarks, and the RBC/UKQCD 2+1 flavor gauge field configurations with domain-wall fermions and the Iwasaki gauge action.
        Speaker: Mr Edwin Lizarazo (University of Southampton)
        Slides
      • 307
        $B-\bar B$ mixing with domain-wall light quarks and relativistic $b$-quarks
        Using our well established setup of domain-wall light quarks and nonperturbatively tuned, relativistic $b$-quarks we extend our $B$-physics program to the computation of $B$-meson bag parameters and their ratio $\xi$. We present first results based on the RBC/UKQCD 2+1 flavor gauge field configurations with domain-wall fermions and the Iwasaki gauge action at two lattice spacing of 0.086 fm and 0.11 fm.
        Speaker: Dr Oliver Witzel (University of Edinburgh)
        Slides
      • 308
        Radial distributions of the axial density and the $B^{*\prime}B\pi$ coupling
        We will present our results for the radial distributions of the axial density in the heavy-light B meson in the static limit of HQET: $\langle B(\vec{p}) | A_{\mu}(\vec{r}) | B^{*(\prime)}(\vec{p}^{\ \prime}, \lambda) \rangle$. Here, $r$ is the distance between the static heavy quark and the insertion of the axial current acting like a probe. Using the Generalized Eigenvalue Problem, we are able to isolate the first radial excitation of the vector $B^{*}$ meson and the associated density. After discussing the results, we will show how these distributions allow us to determine the $B^{*\prime}B\pi$ coupling at the kinematical point $q^2=0 \neq m_{B^{*\prime}}-m_B$. The simulations are performed using CLS gauge configurations with Nf=2 non-pertubatively $O(a)$-improved Wilson-Clover fermions and Heavy Quark Effective Theory in the static limit.
        Speaker: Dr Antoine Gérardin (Laboratoire de Physique Corpusculaire (LPC))
        Slides
    • 18:30
      Banquet Ohwada (Kobe Portopia Hotel)

      Ohwada

      Kobe Portopia Hotel

    • Applications Beyond QCD 404

      404

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Yasufumi Araki
      • 309
        Effective action for the Abelian Higgs model for a gauge-invariant implementation on optical lattices
        We present a gauge-invariant effective action for the Abelian Higgs model in 1+1 dimensions. It is constructed by integrating out the gauge field and then using the hopping parameter expansion. The latter is tested with Monte Carlo simulations for small values of the scalar self-coupling. In the opposite limit, at infinitely large self-coupling, the Higgs mode is frozen and the partition function can be written in terms of local tensors and the tensor renormalization group blocking can be applied. The numerical implementation requires truncations and the time continuum limit of the blocked transfer matrix can be obtained numerically. At zero gauge coupling and with a spin-1 truncation, the small volume energy spectrum is identical to the low energy spectrum of a two-species Bose-Hubbard model in the limit of large onsite repulsion. The procedure is extended to finite gauge coupling and we derive a spin-1 approximation of the Hamiltonian which involves terms corresponding to transitions among the two species in the Bose-Hubbard model. An optical lattice implementation involving a ladder structure is proposed.
        Speaker: Dr Alexei Bazavov (University of California, Riverside/University of Iowa)
        Slides
      • 310
        Schwinger Model Mass Anomalous Dimension
        The mass anomalous dimension for several gauge theories with an infrared fixed point has recently been determined using the mode number of the Dirac operator. In order to better understand the sources of systematic error in this method, we apply it to a simpler model, the massive Schwinger model with two flavours of fermions, where analytical results are available for comparison with the lattice data.
        Speaker: Liam Keegan (CERN)
        Slides
      • 311
        Study of the continuum limit of the Schwinger model using Wilson's lattice formulation
        We estimate the critical point of the Schwinger model in the theta vaccum by taking the continuum limit of Wilson's lattice formulation. The decorated tensor renormalization group which has been recently proposed by Dittrich et al. is employed for numerical investigation. We compare our numerical results with those of Byrnes et al., which were derived from Kogut-Susskind's lattice formulation employing the density matrix renormalization group.
        Speaker: Dr Yuya Shimizu (RIKEN Advanced Institute for Computational Science)
        Slides
      • 312
        Critical flavour number of the Thirring model in three dimensions
        The Thirring model is a four fermion theory with vector interaction. We study it in three dimensions, where it is closely related to QED and other models used to describe properties of graphene. In addition it is a good toy model to study chiral symmetry breaking, since a phase with broken chiral symmetry is present for the model with one fermion flavour. On the other hand, there is no such phase in the limit of infinitely many fermion flavours. Thus, a transition at some critical flavour number Nfc is expected, where the broken phase vanishes. The model was already studied with different methods, including Schwinger-Dyson, functional renormalization group and lattice approaches. Most studies agree that there is indeed a phase transition from a chirally symmetric phase to a spontaneously broken phase for a small number of fermion flavours. But there is no agreement on the critical flavour number and further details of the critical behaviour. Values of Nfc found in the literature usually range between 2 and 7. All earlier lattice studies were performed with staggered fermions, where it is questionable if the lattice model has the same chiral symmetry as the continuum version. We present new results from simulations of the Thirring model with SLAC fermions. With this choice, we can be sure to implement the same chiral symmetry of the continuum model. First estimates for the critical flavour number are given in this talk.
        Speaker: Mr Daniel Schmidt (Theoretisch-Physikalisches Institut, Jena University)
        Slides
      • 313
        Asymptotically free lattice gauge theory in five dimensions
        The first lattice formulation of Lifshitz-type gauge theories is presented. While the Lorentz-invariant Yang-Mills theory is not renormalizable in five dimensions, non-Abelian Lifshitz-type gauge theories are renormalizable and asymptotically free. We construct a lattice gauge action and numerically examine the continuum limit and the bulk phase structure.
        Speaker: Dr Takuya Kanazawa (RIKEN iTHES)
        Slides
      • 314
        The gradient flow in simple field theories
        The gradient flow has proved a valuable tool to the lattice community, with a range of applications to a variety of lattice calculations, from scale-setting to renormalisation. In this talk, I will focus on the gradient flow as a tool to suppress power-divergent mixing, a consequence of the hypercubic symmetry of the lattice regulator that is a particular difficulty for calculations of high moments of parton distribution functions. Provided the flow time is kept fixed in physical units, the gradient flow removes power-divergent mixing at the expense of introducing a new physical scale in the continuum. The smeared operator product expansion is a formalism that systematically accounts for this additional scale and connects nonperturbative calculations of flowed operators to continuum physics. Here I investigate the role of the gradient flow in suppressing power-divergent mixing and discuss the smeared operator product expansion for simple field theories.
        Speaker: Dr Christopher Monahan (University of Utah)
        Slides
    • Hadron Structure 403

      403

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Rajan Gupta
      • 315
        Connected contribution to hadron correlation functions
        We present an update of our lattice QCD calculation of the connected contribution to the hadronic vacuum polarization. We discuss different evaluations procedures (including the moment base approach and Pade fits) and the corresponding systematic errors.
        Speaker: Dr Alfonso Sastre (University of Wuppertal)
        Slides
      • 316
        Disconnected contribution to hadron correlation functions
        We report on our lattice QCD calculation of the disconnected contribution to various hadron correlation functions. We discuss different evaluation strategies, optimization of algorithmic parameters and we show first results at a coarse lattice spacing.
        Speaker: Dr Balint Toth (University of Wuppertal)
        Slides
      • 317
        Disconnected quark loop contributions to nucleon observables using $N_f=2$ twisted clover fermions at a physical value of the light quark mass
        We compute the disconnected quark loops contributions entering the determination of nucleon observables, by using an N_f = 2 ensemble of twisted mass fermions with a clover term at a pion mass $m_\pi~130$ MeV. We employ exact deflation and implement all calculations in QUDA, enabling us to achieve large statistics and a good signal.
        Speaker: Dr Alejandro Vaquero Avilés-Casco (INFN Sezione Milano Biccoca)
        Slides
      • 318
        Exploring the effects of open boundary conditions on baryonic observables
        We explore the effects of open boundary conditions, as used in the CLS $N_{\mathrm f}=2+1$ ensembles, on nucleon two- and three-point functions, with a view towards studying nucleon form factors on these ensembles. Particular emphasis is put on controlling systematic effects due to the boundaries and excited states, and to the approach to the physical point.
        Speaker: Dr Tim Harris (Helmholtz Institute Mainz)
        Slides
      • 319
        Improving the lattice axial vector current
        For Wilson and clover fermions traditional formulations of the axial vector current do not respect the continuum Ward identity which relates the divergence of that current to the pseudoscalar density. Here we propose to use an one-link axial vector current whose divergence exactly satisfies a lattice Ward identity, involving the pseudoscalar density and a number of irrelevant operators. We check in one-loop lattice perturbation theory with SLiNC fermion and gauge plaquette action that this is indeed the case including order O(a) effects. With those operators that axial Ward identity remains renormalisation invariant. First preliminary results of a nonperturbative check of the Ward identity are presented.
        Speaker: Dr Arwed Schiller (Leipzig University)
        Slides
      • 320
        A systematic study of excited-state effects on nucleon axial form factors
        We extend our study of excited-state effects on nucleon vector form factors to the case of the axial vector and pseudoscalar form factors. Combining information from a variety of different ratios of two- and three-point functions, we are able to extract the form factors $G_A$ and $G_P$ over a range of momentum transfers Q^2; together with the use of different methods to suppress excited-state contaminations this allows us to systematically study the effect of excited states.
        Speaker: Dr Georg von Hippel (University of Mainz)
        Slides
    • Nonzero Temperature and Density 401

      401

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Urs Heller
      • 321
        On the axial U(1) symmetry at finite temperature
        We report the results of our finite temperature investigation of the axial symmetry restoration at finite temperature. We simulated two flavors of domain-wall fermions at several volumes and lattice spacings. After taking into account the systematic errors from the violation of the Ginsparg-Wilson relation, our results show that in the chiral limit there is a strong suppression of the axial U(1) symmetry breaking measured using meson susceptibilities. This suppression is compatible with a no-breaking scenario at zero quark mass in the chirally symmetric phase. We observed that the contribution of the violations of the GW relation to the meson susceptibilities are much larger that what the residual mass measurement would suggest. We also show some insights on the sources of the violations related to the lowest eigenmodes of the Dirac operator.
        Speaker: Dr Guido Cossu (KEK High Energy Accelerator Research Organization)
        Slides
      • 322
        Study of the U(1)A symmetry restoration in two-flavor QCD at finite temperature with reweighed overlap fermions
        We study the U(1)A anomaly in two-flavor lattice QCD at finite temperature using the Ginsparg-Wilson fermions. The gauge configurations are generated with the Mobius domain-wall fermion at and above the critical temperature on 32^3x8 and 32^3x12 lattices. We apply the reweighting of the fermion determinant to that of domain-wall fermion satisfying the GW relation exactly. The results for low-lying eigenmodes are significantly different, indicating that even a small violation of the GW relation may strongly affect the low-modes. We find a strong suppression of the low-modes after the reweighting, which suggests vanishing U(1)A in the chiral limit.
        Speaker: Dr Akio Tomiya (Osaka university)
        Slides
      • 323
        The $U_A(1)$ anomaly in high temperature QCD with chiral fermions on the lattice
        The magnitude of the $U_A(1)$ symmetry breaking is expected to affect the nature of the $N_f=2$ QCD chiral phase transition. The explicit breaking of chiral symmetry due to realistic light quark mass is small, hence it is important to use chiral fermions on the lattice to understand the effect of $U_A(1)$ near the chiral crossover temperature, $T_c$. I present our latest results for the eigenvalue spectrum of 2+1 flavour QCD with dynamical M\"{o}bius domain wall fermions at finite temperature probed using the overlap fermion operator on a $32^3\times 8$ lattice. We observe that the $U_A(1)$ is broken near $T_c$ and the low-lying eigenvalues primarily contributing to it. We check how sensitive the low-lying eigenvalues are on the sea-light quark mass and the lattice volume. We also present comparison with the earlier independent results with domain wall fermions.
        Speaker: Dr sayantan sharma (BNL)
        Slides
      • 324
        Thermal modification of mesons and restoration of broken symmetries from spatial correlation functions with HISQ
        By using spatial correlation functions we study thermal modifications of meson states and the restoration of broken symmetries at finite temperature in 2+1 flavor QCD with the Highly Improved Staggered Quarks (HISQ) action. The spatial correlation functions provide a direct signal for the thermal modification of meson spectral functions and indicate the restoration of chiral and U_A(1) symmetries at finite temperature. We calculate spatial correlation functions for several meson states on lattices with aspect ratio Ns/Nt=4 in a wide temperature range of 110--170 MeV on Nt=8 and 150--400 MeV on Nt=12 with a physical strange quark mass (ms) and light quark masses (ml) set to ml=ms/20. In pion and kaon sectors we find that the thermal modifications become obvious even below the critical temperature (Tc) and the modification of positive parity states is more significant than those of the negative parity states. We also find that the spatial correlation functions of the vector parity partners degenerate at Tc, whereas those of the scalar parity partners show clear differences at Tc and degenerate only above 1.6Tc. This means that the breaking of the U_A(1) symmetry is still significant at Tc. This fact is also confirmed through simulations at the physical point (ml=ms/27). We also discuss several properties of open-strange and open-charm mesons. This is study on behalf of the BNL-Bielfeld Collaboration.
        Speaker: Dr Yu Maezawa (Yukawa Institute for Theoretical Physics, Kyoto University)
        Slides
      • 325
        Footprint of non-decoupling in chiral phase transition
        We study chiral phase transition of two flavor massless QCD in the presence of a finite breaking of U(1) axial symmetry using the corresponding linear sigma model (LSM), in which half of eight scalar fields acquire the mass proportional to the U(1) breaking. Naively, the model is expected to reduce to ordinary O(4) LSM in the infrared limit if the massive degrees of freedom decouple sufficiently fast. We examine this expectation through the analysis of the renormalization group flows of the model, and discuss possible non-decoupling effects in chiral phase transition.
        Speaker: Mr Tomomi Sato (The Graduate University for Advanced Studies)
        Slides
      • 326
        Determination of U_A (1) restoration from meson screening masses by using the entanglement PNJL model: Toward chiral regime
        We determine temperature (T) dependence of U_A (1) restoration from meson screening masses calculated with 2+1 flavor lattice QCD, using PNJL model with entanglement vertex. The entanglement PNJL (EPNJL) model exhibits the U_A (1) anomaly through the Kobayashi- Maskawa- ’t Hooft (KMT) interaction. T dependence of KMT interaction strength is then determined from the difference between pion and a0 meson screening masses. The strength is strongly suppressed around the pseudocritical temperature of chiral transition. This suppression is much stronger than that predicted by instanton liquid model. Using this T-dependent KMT interaction, we draw the Columbia plot near the physical point. In the light-quark chiral-limit with the strange quark mass (m_s) fixed at the physical value, the chiral transition becomes the second order. A tricritical point appears when m_s is slightly below the physical value. Finally, we calculate T dependence of other meson screening masses with the EPNJL model and compare the model results with lattice QCD ones to check the validity of the present model.
        Speaker: Mr Masahiro Ishii (Kyushu University)
        Slides
    • Nonzero Temperature and Density 406

      406

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Christof Gattringer
      • 327
        Hagedorn spectrum and equation of state of Yang-Mills theories
        We present a novel lattice calculation of the equation of state of SU(2) Yang-Mills theory in the confining phase. We show that a gas of massive, non-interacting glueballs describes remarkably well the results, provided that a bosonic closed-string model is used to derive an exponentially growing Hagedorn spectrum for the heavy glueball states with no free parameters. This effective model can be applied to SU(3) Yang-Mills theory and the theoretical prediction agrees nicely with the lattice results reported by Borsányi et al. in JHEP 07 (2012) 056.
        Speaker: Mr Alessandro Nada (Università di Torino)
        Slides
      • 328
        Higher order net baryon number cumulants in the strong coupling lattice QCD
        In the Beam Energy Scan Phase-1 (BES-1), net proton number fluctuations have been measured in search for the QCD critical point [1]. Finding the critical point is one of the most challenging subjects also in lattice QCD due to the sign problem. In the strong coupling limit of QCD, one can investigate the finite density region by employing the Auxiliary Field Monte-Carlo (AFMC) method [2] or monomer-dimer-polymer simulation [3] in which the sign problem becomes milder. We will report results on the baryon number fluctuations in the chiral limit with AFMC and discuss the influence of the mesonic fluctuations on the critical behavior of the cumulants. [1]L. Adamczyk, et al., [STAR Collaboration], Phys. Rev. Lett. 112 (2014) 032302. [2] T. Ichihara, A. Ohnishi, and T. Z. Nakano, PTEP 2014 (2014) 12, 123D02. [3] W. Unger and P. de Forcrand, J.Phys. G38 (2011) 124190.
        Speaker: Mr Terukazu Ichihara (Kyoto University)
        Slides
      • 329
        Polyakov line actions from SU(3) lattice gauge theory with dynamical fermions: first results via relative weights
        We ​apply the relative weights method to extract an effective Polyakov line action, at finite chemical potential, from an underlying SU(3) lattice gauge theory with dynamical fermions. The center-symmetry breaking terms in the effective theory are fit to a form suggested by the hopping-parameter expansion, and the effective action is solved at finite chemical potential by a mean field approach. We present preliminary results for staggered fermions at ma=0.3.
        Speaker: Roman Höllwieser (NMSU / VUT)
        Slides
      • 330
        Effective Polyakov loop models for QCD-like theories at finite chemical potential
        We study QCD-like theories at finite density with three-dimensional Polyakov theories for heavy quarks. These effective theories are derived by combined strong coupling and hopping expansion techniques as previously been used for QCD. In particular we investigate the cold and dense regimes of the phase diagrams where one either expects Bose-Einstein condensation of bosonic or a liquid-gas transition of fermionic baryons, depending on the theory. We explicitly verify that there is no Silver Blaze problem, i.e. that the density tends to zero in the zero temperature limit for all quark chemical potentials below half the diquark mass in two-color and G2-QCD. As the quark number density increases with larger chemical potentials we generally also observe that the Polyakov loop rises indicating the creation of deconfined quark matter.
        Speaker: Mr Philipp Scior (Technische Universitaet Darmstadt)
        Slides
      • 331
        G(2)-QCD at finite temperature and density
        G(2)-QCD is a QCD-like theory that can be simulated with standard Monte-Carlo methods at finite baryon density. It shares many properties with QCD, especially fermionic baryons. In addition also bosonic baryons are present in the theory. In the talk we review the effect of bosonic and fermionic baryons on the phase diagram at zero temperature and show evidence for a first oder nuclear matter transition. Furthermore we present our latest results on the phase diagram at finite temperature and density.
        Speaker: Dr Bjoern Wellegehausen (Theoretisch-Physikalisches Institut, FSU Jena)
        Slides
      • 332
        Phase diagram of the U(2)xU(2) scalar model in three dimensions
        If the chiral anomaly is effectively restored at finite temperature, the chiral phase transition of 2-flavor QCD with massless quarks is characterised by the symmetry breaking patter U_L(2) x U_R(2) -> U_V(2). Therefore, the order of the phase transition of the U(2) x U(2) scalar model in three dimensions has been extensively studied. The results of the perturbative renormalization group (RG) analyses are controversial. Up to 4th order of the perturbative expansion, there was no stable IR fixed point of the beta functions, therefore, the phase transition had been believed to be of first order. However there are reports that an infrared fixed point emerges at 5th order. Then the existence of the IR fixed point (or second order phase transition) is still under debate. In this talk, we present results of a Monte Carlo simulation of the U(2) x U(2) scalar model on three dimensional lattice. We will make clear that the system undergoes both second and first order phase transition depending on the choice of the model parameters. We represent a phase diagram in the model parameter space from the view point of the order of the phase transition. We make a comparison the phase diagram with the RG flow of the 5th order RG analyses. We will show both results are consistent and the IR fixed point of the RG flow likely to exist.
        Speaker: Dr Kazuhiko Kamikado (Riken)
        Slides
    • Physics Beyond the Standard Model 405

      405

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Noboru KAWAMOTO
      • 333
        Witten index and phase diagram of compactified N=1 supersymmetric Yang-Mills theory on the lattice
        Owing to confinement, the fundamental particles of N=1 Supersymmetric Yang-Mills (SYM) theory, gluons and gluinos, appear only in colourless bound states at zero temperature. Compactifying one space-time dimension with periodic boundary conditions for fermions preserves supersymmetry, and confinement is predicted to persist independently of the length of the compactified dimension. This scenario can be tested non-perturbatively with Monte-Carlo simulations on a lattice. SUSY is, however, broken on the lattice and can be recovered only in the continuum limit. The partition function of compactified N=1 SYM with periodic fermion boundary conditions corresponds to the Witten index. Therefore it can be used to test whether supersymmetry is realized on the lattice. Results of our recent numerical simulations will be presented, supporting the disappearance of the deconfinement transition in the supersymmetric limit and the restoration of SUSY at low energies.
        Speaker: Mr Stefano Piemonte (WWU Muenster)
        Slides
      • 334
        Supermultiplets of the N=1 supersymmetric Yang-Mills theory in the continuum limit
        The spectrum of N=1 supersymmetric Yang-Mills theory, calculated on the lattice, is presented. The masses have been determined on three different lattice spacings and extrapolated towards vanishing gluino mass. We present the extrapolation to the continuum limit which is consistent with the formation of degenerate supermultiplets.
        Speaker: Dr pietro giudice (muenster university)
        Slides
      • 335
        Canonical simulations of supersymmetric SU(N) Yang-Mills quantum mechanics
        The fermion loop formulation naturally separates partition functions into their canonical sectors. Here we discuss various strategies to make use of this for supersymmetric SU(N) Yang-Mills quantum mechanics obtained from dimensional reduction in various dimensions and present numerical results for the separate canonical sectors with fixed fermion numbers. We comment on potential problems due to the sign of the contributions from the fermions and due to flat directions. Finally, we discuss the possibility to extend the simulation strategies to higher dimensional gauge field theories, such as QCD.
        Speaker: Dr Urs Wenger (Albert Einstein Center for fundamental physics)
        Slides
      • 336
        New results from lattice N=4 supersymmetric Yang--Mills
        I will present results from numerical studies of maximally supersymmetric Yang--Mills theory, using a new improved lattice action. I will also summarize the improvement procedure, which modifies the moduli equations in order to lift U(1) flat directions without violating the exact supersymmetry preserved at non-zero lattice spacing by the lattice formulation. The resulting improved action leads to dramatically reduced lattice artifacts and much more rapid approach to the continuum limit, allowing us to investigate stronger couplings on accessible lattice volumes.
        Speaker: David Schaich (Syracuse University)
        Slides
      • 337
        Large-scale computation of the exponentially expanding universe in a simplified Lorentzian IIB matrix model
        The type IIB matrix model is conjectured to be a nonperturbative formulation of superstring theory.In 2011, the Lorentzian version of the model was studied for the first time by Monte Carlo simulation, and it was shown that (3+1)D expanding universe emerges from (9+1)D space-time predicted by superstring theory. Recently, a simplified model that describes the early time behaviors of the original model was studied with matrix size up to N=64, and it was suggested that the expansion of the (3+1)D universe is actually exponential, which is reminiscent of the inflation. In this work, we confirm this exponential expansion by simulating the simplified model with much larger matrix size up to N=512 using a large-scale parallel computer.
        Speaker: Mr Yuta Ito (SOKENDAI)
        Slides
      • 338
        Lattice and string worldsheet in AdS/CFT: a numerical study
        A possible discretization of the Green-Schwarz string in AdS5xS5 background is discussed and applied to the numerical study of the so-called cusp anomaly f(g). The latter, a pure function of the coupling g, when studied assuming the integrability of the underlying model, is an example of non-trivial interpolation smoothly connecting weak and strong coupling regimes and testing the AdS/CFT correspondence at a very deep level.
        Speaker: Dr Valentina Forini (Humboldt University Berlin)
        Slides
    • Theoretical Developments 402

      402

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Goro Ishiki
      • 339
        Lattice gradient flow with tree-level $O(a^4)$ improvement in pure Yang-Mills theory
        Recently, the Yang-Mills gradient flow method has continued to develop remarkably. The most successful application is demonstrated in an accurate determination of a reference scale. However, there is still room for improvement with respect to the lattice gradient flow, where some lattice artifacts are found to be non-negligible. Following a recent paper by Fodor et al. (arXiv:1406.0827), we examine several types of tree-level improvements on the flow action with various gauge actions in order to reduce the lattice discretization error. In this talk, we will present our numerical results for reference scale $t_0$ using tree-level, $O(a^4)$ improved lattice gradient flow including the rectangle term in the flow action, and discuss effects of the discretization error in the lattice gradient flow.
        Speaker: Mr Norihiko Kamata (Tohoku Univ)
        Slides
      • 340
        Renormalization constants of the lattice energy momentum tensor using the gradient flow
        We implement two different strategies for the non perturbative determination of the renormalized energy momentum tensor in SU(3) Yang-Mills on the lattice. Both strategies employ observables built with gauge fields that evolve according to the Yang Mills gradient flow. In the first case, we use observables computed along the flow in order to define suitable lattice Ward Identities from which the renormalization constants of the energy momentum tensor can be measured. In the second one, we show how to compute the renormalized energy momentum tensor using the small flow time expansion of properly chosen flowed observables. We show and discuss the numerical results coming from the application of the first method, as well as some preliminary data from the second one.
        Speaker: Mr Francesco Capponi (Plymouth University)
        Slides
      • 341
        (2+1)-flavor QCD Thermodynamics from the Gradient Flow
        Recently, we proposed a novel method to define and calculate the energy-momentum tensor (EMT) in lattice gauge theory on the basis of the Yang-Mills gradient flow. In this talk, we show the bulk thermodynamic quantities in lattice gauge theory using this method for (2+1)-flavor QCD. The entropy density of (2+1)-flavor QCD at fixed temperature are calculated. We also show the flow time dependence of the EMT including the dynamical fermion contributions. This work is based on a joint-collaboration between FlowQCD and WHOT QCD.
        Speaker: Dr Etsuko Itou (KEK)
        Slides
      • 342
        Background field method in the gradient flow
        We develop a background field method and associated gauge fixing prescription in gradient flow equations. By using this formulation, we carry out non-diagramatic one-loop computations of the small flow time expansion of gauge invariant composite operators, aiming for the application to the energy-momentum tensor and the flavor-singlet axial-vector current on the lattice.
        Speaker: Prof. Hiroshi Suzuki (Kyushu University)
        Slides
      • 343
        Generalized Gradient Flow Equation and Its Applications
        We generalize the gradient flow equation for field theories with nonlinearly realized symmetry. We apply the method to two theories, the super Yang-Mills theory in four dimensions and the O(N) nonlinear sigma model in two dimensions. Firstly, applying the formalism to super Yang-Mills theory, we construct the supersymmetric gradient flow equation. Furthermore, choosing an appropriate modification term to damp the gauge degrees of freedom, we obtain the gradient flow equation which is closed within the Wess-Zumino gauge. Secondly, applying the formalism to the O(N) nonlinear sigma model in two dimensions, we construct the gradient flow equation in the 1/N expansion. Solving this equation in the large N limit, we show non-perturbatively that the two point function at finite flow time is free from UV divergence. We also discuss the solution for four point function in the next-to-leading order of the 1/N expansion.
        Speaker: Dr Kengo Kikuchi (Yukawa Institute for Theoretical Physics, Kyoto University)
        Slides
      • 344
        The Wilson flow in scalar field theory
        A nonperturbative renormalization prescription for the energy-momentum tensor, based on space-time symmetries along the Wilson flow, has been proposed recently in the context of four-dimensional gauge theories. We extend this construction to the case of a scalar field theory, and investigate its numerical feasibility by studying Ward identities in 3-dimensional phi^4 theory. In this talk, we introduce the Wilson flow for the scalar theory, discuss its renormalization properties, and present some preliminary results about the determination of the renormalization constants for the energy-momentum tensor.
        Speaker: Ms Susanne Ehret (University of Edinburgh)
        Slides
    • 11:00
      Lunch Break
    • Plenary Session BF & 1F Main Hall

      BF & 1F Main Hall

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Kazuyuki Kanaya
      • 345
        QCD at non-zero temperature from the lattice
        I will review the status of calculations of thermodynamics quantities, spatial correlation lengths and real-time properties of strongly interacting matter at non-zero temperature. An attempt at a synthesis will be made. Quark number susceptibilities will not be covered here.
        Speaker: Prof. Harvey Meyer (Johannes Gutenberg-University Mainz)
        Slides
      • 346
        Fluctuations of conserved charges at finite temperature and density
        Fluctuations of conserved charges in a grand canonical ensemble can be calculated as derivatives of the free energy with respect to the respective chemical potential. They are directly related to experimentally available observables that describe the hadronization in heavy ion collisions. The same derivatives can be used to extrapolate zero density results to finite chemical potential. We review the recent lattice calculations in the staggered formalism and discuss its implications to phenomenology and resummed perturbation theory.
        Speaker: Dr Szabolcs Borsanyi (University of Wuppertal)
        Slides
      • 347
        The Lefschetz thimble and the sign problem
        In this talk I will review the proposal to formulate quantum field theories (QFTs) on a Lefschetz thimble, which was put forward to enable Monte Carlo simulations of lattice QFTs affected by a sign problem. First I will review the theoretical justification of the approach, and comment on some open issues. Then, I will review the algorithms that have been proposed and are being tested to represent and simulate a lattice QFT on a Lefschetz thimble. In particular, I will review the lessons from the very first models of QFTs that have been studied with this approach.
        Speaker: Dr Luigi Scorzato
        Slides
    • 14:45
      Coffee Break
    • Plenary Session: Conference Closing BF & 1F Main Hall

      BF & 1F Main Hall

      Kobe International Conference Center

      6-9-1 Minatojima-nakamachi, Chuo-ku, Kobe, Hyogo 650-0046, Japan
      Convener: Norman Christ
      • 348
        Recent progress in lattice supersymmetry -- from lattice gauge theory to black hole
        Non-perturbative investigations of supersymmetry and superstring theory through lattice simulations are promising research directions. In this talk, I will review the recent progress in lattice supersymmetry, in particular, the numerical verification of the gauge/gravity duality.
        Speaker: Dr Daisuke Kadoh (KEK)
      • 349
        Beyond the Standard Model: lattice calculations at the energy frontier
        The discovery of the Higgs boson in 2012 completed the Standard Model but many puzzles remain. Composite systems where the Higgs boson is a bound state of some new fermion species are viable models to describe beyond-SM phenomenology, but are most likely strongly coupled and require non-perturbative investigations. In this talk I will review recent lattice calculations that investigate general properties of near-conformal systems and might describe the Higgs boson.
        Speaker: Dr Anna Hasenfratz (University of Colorado)
        Slides
      • 16:30
        Conference Closing